Trends in Biotechnology
Available online 28 September 2024
Journal home page for Trends in Biotechnology

Review
Engineering next-generation oxygen-generating scaffolds to enhance bone regeneration

工程技术TOPEI检索SCI升级版 工程技术1区SCI基础版 工程技术1区IF 14.3SWJTU A++
https://doi.org/10.1016/j.tibtech.2024.09.006 Get rights and content

Highlights

  • Oxygen (O2) signaling in the bone microenvironment is tightly regulated with spatiotemporal precision, and this poses challenges when replicating these conditions for healing bone defects.
  • In bone tissue engineering (BTE), oxygen-generating scaffolds (OGS) provide spatiotemporal control of the O2 supply that is essential for cell viability and tissue regeneration, particularly in poorly vascularized and critical-sized bone defects.
  • Current efforts in BTE trend towards the development of stimulus-responsive 'smart' biomaterials for OGS.
  • Emerging research emphasizes the development of new imaging techniques to monitor the spatiotemporal distribution of O2 and reactive oxygen species (ROS) in vivo.
  • Combining data-driven machine learning (ML) strategies, in vivo imaging techniques, and materials science with biomanufacturing provides a unique opportunity for the design of novel OGS.

Abstract

In bone, an adequate oxygen (O2) supply is crucial during development, homeostasis, and healing. Oxygen-generating scaffolds (OGS) have demonstrated significant potential to enhance bone regeneration. However, the complexity of O2 delivery and signaling in vivo makes it challenging to tailor the design of OGS to precisely meet this biological requirement. We review recent advances in OGS and analyze persisting engineering and translational hurdles. We also discuss the potential of computational and machine learning (ML) models to facilitate the integration of novel imaging data with biological readouts and advanced biomanufacturing technologies. By elucidating how to tackle current challenges using cutting-edge technologies, we provide insights for transitioning from traditional to next-generation OGS to improve bone regeneration in patients.

Keywords

oxygen generation
biomaterials
hypoxia
imaging
computational modeling
machine learning

Oxygen-generating scaffolds in bone tissue engineering: impact and opportunities

Bone tissue engineering (BTE) combines stem cells, bioactive factors, and/or biodegradable carriers and has been developed as a promising therapeutic alternative to autografts (see Glossary) and allografts for treating critical-sized bone defects. A major challenge in developing scaffolds to regenerate bone in large injuries is an inadequate supply of O2 to cells within the regenerating tissue, and this can lead to cell apoptosis, tissue necrosis, and delayed healing due to hypoxic environments [1]. Oxygenated biomaterials can overcome ischemia and stimulate tissue healing and bone regeneration [2]. Within large BTE constructs, oxygenated scaffolds are broadly categorized into O2-carrying scaffolds (OCS) or OGS, both of which enable spatiotemporal control of O2 delivery. In preclinical or animal models, OCS and OGS have resulted in enhanced bone regeneration [3,4]. However, improving the scaffold design requires that the field addresses deficits in three broad areas: (i) biological underpinnings – what are the underlying mechanisms via which O2 enhances bone healing, and what is the optimal delivery profile (dosing and timing) necessary to maximize bone regeneration? (ii) Technical limitations in the ability to manufacture scaffolds with the appropriate features to provide O2 with tailored spatiotemporal profiles based on dynamic physiological demands. (iii) Overcoming the logistic and translational barriers to using OCS and OGS in clinical settings. This review synthesizes recent advances in the design and preclinical applications of OCS and OGS while highlighting ongoing challenges and providing insights into how an integrative framework that incorporates cutting-edge but underexplored technologies can herald the design and clinical translation of next-generation OGS.

O2 as a cellular signaling molecule

O2 is a major signaling molecule that regulates cellular proliferation and differentiation. However, defining normoxia, hypoxia, and hyperoxia for specific biological tissues or organs is complex, as O2 needs vary among different tissues [5]. For instance, the partial O2 pressure (pO2) of arterial blood is only 95 mmHg whereas atmospheric pO2 is 160 mmHg (20%). This is lower in the periosteum, cortical bone, and bone marrow; the latter has a pO2 of <32 mmHg (4%) [6]. Therefore, the level of atmospheric O2 is hyperoxic for the bone microenvironment, increases the production of reactive oxygen species (ROS), and has detrimental effects on osseous cells [7]. Hence, to develop optimized OCS and OGS for bone regeneration, a thorough understanding of variations in the O2 levels in key cell types within the regenerating bone microenvironment is essential.

Effect of O2 on osteogenesis during bone healing

O2 tensions of 6.6–8.6% are physiologically relevant for overall bone formation and homeostasis [8., 9., 10.]. Even so, excessively low pO2 has detrimental effects on bone cells [11]. For example, exposure to 2% O2 resulted in a tenfold decrease in bone formation, nearly halting it when reduced to 0.2% [12]. This reduction may be due to decreased Runx2 expression [13], inhibition of phosphatidylinositol 3-kinase (PI3K)/Akt signaling pathways [14], or impaired collagen crosslinking which disrupts bone matrix mineralization [12]. Furthermore, osteoclast activity and number increase in response to 2% O2, with a tenfold increase in resorption pit formation and a 21-fold increase in osteoclast activity. During fracture healing, which resembles the healing process of artificial implants, osteoprogenitors were recruited early and proliferated during soft callus formation, and their activity was crucial during hard callus formation and bone remodeling [15]. O2 tension was measured to be 0.8% at the 4 day post-fracture hematoma stage and increased to only 3.8% at the 2 week stage of newly formed fibrous bone [7], suggesting that elevating O2 tension for sustained periods following injury may enhance tissue healing (Figure 1).
Figure 1
  1. Download: Download high-res image (345KB)
  2. Download: Download full-size image

Figure 1. Cellular activity profile based on cell types and O2 demand and desired supply for the design of an oxygen-generating scaffold (OGS).

(A) The cellular activity of osteoprogenitors [122,123], neovascularization processes [33], and M1/M2 macrophages [23,26., 27., 28., 29.] are illustrated for the four stages of bone healing (hematoma formation, soft callus formation, hard callus formation, and bone remodeling). (B) The black curve represents the tissue O2 partial pressure (pO2) at different bone-healing stages without OGS intervention [7,17., 18., 19.]. The healthy bone pO2 range of 6.6–8.6% indicates normal O2 levels in bone homeostasis [8., 9., 10.]. Based on the deficient O2 profile, one can predict a desired O2 supply profile (purple) provided by OGS to meet real-time O2 needs during bone healing. Figure created with BioRender.

Effect of O2 on angiogenesis during bone healing

Hypoxia promotes angiogenesis through multifaceted signaling involving a combination of angiogenic factors and inflammatory cytokines. Although acute hypoxia appears to be essential for stimulating wound healing, prolonged or chronic hypoxia lasting for >7 days has been shown to delay the tissue repair process [16]. Despite vascular infiltration, low local O2 tension persists longer with increasing defect size, extending from 2 weeks [17,18] to 10 weeks post-injury [19]. Since angiogenesis is crucial for bone healing and regeneration, a major concern with O2-releasing scaffolds was that they might hinder vascular infiltration and subsequent tissue regeneration. However, O2-generating calcium peroxide (CaO2)/gelatin microspheres developed for treating osteonecrosis demonstrated enhanced angiogenesis, and CD31+ vessels were prevalent in the defect area 4 weeks post-surgery [20]. Similarly, injectable sodium alginate/carboxymethyl chitosan (CMC)/CaO2 hydrogel led to ectopic osteogenesis in rat subcutaneous defects [21], and significantly greater blood vessel volumes were observed in scaffolds with 25% CaO2 compared with control poly(ε-caprolactone) (PCL) scaffolds [22]. Although not statistically significant, a higher density of 'type H' vessels was observed in CaO2PCL scaffold groups, suggesting that these proregenerative vessel subtypes may be characteristic of the oxygenated microenvironment [22].

Effects of O2 on immune responses during bone healing

Hypoxia elevates hypoxia-inducible factor 1α (HIF-1α) expression and induces the recruitment of neutrophils, macrophages, endothelial cells, and subsequently fibroblasts to the wound site [23]. In hypoxic environments, macrophages switch their metabolism from oxidative phosphorylation towards glycolysis, and this biases them toward the M1 phenotype [24,25]. Although the prolonged presence of M1 macrophages causes inflammation and hinders healing, they are crucial for the initiation of bone healing through clearing cell debris and pathogens, and the recruitment of more immune cells [26]. Thus, ideally, acute hypoxia for up to 48 h might benefit bone healing. Subsequent angiogenesis alleviates hypoxia and facilitates the transition of macrophages from M1 to M2, thus promoting tissue repair [27,28]. M2 macrophages continue to be active during the ossification phase and support bone remodeling (Figure 1A) [29]. Oxygenated biomaterials with a transient O2 release profile provide an opportunity to alter immune cell behavior from proinflammatory to proregenerative within a few days post-implantation. Efforts have been directed towards investigating the immunomodulatory potential of ROS-scavenging bone scaffolds which can influence M1 to M2 macrophage polarization [27,30]. O2 release from OGS could facilitate such a transition and stabilize bone repair through M2 polarization.
In summary, O2 generation and/or release from an implant should maintain O2 tension between 6.6 and 8.6%, a level that is crucial for bone formation and homeostasis, until blood vessel infiltration is stabilized – a process that typically takes at least 2 weeks (Figure 1B) [17]. An acute moderately hypoxic environment (1–5% O2 for 24–48 h) immediately post-implantation could enhance immune response and angiogenesis [26,31,32], although this is highly dependent on the cell type and might require further consideration regarding the specific experimental conditions. As the BTE construct size increases to a clinically relevant scale, O2 deficiency would progressively become a problem, necessitating adjustments in the O2 delivery profile to ensure immediate O2 availability. Consequently, the release rate of oxygenated scaffolds designed for bone repair should match the physiological O2 demand within the defect site.

In vivo O2 measurements during bone regeneration

Measurements of O2 in vivo following fractures could provide further insights into the desired biomaterial designs. However, it is challenging to characterize the complex in vivo dynamics of oxygenation during bone-healing processes. Recent advances in in vivo imaging techniques have provided insights into the role of O2 during bone healing and how tissue-engineering methods can be applied more efficaciously (Table 1). The different imaging methods for characterizing O2 in vivo have their own advantages and disadvantages, as summarized in Table 1. For example, high-resolution optical imaging techniques can suffer from poor tissue penetration [18,33., 34., 35., 36., 37., 38.], whereas methods with deeper tissue penetration tend to have lower spatial resolution [39., 40., 41.] or require complex hardware [42] or specialized contrast agents [43., 44., 45.]. Clinically relevant methods for assessing O2 may provide excellent tissue depth coverage, but may lack sensitivity or have limited probe availability [46., 47., 48.].

Table 1. Methods suitable for imaging OGSa

Imaging modalityMeasured parameterTissue penetrationSpatial resolutionAdvantagesLimitationsRefs
Imaging modalities measuring O2-related parameters
IOS imagingSO2<500 μm5–100 μmEasy to implementLow penetration, 2D[33,34]
PA imagingSO2>5 mm50–150 μmProvides both structural and functional informationComplex system[42]
OCTSO21–3 mm1–15 μmHigh resolutionLimited penetration[35., 36., 37.]
EPRpO25–10 mm for high frequency, 8 cm or more for low frequency0.5–5 mmHigh accuracy and precision for O2 measurementsLow resolution[39., 40., 41.]
Imaging modality measuring ROS
UltrasoundROS>10mm20–200 μmCan mediate ROS generationRelatively low resolution, a contrast agent (microbubbles) is needed[43., 44., 45.]
Imaging modalities measuring O2-related parameters and ROS
Two/multiphotonpO2, ROS500 μm–1 mmSub-micronHigh resolution, depth-resolvedLow penetration, small FOV[18,38]
MRIStO2, ROSFull body5–200 μmExcellent imaging penetration depth, clinically relevantLow sensitivity, limited molecular probes[46., 47., 48.]
BLIHypoxia, ROS10mm2–3 mmHigh SNR and selectivityLow resolution[124,125]
PETHypoxia, ROSFull body1–2 mmHigh sensitivity, various probes availableLow resolution[126,127]
a
Abbreviations: BLI, bioluminescence imaging; EPR, electron paramagnetic resonance; IOS, intrinsic optical signal; MRI, magnetic resonance imaging; OCT, optical coherence tomography; PA imaging, photoacoustic imaging; PET, positron emission tomography; PO2, oxygen partial pressure; ROS, reactive oxygen species; SNR, signal to noise ratio; SO2, oxygen saturation; StO2, tissue oxygen saturation.
Given these developments, we foresee that future research on in vivo O2 imaging will be focused on combining the strengths of extant methods while mitigating their inherent weaknesses. This could include the development of multimodality imaging approaches that integrate complementary image contrast mechanisms to provide more comprehensive O2 mapping during the bone-healing cascade [49]. In addition, there is a dearth of imaging techniques tailored for real-time in vivo imaging of OGS designed for bone-healing applications. Although several in vivo imaging methods have been employed to study OGS within wound-healing contexts [50], and other in vivo O2 imaging methods, such as optical coherence tomography (OCT) [51] and electron paramagnetic resonance (EPR) imaging [52], have been utilized for tissue-engineering applications, their use for bone imaging remains limited due to issues related to their widespread availability, tissue penetration, and light scattering. Furthermore, there is a crucial need for novel noninvasive methods with high spatial and temporal resolution to characterize OGS efficacy in vivo. Enhancing the specificity and sensitivity of O2-sensitive imaging probes and sensors, as well as reducing the cost and complexity of imaging hardware, will be essential to advance the field and make such methods accessible and clinically translatable. Looking ahead, we anticipate that technical advances in this area will include the synthesis of novel O2-sensitive imaging probes [53] or sensors, or wearable devices [54] that enable real-time, localized imaging of O2 release within bone.
The methodological advances described earlier are essential for the development of next-generation OGS because they not only can enhance our understanding of the role of O2 in bone healing but can also inform the design of better biomaterials. For example, exciting results from prior imaging studies have guided the design of novel OGS (Figure 2). Multiphoton microscopy has revealed how the early stages of bone repair often involve severe hypoxia, whereas the later stages exhibit O2 fluctuations, indicating the need for targeted oxygenation strategies at the different stages of healing (Figure 2A,B) [18]. Furthermore, continuous monitoring of O2 with multi-wavelength intrinsic optical signal (IOS) imaging revealed intricate changes in intravascular oxygenation in the spatial and temporal domains [33], highlighting the need to factor in spatiotemporal variations in O2 when optimizing OGS design and bone regeneration (Figure 2C,D). In addition, photoacoustic imaging revealed the role of O2 availability for successful healing of union and non-union bone defects, wherein the oxygen saturation (SO2) in non-union defects was much lower compared with union defects at various phases of healing (Figure 2E,F) [19]. Although still in the nascent stage within the context of OGS development, in vivo imaging technologies hold great promise for unraveling the mechanisms underlying bone oxygenation dynamics. The insights derived from these cutting-edge imaging techniques underscore their potential to inform therapeutic tissue-engineering approaches and improve the management of skeletal injuries.
Figure 2
  1. Download: Download high-res image (829KB)
  2. Download: Download full-size image

Figure 2. In vivo imaging provides novel insights into O2 changes during bone healing.

(A,B) Multiphoton microscopy was used to detect changes in O2 tension during calvarial bone healing [18]. (A) PtP-C343, an O2-sensitive probe, was integrated into an electrospun fiber mesh (green) to enable fluorescence-based measurements of O2 tension during calvarial defect (broken-line circle) healing through a cranial window over 32 days. Multiple regions of interest (labeled A–C) were measured within the defect. (B) Quantitative analysis showed a significant decrease in O2 in all the measured regions on the 7th day post-surgery, followed by their recovery over the next few weeks (*P < 0.05). (C,D) In vivo multiwavelength intrinsic optical signal (IOS) imaging reveals intricate changes in intravascular SO2 during calvarial bone defect healing [33]. (C) Intravascular SO2 computed from multiwavelength IOS images indicates that angiogenic vessels near the defect edge (white broken-line circle) exhibited elevated SO2 during the second week of healing (days D10–D16) compared with the rest of the vasculature. (D) Radial plots for the same defect illustrate the spatiotemporal evolution of mean SO2 over 4 weeks. (E,F) Photoacoustic imaging demonstrated differences in O2 saturation during union and non-union bone healing [19]. (E) Hemoglobin O2 saturation was measured in union and non-union mice long bone defects at weeks 2, 5, and 10 after injury. (F) Quantitative analysis revealed significantly lower SO2 in non-union defects compared with union defects at different phases (*P < 0.05 vs. union; a < 0.05 vs. 2 weeks; b < 0.05 vs. 2 and 5 weeks). Image adapted, with permission, from [18,19,33]. Abbreviations: MPSLM, multiphoton laser scanning microscopy; PtP-C343, platinum porphyrin-coumarin-343; SO2, oxygen saturation.

Recent advances in scaffold design

OCS in bone healing

OCS do not generate O2 on their own but serve as vehicles for transporting and delivering O2 at a desired dose to target sites. Examples include hemoglobin-based oxygen carriers (HBOCs), lipid-based O2 microbubbles, O2-laden nanosponges, polymeric microtanks, and perfluorocarbons (PFC). O2 release by these systems is rate-limited by the diffusion of O2 through the polymeric or lipid membrane. For example, O2 was hyperbarically loaded into hollow biodegradable microtanks and incorporated into PCL for subsequent 3D printing into scaffolds. The scaffolds released O2 for up to 8 h – a relatively short period – but enhanced the deposition of extracellular matrix (ECM) in murine calvarial and subcutaneous defects [55], suggesting that even short O2 release at the earliest stages of healing can be beneficial. Owing to their excellent O2 solubility, PFC-based emulsion-loaded hollow microparticles enabled controlled release of O2 for 10 days [56]. They maintained the survival and osteogenic differentiation potential of human periosteal-derived cells under hypoxic conditions and provided a conducive environment for enhanced bone regeneration in miniature pig mandibular defects. However, for PFCs, long-term stability remains a significant concern, as the emulsions often decay during the freeze–thaw storage process or destabilize due to Ostwald ripening [57]. Oxygent, an FDA-approved PFC product, has a half-life of only 12–48 h, which limits its use primarily to temporary blood substitution during surgery [58]. These issues pose significant challenges for the clinical translation of OCS in BTE applications.

OGS in bone healing

Compared with OCS, the maximum O2 load in OGS may be two to four orders of magnitude higher, allowing O2 release lasting from several days to a month. This longer release period potentially provides an advantage in the context of treating critical-sized bone defects. Examples include solid peroxides (sodium percarbonates, manganese oxide, CaO2, etc.), liquid peroxides (hydrogen peroxide, H2O2), and photosynthetic algae [59], of which the solid peroxides are the most promising for BTE applications Sustained O2 release from solid peroxides is achieved by their chemical reaction with the aqueous environment. O2-generating particles are usually encapsulated within micro/nanocarriers such as liposomes, dendrimers, exosomes, nanospheres, nanocapsules, solid–lipid nanoparticles, nanofibers, or polymeric micelles [60], and O2 release from these systems depends on water penetration and polymeric matrix biodegradation.
Among solid peroxides, CaO2 is the most prevalently researched in BTE and has been employed in different formulations such as coatings [4,61], composite microspheres and filaments [62., 63., 64.], electrospun nanofibers [60,65,66], and hydrogel systems [3,20,21,67,68]. Used in hydrogel systems [3,67], near-complete bone regeneration in critical-sized calvaria defects was observed after 12 weeks [3], but hydrogels are not suitable for load-bearing applications. Alternatively, CaO2 has been coated onto robocasted biphasic calcium phosphate (BCP) scaffolds [61]. When implanted into a 15 mm segmental radial defect in rabbits, a nearly twofold increase in bone formation was observed 6 months after surgery compared with uncoated scaffolds [4]. In polymer–CaO2 composites, the polymer matrix serves as a hydrophobic barrier that slows O2 release by limiting CaO2–water interactions. However, in electrospinning systems [60,65,66], sustained O2 release does not exceed 2 weeks, probably due to their high surface area. Nevertheless, CaO2-based OGS have shown efficacy in promoting bone growth. PCL–CaO2 microparticles developed by Suvarnapathaki and coworkers demonstrated high viability of preosteoblasts and reduced apoptosis, as confirmed by minimal in vitro caspase 3/7 activity. In vivo results from the same study also revealed that their scaffolds accelerated new bone formation, with less fibrous tissue and better tissue integration over 12 weeks, suggesting long-term safety and efficacy [3]. Gelatin–CaO2 microparticles gave similar in vivo results over 10 weeks, and showed low immunogenicity and a significant reduction in inflammatory cells from the first to the third week post-implantation [69]. These studies collectively highlight the promise of CaO2-based scaffolds in promoting bone regeneration for future clinical applications, while maintaining biocompatibility and reducing adverse immune responses.
Magnesium peroxide (MgO2), another solid peroxide, exhibits the slowest and most persistent release kinetics of O2 owing to its low decomposition rate [70]. A 3D-printed scaffold composed of PCL, β-tricalcium phosphate (β-TCP), and MgO2 created by Peng and colleagues [70] released O2 over 21 days and facilitated the formation of new bone in a femoral condyle defect model. Mg2+ ions also directly aid bone formation and impede the growth of various tumors through oxidative damage to DNA, making it particularly effective for repairing bone defects after osteosarcoma resection [71]. Similarly, strontium peroxide (SrO2)-based OGS also promoted new bone formation and inhibited bone resorption, leveraging the dual role of strontium in bone metabolism by stimulating osteoblasts and inhibiting osteoclasts [72]. Of note, strontium is recognized not only as a bioactive trace element that supports bone health but also as a component of strontium ranelate, a well-known treatment for postmenopausal osteoporosis.

Smart stimulus-responsive biomaterials for bone regeneration: harnessing and modulating ROS

During the inflammatory and regenerative phases of bone healing, ROS concentrations can increase 100-fold [73., 74., 75.]. This induces oxidative stress in leukocytes and platelets [76], and promotes osteoclast formation and bone resorption [77]. Despite the promise of OGS for bone defect treatment, a significant technical challenge is excessive O2 generation that can exacerbate ROS accumulation [60,72,78]. Addressing this challenge when evaluating OGS performance requires rigorous characterization and quantification of in vivo ROS levels. Various probes are available, such as superoxide, H2O2, peroxynitrite, and reactive halogen species [79], and imaging techniques such as two-photon microscopy, ultrasound imaging, and bioluminescence imaging are widely used for noninvasive monitoring of ROS dynamics in vivo [80]. Despite challenges, in vivo imaging of ROS in bone tissues has been achieved by Xie and colleagues [81] using hydro-indocyanine green injections. Ultimately, the goal of these measurements was to engineer ROS-responsive systems [82,83] to counteract oxidative stress and enhance bone healing. For example, catalase (CAT)-encapsulated particles applied to bone, joints, and tumors with high H2O2 content led to O2 production via the reaction: 2H2O2 → 2H2O + O2 (Figure 3) [84., 85., 86.]. Similarly, MnO2-based nanozymes with CAT-like activity decompose H2O2 [87,88]. Based on the differential spatiotemporal needs of O2 in various preclinical models, smart stimulus-responsive materials have been developed that dynamically respond to chemical changes within the bone-healing environment or to external physical cues such as photo-, acoustic, magnetic field, and electrical stimulation [89]. They combine real-time sensing with the ability to trigger the in situ reconfiguration of the BTE scaffold for therapeutic interventions.
Figure 3
  1. Download: Download high-res image (538KB)
  2. Download: Download full-size image

Figure 3. Representative oxidative stimulus-responsive biomaterials.

(A) ROS scavenging and responsive prolonged O2-generating hydrogels (CPP-L/GelMA) which comprise antioxidant enzyme CAT and ROS-responsive O2-releasing nanoparticles (NPs) (PFC@PLGA/PPS), coloaded liposomes (CCP-Ls), and GelMA hydrogels [84]. (B) H2O2-responsive O2-generating PLGA nanoparticles containing CAT for platinum anticancer agent delivery [86]. (C) pH/H2O2-responsive albumin–MnO2 nanoparticles [88]. Image adapted, with permission, from [84,86,88]. Abbreviations: CAT, catalase; CCP-L, co-loaded liposome; GelMA, gelatin methacryloyl; HMCP, human mandible coronoid process; H2O2, hydrogen peroxide; MnO2, manganese dioxide; NPs, nanoparticles; PFC, perfluorocarbon; PLGA, poly(D,L-lactide-co-glycolic acid); PPS, poly(propylene sulfide); [PtLCl]Cl, 4′-bis(pyridine-2-ylmethyl)amino-2-phenylbenzothiazole; ROS, reactive oxygen species.
ROS-responsive scaffolds have been utilized in drug delivery. Lee and colleagues developed poly(D,L-lactide-co-glycolic acid) (PLGA)-based nanoparticles that managed ischemia/reperfusion injuries via the controlled release of heparin and glutathione in response to H2O2, thereby enhancing bone healing [90]. BTE scaffold coatings using ROS-responsive thioketal-based polymers enabled tunable release of BMP-2 and increased bone regeneration by 50% in rat calvarial defect models [91]. Similarly, a titanium implant coated with a hydrogel responsive to ROS in femoral bone defects released Tβ4 to promote macrophage transformation and osteogenic differentiation, and improve bone healing [92]. Aleemardani and colleagues integrated ROS-mitigating quercetin, a plant flavanol, into silk fibroin–CaO2 hydrogel–electrospun fiber–hydrogel constructs [68]. This novel construct released O2 for 12 days while simultaneously performing ROS scavenging, and increased cell viability by 30% under normoxia and by 10% under hypoxia compared with its quercetin-lacking counterpart. Overall, harnessing and modulating ROS within OGS is a promising strategy for overcoming the challenges associated with balancing O2 supply and demand. Continued developments will include fine-tuning their sensitivity to ROS levels. Future research could explore the incorporation of novel O2-sensing elements, chemical bonds, or O2 carriers within the scaffolds. In addition, the long-term stability and biocompatibility of redox-active materials and their byproducts will necessitate further research to fully understand their interactions with the biological environment [93].

OGS optimization via biophysical modeling and machine learning

Biophysical models and machine learning (ML; Box 1) tools have the potential to address some of the challenges associated with the design of OGS by predicting biological responses [94., 95., 96., 97., 98., 99., 100.] to scaffold design and composition [101., 102., 103., 104., 105., 106., 107., 108., 109., 110., 111., 112., 113., 114., 115., 116.]. The ultimate vision of an integrated, iterative workflow involves using ML to analyze biological and imaging data to enhance OGS design and optimize bone regeneration (Figure 4). However, this modeling of complex biological systems remains challenging due to insufficient mechanistic data (e.g., derived from single-cell or spatial transcriptomic approaches [117., 118., 119.]) or phenomenological data (e.g., derived from the imaging strategies described earlier). To date, computational models have been used to predict biological responses and the mechanical properties of scaffolds as a function of their composition, or to enhance the manufacturing process.
Box 1
Fundamental computational and machine learning (ML) models that can be applied to BTE design
ML is a data-driven artificial intelligence (AI) technology that involves creating predictive models by learning from existing data. It operates by fitting models to a small set of 'training data' with known outcomes to develop its predictive capability, and then the models are applied to 'test data' to make generalized predictions about unseen data [128]. ML can be broadly categorized into supervised learning, unsupervised learning, and reinforcement learning. In supervised learning, each training dataset (scaffold geometry, material properties, etc.) is associated with corresponding outputs (permeability, hydrophobicity, etc.) as its labels. Supervised learning learns the input–output relationships from the training data, and therefore is 'supervised' to predict the outputs of the testing data [129].
(i) A decision tree is a supervised learning algorithm for classification and regression. In a tree-like model of decisions, each internal node represents a test on an attribute, each leaf node indicates a class label or decision based on the attributes, and each branch represents the outcome of a test [130].
(ii) A support vector machine (SVM) is a supervised learning method for classification, regression, and outlier detection. SVMs work by identifying the hyperplane that optimally divides a dataset into distinct classes with high robustness and a low risk of overfitting [130].
Unsupervised learning, by contrast, represents the process of training algorithms with unlabeled data to discover inherent patterns within the input. Finally, reinforcement learning is a method that navigates through uncertain environments by learning from past behaviors to select optimal actions that maximize rewards [129].
Neural networks (NNs) are key architectures in ML, inspired by the biological neural network in the human brain, and are applicable to both supervised and unsupervised learning paradigms. An NN consists of layers of interconnected neurons, which include an input layer, one or more hidden layer(s), and an output layer. The neurons and their connections, associated with different weights, allow each layer to perform transformation on inputs with varying signal strength [130].
Ensemble learning in ML trains and combines multiple models to solve specific problems, and improves model performance by aggregating their predictions, which effectively reduces overfitting. Random forest is a notable example of ensemble learning that combines multiple decision trees when each tree is trained on independent input datasets. It is known for its high accuracy, ability to process unbalanced or missing data, and ability to identify important classification variables [130].
Figure 4
  1. Download: Download high-res image (480KB)
  2. Download: Download full-size image

Figure 4. Framework for next-generation oxygen-generating scaffold (OGS) engineering for bone regeneration.

(1) Scaffold development. Based on our biological understanding of the impact of O2 on bone regeneration, OGS design incorporates empirically derived biomaterial properties and advanced biomanufacturing. Biophysical simulations can predict OGS efficiency and optimize O2 generation for customized biological impact. This step represents the preliminary scaffold optimization before in vitro or in vivo experiments. (2) Biological characterization. In vitro and in vivo behaviors of fabricated OGS are monitored via noninvasive imaging techniques to assess the distribution of O2, reactive oxygen species (ROS), and bone regeneration, while also elucidating the role of O2 in bone regeneration. Transcriptomic analysis could provide comprehensive and versatile RNA profiling to elucidate the molecular mechanisms involved. This step provides essential real-world data for subsequent machine learning (ML)-guided optimization. (3) ML-guided scaffold optimizations. Ideally, ML algorithms will use OGS properties (e.g., biomaterial and biomanufacturing aspects) and preclinical and clinical imaging data as inputs to simulate/predict the biological performance of scaffolds in patients and enable a priori optimization of the scaffold design to maximize bone regeneration. The ML models help to identify the most impactful OGS properties for the desired biological outcome through iterative learning and feature significance analyses. This process enables the focused refinement of scaffold design with high efficiency, leading to accelerated next-generation OGS engineering for optimized bone regeneration. Figure created with BioRender.

Predictive modeling of biological responses

Biophysical models applied to cell-seeded PLGA–sodium percarbonate OGS could predict dynamic spatiotemporal O2 interactions during bone healing by accounting for scaffold size, O2 release, and cellular consumption rates in avascular bone defects [94]. Finite element modeling (FEM) has been used to explore how local O2 tension, coupled with the mechanical environment within an implanted osteochondral scaffold, impacted on tissue differentiation and repair outcomes of mesenchymal stem cells (MSCs) [95]. FEM also simulated cellular-level vascular development by mathematically linking O2 sources and Dll4-Notch signaling in endothelial cells [96,97]. At the tissue level, FEM-based approaches have been used to simulate bone regeneration by integrating mechanical stimuli, scaffold degradation, biological responses to VEGF supplementation, vascularization, and O2 delivery [98]. In critical-sized bone defects, spatial variations in these parameters occur during bone regeneration. For accurate FEM simulations, thousands of representative volume elements are needed, resulting in high computational costs. However, the use of neural networks to predict bone regeneration in a sheep tibia model achieved results as accurate as FEM but at one-fifth of the computational cost. The ML model was validated by longitudinal in vivo X-ray images, and showed consistent bone regeneration at 12 month based on remodeling parameters established from the first 9 months of data [99]. Similarly, Ghosh and colleagues used neural networks to substitute tedious FEM analyses in predicting bone growth over macro-textured bone implant surfaces, further demonstrating the efficiency improvements provided by ML compared with traditional FEM methods [100]. Biological parameters for both FEM and ML models are usually sourced from the literature, in vitro/in vivo data, and empirical estimation, which can be updated via iterations and new experimental data [95,99]. ML tools integrate data from in vitro and in vivo experiments, thereby providing a comprehensive framework to elucidate the underlying mechanisms by which spatiotemporal O2 gradients and bone scaffolds enhance osteogenesis, angiogenesis, and overall bone healing. ML approaches can also guide on-demand OGS design to match a desired optimal cellular response. Integrating noninvasive imaging data with ML techniques has the potential to further enhance the simulation of biological processes such as bone regeneration and revascularization. ML models trained on extensive imaging datasets can provide more accurate and robust predictions of the in vivo O2 distribution. Ultimately, these models would be used clinically to facilitate customized interventions tailored to the patient and their defect based on initial imaging data.

Computation-based customization of scaffold material properties

Biophysical simulations and ML also offer powerful solutions to address the technical limitations of OGS designs by predicting the effects of different materials on O2 generation, diffusion, and distribution. For example, O2 gradients within CaO2-based scaffolds supporting the high O2 demand of islet cells were simulated using FEM [101., 102., 103.]. These models, based on oxygenated reaction, diffusion, and consumption kinetics, successfully guided the optimization of cell loading [103]. There are limited studies highlighting the potential of ML for advancing biomaterial designs for BTE applications. In the Polymer Genome Project, computational methods were used to predict polymer properties, such as glass transition temperature and solubility parameters, thereby enhancing material selection and design [104]. ML accelerated peptide discovery for self-assembling hydrogels by combining Monte Carlo tree search and random forest algorithms with molecular dynamics simulations, and efficiently identified sequences with high self-assembly potential [105]. In addition, ML could predict the material properties and bone-forming ability of metals using artificial neural networks (ANNs) [106] or in ceramics [106,107]. ML has been used to predict scaffold biodegradability, which affects O2 release kinetics, tissue integration, and scaffold longevity [108]. Bioceramics and natural substances (e.g., collagen, fibrin), commonly used as OGS binding materials, exhibit complex degradation patterns owing to their inherent heterogeneity and bioactivity. Predicting these patterns using traditional mathematical principles is challenging. However, random forest algorithms have been trained on histological images to predict collagen-based scaffold biodegradation [109]. In addition, neural networks trained on the relationships between mechanical properties, crosslinking, and swelling characteristics and degradation rates have predicted gelatin-based bone scaffold degradation [110].

Computation-based customization of advanced biomanufacturing

Finally, ML has been used to create bone scaffolds with controlled porosity, Young's modulus, and compression strength, and achieved <5% error from design targets [111,112]. Although integrating FEM and computational fluid dynamics (CFD) helps to simulate mechanical properties and O2 permeability concurrently to iteratively adjust scaffold design [113,114], recent advances have utilized a more efficient ML approach for inverse design, such as supervised learning to generate Voronoi lattices with targeted mechanical properties that were validated through compression tests and simulations [112]. Another study used FEM simulations to create a neural network dataset linking microstructural parameters to the elastic matrix of bone [111]. ML also provides significant benefits to scaffold accuracy and consistency in image-based computer-aided design (CAD) and additive manufacturing (AM)-based biomanufacturing. Creating patient-specific bone scaffolds relies on accurate bone segmentation from computerized tomography (CT) scans. Minnema and colleagues improved this step by using automated convolutional neural networks (CNNs), and achieved a 92% similarity to gold standard segmentation tools in shorter times [115]. During the subsequent bioprinting process, ML models analyzed data signatures from in situ IR sensors in a feedforward process to predict print quality metrics and adjust process parameters in real time, thereby preventing manufacturing flaws and improving print quality [116].
Challenges to ML for OGS include the need for extensive, high-quality data for model training. The variability in the experimental conditions and requirements for preclinical versus clinical applications require collaborative efforts to build comprehensive datasets for OGS materials. Currently, there are limited publicly available biomaterial databases, and most groups still rely on extracting descriptors from multiple sources, which reduces scalability and increases computing power requirements [120]. High-quality, standardized data collection protocols will be needed in the future to ensure accurate and reproducible experimental data. In addition, the 'black-box' nature of many ML models limits the interpretability of ML-generated designs for in vivo applications. Therefore, ML should be considered as an aid, and not as a replacement for empirical trials, to be accompanied by rigorous experimental validation of OGS performance. Enhancing the transparency of ML algorithms through explainable artificial intelligence (AI) techniques will increase their clinical acceptance among medical professionals [121]. Creating effective ML models for OGS requires interdisciplinary expertise in materials science, tissue engineering, and computer science, which would further advance bone scaffold and next-generation OGS design.

Concluding remarks and future perspectives

Overall, computational and ML modeling can predict O2 release, scaffold degradation, and multilevel biological responses to spatiotemporal O2 gradients in bone defects. Integrating these predictions into a collective framework would enhance OGS design and encourage the development of more standardized manufacturing for OGS applications in patients. However, critical-sized craniomaxillofacial bone injuries in patients are more varied and complex than the controlled, predictable injuries in preclinical models. Consequently, the need for varied OGS precludes a one-size-fits-all solution. A major advance would be simulations that identify suitable OGS designs for optimal bone regeneration. Predictions of the design of an OGS for adequate and effective treatment would need to account for multiple parameters including: (i) the number and phenotype of transplanted osteoprogenitors and the impact of pO2 on their survival, proliferation, and differentiation. (ii) The anatomical location of the defect site: the prevailing pO2 and healing outcomes are affected by soft-tissue coverage that facilitates angiogenesis, and the interface with the adjacent host bone that impacts on efficient osteoconduction. (iii) The geometry of the injury: optimizing the shape, size, microstructure, and pore architecture of the scaffold is crucial because it directly influences the mechanical properties as well as the rate and duration of O2 release. Despite advances that enable computer simulations of animal models, a significant barrier to the application of ML or computational models for simulating de novo bone regeneration in patients is the lack of meaningful clinical data (see Outstanding questions). Extrapolating healing dynamics to patients from imaging or transcriptomic data obtained in preclinical models is challenging. Thus, future studies that include longitudinal imaging data on bone healing and revascularization in patients would greatly advance the clinical potential of OGS.
Outstanding questions
What are the physiologically relevant O2 levels for bone healing, given its spatiotemporal variations in vivo?
What innovations are necessary to develop or translate noninvasive imaging techniques for real-time monitoring of in vivo O2 and ROS levels, with the goal of optimizing scaffold performance and successful bone regeneration?
How can ML facilitate the optimization of OGS properties and design for efficacious bone regeneration?
How can we better integrate and leverage ML models in tissue engineering?

Acknowledgments

This work was supported by funding from the National Institute of Dental and Craniofacial Research (NIDCR; grant 1R01DE027957), the Maryland Stem Cell Research Fund (2022-MSCRFV-5782), and the National Cancer Institute (NCI; 5R01CA237597-05, 5R01CA196701-07, and 5R01CA237597-05).

Declaration of interests

W.L.G. hold shares in EpiBone, Inc. The other authors declare no competing interests.

References

Cited by (0)

Glossary

Additive manufacturing (AM)
the process of creating objects by adding material layer by layer based on a digital model.
Allografts
grafts transplanted into a recipient from a genetically non-identical donor of the same species.
Autografts
grafts harvested from and transferred to another site in the same individual of a species as the 'gold standard' for grafting.
Computational fluid dynamics (CFD)
a computational technique that analyzes and predicts fluid flow behaviors in complex systems, often based on Navier–Stokes equations.
Computer-aided design (CAD)
detailed designs and technical drawings of objects generated by computational software to visualize, modify, analyze, and optimize designs before physical production.
Cranial window
a surgical setup in which the skull is thinned or replaced with a synthetic optical interface such as a glass coverslip, allowing optical access of a microscope.
Critical-sized bone defects
bone defects that do not heal spontaneously without medical intervention to facilitate bone healing.
Finite element modeling (FEM)
a numerical technique that simulates and analyze complex structures or systems by dividing them to smaller elements and solving differential equations to predict their responses to applied loads or stimuli.
Hypoxia
tissue O2 levels that are insufficient to maintain adequate homeostasis, which can lead to impaired cellular function, reduced tissue regeneration, increased inflammation, and activation of pathways that can trigger adaptive responses.
Ischemia
shortage of O2 and nutrients needed for cellular metabolism caused by restriction in the blood supply.
Oxygen saturation (SO2)
the ratio of oxyhemoglobin concentration (HbO2) to the concentration of total functional hemoglobin (HbT), which indicates the utilization of the present O2 transport capacity.
Platinum porphyrin-coumarin-343 (PtP-C343)
a phosphorescent O2 probe that can be safely delivered into the body for tissue O2 measurements and imaging.
Reactive oxygen species (ROS)
highly reactive oxygen-containing molecules such as superoxide radicals, hydrogen peroxide, and hydroxyl radicals. ROS are essential for cellular signaling and immune functions, but excessive ROS can cause oxidative damage, contributing to metabolic diseases and cell aging.
Representative volume elements
small representative sections of a structure used in computational modeling for property analysis which capture essential features while minimizing computational complexity.
View Abstract