The key to a robust life system is to ensure that each cell population is maintained in an appropriate state. In this work, a mathematical model is used to investigate the control of the switching between the migrating and non-migrating states of the Bacillus subtilis cell population. In this case, the motile cells and matrix producers are the predominant cell types in the migrating cell population and non-migrating state, respectively, and can be suitably controlled according to the environmental conditions and cell density information. A minimal smooth model consisting of four ordinary differential equations is used as the mathematical model to control the B. subtilis cell types. Furthermore, the necessary and sufficient conditions for the hysteresis, which pertains to the change in the pheromone concentration, are clarified. In general, the hysteretic control of the cell state enables stable switching between the migrating and growth states of the B. subtilis cell population, thereby facilitating the biofilm life cycle. The results of corresponding culture experiments are examined, and the obtained corollaries are used to develop a model to input environmental conditions, especially, the external pH. On this basis, the environmental conditions are incorporated in a simulation model for the cell type control. In combination with a mathematical model of the cell population dynamics, a prediction model for colony growth involving multiple cell states, including concentric circular colonies of B. subtilis, can be established.
{"title":"Necessary and sufficient condition for hysteresis in the mathematical model of the cell type regulation of Bacillus subtilis.","authors":"Sohei Tasaki, Madoka Nakayama, Izumi Takagi, Jun-Ichi Wakita, Wataru Shoji","doi":"10.1007/s00285-025-02316-8","DOIUrl":"https://doi.org/10.1007/s00285-025-02316-8","url":null,"abstract":"<p><p>The key to a robust life system is to ensure that each cell population is maintained in an appropriate state. In this work, a mathematical model is used to investigate the control of the switching between the migrating and non-migrating states of the Bacillus subtilis cell population. In this case, the motile cells and matrix producers are the predominant cell types in the migrating cell population and non-migrating state, respectively, and can be suitably controlled according to the environmental conditions and cell density information. A minimal smooth model consisting of four ordinary differential equations is used as the mathematical model to control the B. subtilis cell types. Furthermore, the necessary and sufficient conditions for the hysteresis, which pertains to the change in the pheromone concentration, are clarified. In general, the hysteretic control of the cell state enables stable switching between the migrating and growth states of the B. subtilis cell population, thereby facilitating the biofilm life cycle. The results of corresponding culture experiments are examined, and the obtained corollaries are used to develop a model to input environmental conditions, especially, the external pH. On this basis, the environmental conditions are incorporated in a simulation model for the cell type control. In combination with a mathematical model of the cell population dynamics, a prediction model for colony growth involving multiple cell states, including concentric circular colonies of B. subtilis, can be established.</p>","PeriodicalId":50148,"journal":{"name":"Journal of Mathematical Biology","volume":"92 1","pages":"4"},"PeriodicalIF":2.3,"publicationDate":"2025-12-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145656230","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-12-01DOI: 10.1007/s00285-025-02317-7
Alexander Hermann, Tobias Köppl, Andreas Wagner, Arman Shojaei, Barbara Wohlmuth, Roland Aydin, Christian J Cyron, Roustem Miftahof
Cerebral blood flow regulation is critical for brain function, and its disruption is implicated in various neurological disorders. Many existing models do not fully capture the complex, multiscale interactions among neuronal activity, astrocytic signaling, and vascular dynamics, especially in key brainstem regions. In this work, we present a 3D-1D-0D multiscale computational framework for modeling the neuro-glial-vascular unit (NGVU) in the dorsal vagal complex (DVC). Our approach integrates a quadripartite synapse model, which captures the dynamic interactions among excitatory and inhibitory neurons, astrocytes, and vascular smooth muscle cells, with a hierarchical description of vascular dynamics that couples a three-dimensional microcirculatory network with a one-dimensional macrocirculatory representation and a zero-dimensional synaptic component. By linking neuronal spiking, astrocytic calcium and gliotransmitter signaling, and vascular tone regulation, our model reproduces key features of neurovascular regulation and elucidates the feedback loops that help maintain cerebral blood flow. Simulation results demonstrate that neurotransmitter release triggers astrocytic responses that modulate vessel radius, thereby influencing local oxygen and nutrient delivery. This integrated framework provides a robust and modular platform for future investigations into the pathophysiology of cerebral blood flow regulation and its role in autonomic control, including the regulation of gastric function.
{"title":"A 3D-1D-0D multiscale model of the neuro-glial-vascular unit for synaptic and vascular dynamics in the dorsal vagal complex.","authors":"Alexander Hermann, Tobias Köppl, Andreas Wagner, Arman Shojaei, Barbara Wohlmuth, Roland Aydin, Christian J Cyron, Roustem Miftahof","doi":"10.1007/s00285-025-02317-7","DOIUrl":"10.1007/s00285-025-02317-7","url":null,"abstract":"<p><p>Cerebral blood flow regulation is critical for brain function, and its disruption is implicated in various neurological disorders. Many existing models do not fully capture the complex, multiscale interactions among neuronal activity, astrocytic signaling, and vascular dynamics, especially in key brainstem regions. In this work, we present a 3D-1D-0D multiscale computational framework for modeling the neuro-glial-vascular unit (NGVU) in the dorsal vagal complex (DVC). Our approach integrates a quadripartite synapse model, which captures the dynamic interactions among excitatory and inhibitory neurons, astrocytes, and vascular smooth muscle cells, with a hierarchical description of vascular dynamics that couples a three-dimensional microcirculatory network with a one-dimensional macrocirculatory representation and a zero-dimensional synaptic component. By linking neuronal spiking, astrocytic calcium and gliotransmitter signaling, and vascular tone regulation, our model reproduces key features of neurovascular regulation and elucidates the feedback loops that help maintain cerebral blood flow. Simulation results demonstrate that neurotransmitter release triggers astrocytic responses that modulate vessel radius, thereby influencing local oxygen and nutrient delivery. This integrated framework provides a robust and modular platform for future investigations into the pathophysiology of cerebral blood flow regulation and its role in autonomic control, including the regulation of gastric function.</p>","PeriodicalId":50148,"journal":{"name":"Journal of Mathematical Biology","volume":"92 1","pages":"3"},"PeriodicalIF":2.3,"publicationDate":"2025-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC12669336/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145650054","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-12-01DOI: 10.1007/s00285-025-02298-7
Vasileios E Papageorgiou
Stochastic epidemic modeling has become increasingly crucial for assessing the severity of infectious diseases, attracting considerable attention in recent years. In this paper, we present three Markov-based epidemic models that incorporate demographic dynamics, including births, deaths, and migration. The inclusion of transition rates associated with these factors defines open-population systems, leading to a time-dependent transition pattern from the susceptible to the infectious phase. Notably, this work is the first to investigate epidemic models with time-varying population sizes within a Markovian framework. Furthermore, we introduce novel computational approaches for estimating stochastic features related to the number of secondary infections originating from an index case and the onset of a hazard (hitting) time associated with the number of susceptible cases in the system. Through extensive sensitivity analysis, we assess the impact of demographic dynamics on these descriptors and, consequently, on the severity of epidemic outbreaks. To validate the effectiveness of the introduced models, we utilize data from the 2022 mpox outbreak in Greece and examine the effect of interventions such as lockdowns on disease severity. This analysis helps health authorities identify optimal initiation periods and more effectively adjust the stringency of restrictive measures.
{"title":"A stochastic Markov-based modeling framework with demography.","authors":"Vasileios E Papageorgiou","doi":"10.1007/s00285-025-02298-7","DOIUrl":"https://doi.org/10.1007/s00285-025-02298-7","url":null,"abstract":"<p><p>Stochastic epidemic modeling has become increasingly crucial for assessing the severity of infectious diseases, attracting considerable attention in recent years. In this paper, we present three Markov-based epidemic models that incorporate demographic dynamics, including births, deaths, and migration. The inclusion of transition rates associated with these factors defines open-population systems, leading to a time-dependent transition pattern from the susceptible to the infectious phase. Notably, this work is the first to investigate epidemic models with time-varying population sizes within a Markovian framework. Furthermore, we introduce novel computational approaches for estimating stochastic features related to the number of secondary infections originating from an index case and the onset of a hazard (hitting) time associated with the number of susceptible cases in the system. Through extensive sensitivity analysis, we assess the impact of demographic dynamics on these descriptors and, consequently, on the severity of epidemic outbreaks. To validate the effectiveness of the introduced models, we utilize data from the 2022 mpox outbreak in Greece and examine the effect of interventions such as lockdowns on disease severity. This analysis helps health authorities identify optimal initiation periods and more effectively adjust the stringency of restrictive measures.</p>","PeriodicalId":50148,"journal":{"name":"Journal of Mathematical Biology","volume":"92 1","pages":"2"},"PeriodicalIF":2.3,"publicationDate":"2025-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145650165","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-11-28DOI: 10.1007/s00285-025-02320-y
Justin Eilertsen, Wylie Stroberg
The linear noise approximation (LNA) describes the random fluctuations from the mean-field concentrations of a chemical reaction network due to intrinsic noise. It is also used as a test probe to determine the accuracy of reduced formulations of the chemical master equation and to understand the relationship between timescale disparity and model reduction in stochastic environments. Although several reduced LNAs have been proposed, they have not been placed into a general theory concerning the accuracy of reduced LNAs derived from center manifold and singular perturbation theory. This has made it difficult to understand why certain reductions of the master or Langevin equations fail or succeed. In this work, we develop a deeper understanding of slow manifold projection in the linear noise regime by answering a straightforward but open question: In the presence of eigenvalue disparity, does the appropriate oblique projection of the LNA onto the slow eigenspace accurately approximate the first and second moments of complete LNA, and if not, why? Although most studies concentrate on the role of eigenvalue disparity arising from the drift matrix, we go further and examine the interplay between disparate 'drift" eigenvalues and the eigenvalues of the diffusion matrix, the latter of which may or may not be disparate. Furthermore, we place the previously established reductions of the LNA into a more general framework and formulate the necessary and sufficient conditions for the projected LNA to accurately approximate the first and second moments of the complete LNA.
{"title":"On the reduction of stochastic chemical reaction networks.","authors":"Justin Eilertsen, Wylie Stroberg","doi":"10.1007/s00285-025-02320-y","DOIUrl":"https://doi.org/10.1007/s00285-025-02320-y","url":null,"abstract":"<p><p>The linear noise approximation (LNA) describes the random fluctuations from the mean-field concentrations of a chemical reaction network due to intrinsic noise. It is also used as a test probe to determine the accuracy of reduced formulations of the chemical master equation and to understand the relationship between timescale disparity and model reduction in stochastic environments. Although several reduced LNAs have been proposed, they have not been placed into a general theory concerning the accuracy of reduced LNAs derived from center manifold and singular perturbation theory. This has made it difficult to understand why certain reductions of the master or Langevin equations fail or succeed. In this work, we develop a deeper understanding of slow manifold projection in the linear noise regime by answering a straightforward but open question: In the presence of eigenvalue disparity, does the appropriate oblique projection of the LNA onto the slow eigenspace accurately approximate the first and second moments of complete LNA, and if not, why? Although most studies concentrate on the role of eigenvalue disparity arising from the drift matrix, we go further and examine the interplay between disparate 'drift\" eigenvalues and the eigenvalues of the diffusion matrix, the latter of which may or may not be disparate. Furthermore, we place the previously established reductions of the LNA into a more general framework and formulate the necessary and sufficient conditions for the projected LNA to accurately approximate the first and second moments of the complete LNA.</p>","PeriodicalId":50148,"journal":{"name":"Journal of Mathematical Biology","volume":"92 1","pages":"1"},"PeriodicalIF":2.3,"publicationDate":"2025-11-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145642110","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-11-25DOI: 10.1007/s00285-025-02319-5
Mengqi Peng, Weihua Jiang
Crohn's disease (CD) is a recurrent chronic autoimmune disease, which is an inflammatory disease of the intestine with epithelial granulomas. The number of patients has been increasing significantly, and its pathogenesis and treatments are arousing hot discussions in the academic community. Taking into account the spatial heterogeneity of lesion distribution and the periodic recurrence, this paper uses a partial functional differential system with the free diffusion of bacteria and immunocytes and immune response latency to model the process of CD, based on the Lauffenburger-Kennedy bacterial infection model. In order to describe the spatial distribution and recurrence, we analyze the stability of the inflammation equilibrium state, and deduce the diffusion-driven Turing bifurcations and delay-driven Hopf bifurcations, drive the critical conditions for occurrence. Furthermore, through the analysis of Turing-Hopf bifurcations, the coupling effect of two factors is explored to obtain spatiotemporal patterns that better reflect clinical manifestations of CD. In addition, both theoretical and numerical results reveal that the motility is a necessary factor in the production of intestinal epithelial granulomas, while the immune response latency is an important factor in the recurrence. A small effective diffusion rate and a large time delay would lead to two spatially non-homogeneous steady states and a stable periodic solution, ultimately giving rise to a pair of stable spatially non-homogeneous periodic solutions through Turing-Hopf bifurcations. Our conclusions may provide some insights into the control mechanisms for Crohn's disease.
{"title":"Dynamical mechanisms of inflammatory spatial distribution and its association with recurrence in Crohn's disease.","authors":"Mengqi Peng, Weihua Jiang","doi":"10.1007/s00285-025-02319-5","DOIUrl":"10.1007/s00285-025-02319-5","url":null,"abstract":"<p><p>Crohn's disease (CD) is a recurrent chronic autoimmune disease, which is an inflammatory disease of the intestine with epithelial granulomas. The number of patients has been increasing significantly, and its pathogenesis and treatments are arousing hot discussions in the academic community. Taking into account the spatial heterogeneity of lesion distribution and the periodic recurrence, this paper uses a partial functional differential system with the free diffusion of bacteria and immunocytes and immune response latency to model the process of CD, based on the Lauffenburger-Kennedy bacterial infection model. In order to describe the spatial distribution and recurrence, we analyze the stability of the inflammation equilibrium state, and deduce the diffusion-driven Turing bifurcations and delay-driven Hopf bifurcations, drive the critical conditions for occurrence. Furthermore, through the analysis of Turing-Hopf bifurcations, the coupling effect of two factors is explored to obtain spatiotemporal patterns that better reflect clinical manifestations of CD. In addition, both theoretical and numerical results reveal that the motility is a necessary factor in the production of intestinal epithelial granulomas, while the immune response latency is an important factor in the recurrence. A small effective diffusion rate and a large time delay would lead to two spatially non-homogeneous steady states and a stable periodic solution, ultimately giving rise to a pair of stable spatially non-homogeneous periodic solutions through Turing-Hopf bifurcations. Our conclusions may provide some insights into the control mechanisms for Crohn's disease.</p>","PeriodicalId":50148,"journal":{"name":"Journal of Mathematical Biology","volume":"91 6","pages":"84"},"PeriodicalIF":2.3,"publicationDate":"2025-11-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145607337","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-11-25DOI: 10.1007/s00285-025-02315-9
Ning Wei, Yoichiro Mori
Cardiovascular disease continues to be the leading cause of death in the United States. A major contributing factor is cardiac arrhythmia, which results from irregular electrical activity in the heart. On a tissue level, cardiac conduction involves the spread of action potentials (AP) across the heart, enabling coordinated contraction of the myocardium. On a cellular level, the transmission of signals between cells is facilitated by low-resistance pathways formed by gap junctions (GJs). Recent experimental studies have sparked discussion on whether GJs play a dominant role in cell communication. Interestingly, research has revealed that GJ knockout mice can still demonstrate signal propagation in the heart, albeit more slowly and discontinuously, indicating the presence of an alternative mechanism for cardiac conduction. Unlike GJ-mediated propagation, ephaptic coupling (EpC) has emerged as a distinct form of electrical transmission, characterized by contactless electrochemical signaling across the narrow intercalated discs (IDs) between cardiomyocytes. Advancements in cardiac research have highlighted the crucial role of EpC in restoring conduction by increasing conduction velocity (CV), reducing conduction block (CB), and terminating reentry arrhythmias, particularly when GJs are impaired. However, most EpC studies are either numerical or experimental, while analytical studies on ephaptic conduction-an equally important aspect of understanding EpC-remain extremely limited. In this paper, we applied asymptotic theory to calculate the CV in the presence of weak EpC. To achieve this, we developed both continuous and discrete models to describe ephaptic conduction along a strand of cells. Ionic dynamics were modeled using the piecewise linear and cubic functions. The resulting system represents a bistable system with weak EpC. We calculated an expression for CV in the presence of weak EpC for both models, and validated our analytical results with numerical simulations. Additionally, we showed that under weak EpC, CV can increase if the distribution of INa is more prominent on the end membrane.
{"title":"Analytical Insights into Ephaptic Coupling and Its Effect on Conduction Velocity.","authors":"Ning Wei, Yoichiro Mori","doi":"10.1007/s00285-025-02315-9","DOIUrl":"10.1007/s00285-025-02315-9","url":null,"abstract":"<p><p>Cardiovascular disease continues to be the leading cause of death in the United States. A major contributing factor is cardiac arrhythmia, which results from irregular electrical activity in the heart. On a tissue level, cardiac conduction involves the spread of action potentials (AP) across the heart, enabling coordinated contraction of the myocardium. On a cellular level, the transmission of signals between cells is facilitated by low-resistance pathways formed by gap junctions (GJs). Recent experimental studies have sparked discussion on whether GJs play a dominant role in cell communication. Interestingly, research has revealed that GJ knockout mice can still demonstrate signal propagation in the heart, albeit more slowly and discontinuously, indicating the presence of an alternative mechanism for cardiac conduction. Unlike GJ-mediated propagation, ephaptic coupling (EpC) has emerged as a distinct form of electrical transmission, characterized by contactless electrochemical signaling across the narrow intercalated discs (IDs) between cardiomyocytes. Advancements in cardiac research have highlighted the crucial role of EpC in restoring conduction by increasing conduction velocity (CV), reducing conduction block (CB), and terminating reentry arrhythmias, particularly when GJs are impaired. However, most EpC studies are either numerical or experimental, while analytical studies on ephaptic conduction-an equally important aspect of understanding EpC-remain extremely limited. In this paper, we applied asymptotic theory to calculate the CV in the presence of weak EpC. To achieve this, we developed both continuous and discrete models to describe ephaptic conduction along a strand of cells. Ionic dynamics were modeled using the piecewise linear and cubic functions. The resulting system represents a bistable system with weak EpC. We calculated an expression for CV in the presence of weak EpC for both models, and validated our analytical results with numerical simulations. Additionally, we showed that under weak EpC, CV can increase if the distribution of INa is more prominent on the end membrane.</p>","PeriodicalId":50148,"journal":{"name":"Journal of Mathematical Biology","volume":"91 6","pages":"85"},"PeriodicalIF":2.3,"publicationDate":"2025-11-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC12647203/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145607271","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-11-19DOI: 10.1007/s00285-025-02313-x
R Gutiérrez, C Rojas-Jara
The identification of psychosocial risk and protective factors associated with substance use during adolescence is crucial for the prevention and treatment of addictive behaviors. However, in the mathematical modeling of addictions, these psychological and social dimensions are often addressed in isolation. By focusing on the immediate socialization environment of adolescents, comprising peers, parents and family, and school, and considering them as risk and protective factors, we mathematically model binge drinking behavior from interaction dynamics among individuals and their predisposition to consumption based on evolutionary and adaptive reasons. Our findings indicate that, in the face of negative peer influence, prosocial behaviors, as well as supportive attitudes from parents and family, and the school environment, have the potential to inhibit long-term binge drinking. This tendency is strengthened when protective situations are consistently promoted in response to risky scenarios, which is in line with both theoretical approaches to prevention and current public policy on drug consumption.
{"title":"The importance of the individual's micro-environment in drug consumption : A novel conceptual framework applied to the modeling of binge drinking.","authors":"R Gutiérrez, C Rojas-Jara","doi":"10.1007/s00285-025-02313-x","DOIUrl":"10.1007/s00285-025-02313-x","url":null,"abstract":"<p><p>The identification of psychosocial risk and protective factors associated with substance use during adolescence is crucial for the prevention and treatment of addictive behaviors. However, in the mathematical modeling of addictions, these psychological and social dimensions are often addressed in isolation. By focusing on the immediate socialization environment of adolescents, comprising peers, parents and family, and school, and considering them as risk and protective factors, we mathematically model binge drinking behavior from interaction dynamics among individuals and their predisposition to consumption based on evolutionary and adaptive reasons. Our findings indicate that, in the face of negative peer influence, prosocial behaviors, as well as supportive attitudes from parents and family, and the school environment, have the potential to inhibit long-term binge drinking. This tendency is strengthened when protective situations are consistently promoted in response to risky scenarios, which is in line with both theoretical approaches to prevention and current public policy on drug consumption.</p>","PeriodicalId":50148,"journal":{"name":"Journal of Mathematical Biology","volume":"91 6","pages":"83"},"PeriodicalIF":2.3,"publicationDate":"2025-11-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145551574","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-11-18DOI: 10.1007/s00285-025-02306-w
Yue Jiao, Xiaoxian Tang, Xiaowei Zeng
Zero-one biochemical reaction networks play key roles in cell signalling such as signalling pathways regulated by protein phosphorylation. Multistability of reaction networks is a crucial dynamics feature enabling decision-making in cells. It is well known that multistability can be lifted from a "subnetwork" (a network with less species and fewer reactions) to large networks. So, we aim to explore the multistability problem of small zero-one networks. In this work, we prove the following main results: 1. any zero-one network with a one-dimensional stoichiometric subspace admits at most one positive steady state (it must be stable), and all the one-dimensional zero-one networks can be classified according to if they indeed admit a stable positive steady state or not; 2. any two-dimensional zero-one network with up to three species either admits only degenerate positive steady states, or admits at most one positive steady state (it must be stable); 3. the smallest zero-one networks (here, by "smallest", we mean these networks contain species as few as possible) that admit nondegenerate multistationarity/multistability contain three species and five/six reactions, and they are three-dimensional. In these proofs, we use the theorems based on the Brouwer degree theory and the theory of real algebraic geometry. Moreover, applying the tools of computational real algebraic geometry, we provide a systematical way for detecting the networks that admit nondegenerate multistationarity/multistability.
{"title":"Multistability of small zero-one reaction networks.","authors":"Yue Jiao, Xiaoxian Tang, Xiaowei Zeng","doi":"10.1007/s00285-025-02306-w","DOIUrl":"10.1007/s00285-025-02306-w","url":null,"abstract":"<p><p>Zero-one biochemical reaction networks play key roles in cell signalling such as signalling pathways regulated by protein phosphorylation. Multistability of reaction networks is a crucial dynamics feature enabling decision-making in cells. It is well known that multistability can be lifted from a \"subnetwork\" (a network with less species and fewer reactions) to large networks. So, we aim to explore the multistability problem of small zero-one networks. In this work, we prove the following main results: 1. any zero-one network with a one-dimensional stoichiometric subspace admits at most one positive steady state (it must be stable), and all the one-dimensional zero-one networks can be classified according to if they indeed admit a stable positive steady state or not; 2. any two-dimensional zero-one network with up to three species either admits only degenerate positive steady states, or admits at most one positive steady state (it must be stable); 3. the smallest zero-one networks (here, by \"smallest\", we mean these networks contain species as few as possible) that admit nondegenerate multistationarity/multistability contain three species and five/six reactions, and they are three-dimensional. In these proofs, we use the theorems based on the Brouwer degree theory and the theory of real algebraic geometry. Moreover, applying the tools of computational real algebraic geometry, we provide a systematical way for detecting the networks that admit nondegenerate multistationarity/multistability.</p>","PeriodicalId":50148,"journal":{"name":"Journal of Mathematical Biology","volume":"91 6","pages":"82"},"PeriodicalIF":2.3,"publicationDate":"2025-11-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145551589","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
The closed fishing season policy plays a crucial role in fishery management by contributing to the restoration and protection of fishery resources, maintaining ecological balance and promoting sustainable development. The population dynamics of fish, particularly marine species, are highly complex. Under the combined effects of ecological mechanisms (such as predation, resource limitations, and competition), fish populations can exhibit multiple stable states. Overfishing increases the vulnerability of fish populations, making them prone to shift from a high-density stable state to a low-density one, and in some cases, leading to the risk of extinction. In this context, developing effective closed fishing season policies to ensure the sustainable development of fishery resources has become a pressing issue. In this paper, we propose a reaction-diffusion model consisting of two sub-equations with multiple stable states and a linear harvesting rate to describe the continuous switching between closed and open fishing seasons. We define a threshold value for the duration of the fishing ban . When , the trivial stable state is globally asymptotically stable. Uniqueness of periodic solutions is generally a mathematically challenging problem. However, employing the comparison theorem, we find that conditions on the uniqueness of periodic solutions to the associated ODE system are also applicable to our model. Specifically, under certain conditions, when , we provide sufficient conditions on the existence of a globally asymptotically stable periodic solution. Finally, we offer discussion and numerical simulations to illustrate our findings.
{"title":"Global dynamics of a reaction-diffusion switching model with multiple stable states and linear harvesting rate.","authors":"Yunfeng Liu, Huaqin Peng, Jianshe Yu, Yuming Chen, Zhiming Guo","doi":"10.1007/s00285-025-02312-y","DOIUrl":"10.1007/s00285-025-02312-y","url":null,"abstract":"<p><p>The closed fishing season policy plays a crucial role in fishery management by contributing to the restoration and protection of fishery resources, maintaining ecological balance and promoting sustainable development. The population dynamics of fish, particularly marine species, are highly complex. Under the combined effects of ecological mechanisms (such as predation, resource limitations, and competition), fish populations can exhibit multiple stable states. Overfishing increases the vulnerability of fish populations, making them prone to shift from a high-density stable state to a low-density one, and in some cases, leading to the risk of extinction. In this context, developing effective closed fishing season policies to ensure the sustainable development of fishery resources has become a pressing issue. In this paper, we propose a reaction-diffusion model consisting of two sub-equations with multiple stable states and a linear harvesting rate to describe the continuous switching between closed and open fishing seasons. We define a threshold value <math> <mmultiscripts><mover><mi>T</mi> <mo>¯</mo></mover> <mrow></mrow> <mrow><mrow></mrow> <mo>∗</mo></mrow> </mmultiscripts> </math> for the duration of the fishing ban <math><mover><mi>T</mi> <mo>¯</mo></mover> </math> . When <math> <mrow><mover><mi>T</mi> <mo>¯</mo></mover> <mo>≤</mo> <mmultiscripts><mover><mi>T</mi> <mo>¯</mo></mover> <mrow></mrow> <mrow><mrow></mrow> <mo>∗</mo></mrow> </mmultiscripts> </mrow> </math> , the trivial stable state is globally asymptotically stable. Uniqueness of periodic solutions is generally a mathematically challenging problem. However, employing the comparison theorem, we find that conditions on the uniqueness of periodic solutions to the associated ODE system are also applicable to our model. Specifically, under certain conditions, when <math> <mrow><mover><mi>T</mi> <mo>¯</mo></mover> <mo>></mo> <mmultiscripts><mover><mi>T</mi> <mo>¯</mo></mover> <mrow></mrow> <mrow><mrow></mrow> <mo>∗</mo></mrow> </mmultiscripts> </mrow> </math> , we provide sufficient conditions on the existence of a globally asymptotically stable periodic solution. Finally, we offer discussion and numerical simulations to illustrate our findings.</p>","PeriodicalId":50148,"journal":{"name":"Journal of Mathematical Biology","volume":"91 6","pages":"81"},"PeriodicalIF":2.3,"publicationDate":"2025-11-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145543738","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-11-12DOI: 10.1007/s00285-025-02308-8
Tom Britton, Andrea Pugliese
We consider a model for the spread of an influenza-like disease in which, between seasons, the virus makes a random genetic drift (reducing immunity) and obtains a new random transmissibility (closely related to ). Given the immunity status at the start of season k, i.e. the community distribution of years since last infection and their associated immunity levels, the outcome of the epidemic season k, characterized by the effective reproduction number and the fractions infected in the different immunity groups , is determined by the random genetic drift and transmissibility. It is shown that the community immunity status of consecutive seasons, is an ergodic Markov chain, which converges to a stationary distribution. More analytical progress is made for the case where immunity only lasts for one season: we then characterize the stationary distribution of the community fraction having partial immunity (from being infected last season) as well as the stationary distribution of , and the conditional distribution of given . The effective reproduction number is closely related to the initial exponential growth rate of the outbreak, a quantity which can be estimated early in the epidemic season. As a consequence, this conditional distribution may be used for predicting the final size of the epidemic based on its initial growth and immunity status.
我们考虑一种流感样疾病的传播模型,在季节之间,病毒进行随机遗传漂变(降低免疫力)并获得新的随机传播力(与r0密切相关)。鉴于k季节开始时的免疫状况,即自上次感染以来的年份群落分布及其相关的免疫水平,以有效繁殖数R e (k)和不同免疫组中感染的分数z (k)为特征的k流行季节的结果由随机遗传漂变和传播性决定。结果表明,连续季节的群体免疫状态是一条遍历马尔可夫链,并收敛于平稳分布。对于免疫仅持续一个季节的情况,我们进行了更多的分析进展:然后我们描述了具有部分免疫(从上一季节感染)的社区分数的平稳分布,以及(R e (k), z (k))的平稳分布,以及给定R e (k)的z (k)的条件分布。有效繁殖数R e (k)与暴发的初始指数增长率ρ (k)密切相关,ρ (k)可在流行季节早期估计。因此,该条件分布可用于根据其初始生长和免疫状态预测该流行病的最终规模。
{"title":"A multi-season epidemic model with random genetic drift and transmissibility.","authors":"Tom Britton, Andrea Pugliese","doi":"10.1007/s00285-025-02308-8","DOIUrl":"10.1007/s00285-025-02308-8","url":null,"abstract":"<p><p>We consider a model for the spread of an influenza-like disease in which, between seasons, the virus makes a random genetic drift (reducing immunity) and obtains a new random transmissibility (closely related to <math><msub><mi>R</mi> <mn>0</mn></msub> </math> ). Given the immunity status at the start of season k, i.e. the community distribution of years since last infection and their associated immunity levels, the outcome of the epidemic season k, characterized by the effective reproduction number <math><msubsup><mi>R</mi> <mi>e</mi> <mrow><mo>(</mo> <mi>k</mi> <mo>)</mo></mrow> </msubsup> </math> and the fractions infected in the different immunity groups <math> <msup><mrow><mi>z</mi></mrow> <mrow><mo>(</mo> <mi>k</mi> <mo>)</mo></mrow> </msup> </math> , is determined by the random genetic drift and transmissibility. It is shown that the community immunity status of consecutive seasons, is an ergodic Markov chain, which converges to a stationary distribution. More analytical progress is made for the case where immunity only lasts for one season: we then characterize the stationary distribution of the community fraction having partial immunity (from being infected last season) as well as the stationary distribution of <math><mrow><mo>(</mo> <msubsup><mi>R</mi> <mi>e</mi> <mrow><mo>(</mo> <mi>k</mi> <mo>)</mo></mrow> </msubsup> <mo>,</mo> <msup><mi>z</mi> <mrow><mo>(</mo> <mi>k</mi> <mo>)</mo></mrow> </msup> <mo>)</mo></mrow> </math> , and the conditional distribution of <math><msup><mi>z</mi> <mrow><mo>(</mo> <mi>k</mi> <mo>)</mo></mrow> </msup> </math> given <math><msubsup><mi>R</mi> <mi>e</mi> <mrow><mo>(</mo> <mi>k</mi> <mo>)</mo></mrow> </msubsup> </math> . The effective reproduction number <math><msubsup><mi>R</mi> <mi>e</mi> <mrow><mo>(</mo> <mi>k</mi> <mo>)</mo></mrow> </msubsup> </math> is closely related to the initial exponential growth rate <math><msup><mi>ρ</mi> <mrow><mo>(</mo> <mi>k</mi> <mo>)</mo></mrow> </msup> </math> of the outbreak, a quantity which can be estimated early in the epidemic season. As a consequence, this conditional distribution may be used for predicting the final size of the epidemic based on its initial growth and immunity status.</p>","PeriodicalId":50148,"journal":{"name":"Journal of Mathematical Biology","volume":"91 6","pages":"80"},"PeriodicalIF":2.3,"publicationDate":"2025-11-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC12612026/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145497418","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}