Free energy landscape and folding mechanism of a beta-hairpin in explicit water: a replica exchange molecular dynamics study

The free energy landscape and the folding mechanism of the C-terminal beta-hairpin of protein G is studied by extensive replica exchange molecular dynamics simulations (40 replicas and 340 ns total simulation time), using the GROMOS96 force field and the SPC explicit water solvent. The study reveals that the system preferentially adopts a beta-hairpin structure at biologically important temperatures, and that the helix content is low at all temperatures studied. Representing the free energy landscape as a function of several types of reaction coordinates, four local minima corresponding to the folded, partially folded, molten globule, and unfolded states are identified. The findings suggest that the folding of the beta-hairpin occurs as the sequence: collapse of hydrophobic core --> formation of H-bond --> formation of the turn. Identifying the folded and molten globule states as the main conformations, the free energy landscape of the beta-hairpin is consistent with a two-state behavior with a broad transition state. The temperature dependence of the folding-unfolding transition is investigated in some detail. The enthalpy and entropy jumps at the folding transition temperature are found to be about three times lower than the experimental estimates, indicating that the folding-unfolding transition in silico is less cooperative than its in vitro counterpart.