Previous Article | Next Article ![]()
Journal of Virology, February 2008, p. 1238-1248, Vol. 82, No. 3
0022-538X/08/$08.00+0 doi:10.1128/JVI.02026-07
Copyright © 2008, American Society for Microbiology. All Rights Reserved.

Department of Biochemistry, Tulane University Health Sciences Center, 1430 Tulane Avenue, New Orleans, Louisiana 70112
Received 13 September 2007/ Accepted 15 November 2007
|
|
|---|
|
|
|---|
Epitope dominance can be based on the frequency of response in a population of individuals or the intensity of response in a particular individual, although dominant epitopes often satisfy both definitions (15, 21). Peptide affinities for class II major histocompatibility (MHC) antigen-presenting proteins have not explained CD4+ epitope immunodominance as well as peptide affinities for class I MHC proteins have explained CD8+ epitope immunodominance (76). The class II MHC protein has been distinguished from the class I MHC protein in that the class II MHC protein is exquisitely adapted for nonspecific binding to peptides (38). The presence or absence of certain class II MHC alleles influences the dominance of particular epitopes (1, 43). However, peptide affinity for the class II MHC protein is a poor predictor of epitope dominance in influenza virus hemagglutinin (HA) (22) and other antigens (40, 46, 52). The weak influence of MHC-peptide affinity on epitope dominance is further indicated by the frequent observation of "promiscuous" or "universal" CD4+ epitopes, which are dominant in multiple MHC backgrounds (17, 20, 28, 39, 50, 54, 60). Thus, other factors exert a major influence on epitope immunogenicity.
An epitope's structural context has been shown to influence CD4+ T-cell epitope dominance at the level of antigen processing (4, 9, 14, 45, 62, 65, 67, 68, 73), but a general relationship between structure and immunodominance has not been described. Processing of CD4+ epitopes occurs in a lysosome-like compartment, where there is no evidence for ATP-dependent protein-unfolding mechanisms other than lysosomal acidification (72). Protease-sensitive sites in proteins have been correlated with conformational flexibility, as reported by high crystallographic B-factors (8, 19, 36, 49), and dominant CD4+ epitopes have been shown to occur adjacent to flexible antigen segments identified by crystallographic B-factors in several model antigens (42), tetanus toxoid (17), Hsp10s (8), and HIV gp120 (15).
In the present study, published epitope maps are compared to profiles of conformational stability in influenza virus HA and eight other well-studied antigens/allergens. The profiles of conformational stability were generated with the COREX algorithm, which evaluates the dominant forces of protein folding, namely, the hydrophobic interaction and conformational entropy (33). The analysis yields information about probable pathways of proteolytic antigen processing and indicates that dominant CD4+ epitopes arise on the flanks of stable antigen segments, which most likely correspond to proteolytic processing intermediates. The relationship of epitope dominance to antigen structure focuses immunity onto conserved antigen sequences but also creates a mechanism for immune evasion by pathogens such as influenza virus.
|
|
|---|
Identification of dominance peaks. (i) HA. Proliferative responses to HA in 30 human subjects reported by Gelder et al. (22-24) were converted to stimulation indices (SI) by dividing the response (cpm) for each peptide pool by the response for the unstimulated control. The response was scored as positive for a given subject if the SI was greater than 4. Epitope scores (expressed as the percentage of subjects responding) were assigned to residues according to the frequency of response to the peptide pool that contained the residue. Each pool of five 15-mer peptides spanned 35 residues, and the 24 pools spanned the complete sequence of HA. The peptide pools partially overlapped, and therefore residues that appeared in two peptide pools were assigned the average of the two scores. Dominance peaks were identified as local maxima in the smoothed profile of epitope scores. Smoothing facilitated the identification of a single amino acid residue at the center of each dominance peak, which was characterized by a broad plateau in the raw profile. Smoothing did not affect the number of dominance peaks. The profile of raw epitope scores was smoothed with a 23-residue moving window average. A "first-derivative" profile was generated by taking the difference between scores for the current residue and the preceding residue, and then this profile was smoothed with a 13-residue moving window average. Dominance peaks were assigned to the residues where the smoothed first derivative changes from positive to negative. For Fig. 1, the responses of H3-infected C57BL/6 mice and of H1-infected HLA-DR1 trangenic mice were assigned to residues according to the average response to the peptides that contained the residue. In the study using H1-infected HLA-DR1 mice, responses had been scored according to the percentage of the maximum, defined by the number of enzyme-linked immunospots (ELISPOTs) obtained for peptide 63 (corresponding to 443 to 465 in the x31 strain) as follows: negative (less than 6%), weak (6 to 24%), moderate (25 to 60%), or strong (greater than 60%). In the present work, peptides were assigned scores corresponding to the average response for the range: negative (3%), weak (15%), moderate (43%), or strong (81%). Residues that appeared in multiple peptides were assigned the average of the scores.
![]() View larger version (28K): [in a new window] |
FIG. 1. Similar profiles of CD4+ epitope dominance for influenza virus HA in mice and humans infected with H3 or H1 strains of virus. (a) Schematic diagram of HA subdomains. F', HA1 portion of fusion domain; E, vestigial esterase; RB, receptor binding; F, HA2 portion of fusion domain. (b) Profile of epitope scores for a population of 30 human subjects. (c) Profile of epitope scores for a group of C57BL/6 mice that had been inflected with an H3N1 strain. IFN , gamma interferon. (d) Profile of epitope scores for a group of HLA-DR1 transgenic mice that had been infected with an H1N1 strain. IL-2, interleukin-2. (e) Profiles of sequence entropy calculated for HA molecules representing five HA subtypes ("Shift") or for 135 HA molecules representing different isolates of viral subtype H3N2 ("Drift"). Vertical lines indicate peaks of epitope dominance that are shared in all three studies (heavy lines) or shared between only two studies (light lines).
|
(iii) Mycobacterium tuberculosis Ag85A. The responses of 10 tuberculin-positive subjects reported by Launois et al. (44) were based on lymphocyte proliferation stimulated by 20-mer peptides spanning Ag85A with 10-residue overlaps. The authors scored the response for a given subject as positive if the SI was greater than 5. Epitope scores (expressed as the percent of subjects responding) were assigned to residues according to the average frequency of response to the peptides that contained the residue. The profile of raw epitope scores was smoothed with a 13-residue moving window average. The first-derivative profile generated as for HA was smoothed with a 13-residue moving window average. Dominance peaks were assigned as for HA.
(iv) M. tuberculosis Mpb70. The responses of 14 BCG-vaccinated subjects reported by Al-Attiyah et al. (2) were based on lymphocyte proliferation stimulated by 25-mer peptides spanning Mpb70 with 10-residue overlaps. The authors scored the response for a given subject as positive if the SI was greater than 3. Epitope scores (expressed as the percentage of subjects responding) were assigned to residues according to the frequency of response to the peptide that contained the residue. Residues that appeared in multiple peptides were assigned the average of the scores. The profile of raw epitope scores was smoothed with a 15-residue moving window average. The first-derivative profile generated as for HA was smoothed with a 15-residue moving window average. Dominance peaks were assigned as for HA.
(v) Mycobacterium bovis Ag85B. The responses of 15 BCG-vaccinated subjects reported by Roche et al. (57) were based on lymphocyte proliferation stimulated by 20-mer peptides spanning Ag85B with 10-residue overlaps. The authors scored the response for a given subject as positive if the SI was greater than 4. Epitope scores (expressed as the percentage of subjects responding) were assigned to residues according to the frequency of response to the peptide that contained the residue. Residues that appeared in multiple peptides were assigned the average of the scores. The profile of raw epitope scores was smoothed with a 17-residue moving window average. The first-derivative profile generated as for HA was smoothed with a 13-residue moving window average. Dominance peaks were assigned as for HA.
(vi) M. tuberculosis Hsp65. The responses of immunized Lewis rats reported by Moudgil et al. (47) were based on lymphocyte proliferation stimulated by 15-mer peptides spanning Hsp65 with 11-residue overlaps. Epitope scores (expressed as thousands of cpm) were assigned to residues according to the response to the peptide that contained the residue. Residues that appeared in multiple peptides were assigned the average of the scores. The profile of raw epitope scores was smoothed with a 15-residue moving window average. The first-derivative profile generated as for HA was smoothed with a five-residue moving window average. Dominance peaks were assigned as for HA.
(vii) Bet v 1. The responses of 22 T-cell clones from nine allergic patients reported by Dormann et al. (18) were based on proliferation stimulated by 12-mer peptides spanning Bet v 1 with 10-residue overlaps. The authors scored the response for a given clone as positive if the SI was greater than 5. Epitope scores were assigned to residues according to the frequency that T-cell clones responded to the peptide that contained the residue. The profile of raw epitope scores was smoothed with a seven-residue moving window average. The first-derivative profile generated as for HA was smoothed with a seven-residue moving window average. Dominance peaks were assigned as for HA.
(viii) Ves v 5. The responses of T-cell lines from 13 allergic patients reported by Bohle et al. (6) were based on proliferation stimulated by 12-mer peptides spanning Ves v 5 with nine-residue overlaps. The authors scored the response for a given clone as positive if the SI was greater than 4. Epitope scores were assigned to residues according to the number of patients whose lines responded to the peptide that contained the residue. The profile of raw epitope scores was smoothed with a 15-residue moving window average. The first-derivative profile generated as for HA was smoothed with a seven-residue moving window average. Dominance peaks were assigned as for HA.
(ix) HCV NS3 helicase domain. The responses of 20 HCV patients reported by Gerlach et al. (25) were based on lymphocyte proliferation stimulated by 20-mer peptides spanning the NS3 helicase with 10-residue overlaps. The authors scored the response for a given patient as positive if the SI was greater than 3. Epitope scores were assigned to residues according to the number of responses to the peptide(s) that contained the residue. Residues that appeared in multiple peptides were assigned the average of the scores. The profile of raw epitope scores was smoothed with a 15-residue moving window average. The first-derivative profile generated as for HA by was smoothed with a 15-residue moving window average. Dominance peaks were assigned as for HA.
Identification of stable and unstable segments. (i) HA. A composite stability profile, including residue stabilities for the neutral-pH globular domain, a portion of the neutral-pH stem, and the entire low-pH fusogenic fragment, was used for the identification of stable and unstable segments. Stability profiles were calculated for the globular domain at neutral pH (residues 73 to 286 from PDB:2VIU, entropy factor 1.0), the stem at neutral pH (residues 25 to 72 and 287 to 524 from PDB:2VIU, entropy factor 0.92), and the fusogenic fragment (residues 27 to 32 and 385 to 507 from PDB:1HTM, chain D, entropy factor 0.80). For the identification of stability minima and maxima, the first-derivative profile was smoothed with a 15-residue moving window average and then the local minima and maxima were assigned to residues where the smoothed first derivative became positive and negative, respectively. For the H1, H5, and H9 stability profiles, COREX profiles were calculated for aligned portions of the corresponding structures. Input structures and entropy factors, respectively, were as follows: for the H1 globular domain, PDB:1RU7 and 1.009, and for the H1 stem, PDB:1RU7 and 0.932; for the H5 globular domain, PDB:1JSM and 1.008, and for the H5 stem, PDB:1JSM and 0.931; and for the H9 globular domain, PDB:1JSD and 0.980, and for the H9 stem, PDB:1JSD and 0.940.
(ii) Phl p 1. A stability profile was calculated for a structure of Phl p 1 that was generated by Swiss-Model (61) (http://swissmodel.expasy.org/) based on the incomplete crystal structures of Phl p 1 (PDB:1N10, chains A and B) and the complete amino acid sequence (accession no. P43213). Automated procedures in Swiss-Model constructed a chemically reasonable structure for a loop (residues 29 to 38 of the mature Phl p 1) that was missing from the experimentally determined structures. The COREX calculation employed an entropy factor of 0.850. For the identification of stability minima and maxima, the first-derivative profile was smoothed with a 15-residue moving window average.
(iii) M. tuberculosis Ag85A. A stability profile was calculated for the crystal structure of Ag85A (PDB:1SFR, chain A), using an entropy factor of 1.058. For the identification of stability minima and maxima, the first-derivative profile was smoothed with a 15-residue moving window average. The M. tuberculosis Ag85A sequence has 76% identity to the M. bovis Ag85B sequence.
(iv) M. tuberculosis Mpb70. A stability profile was calculated for the nuclear magnetic resonance (NMR) structure of Mpb70 (PDB:1NYO, model 1), using an entropy factor of 0.93 and default values for all other parameters. For the identification of stability minima and maxima, the first-derivative profile was smoothed with a 15-residue moving window average.
(v) M. bovis Ag85B. A stability profile was calculated for the crystal structure of M. tuberculosis Ag85B (PDB:1F0N, chain A), using an entropy factor of 1.058 and all other parameters set to default values. For the identification of stability minima and maxima, the first-derivative profile was smoothed with a 15-residue moving window average. The M. bovis Ag85B sequence has 99% identity to the M. tuberculosis Ag85B sequence and 76% identity to the M. tuberculosis Ag85A sequence.
(vi) M. tuberculosis Hsp65. A stability profile was calculated for a "homology model" of M. tuberculosis Hsp65 (accession no. P0A520) that was based on crystal structures of the same protein (PDB:1SJP, chains A and B). The homology model was prepared using Swiss-Model. Automated procedures in Swiss-Model constructed a chemically reasonable structure for a loop (residues 81 to 88) that was missing from the experimentally determined structures. Due to the large size of Hsp65, residue stabilities for the apical and equatorial domains were calculated separately and then combined into a composite stability profile. For the apical domain (residues 188 to 372), the COREX calculation employed an entropy factor of 0.952 and default values for all other parameters. For the equatorial domain (residues 1 to 187 and 373 to 525), the calculation used Monte Carlo sampling at 5,000 samples per partition, an entropy factor of 0.957, and default values for all other parameters. For the identification of stability minima and maxima, the first-derivative profile was smoothed with a 15-residue moving window average.
(vii) Bet v 1. A stability profile was calculated for the NMR structure of Bet v 1 (PDB:1BTV, model 1) using an entropy factor of 0.90 and default values for all other parameters. For the identification of stability minima and maxima, the first-derivative profile was smoothed with a 15-residue moving window average.
(viii) Ves v 5. A stability profile was calculated for the crystal structure of Ves v 5 (PDB:1QNX), using an entropy factor of 0.994 and default values for all other parameters. For the identification of stability minima and maxima, the first-derivative profile was smoothed with a seven-residue moving window average.
(ix) HCV NS3 helicase domain. A stability profile was calculated for a homology model of the NS3 helicase domain (HCV genotype 1a; accession no. P26664) that was based on the following structures (percent identity): PDB:1CU1 (99%), PDB:1HEI (94%), PDB:8OHM (96%), and PDB:2F55 (97%). The homology model was prepared using Swiss-Model. The stability profile was calculated using COREX with Monte Carlo sampling at 10,000 samples per partition, an entropy factor of 0.99, and all other parameters set to default values. For the identification of stability minima and maxima the first-derivative profile was smoothed with a 15-residue moving window average.
|
|
|---|
The most outstanding similarity in the CD4 T-cell epitope profiles is the coincidence of three immunodominant regions (centered at residues 215, 265, and 443) in the human and both mouse profiles (Fig. 1 and 2a). The human profile exhibits additional dominant regions (at residues 333 and 490) that are shared in common with the profile for H3-infected C57BL/6 mice. Humans have a distinct dominant region in the globular domain (at residue 130) that is not shared by either group of mice. The H3-infected C57BL/6 mice have a dominant region in the stem (at residue 390) that is not shared by the H1-infected HLA-DR1 mice but which is encompassed by a large shoulder in the peak of dominance in humans at residue 443. A dominant region in the globular domain (at residue 172) in the H1-infected HLA-DR1 mice appears to be shared with the H3-infected C57BL/6 mice but is shifted to a slightly more N-terminal position in the H3-infected mice.
![]() View larger version (31K): [in a new window] |
FIG. 2. Structural context of structurally stable segments and regions of conserved epitope dominance in the HA precursor and the low-pH conformation of the HA fusogenic fragment. In panel a, regions of above-average epitope score that were observed in humans and mice are shaded (dark gray). Note that the region of dominance centered at residue 443 is not shaded in the stem of the HA precursor. In panel b, regions of above-average residue stability are shaded (dark gray). Note that some stable regions in the stem of the precursor become unstable in the fusogenic fragment and vice versa. Arrows indicate established proteolytic cleavage sites in HA. Lines connect equivalent positions in the neutral-pH structure (PDB:1HA0) and fusogenic fragment (PDB:1HTM).
|
CD4+ epitopes adjacent to flexible, protease-sensitive HA segments. The similarity of HA epitope dominance profiles obtained for diverse subjects and different viral subtypes suggests that the conserved HA structure shapes the dominance profile by influencing HA processing and/or peptide loading. Conserved, dominant epitopes at residues 215, 265, 333, and 490 lie adjacent to established protease-sensitive sites in HA (Fig. 3b). The epitopes at residues 215 and 265 flank the tryptic cleavage site observed following residue 240 in the globular domain (74). This site becomes sensitive to trypsin only after HA undergoes the low-pH conformational change (63), which is expected in an antigen-processing compartment. The epitope at residue 333 lies at the C terminus of HA1, which is produced by a trypsin-like cleavage following residue 344 of the HA precursor in all infected tissues (74). The epitope at residue 490 lies at the C terminus of the HA ectodomain, which is the soluble fragment generated by bromelain cleavage of the membrane-bound protein (74). Local regions of conformational flexibility in the native or acid-induced conformation of HA most likely promote cleavage by trypsin and bromelain, and the same regions of flexibility are expected to promote cleavage by antigen-processing proteases and therefore facilitate presentation of adjacent epitopes.
![]() View larger version (28K): [in a new window] |
FIG. 3. Influenza virus HA epitopes cluster on the flanks of structurally stable HA segments. (a) Schematic diagram of HA subdomains (as in Fig. 1a). (b) Profiles of crystallographic B-factors for the subtype H3 precursor (heavy line) and fusogenic fragment (light line) superimposed on the profile of epitope scores for humans (as in Fig. 1b) and illustrating the locations of known proteolytic cleavage sites (arrows) (74). (c) Profiles of residue stability for the subtype H3 precursor (heavy line) and fusogenic fragment (light line) superimposed on the profile of epitope scores for humans (as in Fig. 1b). (d) Residue stability for four subtypes of HA. Residues are numbered according to the sequence of the HA precursor from the x31 strain of influenza virus (UniProtKB/Swiss-Prot accession no. P03437).
|
The correlations were evaluated on a peptide basis in order to properly weight the sampling frequency: i.e., the average epitope score for a peptide sequence was paired with the average B-factor for the same sequence or a similarly sized sequence at a given offset from the peptide. Using B-factors for the HA precursor, human epitopes were anticorrelated with B-factors, achieving an optimal value (rmax) of –0.60 (P = 0.003) with offset –10, suggesting that epitopes tend to be centered 10 residues C terminal from regions of low flexibility (Fig. 4a, upper panel).
![]() View larger version (22K): [in a new window] |
FIG. 4. Spatial relationship of CD4+ epitopes to stable antigen segments in influenza virus HA. (a) Correlation versus offset using epitope scores from studies of humans or H3-infected mice. The y axis indicates the center of the epitope, and the offset indicates the relative position of the stable/unstable segment being correlated (positive offset indicates C terminal). High B-factors indicate low stability, and thus correlations involving B-factor and residue stability generally have the opposite sign. In panel b, the average distances of mouse and human epitopes from the centers of stable segments are illustrated schematically.
|
Using B-factors from the uncleaved HA precursor, the epitope scores from H3-infected mice achieved values of rmax that indicate both correlation and anticorrelation, depending on the offset (Fig. 4a, lower panel). The rmax of 0.34 (P = 0.0007) at offset +16, suggests that epitopes tend to be located 16 residues N terminal from flexible, protease-sensitive sites. The rmax of –0.26 (P = 0.0123) at offset –20 suggests that epitopes tend to be located 20 residues C terminal from segments of low flexibility. This anticorrelation is similar to the anticorrelation observed for humans but at a larger offset from segments of low flexibility. Since dominance correlates well with nearby regions of low flexibility, these results raise the possibility that the positions of epitopes are related to the distribution of stability within the proteolytic fragments generated during antigen processing. As observed for the correlation with epitope scores from humans, the inclusion of B-factors from the fusogenic fragment did not improve the correlation with epitope scores from H3-infected mice (Fig. 4a, lower panel).
In spite of the statistical correlation of epitope scores and B-factors, it is difficult to discern the relationship by inspection of the overlaid profiles (Fig. 3b). Although B-factors represent the most well-documented structural correlate with proteolytic sensitivity, B-factors often do not reveal conformational fluctuations on time scales of seconds to minutes, which could be important during antigen processing (42).
CD4+ epitopes C terminal from stable HA segments. Amide-group hydrogen-deuterium exchange NMR (HX-NMR) can detect conformational fluctuations of the appropriate size and time scale for prediction of protease-sensitive sites (42, 71); however, HX-NMR data are not available for many proteins, especially large ones such as HA. The COREX algorithm developed by Hilser and Friere calculates a residue stability profile by computing the equilibrium constant for folding, Kf, of all possible segments in the protein (33). The stability (log Kf) profiles resemble the patterns of amide-group hydrogen-exchange protection observed by HX-NMR (31-34, 58). Residue stability profiles were calculated separately for monomers of the HA globular domain, stem, and fusogenic fragment. The residue stability profile for the complete HA monomer calculated at lower precision (due to computational overhead) was similar to the combined profiles for the globular domain and stem, and the profiles calculated for trimers of the stem and fusogenic fragment were similar to those of the monomers (data not shown). The acid-induced structure of the fusogenic fragment may be a better model for processing of the HA stem in the acidic lysosome. Thus, a composite stability profile, including residue stabilities for the neutral-pH globular domain, a portion of the neutral-pH stem, and the entire low-pH fusogenic fragment, also was compared to the epitope profile of HA.
Each peak of epitope dominance overlaps a structurally stable segment of HA (Fig. 3c). Optimum agreement between the number of stable segments and number of epitopes is obtained when HA's low-pH-induced conformational change is taken into account. The stability profile for the stem at neutral pH indicates only one peak in the region of the stem corresponding to the fusogenic fragment (positions 385 to 507), but the dominance profile indicates two peaks in this region. In contrast, the stability profile for the low-pH conformation of the fusogenic fragment indicates two peaks that overlap the epitopes at 443 and 490 (Fig. 3c). Correlation of epitope dominance with residue stability was evaluated for a range of offsets and using the stability profile from the neutral-pH structure or the composite profile that included residue stabilities from the low-pH fusogenic fragment. Epitope dominance in humans and H3-infected mice correlated at a low level with residue stability for the neutral-pH structure (rmax < 0.26) (Fig. 4a, upper and lower panels). The correlation was significantly improved by including residue stabilities from the fusogenic fragment. For dominance in humans, rmax reached 0.47 (P = 0.03) at offset –10 (Fig. 4a, upper panel). For the epitopes in H3-infected mice, the value of rmax reached 0.35 (P = 0.0005) at offset –16 (Fig. 4a, lower panel). The correlations with residue stability at offsets –10 and –16 are similar to the anticorrelations with B-factors at offsets –10 and –20 in humans and mice, respectively. Thus, HA epitopes in both humans and in the H3-infected mice tend to be located on the C-terminal flanks of stable segments, but in humans, the epitopes are shifted closer to the peaks of stability (Fig. 4b). Inclusion of residue stabilities from the fusogenic fragment improves the correlation because the fusogenic fragment has two stable segments, one on each side of the reverse turn formed at low pH, and an epitope lies on the C-terminal flank of each stable segment (Fig. 2b).
A model for antigen processing and peptide loading provides a framework for understanding CD4+ epitope immunodominance in HA (Fig. 5). In step 1, antigen three-dimensional structure is destabilized by acidification of the antigen-processing compartment and an endoprotease(s) cleaves near the center of the unstable segments, which generates polypeptide fragments having structurally frayed terminal segments. In step 2, a frayed C-terminal segment binds to the MHC protein, the stable segment unfolds, and the fragment adjusts its position in the MHC protein. This step may be aided by multiple DM-catalyzed cycles of fragment releasing and rebinding. DM reduces the selectivity of peptide binding to MHC proteins (48, 51), and therefore it may promote presentation of the most abundant proteolytic fragments. The final complex could be the product of an equilibrium that balances the stability of the MHC-fragment complex with the concentration and conformational stability of the free fragment. In step 3, unbound segments of the fragment are proteolytically trimmed and the completed MHC-peptide complex traffics to the surface of the antigen-presenting cell.
|
View larger version (10K): [in a new window] |
FIG. 5. General model for HA processing and presentation that accounts for CD4+ epitope dominance on the C-terminal flanks of stable antigen segments. In step 1, an endoprotease(s) cleaves in the unstable loops between stable domains. In step 2, the C-terminal region of a fragment cycles on and off of the MHC protein (catalyzed by DM), allowing the MHC protein to sample more unfolded conformations of the fragment. In step 3, the stably bound fragment is trimmed by proteases and the MHC-peptide complex traffics to the cell surface.
|
|
View this table: [in a new window] |
TABLE 1. Frequency of CD4+ epitope dominance peaks and structurally stable segments
|
![]() View larger version (35K): [in a new window] |
FIG. 6. Illustrations of CD4+ epitope dominance mapped onto an antigen's stable and unstable segments. (a) Schematic diagram and graph illustrating the frequency of epitope-dominance peaks in stable segments of allergen Phl p 1. In the schematic diagram, the approximate positions of peaks of stability (dark ovals) and epitope score (MHC proteins stylized with four circles) are illustrated. The number of peaks of the epitope score for each stable segment is indicated below the diagram. For the purposes of statistics, stable segments are defined as spanning between stability minima (indicated as broken vertical lines in the graph). (b) Schematic diagram and graph illustrating the locations of epitope dominance within unstable segments of M. tuberculosis Ag85A. The schematic diagram is as in panel a, except the arrows indicate distances (Df) from the stability minimum to peaks of epitope scores within the same unstable segment. Locations of stability minima (broken vertical lines) and stability maxima (solid vertical lines) are illustrated.
|
![]() View larger version (23K): [in a new window] |
FIG. 7. CD4+ epitope dominance mapped onto stable and unstable segments in three mycobacterial antigens. In panels a to c, schematic diagrams indicate the approximate positions of peaks of stability (dark ovals) and epitope scores (MHC proteins stylized with four circles) and graphs indicate residue stability (line graphs) and the frequency of epitopes (area graphs).
|
![]() View larger version (24K): [in a new window] |
FIG. 8. CD4+ epitope dominance mapped onto stable and unstable segments in two allergens and one viral antigen. In panels a to c, schematic diagrams indicate the approximate positions of peaks of stability (dark ovals) and epitope scores (MHC proteins stylized with four circles) and graphs indicate residue stability (line graphs) and the frequency of epitopes (area graphs).
|
In order to study the spatial distribution of epitope dominance relative to putative antigen-processing sites, the unstable segments (spanning between stability maxima) were aligned according to their stability minima. For each dominance peak, a flanking distance (Df) was measured as the distance in amino acid residues from the dominance peak to the stability minimum. For example, M. tuberculosis Ag85a has one dominance peak on the left side of the stability minimum (Df, <0) and five dominance peaks on the right side of the stability minimum (Df, >0) (Fig. 6b). For the 71 dominance peaks representing nine antigens, the distribution of dominance peaks between the left and right sides of the stability minimum is essentially even (34 peaks and 36 peaks, respectively). Thus, the tendency for HA epitopes to occur on the C-terminal flanks of stable segments does not reflect a general feature of the antigen-presentation machinery.
The basis for the tendency of HA epitopes to occur on the C-terminal flanks of stable segments could be related to the effect of the influenza virus M2 protein on endosomal and trans-Golgi vesicle acidification and sorting. The M2 protein is an acid-activated proton channel that selectively inhibits acidification of a subset of transport vesicles and partially blocks the delivery and maturation of lysosomal enzymes (29, 30). The M2 protein does not affect lysosomal pH, as indicated by normal morphology of the lysosomes and normal degradation of epidermal growth factor. Thus, it is possible that M2 alters the composition of the antigen-processing machinery in such a way that it introduces a bias toward loading of the C-terminal flanks of stable segments. The processing compartment could be deficient in a lysosomal protease(s) that is normally involved in degradation of the C-terminal flanks. Alternatively, the activity of DM could be affected. The existence of an HA-like pattern in other influenza virus proteins would support the M2-based explanation for a bias in presenting C-terminal flanks. Unfortunately, there are few data to consult. In the study by Woodland and coworkers in mice, aside from HA, only the acidic polymerase (PA) and nucleoprotein (NP) elicited significant T-cell responses (12), and there are no high-resolution structures available for either protein. Amino acid sequence conservation serves as a reasonably good surrogate for structural stability when analyzing the relationship to CD4+ epitope immunodominance (15). When amino acid sequence variation in PA and NP was evaluated in terms of sequence entropy (a surrogate for flexibility), significant correlations with epitope dominance were found at positive offsets. For PA, the rmax of 0.32 occurred at offset +14, and for NP, the rmax of 0.19 occurred at offset +4 (data not shown). Thus, influenza virus could be using M2 protein to focus the CD4+ responses onto the C-terminal flanks of its protein fragments.
Impact of HA variation in the context of structure-based CD4+ epitope dominance. Results presented here and in previous studies indicate that mechanisms of antigen processing disfavor presentation of epitopes in the middle of flexible antigen segments and favor presentation of epitopes in more stable antigen segments. It is likely that these mechanisms evolved in order to focus CD4+ T-cell immunity on antigen sequences that must be conserved to maintain protein structure and function. The advantage to the host is derived from the cost of immune escape, which is paid by the pathogen in reduced pathogen fitness (3, 5, 66). By extension of that argument, one might expect the processing and presentation machinery to target the most conserved, most stable segments for presentation, for example, by demanding high resistance to proteolysis. However, it is possible that such mechanisms would be too stringent in their requirements for stability, resulting in the presentation of too few epitopes—strong epitope dominance. Strong dominance probably is disadvantageous for the host because immune protection can be defeated by only a few mutations in the dominant epitopes.
The long-standing battle of influenza virus with vertebrate immune systems probably involves strategies to shape the CD4+ response. HA is the target for neutralizing antibodies, and it is a major antigen for both the murine and human CD4+ T-cell responses to influenza virus (74). Antigenic variation in HA arises by two mechanisms that render immunity ineffective against influenza: antigenic drift and antigenic shift (75). In antigenic drift, single amino acid changes accumulate in antibody epitopes of the prevailing viral subtype. In antigenic shift, the reassortment of genes from different viral subtypes produces a virus with more extensive changes in sequence and structure. As discussed below, CD4+ epitope dominance could mitigate or intensify losses in immune protection caused by antigenic drift, but dominance almost certainly intensifies the loss of protection caused by antigenic shift.
Antigenic drift is not expected to alter the CD4+ epitope dominance pattern because the small number of amino acid substitutions would have only minor effects on three-dimensional structure. Thus, the impact of antigenic drift on the recall of CD4+ T-cell responses depends on whether dominance is focused on conserved or divergent sequences. For example, epitopes in humans that are clustered near residue 443 are conserved, and thus responses to these epitopes are likely to be recalled by a new infection. In contrast, epitopes near residue 215 are poorly conserved, and thus dominant responses for that region of HA are less likely to be recalled in a new infection. Since particular CD4+ responses can reduce or exacerbate disease (7, 26, 35), HA is under selective pressure to evade the response to some epitopes more than others, and the dominance pattern could benefit the host or the virus. Mechanisms for epitope-specific CD4+ responses to influence viral clearance are not well understood, but epitope-specific CD8+ responses provide examples that could apply to CD4+ responses. On the one hand, the host benefits from more durable immunity if the dominant epitope is conserved for viral replicative fitness, such as has been observed for CD8+ T-cell epitopes in HIV (3, 5). On the other hand, the virus benefits if the immune response is somehow beneficial to the virus. For example, the CD8+ response to "nonprotective" T-cell epitopes can delay clearance of influenza virus (13).
When an influenza virus having a new HA subtype emerges, the allocation of CD4+ epitope dominance could determine whether the individual realizes any protection from exposure to the previously circulating viral subtype. Residue stability profiles for HA molecules of subtypes H1, H3, H5, and H9 illustrate striking differences in subdomain E and the N-terminal portion of RB (Fig. 3d), which correspond to the poorly conserved dominant regions at residues 130 and 172 (Fig. 1). Shifts of dominance in these regions could prevent the recall of CD4+ T-cell help for B cells specific for neutralizing antibody epitopes, which are most frequently observed in the same region of the protein (63, 74). Likewise, a weak recall of CD4+ T-cell help would reduce the response of CD8+ T cells that depend on CD4+ help (37). Further study of the relationship between HA structure and T-cell immunodominance is likely to reveal new insights on immunity to influenza virus and the feasibility of heterosubtypic immunization.
I thank the New Orleans Protein Folding Intergroup; James Robinson, Robert Garry, and Colin Watts for helpful discussions; and David Woodland, Jay Berzofsky, and Eli Sercarz for critical reviews of the manuscript.
Published ahead of print on 5 December 2007. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»