Ethnic variability in the expression of hepatic drug transporters: absolute quantification by an optimized targeted quantitative proteomic approach

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.

Development of a UPLC-MRM-based targeted proteomic method to profile subcellular organelle marker proteins from human liver tissues

Subcellular organelles have long been an interest in biochemical research and drug development as the isolation of those organelles can help to probe protein functions and elucidate drug disposition within the cell. Usually, the purity of isolated subcellular organelle fractions was determined using immunoblot analysis of subcellular organelle marker proteins, which can be labor-intensive and lack reproducibility due to antibody batch-to-batch variability. As such, a higher throughput and more robust method is needed. Here, a UPLC-MRM-based targeted proteomic method was developed for a panel of human organelle marker proteins and used to profile a series of sucrose fractions isolated from the protein extract of human liver tissues. The method was validated by comparing to the traditional immunoblot and determining subcellular localization of three case study proteins (CYP3A4, FcRn, and β2M) pertaining to the disposition of small molecule and biologic drugs. All three case study proteins were co-enriched with their corresponding subcellular protein marker, and complete recoveries were achieved from isolated fractions. This newly developed MRM method for the panel of human organelle marker proteins can potentially accelerate future intracellular drug disposition analysis and facilitate subcellular organelle quality assessment.

Functional in vitro characterization of SLCO1B1 variants and simulation of the clinical pharmacokinetic impact of impaired OATP1B1 function

Organic Anion Transporting Polypeptide 1B1 is important to the hepatic elimination and distribution of many drugs. If OATP1B1 function is decreased, it can increase plasma exposure of e.g. several statins leading to increased risk of muscle toxicity. First, we examined the impact of three naturally occurring rare variants and the frequent SLCO1B1 c.388A>G variant on in vitro transport activity with cellular uptake assay using two substrates: 2', 7'-dichlorofluorescein (DCF) and rosuvastatin. Secondly, LC-MS/MS based quantitative targeted absolute proteomics measured the OATP1B1 protein abundance in crude membrane fractions of HEK293 cells over-expressing these SNVs. Additionally, we simulated the effect of impaired OATP1B1 function on rosuvastatin pharmacokinetics to estimate the need for genotype-guided dosing. R57Q impaired DCF and rosuvastatin transport significantly yet did not change protein expression considerably, while N130D and N151S did not alter activity but increased protein expression. R253Q did not change protein expression but reduced DCF uptake and increased rosuvastatin K_m. Based on pharmacokinetics simulations, doses of 30 mg (with 50% OATP1B1 function) and 20 mg (with 0% OATP1B1 function) result in plasma exposure similar to 40 mg dose (with 100% OATP1B1 function). Therefore dose reductions might be considered to avoid increased plasma exposure caused by function-impairing OATP1B1 genetic variants, such as R57Q.

Assessment of Liver Function With MRI: Where Do We Stand?

Liver disease and hepatocellular carcinoma (HCC) have become a global health burden. For this reason, the determination of liver function plays a central role in the monitoring of patients with chronic liver disease or HCC. Furthermore, assessment of liver function is important, e.g., before surgery to prevent liver failure after hepatectomy or to monitor the course of treatment. Liver function and disease severity are usually assessed clinically based on clinical symptoms, biopsy, and blood parameters. These are rather static tests that reflect the current state of the liver without considering changes in liver function. With the development of liver-specific contrast agents for MRI, noninvasive dynamic determination of liver function based on signal intensity or using T1 relaxometry has become possible. The advantage of this imaging modality is that it provides additional information about the vascular structure, anatomy, and heterogeneous distribution of liver function. In this review, we summarized and discussed the results published in recent years on this technique. Indeed, recent data show that the T1 reduction rate seems to be the most appropriate value for determining liver function by MRI. Furthermore, attention has been paid to the development of automated tools for image analysis in order to uncover the steps necessary to obtain a complete process flow from image segmentation to image registration to image analysis. In conclusion, the published data show that liver function values obtained from contrast-enhanced MRI images correlate significantly with the global liver function parameters, making it possible to obtain both functional and anatomic information with a single modality.

Inter-ethnic differences in pharmacokinetics-is there more that unites than divides?

Inter-ethnic variability in pharmacokinetics (PK) has been attributed to several factors ranging from genetic to environmental. It is not clear how current teaching in higher education (HE) reflects what published literature suggests on this subject. This study aims to gain insights into current knowledge about inter-ethnic differences in PK based on reports from published literature and current teaching practices in HE. A systematic literature search was conducted on PubMed and Scopus to identify suitable literature to be reviewed. Insights into inter-ethnic differences in PK teaching among educators in HE and industry were determined using a questionnaire. Thirty-one percent of the studies reviewed reported inter-ethnic differences in PK, of these, 37% of authors suggested genetic polymorphism as possible explanation for the inter-ethnic differences observed. Other factors authors proposed included diet and weight differences between ethnicities. Most respondents (80%) who taught inter-ethnic difference in PK attributed inter-ethnic differences to genetic polymorphism. While genetic polymorphism is one source of variability in PK, the teaching of genetic polymorphism is better associated with interindividual variabilities rather than inter-ethnic differences in PK as there are no genes with PK implications specific to any one ethnic group. Nongenetic factors such as diet, weight, and environmental factors, should be highlighted as potential sources of interindividual variation in the PK of drugs.

Physiologically Based Modeling of the Effect of Physiological and Anthropometric Variability on Indocyanine Green Based Liver Function Tests

Accurate evaluation of liver function is a central task in hepatology. Dynamic liver function tests (DLFT) based on the time-dependent elimination of a test substance provide an important tool for such a functional assessment. These tests are used in the diagnosis and monitoring of liver disease as well as in the planning of hepatobiliary surgery. A key challenge in the evaluation of liver function with DLFTs is the large inter-individual variability. Indocyanine green (ICG) is a widely applied test compound used for the evaluation of liver function. After an intravenous administration, pharmacokinetic (PK) parameters are calculated from the plasma disappearance curve of ICG which provide an estimate of liver function. The hepatic elimination of ICG is affected by physiological factors such as hepatic blood flow or binding of ICG to plasma proteins, anthropometric factors such as body weight, age, and sex, or the protein amount of the organic anion-transporting polypeptide 1B3 (OATP1B3) mediating the hepatic uptake of ICG. Being able to account for and better understand these various sources of inter-individual variability would allow to improve the power of ICG based DLFTs and move toward an individualized evaluation of liver function. Within this work we systematically analyzed the effect of various factors on ICG elimination by the means of computational modeling. For the analysis, a recently developed and validated physiologically based pharmacokinetics (PBPK) model of ICG distribution and hepatic elimination was utilized. Key results are (i) a systematic analysis of the variability in ICG elimination due to hepatic blood flow, cardiac output, OATP1B3 abundance, liver volume, body weight and plasma bilirubin level; (ii) the evaluation of the inter-individual variability in ICG elimination via a large in silico cohort of n = 100,000 subjects based on the NHANES cohort with special focus on stratification by age, sex, and body weight; (iii) the evaluation of the effect of various degrees of cirrhosis on variability in ICG elimination. The presented results are an important step toward individualizing liver function tests by elucidating the effects of confounding physiological and anthropometric parameters in the evaluation of liver function via ICG.

PBPK Model of Coproporphyrin I: Evaluation of the Impact of SLCO1B1 Genotype, Ethnicity, and Sex on its Inter-Individual Variability

Coproporphyrin I (CPI) is an endogenous biomarker of OATP1B activity and associated drug-drug interactions. In this study, a minimal physiologically-based pharmacokinetic model was developed to investigate the impact of OATP1B1 genotype (c.521T>C), ethnicity, and sex on CPI pharmacokinetics and interindividual variability in its baseline. The model implemented mechanistic descriptions of CPI hepatic transport between liver blood and liver tissue and renal excretion. Key model parameters (e.g., endogenous CPI synthesis rate, and CPI hepatic uptake clearance) were estimated by fitting the model simultaneously to three independent CPI clinical datasets (plasma and urine data) obtained from white (n = 16, men and women) and Asian-Indian (n = 26, all men) subjects, with c.521 variants (TT, TC, and CC). The optimized CPI model successfully described the observed data using c.521T>C genotype, ethnicity, and sex as covariates. CPI hepatic active was 79% lower in 521CC relative to the wild type and 42% lower in Asian-Indians relative to white subjects, whereas CPI synthesis was 23% higher in male relative to female subjects. Parameter sensitivity analysis showed marginal impact of the assumption of CPI synthesis site (blood or liver), resulting in comparable recovery of plasma and urine CPI data. Lower magnitude of CPI-drug interaction was simulated in 521CC subjects, suggesting the risk of underestimation of CPI-drug interaction without prior OATP1B1 genotyping. The CPI model incorporates key covariates contributing to interindividual variability in its baseline and highlights the utility of the CPI modeling to facilitate the design of prospective clinical studies to maximize the sensitivity of this biomarker.

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.

Quantification of Neonatal Fc Receptor and Beta-2 Microglobulin in Human Liver Tissues by Ultraperformance Liquid Chromatography-Multiple Reaction Monitoring-based Targeted Quantitative Proteomics for Applications in Biotherapeutic Physiologically-based Pharmacokinetic Models

Neonatal Fc receptor (FcRn) and beta-2 microglobulin (β2M) play an important role in transporting maternal IgG to fetuses, maintaining the homeostasis of IgG and albumin in human body, and prolonging the half-life of IgG- or albumin-based biotherapeutics. Little is known about the influence of age, gender and race, and interindividual variability of human FcRn and β2M on the protein level. In this study, an ultraperformance liquid chromatography-multiple reaction monitoring mass spectrometry-based targeted quantitative proteomic method was developed and optimized for the quantification of human FcRn and β2M. Among the 39 human livers studied (age 13-80 years), the mean (±S.D.) concentrations of FcRn and β2M were 147 (±39) and 1250 (±460) pmol/g of liver tissue, respectively. A four-fold interindividual variability (63-243 pmol/g of liver tissue) was observed for the hepatic FcRn concentration. A moderate correlation was found between the hepatic β2M and FcRn expression levels. Influences of age, gender, and race on the hepatic expression of FcRn and β2M were evaluated. The findings from this study may aid the development of physiologically-based pharmacokinetic models that incorporate empirical FcRn tissue concentrations and interindividual variabilities, and the development of personalized dosing of biopharmaceuticals. SIGNIFICANCE STATEMENT: This is the first study to evaluate the influence of age, gender, and race on the expression of neonatal Fc receptor (FcRn) and beta-2 microglobulin (β2M) and their interindividual variability in human livers. This study describes a validated ultraperformance liquid chromatography-multiple reaction monitoring-based targeted quantitative proteomic method for quantifying human FcRn and β2M in biological tissues. Results from this study may aid current development of physiologically-based pharmacokinetic models for biotherapeutics, where FcRn plays a significant role in clearance mechanism, and its expression level and interindividual variability are largely unknown.

Microengineered systems with iPSC-derived cardiac and hepatic cells to evaluate drug adverse effects

Hepatic and cardiac drug adverse effects are among the leading causes of attrition in drug development programs, in part due to predictive failures of current animal or in vitro models. Hepatocytes and cardiomyocytes differentiated from human induced pluripotent stem cells (iPSCs) hold promise for predicting clinical drug effects, given their human-specific properties and their ability to harbor genetically determined characteristics that underlie inter-individual variations in drug response. Currently, the fetal-like properties and heterogeneity of hepatocytes and cardiomyocytes differentiated from iPSCs make them physiologically different from their counterparts isolated from primary tissues and limit their use for predicting clinical drug effects. To address this hurdle, there have been ongoing advances in differentiation and maturation protocols to improve the quality and use of iPSC-differentiated lineages. Among these are in vitro hepatic and cardiac cellular microsystems that can further enhance the physiology of cultured cells, can be used to better predict drug adverse effects, and investigate drug metabolism, pharmacokinetics, and pharmacodynamics to facilitate successful drug development. In this article, we discuss how cellular microsystems can establish microenvironments for these applications and propose how they could be used for potentially controlling the differentiation of hepatocytes or cardiomyocytes. The physiological relevance of cells is enhanced in cellular microsystems by simulating properties of tissue microenvironments, such as structural dimensionality, media flow, microfluidic control of media composition, and co-cultures with interacting cell types. Recent studies demonstrated that these properties also affect iPSC differentiations and we further elaborate on how they could control differentiation efficiency in microengineered devices. In summary, we describe recent advances in the field of cellular microsystems that can control the differentiation and maturation of hepatocytes and cardiomyocytes for drug evaluation. We also propose how future research with iPSCs within engineered microenvironments could enable their differentiation for scalable evaluations of drug effects.

Interindividual Diversity in Expression of Organic Anion Uptake Transporters in Normal and Cirrhotic Human Liver

The liver plays an essential role in removing endogenous and exogenous compounds from the circulation. This function is mediated by specific transporters, including members of the family of organic anion transport proteins (OATPs) and the Na⁺-taurocholate transporting polypeptide (NTCP). In the present study, transporter protein expression was determined in liver samples from patients with cirrhosis or controls without liver disease. Five transporters (OATP1A2, OATP1B1, OATP1B3, OATP2B1, and NTCP) were studied. Transporter content in homogenates of human liver was quantified on western blots probed with transporter-specific antibodies in which a calibrated green fluorescent protein-tagged transporter standard was included. Liver samples from 21 patients with cirrhosis (hepatitis C in 17 and alcohol abuse in 4) and 17 controls without liver disease were analyzed. Expression of each of the transporters had a large spread, varying by an order of magnitude in cirrhotic and control livers. OATP1B1 was the most abundant transporter in controls (P < 0.01) but was significantly lower in cirrhotic livers as was NTCP expression (P < 0.01). There was little difference in transporter expression with respect to age or sex. Despite the large variability in transporter expression within a group, analysis in individuals showed that those with high or low expression of one transporter had a similar magnitude in expression of the others. Conclusion: Differences in transporter expression could explain unanticipated heterogeneity of drug transport and metabolism in individuals with and without liver disease.

Protein Expression and Functional Relevance of Efflux and Uptake Drug Transporters at the Blood-Brain Barrier of Human Brain and Glioblastoma

The knowledge of transporter protein expression and function at the human blood-brain barrier (BBB) is critical to prediction of drug BBB penetration and design of strategies for improving drug delivery to the brain or brain tumor. This study determined absolute transporter protein abundances in isolated microvessels of human normal brain (N = 30), glioblastoma (N = 47), rat (N = 10) and mouse brain (N = 10), and cell membranes of MDCKII cell lines, using targeted proteomics. In glioblastoma microvessels, efflux transporters (ABCB1 and ABCG2), monocarboxylate transporter 1 (MCT1), glucose transporter 1 (GLUT1), sodium-potassium pump (Na/K ATPase), and Claudin-5 protein levels were significantly reduced, while large neutral amino acid transporter 1 (LAT1) was increased and GLU3 remained the same, as compared with human normal brain microvessels. ABCC4, OATP1A2, OATP2B1, and OAT3 were undetectable in microvessels of both human brain and glioblastoma. Species difference in BBB transporter abundances was noted. Cellular permeability experiments and modeling simulations suggested that not a single apical uptake transporter but a vectorial transport system consisting of an apical uptake transporter and basolateral efflux mechanism was required for efficient delivery of poor transmembrane permeability drugs from the blood to brain.

Toward a Consensus on Applying Quantitative Liquid Chromatography-Tandem Mass Spectrometry Proteomics in Translational Pharmacology Research: A White Paper

Quantitative translation of information on drug absorption, disposition, receptor engagement, and drug-drug interactions from bench to bedside requires models informed by physiological parameters that link in vitro studies to in vivo outcomes. To predict in vivo outcomes, biochemical data from experimental systems are routinely scaled using protein quantity in these systems and relevant tissues. Although several laboratories have generated useful quantitative proteomic data using state-of-the-art mass spectrometry, no harmonized guidelines exit for sample analysis and data integration to in vivo translation practices. To address this gap, a workshop was held on September 27 and 28, 2018, in Cambridge, MA, with 100 experts attending from academia, the pharmaceutical industry, and regulators. Various aspects of quantitative proteomics and its applications in translational pharmacology were debated. A summary of discussions and best practices identified by this expert panel are presented in this "White Paper" alongside unresolved issues that were outlined for future debates.

Explaining Ethnic Variability of Transporter Substrate Pharmacokinetics in Healthy Asian and Caucasian Subjects with Allele Frequencies of OATP1B1 and BCRP: A Mechanistic Modeling Analysis

Background

Ethnic variability in the pharmacokinetics of organic anion transporting polypeptide (OATP) 1B1 substrates has been observed, but its basis is unclear. A previous study hypothesizes that, without applying an intrinsic ethnic variability in transporter activity, allele frequencies of transporters cannot explain observed ethnic variability in pharmacokinetics. However, this hypothesis contradicts the data collected from compounds that are OATP1B1 substrates but not breast cancer resistance protein (BCRP) substrates.

Objective

The objective of this study is to evaluate a hypothesis that is physiologically reasonable and more consistent with clinical observations.

Methods

We evaluated if allele frequencies of two transporters (OATP1B1 and BCRP) are key contributors to ethnic variability. In this hypothesis, the same genotype leads to the same activity independent of ethnicity, in contrast to the previous hypothesis of intrinsic ethnic variability in OATP1B1 activity. As a validation, we perform mechanistic pharmacokinetic modeling for SLCO1B1 (encoding OATP1B1) and ABCG2 (encoding BCRP) genotyped pharmacokinetic data from 18 clinical studies with healthy Caucasian and/or Asian subjects.

Results

Simulations based on the current hypothesis reasonably describe SLCO1B1 and ABCG2 genotyped pharmacokinetic time course data for five transporter substrates (atorvastatin, pitavastatin, pravastatin, repaglinide, and rosuvastatin) in Caucasian and Asian populations.

Conclusion

This hypothesis covers the observations that can (e.g., ethnic differences in rosuvastatin pharmacokinetics) or cannot (e.g., lack of differences for pitavastatin pharmacokinetics) be explained by the previous hypothesis. It helps to characterize sources of ethnic variability and provides a foundation for predicting ethnic variability in transporter substrate pharmacokinetics.

Electronic supplementary material

The online version of this article (doi:10.1007/s40262-017-0568-7) contains supplementary material, which is available to authorized users.

A repository of protein abundance data of drug metabolizing enzymes and transporters for applications in physiologically based pharmacokinetic (PBPK) modelling and simulation

Population factors such as age, gender, ethnicity, genotype and disease state can cause inter-individual variability in pharmacokinetic (PK) profile of drugs. Primarily, this variability arises from differences in abundance of drug metabolizing enzymes and transporters (DMET) among individuals and/or groups. Hence, availability of compiled data on abundance of DMET proteins in different populations can be useful for developing physiologically based pharmacokinetic (PBPK) models. The latter are routinely employed for prediction of PK profiles and drug interactions during drug development and in case of special populations, where clinical studies either are not feasible or have ethical concerns. Therefore, the main aim of this work was to develop a repository of literature-reported DMET abundance data in various human tissues, which included compilation of information on sample size, technique(s) involved, and the demographic factors. The collation of literature reported data revealed high inter-laboratory variability in abundance of DMET proteins. We carried out unbiased meta-analysis to obtain weighted mean and percent coefficient of variation (%CV) values. The obtained %CV values were then integrated into a PBPK model to highlight the variability in drug PK in healthy adults, taking lamotrigine as a model drug. The validated PBPK model was extrapolated to predict PK of lamotrigine in paediatric and hepatic impaired populations. This study thus exemplifies importance of the DMET protein abundance database, and use of determined values of weighted mean and %CV after meta-analysis in PBPK modelling for the prediction of PK of drugs in healthy and special populations.

Effects of Food and Gender on Pharmacokinetics of Rosuvastatin in a Chinese Population Based on 4 Bioequivalence Studies

The effects of food and gender on the pharmacokinetics of rosuvastatin in healthy Chinese subjects were investigated from 4 bioequivalence studies. These studies were designed as randomized, open-label, and 2-period crossover in both fasting and fed states. A total of 204 subjects were enrolled, 134 men and 70 women. These subjects received a single oral 10-mg dose of rosuvastatin with a 7-day washout between 2 periods. The plasma concentrations were determined using a validated liquid chromatography tandem mass spectrometry method, and pharmacokinetic parameters were calculated by noncompartmental methods. Compared with the fasting condition, administration after a high-fat and high-calorie meal resulted in an approximately 40% reduction of rosuvastatin exposure and a near 50% decrease in absorption rate. Moreover, the apparent clearance was significantly greater in the fed state than that in the fasting state. It was noted that the adverse events incidence is increased by approximately 30% in the fasting state; however, no serious adverse events were observed. Additionally, small differences in pharmacokinetic characteristics were found between male and female subjects. Food effect might be considered for optimal effectiveness and safety of rosuvastatin therapy.

Protein Abundance of Clinically Relevant Drug Transporters in the Human Liver and Intestine: A Comparative Analysis in Paired Tissue Specimens

Bioavailability of orally administered drugs is partly determined by function of drug transporters in the liver and intestine. Therefore, we explored adenosine triphosphate-binding cassette (ABC) and solute carriers family transporters expression (quantitative polymerase chain reaction) and protein abundance (liquid chromatography tandem mass spectrometry (LC-MS/MS)) in human liver and duodenum, jejunum, ileum, and colon in paired tissue specimens from nine organ donors. The transporter proteins were detected in the liver (permeability-glycoprotein (P-gp), multidrug resistance protein (MRP)2, MRP3, breast cancer resistance protein (BCRP), organic anion-transporting polypeptide (OATP)1B1, OATP1B3, OATP2B1, organic cation transporter (OCT)1, OCT3, organic anion transporter 2, Na+-taurocholate cotransporting polypeptide, monocarboxylate transporter (MCT)1, and multidrug and toxin extrusion 1) and the intestine (P-gp, multidrug-resistance protein (MRP)2, MRP3, MRP4, BCRP, OATP2B1, OCT1, apical sodium-bile acid transporter (only ileum), MCT1, and peptide transporter (PEPT1)). Significantly higher hepatic gene expression and protein abundance of ABCC2/MRP2, SLC22A1/OCT1, and SLCO2B1/OATP2B1 were found, as compared to all intestinal segments. No correlations between hepatic and small intestinal protein levels were observed. These observations provide a description of drug transporters distribution without the impact of interindividual variability bias and may help in construction of superior physiologically based pharmacokinetic and humanized animal models.

Advancing Predictions of Tissue and Intracellular Drug Concentrations Using In Vitro, Imaging and Physiologically Based Pharmacokinetic Modeling Approaches

This white paper examines recent progress, applications, and challenges in predicting unbound and total tissue and intra/subcellular drug concentrations using in vitro and preclinical models, imaging techniques, and physiologically based pharmacokinetic (PBPK) modeling. Published examples, regulatory submissions, and case studies illustrate the application of different types of data in drug development to support modeling and decision making for compounds with transporter-mediated disposition, and likely disconnects between tissue and systemic drug exposure. The goals of this article are to illustrate current best practices and outline practical strategies for selecting appropriate in vitro and in vivo experimental methods to estimate or predict tissue and plasma concentrations, and to use these data in the application of PBPK modeling for human pharmacokinetic (PK), efficacy, and safety assessment in drug development.

Quantification of Hepatic Organic Anion Transport Proteins OAT2 and OAT7 in Human Liver Tissue and Primary Hepatocytes

Organic anion transporter (OAT) 2 and OAT7 were recently shown to be involved in the hepatic uptake of drugs; however, there is limited understanding of the population variability in the expression of these transporters in liver. There is also a need to derive relative expression-based scaling factors (REFs) that can be used to bridge in vitro functional data to the in vivo drug disposition. To this end, we quantified OAT2 and OAT7 surrogate peptide abundance in a large number of human liver tissue samples ( n = 52), as well as several single-donor cryopreserved human hepatocyte lots ( n = 30) by a novel, validated liquid chromatography tandem mass spectrometry (LC-MS/MS) method. The average surrogate peptide expression of OAT2 and OAT7 in the liver samples was 1.52 ± 0.57 and 4.63 ± 1.58 fmol/μg membrane protein, respectively. While we noted statistically significant differences ( p < 0.05) in hepatocyte and liver tissue abundances for both OAT2 and OAT7, the differences were relatively small (1.8- and 1.5-fold difference in median values, respectively). Large interindividual variability was noted in the hepatic expression of OAT2 (16-fold in liver tissue and 23-fold in hepatocytes). OAT7, on the other hand, showed less interindividual variability (4-fold) in the livers, but high variability for the hepatocyte lots (27-fold). A significant positive correlation in OAT2 and OAT7 expression was observed, but expression levels were neither associated with age nor sex. In conclusion, our data suggest marked interindividual variability in the hepatic expression of OAT2/7, which may contribute to the pharmacokinetic variability of their substrates. Because both transporters were less abundant in hepatocytes than livers, a REF-based approach is recommended when scaling in vitro hepatocyte transport data to predict hepatic drug clearance and liver exposure of OAT2/7 substrates.

Points-to-consider documents: Scientific information on the evaluation of genetic polymorphisms during non-clinical studies and phase I clinical trials in the Japanese population

Pharmacotherapy shows striking individual differences in pharmacokinetics and pharmacodynamics, involving drug efficacy and adverse reactions. Recent genetic research has revealed that genetic polymorphisms are important intrinsic factors for these inter-individual differences. This pharmacogenomic information could help develop safer and more effective precision pharmacotherapies and thus, regulatory guidance/guidelines were developed in this area, especially in the EU and US. The Project for the Promotion of Progressive Medicine, Medical Devices, and Regenerative Medicine by the Ministry of Health, Labour and Welfare, performed by Tohoku University, reported scientific information on the evaluation of genetic polymorphisms, mainly on drug metabolizing enzymes and transporters, during non-clinical studies and phase I clinical trials in Japanese subjects/patients. We anticipate that this paper will be helpful in drug development for the regulatory usage of pharmacogenomic information, most notably pharmacokinetics.

Prediction and Quantification of Hepatic Transporter-Mediated Uptake of Pitavastatin Utilizing a Combination of the Relative Activity Factor Approach and Mechanistic Modeling

Quantification of the fraction transported (f_t) by a particular transporter will facilitate more robust estimations of transporter interactions. Using pitavastatin as a model uptake transporter substrate, we investigated the utility of the relative activity factor (RAF) approach and mechanistic modeling to estimate f_t in hepatocytes. The transporters evaluated were organic anion-transporting polypeptides OATP1B1 and OATP1B3 and sodium-taurocholate cotransporting polypeptide. Transporter-expressing human embryonic kidney 293 cells and human hepatocytes were used for determining RAF values, which were then incorporated into the mechanistic model to simulate hepatocyte uptake of pitavastatin over time. There was excellent agreement between simulated and observed hepatocyte uptake of pitavastatin, indicating the suitability of this approach for translation of uptake from individual transporter-expressing cells to more holistic in vitro models. Subsequently, f_t values were determined. The largest contributor to hepatocyte uptake of pitavastatin was OATP1B1, which correlates with what is known about the in vivo disposition of pitavastatin. The f_t values were then used for evaluating in vitro-in vivo correlations of hepatic uptake inhibition with OATP inhibitors rifampicin and cyclosporine. Predictions were compared with previously reported plasma exposure changes of pitavastatin with these inhibitors. Although hepatic uptake inhibition of pitavastatin was 2-3-fold underpredicted, incorporation of scaling factors (SFs) into RAF values significantly improved the predictive ability. We propose that calibration of hepatocytes with standard transporter substrates and inhibitors would allow for determination of system-specific SFs, which could subsequently be used for refining predictions of clinical DDI potential for new chemical entities that undergo active hepatic uptake.

Organic Anion Transporting Polypeptides Contribute to the Disposition of Perfluoroalkyl Acids in Humans and Rats

Perfluoroalkyl sulfonates (PFSAs) such as perfluorohexane sulfonate (PFHxS) and perfluorooctane sulfonate (PFOS) have very long serum elimination half-lives in humans, and preferentially distribute to serum and liver. The enterohepatic circulation of PFHxS and PFOS likely contributes to their extended elimination half-lives. We previously demonstrated that perfluorobutane sulfonate (PFBS), PFHxS, and PFOS are transported into hepatocytes both in a sodium-dependent and a sodium-independent manner. We identified Na+/taurocholate cotransporting polypeptide (NTCP) as the responsible sodium-dependent transporter. Furthermore, we demonstrated that the human apical sodium-dependent bile salt transporter (ASBT) contributes to the intestinal reabsorption of PFOS. However, so far no sodium-independent uptake transporters for PFSAs have been identified in human hepatocytes or enterocytes. In addition, perfluoroalkyl carboxylates (PFCAs) with 8 and 9 carbons were shown to preferentially distribute to the liver of rodents; however, no rat or human liver uptake transporters are known to transport these PFCAs. Therefore, we tested whether PFBS, PFHxS, PFOS, and PFCAs with 7-10 carbons are substrates of organic anion transporting polypeptides (OATPs). We used CHO and HEK293 cells to demonstrate that human OATP1B1, OATP1B3, and OATP2B1 can transport PFBS, PFHxS, PFOS, and the 2 PFCAs (C8 and C9). In addition, we show that rat OATP1A1, OATP1A5, OATP1B2, and OATP2B1 transport all 3 PFSAs. In conclusion, our results suggest that besides NTCP and ASBT, OATPs also are capable of contributing to the enterohepatic circulation and extended human serum elimination half-lives of the tested perfluoroalkyl acids.

Modelling of atorvastatin pharmacokinetics and the identification of the effect of a BCRP polymorphism in the Japanese population

AIM

Ethnicity plays a modulating role in atorvastatin pharmacokinetics (PK), with Asian patients reported to have higher exposure compared with Caucasians. Therefore, it is difficult to safely extrapolate atorvastatin PK data and models across ethnic groups. This work aims to develop a population PK model for atorvastatin and its pharmacologically active metabolites specifically for the Japanese population. Subsequently, it aimed to identify genetic polymorphisms affecting atorvastatin PK in this population.

METHODS

Atorvastatin acid (ATA) and ortho-hydroxy-atorvastatin acid (o-OH-ATA) plasma concentrations, clinical/demographic characteristics and genotypes for 18 (3, 3, 1, 1, 7, 2 and 1 in the ABCB1, ABCG2, CYP3A4, CYP3A5, SLCO1B1, SLCO2B1 and PPARA genes, respectively) genetic polymorphisms were collected from 27 Japanese individuals (taking 10 mg atorvastatin once daily) and analysed using a population PK modelling approach.

RESULTS

The population PK model developed (one-compartment for ATA linked through metabolite formation to an additional compartment describing the disposition of o-OH-ATA) accurately described the observed data and the associated population variability. Our analysis suggested that patients carrying one variant allele for the rs2622604 polymorphism (ABCG2) show a 55% (95% confidence interval: 16-131%) increase in atorvastatin oral bioavailability relative to the value in individuals without the variant allele.

CONCLUSION

The current work reports the identification in the Japanese population of a BCRP polymorphism, not previously associated with the PK of any statin, that markedly increases ATA and o-OH-ATA exposure. The model developed may be of clinical importance to guide dosing recommendations tailored specifically for the Japanese.

Development of a rapid and sensitive multiple reaction monitoring proteomic approach for quantification of transporters in human liver tissue

With increasing knowledge on the role of hepatic transporters in drug disposition, numerous efforts have been described to quantify the expression of human hepatic transporters. However, reported quantitative proteomic approaches often require long analysis times. Additionally, greater assay sensitivity is still necessary for less abundant transporters or limited quantity of samples (e.g. hepatocytes and liver tissue). In the present study, an LC-MS/MS method for rapid and simultaneous quantification of 12 hepatic transporters (BCRP, BSEP, MATE1, MRP2, MRP3, MRP4, NTCP, OATP1B1, 1B3, 2B1, OCT1, and P-gp) was developed. Using a high LC flow rate (1.5mL/min) and fast LC gradient (4min total cycle time), the run time was markedly reduced to 4min, which was much shorter than most previously published assays. Chromatographic separation was achieved using ACE UltraCore SuperC18 50mm×2.1mm 5-μm HPLC column. In addition, greater analytical sensitivity was achieved with both high LC flow rate/fast LC gradient and post-column infusion of ethylene glycol. The on-column LLOQ for signature peptides in this method ranged from 0.194 to 0.846 femtomoles. The impact of five protein solubilizers, including extraction buffer II of ProteoExtract Native Membrane Protein Extraction Kit, 3% (w/v) sodium deoxycholate, 20% (v/v) Invitrosol, 0.2% (w/v) RapiGest SF, and 10% (w/v) formamide on total membrane protein extraction and trypsin digestion was investigated. Sodium deoxycholate was chosen because of good total membrane protein extraction and trypsin digestion efficiency, as well as no significant MS interference. Good precision (within 15% coefficient of variation) and accuracy (within ±15% bias), and inter-day trypsin digestion efficiency (within 28% coefficient of variation) was observed for quality controls. This method can quantify human hepatic transporter expression in a high-throughput manner and due to the increased sensitivity can be used to investigate the down-regulation of hepatic transporter protein (e.g., different ethnic groups and liver disease patients).

Strategies of Drug Transporter Quantitation by LC-MS: Importance of Peptide Selection and Digestion Efficiency

Huge variation of drug transporter abundance was seen in the literature, making PBPK prediction difficult when transporters play a major role. Among multiple factors such as membrane fraction, digestion, and peptide selection that contributed to such variation, peptide selection is the least discussed. Herein, a strategy was established by using a small amount of purified protein standard to select a peptide with near 100% digestion efficiency for quantitation of a transporter protein MDR1. The impact of native membrane protein's tertiary structure on the digestion efficiency of surrogate peptides of MDR1 was investigated. Peptides in more solvent accessible regions are found to be digested much more efficiently than those in large stretches of helical structures. The concentration of peptide EALDESIPPVSFWR(EAL) in the most solvent accessible linker region of MDR1 was found closest to the true protein concentration. When using EAL for MDR1 quantitation, the abundance is over 10 times higher than previously reported, indicating the importance of peptide selection for transporter quantitation. In addition, this study also proposes a screening strategy to select peptides appropriate for relative quantitation for in vitro-in vivo extrapolation in the absence of any protein standard.

Glucuronidation: driving factors and their impact on glucuronide disposition

Glucuronidation is a well-recognized phase II metabolic pathway for a variety of chemicals including drugs and endogenous substances. Although it is usually the secondary metabolic pathway for a compound preceded by phase I hydroxylation, glucuronidation alone could serve as the dominant metabolic pathway for many compounds, including some with high aqueous solubility. Glucuronidation involves the metabolism of parent compound by UDP-glucuronosyltransferases (UGTs) into hydrophilic and negatively charged glucuronides that cannot exit the cell without the aid of efflux transporters. Therefore, elimination of parent compound via glucuronidation in a metabolic active cell is controlled by two driving forces: the formation of glucuronides by UGT enzymes and the (polarized) excretion of these glucuronides by efflux transporters located on the cell surfaces in various drug disposition organs. Contrary to the common assumption that the glucuronides reaching the systemic circulation were destined for urinary excretion, recent evidences suggest that hepatocytes are capable of highly efficient biliary clearance of the gut-generated glucuronides. Furthermore, the biliary- and enteric-eliminated glucuronides participate into recycling schemes involving intestinal microbes, which often prolong their local and systemic exposure, albeit at low systemic concentrations. Taken together, these recent research advances indicate that although UGT determines the rate and extent of glucuronide generation, the efflux and uptake transporters determine the distribution of these glucuronides into blood and then to various organs for elimination. Recycling schemes impact the apparent plasma half-life of parent compounds and their glucuronides that reach intestinal lumen, in addition to prolonging their gut and colon exposure.

Rosuvastatin Pharmacokinetics in Asian and White Subjects Wild Type for Both OATP1B1 and BCRP Under Control and Inhibited Conditions

The Food and Drug Administration recommends rosuvastatin dosage reductions in Asian patients because pharmacokinetic studies have demonstrated an approximate 2-fold increase in median exposure to rosuvastatin in Asian subjects compared with Caucasian controls. Yet, no explanation for this ethnic difference has been confirmed. Here we show that rosuvastatin exposure in Asians and Whites does not differ significantly when all subjects are wild-type carriers for both solute carrier organic anion transporter 1B1 *1a and ATP-binding cassette subfamily G member 2 c.421 transporters in a 2-arm, randomized, cross-over rosuvastatin pharmacokinetics study in healthy white and Asian volunteers. For single rosuvastatin doses, AUC_0-48 were 92.5 (±36.2) and 83.5 (±32.2) ng/mL × h and C_max were 10.0 (±4.1) and 7.6 (±3.0) ng/mL for Asians and Whites, respectively. When transporters were inhibited by intravenous rifampin, rosuvastatin AUC_0-48 and C_max also showed no ethnic differences. Our study suggests that both SLCO1B1 and ABCG2 polymorphisms are better predictors of rosuvastatin exposure than ethnicity alone and could be considered in precision medicine dosing of rosuvastatin.

Interindividual and interethnic variability in drug disposition: polymorphisms in organic anion transporting polypeptide 1B1 (OATP1B1; SLCO1B1)

OATP1B1 (SLCO1B1) is predominantly expressed at the basolateral membrane of hepatocytes and is critically important for the hepatic uptake and clearance of numerous drug substrates and endogenous compounds. In general, the organic anion transporting polypeptides (OATP; SLCO) represent a superfamily of uptake transporters that mediate the sodium-independent transport of a diverse range of amphipathic organic compounds including bile salts, steroid conjugates, thyroid hormones, anionic peptides, numerous drugs and other xenobiotic substances. OATP1B1 is highly polymorphic and a number of relevant and ethnically dependent polymorphisms have been identified and functionally characterized. In particular, the SLCO1B1 521T>C and 388A>G polymorphisms are commonly occurring variants in ethnically diverse populations and numerous in vitro and clinical studies have evaluated the consequences of these variants to interindividual differences in drug disposition and response. OATP1B1 is particularly important for the disposition of HMG-CoA reductase inhibitors, or statins, as it is known to efficiently transport most statins to their site of action within hepatocytes. Many studies have focused on the consequences of OATP1B1 variants to statin disposition in vitro and in vivo and would suggest that genetic variability in SLCO1B1 has important implications for statin pharmacokinetics, risk for statin-induced myopathy, and modulation of statin treatment response. This review describes what is currently known regarding SLCO1B1 genotype, OATP1B1 protein expression and interindividual and interethnic consequences to drug disposition, with particular focus on statin pharmacokinetics and implications for drug response and toxicity.

Pharmacogenetics of Lipid-Lowering Agents: Precision or Indecision Medicine?

Lipid-lowering medications, particularly statins, have been a popular target for pharmacogenetic studies. A handful of genes have shown promise for predicting response to therapy from the perspective of lipid lowering, as well as myopathy. A number of genes have been implicated and have biological plausibility based on their involvement with the pharmacokinetics or pharmacodynamics of statins or other lipid-lowering medications. The level of confidence and replication of these findings varies, although several associations are likely true. Novel classes of lipid-lowering therapy have opened up new possibilities in the treatment of severe inherited forms of dyslipidemia, making the identification of such mutations an important pharmacogenetic predictor of failure of standard therapy, with potential response to novel therapy. Advances in next-generation sequencing technology bring the application of pharmacogenetics even closer to routine clinical practice.

Abundance of Hepatic Transporters in Caucasians: A Meta-Analysis

This study aimed to derive quantitative abundance values for key hepatic transporters suitable for in vitro–in vivo extrapolation within a physiologically based pharmacokinetic modeling framework. A meta-analysis was performed whereby data on abundance measurements, sample preparation methods, and donor demography were collated from the literature. To define values for a healthy Caucasian population, a subdatabase was created whereby exclusion criteria were applied to remove samples from non-Caucasian individuals, those with underlying disease, or those with subcellular fractions other than crude membrane. Where a clinically relevant active genotype was known, only samples from individuals with an extensive transporter phenotype were included. Authors were contacted directly when additional information was required. After removing duplicated samples, the weighted mean, geometric mean, standard deviation, coefficient of variation, and between-study homogeneity of transporter abundances were determined. From the complete database containing 24 transporters, suitable abundance data were available for 11 hepatic transporters from nine studies after exclusion criteria were applied. Organic anion transporting polypeptides OATP1B1 and OATP1B3 showed the highest population abundance in healthy adult Caucasians. For several transporters, the variability in abundance was reduced significantly once the exclusion criteria were applied. The highest variability was observed for OATP1B3 > OATP1B1 > multidrug resistance protein 2 > multidrug resistance gene 1. No relationship was found between transporter expression and donor age. To our knowledge, this study provides the first in-depth analysis of current quantitative abundance data for a wide range of hepatic transporters, with the aim of using these data for in vitro–in vivo extrapolation, and highlights the significance of investigating the background of tissue(s) used in quantitative transporter proteomic studies. Similar studies are now warranted for other ethnicities.

Quantification of Four Efflux Drug Transporters in Liver and Kidney Across Species Using Targeted Quantitative Proteomics by Isotope Dilution NanoLC-MS/MS

PURPOSE

The expression levels of several efflux drug transporters in the liver and kidney were evaluated across species to address potential roles of the transporters in species dependent excretion of drugs and their metabolites.

METHODS

Four efflux transporters, namely MDR1/P-gp, BCRP/Bcrp, MRP2/Mrp2 and MRP3/Mrp3 in liver and kidney in three preclinical species and humans were quantified using targeted quantitative proteomics by isotope dilution nanoLC-MS/MS.

RESULTS

In liver, the level of P-gp was highest in monkey and lowest in rat. The concentration of BCRP/Bcrp was highest in dog followed by monkey. MRP2/Mrp2 level was highest in monkey and rat, whereas MRP3/Mrp3 levels were similar in human, monkey and dog. In the kidney, the concentrations of MDR1/P-gp in human and monkey were roughly 2 to 3-fold higher than in rat and dog. In rat, BCRP/Bcrp concentrations were substantially higher than in any of the other species. MRP2/Mrp2 concentrations were similar across species, whereas expression of MRP3/Mrp3 was highest in rat.

CONCLUSION

Overall, the results indicated that the pattern of hepatic and renal expression of the transporters was quite species dependent. This information should be helpful in the estimation of transport mediated drug and metabolites excretion in liver and kidney across species.

Quantification of Flavin-containing Monooxygenases 1, 3, and 5 in Human Liver Microsomes by UPLC-MRM-Based Targeted Quantitative Proteomics and Its Application to the Study of Ontogeny

Flavin-containing monooxygenases (FMOs) have a significant role in the metabolism of small molecule pharmaceuticals. Among the five human FMOs, FMO1, FMO3, and FMO5 are the most relevant to hepatic drug metabolism. Although age-dependent hepatic protein expression, based on immunoquantification, has been reported previously for FMO1 and FMO3, there is very little information on hepatic FMO5 protein expression. To overcome the limitations of immunoquantification, an ultra-performance liquid chromatography (UPLC)-multiple reaction monitoring (MRM)-based targeted quantitative proteomic method was developed and optimized for the quantification of FMO1, FMO3, and FMO5 in human liver microsomes (HLM). A post-in silico product ion screening process was incorporated to verify LC-MRM detection of potential signature peptides before their synthesis. The developed method was validated by correlating marker substrate activity and protein expression in a panel of adult individual donor HLM (age 39-67 years). The mean (range) protein expression of FMO3 and FMO5 was 46 (26-65) pmol/mg HLM protein and 27 (11.5-49) pmol/mg HLM protein, respectively. To demonstrate quantification of FMO1, a panel of fetal individual donor HLM (gestational age 14-20 weeks) was analyzed. The mean (range) FMO1 protein expression was 7.0 (4.9-9.7) pmol/mg HLM protein. Furthermore, the ontogenetic protein expression of FMO5 was evaluated in fetal, pediatric, and adult HLM. The quantification of FMO proteins also was compared using two different calibration standards, recombinant proteins versus synthetic signature peptides, to assess the ratio between holoprotein versus total protein. In conclusion, a UPLC-MRM-based targeted quantitative proteomic method has been developed for the quantification of FMO enzymes in HLM.

In Vitro-In Vivo Extrapolation Scaling Factors for Intestinal P-Glycoprotein and Breast Cancer Resistance Protein: Part I: A Cross-Laboratory Comparison of Transporter-Protein Abundances and Relative Expression Factors in Human Intestine and Caco-2 Cells

Over the last 5 years the quantification of transporter-protein absolute abundances has dramatically increased in parallel to the expanded use of in vitro-in vivo extrapolation (IVIVE) and physiologically based pharmacokinetics (PBPK)-linked models, for decision-making in pharmaceutical company drug development pipelines and regulatory submissions. Although several research groups have developed laboratory-specific proteomic workflows, it is unclear if the large range of reported variability is founded on true interindividual variability or experimental variability resulting from sample preparation or the proteomic methodology used. To assess the potential for methodological bias on end-point abundance quantification, two independent laboratories, the University of Manchester (UoM) and Bertin Pharma (BPh), employing different proteomic workflows, quantified the absolute abundances of Na/K-ATPase, P-gp, and breast cancer resistance protein (BCRP) in the same set of biologic samples from human intestinal and Caco-2 cell membranes. Across all samples, P-gp abundances were significantly correlated (P = 0.04, Rs = 0.72) with a 2.4-fold higher abundance (P = 0.001) generated at UoM compared with BPh. There was a systematically higher BCRP abundance in Caco-2 cell samples quantified by BPh compared with UoM, but not in human intestinal samples. Consequently, a similar intestinal relative expression factor (REF), derived from distal jejunum and Caco-2 monolayer samples, between laboratories was found for P-gp. However, a 2-fold higher intestinal REF was generated by UoM (2.22) versus BPh (1.11). We demonstrate that differences in absolute protein abundance are evident between laboratories and they probably result from laboratory-specific methodologies relating to peptide choice.