Approximately 2% of mammalian genes encode proteases. Comparative genomics reveals that those involved in immunity and reproduction show the most interspecies diversity and evidence of positive selection during evolution. This is particularly true of granzymes, the cytotoxic proteases of natural killer cells and CD8+ T cells. There are 5 granzyme genes in humans and 10 in mice, and it is suggested that granzymes evolve to meet species-specific immune challenge through gene duplication and more subtle alterations to substrate specificity. We show that mouse and human granzyme B have distinct structural and functional characteristics. Specifically, mouse granzyme B is 30 times less cytotoxic than human granzyme B and does not require Bid for killing but regains cytotoxicity on engineering of its active site cleft. We also show that mouse granzyme A is considerably more cytotoxic than human granzyme A. These results demonstrate that even “orthologous” granzymes have species-specific functions, having evolved in distinct environments that pose different challenges.
Granzymes are cytotoxic serine proteases used by cytotoxic lymphocytes (CLs) to destroy virus-infected and malignant cells. They are delivered into the cytoplasm of the target cell by the pore-forming protein perforin. Once inside the target cell, granzymes cleave specific proteins and trigger apoptosis (for review see Russell and Ley, 2002).
Humans have 5 granzymes (A, B, H, K, and M), whereas mice have 10 granzymes, having duplicated the ancestral granzyme B/H gene. In the absence of reported human granzyme deficiency, understanding of granzymes has been built by experimental approaches involving knockout mice or the use of purified granzymes on cells or cell extracts. For example, loss of granzyme A (GrA) or granzyme B (GrB) in mice is not catastrophic, but loss of both results in a severe immunological defect. Numerous studies using purified human and rodent granzymes have identified substrates, regulators, and trafficking pathways, whereas granzyme- deficient mice are widely used to probe the contribution and importance of granzymes to immunity (Russell and Ley, 2002).
GrA and GrB are the major and most-studied granzymes. GrA cleaves substrates after basic residues and induces caspase-independent apoptosis (for review see Lieberman and Fan, 2003). GrB cleaves after acidic residues, primarily Asp, and induces apoptosis by directly or indirectly activating caspases, and by cleaving other proapoptotic proteins such as Bid or ICAD/DFF45 (for review see Trapani and Sutton, 2003). Studies using CL from knockout mice confirm that GrA and GrB trigger distinct apoptotic pathways and are thus complementary (Pardo et al., 2004).
The high sequence homology and conserved primary cleavage specificity of human and mouse granzymes has led to widespread and interchangeable use of human and mouse enzymes in experimental systems, usually without side-by-side comparisons being made. Because GrA is a unique gene in both species, it is reasonable to assume an orthologous relationship between mouse GrA (mGrA) and human GrA (hGrA), but recent structural data show that they have distinct extended substrate specificities (Bell et al., 2003), raising the possibility that they possess distinct functions. The situation with GrB is even more complicated. hGrB has multiple paralogues in the mouse and rat, suggesting that the roles fulfilled by a single enzyme in humans may be split among several in rodents (Pham et al., 1996; Revell et al., 2005).
Without detailed comparative knowledge of the characteristics of human and mouse granzymes, data from mice cannot be extrapolated with confidence to humans. Here, we show that the cytotoxicity, substrate preferences, and inhibitor interactions of hGrB and mGrB are substantially different. We also extend the study of Bell et al. (2003) to show that hGrA and mGrA have markedly different cytotoxic potential. Such differences probably reflect the two species' response to differing immune evolutionary pressure.
hGrB is more cytotoxic than mGrB
Granzyme cytotoxicity can be studied in vitro by adding purified perforin and granzyme to cultured cells and monitoring cell survival (Fig. S1, available at http://www.jcb.org/cgi/content/full/jcb.200606073/DC1). Keeping the concentration of perforin constant, granzyme cytotoxicity is quantitated using a dose–response approach, in which the amount of granzyme required to kill 50% of the target cells (EC50) is measured (Sun et al., 2004). Initial experiments using this system with either human or mouse cells showed that hGrB is substantially more cytotoxic than mGrB.
This apparent difference in cytotoxicity between mGrB and hGrB was confirmed using several different targets, as cell type–dependent differences in susceptibility to granzymes have been reported (Pardo et al., 2002). In one human cell line and three mouse cell lines examined in detail, hGrB was consistently more cytotoxic than mGrB (Fig. 1). With the exception of YAC-1, cells were 10–60-fold more sensitive to hGrB (Fig. 1). Killing of mouse YAC-1 cells required similarly high levels of hGrB or mGrB, suggesting that a key substrate is low or absent in these cells. mGrB was also 30-fold less effective in killing mouse lymphoma cells (Fig. 2 C).
To test whether poorer uptake of mGrB explains its inefficient killing, we compared accumulation of fluoresceinated mGrB and hGrB in Jurkat cells (Bird et al., 2005). No impairment in the rate or extent of accumulation of mGrB versus hGrB was noted in these cells (Fig. S2, available at http://www.jcb.org/cgi/content/full/jcb.200606073/DC1) or in the mouse P815 and lymphoma cell lines used in this study. The differences in cytotoxicity between mGrB and hGrB also remained when the bacterial porin, streptolysin O (SLO) was used instead of perforin (Fig. 1 C). SLO admits proteins directly to the cytoplasm via transient plasma membrane pores (Walev et al., 2001) and provides a perforin-independent mechanism for granzyme uptake, as shown by its ability to deliver mutant granzymes that cannot be delivered by perforin (Bird et al., 2005). The fact that mGrB killing is as inefficient in SLO as it is in perforin suggests that there is no restriction in perforin-mediated delivery of mGrB. Thus, differences in uptake or delivery of hGrB and mGrB do not explain their markedly different cytotoxic potential.
mGrB-mediated cytotoxicity does not require Bid
To investigate the basis for the differences in cytotoxicity between mGrB and hGrB, we compared their activities on peptide substrates and on the natural substrates procaspase 3 and Bid (Metkar et al., 2003; Waterhouse et al., 2005). Combinatorial peptide substrate analysis has identified the peptide IleGluXxxAsp↓ as an optimal recognition and cleavage sequence for hGrB and rat GrB (Thornberry et al., 1997; Harris et al., 1998). This appears in human and mouse Bid, and both are cleaved at this site by hGrB (Pinkoski et al., 2001). Although mGrB cleaves Bid (Thomas et al., 2001), the cleavage site has not been determined. A similar P4-P1 cleavage sequence is present on procaspase 3, but the downstream (P′) sequence differs from Bid. The nomenclature of cleavage site positions of a protease substrate follows Schechter and Berger (1967). The protease cleaves the peptide bond between P1 and P1′. Adjacent residues in the N-terminal direction are numbered P2, P3, P4, etc. On the carboxyl side of the cleavage site, residues are numbered P2′, P3′, P4′, etc. The P1 residue is accommodated in the S1 pocket of the catalytic cleft, the P2 residue is accommodated in the S2 pocket, and so forth.
mGrB and hGrB both cleaved the small peptide substrate Boc-AlaAlaAsp-thio benzyl ester (AAD-sbzl) with similar efficiency (Fig. 2 A, right), consistent with their known Aspase activities, indicating that their apparent functional distinctions are not due to differences in ability to hydrolyse a P1(Asp)-P1′ bond. As expected, hGrB also efficiently cleaved IleGluThrAsp-pNA (IETD-pNA; Fig. 2 A, left). In contrast, mGrB hydrolyzed IETD-pNA 30-fold less efficiently than hGrB (Fig. 2 A, left), suggesting that it cannot comfortably accommodate this peptide in its substrate binding pocket.
The markedly inferior ability of mGrB to cleave IETD-pNA also suggested that it would not cleave full-length Bid or procaspase 3 efficiently. To test this, we produced 35S-labeled mouse Bid and mouse procaspase 3 and incubated the protein with hGrB or mGrB. Under these conditions, complete cleavage of Bid was achieved in 25 nM hGrB, whereas 2 μM mGrB failed to completely cleave Bid (Fig. 2 B, right). Thus, mGrB cleaves Bid 80–100-fold less efficiently than hGrB. In contrast, the two granzymes cleaved procaspase 3 with similar efficiency (Fig. 2 B, left), suggesting that other factors such as P′ residues also influence cleavage (Table I and Fig. 3).
Previous studies using hGrB have indicated that Bid cleavage is a key event in GrB-mediated apoptosis (Goping et al., 2003; Sutton et al., 2003; Waterhouse et al., 2005). GrB initiates procaspase 3 activation by cleaving it once, but full activation requires a second cleavage, autocatalytically or via another caspase. This is illustrated in Fig. 2 B and Fig. S3 (available at http://www.jcb.org/cgi/content/full/jcb.200606073/DC1) by the inclusion of the caspase inhibitor z-VAD in the reactions. In vivo, this second cleavage step is prevented by inhibitor of apoptosis protein (IAP) bound to caspase 3 and will not occur unless IAP is displaced by Smac/Diablo released from mitochondria by Bid (Goping et al., 2003; Sutton et al., 2003). Thus, in the absence of Bid, hGrB cannot drive full caspase 3 activation.
To examine the importance of Bid in mGrB-mediated apoptosis, we compared the ability of hGrB and mGrB to kill mouse cells lacking Bid. As reported previously (Waterhouse et al., 2005), hGrB efficiently killed mouse lymphoma cells containing Bid (EC50 2.5 U/ml), but was >1,000-fold less cytotoxic toward lymphoma cells lacking Bid (Fig. 2 C, left). In contrast, mGrB killed the Bid+/+ and Bid−/− cells with similar efficiency, indicating that Bid cleavage is not essential for mGrB cytotoxicity (Fig. 2 C, right). These observations cannot be explained by differences in procaspase 3 levels between these lines or by differences in the ability of mGrB or hGrB to cleave procaspase 3 either in vitro or in cell extracts (compare Fig. 2 B and Fig. S3). It is also unlikely that mGrB can circumvent Bid by directly and fully activating procaspase 3, as z-VAD blocks the second cleavage event when procaspase 3 is exposed to mGrB (Fig. 2 B and Fig. S3). Collectively, the aforementioned results indicate that mGrB and hGrB have similar but distinct substrate specificities and that they trigger death differently.
Substrate phage display confirms that hGrB and mGrB have distinct specificities
To compare the substrate specificities of hGrB and mGrB, we used substrate phage display. Phage carry a variable (randomized) nonameric peptide linked to a His tag, fused to a coat protein Phage adsorbed to nickel-agarose are treated with protease, and cleaved phage are analyzed after several rounds of capture and release.
The granzymes were used to probe a library in which the fourth residue of the nonamer had been fixed as Asp. As shown in Fig. 3, phage enriched by hGrB treatment displayed a P4-P3-P2-P1-P1′-P2′-P3′ consensus sequence, [I/V][G/E]ADVLV. This closely matches the P4-P1 sequence previously identified by combinatorial peptide analysis (Thornberry et al., 1997; Harris et al., 1998) and validates the phage display approach. In contrast, treatment with mGrB yielded the consensus sequence [I/L]X[F/Y]DXGV. Specifically, the preferences at the P4, P2, and P2′ positions are different between mGrB and hGrB. mGrB shows an equal preference for Ile or Leu at P4 and an aromatic residue at P2, whereas hGrB strongly prefers Ile at P4 and has a relaxed requirement at P2. mGrB also has a very strong preference for Gly at P2′, whereas hGrB has no obvious preference at this position.
To confirm that mGrB can cleave a substrate with the aforementioned characteristics, we synthesized a quenched-fluorescence substrate peptide based on one of the phage sequences that met the consensus [L/I]X[F/Y]DXGX. This peptide abz-LEYDLGALK(dnp)S was cleaved by both mGrB and hGrB, but kinetic analysis showed that it is a 370-fold better substrate for mGrB (Table I, compare specificity constants). These results show that the consensus substrate sequences for hGrB and mGrB differ at the key P4, P2, and P2′ positions and indicate structural differences in the corresponding pockets in the catalytic clefts of the two enzymes.
Efficient cleavage of Bid by mGrB requires a P2′ Gly
As shown in Fig. 2, mGrB cleaves Bid poorly but procaspase 3 relatively well. The substrate phage display results suggest that mGrB requires Gly at the P2′ to efficiently cleave substrates, but both mouse and human Bid have Glu at this position. To determine whether P2′ Gly enhances cleavage by mGrB, we synthesized a peptide substrate based on the GrB cleavage site in mouse Bid (IEPD↓SESQ) and a derivative with P2′ Gly (IEPD↓SGSQ). Both hGrB and mGrB cleaved the wild-type (wt) Bid substrate, but hGrB cleaved it >40-fold more efficiently (Table I). This is consistent with the 30-fold difference in the ability of the two proteases to cleave IETD-pNA (Fig. 2 A). In contrast, mGrB cleaved the P2′ Gly substrate sevenfold better than the wt Bid substrate, whereas hGrB showed a threefold improvement over cleavage of the wt sequence, consistent with the slight overrepresentation of Gly at P2′ in the phage display data (Fig. 3). Thus, mGrB prefers P2′ Gly in natural substrates, and its poor cleavage of Bid can be explained partly by lack of Gly at this position. Given the conservation of the P4-P1 sequence in Bid and procaspase 3 (IEXD), these results also explain why procaspase 3 is apparently a better substrate of mGrB than Bid (Fig. 2 B): mouse procaspase 3 has a P2′ Gly at its GrB cleavage site.
mGrB can be “humanized” by alteration of its S4/S3 subsite
The inability of mGrB to cleave after IETD suggested differences between mGrB and hGrB in the substrate binding cleft, presumably in the pockets that interact with the side chains of the P4, P3, or P2 residues of a substrate. Inspection of the structure of hGrB bound to the tetrapeptide aldehyde inhibitor IEPD-CHO (Rotonda et al., 2001) revealed that Asn218 in the S3 pocket forms hydrogen bonds with the P3 Glu of the substrate, whereas Tyr174 in the S4 pocket undergoes hydrophobic interactions with P4 Ile of the substrate (Fig. 4 A). These two residues are conserved in rat GrB and have been shown to influence substrate specificity (Ruggles et al., 2004). In contrast, they are not conserved in mGrB (Fig. 4 B, Arg174 and Lys218). Thus, we predicted that mutating mGrB toward hGrB by replacing Arg174 with Tyr and Lys218 with Asn would enable it to more efficiently cleave IETD and Bid.
As shown in Fig. 5 A, mutated mGrB retained the ability to cleave AAD-sbzl and showed substantially enhanced activity on IETD-pNA, with kcat/Km values on both substrates similar to those of hGrB (compare Fig. 5 A with Fig. 2 A). In addition, the mutant mGrB showed enhanced ability to cleave Bid, together with a slight increase in the ability to cleave procaspase 3 (Fig. 5 B). This translated into an increased ability to kill the Bid+/+ lymphomas (Fig. 5 C). Interestingly, the mutated mGrB also showed an increased ability to kill the Bid−/− lymphomas (Fig. 5 C), suggesting that rather than being simply converted to hGrB selectivity, it has gained the ability to cleave hGrB-specific targets (e.g., Bid) while retaining the ability to cleave mGrB-specific targets. Finally, mutant mGrB showed 100-fold enhanced ability to kill P815 cells, with similar efficiency to hGrB (Fig. 5 D, left).
At present, the reason for the difference in mutant mGrB killing of the lymphoma cells (two- to fourfold increase) compared with P815 cells (100-fold increase) is unclear, but it is possible that the altered mGrB cleaves a target in P815 cells that is absent from the lymphoma cells. Nevertheless, collectively, the aforementioned results confirm that mGrB has a substrate binding cleft that is structurally distinct from hGrB, whichexplains its different substrate specificity profile and killing mechanism.
Serpins that regulate hGrB and mGrB cytotoxicity are not conserved
hGrB is regulated by the intracellular protease inhibitor, PI-9, a member of the serpin superfamily (Bird et al., 1998; Hirst et al., 2003). The serpin inhibitory mechanism requires recognition of an exposed reactive center loop (RCL) as a potential substrate by a cognate protease, followed by cleavage between the P1 and P1′ residues. This triggers conformational change in the serpin, distorting the active site of the protease and trapping an intermediate consisting of the protease bound to the serpin (for review see Gettins, 2002). hGrB cleaves the PI-9 RCL after Glu, consistent with its preference for acidic residues, and is very rapidly trapped (Sun et al., 2001).
Mice have seven PI-9 homologues (Kaiserman et al., 2002), but the counterpart of PI-9 is not immediately obvious. The closest homologue, SPI6, has not been tested with mGrB but inhibits hGrB 20-fold slower than PI-9 (Sun et al., 1997). When overexpressed, SPI6 protects dendritic cells (Medema et al., 2001) and T cells (Phillips et al., 2004) from CL attack and increases the survival of CD8+ memory cells (Phillips et al., 2004). Loss of SPI6 results in diminished survival of CTL (Zhang et al., 2006). However, none of these studies has demonstrated that SPI6 directly or efficiently inhibits mGrB. To determine whether SPI6 is a mouse GrB inhibitor, we identified the cleavage site for mGrB in the SPI6 RCL and measured the kinetics of the interaction, assessed the ability of SPI6 and mGrB to interact in CL extracts, and tested whether SPI6 protects cells against death mediated by purified mGrB and perforin.
The hGrB–PI-9 interaction requires cleavage of the PI-9 RCL at Glu340 (Sun et al., 2001; Fig. 6 A), and mGrB also cleaves PI-9 at this point (unpublished data). N-terminal sequencing of the 6-kD fragment generated by treatment of SPI6 with mGrB or hGrB showed that cleavage occurs at Cys339, adjacent to the predicted P1 Cys340 (Fig. 6 A). This is a previously unrecognized cleavage specificity of GrB and broadens the range of potential substrates.
Although cleavage at Cys339 might allow capture of some mGrB by SPI6, the kinetics of the interaction will indicate whether mGrB is effectively regulated by SPI6 in vivo. Accordingly, we measured the rate constant (kass) and stoichiometry of inhibition of SPI6 and mGrB complex formation (Table II). The results strongly suggest that SPI6 is a poor inhibitor of mGrB. The stoichiometry of inhibition indicates that fewer than one out of five contacts results in mGrB trapping, and the low kass indicates that complex formation is very slow. SPI6 is even less effective in inhibiting hGrB (Table II). To put these results into context, PI-9 inhibits hGrB with equimolar stoichiometry and a rate of 1.2 × 106 M−1s−1 (Sun et al., 2001). Furthermore, the inhibition of mGrB by SPI6 is slower than the inhibition of hGrB by the viral serpin CrmA (Quan et al., 1995), which does not block GrB-mediated apoptosis in vivo.
Another signature of a serpin–protease interaction is the formation of SDS-stable complexes. For example, when human CLs are lysed, GrB is released from granules and rapidly binds cytosolic PI-9, and because GrB is in excess, essentially all of the PI-9 moves into complex (Fig. 6 B, left; Hirst et al., 2003). Similarly, in extracts of mouse CLs probed with anti-mGrB or anti-SPI6 antibodies, a complex containing mGrB and SPI6 was evident (Fig. 6 B, right). However, less than half of the SPI6 moved into complex, consistent with a less efficient interaction between this serpin and GrB, and indicating that the slow in vitro kinetics are not due to lack of a SPI6 cofactor present in CLs. To confirm that complex formation is GrB dependent, extracts were prepared from CLs taken from mice lacking GrB. No complex was evident in these cells (unpublished data).
If mouse and human GrB are orthologous, PI-9 and SPI6 should provide cross-species protection from GrB-mediated cell death. To test this, we challenged mouse P815 cells expressing SPI6 with hGrB and human SKW6.4 cells expressing PI-9 with mGrB (Fig. 6 C). As expected, PI-9 protected SKW6.4 cells against hGrB (Fig. 6 C, top left), increasing the EC50 ∼20-fold. SPI6 protected P815 cells against mGrB (Fig. 6 C, bottom right), increasing the EC50 >14-fold. Furthermore, the inactive SPI6T327R hinge mutant did not prevent mGrB-mediated death, indicating that the classic serpin inhibitory mechanism is involved in SPI6 cytoprotection (unpublished data). However, PI-9 did not protect cells against mGrB (Fig. 6 C, top right), nor did SPI6 effectively protect cells against hGrB (Fig. 6 C, bottom left).
Finally, the humanized mutant mGrB almost completely lost its ability to be regulated by SPI6, exhibiting a similar cytotoxicity profile to hGrB on SPI6-expressing cells (Fig. 5 D, right). This is consistent with the alterations to its active site cleft and substrate interactions and emphasizes the high specificity of a cognate serpin–protease interaction. From the aforementioned experiments, we concluded that mGrB is indirectly inhibited by SPI6 and that PI-9 and SPI6 are not functionally interchangeable. These results confirm that mGrB and hGrB have distinct substrate specificities and indicate that they are regulated differently in vivo.
Mouse and human GrA have significantly different cytotoxic potential
The functional differences that we identified between mGrB and hGrB raised the possibility that such interspecies differences would also be evident in other granzymes. Although mouse and human GrA cleave the peptide substrate benzyloxy-carbonyl-Lys-thiobenzyl ester with similar efficiency (Fig. 7 A; Odake et al., 1991), structural and substrate preference analysis of GrA has identified differences between the mouse and human enzymes in the active site cleft (Bell et al., 2003). This study identified an optimal P4-P1 substrate for hGrA as VANR↓, whereas the optimal mGrA substrate was identified as GYFR↓. Mutation of the hGrA S4 and S3 pockets converted it to a mouse-like GrA (Bell et al., 2003).
To determine whether these structural differences translate into differences in cytotoxic potential, we compared the ability of mGrA and hGrA to kill human or mouse cells in the presence of perforin. We used two different death assays: our standard MTT (3-[4,5-Dimethylthiazol-2-yl]-2,5-diphenyltetrazolium bromide) test, which measures loss of mitochondrial respiration (Fig. 7 C), and 51Cr release, which measures loss of membrane integrity (Fig. 7 D). As shown in Fig. 7 (C and D), mGrA efficiently killed mouse and human cells, but hGrA did not. mGrA-induced death was perforin dependent, as cells exposed to mGrA alone survived (Fig. 7 C). The inability of hGrA to kill cells was not related to the source of the protease, as both recombinant and native forms were used. All of the hGrA used was in the physiological 60-kD dimeric form, which could be reduced to the 30-kD monomeric form (Fig. 7 B). Monomeric hGrA also failed to kill cells (unpublished data).
GrA and GrB are believed to be the primary cytotoxins in human and mouse CLs. The evidence presented here clearly demonstrates substantial interspecies distinctions that suggest caution when using mouse models to elucidate the molecular and pathophysiological functions of these proteases in humans. For example, loss of GrB may be much more devastating to humans than the mild phenotype of GrB-deficient mice would suggest, given that hGrB is more efficient than mGrB and that hGrA is apparently nonlethal and could not compensate for GrB absence.
Such interspecies differences between granzymes are illustrated by the GrB–Bid interaction, which indicates that distinct death pathways are initiated by hGrB and mGrB (Fig. 8 A). Our results show that hGrB is a much more effective killer than mGrB because it can cleave and activate Bid, whereas mGrB cleaves Bid poorly and does not require it for initiating apoptosis. The requirement for Bid in hGrB-mediated killing is not an artifact arising from its use on mouse target cells, as Bcl-2 overexpression in human cells blocks hGrB cytotoxicity but can be relieved by Bid coexpression (Waterhouse et al., 2006). The failure of hGrB to kill in the absence of Bid indicates that it lacks the capacity to independently and fully activate effector caspases, as suggested previously (Goping et al., 2003; Sutton et al., 2003). Its ability to cleave other substrates must also be insufficient to cause death. In contrast, mGrB probably directly activates another effector caspase, such as procaspase 7, or it relieves IAP inhibition of caspase 3 independently of Bid. It may also cause caspase-independent death by cleaving substrates such as ICAD or as-yet-undefined mitochondrial targets (Thomas et al., 2001).
The differences between hGrB and mGrB are also emphasized by our ability to humanize mGrB and improve its capacity to cleave Bid. However, humanized mGrB remains relatively inefficient compared with hGrB, indicating that additional changes to the mGrB substrate binding pocket are required to fully accommodate Bid. Supporting evidence for differences between mGrB and hGrB can also be found in previous studies. For example, synthetic catalytic inhibitors designed to hGrB do not inhibit mGrB (Willoughby et al., 2002); the human adenoviral protein L4-100K inhibits hGrB but not mGrB (Andrade et al., 2003); mGrB does not process human Bid (Adrain et al., 2005); and mGrB does not require Bax or Bak for cytotoxicity, whereas hGrB does (Thomas et al., 2001; Cartron et al., 2003).
Similarly, hGrA and mGrA do not appear to be functionally equivalent, as we have shown that mGrA is cytotoxic and that hGrA is not. This is supported by several other lines of evidence. Modeling of mGrA on the structure of hGrA has identified differences in their active site clefts, substantiated by differences in peptide substrate preference and the conversion of hGrA to an mGrA-like protease by mutation of the S4/S3 pockets (Bell et al., 2003). Early work showed that mGrA is cytotoxic when produced with perforin in rat basophilic leukemia cells (Shiver et al., 1992), and many studies on the mechanisms of GrA-mediated death have used CTL and target cells from mice (Pardo et al., 2002, 2004). In contrast, recent studies using specific inhibitors show that hGrA is a minor contributor to CL cytotoxicity (Mahrus and Craik, 2005). Such findings do not necessarily overturn the idea that hGrA engages apoptotic pathways by cleaving specific subcellular proteins and generating reactive oxygen species (Pardo et al., 2004; Martinvalet et al., 2005), but they suggest that it does so much less efficiently than mGrA and that it is not sufficient to induce cell death.
The differences between hGrB and mGrB also extend to their regulation, as indicated by differences between PI-9 and SPI6. PI-9 is an efficient hGrB inhibitor that protects CLs and bystanders from misdirected GrB (Sun et al., 1996; Bird et al., 1998; Hirst et al., 2003), and SPI6 is a putative counterpart of PI-9 (Sun et al., 1997; Medema et al., 2001; Phillips et al., 2004). Mice lacking SPI6 have fewer CTL, and remaining CTL have fewer granules and less evidence of cytoplasmic mGrB (Zhang et al., 2006). Crossing SPI6-deficient mice with mGrB-deficient mice rescues the CTL defect, but because the latter mice also lack several mGrB paralogues (Pham et al., 1996), it cannot be assumed that SPI6 directly controls mGrB. Although the accumulating evidence is consistent with a cytoprotective role for both PI-9 and SPI6 in CLs, our present work strongly suggests that the question of the mechanism of cytoprotection by SPI6 remains open, as it is not an effective mGrB inhibitor. Given that other mGrB paralogues are cytotoxic, four possibilities exist: the primary target of SPI6 is a cytotoxic mGrB paralogue; SPI6 simultaneously controls multiple mGrB paralogues through a one-size-fits-all RCL; mGrB or paralogues activate a downstream apoptotic protease, which is efficiently inhibited by SPI6 (and not activatable by hGrB); or SPI6 protects granule integrity independently of granzyme inhibition (in which case, PI-9 would be predicted to control GrB and independently protect granule integrity). We consider the last two possibilities more likely because they explain why the rat, which has a more human-like GrB (Harris et al., 1998), has an orthologue of SPI6 rather than PI-9.
Among mammals, the multiple granzyme and serpin genes found in rodents appear to be the exception rather than the rule (Fig. 8 B), and it is evident that GrB and PI-9 paralogues emerged after divergence of Lagomorpha and Rodentia. A priori, the disappearance of the PI-9 RCL sequence from rodents suggests a fundamental change to the structure of mGrB, which we have confirmed by demonstrating that mGrB is catalytically distinct from hGrB and is not controlled by PI-9. Conversely, we have shown that hGrB is not controlled by SPI6. Although mGrB has changed, the GrB cleavage site in mouse Bid is conserved and is accessible to hGrB. Thus, it can be surmised that mGrB has lost cytotoxicity rather than hGrB gaining it. Conservation of the mouse Bid cleavage site suggests that another protease attacks it or that it serves another function in the molecule.
Why have these changes occurred in rodents? Given its role in immunity, the simplest explanation is that mGrB has evolved to target a virulence factor from a particular pathogen or to avoid inhibition by a specific pathogen protein. Specific targeting of virulence factors has been demonstrated for another leukocyte protease, neutrophil elastase (Weinrauch et al., 2002). Duplication of mGrB to generate paralogues may have allowed compensating cytotoxic mechanisms to evolve and/or multiple virulence factors to be targeted. The former possibility is supported by the reported cytotoxic functions of GrC and GrF (Revell et al., 2005). The appearance of cytotoxic mGrB paralogues may have then driven the evolution of cognate PI-9 paralogues.
Although functioning primarily as a cytotoxin, GrB probably has other roles. For example, in humans it is found in testis and has been implicated in reproduction (Hirst et al., 2001). It may also be involved in extracellular matrix remodeling (Buzza et al., 2005). With multiple substrates and roles, it is easy to envisage that particular functions of GrB may have devolved to separate mGrB paralogues in rodents. For example, mouse granyzme N appears to be restricted to the testis and may have assumed the testicular role of hGrB (Takano et al., 2004).
In closing, we have shown that mGrB is structurally and enzymatically distinct from hGrB and has a different killing mechanism. We have also shown that human GrA is not cytotoxic, whereas mGrA is substantially cytotoxic. Such species differences may also extend to granzymes M and K, but appropriate comparisons have yet to be made. It may be that the overall killing efficiency of human and mouse CLs is comparable and that the loss of GrB cytotoxicity in the mouse has been offset by a gain in GrA cytotoxicity or the evolution of GrB paralogues. Nevertheless, use of a granzyme from one species on cells or extracts from another species may produce misleading results because key substrates are missing or are inefficiently cleaved because recognition sites are not conserved. Use of mice to test inhibitors of human granzymes will fail because differences in active site topology of the mouse enzymes will prevent inhibitor binding. Likewise, use of mice to study viral infection may be compromised, as conserved virus inhibitors (e.g., SPI2/CrmA, Bcl-2 homologues, and AdL4-100K homologues) will interact differently with or fail to control mouse granzymes. Finally, use of mouse models to study the impact of granzymes on tumorigenesis is also problematic because the species-specific death pathways and substrates outlined here imply that tumors will subvert the immune system differently in mice and humans by up- or down-regulating different proteins.
Materials And Methods
Granzymes and serpins
Recombinant serpins were produced in Pichia pastoris as described previously (Sun et al., 1995, 1997). Recombinant granzymes were also produced from P. pastoris (Sun et al., 1999, 2004). To improve expression levels, the wt mGrB cDNA was modified by QuikChange site-directed mutagenesis (Stratagene) to remove a cryptic polyadenylation site, using the oligonucleotide 5′-CACTGTGAAGGAAGTATAATCAATGTCACTTTGGGGGC-3′ and its complement. The S4/S3 subsites in mGrB were mutated using the oligos 5′-GAGTCCTACTTTAAAAATTATTACAACAAAACCAATCAG-3′ and its complement (R174Y) and 5′-TCCTATGGATATAATGATGGTTCACCTCC-3′ and its complement (K218N). Activity of recombinant GrB was assessed by cleavage of the peptide substrate benzyloxycarbonyl- AlaAlaAsp-thiobenzyl ester (Sun et al., 1999, 2004). The specific activity of hGrB was 90 U/μg, and the specific activity of mGrB was 120 U/μg. The proportion of active enzyme in each batch of purified protein was estimated by the ability to complex with PI-9 or SPI-6. Routinely, >90% of granzyme in a batch could complex with the appropriate serpin. Native human GrA was from lymphokine-activated killer cells. Activity of GrA was assessed by cleavage of the substrate benzyloxycarbonyl- Lys-thiobenzyl ester.
Kinetic analysis of substrate cleavage by granzymes and inhibition by serpins
Quenched fluorescence substrates were synthesized and used as described previously (Sun et al., 2001). Activity of granzymes on the peptide substrates Boc-AAD-thiobenzyl ester (BIOMOL Research Laboratories, Inc.) and IETD-pNA (Calbiochem) was assessed as described previously (Sun et al., 2001). Inhibition of granzymes by serpins was measured using standard approaches (Hopkins and Stone, 1995; Sun et al., 2001). Granzyme cleavage sites in SPI6 and PI-9 were identified by N-terminal sequencing (Sun et al., 2001).
Production of granzyme substrates via in vitro transcription and translation
35S-labeled mouse procaspase 3 or mouse Bid was produced from cDNAs in the expression vector pSVTf via in vitro transcription and translation (Sun et al., 2004). Cleavage by granzymes was assessed by adding a specified amount of protease directly to the translation mix and incubating for 30 min at 37°C. The caspase inhibitor z-VAD was obtained from Calbiochem. Products were separated by SDS-PAGE and visualized by fluorography.
Analysis of granzyme substrate specificity by phage display
Construction of substrate libraries based on bacteriophage T7 used the T7Select1-1b phage display system (Novagen). Each library was made by synthesizing a degenerate oligonucleotide, annealing it to complementary half-site oligonucleotides, ligating the resulting heteroduplex to vector arms, and adding it to a T7 phage packaging extract. The half-site oligonucleotides were 5′-GCCGCCTGGAGTGAGAG-3′ and 5′-AGCTTAGTGATGGTGATGGTGATG-3′. One library was constructed using the degenerate oligonucleotide 5′-AATTCTCTCACTCCAGGCGGC(NNK)9CATCACCATCACCATCACA-3′ (where N represents any nucleotide and K represents T/C). This added a randomized nonameric peptide and a hexahistidine tag to the C terminus of the 10B coat protein. The complexity of this “random” library was 7 × 106 plaque-forming units. Another library was constructed using 5′-AATTCTCTCACTCCAGGCGGC(NNK)3GAC(NNK)5CATCACCATCACCATCACA-3′. This encoded an aspartate near the middle of the randomized peptide. The complexity of this “fixed” library was 4 × 106 plaque-forming units. Each library was amplified to >1010 plaque-forming units. About 109 plaque-forming units of amplified phage were bound to nickel-chelated Sepharose beads at 4°C. Unbound phage was removed by washing the beads in phosphate-buffered saline containing 850 mM NaCl and 0.1% (vol/vol) Tween 20. After two further washes in phosphate-buffered saline containing 1 mM MgSO4, the suspension was split into two parts (treatment and control; control treatments were performed to assess the extent of phage release from beads in the absence of protease). 200 nM granzyme was added to the treatment tube, and both tubes were incubated overnight at 37°C. Plaque-forming units in the supernatants representing cleaved or released phage were then counted and amplified to form sublibraries for the next round of selection. Phage remaining bound to the beads was eluted with 0.5 M imidazole, and plaque-forming units were counted to assess cleavage efficiency. After several selection rounds, individual plaques were chosen at random for sequence analysis. Phage DNA was amplified by PCR using dedicated primers (T7Select cloning kit; Novagen). Sequencing of PCR products using the same primers was performed using the Big Dye 3.1 kit (GE Healthcare).
The sequencing results were analyzed to determine the statistical distribution of each amino acid at each position of the nonamer (Matthews et al., 1994), allowing for the redundancy of the code, the fact that only 32 out of 64 codons are represented by NNK, and the exclusion from the analysis of any sequences encoding a stop codon in the nonamer. In a binomial distribution of amino acids, Δσ yields the difference of the observed frequency from the expected frequency in terms of standard deviations:where Obs(x) is the number of times amino acid x occurs in the selected sequences, P(x) is the theoretical probability of amino acid x occurring, and n is the total number of sequences analyzed.
Cells, cell culture, and transfections
Human YT and Jurkat cells were maintained as described previously (Bird et al., 1998). Eμ-myc/bid−/− and Eμ-myc/bid+/+ mouse B cell lymphomas (Waterhouse et al., 2005) were obtained from R. Johnstone (Peter MacCallum Cancer Institute, Melbourne, Australia). Human SKW6.4 and mouse EL4 and YAC-1 cells were maintained in Dulbecco's modified Eagle medium containing 10% heat-inactivated fetal calf serum, 2 mM glutamine, and 55 μM β-mercaptoethanol. P815 cells were maintained in RPMI 1640 medium containing 10% heat-inactivated fetal calf serum, 2 mM glutamine, and 55 μM β-mercaptoethanol. Isolation and culture of mouse splenocytes followed standard procedures.
To generate SKW6.4 cells stably expressing PI-9, a PI-9 expression vector (Bird et al., 1998) was cotransfected by electroporation with a plasmid expressing neomycin resistance, and selection was performed in 1 mg/ml G418 (Sigma-Aldrich). Clones were established by limiting dilution and screened for PI-9 expression by indirect immunofluorescence and immunoblotting. FLAG-tagged SPI6 and SPI6T327R (containing an inactivating mutation in the proximal hinge) cDNAs were generated by PCR from the appropriate templates using the primer pair 5′-CGGGATCCATGGACTACAAAGACGATGACGATAAAGGGAATACTCTGTCTGAAG-3′ and 5′-GCTCTAGATTATGGAGATGAGAACCTGCCACA-3′. Products were cloned into the BamHI and XbaI sites of the expression vector pEF-PGKpuropA (Huang et al., 1997). 1.25 × 107 P815 cells were electroporated at 330 μF and 350 V with either pEFpuro/FlagSPI6 or pEFpuro/FlagSPI6T327R vectors linearized with SalI. Transfected pools were grown for 3 d in 1.1 μg/ml puromycin, and subsequently 45 SPI6 and 28 SPI6T327R lines were cloned by limiting dilution in puromycin selection. Clones were scored for relative SPI6 expression by indirect immunofluorescence using the α-FLAG M2 antibody. Expression was confirmed by immunoblotting using rabbit antiserum (R25) raised to recombinant SPI6. This antiserum does not recognize the highly related mouse serpins (Serpinb9b, Serpinb9e, Serpinb6, Serpinb6b, and Serpinb1) and can be considered specific for SPI6.
Cell death mediated by granzymes and recombinant perforin was assessed as previously described (Sun et al., 2004; Bird et al., 2005). Recombinant perforin was produced using a baculovirus expression system (Liu et al., 1994). Recombinant SLO was prepared and used according to Walev et al. (2001).
Preparation of cell extracts, antibodies, and immunoblotting
For analysis of serpin complexes, cells were lysed in either 1 volume of Laemmli sample buffer (LSB) or 1/2 volume of NP-40 lysis buffer (50 mM Tris HCl, pH 8.0, 10 mM EDTA, and 1% [vol/vol] Nonidet P40), followed by incubation at 37°C for 10 min, after which 1/2 volume of 2× LSB was added. Viscosity was reduced by mechanically shearing the DNA using a needle and syringe. For analysis of caspase 3 cleavage, cells were lysed in NP-40 lysis buffer containing 1 μg/ml pepstatin, 1 μg/ml leupeptin, 1 μg/ml aprotinin, and 10 μg/ml PMSF, and a postnuclear supernatant was prepared by centrifugation at 16,000 g. Antibodies to PI-9 (R15) and hGrB (2C5 hybridoma supernatant) were used for immunoblotting at 1:2,000 and 1:100, respectively. A rat monoclonal antibody against mGrB was obtained from eBioscience (clone 16G6) and was used at 1:1,000. Goat anti-actin sera (Santa Cruz Biotechnology, Inc.) was used at 1:1,000. Rabbit antisera to SPI6 (R25) and hGrA (R045) was raised against recombinant protein purified from P. pastoris following standard procedures and was used at 1:1,000. Rabbit antisera to caspase 3 (Cell Signaling Technology) was used at 1:100. HRP-conjugated secondary antibodies was used at 1:5,000 and detected using enhanced chemiluminescence.
Image acquisition and manipulation
Images of blots were initially captured on x-ray films, which were subsequently scanned into Photopaint (Corel). Adjustments to brightness or contrast of digital images were applied to the whole image. No nonlinear adjustments were made.
Online supplemental material
Fig. S1 shows a comparison of methods used to assess cell survival after GrB treatment. Fig. S2 shows a comparison of uptake of mGrB and hGrB into target cells. Fig. S3 shows GrB-mediated procaspase 3 cleavage in Bid+/+ or Bid−/− mouse B cell lymphoma extracts. Online supplemental material is available at http://www.jcb.org/cgi/content/full/jcb.200606073/DC1.
We thank K. Sedelies for assistance with 51Cr assays, Q. Zhou for help with construction of mutant mGrB, Dr. R. Johnstone for providing the Eμ-myc cells, and Dr. N. Waterhouse (Peter McCallum Cancer Institute) for advice and discussions.
This work was supported by program grants to P.I. Bird, J.C. Whisstock (284233), and J.A. Trapani from the National Health and Medical Research Council (Australia).
D. Kaiserman and C.H. Bird contributed equally to this paper.
Abbreviations used in this paper: CL, cytotoxic lymphocyte; IAP, inhibitor of apoptosis protein; MTT, 3-(4,5-Dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide; RCL, reactive center loop; SLO, streptolysin O; wt, wild-type.
- Submitted: 14 June 2006
- Accepted: 19 October 2006