Unfolding of titin immunoglobulin domains by steered molecular dynamics simulation

Titin, a 1-mu m-long protein found in striated muscle myofibrils, possesses unique elastic and extensibility properties in its I-band region, which is largely composed of a PEVK region (70% proline, glutamic acid, valine, and lysine residue) and seven-strand beta-sandwich immunoglobulin-like (Ig) domains. The behavior of titin as a multistage entropic spring has been shown in atomic force microscope and optical tweezer experiments to partially depend on the reversible unfolding of individual Ig domains. We performed steered molecular dynamics simulations to stretch single titin Ig domains in solution with pulling speeds of 0.5 and 1.0 Angstrom/ps. Resulting force-extension profiles exhibit a single dominant peak for each Ig domain unfolding, consistent with the experimentally observed sequential, as opposed to concerted, unfolding of Ig domains under external stretching forces. This force peak can be attributed to an initial burst of backbone hydrogen bonds, which takes place between antiparallel beta-strands A and B and between parallel p-strands A' and G. Additional features of the simulations, including the position of the force peak and relative unfolding resistance of different Ig domains, can be related to experimental observations.