Microporosity Clustering Assessment in Calcium Phosphate Bioceramic Particles
Raúl Vallejos Baier1 
Isabel Benjumeda Wijnhoven2 Víctor Irribarra del Valle1 Carola Millán Giovanetti2 Juan F. Vivanco1*- 1Faculty of Engineering and Sciences, Adolfo Ibáñez University, Viña del Mar, Chile
- 2Faculty of Liberal Arts, Adolfo Ibáñez University, Santiago, Chile
There has been an increase in the application of different biomaterials to repair hard tissues. Within these biomaterials, calcium phosphate (CaP) bioceramics are suitable candidates, since they can be biocompatible, biodegradable, osteoinductive, and osteoconductive. Moreover, during sintering, bioceramic materials are prone to form micropores and undergo changes in their surface topographical features, which influence cellular physiology and bone ingrowth. In this study, five geometrical properties from the surface of CaP bioceramic particles and their micropores were analyzed by data mining techniques, driven by the research question: what are the geometrical properties of individual micropores in a CaP bioceramic, and how do they relate to each other? The analysis not only shows that it is feasible to determine the existence of micropore clusters, but also to quantify their geometrical properties. As a result, these CaP bioceramic particles present three groups of micropore clusters distinctive by their geometrical properties. Consequently, this new methodological clustering assessment can be applied to advance the knowledge about CaP bioceramics and their role in bone tissue engineering.
Introduction
It is estimated that ~1.5 million devices are implanted worldwide per year to heal musculoskeletal diseases (Holzapfel et al., 2013), these procedures present difficulties that have economic implications for health care providers. Likewise, given the complications in the access or usage of grafts, there has been an increase in the application of different biomaterials in order to repair hard tissues, such as bone and teeth (Giannoudis et al., 2005).
One type of biomaterials currently used for hard tissue regeneration are calcium phosphate (CaP) based bioceramics, which not only have a similar composition to the mineral phase of bone (Dorozhkin, 2010), but can be biocompatible, biodegradable, osteoinductive, and osteoconductive (Vivanco et al., 2012; Denry and Kuhn, 2016; Kim et al., 2017). In addition, CaP bioceramics promote rapid bone formation and may assure bone healing within a year (Habraken et al., 2016). Moreover, it has been shown that calcium and phosphate ions trigger an osteoinductive response during bone regeneration, being resorbed by cell-mediated processes, therefore controlling the potential toxicity of degradation products (Habraken et al., 2016).
Recent advances have started to focus on using CaP bioceramics for bone tissue regeneration, given that they can be inserted in the defect site with a minimally invasive surgery, preventing the risks of infections, surgical scars, and blood loss (Low et al., 2010; Uswatta et al., 2016). Among CaP bioceramics, hydroxyapatite (HA) and tricalcium phosphate (TCP) have been commonly used in clinical applications and in in vivo studies. For example, some studies have shown that 95% of these CaP bioceramics are resorbed in 26–28 weeks (Knaack et al., 1998; Wiltfang et al., 2002), with crystalline TCP having a higher degradation rate than crystalline HA (Vicente et al., 1996; El-Ghannam, 2004). This degradation profile is a desirable property, given that it allows the replacement of the bioceramic material with newly synthetized bone.
Previous evidence suggests that the structural and material properties of a CaP bioceramic strongly influence its capacity to induce new bone formation (Habibovic and de Groot, 2007). Additionally, it is not only the material per se or its fabrication parameters which influence biological processes (Habraken et al., 2016), but also its final microarchitecture, such as pore size and pore distribution (Gauthier et al., 1998; Mastrogiacomo et al., 2006; Novotna et al., 2019). Correspondingly, the porosity of a bioceramic is an important requirement for vascularization (Karageorgiou and Kaplan, 2005).
Although there is a general agreement that pore size is a key factor affecting cell ingrowth and bone formation in bioceramics (Hutmacher, 2000), there is no conclusive data regarding optimal pore size. For example, a minimal pore diameter of 100 μm has been proposed to influence cell ingrowth (Karageorgiou and Kaplan, 2005), while a diameter of 200 μm or more has been proposed to support new bone formation (Gauthier et al., 1998; Flautre et al., 2001; Galois and Mainard, 2004). Interestingly, some in vivo studies have suggested that there is no significant difference in bone regeneration for pore sizes in the range of 400–1,200 μm (Hollister et al., 2005; Schek et al., 2006). However, there is evidence for bone formation in structures with interconnected micro-pores of <10 μm in size, while also having a macro-porosity of more than 100 μm (Lan Levengood et al., 2010a,b).
It is well-known that microporosity affects the process of osteogenesis (Zhang et al., 2018; Rustom et al., 2019), for instance, the increased specific surface areas by microporosity in biomaterials can offer more protein adsorption sites. Additionally, the capillary force generated by this microporosity can improve the attachment and immigration of bone-related cells on the biomaterial surface, with these cells penetrating them even if these micropores are smaller (Polak et al., 2013; Zhang et al., 2018). The contribution of capillarity induced by microporosity on biphasic CaP biomaterials has also been tested in vivo, by experiments with pig mandibles, which showed that capillarity induced by microporosity improved the homogeneity of bone distribution in these biomaterials (Rustom et al., 2016).
Although the above-mentioned features of a biomaterial, such as the specific surface area and the capillary force generated by its microporosity, have been shown to play a role in osteogenesis, detailed analyses of each of these features might not be yet fully quantified. Likewise, while properties such as pore size, localization, or gradients have been studied in relation to cell growth and bone formation, more global properties related to their organization, such as pore clustering, have not been thoroughly studied.
The process of bioceramic sintering can be prone to form a complex microarchitecture with high surface roughness and microporosity (Wilson et al., 2004; Champion, 2013). For example, the creation and distribution of micropores during sintering can create clusters of them with amounts in the order of 103, which could be difficult to quantify with the traditional tools used in tissue engineering. Thus, a necessity appears to develop techniques to analyse these big data sets. These techniques can be used to relate these large micropore datasets and their influence on the mechanical and biological performance of bioceramics.
In this work, a novel method is used to asses pore clustering and to study the microarchitecture of CaP bioceramic particles as a model that presents both micro and macroporosity, driven by the research question: what are the geometrical properties of individual micropores in a CaP bioceramic, and how do they relate to each other? Material properties were studied by crystallographic phases and microstructural morphology. Hence, this study adds to the knowledge of CaP bioceramics and their use in regenerative medicine.
Materials and Methods
Bioceramic Material
In this work, the bioceramics used as a model were commercial Megagen Boneplus Eagle Eye Synthetic Bone Graft (http://www.imegagen.com).
These CaP bioceramics were used as received, which come in particles of a torus shape with a diameter of max 1 mm. According to the manufacturer, this material is made of a synthetic HA/β-TCP composite (60–40%), and each particle presents an interconnected micropore channel structure, while also having a central macro-pore (as shown in Figure 1).
Figure 1. CaP bioceramic particle microarchitecture. (A) SEM micrograph of the bioceramic surface. (B) Area of interest marked in (A), showing the microporosity of the bioceramic, the yellow square indicates a representative micropore. (C) Diagram of the five attributes measured for each micropore, based on Rueden et al. (2017).
Microstructural Morphology and Phase Identification
The microstructural morphology of the bioceramic surface was analyzed by scanning electron microscopy (SEM). Samples were mounted on aluminum stubs with double-sided carbon tape and sputter coated with gold for 30 s at 45 mA (Denton Vacuum Desk V). Subsequently, SEM images were obtained (Jeol JSM IT300LV).
The crystallographic phase was determined by means of X-ray diffraction (XRD) (Bruker D8 Advance), with Cu Kα1 radiation (λ = 1.5406 Å) operated at 40 kV and 30 mA. The XRD patterns were recorded in the 2θ range of 10° to 80°, with a step size of 0.001° and step duration of 34 s. The recorded XRD spectra were identified by matching the spectra with ones based on the structural data of similar apatite bioceramics available in the Crystallography Open Database (COD). The phase composition analysis was done using the Rietveld Method (Bish and Howard, 1988).
Assessment of Micropore Clusters of Cap Bioceramic Particles
A total of five (n = 5) bioceramic particles were analyzed by SEM imaging, obtaining more than 5 × 103 micropores to be analyzed. The images were processed using ImageJ software (NIH, USA, version 1.51j8), with the aim of isolating only the micropores from the surface of the bioceramic particles, to perform the measurements and the subsequent analysis of clusters.
In order to isolate the micropores of each bioceramic particle, the threshold tool from ImageJ was used, in which the SEM microphotographs were converted to their corresponding binary images, where black pixels represent empty spaces (pores), following a similar protocol used previously (Xu and Chan, 2008). The validation of this porosity measurement was done using an accepted methodology (Xu et al., 2007).
For each of the more than 5 × 103 isolated micropores (observations) the values of five geometrical attributes (variables) were determined: area, perimeter, circularity, roundness, and Feret's diameter (Walton, 1948). The area is expressed in μm2 and the perimeter in μm. The circularity and roundness are expressed in values between 0 (one line) and 1 (one circle and sphere, respectively). Feret's diameter represents the distance between two parallel lines that are tangential to the contour of the projection of a particle, in this case micropores, which allows having a measurement in μm of the diameter of the micropores that have an irregular shape (Merkus, 2009). Details of these parameters are shown in Figure 1.
Clustering Analysis was conducted by using the K-means clustering algorithm, which is a method that finds K vectors that represent a complete data set and aims to minimize the sum of the squared distances between all the points and the center of the cluster (Gonzalez and Tou, 1974; Pham et al., 2005). K-means was applied to a series of groups from 1 to 10, and inertia values were obtained (sum of the square distances of each cluster micropore to its centroid), which were used to determine the number of groups by means of the “Elbow” method. This method allows representing in a linear graph the inertia with respect to the number of clusters, indicating an inflection point and evaluating the maximal inflection. This point (elbow) indicates the optimal number of groups for the data set (Kaufman and Rousseeuw, 1990).
The validation of the number of clusters obtained was done by the Davies-Bouldin index (DBI), which is based on the fact that those algorithms that produce clusters with low intracluster distances (high intracluster similarity) and high intercluster distances (low intercluster similarity) will have a low DBI. Consequently, the clustering algorithm that generates a collection of clusters with the lowest value of this index was considered the best algorithm (Davies and Bouldin, 1979).
Finally, the analysis of clusters was deepened by relating the Feret's diameter of the micropores with their circularity (RStudio Inc., version 1.0.153), since it is known that the size and shape of the pores of a biomaterial is critical for cellular behavior (Akay et al., 2004; Zadpoor, 2015).
Theory/Calculation
Determination of Davies-Bouldin Index
The validation of the number of clusters obtained was done by means of the DBI (Davies and Bouldin, 1979), using the following equation:
Where n is the number of clusters, cx is the centroid of cluster x, σx is the mean distance of all observations in cluster x to the centroid cx, and d(ci,cj) is the distance between the centroids ci and cj.
Results
Analysis of Cap Bioceramic Particle Microarchitecture
Since the architecture and microporosity of CaP biomaterials can influence its mechanical and biological performance (Rustom et al., 2019), an analysis of the bioceramic microarchitecture was carried out from SEM images. These data were processed using ImageJ, with the aim of isolating only the micropores from the surface of the bioceramic particles, to perform the measurements and the subsequent cluster analysis.
The micropores were individually analyzed, obtaining 5,338 observations from 5 bioceramic particles, in which the area, perimeter, circularity, roundness, and Feret's diameter were measured (Figure 1).
Microporosity Cluster Determination of Cap Bioceramic Particles
Clustering is a data exploration technique that allows objects with similar characteristics to be grouped, in order to facilitate their subsequent processing. In this particular case, the K-means method was used.
K-means was applied to a number of clusters from 1 to 10, and the inertia values were obtained. Then, after using the Elbow method, the suggested number of clusters in these particles was three (K = 3). In order to confirm this result, the data were validated by means of the DBI, which showed that when comparing a K of 2, 3, and 4 (using the closest neighbors), a K = 3 gives the lowest DBI, and thus validates this result (Figure 2).
Figure 2. The CaP bioceramic particles show three micropore clusters. (A) Elbow method to determine the number of clusters in the 5,338 micropores data set, the “elbow” marked with the red asterisk suggests that the number of micropore clusters is three. (B) Davies-Bouldin Index (DBI) calculated considering the five attributes of the 5,338 micropores isolated from the five bioceramic particles, for which the K-value obtained by the elbow method, together with its closest neighbors, was used as an approximation. The lowest value of the DBI indicates the number of suitable clusters, in this case three.
The analysis of clusters was deepened by relating the Feret's diameter of the micropores with their circularity. This analysis indicates that within these three clusters of micropores, the first cluster represents <3% of the micropores, has a Feret's diameter of more than 100 μm but a circularity generally lower than 0.1, being the least circular pores. The second cluster represents <20% of the micropores, has a Feret's diameter of about 50 μm and circularity near 0.2. While the third cluster represents most of the micropores (more than 80%), has a Feret's diameter of <50 μm and a wide range of circularity up to 1, being the most circular pores (Figure 3).
Figure 3. The three micropores clusters based on Feret's diameter and circularity. (A) The three clusters of micropores are highlighted in black (1), red (2), and green (3). The yellow dots indicate the centers of each cluster. (B) Table showing the Feret's diameter and the circularity of the centers of each cluster (yellow dots in A). It also indicates the number of micropores that belong to each cluster and their representative percentage within the bioceramic particles.
Phase Identification of Cap Bioceramic Particles
In order to corroborate this CaP bioceramic crystallographic phases, their XRD pattern was indexed using standard cards by the COD. This bioceramic showed two phases, which matched with the card number PDF 86-0740, correlated to hydroxyapatite, and PDF 03-0713, correlated to Whitlockite, which corresponds to β-tricalcium phosphate (Jarcho et al., 1979).
The Rietveld composition analysis showed that this bioceramic is composed of 53.74% hydroxyapatite (Ca5(PO4)3OH) and 39.02% β-tricalcium phosphate (Ca3(PO4)2), with the rest of the components being other calcium phosphates or calcium oxide: 3.57% of CaP4O11 and 3.68% of CaO. Therefore, the main crystallographic phases of these bioceramic particles are indeed hydroxyapatite and β-tricalcium phosphate (Figure 4).
Figure 4. X-Ray diffraction pattern for the bioceramic particle used. The peaks corresponding to Hydroxyapatite and β-TCP are shown by blue squares and red circles, respectively. Peaks identification based on Lee et al. (2013).
Discussion
This research analysis shows that in these CaP bioceramic particles, three groups of micropores are distinguished by their area, perimeter, circularity, Feret's diameter and roundness. To the best of the author's knowledge, this is the first report on the analysis of microporosity clusters of CaP bioceramics particles with the described method, driven by the research question: what are the geometrical properties of individual micropores in a CaP bioceramic, and how do they relate to each other?
Given that on the one hand, the fabrication process of bioceramics leads to the generation of a complex microarchitecture (Wilson et al., 2004), which have an effect on their biological (Polak et al., 2013; Rustom et al., 2016) and mechanical properties (Pecqueux et al., 2010). And on the other hand, the generation of these large amounts of micropores requires techniques to analyse these big data sets. This study used CaP bioceramic particles as a model, to validate a novel method to analyse such microarchitecture in terms of micropore clustering, which has not been previously studied on this bioceramic (Conserva et al., 2011).
The microstructural and morphological properties of CaP bioceramics are related to biological and cellular performance, for example, it has been shown that the size and shape of HA microparticles have an influence on the osteogenic differentiation of MC3T3-E1 preosteoblasts (Xu et al., 2018). Furthermore, microporosity (pores of < 50 μm) plays a role in bone ingrowth and in the improvement of a scaffold mechanical properties (Rustom et al., 2019).
In this study, almost 20% of the micropores (clusters 1 and 2) have a Feret's diameter of 50 μm or more. In a study that used biphasic CaP bioceramics in vivo, these large pores have been shown to promote natural bone deposition and resorption processes (Lan Levengood et al., 2010b). The remaining 80% (cluster 3), corresponds to micropores of smaller size (20 μm or less). These small pores can have a higher capability to adsorb proteins in vitro, which can influence the adhesion and proliferation of human bone cells (Rouahi et al., 2006). Additionally, it has been shown that an increase in surface roughness significantly improves cell attachment and proliferation (Li et al., 2005), as well as improving the homogeneity of bone distribution by micropore-induced capillarity (Rustom et al., 2016).
The abovementioned pore properties, added to the existence of a central macropore in each of these particles, which could allow vasculogenesis, as shown for the macropores of other CaP bioceramics (Lan Levengood et al., 2011), make these particles suitable for bone tissue regeneration, with proper oxygen and blood supply. Besides, it is known that biphasic CaP bioceramics create a unique microenvironment that favors bone regeneration (Lobo and Arinzeh, 2010).
By using the presented methodology, future studies will be able to relate micropore geometry and clustering with biological performance of different CaP bioceramics in terms of adhesion, proliferation, or differentiation, among others. Consequently, an optimal set of microarchitectural parameters can be found, which will lead to a desired response in clinical therapies. Furthermore, this methodology could be applied to structures made of other materials, and fabricated by processes that could allow for major control of their microarchitecture for bone tissue engineering applications (Turnbull et al., 2017).
Within the limitations of this study, it is the surface microarchitecture of the CaP bioceramic and not the internal microarchitecture of the particle that was examined, which in future studies should be investigated to determine if there are variations in the characteristics of the micropores and, consequently, of the clusters. Likewise, it is necessary to analyze a larger number of particles, different manufacturing lots, and incorporate a more comprehensive and multi-scale analysis in order to integrate biology with the materials science paradigm (Roeder, 2010).
Conclusions
In this study, a novel method based on data mining techniques was used to analyze micropore clustering in CaP bioceramic particles, and five geometrical parameters were assessed in each micropore (area, perimeter, circularity, Feret's diameter, and roundness). As a result, three clusters of micropores are described, in which almost 20% of the micropores (clusters 1 and 2) have a Feret's diameter of 50 μm or more, with the remaining 80% (cluster 3), having micropores of smaller diameter (20 μm or less).
This new methodology can be applied to advance the knowledge about CaP bioceramics by analyzing global micropore properties, such as pore clustering, and their relation with other properties of interest, such as mechanical or biological performance, for bone tissue engineering.
Data Availability Statement
The datasets generated for this study are available on request to the corresponding author.
Ethics Statement
The studies involving human participants were reviewed and approved by Universidad Adolfo Ibáñez. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
CM and JV planned the study. IB, VI, and JV performed the experiments. All authors contributed to the results analysis and RV, IB, CM, and JV contributed to the manuscript writing.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
The authors thank Megagen for in-kind donation of bioceramic materials for analysis. RV thanks Universidad Adolfo Ibáñez for providing postdoctoral funding. The authors also thank Rocío Orellana for technical assistance with SEM. Lastly, the authors thank FONDECYT Inicio (N° 11170957) for providing funding for this research.
References
Akay, G., Birch, M. A., and Bokhari, M. A. (2004). Microcellular polyHIPE polymer supports osteoblast growth and bone formation in vitro. Biomaterials 25, 3991–4000. doi: 10.1016/j.biomaterials.2003.10.086
Bish, D. L., and Howard, S. A. (1988). Quantitative phase analysis using the Rietveld method. J. Appl. Crystallogr. 21, 86–91. doi: 10.1107/S0021889887009415
Champion, E. (2013). Sintering of calcium phosphate bioceramics. Acta Biomater. 9, 5855–5875. doi: 10.1016/j.actbio.2012.11.029
Conserva, E., Foschi, F., Cancedda, R., and Mastrogiacomo, M. (2011). In vitro and in vivo osteoinductive and osteoconductive properties of a synthetic bone substitute. Int. J. Oral Maxillofacial Implants 28, e432–e439. doi: 10.11607/jomi.te23
Davies, D. L., and Bouldin, D. W. (1979). A cluster separation measure. IEEE Trans. Pattern Anal. Mach. Intell. 224–227. doi: 10.1109/TPAMI.1979.4766909
Denry, I., and Kuhn, L. T. (2016). Design and characterization of calcium phosphate ceramic scaffolds for bone tissue engineering. Dental Mater. 32, 43–53. doi: 10.1016/j.dental.2015.09.008
Dorozhkin, S. (2010). Bioceramics of calcium orthophosphates. Biomaterials 31, 1465–1485. doi: 10.1016/j.biomaterials.2009.11.050
El-Ghannam, A. R. (2004). Advanced bioceramic composite for bone tissue engineering: design principles and structure-bioactivity relationship. J. Biomed. Mater. Res. A 69, 490–501. doi: 10.1002/jbm.a.30022
Flautre, B., Descamps, M., Delecourt, C., Blary, M. C., and Hardouin, P. (2001). Porous HA ceramic for bone replacement: role of the pores and interconnections - experimental study in the rabbit. J. Mater. Sci. Mater. Med. 12, 679–682. doi: 10.1023/A:1011256107282
Galois, L., and Mainard, D. (2004). Bone ingrowth into two porous ceramics with different pore sizes: an experimental study. Acta Orthop. Belg. 70, 598–603.
Gauthier, O., Bouler, J.-M., Aguado, E., Pilet, P., and Daculsi, G. (1998). Macroporous biphasic calcium phosphate ceramics: influence of macropore diameter and macroporosity percentage on bone ingrowth. Biomaterials 19, 133–139. doi: 10.1016/S0142-9612(97)00180-4
Giannoudis, P. V., Dinopoulos, H., and Tsiridis, E. (2005). Bone substitutes: an update. Injury 36, S20–S27. doi: 10.1016/j.injury.2005.07.029
Gonzalez, R. C., and Tou, J. T. (1974). Pattern Recognition Principles. Applied Mathematics and Computation. Reading, MA: Addison-Wesley.
Habibovic, P., and de Groot, K. (2007). Osteoinductive biomaterials—properties and relevance in bone repair. J. Tissue Eng. Regener. Med. 1, 25–32. doi: 10.1002/term.5
Habraken, W., Habibovic, P., Epple, M., and Bohner, M. (2016). Calcium phosphates in biomedical applications: materials for the future? Mater. Tdy. 19, 69–87. doi: 10.1016/j.mattod.2015.10.008
Hollister, S. J., Lin, C. Y., Saito, E., Lin, C. Y., Schek, R. D., Taboas, J. M., et al. (2005). Engineering craniofacial scaffolds. Orthod. Craniofac. Res. 8, 162–173. doi: 10.1111/j.1601-6343.2005.00329.x
Holzapfel, B. M., Reichert, J. C., Schantz, J. T., Gbureck, U., Rackwitz, L., Nöth, U., et al. (2013). How smart do biomaterials need to be? A translational science and clinical point of view. Adv. Drug Deliv. Rev. 65, 581–603. doi: 10.1016/j.addr.2012.07.009
Hutmacher, D. W. (2000). “Scaffolds in tissue engineering bone and cartilage,” in ed D. F. B. T.-T. B. S. J. C. Williams (Oxford: Elsevier Science), 175–189. doi: 10.1016/B978-008045154-1.50021-6
Jarcho, M., Salsbury, R. L., Thomas, M. B., and Doremus, R. H. (1979). Synthesis and fabrication of β-tricalcium phosphate (whitlockite) ceramics for potential prosthetic applications. J. Mater. Sci. 14, 142–150. doi: 10.1007/BF01028337
Karageorgiou, V., and Kaplan, D. (2005). Porosity of 3D biomaterial scaffolds and osteogenesis. Biomaterials 26, 5474–5491. doi: 10.1016/j.biomaterials.2005.02.002
Kaufman, L., and Rousseeuw, P. J. (1990). Finding Groups in Data: An Introduction to Cluster Analysis. New York, NY: Wiley. doi: 10.1002/9780470316801
Kim, H. D., Amirthalingam, S., Kim, S. L., Lee, S. S., Rangasamy, J., and Hwang, N. S. (2017). Biomimetic materials and fabrication approaches for bone tissue engineering. Adv. Healthcare Mater. 6:1700612. doi: 10.1002/adhm.201700612
Knaack, D., Goad, M. E., Aiolova, M., Rey, C., Tofighi, A., Chakravarthy, P., et al. (1998). Resorbable calcium phosphate bone substitute. J. Biomed. Mater. Res. 43, 399–409.
Lan Levengood, S. K., Poellmann, M. J., Clark, S. G., Ingram, D. A., Yoder, M. C., and Wagoner Johnson, A. J. (2011). Human endothelial colony forming cells undergo vasculogenesis within biphasic calcium phosphate bone tissue engineering constructs. Acta Biomater. 7, 4222–4228. doi: 10.1016/j.actbio.2011.07.006
Lan Levengood, S. K., Polak, S. J., Poellmann, M. J., Hoelzle, D. J., Maki, A. J., Clark, S. G., et al. (2010a). The effect of BMP-2 on micro- and macroscale osteointegration of biphasic calcium phosphate scaffolds with multiscale porosity. Acta Biomater. 6, 3283–3291. doi: 10.1016/j.actbio.2010.02.026
Lan Levengood, S. K., Polak, S. J., Wheeler, M. B., Maki, A. J., Clark, S. G., Jamison, R. D., et al. (2010b). Multiscale osteointegration as a new paradigm for the design of calcium phosphate scaffolds for bone regeneration. Biomaterials 31, 3552–3563. doi: 10.1016/j.biomaterials.2010.01.052
Lee, D. S. H., Pai, Y., and Chang, S. (2013). Effect of thermal treatment of the hydroxyapatite powders on the micropore and microstructure of porous biphasic calcium phosphate composite granules. J. Biomater. Nanobiotechnol. 4, 114–118. doi: 10.4236/jbnb.2013.42015
Li, J., Yun, H., Gong, Y., Zhao, N., and Zhang, X. (2005). Effects of surface modification of poly (3-hydroxybutyrate-co-3-hydroxyhexanoate) (PHBHHx) on physicochemical properties and on interactions with MC3T3-E1 cells. J. Biomed. Mater. Res. A 75, 985–998. doi: 10.1002/jbm.a.30504
Lobo, S. E., and Arinzeh, T. L. (2010). Biphasic calcium phosphate ceramics for bone regeneration and tissue engineering applications. Materials 3, 815–826. doi: 10.3390/ma3020815
Low, K. L., Tan, S. H., Zein, S. H., Roether, J. A., Mouriño, V., and Boccaccini, A. R. (2010). Calcium phosphate-based composites as injectable bone substitute materials. J. Biomed. Mater. Res. B Appl. Biomater. 94, 273–286. doi: 10.1002/jbm.b.31619
Mastrogiacomo, M., Scaglione, S., Martinetti, R., Dolcini, L., Beltrame, F., Cancedda, R., et al. (2006). Role of scaffold internal structure on in vivo bone formation in macroporous calcium phosphate bioceramics. Biomaterials 27, 3230–3237. doi: 10.1016/j.biomaterials.2006.01.031
Merkus, H. G. (2009). Particle Size Measurements: Fundamentals, Practice, Quality, Vol. 17. Springer Science & Business Media. Available online at: https://books.google.cl/books?hl=es&lr=&id=lLx4GzA-7AUC&oi=fnd&pg=PP5&dq=Merkus,+H.+G.+(2009).+Particle+Size+Measurements:+Fundamentals,+Practice,+Quality,+Vol.+17.+Springer+Science+%26+Business+Media.&ots=VbGX6I124i&sig=iQieZHJI2rSTJFjZZ9TpMCkqNvU#v=onepage&q=Merkus%2C%20H.%20G.%20(2009).%20Particle%20Size%20Measurements%3A%20Fundamentals%2C%20Practice%2C%20Quality%2C%20Vol.%2017.%20Springer%20Science%20%26%20Business%20Media.&f=false
Novotna, L., Kucera, L., Hampl, A., Drdlik, D., and Cihlar, J. (2019). Biphasic calcium phosphate scaffolds with controlled pore size distribution prepared by in-situ foaming. Mater. Sci. Eng. C 95, 363–370. doi: 10.1016/j.msec.2018.03.022
Pecqueux, F., Tancret, F., Payraudeau, N., and Bouler, J. M. (2010). Influence of microporosity and macroporosity on the mechanical properties of biphasic calcium phosphate bioceramics: modelling and experiment. J. Eur. Ceramic Soc. 30, 819–829. doi: 10.1016/j.jeurceramsoc.2009.09.017
Pham, D. T., Dimov, S. S., and Nguyen, C. D. (2005). Selection of K in K-means clustering. Proc. Instit. Mech. Eng. C J. Mech. Eng. Sci. 219, 103–119. doi: 10.1243/095440605X8298
Polak, S. J., Rustom, L. E., Genin, G. M., Talcott, M., and Wagoner Johnson, A. J. (2013). A mechanism for effective cell-seeding in rigid, microporous substrates. Acta Biomater. 9, 7977–7986. doi: 10.1016/j.actbio.2013.04.040
Roeder, R. K. (2010). A paradigm for the integration of biology in materials science and engineering. JOM 62, 49–55. doi: 10.1007/s11837-010-0108-5
Rouahi, M., Gallet, O., Champion, E., Dentzer, J., Hardouin, P., and Anselme, K. (2006). Influence of hydroxyapatite microstructure on human bone cell response. J. Biomed. Mater. Res. A 78, 222–235. doi: 10.1002/jbm.a.30682
Rueden, C. T., Schindelin, J., Hiner, M. C., DeZonia, B. E., Walter, A. E., Arena, E. T., et al. (2017). ImageJ2: ImageJ for the next generation of scientific image data. BMC Bioinformatics 18, 1–26. doi: 10.1186/s12859-017-1934-z
Rustom, L. E., Boudou, T., Lou, S., Pignot-Paintrand, I., Nemke, B. W., Lu, Y., et al. (2016). Micropore-induced capillarity enhances bone distribution in vivo in biphasic calcium phosphate scaffolds. Acta Biomater. 44, 144–154. doi: 10.1016/j.actbio.2016.08.025
Rustom, L. E., Poellmann, M. J., and Wagoner Johnson, A. J. (2019). Mineralization in micropores of calcium phosphate scaffolds. Acta Biomater. 83, 435–455. doi: 10.1016/j.actbio.2018.11.003
Schek, R. M., Wilke, E. N., Hollister, S. J., and Krebsbach, P. H. (2006). Combined use of designed scaffolds and adenoviral gene therapy for skeletal tissue engineering. Biomaterials 27, 1160–1166. doi: 10.1016/j.biomaterials.2005.07.029
Turnbull, G., Clarke, J., Picard, F., Riches, P., Jia, L., Han, F., et al. (2017). 3D bioactive composite scaffolds for bone tissue engineering. Bioactive Mater. 3, 278–314. doi: 10.1016/j.bioactmat.2017.10.001
Uswatta, S. P., Okeke, I. U., and Jayasuriya, A. C. (2016). Injectable porous nano-hydroxyapatite/chitosan/tripolyphosphate scaffolds with improved compressive strength for bone regeneration. Mater. Sci. Eng. C 69, 505–512. doi: 10.1016/j.msec.2016.06.089
Vicente, V., Meseguer, L., Martinez, F., Galian, A., Rodriguez, J., Alcaraz, M., et al. (1996). Ultrastructural study of the osteointegration of bioceramics (Whitlockite and composite β-TCP + collagen) in rabbit bone. Ultrastruct. Pathol. 20, 179–188. doi: 10.3109/01913129609016313
Vivanco, J., Aiyangar, A., Araneda, A., and Ploeg, H. L. (2012). Mechanical characterization of injection-molded macro porous bioceramic bone scaffolds. J. Mech. Behav. Biomed. Mater. 9, 137–152. doi: 10.1016/j.jmbbm.2012.02.003
Walton, W. H. (1948). Feret's statistical diameter as a measure of particle size. Nature 162, 329. doi: 10.1038/162329b0
Wilson, C. E., De Bruijn, J. D., Van Blitterswijk, C. A., Verbout, A. J., and Dhert, W. J. (2004). Design and fabrication of standardized hydroxyapatite scaffolds with a defined macro-architecture by rapid prototyping for bone-tissue-engineering research. J. Biomed. Mater. Res. A 68, 123–132. doi: 10.1002/jbm.a.20015
Wiltfang, J., Merten, H. A., Schlegel, K. A., Schultze-Mosgau, S., Kloss, F. R., Rupprecht, S., et al. (2002). Degradation characteristics of α and β tri-calcium-phosphate (TCP) in minipigs. J. Biomed. Mater. Res. 63, 115–121. doi: 10.1002/jbm.10084
Xu, C. C., and Chan, R. W. (2008). Pore architecture of a bovine acellular vocal fold scaffold. Tissue Eng. A 14, 1893–1903. doi: 10.1089/ten.tea.2007.0243
Xu, C. C., Chan, R. W., and Tirunagari, N. (2007). A biodegradable, acellular xenogeneic scaffold for regeneration of the vocal fold lamina propria. Tissue Eng. 13, 551–566. doi: 10.1089/ten.2006.0169
Xu, W., Xu, C., Yi, J., and Dai, H. (2018). The effect of different hydroxyapatite microparticles on the osteogenic differentiation of MC3T3-E1 preosteoblasts. J. Mater. Chem B. 6, 5234–5242. doi: 10.1039/C8TB01352G
Zadpoor, A. A. (2015). Bone tissue regeneration: the role of scaffold geometry. Biomater. Sci. 3, 231–245. doi: 10.1039/C4BM00291A
Keywords: calcium phosphate, bioceramic particle, microporosity, data mining, K-means clustering
Citation: Vallejos Baier R, Benjumeda Wijnhoven I, Irribarra del Valle V, Millán Giovanetti C and Vivanco JF (2019) Microporosity Clustering Assessment in Calcium Phosphate Bioceramic Particles. Front. Bioeng. Biotechnol. 7:281. doi: 10.3389/fbioe.2019.00281
Received: 13 August 2019; Accepted: 04 October 2019;
Published: 18 October 2019.
Edited by:
Francesco Baino, Politecnico di Torino, ItalyReviewed by:
Francesca Diomede, Università degli Studi G. d'Annunzio Chieti e Pescara, ItalyD'Arros Cyril, INSERM U1229 Regenerative Medicine and Skeleton, France
Sergey V. Dorozhkin, Independent Researcher, Moscow, Russia
Copyright © 2019 Vallejos Baier, Benjumeda Wijnhoven, Irribarra del Valle, Millán Giovanetti and Vivanco. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Juan F. Vivanco, juan.vivanco@uai.cl