Analysis of peptide-protein binding using amino acid descriptors: prediction and experimental verification for human histocompatibility complex HLA-A0201.

Amino acid descriptors are often used in quantitative structure-activity relationship (QSAR) analysis of proteins and peptides. In the present study, descriptors were used to characterize peptides binding to the human MHC allele HLA-A0201. Two sets of amino acid descriptors were chosen: 93 descriptors taken from the amino acid descriptor database AAindex and the z descriptors defined by Wold and Sandberg. Variable selection techniques (SIMCA, genetic algorithm, and GOLPE) were applied to remove redundant descriptors. Our results indicate that QSAR models generated using five z descriptors had the highest predictivity and explained variance (q2 between 0.6 and 0.7 and r2 between 0.6 and 0.9). Further to the QSAR analysis, 15 peptides were synthesized and tested using a T2 stabilization assay. All peptides bound to HLA-A0201 well, and four peptides were identified as high-affinity binders.

[1]  A. Sette,et al.  Cutting Edge: The Conversion of Arginine to Citrulline Allows for a High-Affinity Peptide Interaction with the Rheumatoid Arthritis-Associated HLA-DRB1*0401 MHC Class II Molecule1 , 2003, The Journal of Immunology.

[2]  A Sette,et al.  Improved induction of melanoma-reactive CTL with peptides from the melanoma antigen gp100 modified at HLA-A*0201-binding residues. , 1996, Journal of immunology.

[3]  Robert C. Glen,et al.  Novel Methods for the Prediction of logP, pKa, and logD , 2002, J. Chem. Inf. Comput. Sci..

[4]  J. Sidney,et al.  Prominent role of secondary anchor residues in peptide binding to HLA-A2.1 molecules , 1993, Cell.

[5]  Hongmao Sun,et al.  A Universal Molecular Descriptor System for Prediction of LogP, LogS, LogBB, and Absorption , 2004, J. Chem. Inf. Model..

[6]  G. Cruciani,et al.  Comparative molecular field analysis using GRID force-field and GOLPE variable selection methods in a study of inhibitors of glycogen phosphorylase b. , 1994, Journal of medicinal chemistry.

[7]  A Sette,et al.  Identification of a shared HLA-A*0201-restricted T-cell epitope from the melanoma antigen tyrosinase-related protein 2 (TRP2). , 1998, Cancer research.

[8]  S. Wold,et al.  Peptide QSAR on substance P analogues, enkephalins and bradykinins containing L- and D-amino acids. , 1990, Acta Chemica Scandinavica.

[9]  Lunli Yuan,et al.  Comparison of cytotoxic T lymphocyte responses induced by peptide or DNA immunization: Implications on immunogenicity and immunodominance , 1997, European journal of immunology.

[10]  S. Wold,et al.  The prediction of bradykinin potentiating potency of pentapeptides. An example of a peptide quantitative structure-activity relationship. , 1986, Acta chemica Scandinavica. Series B: Organic chemistry and biochemistry.

[11]  Xiaoyang Xia,et al.  High-throughput logP measurement using parallel liquid chromatography/ultraviolet/mass spectrometry and sample-pooling. , 2002, Rapid communications in mass spectrometry : RCM.

[12]  D. Flower,et al.  A comparative molecular similarity indices (CoMSIA) study of peptide binding to the HLA-A3 superfamily. , 2003, Bioorganic & medicinal chemistry.

[13]  Rebecca Harris,et al.  Genetic algorithms and self-organizing maps: a powerful combination for modeling complex QSAR and QSPR problems , 2004, J. Comput. Aided Mol. Des..

[14]  Irini Doytchinova,et al.  The HLA-A2-supermotif: a QSAR definition. , 2003, Organic & biomolecular chemistry.

[15]  P. Prathipati,et al.  2D-QSAR in hydroxamic acid derivatives as peptide deformylase inhibitors and antibacterial agents. , 2002, Bioorganic & medicinal chemistry.

[16]  M. Oldstone,et al.  Optimal lymphocytic choriomeningitis virus sequences restricted by H-2Db major histocompatibility complex class I molecules and presented to cytotoxic T lymphocytes , 1995, Journal of virology.

[17]  Irini A. Doytchinova,et al.  JenPep: A Novel Computational Information Resource for Immunobiology and Vaccinology , 2003, J. Chem. Inf. Comput. Sci..

[18]  D. Wiley,et al.  The antigenic identity of peptide-MHC complexes: A comparison of the conformations of five viral peptides presented by HLA-A2 , 1993, Cell.

[19]  Peter A. Hunt.  QSAR using 2D descriptors and TRIPOS' SIMCA , 1999, J. Comput. Aided Mol. Des..

[20]  K. Siebert,et al.  Quantitative structure-activity relationship modeling of peptide and protein behavior as a function of amino acid composition. , 2001, Journal of agricultural and food chemistry.

[21]  S. Park,et al.  Synthetic peptides of human papillomavirus type 18 E6 harboring HLA-A2.1 motif can induce peptide-specific cytotoxic T-cells from peripheral blood mononuclear cells of healthy donors. , 1998, Virus research.

[22]  J. Drijfhout,et al.  Detailed motifs for peptide binding to HLA-A*0201 derived from large random sets of peptides using a cellular binding assay. , 1995, Human immunology.

[23]  Ulf Norinder,et al.  Theoretical amino acid descriptors. Application to bradykinin potentiating peptides , 1991, Peptides.

[24]  Darren R Flower,et al.  Coupling In Silico and In Vitro Analysis of Peptide-MHC Binding: A Bioinformatic Approach Enabling Prediction of Superbinding Peptides and Anchorless Epitopes , 2004, The Journal of Immunology.

[25]  A Sette,et al.  Role of HLA-A motifs in identification of potential CTL epitopes in human papillomavirus type 16 E6 and E7 proteins. , 1994, Journal of immunology.

[26]  Channa K. Hattotuwagama,et al.  New horizons in mouse immunoinformatics: reliable in silico prediction of mouse class I histocompatibility major complex peptide binding affinity. , 2004, Organic & biomolecular chemistry.

[27]  M. Kanehisa,et al.  Cluster analysis of amino acid indices for prediction of protein structure and function. , 1988, Protein engineering.

[28]  J. Sacchettini,et al.  The three-dimensional structure of H-2Db at 2.4 Å resolution: Implications for antigen-determinant selection , 1994, Cell.

[29]  Pingping Guan,et al.  HLA-A3 supermotif defined by quantitative structure-activity relationship analysis. , 2003, Protein engineering.

[30]  N. Viner,et al.  T cell recognition of the dominant I-A(k)-restricted hen egg lysozyme epitope: critical role for asparagine deamidation. , 2001, The Journal of experimental medicine.

[31]  A Sette,et al.  Induction of tumor-reactive CTL from peripheral blood and tumor-infiltrating lymphocytes of melanoma patients by in vitro stimulation with an immunodominant peptide of the human melanoma antigen MART-1. , 1995, Journal of immunology.

[32]  M. Oldstone,et al.  Relative implication of peptide residues in binding to major histocompatibility complex class I H-2Db: application to the design of high-affinity, allele-specific peptides. , 1995, Molecular immunology.

[33]  M. Oldstone,et al.  Binding of Viral Antigens to Major Histocompatibility Complex Class I H-2Db Molecules Is Controlled by Dominant Negative Elements at Peptide Non-anchor Residues , 1996, The Journal of Biological Chemistry.

[34]  Irini A. Doytchinova,et al.  A Comparative Molecular Similarity Index Analysis (CoMSIA) study identifies an HLA-A2 binding supermotif , 2002, J. Comput. Aided Mol. Des..

[35]  O. Larsson,et al.  Carcinomas and Melanomas CTL Against Autologous and Allogeneic Peptide Epitopes That Can Elicit Specific Derived neu-Identification of New HER 2 / , 1999 .

[36]  L. Nádasdi,et al.  A study of the applicability of QSAR calculation for peptide hormones. , 1981, Biochemical and biophysical research communications.

[37]  L. Kier,et al.  Amino acid side chain parameters for correlation studies in biology and pharmacology. , 2009, International journal of peptide and protein research.

[38]  P. Sneath.  Relations between chemical structure and biological activity in peptides. , 1966, Journal of theoretical biology.

[39]  S. Wold,et al.  Peptide quantitative structure-activity relationships, a multivariate approach. , 1987, Journal of medicinal chemistry.

[40]  R. Holmdahl,et al.  Identification of the minimal glycopeptide core recognized by T cells in a model for rheumatoid arthritis. , 2005, Bioorganic & medicinal chemistry.

[41]  Hans-Georg Rammensee,et al.  MHC ligands and peptide motifs: first listing , 2004, Immunogenetics.

[42]  Erik Johansson,et al.  Multivariate design and modeling in QSAR , 1996 .

[43]  Bruce L. Bush,et al.  Sample-distance partial least squares: PLS optimized for many variables, with application to CoMFA , 1993, J. Comput. Aided Mol. Des..

[44]  S. Wold,et al.  Minimum analogue peptide sets (MAPS) for quantitative structure-activity relationships. , 2009, International journal of peptide and protein research.

[45]  R. Rees,et al.  Identification of peptide epitopes of MAGE-1, -2, -3 that demonstrate HLA-A3-specific binding , 1996, Cancer Immunology, Immunotherapy.

[46]  K Hasegawa,et al.  Partial Least Squares Modeling and Genetic Algorithm Optimization in Quantitative Structure-Activity Relationships , 2000, SAR and QSAR in environmental research.

[47]  O. Rötzschke,et al.  Consensus motifs and peptide ligands of MHC class I molecules. , 1993, Seminars in immunology.

[48]  S. Wold,et al.  New chemical descriptors relevant for the design of biologically active peptides. A multivariate characterization of 87 amino acids. , 1998, Journal of medicinal chemistry.

[49]  D. Flower,et al.  Additive method for the prediction of protein-peptide binding affinity. Application to the MHC class I molecule HLA-A*0201. , 2002, Journal of proteome research.

[50]  Channa K. Hattotuwagama,et al.  AntiJen: a quantitative immunology database integrating functional, thermodynamic, kinetic, biophysical, and cellular data , 2005, Immunome research.

[51]  T. Elliott,et al.  Phosphorylated peptides can be transported by TAP molecules, presented by class I MHC molecules, and recognized by phosphopeptide-specific CTL. , 1999, Journal of immunology.

[52]  Minoru Kanehisa,et al.  AAindex: Amino Acid index database , 2000, Nucleic Acids Res..

[53]  P. Cresswell,et al.  Impaired assembly and transport of HLA‐A and ‐B antigens in a mutant TxB cell hybrid. , 1986, The EMBO journal.

[54]  H. Scheraga,et al.  Statistical analysis of the physical properties of the 20 naturally occurring amino acids , 1985 .

[55]  G Hermanson,et al.  Binding of a peptide antigen to multiple HLA alleles allows definition of an A2-like supertype. , 1995, Journal of immunology.

[56]  H. Rammensee,et al.  Allele-specific motifs revealed by sequencing of self-peptides eluted from MHC molecules , 1991, Nature.

[57]  Xiaomin Luo,et al.  Molecular docking and 3D-QSAR studies on gag peptide analogue inhibitors interacting with human cyclophilin A. , 2002, Journal of medicinal chemistry.

[58]  Fumiyoshi Yamashita,et al.  Two‐ and three‐dimensional QSAR of carrier‐mediated transport of β‐lactam antibiotics in Caco‐2 cells , 2004 .

[59]  W. Dunn,et al.  Amino acid side chain descriptors for quantitative structure-activity relationship studies of peptide analogues. , 1995, Journal of medicinal chemistry.

[60]  J. Sidney,et al.  Identification of subdominant CTL epitopes of the GP100 melanoma-associated tumor antigen by primary in vitro immunization with peptide-pulsed dendritic cells. , 1997, Journal of immunology.

[61]  A. Vitiello,et al.  The relationship between class I binding affinity and immunogenicity of potential cytotoxic T cell epitopes. , 1994, Journal of immunology.

[62]  S Wold,et al.  Multivariate parametrization of coded and non-coded amino acids by thin layer chromatography. , 1989, Progress in clinical and biological research.