ODE Models I: Biological Systems

2018 ◽  
Author(s):  
Yasmin Z. Paterson ◽  
David Shorthouse ◽  
Markus W. Pleijzier ◽  
Nir Piterman ◽  
Claus Bendtsen ◽  
...  

ABSTRACTIn an age where the volume of data regarding biological systems exceeds our ability to analyse it, many researchers are looking towards systems biology and computational modelling to help unravel the complexities of gene and protein regulatory networks. In particular, the use of discrete modelling allows generation of signalling networks in the absence of full quantitative descriptions of systems, which are necessary for ordinary differential equation (ODE) models. In order to make such techniques more accessible to mainstream researchers, tools such as the BioModelAnalyzer (BMA) have been developed to provide a user-friendly graphical interface for discrete modelling of biological systems. Here we use the BMA to build a library of discrete target functions of known canonical molecular interactions, translated from ordinary differential equations (ODEs). We then show that these BMA target functions can be used to reconstruct complex networks, which can correctly predict many known genetic perturbations. This new library supports the accessibility ethos behind the creation of BMA, providing a toolbox for the construction of complex cell signalling models without the need for extensive experience in computer programming or mathematical modelling, and allows for construction and simulation of complex biological systems with only small amounts of quantitative data.AUTHOR SUMMARYOrdinary differential equation (ODE) based models are a popular approach for modelling biological networks. A limitation of ODE models is that they require complete networks and detailed kinetic parameterisation. An alternative is the use of discrete, executable models, in which nodes are assigned discrete value ranges, and the relationship between them defined with simple mathematical operations. One tool for constructing such models is the BioModelAnalyzer (BMA), an open source and publicly available (https://www.biomodelanalyzer.org) software, aimed to be fully usable by researchers without extensive computational or mathematical experience. A fundamental question for executable models is whether the high level of abstraction substantially reduces expressivity relative to continuous approaches. Here, we present a canonical library of biological signalling motifs, initially defined by Tyson et al (2003), translated for the first time into the BMA. We show that; 1) these motifs are easily and fully translatable from continuous to discrete models, 2) Combining these motifs in a computationally naïve way generates a fully functional and predictive model of the yeast cell cycle.


2021 ◽  
Author(s):  
T.J. Sego ◽  
Josua O. Aponte-Serrano ◽  
Juliano F. Gianlupi ◽  
James A. Glazier

AbstractThe biophysics of an organism span scales from subcellular to organismal and include spatial processes like diffusion of molecules, cell migration, and flow of intravenous fluids. Mathematical biology seeks to explain biophysical processes in mathematical terms at, and across, all relevant spatial and temporal scales. While non-spatial, ordinary differential equation (ODE) models are often used and readily calibrated to experimental data, they do not explicitly represent spatial and stochastic features of a biological system, limiting their insights and applications. Spatial models describe biological systems with spatial information but are mathematically complex and computationally expensive, which limits the ability to calibrate and deploy them. In this work we develop a formal method for deriving cell-based, spatial, multicellular models from ODE models of population dynamics in biological systems, and vice-versa. We provide examples of generating spatiotemporal, multicellular models from ODE models of viral infection and immune response. In these models the determinants of agreement of spatial and non-spatial models are the degree of spatial heterogeneity in viral production and rates of extracellular viral diffusion and decay. We show how ODE model parameters can implicitly represent spatial parameters, and cell-based spatial models can generate uncertain predictions through sensitivity to stochastic cellular events, which is not a feature of ODE models. Using our method, we can test ODE models in a multicellular, spatial context and translate information to and from non-spatial and spatial models, which help to employ spatiotemporal multicellular models using calibrated ODE model parameters, investigate objects and processes implicitly represented by ODE model terms and parameters, and improve the reproducibility of spatial, stochastic models. We hope to employ our method to generate new ODE model terms from spatiotemporal, multicellular models, recast popular ODE models on a cellular basis, and generate better models for critical applications where spatial and stochastic features affect outcomes.Statement of SignificanceOrdinary differential equations (ODEs) are widely used to model and efficiently simulate multicellular systems without explicit spatial information, while spatial models permit explicit spatiotemporal modeling but are mathematically complicated and computationally expensive. In this work we develop a method to generate stochastic, agent-based, multiscale models of multicellular systems with spatial resolution at the cellular level according to non-spatial ODE models. We demonstrate how to directly translate model terms and parameters between ODE and spatial models and apply non-spatial model terms to boundary conditions using examples of viral infection modeling, and show how spatial models can interrogate implicitly represented biophysical mechanisms in non-spatial models. We discuss strategies for co-developing spatial and non-spatial models and reconciling disagreements between them.


2015 ◽  
Author(s):  
Louis Yang ◽  
Ming Yang

Sustained oscillations are frequently observed in biological systems consisting of a negative feedback loop, but a mathematical model with two ordinary differential equations (ODE) that has a negative feedback loop structure fails to produce sustained oscillations. Only when a time delay is introduced into the system by expanding to a three-ODE model, transforming to a two-DDE model, or introducing a bistable trigger do stable oscillations present themselves. In this study, we propose another mechanism for producing sustained oscillations based on periodic reaction pauses of chemical reactions in a negative feedback system. We model the oscillatory system behavior by allowing the coefficients in the two-ODE model to be periodic functions of time – called pulsate functions – to account for reactions with go-stop pulses. We find that replacing coefficients in the two-ODE system with pulsate functions with micro-scale (several seconds) pauses can produce stable system-wide oscillations that have periods of approximately one to several hours long. We also compare our two-ODE and three-ODE models with the two-DDE, three-ODE, and three-DDE models without the pulsate functions. Our numerical experiments suggest that sustained long oscillations in biological systems with a negative feedback loop may be an intrinsic property arising from the slow diffusion-based pulsate behavior of biochemical reactions.


BMC Biology ◽  
2021 ◽  
Vol 19 (1) ◽  
Author(s):  
T. J. Sego ◽  
Josua O. Aponte-Serrano ◽  
Juliano F. Gianlupi ◽  
James A. Glazier

Abstract Background The biophysics of an organism span multiple scales from subcellular to organismal and include processes characterized by spatial properties, such as the diffusion of molecules, cell migration, and flow of intravenous fluids. Mathematical biology seeks to explain biophysical processes in mathematical terms at, and across, all relevant spatial and temporal scales, through the generation of representative models. While non-spatial, ordinary differential equation (ODE) models are often used and readily calibrated to experimental data, they do not explicitly represent the spatial and stochastic features of a biological system, limiting their insights and applications. However, spatial models describing biological systems with spatial information are mathematically complex and computationally expensive, which limits the ability to calibrate and deploy them and highlights the need for simpler methods able to model the spatial features of biological systems. Results In this work, we develop a formal method for deriving cell-based, spatial, multicellular models from ODE models of population dynamics in biological systems, and vice versa. We provide examples of generating spatiotemporal, multicellular models from ODE models of viral infection and immune response. In these models, the determinants of agreement of spatial and non-spatial models are the degree of spatial heterogeneity in viral production and rates of extracellular viral diffusion and decay. We show how ODE model parameters can implicitly represent spatial parameters, and cell-based spatial models can generate uncertain predictions through sensitivity to stochastic cellular events, which is not a feature of ODE models. Using our method, we can test ODE models in a multicellular, spatial context and translate information to and from non-spatial and spatial models, which help to employ spatiotemporal multicellular models using calibrated ODE model parameters. We additionally investigate objects and processes implicitly represented by ODE model terms and parameters and improve the reproducibility of spatial, stochastic models. Conclusion We developed and demonstrate a method for generating spatiotemporal, multicellular models from non-spatial population dynamics models of multicellular systems. We envision employing our method to generate new ODE model terms from spatiotemporal and multicellular models, recast popular ODE models on a cellular basis, and generate better models for critical applications where spatial and stochastic features affect outcomes.


Author(s):  
Louis Yang ◽  
Ming Yang

Sustained oscillations are frequently observed in biological systems consisting of a negative feedback loop, but a mathematical model with two ordinary differential equations (ODE) that has a negative feedback loop structure fails to produce sustained oscillations. Only when a time delay is introduced into the system by expanding to a three-ODE model, transforming to a two-DDE model, or introducing a bistable trigger do stable oscillations present themselves. In this study, we propose another mechanism for producing sustained oscillations based on periodic reaction pauses of chemical reactions in a negative feedback system. We model the oscillatory system behavior by allowing the coefficients in the two-ODE model to be periodic functions of time – called pulsate functions – to account for reactions with go-stop pulses. We find that replacing coefficients in the two-ODE system with pulsate functions with micro-scale (several seconds) pauses can produce stable system-wide oscillations that have periods of approximately one to several hours long. We also compare our two-ODE and three-ODE models with the two-DDE, three-ODE, and three-DDE models without the pulsate functions. Our numerical experiments suggest that sustained long oscillations in biological systems with a negative feedback loop may be an intrinsic property arising from the slow diffusion-based pulsate behavior of biochemical reactions.


Author(s):  
Henry S. Slayter

Electron microscopic methods have been applied increasingly during the past fifteen years, to problems in structural molecular biology. Used in conjunction with physical chemical methods and/or Fourier methods of analysis, they constitute powerful tools for determining sizes, shapes and modes of aggregation of biopolymers with molecular weights greater than 50, 000. However, the application of the e.m. to the determination of very fine structure approaching the limit of instrumental resolving power in biological systems has not been productive, due to various difficulties such as the destructive effects of dehydration, damage to the specimen by the electron beam, and lack of adequate and specific contrast. One of the most satisfactory methods for contrasting individual macromolecules involves the deposition of heavy metal vapor upon the specimen. We have investigated this process, and present here what we believe to be the more important considerations for optimizing it. Results of the application of these methods to several biological systems including muscle proteins, fibrinogen, ribosomes and chromatin will be discussed.


Author(s):  
Nicholas J Severs

In his pioneering demonstration of the potential of freeze-etching in biological systems, Russell Steere assessed the future promise and limitations of the technique with remarkable foresight. Item 2 in his list of inherent difficulties as they then stood stated “The chemical nature of the objects seen in the replica cannot be determined”. This defined a major goal for practitioners of freeze-fracture which, for more than a decade, seemed unattainable. It was not until the introduction of the label-fracture-etch technique in the early 1970s that the mould was broken, and not until the following decade that the full scope of modern freeze-fracture cytochemistry took shape. The culmination of these developments in the 1990s now equips the researcher with a set of effective techniques for routine application in cell and membrane biology.Freeze-fracture cytochemical techniques are all designed to provide information on the chemical nature of structural components revealed by freeze-fracture, but differ in how this is achieved, in precisely what type of information is obtained, and in which types of specimen can be studied.


2019 ◽  
Vol 3 (5) ◽  
pp. 435-443 ◽  
Author(s):  
Addy Pross

Despite the considerable advances in molecular biology over the past several decades, the nature of the physical–chemical process by which inanimate matter become transformed into simplest life remains elusive. In this review, we describe recent advances in a relatively new area of chemistry, systems chemistry, which attempts to uncover the physical–chemical principles underlying that remarkable transformation. A significant development has been the discovery that within the space of chemical potentiality there exists a largely unexplored kinetic domain which could be termed dynamic kinetic chemistry. Our analysis suggests that all biological systems and associated sub-systems belong to this distinct domain, thereby facilitating the placement of biological systems within a coherent physical/chemical framework. That discovery offers new insights into the origin of life process, as well as opening the door toward the preparation of active materials able to self-heal, adapt to environmental changes, even communicate, mimicking what transpires routinely in the biological world. The road to simplest proto-life appears to be opening up.


Sign in / Sign up

Export Citation Format

Share Document