Abstract
Drug-drug interactions (DDIs) via cytochrome P450 (P450) induction are one clinical problem leading to increased risk of adverse effects and the need for dosage adjustments and additional therapeutic monitoring. In silico models for predicting P450 induction are useful for avoiding DDI risk. In this study, we have established regression models for CYP3A4 and CYP2B6 induction in human hepatocytes using several physicochemical parameters for a set of azole compounds with different P450 induction as characteristics as model compounds. To obtain a well-correlated regression model, the compounds for CYP3A4 or CYP2B6 induction were independently selected from the tested azole compounds using principal component analysis with fold-induction data. Both of the multiple linear regression models obtained for CYP3A4 and CYP2B6 induction are represented by different sets of physicochemical parameters. The adjusted coefficients of determination for these models were of 0.8 and 0.9, respectively. The fold-induction of the validation compounds, another set of 12 azole–containing compounds, were predicted within twofold limits for both CYP3A4 and CYP2B6. The concordance for the prediction of CYP3A4 induction was 87% with another validation set, 23 marketed drugs. However, the prediction of CYP2B6 induction tended to be overestimated for these marketed drugs. The regression models show that lipophilicity mostly contributes to CYP3A4 induction, whereas not only the lipophilicity but also the molecular polarity is important for CYP2B6 induction. Our regression models, especially that for CYP3A4 induction, might provide useful methods to avoid potent CYP3A4 or CYP2B6 inducers during the lead optimization stage without performing induction assays in human hepatocytes.
Introduction
Drug-drug interactions (DDIs) can affect the pharmacokinetics of a coadministered drug or its metabolites, which consequently increases the risk of adverse effects from the drug and causes the need for dosage adjustment and additional therapeutic monitoring. Therefore, it is important to predict DDIs during the drug discovery process. In the past several years, the pharmaceutical regulatory agencies in the United States, Europe, and Japan have issued guideline/draft guidance documents for in vitro and in vivo DDI studies that need to be conducted during the development of new drug candidates (http://www.fda.gov/downloads/Drugs/GuidanceComplianceRegulatoryInformation/Guidances/UCM292362.pdf, http://www.ema.europa.eu/docs/en_GB/document_library/Scientific_guideline/2012/07/WC500129606.pdf, http://www.nihs.go.jp/mhlw/20131488.pdf). To reduce DDI risks, several screening methods are conducted to design safer and more effective drug molecules.
The inhibition and/or induction of cytochrome P450s (P450s) are the main mechanisms of DDIs. Several clinically used drugs cause P450 induction, including the antibiotic rifampicin (Niemi et al., 2003), the acyl-CoA/cholesterol acyltransferase inhibitor avasimibe (Sahi et al., 2003), and the dual endothelin receptor antagonist bosentan (Dingemanse and van Giersbergen, 2004). These drugs reduce the plasma concentration of midazolam and/or those of themselves. General induction mechanisms involve the activation of nuclear receptor-type transcription factors. Inducers directly bind to or indirectly activate one or more of the following receptors: aryl hydrocarbon receptor, pregnane X receptor (PXR), and constitutive androstane receptor (CAR). Among them, PXR plays a key role in the drug-mediated induction of CYP3A4, a major human drug metabolizing P450 (LeCluyse, 2001). Therefore, PXR activation is a good prediction marker for CYP3A4 induction (Luo et al., 2002). However, many of the ligands that activate PXR can also activate CAR, and there is an overlap in the target genes of these two receptors (Maglich et al., 2002, 2003; Wei et al., 2002), and there are uncertainties about CYP3A4 induction with CAR activators and/or PXR/CAR dual agonists because there is no perfect reporter assay system of CAR (Imai, et al., 2013).
Although several prediction systems have been developed for the substrates and inhibitors of P450 enzymes (Wanchana et al., 2003; Mishra et al., 2010), systems for P450 induction are limited. As an in silico approach to predict enzyme induction, computational ligand docking, quantitative structure-activity relationship approach using machine-learning methods, and classifications using recursive partitioning methods have been developed for PXR that plays a key role in P450 induction (Khandelwal et al., 2008; Xiao et al., 2011; Dybdahl et al., 2012; Wu et al., 2013; Handa et al., 2015). Docking studies provide important interactions between ligands and proteins but have limitations in quantitative predictions. Although general quantitative structure-activity relationship models can provide a quantitative prediction, they require a large number of agonists to establish a reliable model. Therefore, a quantitative prediction model conducted with a moderate number of compounds and commonly used physicochemical parameters would be a useful tool in lead-optimizing and developing processes.
Conazoles are a class of N-substituted azole antifungal drugs, which include two major classes, namely, imidazole- and triazole-containing compounds. These azole-containing antifungals are also used as environmental agents to control fungal growth. Many conazoles inhibit ergosterol synthesis through the inhibition of fungal CYP51 (lanosterol 14α-demethylase) activity. Furthermore, these conazoles inhibit mammalian hepatic P450s in the clinic (Niwa et al., 2014), especially because imidazole antimycotics are known to show extremely high affinity for the heme iron atom of P450. In addition, a number of azole-containing compounds, including cyproconazole, fenbuconazole, propiconazole, fluconazole, triadimefon, and myclobutanil, have been reported with abilities to induce CYP1A, CYP2B, CYP2C, and CYP3A enzymes in mouse livers (Goetz et al., 2006; Juberg et al., 2006; Sun et al., 2006, 2007, Tamura et al. 2013). Voriconazole, a triazole-containing conazole, shows PXR- and CAR-mediated autoinduction in mice (Ohbuchi et al., 2013). Clotrimazole, also known as a CYP3A4 inducer, and its analogs were potent ligands of human PXR (Matsuura et al., 1991; Luo et al., 2002; Sahi et al., 2009). Therefore, the azole-containing compounds might be good model compounds for a P450 induction prediction model due to their variable abilities to induce P450s.
In this study, we used a group of azole compounds for a learning set and established regression models with several physicochemical parameters for CYP3A4 and CYP2B6 induction in human hepatocytes. The accuracy and limitations of the models were clarified using another set of azole compounds and non-azole compounds (commercially used pharmaceutical drugs).
Materials and Methods
Materials.
Climbazole, hexaconazole, triflumizole, and uniconazole were purchased from Abcam Biochemicals (Tokyo, Japan). Fenbuconazole, imazalil, and lanoconazole were obtained from Santa Cruz Biotechnology (Santa Cruz, CA), AK Scientific (Union City, CA), and Sequoia Research Products (Pangbourne, UK), respectively. Bifonazole, clotrimazole, econazole, miconazole, avasimibe, flumazenil, nifedipine, rifampicin, 3-methylcholanthrene (3-MC), and 6-(4-chlorophenyl)imidazo[2,1-b][1,3]thiazole-5-carbaldehyde O-(3,4-dichlorobenzyl)oxime (CITCO) were from Sigma-Aldrich (St. Louis, MO). Fluconazole, isoconazole, voriconazole, ethinyl estradiol, fluoxetine, ranitidine, and topiramate were purchased from Tokyo Chemical Industry (Tokyo, Japan). Myclobutanil, oxiconazole, propiconazole, tebuconazole, tioconazole, carbamazepine, gatifloxacin, leflunomide, nafcillin, phenobarbital, phenytoin, pioglitazone, and a metabolite of miconazole (1-(2,4-dichlorophenyl)-2-(1H-imidazol-1-yl)ethanol) were obtained from Wako Pure Chemical Industries (Osaka, Japan). Luliconazole, a metabolite of voriconazole [1-(2,4-difluorophenyl)-2-(1H-1,2,4-triazole-1-yl)ethanol], and the validation compounds containing imidazole moieties (WO 2008/156092) were produced at Kaken Pharmaceutical (Kyoto, Japan). Aprepitant, bosentan, efavirenz, fluvoxamine, pleconaril, probenecid, and troglitazone were from Axon Medchem (Reston, VA), Pharmten Chemical (Leshan, China), Tronto Research Chemicals (Tronto, Canada), Tocris Bioscience (Bristol, UK), Cheminstock (Shanghai, China), Nacalai Tesque (Kyoto, Japan), and Cayman Chemical (Ann Arbor, MI), respectively. Nevirapine and rosiglitazone were obtained from Chem Pacific (Zhejiang, China). Cryopreserved human hepatocytes, hepatocyte thawing, plating, incubation media, antibiotics, and Geltrex were purchased from Life Technologies (Carlsbad, CA). All other chemicals and reagents were of reagent grade or better and were obtained from commercial sources.
Induction Experiments Using Cryopreserved Human Hepatocytes.
Cryopreserved human hepatocytes (lot no. Hu1423, Caucasian, male, 61 years old, and Hu8125, Caucasian, male, 44 years old) were thawed with thawing medium, according to the manufacturer’s protocol. The cells were seeded in collagen I–precoated 96-well plates (Life Technologies) at a density of 0.5 × 105 viable cells/well. Viability was determined by trypan blue exclusion of 90% or better for this study. Cells were maintained in fetal bovine serum (Life Technologies)–added Williams’ E medium (Life Technologies) at 37°C in a humidified incubator with an atmosphere of 5% CO2–95% air for 4–6 hours. Subsequently, the medium was changed to serum-free Williams’ E medium (incubation medium) containing 0.35 mg/mL Geltrex. The following day, the medium was changed to fresh incubation medium, and the hepatocytes were incubated at 37°C. The medium was removed 24 hours after plating, and the hepatocytes were treated with medium containing the test compounds for 48 hours. Fresh medium containing the test compounds was supplied daily. The test compounds were used at a single concentration of 10 μM, and the reference compounds, 3-MC (1 μM) for CYP1A2, phenobarbital (750 μM) and CITCO (1 μM) for CYP2B6, rifampicin (10 μM) for CYP3A4, and probenecid (10 μM) as a negative control, were used to confirm the performance of hepatocyte preparations. Both test compounds and reference inducers were dissolved in dimethyl sulfoxide, and the final concentration of dimethyl sulfoxide in the incubation medium was 0.1%.
mRNA Level Determination.
A QuantiGene Plex 2.0 Assay Kit (Panel 12117) from Affymetrix (Santa Clara, CA) was used to determine mRNA levels. After a 48-hour incubation period, the hepatocytes were washed with Krebs-Henseleit buffer and then lysed. The mRNA level determination was performed according to the manufacturer’s protocol. CYP1A2, CYP2B6, CYP2C9, and CYP3A4 mRNA levels were normalized against the geometric mean value of two housekeeping genes: PGK1 and HPRT1.
Data Processing and Model Development.
The optimized three-dimensional conformations of the test compounds were generated using ChemBio3D Ultra 12.0 (PerkinElmer, Waltham, MA), and their physicochemical parameters were calculated with QikProp (Schrödinger, Mannheim, Germany) and clogP values were calculated with ChemBioDraw Ultra 12.0 (PerkinElmer). Regression models for P450 induction were established using the set of compounds shown in Fig. 1. Multiple linear regression modeling was performed using JMP 7 (SAS Institute Japan, Tokyo, Japan) with physicochemical parameters. The goodness of fit for the models was assessed using the adjusted coefficient of determination (R2) and the root mean square error (RMSE):(1)The selection of variables was carried out by a stepwise forward-backward selection. To narrow the test compounds down based on their P450 induction profiles, principal component analysis (PCA) was performed using JMP 7 with fold-induction values of CYP1A2, CYP2B6, CYP2C9, and CYP3A4, and a ratio of CYP3A4/CYP2B6.
Chemical structures of the test azole compounds.
Prediction of CYP3A4 and CYP2B6 Induction Using the Regression Models.
Two types of validation compound sets were used to assess the predictive applicability: another set of azole compounds and a set of non-azole compounds [commercially marketed drugs used in the report by Fahmi et al. (2010)]. For CYP3A4 induction of the marketed drugs, the predictive capability of the regression models established with the test compounds was assessed with the classification system, including sensitivity, specificity, and concordance, which were defined as follows:(2)
(3)
(4)To assess the applicability of the models for CYP3A4 and CYP2B6 induction using the different lots of hepatocytes, the fold-induction values of the test azole compounds in lot Hu8125 were predicted with the regression models established by lot Hu1423. Because the fold-induction values of non-inducers are always approximately 1, in any lots hepatocytes and the fold-induction values of inducers are usually different based on the characteristics of hepatocyte lots used, the potency of P450 induction was normalized using the fold-induction values of positive controls: rifampicin for CYP3A4 and CITCO for CYP2B6. The normalization was conducted using the following equation:
(5)x, y, A, and B were the normalized predicted fold-induction value in lot Hu1423, the predicted fold-induction value in lot Hu8125, the slope, and the intercept of the straight line used for the normalization, respectively.
The equation contained two points, (x, y) = (1, 1) and (an observed fold-induction value of positive control in lot Hu1423, an observed fold-induction value of positive control in lot Hu8125).
Results
P450 Inductions in Human Hepatocytes.
Most of the imidazole- or triazole-containing compounds (Fig. 1) increased the mRNA levels of CYP3A4 and CYP2B6, whereas only a few compounds, such as bifonazole and oxiconazole, increased the levels of CYP1A2 more strongly than CYP3A4 and CYP2B6 (Table 1). The results indicated that climbazole, etc., increased CYP3A4 mRNA levels more than CYP2B6 mRNA levels as rifampicin, and that lanoconazole, etc., increased CYP2B6 mRNA levels more than CYP3A4 mRNA levels as CITCO. The 3-MC showed high selectivity to CYP1A2, and probenecid, which was used as a negative control, increased the CYP1A2, CYP2B6, CYP2C9, and CYP3A4 mRNA levels by less than twofold (Table 1).
Fold-induction of CYP1A2, CYP2B6, CYP2C9, and CYP3A4 for test azole and reference compounds
Regression Models for CYP3A4 and CYP2B6 Induction.
The physicochemical parameters of the test compounds calculated are shown in Table 2 and Supplemental Table 1. The CYP3A4 and CYP2B6 fold-induction values showed no significant correlation with clogP or mol. wt., which are commonly used molecular descriptors (Supplemental Fig. 1). Therefore, the regression models using all of the test azole compounds did not predict the CYP3A4 and CYP2B6 induction with enough accuracy (R2 = 0.4984 and RMSE = 2.753 for CYP3A4, and R2 = 0.6098 and RMSE = 1.709 for CYP2B6) (Supplemental Fig. 2).
Physicochemical parameters of the test azole compounds
To improve the regression models, we narrowed the list of test compounds using PCA with induction data. The loading plots and score plots are shown in Fig. 2. In the score plots, the preferential CYP3A4 inducers were plotted in the shadow area of PC1 versus PC2 (Fig. 2C), and CYP2B6 inducers were plotted in the shadow area of PC1 versus PC3 (Fig. 2D). Therefore, these compounds in the shadow areas were selected and subjected to multiple linear regression modeling for each CYP3A4 and CYP2B6 induction. The physicochemical parameters used in the regression model were selected by stepwise forward-backward methods within the criteria of P values less than 0.05. The obtained regression models for CYP3A4 and CYP2B6 are as follows:
Loading and score plots of the test azole compounds from PCA. (A) Loading plots of PC1 versus PC2. (B) Loading plots of PC1 versus PC3. (C) Score plots of PC1 versus PC2. (D) Score plots of PC1 versus PC3. Compounds in shadow areas were considered preferential CYP3A4 inducers (C) and potential CYP2B6 inducers (D).


With these regression models, fold-induction values were well predicted, as shown in Fig. 3 (R2 = 0.8058 and RMSE = 2.092 for CYP3A4, and R2 = 0.9219 and RMSE = 0.890 for CYP2B6). In the regression models, lipophilicity is a main contributor to CYP3A4 induction, and not only lipophilicity but also molecular polarity plays important roles for CYP2B6 induction.
Correlation of the predicted and observed fold-induction values of CYP3A4 (A) and CYP2B6 (B) induction by selected compounds. Correlation analyses were performed on the compounds selected by PCA in Fig. 2. The solid line in each graph is the line of unity, and the dotted lines indicate a twofold limit. The observed fold-induction is shown as the mean ± SD (n = 3–6; n refers to the number of wells in a single hepatocyte experiment).
Verification of the Regression Models.
To verify the regression models, induction assays with cryopreserved human hepatocytes were performed using another set of imidazole-containing compounds (validation compounds: Supplemental Fig. 3). The physicochemical parameters of the validation compounds were obtained with QikProp and ChemBioDraw (Supplemental Table 2; Table 3), and the fold-induction values of the validation compounds were predicted using the established regression models for CYP3A4 and CYP2B6 induction. The predicted fold-induction values for both CYP3A4 and CYP2B6 were within twofold of the fold-inductions observed (Fig. 4; Table 4).
Physicochemical parameters and predicted and observed fold-induction values of azole-containing validation compounds
Correlation between the predicted and observed fold-induction values of CYP3A4 (A) and CYP2B6 (B) induction for the azole-containing validation compounds. The solid line in each graph is the line of unity, and the dotted lines indicate a twofold limit. The observed fold-induction is shown as the mean ± SD (n = 3; n refers to the number of wells in a single hepatocyte experiment).
Physicochemical parameters and predicted and observed fold-induction values of marketed drugs
Applicability of the Regression Models.
To know the validity of the regression models for CYP3A4 and CYP2B6 induction, in vitro induction assays were performed using 23 marketed drugs, which were previously reported by Fahmi et al. (2010) (Table 4). CYP3A4 induction by these drugs was predicted with the model, and the predicted values were within a twofold limit for 17 of the marketed drugs (73.9%) (Fig. 5; Table 4). Because various cutoff criteria were reported to define the inducers versus non-inducers with in vitro induction assays, the accuracy of the prediction for CYP3A4 induction with the models was evaluated using two cutoff categories, as follows: 1) a fold-induction of CYP3A4 mRNA levels (two- and fourfold) and 2) a fold induction achieving 10% and 20% that for rifampicin (10 μM) (Table 5). The best concordance was achieved when 10% of rifampicin was used as a cutoff criterion. Unlike CYP3A4 induction, CYP2B6 induction by these marketed drugs was not well predicted with the established models, and in most cases the fold-induction was overestimated (Fig. 5; Supplemental Table 3; Table 4).
Correlation between the predicted and observed fold-induction values of CYP3A4 (A) and CYP2B6 (B) induction for the marketed drugs. The solid line in each graph is the line of unity, and the dotted lines indicate a twofold limit. When the predicted fold-induction values were under 1, they were plotted at 1. The observed fold-induction values from one experiment are shown.
Summary of the prediction for CYP3A4 induction using marketed drugs
Finally, the applicability of the CYP3A4 and CYP2B6 regression models in different lots of hepatocytes was investigated. CYP3A4 and CYP2B6 induction by the 22 test azole compounds was predicted with the models, and the induction of CYP3A4 and CYP2B6 by these azole compounds was assessed in another lot of hepatocytes. The predicted values of CYP3A4 and CYP2B6 induction were within a twofold limit for 18 (81.8%) and 19 (86.4%) compounds, respectively (Fig. 6).
Correlation between the predicted and observed fold-induction values of CYP3A4 (A) and CYP2B6 (B) induction for the test azole compounds in a different lot of cryopreserved human hepatocytes. The solid line in each graph is the line of unity, and the dotted lines indicate a twofold limit. When the predicted fold-induction values were under 1, they were plotted at 1. The observed fold induction is shown as the mean ± SD (n = 3–6; n refers to the number of wells in a single hepatocyte experiment).
Discussion
In this study, we established the regression models for CYP3A4 and CYP2B6 inductions in human hepatocytes using physicochemical parameters calculated by commonly used software. As shown in Fig. 3, the predicted fold-induction values were strongly correlated with the observed values for both CYP3A4 and CYP2B6 when the multiple regression modeling was conducted using the compounds selected by PCA. Although the test compounds used in this study were congeneric compounds, the regression model for CYP3A4 induction could be applicable to structurally diverse compounds.
Because P450 induction is observed in cultured primary human hepatocytes and various methods have been reported to predict the clinical CYP3A4 induction from in vitro studies (Fahmi et al., 2008, 2009; Kaneko et al., 2009), the induction assays using human hepatocytes and/or immortalized cell lines are usually performed in drug discovery and/or development processes to estimate the induction risk of new drug candidates (Fahmi et al., 2010). In addition, the evaluation in human hepatocytes is recognized as the most accepted approach by regulatory authorities, such as the U.S. Food and Drug Administration. However, because of the physicochemical properties (e.g., insufficient solubility) and/or cytotoxicity of test compounds, it is difficult to assess P450 induction for all of the compounds using human hepatocytes in early drug discovery and lead-optimizing stages. Therefore, the prediction of P450 induction used in silico models is useful and desirable to evaluate the DDI risks of such compounds as well as those with synthetic difficulties, such as metabolites. Because the prediction models established in this study for P450 induction consist of only a few physicochemical parameters calculated with commercially available software, our models enable us to determine the potent inducers within a large number of compounds.
The multiple regression models first obtained using all of the azole compounds were unable to sufficiently predict CYP3A4 and CYP2B6 induction. Many ligands that activate PXR can also activate CAR, and there is overlap in the target genes of these two receptors (Maglich et al., 2002, 2003). In addition, the test azole compounds showed different abilities of inducing several P450s, suggesting that these azole compounds might induce P450 enzymes through multiple mechanisms. Therefore, we selected the test azole compounds to establish good prediction models for P450 induction using PCA, which is commonly used to characterize and categorize data sets. The PCA results suggested that PC1 might represent the potency of CYP3A4 induction, PC2 might show the selectivity for CYP1A2 induction, and PC3 might be the selectivity for CYP2B6 induction. Based on these, predominant CYP3A4 inducers and potent CYP2B6 inducers were selected for modeling. This process greatly improved the accuracy of the regression models for both CYP3A4 and CYP2B6 induction; the adjusted coefficients of determination (R2) were 0.8 and 0.9, respectively.
The variability in P450 induction among different batches of hepatocytes is normally great. However, our regression models showed that more than 80% of compounds could be predicted within twofold limited errors both CYP3A4 and CYP2B6 induction in different lots of cryopreserved primary hepatocytes after normalization with the fold-induction values of positive control compounds. Therefore, our regression models might be generally used to predict the potency of P450 induction independently of hepatocyte lots.
Our final regression model for CYP3A4 induction contains lipophilicity as the main determinant. In fact, the validation compounds E, F, and G have low clogP values compared with compound D, and the predicted and observed fold-induction values of E, F, and G were lower than compound D. Conversely, compounds H, I, J, K, and L, which have high clogP values compared with compound D, showed higher fold-induction values than compound D. Similar results were obtained for non-azole compounds; pleconaril, a known CYP3A4 inducer, clearly increased CYP3A4 mRNA levels. Several marketed drugs, including aprepitant, avasimibe, bosentan, efavirenz, terbinafine, and troglitazone, having high clogP values increased the CYP3A4 mRNA levels much more than pleconaril did. A chlorine atom is often used to increase the lipophilicity and molecular size in the lead-optimizing process. Econazole has a similar structure to isoconazole and miconazole but contains fewer chlorine atoms, suggesting that CYP3A4 induction by econazole is weaker than those by isoconazole and miconazole. The results of the induction assays demonstrated that this prediction was correct. Considering the induction mechanism, Yoshida et al. (2012) reported that PXR agonists tend to be more lipophilic and larger in molecular size than non-agonists. In addition, Ung et al. (2007) reported that PXR activators contain more halogen atoms, especially chlorine atoms, than nonactivators. In summary, our regression model for CYP3A4 induction agrees with these reports, which demonstrates the reliability and applicability of our model.
In contrast to the prediction of CYP3A4 induction, our regression model for CYP2B6 induction seems to have structural limitations. For compounds containing imidazole and triazole moieties, the predicted CYP2b6 fold-induction values were within twofold of the observed values, whereas the predicted values of the marketed drugs, non-azole compounds, tended to be overestimated. Currently, the majority of identified CYP3A4 and CYP2B6 inducers are known activators of PXR, but not CAR (Faucette et al., 2004). Only a limited number of compounds, including CITCO and phenytoin, are shown to induce CYP3A4 and/or CYP2B6 preferentially through CAR (Maglich et al., 2003; Wang et al., 2004), with stronger induction of CYP2B6 than CYP3A4. For 52% of the marketed drugs used, the fold-induction ratio of CYP3A4/CYP2B6 was higher than that of phenobarbital (CYP3A4/CYP2B6: 1.37) (data not shown). Therefore, most of the marketed drugs used in this study might be preferential CYP3A4 inducers through PXR activation.
In addition, CAR is activated through two different mechanisms, as follows: 1) activation by direct binding of agonists, such as CITCO for human CAR (Maglich et al., 2003), and 2) indirect activation through phosphorylation/dephosphorylation signals (Sueyoshi and Negishi, 2001; Mutoh et al., 2013). Phenobarbital, which is commonly used as a positive control for CYP2B6 induction, activates CAR through the latter mechanism (Mutoh et al., 2013). Generally, higher concentrations are needed to achieve significant CYP2B6 induction for these indirect activators compared with direct activators or ligands. Therefore, the reason for the poor applicability of the regression model for CYP2B6 induction might result from the use of fold-induction values at a fixed concentration of 10 μM in this study. In fact, when we used Emax values of CYP2B6 induction, the maximum fold-induction value observed in the concentration range used, instead of the fold-induction values at 10 μM as the observed values for the marketed drugs, the correlation between the predicted and observed values improved; the prediction of 14 compounds (60%) was within a twofold limit (data not shown). These results suggest that the marketed drugs used in this study need higher concentrations to achieve Emax compared with the azole-containing compounds. Thus, the regression model for CYP2B6 induction established by azole compounds probably has structural and/or mechanistic limitations.
Conazoles are widely used as antifungal drugs, although their potent abilities of P450 inhibition and induction limit their clinical application by systematic treatment. First-generation conazoles, such as miconazole, econazole, and isoconazole, are mainly used to treat superficial mycoses, and second/third-generation conazoles, such as fluconazole and voriconazole, are applied to treat emerging invasive fungal infections (Heeres et al., 2010). Because there is a need for more powerful and easy-to-use antifungal drugs for invasive infections, some conazoles, such as albaconazole, ravuconazole, isavuconazole, and pramiconazole, have been developed and/or marketed. Our regression model could aid in developing a conazole with an advantage in pharmacokinetics and DDIs showing no or less potency for CYP3A4 and CYP2B6 induction.
In conclusion, the quantitative regression models for CYP3A4 and CYP2B6 induction were established with a few physicochemical parameters, even though there were limitations with the CYP2B6 induction model. Because CYP3A enzymes, including CYP3A4, account for 30% of the P450s in human livers and mediate the metabolism of more than 60% of marketed drugs (di Masi et al., 2009), the CYP3A4-associated DDIs should be avoided at first. Our regression models for P450 induction provide a useful and promising method to obtain compounds with less DDI risk during the drug-screening process without chemical synthesis.
Authorship Contributions
Participated in research design: Nagai, Konno, Satsukawa, Yamashita, Yoshinari.
Conducted experiments: Nagai, Yoshinari.
Performed data analysis: Nagai, Yoshinari.
Wrote or contributed to the writing of the manuscript: Nagai, Yoshinari.
Footnotes
- Received December 19, 2015.
- Accepted May 18, 2016.
↵
This article has supplemental material available at dmd.aspetjournals.org.
Abbreviations
- CAR
- constitutive androstane receptor
- CITCO
- 6-(4-chlorophenyl)imidazo[2,1-b][1,3]thiazole-5-carbaldehyde O-(3,4-dichlorobenzyl)oxime
- DDI
- drug-drug interaction
- 3-MC
- 3-methylcholanthrene
- PCA
- principal component analysis
- P450
- cytochrome P450
- PXR
- pregnane X receptor
- RMSE
- root mean square error
- Copyright © 2016 by The American Society for Pharmacology and Experimental Therapeutics