Structural Mapping of the Active Site Specificity Determinants of Human Tissue-type Plasminogen Activator
IMPLICATIONS FOR THE DESIGN OF LOW MOLECULAR WEIGHT SUBSTRATES AND INHIBITORS*

(Received for publication, April 2, 1997, and in revised form, May 12, 1997)

Martin Renatus Dagger §, Wolfram Bode Dagger §, Robert Huber Dagger , Jörg Stürzebecher , Dagmar Prasa , Stephan Fischer par , Ulrich Kohnert par and Milton T. Stubbs Dagger **

From the Dagger  Max Planck Institute of Biochemistry, Department of Structural Research, D-82152 Martinsried,  Center for Vascular Biology and Medicine, Nordhäuserstr. 78, University of Jena, D-99089 Erfurt, and par  Biochemical Research Center, Boehringer Mannheim GmbH, Nonnenwald 2, D-82372 Penzberg, Germany

ABSTRACT
INTRODUCTION
EXPERIMENTAL PROCEDURES
RESULTS
DISCUSSION
FOOTNOTES
ACKNOWLEDGEMENTS
REFERENCES


ABSTRACT

The recent structure determination of the catalytic domain of tissue-type plasminogen activator (tPA) suggested residue Arg174 could play a role in P3/P4 substrate specificity. Six synthetic chromogenic tPA substrates of the type R-Xaa-Gly-Arg-p-nitroanilide, in which R is an N-terminal protection group, were synthesized to test this property. Although changing the residue Xaa (in its L or D form) at position P3 from the hydrophobic Phe to an acidic residue, Asp or Glu, gave no improvement in catalytic efficiency, comparative analysis of the substrates indicated a preference for an acidic substituent occupying the S3 site when the S4 site contains a hydrophobic or basic moiety. The 2.9 Å structure determination of the catalytic domain of human tPA in complex with the bis-benzamidine inhibitor 2,7-bis-(4-amidinobenzylidene)-cycloheptan-1-one reveals a three-site interaction, salt bridge formation of the proximal amidino group of the inhibitor with Asp189 in the primary specificity pocket, extensive hydrophobic surface burial, and a weak electrostatic interaction between the distal amidino group of the inhibitor and two carbonyl oxygens of the protein. The latter position was previously occupied by the guanidino group of Arg174, which swings out to form the western edge of the S3 pocket. These data suggest that the side chain of Arg174 is flexible, and does not play a major role in the S4 specificity of tPA. On the other hand, this residue would modulate S3 specificity, and may be exploited to fine tune the specificity and selectivity of tPA substrates and inhibitors.


INTRODUCTION

The 60-kDa multidomain tissue-type plasminogen activator (tPA)1 catalyzes the conversion of the zymogen plasminogen into the active enzyme plasmin, the rate-limiting step in the endogenous fibrinolytic cascade (1). Recombinant tPA is used therapeutically as a fibrinolytic agent in the treatment of acute myocardial infarction and pulmonary embolism. Its fibrin dependent activity has attracted particular interest as it allows targeting of enzymatic activity to its natural substrate plasminogen (2, 3).

The recent structure determination of the catalytic domain of tPA (4) in complex with benzamidine (hereafter termed b-tPA) revealed a strong structural similarity to other trypsin-like serine proteases. 187 residues of tPA have topological equivalents in alpha -chymotrypsin, forming the basis for the chymotrypsinogen numbering of the tPA catalytic domain used in this report (4). Compared with chymotrypsin, tPA contains a deletion of three residues at the C terminus of the A-chain, two single residue deletions in the B-chain, and insertions at six different positions totaling 24 residues. Five of these insertions are noteworthy: one forms a helix, the "intermediate helix" (Pro164-Leu171),2 and four form surface loops, referred to as the 37-loop (Lys36-Arg39), 60-loop (Phe59-His62), 110-loop (Ser110-Cys111), and 186-loop (Asp185-Ala186G) (see Fig. 1). Most of these surface loops cluster around the active site cleft and are involved in specific interactions of the molecule with its substrates and inhibitors.


Fig. 1. Stereo ribbon plot of the catalytic domain of tc-tPA in complex with the bis-benzamidine inhibitor shown in Fig. 2. The active site residues Ser195, His57, Asp102, Asp189 at the bottom of the S1 specificity pocket and Tyr99, Trp215, Arg174 bordering the S4 pocket are highlighted.
[View Larger Version of this Image (41K GIF file)]

tPA exhibits remarkable specificity, a single bond of plasminogen (Arg560-Val561) is the only known substrate cleavage site of tPA in vivo. While the role of the specificity pocket S1 of tPA (and trypsin-like serine proteinases in general) in determining its preference for P1 Arg residues is known, both the location and mechanism of additional specificity determinants remain uncertain (5). In contrast to its high in vivo activity and specificity toward plasminogen, tPA shows low activity toward small peptides, peptide substrates modelled according to the cleavage sequence of plasminogen are cleaved by trypsin with much higher activity (>104-fold) (5). Optimization of the most labile peptide using phage display (6) emphasized the limited reactivity of tPA compared with trypsin and its high sensitivity to changes in residues at distinct positions. This suggests that tPA recognizes complex or multiple elements on the surface of plasminogen distant from the cleavage site. This might not only represent an increased binding affinity, but also a subtle alteration of the catalytic machinery of the protease domain of tPA. On the other hand, tPA is able to cleave the more labile amide bond of p-nitroanilide substrates, for example MeSO2-D-HHT-Gly-Arg-pNA (spectrozyme tPA) and Boc-Leu-Gly-Arg-pNA with a catalytic efficiency (kcat/Km) only 6.5-18-fold lower than that of trypsin (5).

Few protein or synthetic tPA-specific inhibitors are known. Plasminogenolytic activity in vivo is controlled by the highly specific and fast acting serpin PAI-1 (7). Structural aspects of the interaction of PAI-1 and tPA have been described elsewhere (4, 8-11). One of the few non-serpin proteinase inhibitors that block tPA with high affinity is the Erythrina trypsin inhibitor (ETI) (12-14), for which docking studies have been described recently (4). Synthetic inhibitors for enzymes of the trypsin family such as tPA have largely been based on Arg and Lys derivatives and the structurally related benzamidines. Remarkable variations are seen in the structure-activity relationships for the inhibition of trypsin family members with benzamidine-type inhibitors (15).

The three-dimensional structure of tPA (4) reveals some special features of the active site region, suggesting a strategy for obtaining more selective p-nitroanilide substrates and inhibitors. Adjacent to the catalytic triad, His57 (322), Asp102 (371), Ser195 (478) (the tPA numbering for the residues of the catalytic triad are given in parentheses), and the oxyanion hole, the specificity pocket S1 is bordered by the segment Ile213-Cys220, including the residue Trp215 (the entrance of the pocket), Asp189-Ser195 (the base of the pocket), Pro225-Tyr229 (the back of the pocket), and the disulfide bridge Cys191-Cys220 (the south of the pocket). Asp189 at the base of the S1 pocket allows salt bridge formation with basic residues. tPA shows a preference for S1 Arg over Lys, explained by the size of the S1 pocket and its lack of other polar residues (4, 16). The S2 subsite of tPA is restricted by the imposing phenolic side chain of Tyr99, leading to a preference for small residues such as Gly at position P2 as seen for factor Xa (17). Next to the S2 site is a hydrophobic surface, the S4 pocket, with Trp215 at the base and Tyr99 to the east. The pocket is delimited to the west by the side chain of Arg174 and to the north by the main chain atoms of Asp97, Thr98, and Thr175.

The side chain of Arg174 in the b-tPA structure (4), which points toward the S4 pocket and partially occupies it, has drawn our attention. Its position suggested that it could influence the binding of substrates containing acidic P4 residues for all L-amino acid peptides, or acidic P3 peptides when the P3 residue is in the D-configuration (18, 19). In the following, we shall use the terms "S3" and "S4" to denote the sites on tPA that would be occupied by a natural extended substrate. For a peptide of the form D-Xaa-Gly-Arg, the P3 side chain would occupy the S4 site, while an amino-terminal substituent would occupy the S3 site. We have constructed and synthesized a series of substrates of the type R-Xaa-Gly-Arg-pNA. In addition, we have solved the structure of the catalytic domain of two-chain tPA (Fig. 1) in complex with 2,7-bis-(4-amidinobenzylidene)-cycloheptan-1-one, a bis-benzamidine derived inhibitor (see Fig. 2).


Fig. 2. Chemical structure of the 2,7-bis-(4-amidinobenzylidene)-cycloheptan-1-one inhibitor. The inhibitor is represented in its trans/trans conformation as found in the tPA/bis-benzamidine complex. The Ki values for the inhibition of full length sc-tPA and for factor Xa are 0.5 µM and 0.013 µM respectively.
[View Larger Version of this Image (6K GIF file)]


EXPERIMENTAL PROCEDURES

Synthesis of the Inhibitor and Substrates

The inhibitor 2,7-bis-(4-amidinobenzylidene)-cycloheptan-1-one (Pefabloc Xa/tPA, Fig. 2) was kindly provided by the Pentapharm Ltd., Basel, Switzerland. The inhibitory activity was characterized as described elsewhere (20-22). Enzyme assays were carried out following standard procedures. Six substrates of the type Xaa-Gly-Arg-pNA with D- and L-Phe, D- and L-Glu and D- and L-Asp as Xaa residue where synthesized as described previously (23). The substrates Boc-Gly-Gly-Arg-pNA, Cbo-L-Val-Gly-Arg-pNA (chromozyme Try), Cmo-D-Val-Gly-Arg-pNA, Cmo-D-Nle-Gly-Arg-pNA (chromozyme X), Cmo-D-Leu-Gly-Arg-pNA, MeSO2-D-Phe-Gly-Arg-pNA (chromozyme tPA), H-CHG-Gly-Arg-pNA, H-D-CHG-Gly-Arg-pNA, Cmo-D-CHG-Gly-Arg-pNA (spectrozyme Xa), H-D-HHT-Gly-Arg-pNA (chromozyme XII), Cmo-D-HHT-Gly-Arg-pNA (spectrozyme LAL), Cmo-D-CHA-Gly-Arg-pNA (pefachrome Xa) and H-D-Lys-Gly-Arg-pNA were kindly provided by Lars Svendsen, Pentapharm Ltd., Basel, Switzerland. Cbo-D-Arg-Gly-Arg-pNA (S-2765) and pGlu-Gly-Arg-pNA (S-2444) were obtained from Hemochrom Diagnostica GmbH, Essen, Germany.

Hydrolysis of the Chromogenic Substrates

Single-chain CHO-tPA, used for kinetic studies, was purified from CHO cells by affinity chromatography on red Sepharose and lysine-Sepharose (24-26). The enzyme concentration (in mg/ml) was determined from the absorbance at 280 nm with epsilon  = 1.81 ml/mg × cm. The molar concentration was calculated with a molecular weight of 60 kDa.

The substrates were dissolved in distilled water or dimethyl sulfoxide as necessary. Fifty microliters of the substrate solutions (six concentrations between 20 and 1 mM) were mixed with 0.2 ml HEPES buffer, pH 8.0 (0.1 mM, 0.154 NaCl, 0.1% human serum albumin), and prewarmed to 25 °C. The assays were started by the addition of 25 µl of tPA (34 µg/ml) prewarmed to 25 °C. Determination of the substrate cleavage was performed on microplates using a MR 5000 reader (Dynatech, Denkendorf, Germany). The amidolytic activity was calculated from the absorbance at 405 nm (epsilon  = 6.480 l/mmol × s). The velocity was plotted against the substrate concentration, and the data were fitted to the Michaelis-Menten equation with a nonlinear regression analysis program. The kcat was calculated from Vmax/enzyme concentration.

Crystallization and Structure Analysis

The des(Val4-Cys261)tPA variant (27) used for crystallization studies was expressed as inclusion bodies in Escherichia coli (strain K12 C600+, containing the helper plasmid pUBS520 with the DNA Y gene). Isolation and solubilization of the inclusion bodies, derivatization and refolding were transformed as described elsewhere (28-30). Briefly, E. coli cells were lysed by the addition of lysozyme and subsequent sonification, and the inclusion bodies isolated by centrifugation. Solubilization was achieved by incubation with M guanidium, HCl, and 0.1 M dithioerythriol at pH 8.5. After dialysis against 6 M guanidium HCl, the thiol groups were derivatized with glutathione by incubating the solubilized protein in 0.05 M Tris/HCl (pH 7.5), 6 M guanidine, and 0.1 M oxidized glutathione. After dialysis against 3 M guanidine, refolding was performed by stepwise addition of 300 ml of protein solution (mixed disulfide) to 10 liter of refolding buffer (0.7 M L-arginine/HCl, pH 8.6, 2 mM GSH, 1 mM EDTA) (31, 32). The refolded protein was purified by affinity chromatography on an Erythrina trypsin inhibitor-Sepharose column. The two-chain form was prepared from the refolded and purified tPA protein sample (more than 95% in single-chain form) by incubation with plasmin-Sepharose at room temperature under gentle agitation (27). The plasmin-Sepharose was removed by centrifugation.

Crystals were grown using the sitting-drop vapor-diffusion method. 3 µl of the protein (2.8 mg/ml in 10 mM acetate buffer, pH 5, and 0.8 mM inhibitor) and 1.5 µl of the reservoir solution (0.1 M acetate buffer, pH 5.0, 0.1 M (NH4)2SO4, 25% polyethylene glycol 4000) were mixed and equilibrated at 23 °C against 1 ml of reservoir solution. Plate shaped crystal bushels grew within 1 week. A crystal of 0.2×0.15×0.04 mm3 was mounted and belonged to the monoclinic space group C2 with cell constants a = 151.83 Å, b = 60.50 Å, c = 62.61 Å, alpha  = 90.0°, beta  = 110.5°, gamma  = 90.0°. Estimation of the solvent content indicated that the crystals contain 2 molecules per asymmetric unit with a fractional solvent content of 45% (Vm = 2.24 Å3/dalton). Data to 2.8 Å were collected on a multiwire Centronix area detector (Siemens) mounted on a Rigaku rotating anode x-ray generator and processed using XDS (33) (see Table I for crystallographic data). The structure was solved using Patterson search techniques with the program AMoRe (34) (integration radius 30 Å; resolution range, 15.0-3.5 Å). The search model consisted of the catalytic domain of tPA (4). The rotation function gave two peaks with 19.8 and 18.3 correlation, the next best angle triplet was less than 50% of the highest correlation value and not included in further translational searches. The translation function corresponding to the two highest peaks yielded correlation values of 25.5 (20.4 for next highest peak) and 27.5 (18.8 for next highest peak). After rigid body refinement with both molecules in the asymmetric unit, the correlation was 60.0 and the R-factor 35.7%.

Table I. Crystal data and refinement parameters for the bis-benzamidine-tPA complex


Spacegroup C2
Cell constants (Å)
  a 151.83
  b 60.50
  c 62.61
  beta 110.5°
Significant measurements 25355
Independent reflections 12586
Limiting resolution (Å) 2.8
Completeness 94.4% (infinity -2.8 Å)
  outermost shell 86.2% (3.1-2.8 Å)
Rmergea 8.8%
  outermost shell 25.2% (3.1-2.8 Å)
Number of atoms per asymmetric unit
  Non-hydrogen protein atoms 4066
  Non-hydrogen inhibitor atoms 56
  Solvent molecules 91
Resolution range (Å) 7.0-2.9
Reflections used for refinement 9017
Completeness (2sigma cutoff) 81.2%
  outermost shell 60.4% (3.0-2.9 Å)
Rvalueb 19.8%
Root mean square standard deviation
  Bond lengths (Å) 0.008
  Bond angles (degrees) 1.428
  RMSB (Å2)c 2.202

a Rmerge = Sigma (I - < I> )/Sigma I.
b Rvalue = Sigma (|Fobs- |Fcalc|)/Sigma |Fobs|.
c RMSB = root mean square deviation of the B-factor of bonded atoms.

Conventional crystallographic refinement (rigid body, positional and temperature factor) were carried out with XPLOR (35) using the parameters of Engh and Huber (36). In the latter stages restrained individual B-factor refinement was applied. Except for residues belonging to partially undefined surface loops, all atoms were restrained to obey strict noncrystallographic symmetry. Disordered residues (Tyr6, Gly36A-Gly37E, Asp110A-Ser110C and Arg243-Pro244 of both molecules in the asymmetric unit) and side-chains were excluded from phasing but included in the coordinate set (except the N-terminal and the C-terminal residues of the A-chain, Ser1E-Ser1B and Ser7-Phe11 of both molecules in the asymmetric unit). The molecular model of the inhibitor was constructed using InsightII (37) and model building was performed using Main (38). Averaging over the noncrystallographic symmetry elements using CCP4 (39) provided no significant improvement of the electron density. The program PROCHECK (40) indicates that 85% of the residues fall within the most favored region of the Ramachandran plot and that no residues are in the "disallowed regions." Crystallographic and refinement data are given in Table I. All figures showing the three dimensional structure of bisb-PA where prepared using either MAIN (38) or MOLSCRIPT (41).


RESULTS

Substrate Design

Kinetic analysis of the the hydrolysis of peptidic substrates (5, 6, 16, 42) and the three-dimensional structure of b-tPA (4) have shown that Arg and Gly are the preferred P1 and P2 residues for tPA substrates, respectively. One striking feature of b-tPA's active site region is Arg174, which suggests that acidic residues at P3 or P4 might provide additional stabilizing energies. To test this hypothesis, the P3 residue of p-nitroanilide substrates of the type R-Xaa-Gly-Arg-pNA was varied, and the kinetic constants for their cleavage by single chain tPA were determined and compared with known chromogenic substrates of this type (see Table II). In general, a clear preference is observed for D-configured amino acids at P3. Substitution of Xaa for the acidic residues D-Asp or D-Glu failed to produce the desired improvement in catalytic efficiency; indeed, these substrates were among the worst tested. On the other hand, the P3 residue L-Asp exhibited the lowest Km value of all L-configured P3 residues. In this case, however, catalysis by tPA was impaired, suggesting that L-Asp-Gly-Arg-pNA does not interact with the active site region of tPA in a suitable conformation for cleavage.

Table II. Kinetic data for the hydrolysis of p-nitroanilide substrates of the type R-Xaa-Gly-Arg-pNA, where R is an N-terminal protection group, by full-length single-chain tPA


R P3 Km kcat kcat/Km

mM l/s l/mmol × s
Boc Gly 3.50 1.29 0.37
Cbo Val 1.36 1.69 1.24
Cmo D-Val 0.80a 5.77 7.22
Cmo D-Nle 0.99 6.82 6.89
Cmo D-Leu 0.62 6.92 11.2
H Phe 2.00 2.79 1.40
H D-Phe 0.70 7.46 10.7
Cbo D-Phe 0.30 8.96 19.9
MeSO2 D-Phe 0.10 7.36 73.6
H CHG 5.32 1.36 0.26
H D-CHG 1.53 6.74 4.40
Cmo D-CHG 0.42 5.02 11.9
H D-HHT 0.83 10.6 12.8
Cmo D-HHT 0.39 8.31 21.3
Cmo D-CHA 0.34 7.76 22.8
H Asp 0.86 0.28 0.32
Cbo Asp 0.72 0.35 0.49
Cbo Asp(O-t-BU) 2.30 1.84 0.80
H D-Asp 1.70 0.43 0.25
Cbo D-Asp(O-t-BU) 1.80 4.98 2.77
H Glu 3.20 0.13 0.041
Cbo Glu 2.20 0.75 0.34
Cbo Glu(O-t-BU) 0.90 1.29 1.43
H D-Glu 2.90 3.98 1.37
Cbo D-Glu 2.60 7.46 2.87
Cbo D-Glu(O-t-BU) 1.20 4.48 3.73
pGlu 3.15 4.04 1.28
H D-Lys 2.10 6.32 3.01
Cbo D-Arg 0.52 8.53 16.4
Ac Acetyl -CO-CH3
Boc tert-Butoxycarbonyl -CO-O-C(CH3)
Cbo Benzyloxycarbonyl -CO-O-CH2-C6H5
CHA Cyclohexylalanyl -NH-CH-(CH2-C6H11)-CO-
CHG Cyclohexylglycyl -NH-CH-(C6H11)-CO-
HHT Hexahydrotyrosyl -NH-CH-(CH2-C6H10-OH)-CO-
MeSO2 Methylsulfonyl -SO2-CH3
pNA p-Nitroanilide -NH-C6H4-NO2                 
||
pGlu Pyroglutamic acid CH2-CH2-CO-NH-CH-CO-
t-BU tert-Butyl -C(CH3)3

a Km values <1, kcat values >5.0, and kcat/Km values >10 are emphasized in bold.

The substrates with the highest specificity constant (kcat/Km) were of the form R-D-Xaa-Gly-Arg-pNA, where Xaa is an hydrophobic (phenylalanine or hexahydrotyrosine) or basic (Arg or Lys) residue. Among these substrates, a further gradation in catalytic efficiency was observed according to the protection group of the amino terminus (R-). In each case, a free amino terminus was suboptimal. The kinetic constants of all measured chromogenic substrates of the type R-Xaa-Gly-Arg-pNA are given in Table II.

Structure of tPA in Complex with Bis-benzamidine

Many specific features of tPA involve interaction with surface elements of the catalytic domain, in particular the surface loops (see Fig. 1). Their role in substrate recognition, cofactor binding and inhibition under physiological conditions has been discussed comprehensively elsewhere (4). In that study, some of the surface loops were not defined by electron density and therefore termed flexible. The crystal packing of the structure described here is quite different, leading to different environments of the molecular surface and allowing possible visualization of labile parts of the structure. Although the two molecules in the asymmetric unit are related by noncrystallographic symmetry, i.e. they are not identical, no significant differences are observed between them. The overall structure of the bis-benzamidine liganded tPA (bisb-tPA) (Figs. 1 and 2) shows almost complete structural identity with the benzamidine liganded (b-tPA) structure (4).

Most regions that were undefined in b-tPA remain so in bisb-tPA. Thus, residues of the A-chain amino-terminal to Thr1A and carboxyl-terminal to Tyr6 are disordered, as are the six residues Lys36, His37, Arg37A, Arg37B, Ser37C, Pro37D of the 37 insertion loop. The latter has been shown to be of fundamental importance for the interaction with PAI-1 (8-11) and for fibrin specificity and stimulation (43, 44). In contrast to b-tPA, the three residues Asp110A, Ser110B, and Ser110C of the 110 insertion loop are undefined in bisb-tPA; the equivalent loop in uPA (45) also exhibits high mobility. The last two C-terminal residues Arg243 and Pro244, which are in structural proximity to the 110 loop, are also undefined in bisb-tPA.

The "autolysis" loop Tyr141 to Tyr151 to the south of the active site cleft is well defined by electron density, with the same overall conformation as observed in b-tPA. However, the side chain of Glu145 forms a salt bridge with Nzeta Lys17 in bisb-tPA instead of the H-bond to N Leu187 in b-tPA. This is probably a result of different crystal contacts.

The 186-loop to the south of the active site cleft consists of 8 additional residues compared with alpha -chymotrypsinogen. In contrast to b-tPA, all residues are defined by electron density (see Fig. 3). The extended loop projects out from the molecular surface into the solvent, and is stabilized by crystal contacts. The cluster of three residues Arg186A, Gln186F, His188, which point toward the main body of the molecule, might also stabilize this arrangement. It is not clear which is the main determinant for the stabilization of this loop in this structure. It is possible that the different pH conditions used for the crystallization of b-tPA (pH 7.5) and bisb-tPA (pH 5) could have lead to this conformational difference. Thus, the deprotonation of His188 in b-tPA might disfavor formation of the stabilizing cluster Arg186A, Gln186F, His188.


Fig. 3. Experimental density (2Fo - Fc) for the 186-loop of tPA, showing the cluster of Arg186A, Gln186C, His188. Gly19 of the activation peptide is adjacent to this loop (lower right), as is the strand His159-Arg161 (bottom). Orientation as in Fig. 1.
[View Larger Version of this Image (51K GIF file)]

The exact conformation of this loop might not be important, as position 186G is deleted in mouse and rat tPA. Its function is not clear, but it may be involved in stimulation of the plasminogen activator by its cofactor fibrin (46). It is bordered on its eastern side by the N-terminal part of the B-chain, with the side chain of His188 approaching Gly18 and Gly19. In this conformation, the 186-loop is close to the entrance frame of the S1 pocket (Ser214-Gln221A) and the following residues Lys222, Asp223, Val224, which form a turn.

Interaction of Bis-benzamidine with the Active Site of tPA

The inhibitor is well defined by electron density (see Fig. 4). The binding of the bis-benzamidine to tPA is determined by two major interaction sites. The first amidino group, the "proximal" amidino group, binds in the specificity pocket while the second amidino group, the "distal" amidino group, fits in an hydrophobic groove, resulting in an extended binding of the inhibitor. The inhibitor exhibits a trans/trans conformation (see Fig. 2). The proximal amidino group is sandwiched in the S1 pocket between main-chain segments Trp215-Gly216 and Cys191-Gln192 in the same way as amidine in b-tPA. Its amidino function forms a symmetric two O-two N salt bridge with Asp189 and a N-O hydrogen bond with the carbonyl oxygen of Gly219. The proximal amidino group is almost coplanar and coincident with the superimposed amidine from b-tPA, yet its amidino group appears to be more out of the aromatic ring plane than in the latter structure.


Fig. 4. Active site region of human tPA in complex with the bis-benzamidine inhibitor. The proximal benzamidine occupies the specificity pocket, and the distal benzamidine reaches toward an electrophilic pocket formed by the side chain carbonyl atoms of Asp97, Thr98, and Arg174. The benzyl ring of the inhibitor is perpendicular to the side chain of Tyr99; the side chain of Arg174 is displaced by the amidino group. A, experimental density (2Fo - Fc) of the active site region in complex with the bis-benzamidine inhibitor. B, active site region of human tPA in complex with bis-benzamidine. Arg174 from b-tPA is superimposed (thin lines).
[View Larger Version of this Image (84K GIF file)]

The heptanone ring system, which acts as a spacer between the two benzamidino groups, partially occupies the S2 pocket, resulting in a close contact between one carbon atom of the cycloheptanone ring and the OH-group of the Tyr99 (3 Å distance). The restriction of the S2 pocket by the side chain of Tyr99 (see above) does not allow more bulky spacers in bis-benzamidine inhibitors, which correlates with previous published results (42). The rigid nature of the cycloheptanone ring results in a significant distance (3.8 Å) between the carbonyl oxygen of the inhibitor and the backbone of the segment Ile213-Cys220. This distance would be smaller for a cyclopentanone and cyclohexanone spacer, explaining their weaker inhibition of tPA. Gly216 does not contribute to binding by hydrogen bond formation, as seen in some synthetic inhibitor complexes with the related serine proteinase thrombin (47).

The distal benzamidino group penetrates into the hydrophobic groove, displacing the side chain of Arg174. The aliphatic portion of the Arg side chain now makes up the western border of the S4 pocket, making it more hydrophobic and flanking the distal benzamidine. The aromatic benzene ring is parallel to the indole moiety of the Trp215 and edge to face to the phenolic side chain of Tyr99, which may confer additional weak binding energy. Indeed, a complex of the same inhibitor in trypsin,3 indicates a similar overall geometry of the inhibitor, but with the distal amidino group rotated 50° around its molecular axis; the aliphatic side chain of Leu99 in trypsin is unable to achieve the favorable aromatic stacking arrangements seen in tPA. The basic amidino group occupies the same position as the guanidino group of Arg174 in b-tPA (see Fig. 4), approaching an electronegative cavity on the surface of tPA formed by the carbonyl groups of Asp97 and Thr98.


DISCUSSION

A series of p-nitroanilide substrates were synthesized on the premise that the guanidino group of Arg174 partially occupies the S4 site of tPA (4). Substrates containing an acidic D-amino acid at P3, which could form a stable salt bridge to this residue, failed to give the desired improvement in catalytic efficiency. On the other hand, substrates containing an L-Asp residue at P3 exhibited the lowest Km values of all the L-P3 residues, comparable to the Km values of the most specific substrates. In this case, however, cleavage of the substrate is compromised, indicating that binding of P3 L-Asp produces an incorrect peptide conformation at the cleavage site. These results suggest that the conformation of Arg174 as seen in benzamidine-liganded tPA is not fixed.

The 2,7-bis-(4-amidinobenzylidene)-cycloheptan-1-one tPA structure shows that the side chain of this residue can indeed adopt multiple conformations. While the proximal amidino group occupies the S1 pocket as expected, the distal aromatic ring displaces the side chain of Arg174 from its position in b-tPA (see Fig. 4). The aliphatic side chain now rests against that of Leu217, resulting in an hydrophobic western edge to the S3/S4 pocket. The distal amidino group of the inhibitor occupies an electronegative cavity formed by the carbonyl groups of Asp97, Thr98, and Thr175, formerly the location of the guanidino group of Arg174. The cycloheptanone ring system, which makes negligible contacts with the enzyme and acts as a spacer for the two benzamidine moieties, presumably provides an entropic advantage by restricting the conformational degrees of freedom in the unbound state.

This new structural data allows rationalization of the substrate selectivity shown in Table II. The S4 pocket is large and hydrophobic, bounded by Trp215 at the base, Tyr99 to the east, and the aliphatic side chain of Arg174 to the west. This explains the observed preference of tPA for the hydrophobic residues D-Phe and D-hexahydrotyrosine at P3. The presence of an electronegative cavity to the north, formed by the carbonyl groups of Asp97, Thr98, and Thr175, provides a receptor for basic P3 residues such as D-Arg. This electrostatic interaction is presumably weak, as the shorter and less basic D-Lys exhibits a high Km. The protection groups connected to the amino terminus greatly influence the binding constants.

The peptide moieties of Cbo-D-Arg-Gly-Arg-pNA and MeSO2-D-Phe-Gly-Arg-pNA were modeled into the active site of bisb-tPA assuming a "canonical conformation" (48) of the peptides. Although small substrates are in general more flexible than protein substrates and inhibitors, a survey of structures of serine proteinase containing short peptide chloromethyl ketone inhibitors (such as D-Phe-Pro-Arg-chloromethylketone in thrombin) (18) suggests that also small substrates in the Michaelis complex adopt canonical conformation. For the substrates discussed here, we might expect more flexibility associated with the P2 glycyl residues; recent structures of uPA in complex with Glu-Gly-Arg-chloromethylketone (45), single chain tPA in complex with dansyl-Glu-Gly-Arg-chloromethylketone (49), and vampire bat plasminogen activator in complex with Glu-Gly-Arg-chloromethylketone (50) reveal, however, that the P2 Gly residue adopts canonical phi -psi angles in the bound state.

MeSO2-D-Phe-Gly-Arg-pNA

Assuming a short antiparallel beta -sheet between the substrate Gly and Phe residues and the enzyme Trp215 and Gly216, the side chain of P3 D-Phe fits perfectly into the hydrophobic entrance of the S4 pocket (see Fig. 5). The Phe aromatic ring stacks over the indole system of Trp215, with concomitant edge to plane stacking to Tyr99. One oxygen atom of the methylsulfonyl group can form a hydrogen bond with the amide nitrogen of Gly219, while the methyl group could rest against the side chain of Leu217, providing an explanation for the sevenfold decrease in Km for MeSO2-D-Phe-Gly-Arg-pNA compared with H-D-Phe-Gly-Arg-pNA (see Table II).


Fig. 5. Models of the Michaelis complex for two R-Xaa-Gly-Arg-pNA-substrates for tPA with high affinity (see Table II). The peptidic moieties were modeled into the active site of bisb-tPA according to the binding of D-Phe-Pro-Arg chloromethylketone in thrombin (18, 19). The positions of the D-P3 side-chains were modeled based on the bisb-tPA structure. The phenylalanyl side chain of MeSO2-D-Phe-Gly-Arg-pNA (thick black connections) is stabilized by contacts to the aromatic moieties of Tyr99, Trp215, and the aliphatic portion of Arg174. An oxygen atom of the MeSO2 group could make an hydrogen bond to the amide nitrogen of Gly219. The P3 guanidino group of Cbo-D-Arg-Gly-Arg-pNA (gray connections) occupies the electrophilic cavity determined for bisb-tPA. Additional stabilization is afforded by the Cbo group, which covers the hydrophobic surface provided by residues Leu217, Val224. and Leu172.
[View Larger Version of this Image (41K GIF file)]

Cbo-D-Arg-Gly-Arg-pNA

Assuming the same antiparallel beta -sheet as above, the side chain of P3 D-Arg extends deeply into the S4 pocket to reach the electronegative cavity (see Fig. 5). A possible stacking of the aliphatic portions of P3 D-Arg and Arg174 could make a further contribution to the binding energy. The Cbo protecting group points in an extended conformation toward an hydrophobic patch formed by the residues Leu172, Leu217, Gln221A, and Val224, allowing partial burial of this hydrophobic surface by the benzoxy group.

These models suggest possible modifications for improved substrates. Superposition of the two substrates indicate coincidence of P3 D-Phe Czeta (or Cdelta 1) with P3 D-Arg Nepsilon . A P3 D-amidinophenyl group could make dual use of the hydrophobic environment at the entrance at the pocket and the electrophilic environment at the base of the pocket for stabilizing interactions. The large influence of a protection group at the N terminus on catalytic efficiency suggests that an additional hydrophobic residue at position P4 could allow additional stabilization through burial of tPA's hydrophobic surface to the west of the active site.

Our results concerning the specificity sites of tPA are summarized in Fig. 6. Outside the primary specificity pocket, the major determinants for interaction with substrates and inhibitors would appear to be binding in the hydrophobic S4 pocket, together with weak electrostatic interactions at the electronegative cavity. Such an extended binding mode has also been observed for factor Xa specific inhibitors (51, 52). Indeed, the 99-loop forming the north-eastern edge of the S4 pocket is almost identical in factor Xa (Glu97-Thr98-Tyr99). This strongly suggests that both tPA and factor Xa could be inhibited by similar ligands. This raises the possibility that use of synthetic inhibitors of factor Xa as antithrombotics could have the undesired side effect of interfering with fibrinolysis. Extreme care must therefore be taken to manipulate the selectivity of such inhibitors; in this case, this could be achieved by exploiting differences at the western edge of the S4/S3 site.


Fig. 6. Schematic diagram of the active site specificity determinants of human tPA. The relative sizes of the S1, S2, S3, and S4 pockets emphasize their volume ratio in the three-dimensional structure. Asp189 at the bottom of the S1 specificity-pocket results in a negatively charged environment. The S4 pocket, delimited by the side chains of Tyr99 at its eastern side, Trp215 at the bottom and the aliphatic part of Arg174 (represented in the orientation found in bisb-tPA) at its western border, renders the entrance frame hydrophobic. At the north of the S4 pocket, the carbonyl groups of Asp97, Asp98, and Thr175, form a negative cavity.
[View Larger Version of this Image (69K GIF file)]


FOOTNOTES

*   This work was supported by Sonderforschungsbereich (SFB) 207, SFB 469, and European Union Project CHRX-CT94-0535.The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
§   To whom correspondence should be addressed. Tel.: 49-(0)89-8578 2676; Fax: 49-(0)89-8578 3516; E-mail: renatus{at}biochem.mpg.de or bode{at}biochem.mpg.de.
**   Present address: Institute for Pharmaceutical Chemistry, Philipps-University Marburg, Marbacher Weg 6, D-35032 Marburg, Germany.
1   The abbreviations used are: tPA, tissue type plasminogen activator; pNA para-nitroanilide; b-tPA, catalytic domain of activated tissue type plasminogen activator in complex with benzamidine inhibitor; bisb-tPA, catalytic domain of activated tissue type plasminogen activator in complex with the bis-benzamidine inhibitor 2,7-bis-(4-amidinobenzylidene)-cycloheptan-1-one; HHT, hexahydrotyrosyl; CHG, cyclohexylglycyl; CHA, cyclohexylalanyl; t-BU, tert-butyl; Cbo, benzyloxycarbonyl; Boc, tert-butoxycarbonyl; Cmo, methyloxycarbonyl; MeSO2, methylsulfonyl; ETI, Erythrina trypsin inhibitor; CHO, Chinese hamster ovary; dansyl, 5-dimethylaminonaphthalene-1-sulfonyl.
2   Residues in this publication are numbered using the chymotrypsinogen numbering, based on the topological equivalence of the proteolytic domain of human tPA with alpha -chymotrypsinogen (4).
3   M. Renatus, W. Bode, R. Huber, J. Stürzebecher, D. Prasa, and M. T. Stubbs, manuscript in preparation.

ACKNOWLEDGEMENTS

We are grateful to Peter Wikström (Pentapharm Ltd, Basel, Switzerland) for providing the inhibitor for crystallization studies. Assistance in the synthesis of substrates by Lars Svendsen (Pentapharm Ltd, Basel, Switzerland) is gratefully acknowledged.


REFERENCES

  1. Madison, E. L. (1994) Fibrinolysis 8, 221-236
  2. Rijken, D. C., Hoylaerts, M., and Collen, D. (1982) J. Biol. Chem. 257, 2920-2925 [Free Full Text]
  3. Rånby, M. (1982) Biochim. Biophys. Acta 704, 461-469 [Medline] [Order article via Infotrieve]
  4. Lamba, D., Bauer, M., Huber, R., Fischer, S., Rudolph, R., Kohnert, U., and Bode, W. (1996) J. Mol. Biol. 258, 117-135 [CrossRef][Medline] [Order article via Infotrieve]
  5. Madison, E. L., Coombs, G. S., and Corey, D. R. (1995) J. Biol. Chem. 270, 7558-7562 [Abstract/Free Full Text]
  6. Coombs, G. S., Dang, A. T., Madison, E. L., and Corey, D. R. (1996) J. Biol. Chem. 271, 4461-4467 [Abstract/Free Full Text]
  7. Bachmann, F. (1994) Haemostasis and Thrombosis, Vol. 1, pp. 575-613, Churchill Livingstone, New York
  8. Madison, E. L., Goldsmith, E. J., Gerard, R. D., Gething, M.-J. H., and Sambrook, J. F. (1989) Nature 339, 721-724 [CrossRef][Medline] [Order article via Infotrieve]
  9. Madison, E. L., Goldsmith, E. J., Gerard, R. D., Gething, M.-J. H., Sambrook, J. F., and Bassel-Duby, R. S. (1990) Proc. Natl. Acad. Sci. U. S. A. 87, 3530-3533 [Abstract]
  10. Madison, E. L., Goldsmith, E. J., Gething, M.-J. H., Sambrook, J. F., and Gerard, R. D. (1990) J. Biol. Chem. 265, 21423-21426 [Abstract/Free Full Text]
  11. Camani, C., Bachmann, F., and Kruithof, E. K. O. (1994) J. Biol. Chem. 269, 5770-5775 [Abstract/Free Full Text]
  12. Heussen, C., Joubert, F., and Dowdle, E. B. (1984) J. Biol. Chem. 259, 11635-11638 [Abstract/Free Full Text]
  13. Joubert, F. J., and Dowdle, E. B. D. (1987) Thromb. Haemostasis 57, 356-360 [Medline] [Order article via Infotrieve]
  14. Teixeira, A. V., Dowdle, E. B. D., and Botes, D. P. (1994) Biochim. Biophys. Acta 1217, 23-28 [Medline] [Order article via Infotrieve]
  15. Markwardt, F., and Stürzebecher, J. (1989) Design of Enzyme Inhibitors as Drugs, pp. 619-649, Oxford University Press, Oxford, New York
  16. Ding, L., Coombs, G. S., Strandberg, L., Navre, M., Corey, D. R., and Madison, E. L. (1995) Proc. Natl. Acad. Sci. U. S. A. 92, 7627-7631 [Abstract]
  17. Padmanabhan, K., Padmanabhan, K. P., Tulinsky, A., Park, C. H., Bode, W., Huber, R., Blankenship, D. T., Cardin, A. D., and Kisiel, W. (1993) J. Mol. Biol. 232, 947-966 [CrossRef][Medline] [Order article via Infotrieve]
  18. Bode, W., Mayer, I., Baumann, U., Huber, R., Stone, S. R., and Hofsteenge, J. (1989) EMBO J. 8, 3467-3475 [Abstract]
  19. Bode, W., Turk, D., and Karshikov, A. (1992) Protein Sci. 1, 426-471 [Abstract/Free Full Text]
  20. Stürzebecher, J., Markwardt, F., Wagner, G., and Walsmann, P. (1982) Pharmazie 37, 281-283 [Medline] [Order article via Infotrieve]
  21. Stürzebecher, J., Markwardt, F., Wagner, G., and Walsmann, P. (1978) Acta Biol. Med. Ger. 35, 1665-1676
  22. Walsmann, P., Horn, H., Markwardt, F., Richter, P., Stürzebecher, J., Vieweg, H., and Wagner, G. (1976) Acta Biol. Med. Ger. 35, K1-K8 [Medline] [Order article via Infotrieve]
  23. Svendsen, L. G. (1983) European Patent 0 034 122 B1
  24. Macartney, H. W., and Menhardt, J. (1989) European Patent Application 502503
  25. van Zonneveld, A.-J., Veerman, H., and Pannekoek, H. (1986) J. Biol. Chem. 261, 14214-14218 [Abstract/Free Full Text]
  26. Kohnert, U., Rudolph, R., Verheijen, J. H., Weening-Verhoeff, E. J. D., Stern, A., Opitz, U., Martin, U., Lill, H., Prinz, H., Lechner, M., Kresse, G.-B., Buckel, P., and Fischer, S. (1992) Protein Eng. 5, 93-100 [Abstract]
  27. Pannekoek, H., Lijnen, H. R., and Loskutoff, D. J. (1990) Thromb. Haemostasis 64, 600-603 [Medline] [Order article via Infotrieve]
  28. Kohnert, U., Hellerbrand, K., Martin, U., Stern, A., Popp, F., and Fischer, S. (1996) Fibrinolysis 10, 93-102
  29. Rudolph, R. (1990) in Modern Methods in Protein and Nucleic Acid Research (Tschesche, H., ed), pp. 149-172, Walter de Gruyter, Berlin
  30. Stern, A., Kohnert, U., Rudolph, R., Fischer, S., and Martin, U. (1990) European Patent Application 382174
  31. Rudolph, R., and Fischer, S. (1990) United States Patent 4 933 434
  32. Rudolph, R., Fischer, S., and Mattes, R. (1987) European Patent Application 219874
  33. Kabsch, W. (1988) J. Appl. Crystallogr. 21, 916-924 [CrossRef]
  34. Navaza, J. (1994) Acta Crystallogr. A 50,, 157-163
  35. Brünger, A. (1992) X-PLOR (Version 3.1). A System for X-ray Crystallography and NMR, Yale University Press, New Haven, CT
  36. Engh, R. A., and Huber, R. (1991) Acta Crystallogr. A 47,, 392-400 [CrossRef]
  37. Biosym Technologies San Diego (1993) InsightII, Version 2.2.0
  38. ät MünchenTurk, D. (1992) Weiterentwicklung eines Programmes für Molekülgraphik und Elektronendichte-Manipulation und seine Anwendung auf verschiedene Protein-Strukturaufklärungen. Ph.D. thesis, Technische Universität München
  39. Collaborative Computational Project, Number 4 (1994) Acta Crystallogr. Sec. D 50, 760-763 [CrossRef][Medline] [Order article via Infotrieve]
  40. Laskowski, R. A., MacArthur, M. W., Moss, S. D., and Thorton, J. M. (1993) J. Appl. Crystallogr. 26, 283-291 [CrossRef]
  41. Kraulis, P. J. (1991) J. Appl. Crystallogr. 24, 946-950 [CrossRef]
  42. Stürzebecher, J., Neumann, U., Kohnert, U., Kresse, G.-B., and Fischer, S. (1992) Protein Sci. 1, 1007-1013 [Abstract/Free Full Text]
  43. Bennett, W. F., Paoni, N. F., Keyt, B. A., Botstein, D., Jones, A. J. S., Presta, L., Wurm, F. M., and Zoller, M. J. (1991) J. Biol. Chem. 266, 5191-5201 [Abstract/Free Full Text]
  44. Paoni, N. F., Keyt, B. A., Refino, C. J., Chow, A. M., Nguyen, H. V., Bereleau, L. T., Badillo, J., Peña, L. C., Bardy, K., Wurm, F. M., Ogez, J., and Bennett, W. F. (1993) Thromb. Haemostasis 70, 307-312 [Medline] [Order article via Infotrieve]
  45. Spraggon, G., Phillips, C., Novak, U. K., Ponting, C. P., Saunders, D., Dobson, C. M., Stuart, D. I., and Jones, E. Y. (1995) Structure 3, 681-691 [Medline] [Order article via Infotrieve]
  46. Paoni, N. F., Chow, A. M., Peña, L. C., Keyt, B. A., Zoller, M. J., and Bennett, W. F. (1993) Protein Eng. 6, 529-534 [Abstract]
  47. Brandstetter, H., Turk, D., Hoeffken, H. W., Grosse, D., Stürzebecher, J., Martin, P. D., Edwards, B. F. P., and Bode, W. (1992) J. Mol. Biol. 226, 1085-1099 [Medline] [Order article via Infotrieve]
  48. Bode, W., and Huber, R. (1994) Fibrinolysis 8, 161-171
  49. Renatus, M., Engh, R. A., Stubbs, M. T., Huber, R., Fischer, S., Kohnert, U., and Bode, W. (1997) EMBO J., in press
  50. Renatus, M., Stubbs, M. T., Huber, R., Bringmann, P., Donner, P., Schleuning, W.-D., and Bode, W. (1997) Biochemistry, in press
  51. Stubbs, M. T., Huber, R., and Bode, W. (1995) FEBS Lett. 375, 103-107 [CrossRef][Medline] [Order article via Infotrieve]
  52. Brandstetter, H., Kühne, A., Bode, W., Huber, R., von der Saal, W., Wirthensohn, K., and Engh, R. A. (1996) J. Biol. Chem. 271, 29988-29992 [Abstract/Free Full Text]

©1997 by The American Society for Biochemistry and Molecular Biology, Inc.