RESUMEN
Cancer cell-killing by CD8+ T cells demands effective tumor antigen presentation by human leukocyte antigen class I (HLA-I) molecules. Screening and designing highly immunogenic neoantigens require quantitative computations to reliably predict HLA-peptide binding affinities. Here, with all-atom molecular dynamics (MD) simulations and free energy perturbation (FEP) methods, we design a collection of antigenic peptide candidates through in silico mutagenesis studies on immunogenic neoantigens, yielding enhanced binding affinities to HLA-B*44:02. In-depth structural dissection shows that introducing positively charged residues such as arginine to position 6 or lysine to position 7 of the candidates triggers conformational shifts in both peptides and the antigen-binding groove of the HLA, following the "induced-fit" mechanism. Enhancement in binding affinities compared to the wild-type was found in three out of five mutated candidates. The HLA pocket, capable of accommodating positively charged residues in positions from 5 to 7, is designated as the "dynamic pocket". Taken together, we showcase an effective structure-based binding affinity optimization framework for antigenic peptides of HLA-B*44:02 and underscore the importance of dynamic nature of the antigen-binding groove in concert with the anchoring motifs. This work provides structural insights for rational design of favorable HLA-peptide bindings and future developments in neoantigen-based therapeutics.