A robotics-inspired method to sample conformational paths connecting known functionally-relevant structures in protein systems

Characterization of transition trajectories that take a protein between different functional states is an important yet challenging problem in computational biology. Approaches based on Molecular Dynamics can obtain the most detailed and accurate information but at considerable computational cost. To address the cost, sampling-based path planning methods adapted from robotics forego protein dynamics and seek instead conformational paths, operating under the assumption that dynamics can be incorporated later to transform paths to transition trajectories. Existing methods focus either on short peptides or large proteins; on the latter, coarse representations simplify the search space. Here we present a robotics-inspired tree-based method to sample conformational paths that connect known structural states of small- to medium- size proteins. We address the dimensionality of the search space using molecular fragment replacement to efficiently obtain physically-realistic conformations. The method grows a tree in conformational space rooted at a given conformation and biases the growth of the tree to steer it to a given goal conformation. Different bias schemes are investigated for their efficacy. Experiments on proteins up to 214 amino acids long with known functionally-relevant states more than 13ÅA apart show that the method effectively obtains conformational paths connecting significantly different structural states.

[1]  Eva Thulin,et al.  Calcium-induced structural changes and domain autonomy in calmodulin , 1995, Nature Structural Biology.

[2]  Ron Elber,et al.  Pathways of Conformational Transitions in Proteins , 2008 .

[3]  S. LaValle,et al.  Randomized Kinodynamic Planning , 2001 .

[4]  Dahlia R. Weiss,et al.  Can morphing methods predict intermediate structures? , 2009, Journal of molecular biology.

[5]  Lydia E. Kavraki,et al.  Tracing conformational changes in proteins , 2009, BIBM 2009.

[6]  Amarda Shehu,et al.  An Ab-initio tree-based exploration to enhance sampling of low-energy protein conformations , 2009, Robotics: Science and Systems.

[7]  G. Torrie,et al.  Nonphysical sampling distributions in Monte Carlo free-energy estimation: Umbrella sampling , 1977 .

[8]  Amarda Shehu,et al.  Efficient basin hopping in the protein energy surface , 2012, 2012 IEEE International Conference on Bioinformatics and Biomedicine.

[9]  Nancy M. Amato,et al.  A motion planning approach to folding: from paper craft to protein folding , 2001, Proceedings 2001 ICRA. IEEE International Conference on Robotics and Automation (Cat. No.01CH37164).

[10]  Eric Vanden-Eijnden,et al.  Transition pathways in complex systems: Application of the finite-temperature string method to the alanine dipeptide. , 2005, The Journal of chemical physics.

[11]  M. Akke,et al.  Backbone dynamics and energetics of a calmodulin domain mutant exchanging between closed and open conformations. , 1999, Journal of molecular biology.

[12]  B. Brooks,et al.  A super-linear minimization scheme for the nudged elastic band method , 2003 .

[13]  James E. Fitzgerald,et al.  Mimicking the folding pathway to improve homology-free protein structure prediction , 2009, Proceedings of the National Academy of Sciences.

[14]  Jean-Claude Latombe,et al.  Using Stochastic Roadmap Simulation to Predict Experimental Quantities in Protein Folding Kinetics: Folding Rates and Phi-Values , 2007, J. Comput. Biol..

[15]  Haruki Nakamura,et al.  Announcing the worldwide Protein Data Bank , 2003, Nature Structural Biology.

[16]  Bin W. Zhang,et al.  Efficient and verified simulation of a path ensemble for conformational change in a united-residue model of calmodulin , 2007, Proceedings of the National Academy of Sciences.

[17]  Wenjun Zheng,et al.  Identification of dynamical correlations within the myosin motor domain by the normal mode analysis of an elastic network model. , 2005, Journal of molecular biology.

[18]  Gaetano T. Montelione,et al.  3.11 News & Views 031 CDS , 2005 .

[19]  N. Mousseau,et al.  Dynamics of lennard-jones clusters: A characterization of the activation-relaxation technique , 2000 .

[20]  W. E,et al.  Finite temperature string method for the study of rare events. , 2002, Journal of Physical Chemistry B.

[21]  L. Kavraki,et al.  Multiscale characterization of protein conformational ensembles , 2009, Proteins.

[22]  Amarda Shehu,et al.  Guiding Probabilistic Search of the protein conformational Space with Structural Profiles , 2012, J. Bioinform. Comput. Biol..

[23]  J. Evanseck,et al.  Functional Significance of Hierarchical Tiers in Carbonmonoxy Myoglobin: Conformational Substates and Transitions Studied by Conformational Flooding Simulations , 2000 .

[24]  Juan Cortés,et al.  Randomized tree construction algorithm to explore energy landscapes , 2011, J. Comput. Chem..

[25]  Amarda Shehu,et al.  Guiding the Search for Native-like Protein Conformations with an Ab-initio Tree-based Exploration , 2010, Int. J. Robotics Res..

[26]  B. Faverjon,et al.  Probabilistic Roadmaps for Path Planning in High-Dimensional Con(cid:12)guration Spaces , 1996 .

[27]  TWO-WEEK Loan COpy,et al.  University of California , 1886, The American journal of dental science.

[28]  Jean-Claude Latombe,et al.  Randomized Kinodynamic Motion Planning with Moving Obstacles , 2002, Int. J. Robotics Res..

[29]  I. Bahar,et al.  Coarse-grained normal mode analysis in structural biology. , 2005, Current opinion in structural biology.

[30]  E. Vanden-Eijnden,et al.  On-the-fly string method for minimum free energy paths calculation , 2007 .

[31]  Lydia E. Kavraki,et al.  Probabilistic roadmaps for path planning in high-dimensional configuration spaces , 1996, IEEE Trans. Robotics Autom..

[32]  Wenjun Zheng,et al.  Predicting order of conformational changes during protein conformational transitions using an interpolated elastic network model , 2010, Proteins.

[33]  Pemra Doruker,et al.  Conformational transition pathways explored by Monte Carlo simulation integrated with collective modes. , 2008, Biophysical journal.

[34]  Richard Bonneau,et al.  De novo prediction of three-dimensional structures for major protein families. , 2002, Journal of molecular biology.

[35]  He Huang,et al.  A Comparison of Three Perturbation Molecular Dynamics Methods for Modeling Conformational Transitions. , 2009, Journal of chemical theory and computation.

[36]  G. Chirikjian,et al.  Iterative cluster‐NMA: A tool for generating conformational transitions in proteins , 2009, Proteins.

[37]  J. Skolnick,et al.  Local energy landscape flattening: Parallel hyperbolic Monte Carlo sampling of protein folding , 2002, Proteins.

[38]  Michael W Deem,et al.  Parallel tempering: theory, applications, and new perspectives. , 2005, Physical chemistry chemical physics : PCCP.

[39]  Elad Project,et al.  A molecular dynamics study of the effect of Ca2+ removal on calmodulin structure. , 2006, Biophysical journal.

[40]  G. Chirikjian,et al.  Efficient generation of feasible pathways for protein conformational transitions. , 2002, Biophysical journal.

[41]  Oliver Brock,et al.  Guiding conformation space search with an all‐atom energy potential , 2008, Proteins.

[42]  T. Siméon,et al.  An NMA‐guided path planning approach for computing large‐amplitude conformational changes in proteins , 2007, Proteins.

[43]  F. Ding,et al.  Ab initio folding of proteins with all-atom discrete molecular dynamics. , 2008, Structure.

[44]  Emilio Gallicchio,et al.  Conformational equilibria and free energy profiles for the allosteric transition of the ribose-binding protein. , 2005, Journal of molecular biology.

[45]  Amarda Shehu,et al.  In Search of the protein Native State with a Probabilistic Sampling Approach , 2011, J. Bioinform. Comput. Biol..

[46]  Eric Vanden-Eijnden,et al.  Simplified and improved string method for computing the minimum energy paths in barrier-crossing events. , 2007, The Journal of chemical physics.

[47]  Lydia E. Kavraki,et al.  Motion Planning in the Presence of Drift, Underactuation and Discrete System Changes , 2005, Robotics: Science and Systems.

[48]  R. Swendsen,et al.  THE weighted histogram analysis method for free‐energy calculations on biomolecules. I. The method , 1992 .

[49]  Wayne A Hendrickson,et al.  Computation of conformational transitions in proteins by virtual atom molecular mechanics as validated in application to adenylate kinase , 2009, Proceedings of the National Academy of Sciences.

[50]  Wilfred F. van Gunsteren,et al.  Biomolecular Modeling: Goals, Problems, Perspectives , 2006 .

[51]  P. Bradley,et al.  Toward High-Resolution de Novo Structure Prediction for Small Proteins , 2005, Science.

[52]  P. Wolynes,et al.  Water in protein structure prediction. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[53]  C. Brooks,et al.  Large-scale allosteric conformational transitions of adenylate kinase appear to involve a population-shift mechanism , 2007, Proceedings of the National Academy of Sciences.

[54]  Thierry Siméon,et al.  A path planning approach for computing large-amplitude motions of flexible molecules , 2005, ISMB.

[55]  Howie Choset,et al.  Principles of Robot Motion: Theory, Algorithms, and Implementation ERRATA!!!! 1 , 2007 .