Genetic predictors of interindividual variability in hepatic CYP3A4 expression

Pharmacokinetics, Pharmacodynamics, and Side Effects of Midazolam: A Review and Case Example

Midazolam, a short-acting benzodiazepine, is widely used to alleviate patient anxiety, enhance compliance, and aid in anesthesia. While its side effects are typically dose-dependent and manageable with vigilant perioperative monitoring, serious cardiorespiratory complications, including fatalities and permanent neurological impairment, have been documented. Prolonged exposure to benzodiazepines, such as midazolam, has been associated with neurological changes in infants. Despite attempts to employ therapeutic drug monitoring for optimal sedation dosing, its efficacy has been limited. Consequently, efforts are underway to identify alternative predictive markers to guide individualized dosing and mitigate adverse effects. Understanding these factors is crucial for determining midazolam's suitability for future administration, particularly after a severe adverse reaction. This article aims to elucidate the factors influencing midazolam's pharmacokinetics and pharmacodynamics, potentially leading to adverse events. Finally, a case study is presented to exemplify the complex investigation into the causative factors of midazolam-related adverse events.

Circulating Biomarkers Instead of Genotyping to Establish Metabolizer Phenotypes

Pharmacogenomics (PGx) enables personalized treatment for the prediction of drug response and to avoid adverse drug reactions. Currently, PGx mainly relies on the genetic information of absorption, distribution, metabolism, and excretion (ADME) targets such as drug-metabolizing enzymes or transporters to predict differences in the patient's phenotype. However, there is evidence that the phenotype-genotype concordance is limited. Thus, we discuss different phenotyping strategies using exogenous xenobiotics (e.g., drug cocktails) or endogenous compounds for phenotype prediction. In particular, minimally invasive approaches focusing on liquid biopsies offer great potential to preemptively determine metabolic and transport capacities. Early studies indicate that ADME phenotyping using exosomes released from the liver is reliable. In addition, pharmacometric modeling and artificial intelligence improve phenotype prediction. However, further prospective studies are needed to demonstrate the clinical utility of individualized treatment based on phenotyping strategies, not only relying on genetics. The present review summarizes current knowledge and limitations.

Sexual Dimorphism in the Expression of Cytochrome P450 Enzymes in Rat Heart, Liver, Kidney, Lung, Brain, and Small Intestine

Cytochrome P450 (P450) enzymes are monooxygenases that are expressed hepatically and extrahepatically and play an essential role in xenobiotic metabolism. Substantial scientific evidence indicates sex-specific differences between males and females in disease patterns and drug responses, which could be attributed, even partly, to differences in the expression and/or activity levels of P450 enzymes in different organs. In this study, we compared the mRNA and protein expression of P450 enzymes in different organs of male and female Sprague-Dawley rats by real-time polymerase chain reaction and western blot techniques. We found significant sex- and organ-specific differences in several enzymes. Hepatic Cyp2c11, Cyp2c13, and Cyp4a2 showed male-specific expression, whereas Cyp2c12 showed female-specific expression. Cyp2e1 and Cyp4f enzymes demonstrated higher expression in the female heart and kidneys compared with males; however, they showed no significant sexual dimorphism in the liver. Male rats showed higher hepatic and renal Cyp1b1 levels. All assessed enzymes were found in the liver, but some were not expressed in other organs. At the protein expression level, CYP1A2, CYP3A, and CYP4A1 demonstrated higher expression levels in the females in several organs, including the liver. Elucidating sex-specific differences in P450 enzyme levels could help better understand differences in disease pathogeneses and drug responses between males and females and thus improve treatment strategies. SIGNIFICANCE STATEMENT: This study characterized the differences in the mRNA and protein expression levels of different cytochrome P450 (P450) enzymes between male and female rats in the heart, liver, lung, kidney, brain, and small intestine. It demonstrated unique sex-specific differences in the different organs. This study is considered a big step towards elucidating sex-specific differences in P450 enzyme levels, which is largely important for achieving a better understanding of the differences between males and females in the disease's processes and treatment outcomes.

Genome-wide microRNA profiles identify miR-107 as a top miRNA associating with expression of the CYP3As and other drug metabolizing cytochrome P450 enzymes in the liver

Cytochrome P450 (CYP) drug metabolizing enzymes are responsible for the metabolism of over 70% of currently used medications with the CYP3A family being the most important CYP enzymes in the liver. Large inter-person variability in expression/activity of the CYP3As greatly affects drug exposure and treatment outcomes, yet the cause of such variability remains elusive. Micro-RNAs (miRNAs) are small noncoding RNAs that negatively regulate gene expression and are involved in diverse cellular processes including metabolism of xenobiotics and therapeutic outcomes. Target prediction and in vitro functional assays have linked several miRNAs to the control of CYP3A4 expression. Yet, their co-expression with CYP3As in the liver remain unclear. In this study, we used genome-wide miRNA profiling in liver samples to identify miRNAs associated with the expression of the CYP3As. We identified and validated both miR-107 and miR-1260 as strongly associated with the expression of CYP3A4, CYP3A5, and CYP3A43. Moreover, we found associations between miR-107 and nine transcription factors (TFs) that regulate CYP3A expression, with estrogen receptor alpha (ESR1) having the largest effect size. Including ESR1 and the other TFs in the regression model either diminished or abolished the associations between miR-107 and the CYP3As, indicating that the role of miR-107 in CYP3A expression may be indirect and occur through these key TFs. Indeed, testing the other nine CYPs previously shown to be regulated by ESR1 identified similar miR-107 associations that were dependent on the exclusion of ESR1 and other key TFs in the regression model. In addition, we found significant differences in miRNA expression profiles in liver samples between race and sex. Together, our results identify miR-107 as a potential epigenetic regulator that is strongly associated with the expression of many CYPs, likely via impacting the CYP regulatory network controlled by ESR1 and other key TFs. Therefore, both genetic and epigenetic factors that alter the expression of miR-107 may have a broad influence on drug metabolism.

Why We Need to Take a Closer Look at Genetic Contributions to CYP3A Activity

Cytochrome P450 3A (CYP3A) subfamily enzymes are involved in the metabolism of 40% of drugs in clinical use. Twin studies have indicated that 66% of the variability in CYP3A4 activity is hereditary. Yet, the complexity of the CYP3A locus and the lack of distinct drug metabolizer phenotypes has limited the identification and clinical application of CYP3A genetic variants compared to other Cytochrome P450 enzymes. In recent years evidence has emerged indicating that a substantial part of the missing heritability is caused by low frequency genetic variation. In this review, we outline the current pharmacogenomics knowledge of CYP3A activity and discuss potential future directions to improve our genetic knowledge and ability to explain CYP3A variability.

An Integrated Clinical and Genetic Prediction Model for Tacrolimus Levels in Pediatric Solid Organ Transplant Recipients

BACKGROUND

There are challenges in achieving and maintaining therapeutic tacrolimus levels after solid organ transplantation (SOT). The purpose of this genome-wide association study was to generate an integrated clinical and genetic prediction model for tacrolimus levels in pediatric SOT.

METHODS

In a multicenter prospective observational cohort study (2015-2018), children <18 years old at their first SOT receiving tacrolimus as maintenance immunosuppression were included (455 as discovery cohort; 322 as validation cohort). Genotyping was performed using a genome-wide single nucleotide polymorphism (SNP) array and analyzed for association with tacrolimus trough levels during 1-y follow-up.

RESULTS

Genome-wide association study adjusted for clinical factors identified 25 SNPs associated with tacrolimus levels; 8 were significant at a genome-wide level (P < 1.025 × 10-7). Nineteen SNPs were replicated in the validation cohort. After removing SNPs in strong linkage disequilibrium, 14 SNPs remained independently associated with tacrolimus levels. Both traditional and machine learning approaches selected organ type, age at transplant, rs776746, rs12333983, and rs12957142 SNPs as the top predictor variables for dose-adjusted 36- to 48-h posttacrolimus initiation (T1) levels. There was a significant interaction between age and organ type with rs776476*1 SNP (P < 0.05). The combined clinical and genetic model had lower prediction error and explained 30% of the variation in dose-adjusted T1 levels compared with 18% by the clinical and 12% by the genetic only model.

CONCLUSIONS

Our study highlights the importance of incorporating age, organ type, and genotype in predicting tacrolimus levels and lays the groundwork for developing an individualized age and organ-specific genotype-guided tacrolimus dosing algorithm.

CYP 450 enzymes influence (R,S)-ketamine brain delivery and its antidepressant activity

Esketamine, the S-stereoisomer of (R,S)-ketamine was recently approved by drug agencies (FDA, EMA), as an antidepressant drug with a new mechanism of action. (R,S)-ketamine is a N-methyl-d-aspartate receptor (NMDA-R) antagonist putatively acting on GABAergic inhibitory synapses to increase excitatory synaptic glutamatergic neurotransmission. Unlike monoamine-based antidepressants, (R,S)-ketamine exhibits rapid and persistent antidepressant activity at subanesthetic doses in preclinical rodent models and in treatment-resistant depressed patients. Its major brain metabolite, (2R,6R)-hydroxynorketamine (HNK) is formed following (R,S)-ketamine metabolism by various cytochrome P450 enzymes (CYP) mainly activated in the liver depending on routes of administration [e.g., intravenous (largely used for a better bioavailability), intranasal spray, intracerebral, subcutaneous, intramuscular or oral]. Experimental or clinical studies suggest that (2R,6R)-HNK could be an antidepressant drug candidate. However, questions still remain regarding its molecular and cellular targets in the brain and its role in (R,S)-ketamine's fast-acting antidepressant effects. The purpose of the present review is: 1) to review (R,S)-ketamine pharmacokinetic properties in humans and rodents and its metabolism by CYP enzymes to form norketamine and HNK metabolites; 2) to provide a summary of preclinical strategies challenging the role of these metabolites by modifying (R,S)-ketamine metabolism, e.g., by administering a pre-treatment CYP inducers or inhibitors; 3) to analyze the influence of sex and age on CYP expression and (R,S)-ketamine metabolism. Importantly, this review describes (R,S)-ketamine pharmacodynamics and pharmacokinetics to alert clinicians about possible drug-drug interactions during a concomitant administration of (R,S)-ketamine and CYP inducers/inhibitors that could enhance or blunt, respectively, (R,S)-ketamine's therapeutic antidepressant efficacy in patients.

Trans-eQTLs of the CYP3A4 and CYP3A5 associated with tacrolimus trough blood concentration in Chinese renal transplant patients

This study aimed to systematically investigate trans-eQTLs of CYP3A4 and CYP3A5 affecting tacrolimus trough blood concentrations in Chinese renal transplant patients. We used Plink v1.90 to perform data quality control and linear regression analysis on GTEx v8 data. SNPs with p-value < 0.05 were selected and the GTEx eQTL Calculator was used to further prioritize the eQTLs of CYP3A4 and CYP3A5 in the liver and small intestine. The eQTLs with a p-value < 5 × 10^-5 and MAF≥ 0.05 in the CHB population were selected as candidate eQTLs. The genotyping of candidate eQTLs was performed using high-resolution melting (HRM) assays and Sanger DNA sequencing. This study included 845 Chinese renal transplant patients who received tacrolimus as an immunosuppressive agent. Association between 103 candidate eQTLs and log-transformed tacrolimus concentration/dose ratio (log (C₀/D)) in this cohort was conducted using the SNPassoc package of R software. In the end, a total of 75,632 liver eQTLs of CYP3A4, 69,558 liver eQTLs of CYP3A5, 48,596 small intestine eQTLs of CYP3A4 and 28,616 small intestine eQTLs of CYP3A5 were obtained using the GTEx v8 eQTL Calculator. Of the 103 candidate eQTLs, rs75727207, rs181294422 and rs28522676 were significantly associated with tacrolimus log(C₀/D) in different genetic models. We discovered a substantial number of novel eQTLs of CYP3A4 and CYP3A5 in liver and small intestine, also found that rs75727207, rs181294422 and rs28522676 may affect tacrolimus trough blood concentrations in Chinese renal transplant patients.

Regulation of PXR Function by Coactivator and Corepressor Proteins: Ligand Binding Is Just the Beginning

The pregnane X receptor (PXR, NR1I2) is a nuclear receptor which exerts its regulatory function by heterodimerization with the retinoid-X-receptor α (RXRα, NR2B1) and binding to the promoter and enhancer regions of diverse target genes. PXR is involved in the regulation of drug metabolism and excretion, metabolic and immunological functions and cancer pathogenesis. PXR activity is strongly regulated by the association with coactivator and corepressor proteins. Coactivator proteins exhibit histone acetyltransferase or histone methyltransferase activity or associate with proteins having one of these activities, thus promoting chromatin decondensation and activation of the gene expression. On the contrary, corepressor proteins promote histone deacetylation and therefore favor chromatin condensation and repression of the gene expression. Several studies pointed to clear cell- and ligand-specific differences in the activation of PXR. In this article, we will review the critical role of coactivator and corepressor proteins as molecular determinants of the specificity of PXR-mediated effects. As already known for other nuclear receptors, understanding the complex mechanism of PXR activation in each cell type and under particular physiological and pathophysiological conditions may lead to the development of selective modulators with therapeutic potential.

Different Inductive Effects of Praziquantel Racemate and its Enantiomers on the Enzyme CYP3A4 Mediated by Pregnane X Receptor and its Variants

BACKGROUND

Praziquantel (PZQ), which possesses an asymmetric center, is classified as a pyrazinoisoquinoline and has been the mainstay in the treatment of schistosomiasis since 1980. PZQ undergoes a pronounced first-pass metabolism in the liver through the CYP450 system which could be mediated by nuclear receptors.

OBJECTIVE

The purpose of this study was to investigate the possible different induction effects of CYP3A4 by PZQ racemate and enantiomers via the pregnane X receptor (PXR) and the effect of PXR polymorphism on the induction potency of PZQs.

METHODS

The dual-luciferase reporter gene systems constructed in HepG2 cells were used to measure the abilities of PZQs to induce CYP3A4 expression mediated by PXR. The mRNA and protein levels of CYP3A4 were evaluated by polymerase chain reaction (PCR) and western blotting, respectively.

RESULTS

In HepG2 cells transfected with PXRwt, PXR158, PXR163, PXR370 or PXR403 expression plasmids, PZQ racemate and its enantiomers up-regulated the luciferase activity in a concentration-dependent manner, while reaching saturation after transfected with PXR379 expression plasmids. The mRNA and protein expression of CYP3A4 was effectively activated in PXR-transfected HepG2 cells. The induction ability of CYP3A4 mediated by PXR activation by PZQ racemate and its enantiomers were statistically different between the same PXR group and different PXR groups.

CONCLUSION

The enantioselective induction effects of PZQs on CYP3A4 were related to the enantioselective activations of PXR by PZQs and were influenced by the PXR gene polymorphism. These findings provide a basis for further understanding the enantiomeric metabolism and the variable efficacy of PZQs.

Co-expression of drug metabolizing cytochrome P450 enzymes and estrogen receptor alpha (ESR1) in human liver: racial differences and the regulatory role of ESR1

OBJECTIVES

The function and expression of cytochrome P450 (CYP) drug metabolizing enzymes is highly variable, greatly affecting drug exposure, and therapeutic outcomes. The expression of these enzymes is known to be controlled by many transcription factors (TFs), including ligand-free estrogen receptor alpha (ESR1, in the absence of estrogen). However, the relationship between the expression of ESR1, other TFs, and CYP enzymes in human liver is still unclear.

METHODS

Using real-time PCR, we quantified the mRNA levels of 12 CYP enzymes and nine TFs in 246 human liver samples from European American (EA, n = 133) and African American (AA, n = 113) donors.

RESULTS

Our results showed higher expression levels of ESR1 and six CYP enzymes in EA than in AA. Partial least square regression analysis showed that ESR1 is the top-ranking TF associating with the expression of eight CYP enzymes, six of which showed racial difference in expression. Conversely, four CYP enzymes without racial difference in expression did not have ESR1 as a top-ranking TF. These results indicate that ESR1 may contribute to variation in CYP enzyme expression between these two ancestral backgrounds.

CONCLUSIONS

These results are consistent with our previous study showing ESR1 as a master regulator for the expression of several CYP enzymes. Therefore, factors affecting ESR1 expression may have broad influence on drug metabolism through altered expression of CYP enzymes.

Relative addictive potential of opioid analgesic agents

Opioid overdoses and deaths continue to be a problem in the USA with a significant portion related to prescribed opioid analgesic agents. The role of pharmacogentic factors in opioid addiction is an active area of research. While all opioid analgesic agents have some addictive potential, it is clear that there are some with greater addictive potential. Oxycodone is the most widely abused opioid analgesic and it appears to predispose to chronic use with high likability by users. Fentanyl and hydromorphone are both very lipophilic allowing rapid penetration into the CNS, but are not rated as highly as other agents. Providers should consider the risk of addiction with the opioids they prescribe and give those with a lower addictive potential.

Gut Microbiome Alters the Activity of Liver Cytochromes P450 in Mice With Sex-Dependent Differences

Sexual differences and the composition/function of the gut microbiome are not considered the most important players in the drug metabolism field; however, from the recent data it is obvious that they may significantly affect the response of the patient to therapy. Here, we evaluated the effect of microbial colonization and sex differences on mRNA expression and the enzymatic activity of hepatic cytochromes P450 (CYPs) in germ-free (GF) mice, lacking the intestinal flora, and control specific-pathogen-free (SPF) mice. We observed a significant increase in the expression of Cyp3a11 in female SPF mice compared to the male group. However, the sex differences were erased in GF mice, and the expression of Cyp3a11 was about the same in both sexes. We have also found higher Cyp2c38 gene expression in female mice compared to male mice in both the SPF and GF groups. Moreover, these changes were confirmed at the level of enzymatic activity, where the female mice exhibit higher levels of functional CYP2C than males in both groups. Interestingly, we observed the same trend as with CYP3A enzymes: a diminished difference between the sexes in GF mice. The presented data indicate that the mouse gut microbiome plays an important role in sustaining sexual dimorphism in terms of hepatic gene expression and metabolism.

Pharmacogenomics as a Tool to Limit Acute and Long-Term Adverse Effects of Chemotherapeutics: An Update in Pediatric Oncology

In the past decades, new cancer treatments have been introduced in pediatric oncology leading to improvement in clinical outcomes and survival rates. However, due to inter-individual differences, some children experience severe chemotherapy-induced toxicities or a poor clinical outcome. An explanation for the diversity in response to chemotherapy is genetic variation, leading to differences in expression and activity of metabolizing and transport enzymes as well as drug targets. Pharmacogenetic testing has emerged as a promising tool to predict and limit acute and long-term adverse effects in patients. However, in pediatric oncology, limited number of patients and a considerable diversity in study results complicate the interpretation of test results and its clinical relevance. With this review, we provide an overview of new developments over the past four years regarding relevant polymorphisms related to toxicity in pediatric oncology. The following chemotherapeutics and associated toxicities are discussed: alkylating agents, anthracyclines, asparaginase, methotrexate, platinum compounds, steroids, thiopurines, topoisomerase inhibitors, and vinca alkaloids. Our review identifies several questions regarding the role of genetic variants in chemotherapy-induced toxicities. Ambiguities in the literature stem from small population sizes, differences in (statistical) interpretation and variations in sequencing technologies as well as different clinical outcome definitions. Standardization of clinical outcome data and toxicity definitions within electronic health records combined with the increased availability of genomic sequence techniques in clinical practice will help to validate these models in upcoming years.

Transcriptional regulation of CYP3A4 by nuclear receptors in human hepatocytes under hypoxia

The human hepatic cytochrome P-450 3A4 (CYP3A4), recognized as a multifunctional enzyme, has a wide range of substrates including commonly used drugs. Previous investigations demonstrated that the expression of CYP3A4 in human hepatocytes could be regulated by some nuclear receptors (NRs) at transcriptional level under diverse situations. The significance of oxygen on CYP3A4-mediated metabolism seems notable while the regulatory mode of CYP3A4 in the particular case still remains elusive. Recently, striking evidence has emerged that both CYP3A4 and its regulator NR could be inhibited by exposure to hypoxia. Therefore, it is of great importance to elucidate whether and how these NRs act in the transcriptional regulation of CYP3A4 in human hepatocytes under hypoxic conditions. In this review, we mainly summarized transcriptional regulation of the pivotal enzyme CYP3A4 by NRs and explored the possible regulatory pathways of CYP3A4 via these major NRs under hypoxia, expecting to provide favorable evidence for further clinical guidance under such pathological situations.

Influence of POR*28 Polymorphisms on CYP3A5*3-Associated Variations in Tacrolimus Blood Levels at an Early Stage after Liver Transplantation

It is well known that the CYP3A5*3 polymorphism is an important marker that correlates with the tacrolimus dose requirement after organ transplantation. Recently, it has been revealed that the POR*28 polymorphism affects the pharmacokinetics of tacrolimus in renal transplant patients. In this study, we examined whether POR*28 as well as CYP3A5*3 polymorphism in Japanese recipients and donors would be another biomarker for the variation of tacrolimus blood levels in the recipients during the first month after living-donor liver transplantation. We enrolled 65 patients treated with tacrolimus, who underwent liver transplantation between July 2016 and January 2019. Genomic DNA was extracted from whole-blood samples, and genotyping was performed to examine the presence of CYP3A5*3 and POR*28 polymorphisms in the recipients and donors. The CYP3A5*3/*3 genotype (defective CYP3A5) of the recipient (standard partial regression coefficient [median C/D ratio of CYP3A5 expressor vs. CYP3A5 non-expressor, p value]: Pod 1–7, β= −0.389 [1.76 vs. 2.73, p < 0.001]; Pod 8–14, β = −0.345 [2.03 vs. 2.83, p < 0.001]; Pod 15–21, β= −0.417 [1.75 vs. 2.94, p < 0.001]; Pod 22–28, β = −0.627 [1.55 vs. 2.90, p < 0.001]) rather than donor (Pod 1–7, β = n/a [1.88 vs. 2.76]; Pod 8–14, β = n/a [1.99 vs. 2.93]; Pod 15–21, β = −0.175 [1.91 vs. 2.94, p = 0.004]; Pod 22–28, β = n/a [1.61 vs. 2.67]) significantly contributed to the increase in the concentration/dose (C/D) ratio of tacrolimus for at least one month after surgery. We found that the tacrolimus C/D ratio significantly decreased from the third week after transplantation when the recipient carried both CYP3A5*1 (functional CYP3A5) and POR*28 (n = 19 [29.2%], median C/D ratio [inter quartile range] = 1.58 [1.39–2.17]), compared with that in the recipients carrying CYP3A5*1 and POR*1/*1 (n = 8 [12.3%], median C/D ratio [inter quartile range] = 2.23 [2.05–3.06]) (p < 0.001). In conclusion, to our knowledge, this is the first report suggesting that the POR*28 polymorphism is another biomarker for the tacrolimus oral dosage after liver transplantation in patients carrying CYP3A5*1 rather than CYP3A5*3/*3.

Identification and analysis of 35 genes associated with vitamin D deficiency: A systematic review to identify genetic variants

Vitamin D deficiency is a public health concern associated with, but not limited to, skeletal anomalies, chronic diseases, immune conditions, and cancer, among others. Hypovitaminosis D is mainly associated with environmental and lifestyle factors that affect sunlight exposure. However, genetic factors also influence 25-hydroxyvitamin D (25[OH]D) serum concentration. Although there is available information of genes with clear biological relevance or markers identified by Genome-Wide Association Studies, an overall view and screening tool to identify known genetic causes of altered serum levels of 25(OH)D is lacking. Moreover, there are no studies including the total genetic evidence associated with abnormal serum concentration of 25(OH)D. Therefore, we conducted a de-novo systematic literature review to propose a set of genes comprehensive of all genetic variants reported to be associated with deficiency of vitamin D. Abstracts retrieved from PubMed search were organized by gene and curated one-by-one using the PubTerm web tool. The genes identified were classified according to the type of genetic evidence associated with serum 25(OH)D levels and were also compared with the few commonly screened genes related to vitamin D status. This strategy allowed the identification of 35 genes associated with serum 25(OH)D concentrations, 27 (75%) of which are not commercially available and are not, therefore, analyzed in clinical practice for genetic counseling, nor are they sufficiently studied for research purposes. Functional analysis of the genes identified confirmed their role in vitamin D pathways and diseases. Thus, the list of genes is an important source to understand the genetic determinants of 25(OH)D levels. To further support our findings, we provide a map of the reported functional variants and SNPs not included in ClinVar, minor allelic frequencies, SNP effect sizes, associated diseases, and an integrated overview of the biological role of the genes. In conclusion, we identified a comprehensive candidate list of genes associated with serum 25(OH)D concentrations, most of which are not commercially available, but would prove of importance in clinical practice in screening for patients that should respond to supplementation because of alterations in absorption, patients that would have little benefit because alterations in the downstream metabolism of vitamin D, and to study non-responsiveness to supplementation with vitamin D.

Ligand-Free Estrogen Receptor α (ESR1) as Master Regulator for the Expression of CYP3A4 and Other Cytochrome P450 Enzymes in the Human Liver

Cytochrome P450 3A4 isoform (CYP3A4) transcription is controlled by hepatic transcription factors (TFs), but how TFs dynamically interact remains uncertain. We hypothesize that several TFs form a regulatory network with nonlinear, dynamic, and hierarchical interactions. To resolve complex interactions, we have applied a computational approach for estimating Sobol's sensitivity indices (SSI) under generalized linear models to existing liver RNA expression microarray data (GSE9588) and RNA-seq data from genotype-tissue expression (GTEx), generating robust importance ranking of TF effects and interactions. The SSI-based analysis identified TFs and interacting TF pairs, triplets, and quadruplets involved in CYP3A4 expression. In addition to known CYP3A4 TFs, estrogen receptor α (ESR1) emerges as key TF with the strongest main effect and as the most frequently included TF interacting partner. Model predictions were validated using small interfering RNA (siRNA)/short hairpin RNA (shRNA) gene knockdown and clustered regularly interspaced short palindromic repeats (CRISPR)-mediated transcriptional activation of ESR1 in biliary epithelial Huh7 cells and human hepatocytes in the absence of estrogen. Moreover, ESR1 and known CYP3A4 TFs mutually regulate each other. Detectable in both male and female hepatocytes without added estrogen, the results demonstrate a role for unliganded ESR1 in CYP3A4 expression consistent with unliganded ESR1 signaling reported in other cell types. Added estrogen further enhances ESR1 effects. We propose a hierarchical regulatory network for CYP3A4 expression directed by ESR1 through self-regulation, cross regulation, and TF-TF interactions. We also demonstrate that ESR1 regulates the expression of other P450 enzymes, suggesting broad influence of ESR1 on xenobiotics metabolism in human liver. Further studies are required to understand the mechanisms underlying role of ESR1 in P450 regulation. SIGNIFICANCE STATEMENT: This study focuses on identifying key transcription factors and regulatory networks for CYP3A4, the main drug metabolizing enzymes in liver. We applied a new computational approach (Sobol's sensitivity analysis) to existing hepatic gene expression data to determine the role of transcription factors in regulating CYP3A4 expression, and used molecular genetics methods (siRNA/shRNA gene knockdown and CRISPR-mediated transcriptional activation) to test these interactions in life cells. This approach reveals a robust network of TFs, including their putative interactions and the relative impact of each interaction. We find that ESR1 serves as a key transcription factor function in regulating CYP3A4, and it appears to be acting at least in part in a ligand-free fashion.

Nuclear receptor AHR-mediated xenobiotic detoxification pathway involves in atrazine-induced nephrotoxicity in quail (Coturnix C. coturnix)

Atrazine (ATR), one of the most widely used pesticides in agricultural production, are gradually concerned due to potential ecosystem and health risks. Further, the induction of ATR nephrotoxicity and detoxification response is still unknown. To evaluate ATR-induced nephrotoxicity, quails were treated with 0, 50, 250 or 500 mg/kg ATR by gavage administration for 45 days. Histopathology indicated that ATR exposure caused renal tubular epithelial cell swelling and endoplasmic reticulum degeneration, suggesting that ATR exposure causes renal impairment even renal diseases. Notably, ATR interfered cytochrome P450 system (CYP450s) homeostasis by enhancing contents or activities of CYP450s (total CYP450, Cyt b5, AH, APND, NCR and ERND) and the expression of CYP450 isoforms (CYP1A, CYP1B, CYP2C and CYP3A). ATR triggered phase II detoxifying reaction, reflected by the elevated GSH level, GST activity and the up-regulation of GST isoforms (GSTa, GSTa3 and GSTt1) and GSH synthetase (GCLC). Moreover, ABC transporters were activated to expel ATR from the body by increasing expression of MRP1 and P-GP gene. Accompanying these alterations, the nuclear receptors (AHR, CAR and PXR) were activated by ATR in a dose-dependent manner. Analysis results of present study demonstrated that the induction of phase II detoxifying enzyme system and ABC transporters could be modulated by nuclear receptors response and CYP450s disturbance in low-dose ATR-treated quail. In conclusion, all data suggested that nuclear receptors AHR-mediated detoxification pathway was involved in ATR-induced nephrotoxicity. These results provided new evidence about the nephrotoxic effects of ATR on the response of biotransformation and detoxification system.

Integrating physiologically based kinetic (PBK) and Monte Carlo modelling to predict inter-individual and inter-ethnic variation in bioactivation and liver toxicity of lasiocarpine

The aim of the present study was to predict the effect of inter-individual and inter-ethnic human kinetic variation on the sensitivity towards acute liver toxicity of lasiocarpine in the Chinese and the Caucasian population, and to derive chemical specific adjustment factors (CSAFs) by integrating variation in the in vitro kinetic constants V_max and K_m, physiologically based kinetic (PBK) modelling and Monte Carlo simulation. CSAFs were derived covering the 90th and 99th percentile of the population distribution of pyrrole glutathione adduct (7-GS-DHP) formation, reflecting bioactivation. The results revealed that in the Chinese population, as compared to the Caucasian population, the predicted 7-GS-DHP formation at the geometric mean, the 90th and the 99th percentile were 2.1-, 3.3- and 4.3-fold lower respectively. The CSAFs obtained using the 99th percentile values were 8.3, 17.0 and 19.5 in the Chinese, the Caucasian population and the two populations combined, respectively, while the CSAFs were generally 3.0-fold lower at the 90th percentile. These results indicate that when considering the formation of 7-GS-DHP the Caucasian population may be more sensitive towards acute liver toxicity of lasiocarpine, and further point out that the default safety factor of 3.16 for inter-individual human kinetic differences may not be sufficiently protective. Altogether, the results obtained demonstrate that integrating PBK modelling with Monte Carlo simulations using human in vitro data is a powerful strategy to quantify inter-individual variations in kinetics, and can be used to refine the human risk assessment of pyrrolizidine alkaloids.

Interindividual Variation in CYP3A Activity Influences Lapatinib Bioactivation

Lapatinib is a dual tyrosine kinase inhibitor associated with rare but potentially severe idiosyncratic hepatotoxicity. We have previously shown that cytochromes P450 CYP3A4 and CYP3A5 quantitatively contribute to lapatinib bioactivation, leading to formation of a reactive, potentially toxic quinone imine. CYP3A5 is highly polymorphic; however, the impact of CYP3A5 polymorphism on lapatinib metabolism has not been fully established. The goal of this study was to determine the effect of CYP3A5 genotype and individual variation in CYP3A activity on the metabolic activation of lapatinib using human-relevant in vitro systems. Lapatinib metabolism was examined using CYP3A5-genotyped human liver microsomes and cryopreserved human hepatocytes. CYP3A and CYP3A5-selective activities were measured in liver tissues using probe substrates midazolam and T-5 (T-1032), respectively, to evaluate the correlation between enzymatic activity and lapatinib metabolite formation. Drug metabolites were measured by high-performance liquid chromatography-tandem mass spectrometry. Further, the relative contributions of CYP3A4 and CYP3A5 to lapatinib O-debenzylation were estimated using selective chemical inhibitors of CYP3A. The results from this study demonstrated that lapatinib O-debenzylation and quinone imine-GSH conjugate formation were highly correlated with hepatic CYP3A activity, as measured by midazolam 1'-hydroxylation. CYP3A4 played a dominant role in lapatinib bioactivation in all liver tissues evaluated. The CYP3A5 contribution to lapatinib bioactivation varied by individual donor and was dependent on CYP3A5 genotype and activity. CYP3A5 contributed approximately 20%-42% to lapatinib O-debenzylation in livers from CYP3A5 expressers. These findings indicate that individual CYP3A activity, not CYP3A5 genotype alone, is a key determinant of lapatinib bioactivation and likely influences exposure to reactive metabolites. SIGNIFICANCE STATEMENT: This study is the first to examine the effect of CYP3A5 genotype, total CYP3A activity, and CYP3A5-selective activity on lapatinib bioactivation in individual human liver tissues. The results of this investigation indicate that lapatinib bioactivation via oxidative O-debenzylation is highly correlated with total hepatic CYP3A activity, and not CYP3A5 genotype alone. These findings provide insight into the individual factors, namely, CYP3A activity, that may affect individual exposure to reactive, potentially toxic metabolites of lapatinib.

Perspectives from the Innovation and Quality Consortium Induction Working Group on Factors Impacting Clinical Drug-Drug Interactions Resulting from Induction: Focus on Cytochrome 3A Substrates

A recent publication from the Innovation and Quality Consortium Induction Working Group collated a large clinical data set with the goal of evaluating the accuracy of drug-drug interaction (DDI) prediction from in vitro data. Somewhat surprisingly, comparison across studies of the mean- or median-reported area under the curve ratio showed appreciable variability in the magnitude of outcome. This commentary explores the possible drivers of this range of outcomes observed in clinical induction studies. While recommendations on clinical study design are not being proposed, some key observations were informative during the aggregate analysis of clinical data. Although DDI data are often presented using median data, individual data would enable evaluation of how differences in study design, baseline expression, and the number of subjects contribute. Since variability in perpetrator pharmacokinetics (PK) could impact the overall DDI interpretation, should this be routinely captured? Maximal induction was typically observed after 5-7 days of dosing. Thus, when the half-life of the inducer is less than 30 hours, are there benefits to a more standardized study design? A large proportion of CYP3A4 inducers were also CYP3A4 inhibitors and/or inactivators based on in vitro data. In these cases, using CYP3A selective substrates has limitations. More intensive monitoring of changes in area under the curve over time is warranted. With selective CYP3A substrates, the net effect was often inhibition, whereas less selective substrates could discern induction through mechanisms not susceptible to inhibition. The latter included oral contraceptives, which raise concerns of reduced efficacy following induction. Alternative approaches for modeling induction, such as applying biomarkers and physiologically based pharmacokinetic modeling (PBPK), are also considered. SIGNIFICANCE STATEMENT: The goal of this commentary is to stimulate discussion on whether there are opportunities to optimize clinical drug-drug interaction study design. The overall aim is to reduce, understand and contextualize the variability observed in the magnitude of induction across reported clinical studies. A large clinical CYP3A induction dataset was collected and further analyzed to identify trends and gaps. Reporting individual victim PK data, characterizing perpetrator PK and including additional PK assessments for mixed-mechanism perpetrators may provide insights into how these factors impact differences observed in clinical outcomes. The potential utility of biomarkers and PBPK modeling are discussed in considering future directions.

Identification of the caffeine to trimethyluric acid ratio as a dietary biomarker to characterise variability in cytochrome P450 3A activity

PURPOSE

Cytochrome P450 (CYP) 3A plays an important role in the metabolism of many clinically used drugs and exhibits substantial between-subject variability (BSV) in activity. Current methods to assess variability in CYP3A activity have limitations and there remains a need for a minimally invasive clinically translatable strategy to define CYP3A activity. The purpose of this study was to evaluate the potential for a caffeine metabolic ratio to describe variability in CYP3A activity.

METHODS

The metabolic ratio 1,3,7-trimethyluric acid (TMU) to caffeine was evaluated as a biomarker to describe variability in CYP3A activity in a cohort (n = 28) of healthy 21 to 35-year-old males. Midazolam, caffeine, and TMU concentrations were assessed at baseline and following dosing of rifampicin (300 mg daily) for 7 days.

RESULTS

At baseline, correlation coefficients for the relationship between apparent oral midazolam clearance (CL/F) with caffeine/TMU ratio measured at 3, 4, and 6 h post dose were 0.82, 0.79, and 0.65, respectively. The strength of correlations was retained post rifampicin dosing; 0.72, 0.87, and 0.82 for the ratios at 3, 4, and 6 h, respectively. Weaker correlations were observed between the change in midazolam CL/F and change in caffeine/TMU ratio post/pre-rifampicin dosing.

CONCLUSION

BSV in CYP3A activity was well described by caffeine/TMU ratios pre- and post-induction. The caffeine/TMU ratio may be a convenient tool to assess BSV in CYP3A activity, but assessment of caffeine/TMU ratio alone is unlikely to account for all sources of variability in CYP3A activity.

Effect of Grapefruit Juice Intake on Serum Level of the Endogenous CYP3A4 Metabolite 4β-Hydroxycholesterol-an Interaction Study in Healthy Volunteers

4β-Hydroxycholesterol (4βOHC) is an endogenous CYP3A4 metabolite. However, it is unclear whether circulating levels of 4βOHC may reflect hepatic CYP3A4 activity or both hepatic and intestinal enzyme activity. The aim of this study was to investigate the effect of grapefruit juice, regarded to be a selective intestinal CYP3A4 inhibitor, on serum 4βOHC levels in healthy volunteers. The participants (n = 22) consumed grapefruit juice twice daily for 3 weeks followed by a 2-week washout period. Blood samples for measurements of 4βOHC and the non-CYP3A4-derived oxysterols 24-hydroxycholesterol (24OHC) and 27-hydroxycholesterol (27OHC), as well as lathosterol and total cholesterol, were drawn on days 0, 7, 21, and 35. Median individual changes (ratios) in cholesterol-corrected 4βOHC levels from baseline to weeks 1, 3, and 5 were 0.94 (P = 0.2), 0.98 (P = 0.3), and 0.97 (P = 0.9), respectively. In comparison, median changes (ratios) in cholesterol-corrected levels of 24OHC at the same points were 1.01 (P = 0.6), 0.98 (P = 0.3), and 0.99 (P = 0.5), and of 27OHC 1.01 (P = 0.8), 0.97 (P = 0.5), and 0.99 (P = 0.2). Surprisingly, serum concentration of cholesterol was significantly reduced by approximately 5% after 1 week (P = 0.03), while median cholesterol-corrected levels of lathosterol increased significantly and persistently by approximately 15% during the whole 5-week period (P < 0.04). In conclusion, the present findings suggest that intestinal CYP3A4 is not relevant for the overall formation of 4βOHC in healthy volunteers. The fact that grapefruit juice altered cholesterol homeostasis should be further investigated.

Identification of genetic variants associated with tacrolimus metabolism in kidney transplant recipients by extreme phenotype sampling and next generation sequencing

An extreme phenotype sampling (EPS) model with targeted next-generation sequencing (NGS) identified genetic variants associated with tacrolimus (Tac) metabolism in subjects from the Deterioration of Kidney Allograft Function (DeKAF) Genomics cohort which included 1,442 European Americans (EA) and 345 African Americans (AA). This study included 48 subjects separated into 4 groups of 12 (AA high, AA low, EA high, EA low). Groups were selected by the extreme phenotype of dose-normalized Tac trough concentrations after adjusting for common genetic variants and clinical factors. NGS spanned >3 Mb of 28 genes and identified 18,661 genetic variants (3,961 previously unknown). A group of 125 deleterious variants, by SIFT analysis, were associated with Tac troughs in EAs (burden test, p=0.008), CYB5R2 was associated with Tac troughs in AAs (SKAT, p=0.00079). In CYB5R2, rs61733057 (increased allele frequency in AAs) was predicted to disrupt protein function by SIFT and PolyPhen2 analysis. The variants merit further validation.

Mathematical Models in the Description of Pregnane X Receptor (PXR)-Regulated Cytochrome P450 Enzyme Induction

The pregnane X receptor (PXR) is a drug/xenobiotic-activated transcription factor of crucial importance for major cytochrome P450 xenobiotic-metabolizing enzymes (CYP) expression and regulation in the liver and the intestine. One of the major target genes regulated by PXR is the cytochrome P450 enzyme (CYP3A4), which is the most important human drug-metabolizing enzyme. In addition, PXR is supposed to be involved both in basal and/or inducible expression of many other CYPs, such as CYP2B6, CYP2C8, 2C9 and 2C19, CYP3A5, CYP3A7, and CYP2A6. Interestingly, the dynamics of PXR-mediated target genes regulation has not been systematically studied and we have only a few mechanistic mathematical and biologically based models describing gene expression dynamics after PXR activation in cellular models. Furthermore, few indirect mathematical PKPD models for prediction of CYP3A metabolic activity in vivo have been built based on compartmental models with respect to drug–drug interactions or hormonal crosstalk. Importantly, several negative feedback loops have been described in PXR regulation. Although current mathematical models propose these adaptive mechanisms, a comprehensive mathematical model based on sufficient experimental data is still missing. In the current review, we summarize and compare these models and address some issues that should be considered for the improvement of PXR-mediated gene regulation modelling as well as for our better understanding of the quantitative and spatial dynamics of CYPs expression.

Nonalcoholic Fatty Liver Disease and Diabetes Are Associated with Decreased CYP3A4 Protein Expression and Activity in Human Liver

Nonalcoholic fatty liver disease (NAFLD) is a major cause of chronic liver disease in the Western population. We investigated the association of nonalcoholic fatty liver disease (NAFLD) and diabetes mellitus on CYP3A4 activity in human liver tissue from brain dead donors ( n = 74). Histopathologically graded livers were grouped into normal ( n = 24), nonalcoholic fatty liver (NAFL, n = 26), and nonalcoholic steatohepatitis (NASH, n = 24) categories. The rate of conversion of midazolam to its 1-hydroxy metabolite was used to assess in vitro CYP3A4 activity in human liver microsomes (HLM). A proteomics approach was utilized to quantify the protein expression of CYP3A4 and related enzymes. Moreover, a physiologically based pharmacokinetic (PBPK) model was developed to allow prediction of midazolam concentration in NAFL and NASH livers. CYP3A4 activity in NAFL and NASH was 1.9- and 3.1-fold ( p < 0.05) lower than normal donors, respectively. Intrinsic clearance (CLint) was 2.7- ( p < 0.05) and 4.1-fold ( p < 0.01) lower in donors with NAFL and NASH, respectively. CYP3A4 protein expression was significantly lower in NAFL and NASH donors ( p < 0.05) and accounted for significant midazolam hydroxylation variability in a multiple linear regression analysis (β = 0.869, r² = 0.762, P < 0.01). Diabetes was also associated with decreased CYP3A4 activity and protein expression. Both midazolam CLint and CYP3A4 protein abundance decreased significantly with increase in hepatic fat accumulation. Age and gender did not exhibit any significant association with the observed alterations. Predicted midazolam exposure was 1.7- and 2.3-fold higher for NAFL and NASH, respectively, which may result in a longer period of sedation in these disease-states. Data suggests that NAFLD and diabetes are associated with the decreased hepatic CYP3A4 activity. Thus, further evaluation of clinical consequences of these findings on the efficacy and safety of CYP3A4 substrates is warranted.

Relationships between the expression of hepatocyte nuclear factors and factors essential for lipoprotein production in a human mesenchymal stem cell line, UE7T-13

To clarify the mechanisms regulating lipoprotein production by hepatocyte nuclear factors (HNFs), we generated four kinds of transfectants in human bone marrow mesenchymal stem cells: UE7T-13, stably expressing FOXA2 (also known as HNF3β), HNF4α, HNF1α or co-expressing HNF4α, and HNF1α (HNF4α/HNF1α). In HNF4α/HNF1α transfectants, cellular contents of triglycerides (TG) and cholesterol were markedly higher than in UE7T-13 cells and comparable to those in human hepatoma HepG2 cells. However, TG and cholesterol, which are secreted from cells as components of lipoproteins, were hardly detected in the medium for any of the transfectants. ApoB100 and MTP, which are essential for the formation and secretion of lipoproteins, were undetectable and detected at low levels, respectively, in HNF4α/HNF1α transfectants. We suggest that enforced co-expression of HNF4α and HNF1α is effective for cellular lipid accumulation, while additional factors are probably required for lipoprotein formation and secretion.

Impact of genetic and nongenetic factors on interindividual variability in 4β-hydroxycholesterol concentration

PURPOSE

Individual variability in the endogenous CYP3A metabolite 4β-hydroxycholesterol (4βOHC) is substantial, but to which extent this is determined by genetic and nongenetic factors remains unclear. The aim of the study was to evaluate the explanatory power of candidate genetic variants and key nongenetic factors on individual variability in 4βOHC levels in a large naturalistic patient population.

METHODS

We measured 4βOHC concentration in serum samples from 655 patients and used multiple linear regression analysis to estimate the quantitative effects of CYP3A4*22, CYP3A5*3, and POR*28 variant alleles, comedication with CYP3A inducers, inhibitors and substrates, sex, and age on individual 4βOHC levels.

RESULTS

4βOHC concentration ranged >100-fold in the population, and the multiple linear regression model explained about one fourth of the variability (R ² = 0.23). Only comedication with inducers or inhibitors, sex, and POR genotype were significantly associated with individual variability in 4βOHC level. The estimated quantitative effects on 4βOHC levels were greatest for inducer comedication (+>313%, P < 0.001), inhibitor comedication (-34%, P = 0.021), and female sex (+30%, P < 0.001), while only a modestly elevated 4βOHC level was observed in carriers vs. noncarriers of POR*28 (+11%, P = 0.023).

CONCLUSIONS

These findings suggest that the CYP3A4*22, CYP3A5*3, and POR*28 variant alleles are of limited importance for overall individual variability in 4βOHC levels compared to nongenetic factors.

Developmental regulation of CYP3A4 and CYP3A7 in Chinese Han population

CYP3A4 and CYP3A7 are generally served as the major adult and fetal liver forms, respectively, and exhibited a developmental switch during liver maturation. The objective of this study was to explore the potential mechanisms associated with the developmental switch of CYP3A4 and CYP3A7 in the Chinese Han population. We analyzed CYP3A4/7, nuclear receptors, and epigenetic modifications in human liver samples. We found that the expression levels of CYP3A4 mRNA in adults were significantly higher than the levels in fetus. In contrast, CYP3A7 mRNA expression reached a maximal level at an estimated gestational age of 25 weeks and then substantially decreased during the first year after birth. We also found that the expression level of hepatocyte nuclear factor 4 alpha (HNF4A) was most associated with CYP3A4 expression in adult liver; whereas the expression level of glucocorticoid receptor (GR) was intensively correlated with CYP3A7 expression in fetal liver. Furthermore, we illustrated the dynamic changes of H3K4me2 and H3K27me3 in the developmental switch of CYP3A7 and CYP3A4. In summary, our data suggested that HNF4A and GR, and epigenetic changes of H3K4me2 and H3K27me3 are associated with the ontogenic expressions of CYP3A4/3A7 in the livers of the Chinese Han population.

PXR variants: the impact on drug metabolism and therapeutic responses

The pregnane X receptor (PXR) plays an important and diverse role in mediating xenobiotic induction of drug-metabolizing enzymes and transporters. Several protein isoforms of PXR exist, and they have differential transcriptional activity upon target genes; transcript variants 3 (PXR3) and 4 (PXR4) do not induce target gene expression, whereas transcript variants 1 (PXR1) and 2 (PXR2) respond to agonist by activating target gene expression. PXR protein variants also display differences in protein-protein interactions; PXR1 interacts with p53, whereas PXR3 does not. Furthermore, the transcript variants of PXR that encode these protein isoforms are differentially regulated by methylation and deletions in the respective promoters of the variants, and their expression differs in various human cancers and also in cancerous tissue compared to adjacent normal tissues. PXR1 and PXR4 mRNA are downregulated by methylation in cancerous tissue and have divergent effects on cellular proliferation when ectopically overexpressed. Additional detailed and comparative mechanistic studies are required to predict the effect of PXR transcript variant expression on carcinogenesis, therapeutic response, and the development of toxicity.

Bioactivation of Trimethoprim to Protein-Reactive Metabolites in Human Liver Microsomes

The formation of drug-protein adducts via metabolic activation and covalent binding may stimulate an immune response or may result in direct cell toxicity. Protein covalent binding is a potentially pivotal step in the development of idiosyncratic adverse drug reactions (IADRs). Trimethoprim (TMP)-sulfamethoxazole (SMX) is a combination antibiotic that commonly causes IADRs. Recent data suggest that the contribution of the TMP component of TMP-SMX to IADRs may be underappreciated. We previously demonstrated that TMP is bioactivated to chemically reactive intermediates that can be trapped in vitro by N-acetyl cysteine (NAC), and we have detected TMP-NAC adducts (i.e., mercapturic acids) in the urine of patients taking TMP-SMX. However, the occurrence and extent of TMP covalent binding to proteins was unknown. To determine the ability of TMP to form protein adducts, we incubated [(14)C]TMP with human liver microsomes in the presence and absence of NADPH. We observed protein covalent binding that was NADPH dependent and increased with incubation time and concentration of both protein and TMP. The estimated covalent binding was 0.8 nmol Eq TMP/mg protein, which is comparable to the level of covalent binding for several other drugs that have been associated with covalent binding-induced toxicity and/or IADRs. NAC and selective inhibitors of CYP2B6 and CYP3A4 significantly reduced TMP covalent binding. These results demonstrate for the first time that TMP bioactivation can lead directly to protein adduct formation, suggesting that TMP has been overlooked as a potential contributor of TMP-SMX IADRs.

MiRNA-Based Regulation of Hemostatic Factors through Hepatic Nuclear Factor-4 Alpha

MiRNAs have been reported as CIS-acting elements of several hemostatic factors, however, their mechanism as TRANS-acting elements mediated by a transcription factor is little known and could have important effects. HNF4α has a direct and important role in the regulation of multiple hepatic coagulation genes. Previous in vitro studies have demonstrated that miR-24-3p and miR-34a-5p regulate HNF4A expression. Here we aimed to investigate the molecular mechanisms of miR-24 and miR-34a on coagulation through HNF4A. Transfections with miR-24 and miR-34a in HepG2 cells decreased not only HNF4A but also F10, F12, SERPINC1, PROS1, PROC, and PROZ transcripts levels. Positive and significant correlations were observed between levels of HNF4A and several hemostatic factors (F5, F8, F9, F11, F12, SERPINC1, PROC, and PROS1) in human liver samples (N = 104). However, miR-24 and miR-34a levels of the low (10th) and high (90th) percentiles of those liver samples were inversely correlated with HNF4A and almost all hemostatic factors expression levels. These outcomes suggest that miR-24 and miR-34a might be two indirect elements of regulation of several hemostatic factors. Additionally, variations in miRNA expression profiles could justify, at least in part, changes in HNF4A expression levels and its downstream targets of coagulation.

Genome-wide association study of 25(OH) Vitamin D concentrations in Punjabi Sikhs: Results of the Asian Indian diabetic heart study

Vitamin D deficiency is implicated in multiple disease conditions and accumulating evidence supports that the variation in serum vitamin D (25(OH)D) levels, including deficiency, is under strong genetic control. However, the underlying genetic mechanism associated with vitamin 25(OH)D concentrations is poorly understood. We earlier reported a very high prevalence of vitamin D deficiency associated with an increased risk for type 2 diabetes and obesity in a Punjabi Sikh diabetic cohort as part of the Asian Indian diabetic heart study (AIDHS). Here we have performed the first genome-wide association study (GWAS) of serum 25(OH)D on 3538 individuals from this Punjabi Sikh population. Our discovery GWAS comprised of 1387 subjects followed by validation of 24 putative SNPs (P<10(-4)) using an independent replication sample (n=2151) from the same population by direct genotyping. A novel locus at chromosome 20p11.21 represented by rs2207173 with minor allele frequency (MAF) 0.29, [β=-0.13, p=4.47×10(-9)] between FOXA2 and SSTR4 was identified to be associated with 25(OH)D levels. Another suggestive association signal at rs11586313 (MAF 0.54) [β=0.90; p=1.36×10(-6)] was found within the regulatory region of the IVL gene on chromosome 1q21.3. Additionally, our study replicated 3 of 5 known GWAS genes associated with 25(OH)D concentrations including GC (p=0.007) and CYP2R1 (p=0.019) reported in Europeans and the DAB1 (p=0.003), reported in Hispanics. Identification of novel association signals in biologically plausible regions with 25(OH)D metabolism will provide new molecular insights on genetic drivers of vitamin D status and its implications in health disparities.

Factors associated with 1,25-dihydroxyvitamin D3 concentrations in liver transplant recipients: a prospective observational longitudinal study

The aim of the study was to identify factors associated with 1,25(OH)2D3 concentrations in liver transplant recipients with emphasis on the renal function and catabolism. We also tested the hypothesis that tacrolimus increases 1,25(OH)2D3 concentrations. Serum 25(OH)D3, 1,25(OH)2D3, and 24,25(OH)2D3 were measured in 41 patients before, at 2 weeks and 3 months after transplantation. Dose-adjusted tacrolimus concentration was used as a surrogate marker of CYP3A4 activity. Factors associated with 1,25(OH)2D3 were identified using multivariate linear regression analysis. The median 1,25(OH)2D3 levels remained stable: 55 versus 46 pg/ml (P = 0.36) despite an increase in 25(OH)D3 from 18 ng/ml at baseline to 26 ng/ml (P = 0.03), serum albumin (34 to 41 g/l, P = 0.02), and comparable eGFR at baseline and month 3 (94 and 92 ml/min, respectively, P = 0.15). At 3 months 19 % of patients had 1,25(OH)2D3 < 25 pg/ml. Low eGFR and a low dose-adjusted tacrolimus concentration were both independently associated with 1,25(OH)2D3 at 3 months. Liver transplant recipients with impaired renal function or a low dose-adjusted tacrolimus concentration suggesting a high CYP3A4 are at risk of low 1,25(OH)2D3 concentrations. The use of tacrolimus does not lead to an increase in 1,25(OH)2D3 concentrations in a clinical setting.

Interindividual Variability in Cytochrome P450-Mediated Drug Metabolism

The cytochrome P450 (P450) enzymes are the predominant enzyme system involved in human drug metabolism. Alterations in the expression and/or activity of these enzymes result in changes in pharmacokinetics (and consequently the pharmacodynamics) of drugs that are metabolized by this set of enzymes. Apart from changes in activity as a result of drug-drug interactions (by P450 induction or inhibition), the P450 enzymes can exhibit substantial interindividual variation in basal expression and/or activity, leading to differences in the rates of drug elimination and response. This interindividual variation can result from a myriad of factors, including genetic variation in the promoter or coding regions, variation in transcriptional regulators, alterations in microRNA that affect P450 expression, and ontogenic changes due to exposure to xenobiotics during the developmental and early postnatal periods. Other than administering a probe drug or cocktail of drugs to obtain the phenotype or conducting a genetic analysis to determine genotype, methods to determine interindividual variation are limited. Phenotyping via a probe drug requires exposure to a xenobiotic, and genotyping is not always well correlated with phenotype, making both methodologies less than ideal. This article describes recent work evaluating the effect of some of these factors on interindividual variation in human P450-mediated metabolism and the potential utility of endogenous probe compounds to assess rates of drug metabolism among individuals.

Transcriptional Regulation of CYP2B6 Expression by Hepatocyte Nuclear Factor 3β in Human Liver Cells

CYP2B6 plays an increasingly important role in xenobiotic metabolism and detoxification. The constitutive androstane receptor (CAR) and the pregnane X receptor (PXR) have been established as predominant regulators for the inductive expression of CYP2B6 gene in human liver. However, there are dramatic interindividual variabilities in CYP2B6 expression that cannot be fully explained by the CAR/PXR-based modulation alone. Here, we show that expression level of CYP2B6 was correlated with that of hepatocyte nuclear factor 3β (HNF3β) in human primary hepatocytes prepared from 35 liver donors. Utilizing recombinant virus-mediated overexpression or knockdown of HNF3β in HepG2 cells, as well as constructs containing serial deletion and site-directed mutation of HNF3β binding motifs in CYP2B6 luciferase reporter assays, we demonstrated that the presence or lack of HNF3β expression markedly correlated with CYP2B6 gene expression and its promoter activity. Novel enhancer modules of HNF3β located upstream of the CYP2B6 gene transcription start site were identified and functionally validated as key elements governing HNF3β-mediated CYP2B6 expression. Chromatin immunoprecipitation assays in human primary hepatocytes and surface plasmon resonance binding affinity experiments confirmed the essential role of these enhancers in the recruitment of HNF3β to the promoter of CYP2B6 gene. Overall, these findings indicate that HNF3β represents a new liver enriched transcription factor that is involved in the transcription of CYP2B6 gene and contributes to the large interindividual variations of CYP2B6 expression in human population.

Dynamics of Cytosine Methylation in the Proximal Promoters of CYP3A4 and CYP3A7 in Pediatric and Prenatal Livers

Members of the human CYP3A family of metabolizing enzymes exhibit developmental changes in expression whereby CYP3A7 is expressed in fetal tissues, followed by a transition to expression of CYP3A4 in the first months of life. Despite knowledge about the general pattern of CYP3A activity in human development, the mechanisms that regulate developmental expression remain poorly understood. Epigenetic changes, including cytosine methylation, have been suggested to play a role in the regulation of CYP3A expression. The objective of this study was to investigate changes in cytosine methylation of the CYP3A4 and CYP3A7 genes in human pediatric and prenatal livers. The methylation status of cytosine-phospho-guanine dinucleotides was determined in 16 pediatric liver samples using methyl-seq and confirmed by bisulfite sequencing of 48 pediatric and 34 prenatal liver samples. Samples were separated by age into five groups (prenatal, < 1 year of age, 1.8-6 years, 7-11 years, and 12-17 years). Methyl-seq anaylsis revealed that cytosines in the proximal promoter of CYP3A7 are hypomethylated in neonates compared with adolescents (P < 0.001). In contrast, a cytosine 383 base pair upstream of CYP3A4 is hypermethylated in liver samples from neonates compared with adolescents (P = 0.00001). Developmental changes in methylation of cytosines in the proximal promoters of CYP3A4 and CYP3A7 in pediatric livers were confirmed by bisulfite sequencing. In addition, the methylation status of cytosine in the CYP3A4 and CYP3A7 proximal promoters correlated with changes in developmental expression of mRNA for the two enzymes.

Population pharmacokinetic approach to evaluate the effect of CYP2D6, CYP3A, ABCB1, POR and NR1I2 genotypes on donepezil clearance

AIMS

A large interindividual variability in plasma concentrations has been reported in patients treated with donepezil, the most frequently prescribed antidementia drug. We aimed to evaluate clinical and genetic factors influencing donepezil disposition in a patient population recruited from a naturalistic setting.

METHODS

A population pharmacokinetic study was performed including data from 129 older patients treated with donepezil. The patients were genotyped for common polymorphisms in the metabolic enzymes CYP2D6 and CYP3A, in the electron transferring protein POR and the nuclear factor NR1I2 involved in CYP activity and expression, and in the drug transporter ABCB1.

RESULTS

The average donepezil clearance was 7.3 l h(-1) with a 30% interindividual variability. Gender markedly influenced donepezil clearance (P < 0.01). Functional alleles of CYP2D6 were identified as unique significant genetic covariate for donepezil clearance (P < 0.01), with poor metabolizers and ultrarapid metabolizers demonstrating, respectively, a 32% slower and a 67% faster donepezil elimination compared with extensive metabolizers.

CONCLUSION

The pharmacokinetic parameters of donepezil were well described by the developed population model. Functional alleles of CYP2D6 significantly contributed to the variability in donepezil disposition in the patient population and should be further investigated in the context of individual dose optimization to improve clinical outcome and tolerability of the treatment.

Co-modulation analysis of gene regulation in breast cancer reveals complex interplay between ESR1 and ERBB2 genes

BACKGROUND

Gene regulation is dynamic across cellular conditions and disease subtypes. From the aspect of regulation under modulation, regulation strength between a pair of genes can be modulated by (dependent on) expression abundance of another gene (modulator gene). Previous studies have demonstrated the involvement of genes modulated by single modulator genes in cancers, including breast cancer. However, analysis of multi-modulator co-modulation that can further delineate the landscape of complex gene regulation is, to our knowledge, unexplored previously. In the present study we aim to explore the joint effects of multiple modulator genes in modulating global gene regulation and dissect the biological functions in breast cancer.

RESULTS

To carry out the analysis, we proposed the Covariability-based Multiple Regression (CoMRe) method. The method is mainly built on a multiple regression model that takes expression levels of multiple modulators as inputs and regulation strength between genes as output. Pairs of genes were divided into groups based on their co-modulation patterns. Analyzing gene expression profiles from 286 breast cancer patients, CoMRe investigated ten candidate modulator genes that interacted and jointly determined global gene regulation. Among the candidate modulators, ESR1, ERBB2, and ADAM12 were found modulating the most numbers of gene pairs. The largest group of gene pairs was composed of ones that were modulated by merely ESR1. Functional annotation revealed that the group was significantly related to tumorigenesis and estrogen signaling in breast cancer. ESR1-ERBB2 co-modulation was the largest group modulated by more than one modulators. Similarly, the group was functionally associated with hormone stimulus, suggesting that functions of the two modulators are performed, at least partially, through modulation. The findings were validated in majorities of patients (> 99%) of two independent breast cancer datasets.

CONCLUSIONS

We have showed CoMRe is a robust method to discover critical modulators in gene regulatory networks, and it is capable of achieving reproducible and biologically meaningful results. Our data reveal that gene regulatory networks modulated by single modulator or co-modulated by multiple modulators play important roles in breast cancer. Findings of this report illuminate complex and dynamic gene regulation under modulation and its involvement in breast cancer.

The Influence of the CYP3A4*22 Polymorphism and CYP2D6 Polymorphisms on Serum Concentrations of Aripiprazole, Haloperidol, Pimozide, and Risperidone in Psychiatric Patients

INTRODUCTION

Cytochrome P450 3A4 (CYP3A4) is involved in the metabolism of greater than 50% of the prescribed drugs. Recently, the CYP3A4*22 allele was reported to be associated with lower CYP3A4 expression and activity. Quetiapine, an antipsychotic metabolized by only CYP3A4, displayed higher serum levels in CYP3A4*22 carriers. Aripiprazole, haloperidol, pimozide, and risperidone are antipsychotics that are metabolized by CYP3A4 and CYP2D6. We investigated to which degree the CYP3A4*22 single-nucleotide polymorphism affects serum concentrations of patients receiving these drugs and compared this with the influence of CYP2D6 polymorphisms.

METHODS

Eight hundred thirty-four adult patients were included in this study, of whom 130 used aripiprazole, 312 used haloperidol, 86 used pimozide, and 396 used risperidone. Serum levels of the drug and, if available, their active metabolites were collected as well as information on dose. Patients were genotyped for CYP3A4*22 using restriction fragment length polymorphism analysis. Genotyping for CYP2D6 was done with allele-specific polymerase chain reaction.

RESULTS

No differences were found in serum (dose-corrected) concentrations of the antipsychotics between CYP3A4*22 wild-type and carrier groups. In contrast, CYP2D6 genotype did affect dose-corrected concentrations of the antipsychotics: for example, median dose-corrected concentrations were 56%, 86%, and 400% higher in predicted poor metabolizers versus extensive metabolizers for aripiprazole (P = 0.004), haloperidol (P > 0.001), and risperidone (P < 0.001), respectively, although a multiple regression analysis showed that only 4% to 17% of the variation in these concentrations could be explained by CYP2D6 status.

CONCLUSIONS

Heterozygous presence of CYP3A4*22 does not increase serum levels of antipsychotics metabolized by both CYP3A4 and CYP2D6, whereas CYP2D6 polymorphisms do affect serum levels to a limited extent.

The Making of a CYP3A Biomarker Panel for Guiding Drug Therapy

CYP3A ranks among the most abundant cytochrome P450 enzymes in the liver, playing a dominant role in metabolic elimination of clinically used drugs. A main member in CYP3A family, CYP3A4 expression and activity vary considerably among individuals, attributable to genetic and non-genetic factors, affecting drug dosage and efficacy. However, the extent of genetic influence has remained unclear. This review assesses current knowledge on the genetic factors influencing CYP3A4 activity. Coding region CYP3A4 polymorphisms are rare and account for only a small portion of inter-person variability in CYP3A metabolism. Except for the promoter allele CYP3A4*1B with ambiguous effect on expression, common CYP3A4 regulatory polymorphisms were thought to be lacking. Recent studies have identified a relatively common regulatory polymorphism, designated CYP3A4*22 with robust effects on hepatic CYP3A4 expression. Combining CYP3A4*22 with CYP3A5 alleles *1, *3 and *7 has promise as a biomarker predicting overall CYP3A activity. Also contributing to variable expression, the role of polymorphisms in transcription factors and microRNAs is discussed.

Nuclear Receptors in Drug Metabolism, Drug Response and Drug Interactions

Orally delivered small-molecule therapeutics are metabolized in the liver and intestine by phase I and phase II drug-metabolizing enzymes (DMEs), and transport proteins coordinate drug influx (phase 0) and drug/drug-metabolite efflux (phase III). Genes involved in drug metabolism and disposition are induced by xenobiotic-activated nuclear receptors (NRs), i.e. PXR (pregnane X receptor) and CAR (constitutive androstane receptor), and by the 1α, 25-dihydroxy vitamin D3-activated vitamin D receptor (VDR), due to transactivation of xenobiotic-response elements (XREs) present in phase 0-III genes. Additional NRs, like HNF4-α, FXR, LXR-α play important roles in drug metabolism in certain settings, such as in relation to cholesterol and bile acid metabolism. The phase I enzymes CYP3A4/A5, CYP2D6, CYP2B6, CYP2C9, CYP2C19, CYP1A2, CYP2C8, CYP2A6, CYP2J2, and CYP2E1 metabolize >90% of all prescription drugs, and phase II conjugation of hydrophilic functional groups (with/without phase I modification) facilitates drug clearance. The conjugation step is mediated by broad-specificity transferases like UGTs, SULTs, GSTs. This review delves into our current understanding of PXR/CAR/VDR-mediated regulation of DME and transporter expression, as well as effects of single nucleotide polymorphism (SNP) and epigenome (specified by promoter methylation, histone modification, microRNAs, long non coding RNAs) on the expression of PXR/CAR/VDR and phase 0-III mediators, and their impacts on variable drug response. Therapeutic agents that target epigenetic regulation and the molecular basis and consequences (overdosing, underdosing, or beneficial outcome) of drug-drug/drug-food/drug-herb interactions are also discussed. Precision medicine requires understanding of a drug's impact on DME and transporter activity and their NR-regulated expression in order to achieve optimal drug efficacy without adverse drug reactions. In future drug screening, new tools such as humanized mouse models and microfluidic organs-on-chips, which mimic the physiology of a multicellular environment, will likely replace the current cell-based workflow.

Regulation of coagulation factor XI expression by microRNAs in the human liver

High levels of factor XI (FXI) increase the risk of thromboembolic disease. However, the genetic and environmental factors regulating FXI expression are still largely unknown. The aim of our study was to evaluate the regulation of FXI by microRNAs (miRNAs) in the human liver. In silico prediction yielded four miRNA candidates that might regulate FXI expression. HepG2 cells were transfected with miR-181a-5p, miR-23a-3p, miR-16-5p and miR-195-5p. We used mir-494, which was not predicted to bind to F11, as a negative control. Only miR-181a-5p caused a significant decrease both in FXI protein and F11 mRNA levels. In addition, transfection with a miR-181a-5p inhibitor in PLC/PRF/5 hepatic cells increased both the levels of F11 mRNA and extracellular FXI. Luciferase assays in human colon cancer cells deficient for Dicer (HCT-DK) demonstrated a direct interaction between miR-181a-5p and 3′untranslated region of F11. Additionally, F11 mRNA levels were inversely and significantly correlated with miR-181a-5p levels in 114 healthy livers, but not with miR-494. This study demonstrates that FXI expression is directly regulated by a specific miRNA, miR-181a-5p, in the human liver. Future studies are necessary to further investigate the potential consequences of miRNA dysregulation in pathologies involving FXI.

Pharmacogenomics of acetaminophen in pediatric populations: a moving target

Acetaminophen (APAP) is widely used as an over-the-counter fever reducer and pain reliever. However, the current therapeutic use of APAP is not optimal. The inter-patient variability in both efficacy and toxicity limits the use of this drug. This is particularly an issue in pediatric populations, where tools for predicting drug efficacy and developmental toxicity are not well established. Variability in toxicity between age groups may be accounted for by differences in metabolism, transport, and the genetics behind those differences. While pharmacogenomics has been revolutionizing the paradigm of pharmacotherapy for many drugs, its application in pediatric populations faces significant challenges given the dynamic ontogenic changes in cellular and systems physiology. In this review we focused on the ontogenesis of the regulatory pathways involved in the disposition of APAP and on the variability between pediatric, adolescent, and adult patients. We also summarize important polymorphisms of the pharmacogenes associated with APAP metabolism. Pharmacogenetic studies in pediatric APAP treatment are also reviewed. We conclude that while a consensus in pharmacogenetic management of APAP in pediatric populations has not been achieved, a systems biology based strategy for comprehensively understanding the ontogenic regulatory pathway as well as the interaction between age and genetic variations are particularly necessary in order to address this question.

Pregnancy and pharmacogenomics in the context of drug metabolism and response

It is well-known that profound physiological and biochemical changes occur throughout the course of pregnancy. At the same time, the role of pharmacogenomics in modulating the metabolism and response profile to numerous medications has been elucidated. Yet, the clinical impact of pharmacogenomics during pregnancy is less well understood. We present an overview of factors modulating the pharmacokinetics and pharmacodynamics of medications throughout the time span of pregnancy while providing insights on how pharmacogenomics may contribute to interindividual variability in drug metabolism and response amongst pregnant women.