Mathematical modeling of GATA-switching for regulating the differentiation of hematopoietic stem cell

Human Stem Cells for Ophthalmology: Recent Advances in Diagnostic Image Analysis and Computational Modelling

Purpose of Review

To explore the advances and future research directions in image analysis and computational modelling of human stem cells (hSCs) for ophthalmological applications.

Recent Findings

hSCs hold great potential in ocular regenerative medicine due to their application in cell-based therapies and in disease modelling and drug discovery using state-of-the-art 2D and 3D organoid models. However, a deeper characterisation of their complex, multi-scale properties is required to optimise their translation to clinical practice. Image analysis combined with computational modelling is a powerful tool to explore mechanisms of hSC behaviour and aid clinical diagnosis and therapy.

Summary

Many computational models draw on a variety of techniques, often blending continuum and discrete approaches, and have been used to describe cell differentiation and self-organisation. Machine learning tools are having a significant impact in model development and improving image classification processes for clinical diagnosis and treatment and will be the focus of much future research.

A robust method for designing multistable systems by embedding bistable subsystems

Although multistability is an important dynamic property of a wide range of complex systems, it is still a challenge to develop mathematical models for realising high order multistability using realistic regulatory mechanisms. To address this issue, we propose a robust method to develop multistable mathematical models by embedding bistable models together. Using the GATA1-GATA2-PU.1 module in hematopoiesis as the test system, we first develop a tristable model based on two bistable models without any high cooperative coefficients, and then modify the tristable model based on experimentally determined mechanisms. The modified model successfully realises four stable steady states and accurately reflects a recent experimental observation showing four transcriptional states. In addition, we develop a stochastic model, and stochastic simulations successfully realise the experimental observations in single cells. These results suggest that the proposed method is a general approach to develop mathematical models for realising multistability and heterogeneity in complex systems.

Reprogramming multistable monotone systems with application to cell fate control

Multistability is a key property of dynamical systems modeling cellular regulatory networks implicated in cell fate decisions, where, different stable steady states usually represent distinct cell phenotypes. Monotone network motifs are highly represented in these regulatory networks. In this paper, we leverage the properties of monotone dynamical systems to provide theoretical results that guide the selection of inputs that trigger a transition, i.e., reprogram the network, to a desired stable steady state. We first show that monotone dynamical systems with bounded trajectories admit a minimum and a maximum stable steady state. Then, we provide input choices that are guaranteed to reprogram the system to these extreme steady states. For intermediate states, we provide an input space that is guaranteed to contain an input that reprograms the system to the desired state. We then provide implementation guidelines for finite-time procedures that search this space for such an input, along with rules to prune parts of the space during search. We demonstrate these results on simulations of two recurrent regulatory network motifs: self-activation within mutual antagonism and self-activation within mutual cooperation. Our results depend uniquely on the structure of the network and are independent of specific parameter values.

A non-linear reverse-engineering method for inferring genetic regulatory networks

Hematopoiesis is a highly complex developmental process that produces various types of blood cells. This process is regulated by different genetic networks that control the proliferation, differentiation, and maturation of hematopoietic stem cells (HSCs). Although substantial progress has been made for understanding hematopoiesis, the detailed regulatory mechanisms for the fate determination of HSCs are still unraveled. In this study, we propose a novel approach to infer the detailed regulatory mechanisms. This work is designed to develop a mathematical framework that is able to realize nonlinear gene expression dynamics accurately. In particular, we intended to investigate the effect of possible protein heterodimers and/or synergistic effect in genetic regulation. This approach includes the Extended Forward Search Algorithm to infer network structure (top-down approach) and a non-linear mathematical model to infer dynamical property (bottom-up approach). Based on the published experimental data, we study two regulatory networks of 11 genes for regulating the erythrocyte differentiation pathway and the neutrophil differentiation pathway. The proposed algorithm is first applied to predict the network topologies among 11 genes and 55 non-linear terms which may be for heterodimers and/or synergistic effect. Then, the unknown model parameters are estimated by fitting simulations to the expression data of two different differentiation pathways. In addition, the edge deletion test is conducted to remove possible insignificant regulations from the inferred networks. Furthermore, the robustness property of the mathematical model is employed as an additional criterion to choose better network reconstruction results. Our simulation results successfully realized experimental data for two different differentiation pathways, which suggests that the proposed approach is an effective method to infer the topological structure and dynamic property of genetic regulations.

Intracellular Energy Variability Modulates Cellular Decision-Making Capacity

Cells generate phenotypic diversity both during development and in response to stressful and changing environments, aiding survival. Functionally vital cell fate decisions from a range of phenotypic choices are made by regulatory networks, the dynamics of which rely on gene expression and hence depend on the cellular energy budget (and particularly ATP levels). However, despite pronounced cell-to-cell ATP differences observed across biological systems, the influence of energy availability on regulatory network dynamics is often overlooked as a cellular decision-making modulator, limiting our knowledge of how energy budgets affect cell behaviour. Here, we consider a mathematical model of a highly generalisable, ATP-dependent, decision-making regulatory network, and show that cell-to-cell ATP variability changes the sets of decisions a cell can make. Our model shows that increasing intracellular energy levels can increase the number of supported stable phenotypes, corresponding to increased decision-making capacity. Model cells with sub-threshold intracellular energy are limited to a singular phenotype, forcing the adoption of a specific cell fate. We suggest that energetic differences between cells may be an important consideration to help explain observed variability in cellular decision-making across biological systems.

Limit-cycle oscillatory coexpression of cross-inhibitory transcription factors: a model mechanism for lineage promiscuity

Lineage switches are genetic regulatory motifs that govern and maintain the commitment of a developing cell to a particular cell fate. A canonical example of a lineage switch is the pair of transcription factors PU.1 and GATA-1, of which the former is affiliated with the myeloid and the latter with the erythroid lineage within the haematopoietic system. On a molecular level, PU.1 and GATA-1 positively regulate themselves and antagonize each other via direct protein-protein interactions. Here we use mathematical modelling to identify a novel type of dynamic behaviour that can be supported by such a regulatory architecture. Guided by the specifics of the PU.1-GATA-1 interaction, we formulate, using the law of mass action, a system of differential equations for the key molecular concentrations. After a series of systematic approximations, the system is reduced to a simpler one, which is tractable to phase-plane and linearization methods. The reduced system formally resembles, and generalizes, a well-known model for competitive species from mathematical ecology. However, in addition to the qualitative regimes exhibited by a pair of competitive species (exclusivity, bistable exclusivity, stable-node coexpression) it also allows for oscillatory limit-cycle coexpression. A key outcome of the model is that, in the context of cell-fate choice, such oscillations could be harnessed by a differentiating cell to prime alternately for opposite outcomes; a bifurcation-theory approach is adopted to characterize this possibility.

Stochastic multistationarity in a model of the hematopoietic stem cell differentiation network

A central issue in the analysis of multi-stable systems is that of controlling the relative size of the basins of attraction of alternative states through suitable choices of system parameters. We are interested here mainly in the stochastic version of this problem, that of shaping the stationary probability distribution of a Markov chain so that various alternative modes become more likely than others. Although many of our results are more general, we were motivated by an important biological question, that of cell differentiation. In the mathematical modeling of cell differentiation, it is common to think of internal states of cells (quanfitied by activation levels of certain genes) as determining the different cell types. Specifically, we study here the "PU.1/GATA-1 circuit" which is involved in the control of the development of mature blood cells from hematopoietic stem cells (HSCs). All mature, specialized blood cells have been shown to be derived from multipotent HSCs. Our first contribution is to introduce a rigorous chemical reaction network model of the PU.1/GATA-1 circuit, which incorporates current biological knowledge. We then find that the resulting ODE model of these biomolecular reactions is incapable of exhibiting multistability, contradicting the fact that differentiation networks have, by definition, alternative stable steady states. When considering instead the stochastic version of this chemical network, we analytically construct the stationary distribution, and are able to show that this distribution is indeed capable of admitting a multiplicity of modes. Finally, we study how a judicious choice of system parameters serves to bias the probabilities towards different stationary states. We remark that certain changes in system parameters can be physically implemented by a biological feedback mechanism; tuning this feedback gives extra degrees of freedom that allow one to assign higher likelihood to some cell types over others.

Computational Analysis of Altering Cell Fate

The notion of reprogramming cell fate is a direct challenge to the traditional view in developmental biology that a cell's phenotypic identity is sealed after undergoing differentiation. Direct experimental evidence, beginning with the somatic cell nuclear transfer experiments of the twentieth century and culminating in the more recent breakthroughs in transdifferentiation and induced pluripotent stem cell (iPSC) reprogramming, have rewritten the rules for what is possible with cell fate transformation. Research is ongoing in the manipulation of cell fate for basic research in disease modeling, drug discovery, and clinical therapeutics. In many of these cell fate reprogramming experiments, there is often little known about the genetic and molecular changes accompanying the reprogramming process. However, gene regulatory networks (GRNs) can in some cases be implicated in the switching of phenotypes, providing a starting point for understanding the dynamic changes that accompany a given cell fate reprogramming process. In this chapter, we present a framework for computationally analyzing cell fate changes by mathematically modeling these GRNs. We provide a user guide with several tutorials of a set of techniques from dynamical systems theory that can be used to probe the intrinsic properties of GRNs as well as study their responses to external perturbations.

Kinetic models of hematopoietic differentiation

As cell and molecular biology is becoming increasingly quantitative, there is an upsurge of interest in mechanistic modeling at different levels of resolution. Such models mostly concern kinetics and include gene and protein interactions as well as cell population dynamics. The final goal of these models is to provide experimental predictions, which is now taking on. However, even without matured predictions, kinetic models serve the purpose of compressing a plurality of experimental results into something that can empower the data interpretation, and importantly, suggesting new experiments by turning "knobs" in silico. Once formulated, kinetic models can be executed in terms of molecular rate equations for concentrations or by stochastic simulations when only a limited number of copies are involved. Developmental processes, in particular those of stem and progenitor cell commitments, are not only topical but also particularly suitable for kinetic modeling due to the finite number of key genes involved in cellular decisions. Stem and progenitor cell commitment processes have been subject to intense experimental studies over the last decade with some emphasis on embryonic and hematopoietic stem cells. Gene and protein interactions governing these processes can be modeled by binary Boolean rules or by continuous-valued models with interactions set by binding strengths. Conceptual insights along with tested predictions have emerged from such kinetic models. Here we review kinetic modeling efforts applied to stem cell developmental systems with focus on hematopoiesis. We highlight the future challenges including multi-scale models integrating cell dynamical and transcriptional models. This article is categorized under: Models of Systems Properties and Processes > Mechanistic Models Developmental Biology > Stem Cell Biology and Regeneration.

Reverse-engineering of gene networks for regulating early blood development from single-cell measurements

BACKGROUND

Recent advances in omics technologies have raised great opportunities to study large-scale regulatory networks inside the cell. In addition, single-cell experiments have measured the gene and protein activities in a large number of cells under the same experimental conditions. However, a significant challenge in computational biology and bioinformatics is how to derive quantitative information from the single-cell observations and how to develop sophisticated mathematical models to describe the dynamic properties of regulatory networks using the derived quantitative information.

METHODS

This work designs an integrated approach to reverse-engineer gene networks for regulating early blood development based on singel-cell experimental observations. The wanderlust algorithm is initially used to develop the pseudo-trajectory for the activities of a number of genes. Since the gene expression data in the developed pseudo-trajectory show large fluctuations, we then use Gaussian process regression methods to smooth the gene express data in order to obtain pseudo-trajectories with much less fluctuations. The proposed integrated framework consists of both bioinformatics algorithms to reconstruct the regulatory network and mathematical models using differential equations to describe the dynamics of gene expression.

RESULTS

The developed approach is applied to study the network regulating early blood cell development. A graphic model is constructed for a regulatory network with forty genes and a dynamic model using differential equations is developed for a network of nine genes. Numerical results suggests that the proposed model is able to match experimental data very well. We also examine the networks with more regulatory relations and numerical results show that more regulations may exist. We test the possibility of auto-regulation but numerical simulations do not support the positive auto-regulation. In addition, robustness is used as an importantly additional criterion to select candidate networks.

CONCLUSION

The research results in this work shows that the developed approach is an efficient and effective method to reverse-engineer gene networks using single-cell experimental observations.

The transcription factor GATA1 regulates NBEAL2 expression through a long-distance enhancer

Gray platelet syndrome is named after the gray appearance of platelets due to the absence of α-granules. It is caused by recessive mutations in NBEAL2, resulting in macrothrombocytopenia and myelofibrosis. Though using the term gray platelets for GATA1 deficiency has been debated, a reduced number of α-granules has been described for macrothrombocytopenia due to GATA1 mutations. We compared platelet size and number of α-granules for two NBEAL2 and two GATA1-deficient patients and found reduced numbers of α-granules for all, with the defect being more pronounced for NBEAL2 deficiency. We further hypothesized that the granule defect for GATA1 is due to a defective control of NBEAL2 expression. Remarkably, platelets from two patients, and Gata1-deficient mice, expressed almost no NBEAL2. The differentiation of GATA1 patient-derived CD34⁺ stem cells to megakaryocytes showed defective proplatelet and α-granule formation with strongly reduced NBEAL2 protein and ribonucleic acid expression. Chromatin immunoprecipitation sequencing revealed 5 GATA binding sites in a regulatory region 31 kb upstream of NBEAL2 covered by a H3K4Me1 mark indicative of an enhancer locus. Luciferase reporter constructs containing this region confirmed its enhancer activity in K562 cells, and mutagenesis of the GATA1 binding sites resulted in significantly reduced enhancer activity. Moreover, DNA binding studies showed that GATA1 and GATA2 physically interact with this enhancer region. GATA1 depletion using small interfering ribonucleic acid in K562 cells also resulted in reduced NBEAL2 expression. In conclusion, we herein show a long-distance regulatory region with GATA1 binding sites as being a strong enhancer for NBEAL2 expression.

Power-Law Modeling of Cancer Cell Fates Driven by Signaling Data to Reveal Drug Effects

Extracellular signals are captured and transmitted by signaling proteins inside a cell. An important type of cellular responses to the signals is the cell fate decision, e.g., apoptosis. However, the underlying mechanisms of cell fate regulation are still unclear, thus comprehensive and detailed kinetic models are not yet available. Alternatively, data-driven models are promising to bridge signaling data with the phenotypic measurements of cell fates. The traditional linear model for data-driven modeling of signaling pathways has its limitations because it assumes that the a cell fate is proportional to the activities of signaling proteins, which is unlikely in the complex biological systems. Therefore, we propose a power-law model to relate the activities of all the measured signaling proteins to the probabilities of cell fates. In our experiments, we compared our nonlinear power-law model with the linear model on three cancer datasets with phosphoproteomics and cell fate measurements, which demonstrated that the nonlinear model has superior performance on cell fates prediction. By in silico simulation of virtual protein knock-down, the proposed model is able to reveal drug effects which can complement traditional approaches such as binding affinity analysis. Moreover, our model is able to capture cell line specific information to distinguish one cell line from another in cell fate prediction. Our results show that the power-law data-driven model is able to perform better in cell fate prediction and provide more insights into the signaling pathways for cancer cell fates than the linear model.

Mathematical Models for Immunology: Current State of the Art and Future Research Directions

The advances in genetics and biochemistry that have taken place over the last 10 years led to significant advances in experimental and clinical immunology. In turn, this has led to the development of new mathematical models to investigate qualitatively and quantitatively various open questions in immunology. In this study we present a review of some research areas in mathematical immunology that evolved over the last 10 years. To this end, we take a step-by-step approach in discussing a range of models derived to study the dynamics of both the innate and immune responses at the molecular, cellular and tissue scales. To emphasise the use of mathematics in modelling in this area, we also review some of the mathematical tools used to investigate these models. Finally, we discuss some future trends in both experimental immunology and mathematical immunology for the upcoming years.

E2F and GATA switches turn off WD repeat domain 77 expression in differentiating cells

WDR77 (WD repeat domain 77) is expressed during earlier lung development when cells are rapidly proliferating, but is absent from adult lung. It is re-activated during lung tumorigenesis and is essential for lung cancer cell proliferation. Signalling pathways/molecules that control WDR77 gene expression are unknown. Promoter mapping, gel shift assay and ChIP revealed that the WDR77 promoter contains bona fide response elements for E2F and GATA transcriptional factors as demonstrated in prostate cancer, lung cancer and erythroid cells, as well as in mouse lung tissues. The WDR77 promoter is transactivated by E2F1, E2F3, GATA1 and GATA6, but suppressed by E2F6, GATA1 and GATA3 in prostate cancer PC3 cells. WDR77 expression is associated with E2F1, E2F3, GATA2 and GATA6 occupancy on the WDR77 gene, whereas, in contrast, E2F6, GATA1 and GATA3 occupancy is associated with the loss of WDR77 expression during erythroid maturation and lung development. More importantly, the loss of WDR77 expression that results from E2F and GATA switches is required for cellular differentiation of erythroid and lung epithelial cells. In contrast, lung cancer cells avoid post-mitotic differentiation by sustaining WDR77 expression. Altogether, the present study provides a novel molecular mechanism by which WDR77 is regulated during erythroid and lung development and lung tumorigenesis.

Emerging Modeling Concepts and Solutions in Stem Cell Research

Modern stem cell research, as well as other fields of contemporary biology involves quantitative sciences in many ways. Identifying candidates for key differentiation or reprogramming factors, tracing global transcriptome changes, or finding drugs is now broadly involves bioinformatics and biostatistics. However, the next key step, understanding the underlying reasons and establishing causal links leading to differentiation or reprogramming requires qualitative and quantitative biological models describing complex biological systems. Currently, quantitative modeling is a challenging science, capable to deliver rather modest results or predictions. What model types are the most popular and what features of stem cell behavior they are capturing? What new insights do we expect from the computational modeling of stem cells in the foreseeable future? Current review attempts to approach these essential questions by considering published quantitative models and solutions emerging in the area of stem cell research.

Generalized logical model based on network topology to capture the dynamical trends of cellular signaling pathways

Background

Cellular responses to extracellular perturbations require signaling pathways to capture and transmit the signals. However, the underlying molecular mechanisms of signal transduction are not yet fully understood, thus detailed and comprehensive models may not be available for all the signaling pathways. In particular, insufficient knowledge of parameters, which is a long-standing hindrance for quantitative kinetic modeling necessitates the use of parameter-free methods for modeling and simulation to capture dynamic properties of signaling pathways.

Results

We present a computational model that is able to simulate the graded responses to degradations, the sigmoidal biological relationships between signaling molecules and the effects of scheduled perturbations to the cells. The simulation results are validated using experimental data of protein phosphorylation, demonstrating that the proposed model is capable of capturing the main trend of protein activities during the process of signal transduction. Compared with existing simulators, our model has better performance on predicting the state transitions of signaling networks.

Conclusion

The proposed simulation tool provides a valuable resource for modeling cellular signaling pathways using a knowledge-based method.

Approximate Bayesian computation schemes for parameter inference of discrete stochastic models using simulated likelihood density

BACKGROUND

Mathematical modeling is an important tool in systems biology to study the dynamic property of complex biological systems. However, one of the major challenges in systems biology is how to infer unknown parameters in mathematical models based on the experimental data sets, in particular, when the data are sparse and the regulatory network is stochastic.

RESULTS

To address this issue, this work proposed a new algorithm to estimate parameters in stochastic models using simulated likelihood density in the framework of approximate Bayesian computation. Two stochastic models were used to demonstrate the efficiency and effectiveness of the proposed method. In addition, we designed another algorithm based on a novel objective function to measure the accuracy of stochastic simulations.

CONCLUSIONS

Simulation results suggest that the usage of simulated likelihood density improves the accuracy of estimates substantially. When the error is measured at each observation time point individually, the estimated parameters have better accuracy than those obtained by a published method in which the error is measured using simulations over the entire observation time period.

A continuous optimization approach for inferring parameters in mathematical models of regulatory networks

BACKGROUND

The advances of systems biology have raised a large number of sophisticated mathematical models for describing the dynamic property of complex biological systems. One of the major steps in developing mathematical models is to estimate unknown parameters of the model based on experimentally measured quantities. However, experimental conditions limit the amount of data that is available for mathematical modelling. The number of unknown parameters in mathematical models may be larger than the number of observation data. The imbalance between the number of experimental data and number of unknown parameters makes reverse-engineering problems particularly challenging.

RESULTS

To address the issue of inadequate experimental data, we propose a continuous optimization approach for making reliable inference of model parameters. This approach first uses a spline interpolation to generate continuous functions of system dynamics as well as the first and second order derivatives of continuous functions. The expanded dataset is the basis to infer unknown model parameters using various continuous optimization criteria, including the error of simulation only, error of both simulation and the first derivative, or error of simulation as well as the first and second derivatives. We use three case studies to demonstrate the accuracy and reliability of the proposed new approach. Compared with the corresponding discrete criteria using experimental data at the measurement time points only, numerical results of the ERK kinase activation module show that the continuous absolute-error criteria using both function and high order derivatives generate estimates with better accuracy. This result is also supported by the second and third case studies for the G1/S transition network and the MAP kinase pathway, respectively. This suggests that the continuous absolute-error criteria lead to more accurate estimates than the corresponding discrete criteria. We also study the robustness property of these three models to examine the reliability of estimates. Simulation results show that the models with estimated parameters using continuous fitness functions have better robustness properties than those using the corresponding discrete fitness functions.

CONCLUSIONS

The inference studies and robustness analysis suggest that the proposed continuous optimization criteria are effective and robust for estimating unknown parameters in mathematical models.