Laminin I mediates resistance to lapatinib in HER2-positive brain metastatic breast cancer cells in vitro.

The role of extracellular matrix (ECM) prevalent in the brain metastatic breast cancer (BMBC) niche in mediating cancer cell growth, survival, and response to therapeutic agents is not well understood. Emerging evidence suggests a vital role of ECM of the primary breast tumor microenvironment (TME) in tumor progression and survival. Possibly, the BMBC cells are also similarly influenced by the ECM of the metastatic niche; therefore, understanding the effect of the metastatic ECM on BMBC cells is imperative. Herein, we assessed the impact of various ECM components (i.e., Tenascin C, Laminin I, Collagen I, Collagen IV, and Fibronectin) on brain metastatic human epidermal growth factor receptor 2 (HER2)-positive and triple negative breast cancer (TNBC) cell lines in vitro. The highly aggressive TNBC cell line was minimally affected by ECM components exhibiting no remarkable changes in viability and morphology. On the contrary, amongst various ECM components tested, the HER2-positive cell line was significantly affected by Laminin I with higher viability and demonstrated a distinct spread morphology. In addition, HER2-positive BMBC cells exhibited resistance to Lapatinib in presence of Laminin I. Mechanistically, Laminin I-induced resistance to Lapatinib was mediated in part by phosphorylation of Erk 1/2 and elevated levels of Vimentin. Laminin I also significantly enhanced the migratory potential and replicative viability of HER2-positive BMBC cells. In sum, our findings show that presence of Laminin I in the TME of BMBC cells imparts resistance to targeted therapeutic agent Lapatinib, while increasing the possibility of its dispersal and clonogenic survival.

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.

Protocol for generating dormant human brain metastatic breast cancer spheroids in vitro.

Here, we present a protocol to generate dormant brain metastatic breast cancer (BMBC) spheroids utilizing hyaluronic acid (HA) hydrogels. We describe the steps for construction of spheroids from human BMBC cell lines MDA-MB-231Br and BT474Br3, HA hydrogel preparation, and spheroid plating on HA hydrogels and in suspension culture. We then detail the impact of HA hydrogel on the dormant phenotype of spheroids by measuring spheroid cross-sectional area, cell numbers, and EdU staining. For complete details on the use and execution of this protocol, please refer to Kondapaneni et al.1.

Varying the RGD concentration on a hyaluronic acid hydrogel influences dormancy versus proliferation in brain metastatic breast cancer cells.

A majority of breast cancer deaths occur due to metastasis of cancer cells to distant organs. In particular, brain metastasis is very aggressive with an extremely low survival rate. Breast cancer cells that metastasize to the brain can enter a state of dormancy, which allows them to evade death. The brain microenvironment provides biophysical, biochemical, and cellular cues, and plays an important role in determining the fate of dormant cancer cells. However, how these cues influence dormancy remains poorly understood. Herein, we employed hyaluronic acid (HA) hydrogels with a stiffness of ~0.4 kPa as an in vitro biomimetic platform to investigate the impact of biochemical cues, specifically alterations in RGD concentration, on dormancy versus proliferation in MDA-MB-231Br brain metastatic breast cancer cells. We applied varying concentrations of RGD peptide (0, 1, 2, or 4 mg/mL) to HA hydrogel surfaces and confirmed varying degrees of surface functionalization using a fluorescently labeled RGD peptide. Post functionalization, ~10,000 MDA-MB-231Br cells were seeded on top of the hydrogels and cultured for 5 days. We found that an increase in RGD concentration led to changes in cell morphology, with cells transitioning from a rounded to spindle-like morphology as well as an increase in cell spreading area. Also, an increase in RGD concentration resulted in an increase in cell proliferation. Cellular dormancy was assessed using the ratio of phosphorylated extracellular signal-regulated kinase 1/2 (p-ERK) to phosphorylated p38 (p-p38) positivity, which was significantly lower in hydrogels without RGD and in hydrogels with lowest RGD concentration compared to hydrogels functionalized with higher RGD concentration. We also demonstrated that the HA hydrogel-induced cellular dormancy was reversible. Finally, we demonstrated the involvement of ß1 integrin in mediating cell phenotype in our hydrogel platform. Overall, our results provide insight into the role of biochemical cues in regulating dormancy versus proliferation in brain metastatic breast cancer cells.

A Biomimetic Hyaluronic Acid Hydrogel Models Mass Dormancy in Brain Metastatic Breast Cancer Spheroids.

Approximately 90% of breast cancer related mortalities are due to metastasis to distant organs. At the metastatic sites, cancer cells are capable of evading death by exhibiting cellular or mass dormancy. However, the mechanisms involved in attaining dormancy at the metastatic site are not well understood. This is partly due to the lack of experimental models to study metastatic site-specific interactions, particularly in the context of brain metastatic breast cancer (BMBC). Herein, an in vitro hyaluronic acid (HA) hydrogel-based model is developed to study mass dormancy in BMBC. HA hydrogels with a stiffness of ≈0.4 kPa are utilized to mimic the brain extracellular matrix. MDA-MB-231Br or BT474Br3 BMBC spheroids are prepared and cultured on top of HA hydrogels or in suspension for 7 days. HA hydrogel induced a near mass dormant state in spheroids by achieving a balance between proliferating and dead cells. In contrast, these spheroids displayed growth in suspension cultures. The ratio of %p-ERK to %p-p38 positive cells is significantly lower in HA hydrogels compared to suspension cultures. Further, it is demonstrated that hydrogel induced mass dormant state is reversible. Overall, such models provide useful tools to study dormancy in BMBC and could be employed for drug screening.

The impact of temozolomide and lonafarnib on the stemness marker expression of glioblastoma cells in multicellular spheroids.

Glioblastoma multiforme (GBM) is a highly malignant brain tumor with a poor prognosis. The GBM microenvironment is highly heterogeneous and is composed of many cell types including astrocytes and endothelial cells (ECs) along with tumor cells, which are responsible for heightened resistance to standard chemotherapeutic drugs such as Temozolomide (TMZ). Here, we investigated how drug treatments impact stemness marker expression of GBM cells in multicellular tumor spheroid (MCTS) models. Co- and tri-culture MCTS constructed using U87-MG GBM cells, astrocytes, and/or ECs were cultured for 7 days. At Day 7, 5 µM lonafarnib (LNF), 100 µM TMZ, or combination of 5 µM LNF + 100 µM TMZ was added and the MCTS were cultured for an additional 48 h. We assessed the spheroid sizes and expression of stemness markers- NESTIN, SOX2, CD133, NANOG, and OCT4- through qRT-PCR and immunostaining. Following 48 h treatment with LNF, TMZ or their combination (LNF + TMZ), the spheroid sizes decreased compared to the untreated control. We also observed that the expression of most of the stemness markers significantly increased in the LNF + TMZ treated condition as compared to the untreated condition. These results indicate that drug treatment can influence the stemness marker expression of GBM cells in MCTS models and these aspects must be considered while evaluating therapies. In future, by incorporating other relevant cell types, we can further our understanding of their crosstalk, eventually leading to the development of new therapeutic strategies.

The Impact of Astrocytes and Endothelial Cells on Glioblastoma Stemness Marker Expression in Multicellular Spheroids.

INTRODUCTION: Glioblastoma multiforme (GBM), the most common primary brain tumor in adults, is extremely malignant and lethal. GBM tumors are highly heterogenous, being comprised of cellular and matrix components, which contribute to tumor cell invasion, cancer stem cell maintenance, and drug resistance. Here, we developed a heterotypic 3D spheroid model integrating GBM cells with astrocytes and endothelial cells (ECs) to better simulate the cellular components of the tumor microenvironment and investigate their impact on the stemness marker expression of GBM cells, which has not been previously investigated. METHODS: We used U87 GBM cells, C8-D1A mouse astrocytes, and human umbilical vein ECs to construct co- and tri-culture spheroid models in low-attachment U-well plates. We characterized the expression of known stemness markers NESTIN, SOX2, CD133, NANOG, and OCT4 in these models and compared it to respective mixed monoculture spheroids (control) using qRT-PCR and immunostaining. RESULTS: We incorporated GBM cells and astrocytes/ECs in 1:1, 1:2, 1:4, and 1:9 ratio and observed spontaneous self-assembled spheroids in all coculture conditions. We observed changing spheroid size dynamics over 7 days and an increased expression in stemness markers in GBM-astrocyte and GBM-EC coculture spheroids in 1:4 and 1:9 coculture conditions, respectively. In a triculture model employing GBM cells, astrocytes, and ECs in a 1:4:9 ratio, we found an increased expression of all the stemness markers. CONCLUSIONS: We elucidated the impact of astrocytes and ECs on GBM stemness marker expression. This multicellular spheroid model may provide an important tool for investigating the crosstalk between cell types in GBM. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s12195-021-00691-y.

Hyaluronic acid-based hydrogels to study cancer cell behaviors.

Hyaluronic acid (HA) is a natural polysaccharide and a key component of the extracellular matrix (ECM) in many tissues. Therefore, HA-based biomaterials are extensively utilized to create three dimensional ECM mimics to study cell behaviors in vitro. Specifically, derivatives of HA have been commonly used to fabricate hydrogels with controllable properties. In this review, we discuss the various chemistries employed to fabricate HA-based hydrogels as a tunable matrix to mimic the cancer microenvironment and subsequently study cancer cell behaviors in vitro. These include Michael-addition reactions, photo-crosslinking, carbodiimide chemistry, and Diels-Alder chemistry. The utility of these HA-based hydrogels to examine cancer cell behaviors such as proliferation, migration, and invasion in vitro in various types of cancer are highlighted. Overall, such hydrogels provide a biomimetic material-based platform to probe cell-matrix interactions in cancer cells in vitro and study the mechanisms associated with cancer progression.

Matrix stiffness and cluster size collectively regulate dormancy versus proliferation in brain metastatic breast cancer cell clusters.

Breast cancer cells can metastasize either as single cells or as clusters to distant organs from the primary tumor site. Cell clusters have been shown to possess higher metastatic potential compared to single cells. The organ microenvironment is critical in regulating the ultimate phenotype, specifically, the dormant versus proliferative phenotypes, of these clusters. In the context of breast cancer brain metastasis (BCBM), tumor cell cluster-organ microenvironment interactions are not well understood, in part, due to the lack of suitable biomimetic in vitro models. To address this need, herein, we report a biomaterial-based model, utilizing hyaluronic acid (HA) hydrogels with varying stiffnesses to mimic the brain microenvironment. Cell spheroids were used to mimic cell clusters. Using 100-10 000 MDA-MB-231Br BCBM cells, six different sizes of cell spheroids were prepared to study the impact of cluster size on dormancy. On soft HA hydrogels (∼0.4 kPa), irrespective of spheroid size, all cell spheroids attained a dormant phenotype, whereas on stiff HA hydrogels (∼4.5 kPa), size dependent switch between the dormant and proliferative phenotypes was noted (i.e., proliferative phenotype ≥5000 cell clusters < dormant phenotype), as tested via EdU and Ki67 staining. Furthermore, we demonstrated that the matrix stiffness driven dormancy was reversible. Such biomaterial systems provide useful tools to probe cell cluster-matrix interactions in BCBM.

Three-dimensional culture models to study drug resistance in breast cancer.

Despite recent advances in breast cancer treatment, drug resistance frequently presents as a challenge, contributing to a higher risk of relapse and decreased overall survival rate. It is now generally recognized that the extracellular matrix and cellular heterogeneity of the tumor microenvironment influences the cancer cells' ultimate fate. Therefore, strategies employed to examine mechanisms of drug resistance must take microenvironmental influences, as well as genetic mutations, into account. This review discusses three-dimensional (3D) in vitro model systems which incorporate microenvironmental influences to study mechanisms of drug resistance in breast cancer. These bioengineered models include spheroid-based models, biomaterial-based models such as polymeric scaffolds and hydrogels, and microfluidic chip-based models. The advantages of these model systems over traditionally studied two-dimensional tissue culture polystyrene are examined. Additionally, the applicability of such 3D models for studying the impact of tumor microenvironment signals on drug response and/or resistance is discussed. Finally, the potential of such models for use in the development of strategies to combat drug resistance and determine the most promising treatment regimen is explored.

An in vitro hyaluronic acid hydrogel based platform to model dormancy in brain metastatic breast cancer cells.

Breast cancer cells (BCCs) can remain dormant at the metastatic site, which when revoked leads to formation of metastasis several years after the treatment of primary tumor. Particularly, awakening of dormant BCCs in the brain results in breast cancer brain metastasis (BCBrM) which marks the most advanced stage of the disease with a median survival period of ~4-16 months. However, our understanding of dormancy associated with BCBrM remains obscure, in part, due to the lack of relevant in vitro platforms to model dormancy associated with BCBrM. To address this need, we developed an in vitro hyaluronic acid (HA) hydrogel platform to model dormancy in brain metastatic BCCs via exploiting the bio-physical cues provided by HA hydrogels while bracketing the normal brain and metastatic brain malignancy relevant stiffness range. In this system, we observed that MDA-MB-231Br and BT474Br3 brain metastatic BCCs exhibited a dormant phenotype when cultured on soft (0.4 kPa) HA hydrogel compared to stiff (4.5 kPa) HA hydrogel as characterized by significantly lower EdU and Ki67 positivity. Further, we demonstrated the nuclear localization of p21 and p27 (markers associated with dormancy) in dormant MDA-MB-231Br cells contrary to their cytoplasmic localization in the proliferative population. We also demonstrated that the stiffness-based dormancy in MDA-MB-231Br cells was reversible and was, in part, mediated by focal adhesion kinases and the initial cell seeding density. Finally, RNA sequencing confirmed the dormant phenotype in MDA-MB-231Br cells. This platform could further our understanding of dormancy in BCBrM and could be adapted for anti-metastatic drug screening. STATEMENT OF SIGNIFICANCE: Our understanding of dormancy associated with BCBrM remains obscure, in part, due to the lack of relevant in vitro platforms to model dormancy associated with BCBrM. Herein, we present a HA hydrogel-based platform to model dormancy in brain metastatic BCCs while recapitulating key aspects of brain microenvironment. We demonstrated that the biophysical cues provided the HA hydrogel mediates dormancy in brain metastatic BCCs by assessing both proliferation and cell cycle arrest markers. We also established the role of focal adhesion kinases and initial cell seeding density in the stiffness-mediated dormancy in brain metastatic BCCs. Further, RNA-seq. confirmed the dormant phenotype in brain metastatic BCCs. This platform could be utilized to further our understanding of microenvironmental regulation of dormancy in BCBrM.

Microporous scaffolds loaded with immunomodulatory lentivirus to study the contribution of immune cell populations to tumor cell recruitment in vivo.

Metastases are preceded by stochastic formation of a hospitable microenvironment known as the premetastatic niche, which has been difficult to study. Herein, we employ implantable polycaprolactone scaffolds as an engineered premetastatic niche to independently investigate the role of interleukin-10 (IL10), CXCL12, and CCL2 in recruiting immune and tumor cells and impacting breast cancer cell phenotype via lentiviral overexpression. Lentivirus delivered from scaffolds in vivo achieved sustained transgene expression for 56 days. IL10 lentiviral expression, but not CXCL12 or CCL2, significantly decreased tumor cell recruitment to scaffolds in vivo. Delivery of CXCL12 enhanced CD45+ immune cell recruitment to scaffolds while delivery of IL10 reduced immune cell recruitment. CCL2 did not alter immune cell recruitment. Tumor cell phenotype was investigated using conditioned media from immunomodulated scaffolds, with CXCL12 microenvironments reducing proliferation, and IL10 microenvironments enhancing proliferation. Migration was enhanced with CCL2 and reduced with IL10-driven microenvironments. Multiple linear regression identified populations of immune cells associated with tumor cell abundance. CD45+ immune and CD8+ T cells were associated with reduced tumor cell abundance, while CD11b+Gr1+ neutrophils and CD4+ T cells were associated with enhanced tumor cell abundance. Collectively, biomaterial scaffolds provide a tool to probe the formation and function of the premetastatic niche.

High Frequency Spectral Ultrasound Imaging to Detect Metastasis in Implanted Biomaterial Scaffolds.

For most cancers, metastasis is the point at which disease is no longer curable. Earlier detection of metastasis, when it is undetectable by current clinical methods, may enable better outcomes. We have developed a biomaterial implant that recruits metastatic cancer cells in mouse models of breast cancer. Here, we investigate spectral ultrasound imaging (SUSI) as a non-invasive strategy for detecting metastasis to the implanted biomaterial scaffolds. Our results show that SUSI, which detects parameters related to tissue composition and structure, identified changes at an early time point when tumor cells were recruited to scaffolds in orthotopic breast cancer mouse models. These changes were not associated with acellular components in the scaffolds but were reflected in the cellular composition in the scaffold microenvironment, including an increase in CD31 + CD45-endothelial cell number in tumor bearing mice. In addition, we built a classification model based on changes in SUSI parameters from scaffold measurements to stratify tumor free and tumor bearing status. Combination of a linear discriminant analysis and bagged decision trees model resulted in an area under the curve of 0.92 for receiver operating characteristics analysis. With the potential for early non-invasive detection, SUSI could facilitate clinical translation of the scaffolds for monitoring metastatic disease.

Three-dimensional biomimetic hyaluronic acid hydrogels to investigate glioblastoma stem cell behaviors.

Glioblastoma multiforme (GBM) is the deadliest form of primary brain tumor. GBM tumors are highly heterogeneous, being composed of tumor cells as well as glioblastoma stem cells (GSCs) that contribute to drug resistance and tumor recurrence following treatment. To develop therapeutic strategies, an improved understanding of GSC behavior in their microenvironment is critical. Herein, we have employed three-dimensional (3D) hyaluronic acid (HA) hydrogels that allow the incorporation of brain microenvironmental cues to investigate GSC behavior. U87 cell line and patient-derived D456 cells were cultured as suspension cultures (serum-free) and adherently (in the presence of serum) and were then encapsulated in HA hydrogels. We observed that all the seeded single cells expanded and formed spheres, and the size of the spheres increased with time. Increasing the initial cell seeding density of cells influenced the sphere size distribution. Interestingly, clonal expansion of serum-free grown tumor cells in HA hydrogels was observed. Also, stemness marker expression of serum and/or serum-free grown cells was altered when cultured in HA hydrogels. Finally, we demonstrated that HA hydrogels can support long-term GSC culture (up to 60 days) with retention of stemness markers. Overall, such biomimetic culture systems could further our understanding of the microenvironmental regulation of GSC phenotypes.

Targeting Hyaluronan Interactions for Glioblastoma Stem Cell Therapy.

Even with rigorous treatments, glioblastoma multiforme (GBM) has an abysmal median survival rate, greatly due to the drug-resistant glioblastoma stem cell (GSC) population. GSCs are known to remodel their microenvironment, but the precise role of extracellular matrix components hyaluronic acid (HA) and hyaluronidases (HAases) on the GSC population is still largely unknown. Our objective was to determine how HAase can sensitize GSCs to chemotherapy drugs by disrupting the HA-CD44 signaling. GBM cell line U87-MG and patient-derived D456 cells were grown in GSC-enriching media and treated with HA or HAase. Expressions of GSC markers, HA-related genes, and drug resistance genes were measured via flow cytometry, confocal microscopy, and qRT-PCR. Proliferation after combined HAase and temozolomide (TMZ) treatment was measured via WST-8. HA supplementation promoted the expression of GSC markers and CD44 in GBM cells cultured in serum-free media. Conversely, HAase addition inhibited GSC gene expression while promoting CD44 expression. Finally, HAase sensitized GBM cells to TMZ. We propose a combined treatment of HAase and chemotherapy drugs by disrupting the stemness-promoting HA to target GSCs. This combination therapy shows promise even when temozolomide treatment alone causes resistance.

Role of Surface Chemistry in Mediating the Uptake of Ultrasmall Iron Oxide Nanoparticles by Cancer Cells.

Ultrasmall iron oxide nanoparticles (USIONPs) (<4 nm) have recently attracted significant attention because of their potential as positive T1 magnetic resonance imaging (MRI) contrast agent contrary to larger superparamagnetic iron oxide nanoparticles (>6 nm) which act as negative T2 MRI contrast agents. However, studies on the cellular uptake behavior of these nanoparticles are very limited compared to their counterpart, larger-sized superparamagnetic iron oxide nanoparticles. In particular, the effects of specific nanoparticle parameters on the cellular uptake behavior of USIONPs by various cancer cells are not available. Here, we specifically investigated the role of USIONPs' surface functionalities [tannic acid (TA) and quinic acid (QA)] in mediating cellular uptake behavior of cancer cells pertaining to primary (U87 cells) and metastatic (MDA-MB-231Br cells) brain malignancies. Here, we chose TA and QA as representative capping molecules, wherein TA coating provides a general negatively charged nontargeting surface while QA provides a tumor-targeting surface as QA and its derivatives are known to interact with selectin receptors expressed on tumor cells and tumor endothelium. We observed differential cellular uptake in the case of TA- and QA-coated USIONPs by cancer cells. Both the cell types showed significantly higher cellular uptake of QA-coated USIONPs compared to TA-coated USIONPs at 4, 24, and 72 h. Blocking studies indicated that P-selectin cell surface receptors, in part, mediated the cellular uptake of QA-coated USIONPs. Given that P-selectin is overexpressed in cancer cells, tumor microenvironment, and at the metastatic niche, QA-coated USIONPs hold potential to be utilized as a platform for tumor-targeted drug delivery and in imaging and detection of primary and metastatic tumors.

Bioengineered models to study tumor dormancy.

The onset of cancer metastasis is the defining event in cancer progression when the disease is considered lethal. The ability of metastatic cancer cells to stay dormant for extended time periods and reawaken at later stages leading to disease recurrence makes treatment of metastatic disease extremely challenging. The tumor microenvironment plays a critical role in deciding the ultimate fate of tumor cells, yet the mechanisms by which this occurs, including dormancy, is not well understood. This mini-review discusses bioengineered models inspired from tissue engineering strategies that mimic key aspects of the tumor microenvironment to study tumor dormancy. These models include biomaterial based three dimensional models, microfluidic based models, as well as bioreactor based models that incorporate relevant microenvironmental components such as extracellular matrix molecules, niche cells, or their combination to study microenvironmental regulation of tumor dormancy. Such biomimetic models provide suitable platforms to investigate the dormant niche, including cues that drive the dormant to proliferative transition in cancer cells. In addition, the potential of such model systems to advance research in the field of tumor dormancy is discussed.

A combined compression molding, heating, and leaching process for fabrication of micro-porous poly(ε-caprolactone) scaffolds.

Biomaterial scaffolds have been increasingly used for tissue engineering applications as well as three dimensional (3D) cell culture models. Herein, we report a simple procedure combining compression molding, heating, and leaching methods for the fabrication of 3D micro-porous poly(ε-caprolactone) (PCL) biomaterial scaffolds. In this procedure, PCL micro particles are mixed with NaCl of defined sizes and compression molded, followed by heating and subsequent leaching of NaCl particles. This technique eliminates the gas foaming method, which is commonly used in the fabrication of PCL scaffolds. Process and scaffold parameters (i.e., heating time, NaCl concentration, and NaCl particle size) were varied and analyzed to determine their impact on the overall scaffold structural and mechanical properties. An increase in NaCl particle size led to an increase in pore area but did not significantly impact the mechanical properties of the scaffolds. Additionally, NaCl concentration did not show a significant effect on pore area, but considerably impacted the mechanical properties, water absorption capacity and porosity of the scaffolds. Variations in the heating time did not have an effect in the pore area, porosity, water absorption capacity or mechanical properties of the scaffolds. We also demonstrated the ability of these scaffolds to support the proliferation of breast cancer cells. Overall, these results elucidated structure-property relationships in the fabricated micro-porous PCL scaffolds. Further, this procedure could be potentially scaled up for the fabrication of micro-porous PCL scaffolds.

Biomimetic models to examine microenvironmental regulation of glioblastoma stem cells.

Glioblastoma multiforme (GBM), a malignant brain tumor, is the deadliest form of human cancer with low survival rates because of its highly invasive nature. In recent years, there has been a growing appreciation for the role that glioblastoma stem cells (GSCs) play during tumorigenesis and tumor recurrence of GBM. GSCs are a specialized subset of GBM cells with stem cell-like features that contribute to tumor initiation and therapeutic resistance. Thus, to enhance therapeutic efficiency and improve survival, targeting GSCs and their microenvironmental niche appears to be a promising approach. To develop this approach, understanding GSC-microenvironment interactions is crucial. This review discusses various biomimetic model systems to understand the impact of biophysical, biochemical, and cellular microenvironmental cues on GSC behaviors. These models include two-dimensional or matrix-free environment models, engineered biomaterial-based three-dimensional models, co-culture models, and mouse and rat in vivo models. These systems have been used to study the effects of biophysical factors, modulation of signaling pathways, extracellular matrix components, and culture conditions on the GSC phenotype. The advantages and disadvantages of these model systems and their impact in the field of GSC research are discussed.