Listing of peptide prediction tools available on the web [Accessed November 2017].
Abstract
Positive and negative selection in the thymus relies on T-cell receptor recognition of peptides presented by HLA molecules and determines the repertoire of T cells. Immune competent T-lymphocytes target cells display nonself or pathogenic peptides in complex with their cognate HLA molecule. A peptide passes several selection processes before being presented in the peptide binding groove of an HLA molecule; here the sequence of the HLA molecule’s heavy chain determines the mode of peptide recruitment. During inflammatory processes, the presentable peptide repertoire is obviously altered compared to the healthy state, while the peptide loading pathway undergoes modifications as well. The presented peptides dictate the fate of the HLA expressing cell through their (1) sequence, (2) topology, (3) origin (self/nonself). Therefore, the knowledge about peptide competition and presentation in the context of alloreactivity, infection or pathogenic invasion is of enormous significance. Since in adoptive cellular therapies transferred cells should exclusively target peptide-HLA complexes they are primed for, one of the most crucial questions remains at what stage of viral infection viral peptides are presented preferentially over self-peptides. The systematic analyzation of peptide profiles under healthy or pathogenic conditions is the key to immunological success in terms of personalized therapeutics.
Keywords
- HLA
- peptides
- peptide prediction
- adoptive T-cell therapies
- peptide-vaccination
1. Introduction
The immune system of all species has to be able to discriminate self and foreign (nonself) antigens to combat infections without eliciting autoimmune diseases. The presentation of self and nonself occurs through displaying cellular proteins on the cell surface by proteins of the major histocompatibility complex (MHC) gene cluster. In humans, the MHC locus is termed human leukocyte antigen (HLA) and comprises several gene loci with numerous different alleles for most of the genes [1]. One part of the genes is subsumed as HLA class I (HLA-I) with the gene products being expressed on virtually every nucleated cell in the human body. HLA-I molecules present peptides of intracellular proteins on the cell surface. Cytotoxic T cells (CD8+ T cells) as part of the adaptive immune system can recognize these peptide HLA-I (pHLA-I) complexes by the T-cell receptor (TCR) and scan simultaneously the HLA molecule and the peptide [2, 3] to discriminate between healthy and unhealthy cells, for example, virally infected cells. At the same time, natural killer (NK) cells that are part of the innate immune system scan the cell surface of HLA as well. These cells become activated when HLA-I is missing on the cell surface, for example, on virally infected or tumor cells [4].
Peptide loading on HLA-I is a complex mechanism and determines in addition to the HLA allele which peptides will be presented. The central part of this process is the peptide loading complex that is localized in the endoplasmic reticulum (ER). The HLA-I molecule, consisting of a heavy chain and a microglobulin, is stabilized by several chaperons since the structure is unstable when no peptide is bound. The transporter associated with antigen processing (TAP) imports protein fragments that are degraded in the cytosol by the proteasome into the ER. Depending on the sequence, these peptides are trimmed in different ways in the ER [5, 6]. The bridge between TAP and HLA-I is the chaperon tapasin (TPN) that facilitates peptide binding in the peptide binding groove [7].
Because the HLA gene cluster ranks among the most polymorphic region in the human genome [1] and most of these polymorphisms are located in the peptide binding region (PBR) [8, 9], these polymorphisms result in an abundance of structurally different pHLA entities. In this chapter, we focus on the interplay between HLA alleles, bound peptides and the interaction with immune receptors. It is highlighted that even minor differences in the HLA sequence can impact on the bound ligand or the pHLA structure. Every single peptide changes the overall structure of the HLA molecule. Structural alterations that differ from self-pHLA structures will be recognized by the immune system. Therefore, the last parts of the chapter demonstrate the advantage of established immunopeptidomes for immunotherapies.
2. Peptide selection and presentation
The viability of the immune system is governed by interactions between effector cell receptors and their cognate antigenic ligands. Immune effector cells survey HLA-I molecules on the surface of antigen-presenting cells by indirectly scanning the proteomic content of every single cell. The fundamental role of CD8+ T cells, the elimination of pathogens, is elicited through HLA-I molecules complexed to a peptide of foreign (e.g., viral) origin.
Positive and negative selection of T cells in the thymus is a critical step for the development of a mature functional immune system. Immune cells that have not developed immune tolerance against specific pHLA-I complexes during thymus selection will recognize these antigens as foreign. The allele-specific and patient-specific peptides that are presented on the cell surface shape the individual immune response. Even a single alteration in the peptide sequence can be recognized by immune effectors. Single alterations in the sequence of the HLA heavy chain might not affect which peptide sequences can be bound but could lead to a modified overall pHLA-I structure or might affect the strength of peptide binding resulting in pHLA-I complexes with different half-life times. Besides the influence of amino acid (AA) exchanges within the heavy chain, peptides might undergo competition in patients who carry alleles with the same peptide-binding motif. For those reasons, the presentation of a given peptide is dependent on the HLA type and the health status of the patient. Half-life times of pHLA complexes influence the phenotype and/or functionality of T cells. That has to be taken into account when choosing pHLA-I targets for T-cell therapeutics. Based on these facts the importance of knowledge about allelic peptide specificity becomes obvious. Peptides have to pass several intracellular filters and processing steps before being presented to immune effectors cells. The bottleneck is the peptide loading complex (PLC). Not every available peptide in a cell would be necessarily bound to an HLA molecule and displayed at the cell surface since peptides undergo peptide competition during recruitment through the PLC. The PLC consists of several proteins; each has a specialized function for peptide selectivity and specificity. Proteins of this complex are dedicated targets for viral interference and thus viral immune evasion (Figure 1).
The HLA heavy chain has to adopt a peptide-receptive form and complex with certain proteins of the PLC. It could be demonstrated that certain allelic variants interact differently with proteins of the PLC and thus are more prone to present peptides of malignant origin. Those HLA subtypes can differ from alleles that are strictly dependent on the association with the PLC for peptide loading only by a single amino acid within the heavy chain, altering the structural interface that interacts with the PLC [10, 11, 12]. Especially the interaction of TPN, a protein that mediates the binding of high-affinity peptides into the HLA-I PBR, with the HLA-I molecule, is of exquisite importance to produce stable pHLA-I complexes that persist on the cell surface. However, few allelic HLA variants are able to present peptides without the assistance of TPN. That enables those alleles to continuously present viral peptides at an infectious stage where viral interference with TPN occurs; however, the presentation of self-peptides that did not take part in negative T-cell selection would be facilitated and might lead to uncontrollable autoimmune reactions. HLA-I variants that select and load peptides without the assistance of TPN are likely to present a broad range of low-affinity viral-derived peptides during an infection. However, the presentation of viral peptides during an active viral infection is a rare event since still self-peptides are present intracellularly and would compete with viral peptides to fit into the PBR. The viral peptides that would reach the cell surface complexed to an HLA-I molecule could hardly be predicted by peptide prediction tools.
3. Peptide specificity
Due to the fact that peptides undergo different selection steps before being presented by an individual HLA allele, it has to become clear that the individual HLA profile is the major and most distinguished obstacle. Every HLA allele differs from another by the composition of AAs in and/or outside the peptide-binding groove, resulting in allele-specific profiles of the bound peptides [13, 14, 15]. Therefore, the knowledge of an individual peptide binding profile can be used as precedence for the measurement of permissivity between HLA subtypes. The sequence, length and immunogenicity of a given peptide determine the half-life time of the whole pHLA complex and furthermore the specificity and reactivity of their cognate immune receptor.
Several studies demonstrated the impact of the sequence and feature of HLA-bound peptides on receptor recognition. This observation holds true for T-cell receptors of the adaptive immune system as well as for NK-cell receptors of the innate immune system. We recently could highlight the potential of the nonclassical HLA-I molecule HLA-E to select and present peptides of extraordinary length and their effect on differential NK-cell recognition. For the nonpolymorphic HLA-E molecule only two functional variants exist distinguished by a single AA difference. That would imply that HLA-E is, in regard to its peptide profile, invariant. However, by sequencing their bound ligands, we found both alleles presenting a different set of peptides [16]. Since HLA-E is an intermediate molecule for the adaptive and innate immune system, supporting the non-PLC-dependent presentation of peptides during HLA downregulation episodes, the invariability of this molecule would certainly make biological sense. However, the finding that HLA-E subtypes differ in their immunity was somehow unexpected. Reconstitution of empty HLA-E molecules with the designated peptides on the surface of artificial APCs resulted in a peptide-specific immune recognition [17].
Two types of NK-cell receptors interact with HLA-I molecules, killer cell immunoglobulin-like receptors (KIRs) and C-type lectin-like receptors. The former ones are highly polymorphic and polygenic in the population and recognize HLA-A, -B and -C alleles [18], whereas the latter ones bind among others to HLA-E [19, 20]. CD8+ T cells recognize endogenous HLA alleles and assess their immune status by virtue of the presented peptide. The display of different peptides thus allows for precise monitoring of the immune status of the cell through the adaptive immune system. However, this also means that these cells have to be primed on the recognition of specific peptides that are usually derived from endogenously expressed proteins [21]. The presented peptide repertoire can be altered by aberrant protein expression as well as the presence of foreign proteins (e.g., viral proteins) in the cell. To counter recognition by CD8+ T cells, many viruses have developed immune evasion strategies that specifically target HLA peptide loading and presentation. For instance, the HCMV protein US6 interferes with peptide translocation at the ER thus depriving the available peptide pool or even directly procures that the HLA heavy chain is degraded through US2 or US11 [22]. However, in the event of such disrupted HLA-I presentation, NK cells become activated. Although NK cells do not recognize the specific HLA allele, the absence of HLA expression on the cell surface triggers NK-cell activation. Nevertheless, NK cells can still recognize certain peptides in the context of the nonclassical HLA molecule HLA-E that presents a very narrow set of peptides derived from the signal sequence of other HLA class I molecules. These peptides are 9 AA in length and are anchored preferably by Met at peptide position p2 and Leu at pΩ, whereas positions p4, p5 and p6 are accessible to the solvent [23, 24]. On NK cells, HLA-E in combination with these peptides is recognized by inactivating the NKG2A/CD94 heterodimeric receptor complex [25]. In the absence of HLA class I molecules or during cell stress, HLA-E was shown to present noncanonical peptides of different length [16, 17], for example, the Hsp60-derived peptide QMRPVSRVL that causes loss of recognition by NKG2A/CD94 [26] or the HIV Gag-derived peptide AISPRTLNA that causes HLA-E upregulation [27]. In the case of HCMV, a peptide from the UL40 protein that closely resembles the sequence of leader peptides from certain HLA-C allotypes is provided to stabilize HLA-E expression in infected cells. However, in individuals negative for these HLA-C allotypes, the UL40-peptide constitutes the presentation of nonself on HLA-E and can thus elicit a CD8+ T-cell response [19, 28]. Additionally, HLA-E in complex with other pathogen-derived peptides was shown to stimulate CD8+ T-cell responses. For instance, the Epstein–Barr virus-derived peptide SQAPLPCVL was shown to be recognized by the αβTCR of a CD8+-CD94/NKG2C+ T-cell clone [29, 30] described HLA-E-restricted
4. Peptide binding prediction, bioinformatic tools
To identify peptides that would be suitable for application in cellular therapeutic strategies, certain properties have to be analyzed: (1) the peptide-binding motif of the HLA allele of choice and (2) the HLA allele-specific features of the bound peptides such as length and topology. There are several bioinformatic tools that enable scientists and clinicians to predict peptides that would be presented by a certain HLA allele, yet, these tools do not consider allele-specific features and immune dominance of peptides. The kinetics of antigen expression and the competition of peptides to be preferentially bound and presented are also not considered by these bioinformatics prediction tools. Most data available in these tools are based on experimental peptide data (Tables 1 and 2). However, it remains unclear if those predicted peptides would ever be naturally presented. Peptides predicted from for example a viral protein would not necessarily be processed, selected and/or presented by the respective patient awaiting T-cell therapy. Therefore, the pathogen- or peptide-specific T cells that would be transplanted might not be able to find their mutual pHLA molecule. An example of the first successful adoptive transfer of virus-specific T cells described the transfer of HCMV-specific T cells and their reconstitution of antiviral immunity in an immune-deficient bone marrow transplant recipient [31]. The technique of adoptive T-cell transfer could be further improved leading to the selection of specific T-cells based on IFN-γ secretion or pHLA multimer staining and selection following antigen stimulation [32, 33, 34]. Both techniques bear the imperative to know which viral peptides are presented on the particular HLA subtype of, for example, HCMV-infected cells. So far, few HLA-restricted peptides have been studied. The majority of peptides are derived from the well-characterized phosphoprotein (pp)65 or the immediately early (IE)1 protein, however, not for every patient responses against these two proteins are immunodominant [35, 36]. Best studied are the pp65-derived peptides NLVPMVATV and TPRVTGGGAM, restricted to HLA-A*02:01 and HLA-B*07:02, respectively. Those peptides are described to induce extremely strong T-cell responses [37, 38, 39, 40, 41, 42]. Yet, these peptides have been computationally predicted [43] but not been isolated from HLA molecules. Thus, it remains unproven if they would ever be naturally presented. That might be an explanation for the failure of long-term T-cell transfers [33, 44, 45].
Name | Application | Methods | Ref. | Number of HLA class I alleles | Number of HLA class II alleles | Peptide length | Other species | URL |
---|---|---|---|---|---|---|---|---|
BIMAS | Predicts half-time of dissociation of peptides from protein sequences | Coefficient tables | [58] | 41 inc. supertypes | 0 | 8–10 | No | https://www-bimas.cit.nih.gov/molbio/hla_bind/ |
EpiJen | Predicts peptide binding from protein sequence (proteasome cleavage, TAP binding and MHC binding) | Multi-step algorithm | [59] | 18 | 0 | 9 | No | http://www.ddg-pharmfac.net/epijen/EpiJen/EpiJen.htm |
hla_a2_smm | Predicts binding affinity of peptides, high affinity HLA-A2 binding peptides from protein sequence and mutated peptides with higher affinity | SMM pair coefficients | [60] | 1 | 0 | 9–10 | No | https://zlab.bu.edu/SMM/ |
IEDB T Cell Epitope Prediction Tools | Predicts T cell epitopes from proteins (MHC binding, processing and immunogenicity) | Several tools can be chosen | [61, 62, 63, 64] | 77 | n/s | Class I:8–14 Class II: n/s |
Chimpanzee, cow, gorilla, macaque, mouse, pig, rat for MHC class I; mouse for MCH class II | http://tools.iedb.org/main/tcell/ |
Mappp | Predicts antigenic peptides to be processed and presented by MHC class I from peptide or protein sequence | Uses BIMAS or SYFPEITHI for binding prediction | [65] | 35 inc. supertypes | 0 | 8–10 | Mouse, cattle |
http://www.mpiib-berlin.mpg.de/MAPPP/index.html
http://www.mpiib-berlin.mpg.de/MAPPP/binding.html |
MHC2MIL | Predicts binding affinity of MHC-II peptides from protein sequence | MIL | [66] | 0 | 26 | 9–25 | No | http://datamining-iip.fudan.edu.cn/service/MHC2MIL/index.html |
MHC2PRED | Prediction of MHC class II binders | SVM | [67] | 0 | 38 inc. supertypes | 9 | Mouse | http://crdd.osdd.net/raghava/mhc2pred/index.html |
MHCBN | Database with information about allele specific MHC binding peptides, MHC nonbinding, TAP binding, TAP nonbinding peptides and T-cell epitopes | Database | [68, 69] | n/s | n/s | n/s | http://crdd.osdd.net/raghava/mhcbn/index.html | |
MHCMIR | Predicts binding affinity and levels of MHC-II peptides from peptide or protein sequence | MIR | [70] | 0 | 13 | All | Mouse | http://ailab.ist.psu.edu/mhcmir/predict.html |
MHCPRED | Predicts binding affinity of peptides to MHC class I and II molecules and to TAP from protein sequence and calculates binding affinity for heteroclitic peptides | Additive method, partial least square regression | [71, 72, 73] | 11 | 3 | 9 | Mouse |
http://www.ddg-pharmfac.net/mhcpred/MHCPred/
http://www.ddg-pharmfac.net/mhcpred/MHCPred/pepLib.html |
MMBPred | Predicts mutated high affinity and promiscuous MHC class-I binding peptides from protein sequence, epitope enhancement, 1–3 AAs mutation of nonamer peptides | QM | [74] | 40 inc. supertypes | 0 | 9 | Rhesus macaque, mouse | http://crdd.osdd.net/raghava/mmbpred/ |
MULTIPRED | Predicts binding of peptides to HLA class I and class II DR supertypes and individual genotypes | Uses NetMHCpan and NetMHCIIpan | [75] | 13 supertypes | 13 supertypes | 8–11 for HLA class I and genotype 9 for HLA class II |
No | http://cvc.dfci.harvard.edu/multipred2/index.php |
NetCTL | Predicts CTL epitopes in protein sequences (Cleavage, TAP transport, HLA class I binding) | ANN | [76] | 12 supertypes | 0 | 9 | No | http://www.cbs.dtu.dk/services/NetCTL/ |
NetMHC | Predicts binding of peptide to MHC class I molecules from peptide or protein sequence | ANN | [19, 77] | 81 (or 12 supertypes) | 0 | 8–14 | Chimpanzee, rhesus macaque, mouse, cuttle, pig | http://www.cbs.dtu.dk/services/NetMHC/ |
NetMHCcons | Predicts binding of peptides to any known MHC class I molecule from peptide or protein sequence | Consensus (NetMHC, NetMHCpan and PickPocket) | [78] | User specified | 0 | 8–15 | Chimpanzee, gorilla, rhesus macaque, mouse, cuttle, pig | http://www.cbs.dtu.dk/services/NetMHCcons/ |
NetMHCII | Predicts binding of peptides to HLA-DR, HLA-DQ, HLA-DP from peptide or protein sequence | ANN | [79, 80] | 0 | 26 | variable | Mouse | http://www.cbs.dtu.dk/services/NetMHCII/ |
NetMHCIIpan | Predicts binding of peptides to HLA-DR, HLA-DQ, HLA-DP from peptide or protein sequence | ANN | [81, 82] | User specified | 0 | Mouse | http://www.cbs.dtu.dk/services/NetMHCIIpan/ | |
NetMHCpan | Predicts binding of peptides to any known MHC class I molecule from peptide or protein sequence | ANN | [19, 83, 84] | User specified | 0 | 8–14 | Chimpanzee, gorilla, rhesus macaque, mouse, cuttle, pig | http://www.cbs.dtu.dk/services/NetMHCpan/ |
nHLAPred: ANNPred | Predicts MHC Class I binding regions in proteins | ANN | [85] | 26 inc. supertypes | 0 | Mouse | http://crdd.osdd.net/raghava/nhlapred/neural.html | |
nHLAPred: ComPred | Predicts MHC Class I binding regions in proteins | ANN/QM | [85] | 59 inc. supertypes | 0 | Rhesus macaque, mouse | http://crdd.osdd.net/raghava/nhlapred/comp.html | |
PREDPEP | Predicts binding of peptides to HLA class I from peptide or protein sequence | Published coefficient tables | [86] | 6 | 0 | 8–10 (dependent on the allele) | Mouse | http://margalit.huji.ac.il/Teppred/mhc-bind/index.html |
ProPred | Predicts MHC Class II binding regions in an antigen sequence | QM | [87] | 51 | 0 | No | http://crdd.osdd.net/raghava/propred/ | |
ProPred I | Predicts MHC Class I binding regions in an antigen sequence | QM | [88] | 39 inc. supertypes | 0 | Mouse, cattle | http://crdd.osdd.net/raghava/propred1/index.html | |
Rankpep | Predicts binding of peptides to MHC class I and class II molecules from peptide or protein sequence | PSSM | [89, 90, 91] | n/s | n/s | Dependent on the allele | http://imed.med.ucm.es/Tools/rankpep.html | |
svmhc | Predicts binding of peptides to MHC class I molecules from peptide or protein sequence | SVM, uses MHCPEP or SYFPEITHI | [92, 93] | 31 | 0 | 8–10 (dependent on the allele) | No | http://svmhc.bioinfo.se/svmhc// |
SYFPEITHI | Database of MHC ligands and peptide motifs and epitope prediction | Matrix/motif-based, published motifs | [94] | 33 | 6 | 8–11 for HLA class I; 15 for HLA class II | No | http://www.syfpeithi.de/ |
TEPITOPEpan | Predicts tissue-specific binding of peptides to MHC class II molecules from peptide or protein sequence | PSSM | [95] | 0 | 50 | 9–25 | No | http://datamining-iip.fudan.edu.cn/service/TEPITOPEpan/index.html |
Name | Underlying database/data source | URL for matrices/training data |
---|---|---|
BIMAS | Coefficient tables deduced from the published literature by Dr. Kenneth Parker, Children’s Hospital Boston |
https://www-bimas.cit.nih.gov/cgi-bin/molbio/hla_coefficient_viewing_page
https://www-bimas.cit.nih.gov/molbio/hla_bind/hla_references.html |
EpiJen | AntiJen [96, 97], SYFPEITHI [94] | http://www.ddg-pharmfac.net/antijen/AntiJen/antijenhomepage.htm |
hla_a2_smm | BIMAS [58], SYFPEITHI [94], data described in Peters, Tong [60] | https://zlab.bu.edu/SMM/ |
IEDB T Cell Epitope Prediction Tools | IEDB [61], Sette lab, Buus lab, uses diverse predictions methods (see webpage) |
http://tools.iedb.org/mhci/download/
http://tools.iedb.org/mhcii/download/ |
Mappp | BIMAS [58], SYFPEITHI [94], coefficient tables deduced from the literature by Kenneth Parker, Children’s Hospital Boston | — |
MHC2MIL | Data by Wang, Sidney [98] | — |
MHC2PRED | JenPep [19], MHCBN [68] | — |
MHCBN | MHCBN [68] | — |
MHCMIR | IEDB [61] | — |
MHCPRED | JenPep [19] | — |
MMBPred | MHCBN [68] | — |
MULTIPRED | See NetMHCpan and NetMHCIIpan | — |
NetCTL | See NetMHC | — |
NetMHC | Trained for 81 HLA alleles including HLA-A, -B, -C and –E, n/s | — |
NetMHCcons | IEDB [61] | — |
NetMHCII | Data by [19] | http://www.cbs.dtu.dk/suppl/immunology/NetMHCII-2.0.php |
NetMHCIIpan | IEDB [61] | http://www.cbs.dtu.dk/suppl/immunology/NetMHCIIpan-3.0/ |
NetMHCpan | IEDB [61], IMGT/HLA database [1] | — |
nHLAPred: ANNPred | MHCBN [68] | — |
nHLAPred: ComPred | MHCBN [68], BIMAS [58] | http://crdd.osdd.net/raghava/nhlapred/matrix.html |
PREDPEP | Pairwise potential table by Miyazawa and Jernigan [99] | — |
ProPred | QMs by Sturniolo, Bono [95] | http://crdd.osdd.net/raghava/propred/page4.html |
ProPred I | BIMAS [58] and matrices by Ruppert, Sidney [100] and Sidney, Southwood [101] | http://crdd.osdd.net/raghava/propred1/matrices/matrix.html |
Rankpep | MHCPEP [102], SYFPEITHI [94], GenBank [103] | — |
svmhc | MHCPEP [102], SYFPEITHI [94] | http://www.cs.cornell.edu/people/tj/svm_light/ |
SYFPEITHI | Published literature | — |
TEPITOPEpan | n/s | http://datamining-iip.fudan.edu.cn/service/TEPITOPEpan/TEPITOPEpan.html |
5. Analysis of naturally presented peptides
The analysis of the individual patient and cell-type-specific immunopeptidome can be realized through sequencing the HLA-bound peptides. It is imperative for all ongoing peptide studies and cellular therapies to find peptides that are (1) naturally presented by the distinct allele, (2) immunogenic for (at best) a public T-cell repertoire and (3) preferentially presented when different peptides are available. A study from Yaciuk et al. showed for example that the peptides isolated from HIV-infected T cells differ from predicted peptides and exhibit different T-cell reactions, factors that have to be considered in designing immunotherapies [46]. That information is only available after immunopeptidome analyses.
In the past, different methods have been applied to answer these questions comprehensively. There are two reliable methods to determine peptide sequences from selected HLA alleles. First, membrane-bound HLA molecules from recombinant single-antigen-presenting cells [47, 48] or from donor cells [49, 50] can be captured by affinity chromatographic methods and the bound peptides isolated and sequenced by mass spectrometry. Second, the most realizable method is the soluble HLA technology [16, 51]. Vectors encoding for soluble forms of HLA molecules (Exon 1–4) are transfected or lentivirally transduced into the cell line of choice. An optional recombinant tag (e.g., V5 tag) engineered at the C-terminus of the protein enables specific purification of the recombinant HLA molecule of choice without the challenge of contamination by cellular-self-HLA molecules. Both methods have been compared by Scull et al. [52] and indicated as an equivalent for the determination of allele-specific peptides. Furthermore, Badrinath et al. [10] could demonstrate that sHLA molecules associate during peptide acquisition with the loading complex as well. These results prove evidence that the use of sHLA technology for understanding allele-specific peptide-binding motifs, the prerequisite for updating peptide prediction databases, is the most time- and cost-efficient implementation.
For the development of tailor-made T-cell-based immunotherapeutic strategies, the identification of tumor-specific HLA ligands is imperative. The production of recombinant sHLA-expressing cells derived from various tissues of malignant origins would guide towards understanding immune dominance through peptide competition. One of the most innovative applications is the peptide fishing from tumor tissue. Immunological tolerance is mediated through T cells that are primed in the thymus by self-peptides. Therefore, the comprehensive knowledge of the HLA immunopeptidome from diseased cells is fundamental for the development of efficient immunotherapeutic strategies. The presentation of peptides depends on the health state of a patient. During infections, the expression of HLA molecules and thus peptide presentation, including presentation of self-peptides, is diminished through an immune escape mechanism of the invasive pathogen.
6. Peptide vaccination
The treatment of cancer represents a great challenge due to the fact that the vast majority of HLA-restricted peptides differs from tissue to tissue and is dependent on the tumor entity. For that reason, it becomes obvious how fundamentally important the knowledge of the tumor-specific peptidome is. For personalized cancer immunotherapies, the knowledge of naturally presented peptides [53] represents the exclusive possibility for therapeutical success. The analysis of the mutanome, the proteomic content of a diseased cell, includes the discovery of neo-antigens or post-translational-modified peptides and the avoidance of targeting self-antigens from healthy tissue. The results of such individual mutanomes might alter during the course of tumor progression [16]. In peptide vaccination trials, the use of multiple peptides in combination [54, 55] represents a useful method for targeting all MHC-presenting cells with the peptide of choice. Yet, since the cell type where the peptide(s) bind to cannot be traced, the rates of antitumor immune responses might differ from patient to patient and certain tumor cells where, for example, low MHC expression rates might remain undetected from the immune system. To achieve a comprehensive and precise analysis of presented tumor antigens, the method of antigen discovery and appropriate T-cell assay for knowledge of immunogenicity of the dedicated antigen for vaccination is the key factor [56, 57].
7. Conclusion
Peptide selection and presentation is an exquisite biological and immunological event. Every single peptide is a mirror of the health state of a distinct cell and determines the outcome of immune recognition and responses. For all cellular therapies, the knowledge of the HLA-subtype specific proteome is crucial for the utilization of ligand prediction tools, which have to be implemented where no experimental data are available, yet.
References
- 1.
Robinson J et al. The IPD and IMGT/HLA database: Allele variant databases. Nucleic Acids Research. 2015; 43 (Database issue):D423-D431 - 2.
Germain RN. MHC-dependent antigen processing and peptide presentation: Providing ligands for T lymphocyte activation. Cell. 1994; 76 (2):287-299 - 3.
Zinkernagel RM, Doherty PC. Immunological surveillance against altered self components by sensitised T lymphocytes in lymphocytic choriomeningitis. Nature. 1974; 251 (5475):547-548 - 4.
Kiessling R et al. Evidence for a similar or common mechanism for natural killer cell activity and resistance to hemopoietic grafts. European Journal of Immunology. 1977; 7 (9):655-663 - 5.
Roelse J et al. Trimming of TAP-translocated peptides in the endoplasmic reticulum and in the cytosol during recycling. The Journal of Experimental Medicine. 1994; 180 (5):1591-1597 - 6.
Mpakali A et al. Structural basis for antigenic peptide recognition and processing by endoplasmic reticulum (ER) aminopeptidase 2. The Journal of Biological Chemistry. 2015; 290 (43):26021-26032 - 7.
Zarling AL et al. Tapasin is a facilitator, not an editor, of class I MHC peptide binding. The Journal of Immunology. 2003; 171 (10):5287-5295 - 8.
Carreno BM et al. The peptide binding specificity of HLA class I molecules is largely allele-specific and non-overlapping. Molecular Immunology. 1992; 29 (9):1131-1140 - 9.
Bjorkman PJ et al. The foreign antigen binding site and T cell recognition regions of class I histocompatibility antigens. Nature. 1987; 329 (6139):512-518 - 10.
Badrinath S et al. Position 156 influences the peptide repertoire and tapasin dependency of human leukocyte antigen B*44 allotypes. Haematologica. 2012; 97 (1):98-106 - 11.
Badrinath S et al. A micropolymorphism altering the residue triad 97/114/156 determines the relative levels of Tapasin independence and distinct peptide profiles for HLA-A(*)24 allotypes. Journal of Immunology Research. 2014; 2014 :298145 - 12.
Manandhar T et al. Understanding the obstacle of incompatibility at residue 156 within HLA-B*35 subtypes. Immunogenetics. 2016; 68 (4):247-260 - 13.
Bade-Doeding C et al. Amino acid 95 causes strong alteration of peptide position Pomega in HLA-B*41 variants. Immunogenetics. 2007; 59 (4):253-259 - 14.
Badrinath S et al. Position 45 influences the peptide binding motif of HLA-B*44:08. Immunogenetics. 2012; 64 (3):245-249 - 15.
Huyton T et al. Residue 81 confers a restricted C-terminal peptide binding motif in HLA-B*44:09. Immunogenetics. 2012; 64 (9):663-668 - 16.
Kraemer T et al. HLA-E: Presentation of a broader peptide repertoire impacts the cellular immune response-implications on HSCT outcome. Stem Cells International. 2015; 2015 :346714 - 17.
Celik AA et al. The diversity of the HLA-E-restricted peptide repertoire explains the immunological impact of the Arg107Gly mismatch. Immunogenetics. 2016; 68 (1):29-41 - 18.
Uhrberg M et al. Human diversity in killer cell inhibitory receptor genes. Immunity. 1997; 7 (6):753-763 - 19.
Hoare HL et al. Structural basis for a major histocompatibility complex class Ib-restricted T cell response. Nature Immunology. 2006; 7 (3):256-264 - 20.
Brooks AG et al. Specific recognition of HLA-E, but not classical, HLA class I molecules by soluble CD94/NKG2A and NK cells. Journal of Immunology. 1999; 162 (1):305-313 - 21.
Bjorkman PJ et al. Structure of the human class I histocompatibility antigen, HLA-A2. Nature. 1987; 329 (6139):506-512 - 22.
Jackson SE, Mason GM, Wills MR. Human cytomegalovirus immunity and immune evasion. Virus Research. 2011; 157 (2):151-160 - 23.
Braud V, Jones EY, McMichael A. The human major histocompatibility complex class Ib molecule HLA-E binds signal sequence-derived peptides with primary anchor residues at positions 2 and 9. European Journal of Immunology. 1997; 27 (5):1164-1169 - 24.
Petrie EJ et al. CD94-NKG2A recognition of human leukocyte antigen (HLA)-E bound to an HLA class I leader sequence. The Journal of Experimental Medicine. 2008; 205 (3):725-735 - 25.
Braud VM et al. HLA-E binds to natural killer cell receptors CD94/NKG2A, B and C. Nature. 1998; 391 (6669):795-799 - 26.
Michaelsson J et al. A signal peptide derived from hsp60 binds HLA-E and interferes with CD94/NKG2A recognition. The Journal of Experimental Medicine. 2002; 196 (11):1403-1414 - 27.
Nattermann J et al. HIV-1 infection leads to increased HLA-E expression resulting in impaired function of natural killer cells. Antiviral Therapy. 2005; 10 (1):95-107 - 28.
Heatley SL et al. Polymorphism in human cytomegalovirus UL40 impacts on recognition of human leukocyte antigen-E (HLA-E) by natural killer cells. The Journal of Biological Chemistry. 2013; 288 (12):8679-8690 - 29.
Garcia P et al. Human T cell receptor-mediated recognition of HLA-E. European Journal of Immunology. 2002; 32 (4):936-944 - 30.
Salerno-Goncalves R et al. Identification of a human HLA-E-restricted CD8+ T cell subset in volunteers immunized with Salmonella enterica serovar Typhi strain Ty21a typhoid vaccine. Journal of Immunology. 2004; 173 (9):5852-5862 - 31.
Riddell SR et al. Restoration of viral immunity in immunodeficient humans by the adoptive transfer of T cell clones. Science. 1992; 257 (5067):238-241 - 32.
Cobbold M et al. Adoptive transfer of cytomegalovirus-specific CTL to stem cell transplant patients after selection by HLA-peptide tetramers. The Journal of Experimental Medicine. 2005; 202 (3):379-386 - 33.
Feuchtinger T et al. Adoptive transfer of pp65-specific T cells for the treatment of chemorefractory cytomegalovirus disease or reactivation after haploidentical and matched unrelated stem cell transplantation. Blood. 2010; 116 (20):4360-4367 - 34.
Schmitt A et al. Adoptive transfer and selective reconstitution of streptamer-selected cytomegalovirus-specific CD8+ T cells leads to virus clearance in patients after allogeneic peripheral blood stem cell transplantation. Transfusion. 2011; 51 (3):591-599 - 35.
Elkington R et al. Ex vivo profiling of CD8+−T-cell responses to human cytomegalovirus reveals broad and multispecific reactivities in healthy virus carriers. Journal of Virology. 2003; 77 (9):5226-5240 - 36.
Sylwester AW et al. Broadly targeted human cytomegalovirus-specific CD4+ and CD8+ T cells dominate the memory compartments of exposed subjects. The Journal of Experimental Medicine. 2005; 202 (5):673-685 - 37.
Ameres S et al. Presentation of an immunodominant immediate-early CD8+ T cell epitope resists human cytomegalovirus immunoevasion. PLoS Pathogens. 2013; 9 (5):e1003383 - 38.
Gibson L et al. Human cytomegalovirus proteins pp65 and immediate early protein 1 are common targets for CD8+ T cell responses in children with congenital or postnatal human cytomegalovirus infection. Journal of Immunology. 2004; 172 (4):2256-2264 - 39.
Kato R et al. Early detection of cytomegalovirus-specific cytotoxic T lymphocytes against cytomegalovirus antigenemia in human leukocyte antigen haploidentical hematopoietic stem cell transplantation. Annals of Hematology. 2015; 94 (10):1707-1715 - 40.
Nguyen TH et al. Cross-reactive anti-viral T cells increase prior to an episode of viral reactivation post human lung transplantation. PLoS One. 2013; 8 (2):e56042 - 41.
van Bockel D et al. Validation of RNA-based molecular clonotype analysis for virus-specific CD8+ T-cells in formaldehyde-fixed specimens isolated from peripheral blood. Journal of Immunological Methods. 2007; 326 (1-2):127-138 - 42.
Yang X et al. Structural basis for clonal diversity of the public T cell response to a dominant human cytomegalovirus epitope. The Journal of Biological Chemistry. 2015; 290 (48):29106-29119 - 43.
Kuzushima K et al. Efficient identification of HLA-A*2402-restricted cytomegalovirus-specific CD8(+) T-cell epitopes by a computer algorithm and an enzyme-linked immunospot assay. Blood. 2001; 98 (6):1872-1881 - 44.
Doubrovina E et al. Adoptive immunotherapy with unselected or EBV-specific T cells for biopsy-proven EBV+ lymphomas after allogeneic hematopoietic cell transplantation. Blood. 2012; 119 (11):2644-2656 - 45.
Gottschalk S et al. An Epstein-Barr virus deletion mutant associated with fatal lymphoproliferative disease unresponsive to therapy with virus-specific CTLs. Blood. 2001; 97 (4):835-843 - 46.
Yaciuk JC et al. Direct interrogation of viral peptides presented by the class I HLA of HIV-infected T cells. Journal of Virology. 2014; 88 (22):12992-13004 - 47.
Krausa P et al. Definition of peptide binding motifs amongst the HLA-A*30 allelic group. Tissue Antigens. 2000; 56 (1):10-18 - 48.
Macdonald WA et al. A naturally selected dimorphism within the HLA-B44 supertype alters class I structure, peptide repertoire, and T cell recognition. The Journal of Experimental Medicine. 2003; 198 (5):679-691 - 49.
Kruger T et al. Lessons to be learned from primary renal cell carcinomas: Novel tumor antigens and HLA ligands for immunotherapy. Cancer Immunology, Immunotherapy. 2005; 54 (9):826-836 - 50.
Weinzierl AO et al. Distorted relation between mRNA copy number and corresponding major histocompatibility complex ligand density on the cell surface. Molecular & Cellular Proteomics. 2007; 6 (1):102-113 - 51.
Kunze-Schumacher H, Blasczyk R, Bade-Doeding C. Soluble HLA technology as a strategy to evaluate the impact of HLA mismatches. Journal of Immunology Research. 2014; 2014 :246171 - 52.
Scull KE et al. Secreted HLA recapitulates the immunopeptidome and allows in-depth coverage of HLA a*02:01 ligands. Molecular Immunology. 2012; 51 (2):136-142 - 53.
Bassani-Sternberg M, Coukos G. Mass spectrometry-based antigen discovery for cancer immunotherapy. Current Opinion in Immunology. 2016; 41 :9-17 - 54.
Slingluff CL Jr. The present and future of peptide vaccines for cancer: Single or multiple, long or short, alone or in combination? Cancer Journal. 2011; 17 (5):343-350 - 55.
Li W et al. Peptide vaccine: Progress and challenges. Vaccines (Basel). 2014; 2 (3):515-536 - 56.
Purcell AW, Croft NP, Tscharke DC. Immunology by numbers: Quantitation of antigen presentation completes the quantitative milieu of systems immunology! Current Opinion in Immunology. 2016; 40 :88-95 - 57.
Caron E et al. Analysis of major histocompatibility complex (MHC) Immunopeptidomes using mass spectrometry. Molecular & Cellular Proteomics. 2015; 14 (12):3105-3117 - 58.
Parker KC, Bednarek MA, Coligan JE. Scheme for ranking potential HLA-A2 binding peptides based on independent binding of individual peptide side-chains. Journal of Immunology. 1994; 152 (1):163-175 - 59.
Doytchinova IA, Guan PP, Flower DR. EpiJen: A server for multistep T cell epitope prediction. BMC Bioinformatics. 2006; 7 :131-142 - 60.
Peters B et al. Examining the independent binding assumption for binding of peptide epitopes to MHC-I molecules. Bioinformatics. 2003; 19 (14):1765-1772 - 61.
Vita R et al. The immune epitope database (IEDB) 3.0. Nucleic Acids Research. 2015; 43 (Database issue):D405-D412 - 62.
Kim Y et al. Immune epitope database analysis resource. Nucleic Acids Research. 2012; 40 (Web Server issue):W525-W530 - 63.
Wang P et al. A systematic assessment of MHC class II peptide binding predictions and evaluation of a consensus approach. PLoS Computational Biology. 2008; 4 (4):e1000048 - 64.
Moutaftsi M et al. A consensus epitope prediction approach identifies the breadth of murine T(CD8+)-cell responses to vaccinia virus. Nature Biotechnology. 2006; 24 (7):817-819 - 65.
Hakenberg J et al. MAPPP: MHC class I antigenic peptide processing prediction. Applied Bioinformatics. 2003; 2 (3):155-158 - 66.
Xu Y et al. MHC2MIL: A novel multiple instance learning based method for MHC-II peptide binding prediction by considering peptide flanking region and residue positions. BMC Genomics. 2014; 15 (Suppl 9):S9 - 67.
Bhasin M, Raghava GP. SVM based method for predicting HLA-DRB1*0401 binding peptides in an antigen sequence. Bioinformatics. 2004; 20 (3):421-423 - 68.
Lata S, Bhasin M, Raghava GP. MHCBN 4.0: A database of MHC/TAP binding peptides and T-cell epitopes. BMC Research Notes. 2009; 2 :61 - 69.
Bhasin M, Singh H, Raghava GP. MHCBN: A comprehensive database of MHC binding and non-binding peptides. Bioinformatics. 2003; 19 (5):665-666 - 70.
EL-Manzalawy Y, Dobbs D, Honavar V. Predicting MHC-II binding affinity using multiple instance regression. IEEE-ACM Transactions on Computational Biology and Bioinformatics. 2011; 8 (4):1067-1079 - 71.
Guan P et al. MHCPred: A server for quantitative prediction of peptide-MHC binding. Nucleic Acids Research. 2003; 31 (13):3621-3624 - 72.
Guan P et al. MHCPred: Bringing a quantitative dimension to the online prediction of MHC binding. Applied Bioinformatics. 2003; 2 (1):63-66 - 73.
Hattotuwagama CK et al. Quantitative online prediction of peptide binding to the major histocompatibility complex. Journal of Molecular Graphics & Modelling. 2004; 22 (3):195-207 - 74.
Bhasin M, Raghava GP. Prediction of promiscuous and high-affinity mutated MHC binders. Hybridoma and Hybridomics. 2003; 22 (4):229-234 - 75.
Zhang GL et al. MULTIPRED2: A computational system for large-scale identification of peptides predicted to bind to HLA supertypes and alleles. Journal of Immunological Methods. 2011; 374 (1-2):53-61 - 76.
Larsen MV et al. Large-scale validation of methods for cytotoxic T-lymphocyte epitope prediction. BMC Bioinformatics. 2007; 8 :424 - 77.
Andreatta M, Nielsen M. Gapped sequence alignment using artificial neural networks: Application to the MHC class I system. Bioinformatics. 2016; 32 (4):511-517 - 78.
Karosiene E et al. NetMHCcons: A consensus method for the major histocompatibility complex class I predictions. Immunogenetics. 2012; 64 (3):177-186 - 79.
Andreatta M et al. NNAlign: A web-based prediction method allowing non-expert end-user discovery of sequence motifs in quantitative peptide data. PLoS One. 2011; 6 (11):e26781 - 80.
Lundegaard C, Lund O, Nielsen M. Accurate approximation method for prediction of class I MHC affinities for peptides of length 8, 10 and 11 using prediction tools trained on 9mers. Bioinformatics. 2008; 24 (11):1397-1398 - 81.
Karosiene E et al. NetMHCIIpan-3.0, a common pan-specific MHC class II prediction method including all three human MHC class II isotypes, HLA-DR, HLA-DP and HLA-DQ. Immunogenetics. 2013; 65 (10):711-724 - 82.
Andreatta M et al. Accurate pan-specific prediction of peptide-MHC class II binding affinity with improved binding core identification. Immunogenetics. 2015; 67 (11-12):641-650 - 83.
Hoof I et al. NetMHCpan, a method for MHC class I binding prediction beyond humans. Immunogenetics. 2009; 61 (1):1-13 - 84.
Jurtz V et al. NetMHCpan-4.0: Improved peptide-MHC class I interaction predictions integrating eluted ligand and peptide binding affinity data. Journal of Immunology. 2017; 199 (9):3360-3368 - 85.
Bhasin M, Raghava GP. A hybrid approach for predicting promiscuous MHC class I restricted T cell epitopes. Journal of Biosciences. 2007; 32 (1):31-42 - 86.
Schueler-Furman O et al. Structure-based prediction of binding peptides to MHC class I molecules: Application to a broad range of MHC alleles. Protein Science. 2000; 9 (9):1838-1846 - 87.
Singh H, Raghava GP. ProPred: Prediction of HLA-DR binding sites. Bioinformatics. 2001; 17 (12):1236-1237 - 88.
Singh H, Raghava GP. ProPred1: Prediction of promiscuous MHC class-I binding sites. Bioinformatics. 2003; 19 (8):1009-1014 - 89.
Reche PA et al. Enhancement to the RANKPEP resource for the prediction of peptide binding to MHC molecules using profiles. Immunogenetics. 2004; 56 (6):405-419 - 90.
Reche PA, Reinherz EL. Prediction of peptide-MHC binding using profiles. Methods in Molecular Biology. 2007; 409 :185-200 - 91.
Reche PA, Glutting JP, Reinherz EL. Prediction of MHC class I binding peptides using profile motifs. Human Immunology. 2002; 63 (9):701-709 - 92.
Donnes P, Elofsson A. Prediction of MHC class I binding peptides, using SVMHC. BMC Bioinformatics. 2002; 3 :25 - 93.
Donnes P, Kohlbacher O. SVMHC: A server for prediction of MHC-binding peptides. Nucleic Acids Research. 2006; 34 (Web Server issue):W194-W197 - 94.
Rammensee H et al. SYFPEITHI: Database for MHC ligands and peptide motifs. Immunogenetics. 1999; 50 (3-4):213-219 - 95.
Sturniolo T et al. Generation of tissue-specific and promiscuous HLA ligand databases using DNA microarrays and virtual HLA class II matrices. Nature Biotechnology. 1999; 17 (6):555-561 - 96.
Blythe MJ, Doytchinova IA, Flower DR. JenPep: A database of quantitative functional peptide data for immunology. Bioinformatics. 2002; 18 (3):434-439 - 97.
McSparron H et al. JenPep: A novel computational information resource for immunobiology and vaccinology. Journal of Chemical Information and Computer Sciences. 2003; 43 (4):1276-1287 - 98.
Wang P et al. Peptide binding predictions for HLA DR, DP and DQ molecules. BMC Bioinformatics. 2010; 11 :568 - 99.
Miyazawa S, Jernigan RL. Residue-residue potentials with a favorable contact pair term and an unfavorable high packing density term, for simulation and threading. Journal of Molecular Biology. 1996; 256 (3):623-644 - 100.
Ruppert J et al. Prominent role of secondary anchor residues in peptide binding to HLA-A2.1 molecules. Cell. 1993; 74 (5):929-937 - 101.
Sidney J et al. Specificity and degeneracy in peptide binding to HLA-B7-like class I molecules. Journal of Immunology. 1996; 157 (8):3480-3490 - 102.
Brusic V, Rudy G, Harrison LC. MHCPEP, a database of MHC-binding peptides: Update 1997. Nucleic Acids Research. 1998; 26 (1):368-371 - 103.
Benson DA et al. GenBank. Nucleic Acids Research. 2003; 31 (1):23-27 - 104.
Ahn K et al. The ER-luminal domain of the HCMV glycoprotein US6 inhibits peptide translocation by TAP. Immunity. 1997; 6 (5):613-621 - 105.
Park B et al. Human cytomegalovirus inhibits tapasin-dependent peptide loading and optimization of the MHC class I peptide cargo for immune evasion. Immunity. 2004; 20 (1):71-85 - 106.
Hill A et al. Herpes simplex virus turns off the TAP to evade host immunity. Nature. 1995; 375 (6530):411-415 - 107.
Vambutas A et al. Interaction of human papillomavirus type 11 E7 protein with TAP-1 results in the reduction of ATP-dependent peptide transport. Clinical Immunology. 2001; 101 (1):94-99 - 108.
Bennett EM et al. Cutting edge: Adenovirus E19 has two mechanisms for affecting class I MHC expression. Journal of Immunology. 1999; 162 (9):5049-5052 - 109.
Horst D et al. Specific targeting of the EBV lytic phase protein BNLF2a to the transporter associated with antigen processing results in impairment of HLA class I-restricted antigen presentation. Journal of Immunology. 2009; 182 (4):2313-2324 - 110.
Beier DC et al. Association of human class I MHC alleles with the adenovirus E3/19K protein. Journal of Immunology. 1994; 152 (8):3862-3872 - 111.
Kerkau T et al. The human immunodeficiency virus type 1 (HIV-1) Vpu protein interferes with an early step in the biosynthesis of major histocompatibility complex (MHC) class I molecules. The Journal of Experimental Medicine. 1997; 185 (7):1295-1305 - 112.
Gruhler A, Peterson PA, Fruh K. Human cytomegalovirus immediate early glycoprotein US3 retains MHC class I molecules by transient association. Traffic. 2000; 1 (4):318-325 - 113.
Coscoy L, Ganem D. Kaposi’s sarcoma-associated herpesvirus encodes two proteins that block cell surface display of MHC class I chains by enhancing their endocytosis. Proceedings of the National Academy of Sciences of the United States of America. 2000; 97 (14):8051-8056 - 114.
Schwartz O et al. Endocytosis of major histocompatibility complex class I molecules is induced by the HIV-1 Nef protein. Nature Medicine. 1996; 2 (3):338-342 - 115.
Wiertz EJ et al. The human cytomegalovirus US11 gene product dislocates MHC class I heavy chains from the endoplasmic reticulum to the cytosol. Cell. 1996; 84 (5):769-779 - 116.
Wiertz EJ et al. Sec61-mediated transfer of a membrane protein from the endoplasmic reticulum to the proteasome for destruction. Nature. 1996; 384 (6608):432-438 - 117.
Gilbert MJ et al. Cytomegalovirus selectively blocks antigen processing and presentation of its immediate-early gene product. Nature. 1996; 383 (6602):720-722 - 118.
Yin Y, Manoury B, Fahraeus R. Self-inhibition of synthesis and antigen presentation by Epstein-Barr virus-encoded EBNA1. Science. 2003; 301 (5638):1371-1374