Investigations of oxygen transfer into Penicillium chrysogenum pellets by microprobe measurements

A microcoaxial needle sensor with a tip diameter of ca. 0.7 mum was used as a microprobe to measure profiles of dissolved oxygen tension (DOT) within fixed pellets of Penicillium chrysogenum as a function of the DOT level around the pellet, in the presence and absence of bulk convective flow and turbulence. The investigations indicate that the oxygen transfer mechanism is complex. The results were interpreted by assuming the penetration convective flow into the entire pellet and penetration of turbulence into the outer range. A model was developed which was able to describe the measured DOT profiles very well. The model takes into account molecular and turbulent diffusion as well as convective flow as transfer mechanisms inside of the pellet. Structures of pellets used for microprobe measurements were evaluated by histological investigations. Considerable variations of mycelial density with radius within the pellets were found.     

A method to evaluate the isothermal effectiveness factor for dynamic oxygen into mycelial pellets in submerged cultures

Several models have been developed simulating O2 transfer in bioreactors, but three limitations are often found: (i) an inadequate kinetic representation of O2 consumption or wrong boundary conditions, (ii) unrealistic parameter values, and (iii) inadequate experimental systems. In our study we minimized those possible sources of error. Oxygen uptake rate, void fraction of the pellet, and external O2 mass transfer coefficient were experimentally obtained from bioreactor studies in which pellets of Gibberella fujikuroi were naturally formed. Michaelis-Menten kinetics and diffusion equations were used to describe the O2 consumption rate and to evaluate the effectiveness factor in dynamic mode. The nonlinear mathematical model proposed was solved by the orthogonal collocation technique. The O2 consumption rate in pellets of G. fujikuroi of 1.7-2.0 mm is only marginally inhibited by diffusion constraints under conditions tested. Simulation analysis showed that the effectiveness factor decreased as the Thiele modulus and pellet diameter increased. The proposed model was applied to experimental data reported for other fungal pellets and allowed to predict optimal conditions for O2 transfer into mycelial pellets.     

Comparative chondrogenesis of human cells in a 3D integrated experimental–computational mechanobiology model

We present an integrated experimental–computational mechanobiology model of chondrogenesis. The response of human articular chondrocytes to culture medium perfusion, versus perfusion associated with cyclic pressurisation, versus non-perfused culture, was compared in a pellet culture model, and multiphysic computation was used to quantify oxygen transport and flow dynamics in the various culture conditions. At 2 weeks of culture, the measured cell metabolic activity and the matrix content in collagen type II and aggrecan were greatest in the perfused+pressurised pellets. The main effects of perfusion alone, relative to static controls, were to suppress collagen type I and GAG contents, which were greatest in the non-perfused pellets. All pellets showed a peripheral layer of proliferating cells, which was thickest in the perfused pellets, and most pellets showed internal gradients in cell density and matrix composition. In perfused pellets, the computed lowest oxygen concentration was 0.075 mM (7.5% tension), the maximal oxygen flux was 477.5 nmol/m²/s and the maximal fluid shear stress, acting on the pellet surface, was 1.8 mPa (0.018 dyn/cm²). In the non-perfused pellets, the lowest oxygen concentration was 0.003 mM (0.3% tension) and the maximal oxygen flux was 102.4 nmol/m²/s. A local correlation was observed, between the gradients in pellet properties obtained from histology, and the oxygen fields calculated with multiphysic simulation. Our results show up-regulation of hyaline matrix protein production by human chondrocytes in response to perfusion associated with cyclic pressurisation. These results could be favourably exploited in tissue engineering applications.     

Oxygen transfer into mycelial pellets

The oxygen uptake rate in mycelial pellets of Aspergillus niger was studied experimentally and theoretically. The specific rate of respiration of mycelial pellets was found to decrease significantly with increasing pellet size. The distribution of respiratory activity in the mycelial pellets was evaluated and the specific rate of respiration of disrupted mycelia showed adaptation to the concentration of oxygen in the medium. The decrease of the specific rate of respiration of the mycelial pellets could be estimated according to diameter, mycelial density, oxygen diffusivity, and adaptation to the concentration of oxygen. Good agreement was found between the theoretical analysis and the experimental data.     

The effect of oxygen tension on human articular chondrocyte matrix synthesis: Integration of experimental and computational approaches

Significant oxygen gradients occur within tissue engineered cartilaginous constructs. Although oxygen tension is an important limiting parameter in the development of new cartilage matrix, its precise role in matrix formation by chondrocytes remains controversial, primarily due to discrepancies in the experimental setup applied in different studies. In this study, the specific effects of oxygen tension on the synthesis of cartilaginous matrix by human articular chondrocytes were studied using a combined experimental‐computational approach in a “scaffold‐free” 3D pellet culture model. Key parameters including cellular oxygen uptake rate were determined experimentally and used in conjunction with a mathematical model to estimate oxygen tension profiles in 21‐day cartilaginous pellets. A threshold oxygen tension (pO₂ ≈ 8% atmospheric pressure) for human articular chondrocytes was estimated from these inferred oxygen profiles and histological analysis of pellet sections. Human articular chondrocytes that experienced oxygen tension below this threshold demonstrated enhanced proteoglycan deposition. Conversely, oxygen tension higher than the threshold favored collagen synthesis. This study has demonstrated a close relationship between oxygen tension and matrix synthesis by human articular chondrocytes in a “scaffold‐free” 3D pellet culture model, providing valuable insight into the understanding and optimization of cartilage bioengineering approaches. Biotechnol. Bioeng. 2014;111: 1876–1885. © 2014 The Authors. Biotechnology and Bioengineering Published by Wiley Periodicals, Inc.     

Diffusion and reaction within porous packing media: a phenomenological model

A phenomenological model has been developed to describe biomass distribution and substrate depletion in porous diatomaceous earth (DE) pellets colonized by Pseudomonas aeruginosa. The essential features of the model are diffusion, attachment and detachment to/from pore walls of the biomass, diffusion of substrate within the pellet, and external mass transfer of both substrate and biomass in the bulk fluid of a packed bed containing the pellets. A bench-scale reactor filled with DE pellets was inoculated with P. aeruginosa and operated in plug flow without recycle using a feed containing glucose as the limiting nutrient. Steady-state effluent glucose concentrations were measured at various residence times, and biomass distribution within the pellet was measured at the lowest residence time. In the model, microorganism/substrate kinetics and mass transfer characteristics were predicted from the literature. Only the attachment and detachment parameters were treated as unknowns, and were determined by fitting biomass distribution data within the pellets to the mathematical model. The rate-limiting step in substrate conversion was determined to be internal mass transfer resistance; external mass transfer resistance and microbial kinetic limitations were found to be nearly negligible. Only the outer 5% of the pellets contributed to substrate conversion. (c) 1993 Wiley & Sons, Inc.     

A model for oxygen supply in fixed bed reactors with immobilized hybridoma cells

In fixed bed reactors with animal cells immobilized in macroporous carriers sufficient oxygen supply is a critical parameter. For modelling of the oxygen consumption and the oxygen profile in a fixed bed oxygen gradients within the porous carriers and along the length of the fixed bed have to be considered. For the complex geometry of the fixed bed a model structure was assumed, that allows the calculation of the oxygen profile. The model for oxygen supply of the immobilized cells included the transport resistance from the bulk fluid into the carriers and diffusion within the carriers. The model was compared with experimental data obtained with a hybridoma cell line for production of monoclonal antibodies. Model calculations and experimental data agree rather well. The mean volume-specific oxygen uptake rate as an indicator for the cell activity increased with the superficial flow velocity of the bulk liquid flow, and did not depend on the length of the fixed bed in the range tested. This indicates, that the convective transport from the bulk liquid flow between the carriers to the outer surface of the carriers is a dominating transport resistance besides the diffusive oxygen supply within the carriers.     

Nutrient utilization by bovine articular chondrocytes: a combined experimental and theoretical approach

A combined experimental-numerical approach was adopted to characterize glucose and oxygen uptake and lactate production by bovine articular chondrocytes in a model system. For a wide range of cell concentrations, cells in agarose were supplemented with either low or high glucose medium. During an initial culture phase of 48 h, oxygen was monitored noninvasively using a biosensor system. Glucose and lactate were determined by medium sampling. In order to quantify glucose and oxygen uptake, a finite element approach was adopted to describe diffusion and uptake in the experimental model. Numerical predictions of lactate, based on simple relations for cell metabolism, were found to agree well for low glucose, but not for high glucose medium. Oxygen did not play a role in either case. Given the close association between chondrocyte energy metabolism and matrix synthesis, a quantifiable prediction of utilization can present a valuable contribution in the optimization of tissue engineering conditions.     

Impact of pellet size on growth and lignin peroxidase activity of Phanerochaete chrysosporium

We investigated the influence of pellet size on the growth and lignin peroxidase (LiP) productivity of Phanerochaete chrysosporium. Different pellet sizes were obtained by varying the vessel diameter under constant shaking conditions. Under these varying conditions the pellet size was in the range of 2–18 mm, while the number of pellets in a single vessel varied from around 1,200 in the Erlenmeyer flask to around 6 in the narrowest vessel. A correlation between the final pellet size and the shear rate was obtained, demonstrating that the pellet size is mainly affected by hydrodynamics. The growth of large pellets was described by a cubic growth model. Despite different pellet sizes, LiP activity appeared in all vessels, but the onset of LiP activity showed a delay based upon the pellet size, while maximal LiP activities varied by only 15%, being around 850 U/l.     

Monoclonal antibody productivity and the metabolic pattern of perfusion cultures under varying oxygen tensions

The metabolic pattern and cell culture kinetics of high-cell-density perfusion cultures were compared under two different oxygen transfer conditions: oxygen limiting and not limiting. When oxygen was a limiting factor during perfusion culture, both specific glucose uptake and lactate production rates increased, compared to non-oxygen-limited condition, by about 60% and 30%, respectively. The specific glutamine uptake rate under oxygen-limited conditions was almost 4.0 times higher than that under non-oxygen-limited conditions. The activity of lactate dehydrogenase (LDH) released into the medium by the dead cells can be used as an indicator for the metabolic and physiological conditions related to oxygen limitation. There was a 3.2 times higher specific rate of LDH activity released by dead cells in oxygen-limited cultures than those in non-oxygen-limited cultures. The specific production rate of monoclonal antibody was not significantly affected by the oxygen transfer conditions during the rapid cell growth period, but it rapidly increased toward the end of perfusion cultures. The higher perfusion rate may have limited further cell growth during high-cell-density perfusion culture, because cell damage was caused by the hydrodynamic shear within a hollow fiber microfiltration cartridge installed to withdraw the spent medium and the waste metabolites. (c) 1993 John Wiley & Sons, Inc.     

Determination of the respiration kinetics for mycelial pellets of Phanerochaete chrysosporium.

In mycelial pellet cultures of the white rot basidiomycete Phanerochaete chrysosporium, low oxygen concentration negatively affects the production of the extracellular lignin peroxidases and manganese peroxidases which are key components of the lignin-degrading system of this organism. To test the hypothesis that oxygen limitation in the pellets is responsible for this effect, oxygen microelectrodes were used to determine oxygen concentration gradients within the mycelial pellets of P. chrysosporium. Pellets were removed from oxygenated cultures, allowed to equilibrate with air, and probed with oxygen microelectrodes. The oxygen profiles were modelled assuming that O2 uptake follows a Michaelis-Menten relationship. The Vmax and Km values for oxygen uptake were 0.76 +/- 0.10 g/m3 of pellet per s and 0.5 +/- 0.3 g/m3, respectively. These kinetic values were used to predict respiration rates in air-flushed cultures, oxygen-flushed cultures, and cultures with large pellets (diameter greater than 6 mm). The predicted respiration rates were independently validated by experimentally measuring the evolution of carbon dioxide from whole cultures.     

Effects of dissolved oxygen tension and mechanical forces on fungal morphology in submerged fermentation

The effects of dissolved oxygen tension and mechanical forces on fungal morphology were both studied in the submerged fermentation of Aspergillus awamori. Pellet size, the hairy length of pellets, and the free filamentous mycelial fraction in the total biomass were found to be a function of the mechanical force intensity and to be independent of the dissolved oxygen tension provided that the dissolved oxygen tension was neither too low (5%) nor too high (330%). When the dissolved oxygen concentration was close to the saturation concentration corresponding to pure oxygen gas, A. awamori formed denser pellets and the free filamentous mycelial fraction was almost zero for a power input of about 1 W/kg. In the case of very low dissolved oxygen tension, the pellets were rather weak and fluffy so that they showed a very different appearance. The amount of biomass per pellet surface area appeared to be affected only by the dissolved oxygen tension and was proportional to the average dissolved oxygen tension to the power of 0.33. From this it was concluded that molecular diffusion was the dominant mechanism for oxygen transfer in the pellets and that convection and turbulent flow in the pellets were negligible in submerged fermentations. The biomass per wet pellet volume increased with the dissolved oxygen tension and decreased with the size of the pellets. This means that the smaller pellets formed under a higher dissolved oxygen tension had a higher intrinsic strength. Correspondingly, the porosity of the pellets was a function of the dissolved oxygen tension and the size of pellets. Within the studied range, the void fraction in the pellets was high and always much more than 50%.     

Determination of the respiration kinetics for mycelial pellets of Phanerochaete chrysosporium.

In mycelial pellet cultures of the white rot basidiomycete Phanerochaete chrysosporium, low oxygen concentration negatively affects the production of the extracellular lignin peroxidases and manganese peroxidases which are key components of the lignin-degrading system of this organism. To test the hypothesis that oxygen limitation in the pellets is responsible for this effect, oxygen microelectrodes were used to determine oxygen concentration gradients within the mycelial pellets of P. chrysosporium. Pellets were removed from oxygenated cultures, allowed to equilibrate with air, and probed with oxygen microelectrodes. The oxygen profiles were modelled assuming that O2 uptake follows a Michaelis-Menten relationship. The Vmax and Km values for oxygen uptake were 0.76 +/- 0.10 g/m3 of pellet per s and 0.5 +/- 0.3 g/m3, respectively. These kinetic values were used to predict respiration rates in air-flushed cultures, oxygen-flushed cultures, and cultures with large pellets (diameter greater than 6 mm). The predicted respiration rates were independently validated by experimentally measuring the evolution of carbon dioxide from whole cultures.     

Adding talc microparticles to Aspergillus terreus ATCC 20542 preculture decreases fungal pellet size and improves lovastatin production

Changing fungal morphology with the use of morphological engineering techniques leads to improving the production of metabolites by filamentous fungi in the submerged culture. Adding mineral microparticles is one such simple method to change fungal pellet size. Here, it was studied for a lovastatin producer, Aspergillus terreus ATCC 20542. The experiments were conducted in shake flasks and 10 μm talc microparticles were added to the preculture. Intrapellet oxygen concentration profiles were determined by an oxygen microprobe. Talc microparticles caused a decrease of A. terreus pellets diameter from about 2000 to 900 μm, dependent on their concentration in the preculture. Smaller pellets produced more lovastatin, whose titre exceeded then 120 mg L^?1, utilising more lactose. The decrease in pellet size resulted in changes of oxygen concentration profiles in the pellets. The estimated critical pellet diameter, at which the non‐oxygenated zone was observed in the centre of the pellets, was 1700 μm. Smaller pellets were fully penetrated by oxygen. To conclude, facilitated diffusion of oxygen into the pellets of smaller diameter and their less dense structure made lactose utilisation by A. terreus more efficient, which ultimately increased lovastatin production in the runs with talc microparticles added, compared to the control runs.     

Studies of cell pellets: I. Electrical properties and porosity.

Cell pellets formed by centrifugation provided a good system to study the osmotic behavior, electroporation, and interaction between cells. Rabbit erythrocyte pellets were used in this study because they were simpler than nucleated cells to model analytically. Structurally, cell pellets possessed properties of porous solid bodies and gels. Electrically, cell pellets were shown to behave as a parallel set of resistance, Rp, and capacitance, Cp. Information on pellet structures was obtained from electric measurements. The pellet resistance reflected the intercellular conductivity (porosity and gap conductivity), whereas the pellet capacitance depended mostly on membrane capacitance. The pellet resistance was more sensitive to experimental conditions. The intercellular gap distance can be derived from pellet porosity measurements, providing the cell volume and surface area were known. Rp increased and relaxed exponentially with time when centrifugation started and stopped; the cycles were reversible. When supernatants were exchanged with solutions containing hypotonic electrolytes or macromolecules (such as PEG) after the pellets were formed, complicated responses to different colloidal osmotic effects were observed. A transient decrease followed by a large increase of Rp was observed after the application of a porating electric pulse, as expected from a momentary membrane breakdown, followed by a limited colloidal-osmotic swelling of pelleted cells. The equilibrium values of Rp, Cp, pellet porosity, and intercellular distances were measured and calculated as functions of cell number, centrifugation force, and ionic strength of the exchanged supernatant. Thus, the structure and properties of cell pellets can be completely characterized by electrical measurements.     

Extracellular Matrix Deposition in Engineered Micromass Cartilage Pellet Cultures: Measurements and Modelling

This article explores possible mechanisms governing extracellular matrix deposition in engineered cartilaginous cell pellets. A theoretical investigation is carried out alongside an experimental study measuring proteoglycan and collagen volume fractions within murine chondrogenic (ATDC-5) cell pellets. The simple mathematical model, which adopts a nutrient-dependent proteoglycan production rate, successfully reproduces the periphery-dominated proteoglycan deposition, characteristic of the growth pattern observed experimentally within pellets after 21 days of culture. The results suggest that this inhomogeneous proteoglycan production is due to nutrient deficiencies at the pellet centre. Our model analysis further indicates that a spatially uniform distribution of proteoglycan matrix could be maintained by initiating the culture process with a smaller-sized pellet. Finally, possible extensions are put forward with an aim to improve the model predictions for the early behaviour, where different mechanisms appear to dominate the matrix production within the pellets.     

Modeling and measurements of fungal growth and morphology in submerged fermentations

Generalizing results from fungal fermentations is difficult due to their high sensitivity toward slight variation in starting conditions, poor reproducibility, and difference in strains. In this study a mathematical model is presented in which oxygen transfer, agitation intensity, dissolved oxygen tension, pellet size, formation of mycelia, the fraction of mycelia in the total biomass, carbohydrate source consumption, and biomass growth are taken into account. Two parameters were estimated from simulation, whereas all others are based on measurements or were taken from literature. Experimental data are obtained from the fermentations in both 2 L and 100 L fermentors at various conditions. Comparison of the simulation with experiments shows that the model can fairly well describe the time course of fungal growth (such as biomass and carbohydrate source concentrations) and fungal morphology (such as pellet size and the fraction of pellets in the total biomass). The model predicts that a stronger agitation intensity leads to a smaller pellet size and a lower fraction of pellets in the total biomass. At the same agitation intensity, pellet size is hardly affected by the dissolved oxygen tension, whereas the fraction of mycelia decreases slightly with an increase of the dissolved oxygen tension in the bulk. All of these are in line with observations at the corresponding conditions.     

Quantitative inference of cellular parameters from microfluidic cell culture systems

Microfluidic cell culture systems offer a convenient way to measure cell biophysical parameters in conditions close to the physiological environment. We demonstrate the application of a mathematical model describing the spatial distribution of nutrient and growth factor concentrations in inferring cellular oxygen uptake rates from experimental measurements. We use experimental measurements of oxygen concentrations in a poly(dimethylsiloxane) (PDMS) microreactor culturing human hepatocellular liver carcinoma cells (HepG2) to infer quantitative information on cellular oxygen uptake rates. We use a novel microchannel design to avoid the parameter correlation problem associated with simultaneous cellular uptake and diffusion of oxygen through the PDMS surface. We find that the cellular uptake of oxygen is dependent on the cell density and can be modeled using a logistic term in the Michaelis–Menten equation. Our results are significant not only for the development of novel assays to quantitatively infer cell response to stimuli, but also for the development, design, and optimization of novel in vitro systems for drug discovery and tissue engineering. Biotechnol. Bioeng. 2009;103: 966–974. © 2009 Wiley Periodicals, Inc.     

Relation between the reaction cytochrome oxidase-oxygen and oxygen uptake in cells in vivo. The role of diffusion

Before reaching the oxidase located inside the cell on the mitochondrial membrane, oxygen may be slowed down by diffusion within the cytoplasm. Diffusion or enzymic activity may predominate and this is related to the size and morphology of the organism, the intracellular diffusion coefficient, and the K_m of the oxidizing terminal enzyme. This K_m may be apparently diminished by inhibitors. Different experimental situations can arise in some plants (where oxygen diffusion in bulk tissue is not important in comparison to that in cell): either oxygen diffusion is limiting, or diffusion and an enzymic reaction compete, or diffusion does not slow down the respiratory rate. A mathematical model relates the oxygen concentration in the external medium, to the rate of oxygen uptake at the oxygen cytochrome-oxidase reaction level.