Cryo-EM structure of a SARS-CoV-2 omicron spike protein ectodomain

Infectivity and structure of SARS-CoV-2 after hydrogen peroxide treatment

Hydrogen peroxide (H2O2) exhibits broad-spectrum antiviral activity and is commonly used as an over-the-counter disinfecting agent. However, its potential activities against SARS-CoV-2 have not been systematically evaluated, and mechanisms of action are not well understood. In this study, we investigate H2O2's antiviral activity against SARS-CoV-2 infection and its impact on the virion's structural integrity as compared to the commonly used fixative agent paraformaldehyde (PFA). We show that H2O2 rapidly and directly inactivates SARS-CoV-2 with a half-maximal inhibitory concentration (IC50) of 0.0015%. Cryogenic electron tomography (cryo-ET) with subtomogram averaging reveals that treatment with PFA induced the viral trimeric spike protein (S) to adopt a post-fusion conformation, and treatment of viral particles with H2O2 locked S in its pre-fusion conformation. Therefore, H2O2 treatment likely has induced modifications, such as oxidation of cysteine residues within the S subunits of the spike trimer that locked them in their pre-fusion conformation. Locking of the meta-stable pre-fusion trimer prevents its transition to the post-fusion conformation, a process essential for viral fusion with host cells and entry into host cells. Together, our cellular, biochemical, and structural studies established that hydrogen peroxide can inactivate SARS-CoV-2 in tissue culture and uncovered its underlying molecular mechanism.IMPORTANCEHydrogen peroxide (H2O2) is the commonly used, over-the-counter antiseptic solution available in pharmacies, but its effect against the SARS-CoV-2 virus has not been evaluated systematically. In this study, we show that H2O2 inactivates the SARS-CoV-2 infectivity and establish the effective concentration of this activity. Cryogenic electron tomography and sub-tomogram averaging reveal a detailed structural understanding of how H2O2 affects the SARS-CoV-2 spike in comparison with that of the commonly used fixative PFA under identical conditions. We found that PFA promoted a post-fusion conformation of the viral spike protein, while H2O2 could potentially lock the spike in its pre-fusion state. Our findings not only substantiate the disinfectant efficacy of H2O2 as a potent agent against SARS-CoV-2 but also lay the groundwork for future investigations into targeted antiviral therapies that may leverage the virus' structural susceptibilities. In addition, this study may have significant implications for developing new antiviral strategies and improving existing disinfection protocols.

Demonstration and structural basis of a therapeutic DNA aptamer for SARS-CoV-2 spike protein detection

At the onset of the COVID-19 pandemic, the absence of a rapid and highly specific diagnostic method for the SARS-CoV-2 virus led to significant delays in detection, adversely affecting clinical outcomes. This shortfall highlights the urgent need for adaptable, scalable, and reusable diagnostic technologies to improve future pandemic responses. To address this challenge, we developed a renewable electrochemical impedance biosensor device employing a synthetic nucleotide-based therapeutic aptamer (termed 'aptasensor') targeting the SARS-CoV-2 spike (S) protein receptor-binding domain (RBD). We demonstrate that our aptasensor can detect the Omicron BA.2 S protein within one hour and possesses concentration-dependent sensitivity at biologically relevant levels. Notably, the aptasensor is reusable after regeneration by a simple pH 2 buffer treatment. Aptamer binding to the S protein was confirmed by immunogold labeling and visualization by negative-stain electron microscopy. We used cryogenic electron microscopy (cryo-EM) to resolve high-resolution maps of the S protein in both the open and closed conformations and characterized aptamer binding to the up RBD in the open conformation. Taken together, these results establish the versatility and scalability of aptamer-based biosensors, presenting them as a potential transformative diagnostic platform for emerging pathogens. This combination of rapid detection, specificity, and renewable capabilities in a single diagnostic solution marks a significant advance in pandemic preparedness.

Designing and bioengineering of CDRs with higher affinity against receptor-binding domain (RBD) of SARS-CoV-2 Omicron variant

The emergence of the SARS-CoV-2 Omicron variant highlights the need for innovative strategies to address evolving viral threats. This study bioengineered three nanobodies H11-H4, C5, and H3 originally targeting the Wuhan RBD, to bind more effectively to the Omicron RBD. A structure-based in silico affinity maturation pipeline was developed to enhance their binding affinities. The pipeline consists of three key steps: high-throughput in silico mutagenesis of complementarity-determining regions (CDRs), protein-protein docking for screening, and molecular dynamics (MD) simulations for assessment of the complex stability. A total of 741, 551, and 684 mutations were introduced in H11-H4, C5, and H3 nanobodies, respectively. Protein-protein docking and MD simulations shortlisted high-affinity mutants for H11-H4(6), C5(5), and H3(6). Further, recombinant production of H11-H4 mutants and Omicron RBD enabled experimental validation through Isothermal Titration Calorimetry (ITC). The H11-H4 mutants R27E, S57D, S107K, D108W, and A110I exhibited improved binding affinities with dissociation constant (KD) values ranging from ~8.8 to ~27 μM, compared to the H11-H4 nanobody KD of ~32 μM, representing a three-fold enhancement. This study demonstrates the potential of the developed in silico affinity maturation pipeline as a rapid, cost-effective method for repurposing nanobodies, aiding the development of robust prophylactic strategies against evolving SARS-CoV-2 variants and other pathogens.

Visualizing the virus world inside the cell by cryo-electron tomography

Structural studies on purified virus have revealed intricate architectures, but there is little structural information on how viruses interact with host cells in situ. Cryo-focused ion beam (FIB) milling and cryo-electron tomography (cryo-ET) have emerged as revolutionary tools in structural biology to visualize the dynamic conformational of viral particles and their interactions with host factors within infected cells. Here, we review the state-of-the-art cryo-ET technique for in situ viral structure studies and highlight exemplary studies that showcase the remarkable capabilities of cryo-ET in capturing the dynamic virus-host interaction, advancing our understanding of viral infection and pathogenesis.

Preclinical characterization of the Omicron XBB.1.5-adapted BNT162b2 COVID-19 vaccine

As SARS-CoV-2 evolves, increasing in potential for greater transmissibility and immune escape, updated vaccines are needed to boost adaptive immunity to protect against COVID-19 caused by circulating strains. Here, we report features of the monovalent Omicron XBB.1.5-adapted BNT162b2 vaccine, which contains XBB.1.5-specific sequence changes, relative to the original BNT162b2 backbone, in the encoded prefusion-stabilized SARS-CoV-2 spike protein (S(P2)). Biophysical characterization of Omicron XBB.1.5 S(P2) demonstrated that it maintains a prefusion conformation and adopts a flexible, predominantly open, state, with high affinity for the human ACE-2 receptor. When administered as a 4th dose in BNT162b2-experienced mice, the monovalent Omicron XBB.1.5 vaccine elicited substantially higher serum neutralizing titers against pseudotyped viruses of Omicron XBB.1.5, XBB.1.16, XBB.1.16.1, XBB.2.3, EG.5.1 and HV.1 sublineages and phylogenetically distant BA.2.86 lineage than the bivalent Wild Type + Omicron BA.4/5 vaccine. Similar trends were observed against Omicron XBB sublineage pseudoviruses when the vaccine was administered as a 2-dose series in naive mice. Strong S-specific Th1 CD4+ and IFNγ+ CD8+ T cell responses were also observed. These findings, together with real world performance of the XBB.1.5-adapted vaccine, suggest that preclinical data for the monovalent Omicron XBB.1.5-adapted BNT162b2 was predictive of protective immunity against dominant SARS-CoV-2 strains.

Conformational dynamics of SARS-CoV-2 Omicron spike trimers during fusion activation at single molecule resolution

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron entry involves spike (S) glycoprotein-mediated fusion of viral and late endosomal membranes. Here, using single-molecule Förster resonance energy transfer (sm-FRET) imaging and biochemical measurements, we directly visualized conformational changes of individual spike trimers on the surface of SARS-CoV-2 Omicron pseudovirions during fusion activation. We observed that the S2 domain of the Omicron spike is a dynamic fusion machine. S2 reversibly interchanges between the pre-fusion conformation and two previously undescribed intermediate conformations. Acidic pH shifts the conformational equilibrium of S2 toward an intermediate conformation and promotes the membrane hemi-fusion reaction. Moreover, we captured conformational reversibility in the S2 domain, which suggests that spike can protect itself from pre-triggering. Furthermore, we determined that Ca2+ directly promotes the S2 conformational change from an intermediate conformation to post-fusion conformation. In the presence of a target membrane, low pH and Ca2+ stimulate the irreversible transition to S2 post-fusion state and promote membrane fusion.

Discovery of Nanosota-9 as anti-Omicron nanobody therapeutic candidate

Omicron subvariants of SARS-CoV-2 continue to pose a significant global health threat. Nanobodies, single-domain antibodies derived from camelids, are promising therapeutic tools against pandemic viruses due to their favorable properties. In this study, we identified a novel nanobody, Nanosota-9, which demonstrates high potency against a wide range of Omicron subvariants both in vitro and in a mouse model. Cryo-EM data revealed that Nanosota-9 neutralizes Omicron through a unique mechanism: two Nanosota-9 molecules crosslink two receptor-binding domains (RBDs) of the trimeric Omicron spike protein, preventing the RBDs from binding to the ACE2 receptor. This mechanism explains its strong anti-Omicron potency. Additionally, the Nanosota-9 binding epitopes on the spike protein are relatively conserved among Omicron subvariants, contributing to its broad anti-Omicron spectrum. Combined with our recently developed structure-guided in vitro evolution approach for nanobodies, Nanosota-9 has the potential to serve as the foundation for a superior anti-Omicron therapeutic.

Molecular and structural insights into SARS-CoV-2 evolution: from BA.2 to XBB subvariants

Due to the incessant emergence of various SARS-CoV-2 variants with enhanced fitness in the human population, controlling the COVID-19 pandemic has been challenging. Understanding how the virus enhances its fitness during a pandemic could offer valuable insights for more effective control of viral epidemics. In this manuscript, we review the evolution of SARS-CoV-2 from early 2022 to the end of 2023-from Omicron BA.2 to XBB descendants. Focusing on viral evolution during this period, we provide concrete examples that SARS-CoV-2 has increased its fitness by enhancing several functions of the spike (S) protein, including its binding affinity to the ACE2 receptor and its ability to evade humoral immunity. Furthermore, we explore how specific mutations modify these functions of the S protein through structural alterations. This review provides evolutionary, molecular, and structural insights into how SARS-CoV-2 has increased its fitness and repeatedly caused epidemic surges during the pandemic.

Dual-role epitope on SARS-CoV-2 spike enhances and neutralizes viral entry across different variants

Grasping the roles of epitopes in viral glycoproteins is essential for unraveling the structure and function of these proteins. Up to now, all identified epitopes have been found to either neutralize, have no effect on, or enhance viral entry into cells. Here, we used nanobodies (single-domain antibodies) as probes to investigate a unique epitope on the SARS-CoV-2 spike protein, located outside the protein's receptor-binding domain. Nanobody binding to this epitope enhances the cell entry of prototypic SARS-CoV-2, while neutralizing the cell entry of SARS-CoV-2 Omicron variant. Moreover, nanobody binding to this epitope promotes both receptor binding activity and post-attachment activity of prototypic spike, explaining the enhanced viral entry. The opposite occurs with Omicron spike, explaining the neutralized viral entry. This study reveals a unique epitope that can both enhance and neutralize viral entry across distinct viral variants, suggesting that epitopes may vary their roles depending on the viral context. Consequently, antibody therapies should be assessed across different viral variants to confirm their efficacy and safety.

Universal subunit vaccine protects against multiple SARS-CoV-2 variants and SARS-CoV

Although Omicron RBD of SARS-CoV-2 accumulates many mutations, the backbone region (truncated RBD) of spike protein is highly conserved. Here, we designed several subunit vaccines by keeping the conserved spike backbone region of SARS-CoV-2 Omicron BA.1 subvariant (S-6P-no-RBD), or inserting the RBD of Delta variant (S-6P-Delta-RBD), Omicron (BA.5) variant (S-6P-BA5-RBD), or ancestral SARS-CoV-2 (S-6P-WT-RBD) to the above backbone construct, and evaluated their ability to induce immune responses and cross-protective efficacy against various SARS-CoV-2 variants and SARS-CoV. Among the four subunit vaccines, S-6P-Delta-RBD protein elicited broad and potent neutralizing antibodies against all SARS-CoV-2 variants tested, including Alpha, Beta, Gamma, and Delta variants, the BA.1, BA.2, BA.2.75, BA.4.6, and BA.5 Omicron subvariants, and the ancestral strain of SARS-CoV-2. This vaccine prevented infection and replication of SARS-CoV-2 Omicron, and completely protected immunized mice against lethal challenge with the SARS-CoV-2 Delta variant and SARS-CoV. Sera from S-6P-Delta-RBD-immunized mice protected naive mice against challenge with the Delta variant, with significantly reduced viral titers and without pathological effects. Protection correlated positively with the serum neutralizing antibody titer. Overall, the designed vaccine has potential for development as a universal COVID-19 vaccine and/or a pan-sarbecovirus subunit vaccine that will prevent current and future outbreaks caused by SARS-CoV-2 variants and SARS-related CoVs.

Single-molecule imaging reveals allosteric stimulation of SARS-CoV-2 spike receptor binding domain by host sialic acid

Conformational dynamics of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike glycoprotein (S) mediate exposure of the binding site for the cellular receptor, angiotensin-converting enzyme 2 (ACE2). The N-terminal domain (NTD) of S binds terminal sialic acid (SA) moieties on the cell surface, but the functional role of this interaction in virus entry is unknown. Here, we report that NTD-SA interaction enhances both S-mediated virus attachment and ACE2 binding. Through single-molecule Förster resonance energy transfer imaging of individual S trimers, we demonstrate that SA binding to the NTD allosterically shifts the S conformational equilibrium, favoring enhanced exposure of the ACE2-binding site. Antibodies that target the NTD block SA binding, which contributes to their mechanism of neutralization. These findings inform on mechanisms of S activation at the cell surface.

Enhancing the understanding of SARS-CoV-2 protein with structure and detection methods: An integrative review

As the rapid and accurate screening of infectious diseases can provide meaningful information for outbreak prevention and control, as well as owing to the existing limitations of the polymerase chain reaction (PCR), it is imperative to have new and validated detection techniques for SARS-CoV-2. Therefore, the rationale for outlining the techniques used to detect SARS-CoV-2 proteins and performing a comprehensive comparison to serve as a practical benchmark for future identification of similar viral proteins is clear. This review highlights the urgent need to strengthen pandemic preparedness by emphasizing the importance of integrated measures. These include improved tools for pathogen characterization, optimized societal precautions, the establishment of early warning systems, and the deployment of highly sensitive diagnostics for effective surveillance, triage, and resource management. Additionally, with an improved understanding of the virus' protein structure, considerable advances in targeted detection, treatment, and prevention strategies are expected to greatly improve our ability to respond to future outbreaks.

Therapeutic nanobodies against SARS-CoV-2 and other pathogenic human coronaviruses

Nanobodies, single-domain antibodies derived from variable domain of camelid or shark heavy-chain antibodies, have unique properties with small size, strong binding affinity, easy construction in versatile formats, high neutralizing activity, protective efficacy, and manufactural capacity on a large-scale. Nanobodies have been arisen as an effective research tool for development of nanobiotechnologies with a variety of applications. Three highly pathogenic coronaviruses (CoVs), SARS-CoV-2, SARS-CoV, and MERS-CoV, have caused serious outbreaks or a global pandemic, and continue to post a threat to public health worldwide. The viral spike (S) protein and its cognate receptor-binding domain (RBD), which initiate viral entry and play a critical role in virus pathogenesis, are important therapeutic targets. This review describes pathogenic human CoVs, including viral structures and proteins, and S protein-mediated viral entry process. It also summarizes recent advances in development of nanobodies targeting these CoVs, focusing on those targeting the S protein and RBD. Finally, we discuss potential strategies to improve the efficacy of nanobodies against emerging SARS-CoV-2 variants and other CoVs with pandemic potential. It will provide important information for rational design and evaluation of therapeutic agents against emerging and reemerging pathogens.

In Silico Analyses Indicate a Lower Potency for Dimerization of TLR4/MD-2 as the Reason for the Lower Pathogenicity of Omicron Compared to Wild-Type Virus and Earlier SARS-CoV-2 Variants

The SARS-CoV-2 Omicron variants have replaced all earlier variants, due to increased infectivity and effective evasion from infection- and vaccination-induced neutralizing antibodies. Compared to earlier variants of concern (VoCs), the Omicron variants show high TMPRSS2-independent replication in the upper airway organs, but lower replication in the lungs and lower mortality rates. The shift in cellular tropism and towards lower pathogenicity of Omicron was hypothesized to correlate with a lower toll-like receptor (TLR) activation, although the underlying molecular mechanisms remained undefined. In silico analyses presented here indicate that the Omicron spike protein has a lower potency to induce dimerization of TLR4/MD-2 compared to wild type virus despite a comparable binding activity to TLR4. A model illustrating the molecular consequences of the different potencies of the Omicron spike protein vs. wild-type spike protein for TLR4 activation is presented. Further analyses indicate a clear tendency for decreasing TLR4 dimerization potential during SARS-CoV-2 evolution via Alpha to Gamma to Delta to Omicron variants.

Single particle cryo-electron microscopy with an enhanced 200 kV cryo-TEM configuration achieves near-atomic resolution

Single particle cryogenic electron microscopy (cryo-EM) as a structural biology methodology has become increasingly attractive and accessible to investigators in both academia and industry as this ever-advancing technology enables successful structural determination of a wide range of protein and nucleic acid targets. Although data for many high resolution cryo-EM structures are still obtained using a 300 kV cryogenic transmission electron microscope (cryo-TEM), a modern 200 kV cryo-TEM equipped with an advanced direct electron detector and energy filter is a cost-effective choice for most single particle applications, routinely achieving sub 3 angstrom (Å) resolution. Here, we systematically evaluate performance of one such high-end configuration - a 200 kV Glacios microscope coupled with a Falcon 4 direct electron detector and Selectris energy filter (Glacios-F4-S). First, we evaluated data quality on the standard benchmarking sample, rabbit muscle aldolase, using three of the most frequently used cryo-EM data collection software: SerialEM, Leginon and EPU, and found that - despite sample heterogeneity - all final reconstructions yield same overall resolutions of 2.6 Å and map quality when using either of the three software. Furthermore, comparison between Glacios-F4-S and a 300 kV cryo-TEM (Titan Krios with Falcon 4) revealed nominal resolution differences in overall reconstructions of a reconstituted human nucleosome core particle, achieving 2.8 and 2.5 Å, respectively. Finally, we performed comparative data analysis on the human RAD51 paralog complex, BCDX2, a four-protein complex of approximately 150 kilodaltons, and found that a small dataset (≤1,000 micrographs) was sufficient to generate a 3.3 Å reconstruction, with sufficient detail to resolve co-bound ligands, AMP-PNP and Mg +2 . In summary, this study provides evidence that the Glacios-F4-S operates equally well with all standard data collection software, and is sufficient to obtain high resolution structural information of novel macromolecular complexes, readily acquiring single particle data rivaling that of 300 kV cryo-TEMs.

Effects of Strain Differences, Humidity Changes, and Saliva Contamination on the Inactivation of SARS-CoV-2 by Ion Irradiation

One of the methods to inactivate viruses is to denature viral proteins using released ions. However, there have been no reports detailing the effects of changes in humidity or contamination with body fluids on the inactivation of viruses. This study investigated the effects of humidity changes and saliva contamination on the efficacy of SARS-CoV-2 inactivation with ions using multiple viral strains. Virus solutions with different infectious titers were dropped onto a circular nitrocellulose membrane and irradiated with ions from 10 cm above the membrane. After the irradiation of ions for 60, 90, and 120 min, changes in viral infectious titers were measured. The effect of ions on virus inactivation under different humidity conditions was also examined using virus solutions containing 90% mixtures of saliva collected from 10 people. A decrease in viral infectivity was observed over time for all strains, but ion irradiation further accelerated the decrease in viral infectivity. Ion irradiation can inactivate all viral strains, but at 80% humidity, the effect did not appear until 90 min after irradiation. The presence of saliva protected the virus from drying and maintained infectiousness for a longer period compared with no saliva. In particular, the Omicron strain retained its infectivity titer longer than the other strains. Ion irradiation demonstrated a consistent reduction in the number of infectious viruses when compared to the control across varying levels of humidity and irradiation periods. This underscores the notable effectiveness of irradiation, even when the reduction effect is as modest as 50%, thereby emphasizing its crucial role in mitigating the rapid dissemination of SARS-CoV-2.

Cryo-electron microscopy-based drug design

Structure-based drug design (SBDD) has gained popularity owing to its ability to develop more potent drugs compared to conventional drug-discovery methods. The success of SBDD relies heavily on obtaining the three-dimensional structures of drug targets. X-ray crystallography is the primary method used for solving structures and aiding the SBDD workflow; however, it is not suitable for all targets. With the resolution revolution, enabling routine high-resolution reconstruction of structures, cryogenic electron microscopy (cryo-EM) has emerged as a promising alternative and has attracted increasing attention in SBDD. Cryo-EM offers various advantages over X-ray crystallography and can potentially replace X-ray crystallography in SBDD. To fully utilize cryo-EM in drug discovery, understanding the strengths and weaknesses of this technique and noting the key advancements in the field are crucial. This review provides an overview of the general workflow of cryo-EM in SBDD and highlights technical innovations that enable its application in drug design. Furthermore, the most recent achievements in the cryo-EM methodology for drug discovery are discussed, demonstrating the potential of this technique for advancing drug development. By understanding the capabilities and advancements of cryo-EM, researchers can leverage the benefits of designing more effective drugs. This review concludes with a discussion of the future perspectives of cryo-EM-based SBDD, emphasizing the role of this technique in driving innovations in drug discovery and development. The integration of cryo-EM into the drug design process holds great promise for accelerating the discovery of new and improved therapeutic agents to combat various diseases.

Computational design and engineering of self-assembling multivalent microproteins with therapeutic potential against SARS-CoV-2

Multivalent drugs targeting homo-oligomeric viral surface proteins, such as the SARS-CoV-2 trimeric spike (S) protein, have the potential to elicit more potent and broad-spectrum therapeutic responses than monovalent drugs by synergistically engaging multiple binding sites on viral targets. However, rational design and engineering of nanoscale multivalent protein drugs are still lacking. Here, we developed a computational approach to engineer self-assembling trivalent microproteins that simultaneously bind to the three receptor binding domains (RBDs) of the S protein. This approach involves four steps: structure-guided linker design, molecular simulation evaluation of self-assembly, experimental validation of self-assembly state, and functional testing. Using this approach, we first designed trivalent constructs of the microprotein miniACE2 (MP) with different trimerization scaffolds and linkers, and found that one of the constructs (MP-5ff) showed high trimerization efficiency, good conformational homogeneity, and strong antiviral neutralizing activity. With its trimerization unit (5ff), we then engineered a trivalent nanobody (Tr67) that exhibited potent and broad neutralizing activity against the dominant Omicron variants, including XBB.1 and XBB.1.5. Cryo-EM complex structure confirmed that Tr67 stably binds to all three RBDs of the Omicron S protein in a synergistic form, locking them in the "3-RBD-up" conformation that could block human receptor (ACE2) binding and potentially facilitate immune clearance. Therefore, our approach provides an effective strategy for engineering potent protein drugs against SARS-CoV-2 and other deadly coronaviruses.

Protein-protein and protein-nucleic acid binding site prediction via interpretable hierarchical geometric deep learning

Identification of protein-protein and protein-nucleic acid binding sites provides insights into biological processes related to protein functions and technical guidance for disease diagnosis and drug design. However, accurate predictions by computational approaches remain highly challenging due to the limited knowledge of residue binding patterns. The binding pattern of a residue should be characterized by the spatial distribution of its neighboring residues combined with their physicochemical information interaction, which yet cannot be achieved by previous methods. Here, we design GraphRBF, a hierarchical geometric deep learning model to learn residue binding patterns from big data. To achieve it, GraphRBF describes physicochemical information interactions by designing an enhanced graph neural network and characterizes residue spatial distributions by introducing a prioritized radial basis function neural network. After training and testing, GraphRBF shows great improvements over existing state-of-the-art methods and strong interpretability of its learned representations. Applying GraphRBF to the SARS-CoV-2 omicron spike protein, it successfully identifies known epitopes of the protein. Moreover, it predicts multiple potential binding regions for new nanobodies or even new drugs with strong evidence. A user-friendly online server for GraphRBF is freely available at http://liulab.top/GraphRBF/server.

Preclinical evaluation of RQ3013, a broad-spectrum mRNA vaccine against SARS-CoV-2 variants

The global emergence of SARS-CoV-2 variants has led to increasing breakthrough infections in vaccinated populations, calling for an urgent need to develop more effective and broad-spectrum vaccines to combat COVID-19. Here we report the preclinical development of RQ3013, an mRNA vaccine candidate intended to bring broad protection against SARS-CoV-2 variants of concern (VOCs). RQ3013, which contains pseudouridine-modified mRNAs formulated in lipid nanoparticles, encodes the spike (S) protein harboring a combination of mutations responsible for immune evasion of VOCs. Here we characterized the expressed S immunogen and evaluated the immunogenicity, efficacy, and safety of RQ3013 in various animal models. RQ3013 elicited robust immune responses in mice, hamsters, and nonhuman primates (NHP). It can induce high titers of antibodies with broad cross-neutralizing ability against the wild-type, B.1.1.7, B.1.351, B.1.617.2, and the newly emerging Omicron variants. In mice and NHP, two doses of RQ3013 protected the upper and lower respiratory tract against infection by SARS-CoV-2 and its variants. Furthermore, our safety assessment of RQ3013 in NHP showed no observable adverse effects. These results provide strong support for the evaluation of RQ3013 in clinical trials and suggest that it may be a promising candidate for broad protection against COVID-19 and its variants.

Evolution of the SARS-CoV-2 Omicron spike

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron variant of concern, first identified in November 2021, rapidly spread worldwide and diversified into several subvariants. The Omicron spike (S) protein accumulated an unprecedented number of sequence changes relative to previous variants. In this review, we discuss how Omicron S protein structural features modulate host cell receptor binding, virus entry, and immune evasion and highlight how these structural features differentiate Omicron from previous variants. We also examine how key structural properties track across the still-evolving Omicron subvariants and the importance of continuing surveillance of the S protein sequence evolution over time.

Cooperative and structural relationships of the trimeric Spike with infectivity and antibody escape of the strains Delta (B.1.617.2) and Omicron (BA.2, BA.5, and BQ.1)

Herein, we conducted simulations of trimeric Spike from several SARS-CoV-2 variants of concern (Delta and Omicron sub-variants BA.2, BA.5, and BQ.1) and investigated the mechanisms by which specific mutations confer resistance to neutralizing antibodies. We observed that the mutations primarily affect the cooperation between protein domains within and between protomers. The substitutions K417N and L452R expand hydrogen bonding interactions, reducing their interaction with neutralizing antibodies. By interacting with nearby residues, the K444T and N460K mutations in the SpikeBQ.1 variant potentially reduces solvent exposure, thereby promoting resistance to antibodies. We also examined the impact of D614G, P681R, and P681H substitutions on Spike protein structure that may be related to infectivity. The D614G substitution influences communication between a glycine residue and neighboring domains, affecting the transition between up- and -down RBD states. The P681R mutation, found in the Delta variant, enhances correlations between protein subunits, while the P681H mutation in Omicron sub-variants weakens long-range interactions that may be associated with reduced fusogenicity. Using a multiple linear regression model, we established a connection between inter-protomer communication and loss of sensitivity to neutralizing antibodies. Our findings underscore the importance of structural communication between protein domains and provide insights into potential mechanisms of immune evasion by SARS-CoV-2. Overall, this study deepens our understanding of how specific mutations impact SARS-CoV-2 infectivity and shed light on how the virus evades the immune system.

Discovery of Nanosota-2, -3, and -4 as super potent and broad-spectrum therapeutic nanobody candidates against COVID-19

IMPORTANCE

The COVID-19 pandemic exposed limitations of conventional antibodies as therapeutics, including high cost, limited potency, ineffectiveness against new viral variants, and primary reliance on injection-only delivery. Nanobodies are single-domain antibodies with therapeutic potentials. We discovered three anti-SARS-CoV-2 nanobodies, named Nanosota-2, -3, and -4, from an immunized alpaca. Nanosota-2 is super potent against prototypic SARS-CoV-2, Nanosota-3 is highly potent against the omicron variant, and Nanosota-4 is effective against both SARS-CoV-1 and SARS-CoV-2. In addition to their super potency and combined broad antiviral spectrum, these nanobodies are cost-effective, can be easily adapted to new viral variants through phage display, and can potentially be administered as inhalers. The Nanosota series are powerful therapeutic candidates to combat circulating SARS-CoV-2 and prepare for possible future coronavirus pandemics.

Lys417 acts as a molecular switch that regulates the conformation of SARS-CoV-2 spike protein

SARS-CoV-2 spike protein plays a key role in mediating viral entry and inducing host immune responses. It can adopt either an open or closed conformation based on the position of its receptor-binding domain (RBD). It is yet unclear what causes these conformational changes or how they influence the spike's functions. Here, we show that Lys417 in the RBD plays dual roles in the spike's structure: it stabilizes the closed conformation of the trimeric spike by mediating inter-spike-subunit interactions; it also directly interacts with ACE2 receptor. Hence, a K417V mutation has opposing effects on the spike's function: it opens up the spike for better ACE2 binding while weakening the RBD's direct binding to ACE2. The net outcomes of this mutation are to allow the spike to bind ACE2 with higher probability and mediate viral entry more efficiently, but become more exposed to neutralizing antibodies. Given that residue 417 has been a viral mutational hotspot, SARS-CoV-2 may have been evolving to strike a balance between infection potency and immune evasion, contributing to its pandemic spread.

Adsorption of pulmonary and exogeneous surfactants on SARS-CoV-2 spike protein

COVID-19 is transmitted by airborne particles containing virions of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Coronavirus virions represent nanoparticles enveloped by a lipid bilayer decorated by a "crown" of Spike protein protrusions. Virus transmission into the cells is induced by binding of Spike proteins with ACE2 receptors of alveolar epithelial cells. Active clinical search is ongoing for exogenous surfactants and biologically active chemicals capable of hindering virion-receptor binding. Here, we explore by using coarse-grained molecular dynamics simulations the physico-chemical mechanisms of adsorption of selected pulmonary surfactants, zwitterionic dipalmitoyl phosphatidyl choline and cholesterol, and exogeneous anionic surfactant, sodium dodecyl sulfate, on the S1-domain of the Spike protein. We show that surfactants form micellar aggregates that selectively adhere to the specific regions of the S1-domain that are responsible for binding with ACE2 receptors. We find distinctly higher cholesterol adsorption and stronger cholesterol-S1 interactions in comparison with other surfactants, that is consistent with the experimental observations of the effects of cholesterol on COVID-19 infection. Distribution of adsorbed surfactant along the protein residue chain is highly specific and inhomogeneous with preferential adsorption around specific amino acid sequences. We observe preferential adsorption of surfactants on cationic arginine and lysine residues in the receptor-binding domain (RBD) that play an important role in ACE2 binding and are present in higher amounts in Delta and Omicron variants, which may lead to blocking direct Spike-ACE2 interactions. Our findings of strong selective adhesion of surfactant aggregates to Spike proteins have important implications for informing clinical search for therapeutic surfactants for curing and preventing COVID-19 caused by SARS-CoV-2 and its variants.

Structural understanding of SARS-CoV-2 virus entry to host cells

Coronavirus disease 2019 (COVID-19), caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is a major global health concern associated with millions of fatalities worldwide. Mutant variants of the virus have further exacerbated COVID-19 mortality and infection rates, emphasizing the urgent need for effective preventive strategies. Understanding the viral infection mechanism is crucial for developing therapeutics and vaccines. The entry of SARS-CoV-2 into host cells is a key step in the infection pathway and has been targeted for drug development. Despite numerous reviews of COVID-19 and the virus, there is a lack of comprehensive reviews focusing on the structural aspects of viral entry. In this review, we analyze structural changes in Spike proteins during the entry process, dividing the entry process into prebinding, receptor binding, proteolytic cleavage, and membrane fusion steps. By understanding the atomic-scale details of viral entry, we can better target the entry step for intervention strategies. We also examine the impacts of mutations in Spike proteins, including the Omicron variant, on viral entry. Structural information provides insights into the effects of mutations and can guide the development of therapeutics and vaccines. Finally, we discuss available structure-based approaches for the development of therapeutics and vaccines. Overall, this review provides a detailed analysis of the structural aspects of SARS-CoV-2 viral entry, highlighting its significance in the development of therapeutics and vaccines against COVID-19. Therefore, our review emphasizes the importance of structural information in combating SARS-CoV-2 infection.

Spike Protein Mutation-Induced Changes in the Kinetic and Thermodynamic Behavior of Its Receptor Binding Domains Explain Their Higher Propensity to Attain Open States in SARS-CoV-2 Variants of Concern

Spike (S) protein opening in SARS-CoV-2 controls the accessibility of its receptor binding domains (RBDs) to host receptors and immune recognition. Along the evolution of SARS-CoV-2 to its variants of concern (VOC)-alpha, beta, gamma, delta, and omicron-their S proteins showed a higher propensity to attain open states. Deciphering how mutations in S protein can shape its conformational dynamics will contribute to the understanding of viral host tropism. Here using microsecond-scale multiple molecular dynamics simulations (MDS), we provide insights into the kinetic and thermodynamic contributions of these mutations to RBD opening pathways in S proteins of SARS-CoV-2 VOCs. Mutational effects were analyzed using atomistic (i) equilibrium MDS of closed and open states of S proteins and (ii) nonequilibrium MDS for closed-to-open transitions. In MDS of closed or open states, RBDs in S proteins of VOCs showed lower thermodynamic stability with higher kinetic fluctuations, compared to S proteins of ancestral SARS-CoV-2. For closed-to-open transitions in S proteins of VOCs, we observed apparently faster RBD opening with a 1.5-2-fold decrease in the thermodynamic free-energy barrier (ΔGclosed→open). Saturation mutagenesis studies highlighted S protein mutations that may control its conformational dynamics and presentation to host receptors.

The Transmembrane Protease Serine 2 (TMPRSS2) Non-Protease Domains Regulating Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) Spike-Mediated Virus Entry

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) enters cells by binding to the angiotensin-converting enzyme 2 (hACE2) receptor. This process is aided by the transmembrane protease serine 2 (TMPRSS2), which enhances entry efficiency and infectiousness by cleaving the SARS-CoV-2 surface glycoprotein (Spike). The cleavage primes the Spike protein, promoting membrane fusion instead of receptor-mediated endocytosis. Despite the pivotal role played by TMPRSS2, our understanding of its non-protease distinct domains remains limited. In this report, we present evidence indicating the potential phosphorylation of a minimum of six tyrosine residues within the cytosolic tail (CT) of TMPRSS2. Via the use of TMPRSS2 CT phospho-mimetic mutants, we observed a reduction in TMPRSS2 protease activity, accompanied by a decrease in SARS-CoV-2 pseudovirus transduction, which was found to occur mainly via the endosomal pathway. We expanded our investigation beyond TMPRSS2 CT and discovered the involvement of other non-protease domains in regulating infection. Our co-immunoprecipitation experiments demonstrated a strong interaction between TMPRSS2 and Spike. We revealed a 21 amino acid long TMPRSS2-Spike-binding region (TSBR) within the TMPRSS2 scavenger receptor cysteine-rich (SRCR) domain that contributes to this interaction. Our study sheds light on novel functionalities associated with TMPRSS2's cytosolic tail and SRCR region. Both of these regions have the capability to regulate SARS-CoV-2 entry pathways. These findings contribute to a deeper understanding of the complex interplay between viral entry and host factors, opening new avenues for potential therapeutic interventions.

Variant-dependent oxidative and cytokine responses of human neutrophils to SARS-CoV-2 spike protein and anti-spike IgG1 antibodies

Introduction

Severe forms of COVID-19, the disease caused by SARS-CoV-2, are characterized by acute respiratory distress syndrome, robust lung inflammation and death in some patients. Strong evidence has been accumulating that polymorphonuclear neutrophilic granulocytes (PMN) play an important role in the pathophysiology of severe COVID-19. SARS-CoV-2 directly induces in vitro PMN activation, mainly the release of neutrophil extracellular traps (NETs). However, the viral components inducing this PMN response remain unclear.

Methods

In this work human PMN responses were assessed in vitro in response to the spike (S) protein of two different SARS-CoV-2 variants, anti-S IgG1 antibodies or immune complexes formed by them. Production of reactive oxygen species (ROS) was measured by Diogenes-based chemiluminescence. Release of myeloperoxidase (MPO) was assessed by ELISA while secretion of a list of cytokines and growth factors was determined by high-performance multiplex cytokine assay.

Results and discussion

We show that the SARS-CoV-2 Omicron variant S protein and anti-spike IgG1, either alone or together, stimulate ROS production in human PMNs. We also observed that the SARS-CoV-2 Wuhan S protein and anti-S IgG1 antibody together trigger MPO release from PMNs. Based on the relevance of SARS-CoV-2 and influenza co-infections, we have also investigated the impact of influenza virus infection on the previous PMN responses to S proteins or anti-S antibodies. We did not detect any significant effect of influenza co-infection on ROS generation in PMNs. Our data also show that PMN stimulation by S proteins induced the release of different chemokines, growth factors, regulatory and proinflammatory cytokines. Overall, our findings show that the SARS-CoV-2 S protein, an anti-spike IgG1 antibody or their immune complex, promote oxidative responses of PMNs in a variant-dependent manner, contributing to a better understanding of the role of PMN responses during SARS-CoV-2 infection.

Comparative Analysis of Conformational Dynamics and Systematic Characterization of Cryptic Pockets in the SARS-CoV-2 Omicron BA.2, BA.2.75 and XBB.1 Spike Complexes with the ACE2 Host Receptor: Confluence of Binding and Structural Plasticity in Mediating Networks of Conserved Allosteric Sites

In the current study, we explore coarse-grained simulations and atomistic molecular dynamics together with binding energetics scanning and cryptic pocket detection in a comparative examination of conformational landscapes and systematic characterization of allosteric binding sites in the SARS-CoV-2 Omicron BA.2, BA.2.75 and XBB.1 spike full-length trimer complexes with the host receptor ACE2. Microsecond simulations, Markov state models and mutational scanning of binding energies of the SARS-CoV-2 BA.2 and BA.2.75 receptor binding domain complexes revealed the increased thermodynamic stabilization of the BA.2.75 variant and significant dynamic differences between these Omicron variants. Molecular simulations of the SARS-CoV-2 Omicron spike full-length trimer complexes with the ACE2 receptor complemented atomistic studies and enabled an in-depth analysis of mutational and binding effects on conformational dynamic and functional adaptability of the Omicron variants. Despite considerable structural similarities, Omicron variants BA.2, BA.2.75 and XBB.1 can induce unique conformational dynamic signatures and specific distributions of the conformational states. Using conformational ensembles of the SARS-CoV-2 Omicron spike trimer complexes with ACE2, we conducted a comprehensive cryptic pocket screening to examine the role of Omicron mutations and ACE2 binding on the distribution and functional mechanisms of the emerging allosteric binding sites. This analysis captured all experimentally known allosteric sites and discovered networks of inter-connected and functionally relevant allosteric sites that are governed by variant-sensitive conformational adaptability of the SARS-CoV-2 spike structures. The results detailed how ACE2 binding and Omicron mutations in the BA.2, BA.2.75 and XBB.1 spike complexes modulate the distribution of conserved and druggable allosteric pockets harboring functionally important regions. The results are significant for understanding the functional roles of druggable cryptic pockets that can be used for allostery-mediated therapeutic intervention targeting conformational states of the Omicron variants.

SARS-CoV-2 proteins structural studies using synchrotron radiation

In the process of the development of structural biology, both the size and the complexity of the determined macromolecular structures have grown significantly. As a result, the range of application areas for the results of structural studies of biological macromolecules has expanded. Significant progress in the development of structural biology methods has been largely achieved through the use of synchrotron radiation. Modern sources of synchrotron radiation allow to conduct high-performance structural studies with high temporal and spatial resolution. Thus, modern techniques make it possible to obtain not only static structures, but also to study dynamic processes, which play a key role in understanding biological mechanisms. One of the key directions in the development of structural research is the drug design based on the structures of biomolecules. Synchrotron radiation offers insights into the three-dimensional time-resolved structure of individual viral proteins and their complexes at atomic resolution. The rapid and accurate determination of protein structures is crucial for understanding viral pathogenicity and designing targeted therapeutics. Through the application of experimental techniques, including X-ray crystallography and small-angle X-ray scattering (SAXS), it is possible to elucidate the structural details of SARS-CoV-2 virion containing 4 structural, 16 nonstructural proteins (nsp), and several accessory proteins. The most studied potential targets for vaccines and drugs are the structural spike (S) protein, which is responsible for entering the host cell, as well as nonstructural proteins essential for replication and transcription, such as main protease (Mpro), papain-like protease (PLpro), and RNA-dependent RNA polymerase (RdRp). This article provides a brief overview of structural analysis techniques, with focus on synchrotron radiation-based methods applied to the analysis of SARS-CoV-2 proteins.

Exploring Conformational Landscapes and Cryptic Binding Pockets in Distinct Functional States of the SARS-CoV-2 Omicron BA.1 and BA.2 Trimers: Mutation-Induced Modulation of Protein Dynamics and Network-Guided Prediction of Variant-Specific Allosteric Binding Sites

A significant body of experimental structures of SARS-CoV-2 spike trimers for the BA.1 and BA.2 variants revealed a considerable plasticity of the spike protein and the emergence of druggable binding pockets. Understanding the interplay of conformational dynamics changes induced by the Omicron variants and the identification of cryptic dynamic binding pockets in the S protein is of paramount importance as exploring broad-spectrum antiviral agents to combat the emerging variants is imperative. In the current study, we explore conformational landscapes and characterize the universe of binding pockets in multiple open and closed functional spike states of the BA.1 and BA.2 Omicron variants. By using a combination of atomistic simulations, a dynamics network analysis, and an allostery-guided network screening of binding pockets in the conformational ensembles of the BA.1 and BA.2 spike conformations, we identified all experimentally known allosteric sites and discovered significant variant-specific differences in the distribution of binding sites in the BA.1 and BA.2 trimers. This study provided a structural characterization of the predicted cryptic pockets and captured the experimentally known allosteric sites, revealing the critical role of conformational plasticity in modulating the distribution and cross-talk between functional binding sites. We found that mutational and dynamic changes in the BA.1 variant can induce the remodeling and stabilization of a known druggable pocket in the N-terminal domain, while this pocket is drastically altered and may no longer be available for ligand binding in the BA.2 variant. Our results predicted the experimentally known allosteric site in the receptor-binding domain that remains stable and ranks as the most favorable site in the conformational ensembles of the BA.2 variant but could become fragmented and less probable in BA.1 conformations. We also uncovered several cryptic pockets formed at the inter-domain and inter-protomer interface, including functional regions of the S2 subunit and stem helix region, which are consistent with the known role of pocket residues in modulating conformational transitions and antibody recognition. The results of this study are particularly significant for understanding the dynamic and network features of the universe of available binding pockets in spike proteins, as well as the effects of the Omicron-variant-specific modulation of preferential druggable pockets. The exploration of predicted druggable sites can present a new and previously underappreciated opportunity for therapeutic interventions for Omicron variants through the conformation-selective and variant-specific targeting of functional sites involved in allosteric changes.

Examining Functional Linkages Between Conformational Dynamics, Protein Stability and Evolution of Cryptic Binding Pockets in the SARS-CoV-2 Omicron Spike Complexes with the ACE2 Host Receptor: Recombinant Omicron Variants Mediate Variability of Conserved Allosteric Sites and Binding Epitopes

In the current study, we explore coarse-grained simulations and atomistic molecular dynamics together with binding energetics scanning and cryptic pocket detection in a comparative examination of conformational landscapes and systematic characterization of allosteric binding sites in the SARS-CoV-2 Omicron BA.2, BA.2.75 and XBB.1 spike full-length trimer complexes with the host receptor ACE2. Microsecond simulations, Markov state models and mutational scanning of binding energies of the SARS-CoV-2 BA.2 and BA.2.75 receptor binding domain complexes revealed the increased thermodynamic stabilization of the BA.2.75 variant and significant dynamic differences between these Omicron variants. Molecular simulations of the SARS-CoV-2 Omicron spike full length trimer complexes with the ACE2 receptor complemented atomistic studies and enabled an in-depth analysis of mutational and binding effects on conformational dynamic and functional adaptability of the Omicron variants. Despite considerable structural similarities, Omicron variants BA.2, BA.2.75 and XBB.1 can induce unique conformational dynamic signatures and specific distributions of the conformational states. Using conformational ensembles of the SARS-CoV-2 Omicron spike trimer complexes with ACE2, we conducted a comprehensive cryptic pocket screening to examine the role of Omicron mutations and ACE2 binding on the distribution and functional mechanisms of the emerging allosteric binding sites. This analysis captured all experimentally known allosteric sites and discovered networks of inter-connected and functionally relevant allosteric sites that are governed by variant-sensitive conformational adaptability of the SARS-CoV-2 spike structures. The results detailed how ACE2 binding and Omicron mutations in the BA.2, BA.2.75 and XBB.1 spike complexes modulate the distribution of conserved and druggable allosteric pockets harboring functionally important regions. The results of are significant for understanding functional roles of druggable cryptic pockets that can be used for allostery-mediated therapeutic intervention targeting conformational states of the Omicron variants.

Markov State Models and Perturbation-Based Approaches Reveal Distinct Dynamic Signatures and Hidden Allosteric Pockets in the Emerging SARS-Cov-2 Spike Omicron Variant Complexes with the Host Receptor: The Interplay of Dynamics and Convergent Evolution Modulates Allostery and Functional Mechanisms

The new generation of SARS-CoV-2 Omicron variants displayed a significant growth advantage and increased viral fitness by acquiring convergent mutations, suggesting that the immune pressure can promote convergent evolution leading to the sudden acceleration of SARS-CoV-2 evolution. In the current study, we combined structural modeling, microsecond molecular dynamics simulations, and Markov state models to characterize conformational landscapes and identify specific dynamic signatures of the SARS-CoV-2 spike complexes with the host receptor ACE2 for the recently emerged highly transmissible XBB.1, XBB.1.5, BQ.1, and BQ.1.1 Omicron variants. Microsecond simulations and Markovian modeling provided a detailed characterization of the functional conformational states and revealed the increased thermodynamic stabilization of the XBB.1.5 subvariant, which can be contrasted to more dynamic BQ.1 and BQ.1.1 subvariants. Despite considerable structural similarities, Omicron mutations can induce unique dynamic signatures and specific distributions of the conformational states. The results suggested that variant-specific changes of the conformational mobility in the functional interfacial loops of the receptor-binding domain in the SARS-CoV-2 spike protein can be fine-tuned through crosstalk between convergent mutations which could provide an evolutionary path for modulation of immune escape. By combining atomistic simulations and Markovian modeling analysis with perturbation-based approaches, we determined important complementary roles of convergent mutation sites as effectors and receivers of allosteric signaling involved in modulation of conformational plasticity and regulation of allosteric communications. This study also revealed hidden allosteric pockets and suggested that convergent mutation sites could control evolution and distribution of allosteric pockets through modulation of conformational plasticity in the flexible adaptable regions.

Probing conformational landscapes of binding and allostery in the SARS-CoV-2 omicron variant complexes using microsecond atomistic simulations and perturbation-based profiling approaches: hidden role of omicron mutations as modulators of allosteric signaling and epistatic relationships

In this study, we systematically examine the conformational dynamics, binding and allosteric communications in the Omicron BA.1, BA.2, BA.3 and BA.4/BA.5 spike protein complexes with the ACE2 host receptor using molecular dynamics simulations and perturbation-based network profiling approaches. Microsecond atomistic simulations provided a detailed characterization of the conformational landscapes and revealed the increased thermodynamic stabilization of the BA.2 variant which can be contrasted with the BA.4/BA.5 variants inducing a significant mobility of the complexes. Using the dynamics-based mutational scanning of spike residues, we identified structural stability and binding affinity hotspots in the Omicron complexes. Perturbation response scanning and network-based mutational profiling approaches probed the effect of the Omicron mutations on allosteric interactions and communications in the complexes. The results of this analysis revealed specific roles of Omicron mutations as conformationally plastic and evolutionary adaptable modulators of binding and allostery which are coupled to the major regulatory positions through interaction networks. Through perturbation network scanning of allosteric residue potentials in the Omicron variant complexes performed in the background of the original strain, we characterized regions of epistatic couplings that are centered around the binding affinity hotspots N501Y and Q498R. Our results dissected the vital role of these epistatic centers in regulating protein stability, efficient ACE2 binding and allostery which allows for accumulation of multiple Omicron immune escape mutations at other sites. Through integrative computational approaches, this study provides a systematic analysis of the effects of Omicron mutations on thermodynamics, binding and allosteric signaling in the complexes with ACE2 receptor.

Mechanism and evolution of human ACE2 binding by SARS-CoV-2 spike

Spike glycoprotein of SARS-CoV-2 mediates viral entry into host cells by facilitating virus attachment and membrane fusion. ACE2 is the main receptor of SARS-CoV-2 and its interaction with spike has shaped the virus' emergence from an animal reservoir and subsequent evolution in the human host. Many structural studies on the spike:ACE2 interaction have provided insights into mechanisms driving viral evolution during the on-going pandemic. This review describes the molecular basis of spike binding to ACE2, outlines mechanisms that have optimised this interaction during viral evolution, and suggests directions for future research.

The effects of amino acid substitution of spike protein and genomic recombination on the evolution of SARS-CoV-2

Over three years' pandemic of 2019 novel coronavirus disease (COVID-19), multiple variants and novel subvariants have emerged successively, outcompeted earlier variants and become predominant. The sequential emergence of variants reflects the evolutionary process of mutation-selection-adaption of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Amino acid substitution/insertion/deletion in the spike protein causes altered viral antigenicity, transmissibility, and pathogenicity of SARS-CoV-2. Early in the pandemic, D614G mutation conferred virus with advantages over previous variants and increased transmissibility, and it also laid a conservative background for subsequent substantial mutations. The role of genomic recombination in the evolution of SARS-CoV-2 raised increasing concern with the occurrence of novel recombinants such as Deltacron, XBB.1.5, XBB.1.9.1, and XBB.1.16 in the late phase of pandemic. Co-circulation of different variants and co-infection in immunocompromised patients accelerate the emergence of recombinants. Surveillance for SARS-CoV-2 genomic variations, particularly spike protein mutation and recombination, is essential to identify ongoing changes in the viral genome and antigenic epitopes and thus leads to the development of new vaccine strategies and interventions.

Self-Association of ACE-2 with Different RBD Amounts: A Dynamic Simulation Perspective on SARS-CoV-2 Infection

Transmissibility of SARS-CoV-2 initially relies on its trimeric Spike-RBDs to tether the ACE-2 on host cells, and enhanced self-association of ACE-2 engaged with Spike facilitates the viral infection. Two primary packing modes of Spike-ACE2 heteroproteins exist potentially due to discrepant amounts of RBDs loading on ACE-2, but the resultant self-association difference is inherently unclear. We used extensive coarse-grained dynamic simulations to characterize the self-association efficiency, the conformation relevance, and the molecular mechanism of ACE-2 with different RBD amounts. It was revealed that the ACE-2 hanging two/full RBDs (Mode-A) rapidly dimerized into the heteroprotein complex in a compact "linear" conformation, while the bare ACE-2 showed weakened self-association and a protein complex. The RBD-tethered ectodomains of ACE-2 presented a more upright conformation relative to the membrane, and the intermolecular ectodomains were predominantly packed by the neck domains, which was obligated to the rapid protein self-association in a compact pattern. Noted is the fact that the ACE-2 tethered by a single RBD (Mode-B) retained considerable self-association efficiency and clustering capability, which unravels the interrelation of ACE-2 colocalization and protein cross-linkage. The molecular perspectives in this study expound the self-association potency of ACE-2 with different RBD amounts and the viral activity implications, which can greatly enhance our comprehension of SARS-CoV-2 infection details.

High-Throughput Screening of Promising Redox-Active Molecules with MolGAT

Redox flow batteries (RFBs) have emerged as a promising option for large-scale energy storage, owing to their high energy density, low cost, and environmental benefits. However, the identification of organic compounds with high redox activity, aqueous solubility, stability, and fast redox kinetics is a crucial and challenging step in developing an RFB technology. Density functional theory-based computational materials prediction and screening is a time-consuming and computationally expensive technique, yet it has a high success rate. To speed up the discovery of new materials with desired properties, machine-learning-based models can be trained on large data sets. Graph neural networks (GNNs) are particularly well-suited for non-Euclidean data and can model complex relationships, making them ideal for accelerating the discovery of novel materials. In this study, a GNN-based model called MolGAT was developed to predict the redox potential of organic molecules using molecular structures, atomic properties, and bond attributes. The model was trained on a data set of over 15,000 compounds with redox potentials ranging from -4.11 to 2.56. MolGAT outperformed other GNN variants, such as the Graph Attention Network, Graph Convolution Network, and AttentiveFP models. The trained model was used to screen a vast chemical data set comprising 581,014 molecules, namely OMDB, QM9, ZINC, CHEMBL, and DELANEY, and identified 23,467 potential redox-active compounds for use in redox flow batteries. Of those, 20,716 molecules were identified as potential catholytes with predicted redox potentials up to 2.87 V, while 2,751 molecules were deemed potential anolytes with predicted redox potentials as low as -2.88 V. This work demonstrates the capabilities of graph neural networks in condensed matter physics and materials science to screen promising redox-active species for further electronic structure calculations and experimental testing.

A large expert-curated cryo-EM image dataset for machine learning protein particle picking

Cryo-electron microscopy (cryo-EM) is a powerful technique for determining the structures of biological macromolecular complexes. Picking single-protein particles from cryo-EM micrographs is a crucial step in reconstructing protein structures. However, the widely used template-based particle picking process is labor-intensive and time-consuming. Though machine learning and artificial intelligence (AI) based particle picking can potentially automate the process, its development is hindered by lack of large, high-quality labelled training data. To address this bottleneck, we present CryoPPP, a large, diverse, expert-curated cryo-EM image dataset for protein particle picking and analysis. It consists of labelled cryo-EM micrographs (images) of 34 representative protein datasets selected from the Electron Microscopy Public Image Archive (EMPIAR). The dataset is 2.6 terabytes and includes 9,893 high-resolution micrographs with labelled protein particle coordinates. The labelling process was rigorously validated through 2D particle class validation and 3D density map validation with the gold standard. The dataset is expected to greatly facilitate the development of both AI and classical methods for automated cryo-EM protein particle picking.

Rapid detection of intact SARS-CoV-2 using designer DNA Nets and a pocket-size smartphone-linked fluorimeter

Rapid, sensitive, and inexpensive point-of-care diagnosis is vital to controlling highly infectious diseases, including COVID-19. Here, we report the design and characterization of a compact fluorimeter called a "Virus Pod" (V-Pod) that enables sensitive self-testing of SARS-CoV-2 viral load in saliva. The rechargeable battery-operated device reads the fluorescence generated by Designer DNA Nanostructures (DDN) when they specifically interact with intact SARS-CoV-2 virions. DDNs are net-shaped self-assembling nucleic acid constructs that provide an array of highly specific aptamer-fluorescent quencher duplexes located at precise positions that match the pattern of spike proteins. The room-temperature assay is performed by mixing the test sample with DNA Net sensor in a conventional PCR tube and placing the tube into the V-Pod. Fluorescent signals are generated when multivalent aptamer-spike binding releases fluorescent quenchers, resulting in rapid (5-min) generation of dose-dependent output. The V-Pod instrument performs laser excitation, fluorescence intensity quantitation, and secure transmission of data to an App via Bluetooth™. We show that the V-Pod and DNA Net assay achieves clinically relevant detection limits of 3.92 × 103 viral-genome-copies/mL for pseudo-typed wild-type SARS-CoV-2 and 1.84 × 104, 9.69 × 104, 6.99 × 104 viral-genome-copies/mL for pathogenic Delta, Omicron, and D614G variants, representing sensitivity similar to laboratory-based PCR. The pocket-sized instrument (∼$294), inexpensive reagent-cost/test ($1.26), single-step, rapid sample-to-answer, and quantitative output represent a capability that is compatible with the needs of frequent self-testing in a consumer-friendly format that can link with medical service systems such as healthcare providers, contact tracing, and infectious disease reporting.

Correlating the differences in the receptor binding domain of SARS-CoV-2 spike variants on their interactions with human ACE2 receptor

Spike glycoprotein of SARS-CoV-2 variants plays a critical role in infection and transmission through its interaction with human angiotensin converting enzyme 2 (hACE2) receptors. Prior findings using molecular docking and biomolecular studies reported varied findings on the difference in the interactions among the spike variants with the hACE2 receptors. Hence, it is a prerequisite to understand these interactions in a more precise manner. To this end, firstly, we performed ELISA with trimeric spike glycoproteins of SARS-CoV-2 variants including Wuhan Hu-1(Wild), Delta, C.1.2 and Omicron. Further, to study the interactions in a more specific manner by mimicking the natural infection, we developed hACE2 receptors expressing HEK-293T cell line, evaluated their binding efficiencies and competitive binding of spike variants with D614G spike pseudotyped virus. In line with the existing findings, we observed that Omicron had higher binding efficiency compared to Delta in both ELISA and Cellular models. Intriguingly, we found that cellular models could differentiate the subtle differences between the closely related C.1.2 and Delta in their binding to hACE2 receptors. Our study using the cellular model provides a precise method to evaluate the binding interactions between spike sub-lineages to hACE2 receptors.

Markov State Models and Perturbation-Based Approaches Reveal Distinct Dynamic Signatures and Hidden Allosteric Pockets in the Emerging SARS-Cov-2 Spike Omicron Variants Complexes with the Host Receptor: The Interplay of Dynamics and Convergent Evolution Modulates Allostery and Functional Mechanisms

The new generation of SARS-CoV-2 Omicron variants displayed a significant growth advantage and the increased viral fitness by acquiring convergent mutations, suggesting that the immune pressure can promote convergent evolution leading to the sudden acceleration of SARS-CoV-2 evolution. In the current study, we combined structural modeling, extensive microsecond MD simulations and Markov state models to characterize conformational landscapes and identify specific dynamic signatures of the SARS-CoV-2 spike complexes with the host receptor ACE2 for the recently emerged highly transmissible XBB.1, XBB.1.5, BQ.1, and BQ.1.1 Omicron variants. Microsecond simulations and Markovian modeling provided a detailed characterization of the conformational landscapes and revealed the increased thermodynamic stabilization of the XBB.1.5 subvariant which is contrasted to more dynamic BQ.1 and BQ.1.1 subvariants. Despite considerable structural similarities, Omicron mutations can induce unique dynamic signatures and specific distributions of conformational states. The results suggested that variant-specific changes of conformational mobility in the functional interfacial loops of the spike receptor binding domain can be fine-tuned through cross-talk between convergent mutations thereby providing an evolutionary path for modulation of immune escape. By combining atomistic simulations and Markovian modeling analysis with perturbation-based approaches, we determined important complementary roles of convergent mutation sites as effectors and receivers of allosteric signaling involved in modulating conformational plasticity at the binding interface and regulating allosteric responses. This study also characterized the dynamics-induced evolution of allosteric pockets in the Omicron complexes that revealed hidden allosteric pockets and suggested that convergent mutation sites could control evolution and distribution of allosteric pockets through modulation of conformational plasticity in the flexible adaptable regions. Through integrative computational approaches, this investigation provides a systematic analysis and comparison of the effects of Omicron subvariants on conformational dynamics and allosteric signaling in the complexes with the ACE2 receptor.

For Table of Contents Use Only

Nanoparticle Spikes Enhance Cellular Uptake via Regulating Myosin IIA Recruitment

Spike-like nanostructures are omnipresent in natural and artificial systems. Although biorecognition of nanostructures to cellular receptors has been indicated as the primary factor for virus infection pathways, how the spiky morphology of DNA-modified nanoparticles affects their cellular uptake and intracellular fate remains to be explored. Here, we design dually emissive gold nanoparticles with varied spikiness (from 0 to 2) to probe the interactions of spiky nanoparticles with cells. We discovered that nanospikes at the nanoparticle regulated myosin IIA recruitment at the cell membrane during cellular uptake, thereby enhancing cellular uptake efficiency, as revealed by dual-modality (plasmonic and fluorescence) imaging. Furthermore, the spiky nanoparticles also exhibited facilitated endocytosis dynamics, as revealed by real-time dark-field microscopy (DFM) imaging and colorimetry-based classification algorithms. These findings highlight the crucial role of the spiky morphology in regulating the intracellular fate of nanoparticles, which may shed light on engineering theranostic nanocarriers.

Balancing Functional Tradeoffs between Protein Stability and ACE2 Binding in the SARS-CoV-2 Omicron BA.2, BA.2.75 and XBB Lineages: Dynamics-Based Network Models Reveal Epistatic Effects Modulating Compensatory Dynamic and Energetic Changes

Evolutionary and functional studies suggested that the emergence of the Omicron variants can be determined by multiple fitness trade-offs including the immune escape, binding affinity for ACE2, conformational plasticity, protein stability and allosteric modulation. In this study, we systematically characterize conformational dynamics, structural stability and binding affinities of the SARS-CoV-2 Spike Omicron complexes with the host receptor ACE2 for BA.2, BA.2.75, XBB.1 and XBB.1.5 variants. We combined multiscale molecular simulations and dynamic analysis of allosteric interactions together with the ensemble-based mutational scanning of the protein residues and network modeling of epistatic interactions. This multifaceted computational study characterized molecular mechanisms and identified energetic hotspots that can mediate the predicted increased stability and the enhanced binding affinity of the BA.2.75 and XBB.1.5 complexes. The results suggested a mechanism driven by the stability hotspots and a spatially localized group of the Omicron binding affinity centers, while allowing for functionally beneficial neutral Omicron mutations in other binding interface positions. A network-based community model for the analysis of epistatic contributions in the Omicron complexes is proposed revealing the key role of the binding hotspots R498 and Y501 in mediating community-based epistatic couplings with other Omicron sites and allowing for compensatory dynamics and binding energetic changes. The results also showed that mutations in the convergent evolutionary hotspot F486 can modulate not only local interactions but also rewire the global network of local communities in this region allowing the F486P mutation to restore both the stability and binding affinity of the XBB.1.5 variant which may explain the growth advantages over the XBB.1 variant. The results of this study are consistent with a broad range of functional studies rationalizing functional roles of the Omicron mutation sites that form a coordinated network of hotspots enabling a balance of multiple fitness tradeoffs and shaping up a complex functional landscape of virus transmissibility.

Probing Conformational Landscapes of Binding and Allostery in the SARS-CoV-2 Omicron Variant Complexes Using Microsecond Atomistic Simulations and Perturbation-Based Profiling Approaches: Hidden Role of Omicron Mutations as Modulators of Allosteric Signaling and Epistatic Relationships

In this study, we systematically examine the conformational dynamics, binding and allosteric communications in the Omicron BA.1, BA.2, BA.3 and BA.4/BA.5 complexes with the ACE2 host receptor using molecular dynamics simulations and perturbation-based network profiling approaches. Microsecond atomistic simulations provided a detailed characterization of the conformational landscapes and revealed the increased thermodynamic stabilization of the BA.2 variant which is contrasted with the BA.4/BA.5 variants inducing a significant mobility of the complexes. Using ensemble-based mutational scanning of binding interactions, we identified binding affinity and structural stability hotspots in the Omicron complexes. Perturbation response scanning and network-based mutational profiling approaches probed the effect of the Omicron variants on allosteric communications. The results of this analysis revealed specific roles of Omicron mutations as "plastic and evolutionary adaptable" modulators of binding and allostery which are coupled to the major regulatory positions through interaction networks. Through perturbation network scanning of allosteric residue potentials in the Omicron variant complexes, which is performed in the background of the original strain, we identified that the key Omicron binding affinity hotspots N501Y and Q498R could mediate allosteric interactions and epistatic couplings. Our results suggested that the synergistic role of these hotspots in controlling stability, binding and allostery can enable for compensatory balance of fitness tradeoffs with conformationally and evolutionary adaptable immune-escape Omicron mutations. Through integrative computational approaches, this study provides a systematic analysis of the effects of Omicron mutations on thermodynamics, binding and allosteric signaling in the complexes with ACE2 receptor. The findings support a mechanism in which Omicron mutations can evolve to balance thermodynamic stability and conformational adaptability in order to ensure proper tradeoff between stability, binding and immune escape.

S:D614G and S:H655Y are gateway mutations that act epistatically to promote SARS-CoV-2 variant fitness

SARS-CoV-2 variants bearing complex combinations of mutations that confer increased transmissibility, COVID-19 severity, and immune escape, were first detected after S:D614G had gone to fixation, and likely originated during persistent infection of immunocompromised hosts. To test the hypothesis that S:D614G facilitated emergence of such variants, S:D614G was reverted to the ancestral sequence in the context of sequential Spike sequences from an immunocompromised individual, and within each of the major SARS-CoV-2 variants of concern. In all cases, infectivity of the S:D614G revertants was severely compromised. The infectivity of atypical SARS-CoV-2 lineages that propagated in the absence of S:D614G was found to be dependent upon either S:Q613H or S:H655Y. Notably, Gamma and Omicron variants possess both S:D614G and S:H655Y, each of which contributed to infectivity of these variants. Among sarbecoviruses, S:Q613H, S:D614G, and S:H655Y are only detected in SARS-CoV-2, which is also distinguished by a polybasic S1/S2 cleavage site. Genetic and biochemical experiments here showed that S:Q613H, S:D614G, and S:H655Y each stabilize Spike on virions, and that they are dispensable in the absence of S1/S2 cleavage, consistent with selection of these mutations by the S1/S2 cleavage site. CryoEM revealed that either S:D614G or S:H655Y shift the Spike receptor binding domain (RBD) towards the open conformation required for ACE2-binding and therefore on pathway for infection. Consistent with this, an smFRET reporter for RBD conformation showed that both S:D614G and S:H655Y spontaneously adopt the conformation that ACE2 induces in the parental Spike. Data from these orthogonal experiments demonstrate that S:D614G and S:H655Y are convergent adaptations to the polybasic S1/S2 cleavage site which stabilize S1 on the virion in the open RBD conformation and act epistatically to promote the fitness of variants bearing complex combinations of clinically significant mutations.

Endogenous piRNAs Can Interact with the Omicron Variant of the SARS-CoV-2 Genome

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which caused the COVID-19 pandemic, can still infect populations in many countries around the globe. The Omicron strain is the most mutated variant of SARS-CoV-2. The high transmissibility of the strain and its ability to evade immunity necessitate a priority study of its properties in order to quickly create effective means of preventing its spread. The current research aimed to examine the in silico interaction between PIWI-interacting RNAs (piRNAs) and the SARS-CoV-2 genome (gRNA) to identify endogenous piRNAs and propose synthetic piRNAs with strong antiviral activity for drug development. This study used validated bioinformatic approaches regarding the interaction of more than eight million piRNAs with the SARS-CoV-2 genome. The piRNAs’ binding sites (BSs) in the 5′UTR were located with overlapping nucleotide sequences termed clusters of BSs. Several BSs clusters have been found in the nsp3, nsp7, RNA-dependent RNA polymerase, endoRNAse, S surface glycoprotein, ORF7a, and nucleocapsid. Sixteen synthetic piRNAs that interact with gRNA have been proposed with free binding energy ranging from −170 kJ/mol to −175 kJ/mol, which can be used to create drugs that suppress the reproduction of SARS-CoV-2.

Coarse-Grained Molecular Simulations and Ensemble-Based Mutational Profiling of Protein Stability in the Different Functional Forms of the SARS-CoV-2 Spike Trimers: Balancing Stability and Adaptability in BA.1, BA.2 and BA.2.75 Variants

Evolutionary and functional studies have suggested that the emergence of Omicron variants can be determined by multiple fitness tradeoffs including immune escape, binding affinity, conformational plasticity, protein stability, and allosteric modulation. In this study, we embarked on a systematic comparative analysis of the conformational dynamics, electrostatics, protein stability, and allostery in the different functional states of spike trimers for BA.1, BA.2, and BA.2.75 variants. Using efficient and accurate coarse-grained simulations and atomistic reconstruction of the ensembles, we examined the conformational dynamics of the spike trimers that agree with the recent functional studies, suggesting that BA.2.75 trimers are the most stable among these variants. A systematic mutational scanning of the inter-protomer interfaces in the spike trimers revealed a group of conserved structural stability hotspots that play a key role in the modulation of functional dynamics and are also involved in the inter-protomer couplings through local contacts and interaction networks with the Omicron mutational sites. The results of mutational scanning provided evidence that BA.2.75 trimers are more stable than BA.2 and comparable in stability to the BA.1 variant. Using dynamic network modeling of the S Omicron BA.1, BA.2, and BA.2.75 trimers, we showed that the key network mediators of allosteric interactions are associated with the major stability hotspots that are interconnected along potential communication pathways. The network analysis of the BA.1, BA.2, and BA.2.75 trimers suggested that the increased thermodynamic stability of the BA.2.75 variant may be linked with the organization and modularity of the residue interaction network that allows for allosteric communications between structural stability hotspots and Omicron mutational sites. This study provided a plausible rationale for a mechanism in which Omicron mutations may evolve by targeting vulnerable sites of conformational adaptability to elicit immune escape while maintaining their control on balancing protein stability and functional fitness through robust allosteric communications with the stability hotspots.

Distribution and Functional Analyses of Mutations in Spike Protein and Phylogenic Diversity of SARS-CoV-2 Variants Emerged during the Year 2021 in India

Introduction

Prolonged COVID-19 pandemic accelerates the emergence and transmissibility of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) variants through the accumulation of adaptive mutations. Particularly, adaptive mutations in spike (S) protein of SARS-CoV-2 leads to increased viral infectivity, severe morbidity and mortality, and immune evasion. This study focuses on the phylodynamic distribution of SARS-CoV-2 variants during the year 2021 in India besides analyzing the functional significance of mutations in S-protein of SARS-CoV-2 variants.

Methods

Whole genome of SARS-CoV-2 sequences (n = 87957) from the various parts of India over the period of January to December 2021 was retrieved from Global Initiative on Sharing All Influenza Data. All the S-protein sequences were subjected to clade analysis, variant calling, protein stability, immune escape potential, structural divergence, Furin cleavage efficiency, and phylogenetic analysis using various in silico tools.

Results

Delta variant belonging to 21A, 21I, and 21J clades was found to be predominant throughout the year 2021 though many variants were also present. A total of 4639 amino acid mutations were found in S-protein. D614G was the most predominant mutation in the S-protein followed by P681R, L452R, T19R, T478K, and D950N. The highest number of mutations was found in the N-terminal domain of S-protein. Mutations in the crucial sites of S-protein impacting pathogenicity, immunogenicity, and fusogenicity were identified. Intralineage diversity analysis showed that certain variants of SARS-CoV-2 possess high diversification.

Conclusions

The study has disclosed the distribution of various variants including the Delta, the predominant variant, in India throughout the year 2021. The study has identified mutations in S-protein of each SARS-CoV-2 variant that can significantly impact the virulence, immune evasion, increased transmissibility, high morbidity, and mortality. In addition, it is found that mutations acquired during each viral replication cycle introduce new sub-lineages as studied by intralineage diversity analysis.