Quantification of Proteins Involved in Drug Metabolism and Disposition in the Human Liver Using Label-Free Global Proteomics

Transporter expressions as part of required scaling factor to support in vitro in vivo extrapolation for blood-brain barrier drug permeability

Access of drugs to the central nervous system is limited by the blood-brain barrier, and this in turn affects drug efficacy/toxicity. To date, most drug discovery optimization paradigms have relied heavily on in vitro transporter assays and preclinical species pharmacokinetic evaluation to provide a qualitative assessment of human brain penetration. Because of the lack of human brain pharmacokinetic data, mechanistic models for preclinical species, combined with in vitro and in silico data, are useful for translation to human. These models require transporter expression data to be measured in both in vitro and in vivo systems. The purpose of this work was to quantify transporter expression and generate scaling factors (SFs) to enable in vitro in vivo extrapolation (IVIVE) of transporter-mediated processes and to support the development of PBPK model of the brain in rats. SF represents the ratio of abundance of the relevant transporter in the tissue relative to transporter expressing cells. Using quantitative proteomics with QconCAT technology, the expression of human and rat P-gp (ABCB1/Abcb1) and BCRP/Bcrp (ABCG2/Abcg2) was measured in rat brain microvessels, mock and transfected cell lines including, Madin-Darby Canine Kidney I (MDCK I), Madin-Darby Canine Kidney II (MDCK II) and Porcine Kidney epithelial cells (LLC-PK1). P-gp expression ranged from 32 to 71 pmol/mg in rat brain microvessels, exceeding literature values of 14.1-25.2 pmol/mg protein. Conversely, Bcrp expression ranged between 0.02-0.27 pmol/mg protein lower than the literature range (2-6.2 pmol/mg protein). P-gp expression in MDCK I and LLC-PK1 cells transfected with rat Mdr1a was similar (within 1.5-fold) as was human P-gp expression in MDR1 transfected LLC-PK1 and MDCK II cells. The generated SFs were 34.4 and 50.4 for brain P-gp (depending on the cell line used) and 0.53 for brain Bcrp. Endogenous P-gp transporter was detected in MDCK II cell lines when protein expression was measured using a surrogate peptide that was shared across species. The current work provides a framework for proteomics-informed translation of in vitro P-gp and BCRP-related kinetics of drugs and supports the development of PBPK models to predict drug disposition in the brain.

Sex-specific UGT expression and function: prevalence, potential mechanisms and significance

INTRODUCTION

Sex and gender influence pharmacotherapy outcomes, including adverse drug effects which are nearly twice as common in women. Sex differences in drug responses involve factors as diverse as body composition, physiology, and prescribing patterns. Many drugs show higher exposure in women, which can be partly attributed to sex-differences in processes that control drug disposition such as metabolism and transport.

AREAS COVERED

This article reviews sex differences in the expression and function of the critical phase II drug-metabolizing enzymes, UDP-glucuronosyltransferases (UGTs). We curate the literature on sex-biased UGT expression in human tissues, describe the evidence for UGT-mediated sex-differences in drug exposure, and critically evaluate whether UGTs contribute to different drug outcomes in males and females. Relevant literature was identified by searching PubMed with terms including UDP-glucuronosyltransferase/UGT, glucuronidation/glucuronide, sex, gender, male, female, men, and women.

EXPERT OPINION

Several examples of sex-biased UGT expression and drug glucuronidation were identified; however, evidence of clinical impact was more limited. Significant data gaps limit our understanding of the prevalence and importance of sex-biased glucuronidation. Novel methodologies for tissue-level metabolite sampling together with increased sex-aware analysis of clinical/preclinical data, could help address gaps and reveal new avenues for enhancing pharmacotherapy outcomes for all genders.

Developmental neurotoxicity evaluation of acrylamide based on in vitro to in vivo extrapolation by pregnancy PBTK modelling

Acrylamide (ACR) is a known neurotoxicant that can pass the placenta and has been detected in breast milk. Some in vivo and in vitro studies indicate that ACR exposure might lead to developmental neurotoxicity (DNT). Here, we have developed a physiologically-based toxicokinetic model for a pregnant human population using PK-Sim. We performed an in vitro to in vivo extrapolation (IVIVE) of data collected from human neuroblastoma SH-SY5Y cells exposed during differentiation to ACR. The developed PBTK model was successfully evaluated and predicted fetal plasma concentrations in the low nM range after exposing the model to an estimated average daily intake for pregnant women. The IVIVE showed that low concentrations of ACR (fM-nM) that induced attenuated differentiation of the SH-SY5Y neuronal cell model, were relevant for human exposure to ACR from oral intake. However, doses estimated in the IVIVE from concentrations in the µM range, were found to be unrealistic by exposure through food intake for an average daily intake. However, in case of exposure due to environmental pollution or occupational exposure, these concentrations may be reached in fetal plasma. The findings in this study raise the concern regarding ACR exposure during pregnancy as well as the relevance of testing concentrations in vitro that are several orders of magnitude higher than the predicted fetal plasma concentrations.

Differential Tissue Abundance of Membrane-Bound Drug Metabolizing Enzymes and Transporter Proteins by Global Proteomics

Protein abundance data of drug-metabolizing enzymes and transporters (DMETs) are useful for scaling in vitro and animal data to humans for accurate prediction and interpretation of drug clearance and toxicity. Targeted DMET proteomics that relies on synthetic stable isotope-labeled surrogate peptides as calibrators is routinely used for the quantification of selected proteins; however, the technique is limited to the quantification of a small number of proteins. Although the global proteomics-based total protein approach (TPA) is emerging as a better alternative for large-scale protein quantification, the conventional TPA does not consider differential sequence coverage by identifying unique peptides across proteins. Here, we optimized the TPA approach by correcting protein abundance data by the sequence coverage, which was applied to quantify 54 DMETs for characterization of 1) differential tissue DMET abundance in the human liver, kidney, and intestine, and 2) interindividual variability of DMET proteins in individual intestinal samples (n = 13). Uridine diphosphate-glucuronosyltransferase 2B7 (UGT2B7), microsomal glutathione S-transferases (MGST1, MGST2, and MGST3) carboxylesterase 2 (CES2), and multidrug resistance-associated protein 2 (MRP2) were expressed in all three tissues, whereas, as expected, four cytochrome P450s (CYP3A4, CYP3A5, CYP2C9, and CYP4F2), UGT1A1, UGT2B17, CES1, flavin-containing monooxygenase 5, MRP3, and P-glycoprotein were present in the liver and intestine. The top three DMET proteins in individual tissues were: CES1>CYP2E1>UGT2B7 (liver), CES2>UGT2B17>CYP3A4 (intestine), and MGST1>UGT1A6>MGST2 (kidney). CYP3A4, CYP3A5, UGT2B17, CES2, and MGST2 showed high interindividual variability in the intestine. These data are relevant for enhancing in vitro to in vivo extrapolation of drug absorption and disposition and can be used to enhance the accuracy of physiologically based pharmacokinetic prediction of systemic and tissue concentration of drugs. SIGNIFICANCE STATEMENT: This study quantified the abundance and compositions of drug-metabolizing enzymes and transporters in pooled human liver, intestine, and kidney microsomes as well as individual intestinal microsomes using an optimized global proteomics approach. The data revealed large intertissue differences in the abundance of these proteins and high intestinal interindividual variability in the levels of cytochrome P450s (e.g., CYP3A4 and CYP3A5), uridine diphosphate-glucuronosyltransferase 2B17, carboxylesterase 2, and microsomal glutathione S-transferase 2. These data are applicable for the prediction of first-pass metabolism and tissue-specific drug clearance.

Complementarity of two proteomic data analysis tools in the identification of drug-metabolising enzymes and transporters in human liver

Several software packages are available for the analysis of proteomic LC-MS/MS data, including commercial (e.g. Mascot/Progenesis LC-MS) and open access software (e.g. MaxQuant). In this study, Progenesis and MaxQuant were used to analyse the same data set from human liver microsomes (n = 23). Comparison focussed on the total number of peptides and proteins identified by the two packages. For the peptides exclusively identified by each software package, distribution of peptide length, hydrophobicity, molecular weight, isoelectric point and score were compared. Using standard cut-off peptide scores, we found an average of only 65% overlap in detected peptides, with surprisingly little consistency in the characteristics of peptides exclusively detected by each package. Generally, MaxQuant detected more peptides than Progenesis, and the additional peptides were longer and had relatively lower scores. Progenesis-specific peptides tended to be more hydrophilic and basic relative to peptides detected only by MaxQuant. At the protein level, we focussed on drug-metabolising enzymes (DMEs) and transporters, by comparing the number of unique peptides detected by the two packages for these specific proteins of interest, and their abundance. The abundance of DMEs and SLC transporters showed good correlation between the two software tools, but ABC showed less consistency. In conclusion, in order to maximise the use of MS datasets, we recommend processing with more than one software package. Together, Progenesis and MaxQuant provided excellent coverage, with a core of common peptides identified in a very robust way.

Comparative Intra-Subject Analysis of Gene Expression and Protein Abundance of Major and Minor Drug Metabolizing Enzymes in Healthy Human Jejunum and Liver

First pass metabolism by phase I and phase II enzymes in the intestines and liver is a major determinant of the oral bioavailability of many drugs. Several studies analyzed expressions of major drug-metabolizing enzymes (DMEs), such as CYP3A4 and UGT1A1 in the human gut and liver. However, there is still a lack of knowledge regarding other DMEs (i.e., "minor" DMEs), although several clinically relevant drugs are affected by those enzymes. Moreover, there is very limited intra-subject data on hepatic and intestinal expression levels of minor DMEs. To fill this gap of knowledge, we analyzed gene expression (quantitative real-time polymerase chain reaction) and protein abundance (targeted proteomics) of 24 clinically relevant DMEs, that is, carboxylesterases (CES), UDP-glucuronosyltransferases (UGT), and cytochrome P450 (CYP)-enzymes. We performed our analysis using jejunum and liver tissue specimens from the same 11 healthy organ donors (8 men and 3 women, aged 19-60 years). Protein amounts of all investigated DMEs, with the exception of CYP4A11, were detected in human liver samples. CES2, CYP2C18, CYP3A4, and UGT2B17 protein abundance was similar or even higher in the jejunum, and all other DMEs were found in higher amounts in the liver. Significant correlations between gene expression and protein levels were observed only for 2 of 15 jejunal, but 13 of 23 hepatic DMEs. Intestinal and hepatic protein amounts only significantly correlated for CYP3A4 and UGT1A3. Our results demonstrated a notable variability between the individuals, which was even higher in the intestines than in the liver. Our intrasubject analysis of DMEs in the jejunum and liver from healthy donors, may be useful for physiologically-based pharmacokinetic-based modeling and prediction in order to improve efficacy and safety of oral drug therapy.

Comparison of 1Beta- and 5Beta-hydroxylation of Deoxycholate and Glycodeoxycholate as In Vitro Index Reactions for Cytochrome P450 3A Activities

Cytochrome P450 3A (CYP3A) participates in the metabolism of more than 30% of clinical drugs. The vast intra- and inter-individual variations in CYP3A activity pose great challenges to drug development and personalized medicine. It has been disclosed that human CYP3A4 and CYP3A7 are exclusively responsible for the tertiary oxidations of deoxycholic acid (DCA) and glycodeoxycholic acid (GDCA) regioselectivity at C-1β and C-5β This work aimed to compare the 1β- and 5β-hydroxylation of DCA and GDCA as potential in vitro CYP3A index reactions in both human liver microsomes and recombinant P450 enzymes. The results demonstrated that the metabolic activity of DCA 1β- and 5β-hydroxylation was 5-10 times higher than that of GDCA, suggesting that 1β-hydroxyglycodeoxycholic acid and 5β-hydroxyglycodeoxycholic acid may originate from DCA oxidation followed by conjugation in humans. Metabolic phenotyping data revealed that DCA 1β-hydroxylation, DCA 5β-hydroxylation, and GDCA 5β-hydroxylation were predominantly catalyzed by CYP3A4 (>80%), while GDCA 1β-hydroxylation had approximately equal contributions from CYP3A4 (41%) and 3A7 (58%). Robust Pearson correlation was established for the intrinsic clearance of DCA 1β- and 5β-hydroxylation with midazolam (MDZ) 1'- and 4-hydroxylation in fourteen single donor microsomes. Although DCA 5β-hydroxylation exhibited a stronger correlation with MDZ oxidation, DCA 1β-hydroxylation exhibited higher reactivity than DCA 5β-hydroxylation. It is therefore suggested that DCA 1β- and 5β-hydroxylations may serve as alternatives to T 6β-hydroxylation as in vitro CYP3A index reactions. SIGNIFICANCE STATEMENT: The oxidation of DCA and GDCA is primarily catalyzed by CYP3A4 and CYP3A7. This work compared the 1β- and 5β-hydroxylation of DCA and GDCA as in vitro index reactions to assess CYP3A activities. It was disclosed that the metabolic activity of DCA 1β- and 5β-hydroxylation was 5-10 times higher than that of GDCA. Although DCA 1β-hydroxylation exhibited higher metabolic activity than DCA 5β-hydroxylation, DCA 5β-hydroxylation demonstrated stronger correlation with MDZ oxidation than DCA 1β-hydroxylation in individual liver microsomes.

Tree shrew cytochrome P450 2E1 is a functional enzyme that metabolises chlorzoxazone and p-nitrophenol

Cytochromes P450 (CYPs or P450s) are important enzymes for drug metabolism. Tree shrews are non-primate animal species used in various fields of biomedical research, including infection (especially hepatitis viruses), depression, and myopia. A recent tree shrew genome analysis indicated that the sequences and the numbers of P450 genes are similar to those of humans; however, P450s have not been adequately identified and analysed in this species.In this study, a novel CYP2E1 was isolated from tree shrew liver and was characterised in comparison with human, dog, and pig CYP2E1. Tree shrew CYP2E1 and human CYP2E1 showed high amino acid sequence identity (83%) and were closely related in a phylogenetic tree.Gene and genome structures of CYP2E1 were generally similar in humans, dogs, pigs, and tree shrews. Tissue expression patterns showed that tree shrew CYP2E1 mRNA was predominantly expressed in liver, just as for dog and pig CYP2E1 mRNAs. In tree shrews, recombinant CYP2E1 protein and liver microsomes metabolised chlorzoxazone and p-nitrophenol, probe substrates of human CYP2E1, just as they do in dogs and pigs.These results suggest that tree shrew CYP2E1 encodes a functional drug-metabolising enzyme that plays a role in the liver, similar to human CYP2E1.

Quantification of drug metabolising enzymes and transporter proteins in the paediatric duodenum via LC-MS/MS proteomics using a QconCAT technique

Characterising the small intestine absorptive membrane is essential to enable prediction of the systemic exposure of oral formulations. In particular, the ontogeny of key intestinal Drug Metabolising Enzymes and Transporter (DMET) proteins involved in drug disposition needs to be elucidated to allow for accurate prediction of the PK profile of drugs in the paediatric cohort. Using pinch biopsies from the paediatric duodenum (n = 36; aged 11 months to 15 years), the abundance of 21 DMET proteins and two enterocyte markers were quantified via LC-MS/MS. An established LCMS nanoflow method was translated to enable analysis on a microflow LC system, and a new stable-isotope-labelled QconCAT standard developed to enable quantification of these proteins. Villin-1 was used to standardise abundancy values. The observed abundancies and ontogeny profiles, agreed with adult LC-MS/MS-based data, and historic paediatric data obtained via western blotting. A linear trend with age was observed for duodenal CYP3A4 and CES2 only. As this work quantified peptides on a pinch biopsy coupled with a microflow method, future studies using a wider population range are very feasible. Furthermore, this DMET ontogeny data can be used to inform paediatric PBPK modelling and to enhance the understanding of oral drug absorption and gut bioavailability in paediatric populations.

Addressing the Clinical Importance of Equilibrative Nucleoside Transporters in Drug Discovery and Development

The US Food and Drug Administration (FDA), European Medicines Agency (EMA), and Pharmaceuticals and Medical Devices Agency (PMDA) guidances on small-molecule drug-drug interactions (DDIs), with input from the International Transporter Consortium (ITC), recommend the evaluation of nine drug transporters. Although other clinically relevant drug uptake and efflux transporters have been discussed in ITC white papers, they have been excluded from further recommendation by the ITC and are not included in current regulatory guidances. These include the ubiquitously expressed equilibrative nucleoside transporters (ENT) 1 and ENT2, which have been recognized by the ITC for their potential role in clinically relevant nucleoside analog drug interactions for patients with cancer. Although there is comparatively limited clinical evidence supporting their role in DDI risk or other adverse drug reactions (ADRs) compared with the nine highlighted transporters, several in vitro and in vivo studies have identified ENT interactions with non-nucleoside/non-nucleotide drugs, in addition to nucleoside/nucleotide analogs. Some noteworthy examples of compounds that interact with ENTs include cannabidiol and selected protein kinase inhibitors, as well as the nucleoside analogs remdesivir, EIDD-1931, gemcitabine, and fialuridine. Consequently, DDIs involving the ENTs may be responsible for therapeutic inefficacy or off-target toxicity. Evidence suggests that ENT1 and ENT2 should be considered as transporters potentially involved in clinically relevant DDIs and ADRs, thereby warranting further investigation and regulatory consideration.

Evaluating the Pharmacokinetics of Fentanyl in the Brain Extracellular Fluid, Saliva, Urine, and Plasma of Newborns from Transplacental Exposure from Parturient Mothers Dosed with Epidural Fentanyl Utilizing PBPK Modeling

BACKGROUND AND OBJECTIVE

Fentanyl can mitigate the mother and newborn complications resulting from labor pain. However, fentanyl shows a narrow therapeutic index between its respiratory depressive and analgesic effects. Thus, prenatally acquired high fentanyl levels in the newborn brain extracellular fluid (bECF) may induce respiratory depression which requires therapeutic drug monitoring (TDM). TDM using saliva and urine in newborns can reduce the possibility of infections and distress associated with TDM using blood. The objective of this study was to develop  a physiologically based pharmacokinetic (PBPK) model to predict fentanyl concentrations in different newborn tissues due to intrauterine exposure.

METHODS

A fentanyl PBPK model in adults after intravenous and epidural administration was built, validated, and scaled to pregnancy and newborn populations. The dose that the newborn received transplacentally at birth was calculated using the pregnancy model. Then, the newborn bECF, saliva, plasma, and urine concentrations after such a dose were predicted using the newborn PBPK model.

RESULTS

After a maternal epidural dose of fentanyl 245 µg, the predicted newborn plasma and bECF levels were below the toxicity thresholds. Furthermore, the salivary threshold levels in newborns for fentanyl analgesic and respiratory depression effects were estimated to be 0.39 and 14.7-18.2 ng/ml, respectively.

CONCLUSION

The salivary TDM of fentanyl in newborns can be useful in newborns exposed to intrauterine exposure from parturient females dosed with epidural fentanyl. However, newborn-specific values of µ-opioid receptors IC50 for respiratory depression are needed.

CYP2E1 C-1054T and 96-bp I/D genetic variations and risk of gestational diabetes mellitus in chinese women: a case-control study

BACKGROUND

Cytochrome P450 2E1 (CYP2E1) plays a key role in the metabolism of xenobiotic and endogenous low-molecular-weight compounds. This study aimed to determine if the genetic variations of 96-bp insertion/deletion (I/D) and C-1054T (rs2031920) in CYP2E1 were associated with the risk of gestational diabetes mellitus (GDM).

METHODS

CYP2E1 polymorphisms were genotyped in a case-control study of 1,134 women with uncomplicated pregnancies and 723 women with GDM. The effects of genotype on the clinical, metabolic, and oxidative stress indices were assessed.

RESULTS

The CYP2E1 C-1054T variant was associated with an increased risk of GDM based on the genotype, recessive, dominant, and allele genetic models (P < 0.05). The TT + CT genotype remained a significant predictive factor for GDM risk after correcting for maternal age and pre-pregnancy body mass index (OR = 1.277, 95% CI: 1.042-1.563, P = 0.018). Moreover, fasting insulin concentrations and homeostatic model assessment of insulin resistance were significantly higher in GDM patients carrying the T allele than in those with the CC genotype (P < 0.05). Furthermore, the combined genotype II + ID/TT + CT of the 96-bp I/D and C-1054T polymorphisms further increased the risk of GDM when the combined genotype DD/CC was set as the reference category (OR = 1.676, 95% CI: 1.182-2.376, P = 0.004).

CONCLUSIONS

The T allele of the C-1054T polymorphism and its combination with the I allele of the 96-bp I/D variation in CYP2E1 are associated with an increased risk of GDM in the Chinese population. The - 1054T allele may be associated with more serious insulin resistance in patients.

Drug-drug interactions that alter the exposure of glucuronidated drugs: Scope, UDP-glucuronosyltransferase (UGT) enzyme selectivity, mechanisms (inhibition and induction), and clinical significance

Drug-drug interactions (DDIs) arising from the perturbation of drug metabolising enzyme activities represent both a clinical problem and a potential economic loss for the pharmaceutical industry. DDIs involving glucuronidated drugs have historically attracted little attention and there is a perception that interactions are of minor clinical relevance. This review critically examines the scope and aetiology of DDIs that result in altered exposure of glucuronidated drugs. Interaction mechanisms, namely inhibition and induction of UDP-glucuronosyltransferase (UGT) enzymes and the potential interplay with drug transporters, are reviewed in detail, as is the clinical significance of known DDIs. Altered victim drug exposure arising from modulation of UGT enzyme activities is relatively common and, notably, the incidence and importance of UGT induction as a DDI mechanism is greater than generally believed. Numerous DDIs are clinically relevant, resulting in either loss of efficacy or an increased risk of adverse effects, necessitating dose individualisation. Several generalisations relating to the likelihood of DDIs can be drawn from the known substrate and inhibitor selectivities of UGT enzymes, highlighting the importance of comprehensive reaction phenotyping studies at an early stage of drug development. Further, rigorous assessment of the DDI liability of new chemical entities that undergo glucuronidation to a significant extent has been recommended recently by regulatory guidance. Although evidence-based approaches exist for the in vitro characterisation of UGT enzyme inhibition and induction, the availability of drugs considered appropriate for use as 'probe' substrates in clinical DDI studies is limited and this should be research priority.

The SLIT/ROBO Pathway in Liver Fibrosis and Cancer

Liver fibrosis is a common outcome of most chronic liver insults/injuries that can develop into an irreversible process of cirrhosis and, eventually, liver cancer. In recent years, there has been significant progress in basic and clinical research on liver cancer, leading to the identification of various signaling pathways involved in tumorigenesis and disease progression. Slit glycoprotein (SLIT)1, SLIT2, and SLIT3 are secreted members of a protein family that accelerate positional interactions between cells and their environment during development. These proteins signal through Roundabout receptor (ROBO) receptors (ROBO1, ROBO2, ROBO3, and ROBO4) to achieve their cellular effects. The SLIT and ROBO signaling pathway acts as a neural targeting factor regulating axon guidance, neuronal migration, and axonal remnants in the nervous system. Recent findings suggest that various tumor cells differ in SLIT/ROBO signaling levels and show varying degrees of expression patterns during tumor angiogenesis, cell invasion, metastasis, and infiltration. Emerging roles of the SLIT and ROBO axon-guidance molecules have been discovered in liver fibrosis and cancer development. Herein, we examined the expression patterns of SLIT and ROBO proteins in normal adult livers and two types of liver cancers: hepatocellular carcinoma and cholangiocarcinoma. This review also summarizes the potential therapeutics of this pathway for anti-fibrosis and anti-cancer drug development.

The Analysis of Pethidine Pharmacokinetics in Newborn Saliva, Plasma, and Brain Extracellular Fluid After Prenatal Intrauterine Exposure from Pregnant Mothers Receiving Intramuscular Dose Using PBPK Modeling

BACKGROUND AND OBJECTIVE

Pethidine (meperidine) can decrease labor pain-associated mother's hyperventilation and high cortisol-induced newborn complications. However, prenatal transplacentally acquired pethidine can cause side effects in newborns. High pethidine concentrations in the newborn brain extracellular fluid (bECF) can cause a serotonin crisis. Therapeutic drug monitoring (TDM) in newborns' blood distresses them and increases infection incidence, which can be overcome by using salivary TDM. Physiologically based pharmacokinetic (PBPK) modeling can predict drug concentrations in newborn plasma, saliva, and bECF after intrauterine pethidine exposure.

METHODS

A healthy adult PBPK model was constructed, verified, and scaled to newborn and pregnant populations after intravenous and intramuscular pethidine administration. The pregnancy PBPK model was used to predict the newborn dose received transplacentally at birth, which was used as input to the newborn PBPK model to predict newborn plasma, saliva, and bECF pethidine concentrations and set correlation equations between them.

RESULTS

Pethidine can be classified as a Salivary Excretion Classification System class II drug. The developed PBPK model predicted that, after maternal pethidine intramuscular doses of 100 mg and 150 mg, the newborn plasma and bECF concentrations were below the toxicity thresholds. Moreover, it was estimated that newborn saliva concentrations of 4.7 µM, 11.4 µM, and 57.7 µM can be used as salivary threshold concentrations for pethidine analgesic effects, side effects, and the risk for serotonin crisis, respectively, in newborns.

CONCLUSION

It was shown that saliva can be used for pethidine TDM in newborns during the first few days after delivery to mothers receiving pethidine.

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.

Toward systems-informed models for biologics disposition: covariates of the abundance of the neonatal Fc Receptor (FcRn) in human tissues and implications for pharmacokinetic modelling

Biologics are a fast-growing therapeutic class, with intertwined pharmacokinetics and pharmacodynamics, affected by the abundance and function of the FcRn receptor. While many investigators assume adequacy of classical models, such as allometry, for pharmacokinetic characterization of biologics, advocates of physiologically-based pharmacokinetics (PBPK) propose consideration of known systems parameters that affect the fate of biologics to enable a priori predictions, which go beyond allometry. The aim of this study was to deploy a systems-informed modelling approach to predict the disposition of Fc-containing biologics. We used global proteomics to quantify the FcRn receptor [p51 and β₂-microglobulin (B2M) subunits] in 167 samples of human tissue (liver, intestine, kidney and skin) and assessed covariates of its expression. FcRn p51 subunit was highest in liver relative to other tissues, and B2M was 1-2 orders of magnitude more abundant than FcRn p51 across all sets. There were no sex-related differences, while higher expression was confirmed in neonate liver compared with adult liver. Trends of expression in liver and kidney indicated a moderate effect of body mass index, which should be confirmed in a larger sample size. Expression of FcRn p51 subunit was approximately 2-fold lower in histologically normal liver tissue adjacent to cancer compared with healthy liver. FcRn mRNA in plasma-derived exosomes correlated moderately with protein abundance in matching liver tissue, opening the possibility of use as a potential clinical tool. Predicted effects of trends in FcRn abundance in healthy and disease (cancer and psoriasis) populations using trastuzumab and efalizumab PBPK models were in line with clinical observations, and global sensitivity analysis revealed endogenous IgG plasma concentration and tissue FcRn abundance as key systems parameters influencing exposure to Fc-conjugated biologics.

The investigation of the complex population-drug-drug interaction between ritonavir-boosted lopinavir and chloroquine or ivermectin using physiologically-based pharmacokinetic modeling

OBJECTIVES

Therapy failure caused by complex population-drug-drug (PDDI) interactions including CYP3A4 can be predicted using mechanistic physiologically-based pharmacokinetic (PBPK) modeling. A synergy between ritonavir-boosted lopinavir (LPVr), ivermectin, and chloroquine was suggested to improve COVID-19 treatment. This work aimed to study the PDDI of the two CYP3A4 substrates (ivermectin and chloroquine) with LPVr in mild-to-moderate COVID-19 adults, geriatrics, and pregnancy populations.

METHODS

The PDDI of LPVr with ivermectin or chloroquine was investigated. Pearson's correlations between plasma, saliva, and lung interstitial fluid (ISF) levels were evaluated. Target site (lung epithelial lining fluid [ELF]) levels of ivermectin and chloroquine were estimated.

RESULTS

Upon LPVr coadministration, while the chloroquine plasma levels were reduced by 30, 40, and 20%, the ivermectin plasma levels were increased by a minimum of 425, 234, and 453% in adults, geriatrics, and pregnancy populations, respectively. The established correlation equations can be useful in therapeutic drug monitoring (TDM) and dosing regimen optimization.

CONCLUSIONS

Neither chloroquine nor ivermectin reached therapeutic ELF levels in the presence of LPVr despite reaching toxic ivermectin plasma levels. PBPK modeling, guided with TDM in saliva, can be advantageous to evaluate the probability of reaching therapeutic ELF levels in the presence of PDDI, especially in home-treated patients.

Protein Abundance of Drug Metabolizing Enzymes in Human Hepatitis C Livers

Hepatic drug metabolizing enzymes (DMEs), whose activity may be affected by liver diseases, are major determinants of drug pharmacokinetics. Hepatitis C liver samples in different functional states, i.e., the Child-Pugh class A (n = 30), B (n = 21) and C (n = 7) were analyzed for protein abundances (LC-MS/MS) and mRNA levels (qRT-PCR) of 9 CYPs and 4 UGTs enzymes. The protein levels of CYP1A1, CYP2B6, CYP2C8, CYP2C9, and CYP2D6 were not affected by the disease. In the Child-Pugh class A livers, a significant up-regulation of UGT1A1 (to 163% of the controls) was observed. The Child-Pugh class B was associated with down-regulation of the protein abundance of CYP2C19 (to 38% of the controls), CYP2E1 (to 54%), CYP3A4 (to 33%), UGT1A3 (to 69%), and UGT2B7 (to 56%). In the Child-Pugh class C livers, CYP1A2 was found to be reduced (to 52%). A significant trend in down-regulation of the protein abundance was documented for CYP1A2, CYP2C9, CYP3A4, CYP2E1, UGT2B7, and UGT2B15. The results of the study demonstrate that DMEs protein abundances in the liver are affected by hepatitis C virus infection and depend on the severity of the disease.

Comparison of P-glycoprotein function in peripheral blood mononuclear cells ex vivo in stable Black and White male and female kidney transplant recipients

Kidney allograft survival remains poorer in Black compared to White recipients due to racial differences in calcineurin inhibitor (CNI) pharmacology. P-glycoprotein (P-gp), an ABC efflux transporter expressed in peripheral blood mononuclear cells (PBMCs), modulates CNI pharmacokinetics and intracellular pharmacology. This study investigated P-gp function in PBMC ex vivo at 0 (trough), 4, 8, and 12 h in stable Black and White male and female kidney transplant recipients (n = 67) receiving tacrolimus and mycophenolic acid. Tacrolimus doses were adjusted to troughs of 4-10 ng/ml. P-gp function was quantified with flow cytometric measurement of cyclosporine (CYA; 2.5 μM)-reversible efflux of P-gp substrate, 3,3'-Diethyloxacarbocyanine iodide by determining the percentage change of mean fluorescent intensity (MFI) with CYA (% ΔMFI). The composite parameter of area under the concentration versus time (AUC)0-12h % ΔMFI estimated P-gp function. Data analysis examined race, sex, and race-sex associations to P-gp function. A secondary aim analyzed ABCB1 genotypes: 1236C>T (rs1128503), 2677G>T/A (rs2032582), 3435C>T (rs1045642), and P-gp function. P-gp function (% ΔMFI) was higher in White patients at troughs (p = 0.031) compared to Black counterparts with similar trends at 4 and 8 h. Reduced AUC0-12h % ΔMFI was noted in Black recipients (N = 32) compared with Whites (N = 35, p = 0.029) with notable pairwise adjusted differences between Black and White women (p = 0.021). Higher AUC0-12h % ΔMFI was associated with ABCB1 2677 TT compared to GG variants (p = 0.035). The AUC0-12h % ΔMFI was greater in White than Black subjects. P-gp function was higher at troughs in White subjects and differed between race-sex groups. P-gp function in PBMC may influence intracellular tacrolimus exposure and inter-relating pharmacodynamic responses which may support race and sex pharmacologic differences.

Proteomic quantification of receptor tyrosine kinases involved in the development and progression of colorectal cancer liver metastasis

Introduction

Alterations in expression and activity of human receptor tyrosine kinases (RTKs) are associated with cancer progression and in response to therapeutic intervention.

Methods

Thus, protein abundance of 21 RTKs was assessed in 15 healthy and 18 cancerous liver samples [2 primary and 16 colorectal cancer liver metastasis (CRLM)] matched with non-tumorous (histologically normal) tissue, by a validated QconCAT-based targeted proteomic approach.

Results

It was demonstrated, for the first time, that the abundance of EGFR, INSR, VGFR3 and AXL, is lower in tumours relative to livers from healthy individuals whilst the opposite is true for IGF1R. EPHA2 was upregulated in tumour compared with histologically normal tissue surrounding it. PGFRB levels were higher in tumours relative to both histologically normal tissue surrounding tumour and tissues taken from healthy individuals. The abundances of VGFR1/2, PGFRA, KIT, CSF1R, FLT3, FGFR1/3, ERBB2, NTRK2, TIE2, RET, and MET were, however, comparable in all samples. Statistically significant, but moderate correlations were observed (Rs > 0.50, p < 0.05) for EGFR with INSR and KIT. FGFR2 correlated with PGFRA and VGFR1 with NTRK2 in healthy livers. In non-tumorous (histologically normal) tissues from cancer patients, there were correlations between TIE2 and FGFR1, EPHA2 and VGFR3, FGFR3 and PGFRA (p < 0.05). EGFR correlated with INSR, ERBB2, KIT and EGFR, and KIT with AXL and FGFR2. In tumours, CSF1R correlated with AXL, EPHA2 with PGFRA, and NTRK2 with PGFRB and AXL. Sex, liver lobe and body mass index of donors had no impact on the abundance of RTKs, although donor age showed some correlations. RET was the most abundant of these kinases in non-tumorous tissues (~35%), while PGFRB was the most abundant RTK in tumours (~47%). Several correlations were also observed between the abundance of RTKs and proteins relevant to drug pharmacokinetics (enzymes and transporters).

Discussion

DiscussionThis study quantified perturbation to the abundance of several RTKs in cancer and the value generated in this study can be used as input to systems biology models defining liver cancer metastases and biomarkers of its progression.

Association between CYP2E1 C-1054T and 96-bp I/D genetic variations and the risk of polycystic ovary syndrome in Chinese women

PURPOSE

To investigate the association of cytochrome P450 2E1 (CYP2E1) C-1054T (rs2031920) and 96-bp I/D genetic variations with the risk of polycystic ovary syndrome (PCOS), and to estimate the effects of genotypes on the clinical, metabolic, hormonal, and oxidative stress indicators.

METHODS

This case-control study included 762 control women and 1034 patients with PCOS. Genotypes were determined using polymerase chain reaction and/or restriction fragment length polymorphism analysis. Clinical and biochemical parameters were also analyzed.

RESULTS

Frequencies of the TT + CT genotype (35.4 vs. 28.9%) and T allele (19.6 vs. 16.0%) of the CYP2E1 C-1054T polymorphism were significantly higher in the PCOS group than in the control group (OR = 1.350, 95% CI 1.103-1.652, P = 0.004 for the dominant model). Genotype TT + CT remained a significant predictor of PCOS in a logistic regression model including age, body mass index (BMI), and recruitment year of participants (OR = 1.345, 95% CI 1.071-1.688, P = 0.011). No statistical differences were found in the genotype and allele frequencies of CYP2E1 96-bp I/D polymorphism. However, the combined genotype DD/TT + CT was related to an increased risk of PCOS when the DD/CC wild-type combined genotype was used as a reference. Patients with the I allele of 96-bp I/D polymorphism had a lower BMI but higher plasma apolipoprotein B and oxidized low-density lipoprotein cholesterol levels than those with the DD genotype.

CONCLUSION

CYP2E1 C-1054T, but not 96-bp I/D, genetic polymorphism is associated with an increased risk of PCOS in Chinese women.

A simple LC-MS/MS method for the simultaneous quantification of drug metabolic enzymes

OBJECTIVE

The aim of this study is to develop a LC-MS/MS method for the quantitation of seven cytochrome P450 (CYP450) enzymes.

METHODS

A high-performance liquid chromatography-tandem mass spectrometry method was developed using multiple reaction monitoring mode with positive electrospray ionization. The method was validated with selectivity, linearity, stability, accuracy and precious. In addition, the abundance of seven CYP450 enzymes in human liver microsomes and CYP3A4 in placenta were determined using the current method.

RESULTS

The linear range for CYP1A2, CYP2B6 and CYP2C8 was 0.036-3.6 nM and for CYP2C9, CYP2C19, CYP2D6 and CYP3A4 was 0.090-9.0 nM. No interference was found between the blank matrix and each specific peptides. The accuracy and precious results were in accord with the requirement of analytical methods for biological samples in Chinese Pharmacopoeia. In addition, the peptides were stable under current stability conditions. The content of CYP3A4 in placenta and the seven CYP450 enzymes in human liver microsomes were accurately quantified.

CONCLUSION

The developed method is sensitive and specific and can be applied to the quantification of enzymes abundance in different human derived samples like placenta and liver microsomes.

Human Hepatic Transporter Signature Peptides for Quantitative Targeted Absolute Proteomics: Selection, Digestion Efficiency, and Peptide Stability

PURPOSE

Quantitative targeted absolute proteomics (QTAP) quantifies proteins by measuring the signature peptides produced from target proteins by trypsin digestion. The selection of signature peptides is critical for reliable peptide quantification. The purpose of this study was to comprehensively assess the digestion efficiency and stability of tryptic peptides and to identify optimal signature peptides for human hepatic transporters and membrane marker proteins.

METHODS

The plasma membrane fraction of the human liver was digested at different time points and the peptides were comprehensively quantified using quantitative proteomics. Transporters and membrane markers were quantified using the signature peptides by QTAP.

RESULTS

Tryptic peptides were classified into clusters with low digestion efficiency, low stability, and high digestion efficiency and stability. Using the cluster information, we found that a proline residue next to the digestion site or the peptide position in or close to the transmembrane domains lowers digestion efficiency. A peptide containing cysteine at the N-terminus or arginine-glycine lowers peptide stability. Based on this information and the time course of peptide quantification, optimal signature peptides were identified for human hepatic transporters and membrane markers. The quantification of transporters with multiple signature peptides yielded consistent absolute values with less than 30% of coefficient variants in human liver microsomes and homogenates.

CONCLUSIONS

The signature peptides selected in the present study enabled the reliable quantification of human hepatic transporters. The QTAP protocol using these optimal signature peptides provides quantitative data on hepatic transporters usable for integrated pharmacokinetic studies.

Altered Bioavailability and Pharmacokinetics in Crohn's Disease: Capturing Systems Parameters for PBPK to Assist with Predicting the Fate of Orally Administered Drugs

Backgrond and Objective

Crohn’s disease (CD) is a chronic inflammatory bowel disease that affects a wide age range. Hence, CD patients receive a variety of drugs over their life beyond those used for CD itself. The changes to the integrity of the intestine and its drug metabolising enzymes and transporters (DMETs) can alter the oral bioavailability of drugs. However, there are other changes in systems parameters determining the fate of drugs in CD, and understanding these is essential for dose adjustment in patients with CD.

Methods

The current analysis gathered all the available clinical data on the kinetics of drugs in CD (by March 2021), focusing on orally administered small molecule drugs. A meta-analysis of the systems parameters affecting oral drug pharmacokinetics was conducted. The systems information gathered on intestine, liver and blood proteins and other physiological parameters was incorporated into a physiologically based pharmacokinetic (PBPK) platform to create a virtual population of CD patients, with a view for guiding dose adjustment in the absence of clinical data in CD.

Results

There were no uniform trends in the reported changes in reported oral bioavailability. The nature of the drug as well as the formulation affected the direction and magnitude of variation in kinetics in CD patients relative to healthy volunteers. Even for the same drug, the reported changes in exposure varied, possibly due to a lack of distinction between the activity states of CD. The highest alteration was seen with S-verapamil and midazolam, 8.7- and 5.3-fold greater exposure, respectively, in active CD patients relative to healthy volunteers. Only one report was available on liver DMETs in CD, and indicated reduced CYP3A4 activity. In a number of reports, mRNA expression of DMETs in the ileum and colon of CD patients was measured, focussing on P-glycoprotein (p-gp) transporter and CYP3A4 enzyme, and showed contradictory results. No data were available on protein expression in duodenum and jejunum despite their dominant role in oral drug absorption.

Conclusion

There are currently inadequate dedicated clinical or quantitative proteomic studies in CD to enable predictive PBPK models with high confidence and adequate verification. The PBPK models for CD with the available systems parameters were able to capture the major physiological influencers and the gaps to be filled by future research. Quantification of DMETs in the intestine and the liver in CD is warranted, alongside well-defined clinical drug disposition studies with a number of index drugs as biomarkers of changes in DMETs in these patients, to avoid large-scale dedicated studies for every drug to determine the effects of disease on the drug’s metabolism and disposition and the consequential safety and therapeutic concerns.

Supplementary Information

The online version contains supplementary material available at 10.1007/s40262-022-01169-4.

Cabozantinib Carries the Risk of Drug-Drug Interactions via Inhibition of UDPglucuronosyltransferase (UGT) 1A9

BACKGROUND

Cabozantinib is a multiple receptor tyrosine kinases inhibitor (TKI) approved to treat progressive, metastatic medullary thyroid cancer, advanced renal cell carcinoma, and hepatocellular carcinoma. Drugdrug interactions (DDIs) for cabozantinib have been identified involving the role of cytochromes P450. Although the previous study reported that cabozantinib showed a slight inhibition of UDP-glucuronosyltransferase (UGT) 1A1 at the highest concentration tested, there are no reports on the potential for UGTs-mediated-DDIs. Hence, the current study aims to address this knowledge gap.

OBJECTIVE

This study aimed to investigate the inhibitory effect of cabozantinib on human UGTs and to quantitatively evaluate the DDI potential via UGT inhibition.

METHODS

The inhibitory effects of cabozantinib on UGTs were determined by measuring the formation rates for 4- methylumbelliferone (4-MU) glucuronide and trifluoperazine N-glucuronide using recombinant human UGT isoforms in the absence or presence of cabozantinib. Inhibition kinetic studies were conducted to determine the type of inhibition of cabozantinib on UGTs and the corresponding inhibition constant (Ki) value. In vitro-in vivo extrapolation (IVIVE) was further employed to predict the potential risk of DDI in vivo.

RESULTS

Cabozantinib displayed potent inhibition of UGT1A1, 1A3, 1A4, 1A6, 1A7, 1A8, 1A9, 1A10, 2B7, and 2B15. Cabozantinib exhibited noncompetitive inhibition towards UGT1A1 and 1A3 and inhibition towards UGT1A7 and 1A9. The Ki,u values (mean ± standard deviation) were calculated to be 2.15±0.11 μM, 0.83±0.05 μM, 0.75±0.04 μM and 0.18 ± 0.10 μM for UGT1A1, 1A3, 1A7 and 1A9, respectively. Co-administration of cabozantinib at the clinically approved dose of 60 mg/day or 140 mg/day may result in approximately a 26% to 60% increase in the systemic exposure of drugs predominantly cleared by UGT1A9, implying a high risk of DDIs.

CONCLUSION

Cabozantinib has the potential to cause DDIs via the inhibition of UGT1A9; therefore, additional attention should be paid to the safety of the combined use of cabozantinib and drugs metabolized by UGT1A9.

LuxS in Lactobacillus plantarum SS-128 Improves the Texture of Refrigerated Litopenaeus vannamei: Mechanism Exploration Using a Proteomics Approach

This study illustrated the texture changes of Shewanella baltica-inoculated Litopenaeus vannamei during refrigerated storage with the exogenous addition of Lactobacillus plantarum SS-128. The group inoculated with SS-128 had an improved texture compared with that inoculated with the luxS-mutant group (ΔluxS). Proteomics were conducted to analyze the protein alterations in L. vannamei and supernatant, respectively. During storage, many texture-related proteins, including myosin heavy chain and beta-actin, were maintained due to luxS. Some endogenous enzymes related to the energy metabolism and hydrolysis of L. vannamei were downregulated. The luxS-induced interaction with S. baltica showed significant changes in the expression of some critical enzymes and pathways. The ATP-dependent zinc metalloprotease FtsH and protease subunit HslV were downregulated, and the oxidative phosphorylation and glycosaminoglycan degradation pathways in S. baltica were inhibited, resulting in the slow deterioration of L. vannamei. By exploring the mechanism underlying SS-128-led manipulation of the metabolism of spoilage bacteria, we clarified the texture maintenance mechanism of luxS in SS-128, providing theoretical evidence for SS-128 application in food preservation.

Recent advances in the translation of drug metabolism and pharmacokinetics science for drug discovery and development

Drug metabolism and pharmacokinetics (DMPK) is an important branch of pharmaceutical sciences. The nature of ADME (absorption, distribution, metabolism, excretion) and PK (pharmacokinetics) inquiries during drug discovery and development has evolved in recent years from being largely descriptive to seeking a more quantitative and mechanistic understanding of the fate of drug candidates in biological systems. Tremendous progress has been made in the past decade, not only in the characterization of physiochemical properties of drugs that influence their ADME, target organ exposure, and toxicity, but also in the identification of design principles that can minimize drug-drug interaction (DDI) potentials and reduce the attritions. The importance of membrane transporters in drug disposition, efficacy, and safety, as well as the interplay with metabolic processes, has been increasingly recognized. Dramatic increases in investments on new modalities beyond traditional small and large molecule drugs, such as peptides, oligonucleotides, and antibody-drug conjugates, necessitated further innovations in bioanalytical and experimental tools for the characterization of their ADME properties. In this review, we highlight some of the most notable advances in the last decade, and provide future perspectives on potential major breakthroughs and innovations in the translation of DMPK science in various stages of drug discovery and development.

Label-Free but Still Constrained: Assessment of Global Proteomic Strategies for the Quantification of Hepatic Enzymes and Transporters

Building and refining pharmacology models require "system" data derived from tissues and in vitro systems analyzed by quantitative proteomics. Label-free global proteomics offers a wide scope of analysis, allowing simultaneous quantification of thousands of proteins per sample. The data generated from such analysis offer comprehensive protein expression profiles that can address existing gaps in models. In this study, we assessed the performance of three widely used label-free proteomic methods, "high N" ion intensity approach (HiN), intensity-based absolute quantification (iBAQ) and total protein approach (TPA), in relation to the quantification of enzymes and transporters in 27 human liver microsomal samples. Global correlations between the three methods were highly significant (R2 > 0.70, P < 0.001, n = 2232 proteins). Absolute abundances of 57 pharmacokinetic targets measured by standard-based label-free methods (HiN and iBAQ) showed good agreement, whereas the TPA overestimated abundances by two- to threefold. Relative abundance distribution of enzymes was similar for the three methods, while differences were observed with TPA in the case of transporters. Variability (CV) was similar across methods, with consistent between-sample relative quantification. The back-calculated amount of protein in the samples based on each method was compared with the nominal protein amount analyzed in the proteomic workflow, revealing overall agreement with data from the HiN method with bovine serum albumin as standard. The findings herein present a critique of label-free proteomic data relevant to pharmacokinetics and evaluate the possibility of retrospective analysis of historic datasets. SIGNIFICANCE STATEMENT: This study provides useful insights for using label-free methods to generate abundance data applicable for populating pharmacokinetic models. The data demonstrated overall correlation between intensity-based label-free proteomic methods (HiN, iBAQ and TPA), whereas iBAQ and TPA overestimated the total amount of protein in the samples. The extent of overestimation can provide a means of normalization to support absolute quantification. Importantly, between-sample relative quantification was consistent (similar variability) across methods.

Reevaluate In Vitro CYP3A Index Reactions of Benzodiazepines and Steroids between Humans and Dogs

Cytochrome P450 3A (CYP3A), the most important class of drug-metabolizing enzymes, participates in the metabolism of half of clinically used drugs. The CYP3A index reactions of dogs, one of the most widely used preclinical nonrodent species, are still poorly understood. This work evaluated the activity and selectivity of 10 CYP3A index reactions, including midazolam (MDZ) 1'- and 4-hydroxylation, alprazolam (APZ) and triazolam (TRZ) α- and 4-hydroxylation, testosterone (T) 6β-hydroxylation, lithocholate (LCA) 6α-hydroxylation, deoxycholate (DCA) 1β- and 5β-hydroxylation, with quantitative reaction phenotyping and kinetic analysis in human and canine recombinant CYP enzymes (rCYPs). In human studies, all reactions are reconfirmed as mixed index reactions of CYP3A with minor contributions from non-CYP3A isoforms. In canine studies, all reactions are also primarily catalyzed by CYP3A12 with lower contributions from CYP3A26. However, the canine CYP2B11 appreciably contributes to the hydroxylation of benzodiazepines except for APZ 4-hydroxylation. The canine CYP3A isoforms have lower activity than human isoforms toward T 6β-hydroxylation and LCA 6α-hydroxylation and both substrates undergo non-CYP3A catalyzed side reactions. DCA 1β- and 5β-hydroxylation are validated as the CYP3A index reactions in both humans and dogs with limited non-CYP3A contributions and side reactions. In conclusion, this work provides a comprehensive overview for the selectivity and activity of in vitro CYP3A index reactions in humans and dogs. The validated CYP3A index reactions between humans and dogs may benefit future practices in drug metabolism and drug interaction studies. SIGNIFICANCE STATEMENT: Dogs are one of the most important nonrodent animals with limited studies of cytochrome P450 enzymes than humans. This work provides the most comprehensive quantitative data to date for the selectivity and activity of CYP3A index reactions in humans and dogs. The canine CYP2B11 was found to appreciably contribute to hydroxylation of midazolam, alprazolam and triazolam, the well-known probes for human CYP3A. Deoxycholate 1β- and 5β-hydroxylation are validated as the CYP3A index reactions in both humans and dogs.

The Development of a PBPK Model for Atomoxetine Using Levels in Plasma, Saliva and Brain Extracellular Fluid in Patients with Normal and Deteriorated Kidney Function

BACKGROUND

Atomoxetine is a treatment for attention-deficit hyperactivity disorder. It inhibits Norepinephrine Transporters (NET) in the brain. Renal impairment can reduce hepatic CYP2D6 activity and atomoxetine elimination which may increase its body exposure. Atomoxetine can be secreted in saliva.

OBJECTIVE

The objective of this work was to test the hypothesis that atomoxetine saliva levels (sATX) can be used to predict ATX brain Extracellular Fluid (bECF) levels and their pharmacological effects in healthy subjects and those with End-Stage Renal Disease (ESRD).

METHODS

The pharmacokinetics of atomoxetine after intravenous administration to rats with chemically induced acute and chronic renal impairments were investigated. A physiologically-based pharmacokinetic (PBPK) model was built and verified in rats using previously published measured atomoxetine levels in plasma and brain tissue. The rat PBPK model was then scaled to humans and verified using published measured atomoxetine levels in plasma, saliva, and bECF.

RESULTS

The rat PBPK model predicted the observed reduced atomoxetine clearance due to renal impairment in rats. The PBPK model predicted atomoxetine exposure in human plasma, sATX and bECF. Additionally, it predicted that ATX bECF levels needed to inhibit NET are achieved at 80 mg dose. In ESRD patients, the developed PBPK model predicted that the previously reported 65% increase in plasma exposure in these patients can be associated with a 63% increase in bECF. The PBPK simulations showed that there is a significant correlation between sATX and bECF in human.

CONCLUSION

Saliva levels can be used to predict atomoxetine pharmacological response.

CYP2E1 plays a suppressive role in hepatocellular carcinoma by regulating Wnt/Dvl2/β-catenin signaling

Objective

Knowledge of the role of CYP2E1 in hepatocarcinogenesis is largely based on epidemiological and animal studies, with a primary focus on the role of CYP2E1 in metabolic activation of procarcinogens. Few studies have directly assessed the effects of CYP2E1 on HCC malignant phenotypes.

Methods

The expression of CYP2E1 in HCC tissues was determined by qRT-PCR, western blotting and immunohistochemistry. Overexpression of CYP2E1 in HCC cell was achieved by lentivirus transfection. The function of CYP2E1 were detected by CCK-8, wound healing, transwell assays, xenograft models and pulmonary metastasis model. TOP/FOPFlash reporter assay, western blotting, functional rescue experiments, Co-immunoprecipitation and reactive oxygen species detection were conducted to reveal the underlying mechanism of the tumor suppressive role of CYP2E1.

Results

CYP2E1 expression is down-regulated in HCC tissues, and this downregulation was associated with large tumor diameter, vascular invasion, poor differentiation, and shortened patient survival time. Ectopic expression of CYP2E1 inhibits the proliferation, invasion and migration and epithelial-to-mesenchymal transition of HCC cells in vitro, and inhibits tumor formation and lung metastasis in nude mice. Mechanistic investigations show that CYP2E1 overexpression significantly inhibited Wnt/β-catenin signaling activity and decreased Dvl2 expression in HCC cells. An increase in Dvl2 expression restored the malignant phenotype of HCC cells. Notably, CYP2E1 promoted the ubiquitin-mediated degradation of Dvl2 by strengthening the interaction between Dvl2 and the E3 ubiquitin ligase KLHL12 in CYP2E1-stable HCC cells. CYP2E1-induced ROS accumulation was a critical upstream event in the Wnt/β-Catenin pathway in CYP2E1-overexpressing HCC cells.

Conclusions

These results provide novel insight into the role of CYP2E1 in HCC and the tumor suppressor role of CYP2E1 can be attributed to its ability to manipulate Wnt/Dvl2/β-catenin pathway via inducing ROS accumulation, which provides a potential target for the prevention and treatment of HCC.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s12967-022-03396-6.

Quantitative Proteomics of Hepatic Drug-Metabolizing Enzymes and Transporters in Patients with Colorectal Cancer Metastasis

The impact of liver cancer metastasis on protein abundance of 22 drug‐metabolizing enzymes (DMEs) and 25 transporters was investigated using liquid chromatography‐tandem accurate mass spectrometry targeted proteomics. Microsomes were prepared from liver tissue taken from 15 healthy individuals and 18 patients with cancer (2 primary and 16 metastatic). Patient samples included tumors and matching histologically normal tissue. The levels of cytochrome P450 (CYPs 2B6, 2D6, 2E1, 3A4, and 3A5) and uridine 5′‐diphospho‐glucuronosyltransferases (UGTs 1A1, 1A6, 1A9, 2B15, 2B4, and 2B7) were lower in histologically normal tissue from patients relative to healthy controls (up to 6.6‐fold) and decreased further in tumors (up to 21‐fold for CYPs and 58‐fold for UGTs). BSEP and MRPs were also suppressed in histologically normal (up to 3.1‐fold) and tumorous tissue (up to 6.3‐fold) relative to healthy individuals. Abundance of OCT3, OAT2, OAT7, and OATPs followed similar trends (up to 2.9‐fold lower in histologically normal tissue and up to 16‐fold lower in tumors). Abundance of NTCP and OCT1 was also lower (up to 9‐fold). Interestingly, monocarboxylate transporter MCT1 was more abundant (3.3‐fold) in tumors, the only protein target to show this pattern. These perturbations could be attributed to inflammation. Interindividual variability was substantially higher in patients with cancer. Proteomics‐informed physiologically‐based pharmacokinetic (PBPK) models of 50 drugs with different attributes and hepatic extraction ratios (Simcyp) showed substantially lower drug clearance with cancer‐specific parameters compared with default parameters. In conclusion, this study provides values for decreased abundance of DMEs and transporters in liver cancer, which enables using population‐specific abundance for these patients in PBPK modeling.

Drug disposition protein quantification in matched human jejunum and liver from donors with obesity

Mathematical models, such as physiologically‐based pharmacokinetic models, are used to predict, for example, drug disposition and toxicity. However, populations differ in the abundance of proteins involved in these processes. To improve the building and refinement of such models, they must take into account these interindividual variabilities. In this study, we used global proteomics to characterize the protein composition of jejunum and liver from 37 donors with obesity enrolled in the COCKTAIL study. Liver protein levels from the 37 donors were further compared with those from donors without obesity. We quantified thousands of proteins and could present the expression of several drug‐metabolizing enzymes, for the first time, in jejunum, many of which belong to the cytochrome P450 (CYP) (e.g., CYP2U1) and the amine oxidase (flavin‐containing) (e.g., monoamine oxidase A (MAOA)) families. Although we show that many metabolizing enzymes had greater expression in liver, others had higher expression in jejunum (such as, MAOA and CES2), indicating the role of the small intestine in extrahepatic drug metabolism. We further show that proteins involved in drug disposition are not correlated in the two donor‐matched tissues. These proteins also do not correlate with physiological factors such as body mass index, age, and inflammation status in either tissue. Furthermore, the majority of these proteins are not differently expressed in donors with or without obesity. Nonetheless, interindividual differences were considerable, with implications for personalized prediction models and systems pharmacology.

A potential implication of UDP-glucuronosyltransferase 2B10 in the detoxification of drugs used in pediatric hematopoietic stem cell transplantation setting: an in silico investigation

Background

Sinusoidal occlusion syndrome (SOS) is a potentially severe complication following hematopoietic stem cell transplantation (HSCT) in pediatric patients. Treatment related risk factors such as intensity of conditioning, hepatotoxic co-medication and patient related factors such as genetic variants predispose individuals to develop SOS. The variant allele for SNP rs17146905 in UDP-glucuronosyl transferase 2B10 (UGT2B10) gene was correlated with the occurrence of SOS in an exome-wide association study. UGT2B10 is a phase II drug metabolizing enzyme involved in the N-glucuronidation of tertiary amine containing drugs.

Methods

To shed light on the functionality of UGT2B10 enzyme in the metabolism of drugs used in pediatric HSCT setting, we performed in silico screening against custom based library of putative ligands. First, a list of potential substrates for in silico analysis was prepared using a systematic consensus-based strategy. The list comprised of drugs and their metabolites used in pediatric HSCT setting. The three-dimensional structure of UGT2B10 was not available from the Research Collaboratory Structural Bioinformatics - Protein Data Bank (RCSB - PDB) repository and thus we predicted the first human UGT2B10 3D model by using multiple template homology modeling with MODELLER Version 9.2 and molecular docking calculations with AutoDock Vina Version 1.2 were implemented to quantify the estimated binding affinity between selected putative substrates or ligands and UGT2B10. Finally, we performed molecular dynamics simulations using GROMACS Version 5.1.4 to confirm the potential UGT2B10 ligands prioritized after molecular docking (exhibiting negative free binding energy).

Results

Four potential ligands for UGT2B10 namely acetaminophen, lorazepam, mycophenolic acid and voriconazole n-oxide intermediate were identified. Other metabolites of voriconazole satisfied the criteria of being possible ligands of UGT2B10. Except for bilirubin and 4-Hydroxy Voriconazole, all the ligands (particularly voriconazole and hydroxy voriconazole) are oriented in substrate binding site close to the co-factor UDP (mean ± SD; 0.72 ± 0.33 nm). Further in vitro screening of the putative ligands prioritized by in silico pipeline is warranted to understand the nature of the ligands either as inhibitors or substrates of UGT2B10.

Conclusions

These results may indicate the clinical and pharmacological relevance UGT2B10 in pediatric HSCT setting. With this systematic computational methodology, we provide a rational-, time-, and cost-effective way to identify and prioritize the interesting putative substrates or inhibitors of UGT2B10 for further testing in in vitro experiments.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12860-021-00402-5.

Ontogeny of Small Intestinal Drug Transporters and Metabolizing Enzymes Based on Targeted Quantitative Proteomics

Most drugs are administered to children orally. An information gap remains on the protein abundance of small intestinal drug-metabolizing enzymes (DMEs) and drug transporters (DTs) across the pediatric age range, which hinders precision dosing in children. To explore age-related differences in DMEs and DTs, surgical leftover intestinal tissues from pediatric and adult jejunum and ileum were collected and analyzed by targeted quantitative proteomics for apical sodium-bile acid transporter, breast cancer resistance protein (BCRP), monocarboxylate transporter 1 (MCT1), multidrug resistance protein 1 (MDR1), multidrug resistance-associated protein (MRP) 2, MRP3, organic anion-transporting polypeptide 2B1, organic cation transporter 1, peptide transporter 1 (PEPT1), CYP2C19, CYP3A4, CYP3A5, UDP glucuronosyltransferase (UGT) 1A1, UGT1A10, and UGT2B7. Samples from 58 children (48 ileums, 10 jejunums, age range: 8 weeks to 17 years) and 16 adults (8 ileums, 8 jejunums) were analyzed. When comparing age groups, BCRP, MDR1, PEPT1, and UGT1A1 abundance was significantly higher in adult ileum as compared with the pediatric ileum. Jejunal BCRP, MRP2, UGT1A1, and CYP3A4 abundance was higher in the adults compared with children 0-2 years of age. Examining the data on a continuous age scale showed that PEPT1 and UGT1A1 abundance was significantly higher, whereas MCT1 and UGT2B7 abundance was lower in adult ileum as compared with the pediatric ileum. Our data contribute to the deeper understanding of the ontogeny of small intestinal drug-metabolizing enzymes and drug transporters and shows DME-, DT-, and intestinal location-specific, age-related changes. SIGNIFICANCE STATEMENT: This is the first study that describes the ontogeny of small intestinal DTs and DMEs in human using liquid chromatography with tandem mass spectrometry-based targeted quantitative proteomics. The current analysis provides a detailed picture about the maturation of DT and DME abundances in the human jejunum and ileum. The presented results supply age-related DT and DME abundance data for building more accurate PBPK models that serve to support safer and more efficient drug dosing regimens for the pediatric population.

UGT1A1 and UGT1A3 activity and inhibition in human liver and intestinal microsomes and a recombinant UGT system under similar assay conditions using selective substrates and inhibitors

In vitro enzyme kinetics and inhibition data was compared for UGT1A1 and UGT1A3 isoforms under similar assay conditions using human liver microsomes (HLM), human intestinal microsomes (HIM) and recombinant UGT (rUGT) enzyme systems.UGT1A1 catalysed β-estradiol 3-β-D-glucuronide formation showed allosteric sigmoidal kinetics in all enzyme systems; while UGT1A3 catalysed CDCA 24-acyl-β-D-glucuronide formation exhibited Michaelis-Menten kinetics in HLM, substrate inhibition kinetics in HIM and rUGT systems. Corresponding K_m or S₅₀ concentrations of β-estradiol and CDCA were employed in the respective UGT inhibition studies.Atazanavir inhibited the production of β-estradiol 3-β-D-glucuronide with IC50 values of 0.54 µM and 0.16 µM in HLM and rUGT1A1, respectively. But its inhibition potential was not observed in HIM, indicating potential cross-talk with other high-affinity intestinal UGT isozymes. On the other hand, zafirlukast, a pan UGT inhibitor, exhibited moderate inhibition in HIM with an IC50 value of 16.70 µM. Lithocholic acid, inhibited the production of CDCA 24-acyl-β-D-glucuronide with IC50 values of 1.68, 1.84, and 12.42 µM in HLM, rUGT1A3, and HIM, respectively.These results indicated that HLM, HIM, and rUGTs may be used as complementary in vitro systems to evaluate hepatic and intestinal UGT mediated DDIs at the screening stage.

Proteomics of colorectal cancer liver metastasis: A quantitative focus on drug elimination and pharmacodynamics effects

Aims

This study aims to quantify drug‐metabolising enzymes, transporters, receptor tyrosine kinases (RTKs) and protein markers (involved in pathways affected in cancer) in pooled healthy, histologically normal and matched cancerous liver microsomes from colorectal cancer liver metastasis (CRLM) patients.

Methods

Microsomal fractionation was performed and pooled microsomes were prepared. Global and accurate mass and retention time liquid chromatography–mass spectrometry proteomics were used to quantify proteins. A QconCAT (KinCAT) for the quantification of RTKs was designed and applied for the first time. Physiologically based pharmacokinetic (PBPK) simulations were performed to assess the contribution of altered abundance of drug‐metabolising enzymes and transporters to changes in pharmacokinetics.

Results

Most CYPs and UGTs were downregulated in histologically normal relative to healthy samples, and were further reduced in cancer samples (up to 54‐fold). The transporters, MRP2/3, OAT2/7 and OATP2B1/1B3/1B1 were downregulated in CRLM. Application of abundance data in PBPK models for substrates with different attributes indicated substantially lower (up to 13‐fold) drug clearance when using cancer‐specific instead of default parameters in cancer population. Liver function markers were downregulated, while inflammation proteins were upregulated (by up to 76‐fold) in cancer samples. Various pharmacodynamics markers (e.g. RTKs) were altered in CRLM. Using global proteomics, we examined proteins in pathways relevant to cancer (such as metastasis and desmoplasia), including caveolins and collagen chains, and confirmed general over‐expression of such pathways.

Conclusion

This study highlights impaired drug metabolism, perturbed drug transport and altered abundance of cancer markers in CRLM, demonstrating the importance of population‐specific abundance data in PBPK models for cancer.

Non-uniformity of Changes in Drug-Metabolizing Enzymes and Transporters in Liver Cirrhosis: Implications for Drug Dosage Adjustment

Liver cirrhosis is a chronic disease that affects the liver structure, protein expression, and overall metabolic function. Abundance data for drug-metabolizing enzymes and transporters (DMET) across all stages of disease severity are scarce. Levels of these proteins are crucial for the accurate prediction of drug clearance in hepatically impaired patients using physiologically based pharmacokinetic (PBPK) models, which can be used to guide the selection of more precise dosing. This study aimed to experimentally quantify these proteins in human liver samples and assess how they can impact the predictive performance of the PBPK models. We determined the absolute abundance of 51 DMET proteins in human liver microsomes across the three degrees of cirrhosis severity (n = 32; 6 mild, 13 moderate, and 13 severe), compared to histologically normal controls (n = 14), using QconCAT-based targeted proteomics. The results revealed a significant but non-uniform reduction in the abundance of enzymes and transporters, from control, by 30–50% in mild, 40–70% in moderate, and 50–90% in severe cirrhosis groups. Cancer and/or non-alcoholic fatty liver disease-related cirrhosis showed larger deterioration in levels of CYP3A4, 2C8, 2E1, 1A6, UGT2B4/7, CES1, FMO3/5, EPHX1, MGST1/3, BSEP, and OATP2B1 than the cholestasis set. Drug-specific pathways together with non-uniform changes of abundance across the enzymes and transporters under various degrees of cirrhosis necessitate the use of PBPK models. As case examples, such models for repaglinide, dabigatran, and zidovudine were successful in recovering disease-related alterations in drug exposure. In conclusion, the current study provides the biological rationale behind the absence of a single dose adjustment formula for all drugs in cirrhosis and demonstrates the utility of proteomics-informed PBPK modeling for drug-specific dose adjustment in liver cirrhosis.

Biopharmaceutical quality control with mass spectrometry

Mass spectrometry (MS) is a powerful technique for protein identification, quantification and characterization that is widely applied in biochemical studies, and which can provide data on the quantity, structural integrity and post-translational modifications of proteins. It is therefore a versatile and widely used analytic tool for quality control of biopharmaceuticals, especially in quantifying host-cell protein impurities, identifying post-translation modifications and structural characterization of biopharmaceutical proteins. Here, we summarize recent advances in MS-based analyses of these key quality attributes of the biopharmaceutical development and manufacturing processes.

Quantitative Proteomic Map of Enzymes and Transporters in the Human Kidney: Stepping Closer to Mechanistic Kidney Models to Define Local Kinetics

The applications of translational modeling of local drug concentrations in various organs had a sharp increase over the last decade. These are part of the model-informed drug development initiative, adopted by the pharmaceutical industry and promoted by drug regulatory agencies. With respect to the kidney, the models serve as a bridge for understanding animal vs. human observations related to renal drug disposition and any consequential adverse effects. However, quantitative data on key drug-metabolizing enzymes and transporters relevant for predicting renal drug disposition are limited. Using targeted and global quantitative proteomics, we determined the abundance of multiple enzymes and transporters in 20 human kidney cortex samples. Nine enzymes and 22 transporters were quantified (8 for the first time in the kidneys). In addition, > 4,000 proteins were identified and used to form an open database. CYP2B6, CYP3A5, and CYP4F2 showed comparable, but generally low expression, whereas UGT1A9 and UGT2B7 levels were the highest. Significant correlation between abundance and activity (measured by mycophenolic acid clearance) was observed for UGT1A9 (Rs = 0.65, P = 0.004) and UGT2B7 (Rs = 0.70, P = 0.023). Expression of P-gp ≈ MATE-1 and OATP4C1 transporters were high. Strong intercorrelations were observed between several transporters (P-gp/MRP4, MRP2/OAT3, and OAT3/OAT4); no correlation in expression was apparent for functionally related transporters (OCT2/MATEs). This study extends our knowledge of pharmacologically relevant proteins in the kidney cortex, with implications on more prudent use of mechanistic kidney models under the general framework of quantitative systems pharmacology and toxicology.

Mass spectrometry-based targeted proteomics method for the quantification of clinically relevant drug metabolizing enzymes in human specimens

Biotransformation by phase I and II metabolizing enzymes represents the major determinant for the oral bioavailability of many drugs. To estimate the pharmacokinetics, data on protein abundance of hepatic and extrahepatic tissues, such as the small intestine, are required. Targeted proteomics assays are nowadays state-of-the-art for absolute protein quantification and several methods for quantification of drug metabolizing enzymes have been published. However, some enzymes remain still uncovered by the analytical spectra of those methods. Therefore, we developed and validated a quantification assay for two carboxylesterases (CES-1, CES-2), 17 cytochrome P450 enzymes (CYP) (CYP1A1, CYP1A2, CYP2A6, CYP2B6, CYP2C8, CYP2C9, CYP2C18, CYP2C19, CYP2D6, CYP2E1, CYP2J2, CYP3A4, CYP3A5, CYP3A7, CYP4F2, CYP4F12, CYP4A11) and five UDP-glucuronosyltransferases (UGTs) (UGT1A1, UGT1A3, UGT2B7, UGT2B15, UGT2B17). Protein quantification was performed by analyzing proteospecific surrogate peptides after tryptic digestion with stable isotope-labelled standards. Chromatographic separation was performed on a Kinetex® 2.6 µm C18 100 Å core-shell column (100 × 2.1 mm) with a gradient elution using 0.1% formic acid and acetonitrile containing 0.1% formic acid with a flow rate of 200 µl/min. Three mass transitions were simultaneously monitored with a scheduled multiple reaction monitoring (sMRM) method for each analyte and standard. The method was partly validated according to current bioanalytical guidelines and met the criteria regarding linearity (0.1-25 nmol/L), within-day and between-day accuracy and precision as well as multiple stability criteria. Finally, the developed method was successfully applied to determine the abundance of the aforementioned enzymes in human intestinal und liver microsomes. Our work offers a new fit for purpose method for the absolute quantification of CES, CYPs and UGTs in various human tissues and can be used for the acquisition of data for physiologically based pharmacokinetic modelling.

Proteomic Quantification of Changes in Abundance of Drug-Metabolizing Enzymes and Drug Transporters in Human Liver Cirrhosis: Different Methods, Similar Outcomes

Model-based assessment of the effects of liver disease on drug pharmacokinetics requires quantification of changes in enzymes and transporters responsible for drug metabolism and disposition. Different proteomic methods are currently used for protein quantification in tissues and in vitro systems, each with specific procedures and requirements. The outcome of quantitative proteomic assays using four different methods (one targeted and three label-free) applied to the same sample set was compared in this study. Three pooled cirrhotic liver microsomal samples corresponding to cirrhosis with nonalcoholic fatty liver disease, biliary disease, or cancer and a control microsomal pool were analyzed using quantification concatemer-based targeted proteomics, the total protein approach (TPA), high three ion intensity (Hi3) approach, and intensity-based absolute quantification (iBAQ) to determine the absolute and relative abundance in disease compared with control. The relative abundance data provided a "disease perturbation factor" (DPF) for each target protein. Absolute and relative abundances generated by standard-based label-free methods (iBAQ and Hi3) showed good agreement with targeted proteomics (limited bias and scatter), but TPA (standard-free method) overestimated absolute abundances by approximately 2-fold. The DPF was consistent between different proteomic methods but varied between enzymes and transporters, indicating discordance of effects of cirrhosis on various metabolism-related proteins. The DPF ranged from no change (e.g., for glucuronosyltransferase-1A6 in nonalcoholic fatty liver disease group) to less than 0.3 (e.g., carboxylesterases-1 in cirrhosis of biliary origin). SIGNIFICANCE STATEMENT: This study demonstrated that relative changes in enzymes and transporters (DPF) are independent of the quantitative proteomic methods used. Standard-based label-free methods, such as high three ion intensity (Hi3) and intensity-based absolute quantification (iBAQ) methods, were less biased and more precise than the total protein approach (TPA) when compared with targeted data. The DPF reconciled differences across proteomic methods observed with absolute levels. Using this approach, differences were revealed in the expression of enzymes/transporters in cirrhosis associated with different etiologies.

Comparative Proteomics Analysis of the Postmitochondrial Supernatant Fraction of Human Lens-Free Whole Eye and Liver

The increasing incidence of ocular diseases has accelerated research into therapeutic interventions needed for the eye. Ocular enzymes play important roles in the metabolism of drugs and endobiotics. Various ocular drugs are designed as prodrugs that are activated by ocular enzymes. Moreover, ocular enzymes have been implicated in the bioactivation of drugs to their toxic metabolites. The key purpose of this study was to compare global proteomes of the pooled samples of the eye (n = 11) and the liver (n = 50) with a detailed analysis of the abundance of enzymes involved in the metabolism of xenobiotics and endobiotics. We used the postmitochondrial supernatant fraction (S9 fraction) of the lens-free whole eye homogenate as a model to allow accurate comparison with the liver S9 fraction. A total of 269 proteins (including 23 metabolic enzymes) were detected exclusively in the pooled eye S9 against 648 proteins in the liver S9 (including 174 metabolic enzymes), whereas 424 proteins (including 94 metabolic enzymes) were detected in both the organs. The major hepatic cytochrome P450 and UDP-glucuronosyltransferases enzymes were not detected, but aldehyde dehydrogenases and glutathione transferases were the predominant proteins in the eye. The comparative qualitative and quantitative proteomics data in the eye versus liver is expected to help in explaining differential metabolic and physiologic activities in the eye. SIGNIFICANCE STATEMENT: Information on the enzymes involved in xenobiotic and endobiotic metabolism in the human eye in relation to the liver is scarcely available. The study employed global proteomic analysis to compare the proteomes of the lens-free whole eye and the liver with a detailed analysis of the enzymes involved in xenobiotic and endobiotic metabolism. These data will help in better understanding of the ocular metabolism and activation of drugs and endobiotics.

Impact of carboxylesterase 1 genetic polymorphism on trandolapril activation in human liver and the pharmacokinetics and pharmacodynamics in healthy volunteers

Trandolapril, an angiotensin-converting enzyme inhibitor prodrug, needs to be activated by carboxylesterase 1 (CES1) in the liver to exert its intended therapeutic effect. A previous in vitro study demonstrated that the CES1 genetic variant G143E (rs71647871) abolished CES1-mediated trandolapril activation in cells transfected with the variant. This study aimed to determine the effect of the G143E variant on trandolapril activation in human livers and the pharmacokinetics (PKs) and pharmacodynamics (PDs) in human subjects. We performed an in vitro incubation study to assess trandolapril activation in human livers (5 G143E heterozygotes and 97 noncarriers) and conducted a single-dose (1 mg) PK and PD study of trandolapril in healthy volunteers (8 G143E heterozygotes and 11 noncarriers). The incubation study revealed that the mean trandolapril activation rate in G143E heterozygous livers was 42% of those not carrying the variant (p = 0.0015). The clinical study showed that, relative to noncarriers, G143E carriers exhibited 20% and 15% decreases, respectively, in the peak concentration (Cmax ) and area under the curve from 0 to 72 h (AUC0-72 h ) of the active metabolite trandolaprilat, although the differences were not statistically significant. Additionally, the average maximum reductions of systolic blood pressure and diastolic blood pressure in carriers were ~ 22% and 23% less than in noncarriers, respectively, but the differences did not reach a statistically significant level. In summary, the CES1 G143E variant markedly impaired trandolapril activation in the human liver under the in vitro incubation conditions; however, this variant had only a modest impact on the PK and PD of trandolapril in healthy human subjects.

Effect of CES1 genetic variation on enalapril steady-state pharmacokinetics and pharmacodynamics in healthy subjects

AIMS

Enalapril is a prodrug and needs to be activated by carboxylesterase 1 (CES1). A previous in vitro study demonstrated the CES1 genetic variant, G143E (rs71647871), significantly impaired enalapril activation. Two previous clinical studies examined the impact of G143E on single-dose enalapril PK (10 mg); however, the results were inconclusive. A prospective, multi-dose, pharmacokinetics and pharmacodynamics (PK/PD) study was conducted to determine the impact of the CES1 G143E variant on enalapril steady-state PK and PD in healthy volunteers.

METHODS

Study participants were stratified to G143E non-carriers (n = 15) and G143E carriers (n = 6). All the carriers were G143E heterozygotes. Study subjects received enalapril 10 mg daily for seven consecutive days prior to a 72 hour PK/PD study. Plasma concentrations of enalapril and its active metabolite enalaprilat were quantified by an established liquid chromatography-tandem mass spectrometry (LC-MS/MS) method.

RESULTS

The CES1 G143E carriers had 30.9% lower enalaprilat C_max (P =  0.03) compared to the non-carriers (38.01 vs. 55.01 ng/mL). The carrier group had 27.5% lower AUC_0-∞ (P =  0.02) of plasma enalaprilat compared to the non-carriers (374.29 vs. 515.91 ng*h/mL). The carriers also had a 32.3% lower enalaprilat-to-enalapril AUC_0-∞ ratio (P =  0.003) relative to the non-carriers. The average maximum reduction of systolic blood pressure in the non-carrier group was approximately 12.4% at the end of the study compared to the baseline (P =  0.001). No statistically significant blood pressure reduction was observed in the G143E carriers.

CONCLUSIONS

The CES1 loss-of-function G143E variant significantly impaired enalapril activation and its systolic blood pressure-lowering effect in healthy volunteers.

Physiologically-based pharmacokinetic model for alectinib, ruxolitinib, and panobinostat in the presence of cancer, renal impairment, and hepatic impairment

Renal (RIP) and hepatic (HIP) impairments are prevalent conditions in cancer patients. They can cause changes in gastric emptying time, albumin levels, hematocrit, glomerular filtration rate, hepatic functional volume, blood flow rates, and metabolic activity that can modify drug pharmacokinetics. Performing clinical studies in such populations has ethical and practical issues. Using predictive physiologically-based pharmacokinetic (PBPK) models in the evaluation of the PK of alectinib, ruxolitinib, and panobinostat exposures in the presence of cancer, RIP, and HIP can help in using optimal doses with lower toxicity in these populations. Verified PBPK models were customized under scrutiny to account for the pathophysiological changes induced in these diseases. The PBPK model-predicted plasma exposures in patients with different health conditions within average 2-fold error. The PBPK model predicted an area under the curve ratio (AUCR) of 1, and 1.8, for ruxolitinib and panobinostat, respectively, in the presence of severe RIP. On the other hand, the severe HIP was associated with AUCR of 1.4, 2.9, and 1.8 for alectinib, ruxolitinib, and panobinostat, respectively, in agreement with the observed AUCR. Moreover, the PBPK model predicted that alectinib therapeutic cerebrospinal fluid levels are achieved in patients with non-small cell lung cancer, moderate HIP, and severe HIP at 1-, 1.5-, and 1.8-fold that of healthy subjects. The customized PBPK models showed promising ethical alternatives for simulating clinical studies in patients with cancer, RIP, and HIP. More work is needed to quantify other pathophysiological changes induced by simultaneous affliction by cancer and RIP or HIP.

Influence of Proteome Profiles and Intracellular Drug Exposure on Differences in CYP Activity in Donor-Matched Human Liver Microsomes and Hepatocytes

Human liver microsomes (HLM) and human hepatocytes (HH) are important in vitro systems for studies of intrinsic drug clearance (CLint) in the liver. However, the CLint values are often in disagreement for these two systems. Here, we investigated these differences in a side-by-side comparison of drug metabolism in HLM and HH prepared from 15 matched donors. Protein expression and intracellular unbound drug concentration (Kpuu) effects on the CLint were investigated for five prototypical probe substrates (bupropion-CYP2B6, diclofenac-CYP2C9, omeprazole-CYP2C19, bufuralol-CYP2D6, and midazolam-CYP3A4). The samples were donor-matched to compensate for inter-individual variability but still showed systematic differences in CLint. Global proteomics analysis outlined differences in HLM from HH and homogenates of human liver (HL), indicating variable enrichment of ER-localized cytochrome P450 (CYP) enzymes in the HLM preparation. This suggests that the HLM may not equally and accurately capture metabolic capacity for all CYPs. Scaling CLint with CYP amounts and Kpuu could only partly explain the discordance in absolute values of CLint for the five substrates. Nevertheless, scaling with CYP amounts improved the agreement in rank order for the majority of the substrates. Other factors, such as contribution of additional enzymes and variability in the proportions of active and inactive CYP enzymes in HLM and HH, may have to be considered to avoid the use of empirical scaling factors for prediction of drug metabolism.

Label-Free Quantitative Proteomics and Substrate-Based Mass Spectrometry Imaging of Xenobiotic Metabolizing Enzymes in Ex Vivo Human Skin and a Human Living Skin Equivalent Model

We report for the first time label-free quantification of xenobiotic metabolizing enzymes (XME), transporters, redox enzymes, proteases, and nucleases in six human skin explants and a three-dimensional living skin equivalent model from LabSkin. We aimed to evaluate the suitability of LabSkin as an alternative to animal testing for the development of topical formulations. More than 2000 proteins were identified and quantified from total cellular protein. Alcohol dehydrogenase 1C, the most abundant phase I XME in human skin, and glutathione S-transferase pi 1, the most abundant phase II XME in human skin, were present in similar abundance in LabSkin. Several esterases were quantified and esterase activity was confirmed in LabSkin using substrate-based mass spectrometry imaging. No cytochrome P450 (P450) activity was observed for the substrates tested, in agreement with the proteomics data, where the cognate P450s were absent in both human skin and LabSkin. Label-free protein quantification allowed insights into other related processes such as redox homeostasis and proteolysis. For example, the most abundant antioxidant enzymes were thioredoxin and peroxiredoxin-1. This systematic determination of functional equivalence between human skin and LabSkin is a key step toward the construction of a representative human in vitro skin model, which can be used as an alternative to current animal-based tests for chemical safety and for predicting dosage of topically administered drugs. SIGNIFICANCE STATEMENT: The use of label-free quantitative mass spectrometry to elucidate the abundance of xenobiotic metabolizing enzymes, transporters, redox enzymes, proteases, and nucleases in human skin enhance our understanding of the skin physiology and biotransformation of topical drugs and cosmetics. This will help to develop mathematical models to predict drug metabolism in human skin and to develop more robust in vitro engineered human skin tissue as alternatives to animal testing.