Dynamical view of membrane binding and complex formation of human factor VIIa and tissue factor.

SUMMARY BACKGROUND: The molecular mechanism of enhancement of the enzymatic activity of factor VIIa by tissue factor (TF) is not fully understood, primarily because of the lack of atomic models for the membrane-bound form of the TF-FVIIa complex. OBJECTIVES: To construct the first membrane-bound model of the TF-FVIIa complex, and to investigate the dynamics of the complex in solution and on the surface of anionic membranes by using large-scale molecular dynamics (MD) simulations in full atomic detail. METHODS: Membrane-bound models of the TF-FVIIa complex and the individual factors were constructed and subjected to MD simulations, in order to characterize protein-protein and protein-lipid interactions, and to investigate the dynamics of TF and FVIIa. RESULTS: The MD trajectories reveal that isolated FVIIa undergoes large structural fluctuation, primarily due to the hinge motions between its domains, whereas soluble TF (sTF) is structurally stable. Upon complex formation, sTF restricts the motion of FVIIa significantly. The results also show that, in the membrane-bound form, sTF directly interacts with the lipid headgroups, even in the absence of FVIIa. CONCLUSION: The first atomic models of membrane-bound sTF-FVIIa, FVIIa and sTF are presented, revealing that sTF forms direct contacts with the lipids, both in the isolated form and in complex with FVIIa. The main effect of sTF binding to FVIIa is spatial stabilization of the catalytic site of FVIIa, which ensures optimal interaction with the substrate, FX.

Protein-membrane interactions: blood clotting on nanoscale bilayers.

The clotting cascade requires the assembly of protease-cofactor complexes on membranes with exposed anionic phospholipids. Despite their importance, protein-membrane interactions in clotting remain relatively poorly understood. Calcium ions are known to induce anionic phospholipids to cluster, and we propose that clotting proteins assemble preferentially on such anionic lipid-rich microdomains. Until recently, there was no way to control the partitioning of clotting proteins into or out of specific membrane microdomains, so experimenters only knew the average contributions of phospholipids to blood clotting. The development of nanoscale membrane bilayers (Nanodiscs) has now allowed us to probe, with nanometer resolution, how local variations in phospholipid composition regulate the activity of key protease-cofactor complexes in blood clotting. Furthermore, exciting new progress in solid-state NMR and large-scale molecular dynamics simulations allow structural insights into interactions between proteins and membrane surfaces with atomic resolution.

Potential energy function for continuous state models of globular proteins.

One of the approaches to protein structure prediction is to obtain energy functions which can recognize the native conformation of a given sequence among a zoo of conformations. The discriminations can be done by assigning the lowest energy to the native conformation, with the guarantee that the native is in the zoo. Well-adjusted functions, then, can be used in the search for other (near-) natives. Here the aim is the discrimination at relatively high resolution (RMSD difference between the native and the closest nonnative is around 1 A) by pairwise energy potentials. The potential is trained using the experimentally determined native conformation of only one protein, instead of the usual large survey over many proteins. The novel feature is that the native structure is compared to a vastly wider and more challenging array of nonnative structures found not only by the usual threading procedure, but by wide-ranging local minimization of the potential. Because of this extremely demanding search, the native is very close to the apparent global minimum of the potential function. The global minimum property holds up for one other protein having 60% sequence identity, but its performance on completely dissimilar proteins is of course much weaker.

Statistical mechanics of protein folding by exhaustive enumeration.

It is hard to construct theories for the folding of globular proteins because they are large and complicated molecules having enormous numbers of nonnative conformations and having native states that are complicated to describe. Statistical mechanical theories of protein folding are constructed around major simplifying assumptions about the energy as a function of conformation and/or simplifications of the representation of the polypeptide chain, such as one point per residue on a cubic lattice. It is not clear how the results of these theories are affected by their various simplifications. Here we take a very different simplification approach where the chain is accurately represented and the energy of each conformation is calculated by a not unreasonable empirical function. However, the set of amino acid sequences and allowed conformations is so restricted that it becomes computationally feasible to examine them all. Hence we are able to calculate melting curves for thermal denaturation as well as the detailed kinetic pathway of refolding. Such calculations are based on a novel representation of the conformations as points in an abstract 12-dimensional Euclidean conformation space. Fast folding sequences have relatively high melting temperatures, native structures with relatively low energies, small kinetic barriers between local minima, and relatively many conformations in the global energy minimum's watershed. In contrast to other folding theories, these models show no necessary relationship between fast folding and an overall funnel shape to the energy surface, or a large energy gap between the native and the lowest nonnative structure, or the depth of the native energy minimum compared to the roughness of the energy landscape.