Prediction of peptide binding to a major histocompatibility complex class I molecule based on docking simulation

Binding between major histocompatibility complex (MHC) class I molecules and immunogenic epitopes is one of the most important processes for cell-mediated immunity. Consequently, computational prediction of amino acid sequences of MHC class I binding peptides from a given sequence may lead to important biomedical advances. In this study, an efficient structure-based method for predicting peptide binding to MHC class I molecules was developed, in which the binding free energy of the peptide was evaluated by two individual docking simulations. An original penalty function and restriction of degrees of freedom were determined by analysis of 361 published X-ray structures of the complex and were then introduced into the docking simulations. To validate the method, calculations using a 50-amino acid sequence as a prediction target were performed. In 27 calculations, the binding free energy of the known peptide was within the top 5 of 166 peptides generated from the 50-amino acid sequence. Finally, demonstrative calculations using a whole sequence of a protein as a prediction target were performed. These data clearly demonstrate high potential of this method for predicting peptide binding to MHC class I molecules.