Molecular mechanism of Be2+-ion binding to HLA-DP2: tetrahedral coordination, conformational changes and multi-ion binding.

The chemistry of beryllium is rather unusual, however, less explored as compared to other main group elements. This is mainly attributed to the high toxicity of beryllium, leading to chronic granulomatous pneumonitis, called chronic beryllium disease (CBD). It has been reported that Be2+-ion binding to the human leukocyte antigen protein (HLA-DP2) and peptide (M2) results in favorable interaction with the T-cell receptor protein (TCR), which initiates immune-mediated toxicity. We have carried out molecular dynamics (MD) simulations combined with quantum mechanical/molecular mechanical (QM/MM) studies to explore the binding nature of Be2+ with a HLA-DP2 protein and M2 peptide. The interaction between the negatively charged M2 peptide and the negatively charged binding cleft of HLA-DP2 is unfavorable. However, this interaction is stabilized by one Be2+ and two Na+-ions bridged by negatively charged carboxyl groups of glutamate residues (ß26E and ß69E) of the ß-chain of HLA-DP2 and one glutamate (p7E) and one aspartate residue (p4D) of the M2 peptide. This multi-ion cavity consists of tetrahedrally coordinated static Be2+ and Na+-ions, as well as one dynamically exchangeable Na+-ion. The smaller size and higher charge of the Be2+-ion as compared to the Na+-ion reduce the distance between the M2 peptide and the ß-chain of HLA-DP2, which results in conformational change suitable for TCR binding. However, the replacement of the Be2+ by the Na+-ion could not generate a suitable binding site for TCR.

Patterns of non-ARD variation in more than 300 full-length HLA-DPB1 alleles.

Our understanding of sequence variation in the HLA-DPB1 gene is largely restricted to the hypervariable antigen recognition domain (ARD) encoded by exon 2. Here, we employed a redundant sequencing strategy combining long-read and short-read data to accurately phase and characterise in full length the majority of common and well-documented (CWD) DPB1 alleles as well as alleles with an observed frequency of at least 0.0006% in our predominantly European sample set. We generated 664 DPB1 sequences, comprising 279 distinct allelic variants. This allows us to present the, to date, most comprehensive analysis of the nature and extent of DPB1 sequence variation. The full-length sequence analysis revealed the existence of two highly diverged allele clades. These clades correlate with the rs9277534 A → G variant, a known expression marker located in the 3'-UTR. The two clades are fully differentiated by 174 fixed polymorphisms throughout a 3.6 kb stretch at the 3'-end of DPB1. The region upstream of this differentiation zone is characterised by increasingly shared variation between the clades. The low-expression A clade comprises 59% of the distinct allelic sequences including the three by far most frequent DPB1 alleles, DPB1*04:01, DPB1*02:01 and DPB1*04:02. Alleles in the A clade show reduced nucleotide diversity with an excess of rare variants when compared to the high-expression G clade. This pattern is consistent with a scenario of recent proliferation of A-clade alleles. The full-length characterisation of all but the most rare DPB1 alleles will benefit the application of NGS for DPB1 genotyping and provides a helpful framework for a deeper understanding of high- and low-expression alleles and their implications in the context of unrelated haematopoietic stem-cell transplantation.

Computational assessment of miRNA binding to low and high expression HLA-DPB1 allelic sequences.

Cell surface expression of HLA-DP is allele specific. SNP rs9277534 (A/G), located in the 3'UTR of the DPB1 gene, has been associated with either low (A) or high (G) expression of DP on the cell surface. Considering the role of miRNAs in the regulation of gene expression, we computationally identified the miRNAs of two BLCLs, PGF and COX, predicted to interact with their corresponding DPB1 transcripts, DPB1 * 04:01:01:01-low expression and DPB1 * 03:01:01:01-high expression. The identified target sequences are located primarily in intron 2 and the 3'UTR. We hypothesize that gene expression may be influenced first by nuclear pre-mRNA events involving intronic regions, followed by the usual 3'UTR-associated events in the cytoplasm. The low DP expression allele was found to interact in silico with a larger number of miRNAs than the high expression allele. This pattern holds when examining either the entire transcript unit or simply the polymorphic sites that differentiate the alleles. Interestingly, the rs9277534 A/G polymorphism appears to be in linkage disequilibrium with polymorphisms targeted by the identified miRNAs. The multiplicity of sites targeted by different miRNAs suggests that the expression of DPB1 may be a dynamic process, influenced by different miRNAs under different states of the cell.

Amino acid residues in five separate HLA genes can explain most of the known associations between the MHC and primary biliary cholangitis.

Primary Biliary Cholangitis (PBC) is a chronic autoimmune liver disease characterised by progressive destruction of intrahepatic bile ducts. The strongest genetic association is with HLA-DQA1*04:01, but at least three additional independent HLA haplotypes contribute to susceptibility. We used dense single nucleotide polymorphism (SNP) data in 2861 PBC cases and 8514 controls to impute classical HLA alleles and amino acid polymorphisms using state-of-the-art methodologies. We then demonstrated through stepwise regression that association in the HLA region can be largely explained by variation at five separate amino acid positions. Three-dimensional modelling of protein structures and calculation of electrostatic potentials for the implicated HLA alleles/amino acid substitutions demonstrated a correlation between the electrostatic potential of pocket P6 in HLA-DP molecules and the HLA-DPB1 alleles/amino acid substitutions conferring PBC susceptibility/protection, highlighting potential new avenues for future functional investigation.

Identification of the new HLA-DPB1 allele, DPB1*647:01, in an Italian patient with leukemia.

A novel HLA-DPB1 allele, named HLA-DPB1*647:01, identified in a leukemia patient.

Metal-specific CD4+ T-cell responses induced by beryllium exposure in HLA-DP2 transgenic mice.

Chronic beryllium disease (CBD) is a granulomatous lung disorder that is associated with the accumulation of beryllium (Be)-specific CD4(+) T cells into the lung. Genetic susceptibility is linked to HLA-DPB1 alleles that possess a glutamic acid at position 69 (ßGlu69), and HLA-DPB1*02:01 is the most prevalent ßGlu69-containing allele. Using HLA-DP2 transgenic (Tg) mice, we developed a model of CBD that replicates the major features of the human disease. Here we characterized the T-cell receptor (TCR) repertoire of Be-responsive CD4(+) T cells derived from the lungs of Be oxide-exposed HLA-DP2 Tg mice. The majority of Be-specific T-cell hybridomas expressed TCR Vß6, and a subset of these hybridomas expressed identical or nearly identical ß-chains that were paired with different α-chains. We delineated mimotopes that bind to HLA-DP2 and form a complex recognized by Be-specific CD4(+) T cells in the absence of Be. These Be-independent peptides possess an arginine at p5 and a tryptophan at p7 that surround the Be-binding site within the HLA-DP2 acidic pocket and likely induce charge and conformational changes that mimic those induced by the Be(2+) cation. Collectively, these data highlight the interplay between peptides and Be in the generation of an adaptive immune response in metal-induced hypersensitivity.

The impact of amino acid variability on alloreactivity defines a functional distance predictive of permissive HLA-DPB1 mismatches in hematopoietic stem cell transplantation.

A major challenge in unrelated hematopoietic stem cell transplantation (HSCT) is the prediction of permissive HLA mismatches, ie, those associated with lower clinical risks compared to their nonpermissive counterparts. For HLA-DPB1, a clinically prognostic model has been shown to be matching for T cell epitope (TCE) groups assigned by cross reactivity of T cells alloreactive to HLA-DPB1∗09:01; however, the molecular basis of this observation is not fully understood. Here, we have mutated amino acids (aa) in 10 positions of HLA-DPB1∗09:01 to other naturally occurring variants, expressed them by lentiviral vectors in B cell lines, and quantitatively measured allorecognition by 17 CD4(+) T cell effectors from 6 unrelated individuals. A significant impact on the median alloresponse was observed for peptide contact positions 9, 11, 35, 55, 69, 76, and 84, but not for positions 8, 56, and 57 pointing away from the groove. A score for the "functional distance" (FD) from HLA-DPB1∗09:01 was defined as the sum of the median impact of polymorphic aa in a given HLA-DPB1 allele on T cell alloreactivity. Established TCE group assignment of 23 alleles correlated with FD scores of ≤0.5, 0.6 to 1.9 and ≥2 for TCE groups 1, 2, and 3, respectively. Based on this, prediction of TCE group assignment will be possible for any given HLA-DPB1 allele, including currently 367 alleles encoding distinct proteins for which T cell cross reactivity patterns are unknown. Experimental confirmation of the in silico TCE group classification was successfully performed for 7 of 7 of these alleles. Our findings have practical implications for the applicability of TCE group matching in unrelated HSCT and provide new insights into the molecular mechanisms underlying this model. The innovative concept of FD opens new potential avenues for risk prediction in unrelated HSCT.

Chronic Beryllium Disease: revealing the role of beryllium ion and small peptides binding to HLA-DP2.

Chronic Beryllium (Be) Disease (CBD) is a granulomatous disorder that predominantly affects the lung. The CBD is caused by Be exposure of individuals carrying the HLA-DP2 protein of the major histocompatibility complex class II (MHCII). While the involvement of Be in the development of CBD is obvious and the binding site and the sequence of Be and peptide binding were recently experimentally revealed [1], the interplay between induced conformational changes and the changes of the peptide binding affinity in presence of Be were not investigated. Here we carry out in silico modeling and predict the Be binding to be within the acidic pocket (Glu26, Glu68 and Glu69) present on the HLA-DP2 protein in accordance with the experimental work [1]. In addition, the modeling indicates that the Be ion binds to the HLA-DP2 before the corresponding peptide is able to bind to it. Further analysis of the MD generated trajectories reveals that in the presence of the Be ion in the binding pocket of HLA-DP2, all the different types of peptides induce very similar conformational changes, but their binding affinities are quite different. Since these conformational changes are distinctly different from the changes caused by peptides normally found in the cell in the absence of Be, it can be speculated that CBD can be caused by any peptide in presence of Be ion. However, the affinities of peptides for Be loaded HLA-DP2 were found to depend of their amino acid composition and the peptides carrying acidic group at positions 4 and 7 are among the strongest binders. Thus, it is proposed that CBD is caused by the exposure of Be of an individual carrying the HLA-DP2*0201 allele and that the binding of Be to HLA-DP2 protein alters the conformational and ionization properties of HLA-DP2 such that the binding of a peptide triggers a wrong signaling cascade.

Structural basis for the specific recognition of the major antigenic peptide from the Japanese cedar pollen allergen Cry j 1 by HLA-DP5.

The major allergen, Cry j 1, was isolated from Japanese cedar Cryptomeria japonica (Cry j) pollen and was shown to react with immunoglobulin E antibodies in the sera from pollinosis patients. We previously reported that the frequency of HLA-DP5 was significantly higher in pollinosis patients and the immunodominant peptides from Cry j 1 bound to HLA-DP5 to activate Th2 cells. In the present study, we determined the crystal structure of the HLA-DP5 heterodimer in complex with a Cry j 1-derived nine-residue peptide, at 2.4Å resolution. The peptide-binding groove recognizes the minimal peptide with 10 hydrogen bonds, including those between the negatively charged P1 pocket and the Lys side chain at the first position in the peptide sequence. We confirmed that HLA-DP5 exhibits the same Cry j 1-binding mode in solution, through pull-down experiments using structure-based mutations of Cry j 1. We also identified the characteristic residues of HLA-DP5 that are responsible for the distinct properties of the groove, by comparing the structure of HLA-DP5 and the previously reported structures of HLA-DP2 in complexes with pDRA of the self-antigen. The comparison revealed that the HLA-DP5·pCry j 1 complex forms several hydrogen bond/salt bridge networks between the receptor and the antigen that were not observed in the HLA-DP2·pDRA complex. Evolutionary considerations have led us to conclude that HLA-DP5 and HLA-DP2 represent two major groups of the HLA-DP family, in which the properties of the P1 and P4 pockets have evolved and acquired the present ranges of epitope peptide-binding specificities.

Structural basis of chronic beryllium disease: linking allergic hypersensitivity and autoimmunity.

T-cell-mediated hypersensitivity to metal cations is common in humans. How the T cell antigen receptor (TCR) recognizes these cations bound to a major histocompatibility complex (MHC) protein and self-peptide is unknown. Individuals carrying the MHCII allele, HLA-DP2, are at risk for chronic beryllium disease (CBD), a debilitating inflammatory lung condition caused by the reaction of CD4 T cells to inhaled beryllium. Here, we show that the T cell ligand is created when a Be(2+) cation becomes buried in an HLA-DP2/peptide complex, where it is coordinated by both MHC and peptide acidic amino acids. Surprisingly, the TCR does not interact with the Be(2+) itself, but rather with surface changes induced by the firmly bound Be(2+) and an accompanying Na(+) cation. Thus, CBD, by creating a new antigen by indirectly modifying the structure of preexisting self MHC-peptide complex, lies on the border between allergic hypersensitivity and autoimmunity.

Description of the novel HLA-DPB1*137:01 allele found in an Italian subject.

In this report, we describe the identification and sequencing of a novel HLA-DPB1 allele, found in an Italian haematological patient. This allele is identical to DPB1*17:01 except for a single nucleotide substitution (GAC→GAG) at position 57, which changes the encoded amino acid from Asp to Glu.

A novel family of human leukocyte antigen class II receptors may have its origin in archaic human species.

HLA class II α and ß chains form receptors for antigen presentation to CD4(+) T cells. Numerous pairings of class II α and ß subunits from the wide range of haplotypes and isotypes may form, but most of these combinations, in particular those produced by isotype mixing, yielded mismatched dimers. It is unclear how selection of functional receptors is achieved. At the atomic level, it is not known which interactions of class II residues regulate selection of matched αß heterodimers and the evolutionary origin of matched isotype mixed dimer formation. In this study we investigated assembly of isotype-mixed HLA class II α and ß heterodimers. Assembly and carbohydrate maturation of various HLA-class II isotype-mixed α and ß subunits was dependent on the groove binding section of the invariant chain (Ii). By mutation of polymorphic DPß sequences, we identified two motifs, Lys-69 and GGPM-(84-87), that are engaged in Ii-dependent assembly of DPß with DRα. We identified five members of a family of DPß chains containing Lys-69 and GGPM 84-87, which assemble with DRα. The Lys/GGPM motif is present in the DPß sequence of the Neanderthal genome, and this ancient sequence is related to the human allele DPB1*0401. By site-directed mutagenesis, we inspected Neanderthal amino acid residues that differ from the DPB1*0401 allele and aimed to determine whether matched heterodimers are formed by assembly of DPß mutants with DRα. Because the *0401 allele is rare in the sub-Saharan population but frequent in the European population, it may have arisen in modern humans by admixture with Neanderthals in Europe.

Quantitative prediction of peptide binding to HLA-DP1 protein.

The exogenous proteins are processed by the host antigen-processing cells. Peptidic fragments of them are presented on the cell surface bound to the major hystocompatibility complex (MHC) molecules class II and recognized by the CD4+ T lymphocytes. The MHC binding is considered as the crucial prerequisite for T-cell recognition. Only peptides able to form stable complexes with the MHC proteins are recognized by the T-cells. These peptides are known as T-cell epitopes. All T-cell epitopes are MHC binders, but not all MHC binders are T-cell epitopes. The T-cell epitope prediction is one of the main priorities of immunoinformatics. In the present study, three chemometric techniques are combined to derive a model for in silico prediction of peptide binding to the human MHC class II protein HLA-DP1. The structures of a set of known peptide binders are described by amino acid z-descriptors. Data are processed by an iterative self-consisted algorithm using the method of partial least squares, and a quantitative matrix (QM) for peptide binding prediction to HLA-DP1 is derived. The QM is validated by two sets of proteins and showed an average accuracy of 86 percent.

Identification of beryllium-dependent peptides recognized by CD4+ T cells in chronic beryllium disease.

Chronic beryllium disease (CBD) is a granulomatous disorder characterized by an influx of beryllium (Be)-specific CD4⁺ T cells into the lung. The vast majority of these T cells recognize Be in an HLA-DP­restricted manner, and peptide is required for T cell recognition. However, the peptides that stimulate Be-specific T cells are unknown. Using positional scanning libraries and fibroblasts expressing HLA-DP2, the most prevalent HLA-DP molecule linked to disease, we identified mimotopes and endogenous self-peptides that bind to MHCII and Be, forming a complex recognized by pathogenic CD4⁺ T cells in CBD. These peptides possess aspartic and glutamic acid residues at p4 and p7, respectively, that surround the putative Be-binding site and cooperate with HLA-DP2 in Be coordination. Endogenous plexin A peptides and proteins, which share the core motif and are expressed in lung, also stimulate these TCRs. Be-loaded HLA-DP2­mimotope and HLA-DP2­plexin A4 tetramers detected high frequencies of CD4⁺ T cells specific for these ligands in all HLADP2+ CBD patients tested. Thus, our findings identify the first ligand for a CD4⁺ T cell involved in metal-induced hypersensitivity and suggest a unique role of these peptides in metal ion coordination and the generation of a common antigen specificity in CBD.

Structural basis of metal hypersensitivity.

Metal hypersensitivity is a common immune disorder. Human immune systems mount the allergic attacks on metal ions through skin contacts, lung inhalation and metal-containing artificial body implants. The consequences can be simple annoyances to life-threatening systemic illness. Allergic hyper-reactivities to nickel (Ni) and beryllium (Be) are the best-studied human metal hypersensitivities. Ni-contact dermatitis affects 10 % of the human population, whereas Be compounds are the culprits of chronic Be disease (CBD). αß T cells (T cells) play a crucial role in these hypersensitivity reactions. Metal ions work as haptens and bind to the surface of major histocompatibility complex (MHC) and peptide complex. This modifies the binding surface of MHC and triggers the immune response of T cells. Metal-specific αß T cell receptors (TCRs) are usually MHC restricted, especially MHC class II (MHCII) restricted. Numerous models have been proposed, yet the mechanisms and molecular basis of metal hypersensitivity remain elusive. Recently, we determined the crystal structures of the Ni and Be presenting human MHCII molecules, HLA-DR52c (DRA*0101, DRB3*0301) and HLA-DP2 (DPA1*0103, DPB1*0201). These structures revealed unusual features of MHCII molecules and shed light on how metal ions are recognized by T cells.

[Polymorphic analysis of Mhc-DPB1 gene exon 2 in Tibetan macaques (Macaca thibetana)].

Major histocompatibility complex (MHC) molecules play an important role in the susceptibility and/or resistance to many diseases. To gain an insight into the MHC background of the Tibetan macaques (Macaca thibetana), and thereby facilitate their protection and application in biomedical research, the second exon of the Mhc-DPB1 genes from 70 Tibetan macaques in Sichuan Province were characterized by PCR, cloning, sequencing, and statistical analysis. A total of 18 Mhc-DPB1 alleles were identified from Tibetan macaques, of which one (Math-DPB1*01:06N) was a pseudogene. Math-DPB1*06:01:01 (67.14%) was the most frequent allele in all the 18 alleles detected, followed by Math-DPB1* 01:03:01 (37.14%), Math-DPB1*09:02 (25.71%), and Math-DPB1*22:01 (15.71%). The alignment of putative amino acid sequences of the 18 Math-DPB1 alleles showed that 5 variable sites were species-specific to Tibetan macaques. A phylogenetic tree constructed using DPB1 alleles in difference species demonstrated that the alleles for Math-DPB1, Mamu-DPB1, and Mafa-DPB1 tended to mix together, rather than cluster into a separate branch in a species-specific fashion, and the Trans-species polymorphism was also observed in the phylogenetic tree. Selection analysis revealed that balancing selection may play an important role in maintaining the polymorphism of Math-DPB1 genes.

Chronic beryllium disease, HLA-DPB1, and the DP peptide binding groove.

Multiple epidemiologic studies demonstrate associations between chronic beryllium disease (CBD), beryllium sensitization (BeS), and HLA-DPB1 alleles with a glutamic acid residue at position 69 (E69). Results suggest that the less-frequent E69 variants (non-*0201/*0202 alleles) might be associated with greater risk of CBD. In this study, we sought to define specific E69-carrying alleles and their amino acid sequences in the DP peptide binding groove, as well as their relationship to CBD and BeS risk, using the largest case control study to date. We enrolled 502 BeS/CBD subjects and 653 beryllium-exposed controls from three beryllium industries who gave informed consent for participation. Non-Hispanic white cases and controls were frequency-matched by industry. HLA-DPB1 genotypes were determined using sequence-specific primer PCR. The E69 alleles were tested for association with disease individually and grouped by amino acid structure using logistic regression. The results show that CBD cases were more likely than controls to carry a non-*02 E69 allele than an *02 E69, with odds ratios (95% confidence interval) ranging from 3.1 (2.1-4.5) to 3.9 (2.6-5.9) (p < 0.0001). Polymorphic amino acids at positions 84 and 11 were associated with CBD: DD versus GG, 2.8 (1.8-4.6), p < 0.0001; GD versus GG, 2.1 (1.5-2.8), p < 0.0001; LL versus GG, 3.2 (1.8-5.6), p < 0.0001; GL versus GG, 2.8 (2.1-3.8), p < 0.0001. Similar results were found within the BeS group and CBD/BeS combined group. We conclude that the less frequent E69 alleles confer more risk for CBD than does *0201. Recent studies examining how the composition and structure of the binding pockets influence peptide binding in MHC genes, as well of studies showing the topology of the TCR to likely bind DPB1 preferentially, give plausible biological rationale for these findings.

HLA-DP2 binding prediction by molecular dynamics simulations.

Major histocompatibility complex (MHC) II proteins bind peptide fragments derived from pathogen antigens and present them at the cell surface for recognition by T cells. MHC proteins are divided into Class I and Class II. Human MHC Class II alleles are grouped into three loci: HLA-DP, HLA-DQ, and HLA-DR. They are involved in many autoimmune diseases. In contrast to HLA-DR and HLA-DQ proteins, the X-ray structure of the HLA-DP2 protein has been solved quite recently. In this study, we have used structure-based molecular dynamics simulation to derive a tool for rapid and accurate virtual screening for the prediction of HLA-DP2-peptide binding. A combinatorial library of 247 peptides was built using the "single amino acid substitution" approach and docked into the HLA-DP2 binding site. The complexes were simulated for 1 ns and the short range interaction energies (Lennard-Jones and Coulumb) were used as binding scores after normalization. The normalized values were collected into quantitative matrices (QMs) and their predictive abilities were validated on a large external test set. The validation shows that the best performing QM consisted of Lennard-Jones energies normalized over all positions for anchor residues only plus cross terms between anchor-residues.

Mutagenesis of beryllium-specific TCRs suggests an unusual binding topology for antigen recognition.

Unconventional Ags, such as metals, stimulate T cells in a very specific manner. To delineate the binding landscape for metal-specific T cell recognition, alanine screens were performed on a set of Be-specific TCRs derived from the lung of a chronic beryllium disease patient. These TCRs are HLA-DP2-restricted and express nearly identical TCR Vß5.1 chains coupled with different TCR α-chains. Site-specific mutagenesis of all amino acids comprising the CDRs of the TCRA and TCRB genes showed a dominant role for Vß5.1 residues in Be recognition, with little contribution from the TCR α-chain. Solvent-exposed residues along the α-helices of the HLA-DP2 α- and ß-chains were also mutated to alanine. Two ß-chain residues, located near the proposed Be binding site of HLA-DP2, played a dominant role in T cell recognition with no contribution from the HLA-DP2 α-chain. These findings suggest that Be-specific T cells recognize Ag using an unconventional binding topology, with the majority of interactions contributed by TCR Vß5.1 residues and the HLA-DP2 ß1-chain. Thus, unusual docking topologies are not exclusively used by autoreactive T cells, but also for the recognition of unconventional metal Ags, such as Be.