Binding free energy differences in a TCR-peptide-MHC complex induced by a peptide mutation: a simulation analysis.

Recognition by the T-cell receptor (TCR) of immunogenic peptides presented by class I major histocompatibility complexes (MHCs) is the determining event in the specific cellular immune response against virus-infected cells or tumor cells. It is of great interest, therefore, to elucidate the molecular principles upon which the selectivity of a TCR is based. These principles can in turn be used to design therapeutic approaches, such as peptide-based immunotherapies of cancer. In this study, free energy simulation methods are used to analyze the binding free energy difference of a particular TCR (A6) for a wild-type peptide (Tax) and a mutant peptide (Tax P6A), both presented in HLA A2. The computed free energy difference is 2.9 kcal/mol, in good agreement with the experimental value. This makes possible the use of the simulation results for obtaining an understanding of the origin of the free energy difference which was not available from the experimental results. A free energy component analysis makes possible the decomposition of the free energy difference between the binding of the wild-type and mutant peptide into its components. Of particular interest is the fact that better solvation of the mutant peptide when bound to the MHC molecule is an important contribution to the greater affinity of the TCR for the latter. The results make possible identification of the residues of the TCR which are important for the selectivity. This provides an understanding of the molecular principles that govern the recognition. The possibility of using free energy simulations in designing peptide derivatives for cancer immunotherapy is briefly discussed.

[1]  J. Fontecilla-Camps,et al.  The three‐dimensional structure of a T‐cell antigen receptor VαVβ heterodimer reveals a novel arrangement of the Vβ domain , 1997 .

[2]  Thomas Simonson,et al.  Free energy of particle insertion , 1993 .

[3]  Michael H. Abraham,et al.  Thermodynamics of solute transfer from water to hexadecane , 1990 .

[4]  M. Karplus,et al.  Active site dynamics in protein molecules: A stochastic boundary molecular‐dynamics approach , 1985, Biopolymers.

[5]  M Karplus,et al.  Modeling of the TCR-MHC-peptide complex. , 2000, Journal of molecular biology.

[6]  M. Karplus,et al.  Locally accessible conformations of proteins: Multiple molecular dynamics simulations of crambin , 1998, Protein science : a publication of the Protein Society.

[7]  P. A. Bash,et al.  Calculation of the relative change in binding free energy of a protein-inhibitor complex. , 1987, Science.

[8]  P. A. Peterson,et al.  Crystal structure of an H-2Kb-ovalbumin peptide complex reveals the interplay of primary and secondary anchor positions in the major histocompatibility complex binding groove. , 1995, Proceedings of the National Academy of Sciences of the United States of America.

[9]  W F van Gunsteren,et al.  Decomposition of the free energy of a system in terms of specific interactions. Implications for theoretical and experimental studies. , 1994, Journal of molecular biology.

[10]  P. Romero,et al.  Induction of Potent Antitumor CTL Responses by Recombinant Vaccinia Encoding a Melan-A Peptide Analogue , 2000, The Journal of Immunology.

[11]  B. Braden,et al.  Crystal Structure of the Vα Domain of a T Cell Antigen Receptor , 1995, Science.

[12]  M. Mezei,et al.  Studies on free energy calculations. II. A theoretical approach to molecular solvation , 1994 .

[13]  M Karplus,et al.  T Cell Recognition of Hapten , 1999, The Journal of Biological Chemistry.

[14]  A. Smolyar,et al.  Atomic structure of an αβ T cell receptor (TCR) heterodimer in complex with an anti‐TCR Fab fragment derived from a mitogenic antibody , 1998, The EMBO journal.

[15]  Alexander D. MacKerell,et al.  All-atom empirical potential for molecular modeling and dynamics studies of proteins. , 1998, The journal of physical chemistry. B.

[16]  J A McCammon,et al.  Theoretical calculations of relative affinities of binding. , 1991, Methods in enzymology.

[17]  D. Speiser,et al.  Optimal activation of tumor-reactive T cells by selected antigenic peptide analogues. , 1999, International immunology.

[18]  M Karplus,et al.  Simulation analysis of the stability mutant R96H of T4 lysozyme. , 1991, Biochemistry.

[19]  Thomas Simonson,et al.  Continuum Treatment of Long-Range Interactions in Free Energy Calculations. Application to Protein−Ligand Binding. , 1997 .

[20]  M. Karplus,et al.  Deformable stochastic boundaries in molecular dynamics , 1983 .

[21]  William L. Jorgensen,et al.  Efficient computation of absolute free energies of binding by computer simulations. Application to the methane dimer in water , 1988 .

[22]  M. Karplus,et al.  Molecular recognition in proteins. Simulation analysis of substrate binding by a tyrosyl-tRNA synthetase mutant. , 1994, Journal of molecular biology.

[23]  A T Brünger,et al.  Thermodynamics of protein-peptide interactions in the ribonuclease-S system studied by molecular dynamics and free energy calculations. , 1992, Biochemistry.

[24]  M Karplus,et al.  Free energy simulations: The meaning of the individual contributions from a component analysis , 1994, Proteins.

[25]  M. Karplus,et al.  Hidden thermodynamics of mutant proteins: a molecular dynamics analysis. , 1989, Science.

[26]  M. Karplus,et al.  CHARMM: A program for macromolecular energy, minimization, and dynamics calculations , 1983 .

[27]  M Karplus,et al.  Solvent effects on protein motion and protein effects on solvent motion. Dynamics of the active site region of lysozyme. , 1989, Journal of molecular biology.

[28]  B M Baker,et al.  Four A6-TCR/peptide/HLA-A2 structures that generate very different T cell signals are nearly identical. , 1999, Immunity.

[29]  P A Kollman,et al.  Absolute and relative binding free energy calculations of the interaction of biotin and its analogs with streptavidin using molecular dynamics/free energy perturbation approaches , 1993, Proteins.

[30]  L R Pease,et al.  Structural basis of plasticity in T cell receptor recognition of a self peptide-MHC antigen. , 1998, Science.

[31]  M. Billeter,et al.  MOLMOL: a program for display and analysis of macromolecular structures. , 1996, Journal of molecular graphics.

[32]  Partho Ghosh,et al.  Structure of the complex between human T-cell receptor, viral peptide and HLA-A2 , 1996, Nature.

[33]  David A. Pearlman,et al.  A Comparison of Alternative Approaches to Free Energy Calculations , 1994 .

[34]  D. Wiley,et al.  Conversion of a T cell antagonist into an agonist by repairing a defect in the TCR/peptide/MHC interface: implications for TCR signaling. , 2000, Immunity.

[35]  D. Wiley,et al.  Two human T cell receptors bind in a similar diagonal mode to the HLA-A2/Tax peptide complex using different TCR amino acids. , 1998, Immunity.

[36]  Thomas Simonson,et al.  Protein‐protein recognition: An experimental and computational study of the R89K mutation in Raf and its effect on Ras binding , 2008, Protein science : a publication of the Protein Society.

[37]  M Karplus,et al.  Specific amino acid recognition by aspartyl-tRNA synthetase studied by free energy simulations. , 1998, Journal of molecular biology.

[38]  Robyn L. Stanfield,et al.  An αβ T Cell Receptor Structure at 2.5 Å and Its Orientation in the TCR-MHC Complex , 1996, Science.

[39]  Stefan Boresch,et al.  THE ROLE OF BONDED TERMS IN FREE ENERGY SIMULATIONS : 1. THEORETICAL ANALYSIS , 1999 .

[40]  M Karplus,et al.  Polar hydrogen positions in proteins: Empirical energy placement and neutron diffraction comparison , 1988, Proteins.

[41]  J W Gratama,et al.  Grafting primary human T lymphocytes with cancer-specific chimeric single chain and two chain TCR , 2000, Gene Therapy.

[42]  Peter A. Kollman,et al.  Free energy calculations on protein stability: Thr-157 .fwdarw. Val-157 mutation of T4 lysozyme , 1989 .

[43]  K. Garcia,et al.  Alanine Scanning Mutagenesis of an αβ T Cell Receptor: Mapping the Energy of Antigen Recognition , 1998 .

[44]  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.

[45]  A. Fersht,et al.  Structural response to mutation at a protein-protein interface. , 1998, Journal of molecular biology.

[46]  S J Wodak,et al.  Contribution of the hydrophobic effect to protein stability: analysis based on simulations of the Ile-96----Ala mutation in barnase. , 1991, Proceedings of the National Academy of Sciences of the United States of America.

[47]  M Karplus,et al.  The meaning of component analysis: decomposition of the free energy in terms of specific interactions. , 1995, Journal of molecular biology.

[48]  R. Mariuzza,et al.  Crystal structure of the beta chain of a T cell antigen receptor. , 1995, Science.

[49]  G. Ogg,et al.  Ex Vivo Staining of Metastatic Lymph Nodes by Class I Major Histocompatibility Complex Tetramers Reveals High Numbers of Antigen-experienced Tumor-specific Cytolytic T Lymphocytes , 1998, The Journal of experimental medicine.

[50]  Peter A. Kollman,et al.  FREE ENERGY CALCULATIONS : APPLICATIONS TO CHEMICAL AND BIOCHEMICAL PHENOMENA , 1993 .