Abstract
Bioprinting continues to expand as a preferred fabrication technology as new biomaterials emerge. These biomaterials must be biocompatible, in addition to exhibiting suitable printing characteristics. Printability is usually gauged by metrics such as extrudability, fiber morphology, and shape fidelity, and optimized via design of experiments or machine learning methods. However, these optimization techniques often overlook dynamic variables such as time-dependent behavior and batch-to-batch variability that can significantly affect the quality of printed structures. In this study, we 3D-printed colloidal hydrogels to investigate the effects of printing parameters and elapsed time on overall hydrogel print quality.
Impact Statement
We introduce colloidal hydrogels, which, unlike conventional systems, these materials incorporate a pressure-dependent variable that temporally shifts the optimal printing landscape. Understanding this evolving parameter space and establishing systematic models to predict and control it will be essential for harnessing the full potential of colloidal hydrogels in additive manufacturing.
Overview of Printability of Biomaterials
Additive manufacturing, or 3D printing, is a pivotal technology applied across domains from health care to electronics. 3D printing builds objects in a controlled, layer-by-layer fashion, minimizing waste and shortening fabrication time as opposed to subtractive and formative manufacturing. This technique enables the creation of intricate structures from polymers, ceramics, metals, and even biological materials. Common techniques include extrusion, inkjet, stereolithography, and laser-assisted printing, each offering distinct advantages in terms of resolutions and material compatibilities. 1 Among these, extrusion-based 3D printing is the most prevalent due to its adaptability to a wide range of materials, especially for tissue engineering applications.2,3 Ideal materials for extrusion-based 3D printing are those that flow under high shear or exhibit shear-thinning behavior and possess sufficient viscosity to retain the printed shape postprinting. 4 The convergence of 3D printing and tissue engineering has led to the emergence of bioprinting, a subfield focused on developing materials that meet these viscoelastic requirements while incorporating bioactive components essential for tissue engineering applications.
Hydrogels are hydrophilic polymer networks capable of absorbing large volumes of water that are well-suited for bioprinting because of their tissue-like swelling behavior. 4 While their biocompatibility has been well characterized in molded or injectable forms, their inherently low shape fidelity has hindered broader 3D printing applications. To overcome this, researchers have modified hydrogels to enable crosslinking, thereby bolstering shape stability postprinting. In our previous work, we demonstrated that naturally derived nanoparticle colloidal hydrogels are a promising bioprinting ink, highlighting the potential of colloidal hydrogel platforms to enhance tissue engineering applications.5,6 Unlike continuous polymer networks, colloidal gels form through the self-assembly of particles in suspension, driven by weak attractive forces, such as Van der Waals forces and electrostatic interactions. 7 The resulting network exhibits shear-thinning flow with self-healing, viscoelastic properties that facilitate extrusion printing and can serve as quantitative predictors of printability.8,9 With this introduction of colloidal gels as well as other novel hydrogel inks as bioprinting inks, there lacks the literature background necessary to determine optimal printing parameters.
Design of experiments (DOE), including full factorial designs, has been used to identify the primary factors affecting printability and to guide the development of related material systems. 10 However, DOE methods require large datasets, consuming significant time and resources. To overcome these challenges, data-driven machine learning models have been employed to optimize printability metrics on materials like gelatin methacryloyl (GelMA), 11 CELLINK start gel, 12 and granular hydrogels, 13 examining variables including composition, temperature, extrusion pressure, and print speed. However, an underlying assumption of these approaches is that a fixed parameter set will yield reproducible results. Here, a colloidal gel was selected as a novel ink platform to show that dynamic variables, such as elapsed time and batch-to-batch variability, also influence print fidelity and should be considered in other hydrogel printing optimization strategies. This work shows how print fidelity can be impacted by time, represented by scaffold count, and batch-to-batch variability, reflected in the range of printing parameters and material response across different ink batches.
Exploring the Printability of Colloidal Gels
Preparation of colloidal gel inks
Following previously published methods, gelatin type B (Sigma G9382, St. Louis, MO) was dissolved in Milli-Q water to create a 5 wt/v% solution at 50°C. 5 Acetone was added for the first desolvation, separating gelatin into lower and higher molecular weights. The precipitate was redissolved in Milli-Q water, freeze-dried, and subjected to a second desolvation with acetone, forming gelatin nanoparticles (GNPs). GNPs were crosslinked overnight with glutaraldehyde (Sigma G6257, St. Louis, MO), then quenched with glycine and washed by centrifugation. This reaction consumes primary amines on the gelatin backbone to form nanoparticles. GNPs used in this study were 250.4 ± 19.9 nm in diameter and exhibited a polydispersity index of 0.04 ± 0.02 as measured by a Zetasizer. For methacryloylation, GNPs were combined with carbonate-bicarbonate buffer and methacrylic anhydride (Sigma 276685, St. Louis, MO) at 50°C, pH 9, for 1 h. The resulting methacryloylated GNPs were dialyzed for 5 days, then freeze-dried after adjusting the pH to 7.4. Methacrylic anhydride reacts with the remaining primary amines on the GNPs, allowing for inter- and intraparticle photocrosslinking.
The dried methacryloylated GNPs and a photoinitiator, Irgacure 2959 (Advanced BioMatrix 5200) at 0.5 wt/v% in phosphate-buffered saline, were combined to create the colloidal gel inks. The combined nanoparticle content was 20 wt/v% and each ink was swollen at 4°C for 72 h with daily mechanical stirring and centrifugation. Inks were transferred to a 5 mL cartridge and were wrapped in aluminum foil to prevent undesirable photocrosslinking before 3D printing. However, the foil was removed during the printing process to expose the ink for visual inspection. No additional mixing was conducted, as air pockets would be introduced to the ink, biasing printing.
3D printing of GNP colloidal gel inks
The BioAssembly Bot 400 (Advanced Solutions, Louisville, KY) was used to 3D print 10 × 10 mm, two-layered gridded scaffolds at ambient temperature. The ink cartridge was loaded in the 3D Syringe pneumatic printhead and awaited printing within a print bay in the BioAssembly Bot. Initial printing pressure was determined based on the amount required to see a fiber extrude through a 20-gauge needle. Speed and pressure were changed and documented by the user if over- or under-extrusion was seen. After six scaffolds were printed within a petri dish (PD), the printhead was returned to the print bay and rested under ambient temperature and pressure for 15 min in between each dish. On average, each PD takes ## min to print. All PDs from one ink batch were printed on after the other on the same day. Images of the ink cartridge were taken before and after each rest period. After printing, scaffolds were imaged on an AmScope microscope (Irvine, CA) for characterization and were freeze-dried afterward.
Characterization of colloidal gel inks and 3D-printed scaffolds
Viscoelastic properties of colloidal gel inks
Viscoelastic properties, including shear-thinning and self-healing of the colloidal gel inks, were evaluated according to published methods.5,14,15 Briefly, these rheological tests were carried out at room temperature using a Discovery HR-2 rheometer (Waters, New Castle, DE) with a 40 mm stainless steel 2DEG SMART-SWA cone geometry. For assessing shear-thinning and self-healing, a flow ramp was performed with an increasing shear rate from 0.1 to 1000 1/s, and a three-step strain sweep was performed going from (I) constant low strain (1%) to (II) strain sweep increasing strain from 0.1% to 1000% and then (III) a return to constant low stain (1%). All tests were repeated on three technical replicates from one ink batch before and after printing all scaffolds for one ink batch. The average storage and loss moduli were recorded from phases I and III and used to compare the same ink before and after printing across four different ink batches (n = 4).
Fiber diameter of 3D-printed scaffolds
Images from the AmScope microscope were used to measure the average fiber diameter of the scaffolds using ImageJ. The first layer printing parameters produced vertical inner fibers, while the horizontal fibers represented the second layer’s printing parameters.
Representative Effects of Colloidal Gel Printability
Necessary changes to printing parameters throughout the printing process
The BioAssembly Bot 400 allows changes to the printing parameters, namely pressure and speed, while printing, which take effect on the next layer. This allows an experienced user to respond to changes in shape fidelity while printing by either adjusting the pressure, speed, or pausing the print to clear a clogged needle. These changes in shape fidelity and the response by the user were documented within annotated videos of various ink batches. Supplementary Video S1 shows the annotated video displaying these changes for one ink. These changes to the print parameters and shape fidelity from Supplementary Video S1 were summarized within Figure 1A. From this ink batch, it is noted that as time or the number of printed scaffolds increased, an increase in printing pressure is needed to maintain the same shape fidelity. Time is represented by the scaffold count based on how long it takes to print six scaffolds or a full PD. The time scale for printing a full PD is on average 25 min per PD, varying based on the print speed used. All scaffolds from one ink batch were printed within the same day, and each petri was printed back-to-back, with 15 min of rest in between. With this set-up, the scaffold count is continuous over time and can serve as a representation of time. Figure 1B further illustrates how scaffolds with different printing parameters can result in the same printing outcome and how the same printing parameters can result in different shape fidelity.

User adjustments to printing parameters (speed and pressure) during the 3D printing process to maintain the fiber diameter for one ink.
Material responses to the 3D printing process
This constant adjustment to the printing parameters has not been previously reported and raised concerns about changes to the material’s printability caused by the printing process. To assess the material’s response to the printing process, images of the ink cartridge were collected, and the ink’s viscoelastic properties were measured before and after printing. Cartridge images revealed a consistent decrease in piston height between printing each PD, indicating uniform extrusion as seen in Figure 2. Notably, the pneumatic print head compresses the ink during printing but may cause the piston to retract above the ink level once printing the PD is completed. The ink appeared visually homogenous throughout within images in Figure 2 and the video seen in Supplementary Video S2, with no evidence of phase separation between the nanoparticles and the suspension medium. Despite this, rheological measurements showed noticeable changes in material properties. As shown in Figure 3A, viscosity increased after printing across multiple ink batches. Additionally, comparison of the change of the storage and loss moduli during the printing revealed batch-dependent variation in differences seen in Figure 3B. While some batches showed minimal difference, others exhibited a significant increase in storage modulus postprinting, suggesting material properties are impacted during the printing process.

Illustrating the change in the colloidal gel ink within the print cartridge after 3D printing 6 scaffolds within 1 petri dish (PD) with a pneumatic extrusion-based system. Line-up of the ink cartridge, showing the change in ink and piston position after printing 30 scaffolds and waiting 15 min in between each PD.

Rheological evaluation of the viscoelastic properties of the GNP colloidal gel inks before and after printing.
Discussion/Perspective
Extrusion-based 3D printing enables the fabrication of complex structures using diverse materials and has become a powerful tool in tissue engineering. Colloidal gels, formed by the self-assembly of nanoparticles, exhibit promising viscoelastic properties, including shear-thinning and self-healing, making them strong candidates for 3D printing. Before printing, it is required to tune the printing parameters of any material to achieve optimal printability by using DOE or machine learning methods. This study explored how the printing process itself influences the properties of colloidal gels and, ultimately, the printability of the gels that is not considered during the initial tuning studies. Maintaining shape fidelity requires continual adjustments to the pressure and speed during the printing process (Fig. 1A) for several separate batches of colloidal gel inks (Supplementary Figs. S3, S4, S5 and S6). In addition, there is notable variation in the range of pressures needed to extrude a uniform fiber between ink batches. For instance, two ink batches printed well at a pressure of 4 psi, while another achieved similar fidelity at only 1.5 psi. Such variability between inks limits the direct transferability of optimized parameters between batches and undermines the assumption commonly made in DOE and machine learning approaches that a single parameter set yields a consistent print outcome. This inconsistency is often simplified as a “printing range” in literature, which highlights a broader challenge across printing colloidal gel platforms.
To investigate the source of this variability, we measured the change in the viscoelastic properties of the colloidal gel inks before and after printing. Results showed a general increase in storage modulus postprinting, indicating a rise in material stiffness. This shift likely explains the need for increased pressure during printing. We hypothesize that this change arises from a pressure-induced phase separation within the GNP colloidal gel ink, leading to a GNP concentration gradient in the ink cartridge. Phase separation during printing may cause a nonuniform distribution of water leaving the syringe, resulting in progressively less water remaining in the syringe. This process creates a nanoparticle gradient, with a higher particle concentration accumulating toward the end of printing potentially causing aggregation and precipitation of GNPs. Such precipitation is expected during 3D printing, given that hydrated gelatin has a higher density (1.2 g/mL) than water, which explains why water exits the syringe more rapidly. The resulting increase in GNP concentration is expected to raise the viscosity and storage modulus of the ink, as previous work has shown both properties to be dependent on GNP concentration. 5 This subtle change in material properties would be visually undetectable to users (Fig. 2) but results in fluctuations in fiber diameter, print uniformity, and storage modulus.
These findings emphasize the importance of documenting real-time parameter adjustments and printing conditions in bioprinting studies. Variations in print performance should not be dismissed as noise but recognized as fundamental behavior of the material and printing system. This work aims to encourage the gathering of diverse data, including dynamic factors and captured data, to further improve the predictive strength of optimization processes. Capturing supporting data such as photographs and videos during printing (Supplementary Videos S1 and S2) can enrich datasets, enabling more accurate models that link material properties and printing parameters with outcome print quality. Thus, a deeper understanding of colloidal biomaterial behavior in extrusion-based 3D printing is essential for developing adaptive, intelligent biofabrication strategies.
Footnotes
Acknowledgments
The authors gratefully acknowledge Dr. Lydia E. Kavraki from the Department of Computer Science at Rice University for their insightful contributions to the conceptual development of this work and for providing guidance in the preparation of this article.
Author Disclosure Statement
No competing financial interests exist.
Funding Information
The authors acknowledge support from a National Science Foundation Graduate Research Fellowship (MRP), the National Institutes of Health (P41 EB023833, AGM), and the Rice Academy of Postdoctoral Fellows (VKK) and the National Institutes of Health Interdisciplinary Translational Postdoctoral Program in Cancer Nanotechnology (T32CA196561) (VKK). Additional support was provided by the Biomaterials Lab at Rice University.
Supplemental Material
Supplemental Material
Supplemental Material
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
