Proteomic quantification of perturbation to pharmacokinetic target proteins in liver disease

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.

Pharmacogene expression during progression of metabolic dysfunction-associated steatotic liver disease: Studies on mRNA and protein levels and their relevance to drug treatment

Metabolic dysfunction-associated steatotic liver disease (MASLD) is common worldwide. Genes and proteins contributing to drug disposition may show altered expression as MASLD progresses. To assess this further, we undertook transcriptomic and proteomic analysis of 137 pharmacogenes in liver biopsies from a large MASLD cohort. We performed sequencing on RNA from 216 liver biopsies (206 MASLD and 10 controls). Untargeted mass spectrometry proteomics was performed on a 103 biopsy subgroup. Selected RNA sequencing signals were replicated with an additional 187 biopsies. Comparison of advanced MASLD (fibrosis score 3/4) with milder disease (fibrosis score 0-2) by RNA sequencing showed significant alterations in expression of certain phase I, phase II and ABC transporters. For cytochromes P450, CYP2C19 showed the most significant decreased expression (30 % of that in mild disease) but significant decreased expression of other CYPs (including CYP2C8 and CYP2E1) also occurred. CYP2C19 also showed a significant decrease comparing the inflammatory form of MASLD (MASH) with non-MASH biopsies. Findings for CYP2C19 were confirmed in the replication cohort. Proteomics on the original discovery cohort confirmed decreased levels of several CYPs as MASLD advanced but this decrease was greatest for CYP2C19 where levels fell to 40 % control. This decrease may result in decreased CYP2C19 activity that could be problematic for prescription of drugs activated or metabolized by CYP2C19 as MASLD advances. More limited decreases for other P450s suggest fewer issues with non-CYP2C19 drug substrates. Negative correlations at RNA level between CYP2C19 and several cytokine genes provided initial insights into the mechanism underlying decreased expression.

In Vitro to In Vivo Scalars for Drug Clearance in Nonalcoholic Fatty Liver and Steatohepatitis

In vitro-in vivo extrapolation (IVIVE) allows prediction of clinical outcomes across populations from in vitro data using specific scalars tailored to the biologic characteristics of each population. This study experimentally determined scalars for patients with varying degrees of nonalcoholic fatty liver disease (NAFLD), ranging from fatty liver to nonalcoholic steatohepatitis (NASH) and cirrhosis. Microsomal, S9, and cytosol fractions were extracted from 36 histologically normal and 66 NAFLD livers (27 nonalcoholic fatty liver [NAFL], 13 NASH, and 26 NASH with cirrhosis). Corrected microsomal protein per gram liver (MPPGL) progressively decreased with disease severity (26.8, 27.4, and 24.3 mg/g in NAFL, NASH, and NASH/cirrhosis, respectively, compared with 35.6 mg/g in normal livers; ANOVA, P < 0.001). Homogenate, S9, and cytosolic protein showed a consistent trend of decline in NASH/cirrhosis relative to normal control (post-hoc t test, P < 0.05). No differences across the groups were observed in homogenate, S9, cytosolic, and microsomal protein content in matched kidney samples. MPPGL-based scalars that combine protein content with liver size revealed that the reduction in MPPGL in NAFL and NASH was compensated by the reported increase in liver size (relative scalar ratios of 0.96 and 0.99, respectively), which was not the case with NASH/cirrhosis (ratio of 0.63), compared with healthy control. Physiologically based pharmacokinetics-informed global sensitivity analysis of the relative contribution of IVIVE scalars (hepatic CYP3A4 abundance, MPPGL, and liver size) to variability in exposure (area under the curve) to three CYP3A substrates (alprazolam, midazolam, and ibrutinib) revealed enzyme abundance as the most significant parameter, followed by MPPGL, whereas liver volume was the least impactful factor. SIGNIFICANCE STATEMENT: Nonalcoholic fatty liver disease-specific scalars necessary for extrapolation from in vitro systems to liver tissue are lacking. These are required in clearance prediction and dose selection in nonalcoholic fatty liver and steatohepatitis populations. Previously reported disease-driven changes have focused on cirrhosis, with no data on the initial stages of liver disease. The authors obtained experimental values for microsomal, cytosolic, and S9 fractions and assessed the relative impact of microsomal scalars on predicted exposure to substrate drugs using physiologically based pharmacokinetics.

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.

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

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

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.

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.