Next Article in Journal
The Importance of Tumor–Host Interactions in Adult B-Cell Leukemias and Lymphomas
Next Article in Special Issue
The Role of PKM2 in Metabolic Reprogramming: Insights into the Regulatory Roles of Non-Coding RNAs
Previous Article in Journal
Classifications of Neurodegenerative Disorders Using a Multiplex Blood Biomarkers-Based Machine Learning Model
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 21
Issue 18
10.3390/ijms21186913
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

3D Bone Morphology Alters Gene Expression, Motility, and Drug Responses in Bone Metastatic Tumor Cells

by
Ushashi C. Dadwal
1,2,3,†,
Alyssa R. Merkel
1,3,4,†,
Jonathan M. Page
2,
Kristin A. Kwakwa
1,3,
Michael Kessler
5 and
Julie A. Rhoades
1,3,4,6,*
1
Department of Veterans Affairs, Tennessee Valley Healthcare System, Nashville, TN 37212, USA
2
Department of Chemical and Biomolecular Engineering, Vanderbilt University, Nashville, TN 37235, USA
3
Center for Bone Biology, Vanderbilt University Medical Center, Nashville, TN 37232, USA
4
Division of Clinical Pharmacology, Department of Medicine, Vanderbilt University Medical Center, Nashville, TN 37232, USA
5
School for Science and Math, Vanderbilt University, Nashville, TN 37235, USA
6
Department of Biomedical Engineering, Vanderbilt University, Nashville, TN 37235, USA
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Int. J. Mol. Sci. 2020, 21(18), 6913; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms21186913
Submission received: 27 August 2020 / Revised: 16 September 2020 / Accepted: 17 September 2020 / Published: 21 September 2020
(This article belongs to the Special Issue Signal Transduction Pathways Involved in Regulation of Normal and Tumour Cells)
Browse Figures
Review Reports Versions Notes

Abstract

:
Patients with advanced skeletal metastases arising from primary cancers including breast, lung, and prostate suffer from extreme pain, bone loss, and frequent fractures. While the importance of interactions between bone and tumors is well-established, our understanding of complex cell–cell and cell–microenvironment interactions remains limited in part due to a lack of appropriate 3D bone models. To improve our understanding of the influence of bone morphometric properties on the regulation of tumor-induced bone disease (TIBD), we utilized bone-like 3D scaffolds in vitro and in vivo. Scaffolds were seeded with tumor cells, and changes in cell motility, proliferation, and gene expression were measured. Genes associated with TIBD significantly increased with increasing scaffold rigidity. Drug response differed when tumors were cultured in 3D compared to 2D. Inhibitors for Integrin β3 and TGF-β Receptor II significantly reduced bone-metastatic gene expression in 2D but not 3D, while treatment with the Gli antagonist GANT58 significantly reduced gene expression in both 2D and 3D. When tumor-seeded 3D scaffolds were implanted into mice, infiltration of myeloid progenitors changed in response to pore size and rigidity. This study demonstrates a versatile 3D model of bone used to study the influence of mechanical and morphometric properties of bone on TIBD.

1. Introduction

Tumor-induced bone disease (TIBD) is a common occurrence among primary breast, lung, and prostate cancer patients, with up to 80% of patients who die from the disease displaying bone metastases upon post-mortem analysis [1]. Bone metastases lead to chronic bone pain, hypercalcemia, osteolytic lesions, pathological fractures, and bone loss, which significantly decrease patient quality of life [2]. Prior studies have demonstrated that in the primary tumor, mechanical forces exerted on tumor cells by the microenvironment promote integral changes in cell morphology, gene expression, cellular metabolism, protein synthesis, and response to cytotoxic therapeutics [3,4]. The unique architecture of trabecular bone contributes to mechanical forces on cells resulting from the mineralized extracellular matrix [5,6], curved trabecular surfaces [7,8,9,10,11], interstitial fluid flow [12,13,14], and dynamic compressive loading [15]. These bone-specific mechanical forces can affect the colonization and establishment of tumor cells in bone. While many research groups have studied the molecular events leading to tumor establishment and bone destruction, the specific cellular and mechanical cues that govern these responses have not been explored in the context of the 3D bone-tumor microenvironment [16,17].
Biomaterials mimicking the bone-tumor microenvironment have been highlighted as critical for understanding the dynamics of tumor progression, recapitulating in vivo conditions at distinct steps of metastasis, and determining how tumor cells integrate mechanical and chemical signals over multiple length scales [18]. Biomimetic 3D tissue-engineered scaffolds have been proposed as alternatives to 2D culture and mouse models for investigating molecular mechanisms of disease progression and for screening drugs [17,18,19,20]. Hydroxyapatite [21], collagen [22,23], poly(lactic-co-glycolic acid) (PLGA) [24], and poly(ε-caprolactone) (PCL) [25] scaffolds have been investigated as 3D models for assessing tumor cell gene expression and drug response in TIBD. However, the ability of electrospun polymeric scaffolds and hydrogels to recapitulate the mechanical and morphometric properties of the bone microenvironment may be limited by their low substrate modulus and nanoscale pore size. Consequently, the effects of mechanical and morphometric properties on bone-metastatic gene expression and drug response have not been systematically investigated for scaffolds with substrate moduli and pore sizes representative of trabecular bone.
In this study, we designed 3D-printed scaffolds to investigate the effects of mechanical forces in the bone microenvironment on progression of TIBD. We previously reported that the substrate modulus of 2D films regulates tumor cell expression of the transcription factor Gli2 and parathyroid hormone-related protein (PTHrP), both of which have been associated with bone destruction [5,6]. However, in the 3D bone microenvironment, cells are exposed to mechanical forces resulting from surface curvature [26] in addition to the substrate modulus, which stimulates mechanically transduced signaling in tumor [6] and bone [27,28] cells. To more accurately recapitulate the mechanical and morphometric properties of the bone microenvironment, we fabricated 3D scaffolds with tunable substrate moduli and pore size using a templated-Fused Deposition Modeling (3D t-FDM) process [28,29]. The bulk substrate modulus ranged from 5 to 266 MPa, which is representative of values reported for collagen fibers (32 MPa) and trabecular bone (93–366 MPa). Pore sizes in the physiologically relevant range that supports bone formation and vascularization (300–600 µm) were also evaluated [5,27,28,29,30]. Scaffolds were seeded with bone-metastatic tumor cells to investigate the effects of substrate modulus and pore size on bone-metastatic gene expression and drug response in vitro using multiple cancer cell lines, and in vivo using a breast cancer xenograft model.

2. Results

2.1. 3D Scaffolds Recapitulate the Mechanical and Morphometric Properties of Trabecular Bone

We fabricated 3D scaffolds designed to recapitulate both bone-specific mechanical forces arising from the rigid extracellular matrix and macropores between curved trabeculae. Two different substrate moduli were used: (i) compliant (C, 5 MPa) scaffolds that are less rigid than a collagen fiber (32 MPa), and (ii) rigid (R, 266 MPa) scaffolds that are representative of trabecular bone (93–365 MPa) [5,28,30]. Pore sizes were selected to be within the range that supports bone formation and vascularization (>300 μm). The smaller pore size (420 μm) is within the optimal range for osteogenesis [31,32,33,34], while the larger pore size (560 μm) is representative of tumor-bearing trabecular bone (Trabecular Separation, Tb.Sp.) ranging from 430 to 930 μm [28,30,32]. Four scaffolds were fabricated: (1) compliant scaffolds with 560 μm pores (560C), (2) compliant scaffolds with 420 μm pores (420C), (3) rigid scaffolds with 560 μm pores (560R), and (4) rigid scaffolds with 420 μm pores (420R). Representative images of the scaffolds are in Figure 1B and scaffold properties are listed in Table 1.

2.2. 3D Scaffolds Support Proliferation of Tumor Cells

Scaffolds were seeded with bone-tropic MDA-MB-231 tumor cells, and cells were evaluated for changes in gene expression, cell growth, and motility in response to varying pore sizes and rigidity. Viability of MDA-MB-231 and patient-derived chondrosarcoma cells assessed by Trypan blue staining demonstrated < 5% cell death after 48 h, indicating that the 3D scaffolds are not cytotoxic (Supplemental Figures S1 and S2). Scanning electron microscopy (SEM) images show that the tumor cells spread on and adhere to the surface of all four scaffold groups (Figure 1C). Previous studies have shown expression of Collagen-1 and Fibronection in MSCs cultured on 3D scaffolds [28]. Here, we demonstrate changes in expression of proteins associated with tumor attachment and adhesion to the environment including Focal Adhesion Kinase (FAK) and Integrin Beta 3 (Iβ3). FAK expression increases on MDA-MB-231 cells grown on rigid scaffolds compared to compliant scaffolds as shown by increased number of focal adhesions per cell (Supplemental Figure S3). We further investigated whether pore size and substrate modulus affected cell proliferation from day 1 to 3. Cell proliferation as measured by increases in cell number (GFP fluorescence expression) or metabolic activity (MTS assay) was not significantly affected by substrate modulus or pore size (Figure 1D,E).

2.3. Substrate Modulus and Pore Size Regulate the Motility of Tumor Cells

Since previous studies have established that the microenvironment can alter cell motility and invasive potential [16,17,35,36], we investigated the effects of substrate modulus and pore size on motility and migration of bone-metastatic MDA-MB-231-GFP tumor cells cultured on 2D and 3D substrates. 2D films with 5 MPa (2D Compliant) or 266 MPa (2D Rigid) bulk substrate moduli were synthesized from the same components as the 3D scaffolds to minimize the effects of chemical composition on cell fate (Figure 2A). Using live cell tracking, we found that tumor cells were significantly more motile on 3D scaffolds compared to 2D films (Supplemental Movies S1–S4). We observed 6- and 7-fold higher cell displacements for the 560R and 420R scaffolds compared to the 2D compliant films (Figure 2B). When tumor cell-seeded scaffolds were placed in serum-free media for a period of 6 h and exposed to 10% FBS media (chemoattractant, CM), cell speed increased significantly (5-fold 560R and 4-fold 420R) with increasing substrate rigidity but was independent of pore size (Figure 2C).
To determine if this increased motility can result in an increase in tumor migration, we placed tumor seeded scaffolds on top of a transwell and measured the number of cells that migrated through. Transwell migration assays demonstrated significantly higher migration potential of cells on 420R (≈3-fold) and 560R (≈2.5-fold) scaffolds compared to compliant scaffolds in complete media (CM), while there were no significant changes in migration potential without a chemoattractant gradient (SFM) (Figure 2D). Taken together, these observations suggest that the rigidity of the 3D microenvironment, but not pore size, can increase cell motility.

2.4. 3D Scaffolds Influence the Expression of Bone-Metastatic Genes In Vitro

To test the effects of substrate modulus and pore size, both parameters that influence mechanical signaling, on gene expression in bone-metastatic tumor cells, MDA-MB-231, RWGT2, and PC3 cells were seeded onto 2D films or 3D scaffolds and cultured for 48 h. Bone-metastatic gene expression was analyzed by qRT-PCR. ITGB3, Gli2, and PTHrP expression were significantly affected by both increasing rigidity and alterations in pore size. ITGB3 expression was 10-fold higher in MDA-MB-231 cells, 5–7-fold higher in RWGT2 cells, and 5–10-fold higher in PC3 cells on rigid compared to compliant scaffolds (Figure 3A). Furthermore, there was a 2-fold significant increase in ITGB3 gene expression in PC3 cells grown in a 560 μM scaffold compared to 460 μM scaffolds or 2D films. Gli2 expression significantly increased almost 2-fold with increasing rigidity and decreasing pore size in MDA-MB-231 cells, while Gli2 expression was highest (2-fold) in 560R scaffolds for RWGT2 and PC3 cells (Figure 3B). PTHrP expression was 10-fold higher in MDA-MB-231 and RWGT2 cells, and 3-fold higher in PC3 cells on rigid compared to compliant scaffolds (Figure 3C). PTHrP increased with decreasing pore size in MDA-MB-231 cells but pore size differences were not observed for RWGT2 and PC3 cells. These data suggest that substrate modulus and pore size regulate expression of genes associated with bone metastasis in breast cancer (MDA-MB-231), lung cancer (RWGT2), and prostate cancer (PC3).

2.5. 3D Scaffolds Influence the Response of Tumor Cells to Therapeutics

To further explore the effect of the 3D bone microenvironment on tumor cell behavior, we tested the drug response of the tumor cell lines to three inhibitors in short-term mono-culture on 2D vs. 3D rigid scaffolds. MDA-MB-231, RWGT2, and PC3 cells were cultured on rigid 2D films or 420R 3D scaffolds and expression of ITGB3, GLI2, and PTHrP was measured by qRT-PCR after 48 h of drug treatment. The TGF-β Receptor I kinase inhibitor (SD-208) and the Integrin inhibitor (Cilengitide, Cil) significantly reduced expression of genes associated with TIBD in 2D films by 2–3-fold; however, these drugs were less or not effective in 3D scaffolds (Figure 4A,B). In contrast, treatment with the Gli antagonist GANT58 both consistently and significantly reduced bone-metastatic gene expression >3-fold in all three cell lines and in both culture conditions. Similar experiments were performed on compliant films and scaffolds; however, low expression caused a decrease in yield. In support, molecular inhibition of these pathways in MDA-MB-231 cells using a mutant TGF-β receptor type II construct (rII Δcyst) [37], shRNA for integrin β3 (shβ3), or a Gli2 repressor construct (Gli2 EnR) reduced gene expression in both 2D and 3D (Supplemental Figure S4). Taken together, these data reveal that the characteristics of the bone microenvironment like rigidity and 3-dimensional structure mediate the efficacy of therapeutics, highlighting the importance of 3D culture systems when studying bone-metastatic tumor cell response.

2.6. Substrate Modulus and Pore Size Mediate Expression of Bone-Metastatic Genes in a Breast Cancer Xenograft Model

In order to determine how the 3D microenvironment regulates tumor growth in vivo, we implanted MDA-MB-231-seeded 3D scaffolds in the abdomen of athymic nude mice. After 3.5 weeks, scaffolds were harvested for histomorphometry analysis or gene expression analysis (Figure 5A). Expression of genes associated with TIBD increased with greater substrate modulus and decreasing pore size (Figure 5B–D). ITGB3, GLI2, and PTHrP mRNA significantly increased 20-, 18-, and 35-fold on 420R compared to 560C scaffolds. Immunohistochemistry for the proliferation marker Ki67 and the apoptosis marker Caspase 3 (Casp3) show increased cell proliferation on 420R scaffolds and little apoptosis activity on all scaffolds (Supplemental Figure S5). We also determined via a Tumor Metastatic panel several genes associated with bone remodeling were highly expressed in the 3D rigid scaffolds when compared to the complaint scaffolds (Supplemental Table S1). These data suggest that these 3D scaffolds could be a powerful tool for studying tumor growth and response to therapeutics in vivo.

2.7. Substrate Modulus Influences Immune Cell Infiltration

Previous studies have emphasized the contributions of the host microenvironment and immune response to tumor behavior during tumor establishment, proliferation, and growth [38]. We investigated whether the host immune cell populations that infiltrated the scaffolds created a pro-tumorigenic niche by examining the myeloid (CD11b+), macrophages (F4/80+), and granulocytic (Ly6G+) immune cell populations that are known to support TIBD [39,40,41]. We confirmed a significant 6-fold increase in the population of CD11b+ cells and a 10-fold increase in F4/80+ cells in 420R compared to 560C 3D scaffolds (Figure 6A–D). Similar to the bone-metastatic gene expression data, the population of CD11b+ and F4/80+ cells increased with increasing substrate modulus and decreasing pore size. We also detected the presence of Ly-6G+ positive cells in all scaffolds (Supplemental Figure S6). This data shows that the 3D scaffolds recapitulate the tumor–immune cell interactions in the bone microenvironment during TIBD, demonstrating that the 3D scaffolds are a robust and reliable model to study components of bone metastatic disease.

3. Discussion

In the present study, we investigated the effects of substrate modulus and pore size (surface curvature) on tumor cell gene expression and drug response using 3D scaffolds with mechanical and morphometric properties comparable to trabecular bone (Figure 7). Expression of genes associated with TIBD increased with increasing substrate modulus and decreasing pore size (increasing surface curvature). Furthermore, while the Gli antagonist GANT58 was effective in both 2D and 3D culture, Iβ3 and TGF-β inhibitors were not as effective on 3D bone-like scaffolds. These findings underscore the importance of screening tumor-targeting therapies in a microenvironment that is representative of host tissue.
Tissue-engineered 3D scaffolds are a promising technology for independently controlling the mechanical properties of the tissue microenvironment [18]. Previous studies using 3D scaffolds to model bone-metastatic breast cancer [5,42,43], Ewing’s sarcoma [25], and bone-metastatic prostate cancer [22,44] have reported significant differences between 2D and 3D culture. Ewing’s sarcoma cells cultured on 3D electrospun PCL scaffolds showed higher expression and activation of insulin-like growth factor-1 receptor (IGF-1R), as well as higher expression of proteins associated with resistance to IGF-1R therapy, compared to 2D mono-culture [25]. In another study, prostate cancer cells cultured in 3D hydroxyapatite (HA)-collagen scaffolds exhibited reduced expression of MMP1 and MMP9 compared to 2D culture [22]. Similarly, we have previously reported that expression of MMP9 by MDA-MB-231 cells cultured on 2D films decreased with substrate rigidity [6]. Coculture of prostate cancer cells with human osteoblasts on poly(ε-caprolactone) (PCL)-tricalcium phosphate (TCP) scaffolds exhibited upregulation of matrix metalloproteinases (MMPs) and prostate specific antigen (PSA) compared to prostate cancer cells alone [44]. While these studies highlight the substantial contribution of the 3D microenvironment to tumor cell fate, the relative contributions of mechanical forces resulting from substrate rigidity and surface curvature on tumor cell gene expression in bone have not been systematically investigated.
The complexity and spatial heterogeneity of the bone microenvironment [28,35,45,46] has made it difficult to predict which tumors will colonize the bone and induce bone destruction. Recently, we developed fabricated tissue-engineered bone constructs (TEBCs) from human trabecular bone templates that are capable of mimicking human bone morphometric properties—trabecular interconnectivity, surface roughness, and mechanical properties to study human mesenchymal stem cell (hMSCs) regulation of proliferation, differentiation, and mineralization [11]. Our previous study demonstrated hMSCs cultured on TEBCs exhibited significantly different metabolic activities, osteogenic differentiation, and mineralization depending on the anatomic site. Although many approaches to 3D models of bone and tumor exist, by using high throughput and tunable 3D scaffolds presented in this study we demonstrate a faster, more readily available, and more affordable 3D scaffold which is user-friendly to molecular biology laboratories.
Our studies demonstrated that these 3D scaffolds in vitro showed an important distinction between 2D and 3D culture. This difference was further exemplified while testing inhibitors that target components of TIBD. We also demonstrate in vivo that the 3D scaffolds reflected critical components of the tumor-bone microenvironment including immune cell infiltration, tumor proliferation, and pro-TIBD gene expression. Although these studies were performed with cells that primarily form osteolytic lesions, the 3D scaffolds described herein also have potential to be used to study tumors such as prostate cancer that form osteoblastic lesions. Based on these studies, while this model does not fully recapitulate bone, it improves upon tissue culture plastic or simple Matrigel studies to be a more bone mimicking in vitro microenvironment. For many in vitro molecular biology and in vivo cellular studies, this simpler approach is sufficient. However, we recognize that other approaches may still require more bone-like scaffolds like the TEBCs.
Trabecular bone comprises a complex mixture of rod and plate-like trabeculae with an average spacing of 600–800 µm apart and substrate modulus 93–365 MPa, which is several orders of magnitude higher than that of soft tissue [5,6,39,47]. Tumor cells in the bone microenvironment are subject to increased mechanical forces compared to soft tissue due to the mineralized extracellular matrix, curved trabeculae, and fluid shear stress. Consistent with our previous study using 2D films [5], expression of ITGB3, GLI2, and PTHRP by tumor cells cultured on 3D scaffolds increased with increasing substrate modulus due to interactions between integrins and growth factors. Similarly, we have reported that osteoblast differentiation of mesenchymal stem cells increased with increasing substrate modulus and decreasing pore size due to physical association of Iβ1 and BMP Receptor I (BMPRI) [30]. These findings are consistent with previous studies reporting that the magnitude of cellular contractile forces increases with increasing surface curvature [48,49]. Since trabecular separation decreases with proximity to the cortex, our findings suggest that tumor cells may become more invasive as they migrate toward the cortical wall.
Previous studies have reported that antitumor drugs are not as effective in the 3D bone microenvironment compared to 2D culture [17,22,25,50]. Therefore, we tested several inhibitors of TIBD in 2D and 3D culture to identify differences in drug response: Cilengitide (Iβ3 inhibitor), SD-208 (TGF-β Receptor I kinase inhibitor), and GANT58 (Gli antagonist). Cilengitide and SD-208 significantly reduced expression of ITGB3, GLI2, and PTHRP in 2D culture, but were less effective at reducing expression of these genes in 3D culture. While TGF-β inhibitors block metastasis to bone and subsequent bone destruction [51,52], they are less effective when given after tumors have established in bone [38]. However, it is well-established that TGF-β inhibitors increase bone volume independent of tumor [53,54] and thus would be beneficial even when given after tumors have established in bone. Iβ3 inhibitors have also been reported to block tumor metastasis to bone [52,55], as well as directly target osteoclasts and prevent bone resorption [56]. Our findings suggest that once tumors are established in bone, they are subjected to mechanical forces significantly greater than in the primary site and this increased signaling activation may reduce their ability to respond to Iβ3 and TGF-β inhibitors. In contrast, GANT58 reduced expression of ITGB3, GLI2, and PTHRP in tumor cells cultured on both 2D films and bone-like 3D scaffolds, which highlights the potential of Gli2 inhibitors for blocking TIBD. It was recently shown that targeting Gli2 inhibition to the tumor-bone microenvironment effectively reduced TIBD in a model of breast cancer metastasis to bone [57,58]. Thus, this study emphasizes the importance of studying tumor cell signaling and response to therapeutics in 3D cultures that accurately capture the physical properties of the bone microenvironment.

4. Materials and Methods

4.1. Cell Lines

A bone metastatic clone of the human breast cancer cell line MDA-MB-231 was derived by our lab as previously described [59]. The human squamous lung cancer cell line, RWGT2, was derived as previously described [60]. The human prostate cancer cell line, PC3, was purchased from ATCC (Manassas, VA, USA). MDA-MB-231-GFP cells were generated by transfecting a bone-metastatic variant of MDA-MB-231 cells with a Green Fluorescence Protein (GFP) overexpression plasmid and passaged through mice to further select for a bone metastatic clone [5,61,62]. MDA-shβ3, MDA-Gli2 EnR, and MDA-rII ΔCyst cells were generated for molecular inhibition experiments. Detailed methods can be found in supplementary methods section. Cells were tested every 6 months using the MycoAlert Mycoplasma Detection Kit (Lonza, Alpharetta, GA, USA) and every 2 years by the Vanderbilt Translational Pathology Shared Resource for multiple infectious agents.

4.2. Synthesis of 2D Poly(Ester Urethane) (PEUR) Films

2D poly(ester urethane) (PEUR) films were synthesized by reacting a mixture of poly(ε-caprolactone-co-glycolide) triol (Mn = 300 or 3000 g mol−1), hexamethylene diisocyanate trimer (HDIt), and stannous octoate catalyst as described previously [5]. To facilitate cell adhesion, fibronectin (Fn) was adsorbed to the surface of the substrates by incubating them in a 4 µg mL−1 solution of Fn in PBS at 4 °C overnight.

4.3. Fabrication and Characterization of 3D Scaffolds by templated-Fused Deposition Modeling (t-FDM)

3D scaffolds were fabricated as described previously (Figure 1A) [28,30]. Briefly, a poly(lactic acid) (PLA) template was printed using a MakerBot Replicator® 2 Fused Deposition Modeling (FDM) (Brooklyn, NY, USA) printer and filled with the PEUR mixture listed above. The PLA template was leached with dichloromethane (DCM) to yield 3D scaffolds with interconnected pores having a diameter of 423 ± 34 or 557 ± 44 μm for the nominal 420 or 560 μm templates, respectively. Pore size and morphology were characterized by scanning electron microscopy (SEM) (Hitachi S-4200, Finchampstead, Royal Bekshire, UK). Scaffolds were sterilized under UV light for 15 min in 70% ethanol and incubated in a solution of 4 μg mL−1 fibronectin overnight at 4 °C. As we have reported previously [28], the bulk substrate modulus (Ks, measured by compression testing) of the 3D scaffolds was controlled by the molecular weight (Mw) of the polyester triol (3000 or 300 g mol−1) to attain values representative of collagen (5 MPa, Compliant (C)) or trabecular bone (266 MPa, Rigid (R)) (Table 1).

4.4. Cell Culture

MDA-MB-231 cells were cultured in DMEM; RWGT2 cells were cultured in α-MEM; PC3 cells were cultured in RPMI growth media. All media was supplemented with 10% FBS (Peak Serum, Wellington, CO, USA) and 1 % Penicillin streptomycin solution (Mediatech, Manassa, VA, USA). 2D films were seeded with 3.6 × 105 tumor cells and 3D scaffolds were seeded with 0.5 × 106 cells/scaffold.

4.5. Cell Viability, Proliferation, and Metabolic Activity

Trypan Blue (Sigma, St. Louis, MO, USA) was used to measure cell viability 48 h after seeding (detailed methods can be found in supplementary methods). Metabolism was measured by MTS assay (CellTiter 96® Aqueous Non-Radioactive Cell Proliferation Assay, Promega, Madison, WI, USA) per the manufacturer’s protocol after 1, 2, and 3 days in culture (Figure 1D). Cell proliferation was measured by fluorescence intensity measurement in MDA-MB-231-GFP cells after 1, 2, and 3 days in culture using a Synergy 2 plate reader (Biotek Instruments, Winooski, VT, USA) (Figure 1E).

4.6. Cell Migration Assays

MDA-MB-231-GFP cells were plated on fibronectin-coated 2D films or 3D scaffolds, placed in a live cell chamber (LiveCell™, Okolab, Ambridge, PA, USA) at 5% CO2 and 37 °C and monitored by light microscopy (Olympus CKX41, Center Valley, PA, USA). Images of the same field were taken every hour for 48 h, and ImageJ software (NIH, Bethesda, MD, USA) [63] was used to analyze the photo series to track single cell (n = 3) movement. Chemotactic migration was quantified using a 24-well plate Boyden chamber transwell assay (8 μm pore size; Corning Costar, Tewksbury, MA, USA). 3D scaffolds were seeded with 4.0 × 104 tumor cells, serum-starved for 4 h, and mounted into the upper chamber of the transwell. Medium in the lower chamber was supplemented with 10% (positive control) or 0% (negative control) FBS. After 30 h, cells present on the upper surface of the filter were removed using a sterile cotton swab, and cells that migrated through the chamber onto the lower surface were fixed and stained with Crystal Violet Dye (Sigma, St. Louis, MO, USA). The number of migrating cells was counted (n = 5 fields) and compared to the negative control. Experiments were performed in triplicate wells and repeated at least three times.

4.7. Scanning Electron Microscopy

3D scaffolds were seeded with tumor cells and fixed with 4% paraformaldehyde for 1 h. Subsequently, samples were further fixed in a solution of 5% glutaraldehyde, followed by 2% osmium tetraoxide. Samples were dehydrated using ethanol washes of increasing concentration and vacuum dried overnight. Once fixed and dried, 3D scaffolds were cut open to expose pores and mounted on a stub using carbon tape. Samples were sputter-coated with gold using a 108 Auto Sputter Coater (Ted Pella Inc., Redding, CA, USA) and imaged on a Zeiss Merlin scanning electron microscope (Carl Zeiss Inc., Thornwood, White Plains, NY, USA).

4.8. Drug Treatments

2D films or 3D scaffolds were seeded with tumor cells, cultured for 24 h, and then treated with 10 μM Cilengitide (Selleck Chemicals, Huston, TX, USA), 10 μM GANT 58 (Santa Cruz Biotechnology, Dallas, TX, USA), 100 nM SD208 (Sigma-Aldrich, St. Louis, MO, USA), or DMSO as a control in serum-free media. After 48 h of drug treatment, RNA was harvested for qRT-PCR analysis.

4.9. Quantitative RT-PCR

To measure changes in gene expression, quantitative mRNA reverse transcription polymerase chain reaction was carried out using a TaqMan® gene expression assay. Briefly, cells were harvested with TRIzol (Invitrogen, Waltham, MA, USA) after 48 h in culture and total RNA was extracted following manufacturer’s instructions. The qScript cDNA supermix (Quanta, Beverly, MA, USA) was used to synthesize cDNA using 1 μg total RNA. Using validated TaqMan® primers, the expression of PTHLH (Hs00174969_m1), Gli2 (Hs01119974_m1), and ITGB3 (Hs01001469_m1) was measured in triplicate by quantitative RT-PCR with the 7500 Real-Time PCR System (Applied Biosciences, Waltham, MA, USA) under the following cycling conditions: 95 °C for 10 min, followed by 40 cycles of 95 °C for 15 s and 60 °C for 1 min. mRNA concentration was calculated against a standard curve using the absolute quantification method. 18SrRNA was used as an internal control.

4.10. In Vivo Studies

All animal procedures were approved by the Vanderbilt University Institutional Animal Care and Use Committee (IACUC, protocol M/06/213, approval date: 22 August 2013) and were conducted according to NIH guidelines. For in vivo studies, 1 × 106 MDA-MB-231-GFP cells were seeded into 3D scaffolds. For subcutaneous implants, female 4-week-old athymic nude mice (Mus musculus Hsd:Athymic Nude-Foxn1nu, Envigo, Indianapolis, IN, USA) were anesthetized by continuous isoflurane and an incision made on the ventral lower abdomen. The left inguinal mammary gland between the fourth and fifth mammary fat pad was transplanted with cell-seeded 3D scaffolds. Mice were imaged weekly for fluorescence using the MAESTRO™ (CRi, Hopkins, MA, USA) imaging system. The collected images were spectrally unmixed to remove background fluorescence. Mice were sacrificed at 3.5 weeks post-tumor cell inoculation and the tumors excised and processed for gene expression analysis or histomorphometry.

4.11. Histomorphometry

Excised tumors were formalin fixed for 48 h and processed for paraffin embedding. Then, 4 µm sections were cut using a Leica RM2255 (Buffalo Grove, IL, USA) and deparaffinized in xylenes and stained with Hematoxylin and Eosin for histological analyses. Immunohistochemical staining was performed to analyze infiltration of immune cells (CD11b, F4/80, Ly6G) and cellular proliferation (Ki67/Casp3). Detailed methods can be found in the supplemental methods section. All quantitative analyses were performed using Metamorph Analysis (Meta Imaging) or ImageJ software (NIH, Bethesda, MD, USA).

4.12. Statistical Analysis

All studies were performed in triplicate with at least n = 3 test substrates per replicate to ensure batch-to-batch reproducibility. ANOVA and the Tukey’s multi comparison post hoc test were used for statistical analyses with significance set at p < 0.05 unless otherwise stated. Cell culture data was tested for significance across time points, across materials, and interactions between time and materials with n = 3 independent biological experiments.

5. Conclusions

This study showed how substrate modulus and pore size affects tumor cell behavior using a tunable 3D culture system. We found that tumor cells grown in a more rigid environment with a smaller pore size were more motile and had amplified expression of genes associated with bone metastasis, and that the drug response in a 3D bone-like environment differed remarkably from 2D monoculture. Furthermore, this 3D system can be used in an in vivo orthotopic model to recapitulate the tumor-bone microenvironment. This 3D bone-like model will potentially enable screening of new therapeutics both in vitro and in vivo prior to preclinical testing.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/1422-0067/21/18/6913/s1, Figure S1. Rigidity and pore size do not affect viability of MDA-MB-231 cells on 3D scaffolds, Figure S2. Rigidity and pore size do not affect viability of patient derived chondrosarcoma cells on 3D scaffolds, Figure S3. Effect of rigidity and pore size on cell adhesion, Figure S4. Effect of TGFβ, Integrin β3, and Gli2 Inhibition on bone-metastatic gene expression in 3D culture, Figure S5 Representative in vivo immunohistochemistry images of scaffolds and quantification of Ki67 and Caspase3 staining, Figure S6 Representative in vivo immunohistochemistry images of scaffolds stained for ly6G by immunohistochemistry, Table S1: Metastatic panel in 3D t-FDM scaffolds in vivo, Movie S1 Live cell imaging on 560C 3D scaffolds, Movie S2 Live cell imaging on 420C 3D scaffolds, Movie S3 Live cell imaging on 560R 3D scaffolds, Movie S4 Live cell imaging on 420R 3D scaffolds.

Author Contributions

U.C.D. and A.R.M. conceptualized the study design, performed experiments and analysis, and wrote the manuscript. J.M.P., K.A.K., and M.K. contributed to data curation and analysis. J.A.R. supervised completion of work, acquired funding for the project, and reviewed and edited the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This publication was made possible by grant number R01 CA163499 (JAS) from the National Cancer Institute (NCI) at the National Institutes of Health and by grant number 1I01BX001957 (JAS) from the Veterans Administration. Its contents are solely the responsibility of the authors and do not necessarily represent the official views of NCI or the VA. UCD also acknowledges financial support from the Schlumberger Faculty for the Future pre-doctoral fellowship.

Acknowledgments

We would like to thank the Vanderbilt Confocal Core. Funding sources for this publication are from the use of the VUMC Cell Imaging Shared Resource (supported by National Institutes of Health grants CA68485, DK20593, DK58404, DK59637 and EY08126), Schlumberger Faculty for the Future pre-doctoral program and National Cancer Institute (part of the National Institutes of Health, under Award Number CA163499). In addition, we thank Scott Guelcher and Sichang Lu (Vanderbilt University) for assistance with material synthesis and characterization.

Conflicts of Interest

The authors declare that they have no conflict of interest.

References

  1. Johnson, R.W.; Schipani, E.; Giaccia, A.J. HIF targets in bone remodeling and metastatic disease. Pharmacol. Ther. 2015, 150, 169–177. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Saad, F.; Lipton, A.; Cook, R.; Chen, Y.-M.; Smith, M.; Coleman, R. Pathologic fractures correlate with reduced survival in patients with malignant bone disease. Cancer 2007, 110, 1860–1867. [Google Scholar] [CrossRef] [PubMed]
  3. Schrader, J.; Gordon-Walker, T.T.; Aucott, R.L.; Van Deemter, M.; Quaas, A.; Walsh, S.; Benten, D.; Forbes, S.J.; Wells, R.G.; Iredale, J.P. Matrix stiffness modulates proliferation, chemotherapeutic response, and dormancy in hepatocellular carcinoma cells. Hepatology 2011, 53, 1192–1205. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Tilghman, R.W.; Blais, E.M.; Cowan, C.R.; Sherman, N.E.; Grigera, P.R.; Jeffery, E.D.; Fox, J.W.; Blackman, B.R.; Tschumperlin, D.J.; A Papin, J.; et al. Matrix Rigidity Regulates Cancer Cell Growth by Modulating Cellular Metabolism and Protein Synthesis. PLoS ONE 2012, 7, e37231. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Page, J.M.; Merkel, A.R.; Ruppender, N.S.; Guo, R.; Dadwal, U.C.; Cannonier, S.A.; Basu, S.; Guelcher, S.A.; Sterling, J.A. Matrix rigidity regulates the transition of tumor cells to a bone-destructive phenotype through integrin β3 and TGF-β receptor type II. Biomaterials 2015, 64, 33–44. [Google Scholar] [CrossRef] [Green Version]
  6. Ruppender, N.S.; Merkel, A.R.; Martin, T.J.; Mundy, G.R.; Sterling, J.A.; Guelcher, S.A. Matrix Rigidity Induces Osteolytic Gene Expression of Metastatic Breast Cancer Cells. PLoS ONE 2010, 5, e15451. [Google Scholar] [CrossRef] [Green Version]
  7. Bianchi, M.; Edreira, E.R.U.; Wolke, J.G.C.; Birgani, Z.T.; Habibović, P.; Jansen, J.A.; Tampieri, A.; Marcacci, M.; Leeuwenburgh, S.; Beucken, J.J.J.P.V.D. Substrate geometry directs the in vitro mineralization of calcium phosphate ceramics. Acta Biomater. 2014, 10, 661–669. [Google Scholar] [CrossRef]
  8. Bidan, C.M.; Kommareddy, K.P.; Rumpler, M.; Kollmannsberger, P.; Bréchet, Y.J.M.; Fratzl, P.; Dunlop, J.W.C. How Linear Tension Converts to Curvature: Geometric Control of Bone Tissue Growth. PLoS ONE 2012, 7, e36336. [Google Scholar] [CrossRef] [Green Version]
  9. Bidan, C.M.; Kommareddy, K.; Rumpler, M.; Kollmannsberger, P.; Fratzl, P.; Dunlop, J.W.C. Geometry as a Factor for Tissue Growth: Towards Shape Optimization of Tissue Engineering Scaffolds. Adv. Healthc. Mater. 2012, 2, 186–194. [Google Scholar] [CrossRef]
  10. Gamsjäger, E.; Bidan, C.; Fischer, F.; Fratzl, P.; Dunlop, J.W.C. Modelling the role of surface stress on the kinetics of tissue growth in confined geometries. Acta Biomater. 2013, 9, 5531–5543. [Google Scholar] [CrossRef]
  11. Vanderburgh, J.; Sterling, J.A.; Guelcher, S.A. 3D Printing of Tissue Engineered Constructs for In Vitro Modeling of Disease Progression and Drug Screening. Ann. Biomed. Eng. 2016, 45, 164–179. [Google Scholar] [CrossRef] [Green Version]
  12. Cowin, S.C.; Cardoso, L. Blood and interstitial flow in the hierarchical pore space architecture of bone tissue. J. Biomech. 2014, 48, 842–854. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Fritton, S.P.; Weinbaum, S. Fluid and Solute Transport in Bone: Flow-Induced Mechanotransduction. Annu. Rev. Fluid Mech. 2009, 41, 347–374. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Bixel, M.G.; Kusumbe, A.P.; Ramasamy, S.; Sivaraj, K.K.; Butz, S.; Vestweber, D.; Adams, R.H. Flow Dynamics and HSPC Homing in Bone Marrow Microvessels. Cell Rep. 2017, 18, 1804–1816. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Harmata, A.J.; Uppuganti, S.; Granke, M.; Guelcher, S.A.; Nyman, J.S. Compressive fatigue and fracture toughness behavior of injectable, settable bone cements. J. Mech. Behav. Biomed. Mater. 2015, 51, 345–355. [Google Scholar] [CrossRef] [Green Version]
  16. Lynch, M.E.; Fischbach, C. Biomechanical forces in the skeleton and their relevance to bone metastasis: Biology and engineering considerations. Adv. Drug Deliv. Rev. 2014, 119–134. [Google Scholar] [CrossRef] [Green Version]
  17. Fischbach, C.; Chen, R.; Matsumoto, T.; Schmelzle, T.; Brugge, J.S.; Polverini, P.J.; Mooney, D.J. Engineering tumors with 3D scaffolds. Nat. Methods 2007, 4, 855–860. [Google Scholar] [CrossRef]
  18. Schuessler, T.K.; Chan, X.Y.; Chen, H.J.; Ji, K.; Park, K.M.; Roshan-Ghias, A.; Sethi, P.; Thakur, A.; Tian, X.; Villasante, A.; et al. Biomimetic tissue-engineered systems for advancing cancer research: NCI Strategic Workshop report. Cancer Res. 2014, 74, 5359–5363. [Google Scholar] [CrossRef] [Green Version]
  19. Infanger, D.W.; Lynch, M.E.; Fischbach, C. Engineered Culture Models for Studies of Tumor-Microenvironment Interactions. Annu. Rev. Biomed. Eng. 2013, 15, 29–53. [Google Scholar] [CrossRef]
  20. Aveic, S.; Davtalab, R.; Vogt, M.; Weber, M.; Buttler, P.; Tonini, G.P.; Fischer, H. Calcium phosphate scaffolds with defined interconnecting channel structure provide a mimetic 3D niche for bone marrow metastasized tumor cell growth. Acta Biomater. 2019, 88, 527–539. [Google Scholar] [CrossRef]
  21. Pathi, S.P.; Kowalczewski, C.; Tadipatri, R.; Fischbach, C. A Novel 3-D Mineralized Tumor Model to Study Breast Cancer Bone Metastasis. PLoS ONE 2010, 5, e8849. [Google Scholar] [CrossRef] [PubMed]
  22. Fitzgerald, K.A.; Guo, J.; Tierney, E.G.; Curtin, C.M.; Malhotra, M.; Darcy, R.; O’Brien, F.J.; O’Driscoll, C.M. The use of collagen-based scaffolds to simulate prostate cancer bone metastases with potential for evaluating delivery of nanoparticulate gene therapeutics. Biomaterials 2015, 66, 53–66. [Google Scholar] [CrossRef] [PubMed]
  23. Graham, J.S.; Vomund, A.N.; Phillips, C.L.; Grandbois, M. Structural changes in human type I collagen fibrils investigated by force spectroscopy. Exp. Cell Res. 2004, 299, 335–342. [Google Scholar] [CrossRef] [PubMed]
  24. Yang, Z.; Zhao, X. A 3D model of ovarian cancer cell lines on peptide nanofiber scaffold to explore the cell–scaffold interaction and chemotherapeutic resistance of anticancer drugs. Int. J. Nanomed. 2011, 6, 303–310. [Google Scholar] [CrossRef] [Green Version]
  25. Fong, E.L.S.; Lamhamedi-Cherradi, S.-E.; Burdett, E.; Ramamoorthy, V.; Lazar, A.J.; Kasper, F.K.; Farach-Carson, M.C.; Vishwamitra, D.; Demicco, E.G.; Menegaz, B.A.; et al. Modeling Ewing sarcoma tumors in vitro with 3D scaffolds. Proc. Natl. Acad. Sci. USA 2013, 110, 6500–6505. [Google Scholar] [CrossRef] [Green Version]
  26. McCoy, R.J.; O’Brien, F.J. Influence of Shear Stress in Perfusion Bioreactor Cultures for the Development of Three-Dimensional Bone Tissue Constructs: A Review. Tissue Eng. Part B: Rev. 2010, 16, 587–601. [Google Scholar] [CrossRef]
  27. Guo, R.; Lu, S.; Merkel, A.R.; Sterling, J.A.; Guelcher, S.A. Substrate Modulus Regulates Osteogenic Differentiation of Rat Mesenchymal Stem Cells through Integrin β1 and BMP Receptor Type IA. J. Mater. Chem. B 2016, 4, 3584–3593. [Google Scholar] [CrossRef] [Green Version]
  28. Guo, R.; Lu, S.; Page, J.M.; Merkel, A.R.; Basu, S.; Sterling, J.A.; Guelcher, S.A. Fabrication of 3D Scaffolds with Precisely Controlled Substrate Modulus and Pore Size by Templated-Fused Deposition Modeling to Direct Osteogenic Differentiation. Adv. Healthc. Mater. 2015, 4, 1826–1832. [Google Scholar] [CrossRef]
  29. Gupta, M.K.; Meyer, T.A.; Nelson, C.E.; Duvall, C.L. Poly(PS-b-DMA) micelles for reactive oxygen species triggered drug release. J. Control. Release 2012, 162, 591–598. [Google Scholar] [CrossRef] [Green Version]
  30. Guo, R.; Merkel, A.; Sterling, J.; Davidson, J.; Guelcher, S.A. Substrate modulus of 3D-printed scaffolds regulates the regenerative response in subcutaneous implants through the macrophage phenotype and Wnt signaling. Biomaterials 2015, 73, 85–95. [Google Scholar] [CrossRef] [Green Version]
  31. Cyster, L.; Grant, D.M.; Howdle, S.; Rose, F.R.A.J.; Irvine, D.; Freeman, D.; Scotchford, C.A.; Shakesheff, K.M. The influence of dispersant concentration on the pore morphology of hydroxyapatite ceramics for bone tissue engineering. Biomaterials 2005, 26, 697–702. [Google Scholar] [CrossRef] [PubMed]
  32. Karageorgiou, V.; Kaplan, D. Porosity of 3D biomaterial scaffolds and osteogenesis. Biomaterials 2005, 26, 5474–5491. [Google Scholar] [CrossRef] [PubMed]
  33. Penk, A.; Forster, Y.; Scheidt, H.A.; Nimptsch, A.; Hacker, M.; Schulz, M.B.; Ahnert, P.; Schiller, J.; Rammelt, S.; Huster, D. The pore size of PLGA bone implants determines the de novo formation of bone tissue in tibial head defects in rats. Magn. Reson. Med. 2012, 70, 925–935. [Google Scholar] [CrossRef] [PubMed]
  34. Rossi, F.; Santoro, M.; Perale, G. Polymeric scaffolds as stem cell carriers in bone repair. J. Tissue Eng. Regen. Med. 2013, 9, 1093–1119. [Google Scholar] [CrossRef]
  35. Guise, T.A. Breast cancer bone metastases: It’s all about the neighborhood. Cell 2013, 154, 957–959. [Google Scholar] [CrossRef] [Green Version]
  36. Zhang, X.H.-F.; Jin, X.; Malladi, S.; Zou, Y.; Wen, Y.H.; Brogi, E.; Smid, M.; Foekens, J.A.; Massagué, J. Selection of bone metastasis seeds by mesenchymal signals in the primary tumor stroma. Cell 2013, 154, 1060–1073. [Google Scholar] [CrossRef] [Green Version]
  37. Yin, J.J.; Selander, K.; Chirgwin, J.M.; Dallas, M.; Grubbs, B.G.; Wieser, R.; Massagué, J.; Mundy, G.R.; A Guise, T. TGF-β signaling blockade inhibits PTHrP secretion by breast cancer cells and bone metastases development. J. Clin. Investig. 1999, 103, 197–206. [Google Scholar] [CrossRef] [Green Version]
  38. Buenrostro, D.; Kwakwa, K.A.; Putnam, N.E.; Merkel, A.R.; Johnson, J.R.; Cassat, J.; Sterling, J.A. Early TGF-β inhibition in mice reduces the incidence of breast cancer induced bone disease in a myeloid dependent manner. Bone 2018, 113, 77–88. [Google Scholar] [CrossRef]
  39. Azarin, S.M.; Yi, J.; Gower, R.M.; Aguado, B.A.; Sullivan, M.E.; Goodman, A.G.; Jiang, E.J.; Rao, S.S.; Ren, Y.; Tucker, S.L.; et al. In vivo capture and label-free detection of early metastatic cells. Nat. Commun. 2015, 6, 8094. [Google Scholar] [CrossRef]
  40. Roca, H.; McCauley, L.K. Inflammation and skeletal metastasis. BoneKEy Rep. 2015, 4, 706. [Google Scholar] [CrossRef] [Green Version]
  41. Yang, L.; Huang, J.; Ren, X.; Gorska, A.E.; Chytil, A.; Aakre, M.; Carbone, D.P.; Matrisian, L.M.; Richmond, A.; Lin, P.C.; et al. Abrogation of TGFβ Signaling in Mammary Carcinomas Recruits Gr-1+CD11b+ Myeloid Cells that Promote Metastasis. Cancer Cell 2008, 13, 23–35. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Mastro, A.M.; Vogler, E.A. A Three-Dimensional Osteogenic Tissue Model for the Study of Metastatic Tumor Cell Interactions with Bone. Cancer Res. 2009, 69, 4097–4100. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Moreau, J.E.; Anderson, K.; Mauney, J.R.; Nguyen, T.; Kaplan, D.L.; Rosenblatt, M. Tissue-Engineered Bone Serves as a Target for Metastasis of Human Breast Cancer in a Mouse Model. Cancer Res. 2007, 67, 10304–10308. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Sieh, S.; Taubenberger, A.; Lehman, M.L.; Clements, J.; Nelson, C.C.; Hutmacher, D.W. Paracrine interactions between LNCaP prostate cancer cells and bioengineered bone in 3D in vitro culture reflect molecular changes during bone metastasis. Bone 2014, 63, 121–131. [Google Scholar] [CrossRef] [PubMed]
  45. Sterling, J.A.; Edwards, J.R.; Martin, T.J.; Mundy, G.R. Advances in the biology of bone metastasis: How the skeleton affects tumor behavior. Bone 2011, 48, 6–15. [Google Scholar] [CrossRef]
  46. Weilbaecher, K.N.; Guise, T.A.; McCauley, L.K. Cancer to bone: A fatal attraction. Nat. Rev. Cancer 2011, 11, 411–425. [Google Scholar] [CrossRef]
  47. E Lynch, M.; Brooks, D.; Mohanan, S.; Lee, M.J.; Polamraju, P.; Dent, K.; Bonassar, L.J.; Van Der Meulen, M.C.H.; Fischbach, C. In vivo tibial compression decreases osteolysis and tumor formation in a human metastatic breast cancer model. J. Bone Miner. Res. 2013, 28, 2357–2367. [Google Scholar] [CrossRef]
  48. Rumpler, M.; Woesz, A.; Dunlop, J.W.C.; Van Dongen, J.T.; Fratzl, P. The effect of geometry on three-dimensional tissue growth. J. R. Soc. Interface 2008, 5, 1173–1180. [Google Scholar] [CrossRef] [Green Version]
  49. Sanz-Herrera, J.; Moreo, P.; García-Aznar, J.M.; Doblaré, M. On the effect of substrate curvature on cell mechanics. Biomaterials 2009, 30, 6674–6686. [Google Scholar] [CrossRef]
  50. Pickl, M.; Ries, C.H. Comparison of 3D and 2D tumor models reveals enhanced HER2 activation in 3D associated with an increased response to trastuzumab. Oncogene 2008, 28, 461–468. [Google Scholar] [CrossRef] [Green Version]
  51. Mohammad, K.S.; Javelaud, D.; Fournier, P.G.; Niewolna, M.; McKenna, C.R.; Peng, X.H.; Duong, V.; Dunn, L.K.; Mauviel, A.; Guise, T.A. TGF-beta-RI kinase inhibitor SD-208 reduces the development and progression of melanoma bone metastases. Cancer Res. 2010, 71, 175–184. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Nyman, J.S.; Merkel, A.R.; Uppuganti, S.; Nayak, B.; Rowland, B.; Makowski, A.J.; Oyajobi, B.O.; Sterling, J.A. Combined treatment with a transforming growth factor beta inhibitor (1D11) and bortezomib improves bone architecture in a mouse model of myeloma-induced bone disease. Bone 2016, 91, 81–91. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Biswas, S.; Nyman, J.S.; Alvarez, J.; Chakrabarti, A.; Ayres, A.; Sterling, J.; Edwards, J.; Rana, T.; Johnson, R.; Perrien, D.S.; et al. Anti-Transforming Growth Factor ß Antibody Treatment Rescues Bone Loss and Prevents Breast Cancer Metastasis to Bone. PLoS ONE 2011, 6, e27090. [Google Scholar] [CrossRef] [PubMed]
  54. Edwards, J.R.; Nyman, J.S.; Lwin, S.T.; Moore, M.M.; Esparza, J.; O’Quinn, E.C.; Hart, A.J.; Biswas, S.; A Patil, C.; Lonning, S.; et al. Inhibition of TGF-β signaling by 1D11 antibody treatment increases bone mass and quality in vivo. J. Bone Miner. Res. 2010, 25, 2419–2426. [Google Scholar] [CrossRef]
  55. Zhao, Y.; Bachelier, R.; Treilleux, I.; Pujuguet, P.; Peyruchaud, O.; Baron, R.; Clément-Lacroix, P.; Clezardin, P. Tumor v 3 Integrin Is a Therapeutic Target for Breast Cancer Bone Metastases. Cancer Res. 2007, 67, 5821–5830. [Google Scholar] [CrossRef] [Green Version]
  56. Nakamura, I.; Pilkington, M.F.; Lakkakorpi, P.T.; Lipfert, L.; Sims, S.M.; Dixon, S.J.; A Rodan, G.; Duong, L.T. Role of alpha(v)beta(3) integrin in osteoclast migration and formation of the sealing zone. J. Cell Sci. 1999, 112, 3985–3993. [Google Scholar]
  57. Vanderburgh, J.P.; Kwakwa, K.A.; Werfel, T.A.; Merkel, A.R.; Gupta, M.K.; Johnson, R.W.; Guelcher, S.A.; Duvall, C.L.; Rhoades, J.A. Systemic delivery of a Gli inhibitor via polymeric nanocarriers inhibits tumor-induced bone disease. J. Control. Release 2019, 311-312, 257–272. [Google Scholar] [CrossRef]
  58. Vanderburgh, J.; Hill, J.L.; Gupta, M.K.; Kwakwa, K.A.; Wang, S.K.; Moyer, K.; Bedingfield, S.K.; Merkel, A.R.; D’Arcy, R.; Guelcher, S.A.; et al. Tuning Ligand Density To Optimize Pharmacokinetics of Targeted Nanoparticles for Dual Protection against Tumor-Induced Bone Destruction. ACS Nano 2020, 14, 311–327. [Google Scholar] [CrossRef]
  59. Sterling, J.A.; Oyajobi, B.O.; Grubbs, B.; Padalecki, S.S.; Munoz, S.A.; Gupta, A.; Story, B.; Zhao, M.; Mundy, G.R. The Hedgehog Signaling Molecule Gli2 Induces Parathyroid Hormone-Related Peptide Expression and Osteolysis in Metastatic Human Breast Cancer Cells. Cancer Res. 2006, 66, 7548–7553. [Google Scholar] [CrossRef] [Green Version]
  60. Guise, T.A.; Yoneda, T.; Yates, A.J.; Mundy, G.R. The combined effect of tumor-produced parathyroid hormone-related protein and transforming growth factor-alpha enhance hypercalcemia in vivo and bone resorption in vitro. J. Clin. Endocrinol. Metab. 1993, 77, 40–45. [Google Scholar]
  61. Johnson, R.W.; Nguyen, M.P.; Padalecki, S.S.; Grubbs, B.G.; Merkel, A.R.; Oyajobi, B.O.; Matrisian, L.M.; Mundy, G.R.; Sterling, J.A. TGF-beta promotion of Gli2-induced expression of parathyroid hormone-related protein, an important osteolytic factor in bone metastasis, is independent of canonical Hedgehog signaling. Cancer Res. 2010, 71, 822–831. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Campbell, J.P.; Merkel, A.R.; Masood-Campbell, S.K.; Elefteriou, F.; Sterling, J.A. Models of Bone Metastasis. J. Vis. Exp. 2012, 10, e4260. [Google Scholar] [CrossRef] [PubMed]
  63. Schneider, C.A.; Rasband, W.S.; Eliceiri, K.W. NIH Image to ImageJ: 25 years of image analysis. Nat. Methods 2012, 9, 671–675. [Google Scholar] [CrossRef] [PubMed]
Ijms 21 06913 g001 550
Figure 1. Fabrication schematic for 3D t-FDM scaffolds and cellular activity on the scaffolds. (A) Hexamethylene diisocyanate trimer (HDIt) was mixed with a poly(ε-caprolactone-co-glycolide) triol and stannous octoate catalyst and poured into a poly(lactic acid) (PLA) mold fabricated by a MakerBot Replicator 2 fused-deposition modeling (FDM) 3D printer. After curing overnight, the PLA mold was extracted via solvent leaching in dichloromethane for 24 h. The resulting scaffold was washed with DI water and prepared for cell culture. (B,C) Representative scanning electron microscopy (SEM) image of the 3D t-FDM scaffolds (B) and of cell attachment on 3D scaffolds (C)—560 µm compliant (560C), 420 µm compliant (420C), 560 µm rigid (560R), and 420 µm rigid (420R). (D) Metabolic activity of MDA-MB-231 tumor cells cultured on 3D scaffolds measured by the MTS assay corroborates no variation of metabolic expression over a time period (days- D1, D2 and D3). (E) Proliferation of MDA-MB-231 tumor cells on 3D scaffolds measured as fluorescence expression of GFP, demonstrates no statistical differences between 3D scaffold groups over a time period (days- D1, D2 and D3).
Figure 1. Fabrication schematic for 3D t-FDM scaffolds and cellular activity on the scaffolds. (A) Hexamethylene diisocyanate trimer (HDIt) was mixed with a poly(ε-caprolactone-co-glycolide) triol and stannous octoate catalyst and poured into a poly(lactic acid) (PLA) mold fabricated by a MakerBot Replicator 2 fused-deposition modeling (FDM) 3D printer. After curing overnight, the PLA mold was extracted via solvent leaching in dichloromethane for 24 h. The resulting scaffold was washed with DI water and prepared for cell culture. (B,C) Representative scanning electron microscopy (SEM) image of the 3D t-FDM scaffolds (B) and of cell attachment on 3D scaffolds (C)—560 µm compliant (560C), 420 µm compliant (420C), 560 µm rigid (560R), and 420 µm rigid (420R). (D) Metabolic activity of MDA-MB-231 tumor cells cultured on 3D scaffolds measured by the MTS assay corroborates no variation of metabolic expression over a time period (days- D1, D2 and D3). (E) Proliferation of MDA-MB-231 tumor cells on 3D scaffolds measured as fluorescence expression of GFP, demonstrates no statistical differences between 3D scaffold groups over a time period (days- D1, D2 and D3).
Ijms 21 06913 g001
Ijms 21 06913 g002 550
Figure 2. Rigidity and pore size affect migration of tumor cells on 3D scaffolds. (A) Schematic of experimental procedure (B) Live cell imaging using bright field microscopy demonstrates significantly higher cell motility on 560R (>6-fold) and 420R (>7-fold) scaffolds as compared to compliant films (2D Compliant, white column). (C) Motility assay using confocal imaging shows that tumor cells have significantly increased speed on the 560R (>8-fold) and 420R (>7-fold) scaffolds compared to compliant scaffolds. (D) 3D migration transwell assay with either 10% FBS media (CM) or serum-free medium (SFM) determined tumor cells have higher migratory potential on 560R (3-fold) and 420R (>2.5-fold) 3D scaffolds as compared to 560C and 420C scaffolds. 2way ANOVA. Compliant vs. Rigid, **** p < 0.0001. 2D vs. 3D, #### p < 0.0001. 560 vs. 420, + p < 0.05.
Figure 2. Rigidity and pore size affect migration of tumor cells on 3D scaffolds. (A) Schematic of experimental procedure (B) Live cell imaging using bright field microscopy demonstrates significantly higher cell motility on 560R (>6-fold) and 420R (>7-fold) scaffolds as compared to compliant films (2D Compliant, white column). (C) Motility assay using confocal imaging shows that tumor cells have significantly increased speed on the 560R (>8-fold) and 420R (>7-fold) scaffolds compared to compliant scaffolds. (D) 3D migration transwell assay with either 10% FBS media (CM) or serum-free medium (SFM) determined tumor cells have higher migratory potential on 560R (3-fold) and 420R (>2.5-fold) 3D scaffolds as compared to 560C and 420C scaffolds. 2way ANOVA. Compliant vs. Rigid, **** p < 0.0001. 2D vs. 3D, #### p < 0.0001. 560 vs. 420, + p < 0.05.
Ijms 21 06913 g002
Ijms 21 06913 g003 550
Figure 3. Effects of substrate modulus and pore size on gene expression of bone-metastatic tumor cells. The breast cancer cell line, MDA-MB-231 (black), the lung cancer cell line, RWGT2 (red), and the prostate cancer cell line, PC3 (blue), were seeded on 2D films or 3D scaffolds, cultured for 48 h and analyzed for changes in gene expression. Expression of (A) ITGB3, (B) Gli2, and (C) PTHrP were significantly increased for all cell types analyzed with respect to changes in both pore size and rigidity. Data presented as fold change over 2D compliant. Two-way ANOVA. Compliant vs. rigid, * p < 0.05, ** p < 0.01, *** p < 0.001, **** p < 0.0001. 560 vs. 420, + p < 0.05, ++++ p < 0.0001. 2D vs. 3D, # p < 0.05, ## p < 0.01, #### p < 0.0001.
Figure 3. Effects of substrate modulus and pore size on gene expression of bone-metastatic tumor cells. The breast cancer cell line, MDA-MB-231 (black), the lung cancer cell line, RWGT2 (red), and the prostate cancer cell line, PC3 (blue), were seeded on 2D films or 3D scaffolds, cultured for 48 h and analyzed for changes in gene expression. Expression of (A) ITGB3, (B) Gli2, and (C) PTHrP were significantly increased for all cell types analyzed with respect to changes in both pore size and rigidity. Data presented as fold change over 2D compliant. Two-way ANOVA. Compliant vs. rigid, * p < 0.05, ** p < 0.01, *** p < 0.001, **** p < 0.0001. 560 vs. 420, + p < 0.05, ++++ p < 0.0001. 2D vs. 3D, # p < 0.05, ## p < 0.01, #### p < 0.0001.
Ijms 21 06913 g003
Ijms 21 06913 g004 550
Figure 4. Effects of 3D culture on tumor cell drug response. Effects of 3D culture on inhibition of MDA-MB-231 (white—2D scaffolds, black—3D scaffolds), RWGT2 (pink—2D scaffolds, red—3D scaffolds), and PC3 (light blue—2D scaffolds, dark blue—3D scaffolds) tumor cells was tested for (A) Integrin β3 inhibition using the inhibitor Cilengitide (Cil), (B) TGF-β inhibition using an RI kinase inhibitor (SD208), and (C) Gli2 inhibition using the Gli antagonist (Gant58). Cells were grown on 2D rigid films (pale bars) or 420R 3D scaffolds (dark bars) and mRNA expression was analyzed by qRT-PCR after 48 h in culture. ITGB3, Gli2, and PTHrP mRNA expression decreased in response to SD208 and Cilengitide treatment in 2D films, but to a lesser extent in 3D scaffolds. Gli inhibition was effective in both 2D and 3D and all three tumor cell lines. Data presented as fold change over untreated control. Two-way ANOVA. 3D scaffolds/2D scaffolds treatment vs. control, * p < 0.05, ** p < 0.01, *** p < 0.001, **** p < 0.0001.
Figure 4. Effects of 3D culture on tumor cell drug response. Effects of 3D culture on inhibition of MDA-MB-231 (white—2D scaffolds, black—3D scaffolds), RWGT2 (pink—2D scaffolds, red—3D scaffolds), and PC3 (light blue—2D scaffolds, dark blue—3D scaffolds) tumor cells was tested for (A) Integrin β3 inhibition using the inhibitor Cilengitide (Cil), (B) TGF-β inhibition using an RI kinase inhibitor (SD208), and (C) Gli2 inhibition using the Gli antagonist (Gant58). Cells were grown on 2D rigid films (pale bars) or 420R 3D scaffolds (dark bars) and mRNA expression was analyzed by qRT-PCR after 48 h in culture. ITGB3, Gli2, and PTHrP mRNA expression decreased in response to SD208 and Cilengitide treatment in 2D films, but to a lesser extent in 3D scaffolds. Gli inhibition was effective in both 2D and 3D and all three tumor cell lines. Data presented as fold change over untreated control. Two-way ANOVA. 3D scaffolds/2D scaffolds treatment vs. control, * p < 0.05, ** p < 0.01, *** p < 0.001, **** p < 0.0001.
Ijms 21 06913 g004
Ijms 21 06913 g005 550
Figure 5. Effects of substrate modulus and pore size on bone metastatic gene expression in a xenograft model. (A) Schematic of animal model. Scaffolds were seeded with 1 × 106 MDA-MB-231-GFP cells and implanted subcutaneously in the abdomen. After 3.5 weeks, scaffolds were harvested for qRT-PCR analysis to assess (B) ITGB3, (C) Gli2, and (D) PTHrP mRNA expression. mRNA expression was significantly upregulated on 420R scaffolds as compared to 560C 3D scaffolds. Data presented as fold change over 560C. Two-way ANOVA. Compliant vs. rigid, ** p < 0.01, *** p < 0.001, **** p < 0.0001. 560 vs. 420, ++++ p < 0.0001.
Figure 5. Effects of substrate modulus and pore size on bone metastatic gene expression in a xenograft model. (A) Schematic of animal model. Scaffolds were seeded with 1 × 106 MDA-MB-231-GFP cells and implanted subcutaneously in the abdomen. After 3.5 weeks, scaffolds were harvested for qRT-PCR analysis to assess (B) ITGB3, (C) Gli2, and (D) PTHrP mRNA expression. mRNA expression was significantly upregulated on 420R scaffolds as compared to 560C 3D scaffolds. Data presented as fold change over 560C. Two-way ANOVA. Compliant vs. rigid, ** p < 0.01, *** p < 0.001, **** p < 0.0001. 560 vs. 420, ++++ p < 0.0001.
Ijms 21 06913 g005
Ijms 21 06913 g006 550
Figure 6. Effects of substrate modulus and pore size on immune infiltration in a xenograft model. Representative immunohistochemistry images show scaffolds with pores outlined by dashed lines. Tissue infiltrated inside pores is shown by Hematoxylin staining, and positive immunohistochemical staining is represented by dark brown staining. (A) Representative 20× images and (B) quantification of immunohistochemistry staining for CD11b in 3D scaffolds show increased CD11b+ cell populations in 560R (5-fold) and 420R (6-fold) scaffolds compared to compliant scaffolds. (C) Representative 20× images and (D) quantification of immunohistochemistry staining for F4/80 in 3D scaffolds reveal increased populations of F4/80+ cells in 560R (7-fold) and 420R (10-fold) as well as a noticeable increase in staining with decreasing pore size. All data was normalized to 560C. Compliant vs. rigid * p < 0.05, ** p < 0.01. 560 vs. 420, # p < 0.05.
Figure 6. Effects of substrate modulus and pore size on immune infiltration in a xenograft model. Representative immunohistochemistry images show scaffolds with pores outlined by dashed lines. Tissue infiltrated inside pores is shown by Hematoxylin staining, and positive immunohistochemical staining is represented by dark brown staining. (A) Representative 20× images and (B) quantification of immunohistochemistry staining for CD11b in 3D scaffolds show increased CD11b+ cell populations in 560R (5-fold) and 420R (6-fold) scaffolds compared to compliant scaffolds. (C) Representative 20× images and (D) quantification of immunohistochemistry staining for F4/80 in 3D scaffolds reveal increased populations of F4/80+ cells in 560R (7-fold) and 420R (10-fold) as well as a noticeable increase in staining with decreasing pore size. All data was normalized to 560C. Compliant vs. rigid * p < 0.05, ** p < 0.01. 560 vs. 420, # p < 0.05.
Ijms 21 06913 g006
Ijms 21 06913 g007 550
Figure 7. Schematic summary of the 3D bone-like culture system. Bone-like 3D culture of metastatic tumor cells increases signaling activation of genes associated with TIBD, increases cell motility, and influences immune cell infiltration in an in vivo orthotopic model.
Figure 7. Schematic summary of the 3D bone-like culture system. Bone-like 3D culture of metastatic tumor cells increases signaling activation of genes associated with TIBD, increases cell motility, and influences immune cell infiltration in an in vivo orthotopic model.
Ijms 21 06913 g007
Table
Table 1. 3D scaffold characterization and in vitro study design. Material characterization and nomenclature of 2D films and 3D scaffolds.
Table 1. 3D scaffold characterization and in vitro study design. Material characterization and nomenclature of 2D films and 3D scaffolds.
Treatment GroupLabelPore Size (μm)Bulk Substrate Modulus
Ks (MPa)
2D compliant film2DCN/A5 ± 0.4
2D rigid film2DRN/A266 ± 27
3D compliant scaffolds (small pores)420C423 ± 345 ± 0.4
3D compliant scaffolds (large pores)560C557 ± 445 ± 0.4
3D rigid scaffolds (small pores)420R423 ± 34266 ± 27
3D rigid scaffolds (large pores)560R557 ± 44266 ± 27

Share and Cite

MDPI and ACS Style

Dadwal, U.C.; Merkel, A.R.; Page, J.M.; Kwakwa, K.A.; Kessler, M.; Rhoades, J.A. 3D Bone Morphology Alters Gene Expression, Motility, and Drug Responses in Bone Metastatic Tumor Cells. Int. J. Mol. Sci. 2020, 21, 6913. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms21186913

AMA Style

Dadwal UC, Merkel AR, Page JM, Kwakwa KA, Kessler M, Rhoades JA. 3D Bone Morphology Alters Gene Expression, Motility, and Drug Responses in Bone Metastatic Tumor Cells. International Journal of Molecular Sciences. 2020; 21(18):6913. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms21186913

Chicago/Turabian Style

Dadwal, Ushashi C., Alyssa R. Merkel, Jonathan M. Page, Kristin A. Kwakwa, Michael Kessler, and Julie A. Rhoades. 2020. "3D Bone Morphology Alters Gene Expression, Motility, and Drug Responses in Bone Metastatic Tumor Cells" International Journal of Molecular Sciences 21, no. 18: 6913. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms21186913

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop