Structure change of β-hairpin induced by turn optimization: An enhanced sampling molecular dynamics simulation study

Abstract

Our previous study showed that for the tested polypeptides which have similar β-hairpin structures but different sequences, their folding free energy pathways are dominantly determined by the turn conformational propensity. In this study, we study how the turn conformational propensity affects the structure of hairpins. The folding of two mutants of GB1p peptide (GB1m2 and GB1m3), which have the optimized turn sequence ((6)DDATK(11)T → (6)NPATG(11)K) with native structures unsolved, were simulated using integrated tempering sampling molecular dynamics simulations and the predicted stable structures were compared to wild-type GB1p. It was observed that the turn optimization of GB1p generates a more favored 5-residue type I(') turn in addition to the 6-residue type I turn in wild-type GB1p. As a result two distinctly different hairpin structures are formed corresponding to the "misfolded" (M) and the "folded" (F) states. M state is a one-residue-shifted asymmetric β-hairpin structure whereas F state has the similar symmetric hairpin structure as wild-type GB1p. The formation of the favored type I(') turn has a small free energy barrier and leads to the shifted β-hairpin structure, following the modified "zipping" model. The presence of disfavored type I turn structure makes the folding of a β-hairpin consistent with the "hydrophobic-core-centric" model. On the other hand, the folding simulations on other two GB1p mutants (GB1r1 and GBr2), which have the position of the hydrophobic core cluster further away from the turn compared to wild-type GB1p, showed that moving the hydrophobic core cluster away from the turn region destabilizes but does not change the hairpin structure. Therefore, the present study showed that the turn conformational propensity is a key factor in affecting not only the folding pathways but also the stable structure of β-hairpins, and the turn conformational change induced by the turn optimization leads to significant changes of β-hairpin structure.