Next Article in Journal
TRAIL/NF-κB/CX3CL1 Mediated Onco-Immuno Crosstalk Leading to TRAIL Resistance of Pancreatic Cancer Cell Lines
Next Article in Special Issue
Inter-Individual Variability in Acute Toxicity of R-Pulegone and R-Menthofuran in Human Liver Slices and Their Influence on miRNA Expression Changes in Comparison to Acetaminophen
Previous Article in Journal
Cx43 Channel Gating and Permeation: Multiple Phosphorylation-Dependent Roles of the Carboxyl Terminus
Previous Article in Special Issue
Role of 3-Hydroxy Fatty Acid-Induced Hepatic Lipotoxicity in Acute Fatty Liver of Pregnancy
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 19
Issue 6
10.3390/ijms19061660
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Interaction of Nevirapine with the Peptide Binding Groove of HLA-DRB1*01:01 and Its Effect on the Conformation of HLA-Peptide Complex

by
Makoto Hirasawa
1,*,
Katsunobu Hagihara
2,
Koji Abe
1,
Osamu Ando
1 and
Noriaki Hirayama
3
1
Drug Metabolism & Pharmacokinetics Research Laboratories, Daiichi Sankyo Co., Ltd., 1-2-58 Hiromachi, Shinagawa-ku, Tokyo 140-8710, Japan
2
Biomarker Department, Daiichi Sankyo Co., Ltd., 1-2-58 Hiromachi, Shinagawa-ku, Tokyo 140-8710, Japan
3
Institute of Advanced Biosciences, Tokai University, 4-1-1 Kitakaname, Hiratsuka-shi, Kanagawa 259-1292, Japan
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2018, 19(6), 1660; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms19061660
Submission received: 14 May 2018 / Revised: 28 May 2018 / Accepted: 29 May 2018 / Published: 4 June 2018
(This article belongs to the Special Issue Hepatotoxicity: Molecular Mechanisms and Pathophysiology)
Browse Figures
Versions Notes

Abstract

:
Human leukocyte antigen (HLA)-DRB1*01:01 has been shown to be involved in nevirapine-induced hepatic hypersensitivity reactions. In the present study, in silico docking simulations and molecular dynamics simulations were performed to predict the interaction mode of nevirapine with the peptide binding groove of HLA-DRB1*01:01 and its possible effect on the position and orientation of the ligand peptide derived from hemagglutinin (HA). In silico analyses suggested that nevirapine interacts with HLA-DRB1*01:01 around the P4 pocket within the peptide binding groove and the HA peptide stably binds on top of nevirapine at the groove. The analyses also showed that binding of nevirapine at the groove will significantly change the inter-helical distances of the groove. An in vitro competitive assay showed that nevirapine (1000 μM) increases the binding of the HA peptide to HLA-DRB1*01:01 in an allele-specific manner. These results indicate that nevirapine might interact directly with the P4 pocket and modifies its structure, which could change the orientation of loaded peptides and the conformation of HLA-DRB1*01:01; these changes could be distinctively recognized by T-cell receptors. Through this molecular mechanism, nevirapine might stimulate the immune system, resulting in hepatic hypersensitivity reactions.

Graphical Abstract

1. Introduction

Nevirapine (Viramune®) is an orally available antiretroviral agent, which inhibits reverse transcriptase. It is used in combination with other antiretroviral agents for the treatment of HIV/AIDS infections and is included in the WHO Model List of Essential Medicines. Although it is generally well-tolerated and effective, it has been associated with hypersensitivity reactions (HSRs) in approximately 5% of patients with various clinical phenotypes, including cutaneous, hepatic, and systemic symptoms [1,2,3,4]. Several features of nevirapine HSRs, including the delayed onset of reaction, more rapid and severe development of symptoms after rechallenge, and association with the percentage of pretreatment CD4 T-cells, suggest that nevirapine HSRs are immune-mediated [3,5].
Various genetic association studies have shown that cytochrome P450 (CYP) 2B6 genotypes and human leukocyte antigen (HLA) alleles are predisposed to the development of nevirapine HSRs [4,5,6,7,8]; however, these associations were very complicated. For example, CYP2B6 G516T polymorphism increased nevirapine plasma exposure [9] and was associated with its cutaneous HSRs [6], but not with hepatic ones [6,10]. In addition, unlike the cases of HLA-B*57:01 in abacavir HSRs [11] and HLA-B*58:01 in allopurinol-induced Stevens-Johnson syndrome (SJS) and toxic epidermal necrolysis (TEN) [12], multiple class I and class II HLA alleles, which have variable distribution and risk across ethnic groups, have been associated with nevirapine HSRs [6,13,14]. Interestingly, associations with protective HLA alleles have also been reported [7,13,15]. Furthermore, the associations with predisposing HLA alleles tend to be phenotype-specific [6,14,16], where HLA-C*04 has been associated with cutaneous HSRs, in particular SJS/TEN [6,7,8,13,14], but HLA-DRB1*01 has been associated with hepatic HSRs. The exclusion of individuals with coincident cutaneous and hepatic HSRs resulted in a significant increase of the odds ratios [5,6]. These results suggest that cutaneous and hepatic HSRs might have fundamentally different mechanisms.
We previously reported in silico and in vitro interactions of lapatinib [17] and ximelagatran [18] with HLA-DRB1*07:01, which predisposed to the development of idiosyncratic hepatotoxicity by both drugs, and indicated that these interactions might trigger immune responses, leading to drug-related toxicity. In the present study, we performed in silico docking simulations to predict the interaction mode of nevirapine (Figure 1) with the peptide binding groove of HLA-DRB1*01:01, as well as molecular dynamics (MD) simulations to predict the possible effect of nevirapine on the structure of HLA-DRB1*01:01 and the position and orientation of the bound peptide derived from hemagglutinin (HA). An in vitro competitive assay was performed to evaluate the effect of nevirapine on the binding of ligand peptides to three HLA-DR alleles.

2. Results

2.1. Docking Simulations

In silico docking simulations were performed to predict the interaction mode of nevirapine with the peptide binding groove of HLA-DRB1*01:01. The binding affinity of nevirapine to HLA-DRB1*01:01 was judged by a scoring function of generalized-Born volume integral/weighted surface area (GBVI/WSA_dG). The lowest GBVI/WSA_dG value of the complex between nevirapine and HLA-DRB1*01:01 was −5.18 kcal/mol. This value is higher than that of the complex between HLA-B*14:02 and nevirapine (−6.40 kcal/mol) simulated by the same software system [19], indicating that the interaction of nevirapine with HLA-DRB1*01:01 was weaker than that with HLA-B*14:02. The binding mode of nevirapine at the peptide binding groove of HLA-DRB1*01:01 with the lowest GBVI/WSA_dG value is shown in Figure 2. Nevirapine sprawls out across the pocket P4–P6.

2.2. Molecular Dynamics (MD) Simulations

MD simulations were performed to evaluate the possible effect of nevirapine on the structure of the peptide binding groove of HLA-DRB1*01:01 and the position and orientation of bound HA peptide. Each of the total energies of MD simulations stabilized within approximately 0.5 ns (Figure 3a) indicating that the simulations were energetically stable. The flexibility of HLA-DRB1*01:01 molecule was very similar in all simulations (Figure 3b). As shown in Figure 3c, the root mean square deviation (RMSD) values stabilized after 1–2 ns in all simulations, which indicated that the systems relaxed quickly.
The interaction of nevirapine with HLA-DRB1*01:01 resulted in the closing motion of the peptide binding groove (Figure 3d, Table 1), which differs from the tightly closed binding groove observed with lapatinib-DRB1*07:01 [17] and the open groove observed with lapatinib-DRB1*01:01, lapatinib-DRB1*15:01, [17] and ximelagatran-DRB1*07:01 [18]. The mean inter-helical distance of the peptide binding groove of HLA-DRB1*01:01 in the presence of nevirapine (13.1 ± 0.4 Å; Table 1) was very similar to that observed in apo simulation (13.0 ± 0.4 Å), which reflected the transition of HLA-DR molecules from a peptide-receptive conformation to a non-receptive one in the absence of any ligand molecules [20]. This is reasonable considering the much smaller molecular size of nevirapine than that of lapatinib and ximelagatran.
Figure 4 shows structural analyses of five sets of MD simulations presented in three panels for each. In the absence of HA peptide, nevirapine bound to HLA-DRB1*01:01 in or near the P4 pocket in a stable conformation (Figure 4c) and had a small impact on the overall configuration of the α and β chains of HLA-DRB1*01:01, compared to that of lapatinib [17]. Regarding HA peptide-bound simulations, three starting frames (2, 3, and 4) were considered, and frame 3 was the most likely because it was stable and exhibited good contacts with both sides of the groove. The conformation was the closest to the X-ray crystal structure (Figure 4d–f). Frame 2 was also considered for the trimer (nevirapine, HLA-DRB1*01:01, and HA peptide) simulations because of its total energy, although the binding mode and closing motion of the peptide binding groove were somehow unusual (Figure 4g–i, Table 1). In the trimer simulations, conformational changes of both nevirapine and HA peptide during 5-ns simulations were slight in both frames. Nevirapine seems to interact mainly with the peptide binding groove around the P4 pocket, whereas HA peptide bound on top of nevirapine. Slight conformational changes were observed in the α and β chains of HLA-DRB1*01:01 with the α chain helix being maintained. Considering the high root mean square fluctuation (RMSF) values of HA peptide ends in frame 2, the HA peptide in frame 3 appears more stably bound to HLA-DRB1*01:01 in the trimer simulations.
Figure 5 shows the differences in the positions and orientations of ligand peptides between representative structures of trimer and HLA-peptide dimer simulations. Nevirapine had much smaller effects on the position and orientation of the ligand peptide than lapatinib and ximelagatran.

2.3. HLA Class II Competitive Assay

To evaluate the effect of nevirapine on the binding of ligand peptides to three HLA-DR molecules, a DELFIA-based competitive assay was conducted. Nevirapine increased the binding of HA peptide to soluble HLA-DRB1*01:01 molecule only at 1000 μM (p < 0.0001), whereas there were no significant or concentration-dependent effects on the binding of the peptide derived from tetanus toxoid (TT) to HLA-DRB1*07:01 and that derived from myelin basic protein (MBP) to HLA-DRB1*15:01 (Table 2).
The concentration of nevirapine that possibly bound to HLA-DR molecules in the competitive assay sample, where nominally 1000 μM of nevirapine was applied, was measured by liquid chromatography-tandem mass spectrometry (LC-MS/MS). However, a much higher concentration than the theoretically highest one (4.7 nM) with a large variance was detected in the absence of HLA-DR molecules, which indicated high non-specific binding of nevirapine to the plastic tubes and plates in the absence of HLA-DR molecules. In addition, the nevirapine concentration in the presence of both the ligand peptide and HLA-DR was lower than that in the absence of the ligand peptide for all three HLA-DR alleles, which indicated that the presence of the ligand peptide also attenuated the non-specific binding of nevirapine to tubes and plates. Taken together, the nevirapine bound to HLA-DR molecules could not be detected in this experiment owing to its high non-specific binding.

3. Discussion

Associations of multiple phenotypes of nevirapine HSRs with various HLA class I and class II alleles across several ethnic groups have been reported [4,5,6,7,8,13,14,16]. Among the multiple HLA alleles associated with nevirapine HSRs, HLA-C*04:01 is most often associated and commonly carried across ethnic groups [13]. Several docking studies recently have shown the binding of nevirapine in the B and/or F pockets [8,13] within the peptide binding groove of HLA-C*04:01. Pavlos et al. [13] reported that HLA-C risk alleles for nevirapine cutaneous HSRs (*04:01, *05:01, and *18:01) shared the unique F pocket motif, as well as Arg156. They proposed that the disease-causing peptides were anchored in the F pocket together with nevirapine and stabilized by Arg156 in the central portion (P3-P5-P6), which could propagate T-cell mediated responses. In addition, they reported that the P4 pocket motif was also shared by HLA-DRB1 risk alleles (*01:(01/02/03) and *04:(04/05/08/10)).
Present in silico studies have shown that nevirapine sprawls out across the P4–P6 pockets in the docking simulations (Figure 2) and binds to the P4 pocket in a stable conformation both in the presence and absence of HA peptide in the MD simulations (Figure 4). These findings are in accordance with the reported predisposing effect of the shared P4 pocket motif. It is noteworthy that the present docking simulations indicated that the interaction of nevirapine with the peptide binding groove of HLA-DRB1*01:01 was weaker than that of other idiosyncratic drug toxicity-causing drugs, such as ximelagatran [18] and allopurinol [21].
Moreover, Pavlos et al. reported that nevirapine did not affect the repertoire of peptides presented on HLA-DRB1*01:01 in L2 cells, and any peptide eluted from nevirapine-treated cells did not show a significantly higher binding affinity to HLA-DRB1*01:01 in the presence of nevirapine than in the absence of it. These results were unexpected considering the relatively small molecular size and the simulated direct interaction mode of nevirapine with the P4 pocket of HLA-DRB1*01:01, which indicates the likelihood of the “altered-repertoire” mechanism of nevirapine-induced immune stimulation, similar to abacavir hypersensitivity [22].
However, in the present in vitro competitive assay, nevirapine increased the binding of HA peptide to HLA-DRB1*01:01 in an allele-specific manner at 1000 μM (Table 2), although its concentration dependency was not clear, in comparison to that of lapatinib, which increased the binding of the TT peptide to HLA-DRB1*07:01 in a concentration-dependent manner [17]. The lack of effect at lower concentrations of nevirapine might be consistent with its relatively low affinity to HLA-DRB1*01:01, as shown in the docking simulations. This peculiar concentration dependency might account for the negative results of the elution studies conducted by Pavlos et al. [13], where nominally 100 μg/mL (375.5 μM) of nevirapine was applied to the cells.
Additionally, nevirapine might not influence the binding affinity of the ligand peptides to HLA-DRB1*01:01 but could appreciably affect their conformation and orientation at lower concentrations, which are undetectable in the present competitive assay and probably in the reported elution studies. This hypothesis is supported by the present MD simulation results, which suggested that HA peptide can bind to HLA-DRB1*01:01 on top of nevirapine bound to the P4 pocket (Figure 4) but the effect of nevirapine on the position and orientation of the HA peptide on HLA-DRB1*01:01 was much less than those of ximelagatran and lapatinib (Figure 5). This indicates that nevirapine might neither inhibit nor enhance the binding of the ligand peptides to HLA-DRB1*01:01 robustly, which is unlike the results observed with HLA-DRB1*07:01-ximelagatran [18] and HLA-DRB1*07:01-lapatinib [17].
The binding mode of HA peptide in the trimer simulations was seemingly contradictory to the results of molecular modeling reported by Pavlos et al. [13], where the leucine residue of the sequence logo peptide (YRSLKAQRV), which was identified from peptides bound to HLA-DRB1*01:01 in nevirapine-treated cells, occupied the P4 pocket, and thus nevirapine was considered unlikely to bind directly to HLA-DRB1*01:01 in the presence of the peptide. However, in the present MD simulations (frame 3 trimer), P1–P5 (Tyr-Val-Lys-Gln-Asn) moved slightly to cover nevirapine with Tyr and Lys, having a non-bonded interaction with nevirapine. Additionally, the α and β helices slightly changed their conformations to accommodate the binding of the HA peptide. This finding showed the advantage of MD simulations in predicting the possible effects of small-molecule drugs on the position and orientation of the ligand peptides in the binding groove of HLA class II molecules.
As shown in Figure 5, the degree of nevirapine-induced changes in the position and orientation of the HA peptide in HLA-DRB1*01:01 was much smaller than that of lapatinib- and ximelagatran-induced changes in the position of the TT peptide in HLA-DRB1*07:01. The relevant small conformational change in the HLA-peptide complex was considered to be insufficient to trigger the immune response. However, the mechanism of immune stimulation induced by small structural changes in the surface of HLA-ligand complexes was clearly shown for chronic beryllium disease (CBD) associated with HLA-DP2 [23]. Petukh et al. [24] analyzed the conformational changes in HLA-DP2 molecule induced by the binding of the ligand peptides in the presence and absence of Be2+ cation by comparing the inter-helical distances in HLA-DP2 using MD simulations. The simulations indicated that the conformational changes induced by the binding of the ligand peptides in the presence of Be2+ were distinctly different from those caused by natural peptides, which could trigger immune responses in the absence of Be2+, and these conformational changes were expected to trigger a completely different signaling cascade of immune reactions, resulting in CBD.
In the present MD simulations, nevirapine induced significant differences in the mean inter-helical distances of the peptide binding groove of HLA-DRB1*01:01 (Table 1). Accordingly, the interaction of nevirapine with HLA-DRB1*01:01 might induce conformational changes not only in the ligand peptide but also in HLA-DRB1*01:01. These changes might be recognized differently by T cell receptors (TCRs) and trigger nevirapine-induced hepatic HSRs.

4. Materials and Methods

4.1. Docking Simulations

A crystal structure of HLA-DRB1*01:01 (PDB ID: 3PDO) deposited at the Protein Data Bank [25] was used in this study. The binding mode and affinity of nevirapine to the peptide binding groove of HLA-DRB1*01:01 were determined by docking simulations using the program ASEDock [26]. The scoring function of GBVI/WSA_dG [27] was used to evaluate the binding affinity. The software system molecular operating environment (MOE) [28] was used.

4.2. MD Simulations

A crystal structure of HLA-DRB1*01:01 (PDB ID: 1AQD) was used in this study. The top scoring mode of nevirapine with HLA-DRB1*01:01 in the docking simulations performed using AutoDock Vina [29] was selected as a starting conformation for MD simulation. HA peptide (HA 306–318; PRYVKQNTLKLAT) was used as a ligand peptide for HLA-DRB1*01:01. The structure of the HA peptide was generated by homology modeling using a template structure (PDB ID: 1DLH, chain C; PKYVKQNTLKLAT). To select the correct peptide registration, we conducted a short translational scan by sliding the HA peptide along the axis of the binding groove and running additional MD simulations. After modeling a series of peptide registration frames, frame 2 and frame 3 were chosen as starting points for MD simulations, considering the energetic stability, contact with α and β chains, and EpiMatrix-predicted registration. MD simulations were carried out as previously described [17]. Briefly, each system was protonated at pH 7.4 and then solvated in a periodic rectangular box of TIP3P water with 10-Å padding between the edge of the box and the nearest solute atom. Sodium ions were added to neutralize the system. The system was relaxed and then gradually heated to 300 K. Production simulations were performed in the canonical constant-temperature, constant-volume (NVT) ensemble. An integration time step of 4 fs was used. The “representative structure” of each simulation was selected as the conformation, which was the most representative of the conformations assumed over the final 1000 frames of the simulations. MD simulations of nevirapine-bound HLA-DRB1*01:01 in the absence and presence of HA peptide lasted 5 ns for each.

4.3. HLA Class II Competitive Assay

The competitive assay was conducted as previously reported with minor modifications [30]. Briefly, soluble HLA-DR molecules (20 nM), provided by Benaroya Research Institute (Seattle, WA, USA), were incubated in quadruplicate with nevirapine at concentrations ranging between 8 and 1000 μM and biotinylated control ligand peptides (HA 306–318: PRYVKQNTLKLAT for HLA-DRB1*01:01, TT 830–844: QYIKANSKFIGITEL for HLA-DRB1*07:01, and MBP 94–112: NPVVHFFKNIVTPRTPPPS for HLA-DRB1*15:01) at concentrations lower than their Michaelis-Menten constants for 48 h at 37 °C to reach equilibrium. HLA-peptide complexes were then captured onto 96-well plates coated with L243 anti-HLA-DR antibody (BioLegend, San Diego, CA, USA). The plates were washed and incubated with Europium-labeled streptavidin (Perkin-Elmer, Waltham, MA, USA) for 1 h at room temperature. The plates were washed again, and DELFIA enhancement solution (Perkin-Elmer) was added for 20 min at room temperature to develop the plates before they were read using an EnVision plate reader. An aliquot of the incubation sample containing 1000 μM of nevirapine was also captured onto other L243-coated plates. The plates were washed, and nevirapine bound to HLA-DR molecules was extracted with 150 μL of acetonitrile. The concentrations of nevirapine in the extracted samples were determined by LC-MS/MS using a Shim-pack XR-ODS (30 × 2 mm i.d., 2.2 μm) column and Prominence LC-20A system (Shimadzu Corp., Kyoto, Japan) coupled with a Sciex API 4000 triple quadrupole mass spectrometer (Applied Biosystems, Foster City, CA, USA). The aqueous mobile phase (A) was H2O/100 mM CH3CO2NH4/CH3CN (900/50/50, v/v/v), and the organic mobile phase (B) was CH3CN/100 mM CH3CO2NH4 (1000/50, v/v). Gradient elution was used as follows: 50% B for the first 0.5 min, increased to 100% B from 0.5 to 1 min, and maintained at 100% B for additional 0.75 min. The total run time was 1.75 min at a flow rate of 0.75 mL/min. Ionization was conducted in the positive ion mode using the transition m/z 267.1 to m/z 226.1. The retention time of nevirapine was 0.96 min. The lower and upper limits of quantification of the assay were 0.03 and 300 nM, respectively.

4.4. Statistical Analysis

A two-tailed Dunnett’s test was carried out using SAS System Release 9.2 (SAS Institute Inc., Cary, NC, USA) to compare the relative fluorescence counts (% of dimethyl sulfoxide control for each HLA-DR allele) in the presence and absence of nevirapine. A p-value < 0.05 was considered statistically significant.

5. Conclusions

The present study suggested that nevirapine might directly interact with the peptide binding groove of HLA-DRB1*01:01 around the P4 pocket, which could induce conformational changes in the ligand peptide and HLA-DRB1*01:01. These conformational changes recognized by TCR differently from the natural peptide-HLA-DRB1*01:01 complex would stimulate the immune system and trigger nevirapine-induced hepatic HSRs.

Author Contributions

M.H., K.H. and N.H. conceived and designed the experiments; M.H. and N.H. performed the experiments; M.H. and N.H. analyzed the data; M.H., K.A., O.A. and N.H. contributed to reagents/materials/analysis tools; and M.H. wrote this manuscript. All authors read and approved the final manuscript.

Acknowledgments

We would like to thank Anne S. De Groot and Frances Terry (EpiVax, Inc.) for their helpful comments and discussion on this study.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

HLAhuman leukocyte antigen
HAhemagglutinin
HSRhypersensitivity reaction
CYPcytochrome P450
SJSStevens-Johnson syndrome
TENtoxic epidermal necrolysis
MDmolecular dynamics
GBVI/WSA_dGgeneralized-Born volume integral/weighted surface area
RMSDroot mean square deviation
RMSFroot mean square fluctuation
SDstandard deviation
TTtetanus toxoid
MBPmyelin basic protein
LC-MS/MSliquid chromatography-tandem mass spectrometry
DMSOdimethyl sulfoxide
CBDchronic beryllium disease
TCRT-cell receptor

References

  1. Dieterich, D.T.; Robinson, P.A.; Love, J.; Stern, J.O. Drug-induced liver injury associated with the use of nonnucleoside reverse-transcriptase inhibitors. Clin. Infect. Dis. 2004, 38 (Suppl. 2), S80–S89. [Google Scholar] [CrossRef] [PubMed]
  2. Wit, F.W.; Kesselring, A.M.; Gras, L.; Richter, C.; van der Ende, M.E.; Brinkman, K.; Lange, J.M.; de Wolf, F.; Reiss, P. Discontinuation of nevirapine because of hypersensitivity reactions in patients with prior treatment experience, compared with treatment-naive patients: The ATHENA cohort study. Clin. Infect. Dis. 2008, 46, 933–940. [Google Scholar] [CrossRef] [PubMed]
  3. Stern, J.O.; Robinson, P.A.; Love, J.; Lanes, S.; Imperiale, M.S.; Mayers, D.L. A comprehensive hepatic safety analysis of nevirapine in different populations of HIV infected patients. J. Acquir. Immune Defic. Syndr. 2003, 34 (Suppl. 1), S21–S33. [Google Scholar] [CrossRef] [PubMed]
  4. Keane, N.M.; Pavlos, R.K.; McKinnon, E.; Lucas, A.; Rive, C.; Blyth, C.C.; Dunn, D.; Lucas, M.; Mallal, S.; Phillips, E. HLA Class I restricted CD8+ and Class II restricted CD4+ T cells are implicated in the pathogenesis of nevirapine hypersensitivity. AIDS 2014, 28, 1891–1901. [Google Scholar] [CrossRef] [PubMed]
  5. Martin, A.M.; Nolan, D.; James, I.; Cameron, P.; Keller, J.; Moore, C.; Phillips, E.; Christiansen, F.T.; Mallal, S. Predisposition to nevirapine hypersensitivity associated with HLA-DRB1*0101 and abrogated by low CD4 T-cell counts. AIDS 2005, 19, 97–99. [Google Scholar] [CrossRef] [PubMed]
  6. Yuan, J.; Guo, S.; Hall, D.; Cammett, A.M.; Jayadev, S.; Distel, M.; Storfer, S.; Huang, Z.; Mootsikapun, P.; Ruxrungtham, K.; et al. Toxicogenomics of nevirapine-associated cutaneous and hepatic adverse events among populations of African, Asian, and European descent. AIDS 2011, 25, 1271–1280. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Carr, D.F.; Chaponda, M.; Jorgensen, A.L.; Castro, E.C.; van Oosterhout, J.J.; Khoo, S.H.; Lalloo, D.G.; Heyderman, R.S.; Alfirevic, A.; Pirmohamed, M. Association of human leukocyte antigen alleles and nevirapine hypersensitivity in a Malawian HIV-infected population. Clin. Infect. Dis. 2013, 56, 1330–1339. [Google Scholar] [CrossRef] [PubMed]
  8. Carr, D.F.; Bourgeois, S.; Chaponda, M.; Takeshita, L.Y.; Morris, A.P.; Castro, E.M.; Alfirevic, A.; Jones, A.R.; Rigden, D.J.; Haldenby, S.; et al. Genome-wide association study of nevirapine hypersensitivity in a sub-Saharan African HIV-infected population. J. Antimicrob. Chemother. 2017, 72, 1152–1162. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Penzak, S.R.; Kabuye, G.; Mugyenyi, P.; Mbamanya, F.; Natarajan, V.; Alfaro, R.M.; Kityo, C.; Formentini, E.; Masur, H. Cytochrome P450 2B6 (CYP2B6) G516T influences nevirapine plasma concentrations in HIV-infected patients in Uganda. HIV Med. 2007, 8, 86–91. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Haas, D.W.; Bartlett, J.A.; Andersen, J.W.; Sanne, I.; Wilkinson, G.R.; Hinkle, J.; Rousseau, F.; Ingram, C.D.; Shaw, A.; Lederman, M.M.; et al. Pharmacogenetics of nevirapine-associated hepatotoxicity: An Adult AIDS Clinical Trials Group collaboration. Clin. Infect. Dis. 2006, 43, 783–786. [Google Scholar] [CrossRef] [PubMed]
  11. Chessman, D.; Kostenko, L.; Lethborg, T.; Purcell, A.W.; Williamson, N.A.; Chen, Z.; Kjer-Nielsen, L.; Mifsud, N.A.; Tait, B.D.; Holdsworth, R.; 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, 822–832. [Google Scholar] [CrossRef] [PubMed]
  12. Somkrua, R.; Eickman, E.E.; Saokaew, S.; Lohitnavy, M.; Chaiyakunapruk, N. Association of HLA-B*5801 allele and allopurinol-induced Stevens Johnson syndrome and toxic epidermal necrolysis: A systematic review and meta-analysis. BMC Med. Genet. 2011, 12. [Google Scholar] [CrossRef] [PubMed]
  13. Pavlos, R.; McKinnon, E.J.; Ostrov, D.A.; Peters, B.; Buus, S.; Koelle, D.; Chopra, A.; Schutte, R.; Rive, C.; Redwood, A.; et al. Shared peptide binding of HLA Class I and II alleles associate with cutaneous nevirapine hypersensitivity and identify novel risk alleles. Sci. Rep. 2017, 7. [Google Scholar] [CrossRef] [PubMed]
  14. Cornejo Castro, E.M.; Carr, D.F.; Jorgensen, A.L.; Alfirevic, A.; Pirmohamed, M. HLA-allelotype associations with nevirapine-induced hypersensitivity reactions and hepatotoxicity: A systematic review of the literature and meta-analysis. Pharmacogenet. Genom. 2015, 25, 186–198. [Google Scholar] [CrossRef] [PubMed]
  15. Umapathy, S.; Pawar, A.; Bajpai, S.; Pazare, A.R.; Ghosh, K. HLA involvement in nevirapine-induced dermatological reaction in antiretroviral-treated HIV-1 patients. J. Pharmacol. Pharmacother. 2011, 2, 114–115. [Google Scholar] [CrossRef] [PubMed]
  16. Pavlos, R.; Mallal, S.; Phillips, E. HLA and pharmacogenetics of drug hypersensitivity. Pharmacogenomics 2012, 13, 1285–1306. [Google Scholar] [CrossRef] [PubMed]
  17. Hirasawa, M.; Hagihara, K.; Okudaira, N.; Izumi, T. The possible mechanism of idiosyncratic lapatinib-induced liver injury in patients carrying human leukocyte antigen-DRB1*07:01. PLoS ONE 2015, 10. [Google Scholar] [CrossRef] [PubMed]
  18. Hirasawa, M.; Hagihara, K.; Abe, K.; Ando, O.; Hirayama, N. In Silico and In Vitro analysis of interaction between ximelagatran and human leukocyte antigen (HLA)-DRB1*07:01. Int. J. Mol. Sci. 2017, 18. [Google Scholar] [CrossRef] [PubMed]
  19. Isogai, H.; Hirayama, N. In Silico analysis of interactions between nevirapine-related compounds, HLA-B*14:02 and T-cell receptor. Chem-Bio Inf. J. 2016, 16, 9–12. [Google Scholar]
  20. Rabinowitz, J.D.; Vrljic, M.; Kasson, P.M.; Liang, M.N.; Busch, R.; Boniface, J.J.; Davis, M.M.; McConnell, H.M. Formation of a Highly peptide-receptive state of class II MHC. Immunity 1998, 9, 699–709. [Google Scholar] [CrossRef]
  21. Osabe, M.; Tohkin, M.; Hirayama, N. In Silico analysis of interactions between HLA-B*58:01 and allopurinol-related compounds. Chem-Bio Inf. J. 2016, 16, 1–4. [Google Scholar] [CrossRef]
  22. Illing, P.T.; Vivian, J.P.; Dudek, N.L.; Kostenko, L.; Chen, Z.; Bharadwaj, M.; Miles, J.J.; Kjer-Nielsen, L.; Gras, S.; Williamson, N.A.; et al. Immune self-reactivity triggered by drug-modified HLA-peptide repertoire. Nature 2012, 486, 554–558. [Google Scholar] [CrossRef] [PubMed]
  23. Clayton, G.M.; Wang, Y.; Crawford, F.; Novikov, A.; Wimberly, B.T.; Kieft, J.S.; Falta, M.T.; Bowerman, N.A.; Marrack, P.; Fontenot, A.P.; et al. Structural basis of chronic beryllium disease: Linking allergic hypersensitivity and autoimmunity. Cell 2014, 158, 132–142. [Google Scholar] [CrossRef] [PubMed]
  24. Petukh, M.; Wu, B.; Stefl, S.; Smith, N.; Hyde-Volpe, D.; Wang, L.; Alexov, E. Chronic beryllium disease: Revealing the role of beryllium ion and small peptides binding to HLA-DP2. PLoS ONE 2014, 9. [Google Scholar] [CrossRef] [PubMed]
  25. Bernstein, F.C.; Koetzle, T.F.; Williams, G.J.; Meyer, E.F., Jr.; Brice, M.D.; Rodgers, J.R.; Kennard, O.; Shimanouchi, T.; Tasumi, M. The protein data bank: A computer-based archival file for macromolecular structures. J. Mol. Biol. 1977, 112, 535–542. [Google Scholar] [CrossRef]
  26. Goto, J.; Kataoka, R.; Muta, H.; Hirayama, N. ASEDock-docking based on alpha spheres and excluded volumes. J. Chem. Inf. Model. 2008, 48, 583–590. [Google Scholar] [CrossRef] [PubMed]
  27. Corbeil, C.R.; Williams, C.I.; Labute, P.; Spraggs, C.F. Variability in docking success rates due to dataset preparation. J. Comput. Aided Mol. Des. 2012, 26, 775–786. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Chemical Computing Group Inc. Molecular Operating Environment (MOE), 2013.08; Chemical Computing Group Inc.: Montreal, QC, Canada, 2017. [Google Scholar]
  29. Trott, O.; Olson, A.J. AutoDock Vina: Improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. J. Comput. Chem. 2010, 31, 455–461. [Google Scholar] [CrossRef] [PubMed]
  30. Steere, A.C.; Klitz, W.; Drouin, E.E.; Falk, B.A.; Kwok, W.W.; Nepom, G.T.; Baxter-Lowe, L.A. Antibiotic-refractory Lyme arthritis is associated with HLA-DR molecules that bind a Borrelia burgdorferi peptide. J. Exp. Med. 2006, 203, 961–971. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Ijms 19 01660 g001 550
Figure 1. Chemical structure of nevirapine.
Figure 1. Chemical structure of nevirapine.
Ijms 19 01660 g001
Ijms 19 01660 g002 550
Figure 2. Binding mode of nevirapine to human leukocyte antigen (HLA)-DRB1*01:01 in docking simulations, (a) side view and (b) top view. The structures of HLA-DRB1*01:01, nevirapine, and amino acid residues of HLA-DRB1*01:01 in the vicinity of nevirapine are depicted in cartoon mode, ball-and-stick model, and wire model, respectively.
Figure 2. Binding mode of nevirapine to human leukocyte antigen (HLA)-DRB1*01:01 in docking simulations, (a) side view and (b) top view. The structures of HLA-DRB1*01:01, nevirapine, and amino acid residues of HLA-DRB1*01:01 in the vicinity of nevirapine are depicted in cartoon mode, ball-and-stick model, and wire model, respectively.
Ijms 19 01660 g002
Ijms 19 01660 g003 550
Figure 3. Parameter profiles of molecular dynamics (MD) simulations of nevirapine-bound HLA-DRB1*01:01. (a) Calculated energies vs. time plot; (b) RMSF values of polypeptide backbone. The location of α and β chains and residues that comprise the peptide binding groove helices are indicated by the solid and dashed lines running just above the x-axis, respectively. The α chain (solid black), α chain helix (solid purple), β chain (dashed black), and β chain helix (dashed purple); (c) Root mean square deviation (RMSD) values of polypeptide backbone vs. time plot; and (d) The average distance between each Cα in the helix of the α chain and the closest Cα atom in the helix of the β chain.
Figure 3. Parameter profiles of molecular dynamics (MD) simulations of nevirapine-bound HLA-DRB1*01:01. (a) Calculated energies vs. time plot; (b) RMSF values of polypeptide backbone. The location of α and β chains and residues that comprise the peptide binding groove helices are indicated by the solid and dashed lines running just above the x-axis, respectively. The α chain (solid black), α chain helix (solid purple), β chain (dashed black), and β chain helix (dashed purple); (c) Root mean square deviation (RMSD) values of polypeptide backbone vs. time plot; and (d) The average distance between each Cα in the helix of the α chain and the closest Cα atom in the helix of the β chain.
Ijms 19 01660 g003
Ijms 19 01660 g004 550
Figure 4. Simulated representative structures of nevirapine-bound HLA-DRB1*01:01 (a–c), hemagglutinin (HA) peptide-bound HLA-DRB1*01:01 in frame 2 (d–f) and frame 3 (g–i), and nevirapine-bound HLA-DRB1*01:01-HA peptide in frame 2 (j–l) and frame 3 (m–o). (a,d,g,j,m): Alignment of initial structure (red) and representative structure (blue); (b,e,h,k,n): Sausage plot of the structure, where the color and thickness of HLA-DRB1*01:01 are proportional to the RMSF of Cα; and (c,f,i,l,o): The volume occupied by nevirapine and HA peptide. The blue envelope outlines the 50% occupancy volume of nevirapine (red) and HA peptide (green).
Figure 4. Simulated representative structures of nevirapine-bound HLA-DRB1*01:01 (a–c), hemagglutinin (HA) peptide-bound HLA-DRB1*01:01 in frame 2 (d–f) and frame 3 (g–i), and nevirapine-bound HLA-DRB1*01:01-HA peptide in frame 2 (j–l) and frame 3 (m–o). (a,d,g,j,m): Alignment of initial structure (red) and representative structure (blue); (b,e,h,k,n): Sausage plot of the structure, where the color and thickness of HLA-DRB1*01:01 are proportional to the RMSF of Cα; and (c,f,i,l,o): The volume occupied by nevirapine and HA peptide. The blue envelope outlines the 50% occupancy volume of nevirapine (red) and HA peptide (green).
Ijms 19 01660 g004
Ijms 19 01660 g005 550
Figure 5. Drug-induced peptide displacements in (a) nevirapine-bound HLA-DRB1*01:01-HA peptide in frame 2 and (b) frame 3; (c) lapatinib-bound HLA-DRB1*01:01-HA peptide in frame 2 (d) and frame 3; (e) lapatinib-bound HLA-DRB1*07:01-peptide from tetanus toxoid (TT) in frame 3; and (f) ximelagatran-bound HLA-DRB1*07:01-TT peptide in frame 3. Each residue in the ligand peptides is colored according to the displacement of the Cα in the trimer simulations, relative to its position in the HLA-peptide dimer simulations. The same color scale is used and the range is from 0 (blue) to 10 Å (red). Small values mean that the drug has little effect on the position of ligand peptides in the binding groove. Drugs are shown using grey stick representation.
Figure 5. Drug-induced peptide displacements in (a) nevirapine-bound HLA-DRB1*01:01-HA peptide in frame 2 and (b) frame 3; (c) lapatinib-bound HLA-DRB1*01:01-HA peptide in frame 2 (d) and frame 3; (e) lapatinib-bound HLA-DRB1*07:01-peptide from tetanus toxoid (TT) in frame 3; and (f) ximelagatran-bound HLA-DRB1*07:01-TT peptide in frame 3. Each residue in the ligand peptides is colored according to the displacement of the Cα in the trimer simulations, relative to its position in the HLA-peptide dimer simulations. The same color scale is used and the range is from 0 (blue) to 10 Å (red). Small values mean that the drug has little effect on the position of ligand peptides in the binding groove. Drugs are shown using grey stick representation.
Ijms 19 01660 g005
Table
Table 1. Mean molecular dynamics (MD) simulations parameters.
Table 1. Mean molecular dynamics (MD) simulations parameters.
Hemagglutinin PeptideNevirapineEnergies (kcal/mol)RMSF (Å)Inter-Helical Distance (Å)Slope of Inter-Helical Distance (Å/ns)
-−−106715 ± 309.91.0 ± 0.413.0 ± 0.4−0.12
-+−106370 ± 533.50.8 ± 0.313.1 ± 0.4−0.14
Frame 2−−111392 ± 579.80.9 ± 0.412.7 ± 0.2−0.16
Frame 3−−106078 ± 323.10.9 ± 0.314.2 ± 0.30.04
Frame 2+−120550 ± 243.60.8 ± 0.313.5 ± 0.20.02
Frame 3+−112166 ± 230.50.9 ± 0.312.6 ± 0.2−0.13
The total energies, average root mean square fluctuation (RMSF), distance between the α-helices flanking the binding groove, and slope of inter-helical distance curve are listed. Values represent the mean ± standard deviation (SD) during 5-ns simulations.
Table
Table 2. The effect of nevirapine on the binding of ligand peptides to HLA-DR molecules.
Table 2. The effect of nevirapine on the binding of ligand peptides to HLA-DR molecules.
HLA AlleleDRB1*01:01DRB1*07:01DRB1*15:01
Concentration of nevirapine (μM)1000263.9 ± 15.5 #62.6 ± 3.9101.3 ± 14.9
20099.7 ± 11.3110.4 ± 5.594.6 ± 10.6
40106.2 ± 6.576.9 ± 4.693 ± 10.4
8101.9 ± 8.1132.2 ± 6498.9 ± 3.4
The effects of nevirapine on the binding of ligand peptides to three HLA-DRs are expressed as % binding of ligand peptides in the presence of nevirapine to that in dimethyl sulfoxide (DMSO) control (n = 8). Values represent the mean ± SD of quadruplicate. # p < 0.001, compared with the DMSO control.

Share and Cite

MDPI and ACS Style

Hirasawa, M.; Hagihara, K.; Abe, K.; Ando, O.; Hirayama, N. Interaction of Nevirapine with the Peptide Binding Groove of HLA-DRB1*01:01 and Its Effect on the Conformation of HLA-Peptide Complex. Int. J. Mol. Sci. 2018, 19, 1660. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms19061660

AMA Style

Hirasawa M, Hagihara K, Abe K, Ando O, Hirayama N. Interaction of Nevirapine with the Peptide Binding Groove of HLA-DRB1*01:01 and Its Effect on the Conformation of HLA-Peptide Complex. International Journal of Molecular Sciences. 2018; 19(6):1660. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms19061660

Chicago/Turabian Style

Hirasawa, Makoto, Katsunobu Hagihara, Koji Abe, Osamu Ando, and Noriaki Hirayama. 2018. "Interaction of Nevirapine with the Peptide Binding Groove of HLA-DRB1*01:01 and Its Effect on the Conformation of HLA-Peptide Complex" International Journal of Molecular Sciences 19, no. 6: 1660. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms19061660

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop