Epistasis at the SARS-CoV-2 RBD Interface and the Propitiously Boring Implications for Vaccine Escape

At the time of this writing, December 2021, potential emergence of vaccine escape variants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a grave global concern. The interface between the receptor-binding domain (RBD) of SARS-CoV-2 spike (S) protein and the host receptor (ACE2) overlap with the binding site of principal neutralizing antibodies (NAb), limiting the repertoire of viable mutations. Nonetheless, variants with multiple mutations in the RBD have rose to dominance. Non-additive, epistatic relationships among RBD mutations are apparent, and assessing the impact of such epistasis on the mutational landscape is crucial. Epistasis can substantially increase the risk of vaccine escape and cannot be completely characterized through the study of the wild type (WT) alone. We employed protein structure modeling using Rosetta to compare the effects of all single mutants at the RBD-NAb and RBD-ACE2 interfaces for the WT, Delta, Gamma, and Omicron variants. Overall, epistasis at the RBD interface appears to be limited and the effects of most multiple mutations are additive. Epistasis at the Delta variant interface weakly stabilizes NAb interaction relative to ACE2 interaction, whereas in the Gamma variant, epistasis more substantially destabilizes NAb interaction. Although a small, systematic trend towards NAb destabilization not observed for Delta or Gamma was detected for Omicron, and despite bearing significantly more RBD mutations, the epistatic landscape of the Omicron variant closely resembles that of Gamma. These results suggest that, although Omicron poses new risks not observed with Delta, structural constraints on the RBD hamper continued evolution towards more complete vaccine escape. The modest ensemble of mutations relative to the WT that are currently known to reduce vaccine efficacy is likely to comprise the majority of all possible escape mutations for future variants, predicting continued efficacy of the existing vaccines. SignificanceEmergence of vaccine escape variants of SARS-CoV-2 is arguably the most pressing problem during the COVID-19 pandemic as vaccines are distributed worldwide. We employed a computational approach to assess the risk of antibody escape resulting from mutations in the receptor-binding domain of the spike protein of the wild type SARS-CoV-2 virus as well as the Delta, Gamma, and Omicron variants. At the time of writing, December, 2021, Omicron is poised to replace Delta as the dominant variant worldwide. The efficacy of the existing vaccines against Omicron could be substantially reduced relative to the WT and the potential for vaccine escape is of grave concern. Our results suggest that although Omicron poses new evolutionary risks not observed for the Delta variant, structural constraints on the RBD make continued evolution towards more complete vaccine escape unlikely. The modest set of escape-enhancing mutations already identified for the wild type likely include the majority of all possible mutations with this effect.

Insertions in SARS-CoV-2 genome caused by template switch and duplications give rise to new variants of potential concern

The appearance of multiple new SARS-CoV-2 variants during the winter of 2020-2021 is a matter of grave concern. Some of these new variants, such as B.1.617.2, B.1.1.7, and B.1.351, manifest higher infectivity and virulence than the earlier SARS-CoV-2 variants, with potential dramatic effects on the course of the COVID-19 pandemic. So far, analysis of new SARS-CoV-2 variants focused primarily on point nucleotide substitutions and short deletions that are readily identifiable by comparison to consensus genome sequences. In contrast, insertions have largely escaped the attention of researchers although the furin site insert in the spike protein is thought to be a determinant of SARS-CoV-2 virulence and other inserts might have contributed to coronavirus pathogenicity as well. Here, we investigate insertions in SARS-CoV-2 genomes and identify 347 unique inserts of different lengths. We present evidence that these inserts reflect actual virus variance rather than sequencing errors. Two principal mechanisms appear to account for the inserts in the SARS-CoV-2 genomes, polymerase slippage and template switch that might be associated with the synthesis of subgenomic RNAs. We show that inserts in the Spike glycoprotein can affect its antigenic properties and thus merit monitoring. At least, three inserts in the N-terminal domain of the Spike (ins245IME, ins246DSWG, and ins248SSLT) that were first detected in 2021 are predicted to lead to escape from neutralizing antibodies, whereas other inserts might result in escape from T-cell immunity.

Lost in translation: codon optimization inactivates SARS-CoV-2 RdRp

The catalytic subunit of SARS-CoV-2 RNA-dependent RNA polymerase (RdRp), Nsp12, has a unique NiRAN domain that transfers nucleoside monophosphates to the Nsp9 protein. The NiRAN and RdRp modules form a dynamic interface distant from their catalytic sites and both activities are essential for viral replication. We report that codon-optimized (for the pause-free translation) Nsp12 exists in inactive state in which NiRAN/RdRp interactions are broken, whereas translation by slow ribosomes and incubation with accessory Nsp7/8 subunits or NTPs partially rescue RdRp activity. Our data show that adenosine and remdesivir triphosphates promote synthesis of A-less RNAs, as does ppGpp, while amino acid substitutions at the NiRAN/RdRp interface augment activation, suggesting that ligand binding to the NiRAN catalytic site modulates RdRp activity. The existence of allosterically-linked nucleotidyl transferase sites that utilize the same substrates has important implications for understanding the mechanism of SARS-CoV-2 replication and design of its inhibitors. HighlightsO_LICodon-optimization of Nsp12 triggers misfolding and activity loss C_LIO_LISlow translation, accessory Nsp7 and Nsp8 subunits, and NTPs rescue Nsp12 C_LIO_LINon-substrate nucleotides activate RNA chain synthesis, likely via NiRAN domain C_LIO_LICrosstalk between two Nsp12 active sites that bind the same ligands C_LI

Evolution of Human Respiratory Virus Epidemics

BackgroundWhile pathogens often evolve towards reduced virulence, many counterexamples are evident. When faced with a new pathogen, such as SARS-CoV-2, it is highly desirable to be able to forecast the case fatality rate (CFR) into the future. Considerable effort has been invested towards the development of a mathematical framework for predicting virulence evolution. Although these approaches accurately recapitulate some complex outcomes, most rely on an assumed trade-off between mortality and infectivity. It is often impractical to empirically validate this constraint for human pathogens. ResultsUsing a compartment model with parameters tuning the degree to which symptomatic individuals are isolated and the duration of immunity, we reveal kinetic constraints where the variation of multiple parameters in concert leads to decreased virulence and increased pathogen fitness, whereas independent variation of the parameters decreases pathogen fitness. Smallpox, SARS-CoV-2, and Influenza are analyzed as diverse representatives of human respiratory viruses. We show that highly virulent viruses, such as Smallpox, are likely often constrained by host behavior, whereas moderately virulent viruses, such as SARS-CoV-2, appear to be typically constrained by the relationship between the duration of immunity and CFR. ConclusionsThe evolution of human respiratory epidemics appears to be often kinetically constrained and a reduction in virulence should not be assumed. Our findings imply that, without continued public health intervention, SARS-CoV-2 is likely to continue presenting a substantial disease burden. The existence of a parameter regime admitting endemic equilibrium suggests that herd immunity is unachievable. However, we demonstrate that even partial isolation of symptomatic individuals can have a major effect not only by reducing the number of fatalities in the short term but also by potentially changing the evolutionary trajectory of the virus towards reduced virulence.

Ongoing Adaptive Evolution and Globalization of Sars-Cov-2

Understanding the trends in SARS-CoV-2 evolution is paramount to control the COVID- 19 pandemic. We analyzed more than 300,000 high quality genome sequences of SARS-CoV-2 variants available as of January 2021. The results show that the ongoing evolution of SARS-CoV-2 during the pandemic is characterized primarily by purifying selection, but a small set of sites appear to evolve under positive selection. The receptor-binding domain of the spike protein and the nuclear localization signal (NLS) associated region of the nucleocapsid protein are enriched with positively selected amino acid replacements. These replacements form a strongly connected network of apparent epistatic interactions and are signatures of major partitions in the SARS-CoV-2 phylogeny. Virus diversity within each geographic region has been steadily growing for the entirety of the pandemic, but analysis of the phylogenetic distances between pairs of regions reveals four distinct periods based on global partitioning of the tree and the emergence of key mutations. The initial period of rapid diversification into region- specific phylogenies that ended in February 2020 was followed by a major extinction event and global homogenization concomitant with the spread of D614G in the spike protein, ending in March 2020. The NLS associated variants across multiple partitions rose to global prominence in March-July, during a period of stasis in terms of inter- regional diversity. Finally, beginning July 2020, multiple mutations, some of which have since been demonstrated to enable antibody evasion, began to emerge associated with ongoing regional diversification, which might be indicative of speciation. SignificanceUnderstanding the ongoing evolution of SARS-CoV-2 is essential to control and ultimately end the pandemic. We analyzed more than 300,000 SARS-CoV-2 genomes available as of January 2021 and demonstrate adaptive evolution of the virus that affects, primarily, multiple sites in the spike and nucleocapsid protein. Selection appears to act on combinations of mutations in these and other SARS-CoV-2 genes. Evolution of the virus is accompanied by ongoing adaptive diversification within and between geographic regions. This diversification could substantially prolong the pandemic and the vaccination campaign, in which variant-specific vaccines are likely to be required.

Prediction of the virus incubation period for COVID-19 and future outbreaks

A crucial factor in mitigating respiratory viral outbreaks is early determination of the duration of the incubation period and, accordingly, the required quarantine time for potentially exposed individuals. Here, we explore different genomic features of RNA viruses that correlate with the incubation times and provide a predictive model that accurately estimates the upper limit incubation time for diverse viruses including SARS-CoV-2, and thus, could help control future outbreaks.

Genomic determinants of pathogenicity in SARS-CoV-2 and other human coronaviruses

SARS-CoV-2 poses an immediate, major threat to public health across the globe. Here we report an in-depth molecular analysis to reconstruct the evolutionary origins of the enhanced pathogenicity of SARS-CoV-2 and other coronaviruses that are severe human pathogens. Using integrated comparative genomics and machine learning techniques, we identify key genomic features that differentiate SARS-CoV-2 and the viruses behind the two previous deadly coronavirus outbreaks, SARS-CoV and MERS-CoV, from less pathogenic coronaviruses. These features include enhancement of the nuclear localization signals in the nucleocapsid protein and distinct inserts in the spike glycoprotein that appear to be associated with high case fatality rate of these coronaviruses as well as the host switch from animals to humans. The identified features could be crucial elements of coronavirus pathogenicity and possible targets for diagnostics, prognostication and interventions.