Nathan J Cherrington
Associate Dean, Research and Graduate Studies - College of Pharmacy
Director, Southwest Environmental Health Science Center
Professor, BIO5 Institute
Professor, Clinical Translational Sciences
Professor, Pharmacology and Toxicology
Professor, Public Health
Primary Department
Department Affiliations
(520) 626-0219
Research Interest
Numerous drug-induced and environmental exposure-related toxicities are the result of inter-individual variation in the ADME processes of absorption, distribution, metabolism and elimination that control the fate of these compounds from the body. Alterations in these processes provide the mechanistic basis for individual variability in response to drugs and environmental exposures. A common perception is that variability in response is due to genetic polymorphisms within the drug metabolizing enzyme and transporter genes. While there are numerous examples of these differences that play a major role in the susceptibility of genetic subpopulations for specific toxicities, the potential for transient phenotypic conversion due to temporary environmental changes, such as inflammation and disease, are often overlooked.Due to the ensuing liver damage caused by the progressive stages of NAFLD, gene expression patterns can change dramatically resulting in a phenoconversion resembling genetic polymorphisms. Because the liver plays such a key role in the metabolism and disposition of xenobiotics, this temporary phenoconversion could lead to the inability of patients to properly metabolize and excrete drugs and environmental toxicants, increasing the risk of some adverse drug reactions and environmental toxicities.


Han, J., Dzierlenga, A. L., Lu, Z., Billheimer, D. D., Torabzadeh, E., Lake, A. D., Li, H., Novak, P., Shipkova, P., Aranibar, N., Robertson, D., Reily, M. D., Lehman-McKeeman, L. D., & Cherrington, N. J. (2017). Metabolomic profiling distinction of human nonalcoholic fatty liver disease progression from a common rat model. Obesity (Silver Spring, Md.), 25(6), 1069-1076.
BIO5 Collaborators
Dean Billheimer, Nathan J Cherrington

Characteristic pathological changes define the progression of steatosis to nonalcoholic steatohepatitis (NASH) and are correlated to metabolic pathways. A common rodent model of NASH is the methionine and choline deficient (MCD) diet. The objective of this study was to perform full metabolomic analyses on liver samples to determine which pathways are altered most pronouncedly in this condition in humans, and to compare these changes to rodent models of nonalcoholic fatty liver disease (NAFLD).

Lake, A. D., Novak, P., Fisher, C. D., Jackson, J. P., Hardwick, R. N., Billheimer, D. D., Klimecki, W. T., & Cherrington, N. J. (2011). Analysis of global and absorption, distribution, metabolism, and elimination gene expression in the progressive stages of human nonalcoholic fatty liver disease. Drug Metabolism and Disposition, 39(10), 1954-1960.
BIO5 Collaborators
Dean Billheimer, Nathan J Cherrington, Walter Klimecki

PMID: 21737566;PMCID: PMC3186211;Abstract:

Nonalcoholic fatty liver disease (NAFLD) is characterized by a series of pathological changes that range from simple fatty liver to nonalcoholic steatohepatitis (NASH). The objective of this study is to describe changes in global gene expression associated with the progression of human NAFLD. This study is focused on the expression levels of genes responsible for the absorption, distribution, metabolism, and elimination (ADME) of drugs. Differential gene expression between three clinically defined pathological groups - normal, steatosis, and NASH - was analyzed. Genome-wide mRNA levels in samples of human liver tissue were assayed with Affymetrix GeneChip Human 1.0ST arrays. A total of 11,633 genes exhibited altered expression out of 33,252 genes at a 5% false discovery rate. Most gene expression changes occurred in the progression from steatosis to NASH. Principal component analysis revealed that hepatic disease status was the major determinant of differential ADME gene expression rather than age or sex of sample donors. Among the 515 drug transporters and 258 drug-metabolizing enzymes (DMEs) examined, uptake transporters but not efflux transporters or DMEs were significantly over-represented in the number of genes down-regulated. These results suggest that uptake transporter genes are coordinately targeted for down-regulation at the global level during the pathological development of NASH and that these patients may have decreased drug uptake capacity. This coordinated regulation of uptake transporter genes is indicative of a hepatoprotective mechanism acting to prevent accumulation of toxic intermediates in disease- compromised hepatocytes. Copyright © 2011 by The American Society for Pharmacology and Experimental Therapeutics.

Li, H., Canet, M. J., Clarke, J. D., Billheimer, D., Xanthakos, S. A., Lavine, J. E., Erickson, R. P., & Cherrington, N. J. (2017). Pediatric Cytochrome P450 Activity Alterations in Nonalcoholic Steatohepatitis. Drug metabolism and disposition: the biological fate of chemicals, 45(12), 1317-1325.
BIO5 Collaborators
Dean Billheimer, Nathan J Cherrington

Variable drug responses depend on individual variation in the activity of drug-metabolizing enzymes, including cytochrome P450 enzymes (CYP). As the most common chronic liver disease in children and adults, nonalcoholic steatohepatitis (NASH) has been identified as a source of significant interindividual variation in hepatic drug metabolism. Compared with adults, children present age-related differences in pharmacokinetics and pharmacodynamics. The purpose of this study was to determine the impact of fatty liver disease severity on the activity of a variety of CYP enzymes in children and adolescents. Healthy and nonalcoholic fatty liver disease pediatric subjects aged 12-21 years inclusive received an oral cocktail of four probe drugs: caffeine (CYP1A2, 100 mg), omeprazole (CYP2C19, 20 mg), losartan (CYP2C9, 25 mg), and midazolam (CYP3A4, 2 mg). Venous blood and urine were collected before administration and 1, 2, 4, and 6 hours after administration. Concentrations of the parent drugs and CYP-specific metabolites were quantified in plasma and urine using liquid chromatography with tandem mass spectrometry. In plasma, the decreased metabolic area under the curve (AUC) ratio, defined as the metabolite AUC to parent AUC, of omeprazole indicated significant decreases of CYP2C19 (P = 0.002) enzymatic activities in NASH adolescents, while the urine analyses did not show significant differences and were highly variable. A comparison between the present in vivo pediatric studies and a previous ex vivo study in adults indicates distinct differences in the activities of CYP1A2 and CYP2C9. These data demonstrate that pediatric NASH presents an altered pattern of CYP activity and NASH should be considered as a confounder of drug metabolism for certain CYP enzymes. These differences could lead to future investigations that may reveal unexpected variable drug responses that should be considered in pediatric dosage recommendations.

Beilke, L. D., Besselsen, D. G., Cheng, Q., Kulkarni, S., Slitt, A. L., & Cherrington, N. J. (2008). Minimal role of hepatic transporters in the hepatoprotection against LCA-induced intrahepatic cholestasis. Toxicological sciences : an official journal of the Society of Toxicology, 102(1), 196-204.
BIO5 Collaborators
David G Besselsen, Nathan J Cherrington

The multidrug resistance-associated proteins (Mrps) are a family of adenosine triphosphate-dependent transporters that facilitate the movement of various compounds, including bile acids, out of hepatocytes. The current study was conducted to determine whether induction of these transporters alters bile acid disposition as a means of hepatoprotection during bile acid-induced cholestasis. Lithocholic acid (LCA) was used to induce intrahepatic cholestasis. C57BL/6 mice were pretreated with corn oil (CO) or known transporter inducers, phenobarbital (PB), oltipraz (OPZ), or TCPOBOP (TC) for 3 days prior to cotreatment with LCA and inducer for 4 days. Histopathology revealed that PB and TC pretreatments provide a protective effect from LCA-induced toxicity, whereas OPZ pretreatment did not. Both PB/LCA and TC/LCA cotreatment groups also had significantly lower alanine aminotransferase values than the LCA-only group. In TC/LCA cotreated mice compared with LCA only, messenger RNA (mRNA) expression of uptake transporters Ntcp and Oatp4 was significantly increased, as were sinusoidal efflux transporters Mrp3 and Mrp4. Although in PB/LCA cotreated mice, the only significant change compared with LCA-only treatment was an increase in uptake transporter Oatp4. Oatp1 was reduced in all groups compared with CO controls. No significant changes in mRNA expression were observed in Oatp2, Bsep, Mrp2, Bcrp, Mrp1, Mrp5, or Mrp6. Mrp4 protein expression was induced in the OPZ/LCA and TC/LCA cotreated groups, whereas Mrp3 protein levels remained unchanged between groups. Protein expression of Mrp1 and Mrp5 was increased in the unprotected LCA-only and OPZ/LCA mice. Thus, transporter expression did not correlate with histologic hepatoprotection, however, there was a correlation between hepatoprotection and significantly reduced total liver bile acids in the PB/LCA and TC/LCA cotreated mice compared with LCA only. In conclusion, changes in transporter expression did not correlate with hepatoprotection, and therefore, transport may not play a critical role in the observed hepatoprotection from LCA-induced cholestasis in the C57BL/6 mouse.

Beilke, L. D., Aleksunes, L. M., Holland, R. D., Besselsen, D. G., Beger, R. D., Klaassen, C. D., & Cherrington, N. J. (2009). Constitutive androstane receptor-mediated changes in bile acid composition contributes to hepatoprotection from lithocholic acid-induced liver injury in mice. Drug metabolism and disposition: the biological fate of chemicals, 37(5), 1035-45.
BIO5 Collaborators
David G Besselsen, Nathan J Cherrington

Pharmacological activation of the constitutive androstane receptor (CAR) protects the liver during cholestasis. The current study evaluates how activation of CAR influences genes involved in bile acid biosynthesis as a mechanism of hepatoprotection during bile acid-induced liver injury. CAR activators phenobarbital (PB) and 1,4-bis[2-(3,5-dichloropyridyloxy)]benzene (TCPOBOP) or corn oil (CO) were administered to C57BL/6 wild-type (WT) and CAR knockout (CAR-null) mice before and during induction of intrahepatic cholestasis using the secondary bile acid, lithocholic acid (LCA). In LCA-treated WT and all the CAR-null groups (excluding controls), histology revealed severe multifocal necrosis. This pathology was absent in WT mice pretreated with PB and TCPOBOP, indicating CAR-dependent hepatoprotection. Decreases in total hepatic bile acids and hepatic monohydroxy, dihydroxy, and trihydroxy bile acids in PB- and TCPOBOP-pretreated WT mice correlated with hepatoprotection. In comparison, concentrations of monohydroxylated and dihydroxylated bile acids were increased in all the treated CAR-null mice compared with CO controls. Along with several other enzymes (Cyp7b1, Cyp27a1, Cyp39a1), Cyp8b1 expression was increased in hepatoprotected mice, which could be suggestive of a shift in the bile acid biosynthesis pathway toward the formation of less toxic bile acids. In CAR-null mice, these changes in gene expression were not different among treatment groups. These results suggest CAR mediates a shift in bile acid biosynthesis toward the formation of less toxic bile acids, as well as a decrease in hepatic bile acid concentrations. We propose that these combined CAR-mediated effects may contribute to the hepatoprotection observed during LCA-induced liver injury.