Impact Factor 3.258 | CiteScore 2.7
More on impact ›

Original Research ARTICLE

Front. Genet., 21 August 2020 | https://doi.org/10.3389/fgene.2020.00944

Inherited and Acquired Determinants of Hepatic CYP3A Activity in Humans

Johannes Matthaei1, Wagner Hugo Bonat2, Reinhold Kerb3, Mladen Vassilev Tzvetkov1, Jakob Strube1, Stefanie Brunke1, Cordula Sachse-Seeboth1, Daniel Sehrt1, Ute Hofmann3, Jacob von Bornemann Hjelmborg2, Matthias Schwab3,4,5 and Jürgen Brockmöller1*
  • 1Institute for Clinical Pharmacology, University Medical Center Göttingen, Georg-August University, Göttingen, Germany
  • 2Department of Epidemiology, Biostatistics and Biodemography, University of Southern Denmark, Odense, Denmark
  • 3Dr. Margarete Fischer-Bosch Institute of Clinical Pharmacology and University of Tübingen, Stuttgart, Germany
  • 4Department of Clinical Pharmacology, University Hospital Tübingen, Tübingen, Germany
  • 5Department of Pharmacy and Biochemistry, University of Tübingen, Tübingen, Germany

Human CYP3A enzymes (including CYP3A4 and CYP4A5) metabolize about 40% of all drugs and numerous other environmental and endogenous substances. CYP3A activity is highly variable within and between humans. As a consequence, therapy with standard doses often results in too low or too high blood and tissue concentrations resulting in therapeutic failure or dose-related adverse reactions. It is an unanswered question how much of the big interindividual variation in CYP3A activity is caused by genetic or by environmental factors. This question can be answered by the twin study approach. Using midazolam as CYP3A probe drug, we studied 43 monozygotic and 14 dizygotic twins and measured midazolam and its metabolite 1-OH-midazolam. In addition, endogenous biomarkers of CYP3A activity, 4ß-OH-cholesterol and 6ß-OH-cortisol, were analyzed. Additive genetic effects accounted for only 15% of the variation in midazolam AUC, whereas 48% was attributed to common environmental factors. In contrast, 73, 56, and 31% of 1-OH-midazolam, 4ß-OH-cholesterol and 6ß-OH-cortisol variation was due to genetic effects. There was a low phenotypic correlation between the four CYP3A biomarkers. Only between midazolam and its 1-OH-metabolite, and between midazolam and 6ß-OH-cortisol we found significant bivariate genetic correlations. Midazolam AUC differed depending on the CYP3A422 variant (p = 0.001) whereas plasma 4ß-OH-cholesterol was significantly lower in homozygous carriers of CYP3A53 (p = 0.02). Apparently, non-genomic factors played a dominant role in the inter-individual variation of the CYP3A probe drug midazolam. A small intra-individual pharmacokinetic variation after repeated administration of midazolam was rated earlier as indication of high heritability of CYP3A activity, but according to present data that could also largely be due to constant environmental factors and/or heritability of liver blood flow. The higher heritabilities of 4ß-OH-cholesterol and of 1-OH-midazolam may deserve further research on the underlying factors beyond CYP3A genes.

Clinical Trial Registration: ClinicalTrials.gov: NCT01845194 and EUDRA-CT: 2008-006223-31.

Introduction

The cytochrome P450 (CYP) enzymes 3A4 and 3A5 play a key role in the biotransformation of about 40% of currently used drugs (Evans and Relling, 1999). CYP3A enzymes also catalyze the biotransformation of numerous exogenous and endogenous substances, notably cholesterol, bile acids, and many steroid hormones (Patki et al., 2003; Bodin et al., 2005; Klein and Zanger, 2013). CYP3A4 is the most prominent CYP3A subtype in adult humans. The homologous isozyme CYP3A5 has partially overlapping substrate specificity with CYP3A4 but may significantly contribute to CYP3A metabolic activity with some drugs like tacrolimus. Only about 10–15% of persons with European ancestry express CYP3A5 and that can be easily analyzed by genotyping. However, only a small fraction of the variation in CYP3A4 expression and activity can thus far be explained by defined genotypes.

Expression and activity of human CYP3A enzymes in the gut epithelia and liver cells is highly variable both inter- as well as intra-individually (Koch et al., 2002). CYP3A activity may be reduced to 10% by strong inhibitors like itraconazole, ritonavir, or grapefruit juice (Kim et al., 2006; Greenblatt and Harmatz, 2015). It could also be increased up to 10-fold by strong transcriptional inducers like rifampin, carbamazepin or phenobarbital (Bodin et al., 2001; Bui et al., 2016). However, even when excluding such strong environmental modulators of CYP3A activity and expression, there remains a very wide and mostly unexplained variation between individuals (Gorski et al., 1998; Gashaw et al., 2003). By extensive search for genetic polymorphisms in the CYP3A gene locus that could account for inter-individual variation in CYP3A activity only a few variants were found that explained overall only a small fraction of that variation (Klein and Zanger, 2013). Also, polymorphisms in nuclear transcription factors and other regulators of expression appear to contribute to the variation to a minor extent only. A better understanding of the inter-individual variation in CYP3A activity may result in molecular biomarkers that could facilitate individually more precise dosing of drugs metabolized by CYP3A.

Our understanding of CYP3A variation depends on how the CYP3A phenotype was measured. As human liver tissue or intestinal tissue is mostly not available in clinical research, in vivo probe drugs are the most important tools for the study of individual CYP3A activity. Numerous CYP3A probe drugs have been applied. Currently, midazolam is the most widely used (Thummel et al., 1994; Ozdemir et al., 2000; Tateishi et al., 2001; Eap et al., 2004a, b; Kirby et al., 2006; Fromm et al., 2007; Fuhr et al., 2007; Tomalik-Scharte et al., 2014), mostly because of its high sensitivity in reflecting even moderate enzyme inhibition and induction. Approximately 75% of midazolam is eliminated as 1-OH-midazolam (Heizmann et al., 1983), and this metabolite is formed by CYP3A (Kronbach et al., 1989). As an alternative to administering probe drugs, individual CYP3A activity could also be determined using endogenous biomarkers, such as 4ß-OH-cholesterol (Tomalik-Scharte et al., 2009; Diczfalusy et al., 2011) and 6ß-OH-cortisol (Galteau and Shamsa, 2003) that are generated by CYP3A-mediated oxidation of cholesterol and cortisol.

Upon repeated drug administration of CYP3A substrates, a very high intra-individual constancy was reported for midazolam (Ozdemir et al., 2000). From such data, it was interpreted that up to 90% of variation in human CYP3A activity is caused by genetic factors (Ozdemir et al., 2000). However, individual constancy may also be due to constant environmental effects or epigenetic factors.

An analysis of genetic factors determining enzyme induction by St. John’s wort showed that ca. 66% of the inter-individual variation in induced CYP3A4 activity (measured with quinine as probe drug) may be due to additive genetic factors (Rahmioglu et al., 2011). However, induced CYP3A activity is not only reflecting the heritable factors influencing the pharmacokinetics of the CYP3A probe drug, but also by those that affect pharmacokinetics of the inducer (Rahmioglu et al., 2011).

Thus, at present, it is not clear how much of the inter-individual variation in CYP3A4 activity is in fact genetically determined. Twin studies performed in monozygotic (MZ) and dizygotic (DZ) twins allow to distinguish between additive (A) and dominant (D) genetic factors on the one hand, and constant (C) and unique environmental (E) factors on the other hand. We used this approach with midazolam as CYP3A4 probe drug and endogenous biomarkers of CYP3A activity. To allow comparison with the repeated drug administration approach (Ozdemir et al., 2000), each CYP3A test was performed three times in each individual. In conclusion, this study should contribute to a better understanding of the inter-individual variation of human CYP3A activity, which plays a relevant role in up to 40% of current drug therapies as well as in the metabolism of many endogenous and exogenous substances with biomedical relevance.

Materials and Methods

Study Design

The study was a repeated dose pharmacokinetic/pharmaco- genetic study on monozygotic (MZ) and dizygotic (DZ) twin pairs. The ethics committee of the University of Göttingen and the German Federal Drug Administration (BfArM) approved the study. It was registered at ClinicalTrials.gov (NCT01845194) and in the European clinical trials database (EUDRA-CT number: 2008-006223-31).

In total, 116 subjects including 44 MZ and 14 same-sex DZ twin pairs were included in the study. One MZ twin pair was excluded from the analyses of the intravenous midazolam administration because of erroneous blood sampling. All subjects were healthy according to their medical history, medical examination, electrocardiogram, urine status, and clinical chemistry analyses. The latter included sodium, potassium, calcium, aspartate aminotransferase, creatinine, total bilirubin, hemoglobin, erythrocyte, thrombocyte, and leucocyte counts. Written informed consent was obtained from each subject before participation in the study. Except for hormonal contraceptives, no other drugs or potentially interacting nutrients (grapefruits, grapefruit juice, St. John’s wort) were allowed for 1 week prior to the study days and until the last blood and urine sampling was completed. Smoking was allowed if both siblings of each pair smoked to a similar extent. Subjects with history of alcohol or drug addiction were excluded prior to any study related procedures. Additionally, subjects were not allowed to drink alcohol from 48 h prior to 72 h after midazolam application and alcohol breath test was absolved in the morning of each study day. Other dietary habits were recorded in standardized nutrition questionnaires (Winkler and Doring, 1998). Lean body weight was calculated from a formula based on body weight and height (Hume, 1966).

Midazolam was administered as a 2 min intravenous infusion with a total dose of 0.2 mg. The phenotyping test with midazolam was performed on three different occasions with a time interval of at least 1 week and maximally 3 months. Blood samples were drawn at 0 (before infusion), 0.25, 0.5, 0.75, 1, 1.5, 2, 2.5, 3, 4, 5, 6, 7, 8, and 24 h after the midazolam infusion.

Cholesterol and 4ß-OH cholesterol were measured in the baseline blood samples before drug administration. All samples were centrifuged for 10 min (4000 rpm) at 4°C and the plasma was stored at −20°C. To avoid possible circadian variation of the cortisol to 6ß-OH-cortisol metabolism (Ohno et al., 2000), a 24 h urine sampling for quantification of cortisol and 6ß-OH cortisol was performed on each of the three study days.

Bioanalytical Methods

The concentrations of midazolam, 1-OH-midazolam, 6-ß-hydroxycortisol and cortisol were determined by LC-MS-MS analysis on an Agilent 6460 triple quadrupole mass spectrometer (Agilent, Waldbronn, Germany). Sample work-up and LC-MS-MS analysis of midazolam was performed as described previously (Matthaei et al., 2015; Thiel et al., 2015) using midazolam-d6 (Roche Nr. RO0213981-010) as internal standard. For 1-OH-midazolam, 1 ml of plasma was extracted with diethyl ether as described (Becquemont et al., 2006) using 1-OH-midazolam-d5 (Toronto Research Chemicals, no. H948423) as internal standard. Mass spectrometric conditions and transitions in multiple reaction monitoring (MRM) mode were as described previously (Thiel et al., 2015). Standardization of the analytical assays was performed with calibration samples prepared in plasma in the concentration range from 0.04 to 20 nM for 1′-hydroxymidazolam, and 0.2–100 nM for midazolam. Precision of the analytical method at the limit of quantification, which was 0.04 nM for 1′-hydroxymidazolam and 0.2 nM for midazolam, could be characterized with a coefficient of variation (CV) of 6.6 and 3.7% and a bias of 6.1 or 1.5% for 1′-hydroxymidazolam and midazolam, respectively. At higher concentrations (0.1, 1, and 15 nM for 1′-hydroxymidazolam and 0.5, 5, and 75 nM for midazolam) the CV ranged from 1.3 to 8.1%, and the bias was between −3.8 and 4.2%.

Urine samples were centrifuged after thawing, and 10 μl were diluted with 90 μl of water for determination of cortisol and 6ß-OH-cortisol. After addition of internal standard mixture (10 pmol cortisol-d4 and 100 pmol 6ß-OH-cortisol-d4) the samples were centrifuged, and 10 μl of the supernatant were used for LC-MS-MS analysis. HPLC separation was achieved on a reversed phase column (SynergiTM 4 μm Polar RP 80Å, 150 × 2 mm, Phenomenex, Aschaffenburg, Germany) with a gradient of acetonitrile and water with 0.1% formic acid.

Cholesterol and 4ß-hydroxycholesterol were determined by GC-MS in EI mode, on a 5975 XL MSD coupled to a 7890A GC (Agilent, Waldbronn, Germany). Cholesterol quantification was performed as described previously (Lutjohann et al., 1995) with minor modifications. Briefly, 10 μl of plasma were spiked with 10 μg of cholesterol-d5 as internal standard. After saponification with 0.5 ml of 1 M NaOH in 90% ethanol at 70°C for 1 h, 250 μl of water were added and the samples extracted with 2 ml of n-hexane. An aliquot of 50 μl of the extract was evaporated to dryness and derivatized with 20 μl of N,O-bis(trimethylsilyl)trifluoroacetamide (BSTFA) for 30 min at room temperature.

Sample preparation for 4ß-hydroxycholesterol was performed according to a previously published method (Dzeletovic et al., 1995) with minor modifications. Briefly, 1 ml of plasma was spiked with 10 μg of BHT and 50 ng of the internal standard 4ß-hydroxycholesterol-d4 and saponified with 2.5 ml of 1 M NaOH in 90% ethanol at 70°C for 1 h under argon and then extracted with 1 ml of water and 5 ml of chloroform. The chloroform phase was evaporated to dryness in a stream of nitrogen and the residue dissolved in 1 ml of toluene. Samples were purified by solid phase extraction on silica cartridges (Isolute Si 100 mg, 3 ml, Biotage, Uppsala, Sweden) preconditioned with hexane. The cartridges were washed with 1 ml of hexane and 10 ml of 2-propanol in hexane (0.5% v/v) and then eluted with 2 ml of 2-propanol in hexane (30% v/v). The eluate was evaporated to dryness in a stream of nitrogen and derivatized with 20 μl of BSTFA for 30 min at room temperature. GC was performed on a J&W DB-5MS column (25 m, 0.2 mm i.d., 0.33 μm film thickness; Agilent) in the splitless mode. For the analysis of cholesterol, the GC oven program started at 150°C and was held for 1 min. Temperature was increased with 20°C/min to 300°C, with a total run time of 18.5 min. The trimethylsilyl derivatives of cholesterol and the internal standard cholesterol-d5 were detected in SIM mode at m/z 458 and 463, respectively. For the analysis of 4ß-hydroxycholesterol the GC oven program startet at 150°C and was held for 1 min. Temperature was increased with 10°C/min to 250°C, then with 30°C/min to 300°C. The trimethylsilyl derivatives of 4ß-hydroxycholesterol and the internal standard 4ß-hydroxycholesterol-d4 were quantified in SIM mode at m/z 366 and 370, respectively using m/z 456 and 460 as qualifier ions.

Calibration samples for cholesterol were prepared in isooctane with 10% 2-propanol in the concentration range from 1 to 30 μg. Calibration samples for 4ß-hydroxycholesterol were prepared in isooctane with 0.8% 2-propanol from 5 to 100 ng. Calibration samples for cortisol and 6ß-hydroxycortisol were prepared in water from 0.025 to 10 pmol, or 0.25 to 100 pmol, respectively. Calibration samples were worked up as the samples, and analyzed together with the unknown samples. Calibration curves based on internal standard calibration were obtained by weighted (1/x) linear regression for the peak area ratio of the analyte to the respective internal standard against the amount of the analyte. The concentration in unknown samples was obtained from the regression line.

Genotyping

DNA was isolated from whole venous blood by automated solid phase extraction with the EZ1TM DNA Blood 350 μl Kit using the Bio-robot EZ1TM (both Qiagen, Hilden, Germany). The CYP3A41b and 22 alleles as well as the CYP3A51/3 polymorphism were genotyped by DNA sequencing using fluorescence-labeled dideoxynucleotides and detection via Gene Mapper v3.7 Software® (Applied Biosystems®, Foster City, United States). Zygosity was assessed by analyzing variants in 23 polymorphic genes. Monozygosity was concluded when there were no differences regarding the polymorphisms in both siblings of one twin pair.

Pharmacokinetic Analysis

Pharmacokinetic parameters were estimated by non-compartmental analyses using the WinNonlin software (Pharsight Corporation, Mountain View, United States). The AUC from zero to infinity (AUC), was calculated by the linear/log trapezoidal rule and using the terminal elimination rate constant (lambda z). AUC420 was calculated from time of dose to 420 min after dosage. Total plasma clearance (Cl) after administration of midazolam was calculated as the ratio dose over AUC. Time (tmax) and amount of maximum plasma concentration (Cmax) were given as measured and the terminal half-life (t1/2) was calculated as t1/2 = ln(2)/lambda z.

Statistics and Analysis of Heritability

Predefined primary endpoint was midazolam AUC. All concentrations and pharmacokinetic parameters were determined on three independent occasions to compare intra- versus interindividual variation. The genetic component (rGC) concerning intra- (SDW) and inter-individual (SDB) standard deviation was calculated as described by Kalow et al. (1998, 1999). In structural equation modeling of heritability and in the other multifactorial analyses the individual means from the measurements or parameters obtained on the three occasions were used. Using structural equation analysis, variation was separated into additive genetic (A), common environmental (C), and unique environmental (E) effects. The heritability is the proportion of the variance resulting from the genetic effect (A). In order to investigate the genetic, common environmental, and unique environmental effects of the different traits taken into account, the possible correlation between traits were fitted in multivariate ACE, AE, and CE biometrical models using the mcglm package for multivariate regression analysis in R. As an exploratory analysis, we computed the correlation between MZ and DZ twins based on an unstructured (equal mean and variance for MZ and DZ twins) model. In all models, the mean structure for each trait was specified as a linear function of the covariates age and sex.

Results

We studied 44 monozygotic and 14 dizygotic twin pairs. There were no significant differences in demographic data between MZ and DZ twin pairs (Table 1). No serious adverse events related to the drugs or procedures of this study were observed.

TABLE 1
www.frontiersin.org

Table 1. Demographic data of study participants.

Midazolam was administered intravenously on three independent occasions under supine and resting conditions. As illustrated in Figure 1, there was a significant variation in midazolam blood concentrations between the individuals and an even bigger variation in 1-OH-midazolam blood concentrations. Inter-individually, there was a 3.1-fold variation in midazolam AUC (288 – 897.4 μgmin/L) and a 10.4-fold variation of 1-OH-midazolam AUC (30.95 – 322.8). There were no statistically significant differences in the pharmacokinetic parameters between the monozygotic and dizygotic twins (Table 2). Intra-individually (i.e., between-days), the relative variation was similar for midazolam and for 1-OH-midazolam with mean coefficients of variation of 11.3 and 12.3%.

FIGURE 1
www.frontiersin.org

Figure 1. Concentration time curves of midazolam (A) and its major metabolite 1-OH-midazolam (B) after intravenous administration of 0.2 mg midazolam. Data for all subjects of one study day are shown. Concentrations of midazolam varied about 10-fold and that of the CYP3A dependent metabolite even higher. Midazolam AUC (C) and 1-OH-midazolam AUC (D) is shown below illustrating the intra-individual variation and the variation between the three study days. The genetic component (Kalow et al., 1998) was determined from this data independently of the twin status, with values of 0.57 (95% CI = 0.54 – 0.60) for the AUC of midazolam and 0.93 (95% CI = 0.93 – 0.94) for the AUC of 1-OH-midazolam.

TABLE 2
www.frontiersin.org

Table 2. Pharmacokinetic parameters of midazolam and 1-OH-midazolam.

As illustrated in the lower part of Figure 1, in some, but not in all volunteers a high intraindividual constancy of midazolam and 1-OH-midazolam AUC was found. From comparing the intra- and interindividual variation, the so-called genetic component was calculated. The genetic component is an alternative approach suggested as an indicator of heritability (Kalow et al., 1998, 1999). This parameter is determined completely independently of the twin status and solely based on intra- and inter-individual variability. The genetic component was 0.612 for midazolam and 0.931 for 1-OH-midazolam (Table 3). That means, that up to 61% of the variation in midazolam and up to 93% of the variation in 1-OH-midazolam may be due to genetic factors. This parameter was also calculated for the endogenous biomarkers plasma 4ß-OH-cholesterol and 24-h urinary 6ß-OH-cortisol with a quite high genetic component of 0.812 for the cholesterol metabolite. However, in interpretation of this high value it has to be taken into account that 4ß-OH-cholesterol has a long half-life of more than 60 h (Bodin et al., 2002), and thus some of the individual constancy is simply due to that. As illustrated in Table 4, there were no significant differences in the endogenous CYP3A biomarkers between monozygotic and dizygotic twins.

TABLE 3
www.frontiersin.org

Table 3. Within- and between subject variance and genetic component of the 4 biomarkers.

TABLE 4
www.frontiersin.org

Table 4. Endogenous biomarkers of CYP3A activity.

For all four indicators of CYP3A activity, i.e., midazolam, 1-OH-midazolam, 4ß-OH-cholesterol, and 6ß-OH-cortisol, there was a significant correlation between both, the MZ and the DZ pairs (Table 5 and Figure 2). A relatively high correlation of all 4 biomarkers even in the dizygotic twins might indicate that the respective parameters are modulated not only be heritable factors, but also by common environmental factors. But with all four biomarkers, the correlations between MZ twins were closer than between DZ twins, which supports the hypothesis that the biomarkers are also modulated by inherited effects (Table 5). Using structural equation analysis, based on the ACE model including additive genetic (A), common environment (C), and unique environmental (E) effects, the study showed a significant and high heritability only for 1-OH-midazolam (p = 0.03) and 4ß-hydroxy-cholesterol (p < 0.01) (Table 5, ACE columns). According to this data, 73% of the variation in 1-OH-midazolam and 56% of variation in plasma 4ß-OH-cholesterol were due to inherited factors. In contrast, variations of midazolam AUC and of 24-h urine 6ß-OH-cortisol were less strongly determined by genetic factors but more variation was explained by common environmental factors. For comparison, the two more parsimonious models are also given in Table 5, attributing all variation either to heritability (AE model) or to common environmental effects (CE model).

TABLE 5
www.frontiersin.org

Table 5. Heritable and environmental effects on four known or presumed probe substrates of CYP3A activity.

FIGURE 2
www.frontiersin.org

Figure 2. Correlation of the biomarkers for CYP3A activity between monozygotic and dizygotic twins. The upper lane shows midazolam AUC. Below is the AUC of 1-OH-midazolam, 4ß-OH-cholesterol in plasma, and 6ß-OH-cortisol in urine over 24 h. Data represents the mean of three measurements performed on separate days.

Correlation Between the Biomarkers of CY3A Activity

As illustrated in Figure 3, correlations between the four biomarkers were surprisingly low or in some comparisons even absent. Notably, there was no correlation between 1-OH-midazolam and the two endogenous biomarkers of CYP3A. The simple phenotypic correlation as illustrated in Figure 3, does not indicate the origin of the correlation. However, in multivariate structural equation modeling (Table 6), one may employ either genetic and environmental correlations or phenotypically standardized covariances to assess the structure of the genetic and environmental influences. The genetic correlation is an estimate of the additive genetic effect that is shared between a pair of traits. It is obtained as a standardization of the additive genetic covariance between traits by their respective additive variances. As shown in Table 6, that parameter was statistically significant with coefficients of 0.4 only for midazolam and its metabolite (which may appear self-evident because of the drug-metabolite relationship), and, interestingly, for midazolam and 6ß-OH-cortisol, reassuring concerning the use of the latter two biomarkers as indicators of partially heritable CYP3A activity.

FIGURE 3
www.frontiersin.org

Figure 3. Correlations between the four indicators of CYP3A activity: plasma AUC of midazolam and 1-OH-midazolam, plasma concentration of 4ß-OH-cholesterol and the amount of 6ß-OH-cortisol excreted in 24-h urine. Linear regression coefficients are shown in the figure (*p < 0.05, **p < 0.01). As seen, there was a relatively high correlation between midazolam and 6ß-OH-cortisol, whereas the CYP3A generated metabolite 1-OH-midazolam did not correlate at all with the alternative endogenous in vivo biomarkers of CYP3A activity. Only frequent genetic polymorphisms with unequivocal functional relevance were considered, the CYP3A4*22 variant conferring low CYP3A4 activity and the CYP3A5*3 variant conferring loss of CYP3A5 activity.

TABLE 6
www.frontiersin.org

Table 6. Phenotypic correlation, genetic correlation, non-shared correlation, and bivariate heritability of 4 biomarkers of CYP3A activity.

Variance Explained by Environmental Factors and Molecularly Defined Polymorphisms

Midazolam pharmacokinetics were significantly influenced by sex, lean body weight, type of nutrition, usage of hormonal contraceptives, and by the CYP3A422 polymorphism (Table 7). Median (minimum – maximum) midazolam AUC was 8% higher in male than in female subjects [466.0 μg min/L (290.6 – 668.1 μg min/L) versus 432.1 μg min/L (287.8 – 897.4 μg min/L)]. Higher lean body mass was associated with low AUC according to linear regression analysis. The median midazolam AUC in heterozygous carriers of the CYP3A422 polymorphism (n = 4) was 35% higher than in subjects with the 1/1 and 1/1B genotypes (n = 110) (595.1 μg min/L versus 439.1 μg min/L).

TABLE 7
www.frontiersin.org

Table 7. Genetic and environmental effects on the CYP3A probe drugs and biomarkers.

The AUC of 1-OH-midazolam was significantly lower with increasing age (p = 0.012), but there were no statistically significant effects of the CYP3A genotypes (Table 7). Interestingly, plasma 4ß-OH-cholesterol was significantly lower in carriers of the CYP3A53 allele (p = 0.02), with a median 4ß-OH-cholesterol level of 23.01 ng/ml in homozygous carriers of the CYP3A53 allele and 30.69 ng/ml in heterozygous carriers of the active CYP3A51 allele. Our study included 16 carriers of the CYP3A51 genotype, and from the genotype-related differences (Table 3) one may conclude that in vivo CYP3A5 contributes in a relevant manner to cholesterol hydroxylation whereas CYP3A5 does not appear to contribute in a relevant manner to midazolam or cortisol hydroxylation. However, a minor contribution cannot be ruled out since our study was not designed to quantify the effects of specific CYP3A4 or CYP3A5 genetic variants.

Nutrition may affect CYP3A activity and expression. However, the analysis of nutritional effects was not a primary aim of the present study, and the data should thus be considered as explorative. High protein consumption was associated in a dose-dependent manner with low 1-OH-midazolam AUC and low 6ß-OH-cortisol. However, no effects of diet on 4ß-OH-cholesterol were found. High fruit intake was apparently associated with a high midazolam clearance, with a median AUC of 425.4, 446.9, 519.4, and 558.9 μg min/L in subjects with daily, thrice weekly, rare and very rare fruit consumption, respectively. Alcohol consumption may also modulate CYP3A activity but our data could not confirm that because, according to the standardized nutrition questionnaires, only three participants stated that they drank alcohol 1–3 times a week, and all the others drank less alcohol according to the questionnaires. According to alcohol breath tests and according to the liver function tests, there was also no indication of alcohol consumption.

Discussion

Human CYP3A activity is highly variable between subjects. Some data published earlier indicated that up to 98% of that variation may be heritable (Ozdemir et al., 2000; Klein and Zanger, 2013). Up to 66% of induced CYP3A activity may be heritable, according to a twin study using a herbal product, St. John’s wort, for enzyme induction (Rahmioglu et al., 2011). However, there is scarce data from twin studies on heritability of not induced CYP3A activity. Because much of the current thinking about a high heritability in CYP3A activity is based on the high intra-individual constancy of CYP3A in vivo indicators, we designed our study to allow within the study a comparison between that repeated dose approach or genetic component approach (Ozdemir et al., 2000) and the twins study approach specifically reflecting heritabilty. As show in Table 3, there was significant intraindividual constancy but also some variation particularly in midazolam AUC (Figure 1C). Thus, between the study days environmental effects currently not classified as strong CYP3A inducers or inhibitors apparently may have had modified enzyme expression and activity and further research on such possibly factors (e.g., specific foods or changes in the intestinal microbiome) that may reveal interesting new insights. Our study participants were instructed to avoid any known inducer or inhibitor of CYP3A enzymes. Nevertheless, not yet identified inducers or inhibitors from nutrition, from the microbiome and from endogenous metabolism may have modulated CYP3A expression and activity in our study.

With a relatively small number of dizygotic twins, the differentiation between additive genetic effects (A) and constantly acting environmental (including epigenetic) effects (C) was not very precise as given in Table 5. But the sum of both, A and C (Table 5), was about 60% and that surprisingly well corresponded to the independent analytical approach based in intra- versus interindividual variation (Table 3).

A high intra-individual constancy (Figure 1D and Table 3) and a high heritability was indeed found for 1-OH-midazolam AUC (Table 5), but apparently this was due to a high heritability of glucuronyltransferases or other processes relevant for elimination of 1-OH-midazolam. In general, the low to moderate correlation between all four markers indicates that each biomarker also reflects processes beside CYP3A activity. The twin study design with multivariate structural equation analyses allowed exploring more about the causes. A low correlation of 4ß-OH-cholesterol with the other parameters may be due to the long half-life of 4ß-OH-cholesterol. Therefore, this biomarker does more reflect the overall integrated effects of CYP3A inducing and inhibiting factors over a longer time span, whereas the other three biomarkers better reflect the metabolic status at the day of testing. Indeed, the highest genetic correlation of −0.42 was found between midazolam and 6ß-OH-cortisol (Table 6) consistent with the assumption that both reflect a common genetic component in the CYP3A genes or in other genes regulating CYP3A genes.

Due to its high sensitivity to inhibitors and inducers, midazolam is currently a preferred CYP3A probe drug in drug research, but data shown here indicate that midazolam may also reflect other processes and activities besides CYP3A. With a clearance of about 0.5 L/min (Table 2), midazolam is a so-called high clearance drug (Foster et al., 2011; Burke et al., 2014) and all high-clearance drugs do depend in their elimination strongly on liver blood flow (Nies et al., 1976). However, the factors known to affect liver blood flow like physical activity or food intake, were highly controlled here on all three study days. Nevertheless, there was an apparently significant intra-individual variation in several volunteers between the three study days (Figure 1). With a mean coefficient of variation of 11.3% for intra-individual variation, the intra-individual variation in our study was overall not significantly higher than that observed in an earlier study also with repeated intravenous administration of midazolam (Kashuba et al., 1998). A moderately bigger intra-individual variation in our study may simply be due to the fact that in our study the time between the measurements was longer (on average 1 month). Ideally, time between measurements should not be relevant for truly only genetically determined factors, but the genetic component (Kalow et al., 1998) does not differentiate between A and C (additive genetic effects and constant environmental effects, Table 5). This intraindividual variation together with the common genetic correlation of only 0.42 between 6ß-OH-cortisol and midazolam indicates that there were significant random and individually constant effects modulating the midazolam pharmacokinetics in humans. Such constant environmental factors may include long-term epigenetic imprinting but also constant factors in lifestyle and nutrition. Subjects in our study were not allowed to consume known inducers or inhibitors of CYP3A like grapefruits or St. John’s wort. However, in an explorative analysis based on standardized nutrition questionnaires high fruit intake was associated with high midazolam clearance. This may support the hypothesis that there are CYP3A inducers in some fruits and vegetables, but the food questionnaire was not specific enough and sample size was not sufficient to analyze that in more detail. But we think that this is an interesting point of further research, possibly also including the interactions between CYP3A activity, food and the intestinal microbiome.

In vitro, one often prefers measuring enzyme activity based on the metabolite formation. In vivo, in our study, there was a high intra-individual constancy and significant heritability of the main midazolam metabolite 1-OH-midazolam (Figure 1 and Tables 3, 5). However, the 1-OH-midazolam AUC is determined by both metabolite formation and elimination. While metabolite formation is mostly mediated by CYP3A4, the AUC of 1-OH-midazolam was not different for the functionally relevant CYP3A422 polymorphisms and correlated only poorly with the other CYP3A biomarkers. This indicates that the variation in elimination of 1-OH-midazolam may be more important than the variation in its formation for our understanding of the 1-OH-midazolam plasma concentrations. That is also strongly supported by another finding: As seen in Figure 3, there was a positive correlation between midazolam and 1-OH-midazolam. With the substrate-product relationship between midazolam and 1-OH-midazolam that correlation should have been negative. This paradox is best explained by the assumption that plasma 1-OH-midazolam concentrations were much more dependent on factors other than CYP3A activity. That corresponds to earlier findings that under rifampin induction there was less plasma 1-OH-midazolam contrary to the expectations and only with analysis of the 1-OH-midazolam glucuronide there was more after rifampin induction (Link et al., 2008). Unfortunately, in our study the 1-OH-midazolam glucuronide was not analyzed. Elimination of 1-OH-midazolam may mostly be determined by glucuronidation via uridine 5’-diphospho-glucuronosyltransferase (UGT) enzymes 1A4, 2B6, and 2B7 (Seo et al., 2010; Stingl et al., 2014). The apparent high heritability in 1-OH-midazolam AUC may reflect, to a great extent, a high heritability of the corresponding UGT enzyme. Alternatively, membrane transporters relevant for 1-OH-midazolam and its glucuronide may be highly heritable. This remains speculative without further experimental identification, but at least it may teach us two points: firstly, drug-to-metabolite ratios may sometimes be misleading (or at least not increase specificity), and, secondly, highly predictive genetic polymorphisms may exist in the relevant UGT enzymes or in the glucuronide membrane transport proteins which should be further studied.

Our data on the poor correlation between 4ß-OH-cholesterol and midazolam AUC is in line with results published by others: For instance, the correlation between midazolam total clearance and 4ß-OH-cholesterol plasma concentration was only 0.24 (p = 0.048) in 50 healthy volunteers (Tomalik-Scharte et al., 2009). This is similar to the results obtained here (Figure 3). One explanation lies in the large differences in pharmacokinetics. While midazolam has an elimination half-life of about 2 h (Table 2), the half-life of 4ß-OH-cholesterol was estimated to be 60 h by application of radiolabeled substance to humans (Bodin et al., 2002) but other studies estimated that half-life even as 17 days or longer (Diczfalusy et al., 2008). Thus, 4ß-OH-cholesterol is an indicator for the mean average CYP3A activity over weeks whereas midazolam reflects CYP3A activity within a time span of a few hours only. There may exist moderate circadian variation in CYP3A activity, but this can only explain part of the missing correlation. Again, the pathways mediating further elimination of 4ß-OH-cholesterol by membrane transport and metabolism may play a role, and thus, the high heritability seen with 4ß-OH-cholesterol may not only be heritability related to CYP3A activity.

The cortisol metabolite 6ß-OH-cortisol has been in use in clinical pharmacology for more than 50 years (Roots et al., 1975). Under baseline conditions excluding any of the known strong inducers, there was only a moderate correlation between midazolam and 6ß-OH-cortisol in our study (Figure 3). Other studies found even no correlation at all between midazolam systemic clearance and the 6ß-OH-cortisol urinary ratio under such baseline conditions (Eeckhoudt et al., 2001; Chan et al., 2016). However, a high correlation was found between midazolam clearance and the 6ß-OH-cortisol ratio following hepatic CYP3A inhibition and induction (Shin et al., 2016). The endogenous plasma concentrations of 4ß-OH-cholesterol were similarly identified as a suitable biomarker correlating with midazolam after long-term CYP3A induction (Kasichayanula et al., 2014). However, with enzyme induction, probably all four in vivo markers of CYP3A activity used here most likely would have correlated much stronger than found here. Heritability of induced CYP3A activity was indeed found to be about 66% after a 14 days treatment with St. John’s wort (Rahmioglu et al., 2011), but the heritability of baseline (not significantly induced) CYP3A activity was not studied there.

Conclusion

In a population not exposed to known inducers or inhibitors, heritability of CYP3A activity was found to be relatively small, particularly when using midazolam as biomarker. Considering the known very strong susceptibility of CYP3A4 and CYP3A5 to enzyme induction and enzyme inhibition, this finding may not appear unexpected. In our study, the estimates of CYP3A heritability varied depending on the biomarker used, which could partially be explained by the specific features of each probe substance, but it clearly indicates that the one and only correct in vivo biomarker of CYP3A activity does not exist. The multivariate twin study design allowed searching for the common hidden factor using structural equation modeling, but this did not identify a strong genetic correlation between all four presumed CYP3A biomarkers. The relatively poor correlation between different biomarkers of CYP3A activity may not generally devalue them, but we should keep in mind that CYP3A biomarkers may reflect to a significant extent more of other processes than commonly presumed. The fraction of variation in CYP3A4 activity currently explained by defined genetic polymorphisms is relatively small and according to our study, search for further genomic biomarkers to predict reliably the individual CYP3A activity may not be very promising. Instead, functional CYP3A biomarkers may be more suitable for therapy individualization and further search for epigenetic factors modulating intermediate- and long-term CYP3A activity may be promising. However, the high heritability of 1-OH-midazolam might stimulate further research on the underlying genetic polymorphisms relevant for this metabolite.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics Statement

The studies involving human participants were reviewed and approved by the Ethikkommission der Universitätsmedizin Göttingen, Georg-August-Universität Göttingen. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

JM, WB, and JB wrote the manuscript. JB, DS, MS, and RK designed the research. JM, JB, MT, JS, DS, CS-S, and UH performed the research. WB, JB, JM, and JB analyzed the data. All authors contributed to the article and approved the submitted version.

Funding

This study was funded by the Georg-August University Göttingen. Additional external funding was received from AstraZeneca and from the German Ministry of Education and Research, BMBF, project “Virtual Liver” grant 0315755 and from the Robert-Bosch Foundation (Stuttgart, Germany).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We gratefully acknowledge the contributions of Josephine Salaw, Valerie Gal, and Ulrike Schöning to the clinical study conduct and Karoline Jobst for the laboratory analyses.

References

Becquemont, L., Glaeser, H., Drescher, S., Hitzl, M., Simon, N., Murdter, T. E., et al. (2006). Effects of ursodeoxycholic acid on P-glycoprotein and cytochrome P450 3A4-dependent pharmacokinetics in humans. Clin. Pharmacol. Therap. 79, 449–460. doi: 10.1016/j.clpt.2006.01.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Bodin, K., Andersson, U., Rystedt, E., Ellis, E., Norlin, M., Pikuleva, I., et al. (2002). Metabolism of 4 beta -hydroxycholesterol in humans. J. Biol. Chem. 277, 31534–31540.

Google Scholar

Bodin, K., Bretillon, L., Aden, Y., Bertilsson, L., Broome, U., Einarsson, C., et al. (2001). Antiepileptic drugs increase plasma levels of 4beta-hydroxycholesterol in humans: evidence for involvement of cytochrome p450 3A4. J. Biol. Chem. 276, 38685–38689. doi: 10.1074/jbc.m105127200

PubMed Abstract | CrossRef Full Text | Google Scholar

Bodin, K., Lindbom, U., and Diczfalusy, U. (2005). Novel pathways of bile acid metabolism involving CYP3A4. Biochim. Biophys. Acta 1687, 84–93. doi: 10.1016/j.bbalip.2004.11.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Bui, K., Zhou, D., Sostek, M., She, F., and Al-Huniti, N. (2016). Effects of CYP3A modulators on the pharmacokinetics of Naloxegol. J. Clin. Pharmacol. 56, 1019–1027. doi: 10.1002/jcph.693

PubMed Abstract | CrossRef Full Text | Google Scholar

Burke, R. S., Somasuntharam, I., Rearden, P., Brown, D., Deshmukh, S. V., Dipietro, M. A., et al. (2014). siRNA-mediated knockdown of P450 oxidoreductase in rats: a tool to reduce metabolism by CYPs and increase exposure of high clearance compounds. Pharm Res 31, 3445–3460. doi: 10.1007/s11095-014-1433-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Chan, S. W., Xiao, Y., Hu, M., Yin, O. Q., Chu, T. T., Fok, B. S., et al. (2016). Associations of the CYP3A53 and CYP3A41G polymorphisms with the pharmacokinetics of oral midazolam and the urinary 6beta-hydroxycortisol/cortisol ratio as markers of CYP3A activity in healthy male Chinese. J. Clin. Pharm. Ther. 41, 552–558. doi: 10.1111/jcpt.12433

PubMed Abstract | CrossRef Full Text | Google Scholar

Diczfalusy, U., Miura, J., Roh, H. K., Mirghani, R. A., Sayi, J., Larsson, H., et al. (2008). 4Beta-hydroxycholesterol is a new endogenous CYP3A marker: relationship to CYP3A5 genotype, quinine 3-hydroxylation and sex in Koreans, Swedes and Tanzanians. Pharmacogenet. Genom. 18, 201–208. doi: 10.1097/fpc.0b013e3282f50ee9

PubMed Abstract | CrossRef Full Text | Google Scholar

Diczfalusy, U., Nylen, H., Elander, P., and Bertilsson, L. (2011). 4beta-Hydroxycholesterol, an endogenous marker of CYP3A4/5 activity in humans. Br. J. Clin. Pharmacol. 71, 183–189. doi: 10.1111/j.1365-2125.2010.03773.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Dzeletovic, S., Breuer, O., Lund, E., and Diczfalusy, U. (1995). Determination of cholesterol oxidation products in human plasma by isotope dilution-mass spectrometry. Analyt. Biochem. 225, 73–80. doi: 10.1006/abio.1995.1110

PubMed Abstract | CrossRef Full Text | Google Scholar

Eap, C. B., Buclin, T., Cucchia, G., Zullino, D., Hustert, E., Bleiber, G., et al. (2004a). Oral administration of a low dose of midazolam (75 microg) as an in vivo probe for CYP3A activity. Eur. J. Clin. Pharmacol. 60, 237–246.

Google Scholar

Eap, C. B., Buclin, T., Hustert, E., Bleiber, G., Golay, K. P., Aubert, A. C., et al. (2004b). Pharmacokinetics of midazolam in CYP3A4- and CYP3A5-genotyped subjects. Eur. J. Clin. Pharmacol. 60, 231–236.

Google Scholar

Eeckhoudt, S. L., Desager, J. P., Robert, A. R., Leclercq, I., Verbeeck, R. K., and Horsmans, Y. (2001). Midazolam and cortisol metabolism before and after CYP3A induction in humans. Int. J. Clin. Pharmacol. Ther. 39, 293–299. doi: 10.5414/cpp39293

PubMed Abstract | CrossRef Full Text | Google Scholar

Evans, W. E., and Relling, M. V. (1999). Pharmacogenomics: translating functional genomics into rational therapeutics. Science 286, 487–491. doi: 10.1126/science.286.5439.487

PubMed Abstract | CrossRef Full Text | Google Scholar

Foster, J. A., Houston, J. B., and Hallifax, D. (2011). Comparison of intrinsic clearances in human liver microsomes and suspended hepatocytes from the same donor livers: clearance-dependent relationship and implications for prediction of in vivo clearance. Xenobiotica 41, 124–136. doi: 10.3109/00498254.2010.530700

PubMed Abstract | CrossRef Full Text | Google Scholar

Fromm, M. F., Schwilden, H., Bachmakov, I., Konig, J., Bremer, F., and Schuttler, J. (2007). Impact of the CYP3A5 genotype on midazolam pharmacokinetics and pharmacodynamics during intensive care sedation. Eur. J. Clin. Pharmacol. 63, 1129–1133. doi: 10.1007/s00228-007-0365-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Fuhr, U., Jetter, A., and Kirchheiner, J. (2007). Appropriate phenotyping procedures for drug metabolizing enzymes and transporters in humans and their simultaneous use in the “cocktail” approach. Clin. Pharmacol. Ther. 81, 270–283. doi: 10.1038/sj.clpt.6100050

PubMed Abstract | CrossRef Full Text | Google Scholar

Galteau, M. M., and Shamsa, F. (2003). Urinary 6beta-hydroxycortisol: a validated test for evaluating drug induction or drug inhibition mediated through CYP3A in humans and in animals. Eur. J. Clin. Pharmacol. 59, 713–733. doi: 10.1007/s00228-003-0690-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Gashaw, I., Kirchheiner, J., Goldammer, M., Bauer, S., Seidemann, J., Zoller, K., et al. (2003). Cytochrome p450 3A4 messenger ribonucleic acid induction by rifampin in human peripheral blood mononuclear cells: correlation with alprazolam pharmacokinetics. Clin. Pharmacol. Ther. 74, 448–457. doi: 10.1016/s0009-9236(03)00237-6

CrossRef Full Text | Google Scholar

Gorski, J. C., Jones, D. R., Haehner-Daniels, B. D., Hamman, M. A., and O’mara, E. M. Jr. (1998). The contribution of intestinal and hepatic CYP3A to the interaction between midazolam and clarithromycin. Clin. Pharmacol. Ther. 64, 133–143. doi: 10.1016/s0009-9236(98)90146-1

CrossRef Full Text | Google Scholar

Greenblatt, D. J., and Harmatz, J. S. (2015). Ritonavir is the best alternative to ketoconazole as an index inhibitor of cytochrome P450-3A in drug-drug interaction studies. Br. J. Clin. Pharmacol. 80, 342–350. doi: 10.1111/bcp.12668

PubMed Abstract | CrossRef Full Text | Google Scholar

Heizmann, P., Eckert, M., and Ziegler, W. H. (1983). Pharmacokinetics and bioavailability of midazolam in man. Br. J. Clin. Pharmacol. 16(Suppl. 1), 43S–49S. doi: 10.1111/j.1365-2125.1983.tb02270.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Hume, R. (1966). Prediction of lean body mass from height and weight. J. Clin. Pathol. 19, 389–391. doi: 10.1136/jcp.19.4.389

PubMed Abstract | CrossRef Full Text | Google Scholar

Kalow, W., Endrenyi, L., and Tang, B. (1999). Repeat administration of drugs as a means to assess the genetic component in pharmacological variability. Pharmacology 58, 281–284. doi: 10.1159/000028292

PubMed Abstract | CrossRef Full Text | Google Scholar

Kalow, W., Tang, B. K., and Endrenyi, L. (1998). Hypothesis: comparisons of inter- and intra-individual variations can substitute for twin studies in drug research. Pharmacogenetics 8, 283–289. doi: 10.1097/00008571-199808000-00001

PubMed Abstract | CrossRef Full Text | Google Scholar

Kashuba, A. D., Bertino, J. S. Jr., Rocci, M. L. Jr., and Kulawy, R. W. (1998). Quantification of 3-month intraindividual variability and the influence of sex and menstrual cycle phase on CYP3A activity as measured by phenotyping with intravenous midazolam. Clin. Pharmacol. Ther. 64, 269–277. doi: 10.1016/s0009-9236(98)90175-8

CrossRef Full Text | Google Scholar

Kasichayanula, S., Boulton, D. W., Luo, W. L., Rodrigues, A. D., Yang, Z., Goodenough, A., et al. (2014). Validation of 4beta-hydroxycholesterol and evaluation of other endogenous biomarkers for the assessment of CYP3A activity in healthy subjects. Br. J. Clin. Pharmacol. 78, 1122–1134. doi: 10.1111/bcp.12425

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, H., Yoon, Y. J., Shon, J. H., Cha, I. J., Shin, J. G., and Liu, K. H. (2006). Inhibitory effects of fruit juices on CYP3A activity. Drug Metab. Dispos. 34, 521–523. doi: 10.1124/dmd.105.007930

PubMed Abstract | CrossRef Full Text | Google Scholar

Kirby, B., Kharasch, E. D., Thummel, K. T., Narang, V. S., Hoffer, C. J., and Unadkat, J. D. (2006). Simultaneous measurement of in vivo P-glycoprotein and cytochrome P450 3A activities. J. Clin. Pharmacol. 46, 1313–1319. doi: 10.1177/0091270006292625

PubMed Abstract | CrossRef Full Text | Google Scholar

Klein, K., and Zanger, U. M. (2013). Pharmacogenomics of cytochrome P450 3A4: recent progress toward the “missing heritability”. Problem. Front. Genet. 4:12. doi: 10.3389/fgene.2013.00012

PubMed Abstract | CrossRef Full Text | Google Scholar

Koch, I., Weil, R., Wolbold, R., Brockmoller, J., Hustert, E., Burk, O., et al. (2002). Interindividual variability and tissue-specificity in the expression of cytochrome P450 3A mRNA. Drug Metab. Dispos. 30, 1108–1114. doi: 10.1124/dmd.30.10.1108

PubMed Abstract | CrossRef Full Text | Google Scholar

Kronbach, T., Mathys, D., Umeno, M., Gonzalez, F. J., and Meyer, U. A. (1989). Oxidation of midazolam and triazolam by human liver cytochrome P450IIIA4. Mol. Pharmacol. 36, 89–96.

Google Scholar

Link, B., Haschke, M., Grignaschi, N., Bodmer, M., Aschmann, Y. Z., Wenk, M., et al. (2008). Pharmacokinetics of intravenous and oral midazolam in plasma and saliva in humans: usefulness of saliva as matrix for CYP3A phenotyping. Br. J. Clin. Pharmacol. 66, 473–484. doi: 10.1111/j.1365-2125.2008.03201.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Lutjohann, D., Bjorkhem, I., Beil, U. F., and Von Bergmann, K. (1995). Sterol absorption and sterol balance in phytosterolemia evaluated by deuterium-labeled sterols: effect of sitostanol treatment. J. Lipid Res. 36, 1763–1773.

Google Scholar

Matthaei, J., Brockmoller, J., Tzvetkov, M. V., Sehrt, D., Sachse-Seeboth, C., Hjelmborg, J. B., et al. (2015). Heritability of metoprolol and torsemide pharmacokinetics. Clin. Pharmacol. Ther. 98, 611–621. doi: 10.1002/cpt.258

PubMed Abstract | CrossRef Full Text | Google Scholar

Nies, A. S., Shand, D. G., and Wilkinson, G. R. (1976). Altered hepatic blood flow and drug disposition. Clin. Pharmacokinet. 1, 135–155. doi: 10.2165/00003088-197601020-00005

PubMed Abstract | CrossRef Full Text | Google Scholar

Ohno, M., Yamaguchi, I., Ito, T., Saiki, K., Yamamoto, I., and Azuma, J. (2000). Circadian variation of the urinary 6beta-hydroxycortisol to cortisol ratio that would reflect hepatic CYP3A activity. Eur. J. Clin. Pharmacol. 55, 861–865. doi: 10.1007/s002280050708

PubMed Abstract | CrossRef Full Text | Google Scholar

Ozdemir, V., Kalow, W., Tang, B. K., Paterson, A. D., Walker, S. E., Endrenyi, L., et al. (2000). Evaluation of the genetic component of variability in CYP3A4 activity: a repeated drug administration method. Pharmacogenetics 10, 373–388.

Google Scholar

Patki, K. C., Von Moltke, L. L., and Greenblatt, D. J. (2003). In vitro metabolism of midazolam, triazolam, nifedipine, and testosterone by human liver microsomes and recombinant cytochromes p450: role of cyp3a4 and cyp3a5. Drug Metab. Dispos. 31, 938–944. doi: 10.1124/dmd.31.7.938

PubMed Abstract | CrossRef Full Text | Google Scholar

Rahmioglu, N., Heaton, J., Clement, G., Gill, R., Surdulescu, G., Zlobecka, K., et al. (2011). Genetic epidemiology of induced CYP3A4 activity. Pharmacogenet. Genom. 21, 642–651. doi: 10.1097/fpc.0b013e3283498ecf

PubMed Abstract | CrossRef Full Text | Google Scholar

Roots, I., Saalfrank, K., and Hildebrandt, A. G. (1975). Comparison of methods to study enzyme induction in man. Adv. Exp. Med. Biol. 58, 485–502. doi: 10.1007/978-1-4615-9026-2_33

CrossRef Full Text | Google Scholar

Seo, K. A., Bae, S. K., Choi, Y. K., Choi, C. S., Liu, K. H., and Shin, J. G. (2010). Metabolism of 1’- and 4-hydroxymidazolam by glucuronide conjugation is largely mediated by UDP-glucuronosyltransferases 1A4, 2B4, and 2B7. Drug Metab. Dispos. 38, 2007–2013. doi: 10.1124/dmd.110.035295

PubMed Abstract | CrossRef Full Text | Google Scholar

Shin, K. H., Ahn, L. Y., Choi, M. H., Moon, J. Y., Lee, J., Jang, I. J., et al. (2016). Urinary 6beta-Hydroxycortisol/cortisol ratio most highly correlates with midazolam clearance under hepatic CYP3A inhibition and induction in females: a pharmacometabolomics approach. AAPS J. 18, 1254–1261. doi: 10.1208/s12248-016-9941-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Stingl, J. C., Bartels, H., Viviani, R., Lehmann, M. L., and Brockmoller, J. (2014). Relevance of UDP-glucuronosyltransferase polymorphisms for drug dosing: a quantitative systematic review. Pharmacol. Ther. 141, 92–116. doi: 10.1016/j.pharmthera.2013.09.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Tateishi, T., Watanabe, M., Nakura, H., Asoh, M., Shirai, H., Mizorogi, Y., et al. (2001). CYP3A activity in European American and Japanese men using midazolam as an in vivo probe. Clin. Pharmacol. Ther. 69, 333–339. doi: 10.1067/mcp.2001.115447

PubMed Abstract | CrossRef Full Text | Google Scholar

Thiel, C., Schneckener, S., Krauss, M., Ghallab, A., Hofmann, U., Kanacher, T., et al. (2015). A systematic evaluation of the use of physiologically based pharmacokinetic modeling for cross-species extrapolation. J. Pharm. Sci. 104, 191–206. doi: 10.1002/jps.24214

PubMed Abstract | CrossRef Full Text | Google Scholar

Thummel, K. E., Shen, D. D., Podoll, T. D., Kunze, K. L., Trager, W. F., Hartwell, P. S., et al. (1994). Use of midazolam as a human cytochrome P450 3A probe: I. In vitro-in vivo correlations in liver transplant patients. J. Pharmacol. Exp. Ther. 271, 549–556.

Google Scholar

Tomalik-Scharte, D., Lutjohann, D., Doroshyenko, O., Frank, D., Jetter, A., and Fuhr, U. (2009). Plasma 4beta-hydroxycholesterol: an endogenous CYP3A metric? Clin. Pharmacol. Ther. 86, 147–153. doi: 10.1038/clpt.2009.72

PubMed Abstract | CrossRef Full Text | Google Scholar

Tomalik-Scharte, D., Suleiman, A. A., Frechen, S., Kraus, D., Kerkweg, U., Rokitta, D., et al. (2014). Population pharmacokinetic analysis of circadian rhythms in hepatic CYP3A activity using midazolam. J. Clin. Pharmacol. 54, 1162–1169. doi: 10.1002/jcph.318

PubMed Abstract | CrossRef Full Text | Google Scholar

Winkler, G., and Doring, A. (1998). Validation of a short qualitative food frequency list used in several German large scale surveys. Z. Ernahrungswiss 37, 234–241. doi: 10.1007/pl00007377

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: CYP3A, CYP3A4, CYP3A5, midazolam, 4ß-OH-cholesterol, 6ß-OH-cortisol, heritability, twin study

Citation: Matthaei J, Bonat WH, Kerb R, Tzvetkov MV, Strube J, Brunke S, Sachse-Seeboth C, Sehrt D, Hofmann U, von Bornemann Hjelmborg J, Schwab M and Brockmöller J (2020) Inherited and Acquired Determinants of Hepatic CYP3A Activity in Humans. Front. Genet. 11:944. doi: 10.3389/fgene.2020.00944

Received: 10 June 2019; Accepted: 28 July 2020;
Published: 21 August 2020.

Edited by:

Martin A. Kennedy, University of Otago, New Zealand

Reviewed by:

Janne Hukkanen, University of Oulu, Finland
Eng Wee Chua, National University of Malaysia, Malaysia

Copyright © 2020 Matthaei, Bonat, Kerb, Tzvetkov, Strube, Brunke, Sachse-Seeboth, Sehrt, Hofmann, von Bornemann Hjelmborg, Schwab and Brockmöller. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Jürgen Brockmöller, jbrockm@gwdg.de