3D modelling of angiogenesis and vascular tumour growth

This paper was presented at the 3rd Micro and Nano Flows Conference (MNF2011), which was held at the Makedonia Palace Hotel, Thessaloniki in Greece. The conference was organised by Brunel University and supported by the Italian Union of Thermofluiddynamics, Aristotle University of Thessaloniki, University of Thessaly, IPEM, the Process Intensification Network, the Institution of Mechanical Engineers, the Heat Transfer Society, HEXAG - the Heat Exchange Action Group, and the Energy Institute

Multiscale modelling of vascular tumour growth in 3D: the roles of domain size & boundary condition

We investigate a three-dimensional multiscale model of vascular tumour growth, which couples blood flow, angiogenesis, vascular remodelling, nutrient/growth factor transport, movement of, and interactions between, normal and tumour cells, and nutrient-dependent cell cycle dynamics within each cell. In particular, we determine how the domain size, aspect ratio and initial vascular network influence the tumour's growth dynamics and its long-time composition. We establish whether it is possible to extrapolate simulation results obtained for small domains to larger ones, by constructing a large simulation domain from a number of identical subdomains, each subsystem initially comprising two parallel parent vessels, with associated cells and diffusible substances. We find that the subsystem is not representative of the full domain and conclude that, for this initial vessel geometry, interactions between adjacent subsystems contribute to the overall growth dynamics. We then show that extrapolation of results from a small subdomain to a larger domain can only be made if the subdomain is sufficiently large and is initialised with a sufficiently complex vascular network. Motivated by these results, we perform simulations to investigate the tumour's response to therapy and show that the probability of tumour elimination in a larger domain can be extrapolated from simulation results on a smaller domain. Finally, we demonstrate how our model may be combined with experimental data, to predict the spatio-temporal evolution of a vascular tumour

High-Throughput Analysis of Calcium Signalling Kinetics in Astrocytes Stimulated with Different Neurotransmitters

Astrocytes express a wide range of receptors for neurotransmitters and hormones that are coupled to increases in intracellular Ca2+ concentration, enabling them to detect activity in both neuronal and vascular networks. There is increasing evidence that astrocytes are able to discriminate between different Ca2+-linked stimuli, as the efficiency of some Ca2+ dependent processes – notably release of gliotransmitters – depends on the stimulus that initiates the Ca2+ signal. The spatiotemporal complexity of Ca2+ signals is substantial, and we here tested the hypothesis that variation in the kinetics of Ca2+ responses could offer a means of selectively engaging downstream targets, if agonists exhibited a “signature shape” in evoked Ca2+ response. To test this, astrocytes were exposed to three different receptor agonists (ATP, glutamate and histamine) and the resultant Ca2+ signals were analysed for systematic differences in kinetics that depended on the initiating stimulus. We found substantial heterogeneity between cells in the time course of Ca2+ responses, but the variation did not correlate with the type or concentration of the stimulus. Using a simple metric to quantify the extent of difference between populations, it was found that the variation between agonists was insufficient to allow signal discrimination. We conclude that the time course of global intracellular Ca2+ signals does not offer the cells a means for distinguishing between different neurotransmitters

REST mediates resolution of HIF-dependent gene expression in prolonged hypoxia

The hypoxia-inducible factor (HIF) is a key regulator of the cellular response to hypoxia which promotes oxygen delivery and metabolic adaptation to oxygen deprivation. However, the degree and duration of HIF-1α expression in hypoxia must be carefully balanced within cells in order to avoid unwanted side effects associated with excessive activity. The expression of HIF-1α mRNA is suppressed in prolonged hypoxia, suggesting that the control of HIF1A gene transcription is tightly regulated by negative feedback mechanisms. Little is known about the resolution of the HIF-1α protein response and the suppression of HIF-1α mRNA in prolonged hypoxia. Here, we demonstrate that the Repressor Element 1-Silencing Transcription factor (REST) binds to the HIF-1α promoter in a hypoxia-dependent manner. Knockdown of REST using RNAi increases the expression of HIF-1α mRNA, protein and transcriptional activity. Furthermore REST knockdown increases glucose consumption and lactate production in a HIF-1α- (but not HIF-2α-) dependent manner. Finally, REST promotes the resolution of HIF-1α protein expression in prolonged hypoxia. In conclusion, we hypothesize that REST represses transcription of HIF-1α in prolonged hypoxia, thus contributing to the resolution of the HIF-1α response

Where It’s at Really Matters: In Situ In Vivo Vascular Endothelial Growth Factor Spatially Correlates with Electron Paramagnetic Resonance pO2 Images in Tumors of Living Mice

Purpose: Tumor microenvironments show remarkable tumor pO_{2} heterogeneity, as seen in prior EPR pO_{2} images (EPROI). pO_{2} correlation with hypoxia response proteins is frustrated by large rapid pO2 changes with position. Procedures: To overcome this limitation, biopsies stereotactically located in the EPROI were used to explore the relationship between vascular endothelial growth factor A (VEGF) concentrations in living mouse tumors and the local EPROI pO_{2}. Results: Quantitative ELISA VEGF concentrations correlated (p = 0.0068 to 0.019) with mean pO_{2}, median pO_{2}, and the fraction of voxels in the biopsy volume with pO_{2} less than 3, 6, and 10 Torr. Conclusions: This validates EPROI hypoxic fractions at the molecular level and provides a new paradigm for the assessment of the relationship, in vivo, between hypoxia and hypoxia response proteins. When translated to human subjects, this will enhance understanding of human tumor pathophysiology and cancer response to therapy

Quantitative Proteomic and Interaction Network Analysis of Cisplatin Resistance in HeLa Cells

Cisplatin along with other platinum based drugs are some of the most widely used chemotherapeutic agents. However drug resistance is a major problem for the successful chemotherapeutic treatment of cancer. Current evidence suggests that drug resistance is a multifactorial problem due to changes in the expression levels and activity of a wide number of proteins. A majority of the studies to date have quantified mRNA levels between drug resistant and drug sensitive cell lines. Unfortunately mRNA levels do not always correlate with protein expression levels due to post-transcriptional changes in protein abundance. Therefore global quantitative proteomics screens are needed to identify the protein targets that are differentially expressed in drug resistant cell lines. Here we employ a quantitative proteomics technique using stable isotope labeling with amino acids in cell culture (SILAC) coupled with mass spectrometry to quantify changes in protein levels between cisplatin resistant (HeLa/CDDP) and sensitive HeLa cells in an unbiased fashion. A total of 856 proteins were identified and quantified, with 374 displaying significantly altered expression levels between the cell lines. Expression level data was then integrated with a network of protein-protein interactions, and biological pathways to obtain a systems level view of proteome changes which occur with cisplatin resistance. Several of these proteins have been previously implicated in resistance towards platinum-based and other drugs, while many represent new potential markers or therapeutic targets

Global gene expression analysis of canine osteosarcoma stem cells reveals a novel role for COX-2 in tumour initiation

Osteosarcoma is the most common primary bone tumour of both children and dogs. It is an aggressive tumour in both species with a rapid clinical course leading ultimately to metastasis. In dogs and children distant metastasis occurs in >80% of individuals treated by surgery alone. Both canine and human osteosarcoma has been shown to contain a sub-population of cancer stem cells (CSCs), which may drive tumour growth, recurrence and metastasis, suggesting that naturally occurring canine osteosarcoma could act as a preclinical model for the human disease. Here we report the successful isolation of CSCs from primary canine osteosarcoma, as well as established cell lines. We show that these cells can form tumourspheres, and demonstrate relative resistance to chemotherapy. We demonstrate similar results for the human osteosarcma cell lines, U2OS and SAOS2. Utilizing the Affymetrix canine microarray, we are able to definitively show that there are significant differences in global gene expression profiles of isolated osteosarcoma stem cells and the daughter adherent cells. We identified 13,221 significant differences (p = 0.05), and significantly, COX-2 was expressed 141-fold more in CSC spheres than daughter adherent cells. To study the role of COX-2 expression in CSCs we utilized the COX-2 inhibitors meloxicam and mavacoxib. We found that COX-2 inhibition had no effect on CSC growth, or resistance to chemotherapy. However inhibition of COX-2 in daughter cells prevented sphere formation, indicating a potential significant role for COX-2 in tumour initiation

Lactic Acidosis Triggers Starvation Response with Paradoxical Induction of TXNIP through MondoA

Although lactic acidosis is a prominent feature of solid tumors, we still have limited understanding of the mechanisms by which lactic acidosis influences metabolic phenotypes of cancer cells. We compared global transcriptional responses of breast cancer cells in response to three distinct tumor microenvironmental stresses: lactic acidosis, glucose deprivation, and hypoxia. We found that lactic acidosis and glucose deprivation trigger highly similar transcriptional responses, each inducing features of starvation response. In contrast to their comparable effects on gene expression, lactic acidosis and glucose deprivation have opposing effects on glucose uptake. This divergence of metabolic responses in the context of highly similar transcriptional responses allows the identification of a small subset of genes that are regulated in opposite directions by these two conditions. Among these selected genes, TXNIP and its paralogue ARRDC4 are both induced under lactic acidosis and repressed with glucose deprivation. This induction of TXNIP under lactic acidosis is caused by the activation of the glucose-sensing helix-loop-helix transcriptional complex MondoA:Mlx, which is usually triggered upon glucose exposure. Therefore, the upregulation of TXNIP significantly contributes to inhibition of tumor glycolytic phenotypes under lactic acidosis. Expression levels of TXNIP and ARRDC4 in human cancers are also highly correlated with predicted lactic acidosis pathway activities and associated with favorable clinical outcomes. Lactic acidosis triggers features of starvation response while activating the glucose-sensing MondoA-TXNIP pathways and contributing to the “anti-Warburg” metabolic effects and anti-tumor properties of cancer cells. These results stem from integrative analysis of transcriptome and metabolic response data under various tumor microenvironmental stresses and open new paths to explore how these stresses influence phenotypic and metabolic adaptations in human cancers