1. Introduction:
Bone is a versatile tissue that has a significant capacity for regeneration. However, the successful regeneration of complex bone defects is currently one of the biggest challenges faced by reconstructive surgeons [
1]. Smaller bone defects can usually heal without intervention, as is the case with many types of fractures and small cranial bone defects [
2]. Large or complex defects, termed critical-sized defects, exceed the natural capacity for regeneration, requiring further intervention [
2]. Currently, the most common methods used for repairing large bone defects are autografting (transplanting bone from another site on a single patient) and allografting (transplanting cadaver bone) [
2]. These structures are often hand-carved to fit the size of the defect [
2]. Despite being considered the current gold standard, autografting and allografting are associated with significant drawbacks. Autografting can lead to donor site morbidity and requires more complex surgery. Allografting may be limited by the supply of donors, the potential for immunogenicity, and the risk of disease transmission [
1]. Although conflicts of interest should be noted, studies on cell-infused cellular bone matrices have shown high success rates and low complication rates in spinal fusion [
3].
In the past decade, tissue engineering of bone through 3-Dimensional Bioprinting (3DBP) has emerged as a possible solution for the challenges associated with bone autografting and allografting [
4]. 3DBP is the fabrication of a tissue-like construct through the deposition of living cells often suspended in a hydrogel. When the cells and hydrogel are combined, they make up a liquid “bioink” that can be extruded from a bioprinter in layers to produce a cellular scaffold [
5,
6]. This scaffold can be subjected to covalent photocrosslinking during or after printing, producing a hardened construct that resembles the bone microenvironment. A significant advantage of 3DBP over allografting and autografting is that patient-derived autologous mesenchymal stromal cells (hMSCs) could be used in the construct, reducing the potential for graft immunogenicity and donor site morbidity. Additional advantages of 3DBP over allografting are decreased risk of infection, the ability to easily create large complex shapes, and the use of organic and biological additives to aid osteogenesis [
5,
7].
A significant challenge in 3DBP bone-like constructs is formulating a bioink that closely mimics the bone microenvironment. Natural bone is typically 20-40% collagen, 60-80% mineral, and 10-20% water [
8]. The young’s modulus of cortical and trabecular bone is ~19.9 GPa and ~18.0 GPa, respectively [
9]. While the modulus of pre-calcified bone matrix housing bone marrow is around 64 kPa [
10]. Bone bioinks currently used for 3DBP generally consist of three components: the base material, living cells, and various natural or synthetic additives. The base materials (most commonly alginate, gelatin methacrylate, or hyaluronic acid) in the bioink can be cross-linked after printing to increase the mechanical stiffness of the construct significantly [
7,
11].
Gelatin (a hydrolyzed product of collagen) with methacryloyl side groups (GelMA) is a popular base material in bone 3DBP because of its ability to be cross-linked in the presence of a photoinitiator [
12,
13,
14]. Low concentrations of pure GelMA (5% w/v) can be extruded more easily and generate higher cell viability [
15]. However, low concentrations of GelMA lack the stiffness needed to produce a construct that resembles the rigidity and load-bearing capacity of bone. Higher concentrations of GelMA (15% w/v) require higher pressure to be extruded and may decrease cell viability, although this can be overcome by warming the material to a liquid before extrusion [
15]. To maximize printability and cell viability, most groups have used GelMA concentrations between 5% and 15% w/v for 3DBP [
7]. In addition, low concentrations of pure gelatin can be added to GelMA hydrogels to increase the viscosity and printability of the material, without significantly affecting the material characteristics [
16].
Natural and synthetic additives such as hydroxyapatite, bioactive glass, and graphene oxide have a significant impact on the structural stability of a bioink [
17]. Additives can also increase cell viability within the hydrogel in low concentrations and promote differentiation into target tissues like bone or cartilage [
18]. Hydroxyapatite (HA) is the primary inorganic material found in bone tissue which makes it a desirable additive for bone tissue engineering [
17]. When used as an additive with pure gelatin hydrogels, HA has been shown to improve the mechanical stiffness of engineered bone tissue by up to 70% (42 kPa to 72 kPa) [
19]. However, higher concentrations of HA may decrease cell viability due to its high calcium content [
20].
Some of the most important characteristics of a bioprinted structure are printability, print fidelity, mechanical strength, and cell viability [
21]. Printability, also called extrudability in extrusion bioprinting, is defined as the printing parameters necessary to print a straight line of defined width. Printability is influenced by pressure, temperature, speed, acceleration, and the diameter of the needle used for printing [
21]. Determining printability is an important step in evaluating the feasibility of printing a bioink. Print fidelity is the degree to which the bioprinted structure holds its shape and it is related to the viscosity of the bioink [
22]. 3D bioprinted structures often require a high print fidelity and it is crucial to balance print fidelity with printability to enable custom, defect-specific implants based on high resolution clinical imaging [
7]. The mechanical stiffness of the printed product is particularly important as bioprinted bone structures should have significant load-bearing capabilities to maximize bone healing when implanted [
23]. Mechanical properties can be determined by stress-testing the material. In this study we will calculate the equilibrium modulus to determine the mechanical properties of the materials [
24]. Lastly, the high pressures and cross-linking procedures required for 3DBP bone tissues with high mechanical strength can have a deleterious effect on cell viability [
15]. Evaluating the viability of human cells in each printable bioink after the printing process is an important step in determining their suitability as bone replacement implants.
A previous study used GelMA and HA for 3DBP and determined that the addition of HA in low concentrations (≤2%) to a 10% GelMA hydrogel maintained equal cell viability and decreased hydrogel swelling in comparison to constructs without HA [
25]. The ideal combination of these materials for optimizing printability, print fidelity, mechanical strength, and cell viability remains largely unknown.
Testing bioprinting with a wide range of material concentrations with more than one material is challenging. Design of Experiments software (DoE) can aid in the experimental planning and analysis [
13,
14,
26]. The software allows researchers to determine optimal values for further testing from fewer experiments based on desired outcomes.
For example, if the software is given maximum and minimum values of two material concentrations within a hydrogel, it will determine a subset of combinations within that range that will produce the greatest coverage of the design space with the fewest runs in an experiment. After the experiment, the results can be analyzed by the software, which can then predict the ideal concentrations of both materials in a hydrogel for a desired outcome [
13,
14,
26].
The objective of this study was therefore to determine the ideal GelMA/HA bioink combination with added gelatin for extrusion-based 3DBP of bone-like constructs laden with hMSCs. To accomplish this, we experimentally evaluated the effects of various combinations of GelMA and HA with 2% Gelatin on printability, print fidelity, mechanical properties, and cell viability.
We hypothesized there would be a direct relationship between GelMA/HA concentrations and equilibrium modulus, but an inverse relationship between GelMA/HA concentrations and printability, print fidelity, and cell viability.
2. Materials and Methods
2.1. Cell Culture
Human bone-marrow-derived mesenchymal stromal cells (hMSCs) were thawed at passage 2 from a frozen stock of hMSCs isolated under Institutional Review Board-approved non-human research protocols (University of Central Florida College of Medicine STUDY00001124) and cultured until passage 5, at which point the experiments were performed. The cells were cultured in growth medium containing Dulbecco’s Modified Eagle Medium (DMEM; Gibco Inc.) supplemented with 10% fetal bovine serum (FBS, Atlanta Biologicals Inc.), and 1% penicillin/streptomycin (Gibco). The cells were incubated in a humidified environment at 37°C with 5% CO2. At ~90% confluence, the cells were trypsinized (0.25% trypsin/EDTA, Gibco), neutralized with growth media, centrifuged (500 RCF, 5 min), and resuspended in growth media. Live cells were counted (Cytosmart) using trypan blue exclusion. Cells were combined with the bioinks at a final concentration of ~1 x 106 cells/mL.
2.2. Design of Experiments Parameters
Experimental design software (DesignExpert, StatEase Inc.) was utilized to determine the hydrogel formulations for testing based on GelMA concentrations between 5% and 15% w/v and HA concentrations between 1% and 33%. The software produced 17 experimental runs with various GelMA/HA bioink formulations at 2% Gelatin, some of which are repeated to test experimental variability (
Table 1).
The 17 hydrogels were first assessed based on their pipettability and equilibrium modulus of a cast cylinder. The bioinks with the ideal printability and mechanical properties were later 3D bioprinted and evaluated for their print fidelity and cell viability.
2.3. Bioink formulation
The 17 experimental bioinks (
Table 1) were prepared using GelMA at a methacrylation level of 45-55% (Rousselot Inc.) in PBS containing 0.05% lithium phenyl-2,4,6-trimethyl-benzoyl phosphinate (LAP) as the photoinitiator (Cellink Inc.) with 2% 300 bloom Gelatin (Electron Microscopy Sciences Inc.). The solution was warmed at 37°C on a hot plate and vortexed at 1500 revolutions per minute until dissolved. Various concentrations of HA powder (Sigma-Aldrich) were added to the solution and stirred until suspended. The hydrogel was combined with the cultured hMSCs at a concentration of ~1x10
6 cells/mL [
27].
2.4. Printability
For the purposes of this study, printability was defined as the ability of a material to be pipetted and extruded as a gel at a temperature under 37 ℃. A 0.2 mL sample of each of the 17 composite hydrogels (no cells) was warmed to 37 ℃ and pipetted. Bioinks that could not be pipetted or extruded were deemed unprintable. A composite score was developed for each bioink from the results, where printability is a dichotomous variable (1 = printable, 0 = not printable).
2.5. Mechanical testing
Mechanical testing was performed using a TA.XTPlusC mechanical testing device (Stable Micro Systems Inc.). Each of the 17 acellular composite hydrogels were heated to 37°C and pipetted into circular casts made using a biopsy punch in a silicone mat giving discs 5 mm in diameter by 0.5 mm in height. The discs and casts were held with a microscope slide and coverslip. They were then photo-crosslinked using a 300mW laser with visible light at 405 nm (Luck Laser Inc.). The laser was positioned 3.5cm above the hydrogels which were treated on both sides for 60 seconds.
A Texture-Analyzer 12.7 mm diameter cylindrical probe (TA-10) with a 5kg load cell was used for compression. Uniaxial compression was initiated at a speed of 0.1 mm/s and trigger force of 0.1g. Cylinders were compressed to 5, 10, 15, 20 and 25% strain and allowed to equilibrate for 5, 10, 15, 30 and 45 mins respectively. The materials were tested at room temperature in PBS. The output was plotted using the Exponent Connect software (Stable Micro Systems Inc.) to produce the force/time curve. The slope of the curve fitted to those stress/strain points was taken as the equilibrium modulus.
2.6. Bioprinting
Bioprinting was performed with a BioAssemblyBot 3D bioprinter (Advanced Solutions Inc.) using pneumatic extrusion bioprinting (
Figure 1). Each bioink was loaded into a disposable, UV-blocking amber cartridge (Nordson Inc.) with 20-gauge 0.5 inch straight needle dispensing tips (Nordson Inc.). Tissue Structure Information Modeling (TSIM, Advanced Solutions Inc.) software was used for 3D modeling. After printing, constructs were photo-crosslinked using a 200mW laser and visible light at 405nm positioned 3.5cm above the construct (Luck Laser Inc.).
2.7. Development of a Laser Adapter on Hot Tool
The hot tool available when we started these experiments did not have the ability to photocrosslink, and the separate crosslinking tool was a 365nm light. We therefore sought to develop an on-tool laser to enable photocrosslinking while printing. A 405nm laser (200mW, Luck Laser) was chosen as this wavelength does less damage to DNA than the 365nm wavelength [
28] and can photocrosslink using the LAP photoinitiator [
13,
14]. A holder arm was 3D printed in PLA to attach the laser, battery packs and Raspberry Pi Pico (
Figure 1D) [
27]. This enables the printing of thicker constructs as crosslinking may not be possible after the print is completed.
2.8. Print Settings
The ideal printing pressures and temperatures for the composite hydrogels were determined using a systematic approach where the pressure was gradually increased in increments of 5 psi, while the temperature was increased from 30°C to 37°C at each interval until a continuous filament was achieved.
2.9. Print Fidelity
Print fidelity is defined as the degree to which a bioprinted structure holds its shape [
22]. The bioinks were bioprinted into a square zig zag pattern (
Figure 2) [
13]. The constructs were crosslinked via the above-mentioned procedure and imaged (VHX-7000, Keyence Inc.).
Photographs were taken and overlayed onto the TSIM model. The area was measured using an area calculator [
29]. These measurements were used to calculate the print fidelity according to the equation below, where
Ath is the theoretical pore area and
Ae is the experimental pore area. A value of 100% is perfect print fidelity, a value below 100% indicates underprinting, and a value above 100% represents overprinting.
2.10. Viability
Viability testing was performed using a LIVE/DEAD Viability/Cytotoxicity Kit (ThermoFisher Inc.). Bioinks were both pipetted and printed onto a well plate in triplicate. The bioinks were crosslinked and incubated in growth medium for 2 days and 7 days at 37 °C and 5% CO2. After each time point, growth medium was removed, and the constructs were washed in PBS. A solution containing 2µM Calcein AM (Invitrogen) and 4 µM Ethidium homodimer-1 (Invitrogen) in PBS was added to the wells and incubated for 25 minutes at 37 °C and 5% CO2. The dye was then removed, and PBS was added to the wells. The constructs were then imaged and overlayed (Keyence BZ-X810 Microscope). The images were processed and analyzed using the ImageJ/Fiji software with the Stardist plugin.
2.11. Statistical analysis
The design space covered by the DoE combinations was analyzed (StatEase). The z-axis is the standard error of design, and the x/y-axes are GelMA/HA concentrations as pictured (
Figure 3). Following testing, ANOVA was performed using the DoE software to demonstrate the effects of GelMA and HA concentrations on the equilibrium modulus and printability and the overall fit of the data to the model. Based on the compression and printability results, we used the DoE software to determine the ideal concentrations of GelMA and HA to optimize mechanical strength. The constraints used were that GelMA had to lie between 5 and 15%, HA between 1 and 33%, pipetability was targeted to be = 0.999 (as close to 1 as possible) with a lower limit of 0.8, and the goal to maximize the equilibrium modulus with a lower limit of 160 kPa and upper limit of 354 kPa.
Statistical analysis and data illustration for cell viability were performed using GraphPad Prism software. A two-way ANOVA test followed by Sidak’s correction for multiple comparisons in the post-hoc test was used to determine statistical significance. Where applicable, data is shown as mean +/- the standard deviation.
3. Discussion
An initial screen of hydrogel mixtures was used to make a DoE prediction for the stiffest printable gel. We successfully achieved a well correlated prediction with Gel 1, composed of 12.3% GelMA, 15.7% hydroxyapatite, and 2% gelatin giving an equilibrium modulus ~170 kPa vs. predicted 170 kPa. Gel 2, composed of 12.5% GelMA, 9% hydroxyapatite, and 2% gelatin had a slightly lower equilibrium modulus (~146 kPa) but yielded significantly higher cell viability (~82%) vs. Gel 1 (~60%). Gel 3, composed of 10% GelMA, 2% hydroxyapatite, and 5% gelatin had the highest print fidelity. Thus, we have achieved an increase in the mineral content of bioink to 9% HA vs. 2% in Gel 3 with good viability. However, native bone has a mineral content of ~60-80% [
8] so we fall significantly short of that target. It is expected that this scaffold would increase in mineralization with time either in culture or in vivo [
30].
Cell viability decreases at temperatures above 37 ℃, thus a cutoff of 37 ℃ was chosen for printability. If a gel is not printable at 37 ℃ or below, it cannot sustain cell viability. It was also found that at temperatures above 37℃, the hydrogels slowly weaken (i.e. lose compressive moduli, data not shown), capping the max temperature at 37℃.
There are multiple possible explanations for the drop in printability at high HA concentrations. For example, high amounts of HA may not have remained suspended within the hydrogels, which then formed a plug at the bottom of the syringe preventing extrusion.
Bone has the greatest stiffness of any tissue in the human body and can sustain immense amounts of pressure. As previously mentioned, the young’s modulus of cortical and trabecular bone is ~19.9 GPa and ~18.0 GPa, respectively [
9]. While this number may be impossible to achieve using hydrogels, a 3D bioprinted bone construct with an equilibrium modulus >64 kPa was achieved, making it similar to the hardness of the pre-calcified bone matrix [
10]. In addition, the hMSCs may increase the hardness of the material over time due to the release of osteogenic molecules that increase mechanical stiffness [
29].
The increasing opacity of high % HA may have also affected the mechanical stiffness by preventing the laser from fully penetrating the material, creating non-uniform crosslinking. This was addressed by simultaneously printing and crosslinking. One drawback to this is that each part of the gel is only crosslinked for a moment as the laser passes through.
The two print settings that were quantified were temperature and pressure (
Figure 6B). Both had a significant effect on the quality of the extruded filament. As expected, hydrogels containing higher GelMA and HA concentrations required a higher temperature and pressure to achieve a continuous filament. One of the challenges when heating the gel within the printer was a potential lack of uniformity of temperature within the tube itself. This resulted in regions of varying temperature and while printing resulting in some patterns of over-extrusion followed by under-extrusion. The temperature near the tip of the heated syringe was also likely lower as it was outside of the heating apparatus. As a result, higher pressures were needed in order to start extrusion and then the pressure would be lowered to maintain the print. Another factor that was less optimized was the print speed and bed temperature, which can have a significant effect on the fidelity.
Print fidelity is crucial for assessing the ability of a bioink to be engineered for patient-specific implants and depends on many factors such as gelation time, temperature, print speed and pressure. Bioinks with high print fidelity will produce a value close to 100%, meaning they hold their shape to a significant degree. Bioinks with greater than 100% print fidelity will be overprinted and collapse on the central pore. Bioinks with lower than 100% print fidelity will be underprinted and produce a central pore area that is larger than that of the model.
With the print settings shown in
Figure 6B, Gel 3 produced the most optimal print fidelity. Although all three gels produced constructs that qualitatively modeled the computer aided design, all were slightly overprinted. While some of the print fidelity results can be attributed to the characteristics of the gel, it is difficult to determine the degree to which the bioink itself influences print fidelity. The material concentrations influence the time required for the material to become a solid after printing. However, this effect could be controlled by performing photocrosslinking during the bioprinting process. Gelation time, temperature, print speed and pressure can all heavily influence print fidelity. It is likely that with more precise control over temperature and pressure, and with immediate photocrosslinking, all 3 gels could produce a more optimal print fidelity. Therefore, although Gel 3 produced the print fidelity closest to 100%, this may not mean that it is significantly better at modeling the computer aided design than Gel 1 or Gel 2. Despite these limitations, we determined that Gel 3 performed the best with the print settings shown in
Figure 6B.
All three materials sustained adequate cell viability, indicating that hMSCs can survive within the hydrogels through the photocrosslinking process. The results also showed that all three hydrogels can sustain hMSC viability for at least 7 days in culture. Increasing concentrations of HA may have affected cell viability in gel 1 which had an HA concentration of 15.7%. This is significantly higher than the highest HA concentration previously tested in the literature, which was 2% HA [
16]. We predicted that a printed gel may have lower cell viability as a result of high extrusion pressures. This however was not the case as the cell viability in printed gel 1 was not significantly different from pipetted gel 1. Cell viability did not change with intermittent pressures of up to 40psi. Of the three gels tested, Gel 2 seemed to exhibit a trend toward consistently higher viability, but this difference was not statistically significant when compared to gel 3. These results indicate that although all three hydrogels are suitable for sustaining hMSCs for an extended period, gel 2 may be the most optimal given the combination of stiffness and viability for 3D bioprinted bone implants.
The cell viability in printed gel 1 did not change significantly when compared to pipetted gel 1. Analysis of viability in printed gel 2 and gel 3 was not possible because the data was not quantifiable (
Supplemental Figure 1). Given that the higher pressure necessary to extrude gel 1 (20 psi vs 14-18 psi for gels 2 and 3) did not seem to significantly affect viability, cell viability likely would not have been different when comparing printed and pipetted gel 2 and 3.