Functional and structural characteristics of HLA-B*13:01-mediated specific T cells reaction in dapsone-induced drug hypersensitivity
Journal of Biomedical Science volume 29, Article number: 58 (2022)
Severe cutaneous adverse drug reactions (SCARs) are a group of serious clinical conditions caused by immune reaction to certain drugs. The allelic variance of human leukocyte antigens of HLA-B*13:01 has been strongly associated with hypersensitivities induced by dapsone (DDS). T-cell receptor mediated activation of cytotoxic T lymphocytes (CTLs) has also been suggested to play an essential role in pathogenesis of SCARs. However, HLA-B*13:01-DDS-TCR immune synapse that plays role in drug-induced hypersensitivity syndrome (DIHS) associated T cells activation remains uncharacterized.
To investigate the molecular mechanisms for HLA-B*13:01 in the pathogenesis of Dapsone-induced drug hypersensitivity (DDS-DIHS), we performed crystallization and expanded drug-specific CTLs to analyze the pathological role of DDS-DIHS.
Results showed the crystal structure of HLA-B*13:01-beta-2-microglobulin (β2M) complex at 1.5 Å resolution and performed mutation assays demonstrating that I118 or I119, and R121 of HLA-B*13:01 were the key residues that mediate the binding of DDS. Subsequent single-cell TCR and RNA sequencing indicated that TCRs composed of paired TRAV12-3/TRBV28 clonotype with shared CDR3 region specifically recognize HLA-B*13:01-DDS complex to trigger inflammatory cytokines associated with DDS-DIHS.
Our study identified the novel p-i-HLA/TCR as the model of interaction between HLA-B*13:01, DDS and the clonotype-specific TCR in DDS-DIHS.
Severe cutaneous adverse drug reactions (SCARs), including drug-induced hypersensitivity syndrome or drug reaction with eosinophilia and systemic symptoms (DRESS), Stevens-Johnson syndrome (SJS) and toxic epidermal necrolysis (TEN), are a group of potentially life-threatening adverse drug reactions with different clinical presentations and pathogenesis [1, 2]. SCARs are delayed-type hypersensitivity reactions (i.e., drug-induced immune reactions mediated by T cells that occur after a minimum of 3–4 days of exposure to the causative drug). Naive T cells are primed upon initial treatment with the causative drug, following which a pool of memory T cells is gradually formed by re-stimulation during continued exposure. Accordingly, SCARs resolve when treatment with the causative drug is discontinued and develop more rapidly upon renewal of its administration [3, 4].
Individuals with predisposition to develop SCARs appear to possess specific human leukocyte antigen (HLA) alleles, which have been found to trigger the condition by mediating drug-induced T cell activation [5,6,7,8,9]. Drug-specific T-cell receptors (TCRs) have been observed in oligoclonal T cell populations isolated from skin lesions and blood samples of allopurinol-induced SCAR patients with HLA-B*58:01 allele  and carbamazepine-induced SJS/TEN patients with HLA-B*15:02 allele . Several theories, such as hapten (and pro-hapten), p-i (pharmacological interaction with immune receptors) [12, 13], peptide repertoire alteration  and altered TCR repertoire model, have been proposed to explain the HLA-associated drug hypersensitivity [15, 16]. More recently, a unique population of T cells, termed skin-resident memory T cells, have been identified as possible SCAR mediators [10, 16]. In addition, high-resolution crystal structures of self-peptides presented by abacavir-modified HLA-B*57:01 have been reported, providing valuable insight into HLA-associated drug hypersensitivity on molecular level [11, 17]. Besides, the structural bioinformatics approach, which has seen remarkable advance in recent years, is a promising way to computationally elucidate the structural basis, dynamics and energetic properties of HLA-drug association as well as its pathological role in SCARs.
The HLA-B*13:01 allele has been identified as an important risk allele for SCARs induced by a number of drugs, including dapsone (diaminodiphenyl sulfone, DDS) [18,19,20,21,22,23,24]. DDS-induced hypersensitivity syndrome (DDS-DIHS) is one of the major causes of mortality amongst leprosy patients possessing the HLA-B*13:01 allele (i.e., about 2% of leprosy patients with DDS therapy, mortality rate 12.5%)[25,26,27]. DDS has also been used to treat various dermatologic and non-dermatologic diseases, including Pneumocystis jirovecii pneumonia , dermatitis herpetiformis , immune thrombocytopaenia  and eosinophilic fasciitis . Previous studies found that DDS-specific cytotoxic T cells were significantly activated when co-cultured with HLA-B*13:01-expressing antigen presenting cells (APCs) in the presence of DDS [32, 33]. However, HLA-B*13:01-DDS-TCR immune synapse that plays role in DDS-DIHS associated T cells activation remains uncharacterized, which makes prevention and treatment efforts more difficult.
In this study, we used structural biology, bioinformatics and mutation assays to demonstrate that extra-protein crystallization of HLA-B*13:01 and DDS binds directly to HLA-B*13:01 and to determine key HLA-B*13:01 residues mediating the binding. These key residues also explain why DDS-DIHS is specifically associated with HLA-B*13:01 and not the highly homologous HLA-B*13:02. Moreover, we carried out RNA-sequencing and single-cell TCR sequencing to identify DDS-specific TCR clonotypes from the blister cells and peripheral blood mononuclear cells (PBMCs) of DDS-DIHS and DDS-SJS patients. Using in silico structure modeling, we show that DDS interacts directly both with HLA-B*13:01 and the DDS-specific TCR, indicating that DIHS-associated TCL activation is mediated through the p-i -mediated drug interaction with HLA and TCR (p-i HLA/TCR model). The results of this study provide new insights for prevention and treatment of DDS-DIHS, and highlight the p-i HLA/TCR model as one of possible pathomechanisms for HLA-associated drug hypersensitivities.
Materials and methods
Study design and patient recruitment
From 2012 to 2020, we collected twelve clinical information, as well as eleven blood cells samples of DRESS and one blister cells sample of SJS/TEN from patients diagnosed with DDS-induced hypersensitivity reactions (DDS-DIHS), at the Institute of Dermatology, Chinese Academy of Medical Sciences, Nanjing and the Chang Gung Memorial Hospitals, Taiwan. DRESS was diagnosed according to the RegiSCAR criteria [34, 35] and was characterized by extensive erythematous maculopapular rash with periorbital edema and desquamation. The diagnosis of DRESS using the criteria and scoring system of the RegiSCAR group include cutaneous involvement with typical skin eruptions (e.g., exfoliative dermatitis, generalized maculopapular exanthema), fever (≧ 38.5 °C), enlarged lymph nodes (≧ 2 sites, ≧ 1 cm), presence of atypical lymphocytes and eosinophilia, systemic involvement (e.g., liver, kidney, and lung), time of resolution, and the evaluation of other potential causes. The clinical course, dosage and duration of dapsone, systemic involvement, and mortality were analyzed. The drug causality of SJS and DIHS was determined by the algorithm of drug causality for epidermal necrolysis (ALDEN) and Naranjo algorithm following the guidelines for assessment of drug causality published by the RegiSCAR group [35,36,37,38]. Only the cases with probable or certain causal association with DDS (ALDEN score greater than 4 or Naranjo algorithm greater than 5 were recruited. All DDS-DIHS patients and one DDS-SJS patient were HLA-B*13:01-positive, and were assessed by two dermatologists through review of photographs, pathology tissue slides, and medical records. We also recruited ten DDS-tolerant control subjects who had been administered DDS for at least 3 months without any cutaneous adverse reactions. Among them, three were HLA-B*13:01-positive, and seven were HLA-B*13:01-negative leprosy patient (see Table 1). In addition, twelve healthy control individuals were enrolled. Approvals were obtained from the institutional review board (IRB No. 100-4657A3, 201601761B0, 201902171A3 and 202001645B0), and informed consent was obtained from each participant.
Chemicals, vectors and cell lines
DDS, the DDS metabolite N-acetyl DDS (NAD), sulfadiazine and trichloroethylene were purchased from Sigma-Aldrich (St Louis, Mo). EBV-transformed autologous B-cell lines (B-LCLs), and HLA class I-deficient lymphoblastoid B-cell line (Hmy2.C1R), which were used as antigen-presenting cells (APC), were obtained from ATCC (Maryland, USA). For mutation assays, full-length open-reading frames encoding HLA-B*13:01, HLA-B*13:02 and β2M cloned into VP64-GFP vector. The vectors encoding wild-type HLA-B*13:01 and HLA-B*13:02 were then used as templates for generation of single-site and two-site HLA-B*13:01 and HLA-B*13:02 mutant variants with site-directed mutagenesis (Table 2). The VP64-GFP vectors encoding wild-type or one of mutant HLA-B*13:01 or HLA-B*13:02 variants were then transfected along with VP64-GFP-β2M, pspax2 and pMD2.G to Hmy2. C1R cells (HLA-B-deficient cells) and sorted by flow cytometry with GFP fluorescent protein as a reporter. Stable expression of the C1R transmembrane HLA-B*13:01 (C1R-HLA-B*13:01) clone was determined using flow cytometry by determining protein expression with anti-HLA antibody (w6/32, eBioscience).
Expression and purification of recombinant proteins for crystallographic experiments
Gene sequences for HLA-B*13:01 and beta-2-microglobulin (GenBank ID: AAA59627.1 and AAA51811.1) were synthesized at Sangon Biotech (China). The genes were cloned into the pET21a vector with an C-terminal 6 × His tag and expressed in Escherichia coli BL21 (DE3) cells. The cells were grown in LB medium at 37 °C until the OD600 nm reached 0.6, after which the expression of the recombinant proteins was induced at 16 °C for 20 h with the final concentration of isopropyl-B-D-1-thiogalactopyranoside (IPTG) of 0.1 mM. The cells were then harvested by centrifugation at 12,000 × g for 20 min, resuspended in lysis buffer containing 6 M guanidine-HCl and ultrasonicated. Recombinant proteins (180 mg HLA-B*13:01 and 60 mg β2M) were renatured overnight at 4 °C in dialysis buffer (0.1 M Tris pH = 8.0, 2 mM EDTA, 400 mM l-arginine-HCl, 0.5 mM oxidized glutathione), following which the precipitate was removed by centrifugation. The clear supernatant was loaded onto a 5 mL nickel-nitrilotriacetic acid (Ni–NTA) resin gravity column (Qiagen) pre-equilibrated with binding buffer (20 mM Tris–HCl pH = 8.0, 150 mM NaCl). The column was successively washed with 50 mL binding buffer containing 20 mM imidazole and 100 mM imidazole, and target protein was eluted using binding buffer containing 500 mM imidazole. The target proteins were then concentrated and loaded onto a Superdex G200 size-exclusion chromatography (SEC) column (120 mL, GE Healthcare Life Sciences, USA) pre-equilibrated with SEC buffer (20 mM Tris–HCl pH 8.0, 150 mM NaCl and 2 mM DTT). As initial SEC results indicated that the majority of HLA-B*13:01 did not bind to β2M, we used different oligopeptides predicted by NetMHC 4.0 Server (http://www.cbs.dtu.dk/services/NetMHC/) to stabilize the HLA-B*13:01-β2M complex. After purification of HLA-B*13:01 and β2M, the predicted peptides were added respectively. The formation of HLA-B*13:01-β2M complex was verified by SEC combined with SDS-PAGE, and finally found that the peptide (RQDILDLWI) could promote the formation of the complex.
Initial crystallization screens for the HLA-B*13:01-β2M-peptide complex (concentrated to 5–7 µg/mL) were performed with commercial crystallization screening kits using the sitting-drop vapor diffusion method at 16 °C. Crystallization drops contained 0.5 μL of the protein complex solution mixed with 0.5 μL of reservoir solution. Diffraction quality crystals were grown in 0.1 M HEPES (pH 6.5–7.5), 20% (v/v) PEG 4000, and 0.2 M NaCl. All crystals were cryoprotected in the reservoir solution supplemented with 20–25% (v/v) glycerol and flash-frozen in liquid nitrogen.
Data collection and structure determination
X-ray diffraction data for the HLA-B*13:01-β2M-peptide complex crystal was collected at the beamline BL-17U1 of the Shanghai Synchrotron Radiation Facility (SSRF). All data were indexed and scaled using HKL2000 software. The structure of HLA-B*13:01-β2M-peptide complex was determined by molecular replacement with MOLREP in PHENIX  using the structure of HLA-B*44:03 in complex with Epstein-Barr virus BZLF1-derived peptide (PDB ID: 4JQX) as the search model. AutoBuild in PHENIX was used to automatically build the structure model. After several rounds of positional and B-factor refinement using Phenix. Refine with TLS parameters alternated by manual model revision in Coot , model quality was verified using the PROCHECK program (https://www.ebi.ac.uk/thornton-srv/software/PROCHECK). The quality of the final model was validated with MolProbity . Structures were analyzed with the PDBePISA (Protein Interfaces, Surfaces, and Assemblies) tool on Dali server (http://ekhidna2.biocenter.helsinki.fi/dali). Detailed information on data collection and refinement statistics is listed in the Table 2. All structural figures were prepared in PyMOL (http://www.pymol.org).
Surface plasmon resonance
Biacore T200 surface plasmon resonance (SPR) system (GE Healthcare, Piscataway, NJ) was used to analyze the interaction between HLA-B*13:01 and different drugs. The recombinant HLA-B*13:01 and β2-microglobulin (Chemicon International, Temecula, Calif) were incubated together in PBS buffer and then immobilized on the chips via amine coupling reaction. The drugs used for the analysis were dissolved in PBS supplemented with 5% dimethyl sulfoxide (DMSO and flowed through the solid phase. Responses of the interaction were reference subtracted and corrected with a standard curve for the DMSO effects. BIA Evaluation software Version 1.0 was used for data analysis. The Kd value were determined by three independent experiments.
HLA-B*13:01 fix and pulse and binding assay
We investigated the DDS/HLA-B*13:01 interaction based on cytotoxic T lymphocytes (CTLs) activity. Specific CTLs reacted even if the APCs were fixed by 0.25% paraformaldehyde, excluding that either processing or intracellular metabolism is involved. Pulsed C1R-HLA-B*13:01 were incubated with the DDS for 1 h followed by two washing steps, no stimulation of drug-specific T-cell clones were observed for DDS. A C1R-HLA-B*13:01 stable clone was generated by means of transfection of the full-length cDNA plasmid encoding HLA-B*13:01. The endogenous peptides in HLA-B*13:01 were removed by using cold mild citric acid. After neutralizing with culture medium, the cells were incubated with β2M (4 µg/mL), GolgiStop (1 µg/mL; BD Biosciences, Calif), and the synthetic peptides for 3 h at room temperature.
HLA-B*13:01 protein extraction, co-immunoprecipitation (Co-IP) and peptides identification by liquid chromatograph-mass spectrometer (LC–MS)
C1R and C1R-HLA-B*13:01 were cultured with DDS (50 µg/mL) for 6 h and collected cell pellet to add RIPA lysis buffer. The supernatants were collected by centrifuge at 10,000 rcf for 10 min and transferred to a 3 kD ultrafiltration tube and centrifuge at 12,000 rcf for 10 min at 4 °C, then collected less than 3 kD peptides. The protein samples were testing by HLA-ABC antibody immobilization and Co-IP assay and analyzed by Easy-nLC1000 and Q Exactive™ Hybrid Quadrupole-Orbitrap™ Mass Spectrometer (Thermo Fisher Scientific, USA). The raw MS files were analyzed and searched against protein database based on the species of the samples using PEAKS Studio. The parameters were set as follows: the protein modifications were carbamidomethylation (C) (fixed), oxidation (M) (variable), Acetyl (Protein N-term) (variable); the enzyme specificity was set to Unspecific; the maximum missed cleavages were set to 2; the precursor ion mass tolerance was set to 20 ppm, and MS/MS tolerance was 20 ppm. Only high confident identified peptides were chosen for downstream protein identification analysis.
Growth of DDS-specific T cells
Peripheral blood mononuclear cells (PBMCs) were isolated from whole blood samples using Ficoll-Paque (GE, Life Sciences, USA) density gradient centrifugation. The PBMCs (5 × 10e5 per well) of patients with DDS-induced SCAR and other control individuals were cultured in 96-well microplates in TexMACS medium (miltenyi Biotec, Germany) supplemented with 5% human AB serum  (Sigma-Aldrich, Darmastadt, Germany), penicillin–streptomycin (Gibco Invitrogen, USA), 200U IL-2 (R&D systems, USA), and 50 µg/mL DDS (Sigma-Aldrich, St. Louis, MO) for 1 week at 37 °C in 5% CO2. DDS specific T-cell lines (TCLs) were obtained by culturing the patients’ PBMCs with DDS (50 µg/mL) for 12–14 days, and the expanded T cells were then restimulated with irradiated (50 Gy) autologous B-LCLs and DDS for approximately 4 to 5 cycles. The T-cell clones were obtained by means of serial dilution. Ten CTL TCLs were sorted by using FACSAria (BD, Franklin Lakes, NJ). In addition, dimethyl sulfoxide (DMSO) was added to the medium as the solvent control, and phytohemagglutinin (i.e., PHA) at a concentration of 10 µg/mL was used as the positive control. Proliferative clonotypes of T cells were performed to obtain for T-cell receptors (TCR) study.
TCR repertoire analysis by high-throughput next generation sequencing and single-cell sequencing
Proliferative clonotypes of T cells from DDS-induced SCAR patients were performed for (1) normal TCR sequencing based approach, that RNA was purified from each clone using Qiagen RNeasy, and 5′RACE was performed using BD SmartRace reagents and protocol, using the universal 5′primer, and a 3′gene-specific primer for the TCR constant region. (2) For experiments using the 10 × Genomics platform, Human T Cell (PN-1000005) were used according to the manufacturer’s instructions in the Single Cell V(D)J Reagent Kits User Guide. Cell number and concentration was confirmed with TC20™ Automated Cell Counter. Approximately Number (5000, etc.) cells were subjected immediately onto the 10XGenomics Chromium Controller machine for Gel Beads-in-Emulsion (GEM) generation. Full-length V(D)J segments from either T cell transcripts are enriched from first-strand cDNA via PCR amplification with primers specific to either the TCR or Ig constant regions prior to library construction. Library quality and concentration were assessed using Agilent Bioanalyzer 2100. Libraries were run on the Hiseq X or Novaseq for Illumina PE150 sequencing. Post-processing and quality control were performed by Novogene using the 10X Cell Ranger package (v2.2.0, 10X Genomics). Reads were aligned to mm10 (or GRCh38, etc.) reference assembly (v2.2.0, 10X Genomics).
Generation of single chain TCRα/TCRβ expression constructs and cells
Retroviruses encoding TCR genes carrying human TCR constant regions had the format TCRα-F2A-TCRβ and were produced in HEK-293T cells by transient transfection of retroviral-based plasmids and their packaging vectors (psPAX2 and pMD2.G) using lipofectamine 3000 (Thermo, Life, USA) according to the manufacturer’s protocol. At 48 h after transfection, the virus was collected, filtered through a 0.45-µm syringe filter, and used for infection. The T cells from PBMC were spin-infected with viral supernatant supplemented with 10 µg/mL Polybrene at 2500 r.p.m. at 30 °C for 90 min. On day 3 post-infection, TCR expressing cells were sorted by flow cytometry to establish derivative cell lines as indicated.
Flow cytometry was carried out using distinct fluorochrome-conjugated mAbs that recognize human CD3, CD4, and CD8 (BD Biosciences), human TCRα, TCRβ (BioLegend). These mAbs were labeled with Alexa Fluor 488, phycoerythrin (PE), phycoerythrin–cyanine 5 (PC5), phycoerythrin-cyanin 7 (PC7) or allophycocyanine. The cells were examined by means of multicolor flow cytometry on the BD FACSVerse flow cytometer (BD Biosciences), and data were analyzed with the Flowjo™ V10 software (BD Biosciences).
T cell activation and transduction
To transduce primary human T cells, PBMC (2 × 10e6 cells/mL) were activated in 24-well plates coated with 1 µg/mL antibody to CD3 (clone OKT3, eBioscience) in the presence of 1 µg/mL soluble antibody to CD28 (clone CD28.2, eBioscience) and 300U/mL interleukin-2 (IL-2; Fisher Scientific) in TexMACS medium supplemented with 5% (v/v) human AB serum and antibiotics (penicillin–streptomycin). After 48 h of activation, the cells were spin-infected with viral supernatant supplemented with 10 µg/mL Polybrene for 90 min at 2500 rpm at 30 °C. After spin-infection, the retroviral supernatant was replaced with fresh T cell medium containing 300 U/mL IL-2 and 1 µg/mL anti-CD28. The transduced primary T cells were cultured for 48 h and then used for cytotoxicity assays.
PBMCs isolated from DDS-SCAR patients were incubated with DDS for 1 week. The DDS-specific T cells were then expanded by adding 50 µg/mL DDS, 30 nanogram/mL CD28, 30 nanogram/mL OKT3 (anti-CD3 antibody), and 200 U/mL IL-2 for 2 weeks. The expanded T cells were reset for 3 days, and CTLs were isolated using FACS with anti-human CD45, CD3, and CD8 antibodies as markers. The DDS-specific CTLs were detected by quantifying the IFN-γ release with ELISpot assay following the instructions of the human IFN-γ ELISpot kit (BD Biosciences, USA). To perform ELISpot analysis, 5 × 10e4 DDS-specific CTLs or TCR transfectants were incubated with 1 × 10e4 C1R, C1R-HLA-B*13:01 or EBV-B-LCL-13:01 APC at 37 °C for 24 h. Data was analyzed using the CTL ELISpot Reader (CTL, USA).
Real-time RT-PCR assay
Total RNA extracted from patients’ PBMCs using Qiagen RNA Extraction Kit and was reverse transcribed using the Superscript II reverse transcriptase enzyme (Life Technologies). RT-PCR was performed in 20 µL reactions containing 0.8 µM forward and reverse primers and 5 µL of extracted RNA on an ABI 7300 (Thermo Fisher) using the AgPath-ID One-Step RT-PCR kit (Thermo Fisher). Thermocycle conditions included initial denaturation at 50 °C and 95 °C (10 min each), followed by 40 cycles at 95 °C (15 s) and 60 °C (1 min). The TCR primers used in quantitative-PCR amplification were showed in Additional file 3: Table S1. To verify primer specificities, melting curve analyses and PCR product sequencing were performed. Ct value were calculated via 2-ΔΔct method after normalization to GAPDH.
Modeling TCR/peptide/HLA and DDS
The 3D structures of HLA-B*13:01 and a TCR clonotype were first constructed with homology modeling method Modeller independently . And next the complex structure was assembled by referring the existing crystal structure (PDB code: 4EUP) using structure alignment program TMalign . This complex was further refined by a protein side-chain repacking method CISRR  as well as energy minimization based on CHARMM22 force field . The binding modes of DDS with the HLA-B*1301 were generated by protein–ligand blind docking tool CB-Dock  and superposed onto the complex structure of HLA-B*13:01-TCR. Images were generated with PyMOL (PyMOL Molecular Graphics System, Version 1.2, Schrodinger LLC, New York, NY).
Determination of inflammatory cytokine and cytotoxin production levels
To quantify IFN-γ, IL-3, IL-13, IL-5, IL-21, granulysin and granzyme B production levels, B-LCL (HLA-B*13:01) were plated in 24-well plates (10e5 cells/well) in co-culture with CTL (5 × 10e5 cells/well) and stimulated with DDS (50 µg/mL) in 1 mL complete RPMI1640 media overnight. Afterward, the supernatant concentrations of the inflammatory cytokines and cytotoxins mentioned above were quantified with ELISA (R&D, USA).
Statistical analyses were performed using SPSS for Windows, version 21.0 (IBM, Armonk, NY) and GraphPad prism 8.0 (San Diego, CA). Heatmaps were created using the built-in R heatmap function in stats package, and the Circos plots were generated using the VDJtools software (MiLaboratory). The entire experiment was repeated thrice. P-values for ratio estimates were calculated using a two-sided test. For analysis of ELISpot results, statistical comparison between two variables was performed by two-tailed Student’s t-test. Differences were considered statistically significant at P-values of less than 0.05.
All sequencing data obtained in the study have been deposited in the NCBI sequence read archive (SRA) database with links to BioProject accession ID “PRJNA867422” (it will be uploaded after manuscript accepted). The atomic coordinates and structure factors of the HLA-B*13:01-β2M-peptide complex was deposited in the Protein Data Bank under the accession number 7ER5. All other data are available from the authors upon reasonable request.
Twelve HLA-B*13:01-positive DDS SCAR patients were enrolled in this study. The underlying conditions caused by treatment with DDS were mostly diagnosed as chronic inflammatory dermatoses, including exfoliative dermatitis, erythroderma and eosinophilia (Table 1). The average daily DDS dose was 90.0 ± 20.4 mg, and the average onset of cutaneous manifestations was 29.2 ± 15.7 days after exposure to DDS. In addition, ten DDS-tolerant individuals (three HLA-B*13:01-positive and seven HLA-B*13:01-negative) who had been administered DDS for at least 3 months without any cutaneous adverse drug reaction, and twelve healthy controls were also recruited. All study participants were of Chinese descent.
DDS-specific cytotoxic T lymphocytes (CTLs) are activated by DDS in HLA-B*13:01-dependent manner
We isolated PBMCs from all DDS-DIHS patients and tolerant individuals, performed in vitro T cell expansion in the presence of 50 µg/mL DDS. After approximately three to four cycles of co-culture, CD8+CTLs constituted the majority of the T cell population in samples originating from the DDS-DIHS patients, but were largely absent in the T cell populations of tolerant individuals (Additional file 2: Fig. S1A–C). Thus, DDS-specific CD8+CTL clones could be obtained from DDS-DIHS patients and maintained in vitro. We then verified the dependence of DDS-specific CTL activation on the presence of different HLA-B alleles and stimulation by DDS or other drugs. The activation was detected using IFN-γ ELISPOT assay in the manner dependent on DDS concentration (Fig. 1A, Additional file 2: Fig. S2A). Under the conditions of DDS stimulation, DDS-specific CTLs were activated when co-cultured with autologous HLA-B*13:01-expressing B-lymphoblastoid cell line (B-LCL), but not with allogeneic HLA-B*15:02- or HLA-B*13:02- expressing B-LCL (Fig. 1B). To further confirm HLA-B*13:01-specific activation of CTLs in DDS-DIHS, we performed HLA blocking assay using DDS-specific T-cells obtained at the end of the second cycle (14 days) of clonal expansion. The results indicated that CTL-mediated IFN-γ release could be blocked by anti-HLA class I antibodies but not by anti–HLA class II antibodies (Fig. 1C, Additional file 2: Fig. S2B). Moreover, the cytotoxicity effect of DDS-specific CTLs by flow cytometric with PI/annexin V was dependent on both DDS and HLA-B*13:01 (Fig. 1D), and was significantly intensified with the increase of the effector/target cell ratios at 20:1 (Fig. 1E). Furthermore, the IFN-γ release from CTLs of DDS-DIHS patients was observed in a DDS concentration-dependent manner, but not in sulfadiazine (SFD), the DDS metabolite N-acetyl DDS (NAD), trichloroethylene (TCE) or in solvent control (DMSO) (Fig. 1F). Further investigated the cross-reactivity of dapsone and other sulfa drugs in 6 dapsone-induced SCAR patients carried HLA-B*13:01, the result showed that dapsone-induced SCAR patients carried HLA-B*13:01 had a high cross-reactivity to sulfasalazine and sulfamethoxazole (sensitivity 66.7% and 83.3%, respectively), but not to probenecid, gliclazide and glimepiride (Additional file 2: Fig. S3A and B). This result was consistent with the clinal and genetic finding. We then evaluated antigen presentation in DDS-DIHS patients using APC-fixing assays. DDS-elicited activation of CTLs could not be completely abolished by fixing APCs with 0.25% paraformaldehyde (Fig. 1G), but was abolished by a pulsed procedure (Fig. 1H). Altogether, these results demonstrated that expanded CTLs derived from DDS-DIHS patients are activated by DDS in HLA-B*13:01-dependent manner.
HLA-B*13:01 binds DDS directly
To determine how HLA-B*13:01 and DDS interact, we first expressed recombinant HLA-B*13:01 protein in Escherichia coli and purified it for surface plasmon resonance experiments, and then with the β2-microglobulin form the complex (Additional file 2: Fig. S4A). β2-microglobulin could bound to HLA-B*13:01 protein based on IP assay (Additional file 2: Fig. S4B). The results showed that the recombinant HLA-B*13:01 directly binds DDS with an estimated binding affinity in the micromolar range (Kd = 68.54 µM) (Additional file 2: Fig. S4C). To further investigate whether the endogenous peptides were involved in the HLA-B*13:01-mediated presentation of DDS, HLA-B*13:01 was expressed in the C1R B-cells (HLA-B-deficient cells) by stable transfection. Afterwards, the lysates were collected from HLA-B*13:01-expressing C1R cells and wild-type C1R cells to identify HLA-B*13:01-binding peptides with liquid chromatography-mass spectrometry (LC–MS). Considering that HLA-B*13:01 belongs to HLA class I receptors that preferentially bind 9–12-mer peptides, we used ultrafiltration tubes with the cut-off molecular weight of 3 kDa to collect the endogenous peptides. A total of seven peptides were detected in samples from HLA-B*13:01-expressing C1R cells significantly high expression than wild-type C1R cells. Among these peptides, VGFIGAGQLAF and VGFIGAGQLAY were expressed highly in HLA-B*13:01-expressing C1R cells in the presence of DDS (Additional file 2: Fig. S5A). These two abundant endogenous peptides elicited IFN-γ release from HLA-B*13:01-expressing cells in the presence of DDS, but not wild-type C1R cells (Additional file 2: Fig. S5B). However, treatment with excessive amounts of either of these two peptides did not strengthen the T cell activation when compared to DDS stimulation alone (Additional file 2: Fig. S5B). These results reveal that HLA-B*13:01 directly binds DDS, and that the presentation of endogenous peptides is suggested to be dispensable for binding of DDS to HLA-B*13:01 itself.
Key residues in the DDS-binding groove of HLA-B*13:01
To gain deeper insight into the mechanisms by which HLA-B*13:01 binds DDS, we had tried to crystallize a complex of HLA-B*13:01-β2M with peptide and DDS, but most of the proteins precipitated after adding DDS and no crystals of the complex was obtained. But we determined the crystal structure of HLA-B*13:01-beta-2-microglobulin (β2M)-peptide complex at 1.50 Å resolution (Additional file 1: PDB ID 7YG3). The peptide which predicted by the NetMHC 4.0 server, was used to stabilize the structure of HLA-B*13:01-β2M and facilitate crystallization but not participate in DDS interaction. HLA molecules are heterodimers composed of α and β chains, and the main differences lie in the composition of peptide binding region. The peptide binding region of HLA class I consists of α1 and α2 domains of α chain, while this region of HLA class II consists of α1 domain of α chain and β1 domain of β chain. The HLA-B*13:01-β2M-peptide complex is typical of HLA class I complex (Fig. 2A). We further showed that DDS potentially binds to the F pocket of binding groove (termed site 1), which appeared to be stable due to a constant root-mean-square deviation (RMSDs) value (Fig. 2B). For in-depth analysis of binding of DDS to the site 1, we conducted molecular modeling using Accelrys Discovery Studio. The results revealed that DDS is buried in a hydrophobic pocket formed by 108Y, 109Y, 118I, 119I, 142Y, 147Y and 171W of HLA-B*13:01 (Fig. 2C).
HLA-B*13:01 shares 98% of sequence identity with HLA-B*13:02, and only three amino acid residues differ (118IB13:01/118TB13:02, 119IB13:01/119WB13:02, 121RB13:01/121TB13:02) (Fig. 3A). Nevertheless, DDS-DIHS is strongly associated with HLA-B*13:01 but not HLA-B*13:02, implying that the identity of these three residues plays an important role in DDS binding. Indeed, we found that the side chain of 119WB13:02 creates steric hindrance with DDS (Fig. 3B). Besides, the van der Waals interactions that 118TB13:02 forms with adjacent amino acids, including 119WB13:02, are stronger than those formed by 118IB13:01. Thus, 118TB13:02 and 119WB13:02 work together to obstruct DDS binding by changing the shape of the binding pocket (Figs. 2C, 3B). While there is no significant difference between conformations of 121RB13:01 and 121TB13:02, these two residues differently affect the electrostatic properties of the DDS-binding pocket, indicating that 121RB13:01 may also play an important role in DDS binding (Fig. 3B).
To verify the importance of the aforementioned residues for DDS binding, we generated C1R cell lines stably expressing single-site and double-site HLA-B*13:01 mutants (HLA-B*13:01I118T, HLA-B*13:01I119W, HLA-B*13:01R121T, HLA-B*13:01I118T&I119W, HLA-B*13:01I118T&R121T, and HLA-B*13:01I119W&R121T) and HLA-B*13:02 mutants (HLA-B*13:02T118I, HLA-B*13:02W119I, HLA-B*13:02T121R, HLA-B*13:02T118I&W119I, HLA-B*13:02T118I&T121R, and HLA-B*13:02 I119W&T121R), and used them to perform in vitro co-culture assays with DDS-specific CTLs expanded for two weeks. The CTL-mediated IFN-γ release was measured for each mutant variant by ELISpot assay (Fig. 3C–F). As expected, the levels of IFN-γ released from DDS-specific CTLs were significantly elevated upon co-culture with HLA-B*13:01-expressing APCs and stimulation with DDS (ninefold increase compared with co-culture with wild-type C1R cells, P < 0.01, Fig. 3C). Interestingly, CTL activation was not significantly altered when co-cultured with HLA-B*13:01I118T, HLA-B*13:01I119W, or HLA-B*13:01R121T (Fig. 3C), but was greatly decreased in co-cultures with the double-site mutants HLA-B*13:01I118T&R121T and HLA-B*13:01I119W&R121T (Fig. 3D). Correspondingly, the single-site mutants HLA-B*13:02T118I and HLA-B*13:02W119I did not noticeably change CTL activation (Fig. 3E), whereas the single-site mutant HLA-B*13:02T121R and the two-site mutants HLA-B*13:02T118I&W119I, HLA-B*13:02T118I&T121R, and HLA-B*13:02I119W&T121R significantly increased CTL activation (Fig. 3E, F). Altogether, these results demonstrate that I118 or I119 and R121 are key residues for DDS binding by HLA-B*13:01.
TCR sequencing confirms poly- and mono-clonality of DDS-specific T cells
TCRs are essential for recognition of HLA-drug complex and T cell activation. To analyze 5′ RACE-based TCR sequencing results in parallel with UMI-corrected data, we firstly characterized the distribution of TCRα variable (TRAV) and TCRβ variable (TRBV) genes across blister cells and active T cells of DDS-induced SJS, as well as T cells derived from DDS-DIHS patients (Fig. 4A, B), and CD8+ T cell populations were significantly overexpressed over time for 21 or 35 days (Additional file 2: Fig. S1A–C). The most abundant transcripts included TRAV12-3, TRBV20-1, and TRBV28, which were all strongly correlated with increases in the blister cells, active PBMCs, as well as DDS-expanded T cells from DDS-SCAR patients, but not with the PBMCs from twelve healthy controls (Fig. 4A, B). We then mapped the circos plot for TCR CDR3 assemblies, and found a high level of homology and shared sequence identity within TRA/TRB CDR3 regions of DDS-DIHS patients (Fig. 4C, D). The shared specific TCR CDR3 clonotypes were oligoclonally expanded after 21 days of DDS stimulation, but were absent or present at very low frequencies (< 0.001%) in PBMCs of twelve healthy controls (Fig. 4E, F). Although the clonotype frequency varied, DDS-DIHS-1, DDS-DIHS-2, and DDS-DIHS-3 patients shared the TRA CDR3 clonotypes ‘CAIGAGNNRKLIW’, ‘CGTLSSYNTDKLIF’, ‘CAASLQGGSEKLVF’, ‘CADLDTGRRALTF’, ‘CILRNYNQGG KLIF’, ‘CAARENYGQNFVF’ (Fig. 4E), and the TRB CDR3 clonotypes ‘CASSFSGTGYFNEQFF’, ‘CASSPQGSYEQYF’, ‘CASSVQNGELFF’, ‘CSVVGEAEAFF’, ‘CASSGGRFNEKLFF’ (Fig. 4F). In addition, the TRA CDR3 clonotypes ‘CATLDNYGQNFVF’ and ‘CAASRGSYIPTF’ (Fig. 4E), and the TRB CDR3 clonotype ‘CASSPAGANVLTF’ (Fig. 4F) were expanded in DDS-DIHS-4 patient and were also detected in blister cells of DDS-induced SJS.
To investigate the gene expression profile associated with the aforementioned TCR clonotypes, we performed RNA-Seq and analyzed the differentially expressed genes (DEGs) of DDS-expanded T cells from five DDS-DIHS patients in comparison to unexpanded T cells of five healthy controls (Fig. 4G). The expression levels of IFN-γ, IL-3, IL-13, IL-5, IL-21, granulysin and granzyme B, all of which have been associated with pathogenesis of SJS and DDS-DIHS [48,49,50,51,52], were elevated in the DDS-expanded T cells expressing the aforementioned TCR clonotypes (Fig. 4G), further suggesting that these TCR clonotypes are involved in the pathogenesis of DDS-DIHS.
We next performed single-cell TCR sequencing to determine TRAV-TRBV pairing in T cells isolated from one DDS-DIHS patient (Fig. 4H). The results indicated that 98% of DDS-expanded T cells expressed the paired TRAV12-3/TRBV28 with the TRA CDR3 “CAASRGSYIPTF” clonotype (Fig. 4H). In contrast, only up to 0.3% of T cells from five healthy controls expressed the paired TRAV12-3/TRBV28, and none expressed the TRA CDR3 “CAASRGSYIPTF” clonotype. This finding indicated that T cells expressing the TCR composed of paired TRAV12-3/TRBV28 with TRA CDR3 ‘CAASRGSYIPTF’ clonotype are responsible for T cell activation and hypersensitivity reactions in DDS-DIHS.
Determination of the DDS-specific immunodominant TCR V clonotype and its role in pathogenesis of DDS-DIHS
We next determined the levels of inflammatory cytokines in DDS- and B-LCL (HLA-B*13:01)-stimulated CTLs. The results showed that IFN-γ, IL-3, IL-13, IL-5, IL-21 and granzyme B levels in the DDS- and B-LCL (HLA-B*13:01)-stimulated CTLs were significantly higher than in the mock group (Fig. 5A–F). We then performed qPCR assay to examine the expression levels of TRAV and TRBV clonotypes of T cells isolated from DDS-DIHS patients. Compared to ten healthy controls and ten DDS-tolerant controls, TRAV12-3 and TRBV28 expression levels in T cells from ten DDS-DIHS patients were approximately twofold and ninefold higher on the days 21 and 35 of in vitro expansion under DDS stimulation, respectively (Fig. 5G, H). Moreover, TRAV12-3 and TRBV28 clonotypes were found in HLA-B*13:01- positive DDS-DIHS patients, but were absent in all three DDS-tolerant HLA-B*13:01-positive control subjects. Thus, the absence of the TRAV12-3 and TRBV28 clonotypes in DDS-tolerant control subjects might explain why DDS-DIHS does not always develop in all HLA-B*13:01-positive persons taking DDS.
To further verify the role of the TRAV12-3 and TRBV28 clonotypes in pathogenesis of DDS-DIHS, we overexpressed the paired TRAV12-3/TRBV28, TRAV13-1/TRBV5-1 and TRAV12-3/TRBV30 in CD3 + T cells sorted after 48 h of expansion of PBMCs isolated from DDS-tolerant controls lacking these TCR clonotypes (Additional file 2: Fig. S6). Afterward, we performed in vitro co-culture assay and found that only the co-cultures of the TRAV12-3/TRBV28 T cells with C1R-HLA-B*13:01 or B-LCL-HLA-B*13:01 produced IFN-γ upon DDS stimulation (Fig. 5I). Together, these results further confirm that the TCR with paired TRAV12-3/TRBV28 clonotype is required for formation of immune synapse with HLA-B*13:01-DDS complex and DDS-DIHS-associated T cell activation.
The modeling of the HLA-B*13:01-DDS-TCR immune synapse
According to the results described above, DDS interacts with HLA-B*13:01 directly and is required for the formation of immune synapse with the DDS-DIHS-associated TCR clonotype. However, how does DDS bind to HLA and mediate the formation of the immune synapse with TCR without a peptide remain unclear. To answer the first question, we carried out independent in silico structure modeling of HLA-B*13:01 and TCR, followed by blind docking of DDS and HLA-B*13:01 and assembly based on the existing crystal structure (PDB ID: 4EUP) (Fig. 6A). The obtained HLA-B*13:01-DDS-TCR complex shows that the terminal amino-group of DDS contacts the residue S95 from the H-CDR3 region of TCR, possibly forming an intermolecular hydrogen bond. Concurrently, TCR forms multiple salt bridges with HLA-B*13:01 (Fig. 6B). This implies that both HLA-B*13:01 and DDS contribute to interaction with TCR and that, despite being significantly smaller than a 9-mer peptide, DDS mimics antigenic peptide due to dual interaction with TCR and HLA-B*13:01. We then performed blind docking of NAD, SFD and TCE with HLA-B*13:01, respectively. Among the three molecules, SFD and TCE showed favorable binding modes similar to DDS (Fig. 6E, F), whereas NAD was unable to bind to the F pocket of HLA-B*13:01 (Fig. 6D). Superimposition of the 3D structure of NAD over DDS in the F pocket indicated that Trp171 of HLA-B*13:01 clashes with the N-acetyl moiety of NAD (Fig. 6D), explaining the lack of interaction. Conversely, SFD showed almost the same binding mode as DDS but exposed its hydrophobic aromatic terminal to the CDR3 of TCR (Fig. 6E). This would make the interaction between SFD and TCR energetically unfavorable because the related CDR3 region consists of multiple hydrophilic residues, such as serines. Furthermore, TCE was also too small to form interactions with TCR (Fig. 6F).
Drugs can potentially be recognized as foreign antigens and trigger cell-mediated adaptive immune responses such as SCARs [3, 12]. Although HLA is known as an important immune receptor involved in drug presentation, few studies have investigated the structural basis for interactions between HLA and SCAR-associated drugs. Recently reported the promiscuous immune response induced by HLA Class-II‒restricted T cells in DDS-DIHS patients , but the detailed interaction and mechanism of HLA-B*13:01/DDS-restricted CTL responses remain not fully understood. Furthermore, the previously study used homology modeling to explore the interaction between HLA-B*13:01 and DDS . Our finding demonstrate that DDS binds directly to both HLA-B*13:01 and the DDS-specific TCR, and validate key HLA-B*13:01 residues mediating DDS binding. Thus, our findings show that the HLA-B*13:01-DDS-TCR interaction preferred to the novel p-i HLA/TCR model, in which drug directly interact with HLA and TCR, and peptide does not directly participate in the HLA-drug interaction, as the classic p-i model only focused on HLA and drug interaction in antigen-processing pathway [54, 55]. It suggested that TCR also plays an important role and can directly participate in HLA-drug interaction without any peptide involvement.
According to data presented in this study, DDS is almost entirely buried within a hydrophobic pocket located in the peptide-binding groove of HLA-B*13:01. Such binding mode differs from the way by which abacavir binds to HLA-B*57:01 , but is comparable to the interaction between allopurinol/oxypurinol and HLA-B*58:01 [8, 57]. Furthermore, the structure of DDS appears to be critical for its interaction with HLA-B*13:01 and TCR, as structural differences displayed by the highly similar NAD and SFD are sufficient to hinder interaction with one or both proteins, thereby failing to trigger CTL activation. While HLA-B*13:01 has been previously known as a genetic marker for DDS-DIHS , our study also demonstrated that DDS-specific CTLs are characterized by TCRs of specific TRAV/TRBV and CDR3 clonotypes. Importantly, we found that these clonotypes can be detected in unstimulated PBMCs isolated directly from HLA-B*13:01-positive DDS-DIHS patients. This is in line with clinical observations and may explain why certain HLA-B*13:01 carriers are tolerant to DDS.
In addition, our in vitro co-culture assays with TCR transfectants further validated the pathogenic role of immunodominant TCR clonotype (i.e., TRAV12-3/TRBV28 pair and specific CDR3 region) in DDS-DIHS patients, and the CTLs expressing this TCR clonotype were demonstrated to produce high levels of inflammatory cytokines and cytotoxins in DDS- and HLA-B*13:01-dependent manner. Similarly, recent studies found shared TCR clonotypes in the blister cells  and in vitro expanded T cells  from PBMCs of carbamazepine-induced SJS/TEN patients. Apart from susceptible HLA alleles, the specific TCR clonotypes were deficient for drug metabolism/clearance that has also been linked to SCARs . Thus, the risk factors for DDS-induced hypersensitivity syndrome include HLA-B*13:01 allele, specific TCR clonotypes, and impaired DDS clearance.
In this study, we propose that the existence of both HLA-B*13:01 and T cells with specific TCR clonotypes are required to trigger immune responses in patients with DDS-induced hypersensitivity syndrome. The HLA-B*13:01-DDS-TCR immune synapse is formed according to the novel p-i HLA/TCR model, in which DDS directly interacts both with HLA-B*13:01 and the TCR with specific clonotype. These findings not only reveal the general molecular mechanism for HLA-restricted drug hypersensitivity reactions, but also provide a framework for prevention of DDS-DIHS and development of susceptibility gene diagnosis (detection of HLA-B*13:01 and TCR clonotype) and T cell immunotherapy (antibody therapy for TCR clonotype) for treatment of its symptoms.
Availability of data and materials
The datasets supporting the conclusions of this article are included in this published article (and its Additional files).
Algorithm of drug causality for epidermal necrolysis
Antigen presenting cell
Third complementarity-determining region
Cytotoxic T lymphocytes
Differentially expressed genes
Drug-induced hypersensitivity syndrome
Dapsone-induced hypersensitivity syndrome
Drug reaction with eosinophilia and systemic symptoms
Delayed-type hypersensitivity reactions
Enzyme-linked immunosorbent assays
Human leukocyte antigens
Peripheral blood mononuclear cells
Pharmacological interaction with immune receptors
Severe cutaneous adverse reactions
Single-wavelength anomalous dispersion
Shanghai Synchrotron Radiation Facility
Surface plasmon resonance
T cell receptor
Toxic epidermal necrolysis
- TRAV genes:
T cell receptor alpha variable genes
- TRBV genes:
T cell receptor beta variable genes
Duong TA, Valeyrie-Allanore L, Wolkenstein P, et al. Severe cutaneous adverse reactions to drugs. Lancet. 2017;390(10106):1996–2011.
Cacoub P, Musette P, Descamps V, et al. The DRESS syndrome: a literature review. Am J Med. 2011;124(7):588–97.
Pichler W, Yawalkar N, Schmid S, et al. Pathogenesis of drug-induced exanthems. Allergy. 2002;57(10):884–93.
Hetherington S, McGuirk S, Powell G, et al. Hypersensitivity reactions during therapy with the nucleoside reverse transcriptase inhibitor abacavir. Clin Ther. 2001;23(10):1603–14.
Chessman D, Kostenko L, Lethborg T, et al. Human leukocyte antigen class I-restricted activation of CD8+ T cells provides the immunogenetic basis of a systemic drug hypersensitivity. Immunity. 2008;28(6):822–32.
Chung WH, Hung SI, Hong HS, et al. Medical genetics: a marker for Stevens-Johnson syndrome. Nature. 2004;428(6982):486.
Hung SI, Chung WH, Liou LB, et al. HLA-B*5801 allele as a genetic marker for severe cutaneous adverse reactions caused by allopurinol. Proc Natl Acad Sci U S A. 2005;102(11):4134–9.
Lin CH, Chen JK, Ko TM, et al. Immunologic basis for allopurinol-induced severe cutaneous adverse reactions: HLA-B*58:01-restricted activation of drug-specific T cells and molecular interaction. J Allergy Clin Immunol. 2015;135(4):1063–5.
Wei CY, Chung WH, Huang HW, et al. Direct interaction between HLA-B and carbamazepine activates T cells in patients with Stevens-Johnson syndrome. J Allergy Clin Immunol. 2012;129(6):1562–9.
White KD, Chung WH, Hung SI, et al. Evolving models of the immunopathogenesis of T cell-mediated drug allergy: the role of host, pathogens, and drug response. J Allergy Clin Immunol. 2015;136(2):219–34.
Zhou P, Zhang S, Wang Y, et al. Structural modeling of HLA-B*1502/peptide/carbamazepine /T-cell receptor complex architecture: implication for the molecular mechanism of carbamazepine-induced Stevens-Johnson syndrome/toxic epidermal necrolysis. J Biomol Struct Dyn. 2016;34(8):1806–17.
Ostrov DA, Grant BJ, Pompeu YA, et al. Drug hypersensitivity caused by alteration of the MHC-presented self-peptide repertoire. Proc Natl Acad Sci U S A. 2012;109(25):9959–64.
Pichler WJ. Pharmacological interaction of drugs with antigen-specific immune receptors: the p-i concept. Curr Opin Allergy Clin Immunol. 2002;2(4):301–5.
Illing PT, Vivian JP, Dudek NL, et al. Immune self-reactivity triggered by drug-modified HLA-peptide repertoire. Nature. 2012;486(7404):554–8.
Chung WH, Wang CW, Dao RL. Severe cutaneous adverse drug reactions. J Dermatol. 2016;43(7):758–66.
Pan RY, Chu MT, Wang CW, et al. Identification of drug-specific public TCR driving severe cutaneous adverse reactions. Nat Commun. 2019;10(1):3569.
Pompeu YA, Stewart JD, Mallal S, et al. The structural basis of HLA-associated drug hypersensitivity syndromes. Immunol Rev. 2012;250(1):158–66.
Tangamornsuksan W, Lohitnavy M. Association between HLA-B*1301 and dapsone-induced cutaneous adverse drug reactions: a systematic review and meta-analysis. JAMA Dermatol. 2018;154(4):441–6.
Wang H, Yan L, Zhang G, et al. Association between HLA-B*1301 and dapsone-induced hypersensitivity reactions among leprosy patients in China. J Invest Dermatol. 2013;133(11):2642–4.
Zhang FR, Liu H, Irwanto A, et al. HLA-B*13:01 and the dapsone hypersensitivity syndrome. N Engl J Med. 2013;369(17):1620–8.
Zhao Q, Alhilali K, Alzahrani A, et al. Dapsone- and nitroso dapsone-specific activation of T cells from hypersensitive patients expressing the risk allele HLA-B*13:01. Allergy. 2019;74(8):1533–48.
Park HJ, Park JW, Kim SH, et al. The HLA-B*13:01 and the dapsone hypersensitivity syndrome in Korean and Asian populations: genotype- and meta-analyses. Expert Opin Drug Saf. 2020;19(10):1349–56.
Satapornpong P, Pratoomwun J, Rerknimitr P, et al. HLA-B*13:01 is a predictive marker of dapsone-induced severe cutaneous adverse reactions in Thai patients. Front Immunol. 2021;12: 661135.
Shukkoor AA, Thangavelu S, George NE, et al. Dapsone hypersensitivity syndrome (DHS): a detrimental effect of dapsone? A case report. Curr Drug Saf. 2019;14(1):37–9.
Pandey B, Shrestha K, Lewis J, et al. Mortality due to dapsone hypersensitivity syndrome complicating multi-drug therapy for leprosy in Nepal. Trop Doct. 2007;37(3):162–3.
Shen J, Liu M, Zhou M, et al. Causes of death among active leprosy patients in China. Int J Dermatol. 2011;50(1):57–60.
Liu H, Wang Z, Bao F, et al. Evaluation of prospective HLA-B*13:01 screening to prevent dapsone hypersensitivity syndrome in patients with leprosy. JAMA Dermatol. 2019;155(6):666–72.
Girard PM, Landman R, Gaudebout C, et al. Dapsone-pyrimethamine compared with aerosolized pentamidine as primary prophylaxis against Pneumocystis carinii pneumonia and toxoplasmosis in HIV infection. The PRIO Study Group. N Engl J Med. 1993;328(21):1514–20.
Wozel G. The story of sulfones in tropical medicine and dermatology. Int J Dermatol. 1989;28(1):17–21.
Damodar S, Viswabandya A, George B, et al. Dapsone for chronic idiopathic thrombocytopenic purpura in children and adults–a report on 90 patients. Eur J Haematol. 2005;75(4):328–31.
Servy A, Clerici T, Malines C, et al. Eosinophilic fasciitis: a rare skin sclerosis. Patholog Res Int. 2010;2011: 716935.
Chen WT, Wang CW, Lu CW, et al. The function of HLA-B*13:01 involved in the pathomechanism of dapsone-induced severe cutaneous adverse reactions. J Invest Dermatol. 2018;138(7):1546–54.
Zhao Q, Almutairi M, Tailor A, et al. HLA Class-IIRestricted CD8(+) T cells contribute to the promiscuous immune response in dapsone-hypersensitive patients. J Invest Dermatol. 2021;141:2412.
Bastuji-Garin S, Rzany B, Stern RS, et al. Clinical classification of cases of toxic epidermal necrolysis, Stevens-Johnson syndrome, and erythema multiforme. Arch Dermatol. 1993;129(1):92–6.
Kardaun SH, Sekula P, Valeyrie-Allanore L, et al. Drug reaction with eosinophilia and systemic symptoms (DRESS): an original multisystem adverse drug reaction. Results from the prospective RegiSCAR study. Br J Dermatol. 2013;169(5):1071–80.
Naranjo CA, Busto U, Sellers EM, et al. A method for estimating the probability of adverse drug reactions. Clin Pharmacol Ther. 1981;30(2):239–45.
Sassolas B, Haddad C, Mockenhaupt M, et al. ALDEN, an algorithm for assessment of drug causality in Stevens-Johnson Syndrome and toxic epidermal necrolysis: comparison with case-control analysis. Clin Pharmacol Ther. 2010;88(1):60–8.
Wang CW, Yang LY, Chen CB, et al. Randomized, controlled trial of TNF-alpha antagonist in CTL-mediated severe cutaneous adverse reactions. J Clin Invest. 2018;128(3):985–96.
Liebschner D, Afonine PV, Baker ML, et al. Macromolecular structure determination using X-rays, neutrons and electrons: recent developments in Phenix. Acta Crystallogr D Struct Biol. 2019;75(Pt 10):861–77.
Vagin AA, Steiner RA, Lebedev AA, et al. REFMAC5 dictionary: organization of prior chemical knowledge and guidelines for its use. Acta Crystallogr D Biol Crystallogr. 2004;60(Pt 12 Pt 1):2184–95.
Chen VB, Arendall WB 3rd, Headd JJ, et al. MolProbity: all-atom structure validation for macromolecular crystallography. Acta Crystallogr D Biol Crystallogr. 2010;66(Pt 1):12–21.
Li G, Bethune MT, Wong S, et al. T cell antigen discovery via trogocytosis. Nat Methods. 2019;16(2):183–90.
Webb B, Sali A. Comparative protein structure modeling using MODELLER. Curr Protoc Protein Sci. 2016;86:291–2937.
Zhang Y, Skolnick J. TM-align: a protein structure alignment algorithm based on the TM-score. Nucleic Acids Res. 2005;33(7):2302–9.
Cao Y, Song L, Miao Z, et al. Improved side-chain modeling by coupling clash-detection guided iterative search with rotamer relaxation. Bioinformatics. 2011;27(6):785–90.
MacKerell AD, Bashford D, Bellott M, et al. All-atom empirical potential for molecular modeling and dynamics studies of proteins. J Phys Chem B. 1998;102(18):3586–616.
Liu Y, Grimm M, Dai WT, et al. CB-Dock: a web server for cavity detection-guided protein-ligand blind docking. Acta Pharmacol Sin. 2020;41(1):138–44.
Chung WH, Hung SI, Yang JY, et al. Granulysin is a key mediator for disseminated keratinocyte death in Stevens-Johnson syndrome and toxic epidermal necrolysis. Nat Med. 2008;14(12):1343–50.
Su SC, Mockenhaupt M, Wolkenstein P, et al. Interleukin-15 is associated with severity and mortality in stevens-johnson syndrome/toxic epidermal necrolysis. J Invest Dermatol. 2017;137(5):1065–73.
Tsai YG, Liou JH, Hung SI, et al. Increased type 2 innate lymphoid cells in patients with drug reaction with eosinophilia and systemic symptoms syndrome. J Invest Dermatol. 2019;139(8):1722–31.
Stern RS, Divito SJ. Stevens-johnson syndrome and toxic epidermal necrolysis: associations, outcomes, and pathobiology-thirty years of progress but still much to be done. J Invest Dermatol. 2017;137(5):1004–8.
Shiohara T, Mizukawa Y. Drug-induced hypersensitivity syndrome (DiHS)/drug reaction with eosinophilia and systemic symptoms (DRESS): an update in 2019. Allergol Int. 2019;68(3):301–8.
Watanabe H, Watanabe Y, Tashiro Y, et al. A docking model of dapsone bound to HLA-B*13:01 explains the risk of dapsone hypersensitivity syndrome. J Dermatol Sci. 2017;88(3):320–9.
Pichler WJ. The p-i concept: pharmacological interaction of drugs with immune receptors. World Allergy Organ J. 2008;1(6):96–102.
Yun J, Mattsson J, Schnyder K, et al. Allopurinol hypersensitivity is primarily mediated by dose-dependent oxypurinol-specific T cell response. Clin Exp Allergy. 2013;43(11):1246–55.
Lucas A, Nolan D, Mallal S. HLA-B*5701 screening for susceptibility to abacavir hypersensitivity. J Antimicrob Chemother. 2007;59(4):591–3.
Illing PT, Mifsud NA, Purcell AW. Allotype specific interactions of drugs and HLA molecules in hypersensitivity reactions. Curr Opin Immunol. 2016;42:31–40.
Ko TM, Chung WH, Wei CY, et al. Shared and restricted T-cell receptor use is crucial for carbamazepine-induced Stevens-Johnson syndrome. J Allergy Clin Immunol. 2011;128(6):1266–76.
Chu MT, Wang CW, Chang WC, et al. Granulysin-based lymphocyte activation test for evaluating drug causality in antiepileptics-induced severe cutaneous adverse reactions. J Invest Dermatol. 2021;141(6):1461–72.
We also thank drug-specific T cell culture support of Xiaofeng Qin and Heng Yang from Suzhou Institute of Systems Medicine, Chinese Academy of Medical Sciences & Peking Union Medical College, and the genetic TCR analysis support of members of the Whole-Genome Research Core Laboratory of Human Diseases, and the Cancer Vaccine and Immune Cell Therapy Core Laboratory, Chang Gung Memorial Hospital, Taiwan. The authors would like to thank the staff at the beamline BL17U1 of Shanghai Synchrotron Radiation Facility (SSRF) for crystal data collection.
This work was supported by grants from the grants from the National Key R&D Program of China (2019YFE0113300), National Natural Science Foundation of China (Grant 81972950, 82173431, 82103748, 31770948 and 82073595), CAMS Initiative for Innovative Medicine (Grant 2016-I2M-1-005, 2017-I2M-B&R-14), the Ministry of Science and Technology, Taiwan (MOST 109-2326-B-182A-001-, 108-2314-B-182A-104 -MY3, 108-2320-B-182A-024 -MY2, 107-2811-B-182A-509-), Chang Gung Memorial Hospital (CMRRPG3K0562, CMRPG3K2181, CMRPG3I0382, CIRPG3I0022), the Fujian Provincial Department of Science and Technology (2020Y4007 and 2021H0004, SO), and the High-Level Personnel Introduction Grant of Fujian Normal University (Z0210509, SO).
Ethics approval and consent to participate
Approvals were obtained from the institutional review board (IRB No. 100-4657A3, 201601761B0, 201902171A3 and 202001645B0), and informed consent was obtained from each participant.
Consent for publication
All authors agree for publication.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Flow cytometry analysis of DDS-specific T cells in patients with DIHS. A We performed in vitro T cells expansion for day0, day21 and day35 of PBMC from DDS-induced DIHS patients and tolerant donors, as measured the ratios of CD8+ and CD4+ T cell by flow cytometry. B, C The ratios of CD8+ and CD4+ T cells of DDS-induced DIHS patients and tolerant individuals after T cells expansion were shown. The entire experiment was repeated thrice. *P < 0.05, **P < 0.01, ***P < 0.001, by two-tailed Student t test. Figure S2. The results of IFN-γ ELISpot assays. A DDS-specific CTLs activation was detected by IFN-γ ELISpot assays in a DDS-concentration manner. B T cells-mediated were blocked by anti–HLA class I antibody but not by anti-HLA class II antibody. Figure S3. Determination of the cross-reactivity of DDS and other sulfa drugs for dapsone-SCAR patients carried HLA-B*13:01. A Granulysin-based lymphocyte activation test (LAT) was performed in 6 dapsone-induced SCAR patients carried HLA-B*13:01 and 5 tolerant controls. A positive result was defined as a 1.3-fold increase in granulysin release compared to the controls (dashed blue line). Black solid dot, dapsone-SJS patient; blue solid dot, dapsone-DRESS patient, gray solid dot, dapsone-tolerant subject. B the sensitivity and specificity of LAT for dapsone-induced SCAR patients carried HLA-B*13:01 and their cross-reactive to dapsone are showed. Figure S4. SDS-PAGE analysis for β2M and HLA-B*13:01 protein. A 1. Purified HLA-B*13:01; 2. purified β2M; 3. Marker; 4. HLA-B*13:01-β2M complex. B the recombinant HLA-B*13:01 protein bounded to DDS in a DDS-concentration manner, with an estimated low-affinity of the micromolar range (Kd = 68.54 M). Figure S5. Pulsing assays and LC–MS analysis peptides involved in DDS binding with HLA-B*13:01. A Peptides were analyzed on C1R (HLA-B*13:01) with DDS cocultured by Co-IP and LC–MS. B T Lymphocyte activation test was performed from C1R (HLA-B*13:01) cultured with DDS or two peptides by IFN-γ release. Figure S6. Flow cytometry analysis of DDS-specific TCR transfectants. TRA12-3/TRB28, TCRA13-1/TRB5-1 and TCRA12-3/TRB30 were reconstituted and overexpressed into PBMC form DDS tolerant subjects, and the expression levels of TCRαβ and GFP were determined by FACS.
TCR primers used in RT-PCR analysis.
About this article
Cite this article
Jiang, H., Wang, CW., Wang, Z. et al. Functional and structural characteristics of HLA-B*13:01-mediated specific T cells reaction in dapsone-induced drug hypersensitivity. J Biomed Sci 29, 58 (2022). https://doi.org/10.1186/s12929-022-00845-8
- Dapsone-induced hypersensitivity syndrome
- T-cell receptor