Starting from the X-ray structure of a class I majorhistocompatibility complex (MHC)-encoded protein (HLA-B*2705), a naturallypresented self-nonapeptide and two synthetic analogues were simulated in thebinding groove of two human leukocyte antigen (HLA) alleles (B*2703 andB*2705) differing in a single amino acid residue. After 200 ps moleculardynamics simulations of the solvated HLA-peptide pairs, some molecularproperties of the complexes (distances between ligand and protein center ofmasses, atomic fluctuations, buried versus accessible surface areas,hydrogen-bond frequencies) allow a clear discrimination of potent from weakMHC binders. The binding specificity of the three nonapeptides for the twoHLA alleles could be explained by the disruption of one hydrogen-bondingnetwork in the binding pocket of the HLA-B*2705 protein where the singlemutation occurs. Rearrangements of interactions in the B pocket, which bindsthe side chain of peptidic residue 2, and a weakening of interactionsinvolving the C-terminal end of the peptide also took place. In addition,extension of the peptide backbone using a β-Ala analogue did notabolish binding to any of the two HLA-B27 subtypes, but increased theselectivity for B*2703, as expected from the larger peptide binding groovein this subtype. A better understanding of the atomic details involved inpeptide selection by closely related HLA alleles is of crucial importancefor unraveling the molecular features linking particular HLA alleles toautoimmune diseases, and for the identification of antigenic peptidestriggering such pathologie
