1
|
Yadav Y, Subbaroyan A, Martin OC, Samal A. Relative importance of composition structures and biologically meaningful logics in bipartite Boolean models of gene regulation. Sci Rep 2022; 12:18156. [PMID: 36307465 PMCID: PMC9616893 DOI: 10.1038/s41598-022-22654-7] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/14/2022] [Accepted: 10/18/2022] [Indexed: 12/31/2022] Open
Abstract
Boolean networks have been widely used to model gene networks. However, such models are coarse-grained to an extent that they abstract away molecular specificities of gene regulation. Alternatively, bipartite Boolean network models of gene regulation explicitly distinguish genes from transcription factors (TFs). In such bipartite models, multiple TFs may simultaneously contribute to gene regulation by forming heteromeric complexes, thus giving rise to composition structures. Since bipartite Boolean models are relatively recent, an empirical investigation of their biological plausibility is lacking. Here, we estimate the prevalence of composition structures arising through heteromeric complexes. Moreover, we present an additional mechanism where composition structures may arise as a result of multiple TFs binding to cis-regulatory regions and provide empirical support for this mechanism. Next, we compare the restriction in BFs imposed by composition structures and by biologically meaningful properties. We find that though composition structures can severely restrict the number of Boolean functions (BFs) driving a gene, the two types of minimally complex BFs, namely nested canalyzing functions (NCFs) and read-once functions (RoFs), are comparatively more restrictive. Finally, we find that composition structures are highly enriched in real networks, but this enrichment most likely comes from NCFs and RoFs.
Collapse
Affiliation(s)
- Yasharth Yadav
- The Institute of Mathematical Sciences (IMSc), Chennai, 600113, India
| | - Ajay Subbaroyan
- The Institute of Mathematical Sciences (IMSc), Chennai, 600113, India
- Homi Bhabha National Institute (HBNI), Mumbai, 400094, India
| | - Olivier C Martin
- Université Paris-Saclay, CNRS, INRAE, Univ Evry, Institute of Plant Sciences Paris-Saclay (IPS2), 91190, Gif sur Yvette, France.
- Université Paris Cité, CNRS, INRAE, Institute of Plant Sciences Paris-Saclay (IPS2), 91190, Gif sur Yvette, France.
| | - Areejit Samal
- The Institute of Mathematical Sciences (IMSc), Chennai, 600113, India.
- Homi Bhabha National Institute (HBNI), Mumbai, 400094, India.
| |
Collapse
|
2
|
Schnitzer B, Österberg L, Skopa I, Cvijovic M. Multi-scale model suggests the trade-off between protein and ATP demand as a driver of metabolic changes during yeast replicative ageing. PLoS Comput Biol 2022; 18:e1010261. [PMID: 35797415 PMCID: PMC9295998 DOI: 10.1371/journal.pcbi.1010261] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/11/2022] [Revised: 07/19/2022] [Accepted: 05/31/2022] [Indexed: 11/18/2022] Open
Abstract
The accumulation of protein damage is one of the major drivers of replicative ageing, describing a cell's reduced ability to reproduce over time even under optimal conditions. Reactive oxygen and nitrogen species are precursors of protein damage and therefore tightly linked to ageing. At the same time, they are an inevitable by-product of the cell's metabolism. Cells are able to sense high levels of reactive oxygen and nitrogen species and can subsequently adapt their metabolism through gene regulation to slow down damage accumulation. However, the older or damaged a cell is the less flexibility it has to allocate enzymes across the metabolic network, forcing further adaptions in the metabolism. To investigate changes in the metabolism during replicative ageing, we developed an multi-scale mathematical model using budding yeast as a model organism. The model consists of three interconnected modules: a Boolean model of the signalling network, an enzyme-constrained flux balance model of the central carbon metabolism and a dynamic model of growth and protein damage accumulation with discrete cell divisions. The model can explain known features of replicative ageing, like average lifespan and increase in generation time during successive division, in yeast wildtype cells by a decreasing pool of functional enzymes and an increasing energy demand for maintenance. We further used the model to identify three consecutive metabolic phases, that a cell can undergo during its life, and their influence on the replicative potential, and proposed an intervention span for lifespan control.
Collapse
Affiliation(s)
- Barbara Schnitzer
- Department of Mathematical Sciences, Chalmers University of Technology, Gothenburg, Sweden
- Department of Mathematical Sciences, University of Gothenburg, Gothenburg, Sweden
| | - Linnea Österberg
- Department of Mathematical Sciences, Chalmers University of Technology, Gothenburg, Sweden
- Department of Mathematical Sciences, University of Gothenburg, Gothenburg, Sweden
- Department of Biology and Biological Engineering, Chalmers University of Technology, Gothenburg, Sweden
| | - Iro Skopa
- Department of Mathematical Sciences, Chalmers University of Technology, Gothenburg, Sweden
- Department of Mathematical Sciences, University of Gothenburg, Gothenburg, Sweden
| | - Marija Cvijovic
- Department of Mathematical Sciences, Chalmers University of Technology, Gothenburg, Sweden
- Department of Mathematical Sciences, University of Gothenburg, Gothenburg, Sweden
| |
Collapse
|
3
|
Romers J, Thieme S, Münzner U, Krantz M. A scalable method for parameter-free simulation and validation of mechanistic cellular signal transduction network models. NPJ Syst Biol Appl 2020; 6:2. [PMID: 31934349 PMCID: PMC6954118 DOI: 10.1038/s41540-019-0120-5] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/18/2018] [Accepted: 11/20/2019] [Indexed: 11/09/2022] Open
Abstract
The metabolic modelling community has established the gold standard for bottom-up systems biology with reconstruction, validation and simulation of mechanistic genome-scale models. Similar methods have not been established for signal transduction networks, where the representation of complexes and internal states leads to scalability issues in both model formulation and execution. While rule- and agent-based methods allow efficient model definition and execution, respectively, model parametrisation introduces an additional layer of uncertainty due to the sparsity of reliably measured parameters. Here, we present a scalable method for parameter-free simulation of mechanistic signal transduction networks. It is based on rxncon and uses a bipartite Boolean logic with separate update rules for reactions and states. Using two generic update rules, we enable translation of any rxncon model into a unique Boolean model, which can be used for network validation and simulation-allowing the prediction of system-level function directly from molecular mechanistic data. Through scalable model definition and simulation, and the independence of quantitative parameters, it opens up for simulation and validation of mechanistic genome-scale models of signal transduction networks.
Collapse
Affiliation(s)
- Jesper Romers
- Institute for Biology, Humboldt-Universität zu Berlin, Berlin, Germany
| | - Sebastian Thieme
- Institute for Biology, Humboldt-Universität zu Berlin, Berlin, Germany
| | - Ulrike Münzner
- Institute for Biology, Humboldt-Universität zu Berlin, Berlin, Germany
- Bioinformatics Center, Institute for Chemical Research, Kyoto University, Uji, Japan
| | - Marcus Krantz
- Institute for Biology, Humboldt-Universität zu Berlin, Berlin, Germany
| |
Collapse
|
4
|
Abstract
We present a protocol for building, validating, and simulating models of signal transduction networks. These networks are challenging modeling targets due to the combinatorial complexity and sparse data, which have made it a major challenge even to formalize the current knowledge. To address this, the community has developed methods to model biomolecular reaction networks based on site dynamics. The strength of this approach is that reactions and states can be defined at variable resolution, which makes it possible to adapt the model resolution to the empirical data. This improves both scalability and accuracy, making it possible to formalize large models of signal transduction networks. Here, we present a method to build and validate large models of signal transduction networks. The workflow is based on rxncon, the reaction-contingency language. In a five-step process, we create a mechanistic network model, convert it into an executable Boolean model, use the Boolean model to evaluate and improve the network, and finally export the rxncon model into a rule-based format. We provide an introduction to the rxncon language and an annotated, step-by-step protocol for the workflow. Finally, we create a small model of the insulin signaling pathway to illustrate the protocol, together with some of the challenges-and some of their solutions-in modeling signal transduction.
Collapse
Affiliation(s)
- Jesper Romers
- Institute of Biology, Humboldt-Universität zu Berlin, Berlin, Germany
| | - Sebastian Thieme
- Institute of Biology, Humboldt-Universität zu Berlin, Berlin, Germany
| | - Ulrike Münzner
- Institute of Biology, Humboldt-Universität zu Berlin, Berlin, Germany.,Bioinformatics Center, Institute for Chemical Research, Kyoto University, Uji, Japan
| | - Marcus Krantz
- Institute of Biology, Humboldt-Universität zu Berlin, Berlin, Germany.
| |
Collapse
|
5
|
Chaves M, Tournier L. Analysis Tools for Interconnected Boolean Networks With Biological Applications. Front Physiol 2018; 9:586. [PMID: 29896108 PMCID: PMC5987301 DOI: 10.3389/fphys.2018.00586] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2018] [Accepted: 05/02/2018] [Indexed: 11/13/2022] Open
Abstract
Boolean networks with asynchronous updates are a class of logical models particularly well adapted to describe the dynamics of biological networks with uncertain measures. The state space of these models can be described by an asynchronous state transition graph, which represents all the possible exits from every single state, and gives a global image of all the possible trajectories of the system. In addition, the asynchronous state transition graph can be associated with an absorbing Markov chain, further providing a semi-quantitative framework where it becomes possible to compute probabilities for the different trajectories. For large networks, however, such direct analyses become computationally untractable, given the exponential dimension of the graph. Exploiting the general modularity of biological systems, we have introduced the novel concept of asymptotic graph, computed as an interconnection of several asynchronous transition graphs and recovering all asymptotic behaviors of a large interconnected system from the behavior of its smaller modules. From a modeling point of view, the interconnection of networks is very useful to address for instance the interplay between known biological modules and to test different hypotheses on the nature of their mutual regulatory links. This paper develops two new features of this general methodology: a quantitative dimension is added to the asymptotic graph, through the computation of relative probabilities for each final attractor and a companion cross-graph is introduced to complement the method on a theoretical point of view.
Collapse
Affiliation(s)
- Madalena Chaves
- Inria Sophia Antipolis - Méditerranée, Université Côte d'Azur, Valbonne, France
| | | |
Collapse
|
6
|
Bakker E, Tian K, Mutti L, Demonacos C, Schwartz JM, Krstic-Demonacos M. Insight into glucocorticoid receptor signalling through interactome model analysis. PLoS Comput Biol 2017; 13:e1005825. [PMID: 29107989 PMCID: PMC5690696 DOI: 10.1371/journal.pcbi.1005825] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2017] [Revised: 11/16/2017] [Accepted: 10/16/2017] [Indexed: 12/27/2022] Open
Abstract
Glucocorticoid hormones (GCs) are used to treat a variety of diseases because of their potent anti-inflammatory effect and their ability to induce apoptosis in lymphoid malignancies through the glucocorticoid receptor (GR). Despite ongoing research, high glucocorticoid efficacy and widespread usage in medicine, resistance, disease relapse and toxicity remain factors that need addressing. Understanding the mechanisms of glucocorticoid signalling and how resistance may arise is highly important towards improving therapy. To gain insight into this we undertook a systems biology approach, aiming to generate a Boolean model of the glucocorticoid receptor protein interaction network that encapsulates functional relationships between the GR, its target genes or genes that target GR, and the interactions between the genes that interact with the GR. This model named GEB052 consists of 52 nodes representing genes or proteins, the model input (GC) and model outputs (cell death and inflammation), connected by 241 logical interactions of activation or inhibition. 323 changes in the relationships between model constituents following in silico knockouts were uncovered, and steady-state analysis followed by cell-based microarray genome-wide model validation led to an average of 57% correct predictions, which was taken further by assessment of model predictions against patient microarray data. Lastly, semi-quantitative model analysis via microarray data superimposed onto the model with a score flow algorithm has also been performed, which demonstrated significantly higher correct prediction ratios (average of 80%), and the model has been assessed as a predictive clinical tool using published patient microarray data. In summary we present an in silico simulation of the glucocorticoid receptor interaction network, linked to downstream biological processes that can be analysed to uncover relationships between GR and its interactants. Ultimately the model provides a platform for future development both by directing laboratory research and allowing for incorporation of further components, encapsulating more interactions/genes involved in glucocorticoid receptor signalling. Here we present modelling of the glucocorticoid receptor (GR) signalling network. The GR is the effector for a class of drugs known as corticosteroids, which are widely used in medicine for their anti-inflammatory effects and ability to induce apoptosis in leukaemic cells. However, side effects, treatment-related toxicity and glucocorticoid resistance remain and therefore increased understanding of the glucocorticoid receptor mechanism of action may improve therapeutic outcomes. The GEB052 model presented herein has been used to generate predictions for how the network is altered between glucocorticoid-sensitive and glucocorticoid-resistant scenarios, and these predictions have been verified using published gene expression data from established cell lines (for both qualitative and semi-quantitative analysis). The model has also been preliminarily assessed as a predictive clinical tool by correlating model predictions with clinical outcomes of thirteen leukaemia patients. Thus, the GEB052 model demonstrates successful modelling to understand GR function. GEB052 provides accurate predictions and has indicated potential routes through which glucocorticoid resistance may arise. The work presented herein thus demonstrates a proof-of-principle of this modelling approach to furthering GR research, and provides insight into potential mechanisms of corticosteroids resistance.
Collapse
Affiliation(s)
- Emyr Bakker
- Biomedical Research Centre, School of Environment and Life Sciences, University of Salford, Salford, United Kingdom
| | - Kun Tian
- Biomedical Research Centre, School of Environment and Life Sciences, University of Salford, Salford, United Kingdom
| | - Luciano Mutti
- Biomedical Research Centre, School of Environment and Life Sciences, University of Salford, Salford, United Kingdom
| | - Constantinos Demonacos
- Faculty of Biology, Medicine and Health, University of Manchester, Manchester, United Kingdom
| | - Jean-Marc Schwartz
- Faculty of Biology, Medicine and Health, University of Manchester, Manchester, United Kingdom
- * E-mail: (JMS); (MKD)
| | - Marija Krstic-Demonacos
- Biomedical Research Centre, School of Environment and Life Sciences, University of Salford, Salford, United Kingdom
- * E-mail: (JMS); (MKD)
| |
Collapse
|
7
|
Mathematical and Computational Modeling in Complex Biological Systems. BIOMED RESEARCH INTERNATIONAL 2017; 2017:5958321. [PMID: 28386558 PMCID: PMC5366773 DOI: 10.1155/2017/5958321] [Citation(s) in RCA: 31] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 10/16/2016] [Revised: 12/20/2016] [Accepted: 01/16/2017] [Indexed: 12/22/2022]
Abstract
The biological process and molecular functions involved in the cancer progression remain difficult to understand for biologists and clinical doctors. Recent developments in high-throughput technologies urge the systems biology to achieve more precise models for complex diseases. Computational and mathematical models are gradually being used to help us understand the omics data produced by high-throughput experimental techniques. The use of computational models in systems biology allows us to explore the pathogenesis of complex diseases, improve our understanding of the latent molecular mechanisms, and promote treatment strategy optimization and new drug discovery. Currently, it is urgent to bridge the gap between the developments of high-throughput technologies and systemic modeling of the biological process in cancer research. In this review, we firstly studied several typical mathematical modeling approaches of biological systems in different scales and deeply analyzed their characteristics, advantages, applications, and limitations. Next, three potential research directions in systems modeling were summarized. To conclude, this review provides an update of important solutions using computational modeling approaches in systems biology.
Collapse
|
8
|
Abstract
Molecular profiling of proteins and phosphoproteins using a reverse phase protein array (RPPA) platform, with a panel of target-specific antibodies, enables the parallel, quantitative proteomic analysis of many biological samples in a microarray format. Hence, RPPA analysis can generate a high volume of multidimensional data that must be effectively interrogated and interpreted. A range of computational techniques for data mining can be applied to detect and explore data structure and to form functional predictions from large datasets. Here, two approaches for the computational analysis of RPPA data are detailed: the identification of similar patterns of protein expression by hierarchical cluster analysis and the modeling of protein interactions and signaling relationships by network analysis. The protocols use freely available, cross-platform software, are easy to implement, and do not require any programming expertise. Serving as data-driven starting points for further in-depth analysis, validation, and biological experimentation, these and related bioinformatic approaches can accelerate the functional interpretation of RPPA data.
Collapse
Affiliation(s)
- Adam Byron
- Cancer Research UK Edinburgh Centre, Institute of Genetics and Molecular Medicine, University of Edinburgh, Edinburgh, EH4 2XR, UK.
| |
Collapse
|
9
|
Cheng X, Mori T, Qiu Y, Ching WK, Akutsu T. Exact Identification of the Structure of a Probabilistic Boolean Network from Samples. IEEE/ACM TRANSACTIONS ON COMPUTATIONAL BIOLOGY AND BIOINFORMATICS 2016; 13:1107-1116. [PMID: 26661790 DOI: 10.1109/tcbb.2015.2505310] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Abstract
We study the number of samples required to uniquely determine the structure of a probabilistic Boolean network (PBN), where PBNs are probabilistic extensions of Boolean networks. We show via theoretical analysis and computational analysis that the structure of a PBN can be exactly identified with high probability from a relatively small number of samples for interesting classes of PBNs of bounded indegree. On the other hand, we also show that there exist classes of PBNs for which it is impossible to uniquely determine the structure of a PBN from samples.
Collapse
|
10
|
Zhang F, Chen H, Zhao LN, Liu H, Przytycka TM, Zheng J. Generalized logical model based on network topology to capture the dynamical trends of cellular signaling pathways. BMC SYSTEMS BIOLOGY 2016; 10 Suppl 1:7. [PMID: 26818802 PMCID: PMC4895646 DOI: 10.1186/s12918-015-0249-9] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
Abstract
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.
Collapse
Affiliation(s)
- Fan Zhang
- Biomedical Informatics Graduate Lab, School of Computer Engineering, Nanyang Technological University, Singapore 639798, Singapore.
| | - Haoting Chen
- Biomedical Informatics Graduate Lab, School of Computer Engineering, Nanyang Technological University, Singapore 639798, Singapore. .,Department of Industrial Engineering and Operations Research, Columbia University, New York, NY 10027, USA.
| | - Li Na Zhao
- Bioinformatics Institute, Agency for Science, Technology and Research, Singapore 138671, Singapore.
| | - Hui Liu
- Biomedical Informatics Graduate Lab, School of Computer Engineering, Nanyang Technological University, Singapore 639798, Singapore. .,Lab of Information Management, Changzhou University, Changzhou, Jiangsu 213164, China.
| | - Teresa M Przytycka
- National Center for Biotechnology Information, NLM/NIH, Bethesda, MD 20894, USA.
| | - Jie Zheng
- Biomedical Informatics Graduate Lab, School of Computer Engineering, Nanyang Technological University, Singapore 639798, Singapore. .,Complexity Institute, Nanyang Technological University, Singapore 637723, Singapore. .,Genome Institute of Singapore, Agency for Science, Technology and Research, Singapore 138672, Singapore.
| |
Collapse
|