|
|
||||||||
-2 Domain of Two Naturally Occurring Equine MHC Class I Molecules Alters the Recognition of Gag and Rev Epitopes by Equine Infectious Anemia Virus-Specific CTL1
* Department of Veterinary Microbiology and Pathology, Washington State University, Pullman, WA 99164; and
Bioinformatics and Computational Biology Program, Iowa State University, Ames, IA 50011
| Abstract |
|---|
|
|
|---|
-2 domain. This substitution did not affect recognition of Env-RW12, but resulted in more efficient recognition of Rev-QW11. Significantly, CTL recognition of Gag-GW12 was abrogated, despite Gag-GW12 binding to 141. Molecular modeling suggested that conformational changes in the 141/Gag-GW12 complex led to a loss of TCR recognition. These results confirmed that the ELA-A1 haplotype is comprised of functionally distinct alleles, and demonstrated for the first time that naturally occurring MHC class I molecules that vary by only a single amino acid can result in significantly different patterns of epitope recognition by lentivirus-specific CTL. | Introduction |
|---|
|
|
|---|
Equine infectious anemia virus (EIAV)3 is a macrophage-tropic lentivirus that causes persistent infections in horses worldwide. In contrast to HIV-1 infection however, most EIAV-infected horses eventually control plasma viremia and clinical disease and remain life-long inapparently infected carriers (13, 14, 15). The initial as well as the long-term control of viremia and clinical disease is the result of adaptive immune responses, including neutralizing Ab and importantly, CTL (16, 17, 18, 19, 20, 21, 22). Due to these robust viral-specific immune responses that contain viral replication, EIAV infection in horses is a unique and useful model system for the study of lentiviral immune control.
As in HIV-1, Env, Gag, and Rev proteins are important CTL targets in EIAV-infected horses. EIAV Env- and Gag-specific CTL are detected during acute and inapparent infection (17, 23, 24), with Gag-specific CTL responses occurring in the majority of horses tested (25, 26). Although EIAV Env-specific CTL can be immunodominant and may be important in early viral control, viral escape can limit the effectiveness of CTL directed against variable Env epitopes, and Gag-specific CTL frequency can correlate with the ability to control viral load and clinical disease (27, 28). Epitope clusters occur within EIAV Gag proteins that are recognized by CTL (including high-avidity CTL) in horses with disparate MHC class I (MHC I) haplotypes and are likely important in control of clinical disease (29, 30). Additionally, moderate and high-avidity Rev-specific CTL are associated with control of viremia and disease in EIAV-infected nonprogressor horses (27). Significantly, no viral escape from high-avidity Rev-specific CTL has been observed. The Rev epitope recognized by these CTL is highly conserved among other strains of EIAV, and independent studies by other investigators show that this sequence does not change during long-term EIAV infection (31, 32).
Given that CTL epitopes in Env, Gag, and Rev are critical for lentiviral immune control, factors that affect the presentation and recognition of these peptides by CTL in a population of infected individuals are important considerations for vaccine development. One of the most important factors is MHC I polymorphism, because CTL recognize peptide epitopes only in the context of allelic forms of MHC I molecules. The CTL TCR binds to a cell surface MHC I complex consisting of the MHC I molecule, a peptide processed from a viral protein, and β2-microglobulin (β2m) (33). The HLA classical MHC I genes are the products of three loci, and as of July 2006, there were 478 HLA-A, 805 HLA-B, and 256 HLA-C named alleles (
www.ebi.ac.uk/imgt/hla/intro.html
) (34). MHC I polymorphism occurs primarily in the residues that form the peptide-binding domain, therefore influencing the types of peptides that are selectively bound. Similarities in peptide-binding specificity between human MHC I molecules have been identified based on the structure of peptide-binding pockets, peptide-binding assays, and analysis of motifs, and sets of molecules with similar specificities are called supertypes (35, 36, 37, 38). Despite the fact that different MHC I molecules can belong to a supertype and have similar peptide-binding specificities, differences in only a few residues among class I molecules can significantly affect the recognition of the MHC I-peptide complex by CTL. Saturation mutagenesis of a murine MHC I molecule has shown that a single amino acid substitution in the
-1 domain can result in more effective killing of lymphocytic choriomeningitis virus-infected target cells by lymphocytic choriomeningitis virus-specific CTL (39). Importantly, the HLA-B*35 subtypes HLA-B*3503 and HLA-B*3502 (B*35-Px genotype) correlate more strongly with rapid progression to AIDS than does the HLA-B*3501 subtype (B*35-PY genotype), which varies from HLA-B*3503 and HLA-B*3502 by only one and two amino acids, respectively, in the peptide-binding domain (40). It is presumed that this effect is due to differential presentation of epitopes to HIV-1-specific CTL, and it has been demonstrated that a higher frequency of Gag-specific CTL correlate with lower viral loads in individuals with the B*35-PY genotype, whereas no significant relationship exists between CTL activity and viral load in the B*35-Px group (41).
Compared with humans, less is known regarding the numbers and locus assignments of equine MHC I alleles. Serology has been the most widely used MHC I typing method in horses, with 17 equine leukocyte Ag (ELA)-A haplotypes defined (42, 43, 44, 45). More recently, classical MHC I alleles have been identified and sequenced, but it has not been possible to assign alleles to loci (46, 47, 48, 49, 50, 51, 52). Recent work indicates that up to three or four (or more) classical MHC I loci exist, and it is likely that the number of loci is variable and dependent on haplotype (48, 52). Regardless, we have identified the classical equine MHC I molecule 7-6, which is associated with the serologically defined ELA-A1 haplotype (51). The 7-6 molecule presents Env and Gag epitopes (Env-RW12 and Gag-GW12) to distinct populations of EIAV-specific CTL (27, 51). Target cells from some ELA-A1 horses present both epitopes, whereas target cells from other ELA-A1 horses present only Env-RW12 (51), suggesting that subtypes that comprise functionally distinct alleles exist within the ELA-A1 haplotype. In addition, a CTL epitope in EIAV Rev (Rev-QW11) is also restricted by the ELA-A1 haplotype (27), but the molecule that presents Rev-QW11 has not been identified. The purpose of the present study was to determine whether different horses sharing the ELA-A1 haplotype used functionally distinct MHC I molecules to present ELA-A1-restricted CTL epitopes. To test this hypothesis, we sequenced and expressed MHC I genes from three different horses with the ELA-A1 haplotype and determined how these different molecules affected the recognition of important Env, Gag, and Rev epitopes by EIAV-specific CTL.
| Materials and Methods |
|---|
|
|
|---|
Arabian horses A2140, A2150, and A2152 were used in this study. A2152 is a noninfected 7-year-old breeding stallion in which the 7-6 MHC I allele was initially identified (51). A2140 and A2150 are 8- and 7-year-old mares, respectively, that have been infected with EIAVWSU5 for 7 and 6 years, respectively (27). The ELA-A haplotypes were determined serologically by lymphocyte microcytotoxicity (42, 53, 54) using reagents provided by Dr. E. Bailey (University of Kentucky, Lexington, KY). All experiments involving horses were approved by the Washington State University Institutional Animal Care and Use Committee.
Identification of equine MHC I alleles
The full-length 141 gene (GenBank accession no. AY374512) from A2150 was cloned and sequenced as previously described (51). Briefly, PBMC were isolated and cultured for 48 h in RPMI 1640 medium with 10% FBS, 2 mM L-glutamine, and 10 µg/ml gentamicin. Con A (40 µg/ml) was added to up-regulate MHC I expression. mRNA isolated from these cells was used for first-strand synthesis, which was primed with 1 µg of oligo(dT)12–18 and 50 ng of random hexamers. This reaction was incubated at an initial temperature of 37°C for 15 min, followed by an additional hour at 49°C. The second-strand reaction was incubated at 16°C for 2 h using Escherichia coli DNA ligase, polymerase, and RNase H. cDNA was blunt-ended with T4 DNA polymerase, EcoRI adapters were added with T4 ligase, and the adapters were phosphorylated with T4 polynucleotide kinase. The cDNA was then size fractionated, and cDNA between 2 and 3 kb was ligated into EcoRI-digested and dephosphorylated pcDNA3 vector (Invitrogen Life Technologies). E. coli electrocompetent cells (Invitrogen Life Technologies) were transformed with this ligation mixture resulting in a cDNA library. Clones were selected from the library by colony-lift hybridization with a 32P-labeled HindIII-XbaI fragment from horse MHC I gene 8/9 (provided by Dr. D. Antczak, Cornell University, Ithaca, NY) (46). Positive colonies were isolated, and those with inserts of the correct size were sequenced. Sequencing was performed at the Laboratory for Biotechnology and Bioanalysis (Washington State University) using dye-labeled dideoxynucleotide-cycle sequencing with an ABI 377 automated sequencer.
The presence of the 7-6 allele (GenBank accession no. AY225155) was confirmed in A2140 using a RT-PCR as previously described (48) with modifications. Briefly, total RNA was isolated from 2 x 107 equine kidney (EK) cells with an RNeasy mini kit (Qiagen). Ten microliters of RNA and a SuperScript One-Step RT-PCR kit (Invitrogen Life Technologies) was then used with 1 µM each of forward (5'-ATG ATG CCC CCA ACC TTC-3') and reverse (5'-TGA ACA AAT CTT GCA TCA CTT G-3') primers. The reaction conditions were 45°C for 50 min, followed by 94°C for 2 min, then 35 cycles of 94°C for 30 s, 53°C for 30 s, and 72°C for 1 min. The 1116-bp product was gel isolated and extracted (Qiagen), then used for TOPO TA cloning into the pCR4-TOPO vector for sequencing (Invitrogen Life Technologies). Isolated colony minipreps were digested and electrophoresed to screen for inserts and those with inserts of correct size were sequenced by the Laboratory for Biotechnology and Bioanalysis using dye-labeled dideoxynucleotide-cycle sequencing with an ABI 377 automated sequencer.
Although we independently named the alleles in this study, official Immuno Polymorphism Database nomenclature (
www.ebi.ac.uk/ipd/
) will soon be available for equine MHC alleles (55).
Expression of equine MHC I alleles
Two retroviral vectors were constructed as described (51) using the plasmid pLXSN (provided by Dr. A. Dusty Miller, Fred Hutchinson Cancer Research Center, Seattle, WA). The MHC I genes 7-6 and 141 were PCR amplified and ligated into the cloning site of pLXSN downstream of the Moloney murine sarcoma virus long terminal repeat and upstream of the neomycin phosphotransferase gene, which was under the control of the SV40 early promoter (56). The sequences of the inserts and flanking plasmid DNA were determined. To generate vector-producing cell lines, published procedures were used (51, 57). Briefly, an amphotropic packaging cell line, PA317 (CRL-9078; American Type Culture Collection) was transfected with each plasmid. Supernatant from PA317 cells was used to transduce amphotropic PG13 packaging cells (CRL-10686; American Type Culture Collection), which were then selected using 750 µg/ml G-418 sulfate (Invitrogen Life Technologies). Vectors were harvested from the selected PG13 cells and used to transduce CTL target cells. EK cells and human mutant B lymphoblastoid 721.221 cells (58) were transduced with the retroviral vectors expressing 7-6 and 141, pulsed with peptides, and used as CTL targets (51). The 721.221 cells (obtained from Dr. A. Sette, La Jolla Institute for Allergy and Immunology, San Diego, CA) express human β2m, but not HLA-A, -B, or -C class I molecules (58).
PBMC stimulations and CTL assays
PBMC stimulations and CTL assays were performed as described (17, 27, 28, 51) with modifications. Briefly, PBMC were isolated from A2140 and A2150 and stimulated with peptide-pulsed autologous monocytes. EK target cells from A2140 and A2150 and mixed-breed pony H585 were established from kidney tissue obtained by biopsy (17). For stimulation with peptides, 2 µM Env-RW12, Gag-GW12, or Rev-QW11 was added to PBMC in 10% FBS. Peptide and PBMC were incubated for 2 h at 37°C with occasional mixing before centrifugation at 250 x g for 10 min. PBMC were resuspended to 2 x 106/ml in RPMI 1640 medium with 10% FBS, 20 mM HEPES, 10 µg/ml gentamicin, and 10 µM 2-ME. One milliliter was added to each well of a 24-well plate and incubated for 1 wk at 37°C before use in CTL assays. CTL activity was measured using a 51Cr release assay with a 17-h incubation period using EK target cells (17, 27, 28). In addition, human B lymphoblastoid 721.221 target cells transduced with retroviral vectors expressing equine MHC I genes 7-6 or 141 as described above were used in assays with a 5-h incubation period (51). The shorter incubation period was used because 721.221 target cells are less hardy than EK target cells and they develop spontaneous lysis sooner (51). Target cells were pulsed with various amounts of peptide Env-RW12, Gag-GW12, and Rev-QW11 as indicated in the figures. The formula, percent-specific lysis = [(E – S)/(M – S)] x 100, was used, where E is the mean of three test wells, S is the mean spontaneous release from three target cell wells without effector cells, and M is the mean maximal release from three target cell wells with 2% Triton X-100 in distilled water. The E:T ratio was 20:1 or 50:1 as indicated in the figures, and each well contained
30,000 target cells. The 50:1 E:T ratio was used to confirm 20:1 E:T ratio results. Comparisons were only made between assays that used the same E:T ratio. Previous work indicates that these E:T ratios yield consistent results corresponding to the log portion of the killing curve, and that the 50:1 E:T ratio results in the highest percent-specific lysis (17, 19, 24, 25, 26, 27, 28, 29, 30, 51, 57, 59, 60, 61, 62). Assays with similar constant E:T ratios have been used by others (63). Only assays with a spontaneous target cell lysis of <30% were used. The SE of percent-specific lysis was calculated using a formula that accounts for the variability of E, S, and M (64). Significant lysis was defined as the percent-specific lysis of peptide-pulsed target cells that was >10% and also >3 SE above the nonpulsed target cells or above target cells transduced with control vectors and pulsed with the relevant peptide. For comparisons of CTL recognition efficiency, the peptide concentration that resulted in 50% maximal target cell-specific lysis (EC50) was used. The EC50 was calculated after transforming the percent-specific lysis data to percent-maximal lysis (with the lowest percent-specific lysis value set to 0% and the highest percent-specific lysis value set to 100%) and fitting the curve with nonlinear regression using GraphPad Prism version 3.03 (GraphPad). This is an established method to measure and compare CTL recognition efficiency and avidity (27, 30, 65, 66). All CTL assays were performed at least twice, and the results were consistent in each case.
Live cell peptide-binding assay
Peptide binding to equine MHC I molecules 7-6 and 141 was measured as previously described (51), with slight modifications, using the chloramine-T method (67, 68). Briefly, 721.221 cells transduced with MHC I gene 7-6, 141, or with a retroviral vector that did not express an equine gene, were preincubated with human β2m. The cells were washed and resuspended in RPMI 1640 containing β2m, EDTA, PMSF, and N
-p- tosyl-L-lysine chloromethyl ketone. One hundred microliters containing 2 x 106 cells plus 1 µl containing 1.5 x 105 cpm of 125I-labeled Env-RW12 was incubated for 4 h at 22°C in wells of a 96-well U-bottom plate. Cells were then washed three times with serum-free medium, centrifuged through calf serum to remove any remaining unbound radiolabeled peptide, and then washed a final time. The radioactivity of the cell pellet was counted with a gamma scintillation counter (Packard Instrument). Competitive inhibition assays were performed twice in triplicate by adding 10 µl containing sufficient unlabeled peptide competitors to result in final concentrations of 1–1000 nM to the initial mixture of cells before addition of radiolabeled Env-RW12 peptide. Competing peptides were unlabeled Env-RW12, Gag-GW12, Rev-QW11, and control peptide 1b4a (VRVEDVTNTAEY), which does not inhibit the binding of Env-RW12 to 7-6 (51). For each competing peptide, the concentration resulting in 50% inhibition of radiolabeled Env-RW12 binding (IC50) was calculated by fitting the curve with nonlinear regression using GraphPad Prism version 3.03. All peptides used in binding assays had 95% purity and were synthesized by Sigma-Genosys.
Molecular modeling
Because crystal structures of equine MHC I molecules are not available, three-dimensional computer models of 7-6 and 141 were generated based on known structures of human MHC I molecules using MODELLER 8v2 (
www.salilab.org/modeler
) (69). Templates for modeling were chosen based on PSI-BLAST searches of the Brookhaven Protein Data Bank database. The 1XR9A structure was used for 7-6 and the 1ZSDA structure (70) was used for 141. The models were verified using VERIFY3D (
http://nihserver.mbi.ucla.edu/Verify_3D
) (71). Models of the Env-RW12, Gag-GW12, and Rev-QW11 peptides were also generated. To predict the side chain conformations for each peptide, SCWRL3.0 (
http://dunbrack.fccc.edu/SCWRL3.php
) (72) was used. Viral peptides bound to human MHC I molecules served as templates and were chosen from the Brookhaven Protein Data Bank database (1ZHKC for Env-RW12 and Gag-GW12, and 1ZSDC for Rev-QW11). The 1ZSDC structure is an 11-mer peptide (70), like Rev-QW11. Because no structures of 12-mer peptides bound to MHC I molecules were available, the first residue (L) of 1ZHKC, a 13-mer peptide (73), was eliminated before use as the backbone template for Env-RW12 and Gag-GW12. To model the binding of the three peptides to 7-6 and 141, the docking program FTDock2.0 (
www.bmm.icnet.uk/docking
) (74, 75, 76) was used, and the docking score (RPScore) for each complex was determined (76). The interacting residues for each complex were identified, and the interactions by category (hydrophobic, salt bridges, repulsive charged, hydrogen bonds, and aromatic stacking) between atoms of contact residues were determined using the STING Millennium Suite program (
http://trantor.bioc.columbia.edu/SMS/index_m.html
) (77, 78). Finally, the models were visualized using the PyMOL Molecular Graphics System (DeLano Scientific,
www.pymol.org
).
| Results |
|---|
|
|
|---|
Horses A2140, A2152, and A2150 all had the ELA-A1 haplotype as determined serologically, inheriting ELA-A1 from unrelated dams 172, 162, and 169, respectively (Table I).
|
The Rev-QW11 (QAEVLQERLEW) epitope is recognized by CTL from A2150 (27). Although EK target cells from all three horses were capable of presenting Rev-QW11 to Rev-QW11-specific A2150 CTL, these CTL recognized A2150 EK targets more efficiently than A2140 and A2152 targets (Fig. 1a). This observation suggested that the MHC I molecule presenting Rev-QW11 in A2150 was different from the one presenting Rev-QW11 in A2140 and A2152.
|
Because the 7-6 molecule that presents Env-RW12 and Gag-GW12 occurs in A2152 (51), and because Env-RW12- and Gag-GW12-specific CTL display similar recognition of A2152 and A2140 EK target cells (51), it was hypothesized that A2140 also possessed the 7-6 allele. Therefore, RT-PCR was used to amplify MHC I genes from A2140 PBMC. Cloning and sequencing confirmed the presence of the 7-6 allele in A2140. Of the 21 isolates processed, 5 were copies of a pseudogene, 3 were other pseudogenes, 6 were copies of a classical gene, and 7 were other classical genes, which included 7-6 (data not shown).
Because previous work indicates that A2150 EK target cells are recognized differently by Env-RW12- and Gag-GW12-specific CTL (51), and because Rev-QW11-specific CTL recognized A2150 EK target cells more efficiently than A2140 and A2152 targets, it was hypothesized that a MHC I molecule distinct from 7-6 presented these epitopes in horse A2150. A previous study identified partial sequences for three MHC I alleles in A2150, one of which, designated 141, shared the 7-6 sequence except for a single amino acid difference encoded at position 152 (E
V) in the
-2 domain (48). Due to its sequence similarity to 7-6, it was of interest to determine whether 141 had similar functional characteristics. To obtain the full-length 141 gene for expression, a cDNA library from A2150 was screened for MHC I genes and 141 was subsequently identified by sequencing (Fig. 1b).
Previous work indicated that the 141 allele is not present in A2152 and that the 7-6 allele is not present in A2150 (48). In addition, allele 141 was not identified among sequenced MHC I clones in A2140 (data not shown). Therefore, the 141 class I molecule was a likely candidate for presenting Env-RW12 and Rev-QW11 but not Gag-GW12 in horse A2150.
CTL recognized Env-RW12 more efficiently than Gag-GW12 when presented by the 7-6 molecule, and also recognized Rev-QW11 presented by 7-6
PBMC from A2140 were stimulated separately with Env-RW12 and Gag-GW12 peptides, then assayed for CTL activity on 7-6-transduced human 721.221 target cells that were pulsed with increasing concentrations of Env-RW12 and Gag-GW12 peptides, respectively. Human mutant B lymphoblastoid 721.221 cells express human β2m, but not HLA-A, -B, or -C class I molecules (58). Env-RW12 was recognized by CTL more efficiently than Gag-GW12, with 50% maximal lysis (EC50) of 2.8 nM for Env-RW12 vs 220 nM for Gag-GW12 (Fig. 1c).
Initial assays indicated that recognition of 7-6-transduced 721.221 cells pulsed with Rev-QW11 by A2150 Rev-QW11-specific CTL was equivocal (Fig. 2a). Because it was possible that the absence of equine β2m on 721.221 cells contributed to poor recognition of Rev-QW11 by CTL, ELA-A mismatched EK cells from pony H585 (ELA-A6) were transduced with 7-6, pulsed with Rev-QW11, and used as targets. Results indicated that in addition to Env-RW12 and Gag-GW12, 7-6 also presented Rev-QW11 (Fig. 2b).
|
A retroviral vector containing the 141 gene was constructed and used to transduce 721.221 cells and H585 EK cells. Rev-QW11-specific CTL from A2150 recognized both 141-transduced 721.221 and 141-transduced H585 EK target cells pulsed with the Rev-QW11 peptide (Fig. 2, c and d). Similarly, Env-RW12-specific CTL from A2140 recognized both 141-transduced 721.221 and 141-transduced H585 EK target cells pulsed with Env-RW12 (Fig. 2, c and d). In contrast, Gag-GW12-specific CTL from A2140 failed to recognize both 141-transduced 721.221 and 141-transduced H585 EK target cells pulsed with Gag-GW12 (Fig. 2, c and d).
CTL recognized Env-RW12 with similar efficiency when presented by the 7-6 and 141 molecules, whereas CTL recognized Rev-QW11 more efficiently when presented by 141
PBMC from horse A2140 were stimulated with Env-RW12, then assayed for CTL activity on 7-6- and 141-transduced 721.221 target cells that were pulsed with increasing concentrations of Env-RW12. The efficiency of Env-RW12 recognition on 7-6-transduced targets (EC50: 15 nM) was similar to that on 141-transduced targets (EC50: 8.8 nM) (Fig. 3a).
|
V amino acid substitution in class I molecule 141 increased the recognition efficiency of Rev-QW11 by CTL. Env-RW12 and Gag-GW12 bound to 7-6 and 141 with higher affinity than Rev-QW11
Live-cell Env-RW12 peptide-binding inhibition experiments were performed to determine the relative binding affinities of the three peptides to 7-6 and 141. The Env-RW12 peptide was chosen for 125I labeling because CTL recognized Env-RW12 on 7-6- and 141-expressing targets with similar efficiency (Fig. 3a), suggesting that Env-RW12 bound 7-6 and 141 with similar affinity. Moreover, Env-RW12 was the only peptide with a tyrosine residue, necessary for the chloramine-T method used for 125I labeling (68). The relative binding affinities of each peptide were then determined by using unlabeled peptides in competitive binding inhibition assays. Binding of 125I-labeled Env-RW12 to 7-6 was more efficiently inhibited by unlabeled Env-RW12 (IC50: 54 nM) than by Gag-GW12 (IC50: 71 nM) or Rev-QW11 (IC50: 525 nM) (Fig. 3c). For 7-6 binding, the IC50 of the negative control peptide 1b4a was >1000 nM (percent inhibition caused by 1000 nM was 27%; data not shown). These results were in agreement with the CTL recognition data.
Binding of 125I-labeled Env-RW12 to 141 was efficiently inhibited by unlabeled Env-RW12 (IC50: 41 nM) (Fig. 3d), consistent with the observation for 7-6. Surprisingly, binding of Env-RW12 to 141 was also inhibited by Gag-GW12 (IC50: 120 nM). Thus, the lack of Gag-GW12-specific CTL recognition of 141-expressing target cells was not due to the inability of Gag-GW12 to bind 141. However, the IC50 values indicated that Gag-GW12 bound 141 with lower affinity than 7-6 (120 nM vs 71 nM). Also unexpectedly, Rev-QW11 bound to 141 with lower affinity (IC50: 251 nM) than Gag-GW12. Consistent with the CTL results however, Rev-QW11 bound to 141 with higher affinity than it did to 7-6 (251 nM vs 525 nM). For 141 binding, the IC50 of the negative control peptide 1b4a was >1000 nM (percent inhibition caused by 1000 nM was 33%; data not shown). Taken together, these experiments suggested that differences in MHC/peptide-binding affinity were not sufficient to explain the differential CTL recognition of Gag-GW12 and Rev-QW11 on 7-6- and 141-expressing target cells.
Computer modeling and docking of peptides with the 7-6 and 141 molecules
Because crystal structures are not available, three-dimensional molecular modeling was performed to determine the possible structural and functional effects of the 152E
V substitution. Computer models were generated for 7-6 and 141 and the Env-RW12, Gag-GW12, and Rev-QW11 peptides. A docking algorithm was then used to dock each of the three peptides with 7-6 and 141, and docking scores for each complex were calculated. Modeling each of the MHC-peptide complexes indicated that amino acid position 152 occurred in the
-2 helix of 7-6 and 141, within the wall of the peptide-binding cleft, and it appeared that the conformations of the bound peptides were affected differently (Fig. 4). Modeling suggested that the peptides bound 7-6 and 141 in a bulged conformation, with the first and last residues of each peptide anchored in the binding clefts (Fig. 5). The conformation of Env-RW12 was similar when bound to 7-6 and 141 (Fig. 5a), whereas Rev-QW11 was slightly more bulged when bound to 7-6, presumably because of the W11 residue binding less deeply in the peptide-binding cleft of 7-6 (Fig. 5b). Interestingly, the bound conformation of Gag-GW12 was shifted and more sharply bulged in the 141 complex, apparently because the G1 residue of Gag-GW12 bound less deeply in the cleft of 141 as compared with 7-6 (Fig. 5c). Based on an analysis of pair potentials at the interface (76), the docking algorithm predicted the 7-6/Env-RW12 complex as the most favorable (highest RPScore docking score), followed by the 141/Rev-QW11 and 141/Env-RW12 complexes (Table II). The 7-6/Gag-GW12 and 7-6/Rev-QW11 complexes were less favorable and the 141/Gag-GW12 complex was the least favorable, although the docking score differences between these latter three complexes were not profound (Table II). Based on the MHC/peptide-docking models, the 7-6/Env-RW12 complex had a greater number of hydrophobic, salt bridge, hydrogen bond, and aromatic stacking interactions than did the 141/Env-RW12 complex (Table II and Fig. 6). For the 7-6/Rev-QW11 and 141/Rev-QW11 complexes, there were more hydrophobic, hydrogen bond, and aromatic stacking interactions for 7-6/Rev-QW11, but 141/Rev-QW11 had a greater number of salt bridges (Table II and Fig. 6). Importantly, the 7-6/Rev-QW11 complex had two destabilizing repulsive charged interactions that were absent in the 141/Rev-QW11 complex. For the 7-6/Gag-GW12 and 141/Gag-GW12 complexes, the two salt bridges in 7-6/Gag-GW12 were absent in 141/Gag-GW12, and 141/Gag-GW12 had one less hydrogen bond (Table II and Fig. 6). The salt bridges present in 7-6/Gag-GW12 but absent in 141/Gag-GW12 occurred between the 69E residue of 7-6 and the K4 residue of Gag-GW12 (Fig. 6, c and d, and Fig. 7, a and b). In addition, the G1 residue of Gag-GW12 had more interactions with residues of 7-6 (7Y: two hydrophobic, one hydrogen bond; 99Y: two hydrophobic; 159Y: three hydrophobic, two hydrogen bonds) than it did with residues of 141 (159Y: three hydrophobic) (Fig. 6, c and d, and Fig. 7, c and d). Based on the numbers of interactions between contact residues in the docking models for each of the six MHC/peptide complexes, the first two residues and W12 were probable anchor residues for Env-RW12 and Gag-GW12, whereas the first three residues and W11 were probable anchor residues for Rev-QW11 (Fig. 6). In general, the molecular modeling results supported the experimental data and suggested specific mechanisms for the observed differences in peptide binding and CTL recognition.
|
|
|
|
|
| Discussion |
|---|
|
|
|---|
V substitution in the
-2 domain of the 141 molecule enhanced CTL recognition of the Rev-QW11 epitope, but abolished CTL recognition of the Gag-GW12 epitope. Recognition of the Env-RW12 epitope by CTL was not affected. The hypothesis that the ELA-A1 haplotype was comprised of functionally distinct alleles was based on the initial observations that target cells from ELA-A1 horses A2140 and A2152 were recognized by both Env-RW12- and Gag-GW12-specific CTL whereas target cells from ELA-A1 horse A2150 were not recognized by Gag-GW12-specific CTL, and that Rev-QW11-specific CTL recognized A2150 target cells more efficiently than A2140 and A2152 target cells. Although the latter observation was consistent with the conclusion that the 7-6 MHC I molecule was used by A2140 and A2152 to present these epitopes while the 141 molecule was used by A2150, the efficiency of Rev-QW11-specific CTL recognition of A2140 and A2152 target cells was not exactly the same. The reason for the slightly greater recognition efficiency of the A2152 targets is not known. However, in addition to inheriting ELA-A1 haplotypes from different dams, A2140 and A2150 inherited different ELA haplotypes from the sire. It was therefore possible that the heterogeneous complement of other class I molecules on the respective target cells affected the expression of 7-6, or otherwise decreased peptide binding by 7-6 through competition or other unknown mechanisms (79, 80). Nonetheless, later experiments confirmed the hypothesis.
The loss of Gag-GW12 recognition by CTL was the most striking result of the single amino acid difference between the 7-6 and 141 molecules. Live-cell peptide-binding inhibition assays indicated that Gag-GW12 was able to bind 141, albeit with 1.7 times lower affinity than 7-6. Although this lower binding affinity could have contributed to the inability of Gag-GW12-specific CTL to lyse 141-expressing target cells, loss of TCR recognition of the 141/Gag-GW12 complex was probably more important. Others have shown that a single residue difference between two MHC class I molecules can dictate marked conformational differences in the same bound peptide, directly affecting TCR recognition by CTL (81). Although crystal structures are necessary for confirmation, molecular modeling suggested that the absence of salt bridges, along with the paucity of interactions between the G1 residue (a probable anchor residue) of Gag-GW12 and the 141 molecule, resulted in the observed lower-affinity binding, lower calculated docking score, and a more sharply protruding Gag-GW12 bound confirmation as compared with the 7-6/Gag-GW12 complex. These factors could have lowered TCR affinity for the 141/Gag-GW12 complex such that 141-restricted killing by Gag-GW12-specific CTL no longer occurred.
For the 7-6 molecule, Env-RW12 bound with only 1.3 times higher affinity than Gag-GW12, yet Env-RW12 CTL recognized 7-6 targets with 78 times greater efficiency than Gag-GW12 CTL. Therefore, the difference in affinity for the 7-6 MHC class I molecule between Env-RW12 and Gag-GW12 did not account for the difference in CTL-recognition efficiency, again suggesting that TCR affinity for the MHC/peptide complex played a major role in the discrepancy. Earlier work yielded similar results and conclusions (27, 51). Although the observed difference in 7-6 binding affinity for Env-RW12 and Gag-GW12 was small, molecular modeling provided some possible explanations for the difference. Specifically, the 7-6/Gag-GW12 complex had fewer hydrophobic interactions, fewer salt bridges, and a lower docking score than the 7-6/Env-RW12 complex.
Despite the detrimental effects on Gag-GW12-specific CTL recognition, the 152E
V change in 141 enhanced recognition of Rev-QW11 by CTL. Although CTL recognized Rev-QW11 on 7-6-transduced MHC I-mismatched EK cell targets, 7-6-transduced human 721.221 cells presented Rev-QW11 poorly. This type of differential recognition (equine vs human targets) was not observed for the other peptides, suggesting that human β2m did not effectively stabilize the 7-6/Rev-QW11 complex. This explanation is supported by the observation that human and murine β2m have different murine MHC I-stabilizing effects (82). It is not known why the presentation of the other peptides was not negatively affected by human β2m. Binding inhibition assays indicated that Rev-QW11 bound to 7-6 with 2.1 times lower affinity than to 141, and this inefficient binding could have made the stabilizing effects of autologous β2m more important for the 7-6/Rev-QW11 complex. Based on the molecular modeling results, the 7-6/Rev-QW11 complex had the second lowest docking score overall, and had fewer salt bridges and more destabilizing repulsive-charged interactions than did the 141/Rev-QW11 complex. In addition, and likely because of the reasons just listed, Rev-QW11 bulged out from the 7-6 binding cleft more than it did from 141. This conformational change in bound Rev-QW11 could have negatively affected TCR recognition of the 7-6/Rev-QW11 complex by Rev-QW11-specific CTL. Although crystal structures are needed, modeling provided plausible mechanisms for the inefficient binding of Rev-QW11 to 7-6, and for the inefficient Rev-QW11-specific CTL recognition of 7-6-expressing target cells.
Given the efficient Rev-QW11-specific CTL recognition of 141-expressing target cells and the absence of Gag-GW12-specific CTL recognition of 141-expressing target cells, the observation that Rev-QW11 bound 141 with half the affinity of Gag-GW12 was quite unexpected. Despite the lower MHC binding affinity, the 141-bound conformation of Rev-QW11 must have been efficiently recognized by the TCR of Rev-QW11-specific CTL. In contrast, the TCR of Gag-GW12-specific CTL probably could not bind the sharply protruding and highly bulged conformation of Gag-GW12 in the 141/Gag-GW12 complex.
The CTL epitopes (Env-RW12, Gag-GW12, and Rev-QW11) evaluated in this study were similar in that all three have a large aromatic tryptophan at the C terminus. Additionally, both the N-terminal and C-terminal residues are required for MHC binding and/or CTL recognition for all three epitopes (27, 51). For Env-RW12, the V2 residue is also a probable anchor residue based on 7-6 binding inhibition assays using peptides with amino acid substitutions (51). These observations are consistent with the molecular modeling results obtained in the present study. Analysis of the hydrophilic and hydrophobic amino acid residues for each of the three peptides indicate that Env-RW12 and Rev-QW11 differ from Gag-GW12 at positions 1, 2, 8, and 9. Both Env-RW12 and Rev-QW11 have hydrophilic residues at positions 1 and 8, whereas Gag-GW12 has hydrophobic residues at these positions. At positions 2 and 9, both Env-RW12 and Rev-QW11 have hydrophobic residues, whereas Gag-GW12 has hydrophilic residues at these positions. The differences in these residues, which included probable N-terminal anchor residues, likely contributed to the differences in MHC/peptide complex conformations that allowed CTL TCR recognition of all three peptides when presented by 7-6, but recognition of only Env-RW12 and Rev-QW11 when presented by 141.
Interestingly, all three peptides in this study were longer than the 8–10 aa generally considered optimal for MHC I binding. However, MHC I molecules can bind peptides as long as 14 aa in a bulged conformation and elicit dominant CTL responses (83). Importantly, the BZLF1 protein of EBV includes three completely overlapping CTL epitopes of 9, 11, and 13 aa in length (84). Although all three peptides bind well to the HLA-B*3501 molecule, the CTL response in individuals with this allele is directed exclusively toward the 11-mer epitope (84). Of particular interest, individuals with the B*3503 allele, which differs from B*3501 by a single amino acid in the F pocket of the peptide-binding cleft, do not mount CTL responses to these peptides because they do not bind B*3503. However, individuals with B*3508, which differs from B*3501 by a single amino acid in the D pocket of the peptide-binding cleft, develop CTL responses to the 13-mer epitope (84). The crystal structures indicate that the 13-mer binds both B*3508 and B*3501 in a centrally bulged conformation with the N and C termini anchored in the A and F pockets of the peptide-binding cleft (73). The differential CTL response is due to a broader peptide-binding cleft in B*3508, since the narrower binding cleft of the B*3501-peptide complex interacts poorly with the dominant TCR (73). Crystal structures will be required to confirm that the peptides in the present study bind the 7-6 and 141 molecules in a bulged conformation, and whether or not similar mechanisms are involved in the differential recognition by Env-RW12-, Gag-GW12-, and Rev-QW11-specific CTL.
Although crystal structures of equine MHC I molecules are lacking, the sequences of the two equine class I alleles presented here are surprisingly similar to human and murine class I alleles, and many of the residues forming the binding pockets A–F (in human and murine class I molecules) are shared (85, 86, 87). This suggests that the structure of equine MHC I molecules is similar to that of the mouse and human. In the absence of crystallography, molecular modeling provided important insights into the differential recognition of 7-6- and 141-expressing targets by EIAV Gag-GW12- and Rev-QW11-specific CTL. For example, the charged 152E to hydrophobic 152V substitution in the 141 molecule likely affected stabilizing salt bridges for Gag-GW12 binding, as seen when the 13-mer BZLF1 EBV peptide binds to B*3508 and B*3501, which differ only at residue position 156 (charged 156R
hydrophobic 156L) (73). If the modeling is correct, the absence of these stabilizing salt bridges contributed to the conformational change leading to loss of TCR recognition of the 141/Gag-GW12 complex.
This study confirms that a single amino acid difference between naturally occurring MHC I molecules can result in the loss of Gag-specific CTL recognition and enhanced (or diminished) efficiency of Rev-specific CTL recognition. The CTL, peptide-binding, and molecular modeling observations in this study support and suggest molecular mechanisms for the observed differences in disease progression and CTL responses in HIV-1-infected individuals of the B*35-Px and B*35-PY MHC I genotypes, because these genotypes also only differ by as few as one amino acid in the peptide-binding domain (40, 41). The implications of the observed differences in CTL recognition of important EIAV epitopes due to a single amino acid difference between otherwise identical MHC I molecules are important for designing protective lentivirus-specific CTL-inducing vaccines and understanding differential lentivirus disease progression in individuals within a population.
| Acknowledgments |
|---|
| Disclosures |
|---|
|
|
|---|
| Footnotes |
|---|
1 This work was supported in part by U.S. Public Health Service, National Institutes of Health Grants AI058787 (to R.H.M. and T.C.M.), AI067125 (to R.H.M. and T.C.M.), AI060395 (to T.C.M. and R.H.M.), CA97936 (to J.L.), U.S. Department of Agriculture National Research Initiative Grant 2002-35204-12699 (to J.L.), and the Center for Integrated Animal Genomics at Iowa State University (to J.L.). ![]()
2 Address correspondence and reprint requests to Dr. Robert H. Mealey, Department of Veterinary Microbiology and Pathology, Washington State University, Pullman, WA 99164-7040. E-mail address: rhm{at}vetmed.wsu.edu ![]()
3 Abbreviations used in this paper: EIAV, equine infectious anemia virus; MHC I, MHC class I; β2m, β2-microglobulin; ELA, equine leukocyte Ag; EK, equine kidney. ![]()
Received for publication January 17, 2006. Accepted for publication August 30, 2006.
| References |
|---|
|
|
|---|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |