Reconstructing potential energy functions from simulated force-induced unbinding processes

Steered molecular dynamics simulation of the Rieske subunit motion in the cytochrome bc(1) complex

Crystallographic structures of the mitochondrial ubiquinol/cytochrome c oxidoreductase (cytochrome bc(1) complex) suggest that the mechanism of quinol oxidation by the bc(1) complex involves a substantial movement of the soluble head of the Rieske iron-sulfur protein (ISP) between reaction domains in cytochrome b and cytochrome c(1) subunits. In this paper we report the results of steered molecular dynamics simulations inducing, through an applied torque within 1 ns, a 56 degrees rotation of the soluble domain of ISP. For this purpose, a solvated structure of the bc(1) complex in a phospholipid bilayer (a total of 206,720 atoms) was constructed. A subset of 91,061 atoms was actually simulated with 45,131 moving atoms. Point charge distributions for the force field parametrization of heme groups and the Fe(2)S(2) cluster of the Rieske protein included in the simulated complex were determined. The simulations showed that rotation of the soluble domain of ISP is actually feasible. Several metastable conformations of the ISP during its rotation were identified and the interactions stabilizing the initial, final, and intermediate positions of the soluble head of the ISP domain were characterized. A pathway for proton conduction from the Q(o) site to the solvent via a water channel has been identified.

Forced unfolding of the fibronectin type III module reveals a tensile molecular recognition switch

The 10th type III module of fibronectin possesses a beta-sandwich structure consisting of seven beta-strands (A-G) that are arranged in two antiparallel sheets. It mediates cell adhesion to surfaces via its integrin binding motif, Arg78, Gly79, and Asp80 (RGD), which is placed at the apex of the loop connecting beta-strands F and G. Steered molecular dynamics simulations in which tension is applied to the protein's terminal ends reveal that the beta-strand G is the first to break away from the module on forced unfolding whereas the remaining fold maintains its structural integrity. The separation of strand G from the remaining fold results in a gradual shortening of the distance between the apex of the RGD-containing loop and the module surface, which potentially reduces the loop's accessibility to surface-bound integrins. The shortening is followed by a straightening of the RGD-loop from a tight beta-turn into a linear conformation, which suggests a further decrease of affinity and selectivity to integrins. The RGD-loop therefore is located strategically to undergo strong conformational changes in the early stretching stages of the module and thus constitutes a mechanosensitive control of ligand recognition.

Unbinding of retinoic acid from its receptor studied by steered molecular dynamics

Retinoic acid receptor (RAR) is a ligand-dependent transcription factor that regulates the expression of genes involved in cell growth, differentiation, and development. Binding of the retinoic acid hormone to RAR is accompanied by conformational changes in the protein which induce transactivation or transrepression of the target genes. In this paper we present a study of the hormone binding/unbinding process in order to clarify the role of some of the amino acid contacts and identify possible pathways of the all-trans retinoic acid binding/unbinding to/from human retinoic acid receptor (hRAR)-gamma. Three possible pathways were explored using steered molecular dynamics simulations. Unbinding was induced on a time scale of 1 ns by applying external forces to the hormone. The simulations suggest that the hormone may employ one pathway for binding and an alternative "back door" pathway for unbinding.

Unfolding of titin immunoglobulin domains by steered molecular dynamics simulation

Titin, a 1-microm-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 A/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 beta-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.

Binding pathway of retinal to bacterio-opsin: a prediction by molecular dynamics simulations

Formation of bacteriorhodopsin (bR) from apoprotein and retinal has been studied experimentally, but the actual pathway, including the point of entry, is little understood. Molecular dynamics simulations provide a surprisingly clear prediction. A window between bR helices E and F in the transmembrane part of the protein can be identified as an entry point for retinal. Steered molecular dynamics, performed by applying a series of external forces in the range of 200-1000 pN over a period of 0.2 ns to retinal, allows one to extract this chromophore from bR once the Schiff base bond to Lys216 is cleaved. Extraction proceeds until the retinal tail forms a hydrogen bond network with Ala144, Met145, and Ser183 side groups lining the exit/entry window. The manipulation induces a distortion with a fitted root mean square deviation of coordinates (ignoring retinal, water, and hydrogen atoms) of less than 1.9 A by the time the retinal carbonyl reaches the protein surface. The forces needed to extract retinal are due to friction and do not indicate significant potential barriers. The simulations therefore suggest a pathway for the binding of retinal. Water molecules are found to play a crucial role in the binding process.