Cloud computing approaches for prediction of ligand binding poses and pathways

We describe an innovative protocol for ab initio prediction of ligand crystallographic binding poses and highly effective analysis of large datasets generated for protein-ligand dynamics. We include a procedure for setup and performance of distributed molecular dynamics simulations on cloud computing architectures, a model for efficient analysis of simulation data, and a metric for evaluation of model convergence. We give accurate binding pose predictions for five ligands ranging in affinity from 7 nM to > 200 μM for the immunophilin protein FKBP12, for expedited results in cases where experimental structures are difficult to produce. Our approach goes beyond single, low energy ligand poses to give quantitative kinetic information that can inform protein engineering and ligand design.

[1]  Gisbert Schneider,et al.  Virtual screening: an endless staircase? , 2010, Nature Reviews Drug Discovery.

[2]  J. Mccammon,et al.  Exploring the role of receptor flexibility in structure-based drug discovery. , 2014, Biophysical chemistry.

[3]  G. de Fabritiis,et al.  Complete reconstruction of an enzyme-inhibitor binding process by molecular dynamics simulations , 2011, Proceedings of the National Academy of Sciences.

[4]  M. Parrinello,et al.  Polymorphic transitions in single crystals: A new molecular dynamics method , 1981 .

[5]  Wim F Vranken,et al.  ACPYPE - AnteChamber PYthon Parser interfacE , 2012, BMC Research Notes.

[6]  Daniel-Adriano Silva,et al.  A Role for Both Conformational Selection and Induced Fit in Ligand Binding by the LAO Protein , 2011, PLoS Comput. Biol..

[7]  Emilio Gallicchio,et al.  Conformational Transitions and Convergence of Absolute Binding Free Energy Calculations. , 2012, Journal of chemical theory and computation.

[8]  P Burkhard,et al.  The discovery of steroids and other novel FKBP inhibitors using a molecular docking program. , 1999, Journal of molecular biology.

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

[10]  Vijay S. Pande,et al.  Screen Savers of the World Unite! , 2000, Science.

[11]  Albert C. Pan,et al.  Activation mechanism of the β2-adrenergic receptor , 2011, Proceedings of the National Academy of Sciences.

[12]  Ajay N. Jain Surflex-Dock 2.1: Robust performance from ligand energetic modeling, ring flexibility, and knowledge-based search , 2007, J. Comput. Aided Mol. Des..

[13]  Yan Chen Shang,et al.  Mammalian Target of Rapamycin: Hitting the Bull's-Eye for Neurological Disorders , 2010, Oxidative medicine and cellular longevity.

[14]  J. Massagué,et al.  Mechanism of TGFβ receptor inhibition by FKBP12 case of the activin-binding protein follistatin (Nakamura , 2013 .

[15]  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.

[16]  Berk Hess,et al.  LINCS: A linear constraint solver for molecular simulations , 1997, J. Comput. Chem..

[17]  Jon Clardy,et al.  DESIGN, SYNTHESIS, AND KINETIC EVALUATION OF HIGH-AFFINITY FKBP LIGANDS AND THE X-RAY CRYSTAL-STRUCTURES OF THEIR COMPLEXES WITH FKBP12. , 1994 .

[18]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[19]  G. Bowman,et al.  Equilibrium fluctuations of a single folded protein reveal a multitude of potential cryptic allosteric sites , 2012, Proceedings of the National Academy of Sciences.

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

[21]  Diwakar Shukla,et al.  Activation pathway of Src kinase reveals intermediate states as targets for drug design , 2014 .

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

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

[24]  P. Caron,et al.  X-ray structure of calcineurin inhibited by the immunophilin-immunosuppressant FKBP12-FK506 complex , 1995, Cell.

[25]  Hongmin Li,et al.  Structural basis of conformational transitions in the active site and 80′s loop in the FK506-binding protein FKBP12 , 2014, The Biochemical journal.

[26]  R. Altman,et al.  Cloud-based simulations on Google Exacycle reveal ligand-modulation of GPCR activation pathways , 2013, Nature chemistry.

[27]  Eric Vanden-Eijnden,et al.  Transition-path theory and path-finding algorithms for the study of rare events. , 2010, Annual review of physical chemistry.

[28]  Hoover,et al.  Canonical dynamics: Equilibrium phase-space distributions. , 1985, Physical review. A, General physics.

[29]  Kazuhiko Matsuo,et al.  A peptidyl–prolyl isomerase, FKBP12, accumulates in Alzheimer neurofibrillary tangles , 2009, Neuroscience Letters.

[30]  R. Dror,et al.  Improved side-chain torsion potentials for the Amber ff99SB protein force field , 2010, Proteins.

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

[32]  Kyle A. Beauchamp,et al.  Molecular simulation of ab initio protein folding for a millisecond folder NTL9(1-39). , 2010, Journal of the American Chemical Society.

[33]  Thomas J Lane,et al.  MSMBuilder2: Modeling Conformational Dynamics at the Picosecond to Millisecond Scale. , 2011, Journal of chemical theory and computation.

[34]  Acknowledgments , 2006, Molecular and Cellular Endocrinology.

[35]  M. Karplus,et al.  A mechanism for rotamase catalysis by the FK506 binding protein (FKBP). , 1993, Biochemistry.

[36]  Carsten Kutzner,et al.  GROMACS 4:  Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation. , 2008, Journal of chemical theory and computation.

[37]  Michael R. Shirts,et al.  Direct calculation of the binding free energies of FKBP ligands. , 2005, The Journal of chemical physics.

[38]  J J Burbaum,et al.  Improved calcineurin inhibition by yeast FKBP12-drug complexes. Crystallographic and functional analysis. , 1993, The Journal of biological chemistry.

[39]  Jeffrey K Weber,et al.  Characterization and rapid sampling of protein folding Markov state model topologies. , 2011, Journal of chemical theory and computation.

[40]  Lutz Schmitt,et al.  Proteins and Their Ligands: Their Importance and How to Crystallize Them , 2013 .

[41]  Junmei Wang,et al.  Development and testing of a general amber force field , 2004, J. Comput. Chem..

[42]  Vijay S Pande,et al.  Progress and challenges in the automated construction of Markov state models for full protein systems. , 2009, The Journal of chemical physics.

[43]  Randy J. Read,et al.  Acta Crystallographica Section D Biological , 2003 .

[44]  Daniel Picot,et al.  Maltose-neopentyl glycol (MNG) amphiphiles for solubilization, stabilization and crystallization of membrane proteins , 2010, Nature Methods.

[45]  Eric Vanden-Eijnden,et al.  Transition Path Theory for Markov Jump Processes , 2009, Multiscale Model. Simul..