The thrombin epitope recognizing thrombomodulin is a highly cooperative hot spot in exosite I

Ligand binding to anion-binding exosites regulates conformational properties of thrombin

Thrombin participates in coagulation, anticoagulation, and initiation of platelet activation. To fulfill its diverse roles and maintain hemostasis, this serine protease is regulated via the extended active site region and anion-binding exosites (ABEs) I and II. For the current project, amide proton hydrogen-deuterium exchange coupled with MALDI-TOF mass spectrometry was used to characterize ligand binding to individual exosites and to investigate the presence of exosite-active site and exosite-exosite interactions. PAR3(44-56) and PAR1(49-62) were observed to bind to thrombin ABE I and then to exhibit long range effects over to ABE II. By contrast, Hirudin(54-65) focused more on ABE I and did not transmit influences over to ABE II. Although these three ligands were each directed to ABE I, they did not promote the same conformational consequences. D-Phe-Pro-Arg-chloromethyl ketone inhibition at the thrombin active site led to further local and long range consequences to thrombin-ABE I ligand complexes with the autolysis loop often most affected. When Hirudin(54-65) was bound to ABE I, it was still possible to bind GpIbα(269-286) or fibrinogen γ'(410-427) to ABE II. Each ligand exerted its predominant influences on thrombin and also allowed interexosite communication. The results obtained support the proposal that thrombin is a highly dynamic protein. The transmission of ligand-specific local and long range conformational events is proposed to help regulate this multifunctional enzyme.

DBAC: a simple prediction method for protein binding hot spots based on burial levels and deeply buried atomic contacts

BACKGROUND

A protein binding hot spot is a cluster of residues in the interface that are energetically important for the binding of the protein with its interaction partner. Identifying protein binding hot spots can give useful information to protein engineering and drug design, and can also deepen our understanding of protein-protein interaction. These residues are usually buried inside the interface with very low solvent accessible surface area (SASA). Thus SASA is widely used as an outstanding feature in hot spot prediction by many computational methods. However, SASA is not capable of distinguishing slightly buried residues, of which most are non hot spots, and deeply buried ones that are usually inside a hot spot.

RESULTS

We propose a new descriptor called "burial level" for characterizing residues, atoms and atomic contacts. Specifically, burial level captures the depth the residues are buried. We identify different kinds of deeply buried atomic contacts (DBAC) at different burial levels that are directly broken in alanine substitution. We use their numbers as input for SVM to classify between hot spot or non hot spot residues. We achieve F measure of 0.6237 under the leave-one-out cross-validation on a data set containing 258 mutations. This performance is better than other computational methods.

CONCLUSIONS

Our results show that hot spot residues tend to be deeply buried in the interface, not just having a low SASA value. This indicates that a high burial level is not only a necessary but also a more sufficient condition than a low SASA for a residue to be a hot spot residue. We find that those deeply buried atoms become increasingly more important when their burial levels rise up. This work also confirms the contribution of deeply buried interfacial atomic contacts to the energy of protein binding hot spot.

Mutagenesis studies toward understanding allostery in thrombin

The binding of thrombomodulin (TM) to exosite-1 and the binding of Na(+) to 225-loop allosterically modulate the catalytic activity and substrate specificity of thrombin. To determine whether the conformation of these two cofactor-binding loops are energetically linked to each other and to the active site, we rationally designed two thrombin mutants in which either the 70-80 loop of exosite-1 or the 225-loop of the Na(+)-binding site was stabilized by an engineered disulfide bond. This was possible by replacing two residues, Arg-67 and Ile-82, in the first mutant and two residues, Glu-217 and Lys-224, in the second mutant with Cys residues. These mutants were expressed in mammalian cells as monomeric molecules, purified to homogeneity and characterized with respect to their ability to bind TM and Na(+) by kinetic and direct binding approaches. The Cys-67/Cys-82 mutant did not bind TM and exhibited a normal amidolytic activity, however, the activity of Cys-217/Cys-224 was dramatically impaired, though TM interacted with this mutant with >20-fold elevated K(D) to partially restore its activity. Both mutants exhibited approximately 2-3-fold higher K(D) for interaction with Na(+), and neither mutant clotted fibrinogen or activated protein C in the presence of TM. Both mutants interacted with heparin with a normal affinity. These results suggest that, while exosite-2 of thrombin is an independent cofactor binding-site, both Na(+)-binding and exosite-1 are energetically linked. Further studies with the fluorescein labeled Cys-195 mutant of thrombin revealed that the catalytic residue of thrombin is modulated by Na(+), but TM has no effect on the conformation of this residue.

Stabilization of the E* form turns thrombin into an anticoagulant

Previous studies have shown that deletion of nine residues in the autolysis loop of thrombin produces a mutant with an anticoagulant propensity of potential clinical relevance, but the molecular origin of the effect has remained unresolved. The x-ray crystal structure of this mutant solved in the free form at 1.55 A resolution reveals an inactive conformation that is practically identical (root mean square deviation of 0.154 A) to the recently identified E* form. The side chain of Trp(215) collapses into the active site by shifting > 10 A from its position in the active E form, and the oxyanion hole is disrupted by a flip of the Glu(192)-Gly(193) peptide bond. This finding confirms the existence of the inactive form E* in essentially the same incarnation as first identified in the structure of the thrombin mutant D102N. In addition, it demonstrates that the anticoagulant profile often caused by a mutation of the thrombin scaffold finds its likely molecular origin in the stabilization of the inactive E* form that is selectively shifted to the active E form upon thrombomodulin and protein C binding.

Volume-based solvation models out-perform area-based models in combined studies of wild-type and mutated protein-protein interfaces

BACKGROUND

Empirical binding models have previously been investigated for the energetics of protein complexation (DeltaG models) and for the influence of mutations on complexation (i.e. differences between wild-type and mutant complexes, DeltaDeltaG models). We construct binding models to directly compare these processes, which have generally been studied separately.

RESULTS

Although reasonable fit models were found for both DeltaG and DeltaDeltaG cases, they differ substantially. In a dataset curated for the absence of mainchain rearrangement upon binding, non-polar area burial is a major determinant of DeltaG models. However this DeltaG model does not fit well to the data for binding differences upon mutation. Burial of non-polar area is weighted down in fitting of DeltaDeltaG models. These calculations were made with no repacking of sidechains upon complexation, and only minimal packing upon mutation. We investigated the consequences of more extensive packing changes with a modified mean-field packing scheme. Rather than emphasising solvent exposure with relatively extended sidechains, rotamers are selected that exhibit maximal packing with protein. This provides solvent accessible areas for proteins that are much closer to those of experimental structures than the more extended sidechain regime. The new packing scheme increases changes in non-polar burial for mutants compared to wild-type proteins, but does not substantially improve agreement between DeltaG and DeltaDeltaG binding models.

CONCLUSION

We conclude that solvent accessible area, based on modelled mutant structures, is a poor correlate for DeltaDeltaG upon mutation. A simple volume-based, rather than solvent accessibility-based, model is constructed for DeltaG and DeltaDeltaG systems. This shows a more consistent behaviour. We discuss the efficacy of volume, as opposed to area, approaches to describe the energetic consequences of mutations at interfaces. This knowledge can be used to develop simple computational screens for binding in comparative modelled interfaces.

N-glycans and the N terminus of protein C inhibitor affect the cofactor-enhanced rates of thrombin inhibition

Protein C inhibitor (PCI) is a serine protease inhibitor, displaying broad protease specificity, found in blood and other tissues. In blood, it is capable of inhibiting both procoagulant and anticoagulant proteases. Mechanisms that provide specificity to PCI remain largely unrevealed. In this study we have for the first time provided a full explanation for the marked size heterogeneity of blood-derived PCI and identified functional differences between naturally occurring PCI variants. The heterogeneity was caused by differences in N-glycan structures, N-glycosylation occupancy, and the presence of a Delta6-N-cleaved form. Bi-, tri-, and tetra-antennary complex N-glycans were identified. Fucose residues were identified both on the core GlcNAc and as parts of sialyl-Le(a/x) epitopes. Moreover, a glycan with a composition that implied a di-sialyl antenna was observed. PCI was N-glycosylated at all three potential N-glycosylation sites, Asn-230, Asn-243, and Asn-319, but a small fraction of PCI lacked the N-glycan at Asn-243. The overall removal of N-glycans affected the maximal heparin- and thrombomodulin-enhanced rates of thrombin inhibition differently in different solution conditions. In contrast, the Delta6-N-region increased both the heparin- and the thrombomodulin-enhanced rates of thrombin inhibition at all conditions examined. These results thus demonstrate that the N-linked glycans and the N-terminal region of blood-derived PCI in different ways affect the cofactor-enhanced rates of thrombin inhibition and provide information on the mechanisms by which this may be achieved. The findings are medically important, in view of the documented association of PCI with atherosclerotic plaques and the promising effect of PCI on reducing hypercoagulability states.

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.

Crystal structures of murine thrombin in complex with the extracellular fragments of murine protease-activated receptors PAR3 and PAR4

It has been proposed that the cleaved form of protease-activated receptor 3 (PAR3) acts as a cofactor for thrombin cleavage and activation of PAR4 on murine platelets, but the molecular basis of this physiologically important effect remains elusive. X-ray crystal structures of murine thrombin bound to extracellular fragments of the murine receptors PAR3 ((38)SFNGGPQNTFEEFPLSDIE(56)) and PAR4 ((51)KSSDKPNPR downward arrow GYPGKFCANDSDTLELPASSQA(81), downward arrow = site of cleavage) have been solved at 2.0 and 3.5 A resolution, respectively. The cleaved form of PAR3, traced in the electron density maps from Gln-44 to Glu-56, makes extensive hydrophobic and electrostatic contacts with exosite I of thrombin through the fragment (47)FEEFPLSDIE(56). Occupancy of exosite I by PAR3 allosterically changes the conformation of the 60-loop and shifts the position of Trp-60d approximately 10 A with a resulting widening of the access to the active site. The PAR4 fragment, traced entirely in the electron density maps except for five C-terminal residues, clamps Trp-60d, Tyr-60a, and the aryl-binding site of thrombin with Pro-56 and Pro-58 at the P2 and P4 positions and engages the primary specificity pocket with Arg-59. The fragment then leaves the active site with Gly-60 and folds into a short helical turn that directs the backbone away from exosite I and over the autolysis loop. The structures demonstrate that thrombin activation of PAR4 may occur with exosite I available to bind cofactor molecules, like the cleaved form of PAR3, whose function is to promote substrate diffusion into the active site by allosterically changing the conformation of the 60-loop.

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.

Thrombin as procoagulant and anticoagulant

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. This basic regulatory feature of thrombin has fostered the rational engineering of mutants with selectively compromised fibrinogen and PAR1 cleavage. The discovery of the Na(+) effect on thrombin interaction with substrates and the mapping of functional epitopes by Ala scanning mutagenesis have provided a rational and effective strategy for dissociating the procoagulant and anticoagulant activities of the enzyme. Thrombin mutants with selectively compromised activity toward fibrinogen and PAR1 are effective in vivo as anticoagulant and antithrombotic agents.

Thrombin allostery

Thrombin is a Na(+)-activated, allosteric serine protease that plays multiple functional roles in blood pathophysiology. Binding of Na(+) is the major driving force behind the procoagulant, prothrombotic and signaling functions of the enzyme. This review summarizes our current understanding of the molecular basis of thrombin allostery with special emphasis on the kinetic aspects of Na(+) activation. The molecular mechanism of thrombin allostery is a remarkable example of long-range communication that offers a paradigm for many other biological systems.

A control switch for prothrombinase: characterization of a hirudin-like pentapeptide from the COOH terminus of factor Va heavy chain that regulates the rate and pathway for prothrombin activation

Membrane-bound factor Xa alone catalyzes prothrombin activation following initial cleavage at Arg(271) and prethrombin 2 formation (pre2 pathway). Factor Va directs prothrombin activation by factor Xa through the meizothrombin pathway, characterized by initial cleavage at Arg(320) (meizo pathway). We have shown previously that a pentapeptide encompassing amino acid sequence 695-699 from the COOH terminus of the heavy chain of factor Va (Asp-Tyr-Asp-Tyr-Gln, DYDYQ) inhibits prothrombin activation by prothrombinase in a competitive manner with respect to substrate. To understand the mechanism of inhibition of thrombin formation by DYDYQ, we have studied prothrombin activation by gel electrophoresis. Titration of plasma-derived prothrombin activation by prothrombinase, with increasing concentrations of peptide, resulted in complete inhibition of the meizo pathway. However, thrombin formation still occurred through the pre2 pathway. These data demonstrate that the peptide preferentially inhibits initial cleavage of prothrombin by prothrombinase at Arg(320). These findings were corroborated by studying the activation of recombinant mutant prothrombin molecules rMZ-II (R155A/R284A/R271A) and rP2-II (R155A/R284A/R320A) which can be only cleaved at Arg(320) and Arg(271), respectively. Cleavage of rMZ-II by prothrombinase was completely inhibited by low concentrations of DYDYQ, whereas high concentrations of pentapeptide were required to inhibit cleavage of rP2-II. The pentapeptide also interfered with prothrombin cleavage by membrane-bound factor Xa alone in the absence of factor Va increasing the rate for cleavage at Arg(271) of plasma-derived prothrombin or rP2-II. Our data demonstrate that pentapeptide DYDYQ has opposing effects on membrane-bound factor Xa for prothrombin cleavage, depending on the incorporation of factor Va in prothrombinase.

Role of Na+and K+in Enzyme Function

Metal complexation is a key mediator or modifier of enzyme structure and function. In addition to divalent and polyvalent metals, group IA metals Na+and K+play important and specific roles that assist function of biological macromolecules. We examine the diversity of monovalent cation (M+)-activated enzymes by first comparing coordination in small molecules followed by a discussion of theoretical and practical aspects. Select examples of enzymes that utilize M+as a cofactor (type I) or allosteric effector (type II) illustrate the structural basis of activation by Na+and K+, along with unexpected connections with ion transporters. Kinetic expressions are derived for the analysis of type I and type II activation. In conclusion, we address evolutionary implications of Na+binding in the trypsin-like proteases of vertebrate blood coagulation. From this analysis, M+complexation has the potential to be an efficient regulator of enzyme catalysis and stability and offers novel strategies for protein engineering to improve enzyme function.

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.

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.

The Na+ binding channel of human coagulation proteases: novel insights on the structure and allosteric modulation revealed by molecular surface analysis

Thrombovascular diseases result from imbalanced haemostasis and comprise important health problems in the aging population worldwide. The activity of enzymes pertaining to the coagulation cascade of mammalians exhibit several control mechanisms in order to maintain a proper balance between bleeding and thrombosis. For instance, human coagulation serine proteases carrying a F225 or Y225 are allosteric modulated by the binding of Na+ in a water-filled channel connected to the primary specificity pocket (S1 subsite) of these enzymes. We have characterized the structure, topography and lipophilicity of this channel in the ligand-free fast (sodium-bound) and slow (sodium-free) forms of thrombin, in the sole available structure of activated protein C and in several structures of the coagulation factors VIIa, IXa and Xa, differing in the nature of the bound inhibitor and in the occupancy of exosite-I as well as the Ca2+ and Na+ binding sites. Opposite to thrombin, the aqueous channels in all other coagulation enzymes sheltering a Na+ binding site do not have an aperture on the enzyme surface opposite to the S1 subsite entrance. In these enzymes, the lack of the three-residue insertion in loop 1 (183-189) as found in thrombin allied to compensatory mutations in the positions 187-185 and 222 effects a constriction in the water-filled channel that ends up by segregating the ion binding site from the S1 subsite. We also disclosed major topographical changes on the thrombin's surface upon sodium release and transition to the slow form that culminate in the narrowing of the S1 subsite entrance and, strikingly, in the loss of communication between the primary specificity pocket and the exosite-I. Such observation is in accordance with existing experimental data demonstrating thermodynamic linkage between these distant regions on the thrombin surface. Conformational changes in F34, L40, R73 and T74 were the main responsible for this effect. A path by which these changes in the vicinity of exosite-I could be transmitted to the S1 subsite and, consequently, to the sodium binding site is proposed.

Determinants of specificity in coagulation proteases

Proteases play diverse roles in a variety of essential biological processes, both as non-specific catalysts of protein degradation and as highly specific agents that control physiologic events. Here, we review the mechanisms of substrate specificity employed by serine proteases and focus our discussion on coagulation proteases. We dissect the interplay between active site and exosite specificity and how substrate recognition is regulated allosterically by Na+ binding. We also draw attention to a functional polarity that exists in the serine protease fold, which sheds light on the structural linkages between the active site and exosites.

Thrombin functions through its RGD sequence in a non-canonical conformation

Previous studies have suggested that thrombin interacts with integrins in endothelial cells through its RGD (Arg-187, Gly-188, Asp-189) sequence. All existing crystal structures of thrombin show that most of this sequence is buried under the 220-loop and therefore interaction via RGD implies either partial unfolding of the enzyme or its proteolytic digestion. Here, we demonstrate that surface-absorbed thrombin promotes attachment and migration of endothelial cells through interaction with alpha(v)beta(3) and alpha(5)beta(1) integrins. Using site-directed mutants of thrombin we prove that this effect is mediated by the RGD sequence and does not require catalytic activity. The effect is abrogated when residues of the RGD sequence are mutated to Ala and is not observed with proteases like trypsin and tissue-type plasminogen activator, unless the RGD sequence is introduced at position 187-189. The potent inhibitor hirudin does not abrogate the effect, suggesting that thrombin functions through its RGD sequence in a non-canonical conformation. A 1.9-Angstroms resolution crystal structure of free thrombin grown in the presence of high salt (400 mm KCl) shows two molecules in the asymmetric unit, one of which assumes an unprecedented conformation with the autolysis loop shifted 20 Angstroms away from its canonical position, the 220-loop entirely disordered, and the RGD sequence exposed to the solvent.

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.

Energetic and structural consequences of perturbing Gly-193 in the oxyanion hole of serine proteases

The oxyanion hole of serine proteases is formed by the backbone N atoms of the catalytic Ser-195 and Gly-193 and engages the backbone O atom of the P1 residue of substrate in an important H-bonding interaction. The energetic contribution of this interaction in the ground and transition states is presently unknown. Measurements of the individual rate constants defining the catalytic mechanism of substrate hydrolysis for wild-type thrombin and trypsin and their G193A and G193P mutants reveal that Gly-193 is required for optimal substrate binding and acylation. Crystal structures of the G193A and G193P mutants of thrombin bound to the active site inhibitor H-d-Phe-Pro-Arg-CH2Cl document the extent of perturbation induced by the replacement of Gly-193. The Ala mutant weakens the H-bonding interaction of the N atom of residue 193, whereas the Pro substitution abrogates it altogether with additional small shifts of the protein backbone. From the kinetic and structural data, we estimate that the H-bonding interaction in the oxyanion hole contributes a stabilization of the ground and transition states of > 1.5 kcal/mol but < 3.0 kcal/mol. These results shed light on a basic aspect of the enzyme-substrate interaction in the entire family of trypsin-like serine proteases.

A critical role for Gly25 in the B chain of human thrombin

We have recently identified (Akhavan S et al., Thromb Haemost 2000; 84: 989-97) a patient with a mild bleeding diathesis associated to an homozygous mutation in the thrombin B chain (Gly25Ser, chymotrypsinogen numbering, i.e. position 330 in human prothrombin numbering). Transient transfection of wild-type prothrombin (FII-WT) and mutant prothrombin (designated FII-G25(330)S) cDNA in COS-7 cells showed a mild reduction (50%) in FII-G25(330)S production. Recombinant proteins, stably expressed in Chinese hamster ovary cells, were isolated and activated by Taipan snake or Echis carinatus venoms. We show that the G25(330)S mutation results in a decrease in the rate of prothrombin proteolytic activation. The mutation also significantly decreases (i) the catalytic activity of thrombin with a 9-fold reduction in catalytic efficiency of the mutant toward S-2238; (ii) the interaction with benzamidine; (iii) the rate of inhibition by TLCK and antithrombin; and (iv) the rate of hydrolysis of macromolecular substrates (fibrinogen, protein C). In contrast, exosite I does not appear to be affected by the molecular defect. These results, together with molecular modeling and dynamics, indicate that Gly25(330) is important for proper expression and probably proper folding of prothrombin, and also plays a critical role in both the alignment of the catalytic triad and the flexibility of one of the activation segments of prothrombin.

Thrombomodulin changes the molecular surface of interaction and the rate of complex formation between thrombin and protein C

The interaction of thrombin with protein C triggers a key down-regulatory process of the coagulation cascade. Using a panel of 77 Ala mutants, we have mapped the epitope of thrombin recognizing protein C in the absence or presence of the cofactor thrombomodulin. Residues around the Na(+) site (Thr-172, Lys-224, Tyr-225, and Gly-226), the aryl binding site (Tyr-60a), the primary specificity pocket (Asp-189), and the oxyanion hole (Gly-193) hold most of the favorable contributions to protein C recognition by thrombin, whereas a patch of residues in the 30-loop (Arg-35 and Pro-37) and 60-loop (Phe-60h) regions produces unfavorable contributions to binding. The shape of the epitope changes drastically in the presence of thrombomodulin. The unfavorable contributions to binding disappear and the number of residues promoting the thrombin-protein C interaction is reduced to Tyr-60a and Asp-189. Kinetic studies of protein C activation as a function of temperature reveal that thrombomodulin increases >1,000-fold the rate of diffusion of protein C into the thrombin active site and lowers the activation barrier for this process by 4 kcal/mol. We propose that the mechanism of thrombomodulin action is to kinetically facilitate the productive encounter of thrombin and protein C and to allosterically change the conformation of the activation peptide of protein C for optimal presentation to the thrombin active site.

Molecular Mapping of the Thrombin-Heparin Cofactor II Complex

We used 55 Ala-scanned recombinant thrombin molecules to define residues important for inhibition by the serine protease inhibitor (serpin) heparin cofactor II (HCII) in the absence and presence of glycosaminoglycans. We verified the importance of numerous basic residues in anion-binding exosite-1 (exosite-1) and found 4 additional residues, Gln24, Lys65, His66, and Tyr71 (using the thrombin numbering system), that were resistant to HCII inhibition with and without glycosaminoglycans. Inhibition rate constants for these exosite-1 (Q24A, K65A, H66A, Y71A) thrombin mutants (0.02-0.38 x 10(8) m(-1) min(-1) for HCII-heparin when compared with 2.36 x 10(8) m(-1) min(-1) with wild-type thrombin and 0.03-0.53 x 10(8) m(-1) min(-1) for HCII-dermatan sulfate when compared with 5.23 x 10(8) m(-1) min(-1) with wild-type thrombin) confirmed that the structural integrity of thrombin exosite-1 is critical for optimal HCII-thrombin interactions in the presence of glycosaminoglycans. However, our results are also consistent for HCII-glycosaminoglycan-thrombin ternary complex formation. Ten residues surrounding the active site of thrombin were implicated in HCII interactions. Four mutants (Asp51, Lys52, Lys145/Thr147/Trp148, Asp234) showed normal increased rates of inhibition by HCII-glycosaminoglycans, whereas four mutants (Trp50, Glu202, Glu229, Arg233) remained resistant to inhibition by HCII with glycosaminoglycans. Using 11 exosite-2 thrombin mutants with 20 different mutated residues, we saw no major perturbations of HCII-glycosaminoglycan inhibition reactions. Collectively, our results support a "double bridge" mechanism for HCII inhibition of thrombin in the presence of glycosaminoglycans, which relies in part on ternary complex formation but is primarily dominated by an allosteric process involving contact of the "hirudin-like" domain of HCII with thrombin exosite-1.

Zymogenic and enzymatic properties of the 70-80 loop mutants of factor X/Xa

The Ca(2+) binding 70-80 loop of factor X (fX) contains one basic (Arg(71)) and three acidic (Glu(74), Glu(76), and Glu(77)) residues whose contributions to the zymogenic and enzymatic properties of the protein have not been evaluated. We prepared four Ala substitution mutants of fX (R71A, E74A, E76A, and E77A) and characterized their activation kinetics by the factor VIIa and factor IXa in both the absence and presence of cofactors. Factor VIIa exhibited normal activity toward E74A and E76A and less than a twofold impaired activity toward R71A and E77A in both the absence and presence of tissue factor. Similarly, factor IXa in the absence of factor VIIIa exhibited normal activity toward both E74A and E76A; however, its activity toward R71A and E77A was impaired approximately two- to threefold. In the presence of factor VIIIa, factor IX activated all mutants with approximately two- to fivefold impaired catalytic efficiency. In contrast to changes in their zymogenic properties, all mutant enzymes exhibited normal affinities for factor Va, and catalyzed the conversion of prothrombin to thrombin with normal catalytic efficiencies. However, further studies revealed that the affinity of mutant enzymes for interaction with metal ions Na(+) and Ca(2+) was impaired. These results suggest that although charged residues of the 70-80 loop play an insignificant role in fX recognition by the factor VIIa-tissue factor complex, they are critical for the substrate recognition by factor IXa in the intrinsic Xase complex. The results further suggest that mutant residues do not play a specific role in the catalytic function of fXa in the prothrombinase complex.

Molecular dissection of Na+ binding to thrombin

Na(+) binding near the primary specificity pocket of thrombin promotes the procoagulant, prothrombotic, and signaling functions of the enzyme. The effect is mediated allosterically by a communication between the Na(+) site and regions involved in substrate recognition. Using a panel of 78 Ala mutants of thrombin, we have mapped the allosteric core of residues that are energetically linked to Na(+) binding. These residues are Asp-189, Glu-217, Asp-222, and Tyr-225, all in close proximity to the bound Na(+). Among these residues, Asp-189 shares with Asp-221 the important function of transducing Na(+) binding into enhanced catalytic activity. None of the residues of exosite I, exosite II, or the 60-loop plays a significant role in Na(+) binding and allosteric transduction. X-ray crystal structures of the Na(+)-free (slow) and Na(+)-bound (fast) forms of thrombin, free or bound to the active site inhibitor H-d-Phe-Pro-Arg-chloromethyl-ketone, document the conformational changes induced by Na(+) binding. The slow --> fast transition results in formation of the Arg-187:Asp-222 ion pair, optimal orientation of Asp-189 and Ser-195 for substrate binding, and a significant shift of the side chain of Glu-192 linked to a rearrangement of the network of water molecules that connect the bound Na(+) to Ser-195 in the active site. The changes in the water network and the allosteric core explain the thermodynamic signatures linked to Na(+) binding and the mechanism of thrombin activation by Na(+). The role of the water network uncovered in this study establishes a new paradigm for the allosteric regulation of thrombin and other Na(+)-activated enzymes involved in blood coagulation and the immune response.

The Contribution of Amino Acid Region Asp695-Tyr698 of Factor V to Procofactor Activation and Factor Va Function

There is strong evidence that a functionally important cluster of amino acids is located on the COOH-terminal portion of the heavy chain of factor Va, between amino acid residues 680 and 709. To ascertain the importance of this region for cofactor activity, we have synthesized five overlapping peptides representing this amino acid stretch (10 amino acids each, HC1-HC5) and tested them for inhibition of prothrombinase assembly and function. Two peptides, HC3 (spanning amino acid region 690-699) and HC4 (containing amino acid residues 695-704), were found to be potent inhibitors of prothrombinase activity with IC(50) values of approximately 12 and approximately 10 microm, respectively. The two peptides were unable to interfere with the binding of factor Va to active site fluorescently labeled Glu-Gly-Arg human factor Xa, and kinetic analyses showed that HC3 and HC4 are competitive inhibitors of prothrombinase with respect to prothrombin with K(i) values of approximately 6.3 and approximately 5.3 microm, respectively. These data suggest that the peptides inhibit prothrombinase because they interfere with the incorporation of prothrombin into prothrombinase. The shared amino acid motif between HC3 and HC4 is composed of Asp(695)-Tyr-Asp-Tyr-Gln(699) (DYDYQ). A pentapeptide with this sequence inhibited both prothrombinase function with an IC(50) of 1.6 microm (with a K(D) for prothrombin of 850 nm), and activation of factor V by thrombin. Peptides HC3, HC4, and DYDYQ were also found to interact with immobilized thrombin. A recombinant factor V molecule with the mutations Asp(695) --> Lys, Tyr(696) --> Phe, Asp(697) --> Lys, and Tyr(698) --> Phe (factor V(2K2F)) was partially resistant to activation by thrombin but could be readily activated by RVV-V activator (factor Va(RVV)(2K2F)) and factor Xa (factor Va(Xa)(2K2F)). Factor Va(RVV)(2K2F) and factor Va(Xa)(2K2F) had impaired cofactor activity within prothrombinase in a system using purified reagents. Our data demonstrate for the first time that amino acid sequence 695-698 of factor Va heavy chain is important for procofactor activation and is required for optimum prothrombinase function. These data provide functional evidence for an essential and productive contribution of factor Va to the activity of prothrombinase.

Thrombomodulin allosterically modulates the activity of the anticoagulant thrombin

Exosite 1 of thrombin consists of a cluster of basic residues (Arg-35, Lys-36, Arg-67, Lys-70, Arg-73, Arg-75, and Arg-77 in chymotrypsinogen numbering) that play key roles in the function of thrombin. Structural data suggest that the side chain of Arg-35 projects toward the active site pocket of thrombin, but all other residues are poised to interact with thrombomodulin (TM). To study the role of these residues in TM-mediated protein C (PC) activation by thrombin, a charge reversal mutagenesis approach was used to replace these residues with a Glu in separate constructs. The catalytic properties of the mutants toward PC were analyzed in both the absence and presence of TM and Ca2+. It was discovered that, with the exception of the Arg-67 and Lys-70 mutants, all other mutants activated PC with similar maximum rate constants in the presence of a saturating concentration of TM and Ca2+, although their affinity for interaction with TM was markedly impaired. The catalytic properties of the Arg-35 mutant were changed so that PC activation by the mutant no longer required Ca2+ in the presence of TM, but, instead, it was accelerated by EDTA. Moreover, the activity of this mutant toward PC was improved approximately 25-fold independent of TM. These results suggest that Arg-35 is responsible for the Ca2+ dependence of PC activation by the thrombin-TM complex and that a function for TM in the activation complex is the allosteric alleviation of the inhibitory interaction of Arg-35 with the substrate.

Thrombomodulin enhances the reactivity of thrombin with protein C inhibitor by providing both a binding site for the serpin and allosterically modulating the activity of thrombin

Thrombomodulin (TM), or its epidermal growth factor-like domains 456 (TM456), enhances the catalytic efficiency of thrombin toward both protein C and protein C inhibitor (PCI) by 2-3 orders of magnitude. Structural and mutagenesis data have indicated that the interaction of basic residues of the heparin-binding exosite of protein C with the acidic residues of TM4 is partially responsible for the efficient activation of the substrate by the thrombin-TM456 complex. Similar to protein C, PCI has a basic exosite (H-helix) that constitutes the heparin-binding site of the serpin. To determine whether TM accelerates the reactivity of thrombin with PCI by providing a binding site for the H-helix of the serpin, an antithrombin (AT) mutant was constructed in which the H-helix of the serpin was replaced with the same region of PCI (AT-PCIH-helix). Unlike PCI, the H-helix of AT is negatively charged. It was discovered that TM456 slightly (<2-fold) impaired the reactivity of AT with thrombin; however, it enhanced the reactivity of AT-PCIH-helix with the protease by an order of magnitude. Further studies revealed that the substitution of Arg35 of thrombin with an Ala also resulted in an order of magnitude enhancement in reactivity of the protease with both PCI and AT-PCIH-helix independent of TM. We conclude that TM enhances the reactivity of PCI with thrombin by providing both a binding site for the serpin and a conformational modulation of the extended binding pocket of thrombin.

Thrombin interactions

After generation from prothrombin, thrombin plays multiple roles in the blood coagulation cascade that are mediated by interaction with a number of physiologic substrates, effectors, and inhibitors. Structural and mutagenesis studies have helped unravel the molecular basis of thrombin interactions in the context of both well-established and emerging new roles of the enzyme. The functional versatility of thrombin owes much to its evolutionary origin and results from structural determinants and mechanisms that can be exploited by pharmacologic intervention.

Structural requirements for expression of factor Va activity

Thrombin activated factor Va (factor VIIa, residues 1-709 and 1546-2196) has an apparent dissociation constant (Kd,app) for factor Xa within prothrombinase of approximately 0.5 nM. A protease (NN) purified from the venom of the snake Naja nigricollis nigricollis, cleaves human factor V at Asp697, Asp1509, and Asp1514 to produce a molecule (factor VNN) that is composed of a Mr 100,000 heavy chain (amino acid residues 1-696) and a Mr 80,000 light chain (amino acid residues 1509/1514-2196). Factor VNN, has a Kd,app for factor Xa of 4 nm and reduced clotting activity. Cleavage of factor VIIa by NN at Asp697 results in a cofactor that loses approximately 60-80% of its clotting activity. An enzyme from Russell's viper venom (RVV) cleaves human factor V at Arg1018 and Arg1545 to produce a Mr 150,000 heavy chain and Mr 74,000 light chain (factor VRVV, residues 1-1018 and 1546-2196). The RVV species has affinity for factor Xa and clotting activity similar to the thrombin-activated factor Va. Cleavage of factor VNN at Arg1545 by alpha-thrombin (factor VNN/IIa) or RVV (factor VNN/RVV) leads to enhanced affinity of the cofactor for factor Xa (Kd,app approximately 0.5 nM). A synthetic peptide containing the last 13 residues from the heavy chain of factor Va (amino acid sequence 697-709, D13R) was found to be a competitive inhibitor of prothrombinase with respect to prothrombin. The peptide was also found to specifically interact with thrombin-agarose. These data demonstrate that 1) cleavage at Arg1545 and formation of the light chain of factor VIIa is essential for high affinity binding and function of factor Xa within prothrombinase and 2) a binding site for prothrombin is contributed by amino acid residues 697-709 of the heavy chain of the cofactor.

A steady-state competition model describes the modulating effects of thrombomodulin on thrombin inhibition by plasminogen activator inhibitor-1 in the absence and presence of vitronectin

Thrombomodulin (TM) slows down the interaction rate between thrombin and plasminogen activator inhibitor 1 (PAI-1). We now show that the 12-fold reduced inhibition rate in the presence of TM does not result from an altered distribution between PAI-1 cleavage and irreversible complex formation. Surface plasmon resonance (SPR) revealed an over 200-fold reduced affinity of TM for thrombin-VR1tPA as compared to thrombin, demonstrating the importance of the VR1 loop in the interaction of thrombin with both TM and PAI-1. Furthermore, in contrast to ATIII, PAI-1 was not able to bind the thrombin/TM complex demonstrating complete competitive binding between PAI-1 and TM. Kinetic modeling on the inhibitory effect of TM confirms a mechanism that involves complete steric blocking of the thrombin/PAI-1 interaction. Also, it accurately decribes the biphasic inhibition profile resulting from the substantial reduction of the extremely fast rate of reversible Michaelis complex formation, which is essential for efficient inhibition of thrombin by PAI-1. Vitronectin (VN) is shown to partially relieve TM inhibitory action only by vastly increasing the initial rate of interaction between free thrombin and PAI-1. In addition, SPR established that solution-phase PAI-1/VN complexes and non-native VN (extracellular matrix form) bind TM directly via the chondroitin sulphate moiety of TM. Collectively, these results show that VR1 is a subsite of exosite 1 on thrombin's surface, which regulates exclusive binding of either PAI-1 or TM. This competition will be physiologically significant in controlling the mitogenic activity of thrombin during vascular disease.

A novel class of small functional peptides that bind and inhibit human alpha-thrombin isolated by mRNA display

Here we report the in vitro selection of novel small peptide motifs that bind to human alpha-thrombin. We have applied mRNA display to select for thrombin binding peptides from an unbiased library of 1.2 x 10(11) different 35-mer peptides, each containing a random sequence of 15 amino acids. Two clones showed binding affinities ranging from 166 to 520 nM. A conserved motif of four amino acids, DPGR, was identified. Clot formation of human plasma is inhibited by the selected clones, and they downregulate the thrombin-mediated activation of protein C. The identified peptide motifs do not share primary sequence similarities to any of the known natural thrombin binding motifs. As new inhibitors for human thrombin open interesting possibilities in thrombosis research, our newly identified peptides may provide further insights into this field of investigation and may be possible candidates for the development of new anti-thrombotic agents.

Crystal structure of the anticoagulant slow form of thrombin

Using the thrombin mutant R77aA devoid of the site of autoproteolytic degradation at exosite I, we have solved for the first time the structure of thrombin free of any inhibitors and effector molecules and stabilized in the Na(+)-free slow form. The slow form shows subtle differences compared with the currently available structures of the Na(+)-bound fast form that carry inhibitors at the active site or exosite I. The most notable differences are the displacement of Asp-189 in the S1 specificity pocket, a downward shift of the 190-193 strand, a rearrangement of the side chain of Glu-192, and a significant shift in the position of the catalytic Ser-195 that is no longer within H-bonding distance from His-57. The structure of the slow form explains the reduced specificity toward synthetic and natural substrates and suggests a molecular basis for its anticoagulant properties.