Measuring the mechanical stress induced by an expanding multicellular tumor system: a case study

Glioblastoma mechanobiology at multiple length scales

Glioblastoma multiforme (GBM), a primary brain cancer, is one of the most aggressive forms of human cancer, with a very low patient survival rate. A characteristic feature of GBM is the diffuse infiltration of tumor cells into the surrounding brain extracellular matrix (ECM) that provide biophysical, topographical, and biochemical cues. In particular, ECM stiffness and composition is known to play a key role in controlling various GBM cell behaviors including proliferation, migration, invasion, as well as the stem-like state and response to chemotherapies. In this review, we discuss the mechanical characteristics of the GBM microenvironment at multiple length scales, and how biomaterial scaffolds such as polymeric hydrogels, and fibers, as well as microfluidic chip-based platforms have been employed as tissue mimetic models to study GBM mechanobiology. We also highlight how such tissue mimetic models can impact the field of GBM mechanobiology.

Elastocapillary effects determine early matrix deformation by glioblastoma cell spheroids

During cancer pathogenesis, cell-generated mechanical stresses lead to dramatic alterations in the mechanical and organizational properties of the extracellular matrix (ECM). To date, contraction of the ECM is largely attributed to local mechanical stresses generated during cell invasion, but the impact of "elastocapillary" effects from surface tension on the tumor periphery has not been examined. Here, we embed glioblastoma cell spheroids within collagen gels, as a model of tumors within the ECM. We then modulate the surface tension of the spheroids, such that the spheroid contracts or expands. Surprisingly, in both cases, at the far-field, the ECM is contracted toward the spheroids prior to cellular migration from the spheroid into the ECM. Through computational simulation, we demonstrate that contraction of the ECM arises from a balance of spheroid surface tension, cell-ECM interactions, and time-dependent, poroelastic effects of the gel. This leads to the accumulation of ECM near the periphery of the spheroid and the contraction of the ECM without regard to the expansion or contraction of the spheroid. These results highlight the role of tissue-level surface stresses and fluid flow within the ECM in the regulation of cell-ECM interactions.

The role of cell-matrix adhesion and cell migration in breast tumor growth and progression

During breast cancer progression, there is typically increased collagen deposition resulting in elevated extracellular matrix rigidity. This results in changes to cell-matrix adhesion and cell migration, impacting processes such as the epithelial-mesenchymal transition (EMT) and metastasis. We aim to investigate the roles of cell-matrix adhesion and cell migration on breast tumor growth and progression by studying the impacts of different types of extracellular matrices and their rigidities. We embedded MCF7 spheroids within three-dimensional (3D) collagen matrices and agarose matrices. MCF7 cells adhere to collagen but not agarose. Contrasting the results between these two matrices allows us to infer the role of cell-matrix adhesion. We found that MCF7 spheroids exhibited the fastest growth rate when embedded in a collagen matrix with a rigidity of 5.1 kPa (0.5 mg/mL collagen), whereas, for the agarose matrix, the rigidity for the fastest growth rate is 15 kPa (1.0% agarose) instead. This discrepancy is attributable to the presence of cell adhesion molecules in the collagen matrix, which initiates collagen matrix remodeling and facilitates cell migration from the tumor through the EMT. As breast tumors do not adhere to agarose matrices, it is suitable to simulate the cell-cell interactions during the early stage of breast tumor growth. We conducted further analysis to characterize the stresses exerted by the expanding spheroid on the agarose matrix. We identified two distinct MCF7 cell populations, namely, those that are non-dividing and those that are dividing, which exerted low and high expansion stresses on the agarose matrix, respectively. We confirmed this using Western blot which showed the upregulation of proliferating cell nuclear antigen, a proliferation marker, in spheroids grown in the 1.0% agarose (≈13 kPa). By treating the embedded MCF7 spheroids with an inhibitor or activator of myosin contractility, we showed that the optimum spheroids' growth can be increased or decreased, respectively. This finding suggests that tumor growth in the early stage, where cell-cell interaction is more prominent, is determined by actomyosin tension, which alters cell rounding pressure during cell division. However, when breast tumors begin generating collagen into the surrounding matrix, collagen remodeling triggers EMT to promote cell migration and invasion, ultimately leading to metastasis.

Toward Measuring the Mechanical Stresses Exerted by Branching Embryonic Airway Epithelial Explants in 3D Matrices of Matrigel

Numerous organs in the bodies of animals, including the lung, kidney, and mammary gland, contain ramified networks of epithelial tubes. These structures arise during development via a process known as branching morphogenesis. Previous studies have shown that mechanical forces directly impact this process, but the patterns of mechanical stress exerted by branching embryonic epithelia are not well understood. This is, in part, owing to a lack of experimental tools. Traditional traction force microscopy assays rely on the use of compliant hydrogels with well-defined mechanical properties. Isolated embryonic epithelial explants, however, have only been shown to branch in three-dimensional matrices of reconstituted basement membrane protein, or Matrigel, a biomaterial with poorly characterized mechanical behavior, especially in the regime of large deformations. Here, to compute the traction stresses generated by branching epithelial explants, we quantified the finite-deformation constitutive behavior of gels of reconstituted basement membrane protein subjected to multi-axial mechanical loads. We then modified the mesenchyme-free assay for the ex vivo culture of isolated embryonic airway epithelial explants by suspending fluorescent microspheres within the surrounding gel and tracking their motion during culture. Surprisingly, the tracked bead motion was non-zero in regions of the gel far away from the explants, suggestive of passive swelling deformations within the matrix. To compute accurate traction stresses, these swelling deformations must be decomposed from those generated by the branching explants. We thus tracked the motion of beads suspended within cell-free matrices and quantified spatiotemporal patterns of gel swelling. Taken together, these passive swelling data can be combined with the measured mechanical properties of the gel to compute the traction forces exerted by intact embryonic epithelial explants.

Targeting the tumor biophysical microenvironment to reduce resistance to immunotherapy

Immunotherapy based on immune checkpoint inhibitors has evolved into a new pillar of cancer treatment in clinics, but dealing with treatment resistance (either primary or acquired) is a major challenge. The tumor microenvironment (TME) has a substantial impact on the pathological behaviors and treatment response of many cancers. The biophysical clues in TME have recently been considered as important characteristics of cancer. Furthermore, there is mounting evidence that biophysical cues in TME play important roles in each step of the cascade of cancer immunotherapy that synergistically contribute to immunotherapy resistance. In this review, we summarize five main biophysical cues in TME that affect resistance to immunotherapy: extracellular matrix (ECM) structure, ECM stiffness, tumor interstitial fluid pressure (IFP), solid stress, and vascular shear stress. First, the biophysical factors involved in anti-tumor immunity and therapeutic antibody delivery processes are reviewed. Then, the causes of these five biophysical cues and how they contribute to immunotherapy resistance are discussed. Finally, the latest treatment strategies that aim to improve immunotherapy efficacy by targeting these biophysical cues are shared. This review highlights the biophysical cues that lead to immunotherapy resistance, also supplements their importance in related technologies for studying TME biophysical cues in vitro and therapeutic strategies targeting biophysical cues to improve the effects of immunotherapy.

Prostate Surface Distension and Tumor Texture Descriptors From Pre-Treatment MRI Are Associated With Biochemical Recurrence Following Radical Prostatectomy: Preliminary Findings

Objective

To derive and evaluate the association of prostate shape distension descriptors from T2-weighted MRI (T2WI) with prostate cancer (PCa) biochemical recurrence (BCR) post-radical prostatectomy (RP) independently and in conjunction with texture radiomics of PCa.

Methods

This retrospective study comprised 133 PCa patients from two institutions who underwent 3T-MRI prior to RP and were followed up with PSA measurements for ≥3 years. A 3D shape atlas-based approach was adopted to derive prostate shape distension descriptors from T2WI, and these descriptors were used to train a random forest classifier (C_S) to predict BCR. Texture radiomics was derived within PCa regions of interest from T2WI and ADC maps, and another machine learning classifier (C_R) was trained for BCR. An integrated classifier C_S+R was then trained using predictions from C_S and C_R. These models were trained on D₁ (N = 71, 27 BCR+) and evaluated on independent hold-out set D₂ (N = 62, 12 BCR+). C_S+R was compared against pre-RP, post-RP clinical variables, and extant nomograms for BCR-free survival (bFS) at 3 years.

Results

C_S+R resulted in a higher AUC (0.75) compared to C_R (0.70, p = 0.04) and C_S (0.69, p = 0.01) on D₂ in predicting BCR. On univariable analysis, C_S+R achieved a higher hazard ratio (2.89, 95% CI 0.35–12.81, p < 0.01) compared to other pre-RP clinical variables for bFS. C_S+R, pathologic Gleason grade, extraprostatic extension, and positive surgical margins were associated with bFS (p < 0.05). C_S+R resulted in a higher C-index (0.76 ± 0.06) compared to CAPRA (0.69 ± 0.09, p < 0.01) and Decipher risk (0.59 ± 0.06, p < 0.01); however, it was comparable to post-RP CAPRA-S (0.75 ± 0.02, p = 0.07).

Conclusions

Radiomic shape descriptors quantifying prostate surface distension complement texture radiomics of prostate cancer on MRI and result in an improved association with biochemical recurrence post-radical prostatectomy.

Physical Forces in Glioblastoma Migration: A Systematic Review

The invasive capabilities of glioblastoma (GBM) define the cancer’s aggressiveness, treatment resistance, and overall mortality. The tumor microenvironment influences the molecular behavior of cells, both epigenetically and genetically. Current forces being studied include properties of the extracellular matrix (ECM), such as stiffness and “sensing” capabilities. There is currently limited data on the physical forces in GBM—both relating to how they influence their environment and how their environment influences them. This review outlines the advances that have been made in the field. It is our hope that further investigation of the physical forces involved in GBM will highlight new therapeutic options and increase patient survival. A search of the PubMed database was conducted through to 23 March 2022 with the following search terms: (glioblastoma) AND (physical forces OR pressure OR shear forces OR compression OR tension OR torsion) AND (migration OR invasion). Our review yielded 11 external/applied/mechanical forces and 2 tumor microenvironment (TME) forces that affect the ability of GBM to locally migrate and invade. Both external forces and forces within the tumor microenvironment have been implicated in GBM migration, invasion, and treatment resistance. We endorse further research in this area to target the physical forces affecting the migration and invasion of GBM.

Microsphere sensors for characterizing stress fields within three-dimensional extracellular matrix

Stress in the three-dimensional extracellular matrix is one of the key cues in regulating multiscale biological processes. Thus far, noticeable progress in methods and techniques (e.g., micropipette aspiration, AFM, and molecule probes) has been made to quantify stress in cell microenvironment at different length scales. Among them, the microsphere sensor-based method (MSS-based method) has emerged as an advantageous approach over conventional techniques in quantifying stress in situ and in vivo at cellular and supra-cellular scales. This method is implemented by seven sequential steps, including fabrication, modification, characterization, cell adhesion, imaging, displacement field extraction and stress calculation. Precise control of each step and inter-tunning between steps can provide quantitative characterization of stress field. However, detailed procedural information associated with each step and process has been scattered. This review aims to provide a comprehensive overview of MSS-based method, systematically summarizing the principles and research progresses. Firstly, the basic principles are introduced, and the specific experiment and calculation processes of MSS-based method are presented in detail. Then, recent advances and applications of this method are summarized. Finally, perspectives of the limitations and development trends of MSS-based method are discussed. This specific and comprehensive review would provide a guideline for the widespread application of MSS-based method as an advantageous method for in situ and in vivo stress characterization at cellular and supra-cellular scale within three-dimensional extracellular matrix. STATEMENT OF SIGNIFICANCE: In this review, a method based on a microsphere sensor (MSS-based method) as an advantageous approach over conventional techniques in quantifying stress in situ and in vivo at cellular and supra-cellular scales is introduced and discussed. This technique is implemented by seven sequential steps, including fabrication, modification, characterization, cell junction, imaging, displacement field extraction, and stress calculation. Precise control of each step and inter-tunning between steps can provide quantitative stress field. However, detailed procedural information associated with each step has been scattered. Thus, a comprehensive review collating recent advances and perspective discussions is a necessity to introduce a better option for quantifying the stress field in biological processes at the cellular and supra-cellular scales.

Cancer cell migration on straight, wavy, loop and grid microfibre patterns

Cell migration plays an important role in physiological and pathological processes where the fibrillar morphology of extracellular matrices (ECM) could regulate the migration dynamics. To mimic the morphological characteristics of fibrillar matrix structures, low-voltage continuous electrospinning was adapted to construct straight, wavy, looped and gridded fibre patterns made of polystyrene (of fibre diameter ca. 3 μm). Cells were free to explore their different shapes in response to the directly-adhered fibre, as well as to the neighbouring patterns. For all the patterns studied, analysing cellular migration dynamics of MDA-MB-231 (a highly migratory breast cancer cell line) demonstrated two interesting findings: first, although cells dynamically adjust their shapes and migration trajectories in response to different fibrillar environments, their average step speed is minimally affected by the fibre global pattern; secondly, a switch in behaviour was observed when the pattern features approach the upper limit of the cell body's minor axis, reflecting that cells' ability to divert from an existing fibre track is limited by the size along the cell body's minor axis. It is therefore concluded that the upper limit of cell body's minor axis might act as a guide for the design of microfibre patterns for different purposes of cell migration.

Cell – extracellular matrix interaction in glioma growth. In silico model

The study aims to investigate the role of viscoelastic interactions between cells and extracellular matrix (ECM) in avascular tumor growth. Computer simulations of glioma multicellular tumor spheroid (MTS) growth are being carried out for various conditions. The calculations are based on a continuous model, which simulates oxygen transport into MTS; transitions between three cell phenotypes, cell transport, conditioned by hydrostatic forces in cell–ECM composite system, cell motility and cell adhesion. Visco-elastic cell aggregation and elastic ECM scaffold represent two compressible constituents of the composite. Cell–ECM interactions form a Transition Layer on the spheroid surface, where mechanical characteristics of tumor undergo rapid transition. This layer facilitates tumor progression to a great extent. The study demonstrates strong effects of ECM stiffness, mechanical deformations of the matrix and cell–cell adhesion on tumor progression. The simulations show in particular that at certain, rather high degrees of matrix stiffness a formation of distant multicellular clusters takes place, while at further increase of ECM stiffness subtumors do not form. The model also illustrates to what extent mere mechanical properties of cell–ECM system may contribute into variations of glioma invasion scenarios.

Practical Review on Preclinical Human 3D Glioblastoma Models: Advances and Challenges for Clinical Translation

Fifteen years after the establishment of the Stupp protocol as the standard of care to treat glioblastomas, no major clinical advances have been achieved and increasing patient’s overall survival remains a challenge. Nevertheless, crucial molecular and cellular findings revealed the intra-tumoral and inter-tumoral complexities of these incurable brain tumors, and the essential role played by cells of the microenvironment in the lack of treatment efficacy. Taking this knowledge into account, fulfilling gaps between preclinical models and clinical samples is necessary to improve the successful rate of clinical trials. Since the beginning of the characterization of brain tumors initiated by Bailey and Cushing in the 1920s, several glioblastoma models have been developed and improved. In this review, we focused on the most widely used 3D human glioblastoma models, including spheroids, tumorospheres, organotypic slices, explants, tumoroids and glioblastoma-derived from cerebral organoids. We discuss their history, development and especially their usefulness.

The Extracellular Matrix and Biocompatible Materials in Glioblastoma Treatment

During cancer genesis, the extracellular matrix (ECM) in the human brain undergoes important transformations, starting to resemble embryonic brain cell milieu with a much denser structure. However, the stiffness of the tumor ECM does not preclude cancer cells from migration. The importance of the ECM role in normal brain tissue as well as in tumor homeostasis has engaged much effort in trials to implement ECM as a target and an instrument in the treatment of brain cancers. This review provides a detailed analysis of both experimental and applied approaches in combined therapy for gliomas in adults. In general, matrix materials for glioma treatment should have properties facilitating the simplest delivery into the body. Hence, to deliver an artificial implant directly into the operation cavity it should be packed into a gel form, while for bloodstream injections matrix needs to be in the form of polymer micelles, nanoparticles, etc. Furthermore, the delivered material should mimic biomechanical properties of the native tissue, support vital functions, and slow down or stop the proliferation of surrounding cells for a prolonged period. The authors propose a two-step approach aimed, on the one hand, at elimination of remaining cancer cells and on the other hand, at restoring normal brain tissue. Thereby, the first bioartificial matrix to be applied should have relatively low elastic modulus should be loaded with anticancer drugs, while the second material with a higher elastic modulus for neurite outgrowth support should contain specific factors stimulating neuroregeneration.

Computer simulations suggest that prostate enlargement due to benign prostatic hyperplasia mechanically impedes prostate cancer growth

Prostate cancer and benign prostatic hyperplasia are common genitourinary diseases in aging men. Both pathologies may coexist and share numerous similarities, which have suggested several connections or some interplay between them. However, solid evidence confirming their existence is lacking. Recent studies on extensive series of prostatectomy specimens have shown that tumors originating in larger prostates present favorable pathological features. Hence, large prostates may exert a protective effect against prostate cancer. In this work, we propose a mechanical explanation for this phenomenon. The mechanical stress fields that originate as tumors enlarge have been shown to slow down their dynamics. Benign prostatic hyperplasia contributes to these mechanical stress fields, hence further restraining prostate cancer growth. We derived a tissue-scale, patient-specific mechanically coupled mathematical model to qualitatively investigate the mechanical interaction of prostate cancer and benign prostatic hyperplasia. This model was calibrated by studying the deformation caused by each disease independently. Our simulations show that a history of benign prostatic hyperplasia creates mechanical stress fields in the prostate that impede prostatic tumor growth and limit its invasiveness. The technology presented herein may assist physicians in the clinical management of benign prostate hyperplasia and prostate cancer by predicting pathological outcomes on a tissue-scale, patient-specific basis.

The emergence of solid stress as a potent biomechanical marker of tumour progression

Cancer is a disease of dysregulated mechanics which alters cell behaviour, compromises tissue structure, and promotes tumour growth and metastasis. In the context of tumour progression, the most widely studied of biomechanical markers is matrix stiffness as tumour tissue is typically stiffer than healthy tissue. However, solid stress has recently been identified as another marker of tumour growth, with findings strongly suggesting that its role in cancer is distinct from that of stiffness. Owing to the relative infancy of the field which draws from diverse disciplines, a comprehensive knowledge of the relationships between solid stress, tumorigenesis, and metastasis is likely to provide new and valuable insights. In this review, we discuss the micro- and macro-scale biomechanical interactions that give rise to solid stresses, and also examine the techniques developed to quantify solid stress within the tumour environment. Moreover, by reviewing the effects of solid stress on tissues, cancer and stromal cells, and signalling pathways, we also detail its mode of action at each level of the cancer cascade.

Tropomyosin Tpm 2.1 loss induces glioblastoma spreading in soft brain-like environments

INTRODUCTION

The brain is a very soft tissue. Glioblastoma (GBM) brain tumours are highly infiltrative into the surrounding healthy brain tissue and invasion mechanisms that have been defined using rigid substrates therefore may not apply to GBM dissemination. GBMs characteristically lose expression of the high molecular weight tropomyosins, a class of actin-associating proteins and essential regulators of the actin stress fibres and focal adhesions that underpin cell migration on rigid substrates.

METHODS

Here, we investigated how loss of the high molecular weight tropomyosins affects GBM on soft matrices that recapitulate the biomechanical architecture of the brain.

RESULTS

We find that Tpm 2.1 is down-regulated in GBM grown on soft substrates. We demonstrate that Tpm 2.1 depletion by siRNA induces cell spreading and elongation in soft 3D hydrogels, irrespective of matrix composition. Tpm 1.7, a second high molecular weight tropomyosin is also down-regulated when cells are cultured on soft brain-like surfaces and we show that effects of this isoform are matrix dependent, with Tpm 1.7 inducing cell rounding in 3D collagen gels. Finally, we show that the absence of Tpm 2.1 from primary patient-derived GBMs correlates with elongated, mesenchymal invasion.

CONCLUSIONS

We propose that Tpm 2.1 down-regulation facilitates GBM colonisation of the soft brain environment. This specialisation of the GBM actin cytoskeleton organisation that is highly suited to the soft brain-like environment may provide novel therapeutic targets for arresting GBM invasion.

Fingering instabilities in tissue invasion: an active fluid model

Metastatic tumours often invade healthy neighbouring tissues by forming multicellular finger-like protrusions emerging from the cancer mass. To understand the mechanical context behind this phenomenon, we here develop a minimalist fluid model of a self-propelled, growing biological tissue. The theory involves only four mechanical parameters and remains analytically trackable in various settings. As an application of the model, we study the evolution of a two-dimensional circular droplet made of our active and expanding fluid, and embedded in a passive non-growing tissue. This system could be used to model the evolution of a carcinoma in an epithelial layer. We find that our description can explain the propensity of tumour tissues to fingering instabilities, as conditioned by the magnitude of active traction and the growth kinetics. We are also able to derive predictions for the tumour size at the onset of metastasis, and for the number of subsequent invasive fingers. Our active fluid model may help describe a wider range of biological processes, including wound healing and developmental patterning.

Adaptive adhesion systems mediate glioma cell invasion in complex environments

Diffuse brain invasion by glioma cells prevents effective surgical or molecular-targeted therapy and underlies a detrimental outcome. Migrating glioma cells are guided by complex anatomical brain structures but the exact mechanisms remain poorly defined. To identify adhesion receptor systems and matrix structures supporting glioma cell invasion into brain-like environments we used 2D and 3D organotypic invasion assays in combination with antibody-, peptide- and RNA-based interference. Combined interference with β1 and αV integrins abolished the migration of U-251 and E-98 glioma cells on reconstituted basement membrane; however, invasion into primary brain slices or 3D astrocyte-based scaffolds and migration on astrocyte-deposited matrix was only partly inhibited. Any residual invasion was supported by vascular structures, as well as laminin 511, a central constituent of basement membrane of brain blood vessels. Multi-targeted interference against β1, αV and α6 integrins expressed by U-251 and E-98 cells proved insufficient to achieve complete migration arrest. These data suggest that mechanocoupling by integrins is relatively resistant to antibody- or peptide-based targeting, and cooperates with additional, as yet unidentified adhesion systems in mediating glioma cell invasion in complex brain stroma.

Fabrication and Characterization of Chitosan-Hyaluronic Acid Scaffolds with Varying Stiffness for Glioblastoma Cell Culture

The invasive and recurrent nature of glioblastoma multiforme (GBM) is linked to a small subpopulation of cancer cells, which are self-renewing, resistant to standard treatment regimens, and induce formation of new tumors. Matrix stiffness is implicated in the regulation of cell proliferation, drug resistance, and reversion to a more invasive phenotype. Therefore, understanding the relationship between matrix stiffness and tumor cell behavior is vital to develop appropriate in vitro tumor models. Here, chitosan-hyaluronic acid (CHA) polyelectrolyte complex scaffolds are fabricated with statistically significant stiffness variances to characterize the effect of scaffold stiffness on morphology, proliferation, drug resistance, and gene expression in human glioblastoma cells (U-87 MG). All scaffolds support GBM proliferation over a 12-day culture period, yet larger spheroids are observed in scaffolds with higher stiffness. Additionally, GBM cells cultured in stiffer CHA scaffolds prove significantly more resistant to the common chemotherapeutic temozolomide. Moreover, the stiffer 8% CHA scaffolds exhibit an increase in expression of drug resistance and invasion related genes compared to 2D culture. CHA scaffolds present a tunable microenvironment for enhanced tumor cell malignancy and may provide a valuable in vitro microenvironment for studying tumor progression and screening anticancer therapies.

Brain-Tumor-Regenerating 3D Scaffold-Based Primary Xenograft Models for Glioma Stem Cell Targeted Drug Screening

Glioma stem cells (GSC) present a critical therapeutic challenge for glioblastoma multiforme (GBM). Drug screening against GSC demands development of novel in vitro and in vivo platforms that can mimic brain microenvironment and support GSC maintenance and tumorigenesis. Here, we report, a 3-dimensionel (3D) biomimetic macro-porous scaffold developed by incorporating hyaluronic acid, porcine brain extra cellular matrix (ECM) and growth factors that facilitates regeneration of GBM from primary GSCs, ex vivo and in vivo. After characterizing with human and rat GBM cell lines and neurospheres, human GSCs expressing Notch1, Sox-2, Nestin, and CD133 biomarkers were isolated from GBM patients, cultured in the 3D scaffold, and implanted subcutaneously in nude mice to develop patient derived xenograft (PDX) models. Aggressive growth pattern of PDX with formation of intratumoral vascularization was monitored by magnetic resonance imaging (MRI). Histopathological and phenotypial features of the original tumors were retained in the PDX models. We used this regenerated GBM platform to screen novel siRNA nanotherapeutics targeting Notch, Sox-2, FAK signaling for its ability to inhibit the tumorigenic potential of GSCs. Current clinical drug, Temozolomide and an anticancer phytochemical, nanocurcumin, were used as controls. The siRNA nanoparticles showed excellent efficacy in inhibiting tumorigenesis by GSCs in vivo. Our study suggests that the brain-ECM mimicking scaffold can regenerate primary gliomas from GSCs in vitro and in vivo, and the same can be used as an effective platform for screening drugs against glioma stem cells.

The mode and dynamics of glioblastoma cell invasion into a decellularized tissue-derived extracellular matrix-based three-dimensional tumor model

Glioblastoma multiforme (GBM) is the most common brain tumor with very aggressive and infiltrative. Extracellular matrix (ECM) plays pivotal roles in the infiltrative characteristics of GBM. To understand the invasive characteristic of GBM, it is necessary to study cell-ECM interaction in the physiologically relevant biomimetic model that recapitulates the GBM-specific ECM microenvironment. Here, we propose biomimetic GBM-specific ECM microenvironment for studying mode and dynamics of glioblastoma cell invasion. Using tissue decellularization process, we constructed a patient tissue-derived ECM (pdECM)-based three-dimensional in vitro model. In our model, GBM cells exhibited heterogeneous morphology and altered the invasion routes in a microenvironment-adaptive manner. We further elucidate the effects of inhibition of ECM remodeling-related enzymatic activity (Matrix metalloproteinase (MMP) 2/9, hyaluronan synthase (HAS)) on GBM cell invasion. Interestingly, after blocking both enzyme activity, GBM cells underwent morphological transition and switch the invasion mode. Such adaptability could render cell invasion resistant to anti-cancer target therapy. There results provide insight of how organ-specific matrix differentially regulates cancer cell phenotype, and have significant implications for the design of matrix with appropriate physiologically relevant properties for in vitro tumor model.

A mechanically coupled reaction-diffusion model that incorporates intra-tumoural heterogeneity to predict in vivo glioma growth

While gliomas have been extensively modelled with a reaction-diffusion (RD) type equation it is most likely an oversimplification. In this study, three mathematical models of glioma growth are developed and systematically investigated to establish a framework for accurate prediction of changes in tumour volume as well as intra-tumoural heterogeneity. Tumour cell movement was described by coupling movement to tissue stress, leading to a mechanically coupled (MC) RD model. Intra-tumour heterogeneity was described by including a voxel-specific carrying capacity (CC) to the RD model. The MC and CC models were also combined in a third model. To evaluate these models, rats (n = 14) with C6 gliomas were imaged with diffusion-weighted magnetic resonance imaging over 10 days to estimate tumour cellularity. Model parameters were estimated from the first three imaging time points and then used to predict tumour growth at the remaining time points which were then directly compared to experimental data. The results in this work demonstrate that mechanical-biological effects are a necessary component of brain tissue tumour modelling efforts. The results are suggestive that a variable tissue carrying capacity is a needed model component to capture tumour heterogeneity. Lastly, the results advocate the need for additional effort towards capturing tumour-to-tissue infiltration.

Recapitulating in vivo-like plasticity of glioma cell invasion along blood vessels and in astrocyte-rich stroma

Diffuse invasion of glioma cells into the brain parenchyma leads to nonresectable brain tumors and poor prognosis of glioma disease. In vivo, glioma cells can adopt a range of invasion strategies and routes, by moving as single cells, collective strands and multicellular networks along perivascular, perineuronal and interstitial guidance cues. Current in vitro assays to probe glioma cell invasion, however, are limited in recapitulating the modes and adaptability of glioma invasion observed in brain parenchyma, including collective behaviours. To mimic in vivo-like glioma cell invasion in vitro, we here applied three tissue-inspired 3D environments combining multicellular glioma spheroids and reconstituted microanatomic features of vascular and interstitial brain structures. Radial migration from multicellular glioma spheroids of human cell lines and patient-derived xenograft cells was monitored using (1) reconstituted basement membrane/hyaluronan interfaces representing the space along brain vessels; (2) 3D scaffolds generated by multi-layered mouse astrocytes to reflect brain interstitium; and (3) freshly isolated mouse brain slice culture ex vivo. The invasion patterns in vitro were validated using histological analysis of brain sections from glioblastoma patients and glioma xenografts infiltrating the mouse brain. Each 3D assay recapitulated distinct aspects of major glioma invasion patterns identified in mouse xenografts and patient brain samples, including individually migrating cells, collective strands extending along blood vessels, and multicellular networks of interconnected glioma cells infiltrating the neuropil. In conjunction, these organotypic assays enable a range of invasion modes used by glioma cells and will be applicable for mechanistic analysis and targeting of glioma cell dissemination.

Exploration of Dynamic Elastic Modulus Changes on Glioblastoma Cell Populations with Aberrant EGFR Expression as a Potential Therapeutic Intervention Using a Tunable Hyaluronic Acid Hydrogel Platform

Glioblastoma (GBM) is one of most aggressive forms of brain cancer, with a median survival time of 14.6 months following diagnosis. This low survival rate could in part be attributed to the lack of model systems of this type of cancer that faithfully recapitulate the tumor architecture and microenvironment seen in vivo in humans. Therapeutic studies would provide results that could be translated to the clinic efficiently. Here, we assess the role of the tumor microenvironment physical parameters on the tumor, and its potential use as a biomarker using a hyaluronic acid hydrogel system capable of elastic modulus tuning and dynamic elastic moduli changes. Experiments were conducted to assess the sensitivity of glioblastoma cell populations with different mutations to varying elastic moduli. Cells with aberrant epithelial growth factor receptor (EGFR) expression have a predilection for a stiffer environment, sensing these parameters through focal adhesion kinase (FAK). Importantly, the inhibition of FAK or EGFR generally resulted in reversed elastic modulus preference. Lastly, we explore the concept of therapeutically targeting the elastic modulus and dynamically reducing it via chemical or enzymatic degradation, both showing the capability to reduce or stunt proliferation rates of these GBM populations.

Cadherin- and Rigidity-Dependent Growth of Lung Cancer Cells in a Partially Confined Microenvironment

During tumor development, cancer cells constantly confront different types of extracellular barriers. However, fundamental questions like whether tumor cells will continue to grow against confinement or away from it and what key factors govern this process remain poorly understood. To address these issues, here we examined the growth dynamics of human lung epithelial carcinoma A549 cells partially confined in micrometer-sized cylindrical pores with precisely controlled wall stiffness. It was found that, after reaching confluency, the cell monolayer enclosed by a compliant wall was able to keep growing and pushing the boundary, eventually leading to a markedly enlarged pore. In contrast, a much reduced in-plane growth and elevated strain level among cells were observed when the confining wall becomes stiff. Furthermore, under such circumstance, cells switched their growth from within the monolayer to along the out-of-plane direction, resulting in cell stacking. We showed that these observations can be well explained by a simple model taking into account the deformability of the wall and the threshold stress for inhibiting cell growth. Interestingly, cadherins were found to play an important role in the proliferation and stress buildup within the cell monolayer by aggregating at cell-cell junctions. The stiff confinement led to an elevated expression level of cadherins. Furthermore, inhibition of N-cadherin resulted in a significantly suppressed cell growth under the same confining conditions.

Dynamics of cancerous tissue correlates with invasiveness

Two of the classical hallmarks of cancer are uncontrolled cell division and tissue invasion, which turn the disease into a systemic, life-threatening condition. Although both processes are studied, a clear correlation between cell division and motility of cancer cells has not been described previously. Here, we experimentally characterize the dynamics of invasive and non-invasive breast cancer tissues using human and murine model systems. The intrinsic tissue velocities, as well as the divergence and vorticity around a dividing cell correlate strongly with the invasive potential of the tissue, thus showing a distinct correlation between tissue dynamics and aggressiveness. We formulate a model which treats the tissue as a visco-elastic continuum. This model provides a valid reproduction of the cancerous tissue dynamics, thus, biological signaling is not needed to explain the observed tissue dynamics. The model returns the characteristic force exerted by an invading cell and reveals a strong correlation between force and invasiveness of breast cancer cells, thus pinpointing the importance of mechanics for cancer invasion.

Measurement of dynamic cell-induced 3D displacement fields in vitro for traction force optical coherence microscopy

Traction force microscopy (TFM) is a method used to study the forces exerted by cells as they sense and interact with their environment. Cell forces play a role in processes that take place over a wide range of spatiotemporal scales, and so it is desirable that TFM makes use of imaging modalities that can effectively capture the dynamics associated with these processes. To date, confocal microscopy has been the imaging modality of choice to perform TFM in 3D settings, although multiple factors limit its spatiotemporal coverage. We propose traction force optical coherence microscopy (TF-OCM) as a novel technique that may offer enhanced spatial coverage and temporal sampling compared to current methods used for volumetric TFM studies. Reconstructed volumetric OCM data sets were used to compute time-lapse extracellular matrix deformations resulting from cell forces in 3D culture. These matrix deformations revealed clear differences that can be attributed to the dynamic forces exerted by normal versus contractility-inhibited NIH-3T3 fibroblasts embedded within 3D Matrigel matrices. Our results are the first step toward the realization of 3D TF-OCM, and they highlight the potential use of OCM as a platform for advancing cell mechanics research.

Tensile Forces Originating from Cancer Spheroids Facilitate Tumor Invasion

The mechanical properties of tumors and the tumor environment provide important information for the progression and characterization of cancer. Tumors are surrounded by an extracellular matrix (ECM) dominated by collagen I. The geometrical and mechanical properties of the ECM play an important role for the initial step in the formation of metastasis, presented by the migration of malignant cells towards new settlements as well as the vascular and lymphatic system. The extent of this cell invasion into the ECM is a key medical marker for cancer prognosis. In vivo studies reveal an increased stiffness and different architecture of tumor tissue when compared to its healthy counterparts. The observed parallel collagen organization on the tumor border and radial arrangement at the invasion zone has raised the question about the mechanisms organizing these structures. Here we study the effect of contractile forces originated from model tumor spheroids embedded in a biomimetic collagen I matrix. We show that contractile forces act immediately after seeding and deform the ECM, thus leading to tensile radial forces within the matrix. Relaxation of this tension via cutting the collagen does reduce invasion, showing a mechanical relation between the tensile state of the ECM and invasion. In turn, these results suggest that tensile forces in the ECM facilitate invasion. Furthermore, simultaneous contraction of the ECM and tumor growth leads to the condensation and reorientation of the collagen at the spheroid’s surface. We propose a tension-based model to explain the collagen organization and the onset of invasion by forces originating from the tumor.

Lighting Up the Force: Investigating Mechanisms of Mechanotransduction Using Fluorescent Tension Probes

The ability of cells to sense the physical nature of their surroundings is critical to the survival of multicellular organisms. Cellular response to physical cues from adjacent cells and the extracellular matrix leads to a dynamic cycle in which cells respond by remodeling their local microenvironment, fine-tuning cell stiffness, polarity, and shape. Mechanical regulation is important in cellular development, normal morphogenesis, and wound healing. The mechanisms by which these finely balanced mechanotransduction events occur, however, are not well understood. In large part, this is due to the limited availability of tools to study molecular mechanotransduction events in live cells. Several classes of molecular tension probes have been recently developed which are rapidly transforming the study of mechanotransduction. Molecular tension probes are primarily based on fluorescence resonance energy transfer (FRET) and report on piconewton scale tension events in live cells. In this minireview, we describe the two main classes of tension probes, genetically encoded tension sensors and immobilized tension sensors, and discuss the advantages and limitations of each type. We discuss future opportunities to address major biological questions and outline the challenges facing the next generation of molecular tension probes.

Modeling contact guidance and invasion by cancer cells

The first step in the spread of cancer is invasion by malignant cells of the normal tissue surrounding a tumor. There is considerable evidence both in vitro and in vivo that mechanical interactions with the tissue, in particular with the biopolymer network that makes up the extracellular matrix (ECM), are important factors in invasion. The interactions take two forms: (i) contractile cells on the surface of the tumor act on the nearby ECM and remodel it; in some cases, they align the fibers of the biopolymers; (ii) the aligned fibers can enhance invasion via contact guidance, the tendency of motile cells to follow alignment. Here, we give evidence, mainly for in vitro systems, that both effects are important. We discuss how alignment occurs in biopolymers such as collagen-I (a major component of the ECM). We propose a modeling framework for computing alignment and propose phenomenologic models for contact guidance. See all articles in this Cancer Research section, "Physics in Cancer Research."

Force production and mechanical accommodation during convergent extension

Forces generated within the embryo during convergent extension (CE) must overcome mechanical resistance to push the head away from the rear. As mechanical resistance increases more than eightfold during CE and can vary twofold from individual to individual, we have proposed that developmental programs must include mechanical accommodation in order to maintain robust morphogenesis. To test this idea and investigate the processes that generate forces within early embryos, we developed a novel gel-based sensor to report force production as a tissue changes shape; we find that the mean stress produced by CE is 5.0±1.6 Pascal (Pa). Experiments with the gel-based force sensor resulted in three findings. (1) Force production and mechanical resistance can be coupled through myosin contractility. The coupling of these processes can be hidden unless affected tissues are challenged by physical constraints. (2) CE is mechanically adaptive; dorsal tissues can increase force production up to threefold to overcome a stiffer microenvironment. These findings demonstrate that mechanical accommodation can ensure robust morphogenetic movements against environmental and genetic variation that might otherwise perturb development and growth. (3) Force production is distributed between neural and mesodermal tissues in the dorsal isolate, and the notochord, a central structure involved in patterning vertebrate morphogenesis, is not required for force production during late gastrulation and early neurulation. Our findings suggest that genetic factors that coordinately alter force production and mechanical resistance are common during morphogenesis, and that their cryptic roles can be revealed when tissues are challenged by controlled biophysical constraints.

Quantification of collagen contraction in three-dimensional cell culture

Many different cell types including fibroblasts, smooth muscle cells, endothelial cells, and cancer cells exert traction forces on the fibrous components of the extracellular matrix. This can be observed as matrix contraction both macro- and microscopically in three-dimensional (3D) tissues models such as collagen type I gels. The quantification of local contraction at the micron scale, including its directionality and speed, in correlation with other parameters such as cell invasion, local protein or gene expression, can provide useful information to study wound healing, organism development, and cancer metastasis. In this article, we present a set of tools to quantify the flow dynamics of collagen contraction, induced by cells migrating out of a multicellular cancer spheroid into a three-dimensional (3D) collagen matrix. We adapted a pseudo-speckle technique that can be applied to bright-field and fluorescent microscopy time series. The image analysis presented here is based on an in-house written software developed in the Matlab (Mathworks) programming environment. The analysis program is freely available from GitHub following the link: http://dx.doi.org/10.5281/zenodo.10116. This tool provides an automatized technique to measure collagen contraction that can be utilized in different 3D cellular systems.

Engineering strategies to mimic the glioblastoma microenvironment

Glioblastoma multiforme (GBM) is the most common and deadly brain tumor, with a mean survival time of only 21months. Despite the dramatic improvements in our understanding of GBM fueled by recent revolutions in molecular and systems biology, treatment advances for GBM have progressed inadequately slowly, which is due in part to the wide cellular and molecular heterogeneity both across tumors and within a single tumor. Thus, there is increasing clinical interest in targeting cell-extrinsic factors as way of slowing or halting the progression of GBM. These cell-extrinsic factors, collectively termed the microenvironment, include the extracellular matrix, blood vessels, stromal cells that surround tumor cells, and all associated soluble and scaffold-bound signals. In this review, we will first describe the regulation of GBM tumors by these microenvironmental factors. Next, we will discuss the various in vitro approaches that have been exploited to recapitulate and model the GBM tumor microenvironment in vitro. We conclude by identifying future challenges and opportunities in this field, including the development of microenvironmental platforms amenable to high-throughput discovery and screening. We anticipate that these ongoing efforts will prove to be valuable both as enabling tools for accelerating our understanding of microenvironmental regulation in GBM and as foundations for next-generation molecular screening platforms that may serve as a conceptual bridge between traditional reductionist systems and animal or clinical studies.

Longitudinal measurement of extracellular matrix rigidity in 3D tumor models using particle-tracking microrheology

The mechanical microenvironment has been shown to act as a crucial regulator of tumor growth behavior and signaling, which is itself remodeled and modified as part of a set of complex, two-way mechanosensitive interactions. While the development of biologically-relevant 3D tumor models have facilitated mechanistic studies on the impact of matrix rheology on tumor growth, the inverse problem of mapping changes in the mechanical environment induced by tumors remains challenging. Here, we describe the implementation of particle-tracking microrheology (PTM) in conjunction with 3D models of pancreatic cancer as part of a robust and viable approach for longitudinally monitoring physical changes in the tumor microenvironment, in situ. The methodology described here integrates a system of preparing in vitro 3D models embedded in a model extracellular matrix (ECM) scaffold of Type I collagen with fluorescently labeled probes uniformly distributed for position- and time-dependent microrheology measurements throughout the specimen. In vitro tumors are plated and probed in parallel conditions using multiwell imaging plates. Drawing on established methods, videos of tracer probe movements are transformed via the Generalized Stokes Einstein Relation (GSER) to report the complex frequency-dependent viscoelastic shear modulus, G*(ω). Because this approach is imaging-based, mechanical characterization is also mapped onto large transmitted-light spatial fields to simultaneously report qualitative changes in 3D tumor size and phenotype. Representative results showing contrasting mechanical response in sub-regions associated with localized invasion-induced matrix degradation as well as system calibration, validation data are presented. Undesirable outcomes from common experimental errors and troubleshooting of these issues are also presented. The 96-well 3D culture plating format implemented in this protocol is conducive to correlation of microrheology measurements with therapeutic screening assays or molecular imaging to gain new insights into impact of treatments or biochemical stimuli on the mechanical microenvironment.

Alignment localization in nonlinear biological media

Cells imbedded in biopolymer gels are important components of tissue engineering models and cancer tumor microenvironments. In both these cases, contraction of cells attached to the gel is an important phenomenon, and the nonlinear nature of most biopolymers (such as collagen) makes understanding the mechanics of the contraction a challenging problem. Here, we investigate a unique feature of such systems: a point source of contraction leads to substantial deformation of the environment, but large strains and large alignment of the fibers of the gel are confined to a small region surrounding the source. For fibroblasts in collagen-I, we estimate that the radius of this region is of order 90 μ. We investigate this idea using continuum estimates and a finite element code, and we point out experimental manifestations of the effect.

Toward 3D biomimetic models to understand the behavior of glioblastoma multiforme cells

Glioblastoma multiforme (GBM) tumors are one of the most deadly forms of human cancer and despite improved treatments, median survival time for the majority of patients is a dismal 12-15 months. A hallmark of these aggressive tumors is their unique ability to diffusively infiltrate normal brain tissue. To understand this behavior and successfully target the mechanisms underlying tumor progression, it is crucial to develop robust experimental ex vivo disease models. This review discusses current two-dimensional (2D) experimental models, as well as animal-based models used to examine GBM cell migration, including their advantages and disadvantages. Recent attempts to develop three-dimensional (3D) tissue engineering-inspired models and their utility in unraveling the role of microenvironment on tumor cell behaviors are also highlighted. Further, the use of 3D models to bridge the gap between 2D and animal models is explored. Finally, the broad utility of such models in the context of brain cancer research is examined.

Glioblastoma behaviors in three-dimensional collagen-hyaluronan composite hydrogels

Glioblastoma multiforme (GBM) tumors, which arise from glia in the central nervous system (CNS), are one of the most deadly forms of human cancer with a median survival time of ∼1 year. Their high infiltrative capacity makes them extremely difficult to treat, and even with aggressive multimodal clinical therapies, outcomes are dismal. To improve understanding of cell migration in these tumors, three-dimensional (3D) multicomponent composite hydrogels consisting of collagen and hyaluronic acid, or hyaluronan (HA), were developed. Collagen is a component of blood vessels known to be associated with GBM migration; whereas, HA is one of the major components of the native brain extracellular matrix (ECM). We characterized hydrogel microstructural features and utilized these materials to investigate patient tumor-derived, single cell morphology, spreading, and migration in 3D culture. GBM morphology was influenced by collagen type with cells adopting a rounded morphology in collagen-IV versus a spindle-shaped morphology in collagen-I/III. GBM spreading and migration were inversely dependent on HA concentration; with higher concentrations promoting little or no migration. Further, noncancerous astrocytes primarily displayed rounded morphologies at lower concentrations of HA; in contrast to the spindle-shaped (spread) morphologies of GBMs. These results suggest that GBM behaviors are sensitive to ECM mimetic materials in 3D and that these composite hydrogels could be used to develop 3D brain mimetic models for studying migration processes.

Vaccine strategies for glioblastoma: progress and future directions

Recent advances in glioblastoma therapy have led to optimism that more effective therapies will improve outcomes. Immunotherapy is a promising approach that has demonstrated the potential to eradicate cancer cells with cellular-level accuracy while minimizing damage to surrounding healthy tissue. Several vaccination strategies have been evaluated for activity against glioblastoma in clinical trials. These include peptide vaccines, polyvalent dendritic cell vaccines, heat shock protein vaccines and adoptive immunotherapy. In this review, we highlight clinical trials representative of each of these approaches and discuss strategies for integrating these therapies into routine patient care.

Evolution and morphology of microenvironment-enhanced malignancy of three-dimensional invasive solid tumors

The emergence of invasive and metastatic behavior in malignant tumors can often lead to fatal outcomes for patients. The collective malignant tumor behavior resulting from the complex tumor-host interactions and the interactions between the tumor cells is currently poorly understood. In this paper, we employ a cellular automaton (CA) model to investigate microenvironment-enhanced malignant behaviors and morphologies of in vitro avascular invasive solid tumors in three dimensions. Our CA model incorporates a variety of microscopic-scale tumor-host interactions, including the degradation of the extracellular matrix by the malignant cells, nutrient-driven cell migration, pressure buildup due to the deformation of the microenvironment by the growing tumor, and its effect on the local tumor-host interface stability. Moreover, the effects of cell-cell adhesion on tumor growth are explicitly taken into account. Specifically, we find that while strong cell-cell adhesion can suppress the invasive behavior of the tumors growing in soft microenvironments, cancer malignancy can be significantly enhanced by harsh microenvironmental conditions, such as exposure to high pressure levels. We infer from the simulation results a qualitative phase diagram that characterizes the expected malignant behavior of invasive solid tumors in terms of two competing malignancy effects: the rigidity of the microenvironment and cell-cell adhesion. This diagram exhibits phase transitions between noninvasive and invasive behaviors. We also discuss the implications of our results for the diagnosis, prognosis, and treatment of malignant tumors.

Mimicking white matter tract topography using core-shell electrospun nanofibers to examine migration of malignant brain tumors

Glioblastoma multiforme (GBM), one of the deadliest forms of human cancer, is characterized by its high infiltration capacity, partially regulated by the neural extracellular matrix (ECM). A major limitation in developing effective treatments is the lack of in vitro models that mimic features of GBM migration highways. Ideally, these models would permit tunable control of mechanics and chemistry to allow the unique role of each of these components to be examined. To address this need, we developed aligned nanofiber biomaterials via core-shell electrospinning that permit systematic study of mechanical and chemical influences on cell adhesion and migration. These models mimic the topography of white matter tracts, a major GBM migration 'highway'. To independently investigate the influence of chemistry and mechanics on GBM behaviors, nanofiber mechanics were modulated by using different polymers (i.e., gelatin, poly(ethersulfone), poly(dimethylsiloxane)) in the 'core' while employing a common poly(ε-caprolactone) (PCL) 'shell' to conserve surface chemistry. These materials revealed GBM sensitivity to nanofiber mechanics, with single cell morphology (Feret diameter), migration speed, focal adhesion kinase (FAK) and myosin light chain 2 (MLC2) expression all showing a strong dependence on nanofiber modulus. Similarly, modulating nanofiber chemistry using extracellular matrix molecules (i.e., hyaluronic acid (HA), collagen, and Matrigel) in the 'shell' material with a common PCL 'core' to conserve mechanical properties revealed GBM sensitivity to HA; specifically, a negative effect on migration. This system, which mimics the topographical features of white matter tracts, should allow further examination of the complex interplay of mechanics, chemistry, and topography in regulating brain tumor behaviors.

Leading malignant cells initiate collective epithelial cell invasion in a three-dimensional heterotypic tumor spheroid model

Solid tumors consist of genetically and phenotypically diverse subpopulations of cancer cells with unique capacities for growth, differentiation, and invasion. While the molecular and microenvironmental bases for heterogeneity are increasingly appreciated, the outcomes of such intratumor heterogeneity, particularly in the context of tumor invasion and metastasis, remain poorly understood. To study heterotypic cell-cell interactions and elucidate the biological consequences of intratumor heterogeneity, we developed a tissue-engineered multicellular spheroid (MCS) co-culture model that recapitulates the cellular diversity and fully three-dimensional cell-cell and cell-matrix interactions that characterize human carcinomas. We found that "invasion-competent" malignant cells induced the collective invasion of otherwise "invasion-incompetent" epithelial cells, and that these two cell types consistently exhibited distinct leader and follower roles during invasion. Analysis of extracellular matrix (ECM) microarchitecture revealed that malignant cell invasion was accompanied by extensive ECM remodeling including matrix alignment and proteolytic track-making. Inhibition of cell contractility- and proteolysis-mediated matrix reorganization prevented leader-follower behavior and malignant cell-induced epithelial cell invasion. These results indicate that heterogeneous subpopulations within a tumor may possess specialized roles during tumor progression and suggest that complex interactions among the various subpopulations of cancer cells within a tumor may regulate critical aspects of tumor biology and affect clinical outcome.

Mapping of mechanical strains and stresses around quiescent engineered three-dimensional epithelial tissues

Understanding how physical signals guide biological processes requires qualitative and quantitative knowledge of the mechanical forces generated and sensed by cells in a physiologically realistic three-dimensional (3D) context. Here, we used computational modeling and engineered epithelial tissues of precise geometry to define the experimental parameters that are required to measure directly the mechanical stress profile of 3D tissues embedded within native type I collagen. We found that to calculate the stresses accurately in these settings, we had to account for mechanical heterogeneities within the matrix, which we visualized and quantified using confocal reflectance and atomic force microscopy. Using this technique, we were able to obtain traction forces at the epithelium-matrix interface, and to resolve and quantify patterns of mechanical stress throughout the surrounding matrix. We discovered that whereas single cells generate tension by contracting and pulling on the matrix, the contraction of multicellular tissues can also push against the matrix, causing emergent compression. Furthermore, tissue geometry defines the spatial distribution of mechanical stress across the epithelium, which communicates mechanically over distances spanning hundreds of micrometers. Spatially resolved mechanical maps can provide insight into the types and magnitudes of physical parameters that are sensed and interpreted by multicellular tissues during normal and pathological processes.

Diversity of dynamics and morphologies of invasive solid tumors

Complex tumor-host interactions can significantly affect the growth dynamics and morphologies of progressing neoplasms. The growth of a confined solid tumor induces mechanical pressure and deformation of the surrounding microenvironment, which in turn influences tumor growth. In this paper, we generalize a recently developed cellular automaton model for invasive tumor growth in heterogeneous microenvironments [Y. Jiao and S. Torquato, PLoS Comput. Biol.7, e1002314 (2011)] by incorporating the effects of pressure. Specifically, we explicitly model the pressure exerted on the growing tumor due to the deformation of the microenvironment and its effect on the local tumor-host interface instability. Both noninvasive-proliferative growth and invasive growth with individual cells that detach themselves from the primary tumor and migrate into the surrounding microenvironment are investigated. We find that while noninvasive tumors growing in "soft" homogeneous microenvironments develop almost isotropic shapes, both high pressure and host heterogeneity can strongly enhance malignant behavior, leading to finger-like protrusions of the tumor surface. Moreover, we show that individual invasive cells of an invasive tumor degrade the local extracellular matrix at the tumor-host interface, which diminishes the fingering growth of the primary tumor. The implications of our results for cancer diagnosis, prognosis and therapy are discussed.

Emergent behaviors from a cellular automaton model for invasive tumor growth in heterogeneous microenvironments

Understanding tumor invasion and metastasis is of crucial importance for both fundamental cancer research and clinical practice. In vitro experiments have established that the invasive growth of malignant tumors is characterized by the dendritic invasive branches composed of chains of tumor cells emanating from the primary tumor mass. The preponderance of previous tumor simulations focused on non-invasive (or proliferative) growth. The formation of the invasive cell chains and their interactions with the primary tumor mass and host microenvironment are not well understood. Here, we present a novel cellular automaton (CA) model that enables one to efficiently simulate invasive tumor growth in a heterogeneous host microenvironment. By taking into account a variety of microscopic-scale tumor-host interactions, including the short-range mechanical interactions between tumor cells and tumor stroma, degradation of the extracellular matrix by the invasive cells and oxygen/nutrient gradient driven cell motions, our CA model predicts a rich spectrum of growth dynamics and emergent behaviors of invasive tumors. Besides robustly reproducing the salient features of dendritic invasive growth, such as least-resistance paths of cells and intrabranch homotype attraction, we also predict nontrivial coupling between the growth dynamics of the primary tumor mass and the invasive cells. In addition, we show that the properties of the host microenvironment can significantly affect tumor morphology and growth dynamics, emphasizing the importance of understanding the tumor-host interaction. The capability of our CA model suggests that sophisticated in silico tools could eventually be utilized in clinical situations to predict neoplastic progression and propose individualized optimal treatment strategies.

The goal of the present work is to develop an efficient single-cell based cellular automaton (CA) model that enables one to investigate the growth dynamics and morphology of invasive solid tumors. Recent experiments have shown that highly malignant tumors develop dendritic branches composed of tumor cells that follow each other, which massively invade into the host microenvironment and ultimately lead to cancer metastasis. Previous theoretical/computational cancer modeling neither addressed the question of how such chain-like invasive branches form nor how they interact with the host microenvironment and the primary tumor. Our CA model, which incorporates a variety of microscopic-scale tumor-host interactions (e.g., the mechanical interactions between tumor cells and tumor stroma, degradation of the extracellular matrix by the tumor cells and oxygen/nutrient gradient driven cell motions), can robustly reproduce experimentally observed invasive tumor evolution and predict a wide spectrum of invasive tumor growth dynamics and emergent behaviors in various different heterogeneous environments. Further refinement of our CA model could eventually lead to the development of a powerful simulation tool for clinical purposes capable of predicting neoplastic progression and suggesting individualized optimal treatment strategies.

Quantitative microscopy and imaging tools for the mechanical analysis of morphogenesis

The importance of mechanical signals during embryogenesis and development, through both intercellular and extracellular signals, is coming into focus. It is widely hypothesized that physical forces help to guide the shape, cellular differentiation and the patterning of tissues. To test these ideas many classical engineering principles and imaging technologies are being adapted. Recent advances in microscopy, mechanical testing and genetic and pharmacological techniques, alongside computational models are helping to dissect the activity of mechanical signals in development at the cellular and molecular level. These inroads are providing maps of mechanical changes in tissue structure and stiffness, and will permit deeper insights into the role of mechanics in both developmental biology and disease.

Influence of chondroitin sulfate and hyaluronic acid on structure, mechanical properties, and glioma invasion of collagen I gels

To mimic the extracellular matrix surrounding high grade gliomas, composite matrices composed of either acid-solubilized (AS) or pepsin-treated (PT) collagen and the glycosaminoglycans chondroitin sulfate (CS) and hyaluronic acid (HA) are prepared and characterized. The structure and mechanical properties of collagen/CS and collagen/HA gels are studied via confocal reflectance microscopy (CRM) and rheology. CRM reveals that CS induces fibril bundling and increased mesh size in AS collagen but not PT collagen networks. The presence of CS also induces more substantial changes in the storage and loss moduli of AS gels than of PT gels, in accordance with expectation based on network structural parameters. The presence of HA significantly reduces mesh size in AS collagen but has a smaller effect on PT collagen networks. However, both AS and PT collagen network viscoelasticity is strongly affected by the presence of HA. The effects of CS and HA on glioma invasion is then studied in collagen/GAG matrices with network structure both similar to (PT collagen-based gels) and disparate from (AS collagen-based gels) those of the corresponding pure collagen matrices. It is shown that CS inhibits and HA has no significant effect on glioma invasion in 1.0 mg/ml collagen matrices over 3 days. The inhibitory effect of CS on glioma invasion is more apparent in AS than in PT collagen gels, suggesting invasive behavior in these environments is affected by both biochemical and network morphological changes induced by GAGs. This study is among the few efforts to differentiate structural, mechanical and biochemical effects of changes to matrix composition on cell motility in 3D.

Perspectives on biological growth and remodeling

The continuum mechanical treatment of biological growth and remodeling has attracted considerable attention over the past fifteen years. Many aspects of these problems are now well-understood, yet there remain areas in need of significant development from the standpoint of experiments, theory, and computation. In this perspective paper we review the state of the field and highlight open questions, challenges, and avenues for further development.

Gene expression profiles in 3D tumor analogs indicate compressive strain differentially enhances metastatic potential

Non-physiological mechanobiological stimuli typically occur in tumors and are considered to promote cancer spreading. Non-fluid related pressure (solid stress), which arises as tumors grow against adjacent tissues, is among the least studied endogenous stimuli due to challenges in replicating the in vivo environment. To this end, the novel devices well-pressor and the videomicroscopy-compatible optic-pressor were developed to exert precise compressive strain on cells in 3D gels in absence of other mechanical stimuli and soluble gradients. Glioblastoma (U87, HGL21) and breast cancer (MDA-MB-231) cells in 1% agarose hydrogels were exposed to 50% compressive strain for 3 h (0.25-0.05 kPa). Live imaging showed that cells elongate and deflect vertically to the load. This stimulation is shown for the first time to differentially regulate metastasis-associated genes. Furthermore, a group of differentially expressed genes was identified in all cell types, both by microarrays and confirmed by RT-PCR for select genes (caveolin-1, integrin-β1, Rac1), indicating shared response mechanisms. These genes are functionally linked and involved in decreasing cell-cell contact, increasing ECM degradation, and ultimately promoting invasion. Caveolin could orchestrate these responses while the uPA and PI3K/Akt systems could play major roles. Future work will focus on specific molecular partnerships under compression and their impact on cancer progression.

In silico estimates of the free energy rates in growing tumor spheroids

The physics of solid tumor growth can be considered at three distinct size scales: the tumor scale, the cell-extracellular matrix (ECM) scale and the sub-cellular scale. In this paper we consider the tumor scale in the interest of eventually developing a system-level understanding of the progression of cancer. At this scale, cell populations and chemical species are best treated as concentration fields that vary with time and space. The cells have chemo-mechanical interactions with each other and with the ECM, consume glucose and oxygen that are transported through the tumor, and create chemical by-products. We present a continuum mathematical model for the biochemical dynamics and mechanics that govern tumor growth. The biochemical dynamics and mechanics also engender free energy changes that serve as universal measures for comparison of these processes. Within our mathematical framework we therefore consider the free energy inequality, which arises from the first and second laws of thermodynamics. With the model we compute preliminary estimates of the free energy rates of a growing tumor in its pre-vascular stage by using currently available data from single cells and multicellular tumor spheroids.