The Na+ binding site of thrombin

Role of residue Phe225 in the cofactor-mediated, allosteric regulation of the serine protease coagulation factor VIIa

Functional regulation by cofactors is fundamentally important for the highly ordered, consecutive activation of the coagulation cascade. The initiating protease of the coagulation system, factor VIIa (VIIa), retains zymogen-like features after proteolytic cleavage of the activating Arg(15)-Ile(16) peptide bond and requires the binding of the cofactor tissue factor (TF) to stabilize the protease domain in an active enzyme conformation. Structural comparison of TF-bound and free VIIa failed to provide a conclusive mechanism for this catalytic activation. This study provides novel insight into the cofactor-dependent regulation of VIIa by demonstrating that the side chain of Phe(225), an aromatic residue that is common to allosterically regulated serine proteases, is necessary for optimal TF-mediated activation of VIIa's catalytic function. However, mutation of Phe(225) did not abolish the cofactor-induced stabilization of the Ile(16)-Asp(194) salt bridge, previously considered the primary switch mechanism for activating VIIa. Moreover, mutation of other residue side chains in the VIIa protease domain resulted in a reduced level of or no stabilization of the amino-terminal insertion site upon TF binding, with little or no effect on the TF-mediated enhancement of catalysis. This study thus establishes a crucial role for the aromatic Phe(225) residue position in the allosteric network that transmits the activating switch from the cofactor interface to the catalytic cleft, providing insight into the highly specific conformational linkages that regulate serine protease function.

Interaction of the factor XIII activation peptide with alpha -thrombin. Crystal structure of its enzyme-substrate analog complex

The serine protease thrombin proteolytically activates blood coagulation factor XIII by cleavage at residue Arg(37); factor XIII in turn cross-links fibrin molecules and gives mechanical stability to the blood clot. The 2.0-A resolution x-ray crystal structure of human alpha-thrombin bound to the factor XIII-(28-37) decapeptide has been determined. This structure reveals the detailed atomic level interactions between the factor XIII activation peptide and thrombin and provides the first high resolution view of this functionally important part of the factor XIII molecule. A comparison of this structure with the crystal structure of fibrinopeptide A complexed with thrombin highlights several important determinants of thrombin substrate interaction. First, the P1 and P2 residues must be compatible with the geometry and chemistry of the S1 and S2 specificity sites in thrombin. Second, a glycine in the P5 position is necessary for the conserved substrate conformation seen in both factor XIII-(28-37) and fibrinopeptide A. Finally, the hydrophobic residues, which occupy the aryl binding site of thrombin determine the substrate conformation further away from the catalytic residues. In the case of factor XIII-(28-37), the aryl binding site is shared by hydrophobic residues P4 (Val(34)) and P9 (Val(29)). A bulkier residue in either of these sites may alter the substrate peptide conformation.

The preparation and catalytic properties of recombinant human prostate-specific antigen (rPSA)

The serine proteinase prostate-specific antigen (PSA), and its complex with the serine proteinase inhibitor alpha(1)-antichymotrypsin (ACT), have been used as markers for the diagnosis of prostate cancer. PSA prepared from seminal fluid is typically contaminated with the trypsin-like glandular kallikrein (hK2). Here we describe a convenient and reproducible preparation of catalytically active recombinant PSA (rPSA) and demonstrate an overall similarity in the properties of cloned and refolded rPSA to PSA purified from seminal fluid. We also present results that are relevant for increasing the sensitivity of assays of PSA activity in biological fluids, for the putative role of PSA activity in physiologically important processes, including prostate cancer metastasis, and for the design of PSA inhibitors. Specifically, we find that added salts, in particular NaCl, give rise to dramatic increases in rPSA catalytic activity, as does added glycerol. On the other hand, Zn(2+), spermine, and spermidine, each a major component of seminal and prostatic fluid, strongly inhibit rPSA activity, with Zn(2+) being a non-competitive inhibitor while spermine is a competitive inhibitor. Citrate, also a major component of seminal and prostatic fluid, spermine, and spermidine each protect rPSA from Zn(2+) inhibition, presumably via Zn(2+) sequestration. Finally, rPSA efficiently proteolyzes several protein substrates.

Mutation of W215 compromises thrombin cleavage of fibrinogen, but not of PAR-1 or protein C

W215 is a highly conserved residue that shapes the S3 and S4 specificity sites of thrombin and participates in an edge-to-face interaction with residue F8 of the fibrinogen Aalpha chain. Protein C and the platelet receptor PAR-1 carry an acidic residue at P3 and bind to the active site of thrombin without making contact with W215. This suggested that mutation of W215 could dissociate the cleavage of fibrinogen from that of protein C and PAR-1. Replacement of W215 with Phe produces modest effects on thrombin function, whereas the W215Y replacement compromises significantly the catalytic activity toward all chromogenic and natural substrates that are tested. Replacement of W215 with Ala almost obliterates Na(+) binding, reduces the level of fibrinogen cleavage 500-fold, but decreases the levels of protein C activation and PAR-1 cleavage only 3- and 25-fold, respectively. The W215A mutant cleaves PAR-1 with a specificity constant that is more than 13-fold higher than that of fibrinogen and protein C and is the first thrombin derivative to be described that functions as an almost exclusive activator of PAR-1. The environment of W215 influences differentially three physiologically important interactions of thrombin, which should assist in the study of each of these functions separately in vivo.

Synthetic receptors as models for alkali metal cation-pi binding sites in proteins

The alkali metal cations Na(+) and K(+) have several important physiological roles, including modulating enzyme activity. Recent work has suggested that alkali metal cations may be coordinated by pi systems, such as the aromatic amino acid side chains. The ability of K(+) to interact with an aromatic ring has been assessed by preparing a family of synthetic receptors that incorporate the aromatic side chains of phenylalanine, tyrosine, and tryptophan. These receptors are constructed around a diaza-18-crown-6 scaffold, which serves as the primary binding site for an alkali metal cation. The ability of the aromatic rings to coordinate a cation was determined by crystallizing each of the receptors in the presence of K(+) and by solving the solid state structures. In all cases, complexation of K(+) by the pi system was observed. When possible, the structures of the unbound receptors also were determined for comparison. Further proof that the aromatic ring makes an energetically favorable interaction with the cation was obtained by preparing a receptor in which the arene was perfluorinated. Fluorination of the arene reverses the electrostatics, but the aromaticity is maintained. The fluorinated arene rings do not coordinate the cation in the solid state structure of the K(+) complex. Thus, the results of the predicted electrostatic reversal were confirmed. Finally, the biological implications of the alkali metal cation-pi interaction are addressed.

Identification of a binding site for quaternary amines in factor Xa

In the process of characterizing the Na(+)-binding properties of factor Xa, a specific inhibition of this enzyme by quaternary amines was identified, consistent with previous observations. The binding occurs with K(i) in the low millimolar range, with trimethylphenylammonium (TMPA) showing the highest specificity. Binding of TMPA inhibits substrate hydrolysis in a competitive manner, does not inhibit the binding of p-aminobenzamidine to the S1 pocket, and is positively linked to Na(+) binding. Inhibition by TMPA is also seen in thrombin and tissue plasminogen activator (tPA), though to a lesser extent compared to factor Xa. Computer modeling using the crystal structure of factor Xa suggests that TMPA binds to the S2/S3 specificity sites, with its hydrophobic moiety making van der Waals interactions with the side chains of Y99, F174, and W215, and the charged amine coupling electrostatically with the carboxylates of E97. Site-directed mutagenesis of factor Xa, thrombin, and tPA confirms the predictions drawn by docking calculations and reveal a dominant role for residue Y99. Binding of TMPA to factor Xa is drastically (25-fold) reduced by the Y99T replacement. Likewise, the Y99L substitution compromises binding of TMPA to tPA. On the other hand, the affinity of TMPA is enhanced 4-fold in thrombin with the substitution L99Y. The identification of a binding site for quaternary amines in factor Xa has a bearing on the rational design of selective inhibitors of this clotting enzyme.

Thrombin inhibitor design: X-ray and solution studies provide a novel P1 determinant

The crystal structures of proflavin and 6-fluorotryptamine thrombin have been completed showing binding of both ligands at the active site S1 pocket. The structure of proflavin:thrombin was confirmatory, while the structure of 6-fluorotryptamine indicated a novel binding mode at the thrombin active site. Furthermore, speculation that the sodium atom identified in an extended solvent channel beneath the S pocket may play a role in binding of these ligands was investigated by direct proflavin titrations as well as chromogenic activity measurements as a function of sodium concentration at constant ionic strength. These results suggested a linkage between the sodium site and the S1 pocket. This observation could be due to a simple ionic interaction between Asp189 and the sodium ion or a more complicated structural rearrangement of the thrombin S1 pocket. Finally, the unique binding mode of 6-fluorotryptamine provides ideas toward the design of a neutrally charged thrombin inhibitor.

Thermal stability and atomic-resolution crystal structure of the Bacillus caldolyticus cold shock protein

The bacterial cold shock proteins are small compact beta-barrel proteins without disulfide bonds, cis-proline residues or tightly bound cofactors. Bc-Csp, the cold shock protein from the thermophile Bacillus caldolyticus shows a twofold increase in the free energy of stabilization relative to its homolog Bs-CspB from the mesophile Bacillus subtilis, although the two proteins differ by only 12 out of 67 amino acid residues. This pair of cold shock proteins thus represents a good system to study the atomic determinants of protein thermostability. Bs-CspB and Bc-Csp both unfold reversibly in cooperative transitions with T(M) values of 49.0 degrees C and 77.3 degrees C, respectively, at pH 7.0. Addition of 0.5 M salt stabilizes Bs-CspB but destabilizes Bc-Csp. To understand these differences at the structural level, the crystal structure of Bc-Csp was determined at 1.17 A resolution and refined to R=12.5% (R(free)=17.9%). The molecular structures of Bc-Csp and Bs-CspB are virtually identical in the central beta-sheet and in the binding region for nucleic acids. Significant differences are found in the distribution of surface charges including a sodium ion binding site present in Bc-Csp, which was not observed in the crystal structure of the Bs-CspB. Electrostatic interactions are overall favorable for Bc-Csp, but unfavorable for Bs-CspB. They provide the major source for the increased thermostability of Bc-Csp. This can be explained based on the atomic-resolution crystal structure of Bc-Csp. It identifies a number of potentially stabilizing ionic interactions including a cation-binding site and reveals significant changes in the electrostatic surface potential.

Sodium binding site of factor Xa: role of sodium in the prothrombinase complex

The nature of residue 225 on a consensus loop in serine proteases determines whether a protease can bind Na(+). Serine proteases with a Pro at this position are unable to bind Na(+), but those with a Tyr or Phe can bind Na(+). Factor Xa (FXa), the serine protease of the prothrombinase complex, contains a Tyr at this position. Na(+) is also known to stimulate the amidolytic activity of FXa toward cleavage of small synthetic substrates, but the role of Na(+) in the prothrombinase complex has not been investigated. In this study, we engineered a Gla-domainless form of FX (GDFX) in which residue Tyr(225) was replaced with a Pro. We found that Na(+) stimulated the cleavage rate of chromogenic substrates by FXa or GDFXa approximately 8-24-fold with apparent dissociation constants [K(d(app))] of 37 and 182 mM in the presence and absence of Ca(2+), respectively. In contrast, Na(+) minimally affected the cleavage rate of these substrates by the mutant, and no K(d(app)) for Na(+) binding to the mutant could be estimated. Unlike the wild-type enzyme, the reactivity of the mutant with antithrombin was independent of Na(+) and impaired approximately 32-fold. Ca(2+) improved the reactivity of the mutant with antithrombin approximately 5-fold. Affinity of the mutant for binding to factor Va was weakened and its ability to activate prothrombin was severely impaired. Further studies with the wild-type prothrombinase complex revealed that FXa binds to factor Va with a similar K(d(app)) of 1. 1-1.8 nM in the presence of Na(+), K(+), Li(+), Ch(+), and Tris(+) and that the catalytic efficiency of prothrombinase is enhanced less than 1.5-fold by the specific effect of Na(+) in the reaction buffer. These results suggest that (1) the loop including residue 225 (225-loop) is a Na(+) binding site in FXa, (2) the Na(+)- and Ca(2+)-binding loops of FXa are allosterically linked, and (3) the Tyr conformer of the 225-loop is critical for factor Xa function; however, both Na(+)-bound and Na(+)-free forms of factor Xa in the prothrombinase complex can efficiently activate prothrombin.

The role of Glu(192) in the allosteric control of the S(2)' and S(3)' subsites of thrombin

Thrombin is an allosteric protease controlled through exosites flanking the catalytic groove. Binding of a peptide derived from hirudin (Hir(52-65)) and/or of heparin to these opposing exosites alters catalysis. We have investigated the contribution of subsites S(2)' and S(3)' to this allosteric transition by comparing the hydrolysis of two sets of fluorescence-quenched substrates having all natural amino acids at positions P(2)' and P(3)'. Regardless of the amino acids, Hir(52-65) decreased, and heparin increased the k(cat)/K(m) value of hydrolysis by thrombin. Several lines of evidence have suggested that Glu(192) participates in this modulation. We have examined the role of Glu(192) by comparing the catalytic activity of thrombin and its E192Q mutant. Mutation substantially diminishes the selectivity of thrombin. The substrate with the "best" P(2)' residue was cleaved with a k(cat)/K(m) value only 49 times higher than the one having the "least favorable" P(2)' residue (versus 636-fold with thrombin). Mutant E192Q also lost the strong preference of thrombin for positively charged P(3)' residues and its strong aversion for negatively charged P(3)' residues. Furthermore, both Hir(52-65) and heparin increased the k(cat)/K(m) value of substrate hydrolysis. We conclude that Glu(192) is critical for the P(2)' and P(3)' specificities of thrombin and for the allostery mediated through exosite 1.

Thrombin interacts with thrombomodulin, protein C, and thrombin-activatable fibrinolysis inhibitor via specific and distinct domains

A collection of 56 purified thrombin mutants, in which 76 charged or polar surface residues on thrombin were mutated to alanine, was used to identify key residues mediating the interactions of thrombin with thrombomodulin (TM), protein C, and thrombin-activatable fibrinolysis inhibitor (TAFI). Comparison of protein C activation in the presence and absence of TM identified 11 residues mediating the thrombin-TM interaction (Lys(21), Gln(24), Arg(62), Lys(65), His(66), Arg(68), Thr(69), Tyr(71), Arg(73), Lys(77), Lys(106)). Three mutants (E25A, D51A, R89A/R93A/E94A) were found to have decreased ability to activate TAFI yet retained normal protein C activation, whereas three other mutants (R178A/R180A/D183A, E229A, R233A) had decreased ability to activate protein C but maintained normal TAFI activation. One mutant (W50A) displayed decreased activation of both substrates. Mapping of these functional residues on thrombin revealed that the 11 residues mediating the thrombin-TM interaction are all located in exosite I. Residues important in TAFI activation are located above the active-site cleft, whereas residues involved in protein C are located below the active-site cleft. In contrast to the extensive overlap of residues mediating TM binding and fibrinogen clotting, these data show that distinct domains in thrombin mediate its interactions with TM, protein C, and TAFI. These studies demonstrate that selective enzymatic properties of thrombin can be dissociated by site-directed mutagenesis.

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.

Alkali metal ions protect mitochondrial rhodanese against thermal inactivation

Incubation of bovine liver mitochondrial rhodanese in dilute, reducing solutions at temperatures ranging between 30 and 45 degreesC conduced to a rapid loss of enzymatic activity. This inactivation was substantially reduced in the presence of millimolar concentrations of alkali metal ions, divalent cations (including Mg2+, Ca2+, and Ba2+) were ineffective. The extent of protection afforded by monovalent cations was highly dependent on their ionic radii, with K+ and Na+ ions being the most effective protective agents. The protection afforded by a number of anions, including thiosulfate, could be totally ascribed to the presence of the accompanying monovalent cation. The overall results indicate that K+ and Na+, at concentrations and temperatures within the physiological range, substantially contribute to the stabilization of the functional structure of rhodanese.

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.

Anticoagulant Thrombins

Thrombin plays both procoagulant and anticoagulant roles in the blood coagulation cascade. This dual role is influenced allosterically by the binding of Na(+) near the primary specificity pocket of the enzyme. Recent findings demonstrate that it is possible to engineer recombinant thrombins that have practically lost anticoagulant activity but retain their anticoagulant properties. These anticoagulant thrombins bear substitutions in or around the Na(+) binding site, provide important clues on structure-function relations, and offer a promising alternative to current anticoagulant therapies. In addition, they demonstrate the importance of Na(+) as a coagulation factor and broaden our understanding of the function and regulation of all vitamin K-dependent clotting enzymes.

Allosteric effects of a monoclonal antibody against thrombin exosite II

We previously isolated a monoclonal antithrombin IgG from a patient with multiple myeloma [Colwell et al. (1997) Br. J. Haematol. 97, 219-226]. Using a panel of 55 surface mutants of recombinant thrombin, we now show that the epitope for the IgG most likely includes Arg-101, Arg-233, and Lys-236 in exosite II. The IgG affects the rate at which thrombin cleaves various peptide p-nitroanilide substrates with arginine in the P1 position, increasing the kcat for substrates with a P2 glycine residue but generally decreasing the kcat for substrates with a P2 proline. The allosteric effect of the IgG is altered by deletion of Pro-60b, Pro-60c, and Trp-60d from the 60-loop of thrombin, which lies between exosite II and the catalytic triad. The effect of the IgG, however, does not depend on the presence or absence of sodium ions, a known allosteric regulator of thrombin. The IgG does not affect the conformation of thrombin exosite I as determined by binding of a fluorescent derivative of hirudin54-65. These results provide evidence for a direct allosteric linkage between exosite II and the catalytic site of thrombin.

Cluster analysis of consensus water sites in thrombin and trypsin shows conservation between serine proteases and contributions to ligand specificity

Cluster analysis is presented as a technique for analyzing the conservation and chemistry of water sites from independent protein structures, and applied to thrombin, trypsin, and bovine pancreatic trypsin inhibitor (BPTI) to locate shared water sites, as well as those contributing to specificity. When several protein structures are superimposed, complete linkage cluster analysis provides an objective technique for resolving the continuum of overlaps between water sites into a set of maximally dense microclusters of overlapping water molecules, and also avoids reliance on any one structure as a reference. Water sites were clustered for ten superimposed thrombin structures, three trypsin structures, and four BPTI structures. For thrombin, 19% of the 708 microclusters, representing unique water sites, contained water molecules from at least half of the structures, and 4% contained waters from all 10. For trypsin, 77% of the 106 microclusters contained water sites from at least half of the structures, and 57% contained waters from all three. Water site conservation correlated with several environmental features: highly conserved microclusters generally had more protein atom neighbors, were in a more hydrophilic environment, made more hydrogen bonds to the protein, and were less mobile. There were significant overlaps between thrombin and trypsin conserved water sites, which did not localize to their similar active sites, but were concentrated in buried regions including the solvent channel surrounding the Na+ site in thrombin, which is associated with ligand selectivity. Cluster analysis also identified water sites conserved in thrombin but not trypsin, and vice versa, providing a list of water sites that may contribute to ligand discrimination. Thus, in addition to facilitating the analysis of water sites from multiple structures, cluster analysis provides a useful tool for distinguishing between conserved features within a protein family and those conferring specificity.

Synthesis and characterization of more potent analogues of hirudin fragment 1-47 containing non-natural amino acids

Hirudin is the most potent and specific inhibitor of thrombin, a key enzyme in the coagulation process existing in equilibrium between its procoagulant (fast) and anticoagulant (slow) form. In a previous study, we described the solid-phase synthesis of a Trp3 analogue of fragment 1-47 of hirudin HM2, which displayed approximately 5-fold higher thrombin inhibitory potency relative to that of the natural product [De Filippis, V., et al. (1995) Biochemistry 34, 9552-9564]. By combining automated and manual peptide synthesis, here we have produced in high yields seven analogues of fragment 1-47 containing natural and non-natural amino acids. In particular, we have replaced Val1 with tert-butylglycine (tBug), Ser2 with Arg, and Tyr3 with Phe, cyclohexylalanine (Cha), Trp, alpha-naphthylalanine (alphaNal), and beta-naphthylalanine (betaNal). The crude reduced peptides are able to fold almost quantitatively into the disulfide-cross-linked species, whose unique alignment (Cys6-Cys14, Cys16-Cys28, and Cys22-Cys37) has been shown to be identical to that of the natural fragment. The results of conformational characterization provide evidence that synthetic peptides retain the structural features of the natural species, whereas thrombin inhibition data indicate that the synthetic analogues are all more potent inhibitors of thrombin. In particular, Val --> tBug exchange leads to a 3-fold increase in binding, interpreted as arising from a favorable reduction of the entropy of binding, due to the presence of the more symmetric side chain of tBug relative to that of Val. The S2R analogue binds 24- and 125-fold more tightly than the natural fragment to the fast or slow form of thrombin. These results are explained by considering that Arg2 may favorably couple to Glu192, a key residue involved in the slow to fast transition, thus stabilizing the slow form. Replacement of Tyr3 with more hydrophobic residues having different side chain orientations and electronic structures improves binding by 2-40-fold, suggesting that nonpolar interactions and shape-dependent packing effects strongly influence binding at this position. Overall, these results provide new insights for elucidating the mechanism of hirudin-thrombin recognition at the molecular level and highlight new strategies for designing more potent and selective inhibitors 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.

Prothrombin Greenville, Arg517→Gln, Identified in an Individual Heterozygous for Dysprothrombinemia

A 64-year-old white male was referred for evaluation of prolonged prothrombin time (PT) and activated partial thromboplastin time (aPTT) obtained before elective surgery with initial PT and PTT results of 14.9 and 38.4 seconds, respectively, which corrected to normal in 1:1 mixes with normal plasma. Functional prothrombin assay indicated a level of 51% with thromboplastin as an activator. The prothrombin antigen was 102%. This discordance in the functional and immunologic prothrombin levels was evidence for dysprothrombinemia. Western blotting showed that thrombin was formed at a normal rate in diluted plasma consistent with a mutation within the thrombin portion of prothrombin. DNA was isolated from leukocytes and the thrombin exons were amplified by polymerase chain reaction, cloned, and sequenced. For exon 13, eight clones were sequenced with four clones showing a point mutation in the codon for Arg517, which would result in substitution by Gln. Arg517 is part of the Arg-Gly-Asp(RGD) sequence in thrombin and contributes to an ion cluster with aspartic acid residues 552 and 554. Mutation at this residue most probably distorts the structure of the Na+ binding site in thrombin. This is the first report indicating the critical role of Arg517 in the normal physiological interaction of thrombin with fibrinogen. This dysprothrombin is designated Prothrombin Greenville.

Prothrombin Greenville, Arg517→Gln, Identified in an Individual Heterozygous for Dysprothrombinemia

A 64-year-old white male was referred for evaluation of prolonged prothrombin time (PT) and activated partial thromboplastin time (aPTT) obtained before elective surgery with initial PT and PTT results of 14.9 and 38.4 seconds, respectively, which corrected to normal in 1:1 mixes with normal plasma. Functional prothrombin assay indicated a level of 51% with thromboplastin as an activator. The prothrombin antigen was 102%. This discordance in the functional and immunologic prothrombin levels was evidence for dysprothrombinemia. Western blotting showed that thrombin was formed at a normal rate in diluted plasma consistent with a mutation within the thrombin portion of prothrombin. DNA was isolated from leukocytes and the thrombin exons were amplified by polymerase chain reaction, cloned, and sequenced. For exon 13, eight clones were sequenced with four clones showing a point mutation in the codon for Arg517, which would result in substitution by Gln. Arg517 is part of the Arg-Gly-Asp(RGD) sequence in thrombin and contributes to an ion cluster with aspartic acid residues 552 and 554. Mutation at this residue most probably distorts the structure of the Na+ binding site in thrombin. This is the first report indicating the critical role of Arg517 in the normal physiological interaction of thrombin with fibrinogen. This dysprothrombin is designated Prothrombin Greenville.

Conserved water molecules in the specificity pocket of serine proteases and the molecular mechanism of Na+ binding

Conservation of clusters of buried water molecules is a structural motif present throughout the serine protease family. Frequently, these clusters are shaped as water channels forming extensive hydrogen-bonding networks linked to the protein backbone. The most conspicuous example is the water channel present in the specificity pocket of trypsin and thrombin. In thrombin, other vitamin K-dependent proteases, and some complement factors, Na+ binds in this water channel and enhances allosterically the catalytic activity of the enzyme, whereas digestive and fibrinolytic proteases are devoid of such regulation. A comparative analysis of proteases with and without Na+ binding capability reveals the role of the water channel in maintaining the structural organization of the specificity pocket and in Na+ coordination. This enables the formulation of a molecular mechanism for Na+ binding in thrombin and leads to the identification of the structural changes necessary to engineer a functional Na+ site and enhanced catalytic activity in trypsin and other proteases.

Kinetic and chemical mechanisms for the effects of univalent cations on the spectral properties of aromatic amine dehydrogenase

Univalent cations and pH influence the UV-visible absorption spectrum of the tryptophan tryptophylquinone (TTQ) enzyme, aromatic amine dehydrogenase (AADH). Little spectral perturbation was observed when pH was varied in the absence of univalent cations. The addition of alkali metal univalent cations (K+, Na+, Li+, Rb+ or Cs+) to oxidized AADH caused significant changes in its absorption spectrum. The apparent Kd for each cation, determined from titrations of the spectral perturbation, decreased with increasing pH. Transient kinetic studies involving rapid mixing of AADH with cations and pH jump revealed that the rate of the cation-induced spectral changes initially decreased with increasing cation concentration to a minimum value, then increased with increasing cation concentration. A kinetic model was developed to fit these data, determine the true pH-independent Kd values for K+ and Na+, and explain the pH-dependence of the apparent Kd. A chemical reaction mechanism, based on the kinetic data, is presented in which the metallic univalent cation facilitates the chemical modification of the TTQ prosthetic group to form an hydroxide adduct which gives rise to the spectral change. Addition of NH4(+)/NH3 to AADH caused changes in the absorption spectrum which were very different form those caused by addition of the metallic univalent cations. The kinetics of the reaction induced by addition of NH4+/NH3 were also different, being simple saturation kinetics. Another reaction mechanism is proposed for the NH4+/NH3-induced spectral change that involves nucleophilic addition of the unprotonated NH3 to TTQ. The general relevance of these data and models to the physiological reactions of TTQ-dependent enzymes and to the roles of univalent cations in modulating enzyme activity are discussed.

Kinetic pathway for the slow to fast transition of thrombin. Evidence of linked ligand binding at structurally distinct domains

The kinetic pathway for the Na+-induced slow --> fast transition of thrombin was characterized. The slow form was shown to consist of two conformers in a 3:1 ratio (ES2:ES1) at 5 degrees C, pH 7.4, Gamma/2 0.3. ES2 binds Na+ 3 orders of magnitude faster than does ES1. The small molecule active site-directed inhibitor L-371,912, and the exosite I binding ligand hirugen, like Na+, bind selectively to ES2 and induce the slow --> fast conversion of thrombin. The slow --> fast transition is limited by the rate of conversion of ES1 to ES2 (k approximately 28 s-1 at 5 degrees C). Replacement of Arg-221a or Lys-224 at the Na+ binding site with Ala appears to selectively alter the slow form and reduce the apparent affinity of the mutants for Na+ and L-371,912. This replacement, however, has little effect on the affinity for the inhibitor in the presence of saturating concentrations of Na+. The kinetically linked ligand binding at the Na+ binding site, exosite I, and the active site of thrombin characterized in the present study indicates the basis for the plasticity of this important enzyme, and suggests the possibility that the substrate specificity and, therefore, the procoagulant and anticoagulant activities of thrombin may be subject to allosteric regulation by as yet unidentified physiologically important effectors.

Catalytic role of monovalent cations in the mechanism of proton transfer which gates an interprotein electron transfer reaction

Within the methylamine dehydrogenase (MADH)-amicyanin protein complex, long range intermolecular electron transfer (ET) occurs between tryptophan tryptophylquinone (TTQ) of MADH and the type I copper of amicyanin. The reoxidations of two chemically distinct reduced forms of TTQ were studied, a quinol (O-quinol) generated by reduction by dithionite and the physiologically relevant aminoquinol (N-quinol) generated by reduction by methylamine. The latter contains a substrate-derived amino group which displaces the C6 carbonyl oxygen on TTQ. ET from N-quinol MADH to amicyanin is gated by the transfer of a solvent exchangeable proton [Bishop, G. R., & Davidson, V. L. (1995) Biochemistry 34, 12082-12086]. The factors which influence this proton transfer (PT) reaction have been examined. The rate of PT increases with increasing pH and with increasing salt concentration. The salt effect is due to specific monovalent cations and is not a general ionic strength effect. The rate enhancements by pH and cations do not reflect an elimination of the PT step that gates ET. Over the range of pH from 5.5 to 9.0 and with cation concentrations from 0 to 200 mM, the observed rate of the redox reaction is still that of PT. This is proven by kinetic solvent isotope effect studies which show that a primary isotope effect persists even at the highest values of pH and cation concentration. A model is presented to explain how specific cations contribute to catalysis and influence the rate of PT in this reaction. The pH dependence is attributed to an ionizable group that is involved in cation binding. The effect of the cation is stabilization of a negatively charged reaction intermediate that is formed during the deprotonation of the N-quinol, and from which rapid ET to the copper of amicyanin occurs. The relevance of these findings to other enzymes which exhibit reaction rates that are influenced by monovalent cations is also discussed.

The co-crystal structure of unliganded bovine alpha-thrombin and prethrombin-2: movement of the Tyr-Pro-Pro-Trp segment and active site residues upon ligand binding

Unliganded bovine alpha-thrombin and prethrombin-2 have been co-crystallized, in space group P21212, using either ammonium sulfate or polyethylene glycol 2000 (PEG2K), and their structures determined at 2.2 A and 2.3 A, respectively. Initial phases were determined by molecular replacement and refined using XPLOR to final R factors of 0.187 (Rfree = 0.255) and 0.190 (Rfree = 0.282) for the salt and PEG2K models, respectively. The apo-enzyme form of bovine alpha-thrombin shows dramatic shifts in placement for the Tyr-Pro-Pro-Trp segment, for Glu-192, and for the catalytic residues His-57 and Ser-195, when compared to 4 thrombin complexes representing different states of catalysis, namely (1) the Michaelis complex (residues 7-19 of fibrinogen A alpha with a non-cleavable scissile bond), (2) enzyme-inhibitor complex (D-Phe-Pro-Arg chloromethylketone), (3) enzyme product complex (residues 7-16 of fibrinopeptide A), and (4) the exosite complex (residues 53-64 of hirudin). The structures of bovine and human prethrombin-2 are generally similar to one another (RMS deviation of 0.68 A) but differ significantly in the Arg-15/Ile-16 cleavage region and in the three activation domains, which are disordered in bovine prethrombin-2, analogous to that seen for trypsinogen.

Human alpha-thrombin inhibition by the highly selective compounds N-ethoxycarbonyl-D-Phe-Pro-alpha-azaLys p-nitrophenyl ester and N-carbobenzoxy-Pro-alpha-azaLys p-nitrophenyl ester: a kinetic, thermodynamic and X-ray crystallographic study

Kinetics, thermodynamics and structural aspects of human alpha-thrombin (thrombin) inhibition by newly synthesized low molecular weight derivatives of alpha-azalysine have been investigated. The thrombin catalyzed hydrolysis of N-ethoxycarbonyl-D-Phe-Pro-alpha-azaLys p-nitrophenyl ester (Eoc-D-Phe-Pro-azaLys-ONp) and N-carbobenzoxy-Pro-alpha-azaLys p-nitrophenyl ester (Cbz-Pro-azaLys-ONp) was investigated at pH 6.2 and 21.0 degrees C, and analyzed in parallel with that of N-alpha-(N,N-dimethylcarbamoyl)-alpha-azalysine p-nitrophenyl ester (Dmc-azaLys-ONp). Decarboxylation following the enzymatic hydrolysis of these p-nitrophenyl esters gave the corresponding 1-peptidyl-2(4-aminobutyl) hydrazines (peptidyl-Abh) showing properties of thrombin competitive inhibitors. Therefore, thermodynamics for the reversible binding of D-Phe-Pro-Abh, Cbz-Pro-Abh and Dmc-Abh to thrombin was examined. These results are consistent with the minimum four-step catalytic mechanism for product inhibition of serine proteinases. Eoc-D-Phe-Pro-azaLys-ONp and Eoc-D-Phe-Pro-Abh display a sub-micromolar affinity for thrombin together with a high selectivity versus homologous plasmatic and pancreatic serine proteinases acting on cationic substrates. The three-dimensional structures of the reversible non-covalent thrombin:Eoc-D-Phe-Pro-Abh and thrombin:Cbz-Pro-Abh complexes have been determined by X-ray crystallography at 2.0 A resolution (R-factor = 0.169 and 0.179, respectively), and analyzed in parallel with that of the thrombin:Dmc-azaLys acyl-enzyme adduct. Both Eoc-D-Phe-Pro-Abh and Cbz-Pro-Abh competitive inhibitors are accommodated in the thrombin active center, spanning the region between the aryl binding site and the S1 primary specificity subsite.

Energetics of thrombin-thrombomodulin interaction

Temperature and salt dependence studies of thrombin interaction with thrombomodulin, with and without chondroitin sulfate, and two fragments containing the EGF-like domains 4-5 and 4-5-6 reveal the energetic signatures and the mechanism of recognition of this physiologically important cofactor. Binding of thrombomodulin is affected drastically by the particular salt present in solution and is positively linked to Na+ binding to thrombin and the conversion of the enzyme from the slow to the fast form, but is opposed by Cl- binding to the fibrinogen recognition site and especially to the heparin binding site. Binding of thrombomodulin has an unusually large salt dependence (gamma(salt) = -4.8) contributed mostly by the polyelectrolyte-like nature of the chondroitin sulfate moiety that binds to the heparin binding site and increases the affinity of the cofactor by almost 10-fold. On the other hand, the chondroitin sulfate has no effect on the deltaCp of binding, which is determined predominantly by contacts made by the EGF-like domains 5 and 6 with the fibrinogen recognition site. The modest heat capacity change (-0.2 kcal mol(-1) K(-1)) observed when thrombomodulin binds to the fast form suggests a rigid-body association of the cofactor with the enzyme. In the slow form, however, the heat capacity change is significantly more pronounced (-0.5 kcal mol(-1) K(-1)) and signals the presence of a conformational transition of the enzyme linked to binding of the cofactor that mimics the slow-->fast conversion. These results demonstrate that recognition of thrombomodulin by thrombin is steered electrostatically by the highly charged regions of the fibrinogen recognition site and the heparin binding site, to which the chondroitin sulfate moiety binds and enhances the affinity of the interaction. The recognition event also involves conformational changes of the enzyme in the slow form mediated by binding of the EGF-like domains 5-6 to the fibrinogen recognition site. Consistent with this model, binding of thrombomodulin to the fast form has only a small effect on the hydrolysis of nine chromogenic substrates carrying substitutions at P1, P2, and P3 aimed at probing the environment of the specificity sites S1, S2, and S3 of the enzyme. Binding to the slow form, on the other hand, enhances the specificity toward all substrates up to 15-fold. For substrates carrying a Gly at P2, binding of thrombomodulin changes the relative specificity of the slow and fast forms and makes the slow form more specific. Interestingly, these effects are not specific of thrombomodulin and depend solely on binding to the fibrinogen recognition site of the enzyme. In fact, they are also observed with the hirudin C-terminal fragment 55-65. The characterization of the mechanism of thrombin-thrombomodulin interaction and the effects of the cofactor on the hydrolysis of chromogenic substrates probing the interior of the catalytic pocket bear on the thrombomodulin-induced enhancement of protein C cleavage by thrombin. We propose that this enhancement is due predominantly to an effect of thrombomodulin on the bound protein C in the ternary complex. Therefore, thrombomodulin would carry out its physiological function by making protein C a better substrate for thrombin, rather than making thrombin a better enzyme for protein C.

Functional requirements for inhibition of thrombin by antithrombin III in the presence and absence of heparin

Mutation of 79 highly exposed amino acids that comprise approximately 62% of the solvent accessible surface of thrombin identified residues that modulate the inhibition of thrombin by antithrombin III, the principal physiological inhibitor of thrombin. Mutations that decreased the susceptibility of thrombin to inhibition by antithrombin III in the presence and absence of heparin (W50A, E229A, and R233A) also decreased hydrolysis of a small tripeptidyl substrate. These residues were clustered around the active site cleft of thrombin and were predicted to interact directly with the "substrate loop" of antithrombin III. Despite the large size of antithrombin III (58 kDa), no residues outside of the active cleft were identified that interact directly with antithrombin III. Mutations that decreased the susceptibility of thrombin to inhibition by antithrombin III in the presence but not in the absence of heparin (R89A/R93A/E94A, R98A, R245A, K248A, K252A/D255A/Q256A) in general did not also affect hydrolysis of the tripeptidyl substrate. These residues were clustered among a patch of basic residues on a surface of thrombin perpendicular to the face containing the active site cleft and were predicted to interact directly with heparin. Three mutations (E25A, R178A/R180A/D183A, and E202A) caused a slight enhancement of inhibition by antithrombin III.

Critical role of W60d in thrombin allostery

Mutation of residue W60d of thrombin, located 17 A from the Na+ binding site, suppresses Na+ binding and the functional differences between the slow and fast forms. The molecular basis for the long-range effect of this mutation is provided by a conspicuous network of water molecules which connects the Na+ binding environment to the specificity sites S1 and S2 of the enzyme. The mutation appears to stabilize thrombin in a hybrid conformation that is overall similar to the slow form, but with the fibrinogen recognition site functioning as in the fast form. It also affects the switch in specificity from fibrinogen to protein C linked to the release Na+ and the fast-->slow conversion. Under physiological conditions of pH, temperature and NaCl concentration, the W60dS mutant behaves as an anticoagulant. It has a reduced activity toward fibrinogen by 22-fold, while the reduction of protein C activation in the presence of saturating concentrations of thrombomodulin is less than 2-fold. Even more remarkable is the cleavage of fibrin I monomer leading to release of fibrinopeptide B, which is reduced by more than 130-fold. This property is reminiscent of the snake venom ancrod, which only releases fibrinopeptide A, and adds substantially to the anticoagulant potency of the W60dS mutant. In fact, the clotting time in the presence of this mutant is prolonged more than 40-fold compared to the wild-type.

Contribution of lysine 60f to S1' specificity of thrombin

Lys60f has been proposed to limit the S1' substrate binding site specificity of thrombin to small polar P1' residues by occluding the S1' binding pocket, based on the X-ray crystal structure of thrombin. To test this proposal, we prepared a Lys-->Ala (K60fA) mutant of recombinant thrombin and determined whether this mutation enhanced the reactivity of thrombin with a variant inhibitor [antithrombin (AT)-Denver] and a substrate (protein C) containing poorly recognized P1' Leu residues. AT-Denver in the presence of heparin inhibited K60fA thrombin with a second-order association rate constant [k = 4.2 +/- 0.1) x 10(5) M-1 s-1] that was 3.2-fold faster than thrombin [k = (1.3 +/- 0.1) x 10(5) M-1 s-1]. Wildtype AT (P1' Ser) under the same conditions inhibited K60fA thrombin with a 2.5-fold slower rate constant [k = (1.1 +/- 0.1) x 10(7) M-1 s-1] than thrombin [k = (2.8 +/- 0.1) x 10(7) M-1 s-1]. These results indicate an overall 8.3-fold improvement in the recognition of the P1' Leu of AT-Denver by K60fA thrombin over that of wild-type thrombin; i.e., the K60fA mutation partly overcomes the defect in thrombin inhibition produced by the P1' mutation in AT-Denver. Resolution of the two-step reactions of AT and AT-Denver with wild-type and mutant thrombins revealed that the enhanced recognition of P1' Leu in AT-Denver by K60fA thrombin occurs primarily in the second reaction step in which a noncovalent AT-thrombin encounter complex is converted to a stable, covalent complex. Thrombin K60fA activated Gla-domainless protein C (GDPC) approximately 2- and approximately 4-fold faster than thrombin in the presence and absence of thrombomodulin (TM), respectively, consistent with an improved interaction of the Leu P1' residue with the mutant S1' pocket. In contrast, the mutant thrombin clotted fibrinogen (P1' Gly) approximately 3-fold slower than thrombin. Kinetic analysis revealed that the improvement in the catalytic rate of activation of GDPC by K60fA thrombin in the presence of TM was localized in the second reaction step, as reflected by an approximately 2-fold increase in kcat. Direct binding studies showed that the K60fA mutation minimally affected the affinity of thrombin for Na+, indicating that the changes in S1' site-specificity of K60fA thrombin did not result from altering the allosteric transition induced by Na+. We conclude that Lys60f limits the P1' substrate and inhibitor specificity of thrombin by influencing the size and polarity of the S1' site which thereby affects the stability of the transition state for cleavage of the scissile bond in the second reaction step.

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.

The molecular environment of the Na+ binding site of thrombin

When Na+ binds to thrombin, a conformational change is induced that renders the enzyme kinetically faster and more specific in the activation of fibrinogen. Two Na+ binding sites have here been identified crystallographically by exchanging Na+ with Rb+. One is intermolecular, found on the surface between two symmetry-related thrombin molecules. Since it is not present in thrombin crystal structures having different crystal systems, the other Na+ site is the functionally relevant one. The second site has octahedral coordination with the carbonyl oxygen atoms of Arg221A and Lys224 and four conserved water molecules. It is located near Asp189 of the S1 specificity site in an elongated solvent channel (8 x 18 A) formed by four antiparallel beta-strands between Cys182-Cys191 and Val213-Tyr228. This channel, extending from the active site to the opposite surface of the enzyme, was first noted in the hirudin-thrombin structure and contains about 20 conserved water molecules linked together by a hydrogen bonding network that connects to the main chain of thrombin. Although the antiparallel beta-strand interactions of the functional Na+ binding site are the same in prethrombin2, the loops between the strands are very different, so that Asp189 and Arg221A are not positioned properly for either substrate or Na+ binding in prethrombin2. A water molecule with octahedral coordination has also been identified in factor Xa at the topologically equivalent Na+ site position of thrombin. Since activated protein C shows enhanced activity with monovalant cation binding, the same position is probably utilized by Na+. Since thrombin crystals could not be grown in the absence of Na+, the cation was leached from Na(+)-bound thrombin crystals by diffusion/exchange. Although both Na+ and their coordinating water molecules were removed from the Na+ binding sites, the remainder of the thrombin structure was, unexpectedly, the same. The lack of an allosteric change is most likely attributable to crystal packing effects. Thus, the structure of the slow form remains to be established crystallographically.

Identification of surface residues mediating tissue factor binding and catalytic function of the serine protease factor VIIa

Factor VIIa (VIIa), the serine protease that initiates the coagulation pathways, is catalytically activated upon binding to its cell surface receptor and cofactor tissue factor (TF). This study provides a comprehensive analysis of the functional surface of VIIa by alanine scanning mutagenesis of 112 residues. Residue side chains were defined which contribute to TF binding and factor X hydrolysis. Energetically important binding contacts at the interface with TF were identified in the first epidermal growth factor domain of VIIa (Gln-64, Ile-69, Phe-71, Arg-79) and in the protease domain (Arg-277, Met-306, Asp-309). The observed energetic defects are in good agreement with the corresponding residues in TF, suggesting that the VIIa light chain plays a prominent role in high affinity binding of cofactor. Mutation of protease domain interface residues indicated that TF allosterically influences the active site of VIIa. Stabilization of a labile zymogen to enzyme transition could explain the activating effect of TF on VIIa catalytic function. Residues important for factor X hydrolysis were found in three regions of the protease domain: (i) specificity determinants in the catalytic cleft and adjacent loops, (ii) an exosite near the TF binding site, and (iii) a large electronegative exosite which is in a position analogous to the basic exosite I of thrombin. TF regions involved in factor X activation are positioned on the same face of the TF-VIIa complex as the two exosites identified on the protease domain surface, providing evidence for an extended interaction of TF-VIIa with macromolecular substrate.

Enzyme flexibility, solvent and 'weak' interactions characterize thrombin-ligand interactions: implications for drug design

BACKGROUND

The explosive growth in the rate of X-ray determination of protein structures is fuelled largely by the expectation that structural information will be useful for pharmacological and biotechnological applications. For example, there have been intensive efforts to develop orally administrable antithrombotic drugs using information about the crystal structures of blood coagulation factors, including thrombin. Most of the low molecular weight thrombin inhibitors studied so far are based on arginine and benzamidine. We sought to expand the database of information on thrombin-inhibitor binding by studying new classes of inhibitors.

RESULTS

We report the structures of three new inhibitors complexed with thrombin, two based on 4-aminopyridine and one based on naphthamidine. We observe several geometry changes in the protein main chain and side chains which accompany inhibitor binding. The two inhibitors based on 4-aminopyridine bind in notably different ways: one forms a water-mediated hydrogen bond to the active site Ser195, the other induces a rotation of the Ser214-Trp215 peptide plane that is unprecedented in thrombin structures. These binding modes also differ in their 'weak' interactions, including CH-O hydrogen bonds and interactions between water molecules and aromatic pi-clouds. Induced-fit structural changes were also seen in the structure of the naphthamidine inhibitor complex.

CONCLUSIONS

Protein flexibility and variable water structures are essential elements in protein-ligand interactions. Ligand design strategies that fail to take this into account may overlook or underestimate the potential of lead structures. Further, the significance of 'weak' interactions must be considered both in crystallographic refinement and in analysis of binding mechanisms.

Crystal structures of thrombin with thiazole-containing inhibitors: probes of the S1' binding site

Structures of the blood clotting enzyme thrombin complexed with hirugen and two active site inhibitors, RWJ-50353 10080(N-methyl-D-phenylalanyl-N-[5-[(aminoiminomethyl)amino]-1- [[(2-benzothiazolyl)carbonyl]butyl]-L-prolinamide trifluoroacetate hydrate) and RWJ-50215 (N-[4-(aminoiminomethyl)amino-1-[2- (thiazol-2-ylcarbonylethyl)piperidin- 1-ylcarbonyl]butyl]-5-(dimethylamino)naphthalenesulfonamide trifluoroacetate hydrate), were determined by x-ray crystallography. The refinements converged at R values of 0.158 in the 7.0-2.3-A range for RWJ-50353 and 0.155 in the 7.0-1.8-A range for RWJ-50215. Interactions between the protein and the thiazole rings of the two inhibitors provide new valuable information about the S1' binding site of thrombin. The RWJ-50353 inhibitor consists of an S1'-binding benzothiazole group linked to the D-Phe-Pro-Arg chloromethyl ketone motif. Interactions with the S1-S3 sites are similar to the D-phenylalanyl-prolyl-arginyl chloromethylketone structure. In RWJ-50215, a S1'-binding 2-ketothiazole group was added to the thrombin inhibitor-like framework of dansylarginine N-(3-ethyl-1,5-pentanediyl)amide. The geometry at the S1-S3 sites here is also similar to that of the parent compound. The benzothiazole and 2-ketothiazole groups bind in a cavity surrounded by His57, Tyr60A, Trp60D, and Lys60F. This location of the S1' binding site is consistent with previous structures of thrombin complexes with hirulog-3, CVS-995, and hirutonin-2 and -6. The ring nitrogen of the RWJ-50353 benzothiazole forms a hydrogen bond with His57, and Lys60F reorients because of close contacts. The oxygen and nitrogen of the ketothiazole of RWJ-50215 hydrogen bond with the NZ atom of Lys60F.

Cation selective promotion of tubulin polymerization by alkali metal chlorides

A role for charge-based interactions in protein stability at the monomer or dimer level is well known. We show here that such interactions can also be important for the higher-order structures of microtubule assembly. Alkali metal chlorides increase the rate of polymerization of pure tubulin driven by either taxol or dimethyl sulfoxide. The effect is cation selective, exhibiting a sequence Na+ > K+ > Li+ > Cs+, with optimal concentrations for Na+ at approximately 160 mM. Hofmeister anion effects are additive with these rate stimulations. Sodium is less potent than guanidinium ion stimulation reported previously, but produces a larger fraction of normal microtubules. Alkali metal cations lower the critical concentration by a factor of approximately 2, produce cold reversible polymers whose formation is sensitive to podophyllotoxin inhibition, increase the fraction of polymers present as microtubules from approximately 0.9 to 0.99, and reverse or prevent urea-induced depolymerization of microtubules. In the presence of microtubule-associated proteins, the promotion of polymerization is no longer cation selective. In the polymerization of tubulin S, in which the acidic C termini of both monomers have been cleaved, the cation enhancement is markedly decreased, although selective persists. Because the selectivity sequence is similar to that of the coil/helix transition of polyglutamic acid, we suggest that a major part, although not all, of the cation selective enhancement of polymerization results from shielding of the glutamate-rich C termini of the tubulin monomers.

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.

Electrostatic interactions in hirudin-thrombin binding

Hirudin is a good anticoagulant owing to potent inhibition of the serine protease thrombin. An aspartate- and glutamate-rich portion of hirudin plays an important part in its tight binding to thrombin through a ladder of salt bridges, and these residues have previously been mutated to asparagine or glutamine. Detailed calculations of the electrostatic contribution to changes in binding from these mutations have been performed using the finite-difference Poisson-Boltzmann method which include charge--charge interactions, solvation interactions, the residual electrostatic interaction of mutant residues, pKa shifts, and ionic strength. Single mutant effects on binding energy were close to experimental values, except for the D55N mutant whose effect is overestimated, perhaps because of displacement of a bound chloride ion from the site where it binds. Multiple mutation values were generally overestimated. The effect of pKa shifts upon the binding is significant for one hirudin residue E58, but this appears to be due to a poor salt bridge with thrombin caused by crystal contacts. Electrostatic interaction between the acidic residues is unfavorable. However, analysis of experimental multiple mutation/single mutation data shows apparently negative interactions between these residues, from which it is concluded that structural changes can occur in the complex to relieve an unfavorable interaction when more than one acidic residue is mutated. In all cases, there is a loss in stability of the complex from mutations due to loss of favorable charge--charge interactions with thrombin, but this is largely compensated for by reduced unfavorable desolvation interactions, and by residual polar interactions in the Asn/Gln mutants.

Release of fibrinopeptides by the slow and fast forms of thrombin

The release of fibrinopeptides A and B by the slow and fast forms of thrombin was studied over the temperature range from 5 to 45 degrees C and the salt concentration range from 100 to 800 mM. The sequential mechanism for the release of fibrinopeptides originally proposed by Shafer was found to be obeyed under all conditions examined. The origin of preferential binding of fibrinogen and fibrin I to the fast form of thrombin in the transition state is in the second-order rate constant for association, k(l). In the case of fibrinogen, the values of k(l) for interaction with the fast and slow forms at 25 degrees C are 19 +/- 4 and 2.5 +/- 0.3 microM(-1) s(-1), with an activation energy of about 10 kcal/mol in both forms. In the case of fibrin I, the analogous values of k(l) are 9.1 +/- 0.7 and 2.5 +/- 0.2 microM(-1) s(-1), and the activation energy is about 4.5 kcal/mol in both forms. The mechanism of recognition of fibrinogen and fibrin I by thrombin entails a diffusion-controlled step with a small energy barrier. Analysis of the temperature dependence of the coupling free energy for allosteric switching indicates that the preferential interaction of fibrinogen and fibrin I with the fast form of thrombin in the transition state is entropy-driven, signaling a contribution of the hydrophobic effect to the slow-->fast transition. The salt dependence of the release of fibrinopeptides shows a constant coefficient Gamma(salt) = d ln(k(cat)/K(m))/d ln [salt] in the concentration range examined. Interestingly, the value of Gamma(salt) is independent of the salt used (NaCl, ChCl, or NaF) and is -1.5 +/- 0.1 for fibrinopeptide A and -2.5 +/- 0.1 for fibrinopeptide B. Hence, Gamma(salt) reflects predominantly the electrostatic contribution to the formation of the transition state, with a larger contribution seen in the interaction of thrombin with fibrin I. It is concluded that the interaction of thrombin with fibrinogen and fibrin I, leading to the release of fibrinopeptides A and B, is driven by electrostatic forces that presumably favor the correct preorientation of the enzyme and the substrate to form a productive complex in the transition state. This electrostatic-steering effect, also reported for thrombin-hirudin interaction, leads to a diffusion-controlled encounter with a very small energy barrier. Once the complex is formed, the enzyme switches to the fast form as a result of entropic factors presumably linked to water release from a more extended surface of recognition. While the release of fibrinopeptides as a function of salt concentration was being studied, an important observation was made on the role of Cl- in the formation of the fibrin clot. This anion drastically and specifically reduces the thickness of fibrin fibers, as judged by the 10-fold decrease in the equilibrium turbidity of clots developed in NaCl as compared to the turbidity of clots developed in NaF. Hence, the transition from a "coarse" to a "fine" clot induced by an increase in ionic strength as first described by Ferry is, instead, due to the specific binding of Cl- to intermediates in the ensuing polymerization. In fact, no change in the clotting curve is observed when the ionic strength is changed with NaF.

Heparin enhances the catalytic activity of des-ETW-thrombin

The thrombin mutant, des-ETW-thrombin, lacking Glu(146), Thr(147), and Trp(148) within a unique insertion loop located at the extreme end of the primary specificity pocket, has been shown previously to exhibit reduced catalytic activity with respect to macromolecular and synthetic thrombin substrates and reduced or enhanced susceptibility to inhibition. Investigation of the hydrolysis of peptidyl p-nitroanilide substrates by des-ETW-thrombin showed increased activity in the presence of heparin and other sulphated glycosaminoglycans. No effect was observed upon the activity of wild-type thrombin. Heparin was found to decrease the K(m) for cleavage of four thrombin-specific substrates by des-ETW-thrombin by 3-4-fold. Similarly, pentosan polysulphate (PPS) decreased the K(m) with these substrates by 8-10-fold. Heparin also increased the rate of inhibition of des-ETW-thrombin by antithrombin III and D-phenylalanyl-prolyl-arginylchloromethane (PPACK). The inhibition of des-ETW-thrombin by a number of thrombin-specific peptide boronic acids also showed significant reduction in the final K(i) in the presence of heparin, due to reduction in the off-rate. A peptide analogue of a sequence of hirudin which binds thrombin tightly to exosite I (fibrinogen recognition site) potentiated the activity of des-ETW-thrombin against peptide p-nitroanilide substrates in a manner similar to heparin. The K(i) for the inhibition of des-ETW-thrombin by p-aminobenzamidine was decreased by these ligands from 9.7 mM to 7.5 mM, 5.1 mM, and 2.5 mM in the presence of heparin, hirudin peptide and PPS respectively, suggesting the increased catalytic activity is due to enhanced access to the primary specificity pocket. The positive influence of these ligands on des-ETW-thrombin was reversed in the presence of ATP or ADP; the latter has previously been shown to inhibit thrombin activity by blocking initial interaction with fibrinogen at exosite 1. Because the effect of heparin and PPS is similar to that of hirudin peptide, it is proposed that the most likely mechanism is that binding at the heparin-binding site (thrombin exosite 2) facilitates interaction at exosite 1 causing a conformational change which partially corrects the defective ground-state binding of the mutant thrombin. Although no effect was observed upon the activity of wild-type thrombin, our findings do provide further evidence of an allosteric property of thrombin which may control the geometry of, and access to, the primary specificity pocket.

Functionality map analysis of the active site cleft of human thrombin

The Multiple Copy Simultaneous Search methodology has been used to construct functionality maps for an extended region of human thrombin, including the active site. This method allows the determination of energetically favorable positions and orientations for functional groups defined by the user on the three-dimensional surface of a protein. The positions of 10 functional group sites are compared with those of corresponding groups of four thrombin-inhibitor complexes. Many, but not all features, of known thrombin inhibitors are reproduced by the method. The results indicate that certain aspects of the binding modes of these inhibitors are not optimal. In addition, suggestions are made for improving binding by interaction with functional group sites on the thrombin surface that are not used by the thrombin inhibitors.

Theory of allosteric effects in serine proteases

The classical Botts-Morales theory for the action of a modifier on the catalytic properties of an enzyme has been extended to deal with allosteric effects in serine proteases. The exact analytical solution derived for the linkage scheme at steady state provides a rigorous framework for the study of many biologically relevant systems, including enzymes activated by monovalent cations and cofactor-controlled protease-zymogen interactions in blood coagulation. When the enzyme obeys Michaelis-Menten kinetics, the exact solution of the kinetic linkage scheme simplifies considerably. Of particular importance for practical applications is a simple equation expressing the dependence of the specificity constant of the enzyme, kcat/Km, on the concentration of the modifier, from which the equilibrium binding constant for the formation of the enzyme-modifier complex can be estimated. Analysis of the allosteric changes in thrombin activity induced by thrombomodulin and Na+ in terms of this equation yields accurate determinations of the equilibrium binding constants for both effectors.

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

Residues energetically linked to the allosteric transition of thrombin from its anticoagulant slow form to the procoagulant fast form have been identified by site-directed mutagenesis. The energetics of recognition by the two forms of the enzyme were probed by using a synthetic chromogenic substrate, fibrinogen, and hirudin. The thrombin residues E39, W60d, E192, D221, and D222 are linked to the slow-->fast transition and are part of an "allosteric core" through which events originating at the Na+ binding loop propagate to other regions of the enzyme. The thrombin residues Y76, W96, W148, and R173 lie at the periphery of the allosteric core, affect recognition of fibrinogen and hirudin to the same extent in both forms, and are not linked to the slow-->fast transition.