Grouping of Class I HLA Alleles Using Electrostatic Distribution Maps of the Peptide Binding Grooves

Human leukocyte antigen (HLA) molecules involved in immune function by binding to short peptides (8–20 residues) have different sequences in different individuals belonging to distinct ethnic population. Hence, the peptide-binding function of HLA alleles is specific. Class I HLA alleles (alternative forms of a gene) are associated with CD8+ T cells, and their allele-specific sequence information is available at the IMGT/HLA database. The available sequences are one-dimensional (1D), and the peptide-binding functional inference often requires 3-dimensional (3D) structural models of respective alleles. Hence, 3D structures were constructed for 1,000 class I HLA alleles (310 A, 570 B, and 120 C) using MODELLER (a comparative protein modeling program for modeling protein structures). The electrostatic distribution maps were generated for each modeled structure using Deep View (Swiss PDB Viewer Version 3.7). The 1,000 models were then grouped into different categories by visual inspection of their electrostatic distribution maps in the peptide binding grooves. The distribution of the models based on electrostatic distribution was 30% negative (300), 1% positive (12), 8% neutral (84), and 60% (604) mixed (random mixture of negative, positive, and neutral). This grouping provides insight toward the inference for functional overlap among HLA alleles.