Cytochrome P450 2B6 (CYP2B6) belongs to the minor drug metabolizing P450s in human liver. Expression is highly variable both between individuals and within individuals, owing to non-genetic factors, genetic polymorphisms, inducibility, and irreversible inhibition by many compounds. Drugs metabolized mainly by CYP2B6 include artemisinin, bupropion, cyclophosphamide, efavirenz, ketamine, and methadone. CYP2B6 is one of the most polymorphic CYP genes in humans and variants have been shown to affect transcriptional regulation, splicing, mRNA and protein expression, and catalytic activity. Some variants appear to affect several functional levels simultaneously, thus, combined in haplotypes, leading to complex interactions between substrate-dependent and -independent mechanisms. The most common functionally deficient allele is CYP2B6*6 [Q172H, K262R], which occurs at frequencies of 15 to over 60% in different populations. The allele leads to lower expression in liver due to erroneous splicing. Recent investigations suggest that the amino acid changes contribute complex substrate-dependent effects at the activity level, although data from recombinant systems used by different researchers are not well in agreement with each other. Another important variant, CYP2B6*18 [I328T], occurs predominantly in Africans (4–12%) and does not express functional protein. A large number of uncharacterized variants are currently emerging from different ethnicities in the course of the 1000 Genomes Project. The CYP2B6 polymorphism is clinically relevant for HIV-infected patients treated with the reverse transcriptase inhibitor efavirenz, but it is increasingly being recognized for other drug substrates. This review summarizes recent advances on the functional and clinical significance of CYP2B6 and its genetic polymorphism, with particular emphasis on the comparison of kinetic data obtained with different substrates for variants expressed in different recombinant expression systems.

CYP3A4 is the most important drug metabolizing enzyme in adult humans because of its prominent expression in liver and gut and because of its broad substrate specificity, which includes drugs from most therapeutic categories and many endogenous substances. Expression and function of CYP3A4 vary extensively both intra- and interindividually thus contributing to unpredictable drug response and toxicity. A multitude of environmental, genetic, and physiological factors are known to influence CYP3A4 expression and activity. Among the best predictable sources of variation are drug–drug interactions, which are either caused by pregnane X-receptor (PXR), constitutive androstane receptor (CAR) mediated gene induction, or by inhibition through coadministered drugs or other chemicals, including also plant and food ingredients. Among physiological and pathophysiological factors are hormonal status, age, and gender, the latter of which was shown to result in higher levels in females compared to males, as well as inflammatory processes that downregulate CYP3A4 transcription. Despite the influence of these non-genetic factors, the genetic influence on CYP3A4 activity was estimated in previous twin studies and using information on repeated drug administration to account for 66% up to 88% of the interindividual variation. Although many single nucleotide polymorphisms (SNPs) within the CYP3A locus have been identified, genetic association studies have so far failed to explain a major part of the phenotypic variability. The term “missing heritability” has been used to denominate the gap between expected and known genetic contribution, e.g., for complex diseases, and is also used here in analogy. In this review we summarize CYP3A4 pharmacogenetics/genomics from the early inheritance estimations up to the most recent genetic and clinical studies, including new findings about SNPs in CYP3A4 (*22) and other genes (P450 oxidoreductase (POR), peroxisome proliferator-activated receptor alpha (PPARA)) with possible contribution to CYP3A4 variable expression.

Non-alcoholic fatty liver disease (NAFLD), the hepatic manifestation of the metabolic syndrome, is a complex multifactorial disease characterized by metabolic deregulations that include accumulation of lipids in the liver, lipotoxicity, and insulin resistance. The progression of NAFLD to non-alcoholic steatohepatitis and cirrhosis, and ultimately to carcinomas, is governed by interplay of pro-inflammatory pathways, oxidative stress, as well as fibrogenic and apoptotic cues. As the liver is the major organ of biotransformation, deregulations in hepatic signaling pathways have effects on both, xenobiotic and endobiotic metabolism. Several major nuclear receptors involved in the transcription and regulation of phase I and II drug metabolizing enzymes and transporters also have endobiotic ligands including several lipids. Hence, hepatic lipid accumulation in steatosis and NAFLD, which leads to deregulated activation patterns of nuclear receptors, may result in altered drug metabolism capacity in NAFLD patients. On the other hand, genetic and association studies have indicated that a malfunction in drug metabolism can affect the prevalence and severity of NAFLD. This review focuses on the complex interplay between NAFLD pathogenesis and drug metabolism. A better understanding of these relationships is a prerequisite for developing improved drug dosing algorithms for the pharmacotherapy of patients with different stages of NAFLD.

5-Fluorouracil (5-FU) is a commonly used drug for the treatment of malignant cancers. However, approximately 80% of patients undergoing 5-FU treatment suffer from gastrointestinal mucositis. The aim of this report was to identify the drug target for the 5-FU-induced intestinal mucositis. 5-FU-induced intestinal mucositis was established by intraperitoneally administering mice with 100 mg/kg 5-FU. Network analysis of gene expression profile and bioluminescent imaging were applied to identify the critical molecule associated with 5-FU-induced mucositis. Our data showed that 5-FU induced inflammation in the small intestine, characterized by the increased intestinal wall thickness and crypt length, the decreased villus height, and the increased myeloperoxidase activity in tissues and proinflammatory cytokine production in sera. Network analysis of 5-FU-affected genes by transcriptomic tool showed that the expression of genes was regulated by nuclear factor-κB (NF-κB), and NF-κB was the central molecule in the 5-FU-regulated biological network. NF-κB activity was activated by 5-FU in the intestine, which was judged by in vivo bioluminescence imaging and immunohistochemical staining. However, 5-aminosalicylic acid (5-ASA) inhibited 5-FU-induced NF-κB activation and proinflammatory cytokine production. Moreover, 5-FU-induced histological changes were improved by 5-ASA. In conclusion, our findings suggested that NF-κB was the critical molecule associated with the pathogenesis of 5-FU-induced mucositis, and inhibition of NF-κB activity ameliorated the mucosal damage caused by 5-FU.

Background

Organic cation transporters (OCTs) determine not only physiological processes but are also involved in the cellular uptake of anticancer agents. Based on microarray analyses in hepatocellular carcinoma (HCC), SLC22A1/OCT1 mRNA seems to be downregulated, but systematic protein expression data are currently missing. Moreover, the underlying molecular mechanisms responsible for altered SLC22A1 expression in HCC are not fully understood. Therefore, we investigated the role of DNA methylation in the transcriptional regulation of the family members SLC22A1/OCT1, SLC22A2/OCT2 and SLC22A3/OCT3 in HCC.

Methods

Semiquantitative immunohistochemistry of SLC22A1 protein expression was performed in paired HCC and histological normal adjacent liver tissues (n = 71) using tissue microarray analyses, and the results were correlated with clinicopathological features. DNA methylation, quantified by MALDI-TOF mass spectrometry and gene expression of SLC22A1, SLC22A2 and SLC22A3 were investigated using fresh-frozen HCC (n = 22) and non-tumor adjacent liver tissues as well as histologically normal liver samples (n = 120) from a large-scale liverbank.

Results

Based on tissue microarray analyses, we observed a significant downregulation of SLC22A1 protein expression in HCC compared to normal adjacent tissue (P < 0.0001). SLC22A1 expression was significantly inverse correlated with expression of the proliferation marker MIB1/Ki-67 (rs = -0.464, P < 0.0001). DNA methylation of SLC22A1 was significantly higher in HCC compared with non-tumor adjacent liver tissue and was lowest in histologically normal liver tissue. Methylation levels for SLC22A1 in combination with RASSF1A resulted in a specificity of > 90% and a sensitivity of 82% for discriminating HCC and tumor-free liver tissue.

Conclusions

DNA methylation of SLC22A1 is associated with downregulation of SLC22A1 in HCC and might be a new biomarker for HCC diagnosis and prognosis. Moreover, targeting SLC22A1 methylation by demethylating agents may offer a novel strategy for anticancer therapy of HCC.

CYP2D6 belongs to the cytochrome P450 superfamily of enzymes and plays an important role in the metabolism of 20–25% of clinically used drugs including antidepressants. It displays inter-individual and inter-ethnic variability in activity ranging from complete absence to excessive activity which causes adverse drug reactions and toxicity or therapy failure even at normal drug doses. This variability is due to genetic polymorphisms which form poor, intermediate, extensive or ultrarapid metaboliser phenotypes. This study aimed to determine CYP2D6 alleles and their frequencies in the United Arab Emirates (UAE) local population. CYP2D6 alleles and genotypes were determined by direct DNA sequencing in 151 Emiratis with the majority being psychiatric patients on antidepressants. Several new alleles have been identified and in total we identified seventeen alleles and 49 genotypes. CYP2D6*1 (wild type) and CYP2D6*2 alleles (extensive metaboliser phenotype) were found with frequencies of 39.1% and 12.2%, respectively. CYP2D6*41 (intermediate metaboliser) occurred in 15.2%. Homozygous CYP2D6*4 allele (poor metaboliser) was found with a frequency of 2% while homozygous and heterozygous CYP2D6*4 occurred with a frequency of 9%. CYP2D6*2xn, caused by gene duplication (ultrarapid metaboliser) had a frequency of 4.3%. CYP2D6 gene duplication/multiduplication occurred in 16% but only 11.2% who carried more than 2 active functional alleles were considered ultrarapid metabolisers. CYP2D6 gene deletion in one copy occurred in 7.5% of the study group. In conclusion, CYP2D6 gene locus is heterogeneous in the UAE national population and no significant differences have been identified between the psychiatric patients and controls.

Thiopurine methyltransferase (Tpmt) is the primary enzyme responsible for deactivating thiopurine drugs. Thiopurine drugs (i.e., thioguanine [TG], mercaptopurine, azathioprine) are commonly used for the treatment of cancer, organ transplant, and autoimmune disorders. Chronic thiopurine therapy has been linked to the development of brain cancer (most commonly astrocytomas), and Tpmt status has been associated with this risk. Therefore, we investigated whether the level of Tpmt protein activity could predict TG-associated cytotoxicity and DNA damage in astrocytic cells. We found that TG induced cytotoxicity in a dose-dependent manner in Tpmt+/+, Tpmt+/− and Tpmt−/− primary mouse astrocytes and that a low Tpmt phenotype predicted significantly higher sensitivity to TG than did a high Tpmt phenotype. We also found that TG exposure induced significantly more DNA damage in the form of single strand breaks (SSBs) and double strand breaks (DSBs) in primary astrocytes with low Tpmt versus high Tpmt. More interestingly, we found that Tpmt+/− astrocytes had the highest degree of cytotoxicity and genotoxicity (i.e., IC50, SSBs and DSBs) after TG exposure. We then used human glioma cell lines as model astroglial cells to represent high (T98) and low (A172) Tpmt expressers and found that A172 had the highest degree of cytoxicity and SSBs after TG exposure. When we over-expressed Tpmt in the A172 cell line, we found that TG IC50 was significantly higher and SSB's were significantly lower as compared to mock transfected cells. This study shows that low Tpmt can lead to greater sensitivity to thiopurine therapy in astroglial cells. When Tpmt deactivation at the germ-line is considered, this study also suggests that heterozygosity may be subject to the greatest genotoxic effects of thiopurine therapy.

Background

Implication of pharmacogenetic variations and efavirenz pharmacokinetics in concomitant efavirenz based antiviral therapy and anti-tubercular drug induced liver injury (DILI) has not been yet studied. We performed a prospective case-control association study to identify the incidence, pharmacogenetic, pharmacokinetic and biochemical predictors for anti-tubercular and antiretroviral drugs induced liver injury (DILI) in HIV and tuberculosis (TB) co-infected patients.

Methods and Findings

Newly diagnosed treatment naïve TB-HIV co-infected patients (n = 353) were enrolled to receive efavirenz based ART and rifampicin based anti-TB therapy, and assessed clinically and biochemically for DILI up to 56 weeks. Quantification of plasma efavirenz and 8-hydroxyefaviernz levels and genotyping for NAT2, CYP2B6, CYP3A5, ABCB1, UGT2B7 and SLCO1B1 genes were done. The incidence of DILI and identification of predictors was evaluated using survival analysis and the Cox Proportional Hazards Model. The incidence of DILI was 30.0%, or 14.5 per 1000 person-week, and that of severe was 18.4%, or 7.49 per 1000 person-week. A statistically significant association of DILI with being of the female sex (p = 0.001), higher plasma efavirenz level (p = 0.009), efavirenz/8-hydroxyefavirenz ratio (p = 0.036), baseline AST (p = 0.022), ALT (p = 0.014), lower hemoglobin (p = 0.008), and serum albumin (p = 0.007), NAT2 slow-acetylator genotype (p = 0.039) and ABCB1 3435TT genotype (p = 0.001).

Conclusion

We report high incidence of anti-tubercular and antiretroviral DILI in Ethiopian patients. Between patient variability in systemic efavirenz exposure and pharmacogenetic variations in NAT2, CYP2B6 and ABCB1 genes determines susceptibility to DILI in TB-HIV co-infected patients. Close monitoring of plasma efavirenz level and liver enzymes during early therapy and/or genotyping practice in HIV clinics is recommended for early identification of patients at risk of DILI.

Variable warfarin response during treatment initiation poses a significant challenge to providing optimal anticoagulation therapy. We investigated the determinants of initial warfarin response in a cohort of 167 patients. During the first nine days of treatment with pharmacogenetics-guided dosing, S-warfarin plasma levels and international normalized ratio were obtained to serve as inputs to a pharmacokinetic-pharmacodynamic (PK-PD) model. Individual PK (S-warfarin clearance) and PD (Imax) parameter values were estimated. Regression analysis demonstrated that CYP2C9 genotype, kidney function, and gender were independent determinants of S-warfarin clearance. The values for Imax were dependent on VKORC1 and CYP4F2 genotypes, vitamin K status (as measured by plasma concentrations of proteins induced by vitamin K absence, PIVKA-II) and weight. Importantly, indication for warfarin was a major independent determinant of Imax during initiation, where PD sensitivity was greater in atrial fibrillation than venous thromboembolism. To demonstrate the utility of the global PK-PD model, we compared the predicted initial anticoagulation responses with previously established warfarin dosing algorithms. These insights and modeling approaches have application to personalized warfarin therapy.

We characterized relationships between 22 polymorphisms in CYP2B6, ABCB1 and CYP3A5, plasma efavirenz exposure, and/or treatment responses in AIDS Clinical Trials Group protocols 384, A5095, and A5097s. A stepwise logistic regression procedure selected polymorphisms associated with reduced drug clearance adjusted for body mass index and composite CYP2B6 516/983 genotype. Competing risk methodology characterized relationships between selected polymorphisms and treatment responses. Association analyses involved 821 individuals (317 for pharmacokinetics, 643 for treatment response). Models that included CYP2B6 516/983 genotype best predicted pharmacokinetics. Slow metabolizer genotypes were associated with increased central nervous system events among whites, and decreased virologic failure among blacks.

Sex-differences in human liver gene expression were characterized on a genome-wide scale using a large liver sample collection, allowing for detection of small expression differences with high statistical power. 1,249 sex-biased genes were identified, 70% showing higher expression in females. Chromosomal bias was apparent, with female-biased genes enriched on chrX and male-biased genes enriched on chrY and chr19, where 11 male-biased zinc-finger KRAB-repressor domain genes are distributed in six clusters. Top biological functions and diseases significantly enriched in sex-biased genes include transcription, chromatin organization and modification, sexual reproduction, lipid metabolism and cardiovascular disease. Notably, sex-biased genes are enriched at loci associated with polygenic dyslipidemia and coronary artery disease in genome-wide association studies. Moreover, of the 8 sex-biased genes at these loci, 4 have been directly linked to monogenic disorders of lipid metabolism and show an expression profile in females (elevated expression of ABCA1, APOA5 and LDLR; reduced expression of LIPC) that is consistent with the lower female risk of coronary artery disease. Female-biased expression was also observed for CYP7A1, which is activated by drugs used to treat hypercholesterolemia. Several sex-biased drug-metabolizing enzyme genes were identified, including members of the CYP, UGT, GPX and ALDH families. Half of 879 mouse orthologs, including many genes of lipid metabolism and homeostasis, show growth hormone-regulated sex-biased expression in mouse liver, suggesting growth hormone might play a similar regulatory role in human liver. Finally, the evolutionary rate of protein coding regions for human-mouse orthologs, revealed by dN/dS ratio, is significantly higher for genes showing the same sex-bias in both species than for non-sex-biased genes. These findings establish that human hepatic sex differences are widespread and affect diverse cell metabolic processes, and may help explain sex differences in lipid profiles associated with sex differential risk of coronary artery disease.

Objective

To report pediatric cases of paradoxical respiratory adverse drug reactions (ADRs) after exposure to oral mucolytic drugs (carbocysteine, acetylcysteine) that led to the withdrawal of licenses for these drugs for infants in France and then Italy.

Design

The study followed the recommendations of the European guidelines of pharmacovigilance for medicines used in the paediatric population.

Setting

Cases voluntarily reported by physicians from 1989 to 2008 were identified in the national French pharmacovigilance public database and in drug company databases.

Patients

The definition of paradoxical respiratory ADRs was based on the literature. Exposure to mucolytic drugs was arbitrarily defined as having received mucolytic drugs for at least 2 days (>200 mg) and at least until the day before the first signs of the suspected ADR.

Results

The non-exclusive paradoxical respiratory ADRs reported in 59 paediatric patients (median age 5 months, range 3 weeks to 34 months, 98% younger than 2 years old) were increased bronchorrhea or mucus vomiting (n = 27), worsening of respiratory distress during respiratory tract infection (n = 35), dyspnoea (n = 18), cough aggravation or prolongation (n = 11), and bronchospasm (n = 1). Fifty-one (86%) children required hospitalization or extended hospitalization because of the ADR; one patient died of pulmonary oedema after mucus vomiting.

Conclusion

Parents, physicians, pharmacists, and drug regulatory agencies should know that the benefit risk ratio of mucolytic drugs is at least null and most probably negative in infants according to available evidence.

Atorvastatin δ-lactone, a major, pharmacologically inactive metabolite, has been associated with toxicity. In a previous study we showed that polymorphisms of UGT1A3 influence atorvastatin δ-lactone formation. Here we investigated the reverse reaction, atorvastatin δ-lactone hydrolysis, in a human liver bank. Screening of microarray data revealed paraoxonases PON1 and PON3 among 17 candidate esterases. Microsomal δ-lactone hydrolysis was significantly correlated to PON1 and PON3 protein (rs = 0.60; rs = 0.62, respectively; P < 0.0001). PON1 and PON3 were strongly correlated to each other (rs = 0.60) but PON1 was shown to be more extensively glycosylated than PON3. In addition a novel splice-variant of PON3 was identified. Genotyping of 40 polymorphisms within the PON-locus identified PON1 promoter polymorphisms (−108T > C, −832G > A, −1741G > A) and a tightly linked group of PON3 polymorphisms (−4984A > G, −4105G > A, −1091A > G, −746C > T, and F21F) to be associated with changes in atorvastatin δ-lactone hydrolysis and expression of PON1 but not PON3. However, carriers of the common PON1 polymorphisms L55M or Q192R showed no difference in δ-lactone hydrolysis or PON expression. Haplotype analysis revealed decreased δ-lactone hydrolysis in carriers of the most common haplotype *1 compared to carriers of haplotypes *2, *3, *4, and *7. Analysis of non-genetic factors showed association of hepatocellular and cholangiocellular carcinoma with decreased PON1 and PON3 expression, respectively. Increased C-reactive protein and γ-glutamyl transferase levels were associated with decreased protein expression of both enzymes, and increased bilirubin levels, cholestasis, and presurgical exposure to omeprazole or pantoprazole were related to decreased PON3 protein. In conclusion, PON-locus polymorphisms affect PON1 expression whereas non-genetic factors have an effect on PON1 and PON3 expression. This may influence response to therapy or adverse events in statin treatment.

MicroRNAs (miRNA) are small non-coding RNA molecules of ∼22 nucleotides which regulate large numbers of genes by binding to seed sequences at the 3′-untranslated region of target gene transcripts. The target mRNA is then usually degraded or translation is inhibited, although thus resulting in posttranscriptional down regulation of gene expression at the mRNA and/or protein level. Due to the bioinformatic difficulties in predicting functional miRNA binding sites, several publically available databases have been developed that predict miRNA binding sites based on different algorithms. The parallel use of different databases is currently indispensable, but highly uncomfortable and time consuming, especially when working with numerous genes of interest. We have therefore developed a new stand-alone program, termed MIRNA-DISTILLER, which allows to compile miRNA data for given target genes from public databases. Currently implemented are TargetScan, microCosm, and miRDB, which may be queried independently, pairwise, or together to calculate the respective intersections. Data are stored locally for application of further analysis tools including freely definable biological parameter filters, customized output-lists for both miRNAs and target genes, and various graphical facilities. The software, a data example file and a tutorial are freely available at http://www.ikp-stuttgart.de/content/language1/html/10415.asp

Background

Human serum paraoxonase-1 (PON1) prevents oxidation of low density lipoprotein cholesterol (LDL-C) and hydrolyzes the oxidized form, therefore preventing the development of atherosclerosis. The polymorphisms of PON1 gene are known to affect the PON1 activity and thereby coronary artery disease (CAD) risk. As studies are lacking in North-West Indian Punjabi's, a distinct ethnic group with high incidence of CAD, we determined PON1 activity, genotypes and haplotypes in this population and correlated them with the risk of CAD.

Methodology/Principal Findings

350 angiographically proven (≥70% stenosis) CAD patients and 300 healthy controls were investigated. PON1 activity was determined towards paraoxon (Paraoxonase; PONase) and phenylacetate (Arylesterase; AREase) substrates. In addition, genotyping was carried out by using multiplex PCR, allele specific oligonucleotide –PCR and PCR-RFLP methods and haplotyping was determined by PHASE software. The serum PONase and AREase activities were significantly lower in CAD patients as compared to the controls. All studied polymorphisms except L55M had significant effect on PONase activity. However AREase activity was not affected by them. In a logistic regression model, after adjustment for the conventional risk factors for CAD, QR (OR: 2.73 (1.57–4.72)) and RR (OR, 16.24 (6.41–41.14)) genotypes of Q192R polymorphism and GG (OR: 2.07 (1.02–4.21)) genotype of −162A/G polymorphism had significantly higher CAD risk. Haplotypes L-T-G-Q-C (OR: 3.25 (1.72–6.16)) and L-T-G-R-G (OR: 2.82 (1.01–7.80)) were also significantly associated with CAD.

Conclusions

In conclusion this study shows that CAD patients had lower PONase and AREase activities as compared to the controls. The coding Q192R polymorphism, promoter −162A/G polymorphism and L-T-G-Q-C and L-T-G-R-G haplotypes are all independently associated with CAD.

Although the efficacy of methadone maintenance treatment (MMT) in opioid dependence disorder has been well established, the influence of methadone pharmacokinetics in dose requirement and clinical outcome remains controversial. The aim of this study is to analyze methadone dosage in responder and nonresponder patients considering pharmacogenetic and pharmacokinetic factors that may contribute to dosage adequacy. Opioid dependence patients (meeting Diagnostic and Statistical Manual of Mental Disorders, [4th Edition] criteria) from a MMT community program were recruited. Patients were clinically assessed and blood samples were obtained to determine plasma concentrations of (R,S)-, (R) and (S)- methadone and to study allelic variants of genes encoding CYP3A5, CYP2D6, CYP2B6, CYP2C9, CYP2C19, and P-glycoprotein. Responders and nonresponders were defined by illicit opioid consumption detected in random urinalysis. The final sample consisted in 105 opioid dependent patients of Caucasian origin. Responder patients received higher doses of methadone and have been included into treatment for a longer period. No differences were found in terms of genotype frequencies between groups. Only CYP2D6 metabolizing phenotype differences were found in outcome status, methadone dose requirements, and plasma concentrations, being higher in the ultrarapid metabolizers. No other differences were found between phenotype and responder status, methadone dose requirements, neither in methadone plasma concentrations. Pharmacokinetic factors could explain some but not all differences in MMT outcome and methadone dose requirements.

Background

The individual character of pharmacokinetics is of great importance in the risk assessment of new drug leads in pharmacological research. Amongst others, it is severely influenced by the properties and inter-individual variability of the enzymes and transporters of the drug detoxification system of the liver. Predicting individual drug biotransformation capacity requires quantitative and detailed models.

Results

In this contribution we present the de novo deterministic modeling of atorvastatin biotransformation based on comprehensive published knowledge on involved metabolic and transport pathways as well as physicochemical properties. The model was evaluated on primary human hepatocytes and parameter identifiability analysis was performed under multiple experimental constraints. Dynamic simulations of atorvastatin biotransformation considering the inter-individual variability of the two major involved enzymes CYP3A4 and UGT1A3 based on quantitative protein expression data in a large human liver bank (n = 150) highlighted the variability in the individual biotransformation profiles and therefore also points to the individuality of pharmacokinetics.

Conclusions

A dynamic model for the biotransformation of atorvastatin has been developed using quantitative metabolite measurements in primary human hepatocytes. The model comprises kinetics for transport processes and metabolic enzymes as well as population liver expression data allowing us to assess the impact of inter-individual variability of concentrations of key proteins. Application of computational tools for parameter sensitivity analysis enabled us to considerably improve the validity of the model and to create a consistent framework for precise computer-aided simulations in toxicology.

The domestic cat (Felis catus) shows remarkable sensitivity to the adverse effects of phenolic drugs, including acetaminophen and aspirin, as well as structurally-related toxicants found in the diet and environment. This idiosyncrasy results from pseudogenization of the gene encoding UDP-glucuronosyltransferase (UGT) 1A6, the major species-conserved phenol detoxification enzyme. Here, we established the phylogenetic timing of disruptive UGT1A6 mutations and explored the hypothesis that gene inactivation in cats was enabled by minimal exposure to plant-derived toxicants. Fixation of the UGT1A6 pseudogene was estimated to have occurred between 35 and 11 million years ago with all extant Felidae having dysfunctional UGT1A6. Out of 22 additional taxa sampled, representative of most Carnivora families, only brown hyena (Parahyaena brunnea) and northern elephant seal (Mirounga angustirostris) showed inactivating UGT1A6 mutations. A comprehensive literature review of the natural diet of the sampled taxa indicated that all species with defective UGT1A6 were hypercarnivores (>70% dietary animal matter). Furthermore those species with UGT1A6 defects showed evidence for reduced amino acid constraint (increased dN/dS ratios approaching the neutral selection value of 1.0) as compared with species with intact UGT1A6. In contrast, there was no evidence for reduced amino acid constraint for these same species within UGT1A1, the gene encoding the enzyme responsible for detoxification of endogenously generated bilirubin. Our results provide the first evidence suggesting that diet may have played a permissive role in the devolution of a mammalian drug metabolizing enzyme. Further work is needed to establish whether these preliminary findings can be generalized to all Carnivora.

Cynomolgus monkey and rhesus monkey are used in drug metabolism studies due to their evolutionary closeness and physiological resemblance to human. In cynomolgus monkey, we previously identified cytochrome P450 (P450 or CYP) 2C76 that does not have a human ortholog and is partly responsible for species differences in drug metabolism between cynomolgus monkey and human. In this study, we report characterization of CYP2C93 cDNA newly identified in cynomolgus monkey and rhesus monkey. The CYP2C93 cDNA contained an open reading frame of 490 amino acids approximately 84–86% identical to human CYP2Cs. CYP2C93 was located in the genomic region, which corresponded to the intergenic region in the human genome, indicating that CYP2C93 does not correspond to any human genes. CYP2C93 mRNA was expressed predominantly in the liver among 10 tissues analyzed. The CYP2C93 proteins heterologously expressed in Escherichia coli metabolized human CYP2C substrates, diclofenac, flurbiprofen, paclitaxel, S-mephenytoin, and tolbutamide. In addition to a normal transcript (SV1), an aberrantly spliced transcript (SV2) lacking exon 2 was identified, which did not give rise to a functional protein due to frameshift and a premature termination codon. Mini gene assay revealed that the genetic variant IVS2-1G>T at the splice site of intron 1, at least partly, accounted for the exon-2 skipping; therefore, this genotype would influence CYP2C93-mediated drug metabolism. SV1 was expressed in 6 of 11 rhesus monkeys and 1 of 8 cynomolgus monkeys, but the SV1 in the cynomolgus monkey was nonfunctional due to a rare null genotype (c.102T>del). These results suggest that CYP2C93 can play roles as a drug-metabolizing enzyme in rhesus monkeys (not in cynomolgus monkeys), although its relative contribution to drug metabolism has yet to be validated.

Background

We previously showed that blood serum induced cytochrome P450 1A1 (CYP1A1) monooxygenase expression in vitro.

Objective

Our purpose was (i) to identify the molecular mechanism involved and (ii) to characterize the inducer compound(s) in serum involved at least in part.

Methods

Serum was fractionated on hydrophobic columns. PPARα involvement was demonstrated by gene reporter assays, DNA mutagenesis and EMSA. Gene expression was evaluated by qRT-PCR. Serum samples were analyzed using HS-SPME-GC-MS.

Results

The inductive effect of serum did not depend on the AhR pathway and was enhanced by cotransfection of PPARα cDNA. Mutations in the PPAR response elements of the CYP1A1 gene promoter suppressed this effect. One of the PPRE sites appeared highly specific for human PPARα, an unreported PPRE property. A link was found between CYP1A1 inducibility and serum hydrophobic compounds. Characterization of sera showed that hexanal, a metabolite produced by peroxidation of linoleic acid, was involved in CYP1A1 induction by serum, possibly along with other serum entities.

Conclusion

We demonstrate that serum induces CYP1A1 via the PPARα pathway and that hexanal is one of the serum inducers. The two PPRE sites within the CYP1A1 promoter are functional and one of them is specific for PPARα.

We have previously shown that genetic variability in CNR1 is associated with low HDL dyslipidemia in a multigenerational obesity study cohort of Northern European descent (209 families, median  = 10 individuals per pedigree). In order to assess the impact of CNR1 variability on the development of dyslipidemia in the community, we genotyped this locus in all subjects with class III obesity (body mass index >40 kg/m2) participating in a population-based biobank of similar ancestry. Twenty-two haplotype tagging SNPs, capturing the entire CNR1 gene locus plus 15 kb upstream and 5 kb downstream, were genotyped and tested for association with clinical lipid data. This biobank contains data from 645 morbidly obese study subjects. In these subjects, a common CNR1 haplotype (H3, frequency 21.1%) is associated with fasting TG and HDL cholesterol levels (p = 0.031 for logTG; p = 0.038 for HDL-C; p = 0.00376 for log[TG/HDL-C]). The strength of this relationship increases when the data are adjusted for age, gender, body mass index, diet and physical activity. Mean TG levels were 160±70, 155±70, and 120±60 mg/dL for subjects with 0, 1, and 2 copies of the H3 haplotype. Mean HDL-C levels were 45±10, 47±10, and 48±9 mg/dL, respectively. The H3 CNR1 haplotype appears to exert a protective effect against development of obesity-related dyslipidemia.

Background

Tissue-specific expression of CYP450s can regulate the intracellular concentration of drugs and explain inter-subject variability in drug action. The overall objective of our study was to determine in a large cohort of samples, mRNA levels and CYP450 activity expressed in the human heart.

Methodology

CYP450 mRNA levels were determined by RTPCR in left ventricular samples (n = 68) of explanted hearts from patients with end-stage heart failure. Samples were obtained from ischemic and non-ischemic hearts. In some instances (n = 7), samples were available from both the left and right ventricles. A technique for the preparation of microsomes from human heart tissue was developed and CYP450-dependent activity was determined using verapamil enantiomers as probe-drug substrates.

Principal Findings

Our results show that CYP2J2 mRNA was the most abundant isoform in all human heart left ventricular samples tested. Other CYP450 mRNAs of importance were CYP4A11, CYP2E1, CYP1A1 and CYP2C8 mRNAs while CYP2B6 and CYP2C9 mRNAs were present at low levels in only some of the hearts analyzed. CYP450 mRNAs did not differ between ischemic and non-ischemic hearts and appeared to be present at similar levels in the left and right ventricles. Incubation of verapamil with heart microsomes led to the formation of nine CYP450-dependent metabolites: a major finding was the observation that stereoselectivity was reversed compared to human liver microsomes, in which the R-enantiomer is metabolized to a greater extent.

Conclusions

This study determined cardiac mRNA levels of various CYP450 isozymes involved in drug metabolism and demonstrated the prevalent expression of CYP2J2 mRNA. It revealed that cardiomyocytes can efficiently metabolize drugs and that cardiac CYP450s are highly relevant with regard to clearance of drugs in the heart. Our results support the claim that drug metabolism in the vicinity of a drug effector site can modulate drug effects.

The human drug metabolizing cytochrome P450 (CYP) 1A2, is one of the major P450 isoforms contributing by about 5–20% to the hepatic P450 pool and catalyzing oxidative biotransformation of up to 10% of clinically relevant drugs including clozapine and caffeine. CYP1A2 activity is interindividually highly variable and although twin studies have suggested a high heritability, underlying genetic factors are still unknown. Here we adopted a pathway-oriented approach using a large human liver bank (n = 150) to elucidate whether variants in candidate genes of constitutive, ligand-inducible, and pathophysiological inhibitory regulatory pathways may explain different hepatic CYP1A2 phenotypes. Samples were phenotyped for phenacetin O-deethylase activity, and the expression of CYP1A2 protein and mRNA was determined. CYP1A2 expression and function was increased in smokers and decreased in patients with inflammation and cholestasis. Of 169 SNPs in 17 candidate genes including the CYP1A locus, 136 non-redundant SNPs with minor allele frequency >5% were analyzed by univariate and multivariate methods. A total of 13 strong significant associations were identified, of which 10 SNPs in the ARNT, AhRR, HNF1α, IL1β, SRC-1, and VDR genes showed consistent changes for at least two phenotypes by univariate analysis. Multivariate linear modeling indicated that the polymorphisms and non-genetic factors together explained 42, 38, and 33% of CYP1A2 variation at activity, protein and mRNA levels, respectively. In conclusion, we identified novel trans-associations between regulatory genes and hepatic CYP1A2 function and expression, but additional genetic factors must be assumed to explain the full extent of CYP1A2 heritability.

Background

Flavopiridol is a cyclin-dependent kinase inhibitor in phase II clinical development for treatment of various forms of cancer. When administered with a pharmacokinetically (PK)-directed dosing schedule, flavopiridol exhibited striking activity in patients with refractory chronic lymphocytic leukemia. This study aimed to evaluate pharmacogenetic factors associated with inter-individual variability in pharmacokinetics and outcomes associated with flavopiridol therapy.

Methodology/Principal Findings

Thirty-five patients who received single-agent flavopiridol via the PK-directed schedule were genotyped for 189 polymorphisms in genes encoding 56 drug metabolizing enzymes and transporters. Genotypes were evaluated in univariate and multivariate analyses as covariates in a population PK model. Transport of flavopiridol and its glucuronide metabolite was evaluated in uptake assays in HEK-293 and MDCK-II cells transiently transfected with SLCO1B1. Polymorphisms in ABCC2, ABCG2, UGT1A1, UGT1A9, and SLCO1B1 were found to significantly correlate with flavopiridol PK in univariate analysis. Transport assay results indicated both flavopiridol and flavopiridol-glucuronide are substrates of the SLCO1B1/OATP1B1 transporter. Covariates incorporated into the final population PK model included bilirubin, SLCO1B1 rs11045819 and ABCC2 rs8187710. Associations were also observed between genotype and response. To validate these findings, a second set of data with 51 patients was evaluated, and overall trends for associations between PK and PGx were found to be consistent.

Conclusions/Significance

Polymorphisms in transport genes were found to be associated with flavopiridol disposition and outcomes. Observed clinical associations with SLCO1B1 were functionally validated indicating for the first time its relevance as a transporter of flavopiridol and its glucuronide metabolite. A second 51-patient dataset indicated similar trends between genotype in the SLCO1B1 and other candidate genes, thus providing support for these findings. Further study in larger patient populations will be necessary to fully characterize and validate the clinical impact of polymorphisms in SLCO1B1 and other transporter and metabolizing enzyme genes on outcomes from flavopiridol therapy.

Adverse drug reactions (ADR), also known as side-effects, are complex undesired physiologic phenomena observed secondary to the administration of pharmaceuticals. Several phenomena underlie the emergence of each ADR; however, a dominant factor is the drug's ability to modulate one or more biological pathways. Understanding the biological processes behind the occurrence of ADRs would lead to the development of safer and more effective drugs. At present, no method exists to discover these ADR-pathway associations. In this paper we introduce a computational framework for identifying a subset of these associations based on the assumption that drugs capable of modulating the same pathway may induce similar ADRs. Our model exploits multiple information resources. First, we utilize a publicly available dataset pairing drugs with their observed ADRs. Second, we identify putative protein targets for each drug using the protein structure database and in-silico virtual docking. Third, we label each protein target with its known involvement in one or more biological pathways. Finally, the relationships among these information sources are mined using multiple stages of logistic-regression while controlling for over-fitting and multiple-hypothesis testing. As proof-of-concept, we examined a dataset of 506 ADRs, 730 drugs, and 830 human protein targets. Our method yielded 185 ADR-pathway associations of which 45 were selected to undergo a manual literature review. We found 32 associations to be supported by the scientific literature.