Skip to main content

ORIGINAL RESEARCH article

Front. Bioeng. Biotechnol., 21 February 2020
Sec. Tissue Engineering and Regenerative Medicine

Thermal Bioprinting Causes Ample Alterations of Expression of LUCAT1, IL6, CCL26, and NRN1L Genes and Massive Phosphorylation of Critical Oncogenic Drug Resistance Pathways in Breast Cancer Cells

\nAleli CampbellAleli Campbell1Jonathon E. MohlJonathon E. Mohl2Denisse A. GutierrezDenisse A. Gutierrez3Armando Varela-RamirezArmando Varela-Ramirez3Thomas Boland
Thomas Boland1*
  • 1Metallurgical, Materials and Biomedical Engineering, University of Texas at El Paso, El Paso, TX, United States
  • 2Department of Mathematical Sciences and Border Biomedical Research Center, University of Texas at El Paso, El Paso, TX, United States
  • 3Department of Biological Sciences, Border Biomedical Research Center, University of Texas at El Paso, El Paso, TX, United States

Bioprinting technology merges engineering and biological fields and together, they possess a great translational potential, which can tremendously impact the future of regenerative medicine and drug discovery. However, the molecular effects elicited by thermal inkjet bioprinting in breast cancer cells remains elusive. Previous studies have suggested that bioprinting can be used to model tissues for drug discovery and pharmacology. We report viability, apoptosis, phosphorylation, and RNA sequence analysis of bioprinted MCF7 breast cancer cells at separate timepoints post-bioprinting. An Annexin A5-FITC apoptosis stain was used in combination with flow cytometry at 2 and 24 h post-bioprinting. Antibody arrays using a Human phospho-MAPK array kit was performed 24 h post-bioprinting. RNA sequence analysis was conducted in samples collected at 2, 7, and 24 h post-bioprinting. The post-bioprinting cell viability averages were 77 and 76% at 24 h and 48 h, with 31 and 64% apoptotic cells at 2 and 24 h after bioprinting. A total of 21 kinases were phosphorylated in the bioprinted cells and 9 were phosphorylated in the manually seeded controls. The RNA seq analysis in the bioprinted cells identified a total of 12,235 genes, of which 9.7% were significantly differentially expressed. Using a ±2-fold change as the cutoff, 266 upregulated and 206 downregulated genes were observed in the bioprinted cells, with the following 5 genes uniquely expressed NRN1L, LUCAT1, IL6, CCL26, and LOC401585. This suggests that thermal inkjet bioprinting is stimulating large scale gene alterations that could potentially be utilized for drug discovery. Moreover, bioprinting activates key pathways implicated in drug resistance, cell motility, proliferation, survival, and differentiation.

Introduction

In vitro testing for drug discovery keeps making strides, especially with the advancement of genomics, proteomics, pharmacodynamics, bioinformatics, and automated High Throughput Screening (Andrade et al., 2016; Peng et al., 2017). Target-based drug design using appropriate cell assays, has not only transformed the identification of new targets, but it has also been supplemented with virtual testing aka “in silico” testing. With this approach, computer-based methods for drug simulations are utilized (Andrade et al., 2016; Hevener, 2018). In silico methods provide rapid and inexpensive techniques for quick lead test verification which proceed with in vitro cell testing. This method is a critical step in preclinical studies (Swinney and Anthony, 2011; Begley and Ellis, 2012; Peng et al., 2016). Previous studies have suggested that bioprinting can be used to model tissues for drug discovery and pharmacology (Peng et al., 2016, 2017). Peng et al., suggested that 3D bioprinting can help reduce the attrition rate in drug discovery by creating more realistic models. Through manipulation of pattern or anatomical models, it is possible to create permeable structures that ensure adequate delivery of nutrients and vascularization, which is primordial of in vivo environments. By bioprinting realistic models, we mean to generate tissue based on specific targeted characteristics such as lung, bone, cardiac, and even tumors.

While it is important to have a clear insight regarding cell viability and physiological changes of bioprinted (BP) cells, it is critical to understand the molecular changes within these cells in order to identify triggering mechanisms associated with cellular functions and behaviors. To our knowledge, this type of analysis has not been published before. Zhao et al., tested a 3D extrusion based bioprinted model of HeLa cells and found morphological differences, increased matrix metalloproteinase protein expression and higher cell proliferation when compared to the 2D standard cell culture. It is important to comprehend the gross anatomical structure as well as intra-cellular alterations to be able to model external stimuli, either of biological or synthetic nature. However, the comprehensive cellular response of bioprinted MCF7 breast cancer cells (BCC) or any other cells at the molecular level has not been published, yet it is crucial to determine whether bioprinted cancer models can potentially be used to predict drug efficacy, toxicity, and safety.

It has been widely suggested in the literature that bioprinting technology could lead to the pivotal discoveries of tissue engineered products which can be used for a range of clinical applications, e.g., skin grafting, tissue regeneration, cartilage repair, and others (Cui et al., 2012a, Yanez et al., 2014; Gudapati et al., 2016; Miri et al., 2019; Yerneni et al., 2019). However, this approach has not been used to develop tumor models in vitro for drug discovery. Recently Chen et al. and Phamduy et al. developed a bioprinting system where mass spectrometry was used in single printed cells. The authors (Phamduy et al., 2015) used laser direct-write cell bioprinting to bioprint MDA-MB-231 and MCF7s directly onto ex vivo rat mesentery tissue. They were able to monitor cell viability, proliferative and migratory properties and observed cell attachment and cell invasion within 2–5 days. Yet analyzing molecular and physiological changes in BP MCF7 BCCs has not been done and is long overdue. Here we report viability, apoptosis (programmed cell death), kinase phosphorylation, and RNA sequencing (RNA seq) analysis of BP MCF7 BCCs.

Experimental Procedures

Cell Culture

In this study, MCF7 (ATCC® HTB-22™) breast cancer cells were used for the in vitro experiments. Eagle's minimum essential medium (EMEM), supplemented with 0.01 mg/L Human recombinant insulin and 10% fetal bovine serum (referred to as media), trypsin 0.25% EDTA, and sterile phosphate-buffered saline (PBS) solution were used. Briefly, MCF7 cells were cultured per ATCC's cell protocol, a 75 cm flask with 8–10 ml of media were incubated in a humidified incubator maintained at 37°C with 5% CO2 until 80–90% confluency was reached cells were split and passaged to ensure cell stability.

Bioprinting Process

In preparation for bioprinting, MCF7 cells growing around 80–90% of confluency were gently rinsed with PBS (to remove dead cells and debris) and detached with trypsin harvested and centrifuged at 800 rpm for 5 min. The supernatant was removed, and the cell pellet was next rinsed in PBS. Then the cells were counted with trypan blue and a hemocytometer. Modified inkjet cartridges adapted for a modified HP thermal inkjet printer were used to for bioprinting purposes (Wilson and Boland, 2003; De Maria et al., 2013). Next, 100 μL of PBS-cell solution (~1.6 × 106 cells/mL—bioink) was added and printed into a tissue culture treated petri dishes (100 × 15 mm) or a Falcon, 96-well black/clear, tissue culture treated plate, flat bottom with lid. Following the bioprinting process, the cell viability analysis was conducted at 24 and 48 h utilizing two approaches: by using the Invitrogen™ Countess™ automated cell counter and through manual count by two different lab members using trypan blue and a hemocytometer. In these series of experiments, three independent measurements were accomplished, each performed in triplicate. Collected data are depicted as an average and standard deviation.

Examination of the apoptosis/necrosis pathway via flow cytometer: MCF7 Cells were bioprinted at a density of ~600,000 cells/dish in 6 ml of media. After 24 and 48 h post-bioprinting, cells were collected as above and double-stained by using the Annexin A5- FITC/propidium iodide (PI) staining kit, which is typically used to discern whether the cells are dying via apoptosis or necrosis pathway in the Beckman Coulter flow cytometer. The total percentage of apoptotic cells is depicted as the sum of both early and late stages of apoptosis (Annexin A5-FITC positives), whereas the cells stained only with PI, were considered as the necrotic cell population (Robles-Escajeda et al., 2016). Data acquisition and analysis were performed by using the Gallios flow cytometer (Gallios Beckman Coulter: Miami, FL) and the Kaluza software version 3.1. (Beckman Coulter) as previously detailed (Santiago-Vázquez et al., 2016).

Stain Process

Morphological characteristics of the BP and MS MCF7 cells were evaluated with laser-confocal microscopy. Twenty-four hours post-bioprinting, cells were stained as explained elsewhere; (Iglesias-Figueroa et al., 2019) briefly, cells were fixed in 4% formalin for 20 min, washed and incubated in 0.1% v/v Tween 20 in PBS for 10 min at room temperature (RT), washed twice more with permeabilizing solution, incubated again in 200 μl of 5% w/v in bovine serum albumin (BSA; sigma) dissolved in tris-buffered saline solution containing 0.5% v/v Tween 20 for 1 h at RT on a rocker platform. The cells were next stained with the primary antibody Neu (sc-33684, dilution: 1:50) overnight in a rocker platform at 5°C. The cells were next washed three times with permeabilizing solution and a 1:50 v/v dilution of secondary antibody goat anti-mouse IgG conjugated with Alexa Fluor™ 568 (Invitrogen) was added and then incubated for 1 h on a rocker platform at RT. They were rinsed three more times with permeabilizing buffer and posteriorly they were co-stained with 0.165 μM of Phalloidin Alexa Fluor™ 488 (Invitrogen), and 5 μg/ml of 4′,6-diamidine-2′-phenylindole, dihydrochloride (DAPI, Invitrogen) for 1 h on a rocker platform; rinsed three more times with 0.1% v/v Tween 20 in PBS and leaving 200 μl in the wells at the end. Finally, using an inverted confocal-laser-scanning microscope (model LSM 700: Zeiss; New York, NY) assisted by the Zen 2009 software (Zeiss), to acquire high-quality digital images in three fluorescence channels (Alexa 568, Alexa 488, and DAPI). The enhanced contrast Plan-Neofluar 40x/1.3 oil immersion, differential interference contrast objective was used. The parameters to capture the high resolution images were consistently maintained; single-plane images were consecutively scanned for each fluorescence channel setting the pinhole at 1 Airy Unit.

Phospho-MAPK Antibody Array (BP Cells) Process

The antibody array phospho-mitogen-activated protein kinase (phospho-MAPK) (R&D Systems #ARY002B) was used to analyze the phosphorylation levels of 26 kinases, which includes 9 MAPKs, ERK1/2, JNK1-3, and p28 isoforms in order to understand how the relative phosphorylation levels are affected by the bioprinting process. A collective sample of BP MCF7 BCC was harvested at 12 h, and 24 h post-bioprinting. Samples were processed as per the kit's protocol. Briefly, cells were solubilized in lysis buffer, diluted, mixed with detection antibodies, and incubated overnight in a rocking platform at 5°C with the phospho-MAPK array membranes. The membranes were washed the next day and detection reagents were applied. Upon completion of the reagent detection step, the iBright™ FL1000 Imaging System (Invitrogen Thermo Fisher Scientific) was used to develop the membranes. Microarray data was obtained by pixel density values, which were normalized, analyzed, and quality controlled by using Invitrogen™ iBright™ imaging software v.3.0 (Thermo Fisher).

RNA Extraction and Sequencing

RNA sequencing was conducted in BP MCF7 BCCs at 2, 12, and 24 h post-bioprinting with the intention of identifying genes that were upregulated or downregulated by this process. Samples for RNA sequencing were prepared following the same bioprinting protocol as mentioned above. BP MCF7 BCCs were next gently detached with a cell sweeper to avoid exposure to influencing solutions, centrifuged at 800 rpm for 5 min and collected at 2, 12, and 24 h. RNA extraction was conducted for each sample with the PureLink RNA mini kit from Thermo Fisher and used per manufacturer's protocol. RNA concentration was evaluated with a Nano-drop 2000 spectrometer (Thermo Fisher Scientific). RNA Integrity Number equivalent (RINe) was assessed for each sample and ranged from eight to ten. RNA seq data was analyzed for data summarization, normalization and quality control using Trimmomatic (version 0.36), Bowtie2 (version 2.2.5), and Cufflinks (version 2.2.1) (Langmead and Salzberg, 2012; Trapnell et al., 2012; Bolger et al., 2014). Differentially expressed genes were selected by using threshold values of >2 and <-2-fold change and a q ≤ 0.05. The q-value is used in genome-wide expression data, this statistical method is used to filter the proportion of false positives from a collection of p < 0.05.

STRING was used to map pathways of upregulated and downregulated differentially expressed genes. Cytoscape v.3.7.1., was next used to compare both network samples. Finally, network analysis to identify protein-protein interactions (edges or protein-protein connections) was extracted from Cytoscape. The number of connections or interactions are important parameters when targeting specific proteins associated with a disease. Genes with a >2-fold and <-2-fold (upregulated and downregulated genes) and a q < 0.05 were further classified by their gene ontology (GO) terminology. Upregulated and downregulated genes in the BP cells were divided into 3 groups, by their associated biological processes (bp), molecular functions (mf), and cellular components (cc). We also used STRING to create a protein network and compared the results of the MS to the BP MCF7 cells. We further mapped the protein pathways associated with genes expressed only in the BP cells. RNA and DNA information/analysis were compiled from: STRING, Panther, and Cytoscape 3.7.1 the following websites were also used to extract protein/gene connections: SMART, gene ontology, UniProt, and the Online Mendelian Inheritance in Man®.

Statistical Analysis

Continuous variables for cell response to bioprinting were summarized using means and standard deviation (SD) of triplicate samples. In this experiment, 3 independent samples at 2 different time frames (24 and 48 h) were used to calculate viability and a one-way ANOVA analysis was conducted to compare the means for each time frame of the BP and manually seeded MCF7 BCCs, negative controls were also conducted. For the antibody array assay, data was centered, normalized, and clustered utilizing the iBright™ Analysis software version 3.0. All statistical analysis were completed using Minitab 18, IBM SPSS Statistics 25 and excel 2013.

Results

Post-bioprinting Cell Viability and Apoptosis

The mean cell viability from the manual cell count process of the BP cells was 73.5% tested at 24 and 48 h and the average percentage of cell viability for MS (non-printed) cells was 98.6%. The average percentage cell viability from the automated cell count of the BP cells at 24 and 48 h was 76.7 and 72.8%, respectively. In the flow cytometer, samples collected at 2, 24 h and 7 days post-bioprinting, showed a percentage rate of live cells of 70, 30, and 47%, respectively. Apoptosis rates at these time frames were 30, 69, and 29%, respectively, and are shown in Figure 1A.

FIGURE 1
www.frontiersin.org

Figure 1. Percentage of Apoptotic MCF7 breast cancer cells post-bioprinting. The annexin A5-FITC kit was used for this analysis. Cells were collected post-bioprinting at 3 time frames, 2, 24 h, and 7 days Cells were BP in a modified HP thermal inkjet printer, onto a petri dish with EMEM, and the cells were incubated immediately. Graph (A) depicts the total percentage of apoptotic cell population is expressed as the sum of early and late apoptosis percentage (F2+F4), green bars. Dark bars depict necrotic cells, those are the cells stained with PI. Each bar represents the average of triplicates. Controls used were live cells stained, unstained, and dead cells (not shown). Error bars represent the standard deviations. (B) Dot plot graphs depict results at 2 h post-bioprinting. Graph (C) displays results at 24 h post-bioprinting. Graph (D) are the results from non-printed (manually seeded) MCF7s cells. P-values are from a two-tailed student t-test for independent tests (shown over the bars in A). 10,000 events were obtained per sample. The Kaluza Analysis v.1.3 software (Beckman Coulter) was used to extract the data and graphs.

Cell Morphology Results

Comparing the morphology of both, BP and MS MCF7 cells, no major changes were observed between the two cell samples. Additionally, the fluorescence intensity of the Neu antibody was measured to quantify protein expression on both cells samples, the fluorescence intensity measurements were significantly higher by 3-fold in MS cells than in the BP cells (p < 0.001). Fluorescent parameter settings in that channel, Alexa 568, was the same for both samples. Visually, fluorescence intensity in MS cell samples was markedly more evident than BP cells as seen in Figure 2.

FIGURE 2
www.frontiersin.org

Figure 2. (A) Bioprinted MCF7 cells 24 h post-BP. Image showing in all three channels. (B) Manually Seeded MCF7 fixed at 24 h post seeding. (C) Bioprinted MCF7 Cells and (D) manually seeded MCF7 cells, single channel image, stained with Neu (sc-33684) primary antibody and goat IgG anti-mouse secondary conjugated with Alexa 568 channel. The same parameter settings were used for both cell samples. Fluorescent Intensity measurements were significantly different between the two samples, cytosol and nucleus were measured separately. In (E) Mean intensity measurements for the cytosol and the nucleus of the BP cells were 28.9 (1.6) and 24.3 (1.7), respectively. The mean intensity measurements for the cytosol and nucleus of MS cells were 87.0 (7.4) and 183.6 (9.8), respectively (BP Cyto = BP Cytosol, MS Cyto = Manually Seeded Cytosol, MS Nuc = Manually Seeded Nucleus, BP Nuc = BP Nucleus).

Phospho-MAPK Array Results

The BP MCF7 cells phospho-MAPK array shown in Figure 3 revealed a total of 21 highly phosphorylated sites, whereas the MS cells revealed 11 phosphorylated sites. Of the phosphorylated targets in the MS cell samples, 6 analytes displayed 1.6-fold stronger phosphorylation levels than the BP cells (targets: p70S6 Kinase, CREB, ERK2, ERK1, Akt1, and Akt2), though no statistical significance was observed (p = 0.272). On the other hand, the BP MCF7s displayed 21 phosphorylated sites of which 12 did not show in the MS cells. These sites were: RSK1, HSP27, p38δ, p38β, MSK2, p53, MKK6, TOR, MKK3, p38γ, RSK2, and JNK2. Phosphorylated targets above the threshold were further investigated to associate targets with key biological processes, molecular functions, or the cellular components in the cells and to correlate with RNA sequencing results. In the network analysis shown in Figures 4, 5, a total of 24 targets were observed from the phosphorylated kinases in BP cells whereas for MS cells, 6 phosphorylated sites were observed in this network. This analysis showed MAPK1, CREB1, TP53, MAPK3, MAPK8, AKT1, HSPB1, MAP2K3, and AKT2 with >10 direct edges in the BP cell sample group.

FIGURE 3
www.frontiersin.org

Figure 3. Activation of cellular kinases by thermal inkjet bioprinting. Chemiluminescent images in iBright FL1000 of a Proteome Profiler Human Phospho-MAPK Array (Catalog # ARY002B). (A) Membranes of manually seeded and BP MCF7 breast cancer cells. Membrane arrangement: LEFT = manually seeded (MS) MCF7 cells, RIGHT = BP MCF7. Signal for each kinase is represented by a pair of duplicate spots; three reference pairs are shown in three upper/lower corners. (B) Histogram profiles for selected analytes were generated by quantifying the mean spot pixel density exposure in the iBright FL1000 and Invitrogen™ iBright™ Analysis Software v.3.0. Kinases that show increased levels of phosphorylation are identified. Mean pixel density for the analytes is shown in the bar graphs. Twenty one kinases appeared phosphorylated in the BP samples, whereas 10 kinases showed in the manually seeded. Of the MS cells 6 kinases were strongly phosphorylated by >1.6-fold as compared to the BP cells.

FIGURE 4
www.frontiersin.org

Figure 4. Network of analytes phosphorylated in BP MCF7 breast cancer cells. Functions selected in this network were regulator of apoptosis (green, 8), response to stress (red, 17), intracellular signal transduction (yellow, 21), and signal regulators (blue, 16). This network depicts functional interactions among BP BC predisposed genes. In this network of phosphorylated sites, there are significantly more interactions than expected (p ≤ 0.001).

FIGURE 5
www.frontiersin.org

Figure 5. (A) Schematic summary of the intracellular pathways activated by bioprinting MCF7 cells. Kinases in green ovals not included in the phosphorylated kinases but they represent critical kinases, pink ovals represent overlapping kinases, present in both, manually seeded, and bioprinted cells. Blue ovals represent kinases phosphorylated only in BP cells. Red connector lines ending in a perpendicular line mean the originating kinase blocking signaling activity or deactivation. Representation of the critical kinases and their pathways in the bioprinted cells. (B) Network with links classified based on curated evidence in STRING. Biological functions selected are p53 binding process (red, 4), regulation of phosphorylation (blue, 14), stress activated protein kinase signaling cascade (pink, 7), regulation of cellular response to heat (dark green, 4), and cellular response to stress (light green, 16). These nodes have a significant number of interactions, as expected due to phosphorylated targets were extracted. The interaction score was set at 0.7 with k-mean clustering set at 5, thus only links that have a high confidence probability are displayed (P < 0.001). A total of 24 targets were observed from the phosphorylated targets in Bioprinted cells whereas for manually seeded cells, 6 phosphorylated sites were observed in this network. Despite the complexity of this network, we observed that MAPK1, TP53, CREB1, MAPK3, MAPK8, AKT1, HSPB1, AKT2, and MAP2K3 proteins display more than 10 protein-protein interactions.

RNA Sequencing Results

RNA seq analysis conducted in BP MCF7 cells produced a total of 12,235 genes, of which 1,187 (9.7%) were significant (q ≤ 0.05). Using a cutoff value of a ±2-fold change there were 266 upregulated and 206 downregulated genes shown in Supplemental Information II. Furthermore, in the Log2 transformation of expression values, a total of 65 genes were significant (q ≤ 0.05), with five genes expressed only in the BP cells (NRN1L, LUCAT1, IL6, CCL26, and LOC401585). From the downregulated genes, every gene that showed in the MS cells also appeared in the BP cells. Of the five genes expressed in the BP cells, four gene ontologies (GO) in the molecular function were identified: protein binding, catalytic activity, molecular function regulator, and molecular transducer function (GO:0005488, GO:0098772, GO:0003824, and GO:0060089), four in the biological process: cellular response to stimulus, regulation, localization and processes (GO:0065007, GO:0009987, GO:0051179, and GO:0050896), and two in the cellular component ontology: the extracellular region and plasma membrane, which includes other external encapsulating systems such as the cell envelop (GO:0005623 and GO:0005576).

The following 4 genes, EFCAB11, FAM117B, FAM46B, and RBM43 were downregulated <-2-fold and the following 3 genes, GCLM, MOB3C, and VDR were upregulated >2-fold. No significant connections were found among those genes. EFCAB11 contains 11 nodes and 11 edges, though no statistical significance was observed, with an interaction score set at 0.4. Recent publications have associated this gene with hepatocellular carcinoma. Of the genes that were upregulated and have more than 10 protein to protein interactions are: CYP1A1, IL6, UGT1A6, EGFR, and CYP1B1 (see Table 1). Cytochrome P450 1A1 (CYP1A1) was upregulated by more than 300-fold and CYP1B1 by 10-fold in the BP cells. Cytochromes P450 are a group of enzymes involved in catalyzing endogenous substrates, including steroids, fatty acids, xenobiotics, and others (Spink et al., 1998). These CYP enzymes are accountable for DNA damage, which leads to tumor inception (Mikhailova et al., 2005). On the other hand, downregulated genes in BP cells with ≥10 protein interactions are: TP53, FOS, JUN, EGR1, HIST2H2AC, and FOSB (see Table 2). TP53, aka tumor suppressor gene, is further repressed in these cells; this explains the added resolute characteristics hence death evasion (Campbell et al., 2019). A Complete list of gene analysis is available in Supplemental Information II.

TABLE 1
www.frontiersin.org

Table 1. List of genes extracted from the upregulated genes (>2 fold) in the bioprinted MCF7 cells RNA seq analysis, this table contains only genes with ≥5 number of edges, network diagram depicted in Supplement Information I, Figure S2.

TABLE 2
www.frontiersin.org

Table 2. List of genes extracted from the downregulated genes (<-2 fold) in the bioprinted MCF7 cells RNA seq analysis, this table contains only genes with ≥5 number of edges, network diagram depicted in Supplement Information I, Figure S5.

Discussion

The present study evaluated morphology, viability, apoptosis, and necrosis, phospho-MAPK array, and RNA seq analysis of BP MCF7 BCCs. The cellular morphology of BP cells was compared to the MS cells. The BP cells appeared to have blebs, and non-spherical and distorted bulky features when compared to the MS cells (see images in Supplemental Information I). These protrusions may be caused by shear and mechanical stresses during bioprinting. Wickmann et al. identified blebbing as one of the distinctive characteristic of apoptosis (Wickman et al., 2013). This apoptotic morphological feature was confirmed by performing the Annexin A5-FITC necrosis/apoptosis assay, where after several hours following bioprinting, a number of cells appeared in the early and late apoptosis phase (Figure 1).

BP cells were interrogated with the neu antibody, neu is a glycoprotein expressed in ~30% of the cancers, yet not highly expressed in MCF7 cells. Neu (HER2) was probed in the BP and MS MCF7s to test the hypothesis that bioprinting would stimulate expression of the HER2/neu protein. However, results were contradicting, with MS cells showing significantly higher neu expression (Figure 2, p < 0.0001. This response is attributed to induced stress due to bioprinting, thus decreasing neu protein synthesis in the entire cell. We further compared the morphology of bioprinted and manually seeded cells and observed that the BP cells appeared smaller and flatter.

Viability, Apoptosis, Necrosis

Viability results were comparable to other published results which indicate a viability greater than 70% (Xu et al., 2013; Kolesky et al., 2016) but in some systems it can be as high as 90–98% (Cui et al., 2012b; Xu et al., 2013). Twenty-four hours post bioprinting (Ma et al., 2015). Chang et al., reported that BP HepG2 cells' viability in a multi-syringe nozzle ranged between 70 and 80% (Chang et al., 2008). Though, in that study, the bioink used to bioprint cells was sodium alginate crosslinked with CaCl2. In our experiment, the viability of the BP cells ranged from 60 to 87% tested at both, 24 and 48 h. Additionally, there were no media replenishments, which may have inflated viabilities in other studies. In brief, during the thermal BP process, cells are discharged from the modified print cartridge at ~1 m/s (Kim et al., 2016), they pass through a heated filament where the temperature is raised close to 300°C for a few microseconds (Calvert, 2001; Roth et al., 2004). Therefore, BP cells in a thermal inkjet printer are briefly heated as they pass through the printing filament before landing in a petri dish with pre-incubated cell media at a pH of 7.2 (Mironov et al., 2003; Roth et al., 2004; Singh et al., 2010). It is worth noting that during cell count preparation, most of the cells that did not survive bioprinting may have been discarded during the cell count preparation process because those cells were already detached from the petri dish, thus when centrifuged and supernatant was gently decanted some of the floating cells might have been removed, leaving mostly live cells in the pellet. Previous reports of BP cells have indicated viability rates ranging from 70 to 90% (Cui et al., 2012a; Peng et al., 2017).

Moroi et al. observed that tumor cells exposed to a 43°C water bath for 30 min quickly entered apoptosis depending on the localization of the tumor cells (Moroi et al., 1996). In addition, they observed that necrosis did not change immediately but it gradually increased within 3–6 h (Moroi et al., 1996). Measurements of results from the Annexin A5 analysis at 2 h post-bioprinting, resulted in an apoptotic rate of 30.5% with a necrotic rate of 0.1%. The apoptotic rate increased by 15% and the necrotic rate by 5.9% at 24 h following bioprinting. The percentage of live cells in this assay was 69 and 30% at 2 and 24 h, respectively, suggesting that bioprinting causes damage to the cells. On average, 12 and 19% of the BP cells evaluated at 2 and 24 h post-bioprinting appeared in the early apoptotic cycle. As sampling collection time was increased post-bioprinting, a 2.4-fold increase in the late apoptosis stage was observed 24 h later and a 1.15-fold increase 48 h. In another study, Catros et al. reported a mortality rate of endothelial cells embedded in sodium alginate solution BP with a laser-assisted bioprinter between 60 and 37%, which was dependent upon the thickness of the Matrigel™ bioink used and the laser energy applied, which varied between 20 to 100 μm and from 8 to 24 μJ, respectively (Catros et al., 2011). They observed a lower mortality rate in thicker gels that were 80–100 μm thick when higher laser energy was applied (Catros et al., 2011). However, the actual mortality rate was not directly measured, it was concluded from observed images. Bioprinting may be causing thermal and mechanical stress to the cells resulting in a high number of stressed and weak cells that when subjected to the Annexin A5-FITC assay, became apoptotic, instead of fully convalescing. This is important when considering drug tests, BP cells shall be allowed recovery time with at least one media change prior to testing any drugs (Campbell et al., 2019).

Phospho MAPK Array—Phosphorylated Analytes in BP MC7 BCCs

Glycogen synthase kinase-3 (GSK3 or GSK-3β), regulated by the Ras signaling pathway, stimulates cell proliferation (Milisav et al., 2017). GSK-3α, a Ser/Thr kinase, was identified as a deactivator of glycogen synthase and a regulator of other cellular functions, like cancer cell survival and proliferation (Ugolkov et al., 2016). In another arrangement, when GSK-3β is active, it phosphorylates β-catenin causing it to degrade (Brembeck et al., 2006; Weinberg, 2013). If however, GSK-3β is deactivated, β-catenin translocates to the nucleus where it binds Tcf/Lft Transcription Factors (TF) and stimulates expression of various genes related to cell proliferation amongst others (Weinberg, 2013). GSK-3β is deactivated via phosphorylation by the Akt/PKB pathway; this inhibition prevents cells from entering apoptosis (Brembeck et al., 2006; Weinberg, 2013). Inhibition of drug induced apoptosis may be triggered by GSK-3β and therefore, the bioprinted cell model may potentially mimic this process.

GSK-3 works as a multifunctional downstream switch that influences the output of several signaling pathways (Woodgett, 2005; Arcaro and Guerreiro, 2007) and dysregulated GSK-3 has been implicated in cancer. GSK-3α/β is highly phosphorylated in both cell sample conditions, yet in the BP cells, phosphorylation was 1.6-fold higher than in MS cells. GSK3α/β are encoded by GSK3A and GSK3B genes, respectively. In the RNA seq analysis, expression of this gene was upregulated in the BP cells by 1.2-fold (q = 0.1, 0.5) for both. Through the protein-protein network analysis we observed that this protein contains 9 interactions, and mutated GSK-3 can activate or de-activate other signaling pathways and collectively stimulate cancer progression.

Heat shock proteins (HSPs) have been identified to belong to a large family of stress response proteins and they function as molecular chaperones by assisting with the folding/unfolding of other cellular proteins and keeping substrate aggregation (Kato et al., 1994). HSPs are normally expressed at low levels but under heat shock they are highly phosphorylated (Rogalla et al., 1999; Almeida-Souza et al., 2010). Elevated levels of HSPs has been found in ischemia/reperfusion, cancer, and have been linked to cancer angiogenesis (Kato et al., 1994; Vos et al., 2009). HSP27 also functions as an anti-apoptotic molecule, regulating apoptosis through direct interaction of key components of the apoptotic pathway, this occurs in cooperation with HSP70 (Bryantsev et al., 2007; Vos et al., 2009). Phosphorylation for this analyte was observed in BP cells only, and it was not shown in MS cells. In the RNA seq analysis of the MCF7 cells, expression levels is higher than in MS MCF7 cells although no statistical significance was found (q = 0.37). Thermal inkjet bioprinting exposes cells to heat pulses during the bioprinting process, which would explain the high phosphorylation levels found in surviving BP cells (Mironov et al., 2003; Xu et al., 2013). This confirms one of our hypotheses that thermally bioprinted BCC stimulate HSPs, in this case HSP27. Increased levels of heat shock proteins may trigger a cell survival/resistance mechanism in cancer cells. Additionally, HSPs activate the angiogenic pathway that leads to vascular endothelial growth factor (VGEF) expression, vessel formation, and accelerated tumor formation in host tissues (Chatterjee et al., 2018; Plimpton et al., 2018; Solis et al., 2019).

Jun N-terminal Kinases (JNKs) are small proteins, 45–55 kDa, and product of three genes which, through different splicing, create up to 10 isoforms (Alvarez et al., 1991; Hibi et al., 1993). When activated, JNKs translocate to the nucleus where they regulate the process of various transcription factors, including the activator protein 1 (AP1) (Okazaki and Sagata, 1995; Weinberg, 2013). AP1 is implicated in antiestrogen resistance in breast cancer, which results in poor therapy response to hormonal therapy. This kinase has also been implicated in several biological processes such as cell migration, differentiation, proliferation, transformation, and apoptosis (Dai et al., 2000; Weinberg, 2013). It has been shown that environmental stresses, like heat shock and osmotic shock can activate this pathway (Kurokawa et al., 2000), which is suspected to be the case in BP MCF7s. Phosphorylation levels of JNK in BP cells was 4.4-fold higher than in the MS cells, even though expression levels of this gene changed only slightly.

Ribosomal S6 kinase (RSK1) is a broadly expressed component of the RSK consort of growth factor-regulated ser/thr kinases (Robinson and Cobb, 1997), which stimulate cell proliferation and differentiation (Hu et al., 2004). RSK1/2 appeared highly phosphorylated in BP cancer cells but not in MS cells (see Figure 3). RPS6KA1 is highly upregulated in BP MCF7 cells by 1.6-fold compared to MS cells, and a statistical significance was observed (q = 0.001). This suggests that MCF7 cells surviving the bioprinting process may become more resilient and this characteristic makes their use as in vitro models a great option for use during drug discovery.

P38 MAPKs are stimulated by proinflammatory cytokines, environmental stresses, hypoxia, and osmotic shock (Sayed et al., 2000; Tamura et al., 2000). Once stimulated, p38α phosphorylates several targets, which may act as tumor suppressors and the activation of those is important in cell proliferation and cell survival. Several compounds that inhibit p38α continue to be screened as potential anti-cancer therapies (Igea and Nebreda, 2015). P38α activation has been linked to tumor cell survival of chemotherapeutic treatments. In BP MCF7 cells, p38 was highly phosphorylated whereas in MS cells no phosphorylation was observed, while the expression levels in RNA seq were not significantly higher. The activation of p38 by bioprinting is a process that shows promise when used in tumor models as an anti p38 drug candidate.

MSK2 is activated in stress-related signaling by P38a/MAPK (Pierrat et al., 1998) and is phosphorylated in BP cells but inactivated in MS cells. Cancer cells are known to metabolize large amounts of their glucose through the glycolysis cycle, rather than the citric cycle meaning that activation of MSK2 in the BP cells may assist with regulation of glycogen metabolism that promotes cell survival.

Phosphorylated Sites in Both Manually Seeded and Bioprinted MCF7 BCCs

Seventy kilodalton ribosomal protein S6 kinase (p70S6) controls the activation of cell growth and is generally expressed in adult human tissues (Ferrari and Thomas, 1994; Cross et al., 1995). This activity increases during the G1 phase by 20-fold when released from the G0 phase of the cell cycle (Edelmann et al., 1996; Knowlton et al., 2011). P70S6K is highly phosphorylated in MS cells; however, phosphorylation levels decrease by 0.6-fold in BP MCF7 cells which suggest cell cycle arrest at G0. This explains low cell survivability levels post-bioprinting.

The cyclic AMP Response Element-Binding Protein (CREB) is part of the bZIP family of TFs. CREB is associated with cancer growth and poor clinical outcomes in several types of cancer, including breast cancer (Son et al., 2010). In MS cells, CREB and extracellular signal-regulated kinases 1/2 (ERK1/2) were constitutively phosphorylated by 2 and 3.5-fold, respectively, which is higher than in BP MCF7 BCCs but the gene expression levels do not differ significantly. The lower phosphorylation levels in the BP cells suggests that CREB may be regulating cell cycle progression in the MS cells yet not the BP cells; however, this could mean that CREB may be inactivated temporarily in the BP cells yet when these cells are allowed to recuperate from the bioprinting process, CREB levels may rise again to levels typical of cancer cells.

p53 is known as a tumor suppressor protein, and it has been found mutated in many types of tumors (Bond et al., 2005; Vaseva et al., 2012). A number of events can stimulate increased p53 levels, such as exposure to acidic environments or to nitrous oxide, low oxygen tension, lack of an intracellular group of nucleotides, blockage of DNA or RNA synthesis and other insults that may occur in the cells (Ghosh et al., 2004; Vaseva et al., 2012). P53 phosphorylation levels were high in the BP cells however, it was completely deactivated in the MS cells. This suggests that p53 is highly dysregulated in these cells, yet it may also suggest that p53 could have been activated through the JNK signaling pathway as described above. However, in the RNA seq expression analysis, P53 were two-fold higher (q = 0.001) in the MS cells than in the BP cells, suggesting that p53 is further suppressed in the BP cancer cells, though the phosphorylation of p53 remains unknown at this point. It may be possible that p53 is deactivated in BP cells indirectly via the Akt signaling cascade.

Mitogen Activated Protein Kinase Kinase 6 (AP2K6 or MKK6) has a critical role in the MAPK signal transduction pathway (Raingeaud et al., 1996), as it can activate p38 (Moriguchi et al., 1996; Raingeaud et al., 1996; Stein et al., 1996) in which then leads toward the specific activation of transcription factors, such as activating transcription factor 2 (ATF2) and ETS-like 1 (Elk1) (Goedert et al., 1997). ATF2 dysfunction has been linked to cancer metastasis (Watson et al., 2017) These proteins are activated by several insults, such as heat shock, UV rays, osmotic shock, cellular DNA damage and their regulation remains obscure (Goedert et al., 1997). The phosphorylation level of MKK6/MKK3 was observed in the BP cells yet it was not seen in the MS MCF7 cells. The MAP2K3 gene showed no statistical significance in expressions but MAP2K6 expression levels were significantly lower (0.35-fold) than in MS cells (q = 0.027). Hyper-phosphorylation of these kinases in BP cells is attributed to the mechanical, shear, and heat traumas triggered by bioprinting; this feature may be beneficial for use in tumor models due to the high number of phosphorylated sites, which could mimic in vivo stresses.

Akts participate in cell survival and proliferation, cycle progression, vesicle trafficking, glucose transport, metabolism, and biological processes (Noguchi et al., 2014; Ruiz-Medina et al., 2019). Akt's activation is believed to have a key function in tumorigenesis (Martini et al., 2014). Akt1 and Akt2 appeared phosphorylated in MS MCF7 cells. Phosphorylation of Akt1 is not observed in the BP cells. However, phosphorylation of Akt2 was 1.7-fold higher in the MS BCCs. Lower phosphorylation levels in BP cells suggest that its role as glucose uptake regulator is limited temporarily due to the effect of bioprinting. Akts, aka proto oncoproteins, have several downstream effects when activated, two of which are the binding of B-cell lymphoma 2, BCL2-associated X (BAX) protein, producing an anti-apoptotic effect when BAX blocks the mitochondria. Another downstream effect is protein synthesis via Rheb, resulting in the activation of mTOR thus interacting and activating S6K to initiate mRNA translation (Tabatabaian et al., 2010).

The extracellular signal-regulated kinases 1/2 (ERK1/2), once phosphorylated, can regulate several cellular processes such as protein synthesis, chromatin remodeling, and transcription (Weinberg, 2013). Mutations in proteins along the ERK pathway are common in cancers. In MCF7 cells, the ERK pathway is activated by estradiol (Robinson and Cobb, 1997). While ERK1/2s were highly phosphorylated in MS cells, the phosphorylation levels in BP cells is lower by 0.3 and 0.5-fold and the RNA expression levels for both ERK1 and ERK2 are not significantly different. Lower phosphorylation levels in the BP cells indicate some signaling is still evident and this may be a transitional event for ERK1/2 can be activated via the Akt/MEK pathway.

Mammalian Target of Rapamycin (mTOR) a ser/thr protein kinase acts as the main regulator related to energy and stress signaling, cellular metabolism, growth factors and nutrients, growth and survival in response to hormones, autophagy and cytoskeletal re-organization (Kim et al., 2002; Park et al., 2002; Inoki et al., 2003; Brugarolas et al., 2004). mTOR was found to prevent autophagy through phosphorylation both directly through the phosphorylation of ULK12 and indirectly through DAP1 (Koren et al., 2010). This analyte, is phosphorylated in the BP MCF7 cells but not in the MS MCF7 cells, which suggests that newly bioprinted cells may be inducing cell proliferation. Expression levels for the TOR1A and mTOR were not statistically significant. Getting to know the phosphorylated pathways of bioprinted cancer cells at the molecular level is critical when considering this process in order to use these tumor models in drug discovery.

RNA Sequence Analysis

For decades, preclinical and clinical evidence have shown that cancer cells exhibit significantly higher sensitivity to hyperthermia than normal cells. Intratumoral cell temperatures range from 42 to 45°C (Habash et al., 2006). A cellular exposure of 45°C for 30 min specifically, altered the expression of key mitotic regulators and altered G2/M phase progression in breast cancer cells, including the MCF7 cell line (Amaya et al., 2014). Thus through RNA seq analysis, we sought to explore the differences in gene expression of the manually seeded and bioprinted MCF7 cells as the bioprinting process further exposes cells to high levels of heat.

The non-protein coding gene Lung Cancer Associated Transcript 1 (LUCAT1) was among the genes found differentially and significantly overexpressed in the BP cells as compared to the MS MCF7 cells; yet expression levels found in breast cancer stem cells (BCSCs) is higher than in normal BCCs (Zheng et al., 2019). Upregulated LUCAT1 has been implicated with BCC proliferation and shorter overall survival and progression-free survival (Yoon et al., 2018). LUCAT1 is a member of the long intervening noncoding RNAs (lincRNAs) class that seems to function downstream on the nuclear factor erythroid 2–related factor 2 (NRF2) (Tonelli et al., 2018). NRF2 controls gene expression that regulates oxidative stress protection in airway epithelial cells (Thai et al., 2013). Thus, it is likely that the LUCAT1's overexpression is an important contribution to the augmented resistance to death of the BP MCF7 cells via NRF2, which is a developing regulator of cellular resistance to oxidants (Ma, 2013). LUCAT1, expressed only in BP MCF7 BCCs, suggests that bioprinting may be differentially selecting more specifically BCSCs. It is also possible that the genes associated with stem cell phenotype are already expressed in few of the cells beforehand and when the cells were bioprinted, it induced cell proliferation of more cells with the stem cell phenotype. LUCAT1 has been identified as a potential target for drug discovery (Zheng et al., 2019).

The C-C motif Chemokine Ligand 26 (CCL26) was also found to be overexpressed in the BP MCF7 cells, as compared to the MS cells. CCL26, also known as Eotaxin-3, is a chemotactic cytokine for eosinophils and basophils and exerts its effect by binding to its receptor CCR3 (Butterfield et al., 1993; Kitaura et al., 1996; Shinkai et al., 1999). It was reported that CCL26 is regulating expression of cancer-associated genes during airway inflammation (Marijani et al., 2011). Furthermore, the upregulation of this gene, as well as CCL2, IL6, and LOXL2, has been implicated as part of the effects of cancer-associated fibroblasts in promoting progression of hepatocellular carcinoma cells (Lin et al., 2012). Thus, the CCL26-CCR3 ligand-receptor structure appears to be involved in inflammatory processes and it is upregulated in cancer progression, which highlights the suitability as an in vitro model for anti-cancer drug discovery initiatives.

Of the genes that were downregulated in the BP MCF7s, in the molecular function ontology EF-hand calcium-binding domain-containing protein 11 (EFCAB11) is a protein coding gene and it has been implicated in hepatocellular carcinoma, and FAM117B (under biological process) has been found to negatively regulate activation-induced death in T-cells and to respond to osmotic stress. Upregulated, the glutamate-cysteine ligase modifier (GCLM) subunit is linked to apoptotic mitochondrial changes and it is a negative regulator of the extrinsic signaling pathway. This pathway utilizes a ligand binding process to induce apoptosis, which initiates in the extracellular milieu (Weinberg, 2013).

Under normal circumstances, cancer cells can bypass the in vivo removal tactic of the immune system by implementing diverse strategies such as evasion of apoptosis (Kobayashi et al., 1993). Those cellular eluding mechanisms include mutations that provide the ability to counteract the immune system, which fails to detect and eliminate. Cancer cells acquire the ability to undergo rapid proliferation without regard to the normal homeostatic system; uncontrolled proliferation is a distinctive characteristic of cancer development and progression (metastasis) (Feitelson et al., 2015). Furthermore, RNA seq analysis performed in bioprinted MCF7 cells provided significant information by identifying genes involved in response to the bioprinting process. These data clearly demonstrated that thermal inkjet bioprinting is triggering significant number of gene alterations that could potentially be used for target drug discovery.

We believe that bioprinting may stimulate other conditions like exacerbating pathways implicated in drug immunity, cell motility, proliferation, survival, and differentiation such as the expression of NRN1L, LUCAT1, IL6, and CCL26, which have been implicated with numerous diseases. This underscores the need to use these thermal inkjet bioprinting tumor models to simulate in vivo conditions where, in some cases, tumors develop immunity to chemotherapeutic drugs. Additional to these mutations, we have observed that cancer cells surviving the bioprinting process, elicit further phenotypic changes, which we believe makes them even more resistant to cancer drugs. Thus we are proposing to use bioprinted cancer cells, which become resilient and hard to kill, to discover a new class of more efficacious and potent cancer drugs. Collectively, these findings contribute with relevant information to affirm that this BP model can mimic a novel tumor model that should be studied further and applied in preclinical studies.

Conclusion

Bioprinted MCF7 cells showed increased levels of phosphorylation in analytes that have been identified as key players in activating critical pathways that when dysregulated, are associated with biological aggressive oncogenic properties. Results of bioprinted MCF7 BCC at the physiological and molecular levels were evaluated at different time frames following bioprinting in a modified thermal inkjet printer. Despite losing up to 37% of the cell population 24 h post bioprinting, it is evident that this process induces ample phosphorylation at several critical sites. Several key signaling cascades were activated, as was observed in the phosphorylation of MKK6/MKK3 and RSK1/2, which have been implicated with cancer metastasis. Collectively, GSK-3α/β, ERK1, Akt1/2, JNK, RSK1/2, HSP27, p38, MSK2, p53, MKK3/6, and TOR when mutated or hyper-phosphorylated are implicated in biologically aggressive behaviors which could be used to develop an in vitro tumor model that could be used to explore drug discovery. Unregulated levels of heat shock proteins may also be triggering a cell survival mechanism in cancer cells like a conglomeration of receptors in the ECM. In this experiment, bioprinted MCF7 cells showed increased levels of chaperone protein HSP27. It is possible that bioprinting may be stimulating an increased amount of ligand-independent receptors in the cell surface. This occurrence may cause receptor collisions which lead to receptor dimerization thereby making them more sensitive in overexpressed genes associated with cellular resistance to intrinsic biological processes and cellular functions.

The RNA seq analysis of BP MCF7 cells indicates that BP MCF7 cells may be causing several mutations which enable cells to become more robust encouraging MCF7 cells to activate key kinases implicated in cancer development, proliferation, and metastasis. Furthermore, through the phospho-MAPK array, we confirmed that bioprinting activates signaling pathways associated with cellular response to apoptosis, mitosis, cell migration, transcription, and other cellular functions. Additionally, it may also be triggering other undiscovered functions in the BP BCCs, such as activation of critical pathways implicated in drug immunity or inducing cell motility, proliferation, survival, and differentiation. Moreover, these data suggest that thermal inkjet bioprinting is stimulating large scale gene alterations that could potentially be used with autologous drug tests to confirm drug efficacy prior to initiating cancer therapy. Insights into the cell response after bioprinting have demonstrated that BP cells can potentially improve the in vitro models for drug discovery.

Data Availability Statement

The datasets generated for this study are available on request to the corresponding author.

Author Contributions

AC, AV-R, and TB conceived and designed the study. AC, JM, and DG collected the data. AC, JM, AV-R, and TB analyzed and interpreted the data. AC, JM, AV-R, and DG wrote the manuscript. TB wrote parts of the manuscript and edited it to the final version.

Funding

This study was supported by the National Institute of General Medical Sciences of the National Institutes of Health under linked Award Numbers RL5GM118969, TL4GM118971, and UL1GM118970, and NIH-NIMHD-RCMI Grant No. 5G12MD007592 and Grant No. 2U54MD007592 from National Institute on Minority Health and Health Disparities (NIMHD), a component of the National Institutes of Health (NIH).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

The authors wish to acknowledge Clarissa Amaya and Ana Betancourt for their efforts with manuscript preparation. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fbioe.2020.00082/full#supplementary-material

References

Almeida-Souza, L., Goethals, S., de Winter, V., Dierick, I., Gallardo, R., Van Durme, J., et al. (2010). Increased monomerization of mutant HSPB1 leads to protein hyperactivity in Charcot-Marie-Tooth neuropathy. J. Biol. Chem. 285, 12778–12786. doi: 10.1074/jbc.M109.082644

PubMed Abstract | CrossRef Full Text | Google Scholar

Alvarez, E., Northwood, I. C., Gonzalez, F. A., Latour, D. A., Seth, A., Abate, C., et al. (1991). Pro-Leu-Ser/Thr-Pro is a consensus primary sequence for substrate protein phosphorylation. Characterization of the phosphorylation of c-myc and c-jun proteins by an epidermal growth factor receptor threonine 669 protein kinase. J. Biol. Chem. 266, 15277–15285

PubMed Abstract | Google Scholar

Amaya, C., Kurisetty, V., Stiles, J., Nyakeriga, A. M., Arumugam, A., Lakshmanaswamy, R., et al. (2014). A genomics approach to identify susceptibilities of breast cancer cells to “fever-range” hyperthermia. BMC Cancer 14:81. doi: 10.1186/1471-2407-14-81

PubMed Abstract | CrossRef Full Text | Google Scholar

Andrade, E. L., Bento, A. F., Cavalli, J., Oliveira, S. K., Freitas, C. S., Marcon, R., et al. (2016). Non-clinical studies required for new drug development-Part I: early in silico and in vitro studies, new target discovery and validation, proof of principles and robustness of animal studies. Brazil. J. Med. Biol. Res. 49:e5644. doi: 10.1590/1414-431x20165644

PubMed Abstract | CrossRef Full Text | Google Scholar

Arcaro, A., and Guerreiro, A. S. (2007). The phosphoinositide 3-kinase pathway in human cancer: genetic alterations and therapeutic implications. Curr. Genom. 8, 271–306. doi: 10.2174/138920207782446160

PubMed Abstract | CrossRef Full Text | Google Scholar

Begley, C. G., and Ellis, L. M. (2012). Drug development: raise standards for preclinical cancer research. Nature 483:531. doi: 10.1038/483531a

PubMed Abstract | CrossRef Full Text | Google Scholar

Bolger Anthony, M, Lohse, M., and Usadel, B. (2014). Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 30, 2114–2120. doi: 10.1093/bioinformatics/btu170

PubMed Abstract | CrossRef Full Text | Google Scholar

Bond, G. L., Hu, W., and Levine, A. J. (2005). MDM2 is a central node in the p53 pathway: 12 years and counting. Curr. Cancer Drug Targets 5, 3–8. doi: 10.2174/1568009053332627

PubMed Abstract | CrossRef Full Text | Google Scholar

Brembeck, F. H., Rosário, M., and Birchmeier, W. (2006). Balancing cell adhesion and Wnt signaling, the key role of β-catenin. Curr. Opin. Genet. Dev. 16, 51–59. doi: 10.1016/j.gde.2005.12.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Brugarolas, J., Lei, K., Hurley, R. L., Manning, B. D., Reiling, J. H., Hafen, E., et al. (2004). Regulation of mTOR function in response to hypoxia by REDD1 and the TSC1/TSC2 tumor suppressor complex. Genes Dev. 18, 2893–2904. doi: 10.1101/gad.1256804

PubMed Abstract | CrossRef Full Text | Google Scholar

Bryantsev, A. L., Kurchashova, S. Y., Golyshev, S. A., Polyakov, V. Y., Wunderink, H. F., Kanon, B., et al. (2007). Regulation of stress-induced intracellular sorting and chaperone function of Hsp27 (HspB1) in mammalian cells. Biochem. J. 407, 407–417. doi: 10.1042/BJ20070195

PubMed Abstract | CrossRef Full Text | Google Scholar

Butterfield, J., Leiferman, K., and Gleich, G. (1993). Eosinophil-associated diseases. Immunopharmacol. Eosinophils 152–92

PubMed Abstract

Calvert, P. (2001). Inkjet printing for materials and devices. Chem. Mater. 13, 3299–3305. doi: 10.1021/cm0101632

CrossRef Full Text | Google Scholar

Campbell, A., Philipovskiy, A., Heydarian, R., Varela-Ramirez, A., Gutierrez, D. A., Solis, L. H., et al. (2019). 2D and 3D thermally bioprinted human MCF-7 breast cancer cells: a promising model for drug discovery. J. Clin. Oncol. 37, 2605–2605. doi: 10.1200/JCO.2019.37.15_suppl.2605

CrossRef Full Text | Google Scholar

Catros, S., Guillotin, B., Bačáková, M., Fricain, J.-C., and Guillemot, F. (2011). Effect of laser energy, substrate film thickness and bioink viscosity on viability of endothelial cells printed by Laser-Assisted Bioprinting. Appl. Surf. Sci. 257, 5142–5147. doi: 10.1016/j.apsusc.2010.11.049

CrossRef Full Text | Google Scholar

Chang, R., Nam, J., and Sun, W. (2008). Direct cell writing of 3D microorgan for in vitro pharmacokinetic model. Tissue Eng. Part C Methods 14, 157–166. doi: 10.1089/ten.tec.2007.0392

PubMed Abstract | CrossRef Full Text | Google Scholar

Chatterjee, B. K, Puri, S, Sharma, A, Pastor, A, and Chaudhuri, T. K. (2018). “Molecular chaperones: structure-function relationship and their role in protein folding,” in Regulation of Heat Shock Protein Responses. Heat Shock Proteins, Vol. 13, eds A. Asea and P. Kaur (Cham: Springer). doi: 10.1007/978-3-319-74715-6_8

CrossRef Full Text | Google Scholar

Cross, D. A., Alessi, D. R., Cohen, P., Andjelkovich, M., and Hemmings, B. A. (1995). Inhibition of glycogen synthase kinase-3 by insulin mediated by protein kinase B. Nature 378:785. doi: 10.1038/378785a0

PubMed Abstract | CrossRef Full Text | Google Scholar

Cui, X., Breitenkamp, K., Finn, M., Lotz, M., and D'Lima, D. D. (2012a). Direct human cartilage repair using three-dimensional bioprinting technology. Tissue Eng. Part A 18, 1304–1312. doi: 10.1089/ten.tea.2011.0543

PubMed Abstract | CrossRef Full Text | Google Scholar

Cui, X., Breitenkamp, K., Lotz, M., and D'Lima, D. (2012b). Synergistic action of fibroblast growth factor-2 and transforming growth factor-beta1 enhances bioprinted human neocartilage formation. Biotechnol. Bioeng. 109, 2357–2368. doi: 10.1002/bit.24488

PubMed Abstract | CrossRef Full Text | Google Scholar

Dai, R., Frejtag, W., He, B., Zhang, Y., and Mivechi, N. F. (2000). c-Jun NH2-terminal kinase targeting and phosphorylation of heat shock factor-1 suppress its transcriptional activity. J. Biol. Chem. 275, 18210–18218. doi: 10.1074/jbc.M000958200

PubMed Abstract | CrossRef Full Text | Google Scholar

De Maria, C., Rincon, J., Duarte, A. A., Vozzi, G., and Boland, T. (2013). A new approach to fabricate agarose microstructures. Polym. Adv. Technol. 24, 895–902. doi: 10.1002/pat.3162

CrossRef Full Text | Google Scholar

Edelmann, H. M., Kühne, C., Petritsch, C., and Ballou, L. M. (1996). Cell cycle regulation of p70 S6 kinase and p42/p44 mitogen-activated protein kinases in Swiss mouse 3T3 fibroblasts. J. Biol. Chem. 271, 963–971. doi: 10.1074/jbc.271.2.963

PubMed Abstract | CrossRef Full Text | Google Scholar

Feitelson, M. A., Arzumanyan, A., Kulathinal, R. J., Blain, S. W., Holcombe, R. F., Mahajna, J., et al. (2015). “Sustained proliferation in cancer: mechanisms and novel therapeutic targets,” in Seminars in Cancer Biology, Vol. 35 (Academic Press), S25–S54. doi: 10.1016/j.semcancer.2015.02.006

CrossRef Full Text | Google Scholar

Ferrari, S., and Thomas, G. (1994). S6 Phosphorylation and the p70s6k/p85s6k. Crit. Rev. Biochem. Mol. Biol. 29, 385–413. doi: 10.3109/10409239409083485

PubMed Abstract | CrossRef Full Text | Google Scholar

Ghosh, A., Stewart, D., and Matlashewski, G. (2004). Regulation of human p53 activity and cell localization by alternative splicing. Mol. Cell. Biol. 24, 7987–7997. doi: 10.1128/MCB.24.18.7987-7997.2004

PubMed Abstract | CrossRef Full Text | Google Scholar

Goedert, M., Cuenda, A., Craxton, M., Jakes, R., and Cohen, P. (1997). Activation of the novel stress-activated protein kinase SAPK4 by cytokines and cellular stresses is mediated by SKK3 (MKK6); comparison of its substrate specificity with that of other SAP kinases. EMBO J. 16, 3563–3571. doi: 10.1093/emboj/16.12.3563

PubMed Abstract | CrossRef Full Text | Google Scholar

Gudapati, H., Dey, M., and Ozbolat, I. (2016). A comprehensive review on droplet-based bioprinting: past, present and future. Biomaterials 102, 20–42. doi: 10.1016/j.biomaterials.2016.06.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Habash, R. W. Y, Bansal, R., Krewski, D., and Alhafid, H. T. (2006). Thermal therapy, part 2: hyperthermia techniques. Crit. Rev. Biomed. Eng. 34, 491–542. doi: 10.1615/critrevbiomedeng.v34.i6.30

PubMed Abstract | CrossRef Full Text | Google Scholar

Hevener, K. E. (2018). “Computational Toxicology Methods in Chemical Library Design and High-Throughput Screening Hit Validation,” in Computational Toxicology. Methods in Molecular Biology, Vol. 1800, ed O. Nicolotti (New York, NY: Humana Press). doi: 10.1007/978-1-4939-7899-1_13

PubMed Abstract | CrossRef Full Text | Google Scholar

Hibi, M., Lin, A., Smeal, T., Minden, A., and Karin, M. (1993). Identification of an oncoprotein-and UV-responsive protein kinase that binds and potentiates the c-Jun activation domain. Genes Dev. 7, 2135–2148. doi: 10.1101/gad.7.11.2135

PubMed Abstract | CrossRef Full Text | Google Scholar

Hu, Y., Fang, X., Dunham, S. M., Prada, C., Stachowiak, E. K., and Stachowiak, M. K. (2004). 90-kDa ribosomal S6 kinase is a direct target for the nuclear fibroblast growth factor receptor 1 (FGFR1) ROLE IN FGFR1 SIGNALING. J. Biol. Chem. 279, 29325–29335. doi: 10.1074/jbc.M311144200

PubMed Abstract | CrossRef Full Text | Google Scholar

Igea, A., and Nebreda, A. R. (2015). The stress kinase p38α as a target for cancer therapy. Cancer Res. 75, 3997–4002. doi: 10.1158/0008-5472.CAN-15-0173

PubMed Abstract | CrossRef Full Text | Google Scholar

Iglesias-Figueroa, B. F., Siqueiros-Cendón, T. S., Gutierrez, D. A., Aguilera, R. J., Espinoza-Sánchez, E. A., Arévalo-Gallegos, S., et al. (2019). Recombinant human lactoferrin induces apoptosis, disruption of F-actin structure and cell cycle arrest with selective cytotoxicity on human triple negative breast cancer cells. Apoptosis 24, 562–577. doi: 10.1007/s10495-019-01539-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Inoki, K., Zhu, T., and Guan, K. L. (2003). TSC2 mediates cellular energy response to control cell growth and survival. Cell 115, 577–590. doi: 10.1016/S0092-8674(03)00929-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Kato, K., Hasegawa, K., Goto, S., and Inaguma, Y. (1994). Dissociation as a result of phosphorylation of an aggregated form of the small stress protein, hsp27. J. Biol. Chem. 269, 11274–11278.

PubMed Abstract | Google Scholar

Kim, D. H., Sarbassov, D. D., Ali, S. M., King, J. E., Latek, R. R., Erdjument-Bromage, H., et al. (2002). mTOR interacts with raptor to form a nutrient-sensitive complex that signals to the cell growth machinery. Cell 110, 163–175. doi: 10.1016/S0092-8674(02)00808-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, Y. K., Park, J. A., Yoon, W. H., Kim, J., and Jung, S. (2016). Drop-on-demand inkjet-based cell printing with 30-μ m nozzle diameter for cell-level accuracy. Biomicrofluidics 10:064110. doi: 10.1063/1.4968845

CrossRef Full Text | Google Scholar

Kitaura, M., Nakajima, T., Imai, T., Harada, S., Combadiere, C., Tiffany, H. L., et al. (1996). Molecular cloning of human eotaxin, an eosinophil-selective CC chemokine, and identification of a specific eosinophil eotaxin receptor, CC chemokine receptor 3. J. Biol. Chem. 271, 7725–7730. doi: 10.1074/jbc.271.13.7725

PubMed Abstract | CrossRef Full Text | Google Scholar

Knowlton, D. L., Tang, K., Henstock, P. V., and Subramanian, R. R. (2011). miRNA alterations modify kinase activation in the IGF-1 pathway and correlate with colorectal cancer stage and progression in patients. J. Cancer 2:490. doi: 10.7150/jca.2.490

PubMed Abstract | CrossRef Full Text | Google Scholar

Kobayashi, R., Rassenti, L. Z., Meisenholder, G., Carson, D. A., and Kipps, T. J. (1993). Autoantigen inhibits apoptosis of a human B cell leukemia that produces pathogenic rheumatoid factor. J. Immunol. 151, 7273–7283.

PubMed Abstract | Google Scholar

Kolesky, D. B., Homan, K. A., Skylar-Scott, M. A., and Lewis, J. A. (2016). Three-dimensional bioprinting of thick vascularized tissues. Proc. Natl. Acad. Sci. U.S.A. 113, 3179–3184. doi: 10.1073/pnas.1521342113

PubMed Abstract | CrossRef Full Text | Google Scholar

Koren, I., Reem, E., and Kimchi, A. (2010). DAP1, a novel substrate of mTOR, negatively regulates autophagy. Curr. Biol. 20, 1093–1098. doi: 10.1016/j.cub.2010.04.041

PubMed Abstract | CrossRef Full Text | Google Scholar

Kurokawa, M., Mitani, K., Yamagata, T., Takahashi, T., Izutsu, K., Ogawa, S., et al. (2000). The evi-1 oncoprotein inhibits c-Jun N-terminal kinase and prevents stress-induced cell death. EMBO J. 19, 2958–2968. doi: 10.1093/emboj/19.12.2958

PubMed Abstract | CrossRef Full Text | Google Scholar

Langmead, B, and Salzberg, S. L. (2012). Fast gapped-read alignment with Bowtie 2. Nat Methods. 9:357. doi: 10.1038/nmeth.1923

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, Z.-Y., Chuang, Y.-H., and Chuang, W.-L. (2012). Cancer-associated fibroblasts up-regulate CCL2, CCL26, IL6 and LOXL2 genes related to promotion of cancer progression in hepatocellular carcinoma cells. Biomed. Pharmacother. 66, 525–529. doi: 10.1016/j.biopha.2012.02.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Ma, Q. (2013). Role of nrf2 in oxidative stress and toxicity. Annu. Rev. Pharmacol. Toxicol. 53, 401–426. doi: 10.1146/annurev-pharmtox-011112-140320

PubMed Abstract | CrossRef Full Text | Google Scholar

Ma, Y., Ji, Y., Huang, G., Ling, K., Zhang, X., and Xu, F. (2015). Bioprinting 3D cell-laden hydrogel microarray for screening human periodontal ligament stem cell response to extracellular matrix. Biofabrication 7:044105. doi: 10.1088/1758-5090/7/4/044105

PubMed Abstract | CrossRef Full Text | Google Scholar

Marijani, R., Mathews, T., McGann, H., and Abonyo, B. (2011). Regulatory effect of CCL26 on cancer related genes during inflammation. The FASEB Journal. 25, 561–7.

Google Scholar

Martini, M., De Santis, M. C., Braccini, L., Gulluni, F., and Hirsch, E. (2014). PI3K/AKT signaling pathway and cancer: an updated review. Ann. Med. 46, 372–383. doi: 10.3109/07853890.2014.912836

PubMed Abstract | CrossRef Full Text | Google Scholar

Mikhailova, O. N., Gulyaeva, L. F., and Filipenko, M. L. (2005). Gene expression of drug metabolizing enzymes in adult and aged mouse liver: a modulation by immobilization stress. Toxicology. 210, 189–196. doi: 10.1016/j.tox.2005.01.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Milisav, I., Poljšak, B., and Ribarič, S. (2017). Reduced risk of apoptosis: mechanisms of stress responses. Apoptosis 22, 265–283. doi: 10.1007/s10495-016-1317-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Miri, A. K., Khalilpour, A., Cecen, B., Maharjan, S., Shin, S. R., and Khademhosseini, A. (2019). Multiscale bioprinting of vascularized models. Biomaterials 198, 204–216. doi: 10.1016/j.biomaterials.2018.08.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Mironov, V., Boland, T., Trusk, T., Forgacs, G., and Markwald, R. R. (2003). Organ printing: computer-aided jet-based 3D tissue engineering. Trends Biotechnol. 21, 157–161. doi: 10.1016/S0167-7799(03)00033-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Moriguchi, T., Kuroyanagi, N., Yamaguchi, K., Gotoh, Y., Irie, K., Kano, T., et al. (1996). A novel kinase cascade mediated by mitogen-activated protein kinase kinase 6 and MKK3. J. Biol. Chem. 271, 13675–13679. doi: 10.1074/jbc.271.23.13675

PubMed Abstract | CrossRef Full Text | Google Scholar

Moroi, J., Kashiwagi, S., Kim, S., Urakawa, M., Ito, H., and Yamaguchi, K. (1996). Regional differences in apoptosis in murine gliosarcoma (T9) induced by mild hyperthermia. Int. J. Hyperthermia 12, 345–354. doi: 10.3109/02656739609022523

PubMed Abstract | CrossRef Full Text | Google Scholar

Noguchi, M., Hirata, N., and Suizu, F. (2014). The links between AKT and two intracellular proteolytic cascades: ubiquitination and autophagy. Biochim. Biophys. Acta 1846, 342–352. doi: 10.1016/j.bbcan.2014.07.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Okazaki, K., and Sagata, N. (1995). The Mos/MAP kinase pathway stabilizes c-Fos by phosphorylation and augments its transforming activity in NIH 3T3 cells. EMBO J. 14, 5048–5059. doi: 10.1002/j.1460-2075.1995.tb00187.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Park, I. H., Bachmann, R., Shirazi, H., and Chen, J. (2002). Regulation of ribosomal S6 kinase 2 by mammalian target of rapamycin. J. Biol. Chem. 277, 31423–31429. doi: 10.1074/jbc.M204080200

PubMed Abstract | CrossRef Full Text | Google Scholar

Peng, W., Datta, P., Ayan, B., Ozbolat, V., Sosnoski, D., and Ozbolat, I. T. (2017). 3D bioprinting for drug discovery and development in pharmaceutics. Acta Biomater. 57, 26–46. doi: 10.1016/j.actbio.2017.05.025

PubMed Abstract | CrossRef Full Text | Google Scholar

Peng, W., Unutmaz, D., and Ozbolat, I. T. (2016). Bioprinting towards physiologically relevant tissue models for pharmaceutics. Trends Biotechnol. 34, 722–732. doi: 10.1016/j.tibtech.2016.05.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Phamduy, T. B., Sweat, R. S., Azimi, M. S., Burow, M. E., Murfee, W. L., and Chrisey, D. B. (2015). Printing cancer cells into intact microvascular networks: a model for investigating cancer cell dynamics during angiogenesis. Integr. Biol. 7, 1068–1078. doi: 10.1039/C5IB00151J

PubMed Abstract | CrossRef Full Text | Google Scholar

Pierrat, B., Correia, J. S., Mary, J. L., Tomás-Zuber, M., and Lesslauer, W. (1998). RSK-B, a novel ribosomal S6 kinase family member, is a CREB kinase under dominant control of p38α mitogen-activated protein kinase (p38αMAPK). J. Biol. Chem. 273, 29661–29671. doi: 10.1074/jbc.273.45.29661

PubMed Abstract | CrossRef Full Text | Google Scholar

Plimpton, R. L., Valpuesta, J. M., and Willardson, B. M. (2018). Mechanisms of protein folding by type II chaperonins. Role Mol. Chaperones Struct. Fold. Biol. Funct. Drug Interact. Client Proteins 1, 190–213. doi: 10.2174/9781681086156118010008

CrossRef Full Text | Google Scholar

Raingeaud, J., Whitmarsh, A. J., Barrett, T., Dérijard, B., and Davis, R. J. (1996). MKK3-and MKK6-regulated gene expression is mediated by the p38 mitogen-activated protein kinase signal transduction pathway. Mol. Cell. Biol. 16, 1247–1255. doi: 10.1128/MCB.16.3.1247

PubMed Abstract | CrossRef Full Text | Google Scholar

Robinson, M. J., and Cobb, M. H. (1997). Mitogen-activated protein kinase pathways. Curr. Opin. Cell Biol. 9, 180–186. doi: 10.1016/S0955-0674(97)80061-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Robles-Escajeda, E., Das, U., Ortega, N. M., Parra, K., Francia, G., Dimmock, J. R., et al. (2016). A novel curcumin-like dienone induces apoptosis in triple-negative breast cancer cells. Cell Oncol (Dordr). 39, 265–77. doi: 10.1007/s13402-016-0272-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Rogalla, T., Ehrnsperger, M., Preville, X., Kotlyarov, A., Lutsch, G., Ducasse, C., et al. (1999). Regulation of Hsp27 oligomerization, chaperone function, and protective activity against oxidative stress/tumor necrosis factor α by phosphorylation. J. Biol. Chem. 274, 18947–18956. doi: 10.1074/jbc.274.27.18947

PubMed Abstract | CrossRef Full Text | Google Scholar

Roth, E. A., Xu, T., Das, M., Gregory, C., Hickman, J. J., and Boland, T. (2004). Inkjet printing for high-throughput cell patterning. Biomaterials 25, 3707–3715. doi: 10.1016/j.biomaterials.2003.10.052

PubMed Abstract | CrossRef Full Text | Google Scholar

Ruiz-Medina, B. E., Lerma, D., Hwang, M., Ross, J. A., Skouta, R., Aguilera, R. J., et al. (2019). Green barley mitigates cytotoxicity in human lymphocytes undergoing aggressive oxidative stress, via activation of both the Lyn/PI3K/Akt and MAPK/ERK pathways. Sci. Rep. 9:6005. doi: 10.1038/s41598-019-42228-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Santiago-Vázquez, Y., Das, U., Varela-Ramirez, A., Baca, S. T., Ayala-Marin, Y., Lema, C., et al. (2016). Tumor-selective cytotoxicity of a novel pentadiene analogue on human leukemia/lymphoma cells. Clin. Cancer Drugs 3, 138–146. doi: 10.2174/2212697X03666160830165250

PubMed Abstract | CrossRef Full Text | Google Scholar

Sayed, M., Kim, S. O., Salh, B. S., Issinger, O. G., and Pelech, S. L. (2000). Stress-induced activation of protein kinase CK2 by direct interaction with p38 mitogen-activated protein kinase. J. Biol. Chem. 275, 16569–16573. doi: 10.1074/jbc.M000312200

PubMed Abstract | CrossRef Full Text | Google Scholar

Shinkai, A., Yoshisue, H., Koike, M., Shoji, E., Nakagawa, S., Saito, A., et al. (1999). A novel human CC chemokine, eotaxin-3, which is expressed in IL-4-stimulated vascular endothelial cells, exhibits potent activity toward eosinophils. J. Immunol. 163, 1602–1610

PubMed Abstract | Google Scholar

Singh, M., Haverinen, H. M., Dhagat, P., and Jabbour, G. E. (2010). Inkjet printing—process and its applications. Adv. Mater. 22, 673–685. doi: 10.1002/adma.200901141

PubMed Abstract | CrossRef Full Text | Google Scholar

Solis, L. H., Ayala, Y., Portillo, S., Varela-Ramirez, A., Aguilera, R., and Boland, T. (2019). Thermal inkjet bioprinting triggers the activation of the VEGF pathway in human microvascular endothelial cells in vitro. Biofabrication. doi: 10.1088/1758-5090/ab25f9

PubMed Abstract | CrossRef Full Text | Google Scholar

Son, J., Lee, J. H., Kim, H. N., Ha, H., and Lee, Z. H. (2010). cAMP-response-element-binding protein positively regulates breast cancer metastasis and subsequent bone destruction. Biochem. Biophys. Res. Commun. 398, 309–314. doi: 10.1016/j.bbrc.2010.06.087

PubMed Abstract | CrossRef Full Text | Google Scholar

Spink David, C., Spink, B. C., Cao, J. Q., DePasquale, J. A., Pentecost, B. T., Fasco, M. J., et al. (1998). Differential expression of CYP1A1 and CYP1B1 in human breast epithelial cells and breast tumor cells. Carcinogenesis 19, 291–298. doi: 10.1093/carcin/19.2.291

CrossRef Full Text | Google Scholar

Stein, B., Brady, H., Yang, M. X., Young, D. B., and Barbosa, M. S. (1996). Cloning and characterization of MEK6, a novel member of the mitogen-activated protein kinase kinase cascade. J. Biol. Chem. 271, 11427–11433. doi: 10.1074/jbc.271.19.11427

PubMed Abstract | CrossRef Full Text | Google Scholar

Swinney, D. C., and Anthony, J. (2011). How were new medicines discovered? Nat. Rev. Drug Discov. 10:507. doi: 10.1038/nrd3480

PubMed Abstract | CrossRef Full Text | Google Scholar

Tabatabaian, F., Dougherty, K., Di Fulvio, M., and Gomez-Cambronero, J. (2010). Mammalian target of rapamycin (mTOR) and S6 kinase down-regulate phospholipase D2 basal expression and function. J. Biol. Chem. 285, 18991–19001. doi: 10.1074/jbc.M110.111542

PubMed Abstract | CrossRef Full Text | Google Scholar

Tamura, K., Sudo, T., Senftleben, U., Dadak, A. M., Johnson, R., and Karin, M. (2000). Requirement for p38α in erythropoietin expression: a role for stress kinases in erythropoiesis. Cell 102, 221–231. doi: 10.1016/S0092-8674(00)00027-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Thai, P., Statt, S., Chen, C. H., Liang, E., Campbell, C., and Wu, R. (2013). Characterization of a novel long noncoding RNA, SCAL1, induced by cigarette smoke and elevated in lung cancer cell lines. Am. J. Respir. Cell. Mol. Biol. 49, 204–211. doi: 10.1165/rcmb.2013-0159RC

PubMed Abstract | CrossRef Full Text | Google Scholar

Tonelli, C., Christine Chio, I. I., and Tuveson, D. A. (2018). Transcriptional regulation by Nrf2. Antioxid Redox. Signal. 29, 1727–1745. doi: 10.1089/ars.2017.7342

PubMed Abstract | CrossRef Full Text | Google Scholar

Trapnell, C., Roberts, A., Goff, L., Pertea, G., Kim, D., Kelley, D. R., et al. (2012). Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat. Protoc. 7, 562–578. doi: 10.1038/nprot.2012.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Ugolkov, A., Gaisina, I., Zhang, J. S., Billadeau, D. D., White, K., Kozikowski, A., et al. (2016). GSK-3 inhibition overcomes chemoresistance in human breast cancer. Cancer Lett. 380, 384–392. doi: 10.1016/j.canlet.2016.07.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Vaseva, A. V., Marchenko, N. D., Ji, K., Tsirka, S. E., Holzmann, S., and Moll, U. M. (2012). p53 opens the mitochondrial permeability transition pore to trigger necrosis. Cell 149, 1536–1548. doi: 10.1016/j.cell.2012.05.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Vos, M. J., Kanon, B., and Kampinga, H. H. (2009). HSPB7 is a SC35 speckle resident small heat shock protein. Biochim. Biophys. Acta 1793, 1343–1353. doi: 10.1016/j.bbamcr.2009.05.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Watson, G., Ronai, Z. A., and Lau, E. (2017). ATF2, a paradigm of the multifaceted regulation of transcription factors in biology and disease. Pharmacol. Res. 119, 347–357. doi: 10.1016/j.phrs.2017.02.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Weinberg, R. (2013). The Biology of Cancer. New York, NY: Garland Science.

Google Scholar

Wickman, G. R., Julian, L., Mardilovich, K., Schumacher, S., Munro, J., Rath, N., et al. (2013). Blebs produced by actin–myosin contraction during apoptosis release damage-associated molecular pattern proteins before secondary necrosis occurs. Cell Death Different. 20:1293. doi: 10.1038/cdd.2013.69

PubMed Abstract | CrossRef Full Text | Google Scholar

Wilson, W. C., and Boland, T. (2003). Cell and organ printing 1: protein and cell printers. Anat. Rec. A Discov. Mol. Cell. Evol. Biol. 272, 491–496. doi: 10.1002/ar.a.10057

PubMed Abstract | CrossRef Full Text | Google Scholar

Woodgett, J. R. (2005). Recent advances in the protein kinase B signaling pathway. Curr. Opin. Cell Biol. 17, 150–157. doi: 10.1016/j.ceb.2005.02.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, T., Zhao, W., Zhu, J. M., Albanna, M. Z., Yoo, J. J., and Atala, A. (2013). Complex heterogeneous tissue constructs containing multiple cell types prepared by inkjet printing technology. Biomaterials 34, 130–139. doi: 10.1016/j.biomaterials.2012.09.035

PubMed Abstract | CrossRef Full Text | Google Scholar

Yanez, M., Rincon, J., Dones, A., De Maria, C., Gonzales, R., and Boland, T. (2014). in vivo assessment of printed microvasculature in a bilayer skin graft to treat full-thickness wounds. Tissue Eng. Part A 21, 224–233. doi: 10.1089/ten.tea.2013.0561

PubMed Abstract | CrossRef Full Text | Google Scholar

Yerneni, S. S., Whiteside, T. L., Weiss, L. E., and Campbell, P. G. (2019). Bioprinting exosome-like extracellular vesicle microenvironments. Bioprinting 13:e00041. doi: 10.1016/j.bprint.2019.e00041

CrossRef Full Text | Google Scholar

Yoon, J. H., You, B. H., Park, C. H., Kim, Y. J., Nam, J. W., and Lee, S. K. (2018). The long noncoding RNA LUCAT1 promotes tumorigenesis by controlling ubiquitination and stability of DNA methyltransferase 1 in esophageal squamous cell carcinoma. Cancer Lett. 417, 47–57. doi: 10.1016/j.canlet.2017.12.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Zheng, A., Song, X., Zhang, L., Zhao, L., Mao, X., Wei, M., et al. (2019). Long non-coding RNA LUCAT1/miR-5582-3p/TCF7L2 axis regulates breast cancer stemness via Wnt/β-catenin pathway. J. Exp. Clin. Cancer Res. 38:305. doi: 10.1186/s13046-019-1315-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: bioprinting, tumor model, molecular properties, drug discovery, kinase phosphorylation

Citation: Campbell A, Mohl JE, Gutierrez DA, Varela-Ramirez A and Boland T (2020) Thermal Bioprinting Causes Ample Alterations of Expression of LUCAT1, IL6, CCL26, and NRN1L Genes and Massive Phosphorylation of Critical Oncogenic Drug Resistance Pathways in Breast Cancer Cells. Front. Bioeng. Biotechnol. 8:82. doi: 10.3389/fbioe.2020.00082

Received: 27 September 2019; Accepted: 29 January 2020;
Published: 21 February 2020.

Edited by:

Eric Farrell, Erasmus University Rotterdam, Netherlands

Reviewed by:

Daniela Duarte Campos, Stanford University, United States
Roberta Tasso, University of Genoa, Italy

Copyright © 2020 Campbell, Mohl, Gutierrez, Varela-Ramirez and Boland. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Thomas Boland, tboland@utep.edu

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.