Evaluating the impact of anatomical and physiological variability on human equivalent doses using PBPK models

Addressing human anatomical and physiological variability is a crucial component of human health risk assessment of chemicals. Experts have recommended probabilistic chemical risk assessment paradigms in which distributional adjustment factors are used to account for various sources of uncertainty and variability, including variability in the pharmacokinetic behavior of a given substance in different humans. In practice, convenient assumptions about the distribution forms of adjustment factors and human equivalent doses (HEDs) are often used. Parameters such as tissue volumes and blood flows are likewise often assumed to be lognormally or normally distributed without evaluating empirical data for consistency with these forms. In this work, we performed dosimetric extrapolations using physiologically based pharmacokinetic (PBPK) models for dichloromethane (DCM) and chloroform that incorporate uncertainty and variability to determine if the HEDs associated with such extrapolations are approximately lognormal and how they depend on the underlying distribution shapes chosen to represent model parameters. We accounted for uncertainty and variability in PBPK model parameters by randomly drawing their values from a variety of distribution types. We then performed reverse dosimetry to calculate HEDs based on animal points of departure (PODs) for each set of sampled parameters. Corresponding samples of HEDs were tested to determine the impact of input parameter distributions on their central tendencies, extreme percentiles, and degree of conformance to lognormality. This work demonstrates that the measurable attributes of human variability should be considered more carefully and that generalized assumptions about parameter distribution shapes may lead to inaccurate estimates of extreme percentiles of HEDs.

Sex and the Kidney Drug-Metabolizing Enzymes and Transporters: Are Preclinical Drug Disposition Data Translatable to Humans?

Cross-species differences in drug transport and metabolism are linked to poor translation of preclinical pharmacokinetic and toxicology data to humans, often resulting in the failure of new chemical entities (NCEs) during clinical drug development. Specifically, inaccurate prediction of renal clearance and renal accumulation of NCEs due to differential abundance of enzymes and transporters in kidneys can lead to differences in pharmacokinetics and toxicity between experimental animals and humans. We carried out liquid chromatography-tandem mass spectrometry (LC-MS/MS)-based protein quantification of 78 membrane drug-metabolizing enzymes and transporters (DMETs) in the kidney membrane fractions of humans, rats, and mice for characterization of cross-species and sex-dependent differences. In general, majority of DMET proteins were higher in rodents than in humans. Significant cross-species differences were observed in 30 out of 33 membrane DMET proteins quantified in all three species. Although no significant sex-dependent differences were observed in humans, the abundance of 28 and 46 membrane proteins showed significant sex dependence in rats and mice, respectively. These cross-species and sex-dependent quantitative abundance data are valuable for gaining a mechanistic understanding of drug renal disposition and accumulation. Further, these data can also be integrated into systems pharmacology tools, such as physiologically based pharmacokinetic models, to enhance the interpretation of preclinical pharmacokinetic and toxicological data.

Global Proteomics Analysis of Lysophosphatidic Acid Signaling in PC-3 Human Prostate Cancer Cells: Role of CCN1

Cysteine-rich angiogenic factor 61 (CCN1/Cyr61) is a matricellular protein that is induced and secreted in response to growth factors. Our previous work showed that 18:1-lysophosphatidic acid (LPA), which activates the G protein-coupled receptor LPAR1, induces CCN1 between 2-4 h in PC-3 human prostate cancer cells in a manner than enhances cell-substrate adhesion. While the time course of induction suggests that CCN1 contributes to intermediate events in LPA action, the roles of CCN1 in LPA-mediated signal transduction have not been fully elucidated. This study utilized a comprehensive global proteomics approach to identify proteins up- or down-regulated in response to treatment of PC-3 cells with LPA for three hours, during the time of peak CCN1 levels. In addition, the effects of siRNA-mediated CCN1 knockdown on LPA responses were analyzed. The results show that, in addition to CCN1, LPA increased the levels of multiple proteins. Proteins up-regulated by LPA included metastasis-associated in colon cancer protein 1 (MACC1) and thrombospondin-1 (TSP1/THBS1); both MACC1 and TSP1 regulated cancer cell adhesion and motility. LPA down-regulated thioredoxin interacting protein (TXNIP). CCN1 knockdown suppressed the LPA-induced up-regulation of 30 proteins; these included MACC1 and TSP1, as confirmed by immunoblotting. Gene ontology and STRING analyses revealed multiple pathways impacted by LPA and CCN1. These results indicate that CCN1 contributes to LPA signaling cascades that occur during the intermediate phase after the initial stimulus. The study provides a rationale for the development of interventions to disrupt the LPA-CCN1 axis.

Promiscuity and Quantitative Contribution of UGT2B17 in Drug and Steroid Metabolism Determined by Experimental and Computational Approaches

Uridine 5'-diphospho-glulcuronosyltransferase 2B17 (UGT2B17) is important in the metabolism of steroids and orally administered drugs due to its high interindividual variability. However, the structural basis governing the substrate selectivity or inhibition of UGT2B17 remains poorly understood. This study investigated 76 FDA-approved drugs and 20 steroids known to undergo glucuronidation for their metabolism by UGT2B17. Specifically, we assessed the substrate selectivity for UGT2B17 over other UGT enzymes using recombinant human UGT2B17 (rUGT2B17), human intestinal microsomes, and human liver microsomes. The quantitative contribution of intestinal UGT2B17 in the glucuronidation of these compounds was characterized using intestinal microsomes isolated from UGT2B17 expressors and nonexpressors. In addition, a structure-based pharmacophore model for UGT2B17 substrates was built and validated using the studied pool of substrates and nonsubstrates. The results show that UGT2B17 could metabolize 23 out of 96 compounds from various chemical classes, including alcohols and carboxylic acids, particularly in the intestine. Interestingly, amines were less susceptible to UGT2B17 metabolism, though they could inhibit the enzyme. Three main pharmacophoric features of UGT2B17 substrates include (1) the presence of an accessible -OH or -COOH group near His35 residue, (2) a hydrophobic functional group at ∼4.5-5 Å from feature 1, and (3) an aromatic ring ∼5-7 Å from feature 2. Most of the studied compounds inhibited UGT2B17 activity irrespective of their substrate potential, indicating the possibility of multiple mechanisms. These data suggest that UGT2B17 is promiscuous in substrate selectivity and inhibition and has a high potential to produce significant variability in the absorption and disposition of orally administered drugs.

Effect of Cimetidine on Metformin Pharmacokinetics and Endogenous Metabolite Levels in Rats

Tubular secretion is a primary mechanism along with glomerular filtration for renal elimination of drugs and toxicants into urine. Organic cation transporters (OCTs) and multidrug and toxic extrusion (MATE) transporters facilitate the active secretion of cationic substrates, including drugs such as metformin and endogenous cations. We hypothesized that administration of cimetidine, an Oct/Mate inhibitor, will result in increased plasma levels and decreased renal clearance of metformin and endogenous Oct/Mate substrates in rats. A paired rat pharmacokinetic study was carried out in which metformin (5 mg/kg, intravenous) was administered as an exogenous substrate of Oct/Mate transporters to six Sprague-Dawley rats with and without cimetidine (100 mg/kg, intraperitoneal). When co-administered with cimetidine, metformin area under the curve increased significantly by 3.2-fold, and its renal clearance reduced significantly by 73%. Untargeted metabolomics was performed to investigate the effect of cimetidine on endogenous metabolome in the blood and urine samples. Over 8,000 features (metabolites) were detected in the blood, which were shortlisted using optimized criteria, i.e., a significant increase (P value < 0.05) in metabolite peak intensity in the cimetidine-treated group, reproducible retention time, and quality of chromatogram peak. The metabolite hits were classified into three groups that can potentially distinguish inhibition of i) extra-renal uptake transport or catabolism, ii) renal Octs, and iii) renal efflux transporters or metabolite formation. The metabolomics approach identified novel putative endogenous substrates of cationic transporters that could be tested as potential biomarkers to predict Oct/Mate transporter mediated drug-drug interactions in the preclinical stages. SIGNIFICANCE STATEMENT: Endogenous substrates of renal transporters in animal models could be used as potential biomarkers to predict renal drug-drug interactions in early drug development. Here we demonstrated that cimetidine, an inhibitor of organic cation transporters (Oct/Mate), could alter the pharmacokinetics of metformin and endogenous cationic substrates in rats. Several putative endogenous metabolites of Oct/Mate transporters were identified using metabolomics approach, which could be tested as potential transporter biomarkers to predict renal drug-drug interaction of Oct/Mate substrates.

Characterization of Gla proteoforms and non-Gla peptides of gamma carboxylated proteins: Application to quantification of prothrombin proteoforms in human plasma

Gamma (γ) carboxylation is an essential post-translational modification in vitamin K-dependent proteins (VKDPs), involved in maintaining critical biological homeostasis. Alterations in the abundance or activity of these proteins have pharmacological and pathological consequences. Importantly, low levels of fully γ-carboxylated clotting factors increase plasma des-γ-carboxy precursors resulting in little or no biological activity. Therefore, it is important to characterize the levels of γ-carboxylation that reflect the active state of these proteins. The conventional enzyme-linked immunosorbent assay for protein induced by vitamin K absence or antagonist II (PIVKA-II) quantification uses an antibody that is not applicable to distinguish different γ-carboxylation states. Liquid chromatography-mass spectrometry (LC-MS) approaches have been utilized to distinguish different γ-carboxylated proteoforms, however, these attempts were impeded by poor sensitivity due to spontaneous neutral loss of CO2 and simultaneous cleavage of the backbone bond in the collision cell. In this study, we utilized an alkaline mobile phase in combination with polarity switching (positive and negative ionization modes) to simultaneously identify and quantify γ-carboxylated VKDPs. The method was applied to compare Gla proteomics of prothrombin (FII) in 10 μL plasma samples of healthy control and warfarin-treated adults. We also identified surrogate non-Gla peptides for seven other VKDPs to quantify total (active plus inactive) protein levels. The total protein approach (TPA) was used to quantify absolute levels of the VKDPs in human plasma.

The Plasma Membrane Monoamine Transporter is Highly Expressed in Neuroblastoma and Functions as an mIBG Transporter

Neuroblastoma (NB) is a pediatric cancer with low survival rates in high-risk patients. 131I-mIBG has emerged as a promising therapy for high-risk NB and kills tumor cells by radiation. Consequently, 131I-mIBG tumor uptake and retention are major determinants for its therapeutic efficacy. mIBG enters NB cells through the norepinephrine transporter (NET), and accumulates in mitochondria through unknown mechanisms. Here we evaluated the expression of monoamine and organic cation transporters in high-risk NB tumors and explored their relationship with MYCN amplification and patient survival. We found that NB mainly expresses NET, the plasma membrane monoamine transporter (PMAT), and the vesicular membrane monoamine transporter 1/2 (VMAT1/2), and that the expression of these transporters is significantly reduced in MYCN-amplified tumor samples. PMAT expression is the highest and correlates with overall survival in high-risk NB patients without MYCN amplification. Immunostaining showed that PMAT resides intracellularly in NB cells and co-localizes with mitochondria. Using cells expressing PMAT, mIBG was identified as a PMAT substrate. In mitochondria isolated from NB cell lines, mIBG uptake was reduced by ∼50% by a PMAT inhibitor. Together, our data suggest that PMAT is a previously unrecognized transporter highly expressed in NB and could impact intracellular transport and therapeutic response to 131I-mIBG. SIGNIFICANCE STATEMENT: This study identified that plasma membrane monoamine transporter (PMAT) is a novel transporter highly expressed in neuroblastoma and its expression level is associated with overall survival rate in high-risk patients without MYCN amplification. PMAT is expressed intracellularly in neuroblastoma cells, transports meta-iodobenzylguanidine (mIBG) and thus could impact tumor retention and response to 131I-mIBG therapy. These findings have important clinical implications as PMAT could represent a novel molecular marker to help inform disease prognosis and predict response to 131I-mIBG therapy.

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

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

Effect of probenecid on blood levels and renal elimination of furosemide and endogenous compounds in rats: Discovery of putative organic anion transporter biomarkers

Transporter-mediated drug-drug interactions (DDIs) are assessed using probe drugs and in vitro and in vivo models during drug development. The utility of endogenous metabolites as transporter biomarkers is emerging for prediction of DDIs during early phases of clinical trials. Endogenous metabolites such as pyridoxic acid and kynurenic acid have shown potential to predict DDIs mediated by organic anion transporters (OAT1 and OAT3). However, these metabolites have not been assessed in rats as potential transporter biomarkers. We carried out a rat pharmacokinetic DDI study using probenecid and furosemide as OAT inhibitor and substrate, respectively. Probenecid administration led to a 3.8-fold increase in the blood concentrations and a 3-fold decrease in renal clearance of furosemide. High inter-individual and intra-day variability in pyridoxic acid and kynurenic acid, and no or moderate effect of probenecid administration on these metabolites suggest their limited utility for prediction of Oat-mediated DDI in rats. Therefore, rat blood and urine samples were further analysed using untargeted metabolomics. Twenty-one m/z features (out of >8000 detected features) were identified as putative biomarkers of rat Oat1 and Oat3 using a robust biomarker qualification approach. These m/z features belong to metabolic pathways such as fatty acid analogues, peptides, prostaglandin analogues, bile acid derivatives, flavonoids, phytoconstituents, and steroids, and can be used as a panel to decrease variability caused by processes other than Oats. When validated, these putative biomarkers will be useful in predicting DDIs caused by Oats in rats.

Analysis of the interplay of physiological response to food intake and drug properties in food-drug interactions

The effect of food on oral drug absorption is determined by the complex interplay among gut physiological factors and drug properties. The currently used dissolution testing and classification systems (biopharmaceutics classification system, BCS or biopharmaceutics drug disposition classification system, BDDCS) do not account for dynamic changes in gastrointestinal physiology caused by food intake. This study aimed to identify key drug properties that influence food effect (FE) using supervised machine learning approaches. The analysis showed that drugs with high logP, dose number, and extraction ratio have a higher probability of positive FE, while drugs with low permeability and high efflux saturation index have a greater likelihood of negative FE. Weakly acidic drugs also showed a greater probability of positive FE, particularly at pKa >4.3. The importance of drug properties in predicting FE was ranked as logP, dose number, extraction ratio, pKa, and permeability. The accuracy of FE prediction using the models was compared with BCS and extended clearance classification system (ECCS). Overall, the likelihood or magnitude of FE depends on physiological changes to food intake such as altered bile acid secretion rate, intestinal metabolism, transport kinetics, and gastric emptying time, which should be considered along with drug properties (e.g., solubility, logP, and ionization) in predicting FE of orally administered drugs.

Impact of CYP2C:TG Haplotype on CYP2C19 Substrates Clearance In Vivo, Protein Content, and In Vitro Activity

A novel haplotype composed of two non-coding variants, CYP2C18 NM_000772.3:c.*31T (rs2860840) and NM_000772.2:c.819+2182G (rs11188059), referred to as "CYP2C:TG," was recently associated with ultrarapid metabolism of various CYP2C19 substrates. As the underlying mechanism and clinical relevance of this effect remain uncertain, we analyzed existing in vivo and in vitro data to determine the magnitude of the CYP2C:TG haplotype effect. We assessed variability in pharmacokinetics of CYP2C19 substrates, including citalopram, sertraline, voriconazole, omeprazole, pantoprazole, and rabeprazole in 222 healthy volunteers receiving one of these six drugs. We also determined its impact on CYP2C8, CYP2C9, CYP2C18, and CYP2C19 protein abundance in 135 human liver tissue samples, and on CYP2C18/CYP2C19 activity in vitro using N-desmethyl atomoxetine formation. No effects were observed according to CYP2C:TG haplotype or to CYP2C19*1+TG alleles (i.e., CYP2C19 alleles containing the CYP2C:TG haplotype). In contrast, CYP2C19 intermediate (e.g., CYP2C19*1/*2) and poor metabolizers (e.g., CYP2C19*2/*2) showed significantly higher exposure in vivo, lower CYP2C19 protein abundance in human liver microsomes, and lower activity in vitro compared with normal, rapid (i.e., CYP2C19*1/*17), and ultrarapid metabolizers (i.e., CYP2C19*17/*17). Moreover, a tendency toward lower exposure was observed in ultrarapid metabolizers compared with rapid metabolizers and normal metabolizers. Furthermore, when the CYP2C19*17 allele was present, CYP2C18 protein abundance was increased suggesting that genetic variation in CYP2C19 may be relevant to the overall metabolism of certain drugs by regulating not only its expression levels, but also those of CYP2C18. Considering all available data, we conclude that there is insufficient evidence supporting clinical CYP2C:TG testing to inform drug therapy.

Dissecting Parameters Contributing to the Underprediction of Aldehyde Oxidase-Mediated Metabolic Clearance of Drugs

We investigated the effect of variability and instability in aldehyde oxidase (AO) content and activity on the scaling of in vitro metabolism data. AO content and activity in human liver cytosol (HLC) and five recombinant human AO preparations (rAO) were determined using targeted proteomics and carbazeran oxidation assay, respectively. AO content was highly variable as indicated by the relative expression factor (REF; i.e., HLC to rAO content) ranging from 0.001 to 1.7 across different in vitro systems. The activity of AO in HLC degrades at a 10-fold higher rate in the presence of the substrate as compared with the activity performed after preincubation without substrate. To scale the metabolic activity from rAO to HLC, a protein-normalized activity factor (pnAF) was proposed wherein the activity was corrected by AO content, which revealed up to sixfold higher AO activity in HLC versus rAO systems. A similar value of pnAF was observed for another substrate, ripasudil. Physiologically based pharmacokinetic (PBPK) modeling revealed a significant additional clearance (CL; 66%), which allowed for the successful prediction of in vivo CL of four other substrates, i.e., O-benzyl guanine, BIBX1382, zaleplon, and zoniporide. For carbazeran, the metabolite identification study showed that the direct glucuronidation may be contributing to around 12% elimination. Taken together, this study identified differential protein content, instability of in vitro activity, role of additional AO clearance, and unaccounted metabolic pathways as plausible reasons for the underprediction of AO-mediated drug metabolism. Consideration of these factors and integration of REF and pnAF in PBPK models will allow better prediction of AO metabolism. SIGNIFICANCE STATEMENT: This study elucidated the plausible reasons for the underprediction of aldehyde oxidase (AO)-mediated drug metabolism and provided recommendations to address them. It demonstrated that integrating protein content and activity differences and accounting for the loss of AO activity, as well as consideration of extrahepatic clearance and additional pathways, would improve the in vitro to in vivo extrapolation of AO-mediated drug metabolism using physiologically based pharmacokinetic modeling.

High-Throughput Assay of Cytochrome P450-Dependent Drug Demethylation Reactions and Its Use to Re-Evaluate the Pathways of Ketamine Metabolism

In a search for a reliable, inexpensive, and versatile technique for high-throughput kinetic assays of drug metabolism, we elected to rehire an old-school approach based on the determination of formaldehyde (FA) formed in cytochrome P450-dependent demethylation reactions. After evaluating several fluorometric techniques for FA detection, we chose the method based on the Hantzsch reaction with acetoacetanilide as the most sensitive, robust, and adaptable to high-throughput implementation. Here we provide a detailed protocol for using our new technique for automatized assays of cytochrome P450-dependent drug demethylations and discuss its applicability for high-throughput scanning of drug metabolism pathways in the human liver. To probe our method further, we applied it to re-evaluating the pathways of metabolism of ketamine, a dissociative anesthetic and potent antidepressant increasingly used in the treatment of alcohol withdrawal syndrome. Probing the kinetic parameters of ketamine demethylation by ten major cytochrome P450 (CYP) enzymes, we demonstrate that in addition to CYP2B6 and CYP3A enzymes, which were initially recognized as the primary metabolizers of ketamine, an important role is also played by CYP2C19 and CYP2D6. At the same time, the involvement of CYP2C9 suggested in the previous reports was deemed insignificant.

Involvement of the p38/MK2 Pathway in MCLR Hepatotoxicity Revealed through MAPK Pharmacological Inhibition and Phosphoproteomics in HepaRG Cells

Microcystin-leucine arginine (MCLR) is one of the most common and toxic microcystin variants, a class of cyanotoxins produced by cyanobacteria. A major molecular mechanism for MCLR-elicited liver toxicity involves the dysregulation of protein phosphorylation through protein phosphatase (PP) inhibition and mitogen-activated protein kinase (MAPK) modulation. In this study, specific pharmacological MAPK inhibitors were used in HepaRG cells to examine the pathways associated with MCLR cytotoxicity. SB203580 (SB), a p38 inhibitor, rescued HepaRG cell viability, whereas treatment with SP600125 (JNK inhibitor), MK2206 (AKT inhibitor), or N-acetylcysteine (reactive oxygen species scavenger) did not. Phosphoproteomic analysis revealed that phosphosites-which were altered by the addition of SB compared to MCLR treatment alone-included proteins involved in RNA processing, cytoskeletal stability, DNA damage response, protein degradation, and cell death. A closer analysis of specific proteins in some of these pathways indicated that SB reversed the MCLR-mediated phosphorylation of the necroptosis-associated proteins, the mixed lineage kinase domain-like protein (MLKL), receptor-interacting serine/threonine kinase 1 (RIP1), DNA damage response proteins, ataxia telangiectasia and Rad3-related kinase (ATR), and checkpoint kinase 1 (CHK1). Overall, these data implicate p38/MK2, DNA damage, and necroptosis in MCLR-mediated hepatotoxicity, and suggest these pathways may be targets for prevention prior to, or treatment after, MCLR toxicity.

Characterization of xenobiotic and steroid disposition potential of human placental tissue and cell lines (BeWo, JEG-3, JAR, and HTR-8/SVneo) by quantitative proteomics

The placenta is a fetal organ that performs critical functions to maintain pregnancy and support fetal development, including metabolism and transport of xenobiotics and steroids between the maternal-fetal unit. In vitro placenta models are used to study xenobiotic and steroid disposition, but how well these models recapitulate the human placenta is not well understood. We first characterized the abundance of proteins involved in xenobiotic and steroid disposition in human placental tissue. In pooled human placenta, the following xenobiotic and steroid disposition proteins were detected (highest to lowest), i) enzymes: glutathione S-transferase P, carbonyl reductase 1, aldo-keto reductase 1B1, hydroxysteroid dehydrogenases (HSD3B1 and HSD11B1), aromatase, epoxide hydrolase 1 (EPHX1) and steryl-sulfatase, and ii) transporters: monocarboxylate transporters (MCT1 and 4), organic anion transporting polypeptide 2B1, organic anion transporter 4, and breast cancer resistance protein (BCRP). Then, the tissue proteomics data was compared with four placental cell lines (BeWo, JEG-3, JAR, and HTR-8/SVneo). The differential global proteomics analysis revealed that the tissue and cell lines shared 1420 cytosolic and 1186 membrane proteins. Although extravillous trophoblast and cytotrophoblast marker proteins were detected in all cell lines, only BeWo and JEG-3 cells expressed the syncytiotrophoblast marker, chorionic somatomammotropin hormone 1. BeWo and JEG-3 cells expressed most target proteins including aromatase, HSDs, EPHX1, MCT1, and BCRP. JEG-3 cells treated with phthalates commonly detected in human biofluids, showed dysregulation of steroid pathways. The data presented here showed that BeWo and JEG-3 cells are closer to the placental tissue for studying xenobiotic and steroid disposition. Significance Statement This is the first study to compare proteomics data of human placental tissue and cell lines (BeWo, JAR, JEG-3, and HTR-8/SVneo). The placental cell line and tissue proteomes are vastly different, but BeWo and JEG-3 cells showed greater resemblance to the tissue in the expression of xenobiotic and steroid disposition proteins. These data will assist researchers to select an optimum cell model for mechanistic investigations on xenobiotic and steroid disposition in the placenta.

Intestinal metabolism of diclofenac by polymorphic UGT2B17 correlates with its highly variable pharmacokinetics and safety across populations

Although the United States and Europe have shifted to the prescription use of diclofenac due to several serious incidences of cardiotoxicity, it is one of the most commonly used over-the-counter (OTC) medicine in major parts of the world. We elucidated the quantitative and tissue-specific contribution of uridine diphosphate-glucuronosyltransferases 17 (UGT2B17) in diclofenac metabolism and pharmacokinetics. UGT2B17 is one of most deleted genes in humans with the gene deletion frequency ranging from ~20% in Caucasians to 90% in Japanese. The human intestinal and liver microsomes isolated from the high-UGT2B17 expressing individuals showed 21- and 4-fold greater rate of diclofenac glucuronide (DG) formation than in the null-UGT2B17 carriers, respectively. The greater contribution of intestinal UGT2B17 was confirmed by a strong correlation (R= 0.78, p<0.001) between UGT2B17 abundance and DG formation in individual intestinal microsomes (n=14). However, because UGT2B17 is a minor UGT isoform in the liver, DG formation rate correlated better with the expression of UGT2B7. The proteomics-informed physiologically-based pharmacokinetic model explains the reported higher exposure of diclofenac in females consistent with ~3-fold lower expression of UGT2B17. Similarly, our predictions also corroborate with the reported higher exposure and lower standard clinical dose of diclofenac in Japanese population. Therefore, UGT2B17 mediated metabolism of oral diclofenac is a concern, especially in the developing countries where it is still used as an OTC drug. The ontogeny data of UGTs in human hepatocytes can be utilized in developing PBPK models for predicting pharmacokinetics in pediatric population.

Relevance of Human Aldoketoreductases and Microbial β-Glucuronidases in Testosterone Disposition

Testosterone exhibits high variability in pharmacokinetics and glucuronidation after oral administration. Although testosterone metabolism has been studied for decades, the impact of UGT2B17 gene deletion and the role of gut bacterial β-glucuronidases on its disposition are not well characterized. We first performed an exploratory study to investigate the effect of UGT2B17 gene deletion on the global liver proteome, which revealed significant increases in proteins from multiple biological pathways. The most upregulated liver proteins were aldoketoreductases [AKR1D1, AKR1C4, AKR7A3, AKR1A1, and 7-dehydrocholesterol reductase (DHCR7)] and alcohol or aldehyde dehydrogenases (ADH6, ADH1C, ALDH1A1, ALDH9A1, and ALDH5A). In vitro assays revealed that AKR1D1 and AKR1C4 inactivate testosterone to 5β-dihydrotestosterone (5β-DHT) and 3α,5β-tetrahydrotestosterone (3α,5β-THT), respectively. These metabolites also appeared in human hepatocytes treated with testosterone and in human serum collected after oral testosterone dosing in men. Our data also suggest that 5β-DHT and 3α, 5β-THT are then eliminated through glucuronidation by UGT2B7 in UGT2B17 deletion individuals. Second, we evaluated the potential reactivation of testosterone glucuronide (TG) after its secretion into the intestinal lumen. Incubation of TG with purified gut microbial β-glucuronidase enzymes and with human fecal extracts confirmed testosterone reactivation into testosterone by gut bacterial enzymes. Both testosterone metabolic switching and variable testosterone activation by gut microbial enzymes are important mechanisms for explaining the disposition of orally administered testosterone and appear essential to unraveling the molecular mechanisms underlying UGT2B17-associated pathophysiological conditions. SIGNIFICANCE STATEMENT: This study investigated the association of UGT2B17 gene deletion and gut bacterial β-glucuronidases with testosterone disposition in vitro. The experiments revealed upregulation of AKR1D1 and AKR1C4 in UGT2B17 deletion individuals, and the role of these enzymes to inactivate testosterone to 5β-dihydrotestosterone and 3α, 5β-tetrahydrotestosterone, respectively. Key gut bacterial species responsible for testosterone glucuronide activation were identified. These data are important for explaining the disposition of exogenously administered testosterone and appear essential to unraveling the molecular mechanisms underlying UGT2B17-associated pathophysiological conditions.

Quantitative Characterization of Clinically Relevant Drug-Metabolizing Enzymes and Transporters in Rat Liver and Intestinal Segments for Applications in PBPK Modeling

Rats are extensively used as a preclinical model for assessing drug pharmacokinetics (PK) and tissue distribution; however, successful translation of the rat data requires information on the differences in drug metabolism and transport mechanisms between rats and humans. To partly fill this knowledge gap, we quantified clinically relevant drug-metabolizing enzymes and transporters (DMETs) in the liver and different intestinal segments of Sprague-Dawley rats. The levels of DMET proteins in rats were quantified using the global proteomics-based total protein approach (TPA) and targeted proteomics. The abundance of the major DMET proteins was largely comparable using quantitative global and targeted proteomics. However, global proteomics-based TPA was able to detect and quantify a comprehensive list of 66 DMET proteins in the liver and 37 DMET proteins in the intestinal segments of SD rats without the need for peptide standards. Cytochrome P450 (Cyp) and UDP-glycosyltransferase (Ugt) enzymes were mainly detected in the liver with the abundance ranging from 8 to 6502 and 74 to 2558 pmol/g tissue. P-gp abundance was higher in the intestine (124.1 pmol/g) as compared to that in the liver (26.6 pmol/g) using the targeted analysis. Breast cancer resistance protein (Bcrp) was most abundant in the intestinal segments, whereas organic anion transporting polypeptides (Oatp) 1a1, 1a4, 1b2, and 2a1 and multidrug resistance proteins (Mrp) 2 and 6 were predominantly detected in the liver. To demonstrate the utility of these data, we modeled digoxin PK by integrating protein abundance of P-gp and Cyp3a2 into a physiologically based PK (PBPK) model constructed using PK-Sim software. The model was able to reliably predict the systemic as well as tissue concentrations of digoxin in rats. These findings suggest that proteomics-informed PBPK models in preclinical species can allow mechanistic PK predictions in animal models including tissue drug concentrations.

Quantification of Accurate Composition and Total Abundance of Homologous Proteins by Conserved-Plus-Surrogate Peptide Approach: Quantification of UDP Glucuronosyltransferases in Human Tissues

Characterization of accurate compositions and total abundance of homologous drug-metabolizing enzymes, such as UDP glucuronosyltransferases (UGTs), is important for predicting the fractional contribution of individual isoforms involved in the metabolism of a drug for applications in physiologically based pharmacokinetic (PBPK) modeling. Conventional targeted proteomics utilizes surrogate peptides, which often results in high technical and interlaboratory variability due to peptide-specific digestion leading to data inconsistencies. To address this problem, we developed a novel conserved-plus-surrogate peptide (CPSP) approach for determining the accurate compositions and total or cumulative abundance of homologous UGTs in commercially available pooled human liver microsomes (HLM), human intestinal microsomes (HIM), human kidney microsomes (HKM), and human liver S9 (HLS9) fraction. The relative percent composition of UGT1A and UGT2B isoforms in the human liver was 35:5:36:11:13 for UGT1A1:1A3:1A4:1A6:1A9 and 20:32:22:21:5 for UGT2B4:2B7:2B10:2B15:2B17. The human kidney and intestine also showed unique compositions of UGT1As and UGT2Bs. The reproducibility of the approach was validated by assessing correlations of UGT compositions between HLM and HLS9 (R²> 0.91). The analysis of the conserved peptides also provided the abundance for individual UGT isoforms included in this investigation as well as the total abundance (pmol/mg protein) of UGT1As and UGT2Bs across tissues, i.e., 268 and 342 (HLM), 21 and 92 (HIM), and 138 and 99 (HKM), respectively. The CPSP approach could be used for applications in the in-vitro-to-in-vivo extrapolation of drug metabolism and PBPK modeling. SIGNIFICANCE STATEMENT: We quantified the absolute compositions and total abundance of UDP glucuronosyltransferases (UGTs) in pooled human liver, intestine, and kidney microsomes using a novel conserved-plus-surrogate peptide (CPSP) approach. The CPSP approach addresses the surrogate peptide-specific variability in the determination of the absolute composition of UGTs. The data presented in this manuscript are applicable for the estimation of the fraction metabolized by individual UGTs towards better in vitro-to-in vivo extrapolation of UGT-mediated drug metabolism.

Interplay of Breast Cancer Resistance Protein (Bcrp/Abcg2), Sex, and Fed State in Oral Pharmacokinetic Variability of Furosemide in Rats

Poor and variable oral bioavailability of furosemide (FUR) presents critical challenges in pharmacotherapy. We investigated the interplay of breast cancer resistance protein (Bcrp)-mediated transport, sex, and fed state on FUR pharmacokinetics (PK) in rats. A crossover PK study of FUR (5 mg/kg, oral) was performed in Sprague-Dawley rats (3 males and 3 females), alone or with a Bcrp inhibitor, novobiocin (NOV) (20 mg/kg, oral), in both fed and fasted states. Co-administration of NOV significantly increased FUR extent (AUC) and rate (C_max) of exposure by more than two-fold, which indicates efficient Bcrp inhibition in the intestine. The female rats showed two-fold higher AUC and C_max, and two-fold lower renal clearance of FUR compared to the male rats. The latter was correlated with higher renal abundance of Bcrp and organic anion transporters (Oats) in the male rats compared to age-matched female rats. These findings suggest that the PK of Bcrp and/or Oat substrates could be sex-dependent in rats. Moreover, allometric scaling of rat PK and toxicological data of Bcrp substrates should consider species and sex differences in Bcrp and Oat abundance in the kidney. Considering that Bcrp is abundant in the intestine of rats and humans, a prospective clinical study is warranted to evaluate the effect of Bcrp inhibition on FUR PK. The potential confounding effect of the Bcrp transporter should be considered when FUR is used as a clinical probe of renal organic anion transporter-mediated drug-drug interactions. Unlike human data, no food-effect was observed on FUR PK in rats.

An optimized proteomics-based approach to estimate blood contamination and cellular heterogeneity of frozen placental tissue

INTRODUCTION

Archived human placental tissue specimens are vital for studying placenta pathophysiology and toxicology. Proteomics analysis of placental tissue provides mechanistic and translational information, but the highly perfused and heterogenous nature of the placenta creates confounding technical variability. In this study, we developed an optimized proteomics-based approach to address the technical variability of proteomics data by normalizing blood contamination and cellular heterogeneity of archived placenta samples.

METHODS

Placenta samples (n = 99) were homogenized, digested using trypsin, and analyzed by liquid chromatography mass-spectrometry. Label-free quantification (LFQ) intensities of the proteins were analyzed for their correlation with blood (albumin) and placenta (aromatase) markers. Proteins that positively correlated with albumin and negatively correlated with aromatase or vice versa were considered blood and placental proteins, respectively. Next, the cellular heterogeneity of individual placenta samples was evaluated by quantifying specific cellular markers of cytotrophoblasts, syncytiotrophoblasts, extravillous trophoblasts, fibroblasts, Hofbauer cells, and decidual cells.

RESULTS

We found that placental proteins were contaminated by 41 to 85% blood proteins. Analysis of cellular markers confirmed syncytiotrophoblasts as the major cell type in placenta (i.e., 41 ± 9% of all cell types). Two samples showed distinct cell compositions with higher levels of the extravillous trophoblasts and decidual cells.

DISCUSSION

In summary, the optimized proteomics-based approach to estimate blood contamination and cellular heterogeneity of placental tissues has the potential to address technical variability in placenta proteomics analysis, which can be extended to other highly perfused and heterogenous tissues.

Cardiac Disease Alters Myocardial Tissue Levels of Epoxyeicosatrienoic Acids and Key Proteins Involved in Their Biosynthesis and Degradation

CYP2J2 is the main epoxygenase in the heart that is responsible for oxidizing arachidonic acid to cis-epoxyeicosatrienoic acids (EETs). Once formed, EETs can then be hydrolyzed by soluble epoxide hydrolase (sEH, encoded by EPHX2) or re-esterified back to the membrane. EETs have several cardioprotective properties and higher levels are usually associated with better cardiac outcomes/prognosis. This study investigates how cardiovascular disease (CVD) can influence total EET levels by altering protein expression and activity of enzymes involved in their biosynthesis and degradation. Diseased ventricular cardiac tissues were collected from patients receiving Left Ventricular Assist Device (LVAD) or heart transplants and compared to ventricular tissue from controls free of CVD. EETs, and enzymes involved in EETs biosynthesis and degradation, were measured using mass spectrometric assays. Terfenadine hydroxylation was used to probe CYP2J2 activity. Significantly higher cis- and trans-EET levels were observed in control cardiac tissue (n = 17) relative to diseased tissue (n = 24). Control cardiac tissue had higher CYP2J2 protein levels, which resulted in higher rate of terfenadine hydroxylation, compared to diseased cardiac tissues. In addition, levels of both NADPH-Cytochrome P450 oxidoreductase (POR) and sEH proteins were significantly higher in control versus diseased cardiac tissue. Overall, alterations in protein and activity of enzymes involved in the biosynthesis and degradation of EETs provide a mechanistic understanding for decreased EET levels in diseased tissues.

Dimethandrolone (DMA), a Potential Male Contraceptive Pill, is Primarily Metabolized by the Highly Polymorphic UGT2B17 Enzyme in Human Intestine and Liver

Dimethandrolone undecanoate (DMAU), an oral investigational male hormonal contraceptive is a prodrug, which is rapidly converted to its active metabolite, dimethandrolone (DMA). Poor and variable oral bioavailability of DMA after DMAU dosing is a critical challenge to develop it as an oral drug. The objective of our study was to elucidate the mechanisms of variable pharmacokinetics (PK) of DMA. We first identified DMA metabolites formed in vitro and in vivo in human hepatocyte incubation and serum samples following oral DMAU administration in men, respectively. The metabolite identification study revealed two metabolites, DMA-glucuronide (DMA-G; major) and the androstenedione analogue of DMA (minor) in the hepatocyte incubations. After oral DMAU administration, only DMA-G was detected in serum, which was >100-fold compared to DMA levels, supporting glucuronidation as the major elimination mechanism for DMA. Next, thirteen clinically- relevant UDP-glucuronosyltransferase (UGT) enzymes were tested for their involvement in DMA-G formation, which revealed a major role of UGT2B17 isoform with a smaller contribution of UGT1A9 in DMA-G formation. These data were confirmed by dramatically higher DMA glucuronidation rates (>200- and >7-fold) in the high versus the null UGT2B17-expressing human intestinal and liver microsomes, respectively. Since human UGT2B17 is a highly variable enzyme with a 20-80% gene deletion frequency, the in vitro data suggest a major role of UGT2B17 polymorphism on the first-pass metabolism of DMA. Further, considering DMA is a selective and sensitive UGT2B17 substrate, it could be used as a clinical probe of UGT2B17 activity. Significance Statement Dimethandrolone (DMA) is an active metabolite of dimethandrolone undecanoate (DMAU), an investigational male hormonal contraceptive. Previous studies have indicated poor and inconsistent bioavailability of DMAU following oral administration. Our study found that UGT2B17-mediated high intestinal first-pass metabolism is the key mechanism of variable DMA bioavailability.

Can analysis of serum androgens aid in the diagnosis of polycystic ovary syndrome (PCOS) in adolescents?

INTRODUCTION

Polycystic ovary syndrome (PCOS) is a female metabolic disorder that is characterized by ovulatory dysfunction, elevated serum androgen concentrations, and polycystic ovarian morphology (PCOM). However, diagnosis of PCOS in adolescents is challenging.

AREAS COVERED

The mechanisms of PCOS pathophysiology are discussed that include: i) dysregulation of the levels of steroidal enzymes ii) abnormalities in the secretion of gonadotropin releasing hormone, luteinizing hormone, and follicle stimulating hormone , and iii) abnormalities in ovarian Thecal and Granulosa cell function. Current clinical diagnosis protocols for PCOS in women are covered. The challenges in diagnosis of PCOS particularly in adolescents are highlightedWe highlighted an important unmet need for an accurate serum test for the early diagnosis of adolescent girls with PCOS.

EXPERT OPINION

Steroid metabolite profiling that captures hyperandrogenism has shown some early promise to serve as a biomarker for early diagnosis of PCOS in women, something that would be especially useful in adolescents.

Differential proteomics analysis of JEG-3 and JAR placental cell models and the effect of androgen treatment

The placenta is a vital fetal organ that plays an important role in maintaining fetal sex hormone homeostasis. Xenobiotics can alter placental sex-steroidogenic enzymes and transporters, including enzymes such as aromatase (CYP19A1) and the hydroxysteroid dehydrogenases (HSDs) but studying how compounds disrupt in vivo placental metabolism is complex. Utilizing high-throughput in vitro models is critical to predict the disruption of placental sex-steroidogenic enzymes and transporters, particularly by drug candidates in the early stages of drug discovery. JAR and JEG-3 cells are the most common, simple, and cost-effective placental cell models that are capable of high-throughput screening, but how well they express the sex-steroidogenic enzymes and transporters is not well known. Here, we compared the proteomes of JAR and JEG-3 cells in the presence and absence of physiologically relevant concentrations of dehydroepiandrosterone (DHEA, 8 µM) and testosterone (15 nM) to aid the characterization of sex-steroidogenic enzymes and transporters in these cell models. Global proteomics analysis detected 2931 and 3449 proteins in JAR cells and JEG-3 cells, respectively. However, dramatic differences in sex-steroidogenic enzymes and transporters were observed between these cells. In particular, the basal expression of steroid sulfatase (STS), HSD17B1, and HSD17B7 were unique to JEG-3 cells. JEG-3 cells also showed significantly higher protein levels of aldo-keto reductase (AKR) 1A1 and AKR1B1, while JAR cells showed significantly higher levels of HSD17B4 and HSDB12. Aldehyde dehydrogenase (ALDH) 3A2 and HSD17B11 enzymes as well as the transporters sterol O-acyltransferase (SOAT) 1 and ATP binding cassette subfamily G2 (ABCG2) were comparable between the cell lines, whereas sulfotransferases (SULTs) were uniquely present within JAR cells. Androgen treatments significantly lowered HSD17B11, HSD17B4, HSD17B12, and ALDH3A2 levels in JAR cells. DHEA treatment significantly raised the level of HSD17B1 by 51 % in JEG-3 cells, whereas CYP19A1 was increased to significant levels in both JAR and JEG-3 cells after androgen treatments. The proteomics data were supported by a complementary targeted metabolomics analysis of culture media in the DHEA (8 µM) and testosterone (15 nM) treated groups. This study has indicated that untreated JEG-3 cells express more sex-steroidogenic enzymes and transporters. Nevertheless, JEG-3 and JAR cells are unique and their respective proteomics data can be used to select the best model depending on the hypothesis.

Quantitative Proteomics in Translational Absorption, Distribution, Metabolism, and Excretion and Precision Medicine

A reliable translation of in vitro and preclinical data on drug absorption, distribution, metabolism, and excretion (ADME) to humans is important for safe and effective drug development. Precision medicine that is expected to provide the right clinical dose for the right patient at the right time requires a comprehensive understanding of population factors affecting drug disposition and response. Characterization of drug-metabolizing enzymes and transporters for the protein abundance and their interindividual as well as differential tissue and cross-species variabilities is important for translational ADME and precision medicine. This review first provides a brief overview of quantitative proteomics principles including liquid chromatography-tandem mass spectrometry tools, data acquisition approaches, proteomics sample preparation techniques, and quality controls for ensuring rigor and reproducibility in protein quantification data. Then, potential applications of quantitative proteomics in the translation of in vitro and preclinical data as well as prediction of interindividual variability are discussed in detail with tabulated examples. The applications of quantitative proteomics data in physiologically based pharmacokinetic modeling for ADME prediction are discussed with representative case examples. Finally, various considerations for reliable quantitative proteomics analysis for translational ADME and precision medicine and the future directions are discussed. SIGNIFICANCE STATEMENT: Quantitative proteomics analysis of drug-metabolizing enzymes and transporters in humans and preclinical species provides key physiological information that assists in the translation of in vitro and preclinical data to humans. This review provides the principles and applications of quantitative proteomics in characterizing in vitro, ex vivo, and preclinical models for translational research and interindividual variability prediction. Integration of these data into physiologically based pharmacokinetic modeling is proving to be critical for safe, effective, timely, and cost-effective drug development.

Clinical Implications of Altered Drug Transporter Abundance/Function and PBPK Modeling in Specific Populations: An ITC Perspective

The role of membrane transporters on pharmacokinetics (PKs), drug-drug interactions (DDIs), pharmacodynamics (PDs), and toxicity of drugs has been broadly recognized. However, our knowledge of modulation of transporter expression and/or function in the diseased patient population or specific populations, such as pediatrics or pregnancy, is still emerging. This white paper highlights recent advances in studying the changes in transporter expression and activity in various diseases (i.e., renal and hepatic impairment and cancer) and some specific populations (i.e., pediatrics and pregnancy) with the focus on clinical implications. Proposed alterations in transporter abundance and/or activity in diseased and specific populations are based on (i) quantitative transporter proteomic data and relative abundance in specific populations vs. healthy adults, (ii) clinical PKs, and emerging transporter biomarker and/or pharmacogenomic data, and (iii) physiologically-based pharmacokinetic modeling and simulation. The potential for altered PK, PD, and toxicity in these populations needs to be considered for drugs and their active metabolites in which transporter-mediated uptake/efflux is a major contributor to their absorption, distribution, and elimination pathways and/or associated DDI risk. In addition to best practices, this white paper discusses current challenges and knowledge gaps to study and quantitatively predict the effects of modulation in transporter activity in these populations, together with the perspectives from the International Transporter Consortium (ITC) on future directions.

The disposition of polychlorinated biphenyls (PCBs) differs between germ-free and conventional mice

The disposition of toxicants, such as polychlorinated biphenyls (PCBs), in germ-free (GF) vs. conventional (CV) mice has received little attention to date. Here, we investigate PCB levels in three-month-old female CV and GF mice exposed orally daily for 3 days to 0, 6, or 30 mg/kg body weight of the Fox River Mixture (FRM), an environmental PCB mixture. We euthanized animals 24 h after the final dose. PCB profiles in tissues differed from the FRM profile but were similar in tissues across all 4 PCB exposure groups. PCB levels in CV but not GF mice followed the difference in PCB dose. Importantly, PCB levels were higher in CV than GF mice exposed to the same dose. Hepatic cytochrome P450 enzyme or lipid levels did not explain these trends in PCB tissue levels. Thus, toxicity studies with CV and GF animals need to assess the toxicokinetics of the toxicant investigated. CAPSULE: PCB levels are typically higher in conventional than germ-free mice exposed to the same dose of PCBs.

Exploration of the Plausible Mechanism of Ethambutol Induced Ocular Toxicity by Using Proteomics Informed Physiologically Based Pharmacokinetic (PBPK) Modeling

PURPOSE

Ethambutol (EMB) is a first-line anti-tubercular drug that is known to cause optic neuropathy. The exact mechanism of its eye toxicity is unknown; however, proposition is metal chelating effect of both EMB and its metabolite 2,2'-(ethylenediamino)-dibutyric acid (EDBA). The latter is formed by sequential metabolism of EMB by alcohol dehydrogenases (ADHs) and aldehyde dehydrogenases (ALDHs). The purpose of this study was to predict the levels of drug and EDBA in the eye using physiologically based pharmacokinetic (PBPK) modeling.

METHODS

The PBPK model of EMB was developed using GastroPlus. The intrinsic hepatic clearance of ALDH, calculated by the model, was scaled down using proteomics data to estimate the rate of formation of EDBA in the eye. Additionally, the comparative permeability of EMB and EDBA was assessed by employing in silico and in vitro approaches. The rate of formation of EDBA in the eye and permeability data were then incorporated in a compartmental model to predict the ocular levels of EMB and EDBA.

RESULTS

The simulation results of compartmental model highlighted that there was an on-site formation of EDBA upon metabolism of EMB. Furthermore, in silico and in vitro studies revealed that EDBA possessed much lower permeability than EMB. These observations meant that once EDBA was formed in the eye, it was not permeated out and hence achieved higher ocular concentration.

CONCLUSION

The on-site formation of EDBA in the eye, its higher local concentration due to lower ocular clearance and its pre-known characteristic to chelate metal species better explains the ocular toxicity shown by EMB.

Interindividual Variability and Differential Tissue Abundance of Mitochondrial Amidoxime Reducing Component Enzymes in Humans

Mitochondrial amidoxime-reducing component (mARC) enzymes are molybdenum-containing proteins that metabolize a number of endobiotics and xenobiotics. The interindividual variability and differential tissue abundance of mARC1 and mARC2 were quantified using targeted proteomics in three types of tissue fractions: 1) pediatric liver tissue homogenates, 2) total membrane fraction of the paired liver and kidney samples from pediatric and adult donors, and 3) pooled S9 fractions of the liver, intestine, kidney, lung, and heart. The absolute levels of mARC1 and mARC2 in the pediatric liver homogenate were 40.08 ± 4.26 and 24.58 ± 4.02 pmol/mg homogenate protein, respectively, and were independent of age and sex. In the total membrane fraction of the paired liver and kidney samples, the abundance of hepatic mARC1 and mARC2 was comparable, whereas mARC2 abundance in the kidney was approximately 9-fold higher in comparison with mARC1. The analysis of the third set of samples (i.e., S9 fraction) revealed that mARC1 abundance in the kidney, intestine, and lung was 5- to 13-fold lower than the liver S9 abundance, whereas mARC2 abundance was approximately 3- and 16-fold lower in the intestine and lung than the liver S9, respectively. In contrast, the kidney mARC2 abundance in the S9 fraction was approximately 2.5-fold higher as compared with the hepatic mARC2 abundance. The abundance of mARC enzymes in the heart was below the limit of quantification (∼0.6 pmol/mg protein). The mARC enzyme abundance data presented here can be used to develop physiologically based pharmacokinetic models for the prediction of in vivo pharmacokinetics of mARC substrates. SIGNIFICANCE STATEMENT: A precise targeted quantitative proteomics method was developed and applied to quantify newly discovered drug-metabolizing enzymes, mARC1 and mARC2, in pediatric and adult tissue samples. The data suggest that mARC enzymes are ubiquitously expressed in an isoform-specific manner in the human liver, kidney, intestine, and lung, and the enzyme abundance is not associated with age and sex. These data are important for developing physiologically based pharmacokinetic models for the prediction of in vivo pharmacokinetics of mARC substrates.

Exploring the Interactome of Cytochrome P450 2E1 in Human Liver Microsomes with Chemical Crosslinking Mass Spectrometry

Aiming to elucidate the system-wide effects of the alcohol-induced increase in the content of cytochrome P450 2E1 (CYP2E1) on drug metabolism, we explored the array of its protein-protein interactions (interactome) in human liver microsomes (HLM) with chemical crosslinking mass spectrometry (CXMS). Our strategy employs membrane incorporation of purified CYP2E1 modified with photoreactive crosslinkers benzophenone-4-maleimide and 4-(N-succinimidylcarboxy)benzophenone. Exposure of bait-incorporated HLM samples to light was followed by isolating the His-tagged bait protein and its crosslinked aggregates on Ni-NTA agarose. Analyzing the individual bands of SDS-PAGE slabs of thereby isolated protein with the toolset of untargeted proteomics, we detected the crosslinked dimeric and trimeric complexes of CYP2E1 with other drug-metabolizing enzymes. Among the most extensively crosslinked partners of CYP2E1 are the cytochromes P450 2A6, 2C8, 3A4, 4A11, and 4F2, UDP-glucuronosyltransferases (UGTs) 1A and 2B, fatty aldehyde dehydrogenase (ALDH3A2), epoxide hydrolase 1 (EPHX1), disulfide oxidase 1α (ERO1L), and ribophorin II (RPN2). These results demonstrate the exploratory power of the proposed CXMS strategy and corroborate the concept of tight functional integration in the human drug-metabolizing ensemble through protein-protein interactions of the constituting enzymes.

Ontogeny of Scaling Factors for Pediatric Physiology-Based Pharmacokinetic Modeling and Simulation: Microsomal Protein Per Gram of Liver

Microsomal protein per gram of liver (MPPGL) is an important scaling factor for bottom-up physiology-based pharmacokinetic modeling and simulation, but data in pediatrics are limited. Therefore, MPPGL was determined in 160 liver samples from pediatric (n = 129) and adult (n = 31) donors obtained from four sources: the University of Maryland Brain and Tissue Bank (UMBTB), tissue retrieval services at the University of Minnesota and University of Pittsburgh, and Sekisui-Xenotech. Tissues were homogenized and subjected to differential centrifugation to prepare microsomes, and cytochrome c reductase activities in tissue homogenates and microsomes were used to estimate cytochrome P450 reductase (POR) activity as a marker of microsomal recovery; microsomal POR content was also assessed by quantitative proteomics. MPPGL values varied 5- to 10-fold within various age groups/developmental stages, and tissue source was identified as a contributing factor. Using a "trimmed" dataset comprised of samples ranging from 3 to 18 years of age common to the four sources, POR protein abundance and activity in microsomes and POR activity in homogenates was lower in UMBTB samples (autopsy) compared with other sources (perfused/flash-frozen). Regression analyses revealed that the UMBTB samples were driving an apparent age effect as no effect of age on log-transformed MPPGL values was observed when the UMBTB samples were excluded. We conclude that a mean±SD MPPGL value of 30.4±1.7 mg/g is representative between one month postnatal age and early adulthood. Potential source effects should be considered for studies involving tissue samples from multiple sources with different procurement and processing procedures. SIGNIFICANCE STATEMENT: Microsomal protein per gram of liver (MPPGL) is an important scaling factor for bottom up PBPK modeling and simulation, but data in pediatrics are limited. Although MPPGL varies 5- to 10-fold at a given developmental stage, a value of 30.4 ± 1.7 mg/g (mean ± SD) is representative between one month postnatal age and early adulthood. However, when tissue samples are obtained from multiple sources, different procurement and processing procedures may influence the results and should be taken into consideration.

Corrigendum: Ets-2 Propagates IL-6 Trans-Signaling Mediated Osteoclast-Like Changes in Human Rheumatoid Arthritis Synovial Fibroblast

[This corrects the article DOI: 10.3389/fimmu.2021.746503.].

Ets-2 Propagates IL-6 Trans-Signaling Mediated Osteoclast-Like Changes in Human Rheumatoid Arthritis Synovial Fibroblast

Rheumatoid arthritis synovial fibroblasts (RASFs) contribute to synovial inflammation and bone destruction by producing a pleiotropic cytokine interleukin-6 (IL-6). However, the molecular mechanisms through which IL-6 propels RASFs to contribute to bone loss are not fully understood. In the present study, we investigated the effect of IL-6 and IL-6 receptor (IL-6/IL-6R)-induced trans-signaling in human RASFs. IL-6 trans-signaling caused a significant increase in tartrate-resistant acid phosphatase (TRAP)-positive staining in RASFs and enhanced pit formation by ~3-fold in the osteogenic surface in vitro. IL-6/IL-6R caused dose-dependent increase in expression and nuclear translocation of transcription factor Ets2, which correlated with the expression of osteoclast-specific signature proteins RANKL, cathepsin B (CTSB), and cathepsin K (CTSK) in RASFs. Chromatin immunoprecipitation (ChIP) analysis of CTSB and CTSK promoters showed direct Ets2 binding and transcriptional activation upon IL-6/IL-6R stimulation. Knockdown of Ets2 significantly inhibited IL-6/IL-6R-induced RANKL, CTSB, and CTSK expression and TRAP staining in RASFs and suppressed markers of RASF invasive phenotype such as Thy1 and podoplanin (PDPN). Mass spectrometry analysis of the secretome identified 113 proteins produced by RASFs uniquely in response to IL-6/IL-6R that bioinformatically predicted its impact on metabolic reprogramming towards an osteoclast-like phenotype. These findings identified the role of Ets2 in IL-6 trans-signaling induced molecular reprogramming of RASFs to osteoclast-like cells and may contribute to RASF heterogeneity.

A Tribute to Professor Saranjit Singh - A Critical Thinker, Innovator, Mentor, and Educator

This commentary presents contributions and accomplishments of Professor Saranjit Singh, National Institute of Pharmaceutical Education and Research (NIPER), SAS Nagar, India, to pharmaceutical research and education. Prof. Singh completed his successful tenure in October 2021. Over his 40+ years of illustrious academic career, he trained 147 Masters and 15 PhD students in the fields of drug stability testing, degradation chemistry, impurity and metabolite characterization, and advanced analytical technologies. He has published ∼250 research articles, reviews, editorials, patent, book, and book chapters, and received numerous awards, including the Professor M.L. Khorana Memorial Lecture Award from the Indian Pharmaceutical Association (IPA) and the Outstanding Analyst and Eminent Analyst awards from the Indian Drug Manufacturers' Association (IDMA). This commentary highlights Prof. Singh's inspiring personal and renowned professional journey, including early life, education, career, accomplishments, as well as his services to academia, industry, and regulatory. By sharing the contributions and accomplishments of Prof. Singh, we strongly believe that his story will inspire the next generation of scientists to continue his legacy to advance the field.

The Effect of Daily Methylsulfonylmethane (MSM) Consumption on High-Density Lipoprotein Cholesterol in Healthy Overweight and Obese Adults: A Randomized Controlled Trial

Interventions to decrease inflammation and improve metabolic function hold promise for the prevention of obesity-related diseases. Methylsulfonylmethane (MSM) is a naturally occurring compound that demonstrates antioxidant and anti-inflammatory effects. Improvements in measures of metabolic health have been observed in mouse models of obesity and diabetes following MSM treatment. However, the effects of MSM on obesity-related diseases in humans have not been investigated. Therefore, the purpose of this investigation was to determine whether MSM supplementation improves cardiometabolic health, and markers of inflammation and oxidative status. A randomized, double-blind, placebo-controlled design was utilized with a total of 22 overweight or obese adults completing the study. Participants received either a placebo (white rice flour) or 3 g MSM daily for 16 weeks. Measurements occurred at baseline and after 4, 8, and 16 weeks. Outcome measures included fasting glucose, insulin, blood lipids, blood pressure, body composition, metabolic rate, and markers of inflammation and oxidative status. The primary finding of this work shows that high-density lipoprotein cholesterol was elevated at 8 and 16 weeks of daily MSM consumption compared to baseline, (p = 0.008, p = 0.013). Our findings indicate that MSM supplementation may improve the cholesterol profile by resulting in higher levels of high-density lipoprotein cholesterol.

Assembling the P450 puzzle: on the sources of nonadditivity in drug metabolism

There is an increasing number of indications of an oversimplification in the premise that the cumulative properties of the human drug-metabolizing ensemble represent a simple aggregate of the properties of the constituting enzymes. Recent studies of the functional effects of hetero-association of multiple cytochrome P450 species and their interactions with metabolically related enzymes revealed a tight integration in the drug-metabolizing ensemble. In our opinion, the sources of interindividual variability in drug metabolism can be elucidated only when considering this ensemble as a multienzyme system, the functional parameters of which are determined by interactions between its constituents. In this article, we present a conceptual model providing a mechanistic explanation for the functional effects of the interactions between multiple P450 species and propose a clue to understanding the nonadditive behavior of the drug-metabolizing ensemble.

Probing functional interactions between cytochromes P450 with principal component analysis of substrate saturation profiles and targeted proteomics

We investigated the correspondence between drug metabolism routes and the composition of the P450 ensemble in human liver microsomes (HLM). As a probe, we used Coumarin 152 (C152), a fluorogenic substrate metabolized by multiple P450 species. Studying the substrate-saturation profiles (SSP) in seven pooled HLM preparations, we sought to correlate them with the P450 pool's composition characterized by targeted proteomics. This analysis, complemented with the assays with specific inhibitors of CYP3A4 and CYP2C19, the primary C152 metabolizers, demonstrated a significant contrast between different HLM samples. To unveil the source of these differences, we implemented Principal Component Analysis (PCA) of the SSP series obtained with HLM samples with a known composition of the P450 pool. Our analysis revealed that the parameters of C152 metabolism are primarily determined by the content of CYP2A6, CYP2B6, CYP2C8, CYP2E1, and CYP3A5 of those only CYP2B6 and CYP3A5 can metabolize C152. To validate this finding, we studied the effect of enriching HLM with CYP2A6, CYP2E1, and CYP3A5. The incorporation of CYP3A5 into HLM decreases the rate of C152 metabolism while increasing the role of CYP2B6 in its turnover. In contrast, incorporation of CYP2A6 and CYP2E1 reroutes the C152 demethylation towards some P450 enzyme with a moderate affinity to the substrate, most likely CYP3A4. Our results reveal a sharp non-additivity of the individual P450 properties and suggest a pivotal role of P450-P450 interactions in determining drug metabolism routes. This study demonstrates the high potential of our new PCA-based approach in unveiling functional interrelationships between different P450 species.

Functional characterization of novel rare CYP2A6 variants and potential implications for clinical outcomes

CYP2A6 activity, phenotyped by the nicotine metabolite ratio (NMR), is a predictor of several smoking behaviors, including cessation and smoking‐related disease risk. The heritability of the NMR is 60–80%, yet weighted genetic risk scores (wGRSs) based on common variants explain only 30–35%. Rare variants (minor allele frequency <1%) are hypothesized to explain some of this missing heritability. We present two targeted sequencing studies where rare protein‐coding variants are functionally characterized in vivo, in silico, and in vitro to examine this hypothesis. In a smoking cessation trial, 1687 individuals were sequenced; characterization measures included the in vivo NMR, in vitro protein expression, and metabolic activity measured from recombinant proteins. In a human liver bank, 312 human liver samples were sequenced; measures included RNA expression, protein expression, and metabolic activity from extracted liver tissue. In total, 38 of 47 rare coding variants identified were novel; characterizations ranged from gain‐of‐function to loss‐of‐function. On a population level, the portion of NMR variation explained by the rare coding variants was small (~1%). However, upon incorporation, the accuracy of the wGRS was improved for individuals with rare protein‐coding variants (i.e., the residuals were reduced), and approximately one‐third of these individuals (12/39) were re‐assigned from normal to slow metabolizer status. Rare coding variants can alter an individual’s CYP2A6 activity; their integration into wGRSs through precise functional characterization is necessary to accurately assess clinical outcomes and achieve precision medicine for all. Investigation into noncoding variants is warranted to further explain the missing heritability in the NMR.

A Comprehensive Assessment of Associations between Prenatal Phthalate Exposure and the Placental Transcriptomic Landscape

BACKGROUND

Phthalates are commonly used endocrine-disrupting chemicals that are ubiquitous in the general population. Prenatal phthalate exposure may alter placental physiology and fetal development, leading to adverse perinatal and childhood health outcomes.

OBJECTIVE

We examined associations between prenatal phthalate exposure in the second and third trimesters and the placental transcriptome at birth, including genes and long noncoding RNAs (lncRNAs), to gain insight into potential mechanisms of action during fetal development.

METHODS

The ECHO PATHWAYs consortium quantified 21 urinary phthalate metabolites from 760 women enrolled in the CANDLE study (Shelby County, TN) using high-performance liquid chromatography-tandem mass spectrometry. Placental transcriptomic data were obtained using paired-end RNA sequencing. Linear models were fitted to estimate separate associations between maternal urinary phthalate metabolite concentration during the second and third trimester and placental gene expression at birth, adjusted for confounding variables. Genes were considered differentially expressed at a Benjamini-Hochberg false discovery rate (FDR) p<0.05. Associations between phthalate metabolites and biological pathways were identified using self-contained gene set testing and considered significantly altered with an FDR-adjusted p<0.2.

RESULTS

We observed significant associations between second-trimester phthalate metabolites mono (carboxyisooctyl) phthalate (MCIOP), mono-2-ethyl-5-carboxypentyl phthalate, and mono-2-ethyl-5-oxohexyl phthalate and 18 genes in total, including four lncRNAs. Specifically, placental expression of NEAT1 was associated with multiple phthalate metabolites. Third-trimester MCIOP and mono-isobutyl phthalate concentrations were significantly associated with placental expression of 18 genes and two genes, respectively. Expression of genes within 27 biological pathways was associated with mono-methyl phthalate, MCIOP, and monoethyl phthalate concentrations.

DISCUSSION

To our knowledge, this is the first genome-wide assessment of the relationship between the placental transcriptome at birth and prenatal phthalate exposure in a large and diverse birth cohort. We identified numerous genes and lncRNAs associated with prenatal phthalate exposure. These associations mirror findings from other epidemiological and in vitro analyses and may provide insight into biological pathways affected in utero by phthalate exposure. https://doi.org/10.1289/EHP8973.