1
|
Olsen LF, Lunding A. On the coupling of intracellular K + ${{\rm{K}}}^{+}$ to glycolytic oscillations in yeast. Yeast 2024; 41:486-498. [PMID: 39031655 DOI: 10.1002/yea.3972] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2024] [Revised: 05/06/2024] [Accepted: 06/05/2024] [Indexed: 07/22/2024] Open
Abstract
We have investigated the interplay between glycolytic oscillations and intracellularK + ${{\rm{K}}}^{+}$ concentration in the yeast Saccharomyces cerevisiae. IntracellularK + ${{\rm{K}}}^{+}$ concentration was measured using the fluorophore potassium-binding benzofuranisophthalate (PBFI). We found thatK + ${{\rm{K}}}^{+}$ is an essential ion for the occurrence of glycolytic oscillations and that intracellularK + ${{\rm{K}}}^{+}$ concentration oscillates synchronously with other variables such as nicotinamide adenine dinucleotide hydride (NADH), intracellular adenosine triphosphate (ATP), and mitochondrial membrane potential. We also investigated if glycolysis and intracellularK + ${{\rm{K}}}^{+}$ concentration oscillate in a number of yeast strains with mutations inK + ${{\rm{K}}}^{+}$ transporters in the plasma membrane, mitochondrial membrane and in the vacuolar membrane. Most of these strains are still capable of showing glycolytic oscillations, but two strains are not: (i) a strain with a deletion in the mitochondrial Mdm38pK + ∕ H + ${{\rm{K}}}^{+}\unicode{x02215}{{\rm{H}}}^{+}$ transporter and (ii) a strain with deletion of the late endosomal Nhx1pK + ∕ H + ${{\rm{K}}}^{+}\unicode{x02215}{{\rm{H}}}^{+}$ (Na + ∕ H + ${\text{Na}}^{+}\unicode{x02215}{{\rm{H}}}^{+}$ ) transporter. In these two mutant strains intracellularK + ${{\rm{K}}}^{+}$ concentration seems to be low, indicating that the two transporters may be involved in transport ofK + ${{\rm{K}}}^{+}$ into the cytosol. In the strain, Mdm38pΔ ${\rm{\Delta }}$ oscillations in glycolysis could be restored by addition of theK + ∕ H + ${{\rm{K}}}^{+}\unicode{x02215}{{\rm{H}}}^{+}$ exchange ionophore nigericin. Furthermore, in two nonoscillating mutant strains with a defective V-ATPase and deletion of the Arp1p protein the intracellularK + ${{\rm{K}}}^{+}$ is relatively high, suggesting that the V-ATPase is essential for transport ofK + ${{\rm{K}}}^{+}$ out of the cytosol and that the cytoskeleton may be involved in bindingK + ${{\rm{K}}}^{+}$ to reduce the concentration of free ion in the cytosol. Analyses of the time series of oscillations of NADH, ATP, mitochondrial membrane potential, and potassium concentration using data-driven modeling corroborate the conjecture thatK + ${{\rm{K}}}^{+}$ ion is essential for the emergence of oscillations and support the experimental findings using mutant strains.
Collapse
Affiliation(s)
- Lars F Olsen
- PhyLife, Institute of Biochemistry and Molecular Biology, University of Southern Denmark, Odense M, Denmark
| | - Anita Lunding
- PhyLife, Institute of Biochemistry and Molecular Biology, University of Southern Denmark, Odense M, Denmark
| |
Collapse
|
2
|
Fung SYS, Xǔ XJ, Wu M. Nonlinear dynamics in phosphoinositide metabolism. Curr Opin Cell Biol 2024; 88:102373. [PMID: 38797149 PMCID: PMC11186694 DOI: 10.1016/j.ceb.2024.102373] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2024] [Revised: 04/26/2024] [Accepted: 05/01/2024] [Indexed: 05/29/2024]
Abstract
Phosphoinositides broadly impact membrane dynamics, signal transduction and cellular physiology. The orchestration of signaling complexity by this seemingly simple metabolic pathway remains an open question. It is increasingly evident that comprehending the complexity of the phosphoinositides metabolic network requires a systems view based on nonlinear dynamics, where the products of metabolism can either positively or negatively modulate enzymatic function. These feedback and feedforward loops may be paradoxical, leading to counterintuitive effects. In this review, we introduce the framework of nonlinear dynamics, emphasizing distinct dynamical regimes such as the excitable state, oscillations, and mixed-mode oscillations-all of which have been experimentally observed in phosphoinositide metabolisms. We delve into how these dynamical behaviors arise from one or multiple network motifs, including positive and negative feedback loops, coherent and incoherent feedforward loops. We explore the current understanding of the molecular circuits responsible for these behaviors. While mapping these circuits presents both conceptual and experimental challenges, redefining cellular behavior based on dynamical state, lipid fluxes, time delay, and network topology is likely essential for a comprehensive understanding of this fundamental metabolic network.
Collapse
Affiliation(s)
- Suet Yin Sarah Fung
- Department of Cell Biology, Yale University School of Medicine, 333 Cedar Street, New Haven, CT, 06520-8002, USA
| | - X J Xǔ
- Department of Physics, Yale University, New Haven, CT, 06511, USA
| | - Min Wu
- Department of Cell Biology, Yale University School of Medicine, 333 Cedar Street, New Haven, CT, 06520-8002, USA.
| |
Collapse
|
3
|
Del Olmo M, Legewie S, Brunner M, Höfer T, Kramer A, Blüthgen N, Herzel H. Network switches and their role in circadian clocks. J Biol Chem 2024; 300:107220. [PMID: 38522517 PMCID: PMC11044057 DOI: 10.1016/j.jbc.2024.107220] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2023] [Revised: 03/07/2024] [Accepted: 03/18/2024] [Indexed: 03/26/2024] Open
Abstract
Circadian rhythms are generated by complex interactions among genes and proteins. Self-sustained ∼24 h oscillations require negative feedback loops and sufficiently strong nonlinearities that are the product of molecular and network switches. Here, we review common mechanisms to obtain switch-like behavior, including cooperativity, antagonistic enzymes, multisite phosphorylation, positive feedback, and sequestration. We discuss how network switches play a crucial role as essential components in cellular circadian clocks, serving as integral parts of transcription-translation feedback loops that form the basis of circadian rhythm generation. The design principles of network switches and circadian clocks are illustrated by representative mathematical models that include bistable systems and negative feedback loops combined with Hill functions. This work underscores the importance of negative feedback loops and network switches as essential design principles for biological oscillations, emphasizing how an understanding of theoretical concepts can provide insights into the mechanisms generating biological rhythms.
Collapse
Affiliation(s)
- Marta Del Olmo
- Institute for Theoretical Biology, Humboldt Universität zu Berlin and Charité Universitätsmedizin Berlin, Berlin, Germany.
| | - Stefan Legewie
- Department of Systems Biology, Institute for Biomedical Genetics (IBMG), University of Stuttgart, Stuttgart, Germany; Stuttgart Research Center for Systems Biology (SRCSB), University of Stuttgart, Stuttgart, Germany
| | - Michael Brunner
- Biochemistry Center, Universität Heidelberg, Heidelberg, Germany
| | - Thomas Höfer
- Division of Theoretical Systems Biology, German Cancer Research Center (DKFZ), Universität Heidelberg, Heidelberg, Germany
| | - Achim Kramer
- Laboratory of Chronobiology, Institute for Medical Immunology, Charité Universitätsmedizin Berlin, Berlin, Germany
| | - Nils Blüthgen
- Institute for Theoretical Biology, Humboldt Universität zu Berlin and Charité Universitätsmedizin Berlin, Berlin, Germany; Institute of Pathology, Charité Universitätsmedizin Berlin, Berlin, Germany
| | - Hanspeter Herzel
- Institute for Theoretical Biology, Humboldt Universität zu Berlin and Charité Universitätsmedizin Berlin, Berlin, Germany.
| |
Collapse
|
4
|
van Niekerk DD, van Wyk M, Kouril T, Snoep JL. Kinetic modelling of glycolytic oscillations. Essays Biochem 2024; 68:15-25. [PMID: 38206647 DOI: 10.1042/ebc20230037] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2023] [Revised: 12/18/2023] [Accepted: 12/19/2023] [Indexed: 01/12/2024]
Abstract
Glycolytic oscillations have been studied for well over 60 years, but aspects of their function, and mechanisms of regulation and synchronisation remain unclear. Glycolysis is amenable to mechanistic mathematical modelling, as its components have been well characterised, and the system can be studied at many organisational levels: in vitro reconstituted enzymes, cell free extracts, individual cells, and cell populations. In recent years, the emergence of individual cell analysis has opened new ways of studying this intriguing system.
Collapse
Affiliation(s)
- David D van Niekerk
- Department of Biochemistry, Stellenbosch University, Matieland, South Africa
| | - Morne van Wyk
- Department of Biochemistry, Stellenbosch University, Matieland, South Africa
| | - Theresa Kouril
- Department of Biochemistry, Stellenbosch University, Matieland, South Africa
| | - Jacky L Snoep
- Department of Biochemistry, Stellenbosch University, Matieland, South Africa
- Molecular Cell Biology, Vrije Universiteit, Amsterdam, The Netherlands
| |
Collapse
|
5
|
Prokop B, Gelens L. From biological data to oscillator models using SINDy. iScience 2024; 27:109316. [PMID: 38523784 PMCID: PMC10959654 DOI: 10.1016/j.isci.2024.109316] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2023] [Revised: 01/18/2024] [Accepted: 02/18/2024] [Indexed: 03/26/2024] Open
Abstract
Periodic changes in the concentration or activity of different molecules regulate vital cellular processes such as cell division and circadian rhythms. Developing mathematical models is essential to better understand the mechanisms underlying these oscillations. Recent data-driven methods like SINDy have fundamentally changed model identification, yet their application to experimental biological data remains limited. This study investigates SINDy's constraints by directly applying it to biological oscillatory data. We identify insufficient resolution, noise, dimensionality, and limited prior knowledge as primary limitations. Using various generic oscillator models of different complexity and/or dimensionality, we systematically analyze these factors. We then propose a comprehensive guide for inferring models from biological data, addressing these challenges step by step. Our approach is validated using glycolytic oscillation data from yeast.
Collapse
Affiliation(s)
- Bartosz Prokop
- Laboratory of Dynamics in Biological Systems, Department of Cellular and Molecular Medicine, KU Leuven, Herestraat 49, 3000 Leuven, Belgium
| | - Lendert Gelens
- Laboratory of Dynamics in Biological Systems, Department of Cellular and Molecular Medicine, KU Leuven, Herestraat 49, 3000 Leuven, Belgium
| |
Collapse
|
6
|
Khatri N, Kapral R. Clustering of chemically propelled nanomotors in chemically active environments. CHAOS (WOODBURY, N.Y.) 2024; 34:033103. [PMID: 38427933 DOI: 10.1063/5.0188624] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/22/2023] [Accepted: 02/05/2024] [Indexed: 03/03/2024]
Abstract
Synthetic nanomotors powered by chemical reactions have been designed to act as vehicles for active cargo transport, drug delivery, and a variety of other uses. Collections of such motors, acting in consort, can self-assemble to form swarms or clusters, providing opportunities for applications on various length scales. While such collective behavior has been studied when the motors move in a chemically inactive fluid environment, when the medium in which they move is a chemical network that supports complex spatial and temporal patterns, through simulation and theoretical analysis we show that collective behavior changes. Spatial patterns in the environment can guide and control motor collective states, and interactions of the motors with their environment can give rise to distinctive spatiotemporal motor patterns. The results are illustrated by studies of the motor dynamics in systems that support Turing patterns and spiral waves. This work is relevant for potential applications that involve many active nanomotors moving in complex chemical or biological environments.
Collapse
Affiliation(s)
- Narender Khatri
- Chemical Physics Theory Group, Department of Chemistry, University of Toronto, Toronto, Ontario M5S 3H6, Canada
| | - Raymond Kapral
- Chemical Physics Theory Group, Department of Chemistry, University of Toronto, Toronto, Ontario M5S 3H6, Canada
| |
Collapse
|
7
|
Sarkar P, Ray DS. Tuning limit cycles with a noise: Survival and collapse. Phys Rev E 2024; 109:034209. [PMID: 38632777 DOI: 10.1103/physreve.109.034209] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/16/2023] [Accepted: 02/15/2024] [Indexed: 04/19/2024]
Abstract
We consider a general class of limit cycle oscillators driven by an additive Gaussian white noise. Based on the separation of timescales, we construct the equation of motion for slow dynamics after appropriate averaging over the fast motion. The equation for slow motion whose coefficients are modified by noise characteristics is solved to obtain the analytic solution in the long time limit. We show that with increase of noise strength, the loop area of the limit cycle decreases until a critical value is reached, beyond which the limit cycle collapses. We determine the noise threshold from the condition for removal of secular divergence of the perturbation series and work out two explicit examples of Van der Pol and Duffing-Van der Pol oscillators for corroboration between the theory and numerics.
Collapse
Affiliation(s)
- Prasun Sarkar
- Indian Association for the Cultivation of Science, Jadavpur, Kolkata-700032, India
| | - Deb Shankar Ray
- Indian Association for the Cultivation of Science, Jadavpur, Kolkata-700032, India
| |
Collapse
|
8
|
Peng Z, Kapral R. Self-organization of active colloids mediated by chemical interactions. SOFT MATTER 2024; 20:1100-1113. [PMID: 38221884 DOI: 10.1039/d3sm01272g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/16/2024]
Abstract
Self-propelled colloidal particles exhibit rich non-equilibrium phenomena and have promising applications in fields such as drug delivery and self-assembled active materials. Previous experimental and theoretical studies have shown that chemically active colloids that consume or produce a chemical can self-organize into clusters with diverse characteristics depending on the effective phoretic interactions. In this paper, we investigate self-organization in systems with multiple chemical species that undergo a network of reactions and multiple colloidal species that participate in different reactions. Active colloids propelled by complex chemical reactions with potentially nonlinear kinetics can be realized using enzymatic reactions that occur on the surface of enzyme-coated particles. To demonstrate how the self-organizing behavior depends on the chemical reactions active colloids catalyze and their chemical environment, we consider first a single type of colloid undergoing a simple catalytic reaction, and compare this often-studied case with self-organization in binary mixtures of colloids with sequential reactions, and binary mixtures with nonlinear autocatalytic reactions. Our results show that in general active colloids at low particle densities can form localized clusters in the presence of bulk chemical reactions and phoretic attractions. The characteristics of the clusters, however, depend on the reaction kinetics in the bulk and on the particles and phoretic coefficients. With one or two chemical species that only undergo surface reactions, the space for possible self-organizations are limited. By considering the additional system parameters that enter the chemical reaction network involving reactions on the colloids and in the fluid, the design space of colloidal self-organization can be enlarged, leading to a variety of non-equilibrium structures.
Collapse
Affiliation(s)
- Zhiwei Peng
- Chemical Physics Theory Group, Department of Chemistry, University of Toronto, Toronto, ON, M5S 3H6, Canada.
| | - Raymond Kapral
- Chemical Physics Theory Group, Department of Chemistry, University of Toronto, Toronto, ON, M5S 3H6, Canada.
| |
Collapse
|
9
|
Banaji M, Boros B, Hofbauer J. Oscillations in three-reaction quadratic mass-action systems. STUDIES IN APPLIED MATHEMATICS (CAMBRIDGE, MASS.) 2024; 152:249-278. [PMID: 38516583 PMCID: PMC10953466 DOI: 10.1111/sapm.12639] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 04/08/2023] [Revised: 06/27/2023] [Accepted: 08/22/2023] [Indexed: 03/23/2024]
Abstract
It is known that rank-two bimolecular mass-action systems do not admit limit cycles. With a view to understanding which small mass-action systems admit oscillation, in this paper we study rank-two networks with bimolecular source complexes but allow target complexes with higher molecularities. As our goal is to find oscillatory networks of minimal size, we focus on networks with three reactions, the minimum number that is required for oscillation. However, some of our intermediate results are valid in greater generality. One key finding is that an isolated periodic orbit cannot occur in a three-reaction, trimolecular, mass-action system with bimolecular sources. In fact, we characterize all networks in this class that admit a periodic orbit; in every case, all nearby orbits are periodic too. Apart from the well-known Lotka and Ivanova reactions, we identify another network in this class that admits a center. This new network exhibits a vertical Andronov-Hopf bifurcation. Furthermore, we characterize all two-species, three-reaction, bimolecular-sourced networks that admit an Andronov-Hopf bifurcation with mass-action kinetics. These include two families of networks that admit a supercritical Andronov-Hopf bifurcation and hence a stable limit cycle. These networks necessarily have a target complex with a molecularity of at least four, and it turns out that there are exactly four such networks that are tetramolecular.
Collapse
Affiliation(s)
- Murad Banaji
- Mathematical InstituteUniversity of OxfordOxfordUK
| | - Balázs Boros
- Department of MathematicsUniversity of ViennaViennaAustria
| | - Josef Hofbauer
- Department of MathematicsUniversity of ViennaViennaAustria
| |
Collapse
|
10
|
Jin J, Xu F, Liu Z, Qi H, Yao C, Shuai J, Li X. Biphasic amplitude oscillator characterized by distinct dynamics of trough and crest. Phys Rev E 2023; 108:064412. [PMID: 38243441 DOI: 10.1103/physreve.108.064412] [Citation(s) in RCA: 17] [Impact Index Per Article: 17.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2023] [Accepted: 11/28/2023] [Indexed: 01/21/2024]
Abstract
Biphasic amplitude dynamics (BAD) of oscillation have been observed in many biological systems. However, the specific topology structure and regulatory mechanisms underlying these biphasic amplitude dynamics remain elusive. Here, we searched all possible two-node circuit topologies and identified the core oscillator that enables robust oscillation. This core oscillator consists of a negative feedback loop between two nodes and a self-positive feedback loop of the input node, which result in the fast and slow dynamics of the two nodes, thereby achieving relaxation oscillation. Landscape theory was employed to study the stochastic dynamics and global stability of the system, allowing us to quantitatively describe the diverse positions and sizes of the Mexican hat. With increasing input strength, the size of the Mexican hat exhibits a gradual increase followed by a subsequent decrease. The self-activation of input node and the negative feedback on input node, which dominate the fast dynamics of the input node, were observed to regulate BAD in a bell-shaped manner. Both deterministic and statistical analysis results reveal that BAD is characterized by the linear and nonlinear dependence of the oscillation trough and crest on the input strength. In addition, combining with computational and theoretical analysis, we addressed that the linear response of trough to input is predominantly governed by the negative feedback, while the nonlinear response of crest is jointly regulated by the negative feedback loop and the self-positive feedback loop within the oscillator. Overall, this study provides a natural and physical basis for comprehending the occurrence of BAD in oscillatory systems, yielding guidance for the design of BAD in synthetic biology applications.
Collapse
Affiliation(s)
- Jun Jin
- Department of Physics, Xiamen University, Xiamen, Fujian 361005, China
| | - Fei Xu
- Department of Physics, Anhui Normal University, Wuhu, Anhui 241002, China
| | - Zhilong Liu
- Department of Physics, Xiamen University, Xiamen, Fujian 361005, China
| | - Hong Qi
- Complex Systems Research Center, Shanxi University, Taiyuan, Shanxi 030006, China
| | - Chenggui Yao
- College of Data Science, Jiaxing University, Jiaxing, Zhejiang 314000, China
| | - Jianwei Shuai
- Department of Physics, Xiamen University, Xiamen, Fujian 361005, China
- Oujiang Laboratory (Zhejiang Lab for Regenerative Medicine, Vision and Brain Health) and Wenzhou Institute, University of Chinese Academy of Sciences, Wenzhou, Zhejiang 325001, China
| | - Xiang Li
- Department of Physics, Xiamen University, Xiamen, Fujian 361005, China
| |
Collapse
|
11
|
Anderson CNK, Transtrum MK. Sloppy model analysis identifies bifurcation parameters without normal form analysis. Phys Rev E 2023; 108:064215. [PMID: 38243539 DOI: 10.1103/physreve.108.064215] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/23/2023] [Accepted: 11/30/2023] [Indexed: 01/21/2024]
Abstract
Bifurcation phenomena are common in multidimensional multiparameter dynamical systems. Normal form theory suggests that bifurcations are driven by relatively few combinations of parameters. Models of complex systems, however, rarely appear in normal form, and bifurcations are controlled by nonlinear combinations of the bare parameters of differential equations. Discovering reparameterizations to transform complex equations into a normal form is often very difficult, and the reparameterization may not even exist in a closed form. Here we show that information geometry and sloppy model analysis using the Fisher information matrix can be used to identify the combination of parameters that control bifurcations. By considering observations on increasingly long timescales, we find those parameters that rapidly characterize the system's topological inhomogeneities, whether the system is in normal form or not. We anticipate that this novel analytical method, which we call time-widening information geometry (TWIG), will be useful in applied network analysis.
Collapse
Affiliation(s)
| | - Mark K Transtrum
- Department of Physics and Astronomy Brigham Young University and Provo, Utah 84604, USA
| |
Collapse
|
12
|
Huffman A, Shum H. Boundary-bound reactions: Pattern formation with and without hydrodynamics. Phys Rev E 2023; 108:055103. [PMID: 38115506 DOI: 10.1103/physreve.108.055103] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/02/2023] [Accepted: 10/13/2023] [Indexed: 12/21/2023]
Abstract
We study chemical pattern formation in a fluid between two flat plates and the effect of such patterns on the formation of convective cells. This patterning is made possible by assuming the plates are chemically reactive or release reagents into the fluid, both of which we model as chemical fluxes. We consider this as a specific example of boundary-bound reactions. In the absence of coupling with fluid flow, we show that the two-reagent system with nonlinear reactions admits chemical instabilities equivalent to diffusion-driven Turing instabilities. In the other extreme, when chemical fluxes at the two bounding plates are constant, diffusion-driven instabilities do not occur but hydrodynamic phenomena analogous to Rayleigh-Bénard convection are possible. Assuming we can influence the chemical fluxes along the domain and select suitable reaction systems, this presents a mechanism for the control of chemical and hydrodynamic instabilities and pattern formation. We study a generic class of models and find necessary conditions for a bifurcation to pattern formation. Afterwards, we present two examples derived from the Schnakenberg-Selkov reaction. Unlike the classical Rayleigh-Bénard instability, which requires a sufficiently large unstable density gradient, a chemohydrodynamic instability based on Turing-style pattern formation can emerge from a state that is uniform in density. We also find parameter combinations that result in the formation of convective cells whether gravity acts upwards or downwards relative to the reactive plate. The wave number of the cells and the direction of the flow at regions of high/low concentration depend on the orientation, hence, different patterns can be elicited by simply inverting the device. More generally, our results suggest methods for controlling pattern formation and convection by tuning reaction parameters. As a consequence, we can drive and alter fluid flow in a chamber without mechanical pumps by influencing the chemical instabilities.
Collapse
Affiliation(s)
- Aiden Huffman
- Department of Applied Mathematics, University of Waterloo, 200 University Ave W, Waterloo, ON N2L 3G1, Canada
| | - Henry Shum
- Department of Applied Mathematics, University of Waterloo, 200 University Ave W, Waterloo, ON N2L 3G1, Canada
| |
Collapse
|
13
|
Lauber N, Tichacek O, Narayanankutty K, De Martino D, Ruiz-Mirazo K. Collective catalysis under spatial constraints: Phase separation and size-scaling effects on mass action kinetics. Phys Rev E 2023; 108:044410. [PMID: 37978605 DOI: 10.1103/physreve.108.044410] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2023] [Accepted: 09/07/2023] [Indexed: 11/19/2023]
Abstract
Chemical reactions are usually studied under the assumption that both substrates and catalysts are well-mixed (WM) throughout the system. Although this is often applicable to test-tube experimental conditions, it is not realistic in cellular environments, where biomolecules can undergo liquid-liquid phase separation (LLPS) and form condensates, leading to important functional outcomes, including the modulation of catalytic action. Similar processes may also play a role in protocellular systems, like primitive coacervates, or in membrane-assisted prebiotic pathways. Here we explore whether the demixing of catalysts could lead to the formation of microenvironments that influence the kinetics of a linear (multistep) reaction pathway, as compared to a WM system. We implemented a general lattice model to simulate LLPS of a collection of different catalysts and extended it to include diffusion and a sequence of reactions of small substrates. We carried out a quantitative analysis of how the phase separation of the catalysts affects reaction times depending on the affinity between substrates and catalysts, the length of the reaction pathway, the system size, and the degree of homogeneity of the condensate. A key aspect underlying the differences reported between the two scenarios is that the scale invariance observed in the WM system is broken by condensation processes. The main theoretical implications of our results for mean-field chemistry are drawn, extending the mass action kinetics scheme to include substrate initial "hitting times" to reach the catalysts condensate. We finally test this approach by considering open nonlinear conditions, where we successfully predict, through microscopic simulations, that phase separation inhibits chemical oscillatory behavior, providing a possible explanation for the marginal role that this complex dynamic behavior plays in real metabolisms.
Collapse
Affiliation(s)
- Nino Lauber
- Biofisika Institute (CSIC, UPV/EHU), 48940 Leioa, Spain
- Donostia International Physics Center, 20018 Donostia-San Sebastian, Spain
- Department of Philosophy, University of the Basque Country, 20018 Donostia-San Sebastian, Spain
| | - Ondrej Tichacek
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, 160 00 Prague, Czech Republic
| | - Krishnadev Narayanankutty
- Biofisika Institute (CSIC, UPV/EHU), 48940 Leioa, Spain
- Department of Molecular Biology and Biochemistry, University of the Basque Country, 48940 Leioa, Spain
| | - Daniele De Martino
- Biofisika Institute (CSIC, UPV/EHU), 48940 Leioa, Spain
- Ikerbasque Foundation, 48009 Bilbao, Spain
| | - Kepa Ruiz-Mirazo
- Biofisika Institute (CSIC, UPV/EHU), 48940 Leioa, Spain
- Department of Philosophy, University of the Basque Country, 20018 Donostia-San Sebastian, Spain
| |
Collapse
|
14
|
Gagnon JS, Hochberg D. Conditions for the origin of homochirality in primordial catalytic reaction networks. Sci Rep 2023; 13:9885. [PMID: 37336897 DOI: 10.1038/s41598-023-36852-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/28/2023] [Accepted: 06/11/2023] [Indexed: 06/21/2023] Open
Abstract
We study the generation of homochirality in a general chemical model (based on the homogeneous, fully connected Smoluchowski aggregation-fragmentation model) that obeys thermodynamics and can be easily mapped onto known origin of life models (e.g. autocatalytic sets, hypercycles, etc.), with essential aspects of origin of life modeling taken into consideration. Using a combination of theoretical modeling and numerical simulations, we look for minimal conditions for which our general chemical model exhibits spontaneous mirror symmetry breaking. We show that our model spontaneously breaks mirror symmetry in various catalytic configurations that only involve a small number of catalyzed reactions and nothing else. Of particular importance is that mirror symmetry breaking occurs in our model without the need for single-step autocatalytis or mutual inhibition, which may be of relevance for prebiotic chemistry.
Collapse
Affiliation(s)
| | - David Hochberg
- Department of Molecular Evolution, Centro de Astrobiología (CSIC-INTA), Ctra. Ajalvir Km. 4, 28850, Torrejón de Ardóz, Madrid, Spain
| |
Collapse
|
15
|
Protasov E, Martinov M, Sinauridze E, Vitvitsky V, Ataullakhanov F. Prediction of Oscillations in Glycolysis in Ethanol-Consuming Erythrocyte-Bioreactors. Int J Mol Sci 2023; 24:10124. [PMID: 37373271 DOI: 10.3390/ijms241210124] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/27/2023] [Revised: 06/03/2023] [Accepted: 06/11/2023] [Indexed: 06/29/2023] Open
Abstract
A mathematical model of energy metabolism in erythrocyte-bioreactors loaded with alcohol dehydrogenase and acetaldehyde dehydrogenase was constructed and analyzed. Such erythrocytes can convert ethanol to acetate using intracellular NAD and can therefore be used to treat alcohol intoxication. Analysis of the model revealed that the rate of ethanol consumption by the erythrocyte-bioreactors increases proportionally to the activity of incorporated ethanol-consuming enzymes until their activity reaches a specific threshold level. When the ethanol-consuming enzyme activity exceeds this threshold, the steady state in the model becomes unstable and the model switches to an oscillation mode caused by the competition between glyceraldehyde phosphate dehydrogenase and ethanol-consuming enzymes for NAD. The amplitude and period of metabolite oscillations first increase with the increase in the activity of the encapsulated enzymes. A further increase in these activities leads to a loss of the glycolysis steady state, and a permanent accumulation of glycolytic intermediates. The oscillation mode and the loss of the steady state can lead to the osmotic destruction of erythrocyte-bioreactors due to an accumulation of intracellular metabolites. Our results demonstrate that the interaction of enzymes encapsulated in erythrocyte-bioreactors with erythrocyte metabolism should be taken into account in order to achieve the optimal efficacy of these bioreactors.
Collapse
Affiliation(s)
- Evgeniy Protasov
- Laboratory of Biophysics, Dmitriy Rogachev National Medical Research Center of Pediatric Hematology, Oncology, and Immunology, Ministry of Healthcare, Samora Mashel Str., 1, GSP-7, Moscow 117198, Russia
- Laboratory of Physiology and Biophysics of the Cell, Center for Theoretical Problems of Physicochemical Pharmacology, Russian Academy of Sciences, Srednyaya Kalitnikovskaya Str., 30, Moscow 109029, Russia
| | - Michael Martinov
- Laboratory of Physiology and Biophysics of the Cell, Center for Theoretical Problems of Physicochemical Pharmacology, Russian Academy of Sciences, Srednyaya Kalitnikovskaya Str., 30, Moscow 109029, Russia
| | - Elena Sinauridze
- Laboratory of Biophysics, Dmitriy Rogachev National Medical Research Center of Pediatric Hematology, Oncology, and Immunology, Ministry of Healthcare, Samora Mashel Str., 1, GSP-7, Moscow 117198, Russia
- Laboratory of Physiology and Biophysics of the Cell, Center for Theoretical Problems of Physicochemical Pharmacology, Russian Academy of Sciences, Srednyaya Kalitnikovskaya Str., 30, Moscow 109029, Russia
| | - Victor Vitvitsky
- Laboratory of Physiology and Biophysics of the Cell, Center for Theoretical Problems of Physicochemical Pharmacology, Russian Academy of Sciences, Srednyaya Kalitnikovskaya Str., 30, Moscow 109029, Russia
| | - Fazoil Ataullakhanov
- Laboratory of Biophysics, Dmitriy Rogachev National Medical Research Center of Pediatric Hematology, Oncology, and Immunology, Ministry of Healthcare, Samora Mashel Str., 1, GSP-7, Moscow 117198, Russia
- Laboratory of Physiology and Biophysics of the Cell, Center for Theoretical Problems of Physicochemical Pharmacology, Russian Academy of Sciences, Srednyaya Kalitnikovskaya Str., 30, Moscow 109029, Russia
- Department of Molecular and Translational Medicine, Moscow Institute of Physics and Technology, Institutskiy Per., 9, Dolgoprudny 141701, Russia
- Perelman School of Medicine, University of Pennsylvania, 3400 Civic Center Blvd., Philadelphia, PA 19104, USA
| |
Collapse
|
16
|
Murray PJ. Autoregulation of Transcription and Translation: A Qualitative Analysis. Bull Math Biol 2023; 85:57. [PMID: 37233955 DOI: 10.1007/s11538-023-01143-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2022] [Accepted: 03/09/2023] [Indexed: 05/27/2023]
Abstract
The regulation of both mRNA transcription and translation by down-stream gene products allows for a range of rich dynamical behaviours (e.g. homeostatic, oscillatory, excitability and intermittent solutions). Here, qualitative analysis is applied to an existing model of a gene regulatory network in which a protein dimer inhibits its own transcription and upregulates its own translation rate. It is demonstrated that the model possesses a unique steady state, conditions are derived under which limit cycle solutions arise and estimates are provided for the oscillator period in the limiting case of a relaxation oscillator. The analysis demonstrates that oscillations can arise only if mRNA is more stable than protein and the effect of nonlinear translation inhibition is sufficiently strong. Moreover, it is shown that the oscillation period can vary non-monotonically with transcription rate. Thus the proposed framework can provide an explanation for observed species-specific dependency of segmentation clock period on Notch signalling activity. Finally, this study facilitates the application of the proposed model to more general biological settings where post transcriptional regulation effects are likely important.
Collapse
|
17
|
Sahin A, Weilandt DR, Hatzimanikatis V. Optimal enzyme utilization suggests that concentrations and thermodynamics determine binding mechanisms and enzyme saturations. Nat Commun 2023; 14:2618. [PMID: 37147292 PMCID: PMC10162984 DOI: 10.1038/s41467-023-38159-4] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2022] [Accepted: 04/19/2023] [Indexed: 05/07/2023] Open
Abstract
Deciphering the metabolic functions of organisms requires understanding the dynamic responses of living cells upon genetic and environmental perturbations, which in turn can be inferred from enzymatic activity. In this work, we investigate the optimal modes of operation for enzymes in terms of the evolutionary pressure driving them toward increased catalytic efficiency. We develop a framework using a mixed-integer formulation to assess the distribution of thermodynamic forces and enzyme states, providing detailed insights into the enzymatic mode of operation. We use this framework to explore Michaelis-Menten and random-ordered multi-substrate mechanisms. We show that optimal enzyme utilization is achieved by unique or alternative operating modes dependent on reactant concentrations. We find that in a bimolecular enzyme reaction, the random mechanism is optimal over any other ordered mechanism under physiological conditions. Our framework can investigate the optimal catalytic properties of complex enzyme mechanisms. It can further guide the directed evolution of enzymes and fill in the knowledge gaps in enzyme kinetics.
Collapse
Affiliation(s)
- Asli Sahin
- Laboratory of Computational Systems Biotechnology, Ecole Polytechnique Federale de Lausanne (EPFL), 1015, Lausanne, Switzerland
| | - Daniel R Weilandt
- Laboratory of Computational Systems Biotechnology, Ecole Polytechnique Federale de Lausanne (EPFL), 1015, Lausanne, Switzerland
- Department of Chemistry and Lewis-Sigler Institute for Integrative Genomics, Princeton University, Princeton, NJ, USA
| | - Vassily Hatzimanikatis
- Laboratory of Computational Systems Biotechnology, Ecole Polytechnique Federale de Lausanne (EPFL), 1015, Lausanne, Switzerland.
| |
Collapse
|
18
|
Kumar P, Gangopadhyay G. Glycolytic Wave Patterns in a Simple Reaction-diffusion System with Inhomogeneous Influx: Dynamic Transitions. Chemphyschem 2023; 24:e202200643. [PMID: 36478341 DOI: 10.1002/cphc.202200643] [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: 08/26/2022] [Revised: 11/21/2022] [Indexed: 12/12/2022]
Abstract
An inhomogeneous profile of chemostatted species generates a rich variety of patterns in glycolytic waves depicted in a Selkov reaction-diffusion framework here. A key role played by diffusion amplitude and symmetry in the chemostatted species profile in dictating the fate of local spatial dynamics involving periodic, quasiperiodic, and chaotic patterns and transitions among them are investigated systematically. More importantly, various dynamic transitions, including wave propagation direction changes, are illustrated in interesting situations. Besides numerical results, our analytical formulation of the amplitude equation connecting complex Ginzburg-Landau and Lambda-omega representation shed light on the phase dynamics of the system. This systematic study of the glycolytic reaction-diffusion wave is in line with previous experimental results in open spatial reactor and will provide a knowledge about the dynamics that shape and control biological information processing and related phenomena.
Collapse
Affiliation(s)
- Premashis Kumar
- S. N. Bose National Centre For Basic Sciences, Block-JD, Sector-III, Salt Lake, Kolkata, 700 106, India
| | - Gautam Gangopadhyay
- S. N. Bose National Centre For Basic Sciences, Block-JD, Sector-III, Salt Lake, Kolkata, 700 106, India
| |
Collapse
|
19
|
Muñuzuri AP, Gagnon JS, Pérez-Mercader J. Numerical and renormalization group analysis of the phase diagram of a stochastic cubic autocatalytic reaction-diffusion system. Phys Rev E 2023; 107:034213. [PMID: 37073065 DOI: 10.1103/physreve.107.034213] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/08/2022] [Accepted: 03/14/2023] [Indexed: 04/20/2023]
Abstract
The renormalization group is a set of tools that can be used to incorporate the effect of fluctuations in a dynamical system as a rescaling of the system's parameters. Here, we apply the renormalization group to a pattern-forming stochastic cubic autocatalytic reaction-diffusion model and compare its predictions with numerical simulations. Our results demonstrate a good agreement within the range of validity of the theory and show that external noise can be used as a control parameter in such systems.
Collapse
Affiliation(s)
- Alberto P Muñuzuri
- Department of Earth and Planetary Sciences, Harvard University, Cambridge, Massachusetts 02138-1204, USA
- Group of Nonlinear Physics, University of Santiago de Compostela, 15706 Santiago de Compostela, Spain
- Galician Center for Mathematical Research and Technology (CITMAga), 15782 Santiago de Compostela, Spain
| | - Jean-Sébastien Gagnon
- Department of Earth and Planetary Sciences, Harvard University, Cambridge, Massachusetts 02138-1204, USA
- Department of Physics, Norwich University, Northfield, Vermont 05663, USA
| | - Juan Pérez-Mercader
- Department of Earth and Planetary Sciences, Harvard University, Cambridge, Massachusetts 02138-1204, USA
- Santa Fe Institute, Santa Fe, New Mexico 87501, USA
| |
Collapse
|
20
|
Gagrani P, Smith E. Action functional gradient descent algorithm for estimating escape paths in stochastic chemical reaction networks. Phys Rev E 2023; 107:034305. [PMID: 37072963 DOI: 10.1103/physreve.107.034305] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2022] [Accepted: 02/21/2023] [Indexed: 04/20/2023]
Abstract
We first derive the Hamilton-Jacobi theory underlying continuous-time Markov processes, and then we use the construction to develop a variational algorithm for estimating escape (least improbable or first passage) paths for a generic stochastic chemical reaction network that exhibits multiple fixed points. The design of our algorithm is such that it is independent of the underlying dimensionality of the system, the discretization control parameters are updated toward the continuum limit, and there is an easy-to-calculate measure for the correctness of its solution. We consider several applications of the algorithm and verify them against computationally expensive means such as the shooting method and stochastic simulation. While we employ theoretical techniques from mathematical physics, numerical optimization and chemical reaction network theory, we hope that our work finds practical applications with an inter-disciplinary audience including chemists, biologists, optimal control theorists and game theorists.
Collapse
Affiliation(s)
- Praful Gagrani
- Department of Physics, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
- Wisconsin Institute for Discovery, University of Wisconsin-Madison, Madison, Wisconsin 53715, USA
| | - Eric Smith
- Department of Physics, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
- Department of Biology, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
- Earth-Life Science Institute, Tokyo Institute of Technology, Tokyo 152-8550, Japan
- Santa Fe Institute, 1399 Hyde Park Road, Santa Fe, New Mexico 87501, USA
- Ronin Institute, 127 Haddon Place, Montclair, New Jersey 07043, USA
| |
Collapse
|
21
|
Systematic analysis of negative and positive feedback loops for robustness and temperature compensation in circadian rhythms. NPJ Syst Biol Appl 2023; 9:5. [PMID: 36774353 PMCID: PMC9922291 DOI: 10.1038/s41540-023-00268-7] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2022] [Accepted: 01/30/2023] [Indexed: 02/13/2023] Open
Abstract
Temperature compensation and robustness to biological noise are two key characteristics of the circadian clock. These features allow the circadian pacemaker to maintain a steady oscillation in a wide range of environmental conditions. The presence of a time-delayed negative feedback loop in the regulatory network generates autonomous circadian oscillations in eukaryotic systems. In comparison, the circadian clock of cyanobacteria is controlled by a strong positive feedback loop. Positive feedback loops with substrate depletion can also generate oscillations, inspiring other circadian clock models. What makes a circadian oscillatory network robust to extrinsic noise is unclear. We investigated four basic circadian oscillators with negative, positive, and combinations of positive and negative feedback loops to explore network features necessary for circadian clock resilience. We discovered that the negative feedback loop system performs the best in compensating temperature changes. We also show that a positive feedback loop can reduce extrinsic noise in periods of circadian oscillators, while intrinsic noise is reduced by negative feedback loops.
Collapse
|
22
|
Mauksch M. Spontaneous emergence of enantioenriched chiral aldol reaction products from Achiral precursors in solution and origin of biological homochirality of sugars: a first-principles study. Phys Chem Chem Phys 2023; 25:1734-1754. [PMID: 36594779 DOI: 10.1039/d2cp04285a] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
Experimental reports about observation of spontaneous mirror symmetry breaking and chiral amplification in stereoselective Mannich and aldol reactions, run under fully achiral initial conditions, have drawn a lot of attention, fuelled partly by the role these reactions could have played in chemical evolution as a cause for still puzzling observed homochirality of biomolecules, often considered a prerequisite for the origin of life. We have now revisited this still unresolved problem, using DFT computation of all combinatorially possible transition states and numerical solution of complete set of resulting coupled kinetic rate equations to model the aldol reaction rigorously "from the first principles" and without making any a priori assumptions. Spontaneous mirror symmetry breaking in this autocatalytic, reversible, closed and homogenous system is explained by a supercritical pitchfork bifurcation, occurring in concentrations of enantiomers due to time-delayed kinetic instability of racemic composition of reaction mixture, when reactants are initially provided in non-stoichiometric quantities. Same process, taking place under similar conditions in primordial "soup" of chemicals, might conceivably explain origin of biological homochirality of sugar molecules on early earth billions of years ago. Our results suggest that seemingly innocuous chemical reactions could exhibit unexpected and counter-intuitive emergent behaviour, when initial conditions are appropriately chosen. Chiral amplification in self-catalyzed aldol reaction occurs during approach of thermodynamic equilibrium in accord with principle of microscopic reversibility and second law of thermodynamics.
Collapse
Affiliation(s)
- Michael Mauksch
- Department of Chemistry and Pharmacy, Institute of Theoretical Chemistry, Computer Chemistry Center, Nägelsbachstrasse 25a, 91052 Erlangen, Germany.
| |
Collapse
|
23
|
Nath SS, Villadsen J. Modeling dynamics of chemical reaction networks using electrical analogs: Application to autocatalytic reactions. CHEMICAL ENGINEERING JOURNAL ADVANCES 2022. [DOI: 10.1016/j.ceja.2022.100374] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022] Open
|
24
|
Joshi B, Craciun G. Foundations of static and dynamic absolute concentration robustness. J Math Biol 2022; 85:53. [PMID: 36243796 DOI: 10.1007/s00285-022-01823-2] [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: 05/17/2021] [Revised: 09/14/2022] [Accepted: 10/03/2022] [Indexed: 10/17/2022]
Abstract
Absolute Concentration Robustness (ACR) was introduced by Shinar and Feinberg (Science 327:1389-1391, 2010) as robustness of equilibrium species concentration in a mass action dynamical system. Their aim was to devise a mathematical condition that will ensure robustness in the function of the biological system being modeled. The robustness of function rests on what we refer to as empirical robustness-the concentration of a species remains unvarying, when measured in the long run, across arbitrary initial conditions. Even simple examples show that the ACR notion introduced in Shinar and Feinberg (Science 327:1389-1391, 2010) (here referred to as static ACR) is neither necessary nor sufficient for empirical robustness. To make a stronger connection with empirical robustness, we define dynamic ACR, a property related to long-term, global dynamics, rather than only to equilibrium behavior. We discuss general dynamical systems with dynamic ACR properties as well as parametrized families of dynamical systems related to reaction networks. We find necessary and sufficient conditions for dynamic ACR in complex balanced reaction networks, a class of networks that is central to the theory of reaction networks.
Collapse
Affiliation(s)
- Badal Joshi
- Department of Mathematics, California State University San Marcos, San Marcos, USA.
| | - Gheorghe Craciun
- Departments of Mathematics and Biomolecular Chemistry, University of Wisconsin-Madison, Madison, USA
| |
Collapse
|
25
|
Casani-Galdon P, Garcia-Ojalvo J. Signaling oscillations: Molecular mechanisms and functional roles. Curr Opin Cell Biol 2022; 78:102130. [PMID: 36130445 DOI: 10.1016/j.ceb.2022.102130] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/14/2022] [Revised: 07/21/2022] [Accepted: 08/19/2022] [Indexed: 01/31/2023]
Abstract
Mounting evidence shows that oscillatory activity is widespread in cell signaling. Here, we review some of this recent evidence, focusing on both the molecular mechanisms that potentially underlie such dynamical behavior, and the potential advantages that signaling oscillations might have in cell function. The biological processes considered include cellular differentiation and tissue maintenance, intermittent responses in pluripotent stem cells, and collective cell migration during wound healing. With the aid of mathematical modeling, we review recent examples in which delayed negative feedback has been seen to act as a unifying principle that underpins this wide variety of phenomena.
Collapse
Affiliation(s)
- Pablo Casani-Galdon
- Department of Medicine and Life Sciences, Universitat Pompeu Fabra, Barcelona Biomedical Research Park, Dr. Aiguader 88, 08003 Barcelona, Spain
| | - Jordi Garcia-Ojalvo
- Department of Medicine and Life Sciences, Universitat Pompeu Fabra, Barcelona Biomedical Research Park, Dr. Aiguader 88, 08003 Barcelona, Spain.
| |
Collapse
|
26
|
Bashkirtseva I, Pankratov A, Ryashko L. Noise-induced formation of heterogeneous patterns in the Turing stability zones of diffusion systems. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2022; 34:444001. [PMID: 36001986 DOI: 10.1088/1361-648x/ac8c77] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/31/2022] [Accepted: 08/24/2022] [Indexed: 06/15/2023]
Abstract
We study a phenomenon of stochastic generation of waveform patterns for reaction-diffusion systems in the Turing stability zone where the homogeneous equilibrium is a single attractor. In this analysis, we use a distributed variant of the Selkov glycolytic model with diffusion and random forcing. It is shown that in the Turing stability zone, random disturbances can induce a diversity of metastable spatial patterns with different waveforms. We carry out the parametric analysis of statistical characteristics of evolution of these patterns, and reveal the dominant patterns in the stochastic flow of mixed spatial structures.
Collapse
Affiliation(s)
| | | | - Lev Ryashko
- Ural Federal University, Lenina, 51, 620000 Ekaterinburg, Russia
| |
Collapse
|
27
|
Nijholt E, DeVille L. Dynamical systems defined on simplicial complexes: Symmetries, conjugacies, and invariant subspaces. CHAOS (WOODBURY, N.Y.) 2022; 32:093131. [PMID: 36182350 DOI: 10.1063/5.0093842] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/30/2022] [Accepted: 08/15/2022] [Indexed: 06/16/2023]
Abstract
We consider the general model for dynamical systems defined on a simplicial complex. We describe the conjugacy classes of these systems and show how symmetries in a given simplicial complex manifest in the dynamics defined thereon, especially with regard to invariant subspaces in the dynamics.
Collapse
Affiliation(s)
- Eddie Nijholt
- ICMC São Carlos, Universidade de São Paulo, Av. Trab. São Carlense 400, São Carlos SP 13566-590, Brasil
| | - Lee DeVille
- Department of Mathematics, University of Illinois, 1409 W. Green St., Urbana, Illinois 61801, USA
| |
Collapse
|
28
|
Bartlett S, Louapre D. Provenance of life: Chemical autonomous agents surviving through associative learning. Phys Rev E 2022; 106:034401. [PMID: 36266823 DOI: 10.1103/physreve.106.034401] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2022] [Accepted: 06/21/2022] [Indexed: 06/16/2023]
Abstract
We present a benchmark study of autonomous, chemical agents exhibiting associative learning of an environmental feature. Associative learning systems have been widely studied in cognitive science and artificial intelligence but are most commonly implemented in highly complex or carefully engineered systems, such as animal brains, artificial neural networks, DNA computing systems, and gene regulatory networks, among others. The ability to encode environmental information and use it to make simple predictions is a benchmark of biological resilience and underpins a plethora of adaptive responses in the living hierarchy, spanning prey animal species anticipating the arrival of predators to epigenetic systems in microorganisms learning environmental correlations. Given the ubiquitous and essential presence of learning behaviors in the biosphere, we aimed to explore whether simple, nonliving dissipative structures could also exhibit associative learning. Inspired by previous modeling of associative learning in chemical networks, we simulated simple systems composed of long- and short-term memory chemical species that could encode the presence or absence of temporal correlations between two external species. The ability to learn this association was implemented in Gray-Scott reaction-diffusion spots, emergent chemical patterns that exhibit self-replication and homeostasis. With the novel ability of associative learning, we demonstrate that simple chemical patterns can exhibit a broad repertoire of lifelike behavior, paving the way for in vitro studies of autonomous chemical learning systems, with potential relevance to artificial life, origins of life, and systems chemistry. The experimental realization of these learning behaviors in protocell or coacervate systems could advance a new research direction in astrobiology, since our system significantly reduces the lower bound on the required complexity for autonomous chemical learning.
Collapse
Affiliation(s)
- Stuart Bartlett
- Division of Geological and Planetary Sciences, California Institute of Technology, Pasadena, California 91125, USA and Earth-Life Science Institute, Tokyo Institute of Technology, Tokyo 152-8550, Japan
| | - David Louapre
- Ubisoft Entertainment, 94160 Saint-Mandé, France and Science Étonnante, 75014 Paris, France†
| |
Collapse
|
29
|
Cisneros-Ake LA, Gonzalez-Rodriguez JC, González-Ramírez LR. Turing Instabilities and Rotating Spiral Waves in Glycolytic Processes. Bull Math Biol 2022; 84:100. [PMID: 35951127 DOI: 10.1007/s11538-022-01060-0] [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: 01/28/2022] [Accepted: 07/25/2022] [Indexed: 11/29/2022]
Abstract
We study single-frequency oscillations and pattern formation in the glycolytic process modeled by a reduction in the well-known Sel'kov's equations (Sel'kov in Eur J Biochem 4:79, 1968), which describe, in the whole cell, the phosphofructokinase enzyme reaction. By using averaging theory, we establish the existence conditions for limit cycles and their limiting average radius in the kinetic reaction equations. We analytically establish conditions on the model parameters for the appearance of unstable nonlinear modes seeding the formation of two-dimensional patterns in the form of classical spots and stripes. We also establish the existence of a Hopf bifurcation, which characterizes the reaction dynamics, producing glycolytic rotating spiral waves. We numerically establish parameter regions for the existence of these spiral waves and address their linear stability. We show that as the model tends toward a suppression of the relative source rate, the spiral wave solution loses stability. All our findings are validated by full numerical simulations of the model equations. Finally, we discuss in vitro evidence of spatiotemporal activity patterns found in glycolytic experiments, and propose plausible biological implications of our model results.
Collapse
Affiliation(s)
- Luis A Cisneros-Ake
- Departamento de Matemáticas, ESFM, Instituto Politécnico Nacional, Unidad Profesional Adolfo López Mateos Edificio 9, 07738, Mexico, Mexico.
| | - Juan C Gonzalez-Rodriguez
- Posgrado en Ciencias Fisicomatemáticas, ESFM, Instituto Politécnico Nacional, Unidad Profesional Adolfo López Mateos Edificio 9, 07738, Mexico, Mexico
| | - Laura R González-Ramírez
- Departamento de Matemáticas, ESFM, Instituto Politécnico Nacional, Unidad Profesional Adolfo López Mateos Edificio 9, 07738, Mexico, Mexico
| |
Collapse
|
30
|
Lo-Thong-Viramoutou O, Charton P, Cadet XF, Grondin-Perez B, Saavedra E, Damour C, Cadet F. Non-linearity of Metabolic Pathways Critically Influences the Choice of Machine Learning Model. Front Artif Intell 2022; 5:744755. [PMID: 35757298 PMCID: PMC9226554 DOI: 10.3389/frai.2022.744755] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2021] [Accepted: 04/29/2022] [Indexed: 11/13/2022] Open
Abstract
The use of machine learning (ML) in life sciences has gained wide interest over the past years, as it speeds up the development of high performing models. Important modeling tools in biology have proven their worth for pathway design, such as mechanistic models and metabolic networks, as they allow better understanding of mechanisms involved in the functioning of organisms. However, little has been done on the use of ML to model metabolic pathways, and the degree of non-linearity associated with them is not clear. Here, we report the construction of different metabolic pathways with several linear and non-linear ML models. Different types of data are used; they lead to the prediction of important biological data, such as pathway flux and final product concentration. A comparison reveals that the data features impact model performance and highlight the effectiveness of non-linear models (e.g., QRF: RMSE = 0.021 nmol·min-1 and R2 = 1 vs. Bayesian GLM: RMSE = 1.379 nmol·min-1 R2 = 0.823). It turns out that the greater the degree of non-linearity of the pathway, the better suited a non-linear model will be. Therefore, a decision-making support for pathway modeling is established. These findings generally support the hypothesis that non-linear aspects predominate within the metabolic pathways. This must be taken into account when devising possible applications of these pathways for the identification of biomarkers of diseases (e.g., infections, cancer, neurodegenerative diseases) or the optimization of industrial production processes.
Collapse
Affiliation(s)
- Ophélie Lo-Thong-Viramoutou
- University of Paris, BIGR—Biologie Intégrée du Globule Rouge, Inserm, UMR_S1134, Paris, France
- Laboratory of Excellence GR-Ex, Paris, France
- Laboratory DSIMB, UMR_S1134, BIGR, Inserm, Faculty of Sciences and Technology, University of La Reunion, Saint-Denis, France
| | - Philippe Charton
- University of Paris, BIGR—Biologie Intégrée du Globule Rouge, Inserm, UMR_S1134, Paris, France
- Laboratory of Excellence GR-Ex, Paris, France
- Laboratory DSIMB, UMR_S1134, BIGR, Inserm, Faculty of Sciences and Technology, University of La Reunion, Saint-Denis, France
| | | | - Brigitte Grondin-Perez
- EnergyLab, EA 4079, Faculty of Sciences and Technology, University of La Reunion, Saint-Denis, France
| | - Emma Saavedra
- Departamento de Bioquímica, Instituto Nacional de Cardiología Ignacio Chávez, Mexico City, Mexico
| | - Cédric Damour
- EnergyLab, EA 4079, Faculty of Sciences and Technology, University of La Reunion, Saint-Denis, France
| | - Frédéric Cadet
- University of Paris, BIGR—Biologie Intégrée du Globule Rouge, Inserm, UMR_S1134, Paris, France
- Laboratory of Excellence GR-Ex, Paris, France
- Laboratory DSIMB, UMR_S1134, BIGR, Inserm, Faculty of Sciences and Technology, University of La Reunion, Saint-Denis, France
| |
Collapse
|
31
|
Remlein B, Weissmann V, Seifert U. Coherence of oscillations in the weak-noise limit. Phys Rev E 2022; 105:064101. [PMID: 35854576 DOI: 10.1103/physreve.105.064101] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/14/2022] [Accepted: 05/13/2022] [Indexed: 06/15/2023]
Abstract
In a noisy environment, oscillations lose their coherence, which can be characterized by a quality factor. We determine this quality factor for oscillations arising from a driven Fokker-Planck dynamics along a periodic one-dimensional potential analytically in the weak-noise limit. With this expression, we can prove for this continuum model the analog of an upper bound that has been conjectured for the coherence of oscillations in discrete Markov network models. We show that our approach can also be adapted to motion along a noisy two-dimensional limit cycle. Specifically, we apply our scheme to the noisy Stuart-Landau oscillator and the thermodynamically consistent Brusselator as a simple model for a chemical clock. Our approach thus complements the fairly sophisticated extant general framework based on techniques from Hamilton-Jacobi theory with which we compare our results numerically.
Collapse
Affiliation(s)
- Benedikt Remlein
- II. Institut für Theoretische Physik, Universität Stuttgart, 70550 Stuttgart, Germany
| | - Volker Weissmann
- II. Institut für Theoretische Physik, Universität Stuttgart, 70550 Stuttgart, Germany
| | - Udo Seifert
- II. Institut für Theoretische Physik, Universität Stuttgart, 70550 Stuttgart, Germany
| |
Collapse
|
32
|
Khan E, Saghafi S, Diekman CO, Rotstein HG. The emergence of polyglot entrainment responses to periodic inputs in vicinities of Hopf bifurcations in slow-fast systems. CHAOS (WOODBURY, N.Y.) 2022; 32:063137. [PMID: 35778129 DOI: 10.1063/5.0079198] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/18/2021] [Accepted: 04/04/2022] [Indexed: 06/15/2023]
Abstract
Several distinct entrainment patterns can occur in the FitzHugh-Nagumo (FHN) model under external periodic forcing. Investigating the FHN model under different types of periodic forcing reveals the existence of multiple disconnected 1:1 entrainment segments for constant, low enough values of the input amplitude when the unforced system is in the vicinity of a Hopf bifurcation. This entrainment structure is termed polyglot to distinguish it from the single 1:1 entrainment region (monoglot) structure typically observed in Arnold tongue diagrams. The emergence of polyglot entrainment is then explained using phase-plane analysis and other dynamical system tools. Entrainment results are investigated for other slow-fast systems of neuronal, circadian, and glycolytic oscillations. Exploring these models, we found that polyglot entrainment structure (multiple 1:1 regions) is observed when the unforced system is in the vicinity of a Hopf bifurcation and the Hopf point is located near a knee of a cubic-like nullcline.
Collapse
Affiliation(s)
- Emel Khan
- Department of Mathematical Sciences, New Jersey Institute of Technology, Newark, New Jersey 07102, USA
| | - Soheil Saghafi
- Department of Mathematical Sciences, New Jersey Institute of Technology, Newark, New Jersey 07102, USA
| | - Casey O Diekman
- Department of Mathematical Sciences, New Jersey Institute of Technology, Newark, New Jersey 07102, USA
| | - Horacio G Rotstein
- Federated Department of Biological Sciences, New Jersey Institute of Technology & Rutgers University, Newark, New Jersey 07102, USA
| |
Collapse
|
33
|
Unified representation of Life's basic properties by a 3-species Stochastic Cubic Autocatalytic Reaction-Diffusion system of equations. Phys Life Rev 2022; 41:64-83. [DOI: 10.1016/j.plrev.2022.03.003] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2022] [Accepted: 03/29/2022] [Indexed: 01/01/2023]
|
34
|
Dhatt S, Chaudhury P. Study of oscillatory dynamics in a Selkov glycolytic model using sensitivity analysis. INDIAN JOURNAL OF PHYSICS 2022; 96:1649-1654. [DOI: 10.1007/s12648-021-02102-4] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/28/2020] [Accepted: 04/05/2021] [Indexed: 07/19/2023]
|
35
|
Saha S, Gangopadhyay G, Ray DS. Universality in bio-rhythms: A perspective from nonlinear dynamics. J Biosci 2022. [DOI: 10.1007/s12038-021-00249-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
|
36
|
On the Dynamics of Higgins–Selkov, Selkov and Brusellator Oscillators. Symmetry (Basel) 2022. [DOI: 10.3390/sym14030438] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/10/2022] Open
Abstract
A complete algebraic characterization of the first integrals of the Higgins–Selkov, Selkov and Brusellator oscillators is given here. The existence of symmetries sometimes forces the existence of such first integrals. The nonexistence of centers for such oscillators is also proved. In order to determine the Puiseux integrability of such systems, the multiple Puiseux solutions are also studied.
Collapse
|
37
|
Lövfors W, Ekström J, Jönsson C, Strålfors P, Cedersund G, Nyman E. A systems biology analysis of lipolysis and fatty acid release from adipocytes in vitro and from adipose tissue in vivo. PLoS One 2021; 16:e0261681. [PMID: 34972146 PMCID: PMC8719686 DOI: 10.1371/journal.pone.0261681] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/16/2021] [Accepted: 12/07/2021] [Indexed: 12/03/2022] Open
Abstract
Lipolysis and the release of fatty acids to supply energy fuel to other organs, such as between meals, during exercise, and starvation, are fundamental functions of the adipose tissue. The intracellular lipolytic pathway in adipocytes is activated by adrenaline and noradrenaline, and inhibited by insulin. Circulating fatty acids are elevated in type 2 diabetic individuals. The mechanisms behind this elevation are not fully known, and to increase the knowledge a link between the systemic circulation and intracellular lipolysis is key. However, data on lipolysis and knowledge from in vitro systems have not been linked to corresponding in vivo data and knowledge in vivo. Here, we use mathematical modelling to provide such a link. We examine mechanisms of insulin action by combining in vivo and in vitro data into an integrated mathematical model that can explain all data. Furthermore, the model can describe independent data not used for training the model. We show the usefulness of the model by simulating new and more challenging experimental setups in silico, e.g. the extracellular concentration of fatty acids during an insulin clamp, and the difference in such simulations between individuals with and without type 2 diabetes. Our work provides a new platform for model-based analysis of adipose tissue lipolysis, under both non-diabetic and type 2 diabetic conditions.
Collapse
Affiliation(s)
- William Lövfors
- Department of Biomedical Engineering, Linköping University, Linköping, Sweden
- Department of Mathematics, Linköping University, Linköping, Sweden
| | - Jona Ekström
- Department of Biomedical Engineering, Linköping University, Linköping, Sweden
| | - Cecilia Jönsson
- Department of Biomedical and Clinical Sciences, Linköping University, Linköping, Sweden
| | - Peter Strålfors
- Department of Biomedical and Clinical Sciences, Linköping University, Linköping, Sweden
| | - Gunnar Cedersund
- Department of Biomedical Engineering, Linköping University, Linköping, Sweden
- Center for Medical Image Science and Visualization (CMIV), Linköping University, Linköping, Sweden
| | - Elin Nyman
- Department of Biomedical Engineering, Linköping University, Linköping, Sweden
| |
Collapse
|
38
|
Gomez D, Iyaniwura S, Paquin-Lefebvre F, Ward MJ. Pattern forming systems coupling linear bulk diffusion to dynamically active membranes or cells. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2021; 379:20200276. [PMID: 34743601 DOI: 10.1098/rsta.2020.0276] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/03/2023]
Abstract
Some analytical and numerical results are presented for pattern formation properties associated with novel types of reaction-diffusion (RD) systems that involve the coupling of bulk diffusion in the interior of a multi-dimensional spatial domain to nonlinear processes that occur either on the domain boundary or within localized compartments that are confined within the domain. The class of bulk-membrane system considered herein is derived from an asymptotic analysis in the limit of small thickness of a thin domain that surrounds the bulk medium. When the bulk domain is a two-dimensional disk, a weakly nonlinear analysis is used to characterize Turing and Hopf bifurcations that can arise from the linearization around a radially symmetric, but spatially non-uniform, steady-state of the bulk-membrane system. In a singularly perturbed limit, the existence and linear stability of localized membrane-bound spike patterns is analysed for a Gierer-Meinhardt activator-inhibitor model that includes bulk coupling. Finally, the emergence of collective intracellular oscillations is studied for a class of PDE-ODE bulk-cell model in a bounded two-dimensional domain that contains spatially localized, but dynamically active, circular cells that are coupled through a linear bulk diffusion field. Applications of such coupled bulk-membrane or bulk-cell systems to some biological systems are outlined, and some open problems in this area are discussed. This article is part of the theme issue 'Recent progress and open frontiers in Turing's theory of morphogenesis'.
Collapse
Affiliation(s)
- D Gomez
- Dept. of Mathematics, Univ. of British Columbia, Vancouver, British Columbia, Canada V6T 1Z2
| | - S Iyaniwura
- Dept. of Mathematics, Univ. of British Columbia, Vancouver, British Columbia, Canada V6T 1Z2
| | - F Paquin-Lefebvre
- Dept. of Mathematics, Univ. of British Columbia, Vancouver, British Columbia, Canada V6T 1Z2
| | - M J Ward
- Dept. of Mathematics, Univ. of British Columbia, Vancouver, British Columbia, Canada V6T 1Z2
| |
Collapse
|
39
|
Al Saadi F, Champneys A. Unified framework for localized patterns in reaction-diffusion systems; the Gray-Scott and Gierer-Meinhardt cases. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2021; 379:20200277. [PMID: 34743600 DOI: 10.1098/rsta.2020.0277] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/03/2023]
Abstract
A recent study of canonical activator-inhibitor Schnakenberg-like models posed on an infinite line is extended to include models, such as Gray-Scott, with bistability of homogeneous equilibria. A homotopy is studied that takes a Schnakenberg-like glycolysis model to the Gray-Scott model. Numerical continuation is used to understand the complete sequence of transitions to two-parameter bifurcation diagrams within the localized pattern parameter regime as the homotopy parameter varies. Several distinct codimension-two bifurcations are discovered including cusp and quadruple zero points for homogeneous steady states, a degenerate heteroclinic connection and a change in connectedness of the homoclinic snaking structure. The analysis is repeated for the Gierer-Meinhardt system, which lies outside the canonical framework. Similar transitions are found under homotopy between bifurcation diagrams for the case where there is a constant feed in the active field, to it being in the inactive field. Wider implications of the results are discussed for other pattern-formation systems arising as models of natural phenomena. This article is part of the theme issue 'Recent progress and open frontiers in Turing's theory of morphogenesis'.
Collapse
Affiliation(s)
- Fahad Al Saadi
- Department of Engineering Mathematics, University of Bristol, Bristol BS8 1UB, UK
- Department of Systems Engineering, Military Technological College, Muscat, Oman
| | - Alan Champneys
- Department of Engineering Mathematics, University of Bristol, Bristol BS8 1UB, UK
| |
Collapse
|
40
|
Lakhova TN, Kazantsev FV, Lashin SA, Matushkin YG. The finding and researching algorithm for potentially oscillating enzymatic systems. Vavilovskii Zhurnal Genet Selektsii 2021; 25:318-330. [PMID: 34901728 PMCID: PMC8627878 DOI: 10.18699/vj21.035] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/13/2020] [Revised: 03/15/2021] [Accepted: 03/18/2021] [Indexed: 11/19/2022] Open
Abstract
Many processes in living organisms are subject to periodic oscillations at different hierarchical levels of their organization: from molecular-genetic to population and ecological. Oscillatory processes are responsible for cell cycles in both prokaryotes and eukaryotes, for circadian rhythms, for synchronous coupling of respiration with cardiac contractions, etc. Fluctuations in the numbers of organisms in natural populations can be caused by the populations' own properties, their age structure, and ecological relationships with other species. Along with experimental approaches, mathematical and computer modeling is widely used to study oscillating biological systems. This paper presents classical mathematical models that describe oscillatory behavior in biological systems. Methods for the search for oscillatory molecular-genetic systems are presented by the example of their special case - oscillatory enzymatic systems. Factors influencing the cyclic dynamics in living systems, typical not only of the molecular-genetic level, but of higher levels of organization as well, are considered. Application of different ways to describe gene networks for modeling oscillatory molecular-genetic systems is considered, where the most important factor for the emergence of cyclic behavior is the presence of feedback. Techniques for finding potentially oscillatory enzymatic systems are presented. Using the method described in the article, we present and analyze, in a step-by-step manner, first the structural models (graphs) of gene networks and then the reconstruction of the mathematical models and computational experiments with them. Structural models are ideally suited for the tasks of an automatic search for potential oscillating contours (linked subgraphs), whose structure can correspond to the mathematical model of the molecular-genetic system that demonstrates oscillatory behavior in dynamics. At the same time, it is the numerical study of mathematical models for the selected contours that makes it possible to confirm the presence of stable limit cycles in them. As an example of application of the technology, a network of 300 metabolic reactions of the bacterium Escherichia coli was analyzed using mathematical and computer modeling tools. In particular, oscillatory behavior was shown for a loop whose reactions are part of the tryptophan biosynthesis pathway.
Collapse
Affiliation(s)
- T N Lakhova
- Kurchatov Genomics Center of ICG SB RAS, Novosibirsk, Russia
| | - F V Kazantsev
- Kurchatov Genomics Center of ICG SB RAS, Novosibirsk, Russia
| | - S A Lashin
- Kurchatov Genomics Center of ICG SB RAS, Novosibirsk, Russia Novosibirsk State University, Novosibirsk, Russia
| | - Yu G Matushkin
- Institute of Cytology and Genetics of the Siberian Branch of the Russian Academy of Sciences, Novosibirsk, Russia Novosibirsk State University, Novosibirsk, Russia
| |
Collapse
|
41
|
Zhukov KV, Vetcher AA, Gasparuan BA, Shishonin AY. Alteration of Relative Rates of Biodegradation and Regeneration of Cervical Spine Cartilage through the Restoration of Arterial Blood Flow Access to Rhomboid Fossa: A Hypothesis. Polymers (Basel) 2021; 13:polym13234248. [PMID: 34883749 PMCID: PMC8659970 DOI: 10.3390/polym13234248] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/03/2021] [Revised: 11/22/2021] [Accepted: 11/27/2021] [Indexed: 12/12/2022] Open
Abstract
We found the logical way to prove the existence of the mechanism that maintains the rates of biodegradation and regeneration of cervical spine cartilage. We demonstrate, that after we restore access to arterial blood flow through cervical vertebral arteries to rhomboid fossa it causes the prevalence of regeneration over biodegradation. This is in the frames of consideration of the human body as a dissipative structure. Then the recovery of the body should be considered as a reduction of the relative rates of decay below the regeneration ones. Then the recovery of cervical spine cartilage through redirecting of inner dissipative flow depends on the information about oxygen availability that is provided from oxygen detectors in the rhomboid fossa to the cerebellum. Our proposed approach explains already collected data, which satisfies all the scientific requirements. This allows us to draw conclusions that permit reconsidering the way of dealing with multiple chronic diseases.
Collapse
Affiliation(s)
- Kirill V. Zhukov
- Complementary and Integrative Health Clinic of Dr. Shishonin, 5 Yasnogorskaya Str., 117588 Moscow, Russia; (K.V.Z.); (B.A.G.); (A.Y.S.)
| | - Alexandre A. Vetcher
- Complementary and Integrative Health Clinic of Dr. Shishonin, 5 Yasnogorskaya Str., 117588 Moscow, Russia; (K.V.Z.); (B.A.G.); (A.Y.S.)
- Peoples’ Friendship University of Russia (RUDN), 6 Miklukho-Maklaya Str., 117198 Moscow, Russia
- Correspondence:
| | - Bagrat A. Gasparuan
- Complementary and Integrative Health Clinic of Dr. Shishonin, 5 Yasnogorskaya Str., 117588 Moscow, Russia; (K.V.Z.); (B.A.G.); (A.Y.S.)
| | - Alexander Y. Shishonin
- Complementary and Integrative Health Clinic of Dr. Shishonin, 5 Yasnogorskaya Str., 117588 Moscow, Russia; (K.V.Z.); (B.A.G.); (A.Y.S.)
| |
Collapse
|
42
|
Synchronization of the Glycolysis Reaction-Diffusion Model via Linear Control Law. ENTROPY 2021; 23:e23111516. [PMID: 34828214 PMCID: PMC8622022 DOI: 10.3390/e23111516] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 10/14/2021] [Revised: 11/12/2021] [Accepted: 11/13/2021] [Indexed: 11/17/2022]
Abstract
The Selkov system, which is typically employed to model glycolysis phenomena, unveils some rich dynamics and some other complex formations in biochemical reactions. In the present work, the synchronization problem of the glycolysis reaction-diffusion model is handled and examined. In addition, a novel convenient control law is designed in a linear form and, on the other hand, the stability of the associated error system is demonstrated through utilizing a suitable Lyapunov function. To illustrate the applicability of the proposed schemes, several numerical simulations are performed in one- and two-spatial dimensions.
Collapse
|
43
|
A trajectory-based loss function to learn missing terms in bifurcating dynamical systems. Sci Rep 2021; 11:20394. [PMID: 34650131 PMCID: PMC8516982 DOI: 10.1038/s41598-021-99609-x] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2021] [Accepted: 09/21/2021] [Indexed: 11/15/2022] Open
Abstract
Missing terms in dynamical systems are a challenging problem for modeling. Recent developments in the combination of machine learning and dynamical system theory open possibilities for a solution. We show how physics-informed differential equations and machine learning—combined in the Universal Differential Equation (UDE) framework by Rackauckas et al.—can be modified to discover missing terms in systems that undergo sudden fundamental changes in their dynamical behavior called bifurcations. With this we enable the application of the UDE approach to a wider class of problems which are common in many real world applications. The choice of the loss function, which compares the training data trajectory in state space and the current estimated solution trajectory of the UDE to optimize the solution, plays a crucial role within this approach. The Mean Square Error as loss function contains the risk of a reconstruction which completely misses the dynamical behavior of the training data. By contrast, our suggested trajectory-based loss function which optimizes two largely independent components, the length and angle of state space vectors of the training data, performs reliable well in examples of systems from neuroscience, chemistry and biology showing Saddle-Node, Pitchfork, Hopf and Period-doubling bifurcations.
Collapse
|
44
|
Hauser MJB. Synchronisation of glycolytic activity in yeast cells. Curr Genet 2021; 68:69-81. [PMID: 34633492 DOI: 10.1007/s00294-021-01214-y] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/26/2021] [Revised: 09/26/2021] [Accepted: 09/27/2021] [Indexed: 11/28/2022]
Abstract
Glycolysis is the central metabolic pathway of almost every cell and organism. Under appropriate conditions, glycolytic oscillations may occur in individual cells as well as in entire cell populations or tissues. In many biological systems, glycolytic oscillations drive coherent oscillations of other metabolites, for instance in cardiomyocytes near anorexia, or in pancreas where they lead to a pulsatile release of insulin. Oscillations at the population or tissue level require the cells to synchronize their metabolism. We review the progress achieved in studying a model organism for glycolytic oscillations, namely yeast. Oscillations may occur on the level of individual cells as well as on the level of the cell population. In yeast, the cell-to-cell interaction is realized by diffusion-mediated intercellular communication via a messenger molecule. The present mini-review focuses on the synchronisation of glycolytic oscillations in yeast. Synchronisation is a quorum-sensing phenomenon because the collective oscillatory behaviour of a yeast cell population ceases when the cell density falls below a threshold. We review the question, under which conditions individual cells in a sparse population continue or cease to oscillate. Furthermore, we provide an overview of the pathway leading to the onset of synchronized oscillations. We also address the effects of spatial inhomogeneities (e.g., the formation of spatial clusters) on the collective dynamics, and also review the emergence of travelling waves of glycolytic activity. Finally, we briefly review the approaches used in numerical modelling of synchronized cell populations.
Collapse
Affiliation(s)
- Marcus J B Hauser
- Faculty of Natural Science, Otto-Von-Guericke-Universität Magdeburg, 39106, Magdeburg, Germany.
| |
Collapse
|
45
|
Kumar P, Gangopadhyay G. Nonequilibrium thermodynamics of glycolytic traveling wave: Benjamin-Feir instability. Phys Rev E 2021; 104:014221. [PMID: 34412344 DOI: 10.1103/physreve.104.014221] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/13/2021] [Accepted: 07/08/2021] [Indexed: 11/07/2022]
Abstract
Evolution of the nonequilibrium thermodynamic entities corresponding to dynamics of the Hopf instabilities and traveling waves at a nonequilibrium steady state of a spatially extended glycolysis model is assessed here by implementing an analytically tractable scheme incorporating a complex Ginzburg-Landau equation (CGLE). In the presence of self and cross diffusion, a more general amplitude equation exploiting the multiscale Krylov-Bogoliubov averaging method serves as an essential tool to reveal the various dynamical instability criteria, especially Benjamin-Feir (BF) instability, to estimate the corresponding nonlinear dispersion relation of the traveling wave pattern. The critical control parameter, wave-number selection criteria, and magnitude of the complex amplitude for traveling waves are modified by self- and cross-diffusion coefficients within the oscillatory regime, and their variabilities are exhibited against the amplitude equation. Unlike the traveling waves, a low-amplitude broad region appears for the Hopf instability in the concentration dynamics as the system phase passes through minima during its variation with the control parameter. The total entropy production rate of the uniform Hopf oscillation and glycolysis wave not only qualitatively reflects the global dynamics of concentrations of intermediate species but almost quantitatively. Despite the crucial role of diffusion in generating and shaping the traveling waves, the diffusive part of the entropy production rate has a negligible contribution to the system's total entropy production rate. The Hopf instability shows a more complex and colossal change in the energy profile of the open nonlinear system than in the traveling waves. A detailed analysis of BF instability shows a contrary nature of the semigrand Gibbs free energy with discrete and continuous wave numbers for the traveling wave. We hope the Hopf and traveling wave pattern around the BF instability in terms of energetics and dissipation will open up new applications of such dynamical phenomena.
Collapse
Affiliation(s)
- Premashis Kumar
- S. N. Bose National Centre For Basic Sciences, Block-JD, Sector-III, Salt Lake, Kolkata 700 106, India
| | - Gautam Gangopadhyay
- S. N. Bose National Centre For Basic Sciences, Block-JD, Sector-III, Salt Lake, Kolkata 700 106, India
| |
Collapse
|
46
|
Triggering avalanches by transverse perturbations in a rotating drum. Sci Rep 2021; 11:13936. [PMID: 34230549 PMCID: PMC8260778 DOI: 10.1038/s41598-021-93422-2] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/15/2021] [Accepted: 06/18/2021] [Indexed: 11/16/2022] Open
Abstract
We study the role of small-scale perturbations in the onset of avalanches in a rotating drum in the stick-slip regime. By vibrating the system along the axis of rotation with an amplitude orders of magnitude smaller than the particles’ diameter, we found that the order parameter that properly describes the system is the kinetic energy. We also show that, for high enough frequencies, the onset of the avalanche is determined by the amplitude of the oscillation, contrary to previous studies that showed that either acceleration or velocity was the governing parameter. Finally, we present a theoretical model that explains the transition between the continuous and discrete avalanche regimes as a supercritical Hopf bifurcation.
Collapse
|
47
|
Autocatalytic and oscillatory reaction networks that form guanidines and products of their cyclization. Nat Commun 2021; 12:2994. [PMID: 34016981 PMCID: PMC8138026 DOI: 10.1038/s41467-021-23206-9] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/12/2020] [Accepted: 03/26/2021] [Indexed: 11/08/2022] Open
Abstract
Autocatalytic and oscillatory networks of organic reactions are important for designing life-inspired materials and for better understanding the emergence of life on Earth; however, the diversity of the chemistries of these reactions is limited. In this work, we present the thiol-assisted formation of guanidines, which has a mechanism analogous to that of native chemical ligation. Using this reaction, we designed autocatalytic and oscillatory reaction networks that form substituted guanidines from thiouronium salts. The thiouronium salt-based oscillator show good stability of oscillations within a broad range of experimental conditions. By using nitrile-containing starting materials, we constructed an oscillator where the concentration of a bicyclic derivative of dihydropyrimidine oscillates. Moreover, the mixed thioester and thiouronium salt-based oscillator show unique responsiveness to chemical cues. The reactions developed in this work expand our toolbox for designing out-of-equilibrium chemical systems and link autocatalytic and oscillatory chemistry to the synthesis of guanidinium derivatives and the products of their transformations including analogs of nucleobases. So far, only a few chemical oscillators based on organic reactions have been developed. Here, the authors report both autocatalytic and oscillatory reaction networks that form substituted guanidines from thiouronium salts; when coupled to cascade cyclization, this reaction network produces oscillations in the production of pyrimidine-based heterocycles.
Collapse
|
48
|
The Emergence of the Bilateral Symmetry in Animals: A Review and a New Hypothesis. Symmetry (Basel) 2021. [DOI: 10.3390/sym13020261] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022] Open
Abstract
Most biological organisms exhibit different kinds of symmetry; an Animal (Metazoa), which is our Darwinist ancestor, has bilateral symmetry, and many plants exhibit rotational symmetry. It raises some questions: I. How can the evolution from an undifferentiated cell without bilateral symmetry to a complex biological organism with symmetry, which is based on asymmetric DNA and enzymes, lead to the bilateral symmetry? II. Is this evolution to an organism with bilateral symmetry obtained by other factors than DNA and enzymatic reactions? The existing literature about the evolution of the bilateral symmetry has been reviewed, and a new hypothesis has been formulated based on these reviews. The hypothesis is that the morphogenesis of biosystems is connected with the metabolism and that the oscillating kinetics in the Glycolysis have played a role in the polarity of the biological cells and in the establishment of the bilateral symmetry in Animals.
Collapse
|
49
|
Brechmann P, Rendall AD. Unbounded solutions of models for glycolysis. J Math Biol 2021; 82:1. [PMID: 33475794 PMCID: PMC7819955 DOI: 10.1007/s00285-021-01560-y] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/23/2020] [Revised: 09/23/2020] [Accepted: 11/18/2020] [Indexed: 11/26/2022]
Abstract
The Selkov oscillator, a simple description of glycolysis, is a system of two ordinary differential equations with mass action kinetics. In previous work the authors established several properties of the solutions of this system. In the present paper we extend this to prove that this system has solutions which diverge to infinity in an oscillatory manner at late times. This is done with the help of a Poincaré compactification of the system and a shooting argument. This system was originally derived from another system with Michaelis-Menten kinetics. A Poincaré compactification of the latter system is carried out and this is used to show that the Michaelis-Menten system, like that with mass action, has solutions which diverge to infinity in a monotone manner. It is also shown to admit subcritical Hopf bifurcations and thus unstable periodic solutions. We discuss to what extent the unbounded solutions cast doubt on the biological relevance of the Selkov oscillator and compare it with other models for the same biological system in the literature.
Collapse
Affiliation(s)
- Pia Brechmann
- Institut für Mathematik Johannes Gutenberg-Universität, Staudingerweg 9, 55099, Mainz, Germany
| | - Alan D Rendall
- Institut für Mathematik Johannes Gutenberg-Universität, Staudingerweg 9, 55099, Mainz, Germany.
| |
Collapse
|
50
|
Maria G. In silico Determination of Some Conditions Leading to Glycolytic Oscillations and Their Interference With Some Other Processes in E. coli Cells. Front Chem 2020; 8:526679. [PMID: 33195042 PMCID: PMC7655968 DOI: 10.3389/fchem.2020.526679] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2020] [Accepted: 09/23/2020] [Indexed: 01/05/2023] Open
Abstract
Autonomous oscillations of species levels in the glycolysis express the self-control of this essential cellular pathway belonging to the central carbon metabolism (CCM), and this phenomenon takes place in a large number of bacteria. Oscillations of glycolytic intermediates in living cells occur according to the environmental conditions and to the cell characteristics, especially the adenosine triphosphate (ATP) recovery system. Determining the conditions that lead to the occurrence and maintenance of the glycolytic oscillations can present immediate practical applications. Such a model-based analysis allows in silico (model-based) design of genetically modified microorganisms (GMO) with certain characteristics of interest for the biosynthesis industry, medicine, etc. Based on our kinetic model validated in previous works, this paper aims to in silico identify operating parameters and cell factors leading to the occurrence of stable glycolytic oscillations in the Escherichia coli cells. As long as most of the glycolytic intermediates are involved in various cellular metabolic pathways belonging to the CCM, evaluation of the dynamics and average level of its intermediates is of high importance for further applicative analyses. As an example, by using a lumped kinetic model for tryptophan (TRP) synthesis from literature, and its own kinetic model for the oscillatory glycolysis, this paper highlights the influence of glycolytic oscillations on the oscillatory TRP synthesis through the PEP (phosphoenolpyruvate) glycolytic node shared by the two oscillatory processes. The numerical analysis allows further TRP production maximization in a fed-batch bioreactor (FBR).
Collapse
Affiliation(s)
- Gheorghe Maria
- Department of Chemical and Biochemical Engineering, University POLITEHNICA of Bucharest, Bucharest, Romania.,Chemical Sciences Section, Romanian Academy, Bucharest, Romania
| |
Collapse
|