The ADME characteristics of siRNA therapeutics and the opportunity to predict disposition in pregnant women.

Small interfering RNA (siRNA) therapeutics represent an emerging class of pharmacotherapy with the potential to address previously hard-to-treat diseases. Currently approved siRNA therapeutics include LNP-encapsulated siRNA and triGalNAc-conjugated siRNA. These siRNA therapeutics exhibit distinct pharmacokinetic characteristics and unique absorption, distribution, metabolism, and elimination (ADME) properties. As a new drug modality, limited clinical data are available for siRNA therapeutics in specific populations, including pediatrics, geriatrics, individuals with renal or hepatic impairment, and pregnant women, making dosing challenging. In this review, a mechanistic overview of the ADME properties of the five currently approved siRNA therapeutics is presented. A concise overview of the clinical data available for therapeutic siRNAs in special populations, focusing on the potential impact of physiological changes during pregnancy on siRNA disposition is provided. The utility of physiologically based pharmacokinetic (PBPK) modeling as a tool to elucidate the characteristics and disposition of siRNA therapeutics in pregnant women is explored. Additionally, opportunities to integrate known physiological alterations induced by pregnancy into PBPK models that incorporate siRNA ADME mechanisms to predict the effects of pregnancy on siRNA disposition are discussed. Clinical data regarding the use of therapeutic siRNA in special populations remains limited. Data for precise parameterization of maternal-fetal siRNA PBPK models is lacking presently and underscores the need for further research in this area. Addressing this gap in knowledge will not only enhance our understanding of siRNA pharmacokinetics during pregnancy but also advance possible development of siRNA therapeutics to treat pregnancy related conditions. Significance Statement This review proposes a framework on how siRNA disposition can be predicted in pregnancy based on mechanistic ADME information using physiologically based pharmacokinetic (PBPK) modeling. The mechanistic ADME information and available clinical data in special populations of currently FDA approved siRNA therapeutics are summarized. A detailed discussion on how physiological changes during pregnancy may affect siRNA disposition in pregnant women and on the opportunities to project siRNA disposition in pregnant women using PBPK modeling is provided.

CYP2C9, CYP3A and CYP2C19 metabolize Δ9-tetrahydrocannabinol to multiple metabolites but metabolism is affected by human liver fatty acid binding protein (FABP1).

Δ9-tetrahydrocannabinol (THC) is the psychoactive constituent of cannabis. It is cleared predominantly via metabolism. Metabolism to 11-OH-THC by cytochrome P450 (CYP) 2C9 has been proposed as the main clearance pathway of THC, with the estimated fraction metabolized (fm) about 70%. The remaining clearance pathways are not well established, and it is unknown how THC is eliminated in individuals with reduced CYP2C9 activity. The goal of this study was to systematically identify the CYP enzymes contributing to THC clearance and characterize the metabolites formed. Further, this study aimed to characterize the impact of liver fatty acid binding protein (FABP1) on THC metabolism by human CYPs. THC was metabolized to at least four different metabolites including 11-OH-THC in human liver microsomes (HLMs) and with recombinant CYPs. 11-OH-THC was formed by recombinant CYP2C9 (Km,u = 0.77 nM, kcat = 12 min-1) and by recombinant CYP2C19 (Km,u = 2.2 nM, kcat = 14 min-1). The other three major metabolites were likely hydroxylations in the cyclohexenyl ring and were formed mainly by recombinant CYP3A4/5 (Km,u > 10 nM). HLM experiments confirmed the contributions of CYP2C9, CYP2C19 and CYP3A to THC metabolism. The presence of FABP1 and THC binding to FABP1 altered THC metabolism by recombinant CYPs and HLMs in an enzyme and metabolite specific manner. This suggests that FABP1 may interact with CYP enzymes and alter the fm by CYPs towards THC metabolism. In conclusion, this study is the first to systematically establish the metabolic profile of THC by human CYPs and characterize how FABP1 binding alters CYP mediated THC metabolism.

Drugs Form Ternary Complexes with Human Liver Fatty Acid Binding Protein 1 (FABP1) and FABP1 Binding Alters Drug Metabolism.

Liver fatty acid binding protein 1 (FABP1) binds diverse endogenous lipids and is highly expressed in the human liver. Binding to FABP1 alters the metabolism and homeostasis of endogenous lipids in the liver. Drugs have also been shown to bind to rat FABP1, but limited data are available for human FABP1 (hFABP1). FABP1 has a large binding pocket, and up to two fatty acids can bind to FABP1 simultaneously. We hypothesized that drug binding to hFABP1 results in formation of ternary complexes and that FABP1 binding alters drug metabolism. To test these hypotheses, native protein mass spectrometry (MS) and fluorescent 11-(dansylamino)undecanoic acid (DAUDA) displacement assays were used to characterize drug binding to hFABP1, and diclofenac oxidation by cytochrome P450 2C9 (CYP2C9) was studied in the presence and absence of hFABP1. DAUDA binding to hFABP1 involved high (Kd,1 = 0.2 µM) and low (Kd,2 > 10 µM) affinity binding sites. Nine drugs bound to hFABP1 with equilibrium dissociation constant (Kd) values ranging from 1 to 20 µM. None of the tested drugs completely displaced DAUDA from hFABP1, and fluorescence spectra showed evidence of ternary complex formation. Formation of DAUDA-hFABP1-diclofenac ternary complex was verified with native MS. Docking predicted diclofenac binding in the portal region of FABP1 with DAUDA in the binding cavity. The catalytic rate constant of diclofenac hydroxylation by CYP2C9 was decreased by ∼50% (P < 0.01) in the presence of FABP1. Together, these results suggest that drugs form ternary complexes with hFABP1 and that hFABP1 binding in the liver will alter drug metabolism and clearance. SIGNIFICANCE STATEMENT: Many commonly prescribed drugs bind fatty acid binding protein 1 (FABP1), forming ternary complexes with FABP1 and the fluorescent fatty acid 11-(dansylamino)undecanoic acid. These findings suggest that drugs will bind to apo-FABP1 and fatty acid-bound FABP1 in the human liver. The high expression of FABP1 in the liver, together with drug binding to FABP1, may alter drug disposition processes in vivo.

LX-2 Stellate Cells Are a Model System for Investigating the Regulation of Hepatic Vitamin A Metabolism and Respond to Tumor Necrosis Factor α and Interleukin 1ß.

Hepatic stellate cells (HSCs) are the major site of vitamin A (retinol) esterification and subsequent storage as retinyl esters within lipid droplets. However, retinyl esters become depleted in many pathophysiological states, including acute and chronic liver injuries. Recently, using a liver slice culture system as a model of acute liver injury and fibrogenesis, a time-dependent increase and decrease in the apparent formation of the bioactive retinoid all-trans-retinoic acid (atRA) and retinyl palmitate was measured, respectively. This coincided with temporal changes in the gene expression of retinoid-metabolizing enzymes and binding proteins, that preceded HSC activation. However, the underlying mechanisms that promote early changes in retinoid metabolism remain unresolved. We hypothesized that LX-2 cells could be applied to investigate differences in quiescent and activated HSC retinoid metabolism. We demonstrate that the hypermetabolic state of activated stellate cells relative to quiescent stellate cells may be attributed to induction of STRA6, RBP4, and CYP26A1, thereby reducing intracellular concentrations of atRA. We further hypothesized that paracrine and autocrine cytokine signaling regulates HSC vitamin A metabolism in both quiescent and activated cells. In quiescent cells, tumor necrosis factor α dose-dependently downregulated LRAT and CRBP1 mRNA, with EC50 values of 30-50 pg/mL. Likewise, interleukin-1ß decreased LRAT and CRBP1 gene expression but with less potency. In activated stellate cells, multiple enzymes were downregulated, suggesting that the full effects of altered hepatic vitamin A metabolism in chronic conditions require both paracrine and autocrine signaling events. Further, this study suggests the potential for cell type-specific autocrine effects in hepatic retinoid signaling. SIGNIFICANCE STATEMENT: HSCs are the major site of vitamin A storage and important determinants of retinol metabolism during liver fibrogenesis. Here, two LX-2 culture methods were applied as models of hepatic retinoid metabolism to demonstrate the effects of activation status and dose-dependent cytokine exposure on the expression of genes involved in retinoid metabolism. This study suggests that compared to quiescent cells, activated HSCs are hypermetabolic and have reduced apparent formation of retinoic acid, which may alter downstream retinoic acid signaling.

Physiologically Based Pharmacokinetic (PBPK) Modeling of small molecules: How Much Progress Have We Made?

Physiologically based pharmacokinetic (PBPK) models of small molecules have become mainstream in drug development and academic research. The use of PBPK models is continuously expanding with the majority of work now focusing on predictions of drug-drug interactions, drug-disease interactions, and changes in drug disposition across lifespan. Recently, publications that use PBPK modeling to predict drug disposition during pregnancy and in organ impairment have increased reflecting the advances in incorporating diverse physiological changes into the models. Due to the expanding computational power and diversity of modeling platforms available, the complexity of PBPK models has also increased. Academic efforts have provided clear advances in better capturing human physiology in PBPK models and incorporating more complex mathematical concepts into PBPK models. Examples of such advances include the segregated gut model with a series gut compartments allowing modeling of physiological blood flow distribution within an organ and zonation of metabolic enzymes, and series compartment liver models allowing simulations of hepatic clearance for high extraction drugs. Despite these advances in academic research, the progress in assessing model quality and defining model acceptance criteria based on the intended use of the models has not kept pace. This review suggests that awareness of the need for predefined criteria for model acceptance has increased but many manuscripts still lack description of scientific justification and/or rationale for chosen acceptance criteria. As artificial intelligence and machine learning approaches become more broadly accepted, these tools offer promise for development of comprehensive assessment for existing observed data and analysis of model performance. Significance Statement PBPK modeling has become a mainstream application in academic literature and is broadly used for predictions, analysis and evaluation of pharmacokinetic data. Many significant advances have been made in developing advanced PBPK models that better capture human physiology but oftentimes sufficient justification for the chosen model acceptance criterion and model structure is still missing. This review provides a summary of the current landscape of PBPK applications used and highlights the needs for advancing PBPK modeling science and training in academia.

The impact of pregnancy and associated hormones on the pharmacokinetics of Δ9-tetrahydrocannabinol.

INTRODUCTION: (-)-Δ9-tetrahydrocannabinol (THC) is the main psychoactive component of cannabis. Cannabis is the most widely used drug of abuse by pregnant individuals, but its maternal-fetal safety is still unclear. The changes in THC disposition during pregnancy may affect THC safety and pharmacology. AREAS COVERED: This review summarizes the current literature on THC metabolism and pharmacokinetics in humans. It provides an analysis of how hormonal changes during pregnancy may alter the expression of cannabinoid metabolizing enzymes and THC and its metabolite pharmacokinetics. THC is predominately (>70%) cleared by hepatic metabolism to its psychoactive active metabolite, 11-OH-THC by cytochrome P450 (CYP) 2C9 and to other metabolites (<30%) by CYP3A4. Other physiological processes that change during pregnancy and may alter cannabinoid disposition are also reviewed. EXPERT OPINION: THC and its metabolites disposition likely change during pregnancy. Hepatic CYP2C9 and CYP3A4 are induced in pregnant individuals and in vitro by pregnancy hormones. This induction of CYP2C9 and CYP3A4 is predicted to lead to altered THC and 11-OH-THC disposition and pharmacodynamic effects. More in vitro studies of THC metabolism and induction of the enzymes metabolizing cannabinoids are necessary to improve the prediction of THC pharmacokinetics in pregnant individuals.

Drugs Form Ternary Complexes with Human Liver Fatty Acid Binding Protein (FABP1) and FABP1 Binding Alters Drug Metabolism.

Liver fatty acid binding protein (FABP1) binds diverse endogenous lipids and is highly expressed in the human liver. Binding to FABP1 alters the metabolism and homeostasis of endogenous lipids in the liver. Drugs have also been shown to bind to rat FABP1, but limited data is available for human FABP1 (hFABP1). FABP1 has a large binding pocket and multiple fatty acids can bind to FABP1 simultaneously. We hypothesized that drug binding to hFABP1 results in formation of ternary complexes and that FABP1 binding alters drug metabolism. To test these hypotheses native protein mass spectrometry (MS) and fluorescent 11-(dansylamino)undecanoic acid (DAUDA) displacement assays were used to characterize drug binding to hFABP1 and diclofenac oxidation by cytochrome P450 2C9 (CYP2C9) was studied in the presence and absence of hFABP1. DAUDA binding to hFABP1 involved high (Kd,1=0.2 µM) and low affinity (Kd,2 >10 µM) binding sites. Nine drugs bound to hFABP1 with Kd values ranging from 1 to 20 µM. None of the tested drugs completely displaced DAUDA from hFABP1 and fluorescence spectra showed evidence of ternary complex formation. Formation of DAUDA-diclofenac-hFABP1 ternary complex was verified with native MS. Docking placed diclofenac in the portal region of FABP1 with DAUDA in the binding cavity. Presence of hFABP1 decreased the kcat and Km,u of diclofenac with CYP2C9 by ~50% suggesting that hFABP1 binding in the liver will alter drug metabolism and clearance. Together, these results suggest that drugs form ternary complexes with hFABP1 and that hFABP1 interacts with CYP2C9.

Predicting complex kidney drug handling using a physiologically-based pharmacokinetic model informed by biomarker-estimated secretory clearance and blood flow.

Kidney function-adjusted drug dosing is currently based solely on the estimated glomerular filtration rate (GFR), however, kidney drug handling is accomplished by a combination of filtration, tubular secretion, and re-absorption. Mechanistic physiologically-based pharmacokinetic (PBPK) models recapitulate anatomic compartments to predict elimination from estimated perfusion, filtration, secretion, and re-absorption, but clinical applications are limited by a lack of empiric individual-level measurements of these functions. We adapted and validated a PBPK model to predict drug clearance from individual biomarker-based estimates of kidney perfusion and secretory clearance. We estimated organic anion transporter-mediated secretion via kynurenic acid clearance and kidney blood flow (KBF) via isovalerylglycine clearance in human participants, incorporating these measurements with GFR into the model to predict kidney drug clearance. We compared measured and model-predicted clearances of administered tenofovir and oseltamivir, which are cleared by both filtration and secretion. There were 27 outpatients (age 55 ± 15 years, mean iohexol-GFR [iGFR] 76 ± 31 mL/min/1.73 m2 ) in this drug clearance study. The mean observed and mechanistic model-predicted tenofovir clearances were 169 ± 102 mL/min and 163 ± 80 mL/min, respectively; estimated mean error of the mechanistic model was 37.1 mL/min (95% confidence interval [CI]: 24-52.9), compared to a mean error of 41.8 mL/min (95% CI: 25-61.6) from regression model. The mean observed and model-predicted oseltamivir carboxylate clearances were 183 ± 104 mL/min and 179 ± 89 mL/min, respectively; estimated mean error of the mechanistic model was 42.9 mL/min (95% CI: 29.7-56.4), versus error of 48.1 mL/min (95% CI: 31.2-67.3) from the regression model. Individualized estimates of tubular secretion and KBF improved the accuracy of PBPK model-predicted tenofovir and oseltamivir kidney clearances, suggesting the potential for biomarker-informed measures of kidney function to refine personalized drug dosing.

Effect of isotretinoin on CYP2D6 and CYP3A activity in patients with severe acne.

AIMS: Previously, retinoids have decreased CYP2D6 mRNA expression in vitro and induced CYP3A4 in vitro and in vivo. This study aimed to determine whether isotretinoin administration changes CYP2D6 and CYP3A activities in patients with severe acne. METHODS: Thirty-three patients (22 females and 11 males, 23.5 ± 6.0 years old) expected to receive isotretinoin treatment completed the study. All participants were genotyped for CYP2D6 and CYP3A5. Participants received dextromethorphan (DM) 30 mg orally as a dual-probe substrate of CYP2D6 and CYP3A activity at two study timepoints: pre-isotretinoin treatment and with isotretinoin for at least 1 week. The concentrations of isotretinoin, DM and their metabolites were measured in 2-h postdose plasma samples and in cumulative 0-4-h urine collections using liquid chromatography-mass spectrometry. RESULTS: In CYP2D6 extensive metabolizers, the urinary dextrorphan (DX)/DM metabolic ratio (MR) (CYP2D6 activity marker) was numerically, but not significantly, lower with isotretinoin administration compared to pre-isotretinoin (geometric mean ratio [GMR] [90% confidence interval (CI)] 0.78 [0.55, 1.11]). The urinary 3-hydroxymorphinan (3HM)/DX MR (CYP3A activity marker) was increased (GMR 1.18 [1.03, 1.35]) and the urinary DX-O-glucuronide/DX MR (proposed UGT2B marker) was increased (GMR 1.22 [1.06, 1.39]) with isotretinoin administration compared to pre-isotretinoin. CONCLUSIONS: Administration of isotretinoin did not significantly reduce CYP2D6 activity in extensive metabolizers, suggesting that the predicted downregulation of CYP2D6 based on in vitro data does not translate into humans. We observed a modest increase in CYP3A activity (predominantly CYP3A4) with isotretinoin treatment. The data also suggest that DX glucuronidation is increased following isotretinoin administration.

Raising NAD+ Level Stimulates Short-Chain Dehydrogenase/Reductase Proteins to Alleviate Heart Failure Independent of Mitochondrial Protein Deacetylation.

BACKGROUND: Strategies to increase cellular NAD+ (oxidized nicotinamide adenine dinucleotide) level have prevented cardiac dysfunction in multiple models of heart failure, but molecular mechanisms remain unclear. Little is known about the benefits of NAD+-based therapies in failing hearts after the symptoms of heart failure have appeared. Most pretreatment regimens suggested mechanisms involving activation of sirtuin, especially Sirt3 (sirtuin 3), and mitochondrial protein acetylation. METHODS: We induced cardiac dysfunction by pressure overload in SIRT3-deficient (knockout) mice and compared their response with nicotinamide riboside chloride treatment with wild-type mice. To model a therapeutic approach, we initiated the treatment in mice with established cardiac dysfunction. RESULTS: We found nicotinamide riboside chloride improved mitochondrial function and blunted heart failure progression. Similar benefits were observed in wild-type and knockout mice. Boosting NAD+ level improved the function of NAD(H) redox-sensitive SDR (short-chain dehydrogenase/reductase) family proteins. Upregulation of Mrpp2 (mitochondrial ribonuclease P protein 2), a multifunctional SDR protein and a subunit of mitochondrial ribonuclease P, improves mitochondrial DNA transcripts processing and electron transport chain function. Activation of SDRs in the retinol metabolism pathway stimulates RXRα (retinoid X receptor α)/PPARα (proliferator-activated receptor α) signaling and restores mitochondrial oxidative metabolism. Downregulation of Mrpp2 and impaired mitochondrial ribonuclease P were found in human failing hearts, suggesting a shared mechanism of defective mitochondrial biogenesis in mouse and human heart failure. CONCLUSIONS: These findings identify SDR proteins as important regulators of mitochondrial function and molecular targets of NAD+-based therapy. Furthermore, the benefit is observed regardless of Sirt3-mediated mitochondrial protein deacetylation, a widely held mechanism for NAD+-based therapy for heart failure. The data also show that NAD+-based therapy can be useful in pre-existing heart failure.

Incorporating Uremic Solute-mediated Inhibition of OAT1/3 Improves PBPK Prediction of Tenofovir Renal and Systemic Disposition in Patients with Severe Kidney Disease.

BACKGROUND: Dose modification of renally secreted drugs in patients with chronic kidney disease (CKD) has relied on serum creatinine concentration as a biomarker to estimate glomerular filtration (GFR) under the assumption that filtration and secretion decline in parallel. A discrepancy between actual renal clearance and predicted renal clearance based on GFR alone is observed in severe CKD patients with tenofovir, a compound secreted by renal OAT1/3. Uremic solutes that inhibit OAT1/3 may play a role in this divergence. METHODS: To examine the impact of transporter inhibition by uremic solutes on tenofovir renal clearance, we determined the inhibitory potential of uremic solutes hippuric acid, indoxyl sulfate, and p-cresol sulfate. The inhibition parameters (IC50) were incorporated into a previously validated mechanistic kidney model; simulated renal clearance and plasma PK profile were compared to data from clinical studies. RESULTS: Without the incorporation of uremic solute inhibition, the PBPK model failed to capture the observed data with an absolute average fold error (AAFE) > 2. However, when the inhibition of renal uptake transporters and uptake transporters in the slow distribution tissues were included, the AAFE value was within the pre-defined twofold model acceptance criterion, demonstrating successful model extrapolation to CKD patients. CONCLUSION: A PBPK model that incorporates inhibition by uremic solutes has potential to better predict renal clearance and systemic disposition of secreted drugs in patients with CKD. Ongoing research is warranted to determine if the model can be expanded to include other OAT1/3 substrate drugs and to evaluate how these findings can be translated to clinical guidance for drug selection and dose optimization in patients with CKD.

Application of a physiologically based pharmacokinetic model to predict isoniazid disposition during pregnancy.

Pregnancy can increase the risk of latent tuberculosis infection (LTBI) progression to tuberculosis (TB) disease. Isoniazid (INH) is the preferred preventative treatment for LTBI in pregnancy. INH is mainly cleared by N-acetyltransferase 2 (NAT2) but the pharmacokinetics (PK) of INH in different NAT2 phenotypes during pregnancy is not well characterized. To address this knowledge gap, we used physiologically based pharmacokinetic (PBPK) modeling to evaluate NAT2 phenotype-specific effects of pregnancy on INH disposition. A whole-body PBPK model for INH was developed and verified for non-pregnant NAT2 fast (FA), intermediate (IA), and slow (SA) acetylators. Model predictive performance was assessed using a drug-specific model acceptance criterion for mean plasma area under the curve (AUC) and peak plasma concentration (Cmax ), and the absolute average fold error (AAFE) for individual plasma concentrations. The verified model was extended to simulate INH disposition during pregnancy in NAT2 SA, IA, and FA populations. A sensitivity analysis was conducted using the verified PBPK model and known changes in INH disposition during pregnancy to determine whether NAT2 activity changes during pregnancy or other INH clearance pathways are altered. This analysis suggested that NAT2 activity is unchanged while other INH clearance pathways increase by ~80% during pregnancy. The model was applied to explore the effect of pregnancy on INH disposition in two ethnic populations with different NAT2 phenotype distributions and with high TB burden. Our PBPK model can be used to predict INH disposition during pregnancy in diverse populations and expanded to other drugs cleared by NAT2 during pregnancy.

Effects of Pregnancy on Plasma Sphingolipids Using a Metabolomic and Quantitative Analysis Approach.

Changes in the maternal metabolome, and specifically the maternal lipidome, that occur during pregnancy are relatively unknown. The objective of this investigation was to evaluate the effects of pregnancy on sphingolipid levels using metabolomics analysis followed by confirmational, targeted quantitative analysis. We focused on three subclasses of sphingolipids: ceramides, sphingomyelins, and sphingosines. Forty-seven pregnant women aged 18 to 50 years old participated in this study. Blood samples were collected on two study days for metabolomics analysis. The pregnancy samples were collected between 25 and 28 weeks of gestation and the postpartum study day samples were collected ≥3 months postpartum. Each participant served as their own control. These samples were analyzed using a Ultra-performance liquid chromatography/mass spectroscopy/mass spectroscopy (UPLC/MS/MS) assay that yielded semi-quantitative peak area values that were used to compare sphingolipid levels between pregnancy and postpartum. Following this lipidomic analysis, quantitative LC/MS/MS targeted/confirmatory analysis was performed on the same study samples. In the metabolomic analysis, 43 sphingolipid metabolites were identified and their levels were assessed using relative peak area values. These profiled sphingolipids fell into three categories: ceramides, sphingomyelins, and sphingosines. Of the 43 analytes measured, 35 were significantly different during pregnancy (p < 0.05) (including seven ceramides, 26 sphingomyelins, and two sphingosines) and 32 were significantly higher during pregnancy compared to postpartum. Following metabolomics, a separate quantitative analysis was performed and yielded quantified concentration values for 23 different sphingolipids, four of which were also detected in the metabolomics study. Quantitative analysis supported the metabolomics results with 17 of the 23 analytes measured found to be significantly different during pregnancy including 11 ceramides, four sphingomyelins, and two sphingosines. Fourteen of these were significantly higher during pregnancy. Our data suggest an overall increase in plasma sphingolipid concentrations with possible implications in endothelial function, gestational diabetes mellitus (GDM), intrahepatic cholestasis of pregnancy, and fetal development. This study provides evidence for alterations in maternal sphingolipid metabolism during pregnancy.

CYP2D6 Activity Is Correlated with Changes in Plasma Concentrations of Taurocholic Acid during Pregnancy and Postpartum in CYP2D6 Extensive Metabolizers.

Cytochrome P450 2D6 (CYP2D6) is involved in the metabolism of >20% of marketed drugs. CYP2D6 expression and activity exhibit high interindividual variability and is induced during pregnancy. The farnesoid X receptor (FXR) is a transcriptional regulator of CYP2D6 that is activated by bile acids. In pregnancy, elevated plasma bile acid concentrations are associated with maternal and fetal risks. However, modest changes in bile acid concentrations may occur during healthy pregnancy, thereby altering FXR signaling. A previous study demonstrated that hepatic tissue concentrations of bile acids positively correlated with the hepatic mRNA expression of CYP2D6. This study sought to characterize the plasma bile acid metabolome in healthy women (n = 47) during midpregnancy (25-28 weeks gestation) and ≥3 months postpartum and to determine if plasma bile acids correlate with CYP2D6 activity. It is hypothesized that during pregnancy, plasma bile acids would favor less hydrophobic bile acids (cholic acid vs. chenodeoxycholic acid) and that plasma concentrations of cholic acid and its conjugates would positively correlate with the urinary ratio of dextrorphan/dextromethorphan. At 25-28 weeks gestation, taurine-conjugated bile acids comprised 23% of the quantified serum bile acids compared with 7% ≥3 months postpartum. Taurocholic acid positively associated with the urinary ratio of dextrorphan/dextromethorphan, a biomarker of CYP2D6 activity. Collectively, these results confirm that the bile acid plasma metabolome differs between pregnancy and postpartum and provide evidence that taurocholic acid may impact CYP2D6 activity during pregnancy. SIGNIFICANCE STATEMENT: Bile acid homeostasis is altered in pregnancy, and plasma concentrations of taurocholic acid positively correlate with CYP2D6 activity. Differences between plasma and/or tissue concentrations of farnesoid X receptor ligands such as bile acids may contribute to the high interindividual variability in CYP2D6 expression and activity.

The Utility of Mixed Effects Models in the Evaluation of Complex Genomic Traits In Vitro.

In pharmacogenomic studies, the use of human liver microsomes as a model system to evaluate the impact of complex genomic traits (i.e., linkage-disequilibrium patterns, coding, and non-coding variation, etc.) on efficiency of drug metabolism is challenging. To accurately predict the true effect size of genomic traits requires large richly sampled datasets representative of the study population. Moreover, the acquisition of this data can be labor-intensive if the study design or bioanalytical methods are not high throughput, and it is potentially unfeasible if the abundance of sample needed for experiments is limited. To overcome these challenges, we developed a novel strategic approach using non-linear mixed effects models (NLME) to determine enzyme kinetic parameters for individual liver specimens using sparse data. This method can facilitate evaluation of the impact that complex genomic traits have on the metabolism of xenobiotics in vitro when tissue and other resources are limited. In addition to facilitating the accrual of data, it allows for rigorous testing of covariates as sources of kinetic parameter variability. In this in silico study, we present a practical application of such an approach using previously published in vitro cytochrome P450 (CYP) 2D6 data and explore the impact of sparse sampling, and experimental error on known kinetic parameter estimates of CYP2D6 mediated formation of 4-hydroxy-atomoxetine in human liver microsomes. SIGNIFICANCE STATEMENT: This study presents a novel non-linear mixed effects model (NLME)-based framework for evaluating the impact of complex genomic traits on saturable processes described by a Michaelis-Menten kinetics in vitro using sparse data. The utility of this approach extends beyond gene variant associations, including determination of covariate effects on in vitro kinetic parameters and reduced demand for precious experimental material.

Randomized controlled phase IIa clinical trial of safety, pharmacokinetics and pharmacodynamics of tenofovir and tenofovir plus levonorgestrel releasing intravaginal rings used by women in Kenya.

Introduction: Globally, many young women face the overlapping burden of HIV infection and unintended pregnancy. Protection against both may benefit from safe and effective multipurpose prevention technologies. Methods: Healthy women ages 18-34 years, not pregnant, seronegative for HIV and hepatitis B surface antigen, not using hormonal contraception, and at low risk for HIV were randomized 2:2:1 to continuous use of a tenofovir/levonorgestrel (TFV/LNG), TFV, or placebo intravaginal ring (IVR). In addition to assessing genital and systemic safety, we determined TFV concentrations in plasma and cervicovaginal fluid (CVF) and LNG levels in serum using tandem liquid chromatography-mass spectrometry. We further evaluated TFV pharmacodynamics (PD) through ex vivo CVF activity against both human immunodeficiency virus (HIV)-1 and herpes simplex virus (HSV)-2, and LNG PD using cervical mucus quality markers and serum progesterone for ovulation inhibition. Results: Among 312 women screened, 27 were randomized to use one of the following IVRs: TFV/LNG (n = 11); TFV-only (n = 11); or placebo (n = 5). Most screening failures were due to vaginal infections. The median days of IVR use was 68 [interquartile range (IQR), 36-90]. Adverse events (AEs) were distributed similarly among the three arms. There were two non-product related AEs graded >2. No visible genital lesions were observed. Steady state geometric mean amount (ssGMA) of vaginal TFV was comparable in the TFV/LNG and TFV IVR groups, 43,988 ng/swab (95% CI, 31,232, 61,954) and 30337 ng/swab (95% CI, 18,152, 50,702), respectively. Plasma TFV steady state geometric mean concentration (ssGMC) was <10 ng/ml for both TFV IVRs. In vitro, CVF anti-HIV-1 activity showed increased HIV inhibition over baseline following TFV-eluting IVR use, from a median of 7.1% to 84.4% in TFV/LNG, 15.0% to 89.5% in TFV-only, and -27.1% to -20.1% in placebo participants. Similarly, anti-HSV-2 activity in CVF increased >50 fold after use of TFV-containing IVRs. LNG serum ssGMC was 241 pg/ml (95% CI 185, 314) with rapid rise after TFV/LNG IVR insertion and decline 24-hours post-removal (586 pg/ml [95% CI 473, 726] and 87 pg/ml [95% CI 64, 119], respectively). Conclusion: TFV/LNG and TFV-only IVRs were safe and well tolerated among Kenyan women. Pharmacokinetics and markers of protection against HIV-1, HSV-2, and unintended pregnancy suggest the potential for clinical efficacy of the multipurpose TFV/LNG IVR. Clinical Trial Registration: NCT03762382 [https://clinicaltrials.gov/ct2/show/NCT03762382].

Cytosolic Enzymes Generate Cannabinoid Metabolites 7-Carboxycannabidiol and 11-Nor-9-carboxytetrahydrocannabinol.

The cannabinoids cannabidiol (CBD) and delta-9-tetrahydrocannabinol (THC) undergo extensive oxidative metabolism in the liver. Although cytochromes P450 form the primary, pharmacologically active, hydroxylated metabolites of CBD and THC, less is known about the enzymes that generate the major in vivo circulating metabolites of CBD and THC, 7-carboxy-CBD and 11-carboxy-THC, respectively. The purpose of this study was to elucidate the enzymes involved in forming these metabolites. Cofactor dependence experiments with human liver subcellular fractions revealed that 7-carboxy-CBD and 11-carboxy-THC formation is largely dependent on cytosolic NAD+-dependent enzymes, with lesser contributions from NADPH-dependent microsomal enzymes. Experiments with chemical inhibitors provided evidence that 7-carboxy-CBD formation is mainly dependent on aldehyde dehydrogenases and 11-carboxy-THC formation is mediated also in part by aldehyde oxidase. This study is the first to demonstrate the involvement of cytosolic drug-metabolizing enzymes in generating major in vivo metabolites of CBD and THC and addresses a knowledge gap in cannabinoid metabolism.

Impact of Intracellular Lipid Binding Proteins on Endogenous and Xenobiotic Ligand Metabolism and Disposition.

The family of intracellular lipid binding proteins (iLBPs) is comprised of 16 members of structurally related binding proteins that have ubiquitous tissue expression in humans. iLBPs collectively bind diverse essential endogenous lipids and xenobiotics. iLBPs solubilize and traffic lipophilic ligands through the aqueous milieu of the cell. Their expression is correlated with increased rates of ligand uptake into tissues and altered ligand metabolism. The importance of iLBPs in maintaining lipid homeostasis is well established. Fatty acid binding proteins (FABPs) make up the majority of iLBPs and are expressed in major organs relevant to xenobiotic absorption, distribution, and metabolism. FABPs bind a variety of xenobiotics including nonsteroidal anti-inflammatory drugs, psychoactive cannabinoids, benzodiazepines, antinociceptives, and peroxisome proliferators. FABP function is also associated with metabolic disease, making FABPs currently a target for drug development. Yet the potential contribution of FABP binding to distribution of xenobiotics into tissues and the mechanistic impact iLBPs may have on xenobiotic metabolism are largely undefined. This review examines the tissue-specific expression and functions of iLBPs, the ligand binding characteristics of iLBPs, their known endogenous and xenobiotic ligands, methods for measuring ligand binding, and mechanisms of ligand delivery from iLBPs to membranes and enzymes. Current knowledge of the importance of iLBPs in affecting disposition of xenobiotics is collectively described. SIGNIFICANCE STATEMENT: The data reviewed here show that FABPs bind many drugs and suggest that binding of drugs to FABPs in various tissues will affect drug distribution into tissues. The extensive work and findings with endogenous ligands suggest that FABPs may also alter the metabolism and transport of drugs. This review illustrates the potential significance of this understudied area.

The Effects of Pregnancy on Amino Acid Levels and Nitrogen Disposition.

Limited data are available on the effects of pregnancy on the maternal metabolome. Therefore, the objective of this study was to use metabolomics analysis to determine pathways impacted by pregnancy followed by targeted confirmatory analysis to provide more powerful conclusions about metabolic alterations during pregnancy. Forty-seven pregnant women, 18-50 years of age were included in this study, with each subject serving as their own control. Plasma samples were collected between 25 and 28 weeks gestation and again ≥3 months postpartum for metabolomics analysis utilizing an HILIC/UHPLC/MS/MS assay with confirmatory targeted specific concentration analysis for 10 of the significantly altered amino acids utilizing an LC/MS assay. Principle component analysis (PCA) on metabolomics data clearly separated pregnant and postpartum groups and identified outliers in a preliminary assessment. Of the 980 metabolites recorded, 706 were determined to be significantly different between pregnancy and postpartum. Pathway analysis revealed three significantly impacted pathways, arginine biosynthesis (p = 2 × 10-5 and FDR = 1 × 10-3), valine, leucine, and isoleucine metabolism (p = 2 × 10-5 and FDR = 2 × 10-3), and xanthine metabolism (p = 4 × 10-5 and FDR = 4 × 10-3). Of these we focused analysis on arginine biosynthesis and branched-chain amino acid (BCAA) metabolism due to their clinical importance and interconnected roles in amino acid metabolism. In the confirmational analysis, 7 of 10 metabolites were confirmed as significant and all 10 confirmed the direction of change of concentrations observed in the metabolomics analysis. The data support an alteration in urea nitrogen disposition and amino acid metabolism during pregnancy. These changes could also impact endogenous nitric oxide production and contribute to diseases of pregnancy. This study provides evidence for changes in both the ammonia-urea nitrogen and the BCAA metabolism taking place during pregnancy.

Temporal changes in the systemic concentrations of retinoids in pregnant and postpartum women.

Retinoids and vitamin A are essential for multiple biological functions, including vision and immune responses, as well as the development of an embryo during pregnancy. Despite its importance, alterations in retinoid homeostasis during normal human pregnancy are incompletely understood. We aimed to characterize the temporal changes in the systemic retinoid concentrations across pregnancy and postpartum period. Monthly blood samples were collected from twenty healthy pregnant women, and plasma concentrations of retinol, all-trans-retinoic acid (atRA), 13-cis-retinoic acid (13cisRA), and 4-oxo-retinoic acids were measured using liquid chromatography-tandem mass spectrometry. Significant decreases in 13cisRA concentrations over the pregnancy were observed, with rebound increases in retinol and 13cisRA levels after delivery. Of note, atRA concentrations exhibited a unique temporal pattern with levels peaking at mid-pregnancy. While the 4-oxo-atRA concentration was below the limit of quantification, 4-oxo-13cisRA was readily detectable, and its temporal change mimicked that of 13cisRA. The time profiles of atRA and 13cisRA remained similar after correction by albumin levels for plasma volume expansion adjustment. Together, the comprehensive profiling of systemic retinoid concentrations over the course of pregnancy provides insights into pregnancy-mediated changes in retinoid disposition to maintain its homeostasis.