ABSTRACT
Immunogenicity of gliadin peptides in celiac disease (CD) is majorly determined by the pattern of molecular interactions with HLA-DQ and T-cell receptors (TCR). Investigation of the interactions between immune-dominant gliadin peptides, DQ protein, and TCR are warranted to unravel the basis of immunogenicity and variability contributed by the genetic polymorphisms. Homology modeling of HLA and TCR done using Swiss Model and iTASSER, respectively. Molecular interactions of eight common deamidated immune-dominant gliadin with HLA-DQ allotypes and specific TCR gene pairs were evaluated. Docking of the three structures was performed with ClusPro2.0 and ProDiGY was used to predict binding energies. Effects of known allelic polymorphisms and reported susceptibility SNPs were predicted on protein-protein interactions. CD susceptible allele, HLA-DQ2.5 was shown to have considerable binding affinity to 33-mer gliadin (ΔG = - 13.9; Kd = 1.5E - 10) in the presence of TRAV26/TRBV7. Higher binding affinity was predicted (ΔG = - 14.3, Kd = 8.9E - 11) when TRBV28 was replaced with TRBV20 paired with TRAV4 suggesting its role in CD predisposition. SNP rs12722069 at HLA-DQ8 that codes Arg76α forms three H-bonds with Glu12 and two H-bonds with Asn13 of DQ2 restricted gliadin in the presence of TRAV8-3/TRBV6. None of the HLA-DQ polymorphisms was found to be in linkage disequilibrium with reported CD susceptibility markers. Haplotypic presentations of rs12722069-G, rs1130392-C, rs3188043-C and rs4193-A with CD reported SNPs were observed in sub-ethnic groups. Highly polymorphic sites of HLA alleles and TCR variable regions could be utilized for better risk prediction models in CD. Therapeutic strategies by identifying inhibitors or blockers targeting specific gliadin:HLA-DQ:TCR binding sites could be investigated.
