Free energy landscapes of peptides by enhanced conformational sampling

The free energy landscapes of peptide conformations in water have been observed by the enhanced conformational sampling method, applying the selectively enhanced multicanonical molecular dynamics simulations. The conformations of the peptide dimers, -Gly-Gly-, -Gly-Ala-, -Gly-Ser-, -Ala-Gly-, -Asn-Gly-, -Pro-Gly-, -Pro-Ala-, and -Ala-Ala-, which were all blocked with N-terminal acetyl and C-terminal N-methyl groups, were individually sampled with the explicit TIP3P water molecules. From each simulation trajectory, we obtained the canonical ensemble at 300 K, from which the individual three-dimensional landscape was drawn by the potential of mean force using the three reaction coordinates: the backbone dihedral angle, psi, of the first amino acid, the backbone dihedral angle, phi, of the second amino acid, and the distance between the carbonyl oxygen of the N-terminal acetyl group and the C-terminal amide proton. The most stable state and several meta-stable states correspond to extended conformations and typical beta-turn conformations, and their free energy values were accounted for from the potentials of mean force at the states. In addition, the contributions from the intra-molecular energies of peptides and those from the hydration effects were analyzed. Consequently, the stable beta-turn conformations in the free energy landscape were consistent with the empirically preferred beta-turn types for each amino acid sequence. The thermodynamic values for the hydration effect were decomposed and they correlated well with the empirical values estimated from the solvent accessible surface area of each molecular conformation during the trajectories. The origin of the architecture of protein local fragments was analyzed from the viewpoint of the free energy and its decomposed factors.