Random Parametric Perturbations of Gene Regulatory Circuit Uncover State Transitions in Cell Cycle

Data-driven modeling of core gene regulatory network underlying leukemogenesis in IDH mutant AML

Acute myeloid leukemia (AML) is characterized by uncontrolled proliferation of poorly differentiated myeloid cells, with a heterogenous mutational landscape. Mutations in IDH1 and IDH2 are found in 20% of the AML cases. Although much effort has been made to identify genes associated with leukemogenesis, the regulatory mechanism of AML state transition is still not fully understood. To alleviate this issue, here we develop a new computational approach that integrates genomic data from diverse sources, including gene expression and ATAC-seq datasets, curated gene regulatory interaction databases, and mathematical modeling to establish models of context-specific core gene regulatory networks (GRNs) for a mechanistic understanding of tumorigenesis of AML with IDH mutations. The approach adopts a new optimization procedure to identify the top network according to its accuracy in capturing gene expression states and its flexibility to allow sufficient control of state transitions. From GRN modeling, we identify key regulators associated with the function of IDH mutations, such as DNA methyltransferase DNMT1, and network destabilizers, such as E2F1. The constructed core regulatory network and outcomes of in-silico network perturbations are supported by survival data from AML patients. We expect that the combined bioinformatics and systems-biology modeling approach will be generally applicable to elucidate the gene regulation of disease progression.

Constructing maps between distinct cell fates and parametric conditions by systematic perturbations

MOTIVATION

Cell fate transitions are common in many developmental processes. Therefore, identifying the mechanisms behind them is crucial. Traditionally, due to complexity of networks and existence of plenty of kinetic parameters, dynamical analysis of biomolecular networks can only be performed by simultaneously perturbing a small number of parameters. Although many efforts have focused on how cell states change under specific perturbations, conversely, how to infer parametric conditions underlying distinct cell fates by systematic perturbations is less clear and needs to be further investigated.

RESULTS

In this article, we present a general computational method by integrating systematic perturbations, unsupervised clustering, principal component analysis, and fitting analysis. The method can be used to to construct maps between distinct cell fates and parametric conditions by systematic perturbations. In particular, there are no needs of accurate parameter measurements and occurrence of bifurcations to establish the maps. To validate feasibility and inference performance of the method, we use toggle switch, inner cell mass, and epithelial mesenchymal transition as model systems to show how the maps are constructed and how system parameters encode essential information on cell fates. The maps tell us how systematic perturbations drive cell fate decisions and transitions, and allow us to purposefully predict, manipulate, and even control cell states. The approach is especially helpful in understanding crucial roles of certain parameter combinations during fate transitions. We hope that the approach can provide us valuable information on parametric or perturbation conditions so some specific targets, e.g. directional differentiation, can be realized.

AVAILABILITY AND IMPLEMENTATION

No public data are used. The data we used are generated by randomly chosen values of model parameters in certain ranges, and the corresponding parameters are already attached in supplementary materials.

Data-driven modeling of core gene regulatory network underlying leukemogenesis in IDH mutant AML

Acute myeloid leukemia (AML) is characterized by uncontrolled proliferation of poorly differentiated myeloid cells, with a heterogenous mutational landscape. Mutations in IDH1 and IDH2 are found in 20% of the AML cases. Although much effort has been made to identify genes associated with leukemogenesis, the regulatory mechanism of AML state transition is still not fully understood. To alleviate this issue, here we develop a new computational approach that integrates genomic data from diverse sources, including gene expression and ATAC-seq datasets, curated gene regulatory interaction databases, and mathematical modeling to establish models of context-specific core gene regulatory networks (GRNs) for a mechanistic understanding of tumorigenesis of AML with IDH mutations. The approach adopts a novel optimization procedure to identify the optimal network according to its accuracy in capturing gene expression states and its flexibility to allow sufficient control of state transitions. From GRN modeling, we identify key regulators associated with the function of IDH mutations, such as DNA methyltransferase DNMT1, and network destabilizers, such as E2F1. The constructed core regulatory network and outcomes of in-silico network perturbations are supported by survival data from AML patients. We expect that the combined bioinformatics and systems-biology modeling approach will be generally applicable to elucidate the gene regulation of disease progression.

Computational analysis of synergism in small networks with different logic

Cell fate decision processes are regulated by networks which contain different molecules and interactions. Different network topologies may exhibit synergistic or antagonistic effects on cellular functions. Here, we analyze six most common small networks with regulatory logic AND or OR, trying to clarify the relationship between network topologies and synergism (or antagonism) related to cell fate decisions. We systematically examine the contribution of both network topologies and regulatory logic to the cell fate synergism by bifurcation and combinatorial perturbation analysis. Initially, under a single set of parameters, the synergism of three types of networks with AND and OR logic is compared. Furthermore, to consider whether these results depend on the choices of parameter values, statistics on the synergism of five hundred parameter sets is performed. It is shown that the results are not sensitive to parameter variations, indicating that the synergy or antagonism mainly depends on the network topologies rather than the choices of parameter values. The results indicate that the topology with "Dual Inhibition" shows good synergism, while the topology with "Dual Promotion" or "Hybrid" shows antagonism. The results presented here may help us to design synergistic networks based on network structure and regulation combinations, which has promising implications for cell fate decisions and drug combinations.

A quantitative evaluation of topological motifs and their coupling in gene circuit state distributions

One of the major challenges in biology is to understand how gene interactions collaborate to determine overall functions of biological systems. Here, we present a new computational framework that enables systematic, high-throughput, and quantitative evaluation of how small transcriptional regulatory circuit motifs, and their coupling, contribute to functions of a dynamical biological system. We illustrate how this approach can be applied to identify four-node gene circuits, circuit motifs, and motif coupling responsible for various gene expression state distributions, including those derived from single-cell RNA sequencing data. We also identify seven major classes of four-node circuits from clustering analysis of state distributions. The method is applied to establish phenomenological models of gene circuits driving human neuron differentiation, revealing important biologically relevant regulatory interactions. Our study will shed light on a better understanding of gene regulatory mechanisms in creating and maintaining cellular states.

A data-driven optimization method for coarse-graining gene regulatory networks

One major challenge in systems biology is to understand how various genes in a gene regulatory network (GRN) collectively perform their functions and control network dynamics. This task becomes extremely hard to tackle in the case of large networks with hundreds of genes and edges, many of which have redundant regulatory roles and functions. The existing methods for model reduction usually require the detailed mathematical description of dynamical systems and their corresponding kinetic parameters, which are often not available. Here, we present a data-driven method for coarse-graining large GRNs, named SacoGraci, using ensemble-based mathematical modeling, dimensionality reduction, and gene circuit optimization by Markov Chain Monte Carlo methods. SacoGraci requires network topology as the only input and is robust against errors in GRNs. We benchmark and demonstrate its usage with synthetic, literature-based, and bioinformatics-derived GRNs. We hope SacoGraci will enhance our ability to model the gene regulation of complex biological systems.

NetAct: a computational platform to construct core transcription factor regulatory networks using gene activity

A major question in systems biology is how to identify the core gene regulatory circuit that governs the decision-making of a biological process. Here, we develop a computational platform, named NetAct, for constructing core transcription factor regulatory networks using both transcriptomics data and literature-based transcription factor-target databases. NetAct robustly infers regulators' activity using target expression, constructs networks based on transcriptional activity, and integrates mathematical modeling for validation. Our in silico benchmark test shows that NetAct outperforms existing algorithms in inferring transcriptional activity and gene networks. We illustrate the application of NetAct to model networks driving TGF-β-induced epithelial-mesenchymal transition and macrophage polarization.

What Makes a Functional Gene Regulatory Network? A Circuit Motif Analysis

One of the key questions in systems biology is to understand the roles of gene regulatory circuits in determining cellular states and their functions. In previous studies, some researchers have inferred large gene networks from genome wide genomics/transcriptomics data using the top-down approach, while others have modeled core gene circuits of small sizes using the bottom-up approach. Despite many existing systems biology studies, there is still no general rule on what sizes of gene networks and what types of circuit motifs a system would need to achieve robust biological functions. Here, we adopt a gene circuit motif analysis to discover four-node circuits responsible for multiplicity (rich in dynamical behavior), flexibility (versatile to alter gene expression), or both. We identify the most reoccurring two-node circuit motifs and the co-occurring motif pairs. Furthermore, we investigate the contributing factors of multiplicity and flexibility for large gene networks of different types and sizes. We find that gene networks of intermediate sizes tend to have combined high levels of multiplicity and flexibility. Our study will contribute to a better understanding of the dynamical mechanisms of gene regulatory circuits and provide insights into rational designs of robust gene circuits in synthetic and systems biology.

Quantifying Mutational Response to Track the Evolution of SARS-CoV-2 Spike Variants: Introducing a Statistical-Mechanics-Guided Machine Learning Method

The emergence of SARS-CoV-2 and its variants that critically affect global public health requires characterization of mutations and their evolutionary pattern from specific Variants of Interest (VOIs) to Variants of Concern (VOCs). Leveraging the concept of equilibrium statistical mechanics, we introduce a new responsive quantity defined as "Mutational Response Function (MRF)" aptly quantifying domain-wise average entropy-fluctuation in the spike glycoprotein sequence of SARS-CoV-2 based on its evolutionary database. As the evolution transits from a specific variant to VOC, we find that the evolutionary crossover is accompanied by a dramatic change in MRF, upholding the characteristic of a dynamic phase transition. With this entropic information, we have developed an ancestral-based machine learning method that helps predict future domain-specific mutations. The feedforward binary classification model pinpoints possible residues prone to future mutations that have implications for enhanced fusogenicity and pathogenicity of the virus. We believe such MRF analyses followed by a statistical mechanics augmented ML approach could help track different evolutionary stages of such species and identify a critical evolutionary transition that is alarming.

Gene Circuit Explorer (GeneEx): an interactive web-app for visualizing, simulating and analyzing gene regulatory circuits

SUMMARY

GeneEx is an interactive web-app that uses an ODE-based mathematical modeling approach to simulate, visualize and analyze gene regulatory circuits (GRCs) for an explicit kinetic parameter set or for a large ensemble of random parameter sets. GeneEx offers users the freedom to modify many aspects of the simulation such as the parameter ranges, the levels of gene expression noise and the GRC network topology itself. This degree of flexibility allows users to explore a variety of hypotheses by providing insight into the number and stability of attractors for a given GRC. Moreover, users have the option to upload, and subsequently compare, experimental gene expression data to simulated data generated from the analysis of a built or uploaded custom circuit. Finally, GeneEx offers a curated database that contains circuit motifs and known biological GRCs to facilitate further inquiry into these. Overall, GeneEx enables users to investigate the effects of parameter variation, stochasticity and/or topological changes on gene expression for GRCs using a systems-biology approach.

AVAILABILITY AND IMPLEMENTATION

GeneEx is available at https://geneex.jax.org. This web-app is released under the MIT license and is free and open to all users and there is no mandatory login requirement.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.