1
|
Ojha AA, Votapka LW, Amaro RE. QMrebind: incorporating quantum mechanical force field reparameterization at the ligand binding site for improved drug-target kinetics through milestoning simulations. Chem Sci 2023; 14:13159-13175. [PMID: 38023523 PMCID: PMC10664576 DOI: 10.1039/d3sc04195f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/11/2023] [Accepted: 10/22/2023] [Indexed: 12/01/2023] Open
Abstract
Understanding the interaction of ligands with biomolecules is an integral component of drug discovery and development. Challenges for computing thermodynamic and kinetic quantities for pharmaceutically relevant receptor-ligand complexes include the size and flexibility of the ligands, large-scale conformational rearrangements of the receptor, accurate force field parameters, simulation efficiency, and sufficient sampling associated with rare events. Our recently developed multiscale milestoning simulation approach, SEEKR2 (Simulation Enabled Estimation of Kinetic Rates v.2), has demonstrated success in predicting unbinding (koff) kinetics by employing molecular dynamics (MD) simulations in regions closer to the binding site. The MD region is further subdivided into smaller Voronoi tessellations to improve the simulation efficiency and parallelization. To date, all MD simulations are run using general molecular mechanics (MM) force fields. The accuracy of calculations can be further improved by incorporating quantum mechanical (QM) methods into generating system-specific force fields through reparameterizing ligand partial charges in the bound state. The force field reparameterization process modifies the potential energy landscape of the bimolecular complex, enabling a more accurate representation of the intermolecular interactions and polarization effects at the bound state. We present QMrebind (Quantum Mechanical force field reparameterization at the receptor-ligand binding site), an ORCA-based software that facilitates reparameterizing the potential energy function within the phase space representing the bound state in a receptor-ligand complex. With SEEKR2 koff estimates and experimentally determined kinetic rates, we compare and interpret the receptor-ligand unbinding kinetics obtained using the newly reparameterized force fields for model host-guest systems and HSP90-inhibitor complexes. This method provides an opportunity to achieve higher accuracy in predicting receptor-ligand koff rate constants.
Collapse
Affiliation(s)
- Anupam Anand Ojha
- Department of Chemistry and Biochemistry, University of California San Diego La Jolla California 92093 USA
| | - Lane William Votapka
- Department of Chemistry and Biochemistry, University of California San Diego La Jolla California 92093 USA
| | - Rommie Elizabeth Amaro
- Department of Molecular Biology, University of California San Diego La Jolla California 92093 USA
| |
Collapse
|
2
|
Iida S, Kameda T. Dissociation Rate Calculation via Constant-Force Steered Molecular Dynamics Simulation. J Chem Inf Model 2023. [PMID: 37188657 DOI: 10.1021/acs.jcim.2c01529] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/17/2023]
Abstract
Steered molecular dynamics (SMD) simulations are used to study molecular dissociation events by applying a harmonic force to the molecules and pulling them at a constant velocity. Instead of constant-velocity pulling, we use a constant force: the constant-force SMD (CF-SMD) simulation. The CF-SMD simulation employs a constant force to reduce the activation barrier of molecular dissociation, thereby enhancing the dissociation event. Here, we present the capability of the CF-SMD simulation to estimate the dissociation time at equilibrium. We performed all-atom CF-SMD simulations for NaCl and protein-ligand systems, producing dissociation time at various forces. We extrapolated these values to the dissociation rate without a constant force using Bell's model or the Dudko-Hummer-Szabo model. We demonstrate that the CF-SMD simulations with the models predicted the dissociation time in equilibrium. A CF-SMD simulation is a powerful tool for estimating the dissociation rate in a direct and computationally efficient manner.
Collapse
Affiliation(s)
- Shinji Iida
- Artificial Intelligence Research Center, National Institute of Advanced Industrial Science and Technology (AIST), 2-4-7 Aomi, Koto-ku, Tokyo 135-0064, Japan
| | - Tomoshi Kameda
- Artificial Intelligence Research Center, National Institute of Advanced Industrial Science and Technology (AIST), 2-4-7 Aomi, Koto-ku, Tokyo 135-0064, Japan
| |
Collapse
|
3
|
Muñiz‐Chicharro A, Votapka LW, Amaro RE, Wade RC. Brownian dynamics simulations of biomolecular diffusional association processes. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2022. [DOI: 10.1002/wcms.1649] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Affiliation(s)
- Abraham Muñiz‐Chicharro
- Molecular and Cellular Modeling Group Heidelberg Institute for Theoretical Studies (HITS) Heidelberg Germany
- Faculty of Biosciences and Heidelberg Graduate School of Mathematical and Computational Methods for the Sciences (HGS MathComp) Heidelberg University Heidelberg Germany
| | | | | | - Rebecca C. Wade
- Molecular and Cellular Modeling Group Heidelberg Institute for Theoretical Studies (HITS) Heidelberg Germany
- Center for Molecular Biology (ZMBH), DKFZ‐ZMBH Alliance, and Interdisciplinary Center for Scientific Computing (IWR) Heidelberg University Heidelberg Germany
| |
Collapse
|
4
|
Cardenas AE, Hunter A, Wang H, Elber R. ScMiles2: A Script to Conduct and Analyze Milestoning Trajectories for Long Time Dynamics. J Chem Theory Comput 2022; 18:6952-6965. [PMID: 36191005 PMCID: PMC10336853 DOI: 10.1021/acs.jctc.2c00708] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Milestoning is a theory and an algorithm that computes kinetics and thermodynamics at long time scales. It is based on partitioning the (phase) space into cells and running a large number of short trajectories between the boundaries of the cells. The termination points of the trajectories are analyzed with the Milestoning theory to obtain kinetic and thermodynamic information. Managing the tens to hundreds of thousands of Milestoning trajectories is a challenge, which we handle with a python script, ScMiles. Here, we introduce a new version of the python script ScMiles2 to conduct Milestoning simulations. Major enhancements are: (i) post analysis of Milestoning trajectories to obtain the free energy, mean first passage time, the committor function, and exit times; (ii) similar to (i) but the post analysis is for a single long trajectory; (iii) we support the use of the GROMACS software in addition to NAMD; (iv) a restart option; (v) the automated finding, sampling, and launching trajectories from new milestones that are found on the fly; and (vi) support Milestoning calculations with several coarse variables and for complex reaction coordinates. We also evaluate the simulation parameters and suggest new algorithmic features to enhance the rate of convergence of observables. We propose the use of an iteration-averaged kinetic matrix for a rapid approach to asymptotic values. Illustrations are provided for small systems and one large example.
Collapse
Affiliation(s)
- Alfredo E. Cardenas
- The Oden Institute, University of Texas at Austin, Austin, Texas, 78712, USA
| | - Allison Hunter
- The Oden Institute, University of Texas at Austin, Austin, Texas, 78712, USA
| | - Hao Wang
- The Oden Institute, University of Texas at Austin, Austin, Texas, 78712, USA
- Qingdao Institute for Theoretical and Computational Sciences, Institute of Frontier and Interdisciplinary Science, Shandong University, Qingdao, Shandong 266237, China
| | - Ron Elber
- The Oden Institute, University of Texas at Austin, Austin, Texas, 78712, USA
- Department of Chemistry, University of Texas at Austin, Austin, Texas, 78712, USA
| |
Collapse
|
5
|
Shekhar M, Smith Z, Seeliger MA, Tiwary P. Protein Flexibility and Dissociation Pathway Differentiation Can Explain Onset of Resistance Mutations in Kinases. Angew Chem Int Ed Engl 2022; 61:e202200983. [PMID: 35486370 DOI: 10.1002/anie.202200983] [Citation(s) in RCA: 15] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2022] [Indexed: 12/14/2022]
Abstract
Understanding how mutations render a drug ineffective is a problem of immense relevance. Often the mechanism through which mutations cause drug resistance can be explained purely through thermodynamics. However, the more perplexing situation is when two proteins have the same drug binding affinities but different residence times. In this work, we demonstrate how all-atom molecular dynamics simulations using recent developments grounded in statistical mechanics can provide a detailed mechanistic rationale for such variances. We discover dissociation mechanisms for the anti-cancer drug Imatinib (Gleevec) against wild-type and the N368S mutant of Abl kinase. We show how this point mutation triggers far-reaching changes in the protein's flexibility and leads to a different, much faster, drug dissociation pathway. We believe that this work marks an efficient and scalable approach to obtain mechanistic insight into resistance mutations in biomolecular receptors that are hard to explain using a structural perspective.
Collapse
Affiliation(s)
- Mrinal Shekhar
- Center for Development of Therapeutics, Broad Institute of MIT and Harvard, Cambridge, MA 02142, USA
| | - Zachary Smith
- Biophysics Program and Institute for Physical Science and Technology, University of Maryland, College Park, MD 20742, USA
| | - Markus A Seeliger
- Department of Pharmacological Sciences, Stony Brook University, Stony Brook, NY 11794-8651, USA
| | - Pratyush Tiwary
- Department of Chemistry and Biochemistry and Institute for Physical Science and Technology, University of Maryland, College Park, MD 20742, USA
| |
Collapse
|
6
|
Votapka LW, Stokely AM, Ojha AA, Amaro RE. SEEKR2: Versatile Multiscale Milestoning Utilizing the OpenMM Molecular Dynamics Engine. J Chem Inf Model 2022; 62:3253-3262. [PMID: 35759413 PMCID: PMC9277580 DOI: 10.1021/acs.jcim.2c00501] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
We present SEEKR2
(simulation-enabled estimation of kinetic rates
version 2)—the latest iteration in the family of SEEKR programs
for using multiscale simulation methods to computationally estimate
the kinetics and thermodynamics of molecular processes, in particular,
ligand-receptor binding. SEEKR2 generates equivalent, or improved,
results compared to the earlier versions of SEEKR but with significant
increases in speed and capabilities. SEEKR2 has also been built with
greater ease of usability and with extensible features to enable future
expansions of the method. Now, in addition to supporting simulations
using NAMD, calculations may be run with the fast and extensible OpenMM
simulation engine. The Brownian dynamics portion of the calculation
has also been upgraded to Browndye 2. Furthermore, this version of
SEEKR supports hydrogen mass repartitioning, which significantly reduces
computational cost, while showing little, if any, loss of accuracy
in the predicted kinetics.
Collapse
Affiliation(s)
- Lane W Votapka
- University of California, San Diego, 9500 Gilman Dr., La Jolla, California 92093, United States
| | - Andrew M Stokely
- University of California, San Diego, 9500 Gilman Dr., La Jolla, California 92093, United States
| | - Anupam A Ojha
- University of California, San Diego, 9500 Gilman Dr., La Jolla, California 92093, United States
| | - Rommie E Amaro
- University of California, San Diego, 9500 Gilman Dr., La Jolla, California 92093, United States
| |
Collapse
|
7
|
Kahana A, Lancet D, Palmai Z. Micellar Composition Affects Lipid Accretion Kinetics in Molecular Dynamics Simulations: Support for Lipid Network Reproduction. Life (Basel) 2022; 12:955. [PMID: 35888044 PMCID: PMC9325298 DOI: 10.3390/life12070955] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/14/2022] [Revised: 06/02/2022] [Accepted: 06/21/2022] [Indexed: 11/25/2022] Open
Abstract
Mixed lipid micelles were proposed to facilitate life through their documented growth dynamics and catalytic properties. Our previous research predicted that micellar self-reproduction involves catalyzed accretion of lipid molecules by the residing lipids, leading to compositional homeostasis. Here, we employ atomistic Molecular Dynamics simulations, beginning with 54 lipid monomers, tracking an entire course of micellar accretion. This was done to examine the self-assembly of variegated lipid clusters, allowing us to measure entry and exit rates of monomeric lipids into pre-micelles with different compositions and sizes. We observe considerable rate-modifications that depend on the assembly composition and scrutinize the underlying mechanisms as well as the energy contributions. Lastly, we describe the measured potential for compositional homeostasis in our simulated mixed micelles. This affirms the basis for micellar self-reproduction, with implications for the study of the origin of life.
Collapse
Affiliation(s)
| | | | - Zoltan Palmai
- Department of Molecular Genetics, Weizmann Institute of Science, Rehovot 761001, Israel; (A.K.); (D.L.)
| |
Collapse
|
8
|
Shekhar M, Smith Z, Seeliger M, Tiwary P. Protein Flexibility and Dissociation Pathway Differentiation Can Explain Onset Of Resistance Mutations in Kinases. Angew Chem Int Ed Engl 2022. [DOI: 10.1002/ange.202200983] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Affiliation(s)
- Mrinal Shekhar
- Broad Institute Center for Development of Therapeutics UNITED STATES
| | - Zachary Smith
- University of Maryland at College Park Institute for Physical Science and Technology UNITED STATES
| | - Markus Seeliger
- Stony Brook University Department of Pharmacological Sciences UNITED STATES
| | - Pratyush Tiwary
- university of maryland chemistry and biochemistry university of maryland 20740 college park UNITED STATES
| |
Collapse
|
9
|
Ge Y, Voelz VA. Estimation of binding rates and affinities from multiensemble Markov models and ligand decoupling. J Chem Phys 2022; 156:134115. [PMID: 35395889 PMCID: PMC8993428 DOI: 10.1063/5.0088024] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Accurate and efficient simulation of the thermodynamics and kinetics of protein-ligand interactions is crucial for computational drug discovery. Multiensemble Markov Model (MEMM) estimators can provide estimates of both binding rates and affinities from collections of short trajectories but have not been systematically explored for situations when a ligand is decoupled through scaling of non-bonded interactions. In this work, we compare the performance of two MEMM approaches for estimating ligand binding affinities and rates: (1) the transition-based reweighting analysis method (TRAM) and (2) a Maximum Caliber (MaxCal) based method. As a test system, we construct a small host-guest system where the ligand is a single uncharged Lennard-Jones (LJ) particle, and the receptor is an 11-particle icosahedral pocket made from the same atom type. To realistically mimic a protein-ligand binding system, the LJ ϵ parameter was tuned, and the system was placed in a periodic box with 860 TIP3P water molecules. A benchmark was performed using over 80 µs of unbiased simulation, and an 18-state Markov state model was used to estimate reference binding affinities and rates. We then tested the performance of TRAM and MaxCal when challenged with limited data. Both TRAM and MaxCal approaches perform better than conventional Markov state models, with TRAM showing better convergence and accuracy. We find that subsampling of trajectories to remove time correlation improves the accuracy of both TRAM and MaxCal and that in most cases, only a single biased ensemble to enhance sampled transitions is required to make accurate estimates.
Collapse
Affiliation(s)
- Yunhui Ge
- Department of Pharmaceutical Sciences, University of California, Irvine, California 92697, USA
| | - Vincent A Voelz
- Department of Chemistry, Temple University, Philadelphia, Pennsylvania 19122, USA
| |
Collapse
|
10
|
Ray D, Stone SE, Andricioaei I. Markovian Weighted Ensemble Milestoning (M-WEM): Long-Time Kinetics from Short Trajectories. J Chem Theory Comput 2021; 18:79-95. [PMID: 34910499 DOI: 10.1021/acs.jctc.1c00803] [Citation(s) in RCA: 17] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We introduce a rare-event sampling scheme, named Markovian Weighted Ensemble Milestoning (M-WEM), which inlays a weighted ensemble framework within a Markovian milestoning theory to efficiently calculate thermodynamic and kinetic properties of long-time-scale biomolecular processes from short atomistic molecular dynamics simulations. M-WEM is tested on the Müller-Brown potential model, the conformational switching in alanine dipeptide, and the millisecond time-scale protein-ligand unbinding in a trypsin-benzamidine complex. Not only can M-WEM predict the kinetics of these processes with quantitative accuracy but it also allows for a scheme to reconstruct a multidimensional free-energy landscape along additional degrees of freedom, which are not part of the milestoning progress coordinate. For the ligand-receptor system, the experimental residence time, association and dissociation kinetics, and binding free energy could be reproduced using M-WEM within a simulation time of a few hundreds of nanoseconds, which is a fraction of the computational cost of other currently available methods, and close to 4 orders of magnitude less than the experimental residence time. Due to the high accuracy and low computational cost, the M-WEM approach can find potential applications in kinetics and free-energy-based computational drug design.
Collapse
Affiliation(s)
- Dhiman Ray
- Department of Chemistry, University of California Irvine, Irvine, California 92697, United States
| | - Sharon Emily Stone
- Department of Chemistry, University of California Irvine, Irvine, California 92697, United States
| | - Ioan Andricioaei
- Department of Chemistry, University of California Irvine, Irvine, California 92697, United States.,Department of Physics and Astronomy, University of California Irvine, Irvine, California 92697, United States
| |
Collapse
|
11
|
Kasahara K, Masayama R, Okita K, Matubayasi N. Atomistic description of molecular binding processes based on returning probability theory. J Chem Phys 2021; 155:204503. [PMID: 34852475 DOI: 10.1063/5.0070308] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The efficiency of molecular binding such as host-guest binding is commonly evaluated in terms of kinetics, such as rate coefficients. In general, to compute the coefficient of the overall binding process, we need to consider both the diffusion of reactants and barrier crossing to reach the bound state. Here, we develop a methodology of quantifying the rate coefficient of binding based on molecular dynamics simulation and returning probability (RP) theory proposed by Kim and Lee [J. Chem. Phys. 131, 014503 (2009)]. RP theory provides a tractable formula of the rate coefficient in terms of the thermodynamic stability and kinetics of the intermediate state on a predefined reaction coordinate. In this study, the interaction energy between reactants is utilized as the reaction coordinate, enabling us to effectively describe the reactants' relative position and orientation on one-dimensional space. Application of this method to the host-guest binding systems, which consist of β-cyclodextrin and small guest molecules, yields the rate coefficients consistent with the experimental results.
Collapse
Affiliation(s)
- Kento Kasahara
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| | - Ren Masayama
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| | - Kazuya Okita
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| | - Nobuyuki Matubayasi
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| |
Collapse
|
12
|
Sadiq SK, Muñiz Chicharro A, Friedrich P, Wade RC. Multiscale Approach for Computing Gated Ligand Binding from Molecular Dynamics and Brownian Dynamics Simulations. J Chem Theory Comput 2021; 17:7912-7929. [PMID: 34739248 DOI: 10.1021/acs.jctc.1c00673] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
We develop an approach to characterize the effects of gating by a multiconformation protein consisting of macrostate conformations that are either accessible or inaccessible to ligand binding. We first construct a Markov state model of the apo-protein from atomistic molecular dynamics simulations from which we identify macrostates and their conformations, compute their relative macrostate populations and interchange kinetics, and structurally characterize them in terms of ligand accessibility. We insert the calculated first-order rate constants for conformational transitions into a multistate gating theory from which we derive a gating factor γ that quantifies the degree of conformational gating. Applied to HIV-1 protease, our approach yields a kinetic network of three accessible (semi-open, open, and wide-open) and two inaccessible (closed and a newly identified, "parted") macrostate conformations. The parted conformation sterically partitions the active site, suggesting a possible role in product release. We find that the binding kinetics of drugs and drug-like inhibitors to HIV-1 protease falls in the slow gating regime. However, because γ = 0.75, conformational gating only modestly slows ligand binding. Brownian dynamics simulations of the diffusional association of eight inhibitors to the protease─having a wide range of experimental association constants (∼104-1010 M-1 s-1)─yields gated rate constants in the range of ∼0.5-5.7 × 108 M-1 s-1. This indicates that, whereas the association rate of some inhibitors could be described by the model, for many inhibitors either subsequent conformational transitions or alternate binding mechanisms may be rate-limiting. For systems known to be modulated by conformational gating, the approach could be scaled computationally efficiently to screen association kinetics for a large number of ligands.
Collapse
Affiliation(s)
- S Kashif Sadiq
- Molecular and Cellular Modeling Group, Heidelberg Institute for Theoretical Studies (HITS), Schloss-Wolfsbrunnenweg 35, 69118 Heidelberg, Germany.,Genome Biology Unit, European Molecular Biology Laboratory, Meyerhofstrasse 1, 69117 Heidelberg, Germany.,Infection Biology Unit, Universitat Pompeu Fabra, Barcelona Biomedical Research Park (PRBB), C/Doctor Aiguader 88, 08003 Barcelona, Spain
| | - Abraham Muñiz Chicharro
- Molecular and Cellular Modeling Group, Heidelberg Institute for Theoretical Studies (HITS), Schloss-Wolfsbrunnenweg 35, 69118 Heidelberg, Germany.,Faculty of Biosciences, Heidelberg University, Im Neuenheimer Feld 234, 69120 Heidelberg, Germany
| | - Patrick Friedrich
- Molecular and Cellular Modeling Group, Heidelberg Institute for Theoretical Studies (HITS), Schloss-Wolfsbrunnenweg 35, 69118 Heidelberg, Germany
| | - Rebecca C Wade
- Molecular and Cellular Modeling Group, Heidelberg Institute for Theoretical Studies (HITS), Schloss-Wolfsbrunnenweg 35, 69118 Heidelberg, Germany.,Center for Molecular Biology (ZMBH), DKFZ-ZMBH Alliance, Heidelberg University, Im Neuenheimer Feld 282, 69120 Heidelberg, Germany.,Interdisciplinary Center for Scientific Computing (IWR), Heidelberg University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| |
Collapse
|
13
|
Maximova E, Postnikov EB, Lavrova AI, Farafonov V, Nerukh D. Protein-Ligand Dissociation Rate Constant from All-Atom Simulation. J Phys Chem Lett 2021; 12:10631-10636. [PMID: 34704768 DOI: 10.1021/acs.jpclett.1c02952] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
Dissociation of a ligand isoniazid from a protein catalase was investigated using all-atom molecular dynamics (MD) simulations. Random acceleration MD (τ-RAMD) was used, in which a random artificial force applied to the ligand facilitates its dissociation. We have suggested a novel approach to extrapolate such obtained dissociation times to the zero-force limit assuming never before attempted universal exponential dependence of the bond strength on the applied force, allowing direct comparison with experimentally measured values. We have found that our calculated dissociation time was equal to 36.1 s with statistically significant values distributed in the interval of 0.2-72.0 s, which quantitatively matches the experimental value of 50 ± 8 s despite the extrapolation over 9 orders of magnitude in time.
Collapse
Affiliation(s)
- Ekaterina Maximova
- Department of Nanobiotechnology, Alferov University, Khlopina Street, 8/3 A, 194021 Saint Petersburg, Russia
- Center for Molecular and Cellular Bioengineering, Technische Universität Dresden, Helmholtzstr. 10, 01069 Dresden, Germany
| | - Eugene B Postnikov
- Department of Theoretical Physics, Kursk State University, Radishcheva Street, 33, 305000 Kursk, Russia
| | - Anastasia I Lavrova
- Saint-Petersburg State University, 7/9 Universitetskaya Emb., 199034 Saint Petersburg, Russia
- Saint-Petersburg State Research Institute of Phthisiopulmonology, 2-4 Ligovskiy Avenue, 194064 Saint-Petersburg, Russia
| | - Vladimir Farafonov
- V. N. Karazin Kharkiv National University, 4 Svobody sq., Kharkiv 61022, Ukraine
| | - Dmitry Nerukh
- Department of Mathematics, Aston University, Birmingham B4 7ET, U.K
| |
Collapse
|
14
|
Edwards T, Foloppe N, Harris SA, Wells G. The future of biomolecular simulation in the pharmaceutical industry: what we can learn from aerodynamics modelling and weather prediction. Part 1. understanding the physical and computational complexity of in silico drug design. Acta Crystallogr D Struct Biol 2021; 77:1348-1356. [PMID: 34726163 PMCID: PMC8561735 DOI: 10.1107/s2059798321009712] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/25/2020] [Accepted: 09/17/2021] [Indexed: 02/04/2023] Open
Abstract
The predictive power of simulation has become embedded in the infrastructure of modern economies. Computer-aided design is ubiquitous throughout industry. In aeronautical engineering, built infrastructure and materials manufacturing, simulations are routinely used to compute the performance of potential designs before construction. The ability to predict the behaviour of products is a driver of innovation by reducing the cost barrier to new designs, but also because radically novel ideas can be piloted with relatively little risk. Accurate weather forecasting is essential to guide domestic and military flight paths, and therefore the underpinning simulations are critical enough to have implications for national security. However, in the pharmaceutical and biotechnological industries, the application of computer simulations remains limited by the capabilities of the technology with respect to the complexity of molecular biology and human physiology. Over the last 30 years, molecular-modelling tools have gradually gained a degree of acceptance in the pharmaceutical industry. Drug discovery has begun to benefit from physics-based simulations. While such simulations have great potential for improved molecular design, much scepticism remains about their value. The motivations for such reservations in industry and areas where simulations show promise for efficiency gains in preclinical research are discussed. In this, the first of two complementary papers, the scientific and technical progress that needs to be made to improve the predictive power of biomolecular simulations, and how this might be achieved, is firstly discussed (Part 1). In Part 2, the status of computer simulations in pharma is contrasted with aerodynamics modelling and weather forecasting, and comments are made on the cultural changes needed for equivalent computational technologies to become integrated into life-science industries.
Collapse
Affiliation(s)
- Tom Edwards
- School of Molecular and Cellular Biology, University of Leeds, Leeds, United Kingdom
- Astbury Centre for Structural and Molecular Biology, University of Leeds, Leeds, United Kingdom
| | | | - Sarah Anne Harris
- Astbury Centre for Structural and Molecular Biology, University of Leeds, Leeds, United Kingdom
- School of Physics and Astronomy, University of Leeds, Leeds, United Kingdom
| | - Geoff Wells
- School of Pharmacy, University College London, London, United Kingdom
| |
Collapse
|
15
|
Del Razo MJ, Dibak M, Schütte C, Noé F. Multiscale molecular kinetics by coupling Markov state models and reaction-diffusion dynamics. J Chem Phys 2021; 155:124109. [PMID: 34598578 DOI: 10.1063/5.0060314] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A novel approach to simulate simple protein-ligand systems at large time and length scales is to couple Markov state models (MSMs) of molecular kinetics with particle-based reaction-diffusion (RD) simulations, MSM/RD. Currently, MSM/RD lacks a mathematical framework to derive coupling schemes, is limited to isotropic ligands in a single conformational state, and lacks multiparticle extensions. In this work, we address these needs by developing a general MSM/RD framework by coarse-graining molecular dynamics into hybrid switching diffusion processes. Given enough data to parameterize the model, it is capable of modeling protein-protein interactions over large time and length scales, and it can be extended to handle multiple molecules. We derive the MSM/RD framework, and we implement and verify it for two protein-protein benchmark systems and one multiparticle implementation to model the formation of pentameric ring molecules. To enable reproducibility, we have published our code in the MSM/RD software package.
Collapse
Affiliation(s)
- Mauricio J Del Razo
- Van 't Hoff Institute for Molecular Sciences, University of Amsterdam, Amsterdam, The Netherlands
| | - Manuel Dibak
- Department of Mathematics and Computer Science, Freie Universität Berlin, Berlin, Germany
| | | | - Frank Noé
- Department of Physics, Freie Universität Berlin, Berlin, Germany
| |
Collapse
|
16
|
Castelli M, Serapian SA, Marchetti F, Triveri A, Pirota V, Torielli L, Collina S, Doria F, Freccero M, Colombo G. New perspectives in cancer drug development: computational advances with an eye to design. RSC Med Chem 2021; 12:1491-1502. [PMID: 34671733 PMCID: PMC8459323 DOI: 10.1039/d1md00192b] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/09/2021] [Accepted: 07/06/2021] [Indexed: 02/06/2023] Open
Abstract
Computational chemistry has come of age in drug discovery. Indeed, most pharmaceutical development programs rely on computer-based data and results at some point. Herein, we discuss recent applications of advanced simulation techniques to difficult challenges in drug discovery. These entail the characterization of allosteric mechanisms and the identification of allosteric sites or cryptic pockets determined by protein motions, which are not immediately evident in the experimental structure of the target; the study of ligand binding mechanisms and their kinetic profiles; and the evaluation of drug-target affinities. We analyze different approaches to tackle challenging and emerging biological targets. Finally, we discuss the possible perspectives of future application of computation in drug discovery.
Collapse
Affiliation(s)
- Matteo Castelli
- Department of Chemistry, University of Pavia via Taramelli 12 27100 Pavia Italy
| | - Stefano A Serapian
- Department of Chemistry, University of Pavia via Taramelli 12 27100 Pavia Italy
| | - Filippo Marchetti
- Department of Chemistry, University of Pavia via Taramelli 12 27100 Pavia Italy
| | - Alice Triveri
- Department of Chemistry, University of Pavia via Taramelli 12 27100 Pavia Italy
| | - Valentina Pirota
- Department of Chemistry, University of Pavia via Taramelli 12 27100 Pavia Italy
| | - Luca Torielli
- Department of Drug Sciences, Medicinal Chemistry and Pharmaceutical Technology Section, University of Pavia via Taramelli 12 27100 Pavia Italy
| | - Simona Collina
- Department of Drug Sciences, Medicinal Chemistry and Pharmaceutical Technology Section, University of Pavia via Taramelli 12 27100 Pavia Italy
| | - Filippo Doria
- Department of Chemistry, University of Pavia via Taramelli 12 27100 Pavia Italy
| | - Mauro Freccero
- Department of Chemistry, University of Pavia via Taramelli 12 27100 Pavia Italy
| | - Giorgio Colombo
- Department of Chemistry, University of Pavia via Taramelli 12 27100 Pavia Italy
| |
Collapse
|
17
|
Narayan B, Buchete NV, Elber R. Computer Simulations of the Dissociation Mechanism of Gleevec from Abl Kinase with Milestoning. J Phys Chem B 2021; 125:5706-5715. [PMID: 33930271 DOI: 10.1021/acs.jpcb.1c00264] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Gleevec (a.k.a., imatinib) is an important anticancer (e.g., chronic myeloid leukemia) chemotherapeutic drug due to its inhibitory interaction with the Abl kinase. Here, we use atomically detailed simulations within the Milestoning framework to study the molecular dissociation mechanism of Gleevec from Abl kinase. We compute the dissociation free energy profile, the mean first passage time for unbinding, and explore the transition state ensemble of conformations. The milestones form a multidimensional network with average connectivity of about 2.93, which is significantly higher than the connectivity for a one-dimensional reaction coordinate. The free energy barrier for Gleevec dissociation is estimated to be ∼10 kcal/mol, and the exit time is ∼55 ms. We examined the transition state conformations using both, the committor and transition function. We show that near the transition state the highly conserved salt bridge K217 and E286 is transiently broken. Together with the calculated free energy profile, these calculations can advance the understanding of the molecular interaction mechanisms between Gleevec and Abl kinase and play a role in future drug design and optimization studies.
Collapse
Affiliation(s)
- Brajesh Narayan
- School of Physics, University College Dublin, Belfield, Dublin 4, Ireland.,Institute for Discovery, University College Dublin, Belfield, Dublin 4, Ireland
| | - Nicolae-Viorel Buchete
- School of Physics, University College Dublin, Belfield, Dublin 4, Ireland.,Institute for Discovery, University College Dublin, Belfield, Dublin 4, Ireland
| | - Ron Elber
- Oden Institute for Computational Engineering and Science, Department of Chemistry, University of Texas at Austin, Austin Texas 78712, United States
| |
Collapse
|
18
|
Abstract
Molecular dynamics simulations can now routinely access the microsecond timescale, making feasible direct sampling of ligand association events. While Markov State Model (MSM) approaches offer a useful framework for analyzing such trajectory data to gain insight into binding mechanisms, accurate modeling of ligand association pathways and kinetics must be done carefully. We describe methods and good practices for constructing MSMs of ligand binding from unbiased trajectory data and discuss how to use time-lagged independent component analysis (tICA) to build informative models, using as an example recent simulation work to model the binding of phenylalanine to the regulatory ACT domain dimer of phenylalanine hydroxylase. We describe a variety of methods for estimating association rates from MSMs and discuss how to distinguish between conformational selection and induced-fit mechanisms using MSMs. In addition, we review some examples of MSMs constructed to elucidate the mechanisms by which p53 transactivation domain (TAD) and related peptides bind the oncoprotein MDM2.
Collapse
Affiliation(s)
- Yunhui Ge
- Department of Chemistry, Temple University, Philadelphia, PA, USA
| | - Vincent A Voelz
- Department of Chemistry, Temple University, Philadelphia, PA, USA.
| |
Collapse
|
19
|
Elber R, Fathizadeh A, Ma P, Wang H. Modeling molecular kinetics with Milestoning. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2020. [DOI: 10.1002/wcms.1512] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Affiliation(s)
- Ron Elber
- Department of Chemistry, The Oden Institute for Computational Engineering and Sciences University of Texas at Austin Austin Texas USA
| | - Arman Fathizadeh
- The Oden Institute for Computational Engineering and Sciences University of Texas at Austin Austin Texas USA
| | - Piao Ma
- Department of Chemistry University of Texas at Austin Austin Texas USA
| | - Hao Wang
- The Oden Institute for Computational Engineering and Sciences University of Texas at Austin Austin Texas USA
| |
Collapse
|
20
|
Ahn SH, Jagger BR, Amaro RE. Ranking of Ligand Binding Kinetics Using a Weighted Ensemble Approach and Comparison with a Multiscale Milestoning Approach. J Chem Inf Model 2020; 60:5340-5352. [PMID: 32315175 DOI: 10.1021/acs.jcim.9b00968] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
To improve lead optimization efforts in finding the right ligand, pharmaceutical industries need to know the ligand's binding kinetics, such as binding and unbinding rate constants, which often correlate with the ligand's efficacy in vivo. To predict binding kinetics efficiently, enhanced sampling methods, such as milestoning and the weighted ensemble (WE) method, have been used in molecular dynamics (MD) simulations of these systems. However, a comparison of these enhanced sampling methods in ranking ligands has not been done. Hence, a WE approach called the concurrent adaptive sampling (CAS) algorithm that uses MD simulations was used to rank seven ligands for β-cyclodextrin, a system in which a multiscale milestoning approach called simulation enabled estimation of kinetic rates (SEEKR) was also used, which uses both MD and Brownian dynamics simulations. Overall, the CAS algorithm can successfully rank ligands using the unbinding rate constant koff values and binding free energy ΔG values, as SEEKR did, with reduced computational cost that is about the same as SEEKR. We compare the CAS algorithm simulations with different parameters and discuss the impact of parameters in ranking ligands and obtaining rate constant and binding free energy estimates. We also discuss similarities and differences and advantages and disadvantages of SEEKR and the CAS algorithm for future use.
Collapse
Affiliation(s)
- Surl-Hee Ahn
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Benjamin R Jagger
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Rommie E Amaro
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| |
Collapse
|
21
|
Ray D, Gokey T, Mobley DL, Andricioaei I. Kinetics and free energy of ligand dissociation using weighted ensemble milestoning. J Chem Phys 2020; 153:154117. [PMID: 33092382 DOI: 10.1063/5.0021953] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We consider the recently developed weighted ensemble milestoning (WEM) scheme [D. Ray and I. Andricioaei, J. Chem. Phys. 152, 234114 (2020)] and test its capability of simulating ligand-receptor dissociation dynamics. We performed WEM simulations on the following host-guest systems: Na+/Cl- ion pair and 4-hydroxy-2-butanone ligand with FK506 binding protein. As a proof of principle, we show that the WEM formalism reproduces the Na+/Cl- ion pair dissociation timescale and the free energy profile obtained from long conventional MD simulation. To increase the accuracy of WEM calculations applied to kinetics and thermodynamics in protein-ligand binding, we introduced a modified WEM scheme called weighted ensemble milestoning with restraint release (WEM-RR), which can increase the number of starting points per milestone without adding additional computational cost. WEM-RR calculations obtained a ligand residence time and binding free energy in agreement with experimental and previous computational results. Moreover, using the milestoning framework, the binding time and rate constants, dissociation constants, and committor probabilities could also be calculated at a low computational cost. We also present an analytical approach for estimating the association rate constant (kon) when binding is primarily diffusion driven. We show that the WEM method can efficiently calculate multiple experimental observables describing ligand-receptor binding/unbinding and is a promising candidate for computer-aided inhibitor design.
Collapse
Affiliation(s)
- Dhiman Ray
- Department of Chemistry, University of California Irvine, Irvine, California 92697, USA
| | - Trevor Gokey
- Department of Chemistry, University of California Irvine, Irvine, California 92697, USA
| | - David L Mobley
- Department of Chemistry, University of California Irvine, Irvine, California 92697, USA
| | - Ioan Andricioaei
- Department of Chemistry, University of California Irvine, Irvine, California 92697, USA
| |
Collapse
|
22
|
Jagger BR, Ojha AA, Amaro RE. Predicting Ligand Binding Kinetics Using a Markovian Milestoning with Voronoi Tessellations Multiscale Approach. J Chem Theory Comput 2020; 16:5348-5357. [DOI: 10.1021/acs.jctc.0c00495] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
- Benjamin R. Jagger
- Department of Chemistry and Biochemistry, University of California San Diego, 9500 Gilman Drive, La Jolla, California 92093, United States
| | - Anupam A. Ojha
- Department of Chemistry and Biochemistry, University of California San Diego, 9500 Gilman Drive, La Jolla, California 92093, United States
| | - Rommie E. Amaro
- Department of Chemistry and Biochemistry, University of California San Diego, 9500 Gilman Drive, La Jolla, California 92093, United States
| |
Collapse
|
23
|
Wang H, Huang N, Dangerfield T, Johnson KA, Gao J, Elber R. Exploring the Reaction Mechanism of HIV Reverse Transcriptase with a Nucleotide Substrate. J Phys Chem B 2020; 124:4270-4283. [PMID: 32364738 PMCID: PMC7260111 DOI: 10.1021/acs.jpcb.0c02632] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
Abstract
Enzymatic reactions consist of several steps: (i) a weak binding event of the substrate to the enzyme, (ii) an induced fit or a protein conformational transition upon ligand binding, (iii) the chemical reaction, and (iv) the release of the product. Here we focus on step iii of the reaction of a DNA polymerase, HIV RT, with a nucleotide. We determine the rate and the free energy profile for the addition of a nucleotide to a DNA strand using a combination of a QM/MM model, the string method, and exact Milestoning. The barrier height and the time scale of the reaction are consistent with experiment. We show that the observables (free energies and mean first passage time) converge rapidly, as a function of the Milestoning iteration number. We also consider the substitution of an oxygen of the incoming nucleotide by a nonbridging sulfur atom and its impact on the enzymatic reaction. This substitution has been suggested in the past as a tool to examine the influence of the chemical step on the overall rate. Our joint computational and experimental study suggests that the impact of the substitution is small. Computationally, the differences between the two are within the estimated error bars. Experiments suggest a small difference. Finally, we examine step i, the weak binding of the nucleotide to the protein surface. We suggest that this step has only a small contribution to the selectivity of the enzyme. Comments are made on the impact of these steps on the overall mechanism.
Collapse
Affiliation(s)
- Hao Wang
- Oden Institute for Computational Engineering and Sciences, University of Texas at Austin, Austin TX 78712
| | - Nathan Huang
- Department of Molecular Biosciences, University of Texas at Austin, Austin, TX 78712
| | - Tyler Dangerfield
- Department of Molecular Biosciences, University of Texas at Austin, Austin, TX 78712
| | - Kenneth A. Johnson
- Department of Molecular Biosciences, University of Texas at Austin, Austin, TX 78712
| | - Jiali Gao
- Department of Chemistry, University of Minnesota, Minneapolis, MN, 55455-0431
| | - Ron Elber
- Oden Institute for Computational Engineering and Sciences, University of Texas at Austin, Austin TX 78712
- Department of Chemistry, University of Texas at Austin, Austin, TX 78712
| |
Collapse
|
24
|
Jagger BR, Kochanek SE, Haldar S, Amaro RE, Mulholland AJ. Multiscale simulation approaches to modeling drug-protein binding. Curr Opin Struct Biol 2020; 61:213-221. [PMID: 32113133 DOI: 10.1016/j.sbi.2020.01.014] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/15/2019] [Revised: 01/16/2020] [Accepted: 01/21/2020] [Indexed: 01/19/2023]
Abstract
Simulations can provide detailed insight into the molecular processes involved in drug action, such as protein-ligand binding, and can therefore be a valuable tool for drug design and development. Processes with a large range of length and timescales may be involved, and understanding these different scales typically requires different types of simulation methodology. Ideally, simulations should be able to connect across scales, to analyze and predict how changes at one scale can influence another. Multiscale simulation methods, which combine different levels of treatment, are an emerging frontier with great potential in this area. Here we review multiscale frameworks of various types, and selected applications to biomolecular systems with a focus on drug-ligand binding.
Collapse
Affiliation(s)
- Benjamin R Jagger
- Department of Chemistry and Biochemistry, University of California San Diego, 9500 Gilman Drive, La Jolla, CA 92093, United States
| | - Sarah E Kochanek
- Department of Chemistry and Biochemistry, University of California San Diego, 9500 Gilman Drive, La Jolla, CA 92093, United States
| | - Susanta Haldar
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol, BS8 1TS, UK
| | - Rommie E Amaro
- Department of Chemistry and Biochemistry, University of California San Diego, 9500 Gilman Drive, La Jolla, CA 92093, United States.
| | - Adrian J Mulholland
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol, BS8 1TS, UK.
| |
Collapse
|
25
|
Tang Z, Chen SH, Chang CEA. Transient States and Barriers from Molecular Simulations and the Milestoning Theory: Kinetics in Ligand-Protein Recognition and Compound Design. J Chem Theory Comput 2020; 16:1882-1895. [PMID: 32031801 DOI: 10.1021/acs.jctc.9b01153] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
This study presents a novel computational approach to study molecular recognition and binding kinetics for drug-like compounds dissociating from a flexible protein system. The intermediates and their free energy profile during ligand association and dissociation processes control ligand-protein binding kinetics and bring a more complete picture of ligand-protein binding. The method applied the milestoning theory to extract kinetics and thermodynamics information from running short classical molecular dynamics (MD) simulations for frames from a given dissociation path. High-dimensional ligand-protein motions (3N-6 degrees of freedom) during ligand dissociation were reduced by use of principal component modes for assigning more than 100 milestones, and classical MD runs were allowed to travel multiple milestones to efficiently obtain ensemble distribution of initial structures for MD simulations and estimate the transition time and rate during ligand traveling between milestones. We used five pyrazolourea ligands and cyclin-dependent kinase 8 with cyclin C (CDK8/CycC) as our model system as well as metadynamics and a pathway search method to sample dissociation pathways. With our strategy, we constructed the free energy profile for highly mobile biomolecular systems. The computed binding free energy and residence time correctly ranked the pyrazolourea ligand series, in agreement with experimental data. Guided by a barrier of a ligand passing an αC helix and activation loop, we introduced one hydroxyl group to parent compounds to design our ligands with increased residence time and validated our prediction by experiments. This work provides a novel and robust approach to investigate dissociation kinetics of large and flexible systems for understanding unbinding mechanisms and designing new small-molecule drugs with desired binding kinetics.
Collapse
Affiliation(s)
- Zhiye Tang
- Department of Chemistry, University of California Riverside, Riverside, California 92521, United States
| | - Si-Han Chen
- Department of Chemistry, University of California Riverside, Riverside, California 92521, United States
| | - Chia-En A Chang
- Department of Chemistry, University of California Riverside, Riverside, California 92521, United States
| |
Collapse
|
26
|
Wei W, Elber R. ScMile: A Script to Investigate Kinetics with Short Time Molecular Dynamics Trajectories and the Milestoning Theory. J Chem Theory Comput 2020; 16:860-874. [PMID: 31922745 PMCID: PMC7031965 DOI: 10.1021/acs.jctc.9b01030] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Abstract
Studies of complex and rare events in condensed phase systems continue to attract considerable attention. Milestoning is a useful theory and algorithm to investigate the long-time dynamics of activated molecular events. It is based on launching a large number of short trajectories and statistical analysis of the outcome. The implementation of the theory in a computer script is described that enables more efficient Milestoning calculation, reducing user time and errors, and automating a significant fraction of the algorithm. The script exploits a molecular dynamics engine, which at present is NAMD, to run the short trajectories. However, since the script is external to the engine, the script can be easily adapted to different molecular dynamics codes. The outcomes of the short trajectories are analyzed to obtain a kinetic and thermodynamic description of the entire process. While many examples of Milestoning were published in the past, we provide two simple examples (a conformational transition of alanine dipeptide in a vacuum and aqueous solution) to illustrate the use of the script.
Collapse
Affiliation(s)
- Wei Wei
- Oden Institute for Computational Engineering and Sciences , University of Texas at Austin , Austin , Texas 78712 , United States
| | - Ron Elber
- Oden Institute for Computational Engineering and Sciences , University of Texas at Austin , Austin , Texas 78712 , United States
- Department of Chemistry , The University of Texas at Austin , Austin , Texas 78712 , United States
| |
Collapse
|
27
|
Abstract
Brownian dynamics (BD) is a technique for carrying out computer simulations of physical systems that are driven by thermal fluctuations. Biological systems at the macromolecular and cellular level, while falling in the gap between well-established atomic-level models and continuum models, are especially suitable for such simulations. We present a brief history, examples of important biological processes that are driven by thermal motion, and those that have been profitably studied by BD. We also present some of the challenges facing developers of algorithms and software, especially in the attempt to simulate larger systems more accurately and for longer times.
Collapse
Affiliation(s)
- Gary A Huber
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, CA 92093-0340, USA.,Department of Pharmocology, University of California San Diego, La Jolla, CA 92093-0636, USA
| | - J Andrew McCammon
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, CA 92093-0340, USA.,Department of Pharmocology, University of California San Diego, La Jolla, CA 92093-0636, USA
| |
Collapse
|
28
|
Bernetti M, Masetti M, Recanatini M, Amaro RE, Cavalli A. An Integrated Markov State Model and Path Metadynamics Approach To Characterize Drug Binding Processes. J Chem Theory Comput 2019; 15:5689-5702. [DOI: 10.1021/acs.jctc.9b00450] [Citation(s) in RCA: 32] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Affiliation(s)
- Mattia Bernetti
- Department of Pharmacy and Biotechnology, Alma Mater Studiorum, Università di Bologna, Via Belmeloro 6, I-40126 Bologna, Italy
| | - Matteo Masetti
- Department of Pharmacy and Biotechnology, Alma Mater Studiorum, Università di Bologna, Via Belmeloro 6, I-40126 Bologna, Italy
| | - Maurizio Recanatini
- Department of Pharmacy and Biotechnology, Alma Mater Studiorum, Università di Bologna, Via Belmeloro 6, I-40126 Bologna, Italy
| | - Rommie E. Amaro
- Department of Chemistry and Biochemistry, University of California San Diego, 9500 Gilman Drive, La Jolla, California 92093-0340, United States
| | - Andrea Cavalli
- Department of Pharmacy and Biotechnology, Alma Mater Studiorum, Università di Bologna, Via Belmeloro 6, I-40126 Bologna, Italy
- Computational & Chemical Biology, Istituto Italiano di Tecnologia, Via Morego 30, I-16163 Genova, Italy
| |
Collapse
|
29
|
Variational implicit-solvent predictions of the dry-wet transition pathways for ligand-receptor binding and unbinding kinetics. Proc Natl Acad Sci U S A 2019; 116:14989-14994. [PMID: 31270236 DOI: 10.1073/pnas.1902719116] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Ligand-receptor binding and unbinding are fundamental biomolecular processes and particularly essential to drug efficacy. Environmental water fluctuations, however, impact the corresponding thermodynamics and kinetics and thereby challenge theoretical descriptions. Here, we devise a holistic, implicit-solvent, multimethod approach to predict the (un)binding kinetics for a generic ligand-pocket model. We use the variational implicit-solvent model (VISM) to calculate the solute-solvent interfacial structures and the corresponding free energies, and combine the VISM with the string method to obtain the minimum energy paths and transition states between the various metastable ("dry" and "wet") hydration states. The resulting dry-wet transition rates are then used in a spatially dependent multistate continuous-time Markov chain Brownian dynamics simulation and the related Fokker-Planck equation calculations of the ligand stochastic motion, providing the mean first-passage times for binding and unbinding. We find the hydration transitions to significantly slow down the binding process, in semiquantitative agreement with existing explicit-water simulations, but significantly accelerate the unbinding process. Moreover, our methods allow the characterization of nonequilibrium hydration states of pocket and ligand during the ligand movement, for which we find substantial memory and hysteresis effects for binding vs. unbinding. Our study thus provides a significant step forward toward efficient, physics-based interpretation and predictions of the complex kinetics in realistic ligand-receptor systems.
Collapse
|
30
|
Abstract
The kinetics of drug binding and unbinding is assuming an increasingly crucial role in the long, costly process of bringing a new medicine to patients. For example, the time a drug spends in contact with its biological target is known as residence time (the inverse of the kinetic constant of the drug-target unbinding, 1/ koff). Recent reports suggest that residence time could predict drug efficacy in vivo, perhaps even more effectively than conventional thermodynamic parameters (free energy, enthalpy, entropy). There are many experimental and computational methods for predicting drug-target residence time at an early stage of drug discovery programs. Here, we review and discuss the methodological approaches to estimating drug binding kinetics and residence time. We first introduce the theoretical background of drug binding kinetics from a physicochemical standpoint. We then analyze the recent literature in the field, starting from the experimental methodologies and applications thereof and moving to theoretical and computational approaches to the kinetics of drug binding and unbinding. We acknowledge the central role of molecular dynamics and related methods, which comprise a great number of the computational methods and applications reviewed here. However, we also consider kinetic Monte Carlo. We conclude with the outlook that drug (un)binding kinetics may soon become a go/no go step in the discovery and development of new medicines.
Collapse
Affiliation(s)
- Mattia Bernetti
- Department of Pharmacy and Biotechnology, University of Bologna, I-40126 Bologna, Italy
| | - Matteo Masetti
- Department of Pharmacy and Biotechnology, University of Bologna, I-40126 Bologna, Italy
| | - Walter Rocchia
- CONCEPT Laboratory, Istituto Italiano di Tecnologia, I-16163 Genova, Italy
| | - Andrea Cavalli
- Department of Pharmacy and Biotechnology, University of Bologna, I-40126 Bologna, Italy
- Computational Sciences Domain, Istituto Italiano di Tecnologia, I-16163 Genova, Italy
| |
Collapse
|
31
|
Schuetz DA, Bernetti M, Bertazzo M, Musil D, Eggenweiler HM, Recanatini M, Masetti M, Ecker GF, Cavalli A. Predicting Residence Time and Drug Unbinding Pathway through Scaled Molecular Dynamics. J Chem Inf Model 2018; 59:535-549. [DOI: 10.1021/acs.jcim.8b00614] [Citation(s) in RCA: 44] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Doris A. Schuetz
- Department of Pharmaceutical Chemistry, University of Vienna, UZA 2, Althanstrasse 14, 1090 Vienna, Austria
| | - Mattia Bernetti
- Department of Pharmacy and Biotechnology, Alma Mater Studiorum—Università di Bologna, via Belmeloro 6, I-40126 Bologna, Italy
| | - Martina Bertazzo
- Department of Pharmacy and Biotechnology, Alma Mater Studiorum—Università di Bologna, via Belmeloro 6, I-40126 Bologna, Italy
- Computational Sciences, Istituto Italiano di Tecnologia, via Morego 30, 16163 Genova, Italy
| | - Djordje Musil
- Discovery Technologies, Merck KGaA, Frankfurter Straße 250, 64293 Darmstadt, Germany
| | | | - Maurizio Recanatini
- Department of Pharmacy and Biotechnology, Alma Mater Studiorum—Università di Bologna, via Belmeloro 6, I-40126 Bologna, Italy
| | - Matteo Masetti
- Department of Pharmacy and Biotechnology, Alma Mater Studiorum—Università di Bologna, via Belmeloro 6, I-40126 Bologna, Italy
| | - Gerhard F. Ecker
- Department of Pharmaceutical Chemistry, University of Vienna, UZA 2, Althanstrasse 14, 1090 Vienna, Austria
| | - Andrea Cavalli
- Department of Pharmacy and Biotechnology, Alma Mater Studiorum—Università di Bologna, via Belmeloro 6, I-40126 Bologna, Italy
- Computational Sciences, Istituto Italiano di Tecnologia, via Morego 30, 16163 Genova, Italy
| |
Collapse
|