Oxygen Consumption in a Hollow Fiber Bioartificial Liver-Revisited

Size-related variability of oxygen consumption rates in individual human hepatic cells

Accurate descriptions of the variability in single-cell oxygen consumption and its size-dependency are key to establishing more robust tissue models. By combining microfabricated devices with multiparameter identification algorithms, we demonstrate that single human hepatocytes exhibit an oxygen level-dependent consumption rate and that their maximal oxygen consumption rate is significantly lower than that of typical hepatic cell cultures. Moreover, we found that clusters of two or more cells competing for a limited oxygen supply reduced their maximal consumption rate, highlighting their ability to adapt to local resource availability and the presence of nearby cells. We used our approach to characterize the covariance of size and oxygen consumption rate within a cell population, showing that size matters, since oxygen metabolism covaries lognormally with cell size. Our study paves the way for linking the metabolic activity of single human hepatocytes to their tissue- or organ-level metabolism and describing its size-related variability through scaling laws.

A sense of proximity: Cell packing modulates oxygen consumption

Accurately modeling oxygen transport and consumption is crucial to predict metabolic dynamics in cell cultures and optimize the design of tissue and organ models. We present a methodology to characterize the Michaelis-Menten oxygen consumption parameters in vitro, integrating novel experimental techniques and computational tools. The parameters were derived for hepatic cell cultures with different dimensionality (i.e., 2D and 3D) and with different surface and volumetric densities. To quantify cell packing regardless of the dimensionality of cultures, we devised an image-based metric, referred to as the proximity index. The Michaelis-Menten parameters were related to the proximity index through an uptake coefficient, analogous to a diffusion constant, enabling the quantitative analysis of oxygen dynamics across dimensions. Our results show that Michaelis-Menten parameters are not constant for a given cell type but change with dimensionality and cell density. The maximum consumption rate per cell decreases significantly with cell surface and volumetric density, while the Michaelis-Menten constant tends to increase. In addition, the dependency of the uptake coefficient on the proximity index suggests that the oxygen consumption rate of hepatic cells is superadaptive, as they modulate their oxygen utilization according to its local availability and to the proximity of other cells. We describe, for the first time, how cells consume oxygen as a function of cell proximity, through a quantitative index, which combines cell density and dimensionality. This study enhances our understanding of how cell-cell interaction affects oxygen dynamics and enables better prediction of aerobic metabolism in tissue models, improving their translational value.

Closed-loop vasculature network design for bioprinting large, solid tissue scaffolds

Vascularization is an indispensable requirement for fabricating large solid tissues and organs. The natural vasculature derived from medical imaging modalities for large tissues and organs are highly complex and convoluted. However, the present bioprinting capabilities limit the fabrication of such complex natural vascular networks. Simplified bioprinted vascular networks, on the other hand, lack the capability to sustain large solid tissues. This work proposes a generalized and adaptable numerical model to design the vasculature by utilizing the tissue/organ anatomy. Starting with processing the patient's medical images, organ structure, tissue-specific cues, and key vasculature tethers are determined. An open-source abdomen magnetic resonance image dataset was used in this work. The extracted properties and cues are then used in a mathematical model for guiding the vascular network formation comprising arterial and venous networks. Next, the generated three-dimensional networks are used to simulate the nutrient transport and consumption within the organ over time and the regions deprived of the nutrients are identified. These regions provide cues to evolve and optimize the vasculature in an iterative manner to ensure the availability of the nutrient transport throughout the bioprinted scaffolds. The mass transport of six components of cell culture media-glucose, glycine, glutamine, riboflavin, human serum albumin, and oxygen was studied within the organ with designed vasculature. As the vascular structure underwent iterations, the organ regions deprived of these key components decreased significantly highlighting the increase in structural complexity and efficacy of the designed vasculature. The numerical method presented in this work offers a valuable tool for designing vascular scaffolds to guide the cell growth and maturation of the bioprinted tissues for faster regeneration post bioprinting.

A combined experimental and computational framework to evaluate the behavior of therapeutic cells for peripheral nerve regeneration

Recent studies have explored the potential of tissue‐mimetic scaffolds in encouraging nerve regeneration. One of the major determinants of the regenerative success of cellular nerve repair constructs (NRCs) is the local microenvironment, particularly native low oxygen conditions which can affect implanted cell survival and functional performance. In vivo, cells reside in a range of environmental conditions due to the spatial gradients of nutrient concentrations that are established. Here we evaluate in vitro the differences in cellular behavior that such conditions induce, including key biological features such as oxygen metabolism, glucose consumption, cell death, and vascular endothelial growth factor secretion. Experimental measurements are used to devise and parameterize a mathematical model that describes the behavior of the cells. The proposed model effectively describes the interactions between cells and their microenvironment and could in the future be extended, allowing researchers to compare the behavior of different therapeutic cells. Such a combinatorial approach could be used to accelerate the clinical translation of NRCs by identifying which critical design features should be optimized when fabricating engineered nerve repair conduits.

Establishment of physiologically relevant oxygen gradients in microfluidic organ chips

In vitro models of human organs must accurately reconstitute oxygen concentrations and gradients that are observed in vivo to mimic gene expression, metabolism, and host-microbiome interactions. Here we describe a simple strategy to achieve physiologically relevant oxygen tension in a two-channel human small intestine-on-a-chip (Intestine Chip) lined with primary human duodenal epithelium and intestinal microvascular endothelium in parallel channels separated by a porous membrane while both channels are perfused with oxygenated medium. This strategy was developed using computer simulations that predicted lowering the oxygen permeability of poly-dimethylsiloxane (PDMS) chips in specified locations using a gas impermeable film will allow the cells to naturally decrease the oxygen concentration through aerobic respiration and reach steady-state oxygen levels <36 mm Hg (<5%) within the epithelial lumen. The approach was experimentally confirmed using chips with embedded oxygen sensors that maintained this stable oxygen gradient. Furthermore, Intestine Chips cultured with this approach supported formation of a villus epithelium interfaced with a continuous endothelium and maintained intestinal barrier integrity for 72 h. This strategy recapitulates in vivo functionality in an efficient, inexpensive, and scalable format that improves the robustness and translatability of Organ Chip technology for studies on microbiome as well as oxygen sensitivity.

Scaling of joint mass and metabolism fluctuations in in silico cell-laden spheroids

Allometric scaling has many applications, from the prediction of pharmacokinetics in animals and humans to the probing of ecosystem dynamics. Most studies have neglected to account for variations and fluctuations, although they are intrinsic features of all biological systems. To understand how metabolic scaling emerges in the presence of variations, we developed computer-generated models of cell-laden spheroids to define the experimental size range of cell cultures with quantifiable similitudes in terms of fluctuations and metabolic scaling with living organisms. We show that the estimates of scaling exponents may change with increasing variability in both mass and metabolic rate. The computational pipeline described underpins the sound design of statistically meaningful cell-based models, with impacts in both biomedical science and ecology.

Variations and fluctuations are characteristic features of biological systems and are also manifested in cell cultures. Here, we describe a computational pipeline for identifying the range of three-dimensional (3D) cell-aggregate sizes in which nonisometric scaling emerges in the presence of joint mass and metabolic rate fluctuations. The 3D cell-laden spheroids with size and single-cell metabolic rates described by probability density functions were randomly generated in silico. The distributions of the resulting metabolic rates of the spheroids were computed by modeling oxygen diffusion and reaction. Then, a method for estimating scaling exponents of correlated variables through statistically significant data collapse of joint probability distributions was developed. The method was used to identify a physiologically relevant range of spheroid sizes, where both nonisometric scaling and a minimum oxygen concentration (0.04 mol⋅m⁻³) is maintained. The in silico pipeline described enables the prediction of the number of experiments needed for an acceptable collapse and, thus, a consistent estimate of scaling parameters. Using the pipeline, we also show that scaling exponents may be significantly different in the presence of joint mass and metabolic-rate variations typically found in cells. Our study highlights the importance of incorporating fluctuations and variability in size and metabolic rates when estimating scaling exponents. It also suggests the need for taking into account their covariations for better understanding and interpreting experimental observations both in vitro and in vivo and brings insights for the design of more predictive and physiologically relevant in vitro models.

A bioartificial liver support system integrated with a DLM/GelMA-based bioengineered whole liver for prevention of hepatic encephalopathy via enhanced ammonia reduction

Although bioartificial liver support systems (BLSSs) play an essential role in maintaining partial liver functions and detoxification for liver failure patients, hepatocytes are unanimously seeded in biomaterials, which lack the hierarchal structures and mechanical cues of native liver tissues. To address this challenge, we developed a new BLSS by combining a decellularized liver matrix (DLM)/GelMA-based bioengineered whole liver and a perfusion-based, oxygenated bioreactor. The novel bioengineered whole liver was fabricated by integrating photocrosslinkable gelatin (GelMA) and hepatocytes into a DLM. The combination of GelMA and the DLM not only provided a biomimetic extracellular microenvironment (ECM) for enhanced cell immobilization and growth with elevated hepatic functions (e.g., albumin secretion and CYP activities), but also provided biomechanical support to maintain the native structure of the liver. In addition, the perfusion-based, oxygenated bioreactor helped deliver oxygen to the interior tissues of the bioengineered liver, which was of importance for long-term culture. Most importantly, this new bioengineered whole liver decreased ammonia concentration by 45%, whereas direct seeding of hepatocytes in a naked DLM showed no significant reduction. Thus, the developed BLSS integrated with the DLM/GelMA-based bioengineered whole liver can potentially help elevate liver functions and prevent HE in liver failure patients while waiting for liver transplantation.

Determining the parameter space for effective oxygen depletion for FLASH radiation therapy

There has been a recent revival of interest in the FLASH effect, after experiments have shown normal tissue sparing capabilities of ultra-high-dose-rate radiation with no compromise on tumour growth restraint. A model has been developed to investigate the relative importance of a number of fundamental parameters considered to be involved in the oxygen depletion paradigm of induced radioresistance. An example eight-dimensional parameter space demonstrates the conditions under which radiation may induce sufficient depletion of oxygen for a diffusion-limited hypoxic cellular response. Initial results support experimental evidence that FLASH sparing is only achieved for dose rates on the order of tens of Gy/s or higher, for a sufficiently high dose, and only for tissue that is slightly hypoxic at the time of radiation. We show that the FLASH effect is the result of a number of biological, radiochemical and delivery parameters. Also, the threshold dose for a FLASH effect occurring would be more prominent when the parameterisation was optimised to produce the maximum effect. The model provides a framework for further FLASH-related investigation and experimental design. An understanding of the mechanistic interactions producing an optimised FLASH effect is essential for its translation into clinical practice.

Oxygen Consumption Characteristics in 3D Constructs Depend on Cell Density

Oxygen is not only crucial for cell survival but also a determinant for cell fate and function. However, the supply of oxygen and other nutrients as well as the removal of toxic waste products often limit cell viability in 3-dimensional (3D) engineered tissues. The aim of this study was to determine the oxygen consumption characteristics of 3D constructs as a function of their cell density. The oxygen concentration was measured at the base of hepatocyte laden constructs and a tightly controlled experimental and analytical framework was used to reduce the system geometry to a single coordinate and enable the precise identification of initial and boundary conditions. Then dynamic process modeling was used to fit the measured oxygen vs. time profiles to a reaction and diffusion model. We show that oxygen consumption rates are well-described by Michaelis-Menten kinetics. However, the reaction parameters are not literature constants but depend on the cell density. Moreover, the average cellular oxygen consumption rate (or OCR) also varies with density. We discuss why the OCR of cells is often misinterpreted and erroneously reported, particularly in the case of 3D tissues and scaffolds.

Large-scale high-density culture of hepatocytes in a liver microsystem with mimicked sinusoid blood flow

In vitro engineering of liver tissue is a rapidly developing field for various biomedical applications. However, liver tissue culture is currently performed on only a small scale with a low density of hepatocytes. In this study, a simple design was introduced in a liver microsystem to enhance the transport of nutrients (e.g., oxygen and glucose) for the three-dimensional large-scale, high-density culture of hepatocytes. In this design, convection across the cell culture zone was generated to mimic sinusoid blood flow (SBF) based on the pressure difference between two fluids flowing in a countercurrent manner on either side of the cell culture zone. First, the distributions of living and dead cells in different culture subzones under various perfusion flow rates were observed, analysed, and compared. Then, the enhanced transport of nutrients was experimentally validated in relation to the viability of cells and theoretically explained by comparing the fluid velocity and oxygen concentration distribution in the cell culture zone in counterflow and coflow modes. Finally, the functions of the SBF-mimicked liver microsystem were assessed on the basis of specific metabolites, synthesized proteins, and bilirubin detoxification of hepatocytes, with collagen and alginate as extracellular matrices. Under this design, the density of hepatocytes cultured at the 3-mm-thickness scale reached ~7 × 10⁷  cells/ml on Day 7, and the metabolism and detoxification functions of the cells worked well. In addition, a liver rope-like structure and sphere-like clusters of cells were observed. This work provides insight for the design of a bionic liver microsystem.

Simulation of blood and oxygen distributions in a hepatic lobule with sinusoids obstructed by cancer cells

The liver is one of the common metastatic sites for many cancers. The obstruction of sinusoids by circulating tumor cells changes liver microenvironments and is thus considered a source of hepatic metastases. To date, few studies provide detailed information, either experimentally or theoretically, concerning the changes in blood and oxygen distributions induced by the obstruction of sinusoids. In this study, we utilized a 3D porous medium-vascular tree geometric structure to mimic the hepatic lobule and studied theoretical blood flow and oxygen transport in the lobule. The simulation was validated with data from the literature. Then, the distributions of blood and oxygen in the presence of the obstruction by cancer cells were simulated. The area and degree of the liver damage induced by the obstruction were analyzed by comparing the difference of liver microenvironments between physiological (non-blocked sinusoid) and pathological (fully or partially blocked sinusoid) conditions and the minimum cancer cell sizes causing liver damage for various obstruction positions were obtained. The work presented in this study can be used to predict the degree of liver damage induced by the local ischemia caused by the obstruction of sinusoids and to characterize the relationship between hepatic metastases and liver microenvironments.

Human, Tissue-Engineered, Skeletal Muscle Myobundles to Measure Oxygen Uptake and Assess Mitochondrial Toxicity

Mitochondrial dysfunction is responsible for the toxicity of a number of drugs. Current isolated mitochondria or cellular monoculture mitochondrial respiration measurement systems lack physiological relevance. Using a tissue engineering rather than cell- or mitochondria-based approach enables a more physiologically relevant detection of drug-induced mitochondrial impairment. To probe oxygen consumption and mitochondrial health, we assayed the bioenergetic profile of engineered three-dimensional human skeletal muscle myobundles derived from primary myoblasts. Through experimental and computational techniques, we did not find external or internal oxygen transport limiting the engineered myobundles in the commercial O2k system to measure oxygen consumption. In response to the complex I inhibitor rotenone, myobundle basal respiration decreased dose dependently with an IC₅₀ of 9.24 ± 0.03 nM. At a 20 nM concentration of rotenone, myobundle maximal respiration decreased by 44.4% ± 9.8%. Respiratory depression by rotenone suggests that cultured myobundles rely heavily on the complex I pathway for ATP synthesis during times of both basal and increased energy demand. To address whether these decrements in mitochondrial function corresponded to alterations in physiological muscle function, we determined fatigue susceptibility that revealed a 46.0% ± 7.4% depression at 20 nM rotenone. The bioenergetic health index, which is a measure of normal oxidative mitochondrial function, was inversely correlated with the extent of fatigue. The human myobundles reproduce normal muscle metabolism under both basal and maximal energy demand conditions enabling the detection of drug-induced mitochondrial toxicity.

Allometric scaling in-vitro

About two decades ago, West and coworkers established a model which predicts that metabolic rate follows a three quarter power relationship with the mass of an organism, based on the premise that tissues are supplied nutrients through a fractal distribution network. Quarter power scaling is widely considered a universal law of biology and it is generally accepted that were in-vitro cultures to obey allometric metabolic scaling, they would have more predictive potential and could, for instance, provide a viable substitute for animals in research. This paper outlines a theoretical and computational framework for establishing quarter power scaling in three-dimensional spherical constructs in-vitro, starting where fractal distribution ends. Allometric scaling in non-vascular spherical tissue constructs was assessed using models of Michaelis Menten oxygen consumption and diffusion. The models demonstrate that physiological scaling is maintained when about 5 to 60% of the construct is exposed to oxygen concentrations less than the Michaelis Menten constant, with a significant concentration gradient in the sphere. The results have important implications for the design of downscaled in-vitro systems with physiological relevance.

Microengineered cell and tissue systems for drug screening and toxicology applications: Evolution of in-vitro liver technologies

The liver performs many key functions, the most prominent of which is serving as the metabolic hub of the body. For this reason, the liver is the focal point of many investigations aimed at understanding an organism's toxicological response to endogenous and exogenous challenges. Because so many drug failures have involved direct liver toxicity or other organ toxicity from liver generated metabolites, the pharmaceutical industry has constantly sought superior, predictive in-vitro models that can more quickly and efficiently identify problematic drug candidates before they incur major development costs, and certainly before they are released to the public. In this broad review, we present a survey and critical comparison of in-vitro liver technologies along a broad spectrum, but focus on the current renewed push to develop "organs-on-a-chip". One prominent set of conclusions from this review is that while a large body of recent work has steered the field towards an ever more comprehensive understanding of what is needed, the field remains in great need of several key advances, including establishment of standard characterization methods, enhanced technologies that mimic the in-vivo cellular environment, and better computational approaches to bridge the gap between the in-vitro and in-vivo results.

Hollow fibre membrane bioreactors for tissue engineering applications

Hollow fibre membrane bioreactors (HFB) provide a novel approach towards tissue engineering applications in the field of regenerative medicine. For adherent cell types, HFBs offer an in vivo-like microenvironment as each fibre replicates a blood capillary and the mass transfer rate across the wall is independent from the shear stresses experienced by the cell. HFB also possesses the highest surface area to volume ratio of all bioreactor configurations. In theory, these factors enable a high quantity of the desired cellular product with less population variation, and favourable operating costs. Experimental analyses of different cell types and bioreactor designs show encouraging steps towards producing a clinically relevant device. This review discusses the basic HFB design for cell expansion and in vitro models; compares data produced on commercially available systems and addresses the operational differences between theory and practice. HFBs are showing some potential for mammalian cell culture but further work is needed to fully understand the complexities of cell culture in HFBs and how best to achieve the high theoretical cell yields.

Effect of membranes on oxygen transfer rate and consumption within a newly developed three-compartment bioartificial liver device: Advanced experimental and theoretical studies

A mathematical model is developed to predict oxygen transfer in the fiber-in-fiber (FIF) bioartificial liver device. The model parameters are taken from the constructed and tested FIF modules. We extended the Krogh cylinder model by including one more zone for oxygen transfer. Cellular oxygen uptake was based on Michaelis-Menten kinetics. The effect of varying a number of important model parameters is investigated, including (1) oxygen partial pressure at the inlet, (2) the hydraulic permeability of compartment B (cell region), (3) the hydraulic permeability of the inner membrane, and (4) the oxygen diffusivity of the outer membrane. The mathematical model is validated by comparing its output against the experimentally acquired values of an oxygen transfer rate and the hydrostatic pressure drop. Three governing simultaneous linear differential equations are derived to predict and validate the experimental measurements, e.g., the flow rate and the hydrostatic pressure drop. The model output simulated the experimental measurements to a high degree of accuracy. The model predictions show that the cells in the annulus can be oxygenated well even at high cell density or at a low level of gas phase PG if the value of the oxygen diffusion coefficient Dm is 16 × 10(-5) . The mathematical model also shows that the performance of the FIF improves by increasing the permeability of polypropylene membrane (inner fiber). Moreover, the model predicted that 60% of plasma has access to the cells in the annulus within the first 10% of the FIF bioreactor axial length for a specific polypropylene membrane permeability and can reach 95% within the first 30% of its axial length.

Oxygen diffusion and consumption in extracellular matrix gels: implications for designing three-dimensional cultures

Three-dimensional (3D) cultures are increasingly used as tissue surrogates to study many physiopathological processes. However, to what extent current 3D culture protocols provide physiologic oxygen tension conditions remains ill defined. To address this limitation, oxygen tension was measured in a panel of acellular or cellularized extracellular matrix (ECM) gels with A549 cells, and analyzed in terms of oxygen diffusion and consumption. Gels included reconstituted basement membrane, fibrin and collagen. Oxygen diffusivity in acellular gels was up to 40% smaller than that of water, and the lower values were observed in the denser gels. In 3D cultures, physiologic oxygen tension was achieved after 2 days in dense (≥3 mg/mL) but not sparse gels, revealing that the latter gels are not suitable tissue surrogates in terms of oxygen distribution. In dense gels, we observed a dominant effect of ECM composition over density in oxygen consumption. All diffusion and consumption data were used in a simple model to estimate ranges for gel thickness, seeding density and time-window that may support physiologic oxygen tension. Thus, we identified critical variables for oxygen tension in ECM gels, and introduced a model to assess initial values of these variables, which may short-cut the optimization step of 3D culture studies.

Secretome analysis using a hollow fiber culture system for cancer biomarker discovery

Secreted proteins, collectively referred to as the secretome, were suggested as valuable biomarkers in disease diagnosis and prognosis. However, some secreted proteins from cell cultures are difficult to detect because of their intrinsically low abundance; they are frequently masked by the released proteins from lysed cells and the substantial amounts of serum proteins used in culture medium. The hollow fiber culture (HFC) system is a commercially available system composed of small fibers sealed in a cartridge shell; cells grow on the outside of the fiber. Recently, because this system can help cells grow at a high density, it has been developed and applied in a novel analytical platform for cell secretome collection in cancer biomarker discovery. This article focuses on the advantages of the HFC system, including the effectiveness of the system for collection of secretomes, and reviews the process of cell secretome collection by the HFC system and proteomic approaches to discover cancer biomarkers. The HFC system not only provides a high-density three-dimensional (3D) cell culture system to mimic tumor growth conditions in vivo but can also accommodate numerous cells in a small volume, allowing secreted proteins to be accumulated and concentrated. In addition, cell lysis rates can be greatly reduced, decreasing the amount of contamination by abundant cytosolic proteins from lysed cells. Therefore, the HFC system is useful for preparing a wide range of proteins from cell secretomes and provides an effective method for collecting higher amounts of secreted proteins from cancer cells. This article is part of a Special Issue entitled: An Updated Secretome.

Finite volume scheme for double convection-diffusion exchange of solutes in bicarbonate high-flux hollow-fiber dialyzer therapy

The efficiency of a high-flux dialyzer in terms of buffering and toxic solute removal largely depends on the ability to use convection-diffusion mechanism inside the membrane. A two-dimensional transient convection-diffusion model coupled with acid-base correction term was developed. A finite volume technique was used to discretize the model and to numerically simulate it using MATLAB software tool. We observed that small solute concentration gradients peaked and were large enough to activate solute diffusion process in the membrane. While CO(2) concentration gradients diminished from their maxima and shifted toward the end of the membrane, HCO(3)(-) concentration gradients peaked at the same position. Also, CO(2) concentration decreased rapidly within the first 47 minutes while optimal HCO(3)(-) concentration was achieved within 30 minutes of the therapy. Abnormally high diffusion fluxes were observed near the blood-membrane interface that increased diffusion driving force and enhanced the overall diffusive process. While convective flux dominated total flux during the dialysis session, there was a continuous interference between convection and diffusion fluxes that call for the need to seek minimal interference between these two mechanisms. This is critical for the effective design and operation of high-flux dialyzers.

A theoretical approach to zonation in a bioartificial liver

Bioartificial livers have yet to gain clinical acceptance. In a previous study, a theoretical model was utilized to create operating region charts that graphically illustrated viable bioartificial liver configurations. On this basis a rationale for the choice of operating and design parameters for the device was created. The concept is extended here to include aspects of liver zonation for further design optimization. In vivo, liver cells display heterogeneity with respect to metabolic activity according to their position in the liver lobule. It is thought that oxygen tension is a primary modulator of this heterogeneity and on this assumption a theoretical model to describe the metabolic zonation within an in vitro bioartificial liver device has been adopted. The distribution of the metabolic zones under varying design and operating parameters is examined. In addition, plasma flow rates are calculated that give rise to an equal distribution of the metabolic zones. The results show that when a clinically relevant number of cells are contained in the BAL (10 billion), it is possible to constrain each of the three metabolic zones to approximately one-third of the cell volume. This is the case for a number of different bioreactor designs. These considerations allow bioartificial liver design to be optimized.

Current development of bioreactors for extracorporeal bioartificial liver (Review)

The research and development of extracorporeal bioartificial liver is gaining pace in recent years with the introduction of a myriad of optimally designed bioreactors with the ability to maintain long-term viability and liver-specific functions of hepatocytes. The design considerations for bioartificial liver are not trivial; it needs to consider factors such as the types of cell to be cultured in the bioreactor, the bioreactor configuration, the magnitude of fluid-induced shear stress, nutrients' supply, and wastes' removal, and other relevant issues before the bioreactor is ready for testing. This review discusses the exciting development of bioartificial liver devices, particularly the various types of cell used in current reactor designs, the state-of-the-art culturing and cryopreservation techniques, and the comparison among many today's bioreactor configurations. This review will also discuss in depth the importance of maintaining optimal mass transfer of nutrients and oxygen partial pressure in the bioreactor system. Finally, this review will discuss the commercially available bioreactors that are currently undergoing preclinical and clinical trials.

A strategy to determine operating parameters in tissue engineering hollow fiber bioreactors

The development of tissue engineering hollow fiber bioreactors (HFB) requires the optimal design of the geometry and operation parameters of the system. This article provides a strategy for specifying operating conditions for the system based on mathematical models of oxygen delivery to the cell population. Analytical and numerical solutions of these models are developed based on Michaelis–Menten kinetics. Depending on the minimum oxygen concentration required to culture a functional cell population, together with the oxygen uptake kinetics, the strategy dictates the model needed to describe mass transport so that the operating conditions can be defined. If c_min ≫ K_m we capture oxygen uptake using zero-order kinetics and proceed analytically. This enables operating equations to be developed that allow the user to choose the medium flow rate, lumen length, and ECS depth to provide a prescribed value of c_min. When , we use numerical techniques to solve full Michaelis–Menten kinetics and present operating data for the bioreactor. The strategy presented utilizes both analytical and numerical approaches and can be applied to any cell type with known oxygen transport properties and uptake kinetics.

A low shear stress modular bioreactor for connected cell culture under high flow rates

A generic "system on a plate" modular multicompartmental bioreactor array which enables microwell protocols to be transferred directly to the bioreactor modules, without redesign of cell culture experiments or protocols is described. The modular bioreactors are simple to assemble and use and can be easily compared with standard controls since cell numbers and medium volumes are quite similar. Starting from fluid dynamic and mass transport considerations, a modular bioreactor chamber was first modeled and then fabricated using "milli-molding," a technique adapted from soft lithography. After confirming that the shear stress was extremely low in the system in the range of useful flow rates, the bioreactor chambers were tested using hepatocytes. The results show that the bioreactor chambers can increase or maintain cell viability and function when the flow rates are below 500 microL/min, corresponding to wall shear stresses of 10(-5) Pa or less at the cell culture surface.

Windows of operation for bioreactor design for the controlled formation of tissue-engineered arteries

The availability of large numbers of units of artificial arteries would offer significant benefits to the clinical management of bypass surgery. Tissue engineering offers the potential of providing vessels that can mimic the morphology, function, and physiological environment of native vessels. Ideally this would involve culturing stem cells in vitro within a biodegradable tubular scaffold so as to construct tissue for implantation. Essential to establishing a robust process for the production of tissue-engineered arteries is the understanding of the impact of changes in the operating conditions and bioreactor design on the construct formation. In this article, models of transport phenomena were developed to predict the critical flow rates and mass transfer requirements of a prototype bioreactor for the formation of tissue-engineered arteries. The impact of the cell concentration, tube geometry, oxygen effective diffusivity in alginate, substrate and metabolite concentration levels, feed rate, and recycle rate on the design of the bioreactor was visualized using windows of operation and contour plots. The result of this analysis determined the best configuration of the bioreactor that meets the cellular transport requirements as well as being reliable in performance while seeking to reduce the amount of nutrients to be used.

Transport advances in disposable bioreactors for liver tissue engineering

Acute liver failure (ALF) is a devastating diagnosis with an overall survival of approximately 60%. Liver transplantation is the therapy of choice for ALF patients but is limited by the scarce availability of donor organs. The prognosis of ALF patients may improve if essential liver functions are restored during liver failure by means of auxiliary methods because liver tissue has the capability to regenerate and heal. Bioartificial liver (BAL) approaches use liver tissue or cells to provide ALF patients with liver-specific metabolism and synthesis products necessary to relieve some of the symptoms and to promote liver tissue regeneration. The most promising BAL treatments are based on the culture of tissue engineered (TE) liver constructs, with mature liver cells or cells that may differentiate into hepatocytes to perform liver-specific functions, in disposable continuous-flow bioreactors. In fact, adult hepatocytes perform all essential liver functions. Clinical evaluations of the proposed BALs show that they are safe but have not clearly proven the efficacy of treatment as compared to standard supportive treatments. Ambiguous clinical results, the time loss of cellular activity during treatment, and the presence of a necrotic core in the cell compartment of many bioreactors suggest that improvement of transport of nutrients, and metabolic wastes and products to or from the cells in the bioreactor is critical for the development of therapeutically effective BALs. In this chapter, advanced strategies that have been proposed over to improve mass transport in the bioreactors at the core of a BAL for the treatment of ALF patients are reviewed.

Convection and hemoglobin-based oxygen carrier enhanced oxygen transport in a hepatic hollow fiber bioreactor

Hepatic hollow fiber bioreactors are a promising class of bioartificial liver assist device (BLAD). The development of this type of device is currently hindered by limited oxygen transport to cultured hepatocytes, due to low solubility of oxygen in aqueous media. In order to increase the oxygen spectrum to cultured hepatocytes housed within a hollow fiber bioreactor, several different engineering strategies were explored in this study. These included: supplementing the circulating media stream of the hollow fiber bioreactor with a hemoglobin-based oxygen carrier (bovine red blood cells) with defined oxygen binding and release kinetics and operating the bioreactor with media flow through the hollow fiber membrane into the extracapillary space (ECS). We hypothesize that these two strategies can be used to improve hepatocyte oxygenation and possibly attain an in vivo-like pO(2) spectrum, similar to that observed in vivo in the liver sinusoid. This work is significant, since provision of an in vivo-like pO(2) spectrum should create a fully functional BLAD that could potentially bridge thousands of liver failure patients towards native liver regeneration of damaged tissue or, if necessary, orthotopic liver transplantation.

Targeted Oxygen Delivery within Hepatic Hollow Fiber Bioreactors via Supplementation of Hemoglobin-Based Oxygen Carriers

Hepatic hollow fiber bioreactors are considered a promising class of bioartificial liver assist device (BLAD). Unfortunately, limited oxygen (O(2)) transport to hepatocytes within this device hinders further development. Hepatocytes in vivo (in the liver sinusoid) experience a wide range of oxygen tensions (pO(2) = 25-70 mmHg), which is important for development of proper differentiated function (zonation). Previously, we observed that bovine red blood cell (bRBC) supplementation of the circulating media stream enhanced oxygenation of cultured C3A hepatoma cells compared to a culture with no O(2) carrier (Gordon, J.; Palmer, A. F. Artif. Cells, BloodSubstitutes, Biotechnol. 2006, 33 (3), 297-306). Despite this success, the cells were not exposed to the desired in vivo O(2) spectrum (Sullivan, J.; Gordon, J.; Palmer, A. Biotechnol. Bioeng. 2006, 93 (2) 306-317). We hypothesize that altering the kinetics of O(2) binding/release to/from hemoglobin-based O(2) carriers (HBOCs) could potentially target O(2) delivery to cell cultures. High P(50) (low O(2) affinity) HBOCs preferentially targeted O(2) delivery at high inlet pO(2) values. Conversely, low P(50) (high O(2) affinity) HBOCs targeted O(2) delivery at low inlet pO(2) values. Additionally, inlet pO(2), flow rate, and HBOC concentration were varied to find optimal bioreactor operating conditions. Our results demonstrate that HBOCs can enhance O(2) delivery to cultured hepatocytes, while exposing them to in vivo-like O(2) tensions, which is critical to create a fully functional BLAD.

Enhanced oxygen delivery to primary hepatocytes within a hollow fiber bioreactor facilitated via hemoglobin-based oxygen carriers

The production of a fully functional bioartificial liver assist device (BLAD) would greatly enhance available treatment options for patients suffering from acute liver failure. Currently, inadequate oxygen provision to hepatocytes seeded within hollow fiber bioreactors hampers development of a viable hollow fiber-based BLAD. Experimentally, oxygen provision to primary rat hepatocytes cultured within hollow fiber bioreactors was measured, it was observed that supplementation with an oxygen carrier (bovine red blood cells at approximately 2% human hematocrit) did not significantly improve oxygenation compared to the absence of an oxygen carrier. Therefore, an oxygen transport model of an individual hollow fiber within the bioreactor was developed and simulated (up to approximately 10% human hematocrit) to more fully examine the effect of oxygen carrier supplementation on oxygenation within the bioreactor. The modeling analysis, supported via the experimental results, was utilized to predict optimal bioreactor operating conditions for the delivery of in vivo-like oxygen gradients to cultured hepatocytes in clinically relevant settings.

Mass transfer and metabolic reactions in hepatocyte spheroids cultured in rotating wall gas-permeable membrane system

Isolated hepatocytes in spheroid configuration exhibit a high degree of cell-cell contacts, which are important in the maintenance of viability and liver specific functions. In the absence of a vascular network, the cells in a large spheroid size experience mass transfer limitations of metabolites and oxygen in the core of aggregates. In this paper transport phenomena related to the diffusion and reaction of oxygen, glucose and lactate are mathematically described and experimentally verified for hepatocyte spheroids cultured in a rotating-wall polystyrene system (RWPS) not permeable for gases and in a rotating-wall membrane system (RWMS) with oxygen-permeable membrane. The concentration profiles of glucose, oxygen and lactate in the hepatocyte spheroids were estimated for different diameters of aggregates by solving the mass transfer equations for simultaneous diffusion and reaction, by finite element method. Simulation results evidenced that, for aggregates with size lower than 300 microm cultured in both RWPS and RWMS systems, the concentration profiles of glucose and lactate towards the core of spheroids (effective diffusion coefficients in the order of 10(-11)m(2)/s) are not significantly affected by the metabolic rate (c.a 10(-6)microg/mm(3)/s for glucose and about one order of magnitude less for lactate). On the contrary, the transport of oxygen (diffusion coefficient: 3.4 x 10(-10)m(2)/s, reaction rate: 1.5 x 10(-5)microg/mm(3)/s) is critically affected by the size of the multicellular spheroids and significant gradients in oxygen concentration may develop in spheroids. Aggregates with a size greater than 200 microm suffer severe oxygen limitation in the most part of its size attaining the lowest partial pressure in the centre. The improved viability predicted by the model culturing hepatocyte spheroids in the RWMS, characterized by a higher O(2) permeability with respect to RWPS, was experimentally confirmed. The results demonstrated that the mathematical model used in this study represents a useful support to experimental procedures in order to obtain hepatocyte spheroids with optimal size.

Hemoglobin-based O2 carrier O2 affinity and capillary inlet pO2 are important factors that influence O2 transport in a capillary

Hemopure (Biopure; Cambridge, MA) and PolyHeme (Northfield Laboratories; Evanston, IL) are two acellular hemoglobin-based O2 carriers (HBOCs) currently in phase III clinical trials for use as red blood cell substitutes. The most common adverse side effect that these HBOCs exhibit is increased vasoconstriction. Autoregulatory theory has been presented as a possible explanation for this physiological effect, where it is hypothesized that low-affinity HBOCs over-deliver O2 to tissues surrounding arterioles, thereby eliciting vasoconstriction. In this paper, we wanted to investigate HBOC oxygenation of tissue surrounding a capillary, which is the smallest element of the circulatory system. An a priori model has been developed in which the performance of mixtures of acellular HBOCs (synthesized by our group and others) and human red blood cells (hRBCs) has been simulated using a Krogh tissue cylinder model (KTCM) comprising a capillary surrounded by a capillary membrane and skeletal muscle tissue in cylindrical coordinates with specified tissue O2 consumption rates and Michaelis-Menten kinetics. In this study, the total hemoglobin (hRBCs and HBOCs) concentration was kept constant. The HBOCs studied possessed O2 affinities that were higher and lower compared to hRBCs (P50's spanned 5-55 mmHg), and the equilibrium binding/release of oxygen to/from the HBOCs was modeled using the Adair equation. At normoxic inlet pO2's, there was no correlation between O2 flux out of the capillary and the O2 affinity of the HBOC. However, a correlation was found between the average pO2 tension in the capillary and the O2 affinity of the HBOC. Additionally, we studied the change in the O2 equilibrium curve of HBOCs with different O2 affinities over a wide range of inlet pO2's and found that changing the inlet pO2 greatly affected which HBOC, having a unique O2 affinity, best delivered O2 to the surrounding tissue. The analysis of oxygen transport presented could lead to a better prediction of which acellular HBOC is best suited for a specific transfusion application that many times depends on the capillary inlet pO2 tension.

Oxygen uptake rates and liver-specific functions of hepatocyte and 3T3 fibroblast co-cultures

Bioartificial liver (BAL) devices have been developed to treat patients undergoing acute liver failure. One of the most important parameters to consider in designing these devices is the oxygen consumption rate of the seeded hepatocytes which are known to have oxygen consumption rates 10 times higher than most other cell types. Hepatocytes in various culture configurations have been tested in BAL devices including those formats that involve co-culture of hepatocytes with other cell types. In this study, we investigated, for the first time, oxygen uptake rates (OUR)s of hepatocytes co-cultured with 3T3-J2 fibroblasts at various hepatocyte to fibroblast seeding ratios. OURs were determined by measuring the rate of oxygen disappearance using a ruthenium-coated optical probe after closing and sealing the culture dish. Albumin and urea production rates were measured to assess hepatocyte function. Lower hepatocyte density co-cultures demonstrated significantly higher OURs (2 to 3.5-fold) and liver- specific functions (1.6-fold for albumin and 4.5-fold for urea production) on a per cell basis than those seeded at higher densities. Increases in OUR correlated well with increased liver-specific functions. OURs (V(m)) were modeled by fitting Michaelis-Menten kinetics and the model predictions closely correlated with the experimental data. This study provides useful information for predicting BAL design parameters that will avoid oxygen limitations, as well as maximize metabolic functions.

Metabolic rate does not scale with body mass in cultured mammalian cells

The allometric scaling of metabolic rate with organism body mass can be partially accounted for by differences in cellular metabolic rates. For example, hepatocytes isolated from horses consume almost 10-fold less oxygen per unit time as mouse hepatocytes [Porter and Brand, Am J Physiol Regul Integr Comp Physiol 269: R226–R228, 1995]. This could reflect a genetically programmed, species-specific, intrinsic metabolic rate set point, or simply the adaptation of individual cells to their particular in situ environment (i.e., within the organism). We studied cultured cell lines derived from 10 mammalian species with donor body masses ranging from 5 to 600,000 g to determine whether cells propagated in an identical environment (media) exhibited metabolic rate scaling. Neither metabolic rate nor the maximal activities of key enzymes of oxidative or anaerobic metabolism scaled significantly with donor body mass in cultured cells, indicating the absence of intrinsic, species-specific, cellular metabolic rate set points. Furthermore, we suggest that changes in the metabolic rates of isolated cells probably occur within 24 h and involve a reduction of cellular metabolism toward values observed in lower metabolic rate organisms. The rate of oxygen delivery has been proposed to limit cellular metabolic rates in larger organisms. To examine the effect of oxygen on steady-state cellular respiration rates, we grew cells under a variety of physiologically relevant oxygen regimens. Long-term exposure to higher medium oxygen levels increased respiration rates of all cells, consistent with the hypothesis that higher rates of oxygen delivery in smaller mammals might increase cellular metabolic rates.

Three-dimensional numerical modeling and computational fluid dynamics simulations to analyze and improve oxygen availability in the AMC bioartificial liver

A numerical model to investigate fluid flow and oxygen (O(2)) transport and consumption in the AMC-Bioartificial Liver (AMC-BAL) was developed and applied to two representative micro models of the AMC-BAL with two different gas capillary patterns, each combined with two proposed hepatocyte distributions. Parameter studies were performed on each configuration to gain insight in fluid flow, shear stress distribution and oxygen availability in the AMC-BAL. We assessed the function of the internal oxygenator, the effect of changes in hepatocyte oxygen consumption parameters in time and the effect of the change from an experimental to a clinical setting. In addition, different methodologies were studied to improve cellular oxygen availability, i.e. external oxygenation of culture medium, culture medium flow rate, culture gas oxygen content (pO(2)) and the number of oxygenation capillaries. Standard operating conditions did not adequately provide all hepatocytes in the AMC-BAL with sufficient oxygen to maintain O(2) consumption at minimally 90% of maximal uptake rate. Cellular oxygen availability was optimized by increasing the number of gas capillaries and pO(2) of the oxygenation gas by a factor two. Pressure drop over the AMC-BAL and maximal shear stresses were low and not considered to be harmful. This information can be used to increase cellular efficiency and may ultimately lead to a more productive AMC-BAL.

Evaluation of a bioartificial liver based on a nonwoven fabric bioreactor with porcine hepatocytes in pigs

BACKGROUND/AIMS

We developed a bioartificial liver (BAL) based on a direct hemoperfusion typed nonwoven fabric bioreactor containing porcine hepatocytes. In this study, the efficacy of our BAL was evaluated with a pig fulminant hepatic failure (FHF) model.

METHODS

FHF was induced with intravenous administration of D-galactosamine (1.3 g/kg) in each pig. Twelve hours post D-galactosamine injection, fifteen pigs were divided into: a BAL group (n = 5), in which pigs received the BAL treatment with 1.0 to 1.3 x 10(9) hepatocytes for 6 h, a sham BAL group (n = 5), in which pigs received the BAL treatment without hepatocytes, and a FHF group (n = 5), in which pigs only received intensive care. Parameters related to liver function and animal survival up to 168 h were determined.

RESULTS

In the BAL group, blood ammonia and plasma lactate levels were lower, and serum glucose levels and Fischer index were higher than those in the other two groups. Survival time of pigs in the BAL group was significantly prolonged as compared with the sham BAL and the FHF group.

CONCLUSIONS

The BAL based on a nonwoven fabric bioreactor containing porcine hepatocytes appears to be effective in the treatment of FHF in pigs.

Effects of oxygen transport on 3-d human mesenchymal stem cell metabolic activity in perfusion and static cultures: experiments and mathematical model

Human mesenchymal stem cells (hMSCs) have unique potential to develop into functional tissue constructs to replace a wide range of tissues damaged by disease or injury. While recent studies have highlighted the necessity for 3-D culture systems to facilitate the proper biological, physiological, and developmental processes of the cells, the effects of the physiological environment on the intrinsic tissue development characteristics in the 3-D scaffolds have not been fully investigated. In this study, experimental results from a 3-D perfusion bioreactor system and the static culture are combined with a mathematical model to assess the effects of oxygen transport on hMSC metabolism and proliferation in 3-D constructs grown in static and perfusion conditions. Cells grown in the perfusion culture had order of magnitude higher metabolic rates, and the perfusion culture supports higher cell density at the end of cultivation. The specific oxygen consumption rate for the constructs in the perfusion bioreactor was found to decrease from 0.012 to 0.0017 micromol/10(6) cells/h as cell density increases, suggesting intrinsic physiological change at high cell density. BrdU staining revealed the noneven spatial distribution of the proliferating cells in the constructs grown under static culture conditions compared to the cells that were grown in the perfusion system. The hypothesis that the constructs in static culture grow under oxygen limitation is supported by higher Y(L/G) in static culture. Modeling results show that the oxygen tension in the static culture is lower than that of the perfusion unit, where the cell density was 4 times higher. The experimental and modeling results show the dependence of cell metabolism and spatial growth patterns on the culture environment and highlight the need to optimize the culture parameters in hMSC tissue engineering.

Simulation of oxygen carrier mediated oxygen transport to C3A hepatoma cells housed within a hollow fiber bioreactor

A priori knowledge of the dissolved oxygen (O2) concentration profile within a hepatic hollow fiber (HF) bioreactor is important in developing an effective bioartificial liver assist device (BLAD). O2 provision is limiting within HF bioreactors and we hypothesize that supplementing a hepatic HF bioreactor's circulating media with bovine red blood cells (bRBCs), which function as an O2 carrier, will improve oxygenation. The dissolved O2 concentration profile within a single HF (lumen, membrane, and representative extra capillary space (ECS)) was modeled with the finite element method, and compared to experimentally measured data obtained on an actual HF bioreactor with the same dimensions housing C3A hepatoma cells. Our results (experimental and modeling) indicate bRBC supplementation of the circulating media leads to an increase in O2 consumed by C3A cells. Under certain experimental conditions (pO2,IN) = 95 mmHg, Q = 8.30 mL/min), the addition of bRBCs at 5% of the average in vivo human red blood cell concentration (% hRBC) results in approximately 50% increase in the O2 consumption rate (OCR). By simply adjusting the operating conditions (pO2,IN) = 25 mmHg, Q = 1.77 mL/min) and increasing bRBC concentration to 25% hRBC the OCR increase is approximately 10-fold. However, the improved O2 concentration profile experienced by the C3A cells could not duplicate the full range of in vivo O2 tensions (25-70 mmHg) typically experienced within the liver sinusoid with this particular HF bioreactor. Nonetheless, we demonstrate that the O2 transport model accurately predicts O2 consumption within a HF bioreactor, thus setting up the modeling framework for improving the design of future hepatic HF bioreactors.