Characterizing Energy Landscapes of Peptides Using a Combination of Stochastic Algorithms

Obtaining accurate representations of energy landscapes of biomolecules such as proteins and peptides is central to the study of their physicochemical properties and biological functions. Peptides are particularly interesting, as they exploit structural flexibility to modulate their biological function. Despite their small size, peptide modeling remains challenging due to the complexity of the energy landscape of such highly-flexible dynamic systems. Currently, only stochastic sampling-based methods can efficiently explore the conformational space of a peptide. In this paper, we suggest to combine two such methods to obtain a full characterization of energy landscapes of small yet flexible peptides. First, we propose a simplified version of the classical Basin Hopping algorithm to reveal low-energy regions in the landscape, and thus to identify the corresponding meta-stable structural states of a peptide. Then, we present several variants of a robotics-inspired algorithm, the Transition-based Rapidly-exploring Random Tree, to quickly determine transition path ensembles, as well as transition probabilities between meta-stable states. We demonstrate this combined approach on met-enkephalin.

[1]  A. Pohorille,et al.  The development/application of a ‘minimalist’ organic/biochemical molecular mechanic force field using a combination of ab initio calculations and experimental data , 1997 .

[2]  Robert D Skeel,et al.  Maximum Flux Transition Paths of Conformational Change. , 2009, Journal of chemical theory and computation.

[3]  Julius Jellinek,et al.  Energy Landscapes: With Applications to Clusters, Biomolecules and Glasses , 2005 .

[4]  Amarda Shehu,et al.  Evolutionary-inspired probabilistic search for enhancing sampling of local minima in the protein energy surface , 2012, Proteome Science.

[5]  Steven M. LaValle,et al.  Rapidly-Exploring Random Trees: Progress and Prospects , 2000 .

[6]  Didier Devaurs,et al.  Parallelizing RRT on Large-Scale Distributed-Memory Architectures , 2013, IEEE Transactions on Robotics.

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

[8]  Rahul Banerjee,et al.  Transition paths of Met-enkephalin from Markov state modeling of a molecular dynamics trajectory. , 2014, The journal of physical chemistry. B.

[9]  H A Scheraga,et al.  Prediction of the native conformation of a polypeptide by a statistical‐mechanical procedure. III. Probable and average conformations of enkephalin , 1987, Biopolymers.

[10]  P. Senet,et al.  Reconstructing the free-energy landscape of Met-enkephalin using dihedral principal component analysis and well-tempered metadynamics. , 2012, The Journal of chemical physics.

[11]  Seung-Hoon Choi,et al.  Density-Based Clustering of Small Peptide Conformations Sampled from a Molecular Dynamics Simulation , 2009, J. Chem. Inf. Model..

[12]  Woo Youn Kim,et al.  Efficient Basin-Hopping Sampling of Reaction Intermediates through Molecular Fragmentation and Graph Theory. , 2014, Journal of chemical theory and computation.

[13]  D Picone,et al.  Conformational sampling of bioactive conformers: a low‐temperature NMR study of 15N‐Leu–enkephalin , 1998, Journal of peptide science : an official publication of the European Peptide Society.

[14]  J. Straub,et al.  The MaxFlux algorithm for calculating variationally optimized reaction paths for conformational transitions in many body systems at finite temperature , 1997 .

[15]  Lydia E Kavraki,et al.  Computational models of protein kinematics and dynamics: beyond simulation. , 2012, Annual review of analytical chemistry.

[16]  Kentaro Shimizu,et al.  Folding free‐energy landscape of a 10‐residue mini‐protein, chignolin , 2006, FEBS letters.

[17]  Thierry Siméon,et al.  Sampling-Based Path Planning on Configuration-Space Costmaps , 2010, IEEE Transactions on Robotics.

[18]  J. Doye,et al.  Global Optimization by Basin-Hopping and the Lowest Energy Structures of Lennard-Jones Clusters Containing up to 110 Atoms , 1997, cond-mat/9803344.

[19]  J. Straub,et al.  Direct computation of long time processes in peptides and proteins: Reaction path study of the coil‐to‐helix transition in polyalanine , 1999, Proteins.

[20]  Berend Smit,et al.  Understanding molecular simulation: from algorithms to applications , 1996 .

[21]  Didier Devaurs,et al.  A Multi-Tree Approach to Compute Transition Paths on Energy Landscapes , 2013, AAAI 2013.

[22]  Berend Smit,et al.  Understanding Molecular Simulations: from Algorithms to Applications , 2002 .

[23]  Nirmal Singh,et al.  Extending pharmacological spectrum of opioids beyond analgesia: Multifunctional aspects in different pathophysiological states , 2011, Neuropeptides.

[24]  Didier Devaurs,et al.  Sampling-based methods for a full characterization of energy landscapes of small peptides , 2014, 2014 IEEE International Conference on Bioinformatics and Biomedicine (BIBM).

[25]  Thierry Siméon,et al.  Motion planning algorithms for molecular simulations: A survey , 2012, Comput. Sci. Rev..