An in silico model predicts the impact of scaffold design in large bone defect regeneration

Graded or random - Effect of pore distribution in 3D titanium scaffolds on corrosion performance and response of hMSCs

Researchers agree that the ideal scaffold for tissue engineering should possess a 3D and highly porous structure, biocompatibility to encourage cell/tissue growth, suitable surface chemistry for cell attachment and differentiation, and mechanical properties that match those of the surrounding tissues. However, there is no consensus on the optimal pore distribution. In this study, we investigated the effect of pore distribution on corrosion resistance and performance of human mesenchymal stem cells (hMSC) using titanium scaffolds fabricated by laser beam powder bed fusion (PBF-LB). We designed two scaffold architectures with the same porosities (i.e., 75 %) but different distribution of pores of three sizes (200, 500, and 700 μm). The pores were either grouped in three zones (graded, GRAD) or distributed randomly (random, RAND). Microfocus X-ray computed tomography revealed that the chemically polished scaffolds had the porosity of 69 ± 4 % (GRAD) and 71 ± 4 % (RAND), and that the GRAD architecture had the higher surface area (1580 ± 101 vs 991 ± 62 mm2) and the thinner struts (221 ± 37 vs 286 ± 14 μm). The electrochemical measurements demonstrated that the apparent corrosion rate of chemically polished GRAD scaffold decreased with the immersion time extension, while that for polished RAND was increased. The RAND architecture outperformed the GRAD one with respect to hMSC proliferation (over two times higher although the GRAD scaffolds had 85 % higher initial cell retention) and migration from a monolayer. Our findings demonstrate that the pore distribution affects the biological properties of the titanium scaffolds for bone tissue engineering.

Influence of bite force and implant elastic modulus on mandibular reconstruction with particulate-cancellous bone marrow grafts healing: An in silico investigation

This study aims to investigate tissue differentiation during mandibular reconstruction with particulate cancellous bone marrow (PCBM) graft healing using biphasic mechanoregulation theory under four bite force magnitudes and four implant elastic moduli to examine its implications on healing rate, implant stress distribution, new bone elastic modulus, mandible equivalent stiffness, and load-sharing progression. The finite element model of a half Canis lupus mandible, symmetrical about the midsagittal plane, with two marginal defects filled by PCBM graft and stabilized by porous implants, was simulated for 12 weeks. Eight different scenarios, which consist of four bite force magnitudes and four implant elastic moduli, were tested. It was found that the tissue differentiation pattern corroborates the experimental findings, where the new bone propagates from the superior side and the buccal and lingual sides in contact with the native bone, starting from the outer regions and progressing inward. Faster healing and quicker development of bone graft elastic modulus and mandible equivalent stiffness were observed in the variants with lower bite force magnitude and or larger implant elastic modulus. A load-sharing condition was found as the healing progressed, with M3 (Ti6Al4V) being better than M4 (stainless steel), indicating the higher stress shielding potentials of M4 in the long term. This study has implications for a better understanding of mandibular reconstruction mechanobiology and demonstrated a novel in silico framework that can be used for post-operative planning, failure prevention, and implant design in a better way.

Impacts of permeability and effective diffusivity of porous scaffolds on bone ingrowth: In silico and in vivo analyses

The permeability and the effective diffusivity of a porous scaffold are critical in the bone-ingrowth process. However, design guidelines for porous structures are still lacking due to inadequate understanding of the complex physiological processes involved. In this study, a model integrating the fundamental biological processes of bone regeneration was constructed to investigate the roles of permeability and effective diffusivity in regulating bone deposition in scaffolds. The in silico analysis results were confirmed in vivo by examining bone depositions in three diamond lattice scaffolds manufactured using selective laser melting. The findings show that the scaffolds with better permeability and effective diffusivity had deeper bone ingrowth and greater bone volume. Compared to permeability, effective diffusivity exhibited greater sensitivity to the orientation of porous structures, and bone ingrowth was deeper in the directions with higher effective diffusivity in spite of identical pore size. A 4.8-fold increase in permeability and a 1.6-fold increase in effective diffusivity by changing the porous structure led to a 1.5-fold increase in newly formed bone. The effective diffusivity of the porous scaffold affects the distribution of osteogenic growth factor, which in turn impacts cell migration and bone deposition through chemotaxis effects. Therefore, effective diffusivity may be a more suitable indicator for porous scaffolds because our study shows changes in this parameter determine changes in bone distribution and bone volume.

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.

Computational design and evaluation of the mechanical and electrical behavior of a piezoelectric scaffold: a preclinical study

Piezoelectric scaffolds have been recently developed to explore their potential to enhance the bone regeneration process using the concept of piezoelectricity, which also inherently occurs in bone. In addition to providing mechanical support during bone healing, with a suitable design, they are supposed to produce electrical signals that ought to favor the cell responses. In this study, using finite element analysis (FEA), a piezoelectric scaffold was designed with the aim of providing favorable ranges of mechanical and electrical signals when implanted in a large bone defect in a large animal model, so that it could inform future pre-clinical studies. A parametric analysis was then performed to evaluate the effect of the scaffold design parameters with regard to the piezoelectric behavior of the scaffold. The designed scaffold consisted of a porous strut-like structure with piezoelectric patches covering its free surfaces within the scaffold pores. The results showed that titanium or PCL for the scaffold and barium titanate (BT) for the piezoelectric patches are a promising material combination to generate favorable ranges of voltage, as reported in experimental studies. Furthermore, the analysis of variance showed the thickness of the piezoelectric patches to be the most influential geometrical parameter on the generation of electrical signals in the scaffold. This study shows the potential of computer tools for the optimization of scaffold designs and suggests that patches of piezoelectric material, attached to the scaffold surfaces, can deliver favorable ranges of electrical stimuli to the cells that might promote bone regeneration.

Virtual Design of 3D-Printed Bone Tissue Engineered Scaffold Shape Using Mechanobiological Modeling: Relationship of Scaffold Pore Architecture to Bone Tissue Formation

Large bone defects are clinically challenging, with up to 15% of these requiring surgical intervention due to non-union. Bone grafts (autographs or allografts) can be used but they have many limitations, meaning that polymer-based bone tissue engineered scaffolds (tissue engineering) are a more promising solution. Clinical translation of scaffolds is still limited but this could be improved by exploring the whole design space using virtual tools such as mechanobiological modeling. In tissue engineering, a significant research effort has been expended on materials and manufacturing but relatively little has been focused on shape. Most scaffolds use regular pore architecture throughout, leaving custom or irregular pore architecture designs unexplored. The aim of this paper is to introduce a virtual design environment for scaffold development and to illustrate its potential by exploring the relationship of pore architecture to bone tissue formation. A virtual design framework has been created utilizing a mechanical stress finite element (FE) model coupled with a cell behavior agent-based model to investigate the mechanobiological relationships of scaffold shape and bone tissue formation. A case study showed that modifying pore architecture from regular to irregular enabled between 17 and 33% more bone formation within the 4-16-week time periods analyzed. This work shows that shape, specifically pore architecture, is as important as other design parameters such as material and manufacturing for improving the function of bone tissue scaffold implants. It is recommended that future research be conducted to both optimize irregular pore architectures and to explore the potential extension of the concept of shape modification beyond mechanical stress to look at other factors present in the body.

The Concept of Scaffold-Guided Bone Regeneration for the Treatment of Long Bone Defects: Current Clinical Application and Future Perspective

The treatment of bone defects remains a challenging clinical problem with high reintervention rates, morbidity, and resulting significant healthcare costs. Surgical techniques are constantly evolving, but outcomes can be influenced by several parameters, including the patient's age, comorbidities, systemic disorders, the anatomical location of the defect, and the surgeon's preference and experience. The most used therapeutic modalities for the regeneration of long bone defects include distraction osteogenesis (bone transport), free vascularized fibular grafts, the Masquelet technique, allograft, and (arthroplasty with) mega-prostheses. Over the past 25 years, three-dimensional (3D) printing, a breakthrough layer-by-layer manufacturing technology that produces final parts directly from 3D model data, has taken off and transformed the treatment of bone defects by enabling personalized therapies with highly porous 3D-printed implants tailored to the patient. Therefore, to reduce the morbidities and complications associated with current treatment regimens, efforts have been made in translational research toward 3D-printed scaffolds to facilitate bone regeneration. Three-dimensional printed scaffolds should not only provide osteoconductive surfaces for cell attachment and subsequent bone formation but also provide physical support and containment of bone graft material during the regeneration process, enhancing bone ingrowth, while simultaneously, orthopaedic implants supply mechanical strength with rigid, stable external and/or internal fixation. In this perspective review, we focus on elaborating on the history of bone defect treatment methods and assessing current treatment approaches as well as recent developments, including existing evidence on the advantages and disadvantages of 3D-printed scaffolds for bone defect regeneration. Furthermore, it is evident that the regulatory framework and organization and financing of evidence-based clinical trials remains very complex, and new challenges for non-biodegradable and biodegradable 3D-printed scaffolds for bone regeneration are emerging that have not yet been sufficiently addressed, such as guideline development for specific surgical indications, clinically feasible design concepts for needed multicentre international preclinical and clinical trials, the current medico-legal status, and reimbursement. These challenges underscore the need for intensive exchange and open and honest debate among leaders in the field. This goal can be addressed in a well-planned and focused stakeholder workshop on the topic of patient-specific 3D-printed scaffolds for long bone defect regeneration, as proposed in this perspective review.

Histological and Histomorphometric Evaluation of Implanted Photodynamic Active Biomaterials for Periodontal Bone Regeneration in an Animal Study

Recently, our group developed two different polymeric biomaterials with photodynamic antimicrobial surface activity for periodontal bone regeneration. The aim of the present study was to analyze the biocompatibility and osseointegration of these materials in vivo. Two biomaterials based on urethane dimethacrylate (BioM1) and tri-armed oligoester-urethane methacrylate (BioM2) that additionally contained ß-tricalcium phosphate and the photosensitizer mTHPC (meso-tetra(hydroxyphenyl)chlorin) were implanted in non-critical size bone defects in the femur (n = 16) and tibia (n = 8) of eight female domestic sheep. Bone specimens were harvested and histomorphometrically analyzed after 12 months. BioM1 degraded to a lower extent which resulted in a mean remnant square size of 17.4 mm², while 12.2 mm² was estimated for BioM2 (p = 0.007). For BioM1, a total percentage of new formed bone by 30.3% was found which was significant higher compared to BioM2 (8.4%, p < 0.001). Furthermore, BioM1 was afflicted by significant lower soft tissue formation (3.3%) as compared to BioM2 (29.5%). Additionally, a bone-to-biomaterial ratio of 81.9% was detected for BioM1, while 8.5% was recorded for BioM2. Implantation of BioM2 caused accumulation of inflammatory cells and led to fibrous encapsulation. BioM1 (photosensitizer-armed urethane dimethacrylate) showed favorable regenerative characteristics and can be recommended for further studies.

Accelerating bone regeneration using poly(lactic-co-glycolic acid)/hydroxyapatite scaffolds containing duck feet-derived collagen

Collagen, with low antigenicity and excellent cell adhesion, is a biomaterial mainly used for regenerating bone, cartilage, and skin, owing to its biocompatibility and biodegradability. Results from a previous study confirmed that a scaffold mixed with duck feet-derived collagen (DC) and Poly(lactic-co-glycolic acid) (PLGA) reduced inflammatory reaction and increased bone regeneration. To develop an optimal bone substitute we included hydroxyapatite (HAp), a key osteoconductive material, in a DC and PLGA mixture. We fabricated 0, 10, 20, 40, 60, and 80 wt% DC/PLGA/HAp scaffolds and studied their potential for bone tissue engineering. Characteristic analysis of the scaffold and seeding of rabbit bone marrow mesenchymal stem cells (rBMSCs) on the scaffold were conducted to investigate cell proliferation, osteogenic differentiation, and bone formation. We confirmed that increasing DC concentration not only improved the compressive strength of the DC/PLGA/HAp scaffold but also cell proliferation and osteogenic differentiation. It was found through comparison with previous studies that including HAp in the scaffold also promotes osteogenic differentiation. Our study thus shows through in vivo results that the 80 wt% DC/PLGA/HAp scaffold promotes bone mineralization and collagen deposition while reducing the inflammatory response. Hence, 80 wt% DC/PLGA/HAp has excellent potential as a biomaterial for bone regeneration applications.

Trabecular bone remodeling in the aging mouse: A micro-multiphysics agent-based in silico model using single-cell mechanomics

Bone remodeling is regulated by the interaction between different cells and tissues across many spatial and temporal scales. Notably, in silico models are regarded as powerful tools to further understand the signaling pathways that regulate this intricate spatial cellular interplay. To this end, we have established a 3D multiscale micro-multiphysics agent-based (micro-MPA) in silico model of trabecular bone remodeling using longitudinal in vivo data from the sixth caudal vertebra (CV6) of PolgA^{(D257A/D257A)} mice, a mouse model of premature aging. Our in silico model includes a variety of cells as single agents and receptor-ligand kinetics, mechanomics, diffusion and decay of cytokines which regulate the cells' behavior. We highlighted its capabilities by simulating trabecular bone remodeling in the CV6 of five mice over 4 weeks and we evaluated the static and dynamic morphometry of the trabecular bone microarchitecture. Based on the progression of the average trabecular bone volume fraction (BV/TV), we identified a configuration of the model parameters to simulate homeostatic trabecular bone remodeling, here named basal. Crucially, we also produced anabolic, anti-anabolic, catabolic and anti-catabolic responses with an increase or decrease by one standard deviation in the levels of osteoprotegerin (OPG), receptor activator of nuclear factor kB ligand (RANKL), and sclerostin (Scl) produced by the osteocytes. Our results showed that changes in the levels of OPG and RANKL were positively and negatively correlated with the BV/TV values after 4 weeks in comparison to basal levels, respectively. Conversely, changes in Scl levels produced small fluctuations in BV/TV in comparison to the basal state. From these results, Scl was deemed to be the main driver of equilibrium while RANKL and OPG were shown to be involved in changes in bone volume fraction with potential relevance for age-related bone features. Ultimately, this micro-MPA model provides valuable insights into how cells respond to their local mechanical environment and can help to identify critical pathways affected by degenerative conditions and ageing.

Wall Shear Stress Analysis and Optimization in Tissue Engineering TPMS Scaffolds

When designing scaffolds for bone tissue engineering (BTE), the wall shear stress (WSS), due to the fluid flow inside the scaffold, is an important factor to consider as it influences the cellular process involved in new tissue formation. The present work analyzed the average WSS in Schwartz diamond (SD) and gyroid (SG) scaffolds with different surface topologies and mesh elements using computational fluid dynamics (CFD) analysis. It was found that scaffold meshes with a smooth surface topology with tetrahedral elements had WSS levels 35% higher than the equivalent scaffold with a non-smooth surface topology with hexahedral elements. The present work also investigated the possibility of implementing the optimization algorithm simulated annealing to aid in the design of BTE scaffolds with a specific average WSS, with the outputs showing that the algorithm was able to reach WSS levels in the vicinity of 5 mPa (physiological range) within the established limit of 100 iterations. This proved the efficacy of combining CFD and optimization methods in the design of BTE scaffolds.

PCL strut-like scaffolds appear superior to gyroid in terms of bone regeneration within a long bone large defect: An in silico study

The treatment of large bone defects represents a major clinical challenge. 3D printed scaffolds appear as a promising strategy to support bone defect regeneration. The 3D design of such scaffolds impacts the healing path and thus defect regeneration potential. Among others, scaffold architecture has been shown to influence the healing outcome. Gyroid architecture, characterized by a zero mean surface curvature, has been discussed as a promising scaffold design for bone regeneration. However, whether gyroid scaffolds are favourable for bone regeneration in large bone defects over traditional strut-like architecture scaffolds remains unknown. Therefore, the aim of this study was to investigate whether gyroid scaffolds present advantages over more traditional strut-like scaffolds in terms of their bone regeneration potential. Validated bone defect regeneration principles were applied in an in silico modeling approach that allows to predict bone formation in defect regeneration. Towards this aim, the mechano-biological bone regeneration principles were adapted to allow simulating bone regeneration within both gyroid and strut-like scaffolds. We found that the large surface curvatures of the gyroid scaffold led to a slower tissue formation dynamic and conclusively reduced bone regeneration. The initial claim, that an overall reduced zero mean surface curvature would enhance bone formation, could not be confirmed. The here presented approach illustrates the potential of in silico tools to evaluate in pre-clinical studies scaffold designs and eventually lead to optimized architectures of 3D printed implants for bone regeneration.

A mechanobiological computer optimization framework to design scaffolds to enhance bone regeneration

The treatment of large bone defects is a clinical challenge. 3D printed scaffolds are a promising treatment option for such critical-size defects. However, the design of scaffolds to treat such defects is challenging due to the large number of variables impacting bone regeneration; material stiffness, architecture or equivalent scaffold stiffness—due it specific architecture—have all been demonstrated to impact cell behavior and regeneration outcome. Computer design optimization is a powerful tool to find optimal design solutions within a large parameter space for given anatomical constraints. Following this approach, scaffold structures have been optimized to avoid mechanical failure while providing beneficial mechanical stimulation for bone formation within the scaffold pores immediately after implantation. However, due to the dynamics of the bone regeneration process, the mechanical conditions do change from immediately after surgery throughout healing, thus influencing the regeneration process. Therefore, we propose a computer framework to optimize scaffold designs that allows to promote the final bone regeneration outcome. The framework combines a previously developed and validated mechanobiological bone regeneration computer model, a surrogate model for bone healing outcome and an optimization algorithm to optimize scaffold design based on the level of regenerated bone volume. The capability of the framework is verified by optimization of a cylindrical scaffold for the treatment of a critical-size tibia defect, using a clinically relevant large animal model. The combined framework allowed to predict the long-term healing outcome. Such novel approach opens up new opportunities for sustainable strategies in scaffold designs of bone regeneration.