Coronaviruses as DNA wannabes: a new model for the regulation of RNA virus replication fidelity

The NSP6-L260F substitution in SARS-CoV-2 BQ.1.1 and XBB.1.16 lineages compensates for the reduced viral polymerase activity caused by mutations in NSP13 and NSP14

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron variants emerged at the end of 2021, and their subvariants are still circulating worldwide. While changes in the S protein of these variants have been extensively studied, the roles of amino acid substitutions in non-structural proteins have not been fully revealed. In this study, we found that SARS-CoV-2 bearing the NSP6-L260F substitution emerged repeatedly when we generated several SARS-CoV-2 variants by reverse genetics or when we passaged SARS-CoV-2 isolated from clinical samples and that it was selected under cell culture conditions. Although this substitution has been detected in BQ.1.1 and XBB.1.16 that circulated in nature, its effect on viral properties is unclear. Here, we generated SARS-CoV-2 with or without the NSP6-L260F by reverse genetics and found that NSP6-L260F promotes virus replication in vitro and in vivo by increasing viral polymerase activity and enhancing virus pathogenicity in hamsters. We also identified disadvantageous substitutions, NSP13-M233I and NSP14-D222Y, that reduced BQ.1.1 and XBB.1.16 replication, respectively. These adverse effects were compensated for by NSP6-L260F. Our findings suggest the importance of NSP6-L260F for virus replication and pathogenicity and reveal part of the evolutionary process of Omicron variants.IMPORTANCEAlthough the properties of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron variants continue to change through the acquisition of various amino acid substitutions, the roles of the amino acid substitutions in the non-structural proteins have not been fully explored. In this study, we found that the NSP6-L260F substitution enhances viral polymerase activity and is important for viral replication and pathogenicity. In addition, we found that the NSP13-M233I substitution in the BQ.1.1 lineage and the NSP14-D222Y substitution in the XBB.1.16 lineage reduce viral polymerase activity, and this adverse effect is compensated for by the NSP6-L260F substitution. Our results provide insight into the evolutionary process of SARS-CoV-2.

Hidden evolutionary constraints dictate the retention of coronavirus accessory genes

Coronaviruses exhibit many mechanisms of genetic innovation, including the acquisition of accessory genes that originate by capture of cellular genes or through duplication of existing viral genes. Accessory genes influence viral host range and cellular tropism, but little is known about how selection acts on these variable regions of virus genomes. We used experimental evolution of mouse hepatitis virus (MHV) encoding a cellular AKAP7 phosphodiesterase and an inactive native phosphodiesterase, NS2, to model the evolutionary fate of accessory genes. After courses of serial infection, the gene encoding inactive NS2, ORF2, unexpectedly remained intact, suggesting it is under cryptic constraint uncoupled from the function of NS2. By contrast, AKAP7 was retained under strong selection but rapidly lost under relaxed selection. Experimental evolution also led to altered viral replication in a cell-type-specific manner and changed the relative proportions of subgenomic viral RNA in plaque-purified viral isolates, revealing additional mechanisms of adaptation. Guided by the retention of MHV ORF2 and similar patterns in related betacoronaviruses, we analyzed ORF8 of SARS-CoV-2, which is proposed to have arisen via gene duplication and contains premature stop codons in several globally successful lineages. As with MHV ORF2, the coding-defective SARS-CoV-2 ORF8 gene remained largely intact in these lineages, mirroring patterns observed during MHV experimental evolution, challenging assumptions on the dynamics of gene loss in virus genomes, and extending these findings to viruses currently adapting to humans.

Inactivation of β-coronavirus MHV-A59 by 2.8 GHz microwave

From the severe acute respiratory syndrome coronavirus in 2003 to the severe acute respiratory syndrome coronavirus 2 in 2019, coronavirus has seriously threatened human health. Electromagnetic waves not only own high penetration and low pollution but also can physically resonate with the virus. Several studies have demonstrated that electromagnetic waves can inactivate viruses efficiently. However, there is still a lack of systemic studies to analyze the potential factors closely associated with the effectiveness of inactivation, such as pH, temperature, and so on. In this study, we evaluated the inactivation ability of a 2.8 GHz microwave (MW) on MHV-A59, a substitute virus for coronavirus. Moreover, the influences of environmental pH and temperature on inactivation abilities were also discussed. The results showed that the viral morphology was destroyed, and the infectivity of MHV-A59 was significantly decreased after exposure to a 2.8 GHz MW at a density of 100 mW/cm2. Furthermore, alteration of pH 8 could produce synergistic effects with MW on virus inactivation. And, it was also proved that MWs could inactivate viruses better at room temperature than that under lower environmental temperatures. These results suggested that electromagnetic wave has great promise to become an effective tool to eliminate coronavirus.

Error thresholds in the presence of epistatic interactions

Models for viral populations with high replication error rates (such as RNA viruses) rely on the quasispecies concept, in which mutational pressure beyond the so-called "error threshold" leads to a loss of essential genetic information and population collapse, an effect known as the "error catastrophe." We explain how crossing this threshold, as a result of increasing mutation rates, can be understood as a second-order phase transition, even in the presence of lethal mutations. In particular, we show that, in fitness landscapes with a single peak, this collapse is equivalent to a ferroparamagnetic transition, where the back-mutation rate plays the role of the external magnetic field. We then generalize this framework to rugged fitness landscapes, like the ones that arise from epistatic interactions, and provide numerical evidence that there is a transition from a high average fitness regime to a low average fitness one, similarly to single-peaked landscapes. The onset of the transition is heralded by a sudden change in the susceptibility to variations in the mutation rate. We use insight from replica symmetry breaking mechanisms in spin glasses, in particular by considering the fluctuations of the genotype similarity distribution as the order parameter.

Hidden evolutionary constraints dictate the retention of coronavirus accessory genes

Coronaviruses exhibit many mechanisms of genetic innovation, including the acquisition of accessory genes that originate by capture of cellular genes or through duplication of existing viral genes. Accessory genes influence viral host range and cellular tropism, but little is known about how selection acts on these variable regions of virus genomes. We used experimental evolution of mouse hepatitis virus (MHV) encoding a cellular AKAP7 phosphodiesterase and an inactive native phosphodiesterase, NS2 to model the evolutionary fate of accessory genes. After courses of serial infection, the gene encoding inactive NS2, ORF2, unexpectedly remained intact, suggesting it is under cryptic constraint uncoupled from the function of NS2. In contrast, AKAP7 was retained under strong selection but rapidly lost under relaxed selection. Experimental evolution also led to altered viral replication in a cell type-specific manner and changed the relative proportions of subgenomic viral RNA in plaque-purified viral isolates, revealing additional mechanisms of adaptation. Guided by the retention of ORF2 and similar patterns in related betacoronaviruses, we analyzed ORF8 of SARS-CoV-2, which arose via gene duplication and contains premature stop codons in several globally successful lineages. As with MHV ORF2, the coding-defective SARS-CoV-2 ORF8 gene remains largely intact, mirroring patterns observed during MHV experimental evolution, challenging assumptions on the dynamics of gene loss in virus genomes and extending these findings to viruses currently adapting to humans.

A ~40-kb flavi-like virus does not encode a known error-correcting mechanism

It is commonly held that there is a fundamental relationship between genome size and error rate, manifest as a notional "error threshold" that sets an upper limit on genome sizes. The genome sizes of RNA viruses, which have intrinsically high mutation rates due to a lack of mechanisms for error correction, must therefore be small to avoid accumulating an excessive number of deleterious mutations that will ultimately lead to population extinction. The proposed exceptions to this evolutionary rule are RNA viruses from the order Nidovirales (such as coronaviruses) that encode error-correcting exonucleases, enabling them to reach genome lengths greater than 40 kb. The recent discovery of large-genome flavi-like viruses (Flaviviridae), which comprise genomes up to 27 kb in length yet seemingly do not encode exonuclease domains, has led to the proposal that a proofreading mechanism is required to facilitate the expansion of nonsegmented RNA virus genomes above 30 kb. Herein, we describe a ~40 kb flavi-like virus identified in a Haliclona sponge metatranscriptome that does not encode a known exonuclease. Structural analysis revealed that this virus may have instead captured cellular domains associated with nucleic acid metabolism that have not been previously found in RNA viruses. Phylogenetic inference placed this virus as a divergent pesti-like lineage, such that we have provisionally termed it "Maximus pesti-like virus." This virus represents an instance of a flavi-like virus achieving a genome size comparable to that of the Nidovirales and demonstrates that RNA viruses have evolved multiple solutions to overcome the error threshold.

Evolutionary deletions within the SARS-CoV-2 genome as signature trends for virus fitness and adaptation

Coronaviruses are large RNA viruses that can infect and spread among humans and animals. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), responsible for coronavirus disease 2019, has evolved since its first detection in December 2019. Deletions are a common occurrence in SARS-CoV-2 evolution, particularly in specific genomic sites, and may be associated with the emergence of highly competent lineages. While deletions typically have a negative impact on viral fitness, some persist and become fixed in viral populations, indicating that they may confer advantageous benefits for the virus's adaptive evolution. This work presents a literature review and data analysis on structural losses in the SARS-CoV-2 genome and the potential relevance of specific signatures for enhanced viral fitness and spread.

Puzzles, challenges, and information reservoir of SARS-CoV-2 quasispecies

Upon the emergence of SARS-CoV-2 in the human population, it was conjectured that for this coronavirus the dynamic intra-host heterogeneity typical of RNA viruses would be toned down. Nothing of this sort is observed. Here we review the main observations on the complexity and diverse composition of SARS-CoV-2 mutant spectra sampled from infected patients, within the framework of quasispecies dynamics. The analyses suggest that the information provided by myriads of genomic sequences within infected individuals may have a predictive value of the genomic sequences that acquire epidemiological relevance. Possibilities to reconcile the presence of broad mutant spectra in the large RNA coronavirus genome with its encoding a 3' to 5' exonuclease proofreading-repair activity are considered. Indeterminations in the behavior of individual viral genomes provide a benefit for the survival of the ensemble. We propose that this concept falls in the domain of "stochastic thinking," a notion that applies also to cellular processes, as a means for biological systems to face unexpected needs.

Energetic and frustration analysis of SARS-CoV-2 nucleocapsid protein mutations

The ongoing COVID-19 spreads worldwide with the ability to evolve in diverse human populations. The nucleocapsid (N) protein is one of the mutational hotspots in the SARS-CoV-2 genome. The N protein is an abundant RNA-binding protein critical for viral genome packaging. It comprises two large domains including the N-terminal domain (NTD) and the C-terminal domain (CTD) linked by the centrally located linker region. Mutations in N protein have been reported to increase the severity of disease by modulating viral transmissibility, replication efficiency as well as virulence properties of the virus in different parts of the world. To study the effect of N protein missense mutations on protein stability, function, and pathogenicity, we analyzed 228 mutations from each domain of N protein. Further, we have studied the effect of mutations on local residual frustration changes in N protein. Out of 228 mutations, 11 mutations were predicted to be deleterious and destabilized. Among these mutations, R32C, R191C, and R203 M mutations fall into disordered regions and show significant change in frustration state. Overall, this work reveals that by altering the energetics and residual frustration, N protein mutations might affect the stability, function, and pathogenicity of the SARS-CoV-2.

Metatranscriptomics analysis reveals a novel transcriptional and translational landscape during Middle East respiratory syndrome coronavirus infection

Among all RNA viruses, coronavirus RNA transcription is the most complex and involves a process termed "discontinuous transcription" that results in the production of a set of 3'-nested, co-terminal genomic and subgenomic RNAs during infection. While the expression of the classic canonical set of subgenomic RNAs depends on the recognition of a 6- to 7-nt transcription regulatory core sequence (TRS), here, we use deep sequence and metagenomics analysis strategies and show that the coronavirus transcriptome is even more vast and more complex than previously appreciated and involves the production of leader-containing transcripts that have canonical and noncanonical leader-body junctions. Moreover, by ribosome protection and proteomics analyses, we show that both positive- and negative-sense transcripts are translationally active. The data support the hypothesis that the coronavirus proteome is much vaster than previously noted in the literature.

Perspective Chapter: Real-Time Genomic Surveillance for SARS-CoV-2 on Center Stage

The course of the COVID-19 pandemic depends not only on how the SARS-CoV-2 virus mutates but on the actions taken to respond to it. Important public health decisions can only be taken if we know viral dynamics, viral variants distribution, and whether new variants are emerging that may be more transmissible or/and more virulent, displaying evasion to vaccines or antiviral treatments. This situation has put the use of different approaches, such as molecular techniques and real-time genomic sequencing, to support public health decision-making on center stage. To achieve this, robust programs based on: (i) diagnostic capacity; (ii) high-throughput sequencing technologies; and (iii) high-performance bioinformatic resources, need to be established. This chapter focuses on how SARS-CoV-2 evolved since its discovery and it summarizes the scientific efforts to obtain genomic data as the virus spread throughout the globe.

Structural evolution of Delta lineage of SARS-CoV-2

One of the main obstacles in prevention and treatment of COVID-19 is the rapid evolution of the SARS-CoV-2 Spike protein. Given that Spike is the main target of common treatments of COVID-19, mutations occurring at this virulent factor can affect the effectiveness of treatments. The B.1.617.2 lineage of SARS-CoV-2, being characterized by many Spike mutations inside and outside of its receptor-binding domain (RBD), shows high infectivity and relative resistance to existing cures. Here, utilizing a wide range of computational biology approaches, such as immunoinformatics, molecular dynamics (MD), analysis of intrinsically disordered regions (IDRs), protein-protein interaction analyses, residue scanning, and free energy calculations, we examine the structural and biological attributes of the B.1.617.2 Spike protein. Furthermore, the antibody design protocol of Rosetta was implemented for evaluation the stability and affinity improvement of the Bamlanivimab (LY-CoV55) antibody, which is not capable of interactions with the B.1.617.2 Spike. We observed that the detected mutations in the Spike of the B1.617.2 variant of concern can cause extensive structural changes compatible with the described variation in immunogenicity, secondary and tertiary structure, oligomerization potency, Furin cleavability, and drug targetability. Compared to the Spike of Wuhan lineage, the B.1.617.2 Spike is more stable and binds to the Angiotensin-converting enzyme 2 (ACE2) with higher affinity.

COVID-19: Challenges of Viral Variants

The COVID-19 pandemic has been accompanied by SARS-CoV-2 evolution and emergence of viral variants that have far exceeded initial expectations. Five major variants of concern (Alpha, Beta, Gamma, Delta, and Omicron) have emerged, each having both unique and overlapping amino acid substitutions that have affected transmissibility, disease severity, and susceptibility to natural or vaccine-induced immune responses and monoclonal antibodies. Several of the more recent variants appear to have evolved properties of immune evasion, particularly in cases of prolonged infection. Tracking of existing variants and surveillance for new variants are critical for an effective pandemic response.

Pseudotyped Viruses for Coronaviruses

Seven coronaviruses have been identified that can infect humans, four of which usually cause mild symptoms, including HCoV-229E, HCoV-NL63, HCoV-OC43, and HCoV-HKU1, three of which are lethal coronaviruses, named severe acute respiratory syndrome coronavirus, Middle East respiratory syndrome coronavirus, and severe acute respiratory syndrome coronavirus 2. Pseudotyped virus is an important tool in the field of human coronavirus research because it is safe, easy to prepare, easy to detect, and highly modifiable. In addition to the application of pseudotyped viruses in the study of virus infection mechanism, vaccine, and candidate antiviral drug or antibody evaluation and screening, pseudotyped viruses can also be used as an important platform for further application in the prediction of immunogenicity and antigenicity after virus mutation, cross-species transmission prediction, screening, and preparation of vaccine strains with better broad spectrum and antigenicity. Meanwhile, as clinical trials of various types of vaccines and post-clinical studies are also being carried out one after another, the establishment of a high-throughput and fully automated detection platform based on SARS-CoV-2 pseudotyped virus to further reduce the cost of detection and manual intervention and improve the efficiency of large-scale detection is also a demand for the development of SARS-CoV-2 pseudotyped virus.

A Computational Study on the Interaction of NSP10 and NSP14: Unraveling the RNA Synthesis Proofreading Mechanism in SARS-CoV-2, SARS-CoV, and MERS-CoV

The interaction of exoribonuclease (ExoN) nonstructural protein (NSP14) with NSP10 co-factors is crucial for high-fidelity proofreading activity of coronavirus replication and transcription. Proofreading function is critical for maintaining the large genomes to ensure replication proficiency; therefore, while maintaining the viral replication fitness, quick resistance has been reported to the nucleotide analogue (NA) drugs. Therefore, targeting the NSP14 and NSP10 interacting interface with small molecules or peptides could be a better strategy to obstruct replication processes of coronaviruses (CoVs). A comparative study on the binding mechanism of NSP10 with the NSP14 ExoN domain of SARS-CoV-2, SARS-CoV, MERS-CoV, and four SARS-CoV-2 NSP14mutant complexes has been carried out. Protein-protein interaction (PPI) dynamics, per-residue binding free energy (BFE) analyses, and the identification of interface hotspot residues have been studied using molecular dynamics simulations and various computational tools. The BFE of the SARS-CoV NSP14-NSP10 complex was higher when compared to novel SARS-CoV-2 and MERS. However, SARS-CoV-2 NSP14mutant systems display a higher BFE as compared to the wild type (WT) but lower than SARS-CoV and MERS. Despite the high BFE, the SARS-CoV NSP14-NSP10 complex appears to be structurally more flexible in many regions especially the catalytic site, which is not seen in SARS-CoV-2 and its mutant or MERS complexes. The significantly high residue energy contribution of key interface residues and hotspots reveals that the high binding energy between NSP14 and NSP10 may enhance the functional activity of the proofreading complex, as the NSP10-NSP14 interaction is essential in maintaining the stability of the ExoN domain for the replicative fitness of CoVs. The factors discussed for SARS-CoV-2 complexes may be responsible for NSP14 ExoN having a high replication proficiency, significantly leading to the evolution of new variants of SARS-CoV-2. The NSP14 residues V66, T69, D126, and I201and eight residues of NSP10 (L16, F19, V21, V42, M44, H80, K93, and F96) are identified as common hotspots. Overall, the interface area, hotspot locations, bonded/nonbonded contacts, and energies between NSP14 and NSP10 may pave a way in designing potential inhibitors to disrupt NSP14-NSP10 interactions of CoVs especially SARS-CoV-2.

Refolding of lid subdomain of SARS-CoV-2 nsp14 upon nsp10 interaction releases exonuclease activity

During RNA replication, coronaviruses require proofreading to maintain the integrity of their large genomes. Nsp14 associates with viral polymerase complex to excise the mismatched nucleotides. Aside from the exonuclease activity, nsp14 methyltransferase domain mediates cap methylation, facilitating translation initiation and protecting viral RNA from recognition by the innate immune sensors. The nsp14 exonuclease activity is modulated by a protein co-factor nsp10. While the nsp10/nsp14 complex structure is available, the mechanistic basis for nsp10-mediated modulation remains unclear in the absence of the nsp14 structure. Here, we provide a crystal structure of nsp14 in an apo-form. Comparative analysis of the apo- and nsp10-bound structures explain the modulatory role of the co-factor protein and reveal the allosteric nsp14 control mechanism essential for drug discovery. Further, the flexibility of the N-terminal lid of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) nsp14 structure presented in this study rationalizes the recently proposed idea of nsp14/nsp10/nsp16 ternary complex.

Molecular dynamics simulations of the SARS-CoV-2 Spike protein and variants of concern: structural evidence for convergent adaptive evolution

The Spike protein's structure of the SARS-CoV-2 provides a unique opportunity to consider perturbations at the atomic level. We used the cryo-electron microscopy structure of the open conformation of the Spike protein to assess the impact of the mutations observed in the variants of concern at the molecular level. Molecular dynamics were subsequently performed with both the wt and the mutated forms to compare the flexibility and variation data for each residue of the three-dimensional fluctuations in the region associated with each alpha carbon. Additionally, protein-protein docking was used to investigate the interaction of each mutated profile with the ACE-2 receptor. After the molecular dynamics, the results show that the mutations increased the stability of the trimeric protein, with greater stability observed in the Gamma variant harboring the 10 characteristic mutations. The results of molecular dynamics, as shown by RMSF demonstrated for the residues that comprise the binding domain receptor (RBD), exhibited a reduction in flexibility, which was more pronounced in the Gamma variant. Finally, protein-protein docking experiments revealed an increase in the number of hydrophobic interactions and hydrogen bonds in the Gamma variant against the ACE-2 receptor, as opposed to the other variants. Taken together, these in silico experiments suggest that the evolution of the mutations favored the increased stability of Spike protein while potentially improving its interaction with the ACE-2 receptor, which in turn may indicate putative structural outcomes of the selection of these mutations in the convergent adaptive evolution as it has been observed for SARS-CoV-2.Communicated by Ramaswamy H. Sarma.

Viral Nucleases from Herpesviruses and Coronavirus in Recombination and Proofreading: Potential Targets for Antiviral Drug Discovery

In this review, we explore recombination in two very different virus families that have become major threats to human health. The Herpesviridae are a large family of pathogenic double-stranded DNA viruses involved in a range of diseases affecting both people and animals. Coronaviridae are positive-strand RNA viruses (CoVs) that have also become major threats to global health and economic stability, especially in the last two decades. Despite many differences, such as the make-up of their genetic material (DNA vs. RNA) and overall mechanisms of genome replication, both human herpes viruses (HHVs) and CoVs have evolved to rely heavily on recombination for viral genome replication, adaptation to new hosts and evasion of host immune regulation. In this review, we will focus on the roles of three viral exonucleases: two HHV exonucleases (alkaline nuclease and PolExo) and one CoV exonuclease (ExoN). We will review the roles of these three nucleases in their respective life cycles and discuss the state of drug discovery efforts against these targets.

SARS-CoV-2 Mutant Spectra at Different Depth Levels Reveal an Overwhelming Abundance of Low Frequency Mutations

Populations of RNA viruses are composed of complex and dynamic mixtures of variant genomes that are termed mutant spectra or mutant clouds. This applies also to SARS-CoV-2, and mutations that are detected at low frequency in an infected individual can be dominant (represented in the consensus sequence) in subsequent variants of interest or variants of concern. Here we briefly review the main conclusions of our work on mutant spectrum characterization of hepatitis C virus (HCV) and SARS-CoV-2 at the nucleotide and amino acid levels and address the following two new questions derived from previous results: (i) how is the SARS-CoV-2 mutant and deletion spectrum composition in diagnostic samples, when examined at progressively lower cut-off mutant frequency values in ultra-deep sequencing; (ii) how the frequency distribution of minority amino acid substitutions in SARS-CoV-2 compares with that of HCV sampled also from infected patients. The main conclusions are the following: (i) the number of different mutations found at low frequency in SARS-CoV-2 mutant spectra increases dramatically (50- to 100-fold) as the cut-off frequency for mutation detection is lowered from 0.5% to 0.1%, and (ii) that, contrary to HCV, SARS-CoV-2 mutant spectra exhibit a deficit of intermediate frequency amino acid substitutions. The possible origin and implications of mutant spectrum differences among RNA viruses are discussed.

Next-Generation Sequencing for Confronting Virus Pandemics

Virus pandemics have happened, are happening and will happen again. In recent decades, the rate of zoonotic viral spillover into humans has accelerated, mirroring the expansion of our global footprint and travel network, including the expansion of viral vectors and the destruction of natural spaces, bringing humans closer to wild animals. Once viral cross-species transmission to humans occurs, transmission cannot be stopped by cement walls but by developing barriers based on knowledge that can prevent or reduce the effects of any pandemic. Controlling a local transmission affecting few individuals is more efficient that confronting a community outbreak in which infections cannot be traced. Genetic detection, identification, and characterization of infectious agents using next-generation sequencing (NGS) has been proven to be a powerful tool allowing for the development of fast PCR-based molecular assays, the rapid development of vaccines based on mRNA and DNA, the identification of outbreaks, transmission dynamics and spill-over events, the detection of new variants and treatment of vaccine resistance mutations, the development of direct-acting antiviral drugs, the discovery of relevant minority variants to improve knowledge of the viral life cycle, strengths and weaknesses, the potential for becoming dominant to take appropriate preventive measures, and the discovery of new routes of viral transmission.

How hydrolytic exoribonucleases impact human disease: Two sides of the same story

RNAs are extremely important molecules inside the cell which perform many different functions. For example, messenger RNAs, transfer RNAs, and ribosomal RNAs are involved in protein synthesis, whereas non-coding RNAs have numerous regulatory roles. Ribonucleases are the enzymes responsible for the processing and degradation of all types of RNAs, having multiple roles in every aspect of RNA metabolism. However, the involvement of RNases in disease is still not well understood. This review focuses on the involvement of the RNase II/RNB family of 3'-5' exoribonucleases in human disease. This can be attributed to direct effects, whereby mutations in the eukaryotic enzymes of this family (Dis3 (or Rrp44), Dis3L1 (or Dis3L), and Dis3L2) are associated with a disease, or indirect effects, whereby mutations in the prokaryotic counterparts of RNase II/RNB family (RNase II and/or RNase R) affect the physiology and virulence of several human pathogens. In this review, we will compare the structural and biochemical characteristics of the members of the RNase II/RNB family of enzymes. The outcomes of mutations impacting enzymatic function will be revisited, in terms of both the direct and indirect effects on disease. Furthermore, we also describe the SARS-CoV-2 viral exoribonuclease and its importance to combat COVID-19 pandemic. As a result, RNases may be a good therapeutic target to reduce bacterial and viral pathogenicity. These are the two perspectives on RNase II/RNB family enzymes that will be presented in this review.

Structure and dynamics of SARS-CoV-2 proofreading exoribonuclease ExoN

High-fidelity replication of the large RNA genome of coronaviruses (CoVs) is mediated by a 3'-to-5' exoribonuclease (ExoN) in nonstructural protein 14 (nsp14), which excises nucleotides including antiviral drugs misincorporated by the low-fidelity viral RNA-dependent RNA polymerase (RdRp) and has also been implicated in viral RNA recombination and resistance to innate immunity. Here, we determined a 1.6-Å resolution crystal structure of severe acute respiratory syndrome CoV 2 (SARS-CoV-2) ExoN in complex with its essential cofactor, nsp10. The structure shows a highly basic and concave surface flanking the active site, comprising several Lys residues of nsp14 and the N-terminal amino group of nsp10. Modeling suggests that this basic patch binds to the template strand of double-stranded RNA substrates to position the 3' end of the nascent strand in the ExoN active site, which is corroborated by mutational and computational analyses. We also show that the ExoN activity can rescue a stalled RNA primer poisoned with sofosbuvir and allow RdRp to continue its extension in the presence of the chain-terminating drug, biochemically recapitulating proofreading in SARS-CoV-2 replication. Molecular dynamics simulations further show remarkable flexibility of multidomain nsp14 and suggest that nsp10 stabilizes ExoN for substrate RNA binding to support its exonuclease activity. Our high-resolution structure of the SARS-CoV-2 ExoN-nsp10 complex serves as a platform for future development of anticoronaviral drugs or strategies to attenuate the viral virulence.

Reduction of Cell Fusion by Deletion in the Hypervariable Region of the Spike Protein of Mouse Hepatitis Virus

Deletions in the spike gene of mouse hepatitis virus (MHV) produce several variants with diverse biological characteristics, highlighting the significance of the spike gene in viral pathogenesis. In this study, we characterized the JHM-X strain, which has a deletion in the hypervariable region (HVR) of the spike gene, compared with the cl-2 strain, which has a full spike gene. Cytopathic effects (CPEs) induced by the two strains revealed that the size of the CPE produced by cl-2 is much greater than that produced by JHM-X in delayed brain tumor (DBT) cells. Thus, this finding explains the greater fusion activity of cl-2 than JHM-X in cultured cells, and we speculate that the deletion region of the spike protein is involved in the fusion activity differences. In contrast with the fusion activity, a comparison of the virus growth kinetics revealed that the titer of JHM-X was approximately 100 times higher than that of cl-2. We found that the deletion region of the spike protein was involved in fusion activity differences, whereas cl-2 produced significantly higher luciferase activity than JHM-X upon similar expression levels of the spike protein. However, the reason behind the growth difference is still unknown. Overall, we discovered that deletion in the HVR of the spike gene could be involved in the fusion activity differences between the two strains.

The antigenicity of SARS-CoV-2 Delta variants aggregated 10 high-frequency mutations in RBD has not changed sufficiently to replace the current vaccine strain

Emerging SARS-CoV-2 variants are the most serious problem for COVID-19 prophylaxis and treatment. To determine whether the SARS-CoV-2 vaccine strain should be updated following variant emergence like seasonal flu vaccine, the changed degree on antigenicity of SARS-CoV-2 variants and H3N2 flu vaccine strains was compared. The neutralization activities of Alpha, Beta and Gamma variants' spike protein-immunized sera were analysed against the eight current epidemic variants and 20 possible variants combining the top 10 prevalent RBD mutations based on the Delta variant, which were constructed using pseudotyped viruses. Meanwhile, the neutralization activities of convalescent sera and current inactivated and recombinant protein vaccine-elicited sera were also examined against all possible Delta variants. Eight HA protein-expressing DNAs elicited-animal sera were also tested against eight pseudotyped viruses of H3N2 flu vaccine strains from 2011-2019. Our results indicate that the antigenicity changes of possible Delta variants were mostly within four folds, whereas the antigenicity changes among different H3N2 vaccine strains were approximately 10-100-fold. Structural analysis of the antigenic characterization of the SARS-CoV-2 and H3N2 mutations supports the neutralization results. This study indicates that the antigenicity changes of the current SARS-CoV-2 may not be sufficient to require replacement of the current vaccine strain.

Within and Beyond the Nucleotide Addition Cycle of Viral RNA-dependent RNA Polymerases

Nucleotide addition cycle (NAC) is a fundamental process utilized by nucleic acid polymerases when carrying out nucleic acid biosynthesis. An induced-fit mechanism is usually taken by these polymerases upon NTP/dNTP substrate binding, leading to active site closure and formation of a phosphodiester bond. In viral RNA-dependent RNA polymerases, the post-chemistry translocation is stringently controlled by a structurally conserved motif, resulting in asymmetric movement of the template-product duplex. This perspective focuses on viral RdRP NAC and related mechanisms that have not been structurally clarified to date. Firstly, RdRP movement along the template strand in the absence of catalytic events may be relevant to catalytic complex dissociation or proofreading. Secondly, pyrophosphate or non-cognate NTP-mediated cleavage of the product strand 3′-nucleotide can also play a role in reactivating paused or arrested catalytic complexes. Furthermore, non-cognate NTP substrates, including NTP analog inhibitors, can not only alter NAC when being misincorporated, but also impact on subsequent NACs. Complications and challenges related to these topics are also discussed.

Discovery of Highly Potent Fusion Inhibitors with Potential Pan-Coronavirus Activity That Effectively Inhibit Major COVID-19 Variants of Concern (VOCs) in Pseudovirus-Based Assays

We report the discovery of several highly potent small molecules with low-nM potency against severe acute respiratory syndrome coronavirus (SARS-CoV; lowest half-maximal inhibitory concentration (IC₅₀: 13 nM), SARS-CoV-2 (IC₅₀: 23 nM), and Middle East respiratory syndrome coronavirus (MERS-CoV; IC₅₀: 76 nM) in pseudovirus-based assays with excellent selectivity index (SI) values (>5000), demonstrating potential pan-coronavirus inhibitory activities. Some compounds showed 100% inhibition against the cytopathic effects (CPE; IC₁₀₀) of an authentic SARS-CoV-2 (US_WA-1/2020) variant at 1.25 µM. The most active inhibitors also potently inhibited variants of concern (VOCs), including the UK (B.1.1.7) and South African (B.1.351) variants and the Delta variant (B.1.617.2) originally identified in India in pseudovirus-based assay. Surface plasmon resonance (SPR) analysis with one potent inhibitor confirmed that it binds to the prefusion SARS-CoV-2 spike protein trimer. These small-molecule inhibitors prevented virus-mediated cell-cell fusion. The absorption, distribution, metabolism, and excretion (ADME) data for one of the most active inhibitors, NBCoV1, demonstrated drug-like properties. An in vivo pharmacokinetics (PK) study of NBCoV1 in rats demonstrated an excellent half-life (t_1/2) of 11.3 h, a mean resident time (MRT) of 14.2 h, and oral bioavailability. We expect these lead inhibitors to facilitate the further development of preclinical and clinical candidates.

Microorganisms as Shapers of Human Civilization, from Pandemics to Even Our Genomes: Villains or Friends? A Historical Approach

Universal history is characterized by continuous evolution, in which civilizations are born and die. This evolution is associated with multiple factors, among which the role of microorganisms is often overlooked. Viruses and bacteria have written or decisively contributed to terrible episodes of history, such as the Black Death in 14th century Europe, the annihilation of pre-Columbian American civilizations, and pandemics such as the 1918 Spanish flu or the current COVID-19 pandemic caused by the coronavirus SARS-CoV-2. Nevertheless, it is clear that we could not live in a world without these tiny beings. Endogenous retroviruses have been key to our evolution and for the regulation of gene expression, and the gut microbiota helps us digest compounds that we could not otherwise process. In addition, we have used microorganisms to preserve or prepare food for millennia and more recently to obtain drugs such as antibiotics or to develop recombinant DNA technologies. Due to the enormous importance of microorganisms for our survival, they have significantly influenced the population genetics of different human groups. This paper will review the role of microorganisms as "villains" who have been responsible for tremendous mortality throughout history but also as "friends" who help us survive and evolve.

Host-dependent editing of SARS-CoV-2 in COVID-19 patients

A common trait among RNA viruses is their high capability to acquire genetic variability due to viral and host mechanisms. Next-generation sequencing (NGS) analysis enables the deep study of the viral quasispecies in samples from infected individuals. In this study, the viral quasispecies complexity and single nucleotide polymorphisms of the SARS-CoV-2 spike gene of coronavirus disease 2019 (COVID-19) patients with mild or severe disease were investigated using next-generation sequencing (Illumina platform). SARS-CoV-2 spike variability was higher in patients with long-lasting infection. Most substitutions found were present at frequencies lower than 1%, and had an A → G or T → C pattern, consistent with variants caused by adenosine deaminase acting on RNA-1 (ADAR1). ADAR1 affected a small fraction of replicating genomes, but produced multiple, mainly non-synonymous mutations. ADAR1 editing during replication rather than the RNA-dependent RNA polymerase (nsp12) was the predominant mechanism generating SARS-CoV-2 genetic variability. However, the mutations produced are not fixed in the infected human population, suggesting that ADAR1 may have an antiviral role, whereas nsp12-induced mutations occurring in patients with high viremia and persistent infection are the main source of new SARS-CoV-2 variants.

The seasonal behaviour of COVID-19 and its galectin-like culprit of the viral spike

Seasonal behaviour is an attribute of many viral diseases. Like other 'winter' RNA viruses, infections caused by the causative agent of COVID-19, SARS-CoV-2, appear to exhibit significant seasonal changes. Here we discuss the seasonal behaviour of COVID-19, emerging viral phenotypes, viral evolution, and how the mutational landscape of the virus affects the seasonal attributes of the disease. We propose that the multiple seasonal drivers behind infectious disease spread (and the spread of COVID-19 specifically) are in 'trade-off' relationships and can be better described within a framework of a 'triangle of viral persistence' modulated by the environment, physiology, and behaviour. This 'trade-off' exists as one trait cannot increase without a decrease in another. We also propose that molecular components of the virus can act as sensors of environment and physiology, and could represent molecular culprits of seasonality. We searched for flexible protein structures capable of being modulated by the environment and identified a galectin-like fold within the N-terminal domain of the spike protein of SARS-CoV-2 as a potential candidate. Tracking the prevalence of mutations in this structure resulted in the identification of a hemisphere-dependent seasonal pattern driven by mutational bursts. We propose that the galectin-like structure is a frequent target of mutations because it helps the virus evade or modulate the physiological responses of the host to further its spread and survival. The flexible regions of the N-terminal domain should now become a focus for mitigation through vaccines and therapeutics and for prediction and informed public health decision making.

In depth analysis of Cyprus-specific mutations of SARS-CoV-2 strains using computational approaches

BACKGROUND

This study aims to characterize SARS-CoV-2 mutations which are primarily prevalent in the Cypriot population. Moreover, using computational approaches, we assess whether these mutations are associated with changes in viral virulence.

METHODS

We utilize genetic data from 144 sequences of SARS-CoV-2 strains from the Cypriot population obtained between March 2020 and January 2021, as well as all data available from GISAID. We combine this with countries' regional information, such as deaths and cases per million, as well as COVID-19-related public health austerity measure response times. Initial indications of selective advantage of Cyprus-specific mutations are obtained by mutation tracking analysis. This entails calculating specific mutation frequencies within the Cypriot population and comparing these with their prevalence world-wide throughout the course of the pandemic. We further make use of linear regression models to extrapolate additional information that may be missed through standard statistical analysis.

RESULTS

We report a single mutation found in the ORF1ab gene (nucleotide position 18,440) that appears to be significantly enriched within the Cypriot population. The amino acid change is denoted as S6059F, which maps to the SARS-CoV-2 NSP14 protein. We further analyse this mutation using regression models to investigate possible associations with increased deaths and cases per million. Moreover, protein structure prediction tools show that the mutation infers a conformational change to the protein that significantly alters its structure when compared to the reference protein.

CONCLUSIONS

Investigating Cyprus-specific mutations for SARS-CoV-2 can lead to a better understanding of viral pathogenicity. Researching these mutations can generate potential links between viral-specific mutations and the unique genomics of the Cypriot population. This can not only lead to important findings from which to battle the pandemic on a national level, but also provide insights into viral virulence worldwide.

Tiled-ClickSeq for targeted sequencing of complete coronavirus genomes with simultaneous capture of RNA recombination and minority variants

High-throughput genomics of SARS-CoV-2 is essential to characterize virus evolution and to identify adaptations that affect pathogenicity or transmission. While single-nucleotide variations (SNVs) are commonly considered as driving virus adaption, RNA recombination events that delete or insert nucleic acid sequences are also critical. Whole genome targeting sequencing of SARS-CoV-2 is typically achieved using pairs of primers to generate cDNA amplicons suitable for next-generation sequencing (NGS). However, paired-primer approaches impose constraints on where primers can be designed, how many amplicons are synthesized and requires multiple PCR reactions with non-overlapping primer pools. This imparts sensitivity to underlying SNVs and fails to resolve RNA recombination junctions that are not flanked by primer pairs. To address these limitations, we have designed an approach called ‘Tiled-ClickSeq’, which uses hundreds of tiled-primers spaced evenly along the virus genome in a single reverse-transcription reaction. The other end of the cDNA amplicon is generated by azido-nucleotides that stochastically terminate cDNA synthesis, removing the need for a paired-primer. A sequencing adaptor containing a Unique Molecular Identifier (UMI) is appended to the cDNA fragment using click-chemistry and a PCR reaction generates a final NGS library. Tiled-ClickSeq provides complete genome coverage, including the 5’UTR, at high depth and specificity to the virus on both Illumina and Nanopore NGS platforms. Here, we analyze multiple SARS-CoV-2 isolates and clinical samples to simultaneously characterize minority variants, sub-genomic mRNAs (sgmRNAs), structural variants (SVs) and D-RNAs. Tiled-ClickSeq therefore provides a convenient and robust platform for SARS-CoV-2 genomics that captures the full range of RNA species in a single, simple assay.

Neutral evolution test of the spike protein of SARS-CoV-2 and its implications in the binding to ACE2

As the SARS-CoV-2 has spread and the pandemic has dragged on, the virus continued to evolve rapidly resulting in the emergence of new highly transmissible variants that can be of public health concern. The evolutionary mechanisms that drove this rapid diversity are not well understood but neutral evolution should open the first insight. The neutral theory of evolution states that most mutations in the nucleic acid sequences are random and they can be fixed or disappear by purifying selection. Herein, we performed a neutrality test to better understand the selective pressures exerted over SARS-CoV-2 spike protein from homologue proteins of Betacoronavirus, as well as to the spikes from human clinical isolates of the virus. Specifically, Tyr and Asn have higher occurrence rates on the Receptor Binding Domain (RBD) and in the overall sequence of spike proteins of Betacoronavirus, whereas His and Arg have lower occurrence rates. The in vivo evolutionary phenomenon of SARS-CoV-2 shows that Glu, Lys, Phe, and Val have the highest probability of occurrence in the emergent viral particles. Amino acids that have higher occurrence than the expected by the neutral control, are favorable and are fixed in the sequence while the ones that have lower occurrence than expected, influence the stability and/or functionality of the protein. Our results show that most unique mutations either for SARS-CoV-2 or its variants of health concern are under selective pressures, which could be related either to the evasion of the immune system, increasing the virus' fitness or altering protein - protein interactions with host proteins. We explored the consequences of those selected mutations in the structure and protein - protein interaction with the receptor. Altogether all these forces have shaped the spike protein and the continually evolving variants.

Mutation Rates, Mutation Frequencies, and Proofreading-Repair Activities in RNA Virus Genetics

The error rate displayed during template copying to produce viral RNA progeny is a biologically relevant parameter of the replication complexes of viruses. It has consequences for virus-host interactions, and it represents the first step in the diversification of viruses in nature. Measurements during infections and with purified viral polymerases indicate that mutation rates for RNA viruses are in the range of 10-3 to 10-6 copying errors per nucleotide incorporated into the nascent RNA product. Although viruses are thought to exploit high error rates for adaptation to changing environments, some of them possess misincorporation correcting activities. One of them is a proofreading-repair 3' to 5' exonuclease present in coronaviruses that may decrease the error rate during replication. Here we review experimental evidence and models of information maintenance that explain why elevated mutation rates have been preserved during the evolution of RNA (and some DNA) viruses. The models also offer an interpretation of why error correction mechanisms have evolved to maintain the stability of genetic information carried out by large viral RNA genomes such as the coronaviruses.

Tiled-ClickSeq for targeted sequencing of complete coronavirus genomes with simultaneous capture of RNA recombination and minority variants

High-throughput genomics of SARS-CoV-2 is essential to characterize virus evolution and to identify adaptations that affect pathogenicity or transmission. While single-nucleotide variations (SNVs) are commonly considered as driving virus adaption, RNA recombination events that delete or insert nucleic acid sequences are also critical. Whole genome targeting sequencing of SARS-CoV-2 is typically achieved using pairs of primers to generate cDNA amplicons suitable for Next-Generation Sequencing (NGS). However, paired-primer approaches impose constraints on where primers can be designed, how many amplicons are synthesized and requires multiple PCR reactions with non-overlapping primer pools. This imparts sensitivity to underlying SNVs and fails to resolve RNA recombination junctions that are not flanked by primer pairs. To address these limitations, we have designed an approach called 'Tiled-ClickSeq', which uses hundreds of tiled-primers spaced evenly along the virus genome in a single reverse-transcription reaction. The other end of the cDNA amplicon is generated by azido-nucleotides that stochastically terminate cDNA synthesis, removing the need for a paired-primer. A sequencing adaptor containing a Unique Molecular Identifier (UMI) is appended to the cDNA fragment using click-chemistry and a PCR reaction generates a final NGS library. Tiled-ClickSeq provides complete genome coverage, including the 5'UTR, at high depth and specificity to the virus on both Illumina and Nanopore NGS platforms. Here, we analyze multiple SARS-CoV-2 isolates and clinical samples to simultaneously characterize minority variants, sub-genomic mRNAs (sgmRNAs), structural variants (SVs) and D-RNAs. Tiled-ClickSeq therefore provides a convenient and robust platform for SARS-CoV-2 genomics that captures the full range of RNA species in a single, simple assay.

Binding of a pyrimidine RNA base-mimic to SARS-CoV-2 nonstructural protein 9

The coronaviral nonstructural protein 9 (Nsp9) is essential for viral replication; it is the primary substrate of Nsp12's pseudokinase domain within the viral replication transcription complex, an association that also recruits other components during different stages of RNA reproduction. In the unmodified state, Nsp9 forms an obligate homodimer via an essential GxxxG protein-interaction motif, but its ssRNA-binding mechanism remains unknown. Using structural biological techniques, here we show that a base-mimicking compound identified from a small molecule fragment screen engages Nsp9 via a tetrameric Pi-Pi stacking interaction that induces the formation of a parallel trimer-of-dimers. This oligomerization mechanism allows an interchange of "latching" N-termini, the charges of which contribute to a series of electropositive channels that suggests a potential interface for viral RNA. The identified pyrrolo-pyrimidine compound may also serve as a potential starting point for the development of compounds seeking to probe Nsp9's role within SARS-CoV-2 replication.

Immunopathogenesis of Different Emerging Viral Infections: Evasion, Fatal Mechanism, and Prevention

Different emerging viral infections may emerge in different regions of the world and pose a global pandemic threat with high fatality. Clarification of the immunopathogenesis of different emerging viral infections can provide a plan for the crisis management and prevention of emerging infections. This perspective article describes how an emerging viral infection evolves from microbial mutation, zoonotic and/or vector-borne transmission that progresses to a fatal infection due to overt viremia, tissue-specific cytotropic damage or/and immunopathology. We classified immunopathogenesis of common emerging viral infections into 4 categories: 1) deficient immunity with disseminated viremia (e.g., Ebola); 2) pneumocytotropism with/without later hyperinflammation (e.g., COVID-19); 3) augmented immunopathology (e.g., Hanta); and 4) antibody-dependent enhancement of infection with altered immunity (e.g., Dengue). A practical guide to early blocking of viral evasion, limiting viral load and identifying the fatal mechanism of an emerging viral infection is provided to prevent and reduce the transmission, and to do rapid diagnoses followed by the early treatment of virus neutralization for reduction of morbidity and mortality of an emerging viral infection such as COVID-19.

The enzymes for genome size increase and maintenance of large (+)RNA viruses

With sizes <50 kb, viral RNA genomes are at the crossroads of genetic, biophysical, and biochemical stability in their host cell. Here, we analyze the enzymatic assets accompanying large RNA genome viruses, mostly based on recent scientific advances in Coronaviridae. We argue that, in addition to the presence of an RNA exonuclease (ExoN), two markers for the large size of viral RNA genomes are (i) the presence of one or more RNA methyltransferases (MTases) and (ii) a specific architecture of the RNA-dependent RNA polymerase active site. We propose that RNA genome expansion and maintenance are driven by an evolutionary ménage-à-trois made of fast and processive RNA polymerases, RNA repair ExoNs, and RNA MTases that relates to the transition between RNA- to DNA-based life.

Severe Acute Respiratory Syndrome Coronavirus 2: The Emergence of Important Genetic Variants and Testing Options for Clinical Laboratories

Monitoring the spread of emerging severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants relies on rapid genetic testing of the viral genome. The sequencing method commonly called next-generation sequencing can identify virus variants. At times, for target-specific mutation detection, reverse transcriptase polymerase chain reaction is used to identify specific variants. The Centers for Disease Control and Prevention's national SARS-CoV-2 Strain Surveillance Program is a comprehensive, population-based U.S. surveillance system to monitor SARS-CoV-2 genes, identifying emerging SARS-CoV-2 variants to determine implications for coronavirus disease 2019 (COVID-19) diagnostics, therapy, and vaccines. This review describes the main viral variants of concern and their potential impacts and briefly describes testing strategies.

SARS-CoV-2 501Y.V2 variants lack higher infectivity but do have immune escape

The 501Y.V2 variants of SARS-CoV-2 containing multiple mutations in spike are now dominant in South Africa and are rapidly spreading to other countries. Here, experiments with 18 pseudotyped viruses showed that the 501Y.V2 variants do not confer increased infectivity in multiple cell types except for murine ACE2-overexpressing cells, where a substantial increase in infectivity was observed. Notably, the susceptibility of the 501Y.V2 variants to 12 of 17 neutralizing monoclonal antibodies was substantially diminished, and the neutralization ability of the sera from convalescent patients and immunized mice was also reduced for these variants. The neutralization resistance was mainly caused by E484K and N501Y mutations in the receptor-binding domain of spike. The enhanced infectivity in murine ACE2-overexpressing cells suggests the possibility of spillover of the 501Y.V2 variants to mice. Moreover, the neutralization resistance we detected for the 501Y.V2 variants suggests the potential for compromised efficacy of monoclonal antibodies and vaccines.

In silico approach for identifying natural lead molecules against SARS-COV-2

The life challenging COVID-19 disease caused by the SARS-CoV-2 virus has greatly impacted smooth survival worldwide since its discovery in December 2019. Currently, it is one of the major threats to humanity. Moreover, any specific drug or vaccine unavailability against COVID-19 forces to discover a new drug on an urgent basis. Viral cycle inhibition could be one possible way to prevent the further genesis of this viral disease, which can be contributed by drug repurposing techniques or screening of small bioactive natural molecules against already validated targets of COVID-19. The main protease (M^pro) responsible for producing functional proteins from polyprotein is an important key step for SARS-CoV-2 virion replication. Natural product or herbal based formulations are an important platform for potential therapeutics and lead compounds in the drug discovery process. Therefore, here we have screened >53,500 bioactive natural molecules from six different natural product databases against M^pro (PDB ID: 6LU7) of COVID-19 through computational study. Further, the top three molecules were subjected to pharmacokinetics evaluation, which is an important factor that reduces the drug failure rate. Moreover, the top three screened molecules (C00014803, C00006660, ANLT0001) were further validated by a molecular dynamics study under a condition similar to the physiological one. Relative binding energy analysis of three lead molecules indicated that C00014803 possess highest binding affinity among all three hits. These extensive studies can be a significant foundation for developing a therapeutic agent against COVID-19 through vet lab studies.

The variant gambit: COVID-19's next move

More than a year after its emergence, COVID-19, the disease caused by SARS-CoV-2, continues to plague the world and dominate our daily lives. Even with the development of effective vaccines, this coronavirus pandemic continues to cause a fervor with the identification of major new variants hailing from the United Kingdom, South Africa, Brazil, and California. Coupled with worries over a distinct mink strain that has caused human infections and potential for further mutations, SARS-CoV-2 variants bring concerns for increased spread and escape from both vaccine and natural infection immunity. Here, we outline factors driving SARS-CoV-2 variant evolution, explore the potential impact of specific mutations, examine the risk of further mutations, and consider the experimental studies needed to understand the threat these variants pose. In this review, Plante et al. examine SARS-CoV-2 variants including B.1.1.7 (UK), B.1.351 (RSA), P.1 (Brazil), and B.1.429 (California). They focus on what factors contribute to variant emergence, mutations in and outside the spike protein, and studies needed to understand the impact of variants on infection, transmission, and vaccine efficacy.

Structure and dynamics of SARS-CoV-2 proofreading exoribonuclease ExoN

High-fidelity replication of the large RNA genome of coronaviruses (CoVs) is mediated by a 3'-to-5' exoribonuclease (ExoN) in non-structural protein 14 (nsp14), which excises nucleotides including antiviral drugs mis-incorporated by the low-fidelity viral RNA-dependent RNA polymerase (RdRp) and has also been implicated in viral RNA recombination and resistance to innate immunity. Here we determined a 1.6-Å resolution crystal structure of SARS-CoV-2 ExoN in complex with its essential co-factor, nsp10. The structure shows a highly basic and concave surface flanking the active site, comprising several Lys residues of nsp14 and the N-terminal amino group of nsp10. Modeling suggests that this basic patch binds to the template strand of double-stranded RNA substrates to position the 3' end of the nascent strand in the ExoN active site, which is corroborated by mutational and computational analyses. Molecular dynamics simulations further show remarkable flexibility of multi-domain nsp14 and suggest that nsp10 stabilizes ExoN for substrate RNA-binding to support its exoribonuclease activity. Our high-resolution structure of the SARS-CoV-2 ExoN-nsp10 complex serves as a platform for future development of anti-coronaviral drugs or strategies to attenuate the viral virulence.

Theoretical Analysis Reveals the Cost and Benefit of Proofreading in Coronavirus Genome Replication

Severe acute respiratory syndrome coronaviruses have unusually large RNA genomes replicated by a multiprotein complex containing an RNA-dependent RNA polymerase (RdRp). Exonuclease activity enables the RdRp complex to remove wrongly incorporated bases via proofreading, a process not utilized by other RNA viruses. However, it is unclear why the RdRp complex needs proofreading and what the associated trade-offs are. Here we investigate the interplay among the accuracy, speed, and energetic cost of proofreading in the RdRp complex using a kinetic model and bioinformatics analysis. We find that proofreading nearly optimizes the rate of functional virus production. However, we find that further optimization would lead to a significant increase in the proofreading cost. Unexpected importance of the cost minimization is further supported by other global analyses. We speculate that cost optimization could help avoid cell defense responses. Thus, proofreading is essential for the production of functional viruses, but its rate is limited by energy costs.

The Worldwide Search for the New Mutations in the RNA-Directed RNA Polymerase Domain of SARS-CoV-2

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an RNA virus, responsible for the current pandemic outbreak. In total, 200 genomes of the SARS-CoV-2 strains from four host organisms have been analyzed. To investigate the presence of the new mutations in the RNA-directed RNA Polymerase (RdRp) of SARS-CoV-2, we analyzed sequences isolated from different hosts, with particular emphasis on human isolates. We performed a search for the new mutations of the RdRp proteins and study how those newly identified mutations could influence RdRp protein stability. Our results revealed 25 mutations in Rhinolophus sinicus, 1 in Mustela lutreola, 6 in Homo sapiens, and none in Mus musculus RdRp proteins of the SARS-CoV-2 isolates. We found that P323L is the most common stabilising radical mutation in human isolates. Also, we described several unique mutations, specific for studied hosts. Therefore, our data suggest that new and emerging variants of the SARS-CoV-2 RdRp have to be considered for the development of effective therapeutic agents and treatments.

Identification of SARS-CoV-2 spike mutations that attenuate monoclonal and serum antibody neutralization

Neutralizing antibodies against the SARS-CoV-2 spike (S) protein are a goal of COVID-19 vaccines and have received emergency use authorization as therapeutics. However, viral escape mutants could compromise efficacy. To define immune-selected mutations in the S protein, we exposed a VSV-eGFP-SARS-CoV-2-S chimeric virus, in which the VSV glycoprotein is replaced with the S protein, to 19 neutralizing monoclonal antibodies (mAbs) against the receptor-binding domain (RBD) and generated 50 different escape mutants. Each mAb had a unique resistance profile, although many shared residues within an epitope of the RBD. Some variants (e.g., S477N) were resistant to neutralization by multiple mAbs, whereas others (e.g., E484K) escaped neutralization by convalescent sera. Additionally, sequential selection identified mutants that escape neutralization by antibody cocktails. Comparing these antibody-mediated mutations with sequence variation in circulating SARS-CoV-2 revealed substitutions that may attenuate neutralizing immune responses in some humans and thus warrant further investigation.

Identification of SARS-CoV-2 spike mutations that attenuate monoclonal and serum antibody neutralization

Neutralizing antibodies against the SARS-CoV-2 spike (S) protein are a goal of COVID-19 vaccines and have received emergency use authorization as therapeutics. However, viral escape mutants could compromise efficacy. To define immune-selected mutations in the S protein, we exposed a VSV-eGFP-SARS-CoV-2-S chimeric virus, in which the VSV glycoprotein is replaced with the S protein, to 19 neutralizing monoclonal antibodies (mAbs) against the receptor-binding domain (RBD) and generated 50 different escape mutants. Each mAb had a unique resistance profile, although many shared residues within an epitope of the RBD. Some variants (e.g., S477N) were resistant to neutralization by multiple mAbs, whereas others (e.g., E484K) escaped neutralization by convalescent sera. Additionally, sequential selection identified mutants that escape neutralization by antibody cocktails. Comparing these antibody-mediated mutations with sequence variation in circulating SARS-CoV-2 revealed substitutions that may attenuate neutralizing immune responses in some humans and thus warrant further investigation.

High Throughput Sequencing for the Detection and Characterization of RNA Viruses

This review aims to assess and recommend approaches for targeted and agnostic High Throughput Sequencing of RNA viruses in a variety of sample matrices. HTS also referred to as deep sequencing, next generation sequencing and third generation sequencing; has much to offer to the field of environmental virology as its increased sequencing depth circumvents issues with cloning environmental isolates for Sanger sequencing. That said however, it is important to consider the challenges and biases that method choice can impart to sequencing results. Here, methodology choices from RNA extraction, reverse transcription to library preparation are compared based on their impact on the detection or characterization of RNA viruses.

Potential Pathogenicity Determinants Identified from Structural Proteomics of SARS-CoV and SARS-CoV-2

Despite SARS-CoV and SARS-CoV-2 being equipped with highly similar protein arsenals, the corresponding zoonoses have spread among humans at extremely different rates. The specific characteristics of these viruses that led to such distinct outcomes remain unclear. Here, we apply proteome-wide comparative structural analysis aiming to identify the unique molecular elements in the SARS-CoV-2 proteome that may explain the differing consequences. By combining protein modeling and molecular dynamics simulations, we suggest nonconservative substitutions in functional regions of the spike glycoprotein (S), nsp1, and nsp3 that are contributing to differences in virulence. Particularly, we explain why the substitutions at the receptor-binding domain of S affect the structure–dynamics behavior in complexes with putative host receptors. Conservation of functional protein regions within the two taxa is also noteworthy. We suggest that the highly conserved main protease, nsp5, of SARS-CoV and SARS-CoV-2 is part of their mechanism of circumventing the host interferon antiviral response. Overall, most substitutions occur on the protein surfaces and may be modulating their antigenic properties and interactions with other macromolecules. Our results imply that the striking difference in the pervasiveness of SARS-CoV-2 and SARS-CoV among humans seems to significantly derive from molecular features that modulate the efficiency of viral particles in entering the host cells and blocking the host immune response.

SARS-CoV-2 infection remodels the host protein thermal stability landscape

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a global threat to human health and has compromised economic stability. In addition to the development of an effective vaccine, it is imperative to understand how SARS-CoV-2 hijacks host cellular machineries on a system-wide scale so that potential host-directed therapies can be developed. In situ proteome-wide abundance and thermal stability measurements using thermal proteome profiling (TPP) can inform on global changes in protein activity. Here we adapted TPP to high biosafety conditions amenable to SARS-CoV-2 handling. We discovered pronounced temporal alterations in host protein thermostability during infection, which converged on cellular processes including cell cycle, microtubule and RNA splicing regulation. Pharmacological inhibition of host proteins displaying altered thermal stability or abundance during infection suppressed SARS-CoV-2 replication. Overall, this work serves as a framework for expanding TPP workflows to globally important human pathogens that require high biosafety containment and provides deeper resolution into the molecular changes induced by SARS-CoV-2 infection.

COVID-19: A Review on the Novel Coronavirus Disease Evolution, Transmission, Detection, Control and Prevention

Three major outbreaks of the coronavirus, a zoonotic virus known to cause respiratory disease, have been reported since 2002, including SARS-CoV, MERS-CoV and the most recent 2019-nCoV, or more recently known as SARS-CoV-2. Bats are known to be the primary animal reservoir for coronaviruses. However, in the past few decades, the virus has been able to mutate and adapt to infect humans, resulting in an animal-to-human species barrier jump. The emergence of a novel coronavirus poses a serious global public health threat and possibly carries the potential of causing a major pandemic outbreak in the naïve human population. The recent outbreak of COVID-19, the disease caused by SARS-CoV-2, in Wuhan, Hubei Province, China has infected over 36.5 million individuals and claimed over one million lives worldwide, as of 8 October 2020. The novel virus is rapidly spreading across China and has been transmitted to 213 other countries/territories across the globe. Researchers have reported that the virus is constantly evolving and spreading through asymptomatic carriers, further suggesting a high global health threat. To this end, current up-to-date information on the coronavirus evolution and SARS-CoV-2 modes of transmission, detection techniques and current control and prevention strategies are summarized in this review.