Publications by Colleges and Departments (MSU - Bozeman)
Permanent URI for this communityhttps://scholarworks.montana.edu/handle/1/3
Browse
38 results
Search Results
Item Ribosome Abundance Control in Prokaryotes(Springer Science and Business Media LLC, 2023-10) Shea, Jacob; Davis, Lisa; Quaye, Bright; Gedeon, TomasCell growth is an essential phenotype of any unicellular organism and it crucially depends on precise control of protein synthesis. We construct a model of the feedback mechanisms that regulate abundance of ribosomes in E. coli, a prototypical prokaryotic organism. Since ribosomes are needed to produce more ribosomes, the model includes a positive feedback loop central to the control of cell growth. Our analysis of the model shows that there can be only two coexisting equilibrium states across all 23 parameters. This precludes the existence of hysteresis, suggesting that the ribosome abundance changes continuously with parameters. These states are related by a transcritical bifurcation, and we provide an analytic formula for parameters that admit either state.Item Resource allocation accounts for the large variability of rate-yield phenotypes across bacterial strains(eLife Sciences Publications, Ltd, 2023-05) Baldazzi, Valentina; Ropers, Delphine; Gouzé, Jean-Luc; Gedeon, Tomas; de Jong, HiddeDifferent strains of a microorganism growing in the same environment display a wide variety of growth rates and growth yields. We developed a coarse-grained model to test the hypothesis that different resource allocation strategies, corresponding to different compositions of the proteome, can account for the observed rate-yield variability. The model predictions were verified by means of a database of hundreds of published rate-yield and uptake-secretion phenotypes of Escherichia coli strains grown in standard laboratory conditions. We found a very good quantitative agreement between the range of predicted and observed growth rates, growth yields, and glucose uptake and acetate secretion rates. These results support the hypothesis that resource allocation is a major explanatory factor of the observed variability of growth rates and growth yields across different bacterial strains. An interesting prediction of our model, supported by the experimental data, is that high growth rates are not necessarily accompanied by low growth yields. The resource allocation strategies enabling high-rate, high-yield growth of E. coli lead to a higher saturation of enzymes and ribosomes, and thus to a more efficient utilization of proteomic resources. Our model thus contributes to a fundamental understanding of the quantitative relationship between rate and yield in E. coli and other microorganisms. It may also be useful for the rapid screening of strains in metabolic engineering and synthetic biology.Item Extending Combinatorial Regulatory Network Modeling to Include Activity Control and Decay Modulation(Society for Industrial & Applied Mathematics, 2022-09) Cummins, Bree; Gameiro, Marcio; Gedeon, Tomas; Kepley, Shane; Mischaikow, Konstantin; Zhang, LunUnderstanding how the structure of within-system interactions affects the dynamics of the system is important in many areas of science. We extend a network dynamics modeling platform DSGRN, which combinatorializes both dynamics and parameter space to construct finite but accurate summaries of network dynamics, to new types of interactions. While the standard DSGRN assumes that each network edge controls the rate of abundance of the target node, the new edges may control either activity level or a decay rate of its target. While motivated by processes of post-transcriptional modification and ubiquitination in systems biology, our extension is applicable to the dynamics of any signed directed network.Item Stability and Bifurcations of Equilibria in Networks with Piecewise Linear Interactions(World Scientific Pub Co Pte Lt, 2021-09) Duncan, William; Gedeon, TomasIn this paper, we study equilibria of differential equation models for networks. When interactions between nodes are taken to be piecewise constant, an efficient combinatorial analysis can be used to characterize the equilibria. When the piecewise constant functions are replaced with piecewise linear functions, the equilibria are preserved as long as the piecewise linear functions are sufficiently steep. Therefore the combinatorial analysis can be leveraged to understand a broader class of interactions. To better understand how broad this class is, we explicitly characterize how steep the piecewise linear functions must be for the correspondence between equilibria to hold. To do so, we analyze the steady state and Hopf bifurcations which cause a change in the number or stability of equilibria as slopes are decreased. Additionally, we show how to choose a subset of parameters so that the correspondence between equilibria holds for the smallest possible slopes when the remaining parameters are fixed.Item Genetic networks encode secrets of their past(Elsevier BV, 2022-03) Crawford-Kahrl, Peter; Nerem, Robert R.; Cummins, Bree; Gedeon, TomasResearch shows that gene duplication followed by either repurposing or removal of duplicated genes is an important contributor to evolution of gene and protein interaction networks. We aim to identify which characteristics of a network can arise through this process, and which must have been produced in a different way. To model the network evolution, we postulate vertex duplication and edge deletion as evolutionary operations on graphs. Using the novel concept of an ancestrally distinguished subgraph, we show how features of present-day networks require certain features of their ancestors. In particular, ancestrally distinguished subgraphs cannot be introduced by vertex duplication. Additionally, if vertex duplication and edge deletion are the only evolutionary mechanisms, then a graph’s ancestrally distinguished subgraphs must be contained in all of the graph’s ancestors. We analyze two experimentally derived genetic networks and show that our results accurately predict lack of large ancestrally distinguished subgraphs, despite this feature being statistically improbable in associated random networks. This observation is consistent with the hypothesis that these networks evolved primarily via vertex duplication. The tools we provide open the door for analyzing ancestral networks using current networks. Our results apply to edge-labeled (e.g. signed) graphs which are either undirected or directed.Item Using extremal events to characterize noisy time series(2020-02) Berry, Eric; Cummins, Bree; Nerem, Robert R.; Smith, Lauren M.; Haase, Steven B.; Gedeon, TomasExperimental time series provide an informative window into the underlying dynamical system, and the timing of the extrema of a time series (or its derivative) contains information about its structure. However, the time series often contain significant measurement errors. We describe a method for characterizing a time series for any assumed level of measurement error 𝜀 by a sequence of intervals, each of which is guaranteed to contain an extremum for any function that 𝜀-approximates the time series. Based on the merge tree of a continuous function, we define a new object called the normalized branch decomposition, which allows us to compute intervals for any level 𝜀. We show that there is a well-defined total order on these intervals for a single time series, and that it is naturally extended to a partial order across a collection of time series comprising a dataset. We use the order of the extracted intervals in two applications. First, the partial order describing a single dataset can be used to pattern match against switching model output (Cummins et al. in SIAM J Appl Dyn Syst 17(2):1589–1616, 2018), which allows the rejection of a network model. Second, the comparison between graph distances of the partial orders of different datasets can be used to quantify similarity between biological replicates.Item Multi-parameter exploration of dynamics of regulatory networks(2020-04) Gedeon, TomasOver the last twenty years advances in systems biology have changed our views on microbial communities and promise to revolutionize treatment of human diseases. In almost all scientific breakthroughs since time of Newton, mathematical modeling has played a prominent role. Regulatory networks emerged as preferred descriptors of how abundances of molecular species depend on each other. However, the central question on how cellular phenotypes emerge from dynamics of these network remains elusive. The principal reason is that differential equation models in the field of biology (while so successful in areas of physics and physical chemistry), do not arise from first principles, and these models suffer from lack of proper parameterization. In response to these challenges, discrete time models based on Boolean networks have been developed. In this review, we discuss an emerging modeling paradigm that combines ideas from differential equations and Boolean models, and has been developed independently within dynamical systems and computer science communities. The result is an approach that can associate a range of potential dynamical behaviors to a network, arrange the descriptors of the dynamics in a searchable database, and allows for multi-parameter exploration of the dynamics akin to bifurcation theory. Since this approach is computationally accessible for moderately sized networks, it allows, perhaps for the first time, to rationally compare different network topologies based on their dynamics.Item Comparison of Combinatorial Signatures of Global Network Dynamics Generated by Two Classes of ODE Models(2019-04) Crawford-Kahrl, Peter; Cummins, Bree; Gedeon, TomasModeling the dynamics of biological networks introduces many challenges, among them the lack of first principle models, the size of the networks, and difficulties with parameterization. Discrete time Boolean networks and related continuous time switching systems provide a computationally accessible way to translate the structure of the network to predictions about the dynamics. Recent work has shown that the parameterized dynamics of switching systems can be captured by a combinatorial object, called a Dynamic Signatures Generated by Regulatory Networks (DSGRN) database, that consists of a parameter graph characterizing a finite parameter space decomposition, whose nodes are assigned a Morse graph that captures global dynamics for all corresponding parameters. We show that for a given network there is a way to associate the same type of object by considering a continuous time ODE system with a continuous right-hand side, which we call an L-system. The main goal of this paper is to compare the two DSGRN databases for the same network. Since the L-systems can be thought of as perturbations (not necessarily small) of the switching systems, our results address the correspondence between global parameterized dynamics of switching systems and their perturbations. We show that, at corresponding parameters, there is an order preserving map from the Morse graph of the switching system to that of the L-system that is surjective on the set of attractors and bijective on the set of fixed-point attractors. We provide important examples showing why this correspondence cannot be strengthened.Item DSGRN: Examining the Dynamics of Families of Logical Models(2018-06) Cummins, Bree; Gedeon, Tomas; Harker, Shaun; Mischaikow, KonstantinWe present a computational tool DSGRN for exploring the dynamics of a network by computing summaries of the dynamics of switching models compatible with the network across all parameters. The network can arise directly from a biological problem, or indirectly as the interaction graph of a Boolean model. This tool computes a finite decomposition of parameter space such that for each region, the state transition graph that describes the coarse dynamical behavior of a network is the same. Each of these parameter regions corresponds to a different logical description of the network dynamics. The comparison of dynamics across parameters with experimental data allows the rejection of parameter regimes or entire networks as viable models for representing the underlying regulatory mechanisms. This in turn allows a search through the space of perturbations of a given network for networks that robustly fit the data. These are the first steps toward discovering a network that optimally matches the observed dynamics by searching through the space of networks.Item Global dynamics for switching systems and their extensions by linear differential equations(2018-11) Huttinga, Zane; Cummins, Bree; Gedeon, Tomas; Mischaikow, KonstantinSwitching systems use piecewise constant nonlinearities to model gene regulatory networks. This choice provides advantages in the analysis of behavior and allows the global description of dynamics in terms of Morse graphs associated to nodes of a parameter graph. The parameter graph captures spatial characteristics of a decomposition of parameter space into domains with identical Morse graphs. However, there are many cellular processes that do not exhibit threshold-like behavior and thus are not well described by a switching system. We consider a class of extensions of switching systems formed by a mixture of switching interactions and chains of variables governed by linear differential equations. We show that the parameter graphs associated to the switching system and any of its extensions are identical. For each parameter graph node, there is an order-preserving map from the Morse graph of the switching system to the Morse graph of any of its extensions. We provide counterexamples that show why possible stronger relationships between the Morse graphs are not valid.