scholarly journals Millisecond-scale molecular dynamics simulation of spike RBD structure reveals evolutionary adaption of SARS-CoV-2 to stably bind ACE2

2020 ◽  
Author(s):  
Gard Nelson ◽  
Oleksandr Buzko ◽  
Aaron Bassett ◽  
Patricia R Spilman ◽  
Kayvan Niazi ◽  
...  

The Receptor Binding Domain (RBD) of the SARS-CoV-2 surface spike (S) protein interacts with host angiotensin converting enzyme 2 (ACE2) to gain entry to host cells and initiate infection 1-3. Detailed, accurate understanding of key interactions between S RBD and ACE2 provides critical information that may be leveraged in the development of strategies for the prevention and treatment of COVID-19. Utilizing the published sequences and cryo-EM structures of both the viral S RBD and ACE2 4,5, we performed in silico molecular dynamics (MD) simulations of free S RBD and of its interaction with ACE2 over the exceptionally long durations of 2.9 and 2 milliseconds, respectively, to elucidate the nature and relative affinity of S RBD surface residues for the ACE2 binding region. Our findings reveal that free S RBD has assumed an optimized ACE2 binding-ready conformation, incurring little entropic penalty for binding, an evolutionary adaptation that contributes to its high affinity for the receptor 6. We further identified high probability molecular binding interactions that inform both vaccine design and therapeutic development, which may include recombinant ACE2-based spike decoys 7 and/or allosteric S RBD-ACE2 binding inhibitors 8,9 to prevent or arrest infection and thus disease.

Processes ◽  
2021 ◽  
Vol 9 (10) ◽  
pp. 1749
Author(s):  
Yuqi Zhang ◽  
Li Chen ◽  
Xiaoyu Wang ◽  
Yanyan Zhu ◽  
Yongsheng Liu ◽  
...  

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) is a type of Ribonucleic Acid (RNA) coronavirus and it has infected and killed many people around the world. It is reported that the receptor binding domain of the spike protein (S_RBD) of the SARS-CoV-2 virus is responsible for attachment to human angiotensin converting enzyme II (ACE2). Many researchers are attempting to search potential inhibitors for fighting SARS-CoV-2 infection using theoretical or experimental methods. In terms of experimental and theoretical research, Cefuroxime, Erythromycin, Lincomycin and Ofloxacin are the potential inhibitors of SARS-CoV-2. However, the interactive mechanism of the protein SARS-CoV-2 and the inhibitors are still elusive. Here, we investigated the interactions between S_RBD and the inhibitors using molecular dynamics (MD) simulations. Interestingly, we found that there are two binding sites of S_RBD for the four small molecules. In addition, our analysis also illustrated that hydrophobic and π-π stacking interactions play crucial roles in the interactions between S_RBD and the small molecules. In our work, we also found that small molecules with glycosyl group have more effect on the conformation of S_RBD than other inhibitors, and they are also potential inhibitors for the genetic variants of SARS-CoV-2. This study provides in silico-derived mechanistic insights into the interactions of S_RBD and inhibitors, which may provide new clues for fighting SARS-CoV-2 infection.


2020 ◽  
Vol 23 (1) ◽  
pp. 41-54 ◽  
Author(s):  
Nouman Rasool ◽  
Waqar Hussain

Background: Human Immunodeficiency Virus 1 (HIV-1) is a lentivirus, which causes various HIV-associated infections. The HIV-1 core dissociation is essential for viral cDNA synthesis and phosphorylation of HIV-1 capsid protein (HIV-1 CA) plays an important role in it. Objective: The aim of this study was to explicate the role of three phosphoserine sites i.e. Ser109, Ser149 and Ser178 in the structural stability of HIV-1 CA, and it’s binding with GS-CA1, a novel potent inhibitor. Method: Eight complexes were analyzed and Molecular Dynamics (MD) simulations were performed to observe the stability of HIV-1 CA in the presence and absence of phosphorylation of serine residues at four different temperatures i.e. 300K, 325K, 340K and 350K, along with molecular docking and DFT analysis. Results: The structures showed maximum stability in the presence of phosphorylated serine residue. However, GS-CA1 docked most strongly with the native structure of HIV-1 CA i.e. binding affinity was -8.5 kcal/mol (Ki = 0.579 µM). Conclusion: These results suggest that the phosphorylation of these three serine residues weakens the binding of GS-CA1 with CA and casts derogatory effect on inhibition potential of this inhibitor, but it supports the stability of HIV-1 CA structure that can enhance regulation and replication of HIV-1 in host cells.


2021 ◽  
Author(s):  
Mert Golcuk ◽  
Aysima Hacisuleyman ◽  
Sema Zeynep Yilmaz ◽  
Elhan Taka ◽  
Ahmet Yildiz ◽  
...  

The Delta variant spreads more rapidly than previous variants of SARS-CoV-2. This variant comprises several mutations on the receptor-binding domain (RBD_Delta) of its spike (S) glycoprotein, which binds to the peptidase domain (PD) of angiotensin-converting enzyme 2 (ACE2) receptors in host cells. The RBD-PD interaction has been targeted by antibodies and nanobodies to prevent viral infection, but their effectiveness against the Delta variant remains unclear. Here, we investigated RBD_Delta-PD interactions in the presence and absence of nanobodies H11-H4, H11-D4, and Ty1 by performing in a total of 19 µs all-atom molecular dynamics (MD) simulations. Unbiased simulations revealed that Delta variant mutations strengthen RBD binding to ACE2 by increasing the hydrophobic interactions and salt bridge formation, but weaken interactions with H11-H4, H11-D4, and Ty1. Consequently, these nanobodies are unable to dislocate ACE2 from RBD_Delta. Steered MD simulations at comparable loading rates to atomic force microscopy (AFM) experiments estimated lower rupture forces of the nanobodies from RBD_Delta compared to ACE2. Our results suggest that existing nanobodies are less effective to inhibit RBD_Delta-PD interactions and a new generation of nanobodies will be needed to neutralize the Delta variant.


Sensors ◽  
2021 ◽  
Vol 21 (8) ◽  
pp. 2621
Author(s):  
Seunghwa Yang

Here, we systematically interrogate the effects of grafting single-walled (SWNT) and multi-walled carbon nanotubes (MWNT) to polymer matrices by using molecular dynamics (MD) simulations. We specifically investigate key material properties that include interfacial load transfer, alteration of nanotube properties, and dispersion of nanotubes in the polymer matrix. Simulations are conducted on a periodic unit cell model of the nanocomposite with a straight carbon nanotube and an amorphous polyethylene terephthalate (PET) matrix. For each type of nanotube, either 0%, 1.55%, or 3.1% of the carbon atoms in the outermost nanotubes are covalently grafted onto the carbon atoms of the PET matrix. Stress-strain curves and the elastic moduli of nanotubes and nanocomposites are determined based on the density of covalent grafting. Covalent grafting promotes two rivalling effects with respect to altering nanotube properties, and improvements in interfacial load transfer in the nanocomposites are clearly observed. The enhanced interface enables external loads applied to the nanocomposites to be efficiently transferred to the grafted nanotubes. Covalent functionalization of the nanotube surface with PET molecules can alter the solubility of nanotubes and improve dispersibility. Finally, we discuss the current limitations and challenges in using molecular modelling strategies to accurately predict properties on the nanotube and polymers systems studied here.


2008 ◽  
Vol 32 ◽  
pp. 255-258
Author(s):  
Bohayra Mortazavi ◽  
Akbar Afaghi Khatibi

Molecular Dynamics (MD) are now having orthodox means for simulation of matter in nano-scale. It can be regarded as an accurate alternative for experimental work in nano-science. In this paper, Molecular Dynamics simulation of uniaxial tension of some face centered cubic (FCC) metals (namely Au, Ag, Cu and Ni) at nano-level have been carried out. Sutton-Chen potential functions and velocity Verlet formulation of Noise-Hoover dynamic as well as periodic boundary conditions were applied. MD simulations at different loading rates and temperatures were conducted, and it was concluded that by increasing the temperature, maximum engineering stress decreases while engineering strain at failure is increasing. On the other hand, by increasing the loading rate both maximum engineering stress and strain at failure are increasing.


2014 ◽  
Vol 1700 ◽  
pp. 61-66
Author(s):  
Guttormur Arnar Ingvason ◽  
Virginie Rollin

ABSTRACTAdding single walled carbon nanotubes (SWCNT) to a polymer matrix can improve the delamination properties of the composite. Due to the complexity of polymer molecules and the curing process, few 3-D Molecular Dynamics (MD) simulations of a polymer-SWCNT composite have been run. Our model runs on the Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS), with a COMPASS (Condensed phase Optimized Molecular Potential for Atomistic Simulations Studies) potential. This potential includes non-bonded interactions, as well as bonds, angles and dihedrals to create a MD model for a SWCNT and EPON 862/DETDA (Diethyltoluenediamine) polymer matrix. Two simulations were performed in order to test the implementation of the COMPASS parameters. The first one was a tensile test on a SWCNT, leading to a Young’s modulus of 1.4 TPa at 300K. The second one was a pull-out test of a SWCNT from an originally uncured EPON 862/DETDA matrix.


2021 ◽  
Vol 12 ◽  
Author(s):  
Trina Ekawati Tallei ◽  
Fatimawali ◽  
Afriza Yelnetty ◽  
Rinaldi Idroes ◽  
Diah Kusumawaty ◽  
...  

The rapid spread of a novel coronavirus known as SARS-CoV-2 has compelled the entire world to seek ways to weaken this virus, prevent its spread and also eliminate it. However, no drug has been approved to treat COVID-19. Furthermore, the receptor-binding domain (RBD) on this viral spike protein, as well as several other important parts of this virus, have recently undergone mutations, resulting in new virus variants. While no treatment is currently available, a naturally derived molecule with known antiviral properties could be used as a potential treatment. Bromelain is an enzyme found in the fruit and stem of pineapples. This substance has been shown to have a broad antiviral activity. In this article, we analyse the ability of bromelain to counteract various variants of the SARS-CoV-2 by targeting bromelain binding on the side of this viral interaction with human angiotensin-converting enzyme 2 (hACE2) using molecular docking and molecular dynamics simulation approaches. We have succeeded in making three-dimensional configurations of various RBD variants using protein modelling. Bromelain exhibited good binding affinity toward various variants of RBDs and binds right at the binding site between RBDs and hACE2. This result is also presented in the modelling between Bromelain, RBD, and hACE2. The molecular dynamics (MD) simulations study revealed significant stability of the bromelain and RBD proteins separately up to 100 ns with an RMSD value of 2 Å. Furthermore, despite increases in RMSD and changes in Rog values of complexes, which are likely due to some destabilized interactions between bromelain and RBD proteins, two proteins in each complex remained bonded, and the site where the two proteins bind remained unchanged. This finding indicated that bromelain could have an inhibitory effect on different SARS-CoV-2 variants, paving the way for a new SARS-CoV-2 inhibitor drug. However, more in vitro and in vivo research on this potential mechanism of action is required.


1998 ◽  
Vol 540 ◽  
Author(s):  
J.M. Perlado ◽  
L. Malerba ◽  
T. Diaz De La Rubia

AbstractMolecular Dynamics (MD) simulations of neutron damage in β-SiC have been performed using a modified version of the Tersoff potential. The Threshold Displacement Energy (TDE) for Si and C atoms at 300 K has been determined along directions [001], [110], [111] and [ 1 1 1 ]. The existence of recombination barriers, which allow the formation of metastable, temperature-sensitive defects even below the threshold, has been observed. Displacement cascades produced by both C- and Si-recoils of energies spanning from 0.5 keV up to, respectively, 5 keV and 8 keV have also been simulated at 300 K and 1300 K. Their analysis, together with the analysis of damage accumulation (∼3.4×10-3 DPA) at 1300 K, reveals that the two sub-lattices exhibit opposite responses to irradiation: whereas only a little damage is produced on the “ductile” Si sub-lattice, many point-defects accumulate on the much more “fragile” C sub-lattice. A preliminary study of the nature and clustering tendency of these defects is performed. The possibility of disorder-induced amorphization is considered and the preliminary result is that no amorphization takes place at the dose and temperature simulated.


2015 ◽  
Vol 817 ◽  
pp. 797-802 ◽  
Author(s):  
Cai Jiang ◽  
Jian Wei Zhang ◽  
Shao Feng Lin ◽  
Su Ju ◽  
Da Zhi Jiang

Molecular dynamics (MD) simulations on three single walled carbon nanotube (SWCNT) reinforced epoxy resin composites were conducted to study the influence of SWCNT type on the glass transition temperature (Tg) of the composites. The composite matrix is cross-linked epoxy resin based on the epoxy monomers bisphenol A diglycidyl ether (DGEBA) cured by diaminodiphenylmethane (DDM). MD simulations of NPT (constant number of particles, constant pressure and constant temperature) dynamics were carried out to obtain density as a function of temperature for each composite system. The Tg was determined as the temperature corresponding to the discontinuity of plot slopes of the densityvsthe temperature. In order to understand the motion of polymer chain segments above and below the Tg, various energy components and the MSD at various temperatures of the composites were investigated and their roles played in the glass transition process were analyzed. The results show that the Tg of the composites increases with increasing aspect ratio of the embedded SWCNT


2014 ◽  
pp. 466-477
Author(s):  
Nabila Tahreen ◽  
K. M. Masud

In recent years, polymer/carbon nanotube composites have attracted increased attention because the polymer properties have significantly improved. In this paper, a single walled carbon nanotube (SWCNT) is used to reinforce polystyrene matrix. Molecular dynamics (MD) simulations are used to study two periodic systems - a long CNT-reinforced polystyrene composite and amorphous polystyrene matrix itself. The axial and transverse elastic moduli of the amorphous polystyrene matrix and nanocomposites are evaluated using constant-strain energy minimization method. The results from MD simulations are compared with corresponding rule-of-mixture predictions. The simulation results show that CNTs significantly improve the stiffness of polystyrene/CNT composite, especially in the longitudinal direction of the nanotube. Polystyrene posses a strong attractive interaction with the surface of the SWCNT and therefore play an important role in providing effective adhesion. The conventional rule-of-mixture predicts a smaller value than MD simulation where there are strong interfacial interactions. Here the authors report a study on the interfacial characteristics of a CNT-PS composite system through MD simulations and continuum mechanics.


Sign in / Sign up

Export Citation Format

Share Document