Premium
Predicting peptide binding to MHC pockets via molecular modeling, implicit solvation, and global optimization
Author(s) -
Schafroth Heather D.,
Floudas Christodoulos A.
Publication year - 2003
Publication title -
proteins: structure, function, and bioinformatics
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.699
H-Index - 191
eISSN - 1097-0134
pISSN - 0887-3585
DOI - 10.1002/prot.10608
Subject(s) - solvation , peptide , amino acid , solvent models , molecular dynamics , chemistry , binding energy , computational chemistry , molecule , poisson–boltzmann equation , implicit solvation , major histocompatibility complex , biochemistry , physics , organic chemistry , nuclear physics , gene , ion
Development of a computational prediction method based on molecular modeling, global optimization, and implicit solvation has produced accurate structure and relative binding affinity predictions for peptide amino acids binding to five pockets of the MHC molecule HLA‐DRB1*0101. Because peptide binding to MHC molecules is essential to many immune responses, development of such a method for understanding and predicting the forces that drive binding is crucial for pharmaceutical design and disease treatment. Underlying the development of this prediction method are two hypotheses. The first is that pockets formed by the peptide binding groove of MHC molecules are independent, separating the prediction of peptide amino acids that bind within individual pockets from those that bind between pockets. The second hypothesis is that the native state of a system composed of an amino acid bound to a protein pocket corresponds to the system's lowest free energy. The prediction method developed from these hypotheses uses atomistic‐level modeling, deterministic global optimization, and three methods of implicit solvation: solvent‐accessible area, solvent‐accessible volume, and Poisson–Boltzmann electrostatics. The method predicts relative binding affinities of peptide amino acids for pockets of HLA‐DRB1*0101 by determining computationally an amino acid's global minimum energy conformation. Prediction results from the method are in agreement with X‐ray crystallography data and experimental binding assays. Proteins 2004;54:000–000. © 2003 Wiley‐Liss, Inc.