A dynamic model of the hypoxia-inducible factor 1α (HIF-1α) network

Integration of superoxide formation and cristae morphology for mitochondrial redox signaling

The mitochondrial network provides the central cell's energetic and regulatory unit, which besides ATP and metabolite production participates in cellular signaling through regulated reactive oxygen species (ROS) production and various protein/ion fluxes. The inner membrane forms extensive folds, called cristae, i.e. cavities enfolded from and situated perpendicularly to its inner boundary membrane portion, which encompasses an inner cylinder within the outer membrane tubule. Mitochondrial cristae ultramorphology reflects various metabolic, physiological or pathological states. Since the mitochondrion is typically a predominant superoxide source and generated ROS also serve for the creation of information redox signals, we review known relationships between ROS generation within the respiratory chain complexes of cristae and cristae morphology. Notably, it is emphasized that cristae shape is governed by ATP-synthase dimers, MICOS complexes, OPA1 isoforms and the umbrella of their regulation, and also dependent on local protonmotive force (electrical potential component) in cristae. Cristae are also affected by redox-sensitive kinases/phosphatases or p66SHC. ATP-synthase dimers decrease in the inflated intracristal space, diminishing pH and hypothetically having minimal superoxide formation. Matrix-released signaling superoxide/H₂O₂ is predominantly integrated along mitochondrial tubules, whereas the diffusion of intracristal signaling ROS species is controlled by crista junctions, the widening of which enables specific retrograde redox signaling such as during hypoxic cell adaptation. Other physiological cases of H₂O₂ release from the mitochondrion include the modulation of insulin release in pancreatic β-cells, enhancement of insulin signaling in peripheral tissues, signaling by T-cell receptors, retrograde signaling during the cell cycle and cell differentiation, specifically that of adipocytes.

REST is a hypoxia-responsive transcriptional repressor

Cellular exposure to hypoxia results in altered gene expression in a range of physiologic and pathophysiologic states. Discrete cohorts of genes can be either up- or down-regulated in response to hypoxia. While the Hypoxia-Inducible Factor (HIF) is the primary driver of hypoxia-induced adaptive gene expression, less is known about the signalling mechanisms regulating hypoxia-dependent gene repression. Using RNA-seq, we demonstrate that equivalent numbers of genes are induced and repressed in human embryonic kidney (HEK293) cells. We demonstrate that nuclear localization of the Repressor Element 1-Silencing Transcription factor (REST) is induced in hypoxia and that REST is responsible for regulating approximately 20% of the hypoxia-repressed genes. Using chromatin immunoprecipitation assays we demonstrate that REST-dependent gene repression is at least in part mediated by direct binding to the promoters of target genes. Based on these data, we propose that REST is a key mediator of gene repression in hypoxia.

Unveiling Hidden Dynamics of Hippo Signalling: A Systems Analysis

The Hippo signalling pathway has recently emerged as an important regulator of cell apoptosis and proliferation with significant implications in human diseases. In mammals, the pathway contains the core kinases MST1/2, which phosphorylate and activate LATS1/2 kinases. The pro-apoptotic function of the MST/LATS signalling axis was previously linked to the Akt and ERK MAPK pathways, demonstrating that the Hippo pathway does not act alone but crosstalks with other signalling pathways to coordinate network dynamics and cellular outcomes. These crosstalks were characterised by a multitude of complex regulatory mechanisms involving competitive protein-protein interactions and phosphorylation mediated feedback loops. However, how these different mechanisms interplay in different cellular contexts to drive the context-specific network dynamics of Hippo-ERK signalling remains elusive. Using mathematical modelling and computational analysis, we uncovered that the Hippo-ERK network can generate highly diverse dynamical profiles that can be clustered into distinct dose-response patterns. For each pattern, we offered mechanistic explanation that defines when and how the observed phenomenon can arise. We demonstrated that Akt displays opposing, dose-dependent functions towards ERK, which are mediated by the balance between the Raf-1/MST2 protein interaction module and the LATS1 mediated feedback regulation. Moreover, Ras displays a multi-functional role and drives biphasic responses of both MST2 and ERK activities; which are critically governed by the competitive protein interaction between MST2 and Raf-1. Our study represents the first in-depth and systematic analysis of the Hippo-ERK network dynamics and provides a concrete foundation for future studies.

Hypoxia Stress Response Pathways: Modeling and Targeted Therapy

Hypoxia is a consequence of the decrease in the oxygen reaching the tissues of the body. It is a prominent feature of most solid tumors and is known to promote malignant progression, metastatic capacity, resistance to chemotherapy, and leads to poor patient prognosis. When a cell is under hypoxic stress, a cascade of cell signals is initiated through a family of transcription factors named as hypoxia inducible factors (HIFs). During hypoxia, HIF stabilizes and enters the nucleus and binds to the DNA via the hypoxia response element (HRE) and leads to the translation of downstream genes. The decision of adaptation or cell death depends on the extent of hypoxic stress faced by the cells. Proper understanding of hypoxic stress response is critical for understanding the mechanism of tumor cell adaptation to hypoxia and to develop efficient therapeutic interventions. In this paper, we develop a Boolean network model with targeted drug intervention in a cell that mimics persistent hypoxia. This hypoxic pathway is combined with pathways that help the cell adapt to the situation or undergo cell death. It is linked to apoptosis, cell survival, and energy production via the p53/Mdm2, PI3k/Akt/mTOR, and Glycolysis/TCA cycle pathways, respectively. In this model, we have incorporated eight known anticancer drugs that target these pathways. Through simulations, we have identified drug combinations that provided overall benefits to the cell in comparison to the no intervention case. Where applicable, the behavior predicted by this model is in agreement with experimental observations from the published literature.

Poor prognosis in human papillomavirus-positive oropharyngeal squamous cell carcinomas that overexpress hypoxia inducible factor-1α

BACKGROUND

Hypoxia induces stabilization of the transcription factor HIF-1alpha (HIF-1α), associated with (chemo-)radiotherapy resistance in oropharyngeal squamous cell carcinoma (SCC). We investigated the effect of HIF-1α expression on survival in relation to human papillomavirus (HPV) status in oropharyngeal SCC.

METHODS

We conducted an immunohistochemical analysis of HIF-1α protein expression and downstream targets carbonic anhydrase-IX (CA-IX) and glucose transporter-1 (GLUT-1) in 274 patients with oropharyngeal SCC. Overall survival (OS) was analyzed in total and stratified for HPV status and treatment.

RESULTS

In HPV-positive tumors (n = 44), HIF-1α overexpression predicted worse OS (hazard ratio [HR] = 6.23; p = .012), whereas TNM classification or treatment modality did not. In HPV-negative tumors (n = 218), advanced T and N classification and HIF-1α overexpression all independently predicted worse OS. However, the effect of HIF-1α overexpression on OS was lower in HPV-negative (HR = 1.50; p = .024) than in HPV-positive tumors.

CONCLUSION

HIF-1α overexpression is associated with worse OS and characterized a subgroup of patients with HPV-positive oropharyngeal SCC with poor prognosis. Possibly, patients with HIF-1α overexpressing HPV-positive tumors should not be eligible for treatment dose deescalation. © 2016 Wiley Periodicals, Inc. Head Neck 38: 1338-1346, 2016.

Substrate-Trapped Interactors of PHD3 and FIH Cluster in Distinct Signaling Pathways

Amino acid hydroxylation is a post-translational modification that regulates intra- and inter-molecular protein-protein interactions. The modifications are regulated by a family of 2-oxoglutarate- (2OG) dependent enzymes and, although the biochemistry is well understood, until now only a few substrates have been described for these enzymes. Using quantitative interaction proteomics, we screened for substrates of the proline hydroxylase PHD3 and the asparagine hydroxylase FIH, which regulate the HIF-mediated hypoxic response. We were able to identify hundreds of potential substrates. Enrichment analysis revealed that the potential substrates of both hydroxylases cluster in the same pathways but frequently modify different nodes of signaling networks. We confirm that two proteins identified in our screen, MAPK6 (Erk3) and RIPK4, are indeed hydroxylated in a FIH- or PHD3-dependent mechanism. We further determined that FIH-dependent hydroxylation regulates RIPK4-dependent Wnt signaling, and that PHD3-dependent hydroxylation of MAPK6 protects the protein from proteasomal degradation.

Role of hypoxia-related proteins in invasion of ameloblastoma cells: crosstalk between NOTCH1, hypoxia-inducible factor 1α, a disintegrin and metalloproteinase 12, and heparin-binding epidermal growth factor

AIMS

Ameloblastoma AME is a benign tumour characterized by local invasiveness, high recurrence rates, and diverse histological patterns. The oxygen concentration is reduced in specific areas of the tumour microenvironment, which leads to intratumoral hypoxia. Crosstalk between NOTCH1, a disintegrin and metalloproteinase 12 (ADAM-12), hypoxia-inducible factor 1α (HIF-1α) and heparin-binding epidermal growth factor (HB-EGF) under hypoxic conditions has been implicated in invadopodia formation, tumour invasiveness, and metastasis development. The aim of this study was to analyse the expression of these proteins, in order to further elucidate the mechanisms underlying AME invasiveness.

METHODS AND RESULTS

Twenty cases of AME, eight calcifying cystic odontogenic tumours CCOTs and 10 samples of dental follicle were used to investigate the expression of these proteins by immunohistochemistry with the primary antibodies anti-NOTCH1, anti-ADAM-12, anti-HIF-1α, and anti-HB-EGF. Immunostaining results were expressed as the percentage of stained area in images acquired in an AxioScope microscope equipped with an AxioCamHRc camera and a × 40 objective. The results showed that immunoexpression of all proteins was higher in the AME samples than in the CCOT and dental follicle samples (P < 0.05).

CONCLUSIONS

AME showed an increased presence of proteins associated with tumour invasiveness, which indicates a possible role of these proteins in the biological behaviour of this tumour.

Understanding complexity in the HIF signaling pathway using systems biology and mathematical modeling

Hypoxia is a common micro-environmental stress which is experienced by cells during a range of physiologic and pathophysiologic processes. The identification of the hypoxia-inducible factor (HIF) as the master regulator of the transcriptional response to hypoxia transformed our understanding of the mechanism underpinning the hypoxic response at the molecular level and identified HIF as a potentially important new therapeutic target. It has recently become clear that multiple levels of regulatory control exert influence on the HIF pathway giving the response a complex and dynamic activity profile. These include positive and negative feedback loops within the HIF pathway as well as multiple levels of crosstalk with other signaling pathways. The emerging model reflects a multi-level regulatory network that affects multiple aspects of the physiologic response to hypoxia including proliferation, apoptosis, and differentiation. Understanding the interplay between the molecular mechanisms involved in the dynamic regulation of the HIF pathway at a systems level is critically important in defining new appropriate therapeutic targets for human diseases including ischemia, cancer, and chronic inflammation. Here, we review our current knowledge of the regulatory circuits which exert influence over the HIF response and give examples of in silico model-based predictions of the dynamic behaviour of this system.

Stochastic Petri Net Modeling of Hypoxia Pathway Predicts a Novel Incoherent Feed-Forward Loop Controlling SDF-1 Expression in Acute Kidney Injury

Homing of stem cells to the sites of injury is crucial for tissue regeneration. Stromal derived factor 1 (SDF-1) is among the most important chemokines recruiting these cells. Unexpectedly, our previous experimental data on mouse models of acute kidney injury showed that SDF-1 has a declining trend following ischemic kidney insult. To describe this unforeseen observation, a stochastic Petri net model of SDF-1 regulation in the hypoxia pathway was constructed based on main related components extracted from literature. Using this strategy, predictions regarding the underlying mechanisms of SDF-1 kinetics are generated and a novel incoherent feed forward loop regulating SDF-1 expression is proposed. The computational approach suggested here can be exploited to propose novel therapies for debilitating disorders such as kidney injury.

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.

Secretome identification of immune cell factors mediating metastatic cell homing

Metastatic cell homing is a complex process mediated in part by diffusible factors secreted from immune cells found at a pre-metastatic niche. We report on connecting secretomics and TRanscriptional Activity CEll aRray (TRACER) data to identify functional paracrine interactions between immune cells and metastatic cells as novel mediators of homing. Metastatic breast cancer mouse models were used to generate a diseased splenocyte conditioned media (D-SCM) containing immune cell secreted factors. MDA-MB-231 metastatic cell activity including cell invasion, migration, transendothelial migration, and proliferation were increased in D-SCM relative to control media. Our D-SCM secretome analysis yielded 144 secreted factor candidates that contribute to increased metastatic cell activity. The functional mediators of homing were identified using MetaCore software to determine interactions between the immune cell secretome and the TRACER-identified active transcription factors within metastatic cells. Among the 5 candidate homing factors identified, haptoglobin was selected and validated in vitro and in vivo as a key mediator of homing. Our studies demonstrate a novel systems biology approach to identify functional signaling factors associated with a cellular phenotype, which provides an enabling tool that complements large-scale protein identification provided by proteomics.

Computational Model of MicroRNA Control of HIF-VEGF Pathway: Insights into the Pathophysiology of Ischemic Vascular Disease and Cancer

HRMs (hypoxia-responsive miRNAs) are a specific group of microRNAs that are regulated by hypoxia. Recent studies revealed that several HRMs including let-7 family miRNAs were highly induced in response to HIF (hypoxia-inducible factor) stabilization in hypoxia, and they potently participated in angiogenesis by targeting AGO1 (argonaute 1) and upregulating VEGF (vascular endothelial growth factor). Here we constructed a novel computational model of microRNA control of HIF-VEGF pathway in endothelial cells to quantitatively investigate the role of HRMs in modulating the cellular adaptation to hypoxia. The model parameters were optimized and the simulations based on these parameters were validated against several published in vitro experimental data. To advance the mechanistic understanding of oxygen sensing in hypoxia, we demonstrated that the rate of HIF-1α nuclear import substantially influences its stabilization and the formation of HIF-1 transcription factor complex. We described the biological feedback loops involving let-7 and AGO1 in which the impact of external perturbations were minimized; as a pair of master regulators when low oxygen tension was sensed, they coordinated the critical process of VEGF desuppression in a controlled manner. Prompted by the model-motivated discoveries, we proposed and assessed novel pathway-specific therapeutics that modulate angiogenesis by adjusting VEGF synthesis in tumor and ischemic cardiovascular disease. Through simulations that capture the complex interactions between miRNAs and miRNA-processing molecules, this model explores an innovative perspective about the distinctive yet integrated roles of different miRNAs in angiogenesis, and it will help future research to elucidate the dysregulated miRNA profiles found in cancer and various cardiovascular diseases.

Cells living in a hypoxic environment secrete signals to stimulate new blood vessel growth, a process termed angiogenesis, to acquire more oxygen and nutrients. Hypoxia-inducible factor 1 (HIF-1) accumulates in hypoxia and expedites the release of pro-angiogenic cytokines such as vascular endothelial growth factor (VEGF), a prime inducer of angiogenesis. The intermediate signaling events connecting HIF-1 and VEGF are tightly controlled by microRNAs (miRs), which are endogenous, non-coding RNA molecules and powerful regulators in cancer and cardiovascular disease. Given the importance of angiogenesis in tumor development and post-ischemia reperfusion, it holds great basic research and therapeutic value to investigate how miRs modulate intracellular VEGF synthesis to control angiogenesis in hypoxia. We present a computational model that details the interactions between miRs and other key molecules which make up different hierarchies in HIF-miR-VEGF pathway. Based on simulation analysis, new potential therapies are introduced and tested in silico, from which the strategies that most effectively reduce VEGF synthesis in cancer, or enhance VEGF release in ischemic vascular disease are identified. We conclude that in hypoxia different miRs work consonantly to fine-tune the cellular adaptations; when a master miR alters its expression, dynamics of other miRs vary accordingly which together contribute to aberrant RNA/protein profiles observed in the pathophysiology of multiple diseases.

Imaging of oxygen gradients in giant umbrella cells: an ex vivo PLIM study

O2 plays a pivotal role in aerobic metabolism and regulation of cell and tissue function. Local differences and fluctuations in tissue O2 levels are well documented; however, the physiological significance of O2 microgradients, particularly at the subcellular level, remains poorly understood. Using the cell-penetrating phosphorescent O2 probe Pt-Glc and confocal fluorescence microscopy, we visualized O2 distribution in individual giant (>100-μm) umbrella cells located superficially in the urinary bladder epithelium. We optimized conditions for in vivo phosphorescent staining of the inner surface of the mouse bladder and subsequent ex vivo analysis of excised live tissue. Imaging experiments revealed significant (≤85 μM) and heterogeneous deoxygenation within respiring umbrella cells, with radial O2 gradients of up to 40 μM across the cell, or ∼0.6 μM/μm. Deeply deoxygenated (5–15 μM O2) regions were seen to correspond to the areas enriched with polarized mitochondria. Pharmacological activation of mitochondrial respiration decreased oxygenation and O2 gradients in umbrella cells, while inhibition with antimycin A dissipated the gradients and caused gradual reoxygenation of the tissue to ambient levels. Detailed three-dimensional maps of O2 distribution potentially can be used for the modeling of intracellular O2-dependent enzymatic reactions and downstream processes, such as hypoxia-inducible factor signaling. Further ex vivo and in vivo studies on intracellular and tissue O2 gradients using confocal imaging can shed light on the molecular mechanisms regulating O2-dependent (patho)physiological processes in the bladder and other tissues.

Modeling dioxygenase enzyme kinetics in familial paraganglioma

Hypoxia inducible factors (HIFs) play vital roles in cellular maintenance of oxygen homeostasis. These transcription factors are responsible for the expression of genes involved in angiogenesis, metabolism, and cell proliferation. Here, we generate a detailed mathematical model for the enzyme kinetics of α-ketoglutarate-dependent HIF prolyl 4-hydroxylase domain (PHD) dioxygenases to simulate our in vitro data showing synergistic PHD inhibition by succinate and hypoxia in experimental models of succinate dehydrogenase loss, which phenocopy familial paraganglioma. Our mathematical model confirms the inhibitory synergy of succinate and hypoxia under physiologically-relevant conditions. In agreement with our experimental data, the model predicts that HIF1α is not stabilized under atmospheric oxygen concentrations, as observed. Further, the model confirms that addition of α-ketoglutarate can reverse PHD inhibition by succinate and hypoxia in SDH-deficient cells.

Modeling the interplay between the HIF-1 and p53 pathways in hypoxia

Both the hypoxia-inducible factor-1 (HIF-1) and tumor suppressor p53 are involved in the cellular response to hypoxia. How the two transcription factors interact to determine cell fates is less well understood. Here, we developed a network model to characterize crosstalk between the HIF-1 and p53 pathways, taking into account that HIF-1α and p53 are targeted for proteasomal degradation by Mdm2 and compete for binding to limiting co-activator p300. We reported the network dynamics under various hypoxic conditions and revealed how the stabilization and transcriptional activities of p53 and HIF-1α are modulated to determine the cell fate. We showed that both the transrepression and transactivation activities of p53 promote apoptosis induction. This work provides new insight into the mechanism for the cellular response to hypoxia.

Systems biology of the microvasculature

The vascular network carries blood throughout the body, delivering oxygen to tissues and providing a pathway for communication between distant organs. The network is hierarchical and structured, but also dynamic, especially at the smaller scales. Remodeling of the microvasculature occurs in response to local changes in oxygen, gene expression, cell-cell communication, and chemical and mechanical stimuli from the microenvironment. These local changes occur as a result of physiological processes such as growth and exercise, as well as acute and chronic diseases including stroke, cancer, and diabetes, and pharmacological intervention. While the vasculature is an important therapeutic target in many diseases, drugs designed to inhibit vascular growth have achieved only limited success, and no drug has yet been approved to promote therapeutic vascular remodeling. This highlights the challenges involved in identifying appropriate therapeutic targets in a system as complex as the vasculature. Systems biology approaches provide a means to bridge current understanding of the vascular system, from detailed signaling dynamics measured in vitro and pre-clinical animal models of vascular disease, to a more complete picture of vascular regulation in vivo. This will translate to an improved ability to identify multi-component biomarkers for diagnosis, prognosis, and monitoring of therapy that are easy to measure in vivo, as well as better drug targets for specific disease states. In this review, we summarize systems biology approaches that have advanced our understanding of vascular function and dysfunction in vivo, with a focus on computational modeling.

Clinical implications of hypoxia biomarker expression in head and neck squamous cell carcinoma: a systematic review

Awareness increases that the tumor biology influences treatment outcome and prognosis in cancer. Tumor hypoxia is thought to decrease sensitivity to radiotherapy and some forms of chemotherapy. Presence of hypoxia may be assessed by investigating expression of endogenous markers of hypoxia (EMH) using immunohistochemistry (IHC). In this systematic review we investigated the effect of EMH expression on local control and survival according to treatment modality in head and neck cancer (head and neck squamous cell carcinoma [HNSCC]). A search was performed in MEDLINE and EMBASE. Studies were eligible for inclusion that described EMH expression in relation to outcome in HNSCC patients. Quality was assessed using the Quality in Prognosis Studies (QUIPS) tool. Hazard ratios for locoregional control and survival were extracted. Forty studies of adequate quality were included. HIF-1a, HIF-2a, CA-IX, GLUT-1, and OPN were identified as the best described EMHs. With exception of HIF-2a, all EMHs were significantly related to adverse outcome in multiple studies, especially in studies where patients underwent single-modality treatment. Positive expression was often correlated with adverse clinical characteristics, including disease stage and differentiation grade. In summary, EMH expression was common in HNSCC patients and negatively influenced their prognosis. Future studies should investigate the effect of hypoxia-modified treatment schedules in patients with high In summary, EMH expression. These may include ARCON, treatment with nimorazole, or novel targeted therapies directed at hypoxic tissue. Also, the feasibility of surgical removal of the hypoxic tumor volume prior to radiotherapy should be investigated.

Nonlinear signalling networks and cell-to-cell variability transform external signals into broadly distributed or bimodal responses

We show theoretically and experimentally a mechanism behind the emergence of wide or bimodal protein distributions in biochemical networks with nonlinear input-output characteristics (the dose-response curve) and variability in protein abundance. Large cell-to-cell variation in the nonlinear dose-response characteristics can be beneficial to facilitate two distinct groups of response levels as opposed to a graded response. Under the circumstances that we quantify mathematically, the two distinct responses can coexist within a cellular population, leading to the emergence of a bimodal protein distribution. Using flow cytometry, we demonstrate the appearance of wide distributions in the hypoxia-inducible factor-mediated response network in HCT116 cells. With help of our theoretical framework, we perform a novel calculation of the magnitude of cell-to-cell heterogeneity in the dose-response obtained experimentally.

Role of compartmentalization on HiF-1α degradation dynamics during changing oxygen conditions: a computational approach

HiF-1α is the central protein driving the cellular response to hypoxia. Its accumulation in cancer cells is linked to the appearance of chemoresistant and aggressive tumor phenotypes. As a consequence, understanding the regulation of HiF-1α dynamics is a major issue to design new anti-cancer therapies. In this paper, we propose a model of the hypoxia pathway, involving HiF-1α and its inhibitor pVHL. Based on data from the literature, we made the hypothesis that the regulation of HiF-1α involves two compartments (nucleus and cytoplasm) and a constitutive shuttle of the pVHL protein between them. We first show that this model captures correctly the main features of HiF-1α dynamics, including the bi-exponential degradation profile in normoxia, the kinetics of induction in hypoxia, and the switch-like accumulation. Second, we simulated the effects of a hypoxia/reoxygenation event, and show that it generates a strong instability of HiF-1α. The protein concentration rapidly increases 3 hours after the reoxygenation, and exhibits an oscillating pattern. This effect vanishes if we do not consider compartmentalization of HiF-1α. This result can explain various counter-intuitive observations about the specific molecular and cellular response to the reoxygenation process. Third, we simulated the HiF-1α dynamics in the tumor case. We considered different types of mutations associated with tumorigenesis, and we compared their consequences on HiF-1α dynamics. Then, we tested different therapeutics strategies. We show that a therapeutic decrease of HiF-1α nuclear level is not always correlated with an attenuation of reoxygenation-induced instabilities. Thus, it appears that the design of anti-HiF-1α therapies have to take into account these two aspects to maximize their efficiency.

Differential contribution of key metabolic substrates and cellular oxygen in HIF signalling

Changes in availability and utilisation of O2 and metabolic substrates are common in ischemia and cancer. We examined effects of substrate deprivation on HIF signalling in PC12 cells exposed to different atmospheric O2. Upon 2-4h moderate hypoxia, HIF-α protein levels were dictated by the availability of glutamine and glucose, essential for deep cell deoxygenation and glycolytic ATP flux. Nuclear accumulation of HIF-1α dramatically decreased upon inhibition of glutaminolysis or glutamine deprivation. Elevation of HIF-2α levels was transcription-independent and associated with the activation of Akt and Erk1/2. Upon 2h anoxia, HIF-2α levels strongly correlated with cellular ATP, produced exclusively via glycolysis. Without glucose, HIF signalling was suppressed, giving way to other regulators of cell adaptation to energy crisis, e.g. AMPK. Consequently, viability of cells deprived of O2 and glucose decreased upon inhibition of AMPK with dorsomorphin. The capacity of cells to accumulate HIF-2α decreased after 24h glucose deprivation. This effect, associated with increased AMPKα phosphorylation, was sensitive to dorsomorphin. In chronically hypoxic cells, glutamine played no major role in HIF-2α accumulation, which became mainly glucose-dependent. Overall, the availability of O2 and metabolic substrates intricately regulates HIF signalling by affecting cell oxygenation, ATP levels and pathways involved in production of HIF-α.

Hypoxia-sensitive pathways in inflammation-driven fibrosis

Tissue injury can occur for a variety of reasons, including physical damage, infection, and ischemia. The ability of tissues to effectively recover from injury is a cornerstone of human health. The healing response in tissues is conserved across organs and typically involves distinct but overlapping inflammatory, proliferative, and maturation/resolution phases. If the inflammatory phase is not successfully controlled and appropriately resolved, an excessive healing response characterized by scar formation can lead to tissue fibrosis, a major clinical complication in disorders such as Crohn's disease (CD). As a result of enhanced metabolic and inflammatory processes during chronic inflammation, profound changes in tissue oxygen levels occur leading to localized tissue hypoxia. Therefore, inflammation, fibrosis, and hypoxia are coincidental events during inflammation-driven fibrosis. Our current understanding of the mechanism(s) underpinning fibrosis is limited as are the therapeutic options available. In this review, we discuss what is known about the cellular and molecular mechanisms underpinning inflammation-driven fibrosis and how hypoxia may play a role in shaping this process.

Computational systems biology approaches to anti-angiogenic cancer therapeutics

Angiogenesis is an exquisitely regulated process that is required for physiological processes and is also important in numerous diseases. Tumors utilize angiogenesis to generate the vascular network needed to supply the cancer cells with nutrients and oxygen, and many cancer drugs aim to inhibit tumor angiogenesis. Anti-angiogenic therapy involves inhibiting multiple cell types, molecular targets, and intracellular signaling pathways. Computational tools are useful in guiding treatment strategies, predicting the response to treatment, and identifying new targets of interest. Here, we describe progress that has been made in applying mathematical modeling and bioinformatics approaches to study anti-angiogenic therapeutics in cancer.

The ErbB4 CYT2 variant protects EGFR from ligand-induced degradation to enhance cancer cell motility

The epidermal growth factor receptor (EGFR) is a member of the ErbB family that can promote the migration and proliferation of breast cancer cells. Therapies that target EGFR can promote the dimerization of EGFR with other ErbB receptors, which is associated with the development of drug resistance. Understanding how interactions among ErbB receptors alter EGFR biology could provide avenues for improving cancer therapy. We found that EGFR interacted directly with the CYT1 and CYT2 variants of ErbB4 and the membrane-anchored intracellular domain (mICD). The CYT2 variant, but not the CYT1 variant, protected EGFR from ligand-induced degradation by competing with EGFR for binding to a complex containing the E3 ubiquitin ligase c-Cbl and the adaptor Grb2. Cultured breast cancer cells overexpressing both EGFR and ErbB4 CYT2 mICD exhibited increased migration. With molecular modeling, we identified residues involved in stabilizing the EGFR dimer. Mutation of these residues in the dimer interface destabilized the complex in cells and abrogated growth factor-stimulated cell migration. An exon array analysis of 155 breast tumors revealed that the relative mRNA abundance of the ErbB4 CYT2 variant was increased in ER+ HER2- breast cancer patients, suggesting that our findings could be clinically relevant. We propose a mechanism whereby competition for binding to c-Cbl in an ErbB signaling heterodimer promotes migration in response to a growth factor gradient.

Metabolic requirements for the maintenance of self-renewing stem cells

A distinctive feature of stem cells is their capacity to self-renew to maintain pluripotency. Studies of genetically-engineered mouse models and recent advances in metabolomic analysis, particularly in haematopoietic stem cells, have deepened our understanding of the contribution made by metabolic cues to the regulation of stem cell self-renewal. Many types of stem cells heavily rely on anaerobic glycolysis, and stem cell function is also regulated by bioenergetic signalling, the AKT-mTOR pathway, Gln metabolism and fatty acid metabolism. As maintenance of a stem cell pool requires a finely-tuned balance between self-renewal and differentiation, investigations into the molecular mechanisms and metabolic pathways underlying these decisions hold great therapeutic promise.

Tight control of hypoxia-inducible factor-α transient dynamics is essential for cell survival in hypoxia

Intracellular signaling involving hypoxia-inducible factor (HIF) controls the adaptive responses to hypoxia. There is a growing body of evidence demonstrating that intracellular signals encode temporal information. Thus, the dynamics of protein levels, as well as protein quantity and/or localization, impacts on cell fate. We hypothesized that such temporal encoding has a role in HIF signaling and cell fate decisions triggered by hypoxic conditions. Using live cell imaging in a controlled oxygen environment, we observed transient 3-h pulses of HIF-1α and -2α expression under continuous hypoxia. We postulated that the well described prolyl hydroxylase (PHD) oxygen sensors and HIF negative feedback regulators could be the origin of the pulsatile HIF dynamics. We used iterative mathematical modeling and experimental analysis to scrutinize which parameter of the PHD feedback could control HIF timing and we probed for the functional redundancy between the three main PHD proteins. We identified PHD2 as the main PHD responsible for HIF peak duration. We then demonstrated that this has important consequences, because the transient nature of the HIF pulse prevents cell death by avoiding transcription of p53-dependent pro-apoptotic genes. We have further shown the importance of considering HIF dynamics for coupling mathematical models by using a described HIF-p53 mathematical model. Our results indicate that the tight control of HIF transient dynamics has important functional consequences on the cross-talk with key signaling pathways controlling cell survival, which is likely to impact on HIF targeting strategies for hypoxia-associated diseases such as tumor progression and ischemia.

Monitoring of transcriptional dynamics of HIF and NFκB activities

Genetic experiments over the last few decades have identified many regulatory proteins critical for DNA transcription. The dynamics of their transcriptional activities shape the differential expression of the genes they control. Here we describe a simple method, based on the secreted luciferase, to measure the activities of two transcription factors NFκB and HIF. This technique can effectively monitor dynamics of transcriptional events in a population of cells and be up-scaled for high-throughput screening and promoter analysis, making it ideal for data-demanding applications such as mathematical modelling.

FoxO3a modulates hypoxia stress induced oxidative stress and apoptosis in cardiac microvascular endothelial cells

Cardiac microvascular endothelial cells (CMECs) dysfunction induced by hypoxia is an important pathophysiological event in myocardium ischemic injury, whereas, the underlying mechanism is not fully clarified. FoxO transcription factors regulate target genes involved in apoptosis and cellular reactive oxygen species (ROS) production. Therefore, the present study was designed to elucidate the potential role of FoxOs on the hypoxia-induced ROS formation and apoptosis in CMECs. Exposure to low oxygen tension stimulated ROS accumulation and increased apoptosis in CMECs within 6–24 h. Hypoxia also significantly increased the expressions of HIF-1α and FoxO3a. However, hypoxia decreased the phosphorylation of Akt and FoxO3a, correlated with increased nuclear accumulation. Conversely, the expression of FoxO1 was not significantly altered by hypoxia. After inhibition of HIF-1α by siRNA, we observed that hypoxia-induced ROS accumulation and apoptosis of CMECs were decreased. Meanwhile, knockdown of HIF-1α also inhibited hypoxia induced FoxO3a expression in CMECs, but did not affect FoxO1 expression. Furthermore, hypoxia-induced ROS formation and apoptosis in CMECs were correlated with the disturbance of Bcl-2 family proteins, which were abolished by FoxO3a silencing with siRNA. In conclusion, our data provide evidence that FoxO3a leads to ROS accumulation in CMECs, and in parallel, induces the disturbance of Bcl-2 family proteins which results in apoptosis.

A systems biology view of blood vessel growth and remodelling

Blood travels throughout the body in an extensive network of vessels – arteries, veins and capillaries. This vascular network is not static, but instead dynamically remodels in response to stimuli from cells in the nearby tissue. In particular, the smallest vessels – arterioles, venules and capillaries – can be extended, expanded or pruned, in response to exercise, ischaemic events, pharmacological interventions, or other physiological and pathophysiological events. In this review, we describe the multi-step morphogenic process of angiogenesis – the sprouting of new blood vessels – and the stability of vascular networks in vivo. In particular, we review the known interactions between endothelial cells and the various blood cells and plasma components they convey. We describe progress that has been made in applying computational modelling, quantitative biology and high-throughput experimentation to the angiogenesis process.

Regulation of IL-1β-induced NF-κB by hydroxylases links key hypoxic and inflammatory signaling pathways

Hypoxia is a prominent feature of chronically inflamed tissues. Oxygen-sensing hydroxylases control transcriptional adaptation to hypoxia through the regulation of hypoxia-inducible factor (HIF) and nuclear factor κB (NF-κB), both of which can regulate the inflammatory response. Furthermore, pharmacologic hydroxylase inhibitors reduce inflammation in multiple animal models. However, the underlying mechanism(s) linking hydroxylase activity to inflammatory signaling remains unclear. IL-1β, a major proinflammatory cytokine that regulates NF-κB, is associated with multiple inflammatory pathologies. We demonstrate that a combination of prolyl hydroxylase 1 and factor inhibiting HIF hydroxylase isoforms regulates IL-1β-induced NF-κB at the level of (or downstream of) the tumor necrosis factor receptor-associated factor 6 complex. Multiple proteins of the distal IL-1β-signaling pathway are subject to hydroxylation and form complexes with either prolyl hydroxylase 1 or factor inhibiting HIF. Thus, we hypothesize that hydroxylases regulate IL-1β signaling and subsequent inflammatory gene expression. Furthermore, hydroxylase inhibition represents a unique approach to the inhibition of IL-1β-dependent inflammatory signaling.

Hypoxia-inducible factor (HIF) network: insights from mathematical models

Oxygen is a crucial molecule for cellular function. When oxygen demand exceeds supply, the oxygen sensing pathway centred on the hypoxia inducible factor (HIF) is switched on and promotes adaptation to hypoxia by up-regulating genes involved in angiogenesis, erythropoiesis and glycolysis. The regulation of HIF is tightly modulated through intricate regulatory mechanisms. Notably, its protein stability is controlled by the oxygen sensing prolyl hydroxylase domain (PHD) enzymes and its transcriptional activity is controlled by the asparaginyl hydroxylase FIH (factor inhibiting HIF-1).To probe the complexity of hypoxia-induced HIF signalling, efforts in mathematical modelling of the pathway have been underway for around a decade. In this paper, we review the existing mathematical models developed to describe and explain specific behaviours of the HIF pathway and how they have contributed new insights into our understanding of the network. Topics for modelling included the switch-like response to decreased oxygen gradient, the role of micro environmental factors, the regulation by FIH and the temporal dynamics of the HIF response. We will also discuss the technical aspects, extent and limitations of these models. Recently, HIF pathway has been implicated in other disease contexts such as hypoxic inflammation and cancer through crosstalking with pathways like NFκB and mTOR. We will examine how future mathematical modelling and simulation of interlinked networks can aid in understanding HIF behaviour in complex pathophysiological situations. Ultimately this would allow the identification of new pharmacological targets in different disease settings.