The Na+ binding site of thrombin

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

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

Hirudin Binding Reveals Key Determinants of Thrombin Allostery

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

Allosteric modulation of monomeric proteins*

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

Exosite-interactive regions in the A1 and A2 domains of factor VIII facilitate thrombin-catalyzed cleavage of heavy chain

Thrombin catalyzes the proteolytic activation of factor VIII, cleaving two sites in the heavy chain and one site in the light chain of the procofactor. Evaluation of thrombin binding the reaction products from heavy chain cleavage by steady state fluorescence energy transfer using a fluorophore-labeled, active site-modified thrombin as well as by solid phase binding assays using a thrombin Ser(205) --> Ala mutant indicated a high affinity site in the A1 subunit (K(d) approximately 5 nm) that was dependent upon the Na(+)-bound form of thrombin, whereas a moderate affinity site in the A2 subunit (K(d) approximately 100 nm) was observed for both Na(+)-bound and -free forms. The solid phase assay also indicated that hirudin blocked thrombin interaction with the A1 subunit and had little, if any, effect on its interaction with the A2 subunit. Conversely, heparin blocked thrombin interaction with the A2 subunit and showed a marginal effect on A1 binding. Evaluation of the A2 sequence revealed two regions rich in acidic residues that are localized close to the N and C termini of this domain. Peptides encompassing these clustered acidic regions, residues 373-395 and 719-740, blocked thrombin cleavage of the isolated heavy chain at Arg(372) and Arg(740) and inhibited A2 binding to thrombin Ser(205) --> Ala, suggesting that both A2 domain regions potentially support interaction with thrombin. A B-domainless, factor VIII double mutant Asp(392) --> Ala/Asp(394) --> Ala was constructed, expressed, and purified and possessed specific activity equivalent to a severe hemophilia phenotype. This mutant was resistant to cleavage at Arg(740), whereas cleavage at Arg(372) was not affected. These data suggest the acidic region comprising residues 389-394 in factor VIII A2 domain interacts with thrombin via its heparin-binding exosite and facilitates cleavage at Arg(740) during procofactor activation.

Allosteric regulation of the cofactor-dependent serine protease coagulation factor VIIa

The integration of structure and function analysis of the tissue factor-factor VIIa complex has provided a detailed view of the functional surface of the extrinsic activation complex. An incomplete zymogen to enzyme transition is responsible for the strict cofactor dependence of catalytic function of factor VIIa. The mutational analysis demonstrates that factor VIIa is allosterically regulated by specific conformational linkages that involve the cofactor binding site, the catalytic cleft, and the macromolecular substrate exosite. Regions of the flexible activation domain appear to play an important role in the allosteric regulation of this cofactor-dependent coagulation serine protease.

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

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

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.

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

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

PDBSiteScan: a program for searching for active, binding and posttranslational modification sites in the 3D structures of proteins

PDBSiteScan is a web-accessible program designed for searching three-dimensional (3D) protein fragments similar in structure to known active, binding and posttranslational modification sites. A collection of known sites we designated as PDBSite was set up by automated processing of the PDB database using the data on site localization in the SITE field. Additionally, protein-protein interaction sites were generated by analysis of atom coordinates in heterocomplexes. The total number of collected sites was more than 8100; they were assigned to more than 80 functional groups. PDBSiteScan provides automated search of the 3D protein fragments whose maximum distance mismatch (MDM) between N, Calpha and C atoms in a fragment and a functional site is not larger than the MDM threshold defined by the user. PDBSiteScan requires perfect matching of amino acids. PDBSiteScan enables recognition of functional sites in tertiary structures of proteins and allows proteins with functional information to be annotated. The program PDBSiteScan is available at http://wwwmgs.bionet.nsc.ru/mgs/systems/fastprot/pdbsitescan.html.

Proton inventory studies of alpha-thrombin-catalyzed reactions of substrates with selected P and P' sites

Deuterium kinetic solvent isotope effects for the human alpha-thrombin-catalyzed hydrolysis of (1) substrates with selected P(1)-P(3) sites, Z-Pro-Arg-7-amido-4-methylcoumarin (7-AMC), N-t-Boc-Val-Pro-Arg-7-AMC, Bz-Phe-Val-Arg-4-nitroanilide (pNA), and H-D-Phe-L-Pip-Arg-pNA, are (DOD)k(cat) = (2.8-3.3) +/- 0.1 and (DOD)(k(cat)/K(m)) = (0.8-2.1) +/- 0.1 and (2) internally fluorescence-quenched substrates (a) (AB)Val-Phe-Pro-Arg-Ser-Phe-Arg-Leu-Lys(DNP)-Asp-OH, an optimal sequence, and (b) (AB)Val-Ser-Pro-Arg-Ser-Phe-Gln-Lys(DNP)-Asp-OH, recognition sequence for factor VIII, are (DOD)k(cat) = 2.2 +/- 0.2 and (DOD)(k(cat)/K(m)) = (0.8-0.9) +/- 0.1, at the pL (L = H, D) maximum, 8.4-9.0, and (25.0-26.0) +/- 0.1 degrees C. The most plausible models fitting the partial isotope effect (proton inventory) data have been selected on the basis of lowest values of the reduced chi squared and consistency of fractionation factors at all substrate concentrations, assuming rate-determining acylation. The data for Z-Pro-Arg-7-AMC are consistent with a single-proton bridge at the transition state phi(TS) = 0.39 +/- 0.05 and components for solvent reorganization phi(S) = 0.8 +/- 0.1 and phi(S) = 1.22 for k(cat) and k(cat)/K(m), respectively. The data for tripeptide amides fit bowl-shaped curves; an example is N-t-Boc-Val-Pro-Arg-7-AMC: phi(TS)(1) = phi(TS)(2) = 0.57 +/- 0.01 and phi(S) = 1 for k(cat) and 1.6 +/- 0.1 for k(cat)/K(m). Proton inventories for the nonapeptide (2b) are linear. The data for k(cat) for H-D-Phe-L-Pip-Arg-pNA and the decapeptide (2a) are most consistent with two identical fractionation factors for catalytic proton bridging, phi(TS)(1) = phi(TS)(2) = 0.68 +/- 0.02 and a large inverse component (phi(S) = 3.1 +/- 0.5) for the latter, indicative of substantial solvent reorganization upon leaving group departure. Proton inventory curves for k(cat)/K(m) for nearly all substrates are dome-shaped with an inverse isotope effect component (phi(S) = 1.2-2.4) originating from solvent reorganization during association of thrombin with substrate. These large contributions from medium effects are in full accord with the conformational adjustments required for the fulfillment of the dual, hemostatic and thrombolytic, functions of thrombin.

Structure of the antithrombin-thrombin-heparin ternary complex reveals the antithrombotic mechanism of heparin

The maintenance of normal blood flow depends completely on the inhibition of thrombin by antithrombin, a member of the serpin family. Antithrombin circulates at a high concentration, but only becomes capable of efficient thrombin inhibition on interaction with heparin or related glycosaminoglycans. The anticoagulant properties of therapeutic heparin are mediated by its interaction with antithrombin, although the structural basis for this interaction is unclear. Here we present the crystal structure at a resolution of 2.5 A of the ternary complex between antithrombin, thrombin and a heparin mimetic (SR123781). The structure reveals a template mechanism with antithrombin and thrombin bound to the same heparin chain. A notably close contact interface, comprised of extensive active site and exosite interactions, explains, in molecular detail, the basis of the antithrombotic properties of therapeutic heparin.

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

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

Crystal Structure of the C-terminal Peptidoglycan-binding Domain of Human Peptidoglycan Recognition Protein Iα

Peptidoglycan recognition proteins (PGRPs) are pattern recognition receptors of the innate immune system that bind, and in some cases hydrolyze, peptidoglycans (PGNs) on bacterial cell walls. These molecules, which are highly conserved from insects to mammals, participate in host defense against both Gram-positive and Gram-negative bacteria. We report the crystal structure of the C-terminal PGN-binding domain of human PGRP-Ialpha in two oligomeric states, monomer and dimer, to resolutions of 2.80 and 1.65 A, respectively. In contrast to PGRPs with PGN-lytic amidase activity, no zinc ion is present in the PGN-binding site of human PGRP-Ialpha. The structure reveals that PGRPs exhibit extensive topological variability in a large hydrophobic groove, located opposite the PGN-binding site, which may recognize host effector proteins or microbial ligands other than PGN. We also show that full-length PGRP-Ialpha comprises two tandem PGN-binding domains. These domains differ at most potential PGN-contacting positions, implying different fine specificities. Dimerization of PGRP-Ialpha, which occurs through three-dimensional domain swapping, is mediated by specific binding of sodium ions to a flexible hinge loop, stabilizing the conformation found in the dimer. We further demonstrate sodium-dependent dimerization of PGRP-Ialpha in solution, suggesting a possible mechanism for modulating PGRP activity through the formation of multivalent adducts.

Structural and functional basis of the serine protease-like hepatocyte growth factor beta-chain in Met binding and signaling

Hepatocyte growth factor (HGF), a plasminogen-related growth factor, is the ligand for Met, a receptor tyrosine kinase implicated in development, tissue regeneration, and invasive tumor growth. HGF acquires signaling activity only upon proteolytic cleavage of single-chain HGF into its alpha/beta heterodimer, similar to zymogen activation of structurally related serine proteases. Although both chains are required for activation, only the alpha-chain binds Met with high affinity. Recently, we reported that the protease-like HGF beta-chain binds to Met with low affinity (Stamos, J., Lazarus, R. A., Yao, X., Kirchhofer, D., and Wiesmann, C. (2004) EMBO J. 23, 2325-2335). Here we demonstrate that the zymogen-like form of HGF beta also binds Met, albeit with 14-fold lower affinity than the protease-like form, suggesting optimal interactions result from conformational changes upon cleavage of the single-chain form. Extensive mutagenesis of the HGF beta region corresponding to the active site and activation domain of serine proteases showed that 17 of the 38 purified two-chain HGF mutants resulted in impaired cell migration or Met phosphorylation but no loss in Met binding. However, reduced biological activities were well correlated with reduced Met binding of corresponding mutants of HGF beta itself in assays eliminating dominant alpha-chain binding contributions. Moreover, the crystal structure of HGF beta determined at 2.53 A resolution provides a structural context for the mutagenesis data. The functional Met binding site is centered on the "active site region" including "triad" residues Gln(534) [c57], Asp(578) [c102], and Tyr(673) [c195] and neighboring "activation domain" residues Val(692), Pro(693), Gly(694), Arg(695), and Gly(696) [c214-c219]. Together they define a region that bears remarkable resemblance to substrate processing regions of serine proteases. Models of HGF-dependent Met receptor activation are discussed.

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.

Crystal structure of anticoagulant thrombin variant E217K provides insights into thrombin allostery

Thrombin is the ultimate protease of the blood clotting cascade and plays a major role in its own regulation. The ability of thrombin to exhibit both pro- and anti-coagulant properties has spawned efforts to turn thrombin into an anticoagulant for therapeutic purposes. This quest culminated in the identification of the E217K variant through scanning and saturation mutagenesis. The antithrombotic properties of E217K thrombin are derived from its inability to convert fibrinogen to a fibrin clot while maintaining its thrombomodulin-dependent ability to activate the anticoagulant protein C pathway. Here we describe the 2.5-A crystal structure of human E217K thrombin, which displays a dramatic restructuring of the geometry of the active site. Of particular interest is the repositioning of Glu-192, which hydrogen bonds to the catalytic Ser-195 and which results in the complete occlusion of the active site and the destruction of the oxyanion hole. Substrate binding pockets are further blocked by residues previously implicated in thrombin allostery. We have concluded that the E217K mutation causes the allosteric inactivation of thrombin by destabilizing the Na(+) binding site and that the structure thus may represent the Na(+)-free, catalytically inert "slow" form.

Residue Asp-189 Controls both Substrate Binding and the Monovalent Cation Specificity of Thrombin

Residue Asp-189 plays an important dual role in thrombin: it defines the primary specificity for Arg side chains and participates indirectly in the coordination of Na(+). The former role is shared by other proteases with trypsin-like specificity, whereas the latter is unique to Na(+)-activated proteases in blood coagulation and the complement system. Replacement of Asp-189 with Ala, Asn, Glu, and Ser drastically reduces the specificity toward substrates carrying Arg or Lys at P1, whereas it has little or no effect toward the hydrolysis of substrates carrying Phe at P1. These findings confirm the important role of Asp-189 in substrate recognition by trypsin-like proteases. The substitutions also affect significantly and unexpectedly the monovalent cation specificity of the enzyme. The Ala and Asn mutations abrogate monovalent cation binding, whereas the Ser and Glu mutations change the monovalent cation preference from Na(+) to the smaller cation Li(+) or to the larger cation Rb(+), respectively. The observation that a single amino acid substitution can alter the monovalent cation specificity of thrombin from Na(+) (Asp-189) to Li(+) (Ser-189) or Rb(+) (Glu-189) is unprecedented in the realm of monovalent cation-activated enzymes.

A Natural Prothrombin Mutant Reveals an Unexpected Influence of A-chain Structure on the Activity of Human α-Thrombin

We have recently identified in two unrelated patients with bleeding tendency a homozygous mutation causing a deletion of one of the two contiguous Lys(9)/Lys(10) residues in the A-chain of alpha-thrombin (DeltaK9). We used in vitro expression analysis to clarify the role of the deletion of Lys(9) or Lys(10) in the thrombin function. The k(cat)/K(m) value of the hydrolysis by DeltaK9 of the synthetic substrate Phe-Pip-Arg-p-nitroanilide (where Pip represents l-pipecolyl) and fibrinopeptide A was 18- and 60-fold lower, respectively, compared with wild type (WT). Interaction with antithrombin was also reduced in the mutant, the association rate being about 20-fold lower than in the WT thrombin. The sensitivity to sodium ion of DeltaK9 was found significantly attenuated compared with the WT form. DeltaK9 has a very weak platelet-activating capacity, attributed to a severely defective PAR1 interaction, whereas the binding to the platelet glycoprotein Ibalpha was unaffected. Likewise, the interaction with protein C was severely impaired, whereas interaction with thrombomodulin had a normal K(d) value. At variance with these findings, both low affinity (basic pancreatic trypsin inhibitor) and high affinity (N-alpha-[2-naphthylsulfonyl-glycyl]-4-amidinophenylalanine-piperidide) thrombin inhibitors displayed a better binding to DeltaK9 than to the WT form, indicating a better accommodation of these inhibitors into the catalytic pocket of DeltaK9. A molecular dynamics simulation of the DeltaK9 thrombin in full explicit water solvent provided support to the role of the A-chain in affecting conformation and catalytic properties of the B-chain, especially in some insertion loops of the enzyme, such as the 60-loop, as well as in the geometry of the catalytic triad residues.

Molecular mechanism of enzyme inhibition: prediction of the three-dimensional structure of the dimeric trypsin inhibitor from Leucaena leucocephala by homology modelling

Serine proteinase inhibitors are widely distributed in nature and inhibit the activity of enzymes like trypsin and chymotrypsin. These proteins interfere with the physiological processes such as germination, maturation and form the first line of defense against the attack of seed predator. The most thoroughly examined plant serine proteinase inhibitors are found in the species of the families Leguminosae, Graminae, and Solanaceae. Leucaena leucocephala belongs to the family Leguminosae. It is widely used both as an ornamental tree as well as cattle food. We have constructed a three-dimensional model of a serine proteinase inhibitor from L. leucocephala seeds (LTI) complexed with trypsin. The model was built based on its comparative homology with soybean trypsin inhibitor (STI) using the program, MODELLER6. The quality of the model was assessed stereochemically by PROCHECK. LTI shows structural features characteristic of the Kunitz type trypsin inhibitor and shows 39% residue identity with STI. LTI consists of 172 amino acid residues and is characterized by two disulfide bridges. The protein is a dimer with the two chains being linked by a disulfide bridge. Despite the high similarity in the overall tertiary structure, significant differences exist at the active site between STI and LTI. The present study aims at analyzing these interactions based on the available amino acid sequences and structural data. We have also studied some functional sites such as phosphorylation, myristoylation, which can influence the inhibitory activity or complexation with other molecules. Some of the differences observed at the active site and functional sites can explain the unique features of LTI.

Thrombin activation of factor XI on activated platelets requires the interaction of factor XI and platelet glycoprotein Ib alpha with thrombin anion-binding exosites I and II, respectively

Activation of factor XI (FXI) by thrombin on stimulated platelets plays a physiological role in hemostasis, providing additional thrombin generation required in cases of severe hemostatic challenge. Using a collection of 53 thrombin mutants, we identified 16 mutants with <50% of the wild-type thrombin FXI-activating activity in the presence of dextran sulfate. These mutants mapped to anion-binding exosite (ABE) I, ABE-II, the Na+-binding site, and the 50-insertion loop. Only the ABE-II mutants showed reduced binding to dextran sulfate-linked agarose. Selected thrombin mutants in ABE-I (R68A, R70A, and R73A), ABE-II (R98A, R245A, and K248A), the 50-insertion loop (W50A), and the Na+-binding site (E229A and R233A) with <10% of the wild-type activity also showed a markedly reduced ability to activate FXI in the presence of stimulated platelets. The ABE-I, 50-insertion loop, and Na+-binding site mutants had impaired binding to FXI, but normal binding to glycocalicin, the soluble form of glycoprotein Ibalpha (GPIb alpha). In contrast, the ABE-II mutants were defective in binding to glycocalicin, but displayed normal binding to FXI. Our data support a quaternary complex model of thrombin activation of FXI on stimulated platelets. Thrombin bound to one GPIb alpha molecule, via ABE-II on its posterior surface, is properly oriented for its activation of FXI bound to a neighboring GPI alpha molecule, via ABE-I on its anterior surface. GPIb alpha plays a critical role in the co-localization of thrombin and FXI and the resultant efficient activation of FXI.

Redesigning the monovalent cation specificity of an enzyme

Monovalent-cation-activated enzymes are abundantly represented in plants and in the animal world. Most of these enzymes are specifically activated by K+, whereas a few of them show preferential activation by Na+. The monovalent cation specificity of these enzymes remains elusive in molecular terms and has not been reengineered by site-directed mutagenesis. Here we demonstrate that thrombin, a Na+-activated allosteric enzyme involved in vertebrate blood clotting, can be converted into a K+-specific enzyme by redesigning a loop that shapes the entrance to the cation-binding site. The conversion, however, does not result into a K+-activated enzyme.

A common mutation, Arg457-->Gln, links prothrombin deficiencies in the Puerto Rican population

Five unrelated families with Puerto Rican ancestry were identified as having at least one member with bleeding due to a prothrombin deficiency. Genetic prothrombin deficiencies are extremely rare, but at the University of Puerto Rico Hemophilia Center, prothrombin deficiency is the third most common congenital coagulation factor deficiency. Because Puerto Rico is relatively isolated, there was a reasonable expectation of a founder effect. Prothrombin genes from probands and their parents were directly sequenced from PCR amplified exons using forward and reverse primers. Four novel prothrombin mutations were identified. The first, a G-->A substitution at DNA position 10150 predicting an Arg457-->Gln (R457Q) replacement, is common to all five families. In two of the families, the proband children are homozygous for R457Q. In the other three families, the probands are compound heterozygotes for R457Q and one of the other three mutations, which include another point mutation (gamma16Q), a deletion and a splice junction mutation. The two point mutations have been designated Puerto Rico I and Puerto Rico II. The crystal structure of alpha-thrombin predicts that the R457Q mutation removes a salt bridge that links the A- and B-chains of thrombin. The primary effect of this defect appears to be destabilization of the circulating prothrombin, creating a moderate hypoprothrombinemia. However, prothrombin antigen/activity ratios indicate a dysprothrombinemia as well, most likely due to the inability of R457Q prothrombin to activate fully to 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.

Roles of low specificity and cofactor interaction sites on thrombin during factor XIII activation. Competition for cofactor sites on thrombin determines its fate

Factor XIII is activated by thrombin, and this reaction is enhanced by the presence of fibrin(ogen). Using a substrate-based screening assay for factor XIII activity complemented by kinetic analysis of activation peptide cleavage, we show by using thrombin mutants of surface-exposed residues that Arg-178, Arg-180, Asp-183, Glu-229, Arg-233, and Trp-50 of thrombin are necessary for direct activation of factor XIII. These residues define a low specificity site known to be important also for both protein C activation and for inhibition of thrombin by antithrombin. The enhancing effect of fibrinogen occurs as a consequence of its conversion to fibrin and subsequent polymerization. Surface residues of thrombin further involved in high specificity fibrin-enhanced factor XIII activation were identified as His-66, Tyr-71, and Asn-74. These residues represent a distinct interaction site on thrombin (within exosite I) also employed by thrombomodulin in its cofactor-enhanced activation of protein C. In competition experiments, thrombomodulin inhibited fibrin-enhanced factor XIII activation. Based upon these and prior published results, we propose that the polymerization process forms a fibrin cofactor that acts to approximate thrombin and factor XIII bound to separate and complementary domains of fibrinogen. This enables enhanced factor XIII activation to be localized around the fibrin clot. We also conclude that proximity to and competition for cofactor interaction sites primarily directs the fate of thrombin.

Mutations in autolytic loop-2 and at Asp554 of human prothrombin that enhance protein C activation by meizothrombin

Thrombin acts on many protein substrates during the hemostatic process. Its specificity for these substrates is modulated through interactions at regions remote from the active site of the thrombin molecule, designated exosites. Exosite interactions can be with the substrate, cofactors such as thrombomodulin, or fragments from prothrombin. The relative activity of alpha-thrombin for fibrinogen is 10 times greater than that for protein C. However, the relative activity of meizothrombin for protein C is 14 times greater than that for fibrinogen. Modulation of thrombin specificity is linked to its Na(+)-binding site and residues in autolytic loop-2 that interact with the Na(+)-binding site. Recombinant prothrombins that yield recombinant meizothrombin (rMT) and rMT des-fragment 1 (rMT(desF1)) enable comparisons of the effects of mutations at the Na(+)-binding residue (Asp(554)) and deletion of loop-2 (Glu(466)-Thr(469)) on the relative activity of meizothrombin for several substrates. Hydrolysis of t-butoxycarbonyl-VPR-p-nitroanilide by alpha-thrombin, recombinant alpha-thrombin, or rMT(desF1) was almost identical, but that by rMT was only 40% of that by alpha-thrombin. Clotting of fibrinogen by rMT and rMT(desF1) was 12-16% of that by alpha-thrombin, as already known. Strikingly, however, although meizothrombins modified by substitution of Asp(554) with either Ala or Leu or by deletion of loop-2 had 6-8 and <1%, respectively, of the clotting activity of alpha-thrombin, the activity of these meizothrombins for protein C was increased to >10 times that of alpha-thrombin. It is proposed that interactions within thrombin that involve autolytic loop-2 and the Na(+)-binding site primarily enhance thrombin action on fibrinogen, but impair thrombin action on protein C.

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

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

Three-dimensional models of proteases involved in patterning of the Drosophila Embryo. Crucial role of predicted cation binding sites

Three-dimensional models of the catalytic domains of Nudel (Ndl), Gastrulation Defective (Gd), Snake (Snk), and Easter (Ea), and their complexes with substrate suggest a possible organization of the enzyme cascade controlling the dorsoventral fate of the fruit fly embryo. The models predict that Gd activates Snk, which in turn activates Ea. Gd can be activated either autoproteolytically or by Ndl. The three-dimensional models of each enzyme-substrate complex in the cascade rationalize existing mutagenesis data and the associated phenotypes. The models also predict unanticipated features like a Ca(2+) binding site in Ea and a Na(+) binding site in Ndl and Gd. These binding sites are likely to play a crucial role in vivo as suggested by mutant enzymes introduced into embryos as mRNAs. The mutations in Gd that eliminate Na(+) binding cause an apparent increase in activity, whereas mutations in Ea that abrogate Ca(2+) binding result in complete loss of activity. A mutation in Ea predicted to introduce Na(+) binding results in apparently increased activity with ventralization of the embryo, an effect not observed with wild-type Ea mRNA.

Reconstructing the binding site of factor Xa in trypsin reveals ligand-induced structural plasticity

In order to investigate issues of selectivity and specificity in protein-ligand interactions, we have undertaken the reconstruction of the binding pocket of human factor Xa in the structurally related rat trypsin by site-directed mutagenesis. Three sequential regions (the "99"-, the "175"- and the "190"- loops) were selected as representing the major structural differences between the ligand binding sites of the two enzymes. Wild-type rat trypsin and variants X99rT and X(99/175/190)rT were expressed in yeast, and analysed for their interaction with factor Xa and trypsin inhibitors. For most of the inhibitors studied, progressive loop replacement at the trypsin surface resulted in inhibitory profiles akin to factor Xa. Crystals of the variants were obtained in the presence of benzamidine (3), and could be soaked with the highly specific factor Xa inhibitor (1). Binding of the latter to X99rT results in a series of structural adaptations to the ligand, including the establishment of an "aromatic box" characteristic of factor Xa. In X(99/175/190)rT, introduction of the 175-loop results in a surprising re-orientation of the "intermediate helix", otherwise common to trypsin and factor Xa. The re-orientation is accompanied by an isomerisation of the Cys168-Cys182 disulphide bond, and burial of the critical Phe174 side-chain. In the presence of (1), a major re-organisation of the binding site takes place to yield a geometry identical to that of factor Xa. In all, binding of (1) to trypsin and its variants results in significant structural rearrangements, inducing a binding surface strongly reminiscent of factor Xa, against which the inhibitor was optimised. The structural data reveal a plasticity of the intermediate helix, which has been implicated in the functional cofactor dependency of many trypsin-like serine proteinases. This approach of grafting loops onto scaffolds of known related structures may serve to bridge the gap between structural genomics and drug design.

Interconversion of ATP binding and conformational free energies by tryptophanyl-tRNA synthetase: structures of ATP bound to open and closed, pre-transition-state conformations

Binding ATP to tryptophanyl-tRNA synthetase (TrpRS) in a catalytically competent configuration for amino acid activation destabilizes the enzyme structure prior to forming the transition state. This conclusion follows from monitoring the titration of TrpRS with ATP by small angle solution X-ray scattering, enzyme activity, and crystal structures. ATP induces a significantly smaller radius of gyration at pH=7 with a transition midpoint at approximately 8mM. A non-reciprocal dependence of Trp and ATP dissociation constants on concentrations of the second substrate show that Trp binding enhances affinity for ATP, while the affinity for Trp falls with the square of the [ATP] over the same concentration range ( approximately 5mM) that induces the more compact conformation. Two distinct TrpRS:ATP structures have been solved, a high-affinity complex grown with 1mM ATP and a low-affinity complex grown at 10mM ATP. The former is isomorphous with unliganded TrpRS and the Trp complex from monoclinic crystals. Reacting groups of the two individually-bound substrates are separated by 6.7A. Although it lacks tryptophan, the low-affinity complex has a closed conformation similar to that observed in the presence of both ATP and Trp analogs such as indolmycin, and resembles a complex previously postulated to form in the closely-related TyrRS upon induced-fit active-site assembly, just prior to catalysis. Titration of TrpRS with ATP therefore successively produces structurally distinct high- and low-affinity ATP-bound states. The higher quality X-ray data for the closed ATP complex (2.2A) provide new structural details likely related to catalysis, including an extension of the KMSKS loop that engages the second lysine and serine residues, K195 and S196, with the alpha and gamma-phosphates; interactions of the K111 side-chain with the gamma-phosphate; and a water molecule bridging the consensus sequence residue T15 to the beta-phosphate. Induced-fit therefore strengthens active-site interactions with ATP, substantially intensifying the interaction of the KMSKS loop with the leaving PP(i) group. Formation of this conformation in the absence of a Trp analog implies that ATP is a key allosteric effector for TrpRS. The paradoxical requirement for high [ATP] implies that Gibbs binding free energy is stored in an unfavorable protein conformation and can then be recovered for useful purposes, including catalysis in the case of TrpRS.

Probing the hirudin-thrombin interaction by incorporation of noncoded amino acids and molecular dynamics simulation

Thrombin is a primary target for the development of novel anticoagulants, since it plays two important and opposite roles in hemostasis: procoagulant and anticoagulant. All thrombin functions are influenced by Na+ binding, which triggers the transition of this enzyme from an anticoagulant (slow) form to a procoagulant (fast) form. In previous studies, we have conveniently produced by chemical synthesis analogues of the N-terminal fragment 1-47 of hirudin HM2 containing noncoded amino acids and displaying up to approximately 2700-fold more potent antithrombin activity, comparable to that of full-length hirudin. In the work presented here, we have exploited the versatility of chemical synthesis to probe the structural and energetic properties of the S3 site of thrombin through perturbations introduced in the structure of hirudin fragment 1-47. In particular, we have investigated the effects of systematic replacement of Tyr3 with noncoded amino acids retaining the aromatic nucleus of Tyr, as well as similar hydrophobic and steric properties, but possessing different electronic (e.g., p-fluoro-, p-iodo-, or p-nitro-Phe), charge (p-aminomethyl-Phe), or conformational (homo-Phe) properties. Our results indicate that the affinity of fragment 1-47 for thrombin is proportional to the desolvation free energy change upon complex formation, and is inversely related to the electric dipole moment of the amino acid side chain at position 3 of hirudin. In this study, we have also identified the key features that are responsible for the preferential binding of hirudin to the procoagulant (fast) form of thrombin. Strikingly, shaving at position 3, by Tyr --> Ala exchange, abolishes the differences in the affinity for thrombin allosteric forms, whereas a bulkier side chain (e.g., beta-naphthylalanine) improves binding preferentially to the fast form. These results provide strong, albeit indirect, evidence that the procoagulant (fast) form of thrombin is in a more open and accessible conformation with respect to the less forgiving structure it acquires in the slow form. This view is also supported by the results of molecular dynamics simulations conducted for 18 ns on free thrombin in full explicit water, showing that after approximately 5 ns thrombin undergoes a significant conformational transition, from a more open conformation (which we propose can be related to the fast form) to a more compact and closed one (which we propose can be related to the slow form). This transition mainly involves the Trp148 and Trp60D loop, the S3 site, and the fibrinogen binding site, whereas the S1 site, the Na+-binding site, and the catalytic pocket remain essentially unchanged. In particular, our data indicate that the S3 site of the enzyme is less accessible to water in the putative slow form. This structural picture provides a reasonable molecular explanation for the fact that physiological substrates related to the procoagulant activity of thrombin (fibrinogen, thrombin receptor 1, and factor XIII) orient a bulky side chain into the S3 site of the enzyme. Taken together, our results can have important implications for the design of novel thrombin inhibitors, of practical utility in the treatment of coagulative disorders.

Synthesis and physical characterization of a P1 arginine combinatorial library, and its application to the determination of the substrate specificity of serine peptidases

Serine peptidases are a large, well-studied, and medically important class of peptidases. Despite the attention these enzymes have received, details concerning the substrate specificity of even some of the best known enzymes in this class are lacking. One approach to rapidly characterizing substrate specificity for peptidases is the use of positional scanning combinatorial substrate libraries. We recently synthesized such a library for enzymes with a preference for arginine at P1 and demonstrated the use of this library with thrombin (Edwards et al. Bioorg. Med. Chem. Lett. 2000, 10, 2291). In the present work, we extend these studies by demonstrating good agreement between the theroretical and measured content of portions of this library and by showing that the library permits rapid characterization of the substrate specificity of additional SA clan serine peptidases including factor Xa, tryptase, and trypsin. These results were consistent both with cleavage sites in natural substrates and cleavage of commercially available synthetic substrates. We also demonstrate that pH or salt concentration have a quantitative effect on the rate of cleavage of the pooled library substrates but that correct prediction of optimal substrates for the enzymes studied appeared to be independent of these parameters. These studies provide new substrate specificity data on an important class of peptidases and are the first to provide physical characterization of a peptidase substrate library.

Ser(214) is crucial for substrate binding to serine proteases

Highly conserved amino acids that form crucial structural elements of the catalytic apparatus can be used to account for the evolutionary history of serine proteases and the cascades into which they are organized. One such evolutionary marker in chymotrypsin-like proteases is Ser(214), located adjacent to the active site and forming part of the primary specificity pocket. Here we report the mutation of Ser(214) in thrombin to Ala, Thr, Cys, Asp, Glu, and Lys. None of the mutants seriously compromises active site catalytic function as measured by the kinetic parameter k(cat). However, the least conservative mutations result in large increases in K(m) because of lower rates of substrate diffusion into the active site. Therefore, the role of Ser(214) is to promote the productive formation of the enzyme-substrate complex. The S214C mutant is catalytically inactive, which suggests that during evolution the TCN-->AGY codon transitions for Ser(214) occurred through Thr intermediates.

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.

Structural requirements for the activation of human factor VIII by thrombin

The coagulation factors V (FV) and VIII (FVIII) are important at sites of vascular injury for the amplification of the clotting cascade. Natural variants of these factors frequently lead to severe bleeding disorders. To understand the mechanisms of activation of FVIII by thrombin, we used a bank of mutant thrombins to define residues important for its activation. From the initial screening of 53 mutant thrombins for the activation of human recombinant FVIII, we mapped thrombin mutants with 50% or less activity to anion-binding exosite-I (Lys21Ala, His66Ala, Lys65Ala, Arg68Ala, Arg70Ala, and Tyr71Ala) and anion-binding exosite-II (Arg98Ala), the Na(+)-binding site (Glu229Ala, Arg233Ala, Asp234Ala, and Asp193Ala/Lys196Ala), and the 50-insertion loop (Trp50Ala), which were similar to our results for the activation of FV. The role of these residues for cleavage at Arg372 and Arg1689 was investigated using plasma FVIII. Anion-binding exosite-I appears to be important for cleavage at both sites, whereas the anion-binding exosite-II residue Arg98Ala is important for cleavage at Arg372 alone. The Glu229Ala mutant, which contributes to the Na(+)-binding site, and the 50-insertion loop mutant W50A have severely impaired cleavage at Arg372 and Arg1689. This suggests that the integrity of the active site and the Na(+)-bound form of thrombin are important for its procoagulant activity against FVIII. Detailed mutagenic analysis of thrombin can assist in understanding the pathogenesis of bleeding disorders and may lead to the rational design of selective thrombin inhibitors.

Prothrombinase assembly and S1 site occupation restore the catalytic activity of FXa impaired by mutation at the sodium-binding site

Two loop segments (183-189 and 221-225) in the protease domain of factor Xa contribute to the formation of a Na(+)-binding site. Studies with factor Xa indicate that binding of a single Na(+) ion to this site influences its activity by altering the S1 specificity site, and substitution of Tyr(225) with Pro diminishes sensitivity to Na(+). Using full-length factor Xa(Y225P), the allosteric relationship between the Na(+) site and other structural determinants in factor Xa and prothrombinase was investigated. Direct binding and kinetic measurements with probes that target the S1 specificity pocket indicate that assembly of the mutant in prothrombinase corrected the impaired binding of these probes observed with free factor Xa(Y225P). This appears to result from the apparent allosteric linkage between the factor Va, S1, and Na(+)-binding sites, since binding of the cofactor to membrane-bound factor Xa(Y225P) enhances binding at the S1 site and vice versa. Additional studies revealed that the internal salt bridge (Ile(16)-Asp(194)) of factor Xa(Y225P) is partially destabilized, a process that is reversible upon occupation of the S1 site. The data establish that alterations at the factor Xa Na(+)-binding site shift the zymogen-protease equilibrium to a more zymogen-like state, and as a consequence binding of S1-directed probes and factor Va are adversely affected. Therefore, the zymogen-like characteristics of factor Xa(Y225P) have allowed for the apparent allosteric linkage between the S1, factor Va, and Na(+) sites to become evident and has provided insight into the structural transitions which accompany the conversion of factor X to factor Xa.

Molecular and functional characterization of a natural homozygous Arg67His mutation in the prothrombin gene of a patient with a severe procoagulant defect contrasting with a mild hemorrhagic phenotype

In a patient who presented with a severe coagulation deficiency in plasma contrasting with a very mild hemorrhagic diathesis a homozygous Arg67His mutation was identified in the prothrombin gene. Wild-type (factor IIa [FIIa]-WT) and mutant Arg67His thrombin (FIIa-MT67) had similar amidolytic activity. By contrast, the k(cat)/K(m) value of fibrinopeptide A hydrolysis by FIIa-WT and FIIa-MT67 was equal to 2.1 x 10(7) M(-1)s(-1) and 9 x 10(5) M(-1)s(-1). Decreased activation of protein C (PC) correlated with the 33-fold decreased binding affinity for thrombomodulin (TM; K(d) = 65.3 nM vs 2.1 nM, in FIIa-MT67 and in FIIa-WT, respectively). In contrast, hydrolysis of PC in the absence of TM was normal. The Arg67His mutation had a dramatic effect on the cleavage of protease-activated G protein-coupled receptor 1 (PAR-1) 38-60 peptide (k(cat/)K(m) = 4 x 10(7) M(-1)s(-1) to 1.2 x 10(6) M(-1)s(-1)). FIIa-MT67 showed a weaker platelet activating capacity, attributed to a defective PAR-1 interaction, whereas the interaction with glycoprotein Ib was normal. A drastic decrease (up to 500-fold) of the second-order rate constant pertaining to heparin cofactor II (HCII) interaction, especially in the presence of dermatan sulfate, was found for the FIIa-MT67 compared with FIIa-WT, suggesting a severe impairment of thrombin inhibition by HCII in vivo. Finally, the Arg67His mutation was associated with a 5-fold decrease of prothrombin activation by the factor Xa-factor Va complex, perhaps through impairment of the prothrombin-factor Va interaction. These experiments show that the Arg67His substitution affects drastically both the procoagulant and the anticoagulant functions of thrombin as well as its inhibition by HCII. The mild hemorrhagic phenotype might be explained by abnormalities that ultimately counterbalance each other.

Salt-induced stabilization of apoflavodoxin at neutral pH is mediated through cation-specific effects

Electrostatic contributions to the conformational stability of apoflavodoxin were studied by measurement of the proton and salt-linked stability of this highly acidic protein with urea and temperature denaturation. Structure-based calculations of electrostatic Gibbs free energy were performed in parallel over a range of pH values and salt concentrations with an empirical continuum method. The stability of apoflavodoxin was higher near the isoelectric point (pH 4) than at neutral pH. This behavior was captured quantitatively by the structure-based calculations. In addition, the calculations showed that increasing salt concentration in the range of 0 to 500 mM stabilized the protein, which was confirmed experimentally. The effects of salts on stability were strongly dependent on cationic species: K(+), Na(+), Ca(2+), and Mg(2+) exerted similar effects, much different from the effect measured in the presence of the bulky choline cation. Thus cations bind weakly to the negatively charged surface of apoflavodoxin. The similar magnitude of the effects exerted by different cations indicates that their hydration shells are not disrupted significantly by interactions with the protein. Site-directed mutagenesis of selected residues and the analysis of truncation variants indicate that cation binding is not site-specific and that the cation-binding regions are located in the central region of the protein sequence. Three-state analysis of the thermal denaturation indicates that the equilibrium intermediate populated during thermal unfolding is competent to bind cations. The unusual increase in the stability of apoflavodoxin at neutral pH affected by salts is likely to be a common property among highly acidic proteins.

Factorising ligand affinity: a combined thermodynamic and crystallographic study of trypsin and thrombin inhibition

The binding of a series of low molecular weight ligands towards trypsin and thrombin has been studied by isothermal titration calorimetry and protein crystallography. In a series of congeneric ligands, surprising changes of protonation states occur and are overlaid on the binding process. They result from induced pK(a) shifts depending on the local environment experienced by the ligand and protein functional groups in the complex (induced dielectric fit). They involve additional heat effects that must be corrected before any conclusion on the binding enthalpy (DeltaH) and entropy (DeltaS) can be drawn. After correction, trends in both contributions can be interpreted in structural terms with respect to the hydrogen bond inventory or residual ligand motions. For all inhibitors studied, a strong negative heat capacity change (DeltaC(p)) is detected, thus binding becomes more exothermic and entropically less favourable with increasing temperature. Due to a mutual compensation, Gibbs free energy remains virtually unchanged. The strong negative DeltaC(p) value cannot solely be explained by the removal of hydrophobic surface portions of the protein or ligand from water exposure. Additional contributions must be considered, presumably arising from modulations of the local water structure, changes in vibrational modes or other ordering parameters. For thrombin, smaller negative DeltaC(p) values are observed for ligand binding in the presence of sodium ions compared to the other alkali ions, probably due to stabilising effects on the protein or changes in the bound water structure.

Determinants of thrombin specificity

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

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

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

Using surface-bound rubidium ions for protein phasing

Rubidium is a monovalent metal that can be used as a counterion in protein solutions. X-ray anomalous scattering from rubidium ions bound to the protein surface was used for phasing of the crystal structure of the hsp60 apical domain from Thermus thermophilus. Multiple-wavelength anomalous dispersion (MAD) data were collected from a crystal obtained from a solution containing 0.2 M rubidium salt. One molecule of protein (147 amino acids) binds one well ordered and one poorly ordered Rb atom. Phases calculated with the program SHARP were sufficient for automatic tracing and side-chain assignment using the program ARP/wARP. The data show that bound rubidium ions can be used to determine protein structures and to study the interaction of monovalent metal ions with proteins and other macromolecules.

Replacement of thrombin residue G184 with Lys or Arg fails to mimic Na+ binding

Na+ binding to thrombin enhances the catalytic activity toward numerous synthetic and natural substrates. The bound Na+ is located in a solvent channel 16 A away from the catalytic triad, and connects with D189 in the S1 site through an intervening water molecule. Molecular modeling indicates that the G184K substitution in thrombin positions the protonated epsilon-amino group of the Lys side-chain to replace the bound Na+. Likewise, the G184R substitution positions the guanidinium group of the longer Arg side-chain to replace both the bound Na+ and the connecting water molecule to D189. We explored whether the G184K or G184R substitution would replace the bound Na+ and yield a thrombin derivative stabilized in the highly active fast form. Both the G184K and G184R mutants lost sensitivity to monovalent cations, as expected, but their activity toward a chromogenic substrate was compromised up to 200-fold as a result of impaired diffusion into the S1 site and decreased deacylation rate. Interestingly, both G184K and G184R substitutions compromised cleavage of procoagulant substrates fibrinogen and PAR1 more than that of the anticoagulant substrate protein C. These findings demonstrate that Na+ binding to thrombin is difficult to mimic functionally with residue side-chains, in analogy with results from other systems.

Heparin enhances the specificity of antithrombin for thrombin and factor Xa independent of the reactive center loop sequence. Evidence for an exosite determinant of factor Xa specificity in heparin-activated antithrombin

Heparin activates the primary serpin inhibitor of blood clotting proteinases, antithrombin, both by an allosteric conformational change mechanism that specifically enhances factor Xa inactivation and by a ternary complex bridging mechanism that promotes the inactivation of thrombin and other target proteinases. To determine whether the factor Xa specificity of allosterically activated antithrombin is encoded in the reactive center loop sequence, we attempted to switch this specificity by mutating the P6-P3' proteinase binding sequence excluding P1-P1' to a more optimal thrombin recognition sequence. Evaluation of 12 such antithrombin variants showed that the thrombin specificity of the serpin allosterically activated by a heparin pentasaccharide could be enhanced as much as 55-fold by changing P3, P2, and P2' residues to a consensus thrombin recognition sequence. However, at most 9-fold of the enhanced thrombin specificity was due to allosteric activation, the remainder being realized without activation. Moreover, thrombin specificity enhancements were attenuated to at most 5-fold with a bridging heparin activator. Surprisingly, none of the reactive center loop mutations greatly affected the factor Xa specificity of the unactivated serpin or the several hundred-fold enhancement in factor Xa specificity due to activation by pentasaccharide or bridging heparins. Together, these results suggest that the specificity of both native and heparin-activated antithrombin for thrombin and factor Xa is only weakly dependent on the P6-P3' residues flanking the primary P1-P1' recognition site in the serpin-reactive center loop and that heparin enhances serpin specificity for both enzymes through secondary interaction sites outside the P6-P3' region, which involve a bridging site on heparin in the case of thrombin and a previously unrecognized exosite on antithrombin in the case of factor Xa.

Thermodynamics of Na+ binding to coagulation serine proteases

The sodium binding to serine proteases triggers a conformational change in the proteins that enhances the catalytic activity of the enzymes. The interaction of the cation with the protein is mediated by the hydrogen-bonding network of water molecules that embed the Na+ site. We pointed out the crucial role of the insertion loop 186a-d and the I16-D194 ion pair in the stabilization of sodium binding pocket in thrombin. This paper contributes to better explain the molecular mechanism of sodium binding for different serine proteases leading to the identification of the structural changes necessary to engineer a functional Na+ site and regulate catalytic activity in serine proteases.

Platelet glycoprotein Ib alpha binds to thrombin anion-binding exosite II inducing allosteric changes in the activity of thrombin

The glycoprotein (GP) Ib-IX complex is a platelet surface receptor that binds thrombin as one of its ligands, although the biological significance of thrombin interaction remains unclear. In this study we have used several approaches to investigate the GPIb alpha-thrombin interaction in more detail and to study its effect on the thrombin-induced elaboration of fibrin. We found that both glycocalicin and the amino-terminal fragment of GPIb alpha reduced the release of fibrinopeptide A from fibrinogen by about 50% by a noncompetitive allosteric mechanism. Similarly, GPIb alpha caused in thrombin an allosteric reduction in the rate of turnover of the small peptide substrate d-Phe-Pro-Arg-pNA. The K(d) for the glycocalicin-thrombin interaction was 1 microm at physiological ionic strength but was highly salt-dependent, decreasing to 0.19 microm at 100 mm NaCl (Gamma(salt) = -4.2). The salt dependence was characteristic of other thrombin ligands that bind to exosite II of this enzyme, and we confirmed this as the GPIb alpha-binding site on thrombin by using thrombin mutants and by competition binding studies. R68E or R70E mutations in exosite I of thrombin had little effect on its interaction with GPIb alpha. Both the allosteric inhibition of fibrinogen turnover caused by GPIb alpha binding to these mutants, and the K(d) values for their interactions with GPIb alpha were similar to those of wild-type thrombin. In contrast, R89E and K248E mutations in exosite II of thrombin markedly increased the K(d) values for the interactions of these thrombin mutants with GPIb alpha by 10- and 25-fold, respectively. Finally, we demonstrated that low molecular weight heparin (which binds to thrombin exosite II) but not hirugen (residues 54-65 of hirudin, which binds to exosite I of thrombin) inhibited thrombin binding to GPIb alpha. These data demonstrate that GPIb alpha binds to thrombin exosite II and in so doing causes a conformational change in the active site of thrombin by an allosteric mechanism that alters the accessibility of both its natural substrate, fibrinogen, and the small peptidyl substrate d-Phe-Pro-Arg-pNA.

Preparation of anhydrothrombin and characterization of its interaction with natural thrombin substrates

Thrombin is a serine proteinase that plays a key role in thrombosis and haemostasis through its interaction with several coagulation factors. Anhydrothrombin was prepared from PMSF-inactivated thrombin under alkaline conditions, and the folded anhydrothrombin was successfully recovered after dialysis in the presence of glycerol. Anhydro-derivatives of factor Xa, factor VIIa and activated protein C could also be prepared essentially by the same procedure. Anhydrothrombin retained affinity for various natural substrates of thrombin, including fibrinogen, factor VIII, factor XIII and protein C. In addition, these proteins were bound to anhydrothrombin-agarose in a reversible manner. The K(d) values for factor VIII, fibrinogen, factor XIII and protein C were 1.2x10(-8), 4.4x10(-8), 2.8x10(-7) and 8.1x10(-5) M, respectively. Thus thrombin substrates known to interact with the exosite I of thrombin demonstrated high affinity for anhydrothrombin. Furthermore, in the presence of Na+, substantial enhancement of the association rate constant (k(ass)) was observed for interactions of fibrinogen and factor VIII with anhydrothrombin. These results suggest that anhydrothrombin is useful in the purification of thrombin substrate proteins as well as in the investigation of detailed interactions between thrombin and these substrates in their activation or degradation processes.

Amino acids involved in sodium interaction of murine type II Na(+)-P(i) cotransporters expressed in Xenopus oocytes

Type IIa and IIb Na+-Pi cotransporters are highly conserved proteins expressed in brush border membranes of proximal tubules and small intestine, respectively. The kinetics of IIa and IIb differ significantly: type IIb is saturated at lower concentrations of Na+ and Pi. To define the domain responsible for the difference in Na+ affinity we constructed several mouse IIa-IIb chimeras as well as site-directed mutagenized cotransporters. Pi uptake activity was determined after injection of cRNAs into Xenopus laevis oocytes. From the chimera experiments we concluded that the domain containing part of the second intracellular loop, the fifth transmembrane domain (TD) and part of the third extracellular loop determines the specific Na+ activation properties for both types of cotransporter. Within this domain only a few residues located in the fifth TD are not conserved between type IIa and IIb. Site-directed mutagenesis on non-conserved residues was performed. Substitution of F402 of IIa by the corresponding L418 from IIb yielded a cotransporter that behaved like the IIb. On the other hand, substitution of the specific L418 of IIb by the corresponding F402 of IIa produced a cotransporter with a Na+ activation similar to IIa. (Single letter amino acid nomenclature is used throughout the paper.) These data suggest that the specific Na+ activation properties exhibited by type IIa and type IIb Na+-Pi cotransporters are at least in part due to the presence of a specific amino acid (F402 in IIa, and L418 in IIb) within the fifth TD of the protein.

Enhancement of chymotrypsin-inhibitor/substrate interactions by 3 M NaCl

A series of 16 bovine pancreatic trypsin inhibitor variants mutated at the P(1) position of the binding loop and seven tetrapeptide p-nitroanilide (pNa) substrates of the general formula: suc-Ala-Ala-Pro-Aaa-pNa (where Aaa denotes either: Phe, Arg, Lys, Leu, Met, Nva, Nle) were used to investigate the influence of high salt concentration on the activity of bovine chymotrypsin. The increase of the association constant (K(a)) and the specificity index (k(cat)/K(m)) in the presence of 3 M NaCl highly depends on the chemical nature of the residue at the P(1) position. The highest increase was observed for inhibitors/substrates containing the basic side chains at this site. Surprisingly, for the remaining 13 residues the observed salt effect is not correlated with any side chain properties. In particular, there is a lack of correlation between the accessible non-polar surface area and the magnitude of the salt effect. It suggests that salt-induced increase of the K(a) and k(cat)/K(m) values is not caused by the enhancement of the hydrophobic interactions in chymotrypsin-inhibitor/substrate complex. Moreover, the increase of the K(a) and k(cat)/K(m) values occurs only in the presence of Na(+) ions, while K(+) and Li(+) ions do not change the activity of chymotrypsin. Additionally, the activities of two other proteinases: bovine trypsin and Streptomyces griseus proteinase B were tested in the presence of 3 M NaCl using their specific substrates. The activity of both enzymes was almost not affected by the presence of high NaCl concentration.

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.