Molecular mechanism of the susceptibility difference between HLA-B*27:02/04/05 and HLA-B*27:06/09 to ankylosing spondylitis: substitution analysis, MD simulation, QSAR modelling, and in vitro assay.

The human leukocyte antigen HLA-B27 is directly involved in the disease pathogenesis of ankylosing spondylitis (AS). HLA-B27 has a high degree of genetic polymorphism, with 105 currently known subtypes; the presence of aspartic acid at residue 116 (Asp116) has been found to play an essential role in AS susceptibility. Here, we systematically investigated the molecular mechanism of the susceptibility difference between the AS-associated subtypes HLA-B*27:02/04/05 and AS-unassociated subtypes HLA-B*27:06/09 to AS at sequence, structure, energetic and dynamic levels. In total seven variable residues were identified among the five studied HLA-B27 subtypes, in which Asp116 can be largely stabilized by a spatially vicinal, positively charged His114 through a salt bridge, while five other variable residues seem to have only a marginal effect on AS susceptibility. We also employed a quantitative structure-activity relationship approach to model the statistical correlation between peptide structure and affinity to HLA-B*27:05, a genetic ancestor of all other HLA-B27 subtypes and associated strongly with AS. The built regression predictor was verified rigorously through both internal cross-validation and external blind validation, and was then employed to identify potential HLA-B*27:05 binders from >20,000 cartilage-derived self-peptides. Subsequently, the binding potency of the top five antigenic peptides to HLA-B*27:05 was assayed in vitro using a FACS-based MHC stabilization experiment. Consequently, two (QRVGSDEFK and LRGAGTNEK) out of the five peptides were determined to have high affinity (BL50 = 5.5 and 15.8 nM, respectively) and, as expected, both of them possess positively charged Lys at the C-terminus.