scholarly journals Computational Insights Into the Effects of the R190K and N121Q Mutations on the SARS-CoV-2 Spike Complex With Biliverdin

2021 ◽  
Vol 8 ◽  
Author(s):  
Zhiyuan Qu ◽  
Kaihang Li ◽  
Xiaoju Geng ◽  
Bo Huang ◽  
Jian Gao
Keyword(s):  
Conformational Changes ◽  
Binding Mode ◽  
Dynamics Simulation ◽  
Side Chain ◽  
Energy Calculation ◽  
Upward Movement ◽  
Wild Type ◽  
Generalized Born ◽  
Main Target ◽  
Antibody Design

The SARS-CoV-2 spike has been regarded as the main target of antibody design against COVID-19. Two single-site mutations, R190K and N121Q, were deemed to weaken the binding affinity of biliverdin although the underlying molecular mechanism is still unknown. Meanwhile, the effect of the two mutations on the conformational changes of “lip” and “gate” loops was also elusive. Thus, molecular dynamics simulation and molecular mechanics/generalized Born surface area (MM/GBSA) free energy calculation were conducted on the wild-type and two other SARS-CoV-2 spike mutants. Our simulations indicated that the R190K mutation causes Lys190 to form six hydrogen bonds, guided by Asn99 and Ile101, which brings Lys190 closer to Arg102 and Asn121, thereby weakening the interaction energy between biliverdin and Ile101 as well as Lys190. For the N121Q mutation, Gln121 still maintained a hydrogen bond with biliverdin; nevertheless, the overall binding mode deviated significantly under the reversal of the side chain of Phe175. Moreover, the two mutants would stabilize the lip loop, which would restrain the meaningful upward movement of the lip. In addition, N121Q significantly promoted the gate loop deviating to the biliverdin binding site and compressed the site. This work would be useful in understanding the dynamics binding biliverdin to the SARS-CoV-2 spike.

scholarly journals Based on the Virtual Screening of Multiple Pharmacophores, Docking and Molecular Dynamics Simulation Approaches toward the Discovery of Novel HPPD Inhibitors

2020 ◽  
Vol 21 (15) ◽  
pp. 5546
Author(s):  
Ying Fu ◽  
Tong Ye ◽  
Yong-Xuan Liu ◽  
Jian Wang ◽  
Fei Ye
Keyword(s):  
Coulomb Force ◽  
Docking Study ◽  
Binding Mode ◽  
Homogentisic Acid ◽  
Dynamics Simulation ◽  
Catabolic Pathway ◽  
Molecular Docking Study ◽  
Generalized Born ◽  
Pharmacophore Models ◽  
Important Enzyme

4-Hydroxyphenylpyruvate dioxygenase (HPPD) is an iron-dependent non-heme oxygenase involved in the catabolic pathway of tyrosine, which is an important enzyme in the transformation of 4-hydroxyphenylpyruvic acid to homogentisic acid, and thus being considered as herbicide target. Within this study, a set of multiple structure-based pharmacophore models for HPPD inhibitors were developed. The ZINC and natural product database were virtually screened, and 29 compounds were obtained. The binding mode of HPPD and its inhibitors obtained through molecular docking study showed that the residues of Phe424, Phe381, His308, His226, Gln307 and Glu394 were crucial for activity. Molecular-mechanics-generalized born surface area (MM/GBSA) results showed that the coulomb force, lipophilic and van der Waals (vdW) interactions made major contributions to the binding affinity. These efforts will greatly contribute to design novel and effective HPPD inhibitory herbicides.

scholarly journals Molecular Dynamics of SARS-CoV-2 Nsp1 Structural Changes Associated with Variant Mutations

2021 ◽  
Author(s):  
Shokouh Rezaei ◽  
Yahya Sefidbakht ◽  
Filipe Pereira
Keyword(s):  
Molecular Dynamics ◽  
Conformational Changes ◽  
Structural Changes ◽  
Structure And Function ◽  
Dynamics Simulation ◽  
Structural Protein ◽  
Wild Type ◽  
Deletion Mutations ◽  
And Function

Abstract SARS-CoV-2 non-structural protein 1 (Nsp1) is a virulence factor that inhibits the translation of host mRNAs and interact with viral RNA. Despite the relevance of Nsp1, few studies have been conducted to understand the effect of mutations on Nsp1 structure and function. Here, we provide a molecular dynamics simulation of SARS-CoV-2 Nsp1, wild type and variants. We found that SARS-CoV-2 Nsp1 has a more Rg value than SARS-CoV-1 Nsp1, with indicate an effect on the folding protein. This result suggest that SARS-CoV-2 Nsp1 can more easily approach the active site of the ribosome compared to SARS-CoV-1 Nsp1. In addition, we found that the C-terminal of the SARS-CoV-2 Nsp1, in particular residues 164 to 170, are more flexible than other regions of SARS-CoV-2 Nsp1 and SARS-CoV-1 Nsp1, confirming the role of this region in the interaction with the 40S subunit. Moreover, multiple deletion mutations have been found in the N/C-terminal of the SARS-CoV-2 Nsp1, which seems the effect of SARS-CoV-2 Nsp1 multiple deletions is greater than that of substitutions. Among all deletions, D156-158 and D80-90 may destabilize the protein structure and possibly increase the virulence of the SARS-CoV-2. Overall, our findings reinforce the importance of studying Nsp1 conformational changes in new variants and its effect on virulence of SARS-CoV-2.

scholarly journals The Impact of Mutation L138F/L210F on the Orai Channel: A Molecular Dynamics Simulation Study

2021 ◽  
Vol 8 ◽  
Author(s):  
Xiaoqian Zhang ◽  
Hua Yu ◽  
Xiangdong Liu ◽  
Chen Song
Keyword(s):  
Molecular Dynamics ◽  
Conformational Changes ◽  
Calcium Release ◽  
Hydrophobic Interactions ◽  
Dynamics Simulation ◽  
Wild Type ◽  
Hydrophobic Region ◽  
In The Wild ◽  
Dynamics Simulations ◽  
The Impact

The calcium release-activated calcium channel, composed of the Orai channel and the STIM protein, plays a crucial role in maintaining the Ca2+ concentration in cells. Previous studies showed that the L138F mutation in the human Orai1 creates a constitutively open channel independent of STIM, causing severe myopathy, but how the L138F mutation activates Orai1 is still unclear. Here, based on the crystal structure of Drosophila melanogaster Orai (dOrai), molecular dynamics simulations for the wild-type (WT) and the L210F (corresponding to L138F in the human Orai1) mutant were conducted to investigate their structural and dynamical properties. The results showed that the L210F dOrai mutant tends to have a more hydrated hydrophobic region (V174 to F171), as well as more dilated basic region (K163 to R155) and selectivity filter (E178). Sodium ions were located deeper in the mutant than in the wild-type. Further analysis revealed two local but essential conformational changes that may be the key to the activation. A rotation of F210, a previously unobserved feature, was found to result in the opening of the K163 gate through hydrophobic interactions. At the same time, a counter-clockwise rotation of F171 occurred more frequently in the mutant, resulting in a wider hydrophobic gate with more hydration. Ultimately, the opening of the two gates may facilitate the opening of the Orai channel independent of STIM.

scholarly journals Discovery of the Novel Inhibitor Against New Delhi Metallo-β-Lactamase Based on Virtual Screening and Molecular Modelling

2020 ◽  
Vol 21 (10) ◽  
pp. 3567 ◽  
Author(s):  
Xiyan Wang ◽  
Yanan Yang ◽  
Yawen Gao ◽  
Xiaodi Niu
Keyword(s):  
Virtual Screening ◽  
Bacterial Resistance ◽  
Competitive Inhibition ◽  
Binding Free Energy ◽  
Binding Mode ◽  
Hydrolytic Activity ◽  
Dynamics Simulation ◽  
Energy Calculation ◽  
New Delhi ◽  
The Novel

New Delhi metallo-β-lactamase (NDM-1), one of the metallo-β-lactamases (MBLs), leads to antibiotic resistance in clinical treatments due to the strong ability of hydrolysis to almost all kinds of β-lactam antibiotics. Therefore, there is the urgent need for the research and development of the novel drug-resistant inhibitors targeting NDM-1. In this study, ZINC05683641 was screened as potential NDM-1 inhibitor by virtual screening and the inhibitor mechanism of this compound was explored based on molecular dynamics simulation. The nitrocefin assay showed that the IC50 value of ZINC05683641 was 13.59 ± 0.52 μM, indicating that the hydrolytic activity of NDM-1 can be obviously suppressed by ZINC05683641. Further, the binding mode of ZINC05683641 with NDM-1 was obtained by molecular modeling, binding free energy calculation, mutagenesis assays and fluorescence-quenching assays. As results, ILE-35, MET-67, VAL-73, TRP-93, CYS-208, ASN-220 and HIS-250 played the key roles in the binding of NDM-1 with ZINC05683641. Interestingly, these key residues were exactly located in the catalytic activity region of NDM-1, implying that the inhibitor mechanism of ZINC05683641 against NDM-1 was the competitive inhibition. These findings will provide an available approach to research and develop new drug against NDM-1 and treatment for bacterial resistance.

scholarly journals Structural Basis of Local, pH-Dependent Conformational Changes in Glycoprotein B from Herpes Simplex Virus Type 1

2010 ◽  
Vol 84 (24) ◽  
pp. 12924-12933 ◽  
Author(s):  
Samuel D. Stampfer ◽  
Huan Lou ◽  
Gary H. Cohen ◽  
Roselyn J. Eisenberg ◽  
Ekaterina E. Heldwein
Keyword(s):  
Herpes Simplex ◽  
Conformational Changes ◽  
Fusion Proteins ◽  
Low Ph ◽  
Side Chain ◽  
Wild Type ◽  
Ph Dependent ◽  
Simplex Virus ◽  
Mutant Forms

ABSTRACT Herpesviruses enter cells by membrane fusion either at the plasma membrane or in endosomes, depending on the cell type. Glycoprotein B (gB) is a conserved component of the multiprotein herpesvirus fusion machinery and functions as a fusion protein, with two internal fusion loops, FL1 and FL2. We determined the crystal structures of the ectodomains of two FL1 mutants of herpes simplex virus type 1 (HSV-1) gB to clarify whether their fusion-null phenotypes were due to global or local effects of the mutations on the structure of the gB ectodomain. Each mutant has a single point mutation of a hydrophobic residue in FL1 that eliminates the hydrophobic side chain. We found that neither mutation affected the conformation of FL1, although one mutation slightly altered the conformation of FL2, and we conclude that the fusion-null phenotype is due to the absence of a hydrophobic side chain at the mutated position. Because the ectodomains of the wild-type and the mutant forms of gB crystallized at both low and neutral pH, we were able to determine the effect of pH on gB conformation at the atomic level. For viruses that enter cells by endocytosis, the low pH of the endosome effects major conformational changes in their fusion proteins, thereby promoting fusion of the viral envelope with the endosomal membrane. We show here that upon exposure of gB to low pH, FL2 undergoes a major relocation, probably driven by protonation of a key histidine residue. Relocation of FL2, as well as additional small conformational changes in the gB ectodomain, helps explain previously noted changes in its antigenic and biochemical properties. However, no global pH-dependent changes in gB structure were detected in either the wild-type or the mutant forms of gB. Thus, low pH causes local conformational changes in gB that are very different from the large-scale fusogenic conformational changes in other viral fusion proteins. We propose that these conformational changes, albeit modest, play an important functional role during endocytic entry of HSV.

Binding Modes of Ligands Using Enhanced Sampling (BLUES): Rapid Decorrelation of Ligand Binding Modes Using Nonequilibrium Candidate Monte Carlo

2017 ◽  
Author(s):  
Samuel Gill ◽  
Nathan M. Lim ◽  
Patrick Grinaway ◽  
Ariën S. Rustenburg ◽  
Josh Fass ◽  
...  
Keyword(s):  
Monte Carlo ◽  
Ligand Binding ◽  
Binding Mode ◽  
Free Energy Calculations ◽  
Dynamics Simulation ◽  
Binding Modes ◽  
Enhanced Sampling ◽  
Recent Success ◽  
Multiple Binding ◽  
Proper Sampling

<div>Accurately predicting protein-ligand binding is a major goal in computational chemistry, but even the prediction of ligand binding modes in proteins poses major challenges. Here, we focus on solving the binding mode prediction problem for rigid fragments. That is, we focus on computing the dominant placement, conformation, and orientations of a relatively rigid, fragment-like ligand in a receptor, and the populations of the multiple binding modes which may be relevant. This problem is important in its own right, but is even more timely given the recent success of alchemical free energy calculations. Alchemical calculations are increasingly used to predict binding free energies of ligands to receptors. However, the accuracy of these calculations is dependent on proper sampling of the relevant ligand binding modes. Unfortunately, ligand binding modes may often be uncertain, hard to predict, and/or slow to interconvert on simulation timescales, so proper sampling with current techniques can require prohibitively long simulations. We need new methods which dramatically improve sampling of ligand binding modes. Here, we develop and apply a nonequilibrium candidate Monte Carlo (NCMC) method to improve sampling of ligand binding modes.</div><div><br></div><div>In this technique the ligand is rotated and subsequently allowed to relax in its new position through alchemical perturbation before accepting or rejecting the rotation and relaxation as a nonequilibrium Monte Carlo move. When applied to a T4 lysozyme model binding system, this NCMC method shows over two orders of magnitude improvement in binding mode sampling efficiency compared to a brute force molecular dynamics simulation. This is a first step towards applying this methodology to pharmaceutically relevant binding of fragments and, eventually, drug-like molecules. We are making this approach available via our new Binding Modes of Ligands using Enhanced Sampling (BLUES) package which is freely available on GitHub.</div>

Binding Modes of Ligands Using Enhanced Sampling (BLUES): Rapid Decorrelation of Ligand Binding Modes Using Nonequilibrium Candidate Monte Carlo

2018 ◽  
Author(s):  
Samuel Gill ◽  
Nathan M. Lim ◽  
Patrick Grinaway ◽  
Ariën S. Rustenburg ◽  
Josh Fass ◽  
...  
Keyword(s):  
Monte Carlo ◽  
Ligand Binding ◽  
Binding Mode ◽  
Free Energy Calculations ◽  
Dynamics Simulation ◽  
Binding Modes ◽  
Enhanced Sampling ◽  
Recent Success ◽  
Multiple Binding ◽  
Proper Sampling

<div>Accurately predicting protein-ligand binding is a major goal in computational chemistry, but even the prediction of ligand binding modes in proteins poses major challenges. Here, we focus on solving the binding mode prediction problem for rigid fragments. That is, we focus on computing the dominant placement, conformation, and orientations of a relatively rigid, fragment-like ligand in a receptor, and the populations of the multiple binding modes which may be relevant. This problem is important in its own right, but is even more timely given the recent success of alchemical free energy calculations. Alchemical calculations are increasingly used to predict binding free energies of ligands to receptors. However, the accuracy of these calculations is dependent on proper sampling of the relevant ligand binding modes. Unfortunately, ligand binding modes may often be uncertain, hard to predict, and/or slow to interconvert on simulation timescales, so proper sampling with current techniques can require prohibitively long simulations. We need new methods which dramatically improve sampling of ligand binding modes. Here, we develop and apply a nonequilibrium candidate Monte Carlo (NCMC) method to improve sampling of ligand binding modes.</div><div><br></div><div>In this technique the ligand is rotated and subsequently allowed to relax in its new position through alchemical perturbation before accepting or rejecting the rotation and relaxation as a nonequilibrium Monte Carlo move. When applied to a T4 lysozyme model binding system, this NCMC method shows over two orders of magnitude improvement in binding mode sampling efficiency compared to a brute force molecular dynamics simulation. This is a first step towards applying this methodology to pharmaceutically relevant binding of fragments and, eventually, drug-like molecules. We are making this approach available via our new Binding Modes of Ligands using Enhanced Sampling (BLUES) package which is freely available on GitHub.</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