Dealing With Variable Drug Exposure Due to Variable Hepatic Metabolism: A Proof-of-Concept Application of Liquid Biopsy in Renal Impairment.

Precision dosing strategies require accounting for between-patient variability in pharmacokinetics (PK), affecting drug exposure, and in pharmacodynamics (PD), affecting response achieved at the same drug concentration at the site of action. Although liquid biopsy for assessing different levels of molecular drug targets has yet to be established, individual characterization of drug elimination pathways using liquid biopsy has recently been demonstrated. The feasibility of applying this approach in conjunction with modeling tools to guide individual dosing remains unexplored. In this study, we aimed to individualize physiologically-based pharmacokinetic (PBPK) models based on liquid biopsy measurements in plasma from 25 donors with different grades of renal function who were previously administered oral midazolam as part of a microdose cocktail. Virtual twin models were constructed based on demographics, renal function, and hepatic expression of relevant pharmacokinetic pathways projected from liquid biopsy output. Simulated exposure (AUC) to midazolam was in agreement with observed data (AFE = 1.38, AAFE = 1.78). Simulated AUC variability with three dosing approaches indicated higher variability with uniform dosing (14-fold) and stratified dosing (13-fold) compared with individualized dosing informed by liquid biopsy (fivefold). Further, exosome screening revealed mRNA expression of 532 targets relevant to drug metabolism and disposition (169 enzymes and 361 transporters). Data related to these targets can be used to further individualize PBPK models for pathways relevant to PK of other drugs. This study provides additional verification of liquid biopsy-informed PBPK modeling approaches, necessary to advance strategies that seek to achieve precise dosing from the start of treatment.

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.

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

Introduction: Alterations in expression and activity of human receptor tyrosine kinases (RTKs) are associated with cancer progression and in response to therapeutic intervention. Methods: Thus, protein abundance of 21 RTKs was assessed in 15 healthy and 18 cancerous liver samples [2 primary and 16 colorectal cancer liver metastasis (CRLM)] matched with non-tumorous (histologically normal) tissue, by a validated QconCAT-based targeted proteomic approach. Results: It was demonstrated, for the first time, that the abundance of EGFR, INSR, VGFR3 and AXL, is lower in tumours relative to livers from healthy individuals whilst the opposite is true for IGF1R. EPHA2 was upregulated in tumour compared with histologically normal tissue surrounding it. PGFRB levels were higher in tumours relative to both histologically normal tissue surrounding tumour and tissues taken from healthy individuals. The abundances of VGFR1/2, PGFRA, KIT, CSF1R, FLT3, FGFR1/3, ERBB2, NTRK2, TIE2, RET, and MET were, however, comparable in all samples. Statistically significant, but moderate correlations were observed (Rs > 0.50, p < 0.05) for EGFR with INSR and KIT. FGFR2 correlated with PGFRA and VGFR1 with NTRK2 in healthy livers. In non-tumorous (histologically normal) tissues from cancer patients, there were correlations between TIE2 and FGFR1, EPHA2 and VGFR3, FGFR3 and PGFRA (p < 0.05). EGFR correlated with INSR, ERBB2, KIT and EGFR, and KIT with AXL and FGFR2. In tumours, CSF1R correlated with AXL, EPHA2 with PGFRA, and NTRK2 with PGFRB and AXL. Sex, liver lobe and body mass index of donors had no impact on the abundance of RTKs, although donor age showed some correlations. RET was the most abundant of these kinases in non-tumorous tissues (~35%), while PGFRB was the most abundant RTK in tumours (~47%). Several correlations were also observed between the abundance of RTKs and proteins relevant to drug pharmacokinetics (enzymes and transporters). Discussion: DiscussionThis study quantified perturbation to the abundance of several RTKs in cancer and the value generated in this study can be used as input to systems biology models defining liver cancer metastases and biomarkers of its progression.

Seeing what is behind the smokescreen: A systematic review of methodological aspects of smoking interaction studies over the last three decades and implications for future clinical trials.

Smoking drug interaction studies represent a common approach for the clinical investigation of CYP1A2 induction. Despite this important role, they remain an "orphan topic" in the existing regulatory framework of drug interaction studies, and important methodological aspects remain unaddressed. The University of Washington Drug Interaction Database (DIDB) was used to systematically review the published literature on dedicated smoking pharmacokinetic interaction studies in healthy subjects (1990 to 2021, inclusive). Various methodological aspects of identified studies were reviewed. A total of 51 studies met all inclusion criteria and were included in the analysis. Our review revealed that methods applied in smoking interaction studies are heterogeneous and often fall short of established methodological standards of other interaction trials. Methodological deficiencies included incomplete description of study populations, poor definition and lack of objective confirmation of smoker and nonsmoker characteristics, under-representation of female subjects, small sample sizes, frequent lack of statistical sample size planning, frequent lack of use of existing markers of nicotine exposure and CYP1A2 activity measurements, and frequent lack of control of extrinsic CYP1A2 inducing or inhibiting factors. The frequent quality issues in the assessment and reporting of smoking interaction trials identified in this review call for a concerted effort in this area, if the results of these studies are meant to be followed by actionable decisions on dose optimization, when needed, for the effects of smoking on CYP1A2 victim drugs in smokers.

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

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

Delineating gene-environment effects using virtual twins of patients treated with clozapine.

Studies that focus on individual covariates, while ignoring their interactions, may not be adequate for model-informed precision dosing (MIPD) in any given patient. Genetic variations that influence protein synthesis should be studied in conjunction with environmental covariates, such as cigarette smoking. The aim of this study was to build virtual twins (VTs) of real patients receiving clozapine with interacting covariates related to genetics and environment and to delineate the impact of interacting covariates on predicted clozapine plasma concentrations. Clozapine-treated patients with schizophrenia (N = 42) with observed clozapine plasma concentrations, demographic, environmental, and genotype data were used to construct VTs in Simcyp. The effect of increased covariate virtualization was assessed by performing simulations under three conditions: "low" (demographic), "medium" (demographic and environmental interaction), and "high" (demographic and environmental/genotype interaction) covariate virtualization. Increasing covariate virtualization with interaction improved the coefficient of variation (R2 ) from 0.07 in the low model to 0.391 and 0.368 in the medium and high models, respectively. Whereas R2 was similar between the medium and high models, the high covariate virtualization model had improved accuracy, with systematic bias of predicted clozapine plasma concentration improving from -138.48 ng/ml to -74.65 ng/ml. A high level of covariate virtualization (demographic, environmental, and genotype) may be required for MIPD using VTs.

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

BACKGROND AND OBJECTIVE: Crohn's disease (CD) is a chronic inflammatory bowel disease that affects a wide age range. Hence, CD patients receive a variety of drugs over their life beyond those used for CD itself. The changes to the integrity of the intestine and its drug metabolising enzymes and transporters (DMETs) can alter the oral bioavailability of drugs. However, there are other changes in systems parameters determining the fate of drugs in CD, and understanding these is essential for dose adjustment in patients with CD. METHODS: The current analysis gathered all the available clinical data on the kinetics of drugs in CD (by March 2021), focusing on orally administered small molecule drugs. A meta-analysis of the systems parameters affecting oral drug pharmacokinetics was conducted. The systems information gathered on intestine, liver and blood proteins and other physiological parameters was incorporated into a physiologically based pharmacokinetic (PBPK) platform to create a virtual population of CD patients, with a view for guiding dose adjustment in the absence of clinical data in CD. RESULTS: There were no uniform trends in the reported changes in reported oral bioavailability. The nature of the drug as well as the formulation affected the direction and magnitude of variation in kinetics in CD patients relative to healthy volunteers. Even for the same drug, the reported changes in exposure varied, possibly due to a lack of distinction between the activity states of CD. The highest alteration was seen with S-verapamil and midazolam, 8.7- and 5.3-fold greater exposure, respectively, in active CD patients relative to healthy volunteers. Only one report was available on liver DMETs in CD, and indicated reduced CYP3A4 activity. In a number of reports, mRNA expression of DMETs in the ileum and colon of CD patients was measured, focussing on P-glycoprotein (p-gp) transporter and CYP3A4 enzyme, and showed contradictory results. No data were available on protein expression in duodenum and jejunum despite their dominant role in oral drug absorption. CONCLUSION: There are currently inadequate dedicated clinical or quantitative proteomic studies in CD to enable predictive PBPK models with high confidence and adequate verification. The PBPK models for CD with the available systems parameters were able to capture the major physiological influencers and the gaps to be filled by future research. Quantification of DMETs in the intestine and the liver in CD is warranted, alongside well-defined clinical drug disposition studies with a number of index drugs as biomarkers of changes in DMETs in these patients, to avoid large-scale dedicated studies for every drug to determine the effects of disease on the drug's metabolism and disposition and the consequential safety and therapeutic concerns.

Quantitative Assessment of the Impact of Crohn's Disease on Protein Abundance of Human Intestinal Drug-Metabolising Enzymes and Transporters.

Crohn's disease affects the mucosal layer of the intestine, predominantly ileum and colon segments, with the potential to affect the expression of intestinal enzymes and transporters, and consequently, oral drug bioavailability. We carried out a quantitative proteomic analysis of inflamed and non-inflamed ileum and colon tissues from Crohn's disease patients and healthy donors. Homogenates from samples in each group were pooled and protein abundance determined by liquid chromatography-mass spectrometry (LC-MS). In inflamed Crohn's ileum, CYP3A4, CYP20A1, CYP51A1, ADH1B, ALPI, FOM1, SULT1A2, SULT1B1 and ABCB7 showed ≥10-fold reduction in abundance compared with healthy baseline. By contrast, only MGST1 showed ≥10 fold reduction in inflamed colon. Ileal UGT1A1, MGST1, MGST2, and MAOA levels increased by ≥2 fold in Crohn's patients, while only ALPI showed ≥2 fold increase in the colon. Counter-intuitively, non-inflamed ileum had a higher magnitude of fold change than inflamed tissue when compared with healthy tissue. Marked but non-uniform alterations were observed in the expression of various enzymes and transporters in ileum and colon compared with healthy samples. Modelling will allow improved understanding of the variable effects of Crohn's disease on bioavailability of orally administered drugs.

Effect of Chronic Kidney Disease on the Renal Secretion via Organic Anion Transporters 1/3: Implications for Physiologically-Based Pharmacokinetic Modeling and Dose Adjustment.

There is growing evidence that active tubular secretory clearance (CLs ) may not decline proportionally with the glomerular filtration rate (GFR) in chronic kidney disease (CKD), leading to the overestimation of renal clearance (CLr ) when using solely GFR to approximate disease effect on renal elimination. The clinical pharmacokinetic data of 33 renally secreted OAT1/3 substrates were collated to investigate the impact of mild, moderate, and severe CKD on CLr , tubular secretion and protein binding (fu,p ). The fu,p of the collated substrates ranged from 0.0026 to 1.0 in healthy populations; observed CKD-related increase in the fu,p (up to 2.7-fold) of 8 highly bound substrates (fu,p ≤ 0.2) was accounted for in the analysis. Use of prediction equation based on disease-related changes in albumin resulted in underprediction of the CKD-related increase in fu,p of highly bound substrates, highlighting the necessity to measure protein binding in severe CKD. The critical analysis of clinical data for 33 OAT1/3 probes established that decrease in OAT1/3 activity proportional to the changes in GFR was insufficient to recapitulate effects of severe CKD on unbound tubular secretion clearance. OAT1/3-mediated CLs was estimated to decline by an additional 50% relative to the GFR decline in severe CKD, whereas change in active secretion in mild and moderate CKD was proportional to GFR. Consideration of this additional 50% decline in OAT1/3-mediated CLs is recommended for physiologically-based pharmacokinetic models and dose adjustment of OAT1/3 substrates in severe CKD, especially for substrates with high contribution of the active secretion to CLr .

Proteomic quantification of perturbation to pharmacokinetic target proteins in liver disease.

Model-based assessment of drug pharmacokinetics in liver disease requires quantification of abundance and disease-related changes in hepatic enzymes and transporters. This study aimed to assess performance of three label-free methods [high N (HiN), intensity-based absolute quantification (iBAQ) and total protein approach (TPA)] against QconCAT-based targeted data in healthy and diseased (cancer and cirrhosis) liver tissue. Measurements were compared across methods and disease-to-control ratios provided a 'disease perturbation factor' (DPF) for each protein. Mean label-free measurements of targets correlated well (Pearson's coefficient, r = 0.91-0.98 p < 0.001) and with targeted data (r = 0.65-0.95, p < 0.001). Concordance with targeted data was generally moderate (Lin's concordance coefficient, ρc = 0.46-0.92), depending on methodology. Moderate precision and accuracy were observed for label-free methods (average fold error, AFE = 1.44-1.68; absolute average fold error, AAFE = 2.44-3.23). The DPF reconciled the data and indicated downregulated expression in cancer and cirrhosis, consistent with an inflammatory effect. HiN estimated perturbation consistently with targeted data (AFEHiN = 1.07, AAFEHiN = 1.57), whereas iBAQ overestimated (AFEiBAQ = 0.81, AAFEiBAQ = 1.67) and TPA underestimated (AFETPA = 1.37, AAFETPA = 1.65) disease effect. Progression from mild to severe cirrhosis was consistent with progressive decline in expression, reproduced by HiN but overestimated by iBAQ and underestimated by TPA (AFEHiN = 0.98, AFEiBAQ = 0.60, AFETPA = 1.24). DPF data confirmed non-uniform disease effect on drug-elimination pathways and progressive impact of disease severity. SIGNIFICANCE: This study demonstrated good correlation and moderate concordance between intensity-based label-free proteomic methods (HiN, iBAQ and TPA) and targeted data. Label-free measurements tended to overestimate abundance, but differences were reconciled using a disease perturbation factor (DPF) for each protein. With targeted data as a reference, HiN defined disease perturbation and the impact of disease progression consistently, indicating that the use of 'razor' peptides for quantification against an exogenous standard provides biologically sensible quantitative fingerprints of disease. Disease-driven perturbations in expression relative to healthy baseline are incorporated into drug kinetic models used to predict drug exposure in disease populations where clinical studies may not be feasible.

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

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

Liquid Biopsy for Patient Characterization in Cardiovascular Disease: Verification against Markers of Cytochrome P450 and P-Glycoprotein Activities.

Precision dosing strategies require accounting for between-patient variability in pharmacokinetics together with subsequent pharmacodynamic differences. Liquid biopsy is a valuable new approach to diagnose disease prior to the appearance of clinical signs and symptoms, potentially circumventing invasive tissue biopsies. However, the possibility of quantitative grading of biomarkers, as opposed to simply confirming their presence or absence, is relatively new. In this study, we aimed to verify expression measurements of cytochrome P450 (CYP) enzymes and the transporter P-glycoprotein (P-gp) in liquid biopsy against genotype and activity phenotype (assessed by the Geneva cocktail approach) in 30 acutely ill patients with cardiovascular disease in a hospital setting. After accounting for exosomal shedding, expression in liquid biopsy correlated with activity phenotype for CYP1A2, CYP2B6, CYP2C9, CYP3A, and P-gp (r = 0.44-0.70, P ≤ 0.05). Although genotype offered a degree of stratification, large variability (coefficient of variation (CV)) in activity (up to 157%) and expression in liquid biopsy (up to 117%) was observed within each genotype, indicating a mismatch between genotype and phenotype. Further, exosome screening revealed expression of 497 targets relevant to drug metabolism and disposition (159 enzymes and 336 transporters), as well as 20 molecular drug targets. Although there were no functional data available to correlate against these large-scale measurements, assessment of disease perturbation from healthy baseline was possible. Verification of liquid biopsy against activity phenotype is important to further individualize modeling approaches that aspire to achieve precision dosing from the start of drug treatment without the need for multiple rounds of dose optimization.

A family of QconCATs (Quantification conCATemers) for the quantification of human pharmacological target proteins.

We have developed a family of QconCAT standards for the absolute quantification of pharmacological target proteins in a variety of human tissues. The QconCATs consist of concatenated proteotypic peptides, are designed in silico, and expressed in E. coli in media enriched with [13C6] arginine and [13C6] lysine to generate stable isotope-labeled multiplexed absolute quantification standards. The so-called MetCAT (used to quantify cytochrome P450 (CYP) and glucuronosyltransferase (UGT) enzymes), the liver TransCAT (used to quantify plasma-membrane drug transporters) and the brain TransCAT (used to quantify transporters expressed in the blood-brain barrier) were previously reported. We now report new QconCATs for the quantification of non-UGT non-CYP drug metabolizing enzymes (NuncCAT) and receptor tyrosine kinases (KinCAT). We have also redesigned the liver TransCAT, replacing problematic peptides and the N-terminal tag, for better characterization and expression. All these QconCATs showed high purity, high labelling efficiency with stable 13C isotope (>95%), and high sequence coverage (>88%). They represent a close-knit family of standards for quantifying pharmacokinetic targets, together with a more distant cousin, the KinCAT, used to quantify pharmacodynamic targets. SIGNIFICANCE: Multiplexed determination of absolute protein abundances using quantitative conCATemers (QconCATs) has already been successfully demonstrated in different human tissues. We have previously reported two QconCATs; MetCAT and TransCAT, for the quantification of key enzymes (cytochrome P450 enzymes (CYP) and glucuronosyltransferases (UGT)) and drug transporters. To build on these reports, application of the QconCAT methodology for the determination of non-UGT non-CYP enzymes and receptor tyrosine kinases (RTKs) in human tissue is reported here. This report focuses on development and characterization of two QconCAT constructs for the quantification of 24 enzymes and 21 RTKs. We demonstrate that the developed QconCATs have high purity, high incorporation efficiency and low peptide miscleavage upon proteolysis. Application of these QconCATs for reliable quantification of target proteins was achieved in human liver.

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

AIMS: This study aims to quantify drug-metabolising enzymes, transporters, receptor tyrosine kinases (RTKs) and protein markers (involved in pathways affected in cancer) in pooled healthy, histologically normal and matched cancerous liver microsomes from colorectal cancer liver metastasis (CRLM) patients. METHODS: Microsomal fractionation was performed and pooled microsomes were prepared. Global and accurate mass and retention time liquid chromatography-mass spectrometry proteomics were used to quantify proteins. A QconCAT (KinCAT) for the quantification of RTKs was designed and applied for the first time. Physiologically based pharmacokinetic (PBPK) simulations were performed to assess the contribution of altered abundance of drug-metabolising enzymes and transporters to changes in pharmacokinetics. RESULTS: Most CYPs and UGTs were downregulated in histologically normal relative to healthy samples, and were further reduced in cancer samples (up to 54-fold). The transporters, MRP2/3, OAT2/7 and OATP2B1/1B3/1B1 were downregulated in CRLM. Application of abundance data in PBPK models for substrates with different attributes indicated substantially lower (up to 13-fold) drug clearance when using cancer-specific instead of default parameters in cancer population. Liver function markers were downregulated, while inflammation proteins were upregulated (by up to 76-fold) in cancer samples. Various pharmacodynamics markers (e.g. RTKs) were altered in CRLM. Using global proteomics, we examined proteins in pathways relevant to cancer (such as metastasis and desmoplasia), including caveolins and collagen chains, and confirmed general over-expression of such pathways. CONCLUSION: This study highlights impaired drug metabolism, perturbed drug transport and altered abundance of cancer markers in CRLM, demonstrating the importance of population-specific abundance data in PBPK models for cancer.

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

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

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.

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

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

Liquid Biopsy Enables Quantification of the Abundance and Interindividual Variability of Hepatic Enzymes and Transporters.

Variability in individual capacity for hepatic elimination of therapeutic drugs is well recognized and is associated with variable expression and activity of liver enzymes and transporters. Although genotyping offers some degree of stratification, there is often large variability within the same genotype. Direct measurement of protein expression is impractical due to limited access to tissue biopsies. Hence, determination of variability in hepatic drug metabolism and disposition using liquid biopsy (blood samples) is an attractive proposition during drug development and in clinical practice. This study used a multi-"omic" strategy to establish a liquid biopsy technology intended to assess hepatic capacity for metabolism and disposition in individual patients. Plasma exosomal analysis (n = 29) revealed expression of 533 pharmacologically relevant genes at the RNA level, with 147 genes showing evidence of expression at the protein level in matching liver tissue. Correction of exosomal RNA expression using a novel shedding factor improved correlation against liver protein expression for 97 liver-enriched genes. Strong correlation was demonstrated for 12 key drug-metabolizing enzymes and 4 drug transporters. The developed test allowed reliable patient stratification, and in silico trials demonstrated utility in adjusting drug dose to achieve similar drug exposure between patients with variable hepatic elimination. Accordingly, this approach can be applied in characterization of volunteers prior to enrollment in clinical trials and for patient stratification in clinical practice to achieve more precise individual dosing.

Quantification of Proteins Involved in Intestinal Epithelial Handling of Xenobiotics.

The intestinal epithelium represents a natural barrier against harmful xenobiotics, while facilitating the uptake of nutrients and other substances. Understanding the interaction of chemicals with constituents of the intestinal epithelium and their fate in the body requires quantitative measurement of relevant proteins in in vitro systems and intestinal epithelium. Recent studies have highlighted the mismatch between messenger RNA (mRNA) and protein abundance for several drug-metabolizing enzymes and transporters in the highly dynamic environment of the intestinal epithelium; mRNA abundances cannot therefore be used as a proxy for protein abundances in the gut, necessitating direct measurements. The objective was to determine the expression of a wide range proteins pertinent to metabolism and disposition of chemicals and nutrients in the intestinal epithelium. Ileum and jejunum biopsy specimens were obtained from 16 patients undergoing gastrointestinal elective surgery. Mucosal fractions were prepared and analyzed using targeted and global proteomic approaches. A total of 29 enzymes, 32 transporters, 6 tight junction proteins, 2 adhesion proteins, 1 alkaline phosphatase, 1 thioredoxin, 5 markers, and 1 regulatory protein were quantified-60 for the first time. The global proteomic method identified a further 5,222 proteins, which are retained as an open database for interested parties to explore. This study significantly expands our knowledge of a wide array of proteins important for xenobiotic handling in the intestinal epithelium. Quantitative systems biology models will benefit from the novel systems data generated in the present study and the translation path offered for in vitro to in vivo translation.