Next Article in Journal
Measurement of Hepatic CYP3A4 and 2D6 Activity Using Radioiodine-Labeled O-Desmethylvenlafaxine
Previous Article in Journal
Plasticity towards Rigidity: A Macrophage Conundrum in Pulmonary Fibrosis
Previous Article in Special Issue
Biological Response to Bioinspired Microporous 3D-Printed Scaffolds for Bone Tissue Engineering
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 23
Issue 19
10.3390/ijms231911459
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Macropore Regulation of Hydroxyapatite Osteoinduction via Microfluidic Pathway

by
Feng Shi
1,2,†,
Xin Fang
1,†,
Teng Zhou
1,
Xu Huang
1,
Ke Duan
1,
Jianxin Wang
1,
Shuxin Qu
1,
Wei Zhi
1,* and
Jie Weng
1,*
1
Key Laboratory of Advanced Technologies of Materials (Ministry of Education), School of Materials Science and Engineering, Southwest Jiaotong University, Chengdu 610031, China
2
Collaboration and Innovation Center of Tissue Repair Material Engineering Technology, College of Life Science, China West Normal University, Nanchong 637009, China
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Int. J. Mol. Sci. 2022, 23(19), 11459; https://doi.org/10.3390/ijms231911459
Submission received: 10 August 2022 / Revised: 18 September 2022 / Accepted: 26 September 2022 / Published: 28 September 2022
(This article belongs to the Special Issue Novel Bioinspired Materials for Biomedical and Biotechnological Applications)
Browse Figures
Versions Notes

Abstract

:
Macroporous characteristics have been shown to play a key role in the osteoinductivity of hydroxyapatite ceramics, but the physics underlying the new bone formation and distribution in such scaffolds still remain elusive. The work here has emphasized the osteoinductive capacity of porous hydroxyapatite scaffolds containing different macroporous sizes (200–400 μm, 1200–1500 μm) and geometries (star shape, spherical shape). The assumption is that both the size and shape of a macropore structure may affect the microfluidic pathways in the scaffolds, which results in the different bone formations and distribution. Herein, a mathematical model and an animal experiment were proposed to support this hypothesis. The results showed that the porous scaffolds with the spherical macropores and large pore sizes (1200–1500 μm) had higher new bone production and more uniform new bone distribution than others. A finite element analysis suggested that the macropore shape affected the distribution of the medium–high velocity flow field, while the macropore size effected microfluid speed and the value of the shear stress in the scaffolds. Additionally, the result of scaffolds implanted into the dorsal muscle having a higher new bone mass than the abdominal cavity suggested that the mechanical load of the host tissue could play a key role in the microfluidic pathway mechanism. All these findings suggested that the osteoinduction of these scaffolds depends on both the microfluid velocity and shear stress generated by the macropore size and shape. This study, therefore, provides new insights into the inherent osteoinductive mechanisms of bioceramics, and may offer clues toward a rational design of bioceramic scaffolds with improved osteoinductivity.

1. Introduction

Osteoinductive biomaterials are promising bone grafting materials in orthopedic clinics. Among these materials, calcium phosphate (CaP) bioceramics have received extensive interest because of their chemical similarity to bone minerals [1,2,3,4,5,6,7]. In addition to chemical composition, the macro-pore structure is another important factor in the bone-inducing potential of biomaterials [8]. Pore penetration and macro-pore morphology, such as pore size [9,10], geometry [11,12], orientation, porosity, interconnected pore size, and so on, determine the growth of the blood vessels [8,13,14], and the distribution of nutrients, oxygen, cells and tissues in the scaffolds, as well as the total surface area available for material–cell interactions [15,16,17,18]. Earlier studies have found that porous materials with a pore size of ≥400 μm and the interconnective pore diameter of ≥150 μm, provide an ideal space for bone tissue growth [19,20,21]. However, the mechanisms underlying bone formation and distribution in such scaffolds remain elusive, and guidelines for the optimization of products are lacking.
In vitro bioreactor studies found that the microfluid dynamic environment (mainly including microfluid velocity and fluid shear stress), formed in scaffolds after culture medium perfusion, has an important effect on the cells’ distribution, growth, and differentiation on these scaffolds [22,23,24]. The main factors affecting the microfluid dynamic environment include the perfusion velocity of the fluid medium to the scaffolds and the pore structure of the scaffolds [25,26]. That means that when scaffolds are implanted in vivo, due to the impact of the stress load of implantation sites, body fluid as a fluid medium will form a microfluid dynamic environment in the scaffolds. Therefore, the macro-pore structure properties of the materials are transformed into a biological response regulating signals, and play a role in the osteogenesis. However, due to the lack of suitable micro-sensors, the local microfluid velocity and fluid shear force in the scaffolds cannot be measured [25,27,28]. In this case, it is feasible to study the correlation between the material pore structure, the microfluid environment and the bone formation, with computational fluid dynamics (CFD), and finally establish an accurate quantitative relationship.
In this study, two kinds of hydroxyapatite (HA) scaffolds, hydroxyapatite spheres accumulating scaffolds (HASAs) and hydroxyapatite porogen-pore scaffolds (HAPPs), with different macro-pore shapes and size distributions were fabricated with the sol-gel method. The two kinds of scaffolds had star-shaped pores and spherical pores, respectively, good interconnectivity and similar microstructures. After the two kinds of scaffolds were characterized, they were implanted in the dorsal muscle and abdominal cavities of adult dogs, to compare the effects of different macro-pore shapes and/or different macro-pore sizes on the ectopic bone formation ability of the scaffolds. After implantation for 1 month, 3 months and 6 months, the scaffolds were analyzed for histology, histomorphometry and biomechanics. Finally, mathematical models for 3D scaffolds of the HASAs and HAPPS were developed, and a finite element analysis (FEA) was analyzed for the microfluid dynamic environment using computational fluid dynamics (CFD). The results confirmed the hypothesis that both the size and shape of the macropores affected the microfluidic pathways in the scaffolds, and these may underlie the bone formation and distribution in these scaffolds. The goal was to better understand the role of microfluidic environments in regulating the osteoinduction of macroporous CaP ceramics.

2. Results

2.1. The Macro-Pore Structure Properties of HASAs and HAPPs

HASAs and HAPPs are porous cylinders with a height of 2–2.5 cm and a diameter of 1.5 cm. The macro-pore structures of all the scaffolds were prepared and characterized as reported previously [29]. The HAPPs had higher total porosity and macro-porosity than the HASAs, while the two kinds of scaffolds had similar micro-porosities. Micro-CT reconstruction found the macro-pore shapes of the HASAs to be stellate-shaped, which were formed by the HA spheres accumulating (Figure 1a). That is to say, the macro-pore structure properties of the HASAs were determined by the size of the HA spheres and the arrangement of the adjacent spheres. The arrangements of the S-HASAs (Figure 1b) and L-HASAs (Figure 1c) were depicted by three-ball accumulations, four-ball accumulations and five-ball accumulations (Figure 1d). The minimum macro-pore size existed in the three-ball accumulations, and the maximum macro-pore size was smaller than the diameter of the HA sphere. The micro-CT reconstruction showed the macro-pore shapes of the HAPPs were spherical-shaped (Figure 1f). Each macro-pore had two or more inter-connective pores. The ratio of the sizes of inter-connective pore/macro-pores was ~0.4 (Figure 1g,h). The porous structures of all kinds of scaffolds with excellent inter-connectivity were beneficial to humoral flow, protein adhesion, cell migration and tissue ingrowth.

2.2. General Conditions of Animals

Within 3 days after the operation, the animals resumed activities and normal water/food intake. The surgical wounds healed normally. Three weeks after the operation, there was no swelling, inflammation, nor tissue necrosis around the incisions.

2.3. The Histological Characterization of Ectopic Bone Formation and Blood Vessel Ingrowth

The hematoxylin and eosin (H&E) staining results (Figure 2 and Figure 3) showed the ectopic bone formation and distribution in the scaffolds. One month after implantation, no new bone tissue was observed in the L-HASAs, S-HASAs and S-HAPPs, neither in the dorsal muscle nor in the abdominal cavity (Figure 2a1,b1 and Figure 3a1,b1). However, primary callus was widely distributed in the central and outer areas of the L-HAPPs (Figure 2c1 and Figure 3c1). Even a small amount of new bone tissue appeared in the central area of the L-HAPPs. Three months after implantation, new bone tissue was seen in the outer area of the L-HASAs with a colony and diffused distribution both in dorsal muscles and in the abdominal cavities (Figure 2a2 and Figure 3a2), while new bone tissues were only sporadically present in the outer areas of the S-HASAs (Figure 2b2 and Figure 3b2). New bone tissue was also evident in the central and outer areas of the L-HAPPs (Figure 2c2 and Figure 3c2), while only a small amount of callus was observed in the central and outer areas of the S-HAPPs (Figure 2d2 and Figure 3d2). Six months after implantation, more new bone tissues gradually grew into the central area of the L-HASAs, showing an increased bone maturation and the appearance of a medullary cavity (Figure 2a3 and Figure 3a3). There was a similar trend of new bone formation in the S-HASAs, but with fewer bone tissues and an absence of a medullary cavity (Figure 2b3 and Figure 3b3). Meanwhile, new bone tissues with medullary cavities widely appeared in the central and outer areas of the L-HAPPs (Figure 2c3 and Figure 3c3), whereas a small amount of new bone tissues was observed in the S-HAPPs (Figure 2d3 and Figure 3d3). In comparison, at each time point, the scaffolds implanted in the abdominal cavities exhibited less and later bone formation (vs. in dorsal muscles).
A histomorphometric measurement (Figure 4a,b) found that in the dorsal muscle implantations, regardless of the kind of scaffolds, there was more active new bone formation and more new bone tissue quantity in the L-HASAs and the L-HAPPs than those in the S-HASAs and the S-HAPPs. Interestingly, the new bone formation and new bone tissue quantities of the L-HAPPs were significantly higher than those of the L-HASAs, while those of the S-HAPPs were significantly lower than those of the S-HASAs from 1 month to 6 months’ implantation. In addition, a similar trend was observed in the abdominal cavity implantations.
One month after implantation, new bone formation was only observed in the L-HAPPs, both in the dorsal muscle (4.27 ± 1.85%) and in the abdominal cavity (3.06 ± 0.94%). Six months after implantation, the ratio of new bone formation in the L-HAPPs was 28.31 ± 3.06% in the dorsal muscle and 15.93 ± 3.87% in the abdominal cavity. In comparison, for the L-HASAs, new bone formation was 22.27 ± 2.89% in the dorsal muscle and 10.28 ± 2.06% in the abdominal cavity (Figure 4a,b). For the S-HAPPs, it was 7.61 ± 2.48% in the dorsal muscle and 3.94 ± 1.33% in the abdominal cavity. For the S-HASAs, it was 14.0 ± 3.67% in the dorsal muscle and 5.65 ± 1.17% in the abdominal cavity (Figure 4a,b). Throughout this study, the S-HAPPs exhibited inferior osteoinduction to the other three groups.
The measured blood vessel density in the groups at the early implantation stage showed no significant difference for vascularization of all the scaffolds. Subsequently, the new vessel densities in the four groups followed similar trends to the new bone formations (Figure 4c,d). The amount of vessels/mm2 in the L-HAPPs group (6.44 ± 0.75 mm−2 in the dorsal muscle and 4.38 ± 0.53 mm−2 in the abdominal cavity) was significantly higher than the L-HASAs group (4.98 ± 0.4 mm−2 in the dorsal muscle and 3.17 ± 0.51 mm−2 in the abdominal cavity) at 6 months’ implantation. Six months after implantation, the density in the S-HASAs group (4.22 ± 0.37 mm−2 in the dorsal muscle and 2.87 ± 0.39 mm−2 in the abdominal cavity) was significantly higher than that in the S-HAPPs (3.11 ± 0.74 mm−2 in the dorsal muscle and 1.9 ± 0.63 mm−2 in the abdominal cavity).
These results showed that new bone and vessel formation was better in the L-HAPPs and L-HASAs than those in the S-HAPPs and S-HASAs, as the implantation time extended. In summary (Table 1), the new bone formation in the L-HASAs and L-HAPPs after being implanted in the dorsal muscle and the abdominal cavity were better than those of the S-HASAs and S-HAPPs, respectively. The scaffolds implanted in the abdominal cavities experienced less and later bone formation at each time point than those in the dorsal muscles. Moreover, the new bone distribution in the HASAs mainly appeared in the outer area of the scaffolds, while those in HAPPs mostly existed both in the central and outer areas of the scaffolds since the early stage.

2.4. Compressive Strength of Reconstructed Scaffolds

Figure 5 showed L-HASAs and S-HASAs had a similar strength before implantation, while the L-HAPPs had a significantly lower compressive strength than the S-HAPPs. Following the ingrowth and filling of tissue, and the increase in new bone formation in the scaffolds, the compressive strengths of all scaffolds increased at different degrees. The L-HASAs and L-HAPPs increased more prominently. Six months after implantation, the compressive strength of the L-HASAs implanted in the dorsal muscle and the abdominal cavity increased approximately 16-fold (from 0.17 MPa to 2.73 MPa) and 13.3-fold (from 0.17 MPa to 2.26 MPa), respectively. In comparison, the L-HAPPs implanted in the dorsal muscle and the abdominal cavity increased almost 4-fold (from 0.63 MPa to 2.51 MPa) and 3.2-fold (from 0.63 MPa to 2.02 MPa), respectively. The S-HASAs increased almost 12.5-fold (from 0.21 MPa to 2.63 MPa) and 3.2-fold (from 0.21 MPa to 2.04 MPa), respectively. The S-HAPPs increased almost 2.3-fold (from 0.89 MPa to 2.03 MPa) and 1.5-fold (from 0.89 MPa to 1.33 MPa), respectively. Compared with the other three groups, the S-HAPPs group was obviously lower in compressive strength because of the lesser amount and maturity of new bone formations at 6 months.
Additionally, the results showed that the compressive strengths of the scaffolds implanted in the dorsal muscle at three time points were higher than of those implanted in the abdominal cavity, which were consistent with the histomorphometric findings, and indirectly showed that the new bone formation in the dorsal muscle was better than in the abdominal cavity.

2.5. Numerical Simulation of Fluid Flow in the Scaffolds

A CFD simulation showed that the flow through the scaffold’s microarchitectures was highly non-uniform (Figure 6), the highest flow speeds in the HASAs (star-shape macro-pores) appeared near the edges (red area of Figure 6a1,b1), while the lowest speeds were at the central region (blue area of Figure 6a1,b1). The distribution of local shear stresses on the whole fluid pathway surfaces of the HASAs is shown in Figure 6a2,b2. The highest values for local shear stresses were found on the walls of pores close to the edge of the scaffold (red area of Figure 6a2,b2). In the L-HASAs, the highest surface shear stress in the range from ≈ 1.02 to 1.15 Pa correlated to the highest flow speed from ≈ 0.039 to 0.044 m/s; in the S-HASAs, the highest surface shear stress in the range from ≈ 4.50 to 5.07 Pa correlated to the highest flow speed from ≈ 0.053 to 0.059 m/s. For the HAPPs (sphere-shaped macro-pores), the highest flow speeds were found in the middle of the pores in the central region (yellow or green area of Figure 6c1,d1), while the lowest speeds occurred near the surface of the pores (blue area of Figure 6c1,d1). The flow speeds were also high near the interconnection of the pores. The distribution of local shear stresses on the whole fluid pathway surfaces of the HASAs is shown in Figure 6c2,d2. The highest values for the local shear stresses were found on the surface of the interconnection of the pores (red area of Figure 6c2,d2). In the L-HAPPs, a surface shear stress around an interconnected pore in the range from ≈1.80 to 0.90 Pa correlated to a flow speed from ≈0.025 to 0.034 m/s; in the S-HAPPs, a surface shear stress around an interconnected pore in the range from ≈3.91 to 5.87 Pa correlated to a flow speed from ≈0.038 to 0.044 m/s.
The results suggest that the macro-pore shape of the scaffold affects the distribution of the fluid field and fluid shear stress, while the macro-pore size of the scaffold affects the flow speed and the value of shear stress.

3. Discussion

In this study, four kinds of porous HA scaffolds with the same chemical phase composition and micro/nano structures [29,30,31], but different macro-pore structures (size, shape) were used to ensure the material variables caused by the macro-pore structural variables. The HASAs and HAPPs had two sizes of macro-pore, 0.2–0.4 mm and 1.2–1.5 mm, respectively, which can be seen as two ranges of macro-pore size, small size and large size. Additionally, other studies have reported that the macro-pores of an HA porous scaffold with inter-connectivity network pores with the size of 140–160 μm could be completely filled by bone tissue after implantation in the orthotopic bone tissue [32,33,34,35]. Therefore, the HAPPs used in this study were selected to an inter-connective pore/macro-pore size ratio of ~0.4, to avoid the effect of the size of inter-connective pores.
The results of the histological analysis indicated that the pore structure of the HA porous scaffolds was more beneficial for the growth of cells and tissues, and the formation of new bone, when the size of the macro-pores were >400 μm (Figure 2 and Figure 3). Previous research has found that porous α-TCP bone cement columns with the 150 μm macro-pore size had no new bone formation after 90–180 days’ subcutaneous implantation in goats, and lost the integrity of the scaffolds after 3 months’ implantation [16]. Hence, the small macro-pore size was considered not conducive to new bone formation, and it was believed that the loss of the macro-pore structure of scaffolds would impede the transport of nutrition in the scaffolds, and be not conducive to protein adsorption, cell adhesion and differentiation [16,36]. Although all of the kinds of HA scaffolds had good connectivity and high porosity, the osteoinduction of the L-HASAs and L-HAPPs was obviously superior to the S-HASAs and S-HAPPs. This may be due to the fact that the scaffolds with macro-pore sizes of <400 μm, especially for the HAPPs, are likely to allow the blood or tissue fluid generated during implantation to be trapped in the macro-pores, resulting in a large number of cells and proteins blocking the inter-connective pores of the scaffold, thereby reducing scaffold connectivity. This is not conducive for the vascular network of host tissue to grow into the scaffolds and to provide more adequate oxygen and nutritional support for the cells and tissues grown into the scaffold. This indicates that there may be a threshold for the macro-pore size of the HA porous scaffold to determine osteoinduction. Only when the size of a macro-pore exceeds the threshold, the scaffold would have a strong osteogenic potential.
In this study, we found that the shape of the macro-pores affected the distribution of new bone tissue in the HA porous scaffolds. The bone formation in the HASAs first appeared in the outer region of the scaffolds, and gradually grew into the central region as implantation time was extended. The bone formation in the HAPPs appeared in the central and outer regions simultaneously, and was particularly active in the central region. A previous study has reported that bone formation was observed in the concave of the graft, and had not appeared in the protruding part of the graft when rod- and disc-shaped HA scaffolds with concave surfaces were implanted into the muscle of a baboon [37]. Other studies have reported that the bone formation of porous ceramic scaffolds mainly occurred in the pore’s inner wall, near the opening of the pore after an ectopic implantation in vivo [38,39,40]. Recently, Zhang et al. has investigated the osteogenic properties of triply periodic minimum surfaces structure-based HA scaffolds, and also found that the amount of new bone formation and the distribution area above the scaffolds were significantly different from those of the HA scaffolds with cross-hatch structures [41]. Those results have indicated the influence of the shape of the macro-pores on the osteoinduction of materials [37,42]. The porous gradient distribution was more conducive to osteogenesis and maintained better mechanical properties than the uniform distribution [43]. These data should be related to the macro-pore structure affecting the distribution of blood or liquid fluid in the scaffolds, thus affecting the protein and cell adhesion, nutrient and oxygen transport.
The results of this study have also further revealed that the shape and the size of macro-pores play a synergistic effect, affecting the growth of new bone in the scaffolds. When the size is about 1300–1500 μm, spherical macro-pores are likely to be more conducive to osteogenesis, but when the size is <400 μm, stellate macro-pores may be more conducive to osteogenesis. The reason for this phenomenon may be that there exists a threshold of macro-pore size for the concave curved surface of the spherical pore structure of the HAPPs. When the pore size is below this threshold, the oxygen supply, nutrients diffusion, cell migration, and tissue growth will be hindered after a scaffold implantation. In contrast, the convex curved surface of the stellate pore structure of the HASAs was more conducive to liquid flow, protein adhesion and cell migration than the HAPPs when the size of the macro-pore was <400 μm. However, when the pore size exceeded the threshold, the spherical pore structure was more conducive to tissue growth and new bone growth, as the HAPPs have a greater interconnection of space. Hence, the L-HAPPs had better osteogenic ability than the L-HASAs, while the S-HASAs had better osteogenic ability than the S-HAPPs.
In addition, the stress environments of the implanted sites in vivo also significantly affected and regulated the process of bone formation. Under a higher stress in the dorsal muscle implantation, compared with those in the abdominal cavity, the scaffolds gradually changed to spindle-like shapes, and the new bone formation was significantly higher than those implanted in the abdominal cavity. Compared with the abdominal cavity, muscles are more commonly selected in osteoinduction research. The muscle tissue has adequate blood supply, rich interstitial cells, and a compressive stress environment following a rhythmic contraction–diastole, which is conducive to cell, oxygen and nutrients, followed by tissue fluid or granulation tissue into the scaffolds. The pressure loading on the scaffold can accelerate the degradation of the calcium phosphate ceramic material, thus forming a locally high concentration of calcium and phosphorus ions in the environment around the scaffold. Calcium and phosphorus ions formed a mineralized layer on the surface of the macro-pore wall of the scaffold, which provided a favorable environment for osteogenesis.
The biomechanical results reflected the regulatory effect of the macro-pore structure of the HA scaffolds and the stress load of the implantation sites on the new bone formation and tissue ingrowth. Three months after implantation, the compressive strength of the L-HAPPs and L-HASAs became significantly higher than that of the S-HAPPs and S-HASAs, respectively. Three and six months after implantation, tissue ingrowth and new bone formation were observed in the L-HASAs and L-HAPPs. For the L-HASAs, the HA spheres were encapsulated by new bone tissue and connective tissue as a whole, and the movement between the HA spheres was restricted, which further enhanced the compressive strength of the graft. For the L-HAPPs, excellent tissue ingrowth and wide reconstruction in the macro-pores’ walls increased the compressive strength significantly. On the other hand, since the callus formation was present in the L-HAPPs at the early implantation stage, the compressive strength of the L-HAPPs was significantly higher than that of the L-HASAs after 1 month of implantation. With the implantation time increased, the compressive strength of the S-HASAs was promoted by new bone formation, mainly in the outer region of the graft and connective tissue ingrowth. The compressive strength of the S-HAPPs was increased only moderately because of the limited bone formation and the connective tissue ingrowth. Hence, the compressive strength of the S-HASAs was substantially higher than that of the S-HAPPs 3 months after implantation. In addition, the compressive strength of the scaffolds in the dorsal muscle was significantly higher than that in the abdominal cavity at each implanted time point. The stress load of the scaffold implanted into the dorsal muscle was higher than that of the abdominal cavity, and the scaffold was located in the internal environment rich in blood supply and mesenchymal cells. Therefore, the new bone formation and tissue ingrowth of grafts implanted in the dorsal muscle were better than those implanted in the abdominal cavity. The change of the mechanical properties after the implantation of the scaffold has an important impact on its application. Compared with incorporated composite gels, which are commonly used as bone tissue engineering scaffolds to fill bone defects in non-load-bearing sites, calcium phosphate bioceramics have the potential to be applied in the repairment of load-bearing bone defects [44,45]. The important influence of the shape and size of the macropores on the change of the mechanical strength in vivo undoubtedly makes it a key consideration in the optimization of the design of the scaffolds.
The difference of new bone formation and distribution between the HASAs and HAPPs was explained by the results of the flow simulation calculation. In the microfluid environment of the HASAs, the medium–high velocity flow field was mainly distributed in the outer region. The medium–high velocity of the fluid was conducive to the transmission and diffusion of oxygen and nutrients, the metabolites exchange, the cell migration, and tissue ingrowth. Therefore, in the HASAs, the new bone formation always appeared first in the outer region, and gradually grew into the central region. Furthermore, the difference of the medium–high velocity fluids suggested that there existed a better osteogenic environment in the L-HASAs than in the S-HASAs. While the medium–high velocity flow field of the HAPPs was uniformly distributed in the whole scaffold, and the highest flow velocity was located in the central region, which was beneficial to the uniform distribution of oxygen and nutrients in the scaffolds, and also beneficial for cell migration, survival and tissue ingrowth. Hence the ectopic bone formation often occurred simultaneously in the central region and the outer region of the HAPPs. At the same time, the L-HAPPs were superior to the S-HAPPS, both in the distribution of the medium–high velocity flows and the difference of flow velocities, suggesting that the L-HAPPs had a better osteogenic environment.
In this study, we found that the size and shape of the macro-pores had an obvious effect on the new bone formation and distribution in the HA scaffolds. The flow simulation results showed that this effect was most likely to work through the microfluid environment of the scaffold. That is to say, the microfluid dynamic environment acts as a transition pathway linking the materials’ macrostructure properties to signals regulating the osteoinduction of these materials. Xu et al. have reported that the fluid dynamic environment and macrostructure properties had important effects on the cell growth and distribution in the scaffolds [25]. When the perfusion of porous scaffolds with random macrostructures at a certain velocity continued, the local proliferation rates and distributions of cells on the scaffold were different [25]. Another study has shown that the distribution of the fluid dynamic environment was highly dependent on the macro-pore distribution in the scaffold, and the different macro-pore structures provided different fluid accessibility [46]. The study has also shown that the process of the cells’ differentiation in the fluid dynamic environment was more regulated by the axial perfusion of the scaffold [47]. Porter et al. have reported that the microfluid dynamic environment can enhance the transport of oxygen and nutrients, and the exclusion of metabolites, in in vitro tissue engineering constructs, but have also provided fluid stress to stimulate cell proliferation and differentiation, while the microfluid dynamic environment was determined by the macro-porosity, the size of the macro-pores and the fluid velocity [26]. These studies have suggested close correlations among the scaffold macro-porous structures, microfluid dynamic environments and cell distribution and differentiation.
Although the results of the finite element simulation can well explain the histological results in the in vivo experiment, it was hard to completely match the simulation and animal experiment. In the process of the scaffold preparation, influenced by the process parameters and preparation methods, there were some individual differences and structural imperfections in the obtained scaffolds, thus affecting the in vivo experimental results. In the simulation calculation, the macro-pore structures of the scaffolds were set to complete interconnectivity, the sizes of the macro-pores (in the HAPPs) or the sizes of the spheres (in the HASAs) were set to a uniform value, and the ratio of the inter-connective pore size to the macro-pore size was set to 0.4. However, in the actual scaffold preparation, the sizes of pores or spheres are difficult to be completely consistent, and the ratio of the inter-connective pore sizes to the macro-pore sizes was also difficult to control in an exact value. Although the parameters of the macro-pore structure tend to theoretical configuration, the difference is inevitable.
In addition, in our finite element simulation, only one velocity-inlet value was set. However, the stress loading on the scaffold implanted in the different sites in vivo is not the same, which will affect the velocity-inlet of the fluid flow across the scaffold, and hence affect the microfluid dynamic environment of the scaffold [25]. This was demonstrated in our in vivo experimental results, the new bone formation in the scaffolds implanted in the abdominal cavity with lower stress loading was generally worse than those implanted in the dorsal muscle with higher stress loading.
Bone formation in bioceramic scaffolds is considered a dynamic equilibrium process, involving the growth of blood vessels, supersaturation of calcium ion concentrations, deposition of bioactive apatite layers, and migration, adhesion, proliferation and differentiation of cells, which are necessary for the formation of new bone [48,49,50,51,52]. On this basis, combined with the results of this study, we explained the mechanism of the formation and distribution of the new bone tissues, which were influenced by the microfluid dynamic environment (Figure 7). At the early implantation stage, blood, tissue fluid, mesenchymal cells in the surrounding tissue, and inflammatory factors were driven by the stress load of the implant site into the scaffold. In this process, the microfluid dynamic environment (fluid velocity, flow field distribution, and fluid stress distribution, etc.) within the scaffold determined the distribution of the activity factors and the cells. Subsequently, in the processes of the blood scab formation, calcium phosphate dissolution–redeposition, bone-like apatite formation, protein adhesion, granulation tissue formation and ingrowth, and the microfluid dynamic environment determined the deposition volume and distribution of the bone-like apatite and the blood vessels’ distribution. Thereafter, the mesenchymal cells recruited from the surrounding tissue and blood vessels adhered to the bone-like apatite layer, the cells proliferated and osteogenesis differentiated under chemical, surface micro-morphology and mechanical stimulation. Osteoclasts appeared, fibrous callus formed, and finally new bone appeared afterwards. The microfluid dynamic environment affected the distribution of nutrients and oxygen, the metabolites excretion, the cell adhesion and differentiation, as well as the formation and distribution of new bone in this sequence of events. Throughout these processes, the microfluid dynamic environment was determined by the size and shape of the macro-pore and the stress load at the implant site.

4. Materials and Methods

4.1. HA Scaffolds Preparation

The spherical HA powder used in this experiment was from the National Biomedical Materials Engineering Research Center of Sichuan University. Its main particle size was about 20 μm.

4.1.1. Two Macro-Pore Sizes of HASAs Preparation

The HASAs were prepared by accumulating the HA spheres in porous HA tubes, using the method reported in our group [29,31,53,54]. Briefly, the HA powders were dispersed in a chitin-lithium chloride solution in dimethylformamide (all from Chengdu Kelong Chemical Co., Ltd., Chengdu, China) to form a suspension. The suspension was shaped into spheres, porous tubes, and porous discs, and subsequently sintered into ceramic objects. The sintered HA ceramic tubes and discs were reinforced by immersing in polylactide (PLA, Mw 120 K, Institute of Organic Chemistry, Chinese Academy of Sciences, Chengdu, China)/acetone solution (5% w/w). The HA spheres (diameter: 280–430 μm, 1300–1500 μm) were filled, respectively, into the porous HA tubes (external diameter: 15 mm, internal diameter: 10 mm, height: 20 mm), gently shaken to ensure the spheres were tightly packed, and finally sealed with the porous HA discs (diameter: 15 mm, height: 3 mm) to the final scaffolds (Figure 1). The HASAs with two sphere sizes were labeled the L-HASAs and the S-HASAs.

4.1.2. Two Macro-Pore Sizes of HAPPs Preparation

The HAPPs were prepared by using paraffin spheres with diameters of about 3 mm and about 0.8 mm as the pore-forming agents, respectively, according to the preparation method of our earlier publication [30]. Briefly, HA powders were dispersed in a chitin-lithium chloride solution in dimethylformamide (all from Chengdu Kelong Chemical Co., Ltd., China) to form a suspension. Subsequently, the melting paraffin was dropped into a stirred polyvinyl alcohol (PVA, Chengdu Kelong Chemical Co., Ltd., China) to form the paraffin spheres. These spheres were then accumulated in a cylindrical mold and immersed in xylene (Chengdu Kelong Chemical Co., Ltd., China) and given pressure for a while to control the inter-sphere connectivity. The HA suspension was infused and gelled with water in air. Paraffin was wiped off by immersing in n-hexane (Chengdu Kelong Chemical Co., Ltd., China). The initial scaffolds were sintered into ceramic objects (Figure 2). The two macro-pore sizes of HAPPs were labeled L-HAPPs and S-HAPPs, respectively.
All the scaffolds were sealed in vacuum packages and sterilized with ethylene oxide for subsequent implantation.

4.2. Animal Experiment

All processes involving animal experiments were conducted in accordance with the Chinese local ethics committee and the experimental animal management system. The experiment follows the internationally recognized 3R principle (reduction, replacement, and refinement) for the more ethical use of animals, and conforms to the relevant contents of the guiding opinions on being kind to experimental animals. Twelve adult healthy dogs (1 year old, 15–20 kg) were used in this study for three time points (implantation for 1 month, 3 months, 6 months; four dogs for each time point). Immediately before surgery, sodium pentobarbital was dissolved in physiological saline (W% = 3%), sterilized by filtration, and intraperitoneally injected as a l mL/kg basis.
Two samples of four kinds of scaffolds were implanted, respectively, into the dorsal muscle and abdominal cavity of each animal, as described in detail as follows.
Dorsal muscle implantation: Scaffolds were implanted into the dorsal muscle following a previous study [55]. The dog was placed in a prone position, and the dorsal thoracolumbar area was shaved and disinfected. An incision was made along the dorsal midline from T8 to L5 to expose the para-spinal muscles. Eight intramuscular pockets (four/side) were created in the dorsal muscles parallel to the spine by blunt dissection, and a scaffold was placed in each pocket. Finally, the wound was closed in layers (Figure 3b).
Abdominal cavity implantation: The dog was turned to a supine position. The lower abdominal region was shaved and disinfected. Incisions were made parallel to (and ~2 cm from) the groins to expose the abdominal muscles, which were then blunt-separated to expose the peritoneum. The peritoneum was opened. Eight scaffolds (four/side) were placed into the abdominal cavity and sutured to the parietal peritoneum to ensure a direct contact with the parietal peritoneum and omentum. Finally, the wound was closed in layers (Figure 3c).
After the operations, all animals were given daily intramuscular injections of penicillin (80 × 104 IU) for 3 days. The animals were sacrificed by injecting an overdose of sodium pentobarbital after 1 month, 3 months and 6 months. The grafts were harvested and the tissue around the grafts was carefully removed. The general morphology of the specimens were observed.

4.3. Sample Preparation

Samples were fixed in 4% paraformaldehyde for 2 weeks at room temperature, rinsed with deionized water, dehydrated in an ethanol series and embedded with polymethyl methacrylate (PMMA) as per the method of our previous publications [35]. Then the samples were sectioned using a slicer (Microm HM360E, Thermo Fisher Scientific, Waltham, MA, USA) to yield 8 μm thick slices.

4.4. Histology, Histomorphometry, and Compression Tests

The slices were treated with H&E staining. The new bone formation, vascularization, inflammatory response and material–tissue interactions in the grafts were observed with optical microscopy (Nikon TS-100, Yokohama, Kanagawa, Japan). The percentage of the new bone area and the density of blood vessels were utilized to reveal the effects of the size and shape of the macro-pore structure on the osteoinduction and vascularization of the scaffolds. Vascular density: Each slice was randomly selected from 400 magnified fields of view, and the number of small vessels with a diameter less than 20 μm in each field of view was counted and averaged to represent the number of blood vessels per 0.2 mm2 unit area.
Compression tests (INSTRON 5567 Q4052, Canton, MA, USA) were performed to investigate the mechanical improvement of the scaffolds after the implantation. Referencing the ASTM D1621–2016 and ISO 844 standard test methods [56], the mechanical properties of the scaffolds were evaluated under the compression test at a strain rate of 0.5 mm/min using a 1 N load measurement and 0.005 mm deformation measurement accuracy. Four specimens were used for each scaffold.

4.5. Mathematical Models of the Macroporous Structure of 3D Scaffolds and Microfluidic Pathway

4.5.1. Geometry Reconstruction

A μ-CT scanner (μCT-50, SCANCO Medical AG, Bassersdorf, Switzerland) was used to acquire images of both scaffolds with a voxel resolution of 20 μm. The threshold gray value was determined by matching the material volume obtained from the Archimedes measurements to that of the solid 3D volume calculated from the images in the evaluation (SCANCO Medical AG, Bassersdorf, Switzerland).

4.5.2. CFD Model Creation

Our preliminary study found that two types of scaffolds (the HASAs and HAPPs) had complementary macro-porous structures [17]. The reconstruction images of the HASAs showed that the spheres accumulating type tends to the hexagonal close-accumulation (Figure 1a). Therefore, this study selected the hexagonal close-accumulation mode for the HASAs and its three-dimensional inversion mode for the HAPPs (Figure 1e), to construct two types of calculation models using CAD software SolideWorks 2014 (Dassault Systems S.A, Waltham, MA, USA). According to the macro-pore size and macro-pore shape of the actual scaffolds, the macro-pore sizes or sphere sizes of each model were set to 0.4 mm and 1.5 mm, respectively, and the sizes of the connective pores were set to 0.16 mm and 0.6 mm (the values were set according to the actual ratio of the connective pore size/macro-pore size = 0.4), respectively. Four models simulating the actual scaffold were used for the subsequent finite element analysis.
The fluid was defined as blood (viscosity of 0.033 Pa∙s, density of 1050 kg/m3, and the velocity-inlet of 0.01 m/s) according to previous studies [25,46], which was steady, 3D, laminar, and had a fully developed flow all over the scaffold. The Navier-Stokes equations governing fluid flow were solved by finite element software, Fluent 6.3 (ANSYS, Canonsburg, PA, USA). A no-slip boundary condition was applied on the flow channel wall. The flow speed near the walls was very low. The pressure was free at the flow channel outlets. To ensure the accuracy of flow simulation, each model was processed, in duplicate, four times.

4.5.3. Mesh and Numerical Simulation

With an irregular geometry no analytical solution is known, so shear stress must be estimated from the grids adjacent to the wall with the finite element method. Meshing was implemented in the geometry of the pore structure of the scaffold, and the elements of the model were obtained. The flow channel space was broken into a large number of elements, each comprised a set of cubic grids.
The model simulated a Newtonian fluid, and the Navier-Stokes equations describing the flow inside the pore structure of the scaffolds were solved. After the solution was completed, the postprocessing was undergone for quantifying the wall shear rates, as well as the velocity profiles, along the flow channels.
In order to reflect the variation of the flow fields in the different macro-pore structures, the center cross section of the area where the fluid flowed through the macro-pore and the interconnected pore in CFD models was selected as the flow field analysis contrast section.

4.6. Statistical Analysis

The results were expressed as mean ± standard deviation. The differences between the two groups were analyzed with a t-test. The differences between the three groups were analyzed with an analysis of variance (ANOVA; SPSS 15.0, IBM, Armonk, NY, USA) and a Tukey multiple comparison test. A p-value < 0.05 was considered statistically significant.

5. Conclusions

In this study, the correlation between the size and shape of macro-pores and the osteoinduction of HA porous scaffolds was systematically studied. An animal model of ectopic implantation was used to explore the relationship between the macroscopic pore structure characteristics of the scaffolds and the microfluidic dynamic environment, which was further confirmed by a finite element simulation. The results showed that the macroscopic pore structure is a key factor determining the microfluidic dynamic environment of the scaffold, which in turn affects the bone formation and distribution of the scaffold. The dynamic environment of a microfluid acts as an important intermediate pathway, connecting the macrostructural properties and regulating the osteoinduction of these materials. The results of this study have provided experiment and theoretical bases toward quantitative analysis research of the correlation between the microfluid dynamic environment and the ectopic osteogenesis of calcium phosphate ceramics, and the optimal fabrication of calcium phosphate ceramics with controlled osteoinduction. It should be noted that the scaffold pore model used for the finite element analysis in this study is a scaffold model with ideal structural features, rather than a 3D reconstruction model of the actual scaffold. This caused some discrepancy between the computational simulation results and the actual experimental results. However, this difference does not affect the final conclusion of this study.

Author Contributions

Conceptualization, W.Z. and J.W. (Jie Weng); methodology, F.S., T.Z. and X.F.; validation, F.S., X.F. and W.Z.; formal analysis, F.S., T.Z. and W.Z.; investigation, J.W. (Jianxin Wang) and S.Q.; resources, W.Z.; data curation, X.F. and X.H.; writing—original draft preparation, F.S., X.H. and W.Z.; writing—review and editing, K.D. and W.Z.; visualization, W.Z.; supervision, W.Z. and J.W. (Jie Weng); project administration, W.Z. and J.W. (Jie Weng); funding acquisition, F.S., W.Z. and J.W. (Jie Weng). All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key Research and Development Program of China (2016YFC1102000), the Science and Technology Innovation Seedling Project of Sichuan Province, China (2021057), and the Fundamental Research Funds for the Central Universities (2682021ZTPY003), the Sichuan Science and Technology Program (21MZGC0218), Degree and graduate Education and teaching Reform Project of SWJTU (YJG5-2022-Y018).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Xingdong, Z.; Pin, Z.; Jianguo, Z.; Weiqun, C.; Wu, C. A Study of Hydroxyapatite Ceramics and its Osteogenesis. In Bioceramics and the Human Body; Ravaglioli, A., Krajewski, A., Eds.; Springer: Dordrecht, The Netherlands, 1992; pp. 408–416. [Google Scholar] [CrossRef]
  2. Ripamonti, U. Osteoinduction in porous hydroxyapatite implanted in heterotopic sites of different animal models. Biomaterials 1996, 17, 31–35. [Google Scholar] [CrossRef]
  3. Yuan, H.; Li, Y.; de Bruijn, J.D.; de Groot, K.; Zhang, X. Tissue responses of calcium phosphate cement: A study in dogs. Biomaterials 2000, 21, 1283–1290. [Google Scholar] [CrossRef]
  4. Habibovic, P.; Gbureck, U.; Doillon, C.J.; Bassett, D.C.; van Blitterswijk, C.A.; Barralet, J.E. Osteoconduction and osteoinduction of low-temperature 3D printed bioceramic implants. Biomaterials 2008, 29, 944–953. [Google Scholar] [CrossRef] [PubMed]
  5. Tang, X.; Mao, L.; Liu, J.; Yang, Z.; Zhang, W.; Shu, M.; Hu, N.; Jiang, L.; Fang, B. Fabrication, characterization and cellular biocompatibility of porous biphasic calcium phosphate bioceramic scaffolds with different pore sizes. Ceram. Int. 2016, 42, 15311–15318. [Google Scholar] [CrossRef]
  6. Guo, X.; Li, M.; Qi, W.; Bai, H.; Nie, Z.; Hu, Z.; Xiao, Y.; de Bruijn, J.D.; Bao, C.; Yuan, H. Serial cellular events in bone formation initiated by calcium phosphate ceramics. Acta Biomater. 2021, 134, 730–743. [Google Scholar] [CrossRef]
  7. García-Gareta, E.; Coathup, M.J.; Blunn, G.W. Osteoinduction of bone grafting materials for bone repair and regeneration. Bone 2015, 81, 112–121. [Google Scholar] [CrossRef]
  8. Mehdizadeh, H.; Sumo, S.; Bayrak, E.S.; Brey, E.M.; Cinar, A. Three-dimensional modeling of angiogenesis in porous biomaterial scaffolds. Biomaterials 2013, 34, 2875–2887. [Google Scholar] [CrossRef]
  9. Perez, R.A.; Mestres, G. Role of pore size and morphology in musculo-skeletal tissue regeneration. Mater. Sci. Eng. C 2016, 61, 922–939. [Google Scholar] [CrossRef]
  10. Taniguchi, N.; Fujibayashi, S.; Takemoto, M.; Sasaki, K.; Otsuki, B.; Nakamura, T.; Matsushita, T.; Kokubo, T.; Matsuda, S. Effect of pore size on bone ingrowth into porous titanium implants fabricated by additive manufacturing: An in vivo experiment. Mater. Sci. Eng. C 2016, 59, 690–701. [Google Scholar] [CrossRef] [Green Version]
  11. Wang, X.; Xu, S.; Zhou, S.; Xu, W.; Leary, M.; Choong, P.; Qian, M.; Brandt, M.; Xie, Y.M. Topological design and additive manufacturing of porous metals for bone scaffolds and orthopaedic implants: A review. Biomaterials 2016, 83, 127–141. [Google Scholar] [CrossRef]
  12. Zadpoor, A.A. Bone tissue regeneration: The role of scaffold geometry. Biomater. Sci. 2015, 3, 231–245. [Google Scholar] [CrossRef] [PubMed]
  13. Xiao, X.; Wang, W.; Liu, D.; Zhang, H.; Gao, P.; Geng, L.; Yuan, Y.; Lu, J.; Wang, Z. The promotion of angiogenesis induced by three-dimensional porous beta-tricalcium phosphate scaffold with different interconnection sizes via activation of PI3K/Akt pathways. Sci. Rep. 2015, 5, 9409. [Google Scholar] [CrossRef] [PubMed]
  14. Fu, J.; Wiraja, C.; Muhammad, H.B.; Xu, C.; Wang, D.-A. Improvement of endothelial progenitor outgrowth cell (EPOC)-mediated vascularization in gelatin-based hydrogels through pore size manipulation. Acta Biomater. 2017, 58, 225–237. [Google Scholar] [CrossRef] [PubMed]
  15. Harada, S.-i.; Rodan, G.A. Control of osteoblast function and regulation of bone mass. Nature 2003, 423, 349–355. [Google Scholar] [CrossRef]
  16. Bodde, E.W.H.; Cammaert, C.T.R.; Wolke, J.G.C.; Spauwen, P.H.M.; Jansen, J.A. Investigation as to the osteoinductivity of macroporous calcium phosphate cement in goats. J. Biomed. Mater. Res. Part B Appl. Biomater. 2007, 83B, 161–168. [Google Scholar] [CrossRef]
  17. Wang, H.; Zhi, W.; Lu, X.; Li, X.; Duan, K.; Duan, R.; Mu, Y.; Weng, J. Comparative studies on ectopic bone formation in porous hydroxyapatite scaffolds with complementary pore structures. Acta Biomater 2013, 9, 8413–8421. [Google Scholar] [CrossRef]
  18. Barradas, A.M.; Yuan, H.; van Blitterswijk, C.A.; Habibovic, P. Osteoinductive biomaterials: Current knowledge of properties, experimental models and biological mechanisms. Eur. Cell Mater. 2011, 21, 407–429; discussion 429. [Google Scholar] [CrossRef]
  19. Karageorgiou, V.; Kaplan, D. Porosity of 3D biomaterial scaffolds and osteogenesis. Biomaterials 2005, 26, 5474–5491. [Google Scholar] [CrossRef]
  20. Chang, B.; Song, W.; Han, T.; Yan, J.; Li, F.; Zhao, L.; Kou, H.; Zhang, Y. Influence of pore size of porous titanium fabricated by vacuum diffusion bonding of titanium meshes on cell penetration and bone ingrowth. Acta Biomater. 2016, 33, 311–321. [Google Scholar] [CrossRef]
  21. Choi, D.J.; Park, S.J.; Gu, B.K.; Kim, Y.-J.; Chung, S.; Kim, C.-H. Effect of the pore size in a 3D bioprinted gelatin scaffold on fibroblast proliferation. J. Ind. Eng. Chem. 2018, 67, 388–395. [Google Scholar] [CrossRef]
  22. Yeatts, A.B.; Fisher, J.P. Bone tissue engineering bioreactors: Dynamic culture and the influence of shear stress. Bone 2011, 48, 171–181. [Google Scholar] [CrossRef] [PubMed]
  23. Zhao, F.; Lacroix, D.; Ito, K.; van Rietbergen, B.; Hofmann, S. Changes in scaffold porosity during bone tissue engineering in perfusion bioreactors considerably affect cellular mechanical stimulation for mineralization. Bone Rep. 2020, 12, 100265. [Google Scholar] [CrossRef] [PubMed]
  24. Montorsi, M.; Genchi, G.G.; De Pasquale, D.; De Simoni, G.; Sinibaldi, E.; Ciofani, G. Design, fabrication, and characterization of a multimodal reconfigurable bioreactor for bone tissue engineering. Biotechnol. Bioeng. 2022, 119, 1965–1979. [Google Scholar] [CrossRef]
  25. Shanglong, X.; Pingan, D.; Youzhuan, X.; Yang, Y. Cell distribution in a scaffold with random architectures under the influence of fluid dynamics. J. Biomater. Appl. 2008, 23, 229–245. [Google Scholar] [CrossRef] [PubMed]
  26. Porter, B.; Zauel, R.; Stockman, H.; Guldberg, R.; Fyhrie, D. 3-D computational modeling of media flow through scaffolds in a perfusion bioreactor. J. Biomech. 2005, 38, 543–549. [Google Scholar] [CrossRef] [PubMed]
  27. Zhang, C.; Huang, P.; Weng, J.; Zhi, W.; Hu, Y.; Feng, H.; Yao, Y.; Li, S.; Xia, T. Histomorphological researches on large porous hydroxyapatite cylinder tubes with polylactic acid surface coating in different nonskeletal sites in vivo. J. Biomed. Mater. Res. A 2012, 100, 1203–1208. [Google Scholar] [CrossRef]
  28. Zhang, Y.; Sun, N.; Zhu, M.; Qiu, Q.; Zhao, P.; Zheng, C.; Bai, Q.; Zeng, Q.; Lu, T. The contribution of pore size and porosity of 3D printed porous titanium scaffolds to osteogenesis. Biomater. Adv. 2022, 133, 112651. [Google Scholar] [CrossRef]
  29. Peng, Q.; Jiang, F.; Huang, P.; Zhou, S.; Weng, J.; Bao, C.; Zhang, C.; Yu, H. A novel porous bioceramics scaffold by accumulating hydroxyapatite spherules for large bone tissue engineering in vivo. I. Preparation and characterization of scaffold. J. Biomed. Mater. Res. Part A 2010, 93A, 920–929. [Google Scholar] [CrossRef]
  30. Zhao, J.; Li, J.-Y.; Zhi, W.; Lu, X.; Jia, Z.-B.; Weng, J. Preparation and Optimization of Porous HA Ceramic Scaffolds by Wax Spheres Leaching Method. J. Inorg. Mater. 2013, 28, 74–78. [Google Scholar] [CrossRef]
  31. Zhao, J.; Guo, L.Y.; Yang, X.B.; Weng, J. Preparation of bioactive porous HA/PCL composite scaffolds. Appl. Surf. Sci. 2008, 255, 2942–2946. [Google Scholar] [CrossRef]
  32. Mastrogiacomo, M.; Scaglione, S.; Martinetti, R.; Dolcini, L.; Beltrame, F.; Cancedda, R.; Quarto, R. Role of scaffold internal structure on in vivo bone formation in macroporous calcium phosphate bioceramics. Biomaterials 2006, 27, 3230–3237. [Google Scholar] [CrossRef] [PubMed]
  33. Bai, F.; Wang, Z.; Lu, J.; Liu, J.; Chen, G.; Lv, R.; Wang, J.; Lin, K.; Zhang, J.; Huang, X. The Correlation Between the Internal Structure and Vascularization of Controllable Porous Bioceramic Materials In Vivo: A Quantitative Study. Tissue Eng. Part A 2010, 16, 3791–3803. [Google Scholar] [CrossRef] [PubMed]
  34. Li, J.; Xu, T.; Hou, W.; Liu, F.; Qing, W.; Huang, L.; Ma, G.; Mu, Y.; Weng, J. The response of host blood vessels to graded distribution of macro-pores size in the process of ectopic osteogenesis. Mater. Sci. Eng. C 2020, 109, 110641. [Google Scholar] [CrossRef]
  35. Li, J.; Huang, H.; Xu, T.; Li, J.; Guo, T.; Lu, X.; Ren, J.; Ren, X.; Mu, Y.; Weng, J. Effect of the interconnecting window diameter of hydroxyapatite scaffolds on vascularization and osteoinduction. Ceram. Int. 2022, 48, 25070–25078. [Google Scholar] [CrossRef]
  36. Habibovic, P.; Sees, T.M.; Doel, M.A.V.D.; Blitterswijk, C.A.V.; Groot, K.D. Osteoinduction by biomaterials—Physicochemical and structural influences. J. Biomed. Mater. Res. Part A 2006, 77, 747–762. [Google Scholar] [CrossRef] [PubMed]
  37. Ripamonti, U.; Crooks, J.; Kirkbride, A.N. Sintered porous hydroxyapatites with intrinsic osteoinductive activity: Geometric induction of bone formation. S. Afr. J. Sci. 1999, 95, 335–343. [Google Scholar] [CrossRef]
  38. Yuan, H.; Fernandes, H.; Habibovic, P.; de Boer, J.; Barradas, A.M.C.; de Ruiter, A.; Walsh, W.R.; van Blitterswijk, C.A.; de Bruijn, J.D. Osteoinductive ceramics as a synthetic alternative to autologous bone grafting. Proc. Natl. Acad. Sci. USA 2010, 107, 13614–13619. [Google Scholar] [CrossRef]
  39. Yuan, H.P.; Yang, Z.J.; Li, Y.B.; Zhang, X.D.; De Bruijn, J.D.; De Groot, K. Osteoinduction by calcium phosphate biomaterials. J. Mater. Sci.-Mater. Med. 1998, 9, 723–726. [Google Scholar] [CrossRef]
  40. Yuan, H.; van Blitterswijk, C.A.; de Groot, K.; de Bruijn, J.D. A comparison of bone formation in biphasic calcium phosphate (BCP) and hydroxyapatite (HA) implanted in muscle and bone of dogs at different time periods. J. Biomed. Mater. Res. Part A 2006, 78A, 139–147. [Google Scholar] [CrossRef]
  41. Zhang, Q.; Ma, L.; Ji, X.; He, Y.; Cui, Y.; Liu, X.; Xuan, C.; Wang, Z.; Yang, W.; Chai, M.; et al. High-Strength Hydroxyapatite Scaffolds with Minimal Surface Macrostructures for Load-Bearing Bone Regeneration. Adv. Funct. Mater. 2022, 32, 2204182. [Google Scholar] [CrossRef]
  42. Ripamonti, U.; Van Den Heever, B.; Crooks, J.; Tucker, M.M.; Sampath, T.K.; Rueger, D.C.; Reddi, A.H. Long-term evaluation of bone formation by osteogenic protein 1 in the baboon and relative efficacy of bone-derived bone morphogenetic proteins delivered by irradiated xenogeneic collagenous matrices. J. Bone Miner. Res. 2000, 15, 1798–1809. [Google Scholar] [CrossRef] [PubMed]
  43. Yang, Z.; Wang, C.; Gao, H.; Jia, L.; Zeng, H.; Zheng, L.; Wang, C.; Zhang, H.; Wang, L.; Song, J.; et al. Biomechanical Effects of 3D-Printed Bioceramic Scaffolds With Porous Gradient Structures on the Regeneration of Alveolar Bone Defect: A Comprehensive Study. Front. Bioeng Biotechnol. 2022, 10, 882631. [Google Scholar] [CrossRef] [PubMed]
  44. Suvarnapathaki, S.; Wu, X.; Lantigua, D.; Nguyen, M.A.; Camci-Unal, G. Hydroxyapatite-Incorporated Composite Gels Improve Mechanical Properties and Bioactivity of Bone Scaffolds. Macromol. Biosci. 2020, 20, 2000176. [Google Scholar] [CrossRef]
  45. Wu, X.; Zhang, T.; Hoff, B.; Suvarnapathaki, S.; Lantigua, D.; McCarthy, C.; Wu, B.; Camci-Unal, G. Mineralized Hydrogels Induce Bone Regeneration in Critical Size Cranial Defects. Adv. Healthc. Mater. 2021, 10, 2001101. [Google Scholar] [CrossRef]
  46. Jaasma, M.J.; Plunkett, N.A.; O’Brien, F.J. Design and validation of a dynamic flow perfusion bioreactor for use with compliant tissue engineering scaffolds. J. Biotechnol. 2008, 133, 490–496. [Google Scholar] [CrossRef] [PubMed]
  47. Bonvin, C.; Overney, J.; Shieh, A.C.; Dixon, J.B.; Swartz, M.A. A multichamber fluidic device for 3D cultures under interstitial flow with live imaging: Development, characterization, and applications. Biotechnol. Bioeng. 2010, 105, 982–991. [Google Scholar] [CrossRef]
  48. Nihouannen, D.L.; Daculsi, G.; Saffarzadeh, A.; Gauthier, O.; Delplace, S.; Pilet, P.; Layrolle, P. Ectopic bone formation by microporous calcium phosphate ceramic particles in sheep muscles. Bone 2005, 36, 1086. [Google Scholar] [CrossRef]
  49. Le Huec, J.C.; Clément, D.; Brouillaud, B.; Barthe, N.; Dupuy, B.; Foliguet, B.; Basse-Cathalinat, B. Evolution of the local calcium content around irradiated β-tricalcium phosphate ceramic implants: In vivo study in the rabbit. Biomaterials 1998, 19, 733–738. [Google Scholar] [CrossRef]
  50. Roux, S.; Orcel, P. Bone loss: Factors that regulate osteoclast differentiation—An update. Arthritis Res. Ther. 2000, 2, 451–456. [Google Scholar] [CrossRef] [Green Version]
  51. Yang, Z.J.; Yuan, H.P.; Zou, P.; Tong, W.D.; Qu, S.X.; Zhang, X.D. Osteogenic responses to extraskeletally implanted synthetic porous calcium phosphate ceramics: An early stage histomorphological study in dogs. J. Mater. Sci.-Mater. Med. 1997, 8, 697–701. [Google Scholar] [CrossRef]
  52. Habibovic, P.; Yuan, H.P.; van der Valk, C.M.; Meijer, G.; van Blitterswijk, C.A.; de Groot, K. 3D microenvironment as essential element for osteoinduction by biomaterials. Biomaterials 2005, 26, 3565–3575. [Google Scholar] [CrossRef] [PubMed]
  53. Peng, Q.; Qu, S.X.; Lu, X.Y.; Weng, J. Fabricating Hydroxyapatite Spherulites with Characteristic Microporous Structure by Chitin Sol Emulsified in Oil and Gelatinized In Situ. Key Eng. Mater. 2005, 284–286, 419–422. [Google Scholar] [CrossRef]
  54. Zhao, J.; Lu, X.; Duan, K.; Guo, L.Y.; Zhou, S.B.; Weng, J. Improving mechanical and biological properties of macroporous HA scaffolds through composite coatings. Colloids Surf. B Biointerfaces 2009, 74, 159–166. [Google Scholar] [CrossRef]
  55. Zhi, W.; Zhang, C.; Duan, K.; Li, X.; Qu, S.; Wang, J.; Zhu, Z.; Huang, P.; Xia, T.; Liao, G.; et al. A novel porous bioceramics scaffold by accumulating hydroxyapatite spherulites for large bone tissue engineering in vivo. II. Construct large volume of bone grafts. J. Biomed. Mater. Res. Part A 2014, 102, 2491–2501. [Google Scholar] [CrossRef] [PubMed]
  56. Karimipour-Fard, P.; Jeffrey, M.P.; JonesTaggart, H.; Pop-Iliev, R.; Rizvi, G. Development, processing and characterization of Polycaprolactone/Nano-Hydroxyapatite/Chitin-Nano-Whisker nanocomposite filaments for additive manufacturing of bone tissue scaffolds. J. Mech. Behav. Biomed. Mater. 2021, 120, 104583. [Google Scholar] [CrossRef]
Ijms 23 11459 g001 550
Figure 1. 3D reconstruction using Micro-CT and showing macropore of HASAs and HAPPs.The stellate-shape macropore of HASAs (a). The pore size was decided by the size of HA spheres ((b): S-HASAs; (c): L-HASAs) and the arrangement of the spheres (d). The spherical-shape macropore of HAPPs (e). SEM showed the macropores and interconnective pores of L-HAPPs (f) and S-HAPPs (g).
Figure 1. 3D reconstruction using Micro-CT and showing macropore of HASAs and HAPPs.The stellate-shape macropore of HASAs (a). The pore size was decided by the size of HA spheres ((b): S-HASAs; (c): L-HASAs) and the arrangement of the spheres (d). The spherical-shape macropore of HAPPs (e). SEM showed the macropores and interconnective pores of L-HAPPs (f) and S-HAPPs (g).
Ijms 23 11459 g001
Ijms 23 11459 g002 550
Figure 2. Relationship between ectopic bone formations of HASAs and the macropore sizes implanted in dorsal muscle/abdominal cavity by H&E staining. (a1–a3) L-HASAs implanted 1, 3 and 6 months, (b1–b3) S-HASAs implanted 1, 3 and 6 months, (c1–c3) L-HASAs implanted 1, 3 and 6 months, (d1–d3) S-HASAs implanted 1, 3 and 6 months. With the implantation time extended, new bone tissues were found in all scaffolds, which was better in the L-HASAs than those in the S-HASAs. White triangles show HA materials, yellow stars show new bone. The bar = 2 mm.
Figure 2. Relationship between ectopic bone formations of HASAs and the macropore sizes implanted in dorsal muscle/abdominal cavity by H&E staining. (a1–a3) L-HASAs implanted 1, 3 and 6 months, (b1–b3) S-HASAs implanted 1, 3 and 6 months, (c1–c3) L-HASAs implanted 1, 3 and 6 months, (d1–d3) S-HASAs implanted 1, 3 and 6 months. With the implantation time extended, new bone tissues were found in all scaffolds, which was better in the L-HASAs than those in the S-HASAs. White triangles show HA materials, yellow stars show new bone. The bar = 2 mm.
Ijms 23 11459 g002
Ijms 23 11459 g003 550
Figure 3. Relationship between the ectopic bone formations of HAPPs and the macropore sizes implanted in dorsal muscle/abdominal cavity by H&E staining. (a1–a3) L-HAPPs implanted 1, 3 and 6 months, (b1–b3) S-HASAs implanted 1, 3 and 6 months, (c1–c3) L-HAPPs implanted 1, 3 and 6 months, (d1–d3) S-HASAs implanted 1, 3 and 6 months. With the implantation time extended, new bone tissues were found in all scaffolds, which was better in the L-HAPPs than those in the S-HAPPs. Black triangles show HA materials, yellow stars show new bone. The bar = 2 mm.
Figure 3. Relationship between the ectopic bone formations of HAPPs and the macropore sizes implanted in dorsal muscle/abdominal cavity by H&E staining. (a1–a3) L-HAPPs implanted 1, 3 and 6 months, (b1–b3) S-HASAs implanted 1, 3 and 6 months, (c1–c3) L-HAPPs implanted 1, 3 and 6 months, (d1–d3) S-HASAs implanted 1, 3 and 6 months. With the implantation time extended, new bone tissues were found in all scaffolds, which was better in the L-HAPPs than those in the S-HAPPs. Black triangles show HA materials, yellow stars show new bone. The bar = 2 mm.
Ijms 23 11459 g003
Ijms 23 11459 g004 550
Figure 4. New bone formation and vessel density in HA porous scaffolds with different macropore sizes implanted in dorsal muscle and abdominal cavity using histomorphometry analysis. (a): new bone formation area % of HA scaffolds implanted in dorsal muscle; (b): new bone formation area % of HA scaffolds implanted in abdominal cavity; (c): vessel density of HA scaffolds implanted in dorsal muscle; (d): vessel density of HA scaffolds implanted in abdominal cavity (* p < 0.05).
Figure 4. New bone formation and vessel density in HA porous scaffolds with different macropore sizes implanted in dorsal muscle and abdominal cavity using histomorphometry analysis. (a): new bone formation area % of HA scaffolds implanted in dorsal muscle; (b): new bone formation area % of HA scaffolds implanted in abdominal cavity; (c): vessel density of HA scaffolds implanted in dorsal muscle; (d): vessel density of HA scaffolds implanted in abdominal cavity (* p < 0.05).
Ijms 23 11459 g004
Ijms 23 11459 g005 550
Figure 5. The compressive strength in HA porous scaffolds with different macropore sizes implanted in dorsal muscle (a) and abdominal cavity (b) (* p < 0.05).
Figure 5. The compressive strength in HA porous scaffolds with different macropore sizes implanted in dorsal muscle (a) and abdominal cavity (b) (* p < 0.05).
Ijms 23 11459 g005
Ijms 23 11459 g006 550
Figure 6. CFD fluid simulation analysis of HASAs and HAPPS scaffolds reveal the relationship between distribution of flow field and macropore size. The distribution of flow field on the horizontal plane in: L-HASAs (a1), S-HASAs (b1), L-HAPPs (c1), S-HAPPs (d1). The distribution of flow field on the sagittal plane in: L-HASAs (a2), S-HASAs (b2), L-HAPPs (c2), S-HAPPs (d2). Color bar shows flow velocity.
Figure 6. CFD fluid simulation analysis of HASAs and HAPPS scaffolds reveal the relationship between distribution of flow field and macropore size. The distribution of flow field on the horizontal plane in: L-HASAs (a1), S-HASAs (b1), L-HAPPs (c1), S-HAPPs (d1). The distribution of flow field on the sagittal plane in: L-HASAs (a2), S-HASAs (b2), L-HAPPs (c2), S-HAPPs (d2). Color bar shows flow velocity.
Ijms 23 11459 g006
Ijms 23 11459 g007 550
Figure 7. The osteoinductive mechanism of macropore structure of HA scaffolds.
Figure 7. The osteoinductive mechanism of macropore structure of HA scaffolds.
Ijms 23 11459 g007
Table
Table 1. Histological summary of ectopic bone formations of HA porous scaffolds implanted in dorsal muscle and abdominal cavity.
Table 1. Histological summary of ectopic bone formations of HA porous scaffolds implanted in dorsal muscle and abdominal cavity.
Implant TimeScaffoldDorsal MuscleAbdominal Cavity
New Bone LocationNew Bone Formation (%)Vessel Density
(mm−2)		New Bone LocationNew Bone Formation (%)Vessel Density
(mm−2)
1 ML-HASAs\\2.34 ± 0.55\\1.86 ± 0.31
S-HASAs\\1.60 ± 0.79\\1.43 ± 1.02
L-HAPPsCentral4.27 ± 1.852.67 ± 0.79Central3.06 ± 0.941.83 ± 0.23
S-HAPPs\\1.64 ± 0.21\\1.39 ± 0.83
3 ML-HASAsOuter10.36 ± 1.483.18 ± 0.37Outer5.47 ± 1.932.59 ± 0.41
S-HASAsOuter6.18 ± 2.53.31 ± 0.72Outer1.74 ± 0.472.13 ± 0.47
L-HAPPsAll16.57 ± 2.794.12 ± 0.83All6.38 ± 2.053.15 ± 0.37
S-HAPPsCentral2.14 ± 1.332.58 ± 0.62Central1.28 ± 0.341.75 ± 0.15
6 ML-HASAsAll22.27 ± 2.894.02 ± 0.46Outer10.28 ± 2.063.17 ± 0.51
S-HASAsOuter14.0 ± 3.674.78 ± 0.45Outer5.65 ± 1.172.87 ± 0.39
L-HAPPsAll28.31 ± 3.066.44 ± 0.75All15.93 ± 3.874.38 ± 0.53
S-HAPPsAll7.61 ± 2.483.11 ± 1.03Central3.94 ± 1.331.90 ± 0.63
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Shi, F.; Fang, X.; Zhou, T.; Huang, X.; Duan, K.; Wang, J.; Qu, S.; Zhi, W.; Weng, J. Macropore Regulation of Hydroxyapatite Osteoinduction via Microfluidic Pathway. Int. J. Mol. Sci. 2022, 23, 11459. https://doi.org/10.3390/ijms231911459

AMA Style

Shi F, Fang X, Zhou T, Huang X, Duan K, Wang J, Qu S, Zhi W, Weng J. Macropore Regulation of Hydroxyapatite Osteoinduction via Microfluidic Pathway. International Journal of Molecular Sciences. 2022; 23(19):11459. https://doi.org/10.3390/ijms231911459

Chicago/Turabian Style

Shi, Feng, Xin Fang, Teng Zhou, Xu Huang, Ke Duan, Jianxin Wang, Shuxin Qu, Wei Zhi, and Jie Weng. 2022. "Macropore Regulation of Hydroxyapatite Osteoinduction via Microfluidic Pathway" International Journal of Molecular Sciences 23, no. 19: 11459. https://doi.org/10.3390/ijms231911459

APA Style

Shi, F., Fang, X., Zhou, T., Huang, X., Duan, K., Wang, J., Qu, S., Zhi, W., & Weng, J. (2022). Macropore Regulation of Hydroxyapatite Osteoinduction via Microfluidic Pathway. International Journal of Molecular Sciences, 23(19), 11459. https://doi.org/10.3390/ijms231911459

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop