Forced unfolding of fibronectin type 3 modules: An analysis by biased molecular dynamics simulations

Titin, an important constituent of vertebrate muscles, is a protein of the order of a micrometer in length in the folded state. Atomic force microscopy and laser tweezer experiments have been used to stretch titin molecules to more than ten times their folded lengths. To explain the observed relation between force and extension, it has been suggested that the immunoglobulin and fibronectin domains unfold one at a time in an all-or-none fashion. We use molecular dynamics simulations to study the forced unfolding of two different fibronectin type 3 domains (the ninth, 9Fn3, and the tenth, 10Fn3, from human fibronectin) and of their heterodimer of known structure. An external biasing potential on the N to C distance is employed and the protein is treated in the polar hydrogen representation with an implicit solvation model. The latter provides an adiabatic solvent response, which is important for the nanosecond unfolding simulation method used here. A series of simulations is performed for each system to obtain meaningful results. The two different fibronectin domains are shown to unfold in the same way along two possible pathways. These involve the partial separation of the 'P-sandwich', an essential structural element, and the unfolding of the individual sheets in a stepwise fashion. The biasing potential results are confirmed by constant force unfolding simulations. For the two connected domains, there is complete unfolding of one domain (9Fn3) before major unfolding of the second domain (10Fn3). Comparison of different models for the potential energy function demonstrates that the dominant cohesive element in both proteins is due to the attractive van der Waals interactions; electrostatic interactions play a structural role but appear to make only a small contribution to the stabilization of the domains, in agreement with other studies of β-sheet stability. The unfolding forces found in the simulations are of the order of those observed experimentally, even though the speed of the former is more than six orders of magnitude greater than that used in the latter.