A multiphysics 3D model of tissue growth under interstitial perfusion in a tissue-engineering bioreactor

Cell-Based Modeling of Tissue Developing in the Scaffold Pores of Varying Cross-Sections

In this work, we present a mathematical model of cell growth in the pores of a perfusion bioreactor through which a nutrient solution is pumped. We have developed a 2-D vertex model that allows us to reproduce the microscopic dynamics of the microenvironment of cells and describe the occupation of the pore space with cells. In this model, each cell is represented by a polygon; the number of vertices and shapes may change over time. The model includes mitotic cell division and intercalation. We study the impact of two factors on cell growth. On the one hand, we consider a channel of variable cross-section, which models a scaffold with a porosity gradient. On the other hand, a cluster of cells grows under the influence of a nutrient solution flow, which establishes a non-uniform distribution of shear stresses in the pore space. We present the results of numerical simulation of the tissue growth in a wavy channel. The model allows us to obtain complete microscopic information that includes the dynamics of intracellular pressure, the local elastic energy, and the characteristics of cell populations. As we showed, in a functional-graded scaffold, the distribution of the shear stresses in the pore space has a complicated structure, which implies the possibility of controlling the growth zones by varying the pore geometry.

Recent Advances in "Functional Engineering of Articular Cartilage Zones by Polymeric Biomaterials Mediated with Physical, Mechanical, and Biological/Chemical Cues"

Articular cartilage (AC) plays an unquestionable role in joint movements but unfortunately the healing capacity is restricted due to its avascular and acellular nature. While cartilage tissue engineering has been lifesaving, it is very challenging to remodel the complex cartilage composition and architecture with gradient physio-mechanical properties vital to proper tissue functions. To address these issues, a better understanding of the intrinsic AC properties and how cells respond to stimuli from the external microenvironment must be better understood. This is essential in order to take one step closer to producing functional cartilaginous constructs for clinical use. Recently, biopolymers have aroused much attention due to their versatility, processability, and flexibility because the properties can be tailored to match the requirements of AC. This review highlights polymeric scaffolds developed in the past decade for reconstruction of zonal AC layers including the superficial zone, middle zone, and deep zone by means of exogenous stimuli such as physical, mechanical, and biological/chemical signals. The mimicked properties are reviewed in terms of the biochemical composition and organization, cell fate (morphology, orientation, and differentiation), as well as mechanical properties and finally, the challenges and potential ways to tackle them are discussed.

Challenges in computational fluid dynamics applications for bone tissue engineering

Bone injuries or defects that require invasive surgical treatment are a serious clinical issue, particularly when it comes to treatment success and effectiveness. Accordingly, bone tissue engineering (BTE) has been researching the use of computational fluid dynamics (CFD) analysis tools to assist in designing optimal scaffolds that better promote bone growth and repair. This paper aims to offer a comprehensive review of recent studies that use CFD analysis in BTE. The mechanical and fluidic properties of a given scaffold are coupled to each other via the scaffold architecture, meaning an optimization of one may negatively affect the other. For example, designs that improve scaffold permeability normally result in a decreased average wall shear stress. Linked with these findings, it appears there are very few studies in this area that state a specific application for their scaffolds and those that do are focused on in vitro bioreactor environments. Finally, this review also demonstrates a scarcity of studies that combine CFD with optimization methods to improve scaffold design. This highlights an important direction of research for the development of the next generation of BTE scaffolds.

Characterization of the Microflow Through 3D Synthetic Niche Microenvironments Hosted in a Millifluidic Bioreactor

Background: Development of new medicines is a lengthy process with high risk of failure since drug efficacy measured in vitro is difficult to confirm in vivo. Intended to add a new tool aiding drug discovery, the MOAB-NICHOID device was developed: a miniaturized optically accessible bioreactor (MOAB) housing the 3D engineered scaffold NICHOID. The aim of our study was to characterize the microflow through the 3D nichoid microenvironment hosted in the MOAB-NICHOID device.

Methods: We used computational fluid dynamics (CFD) simulations to compute the flow field inside a very fine grid resembling the scaffold microenvironment.

Results: The microflow inside the multi-array of nichoid blocks is fed and locally influenced by the mainstream flow developed in the perfusion chamber of the device. Here we have revealed a low velocity, complex flow field with secondary, backward, or local recirculation micro-flows induced by the intricate architecture of the nichoid scaffold.

Conclusion: Knowledge of the microenvironment inside the 3D nichoids allows planning of cell experiments, to regulate the transport of cells towards the scaffold substrate during seeding or the spatial delivery of nutrients and oxygen which affects cell growth and viability.

Porous Geometry Guided Micro-mechanical Environment Within Scaffolds for Cell Mechanobiology Study in Bone Tissue Engineering

Mechanobiology research is for understanding the role of mechanics in cell physiology and pathology. It will have implications for studying bone physiology and pathology and to guide the strategy for regenerating both the structural and functional features of bone. Mechanobiological studies in vitro apply a dynamic micro-mechanical environment to cells via bioreactors. Porous scaffolds are commonly used for housing the cells in a three-dimensional (3D) culturing environment. Such scaffolds usually have different pore geometries (e.g. with different pore shapes, pore dimensions and porosities). These pore geometries can affect the internal micro-mechanical environment that the cells experience when loaded in the bioreactor. Therefore, to adjust the applied micro-mechanical environment on cells, researchers can tune either the applied load and/or the design of the scaffold pore geometries. This review will provide information on how the micro-mechanical environment (e.g. fluid-induced wall shear stress and mechanical strain) is affected by various scaffold pore geometries within different bioreactors. It shall allow researchers to estimate/quantify the micro-mechanical environment according to the already known pore geometry information, or to find a suitable pore geometry according to the desirable micro-mechanical environment to be applied. Finally, as future work, artificial intelligent - assisted techniques, which can achieve an automatic design of solid porous scaffold geometry for tuning/optimising the micro-mechanical environment are suggested.

Improving cell seeding efficiency through modification of fiber geometry in 3D printed scaffolds

Cell seeding on 3D scaffolds is a very delicate step in tissue engineering applications, influencing the outcome of the subsequent culture phase, and determining the results of the entire experiment. Thus, it is crucial to maximize its efficiency. To this purpose, a detailed study of the influence of the geometry of the scaffold fibers on dynamic seeding efficiency is presented. 3D printing technology was used to realize PLA porous scaffolds, formed by fibers with a non-circular cross-sectional geometry, named multilobed to highlight the presence of niches and ridges. An oscillating perfusion bioreactor was used to perform bidirectional dynamic seeding of MG63 cells. The fiber shape influences the fluid dynamic parameters of the flow, affecting values of fluid velocity and wall shear stress. The path followed by cells through the scaffold fibers is also affected and results in a larger number of adhered cells in multilobed scaffolds compared to scaffolds with standard pseudo cylindrical fibers. Geometrical and fluid dynamic features can also have an influence on the morphology of adhered cells. The obtained results suggest that the reciprocal influence of geometrical and fluid dynamic features and their combined effect on cell trajectories should be considered to improve the dynamic seeding efficiency when designing scaffold architecture.

Engineering functional microvessels in synthetic polyurethane random-pore scaffolds by harnessing perfusion flow

Recently reported biomaterial-based approaches toward prevascularizing tissue constructs rely on biologically or structurally complex scaffolds that are complicated to manufacture and sterilize, and challenging to customize for clinical applications. In the current work, a prevascularization method for soft tissue engineering that uses a non-patterned and non-biological scaffold is proposed. Human fibroblasts and HUVECs were seeded on an ionomeric polyurethane-based hydrogel and cultured for 14 days under medium perfusion. A flow rate of 0.05 mL/min resulted in a greater lumen density in the constructs relative to 0.005 and 0.5 mL/min, indicating the critical importance of flow magnitude in establishing microvessels. Constructs generated at 0.05 mL/min perfusion flow were implanted in a mouse subcutaneous model and intravital imaging was used to characterize host blood perfusion through the construct after 2 weeks. Engineered microvessels were functional (i.e. perfused with host blood and non-leaky) and neovascularization of the construct by host vessels was enhanced relative to non-prevascularized constructs. We report on the first strategy toward engineering functional microvessels in a tissue construct using non-bioactive, non-patterned synthetic polyurethane materials.

Changes in scaffold porosity during bone tissue engineering in perfusion bioreactors considerably affect cellular mechanical stimulation for mineralization

Bone tissue engineering (BTE) experiments in vitro have shown that fluid-induced wall shear stress (WSS) can stimulate cells to produce mineralized extracellular matrix (ECM). The application of WSS on seeded cells can be achieved through bioreactors that perfuse medium through porous scaffolds. In BTE experiments in vitro, commonly a constant flow rate is used. Previous studies have found that tissue growth within the scaffold will result in an increase of the WSS over time. To keep the WSS in a reported optimal range of 10–30 mPa, the applied external flow rate can be decreased over time. To investigate what reduction of the external flow rate during culturing is needed to keep the WSS in the optimal range, we here conducted a computational study, which simulated the formation of ECM, and in which we investigated the effect of constant fluid flow and different fluid flow reduction scenarios on the WSS. It was found that for both constant and reduced fluid flow scenarios, the WSS did not exceed a critical value, which was set to 60 mPa. However, the constant flow velocity resulted in a reduction of the cell/ECM surface being exposed to a WSS in the optimal range from 50% at the start of culture to 18.6% at day 21. Reducing the fluid flow over time could avoid much of this effect, leaving the WSS in the optimal range for 40.9% of the surface at 21 days. Therefore, for achieving more mineralized tissue, the conventional manner of loading the perfusion bioreactors (i.e. constant flow rate/velocity) should be changed to a decreasing flow over time in BTE experiments. This study provides an in silico tool for finding the best fluid flow reduction strategy.

Fluid flow-induced cell stimulation in bone tissue engineering changes due to interstitial tissue formation in vitro

In tissue engineering experiments in vitro, bioreactors have been used for applying wall shear stress (WSS) on cells to regulate cellular activities. To determine the loading conditions within bioreactors and to design tissue engineering products, in silico models are used. Previous in silico studies in bone tissue engineering (BTE) focused on quantifying the WSS on cells and the influence on appositional tissue growth. However, many BTE experiments also show interstitial tissue formation (i.e., tissue infiltrated in the pores rather than growing on the struts - appositional growth), which has not been considered in previous in silico studies. We hereby used a multiscale fluid-solid interaction model to quantify the WSS and mechanical strain on cells with interstitial tissue formation, taken from a reported BTE experiment. The WSS showed a high variation among different interstitial tissue morphologies. This is different to the situation under appositional tissue growth. It is found that a 35% filling of the pores results (by mineralised bone tissue) when the average WSS increases from 1.530 (day 0) to 5.735 mPa (day 28). Furthermore, the mechanical strain on cells caused by the fluid flow was extremely low (at the level of 10^-14 -10^-15 ), comparing to the threshold in a previous mechanobiological theory of osteogenesis (eg, 10^-2 ). The output from this study offers a significant insight of the WSS changes during interstitial tissue growth under a constant perfusion flow rate in a BTE experiment. It has paved the way for optimising the local micro-fluidic environment for interstitial tissue mineralisation.

Channeling Effect and Tissue Morphology in a Perfusion Bioreactor Imaged by X-Ray Microtomography

BACKGROUND

Perfusion bioreactors for tissue engineering hold great promises. Indeed, the perfusion of culture medium enhances species transport and mechanically stimulates the cells, thereby increasing cell proliferation and tissue formation. Nonetheless, their development is still hampered by a lack of understanding of the relationship between mechanical cues and tissue growth.

METHODS

Combining tissue engineering, three-dimensional visualization and numerical simulations, we analyze the morphological evolution of neo-tissue in a model bioreactor with respect to the local flow pattern. NIH-3T3 cells were grown under perfusion for one, two and three weeks on a stack of 2 mm polyacetal beads. The model bioreactor was then imaged by X-ray micro-tomography and local tissue morphology was analyzed. To relate experimental observations and mechanical stimulii, a computational fluid dynamics model of flow around spheres in a canal was developed and solved using the finite element method.

RESULTS

We observe a preferential tissue formation at the bioreactor periphery, and relate it to a channeling effect leading to regions of higher flow intensity. Additionally, we find that circular crater-like tissue patterns form in narrow channel regions at early culture times. Using computational fluid dynamic simulations, we show that the location and morphology of these patterns match those of shear stress maxima. Finally, the morphology of the tissue is qualitatively described as the tissue grows and reorganizes itself.

CONCLUSION

Altogether, our study points out the key role of local flow conditions on the tissue morphology developed on a stack of beads in perfusion bioreactors and provides new insights for effective design of hydrodynamic bioreactors for tissue engineering using bead packings.

Bioengineering tools to speed up the discovery and preclinical testing of vaccines for SARS-CoV-2 and therapeutic agents for COVID-19

This review provides an update for the international research community on the cell modeling tools that could accelerate the understanding of SARS-CoV-2 infection mechanisms and could thus speed up the development of vaccines and therapeutic agents against COVID-19. Many bioengineering groups are actively developing frontier tools that are capable of providing realistic three-dimensional (3D) models for biological research, including cell culture scaffolds, microfluidic chambers for the culture of tissue equivalents and organoids, and implantable windows for intravital imaging. Here, we review the most innovative study models based on these bioengineering tools in the context of virology and vaccinology. To make it easier for scientists working on SARS-CoV-2 to identify and apply specific tools, we discuss how they could accelerate the discovery and preclinical development of antiviral drugs and vaccines, compared to conventional models.

A Hybrid Model of Cartilage Regeneration Capturing the Interactions Between Cellular Dynamics and Porosity

To accelerate the development of strategies for cartilage tissue engineering, models are necessary to investigate the interactions between cellular dynamics and the local microenvironment. We use a discrete framework to capture the individual behavior of cells, modeling experiments where cells are seeded in a porous scaffold or hydrogel and over the time course of a month, the scaffold slowly degrades while cells divide and synthesize extracellular matrix constituents. The movement of cells and the ability to proliferate is a function of the local porosity, defined as the volume fraction of fluid in the surrounding region. A phenomenological approach is used to capture a continuous profile for the degrading scaffold and accumulating matrix, which will then change the local porosity throughout the construct. We parameterize the model by first matching total cell counts in the construct to chondrocytes seeded in a polyglycolic acid scaffold (Freed et al. in Biotechnol Bioeng 43:597-604, 1994). We investigate the influence of initial scaffold porosity on the total cell count and spatial profiles of cell and ECM in the construct. Cell counts were higher at day 30 in scaffolds of lower initial porosity, and similar cell counts were obtained using different models of scaffold degradation and matrix accumulation (either uniform or cell-specific). Using this modeling framework, we study the interplay between a phenomenological representation of scaffold architecture and porosity as well as the potential continuous application of growth factors. We determine parameter regimes where large cellular aggregates occur, which can hinder matrix accumulation and cellular proliferation. The developed modeling framework can easily be extended and can be used to identify optimal scaffolds and culture conditions that lead to a desired distribution of extracellular matrix and cell counts throughout the construct.

Applications of Computer Modeling and Simulation in Cartilage Tissue Engineering

BACKGROUND

Advances in cartilage tissue engineering have demonstrated noteworthy potential for developing cartilage for implantation onto sites impacted by joint degeneration and injury. To supplement resource-intensive in vivo and in vitro studies required for cartilage tissue engineering, computational models and simulations can assist in enhancing experimental design.

METHODS

Research articles pertinent to cartilage tissue engineering and computer modeling were identified, reviewed, and summarized. Various applications of computer modeling for cartilage tissue engineering are highlighted, limitations of in silico modeling are addressed, and suggestions for future work are enumerated.

RESULTS

Computational modeling can help better characterize shear stresses generated by bioreactor fluid flow, refine scaffold geometry, customize the mechanical properties of engineered cartilage tissue, and model rates of cell growth and dynamics. Thus, results from in silico studies can help resourcefully enhance in vitro and in vivo studies; however, the limitations of these studies, such as the underlying assumptions and simplifications applied in each model, should always be addressed and justified where applicable. In silico models should also seek validation and verification when possible.

CONCLUSION

Future studies may adopt similar approaches to supplement in vitro trials and further investigate effects of mechanical stimulation on chondrocyte and stem cell dynamics. Additionally, as precision medicine, machine learning, and powerful open-source software become more popular and accessible, applications of multi-scale and multiphysics computational models in cartilage tissue engineering are expected to increase.

Continuum Modeling and Simulation in Bone Tissue Engineering

Bone tissue engineering is currently a mature methodology from a research perspective. Moreover, modeling and simulation of involved processes and phenomena in BTE have been proved in a number of papers to be an excellent assessment tool in the stages of design and proof of concept through in-vivo or in-vitro experimentation. In this paper, a review of the most relevant contributions in modeling and simulation, in silico, in BTE applications is conducted. The most popular in silico simulations in BTE are classified into: (i) Mechanics modeling and scaffold design, (ii) transport and flow modeling, and (iii) modeling of physical phenomena. The paper is restricted to the review of the numerical implementation and simulation of continuum theories applied to different processes in BTE, such that molecular dynamics or discrete approaches are out of the scope of the paper. Two main conclusions are drawn at the end of the paper: First, the great potential and advantages that in silico simulation offers in BTE, and second, the need for interdisciplinary collaboration to further validate numerical models developed in BTE.

Lattice and continuum modelling of a bioactive porous tissue scaffold

A contemporary procedure to grow artificial tissue is to seed cells onto a porous biomaterial scaffold and culture it within a perfusion bioreactor to facilitate the transport of nutrients to growing cells. Typical models of cell growth for tissue engineering applications make use of spatially homogeneous or spatially continuous equations to model cell growth, flow of culture medium, nutrient transport and their interactions. The network structure of the physical porous scaffold is often incorporated through parameters in these models, either phenomenologically or through techniques like mathematical homogenization. We derive a model on a square grid lattice to demonstrate the importance of explicitly modelling the network structure of the porous scaffold and compare results from this model with those from a modified continuum model from the literature. We capture two-way coupling between cell growth and fluid flow by allowing cells to block pores, and by allowing the shear stress of the fluid to affect cell growth and death. We explore a range of parameters for both models and demonstrate quantitative and qualitative differences between predictions from each of these approaches, including spatial pattern formation and local oscillations in cell density present only in the lattice model. These differences suggest that for some parameter regimes, corresponding to specific cell types and scaffold geometries, the lattice model gives qualitatively different model predictions than typical continuum models. Our results inform model selection for bioactive porous tissue scaffolds, aiding in the development of successful tissue engineering experiments and eventually clinically successful technologies.

Modeling, simulations, and optimization of smooth muscle cell tissue engineering for the production of vascular grafts

The paper presents a transient, continuum, two-phase model of the tissue engineering in fibrous scaffolds, including transport equations for the flowing culture medium, nutrient and cell concentration with transverse and in-plane diffusion and cell migration, a novel feature of local in-plane transport across a phenomenological pore and innovative layer-by-layer cell filling approach. The model is successfully validated for the smooth muscle cell tissue engineering of a vascular graft using crosslinked, electrospun gelatin fiber scaffolds for both static and dynamic cell culture, the latter in a dynamic bioreactor with a rotating shaft on which the tubular scaffold is attached. Parametric studies evaluate the impact of the scaffold microstructure, cell dynamics, oxygen transport, and static or dynamic conditions on the rate and extent of cell proliferation and depth of oxygen accessibility. An optimized scaffold of 75% dry porosity is proposed that can be tissue engineered into a viable and still fully oxygenated graft of the tunica media of the coronary artery within 2 days in the dynamic bioreactor. Such scaffold also matches the mechanical properties of the tunica media of the human coronary artery and the suture retention strength of a saphenous vein, often used as a coronary artery graft.

Flow rates in perfusion bioreactors to maximise mineralisation in bone tissue engineering in vitro

In bone tissue engineering experiments, fluid-induced shear stress is able to stimulate cells to produce mineralised extracellular matrix (ECM). The application of shear stress on seeded cells can for example be achieved through bioreactors that perfuse medium through porous scaffolds. The generated mechanical environment (i.e. wall shear stress: WSS) within the scaffolds is complex due to the complexity of scaffold geometry. This complexity has so far prevented setting an optimal loading (i.e. flow rate) of the bioreactor to achieve an optimal distribution of WSS for stimulating cells to produce mineralised ECM. In this study, we demonstrate an approach combining computational fluid dynamics (CFD) and mechano-regulation theory to optimise flow rates of a perfusion bioreactor and various scaffold geometries (i.e. pore shape, porosity and pore diameter) in order to maximise shear stress induced mineralisation. The optimal flow rates, under which the highest fraction of scaffold surface area is subjected to a wall shear stress that induces mineralisation, are mainly dependent on the scaffold geometries. Nevertheless, the variation range of such optimal flow rates are within 0.5-5 mL/min (or in terms of fluid velocity: 0.166-1.66 mm/s), among different scaffolds. This approach can facilitate the determination of scaffold-dependent flow rates for bone tissue engineering experiments in vitro, avoiding performing a series of trial and error experiments.

Factors Affecting Mass Transport Properties of Poly(ε-caprolactone) Membranes for Tissue Engineering Bioreactors

High porosity and mass transport properties of microfiltration polymeric membranes benefit nutrients supply to cells when used as scaffolds in interstitial perfusion bioreactors for tissue engineering. High nutrients transport is assumed when pore size and porosity of the membrane are in the micrometric range. The present work demonstrates that the study of membrane fouling by proteins present in the culture medium, though not done usually, should be included in the routine testing of new polymer membranes for this intended application. Two poly(ε-caprolactone) microfiltration membranes presenting similar average pore size (approximately 0.7 µm) and porosity (>80%) but different external surface porosity and pore size have been selected as case studies. The present work demonstrates that a membrane with lower surface pore abundance and smaller external pore size (approximately 0.67 µm), combined with adequate hydrodynamics and tangential flow filtration mode is usually more convenient to guarantee high flux of nutrients. On the contrary, having large external pore size (approximately 1.70 µm) and surface porosity would incur important internal protein fouling that could not be prevented with the operation mode and hydrodynamics of the perfusion system. Additionally, the use of glycerol in the drying protocols of the membranes might cause plasticization and a consequent reduction of mass transport properties due to membrane compaction by the pressure exerted to force perfusion. Therefore, preferentially, drying protocols that omit the use of plasticizing agents are recommended.

Relevant biological processes for tissue development with stem cells and their mechanistic modeling: A review

A potential alternative for tissue transplants is tissue engineering, in which the interaction of cells and biomaterials can be optimized. Tissue development in vitro depends on the complex interaction of several biological processes such as extracellular matrix synthesis, vascularization and cell proliferation, adhesion, migration, death, and differentiation. The complexity of an individual phenomenon or of the combination of these processes can be studied with phenomenological modeling techniques. This work reviews the main biological phenomena in tissue development and their mathematical modeling, focusing on mesenchymal stem cell growth in three-dimensional scaffolds.

Efficient Computational Design of a Scaffold for Cartilage Cell Regeneration

Due to the sensitivity of mammalian cell cultures, understanding the influence of operating conditions during a tissue generation procedure is crucial. In this regard, a detailed study of scaffold based cell culture under a perfusion flow is presented with the aid of mathematical modelling and computational fluid dynamics (CFD). With respect to the complexity of the case study, this work focuses solely on the effect of nutrient and metabolite concentrations, and the possible influence of fluid-induced shear stress on a targeted cell (cartilage) culture. The simulation set up gives the possibility of predicting the cell culture behavior under various operating conditions and scaffold designs. Thereby, the exploitation of the predictive simulation into a newly developed stochastic routine provides the opportunity of exploring improved scaffold geometry designs. This approach was applied on a common type of fibrous structure in order to increase the process efficiencies compared with the regular used formats. The suggested topology supplies a larger effective surface for cell attachment compared to the reference design while the level of shear stress is kept at the positive range of effect. Moreover, significant improvement of mass transfer is predicted for the suggested topology.

Spatiotemporal Image Correlation Analysis for 3D Flow Field Mapping in Microfluidic Devices

Microfluidic devices reproducing 3D networks are particularly valuable for nanomedicine applications such as tissue engineering and active cell sorting. There is however a gap in the possibility to measure how the flow evolves in such 3D structures. We show here that it is possible to map 3D flows in complex microchannel networks by combining wide field illumination to image correlation approaches. For this purpose, we have derived the spatiotemporal image correlation analysis of time stacks of single-plane illumination microscopy images. From the detailed analytical and numerical analysis of the resulting model, we developed a fitting method that allows us to measure, besides the in-plane velocity, the out-of-plane velocity component down to v_z ≅ 65 μm/s. We have applied this method successfully to the 3D reconstruction of flows in microchannel networks with planar and 3D ramifications. These different network architectures have been realized by exploiting the great prototyping ability of a 3D printer, whose precision can reach few tens of micrometers, coupled to poly dimethyl-siloxane soft-printing lithography.

A flow bioreactor system compatible with real-time two-photon fluorescence lifetime imaging microscopy

Bioreactors are essential cell and tissue culture tools that allow the introduction of biophysical signals into in vitro cultures. One major limitation is the need to interrupt experiments and sacrifice samples at certain time points for analyses. To address this issue, we designed a bioreactor that combines high-resolution contact-free imaging and continuous flow in a closed system that is compatible with various types of microscopes. The high throughput fluid flow bioreactor was combined with two-photon fluorescence lifetime imaging microscopy (2P-FLIM) and validated. The hydrodynamics of the bioreactor chamber were characterized using COMSOL. The simulation of shear stress indicated that the bioreactor system provides homogeneous and reproducible flow conditions. The designed bioreactor was used to investigate the effects of low shear stress on human umbilical vein endothelial cells (HUVECs). In a scratch assay, we observed decreased migration of HUVECs under shear stress conditions. Furthermore, metabolic activity shifts from glycolysis to oxidative phosphorylation-dependent mechanisms in HUVECs cultured under low shear stress conditions were detected using 2P-FLIM. Future applications for this bioreactor range from observing cell fate development in real-time to monitoring the environmental effects on cells or metabolic changes due to drug applications.

An experimental-numerical investigation on the effects of macroporous scaffold geometry on cell culture parameters

INTRODUCTION

Perfused bioreactors have been demonstrated to be effective in the delivery of nutrients and in the removal of waste products to and from the interior of cell-populated three-dimensional scaffolds. In this paper, a perfused bioreactor hosting a macroporous scaffold provided with a channel is used to investigate transport phenomena and culture parameters on cell growth.

METHODS

MG63 human osteosarcoma cells were seeded on macroporous poly​(ε-caprolactone) scaffolds provided with a channel. The scaffolds were cultured in a perfused bioreactor and in static conditions for 5 days. Cell viability and growth were assessed while the concentration of oxygen, glucose and lactate were measured. An in silico, multiphysics, numerical model was set up to study the fluid dynamics and the mass transport of the nutrients in the perfused bioreactor hosting different scaffold geometries.

RESULTS

The experimental and numerical results indicated that the specific cell metabolic activity in scaffolds cultured under perfusion was 30% greater than scaffolds cultured under static conditions. In addition, the scaffold provided with a channel enabled the shear stress to be controlled, the initial seeding density to be retained, and adequate mass transport and waste removal.

CONCLUSIONS

We show that the macroporous scaffold provided with a channel cultured in a macroscale bioreactor can be a robust reference experimental model system to systematically investigate and assess crucial culture parameters. We also show that such an experimental model system can be employed as a simplified "representative unit" to improve the performance of both perfused culture systems and hollow, fiber-integrated scaffolds for large-scale tissue engineering.

2D µ-Particle Image Velocimetry and Computational Fluid Dynamics Study Within a 3D Porous Scaffold

Transport properties of 3D scaffolds under fluid flow are critical for tissue development. Computational fluid dynamics (CFD) models can resolve 3D flows and nutrient concentrations in bioreactors at the scaffold-pore scale with high resolution. However, CFD models can be formulated based on assumptions and simplifications. μ-Particle image velocimetry (PIV) measurements should be performed to improve the reliability and predictive power of such models. Nevertheless, measuring fluid flow velocities within 3D scaffolds is challenging. The aim of this study was to develop a μPIV approach to allow the extraction of velocity fields from a 3D additive manufacturing scaffold using a conventional 2D μPIV system. The μ-computed tomography scaffold geometry was included in a CFD model where perfusion conditions were simulated. Good agreement was found between velocity profiles from measurements and computational results. Maximum velocities were found at the centre of the pore using both techniques with a difference of 12% which was expected according to the accuracy of the μPIV system. However, significant differences in terms of velocity magnitude were found near scaffold substrate due to scaffold brightness which affected the μPIV measurements. As a result, the limitations of the μPIV system only permits a partial validation of the CFD model. Nevertheless, the combination of both techniques allowed a detailed description of velocity maps within a 3D scaffold which is crucial to determine the optimal cell and nutrient transport properties.

A mathematical model of tissue-engineered cartilage development under cyclic compressive loading

In this work a coupled model of solute transport and uptake, cell proliferation, extracellular matrix synthesis and remodeling of mechanical properties accounting for the impact of mechanical loading is presented as an advancement of a previously validated coupled model for free-swelling tissue-engineered cartilage cultures. Tissue-engineering constructs were modeled as biphasic with a linear elastic solid, and relevant intrinsic mechanical stimuli in the constructs were determined by numerical simulation for use as inputs of the coupled model. The mechanical dependent formulations were derived from a calibration and parametrization dataset and validated by comparison of normalized ratios of cell counts, total glycosaminoglycans and collagen after 24-h continuous cyclic unconfined compression from another dataset. The model successfully fit the calibration dataset and predicted the results from the validation dataset with good agreement, with average relative errors up to 3.1 and 4.3 %, respectively. Temporal and spatial patterns determined for other model outputs were consistent with reported studies. The results suggest that the model describes the interaction between the simultaneous factors involved in in vitro tissue-engineered cartilage culture under dynamic loading. This approach could also be attractive for optimization of culture protocols, namely through the application to longer culture times and other types of mechanical stimuli.

Mathematical modelling of the degradation behaviour of biodegradable metals

A mathematical model for the biodegradation of magnesium is developed in this study to inspect the corrosion behaviour of biodegradable implants. The aim of this study was to provide a suitable framework for the assessment of the corrosion rate of magnesium which includes the process of formation/dissolution of the protective film. The model is intended to aid the design of implants with suitable geometries. The level-set method is used to follow the changing geometry of the implants during the corrosion process. A system of partial differential equations is formulated based on the physical and chemical processes that occur at the implant-medium boundary in order to simulate the effect of the formation of a protective film on the degradation rate. The experimental data from the literature on the corrosion of a high-purity magnesium sample immersed in simulated body fluid is used to calibrate the model. The model is then used to predict the degradation behaviour of a porous orthopaedic implant. The model successfully reproduces the precipitation of the corrosion products on the magnesium surface and the effect on the degradation rate. It can be used to simulate the implant degradation and the formation of the corrosion products on the surface of biodegradable magnesium implants with complex geometries.

A novel dual-frequency loading system for studying mechanobiology of load-bearing tissue

In mechanobiological research, an appropriate loading system is an essential tool to mimic mechanical signals in a native environment. To achieve this goal, we have developed a novel loading system capable of applying dual-frequency loading including both a low-frequency high-amplitude loading and a high-frequency low-amplitude loading, according to the mechanical conditions experienced by bone and articular cartilage tissues. The low-frequency high-amplitude loading embodies the main force from muscular contractions and/or reaction forces while the high-frequency low-amplitude loading represents an assistant force from small muscles, ligaments and/or other tissue in order to maintain body posture during human activities. Therefore, such dual frequency loading system may reflect the natural characteristics of complex mechanical load on bone or articular cartilage than the single frequency loading often applied during current mechanobiological experiments. The dual-frequency loading system is validated by experimental tests using precision miniature plane-mirror interferometers. The dual-frequency loading results in significantly more solute transport in articular cartilage than that of the low-frequency high-amplitude loading regiment alone, as determined by quantitative fluorescence microscopy of tracer distribution in articular cartilage. Thus, the loading system can provide a new method to mimic mechanical environment in bone and cartilage, thereby revealing the in vivo mechanisms of mechanosensation, mechanotransduction and mass-transport, and improving mechanical conditioning of cartilage and/or bone constructs for tissue engineering.

Optimization of a 3D Dynamic Culturing System for In Vitro Modeling of Frontotemporal Neurodegeneration-Relevant Pathologic Features

Frontotemporal lobar degeneration (FTLD) is a severe neurodegenerative disorder that is diagnosed with increasing frequency in clinical setting. Currently, no therapy is available and in addition the molecular basis of the disease are far from being elucidated. Consequently, it is of pivotal importance to develop reliable and cost-effective in vitro models for basic research purposes and drug screening. To this respect, recent results in the field of Alzheimer’s disease have suggested that a tridimensional (3D) environment is an added value to better model key pathologic features of the disease. Here, we have tried to add complexity to the 3D cell culturing concept by using a microfluidic bioreactor, where cells are cultured under a continuous flow of medium, thus mimicking the interstitial fluid movement that actually perfuses the body tissues, including the brain. We have implemented this model using a neuronal-like cell line (SH-SY5Y), a widely exploited cell model for neurodegenerative disorders that shows some basic features relevant for FTLD modeling, such as the release of the FTLD-related protein progranulin (PRGN) in specific vesicles (exosomes). We have efficiently seeded the cells on 3D scaffolds, optimized a disease-relevant oxidative stress experiment (by targeting mitochondrial function that is one of the possible FTLD-involved pathological mechanisms) and evaluated cell metabolic activity in dynamic culture in comparison to static conditions, finding that SH-SY5Y cells cultured in 3D scaffold are susceptible to the oxidative damage triggered by a mitochondrial-targeting toxin (6-OHDA) and that the same cells cultured in dynamic conditions kept their basic capacity to secrete PRGN in exosomes once recovered from the bioreactor and plated in standard 2D conditions. We think that a further improvement of our microfluidic system may help in providing a full device where assessing basic FTLD-related features (including PRGN dynamic secretion) that may be useful for monitoring disease progression over time or evaluating therapeutic interventions.

In silico regenerative medicine: how computational tools allow regulatory and financial challenges to be addressed in a volatile market

The cell therapy market is a highly volatile one, due to the use of disruptive technologies, the current economic situation and the small size of the market. In such a market, companies as well as academic research institutes are in need of tools to advance their understanding and, at the same time, reduce their R&D costs, increase product quality and productivity, and reduce the time to market. An additional difficulty is the regulatory path that needs to be followed, which is challenging in the case of cell-based therapeutic products and should rely on the implementation of quality by design (QbD) principles. In silico modelling is a tool that allows the above-mentioned challenges to be addressed in the field of regenerative medicine. This review discusses such in silico models and focuses more specifically on the bioprocess. Three (clusters of) examples related to this subject are discussed. The first example comes from the pharmaceutical engineering field where QbD principles and their implementation through the use of in silico models are both a regulatory and economic necessity. The second example is related to the production of red blood cells. The described in silico model is mainly used to investigate the manufacturing process of the cell-therapeutic product, and pays special attention to the economic viability of the process. Finally, we describe the set-up of a model capturing essential events in the development of a tissue-engineered combination product in the context of bone tissue engineering. For each of the examples, a short introduction to some economic aspects is given, followed by a description of the in silico tool or tools that have been developed to allow the implementation of QbD principles and optimal design.

Oxygen measurement in interstitially perfused cellularized constructs cultured in a miniaturized bioreactor

AIMS

The possibility of developing engineered tissue in vitro and maintaining the cell viability and functionality is primarily related to the possibility of controlling key culture parameters such as oxygen concentration and cell-specific oxygen consumption. We measured these parameters in a three-dimensional (3D) cellularized construct maintained under interstitially perfused culture in a miniaturized bioreactor.

METHODS

MG63 osteosarcoma cells were seeded at high density on a 3D polystyrene scaffold. The 3D scaffolds were sensorized with sensor foils made of a polymer, which fluoresce with intensity proportional to the local oxygen tension. Images of the sensor foil in contact with the cellularized construct were acquired with a video camera every four hours for six culture days and were elaborated with analytical imaging software to obtain oxygen concentration maps.

RESULTS

The data collected indicate a globally decreasing oxygen concentration profile, with a total drop of 28% after six days of culture and an average drop of 10.5% between the inlet and outlet of the perfused construct. Moreover, by importing the measured oxygen concentration data and the cell counts in a model of mass transport, we calculated the cell-specific oxygen consumption over the whole culture period. The consumption increased with oxygen availability and ranged from 0.1 to 0.7 µmol/h/106 cells.

CONCLUSIONS

The sensors used here allowed a non-invasive, contamination-free and non-destructive oxygen measurement over the whole culture period. This study is the basis for optimization of the culture parameters involved in oxygen supply, in order to guarantee maintenance of cell viability in our system.

A mathematical model and computational framework for three-dimensional chondrocyte cell growth in a porous tissue scaffold placed inside a bi-directional flow perfusion bioreactor

The in vitro chondrocyte cell culture for cartilage tissue regeneration in a perfusion bioreactor is a complex process. Mathematical modeling and computational simulation can provide important insights into the culture process, which would be helpful for selecting culture conditions to improve the quality of the developed tissue constructs. However, simulation of the cell culture process is a challenging task due to the complicated interaction between the cells and local fluid flow and nutrient transport inside the complex porous scaffolds. In this study, a mathematical model and computational framework has been developed to simulate the three-dimensional (3D) cell growth in a porous scaffold placed inside a bi-directional flow perfusion bioreactor. The model was developed by taking into account the two-way coupling between the cell growth and local flow field and associated glucose concentration, and then used to perform a resolved-scale simulation based on the lattice Boltzmann method (LBM). The simulation predicts the local shear stress, glucose concentration, and 3D cell growth inside the porous scaffold for a period of 30 days of cell culture. The predicted cell growth rate was in good overall agreement with the experimental results available in the literature. This study demonstrates that the bi-directional flow perfusion culture system can enhance the homogeneity of the cell growth inside the scaffold. The model and computational framework developed is capable of providing significant insight into the culture process, thus providing a powerful tool for the design and optimization of the cell culture process.

Cell based advanced therapeutic medicinal products for bone repair: Keep it simple?

The development of cell based advanced therapeutic medicinal products (ATMPs) for bone repair has been expected to revolutionize the health care system for the clinical treatment of bone defects. Despite this great promise, the clinical outcomes of the few cell based ATMPs that have been translated into clinical treatments have been far from impressive. In part, the clinical outcomes have been hampered because of the simplicity of the first wave of products. In response the field has set-out and amassed a plethora of complexities to alleviate the simplicity induced limitations. Many of these potential second wave products have remained "stuck" in the development pipeline. This is due to a number of reasons including the lack of a regulatory framework that has been evolving in the last years and the shortage of enabling technologies for industrial manufacturing to deal with these novel complexities. In this review, we reflect on the current ATMPs and give special attention to novel approaches that are able to provide complexity to ATMPs in a straightforward manner. Moreover, we discuss the potential tools able to produce or predict 'goldilocks' ATMPs, which are neither too simple nor too complex.

Computational modelling of the scaffold-free chondrocyte regeneration: a two-way coupling between the cell growth and local fluid flow and nutrient concentration

The in vitro chondrocyte cell culture process in a perfusion bioreactor provides enhanced nutrient supply as well as the flow-induced shear stress that may have a positive influence on the cell growth. Mathematical and computational modelling of such a culture process, by solving the coupled flow, mass transfer and cell growth equations simultaneously, can provide important insight into the biomechanical environment of a bioreactor and the related cell growth process. To do this, a two-way coupling between the local flow field and cell growth is required. Notably, most of the computational and mathematical models to date have not taken into account the influence of the cell growth on the local flow field and nutrient concentration. The present research aimed at developing a mathematical model and performing a numerical simulation using the lattice Boltzmann method to predict the chondrocyte cell growth without a scaffold on a flat plate placed inside a perfusion bioreactor. The model considers the two-way coupling between the cell growth and local flow field, and the simulation has been performed for 174 culture days. To incorporate the cell growth into the model, a control-volume-based surface growth modelling approach has been adopted. The simulation results show the variation of local fluid velocity, shear stress and concentration distribution during the culture period due to the growth of the cell phase and also illustrate that the shear stress can increase the cell volume fraction to a certain extent.

Influence of the scaffold geometry on the spatial and temporal evolution of the mechanical properties of tissue-engineered cartilage: insights from a mathematical model

The production of tissue-engineered cartilage in vitro with inhomogeneous mechanical properties is a problem yet to be solved. Different geometries have been studied to overcome this caveat; however, the reported measurements are limited to average values of some properties and qualitative measures of spatial distributions. We will apply a coupled model to extend knowledge about the introduction of a macrochannel in a scaffold by calculating spatiotemporal patterns for several interest variables related to the remodeling of the mechanical properties. Model parameters were estimated based on experimental data on the temporal patterns of glycosaminoglycans, collagen and compressive Young's modulus for channel-free constructs. The model reproduced the experimental data trends in both geometries, with experimental-numerical correlations between 0.84 and 0.97. The channel had a higher impact on the reduction in spatial heterogeneities and delay of saturation of core properties than in the improvement of average properties. Despite the possible improvement of cell densities for longer periods than 56 days, it is estimated that it will not cause further significant improvements of the mechanical properties. The degrees of spatial heterogeneity of the Young's modulus and permeability in the channeled geometry are 23 and 27 % of the channel-free values. While the average Young's modulus values are in the range of native cartilage, the permeabilities are one to three degrees of magnitude higher than the native cartilage, suggesting that limiting factors such as scaffold porosity and initial permeability are more relevant than scaffold geometry to effectively decrease the tissue permeability.

Modelling and simulation of the chondrocyte cell growth, glucose consumption and lactate production within a porous tissue scaffold inside a perfusion bioreactor

Mathematical and numerical modelling of the tissue culture process in a perfusion bioreactor is able to provide insight into the fluid flow, nutrients and wastes transport, dynamics of the pH value, and the cell growth rate. Knowing the complicated interdependence of these processes is essential for optimizing the culture process for cell growth. This paper presents a resolved scale numerical simulation, which allows one not only to characterize the supply of glucose inside a porous tissue scaffold in a perfusion bioreactor, but also to assess the overall culture condition and predict the cell growth rate. The simulation uses a simplified scaffold that consists of a repeatable unit composed of multiple strands. The simulation results explore some problematic regions inside the simplified scaffold where the concentration of glucose becomes lower than the critical value for the chondrocyte cell viability and the cell growth rate becomes significantly reduced.

Prediction of cell growth rate over scaffold strands inside a perfusion bioreactor

Mathematical and computational modeling of the dynamic process where tissue scaffolds are cultured in perfusion bioreactors is able to provide insight into the cell and tissue growth which can facilitate the design of tissue scaffolds and selection of optimal operating conditions. To date, a resolved-scale simulation of cell growth in the culture process, by taking account of the influences of the supply of nutrients and fluid shear stress on the cells, is not yet available in the literature. This paper presents such a simulation study specifically on cartilage tissue regeneration by numerically solving the momentum, scalar transport and cell growth equations, simultaneously, based on the lattice Boltzmann method. The simulation uses a simplified scaffold that consists of two circular strands placed in tandem inside a microchannel, with the object of identifying the effect of one strand on the other. The results indicate that the presence of the front strand can reduce the cell growth rate on the surface of the rear strand, depending on the distance between them. As such, the present study allows for investigation into the influence of the scaffold geometry on the cell growth rate within scaffolds, thus providing a means to improve the scaffold design and the culture process.

Regenerative orthopaedics: in vitro, in vivo...in silico

In silico, defined in analogy to in vitro and in vivo as those studies that are performed on a computer, is an essential step in problem-solving and product development in classical engineering fields. The use of in silico models is now slowly easing its way into medicine. In silico models are already used in orthopaedics for the planning of complicated surgeries, personalised implant design and the analysis of gait measurements. However, these in silico models often lack the simulation of the response of the biological system over time. In silico models focusing on the response of the biological systems are in full development. This review starts with an introduction into in silico models of orthopaedic processes. Special attention is paid to the classification of models according to their spatiotemporal scale (gene/protein to population) and the information they were built on (data vs hypotheses). Subsequently, the review focuses on the in silico models used in regenerative orthopaedics research. Contributions of in silico models to an enhanced understanding and optimisation of four key elements-cells, carriers, culture and clinics-are illustrated. Finally, a number of challenges are identified, related to the computational aspects but also to the integration of in silico tools into clinical practice.