Folding mechanism of β‐hairpins studied by replica exchange molecular simulations

Abstract

The folding process of trpzip2 beta-hairpin is studied by the replica exchange molecular dynamics (REMD) and normal MD simulations, aiming to understand the folding mechanism of this unique small, stable, and fast folder, as well as to reveal the general principles in the folding of beta-hairpins. According to our simulations, the TS ensemble is mainly characterized by a largely formed turn and the interaction between the inner pair of hydrophobic core residues. The folding is a zipping up of hydrogen bonds. However, the nascent turn has to be stabilized by the partially formed hydrophobic core to cross the TS. Thus our folding picture is in essence a blend of hydrogen bond-centric and hydrophobic core-centric mechanism. Our simulations provide a direct evidence for a very recent experiment (Du et al., Proc Natl Acad Sci USA 2004;101:15915-15920), which suggests that the turn formation is the rate-limiting step for beta-hairpin folding and the unfolding is mainly determined by the hydrophobic interactions. Besides, the relationship between hydrogen bond stabilities and their relative importance in folding are investigated. It is found that the hydrogen bonds with higher stabilities need not play more important roles in the folding process, and vice versa.