Modeling zymogen protein C

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.

Molecular Insights and Functional Consequences of the Interaction of Heme with Activated Protein C

Aims: In hemolysis, which is accompanied by increased levels of labile redox-active heme and is often associated with hemostatic abnormalities, a decreased activity of activated protein C (APC) is routinely detected. APC is a versatile enzyme that exerts its anticoagulant function through inactivation of clotting factors Va and VIIIa. APC has not been demonstrated to be affected by heme as described for other clotting factors and, thus, is a subject of investigation. Results: We report the interaction of heme with APC and its impact on the protein function by employing spectroscopic and physiologically relevant methods. Binding of heme to APC results in inhibition of its amidolytic and anticoagulant activity, increase of the peroxidase-like activity of heme, and protection of human umbilical vein endothelial cells from heme-induced hyperpermeability. To define the sites that are responsible for heme binding, we mapped the surface of APC for potential heme-binding motifs. T²⁸⁵GWGYHSSR²⁹³ and W³⁸⁷IHGHIRDK³⁹⁵, both located on the basic exosite, turned out as potential heme-binding sites. Molecular docking employing a homology model of full-length APC indicated Tyr²⁸⁹ and His³⁹¹ as the Fe(III)-coordinating amino acids. Innovation: The results strongly suggest that hemolysis-derived heme may directly influence the protein C pathway through binding to APC, conceivably explaining the decreased activity of APC under hemolytic conditions. Further, these results extend our understanding of heme as a multifaceted effector molecule within coagulation and may allow for an improved understanding of disease development in hemostasis under hemolytic conditions. Conclusion: Our study identifies APC as a heme-binding protein and provides insights into the functional consequences.

Zymogen and activated protein C have similar structural architecture

Activated protein C is a trypsin-like protease with anticoagulant and cytoprotective properties that is generated by thrombin from the zymogen precursor protein C in a reaction greatly accelerated by the cofactor thrombomodulin. The molecular details of this activation remain elusive due to the lack of structural information. We now fill this gap by providing information on the overall structural organization of these proteins using single molecule FRET and small angle X-ray scattering. Under physiological conditions, both zymogen and protease adopt a conformation with all domains vertically aligned along an axis 76 Å long and maximal particle size of 120 Å. This conformation is stabilized by binding of Ca2+ to the Gla domain and is affected minimally by interaction with thrombin. Hence, the zymogen protein C likely interacts with the thrombin-thrombomodulin complex through a rigid body association that produces a protease with essentially the same structural architecture. This scenario stands in contrast to an analogous reaction in the coagulation cascade where conversion of the zymogen prothrombin to the protease meizothrombin by the prothrombinase complex is linked to a large conformational transition of the entire protein. The presence of rigid epidermal growth factor domains in protein C as opposed to kringles in prothrombin likely accounts for the different conformational plasticity of the two zymogens. The new structural features reported here for protein C have general relevance to vitamin K-dependent clotting factors containing epidermal growth factor domains, such as factors VII, IX, and X.

Role of the activation peptide in the mechanism of protein C activation

Protein C is a natural anticoagulant activated by thrombin in a reaction accelerated by the cofactor thrombomodulin. The zymogen to protease conversion of protein C involves removal of a short activation peptide that, relative to the analogous sequence present in other vitamin K-dependent proteins, contains a disproportionately high number of acidic residues. Through a combination of bioinformatic, mutagenesis and kinetic approaches we demonstrate that the peculiar clustering of acidic residues increases the intrinsic disorder propensity of the activation peptide and adversely affects the rate of activation. Charge neutralization of the acidic residues in the activation peptide through Ala mutagenesis results in a mutant activated by thrombin significantly faster than wild type. Importantly, the mutant is also activated effectively by other coagulation factors, suggesting that the acidic cluster serves a protective role against unwanted proteolysis by endogenous proteases. We have also identified an important H-bond between residues T176 and Y226 that is critical to transduce the inhibitory effect of Ca2+ and the stimulatory effect of thrombomodulin on the rate of zymogen activation. These findings offer new insights on the role of the activation peptide in the function of protein C.

C-terminal residues of activated protein C light chain contribute to its anticoagulant and cytoprotective activities

BACKGROUND

Activated protein C (APC) is an important homeostatic blood coagulation protease that conveys anticoagulant and cytoprotective activities. Proteolytic inactivation of factors Va and VIIIa facilitated by cofactor protein S is responsible for APC's anticoagulant effects, whereas cytoprotective effects of APC involve primarily the endothelial protein C receptor (EPCR), protease activated receptor (PAR)1 and PAR3.

OBJECTIVE

To date, several binding exosites in the protease domain of APC have been identified that contribute to APC's interaction with its substrates but potential contributions of the C-terminus of the light chain have not been studied in detail.

METHODS

Site-directed Ala-scanning mutagenesis of six positively charged residues within G142-L155 was used to characterize their contributions to APC's anticoagulant and cytoprotective activities.

RESULTS AND CONCLUSIONS

K151 was involved in protein S dependent-anticoagulant activity of APC with some contribution of K150. 3D structural analysis supported that these two residues were exposed in an extended protein S binding site on one face of APC. Both K150 and K151 were important for PAR1 and PAR3 cleavage by APC, suggesting that this region may also mediate interactions with PARs. Accordingly, APC's cytoprotective activity as determined by endothelial barrier protection was impaired by Ala substitutions of these residues. Thus, both K150 and K151 are involved in APC's anticoagulant and cytoprotective activities. The differential contribution of K150 relative to K151 for protein S-dependent anticoagulant activity and PAR cleavage highlights that binding exosites for protein S binding and for PAR cleavage in the C-terminal region of APC's light chain overlap.

Genetic association of PROC variants with pulmonary embolism in Northern Chinese Han population

To evaluate SNPs (single nucleotide polymorphism) in PROC (protein C gene) associated with pulmonary embolism (PE) susceptibility in North Chinese Han population. A case-control study design was used, and patients with PE and healthy participants were enrolled from the Emerging Department of the several hospitals in Weifang, Shandong, China. SNPs in PROC were genotyped using Mass ARRAY system. The allele frequency of rs199469469 was significantly different between PE patients and the control [OR (95 % CI) = 5.00 (1.66-15.12), P = 0.004], and the difference remained significantly after controlling for age and gender [OR (95 % CI) = 5.34 (1.47-19.39), P = 0.011). The G(del)G in the haplotype includes rs1799809|rs199469469|rs2069928 was of a significantly difference (P = 0.016) among PE patients and the controls, and remained significant (P = 0.015) after adjustment for age and sex. Our study reports that PROC SNPs (rs199469469) might be associated with PE susceptibility, with the G allele of rs199469469 serving as the protective factors for incidence of PE. These findings may contribute to the understanding and primary prevention of PE.

The role of Ca(2+) ions in the complex assembling of protein Z and Z-dependent protease inhibitor: A structure and dynamics investigation

We investigated the solution structure and dynamics of the human anti-coagulation protein Z (PZ) in the complex with protein Zdependent protease inhibitor (ZPI) to order to understand key structural changes in the presence and absence of Ca(2+). Structural features of the complete complex of PZ-ZPI are poorly understood due to lack of complete atomic model of the PZ-ZPI complex. We have constructed a model of the complete PZ-ZPI complex and molecular dynamics (MD) simulation of the solvated PZ-ZPI complex with and without Ca(2+) was achieved for 100ns. It is consider that the Ω-loop of GLA domains interacts with negatively charged biological membranes in the presence of Ca(2+) ions. The PZ exerts its role as cofactor in a similar way. However, we used solvent-equilibrated dynamics to show structural features of the PZ-ZPI complex in the presence and the absence of Ca(2+)ions. We observed that the distance between the interacting sites of the ZPI with the PZ and the GLA domain decreases in the presence of Ca(2+) ions. Further, we postulated that the calculated distance between the dominant plane of the Ca(2+) ions and Ser196 of the pseudo-catalytic triad of the PZ is similar to the equivalent distance of FXa. This suggests that the central role of the PZ in the blood coagulation may be to align the inhibitory site of the ZPI with the active site of the FXa, which is depends on the interaction of the calcium bound GLA domain of the PZ with the active membrane.

Genetic background analysis of protein C deficiency demonstrates a recurrent mutation associated with venous thrombosis in Chinese population

Background

Protein C (PC) is one of the most important physiological inhibitors of coagulation proteases. Hereditary PC deficiency causes a predisposition to venous thrombosis (VT). The genetic characteristics of PC deficiency in the Chinese population remain unknown.

Methods

Thirty-four unrelated probands diagnosed with hereditary PC deficiency were investigated. PC activity and antigen levels were measured. Mutation analysis was performed by sequencing the PROC gene. In silico analyses, including PolyPhen-2, SIFT, multiple sequence alignment, splicing prediction, and protein molecular modeling were performed to predict the consequences of each variant identified. One recurrent mutation and its relative risk for thrombosis in relatives were analyzed in 11 families. The recurrent mutation was subsequently detected in a case (VT patients)-control study, and the adjusted odds ratio (OR) for VT risk was calculated by logistic regression analysis.

Results

A total of 18 different mutations, including 12 novel variants, were identified. One common mutation, PROC c.565C>T (rs146922325:C>T), was found in 17 of the 34 probands. The family study showed that first-degree relatives bearing this variant had an 8.8-fold (95%CI = 1.1–71.6) increased risk of venous thrombosis. The case-control (1003 vs. 1031) study identified this mutation in 5.88% patients and in 0.87% controls, respectively. The mutant allele conferred a high predisposition to venous thrombosis (adjusted OR = 7.34, 95%CI = 3.61–14.94). The plasma PC activity and antigen levels in heterozygotes were 51.73±6.92 U/dl and 75.17±4.84 U/dl, respectively.

Conclusions

This is the first study on the genetic background of PC deficiency in the Chinese population. The PROC c.565C>T mutation is the most frequent cause of PC deficiency as well as a prevalent risk factor for VT in Chinese individuals. The inclusion of this variant in routine thrombophilic detection may improve the diagnosis and prevention of venous thrombosis.

Haemostatically active proteins in snake venoms

Snake venom proteins that affect the haemostatic system can cause (a) lowering of blood coagulability, (b) damage to blood vessels, resulting in bleeding, (c) secondary effects of bleeding, e.g. hypovolaemic shock and organ damage, and (d) thrombosis. These proteins may, or may not, be enzymes. We review the data on the most relevant haemostatically active proteinases, phospholipases A₂, L-amino acid oxidases and 5'-nucleotidases from snake venoms. We also survey the non-enzymatic effectors of haemostasis from snake venoms--disintegrins, C-type lectins and three-finger toxins. Medical applications have already been found for some of these snake venom proteins. We describe those that have already been approved as drugs to treat haemostatic disorders or are being used to diagnose such health problems. No clinical applications, however, currently exist for the majority of snake venom proteins acting on haemostasis. We conclude with the most promising potential uses in this respect.

Understanding selectivity of hard and soft metal cations within biological systems using the subvalence concept. I. Application to blood coagulation: direct cation-protein electronic effects vs. indirect interactions through water networks

Following a previous study by de Courcy et al. ((2009) Interdiscip. Sci. Comput. Life Sci. 1, 55-60), we demonstrate in this contribution, using quantum chemistry, that metal cations exhibit a specific topological signature in the electron localization of their density interacting with ligands according to its "soft" or "hard" character. Introducing the concept of metal cation subvalence, we show that a metal cation can split its outer-shell density (the so-called subvalent domains or basins) according to it capability to form a partly covalent bond involving charge transfer. Such behaviour is investigated by means of several quantum chemical interpretative methods encompasing the topological analysis of the Electron Localization Function (ELF) and Bader's Quantum Theory of Atoms in Molecules (QTAIM) and two energy decomposition analyses (EDA), namely the Restricted Variational Space (RVS) and Constrained Space Orbital Variations (CSOV) approaches. Further rationalization is performed by computing ELF and QTAIM local properties such as electrostatic distributed moments and local chemical descriptors such as condensed Fukui Functions and dual descriptors. These reactivity indexes are computed within the ELF topological analysis in addition to QTAIM offering access to non atomic reactivity local index, for example on lone pairs. We apply this "subvalence" concept to study the cation selectivity in enzymes involved in blood coagulation (GLA domains of three coagulation factors). We show that the calcium ions are clearly able to form partially covalent charge transfer networks between the subdomain of the metal ion and the carboxylate oxygen lone pairs whereas magnesium does not have such ability. Our analysis also explains the different role of two groups (high affinity and low affinity cation binding sites) present in GLA domains. If the presence of Ca(II) is mandatory in the central "high affinity" region to conserve a proper folding and a charge transfer network, external sites are better stabilised by Mg(II), rather than Ca(II), in agreement with experiment. The central role of discrete water molecules is also discussed in order to understand the stabilities of the observed X-rays structures of the Gla domain. Indeed, the presence of explicit water molecules generating indirect cation-protein interactions through water networks is shown to be able to reverse the observed electronic selectivity occuring when cations directly interact with the Gla domain without the need of water.

Dopamine D1 receptor agonist and D2 receptor antagonist effects of the natural product (-)-stepholidine: molecular modeling and dynamics simulations

(-)-Stepholidine (SPD), an active ingredient of the Chinese herb Stephania, is the first compound found to have dual function as a dopamine receptor D1 agonist and D2 antagonist. Insights into dynamical behaviors of D1 and D2 receptors and their interaction modes with SPD are crucial in understanding the structural and functional characteristics of dopamine receptors. In this study a computational approach, integrating protein structure prediction, automated molecular docking, and molecular dynamics simulations were employed to investigate the dual action mechanism of SPD on the D1 and D2 receptors, with the eventual aim to develop new drugs for treating diseases affecting the central nervous system such as schizophrenia. The dynamics simulations revealed the surface features of the electrostatic potentials and the conformational "open-closed" process of the binding entrances of two dopamine receptors. Potential binding conformations of D1 and D2 receptors were obtained, and the D1-SPD and D2-SPD complexes were generated, which are in good agreement with most of experimental data. The D1-SPD structure shows that the K-167_EL-2-E-302_EL-3 (EL-2: extracellular loop 2; EL-3: extracellular loop 3) salt bridge plays an important role for both the conformational change of the extracellular domain and the binding of SPD. Based on our modeling and simulations, we proposed a mechanism of the dual action of SPD and a subsequent signal transduction model. Further mutagenesis and biophysical experiments are needed to test and improve our proposed dual action mechanism of SPD and signal transduction model.

Multifunctional specificity of the protein C/activated protein C Gla domain

Activated protein C (APC) has potent anticoagulant and anti-inflammatory properties that are mediated in part by its interactions with its cofactor protein S and the endothelial cell protein C receptor (EPCR). The protein C/APC Gla domain is implicated in both interactions. We sought to identify how the protein C Gla domain enables specific protein-protein interactions in addition to its conserved role in phospholipid binding. The human prothrombin Gla domain, which cannot bind EPCR or support protein S cofactor activity, has 22/45 residues that are not shared with the human protein C Gla domain. We hypothesized that the unique protein C/APC Gla domain residues were responsible for mediating the specific interactions. To assess this, we generated 13 recombinant protein C/APC variants incorporating the prothrombin residue substitutions. Despite anticoagulant activity similar to wild-type APC in the absence of protein S, APC variants APC(PT33-39) (N33S/V34S/D35T/D36A/L38D/A39V) and APC(PT36/38/39) (D36A/L38D/A39V) were not stimulated by protein S, whereas APC(PT35/36) (D35T/D36A) exhibited reduced protein S sensitivity. Moreover, PC(PT8/10) (L8V/H10K) displayed negligible EPCR affinity, despite normal binding to anionic phospholipid vesicles and factor Va proteolysis in the presence and absence of protein S. A single residue variant, PC(PT8), also failed to bind EPCR. Factor VIIa, which also possesses Leu-8, bound soluble EPCR with similar affinity to wild-type protein C, collectively confirming Leu-8 as the critical residue for EPCR recognition. These results reveal the specific Gla domain residues responsible for mediating protein C/APC molecular recognition with both its cofactor and receptor and further illustrate the multifunctional potential of Gla domains.

Genetic background of type I protein C deficiency in Finland

In contrast to other populations the usually rare type II form of protein C deficiency is as common in Finland as type I deficiency. We recently reported that a single mutation explained virtually all cases of type II protein C deficiency in Finland, indicating strong founder effect. We now investigated in the same population the genetic background of type I protein C deficiency. Thirty-eight apparently unrelated families were studied. They represent the vast majority of all families with type I deficiency in Finland. A genetic defect was identified in 23 (61%) families who carried 13 different mutations. Only three of the 13 mutations have been reported in other populations. Unlike in type II deficiency, considerable heterogeneity in mutations was found in type I deficiency. Our results indicate interesting differences in mutational histories of these two different forms of protein C deficiency in Finland.

The functional significance of the autolysis loop in protein C and activated protein C

The autolysis loop of activated protein C (APC) is five residues longer than the autolysis loop of other vitamin K-dependent coagulation proteases. To investigate the role of this loop in the zymogenic and anticoagulant properties of the molecule, a protein C mutant was constructed in which the autolysis loop of the protein was replaced with the corresponding loop of factor X. The protein C mutant was activated by thrombin with approximately 5-fold higher rate in the presence of Ca2+. Both kinetics and direct binding studies revealed that the Ca2+ affinity of the mutant has been impaired approximately 3-fold. The result of a factor Va degradation assay revealed that the anticoagulant function of the mutant has been improved 4-5-fold in the absence but not in the presence of protein S. The improvement was due to a better recognition of both the P1-Arg506 and P1-Arg306 cleavage sites by the mutant protease. However, the plasma half-life of the mutant was markedly shortened due to faster inactivation by plasma serpins. These results suggest that the autolysis loop of protein C is critical for the Ca(2+)-dependence of activation by thrombin. Moreover, a longer autolysis loop in APC is not optimal for interaction with factor Va in the absence of protein S, but it contributes to the lack of serpin reactivity and longer half-life of the protease in plasma.

The functional significance of the autolysis loop in protein C and activated protein C

The autolysis loop of activated protein C (APC) is five residues longer than the autolysis loop of other vitamin K-dependent coagulation proteases. To investigate the role of this loop in the zymogenic and anticoagulant properties of the molecule, a protein C mutant was constructed in which the autolysis loop of the protein was replaced with the corresponding loop of factor X. The protein C mutant was activated by thrombin with approximately 5-fold higher rate in the presence of Ca2+. Both kinetics and direct binding studies revealed that the Ca2+ affinity of the mutant has been impaired approximately 3-fold. The result of a factor Va degradation assay revealed that the anticoagulant function of the mutant has been improved 4-5-fold in the absence but not in the presence of protein S. The improvement was due to a better recognition of both the P1-Arg506 and P1-Arg306 cleavage sites by the mutant protease. However, the plasma half-life of the mutant was markedly shortened due to faster inactivation by plasma serpins. These results suggest that the autolysis loop of protein C is critical for the Ca(2+)-dependence of activation by thrombin. Moreover, a longer autolysis loop in APC is not optimal for interaction with factor Va in the absence of protein S, but it contributes to the lack of serpin reactivity and longer half-life of the protease in plasma.

Molecular docking and 3D QSAR studies on 1-amino-2-phenyl-4-(piperidin-1-yl)-butanes based on the structural modeling of human CCR5 receptor

In the present study, we have used an approach combining protein structure modeling, molecular dynamics (MD) simulation, automated docking, and 3D QSAR analyses to investigate the detailed interactions of CCR5 with their antagonists. Homology modeling and MD simulation were used to build the 3D model of CCR5 receptor based on the high-resolution X-ray structure of bovine rhodopsin. A series of 64 CCR5 antagonists, 1-amino-2-phenyl-4-(piperidin-1-yl)-butanes, were docked into the putative binding site of the 3D model of CCR5 using the docking method, and the probable interaction model between CCR5 and the antagonists were obtained. The predicted binding affinities of the antagonists to CCR5 correlate well with the antagonist activities, and the interaction model could be used to explain many mutagenesis results. All these indicate that the 3D model of antagonist-CCR5 interaction is reliable. Based on the binding conformations and their alignment inside the binding pocket of CCR5, three-dimensional structure-activity relationship (3D QSAR) analyses were performed on these antagonists using comparative molecular field analysis (CoMFA) and comparative molecular similarity analysis (CoMSIA) methods. Both CoMFA and CoMSIA provide statistically valid models with good correlation and predictive power. The q(2)(r(cross)(2)) values are 0.568 and 0.587 for CoMFA and CoMSIA, respectively. The predictive ability of these models was validated by six compounds that were not included in the training set. Mapping these models back to the topology of the active site of CCR5 leads to a better understanding of antagonist-CCR5 interaction. These results suggest that the 3D model of CCR5 can be used in structure-based drug design and the 3D QSAR models provide clear guidelines and accurate activity predictions for novel antagonist design.

The conformation of the activation peptide of protein C is influenced by Ca2+ and Na+ binding

Previous studies have suggested that the conformation of the activation peptide of protein C is influenced by the binding of Ca(2+). To provide direct evidence for the linkage between Ca(2+) binding and the conformation of the activation peptide, we have constructed a protein C mutant in the gamma-carboxyglutamic acid-domainless form in which the P1 Arg(169) of the activation peptide is replaced with the fluorescence reporter Trp. Upon binding of Ca(2+), the intrinsic fluorescence of the mutant decreases approximately 30%, as opposed to only 5% for the wild-type, indicating that Trp(169) is directly influenced by the divalent cation. The K(d) of Ca(2+) binding for the mutant protein C was impaired approximately 4-fold compared with wild-type. Interestingly, the conformation of the activation peptide was also found to be sensitive to the binding of Na(+), and the affinity for Na(+) binding increased approximately 5-fold in the presence of Ca(2+). These findings suggest that Ca(2+) changes the conformation of the activation peptide of protein C and that protein C is also capable of binding Na(+), although with a weaker affinity compared with the mature protease. The mutant protein C can no longer be activated by thrombin but remarkably it can be activated efficiently by chymotrypsin and by the thrombin mutant D189S. Activation of the mutant protein C by chymotrypsin proceeds at a rate comparable to the activation of wild-type protein C by the thrombin-thrombomodulin complex.

Characterization of a thrombomodulin binding site on protein C and its comparison to an activated protein C binding site for factor Va

Activation of the anticoagulant human plasma serine protease zymogen, protein C, by a complex of thrombin and the membrane protein, thrombomodulin, generates activated protein C, a physiologic anti-thrombotic, anti-inflammatory and anti-apoptotic agent. Alanine-scanning site-directed mutagenesis of residues in five surface loops of an extensive basic surface on protein C was used to identify residues that play essential roles in its activation by the thrombin-thrombomodulin complex. Twenty-three residues in the protein C protease domain were mutated to alanine, singly, in pairs or in triple mutation combinations, and mutants were characterized for their effectiveness as substrates of the thrombin-thrombomodulin complex. Three protein C residues, K192, R229, and R230, in two loops, were identified that provided major contributions to interactions with thrombin-thrombomodulin, while six residues, S190, K191, K217, K218, W231, and R312, in four loops, appeared to provide minor contributions. These protein C residues delineated a positively charged area on the molecule's surface that largely overlapped the previously characterized factor Va binding site on activated protein C. Thus, the extensive basic surface of protein C and activated protein C provides distinctly different, though significantly overlapping, binding sites for recognition by thrombin-thrombomodulin and factor Va.

An all-atom solution-equilibrated model for human extrinsic blood coagulation complex (sTF-VIIa-Xa): a protein-protein docking and molecular dynamics refinement study

Tissue factor (TF)-bound factor (F)VIIa plays a critical role in activating FX, an event that rapidly results in blood coagulation. Despite recent advances in the structural information about soluble TF (sTF)-bound VIIa and Xa individually, the atomic details of the ternary complex are not known. As part of our long-term goal to provide a structural understanding of the extrinsic blood coagulation pathway, we built an all atom solution-equilibrated model of the human sTF-VIIa-Xa ternary complex using protein-protein docking and molecular dynamics (MD) simulations. The starting structural coordinates of sTF-VIIa and Xa were derived from dynamically equilibrated solution structures. Due to the flexible nature of the light-chain of the Xa molecule, a three-stage docking approach was employed in which SP (Arg195-Lys448)/EGF2 (Arg86-Arg139), EGF1 (Asp46-Thr85) and GLA (Ala1-Lys45) domains were docked in a sequential manner. The rigid-body docking approach of the FTDOCK method in conjunction with filtering based on biochemical knowledge from experimental site-specific mutagenesis studies provided the strategy. The best complex obtained from the docking experiments was further refined using MD simulations for 3 ns in explicit water. In addition to explaining most of the known experimental site-specific mutagenesis data pertaining to sTF-VIIa, our model also characterizes likely enzyme-binding exosites on FVIIa and Xa that may be involved in the ternary complex formation. According to the equilibrated model, the 140s loop of VIIa serves as the key recognition motif for complex formation. Stable interactions occur between the FVIIa 140s loop and the FXa -strand B2 region near the sodium-binding domain, the 160 s loop and the N-terminal activation loop regions. The helical-hydrophobic stack region that connects the GLA and EGF1 domains of VIIa and Xa appears to play a potential role in the membrane binding region of the ternary complex. The proposed model may serve as a reasonable structural basis for understanding the exosite-mediated substrate recognition of sTF-VIIa and to advance understanding of the TFPI-mediated regulatory pathway of the extrinsic blood coagulation cascade.

A monoclonal antibody against activated protein C allows rapid detection of activated protein C in plasma and reveals a calcium ion dependent epitope involved in factor Va inactivation

Activated protein C (APC) serves as an 'on demand' anticoagulant. Defects in the APC anticoagulant pathway are underlying risk factors for the development of venous and arterial thrombosis. APC has recently been shown to significantly reduce mortality in patients with severe sepsis, presumably by virtue of its ability to down-regulate coagulation as well as inflammation. Our objective was to develop an assay that, for the first time, permits rapid detection of plasma APC. This assay will expedite studies of APC in a variety of vascular disease states including sepsis, severe atherosclerosis, diabetes, and vasculitis. By generating a highly APC-specific monoclonal antibody (HAPC 1555), we have developed an assay that, for the first time, allows rapid detection of plasma APC. The Kd measured for the interaction between APC and HAPC 1555 based on BIAcore studies and binding to immobilized HAPC on microtiter plates is 6.2 +/- 0.9 and 8.8 +/- 1.0 nmol L(-1), respectively. The interaction between HAPC 1555 and APC is Ca2+-dependent, with a Ca2+ concentration of 313 +/- 48 micro mol L(-1) required for half maximal binding. HAPC 1555 interferes with APC-mediated inactivation of factor (F)Va in the presence and absence of phospholipids, suggesting that HAPC 1555 binds to the FVa binding domain of APC. When HAPC 1555 was used in an APC enzyme capture assay, therapeutic APC levels could be measured in 1.5 h, and physiologic levels of APC could be detected between 3 and 19 h. APC levels were also shown to vary markedly in patients with severe sepsis. The rapidity of our APC assay makes APC detection in patients practical clinically. This assay will expedite studies of APC in a variety of vascular disease states including sepsis, severe atherosclerosis, diabetes, and vasculitis.

Molecular dynamics simulations on SDF-1alpha: binding with CXCR4 receptor

Insights into the interacting mode of CXCR4 with SDF-1alpha are crucial in understanding the structural and functional characteristics of CXCR4 receptor. In this paper a computational pipeline, integrating protein structure prediction, molecular dynamics simulations, automated molecular docking, and Brownian dynamics simulations were employed to investigate the dynamic and energetic aspects of CXCR4 associating with SDF-1alpha. The entire simulation revealed the surface distribution feature of electrostatic potentials and conformational "open-close" process of the receptor. The possible binding conformation of CXCR4 was identified, and the CXCR4-SDF-1alpha binding complex was generated. Arg188-Glu277 salt bridge plays an important role for both the extracellular domain conformational change and SDF-1alpha binding. Two binding sites were mapped at the extracellular domain (Site 1) and inside the transmembrane domain (Site 2), which are composed of conserved residues. Sites 1 and 2 contribute approximately 60% and 40% to the binding affinity with SDF-1alpha, respectively. The binding model is in agreement with most of the experimental data. Transmembrane VI has more significant motion in the harmonious conformational transition of CXCR4 during SDF-1alpha binding, which may be possibly associated with signal transduction. Based on the modeling and simulation, a binding mechanism hypothesis between CXCR4 and SDF-1alpha and its relationship to the signal transduction has been proposed.

Addition of missing loops and domains to protein models by x-ray solution scattering

Inherent flexibility and conformational heterogeneity in proteins can often result in the absence of loops and even entire domains in structures determined by x-ray crystallographic or NMR methods. X-ray solution scattering offers the possibility of obtaining complementary information regarding the structures of these disordered protein regions. Methods are presented for adding missing loops or domains by fixing a known structure and building the unknown regions to fit the experimental scattering data obtained from the entire particle. Simulated annealing was used to minimize a scoring function containing the discrepancy between the experimental and calculated patterns and the relevant penalty terms. In low-resolution models where interface location between known and unknown parts is not available, a gas of dummy residues represents the missing domain. In high-resolution models where the interface is known, loops or domains are represented as interconnected chains (or ensembles of residues with spring forces between the C(alpha) atoms), attached to known position(s) in the available structure. Native-like folds of missing fragments can be obtained by imposing residue-specific constraints. After validation in simulated examples, the methods have been applied to add missing loops or domains to several proteins where partial structures were available.

Structure and dynamics of zymogen human blood coagulation factor X

The solution structure and dynamics of the human coagulation factor X (FX) have been investigated to understand the key structural elements in the zymogenic form that participates in the activation process. The model was constructed based on the 2.3-A-resolution x-ray crystallographic structure of active-site inhibited human FXa (PDB:1XKA). The missing gamma-carboxyglutamic acid (GLA) and part of epidermal growth factor 1 (EGF1) domains of the light chain were modeled based on the template of GLA-EGF1 domains of the tissue factor (TF)-bound FVIIa structure (PDB:1DAN). The activation peptide and other missing segments of FX were introduced using homology modeling. The full calcium-bound model of FX was subjected to 6.2 ns of molecular dynamics simulation in aqueous medium using the AMBER6.0 package. We observed significant reorientation of the serine-protease (SP) domain upon activation leading to a compact multi-domain structure. The solution structure of zymogen appears to be in a well-extended conformation with the distance between the calcium ions in the GLA domain and the catalytic residues estimated to be approximately 95 A in contrast to approximately 83 A in the activated form. The latter is in close agreement with fluorescence studies on FXa. The S1-specificity residues near the catalytic triad show significant differences between the zymogen and activated structures.

Predicted solution structure of zymogen human coagulation FVII

A model solution structure for the complete tissue factor-free calcium ion-bound human zymogen FVII (residues 1-406) (FVII) has been constructed to study possible conformational changes associated with the activation process and tissue factor (TF) binding. The initial structure for the present model was constructed using the X-ray crystallographic structure of human coagulation FVIIa/TF complex bound with calcium ions (Banner et al., Nature 1996, 380, 41-46). This model was subsequently subjected to lengthy molecular dynamics simulations. The Amber force field in conjunction with the PME electrostatic summation method was employed. The estimated TF free solution structure was then compared with the currently available X-ray crystal structures of FVIIa (with or without TF, variable inhibitor bound) to estimate the restructuring of FVII due to TF binding and activation. The solution structure of the zymogen FVII in the absence of TF is predicted to be an extended domain structure similar to that of the TF-bound X-ray crystal structure. An additional extension of the serine protease (SP) domain of the zymogen above a reference lipid surface by approximately 7 A was in agreement with experiment. Significant Gla-EGF1 and EGF1-EGF2 interdomain motions in the zymogen were observed. Carbohydrate dimers attached to Ser-52 and Ser-60 did not cause restructuring in this domain. Minimal restructuring of the SP domain is found upon inference of the zymogen from the activated form. The catalytic triad residues maintain the H-bonded network while Lys-341 occupies the S1 specific site in the zymogen.