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

<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>

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

10.26434/chemrxiv.10032218.v1 ◽

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 ◽

<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>

Can molecular dynamics simulations improve the structural accuracy and virtual screening performance of GPCR models?

PLoS Computational Biology ◽

10.1371/journal.pcbi.1008936 ◽

2021 ◽

Vol 17 (5) ◽

pp. e1008936

Author(s):

Jon Kapla ◽

Ismael Rodriguez Espigares ◽

Flavio Ballante ◽

Jana Selent ◽

Jens Carlsson

Keyword(s):

Crystal Structure ◽

Molecular Dynamics ◽

Virtual Screening ◽

Ligand Binding ◽

Md Simulations ◽

Binding Mode ◽

Extracellular Loop ◽

Loop Region ◽

Screening Performance ◽

Virtual Screening Performance

The determination of G protein-coupled receptor (GPCR) structures at atomic resolution has improved understanding of cellular signaling and will accelerate the development of new drug candidates. However, experimental structures still remain unavailable for a majority of the GPCR family. GPCR structures and their interactions with ligands can also be modelled computationally, but such predictions have limited accuracy. In this work, we explored if molecular dynamics (MD) simulations could be used to refine the accuracy of in silico models of receptor-ligand complexes that were submitted to a community-wide assessment of GPCR structure prediction (GPCR Dock). Two simulation protocols were used to refine 30 models of the D3 dopamine receptor (D3R) in complex with an antagonist. Close to 60 μs of simulation time was generated and the resulting MD refined models were compared to a D3R crystal structure. In the MD simulations, the transmembrane helix region of the models generally drifted further away from the crystal structure conformation. However, MD refinement was able to improve the accuracy of the ligand binding mode and the second extracellular loop region. The best refinement protocol improved agreement with the experimentally observed ligand binding mode for a majority of the models. Receptor structures with improved virtual screening performance, which was assessed by molecular docking of ligands and decoys, could also be identified among the MD refined models. Application of weak restraints to the transmembrane helixes in the MD simulations further improved predictions of the ligand binding mode and second extracellular loop. These results provide guidelines for application of MD refinement in prediction of GPCR-ligand complexes and directions for further method development.

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

10.26434/chemrxiv.5406907.v2 ◽

2017 ◽

Author(s):

Samuel Gill ◽

Nathan M. Lim ◽

Patrick Grinaway ◽

Ariën S. Rustenburg ◽

Josh Fass ◽

...

Keyword(s):

Monte Carlo ◽

Ligand Binding ◽

Binding Mode ◽

Dynamics Simulation ◽

Binding Modes ◽

Enhanced Sampling ◽

Recent Success ◽

Multiple Binding ◽

<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

10.26434/chemrxiv.5406907 ◽

2018 ◽

Author(s):

Samuel Gill ◽

Nathan M. Lim ◽

Patrick Grinaway ◽

Ariën S. Rustenburg ◽

Josh Fass ◽

...

Keyword(s):

Monte Carlo ◽

Ligand Binding ◽

Binding Mode ◽

Dynamics Simulation ◽

Binding Modes ◽

Enhanced Sampling ◽

Recent Success ◽

Multiple Binding ◽

<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

Bioinformatics ◽

10.1093/bioinformatics/btaa565 ◽

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.

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

10.26434/chemrxiv.5406907.v5 ◽

2018 ◽

Author(s):

Samuel Gill ◽

Nathan M. Lim ◽

Patrick Grinaway ◽

Ariën S. Rustenburg ◽

Josh Fass ◽

...

Keyword(s):

Monte Carlo ◽

Ligand Binding ◽

Binding Mode ◽

Dynamics Simulation ◽

Binding Modes ◽

Enhanced Sampling ◽

Recent Success ◽

Multiple Binding ◽

<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

10.26434/chemrxiv.5406907.v3 ◽

2017 ◽

Author(s):

Samuel Gill ◽

Nathan M. Lim ◽

Patrick Grinaway ◽

Ariën S. Rustenburg ◽

Josh Fass ◽

...

Keyword(s):

Monte Carlo ◽

Ligand Binding ◽

Binding Mode ◽

Dynamics Simulation ◽

Binding Modes ◽

Enhanced Sampling ◽

Recent Success ◽

Multiple Binding ◽

<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>

Enhancing Sidechain Rotamer Sampling Using Non-Equilibrium Candidate Monte Carlo

10.26434/chemrxiv.7184876 ◽

2019 ◽

Author(s):

Kalistyn H. Burley ◽

Samuel C. Gill ◽

Nathan M. Lim ◽

David Mobley

Keyword(s):

Molecular Dynamics ◽

Monte Carlo ◽

Ligand Binding ◽

Electrostatic Interactions ◽

Molecular Simulations ◽

Adequate Sampling ◽

Binding Pockets ◽

Sidechain Rotamers ◽

Dynamics Simulations ◽

<div>Molecular simulations are a valuable tool for studying biomolecular motions and thermodynamics. However, such motions can be slow compared to simulation timescales, yet critical. Specifically, adequate sampling of sidechain motions in protein binding pockets proves crucial for obtaining accurate estimates of ligand binding free energies from molecular simulations. The timescale of sidechain rotamer flips can range from a few ps to several hundred ns or longer, particularly in crowded environments like the interior of proteins. Here, we apply a mixed non-equilibrium candidate Monte Carlo (NCMC)/molecular dynamics (MD) method to enhance sampling of sidechain rotamers. The NCMC portion of our method applies a switching protocol wherein the steric and electrostatic interactions between target sidechain atoms and the surrounding environment are cycled off and then back on during the course of a move proposal. Between NCMC move proposals, simulation of the system continues via traditional molecular dynamics. Here, we first validate this approach on a simple, solvated valine-alanine dipeptide system and then apply it to a well-studied model ligand binding site in T4 lysozyme L99A. We compute the rate of rotamer transitions for a valine sidechain using our approach and compare it to that of traditional molecular dynamics simulations. Here, we show that our NCMC/MD method substantially enhances sidechain sampling, especially in systems where the torsional barrier to rotation is high (>10 kcal/mol). These barriers can be intrinsic torsional barriers or steric barriers imposed by the environment.</div><div>Overall, this may provide a promising strategy to selectively improve sidechain sampling in molecular simulations.</div>

Enhancing Sidechain Rotamer Sampling Using Non-Equilibrium Candidate Monte Carlo

10.26434/chemrxiv.7184876.v1 ◽

2018 ◽

Author(s):

Kalistyn H. Burley ◽

Samuel C. Gill ◽

Nathan M. Lim ◽

David Mobley

Keyword(s):

Molecular Dynamics ◽

Monte Carlo ◽

Ligand Binding ◽

Electrostatic Interactions ◽

Molecular Simulations ◽

Adequate Sampling ◽

Binding Pockets ◽

Sidechain Rotamers ◽

Dynamics Simulations ◽

<div>Molecular simulations are a valuable tool for studying biomolecular motions and thermodynamics. However, such motions can be slow compared to simulation timescales, yet critical. Specifically, adequate sampling of sidechain motions in protein binding pockets proves crucial for obtaining accurate estimates of ligand binding free energies from molecular simulations. The timescale of sidechain rotamer flips can range from a few ps to several hundred ns or longer, particularly in crowded environments like the interior of proteins. Here, we apply a mixed non-equilibrium candidate Monte Carlo (NCMC)/molecular dynamics (MD) method to enhance sampling of sidechain rotamers. The NCMC portion of our method applies a switching protocol wherein the steric and electrostatic interactions between target sidechain atoms and the surrounding environment are cycled off and then back on during the course of a move proposal. Between NCMC move proposals, simulation of the system continues via traditional molecular dynamics. Here, we first validate this approach on a simple, solvated valine-alanine dipeptide system and then apply it to a well-studied model ligand binding site in T4 lysozyme L99A. We compute the rate of rotamer transitions for a valine sidechain using our approach and compare it to that of traditional molecular dynamics simulations. Here, we show that our NCMC/MD method substantially enhances sidechain sampling, especially in systems where the torsional barrier to rotation is high (>10 kcal/mol). These barriers can be intrinsic torsional barriers or steric barriers imposed by the environment.</div><div>Overall, this may provide a promising strategy to selectively improve sidechain sampling in molecular simulations.</div>