Abstract
Steered molecular dynamics simulation of force-induced titin immunoglobulin domain I27 unfolding led to the discovery of a significant potential energy barrier at an extension of approximately 14 A on the unfolding pathway that protects the domain against stretching. Previous simulations showed that this barrier is due to the concurrent breaking of six interstrand hydrogen bonds (H-bonds) between beta-strands A' and G that is preceded by the breaking of two to three hydrogen bonds between strands A and B, the latter leading to an unfolding intermediate. The simulation results are supported by Angstrom-resolution atomic force microscopy data. Here we perform a structural and energetic analysis of the H-bonds breaking. It is confirmed that H-bonds between strands A and B break rapidly. However, the breaking of the H-bond between strands A' and G needs to be assisted by fluctuations of water molecules. In nanosecond simulations, water molecules are found to repeatedly interact with the protein backbone atoms, weakening individual interstrand H-bonds until all six A'-G H-bonds break simultaneously under the influence of external stretching forces. Only when those bonds are broken can the generic unfolding take place, which involves hydrophobic interactions of the protein core and exerts weaker resistance against stretching than the key event.