# 3D printing of cellulose nanocrystals based composites to build robust biomimetic scaffolds for bone tissue engineering

Kanga Marius N'gatta, Habib Belaid, Joelle El Hayek, Edja Florentin

Assanvo, Marilyn Kajdan, Nathalie Masquelez, David Boa, Vincent Cavaillès,

Mikhael Bechelany, Chrystelle Salameh

► To cite this version:

Kanga Marius N'gatta, Habib Belaid, Joelle El Hayek, Edja Florentin Assanvo, Marilyn Kajdan, et al.. 3D printing of cellulose nanocrystals based composites to build robust biomimetic scaffolds for bone tissue engineering. Scientific Reports, 2022, 12 (1), pp.21244. 10.1038/s41598-022-25652-x. hal-04109785

HAL Id: hal-04109785

https://hal.umontpellier.fr/hal-04109785

Submitted on 30 May 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

OPEN

# 3D printing of cellulose nanocrystals based composites to build robust biomimetic scaffolds for bone tissue engineering

Kanga Marius N'Gatta<fnref n="1" />,², Habib Belaid<fnref n="1" />,³, Joelle El Hayek<fnref n="1" />, Edja Florentin Assanvo², Marilyn Kajdan³, Nathalie Masquelez<fnref n="1" />, David Boa², Vincent Cavaillès³, Mikhael Bechelany<fnref n="1" /> & Chrystelle Salameh<fnref n="1" />✉

Cellulose nanocrystals (CNC) are drawing increasing attention in the fields of biomedicine and healthcare owing to their durability, biocompatibility, biodegradability and excellent mechanical properties. Herein, we fabricated using fused deposition modelling technology 3D composite scaffolds from polylactic acid (PLA) and CNC extracted from *Ficus thonningii*. Scanning electron microscopy revealed that the printed scaffolds exhibit interconnected pores with an estimated average pore size of approximately 400 µm. Incorporating 3% (w/w) of CNC into the composite improved PLA mechanical properties (Young's modulus increased by ~30%) and wettability (water contact angle decreased by ~17%). The mineralization process of printed scaffolds using simulated body fluid was validated and nucleation of hydroxyapatite confirmed. Additionally, cytocompatibility tests revealed that PLA and CNC-based PLA scaffolds are non-toxic and compatible with bone cells. Our design, based on rapid 3D printing of PLA/CNC composites, combines the ability to control the architecture and provide improved mechanical and biological properties of the scaffolds, which opens perspectives for applications in bone tissue engineering and in regenerative medicine.

The development of materials for bone engineering remains a challenge due to the complexity of the natural bone structure and the biomechanical environment. To repair damaged bone tissue, autografts from different bones are harvested and used to replace missing bones. Less available allografts are discouraging<fnref n="1" />⁻³, and artificial junctions used as implants often have to be removed after healing⁴. Recently, new bone repair strategies have emerged, including scaffold-assisted regenerative medicine used to promote bone growth⁵.

An ideal bone scaffold should be a three-dimensional matrix capable of mimicking the complex composition and structure of the bone for cell attachment and proliferation⁶. Therefore, it requires high biocompatibility, biodegradability, non-toxicity, excellent mechanical properties and appropriate architecture in terms of porosity and pore sizes to integrate with native host tissue⁷. The chemical composition and physicochemical characteristics of the scaffold, directly influencing the mechanical and biological performance⁸, are therefore important parameters to study.

Synthetic biopolymers have been widely used in bone tissue engineering due to their biocompatibility and their ability to control the physicochemical properties of the scaffold. They consist of aliphatic polyesters such as polyglycolic acid (PGA), polycaprolactone (PCL) and polylactic acid (PLA)⁹. Unfortunately, they are rather brittle and usually lose their strength due to rapid degradation in vivo. Moreover, their hydrophobic nature hinders the attachment and proliferation of bone cells¹⁰. To overcome these limitations, scaffolds based on synthetic polymers, namely PLA or PCL can be improved by incorporating natural polymer reinforcements such as cellulose¹¹⁻¹⁵, alginate¹⁶, gelatine¹⁷, chitosan¹⁸,¹⁹ or keratin²⁰ known for their interesting characteristics.

<footnote n="1">¹Institut Européen des Membranes, IEM, UMR 5635, Univ Montpellier, ENSCM, CNRS, Montpellier, France. ²Laboratoire de Thermodynamique et de Physico-Chimie du Milieu, UFR SFA, Université Nangui Abrogoua, 02 BP 801, Abidjan 02, Côte d'Ivoire. ³IRCM, Institut de Recherche en Cancérologie de Montpellier, INSERM U1194, Université Montpellier, 34298 Montpellier, France. ✉email: chrystelle.salameh@enscm.fr</footnote>
<figure><img src="image_3.png" /><figcaption>**Figure 1.** SEM images of extruded PLA and PLA/CNCx composite filaments.</figcaption></figure>

Several techniques, including solvent casting and particle leaching, emulsion freeze-drying, phase separation, or electrospinning<sup>21-23</sup>, are used to develop scaffolds for hard tissue engineering. However, they do not allow efficient control of the morphology and porosity.

Additive manufacturing has proven to be a technique of choice for designing and preparing biomimetic bone repair materials. 3D-controlled scaffold architectures significantly affect mechanical properties as well as bone cell adhesion and proliferation<sup>2,24-28</sup>. Therefore, various works have focused on the development of 3D printed scaffolds using various technologies, such as stereolithography, 3D plotting, selective laser sintering, bioprinting, and fused deposition modelling (FDM). FDM is the most widely used additive manufacturing technology. It is a simple and fast technique at low cost offering great possibilities for handling polymers.

Several studies have investigated the reinforcement of PLA with nanomaterials such as graphene oxide<sup>29</sup>, boron nitride (BN) nanofillers<sup>21</sup>, nano-hydroxyapatite<sup>30-32</sup> to improve the mechanical and bioactive properties of the polymeric scaffold. Cellulose is a promising material for biomedical applications, including tissue engineering, stem cell research and regenerative medicine<sup>33-37</sup> due to its excellent biocompatibility, favorable biodegradability, good mechanical properties, high specific surface area, low density and non-abrasive nature<sup>38-41</sup>. In particular, cellulose nanocrystals (CNC) are effective to be applied as reinforcing agents or nanofillers in scaffold design and the results confirmed that CNC significantly improve the mechanical performance and cytocompatibility of scaffolds<sup>38,39,42,43</sup>. In addition, CNC have many hydroxyl groups which can form intermolecular hydrogen bonds with carbonyl groups of PLA to achieve favorable interfacial adhesion. Blending PLA with biopolymers such as cellulose improves the properties of PLA without altering its biocompatibility. PLA/cellulose composites have been widely developed and studies on structural, thermal and mechanical properties have been reported<sup>12,44-45</sup>. However, even though PLA and cellulose nanocomposites are reported in the literature for tissue engineering applications, we believe that there is still room to develop 3D scaffolds for tissue engineering by FDM which is a rapid and low-cost additive manufacturing technique especially to overcome certain limitations encountered in the PLA, namely the mechanical resistance of the biopolymer.

This work therefore aims to unravel the effect of CNC content on PLA/CNC scaffolds processed by FDM and its influence on the final structural, surface and mechanical properties of the biocomposites. Biomineralization, in vitro biodegradability and cytocompatibility tests were performed on the resulting PLA/CNC3 (i.e. with 3% (w/w) CNC incorporated into the composite). The resulting material is a mechanically enhanced scaffold that offers great opportunities for the design and rapid production of biomimetic 3D scaffolds with appropriate biological properties.

# Results and discussion

**Characterization of PLA and PLA/CNCx filaments.** Scanning Electron Microscopy (SEM) images reported in Fig. 1 clearly show pure PLA filaments with a smooth and regular surface while those of the composites are rough due to the formation of aggregates related to the presence of CNC. Indeed the evolution of the PLA surface, after the incorporation of CNC, from a smooth state to a rough state could suggest a good dispersion of CNC over the entire surface of PLA. The incorporation of more than 3% (w/w) of CNC results in the formation of large particles, which could lead to less hydrogen bonds between the matrix (PLA) and the reinforcement (CNC). Indeed, the dispersion of CNC in PLA is low due to the hydrophilic nature of cellulose *versus* the hydrophobic nature of PLA rendering their interaction and uniform dispersion difficult<sup>52</sup>.

The composite filaments, resulting from the hybridization of PLA, exhibit good capacities for improving cell binding and proliferation. In this respect, PLA/CNC composites are good candidates for biomedical applications, particularly in bone regeneration. Indeed, cellulose fillers (CNC) could create important cell attachment and proliferation sites for any implantable scaffold material in hard tissue engineering. Several studies have already been carried out to develop PLA/CNC composite filaments by single-screw extrusion printable by FDM<sup>52-54</sup>. However, the development of PLA/CNC filaments for the production of hybrid scaffolds for biomedical applications remains new to our knowledge.

Fourier-transform infrared spectroscopy (FTIR) analysis of the filaments are shown in Figure S1 in Supplementary Information. Our results show a slight increase in peak intensities with the increase of the CNC loading within the PLA matrix, which is in agreement with the literature<sup>52</sup>. The most significant differences were
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Thermal analysis of PLA and PLA/CNCx extruded filaments. (a) thermogravimetric analysis (TGA, under nitrogen from room temperature to 800 °C at a heating rate of 10 °C min⁻¹), (b) derivative TGA (DTG) and (c) differential scanning calorimetry (DSC curves, from room temperature to 250 °C with a heating rate of 20 °C min⁻¹ under nitrogen atmosphere) showing the T<sub>g</sub>, T<sub>cc</sub> and T<sub>m</sub> of PLA and the nanocomposites.</figcaption></figure>

<table><caption>**Table 1.** Glass transition temperature (Tg), cold crystallization temperature (Tcc), melt temperature (Tm), cold crystallization enthalpy (ΔHcc)) and degree of crystallinity (χc) for PLA and PLA/CNCx composites from DSC.</caption><thead><tr><th>Sample</th><th>T<sub>g</sub> (°C)</th><th>T<sub>cc</sub> (°C)</th><th>T<sub>m</sub> (°C)</th><th>ΔH<sub>cc</sub> (J/g)</th><th>χ<sub>c</sub> (%)</th></tr></thead><tbody><tr><td>PLA</td><td>66</td><td>100</td><td>171</td><td>24</td><td>26</td></tr><tr><td>PLA/CNC1</td><td>57</td><td>88</td><td>168</td><td>20</td><td>22</td></tr><tr><td>PLA/CNC3</td><td>56</td><td>89</td><td>167</td><td>22</td><td>27</td></tr><tr><td>PLA/CNC5</td><td>57</td><td>87</td><td>166</td><td>23</td><td>26</td></tr></tbody></table>

observed at 2946, 1743, and between 1500 and 1000 cm⁻¹ which are respectively attributed to stretching of CH, C=O bonds, carbonyl stretching vibration (C=O) and to CO of PLA⁵⁵.

The decomposition of CNC, PLA filaments and the derived composites was examined from the thermogravimetric analysis (TGA), derivative thermogravimetry (DTG) and differential scanning calorimetry (DSC) curves reported in Fig. 2. A slight weight loss is observed between 100 and approximately 325 °C in all samples due to the dehydration of filaments and depolymerization of PLA as well as the breaking of glycosidic bonds in cellulose. Compared to PLA, the composites exhibit higher weight loss due to the hygroscopic nature of the composites after incorporation of hydrophilic CNC. Therefore, when subjected to thermal stress, composites release the absorbed moisture, resulting in this relatively high weight loss. In addition, a strong weight loss occurs at 400 °C probably related to the degradation of PLA and CNC. The onset decomposition temperatures of CNC and the filaments are observed at 225, 318.8, 311.8, 325.5 and 320.4 °C for CNC, PLA, PLA/CNC1, PLA/CNC3 and PLA/CNC5, respectively. The developed filaments can therefore be processed at temperatures reaching 300 °C without damaging them. Therefore, the printing temperature of 200 °C used in our study, seems to be optimal for developing PLA/CNC scaffolds. However, CNC do not significantly improve the thermal stability of PLA. Moreover, the DTG curves show a drop in the maximum degradation temperature of the composites compared to that of pure PLA. The maximum degradation of PLA/CNC3 occurs over a temperature (T) range of approximately 365–374 °C while that of PLA, PLA/CNC1 and PLA/CNC5 occurs at T=373.39, 370.62 and 369.22 °C, respectively. Indeed, the relatively wide degradation temperature of PLA/CNC3 can be attributed to the better hybridization of the composite reflected by a better interaction between 3% of CNC and PLA.

Furthermore, the analysis of the DSC curves (Fig. 2c) of all the samples successively reveals endothermic stage changes (glass transition), followed by an exothermic peak dedicated to cold crystallization and an endothermic peak due to the fusion. Compared to PLA, the composite filaments show a shift of the glass transition, cold crystallization and melting peaks towards the lowest temperatures. This suggests that the CNC tends to decrease the change-of-state temperatures of PLA (Table 1). Indeed, the glass transition range widens with the increase of CNC loading. This suggests that the addition of cellulose increases crystallinity due to the enhanced mobility of the PLA chains⁵⁶. This is consistent with the calculations of the degree of crystallinity presented in Table 1. The addition of 3% (w/w) CNC into the PLA matrix caused a wider glass transition range with a weaker peak, which indicates better interaction between CNC and PLA and therefore could reflect a good dispersion and good interfacial adhesion.

In the Supplementary Information, Figures S2 and S3 display the thermal analysis of CNC and of 3D printed PLA and PLA/CNC3 respectively, which proves that PLA exhibits no thermodegradability when printed.
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Mechanical properties: (a) Photo of dog bone shaped PLA and PLA/CNCx scaffolds for tensile strength testing. (b) Young's modulus. (c) Elasticity limit tensile. (d&e) stress at break (***p < 0.001, **p < 0.01, ns = not significant).</figcaption></figure>

**Mechanical properties of 3D printed scaffolds.** Implantable scaffolds for biomedical applications require specific attributes to best perform native organ functions and help regenerate damaged tissue. Among their attributes, the mechanical properties remain the most crucial to support the regeneration process. Several authors have investigated scaffolds combining various materials to provide successful substitution or mechanical support needed to promote new tissue growth at defect sites. However, the mechanical properties of hybrid scaffolds depend on the method of manufacture, the content, the orientation of the nanomaterials, the nature of the matrix, the interactions between the polymer matrix and the nanomaterials, etc<sup>21,29,57-60</sup>

Different filament compositions were used for 3D printing of PLA/CNC scaffolds using FDM as described in the experimental section. The analysis of the effect of cellulose on Young's modulus, elastic limit and strain at break of the 3D printed scaffolds was studied and the results are shown in Fig. 3. PLA scaffolds have an average Young's modulus of 2.4 ± 0.1 GPa in accordance with the literature<sup>29,61</sup> and an average elongation at break of 2.4 ± 0.3% while the average Young's modulus and elongation at break of the composites range between 2–3 GPa and 2–3.6%, respectively. This is in good agreement with the properties of PLA/natural fiber composites<sup>62</sup>. It seems that the addition of 3% of CNC enhances the mechanical strength of PLA (*p* < 0.001). This improvement may be due to better hybridization of CNC particles in PLA chains to ensure the best compatibility and absorb any external load. Indeed, at a relatively low concentration of cellulosic reinforcement, the molecular interactions between PLA and CNC provide good rigidity to the composites. This results in an improvement in Young's modulus<sup>63</sup> from 2.4 ± 0.1 GPa for PLA to 3.1 ± 0.2 GPa for PLA/CNC3 because at low concentrations of CNC the particles tend to align in the direction of the orientation of the PLA chains<sup>52</sup>. Consequently, the applied stress is transferred from one fibril to the next within the composite structure, thus allowing a uniform distribution of the stress in the material. The improvement in elongation at break and deformation forces at break of the composites with cellulosic reinforcement ≤ 3% suggests that the presence of CNC tends to improve the elastoplasticity of PLA. This could result from the good fiber-matrix interaction due to a better interlocking of the CNC network in the PLA matrix. The incorporation of 5% of CNC tends to form large particles distributed in the composite, which results in the formation of fewer bonds between the matrix (PLA) and the reinforcement (CNC) and poor compatibility resulting in poor mechanical performance.

Based on these results, we focused the following analyses on PLA/CNC3 only. Indeed, when the CNC content is more than 3% (w/w), the PLA cannot sufficiently wet the surface of the CNC, resulting in strong fiber-fiber interaction and/or fiber agglomeration in the PLA matrix<sup>52,54</sup>. However, Kumar et al.<sup>52</sup> found a 50% increase in the elastic modulus (E) with 1% CNC (4550 MPa) compared to pure PLA (3030 MPa). The elastic modulus of PLA and its composites was better than ours. Nevertheless, there is a decrease in the strain at break observed for 1% CNC (2.8%) compared to the pure PLA sample (8.7%).

**Morphology of 3D printed scaffolds.** The surface morphology of PLA and PLA/CNC3 scaffolds was analyzed by SEM, 3D optical microscopy and topography (Fig. 4). These images clearly show that the PLA scaffold has a smooth surface while that of PLA/CNC3 is rough due to the presence of cellulose. Ideally, the printed scaffolds should have pores of the same size. However, analysis of SEM images (Fig. 4) reveals that the pore size
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Scanning electron micrographs of 3D printed scaffolds: Top view of the morphology of (a) PLA and (c) PLA/CNC3. Cross-section with a zoom in inset of (b) PLA and (d) PLA/CNC3. (g) Statistical analysis of the 3D printed PLA and PLA/CNC diameters observed on top of the samples by SEM (Test Anova, **p**<0.01). Topography of 3D printed scaffolds typically showing the surface roughness of (e) PLA and (f) PLA/CNC3.</figcaption></figure>

of the hybrid scaffold (PLA/CNC3) increased while the width of the rods decreased similarly to what is reported in the literature²¹,²⁹. In our work, the average pore size and the rod width of PLA and PLA/CNC3 scaffolds are respectively evaluated at (362 ± 19 µm and 458 ± 17 µm) and (450 ± 15 µm and 337 ± 30 µm). This variation in the morphology of the scaffold surface can be attributed to the modification of the rheological behavior of the matrix by the addition of CNC. Indeed, the addition of CNC could limit the elastic behavior of PLA and therefore give it stability to retain its shape during deposition. However, the reinforcement produces an instability of the material flow during printing probably affecting the quality and resolution.

Moreover, porosity is one of the most critical parameters to evaluate the efficiency and applicability of bio-
medical scaffolds, particularly in bone tissue engineering, drug delivery, etc. A high porosity promotes the effi-
cient release of biofactors including cells, genes and proteins and provides an environment conducive to nutrient
exchange⁶⁴,⁶⁵. The porosity of PLA and PLA/CNC3 scaffolds is relatively high (58 ± 2% and 67 ± 2%), respectively.
The addition of 3% CNC increases the porosity of PLA scaffolds (**p**<0.01), which can be attributed to the ionic
interaction between the different components creating more voids in the scaffolds. With such porosity profiles,
printed scaffolds could help improve bone ingrowth by allowing optimal vascularization, hydroxyapatite nuclea-
tion and mineral maturation⁶⁶,⁶⁷.

**Contact angle and swelling of 3D printed scaffolds.** Hydrophilicity is a property that significantly affects the adhesion, cell proliferation and rehydration capacity of scaffolds; therefore, it appears crucial for the performance of scaffolds for any application in tissue engineering including bone regeneration. Thus, the reactivity and interaction between the printed scaffolds and the surrounding surface could be predicted by measuring the contact angle and studying the swelling in water. Figure S4 (see Supplementary Information) illustrates the water contact angles of PLA and PLA/CNC3 scaffolds and the values are presented in Table 2. CNC brings a hydrophilic character to the hydrophobic surface of PLA scaffolds due to the rearrangement of CNC in the PLA matrix. To better understand the effect of cellulose filler on the hydrophilicity of PLA, we evaluated the water absorption rate of scaffolds (Table 2). The addition of CNC leads to an improvement in the water uptake of PLA according to the results of the contact angle and surface roughness of the scaffolds. Similar water retention results were observed by *Murphy* et al.⁵⁶ on PLA and PLA/cellulose filaments. They revealed a significant
<table><caption>**Table 2.** Contact angle and swelling index of PLA and PLA/CNC scaffolds.</caption><thead><tr><td>Materials</td><td>PLA</td><td>PLA/CNC3</td></tr></thead><tbody><tr><td>Contact angle (°)</td><td>101 ± 1</td><td>84 ± 2</td></tr><tr><td>Swelling index (%)</td><td>4.1 ± 0.3</td><td>9.8 ± 0.7</td></tr></tbody></table>

<figure><img src="image_7.png" /><figcaption>**Figure 5.** Weight loss of PLA and PLA/CNC3 scaffolds as a function of degradation time (scaffolds immersed in alcalase buffer over a period of 28 days). (\*\*p < 0.01 and \*\*\*p < 0.001 using an unpaired-t-test analysis).</figcaption></figure>

increase in water uptake during the first 24 h followed by a stabilization to gradually reach a plateau (1.2% for the PLA filament and 1.4% for the PLA/cellulose composite filament).

In view of these results, PLA/CNC scaffolds are capable of creating an environment favorable to cell proliferation and therefore appear as the composites of choice for bone regeneration.

**Biodegradation.** The enzymatic degradation of biomedical scaffolds is routinely examined by determining the weight loss. In our study, PLA and PLA/CNC3 scaffolds were subjected to an in vitro degradation study by immersing them in alcalase buffer over a period of 28 days. According to the results shown in Fig. 5, it appears that the PLA scaffold has a slower degradation compared to PLA/CNC3 composite. Indeed, PLA and PLA/CNC3 have degradation rates evaluated respectively at 1.7 ± 0.3% and 3.9 ± 0.1% after 28 days of immersion in the alcalase buffer. They are therefore considered stable. The incorporation of CNC into the PLA matrix relatively improves the degradation of PLA. Compared to studies by Belaid et al.²⁹ on PLA and graphene oxide (PLA/GO) composite scaffolds, CNC loading leads to a moderate improvement in the degradation rate of PLA scaffolds in alcalase buffer. Thus, CNC appear as a suitable reinforcement of PLA. Therefore, PLA/CNC biocomposite can be designed as a potential material to develop scaffolds that can be used as temporary substitutes, especially for bone tissue regeneration. PLA/CNC composites could exhibit favorable transient properties by resorbing over time to provide space for newly grown bone tissue that will replace the scaffold in the body.

**Mineralization assays.** The biomimetic mineralization process of printed scaffolds in SBF was chosen to evaluate the formation of apatite on the surface of PLA and PLA/CNC3 in order to conclude on their bone regeneration and binding capacity. Figure 6 shows SEM images of the surface of PLA and PLA/CNC3 scaffolds after 14 days of incubation in the SBF solution. The process of mineral growth on the surface of the scaffolds was successfully achieved by forming a bioactive coating. For all the samples, the formation of an inhomogeneous mineral layer is clearly observed, reflected by the irregular distribution of aggregates of nanocrystals on the surface, the size and density of which are dependent on the incubation time and the composition of the scaffold. Several authors have observed similar results⁶⁸⁻⁷³. The nucleated mineral particles on the surface of the scaffolds have a rod-like aspect. The mineralization of apatite crystals depends on the pH of the medium, the adsorption and release of ions at the interface and the surface wettability⁷¹. Overall, mineralization increases with the incubation time and is larger and faster for PLA/CNC3 scaffolds. Indeed, the deposition of apatite on the surface of PLA is relatively slow due to its slow hydrolysis in SBF solution⁶⁹. Hydrolysis of PLA in SBF solution can take several weeks to form new carboxyl (–COOH) and hydroxyl (–OH) groups on its surface. Then the carboxyl groups undergo partial dissociation to give carboxylate ions (COO⁻) on the surface causing its negative charge. The strong electrostatic interaction between COO⁻ and Ca²⁺ ions as well as the strong hydrogen bonds between
<figure><img src="image_8.png" /><figcaption>**Figure 6.** Evaluation of the apatite growth on PLA and PLA/CNC3 scaffolds mineralized in SBF at 14 days by scanning electron microscopy at high and low magnification (SEM) to validate the formation of the inorganic layer and by elemental mapping to show the presence and homogeneous distribution of the phosphate and calcium elements that constitute the apatite. (a) SEM images and elemental mapping of PLA, (b) SEM images and elemental mapping of PLA/CNC3.</figcaption></figure>

carboxyl groups and PO₄³⁻ ions will then drive the process of accumulation of calcium Ca²⁺ and phosphate ions PO₄³⁻⁻⁶⁹,⁷⁴.

On the other hand, we might speculate that the addition of CNC to the PLA matrix could modify the mineral layer formation process due to hydrophilicity resulting from the formation of additional hydroxyl groups on the PLA surface. These functional groups play an important role in the ion exchange between the surface and the surrounding liquid and in the deposition of apatite on the scaffold. Indeed, these hydroxide groups bring additional negative charges to the surface of the scaffold to interact with the Ca²⁺ and PO₄³⁻ ions of the SBF thus promoting the adsorption and the nucleation of apatite⁶⁸,⁷¹,⁷⁵,⁷⁶. In view of these results, PLA/CNC composites are good candidates to support apatite formation and therefore contribute to bone regeneration and bonding by providing convenient microenvironment for cell growth and repair function.

**Cytocompatibility tests.** The biocompatibility of the PLA/CNC scaffold was investigated using hFOB1.19 human osteoblasts⁷⁷. Compared to untreated control cells (Ctrl), incubation with PLA and PLA/CNC3 filaments showed no inhibitory effect on the cell viability at day 5 (Fig. 7a) confirming that PLA and CNC are cytocompatible.

To determine that bone cells were able to adhere and proliferate on the corresponding scaffolds, we performed scanning electron microscopy after seeding the hFOB cell onto the scaffolds. As shown in Fig. 7b, hFOB cells fully spread onto the PLA and PLA/CNC3 scaffolds and exhibited normal morphology. Additionally, we also seeded bioluminescent hFOB cells (osteoblastic cell line) onto the scaffolds and monitored their growth over time. After 6 days, the osteoblasts proliferate significantly on both scaffolds and no significant differences were observed on PLA/CNC3 scaffold as compared to PLA (Fig. 7c). Finally, we assess the effect of PLA and PLA/CNC3 on hFOB differentiation by measuring by RT-QPCR the expression of genes (ALP, RUNX2 and OCN) used as markers of bone differentiation. As shown in Fig. 7d, the level of expression of the 3 genes was higher in the presence of PLA/CNC3 than in the presence of PLA alone. Altogether, these data suggest that PLA/CNC3 scaffolds fully support bone cell growth and are biocompatible.

Overall, our approach in designing 3D biomimetic scaffolds for bone tissue engineering is in agreement with previous investigations demonstrating that reinforcing PLA with nanomaterials improves the biological and mechanical properties of the material²¹,²⁹. Nevertheless, the existing nanofillers are synthetic and could trigger undesired effects if they accumulate in the tissues. Here, a natural biopolymer with relevant biological responses, i.e., CNC, incorporated into the PLA matrix at an optimal amount, significantly improves the mechanical strength of the scaffold without altering the tunable porosity and structure. This is a key parameter to achieve optimal functionality of state-of-the-art biomimetic tissue scaffolds.

**Conclusion**

PLA scaffolds reinforced with cellulose nanocrystals have been successfully designed and manufactured using FDM technology. Their structure, surface, morphological and mechanical properties were evaluated as well as their biological performance. The 3D cellular scaffolds are highly porous with evident presence of an interconnected repetitive porous architecture favorable to the enhancement of bone tissue growth by allowing optimal vascularization, hydroxyapatite nucleation and mineral maturation. An incorporation of 3% (w/w) of CNC in
<figure><img src="image_9.png" /><figcaption>**Figure 7.** Cytocompatibility of PLA and PLA/CNC filaments and scaffolds. (a) Effect of PLA or PLA/CNC3 filaments on hFOB human osteoblasts viability at day 5 of culture. (b) Scanning electron microscopy showing hFOB cells grown on the PLA (left panel) and PLA/CNC3 (right panel) scaffolds (scale bar = 200 µm). (c) Growth of luminescent hFOB cells on the PLA and PLA/CNC3 scaffolds at days 1, 2 and 4. (d) Expression of osteogenic differentiation marker genes on hFOB cells grown for 2 days in presence of PLA and PLA/CNC3 eluates. Statistical analysis has been performed using the Mann-Whitney test between PLA and PLA/CNC3 (ns, not significant; ***, *p* < 0.001).</figcaption></figure>

the composite increases the mechanical properties of PLA and the surface roughness, while contributing to the improvement of the wettability by a surface transition of PLA from hydrophobic to hydrophilic. The mineralization process of the printed scaffolds using simulated body fluid and the nucleation of hydroxyapatite were confirmed. Finally, PLA/CNC composites appeared fully compatible with human osteoblasts and allowed their adhesion and proliferation. Our results shed light on the relevance of CNC in PLA-based composites for tissue engineering and open simple and promising avenues to design multifunctional 3D printed scaffolds that are active for future tissue regeneration applications.

## Methods

**Materials.** Polylactic acid (PLA) pellets were purchased from NatureWorks LLC. Dichloromethane (CH₂Cl₂, <99.7%, CAS 75-09-2, Sigma-Aldrich), sodium hydroxide (NaOH, pellets, 98.9%, Sigma-Aldrich), sulfuric acid (H₂SO₄, 95.0-98.0%, CAS 7664-93-9, Sigma-Aldrich), hydrogen peroxide (H₂O₂, 30% (w/w), CAS 7722-84-1, Sigma-Aldrich), ethanol (96% vol, CAS 64-17-5, Sigma-Aldrich), sodium chloride (NaCl, ≥ 99%, CAS 7747-14-5, Sigma-Aldrich), sodium hydrogen carbonate (NaHCO₃, ≥ 99.7%, CAS 144-55-8, Sigma-Aldrich), potassium chloride (KCl, ≥ 99.7%, CAS 7447-40-7, Sigma-Aldrich), di-potassium hydrogen phosphate trihydrate (K₂HPO₄·3H₂O, ≥ 99%, CAS 16788-57-1, Sigma-Aldrich) magnesium chloride hexahydrate (MgCl₂·6H₂O, ≥ 99%, CAS 7791-18-6, Sigma-Aldrich), calcium chloride (CaCl₂, ≥ 97%, CAS 10043-52-4, Sigma-Aldrich), sodium sulfate (Na₂SO₄, ≥ 99.0%, CAS 7757-82-6, Sigma-Aldrich), Tris-hydroxymethyl aminomethane ((HOCH₂)₃CNH₂) (tris), ≥ 99.8%, CAS 1185-53-1, Sigma-Aldrich), hydrochloric acid (HCl, 36% CAS 7647-01-0, Sigma-Aldrich), Bacillus licheniformis (alcalase, CAS 126741 Sigma-Aldrich), L-cysteine (CAS,52-90-4), Azide sodium (CAS 26628-22-8, Sigma-Aldrich), TRIS buffer pH 9.0 (CAS 77-86-1, Alfa Aesar), Dimethyl sulfoxide (DMSO, CAS 23486-297; BDH Prolab), Dulbecco's Modified Eagle Medium α (DMEM/F12, CAS 10565018; Gibco)), 3-(4,5-Dimethylthiazol-2-yl)-2,5-Diphenyltetrazolium Bromide (MTT, 98%, CAS 298-93-1, Sigma-Aldrich), Fetal bovine serum (FBS, CVFSVF00-01; Eurobio), Trypsin-EDTA (CAS 25300-054; Gibco), hexamethyldisilazane (CAS 999-97-3, Sigma-Aldrich), were obtained and used without further purification. hFOB osteoblastic cell line obtained from ATCC was used for cytocompatibility analysis.

**Cellulose nanocrystals synthesis.** Cellulose nanocrystals (CNC) were synthesized from *Ficus thonningii* (FT) by a three-step method based on (1) alkali treatment, (2) bleaching and (3) acid hydrolysis as proposed
<table><caption>**Table 3.** Formulation of PLA/CNC composites.</caption><thead><tr><td>Notation</td><td>Composition</td><td>PLA weight (g)</td><td>CNC weight (g)</td></tr></thead><tbody><tr><td>PLA</td><td>100% PLA + 0% CNC</td><td>30</td><td>0.0</td></tr><tr><td>PLA/CNC1</td><td>99% PLA + 1% CNC</td><td>29.7</td><td>0.3</td></tr><tr><td>PLA/CNC3</td><td>97% PLA + 3% CNC</td><td>29.1</td><td>0.9</td></tr><tr><td>PLA/CNC5</td><td>95% PLA + 5% CNC</td><td>28.5</td><td>1.5</td></tr></tbody></table>

by Dingyuan Zheng<sup>78</sup> with modifications. Briefly, FT raw fibers, previously dried and stored at 70 °C, were cut, crushed and then dispersed in an alkaline solution of sodium hydroxide (NaOH, 0.5 N) with FTFs/NaOH fiber to solution ratio of 3 g/80 mL. They were stirred continuously at 100 °C for 4 h. After the alkali treatment, the fibers were filtered, washed with 18 M MilliQ water until pH neutralization. The resulting fibers were bleached with hydrogen peroxide (H₂O₂) for 1 h at 80 °C. The bleaching sequence was repeated several times to obtain brilliantly white fibers. Finally, nanocellulose synthesis was performed by the acid hydrolysis process followed by ultrasound sonication. In details, 10 g of bleached fibers were dissolved in 200 g of an aqueous solution of sulfuric acid (64%) with magnetic stirring for 1 h 30 min at 45 °C. Then, the suspension obtained was filtered, washed with 18 M MilliQ water several times and was then subjected to dialysis for 7 days (in bags of dialysis from 12,000 to 14,000 Da). After dialysis, the resulting solution was sonicated in an ice bath using an ultrasonic homogenizer (BANDELIN electronic, GM 3100) fitted with a 7 cm probe tip. The operating range of the ultrasonic homogenizer was set at 50%. The resulting suspension was centrifuged at 10,000 rpm for 30 min and the paste was frozen, lyophilized and then ground using a Moulinex mixer (LM935) at 8 000 rpm for 3 min to obtain a white powder of cellulose nanocrystals (CNC). Figure S5 summarizes the steps of the cellulose nanocrystals synthesis from FT stem (See Supplementary Information).

**Filaments extrusion and 3D printing of PLA/CNC scaffolds.** In order to produce homogeneous PLA/CNC composite filaments by single-screw extrusion, we initially proposed to formulate composites by incorporating nanocellulose fillers into the PLA matrix by chemical dissolution. 30 g of PLA were dissolved in 200 mL of dichloromethane (DCM) at room temperature with continuous magnetic stirring followed by sonication in an ice bath for 5 min. Then, the suspension of CNC/DCM was added to that of PLA to give blend slurries of PLA/CNC/DCM with 0, 1, 3 and 5% (w/w) of CNC relative to PLA (Table 3). The PLA/CNC/DCM suspensions obtained were sonicated in an ice bath using an ultrasonic homogenizer for 15 min and kept under continuous magnetic stirring at room temperature for 24 h to ensure good homogenization. They were then left to dry overnight at room temperature, then cut into small pieces and dried under vacuum for 24 h to ensure complete drying. The resulting dried composite polymers were extruded using a Noztek model MHB26234 single screw extruder at a temperature range of 160–165 °C. They were then cooled by conversion at the nozzle outlet using the extruder fan to give PLA/CNC hybrid filaments denoted PLA/CNCx, where x corresponds to the content of CNC. Finally, filaments of approximately 1.75 ± 0.05 mm in diameter were selected for 3D printing. The pure PLA filaments were produced by direct extrusion from PLA scoops under the same conditions as the composite filaments. Figure 8 summarizes the preparation of PLA/CNC composite filaments. The extrusion sequence was repeated 2 to 3 times to obtain PLA/CNCx filaments of homogeneous diameters. The scaffold was modeled by Computer Aided Design (CAD) using Design Spark Mechanical software, translated to STL file and sliced by Prusa3D slicer software to generate the G-code. Finally, the cylindrical designs (porous) (diameter (10 mm) × height (2 mm)) and the rectangular pieces (dense) were designed using an FDM Prusa MK2S Research model printer equipped with a nozzle of 0,4 mm diameter at 200 °C. The filaments were deposited with a straight filling architecture in precise directions of 0° and 90° between two successive layers. The thickness of each layer was 0.15 mm with 80% filling for the cylinders and 100% filling for the rectangular pieces and the printing speed was 30 mm/s.

**Thermal properties.** PLA and PLA/CNCx composites were examined using differential scanning calorimeter (DSC) (Q20, TA instruments), equipped with an RCS90 cooling system (TA instruments). Samples were weighed in an aluminum TA pan and sealed. An empty sealed pan was used as reference. Samples were first cooled to 25 °C and then heated to 250 °C with a heating rate of 20 °C min⁻¹ under nitrogen atmosphere. The thermal properties of PLA and PLA/CNCx composites such as the glass transition temperature (Tg), the cold crystallization temperature (Tcc), the melting temperature (Tm) and the enthalpy of fusion (ΔHm) were evaluated by the heating scans. The degree of crystallinity (χc) was determined for both the cold crystallization peak and the melting peak using the following equations<sup>56</sup>:

$$ \chi_c = \frac{\Delta H_{cc}}{\Delta H_{m0}} * \frac{100}{W} \qquad (1) $$

where ΔH<sub>cc</sub> is the enthalpy of cold crystallization, W is the weight fraction of PLA in the sample and ΔH<sub>m0</sub> is the enthalpy of melting for 100% crystalline PLA material, which was taken as 93 J g<sup>-1</sup>. The thermogravimetric analysis (TGA) was performed using a TGA Q500 device (TA instruments). 10 mg of each sample were heated under nitrogen from room temperature to 800 °C at a heating rate of 10 °C min<sup>-1</sup>.
<figure><img src="image_10.png" /><figcaption>Figure 8. Fabrication steps of PLA/CNCx scaffold (extrusion of the composite filament and 3D printing).</figcaption></figure>

**Structural and morphological characterization.** The Fourier transform infrared (FTIR) spectroscopy was performed in an attenuated total reflectance mode (ATR-FTIR) on a NICOLET NEXUS model spectrometer. All spectra were recorded in a spectral range of 4000–650 cm⁻¹ with an accumulation of 32 scans at a resolution of 4 cm⁻¹.

The morphology of CNC, extruded PLA and composite filaments as well as 3D printed scaffolds were observed using HITACHI S4800 scanning electron microscopy. The samples were sputter-coated for 30 s with platinum using a Polaron SC7620 Mini Sputter Coater for SEM analysis. Image J software was employed to calculate the mean filament diameter by taking average at 20 points and the pore size by taking average of 8 pores, which is denoted as mean ± standard deviation. Topography of the scaffolds' surfaces was analyzed using 3D optical microscopy (Keyence) and a confocal chromatic roughness tester (STIL SA) equipped with a CHR1000 sensor.

The porosity of the scaffold was obtained by a liquid displacement method as reported in the literature⁶⁴. Briefly, the scaffolds were immersed in tubes containing a specific amount of ethanol (W₁) for 30 min. Then, the total weight of immersed scaffolds and ethanol was noted as W₂. After removing the scaffolds, the residual ethanol in the tubes was noted as W₃. The porosity of the scaffolds was measured according to the following equation:

$$ Porosity(\%) = 100 * \frac{W_1 - W_3}{W_2 - W_3} \qquad (2) $$

**Mechanical properties.** The mechanical properties of 3D printed PLA/CNC scaffolds were characterized using a tensile system (Zwick Roell), coupled with a 5KN load cell. The specimens were printed in the form of a dog bone (40 mm long, 4 mm wide and 1.5 mm thick). The specimens were then clamped between dedicated jaws and pulled at a speed of 0.05 mm s⁻¹ until they broke. The Zwick Roell software is able to calculate Young's modulus, maximum force at break and elongation at break. At least 5 specimens are printed under the same conditions in order to perform a statistical measurement.

**Contact angle and swelling.** The contact angle of the 3D printed scaffolds with water was examined to determine their hydrophilicity at room temperature using the sessile drop method. Two scaffolds were analyzed: PLA as the reference & PLA/CNC3 as the composite with the best mechanical performance. The drops were photographed using a Drop Shape analyzer—DSA25 equipped with a monochrome camera B-CAM-21-BW (CCCIR) and an R60 Led lamp (Conrad). Briefly, 6 µL of 18 M MilliQ water was dropped onto the surface of each scaffold and the contact angle at equilibrium (considered at 60 s) was recorded. One Touch Grabber and Image J software were used to calculate the contact angles. The swelling behavior of the scaffolds was determined by gravimetric method. Scaffolds of previously known weight (Wₐ) are immersed in 10 mL of 18 M MilliQ water, then incubated with continuous stirring at 37 °C for an equilibrium time assumed to be 7 days. Then the samples
<table><caption>**Table 4.** Description of the primers used for the gene expression analysis.</caption><thead><tr><td>Gene</td><td>Forward primer</td><td>Reverse primer</td></tr></thead><tbody><tr><td>ALP</td><td>5' GCTGGCAGTGGTCAGATGTT 3'</td><td>5' CTATCCTGGCTCCGTGCTC 3'</td></tr><tr><td>OCN</td><td>5' TTGGACACAAAGGCTGCAC 3'</td><td>5' CTCACACTCCTCGCCCTATT 3'</td></tr><tr><td>RunX2</td><td>5' CCTAAATCACTGAGGCGGTC 3'</td><td>5' CAGTAGATGGACCTCGGGAA 3'</td></tr></tbody></table>

were removed from the water, cleaned with filter paper to remove excess residual water adsorbed on the surface and weighed ($W_h$). Finally, the scaffold swelling index was determined as follows:

$$ Swelling index(\%) = 100 * \frac{W_h - W_d}{W_d} \qquad (3) $$

**Enzymatic degradation.** The enzymatic degradation of PLA and PLA/CNC3 was carried out using alcalase enzyme according to a protocol previously reported in the literature with modifications²⁹,⁷⁹,⁸⁰. Indeed, strips of PLA and PLA/CNC3 tissues (dimensions 5 [W] × 15 [l] × 0.13 [h]) of approximately 0.25 g were printed and immersed in 25 mL of TRIS buffer (pH 9.5, 60 °C) with an optimum concentration of 50% (w/w) (relative to the weight of the tissue) of alcalase required for enzymatic hydrolysis, 3 mM L-cysteine and 0.05% (w/w) of azide sodium (relative to the weight of the TRIS buffer). Degradation was assessed by determining weight loss at 4, 7, 14, 21 and 28 days. The samples were first dried at 105 °C for 90 min, cooled in a desiccator then weighed in a closed weighing bottle. The percentage of weight loss was calculated as follows:

$$ Weight\ loss(\%) = 100 * \frac{W_1 - W_2}{W_1} \qquad (4) $$

where $W_1$ and $W_2$ are the dry weight of the samples before and after biodegradation, respectively.

**Cell viability.** The osteoblastic cell line hFOB 1.19 (ATCC® CRL11372™, ATCC, USA)²⁷ was cultured using DMEM/F12 (Dulbecco's Modified Eagle Medium α) (Gibco 10565018) conditioned media supplemented with 10% (V/V) foetal bovine serum (FBS) (Eurobio CVFSVF00-01). Cells were cultured at 37 °C in 5% CO₂ in a 10 cm diameter petri dish and trypsinized using 0.05% Trypsin-EDTA (Gibco 25300-054). After sterilization with 70% (w/v) ethanol for 30 min and UV irradiation for 1 h, the filaments and printed scaffolds were dried at room temperature and then placed in contact with hFOB cells for 5 days.

Cell viability was analyzed using MTT assay carried out by incubating 100 µL of 0.5 mg/mL of MTT solution on the cells for 3 h. Purple coloured formazan crystals were dissolved using 100 µL of DMSO (BDH Prolab 23486.297) and the absorbance was recorded at 560 nm using Multiskan flat reader (thermos, USA).

**Adhesion assays.** For adhesion assays, hFOB 1.19 cells were cultured on the scaffolds for 4 days (seeding of $2 \times 10^4$ cells per well of 24 well plate containing the scaffold). Cells cultured on scaffolds were fixed with 2.5% glutaraldehyde solution for 2 h at room temperature. Scaffolds with cells were then incubated twice for 15 min at room temperature, using increasing ethanol concentrations (30%, 50%, 70%, 95% (v/v)), followed by absolute ethanol then hexamethyldisilazane overnight. For scaffold growth assays, luciferase expressing hFOB 1.19 cells were cultured on scaffolds for 6 days (seeding of $2 \times 10^4$ cells per well of 24 well plate containing the scaffold). Luciferase activity was quantified using a Luminometer after cell lysis. In order to investigate only the influence of cellulose on PLA biological properties, a 2D model was used.

**Gene expression analysis.** hFOB 1.19 cells were seeded at 500 000 cells/well in 6-well flat-bottom plates and incubated at 37 °C/5% CO₂ for 48 h in 2 ml of medium previously incubated for 24 h with the corresponding scaffolds at 0.2 g/ml. Total RNA was extracted using the EZNA Total RNA Kit I, following the manufacturer's instructions. The cDNA was synthesized from 1 µg of RNA by reverse transcriptase (Mix qScript, Invitrogen Life Technologies). PCR reactions were performed in a thermal cycler (LightCycler® 480) for 30 cycles. The primers used are provided in Table 4.

**Mineralization assays.** To assess bone bioactivity, the ability to form apatite on PLA and PLA/CNC3 scaffolds was studied by mineralized culture in vitro in simulated body fluid (SBF) according to the method described in the literature⁸¹. Briefly, the scaffolds were immersed in 10 mL of fresh SBF and then incubated in a thermoplastic incubator shaker at 37 °C. The SBF solution was changed every other day throughout the 4-weeks study period. After incubation, weekly samples were removed from the fluid and dried in a desiccator. Finally, scanning electron microscopy (SEM) and Energy Dispersive X-ray spectroscopy (Detector: Oxford Instruments X-Max AZTEC, UK; Microscopy: Zeiss EVO HD15, Germany) were carried out to investigate the degree of mineralization of the scaffolds.

**Statistical analysis.** For most of the quantitative characterization method, the ANOVA test was used to evaluate if the data from every group of samples showed a significant difference ($p < 0.05$ for significant). For the
degradation assay, an unpaired-t-test was performed to compare data from PLA to PLA/CNC3. For the biological tests, the Mann-Whitney test was performed for the statistical analysis.

## Data availability

All data generated or analyzed during this study are included in this published article (and its Supplementary Information files).

Received: 13 June 2022; Accepted: 2 December 2022

Published online: 08 December 2022

## References

1. Noyes, F. R. *Noyes' Knee Disorders: Surgery, Rehabilitation, Clinical Outcomes E-Book* (Elsevier Health Sciences, 2016).

2. Kanwar, S. & Vijayavenkataraman, S. Bioprinting design of 3D printed scaffolds for bone tissue engineering: A review. *Bioprinting* **24**, e00167 (2021).

3. Krzysztof, P. & Pokrowiecki, R. Porous titanium implants: A review. *Adv. Mater. Eng.* **20**, 1700648 (2018).

4. Aro, H. T. & Aho, A. J. Clinical use of bone allografts. *Ann. Med.* **25**, 403–412 (1993).

5. Iqbal, N. et al. Recent concepts in biodegradable polymers for tissue engineering paradigms: A critical review. *Int. Mater. Rev.* **64**, 91–126 (2019).

6. Bose, S., Roy, M. & Bandyopadhyay, A. Recent advances in bone tissue engineering scaffolds. *Trends Biotechnol.* **30**, 546–554 (2012).

7. Porter, J. R., Ruckh, T. T. & Popat, K. C. Bone tissue engineering: A review in bone biomimetics and drug delivery strategies. *Biotechnol. Prog.* **25**, 1539–1560 (2009).

8. Obrien, F. J. Biomaterials and scaffolds for tissue engineering. *Mater. Today* **14**, 88–95 (2011).

9. Ma, P. X. Scaffolds for tissue fabrication. *Mater. Today* **7**, 30–40 (2004).

10. Cheung, H.-Y., Lau, K.-T., Lu, T.-P. & Hui, D. A critical review on polymer-based bio-engineered materials for scaffold development. *Compos. Part B Eng.* **38**, 291–300 (2007).

11. Jonoobi, M., Harun, J., Mathew, A. P. & Oksman, K. Mechanical properties of cellulose nanofiber (CNF) reinforced polylactic acid (PLA) prepared by twin screw extrusion. *Compos. Sci. Technol.* **70**, 1742–1747 (2010).

12. Patel, D. K., Deb, S., Hexiu, J., Ganguly, K. & Lim, K. Bioactive electrospun nanocomposite scaffolds of poly(lactic acid)/ cellulose nanocrystals for bone tissue engineering. *Int. J. Biol. Macromol.* **162**, 1429–1441 (2020).

13. Elena, M. et al. Three-dimensional printed polycaprolactone-microcrystalline cellulose scaffolds. *J. Biomed. Mater. Res. Part B* **107B**, 521–528 (2019).

14. Cakmak, A. M. et al. 3D printed polycaprolactone/gelatin/bacterial cellulose/hydroxyapatite composite scaffold for bone tissue engineering. *Polymers (Basel)* **12**, 1962 (2020).

15. Meftahi, A. et al. Nanocelluloses as skin biocompatible materials for skincare, cosmetics, and healthcare: Formulations, regulations, and emerging applications. *Carbohydr. Polym.* **278**, 118956 (2022).

16. Bhattarai, N., Li, Z., Edmondson, D. & Zhang, M. Alginate-based nanofibrous scaffolds: Structural, mechanical, and biological properties. *Adv. Mater.* **18**, 1463–1467 (2006).

17. Feng, B., Tu, H., Yuan, H., Peng, H. & Zhang, Y. Acetic-acid-mediated miscibility toward electrospinning homogeneous composite nanofibers of GT/PCL. *Biomacromol* **13**, 3917–3925 (2012).

18. Bhattarai, N. et al. Natural-synthetic polyblend nanofibers for biomedical applications. *Adv. Mater.* **21**, 2792–2797 (2009).

19. Pon-on, W., Suntornsaratoon, P. & Charoenphandhu, N. Synthesis and investigations of mineral ions-loaded apatite from fish scale and PLA/chitosan composite for bone scaffolds. *Mater. Lett.* **221**, 143–146 (2018).

20. Edwards, A., Jarvis, D., Hopkins, T., Pixley, S. & Bhattarai, N. Poly(ε-caprolactone)/keratin-based composite nanofibers for biomedical applications. *J. Biomed. Mater. Res. Part B Appl. Biomater.* **103**, 21–30 (2015).

21. Belaid, H. et al. Boron nitride based nanobiocomposites: design by 3D printing for bone tissue engineering. *Acs Appl. Bio Mater.* **3**, 1865–1875 (2020).

22. Soundarya, S. P., Menon, A. H., Chandran, S. V. & Selvamurugan, N. Bone tissue engineering: Scaffold preparation using chitosan and other biomaterials with different design and fabrication techniques. *Int. J. Biol. Macromol.* **119**, 1228–1239 (2018).

23. Guo, L. et al. The role of natural polymers in bone tissue engineering. *J. Control. Release* **338**, 571–582 (2021).

24. Martínez-moreno, D., Jim, G. & Marchal, J. A. Pore geometry influences growth and cell adhesion of infrapatellar mesenchymal stem cells in biofabricated 3D thermoplastic scaffolds useful for cartilage tissue engineering. *Mater. Sci. Eng. C* **122**, 111933 (2021).

25. Bobbert, F. S. L. & Zadpoor, A. A. Effects of bone substitute architecture and surface properties on cell response, angiogenesis, and structure of new bone. *Mater. Chem. B* **5**, 6157–6414 (2017).

26. Mirkhalaf, M. et al. Redefining architectural effects in 3D printed scaffolds through rational design for optimal bone tissue regeneration. *Appl. Mater. Today* **25**, 101168 (2021).

27. Shirzad, M., Zolfagharian, A., Matbouei, A. & Bodaghi, M. Design, evaluation, and optimization of 3D printed truss scaffolds for bone tissue engineering. *J. Mech. Behav. Biomed. Mater.* **120**, 104594 (2021).

28. Radhakrishnan, S. et al. Fabrication of 3D printed antimicrobial polycaprolactone scaffolds for tissue engineering applications. *Mater. Sci. Eng. C* **118**, 111525 (2021).

29. Belaid, H. et al. Development of new biocompatible 3D printed graphene oxidebased scaffolds. *Mater. Sci. Eng. C* **110**, 110595 (2020).

30. Wang, W. et al. 3D printing of PLA/n-HA composite scaffolds with customized mechanical properties and biological functions for bone tissue engineering. *Compos. Part B* **224**, 109192 (2021).

31. Mondal, S., Phuoc, T., Pham, V. H. & Hoang, G. Hydroxyapatite nano bioceramics optimized 3D printed poly lactic acid scaffold for bone tissue engineering application. *Ceram. Int.* **46**, 3443–3455 (2020).

32. Zhang, B. et al. 3D printed bone tissue regenerative PLA/HA scaffolds with comprehensive performance optimizations. *Mater. Des.* **201**, 109490 (2021).

33. Du, J. et al. Comparative evaluation of chitosan, cellulose acetate, and polyethersulfone nanofiber scaffolds for neural differentiation. *Carbohydr. Polym.* **99**, 483–490 (2014).

34. Shi, Q. et al. The osteogenesis of bacterial cellulose scaffold loaded with bone morphogenetic protein-2. *Biomaterials* **33**, 6644–6649 (2012).

35. Joseph, B., Sagarika, V. K., Sabu, C., Kalarikkal, N. & Thomas, S. Cellulose nanocomposites: Fabrication and biomedical applications. *J. Bioresour. Bioprod.* **5**, 223–237 (2020).

36. Lin, N. & Dufresne, A. Nanocellulose in biomedicine: Current status and future prospect. *Eur. Polym. J.* **59**, 302–325 (2014).

37. Murizan, N. I. S., Mustafa, N. S., Ngadiman, N. H. A., Mohd Yusof, N. & Idris, A. Review on nanocrystalline cellulose in bone tissue engineering applications. *Polymers (Basel)* **12**, 2818 (2020).

38. Shuai, C., Yuan, X., Yang, W., Peng, S. & Qian, G. Synthesis of a mace-like cellulose nanocrystal @ Ag nanosystem via in-situ growth for antibacterial activities of poly-L-lactide scaffold. *Carbohydr. Polym.* **262**, 117937 (2021).
39. He, Y., Li, H., Fei, X. & Peng, L. Carboxymethyl cellulose/cellulose nanocrystals immobilized silver nanoparticles as an effective coating to improve barrier and antibacterial properties of paper for food packaging applications. *Carbohydr. Polym.* **252**, 117156 (2021).

40. Kumar, A. et al. Effect of crosslinking functionality on microstructure, mechanical properties, and in vitro cytocompatibility of cellulose nanocrystals reinforced poly(vinyl alcohol)/sodium alginate hybrid scaffolds. *Int. J. Biol. Macromol.* **95**, 962–973 (2017).

41. Rashtchian, M., Hivechi, A., Bahrami, S. H., Milan, P. B. & Simorgh, S. Fabricating alginate/poly(caprolactone) nano fibers with enhanced bio-mechanical properties via cellulose nanocrystal incorporation. *Carbohydr. Polym.* **233**, 115873 (2020).

42. Kumar, A., Madhusudana, K. & Soo, S. Development of sodium alginate-xanthan gum based nanocomposite scaffolds reinforced with cellulose nanocrystals and halloysite nanotubes. *Polym. Test.* **63**, 214–225 (2017).

43. Kumar, A., Rao, K. M. & Han, S. S. Mechanically viscoelastic nanoreinforced hybrid hydrogels composed of polyacrylamide, sodium carboxymethylcellulose, graphene oxide, and cellulose nanocrystals. *Carbohydr. Polym.* **193**, 228–238 (2018).

44. Sucinda, E. F. et al. International journal of biological macromolecules development and characterization of packaging film from Napier cellulose nanowhisker reinforced polylactic acid (PLA) bionanocomposites. *Int. J. Biol. Macromol.* **187**, 43–53 (2021).

45. Vorawongsagul, S., Pratumpong, P. & Pechyen, C. Preparation and foaming behavior of poly(lactic acid)/poly(butylene succinate)/ cellulose fiber composite for hot cups packaging application. *Food Packag. Shelf Life* **27**, 100608 (2021).

46. Ghaffari-bohlouli, P., Jafari, H., Khatibi, A., Bakhtiari, M. & Tavana, B. Osteogenesis enhancement using poly (L-lactide-co-D, L-lactide)/poly(vinyl alcohol) nano fibrous scaffolds reinforced by phospho-calcified cellulose nanowhiskers. *Int. J. Biol. Macromol.* **182**, 168–178 (2021).

47. Cui, L. et al. Rheology of PLA/regenerated cellulose nanocomposites prepared by the pickering emulsion process: Network formation and modeling. *Mater. Des.* **206**, 109774 (2021).

48. Chen, J., Zhang, T., Hua, W., Li, P. & Wang, X. 3D Porous poly(lactic acid)/regenerated cellulose composite scaffolds based on electrospun nanofibers for biomineralization. *Colloids Surf. A: Physicochem. Eng. Asp.* **585**, 124048 (2019).

49. Rajeshkumar, G. et al. Environment friendly, renewable and sustainable poly lactic acid (PLA) based natural fiber reinforced composites—A comprehensive review. *J. Clean. Prod.* **310**, 127483 (2021).

50. Mao, D., Li, Q., Bai, N., Dong, H. & Li, D. Porous stable poly(lactic acid)/ethyl cellulose/hydroxyapatite composite scaffolds prepared by a combined method for bone regeneration. *Carbohydr. Polym.* **180**, 104–111 (2018).

51. Yan, D. et al. Surface modified electrospun poly (lactic acid) fibrous scaffold with cellulose nano fibrils and Ag nanoparticles for ocular cell proliferation and antimicrobial application. *Mater. Sci. Eng.* **C** **111**, 110767 (2020).

52. Kumar, S. D., Venkateshwaran, K. & Aravindan, M. K. Fused deposition modelling of PLA reinforced with cellulose nano-crystals. *Mater. Today Proc.* **33**, 868–875 (2020).

53. Wang, Z. et al. Preparation of 3D printable micro/nanocellulose-polylactic acid (MNC/PLA) composite wire rods with high MNC constitution. *Ind. Crop. Prod.* **109**, 889–896 (2017).

54. Wang, Q., Ji, C., Sun, L., Sun, J. & Liu, J. Cellulose nanofibrils filled poly(lactic acid) biocomposite filament for FDM 3D printing. *Molecules* **25**, 2319 (2020).

55. Hafizi, W., Ishak, W., Rosli, N. A. & Ahmad, I. Influence of amorphous cellulose on mechanical, thermal, and hydrolytic degradation of poly(lactic acid) biocomposites. *Sci. Rep.* **10**, 11342 (2020).

56. Murphy, C. A. & Collins, M. N. Microcrystalline cellulose reinforced polylactic acid biocomposite filaments for 3D printing. *Polym. Compos.* **39**, 1311–1320 (2018).

57. Nagarajan, S. et al. Design of boron nitride/gelatin electrospun nanofibers for bone tissue engineering. *ACS Appl. Mater. Interfaces* **9**, 33695–33706 (2017).

58. Cataldi, A., Rigotti, D., Nguyen, V. D. H. & Pegoretti, A. Polyvinyl alcohol reinforced with crystalline nanocellulose for 3D printing application. *Mater. Today Commun.* **15**, 236–244 (2018).

59. Azarudeen, R. S. et al. 3D printable polycaprolactone-gelatin blends characterized for in vitro osteogenic potency Raja. *React. Funct. Polym.* **146**, 104445 (2020).

60. Alam, F., Shukla, V. R., Varadarajan, K. M. & Kumar, S. Microarchitected polylactic acid (PLA) nanocomposite scaffolds for biomedical applications. *J. Mech. Behav. Biomed. Mater.* **103**, 103576 (2020).

61. Stoof, D. & Pickering, K. Fused deposition modelling of natural fibre/polylactic acid composites. *J. Compos. Sci.* **1**, 8 (2017).

62. Mazzanti, V., Malagutti, L. & Mollica, F. FDM 3D printing of polymers containing natural fillers: A review of their mechanical properties. *Polymers (Basel)* **11**, 1094 (2019).

63. Xu, L., Zhao, J., Qian, S., Zhu, X. & Takahashi, J. Green-plasticized poly(lactic acid)/nanofibrillated cellulose biocomposites with high strength, good toughness and excellent heat resistance. *Compos. Sci. Technol.* **203**, 108613 (2021).

64. Tarrahi, R. et al. Development of a cellulose-based scaffold for sustained delivery of curcumin. *Int. J. Biol. Macromol.* **183**, 132–144 (2021).

65. Torgbo, S. & Sukyai, P. Fabrication of microporous bacterial cellulose embedded with magnetite and hydroxyapatite nanocomposite scaffold for bone tissue engineering. *Mater. Chem. Phys.* **237**, 121868 (2019).

66. Krieghoff, J. et al. Increased pore size of scaffolds improves coating efficiency with sulfated hyaluronan and mineralization capacity of osteoblasts. *Biomater. Res.* **23**, 23–26 (2019).

67. Malliappan, S. P., Yetisgin, A. A., Sahin, S. B., Demir, E. & Cetinel, S. Bone tissue engineering: Anionic polysaccharides as promising scaffolds Ponnurengam. *Carbohydr. Polym.* **283**, 119142 (2022).

68. Gorgieva, S., Girandon, L. & Kokol, V. Mineralization potential of cellulose-nano fi brils reinforced gelatine scaffolds for promoted calcium deposition by mesenchymal stem cells. *Mater. Sci. Eng.* **C** **73**, 478–489 (2017).

69. Szustakiewicz, K. The effect of selective mineralization of PLLA in simulated body fluid induced by ArF excimer laser irradiation: Tailored composites with potential in bone tissue engineering. *Compos. Sci. Technol.* **197**, 108279 (2020).

70. Dridi, A., Zlaoui, K. & Somrani, S. Mechanism of apatite formation on a poorly crystallized calcium phosphate in a simulated body fluid (SBF) at 37 °C. *J. Phys. Chem. Solids* **156**, 110122 (2021).

71. Zhang, Q. et al. Mechanical properties and biomineralization of multifunctional nanodiamond-PLLA composites for bone tissue engineering. *Biomaterials* **33**, 5067–5075 (2012).

72. Zou, L., Zhang, Y., Liu, X., Chen, J. & Zhang, Q. Biomimetic mineralization on natural and synthetic polymers to prepare hybrid scaffolds for bone tissue engineering. *Colloids Surf. B: Biointerfaces* **178**, 222–229 (2019).

73. Syahir, M. et al. Proceedings entrapment of collagen on polylactic acid 3D scaffold surface as a potential artificial bone replacement. *Mater. Today Proc.* **46**, 1668–1673 (2021).

74. Zhang, R. & Ma, P. X. Porous poly (L-lactic acid)/apatite composites created by biomimetic process. *J. Biomed. Mater. Resour.* **45**, 285–293 (1999).

75. Ngiam, M. et al. Fabrication of mineralized polymeric nanofibrous composites for bone graft materials. *Tissue Eng.* **15**, 535–546 (2009).

76. Fragal, E. H. et al. Hybrid materials for bone tissue engineering from biomimetic growth of hydroxiapatite on cellulose nanowhiskers. *Carbohydr. Polym.* **152**, 734–746 (2016).

77. Harris, S. A., Enger, R. J., Riggs, B. L. & Spelsberg, T. C. Development and characterization of a conditionally immortalized human fetal osteoblastic cell line. *J. Bone Miner. Res.* **10**, 178–186 (1995).

78. Zheng, D., Zhang, Y., Guo, Y. & Yue, J. Isolation and characterization of nanocellulose with a novel shape from walnut (*Juglans Regia* L.) shell agricultural waste. *Polymers (Basel)* **11**, 1130 (2019).
79. Lee, S. H. & Song, W. S. Modification of polylactic acid fabric by two lipolytic enzyme hydrolysis. *Text. Res. J.* **83**, 229–237 (2013).

80. Lee, S. H., Kim, I. Y. & Song, W. S. Biodegradation of polylactic acid (PLA) fibers using different enzymes. *Macromol. Res.* **22**, 657–663 (2014).

81. Kokubo, T. & Takadama, H. How useful is SBF in predicting in vivo bone bioactivity? *Biomaterials* **27**, 2907–2915 (2006).

## Acknowledgements

Hélène Garay is acknowledged for her support in topography analyses.

## Author contributions

C.S.: Conceptualization, Methodology, Validation, Investigation, Resources, Data Curation, Writing-Review & Editing, Visualization, Supervision, Project administration, Funding acquisition. M.B.: Methodology, Writing-Review & Editing. V.C.: Investigation, Resources, Data Curation, Writing-Review & Editing. D.B.: Funding acquisition, Writing-Review & Editing. N.M.: Investigation, Writing-Review & Editing. M.K.: Investigation, Validation. E.F.A.: Writing-Review & Editing. H.B.: Investigation, validation, Writing-Review & Editing. K.M.N.: Investigation, Validation, Writing—Original draft. J.E.H.: Investigation, Validation.

## Funding

This work was funded by the Ministry of Higher Education and Scientific Research (MESRS) of the Republic of Côte d'Ivoire, the French National Agency (ANR, JCJC program, MONOME ANR-20-CE08-0009) and IEM Axe Energy project PAT (Solid CeraLique).

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-022-25652-x.

Correspondence and requests for materials should be addressed to C.S.

Reprints and permissions information is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_11.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International

License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2022