Abstract
The primary function of the monomorphic MHC class Ib molecule Qa-1b is to present peptides derived from the leader sequences of other MHC class I molecules for recognition by the CD94-NKG2 receptors expressed by NK and T cells. Whereas the mode of peptide presentation by its ortholog HLA-E, and subsequent recognition by CD94-NKG2A, is known, the molecular basis of Qa-1b function is unclear. We have assessed the interaction between Qa-1b and CD94-NKG2A and shown that they interact with an affinity of 17 μM. Furthermore, we have determined the structure of Qa-1b bound to the leader sequence peptide, Qdm (AMAPRTLLL), to a resolution of 1.9 Å and compared it with that of HLA-E. The crystal structure provided a basis for understanding the restricted peptide repertoire of Qa-1b. Whereas the Qa-1b-AMAPRTLLL complex was similar to that of HLA-E, significant sequence and structural differences were observed between the respective Ag-binding clefts. However, the conformation of the Qdm peptide bound by Qa-1b was very similar to that of peptide bound to HLA-E. Although a number of conserved innate receptors can recognize heterologous ligands from other species, the structural differences between Qa-1b and HLA-E manifested in CD94-NKG2A ligand recognition being species specific despite similarities in peptide sequence and conformation. Collectively, our data illustrate the structural homology between Qa-1b and HLA-E and provide a structural basis for understanding peptide repertoire selection and the specificity of the interaction of Qa-1b with CD94-NKG2 receptors.
The MHC class I molecules play a central role in both adaptive and innate immunity. The highly polymorphic class Ia (or classical) MHC molecules are widely expressed in most tissues, where they function to present antigenic peptides to the αβ TCR expressed on the surface of T cells. The polymorphic sites within MHC-Ia molecules are generally focused within the Ag-binding cleft, allowing for the presentation of a large number of antigenic peptides (1). The selection of peptides that bind within the MHC-Ia cleft is usually determined by two pockets that specify primary anchor sites, with one or more secondary anchor sites that fine-tune the binding motifs for individual MHC-Ia allotypes. For example, HLA-B44 has a preference for P2-Glu and hydrophobic residues at PΩ (2), whereas HLA-B8 exhibits a preference for arginine residues at P3 and P5 of the peptide (3, 4). Thus, MHC polymorphism enhances immune defense throughout the population by ensuring a wide variation in the pathogen-derived peptides that are presented for T cell recognition (5).
In contrast to MHC-Ia molecules, MHC-Ib molecules, which include HLA-E, HLA-F, and HLA-G in humans and H2-M3 and Qa-1b in mice, are much less polymorphic and often exhibit a relatively restricted tissue distribution (6, 7). This limited polymorphism results in a marked reduction in the diversity of peptides presented by MHC-Ib proteins across the species. The capacity to present a diverse array of peptides appears further restricted by subtle changes to the peptide-binding groove that impose more stringent requirements for peptide binding. For example, H2-M3 preferentially binds peptides that have an N-formylated Met at P1, a modification that is largely restricted to peptides derived from mitochondrial and bacterial proteins (8). Other MHC-Ib molecules such as HLA-E and HLA-G use additional residues in the peptide to act as MHC-I anchor residues (9–12). HLA-E possesses five anchor sites that engender it being ideally suited to bind peptides derived from the conserved region of the leader sequence of other MHC class I molecules (e.g., amino acid sequence VMAPRTLFL from residues 3–11 of HLA-G). Similarly, HLA-G also possesses five anchor sites, and, although binding a broader repertoire of peptides in comparison with HLA-E, its peptide repertoire appears limited when compared with MHC-Ia molecules (13).
In humans, MHC-Ib molecules play a more prominent role in innate immunity, although they can also stimulate adaptive immune responses (6, 7, 14). For example, HLA-G has been reported to regulate NK cell activation (15, 16), whereas HLA-E also regulates the activation of NK cells by acting as a ligand for the CD94-NKG2 receptors (17, 18). The interaction between HLA-E and CD94-NKG2A represents a central innate mechanism by which NK cells indirectly monitor the expression of other MHC class I molecules within a cell. Recently, the structure of the CD94-NKG2A-HLA-EVMAPRTLFL complex provided a basis for understanding the specificity of this interaction as well as the relative contribution of the CD94 and NKG2A subunits in recognizing the HLA-EVMAPRTLFL complex (19, 20). In addition, HLA-E can also present peptides from pathogens such as CMV and Mycobacterium tuberculosis and stimulate cytotoxic T cell responses (21, 22).
Whereas the mode of HLA-E–mediated Ag presentation and recognition by innate and adaptive immune receptors is well established, much less is known regarding the murine Qa-1b molecule. Early studies suggested that Qa-1b is widely expressed (23) and involved in both innate and adaptive immunity stimulating both αβ (24, 25) and γδ (26) T cells. Recognition of Qa-1b by alloreactive CTL was dependent on the presence of a H-2D–linked gene referred to as the Qa-1 determinant modifier (Qdm) (23). Qdm was subsequently identified as a nonamer peptide (AMAPRTLLL) derived from residues 3–11 of the leader peptide present in most H-2D proteins (27). Peptide elution studies subsequently demonstrated that Qdm was the dominant peptide associated with Qa-1b (28).
Like HLA-E, Qa-1b is recognized by NK cells, specifically by CD94-NKG2A (29), in addition to the activating receptors CD94-NKG2C and CD94-NKG2E (30). Whereas initial studies suggested that HLA-E and Qa-1b are the result of convergent evolution (31), more recent analysis has found that because individual domains of the MHC undergo different rates of evolutionary change, HLA-E and Qa-1b are orthologs (32). As Qa-1b acts in an analogous role to HLA-E, acting as a checkpoint for class I expression by presenting leader peptides from other MHC-I molecules to CD94-NKG2A (18, 29), the evolutionary relationship between these two MHC-Ib molecules is of interest. Our findings explore the structure and function of Qa-1b and compare them with that of HLA-E.
Materials and Methods
Cloning and expression of Qa-1b and murine CD94-NKG2A
RNA from the cell line RMA-S was reverse transcribed and used to amplify a cDNA fragment encoding the extracellular domain of Qa-1b (residues 1–280) that was then cloned into a pET-expression vector. The Qa-1b H chain and murine β2-microglobulin were expressed separately in Escherichia coli, refolded, and purified, essentially as described previously (33, 34). Following reverse transcription of RNA prepared from C57BL/6 splenocytes, cDNAs encoding residues 43–179 and 119–224 of CD94 and NKG2A, respectively, were generated by PCR and cloned into a pET-30–based vector.
The extracellular domains of murine CD94 and NKG2A were expressed individually in E. coli from which inclusion bodies were isolated. These were resuspended in 8 M urea, 0.5 mM Na-EDTA, and 1 mM DTT, and CD94-NKG2A was refolded by flash dilution in a solution containing 5 M urea, 100 mM Tris-HCl (pH 8.0), 400 mM l-arginine-HCl, 5 mM reduced glutathione, and 0.5 mM oxidized glutathione. The refolding solution was then dialyzed against 1 M urea and 10 mM Tris-HCl (pH 8.0), followed by dialysis against 10 mM Tris-HCl (pH 8.0). The refolded CD94-NKG2A heterodimer was purified by anion exchange and size exclusion chromatography, followed by an additional anion–exchange chromatography step. HLA-E and human CD94-NKG2A (hCD94-NKG2A) were prepared, as described previously (20).
Surface plasmon resonance
The interaction between soluble murine and human CD94-NKG2A with Qa-1b or HLA-E was analyzed by surface plasmon resonance (SPR) using a Bio-Rad ProteOn XPR36 instrument (Hercules, CA). All experiments were performed at 25°C in a buffer containing 10 mM HEPES (pH 7.4), 150 mM NaCl, and 0.05% Tween 20 (HBS-T). hCD94-NKG2A and murine CD94-NKG2A (mCD94-NKG2A) were diluted into 10 mM sodium acetate (pH 5) and immobilized on adjacent flow cells of a GLC Sensorchip (Bio-Rad) by amine coupling (∼300 response units of each). An adjacent flow cell to which CD94-NKG2A was not added was activated and quenched in the same manner and served as a control cell. Recombinant Qa-1b and HLA-E were purified using size exclusion chromatography within 24 h of analysis, and the concentration of purified protein was estimated by measuring the absorbance at 280 nm. Qa-1b or HLA-E was then serially diluted in HBS-T (0.03–100 μM) and injected simultaneously over the test and control surfaces at a flow rate of 30 μl/min, and measurements were taken in duplicate. Following subtraction of data from control flow cells, the equilibrium KD values were analyzed using ProteOn Manager software version 3.0.1 (Bio-Rad) and Prism (GraphPad).
Crystallization, data collection, structure determination, and refinement
Crystals of Qa-1b-AMAPRTLLL were grown at 4 mg/ml protein in 6% PEG 4000, 10 mM NiCl2, and 0.1 M sodium acetate (pH 6.3), using the hanging drop vapor diffusion technique at 20°C. The crystals were flash cooled at 100°K using 15% glycerol as a cryoprotectant. A 1.9 Å resolution data set was collected at GM/CA-CAT (Advanced Photon Source, Chicago, IL) with a MarMosaic 300 CCD. The data were processed and scaled with HKL2000 (35). The crystal structure was solved by molecular replacement, as implemented in PHASER (36). The search model was HLA-E (Brookhaven Protein Data Bank [PDB] accession number: 2ESV) with all sequence differences mutated to alanine and the peptide removed. The progress of refinement was monitored by the Rfree value (5% of the data) with neither a σ nor a low-resolution cutoff being applied to the data. The structure was refined with REFMAC (37) and simulated annealing implemented in PHENIX (38) interspersed with rounds of model building in COOT (39). The final structure was validated with MolProbity (40, 41) and has 97.4% of residues within the favored region of the Ramachandran plot with no outliers. For data collection and refinement statistics, see Table I. Coordinates and structure factors are deposited in the PDB under accession number 3VJ6.
Data collection statistics | |
Temperature (°K) | 100 |
X-ray source | GM/CA-CAT, APS |
Detector | MarMosaic 300 CCD |
Space group | P212121 |
Cell dimensions | 59.27, 72.83, 95.84 |
Resolution | 41.45–1.9 (1.97–1.90) |
Total no. of observations | 193,458 |
No. of unique observations | 33,622 |
Multiplicity | 5.8 |
Data completeness (%) | 99.8 (99.0) |
1/σI | 31.5 (3.2) |
Rmerge (%)a | 7.8 (58.6) |
Refinement statistics | |
Nonhydrogen atoms | |
Protein | 3,184 |
Water | 273 |
Nickel | 4 |
Resolution | 1.9 |
Rfactor (%)b | 19.5 |
Rfree (%)b | 24.9 |
r.m.s.d. from ideality | |
Bond lengths (Å) | 0.007 |
Bond angles (°) | 1.188 |
Dihedrals (°) | 16.27 |
Ramachandran plotc | |
Favored regions (%) | 97.4 |
Allowed regions (%) | 2.6 |
B-factors (Å2) | |
Average main chain | 28.91 |
Average side chain | 36.92 |
Average water | 39.64 |
Average nickel | 49.13 |
Data collection statistics | |
Temperature (°K) | 100 |
X-ray source | GM/CA-CAT, APS |
Detector | MarMosaic 300 CCD |
Space group | P212121 |
Cell dimensions | 59.27, 72.83, 95.84 |
Resolution | 41.45–1.9 (1.97–1.90) |
Total no. of observations | 193,458 |
No. of unique observations | 33,622 |
Multiplicity | 5.8 |
Data completeness (%) | 99.8 (99.0) |
1/σI | 31.5 (3.2) |
Rmerge (%)a | 7.8 (58.6) |
Refinement statistics | |
Nonhydrogen atoms | |
Protein | 3,184 |
Water | 273 |
Nickel | 4 |
Resolution | 1.9 |
Rfactor (%)b | 19.5 |
Rfree (%)b | 24.9 |
r.m.s.d. from ideality | |
Bond lengths (Å) | 0.007 |
Bond angles (°) | 1.188 |
Dihedrals (°) | 16.27 |
Ramachandran plotc | |
Favored regions (%) | 97.4 |
Allowed regions (%) | 2.6 |
B-factors (Å2) | |
Average main chain | 28.91 |
Average side chain | 36.92 |
Average water | 39.64 |
Average nickel | 49.13 |
Values in parentheses are for the highest resolution shell.
Rmerge = 100 × Σ |Ihkl − <Ihkl>|/ΣIhkl.
Rfactor = 100 × Σhkl| |Fo| − |Fc| |/Σhkl |Fo| for all data except 5%, which was used for the Rfree calculation.
Molprobity (41).
Results
Qa-1b-AMAPRTLLL binds CD94-NKG2A with low affinity
Although the interaction between HLA-E and hCD94-NKG2A has been well characterized (19, 20, 42, 43), there has been no measurement of the affinity of the interaction between Qa-1b and mCD94-NKG2 receptors. Consequently, cDNAs encoding the extracellular domains of murine CD94 and NKG2A were cloned and expressed in E. coli, and the proteins refolded to generate soluble recombinant mCD94-NKG2A. Similarly, the extracellular domain of Qa-1b was also expressed in E. coli and refolded with β2-microglobulin and the Qdm peptide, and the interaction between refolded Qa-1b and mCD94-NKG2A was assessed by SPR. mCD94-NKG2A bound Qa-1b-AMAPRTLLL with rapid binding kinetics that precluded the accurate determination of both on and off rates, similar to the interaction between hCD94-NKG2A and HLA-E (Fig. 1) (43). Equilibrium-binding analyses yielded a KD of ∼17 μM. Given the capacity of both Qa-1b and HLA-E to bind peptides derived from the leader sequences of MHC-I molecules and the extensive sequence conservation between the mCD94-NKG2A and hCD94-NKG2A receptors, the extent to which Qa-1b and HLA-E could interact with either mCD94-NKG2A or hCD94-NKG2A was also assessed by SPR. hCD94-NKG2A bound HLA-EVMAPRTLLL with an affinity of ∼5 μM, an affinity similar to previously published data (Fig. 1) (20). Comparable affinities were obtained with lower amounts of receptor immobilized (∼65 response units; data not shown). Strikingly, there was no detectable interaction between hCD94-NKG2A and Qa-1b-AMAPRTLLL or between mCD94-NKG2A and HLA-EVMAPRTLLL (Fig. 1B, 1D).
SPR analysis of mCD94-NKG2A and hCD94-NKG2A/ligand interactions. Recombinant murine or hCD94-NKG2A was immobilized on a GLC Sensorchip and recombinant Qa-1b or HLA-E (0.03–100 μM) injected over the surfaces. Sensorgrams show Qa-1b binding to mCD94-NKG2A (A) and hCD94-NKG2A binding to HLA-E (C), both in duplicate. Saturation plots show the interactions of mCD94-NKG2A (B) and hCD94-NKG2A (D) with Qa-1b and HLA-E. Data shown are representative of at least four independent experiments.
SPR analysis of mCD94-NKG2A and hCD94-NKG2A/ligand interactions. Recombinant murine or hCD94-NKG2A was immobilized on a GLC Sensorchip and recombinant Qa-1b or HLA-E (0.03–100 μM) injected over the surfaces. Sensorgrams show Qa-1b binding to mCD94-NKG2A (A) and hCD94-NKG2A binding to HLA-E (C), both in duplicate. Saturation plots show the interactions of mCD94-NKG2A (B) and hCD94-NKG2A (D) with Qa-1b and HLA-E. Data shown are representative of at least four independent experiments.
Structural overview of Qa-1b
To address the structural basis of peptide binding to Qa-1b and the strict species specificity of the interactions between hCD94-NKG2A and mCD94-NKG2A with their MHC-Ib ligands, the crystal structure of Qa-1b in complex with the Qdm (AMAPRTLLL) peptide was determined. Crystals of the refolded Qa-1b-AMAPRTLLL complex formed in the space group P212121 with unit cell dimensions of a = 59.3, b = 72.8, and c = 95.8 Å, and diffracted to a resolution of 1.9 Å. The asymmetric unit comprised one heterodimer, and the structure was solved by molecular replacement using HLA-E as the search probe (34). The Qa-1b-AMAPRTLLL structure was subsequently refined to Rfactor and Rfree values of 19.5 and 24.9%, respectively. The final model contains all residues between 1 and 277 of the H chain and residues 1–99 of β2-microglobulin with an extra N-terminal, vector-derived, methionine residue (Fig. 2A, see Table I). The resultant electron density of the Qdm peptide and peptide-binding residues was unambiguous, allowing for ready interpretation of the interactions between the peptide and the Qa-1b H chain (Fig. 2B).
Overall structure of Qa-1b–Qdm complex. A, Qa-1b is shown in cartoon representation, with the Qdm peptide represented as sticks. The H chain is colored blue, with the α1 domain colored dark blue, the α2 domain colored cyan, and the α3 domain colored gray. The β2-microglobulin is colored orange. B, A side view of the Qdm peptide in the peptide-binding groove orthogonal to A. The α2-helix has been rendered transparent for clarity. A 2Fo–Fc electron density map contoured at 1 σ is shown as a yellow mesh within 1.5 Å of the peptide. The position of the residues within the Qdm peptide is indicated.
Overall structure of Qa-1b–Qdm complex. A, Qa-1b is shown in cartoon representation, with the Qdm peptide represented as sticks. The H chain is colored blue, with the α1 domain colored dark blue, the α2 domain colored cyan, and the α3 domain colored gray. The β2-microglobulin is colored orange. B, A side view of the Qdm peptide in the peptide-binding groove orthogonal to A. The α2-helix has been rendered transparent for clarity. A 2Fo–Fc electron density map contoured at 1 σ is shown as a yellow mesh within 1.5 Å of the peptide. The position of the residues within the Qdm peptide is indicated.
The overall structure of Qa-1b is similar to that of other MHC-I peptide complexes (Fig. 2A). The H chain contains three domains, the α1 and α2 domains that form the peptide binding cleft and the α3 domain that forms the heterodimeric interface with β2-microglobulin. Sequence comparison of Qa-1b with the MHC-Ia molecules H-2Kb (44) and H-2Db (45) and with the murine MHC-Ib molecules Qa-2 (46) and H2-M3 (47) shows significant similarities with sequence identities of 79, 76, 77, and 68%, respectively. The root mean square deviations (r.m.s.d.) upon superposition of Qa-1b and H2-Kb, Qa-2, H2-M3, and H-2Db are 1.0, 1.1, 1.2, and 1.4 Å, respectively (over all main-chain atoms). The majority of the structural differences between Qa-1b and the other murine MHC-I involved side chains clustered around the peptide-binding cleft, highlighting that subtle changes in this region reflect differences in Ag presentation and recognition by receptors expressed primarily by NK cells and T cells.
Peptide-binding specificity
The Qdm peptide is bound in an extended conformation, tethered by an extensive network of hydrogen bonds (Fig. 3, Table II). Contacts between Qa-1b and Qdm include 14 direct hydrogen bonds distributed across the length of the peptide and 8 water-mediated hydrogen bonds centered on the 4 C-terminal peptide residues P6–P9 (Table II). All of the Qa-1b residues involved in hydrogen bonding to Qdm were also conserved across Qa-1a, Qa-1c, and Qa-1d, and thus it is anticipated the mode of peptide binding will be shared across Qa-1 allotypes. With the exception of the P4-Pro and the P6-Thr, the main chain of each peptide residue is stabilized by a direct hydrogen bond with Qa-1b. For example, the amide group of P2-Met forms a hydrogen bond, which is conserved in most MHC-Ia structures, to the terminal carboxyl oxygen of Glu63, whereas Lys66 Nη forms a hydrogen bond to the P2-Met carbonyl group. The N terminus of the peptide is anchored by hydrogen bonds between the P1-Ala and the side-chain hydroxyl groups of Tyr7 and Tyr171 that are conserved in MHC class I–peptide structures (Figs. 3, 4A) (48). Similarly, the C terminus of P9-Leu is anchored by a network of hydrogen bonds to the side-chain hydroxyl groups of Tyr84 and Ser143 and to the side-chain amino group of Lys146, as well as forming water-mediated hydrogen bonds to the side chains of Asn77 and Thr80 (Fig. 3, Table II).
The hydrogen-bond network between Qa-1b and the Qdm peptide. The peptide is shown as a schematic in ball-and-stick representation, with carbons colored gray, oxygens red, nitrogens blue, and sulfurs yellow; the amino acid constituents of the peptide are labeled. Water molecules are represented as blue circles. Qa-1b residues and atoms involved in direct hydrogen bonding are labeled and highlighted in green boxes, and those involved in water-mediated hydrogen bonds are highlighted in blue boxes. The hydrogen bonds are represented as dashed lines, with the length of the bond indicated within the dashed lines in Å. The relative solvent-accessible surface area of the peptide atoms is represented as scaled cyan circles (solvent-accessible surface areas of <8 Å2 are not shown).
The hydrogen-bond network between Qa-1b and the Qdm peptide. The peptide is shown as a schematic in ball-and-stick representation, with carbons colored gray, oxygens red, nitrogens blue, and sulfurs yellow; the amino acid constituents of the peptide are labeled. Water molecules are represented as blue circles. Qa-1b residues and atoms involved in direct hydrogen bonding are labeled and highlighted in green boxes, and those involved in water-mediated hydrogen bonds are highlighted in blue boxes. The hydrogen bonds are represented as dashed lines, with the length of the bond indicated within the dashed lines in Å. The relative solvent-accessible surface area of the peptide atoms is represented as scaled cyan circles (solvent-accessible surface areas of <8 Å2 are not shown).
Peptide/Qa-1b contact residues
Peptide Residue | Qa-1b Residue | Interaction |
Ala1 | ||
Ala | Tyr7, Glu63, Tyr159, Trp167, Tyr171 | van der Waals |
AlaN | Tyr7 OH | H-bond |
Tyr171 OH | H-bond | |
AlaO | Tyr159 OH | H-bond |
Met2 | ||
Met | Tyr7, Ile24, Met45, Glu63, Lys66, Met70, Tyr99, Tyr159 | van der Waals |
MetN | Glu63 Oε2 | H-bond |
MetO | Lys66 Nζ | H-bond |
Ala3 | ||
Ala | Met70, Tyr99, Tyr159 | van der Waals |
AlaN | Tyr99 OH | H-bond |
Pro4 | ||
Pro | Tyr159 | van der Waals |
Arg5 | ||
Arg | Trp97, Val150, Glu152, Gln155, Gln156 | van der Waals |
ArgNε | Glu152 Oε2 | H-bond |
ArgNH2 | Glu152 Oε1 | Salt bridge |
ArgO | Gln156 Nε2 | H-bond |
Thr6 | ||
Thr | Asn73, Phe74, Trp97, Glu116, Glu152 | van der Waals |
ThrOγ1 | Asn77 Nδ2 | Water-mediated H-bond |
ThrOγ1 | Glu116 Oε2 | Water-mediated H-bond |
Leu7 | ||
Leu | Asn77, Glu116, Ile124, Trp133, Glu152 | van der Waals |
LeuO | Asn77 Nδ2 | H-bond |
Glu116 Oε2 | Water-mediated H-bond | |
Leu8 | ||
Leu | Asn77, Lys146 | van der Waals |
LeuN | Ser147 Oγ | Water-mediated H-bond |
Glu152 Oε1 | Water-mediated H-bond | |
LeuO | Ser143 O | Water-mediated H-bond |
Ser147 Oγ | Water-mediated H-bond | |
Leu9 | ||
Leu | Asn77, Thr80, Leu81, Tyr84, Leu95, Glu116, Ser143, Lys146 | van der Waals |
LeuN | Asn77 Oδ1 | H-bond |
LeuO | Tyr84 OH | H-bond |
Ser143 Oγ | ||
LeuOXT | Lys146 Nζ | H-bond |
Asn77 Oδ1 | Water-mediated H-bond | |
Thr80 Oγ1 | Water-mediated H-bond |
Peptide Residue | Qa-1b Residue | Interaction |
Ala1 | ||
Ala | Tyr7, Glu63, Tyr159, Trp167, Tyr171 | van der Waals |
AlaN | Tyr7 OH | H-bond |
Tyr171 OH | H-bond | |
AlaO | Tyr159 OH | H-bond |
Met2 | ||
Met | Tyr7, Ile24, Met45, Glu63, Lys66, Met70, Tyr99, Tyr159 | van der Waals |
MetN | Glu63 Oε2 | H-bond |
MetO | Lys66 Nζ | H-bond |
Ala3 | ||
Ala | Met70, Tyr99, Tyr159 | van der Waals |
AlaN | Tyr99 OH | H-bond |
Pro4 | ||
Pro | Tyr159 | van der Waals |
Arg5 | ||
Arg | Trp97, Val150, Glu152, Gln155, Gln156 | van der Waals |
ArgNε | Glu152 Oε2 | H-bond |
ArgNH2 | Glu152 Oε1 | Salt bridge |
ArgO | Gln156 Nε2 | H-bond |
Thr6 | ||
Thr | Asn73, Phe74, Trp97, Glu116, Glu152 | van der Waals |
ThrOγ1 | Asn77 Nδ2 | Water-mediated H-bond |
ThrOγ1 | Glu116 Oε2 | Water-mediated H-bond |
Leu7 | ||
Leu | Asn77, Glu116, Ile124, Trp133, Glu152 | van der Waals |
LeuO | Asn77 Nδ2 | H-bond |
Glu116 Oε2 | Water-mediated H-bond | |
Leu8 | ||
Leu | Asn77, Lys146 | van der Waals |
LeuN | Ser147 Oγ | Water-mediated H-bond |
Glu152 Oε1 | Water-mediated H-bond | |
LeuO | Ser143 O | Water-mediated H-bond |
Ser147 Oγ | Water-mediated H-bond | |
Leu9 | ||
Leu | Asn77, Thr80, Leu81, Tyr84, Leu95, Glu116, Ser143, Lys146 | van der Waals |
LeuN | Asn77 Oδ1 | H-bond |
LeuO | Tyr84 OH | H-bond |
Ser143 Oγ | ||
LeuOXT | Lys146 Nζ | H-bond |
Asn77 Oδ1 | Water-mediated H-bond | |
Thr80 Oγ1 | Water-mediated H-bond |
Comparison of the contacts mediated by the binding pockets of Qa-1b and HLA-E. Comparison of the peptide conformation within the Qa-1b– and HLA-E–binding groove (PDB accession number: 3BZE). The Qdm peptide is shown in stick representation, with carbons colored green, oxygens red, nitrogen blue, and sulfur yellow. The peptide bound to HLA-E is similarly represented, with carbons colored magenta. The position of the amino acid constituents of the peptides is labeled. The residues of Qa-1b and HLA-E are represented as sticks. Qa-1b carbons are colored in cyan, with HLA-E carbon atoms colored in pink. The residue number and type are labeled. A, View of the P1 amino acid in the binding cleft. B, The binding pocket of the P2 position. C, The binding pocket of the P3 position. D, The binding pocket of the P6 position and the salt-bridge interaction between the P5-Arg and Glu152. E, The binding pocket of the P7 position. F, The binding pocket of the P9 position.
Comparison of the contacts mediated by the binding pockets of Qa-1b and HLA-E. Comparison of the peptide conformation within the Qa-1b– and HLA-E–binding groove (PDB accession number: 3BZE). The Qdm peptide is shown in stick representation, with carbons colored green, oxygens red, nitrogen blue, and sulfur yellow. The peptide bound to HLA-E is similarly represented, with carbons colored magenta. The position of the amino acid constituents of the peptides is labeled. The residues of Qa-1b and HLA-E are represented as sticks. Qa-1b carbons are colored in cyan, with HLA-E carbon atoms colored in pink. The residue number and type are labeled. A, View of the P1 amino acid in the binding cleft. B, The binding pocket of the P2 position. C, The binding pocket of the P3 position. D, The binding pocket of the P6 position and the salt-bridge interaction between the P5-Arg and Glu152. E, The binding pocket of the P7 position. F, The binding pocket of the P9 position.
In addition to the hydrogen bond tethering of the termini of the Qdm peptide, Qa-1b also has five, primarily hydrophobic, binding pockets that also serve to anchor the peptide. These pockets accommodate the P2-Met, P3-Ala, P6-Thr, P7-Leu, and P9-Leu (Fig. 4B–F). The pockets accommodating P2, P7, and P9 are all relatively deep, whereas those in which P3 and P6 sit are significantly shallower (Fig. 5). Previous studies had established that P2-Met served as a dominant anchor residue for peptide binding to Qa-1b (49). This residue sits in the deep B pocket, making a series of van der Waals interactions within a hydrophobic cavity lined by Tyr7, Ile24, Met45, Met70, and Tyr99 (Fig. 4B). The residues Met70, Tyr99, and Tyr159 from the B pocket also combine with Trp97 to shape the shallow P3 pocket (Fig. 4C). Similarly, the P6-Thr also sits in a shallow hydrophobic pocket formed by the side chains of Phe74 and Trp97 and the aliphatic side chain of Asn73 (Fig. 4D). The P7-Leu is located in a deep pocket accommodated by van der Waals interactions with Glu116, Ile124, Trp133, and Glu152 (Fig. 4E). Finally, the C-terminal peptide residue, P9-Leu, is tethered by a series of main-chain hydrogen bonds as well as being anchored via side-chain interactions with Asn77, Thr80, Leu81, Tyr84, Leu95, Glu116, and Ser143 (Fig. 4F).
Comparison of the anchoring pockets in Qa-1b and HLA-E. Space-filling model showing peptide binding to Qa-1b (A) and HLA-E (B). The Qdm peptide in Qa-1b is colored green and the VMAPRTLFL peptide in HLA-E (PDB accession number: 3BZE) is colored magenta. The pockets of the MHCs are contoured about the side-chain atoms of the anchoring residues with a probe radius of 1.4 Å.
Comparison of the anchoring pockets in Qa-1b and HLA-E. Space-filling model showing peptide binding to Qa-1b (A) and HLA-E (B). The Qdm peptide in Qa-1b is colored green and the VMAPRTLFL peptide in HLA-E (PDB accession number: 3BZE) is colored magenta. The pockets of the MHCs are contoured about the side-chain atoms of the anchoring residues with a probe radius of 1.4 Å.
Substitution at the P4 and P8 positions of the Qdm peptide had little effect on binding to Qa-1b (49). From the structure, these positions were not confined to Qa-1b–binding pockets, but were solvent exposed (Figs. 2B, 5). Interestingly, the P4-Pro made no significant interactions with Qa-1b, suggesting it may play a role in constraining the conformation of the peptide, which is consistent with the importance of intramolecular constraints in MHC-Ia peptide presentation (50). Similarly, whereas the side chain of P8-Leu made no significant contacts with Qa-1b, its main chain bound to Qa-1b by a series of water-mediated hydrogen bonds. Interestingly, whereas the P5-Arg side chain was largely solvent exposed, substitution of this residue to Ala impairs peptide binding to Qa-1b (49). This effect appears attributable to the disruption of a salt bridge between P5-Arg and Glu152 on the α2 helix of Qa-1b (Fig. 4D). Collectively, the structure of the Qa-1b-AMAPRTLLL complex confirmed the role of the P2, P3, P6, P7, and P9 positions as anchor residues, which together act to impose considerable constraints on diversity of the peptide repertoire associated with Qa-1b. In contrast, the P4, P5, and P8 side chains are solvent exposed, thus potentially forming part of the ligand recognized by NK and T cells.
Qa-1b is a structural homolog of HLA-E
The H chains of Qa-1b and HLA-E share 74% overall identity, and 71% identity in the peptide-binding region (α1 and α2 domains) (Figs. 6, 7). The overall architecture of this region is structurally well conserved with r.m.s.d. of 0.90 Å over main-chain atoms (residues 2–182), a point reflected in the ability of Qa-1b to bind human MHC-I leader sequences presented by HLA-E (51) (Fig. 6). Similarly, this structural conservation is reflected in the conformation of the Qdm peptide, which shows remarkable similarity to the MHC-I–derived peptides presented by HLA-E with backbone r.m.s.d. of 0. 6 Å over main-chain atoms (Fig. 6).
Overview of the structural conservation of the Qa-1b and HLA peptide-binding groove. A, Structural superposition of the α1 and α2 domains of Qa-1b (blue) with the three HLA-E structures, as follows: PDB accession number 3BZE (10) (red), 3BZF (10) (yellow), and 1KTL (52) (gray), demonstrating the highly conserved structure of the Ag-binding cleft. B, Conformation of peptides bound to Qa-1b and HLA-E, as follows: Qdm (green), VMAPRTLFL (pink), VMAPRTLLL (yellow), and VTAPRTLLL (gray). The α1 and α2 domains were superposed, as in A.
Overview of the structural conservation of the Qa-1b and HLA peptide-binding groove. A, Structural superposition of the α1 and α2 domains of Qa-1b (blue) with the three HLA-E structures, as follows: PDB accession number 3BZE (10) (red), 3BZF (10) (yellow), and 1KTL (52) (gray), demonstrating the highly conserved structure of the Ag-binding cleft. B, Conformation of peptides bound to Qa-1b and HLA-E, as follows: Qdm (green), VMAPRTLFL (pink), VMAPRTLLL (yellow), and VTAPRTLLL (gray). The α1 and α2 domains were superposed, as in A.
Structure-based sequence conservation of Qa-1b with HLA-E. A, Surface representation of Qa-1b–Qdm. Qdm is shown in yellow. The surface of Qa-1b is colored according to the sequence conservation with HLA-E shown in Fig. 4B. Conserved residues are colored green, and nonconserved residues are colored white. Residues of HLA-E that contact the CD94-NKG2A receptor are labeled. Residues contacting the NKG2A are labeled in black, and those that contact CD94 are labeled in pink. Trp65 is labeled in blue, as it contacts both CD94 and NKG2A. B, Structure-based sequence alignment of the α1 and α2 domains of Qa-1b and HLA-E. Strictly conserved residues are highlighted in green boxes.
Structure-based sequence conservation of Qa-1b with HLA-E. A, Surface representation of Qa-1b–Qdm. Qdm is shown in yellow. The surface of Qa-1b is colored according to the sequence conservation with HLA-E shown in Fig. 4B. Conserved residues are colored green, and nonconserved residues are colored white. Residues of HLA-E that contact the CD94-NKG2A receptor are labeled. Residues contacting the NKG2A are labeled in black, and those that contact CD94 are labeled in pink. Trp65 is labeled in blue, as it contacts both CD94 and NKG2A. B, Structure-based sequence alignment of the α1 and α2 domains of Qa-1b and HLA-E. Strictly conserved residues are highlighted in green boxes.
The network of charged interactions between the H chain and the peptide is well conserved between Qa-1b and HLA-E. For example, P5-Arg forms a salt bridge to Glu152 of Qa-1b, a feature also observed in HLA-E/peptide interactions (Fig. 4D) (11). Likewise, with the exception of positions 66 and 99, all residues involved in direct hydrogen bonds with the peptide main chain are conserved (Fig. 4). Lys66 in Qa-1b forms a hydrogen bond with the backbone carbonyl of the P2-Met (Fig. 4B). In HLA-E, position 66 is a Ser and also forms a stabilizing contact to P2, albeit through a bridging water molecule that contacts Glu63. By contrast, the hydrogen bond from Tyr99 of Qa-1b to the peptidyl amide of P3 is not maintained by the His99 substitution in HLA-E (10) (Fig. 4C). Conceivably, this may allow greater flexibility in the positioning of the main chain at P3, which may also impact on the positioning of the P4 residue. Indeed, variation between HLA-E structures is found in the P3 and P4 positions (10, 11, 52) (Fig. 6B). The HLA-EVTAPRTLLL and HLA-EVMAPRTVLL structures adopt a similar conformation at P3/P4, which differs by ∼0.7 Å and ∼1.2 Å, at P3 and P4, respectively, compared with the conformation in the HLA-EVMAPRALLL and HLA-EVMAPRTLFL structures (10, 11, 52). This variation is coincident with a shift of His99Cε1 of ∼1.1 Å.
The dimensions of the anchor residue-binding pockets are well conserved between Qa-1b and HLA-E structures, with subtle differences underlying small variations in peptide conformation. The orientation of residues Tyr7, Met45, Glu63, and Tyr159 that line the deep hydrophobic cavity of the P2-binding pocket is similar between both Qa-1b and HLA-E (Fig. 4B). However, substitution at position 24 of Ile in Qa-1b for the polar Ser in HLA-E and at position 70 of the hydrophobic Met for the polar Thr results in a shallower cavity in Qa-1b and an altered conformation of the P2-Met side chain. The deeper cavity of HLA-E allows the P2-Met to adopt two conformations, as observed in the crystal structures determined to date (10) (3BZE versus 3BZF). The steric hindrance imposed by Met70 in Qa-1b prevents the P2-Met from adopting either of the conformations observed in HLA-E. Instead, the side chain is shifted 1.1 Å toward the face of the pocket lined by residues Tyr7, Ile24, and Met45 (Figs. 4B, 5). Furthermore, the substitution of Tyr99 and Lys66 in Qa-1b for His and Ser, respectively, in HLA-E not only affects hydrogen-bonding interactions with P2 and P3, but also narrows the binding groove about the Cα atom of P2. These substitutions are coincident with the P2-Met being shifted 0.8 Å toward the α1 helix relative to its HLA-E counterpart (Fig. 4B).
Underscoring the strict conservation of Leu at P9 in peptides presented by both Qa-1b and HLA-E, residues Asn77, Thr80, Leu81, Leu95, and Ser147 that line the P9 pocket are conserved (Figs. 4F, 7). Likewise, the orientation of the P9-Leu is structurally conserved between all Qa-1b and HLA-E structures (Fig. 6B). The P7 pocket, by contrast, has subtle structural differences that lead to a comparatively deeper pocket in Qa-1b. The substitution of Glu116 and Ile124 in Qa-1b to Phe116 and Leu124 in HLA-E combined with the rotation of Phe74 opens the P6- and P7-binding pockets. As a result, the P6-Thr and P7-Leu sit 0.8 Å and 1.8 Å deeper in their respective pockets in Qa-1b (Figs. 4E, 6B). Conceivably, the deeper P7 pocket in Qa-1b would allow for greater structural flexibility at this position and the ability to accommodate bulkier amino acids. The deeper positioning of P7 corresponds also with a shift of Trp133 by 1.5 Å. The combined shifts of P7 and Trp133 position the carboxylate of Glu152 2.7 Å closer to the floor of the peptide-binding groove (Fig. 4D, 4E), which is concomitant with a small shift of the guanidinium group of P5-Arg (Fig. 4D).
Structural and mutagenesis studies have shown that CD94-NKG2A recognition of HLA-E is dependent on a cluster of basic residues on the α1 helix, including residues Arg65, Arg75, and Arg79, a hydrophobic patch that includes the P8 side chain that abuts this region together with an acidic region on the α2 helix spanning residues 154–162 (19, 20, 53). Comparison of the corresponding regions in Qa-1b showed a similar charge distribution, with an electropositive patch on the α1 domain and an electronegative patch on the α2 domain (Fig. 8), both of which are also conserved across Qa-1b allotypes. Despite these broad similarities, however, there are numerous differences between Qa-1b and the corresponding residues in HLA-E that have been shown to contact CD94-NKG2A (19), in particular residues 154 and 155, which directly interact with human NKG2A, and residues 72 and 73, the latter of which interacts with P6-Thr, a residue shown to be important for the interaction with hCD94-NKG2A (Fig. 7) (54). Similarly, mutation of Arg72 of Qa-1b has also been shown to disrupt recognition by mCD94-NKG2A (55).
Comparison of the electrostatic potential on the surface of Qa-1b and HLA-E. Electronegative and electropositive patches shown to complement binding regions on CD94 and NKG2A, respectively, are labeled. Likewise, a patch of nonconserved residues shown to be important for the interaction with CD94 is indicated. Although the overall charge distribution between the two proteins is similar, it is likely that differences in electrostatic interactions determine differences in CD94-NKG2A recognition. The potentials are displayed as surfaces contoured at levels of −5 kT/e (red) and +5 kT/e (blue).
Comparison of the electrostatic potential on the surface of Qa-1b and HLA-E. Electronegative and electropositive patches shown to complement binding regions on CD94 and NKG2A, respectively, are labeled. Likewise, a patch of nonconserved residues shown to be important for the interaction with CD94 is indicated. Although the overall charge distribution between the two proteins is similar, it is likely that differences in electrostatic interactions determine differences in CD94-NKG2A recognition. The potentials are displayed as surfaces contoured at levels of −5 kT/e (red) and +5 kT/e (blue).
Discussion
Qa-1b has a number of distinct roles in the immune system serving both to regulate innate immunity via interactions with members of the CD94-NKG2 receptor family and stimulate adaptive immunity presenting peptides for recognition by the αβ TCR (29, 56). The data presented in this work identify the structural features of Qa-1b that allow for the selective binding of the Qdm peptide. The key features of Qdm involved in Qa-1b binding, notably P2-Met, P7-Leu, and P9-Leu, are highly conserved across mammalian MHC class I sequences, consistent with observations that similar peptides derived from other mammalian class I leader sequences also bind Qa-1b (51).
Peptide elution studies suggested that Qa-1b imposes tight specificity on the bound peptide (28). Furthermore, binding of Qdm to Qa-1b was greatly reduced by substitution of P2-Met and P9-Leu with Ala, and the double substitution at P2 and P9 completely abrogated peptide binding (49). These observations are consistent with the structure of Qa-1b-AMAPRTLLL, which revealed deep pockets that could accommodate the P2-Met and P9-Leu side chains of the Qdm determinant. However, the structural data presented also show that residues lining the shallow binding pockets that accommodate the P3 and P6 side chains make critical contacts with Qdm, consistent with the observation that a poly-Gly peptide with P2-Met and P9-Leu did not bind to Qa-1b (51).
Whereas the Qdm peptide is the dominant peptide eluted from Qa-1b purified from healthy cells, there is evidence to suggest that Qa-1b can bind a diverse array of peptides (56, 57). In particular, Qa-1b isolated from TAP-deficient cells was shown to have a markedly broader repertoire of peptides than originally reported (56). These studies showed little evidence of enrichment of any amino acid at P2, suggesting that a variety of amino acids can be accommodated at this position. Similarly, whereas P9 of both the Qdm sequence and its human counterparts (residue 11) is a Leu, the peptides isolated from Qa-1b obtained from TAP-deficient cells were enriched for additional aliphatic amino acids, including Ile, Ala, Met, and Phe at this position. Indeed, a variety of defined peptides that are distinct from MHC-I leader sequences has been shown to bind to Qa-1b. These include peptides from Hsp60 (GMKFDRGYI) (58), Salmonella (GMQFDRGYL) (58), insulin (ALWMRFLPL) (59), and influenza (FLARSALIL) (57). The structural data in this study suggest that it is improbable that the Hsp60- and Salmonella-derived peptides could be accommodated in a conformation similar to Qdm. First, the introduction of long side chains Lys and Gln at the shallow P3 pocket and Arg at the shallow P6 pocket would require conformational changes to the peptide-binding groove. Furthermore, the loss of the P5 salt bridge to Glu152 and the loss of the stabilizing Leu in the deep P7 pocket by the substitution to Gly would impact on peptide stability and conformation. Likewise, modeling of the insulin-derived peptide suggests significant conformational adjustment of the shallow P3 and P6 pockets to accommodate the bulky Trp and Phe side chains. The influenza-derived peptide appears closest to the Qdm peptide in conforming to the needs of Met/Leu at P2 and Leu at P7 and P9 with small side-chain residues at P3 and P6. However, the introduction of Ser at P5 would alter the salt-bridge interaction with Glu152, potentially impacting on the receptor recognition surface of the molecule. Although these peptides have all been shown to stimulate a T cell response, their capacity to interact with CD94-NKG2A is unclear. Modeling based on the Qa-1b-AMAPRTLLL complex suggests that with the possible exception of the influenza peptide, these peptides are unlikely to adopt a conformation similar to that of Qdm.
There are a number of crystal structures of different MHC molecules binding the one peptide (2, 60, 61). Typically, these involve highly related allotypes that exhibit limited sequence differences, which nevertheless often result in biologically significant changes in conformation of the peptide and/or the MHC-I molecule (2, 61–63). In contrast, despite the comparatively low level of sequence identity between Qa-1b and HLA-E, remarkably, they present their respective peptide cargoes in a very similar manner for recognition by CD94-NKG2 receptors. Whereas previous studies have shown that Qdm peptide can bind HLA-E and be recognized by hCD94-NKG2A (42), we show that the Qa-1b-AMAPRTLLL complex does not bind hCD94-NKG2A consistent with previous data showing that HLA-E and Qa-1b tetramers fail to bind murine splenocytes or human PBMC, respectively (64). Thus, we suggest that the inability of hCD94-NKG2A to recognize Qa-1b-AMAPRTLLL most likely results from differences between the Qa-1b and HLA-E H chain rather than changes in orientation and sequence of the peptide per se.
Finally, a notable feature of Qa-1b and indeed its human counterpart, HLA-E, is the relative depth of the C and E pockets that accommodate P6 and P7, respectively. Consequently, the peptide sits deeper in the groove compared with typical pMHC-Ia complexes. Given that the primary function of MHC-Ia molecules is to display a diverse group of peptides for T cell recognition, the elevated positioning of the peptide backbone in MHC-Ia molecules may result in a greater range of ligands for TCR recognition. In contrast, the deep position of the peptide backbone in both Qa-1b and HLA-E, which is also a feature of peptide binding to a number of other nonclassical MHC class I molecules, including HLA-G (9), might reflect their primary role is the presentation of a conserved group of peptides for recognition by essentially invariant receptors that regulate innate immune responses.
Footnotes
This work was supported by grants from the Australian Research Council and the National Health and Medical Research Council; an Australian postgraduate scholarship (to L.Z.); a Career Development award from the National Health and Medical Research Council (to L.C.S.); a Doherty fellowship from the National Health and Medical Research Council (to J.P.V.); a Federation fellowship from the Australian Research Council (to J.R.); and a Queen Elizabeth II fellowship from the Australian Research Council (to C.S.C.).
The coordinates presented in this article have been submitted to the Protein Data Bank (http://www.rcsb.org/pdb/home/home.do) under accession number 3VJ6.
References
Disclosures
The authors have no financial conflicts of interest.