Membrane binding and lipid-protein interaction of the C2 domain from coagulation factor V.

Anchoring of coagulation factors to anionic regions of the membrane involves the C2 domain as a key player. The rate of enzymatic reactions of the coagulation factors is increased by several orders of magnitude upon membrane binding. However, the precise mechanisms behind the rate acceleration remain unclear, primarily because of a lack of understanding of the conformational dynamics of the C2-containing factors and corresponding complexes. We elucidate the membrane-bound form of the C2 domain from human coagulation factor V (FV-C2) by characterizing its membrane binding the specific lipid-protein interactions. Employing all-atom molecular dynamics simulations and leveraging the highly mobile membrane-mimetic (HMMM) model, we observed spontaneous binding of FV-C2 to a phosphatidylserine (PS)-containing membrane within 2-25 ns across twelve independent simulations. FV-C2 interacted with the membrane through three loops (spikes 1-3), achieving a converged, stable orientation. Multiple HMMM trajectories of the spontaneous membrane binding provided extensive sampling and ample data to examine the membrane-induced effects on the conformational dynamics of C2 as well as specific lipid-protein interactions. Despite existing crystal structures representing presumed "open" and "closed" states of FV-C2, our results revealed a continuous distribution of structures between these states, with the most populated structures differing from both "open" and "closed" states observed in crystal environments. Lastly, we characterized a putative PS-specific binding site formed by K23, Q48, and S78 located in the groove enclosed by spikes 1-3 (PS-specificity pocket), suggesting a different orientation of a bound headgroup moiety compared to previous proposals based upon analysis of static crystal structures.

Uncovering Membrane-Bound Models of Coagulation Factors by Combined Experimental and Computational Approaches.

In the life sciences, including hemostasis and thrombosis, methods of structural biology have become indispensable tools for shedding light on underlying mechanisms that govern complex biological processes. Advancements of the relatively young field of computational biology have matured to a point where it is increasingly recognized as trustworthy and useful, in part due to their high space-time resolution that is unparalleled by most experimental techniques to date. In concert with biochemical and biophysical approaches, computational studies have therefore proven time and again in recent years to be key assets in building or suggesting structural models for membrane-bound forms of coagulation factors and their supramolecular complexes on membrane surfaces where they are activated. Such endeavors and the proposed models arising from them are of fundamental importance in describing and understanding the molecular basis of hemostasis under both health and disease conditions. We summarize the body of work done in this important area of research to drive forward both experimental and computational studies toward new discoveries and potential future therapeutic strategies.

Membrane Interaction of the Factor VIIIa Discoidin Domains in Atomistic Detail.

A recently developed membrane-mimetic model was applied to study membrane interaction and binding of the two anchoring C2-like discoidin domains of human coagulation factor VIIIa (FVIIIa), the C1 and C2 domains. Both individual domains, FVIII C1 and FVIII C2, were observed to bind the phospholipid membrane by partial or full insertion of their extruding loops (the spikes). However, the two domains adopted different molecular orientations in their membrane-bound states; FVIII C2 roughly was positioned normal to the membrane plane, while FVIII C1 displayed a multitude of tilted orientations. The results indicate that FVIII C1 may be important in modulating the orientation of the FVIIIa molecule to optimize the interaction with FIXa, which is anchored to the membrane via its γ-carboxyglutamic acid-rich (Gla) domain. Additionally, a structural change was observed in FVIII C1 in the coiled main chain leading the first spike. A tight interaction with one lipid per domain, similar to what has been suggested for the homologous FVa C2, is characterized. Finally, we rationalize known FVIII antibody epitopes and the scarcity of documented hemophilic missense mutations related to improper membrane binding of FVIIIa, based on the prevalent nonspecificity of ionic interactions in the simulated membrane-bound states of FVIII C1 and FVIII C2.

Accelerating membrane insertion of peripheral proteins with a novel membrane mimetic model.

Characterizing atomic details of membrane binding of peripheral membrane proteins by molecular dynamics (MD) has been significantly hindered by the slow dynamics of membrane reorganization associated with the phenomena. To expedite lateral diffusion of lipid molecules without sacrificing the atomic details of such interactions, we have developed a novel membrane representation, to our knowledge, termed the highly mobile membrane-mimetic (HMMM) model to study binding and insertion of various molecular species into the membrane. The HMMM model takes advantage of an organic solvent layer to represent the hydrophobic core of the membrane and short-tailed phospholipids for the headgroup region. We demonstrate that using these components, bilayer structures are formed spontaneously and rapidly, regardless of the initial position and orientation of the lipids. In the HMMM membrane, lipid molecules exhibit one to two orders of magnitude enhancement in lateral diffusion. At the same time, the membrane atomic density profile of the headgroup region produced by the HMMM model is essentially identical to those obtained for full-membrane models, indicating the faithful representation of the membrane surface by the model. We demonstrate the efficiency of the model in capturing spontaneous binding and insertion of peripheral proteins by using the membrane anchor (γ-carboxyglutamic-acid-rich domain; GLA domain) of human coagulation factor VII as a test model. Achieving full insertion of the GLA domain consistently in 10 independent unbiased simulations within short simulation times clearly indicates the robustness of the HMMM model in capturing membrane association of peripheral proteins very efficiently and reproducibly. The HMMM model will provide significant improvements to the current all-atom models by accelerating lipid dynamics to examine protein-membrane interactions more efficiently.

Molecular determinants of phospholipid synergy in blood clotting.

Many regulatory processes in biology involve reversible association of proteins with membranes. Clotting proteins bind to phosphatidylserine (PS) on cell surfaces, but a clear picture of this interaction has yet to emerge. We present a novel explanation for membrane binding by GLA domains of clotting proteins, supported by biochemical studies, solid-state NMR analyses, and molecular dynamics simulations. The model invokes a single "phospho-L-serine-specific" interaction and multiple "phosphate-specific" interactions. In the latter, the phosphates in phospholipids interact with tightly bound Ca(2+) in GLA domains. We show that phospholipids with any headgroup other than choline strongly synergize with PS to enhance factor X activation. We propose that phosphatidylcholine and sphingomyelin (the major external phospholipids of healthy cells) are anticoagulant primarily because their bulky choline headgroups sterically hinder access to their phosphates. Following cell damage or activation, exposed PS and phosphatidylethanolamine collaborate to bind GLA domains by providing phospho-L-serine-specific and phosphate-specific interactions, respectively.

Capturing spontaneous partitioning of peripheral proteins using a biphasic membrane-mimetic model.

Membrane binding of peripheral proteins, mediated by specialized anchoring domains, is a crucial step for their biological function. Computational studies of membrane insertion, however, have proven challenging and largely inaccessible, due to the time scales required for the complete description of the process, mainly caused by the slow diffusion of the lipid molecules composing the membrane. Furthermore, in many cases, the nature of the membrane "anchor", i.e., the part of the protein that inserts into the membrane, is also unknown. Here, we address some of these issues by developing and employing a simplified representation of the membrane by a biphasic solvent model which we demonstrate can be used efficiently to capture and describe the process of hydrophobic insertion of membrane anchoring domains in all-atom molecular dynamics simulations. Applying the model, we have studied the insertion of the anchoring domain of a coagulation protein (the GLA domain of human protein C), starting from multiple initial configurations varying with regard to the initial orientation and height of the protein with respect to the membrane. In addition to efficiently and consistently identifying the "keel" region as the hydrophobic membrane anchor, within a few nanoseconds each configuration simulated showed a convergent height (2.20 ± 1.04 Å) and angle with respect to the interface normal (23.37 ± 12.48°). We demonstrate that the model can produce the same results as those obtained from a full representation of a membrane, in terms of both the depth of penetration and the orientation of the protein in the final membrane-bound form with an order of magnitude decrease in the required computational time compared to previous models, allowing for a more exhaustive search for the correct membrane-bound configuration.

Atomic view of calcium-induced clustering of phosphatidylserine in mixed lipid bilayers.

Membranes play key regulatory roles in biological processes, with bilayer composition exerting marked effects on binding affinities and catalytic activities of a number of membrane-associated proteins. In particular, proteins involved in diverse processes such as vesicle fusion, intracellular signaling cascades, and blood coagulation interact specifically with anionic lipids such as phosphatidylserine (PS) in the presence of Ca(2+) ions. While Ca(2+) is suspected to induce PS clustering in mixed phospholipid bilayers, the detailed structural effects of this ion on anionic lipids are not established. In this study, combining magic angle spinning (MAS) solid-state NMR (SSNMR) measurements of isotopically labeled serine headgroups in mixed lipid bilayers with molecular dynamics (MD) simulations of PS lipid bilayers in the presence of different counterions, we provide site-resolved insights into the effects of Ca(2+) on the structure and dynamics of lipid bilayers. Ca(2+)-induced conformational changes of PS in mixed bilayers are observed in both liposomes and Nanodiscs, a nanoscale membrane mimetic of bilayer patches. Site-resolved multidimensional correlation SSNMR spectra of bilayers containing (13)C,(15)N-labeled PS demonstrate that Ca(2+) ions promote two major PS headgroup conformations, which are well resolved in two-dimensional (13)C-(13)C, (15)N-(13)C, and (31)P-(13)C spectra. The results of MD simulations performed on PS lipid bilayers in the presence or absence of Ca(2+) provide an atomic view of the conformational effects underlying the observed spectra.

Protein-phospholipid interactions in blood clotting.

Most steps of the blood clotting cascade require the assembly of a serine protease with its specific regulatory protein on a suitable phospholipid bilayer. Unfortunately, the molecular details of how blood clotting proteins bind to membrane surfaces remain poorly understood, owing to a dearth of techniques for studying protein-membrane interactions at high resolution. Our laboratories are tackling this question using a combination of approaches, including nanoscale membrane bilayers, solid-state NMR, and large-scale molecular dynamics simulations. These studies are now providing structural insights at atomic resolution into clotting protein-membrane interactions.

Blood clotting reactions on nanoscale phospholipid bilayers.

Blood clotting reactions, such as those catalyzed by the tissue factor:factor VIIa complex (TF:FVIIa), assemble on membrane surfaces containing anionic phospholipids such as phosphatidylserine (PS). In fact, membrane binding is critical for the function of most of the steps in the blood clotting cascade. In spite of this, our understanding of how the membrane contributes to catalysis, or even how these proteins interact with phospholipids, is incomplete. Making matters more complicated, membranes containing mixtures of PS and neutral phospholipids are known to spontaneously form PS-rich membrane microdomains in the presence of plasma concentrations of calcium ions, and it is likely that blood-clotting proteases such as TF:FVIIa partition into these PS-rich microdomains. Unfortunately, little is known about how membrane microdomain composition influences the activity of blood-clotting proteases, which is typically not under experimental control even in "simple" model membranes. Our laboratories have developed and applied new technologies for studying membrane proteins to gain insights into how blood-clotting protease-cofactor pairs assemble and function on membrane surfaces. This includes using a novel, nanoscale bilayer system (Nanodiscs) that permits assembling blood-clotting protease-cofactor pairs on stable bilayers containing from 65 to 250 phospholipid molecules per leaflet. We have used this system to investigate how local (nanometer-scale) changes in phospholipid bilayer composition modulate TF:FVIIa activity. We have also used detailed molecular-dynamics simulations of nanoscale bilayers to provide atomic-scale predictions of how the membrane-binding domain of factor VIIa interacts with PS in membranes.

Distinct structural and adhesive roles of Ca2+ in membrane binding of blood coagulation factors.

The GLA domain, a common membrane-anchoring domain of several serine protease coagulation factors, is a key element in membrane association and activation of these factors in a highly Ca2+-dependent manner. However, the critical role of Ca2+ ions in binding is only poorly understood. Here, we present the atomic model of a membrane-bound GLA domain by using MD simulations of the GLA domain of human factor VIIa and an anionic lipid bilayer. The binding is furnished through a complete insertion of the omega-loop into the membrane and through direct interactions of structurally bound Ca2+ ions and protein side chains with negative lipids. The model suggests that Ca2+ ions play two distinct roles in the process: the four inner Ca2+ ions are primarily responsible for optimal folding of the GLA domain for membrane insertion, whereas the outer Ca2+ ions anchor the protein to the membrane through direct contacts with lipids.

Evaluating the conformational entropy of macromolecules using an energy decomposition approach.

We have developed a novel method to compute the conformational entropy of any molecular system via conventional simulation techniques. This method only requires that the total energy of the system is available and that the Hamiltonian is separable, with individual energy terms for the various degrees of freedom. Consequently the method, which we call the energy decomposition (Edcp) approach, is general and applicable to any large polymer in implicit solvent. Edcp is applied to estimate the entropy differences due to the peptide and ester groups in polyalanine and polyalanil ester. Ensembles over a wide range of temperatures were generated by replica exchange molecular dynamics, and densities of states were estimated using the weighted histogram analysis method. The results are compared with those obtained via evaluating the P ln P integral or employing the quasiharmonic approximation, other approaches widely employed to evaluate the entropy of molecular systems. Unlike the former method, Edcp can accommodate the correlations present between separate degrees of freedom. In addition, the Edcp model assumes no specific form for the underlying fluctuations present in the system, in contrast to the quasiharmonic approximation. For the molecules studied, the quasiharmonic approximation is observed to produce a good estimate of the vibrational entropy, but not of the conformational entropy. In contrast, our energy decomposition approach generates reasonable estimates for both of these entropy terms. We suggest that this approach embodies a simple yet effective solution to the problem of evaluating the conformational entropy of large macromolecules in implicit solvent.

Exploring Flory's isolated-pair hypothesis: statistical mechanics of helix-coil transitions in polyalanine and the C-peptide from RNase A.

To evaluate Flory's isolated-pair hypothesis in the context of helical peptides, we explore equilibrium conformations of alpha-helix-forming polypeptides as a function of temperature by means of replica exchange molecular dynamics in conjunction with the CHARMM/GB implicit solvent force field and the weighted histogram analysis method. From these simulations, Zimm-Bragg parameters, s and sigma, of Ac-Alan-NMe are computed as a function of temperature. The values obtained for s(T) and sigma(T) remain unchanged along the length of the polypeptide except for very short chains and yield results consistent with measurements based on short helix-forming peptides but suggest larger s values than anticipated from polymer-based measurements. From direct estimates of the density of states for Ac-Alan-NMe (n = 3-20) and peptide constructs based on the C peptide from RNase A, the conformational entropy is calculated versus temperature. The calculated S(T) shows a clear proportionality to the chain length over a wide range of temperature. This is observed in polypeptides with both significantly branched and simple methyl (alanine) side chains. These results provide evidence for the validity of Flory's isolated pair hypothesis, at least in the context of helical peptides and helix-to-coil transitions in these peptides.