scholarly journals Molecular Dynamics Studies on the Structural Characteristics for the Stability Prediction of SARS-CoV-2

2021 ◽  
Vol 22 (16) ◽  
pp. 8714
Author(s):  
Kwang-Eun Choi ◽  
Jeong-Min Kim ◽  
JeeEun Rhee ◽  
Ae Kyung Park ◽  
Eun-Jin Kim ◽  
...  
Keyword(s):  
Binding Free Energy ◽  
Structural Characteristics ◽  
3D Structure ◽  
Direct Interaction ◽  
Spike Protein ◽  
Energy Calculation ◽  
Protein Chain ◽  
Complex Forms ◽  
The Stability ◽  
Almost All

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) affects the COVID-19 pandemic in the world. The spike protein of the various proteins encoded in SARS-CoV-2 binds to human ACE2, fuses, and enters human cells in the respiratory system. Spike protein, however, is highly variable, and many variants were identified continuously. In this study, Korean mutants for spike protein (D614G and D614A-C terminal domain, L455F and F456L-RBD, and Q787H-S2 domain) were investigated in patients. Because RBD in spike protein is related to direct interaction with ACE2, almost all researches were focused on the RBD region or ACE2-free whole domain region. The 3D structure for spike protein complexed with ACE2 was recently released. The stability analysis through RBD distance among each spike protein chain and the binding free energy calculation between spike protein and ACE2 were performed using MD simulation depending on mutant types in 1-, 2-, and 3-open-complex forms. D614G mutant of CT2 domain, showing to be the most prevalent in the global pandemic, showed higher stability in all open-complex forms than the wild type and other mutants. We hope this study will provide an insight into the importance of conformational fluctuation in the whole domain, although RBD is involved in the direct interaction with ACE2.

scholarly journals Hydrogen Bond Stability of Quinazoline Derivatives Compounds in Complex against EGFR using Molecular Dynamics Simulation

2019 ◽  
Vol 19 (2) ◽  
pp. 461
Author(s):  
Herlina Rasyid ◽  
Bambang Purwono ◽  
Thomas S Hofer ◽  
Harno Dwi Pranowo
Keyword(s):  
Molecular Dynamics ◽  
Hydrogen Bond ◽  
Free Energy ◽  
Molecular Dynamics Simulation ◽  
Binding Free Energy ◽  
Dynamics Simulation ◽  
Inhibition Mechanism ◽  
Energy Calculation ◽  
Protein Receptor ◽  
The Stability

Lung cancer was a second common cancer case due to the high cigarette smoking activity both in men and women. One of protein receptor which plays an important role in the growth of the tumor is Epidermal Growth Factor Receptor (EGFR). EGFR protein is the most frequent protein mutation in cancer and promising target to inhibit the cancer growth. In this work, the stability of the hydrogen bond as the main interaction in the inhibition mechanism of cancer will be evaluated using molecular dynamics simulation. There were two compounds (A1 and A2) as new potential inhibitors that were complexed against the EGFR protein. The dynamic properties of each complexed were compared with respect to erlotinib against EGFR. The result revealed that both compounds had an interaction in the main catalytic area of protein receptor which is at methionine residue. Inhibitor A1 showed additional interactions during simulation time but the interactions tend to be weak. Inhibitor A2 displayed a more stable interaction. Following dynamics simulation, binding free energy calculation was performed by two scoring techniques MM/GB(PB)SA method and gave a good correlation with the stability of the complex. Furthermore, potential inhibitor A2 had a lower binding free energy as a direct consequence of the stability of hydrogen bond interaction.

scholarly journals An Effective MM/GBSA Protocol for Absolute Binding Free Energy Calculations: A Case Study on SARS-CoV-2 Spike Protein and the Human ACE2 Receptor

Molecules ◽  
2021 ◽  
Vol 26 (8) ◽  
pp. 2383
Author(s):  
Negin Forouzesh ◽  
Nikita Mishra
Keyword(s):  
Free Energy ◽  
Binding Free Energy ◽  
Normal Mode Analysis ◽  
Free Energy Calculations ◽  
Spike Protein ◽  
Upper And Lower Bounds ◽  
Mode Analysis ◽  
Implicit Solvent ◽  
Energy Calculation ◽  
Generalized Born

The binding free energy calculation of protein–ligand complexes is necessary for research into virus–host interactions and the relevant applications in drug discovery. However, many current computational methods of such calculations are either inefficient or inaccurate in practice. Utilizing implicit solvent models in the molecular mechanics generalized Born surface area (MM/GBSA) framework allows for efficient calculations without significant loss of accuracy. Here, GBNSR6, a new flavor of the generalized Born model, is employed in the MM/GBSA framework for measuring the binding affinity between SARS-CoV-2 spike protein and the human ACE2 receptor. A computational protocol is developed based on the widely studied Ras–Raf complex, which has similar binding free energy to SARS-CoV-2/ACE2. Two options for representing the dielectric boundary of the complexes are evaluated: one based on the standard Bondi radii and the other based on a newly developed set of atomic radii (OPT1), optimized specifically for protein–ligand binding. Predictions based on the two radii sets provide upper and lower bounds on the experimental references: −14.7(ΔGbindBondi)<−10.6(ΔGbindExp.)<−4.1(ΔGbindOPT1) kcal/mol. The consensus estimates of the two bounds show quantitative agreement with the experiment values. This work also presents a novel truncation method and computational strategies for efficient entropy calculations with normal mode analysis. Interestingly, it is observed that a significant decrease in the number of snapshots does not affect the accuracy of entropy calculation, while it does lower computation time appreciably. The proposed MM/GBSA protocol can be used to study the binding mechanism of new variants of SARS-CoV-2, as well as other relevant structures.

scholarly journals De novo ssRNA Aptamers against the SARS-CoV-2 Main Protease: In Silico Design and Molecular Dynamics Simulation

2021 ◽  
Vol 22 (13) ◽  
pp. 6874
Author(s):  
Francesco Morena ◽  
Chiara Argentati ◽  
Ilaria Tortorella ◽  
Carla Emiliani ◽  
Sabata Martino
Keyword(s):  
De Novo ◽  
Binding Free Energy ◽  
3D Structure ◽  
Dynamics Simulation ◽  
Stable Complex ◽  
Therapeutic Drug ◽  
Energy Calculation ◽  
Potential Candidate ◽  
Initial Sequence ◽  
Main Protease

Herein, we have generated ssRNA aptamers to inhibit SARS-CoV-2 Mpro, a protease necessary for the SARS-CoV-2 coronavirus replication. Because there is no aptamer 3D structure currently available in the databanks for this protein, first, we modeled an ssRNA aptamer using an entropic fragment-based strategy. We refined the initial sequence and 3D structure by using two sequential approaches, consisting of an elitist genetic algorithm and an RNA inverse process. We identified three specific aptamers against SARS-CoV-2 Mpro, called MAptapro, MAptapro-IR1, and MAptapro-IR2, with similar 3D conformations and that fall in the dimerization region of the SARS-CoV-2 Mpro necessary for the enzymatic activity. Through the molecular dynamic simulation and binding free energy calculation, the interaction between the MAptapro-IR1 aptamer and the SARS-CoV-2 Mpro enzyme resulted in the strongest and the highest stable complex; therefore, the ssRNA MAptapro-IR1 aptamer was selected as the best potential candidate for the inhibition of SARS-CoV-2 Mpro and a perspective therapeutic drug for the COVID-19 disease.

Exploring the Stability of Inhibitor Binding to SIK2 Using Molecular Dynamics Simulation and Binding Free Energy Calculation

2021 ◽  
Author(s):  
Mingsong Shi ◽  
Min Zhao ◽  
Lun Wang ◽  
Kongjun Liu ◽  
Penghui Li ◽  
...  
Keyword(s):  
Binding Free Energy ◽  
Dynamics Simulation ◽  
Free Energy Calculation ◽  
Energy Calculation ◽  
Inhibitor Binding ◽  
Dependent Protein Kinase ◽  
Calcium Calmodulin Dependent ◽  
Biological Phenomena ◽  
Dependent Protein ◽  
The Stability

Salt inducible kinase 2 (SIK2) is a calcium/calmodulin-dependent protein kinase-like kinase that is implicated in a variety of biological phenomena, including cellular metabolism, growth, and apoptosis. SIK2 is the key...

scholarly journals Predicting the Potential Effect of E484K Mutation on the Binding of 28 Antibodies to the Spike Protein of SARS-CoV-2 by Molecular Dynamics Simulation and Free Energy Calculation

2021 ◽  
Author(s):  
Leyun Wu ◽  
Cheng Peng ◽  
Zhijian Xu ◽  
weiliang zhu
Keyword(s):  
Free Energy ◽  
Binding Affinity ◽  
Immune Escape ◽  
Binding Free Energy ◽  
Experimental Studies ◽  
Potential Effect ◽  
Dynamics Simulation ◽  
Spike Protein ◽  
Energy Calculation ◽  
Potential Impact

Vaccines and antibody therapeutic are needed to fight the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) that has spread since 2020. Experimental studies have shown that the E484K variant may escape the neutralization of antibodies. To explore the potential impact of E484K mutation on the antibody binding affinity, we calculated the binding free energy of 28 antibodies to the wild type and K484 mutant of the spike protein of SARS-CoV-2. We found that 71% of the antibodies show lower binding affinity to the E484K mutant, indicating the highly possible immune escape risk of the mutated virus. Further analysis revealed that the other mutations, e.g. F490 and V483, are also likely to cause immune escape.

scholarly journals Exploring the Association Between Sialic Acid and SARS-CoV-2 Spike Protein Through a Molecular Dynamics-Based Approach

2021 ◽  
Vol 2 ◽  
Author(s):  
Leonardo Bò ◽  
Mattia Miotto ◽  
Lorenzo Di Rienzo ◽  
Edoardo Milanetti ◽  
Giancarlo Ruocco
Keyword(s):  
Molecular Dynamics ◽  
Sialic Acid ◽  
Structural Characteristics ◽  
Principal Component ◽  
Dynamics Simulation ◽  
Spike Protein ◽  
Binding Effect ◽  
Terminal Domain ◽  
Correlated Motion ◽  
The Stability

Recent experimental evidence demonstrated the capability of SARS-CoV-2 Spike protein to bind sialic acid molecules, which was a trait not present in SARS-CoV and could shed light on the molecular mechanism used by the virus for the cell invasion. This peculiar feature has been successfully predicted by in-silico studies comparing the sequence and structural characteristics that SARS-CoV-2 shares with other sialic acid-binding viruses, like MERS-CoV. Even if the region of the binding has been identified in the N-terminal domain of Spike protein, so far no comprehensive analyses have been carried out on the spike-sialic acid conformations once in the complex. Here, we addressed this aspect performing an extensive molecular dynamics simulation of a system composed of the N-terminal domain of the spike protein and a sialic acid molecule. We observed several short-lived binding events, reconnecting to the avidic nature of the binding, interestingly occurring in the surface Spike region where several insertions are present with respect to the SARS-CoV sequence. Characterizing the bound configurations via a clustering analysis on the Principal Component of the motion, we identified different possible binding conformations and discussed their dynamic and structural properties. In particular, we analyze the correlated motion between the binding residues and the binding effect on the stability of atomic fluctuation, thus proposing regions with high binding propensity with sialic acid.

scholarly journals Predicting the Potential Effect of E484K Mutation on the Binding of 28 Antibodies to the Spike Protein of SARS-CoV-2 by Molecular Dynamics Simulation and Free Energy Calculation

2021 ◽  
Author(s):  
Leyun Wu ◽  
Cheng Peng ◽  
Zhijian Xu ◽  
weiliang zhu
Keyword(s):  
Free Energy ◽  
Binding Affinity ◽  
Immune Escape ◽  
Binding Free Energy ◽  
Experimental Studies ◽  
Potential Effect ◽  
Dynamics Simulation ◽  
Spike Protein ◽  
Energy Calculation ◽  
Potential Impact

Vaccines and antibody therapeutic are needed to fight the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) that has spread since 2020. Experimental studies have shown that the E484K variant may escape the neutralization of antibodies. To explore the potential impact of E484K mutation on the antibody binding affinity, we calculated the binding free energy of 28 antibodies to the wild type and K484 mutant of the spike protein of SARS-CoV-2. We found that 71% of the antibodies show lower binding affinity to the E484K mutant, indicating the highly possible immune escape risk of the mutated virus. Further analysis revealed that the other mutations, e.g. F490 and V483, are also likely to cause immune escape.

Large-Scale Assessment of Binding Free Energy Calculations in Active Drug Discovery Projects

2020 ◽  
Author(s):  
Christina Schindler ◽  
Hannah Baumann ◽  
Andreas Blum ◽  
Dietrich Böse ◽  
Hans-Peter Buchstaller ◽  
...  
Keyword(s):  
Free Energy ◽  
Drug Discovery ◽  
Large Scale ◽  
Active Drug ◽  
Binding Free Energy ◽  
Free Energy Calculations ◽  
Energy Calculation ◽  
Large Scale Assessment ◽  
New Public ◽  
Binding Free Energy Calculations

Here we present an evaluation of the binding affinity prediction accuracy of the free energy calculation method FEP+ on internal active drug discovery projects and on a large new public benchmark set.<br>

Thermodynamics of the interaction between the spike protein of severe acute respiratory syndrome- coronavirus-2 and the receptor of human angiotensin converting enzyme 2. Effects of possible ligands

2020 ◽  
Author(s):  
Cristina Garcia-Iriepa ◽  
Cecilia Hognon ◽  
Antonio Francés-Monerris ◽  
Isabel Iriepa ◽  
Tom Miclot ◽  
...  
Keyword(s):  
Angiotensin Converting Enzyme ◽  
Molecular Mechanisms ◽  
Molecular Design ◽  
Binding Free Energy ◽  
Transmembrane Protein ◽  
Spike Protein ◽  
Converting Enzyme ◽  
Angiotensin Converting Enzyme 2 ◽  
Rational Molecular Design ◽  
Rational Evaluation

<div><p>Since the end of 2019, the coronavirus SARS-CoV-2 has caused more than 180,000 deaths all over the world, still lacking a medical treatment despite the concerns of the whole scientific community. Human Angiotensin-Converting Enzyme 2 (ACE2) was recently recognized as the transmembrane protein serving as SARS-CoV-2 entry point into cells, thus constituting the first biomolecular event leading to COVID-19 disease. Here, by means of a state-of-the-art computational approach, we propose a rational evaluation of the molecular mechanisms behind the formation of the complex and of the effects of possible ligands. Moreover, binding free energy between ACE2 and the active Receptor Binding Domain (RBD) of the SARS-CoV-2 spike protein is evaluated quantitatively, assessing the molecular mechanisms at the basis of the recognition and the ligand-induced decreased affinity. These results boost the knowledge on the molecular grounds of the SARS-CoV-2 infection and allow to suggest rationales useful for the subsequent rational molecular design to treat severe COVID-19 cases.</p></div>

Application of the ESMACS Binding Free Energy Protocol to a Highly Varied Ligand Dataset: Lactate Dehydogenase A

2019 ◽  
Author(s):  
David Wright ◽  
Fouad Husseini ◽  
Shunzhou Wan ◽  
Christophe Meyer ◽  
Herman Van Vlijmen ◽  
...  
Keyword(s):  
Free Energy ◽  
Surface Area ◽  
Binding Free Energy ◽  
Normal Mode Analysis ◽  
Binding Mode ◽  
Accessible Surface Area ◽  
Solvent Accessible Surface Area ◽  
Mode Analysis ◽  
Energy Calculation ◽  
Accessible Surface

<div>Here, we evaluate the performance of our range of ensemble simulation based binding free energy calculation protocols, called ESMACS (enhanced sampling of molecular dynamics with approximation of continuum solvent) for use in fragment based drug design scenarios. ESMACS is designed to generate reproducible binding affinity predictions from the widely used molecular mechanics Poisson-Boltzmann surface area (MMPBSA) approach. We study ligands designed to target two binding pockets in the lactate dehydogenase A target protein, which vary in size, charge and binding mode. When comparing to experimental results, we obtain excellent statistical rankings across this highly diverse set of ligands. In addition, we investigate three approaches to account for entropic contributions not captured by standard MMPBSA calculations: (1) normal mode analysis, (2) weighted solvent accessible surface area (WSAS) and (3) variational entropy. </div>

Sign in / Sign up

Export Citation Format

Share Document