Premium
Application of nonlinear dimensionality reduction to characterize the conformational landscape of small peptides
Author(s) -
Stamati Hernán,
Clementi Cecilia,
Kavraki Lydia E.
Publication year - 2009
Publication title -
proteins: structure, function, and bioinformatics
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.699
H-Index - 191
eISSN - 1097-0134
pISSN - 0887-3585
DOI - 10.1002/prot.22526
Subject(s) - energy landscape , molecular dynamics , curse of dimensionality , folding (dsp implementation) , basis (linear algebra) , a priori and a posteriori , dimensionality reduction , principal component analysis , nonlinear system , metastability , diffusion map , computer science , similarity (geometry) , protein folding , statistical physics , topology (electrical circuits) , biological system , algorithm , chemistry , computational chemistry , nonlinear dimensionality reduction , physics , artificial intelligence , mathematics , geometry , thermodynamics , quantum mechanics , combinatorics , biochemistry , philosophy , image (mathematics) , engineering , biology , epistemology , electrical engineering
The automatic classification of the wealth of molecular configurations gathered in simulation in the form of a few coordinates that help to explain the main states and transitions of the system is a recurring problem in computational molecular biophysics. We use the recently proposed ScIMAP algorithm to automatically extract motion parameters from simulation data. The procedure uses only molecular shape similarity and topology information inferred directly from the simulated conformations, and is not biased by a priori known information. The automatically recovered coordinates prove as excellent reaction coordinates for the molecules studied and can be used to identify stable states and transitions, and as a basis to build free‐energy surfaces. The coordinates provide a better description of the free energy landscape when compared with coordinates computed using principal components analysis, the most popular linear dimensionality reduction technique. The method is first validated on the analysis of the dynamics of an all‐atom model of alanine dipeptide, where it successfully recover all previously known metastable states. When applied to characterize the simulated folding of a coarse‐grained model of β‐hairpin, in addition to the folded and unfolded states, two symmetric misfolding crossings of the hairpin strands are observed, together with the most likely transitions from one to the other. Proteins 2010. © 2009 Wiley‐Liss, Inc.