The structural basis for the increased immunogenicity of two HIV-reverse transcriptase peptide variant/class I major histocompatibility complexes

Vaccinomics to Design a Multiepitope Vaccine against Legionella pneumophila

Legionella pneumophila is found in the natural aquatic environment and can resist a wide range of environmental conditions. There are around fifty species of Legionella, at least twenty-four of which are directly linked to infections in humans. L. pneumophila is the cause of Legionnaires' disease, a potentially lethal form of pneumonia. By blocking phagosome-lysosome fusion, L. pneumophila lives and proliferates inside macrophages. For this disease, there is presently no authorized multiepitope vaccine available. For the multi-epitope-based vaccine (MEBV), the best antigenic candidates were identified using immunoinformatics and subtractive proteomic techniques. Several immunoinformatics methods were utilized to predict B and T cell epitopes from vaccine candidate proteins. To construct an in silico vaccine, epitopes (07 CTL, 03 HTL, and 07 LBL) were carefully selected and docked with MHC molecules (MHC-I and MHC-II) and human TLR4 molecules. To increase the immunological response, the vaccine was combined with a 50S ribosomal adjuvant. To maximize vaccine protein expression, MEBV was cloned and reverse-translated in Escherichia coli. To prove the MEBV's efficacy, more experimental validation is required. After its development, the resulting vaccine is greatly hoped to aid in the prevention of L. pneumophila infections.

Integrated Core Proteomics, Subtractive Proteomics, and Immunoinformatics Investigation to Unveil a Potential Multi-Epitope Vaccine against Schistosomiasis

Schistosomiasis is a parasitic infection that causes considerable morbidity and mortality in the world. Infections of parasitic blood flukes, known as schistosomes, cause the disease. No vaccine is available yet and thus there is a need to design an effective vaccine against schistosomiasis. Schistosoma japonicum, Schistosoma mansoni, and Schistosoma haematobium are the main pathogenic species that infect humans. In this research, core proteomics was combined with a subtractive proteomics pipeline to identify suitable antigenic proteins for the construction of a multi-epitope vaccine (MEV) against human-infecting Schistosoma species. The pipeline revealed two antigenic proteins-calcium binding and mycosubtilin synthase subunit C-as promising vaccine targets. T and B cell epitopes from the targeted proteins were predicted using multiple bioinformatics and immunoinformatics databases. Seven cytotoxic T cell lymphocytes (CTL), three helper T cell lymphocytes (HTL), and four linear B cell lymphocytes (LBL) epitopes were fused with a suitable adjuvant and linkers to design a 217 amino-acid-long MEV. The vaccine was coupled with a TLR-4 agonist (RS-09; Sequence: APPHALS) adjuvant to enhance the immune responses. The designed MEV was stable, highly antigenic, and non-allergenic to human use. Molecular docking, molecular dynamics (MD) simulations, and molecular mechanics/generalized Born surface area (MMGBSA) analysis were performed to study the binding affinity and molecular interactions of the MEV with human immune receptors (TLR2 and TLR4) and MHC molecules (MHC I and MHC II). The MEV expression capability was tested in an Escherichia coli (strain-K12) plasmid vector pET-28a(+). Findings of these computer assays proved the MEV as highly promising in establishing protective immunity against the pathogens; nevertheless, additional validation by in vivo and in vitro experiments is required to discuss its real immune-protective efficacy.

Crystal structure of the giant panda MHC class I complex: First insights into the viral peptide presentation profile in the bear family

The viral cytotoxic T lymphocyte (CTL) epitope peptides presented by classical MHC-I molecules require the assembly of a peptide-MHC-I-β2m (pMHC-I) trimolecular complex for T cell receptor (TCR) recognition, which is the critical activation link for triggering antiviral T cell immunity. Research on T cell immunology in the Ursidae family, especially structural immunology, is still lacking. In this study, the structure of the key trimolecular complex pMHC-I, which binds a peptide from canine distemper virus, was solved for the first time using giant panda as a representative species of Ursidae. The structural characteristics of the giant panda pMHC-I complex (pAime-128), including the unique pockets in the peptide-binding groove (PBG), were analyzed in detail. Comparing the pAime-128 to others in the bear family and extending the comparison to other mammals revealed distinct features. The interaction between MHC-I and β2m, the features of pAime-128 involved in TCR docking and cluster of differentiation 8 (CD8) binding, the anchor sites in the PBG, and the CTL epitopes of potential viruses that infect pandas were clarified. Unique features of pMHC-I viral antigen presentation in the panda were revealed by solving the three-dimensional (3D) structure of pAime-128. The distinct characteristics of pAime-128 indicate an unusual event that emerged during the evolution of the MHC system in the bear family. These results provide a new platform for research on panda CTL immunity and the design of vaccines for application in the bear family.

T cell receptor interactions with human leukocyte antigen govern indirect peptide selectivity for the cancer testis antigen MAGE-A4

T cell-mediated immunity is governed primarily by T cell receptor (TCR) recognition of peptide-human leukocyte antigen (pHLA) complexes and is essential for immunosurveillance and disease control. This interaction is generally stabilized by interactions between the HLA surface and TCR germline-encoded complementarity-determining region (CDR) loops 1 and 2, whereas peptide selectivity is guided by direct interactions with the TCR CDR3 loops. Here, we solved the structure of a newly identified TCR in complex with a clinically relevant peptide derived from the cancer testis antigen melanoma antigen-A4 (MAGE-A4). The TCR bound pHLA in a position shifted toward the peptide's N terminus. This enabled the TCR to achieve peptide selectivity via an indirect mechanism, whereby the TCR sensed the first residue of the peptide through HLA residue Trp-167, which acted as a tunable gateway. Amino acid substitutions at peptide position 1 predicted to alter the HLA Trp-167 side-chain conformation abrogated TCR binding, indicating that this indirect binding mechanism is essential for peptide recognition. These findings extend our understanding of the molecular rules that underpin antigen recognition by TCRs and have important implications for the development of TCR-based therapies.

Modulation of the strength and character of HIV-specific CD8+ T cell responses with heteroclitic peptides

Chronic infection with human immunodeficiency virus (HIV) causes HIV-specific CD8⁺ T cell dysfunction and exhaustion. The strong association between non-progression and maintenance of HIV-specific CD8⁺ T cell cytokine production and proliferative capacities suggests that invigorating CD8⁺ T cell immune responses would reduce viremia and slow disease progression. A series of studies have demonstrated that sequence variants of native immunogenic peptides can generate more robust CD8⁺ T cell responses and that stimulation with these ‘heteroclitic’ peptides can steer responses away from the phenotypic and functional attributes of exhaustion acquired during chronic HIV infection. Incorporation of heteroclitic peptide stimulation within therapeutic vaccines could favour induction of more effective cellular antiviral responses, and in combination with ‘shock and kill’ strategies, contribute towards HIV cure.

Understanding CD8+ T-cell responses toward the native and alternate HLA-A*02:01-restricted WT1 epitope

The Wilms' tumor 1 (WT1) antigen is expressed in solid and hematological malignancies, but not healthy tissues, making it a promising target for cancer immunotherapies. Immunodominant WT1 epitopes, the native HLA-A2/WT1_126-134 (RMFPNAPYL) (HLA-A2/RMFPNAPYL epitope (WT1A)) and its modified variant YMFPNAPYL (HLA-A2/YMFPNAPYL epitope (WT1B)), can induce WT1-specific CD8⁺ T cells, although WT1B is more stably bound to HLA-A*02:01. Here, to further determine the benefits of those two targets, we assessed the naive precursor frequencies; immunogenicity and cross-reactivity of CD8⁺ T cells directed toward these two WT1 epitopes. Ex vivo naive WT1A- and WT1B-specific CD8⁺ T cells were detected in healthy HLA-A*02:01⁺ individuals with comparable precursor frequencies (1 in 10⁵–10⁶) to other naive CD8⁺ T-cell pools (for example, A2/HIV-Gag_77-85), but as expected, ~100 × lower than those found in memory populations (influenza, A2/M1_58-66; EBV, A2/BMLF1_280-288). Importantly, only WT1A-specific naive precursors were detected in HLA-A2.1 mice. To further assess the immunogenicity and recruitment of CD8⁺ T cells responding to WT1A and WT1B, we immunized HLA-A2.1 mice with either peptide. WT1A immunization elicited numerically higher CD8⁺ T-cell responses to the native tumor epitope following re-stimulation, although both regimens produced functionally similar responses toward WT1A via cytokine analysis and CD107a expression. Interestingly, however, WT1B immunization generated cross-reactive CD8⁺ T-cell responses to WT1A and could be further expanded by WT1A peptide revealing two distinct populations of single- and cross-reactive WT1A⁺CD8⁺ T cells with unique T-cell receptor-αβ gene signatures. Therefore, although both epitopes are immunogenic, the clinical benefits of WT1B vaccination remains debatable and perhaps both peptides may have separate clinical benefits as treatment targets.

Enhancing Human Immunodeficiency Virus-Specific CD8(+) T Cell Responses with Heteroclitic Peptides

Human immunodeficiency virus (HIV)-specific CD8(+) T cells play a critical role in containing HIV replication and delaying disease progression. However, HIV-specific CD8(+) T cells become progressively more "exhausted" as chronic HIV infection proceeds. Symptoms of T cell exhaustion range from expression of inhibitory receptors and selective loss of cytokine production capacity through reduced proliferative potential, impaired differentiation into effector cells and increased susceptibility to apoptosis. While effective combination antiretroviral therapy (cART) durably reduces HIV viremia to undetectable levels, this alone does not restore the full pluripotency of HIV-specific CD8(+) T cells. In a number of studies, a subset of peptide epitope variants categorized as heteroclitic, restimulated more potent cellular immune responses in vitro than did the native, immunizing peptides themselves. This property of heteroclitic peptides has been exploited in experimental cancer and chronic viral infection models to promote clearance of transformed cells and persistent viruses. In this review, we consider the possibility that heteroclitic peptides could improve the efficacy of therapeutic vaccines as part of HIV immunotherapy or eradication strategies. We review literature on heteroclitic peptides and illustrate their potential to beneficially modulate the nature of HIV-specific T cell responses toward those found in the small minority of HIV-infected, aviremic cART-naïve persons termed elite controllers or long-term non-progressors. Our review suggests that the efficacy of HIV vaccines could be improved by identification, testing, and incorporation of heteroclitic variants of native HIV peptide epitopes.

Antigenic peptide molecular recognition by the DRB1–DQB1 haplotype modulates multiple sclerosis susceptibility

DRB1–DQB1 binding affinities in peptide molecular recognition process. (A) In protective haplotype DRB1 allele displays a preferential affinity for MBP peptide, while (B) in predisposing haplotype DQB1 allele displays a preferential affinity for EBNA1 peptide.

Dissecting the dynamics of HIV-1 protein sequence diversity

The rapid mutation of human immunodeficiency virus-type 1 (HIV-1) and the limited characterization of the composition and incidence of the variant population are major obstacles to the development of an effective HIV-1 vaccine. This issue was addressed by a comprehensive analysis of over 58,000 clade B HIV-1 protein sequences reported over at least 26 years. The sequences were aligned and the 2,874 overlapping nonamer amino acid positions of the viral proteome, each a possible core binding domain for human leukocyte antigen molecules and T-cell receptors, were quantitatively analyzed for four patterns of sequence motifs: (1) "index", the most prevalent sequence; (2) "major" variant, the most common variant sequence; (3) "minor" variants, multiple different sequences, each with an incidence less than that of the major variant; and (4) "unique" variants, each observed only once in the alignment. The collective incidence of the major, minor, and unique variants at each nonamer position represented the total variant population for the position. Positions with more than 50% total variants contained correspondingly reduced incidences of index and major variant sequences and increased minor and unique variants. Highly diverse positions, with 80 to 98% variant nonamer sequences, were present in each protein, including 5% of Gag, and 27% of Env and Nef, each. The multitude of different variant nonamer sequences (i.e. nonatypes; up to 68%) at the highly diverse positions, represented by the major, multiple minor, and multiple unique variants likely supported variants function both in immune escape and as altered peptide ligands with deleterious T-cell responses. The patterns of mutational change were consistent with the sequences of individual HXB2 and C1P viruses and can be considered applicable to all HIV-1 viruses. This characterization of HIV-1 protein mutation provides a foundation for the design of peptide-based vaccines and therapeutics.

In silico quantitative prediction of peptides binding affinity to human MHC molecule: an intuitive quantitative structure-activity relationship approach

In this paper, we have handpicked 23 kinds of electronic properties, 37 kinds of steric properties, 54 kinds of hydrophobic properties and 5 kinds of hydrogen bond properties from thousands of amino acid structural and property parameters. Principal component analysis (PCA) was applied on these parameters and thus ten score vectors involving significant nonbonding properties of 20 coded amino acids were yielded, called the divided physicochemical property scores (DPPS) of amino acids. The DPPS descriptor was then used to characterize the structures of 152 HLA-A*0201-restricted CTL epitopes, and significant variables being responsible for the binding affinities were selected by genetic algorithm, and a quantitative structure-activity relationship (QSAR) model by partial least square was established to predict the peptide-HLA-A*0201 molecule interactions. Statistical analysis on the resulted DPPS-based QSAR models were consistent well with experimental exhibits and molecular graphics display. Diversified properties of the different residues in binding peptides may contribute remarkable effect to the interactions between the HLA-A*0201 molecule and its peptide ligands. Particularly, hydrophobicity and hydrogen bond of anchor residues of peptides may have a significant contribution to the interactions. The results showed that DPPS can well represent the structural characteristics of the antigenic peptides and is a promising approach to predict the affinities of peptide binding to HLA-A*0201 in a efficient and intuitive way. We expect that this physical-principle based method can be applied to other protein-peptide interactions as well.

Variation in cytotoxic T-lymphocyte responses to peptides derived from tyrosinase-related protein-2

In this study, we developed three optimized peptide ligands (OPL) that demonstrate increased affinities for HLA-A*0201 compared with wild-type tyrosinase-related protein-2 (TRP-2) peptide. The OPL contain amino acids from TRP-2((180-188)) and preferred primary and auxiliary HLA-A*0201 anchor residues. Cytotoxic T lymphocyte (CTL) lines were generated against wild-type TRP-2 peptide and OPL by multiple rounds of peptide stimulation of peripheral blood mononuclear cells from HLA-A2*0201(+) healthy individuals. CTL reactivity profiles to three different OPL were donor-dependent. Among donors, at least one OPL was particularly stimulatory and elicited high levels of CTL that cross-reacted with wild-type TRP-2 peptide. Cytotoxicity assays using CTL raised on wild-type TRP-2 peptide or OPL demonstrated lysis of HLA-A2-positive glioblastoma cells. Molecular models of TRP-2 and OPL peptides docked with HLA-A*0201 demonstrated that substitution of F for S at position 1 (P1) oriented the peptides favoring a pi-pi aromatic interaction with W 167 of HLA-A*0201. This in turn positions P5 and P8 aromatic rings to face solvent that may promote binding to the T-cell receptor, leading to a robust T-cell activation. The results of this study further substantiate the concept that rational design and testing of multiple peptides for the same T-cell epitope should elicit a broader response among different individuals than single peptide immunization. Our results may partially explain why some patients have better clinical responses to peptide-based immunotherapy, whereas others respond poorly.

Structures of MART-126/27-35 Peptide/HLA-A2 complexes reveal a remarkable disconnect between antigen structural homology and T cell recognition

Small structural changes in peptides presented by major histocompatibility complex (MHC) molecules often result in large changes in immunogenicity, supporting the notion that T cell receptors are exquisitely sensitive to antigen structure. Yet there are striking examples of TCR recognition of structurally dissimilar ligands. The resulting unpredictability of how T cells will respond to different or modified antigens impacts both our understanding of the physical bases for TCR specificity as well as efforts to engineer peptides for immunomodulation. In cancer immunotherapy, epitopes and variants derived from the MART-1/Melan-A protein are widely used as clinical vaccines. Two overlapping epitopes spanning amino acid residues 26 through 35 are of particular interest: numerous clinical studies have been performed using variants of the MART-1 26-35 decamer, although only the 27-35 nonamer has been found on the surface of targeted melanoma cells. Here, we show that the 26-35 and 27-35 peptides adopt strikingly different conformations when bound to HLA-A2. Nevertheless, clonally distinct MART-1(26/27-35)-reactive T cells show broad cross-reactivity towards these ligands. Simultaneously, however, many of the cross-reactive T cells remain unable to recognize anchor-modified variants with very subtle structural differences. These dichotomous observations challenge our thinking about how structural information on unligated peptide/MHC complexes should be best used when addressing questions of TCR specificity. Our findings also indicate that caution is warranted in the design of immunotherapeutics based on the MART-1 26/27-35 epitopes, as neither cross-reactivity nor selectivity is predictable based on the analysis of the structures alone.

Template-Controlled Face-to-Face Stacking of Olefinic and Aromatic Carboxylic Acids in the Solid State

A rigid 1,8-dipyridylnaphthalene template for supramolecular organization of unsaturated dicarboxylic acids in the solid state has been developed. Self-assembly of the template and aromatic dicarboxylic acids generates nondistorted paracyclophane and pyridinophane architectures with perfectly superimposed pi-systems. The dipyridyl template can be used to promote stereoselective solid-state dimerization of olefinic dicarboxylic acids.

A Structure-based, Quantitative Structure?Activity Relationship Approach for Predicting HLA-A*0201-restricted Cytotoxic T Lymphocyte Epitopes

In this investigation, we first constructed four types of non-bonding interaction matrixes by defining direct contacting residue types for HLA-A*0201 protein in interaction with each position of HLA-A*0201-restricted cytotoxic T lymphocyte epitope as well as several formulae calculating ligand/receptor non-bonding interactions. Relative to these studies, a method which we refer to as structure-based, quantitative structure-activity relationship is proposed and utilized for studies on 266 HLA-A*0201-restricted cytotoxic T lymphocyte epitopes. The resulting genetic algorithm-partial least square regression model is consistent with both published studies and molecular graphics analysis. Two non-bonding interactions (i.e. hydrophobic and hydrogen bonding), are found to play important roles in antigen recognition and presentation, especially exerting effects at positions of anchor residues in antigen peptides.

Ab initio prediction of peptide-MHC binding geometry for diverse class I MHC allotypes

Since determining the crystallographic structure of all peptide-MHC complexes is infeasible, an accurate prediction of the conformation is a critical computational problem. These models can be useful for determining binding energetics, predicting the structures of specific ternary complexes with T-cell receptors, and designing new molecules interacting with these complexes. The main difficulties are (1) adequate sampling of the large number of conformational degrees of freedom for the flexible peptide, (2) predicting subtle changes in the MHC interface geometry upon binding, and (3) building models for numerous MHC allotypes without known structures. Whereas previous studies have approached the sampling problem by dividing the conformational variables into different sets and predicting them separately, we have refined the Biased-Probability Monte Carlo docking protocol in internal coordinates to optimize a physical energy function for all peptide variables simultaneously. We also imitated the induced fit by docking into a more permissive smooth grid representation of the MHC followed by refinement and reranking using an all-atom MHC model. Our method was tested by a comparison of the results of cross-docking 14 peptides into HLA-A*0201 and 9 peptides into H-2K(b) as well as docking peptides into homology models for five different HLA allotypes with a comprehensive set of experimental structures. The surprisingly accurate prediction (0.75 A backbone RMSD) for cross-docking of a highly flexible decapeptide, dissimilar to the original bound peptide, as well as docking predictions using homology models for two allotypes with low average backbone RMSDs of less than 1.0 A illustrate the method's effectiveness. Finally, energy terms calculated using the predicted structures were combined with supervised learning on a large data set to classify peptides as either HLA-A*0201 binders or nonbinders. In contrast with sequence-based prediction methods, this model was also able to predict the binding affinity for peptides to a different MHC allotype (H-2K(b)), not used for training, with comparable prediction accuracy.

Identification of interleukin-13 receptor alpha2 peptide analogues capable of inducing improved antiglioma CTL responses

Restricted and high-level expression of interleukin-13 receptor alpha2 (IL-13Ralpha2) in a majority of human malignant gliomas makes this protein an attractive vaccine target. We have previously described the identification of the IL-13Ralpha2(345-353) peptide as a human leukocyte antigen-A2 (HLA-A2)-restricted CTL epitope. However, as it remains unclear how efficiently peptide-based vaccines can induce specific CTLs in patients with malignant gliomas, we have examined whether analogue epitopes could elicit heteroclitic antitumor T-cell responses versus wild-type peptides. We have created three IL-13Ralpha2 analogue peptides by substitutions of the COOH-terminal isoleucine (I) for valine (V) and the NH(2)-terminal tryptophan (W) for either alanine (A), glutamic acid (E), or nonsubstituted (W; designated as 1A9V, 1E9V, and 9V, respectively). In comparison with the native IL-13Ralpha2 epitope, the analogue peptides 9V and 1A9V displayed higher levels of binding affinity and stability in HLA-A2 complexes and yielded an improved stimulatory index for patient-derived, specific CTLs against the native epitope expressed by HLA-A2(+) glioma cells. In HLA-A2-transgenic HHD mice, immunization with the peptides 9V and 1A9V induced enhanced levels of CTL reactivity and protective immunity against an intracranial challenge with IL13Ralpha2-expressing syngeneic tumors when compared with vaccines containing the native IL-13Ralpha2 epitope. These findings indicate highly immunogenic IL-13Ralpha2 peptide analogues may be useful for the development of vaccines capable of effectively expanding IL-13Ralpha2-specific, tumor-reactive CTLs in glioma patients.

Structural prediction of peptides binding to MHC class I molecules

Peptide binding to class I major histocompatibility complex (MHCI) molecules is a key step in the immune response and the structural details of this interaction are of importance in the design of peptide vaccines. Algorithms based on primary sequence have had success in predicting potential antigenic peptides for MHCI, but such algorithms have limited accuracy and provide no structural information. Here, we present an algorithm, PePSSI (peptide-MHC prediction of structure through solvated interfaces), for the prediction of peptide structure when bound to the MHCI molecule, HLA-A2. The algorithm combines sampling of peptide backbone conformations and flexible movement of MHC side chains and is unique among other prediction algorithms in its incorporation of explicit water molecules at the peptide-MHC interface. In an initial test of the algorithm, PePSSI was used to predict the conformation of eight peptides bound to HLA-A2, for which X-ray data are available. Comparison of the predicted and X-ray conformations of these peptides gave RMSD values between 1.301 and 2.475 A. Binding conformations of 266 peptides with known binding affinities for HLA-A2 were then predicted using PePSSI. Structural analyses of these peptide-HLA-A2 conformations showed that peptide binding affinity is positively correlated with the number of peptide-MHC contacts and negatively correlated with the number of interfacial water molecules. These results are consistent with the relatively hydrophobic binding nature of the HLA-A2 peptide binding interface. In summary, PePSSI is capable of rapid and accurate prediction of peptide-MHC binding conformations, which may in turn allow estimation of MHCI-peptide binding affinity.

Structural basis for degenerate recognition of natural HIV peptide variants by cytotoxic lymphocytes

It is well established that even small changes in amino acid side chains of antigenic peptide bound to major histocompatibility complex (MHC) protein may completely abrogate recognition of the peptide-MHC (pMHC) complex by the T cell receptor (TCR). Often, however, several nonconservative substitutions in the peptide antigen are accommodated and do not impair its recognition by TCR. For example, a preponderance of natural sequence variants of the human immunodeficiency virus p17 Gag-derived peptide SLYNTVATL (SL9) are recognized by cytotoxic T lymphocytes, which implies that interactions with SL9 variants are degenerate both with respect to the class I MHC molecule and with respect to TCR. Here we study the molecular basis for this degenerate recognition of SL9 variants. We show that several SL9 variants bind comparably well to soluble HLA-A2 and to a particular soluble TCR and that these variants are active in the cognate cytotoxicity assay. Natural SL9 variation is restricted by its context in the HIV p17 matrix protein. High resolution crystal structures of seven selected SL9 variants bound to HLA-A2 all have remarkably similar peptide conformations and side-chain dispositions outside sites of substitution. This preservation of the peptide conformation despite epitope variations suggests a mechanism for the observed degeneracy in pMHC recognition by TCR and may contribute to the persistence of SL9-mediated immune responses in chronically infected individuals.

Dissection of the interaction of the human cytomegalovirus-derived US2 protein with major histocompatibility complex class I molecules: prominent role of a single arginine residue in human leukocyte antigen-A2

Human cytomegalovirus encodes several proteins that interfere with expression of major histocompatibility complex (MHC) class I molecules on the surface of infected cells. The unique short protein 2 (US2) binds to many MHC class I allomorphs in the endoplasmic reticulum, preventing cell surface expression of the class I molecule in question. The molecular interactions underlying US2 binding to MHC class I molecules and its allele specificity have not been fully clarified. In the present study, we first compared the sequences and the structures of US2 retained versus non-retained human leukocyte antigen (HLA) class I allomorphs to identify MHC residues of potential importance for US2 binding. On the basis of this analysis, 18 individual HLA-A2 mutants were generated and the ability of full-length US2 to bind wild-type and mutated HLA-A2 complexes was assessed. We demonstrate that Arg181 plays a critical role in US2-mediated inhibition of HLA-A2 cell surface expression. The structural comparison of all known crystal structures of HLA-A2 either alone, or in complex with T cell receptor or the CD8 co-receptor, indicates that binding of US2 to HLA-A2 results in a unique, large conformational change of the side chain of Arg181. However, although the presence of Arg181 seems to be a prerequisite for US2 binding to HLA-A2, it is not sufficient for binding to all MHC class I alleles.

Structural prediction of peptides bound to MHC class I

An ab initio structure prediction approach adapted to the peptide-major histocompatibility complex (MHC) class I system is presented. Based on structure comparisons of a large set of peptide-MHC class I complexes, a molecular dynamics protocol is proposed using simulated annealing (SA) cycles to sample the conformational space of the peptide in its fixed MHC environment. A set of 14 peptide-human leukocyte antigen (HLA) A0201 and 27 peptide-non-HLA A0201 complexes for which X-ray structures are available is used to test the accuracy of the prediction method. For each complex, 1000 peptide conformers are obtained from the SA sampling. A graph theory clustering algorithm based on heavy atom root-mean-square deviation (RMSD) values is applied to the sampled conformers. The clusters are ranked using cluster size, mean effective or conformational free energies, with solvation free energies computed using Generalized Born MV 2 (GB-MV2) and Poisson-Boltzmann (PB) continuum models. The final conformation is chosen as the center of the best-ranked cluster. With conformational free energies, the overall prediction success is 83% using a 1.00 Angstroms crystal RMSD criterion for main-chain atoms, and 76% using a 1.50 Angstroms RMSD criterion for heavy atoms. The prediction success is even higher for the set of 14 peptide-HLA A0201 complexes: 100% of the peptides have main-chain RMSD values < or =1.00 Angstroms and 93% of the peptides have heavy atom RMSD values < or =1.50 Angstroms. This structure prediction method can be applied to complexes of natural or modified antigenic peptides in their MHC environment with the aim to perform rational structure-based optimizations of tumor vaccines.

Highly congested nondistorted diheteroarylnaphthalenes: model compounds for the investigation of intramolecular pi-stacking interactions

[structure: see text] The rigid, highly congested structure of 1,8-diacridylnaphthalenes has been studied in solution and in the solid state. The unique geometry of these compounds forces the acridyl rings to undergo face-to-face interactions while rendering T-shaped orientations and face-to-edge interactions impossible. Crystallographic analysis shows that splaying between the heteroaryl rings decreases while twisting between the cofacial rings increases as the acridyl nitrogens of the 1,8-diacridylnaphthalene framework are subsequently oxidized. The peri-acridyl rings are slightly splayed but remain perfectly planar in all cases. The significant decrease in splaying indicates enhanced pi-pi-attraction between the electron-rich acridyl N-oxide moieties, which is in agreement with recently reported symmetry-adapted perturbation theory calculations. The pi-stacking and the molecular geometry between the acridyl rings observed in the solid state have been confirmed through in-solution studies showing characteristic proton NMR upfield shifts and optical properties indicative of static intramolecular arene-arene interactions. Acridyl protons located directly above the adjacent aryl moiety as a consequence of twisting between the heteroaryl rings were identified by COSY NMR measurements and found to intrude into the pi-cloud and diamagnetic ring current of the neighboring acridine. Different shapes and strong red shifts of the fluorescence emission maxima of the diacridylnaphthalenes in comparison to parental acridyl monomers have been attributed to static excimer emission.

Increased immunogenicity of an anchor-modified tumor-associated antigen is due to the enhanced stability of the peptide/MHC complex: implications for vaccine design

The use of "anchor-fixed" altered peptide ligands is of considerable interest in the development of therapeutic vaccines for cancer and infectious diseases, but the mechanism by which successful altered peptide ligands elicit enhanced immunity is unclear. In this study, we have determined the crystallographic structure of a major tumor rejection Ag, gp100(209-217), in complex with the HLA-A*0201 (HLA-A2) molecule, as well as the structure of a modified version of the peptide which substitutes methionine for threonine at position 2 (T2M; gp100(209-2M)). The T2M-modified peptide, which is more immunogenic in vitro and in vivo, binds HLA-A2 with a approximately 9-fold greater affinity and has a approximately 7-fold slower dissociation rate at physiological temperature. Within the limit of the crystallographic data, the T2M substitution does not alter the structure of the peptide/HLA-A2 complex. Consistent with this finding, in peripheral blood from 95 human subjects, we were unable to identify higher frequencies of T cells specific for either the native or modified peptide. These data strongly support the conclusion that the greater immunogenicity of the gp100(209-2M) peptide is due to the enhanced stability of the peptide/MHC complex, validating the anchor-fixing approach for generating therapeutic vaccine candidates. Thermodynamic data suggest that the enhanced stability of the T2M-modified peptide/HLA-A2 complex is attributable to the increased hydrophobicity of the modified peptide, but the gain due to hydrophobicity is offset considerably by the loss of a hydrogen bond made by the native peptide to the HLA-A2 molecule. Our findings have broad implications for the optimization of current vaccine-design strategies.

Quantitative structure-activity relationships and the prediction of MHC supermotifs

The underlying assumption in quantitative structure-activity relationship (QSAR) methodology is that related chemical structures exhibit related biological activities. We review here two QSAR methods in terms of their applicability for human MHC supermotif definition. Supermotifs are motifs that characterise binding to more than one allele. Supermotif definition is the initial in silico step of epitope-based vaccine design. The first QSAR method we review here--the additive method--is based on the assumption that the binding affinity of a peptide depends on contributions from both amino acids and the interactions between them. The second method is a 3D-QSAR method: comparative molecular similarity indices analysis (CoMSIA). Both methods were applied to 771 peptides binding to 9 HLA alleles. Five of the alleles (A*0201, A*0202, A*0203, A*0206 and A*6802) belong to the HLA-A2 superfamily and the other four (A*0301, A*1101, A*3101 and A*6801) to the HLA-A3 superfamily. For each superfamily, supermotifs defined by the two QSAR methods agree closely and are supported by many experimental data.

Strategic Mutations in the Class I Major Histocompatibility Complex HLA-A2 Independently Affect Both Peptide Binding and T Cell Receptor Recognition

Mutational studies of T cell receptor (TCR) contact residues on the surface of the human class I major histocompatibility complex (MHC) molecule HLA-A2 have identified a "functional hot spot" that comprises Arg(65) and Lys(66) and is involved in recognition by most peptide-specific HLA-A2-restricted TCRs. Although there is a significant amount of functional data on the effects of mutations at these positions, there is comparatively little biochemical information that could illuminate their mode of action. Here, we have used a combination of fluorescence anisotropy, functional assays, and Biacore binding experiments to examine the effects of mutations at these positions on the peptide-MHC interaction and TCR recognition. The results indicate that mutations at both position 65 and position 66 influence peptide binding by HLA-A2 to various extents. In particular, mutations at position 66 result in significantly increased peptide dissociation rates. However, these effects are independent of their effects on TCR recognition, and the Arg(65)-Lys(66) region thus represents a true "hot spot" for TCR recognition. We also made the observation that in vitro T cell reactivity does not scale with the half-life of the peptide-MHC complex, as is often assumed. Finally, position 66 is implicated in the "dual recognition" of both peptide and TCR, emphasizing the multiple roles of the class I MHC peptide-binding domain.

A novel MHCp binding prediction model

Many statistical and molecular mechanics models have been developed and tested for major histocompatibility complex peptide (MHCp) binding predictions during the last decade. The statistical model prediction using pooled peptide sequence data and three-dimensional modeling prediction by molecular mechanics calculations have been assessed for efficiency and human leukocyte antigen diversity coverage. We describe a novel predictive model using information gleaned from 29 human MHCp crystal structures. The validation for the new model is performed using four different sets of data: (1) MHCp crystal structures, (2) peptides with known IC(50) binding values, (3) peptides tested positive by tetramer staining, (4) peptides with known binding information at the MHCBN database. The model produces high prediction efficiencies (average 60 %) with good sensitivity (approximately 50%-73%) and specificity (52%-58%) values. The average positive predictive value of the model is 89%, while the average negative predictive value is only 18%. The efficiency is very high in predicting binders and very low in predicting nonbinders. This model is superior to many existing methods because of its potential application to any given MHC allele whose sequence is clearly defined.

Multidimensional frameworks assembled from silver(I) coordination polymers containing flexible bis(thioquinolyl) ligands: role of the intra- and intermolecular aromatic stacking interactions

The two flexible multidentate ligands 1,3-bis(8-thioquinolyl)propane (C3TQ) and 1,4-bis(8-thioquinolyl)butane (C4TQ) were reacted with AgX (X = CF(3)SO(3)(-) or ClO(4)(-)) to give four new complexes: ([Ag(C3TQ)](ClO(4)))(n)() 1, ([Ag(C3TQ)](CF(3)SO(3)))(n)() 2, ([Ag(2)(C4TQ)(CF(3)SO(3))(CH(3)CN)](CF(3)SO(3)))(n)() 3, and ([Ag(C4TQ)](ClO(4)))(n)() 4. All complexes have been characterized by elemental analysis, IR, and (1)H NMR spectroscopy. Single-crystal X-ray analysis showed that chain structures form for all complexes in which the quinoline rings interact via various intra- (1) or intermolecular (2, 3, and 4) pi-pi aromatic stacking interactions, which in the latter cases results in multidimensional structures. Additional weak interactions, such as Ag.O and Ag.S contacts and C-H.O hydrogen bonding, are also present and help form stable, crystalline materials. It was found that the (CH(2))(n) spacers (n = 3 or 4) affect the orientation of the two terminal quinolyl rings, thereby significantly influencing the specific framework structure that forms. If the same ligand is used, on the other hand, then the different counteranions have the greatest effect on the final structure.

Analysis of the pi-pi stacking interactions between the aminoglycoside antibiotic kinase APH(3')-IIIa and its nucleotide ligands

A key contact in the active site of an aminoglycoside phosphotransferase enzyme (APH(3')-IIIa) is a pi-pi stacking interaction between Tyr42 and the adenine ring of bound nucleotides. We investigated the prevalence of similar Tyr-adenine contacts and found that many different protein systems employ Tyr residues in the recognition of the adenine ring. The geometry of these stacking interactions suggests that electrostatics play a role in the attraction between these aromatic systems. Kinetic and calorimetric experiments on wild-type and mutant forms of APH(3')-IIIa yielded further experimental evidence of the importance of electrostatics in the adenine binding region and suggested that the stacking interaction contributes approximately 2 kcal/mol of binding energy. This type of information concerning the forces that govern nucleotide binding in APH(3')-IIIa will facilitate inhibitor design strategies that target the nucleotide binding site of APH-type enzymes.

Dynamic characteristics of a peptide-binding groove of human HLA-A2 class I MHC molecules: normal mode analysis of the antigen peptide-class I MHC complex

Class I major histocompatibility complex (MHC) binds antigen peptides with various sequences. We performed a normal mode analysis of HLA-A2 MHC that binds three peptides with different affinity. HLA-A2 MHC has a peptide-binding groove composed of two alpha-helices (residue 49-84, residue 140-179). Some residues in the center of the groove showed an increase in fluctuations and some residue pairs between two helix groups showed a negative change in correlations by removing the antigen peptide. The extent of the fluctuation and correlation changes correlated well with the experimental ranking of the three peptides in binding affinity. Some definite anti-correlative motions were found between two helix groups in low frequency modes (<50 cm(-1)) by removing the antigen peptide. We propose that the above anti-correlative motions play an important role to bind the antigen peptide, especially in obtaining a "dynamic fit."

Physicochemical explanation of peptide binding to HLA-A*0201 major histocompatibility complex: a three-dimensional quantitative structure-activity relationship study

A three-dimensional quantitative structure-activity relationship method for the prediction of peptide binding affinities to the MHC class I molecule HLA-A*0201 was developed by applying the CoMSIA technique on a set of 266 peptides. To increase the self consistency of the initial CoMSIA model, the poorly predicted peptides were excluded from the training set in a stepwise manner and then included in the study as a test set. The final model, based on 236 peptides and considering the steric, electrostatic, hydrophobic, hydrogen bond donor, and hydrogen bond acceptor fields, had q2 = 0.683 and r2 = 0.891. The stability of this model was proven by cross-validations in two and five groups and by a bootstrap analysis of the non-cross-validated model. The residuals between the experimental pIC50 (-logIC50) values and those calculated by "leave-one-out" cross-validation were analyzed. According to the best model, 63.2% of the peptides were predicted with /residuals/ < or = 0.5 log unit; 29.3% with 1.0 < or = /residuals/ < 0.5; and 7.5% with /residuals/ > 1.0 log unit. The mean /residual/ value was 0.489. The coefficient contour maps identify the physicochemical property requirements at each position in the peptide molecule and suggest amino acid sequences for high-affinity binding to the HLA-A*0201 molecule.

A comparative molecular similarity index analysis (CoMSIA) study identifies an HLA-A2 binding supermotif

The 3D-QSAR CoMSIA technique was applied to a set of 458 peptides binding to the five most widespread HLA-A2-like alleles: A*0201, A*0202, A*0203, A*0206 and A*6802. Models comprising the main physicochemical properties (steric bulk, electron density, hydrophobicity and hydrogen-bond formation abilities) were obtained with acceptable predictivity (q2 ranged from 0.385 to 0.683). The use of coefficient contour maps allowed an A2-supermotif to be identified based on common favoured and disfavoured areas. The CoMSIA definition for the best HLA-A2 binder is as follows: hydrophobic aromatic amino acid at position 1; hydrophobic bulky side chains at positions 2, 6 and 9; non-hydrogen-bond-forming amino acids at position 3; small aliphatic hydrogen-bond donors at position 4; aliphatic amino acids at position 5; small aliphatic side chains at position 7; and small aliphatic hydrophilic and hydrogen-bond forming amino acids at position 8.

Rational design of peptide-based tumor vaccines

Administration of synthetic peptides derived from proteins uniquely or overexpressed in tumor cells (tumor-associated antigens) can elicit tumor-specific immune responses in vivo. This is because cytotoxic T lymphocytes can recognize and lyse tumor cells that display peptides derived from tumor-associated antigens (TAAs) in the context of class I major histocompatibility complex (MHC) molecules. TAA peptides, in contrast to peptides of viral origin, generally bind weakly to the MHC molecule. In many cases, this explains the poor magnitude of T cell response directed at the tumor in vivo. Improving MHC binding as a strategy to upregulate antigen recognition can convert low affinity TAA peptides into useful tools in clinical trial settings. High-resolution structures of class I MHC molecules reported over the past two decades provided the framework for designing peptides that can induce optimal T cell response. This review will discuss the basic and clinical aspects of modifying native TAA peptides as tumor vaccines.

Quantitative approaches to computational vaccinology

This article reviews the newly released JenPep database and two new powerful techniques for T-cell epitope prediction: (i) the additive method; and (ii) a 3D-Quantitative Structure Activity Relationships (3D-QSAR) method, based on Comparative Molecular Similarity Indices Analysis (CoMSIA). The JenPep database is a family of relational databases supporting the growing need of immunoinformaticians for quantitative data on peptide binding to major histocompatibility complexes and to the Transporters associated with Antigen Processing (TAP). It also contains an annotated list of T-cell epitopes. The database is available free via the Internet (http://www.jenner.ac.uk/JenPep). The additive prediction method is based on the assumption that the binding affinity of a peptide depends on the contributions from each amino acid as well as on the interactions between the adjacent and every second side-chain. In the 3D-QSAR approach, the influence of five physicochemical properties (steric bulk, electrostatic potential, local hydrophobicity, hydrogen-bond donor and hydrogen-bond acceptor abilities) on the affinity of peptides binding to MHC molecules were considered. Both methods were exemplified through their application to the well-studied problem of peptides binding to the human class I MHC molecule HLA-A*0201.

Minor structural changes in a mutated human melanoma antigen correspond to dramatically enhanced stimulation of a CD4+ tumor-infiltrating lymphocyte line

While most immunotherapies for cancer have focused on eliciting specific CD8+ cytotoxic T lymphocyte killing of tumor cells, a mounting body of evidence suggests that stimulation of anti-tumor CD4+ T cell help may be required for highly effective therapy. Several MHC class II-restricted tumor antigens that specifically activate such CD4+ helper T lymphocytes have now been identified, including one from a melanoma tumor that is caused by a single base-pair mutation in the glycolytic enzyme triosephosphate isomerase. This mutation results in the conversion of a threonine residue to isoleucine within the antigenic epitope, concomitant with a greater than five log-fold increase in stimulation of a CD4+ tumor-infiltrating lymphocyte line. Here, we present the crystal structures of HLA-DR1 in complex with both wild-type and mutant TPI peptide antigens, the first structures of tumor peptide antigen/MHC class II complexes recognized by CD4+ T cells to be reported. These structures show that very minor changes in the binding surface for T cell receptor correspond to the dramatic differences in T cell stimulation. Defining the structural basis by which CD4+ T cell help is invoked in an anti-tumor immune response will likely aid the design of more effective cancer immunotherapies.

Additive method for the prediction of protein-peptide binding affinity. Application to the MHC class I molecule HLA-A*0201

A method has been developed for prediction of binding affinities between proteins and peptides. We exemplify the method through its application to binding predictions of peptides with affinity to major histocompatibility complex class I molecule HLA-A*0201. The method is named "additive" because it is based on the assumption that the binding affinity of a peptide could be presented as a sum of the contributions of the amino acids at each position and the interactions between them. The amino acid contributions and the contributions of the interactions between adjacent side chains and every second side chain were derived using a partial least squares (PLS) statistical methodology using a training set of 420 experimental IC50 values. The predictive power of the method was assessed using rigorous cross-validation and using an independent test set of 89 peptides. The mean value of the residuals between the experimental and predicted pIC50 values was 0.508 for this test set. The additive method was implemented in a program for rapid T-cell epitope search. It is universal and can be applied to any peptide-protein interaction where binding data is known.

Toward the quantitative prediction of T-cell epitopes: coMFA and coMSIA studies of peptides with affinity for the class I MHC molecule HLA-A*0201

A set of 102 peptides with affinity for the class I MHC HLA-A0201 molecule was subjected to three-dimensional quantitative structure-affinity relationship (3D QSAR) studies using comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). A test set of 50 peptides was used to determine the predictive value of the models. The CoMFA models gave q(2) and r(2)pred below 0.5. The best CoMSIA model has q(2) = 0.542 and r(2)pred = 0.679, and includes hydrophobic, steric, and H-bond donor fields. The hydrophobic interactions play a dominant role in peptide-MHC molecule binding. CoMSIA coefficient contour maps were used to analyze the structural features of the peptides accounting for the affinity in terms of the three positively contributing physicochemical properties: local hydrophobicity, steric bulk and hydrogen-bond-donor ability.

Dendritic cell vaccination induces cross-reactive cytotoxic T lymphocytes specific for wild-type and natural variant human immunodeficiency virus type 1 epitopes in HLA-A*0201/Kb transgenic mice

Dendritic cells (DC) are highly efficient at inducing primary T cell responses. Consequently, DC are being investigated for their potential to prevent and/or treat human immunodeficiency virus type 1 (HIV-1) infection. In the current study, we examined the capacity of DC to elicit CD8+ cytotoxic T lymphocyte (CTL) reactivity against an HLA-A*0201-restricted HIV-1 reverse transcriptase (pol) epitope (residues 476-484) and two naturally occurring variants. Previous work demonstrated that the wild-type pol epitope is recognized by CTLs from HIV-1-infected individuals, whereas the variant pol epitopes are not, despite binding to HLA-A*0201. In agreement with these observations, parenteral administration of wild-type pol peptide induced HLA-A*0201-restricted CTL activity in A2Kb transgenic mice. In contrast, similar treatment with the two variant pol peptides failed to stimulate CTL reactivity, and this lack of immunogenicity correlated with reduced peptide:HLA-A*0201 complex stability. However, CTL responses were induced in A2Kb transgenic mice upon adoptive transfer of syngeneic bone marrow DC pulsed with the variant pol peptides. Furthermore, DC pulsed with the wild-type pol peptide elicited CTLs that cross-reacted with the variant pol epitopes. These results demonstrate that DC effectively expand the T cell repertoire of a given epitope to include cross-reactive T cell clonotypes. Accordingly, DC vaccination may aid in immune recognition of HIV-1 escape variants by broadening the T cell response.

High-resolution structure of HLA-A*0201 in complex with a tumour-specific antigenic peptide encoded by the MAGE-A4 gene

The heterotrimeric complex of the human major histocompatibity complex (MHC) molecule HLA-A*0201, beta2-microglobulin and the decameric peptide GVYDGREHTV derived from the melanoma antigen (MAGE-A4 protein has been determined by X-ray crystallography at 1.4 A resolution. MAGE-A4 belongs to a family of genes that are specifically expressed in a variety of tumours. MAGE-A4-derived peptides are presented by MHC molecules at the cell surface to cytotoxic T-lymphocytes. As the HLA-A*0201:MAGE-A4 complex occurs only on tumour cells, it is considered to be an appropriate target for immunotherapy. The structure presented here reveals potential epitopes specific to the complex and indicates which peptide residues could be recognised by T-cell receptors. In addition, as the structure could be refined anisotropically, it was possible to describe the movements of the bound peptide in more detail.

Class I major histocompatibility complex anchor substitutions alter the conformation of T cell receptor contacts

An immunogenic peptide (GP2) derived from HER-2/neu binds to HLA-A2.1 very poorly. Some altered-peptide ligands (APL) of GP2 have increased binding affinity and generate improved cytotoxic T lymphocyte recognition of GP2-presenting tumor cells, but most do not. Increases in binding affinity of single-substitution APL are not additive in double-substitution APL. A common first assumption about peptide binding to class I major histocompatibility complex is that each residue binds independently. In addition, immunologists interested in immunotherapy frequently assume that anchor substitutions do not affect T cell receptor contact residues. However, the crystal structures of two GP2 APL show that the central residues change position depending on the identity of the anchor residue(s). Thus, it is clear that subtle changes in the identity of anchor residues may have significant effects on the positions of the T cell receptor contact residues.

A structural difference limited to one residue of the antigenic peptide can profoundly alter the biological outcome of the TCR-peptide/MHC class I interaction

The vesicular stomatitis virus (VSV) octapeptide RGYVYQGL binds to H-2K(b) and triggers a cytotoxic T cell response in mice. A variant peptide, RGYVYEGL (E6) with a glutamic acid for glutamine replacement at position 6 of the VSV peptide, elicits a T cell response with features that are quite different from those elicited by the wild-type VSV peptide. The differences found in the nature of the T cells responding to the E6 peptide include changes in both the V beta elements and the sequences of the complementarity-determining region 3 loops of their TCRs. Further experiments found that the E6 peptide can act as an antagonist for VSV-specific T cell hybridomas. To determine whether these differences in V beta usage, complementarity-determining region 3 sequences, and the switch from agonism to antagonism are caused by a conformational change on the MHC, the peptide, or both, we determined the crystal structure of the variant E6 peptide bound to H-2K(b). This structure shows that the only significant structural difference between H-2K(b)/E6 and the previously determined H-2K(b)/VSV is limited to the side chain of position 6 of the peptide, with no differences in the MHC molecule. Thus, a minor conformational change in the peptide can profoundly alter the biological outcome of the TCR-peptide/MHC interaction.