The molecular basis of thrombin allostery revealed by a 1.8 A structure of the "slow" form

Simulations suggest double sodium binding induces unexpected conformational changes in thrombin

Thrombin is a Na[Formula: see text]-activated serine protease existing in two forms targeted to procoagulant and anticoagulant activities, respectively. There is one Na[Formula: see text]-binding site that has been the focus of the study of the thrombin. However, molecular dynamics (MD) simulations suggest that there might be actually two Na[Formula: see text]-binding sites in thrombin and that Na[Formula: see text] ions can even bind to two sites simultaneously. In this study, we performed 12 independent 2-µs all-atom MD simulations for the wild-type (WT) thrombin and we studied the effects of the different Na[Formula: see text] binding modes on thrombin. From the root-mean-square fluctuations (RMSF) for the [Formula: see text]-carbons, we see that the atomic fluctuations mainly change in the 60s, 170s, and 220s loops, and the connection (residue 167 to 170). The correlation matrices for different binding modes suggest regions that may play an important role in thrombin's allosteric response and provide us a possible allosteric pathway for the sodium binding. Amorim-Hennig (AH) clustering tells us how the structure of the regions of interest changes on sodium binding. Principal component analysis (PCA) shows us how the different regions of thrombin change conformation together with sodium binding. Solvent-accessible surface area (SASA) exposes the conformational change in exosite I and catalytic triad. Finally, we argue that the double binding mode might be an inactive mode and that the kinetic scheme for the Na[Formula: see text] binding to thrombin might be a multiple-step mechanism rather than a 2-step mechanism.

Antithrombin Resistance Rescues Clotting Defect of Homozygous Prothrombin-Y510N Dysprothrombinemia

A patient with hematuria in our clinic was diagnosed with urolithiasis. Analysis of the patient's plasma clotting time indicated that both activated partial thromboplastin time (52.6 seconds) and prothrombin time (19.4 seconds) are prolonged and prothrombin activity is reduced to 12.4% of normal, though the patient exhibited no abnormal bleeding phenotype and a prothrombin antigen level of 87.9%. Genetic analysis revealed the patient is homozygous for prothrombin Y510N mutation. We expressed and characterized the prothrombin-Y510N variant in appropriate coagulation assays and found that the specificity constant for activation of the mutant zymogen by factor Xa is impaired approximately fivefold. Thrombin generation assay using patient's plasma and prothrombin-deficient plasma supplemented with either wild-type or prothrombin-Y510N revealed that both peak height and time to peak for the prothrombin mutant are decreased; however, the endogenous thrombin generation potential is increased. Further analysis indicated that the thrombin mutant exhibits resistance to antithrombin and is inhibited by the serpin with approximately 12-fold slower rate constant. Protein C activation by thrombin-Y510N was also decreased by approximately 10-fold; however, thrombomodulin overcame the catalytic defect. The Na⁺-concentration-dependence of the amidolytic activities revealed that the dissociation constant for the interaction of Na⁺ with the mutant has been elevated approximately 20-fold. These results suggest that Y510 (Y184a in chymotrypsin numbering) belongs to network of residues involved in binding Na⁺. A normal protein C activation by thrombin-Y510N suggests that thrombomodulin modulates the conformation of the Na⁺-binding loop of thrombin. The clotting defect of thrombin-Y510N appears to be compensated by its markedly lower reactivity with antithrombin, explaining patient's normal hemostatic phenotype.

Proton Bridging in Catalysis by and Inhibition of Serine Proteases of the Blood Cascade System

Inquiries into the participation of short hydrogen bonds in stabilizing transition states and intermediate states in the thrombin, factor Xa, plasmin and activated protein C-catalyzed reactions revealed that specific binding of effectors at Sn, n = 1-4 and S'n, n = 1-3 and at remote exosites elicit complex patterns of hydrogen bonding and involve water networks. The methods employed that yielded these discoveries include; (1) kinetics, especially partial or full kinetic deuterium solvent isotope effects with short cognate substrates and also with the natural substrates, (2) kinetic and structural probes, particularly low-field high-resolution nuclear magnetic resonance (1H NMR), of mechanism-based inhibitors and substrate-mimic peptide inhibitors. Short hydrogen bonds form at the transition states of the catalytic reactions at the active site of the enzymes as they do with mechanism-based covalent inhibitors of thrombin. The emergence of short hydrogen bonds at the binding interface of effectors and thrombin at remote exosites has recently gained recognition. Herein, I describe our contribution, a confirmation of this discovery, by low-field 1H NMR. The principal conclusion of this review is that proton sharing at distances below the sum of van der Waals radii of the hydrogen and both donor and acceptor atoms contribute to the remarkable catalytic prowess of serine proteases of the blood clotting system and other enzymes that employ acid-base catalysis. Proton bridges also play a role in tight binding in proteins and at exosites, i.e., allosteric sites, of enzymes.

Conformational Plasticity-Rigidity Axis of the Coagulation Factor VII Zymogen Elucidated by Atomistic Simulations of the N-Terminally Truncated Factor VIIa Protease Domain

The vast majority of coagulation factor VII (FVII), a trypsin-like protease, circulates as the inactive zymogen. Activated FVII (FVIIa) is formed upon proteolytic activation of FVII, where it remains in a zymogen-like state and it is fully activated only when bound to tissue factor (TF). The catalytic domains of trypsin-like proteases adopt strikingly similar structures in their fully active forms. However, the dynamics and structures of the available corresponding zymogens reveal remarkable conformational plasticity of the protease domain prior to activation in many cases. Exactly how ligands and cofactors modulate the conformational dynamics and function of these proteases is not entirely understood. Here, we employ atomistic simulations of FVIIa (and variants hereof, including a TF-independent variant and N-terminally truncated variants) to provide fundamental insights with atomistic resolution into the plasticity-rigidity interplay of the protease domain conformations that appears to govern the functional response to proteolytic and allosteric activation. We argue that these findings are relevant to the FVII zymogen, whose structure has remained elusive despite substantial efforts. Our results shed light on the nature of FVII and demonstrate how conformational dynamics has played a crucial role in the evolutionary adaptation of regulatory mechanisms that were not present in the ancestral trypsin. Exploiting this knowledge could lead to engineering of protease variants for use as next-generation hemostatic therapeutics.

Co-evolution networks of HIV/HCV are modular with direct association to structure and function

Mutational correlation patterns found in population-level sequence data for the Human Immunodeficiency Virus (HIV) and the Hepatitis C Virus (HCV) have been demonstrated to be informative of viral fitness. Such patterns can be seen as footprints of the intrinsic functional constraints placed on viral evolution under diverse selective pressures. Here, considering multiple HIV and HCV proteins, we demonstrate that these mutational correlations encode a modular co-evolutionary structure that is tightly linked to the structural and functional properties of the respective proteins. Specifically, by introducing a robust statistical method based on sparse principal component analysis, we identify near-disjoint sets of collectively-correlated residues (sectors) having mostly a one-to-one association to largely distinct structural or functional domains. This suggests that the distinct phenotypic properties of HIV/HCV proteins often give rise to quasi-independent modes of evolution, with each mode involving a sparse and localized network of mutational interactions. Moreover, individual inferred sectors of HIV are shown to carry immunological significance, providing insight for guiding targeted vaccine strategies.

Electrostatic recognition in substrate binding to serine proteases

Serine proteases of the Chymotrypsin family are structurally very similar but have very different substrate preferences. This study investigates a set of 9 different proteases of this family comprising proteases that prefer substrates containing positively charged amino acids, negatively charged amino acids, and uncharged amino acids with varying degree of specificity. Here, we show that differences in electrostatic substrate preferences can be predicted reliably by electrostatic molecular interaction fields employing customized GRID probes. Thus, we are able to directly link protease structures to their electrostatic substrate preferences. Additionally, we present a new metric that measures similarities in substrate preferences focusing only on electrostatics. It efficiently compares these electrostatic substrate preferences between different proteases. This new metric can be interpreted as the electrostatic part of our previously developed substrate similarity metric. Consequently, we suggest, that substrate recognition in terms of electrostatics and shape complementarity are rather orthogonal aspects of substrate recognition. This is in line with a 2‐step mechanism of protein‐protein recognition suggested in the literature.

Dynamic Consequences of Mutation of Tryptophan 215 in Thrombin

Thrombin normally cleaves fibrinogen to promote coagulation; however, binding of thrombomodulin to thrombin switches the specificity of thrombin toward protein C, triggering the anticoagulation pathway. The W215A thrombin mutant was reported to have decreased activity toward fibrinogen without significant loss of activity toward protein C. To understand how mutation of Trp215 may alter thrombin specificity, hydrogen-deuterium exchange experiments (HDXMS), accelerated molecular dynamics (AMD) simulations, and activity assays were carried out to compare the dynamics of Trp215 mutants with those of wild type (WT) thrombin. Variation in NaCl concentration had no detectable effect on the sodium-binding (220s_CT) loop, but appeared to affect other surface loops. Trp215 mutants showed significant increases in amide exchange in the 170s_CT loop consistent with a loss of H-bonding in this loop identified by the AMD simulations. The W215A thrombin showed increased amide exchange in the 220s_CT loop and in the N-terminus of the heavy chain. The AMD simulations showed that a transient conformation of the W215A thrombin has a distorted catalytic triad. HDXMS experiments revealed that mutation of Phe227, which engages in a π-stacking interaction with Trp215, also caused significantly increased amide exchange in the 170s_CT loop. Activity assays showed that only the F227V mutant had wild type catalytic activity, whereas all other mutants showed markedly lower activity. Taken together, the results explain the reduced pro-coagulant activity of the W215A mutant and demonstrate the allosteric connection between Trp215, the sodium-binding loop, and the active site.

Visualizing correlated motion with HDBSCAN clustering

Correlated motion analysis provides a method for understanding communication between and dynamic similarities of biopolymer residues and domains. The typical equal-time correlation matrices-frequently visualized with pseudo-colorings or heat maps-quickly convey large regions of highly correlated motion but hide more subtle similarities of motion. Here we propose a complementary method for visualizing correlations within proteins (or general biopolymers) that quickly conveys intuition about which residues have a similar dynamic behavior. For grouping residues, we use the recently developed non-parametric clustering algorithm HDBSCAN. Although the method we propose here can be used to group residues using correlation as a similarity matrix-the most straightforward and intuitive method-it can also be used to more generally determine groups of residues which have similar dynamic properties. We term these latter groups "Dynamic Domains", as they are based not on spatial closeness but rather closeness in the column space of a correlation matrix. We provide examples of this method across three human proteins of varying size and function-the Nf-Kappa-Beta essential modulator, the clotting promoter Thrombin and the mismatch repair protein (dimer) complex MutS-alpha. Although the examples presented here are from all-atom molecular dynamics simulations, this visualization technique can also be used on correlations matrices built from any ensembles of conformations from experiment or computation.

Uncovering Large-Scale Conformational Change in Molecular Dynamics without Prior Knowledge

As the length of molecular dynamics (MD) trajectories grows with increasing computational power, so does the importance of clustering methods for partitioning trajectories into conformational bins. Of the methods available, the vast majority require users to either have some a priori knowledge about the system to be clustered or to tune clustering parameters through trial and error. Here we present non-parametric uses of two modern clustering techniques suitable for first-pass investigation of an MD trajectory. Being non-parametric, these methods require neither prior knowledge nor parameter tuning. The first method, HDBSCAN, is fast-relative to other popular clustering methods-and is able to group unstructured or intrinsically disordered systems (such as intrinsically disordered proteins, or IDPs) into bins that represent global conformational shifts. HDBSCAN is also useful for determining the overall stability of a system-as it tends to group stable systems into one or two bins-and identifying transition events between metastable states. The second method, iMWK-Means, with explicit rescaling followed by K-Means, while slower than HDBSCAN, performs well with stable, structured systems such as folded proteins and is able to identify higher resolution details such as changes in relative position of secondary structural elements. Used in conjunction, these clustering methods allow a user to discern quickly and without prior knowledge the stability of a simulated system and identify both local and global conformational changes.

Thrombomodulin enhances complement regulation through strong affinity interactions with factor H and C3b-Factor H complex

INTRODUCTION

Coagulation and complement systems are simultaneously activated at sites of tissue injury, leading to thrombin generation and opsonisation with C3b. Thrombomodulin (TM) is a cell-bound regulator of thrombin activation, but can also enhance the regulatory activity of complement factor H (FH), thus accelerating the degradation of C3b into inactive iC3b.

OBJECTIVES

This study sought to determine the biophysical interaction affinities of two recombinant TM analogs with thrombin, FH and C3b in order to analyze their ability to regulate serum complement activity.

METHODS

Surface plasmon resonance (SPR) analysis was used to determine binding affinities of TM analogs with FH and C3b, and compared to thrombin as positive control. The capacity of the two recombinant TM analogs to regulate complement in serum was tested in standard complement hemolytic activity assays.

RESULTS

SPR analysis showed that both TM analogs bind FH and C3b-Factor H with nanomolar and C3b with micromolar affinity; binding affinity for its natural ligand thrombin was several fold higher than for FH. At a physiological relevant concentration, TM inhibits complement hemolytic activity in serum via FH dependent and independent mechanisms.

CONCLUSIONS

TM exhibits significant binding affinity for complement protein FH and C3b-FH complex and its soluble form is capable at physiologically relevant concentrations of inhibiting complement activation in serum.

Releasing the brakes in coagulation Factor IXa by co-operative maturation of the substrate-binding site

Coagulation Factor IX is positioned at the merging point of the intrinsic and extrinsic blood coagulation cascades. Factor IXa (activated Factor IX) serves as the trigger for amplification of coagulation through formation of the so-called Xase complex, which is a ternary complex of Factor IXa, its substrate Factor X and the cofactor Factor VIIIa on the surface of activated platelets. Within the Xase complex the substrate turnover by Factor IXa is enhanced 200000-fold; however, the mechanistic and structural basis for this dramatic enhancement remains only partly understood. A multifaceted approach using enzymatic, biophysical and crystallographic methods to evaluate a key set of activity-enhanced Factor IXa variants has demonstrated a delicately balanced bidirectional network. Essential molecular interactions across multiple regions of the Factor IXa molecule co-operate in the maturation of the active site. This maturation is specifically facilitated by long-range communication through the Ile212–Ile213 motif unique to Factor IXa and a flexibility of the 170-loop that is further dependent on the conformation in the Cys168–Cys182 disulfide bond. Ultimately, the network consists of compensatory brakes (Val16 and Ile213) and accelerators (Tyr99 and Phe174) that together allow for a subtle fine-tuning of enzymatic activity.

Molecular Basis of Enhanced Activity in Factor VIIa-Trypsin Variants Conveys Insights into Tissue Factor-mediated Allosteric Regulation of Factor VIIa Activity

The complex of coagulation factor VIIa (FVIIa), a trypsin-like serine protease, and membrane-bound tissue factor (TF) initiates blood coagulation upon vascular injury. Binding of TF to FVIIa promotes allosteric conformational changes in the FVIIa protease domain and improves its catalytic properties. Extensive studies have revealed two putative pathways for this allosteric communication. Here we provide further details of this allosteric communication by investigating FVIIa loop swap variants containing the 170 loop of trypsin that display TF-independent enhanced activity. Using x-ray crystallography, we show that the introduced 170 loop from trypsin directly interacts with the FVIIa active site, stabilizing segment 215-217 and activation loop 3, leading to enhanced activity. Molecular dynamics simulations and novel fluorescence quenching studies support that segment 215-217 conformation is pivotal to the enhanced activity of the FVIIa variants. We speculate that the allosteric regulation of FVIIa activity by TF binding follows a similar path in conjunction with protease domain N terminus insertion, suggesting a more complete molecular basis of TF-mediated allosteric enhancement of FVIIa activity.

Geometrical Frustration in Interleukin-33 Decouples the Dynamics of the Functional Element from the Folding Transition State Ensemble

Interleukin-33 (IL-33) is currently the focus of multiple investigations into targeting pernicious inflammatory disorders. This mediator of inflammation plays a prevalent role in chronic disorders such as asthma, rheumatoid arthritis, and progressive heart disease. In order to better understand the possible link between the folding free energy landscape and functional regions in IL-33, a combined experimental and theoretical approach was applied. IL-33 is a pseudo- symmetrical protein composed of three distinct structural elements that complicate the folding mechanism due to competition for nucleation on the dominant folding route. Trefoil 1 constitutes the majority of the binding interface with the receptor whereas Trefoils 2 and 3 provide the stable scaffold to anchor Trefoil 1. We identified that IL-33 folds with a three-state mechanism, leading to a rollover in the refolding arm of its chevron plots in strongly native conditions. In addition, there is a second slower refolding phase that exhibits the same rollover suggesting similar limitations in folding along parallel routes. Characterization of the intermediate state and the rate limiting steps required for folding suggests that the rollover is attributable to a moving transition state, shifting from a post- to pre-intermediate transition state as you move from strongly native conditions to the midpoint of the transition. On a structural level, we found that initially, all independent Trefoil units fold equally well until a Q_CA of 0.35 when Trefoil 1 will backtrack in order to allow Trefoils 2 and 3 to fold in the intermediate state, creating a stable scaffold for Trefoil 1 to fold onto during the final folding transition. The formation of this intermediate state and subsequent moving transition state is a result of balancing the difficulty in folding the functionally important Trefoil 1 onto the remainder of the protein. Taken together our results indicate that the functional element of the protein is geometrically frustrated, requiring the more stable elements to fold first, acting as a scaffold for docking of the functional element to allow productive folding to the native state.

Thrombomodulin Binding Selects the Catalytically Active Form of Thrombin

Human α-thrombin is a serine protease with dual functions. Thrombin acts as a procoagulant, cleaving fibrinogen to make the fibrin clot, but when bound to thrombomodulin (TM), it acts as an anticoagulant, cleaving protein C. A minimal TM fragment consisting of the fourth, fifth, and most of the sixth EGF-like domain (TM456m) that has been prepared has much improved solubility, thrombin binding capacity, and anticoagulant activity versus those of previous TM456 constructs. In this work, we compare backbone amide exchange of human α-thrombin in three states: apo, D-Phe-Pro-Arg-chloromethylketone (PPACK)-bound, and TM456m-bound. Beyond causing a decreased level of amide exchange at their binding sites, TM and PPACK both cause a decreased level of amide exchange in other regions including the γ-loop and the adjacent N-terminus of the heavy chain. The decreased level of amide exchange in the N-terminus of the heavy chain is consistent with the historic model of activation of serine proteases, which involves insertion of this region into the β-barrel promoting the correct conformation of the catalytic residues. Contrary to crystal structures of thrombin, hydrogen-deuterium exchange mass spectrometry results suggest that the conformation of apo-thrombin does not yet have the N-terminus of the heavy chain properly inserted for optimal catalytic activity, and that binding of TM allosterically promotes the catalytically active conformation.

Toward understanding allosteric activation of thrombin: a conjecture for important roles of unbound Na(+) molecules around thrombin

We shed light on important roles of unbound Na(+) molecules in enzymatic activation of thrombin. Molecular mechanism of Na(+)-activation of thrombin has been discussed in the context of allostery. However, the recent challenge to redesign K(+)-activated thrombin revealed that the allosteric interaction is insufficient to explain the mechanism. Under these circumstances, we have examined the roles of unbound Na(+) molecule in maximization of thrombin-substrate association reaction rate. We performed all-atomic molecular dynamics (MD) simulations of thrombin in the presence of three different cations; Li(+), Na(+), and Cs(+). Although these cations are commonly observed in the vicinity of the S1-pocket of thrombin, smaller cations are distributed more densely and extensively than larger ones. This suggests the two observation rules: (i) thrombin surrounded by Na(+) is at an advantage in the initial step of association reaction, namely, the formation of an encounter complex ensemble, and (ii) the presence of Na(+) molecules does not necessarily have an advantage in the final step of association reaction, namely, the formation of the stereospecific complex. In conclusion, we propose a conjecture that unbound Na(+) molecules also affect the maximization of rate constant of thrombin-substrate association reaction through optimally forming an encounter complex ensemble.

Ensemble refinement shows conformational flexibility in crystal structures of human complement factor D

Human factor D (FD) is a self-inhibited thrombin-like serine proteinase that is critical for amplification of the complement immune response. FD is activated by its substrate through interactions outside the active site. The substrate-binding, or `exosite', region displays a well defined and rigid conformation in FD. In contrast, remarkable flexibility is observed in thrombin and related proteinases, in which Na(+) and ligand binding is implied in allosteric regulation of enzymatic activity through protein dynamics. Here, ensemble refinement (ER) of FD and thrombin crystal structures is used to evaluate structure and dynamics simultaneously. A comparison with previously published NMR data for thrombin supports the ER analysis. The R202A FD variant has enhanced activity towards artificial peptides and simultaneously displays active and inactive conformations of the active site. ER revealed pronounced disorder in the exosite loops for this FD variant, reminiscent of thrombin in the absence of the stabilizing Na(+) ion. These data indicate that FD exhibits conformational dynamics like thrombin, but unlike in thrombin a mechanism has evolved in FD that locks the unbound native state into an ordered inactive conformation via the self-inhibitory loop. Thus, ensemble refinement of X-ray crystal structures may represent an approach alternative to spectroscopy to explore protein dynamics in atomic detail.

GpIbα interacts exclusively with exosite II of thrombin

Activation of platelets by the serine protease thrombin is a critical event in haemostasis. This process involves the binding of thrombin to glycoprotein Ibα (GpIbα) and cleavage of protease-activated receptors (PARs). The N-terminal extracellular domain of GpIbα contains an acidic peptide stretch that has been identified as the main thrombin binding site, and both anion binding exosites of thrombin have been implicated in GpIbα binding, but it remains unclear how they are involved. This issue is of critical importance for the mechanism of platelet activation by thrombin. If both exosites bind to GpIbα, thrombin could potentially act as a platelet adhesion molecule or receptor dimerisation trigger. Alternatively, if only a single site is involved, GpIbα may serve as a cofactor for PAR-1 activation by thrombin. To determine the involvement of thrombin's two exosites in GpIbα binding, we employed the complementary methods of mutational analysis, binding studies, X-ray crystallography and NMR spectroscopy. Our results indicate that the peptide corresponding to the C-terminal portion of GpIbα and the entire extracellular domain bind exclusively to thrombin's exosite II. The interaction of thrombin with GpIbα thus serves to recruit thrombin activity to the platelet surface while leaving exosite I free for PAR-1 recognition.

•

We analysed interactions of the platelet receptor GpIbα with thrombin using three complementary methods.

•

GpIbα exclusively binds to exosite II of thrombin.

•

Exosite I remains available for binding to other ligands.

•

GpIbα recruits thrombin to the platelet membrane as a cofactor for PAR-1 cleavage.

Bi-factor analysis based on noise-reduction (BIFANR): a new algorithm for detecting coevolving amino acid sites in proteins

Previous statistical analyses have shown that amino acid sites in a protein evolve in a correlated way instead of independently. Even though located distantly in the linear sequence, the coevolved amino acids could be spatially adjacent in the tertiary structure, and constitute specific protein sectors. Moreover, these protein sectors are independent of one another in structure, function, and even evolution. Thus, systematic studies on protein sectors inside a protein will contribute to the clarification of protein function. In this paper, we propose a new algorithm BIFANR (Bi-factor Analysis Based on Noise-reduction) for detecting protein sectors in amino acid sequences. After applying BIFANR on S1A family and PDZ family, we carried out internal correlation test, statistical independence test, evolutionary rate analysis, evolutionary independence analysis, and function analysis to assess the prediction. The results showed that the amino acids in certain predicted protein sector are closely correlated in structure, function, and evolution, while protein sectors are nearly statistically independent. The results also indicated that the protein sectors have distinct evolutionary directions. In addition, compared with other algorithms, BIFANR has higher accuracy and robustness under the influence of noise sites.

Thrombin inhibition by the serpins

Thrombin is the central protease in the blood coagulation network. It has multiple substrates and cofactors, and it appears that four serpins are responsible for inhibiting the thrombin produced in haemostasis and thrombosis. Structural studies conducted over the last 10 years have resolved how thrombin recognises these serpins with the aid of cofactors. Although antithrombin (AT), protein C inhibitor (PCI), heparin cofactor II (HCII) and protease nexin-1 (PN1) all share a common fold and mechanism of protease inhibition, they have evolved radically different mechanisms for cofactor-assisted thrombin recognition. This is likely to be due to the varied environments in which thrombin is found. In this review, I discuss the unusual structural features of thrombin that are involved in substrate and cofactor recognition, the serpin mechanism of protease inhibition and the fate of thrombin in the complex, and how the four thrombin-specific serpins exploit the special features of thrombin to accelerate complex formation.

Structure and behavior of human α-thrombin upon ligand recognition: thermodynamic and molecular dynamics studies

Thrombin is a serine proteinase that plays a fundamental role in coagulation. In this study, we address the effects of ligand site recognition by alpha-thrombin on conformation and energetics in solution. Active site occupation induces large changes in secondary structure content in thrombin as shown by circular dichroism. Thrombin-D-Phe-Pro-Arg-chloromethyl ketone (PPACK) exhibits enhanced equilibrium and kinetic stability compared to free thrombin, whose difference is rooted in the unfolding step. Small-angle X-ray scattering (SAXS) measurements in solution reveal an overall similarity in the molecular envelope of thrombin and thrombin-PPACK, which differs from the crystal structure of thrombin. Molecular dynamics simulations performed with thrombin lead to different conformations than the one observed in the crystal structure. These data shed light on the diversity of thrombin conformers not previously observed in crystal structures with distinguished catalytic and conformational behaviors, which might have direct implications on novel strategies to design direct thrombin inhibitors.

Looking at the proteases from a simple perspective

Proteases have received enormous interest from the research and medical communities because of their significant roles in several human diseases. Some examples include the involvement of thrombin in thrombosis, HIV-1 protease in Acquired Immune Deficiency Syndrome, cruzain in Trypanosoma cruzi infection, and membrane-type 1 matrix metalloproteinase in tumor invasion and metastasis. Many efforts has been undertaken to design effective inhibitors featuring potent inhibitory activity, specificity, and metabolic stability to those proteases involved in such pathologies. Protease inhibitors usually target the active site, but some of them act by other inhibitory mechanisms. The understanding of the structure-function relationships of proteases and inhibitors has an impact on new inhibitor drugs designing. In this paper, the structures of four proteases (thrombin, HIV-protease, cruzain, and a matrix metalloproteinase) are briefly reviewed, and used as examples of the importance of proteases for the development of new treatment strategies, leading to a longer and healthier life.

Molecular basis of thrombomodulin activation of slow thrombin

BACKGROUND

Coagulation is a highly regulated process where the ability to prevent blood loss after injury is balanced against the maintenance of blood fluidity. Thrombin is at the center of this balancing act. It is the critical enzyme for producing and stabilizing a clot, but when complexed with thrombomodulin (TM) it is converted to a powerful anticoagulant. Another cofactor that may play a role in determining thrombin function is the monovalent cation Na(+). Its apparent affinity suggests that half of the thrombin generated is in a Na(+)-free 'slow' state and half is in a Na(+)-coordinated 'fast' state. While slow thrombin is a poor procoagulant enzyme, when complexed to TM it is an effective anticoagulant.

METHODS

To better understand this molecular transformation we solved a 2.4 A structure of thrombin complexed with EGF domains 4-6 of TM in the absence of Na(+) and other cofactors or inhibitors.

RESULTS

We find that TM binds as previously observed, and that the thrombin component resembles structures of the fast form. The Na(+) binding loop is observed in a conformation identical to the Na(+)-bound form, with conserved water molecules compensating for the missing ion. Using the fluorescent probe p-aminobenzamidine we show that activation of slow thrombin by TM principally involves the opening of the primary specificity pocket.

CONCLUSIONS

These data show that TM binding alters the conformation of thrombin in a similar manner as Na(+) coordination, resulting in an ordering of the Na(+) binding loop and an opening of the adjacent S1 pocket. We conclude that other, more subtle subsite changes are unlikely to influence thrombin specificity toward macromolecular substrates.

Protein sectors: evolutionary units of three-dimensional structure

Proteins display a hierarchy of structural features at primary, secondary, tertiary, and higher-order levels, an organization that guides our current understanding of their biological properties and evolutionary origins. Here, we reveal a structural organization distinct from this traditional hierarchy by statistical analysis of correlated evolution between amino acids. Applied to the S1A serine proteases, the analysis indicates a decomposition of the protein into three quasi-independent groups of correlated amino acids that we term "protein sectors." Each sector is physically connected in the tertiary structure, has a distinct functional role, and constitutes an independent mode of sequence divergence in the protein family. Functionally relevant sectors are evident in other protein families as well, suggesting that they may be general features of proteins. We propose that sectors represent a structural organization of proteins that reflects their evolutionary histories.

How Na+ activates thrombin--a review of the functional and structural data

Thrombin is the ultimate coagulation factor; it is the final protease generated in the blood coagulation cascade and is the effector of clot formation. Regulation of thrombin activity is thus of great relevance to determining the correct haemostatic balance, with dysregulation leading to bleeding or thrombosis. One of the most enigmatic and controversial regulators of thrombin activity is the monovalent cation Na+. When bound to Na+, thrombin adopts a 'fast' conformation which cleaves all procoagulant substrates more rapidly, and when free of Na+, thrombin reverts to a 'slow' state which preferentially activates the protein C anticoagulant pathway. Thus, Na+-binding allosterically modulates the activity of thrombin and helps determine the haemostatic balance. Over the last 30 years, there has been much research investigating the structural basis of thrombin allostery. Biochemical and mutagenesis studies established which regions and residues are involved in the slow-->fast conformational change, and recently several crystal structures of the putative slow form have been solved. In this article, the biochemical and crystallographic data are reviewed to see if we are any closer to understanding the conformational basis of the Na+ activation of thrombin.

Role of A-chain in functioning of the active site of human alpha-thrombin

This review summarizes current data suggesting that A-chain of the human alpha-thrombin molecule plays a role of allosteric effector in catalytic reactions with various substrates. Special attention is paid to the relationship between A-chain structure and catalytic activity of thrombin. The existence of this relationship is based on studies of natural mutation of A-chain of the alpha-thrombin molecule. Use of molecular and essential dynamics confirmed the role of A-chain in changes of conformation and catalytic properties of this enzyme; these changes involve residues located in the specificity sites and some inserting loops. Current knowledge on structure and properties of thrombin can be used for the development of new antithrombin agents.

Molecular basis of thrombin recognition by protein C inhibitor revealed by the 1.6-A structure of the heparin-bridged complex

Protein C inhibitor (PCI) is a serpin with many roles in biology, including a dual role as pro- and anticoagulant in blood. The protease specificity and local function of PCI depend on its interaction with cofactors such as heparin-like glycosaminoglycans (GAGs) and thrombomodulin (TM). Both cofactors significantly increase the rate of thrombin inhibition, but GAGs serve to promote the anticoagulant activity of PCI, and TM promotes its procoagulant function. To gain insight into how PCI recognition of thrombin is aided by these cofactors, we determined a crystallographic structure of the Michaelis complex of PCI, thrombin, and heparin to 1.6 A resolution. Thrombin interacts with PCI in an unusual fashion that depends on the length of PCI's reactive center loop (RCL) to align the heparin-binding sites of the two proteins. The principal exosite contact is engendered by movement of thrombin's 60-loop in response to the unique P2 Phe of PCI. This mechanism of communication between the active site of thrombin and its recognition exosite is previously uncharacterized and may relate to other thrombin substrate-cofactor interactions. The cofactor activity of heparin thus depends on the formation of a heparin-bridged Michaelis complex and substrate-induced exosite contacts. We also investigated the cofactor effect of TM, establishing that TM bridges PCI to thrombin through additional direct interactions. A model of the PCI-thrombin-TM complex was built and evaluated by mutagenesis and suggests distinct binding sites for heparin and TM on PCI. These data significantly improve our understanding of the cofactor-dependent roles of PCI in hemostasis.

Thrombin

Thrombin is a Na+-activated, allosteric serine protease that plays opposing functional roles in blood coagulation. Binding of Na+ is the major driving force behind the procoagulant, prothrombotic and signaling functions of the enzyme, but is dispensable for cleavage of the anticoagulant protein C. The anticoagulant function of thrombin is under the allosteric control of the cofactor thrombomodulin. Much has been learned on the mechanism of Na+ binding and recognition of natural substrates by thrombin. Recent structural advances have shed light on the remarkable molecular plasticity of this enzyme and the molecular underpinnings of thrombin allostery mediated by binding to exosite I and the Na+ site. This review summarizes our current understanding of the molecular basis of thrombin function and allosteric regulation. The basic information emerging from recent structural, mutagenesis and kinetic investigation of this important enzyme is that thrombin exists in three forms, E*, E and E:Na+, that interconvert under the influence of ligand binding to distinct domains. The transition between the Na+ -free slow from E and the Na+ -bound fast form E:Na+ involves the structure of the enzyme as a whole, and so does the interconversion between the two Na+ -free forms E* and E. E* is most likely an inactive form of thrombin, unable to interact with Na + and substrate. The complexity of thrombin function and regulation has gained this enzyme pre-eminence as the prototypic allosteric serine protease. Thrombin is now looked upon as a model system for the quantitative analysis of biologically important enzymes.

Exosites in the substrate specificity of blood coagulation reactions

The specificity of blood coagulation proteinases for substrate, inhibitor, and effector recognition is mediated by exosites on the surfaces of the catalytic domains, physically separated from the catalytic site. Some thrombin ligands bind specifically to either exosite I or II, while others engage both exosites. The involvement of different, overlapping constellations of exosite residues enables binding of structurally diverse ligands. The flexibility of the thrombin structure is central to the mechanism of complex formation and the specificity of exosite interactions. Encounter complex formation is driven by electrostatic ligand-exosite interactions, followed by conformational rearrangement to a stable complex. Exosites on some zymogens are in low affinity proexosite states and are expressed concomitant with catalytic site activation. The requirement for exosite expression controls the specificity of assembly of catalytic complexes on the coagulation pathway, such as the membrane-bound factor Xa*factor Va (prothrombinase) complex, and prevents premature assembly. Substrate recognition by prothrombinase involves a two-step mechanism with initial docking of prothrombin to exosites, followed by a conformational change to engage the FXa catalytic site. Prothrombin and its activation intermediates bind prothrombinase in two alternative conformations determined by the zymogen to proteinase transition that are hypothesized to involve prothrombin (pro)exosite I interactions with FVa, which underpin the sequential activation pathway. The role of exosites as the major source of substrate specificity has stimulated development of exosite-targeted anticoagulants for treatment of thrombosis.

Essential thrombin residues for inhibition by protein C inhibitor with the cofactors heparin and thrombomodulin

BACKGROUND

Protein C inhibitor (PCI) and antithrombin (AT) are serine protease inhibitors (serpins) that inhibit a wide array of blood coagulation serine proteases including thrombin.

OBJECTIVE

Fifty-five Ala-scanned recombinant thrombin mutants were used to determine thrombin residues important for inhibition by PCI with and without the cofactors heparin and thrombomodulin (TM) and compared with the prototypical serpin, AT.

RESULTS

Residues around the active site (Tyr50 and Glu202) and the sodium-binding site (Glu229 and Arg233) were required for thrombin inhibition by PCI with and without cofactors. Exosite-2 residues (Arg89, Arg93, Glu94, Arg98, Arg245, Arg248, and Gln251) were critical for heparin-accelerated inhibition of thrombin by PCI. Exosite-1 residues (especially Lys65 and Tyr71) were required for enhanced PCI inhibition of thrombin-TM. Interestingly, we also found that the TM chondroitin sulfate moiety is not required for the approximately 150-fold enhanced rate of thrombin inhibition by PCI. Using the aforementioned thrombin exosite-2 mutants that were essential for heparin-catalyzed PCI-thrombin inhibition reactions we found no change in PCI inhibition rates for thrombin-TM.

CONCLUSIONS

Collectively, these results show that (i) similar thrombin exosite-2 residues are critical for the heparin-catalyzed inhibition by PCI and AT, (ii) PCI and AT are different in their thrombin-TM inhibition properties, and (iii) PCI has a distinct advantage over AT in the regulation of the activity of thrombin-TM.

The central role of thrombin in hemostasis

Following vascular injury, blood loss is controlled by the mechanisms of hemostasis. During this process, the serine proteinase, thrombin, is generated both locally and rapidly at sites of vessel damage. It plays a pivotal role in clot promotion and inhibition, and cell signaling, as well as additional processes that influence fibrinolysis and inflammation. These functions involve numerous cleavage reactions, which must be tightly coordinated. Failure to do so can lead to either bleeding or thrombosis. The crystal structures of thrombin, in combination with biochemical analyses of thrombin mutants, have provided insight into the ways in which thrombin functions, and how its different activities are modulated. Many of the interactions of thrombin are facilitated by exosites on its surface that bind to its substrates and/or cofactors. The use of cofactors not only extends the range of thrombin specificity, but also enhances its catalytic efficiency for different substrates. This explains a paradox (i.e. thrombin is a specific proteinase, and yet one that has multiple, and sometimes opposing, substrate reactions). In this review, we describe the context in which thrombin acts during hemostasis and explain the roles that its exosites and cofactors play in directing thrombin function. Thereafter, we develop the concept of cofactor competition as a means by which the activities of thrombin are controlled.

Solvent accessibility of protein surfaces by amide H/2H exchange MALDI-TOF mass spectrometry

One advantage of detecting amide H/2H exchange by mass spectrometry instead of NMR is that the more rapidly exchanging surface amides are still detectable. In this study, we present quench-flow amide H/2H exchange experiments to probe how rapidly the surfaces of two different proteins exchange. We compared the amide H/2H exchange behavior of thrombin, a globular protein, and IkappaBalpha, a nonglobular protein, to explore any differences in the determinants of amide H/2H exchange rates for each class of protein. The rates of exchange of only a few of the surface amides were as rapid as the "intrinsic" exchange rates measured for amides in unstructured peptides. Most of the surface amides exchanged at a slower rate, despite the fact that they were not seen to be hydrogen bonded to another protein group in the crystal structure. To elucidate the influence of the surface environment on amide H/2H exchange, we compared exchange data with the number of amides participating in hydrogen bonds with other protein groups and with the solvent accessible surface area. The best correlation with amide H/2H exchange was found with the total solvent accessible surface area, including side chains. In the case of the globular protein, the correlation was modest, whereas it was well correlated for the nonglobular protein. The nonglobular protein also showed a correlation between amide exchange and hydrogen bonding. These data suggest that other factors, such as complex dynamic behavior and surface burial, may alter the expected exchange rates in globular proteins more than in nonglobular proteins where all of the residues are near the surface.

Crystal structure of monomeric native antithrombin reveals a novel reactive center loop conformation

The poor inhibitory activity of circulating antithrombin (AT) is critical to the formation of blood clots at sites of vascular damage. AT becomes an efficient inhibitor of the coagulation proteases only after binding to a specific heparin pentasaccharide, which alters the conformation of the reactive center loop (RCL). The molecular basis of this activation event lies at the heart of the regulation of hemostasis and accounts for the anticoagulant properties of the low molecular weight heparins. Although several structures of AT have been solved, the conformation of the RCL in native AT remains unknown because of the obligate crystal contact between the RCL of native AT and its latent counterpart. Here we report the crystallographic structure of a variant of AT in its monomeric native state. The RCL shifted approximately 20 A, and a salt bridge was observed between the P1 residue (Arg-393) and Glu-237. This contact explains the effect of mutations at the P1 position on the affinity of AT for heparin and also the properties of AT-Truro (E237K). The relevance of the observed conformation was verified through mutagenesis studies and by solving structures of the same variant in different crystal forms. We conclude that the poor inhibitory activity of the circulating form of AT is partially conferred by intramolecular contacts that restrain the RCL, orient the P1 residue away from attacking proteases, and additionally block the exosite utilized in protease recognition.

Crystal structure of thrombin in a self-inhibited conformation

The activating effect of Na(+) on thrombin is allosteric and depends on the conformational transition from a low activity Na(+)-free (slow) form to a high activity Na(+)-bound (fast) form. The structures of these active forms have been solved. Recent structures of thrombin obtained in the absence of Na(+) have also documented inactive conformations that presumably exist in equilibrium with the active slow form. The validity of these inactive slow form structures, however, is called into question by the presence of packing interactions involving the Na(+) site and the active site regions. Here, we report a 1.87A resolution structure of thrombin in the absence of inhibitors and salts with a single molecule in the asymmetric unit and devoid of significant packing interactions in regions involved in the allosteric slow --> fast transition. The structure shows an unprecedented self-inhibited conformation where Trp-215 and Arg-221a relocate >10A to occlude the active site and the primary specificity pocket, and the guanidinium group of Arg-187 penetrates the protein core to fill the empty Na(+)-binding site. The extreme mobility of Trp-215 was investigated further with the W215P mutation. Remarkably, the mutation significantly compromises cleavage of the anticoagulant protein C but has no effect on the hydrolysis of fibrinogen and PAR1. These findings demonstrate that thrombin may assume an inactive conformation in the absence of Na(+) and that its procoagulant and anticoagulant activities are closely linked to the mobility of residue 215.

Thrombin-Cofactor Interactions

Precise modulation of thrombin activity throughout the hemostatic response is essential for efficient cessation of bleeding while preventing inappropriate clot growth or dissemination which causes thrombosis. Regulating thrombin activity is made difficult by its ability to diffuse from the surface on which it was generated and its ability to cleave at least 12 substrates. To overcome this challenge, thrombin recognition of substrates is largely controlled by cofactors that act by localizing thrombin to various surfaces, blocking substrate binding to critical exosites, engendering new exosites for substrate recognition and by allosterically modulating the properties of the active site of thrombin. Thrombin cofactors can be classified as either pro- or anticoagulants, depending on how substrate preference is altered. The procoagulant cofactors include glycoprotein Ibα, fibrin, and Na + , and the anticoagulants are heparin and thrombomodulin. Over the last few years, crystal structures have been reported for all of the thrombin-cofactor complexes. The purpose of this article is to summarize the features of these structures and to discuss the mechanisms and physiological relevance of cofactor binding in thrombin regulation.

High resolution crystal structures of free thrombin in the presence of K+ reveal the molecular basis of monovalent cation selectivity and an inactive slow form

Structural biology has recently advanced our understanding of the molecular mechanisms of activation and selectivity in monovalent cation activated enzymes. Here we report a 1.9 Angstrom resolution crystal structure of free thrombin, a Na(+) selective enzyme, in the presence of KCl. There are two molecules in the asymmetric unit, one with the cation site bound to K(+) and the other with this site free. The K(+)-bound form shows key differences compared with the Na(+)-bound structure that explain the different kinetics of activation. The cation-free form, on the other hand, assumes a conformation where the monovalent cation binding site is completely disordered, the S1 pocket is inaccessible to substrate and binding to exosite I is compromised by an unprecedented >20 Angstrom shift in the position of the autolysis loop. This form, named S(*), corresponds to the inactive Na(+)-free slow form identified by early kinetic studies. A simple model of thrombin allostery that incorporates the contribution of S(*) is proposed.

Prothrombin Saint-Denis: a natural variant with a point mutation resulting in Asp to Glu substitution at position 552 in prothrombin

A new prothrombin variant, with a point mutation at nucleotide 20 029 resulting in Asp 552 to Glu substitution (prothrombin numbering), has been identified in a male newborn. Plasma prothrombin level was <3%, 16% and 60% when measured by clotting, chromogenic and immunological assays respectively. The substitution did not affect the rate of prothrombin conversion to thrombin but altered thrombin activity. Amino acid 552 has been reported to be involved in the allosteric transition, which is induced by sodium binding to thrombin. This is the first known amino acid substitution at this site to result in dysprothrombinaemia.

Structure and interaction modes of thrombin

Any vascular injury triggers the burst-like release of the trypsin-like serine proteinase alpha-thrombin. Thrombin, the main executioner of the coagulation cascade, exhibits procoagulant as well as anticoagulant and antifibrinolytic properties, very specifically interacting with a number of protein substrates, receptors, cofactors, inhibitors, carbohydrates, and modulators. A large number of crystal structures of alpha-thrombin have shown that the thrombin surface can be subdivided into several functional regions, which recognize different substrates, inhibitors, and mediators with high specificity.

Crystal structure of wild-type human thrombin in the Na+-free state

Regulation of thrombin activity is critical for haemostasis and the prevention of thrombosis. Thrombin has several procoagulant substrates, including fibrinogen and platelet receptors, and essential cofactors for stimulating its own formation. However, thrombin is also capable of serving an anticoagulant function by activating protein C. The specificity of thrombin is primarily regulated by binding to the cofactor TM (thrombomodulin), but co-ordination of Na+ can also affect thrombin activity. The Na+-free form is often referred to as 'slow' because of reduced rates of cleavage of procoagulant substrates, but the slow form is still capable of rapid activation of protein C in the presence of TM. The molecular basis of the slow proteolytic activity of thrombin has remained elusive, in spite of two decades of solution studies and many published crystallographic structures. In the present paper, we report the first structure of wild-type unliganded human thrombin grown in the absence of co-ordinating Na+. The Na+-binding site is observed in a highly ordered position 6 A (1 A=0.1 nm) removed from that seen in the Na+-bound state. The movement of the Na+ loop results in non-catalytic hydrogen-bonding in the active site and blocking of the S1 and S2 substrate-binding pockets. Similar, if more dramatic, changes were observed in a previous structure of the constitutively slow thrombin variant E217K. The slow behaviour of thrombin in solutions devoid of Na+ can now be understood in terms of an equilibrium between an inert species, represented by the crystal structure described in the present paper, and an active form, where the addition of Na+ populates the active state.

The natural mutation by deletion of Lys9 in the thrombin A-chain affects the pKa value of catalytic residues, the overall enzyme's stability and conformational transitions linked to Na+ binding

The catalytic competence of the natural thrombin mutant with deletion of the Lys9 residue in the A-chain (deltaK9) was found to be severely impaired, most likely due to modification of the 60-loop conformation and catalytic triad geometry, as supported by long molecular dynamics (MD) simulations in explicit water solvent. In this study, the pH dependence of the catalytic activity and binding of the low-molecular mass inhibitor N-alpha-(2-naphthylsulfonyl-glycyl)-4-amidinophenylalanine-piperidine (alpha-NAPAP) to the wild-type (WT) and deltaK9 thrombin forms were investigated, along with their overall structural stabilities and conformational properties. Two ionizable groups were found to similarly affect the activity of both thrombins. The pKa value of the first ionizable group, assigned to the catalytic His57 residue, was found to be 7.5 and 6.9 in ligand-free deltaK9 and WT thrombin, respectively. Urea-induced denaturation studies showed higher instability of the deltaK9 mutant compared with WT thrombin, and disulfide scrambling experiments proved weakening of the interchain interactions, causing faster release of the reduced A-chain in the mutant enzyme. The sodium ion binding affinity was not significantly perturbed by Lys9 deletion, although the linked increase in intrinsic fluorescence was lower in the mutant. Essential dynamics (ED) analysis highlighted different conformational properties of the two thrombins in agreement with the experimental conformational stability data. Globally, these findings enhanced our understanding of the perturbations triggered by Lys9 deletion, which reduces the overall stability of the molecule, weakens the A-B interchain interactions, and allosterically perturbs the geometry and protonation state of catalytic residues of the enzyme.

Mapping the energetics of water-protein and water-ligand interactions with the "natural" HINT forcefield: predictive tools for characterizing the roles of water in biomolecules

The energetics and hydrogen bonding pattern of water molecules bound to proteins were mapped by analyzing structural data (resolution better than 2.3A) for sets of uncomplexed and ligand-complexed proteins. Water-protein and water-ligand interactions were evaluated using hydropatic interactions (HINT), a non-Newtonian forcefield based on experimentally determined logP(octanol/water) values. Potential water hydrogen bonding ability was assessed by a new Rank algorithm. The HINT-derived binding energies and Ranks for second shell water molecules were -0.04 kcal mol(-1) and 0.0, respectively, for first shell water molecules -0.38 kcal mol(-1) and 1.6, for active site water molecules -0.45 kcal mol(-1) and 2.3, for cavity water molecules -0.55 kcal mol(-1) and 3.3, and for buried water molecules -0.56 kcal mol(-1) and 4.4. For the last four classes, similar energies indicate that internal and external water molecules interact with protein almost equally, despite different degrees of hydrogen bonding. The binding energies and Ranks for water molecules bridging ligand-protein were -1.13 kcal mol(-1) and 4.5, respectively. This energetic contribution is shared equally between protein and ligand, whereas Rank favors the protein. Lastly, by comparing the uncomplexed and complexed forms of proteins, guidelines were developed for prediction of the roles played by active site water molecules in ligand binding. A water molecule with high Rank and HINT score is unlikely to make further interactions with the ligand and is largely irrelevant to the binding process, while a water molecule with moderate Rank and high HINT score is available for ligand interaction. Water molecule displaced for steric reasons were characterized by lower Rank and HINT score. These guidelines, tested by calculating HINT score and Rank for 50 water molecules bound in the active site of four uncomplexed proteins (for which the structures of the liganded forms were also available), correctly predicted the ultimate roles (in the complex) for 76% of water molecules. Some failures were likely due to ambiguities in the structural data.

Effect of Na+ binding on the conformation, stability and molecular recognition properties of thrombin

In the present work, the effect of Na+ binding on the conformational, stability and molecular recognition properties of thrombin was investigated. The binding of Na+ reduces the CD signal in the far-UV region, while increasing the intensity of the near-UV CD and fluorescence spectra. These spectroscopic changes have been assigned to perturbations in the environment of aromatic residues at the level of the S2 and S3 sites, as a result of global rigidification of the thrombin molecule. Indeed, the Na+-bound form is more stable to urea denaturation than the Na+-free form by approximately 2 kcal/mol (1 cal identical with 4.184 J). Notably, the effects of cation binding on thrombin conformation and stability are specific to Na+ and parallel the affinity order of univalent cations for the enzyme. The Na+-bound form is even more resistant to limited proteolysis by subtilisin, at the level of the 148-loop, which is suggestive of the more rigid conformation this segment assumes in the 'fast' form. Finally, we have used hirudin fragment 1-47 as a molecular probe of the conformation of thrombin recognition sites in the fast and 'slow' form. From the effects of amino acid substitutions on the affinity of fragment 1-47 for the enzyme allosteric forms, we concluded that the specificity sites of thrombin in the Na+-bound form are in a more open and permissible conformation, compared with the more closed structure they assume in the slow form. Taken together, our results indicate that the binding of Na+ to thrombin serves to stabilize the enzyme into a more open and rigid conformation.

Slow thrombin in solution

As a tool for understanding biological mechanisms, X-ray crystallography possesses unparalleled power to enlighten, resolve controversy and shift a field of study on to a secure new paradigm. Thanks largely to developments in crystallographic methods, the technique has become accessible to the general biochemist and we have thus witnessed an exponential increase in the number of protein structures deposited every year. It is now commonplace for several structures to be published of the same protein under different crystallization conditions, sometimes resulting in conflicting mechanistic interpretations. Such a controversy has arisen over thrombin's conformational response to Na+ co-ordination, and in this issue of Biochemical Journal, De Filippis and colleagues put the two structural models of thrombin allostery to the test by returning to the techniques of solution biochemistry.

Na+ site in blood coagulation factor IXa: effect on catalysis and factor VIIIa binding

During blood coagulation, factor IXa (FIXa) activates factor X (FX) requiring Ca2+, phospholipid, and factor VIIIa (FVIIIa). The serine protease domain of FIXa contains a Ca2+ site and is predicted to contain a Na+ site. Comparative homology analysis revealed that Na+ in FIXa coordinates to the carbonyl groups of residues 184A, 185, 221A, and 224 (chymotrypsin numbering). Kinetic data obtained at several concentrations of Na+ and Ca2+ with increasing concentrations of a synthetic substrate (CH3-SO2-d-Leu-Gly-Arg-p-nitroanilide) were fit globally, assuming rapid equilibrium conditions. Occupancy by Na+ increased the affinity of FIXa for the synthetic substrate, whereas occupancy by Ca2+ decreased this affinity but increased k(cat) dramatically. Thus, Na+-FIXa-Ca2+ is catalytically more active than free FIXa. FIXa(Y225P), a Na+ site mutant, was severely impaired in Na+ potentiation of its catalytic activity and in binding to p-aminobenzamidine (S1 site probe) validating that substrate binding in FIXa is linked positively to Na+ binding. Moreover, the rate of carbamylation of NH2 of Val16, which forms a salt-bridge with Asp194 in serine proteases, was faster for FIXa(Y225P) and addition of Ca2+ overcame this impairment only partially. Further studies were aimed at delineating the role of the FIXa Na+ site in macromolecular catalysis. In the presence of Ca2+ and phospholipid, with or without saturating FVIIIa, FIXa(Y225P) activated FX with similar K(m) but threefold reduced k(cat). Further, interaction of FVIIIa:FIXa(Y225P) was impaired fourfold. Our previous data revealed that Ca2+ binding to the protease domain increases the affinity of FIXa for FVIIIa approximately 15-fold. The present data indicate that occupancy of the Na+ site further increases the affinity of FIXa for FVIIIa fourfold and k(cat) threefold. Thus, in the presence of Ca2+, phospholipid, and FVIIIa, binding of Na+ to FIXa increases its biologic activity by approximately 12-fold, implicating its role in physiologic coagulation.

Snake venom serine proteinases: sequence homology vs. substrate specificity, a paradox to be solved

Snake venom glands synthesize a variety of serine proteinases capable of affecting the haemostatic system. They act on macromolecular substrates of the coagulation, fibrinolytic, and kallikrein-kinin systems, and on platelets to cause an imbalance of the haemostatic system of the prey. In this review we describe their biochemical/biophysical characteristics, biological activities as well as aspects of their evolution and structure-activity relationship.