Journal of Physical Chemistry B, Vol.114, No.8, 3028-3037, 2010
Molecular Dynamics Analysis of the Conformations of beta-Hairpin Miniprotein
Molecular dynamics simulations of a beta-hairpin miniprotein, CLN025, were performed to examine the conformational stability of the peptide in H2O at 278, 300, 333, and 363 K, as well as in TFE, MeOH, and DMSO at 300 K. CLN025 is a variant of the Chignolin miniprotein, in which the terminal Gly residues of Chignolin are replaced with Tyr residues, which leads to a 29.7 K increase in melting temperature. The energy of the intramolecular interactions was calculated using DFT quantum chemical calculations at the BHandHLYP/cc-pVTZ level of theory. CLN025 maintained a beta-hairpin conformation in all environments. The beta-hairpin is stabilized by hydrogen bonds, an electrostatic interaction between the charged termini of the peptide, and weakly polar interactions. The interaction between the backbones of the N and C-terminal strands accounts for -97.32 to -120.87 kcal mol(-1) of the stabilization energy. The energies of the CH-pi interactions between Tyr2 and Pro4 were between -1.80 and -8.9 kcal mol(-1), and the energy of the Tyr2-Trp9 Ar-Ar interaction was between -0.43 and -8.11 kcal mol(-1). Increasing temperature caused the Tyr2-Pro4 CH-pi and the Tyr2-Trp9 and Tyr-2-Tyr-10 Ar-Ar interactions to become less favorable, but the Tyr1-Trp9 interaction became more favorable and played ail important role in stabilizing the beta-hairpin of CLN025 that resulted in the increased melting temperature. Weakly polar interactions play ail important role in the structure and stability of CLN025 and other proteins.