Synthesis and Reactivity of Masked Organic Sulfates

Nature offers a variety of structurally unique, sulfated endobiotics including sulfated glycosaminoglycans, sulfated tyrosine peptides, sulfated steroids/bile acids/catecholamines. Sulfated molecules display a large number of biological activities including antithrombotic, antimicrobial, anticancer, anti-inflammatory, and others, which arise from modulation of intracellular signaling and enhanced in vivo retention of certain hormones. These characteristics position sulfated molecules very favorably as drug-like agents. However, few have reached the clinic. Major hurdles exist in realizing sulfated molecules as drugs. This state-of-the-art has been transformed through recent works on the development of sulfate masking technologies for both alkyl (sulfated carbohydrates, sulfated steroids) and aryl (sTyr-bearing peptides/proteins, sulfated flavonoids) sulfates. This review compiles the literature on different strategies implemented for different types of sulfate groups. Starting from early efforts in protection of sulfate groups to the design of newer SuFEx, trichloroethyl, and gem-dimethyl-based protection technologies, this review presents the evolution and application of concepts in realizing highly diverse, sulfated molecules as candidate drugs and/or prodrugs. Overall, the newer strategies for sulfate masking and demasking are likely to greatly enhance the design and development of sulfated molecules as non-toxic drugs of the future.

Understanding the Impacts of Presystemic Metabolism on the Human Oral Bioavailability of Chemicals

Animal-free new approach methods promote chemical assessments based on the comparison between in vitro bioactivity and human internal concentrations, which necessitates a dependable knowledge of human oral bioavailability, i.e., the fraction of an orally ingested chemical that escapes from presystemic ("first-pass") metabolic processes and eventually enters systemic circulation. Using a physiologically based toxicokinetic model, we show how human oral bioavailability is impacted by presystemic metabolism within the gut lumen, gut wall, and liver and how this impact differs among chemicals with various permeability and stability properties. Our results highlight the gut lumen as a primary site of presystemic metabolism of certain chemicals, such as di-2-ethylhexyl phthalate (DEHP), for which the gut lumen may even exceed the liver in importance of presystemic metabolism due to these metabolic processes occurring in sequence. For chemicals with low transmembrane permeability and low stability, metabolism within the gut lumen is the most remarkable of the three presystemic metabolic processes. Notably, for chemicals that undergo substantial metabolism within the gut lumen, where the metabolites have high permeability, there is a notable discrepancy between the "theoretical bioavailability" (bioavailability of the unchanged parent compound) and the "apparent bioavailability" in measurement practices (bioavailability inferred from measured metabolites). Our work highlights the importance of considering presystemic metabolism, notably within the gut lumen, in human exposure and toxicokinetic modeling.

Physiologically based Kinetic Modeling-Facilitated Quantitative In Vitro to In Vivo Extrapolation to Predict the Effects of Aloe-Emodin in Rats and Humans

Aloe-emodin, a natural hydroxyanthraquinone, exerts both adverse and protective effects. This study aimed at investigating these potential effects of aloe-emodin in humans upon the use of food supplements and herbal medicines using a physiologically based kinetic (PBK) modeling-facilitated quantitative in vitro to in vivo extrapolation (QIVIVE) approach. For this, PBK models in rats and humans were established for aloe-emodin including its active metabolite rhein and used to convert in vitro data on hepatotoxicity, nephrotoxicity, reactive oxidative species (ROS) generation, and Nrf2 induction to corresponding in vivo dose-response curves, from which points of departure (PODs) were derived by BMD analysis. The derived PODs were subsequently compared to the estimated daily intakes (EDIs) resulting from the use of food supplements or herbal medicines. It is concluded that the dose levels of aloe-emodin from food supplements or herbal medicines are unlikely to induce toxicity, ROS generation, or Nrf2 activation in liver and kidney.

CYP2J2-mediated metabolism of arachidonic acid in heart: A review of its kinetics, inhibition and role in heart rhythm control

Cytochrome P450 2 J2 (CYP2J2) is primarily expressed extrahepatically and is the predominant epoxygenase in human cardiac tissues. This highlights its key role in the metabolism of endogenous substrates. Significant scientific interest lies in cardiac CYP2J2 metabolism of arachidonic acid (AA), an omega-6 polyunsaturated fatty acid, to regioisomeric bioactive epoxyeicosatrienoic acid (EET) metabolites that show cardioprotective effects including regulation of cardiac electrophysiology. From an in vitro perspective, the accurate characterization of the kinetics of CYP2J2 metabolism of AA including its inhibition and inactivation by drugs could be useful in facilitating in vitro-in vivo extrapolations to predict drug-AA interactions in drug discovery and development. In this review, background information on the structure, regulation and expression of CYP2J2 in human heart is presented alongside AA and EETs as its endogenous substrate and metabolites. The in vitro and in vivo implications of the kinetics of this endogenous metabolic pathway as well as its perturbation via inhibition and inactivation by drugs are elaborated. Additionally, the role of CYP2J2-mediated metabolism of AA to EETs in cardiac electrophysiology will be expounded.

Ontogeny of Scaling Factors for Pediatric Physiologically Based Pharmacokinetic Modeling and Simulation: Cytosolic Protein Per Gram of Liver

Scaling factors are necessary for translating in vitro drug biotransformation data to in vivo clearance values in physiologically-based pharmacokinetic modeling and simulation. Values for microsomal protein per gram of liver are available from several sources for use as a scaling factor to estimate hepatic clearance from microsomal drug biotransformation data. However, data regarding the distribution of cytosolic protein per gram of liver (CPPGL) values across the lifespan are limited, and sparse pediatric data have been published to date. Thus, CPPGL was determined in 160 liver samples from pediatric (n = 129) and adult (n = 31) donors obtained from multiple sources: the University of Maryland Brain and Tissue Bank, tissue retrieval services at the University of Minnesota and University of Pittsburgh, and Sekisui-XenoTech. Tissues were homogenized and subjected to differential centrifugation to isolate cytosolic fractions. Cytosolic protein content was determined by BCA assay. CPPGL varied from two- to sixfold within each age group/developmental stage. Tissue source and sex did not contribute substantially to variability in protein content. Regression analyses revealed minimal change in CPPGL over the first two decades of life (logCPPGL increases 0.1 mg/g per decade). A mean ± S.D. CPPGL value of 44.4 ± 17.4 mg/g or median 41.0 mg/g is representative of values observed between birth and early adulthood (0-18 years, n = 129). SIGNIFICANCE STATEMENT: Cytosolic protein per gram of liver (CPPGL) is a scaling factor required for physiologically based pharmacokinetic modeling and simulation of drug biotransformation by cytosolic enzymes, but pediatric data are limited. Although CPPGL varies from two- to sixfold within developmental stages, a value of 44.4 ± 17.4 mg/g (mean ± S.D.) is representative of the pediatric period (0-18 years, n = 129).

Challenges and Opportunities for In Vitro-In Vivo Extrapolation of Aldehyde Oxidase-Mediated Clearance: Toward a Roadmap for Quantitative Translation

Underestimation of aldehyde oxidase (AO)-mediated clearance by current in vitro assays leads to uncertainty in human dose projections, thereby reducing the likelihood of success in drug development. In the present study we first evaluated the current drug development practices for AO substrates. Next, the overall predictive performance of in vitro-in vivo extrapolation of unbound hepatic intrinsic clearance (CLint,u) and unbound hepatic intrinsic clearance by AO (CLint,u,AO) was assessed using a comprehensive literature database of in vitro (human cytosol/S9/hepatocytes) and in vivo (intravenous/oral) data collated for 22 AO substrates (total of 100 datapoints from multiple studies). Correction for unbound fraction in the incubation was done by experimental data or in silico predictions. The fraction metabolized by AO (fmAO) determined via in vitro/in vivo approaches was found to be highly variable. The geometric mean fold errors (gmfe) for scaled CLint,u (mL/min/kg) were 10.4 for human hepatocytes, 5.6 for human liver cytosols, and 5.0 for human liver S9, respectively. Application of these gmfe's as empirical scaling factors improved predictions (45%-57% within twofold of observed) compared with no correction (11%-27% within twofold), with the scaling factors qualified by leave-one-out cross-validation. A road map for quantitative translation was then proposed following a critical evaluation on the in vitro and clinical methodology to estimate in vivo fmAO In conclusion, the study provides the most robust system-specific empirical scaling factors to date as a pragmatic approach for the prediction of in vivo CLint,u,AO in the early stages of drug development. SIGNIFICANCE STATEMENT: Confidence remains low when predicting in vivo clearance of AO substrates using in vitro systems, leading to de-prioritization of AO substrates from the drug development pipeline to mitigate risk of unexpected and costly in vivo impact. The current study establishes a set of empirical scaling factors as a pragmatic tool to improve predictability of in vivo AO clearance. Developing clinical pharmacology strategies for AO substrates by utilizing mass balance/clinical drug-drug interaction data will help build confidence in fmAO.

Intestinal organoids as an in vitro platform to characterize disposition, metabolism, and safety profile of small molecules

Intestinal organoids derived from LGR5+ adult stem cells allow for long-term culturing, more closely resemble human physiology than traditional intestinal models, like Caco-2, and have been established for several species. Here we evaluated intestinal organoids for drug disposition, metabolism, and safety applications. Enterocyte-enriched human duodenal organoids were cultured as monolayers to enable bidirectional transport studies. 3D enterocyte-enriched human duodenal and colonic organoids were incubated with probe substrates of major intestinal drug metabolizing enzymes (DMEs). To distinguish human intestinal toxic (high incidence of diarrhea in clinical trials and/or black box warning related to intestinal side effects) from non-intestinal toxic compounds, ATP-based cell viability was used as a readout, and compounds were ranked based on their IC50 values in relation to their 30-times maximal total plasma concentration (Cmax). To assess if rat and dog organoids reproduced the respective in vivo intestinal safety profiles, ATP-based viability was assessed in rat and dog organoids and compared to in vivo intestinal findings when available. Human duodenal monolayers discriminated high and low permeable compounds and demonstrated functional activity for the main efflux transporters Multi drug resistant protein 1 (MDR1, P-glycoprotein P-gp) and Breast cancer resistant protein (BCRP). Human 3D duodenal and colonic organoids also showed metabolic activity for the main intestinal phase I and II DMEs. Organoids derived from specific intestinal segments showed activity differences in line with reported DMEs expression. Undifferentiated human organoids accurately distinguished all but one compound from the test set of non-toxic and toxic drugs. Cytotoxicity in rat and dog organoids correlated with preclinical toxicity findings and observed species sensitivity differences between human, rat, and dog organoids. In conclusion, the data suggest intestinal organoids are suitable in vitro tools for drug disposition, metabolism, and intestinal toxicity endpoints. The possibility to use organoids from different species, and intestinal segment holds great potential for cross-species and regional comparisons.

Physiologically Based Biopharmaceutics Model for Selumetinib Food Effect Investigation and Capsule Dissolution Safe Space - Part I: Adults

OBJECTIVE

A physiologically based biopharmaceutics model (PBBM) was developed to mechanistically investigate the effect of formulation and food on selumetinib pharmacokinetics.

METHODS

Selumetinib is presented as a hydrogen sulfate salt, and in vitro and in vivo data were used to verify the precipitation rate to apply to simulations. Dissolution profiles observed for capsules and granules were used to derive product-particle size distributions for model input. The PBBM incorporated gut efflux and first-pass gut metabolism, based on intravenous and oral pharmacokinetic data, alongside in vitro data for the main enzyme isoform and P-glycoprotein efflux. The PBBM was validated across eight clinical scenarios.

RESULTS

The quality-control dissolution method for selumetinib capsules was found to be clinically relevant through PBBM validation. A safe space for capsule dissolution was established using a virtual batch. The effect of food (low fat vs high fat) on capsules and granules was elucidated by the PBBM. For capsules, a lower amount was dissolved in the fed state due to a pH increase in the stomach followed by higher precipitation in the small intestine. First-pass gut extraction is higher for capsules in the fed state due to drug dilution in the stomach chyme and reduced concentration in the lumen. The enteric-coated granules dissolve more slowly than capsules after stomach emptying, attenuating the difference in first-pass gut extraction between prandial states.

CONCLUSIONS

The PBBM was instrumental in understanding and explaining the different behaviors of the selumetinib formulations. The model can be used to predict the impact of food in humans.

Physiologically-Based Kinetic Modeling Predicts Similar In Vivo Relative Potency of Senecionine N-Oxide for Rat and Human at Realistic Low Exposure Levels

SCOPE

This study aims to determine if previously developed physiologically-based kinetic (PBK) model in rat can be modified for senecionine (SEN) and its N-oxide (SENO), and be used to investigate potential species differences between rat and human in relative potency (REP) of the N-oxide relative to the parent pyrrolizidine alkaloid (PA).

METHODS AND RESULTS

In vitro derived kinetic parameters including the apparent maximum velocities (V_max ) and Michaelis-Menten constants (K_m ) for SENO reduction and SEN clearance are used to define the PBK models. The rat model is validated with published animal data, and the toxicokinetic profiles of SEN from either orally-administered SENO or SEN are simulated. REP values of SENO relative to SEN amount to 0.84 and 0.89 in rat and human, respectively.

CONCLUSION

The REP value can be dose- and species-dependent, with the values for rat and human being comparable at low realistic exposure scenarios. In summary, PBK modeling serves as a valuable New Approach Methodology (NAM) tool for predicting REP values of PA-N-oxides and may actually result in more accurate REP values for human risk assessment than what would be defined using in vivo animal experiments.

Stereoselective Metabolism of Bupropion to Active Metabolites in Cellular Fractions of Human Liver and Intestine

Striking stereoselective disposition of the antidepressant and smoking cessation aid bupropion (BUP) and its active metabolites observed clinically influence patients' response to BUP therapy and its clinically important drug-drug interactions (DDI) with CYP2D6 substrates. However, understanding of the biochemical mechanisms responsible is incomplete. This study comprehensively examined hepatic and extrahepatic stereoselective metabolism of BUP in vitro Racemic-, R-, and S-BUP were incubated separately with pooled cellular fractions of human liver [microsomes (HLMs), S9 fractions (HLS9s), and cytosols (HLCs)] and intestinal [microsomes (HIMs), S9 fractions (HIS9s), and cytosols (HICs)] and cofactors. Formations of diastereomers of 4-hydroxyBUP (OHBUP), threohydroBUP (THBUP), and erythrohydroBUP (EHBUP) were quantified using a novel chiral ultra-high performance liquid chromatography/tandem mass spectrometry method. Racemic BUP (but not R- or S-BUP) was found suitable to determine stereoselective metabolism of BUP; both enantiomers showed complete racemization. Compared with that of RR-THBUP, the in vitro intrinsic clearance (Clint) for the formation of SS-THBUP was 42-, 19-, and 8.3-fold higher in HLMs, HLS9 fractions, and HLCs, respectively; Clint for the formation of SS-OHBUP and RS-EHBUP was also higher (2.7- to 3.9-fold) than their R-derived counterparts. In cellular fractions of human intestine, ≥ 95% of total reduction was accounted by the formation of RR-THBUP. Ours is the first to demonstrate marked stereoselective reduction of BUP in HLCs, HIMs, HIS9 fractions, and HICs, providing the first evidence for tissue- and cellular fraction-dependent stereoselective metabolism of BUP. These data may serve as the first critical step toward understanding factors dictating BUP's stereoselective disposition, effects, and DDI risks. SIGNIFICANCE STATEMENT: This work provides a deeper insight into bupropion (BUP) stereoselective oxidation and reduction to active metabolites in cellular fractions of human liver and intestine tissues. The results demonstrate tissue- and cellular fraction-dependent stereospecific metabolism of BUP. These data may improve prediction of BUP stereoselective disposition and understanding of BUP's effects and CYP2D6-dependent drug-drug interaction in vivo.

Investigating the uncertainty of prediction accuracy for the application of physiologically based pharmacokinetic models to animal-free risk assessment of cosmetic ingredients

Physiologically based pharmacokinetic (PBPK) models are considered useful tools in animal-free risk assessment. To utilize PBPK models for risk assessment, it is necessary to compare their reliability with in vivo data. However, obtaining in vivo pharmacokinetics data for cosmetic ingredients is difficult, complicating the utilization of PBPK models for risk assessment. In this study, to utilize PBPK models for risk assessment without accuracy evaluation, we proposed a novel concept-the modeling uncertainty factor (MUF). By calculating the prediction accuracy for 150 compounds, we established that using in vitro data for metabolism-related parameters and limiting the applicability domain increase the prediction accuracy of a PBPK model. Based on the 97.5th percentile of prediction accuracy, MUF was defined at 10 for the area under the plasma concentration curve and 6 for C_max. A case study on animal-free risk assessment was conducted for bisphenol A using these MUFs. As this study was conducted mainly on pharmaceuticals, further investigation using cosmetic ingredients is pivotal. However, since internal exposure is essential in realizing animal-free risk assessment, our concept will serve as a useful tool to predict plasma concentrations without using in vivo data.

Towards best use and regulatory acceptance of generic physiologically based kinetic (PBK) models for in vitro-to-in vivo extrapolation (IVIVE) in chemical risk assessment

With an increasing need to incorporate new approach methodologies (NAMs) in chemical risk assessment and the concomitant need to phase out animal testing, the interpretation of in vitro assay readouts for quantitative hazard characterisation becomes more important. Physiologically based kinetic (PBK) models, which simulate the fate of chemicals in tissues of the body, play an essential role in extrapolating in vitro effect concentrations to in vivo bioequivalent exposures. As PBK-based testing approaches evolve, it will become essential to standardise PBK modelling approaches towards a consensus approach that can be used in quantitative in vitro-to-in vivo extrapolation (QIVIVE) studies for regulatory chemical risk assessment based on in vitro assays. Based on results of an ECETOC expert workshop, steps are recommended that can improve regulatory adoption: (1) define context and implementation, taking into consideration model complexity for building fit-for-purpose PBK models, (2) harmonise physiological input parameters and their distribution and define criteria for quality chemical-specific parameters, especially in the absence of in vivo data, (3) apply Good Modelling Practices (GMP) to achieve transparency and design a stepwise approach for PBK model development for risk assessors, (4) evaluate model predictions using alternatives to in vivo PK data including read-across approaches, (5) use case studies to facilitate discussions between modellers and regulators of chemical risk assessment. Proof-of-concepts of generic PBK modelling approaches are published in the scientific literature at an increasing rate. Working on the previously proposed steps is, therefore, needed to gain confidence in PBK modelling approaches for regulatory use.

The Role of Uptake and Efflux Transporters in the Disposition of Glucuronide and Sulfate Conjugates

Glucuronidation and sulfation are the most typical phase II metabolic reactions of drugs. The resulting glucuronide and sulfate conjugates are generally considered inactive and safe. They may, however, be the most prominent drug-related material in the circulation and excreta of humans. The glucuronide and sulfate metabolites of drugs typically have limited cell membrane permeability and subsequently, their distribution and excretion from the human body requires transport proteins. Uptake transporters, such as organic anion transporters (OATs and OATPs), mediate the uptake of conjugates into the liver and kidney, while efflux transporters, such as multidrug resistance proteins (MRPs) and breast cancer resistance protein (BCRP), mediate expulsion of conjugates into bile, urine and the intestinal lumen. Understanding the active transport of conjugated drug metabolites is important for predicting the fate of a drug in the body and its safety and efficacy. The aim of this review is to compile the understanding of transporter-mediated disposition of phase II conjugates. We review the literature on hepatic, intestinal and renal uptake transporters participating in the transport of glucuronide and sulfate metabolites of drugs, other xenobiotics and endobiotics. In addition, we provide an update on the involvement of efflux transporters in the disposition of glucuronide and sulfate metabolites. Finally, we discuss the interplay between uptake and efflux transport in the intestine, liver and kidneys as well as the role of transporters in glucuronide and sulfate conjugate toxicity, drug interactions, pharmacogenetics and species differences.

Use of Physiologically Based Pharmacokinetic Modeling to Predict Human Gut Microbial Conversion of Daidzein to S-Equol

A physiologically based pharmacokinetic (PBPK) model was developed for daidzein and its metabolite S-equol. Anaerobic in vitro incubations of pooled fecal samples from S-equol producers and nonproducers allowed definition of the kinetic constants. PBPK model-based predictions for the maximum daidzein plasma concentration (C_max) were comparable to literature data. The predictions also revealed that the C_max of S-equol in producers was only up to 0.22% that of daidzein, indicating that despite its higher estrogenicity, S-equol is likely to contribute to the overall estrogenicity upon human daidzein exposure to a only limited extent. An interspecies comparison between humans and rats revealed that the catalytic efficiency for S-equol formation in rats was 210-fold higher than that of human S-equol producers. The described in vitro-in silico strategy provides a proof-of-principle on how to include microbial metabolism in humans in PBPK modeling as part of the development of new approach methodologies (NAMs).

PBK Model-Based Prediction of Intestinal Microbial and Host Metabolism of Zearalenone and Consequences for its Estrogenicity

SCOPE

The aim of the present study is to develop physiologically-based kinetic (PBK) models for rat and human that include intestinal microbial and hepatic metabolism of zearalenone (ZEN) in order to predict systemic concentrations of ZEN and to obtain insight in the contribution of metabolism by the intestinal microbiota to the overall metabolism of ZEN.

METHODS AND RESULTS

In vitro derived kinetic parameters, apparent maximum velocities (Vmax ) and Michaelis-Menten constants (Km ) for liver and intestinal microbial metabolism of ZEN are included in the PBK models. The models include a sub-model for the metabolite, α-zearalenol (α-ZEL), a metabolite known to be 60-times more potent as an estrogen than ZEN. Integrating intestinal microbial ZEN metabolism into the PBK models revealed that hepatic metabolism drives the formation of α-ZEL. Furthermore, the models predicted that at the tolerable daily intake (TDI) of 0.25 µg kg-1 bw the internal concentration of ZEN and α-ZEL are three-orders of magnitude below concentrations reported to induce estrogenicity in vitro.

CONCLUSION

It is concluded that combining kinetic data on liver and intestinal microbial metabolism in a PBK model facilitates a holistic view on the role of the intestinal microbiota in the overall metabolism of the foodborne xenobiotic ZEN and its bioactivation to α-ZEL.

Physiologically based kinetic modelling predicts the in vivo relative potency of riddelliine N-oxide compared to riddelliine in rat to be dose dependent

Pyrrolizidine alkaloids (PAs) are toxic plant constituents occurring often in their N-oxide form. This raises the question on the relative potency (REP) values of PA-N-oxides compared to the corresponding parent PAs. The present study aims to quantify the in vivo REP value of riddelliine N-oxide compared to riddelliine using physiologically based kinetic (PBK) modelling, taking into account that the toxicity of riddelliine N-oxide depends on its conversion to riddelliine by intestinal microbiota and in the liver. The models predicted a lower C_max and higher T_max for the blood concentration of riddelliine upon oral administration of riddelliine N-oxide compared to the C_max and T_max predicted for an equimolar oral dose of riddelliine. Comparison of the area under the riddelliine concentration–time curve (AUC_RID) obtained upon dosing either the N-oxide or riddelliine itself revealed a ratio of 0.67, which reflects the in vivo REP for riddelliine N-oxide compared to riddelliine, and appeared to closely match the REP value derived from available in vivo data. The models also predicted that the REP value will decrease with increasing dose level, because of saturation of riddelliine N-oxide reduction by the intestinal microbiota and of riddelliine clearance by the liver. It is concluded that PBK modeling provides a way to define in vivo REP values of PA-N-oxides as compared to their parent PAs, without a need for animal experiments.

PBPK Modeling as a Tool for Predicting and Understanding Intestinal Metabolism of Uridine 5'-Diphospho-glucuronosyltransferase Substrates

Uridine 5′-diphospho-glucuronosyltransferases (UGTs) are expressed in the small intestines, but prediction of first-pass extraction from the related metabolism is not well studied. This work assesses physiologically based pharmacokinetic (PBPK) modeling as a tool for predicting intestinal metabolism due to UGTs in the human gastrointestinal tract. Available data for intestinal UGT expression levels and in vitro approaches that can be used to predict intestinal metabolism of UGT substrates are reviewed. Human PBPK models for UGT substrates with varying extents of UGT-mediated intestinal metabolism (lorazepam, oxazepam, naloxone, zidovudine, cabotegravir, raltegravir, and dolutegravir) have demonstrated utility for predicting the extent of intestinal metabolism. Drug–drug interactions (DDIs) of UGT1A1 substrates dolutegravir and raltegravir with UGT1A1 inhibitor atazanavir have been simulated, and the role of intestinal metabolism in these clinical DDIs examined. Utility of an in silico tool for predicting substrate specificity for UGTs is discussed. Improved in vitro tools to study metabolism for UGT compounds, such as coculture models for low clearance compounds and better understanding of optimal conditions for in vitro studies, may provide an opportunity for improved in vitro–in vivo extrapolation (IVIVE) and prospective predictions. PBPK modeling shows promise as a useful tool for predicting intestinal metabolism for UGT substrates.

Identification of human UDP-glucuronosyltransferase and sulfotransferase as responsible for the metabolism of dotinurad, a novel selective urate reabsorption inhibitor

Dotinurad, a novel selective urate reabsorption inhibitor, is used to treat hyperuricemia. In humans, orally administered dotinurad is excreted mainly as glucuronide and sulfate conjugates in urine. To identify the isoforms of UDP-glucuronosyltransferase (UGT) and sulfotransferase (SULT) involved in dotinurad glucuronidation and sulfation, microsome and cytosol fractions of liver, intestine, kidney, and lung tissues (cytosol only) were analyzed along with recombinant human UGT and SULT isoforms. Dotinurad was mainly metabolized to its glucuronide conjugate by human liver microsomes (HLMs), and the glucuronidation followed the two-enzyme Michaelis-Menten equation. Among the recombinant human UGT isoforms expressed in the liver, UGT1A1, UGT1A3, UGT1A9, and UGT2B7 catalyzed dotinurad glucuronidation. Based on inhibition analysis using HLMs, bilirubin, imipramine, and diflunisal decreased glucuronosyltransferase activities by 45.5, 22.3, and 22.2%, respectively. Diflunisal and 3'-azido-3'-deoxythymidine, in the presence of 1% BSA, decreased glucuronosyltransferase activities by 21.1 and 13.4%, respectively. Dotinurad was metabolized to its sulfate conjugate by human liver cytosol (HLC) and human intestinal cytosol (HIC) samples, with the sulfation reaction in HLC samples following the two-enzyme Michaelis-Menten equation and that in HIC samples following the Michaelis-Menten equation. All eight recombinant human SULT isoforms used herein catalyzed dotinurad sulfation. Gavestinel decreased sulfotransferase activity by 15.3% in HLC samples, and salbutamol decreased sulfotransferase activity by 68.4% in HIC samples. These results suggest that dotinurad glucuronidation is catalyzed mainly by UGT1A1, UGT1A3, UGT1A9, and UGT2B7, whereas its sulfation is catalyzed by many SULT isoforms, including SULT1B1 and SULT1A3. Significance Statement The identification of enzymes involved in drug metabolism is important to predicting drug-drug interactions (DDIs) and interindividual variability for safe drug use. The present study revealed that dotinurad glucuronidation is catalyzed mainly by UGT1A1, UGT1A3, UGT1A9, and UGT2B7 and that its sulfation is catalyzed by many SULT isoforms, including SULT1B1 and SULT1A3. Therefore, dotinurad, a selective urate reabsorption inhibitor, is considered safe for use with a small risk of DDIs and low interindividual variability.

Impact of simulated lung fluid components on the solubility of inhaled drugs and predicted in vivo performance

Orally inhaled products (OIPs) are gaining increased attention, as pulmonary delivery is a preferred route for the treatment of various diseases. Yet, the field of inhalation biopharmaceutics is still in development phase. For a successful correlation between various in vitro data obtained during formulation characterization and in vivo performance, it is necessary to understand the impact of parameters such as solubility and dissolution of drugs. In this work, we used in vitro-in silico feedback-feedforward approach to gain a better insight into the biopharmaceutics behavior of inhaled Salbutamol Sulphate (SS) and Budesonide (BUD). The thorough characterization of the in vitro test media and the impact of different in vitro fluid components such as lipids and protein on the solubility of aforementioned drugs was studied. These results were subsequently used as an input into the developed in silico models to investigate potential PK parameter changes in vivo. Results revealed that media comprising lipids and albumin were the most biorelevant and impacted the solubility of BUD the most. On the contrary, no notable impact was seen in case of SS. The use of simple media such as phosphate buffer saline (PBS) might be sufficient to use in solubility studies of the highly soluble and permeable drugs. However, its use for the poorly soluble drugs is limited due to the greater potential for interactions within in vivo environment. The use of in silico tools showed that the model response varies, depending on the used media. Therefore, this work highlights the relevance of carefully selecting the media composition when investigating solubility and dissolution behavior, especially in the early phases of drug development and of poorly soluble drugs.

Contribution of Extrahepatic Aldehyde Oxidase Activity to Human Clearance

Aldehyde oxidase (AOX) is a soluble, cytosolic enzyme that metabolizes various N-heterocyclic compounds and organic aldehydes. It has wide tissue distribution with highest levels found in liver, kidney, and lung. Human clearance projections of AOX substrates by in vitro assessments in isolated liver fractions (cytosol, S9) and even hepatocytes have been largely underpredictive of clinical outcomes. Various hypotheses have been suggested as to why this is the case. One explanation is that extrahepatic AOX expression contributes measurably to AOX clearance and is at least partially responsible for the often observed underpredictions. Although AOX expression has been confirmed in several extrahepatic tissues, activities therein and potential contribution to overall human clearance have not been thoroughly studied. In this work, the AOX enzyme activity using the S9 fractions of select extrahepatic human tissues (kidney, lung, vasculature, and intestine) were measured using carbazeran as a probe substrate. Measured activities were scaled to a whole-body clearance using best-available parameters and compared with liver S9 fractions. Here, the combined scaled AOX clearance obtained from the kidney, lung, vasculature, and intestine is very low and amounted to <1% of liver. This work suggests that AOX metabolism from extrahepatic sources plays little role in the underprediction of activity in human. One of the notable outcomes of this work has been the first direct demonstration of AOX activity in human vasculature. SIGNIFICANCE STATEMENT: This work demonstrates aldehyde oxidase (AOX) activity is measurable in a variety of extrahepatic human tissues, including vasculature, yet activities and potential contributions to human clearance are relatively low and insignificant when compared with the liver. Additionally, the modeling of the tissue-specific in vitro kinetic data suggests that AOX may be influenced by the tissue it resides in and thus show different affinity, activity, and modified activity over time.

Investigating the glucuronidation and sulfation pathways contribution and disposition kinetics of Bisphenol S and its metabolites using LC-MS/MS-based nonenzymatic hydrolysis method

Despite showing serious health consequences and widespread exposure, the toxicokinetic information required to evaluate the health risks of BPS is insufficient. Thus, we aim to describe the comprehensive toxicokinetics of BPS and its glucuronide (BPS-G) and sulfate (BPS-S) metabolites in rats. Simultaneous quantification of BPS and its metabolites (authentic standards) was accomplished using UPLC-MS/MS method. BPS displayed rapid absorption, extensive metabolism and fast elimination after oral administration. Following intravenous administration, BPS exhibited CL (8.8 L/h/kg) higher than the rat hepatic blood flow rate suggesting the likelihood of extrahepatic clearance. The CL value differed from those reported previously (sheep and piglets) and the probable reason could be attributed to dose- and/or interspecies differences. BPS was extensively metabolized and excreted primarily through urine as BPS-G (∼56%). BPS and BPS-S exhibited a high protein binding capacity in comparison to BPS-G. In in vitro metabolic stability study, BPS was predominantly metabolized through glucuronidation. The predicted in vivo hepatic clearance of BPS suggested it to be a high and intermediate clearance chemical in rats and humans, respectively. The significant interspecies difference observed in the clearance of BPS between rats and humans indicated that toxicokinetics of BPS should be considered for health risk assessment in humans.

Hepatic Scaling Factors for In Vitro-In Vivo Extrapolation of Metabolic Drug Clearance in Patients with Colorectal Cancer with Liver Metastasis

In vitro-in vivo extrapolation (IVIVE) linked with physiologically based pharmacokinetics (PBPK) modeling is used to predict the fates of drugs in patients. Ideally, the IVIVE-PBPK models should incorporate systems information accounting for characteristics of the specific target population. There is a paucity of such scaling factors in cancer, particularly microsomal protein per gram of liver (MPPGL) and cytosolic protein per gram of liver (CPPGL). In this study, cancerous and histologically normal liver tissue from 16 patients with colorectal liver metastasis were fractionated to microsomes and cytosol. Protein content was measured in homogenates, microsomes, and cytosol. The loss of microsomal protein during fractionation was accounted for using corrections based on NADPH cytochrome P450 reductase activity in different matrices. MPPGL was significantly lower in cancerous tissue (24.8 ± 9.8 mg/g) than histologically normal tissue (39.0 ± 13.8 mg/g). CPPGL in cancerous tissue was 42.1 ± 12.9 mg/g compared with 56.2 ± 16.9 mg/g in normal tissue. No correlations between demographics (sex, age, and body mass index) and MPPGL or CPPGL were apparent in the data. The generated scaling factors together with assumptions regarding the relative volumes of cancerous versus noncancerous tissue were used to simulate plasma exposure of drugs with different extraction ratios. The PBPK simulations revealed a substantial difference in drug exposure (area under the curve), up to 3.3-fold, when using typical scaling factors (healthy population) instead of disease-related parameters in cancer population. These indicate the importance of using population-specific scalars in IVIVE-PBPK for different disease states. SIGNIFICANCE STATEMENT: Accuracy in predicting the fate of drugs from in vitro data using IVIVE-PBPK depends on using correct scaling factors. The values for two of such scalars, namely microsomal and cytosolic protein per gram of liver, is not known in patients with cancer. This study presents, for the first time, scaling factors from cancerous and matched histologically normal livers. PBPK simulations of various metabolically cleared drugs demonstrate the necessity of population-specific scaling for model-informed precision dosing in oncology.

A simple method to measure sulfonation in man using paracetamol as probe drug

Sulfotransferase enzymes (SULT) catalyse sulfoconjugation of drugs, as well as endogenous mediators, gut microbiota metabolites and environmental xenobiotics. To address the limited evidence on sulfonation activity from clinical research, we developed a clinical metabolic phenotyping method using paracetamol as a probe substrate. Our aim was to estimate sulfonation capability of phenolic compounds and study its intraindividual variability in man. A total of 36 healthy adult volunteers (12 men, 12 women and 12 women on oral contraceptives) received paracetamol in a 1 g-tablet formulation on three separate occasions. Paracetamol and its metabolites were measured in plasma and spot urine samples using liquid chromatography-high resolution mass spectrometry. A metabolic ratio (Paracetamol Sulfonation Index—PSI) was used to estimate phenol SULT activity. PSI showed low intraindividual variability, with a good correlation between values in plasma and spot urine samples. Urinary PSI was independent of factors not related to SULT activity, such as urine pH or eGFR. Gender and oral contraceptive intake had no impact on PSI. Our SULT phenotyping method is a simple non-invasive procedure requiring urine spot samples, using the safe and convenient drug paracetamol as a probe substrate, and with low intraindividual coefficient of variation. Although it will not give us mechanistic information, it will provide us an empirical measure of an individual’s sulfonator status. To the best of our knowledge, our method provides the first standardised in vivo empirical measure of an individual’s phenol sulfonation capability and of its intraindividual variability. EUDRA-CT 2016-001395-29, NCT03182595 June 9, 2017.

Prediction of dose-dependent in vivo acetylcholinesterase inhibition by profenofos in rats and humans using physiologically based kinetic (PBK) modeling-facilitated reverse dosimetry

Organophosphate pesticides (OPs) are known to inhibit acetylcholine esterase (AChE), a critical effect used to establish health-based guidance values. This study developed a combined in vitro–in silico approach to predict AChE inhibition by the OP profenofos in rats and humans. A physiologically based kinetic (PBK) model was developed for both species. Parameter values for profenofos conversion to 4-bromo-2-chlorophenol (BCP) were derived from in vitro incubations with liver microsomes, liver cytosol, and plasma from rats (catalytic efficiencies of 1.1, 2.8, and 0.19 ml/min/mg protein, respectively) and humans (catalytic efficiencies of 0.17, 0.79, and 0.063 ml/min/mg protein, respectively), whereas other chemical-related parameter values were derived using in silico calculations. The rat PBK model was evaluated against literature data on urinary excretion of conjugated BCP. Concentration-dependent inhibition of rat and human AChE was determined in vitro and these data were translated with the PBK models to predicted dose-dependent AChE inhibition in rats and humans in vivo. Comparing predicted dose-dependent AChE inhibition in rats to literature data on profenofos-induced AChE inhibition revealed an accurate prediction of in vivo effect levels. Comparison of rat predictions (BMDL10 of predicted dose–response data of 0.45 mg/kg bw) and human predictions (BMDL10 of predicted dose–response data of 0.01 mg/kg bw) suggests that humans are more sensitive than rats, being mainly due to differences in kinetics. Altogether, the results demonstrate that in vivo AChE inhibition upon acute exposure to profenofos was closely predicted in rats, indicating the potential of this novel approach method in chemical hazard assessment.

Pinto J, Paudel A. Understanding Carrier Performance in Low-Dose Dry Powder Inhalation: An In Vitro-In Silico Approach

The use of physiologically based pharmacokinetic (PBPK) models to support drug product development has become increasingly popular. The in vitro characterization of the materials of the formulation provides valuable descriptors for the in silico prediction of the drug’s pharmacokinetic profile. Thus, the application of an in vitro–in silico framework can be decisive towards the prediction of the in vivo performance of a new medicine. By applying such an approach, this work aimed to derive mechanistic based insights into the potential impact of carrier particles and powder bulk properties on the in vivo performance of a lactose-based dry powder inhaler (DPI). For this, a PBPK model was developed using salbutamol sulphate (SS) as a model drug and the in vitro performance of its low-dose blends (2% w/w) with different types of lactose particles was investigated using different DPI types (capsule versus reservoir) at distinct airflows. Likewise, the influence of various carrier’s particle and bulk properties, device type and airflow were investigated in silico. Results showed that for the capsule-based device, low-dose blends of SS had a better performance, when smaller carrier particles (Dv_0.5 ≈ 50 μm) with about 10% of fines were used. This resulted in a better predicted bioavailability of the drug for all the tested airflows. For the reservoir type DPI, the mean particle size (Dv_0.5) was identified as the critical parameter impacting performance. Shear cell and air permeability or compressibility measurements, particle size distribution by pressure titration and the tensile strength of the selected lactose carrier powders were found useful to generate descriptors that could anticipate the potential in vivo performance of the tested DPI blends.

Physiologically-Based Pharmacokinetic-Pharmacodynamics Model Characterizing CYP2C19 Polymorphisms to Predict Clopidogrel Pharmacokinetics and Its Anti-Platelet Aggregation Effect Following Oral Administration to Coronary Artery Disease Patients With or Without Diabetes

Background and Objective: Clopidogrel (CLOP) is commonly used in coronary artery disease (CAD) patients with or without diabetes (DM), but these patients often suffer CLOP resistance, especially those with diabetes. This study was aimed to develop a physiologically-based pharmacokinetic-pharmacodynamic (PBPK-PD) model to describe the pharmacokinetics and pharmacodynamics of clopidogrel active metabolite (CLOP-AM) in CAD patients with or without DM.

Methods: The PBPK-PD model was first established and validated in healthy subjects and then in CAD patients with or without DM. The influences of CYP2C19, CYP2C9, CYP3A4, carboxylesterase 1 (CES1), gastrointestinal transit rates (K_t,i) and platelets response to CLOP-AM (k_irre) on predicted pharmacokinetics and pharmacodynamics were investigated, followed with their individual and integrated effects on CLOP-AM pharmacokinetics due to changes in DM status.

Results: Most predictions fell within 0.5–2.0 folds of observations, indicating successful predictions. Sensitivity analysis showed that contributions of interested factors to pharmacodynamics were CES1> k_irre> K_t,i> CYP2C19 > CYP3A4> CYP2C9. Mimicked analysis showed that the decreased exposure of CLOP-AM by DM was mainly attributed to increased CES1 activity, followed by decreased CYP2C19 activity.

Conclusion: The pharmacokinetics and pharmacodynamics of CLOP-AM were successfully predicted using the developed PBPK-PD model. Clopidogrel resistance by DM was the integrated effects of altered K_t,i, CYP2C19, CYP3A4, CES1 and k_irre.

Prediction of Drug Clearance from Enzyme and Transporter Kinetics

Accurate estimation of in vivo clearance in human is pivotal to determine the dose and dosing regimen for drug development. In vitro-in vivo extrapolation (IVIVE) has been performed to predict drug clearance using empirical and physiological scalars. Multiple in vitro systems and mathematical modeling techniques have been employed to estimate in vivo clearance. The models for predicting clearance have significantly improved and have evolved to become more complex by integrating multiple processes such as drug metabolism and transport as well as passive diffusion. This chapter covers the use of conventional as well as recently developed methods to predict metabolic and transporter-mediated clearance along with the advantages and disadvantages of using these methods and the associated experimental considerations. The general approaches to improve IVIVE by use of appropriate scalars, incorporation of extrahepatic metabolism and transport and application of physiologically based pharmacokinetic (PBPK) models with proteomics data are also discussed. The chapter also provides an overview of the advantages of using such dynamic mechanistic models over static models for clearance predictions to improve IVIVE.

Parameterization of Microsomal and Cytosolic Scaling Factors: Methodological and Biological Considerations for Scalar Derivation and Validation

Mathematical models that can predict the kinetics of compounds have been increasingly adopted for drug development and risk assessment. Data for these models may be generated from in vitro experimental systems containing enzymes contributing to metabolic clearance, such as subcellular tissue fractions including microsomes and cytosol. Extrapolation from these systems is facilitated by common scaling factors, known as microsomal protein per gram (MPPG) and cytosolic protein per gram (CPPG). Historically, parameterization of MPPG and CPPG has employed the use of recovery factors, commonly benchmarked to cytochromes P450 which work well in some contexts, but could be problematic for other enzymes. Here, we propose absolute quantification of protein content and supplementary assays to evaluate microsomal/cytosolic purity that should be employed. Examples include calculation of microsomal latency by mannose-6-phosphatase activity and immunoblotting of subcellular fractions with fraction-specific markers. Further considerations include tissue source, as disease states can affect enzyme expression and activity, and the methodology used for scalar parameterization. Regional- and organ-specific expression of enzymes, in addition to differences in organ physiology, is another important consideration. Because most efforts have focused on the liver that is, for the most part, homogeneous, derived scalars may not capture the heterogeneity of other major tissues contributing to xenobiotic metabolism including the kidneys and small intestine. Better understanding of these scalars, and how to appropriately derive them from extrahepatic tissues can provide support to the inferences made with physiologically based pharmacokinetic modeling, increase its accuracy in characterizing in vivo drug pharmacokinetics, and improve confidence in go-no-go decisions for clinical trials.

Main contribution of UGT1A1 and CYP2C9 in the metabolism of UR-1102, a novel agent for the treatment of gout

UR-1102, a novel uricosuric agent for treating gout, has been confirmed to exhibit a pharmacological effect in patients. We clarified its metabolic pathway, estimated the contribution of each metabolic enzyme, and assessed the impact of genetic polymorphisms using human in vitro materials. Glucuronide, sulfate and oxidative metabolites of UR-1102 were detected in human hepatocytes. The intrinsic clearance by glucuronidation or oxidation in human liver microsomes was comparable, but sulfation in the cytosol was much lower, indicating that the rank order of contribution was glucuronidation ≥ oxidation > sulfation. Recombinant UGT1A1 and UGT1A3 showed high glucuronidation of UR-1102. We took advantage of a difference in the inhibitory sensitivity of atazanavir to the UGT isoforms and estimated the fraction metabolised (fm) with UGT1A1 to be 70%. Studies using recombinant CYPs and CYP isoform-specific inhibitors showed that oxidation was mediated exclusively by CYP2C9. The effect of UGT1A1 and CYP2C9 inhibitors on UR-1102 metabolism in hepatocytes did not differ markedly between the wild type and variants.

A Translational Physiologically Based Pharmacokinetics/Pharmacodynamics Framework of Target-Mediated Disposition, Target Inhibition and Drug-Drug Interactions of Bortezomib

Bortezomib is a potent 20S proteasome inhibitor approved for the treatment of multiple myeloma and mantle cell lymphoma. Despite the extensive clinical use of bortezomib, the mechanism of the complex time-dependent pharmacokinetics of bortezomib has not been fully investigated in context of its pharmacodynamics (PD) and drug-drug interaction (DDI) profiles. Here, we aimed to develop a mechanistic physiologically based (PB) PK/PD model to project PK, blood target inhibition and DDI of bortezomib in patients. A minimal PBPK/PD model consisting of six compartments was constructed using a bottom-up approach with pre-clinical data and human physiological parameters. Specifically, the target-mediated drug disposition (TMDD) of bortezomib in red blood cells (RBC), which determines target inhibition in blood, was characterized by incorporating the proteasome binding affinity of bortezomib and the proteasome concentration in RBC. The hepatic clearance and fraction metabolized by different CYP isoforms were estimated from in vitro metabolism and phenotyping experiments. The established model adequately characterized the multi-exponential and time-dependent plasma pharmacokinetics, target binding and blood proteasome inhibition of bortezomib. Further, the model was able to accurately predict the impact of a strong CYP3A inducer (rifampicin) and inhibitor (ketoconazole) on bortezomib exposure. In conclusion, the mechanistic PBPK/PD model successfully described the complex pharmacokinetics, target inhibition and DDIs of bortezomib in patients. This study illustrates the importance of incorporating target biology, drug-target interactions and in vitro clearance parameters into mechanistic PBPK/PD models and the utility of such models for pharmacokinetic, pharmacodynamic and DDI predictions.

A matrix dispersion based on phospholipid complex system: preparation, lymphatic transport, and pharmacokinetics

Raloxifene hydrochloride (RH) suffers from low oral bioavailability due to its low water-solubility and first-pass metabolism. Therefore, a novel phospholipid complex of RH (RHPC) and a matrix dispersion based on phospholipid complex (RHPC-MD) were successfully prepared and optimized. Several methods were used to validate the formation of RHPC and RHPC-MD, such as differential scanning calorimetry, X-ray diffraction, scanning electron microscopy, transmission electron microscopy, infrared spectroscopy, particle size, and zeta potential, meanwhile, their octanol-water partition coefficient, solubility, and dissolution in vitro were also evaluated. To investigate the absorption mechanism of RHPC in vivo, the RHPC was administered to the chylomicron flow blockage rat model. Interestingly, as we expected, a significant reduction in RHPC absorption (67%) (**p< .01) in presence of cycloheximide (CXI) inhibitor was observed, thus confirming the RHPC could be absorbed by lymphatic transport in vivo. Pharmacokinetic studies revealed that the relative oral bioavailability of RHPC as well as RHPC-MD was 223% and 329%, respectively, when comparing with the commercial RH tablets. These outcomes suggested that the current study provided an attractive formulation to enhance the oral bioavailability of RH and stimulated to further research the absorption mechanism of RHPC in vivo.

Use of Physiologically Based Kinetic Modeling to Predict Rat Gut Microbial Metabolism of the Isoflavone Daidzein to S-Equol and Its Consequences for ERα Activation

SCOPE

To predict gut microbial metabolism of xenobiotics and the resulting plasma concentrations of metabolites formed, an in vitro-in silico-based testing strategy is developed using the isoflavone daidzein and its gut microbial metabolite S-equol as model compounds.

METHODS AND RESULTS

Anaerobic rat fecal incubations are optimized and performed to derive the apparent maximum velocities (Vmax ) and Michaelis-Menten constants (Km ) for gut microbial conversion of daidzein to dihydrodaidzein, S-equol, and O-desmethylangolensin, which are input as parameters for a physiologically based kinetic (PBK) model. The inclusion of gut microbiota in the PBK model allows prediction of S-equol concentrations and slightly reduced predicted maximal daidzein concentrations from 2.19 to 2.16 µm. The resulting predicted concentrations of daidzein and S-equol are comparable to in vivo concentrations reported.

CONCLUSION

The optimized in vitro approach to quantify kinetics for gut microbial conversions, and the newly developed PBK model for rats that includes gut microbial metabolism, provide a unique tool to predict the in vivo consequences of daidzein microbial metabolism for systemic exposure of the host to daidzein and its metabolite S-equol. The predictions reveal a dominant role for daidzein in ERα-mediated estrogenicity despite the higher estrogenic potency of its microbial metabolite S-equol.

Integrating in vitro testing and physiologically-based pharmacokinetic (PBPK) modelling for chemical liver toxicity assessment-A case study of troglitazone

In vitro to in vivo extrapolation (IVIVE) for next-generation risk assessment (NGRA) of chemicals requires computational modeling and faces unique challenges. Using mitochondria-related toxicity data of troglitazone (TGZ), a prototype drug known for liver toxicity, from HepaRG, HepG2, HC-04, and primary human hepatocytes, we explored inherent uncertainties in IVIVE, including cell models, cellular response endpoints, and dose metrics. A human population physiologically-based pharmacokinetic (PBPK) model for TGZ was developed to predict in vivo doses from in vitro point-of-departure (POD) concentrations. Compared to the 200-800 mg/d dose range of TGZ where liver injury was observed clinically, the predicted POD doses for the mean and top one percentile of the PBPK population were 28-372 and 15-178 mg/d respectively based on C_max dosimetry, and 185-2552 and 83-1010 mg/d respectively based on AUC. In conclusion, although with many uncertainties, integrating in vitro assays and PBPK modeling is promising in informing liver toxicity-inducing TGZ doses.

Biotransformation of Piceatannol, a Dietary Resveratrol Derivative: Promises to Human Health

SCOPE

To evaluate the health-promoting potentials of piceatannol (PIC), a dietary resveratrol derivative, its biotransformation is examined.

METHODS AND RESULTS

The biotransformation is tested in human/rat hepatic microsomes and cytosols; its pharmacokinetic profiles are assessed in rats. Although limited phase I metabolism exists in microsomes, PIC is rapidly converted to two pharmacologically active metabolites, namely rhapontigenin (RHA) and isorhapontigenin (ISO) in cytosols. Such biotransformation is completely blocked by entacapone, a well-known catechol-O-methyltransferase (COMT) inhibitor, demonstrating that the O-methylation is mediated by COMT. Moreover, PIC is identified as a substrate inhibitor of COMT, suggesting its potential benefits in Alzheimer's disease. Due to extensive phase II metabolism including glucuronidation, sulfation, and O-methylation, PIC displays rapid clearance and at least 4.02% ± 0.61% and 17.70% ± 0.91% of PIC is converted to RHA and ISO, respectively, in rats after intravenous administration. Similarly, PIC serves as an effective precursor of ISO upon oral administration.

CONCLUSION

Since PIC and its metabolites possess pleiotropic health-promoting activities, it has emerged as a promising nutraceutical candidate for further development. This study also reinforces the importance of in vivo testing in nutritional researches as the active metabolite(s) may be absent from the in vitro system.

Ontogeny of Hepatic Sulfotransferases and Prediction of Age-Dependent Fractional Contribution of Sulfation in Acetaminophen Metabolism

Cytosolic sulfotransferases (SULTs), including SULT1A, SULT1B, SULT1E, and SULT2A isoforms, play noteworthy roles in xenobiotic and endobiotic metabolism. We quantified the protein abundances of SULT1A1, SULT1A3, SULT1B1, and SULT2A1 in human liver cytosol samples (n = 194) by liquid chromatography-tandem mass spectrometry proteomics. The data were analyzed for their associations by age, sex, genotype, and ethnicity of the donors. SULT1A1, SULT1B1, and SULT2A1 showed significant age-dependent protein abundance, whereas SULT1A3 was invariable across 0-70 years. The respective mean abundances of SULT1A1, SULT1B1, and SULT2A1 in neonatal samples was 24%, 19%, and 38% of the adult levels. Interestingly, unlike UDP-glucuronosyltransferases and cytochrome P450 enzymes, SULT1A1 and SULT2A1 showed the highest abundance during early childhood (1 to <6 years), which gradually decreased by approx. 40% in adolescents and adults. SULT1A3 and SULT1B1 abundances were significantly lower in African Americans compared with Caucasians. Multiple linear regression analysis further confirmed the association of SULT abundances by age, ethnicity, and genotype. To demonstrate clinical application of the characteristic SULT ontogeny profiles, we developed and validated a proteomics-informed physiologically based pharmacokinetic model of acetaminophen. The latter confirmed the higher fractional contribution of sulfation over glucuronidation in the metabolism of acetaminophen in children. The study thus highlights that the ontogeny-based age-dependent fractional contribution (fm) of individual drug-metabolizing enzymes has better potential in prediction of drug-drug interactions and the effect of genetic polymorphisms in the pediatric population.

Use of generally recognized as safe or dietary compounds to inhibit buprenorphine metabolism: potential to improve buprenorphine oral bioavailability

The present study evaluated the potential of five generally recognized as safe (GRAS) or dietary compounds (α-mangostin, chrysin, ginger extract, pterostilbene and silybin) to inhibit oxidative (CYP) and conjugative (UGT) metabolism using pooled human intestinal and liver microsomes. Buprenorphine was chosen as the model substrate as it is extensively metabolized by CYPs to norbuprenorphine and by UGTs to buprenorphine glucuronide. Chrysin, ginger extract, α-mangostin, pterostilbene and silybin were tested for their inhibition of the formation of norbuprenorphine or buprenorphine glucuronide in both intestinal and liver microsomes. Pterostilbene was the most potent inhibitor of norbuprenorphine formation in both intestinal and liver microsomes, with IC50 values of 1.3 and 0.8 μM, respectively, while α-mangostin and silybin most potently inhibited buprenorphine glucuronide formation. The equipotent combination of pterostilbene and ginger extract additively inhibited both pathways in intestinal microsomes. Since pterostilbene and ginger extract showed potent CYP and/or UGT inhibition of buprenorphine metabolism, their equipotent combination was tested to assess the presence of synergistic inhibition. However, because the combination showed additive inhibition, it was not used while performing IVIVE analysis. Based on quantitative in vitro-in vivo extrapolation, pterostilbene (21 mg oral dose) appeared to be most effective in improving the mean predicted Foral and AUC∞ PO of buprenorphine from 3 ± 2% and 340 ± 330 ng*min/ml to 75 ± 8% and 36,000 ± 25,000 ng*min/ml, respectively. At a 10-fold lower dose of pterostilbene, the predicted buprenorphine Foral approximated sublingual bioavailability (~35%) and showed a 2-4 fold reduction in the variability around the predicted AUC∞ PO of buprenorphine. These results demonstrate the feasibility of using various GRAS/dietary compounds to inhibit substantially the metabolism by CYP and UGT enzymes to achieve higher and less variable oral bioavailability. This inhibitor strategy may be useful for drugs suffering from low and variable oral bioavailability due to extensive presystemic oxidative and/or conjugative metabolism.

[Further investigation of 3D culture spheroid models of human hepatocytes]

Three-dimensional (3D) cultured hepatocyte capable of maintaining liver-specific function in an in vivo state over a relatively long period of time have drawn attention as a new method for evaluating the metabolic process, hepatotoxicity and enzyme induction potential of drugs. When human hepatocytes were seeded on a plate for spheroid formation, and cell morphology and albumin secretion were examined, hepatocyte spheroid was stably maintained for at least 21 days after seeding. As a result of drug exposure to this spheroid, sequential metabolic reactions by Phase I and Phase II enzymes and metabolic reactions peculiar to only humans were observed. Moreover, when several drugs were exposed to spheroids and hepatotoxicity was evaluated, stable values were obtained for the 50% inhibitory concentration (IC₅₀) of albumin secretion at 14 and 21 days. The IC₅₀ values of most of the tested drugs were lower than in conventional assays, suggesting that the reported evaluation methods might underestimate hepatotoxicity. Furthermore, examination of mRNA expression level and activity of various cytochrome P450 (CYP) after exposure of typical inducers of CYPs to hepatocyte spheroid resulted in a significant increase in the expression level and activity of each. From these results, it was shown that this 3D hepatocyte spheroid system is suitable for follow-up of metabolic processes, long-term tests of hepatotoxicity and enzyme activity induction potential of drugs.

Targeted LC-MS/MS Proteomics-Based Strategy To Characterize in Vitro Models Used in Drug Metabolism and Transport Studies

Subcellular fractionation of tissue homogenate provides enriched in vitro models (e.g., microsomes, cytosol, or membranes), which are routinely used in the drug metabolism or transporter activity and protein abundance studies. However, batch-to-batch or interlaboratory variability in the recovery, enrichment, and purity of the subcellular fractions can affect performance of in vitro models leading to inaccurate in vitro to in vivo extrapolation (IVIVE) of drug clearance. To evaluate the quality of subcellular fractions, we developed a simple, targeted, and sensitive LC-MS/MS proteomics-based strategy, which relies on determination of protein markers of various cellular organelles, i.e., plasma membrane, cytosol, nuclei, mitochondria, endoplasmic reticulum (ER), lysosomes, peroxisomes, cytoskeleton, and exosomes. Application of the quantitative proteomics method confirmed a significant effect of processing variables (i.e., homogenization method and centrifugation speed) on the recovery, enrichment, and purity of isolated proteins in microsomes and cytosol. Particularly, markers of endoplasmic reticulum lumen and mitochondrial lumen were enriched in the cytosolic fractions as a result of their release during homogenization. Similarly, the relative recovery and composition of the total membrane fraction isolated from cell vs tissue samples was quantitatively different and should be considered in IVIVE. Further, analysis of exosomes isolated from sandwich-cultured hepatocyte media showed the effect of culture duration on compositions of purified exosomes. Therefore, the quantitative proteomics-based strategy developed here can be applied for efficient and simultaneous determination of multiple protein markers of various cellular organelles when compared to antibody- or activity-based assays and can be used for quality control of subcellular fractionation procedures including in vitro model development for drug metabolism and transport studies.

A Strategy to Refine the Phenotyping Approach and Its Implementation to Predict Drug Clearance: A Physiologically Based Pharmacokinetic Simulation Study

The phenotyping approach to predict drug metabolism activity is often hampered by a lack of correlation between the probe and the drug of interest. In this article, we present a strategy to refine the phenotyping approach based on a physiologically based pharmacokinetic simulation (implemented in Simcyp Simulator version 17) using previously published models. The apparent clearance (CL/F) of erlotinib was better predicted by the sum of caffeine and i.v. midazolam CL/F (r²  = 0.60) compared to that of either probe drug alone. The clearance of atorvastatin and repaglinide had a strong correlation (r²  = 0.70 and 0.63, respectively) with that of pitavastatin (a SLCO1B1 probe). Use of multiple probes for drugs that are predominantly metabolized by more than one cytochrome P450 (CYP) enzyme should be considered. In a case in which hepatic uptake transporters play a significant role in the disposition of a drug, the pharmacokinetic of a transporter probe will provide better predictions of the drug clearance.

Assessment of Transporter-Mediated and Passive Hepatic Uptake Clearance Using Rifamycin-SV as a Pan-Inhibitor of Active Uptake

The use of in vitro data for the quantitative prediction of transporter-mediated clearance is critical. Central to this evaluation is the use of hepatocytes, since they contain the full complement of transporters and metabolic enzymes. In general, uptake clearance (CL_uptake) is evaluated by measuring the appearance of compound in the cell. Passive clearance (CL_pd) is often determined by conducting parallel studies at 4 °C or by attempting to saturate uptake pathways. Both approaches have their limitations. Recent studies have proposed the use of Rifamycin-SV (RFV) as a pan-inhibitor of hepatic uptake pathways. In our studies, we confirm that transport activity of all major hepatic uptake transporters is inhibited significantly by RFV at 1 mM (OATP1B1, 1B3, and 2B1 = NTCP (80%), OCT1 (65%), OAT2 (60%)). Under these incubation conditions, we found that the free intracellular concentration of RFV is ∼175 μM and that several major CYPs and UGTs can be reversibly inhibited. Using this approach, we also determined CL_uptake and CL_pd of nine known OATP substrates across three different lots of human hepatocytes. The scaling factors generated for these compounds at 37 °C with RFV and 4 °C were found to be similar. The CL_pd of passively permeable compounds like metoprolol and semagacestat were found to be higher at 37 °C compared to 4 °C, indicating a temperature effect on these compounds. In addition, our data also suggests that incorporation of medium concentrations into CL_uptake and CL_pd calculations may be critical for highly protein bound and highly lipophilic drugs. Overall, our data indicate that RFV, instead of 4 °C, can be reliably used to measure CL_uptake and CL_pd of drugs.

A Semi-Physiologically Based Pharmacokinetic Model Describing the Altered Metabolism of Midazolam Due to Inflammation in Mice

PURPOSE

To investigate influence of inflammation on metabolism and pharmacokinetics (PK) of midazolam (MDZ) and construct a semi-physiologically based pharmacokinetic (PBPK) model to predict PK in mice with inflammatory disease.

METHODS

Glucose-6-phosphate isomerase (GPI)-mediated inflammation was used as a preclinical model of arthritis in DBA/1 mice. CYP3A substrate MDZ was selected to study changes in metabolism and PK during the inflammation. The semi-PBPK model was constructed using mouse physiological parameters, liver microsome metabolism, and healthy animal PK data. In addition, serum cytokine, and liver-CYP (cytochrome P450 enzymes) mRNA levels were examined.

RESULTS

The in vitro metabolite formation rate was suppressed in liver microsomes prepared from the GPI-treated mice as compared to the healthy mice. Further, clearance of MDZ was reduced during inflammation as compared to the healthy group. Finally, the semi-PBPK model was used to predict PK of MDZ after GPI-mediated inflammation. IL-6 and TNF-α levels were elevated and liver-cyp3a11 mRNA was reduced after GPI treatment.

CONCLUSION

The semi-PBPK model successfully predicted PK parameters of MDZ in the disease state. The model may be applied to predict PK of other drugs under disease conditions using healthy animal PK and liver microsomal data as inputs.

The Effect of Piperine Pro-Nano Lipospheres on Direct Intestinal Phase II Metabolism: The Raloxifene Paradigm of Enhanced Oral Bioavailability

Phase II biotransformation reactions have been gaining more attention due to their acknowledged significance in drug bioavailability, drug development, and drug-drug interactions. However, the predominant role of phase I metabolism has always overshadowed phase II metabolism, resulting in insufficient data regarding its mechanisms. In this paper, we investigate the effect of an advanced lipid based formulation on the phase II metabolism process of glucuronidation, occuring in the enterocytes monolayer. The investigated formulation is a self-emulsifying drug delivery system, termed pro-nano lipospheres, which contains the natural absorption enhancer piperine. To evaluate the effect of this formulation on direct glucuronidation we chose the model molecule raloxifene. First, glucuronidation is the main clearance pathway of this compound without involvement of preceding mechanisms. Second, raloxifene's extensive glucuronidation site is primarily at the intestine. Raloxifene's oral bioavailability was determined in a series of pharmacokinetic experiments using the freely moving rat model. In order to test the effect of the formulation on the relevant UGT enzymes reported in the clinic, we used the in vitro method of UGT-Glo Assay. Coadministration of raloxifene and piperine pro-nano lipospheres to rats resulted in a 2-fold increase in the relative oral bioavailability of raloxifene. However, coadministration of raloxifene with blank pro-nano lipospheres had no effect on its oral bioavailability. In contrast to the difference found in vivo between the two vehicles, both formulations extended an inhibitory effect on UGT enzymes in vitro. Ultimately, these findings prove the ability of the formulation to diminish intestinal direct phase II metabolism which serves as an absorption obstacle for many of today's marketed drugs. Pro-nano lipospheres is a formulation that serves as a platform for the simultaneous delivery of the absorption enhancer and a required drug. The discrepancy found between the in vivo and in vitro models demonstrates that the in vitro method may not be sensitive enough to distinguish the difference between the formulations.

Identification of Intestinal UDP-Glucuronosyltransferase Inhibitors in Green Tea (Camellia sinensis) Using a Biochemometric Approach: Application to Raloxifene as a Test Drug via In Vitro to In Vivo Extrapolation

Green tea (Camellia sinensis) is a popular beverage worldwide, raising concern for adverse interactions when co-consumed with conventional drugs. Like many botanical natural products, green tea contains numerous polyphenolic constituents that undergo extensive glucuronidation. As such, the UDP-glucuronosyltransferases (UGTs), particularly intestinal UGTs, represent potential first-pass targets for green tea-drug interactions. Candidate intestinal UGT inhibitors were identified using a biochemometrics approach, which combines bioassay and chemometric data. Extracts and fractions prepared from four widely consumed teas were screened (20-180 μg/ml) as inhibitors of UGT activity (4-methylumbelliferone glucuronidation) in human intestinal microsomes; all demonstrated concentration-dependent inhibition. A biochemometrics-identified fraction rich in UGT inhibitors from a representative tea was purified further and subjected to second-stage biochemometric analysis. Five catechins were identified as major constituents in the bioactive subfractions and prioritized for further evaluation. Of these catechins, (-)-epicatechin gallate and (-)-epigallocatechin gallate showed concentration-dependent inhibition, with IC₅₀ values (105 and 59 μM, respectively) near or below concentrations measured in a cup (240 ml) of tea (66 and 240 μM, respectively). Using the clinical intestinal UGT substrate raloxifene, the K_i values were ∼1.0 and 2.0 μM, respectively. Using estimated intestinal lumen and enterocyte inhibitor concentrations, a mechanistic static model predicted green tea to increase the raloxifene plasma area under the curve up to 6.1- and 1.3-fold, respectively. Application of this novel approach, which combines biochemometrics with in vitro-in vivo extrapolation, to other natural product-drug combinations will refine these procedures, informing the need for further evaluation via dynamic modeling and clinical testing.

Metabolic characteristics of Tanshinone I in human liver microsomes and S9 subcellular fractions

Tanshinone I (TSI) is a lipophilic diterpene in Salvia miltiorrhiza with versatile pharmacological activities. However, metabolic pathway of TSI in human is unknown. In this study, we determined major metabolites of TSI using a preparation of human liver microsomes (HLMs) by HPLC-UV and Q-Trap mass spectrometer. A total of 6 metabolites were detected, which indicated the presence of hydroxylation, reduction as well as glucuronidation. Selective chemical inhibition and purified cytochrome P450 (CYP450) isoform screening experiments revealed that CYP2A6 was primarily responsible for TSI Phase I metabolism. Part of generated hydroxylated TSI was glucuronidated via several glucuronosyltransferase (UGT) isoforms including UGT1A1, UGT1A3, UGT1A7, UGT1A9, as well as extrahepatic expressed isoforms UGT1A8 and UGT1A10. TSI could be reduced to a relatively unstable hydroquinone intermediate by NAD(P)H: quinone oxidoreductase 1 (NQO1), and then immediately conjugated with glucuronic acid by a panel of UGTs, especially UGT1A9, UGT1A1 and UGT1A8. Additionally, NQO1 could also reduce hydroxylated TSI to a hydroquinone intermediate, which was immediately glucuronidated by UGT1A1. The study demonstrated that hydroxylation, reduction as well as glucuronidation were the major pathways for TSI biotransformation, and six metabolites generated by CYPs, NQO1 and UGTs were found in HLMs and S9 subcellular fractions.

Challenges and Opportunities for Increasing the Knowledge Base Related to Drug Biotransformation and Pharmacokinetics during Growth and Development

It is generally acknowledged that there is a need and role for informative pharmacokinetic models to improve predictions and simulation as well as individualization of drug therapy in pediatric populations of different ages and developmental stages. This special issue contains more than 20 papers responding to the challenge of providing new information on scaling factors, ontogeny functions for drug metabolizing enzymes and transporters, the mechanisms underlying the observed developmental trajectories for these gene products, age-dependent changes in physiologic processes affecting drug disposition in children, as well as in vitro and in vivo studies describing the relative contribution of ontogeny and genetic factors as sources of variability in drug disposition in children. Considered together, these contributions serve to illustrate some of the current limitations regarding sample availability, number, and quality, but also provide a framework that allows for the potential value of the results of a given study to be interpreted within the context of these limitations. Among the challenges for the future are improving our understanding of the mechanisms regulating age-dependent changes in factors influencing drug disposition and response, thereby facilitating generalization to systems lacking detailed data, better integrating age-dependent changes in pharmacokinetics with age-dependent changes in pharmacodynamics, and allowing better predictability and individualization of drug disposition and response across the pediatric age spectrum.