Construction of the free energy landscape of biomolecules via dihedral angle principal component analysis.

A systematic approach to construct a low-dimensional free energy landscape from a classical molecular dynamics (MD) simulation is presented. The approach is based on the recently proposed dihedral angle principal component analysis (dPCA), which avoids artifacts due to the mixing of internal and overall motions in Cartesian coordinates and circumvents problems associated with the circularity of angular variables. Requiring that the energy landscape reproduces the correct number, energy, and location of the system's metastable states and barriers, the dimensionality of the free energy landscape (i.

Dihedral angle principal component analysis of molecular dynamics simulations.

It has recently been suggested by Mu et al. [Proteins 58, 45 (2005)] to use backbone dihedral angles instead of Cartesian coordinates in a principal component analysis of molecular dynamics simulations. Dihedral angles may be advantageous because internal coordinates naturally provide a correct separation of internal and overall motion, which was found to be essential for the construction and interpretation of the free energy landscape of a biomolecule undergoing large structural rearrangements.

How complex is the dynamics of Peptide folding?

Classical molecular dynamics simulations of the folding of alanine peptides in aqueous solution are analyzed by constructing a deterministic model of the dynamics, using methods from nonlinear time series analysis. While the dimension of the free energy landscape increases with system size, a Lyapunov analysis shows that the effective dimension of the dynamic system is rather small and even decreases with chain length. The observed reduction of phase space is a nonlinear cooperative effect that is caused by intramolecular hydrogen bonds that stabilize the secondary structure of the peptides.