1. Introduction
Cancer, as a cause of death, is only surpassed by cardiovascular diseases [
1]. According to the World Health Organization, in 2020, there were 19.3 million cases worldwide, and about 10 million people died of cancer. Between 2005 and 2015, the number of oncology patients increased by 33%, mainly due to population aging and population growth. Nevertheless, cancer mortality decreased in the same period in 140 of 195 countries or territories monitored in the Global Burden of Disease 2015 Study, demonstrating that progress is being made in cancer treatment and prevention [
1].
Given the complexity and heterogeneity of cancer, various therapeutic targets are being investigated, including components of the tumor milieu. The tumor microenvironment (TME) consists of several types of cells (tumor-associated fibroblasts, immune cells, mesenchymal stem cells, adipocytes, and vascular cells) embedded in extracellular matrix (ECM) soaked by interstitial fluid rich in soluble factors secreted by cells [
2].
A compelling body of evidence indicates that tumor progression depends on the interaction between the tumor and its microenvironment [
3]. It is known that the effectiveness of anti-cancer therapies is modulated by changes in the TME [
2,
4]. Therefore, extensive research efforts have been devoted to investigating the spatial organization of the native TME [
5] and to build in vitro models of the TME using three-dimensional (3D) bioprinting [
6] and lab-on-a-chip techniques [
7].
Conventional cancer models, such as 2D cultures of cancer cells and stromal cells, cannot mimic the TME, especially in what concerns oxygen and nutrient availability, thereby being limited as biomimetic models for anti-cancer drug testing or fundamental research in cancer biology [
8]. The lack of tumor-stroma interaction is also a weak point in patient-derived tumor organoid models otherwise highly valued for their ability to replicate intra- and intertumor heterogeneity in high-throughput cultures of 3D constructs [
9].
Bioprinting enables the precise positioning of cells and biomaterials in a computer-controlled, layer-by-layer procedure that preserves cell viability [
10]. Therefore, it is considered a promising avenue towards building tissue constructs that recapitulate the compositional and geometric complexity of the TME [
11,
12,
13,
14,
15,
16,
17,
18,
19,
20]. Although none of the model tissues constructed so far have been able to replicate all the features of a TME, most of them provided valuable insights into cancer biology. For example, in a rectangular grid of hydrogel threads uniformly loaded with HeLa cells from a cervical adenocarcinoma cell line, multicellular spheroids formed in less than fivedays of culture, which proved to be more resistant to paclitaxel (an antineoplastic chemotherapy drug) than HeLa cells in 2D culture [
20]. In addition, a bioprinted cancer-on-chip model recapitulated the biochemical and biophysical features of glioblastoma and reproduced the clinically observed, patient-specific variability of the response to concurrent administration of chemotherapy and radiotherapy [
18].
The present work aimed to develop 3D bioprinting methods for building in vitro models of the TME made of cells uniformly distributed in an artificial hydrogel matrix.Our working hypothesis was that the tumor couldbe represented by a cancer cell-laden hydrogel droplet, whereas its microenvironment can be modeled by rings of a hydrogel loaded with peritumoral cells. To our knowledge, this is the first report of a nanocellulose-alginate hydrogel being used for the bioprinting of breast cancer models. In this study, we included tumor-associated fibroblasts (TAFs) and, for the first time, peripheral blood mononuclear cells (PBMCs) in the TME model. The design of this study is represented schematically in the Supplementary Material (
Figure S1).
TAFs are the most abundant stromal cell type in various solid tumors (e.g., breast, pancreatic, and colorectal) [
21,
22]. They synthesize most of the tumor’s ECM and modulate tumor progression by secreting growth factors, cytokines, chemokines, and matrix-degrading enzymes [
23]. In mouse models of pancreatic and lung carcinoma, it has been established that TAFs block intratumoral infiltration of lymphocytes by highly expressing CXC-chemokine ligand 12 [
2]. Hence, TAFs also play an immunosuppressive role. There is solid evidence that TAFs originate from bone marrow-derived MSCs [
24] and present immunophenotypic similarities with them [
25]. Recent works, however, suggest that TAFs emerge from cancer stem cells, a subpopulation of tumor cells capable of self-renewal and differentiation [
26,
27].
PBMCs are a heterogeneous population of immune cells. The co-culture of PBMCs and epithelial tumor organoids allowed for the expansion of tumor-reactive T-cells [
28]. We chose to include PBMCs in the TME tissue constructs because several studies have demonstrated that the clinical progression of cancer depends on the level of immune cell infiltration of tumors. A strong lymphocytic infiltration was found to be associated with longer disease-free survival (after surgery) and/or overall survival in many types of cancer, including melanoma, head and neck, breast, prostate, colon, bladder, and lung cancer [
29]. In lung cancer patients, the gene expressionprofile of PBMCs measured before tumor resection was found to be an independent predictor of patient survival: higher expression of genes associated with protein synthesis was an indicator of better survival, whereas higher expression of genes associated with cell cycle was linked to worse survival [
30].
Besides the characterization of the bioprinted tissue constructs, this work also presents computer simulations of post-printing structure formation. Such simulations are necessary because cells take advantage of their motility to form clusters of optimal stability and function [
31]. Hence, 3D bioprinting merely provides a framework for the directed self-assembly of the cells delivered by the printer.
A notable limitation of the TME model proposed here is the lack of endothelial cells in the peritumoral medium. The motivation of this option is twofold: (i) it allows for the optimization of the co-culture conditions for cancer cells, TAFs and PBMCs, and (ii) it differentiates this study from previous investigations [
11], thereby providing complementary information. From the immunological point of view, the TME model avoids the possible interactions between tumor cells and immune cells (PBMCs) based on non-self recognition, leading to a cytotoxic and destructive mechanism. This was possible due to the selection of TME-comprising cells to be major histocompatibility complex (MHC)-compatible with the tumor cells. This makes the proposed model more reliable for replicating in vivo tumor development.
2. Materials and Methods
2.1. Preparation of Malignant Cells and Peritumoral Cells
2.1.1. Culturing SKBR3 Cells
Tumor cells from SK-BR-3 cell line (ATCC® HTB-30™, Lomianki, Poland) were cultured in McCoy’s 5a medium modified (Gibco BRL, Invitrogen, Carlsbad, CA, USA) supplemented with 10% fetal calf serum (FCS, PromoCell, Heidelberg, Germany) and 1% penicillin/streptomycin antibiotic solution (Pen/Strep, 10,000 IU/mL; PromoCell) using adherent culture flasks and incubated at 37 °C in 5% CO2 atmosphere. When reaching 90% confluence, the cells were detached from the plastic surface with 0.25% (w/v) trypsin EDTA solution (Sigma-Aldrich, St. Louis, MO, USA) and further expanded at a subcultivationratio of 1:2, as recommended by the provider’s protocol.
2.1.2. Harvesting and Culturing Tumor-Associated Fibroblasts (TAFs)
Human TAFs were isolated using the collagenase type IV-S prepared from
Clostridium histolyticum (Sigma-Aldrich), as previously described by Paunescu et al. [
25]. In brief, surgical pieces of approximately 5 cm
2 were obtained from breast cancer female patients diagnosed in different carcinoma stages. Tissue-isolated cells were washed several times with phosphate-buffered saline (PBS) solution, were successively passed through 70/40 µm strainer filters and were cultivated in adherent plastic culture plates. The culture medium contained α-minimum essential medium (α-MEM; Gibco), 10% fetal calf serum (FCS; PromoCell) and 1% antibiotics solution (penn/strep, 10,000 IU/mL;PromoCell), and the tumor-associated fibroblasts were placed in cell culture incubator (37 °C and 5% CO
2).
2.1.3. Harvesting and Culturing PBMCs
PBMCs were obtained from 10 mL venous peripheral blood by centrifugation in a density gradient. The blood sample was harvested from the breast cancer patients before the surgical intervention in EDTA collection tubes. For PBMCs isolation, we used 10 mL Ficoll-Paque PLUS (Sigma-Aldrich), which was placed on the bottom of a 50 mL Falcon tube, while on top, we slowly pipetted the peripheral blood diluted with PBS (Gibco) at a ratio of 1:1. We collected the mononuclear cells ring after centrifugation at 500× g and deceleration 0 for 25 min. After washing the PBMCs twice with PBS (Gibco), the cells were cryopreserved in liquid nitrogen (−196 °C) in a medium containing FCS (PromoCell) and 10% dimethyl sulfoxide (DMSO; Sigma-Aldrich), at a concentration of 106 cells/mL, for further use.
The blood samples were submitted for HLA-typing in an independent study, and we selected only the PBMCs and TAF from patients with HLA type A*11 to be compatible with the SK-BR-3 tumor cell line.
All tissue and biological samples were obtained after signing the informed consent elaborated under a protocol approved by the Ethical Commission of the County Emergency Hospital “Pius Brinzeu” Timisoara, according to the World Medical Association Declaration of Helsinki.
2.2. Three-Dimensional Bioprinting of Model Tissues
In this study, we employed an INKREDIBLE bioprinter (CELLINK, Gothenburg, Sweden) to deliver cell-laden hydrogel strands in a computer-controlled arrangement at 10 µm in-plane resolution, and 100 µm layer resolution. This bioprinter is equipped with two pneumatic print heads, whose extrusion rate depends on the pressure applied to the piston, the geometry of the extrusion nozzle and the rheological properties of the bioink. An external compressor, connected to a different outlet than the bioprinter, provides pressurized air at a pressure of up to 250 kPa. The pressure applied to each print head can be set at a fraction of this pressure via two dials. Using blunt needles of 0.6 mm internal diameter as print nozzles, typical print-head pressures employed in our experiments ranged between 20 kPa and 100 kPa.
2.2.1. Preparation of Cell-Laden Hydrogels
For the preparation of cell-laden hydrogels, the three cellular types were processed individually according to their specific requirements. The adherent cells (TAFs and SK-BR-3 cells) were processed using the trypsinization method. Shortly, when reaching 90% confluence, the culture flasks were washed with warm PBS (Gibco), and 0.25% (w/v) trypsin EDTA solution (Sigma-Aldrich, St. Louis, Mo, USA) was added, 1 mL/1 cm2 of culture surface. After centrifugation (10 min, 500× g), the cells were resuspended in a culture medium and counted. The cryopreserved PBMCs were thawed, centrifuged (10 min, 500× g) and counted. Then, TAF and PBMCs were mixed in the same cellular suspension, as they were printed with the same printing head. The cellular mixture between TAF and PBMCs contained 1.5 × 106 cells of each type, suspended in 300 µL of RPMI (Gibco), 10% FCS (PromoCell) and 1% antibiotics solution (penn/strep, 10,000 IU/mL;PromoCell). The SK-BR-3 cells were also resuspended in 300 µL of the same culture medium. Using the cell mixer kit (CELLINK), 3 mL of CELLINK Bioink (universal bioink, CELLINK) was used for uniformly mixing 300 µL of TAFs and PBMCs (3 × 106 cells), following the provider’s protocol; other 3 mL of CELLINK were mixed with 300 µL of SK-BR-3 cells (3 × 106 cells), following the same protocol. As a result of the mixing procedure, we obtained two printing cartridges of 3 mL containing cell-laden hydrogel: one with TAFs and PBMCs mixed in equal proportions and the other with SK-BR-3 cells. In both cartridges, the cellular density was 106 cells/mL.
2.2.2. Digital Modeling and Bioprinting
The Rhinoceros
®software(Robert McNeeland Associates, Seattle, WA, USA) was used to build digital models of two tissue constructs. The first one, called hereafter the toroidal structure, consists of a spheroid tightly embraced by a donut-shaped hydrogel ring; the second, the so-called triple-layered construct, comprises the toroidal structure sandwiched between two hydrogel layers—a pair of concentric rings at the bottom and one ring at the top (see
Supplementary Materials, Figure S2). The pyramid shape of the triple-layered construct is advantageous because it ensures the stability of the hydrogel rings stacked on top of each other; therefore, this geometry is more reproducible than a cylindrical or spherical one.
We did not rely on a slicer software to generate the G-code needed to control the bioprinter. Instead, we employed the digital model as a template for writing the G-code, thereby retaining full control of the path and movement speed of the print head.
During printer calibration, we ascertained that the first print head nozzle was 0.1 mm above the target surface. By printing straight lines at the feed rate used in the actual bioprinting process, we identified the air pressure applied to the piston of the print head for each type of hydrogel employed in our experiments. The pressure was adjusted until the printer delivered a uniform gel strand about 20% higher in diameter than the inner diameter of the extrusion nozzle.
To deliver a hydrogel ring, we employed the G1 command, substituting circular arcs with the corresponding chords (i.e., a circle was written as a succession of 12 straight lines), but the result was a smooth ring due to bioink viscosity (
Figure 1a). The feed rate was set to 150 mm/minute during extrusion and 900 mm/minute during the repositioning of the print head from one construct to another. Before drawing a circle, a dwell time of 0.3 s was allocated, via the G4 command, to initiate hydrogel extrusion. We used Konix ultrasound transmission gel (Turkuaz, Istanbul, Turkey) in testprints meant to establish the settings of the bioprinter. The tumor was modeled by Konix gel stained with Neutral Red dye (Sigma-Aldrich, Hamburg, Germany). Supplementary Video S1 illustrates the entire printing process.
Figure 1 illustrates the reproducibility of the 3D printing process. As shown in
Figure 1c, the tissue construct is about 3 mm in diameter, and the central spheroid (red) is about 0.8 mm in diameter. The actual bioprinting took place in the 24-well cell culture plate, and the tall walls of the wells hampered visualization (
Figure S3).
2.3. In Vitro Culture of Model Tissues
The 3D bioprinted tissue constructs were placed individually in wells of a 24-well culture plate and grown for 14 days in a medium containing RPMI (Gibco), 10% FCS (PromoCell) and 2% antibiotics solution (penn/strep, 10,000 IU/mL;PromoCell). The medium selection reflected the constructs composition of three cell types that can be cultivated in RPMI, but we increased the amount of the antibiotic as the constructs suffered multiple processing steps.
2.4. Histology
TAF immunocytochemistry staining was performed as previously described [
25] on adherent cells for confirming their phenotype, using the following monoclonal mouse anti-human primary antibodies from Dako company (Carpinteria, CA, USA): alpha-smooth muscle actin (α-SMA), CD105 (endoglin; clone SN6h), vimentin (clone VIM 3B4), and pan-cytokeratin (Pan-CK; AE1/AE3).
The bioprinted tissue constructs were sectioned and evaluated by two histological methods: Hoechst and hematoxylin–eosin (HE) staining of cryosections [
32].
Cryosections were prepared from tissue constructs cultured for 14 days. The samples were covered in tissue freezing medium (Leica Biosystems, Wetzlar, Germany), snap-frozen in liquid nitrogen and sectioned at 5 µm thickness using a Leica CM 3050 S Cryostat (Leica) cryotomy instrument. The cryosections were treated with Hoechst 33342 (Invitrogen) nucleic acid stain, and a cover slide was added using Fluoromountaqueous mounting medium (Sigma-Aldrich). When the HE stain was used, the hematoxylin solution (Sigma-Aldrich) colored the nuclei of the cells for 2 min, and the eosin (Sigma-Aldrich) was placed to stain the cellular cytoplasm for 30 s, followed by distilled water washing procedures. The HE stained sections were mounted with gel mount aqueous mounting medium (Sigma-Aldrich). The sections stained with either Hoechst or HE were visualized using ZEISS Axio Observer microscope (Zeiss, Munich, Germany).
2.5. Computational Modeling of the Evolution of Model Tissues
To simulate the evolution of bioprinted model tissues, we incorporated a new algorithm in the
SIMMMC application [
33]. This algorithm is an extension of the “SIMMMC for bioprinting” module [
34], which describes shape changes of live systems made of two types of cells in an environment that contains a cell culture medium, a hydrogel that can be remodeled by cells, as well as a biomaterial that cannot be penetrated or altered by cells. Besides the Metropolis Monte Carlo (MMC) algorithm [
35,
36], the extended module also includes events that account for cell proliferation.
Our computational model is built on a 3D cubic lattice. Each lattice site is occupied by a point particle that stands for a cell or a similar-sized volume element of one of the three possible embedding materials included in the model. Site occupancy is specified by a particle-type index, , as follows: for medium, for hydrogel, for peritumoral cells, for cancer cells, and for the biomaterial. The latter, however, was not present in the experiments conducted in this study; therefore, the simulations only involved particles of type 0, 1, 2, and 3.
In the model system, interactions are described in terms of works of cohesion,
, and works of adhesion,
(a term used when
); these are defined as the mechanical work needed to break the bond between two adjacent particles whose types are specified by the pair of indices [
37]. The total energy of adhesion can be written as [
38]:
where
= 26 is the number of up to second-nearest neighbors of a given lattice site—called hereafter the neighbors of the site;
is the number of heterotypic bonds between particles of type
and
;
is the number of particles of type
, and
is the interfacial tension of the surface of separation between the two media, of type
and
. In Equation (1), the first term is the total interfacial energy of adhesion, whereas the second term is the total energy of cohesion.
We simulated the evolution of the bioprinted TME models using the following computational algorithm:
Compute the total energy of adhesion , given by Equation (1);
Commence the first Monte Carlo step (MCS) by identifying interfacial cells (i.e., particles of types 2 and 3 with neighbors of a different kind);
Take type 3 (cancer) cells, in a random order, and give them the chance to proliferate with a probability or to move (i.e., swap with a neighbor of a different type);
Pick type 2 (peritumoral) cells in a random sequence, and give them the chance to proliferate with a probability or to move (i.e., swap with a type 1 or a type 0 particle—note that allowing a type 2 particle to also swap with a type 3 particle would lead to a repeated rearrangement of the corresponding interface);
If the chosen event is proliferation, create a new cell of the same type by replacing a randomly picked neighbor of type 1 or 0;
If the event of choice is movement, compute the corresponding change in energy, , and accept the move if ; otherwise, accept it with a probability ;
Increment the number of MCS by one and start a new MCS.
In the expression of the trial move acceptance probability,
, (step 6 of the algorithm),
is the biological analog of the energy of thermal fluctuations from statistical physics, equal to Boltzmann’s constant multiplied by the absolute temperature; this quantity is a surrogate measure of cell motility [
39]. It is convenient to express works of adhesion and interfacial tensions in units of
.
4. Discussion
This work presents 3D printed models of the TME consisting of cancer cells incorporated in a central spheroid wrapped in hydrogel rings loaded with TAFs and PBMCs. Unlike in the toroidal tissue construct, in the triple-layered structure, the cancer cells were surrounded on all sides by peritumoral cells uniformly distributed in the CELLINK hydrogel, giving rise to a more consistent microarchitecture of cell clusters separated by cell-free domains.
The hydrogel used in this work was not prepared in-house. Instead, we used a commercial one—the CELLINK universal bioink (CELLINK, Sweden). According to the manufacturer, this hydrogel contains a mixture of nanofibrillar cellulose, sodium alginate, D-mannitol and HEPES buffer. We considered this bioink for our 3D cell culture because several reports demonstrated that nanocellulose-alginate hydrogels mimic well the extracellular matrix properties while providing low cytotoxicity [
41].
Histological data suggest that the cells present in the TME tissue constructs can remodel the hydrogel’s filament network via traction forces and the action of proteolytic enzymes—matrix-metalloproteinases (MMPs). This observation is consistent with findings related to cervical tumor models built by 3D bioprinting [
20]. In their work, Zhao et al. used a hydrogel made of 10% gelatin, 1% sodium alginate and 2% fibrinogen to incorporate HeLa cells (10
6 cells/mL). Cell spheroid formation was observed, 65% of them being larger than 50 µm in diameter after 5 days of culture. Moreover, cells cultured in 3D hydrogels secreted more than twice as much MMP-2 and MMP-9 asthose cultured in Petri dishes. Our work is similar to ref. [
20] in what concerns the initial cell density of the bioprinted constructs and the presence of sodium alginate in the hydrogel. Despite differences in construct geometry, biomaterials, and cell types, the TME model presented in this work recapitulates several aspects of the structure formation observed by Zhao et al., suggesting that the 3D environment plays a central role in the self-organization of cancer cells. Nevertheless, the differences observed between the toroidal and triple-layered tissue constructs in the present work suggest that the topology of the TME is also important. The proximity of TAFs and PBMCs has facilitated cell proliferation in the triple-layered construct, resulting in higher cell density and more consistent structure formation.
The geometry of our toroidal tissue construct is similar to the glioblastoma-on-a-chip model developed by Yi et al. [
18]. They employed cancer cells and endothelial cells suspended in hydrogel (5 × 10
6 cells/mL) to print a flat, compartmentalized tumor-stroma model: a cancer cell-laden hydrogel disk surrounded by a concentric hydrogel ring comprising human umbilical vein endothelial cells (HUVECs). Importantly, the model tissue was enclosed in a microfluidic device that prevented gas and nutrient exchange from the top and bottom, thereby creating radial concentration gradients and central hypoxia. In our study, bioprinted constructs were bathed by cell culture medium in 24 well cell culture plates, so hypoxia could develop only in the triple-layered construct, in which tumor cells were buried in the bulk of the structure, at depths of the order of 0.8–1.5 mm. In their work [
18], Yi et al. employed two types of hydrogel: collagen and another one prepared from decellularized pig brain extracellular matrix (BdECM). In BdECM, cancer cells had a higher proliferation rate and higher expression of pro-angiogenic factors and ECM-remodeling proteins. In light of these findings, we plan to investigate TME models made of cells suspended in different hydrogels.
To avoid the impact of the embedding hydrogel, Langer et al. investigated tumor heterogeneity using TME models built by scaffold-free bioprinting [
15]. They printed a core of malignant cells wrapped in a shell of stromal cell mix (primary human mammary fibroblasts and HUVECs). A variety of TME models were obtained by using different breast cancer cell lines (MCF-7, SKBR3, HCC1143 and MDA-MB-231). Unlike in the present work, the bioink used in. [
15] included tunable hydrogels to assure mechanical support during tissue construct fabrication, but these were removed during subsequent culture, leaving behind a pure cellular structure. The bioprinted constructs of about 2mm × 2mm × 1mm in size evolved spontaneously into a cancer cell spheroid covered by a shell of stromal cells. MMC simulations of multicellular spheroids composed of two cell populations of different cohesivities show, in accord with experiments [
31], that the system spontaneously forms core–shell structure, with the more cohesive population being wrapped by, the less cohesive one [
38]. A core–shell arrangement was also observed in 2D simulations of a high cell–density bioprinted model of the tumor microenvironment [
42]. In the present study, post-printing self-assembly did not result in a concentric arrangement of the constituent cells, neither in experiments nor in simulations, presumably because of the much smaller cell density of the bioprinted TME constructs.
It is unclear why the cells were less organized in toroidal tissue constructs than in triple-layered ones. The experiments indicated a discrepancy between the post-printing evolution of the two constructs (
Figure 3c,f), but the computational modeling did not reproduce it (
Figure 8b,e). A distinctive feature of the triple-layered structure is that it gives rise to gas and nutrient gradients, akin to the in vivo TME, but these are not considered in the computational model. Further research will be needed to test the impact of oxygen and glucose gradients on the evolution of malignant and peritumoral cell populations. One interesting option in this respect is to include the toroidal structure into a tissue-on-chip device [
18].
The TME plays a critical role in metastasis, too. In a recent study [
43], Kok et al. demonstrated that fibrotic niche generation by metastatic cells is decisive for polyclonal metastasis. When nonmetastatic and metastatic cells were co-transplanted to the spleen, both cell types migrated into the liver within three days and formed colonies there. The metastatic cells activated hepatic stellate cells, inducing fibrotic niche generation. In this environment, nonmetastatic cells continued to proliferate even after the selective depletion of the metastatic cells. Kok et al. suggest that targeting the niche may be an effective strategy for preventing the spread of solid tumors [
43]. Such hypotheses may be tested using model tissues that mimic the TME.
In the computational model employed in this work, the algorithm only enables superficial cells to proliferate (i.e., those cells located in the vicinity of hydrogel or cell culture medium). This requirement is a rough approximation of the experimental fact that cell proliferation mainly occurs close to the tumor periphery due to limited oxygen and nutrient availability within the bulk of the tumor [
18,
42].
The limitations of this study include (i) the relatively low cell density of the bioprinted construct, (ii) the lack of vascular cells in the TME, and (iii) the exclusive focus on cell adhesion in the computational analysis of cellular rearrangements. Extrusion bioprinting of bioinks with high cell density is feasible [
15], although it would require far more cells. We plan to increase cell density by an order of magnitude in future studies. The incorporation of vascular cells is an appealing objective, too. Recent progress in 3D bioprinting [
44,
45,
46] makes it possible to fabricate bulky tissue constructs with branched tubular structures within them; moreover, endothelial cells incorporated in the peritumoral environment self-organize into a microvasculature [
18]. Hence, the TME models could be endowed with a complex vasculature. Cell source options for this task include HUVECs [
11] and human-induced pluripotent stem cell (iPSC)-derived endothelial progenitors [
47]. In this respect, immunological aspects also need to be considered. The computational analysis should be extended to describe more types of peritumoral cells, as well as oxygen and nutrient gradients. A 3D version of the approach of Bustamante et al. [
42] seems to be a good starting point, although the differences between the algorithms employed in the cell-switching MMC model and the cellular Potts model make the implementation far from straightforward. If technical difficulties can be surmounted, the extended model might be able to predict the formation of a necrotic core in the center of large cell clusters, and cell proliferation might be simulated as a function of oxygen and glucose concentration.
5. Conclusions
This work presents tissue constructs fabricated by extrusion bioprinting as models of the tumor microenvironment (TME). The models proposed here differ from previous ones in geometry, matrix composition, and cellular composition. Our results confirm the hypothesis that scaffold-based models of the TME can be fabricated by extrusion bioprinting. Despite limited printer resolution, biomimetic structures emerge as a result of multicellular self-assembly.
The artificial hydrogel used in this study was appropriate for sustaining cell growth and proliferation. Cells migrated within the hydrogel and self-organized into heterotypic aggregates, mostly located on the periphery of the construct. Computer simulations suggest that the formation of aggregates comprising malignant cells and peritumoral cells can be explained by the differential adhesion hypothesis (i.e., by the tendency of the cells to form the largest possible number of strong bonds with one another as well as with the surrounding medium [
31]). The superficial localization of the aggregates, however, was not captured by the simulations. Further investigations will be needed to establish whether it can be explained by taking into account oxygen and nutrient gradients within the bioprinted construct.
In vitro engineered TME models of diverse compositions are promising tools for fundamental research that aims to decompose the key players and mechanisms involved in tumor growth and migration. In addition, TME models built from patient-derived cells may be used to test the effectiveness of anti-cancer therapies, thereby devising personalized treatment plans.