The Answer Lies in the Energy: How Simple Atomistic Molecular Dynamics Simulations May Hold the Key to Epitope Prediction on the Fully Glycosylated SARS-CoV-2 Spike Protein.
Ontology highlight
ABSTRACT:
SUBMITTER: Serapian SA
PROVIDER: S-EPMC7491317 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:The bilayer bending modulus (Kc) is one of the most important physical constants characterizing lipid membranes, but precisely measuring it is a challenge, both experimentally and computationally. Experimental measurements on chemically identical bilayers often differ depending upon the techniques employed, and robust simulation results have previously been limited to coarse-grained models (at varying levels of resolution). This Communication demonstrates the extraction of Kc from fully atomistic molecular dynamics simulations for three different single-component lipid bilayers (DPPC, DOPC, and DOPE). The results agree quantitatively with experiments that measure thermal shape fluctuations in giant unilamellar vesicles. Lipid tilt, twist, and compression moduli are also reported.
| S-EPMC4183605 | biostudies-literature
Project description:SARS-CoV-2 outbreaks worldwide caused COVID-19 pandemic, which is related to several million deaths. In particular, SARS-CoV-2 Spike (S) protein is a major biological target for COVID-19 vaccine design. Unfortunately, recent reports indicated that Spike (S) protein mutations can lead to antibody resistance. However, understanding the process is limited, especially at the atomic scale. The structural change of S protein and neutralizing antibody fragment (FAb) complexes was thus probed using molecular dynamics (MD) simulations. In particular, the backbone RMSD of the 501Y.V2 complex was significantly larger than that of the wild-type one implying a large structural change of the mutation system. Moreover, the mean of Rg, CCS, and SASA are almost the same when compared two complexes, but the distributions of these values are absolutely different. Furthermore, the free energy landscape of the complexes was significantly changed when the 501Y.V2 variant was induced. The binding pose between S protein and FAb was thus altered. The FAb-binding affinity to S protein was thus reduced due to revealing over steered-MD (SMD) simulations. The observation is in good agreement with the respective experiment that the 501Y.V2 SARS-CoV-2 variant can escape from neutralizing antibody (NAb).
| S-EPMC8769535 | biostudies-literature
Project description:The effect of surface chemistry on the adsorption characteristics of a fibronectin fragment (FNIII8⁻10) was investigated using fully atomistic molecular dynamics simulations. Model surfaces were constructed to replicate self-assembled monolayers terminated with methyl, hydroxyl, amine, and carboxyl moieties. It was found that adsorption of FNIII8⁻10 on charged surfaces is rapid, specific, and driven by electrostatic interactions, and that the anchoring residues are either polar uncharged or of opposing charge to that of the targeted surfaces. On charged surfaces the presence of a strongly bound layer of water molecules and ions hinders FNIII8⁻10 adsorption. In contrast, adsorption kinetics on uncharged surfaces are slow and non-specific, as they are driven by van der Waals interactions, and the anchoring residues are polar uncharged. Due to existence of a positively charged area around its cell-binding region, FNIII8⁻10 is available for subsequent cell binding when adsorbed on a positively charged surface, but not when adsorbed on a negatively charged surface. On uncharged surfaces, the availability of the fibronectin fragment's cell-binding region is not clearly distinguished because adsorption is much less specific.
| S-EPMC6275015 | biostudies-literature
Project description:The spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) presents a public health crisis, and the vaccines that can induce highly potent neutralizing antibodies are essential for ending the pandemic. The spike (S) protein on the viral envelope mediates human angiotensin-converting enzyme 2 (ACE2) binding and thus is the target of a variety of neutralizing antibodies. In this work, we built various S trimer-antibody complex structures on the basis of the fully glycosylated S protein models described in our previous work, and performed all-atom molecular dynamics simulations to get insight into the structural dynamics and interactions between S protein and antibodies. Investigation of the residues critical for S-antibody binding allows us to predict the potential influence of mutations in SARS-CoV-2 variants. Comparison of the glycan conformations between S-only and S-antibody systems reveals the roles of glycans in S-antibody binding. In addition, we explored the antibody binding modes, and the influences of antibody on the motion of S protein receptor binding domains. Overall, our analyses provide a better understanding of S-antibody interactions, and the simulation-based S-antibody interaction maps could be used to predict the influences of S mutation on S-antibody interactions, which will be useful for the development of vaccine and antibody-based therapy.
| S-EPMC8132224 | biostudies-literature
Project description:The spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) presents a public health crisis, and the vaccines that can induce highly potent neutralizing antibodies are essential for ending the pandemic. The spike (S) protein on the viral envelope mediates human angiotensin-converting enzyme 2 binding and thus is the target of a variety of neutralizing antibodies. In this work, we built various S trimer-antibody complex structures on the basis of the fully glycosylated S protein models described in our previous work and performed all-atom molecular dynamics simulations to gain insight into the structural dynamics and interactions between S protein and antibodies. Investigation of the residues critical for S-antibody binding allows us to predict the potential influence of mutations in SARS-CoV-2 variants. Comparison of the glycan conformations between S-only and S-antibody systems reveals the roles of glycans in S-antibody binding. In addition, we explored the antibody binding modes and the influences of antibody on the motion of S protein receptor binding domains. Overall, our analyses provide a better understanding of S-antibody interactions, and the simulation-based S-antibody interaction maps could be used to predict the influences of S mutation on S-antibody interactions, which will be useful for the development of vaccine and antibody-based therapy.
| S-EPMC8457324 | biostudies-literature
Project description:This technical study describes all-atom modeling and simulation of a fully-glycosylated full-length SARS-CoV-2 spike (S) protein in a viral membrane. First, starting from PDB:6VSB and 6VXX, full-length S protein structures were modeled using template-based modeling, de-novo protein structure prediction, and loop modeling techniques in GALAXY modeling suite. Then, using the recently-determined most occupied glycoforms, 22 N-glycans and 1 O-glycan of each monomer were modeled using Glycan Reader & Modeler in CHARMM-GUI. These fully-glycosylated full-length S protein model structures were assessed and further refined against the low-resolution data in their respective experimental maps using ISOLDE. We then used CHARMM-GUI Membrane Builder to place the S proteins in a viral membrane and performed all-atom molecular dynamics simulations. All structures are available in CHARMM-GUI COVID-19 Archive (http://www.charmm-gui.org/docs/archive/covid19), so researchers can use these models to carry out innovative and novel modeling and simulation research for the prevention and treatment of COVID-19.
| S-EPMC7263518 | biostudies-literature
Project description:This technical study describes all-atom modeling and simulation of a fully glycosylated full-length SARS-CoV-2 spike (S) protein in a viral membrane. First, starting from PDB: 6VSB and 6VXX, full-length S protein structures were modeled using template-based modeling, de-novo protein structure prediction, and loop modeling techniques in GALAXY modeling suite. Then, using the recently determined most occupied glycoforms, 22 N-glycans and 1 O-glycan of each monomer were modeled using Glycan Reader & Modeler in CHARMM-GUI. These fully glycosylated full-length S protein model structures were assessed and further refined against the low-resolution data in their respective experimental maps using ISOLDE. We then used CHARMM-GUI Membrane Builder to place the S proteins in a viral membrane and performed all-atom molecular dynamics simulations. All structures are available in CHARMM-GUI COVID-19 Archive (http://www.charmm-gui.org/docs/archive/covid19) so that researchers can use these models to carry out innovative and novel modeling and simulation research for the prevention and treatment of COVID-19.
| S-EPMC7341691 | biostudies-literature
Project description:Mixtures of long- and short-tail phosphatidylcholine lipids are known to self-assemble into a variety of aggregates combining flat bilayerlike and curved micellelike features, commonly called bicelles. Atomistic simulations of bilayer ribbons and perforated bilayers containing dimyristoylphosphatidylcholine (DMPC, di-C(14) tails) and dihexanoylphosphatidylcholine (DHPC, di-C(6) tails) have been carried out to investigate the partitioning of these components between flat and curved microenvironments and the stabilization of the bilayer edge by DHPC. To approach equilibrium partitioning of lipids on an achievable simulation timescale, configuration-bias Monte Carlo mutation moves were used to allow individual lipids to change tail length within a semigrand-canonical ensemble. Since acceptance probabilities for direct transitions between DMPC and DHPC were negligible, a third component with intermediate tail length (didecanoylphosphatidylcholine, di-C(10) tails) was included at a low concentration to serve as an intermediate for transitions between DMPC and DHPC. Strong enrichment of DHPC is seen at ribbon and pore edges, with an excess linear density of approximately 3 nm(-1). The simulation model yields estimates for the onset of edge stability with increasing bilayer DHPC content between 5% and 15% DHPC at 300 K and between 7% and 17% DHPC at 323 K, higher than experimental estimates. Local structure and composition at points of close contact between pores suggest a possible mechanism for effective attractions between pores, providing a rationalization for the tendency of bicelle mixtures to aggregate into perforated vesicles and perforated sheets.
| S-EPMC2884235 | biostudies-other
Project description:The celebrated Miller experiments reported on the spontaneous formation of amino acids from a mixture of simple molecules reacting under an electric discharge, giving birth to the research field of prebiotic chemistry. However, the chemical reactions involved in those experiments have never been studied at the atomic level. Here we report on, to our knowledge, the first ab initio computer simulations of Miller-like experiments in the condensed phase. Our study, based on the recent method of treatment of aqueous systems under electric fields and on metadynamics analysis of chemical reactions, shows that glycine spontaneously forms from mixtures of simple molecules once an electric field is switched on and identifies formic acid and formamide as key intermediate products of the early steps of the Miller reactions, and the crucible of formation of complex biological molecules.
| S-EPMC4183268 | biostudies-literature
Project description:Connexins are plasma membrane proteins that associate in hexameric complexes to form channels named connexons. Two connexons in neighboring cells may dock to form a "gap junction" channel, i.e. an intercellular conduit that permits the direct exchange of solutes between the cytoplasm of adjacent cells and thus mediate cell-cell ion and metabolic signaling. The lack of high resolution data for connexon structures has hampered so far the study of the structure-function relationships that link molecular effects of disease-causing mutations with their observed phenotypes. Here we present a combination of modeling techniques and molecular dynamics (MD) to infer side chain positions starting from low resolution structures containing only C alpha atoms. We validated this procedure on the structure of the KcsA potassium channel, which is solved at atomic resolution. We then produced a fully atomistic model of a homotypic Cx32 connexon starting from a published model of the C alpha carbons arrangement for the connexin transmembrane helices, to which we added extracellular and cytoplasmic loops. To achieve structural relaxation within a realistic environment, we used MD simulations inserted in an explicit solvent-membrane context and we subsequently checked predictions of putative side chain positions and interactions in the Cx32 connexon against a vast body of experimental reports. Our results provide new mechanistic insights into the effects of numerous spontaneous mutations and their implication in connexin-related pathologies. This model constitutes a step forward towards a structurally detailed description of the gap junction architecture and provides a structural platform to plan new biochemical and biophysical experiments aimed at elucidating the structure of connexin channels and hemichannels.
| S-EPMC2481295 | biostudies-literature