A Model of Lung Tumor Angiogenesis in a Biomimetic Poly(ethylene glycol)-based Hydrogel System

Tumor angiogenesis is critical to tumor growth and metastasis, yet much is unknown about the role vascular cells play in the tumor microenvironment. A major outstanding challenge associated with studying tumor angiogenesis is that existing preclinical models are limited in their recapitulation of in vivo cellular organization in 3D. This disparity highlights the need for better approaches to study the dynamic interplay of relevant cells and signaling molecules as they are organized in the tumor microenvironment. In this thesis, we combined 3D culture of lung adenocarcinoma cells with adjacent 3D microvascular cell culture in 2-layer cell-adhesive, proteolytically-degradable poly(ethylene glycol) (PEG)-based hydrogels to study tumor angiogenesis and the impacts of neovascularization on tumor cell behavior.

In initial studies, 344SQ cells, a highly metastatic, murine lung adenocarcinoma cell line, were characterized alone in 3D in PEG hydrogels. 344SQ cells formed spheroids in 3D culture and secreted proangiogenic growth factors into the conditioned media that significantly increased with exposure to transforming growth factor beta 1 (TGF-β1), a potent tumor progression-promoting factor. Vascular cells alone in hydrogels formed tubule networks with localized activated TGF-β1. To study cancer cell-vascular cell interactions, the engineered 2-layer tumor angiogenesis model with 344SQ and vascular cell layers was employed. Large, invasive 344SQ clusters developed at the interface between the layers, and were not evident further from the interface or in control hydrogels without vascular cells. A modified model with spatially restricted 344SQ and vascular cell layers confirmed that observed 344SQ cluster morphological changes required close proximity to vascular cells. Additionally, TGF-β1 inhibition blocked endothelial cell-driven 344SQ migration.

Two other lung adenocarcinoma cell lines were also explored in the tumor angiogenesis model: primary tumor-derived metastasis-incompetent, murine 393P cells and primary tumor-derived metastasis-capable human A549 cells. These lung cancer cells also formed spheroids in 3D culture and secreted proangiogenic growth factors into the conditioned media. Epithelial morphogenesis varied for the primary tumor-derived cell lines compared to 344SQ cells, with far less epithelial organization present in A549 spheroids. Additionally, 344SQ cells secreted the highest concentration of two of the three angiogenic growth factors assessed. This finding correlated to 344SQ exhibiting the most pronounced morphological response in the tumor angiogenesis model compared to the 393P and A549 cell lines.

Overall, this dissertation demonstrates the development of a novel 3D tumor angiogenesis model that was used to study vascular cell-cancer cell interactions in lung adenocarcinoma cell lines with varying metastatic capacities. Findings in this thesis have helped to elucidate the role of vascular cells in tumor progression and have identified differences in cancer cell behavior in vitro that correlate to metastatic capacity, thus highlighting the usefulness of this model platform for future discovery of novel tumor angiogenesis and tumor progression-promoting targets.

Mathematical Modeling of Perifusion Cell Culture Experiments

In perifusion cell cultures, the culture medium flows continuously through a chamber containing immobilized cells and the effluent is collected at the end. In our main applications, gonadotropin releasing hormone (GnRH) or oxytocin is introduced into the chamber as the input. They stimulate the cells to secrete luteinizing hormone (LH), which is collected in the effluent. To relate the effluent LH concentration to the cellular processes producing it, we develop and analyze a mathematical model consisting of coupled partial differential equations describing the intracellular signaling and the movement of substances in the cell chamber. We analyze three different data sets and give cellular mechanisms that explain the data. Our model indicates that two negative feedback loops, one fast and one slow, are needed to explain the data and we give their biological bases. We demonstrate that different LH outcomes in oxytocin and GnRH stimulations might originate from different receptor dynamics. We analyze the model to understand the influence of parameters, like the rate of the medium flow or the fraction collection time, on the experimental outcomes. We investigate how the rate of binding and dissociation of the input hormone to and from its receptor influence its movement down the chamber. Finally, we formulate and analyze simpler models that allow us to predict the distortion of a square pulse due to hormone-receptor interactions and to estimate parameters using perifusion data. We show that in the limit of high binding and dissociation the square pulse moves as a diffusing Gaussian and in this limit the biological parameters can be estimated.