MM–GBSA binding free energy decomposition and T cell receptor engineering

Recognition by the T‐cell receptor (TCR) of immunogenic peptides (p) presented by class I major histocompatibility complexes (MHC) is the key event in the immune response against virus infected cells or tumor cells. The major determinant of T cell activation is the affinity of the TCR for the peptide‐MHC complex, though kinetic parameters are also important. A study of the 2C TCR/SIYR/H‐2Kb system using a binding free energy decomposition (BFED) based on the MM–GBSA approach had been performed to assess the performance of the approach on this system. The results showed that the TCR‐p‐MHC BFED including entropic terms provides a detailed and reliable description of the energetics of the interaction (Zoete and Michielin, 2007 ). Based on these results, we have developed a new approach to design sequence modifications for a TCR recognizing the human leukocyte antigen (HLA)‐A2 restricted tumor epitope NY‐ESO‐1. NY‐ESO‐1 is a cancer testis antigen expressed not only in melanoma, but also on several other types of cancers. It has been observed at high frequencies in melanoma patients with unusually positive clinical outcome and, therefore, represents an interesting target for adoptive transfer with modified TCR. Sequence modifications of TCR potentially increasing the affinity for this epitope have been proposed and tested in vitro. T cells expressing some of the proposed TCR mutants showed better T cell functionality, with improved killing of peptide‐loaded T2 cells and better proliferative capacity compared to the wild type TCR expressing cells. These results open the door of rational TCR design for adoptive transfer cancer therapy. Copyright © 2010 John Wiley & Sons, Ltd.

[1]  Peter A. Kollman,et al.  Computational alanine scanning of the 1:1 human growth hormone–receptor complex , 2002, J. Comput. Chem..

[2]  S. Rosenberg,et al.  Adoptive cell transfer: a clinical path to effective cancer immunotherapy , 2008, Nature Reviews Cancer.

[3]  P. Romero,et al.  Vaccination with a Melan-A Peptide Selects an Oligoclonal T Cell Population with Increased Functional Avidity and Tumor Reactivity1 , 2002, The Journal of Immunology.

[4]  W. C. Still,et al.  Semianalytical treatment of solvation for molecular mechanics and dynamics , 1990 .

[5]  U. Aebi,et al.  Crystal structure of the human lamin A coil 2B dimer: implications for the head-to-tail association of nuclear lamins. , 2004, Journal of molecular biology.

[6]  M. Bonneville,et al.  Prevalent Role of TCR α-Chain in the Selection of the Preimmune Repertoire Specific for a Human Tumor-Associated Self-Antigen1 , 2003, The Journal of Immunology.

[7]  S. Jameson,et al.  The impact of duration versus extent of TCR occupancy on T cell activation: a revision of the kinetic proofreading model. , 2001, Immunity.

[8]  S. Tosatto,et al.  Application of MM/PBSA colony free energy to loop decoy discrimination: Toward correlation between energy and root mean square deviation , 2005, Protein science : a publication of the Protein Society.

[9]  P. Kollman,et al.  Calculating structures and free energies of complex molecules: combining molecular mechanics and continuum models. , 2000, Accounts of chemical research.

[10]  Pingping Guan,et al.  Computational vaccinology: quantitative approaches. , 2003, Novartis Foundation symposium.

[11]  J. Kirkwood Statistical Mechanics of Fluid Mixtures , 1935 .

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

[13]  Pedro Alexandrino Fernandes,et al.  Computational alanine scanning mutagenesis—An improved methodological approach , 2007, J. Comput. Chem..

[14]  N. Foloppe,et al.  Towards predictive ligand design with free-energy based computational methods? , 2006, Current medicinal chemistry.

[15]  James Andrew McCammon,et al.  Ligand-receptor interactions , 1984, Comput. Chem..

[16]  Tingjun Hou,et al.  Molecular dynamics and free energy studies on the wild-type and double mutant HIV-1 protease complexed with amprenavir and two amprenavir-related inhibitors: mechanism for binding and drug resistance. , 2007, Journal of medicinal chemistry.

[17]  Olivier Michielin,et al.  Binding free energy differences in a TCR-peptide-MHC complex induced by a peptide mutation: a simulation analysis. , 2002, Journal of molecular biology.

[18]  D. Flower,et al.  Modeling the peptide-T cell receptor interaction by the comparative molecular similarity indices analysis-soft independent modeling of class analogy technique. , 2006, Journal of medicinal chemistry.

[19]  Yi Li,et al.  Directed evolution of human T-cell receptors with picomolar affinities by phage display , 2005, Nature Biotechnology.

[20]  Yi Li,et al.  Directed evolution of human T cell receptor CDR2 residues by phage display dramatically enhances affinity for cognate peptide‐MHC without increasing apparent cross‐reactivity , 2006, Protein science : a publication of the Protein Society.

[21]  S E Ealick,et al.  Ab initio structure determination and refinement of a scorpion protein toxin. , 1997, Acta crystallographica. Section D, Biological crystallography.

[22]  R. Zwanzig High‐Temperature Equation of State by a Perturbation Method. I. Nonpolar Gases , 1954 .

[23]  G L Amidon,et al.  Solubility of nonelectrolytes in polar solvents. , 1972, Journal of pharmaceutical sciences.

[24]  S. Rosenberg,et al.  Cancer Regression in Patients After Transfer of Genetically Engineered Lymphocytes , 2006, Science.

[25]  Olivier Michielin,et al.  Structural assessment of single amino acid mutations: application to TP53 function , 2006, Human mutation.

[26]  D. Case,et al.  Insights into protein-protein binding by binding free energy calculation and free energy decomposition for the Ras-Raf and Ras-RalGDS complexes. , 2003, Journal of molecular biology.

[27]  Richard H. Henchman,et al.  Revisiting free energy calculations: a theoretical connection to MM/PBSA and direct calculation of the association free energy. , 2004, Biophysical journal.

[28]  Trent E Balius,et al.  Quantitative prediction of fold resistance for inhibitors of EGFR. , 2009, Biochemistry.

[29]  Zhiping Weng,et al.  Structure‐based design of a T‐cell receptor leads to nearly 100‐fold improvement in binding affinity for pepMHC , 2009, Proteins.

[30]  Bent K. Jakobsen,et al.  Single and Dual Amino Acid Substitutions in TCR CDRs Can Enhance Antigen-Specific T Cell Functions , 2008, The Journal of Immunology.

[31]  Shulin Zhuang,et al.  Insights into unbinding mechanisms upon two mutations investigated by molecular dynamics study of GSK3β–axin complex: Role of packing hydrophobic residues , 2007, Proteins.

[32]  David A Pearlman,et al.  Evaluating the molecular mechanics poisson-boltzmann surface area free energy method using a congeneric series of ligands to p38 MAP kinase. , 2005, Journal of medicinal chemistry.

[33]  Roland L. Dunbrack Rotamer libraries in the 21st century. , 2002, Current opinion in structural biology.

[34]  C. Brooks,et al.  Novel generalized Born methods , 2002 .

[35]  T. N. Bhat,et al.  The Protein Data Bank , 2000, Nucleic Acids Res..

[36]  Yi Li,et al.  High-Affinity TCRs Generated by Phage Display Provide CD4+ T Cells with the Ability to Recognize and Kill Tumor Cell Lines1 , 2007, The Journal of Immunology.

[37]  Matthew R. Clutter,et al.  High-affinity, peptide-specific T cell receptors can be generated by mutations in CDR1, CDR2 or CDR3. , 2005, Journal of molecular biology.

[38]  Irwin D Kuntz,et al.  Free energy calculations for theophylline binding to an RNA aptamer: Comparison of MM-PBSA and thermodynamic integration methods. , 2003, Biopolymers.

[39]  M Karplus,et al.  The contribution of vibrational entropy to molecular association. The dimerization of insulin. , 1994, Journal of molecular biology.

[40]  Samuel H. Yalkowsky,et al.  Solubility of nonelectrolytes in polar solvents. V. Estimation of the solubility of aliphatic monofunctional compounds in water using a molecular surface area approach , 1975 .

[41]  Pascal Bonnet,et al.  Scoring binding affinity of multiple ligands using implicit solvent and a single molecular dynamics trajectory: application to influenza neuraminidase. , 2005, Journal of molecular graphics & modelling.

[42]  B. Honig,et al.  Calculation of the total electrostatic energy of a macromolecular system: Solvation energies, binding energies, and conformational analysis , 1988, Proteins.

[43]  H. Meirovitch Recent developments in methodologies for calculating the entropy and free energy of biological systems by computer simulation. , 2007, Current opinion in structural biology.

[44]  Annick Dejaegere,et al.  Protein–protein recognition and interaction hot spots in an antigen–antibody complex: Free energy decomposition identifies “efficient amino acids” , 2007, Proteins.

[45]  Charles L. Brooks,et al.  New analytic approximation to the standard molecular volume definition and its application to generalized Born calculations , 2003, J. Comput. Chem..

[46]  Peter A. Kollman,et al.  A Ligand That Is Predicted to Bind Better to Avidin than Biotin: Insights from Computational Fluorine Scanning , 2000 .

[47]  Holger Gohlke,et al.  Change in protein flexibility upon complex formation: Analysis of Ras‐Raf using molecular dynamics and a molecular framework approach , 2004, Proteins.

[48]  Aurélien Grosdidier,et al.  Combined Simulation and Mutagenesis Analyses Reveal the Involvement of Key Residues for Peroxisome Proliferator-activated Receptorα Helix 12 Dynamic Behavior* , 2007, Journal of Biological Chemistry.

[49]  P A Kollman,et al.  Continuum solvent studies of the stability of RNA hairpin loops and helices. , 1998, Journal of biomolecular structure & dynamics.

[50]  Markus Meuwly,et al.  Importance of individual side chains for the stability of a protein fold: Computational alanine scanning of the insulin monomer , 2006, J. Comput. Chem..

[51]  P. Kollman,et al.  Computational Alanine Scanning To Probe Protein−Protein Interactions: A Novel Approach To Evaluate Binding Free Energies , 1999 .

[52]  Robyn L Stanfield,et al.  How TCRs bind MHCs, peptides, and coreceptors. , 2006, Annual review of immunology.

[53]  I. Moreira,et al.  Unraveling the importance of protein-protein interaction: application of a computational alanine-scanning mutagenesis to the study of the IgG1 streptococcal protein G (C2 fragment) complex. , 2006, The journal of physical chemistry. B.

[54]  A. Pohorille,et al.  Free energy calculations : theory and applications in chemistry and biology , 2007 .

[55]  Y. Chien,et al.  A TCR binds to antagonist ligands with lower affinities and faster dissociation rates than to agonists. , 1996, Immunity.

[56]  M. Gilson,et al.  Calculation of protein-ligand binding affinities. , 2007, Annual review of biophysics and biomolecular structure.

[57]  S. Valitutti,et al.  The efficiency of antigen recognition by CD8+ CTL clones is determined by the frequency of serial TCR engagement. , 1998, Journal of immunology.

[58]  M. Davis,et al.  Kinetic discrimination in T-cell activation. , 1996, Proceedings of the National Academy of Sciences of the United States of America.

[59]  Conrad C. Huang,et al.  UCSF Chimera—A visualization system for exploratory research and analysis , 2004, J. Comput. Chem..

[60]  Olivier Michielin,et al.  Comparison between computational alanine scanning and per‐residue binding free energy decomposition for protein–protein association using MM‐GBSA: Application to the TCR‐p‐MHC complex , 2007, Proteins.

[61]  Andrew Sewell,et al.  Structural and kinetic basis for heightened immunogenicity of T cell vaccines , 2005, The Journal of experimental medicine.

[62]  K D Wittrup,et al.  In vitro evolution of a T cell receptor with high affinity for peptide/MHC. , 2000, Proceedings of the National Academy of Sciences of the United States of America.

[63]  M K Gilson,et al.  Energetics of charge–charge interactions in proteins , 1988, Proteins.

[64]  P. Coveney,et al.  Molecular Basis of Peptide Recognition by the TCR: Affinity Differences Calculated Using Large Scale Computing1 , 2005, The Journal of Immunology.

[65]  David D Thomas,et al.  Thermodynamic and structural basis of phosphorylation-induced disorder-to-order transition in the regulatory light chain of smooth muscle myosin. , 2008, Journal of the American Chemical Society.

[66]  Markus Meuwly,et al.  Study of the insulin dimerization: Binding free energy calculations and per‐residue free energy decomposition , 2005, Proteins.

[67]  J. Ladbury Just add water! The effect of water on the specificity of protein-ligand binding sites and its potential application to drug design. , 1996, Chemistry & biology.

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

[69]  B. Kuhlman,et al.  Computational design of affinity and specificity at protein-protein interfaces. , 2009, Current opinion in structural biology.