Engineering 3D parallelized microfluidic droplet generators with equal flow profiles by computational fluid dynamics and stereolithographic printing.

Microfluidic droplet generators excel in generating monodisperse micrometer-sized droplets and particles. However, the low throughput of conventional droplet generators hinders their clinical and industrial translation. Current approaches to parallelize microdevices are challenged by the two-dimensional nature of the standard fabrication methods. Here, we report the facile production of three-dimensionally (3D) parallelized microfluidic droplet generators consisting of stacked and radially multiplexed channel designs. Computational fluid dynamics simulations form the design basis for a microflow distributor that ensures similar flow rates through all droplet generators. Stereolithography is the selected technique to fabricate microdevices, which enables the manufacturing of hollow channels with dimensions as small as 50 µm. The microdevices could be operated up to 4 bars without structural damage, including deformation of channels, or leakage of the on-chip printed Luer-Lok type connectors. The printed microdevices readily enable the production of water-in-oil emulsions, as well as polymer containing droplets that act as templates for both solid and core-shell hydrogel microparticles. The cytocompatibility of the 3D printed device is demonstrated by encapsulating mesenchymal stem cells in hydrogel microcapsules, which results in the controllable formation of stem cell spheroids that remain viable and metabolically active for at least 21 days. Thus, the unique features of stereolithography fabricated microfluidic devices allow for the parallelization of droplet generators in a simple yet effective manner by enabling the realization of (complex) 3D designs.

A PBPK model to study the transfer of α-hexabromocyclododecane (α-HBCDD) to tissues of fast- and slow-growing broilers.

A physiologically based pharmacokinetic (PBPK) model was developed to investigate the production-specific factors involved in the transfer of α-hexabromocyclododecane (α-HBCDD) to broiler meat. The model describes growth and lipid deposition in tissues of fast- (FG) and slow- (SG) growing broilers from hatching to slaughter and simulates the exposure through the ingestion of contaminated feed or expanded polystyrene insulation material. Growth parameters were obtained from the literature while parameters relative to uptake, distribution, and elimination of α-HBCDD were adjusted using results of a previous experiment involving broilers exposed through feed throughout the rearing period or allowed to depurate before slaughter. The model was used to compare the two main edible tissues, breast and leg meat, as well as skin, and to investigate the variability within strain. Between strains and within strain, α-HBCDD assimilation efficiency (AE) is higher when the animals are slaughtered young or heavy. However, increasing slaughter age will lower α-HBCDD concentration in tissues, due to dilution. Based on fresh weight, the concentration of α-HBCDD in breast muscles and skin tends to be lower in SG than in FG broilers (-30 to +10%), while it is 10% to 80% higher in leg muscles. Compared to breast muscles, consuming leg muscles would elicit an exposure 9 and 16 times higher in FG and SG broilers, respectively. The consumption of skin together with muscles would multiply the exposure by up to 36 times compared to breast muscle alone. In case of acute exposure, the α-HBCDD concentration in tissues increased sharply, all the more since the animals are lighter in weight, and then decreased rapidly. In FG broilers, dilution through growth contributed for up to 37%, 28% and 97% to the decontamination of breast muscles, leg muscles and skin, respectively, depending on the duration of depuration before slaughter.

Maximizing neotissue growth kinetics in a perfusion bioreactor: An in silico strategy using model reduction and Bayesian optimization.

In regenerative medicine, computer models describing bioreactor processes can assist in designing optimal process conditions leading to robust and economically viable products. In this study, we started from a (3D) mechanistic model describing the growth of neotissue, comprised of cells, and extracellular matrix, in a perfusion bioreactor set-up influenced by the scaffold geometry, flow-induced shear stress, and a number of metabolic factors. Subsequently, we applied model reduction by reformulating the problem from a set of partial differential equations into a set of ordinary differential equations. Comparing the reduced model results to the mechanistic model results and to dedicated experimental results assesses the reduction step quality. The obtained homogenized model is 105 fold faster than the 3D version, allowing the application of rigorous optimization techniques. Bayesian optimization was applied to find the medium refreshment regime in terms of frequency and percentage of medium replaced that would maximize neotissue growth kinetics during 21 days of culture. The simulation results indicated that maximum neotissue growth will occur for a high frequency and medium replacement percentage, a finding that is corroborated by reports in the literature. This study demonstrates an in silico strategy for bioprocess optimization paying particular attention to the reduction of the associated computational cost.

Computational modelling of local calcium ions release from calcium phosphate-based scaffolds.

A variety of natural or synthetic calcium phosphate (CaP)-based scaffolds are currently produced for dental and orthopaedic applications. These scaffolds have been shown to stimulate bone formation due to their biocompatibility, osteoconductivity and osteoinductivity. The release of the [Formula: see text] ions from these scaffolds is of great interest in light of the aforementioned properties. It can depend on a number of biophysicochemical phenomena such as dissolution, diffusion and degradation, which in turn depend on specific scaffold characteristics such as composition and morphology. Achieving an optimal release profile can be challenging when relying on traditional experimental work alone. Mathematical modelling can complement experimentation. In this study, the in vitro dissolution behaviour of four CaP-based scaffold types was investigated experimentally. Subsequently, a mechanistic finite element method model based on biophysicochemical phenomena and specific scaffold characteristics was developed to predict the experimentally observed behaviour. Before the model could be used for local [Formula: see text] ions release predictions, certain parameters such as dissolution constant ([Formula: see text]) and degradation constant ([Formula: see text]) for each type of scaffold were determined by calibrating the model to the in vitro dissolution data. The resulting model showed to yield release characteristics in satisfactory agreement with those observed experimentally. This suggests that the mathematical model can be used to investigate the local [Formula: see text] ions release from CaP-based scaffolds.

Immersed Boundary Models for Quantifying Flow-Induced Mechanical Stimuli on Stem Cells Seeded on 3D Scaffolds in Perfusion Bioreactors.

Perfusion bioreactors regulate flow conditions in order to provide cells with oxygen, nutrients and flow-associated mechanical stimuli. Locally, these flow conditions can vary depending on the scaffold geometry, cellular confluency and amount of extra cellular matrix deposition. In this study, a novel application of the immersed boundary method was introduced in order to represent a detailed deformable cell attached to a 3D scaffold inside a perfusion bioreactor and exposed to microscopic flow. The immersed boundary model permits the prediction of mechanical effects of the local flow conditions on the cell. Incorporating stiffness values measured with atomic force microscopy and micro-flow boundary conditions obtained from computational fluid dynamics simulations on the entire scaffold, we compared cell deformation, cortical tension, normal and shear pressure between different cell shapes and locations. We observed a large effect of the precise cell location on the local shear stress and we predicted flow-induced cortical tensions in the order of 5 pN/µm, at the lower end of the range reported in literature. The proposed method provides an interesting tool to study perfusion bioreactors processes down to the level of the individual cell's micro-environment, which can further aid in the achievement of robust bioprocess control for regenerative medicine applications.

Spatial optimization in perfusion bioreactors improves bone tissue-engineered construct quality attributes.

Perfusion bioreactors have shown great promise for tissue engineering applications providing a homogeneous and consistent distribution of nutrients and flow-induced shear stresses throughout tissue-engineered constructs. However, non-uniform fluid-flow profiles found in the perfusion chamber entrance region have been shown to affect tissue-engineered construct quality characteristics during culture. In this study a whole perfusion and construct, three dimensional (3D) computational fluid dynamics approach was used in order to optimize a critical design parameter such as the location of the regular pore scaffolds within the perfusion bioreactor chamber. Computational studies were coupled to bioreactor experiments for a case-study flow rate. Two cases were compared in the first instance seeded scaffolds were positioned immediately after the perfusion chamber inlet while a second group was positioned at the computationally determined optimum distance were a steady state flow profile had been reached. Experimental data showed that scaffold location affected significantly cell content and neo-tissue distribution, as determined and quantified by contrast enhanced nanoCT, within the constructs both at 14 and 21 days of culture. However, gene expression level of osteopontin and osteocalcin was not affected by the scaffold location. This study demonstrates that the bioreactor chamber environment, incorporating a scaffold and its location within it, affects the flow patterns within the pores throughout the scaffold requiring therefore dedicated optimization that can lead to bone tissue engineered constructs with improved quality attributes.