On the calculation of binding free energies using continuum methods: Application to MHC class I protein‐peptide interactions
- 1 June 1997
- journal article
- research article
- Published by Wiley in Protein Science
- Vol. 6 (6) , 1293-1301
- https://doi.org/10.1002/pro.5560060617
Abstract
This paper describes a methodology to calculate the binding free energy (ΔG) of a protein‐ligand complex using a continuum model of the solvent. A formal thermodynamic cycle is used to decompose the binding free energy into electrostatic and non‐electrostatic contributions. In this cycle, the reactants are discharged in water, associated as purely nonpolar entities, and the final complex is then recharged. The total electrostatic free energies of the protein, the ligand, and the complex in water are calculated with the finite difference Poisson‐Boltzmann (FDPB) method. The nonpolar (hydrophobic) binding free energy is calculated using a free energy‐surface area relationship, with a single alkane/water surface tension coefficient (γaw). The loss in backbone and side‐chain configurational entropy upon binding is estimated and added to the electrostatic and the nonpolar components of ΔG. The methodology is applied to the binding of the murine MHC class I protein H‐2Kb with three distinct peptides, and to the human MHC class I protein HLA‐A2 in complex with five different peptides. Despite significant differences in the amino acid sequences of the different peptides, the experimental binding free energy differences (ΔΔGexp) are quite small (b and HLA‐A2 complexes, respectively). For each protein, the calculations are successful in reproducing a fairly small range of values for ΔΔGcalc (b and HLA‐A2 are not reproduced. For all protein‐peptide complexes that were treated, it was found that electrostatic interactions oppose binding whereas nonpolar interactions drive complex formation. The two types of interactions appear to be correlated in that larger nonpolar contributions to binding are generally opposed by increased electrostatic contributions favoring dissociation. The factors that drive the binding of peptides to MHC proteins are discussed in light of our results.Keywords
This publication has 68 references indexed in Scilit:
- The interpretation of protein structures: Estimation of static accessibilityPublished by Elsevier ,2004
- Ranking potential binding peptides to MHC molecules by a computational threading approachJournal of Molecular Biology, 1995
- A Continuum Model for Protein–Protein Interactions: Application to the Docking ProblemJournal of Molecular Biology, 1995
- Free R value: a novel statistical quantity for assessing the accuracy of crystal structuresNature, 1992
- Molecular Dynamics Simulation on the Connection MachineMolecular Simulation, 1991
- Hydrogen bond stereochemistry in protein structure and functionJournal of Molecular Biology, 1990
- Polar hydrogen positions in proteins: Empirical energy placement and neutron diffraction comparisonProteins-Structure Function and Bioinformatics, 1988
- CHARMM: A program for macromolecular energy, minimization, and dynamics calculationsJournal of Computational Chemistry, 1983
- Calculation of the electric potential in the active site cleft due to α-helix dipolesJournal of Molecular Biology, 1982
- The protein data bank: A computer-based archival file for macromolecular structuresJournal of Molecular Biology, 1977