Multiscale Estimation of Binding Kinetics Using Brownian Dynamics, Molecular Dynamics and Milestoning

The kinetic rate constants of binding were estimated for four biochemically relevant molecular systems by a method that uses milestoning theory to combine Brownian dynamics simulations with more detailed molecular dynamics simulations. The rate constants found using this method agreed well with experimentally and theoretically obtained values. We predicted the association rate of a small charged molecule toward both a charged and an uncharged spherical receptor and verified the estimated value with Smoluchowski theory. We also calculated the kon rate constant for superoxide dismutase with its natural substrate, O2 −, in a validation of a previous experiment using similar methods but with a number of important improvements. We also calculated the kon for a new system: the N-terminal domain of Troponin C with its natural substrate Ca2+. The kon calculated for the latter two systems closely resemble experimentally obtained values. This novel multiscale approach is computationally cheaper and more parallelizable when compared to other methods of similar accuracy. We anticipate that this methodology will be useful for predicting kinetic rate constants and for understanding the process of binding between a small molecule and a protein receptor.

[1]  Laxmikant V. Kalé,et al.  Scalable molecular dynamics with NAMD , 2005, J. Comput. Chem..

[2]  Eric T. Kim,et al.  How does a drug molecule find its target binding site? , 2011, Journal of the American Chemical Society.

[3]  R C Wade,et al.  Simulation of the diffusional association of barnase and barstar. , 1997, Biophysical journal.

[4]  D. Ermak,et al.  Brownian dynamics with hydrodynamic interactions , 1978 .

[5]  J. Putkey,et al.  The kinetic cycle of cardiac troponin C: Calcium binding and dissociation at site II trigger slow conformational rearrangements , 1998, Protein science : a publication of the Protein Society.

[6]  Adrian H Elcock,et al.  Molecular Simulations of Diffusion and Association in Multimacromolecular Systems , 2004, Numerical Computer Methods, Part D.

[7]  Albert C. Pan,et al.  Pathway and mechanism of drug binding to G-protein-coupled receptors , 2011, Proceedings of the National Academy of Sciences.

[8]  Li Yuan-hui,et al.  Diffusion of ions in sea water and in deep-sea sediments , 1974 .

[9]  G. Wilemski On the derivation of Smoluchowski equations with corrections in the classical theory of Brownian motion , 1976 .

[10]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[11]  Wei Zhang,et al.  A point‐charge force field for molecular mechanics simulations of proteins based on condensed‐phase quantum mechanical calculations , 2003, J. Comput. Chem..

[12]  Ron Elber,et al.  Extending molecular dynamics time scales with milestoning: example of complex kinetics in a solvated peptide. , 2007, The Journal of chemical physics.

[13]  R. Elber,et al.  Unassisted transport of N-acetyl-L-tryptophanamide through membrane: experiment and simulation of kinetics. , 2012, The journal of physical chemistry. B.

[14]  Simulation of bimolecular reactions: synthesis of the encounter and reaction steps , 1993 .

[15]  Ron Elber,et al.  Catch bond-like kinetics of helix cracking: network analysis by molecular dynamics and milestoning. , 2013, The Journal of chemical physics.

[16]  Oliver Beckstein,et al.  MDAnalysis: A toolkit for the analysis of molecular dynamics simulations , 2011, J. Comput. Chem..

[17]  Vijay S. Pande,et al.  Everything you wanted to know about Markov State Models but were afraid to ask. , 2010, Methods.

[18]  J. Mccammon,et al.  Brownian dynamics simulation of diffusion‐influenced bimolecular reactions , 1984 .

[19]  J. Berg,et al.  Molecular dynamics simulations of biomolecules , 2002, Nature Structural Biology.

[20]  William L Jorgensen,et al.  Foundations of Biomolecular Modeling , 2013, Cell.

[21]  R C Wade,et al.  Brownian dynamics simulation of protein-protein diffusional encounter. , 1998, Methods.

[22]  K Schulten,et al.  VMD: visual molecular dynamics. , 1996, Journal of molecular graphics.

[23]  F. Noé,et al.  Calculating kinetics and pathways of protein-ligand association. , 2012, European journal of cell biology.

[24]  Peter A. Kollman,et al.  AMBER, a package of computer programs for applying molecular mechanics, normal mode analysis, molecular dynamics and free energy calculations to simulate the structural and energetic properties of molecules , 1995 .

[25]  Huan‐Xiang Zhou,et al.  BDflex: a method for efficient treatment of molecular flexibility in calculating protein-ligand binding rate constants from brownian dynamics simulations. , 2012, The Journal of chemical physics.

[26]  Wilfred Pinfold,et al.  Proceedings of the Conference on High Performance Computing Networking, Storage and Analysis , 2009, HiPC 2009.

[27]  B. Sykes,et al.  Structural based insights into the role of troponin in cardiac muscle pathophysiology , 2004, Journal of Muscle Research & Cell Motility.

[28]  Ron Elber,et al.  Revisiting and computing reaction coordinates with Directional Milestoning. , 2011, The journal of physical chemistry. A.

[29]  Y Ogawa,et al.  Calcium binding to troponin C and troponin: effects of Mg2+, ionic strength and pH. , 1985, Journal of biochemistry.

[30]  Lane W. Votapka,et al.  DelEnsembleElec: Computing Ensemble-Averaged Electrostatics Using DelPhi , 2013 .

[31]  Lee G. Pedersen,et al.  Molecular Dynamics: Techniques and Applications to Proteins , 2002 .

[32]  Frank Noé,et al.  Markov models of molecular kinetics: generation and validation. , 2011, The Journal of chemical physics.

[33]  Nathan A. Baker,et al.  Electrostatics of nanosystems: Application to microtubules and the ribosome , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[34]  K. Dill,et al.  Automatic discovery of metastable states for the construction of Markov models of macromolecular conformational dynamics. , 2007, The Journal of chemical physics.

[35]  Jean-Claude Spehner,et al.  Fast and robust computation of molecular surfaces , 1995, SCG '95.

[36]  James Andrew McCammon,et al.  Browndye: A software package for Brownian dynamics , 2010, Comput. Phys. Commun..

[37]  S S Hasnain,et al.  Crystallographic structures of bovine copper-zinc superoxide dismutase reveal asymmetry in two subunits: functionally important three and five coordinate copper sites captured in the same crystal. , 1999, Journal of molecular biology.

[38]  J. Andrew McCammon,et al.  Molecular dynamics simulation of substrate-enzyme interactions in the active site channel of superoxide dismutase , 1993 .

[39]  T. J. Murphy,et al.  Brownian Motion of N Interacting Particles. I. Extension of the Einstein Diffusion Relation to the N‐Particle Case , 1972 .

[40]  Ron Elber,et al.  On the assumptions underlying milestoning. , 2008, The Journal of chemical physics.

[41]  Jan H. Jensen,et al.  Very fast prediction and rationalization of pKa values for protein–ligand complexes , 2008, Proteins.

[42]  Eric Vanden-Eijnden,et al.  Full kinetics of CO entry, internal diffusion, and exit in myoglobin from transition-path theory simulations. , 2015, Journal of the American Chemical Society.

[43]  R. Elber,et al.  Atomically detailed simulation of the recovery stroke in myosin by Milestoning , 2010, Proceedings of the National Academy of Sciences.

[44]  Long-timescale molecular dynamics simulations elucidate the dynamics and kinetics of exposure of the hydrophobic patch in troponin C. , 2012, Biophysical journal.

[45]  J. P. Grossman,et al.  Millisecond-scale molecular dynamics simulations on Anton , 2009, Proceedings of the Conference on High Performance Computing Networking, Storage and Analysis.

[46]  B. Luty,et al.  Simulation of the bimolecular reaction between superoxide and superoxide dismutase : synthesis of the encounter and reaction steps , 1993 .

[47]  John L. Klepeis,et al.  Anton, a special-purpose machine for molecular dynamics simulation , 2007, ISCA '07.

[48]  Alexander D. MacKerell,et al.  A simple polarizable model of water based on classical Drude oscillators , 2003 .

[49]  R. Elber,et al.  Computing time scales from reaction coordinates by milestoning. , 2004, The Journal of chemical physics.

[50]  B D Sykes,et al.  Calcium-induced structural transition in the regulatory domain of human cardiac troponin C. , 1997, Biochemistry.

[51]  J. Mccammon,et al.  Point charge distributions and electrostatic steering in enzyme/substrate encounter: Brownian dynamics of modified copper/zinc superoxide dismutases. , 1990, Biochemistry.

[52]  E. Argese,et al.  817—Application of short controlled drop-time polarography to the study of superoxide ion dismutation in aqueous solutions: Determination of the activity of superoxide dismutases , 1984 .

[53]  J. Andrew McCammon,et al.  Diffusional dynamics of ligand-receptor association , 1986 .

[54]  Daniel R Roe,et al.  PTRAJ and CPPTRAJ: Software for Processing and Analysis of Molecular Dynamics Trajectory Data. , 2013, Journal of chemical theory and computation.

[55]  R. Elber A milestoning study of the kinetics of an allosteric transition: atomically detailed simulations of deoxy Scapharca hemoglobin. , 2007, Biophysical journal.

[56]  F. Noé Probability distributions of molecular observables computed from Markov models. , 2008, The Journal of chemical physics.

[57]  Eric Vanden-Eijnden,et al.  Markovian milestoning with Voronoi tessellations. , 2009, The Journal of chemical physics.

[58]  R. Copeland,et al.  Drug–target residence time and its implications for lead optimization , 2007, Nature Reviews Drug Discovery.

[59]  K. Mann,et al.  Studies of the capacity of factor Xa to protect factor Va from inactivation by activated protein C. , 1982, The Journal of biological chemistry.

[60]  P. Fernandes,et al.  Molecular dynamics simulations of the enzyme Cu, Zn superoxide dismutase. , 2006, The journal of physical chemistry. B.

[61]  J. Andrew McCammon,et al.  Sodium chloride ion pair interaction in water: computer simulation , 1984 .

[62]  Christof Schütte,et al.  Markov model theory. , 2014, Advances in experimental medicine and biology.

[63]  Frank Noé,et al.  Markov state models based on milestoning. , 2011, The Journal of chemical physics.

[64]  Ron Elber,et al.  Milestoning without a Reaction Coordinate. , 2010, Journal of chemical theory and computation.

[65]  A. North,et al.  Diffusion-controlled reactions , 1966 .

[66]  J. Kestin,et al.  Viscosity of Liquid Water in the Range - 8 C to 150 C, , 1978 .

[67]  Nathan A. Baker,et al.  Finite element solution of the steady-state Smoluchowski equation for rate constant calculations. , 2004, Biophysical journal.

[68]  Dan S. Tawfik,et al.  The moderately efficient enzyme: evolutionary and physicochemical trends shaping enzyme parameters. , 2011, Biochemistry.

[69]  Stephen Jay Gould,et al.  Erratum: Peroxisomal targeting signal-1 recognition by the TPR domains of human PEX5 (Nature Structural Biology (2002) 7 (1091-1095)) , 2002 .

[70]  J. Mccammon,et al.  Rate constants for ion pair formation and dissociation in water , 1986 .

[71]  Ole Tange,et al.  GNU Parallel: The Command-Line Power Tool , 2011, login Usenix Mag..

[72]  Jeremy C. Smith,et al.  Hierarchical analysis of conformational dynamics in biomolecules: transition networks of metastable states. , 2007, The Journal of chemical physics.

[73]  I. Fridovich,et al.  Electrostatic interactions in the reaction mechanism of bovine erythrocyte superoxide dismutase. , 1982, The Journal of biological chemistry.

[74]  J. Divisek,et al.  Electrochemical generation and reactivity of the superoxide ion in aqueous solutions , 1975 .

[75]  W. C. Still,et al.  Semianalytical treatment of solvation for molecular mechanics and dynamics , 1990 .

[76]  Eric T. Kim,et al.  Oncogenic Mutations Counteract Intrinsic Disorder in the EGFR Kinase and Promote Receptor Dimerization , 2012, Cell.

[77]  S. Hardt The diffusion transit time; A simple derivation , 1981 .

[78]  P. Deuflhard,et al.  A Direct Approach to Conformational Dynamics Based on Hybrid Monte Carlo , 1999 .

[79]  Jonathan P. Davis,et al.  Designing Calcium-sensitizing Mutations in the Regulatory Domain of Cardiac Troponin C* , 2004, Journal of Biological Chemistry.

[80]  I. Fridovich,et al.  Superoxide dismutase. An enzymic function for erythrocuprein (hemocuprein). , 1969, The Journal of biological chemistry.

[81]  Elvira Guàrdia,et al.  Potential of mean force by constrained molecular dynamics: A sodium chloride ion-pair in water , 1991 .

[82]  William Swope,et al.  Describing Protein Folding Kinetics by Molecular Dynamics Simulations. 1. Theory , 2004 .

[83]  G. Hummer,et al.  Coarse master equations for peptide folding dynamics. , 2008, The journal of physical chemistry. B.