A computational model of intracellular oxygen sensing by hypoxia-inducible factor HIF1 alpha

Engineered Feedback Employing Natural Hypoxia-Responsive Factors Enhances Synthetic Hypoxia Biosensors

Inadequate oxygen supply is a feature of multiple acute and chronic diseases, and hypoxia biosensors can be deployed in engineered cells to study or treat disease. Although mediators of hypoxia-responsiveness have been characterized, dynamics of this response are less understood, and there is no general approach for tuning biosensor performance to meet application-specific needs. To address these gaps, we investigated the use of genetic circuits to enhance biosensor performance through feedback, ultimately achieving both low background and amplified hypoxia-induced gene expression. To build insight into the mechanisms by which our circuits modulate performance, we developed an explanatory mathematical model. Our analysis suggests a previously unreported dual regulatory mechanism in the natural hypoxia response, providing new insights into regulatory dynamics in chronic hypoxia. This study exemplifies the potential of using synthetic gene circuits to perturb natural systems in a manner that uniquely enables the elucidation of novel facets of natural regulation.

Systems Biology of the Cancer Cell

Questions in cancer have engaged systems biologists for decades. During that time, the quantity of molecular data has exploded, but the need for abstractions, formal models, and simplifying insights has remained the same. This review brings together classic breakthroughs and recent findings in the field of cancer systems biology, focusing on cancer cell pathways for tumorigenesis and therapeutic response. Cancer cells mutate and transduce information from their environment to alter gene expression, metabolism, and phenotypic states. Understanding the molecular architectures that make each of these steps possible is a long-term goal of cancer systems biology pursued by iterating between quantitative models and experiments. We argue that such iteration is the best path to deploying targeted therapies intelligently so that each patient receives the maximum benefit for their cancer.

A Boolean network model of hypoxia, mechanosensing and TGF-β signaling captures the role of phenotypic plasticity and mutations in tumor metastasis

The tumor microenvironment aids cancer progression by promoting several cancer hallmarks, independent of cancer-related mutations. Biophysical properties of this environment, such as the stiffness of the matrix cells adhere to and local cell density, impact proliferation, apoptosis, and the epithelial to mesenchymal transition (EMT). The latter is a rate-limiting step for invasion and metastasis, enhanced in hypoxic tumor environments but hindered by soft matrices and/or high cell densities. As these influences are often studied in isolation, the crosstalk between hypoxia, biomechanical signals, and the classic EMT driver TGF-β is not well mapped, limiting our ability to predict and anticipate cancer cell behaviors in changing tumor environments. To address this, we built a Boolean regulatory network model that integrates hypoxic signaling with a mechanosensitive model of EMT, which includes the EMT-promoting crosstalk of mitogens and biomechanical signals, cell cycle control, and apoptosis. Our model reproduces the requirement of Hif-1α for proliferation, the anti-proliferative effects of strong Hif-1α stabilization during hypoxia, hypoxic protection from anoikis, and hypoxia-driven mechanosensitive EMT. We offer experimentally testable predictions about the effect of VHL loss on cancer hallmarks, with or without secondary oncogene activation. Taken together, our model serves as a predictive framework to synthesize the signaling responses associated with tumor progression and metastasis in healthy vs. mutant cells. Our single-cell model is a key step towards more extensive regulatory network models that cover damage-response and senescence, integrating most cell-autonomous cancer hallmarks into a single model that can, in turn, control the behavior of in silico cells within a tissue model of epithelial homeostasis and carcinoma.

Progressive Deactivation of Hydroxylases Controls Hypoxia-Inducible Factor-1α-Coordinated Cellular Adaptation to Graded Hypoxia

Graded hypoxia is a common microenvironment in malignant solid tumors. As a central regulator in the hypoxic response, hypoxia-inducible factor-1 (HIF-1) can induce multiple cellular processes including glycolysis, angiogenesis, and necroptosis. How cells exploit the HIF-1 pathway to coordinate different processes to survive hypoxia remains unclear. We developed an integrated model of the HIF-1α network to elucidate the mechanism of cellular adaptation to hypoxia. By numerical simulations and bifurcation analysis, we found that HIF-1α is progressively activated with worsening hypoxia due to the sequential deactivation of the hydroxylases prolyl hydroxylase domain enzymes and factor inhibiting HIF (FIH). Bistable switches control the activation and deactivation processes. As a result, glycolysis, immunosuppression, angiogenesis, and necroptosis are orderly elicited in aggravating hypoxia. To avoid the excessive accumulation of lactic acid during glycolysis, HIF-1α induces monocarboxylate transporter and carbonic anhydrase 9 sequentially to export intracellular hydrogen ions, facilitating tumor cell survival. HIF-1α-induced miR-182 facilitates vascular endothelial growth factor production to promote angiogenesis under moderate hypoxia. The imbalance between accumulation and removal of lactic acid in severe hypoxia may result in acidosis and induce cell necroptosis. In addition, the deactivation of FIH results in the destabilization of HIF-1α in anoxia. Collectively, HIF-1α orchestrates the adaptation of tumor cells to hypoxia by selectively inducing its targets according to the severity of hypoxia. Our work may provide clues for tumor therapy by targeting the HIF-1 pathway.

Preconditioning Strategies for Improving the Outcome of Fat Grafting

Autologous fat grafting is a common procedure in plastic, reconstructive, and aesthetic surgery. However, it is frequently associated with an unpredictable resorption rate of the graft depending on the engraftment kinetics. This, in turn, is determined by the interaction of the grafted adipose tissue with the tissue at the recipient site. Accordingly, preconditioning strategies have been developed following the principle of exposing these tissues in the pretransplantation phase to stimuli inducing endogenous protective and regenerative cellular adaptations, such as the upregulation of stress-response genes or the release of cytokines and growth factors. As summarized in the present review, these stimuli include hypoxia, dietary restriction, local mechanical stress, heat, and exposure to fractional carbon dioxide laser. Preclinical studies show that they promote cell viability, adipogenesis, and angiogenesis, while reducing inflammation, fibrosis, and cyst formation, resulting in a higher survival rate and quality of fat grafts in different experimental settings. Hence, preconditioning represents a promising approach to improve the outcome of fat grafting in future clinical practice. For this purpose, it is necessary to establish standardized preconditioning protocols for specific clinical applications that are efficient, safe, and easy to implement into routine procedures.

Deciphering oxygen distribution and hypoxia profiles in the tumor microenvironment: a data-driven mechanistic modeling approach

Objective. The distribution of hypoxia within tissues plays a critical role in tumor diagnosis and prognosis. Recognizing the significance of tumor oxygenation and hypoxia gradients, we introduce mathematical frameworks grounded in mechanistic modeling approaches for their quantitative assessment within a tumor microenvironment. By utilizing known blood vasculature, we aim to predict hypoxia levels across different tumor types.Approach. Our approach offers a computational method to measure and predict hypoxia using known blood vasculature. By formulating a reaction-diffusion model for oxygen distribution, we derive the corresponding hypoxia profile.Main results. The framework successfully replicates observed inter- and intra-tumor heterogeneity in experimentally obtained hypoxia profiles across various tumor types (breast, ovarian, pancreatic). Additionally, we propose a data-driven method to deduce partial differential equation models with spatially dependent parameters, which allows us to comprehend the variability of hypoxia profiles within tissues. The versatility of our framework lies in capturing diverse and dynamic behaviors of tumor oxygenation, as well as categorizing states of vascularization based on the dynamics of oxygen molecules, as identified by the model parameters.Significance. The proposed data-informed mechanistic method quantitatively assesses hypoxia in the tumor microenvironment by integrating diverse histopathological data and making predictions across different types of data. The framework provides valuable insights from both modeling and biological perspectives, advancing our comprehension of spatio-temporal dynamics of tumor oxygenation.

PICK1 modulates glycolysis and angiogenesis of hypoxic endothelial cells by regulating iron homeostasis

Iron accumulation, which is controlled by transferrin receptor 1 (TfR1), modulates hypoxia-inducible factor-1α (HIF-1α) activation and angiogenesis of hypoxic endothelial cells. The study examined the role of protein interacting with C-kinase 1 (PICK1), a scaffold protein containing PDZ domain, in regulating glycolysis and angiogenesis of hypoxic vascular endothelial cells through its potential effect on TfR1, which features a supersecondary structure that interacts with the PDZ domain. Iron chelator deferoxamine and TfR1 siRNA were employed to assess the impact of iron accumulation on angiogenesis, while the effects of PICK1 siRNA and overexpressing lentivirus on TfR1-mediated iron accumulation were also investigated in hypoxic human umbilical vein vascular endothelial cells (HUVECs). The study found that 72-h hypoxia impaired the proliferation, migration, and tube formation of HUVECs, and reduced the upregulation of vascular endothelial growth factor, HIF-1α, 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase 3, and PICK1, while increasing the expression of TfR1 as compared to 24-h hypoxia. Administration of deferoxamine or TfR1 siRNA reversed these effects and led to increased glycolysis, ATP content, and phosphofructokinase activity, along with increased PICK1 expression. PICK1 overexpression improved glycolysis, enhanced angiogenic capacity, and attenuated TfR1 protein upregulation in hypoxic HUVECs, with higher expression of angiogenic markers, which could be significantly reversed by the PDZ domain inhibitor. PICK1 knockdown exerted opposite effects. The study concluded that PICK1 modulated intracellular iron homeostasis, thereby promoting glycolysis and angiogenesis of HUVECs in response to prolonged hypoxia, at least in part, by regulating TfR1 expression.

Endothelial cells signaling and patterning under hypoxia: a mechanistic integrative computational model including the Notch-Dll4 pathway

Introduction: Several signaling pathways are activated during hypoxia to promote angiogenesis, leading to endothelial cell patterning, interaction, and downstream signaling. Understanding the mechanistic signaling differences between endothelial cells under normoxia and hypoxia and their response to different stimuli can guide therapies to modulate angiogenesis. We present a novel mechanistic model of interacting endothelial cells, including the main pathways involved in angiogenesis. Methods: We calibrate and fit the model parameters based on well-established modeling techniques that include structural and practical parameter identifiability, uncertainty quantification, and global sensitivity. Results: Our results indicate that the main pathways involved in patterning tip and stalk endothelial cells under hypoxia differ, and the time under hypoxia interferes with how different stimuli affect patterning. Additionally, our simulations indicate that Notch signaling might regulate vascular permeability and establish different Nitric Oxide release patterns for tip/stalk cells. Following simulations with various stimuli, our model suggests that factors such as time under hypoxia and oxygen availability must be considered for EC pattern control. Discussion: This project provides insights into the signaling and patterning of endothelial cells under various oxygen levels and stimulation by VEGFA and is our first integrative approach toward achieving EC control as a method for improving angiogenesis. Overall, our model provides a computational framework that can be built on to test angiogenesis-related therapies by modulation of different pathways, such as the Notch pathway.

Angiogenic Microvascular Wall Shear Stress Patterns Revealed Through Three-dimensional Red Blood Cell Resolved Modeling

The wall shear stress (WSS) exerted by blood flowing through microvascular capillaries is an established driver of new blood vessel growth, or angiogenesis. Such adaptations are central to many physiological processes in both health and disease, yet three-dimensional (3D) WSS characteristics in real angiogenic microvascular networks are largely unknown. This marks a major knowledge gap because angiogenesis, naturally, is a 3D process. To advance current understanding, we model 3D red blood cells (RBCs) flowing through rat angiogenic microvascular networks using state-of-the-art simulation. The high-resolution fluid dynamics reveal 3D WSS patterns occurring at sub-endothelial cell (EC) scales that derive from distinct angiogenic morphologies, including microvascular loops and vessel tortuosity. We identify the existence of WSS hot and cold spots caused by angiogenic surface shapes and RBCs, and notably enhancement of low WSS regions by RBCs. Spatiotemporal characteristics further reveal how fluctuations follow timescales of RBC "footprints." Altogether, this work provides a new conceptual framework for understanding how shear stress might regulate EC dynamics in vivo.

Endothelial cells signaling and patterning under hypoxia: a mechanistic integrative computational model including the Notch-Dll4 pathway

Several signaling pathways are activated during hypoxia to promote angiogenesis, leading to endothelial cell patterning, interaction, and downstream signaling. Understanding the mechanistic signaling differences between normoxia and hypoxia can guide therapies to modulate angiogenesis. We present a novel mechanistic model of interacting endothelial cells, including the main pathways involved in angiogenesis. We calibrate and fit the model parameters based on well-established modeling techniques. Our results indicate that the main pathways involved in the patterning of tip and stalk endothelial cells under hypoxia differ, and the time under hypoxia affects how a reaction affects patterning. Interestingly, the interaction of receptors with Neuropilin1 is also relevant for cell patterning. Our simulations under different oxygen concentrations indicate time- and oxygen-availability-dependent responses for the two cells. Following simulations with various stimuli, our model suggests that factors such as period under hypoxia and oxygen availability must be considered for pattern control. This project provides insights into the signaling and patterning of endothelial cells under hypoxia, contributing to studies in the field.

Bacillus pumilus TS1 alleviates Salmonella Enteritidis-induced intestinal injury in broilers

BACKGROUND

In the current context of reduced and limited antibiotic use, several pathogens and stressors cause intestinal oxidative stress in poultry, which leads to a reduced feed intake, slow or stagnant growth and development, and even death, resulting in huge economic losses to the poultry breeding industry. Oxidative stress in animals is a non-specific injury for which no targeted drug therapy is available; however, the health of poultry can be improved by adding appropriate feed additives. Bacillus pumilus, as a feed additive, promotes growth and development and reduces intestinal oxidative stress damage in poultry. Heat shock protein 70 (HSP70) senses oxidative damage and repairs unfolded and misfolded proteins; its protective effect has been widely investigated. Mitogen-activated protein kinase/protein kinase C (MAPK/PKC) and hypoxia inducible factor-1 alpha (HIF-1α) are also common proteins associated with inflammatory response induced by several stressors, but there is limited research on these proteins in the context of poultry intestinal Salmonella Enteritidis (SE) infections. In the present study, we isolated a novel strain of Bacillus pumilus with excellent performance from the feces of healthy yaks, named TS1. To investigate the effect of TS1 on SE-induced enteritis in broilers, 120 6-day-old white-feathered broilers were randomly divided into four groups (con, TS1, SE, TS1 + SE). TS1 and TS1 + SE group chickens were fed with 1.4 × 107 colony-forming units per mL of TS1 for 15 days and intraperitoneally injected with SE to establish the oxidative stress model. Then, we investigated whether TS1 protects the intestine of SE-treated broiler chickens using inflammatory cytokine gene expression analysis, stress protein quantification, antioxidant quantification, and histopathological analysis.

RESULTS

The TS1 + SE group showed lower MDA and higher GSH-Px, SOD, and T-AOC than the SE group. TS1 alleviated the effects of SE on intestinal villus length and crypt depth. Our results suggest that SE exposure increased the expression of inflammatory factors (IL-1β, IL-6, TNF-α, IL-4, and MCP-1), p38 MAPK, and PKCβ and decreased the expression of HSP60, HSP70, and HIF-1α, whereas TS1 alleviated these effects.

CONCLUSIONS

Bacillus pumilus TS1 alleviated oxidative stress damage caused by SE and attenuated the inflammatory response in broilers through MAPK/PKC regulation of HSPs/HIF-1α.

NF-κB mediated regulation of tumor cell proliferation in hypoxic microenvironment

Hypoxia is caused by a cancer-promoting milieu characterized by persistent inflammation. NF-κB and HIF-1α are critical participants in this transition. Tumor development and maintenance are aided by NF-κB, while cellular proliferation and adaptability to angiogenic signals are aided by HIF-1α. Prolyl hydroxylase-2 (PHD-2) has been hypothesized to be the key oxygen-dependent regulator of HIF-1α and NF-transcriptional B's activity. Without low oxygen levels, HIF-1α is degraded by the proteasome in a process dependent on oxygen and 2-oxoglutarate. As opposed to the normal NF-κB activation route, where NF-κB is deactivated by PHD-2-mediated hydroxylation of IKK, this method actually activates NF-κB. HIF-1α is protected from degradation by proteasomes in hypoxic cells, where it then activates transcription factors involved in cellular metastasis and angiogenesis. The Pasteur phenomenon causes lactate to build up inside the hypoxic cells. As part of a process known as lactate shuttle, MCT-1 and MCT-4 cells help deliver lactate from the blood to neighboring, non-hypoxic tumour cells. Non-hypoxic tumour cells use lactate, which is converted to pyruvate, as fuel for oxidative phosphorylation. OXOPHOS cancer cells are characterized by a metabolic switch from glucose-facilitated oxidative phosphorylation to lactate-facilitated oxidative phosphorylation. Although PHD-2 was found in OXOPHOS cells. There is no clear explanation for the presence of NF-kappa B activity. The accumulation of the competitive inhibitor of 2-oxo-glutarate, pyruvate, in non-hypoxic tumour cells is well established. So, we conclude that PHD-2 is inactive in non-hypoxic tumour cells due to pyruvate-mediated competitive suppression of 2-oxo-glutarate. This results in canonical activation of NF-κB. In non-hypoxic tumour cells, 2-oxoglutarate serves as a limiting factor, rendering PHD-2 inactive. However, FIH prevents HIF-1α from engaging in its transcriptional actions. Using the existing scientific literature, we conclude in this study that NF-κB is the major regulator of tumour cell growth and proliferation via pyruvate-mediated competitive inhibition of PHD-2.

Multimodal Potentials of Gold Nanoparticles for Bone Tissue Engineering and Regenerative Medicine: Avenues and Prospects

Osseous tissue repair has advanced due to the introduction of tissue engineering. The key elements required while engineering new tissues involve scaffolds, cells, and bioactive cues. The macrostructural to the nanostructural framework of such complex tissue has engrossed the intervention of nanotechnology for efficient neo-bone formation. Gold nanoparticles (GNPs) have recently gained interest in bone tissue regeneration due to their multimodal functionality. They are proven to modulate the properties of scaffolds and the osteogenic cells significantly. GNPs also influence different metabolic functions within the body, which directly or indirectly govern the mechanism of bone regeneration. Therefore, this review highlights nanoparticle-mediated osteogenic development, focusing on different aspects of GNPs ranging from scaffold modulation to cellular stimulation. The toxic aspects of gold nanoparticles studied so far are critically explicated, while further insight into the advancements and prospects of these nanoparticles in bone regeneration is also highlighted.

Systems biology of angiogenesis signaling: Computational models and omics

Angiogenesis is a highly regulated multiscale process that involves a plethora of cells, their cellular signal transduction, activation, proliferation, differentiation, as well as their intercellular communication. The coordinated execution and integration of such complex signaling programs is critical for physiological angiogenesis to take place in normal growth, development, exercise, and wound healing, while its dysregulation is critically linked to many major human diseases such as cancer, cardiovascular diseases, and ocular disorders; it is also crucial in regenerative medicine. Although huge efforts have been devoted to drug development for these diseases by investigation of angiogenesis‐targeted therapies, only a few therapeutics and targets have proved effective in humans due to the innate multiscale complexity and nonlinearity in the process of angiogenic signaling. As a promising approach that can help better address this challenge, systems biology modeling allows the integration of knowledge across studies and scales and provides a powerful means to mechanistically elucidate and connect the individual molecular and cellular signaling components that function in concert to regulate angiogenesis. In this review, we summarize and discuss how systems biology modeling studies, at the pathway‐, cell‐, tissue‐, and whole body‐levels, have advanced our understanding of signaling in angiogenesis and thereby delivered new translational insights for human diseases.

This article is categorized under:

Cardiovascular Diseases > Computational Models

Cancer > Computational Models

Cycling hypoxia selects for constitutive HIF stabilization

Tumors experience temporal and spatial fluctuations in oxygenation. Hypoxia inducible transcription factors (HIF-α) respond to low levels of oxygen and induce re-supply oxygen. HIF-α stabilization is typically facultative, induced by hypoxia and reduced by normoxia. In some cancers, HIF-α stabilization becomes constitutive under normoxia. We develop a mathematical model that predicts how fluctuating oxygenation affects HIF-α stabilization and impacts net cell proliferation by balancing the base growth rate, the proliferative cost of HIF-α expression, and the mortality from not expressing HIF-α during hypoxia. We compare optimal net cell proliferation rate between facultative and constitutive HIF-α regulation in environments with different oxygen profiles. We find that that facultative HIF-α regulation promotes greater net cell proliferation than constitutive regulation with stochastic or slow periodicity in oxygenation. However, cell fitness is nearly identical for both HIF-α regulation strategies under rapid periodic oxygenation fluctuations. The model thus indicates that cells constitutively expressing HIF-α may be at a selective advantage when the cost of expression is low. In cancer, this condition is known as pseudohypoxia or the “Warburg Effect”. We conclude that rapid and regular cycling of oxygenation levels selects for pseudohypoxia, and that this is consistent with the ecological theory of optimal defense.

HACE1 blocks HIF1α accumulation under hypoxia in a RAC1 dependent manner

Uncovering the mechanisms that underpin how tumor cells adapt to microenvironmental stress is essential to better understand cancer progression. The HACE1 (HECT domain and ankyrin repeat-containing E3 ubiquitin-protein ligase) gene is a tumor suppressor that inhibits the growth, invasive capacity, and metastasis of cancer cells. However, the direct regulatory pathways whereby HACE1 confers this tumor-suppressive effect remain to be fully elucidated. In this report, we establish a link between HACE1 and the major stress factor, hypoxia-inducible factor 1 alpha (HIF1α). We find that HACE1 blocks the accumulation of HIF1α during cellular hypoxia through decreased protein stability. This property is dependent on HACE1 E3 ligase activity and loss of Ras-related C3 botulinum toxin substrate 1 (RAC1), an established target of HACE1 mediated ubiquitinylation and degradation. In vivo, genetic deletion of Rac1 reversed the increased HIF1α expression observed in Hace1^-/- mice in murine KRas^G12D-driven lung tumors. An inverse relationship was observed between HACE1 and HIF1α levels in tumors compared to patient-matched normal kidney tissues, highlighting the potential pathophysiological significance of our findings. Together, our data uncover a previously unrecognized function for the HACE1 tumor suppressor in blocking HIF1α accumulation under hypoxia in a RAC1-dependent manner.

Orchestration of lincRNA-p21 and miR-155 in Modulating the Adaptive Dynamics of HIF-1α

Hypoxia-inducible factor-1 (HIF-1) is the key regulator of cellular adaptive response to hypoxia. Accumulating evidence shows that HIF-1 induces some non-coding RNAs (ncRNAs) including lncRNAs and miRNAs to modulate its own activity, enclosing several feedback loops. How the two classes of ncRNAs are orchestrated in the HIF-1-dependent adaptive response to hypoxia is poorly understood. By selecting lincRNA-p21 and miR-155 as the representatives, we develop an integrated model of the HIF-1 network comprising interlinked positive and negative feedback loops to clarify the interplay between the two ncRNAs in the hypoxic response. By numerical simulations, we find that coordination of lincRNA-p21 and miR-155 shapes the adaptive dynamics of HIF-1α: lincRNA-p21 induction in the early phase stimulates the upregulation of HIF-1α via stabilizing it, while miR-155 induction in the late phase promotes the recovery of HIF-1α via enhancing the degradation of its mRNA. Moreover, HIF-1α-induced PHD2 plays an auxiliary role in the decline of HIF-1α. In addition, lincRNA-p21 and miR-155 modulate each other via regulating HIF-1α activity. Together, lincRNA-p21 and miR-155 coordinate in modulating HIF-1α dynamics, and our work may shed light on the role for ncRNAs in the cellular adaptation to hypoxia.

Computational modeling to determine key regulators of hypoxia effects on the lactate production in the glycolysis pathway

In solid tumors, hypoxia can trigger aberrant expression of transcription factors and genes, resulting in abnormal biological functions such as altered energetic pathways in cancer cells. Glucose metabolism is an important part of this phenomenon, which is associated with changes in the functional expression of transporters and enzymes involved in the glycolysis pathway. The latter phenomenon can finally lead to the lactate accumulation and pH dysregulation in the tumor microenvironment and subsequently further invasion and metastasis of cancer cells. Having capitalized on the computational modeling, in this study, for the first time, we aimed to investigate the effects of hypoxia-induced factor-1 (HIF-1) mediated hypoxia on the magnitude of functional expression of all the enzymes and transporters involved in the glycolysis process. The main objective was to establish a quantitative relationship between the hypoxia intensity and the intracellular lactate levels and determine the key regulators of the glycolysis pathway. This model clearly showed an increase in the lactate concentration during the oxygen depletion. The proposed model also predicted that the phosphofructokinase-1 and phosphoglucomutase enzymes might play the most important roles in the regulation of the lactate production.

Defective coronary vessel organization and reduction of VEGF-A in mouse embryonic hearts with gestational mild hypoxia

BACKGROUND

Vasculature is formed by responding to homeostatic tissue demands including in developing hearts. Hypoxia generally stimulates vascular formation in which vascular endothelial growth factor A (VEGF-A) plays a critical role. Gestational hypoxia increases the risk of low intrauterine growth and low birth weight, both of which are known to increase the risk of the fetus developing cardiovascular defects. In fact, continuous gestational mild hypoxia (14% O₂ ) from the mid-embryonic stage causes cardiac anomalies accompanied by a thinning compact layer in mice in vivo. Because coronary vasculature formation is necessary for compact layers to thicken, we hypothesized that defective coronary vessel organization is related to the thinning compact layer under gestational hypoxia conditions.

RESULTS

Continuous gestational mild hypoxia (14% O₂ ) applied from embryonic day 10.5 (E10.5) reduced the expression of VEGF-A mRNA and proteins by over 60% in E12.5 hearts relative to control normoxic hearts. Formation of CD31-positive vascular plexus, blood islands, and microvessels in embryonic ventricles were stunted by gestational hypoxia compared to control E12.5 hearts.

CONCLUSIONS

Our results suggest that mild hypoxia (14% O₂ ) does not induce coronary vessel organization or VEGF-A expression in developing mouse hearts, opposing the general effects of hypoxia-triggering vascular organization and VEGF-A expression.

Interleukin-15 Signaling in HIF-1α Regulation in Natural Killer Cells, Insights Through Mathematical Models

Natural killer (NK) cells belong to the first line of host defense against infection and cancer. Cytokines, including interleukin-15 (IL-15), critically regulate NK cell activity, resulting in recognition and direct killing of transformed and infected target cells. NK cells have to adapt and respond in inflamed and often hypoxic areas. Cellular stabilization and accumulation of the transcription factor hypoxia-inducible factor-1α (HIF-1α) is a key mechanism of the cellular hypoxia response. At the same time, HIF-1α plays a critical role in both innate and adaptive immunity. While the HIF-1α hydroxylation and degradation pathway has been recently described with the help of mathematical methods, less is known concerning the mechanistic mathematical description of processes regulating the levels of HIF-1α mRNA and protein. In this work we combine mathematical modeling with experimental laboratory analysis and examine the dynamic relationship between HIF-1α mRNA, HIF-1α protein, and IL-15-mediated upstream signaling events in NK cells from human blood. We propose a system of non-linear ordinary differential equations with positive and negative feedback loops for describing the complex interplay of HIF-1α regulators. The experimental design is optimized with the help of mathematical methods, and numerical optimization techniques yield reliable parameter estimates. The mathematical model allows for the investigation and prediction of HIF-1α stabilization under different inflammatory conditions and provides a better understanding of mechanisms mediating cellular enrichment of HIF-1α. Thanks to the combination of in vitro experimental data and in silico predictions we identified the mammalian target of rapamycin (mTOR), the nuclear factor-κB (NF-κB), and the signal transducer and activator of transcription 3 (STAT3) as central regulators of HIF-1α accumulation. We hypothesize that the regulatory pathway proposed here for NK cells can be extended to other types of immune cells. Understanding the molecular mechanisms involved in the dynamic regulation of the HIF-1α pathway in immune cells is of central importance to the immune cell function and could be a promising strategy in the design of treatments for human inflammatory diseases and cancer.

Modeling the regulation of p53 activation by HIF-1 upon hypoxia

As a famous tumor suppressor, p53 is also activated under hypoxic conditions. Hypoxia-inducuble factor 1, HIF-1, is involved in the activation of p53 upon hypoxia. However, how p53 is modulated by the HIF-1 pathway to decide cell fate is less understood. In this work, we developed a network model including p53 and HIF-1 pathways to clarify the mechanism of cell fate decision in response to hypoxia. We found that HIF-1α and p53 are activated under different conditions. Under moderate hypoxia, HIF-1α is activated to induce glycolysis or angiogenesis, and promotes partial accumulation of p53 by inducing PNUTS. Under severe hypoxia, p53 rises to high levels due to ATR-dependent stabilization and promotes Mdm2-dependent HIF-1α degradation. As a result, fully activated p53 triggers apoptosis. Of note, competition for p300 between HIF-1α and p53 plays a key role in regulating their transcriptional activities. This work may advance the understanding of the mechanism for p53 regulation by HIF-1 in the hypoxic response.

Mechanism Sharing Between Genetic and Gestational Hypoxia-Induced Cardiac Anomalies

Background: Cardiac development is a dynamic process both temporally and spatially. These complex processes are often disturbed and lead to congenital cardiac anomalies that affect approximately 1% of live births. Disease-causing variants in several genetic loci lead to cardiac anomalies, with variants in transcription factor NKX2-5 gene being one of the largest variants known. Gestational hypoxia, such as seen in high-altitude pregnancy, has been known to affect cardiac development, yet the incidence and underlying mechanisms are largely unknown.

Methods and Results: Normal wild-type female mice mated with heterozygous Nkx2-5 mutant males were housed under moderate hypoxia (14% O₂) or normoxia (20.9% O₂) conditions from 10.5 days of gestation. Wild-type mice exposed to hypoxia demonstrate excessive trabeculation, ventricular septal defects, irregular morphology of interventricular septum as well as atrial septal abnormalities, which overlap with those seen in heterozygous Nkx2-5 mutant mice. Genome-wide transcriptome done by RNA-seq of a 2-day hypoxic exposure on wild-type embryos revealed abnormal transcriptomes, in which approximately 60% share those from Nkx2-5 mutants without hypoxia. Gestational hypoxia reduced the expression of Nkx2-5 proteins in more than one-half along with a reduction in phosphorylation, suggesting that abnormal Nkx2-5 function is a common mechanism shared between genetic and gestational hypoxia-induced cardiac anomalies, at least at a specific developing stage.

Conclusion: The results of our study provide insights into a common molecular mechanism underlying non-genetic and genetic cardiac anomalies.

Mathematical Model for Tissue-Level Hypoxic Response in Microfluidic Environment

Availability of essential species like oxygen is critical in shaping the dynamics of tumor growth. When the intracellular oxygen level falls below normal, it initiates major cascades in cellular dynamics leading to tumor cell survival. In a cellular block with cells growing away from the blood vessel, the scenario can be aggravated for the cells further inside the block. In this study, the dynamics of intracellular species inside a colony of tumor cells are investigated by varying the cell-block thickness and cell types in a microfluidic cell culture device. The oxygen transport across the cell block is modeled through diffusion, while ascorbate (AS) transport from the extracellular medium is addressed by a concentration-dependent uptake model. The extracellular and intracellular descriptions were coupled through the consumption and traffic of species from the microchannel to the cell block. Our model shows that the onset of hypoxia is possible in HeLa cell within minutes depending on the cell location, although the nutrient supply inside the channel is maintained in normoxic levels. This eventually leads to total oxygen deprivation inside the cell block in the extreme case, representing the development of a necrotic core that maintains a dynamic balance with growing cells and scarce supply. The numerical model reveals that species concentration and hypoxic response are different for HeLa and HelaS3 cells. Results also indicate that the long-term hypoxic response from a microfluidic cellular block stays within 5% of the values of a tissue with the basal layer. The hybrid model can be very useful in designing microfluidic experiments to satisfactorily predict the tissue-level response in cancer research.

Mathematical modelling of interacting mechanisms for hypoxia mediated cell cycle commitment for mesenchymal stromal cells

Background

Existing experimental data have shown hypoxia to be an important factor affecting the proliferation of mesenchymal stromal cells (MSCs), but the contrasting observations made at various hypoxic levels raise the questions of whether hypoxia accelerates proliferation, and how. On the other hand, in order to meet the increasing demand of MSCs, an optimised bioreactor control strategy is needed to enhance in vitro production.

Results

A comprehensive, single-cell mathematical model has been constructed in this work, which combines cellular oxygen sensing with hypoxia-mediated cell cycle progression to predict cell cycle commitment as a proxy to proliferation rate. With oxygen levels defined for in vitro cell culture, the model predicts enhanced proliferation under intermediate (2–8%) and mild (8–15%) hypoxia and cell quiescence under severe (< 2%) hypoxia. Global sensitivity analysis and quasi-Monte Carlo simulation revealed that within a certain range (+/− 100%), model parameters affect (with varying significance) the minimum commitment time, but the existence of a range of optimal oxygen tension could be preserved with the hypothesized effects of Hif2α and reactive oxygen species (ROS). It appears that Hif2α counteracts Hif1α and ROS-mediated protein deactivation under intermediate hypoxia and normoxia (20%), respectively, to regulate the response of cell cycle commitment to oxygen tension.

Conclusion

Overall, this modelling study offered an integrative framework to capture several interacting mechanisms and allowed in silico analysis of their individual and collective roles in shaping the hypoxia-mediated commitment to cell cycle. The model offers a starting point to the establishment of a suitable mechanism that can satisfactorily explain the different existing experimental observations from different studies, and warrants future extension and dedicated experimental validation to eventually support bioreactor optimisation.

Electronic supplementary material

The online version of this article (10.1186/s12918-018-0560-3) contains supplementary material, which is available to authorized users.

Computational Systems Biology in Cancer: Modeling Methods and Applications

In recent years it has become clear that carcinogenesis is a complex process, both at the molecular and cellular levels. Understanding the origins, growth and spread of cancer, therefore requires an integrated or system-wide approach. Computational systems biology is an emerging sub-discipline in systems biology that utilizes the wealth of data from genomic, proteomic and metabolomic studies to build computer simulations of intra and intercellular processes. Several useful descriptive and predictive models of the origin, growth and spread of cancers have been developed in an effort to better understand the disease and potential therapeutic approaches. In this review we describe and assess the practical and theoretical underpinnings of commonly-used modeling approaches, including ordinary and partial differential equations, petri nets, cellular automata, agent based models and hybrid systems. A number of computer-based formalisms have been implemented to improve the accessibility of the various approaches to researchers whose primary interest lies outside of model development. We discuss several of these and describe how they have led to novel insights into tumor genesis, growth, apoptosis, vascularization and therapy.

Hypoxic behavior in cells under controlled microfluidic environment

BACKGROUND

Depleted oxygen levels, known as hypoxia, causes considerable changes in the cellular metabolism. Hypoxia-inducible factors (HIF) act as the major protagonist in orchestrating manifold hypoxic responses by escaping cellular degradation mechanisms. These complex and dynamic intracellular responses are significantly dependent on the extracellular environment. In this study, we present a detailed model of a hypoxic cellular microenvironment in a microfluidic setting involving HIF hydroxylation.

METHODS

We have modeled the induction of hypoxia in a microfluidic chip by an unsteady permeation of oxygen from the microchannel through a porous polydimethylsiloxane channel wall. Extracellular and intracellular interactions were modeled with two different mathematical descriptions. Intracellular space is directly coupled to the extracellular environment through uptake and consumption of oxygen and ascorbate similar to cells in vivo.

RESULTS

Our results indicate a sharp switch in HIF hydroxylation behavior with changing prolyl hydroxylase levels from 0.1 to 4.0μM. Furthermore, we studied the effects of extracellular ascorbate concentration, using a new model, to predict its accumulation inside the cell over a relevant physiological range. In different hypoxic conditions, the cellular environment showed a significant dependence on oxygen levels in resulting intracellular response.

CONCLUSIONS

Change in hydroxylation behavior and nutrient supplementation can have significant potential in designing novel therapeutic interventions in cancer and ischemia/reperfusion injuries.

GENERAL SIGNIFICANCE

The hybrid mathematical model can effectively predict intracellular behavior due to external influences providing valuable directions in designing future experiments.

Transcriptional and Post-Transcriptional Regulation of Thrombospondin-1 Expression: A Computational Model

Hypoxia is an important physiological stress signal that drives angiogenesis, the formation of new blood vessels. Besides an increase in the production of pro-angiogenic signals such as vascular endothelial growth factor (VEGF), hypoxia also stimulates the production of anti-angiogenic signals. Thrombospondin-1 (TSP-1) is one of the anti-angiogenic factors whose synthesis is driven by hypoxia. Cellular synthesis of TSP-1 is tightly regulated by different intermediate biomolecules including proteins that interact with hypoxia-inducible factors (HIFs), transcription factors that are activated by receptor and intracellular signaling, and microRNAs which are small non-coding RNA molecules that function in post-transcriptional modification of gene expression. Here we present a computational model that describes the mechanistic interactions between intracellular biomolecules and cooperation between signaling pathways that together make up the complex network of TSP-1 regulation both at the transcriptional and post-transcriptional level. Assisted by the model, we conduct in silico experiments to compare the efficacy of different therapeutic strategies designed to modulate TSP-1 synthesis in conditions that simulate tumor and peripheral arterial disease microenvironment. We conclude that TSP-1 production in endothelial cells depends on not only the availability of certain growth factors but also the fine-tuned signaling cascades that are initiated by hypoxia.

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.

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.

The Force at the Tip--Modelling Tension and Proliferation in Sprouting Angiogenesis

Sprouting angiogenesis, where new blood vessels grow from pre-existing ones, is a complex process where biochemical and mechanical signals regulate endothelial cell proliferation and movement. Therefore, a mathematical description of sprouting angiogenesis has to take into consideration biological signals as well as relevant physical processes, in particular the mechanical interplay between adjacent endothelial cells and the extracellular microenvironment. In this work, we introduce the first phase-field continuous model of sprouting angiogenesis capable of predicting sprout morphology as a function of the elastic properties of the tissues and the traction forces exerted by the cells. The model is very compact, only consisting of three coupled partial differential equations, and has the clear advantage of a reduced number of parameters. This model allows us to describe sprout growth as a function of the cell-cell adhesion forces and the traction force exerted by the sprout tip cell. In the absence of proliferation, we observe that the sprout either achieves a maximum length or, when the traction and adhesion are very large, it breaks. Endothelial cell proliferation alters significantly sprout morphology, and we explore how different types of endothelial cell proliferation regulation are able to determine the shape of the growing sprout. The largest region in parameter space with well formed long and straight sprouts is obtained always when the proliferation is triggered by endothelial cell strain and its rate grows with angiogenic factor concentration. We conclude that in this scenario the tip cell has the role of creating a tension in the cells that follow its lead. On those first stalk cells, this tension produces strain and/or empty spaces, inevitably triggering cell proliferation. The new cells occupy the space behind the tip, the tension decreases, and the process restarts. Our results highlight the ability of mathematical models to suggest relevant hypotheses with respect to the role of forces in sprouting, hence underlining the necessary collaboration between modelling and molecular biology techniques to improve the current state-of-the-art.

Sprouting angiogenesis—a process by which new blood vessels grow from existing ones—is an ubiquitous phenomenon in health and disease of higher organisms, playing a crucial role in organogenesis, wound healing, inflammation, as well as on the onset and progression of over 50 different diseases such as cancer, rheumatoid arthritis and diabetes. Mathematical models have the ability to suggest relevant hypotheses with respect to the mechanisms of cell movement and rearrangement within growing vessel sprouts. The inclusion of both biochemical and mechanical processes in a mathematical model of sprouting angiogenesis permits to describe sprout extension as a function of the forces exerted by the cells in the tissue. It also allows to question the regulation of biochemical processes by mechanical forces and vice-versa. In this work we present a compact model of sprouting angiogenesis that includes the mechanical characteristics of the vessel and the tissue. We use this model to suggest the mechanism for the regulation of proliferation within sprout formation. We conclude that the tip cell has the role of creating a tension in the cells that follow its lead. On those first cells of the stalk, this tension produces strain and/or empty spaces, inevitably triggering cell proliferation. The new cells occupy the space behind the tip, the tension decreases, and the process restarts. The modelling strategy used, deemed phase-field, permits to describe the evolution of the shape of different domains in complex systems. It is focused on the movement of the interfaces between the domains, and not on an exhaustive description of the transport properties within each domain. For this reason, it requires a reduced number of parameters, and has been used extensively in modelling other biological phenomena such as tumor growth. The coupling of mechanical and biochemical processes in a compact mathematical model of angiogenesis will enable the study of lumen formation and aneurisms in the near future. Also, this framework will allow the study of the action of flow in vessel remodelling, since local forces can readily be coupled with cell movement to obtain the final vessel morphology.

A Systems Pharmacology Model of Erythropoiesis in Mice Induced by Small Molecule Inhibitor of Prolyl Hydroxylase Enzymes

Mammalian erythropoiesis is a conserved process tightly controlled by the hypoxia-inducible factor (HIF1) pathway. In this study, a small molecule inhibitor (PHI-1) of prolyl-hydroxylase-2 (PHD2) enzyme involved in regulating HIF1α levels was orally administered to male BALB/c mice at 10 and 30 mg/kg. A systems pharmacology model was developed based on the measured PHI-1 plasma exposures, kidney HIF1α, kidney erythropoietin (EPO) mRNA, plasma EPO, reticulocyte counts, red blood cells, and hemoglobin levels. The model fit resulted in the estimation of drug potency (IC₅₀: 1.7μM), and systems parameters such as EPO mRNA turnover (k_{deg_EPOmRNA}: 0.43 hr^-1) and mean lifespan of reticulocytes (T_r: 81 hours). The model correctly described the observed 30–40-fold increase in kidney HIF1α protein, ∼1,000 fold increase in EPO mRNA and 2–3-fold increase in the reticulocytes at 30 mg/kg. This study presents the first parsimonious systems model of erythropoiesis to quantitatively describe the in vivo effects of PHD2 inhibition.

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.

Cancer systems biology and modeling: microscopic scale and multiscale approaches

Cancer has become known as a complex and systematic disease on macroscopic, mesoscopic and microscopic scales. Systems biology employs state-of-the-art computational theories and high-throughput experimental data to model and simulate complex biological procedures such as cancer, which involves genetic and epigenetic, in addition to intracellular and extracellular complex interaction networks. In this paper, different systems biology modeling techniques such as systems of differential equations, stochastic methods, Boolean networks, Petri nets, cellular automata methods and agent-based systems are concisely discussed. We have compared the mentioned formalisms and tried to address the span of applicability they can bear on emerging cancer modeling and simulation approaches. Different scales of cancer modeling, namely, microscopic, mesoscopic and macroscopic scales are explained followed by an illustration of angiogenesis in microscopic scale of the cancer modeling. Then, the modeling of cancer cell proliferation and survival are examined on a microscopic scale and the modeling of multiscale tumor growth is explained along with its advantages.

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.

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.

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

Activation of the hypoxia-inducible factor (HIF) pathway is a critical step in the transcriptional response to hypoxia. Although many of the key proteins involved have been characterised, the dynamics of their interactions in generating this response remain unclear. In the present study, we have generated a comprehensive mathematical model of the HIF-1α pathway based on core validated components and dynamic experimental data, and confirm the previously described connections within the predicted network topology. Our model confirms previous work demonstrating that the steps leading to optimal HIF-1α transcriptional activity require sequential inhibition of both prolyl- and asparaginyl-hydroxylases. We predict from our model (and confirm experimentally) that there is residual activity of the asparaginyl-hydroxylase FIH (factor inhibiting HIF) at low oxygen tension. Furthermore, silencing FIH under conditions where prolyl-hydroxylases are inhibited results in increased HIF-1α transcriptional activity, but paradoxically decreases HIF-1α stability. Using a core module of the HIF network and mathematical proof supported by experimental data, we propose that asparaginyl hydroxylation confers a degree of resistance upon HIF-1α to proteosomal degradation. Thus, through in vitro experimental data and in silico predictions, we provide a comprehensive model of the dynamic regulation of HIF-1α transcriptional activity by hydroxylases and use its predictive and adaptive properties to explain counter-intuitive biological observations.

Regulation of SREBP during hypoxia requires Ofd1-mediated control of both DNA binding and degradation

The prolyl hydroxylase Ofd1 enables yeast cells to adapt to hypoxia by regulating the DNA binding and degradation of the hypoxic transcription factor Sre1N. These two regulatory functions are inseparable in vivo. A mathematical model of the Ofd1 system is used to show that both functions are necessary for Ofd1 to work as observed.

Cells adapt to changes in ambient oxygen by changing their gene expression patterns. In fission yeast, the sterol regulatory element–binding protein Sre1 is proteolytically cleaved under low oxygen, and its N-terminal segment (Sre1N) serves as a hypoxic transcription factor. When oxygen is present, the prolyl hydroxylase Ofd1 down-regulates Sre1N activity in two ways: first, by inhibiting its binding to DNA, and second, by accelerating its degradation. Here we use a mathematical model to assess what each of these two regulatory functions contributes to the hypoxic response of the cell. By disabling individual regulatory functions in the model, which would be difficult in vivo, we found that the Ofd1 function of inhibiting Sre1N binding to DNA is essential for oxygen-dependent Sre1N regulation. The other Ofd1 function of accelerating Sre1N degradation is necessary for the yeast to quickly turn off its hypoxic response when oxygen is restored. In addition, the model predicts that increased Ofd1 production at low oxygen plays an important role in the hypoxic response, and the model indicates that the Ofd1 binding partner Nro1 tunes the response to oxygen. This model quantifies our understanding of a novel oxygen-sensing mechanism that is widely conserved.

Beyond the antioxidant: the double life of vitamin C

When considering the history of vitamin C, and the names given to this molecule in early days, the Latin proverb nomen est omen suddenly comes to mind. Around 1920, when Casimir Funk introduced the term Vitamin C to indicate the nutritional factor necessary to prevent the pathological state known as scurvy, the nature of the active molecule was still unknown (Davies MB, Austin J, Partridge DA (1991) Vitamin C: Its chemistry and biochemistry. The Royal Society of Chemistry, Cambridge UK). Almost in the same years, Albert Szent-Giörgyi was striving to identify a new 6-carbon sugar he had obtained in crystal form from oranges, lemons, cabbage and adrenal glands. As humorously described by Szent-Giörgyi himself (Szent-Giörgyi A (1963) Lost in the twentieth century. Annu Rev Biochem 36:1-15), he intended to name this yet unknown carbohydrate "ignose". When this name was rejected by Sir Arthur Harden, editor of the Biochemical Journal, he suggested to name it "godnose", meaning that only God could know the real identity of the molecule. Obviously, also this choice was considered inappropriate by Harden, who suggested the plain name "hexuronic acid". Only later, when the structure of "hexuronic acid" had been completely elucidated, and biological tests performed by Swirbely identified this molecule as the anti-scurvy factor vitamin C, Szent-Giörgyi and Walter Norman Haworth decided to eventually name it ascorbic acid (Szent-Giörgyi A (1963) Lost in the twentieth century. Annu Rev Biochem 36:1-15). "Ascorbic" literally means "against scurvy", and scurvy is known to be mainly due to the inactivation of some important dioxygenases involved in the synthesis of a few key molecules, including different collagen forms (De Tullio MC (2004) How does ascorbic acid prevent scurvy? A survey of the nonantioxidant functions of vitamin C. In: Asard H, May J, Smirnoff N (eds) Vitamin C, its functions and biochemistry in animals and plants. Bios Scientific Publishers, Oxford, UK, pp. 159-172). All this has very little to do with the celebrated role of ascorbic acid (ASC) as an antioxidant. So, if the fate of ASC had to be found in its name, its role in the prevention of scurvy (i.e. beyond the antioxidant function) should be considered its main feature. But, in spite of more than 80 years of extensive research (34,424 hits in a PubMed query on January 6 2007), an unprecedented popularity among the general public, an estimated market of several billion dollars (Hancock RD, Viola R (2005) Improving the nutritional value of crops through enhancement of l-ascorbic acid (vitamin C) content: Rationale and biotechnological opportunities. J Agr Food Chem 53:5248-5257), we should honestly conclude that the fate of vitamin C is still in the first name it received, many years ago: we still ignore much of its actual relevance in cell metabolism, although we are progressively getting aware of the many facets of this fascinating molecule, and its direct involvement in the regulation of apparently unrelated pathways (Arrigoni O, De Tullio MC (2002) Ascorbic acid, much more than just an antioxidant. Biochim Biophys Acta 1569:1-9; De Tullio MC, Arrigoni O (2004) Hopes, disillusions and more hopes from vitamin C. Cell Mol Life Sci 61:209-219; Duarte TL, Lunec J (2005) When is an antioxidant not an antioxidant? A review of novel actions and reactions of vitamin C. Free Rad Res 39:671-686). Recent data on ASC involvement in cell signalling and gene expression open new perspectives, that will be presented and discussed in this chapter.

Molecular mechanisms of chronic intermittent hypoxia and hypertension

Obstructive sleep apnea (OSA) is characterized by episodes of repeated airway obstruction resulting in cessation (apnea) or reduction (hypopnea) in airflow during sleep. These events lead to intermittent hypoxia and hypercapnia, sleep fragmentation, and changes in intrathoracic pressure, and are associated with a marked surge in sympathetic activity and an abrupt increase in blood pressure. Blood pressure remains elevated during wakefulness despite the absence of obstructive events resulting in a high prevalence of hypertension in patients with OSA. There is substantial evidence that suggests that chronic intermittent hypoxia (CIH) leads to sustained sympathoexcitation during the day and changes in vasculature resulting in hypertension in patients with OSA. Mechanisms of sympathoexcitation include augmentation of peripheral chemoreflex sensitivity and a direct effect on central sites of sympathetic regulation. Interestingly, the vascular changes that occur with CIH have been ascribed to the same molecules that have been implicated in the augmented sympathetic tone in CIH. This review will discuss the hypothesized molecular mechanisms involved in the development of hypertension with CIH, will build a conceptual model for the development of hypertension following CIH, and will propose a systems biology approach in further elucidating the relationship between CIH and the development of hypertension.

Oxygen Sensing, Cardiac Ischemia, HIF-1α and Some Emerging Concepts

Oxygen plays a critical role in the perpetuation and propagation of almost all forms of life. The primary site of cellular oxygen consumption is the mitochondrial electron transport chain and in addition, oxygen is also used as a substrate for various enzymes involved in cellular homeostasis. Although our knowledge of the biochemistry and physiology of oxygen transport is century old, recent development of sophisticated tools of biophysical chemistry revealed that tissue oxygenation and oxygen sensing is a highly evolved process, especially in mammals. Perturbation of normal oxygen supply is associated with diseases like tumorigenesis, myocardial infarction and stroke. Available information suggests that when tissue oxygen supply is limited, mitochondria emanate signals involving reactive oxygen species generation which in turn stabilizes oxygen sensing transcription factor HIF-1. Upon stabilization, HIF-1 elicits necessary genetic response to cope with the diminished oxygen level. In view of such critical role of HIF-1 in cellular oxygen sensing, recently there has been a heightened interest in understanding the biology of HIF-1 in the context of cardiovascular system. The following review describes some of the recent advances in this regard.

Astrocyte-neuron lactate shuttle may boost more ATP supply to the neuron under hypoxic conditions--in silico study supported by in vitro expression data

Background

Neuro-glial interactions are important for normal functioning of the brain as well as brain energy metabolism. There are two major working models - in the classical view, both neurons and astrocytes can utilize glucose as the energy source through oxidative metabolism, whereas in the astrocyte-neuron lactate shuttle hypothesis (ANLSH) it is the astrocyte which can consume glucose through anaerobic glycolysis to pyruvate and then to lactate, and this lactate is secreted to the extracellular space to be taken up by the neuron for further oxidative degradation.

Results

In this computational study, we have included hypoxia-induced genetic regulation of these enzymes and transporters, and analyzed whether the ANLSH model can provide an advantage to either cell type in terms of supplying the energy demand. We have based this module on our own experimental analysis of hypoxia-dependent regulation of transcription of key metabolic enzymes. Using this experimentation-supported in silico modeling, we show that under both normoxic and hypoxic conditions in a given time period ANLSH model does indeed provide the neuron with more ATP than in the classical view.

Conclusions

Although the ANLSH is energetically more favorable for the neuron, it is not the case for the astrocyte in the long term. Considering the fact that astrocytes are more resilient to hypoxia, we would propose that there is likely a switch between the two models, based on the energy demand of the neuron, so as to maintain the survival of the neuron under hypoxic or glucose-and-oxygen-deprived conditions.

Module-based multiscale simulation of angiogenesis in skeletal muscle

Background

Mathematical modeling of angiogenesis has been gaining momentum as a means to shed new light on the biological complexity underlying blood vessel growth. A variety of computational models have been developed, each focusing on different aspects of the angiogenesis process and occurring at different biological scales, ranging from the molecular to the tissue levels. Integration of models at different scales is a challenging and currently unsolved problem.

Results

We present an object-oriented module-based computational integration strategy to build a multiscale model of angiogenesis that links currently available models. As an example case, we use this approach to integrate modules representing microvascular blood flow, oxygen transport, vascular endothelial growth factor transport and endothelial cell behavior (sensing, migration and proliferation). Modeling methodologies in these modules include algebraic equations, partial differential equations and agent-based models with complex logical rules. We apply this integrated model to simulate exercise-induced angiogenesis in skeletal muscle. The simulation results compare capillary growth patterns between different exercise conditions for a single bout of exercise. Results demonstrate how the computational infrastructure can effectively integrate multiple modules by coordinating their connectivity and data exchange. Model parameterization offers simulation flexibility and a platform for performing sensitivity analysis.

Conclusions

This systems biology strategy can be applied to larger scale integration of computational models of angiogenesis in skeletal muscle, or other complex processes in other tissues under physiological and pathological conditions.

An insult-inducible vector system activated by hypoxia and oxidative stress for neuronal gene therapy

Gene therapy has demonstrated the protective potential of a variety of genes against stroke. However, conventional gene therapy vectors are limited due to the inability to temporally control their expression, which can sometimes lead to deleterious side effects. Thus, an inducible vector that can be temporally controlled and activated by the insult itself would be advantageous. Using hypoxia responsive elements (HRE) and antioxidant responsive elements (ARE), we have constructed an insult-inducible vector activated by hypoxia and reactive oxygen species (ROS). In COS7 cells, the inducible ARE-HRE-luciferase vectors are highly activated by oxygen deprivation, hydrogen peroxide treatment, and the ROS-induced transcription factor NF-E2-related factor 2 (Nrf2). Using a defective herpes virus, the neuroprotective potential of this inducible vector was tested by over-expressing the transcription factor Nrf2. In primary cortical cultures, expression of the inducible ARE-HRE-Nrf2 protects against oxygen glucose deprivation, similar to that afforded by the constitutively expressed Nrf2. This ARE+HRE vector system is advantageous in that it allows the expression of a transgene to be activated not only during hypoxia but also maintained after reperfusion, thus prolonging the transgene expression during an ischemic insult. This insult-inducible vector system will be a valuable gene therapy tool for activating therapeutic/protective genes in cerebrovascular diseases.

Computational models of VEGF-associated angiogenic processes in cancer

Tumour angiogenesis allows a growing mass of cancer cells to overcome oxygen diffusion limitation and to increase cell survival. The growth of capillaries from pre-existing blood vessels is the result of numerous signalling cascades involving different molecules and of cellular events involving multiple cell and tissue types. Computational models offer insight into the mechanisms governing angiogenesis and provide quantitative information on parameters difficult to assess by experiments alone. In this article, we summarize results from computational models of tumour angiogenic processes with a focus on the molecular-detailed vascular endothelial growth factor-associated models that have been developed in our laboratory, spanning multiple scales from the molecular to whole body.

Effects of polynitrogen compounds on the activity of recombinant human HIF-1α prolyl hydroxylase 3 in E. coli

Hypoxia inducible factor 1α (HIF-1α) becomes an important regulation factor within the histiocyte when it is under the hypoxia condition. Recently, prolyl hydroxylases (PHDs) have been identified to inactivation HIF-lα by hydroxylation. In this study, polynitrogen compounds were screened as HIF-1α PHD3 inhibitors. The coding region of human PHD3 DNA was optimized by using synonymous codons according to the code bias of Escherichia coli. Soluble and active human PHD3 was expressed in the E. coli with a Trx fusion tag under a lower induction temperature of 25°C. Mass spectrometry analysis of the resultant peptide product indicated a mass increase of 16 daltons, consistent with hydroxylation of the proline residue in the HIF-1α (556-574) peptide substrate. Polynitrogen compounds (1-4) inhibited the enzymatic hydroxylation of HIF-1α peptide in a concentration-dependent manner, and the apparent IC(50) values were 29.5, 16.0, 12.8 and 60.4 μM respectively. Double reciprocal (1/V versus 1/[HIF-1α peptide]) plots showed that these compounds are noncompetitive inhibitors of the hydroxylation by recombinant human PHD3 with K(i) values of 67.0, 25.3, 67.3, and 82.1 μM respectively. On the other hand, the metal complexes of these polynitrogen compounds (1-4) cannot inhibit the catalytical activity of PHD3. We hypothesized that the inhibitory mechanism of PHD3 activity by polynitrogen compounds is due to their binding to iron to form stable coordination complexes. Our results in this study indicated that polynitrogen compounds (1-4) could be potential inhibitors of PHD3 to regulate the transcriptional activity of HIF-1α.

Systems biology of pro-angiogenic therapies targeting the VEGF system

Vascular endothelial growth factor (VEGF) is a family of cytokines for which the dysregulation of expression is involved in many diseases; for some, excess VEGF causes pathological hypervascularization, while for others VEGF-induced vascular remodeling may alleviate ischemia and/or hypoxia. Anti-angiogenic therapies attacking the VEGF pathway have begun to live up to their promise for treatment of certain cancers and of age-related macular degeneration. However, the corollary is not yet true: in coronary artery disease and peripheral artery disease, clinical trials of pro-angiogenic VEGF delivery have not, so far, proven successful. The VEGF and VEGF-receptor system is complex, with at least five ligand genes, some encoding multiple protein isoforms and five receptor genes. A systems biology approach for designing pro-angiogenic therapies, using a combination of quantitative experimental approaches and detailed computational models, is essential to deal with this complexity and to understand the effects of drugs targeting the system. This approach allows us to learn from unsuccessful clinical trials and to design and test novel single therapeutics or combinations of therapeutics. Among the parameters that can be varied in order to determine optimal strategy are dosage, timing of multiple doses, route of administration, and the molecular target.