1. Introduction
The most direct way to analyse the human body it’s from its origin. Advances in various fields (e.g., proteomics, imaging, scanners) make it possible to study the human body, but there is too much “noise” to understand the data from observations [
1,
2,
3]. In contrast, in in vitro organ models make it easier to understand the underlying trends because there are considerably fewer variables, which also makes the models incomplete [
4,
5]. Each model has its advantages and disadvantages, and considering the purpose of the study, researchers select the most appropriate model for their studies [
6]. In the beginning, the models were oversimplified by using a single cell line, and the models consistently had paradoxical results [
7,
8,
9,
10]. However, in recent years, more comprehensive preclinical models have been developed [
11,
12]. Therefore, we can consider organ models as the first step toward a better understanding of the processes/mechanisms in the human body.
Assuming that an organ model recreates the in vivo microenvironment accurately, any toxicity assessment in such model should be accurate. Under these assumptions, numerous new organ models have been developed to test the safety of drugs [
13,
14,
15]. Two dimensional hepatocyte models have the advantage that there are relevant and primary human hepatocytes are recognised to be the ‘gold standard’ for in vitro hepatotoxic assays [
13]. Yet, they have a short lifespan, and they lack other cell types present in the liver that are essential for the liver normal functioning [
16]. The liver organ-on-a-chip (Ooc) models showed promising results [
15], but they are designed with a limited number of cells and low cell density in comparison to the in vivo environment [
16]. The spheroid liver models have a greater cell number than the Ooc models and resemble most closely the in vivo liver microenvironment [
14]. But the cells in the middle of the spheroid lack nourishment and the cell density is lower than the in vivo microenvironment [
16]. Thus, there are many different liver in vitro models, but none of the existing models had the characteristics that we considered critical for toxicity testing.
In this sense, we designed an in vitro liver model based on the characteristics previously proposed by Madorran et al. [
16]. In addition, we used preclinical and clinical instruments to measure the physiology and anatomy of the built model to favour the comparison between the values we measured in the model with the clinical data from the literature (clinical cases). Under this premises, we were able to build a cell-based model that have similar values to that observe in clinical cases, therefore, building a liver model that recreated the in vivo liver microenvironment. At the same time, we uncovered promising targets that may be used as liver function early markers in clinical diagnostics.
4. Discussion
Our goal was to build an in vitro liver model to resemble the liver microenvironment and allow clinical marker usage. To that end, we used primary liver cells: HSC, LSEC, Kupffer, and hepatocytes, because they reportedly retain the characteristics of the original cells, unlike cancer cells or immortalized cell lines [
6]. We also used different combinations of cell medium volume and formulation, ECM addition, and cell type proportion to configure a liver model. In this sense, we evaluate the influence of each feature in the liver model.
The influence of cell culture media
We observed that cells cultured with a medium rich on supplements (Hep medium) increased the viability of the cells (
Figure 3a), which is concordant with previous studies [
18,
19]. The addition of insulin to the medium increased the glucose uptake by the cells (
Figure 4), as described in the literature [
20,
21,
22]. In contrast, the net release of triglycerides observed in cells cultured with the same media may be due to the higher amino acid content of the media [
23]. The opposing dynamics of glucose and triglycerides (
Figure 4a) are consistent with in vivo observations [
23]. The addition of amino acids in the Hep medium increased the albumin synthesis significantly (
Figure 4a), which concurs with the literature [
24].
The influence of cell culture volume
In the initial phase, the higher cell density and interaction in samples cultured with 65 µl of cell medium (due to the cell arrangement seen in
Figure 2a) may increase cell proliferation [
25]. Yet, later necrosis of the cells may occur due to nutrient deficiency [
6,
26]. In this sense, and considering the lobule arrangement of the liver [
27], the cells cultured with 85 µl more closely resemble the in vivo situation (
Figure 2b), scattered colonies resembling lobule-like arrangement). There are many different types of stimulus that induce cell migration [
28]. But, considering that the different migrations we observed were related to cell density, we associate the migration to chemotaxis [
29]. This complex process needs to be studied in much more detail to better understand the observed migrations.
Furthermore, the clinical instrumentation is designed to measure the markers at certain range of concentration. However, the concentration of the molecules of interest is different in in vitro models and in vivo environment. Thus, lowering the volume to increase the cell number to cell volume ratio enabled the use of clinical instrumentation.
The influence of ECM
The addition of ECM did not affect the morphology of cells in the model (
Figure 2c). On the contrary, the addition of ECM had a noticeable impact on the viability of the cells. ECM induced the proliferation of hepatocytes (as previously observed by Wang Y. et al. [
30]) but decreased the viability of HSC cells. The designed model had HSC cells responsible for synthesising ECM in the liver [
31]. Thus, ECM production by the existing HSCs was sufficient [
31], and the additional ECM coating did not affect any of the biomarker values.
The influence of the proportion of the cell types
Higher hepatocyte density (80%) induced cell proliferation (
Figure 3), which was previously observed [
25]. This samples (80 % hepatocyte density) had less variable glucose content (
Figure 4a), as expected since hepatocytes are the liver cells that regulate glucose levels [
23]. The samples with higher hepatocyte density also had higher albumin content (
Figure 4a), as albumin is only synthesised in hepatocytes [
32].
After evaluating the influence of each feature in the model, we selected the feature combination that most faithfully recreated the liver microenvironment and expose the liver model to 50 µmol/l RIF, 1 mmol/l IBU or 500 µmol/l 5-FU. The concentrations chosen are based on previous studies by other authors [
33,
34,
35], analysing cell proliferation, IC50, EC50, gene expression, CYP activities and data from our own experiments.
Exposing the liver model to 5-FU
5-FU is a widely used chemotherapeutic drug and one of the most commonly utilized drugs for the treatment of various types of cancers because it inhibits thymidylate synthase [
36]. The hepatotoxicity of this drug is well documented [
34,
36,
37,
38], thus, is interesting as a toxic agent for liver models. The samples exposed to 5-FU had the lowest viability of all the evaluated samples (
Figure 6a) because only 5-FU was cytotoxic to hepatocytes (
Figure 6b). Since hepatocytes account for the vast majority of the population, any change in their cell number induces a higher effect in the total cell number of the sample. Yet, the highest cytotoxic effect was evaluated in Kupffer cells. In contrast, 5-FU induced the proliferation of LSEC cells during the first week, for which there are no previous data in literature). However, given the high cell death observed on day 5 (
Table 1), continuous exposure to 5-FU may also induce cell death in LSEC cells [
39]. Inhibition of thymidylate synthase and accumulation of toxic by-products from 5-FU catabolism (fluorocitrate for instance) may induce this delayed toxic effect [
40]. We also observed that adding 5-FU to the model decreased albumin synthesis. This is more likely due to the reduction of hepatocyte number [
41].
Exposing the liver model to IBU
IBU is a nonsteroidal anti-inflammatory drug and one of the most used drugs worldwide [
42]. There are many animal studies that observed the hepatotoxicity of this drug [
43,
44,
45,
46], and considering its exposure [
42], it is an interesting candidate for the evaluation of our model. The addition of IBU had an ambivalent effect on the cell viability of the samples. On the one hand, IBU significantly decreased the viability of Kupffer and HSC cells (
Figure 6b). On the other hand, samples with added IBU had twice as many LSEC cells as the untreated samples (
Table 1). The proliferation of LSEC was also observed in a previous study [
47], yet the mechanisms involved are unknown to date. As a result of both trends (cytotoxicity and proliferation), samples exposed to IBU had similar viability to the untreated samples, but the cell composition of the liver model changed (
Figure 6b). Yet, the most significant biological response induced by IBU was related to glucose dynamics. The samples treated with IBU had a significantly higher net glucose release than untreated samples. And these higher values were not the result of greater net triglyceride uptake. Thus, the influence of IBU on gluconeogenesis and glycogenolysis was also observed in vivo [
43,
48,
49]and could be the cause of this observation. Notably, there was no albumin content in samples treated with IBU, which may be caused due to the higher net glucose release (and lower availability for albumin synthesis) [
41].
Exposing the liver model to RIF
RIF is a widely used antibiotic to treat tuberculosis and other bacterial infections. Treatment with RIF is effective, but is known to induce drug-metabolizing enzymes in the liver [
50]. Thus, it is an interesting compound to evaluate the toxic assessment of our model. The addition of RIF significantly reduced the viability of Kupffer cells, like in the rest of the treated samples (
Figure 6b). On the contrary, RIF addition had a very distinct effect on the glucose and triglyceride dynamics. Both markers had a significantly higher net uptake by the cells compared to untreated samples (
Figure 7). This is in line with different studies reporting fatty acid accumulation upon RIF treatment [
51]. Moreover, this observation could also indicate an increase in FA accumulation within cells [
23]. This evaluation is consistent with previous studies in which upregulation of the free fatty acid transporter was observed in HepG2 cells treated with 10 µM of RIF [
51]. The albumin synthesis was significantly reduced when exposing the cells to RIF, maybe due to the above mentioned changes in the molecular pathways involved in glucose and triglyceride dynamics [
24]. Although further investigation is needed to corroborate the latter.
Evaluation of the liver model
The model has sufficient cell density to be statistically relevant, unlike Ooc models [
52] and it has a cell arrangement similar to that in spheroids, which reinforces the native physiology of the cells [
14]. It is visible under the microscope, allowing observation of detailed cell morphology. It can also be analysed with clinical instruments, which facilitates comparison with clinical data. In addition, the model enables comprehensive cell analysis with preclinical instruments. Most importantly, the observations in the model suggest that it resembles the liver microenvironment. The liver regulates the dynamics of triglycerides and glucose, which are in an inverse relationship [
23,
49] and these conditions are met in the model we developed (
Figure 4a). Furthermore, the addition of hepatotoxic drugs to the model resulted in similar physiological responses as
in vivo. Albumin synthesis is a key element in any liver model because albumin production is one of the liver’s main functions [
24]. In this sense, there are changes in albumin synthesis associated with various pathologies [
24,
53]. And in our model, we observed the influence on albumin synthesis by the hepatotoxic drugs (
Figure 7).
In the tested liver models, we did not observe any correlation between transaminase values and the viability of the cells. But the influence of transaminase levels in the liver is not clear in clinical studies either. High liver transaminases are found in patients whose liver is proliferating (after the liver recession) [
54]. But various studies have measured similar levels in patients with liver diseases [
55,
56,
57,
58,
59]. In addition, clinical studies have observed patients with abnormal levels who did not have liver disease [
58,
59]. These observations have led other authors to search for alternative markers in recent years [
60]. In this sense, this model may aid in the search for alternative markers, which is a rising concern in recent years [
61,
62,
63]. Moreover, this model may shed further light on some specific pathology related features, since it allows a more concrete analysis using preclinical evaluation methods besides the use of clinical instrumentation. To that end, our study may uncover an underlying trend that was unclear earlier. We observed that all hepatotoxic drugs had a greater impact on the viability of LSEC and Kupffer cells than on HSC and hepatocytes. Previous studies observed transcriptomic shifts between healthy and cirrhotic liver disease (CLD) scenarios [
64,
65]on LSEC cells. But they did not focus on the viability of both cell types. This finding are supported by the function of both cells, since LSEC and Kupffer cells are the liver cells responsible for xenobiotic uptake [
66]. Thus, LSEC and Kupffer cells are the liver cells first affected by exogenous agents [
66]. Therefore, we can use this model (or similar) to monitor molecules related to LSEC and Kupffer injury.
4.1. Limitations of the Model
In this experimental setup we observed the cells in the model with an inverted microscope, but we are working on the next experimental setup where we will observe the model with a confocal microscope. The use of a confocal microscope for further studies should allow us a more detailed study of the anatomy of the cells in this model.
The use of clinical instrumentation to evaluate the model was troublesome due to the scarcity of information on the use of clinical instruments in an in vitro liver model. In this sense, more replicates under similar conditions should be performed to increase confidence in the observed trends. Special focus should be drawn to molecules present in LSEC and Kupffer cells in the model when exposing it to toxic compounds. If these molecules were also present in the peripheral circulation, they may be used as markers of liver health.
With the current settings, the model is not a suitable for pharmacokinetic studies, but key changes in its configuration may adequate it to such studies. The introduction of a rocker or even perfusion is possible and would enable pharmacokinetic studies.
4.2. Highlights
Our in vitro liver model had similar biomarker levels to those observed in vivo, confirming the representativeness of the model.
The biological effects elicited by the selected drugs had similar trends to those observed in clinical studies.
This new approach facilitates the translation between basic and clinical research.
LSEC and Kupffer had a more pronounced biological effect (viability) than the other cell types when exposed to hepatotoxic drugs.
Figure 1.
Schematic presentation of the experimental protocol. a) Cells were cultured with different cell culture medium formulation and volume, ECM, and cell percentage of the different cell types. b) Table summarizing the samples with their feature combinations (each combination was cultured in parallel). c) All the combinations were analysed with the inverted microscope, the biochemical analyser and ISX. d) The feature combination that most resembled the liver microenvironment was selected, and e) was exposed to IBU, RIF and 5-FU. f) The same analysing methods were used to determine the toxic effect of the drugs in the liver model with the selected feature combination.
Figure 1.
Schematic presentation of the experimental protocol. a) Cells were cultured with different cell culture medium formulation and volume, ECM, and cell percentage of the different cell types. b) Table summarizing the samples with their feature combinations (each combination was cultured in parallel). c) All the combinations were analysed with the inverted microscope, the biochemical analyser and ISX. d) The feature combination that most resembled the liver microenvironment was selected, and e) was exposed to IBU, RIF and 5-FU. f) The same analysing methods were used to determine the toxic effect of the drugs in the liver model with the selected feature combination.
Figure 2.
Observations of cell arrangement and morphology with the inverted microscope. a) A sample cultured with 65 µL in which the cells have migrated towards the centre (5x magnification). b) A sample cultured with 85 µL in which the cells have formed scattered colonies (5x magnification). In both cases, a) and b), we observed the transition of the cells from day 1 to day 5 in a 96-well microplate. c) On the left, cells cultured with ECM and on the right, cells cultured without ECM (20x magnification).
Figure 2.
Observations of cell arrangement and morphology with the inverted microscope. a) A sample cultured with 65 µL in which the cells have migrated towards the centre (5x magnification). b) A sample cultured with 85 µL in which the cells have formed scattered colonies (5x magnification). In both cases, a) and b), we observed the transition of the cells from day 1 to day 5 in a 96-well microplate. c) On the left, cells cultured with ECM and on the right, cells cultured without ECM (20x magnification).
Figure 3.
Viability of cells concerning the different combinations of features (p-value > 0.05, < 0.05 *, <0.01 **, < 0.001 ***).
Figure 3.
Viability of cells concerning the different combinations of features (p-value > 0.05, < 0.05 *, <0.01 **, < 0.001 ***).
Figure 4.
Values of the liver markers in the various liver models built with different feature combinations. a) Liver marker concentrations. b) Liver marker values normalised to cell number (molecule biomarker per cell). (p-value > 0.05, < 0.05 *, <0.01 **, < 0.001 ***). We analysed alanine aminotransferase (ALT), aspartate aminotransferase (AST), gamma-glutamyl transferase (GGT), alkaline phosphatase (ALP), glucose, triglycerides, and albumin.
Figure 4.
Values of the liver markers in the various liver models built with different feature combinations. a) Liver marker concentrations. b) Liver marker values normalised to cell number (molecule biomarker per cell). (p-value > 0.05, < 0.05 *, <0.01 **, < 0.001 ***). We analysed alanine aminotransferase (ALT), aspartate aminotransferase (AST), gamma-glutamyl transferase (GGT), alkaline phosphatase (ALP), glucose, triglycerides, and albumin.
Figure 5.
Evaluation of cell morphology with the inverted microscope (20x magnification) on each cell treatment. a) Untreated, b) RIF, c) IBU, d) 5- FU.
Figure 5.
Evaluation of cell morphology with the inverted microscope (20x magnification) on each cell treatment. a) Untreated, b) RIF, c) IBU, d) 5- FU.
Figure 6.
Viability of samples treated with RIF, IBU, and 5-FU. a) Viability of samples relative to the initial seeding number. b) Viability of each cell type normalised to the untreated sample(p-value > 0.05, < 0.05 *, <0.01 **, < 0.001 ***). Rifampicin (RIF), ibuprofen (IBU) and 5-fluorouracil (5-FU), liver sinusoidal endothelial cells (LSEC), human stellate cells (HSC).
Figure 6.
Viability of samples treated with RIF, IBU, and 5-FU. a) Viability of samples relative to the initial seeding number. b) Viability of each cell type normalised to the untreated sample(p-value > 0.05, < 0.05 *, <0.01 **, < 0.001 ***). Rifampicin (RIF), ibuprofen (IBU) and 5-fluorouracil (5-FU), liver sinusoidal endothelial cells (LSEC), human stellate cells (HSC).
Figure 7.
Clinical markers of samples treated with RIF, IBU, and 5-FU. (p-value > 0.05, < 0.05 *, <0.01 **, < 0.001 ***). Rifampicin (RIF), ibuprofen (IBU) and 5-fluorouracil (5-FU), liver sinusoidal endothelial cells (LSEC), human stellate cells (HSC).
Figure 7.
Clinical markers of samples treated with RIF, IBU, and 5-FU. (p-value > 0.05, < 0.05 *, <0.01 **, < 0.001 ***). Rifampicin (RIF), ibuprofen (IBU) and 5-fluorouracil (5-FU), liver sinusoidal endothelial cells (LSEC), human stellate cells (HSC).
Table 1.
Cell death ratio and population share of each cell type 5 days after exposing the cell-based liver model to hepatotoxic drugs.
Table 1.
Cell death ratio and population share of each cell type 5 days after exposing the cell-based liver model to hepatotoxic drugs.
Group |
Share of the population (%) |
Cell death ratio (%) |
Share of the population (%) |
Cell death ratio (%) |
|
LSEC |
HSC |
Control |
4%±1 |
2%±0 |
5%±1 |
2%±1 |
5-FU |
6%±0 |
3%±0 |
5%±1 |
5%±1 |
IBU |
8%±1 |
2%±1 |
3%±1 |
5%±1 |
RIF |
5%±1 |
3%±1 |
5%±0 |
3%±0 |
|
Kupffer |
Hepatocytes |
Control |
4%±0 |
3%±1 |
88%±1 |
3%±1 |
5-FU |
2%±0 |
12%±0 |
87%±1 |
4%±0 |
IBU |
2%±0 |
12%±2 |
88%±1 |
2%±1 |
RIF |
3%±0 |
8%±0 |
88%±0 |
2%±0 |