Developing a Fully Glycosylated Full-Length SARS-CoV-2 Spike Protein Model in a Viral Membrane

A Reflection on the Use of Molecular Simulation to Respond to SARS-CoV-2 Pandemic Threats

Molecular simulations play important roles in understanding the lifecycle of the SARS-CoV-2 virus and contribute to the design and development of antiviral agents and diagnostic tests for COVID. Here, we discuss the insights that such simulations have provided and the challenges involved, focusing on the SARS-CoV-2 main protease (Mpro) and the spike glycoprotein. Mpro is the leading target for antivirals, while the spike glycoprotein is the target for vaccine design. Finally, we reflect on lessons from this pandemic for the simulation community. Data sharing initiatives and collaborations across the international research community contributed to advancing knowledge and should be built on to help in future pandemics and other global challenges such as antimicrobial resistance.

Exploring the effects of N234 and N343 linked glycans to SARS CoV 2 spike protein pocket accessibility using Gaussian accelerated molecular dynamics simulations

The N234 and N343-linked glycans of the SARS-CoV 2 spike protein are known to stabilize the up-conformation of its receptor-binding domains (RBDs), enabling human angiotensin enzyme 2 (hACE2) receptor binding. However, the effect of spike-hACE2 binding on these important glycans remains poorly understood, and these changes could have implications in the development of drugs that inhibit viral entry. In this study, Gaussian accelerated molecular dynamics (GaMD) simulations of the hACE2-free and hACE2-bound spike protein are performed. Biophysical analyses were focused on the accessibility of three previously suggested druggable pockets underneath the three RBD subunits. A shielding effect by N234-linked glycans on the components of their adjacent pockets was observed. Although deshielding of central scaffold residues was observed in the hACE2-bound state, pocket A's accessibility was reduced due to an increase in NTDB-RBDB contacts, restricting entry into the pocket. For pocket B, changes in N234C and N343C expose the central scaffold residues in the bound state, increasing accessibility. In Pocket C, increased shielding due to N234A was found in the bound state, reducing accessibility. Despite these changes, the pockets remain accessible to ligands in both states and are still valid targets for drug development studies.

Quantitative Prediction of Protein-Polyelectrolyte Binding Thermodynamics: Adsorption of Heparin-Analog Polysulfates to the SARS-CoV-2 Spike Protein RBD

Interactions of polyelectrolytes (PEs) with proteins play a crucial role in numerous biological processes, such as the internalization of virus particles into host cells. Although docking, machine learning methods, and molecular dynamics (MD) simulations are utilized to estimate binding poses and binding free energies of small-molecule drugs to proteins, quantitative prediction of the binding thermodynamics of PE-based drugs presents a significant obstacle in computer-aided drug design. This is due to the sluggish dynamics of PEs caused by their size and strong charge-charge correlations. In this paper, we introduce advanced sampling methods based on a force-spectroscopy setup and theoretical modeling to overcome this barrier. We exemplify our method with explicit solvent all-atom MD simulations of the interactions between anionic PEs that show antiviral properties, namely heparin and linear polyglycerol sulfate (LPGS), and the SARS-CoV-2 spike protein receptor binding domain (RBD). Our prediction for the binding free-energy of LPGS to the wild-type RBD matches experimentally measured dissociation constants within thermal energy, k B T, and correctly reproduces the experimental PE-length dependence. We find that LPGS binds to the Delta-variant RBD with an additional free-energy gain of 2.4 k B T, compared to the wild-type RBD, due to the additional presence of two mutated cationic residues contributing to the electrostatic energy gain. We show that the LPGS-RBD binding is solvent dominated and enthalpy driven, though with a large entropy-enthalpy compensation. Our method is applicable to general polymer adsorption phenomena and predicts precise binding free energies and reconfigurational friction as needed for drug and drug-delivery design.

Bioinformatic Selection of Mannose-Specific Lectins from Allium genus as SARS-CoV-2 Inhibitors Analysing Protein-Protein Interaction

Mannose-specific lectins are carbohydrate-binding proteins known for their antiviral potential. This study uses a bioinformatic approach to investigate the possibility of lectins from Allium sativum (garlic) and Allium ursinum (wild garlic) as inhibitors of SARS-CoV-2 entry. The information spectrum method (ISM) identified key interaction frequencies between the SARS-CoV-2 spike protein and these lectins, explicitly targeting the receptor-binding domain (RBD) and glycosylated asparagine residues, including N234. Lectins from Allium species showed a high affinity for oligomannose-type glycans on the spike protein, potentially blocking virus entry by preventing the spike-ACE2 receptor interaction. We propose that Allium lectins are promising candidates for further experimental validation as SARS-CoV-2 inhibitors, offering potential therapeutic applications in managing viral infections.

Sequence of the SARS-CoV-2 Spike Transmembrane Domain Encodes Conformational Dynamics

The homotrimeric SARS-CoV-2 spike protein enables viral infection by undergoing a large conformational transition, which facilitates the fusion of the viral envelope with the host cell membrane. The spike protein is anchored to the SARS-CoV-2 envelope by its transmembrane domain (TMD), composed of three TM helices, each contributed by one of the protomers of spike. Although the TMD is known to be important for viral fusion, whether it is a passive anchor of the spike or actively promotes fusion remains unknown. Specifically, it is unclear if the TMD and its dynamics facilitate the prefusion to postfusion conformational transition of the spike. Here, we computationally study the dynamics and self-assembly of the SARS-CoV-2 spike TMD in homogeneous POPC and cholesterol containing membranes. Atomistic simulations of a long TM helix-containing protomer segment show that the membrane-embedded segment bobs, tilts and gains and loses helicity, locally thinning the membrane. Coarse-grained multimerization simulations using representative TM helix structures from the atomistic simulations exhibit diverse trimer populations whose architecture depends on the structure of the TM helix protomer. While a symmetric conformation reflects the symmetry of the resting spike, an asymmetric TMD conformation could promote membrane fusion through the stabilization of a fusion intermediate. Together, our simulations demonstrate that the sequence and length of the SARS-CoV-2 spike TM segment make it inherently dynamic, that trimerization does not abrogate these dynamics and that the various observed TMD conformations may enable viral fusion.

Phenotypic evolution of SARS-CoV-2 spike during the COVID-19 pandemic

SARS-CoV-2 variants are mainly defined by mutations in their spike. It is therefore critical to understand how the evolutionary trajectories of spike affect virus phenotypes. So far, it has been challenging to comprehensively compare the many spikes that emerged during the pandemic in a single experimental platform. Here we generated a panel of recombinant viruses carrying different spike proteins from 27 variants circulating between 2020 and 2024 in the same genomic background. We then assessed several of their phenotypic traits both in vitro and in vivo. We found distinct phenotypic trajectories of spike among and between variants circulating before and after the emergence of Omicron variants. Spike of post-Omicron variants maintained enhanced tropism for the nasal epithelium and large airways but displayed, over time, several phenotypic traits typical of the pre-Omicron variants. Hence, spike with phenotypic features of both pre- and post-Omicron variants may continue to emerge in the future.

CHARMM at 45: Enhancements in Accessibility, Functionality, and Speed

Since its inception nearly a half century ago, CHARMM has been playing a central role in computational biochemistry and biophysics. Commensurate with the developments in experimental research and advances in computer hardware, the range of methods and applicability of CHARMM have also grown. This review summarizes major developments that occurred after 2009 when the last review of CHARMM was published. They include the following: new faster simulation engines, accessible user interfaces for convenient workflows, and a vast array of simulation and analysis methods that encompass quantum mechanical, atomistic, and coarse-grained levels, as well as extensive coverage of force fields. In addition to providing the current snapshot of the CHARMM development, this review may serve as a starting point for exploring relevant theories and computational methods for tackling contemporary and emerging problems in biomolecular systems. CHARMM is freely available for academic and nonprofit research at https://academiccharmm.org/program.

Oligomerization-driven avidity correlates with SARS-CoV-2 cellular binding and inhibition

Cellular processes are controlled by the thermodynamics of the underlying biomolecular interactions. Frequently, structural investigations use one monomeric binding partner, while ensemble measurements of binding affinities generally yield one affinity representative of a 1:1 interaction, despite the majority of the proteome consisting of oligomeric proteins. For example, viral entry and inhibition in SARS-CoV-2 involve a trimeric spike surface protein, a dimeric angiotensin-converting enzyme 2 (ACE2) cell-surface receptor and dimeric antibodies. Here, we reveal that cooperativity correlates with infectivity and inhibition as opposed to 1:1 binding strength. We show that ACE2 oligomerizes spike more strongly for more infectious variants, while exhibiting weaker 1:1 affinity. Furthermore, we find that antibodies use induced oligomerization both as a primary inhibition mechanism and to enhance the effects of receptor-site blocking. Our results suggest that naive affinity measurements are poor predictors of potency, and introduce an antibody-based inhibition mechanism for oligomeric targets. More generally, they point toward a much broader role of induced oligomerization in controlling biomolecular interactions.

A Peptide Microarray Platform Approach for Discovery of Immunodominant Antibody Epitopes

Immunodominant epitope discovery platforms play an important role in identifying novel biomarkers for effective immunotherapies and diagnostics. Methods to analyze the B-cell repertoire have been improved both experimentally and computationally. We developed an enhanced peptide microarray platform to discover and subsequently screen immunodominant epitopes. We utilized SARS-Cov-2 IgG positive and negative samples as a proof-of-concept to demonstrate the power of these improved peptide microarrays. The method identified significantly discriminant epitopes that classify positive and negative samples with good performance both as single peptides and in combination. We provide the assay conditions and parameters that justify the use of peptide microarrays in the selection of high-affinity epitopes, and we directly compare peptide performance against proteins. The results suggest that this platform can be used to confidently identify immunodominant antiviral epitopes while also serving as a useful tool for high-volume screening.

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.

Comparative study of the interactions between C60 fullerene and SARS-CoV-2, HIV, eukaryotic, and bacterial model membranes: Insights into antimicrobial strategies with C60-peptide hybrids

The neutrophil-derived peptide, indolicidin, and the sphere-shaped carbon nanoparticle, C60, are contemporary components capable of acting as bactericides and virucides, among others. Herein, the coarse-grained molecular dynamics simulation method was used to simulate the interactions of gram-negative bacteria, eukaryotes, human immunodeficiency virus (HIV), and SARS-COV-2 membrane models with indolicidin, C60s, and C60-indolicidin hybrids. Our results demonstrated that the carbon nanoparticle penetrated all membrane models, except the bacterial membrane, which remained impenetrable to both the peptide and C60. Additionally, the membrane thickness did not change significantly. The peptide floated above the membranes, with only the side chains of the tryptophan (Trp)-rich site slightly permeating the membranes. After achieving stable contact between the membrane models and nanoparticles, the infiltrated C60s interacted with the unsaturated tail of phospholipids. The density results showed that C60s stayed close to indolicidin and continued to interact with it even after penetration. Indolicidin, especially its Trp-rich site, exhibited more contact with the head and tail of neutral phospholipids compared to other phospholipids. Moreover, both particles interacted with different kinds of glycosphingolipids located in the eukaryote membrane. This investigation has the potential to advance our knowledge of novel approaches to combat antimicrobial resistance.

Single domain antibodies from camelids in the treatment of microbial infections

Infectious diseases continue to pose significant global health challenges. In addition to the enduring burdens of ailments like malaria and HIV, the emergence of nosocomial outbreaks driven by antibiotic-resistant pathogens underscores the ongoing threats. Furthermore, recent infectious disease crises, exemplified by the Ebola and SARS-CoV-2 outbreaks, have intensified the pursuit of more effective and efficient diagnostic and therapeutic solutions. Among the promising options, antibodies have garnered significant attention due to their favorable structural characteristics and versatile applications. Notably, nanobodies (Nbs), the smallest functional single-domain antibodies of heavy-chain only antibodies produced by camelids, exhibit remarkable capabilities in stable antigen binding. They offer unique advantages such as ease of expression and modification and enhanced stability, as well as improved hydrophilicity compared to conventional antibody fragments (antigen-binding fragments (Fab) or single-chain variable fragments (scFv)) that can aggregate due to their low solubility. Nanobodies directly target antigen epitopes or can be engineered into multivalent Nbs and Nb-fusion proteins, expanding their therapeutic potential. This review is dedicated to charting the progress in Nb research, particularly those derived from camelids, and highlighting their diverse applications in treating infectious diseases, spanning both human and animal contexts.

Computational and experimental validation of phthalocyanine and hypericin as effective SARS-CoV-2 fusion inhibitors

Phthalocyanine and hypericin have been previously identified as possible SARS-CoV-2 Spike glycoprotein fusion inhibitors through a virtual screening procedure. In this paper, atomistic simulations of metal-free phthalocyanines and atomistic and coarse-grained simulations of hypericins, placed around a complete model of the Spike embedded in a viral membrane, allowed to further explore their multi-target inhibitory potential, uncovering their binding to key protein functional regions and their propensity to insert in the membrane. Following computational results, pre-treatment of a pseudovirus expressing the SARS-CoV-2 Spike protein with low compounds concentrations resulted in a strong inhibition of its entry into cells, suggesting the activity of these molecules should involve the direct targeting of the viral envelope surface. The combination of computational and in vitro results hence supports the role of hypericin and phthalocyanine as promising SARS-CoV-2 entry inhibitors, further endorsed by literature reporting the efficacy of these compounds in inhibiting SARS-CoV-2 activity and in treating hospitalized COVID-19 patients.Communicated by Ramaswamy H. Sarma.

SPANA: Spatial decomposition analysis for cellular-scale molecular dynamics simulations

The rapid increase in computational power with the latest supercomputers has enabled atomistic molecular dynamics (MDs) simulations of biomolecules in biological membrane, cytoplasm, and other cellular environments. These environments often contain a million or more atoms to be simulated simultaneously. Therefore, their trajectory analyses involve heavy computations that can become a bottleneck in the computational studies. Spatial decomposition analysis (SPANA) is a set of analysis tools in the Generalized-Ensemble Simulation System (GENESIS) software package that can carry out MD trajectory analyses of large-scale biological simulations using multiple CPU cores in parallel. SPANA applies the spatial decomposition of a large biological system to distribute structural and dynamical analyses into individual CPU cores, which reduces the computational time and the memory size, significantly. SPANA opens new possibilities for detailed atomistic analyses of biomacromolecules as well as solvent water molecules, ions, and metabolites in MD simulation trajectories of very large biological systems containing more than millions of atoms in cellular environments.

Inclusion of deuterated glycopeptides provides increased sequence coverage in hydrogen/deuterium exchange mass spectrometry analysis of SARS-CoV-2 spike glycoprotein

RATIONALE

Hydrogen/deuterium exchange mass spectrometry (HDX-MS) can provide precise analysis of a protein's conformational dynamics across varied states, such as heat-denatured versus native protein structures, localizing regions that are specifically affected by such conditional changes. Maximizing protein sequence coverage provides high confidence that regions of interest were located by HDX-MS, but one challenge for complete sequence coverage is N-glycosylation sites. The deuteration of peptides post-translationally modified by asparagine-bound glycans (glycopeptides) has not always been identified in previous reports of HDX-MS analyses, causing significant sequence coverage gaps in heavily glycosylated proteins and uncertainty in structural dynamics in many regions throughout a glycoprotein.

METHODS

We detected deuterated glycopeptides with a Tribrid Orbitrap Eclipse mass spectrometer performing data-dependent acquisition. An MS scan was used to identify precursor ions; if high-energy collision-induced dissociation MS/MS of the precursor indicated oxonium ions diagnostic for complex glycans, then electron transfer low-energy collision-induced dissociation MS/MS scans of the precursor identified the modified asparagine residue and the glycan's mass. As in traditional HDX-MS, the identified glycopeptides were then analyzed at the MS level in samples labeled with D2 O.

RESULTS

We report HDX-MS analysis of the SARS-CoV-2 spike protein ectodomain in its trimeric prefusion form, which has 22 predicted N-glycosylation sites per monomer, with and without heat treatment. We identified glycopeptides and calculated their average isotopic mass shifts from deuteration. Inclusion of the deuterated glycopeptides increased sequence coverage of spike ectodomain from 76% to 84%, demonstrated that glycopeptides had been deuterated, and improved confidence in results localizing structural rearrangements.

CONCLUSION

Inclusion of deuterated glycopeptides improves the analysis of the conformational dynamics of glycoproteins such as viral surface antigens and cellular receptors.

Cancer cells and viruses share common glycoepitopes: exciting opportunities toward combined treatments

Aberrant glycosylation patterns of glycoproteins and glycolipids have long been recognized as one the major hallmarks of cancer cells that has led to numerous glycoconjugate vaccine attempts. These abnormal glycosylation profiles mostly originate from the lack of key glycosyltransferases activities, mutations, over expressions, or modifications of the requisite chaperone for functional folding. Due to their relative structural simplicity, O-linked glycans of the altered mucin family of glycoproteins have been particularly attractive in the design of tumor associated carbohydrate-based vaccines. Several such glycoconjugate vaccine formulations have generated potent monoclonal anti-carbohydrate antibodies useful as diagnostic and immunotherapies in the fight against cancer. Paradoxically, glycoproteins related to enveloped viruses also express analogous N- and O-linked glycosylation patterns. However, due to the fact that viruses are not equipped with the appropriate glycosyl enzyme machinery, they need to hijack that of the infected host cells. Although the resulting N-linked glycans are very similar to those of normal cells, some of their O-linked glycan patterns often share the common structural simplicity to those identified on tumor cells. Consequently, given that both cancer cells and viral glycoproteins share both common N- and O-linked glycoepitopes, glycoconjugate vaccines could be highly attractive to generate potent immune responses to target both conditions.

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.

Multiscale Molecular Dynamics Simulations of the Homodimer Accessory Protein ORF7b of SARS-CoV-2

The SARS-CoV-2 ORF7b protein has drawn attention for its potential role in viral pathogenesis, but its structural details and lateral membrane associations remain elusive. In this study, we conducted multiscale molecular dynamics simulations to provide detailed molecular insights of the protein's dimerization, which is crucial for unraveling its structural model of protein-protein interface important to regulating cellular immune response. To gain a deeper understanding of homodimer configurations, we employed a machine learning algorithm for structural-based clustering. Clusters were categorized into three distinct groups for both parallel and antiparallel orientations, highlighting the influence of the initial monomer conformation on dimer configurations. Analysis of hydrogen bonding and π-π and π-cation stacking interactions within clusters revealed variations in interactions between clusters. In parallel dimers, weak stacking interactions in the transmembrane (TM) region were observed. In contrast, antiparallel dimers exhibited strong hydrogen bonding and stacking interactions contributing to tight dimeric packing, both within and outside the TM domain. Overall, our study provides a comprehensive view of the structural dynamics of ORF7b homodimerization in both parallel and antiparallel orientations. These findings shed light on the molecular interactions involved in ORF7b dimerization, which are crucial for understanding its potential roles in SARS-CoV-2 pathogenesis. This knowledge could inform future research and therapeutic strategies targeting this viral protein.

Assessing the interaction between hemoglobin and the receptor binding domain of SARS-CoV-2 spike protein through MARTINI coarse-grained molecular dynamics

The emergence of different coronavirus-related diseases in the 2000's (SARS, MERS, and Covid-19) warrants the need of a complete understanding of the pathological, biological, and biochemical behavior of this class of pathogens. Great attention has been paid to the SARS-CoV-2 Spike protein, and its interaction with the human ACE2 has been thoroughly investigated. Recent findings suggested that the SARS-CoV-2 components may interact with different human proteins, and hemoglobin has very recently been demonstrated as a potential target for the Spike protein. Here we have investigated the interaction between either adult or fetal hemoglobin and the receptor binding domain of the Spike protein at molecular level through advanced molecular dynamics techniques and proposed rational binding modes and energy estimations. Our results agree with biochemical data previously reported in literature. We also demonstrated that co-incubation of pulmonary epithelial cells with hemoglobin strongly reduces the pro-inflammatory effects exerted by the concomitant administration of Spike protein.

Unraveling viral drug targets: a deep learning-based approach for the identification of potential binding sites

The coronavirus disease 2019 (COVID-19) pandemic has spurred a wide range of approaches to control and combat the disease. However, selecting an effective antiviral drug target remains a time-consuming challenge. Computational methods offer a promising solution by efficiently reducing the number of candidates. In this study, we propose a structure- and deep learning-based approach that identifies vulnerable regions in viral proteins corresponding to drug binding sites. Our approach takes into account the protein dynamics, accessibility and mutability of the binding site and the putative mechanism of action of the drug. We applied this technique to validate drug targeting toward severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike glycoprotein S. Our findings reveal a conformation- and oligomer-specific glycan-free binding site proximal to the receptor binding domain. This site comprises topologically important amino acid residues. Molecular dynamics simulations of Spike in complex with candidate drug molecules bound to the potential binding sites indicate an equilibrium shifted toward the inactive conformation compared with drug-free simulations. Small molecules targeting this binding site have the potential to prevent the closed-to-open conformational transition of Spike, thereby allosterically inhibiting its interaction with human angiotensin-converting enzyme 2 receptor. Using a pseudotyped virus-based assay with a SARS-CoV-2 neutralizing antibody, we identified a set of hit compounds that exhibited inhibition at micromolar concentrations.

Inconspicuous Yet Indispensable: The Coronavirus Spike Transmembrane Domain

Membrane-spanning portions of proteins' polypeptide chains are commonly known as their transmembrane domains (TMDs). The structural organisation and dynamic behaviour of TMDs from proteins of various families, be that receptors, ion channels, enzymes etc., have been under scrutiny on the part of the scientific community for the last few decades. The reason for such attention is that, apart from their obvious role as an "anchor" in ensuring the correct orientation of the protein's extra-membrane domains (in most cases functionally important), TMDs often actively and directly contribute to the operation of "the protein machine". They are capable of transmitting signals across the membrane, interacting with adjacent TMDs and membrane-proximal domains, as well as with various ligands, etc. Structural data on TMD arrangement are still fragmentary at best due to their complex molecular organisation as, most commonly, dynamic oligomers, as well as due to the challenges related to experimental studies thereof. Inter alia, this is especially true for viral fusion proteins, which have been the focus of numerous studies for quite some time, but have provoked unprecedented interest in view of the SARS-CoV-2 pandemic. However, despite numerous structure-centred studies of the spike (S) protein effectuating target cell entry in coronaviruses, structural data on the TMD as part of the entire spike protein are still incomplete, whereas this segment is known to be crucial to the spike's fusogenic activity. Therefore, in attempting to bring together currently available data on the structure and dynamics of spike proteins' TMDs, the present review aims to tackle a highly pertinent task and contribute to a better understanding of the molecular mechanisms underlying virus-mediated fusion, also offering a rationale for the design of novel efficacious methods for the treatment of infectious diseases caused by SARS-CoV-2 and related viruses.

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.

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.

Experimental and in silico evaluations of the possible molecular interaction between airborne particulate matter and SARS-CoV-2

During the early stage of the COVID-19 pandemic (winter 2020), the northern part of Italy has been significantly affected by viral infection compared to the rest of the country leading the scientific community to hypothesize that airborne particulate matter (PM) could act as a carrier for the SARS-CoV-2. To address this controversial issue, we first verified and demonstrated the presence of SARS-CoV-2 RNA genome on PM2.5 samples, collected in the city of Bologna (Northern Italy) in winter 2021. Then, we employed classical molecular dynamics (MD) simulations to investigate the possible recognition mechanism(s) between a newly modelled PM2.5 fragment and the SARS-CoV-2 Spike protein. The potential molecular interaction highlighted by MD simulations suggests that the glycans covering the upper Spike protein regions would mediate the direct contact with the PM2.5 carbon core surface, while a cloud of organic and inorganic PM2.5 components surround the glycoprotein with a network of non-bonded interactions resulting in up to 4769 total contacts. Moreover, a binding free energy of -207.2 ± 3.9 kcal/mol was calculated for the PM-Spike interface through the MM/GBSA method, and structural analyses also suggested that PM attachment does not alter the protein conformational dynamics. Although the association between the PM and SARS-CoV-2 appears plausible, this simulation does not assess whether these established interactions are sufficiently stable to carry the virus in the atmosphere, or whether the virion retains its infectiousness after the transport. While these key aspects should be verified by further experimental analyses, for the first time, this pioneering study gains insights into the molecular interactions between PM and SARS-CoV-2 Spike protein and will support further research aiming at clarifying the possible relationship between PM abundance and the airborne diffusion of viruses.

Docking and Molecular Dynamics Simulations Clarify Binding Sites for Interactions of Novel Marine Sulfated Glycans with SARS-CoV-2 Spike Glycoprotein

The entry of SARS-CoV-2 into the host cell is mediated by its S-glycoprotein (SGP). Sulfated glycans bind to the SGP receptor-binding domain (RBD), which forms a ternary complex with its receptor angiotensin converting enzyme 2. Here, we have conducted a thorough and systematic computational study of the binding of four oligosaccharide building blocks from novel marine sulfated glycans (isolated from Pentacta pygmaea and Isostichopus badionotus) to the non-glycosylated and glycosylated RBD. Blind docking studies using three docking programs identified five potential cryptic binding sites. Extensive site-targeted docking and molecular dynamics simulations using two force fields confirmed only two binding sites (Sites 1 and 5) for these novel, highly charged sulfated glycans, which were also confirmed by previously published reports. This work showed the structural features and key interactions driving ligand binding. A previous study predicted Site 2 to be a potential binding site, which was not observed here. The use of several molecular modeling approaches gave a comprehensive assessment. The detailed comparative study utilizing multiple modeling approaches is the first of its kind for novel glycan-SGP interaction characterization. This study provided insights into the key structural features of these novel glycans as they are considered for development as potential therapeutics.

SARS-CoV-2 uses CD4 to infect T helper lymphocytes

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the agent of a major global outbreak of respiratory tract disease known as Coronavirus Disease 2019 (COVID-19). SARS-CoV-2 infects mainly lungs and may cause several immune-related complications, such as lymphocytopenia and cytokine storm, which are associated with the severity of the disease and predict mortality. The mechanism by which SARS-CoV-2 infection may result in immune system dysfunction is still not fully understood. Here, we show that SARS-CoV-2 infects human CD4⁺ T helper cells, but not CD8⁺ T cells, and is present in blood and bronchoalveolar lavage T helper cells of severe COVID-19 patients. We demonstrated that SARS-CoV-2 spike glycoprotein (S) directly binds to the CD4 molecule, which in turn mediates the entry of SARS- CoV-2 in T helper cells. This leads to impaired CD4 T cell function and may cause cell death. SARS-CoV-2-infected T helper cells express higher levels of IL-10, which is associated with viral persistence and disease severity. Thus, CD4-mediated SARS-CoV-2 infection of T helper cells may contribute to a poor immune response in COVID-19 patients.

Structural effects of spike protein D614G mutation in SARS-CoV-2

A single mutation from aspartate to glycine at position 614 has dominated all circulating variants of the severe acute respiratory syndrome coronavirus 2. D614G mutation induces structural changes in the spike (S) protein that strengthen the virus infectivity. Here, we use molecular dynamics simulations to dissect the effects of mutation and 630-loop rigidification on S-protein structure. The introduction of the mutation orders the 630-loop structure and thereby induces global structural changes toward the cryoelectron microscopy structure of the D614G S-protein. The ordered 630-loop weakens local interactions between the 614th residue and others in contrast to disordered structures in the wild-type protein. The mutation allosterically alters global interactions between receptor-binding domains, forming an asymmetric and mobile down conformation and facilitating transitions toward up conformation. The loss of salt bridge between D614 and K854 upon the mutation generally stabilizes S-protein protomer, including the fusion peptide proximal region that mediates membrane fusion. Understanding the molecular basis of D614G mutation is crucial as it dominates in all variants of concern, including Delta and Omicron.

The SARS-CoV-2 Spike Protein Mutation Explorer: using an interactive application to improve the public understanding of SARS-CoV-2 variants of concern

Due to the COVID-19 pandemic the virus responsible, SARS-CoV-2, became a source of intense interest for non-expert audiences. The viral spike protein gained particular public interest as the main target for protective immune responses, including those elicited by vaccines. The rapid evolution of SARS-CoV-2 resulted in variations in the spike that enhanced transmissibility or weakened vaccine protection. This created new variants of concern (VOCs). The emergence of VOCs was studied using viral sequence data which was shared through portals such as the online Mutation Explorer of the COVID-19 Genomics UK consortium (COG-UK/ME). This was designed for an expert audience, but the information it contained could be of general interest if suitably communicated. Visualisations, interactivity and animation can improve engagement and understanding of molecular biology topics, and so we developed a graphical educational resource, the SARS-CoV-2 Spike Protein Mutation Explorer (SSPME), which used interactive 3D molecular models and animations to explain the molecular biology underpinning VOCs. User testing showed that the SSPME had better usability and improved participant knowledge confidence and knowledge acquisition compared to COG-UK/ME. This demonstrates how interactive visualisations can be used for effective molecular biology communication, as well as improving the public understanding of SARS-CoV-2 VOCs.

CHARMM-GUI Membrane Builder: Past, Current, and Future Developments and Applications

Molecular dynamics simulations of membranes and membrane proteins serve as computational microscopes, revealing coordinated events at the membrane interface. As G protein-coupled receptors, ion channels, transporters, and membrane-bound enzymes are important drug targets, understanding their drug binding and action mechanisms in a realistic membrane becomes critical. Advances in materials science and physical chemistry further demand an atomistic understanding of lipid domains and interactions between materials and membranes. Despite a wide range of membrane simulation studies, generating a complex membrane assembly remains challenging. Here, we review the capability of CHARMM-GUI Membrane Builder in the context of emerging research demands, as well as the application examples from the CHARMM-GUI user community, including membrane biophysics, membrane protein drug-binding and dynamics, protein-lipid interactions, and nano-bio interface. We also provide our perspective on future Membrane Builder development.

Development of a charged model of the SARS-CoV-2 viral surface

A recent study provided experimental evidence of inactivation of viral activity after radio-frequency (RF) exposures in the 6-12 GHz band that was hypothesized to be caused by vibrations of an acoustic dipole mode in the virus that excited the viral membrane to failure. Here, we develop an atomic-scale molecular dynamics (MD) model of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) viral surface to estimate the electric fields necessary to rupture the viral membrane via dipole shaking of the virus. We computed the absorption spectrum of the system via unbiased MD simulations and found no particular strong absorption in the GHz band. We investigated the mechanical resiliency of the viral membrane by introducing uniaxial strains in the system and observed no pore formation in the membrane for strains up to 50%. Because the computed absorption spectrum was found to be essentially flat, and the strain required to break the viral membrane was >0.5, the field strength associated with rupture of the virus was greater than the dielectric breakdown value of air. Thus, RF disinfection of enveloped viruses would occur only once sufficient heat was transferred to the virus via a thermal mechanism and not by direct action (shaking) of the RF field oscillations on the viral membrane.

Investigation of the effects of N-Acetylglucosamine on the stability of the spike protein in SARS-CoV-2 by molecular dynamics simulations

A lot of effort has been made in developing vaccine and therapeutic agents against the SARS-CoV-2, concentrating on the Spike protein that binds angiotensin-converting enzyme 2 on human cells. Nowadays, some researches study the role of the N-linked glycans as potential targets for vaccines and new agents. Due to the flexibility and diversity of the N-linked glycans, in this work, we focus on the N-Acetylglucosamine moiety, which is the precursor of nearly all eukaryotic glycans. We performed molecular dynamics simulations to study the effects of the N-Acetylglucosamine on the stability of the spike glycoprotein in SARS-CoV-2. After a 100 ns of simulation on the spike proteins without and with the N-Acetylglucosamine molecules, we found that the presence of N-Acetylglucosamine increases the local stability in their vicinity; even though their effect on the full structure is negligible. Thus; it can be inferred that the N-Acetylglucosamine moieties can potentially affect the interaction of the S protein with the ACE2 receptor. We also found that the S1 domain is more flexible than the S2 domain. We propose which of the experimentally observed glycans found on the spike may be more functional than the others. Detailed understanding of glycans is key for the development of new therapeutic strategies.

Identification of Potential Lead Compounds Targeting Novel Druggable Cavity of SARS-CoV-2 Spike Trimer by Molecular Dynamics Simulations

The global pandemic of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has become an urgent public health problem. Spike (S) protein mediates the fusion between the virus and the host cell membranes, consequently emerging as an important target of drug design. The lack of comparisons of in situ full-length S homotrimer structures in different states hinders understanding the structures and revealing the function, thereby limiting the discovery and development of therapeutic agents. Here, the steady-state structures of the in situ full-length S trimer in closed and open states (Sclosed and Sopen) were modeled with the constraints of density maps, associated with the analysis of the dynamic structural differences. Subsequently, we identified various regions with structure and property differences as potential binding pockets for ligands that promote the formation of inactive trimeric protein complexes. By using virtual screening strategy and a newly defined druggable cavity, five ligands were screened with potential bioactivities. Then molecular dynamic (MD) simulations were performed on apo protein structures and ligand bound complexes to reveal the conformational changes upon ligand binding. Our simulation results revealed that sulforaphane (SFN), which has the best binding affinity, could inhibit the conformational changes of S homotrimer that would occur during the viral membrane fusion. Our results could aid in the understanding of the regulation mechanism of S trimer aggregation and the structure-activity relationship, facilitating the development of potential antiviral agents.

Structural remodeling of SARS-CoV-2 spike protein glycans reveals the regulatory roles in receptor-binding affinity

Glycans of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike protein are speculated to play functional roles in the infection processes as they extensively cover the protein surface and are highly conserved across the variants. The spike protein has been the principal target for vaccine and therapeutic development while the exact effects of its glycosylation remain elusive. Analytical reports have described the glycan heterogeneity of the spike protein. Subsequent molecular simulation studies provided a knowledge basis of the glycan functions. However, experimental data on the role of discrete glycoforms on the spike protein pathobiology remains scarce. Building an understanding of their roles in SARS-CoV-2 is important as we continue to develop effective medicines and vaccines to combat the disease. Herein, we used designed combinations of glycoengineering enzymes to simplify and control the glycosylation profile of the spike protein receptor-binding domain (RBD). Measurements of the receptor-binding affinity revealed opposite regulatory effects of the RBD glycans with and without sialylation, which presents a potential strategy for modulating the spike protein behaviors through glycoengineering. Moreover, we found that the reported anti-SARS-CoV-(2) antibody, S309, neutralizes the impact of different RBD glycoforms on the receptor-binding affinity. In combination with molecular dynamics simulation, this work reports the regulatory roles that glycosylation plays in the interaction between the viral spike protein and host receptor, providing new insights into the nature of SARS-CoV-2. Beyond this study, enzymatic glycan remodeling offers the opportunity to understand the fundamental role of specific glycoforms on glycoconjugates across molecular biology.

Molecular architecture and dynamics of SARS-CoV-2 envelope by integrative modeling

Despite tremendous efforts, the exact structure of SARS-CoV-2 and related betacoronaviruses remains elusive. SARS-CoV-2 envelope is a key structural component of the virion that encapsulates viral RNA. It is composed of three structural proteins, spike, membrane (M), and envelope, which interact with each other and with the lipids acquired from the host membranes. Here, we developed and applied an integrative multi-scale computational approach to model the envelope structure of SARS-CoV-2 with near atomistic detail, focusing on studying the dynamic nature and molecular interactions of its most abundant, but largely understudied, M protein. The molecular dynamics simulations allowed us to test the envelope stability under different configurations and revealed that the M dimers agglomerated into large, filament-like, macromolecular assemblies with distinct molecular patterns. These results are in good agreement with current experimental data, demonstrating a generic and versatile approach to model the structure of a virus de novo.

Omicron mutations increase interdomain interactions and reduce epitope exposure in the SARS-CoV-2 spike

Omicron BA.1 is a highly infectious variant of SARS-CoV-2 that carries more than thirty mutations on the spike protein in comparison to the Wuhan wild type (WT). Some of the Omicron mutations, located on the receptor-binding domain (RBD), are exposed to the surrounding solvent and are known to help evade immunity. However, the impact of buried mutations on the RBD conformations and on the mechanics of the spike opening is less evident. Here, we use all-atom molecular dynamics (MD) simulations with metadynamics to characterize the thermodynamic RBD-opening ensemble, identifying significant differences between WT and Omicron. Specifically, the Omicron mutations S371L, S373P, and S375F make more RBD interdomain contacts during the spike's opening. Moreover, Omicron takes longer to reach the transition state than WT. It stabilizes up-state conformations with fewer RBD epitopes exposed to the solvent, potentially favoring immune or antibody evasion.

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.

Towards Quantum-Chemical Level Calculations of SARS-CoV-2 Spike Protein Variants of Concern by First Principles Density Functional Theory

The spike protein (S-protein) is a crucial part of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), with its many domains responsible for binding, fusion, and host cell entry. In this review we use the density functional theory (DFT) calculations to analyze the atomic-scale interactions and investigate the consequences of mutations in S-protein domains. We specifically describe the key amino acids and functions of each domain, which are essential for structural stability as well as recognition and fusion processes with the host cell; in addition, we speculate on how mutations affect these properties. Such unprecedented large-scale ab initio calculations, with up to 5000 atoms in the system, are based on the novel concept of amino acid-amino acid-bond pair unit (AABPU) that allows for an alternative description of proteins, providing valuable information on partial charge, interatomic bonding and hydrogen bond (HB) formation. In general, our results show that the S-protein mutations for different variants foster an increased positive partial charge, alter the interatomic interactions, and disrupt the HB networks. We conclude by outlining a roadmap for future computational research of biomolecular virus-related systems.

Binding of human ACE2 and RBD of Omicron enhanced by unique interaction patterns among SARS-CoV-2 variants of concern

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the virus causing COVID-19, has continued to mutate and spread worldwide despite global vaccination efforts. In particular, the Omicron variant, first identified in South Africa in late November 2021, has become the dominant strain worldwide. Compared to the original strain identified in Wuhan, Omicron features 50 genetic mutations, with 15 mutations in the receptor-binding domain (RBD) of the spike protein, which binds to the human angiotensin-converting enzyme 2 (ACE2) receptor for viral entry. However, it is not completely understood how these mutations alter the interaction and binding strength between the Omicron RBD and ACE2. In this study, we used a combined steered molecular dynamics (SMD) simulation and experimental microscale thermophoresis (MST) approach to quantify the interaction between Omicron RBD and ACE2. We report that the Omicron brings an enhanced RBD-ACE2 interface through N501Y, Q498R, and T478K mutations; the changes further lead to unique interaction patterns, reminiscing the features of previously dominated variants, Alpha (N501Y) and Delta (L452R and T478K). Among the Q493K and Q493R, we report that Q493R shows stronger binding to ACE2 than Q493K due to increased interactions. Our MST data confirmed that the Omicron mutations in RBD are associated with a five-fold higher binding affinity to ACE2 compared to the RBD of the original strain. In conclusion, our results could help explain the Omicron variant's prevalence in human populations, as higher interaction forces or affinity for ACE2 likely promote greater viral binding and internalization, leading to increased infectivity.

Use of simplified models for theoretical prediction of the interactions between available antibodies and the receptor-binding domain of SARS-CoV-2 spike protein

The negative impact of infectious diseases like COVID-19 on public health and the global economy is evident. This pandemic represents a significant challenge for the scientific community to develop new practical analytical methods for accurately diagnosing emerging cases. Due to their selectivity and sensitivity, new methodologies based on antigen/antibody interactions to detect COVID-19 biomarkers are necessary. In this context, the theoretical, computational modeling reduces experimental efforts and saves resources for rational biosensor design. This study proposes using molecular dynamics to predict the interactions between the Receptor Binding Domain (RBD) of the SARS-CoV-2 spike protein simplified model and a set of highly characterized antibodies. The binding free energy of the antigen/antibody complexes was calculated for the simplified models and compared against the complete SARS-CoV-2 ectodomain to validate the methodology. The structural data derived from our molecular dynamics and end-point free energy calculations showed a positive correlation between both approximations, with a 0.82 Pearson correlation coefficient; t = 3.661, df = 3, p-value = 0.03522, with a 95% confident interval. Furthermore, we identified the interfacial residues that could generate covalent bonds with a specific chemical surface without perturbing the binding dynamics to develop highly sensitive and specific diagnostic devices. Communicated by Ramaswamy H. Sarma.

A new approach for extracting information from protein dynamics

Increased ability to predict protein structures is moving research focus towards understanding protein dynamics. A promising approach is to represent protein dynamics through networks and take advantage of well-developed methods from network science. Most studies build protein dynamics networks from correlation measures, an approach that only works under very specific conditions, instead of the more robust inverse approach. Thus, we apply the inverse approach to the dynamics of protein dihedral angles, a system of internal coordinates, to avoid structural alignment. Using the well-characterized adhesion protein, FimH, we show that our method identifies networks that are physically interpretable, robust, and relevant to the allosteric pathway sites. We further use our approach to detect dynamical differences, despite structural similarity, for Siglec-8 in the immune system, and the SARS-CoV-2 spike protein. Our study demonstrates that using the inverse approach to extract a network from protein dynamics yields important biophysical insights.

Revealing the Molecular Interactions between Human ACE2 and the Receptor Binding Domain of the SARS-CoV-2 Wild-Type, Alpha and Delta Variants

After a sudden and first spread of the pandemic caused by the novel SARS-CoV-2 (Severe Acute Respiratory Syndrome-Coronavirus 2) wild-type strain, mutants have emerged which have been associated with increased infectivity, inducing surges in the contagions. The first of the so-called variants of concerns, was firstly isolated in the United Kingdom and later renamed Alpha variant. Afterwards, in the middle of 2021, a new variant appeared called Delta. The latter is characterized by the presence of point mutations in the Spike protein of SARS-CoV-2, especially in the Receptor Binding Domain (RBD). When in its active conformation, the RBD can interact with the human receptor Angiotensin-Converting Enzyme 2 (ACE2) to allow the entry of the virions into cells. In this contribution, by using extended all-atom molecular dynamic simulations, complemented with machine learning post-processing, we analyze the changes in the molecular interaction network induced by these different strains in comparison with the wild-type. On one hand, although relevant variations are evidenced, only limited changes in the global stability indicators and in the flexibility profiles have been observed. On the other hand, key differences were obtained by tracking hydrophilic and hydrophobic molecular interactions, concerning both positioning at the ACE2/RBD interface and formation/disruption dynamic behavior.

Valproate-coenzyme A conjugate blocks opening of receptor binding domains in the spike trimer of SARS-CoV-2 through an allosteric mechanism

The receptor-binding domains (RBDs) of the SARS-CoV-2 spike trimer exhibit "up" and "down" conformations often targeted by neutralizing antibodies. Only in the "up" configuration can RBDs bind to the ACE2 receptor of the host cell and initiate the process of viral multiplication. Here, we identify a lead compound (3-oxo-valproate-coenzyme A conjugate or Val-CoA) that stabilizes the spike trimer with RBDs in the down conformation. Val-CoA interacts with three R408 residues, one from each RBD, which significantly reduces the inter-subunit R408-R408 distance by ∼ 13 Å and closes the central pore formed by the three RBDs. Experimental evidence is presented that R408 is part of a triggering mechanism that controls the prefusion to postfusion state transition of the spike trimer. By stabilizing the RBDs in the down configuration, this and other related compounds can likely attenuate viral transmission. The reported findings for binding of Val-CoA to the spike trimer suggest a new approach for the design of allosteric antiviral drugs that do not have to compete for specific virus-receptor interactions but instead hinder the conformational motion of viral membrane proteins essential for interaction with the host cell. Here, we introduce an approach to target the spike protein by identifying lead compounds that stabilize the RBDs in the trimeric "down" configuration. When these compounds trimerize monomeric RBD immunogens as co-immunogens, they could also induce new types of non-ACE2 blocking antibodies that prevent local cell-to-cell transmission of the virus, providing a novel approach for inhibition of SARS-CoV-2.

IRAA: A statistical tool for investigating a protein-protein interaction interface from multiple structures

Understanding protein-protein interactions (PPIs) is fundamental to infer how different molecular systems work. A major component to model molecular recognition is the buried surface area (BSA), that is, the area that becomes inaccessible to solvent upon complex formation. To date, many attempts tried to connect BSA to molecular recognition principles, and in particular, to the underlying binding affinity. However, the most popular approach to calculate BSA is to use a single (or in some cases few) bound structures, consequently neglecting a wealth of structural information of the interacting proteins derived from ensembles corresponding to their unbound and bound states. Moreover, the most popular method inherently assumes the component proteins to bind as rigid entities. To address the above shortcomings, we developed a Monte Carlo method-based Interface Residue Assessment Algorithm (IRAA), to calculate a combined distribution of BSA for a given complex. Further, we apply our algorithm to human ACE2 and SARS-CoV-2 Spike protein complex, a system of prime importance. Results show a much broader distribution of BSA compared to that obtained from only the bound structure or structures and extended residue members of the interface with implications to the underlying biomolecular recognition. We derive that specific interface residues of ACE2 and of S-protein are consistently highly flexible, whereas other residues systematically show minor conformational variations. In effect, IRAA facilitates the use of all available structural data for any biomolecular complex of interest, extracting quantitative parameters with statistical significance, thereby providing a deeper biophysical understanding of the molecular system under investigation.

Computational analysis of protein stability and allosteric interaction networks in distinct conformational forms of the SARS-CoV-2 spike D614G mutant: reconciling functional mechanisms through allosteric model of spike regulation

In this study, we used an integrative computational approach to examine molecular mechanisms underlying functional effects of the D614G mutation by exploring atomistic modeling of the SARS-CoV-2 spike proteins as allosteric regulatory machines. We combined coarse-grained simulations, protein stability and dynamic fluctuation communication analysis with network-based community analysis to examine structures of the native and mutant SARS-CoV-2 spike proteins in different functional states. Through distance fluctuations communication analysis, we probed stability and allosteric communication propensities of protein residues in the native and mutant SARS-CoV-2 spike proteins, providing evidence that the D614G mutation can enhance long-range signaling of the allosteric spike engine. By combining functional dynamics analysis and ensemble-based alanine scanning of the SARS-CoV-2 spike proteins we found that the D614G mutation can improve stability of the spike protein in both closed and open forms, but shifting thermodynamic preferences towards the open mutant form. Our results revealed that the D614G mutation can promote the increased number of stable communities and allosteric hub centers in the open form by reorganizing and enhancing the stability of the S1-S2 inter-domain interactions and restricting mobility of the S1 regions. This study provides atomistic-based view of allosteric communications in the SARS-CoV-2 spike proteins, suggesting that the D614G mutation can exert its primary effect through allosterically induced changes on stability and communications in the residue interaction networks.Communicated by Ramaswamy H. Sarma.

Is the Stalk of the SARS-CoV-2 Spike Protein Druggable?

The spike protein is key to SARS-CoV-2 high infectivity because it facilitates the receptor binding domain (RBD) encounter with ACE2. As targeting subunit S1 has not yet delivered an ACE2-binding inhibitor, we have assessed the druggability of the conserved segment of the spike protein stalk within subunit S2 by means of an integrated computational approach that combines the molecular docking of an optimized library of fragments with high-throughput molecular dynamics simulations. The high propensity of the spike protein to mutate in key regions that are responsible for the recognition of the human angiotensin-converting enzyme 2 (hACE2) or for the recognition of antibodies, has made subunit S1 of the spike protein difficult to target. Despite the inherent flexibility of the stalk region, our results suggest two hidden interhelical binding sites, whose accessibility is only partially hampered by glycan residues.

Computer Simulation of the Interaction between SARS-CoV-2 Spike Protein and the Surface of Coinage Metals

A prominent feature of the SARS-CoV-2 virus is the presence of a large glycoprotein spike protruding from the virus envelope. The spike determines the interaction of the virus with the environment and the host. Here, we used an all-atom molecular dynamics simulation method to investigate the interaction of up- and down-conformations of the S1 subunit of the SARS-CoV-2 spike with the (100) surface of Au, Ag, and Cu. Our results revealed that the spike protein is adsorbed onto the surface of these metals, with Cu being the metal with the highest interaction with the spike. In our simulations, we considered the spike protein in both its up-conformation S^up (one receptor binding domain exposed) and down-conformation S^down (no exposed receptor binding domain). We found that the affinity of the metals for the up-conformation was higher than their affinity for the down-conformation. The structural changes in the spike in the up-conformation were also larger than the changes in the down-conformation. Comparing the present results for metals with those obtained in our previous MD simulations of S^up with other materials (cellulose, graphite, and human skin models), we see that Au induces the highest structural change in S^up, larger than those obtained in our previous studies.

Mutational escape prevention by combination of four neutralizing antibodies that target RBD conserved regions and stem helix

New variants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) appear rapidly every few months. They have showed powerful adaptive ability to circumvent the immune system. To further understand SARS-CoV-2's adaptability so as to seek for strategies to mitigate the emergence of new variants, herein we investigated the viral adaptation in the presence of broadly neutralizing antibodies and their combinations. First, we selected four broadly neutralizing antibodies, including pan-sarbecovirus and pan-betacoronavirus neutralizing antibodies that recognize distinct conserved regions on receptor-binding domain (RBD) or conserved stem-helix region on S2 subunit. Through binding competition analysis, we demonstrated that they were capable of simultaneously binding. Thereafter, a replication-competent vesicular stomatitis virus pseudotyped with SARS-CoV-2 spike protein was employed to study the viral adaptation. Twenty consecutive passages of the virus under the selective pressure of individual antibodies or their combinations were performed. It was found that it was not hard for the virus to adapt to broadly neutralizing antibodies, even for pan-sarbecovirus and pan-betacoronavirus antibodies. The virus was more and more difficult to escape the combinations of two/three/four antibodies. In addition, mutations in the viral population revealed by high-throughput sequencing showed that under the selective pressure of three/four combinational antibodies, viral mutations were not prone to present in the highly conserved region across betacoronaviruses (stem-helix region), while this was not true under the selective pressure of single/two antibodies. Importantly, combining neutralizing antibodies targeting RBD conserved regions and stem helix synergistically prevented the emergence of escape mutations. These studies will guide future vaccine and therapeutic development efforts and provide a rationale for the design of RBD-stem helix tandem vaccine, which may help to impede the generation of novel variants.

Thirty years of molecular dynamics simulations on posttranslational modifications of proteins

Posttranslational modifications (PTMs) are an integral component to how cells respond to perturbation. While experimental advances have enabled improved PTM identification capabilities, the same throughput for characterizing how structural changes caused by PTMs equate to altered physiological function has not been maintained. In this Perspective, we cover the history of computational modeling and molecular dynamics simulations which have characterized the structural implications of PTMs. We distinguish results from different molecular dynamics studies based upon the timescales simulated and analysis approaches used for PTM characterization. Lastly, we offer insights into how opportunities for modern research efforts on in silico PTM characterization may proceed given current state-of-the-art computing capabilities and methodological advancements.

ACE2 N-glycosylation modulates interactions with SARS-CoV-2 spike protein in a site-specific manner

SARS-CoV-2 has evolved continuously and accumulated spike mutations with each variant having a different binding for the cellular ACE2 receptor. It is not known whether the interactions between such mutated spikes and ACE2 glycans are conserved among different variant lineages. Here, we focused on three ACE2 glycosylation sites (53, 90 and 322) that are geometrically close to spike binding sites and investigated the effect of their glycosylation pattern on spike affinity. These glycosylation deletions caused distinct site-specific changes in interactions with the spike and acted cooperatively. Of note, the particular interaction profiles were conserved between the SARS-CoV-2 parental virus and the variants of concern (VOCs) Delta and Omicron. Our study provides insights for a better understanding of the importance of ACE2 glycosylation on ACE2/SARS-CoV-2 spike interaction and guidance for further optimization of soluble ACE2 for therapeutic use.

Computational analysis of interior mutations of SARS-CoV-2 Spike protein suggest a balance of protein stability and S2: S1 separation propensity

SARS-CoV-2 variants often include surface mutations in the Spike protein that are important for viruses to recognize host receptors and evade antibody neutralization. The Spike protein also has mutations in the interior of the protein likely to affect the Spike protein S1 - S2 subunit's separation propensity, the most important of which is the D614G mutation. Remarkably, the Omicron variant contains a large number of internal mutations at the S2: S1 interface, which have not been investigated yet. In this study, we examined the effects of such interfacial mutations on the S2: S1 and subunit domain interactions and on the subunit's dissociation process. We found that the interaction with S2 is mainly contributed by the three encapsulation domains, named INT, ED1 and ED2 of S1, which are sandwiched between the S1 RBD and N-terminal NTD domain. We found that D614 is the strongest contributor for the S2: S1 interaction which is greatly weakened by the D614G mutation. Surprisingly, we found that, mutations T547K, H655Y, N764K, N856K, N969K, L981F in the Omicron variant largely enhance the S2: ED1 interaction, partially compensating the loss of S2: ED2 interaction due to the D614G mutation. Lastly, these results, together with biological considerations, allow us to suggest that in addition to the binding strength of between the RBD and ACE2, the stability of the Spike protein and the propensity of Spike protein S2: S1 separation are critical factors which likely exist in a balance for a particular infectivity and pathogenicity of the virus.