Variability in Mass Spectrometry-based Quantification of Clinically Relevant Drug Transporters and Drug Metabolizing Enzymes

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.

Utilizing a Dual Human Transporter MDCKII-MDR1-BCRP Cell Line to Assess Efflux at the Blood Brain Barrier

To facilitate the design of drugs readily able to cross the blood brain barrier (BBB), a Madin-Darby canine kidney (MDCK) cell line was established that over expresses both P-glycoprotein (Pgp) and breast cancer resistance protein (BCRP), the main human efflux transporters of the BBB. Proteomics analyses indicate BCRP is expressed at a higher level than Pgp in this cell line. This cell line shows good activity for both transporters [BCRP substrate dantrolene efflux ratio (ER) 16.3 ± 0.9, Pgp substrate quinidine ER 27.5 ± 1.2], and use of selective transporter inhibitors enables an assessment of the relative contributions to overall ERs. The MDCKII-MDR1-BCRP ER negatively correlates with rat unbound brain/unbound plasma ratio, Kpuu Highly brain penetrant compounds with rat Kpuu ≥ 0.3 show ERs ≤ 2 in the MDCKII-MDR1-BCRP assay while compounds predominantly excluded from the brain, Kpuu ≤ 0.05, demonstrate ERs ≥ 20. A subset of compounds with MDCKII-MDR1-BCRP ER < 2 and rat Kpuu < 0.3 were shown to be substrates of rat Pgp using a rat transfected cell line, MDCKII-rMdr1a. These compounds also showed ERs > 2 in the human National Institutes of Health (NIH) MDCKI-MDR1 (high Pgp expression) cell line, which suggests that they are weak human Pgp substrates. Characterization of 37 drugs targeting the central nervous system in the MDCKII-MDR1-BCRP efflux assay show 36 have ERs < 2. In drug discovery, use of the MDCKII-MDR1-BCRP in parallel with the NIH MDCKI-MDR1 cell line is useful for identification of compounds with high brain penetration. SIGNIFICANCE STATEMENT: A single cell line that includes both the major human efflux transporters of the blood brain barrier (MDCKII-MDR1-BCRP) has been established facilitating the rapid identification of efflux substrates and enabling the design of brain penetrant molecules. Efflux ratios using this cell line demonstrate a clear relationship with brain penetration as defined by rat brain Kpuu.

Determination of single-molecule transport activity of OATP2B1 by measuring the number of transporter molecules using electrophysiological approach

Transporter-mediated clearance is determined by two factors, its single-molecule clearance, and expression level. However, no reliable method has been developed to evaluate them separately. This study aimed to develop a reliable method for evaluating the single-molecule activity of membrane transporters, such as organic anion transporting polypeptide (OATP) 2B1. HEK293 cells that co-expressed large conductance calcium-activated potassium (BK) channel and OATP2B1 were established and used for the following experiments. i) BK channel-mediated whole-cell conductance was measured using patch-clamp technique and divided by its unitary conductance to estimate the number of channels on plasma membrane (QI). ii) Using plasma membrane fraction, quantitative targeted absolute proteomics determined the stoichiometric ratio (ρ) of OATP2B1 to BK channel. iii) The uptake of estrone 3-sulfate was evaluated to calculate the Michaelis constant and uptake clearance (CL) per cell. Single-molecule clearance (CLint) was calculated by dividing CL by QI·ρ. QI and ρ values were estimated to be 916 and 2.16, respectively, yielding CLint of 5.23 fL/min/molecule. We successfully developed a novel method to reliably measure the single-molecule activity of a transporter, which could be used to evaluate the influences of factors such as genetic variations and post-translational modifications on the intrinsic activity of transporters.

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.

Epithelial and microbial determinants of colonic drug distribution

A dynamic epithelium and a rich microbiota, separated by multi-layered mucus, make up the complex colonic cellular environment. Both cellular systems are characterized by high inter- and intraindividual differences, but their impact on drug distribution and efficacy remains incompletely understood. This research gap is pressing, as, e.g., inflammatory disorders of the colon are on the rise globally. In an effort to help close this gap, we provide considerations on determining colonic epithelial and microbial cellular parameters, and their impact on drug bioavailability. First, we cover the major cell types found in vivo within the epithelium and microbiota, and discuss how they can be modeled in vitro. We then draw attention to their structural similarities and differences with regard to determinants of drug distribution. Once a drug is solubilized in the luminal fluids, there are two main classes of such determinants: 1) binding processes, and 2) transporters and drug-metabolizing enzymes. Binding lowers the unbound intracellular fraction (f_u,cell), which will, in turn, limit the amount of drug available for transport to desired sites. Transporters and drug metabolizing enzymes are ADME proteins impacting intracellular accumulation (Kp). Across cell types, we point out which processes are likely particularly impactful. Together, f_u,cell and Kp can be used to describe intracellular bioavailability (F_ic), which is a measure of local drug distribution, with consequences for efficacy. Determining these cellular parameters will be beneficial in understanding colonic drug distribution and will advance the field of drug delivery.

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

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

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.

New Approach Methods for Hazard Identification: A Case Study with Azole Fungicides Affecting Molecular Targets Associated with the Adverse Outcome Pathway for Cholestasis

Triazole fungicides such as propiconazole (Pi) or tebuconazole (Te) show hepatotoxicity in vivo, e.g., hypertrophy and vacuolization of liver cells following interaction with nuclear receptors such as PXR (pregnane-X-receptor) and CAR (constitutive androstane receptor). Accordingly, azoles affect gene expression associated with these adverse outcomes in vivo but also in human liver cells in vitro. Additionally, genes indicative of liver cholestasis are affected in vivo and in vitro. We therefore analyzed the capability of Pi and Te to cause cholestasis in an adverse outcome pathway (AOP)-driven approach in hepatic cells of human origin in vitro, considering also previous in vivo studies. Bile salt export pump (BSEP) activity assays confirmed that both azoles are weak inhibitors of BSEP. They alternate the expression of various cholestasis-associated target genes and proteins as well as the mitochondrial membrane function. Published in vivo data, however, demonstrate that neither Pi nor Te cause cholestasis in rodent bioassays. This discrepancy can be explained by the in vivo concentrations of both azoles being well below their EC50 for BSEP inhibition. From a regulatory perspective, this illustrates that toxicogenomics and human in vitro models are valuable tools to detect the potential of a substance to cause a specific type of toxicity. To come to a sound regulatory conclusion on the in vivo relevance of such a finding, results will have to be considered in a broader context also including toxicokinetics in a weight-of-evidence approach.

Innovative approaches and recent advances in the study of ontogeny of drug metabolism and transport

The disposition of a drug is driven by various processes, such as drug metabolism, drug transport, glomerular filtration and body composition. These processes are subject to developmental changes reflecting growth and maturation along the paediatric continuum. However, knowledge gaps exist on these changes and their clinical impact. Filling these gaps may aid better prediction of drug disposition and creation of age-appropriate dosing guidelines. We present innovative approaches to study these developmental changes in relation to drug metabolism and transport. First, analytical methods such as including liquid chromatography-mass spectrometry for proteomic analyses allow quantitation of the expressions of a wide variety of proteins, e.g. membrane transporters, in a small piece of organ tissue. The latter is specifically important for paediatric research, where tissues are scarcely available. Second, innovative study designs using radioactive labelled microtracers allowed study-without risk for the child-of the oral bioavailability of compounds used as markers for certain drug metabolism pathways. Third, the use of modelling and simulation to support dosing recommendations for children is supported by both the European Medicines Agency and the US Food and Drug Administration. This may even do away with the need for a paediatric trial. Physiologically based pharmacokinetics models, which include age-specific physiological information are, therefore, increasingly being used, not only to aid paediatric drug development but also to improve existing drug therapies.

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.

Simultaneous absolute protein quantification of seven cytochrome P450 isoforms in rat liver microsomes by LC-MS/MS-based isotope internal standard method

The cytochrome P450 (CYP) enzymes play a pivotal role in drug metabolism. LC-MS/MS-based targeting technology has been applied to the analysis of CYP enzymes, promoting drug development and drug-drug interaction studies. Rat is one of the most commonly used models for drug metabolism assessment, but LC-MS/MS assay quantifying the abundance of CYP enzymes in rats is rarely reported. Herein, an accurate and stable LC-MS/MS based method was developed and validated for the simultaneous quantification of seven major rat CYP isoforms (CYP1A2, 2B1, 2C6, 2C11, 2D1, 2E1, and 3A1) in liver microsomes. The careful optimization of trypsin digestion and chromatography combined with isotope-labeled peptide as internal standard improved the efficiency and accuracy of the analysis. Highly specific surrogate peptides were obtained by a procedure including trypsin digestion for six hours and separated on a Hypersil Gold C18 column (100 × 2.1 mm, 3 μm) using gradient elution for 15 min with a mobile phase of water containing 0.1% formic acid and acetonitrile. In the method validation, linearity, matrix effect, recovery, stability, accuracy, and precision all meet the requirements. Subsequently, this method was applied to detect seven enzymes in rat liver microsomes from four different sources, and the correlation between the abundance and activity of CYP enzymes was further analyzed. The high-throughput detection method provided in this study will provide support for pertinent pharmaceutical research based on rat models.

Proteomic analysis of hepatic effects of phenobarbital in mice with humanized liver

Activation of the constitutive androstane receptor (CAR) may induce adaptive but also adverse effects in rodent liver, including the induction of drug-metabolizing enzymes, transient hepatocellular proliferation, and promotion of liver tumor growth. Human relevance of CAR-related adverse hepatic effects is controversially debated. Here, we used the chimeric FRG-KO mouse model with livers largely repopulated by human hepatocytes, in order to study human hepatocytes and their response to treatment with the model CAR activator phenobarbital (PB) in vivo. Mice received an intraperitoneal injection with 50 mg/kg body weight PB or saline, and were sacrificed after 72–144 h. Non-repopulated FRG-KO mice were used as additional control. Comprehensive proteomics datasets were generated by merging data obtained by targeted as well as non-targeted proteomics approaches. For the first time, a novel proteomics workflow was established to comparatively analyze the effects of PB on human and murine proteins within one sample. Analysis of merged proteome data sets and bioinformatics data mining revealed comparable responses in murine and human hepatocytes with respect to nuclear receptor activation and induction of xenobiotic metabolism. By contrast, activation of MYC, a key regulator of proliferation, was predicted only for mouse but not human hepatocytes. Analyses of 5-bromo-2′-deoxyuridine incorporation confirmed this finding. In summary, this study for the first time presents a comprehensive proteomic analysis of CAR-dependent effects in human and mouse hepatocytes from humanized FRG-KO mice. The data support the hypothesis that PB does induce adaptive metabolic responses, but not hepatocellular proliferation in human hepatocytes in vivo.

Hepatic Expression of the Na+-Taurocholate Cotransporting Polypeptide Is Independent from Genetic Variation

The hepatic Na+-taurocholate cotransporting polypeptide NTCP/SLC10A1 is important for the uptake of bile salts and selected drugs. Its inhibition results in increased systemic bile salt concentrations. NTCP is also the entry receptor for the hepatitis B/D virus. We investigated interindividual hepatic SLC10A1/NTCP expression using various omics technologies. SLC10A1/NTCP mRNA expression/protein abundance was quantified in well-characterized 143 human livers by real-time PCR and LC-MS/MS-based targeted proteomics. Genome-wide SNP arrays and SLC10A1 next-generation sequencing were used for genomic analyses. SLC10A1 DNA methylation was assessed through MALDI-TOF MS. Transcriptomics and untargeted metabolomics (UHPLC-Q-TOF-MS) were correlated to identify NTCP-related metabolic pathways. SLC10A1 mRNA and NTCP protein levels varied 44-fold and 10.4-fold, respectively. Non-genetic factors (e.g., smoking, alcohol consumption) influenced significantly NTCP expression. Genetic variants in SLC10A1 or other genes do not explain expression variability which was validated in livers (n = 50) from The Cancer Genome Atlas. The identified two missense SLC10A1 variants did not impair transport function in transfectants. Specific CpG sites in SLC10A1 as well as single metabolic alterations and pathways (e.g., peroxisomal and bile acid synthesis) were significantly associated with expression. Inter-individual variability of NTCP expression is multifactorial with the contribution of clinical factors, DNA methylation, transcriptional regulation as well as hepatic metabolism, but not genetic variation.

Novel Insights in Drug Transporter Sciences: the Year 2021 in Review

On behalf of the team I am pleased to present the second annual 'novel insights into drug transporter sciences review' focused on peer-reviewed articles that were published in the year 2021. In compiling the articles for inclusion, preprints available in 2021 but officially published in 2022 were considered to be in scope. To support this review the contributing authors independently selected one or two articles that were thought to be impactful and of interest to the broader research community. A similar approach as published last year was adopted whereby key observations, methods and analysis of each paper is concisely summarized in the synopsis followed by a commentary highlighting the impact of the paper in understanding of drug transporters' role in drug disposition.As the goal of this review is not to provide a comprehensive overview of each paper but rather highlight important findings that are well supported by the data, the reader is encouraged to consult the original articles for additional information. Further, and keeping in line with the goals of this review, it should be noted that all authors actively contributed by writing synopsis and commentary for individual papers and no attempt was made to standardize language or writing styles. In this way, the review article is reflective of not only the diversity of the articles but also that of the contributors. I extend my thanks to the authors for their continued support and also welcome Diane Ramsden and Pallabi Mitra as contributing authors for this issue.

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.

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.

High Throughput Screening of a Prescription Drug Library for Inhibitors of Organic Cation Transporter 3, OCT3

INTRODUCTION

The organic cation transporter 3 (OCT3, SLC22A3) is ubiquitously expressed and interacts with a wide array of compounds including endogenous molecules, environmental toxins and prescription drugs. Understudied as a determinant of pharmacokinetics and pharmacodynamics, OCT3 has the potential to be a major determinant of drug absorption and disposition and to be a target for drug-drug interactions (DDIs).

GOAL

The goal of the current study was to identify prescription drug inhibitors of OCT3.

METHODS

We screened a compound library consisting of 2556 prescription drugs, bioactive molecules, and natural products using a high throughput assay in HEK-293 cells stably expressing OCT3.

RESULTS

We identified 210 compounds that at 20 μM inhibit 50% or more of OCT3-mediated uptake of 4-Di-1-ASP (2 μM). Of these, nine were predicted to inhibit the transporter at clinically relevant unbound plasma concentrations. A Structure-Activity Relationship (SAR) model included molecular descriptors that could discriminate between inhibitors and non-inhibitors of OCT3 and was used to identify additional OCT3 inhibitors. Proteomics of human brain microvessels (BMVs) indicated that OCT3 is the highest expressed OCT in the human blood-brain barrier (BBB).

CONCLUSIONS

This study represents the largest screen to identify prescription drug inhibitors of OCT3. Several are sufficiently potent to inhibit the transporter at therapeutic unbound plasma levels, potentially leading to DDIs or off-target pharmacologic effects.

Data-independent acquisition (DIA): An emerging proteomics technology for analysis of drug-metabolizing enzymes and transporters

Data-independent acquisition (DIA) proteomics is a recently-developed global mass spectrometry (MS)-based proteomics strategy. In a DIA method, precursor ions are isolated into pre-defined isolation windows and fragmented; all fragmented ions in each window are then analyzed by a high-resolution mass spectrometer. DIA proteomics analysis is characterized by a broad protein coverage, high reproducibility, and accuracy, and its combination with advances in other techniques such as sample preparation and computational data analysis could lead to further improvements in assay performances. DIA technology has been increasingly utilized in various proteomics studies, including quantifying drug-metabolizing enzymes and transporters. Quantitative proteomics study of drug-metabolizing enzymes and transporters could lead to a better understanding of pharmacokinetics and pharmacodynamics and facilitate drug development. This review summarizes the application of DIA technology in proteomic analysis of drug-metabolizing enzymes and transporters.

Characterization of hepatic zonation in mice by mass-spectrometric and antibody-based proteomics approaches

Periportal and perivenous hepatocytes show zonal heterogeneity in metabolism and signaling. Here, hepatic zonation in mouse liver was analyzed by non-targeted mass spectrometry (MS) and by the antibody-based DigiWest technique, yielding a comprehensive overview of protein expression in periportal and perivenous hepatocytes. Targeted immunoaffinity-based proteomics were used to substantiate findings related to drug metabolism. 165 (MS) and 82 (DigiWest) zonated proteins were identified based on the selected criteria for statistical significance, including 7 (MS) and 43 (DigiWest) proteins not identified as zonated before. New zonated proteins especially comprised kinases and phosphatases related to growth factor-dependent signaling, with mainly periportal localization. Moreover, the mainly perivenous zonation of a large panel of cytochrome P450 enzymes was characterized. DigiWest data were shown to complement the MS results, substantially improving possibilities to bioinformatically identify zonated biological processes. Data mining revealed key regulators and pathways preferentially active in either periportal or perivenous hepatocytes, with β-catenin signaling and nuclear xeno-sensing receptors as the most prominent perivenous regulators, and several kinase- and G-protein-dependent signaling cascades active mainly in periportal hepatocytes. In summary, the present data substantially broaden our knowledge of hepatic zonation in mouse liver at the protein level.

Damage-Net: A program for DNA repair meta-analysis identifies a network of novel repair genes that facilitate cancer evolution

The advent of genome-wide methods for identifying novel components in biological processes including CRISPR screens and proteomic studies, has transformed the research landscape within the biological sciences. However, each study normally investigates a single aspect of a process without integration of other published datasets. Here, we present Damage-Net, a program with a curated database of published results from a broad range of studies investigating DNA repair, that facilitates simple and quick meta-analysis. Users can incorporate their own datasets for analysis, and query genes of interest in the database. Importantly, this program also allows users to examine the correlation of genes of interest with pan-cancer patient survival and mutational burden effects. Interrogating these datasets revealed a network of genes that associated with cancer progression in adrenocortical carcinoma via facilitating mutational burden, ultimately contributing substantially to adrenocortical carcinoma's poor prognosis. Download at www.damage-net.co.uk.

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

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

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.

Evaluation of Normothermic Machine Perfusion of Porcine Livers as a Novel Preclinical Model to Predict Biliary Clearance and Transporter-Mediated Drug-Drug Interactions Using Statins

There is a lack of translational preclinical models that can predict hepatic handling of drugs. In this study, we aimed to evaluate the applicability of normothermic machine perfusion (NMP) of porcine livers as a novel ex vivo model to predict hepatic clearance, biliary excretion, and plasma exposure of drugs. For this evaluation, we dosed atorvastatin, pitavastatin, and rosuvastatin as model drugs to porcine livers and studied the effect of common drug-drug interactions (DDIs) on these processes. After 120 minutes of perfusion, 0.104 mg atorvastatin (n = 3), 0.140 mg pitavastatin (n = 5), or 1.4 mg rosuvastatin (n = 4) was administered to the portal vein, which was followed 120 minutes later by a second bolus of the statin coadministered with OATP perpetrator drug rifampicin (67.7 mg). After the first dose, all statins were rapidly cleared from the circulation (hepatic extraction ratio > 0.7) and excreted into the bile. Presence of human-specific atorvastatin metabolites confirmed the metabolic capacity of porcine livers. The predicted biliary clearance of rosuvastatin was found to be closer to the observed biliary clearance. A rank order of the DDI between the various systems upon coadministration with rifampicin could be observed: atorvastatin (AUC ratio 7.2) > rosuvastatin (AUC ratio 3.1) > pitavastatin (AUC ratio 2.6), which is in good agreement with the clinical DDI data. The results from this study demonstrated the applicability of using NMP of porcine livers as a novel preclinical model to study OATP-mediated DDI and its effect on hepatic clearance, biliary excretion, and plasma profile of drugs. SIGNIFICANCE STATEMENT: This study evaluated the use of normothermic machine perfusion (NMP) of porcine livers as a novel preclinical model to study hepatic clearance, biliary excretion, plasma (metabolite) profile of statins, and OATP-mediated DDI. Results showed that NMP of porcine livers is a reliable model to study OATP-mediated DDI. Overall, the rank order of DDI severity indicated in these experiments is in good agreement with clinical data, indicating the potential importance of this new ex vivo model in early drug discovery.

Characterization of Hepatic UDP-Glucuronosyltransferase Enzyme Abundance-Activity Correlations and Population Variability Using a Proteomics Approach and Comparison with Cytochrome P450 Enzymes

The expression of ten major drug-metabolizing UDP-glucuronosyltransferase (UGT) enzymes in a panel of 130 human hepatic microsomal samples was measured using a liquid chromatography-tandem mass spectrometry-based approach. Simultaneously, ten cytochromes P450 and P450 reductase were also measured, and activity-expression relationships were assessed for comparison. The resulting data sets demonstrated that, with the exception of UGT2B17, 10th to 90th percentiles of UGT expression spanned 3- to 8-fold ranges. These ranges were small relative to ranges of reported mean UGT enzyme expression across different laboratories. We tested correlation of UGT expression with enzymatic activities using selective probe substrates. A high degree of abundance-activity correlation (Spearman's rank correlation coefficient > 0.6) was observed for UGT1As (1A1, 3, 4, 6) and cytochromes P450. In contrast, protein abundance and activity did not correlate strongly for UGT1A9 and UGT2B enzymes (2B4, 7, 10, 15, and 17). Protein abundance was strongly correlated for UGTs 2B7, 2B10, and 2B15. We suggest a number of factors may contribute to these differences including incomplete selectivity of probe substrates, correlated expression of these UGT2B isoforms, and the impact of splice and polymorphic variants on the peptides used in proteomics analysis, and exemplify this in the case of UGT2B10. Extensive correlation analyses identified important criteria for validating the fidelity of proteomics and enzymatic activity approaches for assessing UGT variability, population differences, and ontogenetic changes. SIGNIFICANCE STATEMENT: Protein expression data allow detailed assessment of interindividual variability and enzyme ontogeny. This study has observed that expression and enzyme activity are well correlated for hepatic UGT1A enzymes and cytochromes P450. However, for the UGT2B family, caution is advised when assuming correlation of expression and activity as is often done in physiologically based pharmacokinetic modeling. This can be due to incomplete probe substrate specificities, but may also be related to presence of inactive UGT protein materials and the effect of splicing variations.

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.

Organic Cation Transporter 1 an Intestinal Uptake Transporter: Fact or Fiction?

Intestinal transporter proteins are known to affect the pharmacokinetics and in turn the efficacy and safety of many orally administered drugs in a clinically relevant manner. This knowledge is especially well-established for intestinal ATP-binding cassette transporters such as P-gp and BCRP. In contrast to this, information about intestinal uptake carriers is much more limited although many hydrophilic or ionic drugs are not expected to undergo passive diffusion but probably require specific uptake transporters. A transporter which is controversially discussed with respect to its expression, localization and function in the human intestine is the organic cation transporter 1 (OCT1). This review article provides an up-to-date summary on the available data from expression analysis as well as functional studies in vitro, animal findings and clinical observations. The current evidence suggests that OCT1 is expressed in the human intestine in small amounts (on gene and protein levels), while its cellular localization in the apical or basolateral membrane of the enterocytes remains to be finally defined, but functional data point to a secretory function of the transporter at the basolateral membrane. Thus, OCT1 should not be considered as a classical uptake transporter in the intestine but rather as an intestinal elimination pathway for cationic compounds from the systemic circulation.

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.

OATP2B1 - The underrated member of the organic anion transporting polypeptide family of drug transporters?

The organic anion transporting polypeptide 2B1 (OATP2B1) was one of the first cloned members of the SLCO family. However, its physiological and pharmacological role is still poorly understood, and object of a current debate on the transporter's relevance. Within this commentary, we summarize the data currently available on the transporter's expression and its substrates and highlight the strength and difficulties of the methods that have been applied to gather these data. The conclusion drawn from these findings was that OATP2B1 due to its intestinal expression is most likely involved in oral drug absorption of its substrate and therefore prone for interactions. This has been tested in in vivo drug interaction and/or pharmacogenetic studies. While some of these support the notion of OATP2B1 being of relevance in drug absorption, the pharmacogenetic findings are rather inconclusive. We will explain our thoughts why OATP2B1 may not influence the general systemic pharmacokinetic of certain substrates, but possibly local distribution processes, like the transfer across the blood-brain-barrier. Besides the pharmacokinetic aspects, there are data on endogenous molecules like coproporphyrins and sulfated steroids. Therefore, we will also highlight possible physiological roles of OATP2B1, which are driven by its expression pattern in the tubular cells of the kidney as well as its expression in the blood brain barrier. Finally we also deal with the advantages and disadvantages in the use of animal models to decipher the role of OATP2B1 in pharmacokinetics of its substrates and beyond.

Ontogeny of Hepatic Transporters and Drug-Metabolizing Enzymes in Humans and in Nonclinical Species

The liver represents a major eliminating and detoxifying organ, determining exposure to endogenous compounds, drugs, and other xenobiotics. Drug transporters (DTs) and drug-metabolizing enzymes (DMEs) are key determinants of disposition, efficacy, and toxicity of drugs. Changes in their mRNA and protein expression levels and associated functional activity between the perinatal period until adulthood impact drug disposition. However, high-resolution ontogeny profiles for hepatic DTs and DMEs in nonclinical species and humans are lacking. Meanwhile, increasing use of physiologically based pharmacokinetic (PBPK) models necessitates availability of underlying ontogeny profiles to reliably predict drug exposure in children. In addition, understanding of species similarities and differences in DT/DME ontogeny is crucial for selecting the most appropriate animal species when studying the impact of development on pharmacokinetics. Cross-species ontogeny mapping is also required for adequate translation of drug disposition data in developing nonclinical species to humans. This review presents a quantitative cross-species compilation of the ontogeny of DTs and DMEs relevant to hepatic drug disposition. A comprehensive literature search was conducted on PubMed Central: Tables and graphs (often after digitization) in original manuscripts were used to extract ontogeny data. Data from independent studies were standardized and normalized before being compiled in graphs and tables for further interpretation. New insights gained from these high-resolution ontogeny profiles will be indispensable to understand cross-species differences in maturation of hepatic DTs and DMEs. Integration of these ontogeny data into PBPK models will support improved predictions of pediatric hepatic drug disposition processes. SIGNIFICANCE STATEMENT: Hepatic drug transporters (DTs) and drug-metabolizing enzymes (DMEs) play pivotal roles in hepatic drug disposition. Developmental changes in expression levels and activities of these proteins drive age-dependent pharmacokinetics. This review compiles the currently available ontogeny profiles of DTs and DMEs expressed in livers of humans and nonclinical species, enabling robust interpretation of age-related changes in drug disposition and ultimately optimization of pediatric drug therapy.

Cross-species analysis of hepatic cytochrome P450 and transport protein expression

Most drugs and xenobiotics are metabolized in the liver. Amongst others, different cytochrome P450 (CYP) enzymes catalyze the metabolic conversion of foreign compounds, and various transport proteins are engaged in the excretion of metabolites from the hepatocytes. Inter-species and inter-individual differences in the hepatic levels and activities of drug-metabolizing enzymes and transporters result from genetic as well as from environmental factors, and play a decisive role in determining the pharmacokinetic properties of a compound in a given test system. To allow for a meaningful comparison of results from metabolism studies, it is, therefore, of utmost importance to know about the specific metabolic properties of the test systems, especially about the levels of metabolic enzymes such as the CYPs. Using a targeted proteomics approach, we, therefore, compared the hepatic levels of important CYP enzymes and transporters in different experimental systems in vivo and in vitro, namely Wistar rats, C57/Bl6 mice, mice humanized for the two xeno-sensing receptors PXR (pregnane-X-receptor) and CAR (constitutive androstane receptor), mice with human hepatocyte-repopulated livers, human HepaRG hepatocarcinoma cells, primary human hepatocytes, and human liver biopsies. In addition, the effects of xenobiotic inducers of drug metabolism on CYP enzymes and transporters were analyzed in selected systems. This study for the first time presents a comprehensive overview of similarities and differences in important drug metabolism-related proteins among the different experimental models.

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

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

Examination of Physiologically-Based Pharmacokinetic Models of Rosuvastatin

Physiologically‐based pharmacokinetic (PBPK) modeling is increasingly used to predict drug disposition and drug–drug interactions (DDIs). However, accurately predicting the pharmacokinetics of transporter substrates and transporter‐mediated DDIs (tDDIs) is still challenging. Rosuvastatin is a commonly used substrate probe in DDI risk assessment for new molecular entities (NMEs) that are potential organic anion transporting polypeptide 1B or breast cancer resistance protein transporter inhibitors, and as such, several rosuvastatin PBPK models have been developed to try to predict the clinical DDI and support NME drug labeling. In this review, we examine five representative PBPK rosuvastatin models, discuss common challenges that the models have come across, and note remaining gaps. These shared learnings will help with the continuing efforts of rosuvastatin model validation, provide more information to understand transporter‐mediated drug disposition, and increase confidence in tDDI prediction.

Mass spectrometry-based abundance atlas of ABC transporters in human liver, gut, kidney, brain and skin

ABC transporters (ATP-binding cassette transporter) traffic drugs and their metabolites across membranes, making ABC transporter expression levels a key factor regulating local drug concentrations in different tissues and individuals. Yet, quantification of ABC transporters remains challenging because they are large and low-abundance transmembrane proteins. Here, we analysed 200 samples of crude and membrane-enriched fractions from human liver, kidney, intestine, brain microvessels and skin, by label-free quantitative mass spectrometry. We identified 32 (out of 48) ABC transporters: ABCD3 was the most abundant in liver, whereas ABCA8, ABCB2/TAP1 and ABCE1 were detected in all tissues. Interestingly, this atlas unveiled that ABCB2/TAP1 may have TAP2-independent functions in the brain and that biliary atresia (BA) and control livers have quite different ABC transporter profiles. We propose that meaningful biological information can be derived from a direct comparison of these data sets.

Proteomic Evaluation of Plasma Membrane Fraction Prepared from a Mouse Liver and Kidney Using a Bead Homogenizer: Enrichment of Drug-Related Transporter Proteins

Quantifying the protein levels of drug transporters in plasma membrane fraction helps elucidate the function of these transporters. In this study, we conducted a proteomic evaluation of enriched drug-related transporter proteins in plasma membrane fraction prepared from mouse liver and kidney tissues using the membrane protein extraction kit and a bead homogenizer. Crude and plasma membrane fractions were prepared using either the Dounce or bead homogenizer, and protein levels were determined using quantitative proteomics. In liver tissues, the plasma membrane fractions were more enriched in transporter proteins than the crude membrane fractions; the average enrichment ratios of plasma-to-crude membrane fractions were 3.31 and 6.93 using the Dounce and bead homogenizers, respectively. The concentrations of transporter proteins in plasma membrane fractions determined using the bead homogenizer were higher than those determined using the Dounce homogenizer. Meanwhile, in kidney tissues, the plasma membrane fractions were enriched in transporters localized in the brush-border membrane to the same degree for both the homogenizers; however, the membrane fractions obtained using either homogenizer were not enriched in Na^+/K⁺-ATPase and transporters localized in the basolateral membrane. These results indicate that fractionation, using the bead homogenizer, yielded transporter-enriched plasma membrane fractions from mouse liver and kidney tissues; however, no enrichment of basolateral transporters was observed in plasma membrane fractions prepared from kidney tissues.

Dilemmas With Absolute Quantification of Pharmacologically Relevant Proteins Using Mass Spectrometry

Determination of abundances of proteins involved in uptake, distribution, metabolism and excretion of xenobiotics is a prerequisite to understand and predict elimination mechanisms in tissue. Mass spectrometry promises simple and accurate measurements of individual proteins in complex mixtures using isotopically labeled peptide standards. However, comparisons of measurements performed in different laboratories have shown considerable discrepancies in the data generated. Even when very similar approaches are compared, the results differ significantly. An alternative method of measuring protein titers is global proteomics. Depending on sample type, this allows quantification of hundreds to thousands of proteins in a single analysis. It enables system-wide insights by providing protein copy numbers and cell sizes. Regardless of differences, the workflows of both the labeled standard-based and the proteomic approach share several steps. Each can be critical. Selection of optimal techniques is the prerequisite for accurate and reproducible protein quantification.

A Tribute to Professor Per Artursson - Scientist, Explorer, Mentor, Innovator, and Giant in Pharmaceutical Research

This issue of the Journal of Pharmaceutical Sciences is dedicated to Professor Per Artursson and the groundbreaking contributions he has made and continues to make in the Pharmaceutical Sciences. Per is one of the most cited researchers in his field, with more than 30,000 citations and an h-index of 95 as of September 2020. Importantly, these citations are distributed over the numerous fields he has explored, clearly showing the high impact the research has had on the discipline. We provide a short portrait of Per, with emphasis on his personality, driving forces and the inspirational sources that shaped his career as a world-leading scientist in the field. He is a curious scientist who deftly moves between disciplines and has continued to innovate, expand boundaries, and profoundly impact the pharmaceutical sciences throughout his career. He has developed new tools and provided insights that have significantly contributed to today's molecular and mechanistic approaches to research in the fields of intestinal absorption, cellular disposition, and exposure-efficacy relationships of pharmaceutical drugs. We want to celebrate these important contributions in this special issue of the Journal of Pharmaceutical Sciences in Per's honor.

Proteomics-Informed Prediction of Rosuvastatin Plasma Profiles in Patients With a Wide Range of Body Weight

Rosuvastatin is a frequently used probe to study transporter‐mediated hepatic uptake. Pharmacokinetic models have therefore been developed to predict transporter impact on rosuvastatin disposition in vivo. However, the interindividual differences in transporter concentrations were not considered in these models, and the predicted transporter impact was compared with historical in vivo data. In this study, we investigated the influence of interindividual transporter concentrations on the hepatic uptake clearance of rosuvastatin in 54 patients covering a wide range of body weight. The 54 patients were given an oral dose of rosuvastatin the day before undergoing gastric bypass or cholecystectomy, and pharmacokinetic (PK) parameters were established from each patient’s individual time‐concentration profiles. Liver biopsies were sampled from each patient and their individual hepatic transporter concentrations were quantified. We combined the transporter concentrations with in vitro uptake kinetics determined in HEK293‐transfected cells, and developed a semimechanistic model with a bottom‐up approach to predict the plasma concentration profiles of the single dose of rosuvastatin in each patient. The predicted PK parameters were evaluated against the measured in vivo plasma PKs from the same 54 patients. The developed model predicted the rosuvastatin PKs within two‐fold error for rosuvastatin area under the plasma concentration versus time curve (AUC; 78% of the patients; average fold error (AFE): 0.96), peak plasma concentration (C_max; 76%; AFE: 1.05), and terminal half‐life (t_1/2; 98%; AFE: 0.89), and captured differences in the rosuvastatin PKs in patients with the OATP1B1 521T<C polymorphism. This demonstrates that hepatic uptake clearance determined in transfected cell lines, together with proteomics scaling, provides a useful tool for prediction models, without the need for empirical scaling factors.

The study of ultrasound and iontophoresis on oxaprozin transdermal penetration using surface-enhanced Raman spectroscopy

The potential for physicochemical driving forces facilitating topical transport of the lipid-soluble drug oxaprozin (OXA) was investigated using surface-enhanced Raman spectroscopy (SERS) in this study. Azone, iontophoresis (IP), and sonophoresis (SP) were combined and performed on mouse skin for the OXA transdermal penetration, and the synergistic effect was analyzed using Raman spectroscopy. The data of characteristic peak intensity were processed with overlapping peak resolving and standard normalization. The results showed that Azone promoted the transdermal penetration of OXA (5.9-fold greater than the OXA concentration of normal penetration); SP enhanced OXA transdermal penetration (5.5-fold); IP enhanced OXA transdermal penetration (4.2-fold); the combined application of Azone and SP (Azone+SP) and SP+IP can improve the enhancement coefficient of OXA transdermal penetration (8.4-fold and 6.1-fold, > 5.9, > 5.5, > 4.2), and their combined application has a synergistic effect; Azone+IP does not have a synergistic effect while the enhancement coefficient of Azone+IP (5.3-fold, < 5.9) and Azone+SP+IP (7.2-fold, < 8.4) was slightly reduced. As for the drug OXA, Azone+SP is an effective method of transdermal penetration.

Expanding the Efflux In Vitro Assay Toolbox: A CRISPR-Cas9 Edited MDCK Cell Line with Human BCRP and Completely Lacking Canine MDR1

The Breast Cancer Resistance Protein (BCRP) is a key transporter in drug efflux and drug-drug interactions. However, endogenous expression of Multidrug Resistance Protein 1 (MDR1) confounds the interpretation of BCRP-mediated transport in in vitro models. Here we used a CRISPR-Cas9 edited Madin-Darby canine kidney (MDCK) II cell line (MDCK^cMDR1-KO) for stable expression of human BCRP (hBCRP) with no endogenous canine MDR1 (cMDR1) expression (MDCK-hBCRP^cMDR1-KO). Targeted quantitative proteomics verified expression of hBCRP, and global analysis of the entire proteome corroborated no or very low background expression of other drug transport proteins or metabolizing enzymes. This new cell line, had similar proteome like MDCK^cMDR1-KO and a previously established, corresponding cell line overexpressing human MDR1 (hMDR1), MDCK-hMDR1^cMDR1-KO. Functional studies with MDCK-hBCRP^cMDR1-KO confirmed high hBCRP activity. The MDCK-hBCRP^cMDR1-KO cell line together with the MDCK-hMDR1^cMDR1-KO easily and accurately identified shared or specific substrates of the hBCRP and the hMDR1 transporters. These cell lines offer new, improved in vitro tools for the assessment of drug efflux and drug-drug interactions in drug development.

Impact of iron oxide nanoparticles on xenobiotic metabolism in HepaRG cells

Iron oxide nanoparticles are used in various industrial fields, as a tool in biomedicine as well as in food colorants, and can therefore reach human metabolism via oral uptake or injection. However, their effects on the human body, especially the liver as one of the first target organs is still under elucidation. Here, we studied the influence of different representative iron oxide materials on xenobiotic metabolism of HepaRG cells. These included four iron oxide nanoparticles, one commercially available yellow food pigment (E172), and non-particulate ionic control FeSO₄. The nanoparticles had different chemical and crystalline structures and differed in size and shape and were used at a concentration of 50 µg Fe/mL. We found that various CYP enzymes were downregulated by some but not all iron oxide nanoparticles, with the Fe₃O₄-particle, both γ-Fe₂O₃-particles, and FeSO₄ exhibiting the strongest effects, the yellow food pigment E172 showing a minor effect and an α-Fe₂O₃ nanoparticle leading to almost no inhibition of phase I machinery. The downregulation was seen at the mRNA, protein expression, and activity levels. Thereby, no dependency on the size or chemical structure was found. This underlines the difficulty of the grouping of nanomaterials regarding their physiological impact, suggesting that every iron oxide nanoparticle species needs to be evaluated in a case-by-case approach.

P-Glycoprotein Inhibition Exacerbates Paclitaxel Neurotoxicity in Neurons and Patients With Cancer

Paclitaxel-induced peripheral neuropathy (PIPN) is a common and dose-limiting adverse event. The role of P-glycoprotein (P-gp) in the neuronal efflux of paclitaxel was assessed using a translational approach. SH-SY5Y cells were differentiated to neurons and paclitaxel toxicity in the absence and presence of a P-gp inhibitor was determined. Paclitaxel caused marked dose-dependent toxicity in SH-SY5Y-derived neurons. Paclitaxel neurotoxicity was exacerbated with concomitant P-gp inhibition by valspodar and verapamil, consistent with increased intracellular accumulation of paclitaxel. Patients with cancer treated with paclitaxel and P-gp inhibitors had a 2.4-fold (95% confidence interval (CI) 1.3-4.3) increased risk of peripheral neuropathy-induced dose modification while a 4.7-fold (95% CI 1.9-11.9) increased risk for patients treated with strong P-gp inhibitors was observed, and a 7.0-fold (95% CI 2.3-21.5) increased risk in patients treated with atorvastatin. Atorvastatin also increased neurotoxicity by paclitaxel in SH-SY5Y-derived neurons. Clinicians should be aware that comedication with P-gp inhibitors may lead to increased risk of PIPN.

Induction and repression effects on CYP and transporter protein abundance by azole mixture uptake in rat liver

Detection of mixture effects is a major challenge in current experimental and regulatory toxicology. Robust markers are needed that are easy to quantify and responsive to chemical stressors in a broad dose range. Several hepatic enzymes and proteins related to drug metabolism like cytochrome-P-450 (CYP) enzymes and transporters have been shown to be responsive to pesticide active substances in a broad dose range and are therefore good candidates to be used as markers for mixture toxicity. Even though they can be well quantified at the mRNA level, quantification on the protein level is challenging because most of these proteins are membrane bound. Here we report the development of mass spectrometry-based assays using triple-x-proteomics (TXP) antibodies in combination with targeted selected ion monitoring (tSIM) to quantify changes of protein levels due to exposure to mixtures of pesticide active substances. Our results indicate that changes on the protein level of CYP1A1, ABCB2, ABCC3 are in line with observations on the mRNA and enzyme activity level and are indicative of mixture effects. Therefore, the tests are promising to reveal effects by chemical mixture effects in toxicological studies in rats.