Abstract: The present invention discloses a method for predicting peptides that are capable of binding to HLA molecules that incorporate the crystal structure of HLA molecules. An improved HLA-specific peptide docking workflow is used to simulate the occupancy of a peptide on the binding pocket of an HLA molecule, and three models are trained to predict the binding of the peptide to HLA molecules. The results show that these models predict HLA-allele specific binding peptides with extremely high accuracy.