Identification of residues linked to the slow-->fast transition of thrombin

Enhancing the anticoagulant profile of meizothrombin

Meizothrombin is an active intermediate generated during the proteolytic activation of prothrombin to thrombin in the penultimate step of the coagulation cascade. Structurally, meizothrombin differs from thrombin because it retains the auxiliary Gla domain and two kringles. Functionally, meizothrombin shares with thrombin the ability to cleave procoagulant (fibrinogen), prothrombotic (PAR1) and anticoagulant (protein C) substrates, although its specificity toward fibrinogen and PAR1 is less pronounced. In this study we report information on the structural architecture of meizothrombin resolved by SAXS and single molecule FRET as an elongated arrangement of its individual domains. In addition, we show the properties of a meizothrombin construct analogous to the anticoagulant thrombin mutant W215A/E217A currently in Phase I for the treatment of thrombotic complications and stroke. The findings reveal new structural and functional aspects of meizothrombin that advance our understanding of a key intermediate of the prothrombin activation pathway.

Histone H4 promotes prothrombin autoactivation

Recent studies have documented the ability of prothrombin to spontaneously convert to the mature protease thrombin when Arg-320 becomes exposed to solvent for proteolytic attack upon mutation of residues in the activation domain. Whether prothrombin autoactivation occurs in the wild-type under conditions relevant to physiology remains unknown. Here, we report that binding of histone H4 to prothrombin under physiological conditions generates thrombin by autoactivation. The effect is abrogated by mutation of the catalytic Ser-525 and requires the presence of the Gla domain. Fluorescence titrations document direct binding of histone H4 to prothrombin with an affinity in the low nm range. Stopped flow data and luminescence resonance energy transfer measurements indicate that the binding mechanism obeys conformational selection. Among the two conformations of prothrombin, collapsed and fully extended, histone H4 binds selectively to the collapsed form and induces a transition toward a new conformation where the distance between Ser-101 in kringle-1 and Ser-210 in kringle-2 increases by 13 Å. These findings confirm the molecular plasticity of prothrombin emerged from recent structural studies and suggest that different conformations of the inter-kringle linker domain determine the functional behavior of prothrombin. The results also broaden our mechanistic understanding of the prothrombotic phenotype observed during cellular damage due to the release of histones in the blood stream. Prothrombin autoactivation induced by histone H4 emerges as a mechanism of pathophysiological relevance through which thrombin is generated independently of activation of the coagulation cascade.

Crystal structure of prothrombin reveals conformational flexibility and mechanism of activation

The zymogen prothrombin is composed of fragment 1 containing a Gla domain and kringle-1, fragment 2 containing kringle-2, and a protease domain containing A and B chains. The prothrombinase complex assembled on the surface of platelets converts prothrombin to thrombin by cleaving at Arg-271 and Arg-320. The three-dimensional architecture of prothrombin and the molecular basis of its activation remain elusive. Here we report the first x-ray crystal structure of prothrombin as a Gla-domainless construct carrying an Ala replacement of the catalytic Ser-525. Prothrombin features a conformation 80 Å long, with fragment 1 positioned at a 36° angle relative to the main axis of fragment 2 coaxial to the protease domain. High flexibility of the linker connecting the two kringles suggests multiple arrangements for kringle-1 relative to the rest of the prothrombin molecule. Luminescence resonance energy transfer measurements detect two distinct conformations of prothrombin in solution, in a 3:2 ratio, with the distance between the two kringles either fully extended (54 ± 2 Å) or partially collapsed (≤34 Å) as seen in the crystal structure. A molecular mechanism of prothrombin activation emerges from the structure. Of the two sites of cleavage, Arg-271 is located in a disordered region connecting kringle-2 to the A chain, but Arg-320 is well defined within the activation domain and is not accessible to proteolysis in solution. Burial of Arg-320 prevents prothrombin autoactivation and directs prothrombinase to cleave at Arg-271 first. Reversal of the local electrostatic potential then redirects prothrombinase toward Arg-320, leading to thrombin generation via the prethrombin-2 intermediate.

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.

Important role of the cys-191 cys-220 disulfide bond in thrombin function and allostery

Little is known on the role of disulfide bonds in the catalytic domain of serine proteases. The Cys-191-Cys-220 disulfide bond is located between the 190 strand leading to the oxyanion hole and the 220-loop that contributes to the architecture of the primary specificity pocket and the Na+ binding site in allosteric proteases. Removal of this bond in thrombin produces an approximately 100-fold loss of activity toward several chromogenic and natural substrates carrying Arg or Lys at P1. Na+ activation is compromised, and no fluorescence change can be detected in response to Na+ binding. A 1.54-A resolution structure of the C191A/C220A mutant in the free form reveals a conformation similar to the Na+-free slow form of wild type. The lack of disulfide bond exposes the side chain of Asp-189 to solvent, flips the backbone O atom of Gly-219, and generates disorder in portions of the 186 and 220 loops defining the Na+ site. This conformation, featuring perturbation of the Na+ site but with the active site accessible to substrate, offers a possible representation of the recently identified E* form of thrombin. Disorder in the 186 and 220 loops and the flip of Gly-219 are corrected by the active site inhibitor H-D-Phe-Pro-Arg-CH(2)Cl, as revealed by the 1.8-A resolution structure of the complex. We conclude that the Cys-191-Cys-220 disulfide bond confers stability to the primary specificity pocket by shielding Asp-189 from the solvent and orients the backbone O atom of Gly-219 for optimal substrate binding. In addition, the disulfide bond stabilizes the 186 and 220 loops that are critical for Na+ binding and activation.

Rapid kinetics of Na+ binding to thrombin

The kinetic mechanism of Na(+) binding to thrombin was resolved by stopped-flow measurements of intrinsic fluorescence. Na(+) binds to thrombin in a two-step mechanism with a rapid phase occurring within the dead time of the spectrometer (<0.5 ms) followed by a single-exponential slow phase whose k(obs) decreases hyperbolically with increasing [Na(+)]. The rapid phase is due to Na(+) binding to the enzyme E to generate the E:Na(+) form. The slow phase is due to the interconversion between E(*) and E, where E(*) is a form that cannot bind Na(+). Temperature studies in the range from 5 to 35 degrees C show significant enthalpy, entropy, and heat capacity changes associated with both Na(+) binding and the E to E(*) transition. As a result, under conditions of physiologic temperature and salt concentrations, the E(*) form is negligibly populated (<1%) and thrombin is almost equally partitioned between the E (40%) and E:Na(+) (60%) forms. Single-site Phe mutations of all nine Trp residues of thrombin enabled assignment of the fluorescence changes induced by Na(+) binding mainly to Trp-141 and Trp-215, and to a lesser extent to Trp-148, Trp-207, and Trp-237. However, the fast phase of fluorescence increase is influenced to different extents by all Trp residues. The distribution of these residues over the entire thrombin surface demonstrates that Na(+) binding induces long-range effects on the structure of the enzyme as a whole, contrary to the conclusions drawn from recent structural studies. These findings elucidate the mechanism of Na(+) binding to thrombin and are relevant to other clotting factors and enzymes allosterically activated by monovalent cations.

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.

The Structural Integrity of Anion Binding Exosite I of Thrombin Is Required and Sufficient for Timely Cleavage and Activation of Factor V and Factor VIII

Alpha-thrombin has two separate electropositive binding exosites (anion binding exosite I, ABE-I and anion binding exosite II, ABE-II) that are involved in substrate tethering necessary for efficient catalysis. Alpha-thrombin catalyzes the activation of factor V and factor VIII following discrete proteolytic cleavages. Requirement for both anion binding exosites of the enzyme has been suggested for the activation of both procofactors by alpha-thrombin. We have used plasma-derived alpha-thrombin, beta-thrombin (a thrombin molecule that has only ABE-II available), and a recombinant prothrombin molecule rMZ-II (R155A/R284A/R271A) that can only be cleaved at Arg(320) (resulting in an enzymatically active molecule that has only ABE-I exposed, rMZ-IIa) to ascertain the role of each exosite for procofactor activation. We have also employed a synthetic sulfated pentapeptide (DY(SO(3)(-))DY(SO(3)(-))Q, designated D5Q1,2) as an exosite-directed inhibitor of thrombin. The clotting time obtained with beta-thrombin was increased by approximately 8-fold, whereas rMZ-IIa was 4-fold less efficient in promoting clotting than alpha-thrombin under similar experimental conditions. Alpha-thrombin readily activated factor V following cleavages at Arg(709), Arg(1018), and Arg(1545) and factor VIII following proteolysis at Arg(372), Arg(740), and Arg(1689). Cleavage of both procofactors by alpha-thrombin was significantly inhibited by D5Q1,2. In contrast, beta-thrombin was unable to cleave factor V at Arg(1545) and factor VIII at both Arg(372) and Arg(1689). The former is required for light chain formation and expression of optimum factor Va cofactor activity, whereas the latter two cleavages are a prerequisite for expression of factor VIIIa cofactor activity. Beta-thrombin was found to cleave factor V at Arg(709) and factor VIII at Arg(740), albeit less efficiently than alpha-thrombin. The sulfated pentapeptide inhibited moderately both cleavages by beta-thrombin. Under similar experimental conditions, membrane-bound rMZ-IIa cleaved and activated both procofactor molecules. Activation of the two procofactors by membrane-bound rMZ-IIa was severely impaired by D5Q1,2. Overall the data demonstrate that ABE-I alone of alpha-thrombin can account for the interaction of both procofactors with alpha-thrombin resulting in their timely and efficient activation. Because formation of meizothrombin precedes that of alpha-thrombin, our findings also imply that meizothrombin may be the physiological activator of both procofactors in vivo in the presence of a procoagulant membrane surface during the early stages of coagulation.

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.

Chemically modified thrombin and anhydrothrombin that differentiate macromolecular substrates of thrombin

BACKGROUND

Thrombin is a primary inducer of thrombus formation by activations of coagulation cascade and platelet aggregation. Hitherto, several types of thrombin inhibitors have been developed for therapeutic purpose.

OBJECTIVES

We prepared modified thrombin (M-thrombin) and modified anhydrothrombin (M-anhydrothrombin) by chemical modification of carboxyl groups of thrombin and anhydrothrombin, respectively, to present a new strategy for a potent antiplatelet-anticoagulant agent and new tools for investigation of thrombin functions.

RESULTS

M-anhydrothrombin retained high affinity for factor VIII (FVIII), but demonstrated lower affinity than anhydrothrombin for fibrinogen and factor V (FV). Both M-anhydrothrombin and anhydrothrombin prolonged activated partial thromboplastin time (APTT) without affecting prothrombin time, and M-anhydrothrombin prolonged APTT much more than anhydrothrombin. M-anhydrothrombin also retained affinity for the recombinant extracellular domain peptide of protease-activated receptor 1 (PAR1). M-thrombin exhibited marginal clotting activity (4% of thrombin), but induced platelet aggregation in platelet-rich plasma without forming a fibrin clot, which was completely suppressed by anti-PAR1 antibody (ATAP2) and by M-anhydrothrombin, but not by anhydrothrombin. These results indicate that M-thrombin induced platelet aggregation through the activation of PAR1, and M-anhydrothrombin inhibited this process completely. In contrast, neither M-anhydrothrombin nor anhydrothrombin apparently inhibited thrombin-induced platelet aggregation. Only in the presence of the Gly-Pro-Arg-Pro (GPRP) peptide that inhibits polymerization of fibrin, M-anhydrothrombin completely inhibited thrombin-induced platelet aggregation.

CONCLUSION

M-thrombin is PAR1-specific and M-anhydrothrombin is FVIII- and PAR1-specific derivatives, and thereby, are new tools as specific agonist and antagonist, respectively, of PAR1. Furthermore, M-anhydrothrombin may be an attractive model for development of a potent anticoagulant-antiplatelet agent.

Hirudin Binding Reveals Key Determinants of Thrombin Allostery

Thrombin exists in two allosteric forms, slow (S) and fast (F), that recognize natural substrates and inhibitors with significantly different affinities. Because under physiologic conditions the two forms are almost equally populated, investigation of thrombin function must address the contribution from the S and F forms and the molecular origin of their differential recognition of ligands. Using a panel of 79 Ala mutants, we have mapped for the first time the epitopes of thrombin recognizing a macromolecular ligand, hirudin, in the S and F forms. Hirudin binding is a relevant model for the interaction of thrombin with fibrinogen and PAR1 and is likewise influenced by the allosteric S-->F transition. The epitopes are nearly identical and encompass two hot spots, one in exosite I and the other in the Na+ site at the opposite end of the protein. The higher affinity of the F form is due to the preferential interaction of hirudin with Lys-36, Leu-65, Thr-74, and Arg-75 in exosite I; Gly-193 in the oxyanion hole; and Asp-221 and Asp-222 in the Na+ site. Remarkably, no correlation is found between the energetic and structural involvements of thrombin residues in hirudin recognition, which invites caution in the analysis of protein-protein interactions in general.

Allosteric modulation of monomeric proteins*

Multimeric proteins (e.g. hemoglobin) are considered to be the prototypes of allosteric enzymes, whereas monomeric proteins (e.g. myoglobin) usually are assumed to be nonallosteric. However, the modulation of the functional properties of monomeric proteins by heterotropic allosteric effectors casts doubts on this assumption. Here, the allosteric properties of sperm whale myoglobin, human serum albumin, and human α-thrombin, generally considered as molecular models of monomeric proteins, are summarized.

Fibrinogen and fibrin

Fibrinogen is a large, complex, fibrous glycoprotein with three pairs of polypeptide chains linked together by 29 disulfide bonds. It is 45 nm in length, with globular domains at each end and in the middle connected by alpha-helical coiled-coil rods. Both strongly and weakly bound calcium ions are important for maintenance of fibrinogen's structure and functions. The fibrinopeptides, which are in the central region, are cleaved by thrombin to convert soluble fibrinogen to insoluble fibrin polymer, via intermolecular interactions of the "knobs" exposed by fibrinopeptide removal with "holes" always exposed at the ends of the molecules. Fibrin monomers polymerize via these specific and tightly controlled binding interactions to make half-staggered oligomers that lengthen into protofibrils. The protofibrils aggregate laterally to make fibers, which then branch to yield a three-dimensional network-the fibrin clot-essential for hemostasis. X-ray crystallographic structures of portions of fibrinogen have provided some details on how these interactions occur. Finally, the transglutaminase, Factor XIIIa, covalently binds specific glutamine residues in one fibrin molecule to lysine residues in another via isopeptide bonds, stabilizing the clot against mechanical, chemical, and proteolytic insults. The gene regulation of fibrinogen synthesis and its assembly into multichain complexes proceed via a series of well-defined steps. Alternate splicing of two of the chains yields common variant molecular isoforms. The mechanical properties of clots, which can be quite variable, are essential to fibrin's functions in hemostasis and wound healing. The fibrinolytic system, with the zymogen plasminogen binding to fibrin together with tissue-type plasminogen activator to promote activation to the active enzyme plasmin, results in digestion of fibrin at specific lysine residues. Fibrin(ogen) also specifically binds a variety of other proteins, including fibronectin, albumin, thrombospondin, von Willebrand factor, fibulin, fibroblast growth factor-2, vascular endothelial growth factor, and interleukin-1. Studies of naturally occurring dysfibrinogenemias and variant molecules have increased our understanding of fibrinogen's functions. Fibrinogen binds to activated alphaIIbbeta3 integrin on the platelet surface, forming bridges responsible for platelet aggregation in hemostasis, and also has important adhesive and inflammatory functions through specific interactions with other cells. Fibrinogen-like domains originated early in evolution, and it is likely that their specific and tightly controlled intermolecular interactions are involved in other aspects of cellular function and developmental biology.

S1 subsite in snake venom thrombin-like enzymes: can S1 subsite lipophilicity be used to sort binding affinities of trypsin-like enzymes to small-molecule inhibitors?

Thrombin-like enzymes isolated from snake venoms comprise a group of serine proteinases responsible for many important coagulation disorders in the envenomed victims. Besides, these proteinases have great biotechnological interest as antithrombotic agents and as diagnostic tools. However, in spite of the recent overflow of snake venom thrombin-like enzymes (SVTLEs) on protein sequence databases, there is a lack of three-dimensional (3D) structural information on this family. Without such 3D structures available many aspects of the biological function and biochemical properties of these enzymes still remain obscure. Therefore, we have gone through a series of computational techniques, which enabled us to identify the set of residues involved in molecular recognition of inhibitors bound to the S1 subsite of snake venom thrombin-like enzymes (SVTLEs) and ultimately conclude that nonpolar (van der Waals) intermolecular interactions and ligand's hydrophobicity are the most important factors affecting binding affinities to the S1 subsite of a SVTLE isolated from the venom of Lachesis muta muta (Lmm-TLE). Consequently, we have proposed that S1 subsite lipophilicity may be used to sort binding affinities of trypsin-like enzymes to small molecules by showing that the inhibitory potency of several S1-directed compounds follows subsite lipophilicity among Lmm-TLE and other three homologous proteases. Noteworthy, in the course of our analyses we determined that thrombin's S1 subsite should, in fact, be considered less lipophilic than that of trypsin if we account for the presence of the sodium-controlled water channel communicating with the S1 subsite in the coagulant enzyme.

Crystal structure of the thrombin mutant D221A/D222K: the Asp222:Arg187 ion-pair stabilizes the fast form

The thrombin mutant D221A/D222K (ARK) does not bind Na+ and has interesting functional properties intermediate between those of the slow and fast forms of wild type. We solved the X-ray crystal structure of ARK bound at exosite I with a fragment of hirudin at 2.4-A resolution. The structure shows a slight collapse of the 186 and 220 loops into the Na+ binding site due to disruption of the Asp222:Arg187 ion-pair. The backbone O atoms of Arg221a and Lys224 are shifted into conformations that eliminate optimal interaction with Na+. A paucity of solvent molecules in the Na+ binding site is also noted, by analogy to what is seen in the structure of the slow form. These findings reinforce the crucial role of the Asp222:Arg187 ion-pair in stabilizing the fast form of thrombin.

Generation of ligand conformations in continuum solvent consistent with protein active site topology: application to thrombin

Using the crystal structure of an inhibitor complexed with the serine protease thrombin (PDB code ) and the functional group definitions contained within the Catalyst software, a representation of the enzyme's active site was produced (structure-based pharmacophore model). A training set of 16 homologous non-peptide inhibitors whose conformations had been generated in continuum solvent (MacroModel) and clustered into conformational families (XCluster) was regressed against this pharmacophore so as to obtain a 3D-QSAR model. To test the robustness of the resulting QSAR model, the synthesis of a series of non-peptide thrombin inhibitors based on arylsuphonyl derivatives of an aminophenol ring linked to a pyridyl-based S1 binding group was undertaken. These compounds served as a test set (20-24). The crystal structure for the novel symmetrical disulfonyl compound 24, in complex with thrombin, has been solved. Its calculated binding mode is in general agreement with the crystallographically observed one, and the predicted K(i) value is in close accord with the experimental value.

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

Thrombin participates in its own positive and negative feedback loops, and its allosteric state helps determine the hemostatic balance. Here we present the 1.8 A crystallographic structure of S195A thrombin in two conformational states: active site occupied and active site free. The active site-occupied form shows how thrombin can accommodate substrates, such as protein C. The active site-free form is in a previously unobserved closed conformation of thrombin, which satisfies all the conditions of the so-called "slow" form. A mechanism of allostery is revealed, which relies on the concerted movement of the disulphide bond between Cys168 and 182 and aromatic residues Phe227, Trp215, and Trp60d. These residues constitute an allosteric switch, which is flipped directly through sodium binding, resulting in the fast form with an open active site.

Allosteric effects potentiating the release of the second fibrinopeptide A from fibrinogen by thrombin

Fibrin formation depends on the release of the two N-terminal fibrinopeptides A (FPA) from fibrinogen, and its formation is accompanied by an intermediate, alpha-profibrin, which lacks only one of the FPA. In this study, we confirm that the maximal levels of alpha-profibrin found over the course of thrombin reactions with human fibrinogen are only half of what would be expected if the first and second FPA were being released independently with equal rate constants. The rapidity of release of the fibrinopeptides by thrombin had been shown to depend on an allosteric transformation that is induced when Na(+) binds to a site defined by the 215-227 residues of thrombin, a transformation that results in the exposure of its fibrinogen-binding exosites transforming the thrombin from a slow to a fast acting form toward fibrinogen. When choline was substituted for sodium to transform thrombin to its slow form, the maximal levels of alpha-profibrin rose to those expected for independent release of the two FPA. Thus, it is only the fast thrombin that releases the second FPA fast, and that fast release only occurs when both FPA are present because of a partial coupling of its release with that of the first FPA. The release of the FPA from purified alpha-profibrin with the first FPA already missing is no faster than the release of any FPA. Surprisingly, we also found that slow thrombin became increasingly transformed to a fast form in the absence of sodium when the fibrinogen was elevated to high concentrations. This potentiation by concentrated fibrinogen also occurs with the recombinant mutant thrombin (Y225P), which is otherwise slow in both the presence and absence of Na(+). The potentiation of thrombin by fibrinogen must be short-lived so that the thrombin reverts to its slow acting form in the interim among encounters with other fibrinogen molecules in dilute fibrinogen solutions lacking Na(+), whereas at high fibrinogen concentrations the thrombin encounters other molecules before it reverts back to the slow form.

Localization of phosphatidylserine binding sites to structural domains of factor Xa

Binding of short chain phosphatidylserine (C6PS) enhances the proteolytic activity of factor X(a) by 60-fold (Koppaka, V., Wang, J., Banerjee, M., and Lentz, B. R. (1996) Biochemistry 35, 7482-7491). In the present study, we locate three C6PS binding sites to different domains of factor X(a) using a combination of activity, circular dichroism, fluorescence, and equilibrium dialysis measurements on proteolytic and biosynthetic fragments of factor X(a). Our results demonstrate that the structural responses of human and bovine factor X(a) to C6PS binding are somewhat different. Despite this difference, data obtained with fragments from both human and bovine factor X(a) are consistent with a common hypothesis for the location of C6PS binding sites to different structural domains. First, the gamma-carboxyglutamic acid (Gla) domain binds C6PS only in the absence of Ca(2+) (k(d) approximately 1 mm), although this PS site does not influence the functional response of factor X(a). Second, a Ca(2+)-dependent binding site is in the epidermal growth factor domains (EGF(NC)) that are linked by Ca(2+) and C6PS binding to the Gla domain. This site appears to be the lipid regulatory site of factor X(a). Third, a Ca(2+)-requiring site seems to be in the EGF(C)-catalytic domain. This site appears not to be a lipid regulatory site but rather to share residues with the substrate recognition site. Finally, the full functional response to C6PS requires linkage of the Gla, EGF(NC), and catalytic domains in the presence of Ca(2+), meaning that PS regulation of factor X(a) involves linkage between widely separated parts of the protein.

Determinants of thrombin specificity

Thrombin recognizes a number of natural substrates that are responsible for important physiologic functions. Its high specificity is controlled by residues within the active site, and by separate recognition sites located on the surface of the enzyme. A number of studies have addressed the question of how thrombin changes its specificity from fibrinogen to protein C, switching from a procoagulant to an anticoagulant enzyme. Site directed mutagenesis studies have revealed important aspects of how this switch takes place. Specifically, residues W215 and E217 have emerged as key residues in controlling the interaction with fibrinogen in that mutation of these residues compromises the procoagulant function of the enzyme up to 500-fold. The loss of fibrinogen clotting reaches 20,000-fold in the double mutant W215A/E217A, whereas protein C activation is compromised less than sevenfold. These findings demonstrate that thrombin specificity can be dissected at the molecular level using Ala-scanning mutagenesis and the procoagulant function of the enzyme can be abrogated rationally and selectively. It is now possible to extend this strategy to the study of other interactions of thrombin, as well as to related serine proteases.

An extensive interaction interface between thrombin and factor V is required for factor V activation

The interaction interface between human thrombin and human factor V (FV), necessary for complex formation and cleavage to generate factor Va, was investigated using a site-directed mutagenesis strategy. Fifty-three recombinant thrombins, with a total of 78 solvent-exposed basic and polar residues substituted with alanine, were used in a two-stage clotting assay with human FV. Seventeen mutants with less than 50% of wild-type (WT) thrombin FV activation were identified and mapped to anion-binding exosite I (ABE-I), anion-binding exosite II (ABE-II), the Leu(45)-Asn(57) insertion loop, and the Na(+) binding loop of thrombin. Three ABE-I mutants (R68A, R70A, and Y71A) and the ABE-II mutant R98A had less than 30% of WT activity. The thrombin Na(+) binding loop mutants, E229A and R233A, and the Leu(45)-Asn(57) insertion loop mutant, W50A, had a major effect on FV activation with 5, 15, and 29% of WT activity, respectively. The K52A mutant, which maps to the S' specificity pocket, had 29% of WT activity. SDS-polyacrylamide gel electrophoresis analysis of cleavage reactions using the thrombin ABE mutants R68A, Y71A, and R98A, the Na(+) binding loop mutant E229A, and the Leu(45)-Asn(57) insertion loop mutant W50A showed a requirement for both ABEs and the Na(+)-bound form of thrombin for efficient cleavage at the FV residue Arg(709). Several basic residues in both ABEs have moderate decreases in FV activation (40-60% of WT activity), indicating a role for the positive electrostatic fields generated by both ABEs in enhancing complex formation with complementary negative electrostatic fields generated by FV. The data show that thrombin activation of FV requires an extensive interaction interface with thrombin. Both ABE-I and ABE-II and the S' subsite are required for optimal cleavage, and the Na(+)-bound form of thrombin is important for its procoagulant activity.

Unexpected crucial role of residue 225 in serine proteases

Residue 225 in serine proteases of the chymotrypsin family is Pro or Tyr in more than 95% of nearly 300 available sequences. Proteases with Y225 (like some blood coagulation and complement factors) are almost exclusively found in vertebrates, whereas proteases with P225 (like degradative enzymes) are present from bacteria to human. Saturation mutagenesis of Y225 in thrombin shows that residue 225 affects ligand recognition up to 60,000-fold. With the exception of Tyr and Phe, all residues are associated with comparable or greatly reduced catalytic activity relative to Pro. The crystal structures of three mutants that differ widely in catalytic activity (Y225F, Y225P, and Y225I) show that although residue 225 makes no contact with substrate, it drastically influences the shape of the water channel around the primary specificity site. The activity profiles obtained for thrombin also suggest that the conversion of Pro to Tyr or Phe documented in the vertebrates occurred through Ser and was driven by a significant gain (up to 50-fold) in catalytic activity. In fact, Ser and Phe are documented in 4% of serine proteases, which together with Pro and Tyr account for almost the entire distribution of residues at position 225. The unexpected crucial role of residue 225 in serine proteases explains the evolutionary selection of residues at this position and shows that the structural determinants of protease activity and specificity are more complex than currently believed. These findings have broad implications in the rational design of enzymes with enhanced catalytic properties.

Identification and characterization of the sodium-binding site of activated protein C

Activated protein C (APC) requires both Ca2+ and Na+ for its optimal catalytic function. In contrast to the Ca2+-binding sites, the Na+-binding site(s) of APC has not been identified. Based on a recent study with thrombin, the 221-225 loop is predicted to be a potential Na+-binding site in APC. The sequence of this loop is not conserved in trypsin. We engineered a Gla domainless form of protein C (GDPC) in which the 221-225 loop was replaced with the corresponding loop of trypsin. We found that activated GDPC (aGDPC) required Na+ (or other alkali cations) for its amidolytic activity with dissociation constant (Kd(app)) = 44.1 +/- 8.6 mM. In the presence of Ca2+, however, the requirement for Na+ by aGDPC was eliminated, and Na+ stimulated the cleavage rate 5-6-fold with Kd(app) = 2.3 +/- 0.3 mM. Both cations were required for efficient factor Va inactivation by aGDPC. In the presence of Ca2+, the catalytic function of the mutant was independent of Na+. Unlike aGDPC, the mutant did not discriminate among monovalent cations. We conclude that the 221-225 loop is a Na+-binding site in APC and that an allosteric link between the Na+ and Ca2+ binding loops modulates the structure and function of this anticoagulant enzyme.

Allosteric modulation of BPTI interaction with human alpha- and zeta-thrombin

In this study, thrombin interaction with the basic pancreatic trypsin inhibitor (BPTI) was investigated in the presence of different allosteric modulators of thrombin, that is the C-terminal hirudin peptide 54-65 (Hir54-65), a recombinant thrombomodulin form (TMEGF4-6) and Na+. BPTI binding to alpha-thrombin is positively linked to Na+. Under low sodium concentration (5 mM Na+) the BPTI affinity for alpha-thrombin was roughly threefold lower than in the presence of 150 mM sodium (Ki = 320 microM vs. 100 microM). The hirudin fragment, which binds to the fibrinogen recognition site (FRS) of thrombin, induced a progressive and saturable decrease (3.6-fold) of alpha-thrombin affinity for BPTI, whereas the thrombomodulin peptide, which binds to a more extended region of FRS, caused a 5.5-fold increase of the enzyme affinity for the inhibitor. The opposite effect exerted by Hir54-65 and TMEGF4-6 was also observed for BPTI interaction with zeta-thrombin, in which the amidic bond between W148 and T149 is cleaved. However, in this case the effect by Hir54-65 and TMEGF4-6, although qualitatively similar to that observed with alpha-thrombin, had a smaller magnitude. Thrombin hydrolysis of Protein C was also differently affected by Hir54-65 and TMEGF4-6 peptides. While the latter enhanced the Protein C activation, the former caused a reduction of both alpha- and zeta-thrombin kcat/K(m)' for Protein C cleavage. These results showed that (a) Na+ facilitates BPTI interaction with thrombin; (b) Hir54-65 and TMEGF4-6, though sharing in part the same binding site at the thrombin FRS, can affect in opposite way thrombin's interaction with BPTI and Protein C; (c) such findings along with the results obtained with zeta-thrombin might be explained by admitting that the thermodynamic linkage between FRS and the critical W60-loop is also controlled by ligation and/or conformational state of the W148 insertion loop.

New general approach for determining the solution structure of a ligand bound weakly to a receptor: structure of a fibrinogen Aalpha-like peptide bound to thrombin (S195A) obtained using NOE distance constraints and an ECEPP/3 flexible docking program

A new approach incorporating flexible docking simulations and NMR data is presented for calculating the bound conformation of a ligand that interacts weakly with an enzyme. This approach consists of sampling directly the conformation of a flexible ligand inside a receptor active site containing surrounding flexible loops. To make this sampling efficient, a ligand-growing procedure has been adopted. Optimization of the ECEPP/3-plus-NOE constraint function is carried out by using a collective variable Monte Carlo minimization technique. Numerous energy minimizations are made possible for such a large system by using a Bezier splines energy grid technique. This new flexible docking approach was applied to determine the structure of a fibrinogen Aalpha-like peptide (7DFLAEGGGVRGPRV20) bound to an active site mutant of thrombin [thrombin(S195A)]. Structure calculations of the bound ligand, using 2D-transferred NOESY distance constraints in the DIANA program, showed that the N-terminal portion of the peptide (D7-R16) involves a chain reversal, whereas the C-terminal portion (G17-V20) adopts a fold that exists in several different orientations. In addition, the ECEPP/3 flexible docking package was used to assess the conformational variability of the ligand and surrounding 60D-insertion loop of thrombin. Amino acid residues (17-20) of the peptide interact with a region of the enzyme that exhibits broad specificity, with a preferred direction between the 60D-insertion loop and Pro37 of thrombin.

Role of P225 and the C136-C201 disulfide bond in tissue plasminogen activator

The protease domain of tissue plasminogen activator (tPA), a key fibrinolytic enzyme, was expressed in Escherichia coli with a yield of 1 mg per liter of media. The recombinant protein was titrated with the Erythrina caraffa trypsin inhibitor (ETI) and characterized in its interaction with plasminogen and the natural inhibitor plasminogen activator inhibitor-1 (PAI-1). Analysis of the catalytic properties of tPA using a library of chromogenic substrates carrying substitutions at P1, P2, and P3 reveals a strong preference for Arg over Lys at P1, unmatched by other serine proteases like thrombin or trypsin. In contrast to these proteases and plasmin, tPA shows little or no preference for Pro over Gly at P2. A specific inhibition of tPA by Cu2+ was discovered. The divalent cation presumably binds to H188 near D189 in the primary specificity pocket and inhibits substrate binding in a competitive manner with a Kd = 19 microM. In an attempt to engineer Na+ binding and enhanced catalytic activity in tPA, P225 was replaced with Tyr, the residue present in Na+-dependent allosteric serine proteases. The P225Y mutation did not result in cation binding, but caused a significant loss of specificity (up to 100-fold) toward chromogenic substrates and plasminogen and considerably reduced the inhibition by PAI-1 and ETI. Interestingly, the P225Y substitution enhanced the ability of Cu2+ to inhibit the enzyme. Elimination of the C136-C201 disulfide bond, that is absent in all Na+-dependent allosteric serine proteases, significantly enhanced the yield (5 mg per liter of media) of expression in E. coli, but caused no changes in the properties of the enzyme whether residue 225 was Pro or Tyr. These findings point out an unanticipated crucial role for residue 225 in controlling the catalytic activity of tPA, and suggest that engineering of a Na+-dependent allosteric enhancement of catalytic activity in this enzyme, must involve substantial changes in the region homologous to the Na+ binding site of allosteric serine proteases.

Evidence for allosteric linkage between exosites 1 and 2 of thrombin

Investigations to date have demonstrated that ligand binding to exosites 1 or 2 on thrombin produces conformational changes at the active site. In this study, we directly compared the effect of ligand binding to exosites 1 and 2 on the structure and function of the active site of thrombin and investigated functional linkage between the two exosites. Binding studies were performed in solution with fluorescein-Phe-Pro-Arg-CH2Cl (FPR)-thrombin. Hirudin-(54-65) and sF2, a synthetic peptide corresponding to residues 63-116 of prothrombin fragment 2, were used as ligands for exosites 1 and 2 of thrombin, respectively. The two ligands produce diametric changes in the fluorescence of fluorescein-FPR-thrombin and also have opposing effects on the rate of thrombin hydrolysis of a number of chromogenic substrates. These results indicate that sF2 and hirudin-(54-65) differentially affect the conformation of the active site. Experiments then were performed to investigate whether both ligands can bind to thrombin simultaneously. When thrombin-bound fluorescein-sF2 is titrated with hirudin-(54-65), complete displacement of fluorescein-sF2 is observed. Likewise, when thrombin-bound fluorescein-hirudin-(54-65) is titrated with sF2, complete displacement occurs. Additional support for reciprocal binding was obtained in fluorescence experiments where both probes were labeled and in experiments monitoring ligand binding to agarose-immobilized thrombin. This mutually exclusive binding of either ligand can be explained by reciprocal, allosteric modulation of ligand affinity between the two exosites. Thus, not only do the two exosites differentially influence the active site, they also affect the binding properties of the opposing exosite.

Site-specific dissection of substrate recognition by thrombin

Current approaches to enzyme specificity focus on the identification of consensus sequences from combinatorial chemistry libraries or phage display. These synthetic substrates can also be used as sensitive probes for the molecular environment of the enzyme specificity sites to determine how they contribute to recognition in the transition state. Libraries constructed to include all relevant species for a site-specific analysis contain a relatively small number of substrates and provide quantitative information on the energetics of recognition that can be exploited in studies of structure-function relations and rational drug design. We have constructed a library of substrates carrying substitutions at P1, P2, and P3 to probe the response of the specificity sites S1, S2, and S3 of thrombin. The library has been used to identify differences between the anticoagulant slow and procoagulant fast forms of thrombin and the structural origin of the effects. The results also offer new guidelines for the design of active-site inhibitors of thrombin.

Selective loss of fibrinogen clotting in a loop-less thrombin

The autolysis loop of thrombin comprises nine residues, from Glu146 to Lys149e, five of which (Ala149a-Lys149e) are inserted relative to trypsin and chymotrypsin. Deletion of the insertion Ala149a-Lys149e causes no significant change in the properties of the enzyme, except for a slight enhancement of protein C activation. Deletion of the entire Glu146-Lys149e loop, however, reduces fibrinogen clotting 240-fold, but decreases protein C activation only 2-fold. This loop-less mutant is de facto an exclusive activator of protein C, having lost the primary procoagulant function of thrombin. Because the autolysis loop affects fibrinogen binding, but not protein C activation, it provides a target for new drugs designed to suppress exclusively the procoagulant activity of thrombin.

Rational engineering of activity and specificity in a serine protease

The discovery of the Na(+)-dependent allosteric regulation in serine proteases makes it possible to control catalytic activity and specificity in this class of enzymes in a way never considered before. We demonstrate that rational site-directed mutagenesis of residues controlling Na+ binding can profoundly after the properties of a serine protease. By suppressing Na+ binding to thrombin, we shift the balance between procoagulant and anticoagulant activities of the enzyme. Those mutants, compared to wild-type, have reduced specificity toward fibrinogen, but enhanced or slightly reduced specificity toward protein C. Because this engineering strategy targets a fundamental regulatory mechanism, it is amenable of extension to other enzymes of biological and pharmacological importance.

Residue 225 determines the Na(+)-induced allosteric regulation of catalytic activity in serine proteases

Residue 225 in serine proteases is typically Pro or Tyr and specifies an important and unanticipated functional aspect of this class of enzymes. Proteases with Y225, like thrombin, are involved in highly specialized functions like blood coagulation and complement that are exclusively found in vertebrates. In these proteases, the catalytic activity is enhanced allosterically by Na+ binding. Proteases with P225, like trypsin, are typically involved in digestive functions and are also found in organisms as primitive as eubacteria. These proteases have no requirement for Na+ or other monovalent cations. The molecular origin of this physiologically important difference is remarkably simple and is revealed by a comparison of the Na+ binding loop of thrombin with the homologous region of trypsin. The carbonyl O atom of residue 224 makes a key contribution to the coordination shell of the bound Na+ in thrombin, but is oriented in a manner incompatible with Na+ binding in trypsin because of constraints imposed by P225 on the protein backbone. Pro at position 225 is therefore incompatible with Na+ binding and is a direct predictor of the lack of allosteric regulation in serine proteases. To directly test this hypothesis, we have engineered the thrombin mutant Y225P. This mutant has lost the ability to bind Na+ and behaves like the allosteric slow (Na(+)-free) form. The Na(+)-induced allosteric regulation also bears on the molecular evolution of serine proteases. A strong correlation exists between residue 225 and the codon used for the active site S195. Proteases with P225 typically use a TCN codon for S195, whereas proteases with Y225 use an AGY codon. It is proposed that serine proteases evolved from two main lineages: (i) TCN/P225 with a trypsin-like ancestor and (ii) AGY/Y225 with a thrombin-like ancestor. We predict that the Na(+)-induced allosteric regulation of catalytic activity can be introduced in the TCN/P225 lineage using the P225Y replacement.