scholarly journals Enhancing Sampling of Water Rearrangements on Ligand Binding: A Comparison of Techniques

2021 ◽  
Author(s):  
Yunhui Ge ◽  
David C. Wych ◽  
Marley L. Samways ◽  
Michael E. Wall ◽  
Jonathan W. Essex ◽  
...  
Keyword(s):  
Monte Carlo ◽  
Binding Site ◽  
Md Simulation ◽  
Direct Analysis ◽  
Md Simulations ◽  
Binding Mode ◽  
Bulk Water ◽  
Free Energy Calculations ◽  
Water Sampling ◽  
Ligand Interactions

Water often plays a key role in protein structure, molecular recognition, and mediating protein-ligand interactions. Thus, free energy calculations must adequately sample water motions, which often proves challenging in typical MD simulation timescales. Thus, the accuracy of methods relying on MD simulations ends up limited by slow water sampling. Particularly, as a ligand is removed or modified, bulk water may not have time to fill or rearrange in the binding site. In this work, we focus on several molecular dynamics (MD) simulation-based methods attempting to help address water motions and occupancies: BLUES, using nonequilibrium candidate Monte Carlo (NCMC); grand, using grand canonical Monte Carlo (GCMC); and normal MD. We assess the accuracy and efficiency of these methods in sampling water motions. We selected a range of systems with varying numbers of waters in the binding site, as well as those where water occupancy is coupled to the identity or binding mode of the ligand. We analyzed water motions and occupancies using both clustering of trajectories and direct analysis of electron density maps. Our results suggest both BLUES and grand enhance water sampling relative to normal MD and grand is more robust than BLUES, but also that water sampling remains a major challenge for all of the methods tested. The lessons we learned for these methods and systems are discussed.

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>

A boosted unbiased molecular dynamics method for predicting ligands binding mechanisms: probing the binding pathway of dasatinib to Src-kinase

2020 ◽  
Vol 36 (18) ◽  
pp. 4714-4720
Author(s):  
Farzin Sohraby ◽  
Mostafa Javaheri Moghadam ◽  
Masoud Aliyar ◽  
Hassan Aryapour
Keyword(s):  
Molecular Dynamics ◽  
Ligand Binding ◽  
Md Simulation ◽  
Src Kinase ◽  
Deep Understanding ◽  
Md Simulations ◽  
Binding Mode ◽  
Supplementary Information ◽  
Computational Resources ◽  
Binding Mechanisms

Abstract Summary Small molecules such as metabolites and drugs play essential roles in biological processes and pharmaceutical industry. Knowing their interactions with biomacromolecular targets demands a deep understanding of binding mechanisms. Dozens of papers have suggested that discovering of the binding event by means of conventional unbiased molecular dynamics (MD) simulation urges considerable amount of computational resources, therefore, only one who holds a cluster or a supercomputer can afford such extensive simulations. Thus, many researchers who do not own such resources are reluctant to take the benefits of running unbiased MD simulation, in full atomistic details, when studying a ligand binding pathway. Many researchers are impelled to be content with biased MD simulations which seek its validation due to its intrinsic preconceived framework. In this work, we have presented a workable stratagem to encourage everyone to perform unbiased (unguided) MD simulations, in this case a protein–ligand binding process, by typical desktop computers and so achieve valuable results in nanosecond time scale. Here, we have described a dynamical binding’s process of an anticancer drug, the dasatinib, to the c-Src kinase in full atomistic details for the first time, without applying any biasing force or potential which may lead the drug to artificial interactions with the protein. We have attained multiple independent binding events which occurred in the nanosecond time scales, surprisingly as little as ∼30 ns. Both the protonated and deprotonated forms of the dasatinib reached the crystallographic binding mode without having any major intermediate state during induction. Availability and implementation The links of the tutorial and technical documents are accessible in the article. Supplementary information Supplementary data are available at Bioinformatics online.

scholarly journals Controlling the Heterodimerisation of the Phytosulfokine Receptor 1 (PSKR1) via Island Loop Modulation

2020 ◽  
Author(s):  
Joao Victor de Souza Cunha ◽  
Matthew Kondal ◽  
Piotr Zaborniak ◽  
Ryland Cairns ◽  
Agnieszka K. Bronowska
Keyword(s):  
Small Molecules ◽  
Binding Site ◽  
Cell Communication ◽  
Md Simulations ◽  
Free Energy Calculations ◽  
Agricultural Applications ◽  
Lrr Domain ◽  
Signalling Cascade ◽  
New Crop ◽  
Resolution Structure

Phytosulfokine (PSK) is a phytohormone responsible for cell-to-cell communication in plants, playing pivotal role in plant development and growth. The binding of PSK to its cognate receptor, PSKR1, is modulated by the formation of a binding site located between leucine-rich repeat (LRR) domain of PSKR1 and the loop located in the receptor’s island domain (ID). The atomic resolution structure of the extracellular PSKR1 bound to PSK has been reported, however, the intrinsic dynamics of PSK binding and the architecture of PSKR1 binding site remain to be understood. In this work, we used atomistic molecular dynamics (MD) simulations and free energy calculations to elucidate how the PSKR1 island domain (ID) loop forms and binds PSK. Moreover, we report a novel “druggable” binding site which could be exploited for the targeted modulation of the PSKR1-PSK binding by small molecules. We expect that our results will open new ways to modulate the PSK signalling cascade via small molecules, which can result in new crop control and agricultural applications.

scholarly journals Enhancing Water Sampling of Buried Binding Sites Using Nonequilibrium Candidate Monte Carlo

2020 ◽  
Author(s):  
Teresa Danielle Bergazin ◽  
Ido Ben-Shalom ◽  
Nathan M. Lim ◽  
Samuel C. Gill ◽  
Michael K. Gilson ◽  
...  
Keyword(s):  
Molecular Dynamics ◽  
Monte Carlo ◽  
Binding Sites ◽  
Water Exchange ◽  
Free Energy Calculations ◽  
Water Molecules ◽  
Water Sampling ◽  
Computer Aided Drug Design ◽  
Practical Applications ◽  
Computer Aided

<div>Water molecules can be found interacting with the surface and within cavities in proteins. However, water exchange between bulk and buried hydration sites can be slow compared to simulation timescales, thus leading to the inefficient sampling of the locations of water. This can pose problems for free energy calculations for computer-aided drug design. Here, we apply a hybrid method that combines nonequilibrium candidate Monte Carlo (NCMC) simulations and molecular dynamics (MD) to enhance sampling of water in specific areas of a system, such as the binding site of a protein. Our approach uses NCMC to gradually remove interactions between a selected water molecule and its environment, then translates the water to a new region, before turning the interactions back on. This approach of gradual removal of interactions, followed by a move and then reintroduction of interactions, allows the environment relax in response to the proposed water translation, improving acceptance of moves and thereby accelerating water exchange and sampling. We validate this approach on several test systems including the ligand-bound MUP-1 and HSP90 proteins with buried crystallographic waters removed. We show that our NCMC/MD method enhances water sampling relative to normal MD when applied to these systems. Thus, this approach provides a strategy to improve water sampling in molecular simulations which may be useful in practical applications in drug discovery and biomolecular design.</div>

scholarly journals Controlling the Heterodimerisation of the Phytosulfokine Receptor 1 (PSKR1) via Island Loop Modulation

2020 ◽  
Author(s):  
Joao Victor de Souza Cunha ◽  
Matthew Kondal ◽  
Piotr Zaborniak ◽  
Ryland Cairns ◽  
Agnieszka K. Bronowska
Keyword(s):  
Small Molecules ◽  
Binding Site ◽  
Cell Communication ◽  
Md Simulations ◽  
Free Energy Calculations ◽  
Agricultural Applications ◽  
Lrr Domain ◽  
Signalling Cascade ◽  
New Crop ◽  
Resolution Structure

Phytosulfokine (PSK) is a phytohormone responsible for cell-to-cell communication in plants, playing pivotal role in plant development and growth. The binding of PSK to its cognate receptor, PSKR1, is modulated by the formation of a binding site located between leucine-rich repeat (LRR) domain of PSKR1 and the loop located in the receptor’s island domain (ID). The atomic resolution structure of the extracellular PSKR1 bound to PSK has been reported, however, the intrinsic dynamics of PSK binding and the architecture of PSKR1 binding site remain to be understood. In this work, we used atomistic molecular dynamics (MD) simulations and free energy calculations to elucidate how the PSKR1 island domain (ID) loop forms and binds PSK. Moreover, we report a novel “druggable” binding site which could be exploited for the targeted modulation of the PSKR1-PSK binding by small molecules. We expect that our results will open new ways to modulate the PSK signalling cascade via small molecules, which can result in new crop control and agricultural applications.

Fragment Pose Prediction Using Non-Equilibrium Candidate Monte Carlo and Molecular Dynamics Simulations

2019 ◽  
Author(s):  
Nathan M. Lim ◽  
Meghan Osato ◽  
Gregory L. Warren ◽  
David L. Mobley
Keyword(s):  
Crystal Structure ◽  
Molecular Dynamics ◽  
Monte Carlo ◽  
Ligand Binding ◽  
Md Simulations ◽  
Binding Mode ◽  
Protein Crystal ◽  
Pose Prediction ◽  
Binding Modes ◽  
Non Equilibrium

<div>Part of early stage drug discovery involves determining how molecules may bind to the target protein. Through understanding where and how molecules bind, chemists can begin to build ideas on how to design improvements to increase binding affinities. In this retrospective study, we compare how computational approaches like docking, molecular dynamics (MD) simulations, and a non-equilibrium candidate Monte Carlo (NCMC) based method (NCMC+MD) perform in predicting binding modes for a set of 12 fragment-like molecules which bind to soluble epoxide hydrolase. We evaluate each method's effectiveness in identifying the dominant binding mode and finding any additional binding modes (if any). Then, we compare our predicted binding modes to experimentally obtained X-ray crystal structures.</div><div>We dock each of the 12 small molecules into the apo-protein crystal structure and then run simulations up to 1 microsecond each. Small and fragment-like molecules likely have smaller energy barriers separating different binding modes by virtue of relatively fewer and weaker interactions relative to drug-like molecules, and thus likely undergo more rapid binding mode transitions. We expect, thus, to see more rapid transitions betweeen binding modes in our study. </div><div><br></div><div>Following this, we build Markov State Models (MSM) to define our stable ligand binding modes. We investigate if adequate sampling of ligand binding modes and transitions between them can occur at the microsecond timescale using traditional MD or a hybrid NCMC+MD simulation approach. Our findings suggest that even with small fragment-like molecules, we fail to sample all the crystallographic binding modes using microsecond MD simulations, but using NCMC+MD we have better success in sampling the crystal structure while obtaining the correct populations.</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>

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>

Targeting Peptidyl-prolyl Cis-trans Isomerase NIMA-interacting 1: A Structure-based Virtual Screening Approach to Find Novel Inhibitors

2020 ◽  
Vol 16 (5) ◽  
pp. 605-617 ◽  
Author(s):  
Kauê Santana da Costa ◽  
João M. Galúcio ◽  
Deivid Almeida de Jesus ◽  
Guelber Cardoso Gomes ◽  
Anderson Henrique Lima e Lima ◽  
...  
Keyword(s):  
Inhibitory Activity ◽  
Binding Free Energy ◽  
Md Simulations ◽  
Binding Mode ◽  
Free Energy Calculations ◽  
Pentacyclic Triterpenoids ◽  
Substrate Peptide ◽  
Molecular Mechanism Of Action ◽  
Structural Insights ◽  
Virtual Screening Approach

Background : Peptidyl-prolyl cis-trans isomerase NIMA-interacting 1 (Pin1) is an enzyme that isomerizes phosphorylated serine or threonine motifs adjacent to proline residues. Pin1 has important roles in several cellular signaling pathways, consequently impacting the development of multiple types of cancers. Methods: Based on the previously reported inhibitory activity of pentacyclic triterpenoids isolated from the gum resin of Boswellia genus against Pin1, we designed a computational experiment using molecular docking, pharmacophore filtering, and structural clustering allied to molecular dynamics (MD) simulations and binding free energy calculations to explore the inhibitory activity of new triterpenoids against Pin1 structure. Results: Here, we report different computational evidence that triterpenoids from neem (Azadirachta indica A. Juss), such as 6-deacetylnimbinene, 6-Oacetylnimbandiol, and nimbolide, replicate the binding mode of the Pin1 substrate peptide, interacting with high affinity with the binding site and thus destabilizing the Pin1 structure. Conclusion: Our results are supported by experimental data, and provide interesting structural insights into their molecular mechanism of action, indicating that their structural scaffolds could be used as a start point to develop new inhibitors against Pin1.

Sign in / Sign up

Export Citation Format

Share Document