Skip to main content
Advertisement

Main menu

  • Home
  • Articles
    • Current Issue
    • Fast Forward
    • Latest Articles
    • Special Sections
    • Archive
  • Information
    • Instructions to Authors
    • Submit a Manuscript
    • FAQs
    • For Subscribers
    • Terms & Conditions of Use
    • Permissions
  • Editorial Board
  • Alerts
    • Alerts
    • RSS Feeds
  • Virtual Issues
  • Feedback
  • Submit
  • Other Publications
    • Drug Metabolism and Disposition
    • Journal of Pharmacology and Experimental Therapeutics
    • Molecular Pharmacology
    • Pharmacological Reviews
    • Pharmacology Research & Perspectives
    • ASPET

User menu

  • My alerts
  • Log in
  • My Cart

Search

  • Advanced search
Drug Metabolism & Disposition
  • Other Publications
    • Drug Metabolism and Disposition
    • Journal of Pharmacology and Experimental Therapeutics
    • Molecular Pharmacology
    • Pharmacological Reviews
    • Pharmacology Research & Perspectives
    • ASPET
  • My alerts
  • Log in
  • My Cart
Drug Metabolism & Disposition

Advanced Search

  • Home
  • Articles
    • Current Issue
    • Fast Forward
    • Latest Articles
    • Special Sections
    • Archive
  • Information
    • Instructions to Authors
    • Submit a Manuscript
    • FAQs
    • For Subscribers
    • Terms & Conditions of Use
    • Permissions
  • Editorial Board
  • Alerts
    • Alerts
    • RSS Feeds
  • Virtual Issues
  • Feedback
  • Submit
  • Visit dmd on Facebook
  • Follow dmd on Twitter
  • Follow ASPET on LinkedIn
Research ArticleArticle

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

Jill Barber, Zubida M. Al-Majdoub, Narciso Couto, Areti-Maria Vasilogianni, Annika Tillmann, Sarah Alrubia, Amin Rostami-Hodjegan and Brahim Achour
Drug Metabolism and Disposition June 2022, 50 (6) 762-769; DOI: https://doi.org/10.1124/dmd.121.000780
Jill Barber
Centre for Applied Pharmacokinetic Research, School of Health Sciences, University of Manchester, Manchester, United Kingdom (J.B., Z.M.A.-M., N.C., A.-M.V., A.T., S.A., A.R.-H., B.A.) Simcyp Division, Certara, Sheffield, United Kingdom (A.R.-H.) and Department of Biomedical and Pharmaceutical Sciences, College of Pharmacy, the University of Rhode Island, Kingston, Rhode Island (B.A.)
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Jill Barber
Zubida M. Al-Majdoub
Centre for Applied Pharmacokinetic Research, School of Health Sciences, University of Manchester, Manchester, United Kingdom (J.B., Z.M.A.-M., N.C., A.-M.V., A.T., S.A., A.R.-H., B.A.) Simcyp Division, Certara, Sheffield, United Kingdom (A.R.-H.) and Department of Biomedical and Pharmaceutical Sciences, College of Pharmacy, the University of Rhode Island, Kingston, Rhode Island (B.A.)
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Zubida M. Al-Majdoub
Narciso Couto
Centre for Applied Pharmacokinetic Research, School of Health Sciences, University of Manchester, Manchester, United Kingdom (J.B., Z.M.A.-M., N.C., A.-M.V., A.T., S.A., A.R.-H., B.A.) Simcyp Division, Certara, Sheffield, United Kingdom (A.R.-H.) and Department of Biomedical and Pharmaceutical Sciences, College of Pharmacy, the University of Rhode Island, Kingston, Rhode Island (B.A.)
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Narciso Couto
Areti-Maria Vasilogianni
Centre for Applied Pharmacokinetic Research, School of Health Sciences, University of Manchester, Manchester, United Kingdom (J.B., Z.M.A.-M., N.C., A.-M.V., A.T., S.A., A.R.-H., B.A.) Simcyp Division, Certara, Sheffield, United Kingdom (A.R.-H.) and Department of Biomedical and Pharmaceutical Sciences, College of Pharmacy, the University of Rhode Island, Kingston, Rhode Island (B.A.)
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Areti-Maria Vasilogianni
Annika Tillmann
Centre for Applied Pharmacokinetic Research, School of Health Sciences, University of Manchester, Manchester, United Kingdom (J.B., Z.M.A.-M., N.C., A.-M.V., A.T., S.A., A.R.-H., B.A.) Simcyp Division, Certara, Sheffield, United Kingdom (A.R.-H.) and Department of Biomedical and Pharmaceutical Sciences, College of Pharmacy, the University of Rhode Island, Kingston, Rhode Island (B.A.)
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Sarah Alrubia
Centre for Applied Pharmacokinetic Research, School of Health Sciences, University of Manchester, Manchester, United Kingdom (J.B., Z.M.A.-M., N.C., A.-M.V., A.T., S.A., A.R.-H., B.A.) Simcyp Division, Certara, Sheffield, United Kingdom (A.R.-H.) and Department of Biomedical and Pharmaceutical Sciences, College of Pharmacy, the University of Rhode Island, Kingston, Rhode Island (B.A.)
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Amin Rostami-Hodjegan
Centre for Applied Pharmacokinetic Research, School of Health Sciences, University of Manchester, Manchester, United Kingdom (J.B., Z.M.A.-M., N.C., A.-M.V., A.T., S.A., A.R.-H., B.A.) Simcyp Division, Certara, Sheffield, United Kingdom (A.R.-H.) and Department of Biomedical and Pharmaceutical Sciences, College of Pharmacy, the University of Rhode Island, Kingston, Rhode Island (B.A.)
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Amin Rostami-Hodjegan
Brahim Achour
Centre for Applied Pharmacokinetic Research, School of Health Sciences, University of Manchester, Manchester, United Kingdom (J.B., Z.M.A.-M., N.C., A.-M.V., A.T., S.A., A.R.-H., B.A.) Simcyp Division, Certara, Sheffield, United Kingdom (A.R.-H.) and Department of Biomedical and Pharmaceutical Sciences, College of Pharmacy, the University of Rhode Island, Kingston, Rhode Island (B.A.)
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Brahim Achour
  • Article
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF + SI
  • PDF
Loading

Abstract

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

SIGNIFICANCE STATEMENT This study provides useful insights for using label-free methods to generate abundance data applicable for populating pharmacokinetic models. The data demonstrated overall correlation between intensity-based label-free proteomic methods (HiN, iBAQ and TPA), whereas iBAQ and TPA overestimated the total amount of protein in the samples. The extent of overestimation can provide a means of normalization to support absolute quantification. Importantly, between-sample relative quantification was consistent (similar variability) across methods.

Introduction

Quantitative proteomics has become a standard method in molecular biology, with the central aim of measuring expression profiles at the protein level. Because of its broad scope of analysis (Wang et al., 2019), label-free (global) proteomics allows measurement of a wide range of proteins that govern drug pharmacokinetics (principally drug-metabolizing enzymes and drug transporters) and their changes in response to pathologic and environmental factors (El-Khateeb et al., 2019; Prasad et al., 2019). The methodology does not depend on the availability of isotopically labeled standards, which are expensive (Al Feteisi et al., 2015) and, when used to quantify low-abundance proteins, require additional care in data analysis (Achour et al., 2018). Mass spectrometry is not, however, an inherently quantitative technique; the relationship between the concentration of the analyte and the intensity of the corresponding signal is complex (Couto et al., 2011), and all label-free methods require assumptions that may not be fully justified (Arike et al., 2012).

Label-free measurement relies on signal intensity either of all native peptides [e.g., the total protein approach (TPA) (Wiśniewski and Rakus, 2014) and intensity-based absolute quantification (iBAQ) (Schwanhäusser et al., 2011)] or a set of unique/razor peptides [e.g., high N ion intensity approach (HiN) (Silva et al., 2006)] assigned to a certain protein. Alternatively, relative quantification can be achieved based on spectral counts [e.g., the exponentially modified protein abundance index (emPAI) (Ishihama et al., 2005)], as a semi-quantitative approach to derive an estimate of protein expression.

During the past 5 years, the pharmacology community has produced approximately 20 publications quantifying human tissue proteomes by global proteomics, and this number is set to increase in the next few years. There are ongoing debates about whether samples should be fractionated and about sample preparation methods (Prasad et al., 2019), although filter-aided sample preparation (FASP) methodology (Wiśniewski et al., 2009) is widely adopted when sample is plentiful. More surprisingly, there is no consensus about methods for analysis of the (typically gigabyte) RAW files obtained from mass spectrometry experiments. Different software packages for data analysis are available that make data processing more streamlined, but these do not always produce completely consistent results (Välikangas et al., 2018). Processing can be done using different reference datasets and following different assumptions, and, above all, using different quantification methodologies (El-Khateeb et al., 2021). These factors are particularly important when modeling the impact of different covariates, such as disease, is the focus of investigation, and therefore validating such data plays a major role in increasing trust in the outcome of predictive pharmacology models.

For the work described here, we used a well-characterized set of 27 human liver microsomal (HLM) samples. As previously described, sample preparation was done by standard FASP methodology and mass spectrometry was carried out on an Orbitrap HF QE instrument (Couto et al., 2019). RAW files so generated have been uploaded to the Proteomics Identification (PRIDE) database and are freely available (Al-Majdoub et al., 2020). The focus of this article is to evaluate the quality of label-free abundance data. In particular, we aimed to use these samples as controls for an investigation of experimentally demanding pediatric samples, analyzed with no standards, and we therefore required that standard-free methods should be as reliable as possible.

Materials and Methods

The proteomic dataset has been deposited to the Proteomics Identification (PRIDE) repository under the identifier PXD020910.

Samples and Proteomic Methodology

The preparation of HLM samples (Supplemental Table 1) and analysis by mass spectrometry are fully described elsewhere (Couto et al., 2019). Briefly, liver membrane fractions were prepared using differential centrifugation, first at low speed (10,000g) to separate cellular debris from the post-mitochondrial fraction, followed by high-speed centrifugation (100,000g) to isolate microsomes. Protein content was measured using the Bradford protein assay and sample preparation of 100 µg of each sample (n = 27) followed the FASP protocol with multienzyme digestion (lysyl endopeptidase and trypsin). Three exogenous protein standards were spiked in the samples: bovine serum albumin (BSA, 0.2 µg), bovine cytochrome c (0.15 µg) and equine myoglobin (0.3 µg). Peptides (1 µg) were analyzed by liquid chromatography-tandem mass spectrometry using an UltiMate 3000 rapid separation liquid chromatography system (Dionex, Surrey, UK) coupled to a Q Exactive HF Hybrid Quadrupole-Orbitrap mass spectrometer (ThermoFisher Scientific, Bremen, Germany).

RAW files were processed using Progenesis version 4.2 as a single batch, and the resulting mgf files were processed using Mascot version 2.7 for protein identification. The human database used was Uniprot 000005640, containing 77,027 protein entries. Files were processed four times with the settings: enzyme trypsin/P, MS tolerance 5 ppm, MS/MS tolerance 0.02 Da. In the initial run, up to one missed cleavage was allowed, carbamidomethyl (cysteine) was set as a fixed modification, and oxidation (methionine) was the only variable modification. In subsequent runs, a second missed cleavage, deamidation (of asparagine and glutamine) and phosphorylation (of serine, threonine and tyrosine) were (separately) permitted. Results were compiled using Progenesis and exported as csv files.

Data Analysis and Label-Free Quantification

Peptides were assigned to proteins based on a bespoke razor as described previously (Al-Majdoub et al., 2020) using Microsoft Excel 365. Assignment prioritized full-length characterized sequences over truncated, uncharacterized, and cDNA sequences. A best-fit analysis was then run to minimize the number of proteins assigned to account for all of the peptides. Deamidated peptides with no corresponding native assignment and those that did not match any protein were deleted. The peptide MS intensities attributed to more than one protein were divided among those proteins based on the ratio of unique or razor (peptides with a single assignment within the current dataset) peptide intensities for each protein, as detailed previously (Al-Majdoub et al., 2020).

Three potential standards (BSA, bovine cytochrome c, and equine myoglobin) were assessed. The equations used in quantification by the HiN, iBAQ, and TPA methods are detailed below (eqs. 1–3). Rearrangement of these equations provides the means to compare the total sample estimated from the total intensity against the total sample analyzed. These act as sanity checks on data analysis.

  • High N (HiN) ion intensity method:

Embedded Image

Where [Protein] is the abundance of a target protein, [standard] is the abundance of the standard protein, both expressed in units of pmol mg−1 total protein, and the fraction refers to the ratio of the average intensity of the n/m highest ion peaks of the target protein relative to the standard (in this case, n = m = 2 or 3). Peptides used for quantification are unique to the target proteins; other selection criteria were according to (Achour et al., 2018).

  • Intensity-based absolute quantification (iBAQ):

Embedded Image

Where the summed intensity of all peptides i from the protein of interest j or the standard k is normalized to T, the number of theoretically observable peptides from a tryptic digest of protein j or standard k.

  • The total protein approach (TPA):

Embedded Image

Where the ratio of the sum of intensity of all peptides i derived from a protein j of interest to the sum of intensity of all peptides (from all proteins) in a sample (expressed in parts per billion) is converted to an abundance value (pmol mg−1) by normalizing to the molecular mass of the protein in Daltons.

Statistical Data Analysis

Data were expressed as mean and standard deviation, and variability was assessed as CV and fold difference (maximum-to-minimum ratio). Abundance and activity correlations were tested by linear regression (R2). Relationships between abundance data and either age or body mass index (BMI) were tested using Pearson correlation (r) to show the direction of trends. Differences between abundance data generated using label-free and targeted methods and between male and female donors were assessed using a t test. Differences across genotypes were assessed using either one-way ANOVA or a t test. A probability cutoff of 0.05 was set for statistical significance.

Results

In this study, HLM samples were analyzed using global proteomic methods. Percent identical peptides reflected high integrity of analyses (86–99%) across replicates and across samples (Supplemental Tables 2 and 3). For the purpose of quantification, we chose to assess three MS intensity-based label-free methods on the basis that they provide more robust protein measurements than spectral counts (Arike et al., 2012). Two methods (HiN and iBAQ) rely on exogenous protein standards at known concentrations, whereas the TPA is applied without the use of a standard. The methods allowed quantification of 2232 proteins, and data describing expression of 23 cytochrome P450 enzymes, 11 glucuronosyltransferases (UGT), 17 ABC transporters, and 6 solute carriers are shown in Supplemental Tables 4-6.

Choice of Standard

For the iBAQ and HiN methods, three potential standards were included in the samples. The total amount of protein used in each experiment was 100 µg, and the amounts of standards were also known: BSA 0.2 µg (28.86 pmol mg−1 total protein), myoglobin 0.3 µg (175.61 pmol mg−1 total protein) and cytochrome c 0.15 µg (32.04 pmol mg−1 total protein). BSA was expected to give the best results because its molecular mass (69 kDa) is close to the average molecular mass of the detected proteins (60 kDa), and it yields a high number of unique peptides. Cytochrome c gave rise to a limited number of unique peptides and was therefore discarded. The HiN method was used to calculate the total analyte protein using BSA and myoglobin (Table 1 and Fig. 1A), with results of 53% for BSA and 207% for myoglobin compared with the nominal amount of protein analyzed (assuming an average molecular mass of 60 kDa for native proteins). In the HiN method, proteins represented by a single peptide as well as those falling below the limit of quantification are ignored, making 53% a reasonable value and 207% a substantial overestimate. Further calculations were therefore carried out using BSA as a standard.

Fig. 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 1.

Comparison of liver proteome measurements using three label-free methods. (A) Percentage of the measured total protein content relative to the nominal content (dashed line) analyzed by mass spectrometry. The amount was determined using the TPA, HiN [based on either myoglobin (MYG) or BSA as standards] and iBAQ (based on BSA). (B) Head-to-head comparison of average concentrations (of 2232 proteins) in 27 samples quantified by iBAQ (BSA) and TPA compared with HiN (BSA). (C) Correlation between mean concentrations of 57 key pharmacokinetic targets measured by TPA or iBAQ and mean abundances measured by HiN. The data show that TPA overestimates protein amounts compared with the HiN method, whereas HiN and iBAQ methods produce comparable results in most cases. The data also indicate that it is possible to estimate absolute abundances by iBAQ or TPA using conversion factors. Abundance is expressed in units of pmol mg−1 total protein.

View this table:
  • View inline
  • View popup
TABLE 1

The total amount of protein in analyte estimated by the label-free quantification methods, averaged (± SD) over 27 samples.

Comparison Between Label-Free Quantification Methods

Although iBAQ and TPA clearly overestimate the total amount of protein (Fig. 1A), both enable estimation of the abundance of proteins that give rise to a single detectable peptide, whereas the HiN method does not. Our previous work (El-Khateeb et al., 2021) indicated that all three methods perform quite well in assessing relative change from healthy baseline. We therefore investigated the correlation between absolute quantification values obtained using the three methods. An overall picture of correlations between measurements (of 2232 proteins) using the different methods is presented in Fig. 1B. The overall correlation between the three methods is strong (R2 > 0.70, P < 0.001). The TPA overestimates by a factor of two- to threefold relative to HiN. The iBAQ and HiN measurements are relatively comparable. Correlations between mean abundances of enzymes and transporters (n = 57) are presented in Fig. 1C. Individual abundance data for these targets are presented in Supplemental Tables 4-6. Specific cases related to a number of proteins of interest to drug pharmacokinetics are shown in Fig. 2. Figure 2 shows that in all cases a straight line can be drawn connecting iBAQ and HiN quantifications (R2 = 0.60–0.97) and so is the case for TPA and HiN but with generally considerably more scatter (R2 = 0.25–0.97). Generally, the TPA gives the highest estimation of the concentration of a protein; for lower abundance proteins (e.g., low abundance transporters), TPA gives the lowest estimation. Relative abundances are presented in Supplemental Fig. 1 for cytochrome P450 enzymes, UGTs, and ABC transporters, reflecting overall agreement, except in the case of relative abundance of ABC transporters determined using the TPA method.

Fig. 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 2.

Correlation between protein concentrations of key pharmacokinetic targets measured by iBAQ (blue) or TPA (orange) relative to HiN method in 27 liver samples. The data show examples of drug-metabolizing enzymes and transporters. BSA was used as a standard for HiN and iBAQ, and abundance was measured in units of pmol mg−1 total protein.

Correlation of Label-Free Data with Functional Activity

Functional activity data were available for several cytochrome P450 and UGT enzymes (Achour et al., 2014, 2017). Correlations between abundance and activity of cytochrome P450 3A4, 2D6, 1A2, 2B6 and 2C19 were moderate to strong across the three methods (R2 = 0.56–0.88 with HiN, R2 = 0.57–0.91 with iBAQ, and R2 = 0.63–0.88 with TPA; Fig. 3A). Correlation with CYP2C9 activity was the exception, with different degrees of correlation across methods (R2 = 0.23 with HiN, R2 = 0.43 with iBAQ, and R2 = 0.71 with TPA). Similarly, weak correlation was previously reported for CYP2C9 activity against diclofenac with targeted data in the set of samples (Achour et al., 2014).

Fig. 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 3.

Correlation of protein concentrations of (A) cytochrome P450 and (B) UGT enzymes measured by HiN (red), iBAQ (blue), and TPA (orange) against functional activity in 27 liver samples. Activity was measured with metabolite formation assays against the substrates: phenacetin (CYP1A2), mephenytoin (CYP2B6), diclofenac (CYP2C9), mephenytoin (CYP2C19), bufuralol (CYP2D6), testosterone (CYP3A), β-estradiol (UGT1A1), chenodeoxycholic acid (UGT1A3), 5-hydroxytryptophol (UGT1A6), propofol (UGT1A9), zidovudine (UGT2B7) and S-oxazepam (UGT2B15). Abundance was measured in units of pmol mg−1 total protein; catalytic activity was measured in units of nmol metabolite min−1 mg−1 total protein.

Correlations between abundance and activity of UGTs 1A1, 1A3, 1A6, and 1A9 were moderate to strong (R2 = 0.34–0.69 with HiN, R2 = 0.42–0.85 with iBAQ, and R2 = 0.52–0.77 with TPA; Fig. 3B). The exception was UGT2B7, with different degrees of correlation across methods (R2 = 0.24 with HiN, R2 = 0.44 with iBAQ, and R2 = 0.59 with TPA), whereas correlations of UGT2B15 were generally weaker (R2 = 0.15, 0.23, and 0.39 with HiN, iBAQ, and TPA, respectively). Moderate correlations were previously reported for UGTs 2B7 and 2B15 activity with targeted data for the set of samples (Achour et al., 2017).

Comparison of Label-Free Quantification with Targeted Data

Label-free measurements were compared with previously reported targeted data for cytochrome P450 and UGT enzymes in the same set of samples (Achour et al., 2014, 2017). For cytochrome P450 enzymes, overall agreement was observed with HiN and iBAQ data (Fig. 4A), reflecting 74% of measurements within twofold of targeted data (Fig. 4C). TPA, however, tended to overestimate measurements (only 53% of the data were within twofold). For UGT enzymes, TPA measurements were closer to targeted data (93% of measurements were within twofold), whereas iBAQ and HiN tended to underestimate, with 77% and 60% of measurements within twofold, respectively.

Fig. 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 4.

Comparison of abundances of (A) cytochrome P450 and (B) UGT enzymes measured using label-free methods (HiN, iBAQ, and TPA) against targeted data. Ratios of label-free measurements relative to targeted data for (C) cytochrome P450 and (D) UGT enzymes. In (A) and (B), the whiskers represent the minimum-to-maximum range, the boxes represent the 25th and 75th percentiles, the lines represent the medians, and + signs represent the means. Comparisons based on a t test against targeted data are shown in black and against HiN measurements are shown in red. * P < 0.05, ** P < 0.01, and *** P < 0.001. In (C) and (D), the dashed lines denote the twofold range, and the percentages are the proportion of label-free measurements within twofold of targeted data. Abundance was measured in units of pmol mg−1 total protein.

Assessment of Variability with Label-Free Methods

Fold difference and CV, related to between-sample variability, was very similar across the three methods for all measured proteins (Fig. 5), indicating robustness of relative quantification regardless of methodology. The calculated CV combines technical and biologic variability. Technical variability assessed using a pool of the same set of samples returned values <30% for the targets across all methods. The calculated variability related to technical error, expressed as fold difference between the 5th and 95th percentiles Embedded Image, was therefore within fourfold, whereas total variability reflected up to 50-fold in abundance across the three methods.

Fig. 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 5.

Comparison of between-sample variability across 27 liver samples in the abundance of enzymes and transporters measured by HiN, iBAQ, and TPA, expressed as (A) fold difference (maximum-to-minimum ratio) and (B) percent CV (CV = SD × 100/mean). BSA was used as a standard for HiN and iBAQ.

Covariates of Protein Expression Assessed by Label-Free Methods

Donors’ demographic and clinical information is summarized in Supplemental Table 1. In addition, CYP2B6, CYP2C9, CYP2C19, CYP2D6, and CYP3A5 genotype data were available for 25 (out of 27) samples. Abundance data were assessed against sex, age, and BMI (Supplemental Table 7). The number of confirmed smokers and alcohol users at the time of donation was small (four smokers and three drinkers), and therefore, the effect of these two factors was not probed. No differences in abundance in samples from male (n = 15) and female (n = 12) donors were observed with all three label-free methods (t test, P > 0.05). Weak, negative correlations with age were revealed for CYP2C18, UGT2B4, UGT2B10, and ABCA1 with borderline significance across the label-free methods (r = −0.42 to −0.39, P = 0.03 to 0.05). The effect of BMI was moderate in the cases of UGT1A3, UGT1A4, and OATP1B3, with lower abundance in overweight and obese donors (r = −0.59 to −0.42, P = 0.001 to 0.03). The effect of genotype was significant in the cases of CYP2D6 (ANOVA, P < 0.05), CYP2C19 (t test, P = 0.01) and CYP3A5 (t test, P < 0.05).

Discussion

This study aimed to assess measurements of hepatic enzymes and transporters by widely used label-free proteomic methods (HiN, iBAQ and TPA). We have previously outlined the use of the Disease Perturbation Factor (DPF) (El-Khateeb et al., 2021), which is essentially a ratio connecting the amount of any given protein in a diseased tissue with the amount of the same protein in healthy control tissue. The DPF has been shown to be independent of the quantification methodology (targeted versus global proteomics, HiN versus iBAQ versus TPA). Furthermore, differences in absolute abundance are explored herein, and a key piece of information—the total analyte protein—which is usually discarded in proteomic data analysis, now allows us to adjudicate between the different methods of label-free quantification and even to estimate conversion factors from one method to another.

Absolute abundance correlated across the three methods, with targeted data and with functional activity. Correlation with targeted proteomics was also previously reported (Vildhede et al., 2018; Wiśniewski et al., 2019; El-Khateeb et al., 2021). Although the HiN method (Silva et al., 2006) generally produces data that seem biologically sensible, especially when BSA is used as a standard for human samples, it has two drawbacks. First, the N in HiN is generally taken to mean two, three, or more, but not one. Thus, we are denied even an estimate of the abundances of proteins represented by a single peptide. Second, a high-quality standard is required, one that should have similar properties to the proteins under study. The choice of a suitable standard in this study was, however, empirical as highlighted by the discrepancy between abundances against BSA and myoglobin, with myoglobin clearly overestimating the total amount of protein. Prospectively, it is, of course, possible to include a protein standard at appropriate concentration in new samples; it is not, however, possible to do this retrospectively. Proteomic work in the drug metabolism and disposition arena involves the use of precious, small human samples. It is imperative, both scientifically and ethically, to derive maximum information from each sample, which means that historical samples, prepared by suboptimal protocols, are still of value (Prasad et al., 2019).

The TPA has proved an excellent method for dealing with such samples. No standards are necessary, and the approach provides broad coverage by allowing the quantification (albeit with low accuracy) of proteins represented by a single peptide. Nonunique peptides can be accommodated within the analysis. We have introduced small modifications to the data analysis so that these nonunique peptides are not overrepresented (Al-Majdoub et al., 2020), but the method is still inclined to overestimate protein concentration relative to other label-free methods. This is not surprising. The normalization in a TPA experiment is based on the total signal intensity, but we know that some signal (that due to proteins falling below the limit of quantification) is not measured. A ‘proteomic ruler’ incorporating MS signal of cellular histones was introduced to make TPA measurements more biologically sensible (Wiśniewski et al., 2014). The DPF (El-Khateeb et al., 2021), being a relative factor, allows for the use of the TPA without too much concern for the systematic overestimate.

Both the iBAQ and TPA methods overestimate the total amount of protein in the sample, and the extent of this overestimation can provide a means of normalization should absolute quantification be required. Robust quantification is biased toward higher abundance proteins, and therefore such normalization approach may only work with enzymes and highly abundant transporters. Relative quantification, both between sample and within sample, is often more pertinent than absolute estimates. Similar variability (CV) across individual measurements recovered by all methods indicates relative quantification is robust, regardless of the quantification approach provided consistent proteomic workflows are used (El-Khateeb et al., 2021; Neuhoff et al., 2021). Relative data are particularly useful for assigning stoichiometry in protein expression (Fabre et al., 2014). The current set reveals a particular example, stoichiometry of TAP1 (ABCB2) and TAP2 (ABCB3). The correlation between abundances generated by the methods for these two proteins was excellent (R2 = 0.93–0.97), indicating strong agreement, and TAP2 to TAP1 ratio (average TAP2:TAP1 > 10) was within threefold across methods. TAP1 and TAP2 form a functional heterodimer that transports peptides for antigen presentation, and one would therefore expect 1:1 expression ratio. The clearly higher abundance of TAP2 may indicate divergence in regulatory mechanisms between TAP1 and TAP2, in support of previous observations (Bahram et al., 1991; Zeidler et al., 1997). Application of relative quantification can be useful to derive changes in abundance of enzymes in a disease population compared with healthy volunteers and assess the implications of such changes for drug–drug interactions.

In conclusion, historical samples without appropriate standards can be subject to quantification using the TPA method, with the expectation that similar results would be achieved by standard-based methods. Normalization, or simply adjustment by a factor of two to three leads to estimates of absolute quantification. Where standards are available, BSA is a good choice from readily available purified proteins. Importantly, relative quantification is robust across methods, which allows consistent assignment of between subject variability.

Acknowledgments

The authors thank Pfizer (Groton, CT) for the provision of microsomal samples along with demographic, clinical, genotype, and activity data. The ChELSI Institute, the University of Sheffield, and the Manchester Institute of Biotechnology, the University of Manchester, provided access to mass spectrometry instrumentation.

Authorship Contributions

Participated in research design: Barber, Achour.

Performed data analysis: Barber, Al-Majdoub, Couto, Achour.

Wrote or contributed to the writing of the manuscript: Barber, Al-Majdoub, Couto, Vasilogianni, Tillmann, Alrubia, Rostami-Hodjegan, Achour.

Footnotes

    • Received November 14, 2021.
    • Accepted March 4, 2022.
  • Supported by the Centre for Applied Pharmacokinetic Research (CAPKR) consortium (Z.M.A.-M., B.A.), Cancer Research UK (A.-M.V.), the Saudi Ministry of Education (S.A.), and DrugTrain funded by the European Union’s Horizon 2020 research and innovation program (A.T.).

  • dx.doi.org/10.1124/dmd.121.000780.

  • ↵Embedded ImageThis article has supplemental material available at dmd.aspetjournals.org.

Abbreviations

BSA
bovine serum albumin
FASP
filter-aided sample preparation
HiN
high N ion intensity method
HLM
human liver microsomes
iBAQ
intensity-based absolute quantification
TPA
total protein approach
  • Copyright © 2022 by The American Society for Pharmacology and Experimental Therapeutics

References

  1. ↵
    1. Achour B,
    2. Dantonio A,
    3. Niosi M,
    4. Novak JJ,
    5. Al-Majdoub ZM,
    6. Goosen TC,
    7. Rostami-Hodjegan A, and
    8. Barber J
    (2018) Data generated by quantitative liquid chromatography-mass spectrometry proteomics are only the start and not the endpoint: optimization of quantitative concatemer-based measurement of hepatic uridine-5′-diphosphate-glucuronosyltransferase enzymes with reference to catalytic activity. Drug Metab Dispos 46:805–812.
    OpenUrlAbstract/FREE Full Text
  2. ↵
    1. Achour B,
    2. Dantonio A,
    3. Niosi M,
    4. Novak JJ,
    5. Fallon JK,
    6. Barber J,
    7. Smith PC,
    8. Rostami-Hodjegan A, and
    9. Goosen TC
    (2017) Quantitative characterization of major hepatic UDP-glucuronosyltransferase enzymes in human liver microsomes: comparison of two proteomic methods and correlation with catalytic activity. Drug Metab Dispos 45:1102–1112.
    OpenUrlAbstract/FREE Full Text
  3. ↵
    1. Achour B,
    2. Russell MR,
    3. Barber J, and
    4. Rostami-Hodjegan A
    (2014) Simultaneous quantification of the abundance of several cytochrome P450 and uridine 5′-diphospho-glucuronosyltransferase enzymes in human liver microsomes using multiplexed targeted proteomics. Drug Metab Dispos 42:500–510.
    OpenUrlAbstract/FREE Full Text
  4. ↵
    1. Al-Majdoub ZM,
    2. Achour B,
    3. Couto N,
    4. Howard M,
    5. Elmorsi Y,
    6. Scotcher D,
    7. Alrubia S,
    8. El-Khateeb E,
    9. Vasilogianni AM,
    10. Alohali N et al.
    (2020) Mass spectrometry-based abundance atlas of ABC transporters in human liver, gut, kidney, brain and skin. FEBS Lett 594:4134–4150.
    OpenUrl
  5. ↵
    1. Al Feteisi H,
    2. Achour B,
    3. Barber J, and
    4. Rostami-Hodjegan A
    (2015) Choice of LC-MS methods for the absolute quantification of drug-metabolizing enzymes and transporters in human tissue: a comparative cost analysis. AAPS J 17:438–446.
    OpenUrlCrossRefPubMed
  6. ↵
    1. Arike L,
    2. Valgepea K,
    3. Peil L,
    4. Nahku R,
    5. Adamberg K, and
    6. Vilu R
    (2012) Comparison and applications of label-free absolute proteome quantification methods on Escherichia coli. J Proteomics 75:5437–5448.
    OpenUrlCrossRefPubMed
  7. ↵
    1. Bahram S,
    2. Arnold D,
    3. Bresnahan M,
    4. Strominger JL, and
    5. Spies T
    (1991) Two putative subunits of a peptide pump encoded in the human major histocompatibility complex class II region. Proc Natl Acad Sci USA 88:10094–10098.
    OpenUrlAbstract/FREE Full Text
  8. ↵
    1. Couto N,
    2. Al-Majdoub ZM,
    3. Achour B,
    4. Wright PC,
    5. Rostami-Hodjegan A, and
    6. Barber J
    (2019) Quantification of proteins involved in drug metabolism and disposition in the human liver using label-free global proteomics. Mol Pharm 16:632–647.
    OpenUrl
  9. ↵
    1. Couto N,
    2. Barber J, and
    3. Gaskell SJ
    (2011) Matrix-assisted laser desorption/ionisation mass spectrometric response factors of peptides generated using different proteolytic enzymes. J Mass Spectrom 46:1233–1240.
    OpenUrlPubMed
  10. ↵
    1. El-Khateeb E,
    2. Al-Majdoub ZM,
    3. Rostami-Hodjegan A,
    4. Barber J, and
    5. Achour B
    (2021) Proteomic quantification of changes in abundance of drug-metabolizing enzymes and drug transporters in human liver cirrhosis: different methods, similar outcomes. Drug Metab Dispos 49:610–618.
    OpenUrlAbstract/FREE Full Text
  11. ↵
    1. El-Khateeb E,
    2. Vasilogianni A-M,
    3. Alrubia S,
    4. Al-Majdoub ZM,
    5. Couto N,
    6. Howard M,
    7. Barber J,
    8. Rostami-Hodjegan A, and
    9. Achour B
    (2019) Quantitative mass spectrometry-based proteomics in the era of model-informed drug development: applications in translational pharmacology and recommendations for best practice. Pharmacol Ther 203:107397.
    OpenUrl
  12. ↵
    1. Fabre B,
    2. Lambour T,
    3. Bouyssié D,
    4. Menneteau T,
    5. Monsarrat B,
    6. Burlet-Schiltz O, and
    7. Bousquet-Dubouch M-P
    (2014) Comparison of label-free quantification methods for the determination of protein complexes subunits stoichiometry. EuPA Open Proteom 4:82–86.
    OpenUrl
  13. ↵
    1. Ishihama Y,
    2. Oda Y,
    3. Tabata T,
    4. Sato T,
    5. Nagasu T,
    6. Rappsilber J, and
    7. Mann M
    (2005) Exponentially modified protein abundance index (emPAI) for estimation of absolute protein amount in proteomics by the number of sequenced peptides per protein. Mol Cell Proteomics 4:1265–1272.
    OpenUrlAbstract/FREE Full Text
  14. ↵
    1. Neuhoff S,
    2. Harwood MD,
    3. Rostami-Hodjegan A, and
    4. Achour B
    (2021) Application of proteomic data in the translation of in vitro observations to associated clinical outcomes. Drug Discov Today Technol 39:13–22.
    OpenUrl
  15. ↵
    1. Prasad B,
    2. Achour B,
    3. Artursson P,
    4. Hop CECA,
    5. Lai Y,
    6. Smith PC,
    7. Barber J,
    8. Wisniewski JR,
    9. Spellman D,
    10. Uchida Y et al.
    (2019) Toward a consensus on applying quantitative liquid chromatography-tandem mass spectrometry proteomics in translational pharmacology research: a white paper. Clin Pharmacol Ther 106:525–543.
    OpenUrl
  16. ↵
    1. Schwanhäusser B,
    2. Busse D,
    3. Li N,
    4. Dittmar G,
    5. Schuchhardt J,
    6. Wolf J,
    7. Chen W, and
    8. Selbach M
    (2011) Global quantification of mammalian gene expression control. Nature 473:337–342.
    OpenUrlCrossRefPubMed
  17. ↵
    1. Silva JC,
    2. Gorenstein MV,
    3. Li G-Z,
    4. Vissers JPC, and
    5. Geromanos SJ
    (2006) Absolute quantification of proteins by LCMSE: a virtue of parallel MS acquisition. Mol Cell Proteomics 5:144–156.
    OpenUrlAbstract/FREE Full Text
    1. Välikangas T,
    2. Suomi T, and
    3. Elo LL
    (2018) A comprehensive evaluation of popular proteomics software workflows for label-free proteome quantification and imputation. Brief Bioinform 19:1344–1355.
    OpenUrlCrossRef
  18. ↵
    1. Vildhede A,
    2. Nguyen C,
    3. Erickson BK,
    4. Kunz RC,
    5. Jones R,
    6. Kimoto E,
    7. Bourbonais F,
    8. Rodrigues AD, and
    9. Varma MVS
    (2018) Comparison of proteomic quantification approaches for hepatic drug transporters: multiplexed global quantitation correlates with targeted proteomic quantitation. Drug Metab Dispos 46:692–696.
    OpenUrlAbstract/FREE Full Text
  19. ↵
    1. Wang D,
    2. Eraslan B,
    3. Wieland T,
    4. Hallström B,
    5. Hopf T,
    6. Zolg DP,
    7. Zecha J,
    8. Asplund A,
    9. Li LH,
    10. Meng C et al.
    (2019) A deep proteome and transcriptome abundance atlas of 29 healthy human tissues. Mol Syst Biol 15:e8503.
    OpenUrlAbstract/FREE Full Text
  20. ↵
    1. Wiśniewski JR,
    2. Hein MY,
    3. Cox J, and
    4. Mann M
    (2014) A “proteomic ruler” for protein copy number and concentration estimation without spike-in standards. Mol Cell Proteomics 13:3497–3506.
    OpenUrlAbstract/FREE Full Text
  21. ↵
    1. Wiśniewski JR and
    2. Rakus D
    (2014) Multi-enzyme digestion FASP and the ‘total protein approach’-based absolute quantification of the Escherichia coli proteome. J Proteomics 109:322–331.
    OpenUrl
  22. ↵
    1. Wiśniewski JR,
    2. Wegler C, and
    3. Artursson P
    (2019) Multiple-enzyme-digestion strategy improves accuracy and sensitivity of label- and standard-free absolute quantification to a level that is achievable by analysis with stable isotope-labeled standard spiking. J Proteome Res 18:217–224.
    OpenUrl
  23. ↵
    1. Wiśniewski JR,
    2. Zougman A,
    3. Nagaraj N, and
    4. Mann M
    (2009) Universal sample preparation method for proteome analysis. Nat Methods 6:359–362.
    OpenUrlCrossRefPubMed
  24. ↵
    1. Zeidler R,
    2. Eissner G,
    3. Meissner P,
    4. Uebel S,
    5. Tampé R,
    6. Lazis S, and
    7. Hammerschmidt W
    (1997) Downregulation of TAP1 in B lymphocytes by cellular and Epstein-Barr virus-encoded interleukin-10. Blood 90:2390–2397.
    OpenUrlAbstract/FREE Full Text
PreviousNext
Back to top

In this issue

Drug Metabolism and Disposition: 50 (6)
Drug Metabolism and Disposition
Vol. 50, Issue 6
1 Jun 2022
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover
  • Index by author
  • Editorial Board (PDF)
  • Front Matter (PDF)
Download PDF
Article Alerts
Sign In to Email Alerts with your Email Address
Email Article

Thank you for sharing this Drug Metabolism & Disposition article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Label-Free but Still Constrained: Assessment of Global Proteomic Strategies for the Quantification of Hepatic Enzymes and Transporters
(Your Name) has forwarded a page to you from Drug Metabolism & Disposition
(Your Name) thought you would be interested in this article in Drug Metabolism & Disposition.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
Research ArticleArticle

Critique of Label-Free Proteomics of Pharmacokinetic Targets

Jill Barber, Zubida M. Al-Majdoub, Narciso Couto, Areti-Maria Vasilogianni, Annika Tillmann, Sarah Alrubia, Amin Rostami-Hodjegan and Brahim Achour
Drug Metabolism and Disposition June 1, 2022, 50 (6) 762-769; DOI: https://doi.org/10.1124/dmd.121.000780

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero

Share
Research ArticleArticle

Critique of Label-Free Proteomics of Pharmacokinetic Targets

Jill Barber, Zubida M. Al-Majdoub, Narciso Couto, Areti-Maria Vasilogianni, Annika Tillmann, Sarah Alrubia, Amin Rostami-Hodjegan and Brahim Achour
Drug Metabolism and Disposition June 1, 2022, 50 (6) 762-769; DOI: https://doi.org/10.1124/dmd.121.000780
Reddit logo Twitter logo Facebook logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Introduction
    • Materials and Methods
    • Results
    • Discussion
    • Acknowledgments
    • Authorship Contributions
    • Footnotes
    • Abbreviations
    • References
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF + SI
  • PDF

Related Articles

Cited By...

More in this TOC Section

  • Functional Characterization of 29 CYP4F2 Variants
  • Exposure-toxicity relation of apatinib
  • ABC phenomenon potentiates anti-HCC efficacy
Show more Articles

Similar Articles

Advertisement
  • Home
  • Alerts
Facebook   Twitter   LinkedIn   RSS

Navigate

  • Current Issue
  • Fast Forward by date
  • Fast Forward by section
  • Latest Articles
  • Archive
  • Search for Articles
  • Feedback
  • ASPET

More Information

  • About DMD
  • Editorial Board
  • Instructions to Authors
  • Submit a Manuscript
  • Customized Alerts
  • RSS Feeds
  • Subscriptions
  • Permissions
  • Terms & Conditions of Use

ASPET's Other Journals

  • Journal of Pharmacology and Experimental Therapeutics
  • Molecular Pharmacology
  • Pharmacological Reviews
  • Pharmacology Research & Perspectives
ISSN 1521-009X (Online)

Copyright © 2023 by the American Society for Pharmacology and Experimental Therapeutics