Free Energy Landscapes from SARS-CoV-2 Spike Glycoprotein Simulations Suggest that RBD Opening Can Be Modulated via Interactions in an Allosteric Pocket

Assessing pH-Dependent Conformational Changes in the Fusion Peptide Proximal Region of the SARS-CoV-2 Spike Glycoprotein

One of the entry mechanisms of the SARS-CoV-2 coronavirus into host cells involves endosomal acidification. It has been proposed that under acidic conditions, the fusion peptide proximal region (FPPR) of the SARS-CoV-2 spike glycoprotein acts as a pH-dependent switch, modulating immune response accessibility by influencing the positioning of the receptor binding domain (RBD). This would provide indirect coupling of RBD opening to the environmental pH. Here, we explored this possible pH-dependent conformational equilibrium of the FPPR within the SARS-CoV-2 spike glycoprotein. We analyzed hundreds of experimentally determined spike structures from the Protein Data Bank and carried out pH-replica exchange molecular dynamics to explore the extent to which the FPPR conformation depends on pH and the positioning of the RBD. A meta-analysis of experimental structures identified alternate conformations of the FPPR among structures in which this flexible regions was resolved. However, the results did not support a correlation between the FPPR conformation and either RBD position or the reported pH of the cryo-EM experiment. We calculated pKa values for titratable side chains in the FPPR region using PDB structures, but these pKa values showed large differences between alternate PDB structures that otherwise adopt the same FPPR conformation type. This hampers the comparison of pKa values in different FPPR conformations to rationalize a pH-dependent conformational change. We supplemented these PDB-based analyses with all-atom simulations and used constant-pH replica exchange molecular dynamics to estimate pKa values in the context of flexibility and explicit water. The resulting titration curves show good reproducibility between simulations, but they also suggest that the titration curves of the different FPPR conformations are the same within the error bars. In summary, we were unable to find evidence supporting the previously published hypothesis of an FPPR pH-dependent equilibrium: neither from existing experimental data nor from constant-pH MD simulations. The study underscores the complexity of the spike system and opens avenues for further exploration into the interplay between pH and SARS-CoV-2 viral entry mechanisms.

Allosteric pathways of SARS and SARS-CoV-2 spike protein identified by neural relational inference

The receptor binding domain (RBD) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike protein must undergo a crucial conformational transition to invade human cells. It is intriguing that this transition is accompanied by a synchronized movement of the entire spike protein. Therefore, it is possible to design allosteric regulators targeting non-RBD but hindering the conformational transition of RBD. To understand the allosteric mechanism in detail, we establish a computational framework by integrating coarse-grained molecular dynamic simulations and a state-of-the-art neural network model called neural relational inference. Leveraging this framework, we have elucidated the allosteric pathway of the SARS-CoV-2 spike protein at the residue level and identified the molecular mechanisms involved in the transmission of allosteric signals. The movement of D614 is coupled with that of Q321. This interaction subsequently influences the movement of K528/K529, ultimately coupling with the movement of RBD during conformational changes. Mutations that weaken the interactions within this pathway naturally block the allosteric signal transmission, thereby modulating the conformational transitions. This observation also offers a rationale for the distinct allosteric patterns observed in the SARS-CoV spike protein. Our result provides a useful method for analyzing the dynamics of potential viral variants in the future.

The Malaria Box molecules: a source for targeting the RBD and NTD cryptic pocket of the spike glycoprotein in SARS-CoV-2

CONTEXT

SARS-CoV-2, responsible for COVID-19, has led to over 500 million infections and more than 6 million deaths globally. There have been limited effective treatments available. The study aims to find a drug that can prevent the virus from entering host cells by targeting specific sites on the virus's spike protein.

METHOD

We examined 13,397 compounds from the Malaria Box library against two specific sites on the spike protein: the receptor-binding domain (RBD) and a predicted cryptic pocket. Using virtual screening, molecular docking, molecular dynamics, and MMPBSA techniques, they evaluated the stability of two compounds. TCMDC-124223 showed high stability and binding energy in the RBD, while TCMDC-133766 had better binding energy in the cryptic pocket. The study also identified that the interacting residues are conserved, which is crucial for addressing various virus variants. The findings provide insights into the potential of small molecules as drugs against the spike protein.

Ligand- and Structure-Based Virtual Screening Identifies New Inhibitors of the Interaction of the SARS-CoV-2 Spike Protein with the ACE2 Host Receptor

The Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) is a fast-spreading viral pathogen and poses a serious threat to human health. New SARS-CoV-2 variants have been arising worldwide; therefore, is necessary to explore more therapeutic options. The interaction of the viral spike (S) protein with the angiotensin-converting enzyme 2 (ACE2) host receptor is an attractive drug target to prevent the infection via the inhibition of virus cell entry. In this study, Ligand- and Structure-Based Virtual Screening (LBVS and SBVS) was performed to propose potential inhibitors capable of blocking the S receptor-binding domain (RBD) and ACE2 interaction. The best five lead compounds were confirmed as inhibitors through ELISA-based enzyme assays. The docking studies and molecular dynamic (MD) simulations of the selected compounds maintained the molecular interaction and stability (RMSD fluctuations less than 5 Å) with key residues of the S protein. The compounds DRI-1, DRI-2, DRI-3, DRI-4, and DRI-5 efficiently block the interaction between the SARS-CoV-2 spike protein and receptor ACE2 (from 69.90 to 99.65% of inhibition) at 50 µM. The most potent inhibitors were DRI-2 (IC50 = 8.8 µM) and DRI-3 (IC50 = 2.1 µM) and have an acceptable profile of cytotoxicity (CC50 > 90 µM). Therefore, these compounds could be good candidates for further SARS-CoV-2 preclinical experiments.

N-Glycosylation as a Modulator of Protein Conformation and Assembly in Disease

Glycosylation, a prevalent post-translational modification, plays a pivotal role in regulating intricate cellular processes by covalently attaching glycans to macromolecules. Dysregulated glycosylation is linked to a spectrum of diseases, encompassing cancer, neurodegenerative disorders, congenital disorders, infections, and inflammation. This review delves into the intricate interplay between glycosylation and protein conformation, with a specific focus on the profound impact of N-glycans on the selection of distinct protein conformations characterized by distinct interactomes-namely, protein assemblies-under normal and pathological conditions across various diseases. We begin by examining the spike protein of the SARS virus, illustrating how N-glycans regulate the infectivity of pathogenic agents. Subsequently, we utilize the prion protein and the chaperone glucose-regulated protein 94 as examples, exploring instances where N-glycosylation transforms physiological protein structures into disease-associated forms. Unraveling these connections provides valuable insights into potential therapeutic avenues and a deeper comprehension of the molecular intricacies that underlie disease conditions. This exploration of glycosylation's influence on protein conformation effectively bridges the gap between the glycome and disease, offering a comprehensive perspective on the therapeutic implications of targeting conformational mutants and their pathologic assemblies in various diseases. The goal is to unravel the nuances of these post-translational modifications, shedding light on how they contribute to the intricate interplay between protein conformation, assembly, and disease.

Development of a Data-Driven Integrative Model of a Bacterial Chromosome

The chromosome of archetypal bacteria E. coli is known for a complex topology with a 4.6 × 106 base pairs (bp) long sequence of nucleotides packed within a micrometer-sized cellular confinement. The inherent organization underlying this chromosome eludes general consensus due to the lack of a high-resolution picture of its conformation. Here we present our development of an integrative model of E. coli at a 500 bp resolution (https://github.com/JMLab-tifrh/ecoli_finer), which optimally combines a set of multiresolution genome-wide experimentally measured data within a framework of polymer based architecture. In particular the model is informed with an intragenome contact probability map at 5000 bp resolution derived via the Hi-C experiment and RNA-sequencing data at 500 bp resolution. Via dynamical simulations, this data-driven polymer based model generates an appropriate conformational ensemble commensurate with chromosome architectures that E. coli adopts. As a key hallmark of the E. coli chromosome the model spontaneously self-organizes into a set of nonoverlapping macrodomains and suitably locates plectonemic loops near the cell membrane. As novel extensions, it predicts a contact probability map simulated at a higher resolution than precedent experiments and can demonstrate segregation of chromosomes in a partially replicating cell. Finally, the modular nature of the model helps us devise control simulations to quantify the individual role of key features in hierarchical organization of the bacterial chromosome.

Mechanisms of allostery at the viral surface through the eyes of molecular simulation

The outermost surface layer of any virus is formed by either a capsid shell or envelope. Such layers have traditionally been thought of as immovable structures, but it is becoming apparent that they cannot be viewed exclusively as static architectures protecting the viral genome. A limited number of proteins on the virion surface must perform a multitude of functions in order to orchestrate the viral life cycle, and allostery can regulate their structures at multiple levels of organization, spanning individual molecules, protomers, large oligomeric assemblies, or entire viral surfaces. Here, we review recent contributions from the molecular simulation field to viral surface allostery, with a particular focus on the trimeric spike glycoprotein emerging from the coronavirus surface, and the icosahedral flaviviral envelope complex. As emerging viral pathogens continue to pose a global threat, an improved understanding of viral dynamics and allosteric regulation will prove crucial in developing novel therapeutic strategies.

Computational approaches to investigate fluoride binding, selectivity and transport across the membrane

The use of molecular dynamics (MD) simulations to study biomolecular systems has proven reliable in elucidating atomic-level details of structure and function. In this chapter, MD simulations were used to uncover new insights into two phylogenetically unrelated bacterial fluoride (F-) exporters: the CLCF F-/H+ antiporter and the Fluc F- channel. The CLCF antiporter, a member of the broader CLC family, has previously revealed unique stoichiometry, anion-coordinating residues, and the absence of an internal glutamate crucial for proton import in the CLCs. Through MD simulations enhanced with umbrella sampling, we provide insights into the energetics and mechanism of the CLCF transport process, including its selectivity for F- over HF. In contrast, the Fluc F- channel presents a novel architecture as a dual topology dimer, featuring two pores for F- export and a central non-transported sodium ion. Using computational electrophysiology, we simulate the electrochemical gradient necessary for F- export in Fluc and reveal details about the coordination and hydration of both F- and the central sodium ion. The procedures described here delineate the specifics of these advanced techniques and can also be adapted to investigate other membrane protein systems.

Evolving spike-protein N-glycosylation in SARS-CoV-2 variants

Since >3 years, SARS-CoV-2 has plunged humans into a colossal pandemic. Henceforth, multiple waves of infection have swept through the human population, led by variants that were able to partially evade acquired immunity. The co-evolution of SARS-CoV-2 variants with human immunity provides an excellent opportunity to study the interaction between viral pathogens and their human hosts. The heavily N-glycosylated spike-protein of SARS-CoV-2 plays a pivotal role in initiating infection and is the target for host immune-response, both of which are impacted by host-installed N-glycans. Using highly-sensitive DeGlyPHER approach, we compared the N-glycan landscape on spikes of the SARS-CoV-2 Wuhan-Hu-1 strain to seven WHO-defined variants of concern/interest, using recombinantly expressed, soluble spike-protein trimers, sharing same stabilizing-mutations. We found that N-glycan processing is conserved at most sites. However, in multiple variants, processing of N-glycans from high mannose- to complex-type is reduced at sites N165, N343 and N616, implicated in spike-protein function.

Structural insights into the modulation of coronavirus spike tilting and infectivity by hinge glycans

Coronavirus spike glycoproteins presented on the virion surface mediate receptor binding, and membrane fusion during virus entry and constitute the primary target for vaccine and drug development. How the structure dynamics of the full-length spikes incorporated in viral lipid envelope correlates with the virus infectivity remains poorly understood. Here we present structures and distributions of native spike conformations on vitrified human coronavirus NL63 (HCoV-NL63) virions without chemical fixation by cryogenic electron tomography (cryoET) and subtomogram averaging, along with site-specific glycan composition and occupancy determined by mass spectrometry. The higher oligomannose glycan shield on HCoV-NL63 spikes than on SARS-CoV-2 spikes correlates with stronger immune evasion of HCoV-NL63. Incorporation of cryoET-derived native spike conformations into all-atom molecular dynamic simulations elucidate the conformational landscape of the glycosylated, full-length spike that reveals a role of hinge glycans in modulating spike bending. We show that glycosylation at N1242 at the upper portion of the stalk is responsible for the extensive orientational freedom of the spike crown. Subsequent infectivity assays implicated involvement of N1242-glyan in virus entry. Our results suggest a potential therapeutic target site for HCoV-NL63.

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.

Structure adaptation in Omicron SARS-CoV-2/hACE2: Biophysical origins of evolutionary driving forces

Since its emergence, the COVID-19 threat has been sustained by a series of transmission waves initiated by new variants of the SARS-CoV-2 virus. Some of these arise with higher transmissivity and/or increased disease severity. Here, we use molecular dynamics simulations to examine the modulation of the fundamental interactions between the receptor binding domain (RBD) of the spike glycoprotein and the host cell receptor (human angiotensin-converting enzyme 2 [hACE2]) arising from Omicron variant mutations (BA.1 and BA.2) relative to the original wild-type strain. Our key findings are that glycans play a vital role at the RBD···hACE2 interface for the Omicrons, and the interplay between glycans and sequence mutations leads to enhanced binding. We find significant structural differences in the complexes, which overall bring the spike protein and its receptor into closer proximity. These are consistent with and attributed to the higher positive charge on the RBD conferred by BA.1 and BA.2 mutations relative to the wild-type. However, further differences between subvariants BA.1 and BA.2 (which have equivalent RBD charges) are also evident: mutations reduce interdomain interactions between the up chain and its clockwise neighbor chain in particular for the latter, resulting in enhanced flexibility for BA.2. Consequently, we see occurrence of additional close contacts in one replica of BA.2, which include binding to hACE2 by a second RBD in addition to the up chain. Although this motif is not seen in BA.1, we find that the Omicrons can directly/indirectly bind a down-RBD to hACE2 through glycans: the role of the glycan on N90 of hACE2 switches from inhibiting to facilitating the binding to Omicron spike protein via glycan-protein lateral interactions. These structural and electrostatic differences offer further insight into the mechanisms by which viral mutations modulate host cell binding and provide a biophysical basis for evolutionary driving forces.

Free Energy Simulations of Receptor-Binding Domain Opening of the SARS-CoV-2 Spike Indicate a Barrierless Transition with Slow Conformational Motions

Infection by sarbecoviruses begins with the attachment of the homotrimeric viral "spike" protein to the angiotensin-converting enzyme 2 receptor on the surface of mammalian cells. This requires one or more receptor-binding domains (RBDs) to be in the open (up) position. Here, we present the results of long molecular dynamics simulations with umbrella sampling (US) to compute a one-dimensional free energy profile of RBD opening/closing and the associated transition times. After ≃3.58μs of simulation time per US window (∼229 μs in total), which was required to approach trajectory decorrelation, the computed free energy profile was found to be without large barriers. This suggests that the RBD diffuses between the open and closed positions without significant energetic hindrance. This interpretation appears consistent with experiments but is at odds with some previous simulations. Modeling the RBD motion as diffusive dynamics along the computed free energy profile, we find that the overall time required for the transition is only about 2 μs, which is 5 orders of magnitude shorter than experimentally measured transition times. We speculate that the most likely reason for the transition time mismatch is our use of very short glycans, which was required to make the simulations performed here feasible. Despite the long simulation times, the final free energy profile is not fully converged with statistical errors of ≃1.16 kcal/mol, which were found to be consistent with the slow time decay in the autocorrelation of the conformational motions of the protein. The simulation lengths that would be required to obtain fully converged results remain unknown, but the present calculations would benefit from at least an order-of-magnitude extension.

Intrinsically Disordered Regions Function as a Cervical Collar to Remotely Regulate the Nodding Dynamics of SARS-CoV-2 Prefusion Spike Heads

The SARS-CoV-2 prefusion spike heads (receptor binding domains, RBDs) frequently nod down and up to interact with host cell receptors. As the spike protein is a trimeric unit of significant size, to understand its large-scale structural dynamics associated with the nodding mechanism and the mutational impact on the same, we develop a topological symmetry-information-loaded coarse-grained structure-based model of a spike trimer using recent cryo-EM structural data. Our study reveals the control of two distant intrinsically disordered regions (IDRs), namely, 630 and FPPR loops, over the nodding dynamics of spike heads. We find that the order-disorder transition of IDRs becomes more evident in the variants of concern (VOCs) that are associated with the characteristic mutation, D614G, in the proximity of these IDRs. In some VOCs, the two other mutations A570D and S982A also show an integral effect. The driver mutation D614G instigates a salt-bridge disruption, altering the order-disorder dynamics of both 630 and FPPR loops and their interaction with the C-terminal domains (CTD1/CTD2). This altered connectivity in these mutants allows the two IDRs to act collectively as a "cervical collar" for the RBD, supporting various spike head postures, consistent with cryo-EM results available for specific cases. The IDRs' control over the spike structure and dynamics presents an exciting opportunity where they can be targeted as remote operational switches to artificially maneuver the nod for effective therapeutic interventions.

Interaction of SARS-CoV-2 with host cells and antibodies: experiment and simulation

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of the devastating global COVID-19 pandemic announced by WHO in March 2020. Through unprecedented scientific effort, several vaccines, drugs and antibodies have been developed, saving millions of lives, but the fight against COVID-19 continues as immune escape variants of concern such as Delta and Omicron emerge. To develop more effective treatments and to elucidate the side effects caused by vaccines and therapeutic agents, a deeper understanding of the molecular interactions of SARS-CoV-2 with them and human cells is required. With special interest in computational approaches, we will focus on the structure of SARS-CoV-2 and the interaction of its spike protein with human angiotensin-converting enzyme-2 (ACE2) as a prime entry point of the virus into host cells. In addition, other possible viral receptors will be considered. The fusion of viral and human membranes and the interaction of the spike protein with antibodies and nanobodies will be discussed, as well as the effect of SARS-CoV-2 on protein synthesis in host cells.

COVID infection in 4 steps: Thermodynamic considerations reveal how viral mucosal diffusion, target receptor affinity and furin cleavage act in concert to drive the nature and degree of infection in human COVID-19 disease

We have developed a mechanistic model of SARS-CoV-2 and SARS-CoV infection, exploring the relationship between the viral diffusion in the mucosa and viral affinity for the angiotensin converting enzyme 2 (ACE2) target. Utilising the structural similarity of SARS-CoV and SARS-CoV-2 and a shared viral target receptor (ACE2), but a dramatic difference in upper or lower respiratory tract infectivity, we were able to generate insights into the linkage of mucosal diffusion and target receptor affinity in determining the pathophysiological pathways of these two viruses. Our analysis reveals that for SARS-CoV-2 the higher affinity of ACE2 binding, the faster and more complete the mucosal diffusion in its transport from the upper airway to the region of the ACE2 target on the epithelium. This diffusional process is essential for the presentation of this virus to the furin catalysed highly efficient entry and infection process in the upper respiratory tract epithelial cells. A failure of SARS-CoV to follow this path is associated with lower respiratory tract infection and decreased infectivity. Thus, our analysis supports the view that through tropism SARS-CoV-2 has evolved a highly efficient membrane entry process that can act in concert with a high binding affinity of this virus and its variants for its ACE2 which in turn promotes enhanced movement of the virus from airway to epithelium. In this way ongoing mutations yielding higher affinities of SARS-CoV-2 for the ACE2 target becomes the basis for higher upper respiratory tract infectivity and greater viral spread. It is concluded that SARS-CoV-2 is constrained in the extent of its activities by the fundamental laws of physics and thermodynamics. Laws that describe diffusion and molecular binding. Moreover it can be speculated that the very earliest contact of this virus with the human mucosa defines the pathogenesis of this infection.

The SARS-CoV-2 spike protein structure: a symmetry tale on distortion trail

A preliminary step in the SARS-CoV-2 human infection process is a conformational change of the receptor binding domain (RBD) of its spike protein, characterized by a significant loss of symmetry. During this process, the residues which later on bind to the human angiotensin converting enzyme 2 (ACE2) receptor, become exposed at the surface of the protein. Symmetry analysis of a data set of 33 protein structures from experimental measurements and 32 structures from molecular dynamics simulation, show that the initial state carries clear indications on the structure of the final state, with respect to the local distortion along the sequence. This surprising finding implies that this type of analysis predicts the mechanism of change. We further show that the level of local distortion at the initial state increases with variant's transmissibility, for the wild type (WT) along with past and present variants of concern (WT ∼ alpha < beta < delta < Omicron BA.1), in accordance with the trend of their evolutionary path. In other words, the initial structure of the variant which is most infectious is also the most distorted, making its path to the final state shorter. It has been claimed that the RBD migration of the spike protein is allosterically controlled. Our analysis provides a quantitative support to a major theorem in this respect - that information about an allosteric process is encoded in the structure itself, suggesting that the path of local distortion is related to an allosteric information network.

Exploring and Learning the Universe of Protein Allostery Using Artificial Intelligence Augmented Biophysical and Computational Approaches

Allosteric mechanisms are commonly employed regulatory tools used by proteins to orchestrate complex biochemical processes and control communications in cells. The quantitative understanding and characterization of allosteric molecular events are among major challenges in modern biology and require integration of innovative computational experimental approaches to obtain atomistic-level knowledge of the allosteric states, interactions, and dynamic conformational landscapes. The growing body of computational and experimental studies empowered by emerging artificial intelligence (AI) technologies has opened up new paradigms for exploring and learning the universe of protein allostery from first principles. In this review we analyze recent developments in high-throughput deep mutational scanning of allosteric protein functions; applications and latest adaptations of Alpha-fold structural prediction methods for studies of protein dynamics and allostery; new frontiers in integrating machine learning and enhanced sampling techniques for characterization of allostery; and recent advances in structural biology approaches for studies of allosteric systems. We also highlight recent computational and experimental studies of the SARS-CoV-2 spike (S) proteins revealing an important and often hidden role of allosteric regulation driving functional conformational changes, binding interactions with the host receptor, and mutational escape mechanisms of S proteins which are critical for viral infection. We conclude with a summary and outlook of future directions suggesting that AI-augmented biophysical and computer simulation approaches are beginning to transform studies of protein allostery toward systematic characterization of allosteric landscapes, hidden allosteric states, and mechanisms which may bring about a new revolution in molecular biology and drug discovery.

Integrated analyses reveal a hinge glycan regulates coronavirus spike tilting and virus infectivity

Coronavirus spike glycoproteins presented on the virion surface mediate receptor binding, and membrane fusion during virus entry and constitute the primary target for vaccine and drug development. How the structure dynamics of the full-length spikes incorporated in viral lipid envelope correlates with the virus infectivity remains poorly understood. Here we present structures and distributions of native spike conformations on vitrified human coronavirus NL63 (HCoV-NL63) virions without chemical fixation by cryogenic electron tomography (cryoET) and subtomogram averaging, along with site-specific glycan composition and occupancy determined by mass spectroscopy. The higher oligomannose glycan shield on HCoV-NL63 spikes than on SARS-CoV-2 spikes correlates with stronger immune evasion of HCoV-NL63. Incorporation of cryoET-derived native spike conformations into all-atom molecular dynamic simulations elucidate the conformational landscape of the glycosylated, full-length spike that reveals a novel role of stalk glycans in modulating spike bending. We show that glycosylation at N1242 at the upper portion of the stalk is responsible for the extensive orientational freedom of the spike crown. Subsequent infectivity assays support the hypothesis that this glycan-dependent motion impacts virus entry. Our results suggest a potential therapeutic target site for HCoV-NL63.

Novel coronavirus mutations: Vaccine development and challenges

The ongoing global pandemic of novel coronavirus pneumonia (COVID-19) caused by the SARS-CoV-2 has a significant impact on global health and economy system. In this context, there have been some landmark advances in vaccine development. Over 100 new coronavirus vaccine candidates have been approved for clinical trials, with ten WHO-approved vaccines including four inactivated virus vaccines, two mRNA vaccines, three recombinant viral vectored vaccines and one protein subunit vaccine on the "Emergency Use Listing". Although the SARS-CoV-2 has an internal proofreading mechanism, there have been a number of mutations emerged in the pandemic affecting its transmissibility, pathogenicity and immunogenicity. Of these, mutations in the spike (S) protein and the resultant mutant variants have posed new challenges for vaccine development and application. In this review article, we present an overview of vaccine development, the prevalence of new coronavirus variants and their impact on protective efficacy of existing vaccines and possible immunization strategies coping with the viral mutation and diversity.

SARS-CoV-2 spike opening dynamics and energetics reveal the individual roles of glycans and their collective impact

The trimeric spike (S) glycoprotein, which protrudes from the SARS-CoV-2 viral envelope, binds to human ACE2, initiated by at least one protomer's receptor binding domain (RBD) switching from a "down" (closed) to an "up" (open) state. Here, we used large-scale molecular dynamics simulations and two-dimensional replica exchange umbrella sampling calculations with more than a thousand windows and an aggregate total of 160 μs of simulation to investigate this transition with and without glycans. We find that the glycosylated spike has a higher barrier to opening and also energetically favors the down state over the up state. Analysis of the S-protein opening pathway reveals that glycans at N165 and N122 interfere with hydrogen bonds between the RBD and the N-terminal domain in the up state, while glycans at N165 and N343 can stabilize both the down and up states. Finally, we estimate how epitope exposure for several known antibodies changes along the opening path. We find that the BD-368-2 antibody's epitope is continuously exposed, explaining its high efficacy.

Energetics of Spike Protein Opening of SARS-CoV-1 and SARS-CoV-2 and Its Variants of Concern: Implications in Host Receptor Scanning and Transmission

The receptor binding domain(s) (RBD) of spike (S) proteins of SARS-CoV-1 and SARS-CoV-2 (severe acute respiratory syndrome coronavirus) undergoes closed to open transition to engage with host ACE2 receptors. In this study, using multi atomistic (equilibrium) and targeted (non-equilibrium) molecular dynamics simulations, we have compared energetics of RBD opening pathways in full-length (modeled from cryo-EM structures) S proteins of SARS-CoV-1 and SARS-CoV-2. Our data indicate that amino acid variations at the RBD interaction interface can culminate into distinct free energy landscapes of RBD opening in these S proteins. We further report that mutations in the S protein of SARS-CoV-2 variants of concern can reduce the protein-protein interaction affinity of RBD(s) with its neighboring domains and could favor its opening to access ACE2 receptors. The findings can also aid in predicting the impact of future mutations on the rate of S protein opening for rapid host receptor scanning.

Frustration-driven allosteric regulation and signal transmission in the SARS-CoV-2 spike omicron trimer structures: a crosstalk of the omicron mutation sites allosterically regulates tradeoffs of protein stability and conformational adaptability

Dissecting the regulatory principles underlying function and activity of the SARS-CoV-2 spike protein at the atomic level is of paramount importance for understanding the mechanisms of virus transmissibility and immune escape. In this work, we introduce a hierarchical computational approach for atomistic modeling of allosteric mechanisms in the SARS-CoV-2 Omicron spike proteins and present evidence of a frustration-based allostery as an important energetic driver of the conformational changes and spike activation. By examining conformational landscapes and the residue interaction networks in the SARS-CoV-2 Omicron spike protein structures, we have shown that the Omicron mutational sites are dynamically coupled and form a central engine of the allosterically regulated spike machinery that regulates the balance and tradeoffs between conformational plasticity, protein stability, and functional adaptability. We have found that the Omicron mutational sites at the inter-protomer regions form regulatory hotspot clusters that control functional transitions between the closed and open states. Through perturbation-based modeling of allosteric interaction networks and diffusion analysis of communications in the closed and open spike states, we have quantified the allosterically regulated activation mechanism and uncover specific regulatory roles of the Omicron mutations. Atomistic reconstruction of allosteric communication pathways and kinetic modeling using Markov transient analysis reveal that the Omicron mutations form the inter-protomer electrostatic bridges that operate as a network of coupled regulatory switches that could control global conformational changes and signal transmission in the spike protein. The results of this study have revealed distinct and yet complementary roles of the Omicron mutation sites as a network of hotspots that enable allosteric modulation of structural stability and conformational changes which are central for spike activation and virus transmissibility.

Posttranslational modifications optimize the ability of SARS-CoV-2 spike for effective interaction with host cell receptors

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike glycoprotein is the prime target for vaccines, diagnostics, and therapeutic antibodies against the virus. While anchored in the viral envelope, for effective virulence, the spike needs to maintain structural flexibility to recognize the host cell surface receptors and bind to them, a property that can heavily depend upon the dynamics of the unresolved domains, most prominently the stalk. Construction of the complete, membrane-bound spike model and the description of its dynamics are critical steps in understanding the inner working of this key element of the viral infection by SARS-CoV-2. Combining homology modeling, protein-protein docking, and molecular dynamics (MD) simulations, we have developed a full spike structure in a native membrane. Multimicrosecond MD simulations of this model, the longest known single trajectory of the full spike, reveal conformational dynamics employed by the protein to explore the surface of the host cell. In agreement with cryogenic electron microscopy (cryo-EM), three flexible hinges in the stalk allow for global conformational heterogeneity of spike in the fully glycosylated system mediated by glycan-glycan and glycan-lipid interactions. The dynamical range of the spike is considerably reduced in its nonglycosylated form, confining the area explored by the spike on the host cell surface. Furthermore, palmitoylation of the membrane domain amplifies the local curvature that may prime the fusion. We show that the identified hinge regions are highly conserved in SARS coronaviruses, highlighting their functional importance in enhancing viral infection, and thereby, provide points for discovery of alternative therapeutics against the virus.

Integrating Conformational Dynamics and Perturbation-Based Network Modeling for Mutational Profiling of Binding and Allostery in the SARS-CoV-2 Spike Variant Complexes with Antibodies: Balancing Local and Global Determinants of Mutational Escape Mechanisms

In this study, we combined all-atom MD simulations, the ensemble-based mutational scanning of protein stability and binding, and perturbation-based network profiling of allosteric interactions in the SARS-CoV-2 spike complexes with a panel of cross-reactive and ultra-potent single antibodies (B1-182.1 and A23-58.1) as well as antibody combinations (A19-61.1/B1-182.1 and A19-46.1/B1-182.1). Using this approach, we quantify the local and global effects of mutations in the complexes, identify protein stability centers, characterize binding energy hotspots, and predict the allosteric control points of long-range interactions and communications. Conformational dynamics and distance fluctuation analysis revealed the antibody-specific signatures of protein stability and flexibility of the spike complexes that can affect the pattern of mutational escape. A network-based perturbation approach for mutational profiling of allosteric residue potentials revealed how antibody binding can modulate allosteric interactions and identified allosteric control points that can form vulnerable sites for mutational escape. The results show that the protein stability and binding energetics of the SARS-CoV-2 spike complexes with the panel of ultrapotent antibodies are tolerant to the effect of Omicron mutations, which may be related to their neutralization efficiency. By employing an integrated analysis of conformational dynamics, binding energetics, and allosteric interactions, we found that the antibodies that neutralize the Omicron spike variant mediate the dominant binding energy hotpots in the conserved stability centers and allosteric control points in which mutations may be restricted by the requirements of the protein folding stability and binding to the host receptor. This study suggested a mechanism in which the patterns of escape mutants for the ultrapotent antibodies may not be solely determined by the binding interaction changes but are associated with the balance and tradeoffs of multiple local and global factors, including protein stability, binding affinity, and long-range interactions.

Ionization of D571 Is Coupled with SARS-CoV-2 Spike Up/Down Equilibrium Revealing the pH-Dependent Allosteric Mechanism of Receptor-Binding Domains

As a type I viral fusion protein, SARS-CoV-2 spike undergoes a pH-dependent switch to mediate the endosomal positioning of the receptor-binding domain to facilitate viral entry into cells and immune evasion. Gaps in our knowledge concerning the conformational transitions and key intramolecular motivations have hampered the development of effective therapeutics against the virus. To clarify the pH-sensitive elements on spike-gating the receptor-binding domain (RBD) opening and understand the details of the RBD opening transition, we performed microsecond-time scale constant pH molecular dynamics simulations in this study. We identified the deeply buried D571 with a clear pKa shift, suggesting a potential pH sensor, and showed the coupling of ionization of D571 with spike RBD-up/down equilibrium. We also computed the free-energy landscape for RBD opening and identified the crucial interactions that influence RBD dynamics. The atomic-level characterization of the pH-dependent spike activation mechanism provided herein offers new insights for a better understanding of the fundamental mechanisms of SARS-CoV-2 viral entry and infection and hence supports the discovery of novel therapeutics for COVID-19.

Molecular simulations of complex carbohydrates and glycoconjugates

Complex carbohydrates (glycans) are the most abundant and versatile biopolymers in nature. The broad diversity of biochemical functions that carbohydrates cover is a direct consequence of the variety of 3D architectures they can adopt, displaying branched or linear arrangements, widely ranging in sizes, and with the highest diversity of building blocks of any other natural biopolymer. Despite this unparalleled complexity, a common denominator can be found in the glycans' inherent flexibility, which hinders experimental characterization, but that can be addressed by high-performance computing (HPC)-based molecular simulations. In this short review, I present and discuss the state-of-the-art of molecular simulations of complex carbohydrates and glycoconjugates, highlighting methodological strengths and weaknesses, important insights through emblematic case studies, and suggesting perspectives for future developments.

Application of Coarse-Grained (CG) Models to Explore Conformational Pathway of Large-Scale Protein Machines

Protein machines are clusters of protein assemblies that function in order to control the transfer of matter and energy in cells. For a specific protein machine, its working mechanisms are not only determined by the static crystal structures, but also related to the conformational transition dynamics and the corresponding energy profiles. With the rapid development of crystallographic techniques, the spatial scale of resolved structures is reaching up to thousands of residues, and the concomitant conformational changes become more and more complicated, posing a great challenge for computational biology research. Previously, a coarse-grained (CG) model aiming at conformational free energy evaluation was developed and showed excellent ability to reproduce the energy profiles by accurate electrostatic interaction calculations. In this study, we extended the application of the CG model to a series of large-scale protein machine systems. The spike protein trimer of SARS-CoV-2, ATP citrate lyase (ACLY) tetramer, and P4-ATPases systems were carefully studied and discussed as examples. It is indicated that the CG model is effective to depict the energy profiles of the conformational pathway between two endpoint structures, especially for large-scale systems. Both the energy change and energy barrier between endpoint structures provide reasonable mechanism explanations for the associated biological processes, including the opening of receptor binding domain (RBD) of spike protein, the phospholipid transportation of P4-ATPase, and the loop translocation of ACLY. Taken together, the CG model provides a suitable alternative in mechanistic studies related to conformational change in large-scale protein machines.

Computational design of stapled peptide inhibitor against SARS‐CoV ‐2 receptor binding domain

Since its first detection in 2019, the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS‐CoV‐2) has been the cause of millions of deaths worldwide. Despite the development and administration of different vaccines, the situation is still worrisome as the virus is constantly mutating to produce newer variants some of which are highly infectious. This raises an urgent requirement to understand the infection mechanism and thereby design therapeutic‐based treatment for COVID‐19. The gateway of the virus to the host cell is mediated by the binding of the receptor binding domain (RBD) of the virus spike protein to the angiotensin‐converting enzyme 2 (ACE2) of the human cell. Therefore, the RBD of SARS‐CoV‐2 can be used as a target to design therapeutics. The α1 helix of ACE2, which forms direct contact with the RBD surface, has been used as a template in the current study to design stapled peptide therapeutics. Using computer simulation, the mechanism and thermodynamics of the binding of six stapled peptides with RBD have been estimated. Among these, the one with two lactam stapling agents has shown binding affinity, sufficient to overcome RBD‐ACE2 binding. Analyses of the mechanistic detail reveal that a reorganization of amino acids at the RBD‐ACE2 interface produces favorable enthalpy of binding whereas conformational restriction of the free peptide reduces the loss in entropy to result higher binding affinity. The understanding of the relation of the nature of the stapling agent with their binding affinity opens up the avenue to explore stapled peptides as therapeutic against SARS‐CoV‐2.

Computer Simulations and Network-Based Profiling of Binding and Allosteric Interactions of SARS-CoV-2 Spike Variant Complexes and the Host Receptor: Dissecting the Mechanistic Effects of the Delta and Omicron Mutations

In this study, we combine all-atom MD simulations and comprehensive mutational scanning of S-RBD complexes with the angiotensin-converting enzyme 2 (ACE2) host receptor in the native form as well as the S-RBD Delta and Omicron variants to (a) examine the differences in the dynamic signatures of the S-RBD complexes and (b) identify the critical binding hotspots and sensitivity of the mutational positions. We also examined the differences in allosteric interactions and communications in the S-RBD complexes for the Delta and Omicron variants. Through the perturbation-based scanning of the allosteric propensities of the SARS-CoV-2 S-RBD residues and dynamics-based network centrality and community analyses, we characterize the global mediating centers in the complexes and the nature of local stabilizing communities. We show that a constellation of mutational sites (G496S, Q498R, N501Y and Y505H) correspond to key binding energy hotspots and also contribute decisively to the key interfacial communities that mediate allosteric communications between S-RBD and ACE2. These Omicron mutations are responsible for both favorable local binding interactions and long-range allosteric interactions, providing key functional centers that mediate the high transmissibility of the virus. At the same time, our results show that other mutational sites could provide a "flexible shield" surrounding the stable community network, thereby allowing the Omicron virus to modulate immune evasion at different epitopes, while protecting the integrity of binding and allosteric interactions in the RBD-ACE2 complexes. This study suggests that the SARS-CoV-2 S protein may exploit the plasticity of the RBD to generate escape mutants, while engaging a small group of functional hotspots to mediate efficient local binding interactions and long-range allosteric communications with ACE2.

Landscape-Based Protein Stability Analysis and Network Modeling of Multiple Conformational States of the SARS-CoV-2 Spike D614G Mutant: Conformational Plasticity and Frustration-Induced Allostery as Energetic Drivers of Highly Transmissible Spike Variants

The structural and functional studies of the SARS-CoV-2 spike protein variants revealed an important role of the D614G mutation that is shared across many variants of concern (VOCs), suggesting the effect of this mutation on the enhanced virus infectivity and transmissibility. The recent structural and biophysical studies provided important evidence about multiple conformational substates of the D614G spike protein. The development of a plausible mechanistic model that can explain the experimental observations from a more unified thermodynamic perspective is an important objective of the current work. In this study, we employed efficient and accurate coarse-grained simulations of multiple structural substates of the D614G spike trimers together with the ensemble-based mutational frustration analysis to characterize the dynamics signatures of the conformational landscapes. By combining the local frustration profiling of the conformational states with residue-based mutational scanning of protein stability and network analysis of allosteric interactions and communications, we determine the patterns of mutational sensitivity in the functional regions and sites of variants. We found that the D614G mutation may induce a considerable conformational adaptability of the open states in the SARS-CoV-2 spike protein without compromising the folding stability and integrity of the spike protein. The results suggest that the D614G mutant may employ a hinge-shift mechanism in which the dynamic couplings between the site of mutation and the interprotomer hinge modulate the interdomain interactions, global mobility change, and the increased stability of the open form. This study proposes that mutation-induced modulation of the conformational flexibility and energetic frustration at the interprotomer interfaces may serve as an efficient mechanism for allosteric regulation of the SARS-CoV-2 spike proteins.

Structural and Computational Studies of the SARS-CoV-2 Spike Protein Binding Mechanisms with Nanobodies: From Structure and Dynamics to Avidity-Driven Nanobody Engineering

Nanobodies provide important advantages over traditional antibodies, including their smaller size and robust biochemical properties such as high thermal stability, high solubility, and the ability to be bioengineered into novel multivalent, multi-specific, and high-affinity molecules, making them a class of emerging powerful therapies against SARS-CoV-2. Recent research efforts on the design, protein engineering, and structure-functional characterization of nanobodies and their binding with SARS-CoV-2 S proteins reflected a growing realization that nanobody combinations can exploit distinct binding epitopes and leverage the intrinsic plasticity of the conformational landscape for the SARS-CoV-2 S protein to produce efficient neutralizing and mutation resistant characteristics. Structural and computational studies have also been instrumental in quantifying the structure, dynamics, and energetics of the SARS-CoV-2 spike protein binding with nanobodies. In this review, a comprehensive analysis of the current structural, biophysical, and computational biology investigations of SARS-CoV-2 S proteins and their complexes with distinct classes of nanobodies targeting different binding sites is presented. The analysis of computational studies is supplemented by an in-depth examination of mutational scanning simulations and identification of binding energy hotspots for distinct nanobody classes. The review is focused on the analysis of mechanisms underlying synergistic binding of multivalent nanobodies that can be superior to single nanobodies and conventional nanobody cocktails in combating escape mutations by effectively leveraging binding avidity and allosteric cooperativity. We discuss how structural insights and protein engineering approaches together with computational biology tools can aid in the rational design of synergistic combinations that exhibit superior binding and neutralization characteristics owing to avidity-mediated mechanisms.

Energy landscape of the SARS-CoV-2 reveals extensive conformational heterogeneity

Cryo-electron microscopy (cryo-EM) has produced a number of structural models of the SARS-CoV-2 spike, already prompting biomedical outcomes. However, these reported models and their associated electrostatic potential maps represent an unknown admixture of conformations stemming from the underlying energy landscape of the spike protein. As with any protein, some of the spike's conformational motions are expected to be biophysically relevant, but cannot be interpreted only by static models. Using experimental cryo-EM images, we present the energy landscape of the glycosylated spike protein, and identify the diversity of low-energy conformations in the vicinity of its open (so called 1RBD-up) state. The resulting atomic refinement reveal global and local molecular rearrangements that cannot be inferred from an average 1RBD-up cryo-EM model. Here we report varied degrees of "openness" in global conformations of the 1RBD-up state, not revealed in the single-model interpretations of the density maps, together with conformations that overlap with the reported models. We discover how the glycan shield contributes to the stability of these low-energy conformations. Five out of six binding sites we analyzed, including those for engaging ACE2, therapeutic mini-proteins, linoleic acid, two different kinds of antibodies, switch conformations between their known apo- and holo-conformations, even when the global spike conformation is 1RBD-up. This apo-to-holo switching is reminiscent of a conformational preequilibrium. We found only one binding site, namely that of AB-C135 remains in apo state within all the sampled free energy-minimizing models, suggesting an induced fit mechanism for the docking of this antibody to the spike.

Allosteric Determinants of the SARS-CoV-2 Spike Protein Binding with Nanobodies: Examining Mechanisms of Mutational Escape and Sensitivity of the Omicron Variant

Structural and biochemical studies have recently revealed a range of rationally engineered nanobodies with efficient neutralizing capacity against the SARS-CoV-2 virus and resilience against mutational escape. In this study, we performed a comprehensive computational analysis of the SARS-CoV-2 spike trimer complexes with single nanobodies Nb6, VHH E, and complex with VHH E/VHH V nanobody combination. We combined coarse-grained and all-atom molecular simulations and collective dynamics analysis with binding free energy scanning, perturbation-response scanning, and network centrality analysis to examine mechanisms of nanobody-induced allosteric modulation and cooperativity in the SARS-CoV-2 spike trimer complexes with these nanobodies. By quantifying energetic and allosteric determinants of the SARS-CoV-2 spike protein binding with nanobodies, we also examined nanobody-induced modulation of escaping mutations and the effect of the Omicron variant on nanobody binding. The mutational scanning analysis supported the notion that E484A mutation can have a significant detrimental effect on nanobody binding and result in Omicron-induced escape from nanobody neutralization. Our findings showed that SARS-CoV-2 spike protein might exploit the plasticity of specific allosteric hotspots to generate escape mutants that alter response to binding without compromising activity. The network analysis supported these findings showing that VHH E/VHH V nanobody binding can induce long-range couplings between the cryptic binding epitope and ACE2-binding site through a broader ensemble of communication paths that is less dependent on specific mediating centers and therefore may be less sensitive to mutational perturbations of functional residues. The results suggest that binding affinity and long-range communications of the SARS-CoV-2 complexes with nanobodies can be determined by structurally stable regulatory centers and conformationally adaptable hotspots that are allosterically coupled and collectively control resilience to mutational escape.

Conformational Flexibility and Local Frustration in the Functional States of the SARS-CoV-2 Spike B.1.1.7 and B.1.351 Variants: Mutation-Induced Allosteric Modulation Mechanism of Functional Dynamics and Protein Stability

Structural and functional studies of the SARS-CoV-2 spike proteins have recently determined distinct functional states of the B.1.1.7 and B.1.351 spike variants, providing a molecular framework for understanding the mechanisms that link the effect of mutations with the enhanced virus infectivity and transmissibility. A detailed dynamic and energetic analysis of these variants was undertaken in the present work to quantify the effects of different mutations on functional conformational changes and stability of the SARS-CoV-2 spike protein. We employed the efficient and accurate coarse-grained (CG) simulations of multiple functional states of the D614G mutant, B.1.1.7 and B.1.351 spike variants to characterize conformational dynamics of the SARS-CoV-2 spike proteins and identify dynamic signatures of the functional regions that regulate transitions between the closed and open forms. By combining molecular simulations with full atomistic reconstruction of the trajectories and the ensemble-based mutational frustration analysis, we characterized how the intrinsic flexibility of specific spike regions can control functional conformational changes required for binding with the host-cell receptor. Using the residue-based mutational scanning of protein stability, we determined protein stability hotspots and identified potential energetic drivers favoring the receptor-accessible open spike states for the B.1.1.7 and B.1.351 spike variants. The results suggested that modulation of the energetic frustration at the inter-protomer interfaces can serve as a mechanism for allosteric couplings between mutational sites and the inter-protomer hinges of functional motions. The proposed mechanism of mutation-induced energetic frustration may result in greater adaptability and the emergence of multiple conformational states in the open form. This study suggested that SARS-CoV-2 B.1.1.7 and B.1.351 variants may leverage the intrinsic plasticity of functional regions in the spike protein for mutation-induced modulation of protein dynamics and allosteric regulation to control binding with the host cell receptor.

Supramolecular Organization of SARS-CoV and SARS-CoV-2 Virions Revealed by Coarse-Grained Models of Intact Virus Envelopes

The coronavirus disease 19 (COVID-19) pandemic is causing a global health crisis and has already caused a devastating societal and economic burden. The pathogen, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has a high sequence and architecture identity with SARS-CoV, but far more people have been infected by SARS-CoV-2. Here, combining the structural data from cryo-electron microscopy and structure prediction, we constructed bottom-up Martini coarse-grained models of intact SARS-CoV and SARS-CoV-2 envelopes. Microsecond molecular dynamics simulations were performed, allowing us to explore their dynamics and supramolecular organization. Both SARS-CoV and SARS-CoV-2 envelopes present a spherical morphology, with structural proteins forming multiple string-like islands in the membrane and clusters between the heads of spike proteins. Critical differences between the SARS-CoV and SARS-CoV-2 envelopes are the interaction pattern between the spike proteins and the flexibility of the spike proteins. Our models provide structural and dynamic insights into the SARS virus envelopes and could be used for further investigation, such as drug design and membrane fusion and fission processes.

Ultraprecise Antigen 10-in-1 Pool Testing by Multiantibodies Transistor Assay

Effective screening of infectious diseases requires a fast, cheap, and population-scale testing. Antigen pool testing can increase the test rate and shorten the screening time, thus being a valuable approach for epidemic prevention and control. However, the overall percent agreement (OPA) with polymerase chain reaction (PCR) is one-half to three-quarters, hampering it from being a comprehensive method, especially pool testing, beyond the gold-standard PCR. Here, a multiantibodies transistor assay is developed for sensitive and highly precise antigen pool testing. The multiantibodies capture SARS-CoV-2 spike S1 proteins with different configurations, resulting in an antigen-binding affinity down to 0.34 fM. The limit of detection reaches 3.5 × 10^-17 g mL^-1SARS-CoV-2 spike S1 protein in artificial saliva, 4-5 orders of magnitude lower than existing transistor sensors. The testing of 60 nasopharyngeal swabs exhibits ∼100% OPA with PCR within an average diagnoses time of 38.9 s. Owing to its highly precise feature, a portable integrated platform is fabricated, which achieves 10-in-1 pooled screening for high testing throughput. This work solves the long-standing problem of antigen pool testing, enabling it to be a valuable tool in precise diagnoses and population-wide screening of COVID-19 or other epidemics in the future.

Distant residues modulate conformational opening in SARS-CoV-2 spike protein

The novel coronavirus (SARS-CoV-2) pandemic resulted in the largest public health crisis in recent times. Significant drug design effort against SARS-CoV-2 is focused on the receptor-binding domain (RBD) of the spike protein, although this region is highly prone to mutations causing therapeutic resistance. We applied deep data analysis methods on all-atom molecular dynamics simulations to identify key non-RBD residues that play a crucial role in spike−receptor binding and infection. Because the non-RBD residues are typically conserved across multiple coronaviruses, they can be targeted by broad-spectrum antibodies and drugs to treat infections from new strains that might appear during future epidemics.

Infection by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) involves the attachment of the receptor-binding domain (RBD) of its spike proteins to the ACE2 receptors on the peripheral membrane of host cells. Binding is initiated by a down-to-up conformational change in the spike protein, the change that presents the RBD to the receptor. To date, computational and experimental studies that search for therapeutics have concentrated, for good reason, on the RBD. However, the RBD region is highly prone to mutations, and is therefore a hotspot for drug resistance. In contrast, we here focus on the correlations between the RBD and residues distant to it in the spike protein. This allows for a deeper understanding of the underlying molecular recognition events and prediction of the highest-effect key mutations in distant, allosteric sites, with implications for therapeutics. Also, these sites can appear in emerging mutants with possibly higher transmissibility and virulence, and preidentifying them can give clues for designing pan-coronavirus vaccines against future outbreaks. Our model, based on time-lagged independent component analysis (tICA) and protein graph connectivity network, is able to identify multiple residues that exhibit long-distance coupling with the RBD opening. Residues involved in the most ubiquitous D614G mutation and the A570D mutation of the highly contagious UK SARS-CoV-2 variant are predicted ab initio from our model. Conversely, broad-spectrum therapeutics like drugs and monoclonal antibodies can target these key distant-but-conserved regions of the spike protein.

SARS-CoV-2 Spike Protein Mutations and Escape from Antibodies: A Computational Model of Epitope Loss in Variants of Concern

The SARS-CoV-2 spike (S) protein is exposed on the viral surface and is the first point of contact between the virus and the host. For these reasons it represents the prime target for Covid-19 vaccines. In recent months, variants of this protein have started to emerge. Their ability to reduce or evade recognition by S-targeting antibodies poses a threat to immunological treatments and raises concerns for their consequences on vaccine efficacy. To develop a model able to predict the potential impact of S-protein mutations on antibody binding sites, we performed unbiased multi-microsecond molecular dynamics of several glycosylated S-protein variants and applied a straightforward structure-dynamics-energy based strategy to predict potential changes in immunogenic regions on each variant. We recover known epitopes on the reference D614G sequence. By comparing our results, obtained on isolated S-proteins in solution, to recently published data on antibody binding and reactivity in new S variants, we directly show that modifications in the S-protein consistently translate into the loss of potentially immunoreactive regions. Our findings can thus be qualitatively reconnected to the experimentally characterized decreased ability of some of the Abs elicited against the dominant S-sequence to recognize variants. While based on the study of SARS-CoV-2 spike variants, our computational epitope-prediction strategy is portable and could be applied to study immunoreactivity in mutants of proteins of interest whose structures have been characterized, helping the development/selection of vaccines and antibodies able to control emerging variants.

A glycan gate controls opening of the SARS-CoV-2 spike protein

SARS-CoV-2 infection is controlled by the opening of the spike protein receptor binding domain (RBD), which transitions from a glycan-shielded 'down' to an exposed 'up' state to bind the human angiotensin-converting enzyme 2 receptor and infect cells. While snapshots of the 'up' and 'down' states have been obtained by cryo-electron microscopy and cryo-electron tomagraphy, details of the RBD-opening transition evade experimental characterization. Here over 130 µs of weighted ensemble simulations of the fully glycosylated spike ectodomain allow us to characterize more than 300 continuous, kinetically unbiased RBD-opening pathways. Together with ManifoldEM analysis of cryo-electron microscopy data and biolayer interferometry experiments, we reveal a gating role for the N-glycan at position N343, which facilitates RBD opening. Residues D405, R408 and D427 also participate. The atomic-level characterization of the glycosylated spike activation mechanism provided herein represents a landmark study for ensemble pathway simulations and offers a foundation for understanding the fundamental mechanisms of SARS-CoV-2 viral entry and infection.

Plant-Produced Glycosylated and In Vivo Deglycosylated Receptor Binding Domain Proteins of SARS-CoV-2 Induce Potent Neutralizing Responses in Mice

The COVID-19 pandemic, caused by SARS-CoV-2, has rapidly spread to more than 222 countries and has put global public health at high risk. The world urgently needs cost-effective and safe SARS-CoV-2 vaccines, antiviral, and therapeutic drugs to control it. In this study, we engineered the receptor binding domain (RBD) of the SARS-CoV-2 spike (S) protein and produced it in the plant Nicotiana benthamiana in a glycosylated and deglycosylated form. Expression levels of both glycosylated (gRBD) and deglycosylated (dRBD) RBD were greater than 45 mg/kg fresh weight. The purification yields were 22 mg of pure protein/kg of plant biomass for gRBD and 20 mg for dRBD, which would be sufficient for commercialization of these vaccine candidates. The purified plant-produced RBD protein was recognized by an S protein-specific monoclonal antibody, demonstrating specific reactivity of the antibody to the plant-produced RBD proteins. The SARS-CoV-2 RBD showed specific binding to angiotensin converting enzyme 2 (ACE2), the SARS-CoV-2 receptor. In mice, the plant-produced RBD antigens elicited high titers of antibodies with a potent virus-neutralizing activity. To our knowledge, this is the first report demonstrating that mice immunized with plant-produced deglycosylated RBD form elicited high titer of RBD-specific antibodies with potent neutralizing activity against SARS-CoV-2 infection. Thus, obtained data support that plant-produced glycosylated and in vivo deglycosylated RBD antigens, developed in this study, are promising vaccine candidates for the prevention of COVID-19.