Molecular dynamics simulations of cis-trans isomerization for a proline-containing tripeptide in solution

The cis-trans isomerization of X-Pro (X is Pro or any other residue) often plays a rate-limiting role in protein folding. In order to study the dynamic properties of X-Pro in water, we have analyzed the trajectory data obtained from constrained molecular dynamics (MD) simulations of a tripeptide Phe-Pro-Ala in solution. Via calculation of autocorrelation functions and their corresponding spectral densities, we have found that the structural fluctuations of the solute are affected by the motions of various degrees of freedom, and the solvent dynamic behavior is bimodal rather than exponential or Gaussian. In addition, according to the energy and torque analysis, we suggest that both the anti/endo and the syn/endo could be the possible transition state for the cis-trans isomerization of the tripeptide.