A multiscale computational fluid dynamics approach to simulate the micro-fluidic environment within a tissue engineering scaffold with highly irregular pore geometry

Synthetic biodegradable microporous hydrogels for in vitro 3D culture of functional human bone cell networks

Generating 3D bone cell networks in vitro that mimic the dynamic process during early bone formation remains challenging. Here, we report a synthetic biodegradable microporous hydrogel for efficient formation of 3D networks from human primary cells, analysis of cell-secreted extracellular matrix (ECM) and microfluidic integration. Using polymerization-induced phase separation, we demonstrate dynamic in situ formation of microporosity (5-20 µm) within matrix metalloproteinase-degradable polyethylene glycol hydrogels in the presence of living cells. Pore formation is triggered by thiol-Michael-addition crosslinking of a viscous precursor solution supplemented with hyaluronic acid and dextran. The resulting microporous architecture can be fine-tuned by adjusting the concentration and molecular weight of dextran. After encapsulation in microporous hydrogels, human mesenchymal stromal cells and osteoblasts spread rapidly and form 3D networks within 24 hours. We demonstrate that matrix degradability controls cell-matrix remodeling, osteogenic differentiation, and deposition of ECM proteins such as collagen. Finally, we report microfluidic integration and proof-of-concept osteogenic differentiation of 3D cell networks under perfusion on chip. Altogether, this work introduces a synthetic microporous hydrogel to efficiently differentiate 3D human bone cell networks, facilitating future in vitro studies on early bone development.

An explainable machine learning-based probabilistic framework for the design of scaffolds in bone tissue engineering

Recently, 3D-printed biodegradable scaffolds have shown great potential for bone repair in critical-size fractures. The differentiation of the cells on a scaffold is impacted among other factors by the surface deformation of the scaffold due to mechanical loading and the wall shear stresses imposed by the interstitial fluid flow. These factors are in turn significantly affected by the material properties, the geometry of the scaffold, as well as the loading and flow conditions. In this work, a numerical framework is proposed to study the influence of these factors on the expected osteochondral cell differentiation. The considered scaffold is rectangular with a 0/90 lay-down pattern and a four-layered strut made of polylactic acid with a 5% steel particle content. The distribution of the different types of cells on the scaffold surface is estimated through a scalar stimulus, calculated by using a mechanobioregulatory model. To reduce the simulation time for the computation of the stimulus, a probabilistic machine learning (ML)-based reduced-order model (ROM) is proposed. Then, a sensitivity analysis is performed using the Shapley additive explanations to examine the contribution of the various parameters to the framework stimulus predictions. In a final step, a multiobjective optimization procedure is implemented using genetic algorithms and the ROM, aiming to identify the material parameters and loading conditions that maximize the percentage of surface area populated by bone cells while minimizing the area corresponding to the other types of cells and the resorption condition. The results of the performed analysis highlight the potential of using ROMs for the scaffold design, by dramatically reducing the simulation time while enabling the efficient implementation of sensitivity analysis and optimization procedures.

In silico model development and optimization of in vitro lung cell population growth

Tissue engineering predominantly relies on trial and error in vitro and ex vivo experiments to develop protocols and bioreactors to generate functional tissues. As an alternative, in silico methods have the potential to significantly reduce the timelines and costs of experimental programs for tissue engineering. In this paper, we propose a methodology to formulate, select, calibrate, and test mathematical models to predict cell population growth as a function of the biochemical environment and to design optimal experimental protocols for model inference of in silico model parameters. We systematically combine methods from the experimental design, mathematical statistics, and optimization literature to develop unique and explainable mathematical models for cell population dynamics. The proposed methodology is applied to the development of this first published model for a population of the airway-relevant bronchio-alveolar epithelial (BEAS-2B) cell line as a function of the concentration of metabolic-related biochemical substrates. The resulting model is a system of ordinary differential equations that predict the temporal dynamics of BEAS-2B cell populations as a function of the initial seeded cell population and the glucose, oxygen, and lactate concentrations in the growth media, using seven parameters rigorously inferred from optimally designed in vitro experiments.

Current Advances in the Use of Tissue Engineering for Cancer Metastasis Therapeutics

Cancer is a leading cause of death worldwide and results in nearly 10 million deaths each year. The global economic burden of cancer from 2020 to 2050 is estimated to be USD 25.2 trillion. The spread of cancer to distant organs through metastasis is the leading cause of death due to cancer. However, as of today, there is no cure for metastasis. Tissue engineering is a promising field for regenerative medicine that is likely to be able to provide rehabilitation procedures to patients who have undergone surgeries, such as mastectomy and other reconstructive procedures. Another important use of tissue engineering has emerged recently that involves the development of realistic and robust in vitro models of cancer metastasis, to aid in drug discovery and new metastasis therapeutics, as well as evaluate cancer biology at metastasis. This review covers the current studies in developing tissue-engineered metastasis structures. This article reports recent developments in in vitro models for breast, prostate, colon, and pancreatic cancer. The review also identifies challenges and opportunities in the use of tissue engineering toward new, clinically relevant therapies that aim to reduce the cancer burden.

Mechanical and Computational Fluid Dynamic Models for Magnesium-Based Implants

Today, mechanical properties and fluid flow dynamic analysis are considered to be two of the most important steps in implant design for bone tissue engineering. The mechanical behavior is characterized by Young's modulus, which must have a value close to that of the human bone, while from the fluid dynamics point of view, the implant permeability and wall shear stress are two parameters directly linked to cell growth, adhesion, and proliferation. In this study, we proposed two simple geometries with a three-dimensional pore network dedicated to a manufacturing route based on a titanium wire waving procedure used as an intermediary step for Mg-based implant fabrication. Implant deformation under different static loads, von Mises stresses, and safety factors were investigated using finite element analysis. The implant permeability was computed based on Darcy's law following computational fluid dynamic simulations and, based on the pressure drop, was numerically estimated. It was concluded that both models exhibited a permeability close to the human trabecular bone and reduced wall shear stresses within the biological range. As a general finding, the proposed geometries could be useful in orthopedics for bone defect treatment based on numerical analyses because they mimic the trabecular bone properties.

Lifelike Transformative Materials for Biohybrid Implants: Inspired by Nature, Driven by Technology

Today's living world is enriched with a myriad of natural biological designs, shaped by billions of years of evolution. Unraveling the construction rules of living organisms offers the potential to create new materials and systems for biomedicine. From the close examination of living organisms, several concepts emerge: hierarchy, pattern repetition, adaptation, and irreducible complexity. All these aspects must be tackled to develop transformative materials with lifelike behavior. This perspective article highlights recent progress in the development of transformative biohybrid systems for applications in the fields of tissue regeneration and biomedicine. Advances in computational simulations and data-driven predictions are also discussed. These tools enable the virtual high-throughput screening of implant design and performance before committing to fabrication, thus reducing the development time and cost of biomimetic and biohybrid constructs. The ongoing progress of imaging methods also constitutes an essential part of this matter in order to validate the computation models and enable longitudinal monitoring. Finally, the current challenges of lifelike biohybrid materials, including reproducibility, ethical considerations, and translation, are discussed. Advances in the development of lifelike materials will open new biomedical horizons, where perhaps what is currently envisioned as science fiction will become a science-driven reality in the future.

Characterization of three-dimensional bone-like tissue growth and organization under influence of directional fluid flow

The transition in the field of bone tissue engineering from bone regeneration to in vitro models has come with the challenge of recreating a dense and anisotropic bone-like extracellular matrix (ECM). Although the mechanism by which bone ECM gains its structure is not fully understood, mechanical loading and curvature have been identified as potential contributors. Here, guided by computational simulations, we evaluated cell and bone-like tissue growth and organization in a concave channel with and without directional fluid flow stimulation. Human mesenchymal stromal cells were seeded on donut-shaped silk fibroin scaffolds and osteogenically stimulated for 42 days statically or in a flow perfusion bioreactor. After 14, 28, and 42 days, constructs were investigated for cell and tissue growth and organization. As a result, directional fluid flow was able to improve organic tissue growth but not organization. Cells tended to orient in the tangential direction of the channel, possibly attributed to its curvature. Based on our results, we suggest that organic ECM production but not anisotropy can be stimulated through the application of fluid flow. With this study, an initial attempt in three-dimensions was made to improve the resemblance of in vitro produced bone-like ECM to the physiological bone ECM.

Unique osteogenic profile of bone marrow stem cells stimulated in perfusion bioreactor is Rho-ROCK-mediated contractility dependent

The fate determination of bone marrow mesenchymal stem/stromal cells (BMSC) is tightly regulated by mechanical cues, including fluid shear stress. Knowledge of mechanobiology in 2D culture has allowed researchers in bone tissue engineering to develop 3D dynamic culture systems with the potential for clinical translation in which the fate and growth of BMSC are mechanically controlled. However, due to the complexity of 3D dynamic cell culture compared to the 2D counterpart, the mechanisms of cell regulation in the dynamic environment remain relatively undescribed. In the present study, we analyzed the cytoskeletal modulation and osteogenic profiles of BMSC under fluid stimuli in a 3D culture condition using a perfusion bioreactor. BMSC subjected to fluid shear stress (mean 1.56 mPa) showed increased actomyosin contractility, accompanied by the upregulation of mechanoreceptors, focal adhesions, and Rho GTPase-mediated signaling molecules. Osteogenic gene expression profiling revealed that fluid shear stress promoted the expression of osteogenic markers differently from chemically induced osteogenesis. Osteogenic marker mRNA expression, type 1 collagen formation, ALP activity, and mineralization were promoted in the dynamic condition, even in the absence of chemical supplementation. The inhibition of cell contractility under flow by Rhosin chloride, Y27632, MLCK inhibitor peptide-18, or Blebbistatin revealed that actomyosin contractility was required for maintaining the proliferative status and mechanically induced osteogenic differentiation in the dynamic culture. The study highlights the cytoskeletal response and unique osteogenic profile of BMSC in this type of dynamic cell culture, stepping toward the clinical translation of mechanically stimulated BMCS for bone regeneration.

The Influence of Fluid Shear Stress on Bone and Cancer Cells Proliferation and Distribution

We investigated the potential correlation between the fluid shear stress and the proliferation of bone prostate cancer cells on the surface of nanoclay-based scaffolds in a perfusion bioreactor. Human mesenchymal stem cells (hMSCs) were seeded on the scaffolds to initiate bone growth. After 23 days, prostate cancer cells (MDAPCa2b) were cultured on top of the osteogenically differentiated hMSCs. The scaffolds were separated into two groups subjected to two distinct conditions: (i) static (no flow); and (ii) dynamic (with flow) conditions to recapitulate bone metastasis of prostate cancer. Based on measured data, Computational Fluid Dynamics (CFD) models were constructed to determine the velocity and shear stress distributions on the scaffold surface. Our experimental results show distinct differences in the growth pattern of hMSCs and MDAPCa2b cells between the static and dynamic conditions. Our computational results further suggest that the dynamic flow leads to drastic change in cell morphology and tumorous distribution. Our work points to a strong correlation between tumor growth and local interstitial flows in bones.

Geometry-Based Computational Fluid Dynamic Model for Predicting the Biological Behavior of Bone Tissue Engineering Scaffolds

The use of biocompatible and biodegradable porous scaffolds produced via additive manufacturing is one of the most common approaches in tissue engineering. The geometric design of tissue engineering scaffolds (e.g., pore size, pore shape, and pore distribution) has a significant impact on their biological behavior. Fluid flow dynamics are important for understanding blood flow through a porous structure, as they determine the transport of nutrients and oxygen to cells and the flushing of toxic waste. The aim of this study is to investigate the impact of the scaffold architecture, pore size and distribution on its biological performance using Computational Fluid Dynamics (CFD). Different blood flow velocities (BFV) induce wall shear stresses (WSS) on cells. WSS values above 30 mPa are detrimental to their growth. In this study, two scaffold designs were considered: rectangular scaffolds with uniform square pores (300, 350, and 450 µm), and anatomically designed circular scaffolds with a bone-like structure and pore size gradient (476-979 µm). The anatomically designed scaffolds provided the best fluid flow conditions, suggesting a 24.21% improvement in the biological performance compared to the rectangular scaffolds. The numerical observations are aligned with those of previously reported biological studies.

Implementing computational modeling in tissue engineering: where disciplines meet

In recent years, the mathematical and computational sciences have developed novel methodologies and insights that can aid in designing advanced bioreactors, microfluidic set-ups or organ-on-chip devices, in optimizing culture conditions, or predicting long-term behavior of engineered tissues in vivo. In this review, we introduce the concept of computational models and how they can be integrated in an interdisciplinary workflow for Tissue Engineering and Regenerative Medicine (TERM). We specifically aim this review of general concepts and examples at experimental scientists with little or no computational modeling experience. We also describe the contribution of computational models in understanding TERM processes and in advancing the TERM field by providing novel insights.

Optimization and Validation of a Custom-Designed Perfusion Bioreactor for Bone Tissue Engineering: Flow Assessment and Optimal Culture Environmental Conditions

Various perfusion bioreactor systems have been designed to improve cell culture with three-dimensional porous scaffolds, and there is some evidence that fluid force improves the osteogenic commitment of the progenitors. However, because of the unique design concept and operational configuration of each study, the experimental setups of perfusion bioreactor systems are not always compatible with other systems. To reconcile results from different systems, the thorough optimization and validation of experimental configuration are required in each system. In this study, optimal experimental conditions for a perfusion bioreactor were explored in three steps. First, an in silico modeling was performed using a scaffold geometry obtained by microCT and an expedient geometry parameterized with porosity and permeability to assess the accuracy of calculated fluid shear stress and computational time. Then, environmental factors for cell culture were optimized, including the volume of the medium, bubble suppression, and medium evaporation. Further, by combining the findings, it was possible to determine the optimal flow rate at which cell growth was supported while osteogenic differentiation was triggered. Here, we demonstrated that fluid shear stress up to 15 mPa was sufficient to induce osteogenesis, but cell growth was severely impacted by the volume of perfused medium, the presence of air bubbles, and medium evaporation, all of which are common concerns in perfusion bioreactor systems. This study emphasizes the necessity of optimization of experimental variables, which may often be underreported or overlooked, and indicates steps which can be taken to address issues common to perfusion bioreactors for bone tissue engineering.

Scaffolding Biomaterials for 3D Cultivated Meat: Prospects and Challenges

Cultivating meat from stem cells rather than by raising animals is a promising solution to concerns about the negative externalities of meat production. For cultivated meat to fully mimic conventional meat's organoleptic and nutritional properties, innovations in scaffolding technology are required. Many scaffolding technologies are already developed for use in biomedical tissue engineering. However, cultivated meat production comes with a unique set of constraints related to the scale and cost of production as well as the necessary attributes of the final product, such as texture and food safety. This review discusses the properties of vertebrate skeletal muscle that will need to be replicated in a successful product and the current state of scaffolding innovation within the cultivated meat industry, highlighting promising scaffold materials and techniques that can be applied to cultivated meat development. Recommendations are provided for future research into scaffolds capable of supporting the growth of high-quality meat while minimizing production costs. Although the development of appropriate scaffolds for cultivated meat is challenging, it is also tractable and provides novel opportunities to customize meat properties.

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.

Effects of the 3D Geometry Reconstruction on the Estimation of 3D Porous Scaffold Permeability. ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. ANNUAL INTERNATIONAL CONFERENCE 2021;2021:4403-4407. [PMID: 34892196 DOI: 10.1109/embc46164.2021.9629664] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

3D scaffolds for tissue engineering typically need to adopt a dynamic culture to foster cell distribution and survival throughout the scaffold. It is, therefore, crucial to know fluids' behavior inside the scaffold architecture, especially for complex porous ones. Here we report a comparison between simulated and measured permeability of a porous 3D scaffold, focusing on different modeling parameters. The scaffold features were extracted by microcomputed tomography (µCT) and representative volume elements were used for the computational fluid-dynamic analyses. The objective was to investigate the sensitivity of the model to the degree of detail of the µCT image and the elements of the mesh. These findings highlight the pros and cons of the modeling strategy adopted and the importance of such parameters in analyzing fluid behavior in 3D scaffolds.

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.

Evaluating Directional Dependency of Selective Laser Sintered Patient Specific Biodegradable Devices to Improve Predictive Modeling and Design Verification

Additive manufacturing, or 3D printing, of the bioresorbable polymer [Formula: see text]-polycaprolactone (PCL) is an emerging tissue engineering solution addressing patient specific anatomies. Predictively modeling the mechanical behavior of 3D printed parts comprised of PCL improves the ability to develop patient specific devices that meet design requirements while reducing the testing of extraneous design variants and development time for emergency devices. Predicting mechanical behavior of 3D-printed devices is limited by the variability of effective material moduli that are determined in part by the 3D printing manufacturing process. Powder fusion methods, specifically laser sintering, are known to produce parts with internal porosity ultimately impacting the mechanical performance of printed devices. This study investigates the role of print direction and part size on the material and structural properties of laser sintered PCL parts. Solid PCL cylinders were printed in the XY (perpendicular to laser) and Z direction (parallel to laser), scanned using microcomputed tomography, and mechanically tested under compression. Compositional, structural, and functional properties of the printed parts were evaluated with differential scanning calorimetry, gel permeation chromatography, microcomputed tomography, and mechanical testing. Computational models of printed and scanned cylinders were fit to experimental data to derive effective moduli. Effective moduli were used to predict the mechanical behavior of splints used for emergency repair of severe tracheobronchomalacia. Laser sintering did not cause significant differences in polymer material properties compared to unmanufactured powder. Effective moduli (Eeff) were greater for larger part sizes (p < 0.01) and for parts oriented in the XY direction compared to the Z direction (p < 0.001). These dependencies were congruent with the differences in void volumes associated with the print direction (p < 0.01) and part size (p < 0.01). Finite element models of splint parallel compression tests utilizing the Eeff dependent on print direction and size agreed with experimental closed compression tests of splints. Evaluating the microstructural properties of printed parts and selecting effective moduli for finite element models based on manufacturing parameters allows accurate prediction of device performance. These findings allow testing of a greater number of device design variants in silico to accomodate patient specific anatomies towards providing higher quality parts while lowering overall time and costs of manufacturing and testing.

A machine learning-based multiscale model to predict bone formation in scaffolds

Computational modeling methods combined with non-invasive imaging technologies have exhibited great potential and unique opportunities to model new bone formation in scaffold tissue engineering, offering an effective alternate and viable complement to laborious and time-consuming in vivo studies. However, existing numerical approaches are still highly demanding computationally in such multiscale problems. To tackle this challenge, we propose a machine learning (ML)-based approach to predict bone ingrowth outcomes in bulk tissue scaffolds. The proposed in silico procedure is developed by correlating with a dedicated longitudinal (12-month) animal study on scaffold treatment of a major segmental defect in sheep tibia. Comparison of the ML-based time-dependent prediction of bone ingrowth with the conventional multilevel finite element (FE2) model demonstrates satisfactory accuracy and efficiency. The ML-based modeling approach provides an effective means for predicting in vivo bone tissue regeneration in a subject-specific scaffolding system.

[Application advances in the computational fluid dynamics in tissue engineering]

OBJECTIVE

To review the advances in the computational fluid dynamics (CFD) in tissue engineering.

METHODS

The latest research of CFD applied to tissue engineering were extensively retrieved and analyzed, the optimization of bioreactor design and the simulation of fluid dynamics and cell growth kinetics during tissue regeneration in vitro were mainly reviewed.

RESULTS

The simulation and predictive capabilities of CFD can provide important guidance for the optimization of bioreactor design, and the cultivation of engineering tissue. The accuracy of model prediction results can be further improved by combining with experimental research.

CONCLUSION

As a new and effective research tool, CFD has its unique advantages in the application of tissue engineering. However, a more comprehensive and accurate simulation of the whole process of tissue regeneration still needs further studies.

Digital Twins for Tissue Culture Techniques—Concepts, Expectations, and State of the Art

Techniques to provide in vitro tissue culture have undergone significant changes during the last decades, and current applications involve interactions of cells and organoids, three-dimensional cell co-cultures, and organ/body-on-chip tools. Efficient computer-aided and mathematical model-based methods are required for efficient and knowledge-driven characterization, optimization, and routine manufacturing of tissue culture systems. As an alternative to purely experimental-driven research, the usage of comprehensive mathematical models as a virtual in silico representation of the tissue culture, namely a digital twin, can be advantageous. Digital twins include the mechanistic of the biological system in the form of diverse mathematical models, which describe the interaction between tissue culture techniques and cell growth, metabolism, and the quality of the tissue. In this review, current concepts, expectations, and the state of the art of digital twins for tissue culture concepts will be highlighted. In general, DT’s can be applied along the full process chain and along the product life cycle. Due to the complexity, the focus of this review will be especially on the design, characterization, and operation of the tissue culture techniques.

Multiphysics simulation of a compression-perfusion combined bioreactor to predict the mechanical microenvironment during bone metastatic breast cancer loading experiments

Incurable breast cancer bone metastasis causes widespread bone loss, resulting in fragility, pain, increased fracture risk, and ultimately increased patient mortality. Increased mechanical signals in the skeleton are anabolic and protect against bone loss, and they may also do so during osteolytic bone metastasis. Skeletal mechanical signals include interdependent tissue deformations and interstitial fluid flow, but how metastatic tumor cells respond to each of these individual signals remains underinvestigated, a barrier to translation to the clinic. To delineate their respective roles, we report computed estimates of the internal mechanical field of a bone mimetic scaffold undergoing combinations of high and low compression and perfusion using multiphysics simulations. Simulations were conducted in advance of multimodal loading bioreactor experiments with bone metastatic breast cancer cells to ensure that mechanical stimuli occurring internally were physiological and anabolic. Our results show that mechanical stimuli throughout the scaffold were within the anabolic range of bone cells in all loading configurations, were homogenously distributed throughout, and that combined high magnitude compression and perfusion synergized to produce the largest wall shear stresses within the scaffold. These simulations, when combined with experiments, will shed light on how increased mechanical loading in the skeleton may confer anti-tumorigenic effects during metastasis.

Design and properties of biomimetic irregular scaffolds for bone tissue engineering

The treatment of sizeable segmental bone defects remains a challenge encountered by surgeons. In addition to bone transplantation, porous scaffolds have become a common option. Although the mechanical and biological properties of porous scaffold have recently been the subject of intense research, pore irregularity as a critical characteristic has been poorly explored. Therefore, this study aimed to design an irregular biomimetic scaffold for use in bone tissue engineering applications. The irregular scaffold was based on the Voronoi tessellation method for similarity with the primary histomorphological indexes of bone (porosity, trabecular thickness, cortical bone thickness, and surface to volume ratio). Moreover, a new gradient method was adopted, in which porosity was maintained constant, and the strut diameter was changed to generate a gradient in the irregular scaffold. The permeability and stress concentration characteristics of the irregular scaffold were compared against three conventional scaffolds (the octet, body-centered cubic, pillar body-centered cubic). The results illustrated that the microstructure of the irregular scaffold could be controlled similarly to that of the cortical/cancellous bone unit. Simultaneously, a broad range of permeability was identified for the irregular scaffold, and gradient irregular scaffolds performed better in terms of both permeability and stress distribution than regular scaffolds. This study describes a novel method for the design of irregular scaffolds, which have good controllability and excellent permeability.

A time-dependent mechanobiology-based topology optimization to enhance bone growth in tissue scaffolds

Scaffold-based bone tissue engineering has been extensively developed as a potential means to treatment of large bone defects. To enhance the biomechanical performance of porous tissue scaffolds, computational design techniques have gained growing popularity attributable to their compelling efficiency and strong predictive features compared with time-consuming trial-and-error experiments. Nevertheless, the mechanical stimulus necessary for bone regeneration, which characterizes dynamic nature due to continuous variation in the bone-scaffold construct system as a result of bone-ingrowth and scaffold biodegradation, is often neglected. Thus, this study proposes a time-dependent mechanobiology-based topology optimization framework for design of tissue scaffolds, thereby developing an ongoing favorable microenvironment and ensuring a long-term outcome for bone regeneration. For the first time, a level-set based topology optimization algorithm and a time-dependent shape derivative are developed to optimize the scaffold architecture. In this study, a large bone defect in a simulated 2D femur model and a partial defect in a 3D femur model are considered to demonstrate the effectiveness of the proposed design method. The results are compared with those obtained from stiffness-based topology optimization, time-independent design and typical scaffold constructs, showing significant advantages in continuing bone ingrowth outcomes.

Numerical Analysis of the Influence of Porosity and Pore Geometry on Functionality of Scaffolds Designated for Orthopedic Regenerative Medicine

BACKGROUND

Scaffolds are vital for orthopedic regenerative medicine. Therefore, comprehensive studies evaluating their functionality with consideration of variable parameters are needed. The research aim was to evaluate pore geometry and scaffold porosity influence on first, cell culture efficiency in a perfusion bioreactor and second, osteogenic cell diffusion after its implantation.

METHODS

For the studies, five pore geometries were selected (triangular prism with a rounded and a flat profile, cube, octagonal prism, sphere) and seven porosities (up to 80%), on the basis of which 70 models were created for finite element analyses. First, scaffolds were placed inside a flow channel to estimate growth medium velocity and wall shear stress. Secondly, scaffolds were placed in a bone to evaluate osteogenic cell diffusion.

RESULTS

In terms of fluid minimal velocity (0.005 m/s) and maximal wall shear stress (100 mPa), only cubic and octagonal pores with 30% porosity and spherical pores with 20% porosity fulfilled the requirements. Spherical pores had the highest osteogenic cell diffusion efficiency for porosities up to 30%. For higher porosities, the octagonal prism's pores gave the best results up to 80%, where no differences were noted.

CONCLUSIONS

The data obtained allows for the appropriate selection of pore geometry and scaffold porosity for orthopedic regenerative medicine.

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.

Orbital seeding of mesenchymal stromal cells increases osteogenic differentiation and bone-like tissue formation

In bone tissue engineering (TE), an efficient seeding and homogenous distribution of cells is needed to avoid cell loss and damage as well as to facilitate tissue development. Dynamic seeding methods seem to be superior to the static ones because they tend to result in a more homogeneous cell distribution by using kinetic forces. However, most dynamic seeding techniques are elaborate or require special equipment and its influence on the final bone tissue-engineered construct is not clear. In this study, we applied a simple, dynamic seeding method using an orbital shaker to seed human bone marrow-derived mesenchymal stromal cells (hBMSCs) on silk fibroin scaffolds. Significantly higher cell numbers with a more homogenous cell distribution, increased osteogenic differentiation, and mineral deposition were observed using the dynamic approach both for 4 and 6 hours as compared to the static seeding method. The positive influence of dynamic seeding could be attributed to both cell density and distribution but also nutrient supply during seeding and shear stresses (0.0-3.0 mPa) as determined by computational simulations. The influence of relevant mechanical stimuli during seeding should be investigated in the future, especially regarding the importance of mechanical cues for bone TE applications. Our results highlight the importance of adequate choice of seeding method and its impact on developing tissue-engineered constructs. The application of this simple seeding technique is not only recommended for bone TE but can also be used for seeding similar porous scaffolds with hBMSCs in other TE fields.