The application of the MM/GBSA method in the binding pose prediction of FGFR inhibitors

2020 ◽  
Vol 22 (17) ◽  
pp. 9656-9663 ◽  
Author(s):  
Yu Chen ◽  
Yongxiang Zheng ◽  
Pedro Fong ◽  
Shengjun Mao ◽  
Qiantao Wang
Keyword(s):  
Free Energy ◽  
Binding Free Energy ◽  
Md Simulations ◽  
Pose Prediction

The correct conformation had lower MM/GBSA binding free energy in longer MD simulations for each FGFR1 inhibitor.

Estimation of the ligand-binding free energy of checkpoint kinase 1 via non-equilibrium MD simulations

2020 ◽  
Vol 100 ◽  
pp. 107648 ◽  
Author(s):  
Nguyen Thi Mai ◽  
Ngo Thi Lan ◽  
Thien Y Vu ◽  
Phuong Thi Mai Duong ◽  
Nguyen Thanh Tung ◽  
...  
Keyword(s):  
Free Energy ◽  
Ligand Binding ◽  
Binding Free Energy ◽  
Md Simulations ◽  
Checkpoint Kinase ◽  
Checkpoint Kinase 1 ◽  
Non Equilibrium

scholarly journals Molecular mechanism and binding free energy of doxorubicin intercalation in DNA

2019 ◽  
Vol 21 (7) ◽  
pp. 3877-3893 ◽  
Author(s):  
Bahaa Jawad ◽  
Lokendra Poudel ◽  
Rudolf Podgornik ◽  
Nicole F. Steinmetz ◽  
Wai-Yim Ching
Keyword(s):  
Free Energy ◽  
Molecular Mechanism ◽  
Binding Free Energy ◽  
Md Simulations ◽  
Intercalation Process

The intercalation process of binding doxorubicin (DOX) in DNA is studied by extensive MD simulations.

scholarly journals Identification of promising drug candidates against Non-Structural Protein 15 (NSP15) from SARS-CoV-2: an in silico assisted drug-repurposing studies

2020 ◽  
Author(s):  
Rameez Jabeer Khan ◽  
Rajat Kumar Jha ◽  
Ekampreet Singh ◽  
Monika Jain ◽  
Gizachew Muluneh Amera ◽  
...  
Keyword(s):  
Free Energy ◽  
In Silico ◽  
Binding Free Energy ◽  
Drug Repurposing ◽  
Md Simulations ◽  
Therapeutic Agents ◽  
Global Dynamics ◽  
Host Immune System ◽  
Structural Protein ◽  
Drug Candidates

<div>The recent COVID-19 pandemic caused by SARS-CoV-2 has recorded a high number of infected people across the globe. The notorious nature of the virus makes it necessary for us to identify promising therapeutic agents in a time-sensitive manner. The current study utilises an <i>in silico</i> based drug repurposing approach to identify potential drug candidates targeting non-structural protein 15 (NSP15), i.e. a uridylate specific endoribonuclease of SARS-CoV-2</div><div>which plays an indispensable role in RNA processing and viral immune evasion from the host immune system. NSP15 was screened against an in-house library of 123 antiviral drugs obtained from the DrugBank database from which three promising drug candidates were identified based on their estimated free energy of binding (<i>ΔG</i>), estimated inhibition constant (<i>Ki</i>), the orientation of drug molecules in the active site and the key interacting residues of</div><div>NSP15. The MD simulations were performed for the selected NSP15-drug complexes along with free protein to mimic on their physiological state. The binding free energies of the selected NSP15-drug complexes were also calculated using the trajectories of MD simulations of NSP15-drug complexes through MM/PBSA (Molecular Mechanics with Poisson-Boltzmann and surface area solvation) approach where NSP15-Simeprevir (-242.559 kJ/mol) and NSP15-Paritaprevir (-149.557 kJ/mol) exhibited the strongest binding affinities. Together with the results of molecular docking, global dynamics, essential dynamics and binding free energy analysis, we propose that Simeprevir and Paritaprevir are promising drug candidates for the inhibition of NSP15 and could act as potential therapeutic agents against SARS-CoV-2.</div>

scholarly journals Possibility of HIV-1 Protease Inhibitors-Clinical Trial Drugs as Repurposed Drugs for SARSCoV-2 Main Protease: A Molecular Docking, Molecular Dynamics and Binding Free Energy Simulation Study

2020 ◽  
Author(s):  
Ancy Iruthayaraj ◽  
Sivanandam Magudeeswaran ◽  
Kumaradhas Poomani
Keyword(s):  
Clinical Trial ◽  
Amino Acids ◽  
Free Energy ◽  
Binding Free Energy ◽  
Drug Repurposing ◽  
Md Simulations ◽  
Binding Mechanism ◽  
Main Protease ◽  
Repurposed Drugs ◽  
Hiv 1

<p>Initially, the SARS-CoV-2 virus was emerged from Wuhan, China and rapidly spreading across the world and urges the scientific community to develop antiviral therapeutic agents. Among several strategies, drug repurposing will help to react immediately to overcome COVID-19 pandemic. In the present study, we have chosen two clinical trial drugs TMB607 and TMC310911 are the inhibitors of HIV-1 protease to use as the inhibitors of SARS-CoV-2 main protease (M<sup>pro</sup>) enzyme. To make use of these two inhibitors as the repurposed drugs for COVID-19, it is essential to know the molecular basis of binding mechanism of these two molecules with the SARS-CoV-2 main protease (M<sup>pro</sup>). Understand the binding mechanism; we performed the molecular docking, molecular dynamics (MD) simulations and binding free energy calculations against the SARS-CoV-2 M<sup>pro</sup>. The docking results indicate that both molecules form intermolecular interactions with the active site amino acids of M<sup>pro</sup> enzyme. However, during the MD simulations, TMB607 forms strong interactions with the key amino acids of M<sup>pro</sup> and remains intact. The RMSD and RMSF values of both complexes were stable throughout the MD simulations. The MM-GBSA binding free energy values of both complexes are -43.7 and -34.9 kcal/mol, respectively. This <i>in silico</i> study proves that the TMB607 molecule binds strongly with the SARS-CoV-2 M<sup>pro</sup> enzyme and it is suitable for the drug repurposing of COVID-19 and further drug designing.</p>

Insights to the Binding of a Selective Adenosine A3 Receptor Antagonist using Molecular Dynamics Simulations, Binding Free Energy Calculations and Mutagenesis

2019 ◽  
Author(s):  
Panagiotis Lagarias ◽  
Kerry Barkan ◽  
Eva Tzortzini ◽  
Eleni Vrontaki ◽  
Margarita Stampelou ◽  
...  
Keyword(s):  
Free Energy ◽  
Computational Model ◽  
Molecular Mechanics ◽  
Surface Area ◽  
Binding Free Energy ◽  
Md Simulations ◽  
Free Energy Calculations ◽  
Binding Energies ◽  
Energy Calculations ◽  
Binding Free Energy Calculations

<p>Adenosine A<sub>3 </sub>receptor (A<sub>3</sub>R), is a promising drug target against cancer cell proliferation. Currently there is no experimentally determined structure of A<sub>3</sub>R. Here, we have investigate a computational model, previously applied successfully for agonists binding to A<sub>3</sub>R, using molecular dynamic (MD) simulations, Molecular Mechanics-Poisson Boltzmann Surface Area (MM-PBSA) and Molecular Mechanics-Generalized Born Surface Area (MM-GBSA) binding free energy calculations. Extensive computations were performed to explore the binding profile of O4-{[3-(2,6-dichlorophenyl)-5-methylisoxazol-4-yl]carbonyl}-2-methyl-1,3-thiazole-4-carbohydroximamide (K18) to A<sub>3</sub>R. K18 is a new specific and competitive antagonist at the orthosteric binding site of A<sub>3</sub>R, discovered using virtual screening and characterized pharmacologically in our previous studies. The most plausible binding conformation for the dichlorophenyl group of K18 inside the A<sub>3</sub>R is oriented towards trans-membrane helices (TM) 5 and 6, according to the MM-PBSA and MM-GBSA binding free energy calculations, and by the previous results obtained by mutating residues of TM5, TM6 to alanine which reduce antagonist potency. The results from 14 site-directed mutagenesis experiments were interpreted using MD simulations and MM-GBSA calculations which show that the relative binding free energies of the mutant A<sub>3</sub>R - K18 complexes compare to the WT A<sub>3</sub>R are in agreement with the effect of the mutations, i.e. the reduction, maintenance or increase of antagonist potency. We show that when the residues V169<sup>5.30</sup>, M177<sup>5.38</sup>, I249<sup>6.54</sup> involved in direct interactions with K18 are mutated to alanine, the mutant A<sub>3</sub>R - K18 complexes reduce potency, increase the RMSD value of K18 inside the binding area and the MM-GBSA binding free energy compared to the WT A<sub>3</sub>R complex. Our computational model shows that other mutant A<sub>3</sub>R complexes with K18, including directly interacting residues, i.e. F168<sup>5.29</sup>A, L246<sup>6.51</sup>A, N250<sup>6.55</sup>A complexes with K18 are not stable. In these complexes of A<sub>3</sub>R mutated in directly interacting residues one or more of the interactions between K18 and these residues are lost. In agreement with the experiments, the computations show that, M174<sup>5.35</sup> a residue which does not make direct interactions with K18 is critical for K18 binding. A striking results is that the mutation of residue V169<sup>5.30</sup> to glutamic acid maintained antagonistic potency. This effect is in agreement with the binding free energy calculations and it is suggested that is due to K18 re-orientation but also to the plasticity of A<sub>3</sub>R binding area. The mutation of direct interacting L90<sup>3.32</sup> in the low region and the non-directly interacting L264<sup>7.35</sup> to alanine in the middle region increases the antagonistic potency, suggesting that chemical modifications of K18 can be applied to augment antagonistic potency. The calculated binding energies Δ<i>G</i><sub>eff</sub> values of K18 against mutant A<sub>3</sub>Rs displayed very good correlation with experimental potencies (pA<sub>2</sub> values). These results further approve the computational model for the description of K18 binding with critical residues of the orthosteric binding area which can have implications for the design of more effective antagonists based on the structure of K18.</p>

scholarly journals Accounting for the Central Role of Interfacial Water in Protein-Ligand Binding Free Energy Calculations

2020 ◽  
Author(s):  
Ido Ben-Shalom ◽  
Zhixiong Lin ◽  
Brian Radak ◽  
Charles Lin ◽  
Woody Sherman ◽  
...  
Keyword(s):  
Free Energy ◽  
Binding Site ◽  
Binding Free Energy ◽  
Md Simulations ◽  
Simulation Software ◽  
Free Energy Calculations ◽  
Computer Hardware ◽  
Water Molecules ◽  
Energy Calculations ◽  
Binding Free Energy Calculations

Rigorous binding free energy methods in drug discovery are growing in popularity due to a combination of methodological advances, improvements in computer hardware, and workflow automation. These calculations typically use molecular dynamics (MD) to sample from the Boltzmann distribution of conformational states. However, when part or all the binding site is inaccessible to bulk solvent, the time needed for water molecules to equilibrate between bulk solvent and the binding site can be well beyond what is practical with standard MD. This sampling limitation is problematic in relative binding free energy calculations, which compute the reversible work of converting Ligand 1 to Ligand 2 within the binding site. Thus, if Ligand 1 is smaller and/or more polar than Ligand 2, the perturbation may allow additional water molecules to occupy a region of the binding site. However, this change in hydration may not be captured by standard MD simulations and may therefore lead to errors in the computed free energy. We recently developed a hybrid Monte Carlo/MD (MC/MD) method, which speeds the equilibration of water between bulk solvent and buried cavities, while sampling from the intended distribution of states. Here, we report on the use of this approach in the context of alchemical binding free energy calculations. We find that using MC/MD markedly improves the accuracy of the calculations and also reduces hysteresis between the forward and reverse perturbations, relative to matched calculations using only MD with or without the crystallographic water molecules. The present method is available for use in the AMBER simulation software.<br>

scholarly journals Corrigendum to ‘Estimation of the ligand-binding free energy of checkpoint kinase 1 via non-equilibrium MD simulations’

2021 ◽  
Vol 103 ◽  
pp. 107805
Author(s):  
Nguyen Thi Mai ◽  
Ngo Thi Lan ◽  
Thien Y Vu ◽  
Phuong Thi Mai Duong ◽  
Nguyen Thanh Tung ◽  
...  
Keyword(s):  
Free Energy ◽  
Ligand Binding ◽  
Binding Free Energy ◽  
Md Simulations ◽  
Checkpoint Kinase ◽  
Checkpoint Kinase 1 ◽  
Non Equilibrium

scholarly journals Protein-Protein Interactions: Insight from Molecular Dynamics Simulations and Nanoparticle Tracking Analysis

Molecules ◽  
2021 ◽  
Vol 26 (18) ◽  
pp. 5696
Author(s):  
Wei Lim Chong ◽  
Koollawat Chupradit ◽  
Sek Peng Chin ◽  
Mai Mai Khoo ◽  
Sook Mei Khor ◽  
...  
Keyword(s):  
Molecular Dynamics ◽  
Free Energy ◽  
Protein Binding ◽  
Protein Interactions ◽  
Binding Free Energy ◽  
Single Point ◽  
Md Simulations ◽  
Nanoparticle Tracking Analysis ◽  
Buffer Solution ◽  
Dynamics Simulations

Protein-protein interaction plays an essential role in almost all cellular processes and biological functions. Coupling molecular dynamics (MD) simulations and nanoparticle tracking analysis (NTA) assay offered a simple, rapid, and direct approach in monitoring the protein-protein binding process and predicting the binding affinity. Our case study of designed ankyrin repeats proteins (DARPins)—AnkGAG1D4 and the single point mutated AnkGAG1D4-Y56A for HIV-1 capsid protein (CA) were investigated. As reported, AnkGAG1D4 bound with CA for inhibitory activity; however, it lost its inhibitory strength when tyrosine at residue 56 AnkGAG1D4, the most key residue was replaced by alanine (AnkGAG1D4-Y56A). Through NTA, the binding of DARPins and CA was measured by monitoring the increment of the hydrodynamic radius of the AnkGAG1D4-gold conjugated nanoparticles (AnkGAG1D4-GNP) and AnkGAG1D4-Y56A-GNP upon interaction with CA in buffer solution. The size of the AnkGAG1D4-GNP increased when it interacted with CA but not AnkGAG1D4-Y56A-GNP. In addition, a much higher binding free energy (∆GB) of AnkGAG1D4-Y56A (−31 kcal/mol) obtained from MD further suggested affinity for CA completely reduced compared to AnkGAG1D4 (−60 kcal/mol). The possible mechanism of the protein-protein binding was explored in detail by decomposing the binding free energy for crucial residues identification and hydrogen bond analysis.

scholarly journals Exploring the Binding Interaction of Raf Kinase Inhibitory Protein With the N-Terminal of C-Raf Through Molecular Docking and Molecular Dynamics Simulation

2021 ◽  
Vol 8 ◽  
Author(s):  
Shraddha Parate ◽  
Shailima Rampogu ◽  
Gihwan Lee ◽  
Jong Chan Hong ◽  
Keun Woo Lee
Keyword(s):  
Molecular Dynamics ◽  
Free Energy ◽  
Binding Sites ◽  
Binding Free Energy ◽  
Md Simulations ◽  
Dynamics Simulation ◽  
Web Servers ◽  
Inhibitory Protein ◽  
Raf Kinase ◽  
Raf Kinase Inhibitory Protein

Protein-protein interactions are indispensable physiological processes regulating several biological functions. Despite the availability of structural information on protein-protein complexes, deciphering their complex topology remains an outstanding challenge. Raf kinase inhibitory protein (RKIP) has gained substantial attention as a favorable molecular target for numerous pathologies including cancer and Alzheimer’s disease. RKIP interferes with the RAF/MEK/ERK signaling cascade by endogenously binding with C-Raf (Raf-1 kinase) and preventing its activation. In the current investigation, the binding of RKIP with C-Raf was explored by knowledge-based protein-protein docking web-servers including HADDOCK and ZDOCK and a consensus binding mode of C-Raf/RKIP structural complex was obtained. Molecular dynamics (MD) simulations were further performed in an explicit solvent to sample the conformations for when RKIP binds to C-Raf. Some of the conserved interface residues were mutated to alanine, phenylalanine and leucine and the impact of mutations was estimated by additional MD simulations and MM/PBSA analysis for the wild-type (WT) and constructed mutant complexes. Substantial decrease in binding free energy was observed for the mutant complexes as compared to the binding free energy of WT C-Raf/RKIP structural complex. Furthermore, a considerable increase in average backbone root mean square deviation and fluctuation was perceived for the mutant complexes. Moreover, per-residue energy contribution analysis of the equilibrated simulation trajectory by HawkDock and ANCHOR web-servers was conducted to characterize the key residues for the complex formation. One residue each from C-Raf (Arg398) and RKIP (Lys80) were identified as the druggable “hot spots” constituting the core of the binding interface and corroborated by additional long-time scale (300 ns) MD simulation of Arg398Ala mutant complex. A notable conformational change in Arg398Ala mutant occurred near the mutation site as compared to the equilibrated C-Raf/RKIP native state conformation and an essential hydrogen bonding interaction was lost. The thirteen binding sites assimilated from the overall analysis were mapped onto the complex as surface and divided into active and allosteric binding sites, depending on their location at the interface. The acquired information on the predicted 3D structural complex and the detected sites aid as promising targets in designing novel inhibitors to block the C-Raf/RKIP interaction.

Sign in / Sign up

Export Citation Format

Share Document