Gene expression variability in human hepatic drug metabolizing enzymes and transporters

Innovations, Opportunities, and Challenges for Predicting Alteration in Drug-Metabolizing Enzyme and Transporter Activity in Specific Populations

Drug-metabolizing enzymes and transporters (DMETs) are key regulators of the pharmacokinetics, efficacy, and toxicity of therapeutics. Over the past two decades, significant advancements in in vitro methodologies, targeted proteomics, in vitro to in vivo extrapolation methods, and integrated computational approaches such as physiologically based pharmacokinetic modeling have unequivocally contributed to improving our ability to quantitatively predict the role of DMETs in absorption, distribution, metabolism, and excretion and drug-drug interactions. However, the paucity of data regarding alterations in DMET activity in specific populations such as pregnant individuals, lactation, pediatrics, geriatrics, organ impairment, and disease states such as, cancer, kidney, and liver diseases and inflammation has restricted our ability to realize the full potential of these recent advancements. We envision that a series of carefully curated articles in a special supplementary issue of Drug Metabolism and Disposition will summarize the latest progress in in silico, in vitro, and in vivo approaches to characterize alteration in DMET activity and quantitatively predict drug disposition in specific populations. In addition, the supplementary issue will underscore the current scientific knowledge gaps that present formidable barriers to fully understand the clinical implications of altered DMET activity in specific populations and highlight opportunities for multistakeholder collaboration to advance our collective understanding of this rapidly emerging area. SIGNIFICANCE STATEMENT: This commentary highlights current knowledge and identifies gaps and key challenges in understanding the role of drug-metabolizing enzymes and transporters (DMETs) in drug disposition in specific populations. With this commentary for the special issue in Drug Metabolism and Disposition, the authors intend to increase interest and invite potential contributors whose research is focused or has aided in expanding the understanding around the role and impact of DMETs in drug disposition in specific populations.

Drug transporter expression and activity in cryopreserved human hepatocytes isolated from chimeric TK-NOG mice with humanized livers

Chimeric mice with humanized liver are thought to represent a sustainable source of isolated human hepatocytes for in vitro studying detoxification of drugs in humans. Because drug transporters are now recognized as key-actors of the hepatic detoxifying process, the present study was designed to characterize mRNA expression and activity of main hepatic drug transporters in cryopreserved human hepatocytes isolated from chimeric TK-NOG mice and termed HepaSH cells. Such cells after thawing were shown to exhibit a profile of hepatic solute carrier (SLC) and ATP-binding cassette (ABC) drug transporter mRNA levels well correlated to those found in cryopreserved primary human hepatocytes or human livers. HepaSH cells used either as suspensions or as 24 h-cultures additionally displayed notable activities of uptake SLCs, including organic anion transporting polypeptides (OATPs), organic anion transporter 2 (OAT2) or sodium-taurocholate co-transporting polypeptide (NTCP). SLC transporter mRNA expression, as well as SLC activities, nevertheless fell in HepaSH cells cultured for 120 h, which may reflect a partial dedifferentiation of these cells with time in culture in the conventional monolayer culture conditions used in the study. These data therefore support the use of cryopreserved HepaSH cells as either suspensions or short-term cultures for drug transport studies.

Administration of Secretome Derived from Human Mesenchymal Stem Cells Induces Hepatoprotective Effects in Models of Idiosyncratic Drug-Induced Liver Injury Caused by Amiodarone or Tamoxifen

Drug-induced liver injury (DILI) is one of the leading causes of acute liver injury. While many factors may contribute to the susceptibility to DILI, obese patients with hepatic steatosis are particularly prone to suffer DILI. The secretome derived from mesenchymal stem cell has been shown to have hepatoprotective effects in diverse in vitro and in vivo models. In this study, we evaluate whether MSC secretome could improve DILI mediated by amiodarone (AMI) or tamoxifen (TMX). Hepatic HepG2 and HepaRG cells were incubated with AMI or TMX, alone or with the secretome of MSCs obtained from human adipose tissue. These studies demonstrate that coincubation of AMI or TMX with MSC secretome increases cell viability, prevents the activation of apoptosis pathways, and stimulates the expression of priming phase genes, leading to higher proliferation rates. As proof of concept, in a C57BL/6 mouse model of hepatic steatosis and chronic exposure to AMI, the MSC secretome was administered endovenously. In this study, liver injury was significantly attenuated, with a decrease in cell infiltration and stimulation of the regenerative response. The present results indicate that MSC secretome administration has the potential to be an adjunctive cell-free therapy to prevent liver failure derived from DILI caused by TMX or AMI.

Deciphering genetic causes for sex differences in human health through drug metabolism and transporter genes

Sex differences have been widely observed in human health. However, little is known about the underlying mechanism behind these observed sex differences. We hypothesize that sex-differentiated genetic effects are contributors of these phenotypic differences. Focusing on a collection of drug metabolism enzymes and transporters (DMET) genes, we discover sex-differentiated genetic regulatory mechanisms between these genes and human complex traits. Here, we show that sex-differentiated genetic effects were present at genome-level and at DMET gene regions for many human complex traits. These sex-differentiated regulatory mechanisms are reflected in the levels of gene expression and endogenous serum biomarkers. Through Mendelian Randomization analysis, we identify putative sex-differentiated causal effects in each sex separately. Furthermore, we identify and validate sex differential gene expression of a subset of DMET genes in human liver samples. We observe higher protein abundance and enzyme activity of CYP1A2 in male-derived liver microsomes, which leads to higher level of an active metabolite formation of clozapine, a commonly prescribed antipsychotic drug. Taken together, our results demonstrate the presence of sex-differentiated genetic effects on DMET gene regulation, which manifest in various phenotypic traits including disease risks and drug responses.

Hepatic Cytochrome P450 Abundance and Activity in the Developing and Adult Göttingen Minipig: Pivotal Data for PBPK Modeling

The Göttingen Minipig is gaining ground as nonrodent species in safety testing of drugs for pediatric indications. Due to developmental changes in pharmacokinetics and pharmacodynamics, physiologically based pharmacokinetic (PBPK) models are built to better predict drug exposure in children and to aid species selection for nonclinical safety studies. These PBPK models require high quality physiological and ADME data such as protein abundance of drug metabolizing enzymes. These data are available for man and rat, but scarce for the Göttingen Minipig. The aim of this study was to assess hepatic cytochrome P450 (CYP) protein abundance in the developing Göttingen Minipig by using mass spectrometry. In addition, sex-related differences in CYP protein abundance and correlation of CYP enzyme activity with CYP protein abundance were assessed. The following age groups were included: gestational day (GD) 84–86 (n = 8), GD 108 (n = 6), postnatal day (PND) 1 (n = 8), PND 3 (n = 8), PND 7 (n = 8), PND 28 (n = 8) and adult (n = 8). Liver microsomes were extracted and protein abundance was compared to that in adult animals. Next, the CYP protein abundance was correlated to CYP enzyme activity in the same biological samples. In general, CYP protein abundance gradually increased during development. However, we observed a stable protein expression over time for CYP4A24 and CYP20A1 and for CYP51A1, a high protein expression during the fetal stages was followed by a decrease during the first month of life and an increase toward adulthood. Sex-related differences were observed for CYP4V2_2a and CYP20A1 at PND 1 with highest expression in females for both isoforms. In the adult samples, sex-related differences were detected for CYP1A1, CYP1A2, CYP2A19, CYP2E1_2, CYP3A22, CYP4V2_2a and CYP4V2_2b with higher values in female compared to male Göttingen Minipigs. The correlation analysis between CYP protein abundance and CYP enzyme activity showed that CYP3A22 protein abundance correlated clearly with the metabolism of midazolam at PND 7. These data are remarkably comparable to human data and provide a valuable step forward in the construction of a neonatal and juvenile Göttingen Minipig PBPK model.

The 1β-Hydroxy-Deoxycholic Acid to Deoxycholic Acid Urinary Metabolic Ratio: Toward a Phenotyping of CYP3A Using an Endogenous Marker?

In this study, we assessed the potential use of the 1β-hydroxy-deoxycholic acid (1β-OH-DCA) to deoxycholic acid (DCA) urinary metabolic ratio (UMR) as a CYP3A metric in ten male healthy volunteers. Midazolam (MDZ) 1 mg was administered orally at three sessions: alone (control session), after pre-treatment with fluvoxamine 50 mg (12 h and 2 h prior to MDZ administration), and voriconazole 400 mg (2 h before MDZ administration) (inhibition session), and after a 7-day pre-treatment with the inducer rifampicin 600 mg (induction session). The 1β-OH-DCA/DCA UMR was measured at each session, and correlations with MDZ metrics were established. At baseline, the 1β-OH-DCA/DCA UMR correlated significantly with oral MDZ clearance (r = 0.652, p = 0.041) and C_max (r = -0.652, p = 0.041). In addition, the modulation of CYP3A was reflected in the 1β-OH-DCA/DCA UMR after the intake of rifampicin (induction ratio = 11.4, p < 0.01). During the inhibition session, a non-significant 22% decrease in 1β-OH-DCA/DCA was observed (p = 0.275). This result could be explained by the short duration of CYP3A inhibitors intake fixed in our clinical trial. Additional studies, particularly involving CYP3A inhibition for a longer period and larger sample sizes, are needed to confirm the 1β-OH-DCA/DCA metric as a suitable CYP3A biomarker.

Use of computational toxicology (CompTox) tools to predict in vivo toxicity for risk assessment

Computational Toxicology tools were used to predict toxicity for three pesticides: propyzamide (PZ), carbaryl (CB) and chlorpyrifos (CPF). The tools used included: a) ToxCast/Tox21 assays (AC₅₀ s μM: concentration 50% maximum activity); b) in vitro-to-in vivo extrapolation (IVIVE) using ToxCast/Tox21 AC₅₀s to predict administered equivalent doses (AED: mg/kg/d) to compare to known in vivo Lowest-Observed-Effect-Level (LOEL)/Benchmark Dose (BMD); c) high throughput toxicokinetics population based (HTTK-Pop) using AC₅₀s for endpoints associated with the mode of action (MOA) to predict age-adjusted AED for comparison with in vivo LOEL/BMDs. ToxCast/Tox21 active-hit-calls for each chemical were predictive of targets associated with each MOA, however, assays directly relevant to the MOAs for each chemical were limited. IVIVE AEDs were predictive of in vivo LOEL/BMD₁₀s for all three pesticides. HTTK-Pop was predictive of in vivo LOEL/BMD₁₀s for PZ and CPF but not for CB after human age adjustments 11-15 (PZ) and 6-10 (CB) or 6-10 and 11-20 (CPF) corresponding to treated rat ages (in vivo endpoints). The predictions of computational tools are useful for risk assessment to identify targets in chemical MOAs and to support in vivo endpoints. Data can also aid is decisions about the need for further studies.

Large-scale analysis of human gene expression variability associates highly variable drug targets with lower drug effectiveness and safety

Motivation

The effectiveness of drugs tends to vary between patients. One of the well-known reasons for this phenomenon is genetic polymorphisms in drug target genes among patients. Here, we propose that differences in expression levels of drug target genes across individuals can also contribute to this phenomenon.

Results

To explore this hypothesis, we analyzed the expression variability of protein-coding genes, and particularly drug target genes, across individuals. For this, we developed a novel variability measure, termed local coefficient of variation (LCV), which ranks the expression variability of each gene relative to genes with similar expression levels. Unlike commonly used methods, LCV neutralizes expression levels biases without imposing any distribution over the variation and is robust to data incompleteness. Application of LCV to RNA-sequencing profiles of 19 human tissues and to target genes of 1076 approved drugs revealed that drug target genes were significantly more variable than protein-coding genes. Analysis of 113 drugs with available effectiveness scores showed that drugs targeting highly variable genes tended to be less effective in the population. Furthermore, comparison of approved drugs to drugs that were withdrawn from the market showed that withdrawn drugs targeted significantly more variable genes than approved drugs. Last, upon analyzing gender differences we found that the variability of drug target genes was similar between men and women. Altogether, our results suggest that expression variability of drug target genes could contribute to the variable responsiveness and effectiveness of drugs, and is worth considering during drug treatment and development.

Availability and implementation

LCV is available as a python script in GitHub (https://github.com/eyalsim/LCV).

Supplementary information

Supplementary data are available at Bioinformatics online.

MicroRNA hsa-miR-1301-3p Regulates Human ADH6, ALDH5A1 and ALDH8A1 in the Ethanol-Acetaldehyde-Acetate Metabolic Pathway

Alcohol dehydrogenases (ADHs) and aldehyde dehydrogenases (ALDHs) are vital enzymes involved in the metabolism of a variety of alcohols. Differences in the expression and enzymatic activity of human ADHs and ALDHs correlate with individual variability in metabolizing alcohols and drugs and in the susceptibility to alcoholic liver disease. MicroRNAs (miRNAs) function as epigenetic modulators to regulate the expression of drug-metabolizing enzymes. To characterize miRNAs that target ADHs and ALDHs in human liver cells, we carried out a systematic bioinformatics analysis to analyze free energies of the interaction between miRNAs and their cognate sequences in ADH and ALDH transcripts and then calculated expression correlations between miRNAs and their targeting ADH and ALDH genes using a public data base. Candidate miRNAs were selected to evaluate bioinformatic predictions using a series of biochemical assays. Our results showed that 11 miRNAs have the potential to modulate the expression of two ADH and seven ALDH genes in the human liver. We found that hsa-miR-1301-3p suppressed the expression of ADH6, ALDH5A1, and ALDH8A1 in liver cells and blocked their induction by ethanol. In summary, our results revealed that hsa-miR-1301-3p plays an important role in ethanol metabolism by regulating ADH and ALDH gene expression. SIGNIFICANCE STATEMENT: Systematic bioinformatics analysis showed that 11 microRNAs might play regulatory roles in the expression of two alcohol dehydrogenase (ADH) and seven aldehyde dehydrogenase (ALDH) genes in the human liver. Experimental evidences proved that hsa-miR-1301-3p suppressed the expression of ADH6, ALDH5A1, and ALDH8A1 in liver cells and decreased their inducibility by ethanol.

Strain-Dependent Variability of Early Discovery Small Molecule Pharmacokinetics in Mice: Does Strain Matter?

Drug discovery programs routinely perform pharmacokinetic (PK) studies in mice to prioritize lead compounds based on anticipated exposure-efficacy and exposure-toxicity relationships. Because of logistical and/or technical issues, the strain of mouse in early discovery PK studies may not always match the strain in toxicity or efficacy studies. This elicits the question do appreciable strain-dependent differences in PK parameters exist to an extent that would warrant conducting PK studies in a strain that matches efficacy and toxicity models? To understand the impact that strain may have on PK parameters, we selected eight marketed drugs with well characterized absorption, distribution, metabolism, and excretion properties and diverse structures to perform PK studies in three common mouse strains (Bagg Albino c, C57BL/6, and CD-1). Some statistical strain-dependent differences were observed; however, we found good general agreement of PK parameters between strains: 88%, 100%, 75%, 76%, 94%, and 88% of compounds were within twofold across strains for clearance, volume of distribution at steady state, t _1/2, C _max, T _max, and oral bioavailability, respectively. Overall, we recommend that an approach using a single strain of mouse is appropriate for discovery screening PK studies, provided that proper caution is exercised. SIGNIFICANCE STATEMENT: The mouse strain in discovery pharmacokinetic (PK) studies may not match the strain in efficacy and toxicology studies. Currently, there is a gap in the literature addressing whether differences in PK parameters across mouse strains exist such that multiple PK studies are warranted. The results from this study indicated that the PK properties of clinically used drugs between mouse strains are within an acceptable range such that single strain PK is appropriate.

The role of hepatic cytochrome P450s in the cytotoxicity of sertraline

Sertraline, an antidepressant, is commonly used to manage mental health symptoms related to depression, anxiety disorders, and obsessive-compulsive disorder. The use of sertraline has been associated with rare but severe hepatotoxicity. Previous research demonstrated that mitochondrial dysfunction, apoptosis, and endoplasmic reticulum stress were involved in sertraline-associated cytotoxicity. In this study, we reported that after a 24-h treatment in HepG2 cells, sertraline caused cytotoxicity, suppressed topoisomerase I and IIα, and damaged DNA in a concentration-dependent manner. We also investigated the role of cytochrome P450 (CYP)-mediated metabolism in sertraline-induced toxicity using our previously established HepG2 cell lines individually expressing 14 CYPs (1A1, 1A2, 1B1, 2A6, 2B6, 2C8, 2C9, 2C18, 2C19, 2D6, 2E1, 3A4, 3A5, and 3A7). We demonstrated that CYP2D6, 2C19, 2B6, and 2C9 metabolize sertraline, and sertraline-induced cytotoxicity was significantly decreased in the cells expressing these CYPs. Western blot analysis demonstrated that the induction of ɣH2A.X (a hallmark of DNA damage) and topoisomerase inhibition were partially reversed in CYP2D6-, 2C19-, 2B6-, and 2C9-overexpressing HepG2 cells. These data indicate that DNA damage and topoisomerase inhibition are involved in sertraline-induced cytotoxicity and that CYPs-mediated metabolism plays a role in decreasing the toxicity of sertraline.

Long noncoding RNA LINC00844-mediated molecular network regulates expression of drug metabolizing enzymes and nuclear receptors in human liver cells

Noncoding RNAs, such as long noncoding RNAs (lncRNAs) and microRNAs (miRNAs), regulate gene expression in many physiological and pathological processes, including drug metabolism. Drug metabolizing enzymes (DMEs) are critical components in drug-induced liver toxicity. In this study, we used human hepatic HepaRG cells treated with 5 or 10 mM acetaminophen (APAP) as a model system and identified LINC00844 as a toxicity-responsive lncRNA. We analyzed the expression profiles of LINC00844 in different human tissues. In addition, we examined the correlations between the levels of LINC00844 and those of key DMEs and nuclear receptors (NRs) for APAP metabolism in humans. Our results showed that lncRNA LINC00844 is enriched in the liver and its expression correlates positively with mRNA levels of CYP3A4, CYP2E1, SULT2A1, pregnane X receptor (PXR), and hepatocyte nuclear factor (HNF) 4α. We demonstrated that LINC00844 regulates the expression of these five genes in HepaRG cells using gain- and loss-of-function assays. Further, we discovered that LINC00844 is localized predominantly in the cytoplasm and acts as an hsa-miR-486-5p sponge, via direct binding, to protect SULT2A1 from miRNA-mediated gene silencing. Our data also demonstrated a functional interaction between LINC00844 and hsa-miR-486-5p in regulating DME and NR expression in HepaRG cells and primary human hepatocytes. We depicted a LINC00844-mediated regulatory network that involves miRNA and NRs and influences DME expression in response to APAP toxicity.

Sequencing XMET genes to promote genotype-guided risk assessment and precision medicine

High-throughput next generation sequencing (NGS) is a shotgun approach applied in a parallel fashion by which the genome is fragmented and sequenced through small pieces and then analyzed either by aligning to a known reference genome or by de novo assembly without reference genome. This technology has led researchers to conduct an explosion of sequencing related projects in multidisciplinary fields of science. However, due to the limitations of sequencing-based chemistry, length of sequencing reads and the complexity of genes, it is difficult to determine the sequences of some portions of the human genome, leaving gaps in genomic data that frustrate further analysis. Particularly, some complex genes are difficult to be accurately sequenced or mapped because they contain high GC-content and/or low complexity regions, and complicated pseudogenes, such as the genes encoding xenobiotic metabolizing enzymes and transporters (XMETs). The genetic variants in XMET genes are critical to predicate inter-individual variability in drug efficacy, drug safety and susceptibility to environmental toxicity. We summarized and discussed challenges, wet-lab methods, and bioinformatics algorithms in sequencing "complex" XMET genes, which may provide insightful information in the application of NGS technology for implementation in toxicogenomics and pharmacogenomics.

The role of hepatic cytochrome P450s in the cytotoxicity of dronedarone

Dronedarone is used to treat patients with cardiac arrhythmias and has been reported to be associated with liver injury. Our previous mechanistic work demonstrated that DNA damage-induced apoptosis contributes to the cytotoxicity of dronedarone. In this study, we examined further the underlying mechanisms and found that after a 24-h treatment of HepG2 cells, dronedarone caused cytotoxicity, G1-phase cell cycle arrest, suppression of topoisomerase II, and DNA damage in a concentration-dependent manner. We also investigated the role of cytochrome P450s (CYPs)-mediated metabolism in the dronedarone-induced toxicity using our previously established HepG2 cell lines expressing individually 14 human CYPs (1A1, 1A2, 1B1, 2A6, 2B6, 2C8, 2C9, 2C18, 2C19, 2D6, 2E1, 3A4, 3A5, and 3A7). We demonstrated that CYP3A4, 3A5, and 2D6 were the major enzymes that metabolize dronedarone, and that CYP3A7, 2E1, 2C19, 2C18, 1A1, and 2B6 also metabolize dronedarone, but to a lesser extent. Our data showed that the cytotoxicity of dronedarone was decreased in CYP3A4-, 3A5-, or 2D6-overexpressing cells compared to the control HepG2 cells, indicating that the parent dronedarone has higher potency than the metabolites to induce cytotoxicity in these cells. In contrast, cytotoxicity was increased in CYP1A1-overexpressing cells, demonstrating that CYP1A1 exerts an opposite effect in dronedarone's toxicity, comparing to CYP3A4, 3A5, or 2D6. We also studied the involvement of topoisomerase II in dronedarone-induced toxicity, and demonstrated that the overexpression of topoisomerase II caused an increase in cell viability and a decrease in γ-H2A.X induction, suggesting that suppression of topoisomerase II may be one of the mechanisms involved in dronedarone-induced liver toxicity.

A functional SNP in the 3'-UTR of TAP2 gene interacts with microRNA hsa-miR-1270 to suppress the gene expression

The transporter associated with antigen processing 2 (TAP2) is involved in the development of multidrug resistance and the etiology of immunological diseases. In this study, we investigated whether the expression of TAP2 can be perturbed by single nucleotide polymorphisms (SNPs) located in 3'-untranslated region (3'-UTR) of the gene via interactions with microRNAs. Using a series of in silico assays, we selected the candidate microRNAs (miRNAs) with the potential to interact with functional SNPs of TAP2. The SNP rs241456-located in the 3'-UTR of TAP2-resides in a potential binding site for hsa-miR-1270 and hsa-miR-620. HEK 293 cells, from a human kidney cell line, were used to characterize the extent of binding of miRNAs to each polymorphic allele of the SNP by a luciferase reporter gene assay. RNA electrophoretic mobility shift assays were used to evaluate the interaction between the miRNAs and each allele sequence of the SNP. We found that hsa-miR-1270 inhibited luciferase activity by binding to the T allele of the SNP in an allele-specific manner. A negative correlation was also found between the expression of hsa-miR-1270 and the T allele of the SNP in kidney tissues. Our findings support the hypothesis that hsa-miR-1270 suppresses the production of TAP2 by binding to this SNP in the 3'-UTR of this gene. Environ. Mol. Mutagen. 59:134-143, 2018. © 2017 Wiley Periodicals, Inc.

Perspective on the Genetic Response to Antiparasitics: A Review Article

BACKGROUND

Drugs' pharmacokinetics and pharmacodynamics can be affected by diverse genetic variations, within which simple nucleotide polymorphisms (SNPs) are the most common. Genetic variability is one of the factors that could explain questions like why a given drug does not have the desired effect or why do adverse drug reactions arise.

METHODS

In this retrospective observational study, literature search limits were set within PubMed database as well as the epidemiological bulletins published by the Mexican Ministry of Health, from Jan 1st 2001 to Mar 31st 2017 (16 years).

RESULTS

Metabolism of antiparasitic drugs and their interindividual responses are mainly modified by variations in cytochrome P450 enzymes. These enzymes show high frequencies of polymorphic variability thus affecting the expression of CYP2C, CYP2A, CYP2A6, CYP2D6, CYP2E6 and CYP2A6 isoforms. Research in this field opens the door to new personalized treatment approaches in medicine.

CONCLUSION

Clinical and pharmacological utility yield by applying pharmacogenetics to antiparasitic treatments is not intended as a mean to improve the prescription process, but to select or exclude patients that could present adverse drug reactions as well as to evaluate genetic alterations which result in a diversity of responses, ultimately seeking to provide a more effective and safe treatment; therefore choosing a proper dose for the appropriate patient and the optimal treatment duration. Furthermore, pharmacogenetics assists in the development of vaccines. In other words, the aim of this discipline is to find therapeutic targets allowing personalized treatments.

Genetic and methylation variation in the CYP2B6 gene is related to circulating p,p'-dde levels in a population-based sample

OBJECTIVES

Since the metabolism of the organochlorine pesticide dichlorodiphenyltrichloroethane (DDT) is not fully known in humans, we evaluated if circulating levels of a major breakdown product of DDT, p,p'-DDE, were related to genome-wide genetic and methylation variation in a population-based sample.

METHODS

In the population-based Prospective Investigation of the Vasculature in Uppsala Seniors (PIVUS) study (1016 subjects all aged 70), circulating levels of p,p'-DDE were analyzed by high-resolution chromatography coupled to high-resolution mass spectrometry (HRGC/HRMS). Genetic variants were genotyped and imputed (1000 Genomes reference, March 2012 release). Methylation sites were assayed using the Illumina HumanMethylation450 array in whole blood. A genome-wide association study (GWAS) approach was applied.

RESULTS

Evidence for genome-wide significant association with p,p'-DDE levels was observed only for a locus at chromosome 19 corresponding to the CYP2B6 gene (lead SNP rs7260538). Subjects being homozygote for the G allele showed a median level of 472ng/g lipid, while the corresponding level for those being homozygote for the T allele was 192ng/g lipid (p=1.5×10^-31). An analysis conditioned on the lead SNP disclosed a distinct signal in the same gene (rs7255374, position chr19:41520351; p=2.2×10^-8). A whole-genome methylation analysis showed one significant relationship vs. p,p'-DDE levels (p=6.2×10^-9) located 7kb downstream the CYP2B6 gene (cg27089200, position chr19:41531976). This CpG-site was also related to the lead SNP (p=3.8×10^-35), but mediated only 4% of the effect of the lead SNP on p,p'-DDE levels.

CONCLUSION

Circulating levels of p,p'-DDE were related to genetic variation in the CYP2B6 gene in the general elderly population. DNA methylation in this gene is not closely linked to the p,p'-DDE levels.

Evaluation of a novel PXR-knockout in HepaRG™ cells

The nuclear pregnane X receptor (PXR) regulates the expression of genes involved in the metabolism, hepatobiliary disposition, and toxicity of drugs and endogenous compounds. PXR is a promiscuous nuclear hormone receptor (NHR) with significant ligand and DNA‐binding crosstalk with the constitutive androstane receptor (CAR); hence, defining the precise role of PXR in gene regulation is challenging. Here, utilising a novel PXR‐knockout (KO) HepaRG cell line, real‐time PCR analysis was conducted to determine PXR involvement for a range of inducers. The selective PXR agonist rifampicin, a selective CAR activator, 6‐(4‐chlorophenyl)imidazo[2,1‐b][1,3]thiazole‐5‐carbaldehyde O‐(3,4‐dichlorobenzyl)oxime (CITCO), and dual activators of CAR and PXR including phenobarbital (PB) were analyzed. HepaRG control cells (5F clone) were responsive to prototypical inducers of CYP2B6 and CYP3A4. No response was observed in the PXR‐KO cells treated with rifampicin. Induction of CYP3A4 by PB, artemisinin, and phenytoin was also much reduced in PXR‐KO cells, while the response to CITCO was maintained. This finding is in agreement with the abolition of functional PXR expression. The apparent EC₅₀ values for PB were in agreement between the cell lines; however, CITCO was ~threefold (0.3 μmol/L vs. 1 μmol/L) lower in the PXR‐KO cells compared with the 5F cells for CYP2B6 induction. Results presented support the application of the novel PXR‐KO cells in the definitive assignment of PXR‐mediated CYP2B6 and CYP3A4 induction. Utilization of such cell lines will allow advancement in composing structure activity relationships rather than relying predominantly on pharmacological manipulations and provide in‐depth mechanistic evaluation.

Minimum information required for a DMET experiment reporting

AIM

To provide pharmacogenomics reporting guidelines, the information and tools required for reporting to public omic databases.

MATERIAL & METHODS

For effective DMET data interpretation, sharing, interoperability, reproducibility and reporting, we propose the Minimum Information required for a DMET Experiment (MIDE) reporting.

RESULTS

MIDE provides reporting guidelines and describes the information required for reporting, data storage and data sharing in the form of XML.

CONCLUSION

The MIDE guidelines will benefit the scientific community with pharmacogenomics experiments, including reporting pharmacogenomics data from other technology platforms, with the tools that will ease and automate the generation of such reports using the standardized MIDE XML schema, facilitating the sharing, dissemination, reanalysis of datasets through accessible and transparent pharmacogenomics data reporting.

The regulation of human hepatic drug transporter expression by activation of xenobiotic-sensing nuclear receptors

INTRODUCTION

If a drug is found to be an inducer of hepatic drug metabolizing enzymes via activation of nuclear receptors such as pregnane X receptor (PXR) or constitutive androstane receptor (CAR), it is likely that drug transporters regulated through these same receptors will be induced as well. This review highlights what is currently known about the molecular mechanisms that regulate transporter expression and where the research is directed. Areas covered: This review is focused on publications that describe the role of activated hepatic nuclear receptors in the subsequent regulation of drug uptake and/or efflux transporters following exposure to xenobiotics. Expert opinion: Many of the published studies on the role of nuclear receptors in the regulation of drug transporters involve non-human test animals. But due to species response differences, these associations are not always applicable to humans. For this reason, some relevant human in vitro models have been developed, such as primary or cryopreserved human hepatocytes, human liver slices, or HepG2 or HuH7 cell lines transiently or stably transfected with PXR expression and reporter constructs as well as in vivo models such as PXR-humanized mice. These human-relevant test systems will continue to be developed and applied for the testing of investigational drugs.

MicroRNA hsa-miR-25-3p suppresses the expression and drug induction of CYP2B6 in human hepatocytes

Cytochrome P450 2B6 (CYP2B6), mainly expressed in the liver and brain, is important for processing a number of widely used drugs. Variations in CYP2B6 expression are associated with decreased drug efficacy or adverse effects in some patients. Although CYP2B6 genetic variants are associated with its differential expression, epigenetic mechanisms affecting CYP2B6 gene regulation have not been established. Sequence analysis identified 29 domains in the CYP2B6 mRNA transcript that could be subject to regulation by microRNAs. Inverse correlations were found in human hepatocytes for the levels of the microRNAs hsa-miR-504-5p and hsa-miR-25-3p compared with CYP2B6 mRNA. Reporter gene assays showed that hsa-miR-25-3p suppresses CYP2B6 expression by targeting a specific sequence in the 3'-untranslated region of the mRNA transcript. Electrophoretic mobility shift assays confirmed that hsa-miR-25-3p forms stable complexes with its cognate mRNA sequence and that it recruits cellular factors, including Ago-4. Transfection of HepaRG cells with hsa-miR-25-3p mimics inhibited expression of the endogenous CYP2B6 gene and it also decreased rifampicin-dependent induction of CYP2B6 at the mRNA and protein levels. In summary, in silico and in vitro analyses show that hsa-miR-25-3p suppresses CYP2B6 expression in human liver cells via an epigenetic mechanism.

The role of CYP 3A4 and 1A1 in amiodarone-induced hepatocellular toxicity

Amiodarone is a widely used potent antiarrhythmic for the treatment of cardiac disease; however, its use is often discontinued due to numerous adverse effects, including hepatotoxicity. To investigate the role of drug metabolism in this liver toxicity, amiodarone and its major metabolite desethylamiodarone were incubated with HepG2 cells overexpressing a series of cytochrome P450 (CYP) isoforms. Significantly higher cytotoxicity of amiodarone was observed in HepG2 cells overexpressing CYP3A4 or CYP1A1, compared with that observed in empty vector transduced control cells. Further, higher levels of the more potent hepatotoxic metabolite desethylamiodarone were detected in CYP3A4 or CYP1A1 expressed cells. The CYP3A4 inhibitor ketoconazole and the CYP1A1 inhibitor α-naphthoflavone drastically inhibited the metabolism of amiodarone to desethylamiodarone. Along with the inhibition of CYP1A1 or CYP3A4, the cytotoxicity of amiodarone was significantly reduced. These data indicate that the metabolism of amiodarone to desethylamiodarone by CYP1A1 or CYP3A4 plays an important role in the hepatocellular toxicity of amiodarone.

Challenges, Solutions, and Quality Metrics of Personal Genome Assembly in Advancing Precision Medicine

Even though each of us shares more than 99% of the DNA sequences in our genome, there are millions of sequence codes or structure in small regions that differ between individuals, giving us different characteristics of appearance or responsiveness to medical treatments. Currently, genetic variants in diseased tissues, such as tumors, are uncovered by exploring the differences between the reference genome and the sequences detected in the diseased tissue. However, the public reference genome was derived with the DNA from multiple individuals. As a result of this, the reference genome is incomplete and may misrepresent the sequence variants of the general population. The more reliable solution is to compare sequences of diseased tissue with its own genome sequence derived from tissue in a normal state. As the price to sequence the human genome has dropped dramatically to around $1000, it shows a promising future of documenting the personal genome for every individual. However, de novo assembly of individual genomes at an affordable cost is still challenging. Thus, till now, only a few human genomes have been fully assembled. In this review, we introduce the history of human genome sequencing and the evolution of sequencing platforms, from Sanger sequencing to emerging "third generation sequencing" technologies. We present the currently available de novo assembly and post-assembly software packages for human genome assembly and their requirements for computational infrastructures. We recommend that a combined hybrid assembly with long and short reads would be a promising way to generate good quality human genome assemblies and specify parameters for the quality assessment of assembly outcomes. We provide a perspective view of the benefit of using personal genomes as references and suggestions for obtaining a quality personal genome. Finally, we discuss the usage of the personal genome in aiding vaccine design and development, monitoring host immune-response, tailoring drug therapy and detecting tumors. We believe the precision medicine would largely benefit from bioinformatics solutions, particularly for personal genome assembly.

Clinical association between pharmacogenomics and adverse drug reactions

Adverse drug reactions (ADRs) are a major public health concern and cause significant patient morbidity and mortality. Pharmacogenomics is the study of how genetic polymorphisms affect an individual's response to pharmacotherapy at the level of a whole genome. This article updates our knowledge on how genetic polymorphisms of important genes alter the risk of ADR occurrence after an extensive literature search. To date, at least 244 pharmacogenes identified have been associated with ADRs of 176 clinically used drugs based on PharmGKB. At least 28 genes associated with the risk of ADRs have been listed by the Food and Drug Administration as pharmacogenomic biomarkers. With the availability of affordable and reliable testing tools, pharmacogenomics looks promising to predict, reduce, and minimize ADRs in selected populations.

Drug metabolism and clearance system in tumor cells of patients with multiple myeloma

Resistance to chemotherapy is a major limitation of cancer treatments with several molecular mechanisms involved, in particular altered local drug metabolism and detoxification process. The role of drug metabolism and clearance system has not been satisfactorily investigated in Multiple Myeloma (MM), a malignant plasma cell cancer for which a majority of patients escapes treatment. The expression of 350 genes encoding for uptake carriers, xenobiotic receptors, phase I and II Drug Metabolizing Enzymes (DMEs) and efflux transporters was interrogated in MM cells (MMCs) of newly-diagnosed patients in relation to their event free survival. MMCs of patients with a favourable outcome have an increased expression of genes coding for xenobiotic receptors (RXRα, LXR, CAR and FXR) and accordingly of their gene targets, influx transporters and phase I/II DMEs. On the contrary, MMCs of patients with unfavourable outcome displayed a global down regulation of genes coding for xenobiotic receptors and the downstream detoxification genes but had a high expression of genes coding for ARNT and Nrf2 pathways and ABC transporters. Altogether, these data suggests ARNT and Nrf2 pathways could be involved in MM primary resistance and that targeting RXRα, PXR, LXR and FXR through agonists could open new perspectives to alleviate or reverse MM drug resistance.

Clinically relevant genetic variants of drug-metabolizing enzyme and transporter genes detected in Thai children and adolescents with autism spectrum disorder

Single-nucleotide polymorphisms (SNPs) among drug-metabolizing enzymes and transporters (DMETs) influence the pharmacokinetic profile of drugs and exhibit intra- and interethnic variations in drug response in terms of efficacy and safety profile. The main objective of this study was to assess the frequency of allelic variants of drug absorption, distribution, metabolism, and elimination-related genes in Thai children and adolescents with autism spectrum disorder. Blood samples were drawn from 119 patients, and DNA was extracted. Genotyping was performed using the DMET Plus microarray platform. The allele frequencies of the DMET markers were generated using the DMET Console software. Thereafter, the genetic variations of significant DMET genes were assessed. The frequencies of SNPs across the genes coding for DMETs were determined. After filtering the SNPs, 489 of the 1,931 SNPs passed quality control. Many clinically relevant SNPs, including CYP2C19*2, CYP2D6*10, CYP3A5*3, and SLCO1B1*5, were found to have frequencies similar to those in the Chinese population. These data are important for further research to investigate the interpatient variability in pharmacokinetics and pharmacodynamics of drugs in clinical practice.

microRNAs as pharmacogenomic biomarkers for drug efficacy and drug safety assessment

Much evidence has documented that microRNAs (miRNAs) play an important role in the modulation of interindividual variability in the production of drug metabolizing enzymes and transporters (DMETs) and nuclear receptors (NRs) through multidirectional interactions involving environmental stimuli/stressors, the expression of miRNA molecules and genetic polymorphisms. MiRNA expression has been reported to be affected by drugs and miRNAs themselves may affect drug metabolism and toxicity. In cancer research, miRNA biomarkers have been identified to mediate intrinsic and acquired resistance to cancer therapies. In drug safety assessment, miRNAs have been found associated with cardiotoxicity, hepatotoxicity and nephrotoxicity. This review article summarizes published studies to show that miRNAs can serve as early biomarkers for the evaluation of drug efficacy and drug safety.

Development of HepG2-derived cells expressing cytochrome P450s for assessing metabolism-associated drug-induced liver toxicity

The generation of reactive metabolites from therapeutic agents is one of the major mechanisms of drug-induced liver injury (DILI). In order to evaluate metabolism-related toxicity and improve drug efficacy and safety, we generated a battery of HepG2-derived cell lines that express 14 cytochrome P450s (CYPs) (1A1, 1A2, 1B1, 2A6, 2B6, 2C8, 2C9, 2C18, 2C19, 2D6, 2E1, 3A4, 3A5 and 3A7) individually using a lentiviral expression system. The expression/production of a specific CYP in each cell line was confirmed by an increased abundance of the CYP at both mRNA and protein levels. Moreover, the enzymatic activities of representative CYPs in the corresponding cell lines were also measured. Using our CYP-expressed HepG2 cells, the toxicity of three drugs that could induce DILI (amiodarone, chlorpromazine and primaquine) was assessed, and all of them showed altered (increased or decreased) toxicity compared to the toxicity in drug-treated wild-type HepG2 cells. CYP-mediated drug toxicity examined in our cell system is consistent with previous reports, demonstrating the potential of these cells for assessing metabolism-related drug toxicity. This cell system provides a practical in vitro approach for drug metabolism screening and for early detection of drug toxicity. It is also a surrogate enzyme source for the enzymatic characterization of a particular CYP that contributes to drug-induced liver toxicity.

Modulation of ALDH5A1 and SLC22A7 by microRNA hsa-miR-29a-3p in human liver cells

Observed variations in drug responses among patients may result from differences in heritable genetic traits or from alterations in the epigenetic regulation of drug metabolizing enzymes and transporters (DMETs). MicroRNAs (miRNAs), a group of small non-coding RNAs, provide an epigenetic mechanism for fine-tuning the expression of targeted DMET genes by regulating the efficiency of protein translation and by decreasing mRNA stability via enhanced degradation. In the current study we systematically screened 374 important genes encoding DMETs for potential response elements to hsa-miR-29a-3p, a highly abundant miRNA in human liver. RNA electrophoresis mobility shift assays displayed direct interactions between hsa-miR-29a-3p and its cognate targets within the mRNA transcripts for the ABCC6, SLC22A7 and ALDH5A1 genes. The expression of luciferase reporter genes containing the 3'-UTRs of SLC22A7 or ALDH5A1 and the expression of endogenous SLC22A7 and ALDH5A1 were each suppressed by transfection with hsa-miR-29a-3p mimics. Importantly, chemically-induced up-regulation of hsa-miR-29a-3p correlated inversely with the expression of SLC22A7 and ALDH5A1. However, our studies failed to detect suppressive effects of hsa-miR-29a-3p on ABCC6 expression, which might be explained by the notion that the interaction of hsa-miR-29a-3p and ABCC6 mRNA was unable to recruit ribonucleoproteins to form a RNA-induced silencing complex.

Source-Related Effects of Wastewater on Transcription Factor (AhR, CAR and PXR)-Mediated Induction of Gene Expression in Cultured Rat Hepatocytes and Their Association with the Prevalence of Antimicrobial-Resistant Escherichia coli

Extracts of wastewater collected from 4 sewage treatment plants (STPs) receiving effluents from different sources in South India were investigated for their levels of transcription factor-mediated gene induction in primary cultured rat hepatocytes. In addition, the relation between gene induction levels and the prevalence of antimicrobial-resistant Escherichia coli (E. coli) in wastewater was examined. STP-3, which treats only hospital wastewater, exhibited significantly greater induction potency of all 6 drug metabolizing cytochrome P450 (CYP) genes examined, CYP1A1, 1A2, 1B1, 2B15, 3A1, and 3A2, whereas the wastewater at STP-1, which exclusively receives domestic sewage, showed significantly diminished levels of induction of 3 CYP genes when compared to the levels of CYP induction at STP-2, which receives mixed wastewater. Samples collected during the monsoon season showed a significantly altered gene induction capacity compared to that of samples from the pre-monsoon period. The data suggest that the toxicity of wastewater in STPs was not significantly diminished during the treatment process. The chemical-gene interaction data predicted that a vast number of chemicals present in the wastewater would stimulate the genes studied in the rat hepatocytes. The multivariable logistic regression analysis demonstrated that the prevalence of isolates resistant to cefotaxime, imipenem and streptomycin was significantly correlated with the levels of induction of at least three CYP-isozymes in STP wastewater. In addition, the resistance of isolates in treatment plants was not altered by the treatment steps, whereas the sampling season did have an impact on the resistance to specific antimicrobials. The identification of receptor-mediated gene regulation capacities offers important data not limited to the (synergistic) physiological role of chemicals in biological systems but may provide new insight into the link between the effects of known/unknown drugs and prevalence of antimicrobial-resistant bacteria in wastewater.

MicroRNA hsa-miR-29a-3p modulates CYP2C19 in human liver cells

Cytochrome P450 2C19 (CYP2C19) is involved in the metabolism of many drugs. Extensive studies have demonstrated that genetic variants and endogenous and environmental factors play important roles in the expression of CYP2C19. However, the role of microRNAs (miRNAs) in controlling CYP2C19 expression has not been investigated completely. In the present study, we performed in silico analysis to rank putative miRNA/CYP2C19 hybrids with regards to the predicted stabilities of their duplexes and then we applied a series of biochemical and molecular assays to elucidate the underlying functional mechanisms for the regulation of CYP2C19 by miRNAs. In silico analysis indicated that hsa-miR-23a-3p and hsa-miR-29a-3p target the coding region of CYP2C19 with hybrid stabilities of -27.5kcal/mol and -23.3kcal/mol, respectively. RNA electrophoresis mobility shift assays showed that both hsa-miR-23a-3p and hsa-miR-29a-3p miRNAs were able to bind directly to their cognate targets in the CYP2C19 transcript. Further, a significant inverse correlation was found between chemically-induced up-regulation of hsa-miR-29a-3p and CYP2C19 expression in HepaRG cells. In addition, inverse correlations were also observed in human liver tissue samples between the level of CYP2C19 mRNA expression and both hsa-miR-23a-3p and hsa-miR-29a-3p levels. All these results demonstrated the suppressing role of hsa-miR-29a-3p on CYP2C19 expression.

Two persistent organic pollutants which act through different xenosensors (alpha-endosulfan and 2,3,7,8 tetrachlorodibenzo-p-dioxin) interact in a mixture and downregulate multiple genes involved in human hepatocyte lipid and glucose metabolism

Individuals, typically, are exposed to mixtures of environmental xenobiotics affecting multiple organs and acting through different xenosensors and pathways in species and cell-type specific manners. 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) and α-endosulfan are Persistent Organic Pollutants (POPs) and endocrine disruptors which act through different xenosensors and accumulate in the liver. Our objective in this HEALS study was to investigate the effects of the mixture of these POPs on gene expression in a human-derived hepatocyte cell line, HepaRG. We found that, in spite of having largely uncorrelated effects, TCDD and α-endosulfan, when mixed, alter the expression of genes. The combined effects of the mixture of the POPs significantly altered the expression of 100 genes (42 up- and 58 down-regulated) whereas the same concentration of either POP alone did not alter significantly the expression of these genes. For 32 other genes, selective inhibitory crosstalk between TCDD and α-endosulfan was observed. One of the POPs inhibited the effect, on gene expression, of the other in the mixture although, when used alone, that POP did not affect expression. The expression of another 82 genes was significantly altered (up- or down-regulated) by a single POP. The addition of the second POP either increased, in the same direction, the effect on gene expression or had no further effect. At low concentrations (0.2 nM TCDD and 1 μM α-endosulfan), the POPs still had significant effects and the levels of expression of the corresponding proteins were found to be affected for some genes. Particularly striking was the 80-90% inhibition, by the mixture, of the expression of a number of genes of several hepatic intermediary metabolic pathways (glycerolipid metabolism, FXR/RXR activation, glycolysis/gluconeogenesis, retinoid and bile acid biosynthesis), whereas each pollutant alone had only a moderate effect.

Expression of drug transporters in human kidney: impact of sex, age, and ethnicity

Background

Differences in expression of drug transporters in human kidney contribute to changes in pharmacokinetics and toxicokinetics of a variety of drug compounds. The basal expression levels of genes involved in drug transport processes in the kidney introduces differences in bioavailability, distribution, and clearance of drugs, possibly influencing drug efficacy and adverse reactions. Sex differences in gene expression of transporters are a key cause of differences in sex-dependent pharmacokinetics, which may characterize many drugs and contribute to individual differences in drug efficacy and toxicity. Therefore, evaluating the expression of drug transporters in normal human kidneys is important to better understand differences in drug bioavailability, distribution, and clearance of drugs in humans. Other factors such as age and ethnicity may also contribute to individual differences in gene expression of drug transporters in the human kidney.

Methods

Quantitative real-time PCR (QRT-PCR) was performed to determine the gene expression of 30 drug transporters in 95 age-matched normal human kidney tissues. Multiple Student’s t-tests (Sidak-Bonferroni correction) and two-way ANOVA (Bonferroni correction) analyses were used to determine statistically significant differences.

Results

In the 30 transporter genes examined, sex, ethnicity, and age differences in gene expression were exhibited in normal human kidney tissue. These changes in expression were not found to be differentially significant. However, sex-age and sex-ethnicity interactions were found to be statistically significant. For sex-age interactions, SCL22A12 was found to be significantly higher expressed in females <50 years compared to males <50 years. Expression levels of SLC22A2, SLC22A12, SLC6A16, and ABCB6 were significantly higher in females <50 years compared to females ≥50 years. In sex-ethnicity interactions, expression levels of ATP7B and KCNJ8 were found to be significantly higher in African American females compared to European American females. Also, the expression of SLC31A2 was significantly higher in European American males compared to European American females.

Conclusions

Sex, age, and ethnic differences impacted the expression of drug transporters in normal human kidneys, which suggests that the analysis of gene expression of drug transporters will aid in improving the usage/dosage of drug therapies influencing personalized medicine and susceptibility to adverse drug reactions.

Suppression of CYP2C9 by microRNA hsa-miR-128-3p in human liver cells and association with hepatocellular carcinoma

Published studies have identified genetic variants, somatic mutations, and changes in gene expression profiles that are associated with hepatocellular carcinoma (HCC), particularly involving genes that encode drug metabolizing enzymes (DMEs). CYP2C9, one of the most abundant and important DMEs, is involved in the metabolism of many carcinogens and drugs and is down-regulated in HCC. To investigate the molecular mechanisms that control CYP2C9 expression, we applied integrative approaches including in silico, in vitro, and in vivo analyses to elucidate the role of microRNA hsa-miR-128-3p in the regulation of CYP2C9 expression and translation. RNA electrophoresis mobility shift assays demonstrated a direct interaction between hsa-miR-128-3p and its cognate target, the CYP2C9 transcript. Furthermore, the expression of a luciferase reporter gene containing the 3′-UTR of CYP2C9 and the endogenous expression of CYP2C9 were suppressed by transfection of hsa-miR-128-3p. Importantly, chemically-induced up- or down-regulation of hsa-miR-128-3p correlated inversely with the expression of CYP2C9. Finally, an association analysis revealed that the expression of hsa-miR-128-3p is inversely correlated with the expression of CYP2C9 in HCC tumor tissues. Altogether, the study helped to elucidate the mechanism of CYP2C9 regulation by hsa-miR-128-3p, and the inverse association in HCC.

Mechanistic basis of altered morphine disposition in nonalcoholic steatohepatitis

Morphine is metabolized in humans to morphine-3-glucuronide (M3G) and the pharmacologically active morphine-6-glucuronide (M6G). The hepatobiliary disposition of both metabolites relies upon multidrug resistance-associated proteins Mrp3 and Mrp2, located on the sinusoidal and canalicular membrane, respectively. Nonalcoholic steatohepatitis (NASH), the severe stage of nonalcoholic fatty liver disease, alters xenobiotic metabolizing enzyme and transporter function. The purpose of this study was to determine whether NASH contributes to the large interindividual variability and postoperative adverse events associated with morphine therapy. Male Sprague-Dawley rats were fed a control diet or a methionine- and choline-deficient diet to induce NASH. Radiolabeled morphine (2.5 mg/kg, 30 µCi/kg) was administered intravenously, and plasma and bile (0-150 or 0-240 minutes), liver and kidney, and cumulative urine were analyzed for morphine and M3G. The antinociceptive response to M6G (5 mg/kg) was assessed (0-12 hours) after direct intraperitoneal administration since rats do not produce M6G. NASH caused a net decrease in morphine concentrations in the bile and plasma and a net increase in the M3G/morphine plasma area under the concentration-time curve ratio, consistent with upregulation of UDP-glucuronosyltransferase Ugt2b1. Despite increased systemic exposure to M3G, NASH resulted in decreased biliary excretion and hepatic accumulation of M3G. This shift toward systemic retention is consistent with the mislocalization of canalicular Mrp2 and increased expression of sinusoidal Mrp3 in NASH and may correlate to increased antinociception by M6G. Increased metabolism and altered transporter regulation in NASH provide a mechanistic basis for interindividual variability in morphine disposition that may lead to opioid-related toxicity.

Cytochrome P450-mediated bioactivation of mefenamic acid to quinoneimine intermediates and inactivation by human glutathione S-transferases

Mefenamic acid (MFA) has been associated with rare but severe cases of hepatotoxicity, nephrotoxicity, gastrointestinal toxicity, and hypersensitivity reactions that are believed to result from the formation of reactive metabolites. Although formation of protein-reactive acylating metabolites by phase II metabolism has been well-studied and proposed to be the cause of these toxic side effects, the oxidative bioactivation of MFA has not yet been competely characterized. In the present study, the oxidative bioactivation of MFA was studied using human liver microsomes (HLM) and recombinant human P450 enzymes. In addition to the major metabolite 3'-OH-methyl-MFA, resulting from the benzylic hydroxylation by CYP2C9, 4'-hydroxy-MFA and 5-hydroxy-MFA were identified as metabolites resulting from oxidative metabolism of both aromatic rings of MFA. In the presence of GSH, three GSH conjugates were formed that appeared to result from GSH conjugation of the two quinoneimines formed by further oxidation of 4'-hydroxy-MFA and 5-hydroxy-MFA. The major GSH conjugate was identified as 4'-OH-5'-glutathionyl-MFA and was formed at the highest activity by CYP1A2 and to a lesser extent by CYP2C9 and CYP3A4. Two minor GSH conjugates resulted from secondary oxidation of 5-hydroxy-MFA and were formed at the highest activity by CYP1A2 and to a lesser extent by CYP3A4. Additionally, the ability of seven human glutathione S-transferases (hGSTs) to catalyze the GSH conjugation of the quinoneimines formed by P450s was also investigated. The highest increase of total GSH conjugation was observed with hGSTP1-1, followed by hepatic hGSTs hGSTA2-2 and hGSTM1-1. The results of this study show that, next to phase II metabolites, reactive quinoneimines formed by oxidative bioactivation might also contribute to the idiosyncratic toxicity of MFA.

DDI-CPI, a server that predicts drug-drug interactions through implementing the chemical-protein interactome

Drug–drug interactions (DDIs) may cause serious side-effects that draw great attention from both academia and industry. Since some DDIs are mediated by unexpected drug–human protein interactions, it is reasonable to analyze the chemical–protein interactome (CPI) profiles of the drugs to predict their DDIs. Here we introduce the DDI-CPI server, which can make real-time DDI predictions based only on molecular structure. When the user submits a molecule, the server will dock user's molecule across 611 human proteins, generating a CPI profile that can be used as a feature vector for the pre-constructed prediction model. It can suggest potential DDIs between the user's molecule and our library of 2515 drug molecules. In cross-validation and independent validation, the server achieved an AUC greater than 0.85. Additionally, by investigating the CPI profiles of predicted DDI, users can explore the PK/PD proteins that might be involved in a particular DDI. A 3D visualization of the drug-protein interaction will be provided as well. The DDI-CPI is freely accessible at http://cpi.bio-x.cn/ddi/.

Female gender as a susceptibility factor for drug-induced liver injury

Adverse drug reactions (ADRs) can involve all tissues and organs, but liver injuries are considered among the most serious. A number of prospective, multicenter studies have confirmed a higher risk of ADRs in general among female subjects compared to a male cohort. Although drug-induced liver injury (DILI) is infrequently encountered, the preponderance of evidence suggests that women appear to be more susceptible than men to fulminate hepatic/acute liver failure especially in response to some anti-infective drugs and to autoimmune-like hepatitis following exposure to certain other therapeutic drugs. A number of hypotheses have been proposed to explain this sex difference in susceptibility to DILI. Collectively, these hypotheses suggest three basic sex-dependent mechanisms that include differences in various aspects of drug pharmacokinetics (PK) or pharmacodynamics following the administration of certain drugs; specific hormonal effects or interactions with immunomodulating agents or signaling molecules; and differences in the adverse response of the immune system to some drugs, reactive drug metabolites, or drug-protein adducts. At the preclinical drug safety stage, there is a need for more research on hormonal effects on drug PK and for additional research on gender differences in aberrant immune responses that may lead to idiosyncratic DILI in some female patients. Because the detection of rare but serious hepatic ADRs requires the exposure of very large patient populations, pharmacovigilance networks will continue to play a key role in the postmarketing surveillance for their detection and reporting.