Spike mutation D614G alters SARS-CoV-2 fitness

A pseudovirus system enables deep mutational scanning of the full SARS-CoV-2 spike

A major challenge in understanding SARS-CoV-2 evolution is interpreting the antigenic and functional effects of emerging mutations in the viral spike protein. Here, we describe a deep mutational scanning platform based on non-replicative pseudotyped lentiviruses that directly quantifies how large numbers of spike mutations impact antibody neutralization and pseudovirus infection. We apply this platform to produce libraries of the Omicron BA.1 and Delta spikes. These libraries each contain ∼7,000 distinct amino acid mutations in the context of up to ∼135,000 unique mutation combinations. We use these libraries to map escape mutations from neutralizing antibodies targeting the receptor-binding domain, N-terminal domain, and S2 subunit of spike. Overall, this work establishes a high-throughput and safe approach to measure how ∼105 combinations of mutations affect antibody neutralization and spike-mediated infection. Notably, the platform described here can be extended to the entry proteins of many other viruses.

Structural dynamics in the evolution of SARS-CoV-2 spike glycoprotein

SARS-CoV-2 spike glycoprotein mediates receptor binding and subsequent membrane fusion. It exists in a range of conformations, including a closed state unable to bind the ACE2 receptor, and an open state that does so but displays more exposed antigenic surface. Spikes of variants of concern (VOCs) acquired amino acid changes linked to increased virulence and immune evasion. Here, using HDX-MS, we identified changes in spike dynamics that we associate with the transition from closed to open conformations, to ACE2 binding, and to specific mutations in VOCs. We show that the RBD-associated subdomain plays a role in spike opening, whereas the NTD acts as a hotspot of conformational divergence of VOC spikes driving immune evasion. Alpha, beta and delta spikes assume predominantly open conformations and ACE2 binding increases the dynamics of their core helices, priming spikes for fusion. Conversely, substitutions in omicron spike lead to predominantly closed conformations, presumably enabling it to escape antibodies. At the same time, its core helices show characteristics of being pre-primed for fusion even in the absence of ACE2. These data inform on SARS-CoV-2 evolution and omicron variant emergence.

Characteristics of Cardiac Abnormalities in Pediatric Patients With Acute COVID-19

Introduction Coronavirus disease 2019 (COVID-19) is known to cause cardiac abnormalities in adults. Cardiac abnormalities are well-described in multisystem inflammatory syndrome in children, but effects in children with acute COVID-19 are less understood. In this multicenter study, we assessed the cardiac effects of acute COVID-19 among hospitalized children (<21 years) admitted to three large healthcare systems in New York City. Methods We performed a retrospective observational study. We examined electrocardiograms, echocardiograms, troponin, or B-type natriuretic peptides. Results Of 317 admitted patients, 131 (41%) underwent cardiac testing with 56 (43%) demonstrating cardiac abnormalities. Electrocardiogram abnormalities were the most common (46/117 patients (39%)), including repolarization abnormalities and QT prolongation. Elevated troponin occurred in 14/77 (18%) patients and B-type natriuretic peptide in 8/39 (21%) patients. Ventricular dysfunction was identified in 5/27 (19%) patients with an echocardiogram, all of whom had elevated troponin. Ventricular dysfunction resolved by first outpatient follow-up. Conclusion Electrocardiogram and troponin can assist clinicians in identifying children at risk for cardiac injury in acute COVID-19.

SARS-CoV-2 Spike Protein Post-Translational Modification Landscape and Its Impact on Protein Structure and Function via Computational Prediction

To elucidate the role of post-translational modifications (PTMs) in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike protein's structure and virulence, we generated a high-resolution map of 87 PTMs using liquid chromatography with tandem mass spectrometry data on the extracted spike protein from SARS-CoV-2 virions and then reconstituted its structure heterogeneity caused by PTMs. Nonetheless, Alphafold2, a high-accuracy artificial intelligence tool to perform protein structure prediction, relies solely on primary amino acid sequence, whereas the impact of PTM, which often modulates critical protein structure and function, is much ignored. To overcome this challenge, we proposed the mutagenesis approach-an in silico, site-directed amino acid substitution to mimic the influence of PTMs on protein structure due to altered physicochemical properties in the post-translationally modified amino acids-and then reconstituted the spike protein's structure from the substituted sequences by Alphafold2. For the first time, the proposed method revealed predicted protein structures resulting from PTMs, a problem that Alphafold2 has yet to address. As an example, we performed computational analyses of the interaction of the post-translationally modified spike protein with its host factors such as angiotensin-converting enzyme 2 to illuminate binding affinity. Mechanistically, this study suggested the structural analysis of post-translationally modified protein via mutagenesis and deep learning. To summarize, the reconstructed spike protein structures showed that specific PTMs can be used to modulate host factor binding, guide antibody design, and pave the way for new therapeutic targets. The code and Supplementary Materials are freely available at https://github.com/LTZHKUSTGZ/SARS-CoV-2-spike-protein-PTM.

A potent and broad neutralization of SARS-CoV-2 variants of concern by DARPins

We report the engineering and selection of two synthetic proteins-FSR16m and FSR22-for the possible treatment of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection. FSR16m and FSR22 are trimeric proteins composed of DARPin SR16m or SR22 fused with a T4 foldon. Despite selection by a spike protein from a now historical SARS-CoV-2 strain, FSR16m and FSR22 exhibit broad-spectrum neutralization of SARS-CoV-2 strains, inhibiting authentic B.1.351, B.1.617.2 and BA.1.1 viruses, with respective IC50 values of 3.4, 2.2 and 7.4 ng ml-1 for FSR16m. Cryo-EM structures revealed that these DARPins recognize a region of the receptor-binding domain (residues 456, 475, 486, 487 and 489) overlapping a critical portion of the angiotensin-converting enzyme 2 (ACE2)-binding surface. K18-hACE2 transgenic mice inoculated with B.1.617.2 and receiving intranasally administered FSR16m showed less weight loss and 10-100-fold lower viral burden in upper and lower respiratory tracts. The strong and broad neutralization potency makes FSR16m and FSR22 promising candidates for the prevention and treatment of infection by SARS-CoV-2.

Global SARS-CoV-2 genomic surveillance: What we have learned (so far)

The COVID-19 pandemic has brought significant challenges for genomic surveillance strategies in public health systems worldwide. During the past thirty-four months, many countries faced several epidemic waves of SARS-CoV-2 infections, driven mainly by the emergence and spread of novel variants. In that line, genomic surveillance has been a crucial toolkit to study the real-time SARS-CoV-2 evolution, for the assessment and optimization of novel diagnostic assays, and to improve the efficacy of existing vaccines. During the pandemic, the identification of emerging lineages carrying lineage-specific mutations (particularly those in the Receptor Binding domain) showed how these mutations might significantly impact viral transmissibility, protection from reinfection and vaccination. So far, an unprecedented number of SARS-CoV-2 viral genomes has been released in public databases (i.e., GISAID, and NCBI), achieving 14 million genome sequences available as of early-November 2022. In the present review, we summarise the global landscape of SARS-CoV-2 during the first thirty-four months of viral circulation and evolution. It demonstrates the urgency and importance of sustained investment in genomic surveillance strategies to timely identify the emergence of any potential viral pathogen or associated variants, which in turn is key to epidemic and pandemic preparedness.

The race to understand immunopathology in COVID-19: Perspectives on the impact of quantitative approaches to understand within-host interactions

The COVID-19 pandemic has revealed the need for the increased integration of modelling and data analysis to public health, experimental, and clinical studies. Throughout the first two years of the pandemic, there has been a concerted effort to improve our understanding of the within-host immune response to the SARS-CoV-2 virus to provide better predictions of COVID-19 severity, treatment and vaccine development questions, and insights into viral evolution and the impacts of variants on immunopathology. Here we provide perspectives on what has been accomplished using quantitative methods, including predictive modelling, population genetics, machine learning, and dimensionality reduction techniques, in the first 26 months of the COVID-19 pandemic approaches, and where we go from here to improve our responses to this and future pandemics.

Positive-strand RNA viruses-a Keystone Symposia report

Positive-strand RNA viruses have been the cause of several recent outbreaks and epidemics, including the Zika virus epidemic in 2015, the SARS outbreak in 2003, and the ongoing SARS-CoV-2 pandemic. On June 18-22, 2022, researchers focusing on positive-strand RNA viruses met for the Keystone Symposium "Positive-Strand RNA Viruses" to share the latest research in molecular and cell biology, virology, immunology, vaccinology, and antiviral drug development. This report presents concise summaries of the scientific discussions at the symposium.

Global landscape of SARS-CoV-2 mutations and conserved regions

BACKGROUND

At the end of December 2019, a novel strain of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) disease (COVID-19) has been identified in Wuhan, a central city in China, and then spread to every corner of the globe. As of October 8, 2022, the total number of COVID-19 cases had reached over 621 million worldwide, with more than 6.56 million confirmed deaths. Since SARS-CoV-2 genome sequences change due to mutation and recombination, it is pivotal to surveil emerging variants and monitor changes for improving pandemic management.

METHODS

10,287,271 SARS-CoV-2 genome sequence samples were downloaded in FASTA format from the GISAID databases from February 24, 2020, to April 2022. Python programming language (version 3.8.0) software was utilized to process FASTA files to identify variants and sequence conservation. The NCBI RefSeq SARS-CoV-2 genome (accession no. NC_045512.2) was considered as the reference sequence.

RESULTS

Six mutations had more than 50% frequency in global SARS-CoV-2. These mutations include the P323L (99.3%) in NSP12, D614G (97.6) in S, the T492I (70.4) in NSP4, R203M (62.8%) in N, T60A (61.4%) in Orf9b, and P1228L (50.0%) in NSP3. In the SARS-CoV-2 genome, no mutation was observed in more than 90% of nsp11, nsp7, nsp10, nsp9, nsp8, and nsp16 regions. On the other hand, N, nsp3, S, nsp4, nsp12, and M had the maximum rate of mutations. In the S protein, the highest mutation frequency was observed in aa 508-635(0.77%) and aa 381-508 (0.43%). The highest frequency of mutation was observed in aa 66-88 (2.19%), aa 7-14, and aa 164-246 (2.92%) in M, E, and N proteins, respectively.

CONCLUSION

Therefore, monitoring SARS-CoV-2 proteomic changes and detecting hot spots mutations and conserved regions could be applied to improve the SARS-CoV-2 diagnostic efficiency and design safe and effective vaccines against emerging variants.

Combined Use of RT-qPCR and NGS for Identification and Surveillance of SARS-CoV-2 Variants of Concern in Residual Clinical Laboratory Samples in Miami-Dade County, Florida

Over the course of the COVID-19 pandemic, SARS-CoV-2 variants of concern (VOCs) with increased transmissibility and immune escape capabilities, such as Delta and Omicron, have triggered waves of new COVID-19 infections worldwide, and Omicron subvariants continue to represent a global health concern. Tracking the prevalence and dynamics of VOCs has clinical and epidemiological significance and is essential for modeling the progression and evolution of the COVID-19 pandemic. Next generation sequencing (NGS) is recognized as the gold standard for genomic characterization of SARS-CoV-2 variants, but it is labor and cost intensive and not amenable to rapid lineage identification. Here we describe a two-pronged approach for rapid, cost-effective surveillance of SARS-CoV-2 VOCs by combining reverse-transcriptase quantitative polymerase chain reaction (RT-qPCR) and periodic NGS with the ARTIC sequencing method. Variant surveillance by RT-qPCR included the commercially available TaqPath COVID-19 Combo Kit to track S-gene target failure (SGTF) associated with the spike protein deletion H69-V70, as well as two internally designed and validated RT-qPCR assays targeting two N-terminal-domain (NTD) spike gene deletions, NTD156-7 and NTD25-7. The NTD156-7 RT-qPCR assay facilitated tracking of the Delta variant, while the NTD25-7 RT-qPCR assay was used for tracking Omicron variants, including the BA.2, BA.4, and BA.5 lineages. In silico validation of the NTD156-7 and NTD25-7 primers and probes compared with publicly available SARS-CoV-2 genome databases showed low variability in regions corresponding to oligonucleotide binding sites. Similarly, in vitro validation with NGS-confirmed samples showed excellent correlation. RT-qPCR assays allow for near-real-time monitoring of circulating and emerging variants allowing for ongoing surveillance of variant dynamics in a local population. By performing periodic sequencing of variant surveillance by RT-qPCR methods, we were able to provide ongoing validation of the results obtained by RT-qPCR screening. Rapid SARS-CoV-2 variant identification and surveillance by this combined approach served to inform clinical decisions in a timely manner and permitted better utilization of sequencing resources.

Interference reduction isothermal nucleic acid amplification strategy for COVID-19 variant detection

Common reference methods for COVID-19 variant diagnosis include viral sequencing and PCR-based methods. However, sequencing is tedious, expensive, and time-consuming, while PCR-based methods have high risk of insensitive detection in variant-prone regions and are susceptible to potential background signal interference in biological samples. Here, we report a loop-mediated interference reduction isothermal nucleic acid amplification (LM-IR-INA) strategy for highly sensitive single-base mutation detection in viral variants. This strategy exploits the advantages of nicking endonuclease-mediated isothermal amplification, luminescent iridium(III) probes, and time-resolved emission spectroscopy (TRES). Using the LM-IR-INA strategy, we established a luminescence platform for diagnosing COVID-19 D796Y single-base substitution detection with a detection limit of 2.01 × 10⁵ copies/μL in a linear range of 6.01 × 10⁵ to 3.76 × 10⁸ copies/μL and an excellent specificity with a variant/wild-type ratio of significantly less than 0.0625%. The developed TRES-based method was also successfully applied to detect D796Y single-base substitution sequence in complicated biological samples, including throat and blood, and was a superior to steady-state technique. LM-IR-INA was also demonstrated for detecting the single-base substitution D614G as well as the multiple-base mutation H69/V70del without mutual interference, indicating that this approach has the potential to be used as a universal viral variant detection strategy.

Ancestral Lineage of SARS-CoV-2 Is More Stable in Human Biological Fluids than Alpha, Beta, and Omicron Variants of Concern

SARS-CoV-2 is a zoonotic virus first identified in 2019, and has quickly spread worldwide. The virus is primarily transmitted through respiratory droplets from infected persons; however, the virus-laden excretions can contaminate surfaces which can serve as a potential source of infection. Since the beginning of the pandemic, SARS-CoV-2 has continued to evolve and accumulate mutations throughout its genome leading to the emergence of variants of concern (VOCs) which exhibit increased fitness, transmissibility, and/or virulence. However, the stability of SARS-CoV-2 VOCs in biological fluids has not been thoroughly investigated. The aim of this study was to determine and compare the stability of different SARS-CoV-2 strains in human biological fluids. Here, we demonstrate that the ancestral strain of the Wuhan-like lineage A was more stable than the Alpha VOC B.1.1.7, and the Beta VOC B.1.351 strains in human liquid nasal mucus and sputum. In contrast, there was no difference in stability among the three strains in dried biological fluids. Furthermore, we also show that the Omicron VOC B.1.1.529 strain was less stable than the ancestral Wuhan-like strain in liquid nasal mucus. These studies provide insight into the effect of the molecular evolution of SARS-CoV-2 on environmental virus stability, which is important information for the development of countermeasures against SARS-CoV-2. IMPORTANCE Genetic evolution of SARS-CoV-2 leads to the continuous emergence of novel virus variants, posing a significant concern to global public health. Five of these variants have been classified to date into variants of concern (VOCs); Alpha, Beta, Gamma, Delta, and Omicron. Previous studies investigated the stability of SARS-CoV-2 under various conditions, but there is a gap of knowledge on the survival of SARS-CoV-2 VOCs in human biological fluids which are clinically relevant. Here, we present evidence that Alpha, Beta, and Omicron VOCs were less stable than the ancestral Wuhan-like strain in human biological fluids. Our findings highlight the potential risk of contaminated human biological fluids in SARS-CoV-2 transmission and contribute to the development of countermeasures against SARS-CoV-2.

SARSMutOnto: An Ontology for SARS-CoV-2 Lineages and Mutations

Mutations allow viruses to continuously evolve by changing their genetic code to adapt to the hosts they infect. It is an adaptive and evolutionary mechanism that helps viruses acquire characteristics favoring their survival and propagation. The COVID-19 pandemic declared by the WHO in March 2020 is caused by the SARS-CoV-2 virus. The non-stop adaptive mutations of this virus and the emergence of several variants over time with characteristics favoring their spread constitute one of the biggest obstacles that researchers face in controlling this pandemic. Understanding the mutation mechanism allows for the adoption of anticipatory measures and the proposal of strategies to control its propagation. In this study, we focus on the mutations of this virus, and we propose the SARSMutOnto ontology to model SARS-CoV-2 mutations reported by Pango researchers. A detailed description is given for each mutation. The genes where the mutations occur and the genomic structure of this virus are also included. The sub-lineages and the recombinant sub-lineages resulting from these mutations are additionally represented while maintaining their hierarchy. We developed a Python-based tool to automatically generate this ontology from various published Pango source files. At the end of this paper, we provide some examples of SPARQL queries that can be used to exploit this ontology. SARSMutOnto might become a 'wet bench' machine learning tool for predicting likely future mutations based on previous mutations.

Arrayed Imaging Reflectometry monitoring of anti-viral antibody production throughout vaccination and breakthrough Covid-19

Immune responses to COVID-19 infection and vaccination are individual and varied. There is a need to understand the timeline of vaccination efficacy against current and yet to be discovered viral mutations. Assessing immunity to SARS-CoV-2 in the context of immunity to other respiratory viruses is also valuable. Here we demonstrate the capability of a fully automated prototype Arrayed Imaging Reflectometry system to perform reliable longitudinal serology against a 34-plex respiratory array. The array contains antigens for respiratory syncytial virus, seasonal influenza, common human coronaviruses, MERS, SARS-CoV-1, and SARS-CoV-2. AIR measures a change in reflectivity due to the binding of serum antibodies to the antigens on the array. Samples were collected from convalescent COVID-19 donors and individuals vaccinated with a two-dose mRNA vaccine regimen. Vaccinated samples were collected prior to the first dose, one week after the first dose, one week after the second dose, and monthly thereafter. Information following booster dose and/or breakthrough infection is included for a subset of subjects. Longitudinal samples of vaccinated individuals demonstrate a rise and fall of SARS-CoV-2 spike antibodies in agreement with general knowledge of the adaptive immune response and other studies. Linear Regression analysis was performed to understand the relationship between antibodies binding to different antigens on the array. Our analysis identified strong correlations between closely related influenza virus strains as well as correlations between SARS-CoV-2, SARS-CoV-1, and human coronavirus 229E. A small test of using diluted whole blood from a fingerstick provided clean arrays with antibody binding comparable to serum. Potential applications include assessing immunity in the context of exposure to multiple respiratory viruses, clinical serology, population monitoring to facilitate public health recommendations, and vaccine development against new viruses and virus mutations.

A Protein Co-Conservation Network Model Characterizes Mutation Effects on SARS-CoV-2 Spike Protein

The emergence of numerous variants of SARS-CoV-2 has presented challenges to the global efforts to control the COVID-19 pandemic. The major mutation is in the SARS-CoV-2 viral envelope spike protein that is responsible for virus attachment to the host, and is the main target for host antibodies. It is critically important to study the biological effects of the mutations to understand the mechanisms of how mutations alter viral functions. Here, we propose a protein co-conservation weighted network (PCCN) model only based on the protein sequence to characterize the mutation sites by topological features and to investigate the mutation effects on the spike protein from a network view. Frist, we found that the mutation sites on the spike protein had significantly larger centrality than the non-mutation sites. Second, the stability changes and binding free energy changes in the mutation sites were positively significantly correlated with their neighbors' degree and the shortest path length separately. The results indicate that our PCCN model provides new insights into mutations on spike proteins and reflects the mutation effects on protein function alternations.

More or less deadly? A mathematical model that predicts SARS-CoV-2 evolutionary direction

SARS-CoV-2 has caused tremendous deaths globally. It is of great value to predict the evolutionary direction of SARS-CoV-2. In this paper, we proposed a novel mathematical model that could predict the evolutionary trend of SARS-CoV-2. We focus on the mutational effects on viral assembly capacity. A robust coarse-grained mathematical model is constructed to simulate the virus dynamics in the host body. Both virulence and transmissibility can be quantified in this model. A delicate equilibrium point that optimizes the transmissibility can be numerically obtained. Based on this model, the virulence of SARS-CoV-2 might further decrease, accompanied by an enhancement of transmissibility. However, this trend is not continuous; its virulence will not disappear but remains at a relatively stable range. A virus assembly model which simulates the virus packing process is also proposed. It can be explained why a few mutations would lead to a significant divergence in clinical performance, both in the overall particle formation quantity and virulence. This research provides a novel mathematical attempt to elucidate the evolutionary driving force in RNA virus evolution.

Intrinsic D614G and P681R/H mutations in SARS-CoV-2 VoCs Alpha, Delta, Omicron and viruses with D614G plus key signature mutations in spike protein alters fusogenicity and infectivity

The SARS-CoV-2 virus has been rapidly evolving over the time and the genetic variation has led to the generation of Variants of Concerns (VoC), which have shown increased fitness. These VoC viruses contain the key mutations in the spike protein which have allowed better survival and evasion of host defense mechanisms. The D614G mutation in the spike domain is found in the majority of VoC; additionally, the P681R/H mutation at the S1/S2 furin cleavage site junction is also found to be highly conserved in major VoCs; Alpha, Delta, Omicron, and its' current variants. The impact of these genetic alterations of the SARS-CoV-2 VoCs on the host cell entry, transmissibility, and infectivity has not been clearly identified. In our study, Delta and D614G + P681R synthetic double mutant pseudoviruses showed a significant increase in the cell entry, cell-to-cell fusion and infectivity. In contrast, the Omicron and P681H synthetic single mutant pseudoviruses showed TMPRSS2 independent cell entry, less fusion and infectivity as compared to Delta and D614G + P681R double mutants. Addition of exogenous trypsin further enhanced fusion in Delta viruses as compared to Omicron. Furthermore, Delta viruses showed susceptibility to both E64d and Camostat mesylate inhibitors suggesting, that the Delta virus could exploit both endosomal and TMPRSS2 dependent entry pathways as compared to the Omicron virus. Taken together, these results indicate that the D614G and P681R/H mutations in the spike protein are pivotal which might be favoring the VoC replication in different host compartments, and thus allowing a balance of mutation vs selection for better long-term adaptation.

Low neutralization of SARS-CoV-2 Omicron BA.2.75.2, BQ.1.1 and XBB.1 by parental mRNA vaccine or a BA.5 bivalent booster

The newly emerged severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron sublineages, including the BA.2-derived BA.2.75.2 and the BA.5-derived BQ.1.1 and XBB.1, have accumulated additional spike mutations that may affect vaccine effectiveness. Here we report neutralizing activities of three human serum panels collected from individuals 23-94 days after dose 4 of a parental mRNA vaccine; 14-32 days after a BA.5 bivalent booster from individuals with 2-4 previous doses of parental mRNA vaccine; or 14-32 days after a BA.5 bivalent booster from individuals with previous SARS-CoV-2 infection and 2-4 doses of parental mRNA vaccine. The results showed that a BA.5 bivalent booster elicited a high neutralizing titer against BA.4/5 measured at 14-32 days after boost; however, the BA.5 bivalent booster did not produce robust neutralization against the newly emerged BA.2.75.2, BQ.1.1 or XBB.1. Previous infection substantially enhanced the magnitude and breadth of BA.5 bivalent booster-elicited neutralization. Our data support a vaccine update strategy that future boosters should match newly emerged circulating SARS-CoV-2 variants.

Sequence analysis of SARS-CoV-2 Delta variant isolated from Makassar, South Sulawesi, Indonesia

Introduction

This study aimed to perform mutation and phylogenetic analyses of the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) Delta variants and analyze the characteristic signs and symptoms of patients infected with SARS-CoV-2 Delta variant originated from Makassar during the Delta outbreak.Methods: We collected samples from patients who were infected with coronavirus disease 2019 (COVID-19) between June and October 2021. We selected the Quantitative Reverse Transcription-Polymerase Chain Reaction (PCR)-positive samples with a cycle threshold value of <30 for whole genome sequencing. Total viral ribonucleic acid (RNA) was isolated from 34 PCR-positive nasopharyngeal swab samples, and whole genome sequencing was performed using the Oxford Nanopore GridlON sequencer. Phylogenetic and maximum clade credibility analyses were performed using the Bayesian Markov chain Monte Carlo method.

Results

It was found that 33 patients were infected with the SARS-CoV-2 Delta variant in this cohort study, among whom 63.6% (21) patients were female. According to the clinical data, 24 (72.7%), 7 (21.2%), and 2 (6.1%) patients had mild, moderate, and severe COVID-19 infections. Phylogenetic analysis based on the spike and RNA-dependent RNA polymerase (RdRp) genes showed that the collected samples were clustered in the main lineage of B.1.617.2 (Delta variant). The Delta variants had a high frequency of distinct mutations in the spike protein region, including T19R (94.12%), L452R (88.23%), T478K (91.17%), D614G (97%), P681R (97%), and D950 N (97%). Other unique mutations found in a smaller frequency in our samples were present in the N-terminal domain, including A27T (2.94%) and A222V (14.70%), and in the receptor-binding domain, including Q414K (5.88%), G446V (2.94%), and T470 N (2.94%).

Conclusion

This study revealed the unique mutations in the S protein region of Delta variants. T19R, L452R, T478K/T478R, D614G, P681R, and D950 N were the most common substitutions in Makassar's Delta variant.

Mathematical Modeling of SARS-CoV-2 Omicron Wave under Vaccination Effects

Over the course of the COVID-19 pandemic millions of deaths and hospitalizations have been reported. Different SARS-CoV-2 variants of concern have been recognized during this pandemic and some of these variants of concern have caused uncertainty and changes in the dynamics. The Omicron variant has caused a large amount of infected cases in the US and worldwide. The average number of deaths during the Omicron wave toll increased in comparison with previous SARS-CoV-2 waves. We studied the Omicron wave by using a highly nonlinear mathematical model for the COVID-19 pandemic. The novel model includes individuals who are vaccinated and asymptomatic, which influences the dynamics of SARS-CoV-2. Moreover, the model considers the waning of the immunity and efficacy of the vaccine against the Omicron strain. This study uses the facts that the Omicron strain has a higher transmissibility than the previous circulating SARS-CoV-2 strain but is less deadly. Preliminary studies have found that Omicron has a lower case fatality rate compared to previous circulating SARS-CoV-2 strains. The simulation results show that even if the Omicron strain is less deadly it might cause more deaths, hospitalizations and infections. We provide a variety of scenarios that help to obtain insight about the Omicron wave and its consequences. The proposed mathematical model, in conjunction with the simulations, provides an explanation for a large Omicron wave under various conditions related to vaccines and transmissibility. These results provide an awareness that new SARS-CoV-2 variants can cause more deaths even if their fatality rate is lower.

Reverse genetic systems of SARS-CoV-2 for antiviral research

Reverse genetic systems are widely used to engineer recombinant viruses with desired mutations. In response to the COVID-19 pandemic, four types of reverse genetic systems have been developed for SARS-CoV-2: (i) a full-length infectious clone that can be used to prepare recombinant SARS-CoV-2 at biosafety level 3 (BSL3), (ii) a trans-complementation system that can be used to produce single-round infectious SARS-CoV-2 at BSL2, (iii) an attenuated SARS-CoV-2 vaccine candidate (with deletions of viral accessory genes) that may be developed for veterinary use as well as for antiviral screening at BSL2, and (iv) replicon systems with deletions of viral structural genes that can be used at BSL2. Each of these genetic systems has its advantages and disadvantages that can be used to address different questions for basic and translational research. Due to the long genomic size and bacteria-toxic sequences of SARS-CoV-2, several experimental approaches have been established to rescue recombinant viruses and replicons, including (i) in vitro DNA ligation, (ii) bacterial artificial chromosome (BAC) system, (iii) yeast artificial chromosome (YAC) system, and (iv) circular polymerase extension reaction (CPER). This review summarizes the current status of SARS-CoV-2 genetic systems and their applications for studying viral replication, pathogenesis, vaccines, and therapeutics.

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.

SARS-CoV-2 evolution influences GBP and IFITM sensitivity

SARS-CoV-2 spike requires proteolytic processing for viral entry. A polybasic furin-cleavage site (FCS) in spike, and evolution toward an optimized FCS by dominant variants of concern (VOCs), are linked to enhanced infectivity and transmission. Here we show interferon-inducible restriction factors Guanylate-binding proteins (GBP) 2 and 5 interfere with furin-mediated spike cleavage and inhibit the infectivity of early-lineage isolates Wuhan-Hu-1 and VIC. By contrast, VOCs Alpha and Delta escape restriction by GBP2/5 that we map to the spike substitution D614G present in these VOCs. Despite inhibition of spike cleavage, these viruses remained sensitive to plasma membrane IFITM1, but not endosomal IFITM2 and 3, consistent with a preference for TMPRSS2-dependent plasma membrane entry. Strikingly, we find that Omicron is unique among VOCs, being sensitive to restriction factors GBP2/5, and also IFITM1, 2, and 3. Using chimeric spike mutants, we map the Omicron phenotype and show that the S1 domain determines Omicron's sensitivity to GBP2/5, whereas the S2' domain determines its sensitivity to endosomal IFITM2/3 and preferential use of TMPRSS2-independent entry. We propose that evolution of SARS-CoV-2 for the D614G substitution has allowed for escape from GBP restriction factors, but the selective pressures on Omicron for spike changes that mediate antibody escape, and altered tropism, have come at the expense of increased sensitivity to innate immune restriction factors that target virus entry.

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.

Unconventional secretion of unglycosylated ORF8 is critical for the cytokine storm during SARS-CoV-2 infection

Coronavirus disease 2019 is a respiratory infectious disease caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Evidence on the pathogenesis of SARS-CoV-2 is accumulating rapidly. In addition to structural proteins such as Spike and Envelope, the functional roles of non-structural and accessory proteins in regulating viral life cycle and host immune responses remain to be understood. Here, we show that open reading frame 8 (ORF8) acts as messenger for inter-cellular communication between alveolar epithelial cells and macrophages during SARS-CoV-2 infection. Mechanistically, ORF8 is a secretory protein that can be secreted by infected epithelial cells via both conventional and unconventional secretory pathways. Conventionally secreted ORF8 is glycosylated and loses the ability to recognize interleukin 17 receptor A of macrophages, possibly due to the steric hindrance imposed by N-glycosylation at Asn78. However, unconventionally secreted ORF8 does not undergo glycosylation without experiencing the ER-Golgi trafficking, thereby activating the downstream NF-κB signaling pathway and facilitating a burst of cytokine release. Furthermore, we show that ORF8 deletion in SARS-CoV-2 attenuates inflammation and yields less lung lesions in hamsters. Our data collectively highlights a role of ORF8 protein in the development of cytokine storms during SARS-CoV-2 infection.

Network analysis uncovers the communication structure of SARS-CoV-2 spike protein identifying sites for immunogen design

The COVID-19 pandemic, caused by the SARS-CoV-2 virus, has triggered myriad efforts to understand the structure and dynamics of this complex pathogen. The spike glycoprotein of SARS-CoV-2 is a significant target for immunogens as it is the means by which the virus enters human cells, while simultaneously sporting mutations responsible for immune escape. These functional and escape processes are regulated by complex molecular-level interactions. Our study presents quantitative insights on domain and residue contributions to allosteric communication, immune evasion, and local- and global-level control of functions through the derivation of a weighted graph representation from all-atom MD simulations. Focusing on the ancestral form and the D614G-variant, we provide evidence of the utility of our approach by guiding the selection of a mutation that alters the spike's stability. Taken together, the network approach serves as a valuable tool to evaluate communication "hot-spots" in proteins to guide design of stable immunogens.

Investigation on the interaction mechanism of different SARS-CoV-2 spike variants with hACE2: insights from molecular dynamics simulations

Since the COVID-19 pandemic caused by Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), SARS-CoV-2 has evolved by acquiring genomic mutations, resulting in the recent emergence of several SARS-CoV-2 variants with improved transmissibility and infectivity relative to the original strain. An underlying mechanism may be the increased ability of the mutants to bind the receptor proteins and infect the host cell. In this work, we implemented all-atom molecular dynamics (MD) simulations to study the binding and interaction of the receptor binding domain (RBD) of the SARS-CoV-2 spike protein singly (D614G), doubly (D614G + L452R and D614G + N501Y), triply (D614G + N501Y + E484K), and quadruply (D614G + N501Y + E484K + K417T) mutated variants with the human angiotensin-converting enzyme 2 (hACE2) receptor protein in the host cell. A combination of multiple analysis approaches elucidated the effects of mutations and the extent of molecular divergence from multiple perspectives, including the dynamic correlated motions, interaction patterns, dominant motions, free energy landscape, and charge distribution on the electrostatic potential surface between the hACE2 and all RBD variants. Moreover, free energy calculations using the MM/PBSA method evaluated the binding affinity between these RBD variants and hACE2. The results showed that the D614G + N501Y + E484K variant possessed the lowest free energy value (highest affinity) compared to the D614G + N501Y + E484K + K417T, D614G + L452R, D614G + N501Y, and D614G mutants. The residue-based energy decomposition also indicated that the energy contribution of residues at the mutation site to the total binding energy was highly variable. The interaction mechanisms between the different RBD variants and hACE2 elucidated in this study will provide some insights into the development of drugs targeting the new SARS-CoV-2 variants.

Mutational analysis of the spike protein of SARS-COV-2 isolates revealed atomistic features responsible for higher binding and infectivity

Introduction: The perpetual appearance of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-COV-2), and its new variants devastated the public health and social fabric around the world. Understanding the genomic patterns and connecting them to phenotypic attributes is of great interest to devise a treatment strategy to control this pandemic. Materials and Methods: In this regard, computational methods to understand the evolution, dynamics and mutational spectrum of SARS-CoV-2 and its new variants are significantly important. Thus, herein, we used computational methods to screen the genomes of SARS-CoV-2 isolated from Pakistan and connect them to the phenotypic attributes of spike protein; we used stability-function correlation methods, protein-protein docking, and molecular dynamics simulation. Results: Using the Global initiative on sharing all influenza data (GISAID) a total of 21 unique mutations were identified, among which five were reported as stabilizing while 16 were destabilizing revealed through mCSM, DynaMut 2.0, and I-Mutant servers. Protein-protein docking with Angiotensin-converting enzyme 2 (ACE2) and monoclonal antibody (4A8) revealed that mutation G446V in the receptor-binding domain; R102S and G181V in the N-terminal domain (NTD) significantly affected the binding and thus increased the infectivity. The interaction pattern also revealed significant variations in the hydrogen bonding, salt bridges and non-bonded contact networks. The structural-dynamic features of these mutations revealed the global dynamic trend and the finding energy calculation further established that the G446V mutation increases the binding affinity towards ACE2 while R102S and G181V help in evading the host immune response. The other mutations reported supplement these processes indirectly. The binding free energy results revealed that wild type-RBD has a TBE of -60.55 kcal/mol while G446V-RBD reported a TBE of -73.49 kcal/mol. On the other hand, wild type-NTD reported -67.77 kcal/mol of TBE, R102S-NTD reported -51.25 kcal/mol of TBE while G181V-NTD reported a TBE of -63.68 kcal/mol. Conclusions: In conclusion, the current findings revealed basis for higher infectivity and immune evasion associated with the aforementioned mutations and structure-based drug discovery against such variants.

Physical-Chemical Regulation of Membrane Receptors Dynamics in Viral Invasion and Immune Defense

Mechanical cues dynamically regulate membrane receptors functions to trigger various physiological and pathological processes from viral invasion to immune defense. These cues mainly include various types of dynamic mechanical forces and the spatial confinement of plasma membrane. However, the molecular mechanisms of how they couple with biochemical cues in regulating membrane receptors functions still remain mysterious. Here, we review recent advances in methodologies of single-molecule biomechanical techniques and in novel biomechanical regulatory mechanisms of critical ligand recognition of viral and immune receptors including SARS-CoV-2 spike protein, T cell receptor (TCR) and other co-stimulatory immune receptors. Furthermore, we provide our perspectives of the general principle of how force-dependent kinetics determine the dynamic functions of membrane receptors and of biomechanical-mechanism-driven SARS-CoV-2 neutralizing antibody design and TCR engineering for T-cell-based therapies.

In-Silico Lead Druggable Compounds Identification against SARS COVID-19 Main Protease Target from In-House, Chembridge and Zinc Databases by Structure-Based Virtual Screening, Molecular Docking and Molecular Dynamics Simulations

Pharmacological strategies to lower the viral load among patients suffering from severe diseases were researched in great detail during the SARS-CoV-2 outbreak. The viral protease M^pro (3CLpro) is necessary for viral replication and is among the main therapeutic targets proposed, thus far. To stop the pandemic from spreading, researchers are working to find more effective M^pro inhibitors against SARS-CoV-2. The 33.8 kDa M^pro protease of SARS-CoV-2, being a nonhuman homologue, has the possibility of being utilized as a therapeutic target against coronaviruses. To develop drug-like compounds capable of preventing the replication of SARS-main CoV-2's protease (M^pro), a computer-aided drug design (CADD) approach is extremely viable. Using MOE, structure-based virtual screening (SBVS) of in-house and commercial databases was carried out using SARS-CoV-2 proteins. The most promising hits obtained during virtual screening (VS) were put through molecular docking with the help of MOE. The virtual screening yielded 3/5 hits (in-house database) and 56/66 hits (commercial databases). Finally, 3/5 hits (in-house database), 3/5 hits (ZINC database), and 2/7 hits (ChemBridge database) were chosen as potent lead compounds using various scaffolds due to their considerable binding affinity with M^pro protein. The outcomes of SBVS were then validated using an analysis based on molecular dynamics simulation (MDS). The complexes' stability was tested using MDS and post-MDS. The most promising candidates were found to exhibit a high capacity for fitting into the protein-binding pocket and interacting with the catalytic dyad. At least one of the scaffolds selected will possibly prove useful for future research. However, further scientific confirmation in the form of preclinical and clinical research is required before implementation.

Differences in New Variant of Concern Replication at Physiological Temperatures In Vitro

Using multiple cell types and isolates of Delta and Omicron variants of SARS-CoV-2, we report differences in virus production, replication, and infectivity in vitro. Ancestral and Delta SARS-CoV-2 variant exhibit reduced virus production and replication at 34°C compared to 37°C while Omicron replication is balanced between temperatures.

Enhanced neutralization resistance of SARS-CoV-2 Omicron subvariants BQ.1, BQ.1.1, BA.4.6, BF.7, and BA.2.75.2

The continued evolution of SARS-CoV-2 has led to the emergence of several new Omicron subvariants, including BQ.1, BQ.1.1, BA.4.6, BF.7, and BA.2.75.2. Here, we examine the neutralization resistance of these subvariants against sera from 3-dose vaccinated healthcare workers, hospitalized BA.1-wave patients, and BA.4/5-wave patients. We found enhanced neutralization resistance in all new subvariants, especially in the BQ.1 and BQ.1.1 subvariants driven by N460K and K444T mutations, as well as the BA.2.75.2 subvariant driven largely by its F486S mutation. All Omicron subvariants maintained their weakened infectivity in Calu-3 cells, with the F486S mutation driving further diminished titer for the BA.2.75.2 subvariant. Molecular modeling revealed the mechanisms of antibody-mediated immune evasion by R346T, K444T, F486S, and D1199N mutations. Altogether, these findings shed light on the evolution of newly emerging SARS-CoV-2 Omicron subvariants.

Mucociliary transport deficiency and disease progression in Syrian hamsters with SARS-CoV-2 infection

Substantial clinical evidence supports the notion that ciliary function in the airways is important in COVID-19 pathogenesis. Although ciliary damage has been observed in both in vitro and in vivo models, the extent or nature of impairment of mucociliary transport (MCT) in in vivo models remains unknown. We hypothesize that SARS-CoV-2 infection results in MCT deficiency in the airways of golden Syrian hamsters that precedes pathological injury in lung parenchyma. Micro-optical coherence tomography was used to quantitate functional changes in the MCT apparatus. Both genomic and subgenomic viral RNA pathological and physiological changes were monitored in parallel. We show that SARS-CoV-2 infection caused a 67% decrease in MCT rate as early as 2 days postinfection (dpi) in hamsters, principally due to 79% diminished airway coverage of motile cilia. Correlating quantitation of physiological, virological, and pathological changes reveals steadily descending infection from the upper airways to lower airways to lung parenchyma within 7 dpi. Our results indicate that functional deficits of the MCT apparatus are a key aspect of COVID-19 pathogenesis, may extend viral retention, and could pose a risk factor for secondary infection. Clinically, monitoring abnormal ciliated cell function may indicate disease progression. Therapies directed toward the MCT apparatus deserve further investigation.

Vaccination with Span, an antigen guided by SARS-CoV-2 S protein evolution, protects against challenge with viral variants in mice

SARS-CoV-2 continues to accumulate mutations to evade immunity, leading to breakthrough infections after vaccination. How researchers can anticipate the evolutionary trajectory of the virus in advance in the design of next-generation vaccines requires investigation. Here, we performed a comprehensive study of 11,650,487 SARS-CoV-2 sequences, which revealed that the SARS-CoV-2 spike (S) protein evolved not randomly but into directional paths of either high infectivity plus low immune resistance or low infectivity plus high immune resistance. The viral infectivity and immune resistance of variants are generally incompatible, except for limited variants such as Beta and Kappa. The Omicron variant has the highest immune resistance but showed high infectivity in only one of the tested cell lines. To provide cross-clade immunity against variants that undergo diverse evolutionary pathways, we designed a new pan-vaccine antigen (S_pan). S_pan was designed by analyzing the homology of 2675 SARS-CoV-2 S protein sequences from the NCBI database before the Delta variant emerged. The refined S_pan protein harbors high-frequency residues at given positions that reflect cross-clade generality in sequence evolution. Compared with a prototype wild-type (S_wt) vaccine, which, when administered to mice, induced serum with decreased neutralization activity against emerging variants, S_pan vaccination of mice elicited broad immunity to a wide range of variants, including those that emerged after our design. Moreover, vaccinating mice with a heterologous S_pan booster conferred complete protection against lethal infection with the Omicron variant. Our results highlight the importance and feasibility of a universal vaccine to fight against SARS-CoV-2 antigenic drift.

Nanomechanical assay for ultrasensitive and rapid detection of SARS-CoV-2 based on peptide nucleic acid

The massive global spread of the COVID-19 pandemic makes the development of more effective and easily popularized assays critical. Here, we developed an ultrasensitive nanomechanical method based on microcantilever array and peptide nucleic acid (PNA) for the detection of severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2) RNA. The method has an extremely low detection limit of 0.1 fM (10⁵ copies/mL) for N-gene specific sequence (20 bp). Interestingly, it was further found that the detection limit of N gene (pharyngeal swab sample) was even lower, reaching 50 copies/mL. The large size of the N gene dramatically enhances the sensitivity of the nanomechanical sensor by up to three orders of magnitude. The detection limit of this amplification-free assay method is an order of magnitude lower than RT-PCR (500 copies/mL) that requires amplification. The non-specific signal in the assay is eliminated by the in-situ comparison of the array, reducing the false-positive misdiagnosis rate. The method is amplification-free and label-free, allowing for accurate diagnosis within 1 h. The strong specificity and ultra-sensitivity allow single base mutations in viruses to be distinguished even at very low concentrations. Also, the method remains sensitive to fM magnitude lung cancer marker (miRNA-155). Therefore, this ultrasensitive, amplification-free and inexpensive assay is expected to be used for the early diagnosis of COVID-19 patients and to be extended as a broad detection tool.

ELECTRONIC SUPPLEMENTARY MATERIAL

Supplementary material (experimental section, N gene sequences and all nucleic acid sequences used in the study, Figs. S1-S6, and Tables S1-S3) is available in the online version of this article at 10.1007/s12274-022-4333-3.

Wastewater pandemic preparedness: Toward an end-to-end pathogen monitoring program

Molecular analysis of public wastewater has great potential as a harbinger for community health and health threats. Long-used to monitor the presence of enteric viruses, in particular polio, recent successes of wastewater as a reliable lead indicator for trends in SARS-CoV-2 levels and hospital admissions has generated optimism and emerging evidence that similar science can be applied to other pathogens of pandemic potential (PPPs), especially respiratory viruses and their variants of concern (VOC). However, there are substantial challenges associated with implementation of this ideal, namely that multiple and distinct fields of inquiry must be bridged and coordinated. These include engineering, molecular sciences, temporal-geospatial analytics, epidemiology and medical, and governmental and public health messaging, all of which present their own caveats. Here, we outline a framework for an integrated, state-wide, end-to-end human pathogen monitoring program using wastewater to track viral PPPs.

Molecular characterisation and tracking of severe acute respiratory syndrome coronavirus 2 in Thailand, 2020-2022

The global COVID-19 pandemic, caused by the novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), was first detected in China in December 2019. To date, there have been approximately 3.4 million reported cases of COVID-19 and over 24,000 deaths in Thailand. In this study, we investigated the molecular characteristics and evolution of SARS-CoV-2 in Thailand from 2020 to 2022. Two hundred sixty-eight SARS-CoV-2 isolates, collected mostly in Bangkok from COVID-19 patients, were characterised by partial genome sequencing. Moreover, the viruses in 5,627 positive SARS-CoV-2 samples were identified as viral variants - B.1.1.7 (Alpha), B.1.617.2 (Delta), B.1.1.529 (Omicron/BA.1), or B.1.1.529 (Omicron/BA.2) - by multiplex real-time reverse transcription polymerase chain reaction (RT-PCR) assays. The results revealed that B.1.36.16 caused the predominant outbreak in the second wave (December 2020-January 2021), B.1.1.7 (Alpha) in the third wave (April-June 2021), B.1.617.2 (Delta) in the fourth wave (July-December 2021), and B.1.1.529 (Omicron) in the fifth wave (January-March 2022). The evolutionary rate of the viral genome was 2.60 × 10^-3 (95% highest posterior density [HPD], 1.72 × 10^-3 to 3.62 × 10^-3) nucleotide substitutions per site per year. Continued molecular surveillance of SARS-CoV-2 is crucial for monitoring emerging variants with the potential to cause new COVID-19 outbreaks.

Wastewater surveillance of SARS-CoV-2 genomic populations on a country-wide scale through targeted sequencing

SARS-CoV-2 surveillance of viral populations in wastewater samples is recognized as a useful tool for monitoring epidemic waves and boosting health preparedness. Next generation sequencing of viral RNA isolated from wastewater is a convenient and cost-effective strategy to understand the molecular epidemiology of SARS-CoV-2 and provide insights on the population dynamics of viral variants at the community level. However, in low- and middle-income countries, isolated groups have performed wastewater monitoring and data has not been extensively shared in the scientific community. Here we report the results of monitoring the co-circulation and abundance of variants of concern (VOCs) of SARS-CoV-2 in Uruguay, a small country in Latin America, between November 2020-July 2021 using wastewater surveillance. RNA isolated from wastewater was characterized by targeted sequencing of the Receptor Binding Domain region within the spike gene. Two computational approaches were used to track the viral variants. The results of the wastewater analysis showed the transition in the overall predominance of viral variants in wastewater from No-VOCs to successive VOCs, in agreement with clinical surveillance from sequencing of nasal swabs. The mutations K417T, E484K and N501Y, that characterize the Gamma VOC, were detected as early as December 2020, several weeks before the first clinical case was reported. Interestingly, a non-synonymous mutation described in the Delta VOC, L452R, was detected at a very low frequency since April 2021 when using a recently described sequence analysis tool (SAM Refiner). Wastewater NGS-based surveillance of SARS-CoV-2 is a reliable and complementary tool for monitoring the introduction and prevalence of VOCs at a community level allowing early public health decisions. This approach allows the tracking of symptomatic and asymptomatic individuals, who are generally under-reported in countries with limited clinical testing capacity. Our results suggests that wastewater-based epidemiology can contribute to improving public health responses in low- and middle-income countries.

Isolation and characterization of a SARS-CoV-2 variant with a Q677H mutation in the spike protein

We isolated 20 SARS-CoV-2 strains from positive clinical samples collected in Columbus, Ohio, and investigated the replication of one pair of isolates: a clade 20G strain and a variant of this strain carrying a Q677H mutation in the spike protein and six other amino acid mutations. The OSU.20G variant replicated to a higher peak infectious titer than the 20G base strain in Vero-E6 cells, but the titers were similar when both strains were grown in Calu-3 cells. These results suggest that the OSU.20G variant has increased replication fitness compared to the 20G base strain. This may have contributed to its emergence in December 2020-January 2021.

Semantic and Population Analysis of the Genetic Targets Related to COVID-19 and Its Association with Genes and Diseases

SARS-CoV-2 is a coronavirus responsible for one of the most serious, modern worldwide pandemics, with lasting and multifaceted effects. By late 2021, SARS-CoV-2 has infected more than 180 million people and has killed more than 3 million. The virus gains entrance to human cells through binding to ACE2 via its surface spike protein and causes a complex disease of the respiratory system, termed COVID-19. Vaccination efforts are being made to hinder the viral spread, and therapeutics are currently under development. Toward this goal, scientific attention is shifting toward variants and SNPs that affect factors of the disease such as susceptibility and severity. This genomic grammar, tightly related to the dark part of our genome, can be explored through the use of modern methods such as natural language processing. We present a semantic analysis of SARS-CoV-2-related publications, which yielded a repertoire of SNPs, genes, and disease ontologies. Population data from the 1000 Genomes Project were subsequently integrated into the pipeline. Data mining approaches of this scale have the potential to elucidate the complex interaction between COVID-19 pathogenesis and host genetic variation; the resulting knowledge can facilitate the management of high-risk groups and aid the efforts toward precision medicine.

Bivalent SARS-CoV-2 mRNA vaccines increase breadth of neutralization and protect against the BA.5 Omicron variant in mice

The emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants in the Omicron lineage has resulted in diminished Coronavirus Disease 2019 (COVID-19) vaccine efficacy and persistent transmission. In this study, we evaluated the immunogenicity and protective efficacy of two, recently authorized, bivalent COVID-19 vaccines that contain two mRNAs encoding Wuhan-1 and either BA.1 (mRNA-1273.214) or BA.4/5 (mRNA-1273.222) spike proteins. As a primary two-dose immunization series in mice, both bivalent vaccines induced greater neutralizing antibody responses against Omicron variants than the parental, monovalent mRNA-1273 vaccine. When administered to mice as a booster at 7 months after the primary vaccination series with mRNA-1273, the bivalent vaccines induced broadly neutralizing antibody responses. Whereas most anti-Omicron receptor binding domain antibodies in serum induced by mRNA-1273, mRNA-1273.214 and mRNA-1273.222 boosters cross-reacted with the antecedent Wuhan-1 spike antigen, the mRNA-1273.214 and mRNA-1273.222 bivalent vaccine boosters also induced unique BA.1-specific and BA.4/5-specific responses, respectively. Although boosting with parental or bivalent mRNA vaccines substantially improved protection against BA.5 compared to mice receiving two vaccine doses, the levels of infection, inflammation and pathology in the lung were lowest in animals administered the bivalent mRNA vaccines. Thus, boosting with bivalent Omicron-based mRNA-1273.214 or mRNA-1273.222 vaccines enhances immunogenicity and confers protection in mice against a currently circulating SARS-CoV-2 strain.

An antibody cocktail with broadened mutational resistance and effective protection against SARS-CoV-2

Neutralizing antibodies have been proven to be highly effective in treating mild and moderate COVID-19 patients, but continuous emergence of SARS-CoV-2 variants poses significant challenges. Antibody cocktail treatments reduce the risk of escape mutants and resistance. In this study, a new cocktail composed of two highly potent neutralizing antibodies (HB27 and H89Y) was developed, whose binding epitope is different from those cocktails that received emergency use authorization. This cocktail showed more potent and balanced neutralizing activities (IC₅₀ 0.9-11.3 ng mL^-1) against a broad spectrum of SARS-CoV-2 variants over individual HB27 or H89Y antibodies. Furthermore, the cocktail conferred more effective protection against the SARS-CoV-2 Beta variant in an aged murine model than monotherapy. It was shown to prevent SARS-CoV-2 mutational escape in vitro and effectively neutralize 61 types of pseudoviruses harbouring single amino acid mutation originated from variants and escape strains of Bamlanivimab, Casirivimab and Imdevimab with IC₅₀ of 0.6-65 ng mL^-1. Despite its breadth of variant neutralization, the HB27+H89Y combo and EUA cocktails lost their potencies against Omicron variant. Our results provide important insights that new antibody cocktails covering different epitopes are valuable tools to counter virus mutation and escape, highlighting the need to search for more conserved epitopes to combat Omicron.

Origin and evolution of SARS-CoV-2

SARS-CoV-2 is a novel coronavirus that emerged in China at the end of 2019 causing the severe disease known as coronavirus disease 2019 (COVID-19). SARS-CoV-2, as to the previously highly pathogenic human coronaviruses named SARS-CoV, the etiological agent of severe acute respiratory syndrome (SARS), has a zoonotic origin, although SARS-CoV-2 precise chain of animal-to-human transmission remains undefined. Unlike the 2002-2003 pandemic caused by SARS-CoV whose extinction from the human population was achieved in eight months, SARS-CoV-2 has been spreading globally in an immunologically naïve population in an unprecedented manner. The efficient infection and replication of SARS-CoV-2 has resulted in the emergence of viral variants that have become predominant posing concerns about their containment as they are more infectious with variable pathogenicity in respect to the original virus. Although vaccine availability is limiting severe disease and death caused by SARS-CoV-2 infection, its extinction is far to be close and predictable. In this regard, the emersion of the Omicron viral variant in November 2021 was characterized by humoral immune escape and it has reinforced the importance of the global monitoring of SARS-CoV-2 evolution. Given the importance of the SARS-CoV-2 zoonotic origin, it will also be crucial to monitor the animal-human interface to be better prepared to cope with future infections of pandemic potential.

Reduced neutralization against Delta, Gamma, Mu, and Omicron BA.1 variants of SARS-CoV-2 from previous non-Omicron infection

The understanding of the host immune response to SARS-CoV-2 variants of concern is critical for improving diagnostics, therapy development, and vaccines. Here, we analyzed the level of neutralizing antibodies against SARS-CoV-2 D614G, Delta, Gamma, Mu, and Omicron variants in D614G infected healthcare workers during a follow-up up to 6 months after recovery. We followed up 76 patients: 60.5% were women and 39.5% men. The 96.1% and 3.9% were symptomatic and asymptomatic, respectively. The most frequent symptoms were headache, myalgia, and cough. The 65.8%, 65.8%, and 92.1% of the infected individuals were positive for neutralizing antibodies against D614G variant at 2, 4, and 6 months of follow-up, respectively. The 26.3%, 48.7% and 65.8% of patients neutralized Delta variant, 19.7%, 32.9% and 52.6% of patients neutralized Gamma, 7.9%, 19.7% and 44.7% of patients neutralized Mu, and 4.0%, 9.2% and 15.8% of patients neutralized Omicron. Low neutralization against Gamma and Mu variants was observed during the follow-up, and very low against the Omicron variant was detected during the same period. The median of neutralizing antibody titers against D614G and Delta variants increased significantly during the follow-up. An association was observed between the levels of neutralizing antibodies against D614G and Delta variants and the severity of the disease. Our results suggest an immune escape from neutralizing antibodies with the Omicron variant because of the many mutations localized in the S protein.

The COVID-19 Pandemic: SARS-CoV-2 Structure, Infection, Transmission, Symptomology, and Variants of Concern

Since it was first detected in December 2019, the COVID-19 pandemic has spread across the world and affected virtually every country and territory. The pathogen driving this pandemic is SARS-CoV-2, a positive-sense single-stranded RNA virus which is primarily transmissible though the air and can cause mild to severe respiratory infections in humans. Within the first year of the pandemic, the situation worsened with the emergence of several SARS-CoV-2 variants. Some of these were observed to be more virulent with varying capacities to escape the existing vaccines and were, therefore, denoted as variants of concern. This chapter provides a general overview of the course of the COVID-19 pandemic up to April 2022 with a focus on the structure, infection, transmission, and symptomology of the SARS-CoV-2 virus. The main objectives were to investigate the effects of the variants of concern on the trajectory of the virus and to highlight a potential pathway for coping with the current and future pandemics.

SARS-CoV-2 infection dysregulates the expression of clinically relevant drug metabolizing enzymes in Vero E6 cells and membrane transporters in human lung tissues

SARS-CoV-2-mediated interactions with drug metabolizing enzymes and membrane transporters (DMETs) in different tissues, especially lung, the main affected organ may limit the clinical efficacy and safety profile of promising COVID-19 drugs. Herein, we investigated whether SARS-CoV-2 infection could dysregulate the expression of 25 clinically relevant DMETs in Vero E6 cells and postmortem lung tissues from COVID-19 patients. Also, we assessed the role of 2 inflammatory and 4 regulatory proteins in modulating the dysregulation of DMETs in human lung tissues. We showed for the first time that SARS-CoV-2 infection dysregulates CYP3A4 and UGT1A1 at the mRNA level, as well as P-gp and MRP1 at the protein level, in Vero E6 cells and postmortem human lung tissues, respectively. We observed that at the cellular level, DMETs could potentially be dysregulated by SARS-CoV-2-associated inflammatory response and lung injury. We uncovered the pulmonary cellular localization of CYP1A2, CYP2C8, CYP2C9, and CYP2D6, as well as ENT1 and ENT2 in human lung tissues, and observed that the presence of inflammatory cells is the major driving force for the discrepancy in the localization of DMETs between COVID-19 and control human lung tissues. Because alveolar epithelial cells and lymphocytes are both sites of SARS-CoV-2 infection and localization of DMETs, we recommend further investigation of the pulmonary pharmacokinetic profile of current COVID-19 drug dosing regimen to improve clinical outcomes.

Attenuation and Degeneration of SARS-CoV-2 Despite Adaptive Evolution

The evolution of severe acute respiratory syndrome coronavirus 2 (SARS‑CoV‑2) has followed similar trends as other RNA viruses, such as human immunodeficiency virus type 1 and the influenza A virus. Rapid initial diversification was followed by strong competition and a rapid succession of dominant variants. Host-initiated RNA editing has been the primary mechanism for introducing mutations. A significant number of mutations detrimental to viral replication have been quickly purged. Fixed mutations are mostly diversifying mutations selected for host adaptation and immune evasion, with the latter accounting for the majority of the mutations. However, immune evasion often comes at the cost of functionality, and thus, optimal functionality is still far from being accomplished. Instead, selection for antibody-escaping variants and accumulation of near-neutral mutations have led to suboptimal codon usage and reduced replicative capacity, as demonstrated in non-respiratory cell lines. Beneficial adaptation of the virus includes reduced infectivity in lung tissues and increased tropism for the upper airway, resulting in shorter incubation periods, milder diseases, and more efficient transmission between people.

A hypothesis on designing strategy of effective RdRp inhibitors for the treatment of SARS-CoV-2

Vaccines are used as one of the major weapons for the eradication of pandemic. However, the rise of different variants of the SARS-CoV-2 virus is creating doubts regarding the end of the pandemic. Hence, there is an urgent need to develop more drug candidates which can be useful for the treatment of COVID-19. In the present research for the scientific hypothesis, emphasis was given on the direct antiviral therapy available for the treatment of COVID-19. In lieu of this, the available molecular targets which include Severe Acute Respiratory Syndrome Chymotrypsin-like Protease (SARS-3CLpro), Papain-Like Cysteine Protease (PLpro), and RNA-Dependent RNA Polymerase (RdRp) were explored. As per the current scientific reports and literature, among all the available molecular targets, RNA-Dependent RNA Polymerase (RdRp) was found to be a crucial molecular target for the treatment of COVID-19. Most of the inhibitors which are reported against this target consisted of the free amine group and carbonyl group which might be playing an important role in the binding interaction with the RdRp protein. Among all the reported RdRp inhibitors, remdesivir, favipiravir, and molnupiravir were found to be the most promising drugs against COVID-19. Overall, the structural features of this RNA-Dependent RNA Polymerase (RdRp) inhibitors proved the importance of pyrrolo-triazine and pyrimidine scaffolds. Previous computational models of these drug molecules indicated that substitution with the polar functional group, hydrogen bond donor, and electronegative atoms on these scaffolds may increase the activity against the RdRp protein. Hence, in line with the proposed hypothesis, in the present research work for the evaluation of the hypothesis, new molecules were designed from the pyrrolo-triazine and pyrimidine scaffolds. Further, molecular docking and MD simulation studies were performed with these designed molecules. All these designed molecules (DM-1, DM-2, and DM-3) showed the results as per the proposed hypothesis. Among all the designed molecules, DM-1 showed promising results against the RdRp protein of SARS-CoV-2. In the future, these structural features can be used for the development of new RdRp inhibitors with improved activity. Also, in the future lead compound DM-1 can be explored against the RdRp protein for the treatment of COVID-19.

RNA structure-altering mutations underlying positive selection on Spike protein reveal novel putative signatures to trace crossing host-species barriers in Betacoronavirus

Similar to other RNA viruses, the emergence of Betacoronavirus relies on cross-species viral transmission, which requires careful health surveillance monitoring of protein-coding information as well as genome-wide analysis. Although the evolutionary jump from natural reservoirs to humans may be mainly traced-back by studying the effect that hotspot mutations have on viral proteins, it is largely unexplored if other impacts might emerge on the structured RNA genome of Betacoronavirus. In this survey, the protein-coding and viral genome architecture were simultaneously studied to uncover novel insights into cross-species horizontal transmission events. We analysed 1,252,952 viral genomes of SARS-CoV, MERS-CoV, and SARS-CoV-2 distributed across the world in bats, intermediate animals, and humans to build a new landscape of changes in the RNA viral genome. Phylogenetic analyses suggest that bat viruses are the most closely related to the time of most recent common ancestor of Betacoronavirus, and missense mutations in viral proteins, mainly in the S protein S1 subunit: SARS-CoV (G > T; A577S); MERS-CoV (C > T; S746R and C > T; N762A); and SARS-CoV-2 (A > G; D614G) appear to have driven viral diversification. We also found that codon sites under positive selection on S protein overlap with non-compensatory mutations that disrupt secondary RNA structures in the RNA genome complement. These findings provide pivotal factors that might be underlying the eventual jumping the species barrier from bats to intermediate hosts. Lastly, we discovered that nearly half of the Betacoronavirus genomes carry highly conserved RNA structures, and more than 90% of these RNA structures show negative selection signals, suggesting essential functions in the biology of Betacoronavirus that have not been investigated to date. Further research is needed on negatively selected RNA structures to scan for emerging functions like the potential of coding virus-derived small RNAs and to develop new candidate antiviral therapeutic strategies.

Mutational profile confers increased stability of SARS-CoV-2 spike protein in Brazilian isolates

Spike (S) protein has been recognized as a promising molecular target for diagnostic, vaccines and antiviral drugs development for COVID-19. In this study, we analyzed the most predominant mutations in the S protein of Brazilian isolates and predicted the effect of these amino acid alterations to protein conformation. A total of 25,924 sequences were obtained from GISAID for five regions of Brazilian territory (Midwest, North, Northeast, South, and Southeast), according to exclusion criteria. Most of the SARS-CoV-2 isolates belongs to the G clade and showed a large occurrence of D614G, N501Y and L18F substitutions. Prediction effects of these amino acid substitutions on the structure dynamics of the spike protein indicated a positive ΔΔG values and negative ΔΔS_Vib in most cases which is associated to structural stabilization and flexibility reduction of the S protein. Mutations E484K, N501Y and K417N belong to several SARS-CoV-2 variants of concern such as Alpha, Beta, Gamma and Delta, and showed high incidence among Brazilian isolates. These mutations have been described to increase RBD affinity to ACE-2 host and abolishment of RBD affinity to potent neutralizing ant-RBD. The increase in rates of infection and reinfection requires continuous genomic surveillance studies in order to characterize emerging mutations and monitor vaccine efficacy, and thus consideration structural data and dynamics in the observed phenotypes.Communicated by Ramaswamy H. Sarma.