Hepatocytes play a pivotal role in human homeostasis. They are essential in regulation of glucose and lipid levels in blood and play a central role in metabolism of amino acids, lipids, drugs and ...xenobiotic-compounds. In addition, hepatocytes produce a major portion of proteins circulating in the blood. Hepatocytes were isolated from liver tissue obtained from surgical resections. Proteins were extracted and processed using filter aided sample preparation protocol and were analyzed by LC–MS/MS using high accuracy mass spectrometry. Proteins were quantified by the ‘Total Protein Approach’ and ‘Proteomic Ruler’. We report a comprehensive proteomic analysis of purified human hepatocytes and the human hepatoma cell line HepG2. The complete dataset comprises 9400 proteins and provides a comprehensive and quantitative depiction of the proteomes of hepatocytes and HepG2 cells at the protein titer and copy number dimensions. We describe basic cell organization and in detail energy metabolism pathways and metabolite transport. We provide quantitative insights into protein synthesis and drug and xenobiotics catabolism. Our data delineate differences between the native human hepatocytes and HepG2 cells by providing for the first time quantitative data at protein concentrations and copy numbers.
Display omitted
•In-depth quantitative proteome analysis of the human hepatocyte•Wide-ranging differences between the human hepatocyte and hepatoma cell line HepG2•Reference datasets for hepatocellular research
The hepatic organic anion transporting polypeptides (OATPs) influence the pharmacokinetics of several drug classes and are involved in many clinical drug–drug interactions. Predicting potential ...interactions with OATPs is, therefore, of value. Here, we developed in vitro and in silico models for identification and prediction of specific and general inhibitors of OATP1B1, OATP1B3, and OATP2B1. The maximal transport activity (MTA) of each OATP in human liver was predicted from transport kinetics and protein quantification. We then used MTA to predict the effects of a subset of inhibitors on atorvastatin uptake in vivo. Using a data set of 225 drug-like compounds, 91 OATP inhibitors were identified. In silico models indicated that lipophilicity and polar surface area are key molecular features of OATP inhibition. MTA predictions identified OATP1B1 and OATP1B3 as major determinants of atorvastatin uptake in vivo. The relative contributions to overall hepatic uptake varied with isoform specificities of the inhibitors.
Differences in the expression and function of the organic anion transporting polypeptide (OATP) transporters contribute to interindividual variability in atorvastatin clearance. However, the ...importance of the bile acid transporter sodium taurocholate cotransporting polypeptide (NTCP, SLC10A1) in atorvastatin uptake clearance (CLupt) is not yet clarified. To elucidate this issue, we investigated the relative contribution of NTCP, OATP1B1, OATP1B3, and OATP2B1 to atorvastatin CLupt in 12 human liver samples. The impact of inhibition on atorvastatin CLupt was also studied, using inhibitors of different isoform specificities. Expression levels of the four transport proteins were quantified by liquid chromatography tandem mass spectrometry. These data, together with atorvastatin in vitro kinetics, were used to predict the maximal transport activity (MTA) and interindividual differences in CLupt of each transporter in vivo. Subsequently, hepatic uptake impairment on coadministration of five clinically interacting drugs was predicted using in vitro inhibitory potencies. NTCP and OATP protein expression varied 3.7- to 32-fold among the 12 sample donors. The rank order in expression was OATP1B1 > OATP1B3 ≈ NTCP ≈ OATP2B1. NTCP was found to be of minor importance in atorvastatin disposition. Instead, OATP1B1 and OATP1B3 were confirmed as the major atorvastatin uptake transporters. The average contribution to atorvastatin uptake was OATP1B1 > OATP1B3 >> OATP2B1 > NTCP, although this rank order varied among individuals. The interindividual differences in transporter expression and CLupt resulted in marked differences in drug-drug interactions due to isoform-specific inhibition. We conclude that this variation should be considered in in vitro to in vivo extrapolations.
3D cultures of primary human hepatocytes (PHH) are emerging as a more in vivo-like culture system than previously available hepatic models. This work describes the characterisation of drug metabolism ...in 3D PHH spheroids. Spheroids were formed from three different donors of PHH and the expression and activities of important cytochrome P450 enzymes (CYP1A2, 2B6, 2C9, 2D6, and 3A4) were maintained for up to 21 days after seeding. The activity of CYP2B6 and 3A4 decreased, while the activity of CYP2C9 and 2D6 increased over time (P < 0.05). For six test compounds, that are metabolised by multiple enzymes, intrinsic clearance (CLint) values were comparable to standard in vitro hepatic models and successfully predicted in vivo CLint within 3-fold from observed values for low clearance compounds. Remarkably, the metabolic turnover of these low clearance compounds was reproducibly measured using only 1–3 spheroids, each composed of 2000 cells. Importantly, metabolites identified in the spheroid cultures reproduced the major metabolites observed in vivo, both primary and secondary metabolites were captured. In summary, the 3D PHH spheroid model shows promise to be used in drug discovery projects to study drug metabolism, including unknown mechanisms, over an extended period of time.
Transporter-mediated hepatic uptake is proven to be the rate-determining step in the systemic clearance of several drugs. Therefore, accurate measurement of active and passive uptake clearances
in ...vitro
is critical to facilitate pharmacokinetics and drug-drug interaction predictions. Here, we evaluated the plated human hepatocytes (PHH) and studied the effect of incubation temperature and inhibitor concentration on uptake measurements, in order to reliably estimate hepatic uptake components. Uptake rates measured using PHH, at 37°C without and with rifamycin SV, were comparable with those obtained from suspension hepatocytes and sandwich-cultured hepatocytes for a set of 10–13 compounds. Apparent permeability across monolayers of low-efflux Madin-Darby canine kidney cells was measured at 4, 10, and 37°C. Of the 23 compounds evaluated, 13 compounds showed >2-fold reduction in passive permeability at 4°C compared to 37°C, inferring that low-temperature incubations may underestimate passive uptake. Inhibition studies using transporter-transfected cells suggested that ∼20 μM rifamycin SV completely inhibited organic anion-transporting polypeptides (OATPs), while no significant inhibition was noted for other hepatic uptake transporters. On the basis of inhibition profiles, the contribution of active
versus
passive and OATP
versus
non-OATP transport to the PHH uptake was discerned for various endogenous substrates and statins. With the exception of fluvastatin, the statins studied were predominantly transported by OATPs in PHH and the non-OATP transporters, such as Na
+
-taurocholate co-transporting polypeptide, played a minimal role. In conclusion, PHH is useful for uptake measurements, and rifamycin SV employed at different concentrations can reliably estimate active and passive uptake and characterize OATP-dependent active uptake.
Human hepatocytes show marked differences in cell size, gene expression, and function throughout the liver lobules, an arrangement termed liver zonation. However, it is not clear if these zonal size ...differences, and the associated phenotypic differences, are retained in isolated human hepatocytes, the “gold standard” for in vitro studies of human liver function. Here, we therefore explored size differences among isolated human hepatocytes and investigated whether separation by size can be used to study liver zonation in vitro. We used counterflow centrifugal elutriation to separate cells into different size fractions and analyzed them with label‐free quantitative proteomics, which revealed an enrichment of 151 and 758 proteins (out of 5163) in small and large hepatocytes, respectively. Further analysis showed that protein abundances in different hepatocyte size fractions recapitulated the in vivo expression patterns of previously described zonal markers and biological processes. We also found that the expression of zone‐specific cytochrome P450 enzymes correlated with their metabolic activity in the different fractions. In summary, our results show that differences in hepatocyte size matches zonal expression patterns, and that our size fractionation approach can be used to study zone‐specific liver functions in vitro.
In this study, we used an elutriation‐based approach to study human liver zonation. We separated human hepatocytes into six size fractions and performed proteomic analysis, which showed that proteins that were enriched in different size fractions largely represented biological processes of known zonal specificity. Together with the strong correlations that we observed between expression and metabolic function of CYP enzymes with known zonation in vivo, we conclude that elutriation‐based size fractionation enables studies of zonated human liver functions in vitro.
Organic anion transporter (OAT) 2 and OAT7 were recently shown to be involved in the hepatic uptake of drugs; however, there is limited understanding of the population variability in the expression ...of these transporters in liver. There is also a need to derive relative expression-based scaling factors (REFs) that can be used to bridge in vitro functional data to the in vivo drug disposition. To this end, we quantified OAT2 and OAT7 surrogate peptide abundance in a large number of human liver tissue samples (n = 52), as well as several single-donor cryopreserved human hepatocyte lots (n = 30) by a novel, validated liquid chromatography tandem mass spectrometry (LC–MS/MS) method. The average surrogate peptide expression of OAT2 and OAT7 in the liver samples was 1.52 ± 0.57 and 4.63 ± 1.58 fmol/μg membrane protein, respectively. While we noted statistically significant differences (p < 0.05) in hepatocyte and liver tissue abundances for both OAT2 and OAT7, the differences were relatively small (1.8- and 1.5-fold difference in median values, respectively). Large interindividual variability was noted in the hepatic expression of OAT2 (16-fold in liver tissue and 23-fold in hepatocytes). OAT7, on the other hand, showed less interindividual variability (4-fold) in the livers, but high variability for the hepatocyte lots (27-fold). A significant positive correlation in OAT2 and OAT7 expression was observed, but expression levels were neither associated with age nor sex. In conclusion, our data suggest marked interindividual variability in the hepatic expression of OAT2/7, which may contribute to the pharmacokinetic variability of their substrates. Because both transporters were less abundant in hepatocytes than livers, a REF-based approach is recommended when scaling in vitro hepatocyte transport data to predict hepatic drug clearance and liver exposure of OAT2/7 substrates.
Targeted protein quantification using liquid chromatography-tandem mass spectrometry with stable isotope-labeled standards is recognized as the gold standard of practice for protein quantification. ...Such assays, however, can only cover a limited number of proteins, and developing targeted methods for larger numbers of proteins requires substantial investment. Alternatively, large-scale global proteomic experiments along with computational methods such as the "total protein approach" (TPA) have the potential to provide extensive protein quantification. In this study, we compared the TPA-based quantitation of seven major hepatic uptake transporters in four human liver tissue samples using global proteomic data obtained from two multiplexed tandem mass tag experiments (performed in two independent laboratories) to the quantitative data from targeted proteomic assays. The TPA-based quantitation of these hepatic transporters sodium-taurocholate cotransporting polypeptide (NTCP/SLC10A1), organic anion transporter 2 (OAT2/SLC22A7), OAT7/SLC22A9, organic anion-transporting polypeptide 1B1 (OATP1B1/SLCO1B1), OATP1B3/SLCO1B3, OATP2B1/SLCO2B1, and organic cation transporter (OCT1/SLC22A1) showed good-to-excellent correlations (Pearson
= 0.74-1.00) to the targeted data. In addition, the values were similar to those measured by targeted proteomics with 71% and 86% of the data sets falling within 3-fold of the targeted data. A comparison of the TPA-based quantifications of enzyme abundances to available literature data showed that the majority of the enzyme quantifications fell within the reference data intervals. In conclusion, these results demonstrate the capability of multiplexed global proteomic experiments to detect differences in protein expression between samples and provide reasonable estimations of protein expression levels.
Physiologically‐based pharmacokinetic (PBPK) models have an important role in drug discovery/development and decision making in regulatory submissions. This is facilitated by predefined PBPK ...platforms with user‐friendly graphical interface, such as Simcyp and PK‐Sim. However, evaluations of platform differences and the potential implications for disposition‐related applications are still lacking. The aim of this study was to assess how PBPK model development, input parameters, and model output are affected by the selection of PBPK platform. This is exemplified via the establishment of simvastatin PBPK models (workflow, final models, and output) in PK‐Sim and Simcyp as representatives of established whole‐body PBPK platforms. The major finding was that the choice of PBPK platform influenced the model development strategy and the final model input parameters, however, the predictive performance of the simvastatin models was still comparable between the platforms. The main differences between the structure and implementation of Simcyp and PK‐Sim were found in the absorption and distribution models. Both platforms predicted equally well the observed simvastatin (lactone and acid) pharmacokinetics (20–80 mg), BCRP and OATP1B1 drug–gene interactions (DGIs), and drug–drug interactions (DDIs) when co‐administered with CYP3A4 and OATP1B1 inhibitors/inducers. This study illustrates that in‐depth knowledge of established PBPK platforms is needed to enable an assessment of the consequences of PBPK platform selection. Specifically, this work provides insights on software differences and potential implications when bridging PBPK knowledge between Simcyp and PK‐Sim users. Finally, it provides a simvastatin model implemented in both platforms for risk assessment of metabolism‐ and transporter‐mediated DGIs and DDIs.