Development of the thermosensitive EXtrusion Patterned Embedded construCT (EXPECT)
We approached the design of the EXPECT hydrogel with two goals in mind. First, we decided to place the cells locally together inside an oriented embedded channel within a 3D hydrogel network, to facilitate cellular communication. Embedded 3D printing offered one viable strategy for achieving this. Previously, groups have reported the writing of sacrificial inks within granular gel media exhibiting yield stress and self-healing properties [37–39]. After 3D printing, the sacrificial ink is removed by perfusion and the resulting embedded hollow channels can be used for a variety of applications, such as organoid fusion or providing enhanced nutrient diffusion [40, 41]. Carbopol® 940 (CP) is one of the most extensively studied granular gel media [42–47] known for its good printing performance and cellular compatibility. The microparticles of CP consist of polyacrylic acid with a high molecular weight, featuring methyl and methylene substituent groups arranged in a syndiotactic manner along the polymer chain [48]. Entangled and close-packed polymer chains of CP form sponge-like microgels, in which chains extend from the core and can interact with other particles. The extended chains strengthen inter-particle interactions and forming a fully entangled microparticle network [49]. In our study, we selected CP as one of the components of EXPECT, which would allow us to align cells within an embedded channel of a ring geometry. Because enclosed rings are a geometry that offers the minimum perimeter for a given enclosed area, we projected that the pattern would enable us to minimize the spatial distance between cells. At the same time, the ring could be extruded in a macroscopic size (10 mm diameter), so we anticipated that migratory movements along the ring would constitute straight ahead, directional movement on the cellular scale.
Our second goal for the design of EXPECT was to provide a means of sustaining cell migratory activities in a forward and backward motion, promoting unidirectional organizational behavior. Although embedded bioprinting has marked a notable progression in biomanufacturing, granular media exhibit a high degree of isotropy, lacking any strong asymmetric cues capable of guiding embedded cells in a specific direction throughout extended in vitro culture periods [35, 50, 51]. To address the limitations of using a hydrogel like CP alone, we introduced an additional component, a copolymer of polyN-isopropylacrylamide-graft-chondroitin sulfate (pNIPAAm-CS), into the hydrogel. This copolymer, developed and characterized previously [28, 29], demonstrates LCST behavior in aqueous solution at approximately 32 ºC. The LCST behavior of pNIPAAm-based systems is well-documented to involve a shift in molecular conformation from a solvated, rod-like state to a hydrophobic, globular state [36, 52, 53]. Based on previous investigations, it is known that culturing cells on 2D substrates containing pNIPAAm above its LCST promotes cell spreading and attachment, while the solvated molecular state below the LCST triggers the immediate release of attached cells and proteins [30–33]. Capitalizing on this knowledge, we explored the possibility of printing cells in an embedded ring channel that incorporates a pNIPAAm-based thermosensitive surface. We hypothesized that cooling below the LCST could effectively reduce the random cell migration out of these channels by inducing periodic detachment from the matrix. Figure 1 illustrates a schematic representation of the approach employed in this study.
Establishing the composition of EXPECT
After determining the main components of the embedding medium, the next step was to design a system that achieved a balance of desired characteristics from each phase. To begin, pNIPAAm-CS copolymers were synthesized in house by free-radical polymerization [28, 29, 54]. After synthesis and purification, we blended aqueous solutions of 1%, 3%, and 5% (w/v) pNIPAAm-CS with 0.8% (w/v) CP. In order to compare the results, we used 0.8% (w/v) CP alone, which is in a well-established composition range for embedded bioprinting [45].
We conducted amplitude sweep tests at 25°C (Supplementary Fig. 1a) and found that the formulations exhibited solid-like behavior within the linear viscoelastic region, below 1% strain. Additionally, the formulations could fluidize when subjected to a specific yield stress unique to each formulation, with increases in pNIPAAm-CS content leading to a reduced flow stress (Supplementary Fig. 1b). This suggested that the presence of pNIPAAm-CS chains reduced cohesive CP-CP inter-particle interactions. Further supporting this concept, print fidelity characterization of the patternable hydrogel was carried out by means of printing a free-floating horizontal grid structure with a low viscosity of 6% (w/v) gelatin bioink supplemented with 0.1% (w/v) Coomassie Blue dye. Printing in a high concentration of 5% pNIPAAm-CS + 0.8% CP resulted in deformation of the printed gid compared to the other embedding formulations (Supplementary Fig. 1c).
Despite possible losses in print fidelity at higher concentrations of pNIPAAm-CS, incorporating this temperature-responsive component into the CP microparticles rendered the hydrogel sensitive to temperature changes between 25°C to 37°C°C, exemplifying an increase in the storage modulus (G') between 66 and 730% (Supplementary Fig. 1d). The magnitude of this increase was proportional to the concentration of pNIPAAm-CS in the hydrogel formulation (Fig. 1e). To strike a balance between printability and thermal response, we decided to proceed with intermediate concentrations of pNIPAAm-CS in the thermosensitive embedding medium (3% pNIPAAm-CS + 0.8% CP). Additionally, we discovered that by blending a small concentration of gelatin (1% (w/v)) with this composition, we could recover key rheological properties at 25°C compared to 0.8% CP alone, summarized in the following sections.
Rheological properties of EXPECT
The shear recovery behavior of EXPECT (3% pNIPAAm-CS + 0.8% CP + 1% gelatin) was comparable to that of 0.8% CP alone. In an oscillatory step strain test (Fig. 2a), we measured the storage and loss moduli (G' and G'', respectively) while alternating the strain between 1% and 250%. Over three cycles, EXPECT maintained self-healing properties, with a mean recovery of 94.6 ± 6.5%, compared to the CP control exhibiting a mean recovery of 96.9 ± 4.0% (p = 0.4). These results indicate that EXPECT can effectively reform its network after nozzle translation in an extrusion-based printing process.
In the amplitude sweep test with oscillatory strain at 25°C (Fig. 2b), we determined that the flow stress of EXPECT, calculated at the point where the curves crossed over (G' = G''), was 366 ± 43 Pa, slightly higher than the CP control (277 ± 4 Pa) (Fig. 2c). This suggests that EXPECT can provide equivalent or improved spatial stability to a bioink extruded within its structure. Print fidelity studies with 6% (w/v) gelatin bioionk containing 0.1% (w/v) Coomassie Blue dye (Fig. 2d) further confirmed this, where multiple measures of printing accuracy, including printability (PR) value, pattern angle, and pattern width, were similar to the same bioink extruded within CP.
When the temperature increased from room temperature (25°C) to physiological temperature (37°C), EXPECT exhibited temperature-responsiveness (Fig. 2e), with G' increasing by 195 ± 65% between these temperatures. In comparison, the control CP showed a slight increase in G' of 5 ± 1% (Fig. 2f, p = 0.007). Macroscopic images of EXPECT at 25°C and 37°C (Fig. 2g) indicate that the system is a self-supporting gel at both temperatures, with the phase transition of pNIPAAm copolymer demonstrated visually by the shift to a white opaque appearance. Overall, our results indicate that the material's temperature responsiveness was successfully implemented without compromising 3D printing accuracy.
Small Angle X-ray Scattering (SAXS) analysis of EXPECT
We performed a structural analysis utilizing SAXS to explore the EXPECT system, aiming to understand and interpret its thermoresponsive characteristics. To achieve this, we collected scattering data at different concentrations and temperatures, encompassing both sub-LCST (25°C) and super-LCST (37°C) conditions, in the q range, which is sensitive to structures ranging from 27 to 140 nm in size (Fig. 3, Supplementary Fig. 2, and Supplementary Fig. 3). From our SAXS measurements, we extracted the Guinier slope (denoted as -n = dm, representing the fractal dimension, as shown in Fig. 3a). This parameter offers insights into the arrangement of polymer chains and surface characteristics within the EXPECT system or its individual components.
As a starting point, we ran the SAXS analysis on pNIPAAm homopolymer at a concentration of 1.5% (w/v). We observed a significant change in dm as a function of temperature, finding 1.5% pNIPAAm transitioning from 1.17 at 25.0°C to 4.09 at 37°C, as shown in Fig. 3b. This shift indicated a transformation from a swollen, disentangled polymeric network to collapsed globules with a nearly smooth surface (as illustrated in Fig. 3a i and ii) [55]. This transition pattern remained consistent for all three tested concentrations of pNIPAAm homopolymers (0.3%, 1.5%, and 3% (w/v), as shown in Supplementary Fig. 2a and 2b), with minor variations in the Guinier slope. Previous SAXS analyses of pNIPAAm microgels have reported similar -n slopes (or dm values) with values of approximately 1.5 below the LCST and around 4 above the LCST [56, 57]. This suggests that SAXS analysis is a suitable method for investigating temperature transitions in this system.
Further SAXS studies were performed on aqueous mixtures of pNIPAAm homopolymers with CP and gelatin. At 25°C, a dm value of 1.04 was extracted for this mixture (Supplementary Fig. 3a, blue), indicating rod-like polymeric chains with weak physical entanglement (Fig. 3a iii). As the temperature increased to 37°C, the dm of the mixture shifted to a value of 3.14 (Supplementary Fig. 3a, red), signifying a transition to a collapsed, entangled globular structure with a rough surface (Fig. 3a iv). In the absence of pNIPAAm, a mixture of CP and gelatin yielded dm values of 1.11 and 1.08 at 25°C and 37°C, respectively (Supplementary Fig. 3b), indicating an extended coil structure at both temperatures (Fig. 3a v and vi). This confirms that the observed temperature responsiveness is attributable to pNIPAAm.
We then proceeded to evaluate the SAXS measurements for pNIPAAm-CS to understand how the addition of chondroitin sulfate affects this transition. At 25°C, the dm values were 2.39 and 2.17 for concentrations of 1.5% and 3% (w/v), respectively (Fig. 3c and Supplementary Fig. 2c). These findings suggest increase in internal entanglement or chain branching (Fig. 3a vii [55]) due to the chondroitin sulfate grafting [54] and conform to the expected structure of the synthesized copolymer in the swollen state at 25°C [58]. Despite this structural change, a transition from dm = 2.39 to 2.66 was still observed for pNIPAAm-CS at 25°C and 37°C (Fig. 3c), indicating that the copolymer retained its temperature sensitivity even after the grafting with chondroitin sulfate (Fig. 3a viii). Furthermore, this transition was consistent for both concentrations (1.5% and 3% (w/v), as shown in Supplementary Fig. 2c and 2d), indicating that the transition in pNIPAAm-CS is not concentration-dependent in the concentration range investigated.
The full EXPECT composition (3% (w/v) pNIPAAm-CS + 0.8% (w/v) CP + 1% (w/v) gelatin) also showed a structural transition with temperature (dm from 2.16 → 2.62, Fig. 3d-e) (Fig. 3a ix and x). Without gelatin, 3% (w/v) pNIPAAm-CS mixed in aqueous solution with 0.8% (w/v) CP underwent transition (dm from 1.55 → 2.85, Fig. 3d-e). Interestingly, pNIPAAm-CS in aqueous solution with 1% (w/v) gelatin did not show any transition (dm from 2.51 → 2.46). However, this does not necessarily indicate a complete absence of transition but rather suggests that a less pronounced molecular re-arrangement occurs, considering that both components are thermoresponsive, adding complexity to the structural events since two different molecules are involved.
The transformation of pNIPAAm-based polymers, shifting from a collapsed globule to a swollen, disentangled polymeric network upon cooling through the LCST, has been extensively documented in connection with a transition from adhesive to non-adhesive behavior towards cells and proteins [30–33]. It is suggested that this process occurs in two stages: initially, cell detachment is triggered by the hydration of polymer chains, followed by a cell shape change that relies on F-actin dynamics, facilitating further detachment [59]. This detachment phenomenon has been observed on 2D surfaces for various cell types, including endothelial cells [60, 61], hepatocytes [60], keratinocytes [62], cardiomyocytes [63, 64], tenogenic and osteogenic progenitor cells [34] and has also been demonstrated for extracellular matrix (ECM) proteins [65, 66]. We anticipate that these thermosensitive properties, as validated through our combined rheological and SAXS data, play a crucial role in orchestrating the aligned cellular movements of cells embedded within EXPECT, described in the following sections.
Pilot cell studies with EXPECT
We conducted two pilot studies using an L929 murine fibroblast line. In the first study, we suspended 1 × 106 cells/mL in a sacrificial 6% (w/v) gelatin bioink and extruded them forming 10 mm ring-patterns embedded in the finalized composition of EXPECT (3% (w/v) pNIPAAm-CS + 0.8% (w/v) CP + 1% (w/v) gelatin). The cultures were maintained under standard conditions of 37°C, 5% CO2, and 90% relative humidity for 14 days to assess the biocompatibility of the embedding medium.
Using fluorescent confocal imaging, we observed TRITC-conjugated phalloidin stained cytoskeleton (red) and DAPI-stained nuclei (blue, Supplementary Fig. 4a). Between days 1 and 7, we noticed a tendency for cell spreading and alignment of the F-actin cytoskeleton. This coincided with an increase in metabolic activity (Supplementary Fig. 4b) and DAPI count (Supplementary Fig. 4c, right axis). However, between days 7 and 14, although the DAPI count did not increase, we observed a notable increase in the disorganization of the patterns, with increased spacing between the cells. To quantify the changes in the pattern width, we converted the DAPI/phalloidin images to grayscale, manually rotated the images to align the major axis of the printed structure to the y-direction and summed the pixel intensity along the column (i.e., the y-axis of the images). From the intensity data, we measured the width of the columnar pattern (i.e., width of the peak) at each time point. We found a statistically significant increase from 212 ± 78 µm at day 7 to 499 ± 70 µm at day 14 (Supplementary Fig. 4c, left axis, p < 0.0001). Also observed in previous reported studies, we postulate that this increase in pattern width over the culture period reflects stochastic movements of the cells in vitro and which gradually worsens the initially placed patterns [7, 14, 16].
To address the issue of pattern maintenance over the culture period, our subsequent study focused on investigating whether the temperature-sensitive behavior of EXPECT could offer a solution. In the next pilot study, we extruded 4 × 106 cell/mL fibroblasts in the same ring pattern. After fabrication, a temperature-actuated study group was designated to undergo a cooling cycle to 25°C for 15 minutes at days 1, 3, 5 and then once every 5 days thereafter, while otherwise maintaining normal culture conditions for a 21-day culture period. As a comparison, a static group was kept at normal culture conditions throughout the study.
Co-images of DAPI and phalloidin from day 21 (Fig. 4a) qualitatively showed that the cells in the temperature-actuated group tended to aggregate more. Quantitative analysis of the imaging data revealed that cell area was slightly lower under temperature-actuated compared to static conditions (Fig. 4b), possibly an indicator of increased cell spreading under static conditions. This behavior aligned with our hypothesis that EXPECT is cell-adhesive at 37°C. Additionally, we segmented the nuclei present in a single image through StarDist plugin for ImageJ [67], extracted the center of mass of each nucleus, and used these coordinates to calculate the nearest neighbor to each particle (Fig. 4c). The results showed that the distance between the nearest neighboring particles was smaller, and the distribution was tighter for the temperature-actuated group compared to the static group. Based on these data, we anticipated that temperature actuation could be useful in minimizing the accumulated migratory distances between cells. We decided to continue our investigations by applying this method with primary bone-marrow derived MSCs to explore if it reduced the migratory distance in this cell type.
Temperature-actuated aligned intercalation of MSC single cells
We next prepared a single cell suspension of 4 × 106 cell/mL MSCs in 6% (w/v) gelatin, 3D printed the same free-floating ring pattern into EXPECT, and maintained the constructs under static or temperature-actuated conditions for 36 days. Temperature actuation was introduced on days 1, 3, 5, and subsequently every 5 days. Co-images, combining DAPI and phalloidin, were captured on days 0, 7, and 36 post-printing. The experimental timeline is illustrated in Fig. 5a.
From fluorescence imaging on day 1, immediately after printing but prior to temperature actuation, the cells were distributed across a pattern width of 390 ± 114 µm. Under static conditions, by day 36, we observed some reductions in cellular organization along the channel. The pattern width significantly increased to 537 ± 76 µm (p = 0.0235) and appeared less regular (Fig. 5b). By adding temperature-actuation to EXPECT, we observed a striking difference in cell behavior. By day 7, the cells tended to self-organize into longitudinal patterns along the embedded channels, and this configuration was maintained up to day 36. To analyze this further, we conducted image analysis to measure the width of the columnar patterns at each time point. As can be observed from the day 7 data depicted in Fig. 5c, the narrower distribution of the curves for temperature-actuated conditions indicated a higher degree of anisotropy compared to static conditions. Additionally, we normalized the maximum width of the columnar pattern (Fig. 5d) at days 7 and 36 for each sample group to day 0, immediately after printing. Interestingly, under temperature-actuated conditions, pattern width was reduced by approximately 50% by day 7, and this reduction was sustained at day 36. In contrast, when examined under static conditions, we noted a 115 ± 20% increase in pattern width by day 7 and a further 137 ± 20% increase by day 36 compared to day 0 (p < 0.001). These findings were consistent in two additional independent experiments, one conducted with an extra donor (see Supplementary Fig. 5).
In summary, we successfully induced and sustained lateral intercalation of cells within the EXPECT system, as evidenced by the reduced width of their columnar patterns. Building upon previous findings in pNIPAAm-based systems [30–33], we posit that that the observed temperature-induced effects are associated with the potential periodic release of cells and their endogenously expressed proteins, including chemoattractants, from the surface of the pNIPAAm-based channel, combining to decrease lateral migration away from the channels. Changes in the viscoelastic properties or fluid dynamics resulting from temperature-induced alterations in the polymer network also may have played a role in guiding and sustaining cells in the longitudinal pattern.
In the realm of tissue engineering, the significance of cell proximity and 3D patterns for tissue functionality is well acknowledged. While various methods, such as non-adhesive micro-patterned surfaces [68, 69], acoustic waves [9, 10], or pressure [70, 71] have been developed to guide initial cell assembly, limited attention has been given to the maintenance of these spatial arrangements. Singular applications or static maintenance of these cues are insufficient to drive cell movements to condense along a predetermined axis of alignment over extended in vitro culture periods beyond 7 days [7, 14, 16, 18, 72]. In contrast, our approach introduces a contactless method, leveraging materials design and minor temperature changes, that can be applied periodically during culture to drive and maintain cellular assembly in desired patterns.
This innovation contributes to the field of tissue engineering by addressing the nuanced challenge of maintaining functional patterning and offering a novel perspective for gaining mechanistic insights into the impact of spatial factors on functional bioassembly across prolonged culture periods. Notably, our approach has also led to the formation of cell distributions with lateral widths smaller than those produced by the 3D printer nozzle. Consequently, we propose this method as an enhancement to improve the resolution of extrusion-based 3D printing, one of the most widely utilized and adaptable techniques in biomanufacturing.
Temperature-actuated convergent extension and bioassembly of MSC spherical aggregates
We further investigated the effect of EXPECT temperature actuation on MSCs assembled into spherical aggregates of approximately 200,000 cells using AggreWell™ plates. Using the same cell density (4 × 106 cell/mL) and ring size, we compared the changes in aggregate roundness, pattern width, and length under static versus temperature-actuated conditions. On the first day after bioprinting, the individual aggregates appeared slightly oblong, probably due to the mechanical pressure during the bioprinting. Over the 36-day culture period, we found no statistically significant changes in aggregate roundness or pattern width under static conditions (p > 0.48). However, we did observe some dissociation among the aggregates during the culture period, as they appeared looser over time (Fig. 6a). There was no measurable evidence of the aggregates migrating towards each other or fusing, as the average longitudinal length of any continuous pattern remained approximately 200 µm throughout the culture duration (p > 0.999). These outcomes are in line with our expectations and correspond with the current findings in controlled assembly [73–76], where aggregates are typically placed adjacent to each other to achieve fusion, whereas, in our study, we spaced them apart at the start of the experiment.
In contrast, under temperature-actuated conditions, aggregate roundness decreased from 0.77 ± 0.14 µm to 0.46 ± 0.10 µm between days 1 and 7 (p < 0.0001). Moreover, pattern width tended to decrease within this same time period from 181 ± 33 µm to 156 ± 46 µm (p = 0.64). These changes in shape descriptors denote the narrowing of the cellular patterns, consistent with what we observed by day 7 in the previous study with single cells. Additionally, the aggregates took on a bipedal morphology by day 7 (Fig. 6a, yellow arrows). Remarkably, with temperature actuation, the shape changes were measured again on day 36, with mean roundness further decreasing to 0.23 ± 0.14 (p < 0.0001) and pattern width tending to decrease to 119 ± 22 µm (p = 0.18). Additionally, we measured statistically significant increases in the longitudinal length of the continuous cellular pattern over the 36-day culture period from 206 ± 44 µm at day 0 to 480 ± 158 µm at day 36 (p < 0.0001). This change is attributable to the elongation of the individual aggregates, which enabled instances of fusion of the aggregates at day 36.
The temperature-actuation of EXPECT exerted a comparable influence on both cell aggregates and individual cells, facilitating organized self-assembly along a specific axis. Under temperature actuation, spheroidal groups of MSCs underwent reshaping, adopting a bipedal morphology akin to what has been documented on 2D gel surfaces—a characteristic associated with directionally migrating spheroids [18]. Importantly, the persistent cell motility observed in the referenced study was maintained for 1–6 hours. The presence of a bipedal morphology in cell aggregates with temperature actuation until day 36, in contrast to the isotropic disorganization observed under static conditions, suggests that temperature actuation sustained the directional migratory behavior of MSCs.
Within the domain of biofabrication, there has been a notable emphasis on employing bottom-up [77, 78] and middle-out [79] assembly strategies that utilize spheroidal aggregates as foundational units for constructing larger musculoskeletal or organ structures. In these approaches, spheroidal cell aggregates, often MSCs [9, 51, 80, 81] are strategically positioned within hydrogels, either through embedded 3D printing or acoustic vibration, with the goal of promoting fusion through spatial proximity. However, achieving precise geometric fusion poses inherent challenges for several reasons. Firstly, it requires a high level of positional accuracy for the aggregates, placing substantial technical demands on the biofabrication process. Additionally, spherical tissue building blocks inherently exhibit isotropic and amorphous qualities and are susceptible to random disassociation over time [16, 18, 72, 82]. As a result, the methodology developed in this work has the potential to enhance anisotropic directional microorganization in bottom-up and middle-out strategies intended to scale up to larger tissues.
The temperature-sensitive characteristics of pNIPAAm have been explored in conjunction with various cell types and diverse copolymer chemistries [30, 34, 35, 83]. As a result, the applications of this method can potentially extend beyond MSCs and fibroblasts, exerting a broader influence in multiple areas of tissue engineering. The ability to achieve directional cell orchestration and maintain spatial axial assembly holds significant promise for enhancing signal transmission efficiency and facilitating the complex maturation of tissues with axial organization, including those involved in tendon, bone, neural, and cardiovascular repair. Going beyond its applications in tissue engineering, our approach involves initiating and dynamically sustaining cellular movements reminiscent of convergent extension [6] and directed cell migration [84]. Consequently, the methodology in this study offers a valuable model framework for exploring developmental mechanisms, particularly those linked to axial elongation. Moreover, it holds promise for investigating diseases or repair mechanisms where the directed chemotaxis of cells is crucial, such as in the case of cancer metastasis.