1
|
Troják M, Šafránek D, Pastva S, Brim L. Rule-based modelling of biological systems using regulated rewriting. Biosystems 2023; 225:104843. [PMID: 36736686 DOI: 10.1016/j.biosystems.2023.104843] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2022] [Revised: 01/03/2023] [Accepted: 01/24/2023] [Indexed: 02/04/2023]
Abstract
In systems biology, models play a crucial role in understanding studied systems. There are many modelling approaches, among which rewriting systems provide a framework for describing systems on a mechanistic level. Describing biochemical processes often requires incorporating knowledge on an abstract level to simplify the system description or substitute the missing details. For this purpose, we present regulation mechanisms, an extension of this formalism with additional controls on the rewriting process. We introduce several regulation mechanisms and apply them to a rule-based language, a notation suitable for modelling biological phenomena. Finally, we demonstrate the usage of such regulations on several case studies from the biochemical domain.
Collapse
Affiliation(s)
- Matej Troják
- Systems Biology Laboratory, Masaryk University, Brno, Czech Republic.
| | - David Šafránek
- Systems Biology Laboratory, Masaryk University, Brno, Czech Republic
| | - Samuel Pastva
- Systems Biology Laboratory, Masaryk University, Brno, Czech Republic
| | - Luboš Brim
- Systems Biology Laboratory, Masaryk University, Brno, Czech Republic
| |
Collapse
|
2
|
Alsharaiah MA, Samarasinghe S, Kulasiri D. Proteins as fuzzy controllers: Auto tuning a biological fuzzy inference system to predict protein dynamics in complex biological networks. Biosystems 2023; 224:104826. [PMID: 36610587 DOI: 10.1016/j.biosystems.2023.104826] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/05/2022] [Revised: 11/30/2022] [Accepted: 01/02/2023] [Indexed: 01/06/2023]
Abstract
Biological systems such as mammalian cell cycle are complex systems consisting of a large number of molecular species interacting in ways that produce complex nonlinear systems dynamics. Discrete models such as Boolean models and continuous models such as Ordinary Differential Equations (ODEs) have been widely used to study these systems. Boolean models are simple and can capture qualitative systems behaviour, but they cannot capture the continuous trends of protein concentrations, while ODE models capture continuous trends but require kinetics parameters that are limited. Further, as systems get larger, complexity of these models becomes an issue for parameterization, analysis and interpretation. Also, molecular systems operate under the conditions of uncertainty and noise and our understanding of molecular processes in general is more at a qualitative level characterised by vagueness, imprecision and ambiguity. Hence, as more data are generated, there is a greater need for simpler data driven methods that can approximate continuous system behaviour while representing vagueness and ambiguity without requiring kinetic parameters. Fuzzy inferencing is one such promising method with the ability to work with qualitative vague/imprecise biological knowledge. In this study, we propose a fuzzy inference system for representing continuous behaviour of proteins and apply to some key proteins in the mammalian cell cycle system. The methods we introduced here is novel to protein interaction systems and cell cycle proteins. Our study proposes a three-stage approach to develop fuzzy protein controllers. In stage one, protein system is studied for interactions. We studied some significant core controllers of mammalian cell cycle and their producers and degraders as presented in a published ODE model. Based on the observations from a dataset generated from it, we developed Fuzzy inference systems (FIS) in the second stage, that involved deriving fuzzy IF-THEN rules and their processing, and manually tuned the FIS to predict the dynamics of individual proteins. In stage three, we employed Particle Swarm Optimisation (PSO) for optimising the FIS to further enhance prediction accuracy. Systems dynamics simulation results of the optimised FIS models were in close agreement with the benchmark ODE model results. The results show that the FIS models provide a close approximation to the comprehensive benchmark model in robustly representing continuous protein dynamics while representing the control of protein behavior in an intuitive and transparent format without requiring kinetic parameters. Therefore, FIS models can be an alternative to ODEs in network modelling. Further, FIS models can be assembled to develop large complex systems without losing information or accuracy.
Collapse
Affiliation(s)
| | - Sandhya Samarasinghe
- Complex Systems, Big Data and Informatics Initiative (CSBII), Lincoln University, Christchurch, New Zealand; Centre for Advanced Computational Solutions, Lincoln University, Christchurch, New Zealand.
| | - Don Kulasiri
- Complex Systems, Big Data and Informatics Initiative (CSBII), Lincoln University, Christchurch, New Zealand; Centre for Advanced Computational Solutions, Lincoln University, Christchurch, New Zealand
| |
Collapse
|
3
|
Abroudi A, Samarasinghe S, Kulasiri D. Towards abstraction of computational modelling of mammalian cell cycle: Model reduction pipeline incorporating multi-level hybrid petri nets. J Theor Biol 2020; 496:110212. [PMID: 32142804 DOI: 10.1016/j.jtbi.2020.110212] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/03/2019] [Revised: 12/13/2019] [Accepted: 02/23/2020] [Indexed: 12/31/2022]
Abstract
Cell cycle is a large biochemical network and it is crucial to simplify it to gain a clearer understanding and insights into the cell cycle. This is also true for other biochemical networks. In this study, we present a model abstraction scheme/pipeline to create a minimal abstract model of the whole mammalian cell cycle system from a large Ordinary Differential Equation model of cell cycle we published previously (Abroudi et al., 2017). The abstract model is developed in a way that it captures the main characteristics (dynamics of key controllers), responses (G1-S and G2-M transitions and DNA damage) and the signalling subsystems (Growth Factor, G1-S and G2-M checkpoints, and DNA damage) of the original model (benchmark). Further, our model exploits: (i) separation of time scales (slow and fast reactions), (ii) separation of levels of complexity (high-level and low-level interactions), (iii) cell-cycle stages (temporality), (iv) functional subsystems (as mentioned above), and (v) represents the whole cell cycle - within a Multi-Level Hybrid Petri Net (MLHPN) framework. Although hybrid Petri Nets is not new, the abstraction of interactions and timing we introduced here is new to cell cycle and Petri Nets. Importantly, our models builds on the significant elements, representing the core cell cycle system, found through a novel Global Sensitivity Analysis on the benchmark model, using Self Organising Maps and Correlation Analysis that we introduced in (Abroudi et al., 2017). Taken the two aspects together, our study proposes a 2-stage model reduction pipeline for large systems and the main focus of this paper is on stage 2, Petri Net model, put in the context of the pipeline. With the MLHPN model, the benchmark model with 61 continuous variables (ODEs) and 148 parameters were reduced to 14 variables (4 continuous (Cyc_Cdks - the main controllers of cell cycle) and 10 discrete (regulators of Cyc_Cdks)) and 31 parameters. Additional 9 discrete elements represented the temporal progression of cell cycle. Systems dynamics simulation results of the MLHPN model were in close agreement with the benchmark model with respect to the crucial metrics selected for comparison: order and pattern of Cyc_Cdk activation, timing of G1-S and G2-M transitions with or without DNA damage, efficiency of the two cell cycle checkpoints in arresting damaged cells and passing healthy cells, and response to two types of global parameter perturbations. The results show that the MLHPN provides a close approximation to the comprehensive benchmark model in robustly representing systems dynamics and emergent properties while presenting the core cell cycle controller in an intuitive, transparent and subsystems format.
Collapse
Affiliation(s)
- Ali Abroudi
- Complex Systems, Big Data and Informatics Initiative (CSBII), Lincoln University, New Zealand
| | - Sandhya Samarasinghe
- Complex Systems, Big Data and Informatics Initiative (CSBII), Lincoln University, New Zealand.
| | - Don Kulasiri
- Complex Systems, Big Data and Informatics Initiative (CSBII), Lincoln University, New Zealand
| |
Collapse
|
4
|
Christ B, Dahmen U, Herrmann KH, König M, Reichenbach JR, Ricken T, Schleicher J, Ole Schwen L, Vlaic S, Waschinsky N. Computational Modeling in Liver Surgery. Front Physiol 2017; 8:906. [PMID: 29249974 PMCID: PMC5715340 DOI: 10.3389/fphys.2017.00906] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2017] [Accepted: 10/25/2017] [Indexed: 12/13/2022] Open
Abstract
The need for extended liver resection is increasing due to the growing incidence of liver tumors in aging societies. Individualized surgical planning is the key for identifying the optimal resection strategy and to minimize the risk of postoperative liver failure and tumor recurrence. Current computational tools provide virtual planning of liver resection by taking into account the spatial relationship between the tumor and the hepatic vascular trees, as well as the size of the future liver remnant. However, size and function of the liver are not necessarily equivalent. Hence, determining the future liver volume might misestimate the future liver function, especially in cases of hepatic comorbidities such as hepatic steatosis. A systems medicine approach could be applied, including biological, medical, and surgical aspects, by integrating all available anatomical and functional information of the individual patient. Such an approach holds promise for better prediction of postoperative liver function and hence improved risk assessment. This review provides an overview of mathematical models related to the liver and its function and explores their potential relevance for computational liver surgery. We first summarize key facts of hepatic anatomy, physiology, and pathology relevant for hepatic surgery, followed by a description of the computational tools currently used in liver surgical planning. Then we present selected state-of-the-art computational liver models potentially useful to support liver surgery. Finally, we discuss the main challenges that will need to be addressed when developing advanced computational planning tools in the context of liver surgery.
Collapse
Affiliation(s)
- Bruno Christ
- Molecular Hepatology Lab, Clinics of Visceral, Transplantation, Thoracic and Vascular Surgery, University Hospital Leipzig, University of Leipzig, Leipzig, Germany
| | - Uta Dahmen
- Experimental Transplantation Surgery, Department of General, Visceral and Vascular Surgery, University Hospital Jena, Jena, Germany
| | - Karl-Heinz Herrmann
- Medical Physics Group, Institute for Diagnostic and Interventional Radiology, University Hospital Jena, Friedrich Schiller University Jena, Jena, Germany
| | - Matthias König
- Department of Biology, Institute for Theoretical Biology, Humboldt University of Berlin, Berlin, Germany
| | - Jürgen R Reichenbach
- Medical Physics Group, Institute for Diagnostic and Interventional Radiology, University Hospital Jena, Friedrich Schiller University Jena, Jena, Germany
| | - Tim Ricken
- Mechanics, Structural Analysis, and Dynamics, TU Dortmund University, Dortmund, Germany
| | - Jana Schleicher
- Experimental Transplantation Surgery, Department of General, Visceral and Vascular Surgery, University Hospital Jena, Jena, Germany.,Department of Bioinformatics, Friedrich Schiller University Jena, Jena, Germany
| | | | - Sebastian Vlaic
- Leibniz Institute for Natural Product Research and Infection Biology, Hans Knöll Institute, Jena, Germany
| | - Navina Waschinsky
- Mechanics, Structural Analysis, and Dynamics, TU Dortmund University, Dortmund, Germany
| |
Collapse
|
5
|
Martin O, Krzywicki A, Zagorski M. Drivers of structural features in gene regulatory networks: From biophysical constraints to biological function. Phys Life Rev 2016; 17:124-58. [DOI: 10.1016/j.plrev.2016.06.002] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2015] [Revised: 03/25/2016] [Accepted: 04/20/2016] [Indexed: 12/23/2022]
|