Liver-on-a-chip devices: the pros and cons of complexity.

Physiologically relevant and broadly applicable liver cell culture platforms are of great importance in both drug development and disease modeling. Organ-on-a-chip systems offer a promising alternative to conventional, static two-dimensional (2-D) cultures, providing much-needed cues such as perfusion, shear stress, and three-dimensional (3-D) cell-cell communication. However, such devices cover a broad range of complexity both in manufacture and in implementation. In this review, we summarize the key features of the human liver that should be reflected in a physiologically relevant liver-on-a-chip model. We also discuss different material properties of importance in producing liver-on-a-chip devices and summarize recent and current progress in the field, highlighting different types of devices at different levels of complexity.

Intercellular communication induces glycolytic synchronization waves between individually oscillating cells.

Many organs have internal structures with spatially differentiated and sometimes temporally synchronized groups of cells. The mechanisms leading to such differentiation and coordination are not well understood. Here we design a diffusion-limited microfluidic system to mimic a multicellular organ structure with peripheral blood flow and test whether a group of individually oscillating yeast cells could form subpopulations of spatially differentiated and temporally synchronized cells. Upon substrate addition, the dynamic response at single-cell level shows glycolytic oscillations, leading to wave fronts traveling through the monolayered population and to synchronized communities at well-defined positions in the cell chamber. A detailed mechanistic model with the architectural structure of the flow chamber incorporated successfully predicts the spatial-temporal experimental data, and allows for a molecular understanding of the observed phenomena. The intricate interplay of intracellular biochemical reaction networks leading to the oscillations, combined with intercellular communication via metabolic intermediates and fluid dynamics of the reaction chamber, is responsible for the generation of the subpopulations of synchronized cells. This mechanism, as analyzed from the model simulations, is experimentally tested using different concentrations of cyanide stress solutions. The results are reproducible and stable, despite cellular heterogeneity, and the spontaneous community development is reminiscent of a zoned cell differentiation often observed in multicellular organs.

Studying Glycolytic Oscillations in Individual Yeast Cells by Combining Fluorescence Microscopy with Microfluidics and Optical Tweezers.

In this unit, we provide a clear exposition of the methodology employed to study dynamic responses in individual cells, using microfluidics for controlling and adjusting the cell environment, optical tweezers for precise cell positioning, and fluorescence microscopy for detecting intracellular responses. This unit focuses on the induction and study of glycolytic oscillations in single yeast cells, but the methodology can easily be adjusted to examine other biological questions and cell types. We present a step-by-step guide for fabrication of the microfluidic device, for alignment of the optical tweezers, for cell preparation, and for time-lapse imaging of glycolytic oscillations in single cells, including a discussion of common pitfalls. A user who follows the protocols should be able to detect clear metabolite time traces over the course of up to an hour that are indicative of dynamics on the second scale in individual cells during fast and reversible environmental adjustments. © 2018 by John Wiley & Sons, Inc.

Phosphofructokinase controls the acetaldehyde-induced phase shift in isolated yeast glycolytic oscillators.

The response of oscillatory systems to external perturbations is crucial for emergent properties such as synchronisation and phase locking and can be quantified in a phase response curve (PRC). In individual, oscillating yeast cells, we characterised experimentally the phase response of glycolytic oscillations for external acetaldehyde pulses and followed the transduction of the perturbation through the system. Subsequently, we analysed the control of the relevant system components in a detailed mechanistic model. The observed responses are interpreted in terms of the functional coupling and regulation in the reaction network. We find that our model quantitatively predicts the phase-dependent phase shift observed in the experimental data. The phase shift is in agreement with an adaptation leading to synchronisation with an external signal. Our model analysis establishes that phosphofructokinase plays a key role in the phase shift dynamics as shown in the PRC and adaptation time to external perturbations. Specific mechanism-based interventions, made possible through such analyses of detailed models, can improve upon standard trial and error methods, e.g. melatonin supplementation to overcome jet-lag, which are error-prone, specifically, since the effects are phase dependent and dose dependent. The models by Gustavsson and Goldbeter discussed in the text can be obtained from the JWS Online simulation database: (https://jjj.bio.vu.nl/models/gustavsson5 and https://jjj.bio.vu.nl/models/goldbeter1).

Applying Microfluidic Systems to Study Effects of Glucose at Single-Cell Level.

Microfluidic systems in combination with microscopy (e.g., fluorescence) can be a powerful tool to study, at single-cell level, the behavior and morphology of biological cells after uptake of glucose. Here, we briefly discuss the advantages of using microfluidic systems. We further describe how microfluidic systems are fabricated and how they are utilized. Finally, we discuss how the large amount of data can be analyzed in a "semi-automatic" manner using custom-made software. In summary, we provide a guide to how to use microfluidic systems in single-cell studies.

Single-cell study links metabolism with nutrient signaling and reveals sources of variability.

BACKGROUND: The yeast AMPK/SNF1 pathway is best known for its role in glucose de/repression. When glucose becomes limited, the Snf1 kinase is activated and phosphorylates the transcriptional repressor Mig1, which is then exported from the nucleus. The exact mechanism how the Snf1-Mig1 pathway is regulated is not entirely elucidated. RESULTS: Glucose uptake through the low affinity transporter Hxt1 results in nuclear accumulation of Mig1 in response to all glucose concentrations upshift, however with increasing glucose concentration the nuclear localization of Mig1 is more intense. Strains expressing Hxt7 display a constant response to all glucose concentration upshifts. We show that differences in amount of hexose transporter molecules in the cell could cause cell-to-cell variability in the Mig1-Snf1 system. We further apply mathematical modelling to our data, both general deterministic and a nonlinear mixed effect model. Our model suggests a presently unrecognized regulatory step of the Snf1-Mig1 pathway at the level of Mig1 dephosphorylation. Model predictions point to parameters involved in the transport of Mig1 in and out of the nucleus as a majorsource of cell to cell variability. CONCLUSIONS: With this modelling approach we have been able to suggest steps that contribute to the cell-to-cell variability. Our data indicate a close link between the glucose uptake rate, which determines the glycolytic rate, and the activity of the Snf1/Mig1 system. This study hence establishes a close relation between metabolism and signalling.

The yeast osmostress response is carbon source dependent.

Adaptation to altered osmotic conditions is a fundamental property of living cells and has been studied in detail in the yeast Saccharomyces cerevisiae. Yeast cells accumulate glycerol as compatible solute, controlled at different levels by the High Osmolarity Glycerol (HOG) response pathway. Up to now, essentially all osmostress studies in yeast have been performed with glucose as carbon and energy source, which is metabolised by glycolysis with glycerol as a by-product. Here we investigated the response of yeast to osmotic stress when yeast is respiring ethanol as carbon and energy source. Remarkably, yeast cells do not accumulate glycerol under these conditions and it appears that trehalose may partly take over the role as compatible solute. The HOG pathway is activated in very much the same way as during growth on glucose and is also required for osmotic adaptation. Slower volume recovery was observed in ethanol-grown cells as compared to glucose-grown cells. Dependence on key regulators as well as the global gene expression profile were similar in many ways to those previously observed in glucose-grown cells. However, there are indications that cells re-arrange redox-metabolism when respiration is hampered under osmostress, a feature that could not be observed in glucose-grown cells.

Design and fabrication of a scalable liver-lobule-on-a-chip microphysiological platform.

The design and fabrication of a very large-scale liver-lobule (VLSLL)-on-a-chip device, providing a microphysiological niche for hepatocytes, is described. The device consists of an integrated network of liver-lobule-like hexagonal tissue-culture chambers constructed in a hybrid layout with a separate seed-feed network. As a key feature, each chamber contains a central outlet mimicking the central vein of a liver lobule. Separating chamber walls located between the culture area and feed network protects cells from the shear force of the convective flow. Arrays of designated passages convey nutrients to the cells by diffusion-dominated mass transport. We simulated the flow velocity, shear stress and diffusion of glucose molecules inside and outside the culture chambers under a continuous flow rate of 1 µl min-1. As proof of concept, human hepatocellular carcinoma cells (HepG2) were cultured for periods of 5 and 14 days and human-induced pluripotent stem cell (hiPSC)-derived hepatocytes for 21 days. Stabilized albumin secretion and urea synthesis were observed in the microfluidic devices and cells maintained morphology and functionality during the culture period. Furthermore, we observed 3D tissue-like structure and bile-canaliculi network formation in the chips. Future applications of the described platform include drug development and toxicity studies, as well as the modeling of patient-specific liver diseases, and integration in multi-organ human-on-a-chip systems.

A Nonlinear Mixed Effects Approach for Modeling the Cell-To-Cell Variability of Mig1 Dynamics in Yeast.

The last decade has seen a rapid development of experimental techniques that allow data collection from individual cells. These techniques have enabled the discovery and characterization of variability within a population of genetically identical cells. Nonlinear mixed effects (NLME) modeling is an established framework for studying variability between individuals in a population, frequently used in pharmacokinetics and pharmacodynamics, but its potential for studies of cell-to-cell variability in molecular cell biology is yet to be exploited. Here we take advantage of this novel application of NLME modeling to study cell-to-cell variability in the dynamic behavior of the yeast transcription repressor Mig1. In particular, we investigate a recently discovered phenomenon where Mig1 during a short and transient period exits the nucleus when cells experience a shift from high to intermediate levels of extracellular glucose. A phenomenological model based on ordinary differential equations describing the transient dynamics of nuclear Mig1 is introduced, and according to the NLME methodology the parameters of this model are in turn modeled by a multivariate probability distribution. Using time-lapse microscopy data from nearly 200 cells, we estimate this parameter distribution according to the approach of maximizing the population likelihood. Based on the estimated distribution, parameter values for individual cells are furthermore characterized and the resulting Mig1 dynamics are compared to the single cell times-series data. The proposed NLME framework is also compared to the intuitive but limited standard two-stage (STS) approach. We demonstrate that the latter may overestimate variabilities by up to almost five fold. Finally, Monte Carlo simulations of the inferred population model are used to predict the distribution of key characteristics of the Mig1 transient response. We find that with decreasing levels of post-shift glucose, the transient response of Mig1 tend to be faster, more extended, and displays an increased cell-to-cell variability.

Entrainment of heterogeneous glycolytic oscillations in single cells.

Cell signaling, gene expression, and metabolism are affected by cell-cell heterogeneity and random changes in the environment. The effects of such fluctuations on cell signaling and gene expression have recently been studied intensively using single-cell experiments. In metabolism heterogeneity may be particularly important because it may affect synchronisation of metabolic oscillations, an important example of cell-cell communication. This synchronisation is notoriously difficult to describe theoretically as the example of glycolytic oscillations shows: neither is the mechanism of glycolytic synchronisation understood nor the role of cell-cell heterogeneity. To pin down the mechanism and to assess its robustness and universality we have experimentally investigated the entrainment of glycolytic oscillations in individual yeast cells by periodic external perturbations. We find that oscillatory cells synchronise through phase shifts and that the mechanism is insensitive to cell heterogeneity (robustness) and similar for different types of external perturbations (universality).

Rapid 3D fluorescence imaging of individual optically trapped living immune cells.

We demonstrate an approach to rapidly characterize living suspension cells in 4 dimensions while they are immobilized and manipulated within optical traps. A single, high numerical aperture objective lens is used to separate the imaging plane from the trapping plane. This facilitates full control over the position and orientation of multiple trapped cells using a spatial light modulator, including directed motion and object rotation, while also allowing rapid 4D imaging. This system is particularly useful in the handling and investigation of the behavior of non-adherent immune cells. We demonstrate these capabilities by imaging and manipulating living, fluorescently stained Jurkat T cells.

Allosteric regulation of phosphofructokinase controls the emergence of glycolytic oscillations in isolated yeast cells.

UNLABELLED: Oscillations are widely distributed in nature and synchronization of oscillators has been described at the cellular level (e.g. heart cells) and at the population level (e.g. fireflies). Yeast glycolysis is the best known oscillatory system, although it has been studied almost exclusively at the population level (i.e. limited to observations of average behaviour in synchronized cultures). We studied individual yeast cells that were positioned with optical tweezers in a microfluidic chamber to determine the precise conditions for autonomous glycolytic oscillations. Hopf bifurcation points were determined experimentally in individual cells as a function of glucose and cyanide concentrations. The experiments were analyzed in a detailed mathematical model and could be interpreted in terms of an oscillatory manifold in a three-dimensional state-space; crossing the boundaries of the manifold coincides with the onset of oscillations and positioning along the longitudinal axis of the volume sets the period. The oscillatory manifold could be approximated by allosteric control values of phosphofructokinase for ATP and AMP. DATABASE: The mathematical models described here have been submitted to the JWS Online Cellular Systems Modelling Database and can be accessed at http://jjj.mib.ac.uk/webMathematica/UItester.jsp?modelName=gustavsson5. [Database section added 14 May 2014 after original online publication].

Yeast AMP-activated protein kinase monitors glucose concentration changes and absolute glucose levels.

Analysis of the time-dependent behavior of a signaling system can provide insight into its dynamic properties. We employed the nucleocytoplasmic shuttling of the transcriptional repressor Mig1 as readout to characterize Snf1-Mig1 dynamics in single yeast cells. Mig1 binds to promoters of target genes and mediates glucose repression. Mig1 is predominantly located in the nucleus when glucose is abundant. Upon glucose depletion, Mig1 is phosphorylated by the yeast AMP-activated kinase Snf1 and exported into the cytoplasm. We used a three-channel microfluidic device to establish a high degree of control over the glucose concentration exposed to cells. Following regimes of glucose up- and downshifts, we observed a very rapid response reaching a new steady state within less than 1 min, different glucose threshold concentrations depending on glucose up- or downshifts, a graded profile with increased cell-to-cell variation at threshold glucose concentrations, and biphasic behavior with a transient translocation of Mig1 upon the shift from high to intermediate glucose concentrations. Fluorescence loss in photobleaching and fluorescence recovery after photobleaching data demonstrate that Mig1 shuttles constantly between the nucleus and cytoplasm, although with different rates, depending on the presence of glucose. Taken together, our data suggest that the Snf1-Mig1 system has the ability to monitor glucose concentration changes as well as absolute glucose levels. The sensitivity over a wide range of glucose levels and different glucose concentration-dependent response profiles are likely determined by the close integration of signaling with the metabolism and may provide for a highly flexible and fast adaptation to an altered nutritional status.

Heterogeneity of glycolytic oscillatory behaviour in individual yeast cells.

There are many examples of oscillations in biological systems and one of the most investigated is glycolytic oscillations in yeast. These oscillations have been studied since the 1950s in dense, synchronized populations and in cell-free extracts, but it has for long been unknown whether a high cell density is a requirement for oscillations to be induced, or if individual cells can oscillate also in isolation without synchronization. Here we present an experimental method and a detailed kinetic model for studying glycolytic oscillations in individual, isolated yeast cells and compare them to previously reported studies of single-cell oscillations. The importance of single-cell studies of this phenomenon and relevant future research questions are also discussed.

Osmostress-induced cell volume loss delays yeast Hog1 signaling by limiting diffusion processes and by Hog1-specific effects.

Signal transmission progresses via a series of transient protein-protein interactions and protein movements, which require diffusion within a cell packed with different molecules. Yeast Hog1, the effector protein kinase of the High Osmolarity Glycerol pathway, translocates transiently from the cytosol to the nucleus during adaptation to high external osmolarity. We followed the dynamics of osmostress-induced cell volume loss and Hog1 nuclear accumulation upon exposure of cells to different NaCl concentrations. While Hog1 nuclear accumulation peaked within five minutes following mild osmotic shock it was delayed up to six-fold under severe stress. The timing of Hog1 nuclear accumulation correlated with the degree of cell volume loss and the cells capacity to recover. Also the nuclear translocation of Msn2, the transcription factor of the general stress response pathway, is delayed upon severe osmotic stress suggesting a general phenomenon. We show by direct measurements that the general diffusion rate of Hog1 in the cytoplasm as well as its rate of nuclear transport are dramatically reduced following severe volume reduction. However, neither Hog1 phosphorylation nor Msn2 nuclear translocation were as much delayed as Hog1 nuclear translocation. Our data provide direct evidence that signaling slows down during cell volume compression, probably as a consequence of molecular crowding. Hence one purpose of osmotic adaptation is to restore optimal diffusion rates for biochemical and cell biological processes. In addition, there may be mechanisms slowing down especially Hog1 nuclear translocation under severe stress in order to prioritize Hog1 cytosolic targets.

Calibration of spatial light modulators suffering from spatially varying phase response.

We present a method for converting the desired phase values of a hologram to the correct pixel addressing values of a spatial light modulator (SLM), taking into account detailed spatial variations in the phase response of the SLM. In addition to thickness variations in the liquid crystal layer of the SLM, we also show that these variations in phase response can be caused by a non-uniform electric drive scheme in the SLM or by local heating caused by the incident laser beam. We demonstrate that the use of a global look-up table (LUT), even in combination with a spatially varying scale factor, generally does not yield sufficiently accurate conversion for applications requiring highly controllable output fields, such as holographic optical trapping (HOT). We therefore propose a method where the pixel addressing values are given by a three-dimensional polynomial, with two of the variables being the (x, y)-positions of the pixels, and the third their desired phase values. The coefficients of the polynomial are determined by measuring the phase response in 8 × 8 sub-sections of the SLM surface; the degree of the polynomial is optimized so that the polynomial expression nearly replicates the measurement in the measurement points, while still showing a good interpolation behavior in between. The polynomial evaluation increases the total computation time for hologram generation by only a few percent. Compared to conventional phase conversion methods, for an SLM with varying phase response, we found that the proposed method increases the control of the trap intensities in HOT, and efficiently prevents the appearance of strong unwanted 0th order diffraction that commonly occurs in SLM systems.

Initiation of the transcriptional response to hyperosmotic shock correlates with the potential for volume recovery.

The control of activity and localization of transcription factors is critical for appropriate transcriptional responses. In eukaryotes, signal transduction components such as mitogen-activated protein kinase (MAPK) shuttle into the nucleus to activate transcription. It is not known in detail how different amounts of nuclear MAPK over time affect the transcriptional response. In the present study, we aimed to address this issue by studying the high osmolarity glycerol (HOG) system in Saccharomyces cerevisiae. We employed a conditional osmotic system, which changes the period of the MAPK Hog1 signal independent of the initial stress level. We determined the dynamics of the Hog1 nuclear localization and cell volume by single-cell analysis in well-controlled microfluidics systems and compared the responses with the global transcriptional output of cell populations. We discovered that the onset of the initial transcriptional response correlates with the potential of cells for rapid adaptation; cells that are capable of recovering quickly initiate the transcriptional responses immediately, whereas cells that require longer time to adapt also respond later. This is reflected by Hog1 nuclear localization, Hog1 promoter association and the transcriptional response, but not Hog1 phosphorylation, suggesting that a presently uncharacterized rapid adaptive mechanism precedes the Hog1 nuclear response. Furthermore, we found that the period of Hog1 nuclear residence affects the amplitude of the transcriptional response rather than the spectrum of responsive genes.

The classical progesterone receptor mediates the rapid reduction of fallopian tube ciliary beat frequency by progesterone.

BACKGROUND: The transport of gametes as well as the zygote is facilitated by motile cilia lining the inside of the fallopian tube. Progesterone reduces the ciliary beat frequency within 30 minutes in both cows and mice. This rapid reduction suggest the involvement of a non-genomic signaling mechanism, although it is not known which receptors that are involved. Here we investigated the possible involvement of the classical progesterone receptor in this process. METHOD: The ciliary beat frequency of mice fallopian tube was measured ex vivo using an inverted bright field microscope and a high speed camera. The effects of the agonists progesterone and promegestone and an antagonist, mifeprestone, were investigated in wildtype mice. The effect of progesterone was also investigated in mice lacking the classical progesterone receptor. RESULTS: Progesterone, as well as the more specific PR agonist promegestone, significantly reduced the CBF at concentrations of 10-100 nanomolar within 10-30 minutes. In the absence of progesterone, the PR antagonist mifeprestone had no effect on the ciliary beat frequency at a concentration of 1 micromolar. When ciliated cells were pre-incubated with 1 micromolar mifeprestone, addition of progesterone did not reduce the ciliary beat frequency. Accordingly, in ciliated cells from mice not expressing the classical progesterone receptor, exposure to 100 nanomolar progesterone did not reduce the ciliary beat frequency. CONCLUSIONS: This is the first study to provide comprehensive evidence that the classical progesterone receptor mediates the rapid reduction of the tubal ciliary beat frequency by progesterone.

Anti-Stokes fluorescence from endogenously formed protoporphyrin IX--implications for clinical multiphoton diagnostics.

Multiphoton imaging based on two-photon excitation is making its way into the clinics, particularly for skin cancer diagnostics. It has been suggested that endogenously formed protoporphyrin IX (PpIX) induced by aminolevulinic acid or methylaminolevulinate can be applied to improve tumor contrast, in connection to imaging of tissue autofluorescence. However, previous reports are limited to cell studies and data from tissue are scarce. No report shows conclusive evidence that endogenously formed PpIX increases tumor contrast when performing multiphoton imaging in the clinical situation. We here demonstrate by spectral analysis that two-photon excitation of endogenously formed PpIX does not provide additional contrast in superficial basal cell carcinomas. In fact, the PpIX signal is overshadowed by the autofluorescent background. The results show that PpIX should be excited at a wavelength giving rise to one-photon anti-Stokes fluorescence, to overcome the autofluorescent background. Thus, this study reports on a plausible method, which can be implemented for clinical investigations on endogenously formed PpIX using multiphoton microscopy.

Two-photon fluorescence correlation spectroscopy as a tool for measuring molecular diffusion within human skin.

There is a need for tools enabling quantitative imaging of biological tissue for pharmaceutical applications. In this study, two-photon fluorescence microscopy (TPM) has been combined with fluorescence correlation spectroscopy (FCS), demonstrating proof-of-principle providing quantitative data of fluorophore concentration and diffusion in human skin. Measurements were performed on excised skin exposed to either rhodamine B (RB) or rhodamine B isothiocyanate (RBITC), chosen based on their similarity in fluorescence yield and molecular weight, but difference in chemical reactivity. The measurements were performed at tissue depths in the range 0 and 20 µm, and the diffusion coefficients at skin depths 5 and 10 µm were found to be significantly different (P<0.05). Overall median values for the diffusion coefficients were found to be 4.0×10(-13) m(2)/s and 2.0×10(-13) m(2)/s for RB and RBITC, respectively. These values correspond to the diffusion of a hard sphere with a volume eight times larger for RBITC compared to RB. This indicates that the RBITC have bound to biomolecules in the skin, and the measured signal is obtained from the RBITC-biomolecule complexes, demonstrating the potential of the TPM-FCS method to track molecular interactions in an intricate biological matrix such as human skin.