Pub Date : 2025-09-09DOI: 10.1016/j.jtbi.2025.112267
Ana S. de Pereda, Jihyun Park, Lily S. Cheung
Transporters play key roles in regulating the movement of molecules into and out of cells. Uniporters, the simplest class of transporters, use facilitated diffusion to translocate molecules across membranes down their concentration gradient. This process can be affected by the presence of additional substrates in the intra- and extracellular environment, which can either increase the net transport rate of a molecule via trans acceleration or decrease it via competitive inhibition. In this study, we derived mathematical models to describe the net transport rate of uniporters in the presence of multiple extracellular substrates or inhibitors. Analyses of these models identified four possible states for the system when two substrates are present, with two states leading to trans acceleration and the other two states resulting in inhibition. Finally, we found that the relation between kinetic constants that controls the fraction of transporters in the inward-facing open state is responsible for these behaviors. Our theoretical results provide a mathematical framework for understanding the dynamic response of uniporters in the presence of multiple substrates and inhibitors, which could have implications for various processes, from nutrient utilization to metabolic engineering.
{"title":"A kinetic study of multi-substrate uniporters","authors":"Ana S. de Pereda, Jihyun Park, Lily S. Cheung","doi":"10.1016/j.jtbi.2025.112267","DOIUrl":"10.1016/j.jtbi.2025.112267","url":null,"abstract":"<div><div>Transporters play key roles in regulating the movement of molecules into and out of cells. Uniporters, the simplest class of transporters, use facilitated diffusion to translocate molecules across membranes down their concentration gradient. This process can be affected by the presence of additional substrates in the intra- and extracellular environment, which can either increase the net transport rate of a molecule via trans acceleration or decrease it via competitive inhibition. In this study, we derived mathematical models to describe the net transport rate of uniporters in the presence of multiple extracellular substrates or inhibitors. Analyses of these models identified four possible states for the system when two substrates are present, with two states leading to trans acceleration and the other two states resulting in inhibition. Finally, we found that the relation between kinetic constants that controls the fraction of transporters in the inward-facing open state is responsible for these behaviors. Our theoretical results provide a mathematical framework for understanding the dynamic response of uniporters in the presence of multiple substrates and inhibitors, which could have implications for various processes, from nutrient utilization to metabolic engineering.</div></div>","PeriodicalId":54763,"journal":{"name":"Journal of Theoretical Biology","volume":"616 ","pages":"Article 112267"},"PeriodicalIF":2.0,"publicationDate":"2025-09-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145042514","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}
Several computational models are available for representing the gene expression process, with each having their advantages and disadvantages. Phenomenological models are widely used as they make appropriate simplifications that aim to find a middle ground between accuracy and complexity. The existing phenomenological models compete in terms of how the transcription initiation process is approximated, to achieve high accuracy while having the lowest complexity possible. However, most current models still suffer from high parameter complexity in the case of complex promoters. Herein, we formally derive a phenomenological approach to model RNA polymerase recruitment, stating approximations on cooperativity between transcription factors that are applicable to promoters requiring multifactorial input, which reduces parameter complexity. We then apply this method to biologically relevant networks of varying complexities to show that the approximations improved predictive ability compared to existing models. In summary, our reduced parameter model (RPM) had lower complexity while maintaining high accuracy, which leads to better scalability for complex networks.
{"title":"Phenomenological modeling of gene transcription by approximating cooperativity of transcription factors improves prediction and reduces complexity in gene regulatory network models","authors":"Thiruvickraman Jothiprakasam, Siddharth Jhunjhunwala","doi":"10.1016/j.jtbi.2025.112264","DOIUrl":"10.1016/j.jtbi.2025.112264","url":null,"abstract":"<div><div>Several computational models are available for representing the gene expression process, with each having their advantages and disadvantages. Phenomenological models are widely used as they make appropriate simplifications that aim to find a middle ground between accuracy and complexity. The existing phenomenological models compete in terms of how the transcription initiation process is approximated, to achieve high accuracy while having the lowest complexity possible. However, most current models still suffer from high parameter complexity in the case of complex promoters. Herein, we formally derive a phenomenological approach to model RNA polymerase recruitment, stating approximations on cooperativity between transcription factors that are applicable to promoters requiring multifactorial input, which reduces parameter complexity. We then apply this method to biologically relevant networks of varying complexities to show that the approximations improved predictive ability compared to existing models. In summary, our reduced parameter model (RPM) had lower complexity while maintaining high accuracy, which leads to better scalability for complex networks.</div></div>","PeriodicalId":54763,"journal":{"name":"Journal of Theoretical Biology","volume":"616 ","pages":"Article 112264"},"PeriodicalIF":2.0,"publicationDate":"2025-09-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145030705","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-09-04DOI: 10.1016/j.jtbi.2025.112263
Shilian Xu , Maoxuan Liu
Chimeric antigen receptor (CAR)-macrophage therapy is a promising approach for tumour treatment due to antigen-specific phagocytosis and tumour clearance. However, the precise impact of tumour burden, dose and dosing regimens on therapeutic outcomes remains poorly understood. We developed ordinary differential equation (ODE) mathematical modelling and utilised parameter inference to analyse in vitro FACS-based phagocytosis assay data testing CD19-positive Raji tumour cell against CAR-macrophage, and revealed that phagocytosing efficiency of CAR-macrophage increases but saturates as both Raji cell and CAR-macrophage concentrations increase. This interaction resulted in bistable Raji cell kinetics; specifically, within a particular range of CAR-macrophage concentration, low tumour burdens are effectively inhibited, while high tumour burdens remain refractory. Furthermore, our model predicted that CAR-macrophage dosages typically suggested by current clinical trials yield favourable therapeutic outcomes only when tumour burden is low. For split CAR-macrophage infusion with fixed total dosage, the first infusion with high CAR-macrophage dose delivers superior treatment outcomes. Finally, we identified alternative infusion regimens: five billion cells administered monthly for three months, or seven billion cells every two months for six months, can efficiently suppress Raji cell replication irrespective of tumour burden. Our findings highlight CAR-macrophage therapeutic outcomes are strongly influenced by both tumour burden and different dosing regimens. This work underscores that reducing tumour burden, increasing CAR-macrophage dose in the first infusion and prolonging CAR-macrophage persistence are key strategies for achieving durable responses.
{"title":"Mathematical model suggests current CAR-macrophage dosage is efficient to low pre-infusion tumour burden but refractory to high tumour burden","authors":"Shilian Xu , Maoxuan Liu","doi":"10.1016/j.jtbi.2025.112263","DOIUrl":"10.1016/j.jtbi.2025.112263","url":null,"abstract":"<div><div>Chimeric antigen receptor (CAR)-macrophage therapy is a promising approach for tumour treatment due to antigen-specific phagocytosis and tumour clearance. However, the precise impact of tumour burden, dose and dosing regimens on therapeutic outcomes remains poorly understood. We developed ordinary differential equation (ODE) mathematical modelling and utilised parameter inference to analyse <em>in vitro</em> FACS-based phagocytosis assay data testing CD19-positive Raji tumour cell against CAR-macrophage, and revealed that phagocytosing efficiency of CAR-macrophage increases but saturates as both Raji cell and CAR-macrophage concentrations increase. This interaction resulted in bistable Raji cell kinetics; specifically, within a particular range of CAR-macrophage concentration, low tumour burdens are effectively inhibited, while high tumour burdens remain refractory. Furthermore, our model predicted that CAR-macrophage dosages typically suggested by current clinical trials yield favourable therapeutic outcomes only when tumour burden is low. For split CAR-macrophage infusion with fixed total dosage, the first infusion with high CAR-macrophage dose delivers superior treatment outcomes. Finally, we identified alternative infusion regimens: five billion cells administered monthly for three months, or seven billion cells every two months for six months, can efficiently suppress Raji cell replication irrespective of tumour burden. Our findings highlight CAR-macrophage therapeutic outcomes are strongly influenced by both tumour burden and different dosing regimens. This work underscores that reducing tumour burden, increasing CAR-macrophage dose in the first infusion and prolonging CAR-macrophage persistence are key strategies for achieving durable responses.</div></div>","PeriodicalId":54763,"journal":{"name":"Journal of Theoretical Biology","volume":"616 ","pages":"Article 112263"},"PeriodicalIF":2.0,"publicationDate":"2025-09-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145008508","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-09-02DOI: 10.1016/j.jtbi.2025.112246
Mingqi He , Sophie Hautphenne , Yao-ban Chan
Phylogenetic trees describe the relationships between species in the evolutionary process, and provide information about the rates of diversification. To understand the mechanisms behind macroevolution, we consider a class of multitype branching processes called Markovian binary trees (MBTs). MBTs allow for trait-based variation in diversification rates, and provide a flexible and realistic probabilistic model for phylogenetic trees. We develop an approximate Bayesian computation (ABC) scheme to infer the rates of MBT parameters by exploiting the information in the shapes of phylogenetic trees. We evaluate the accuracy of this inference method using simulation studies, and find that our method is able to detect variation in the diversification rates, with accuracy comparable to, and generally better than, likelihood-based methods. In an application to a real-life phylogeny of squamata, we reinforce conclusions drawn from earlier studies, in particular supporting the existence of ovi-/viviparity transitions in both directions. Our method demonstrates the potential for more complex models of evolution to be employed in phylogenetic inference, in conjunction with likelihood-free schemes.
{"title":"Approximate Bayesian computation for Markovian binary trees in phylogenetics","authors":"Mingqi He , Sophie Hautphenne , Yao-ban Chan","doi":"10.1016/j.jtbi.2025.112246","DOIUrl":"10.1016/j.jtbi.2025.112246","url":null,"abstract":"<div><div>Phylogenetic trees describe the relationships between species in the evolutionary process, and provide information about the rates of diversification. To understand the mechanisms behind macroevolution, we consider a class of multitype branching processes called Markovian binary trees (MBTs). MBTs allow for trait-based variation in diversification rates, and provide a flexible and realistic probabilistic model for phylogenetic trees. We develop an approximate Bayesian computation (ABC) scheme to infer the rates of MBT parameters by exploiting the information in the shapes of phylogenetic trees. We evaluate the accuracy of this inference method using simulation studies, and find that our method is able to detect variation in the diversification rates, with accuracy comparable to, and generally better than, likelihood-based methods. In an application to a real-life phylogeny of squamata, we reinforce conclusions drawn from earlier studies, in particular supporting the existence of ovi-/viviparity transitions in both directions. Our method demonstrates the potential for more complex models of evolution to be employed in phylogenetic inference, in conjunction with likelihood-free schemes.</div></div>","PeriodicalId":54763,"journal":{"name":"Journal of Theoretical Biology","volume":"616 ","pages":"Article 112246"},"PeriodicalIF":2.0,"publicationDate":"2025-09-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145001904","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-08-31DOI: 10.1016/j.jtbi.2025.112255
Andrés Hernández-Rivera , Pablo Velarde , Ascensión Zafra-Cabeza , José M. Maestre
Stochastic Model Predictive Control (SMPC) is an effective decision-making method in applications where uncertainties play a significant role. This work introduces a non-linear formulation of SMPC specifically designed for cancer therapy. The proposed method considers the stochastic nature of tumor growth, non-linear dynamics, and a potential side effect of the treatment. Through one-year simulations, the results showcase the effectiveness of this strategy in controlling drug dosing.
{"title":"Drug dosing for cancer therapy: A stochastic model predictive control perspective","authors":"Andrés Hernández-Rivera , Pablo Velarde , Ascensión Zafra-Cabeza , José M. Maestre","doi":"10.1016/j.jtbi.2025.112255","DOIUrl":"10.1016/j.jtbi.2025.112255","url":null,"abstract":"<div><div>Stochastic Model Predictive Control (SMPC) is an effective decision-making method in applications where uncertainties play a significant role. This work introduces a non-linear formulation of SMPC specifically designed for cancer therapy. The proposed method considers the stochastic nature of tumor growth, non-linear dynamics, and a potential side effect of the treatment. Through one-year simulations, the results showcase the effectiveness of this strategy in controlling drug dosing.</div></div>","PeriodicalId":54763,"journal":{"name":"Journal of Theoretical Biology","volume":"615 ","pages":"Article 112255"},"PeriodicalIF":2.0,"publicationDate":"2025-08-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144979377","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-08-30DOI: 10.1016/j.jtbi.2025.112254
Anna C. Nelson , Scott A. McKinley , Melissa M. Rolls , Maria-Veronica Ciocanel
Microtubules (MTs) are dynamic protein filaments essential for intracellular organization and transport, particularly in long-lived cells such as neurons. The plus and minus ends of neuronal MTs switch between growth and shrinking phases, and the nucleation of new filaments is believed to be regulated in both healthy and injury conditions. We propose stochastic and deterministic mathematical models to investigate the impact of filament nucleation and length-regulation mechanisms on emergent properties such as MT lengths and numbers in living cells. We expand our stochastic continuous-time Markov chain model of filament dynamics to incorporate MT nucleation and capture realistic stochastic fluctuations in MT numbers and tubulin availability. We also propose a simplified partial differential equation (PDE) model, which allows for tractable analytical investigation into steady-state MT distributions under different nucleation and length-regulating mechanisms. We find that the stochastic and PDE modeling approaches show good agreement in MT length distributions, and that both MT nucleation and the catastrophe rate of large-length MTs regulate MT length distributions. In both frameworks, multiple mechanistic combinations achieve the same average MT length. The models proposed can predict parameter regimes where the system is scarce in tubulin, the building block of MTs, and suggest that low filament nucleation regimes are characterized by high variation in MT lengths, while high nucleation regimes drive high variation in MT numbers. These mathematical frameworks have the potential to improve our understanding of MT regulation in both healthy and injured neurons.
{"title":"Emergent microtubule properties in a model of filament turnover and nucleation","authors":"Anna C. Nelson , Scott A. McKinley , Melissa M. Rolls , Maria-Veronica Ciocanel","doi":"10.1016/j.jtbi.2025.112254","DOIUrl":"10.1016/j.jtbi.2025.112254","url":null,"abstract":"<div><div>Microtubules (MTs) are dynamic protein filaments essential for intracellular organization and transport, particularly in long-lived cells such as neurons. The plus and minus ends of neuronal MTs switch between growth and shrinking phases, and the nucleation of new filaments is believed to be regulated in both healthy and injury conditions. We propose stochastic and deterministic mathematical models to investigate the impact of filament nucleation and length-regulation mechanisms on emergent properties such as MT lengths and numbers in living cells. We expand our stochastic continuous-time Markov chain model of filament dynamics to incorporate MT nucleation and capture realistic stochastic fluctuations in MT numbers and tubulin availability. We also propose a simplified partial differential equation (PDE) model, which allows for tractable analytical investigation into steady-state MT distributions under different nucleation and length-regulating mechanisms. We find that the stochastic and PDE modeling approaches show good agreement in MT length distributions, and that both MT nucleation and the catastrophe rate of large-length MTs regulate MT length distributions. In both frameworks, multiple mechanistic combinations achieve the same average MT length. The models proposed can predict parameter regimes where the system is scarce in tubulin, the building block of MTs, and suggest that low filament nucleation regimes are characterized by high variation in MT lengths, while high nucleation regimes drive high variation in MT numbers. These mathematical frameworks have the potential to improve our understanding of MT regulation in both healthy and injured neurons.</div></div>","PeriodicalId":54763,"journal":{"name":"Journal of Theoretical Biology","volume":"616 ","pages":"Article 112254"},"PeriodicalIF":2.0,"publicationDate":"2025-08-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144979356","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-08-30DOI: 10.1016/j.jtbi.2025.112256
F.E. Cornes , R.H. Barriga Rubio , M. Otero
We present an extension of a previously developed stochastic, stage-structured model of Dalbulus maidis (corn leafhopper), an important pest and vector in maize crops. The extended model introduces nonlinear density-dependent regulation on the nymphal stage, mediated by a carrying capacity that dynamically depends on the leaf area of maize plants. Both insect and host-plant dynamics are explicitly modeled, but the interaction is asymmetric, as the plant is not affected by the insect in the present formulation. Our main objective is to explore how the interplay between temperature-driven development and host-plant dynamics shapes the long-term behavior of the insect population, leading to either extinction or persistence. Using simulations parameterized with laboratory and field data, we analyze how temperature and maize development affect insect dynamics, and assess whether the model can reproduce observed abundance patterns under realistic conditions. This modeling framework provides a biologically grounded and flexible basis for future extensions, including pathogen transmission and bidirectional feedback between the maize and the insect.
{"title":"Extinction and persistence in a temperature-driven, stage-structured stochastic model of Dalbulus maidis dynamics with nonlinear density-dependent regulation","authors":"F.E. Cornes , R.H. Barriga Rubio , M. Otero","doi":"10.1016/j.jtbi.2025.112256","DOIUrl":"10.1016/j.jtbi.2025.112256","url":null,"abstract":"<div><div>We present an extension of a previously developed stochastic, stage-structured model of <em>Dalbulus maidis</em> (corn leafhopper), an important pest and vector in maize crops. The extended model introduces nonlinear density-dependent regulation on the nymphal stage, mediated by a carrying capacity that dynamically depends on the leaf area of maize plants. Both insect and host-plant dynamics are explicitly modeled, but the interaction is asymmetric, as the plant is not affected by the insect in the present formulation. Our main objective is to explore how the interplay between temperature-driven development and host-plant dynamics shapes the long-term behavior of the insect population, leading to either extinction or persistence. Using simulations parameterized with laboratory and field data, we analyze how temperature and maize development affect insect dynamics, and assess whether the model can reproduce observed abundance patterns under realistic conditions. This modeling framework provides a biologically grounded and flexible basis for future extensions, including pathogen transmission and bidirectional feedback between the maize and the insect.</div></div>","PeriodicalId":54763,"journal":{"name":"Journal of Theoretical Biology","volume":"615 ","pages":"Article 112256"},"PeriodicalIF":2.0,"publicationDate":"2025-08-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144979347","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-08-29DOI: 10.1016/j.jtbi.2025.112250
Xuesong Bai, Thomas G. Fai
It has been observed that the growth of the nucleus and the cytoplasm is coordinated during cell growth, resulting in a nearly constant nuclear-to-cell volume ratio (N/C) throughout the cell cycle. Previous studies have shown that the N/C ratio is determined by the ratio between the number of proteins in the nucleus and the total number of proteins in the cell. These observations suggest the importance of the nucleocytoplasmic transport process in nuclear size by regulating protein concentrations in the nucleus and cytoplasm. This paper combines a biophysical model of Ran-mediated nucleocytoplasmic transport and a simple cell growth model to provide insights into several key aspects of the N/C ratio homeostasis in growing cells. Our model shows that the permeability of the nuclear envelope needs to grow in line with the cell to maintain a nearly constant N/C ratio, that several parameters involved in the nucleocytoplasmic transport mechanism and gene translation significantly affect the N/C ratio, and that Ran may potentially compensate for the lack of NTF2 in the nucleocytoplasmic transport mechanism to maintain a viable N/C ratio. However, this compensation is possible only if RanGDP is allowed to translocate through the nuclear envelope independently of NTF2.
{"title":"Mathematical model of nucleocytoplasmic transport and nuclear-to-cell ratio in a growing cell","authors":"Xuesong Bai, Thomas G. Fai","doi":"10.1016/j.jtbi.2025.112250","DOIUrl":"10.1016/j.jtbi.2025.112250","url":null,"abstract":"<div><div>It has been observed that the growth of the nucleus and the cytoplasm is coordinated during cell growth, resulting in a nearly constant nuclear-to-cell volume ratio (N/C) throughout the cell cycle. Previous studies have shown that the N/C ratio is determined by the ratio between the number of proteins in the nucleus and the total number of proteins in the cell. These observations suggest the importance of the nucleocytoplasmic transport process in nuclear size by regulating protein concentrations in the nucleus and cytoplasm. This paper combines a biophysical model of Ran-mediated nucleocytoplasmic transport and a simple cell growth model to provide insights into several key aspects of the N/C ratio homeostasis in growing cells. Our model shows that the permeability of the nuclear envelope needs to grow in line with the cell to maintain a nearly constant N/C ratio, that several parameters involved in the nucleocytoplasmic transport mechanism and gene translation significantly affect the N/C ratio, and that Ran may potentially compensate for the lack of NTF2 in the nucleocytoplasmic transport mechanism to maintain a viable N/C ratio. However, this compensation is possible only if RanGDP is allowed to translocate through the nuclear envelope independently of NTF2.</div></div>","PeriodicalId":54763,"journal":{"name":"Journal of Theoretical Biology","volume":"615 ","pages":"Article 112250"},"PeriodicalIF":2.0,"publicationDate":"2025-08-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144979444","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-08-28DOI: 10.1016/j.jtbi.2025.112253
Gideon A. Ngwa , Bime M. Ghakanyuy , Miranda I. Teboh-Ewungkem , Jacek Banasiak
A deterministic nonlinear ordinary differential equation model for mosquito dynamics in which the mosquitoes can quest for blood either within a human population or within non-human/vertebrate populations is derived and studied. The model captures both the mosquito’s aquatic and terrestrial forms and includes a mechanism to investigate the impact of mating on mosquito dynamics. The model uses a restricted form of homogeneous mixing based on the idea that the mosquito has a blood-feeding habit determined by its blood-feeding preferences and its gonotrophic cycle. This characterisation allows us to compartmentalise the total mosquito population into distinct compartments according to the spatial location of the mosquito (breeding site, resting places and questing places) as well as blood-fed status. Issues of overcrowding and intraspecific competition both within the aquatic and the terrestrial stages of the mosquito’s life forms are addressed and considered in the model. Results show that the inclusion of mating induces bistability, a phenomenon whereby locally stable trivial and non-trivial equilibria co-exist with an unstable non-zero equilibrium. The local nature of the stable equilibria is demonstrated by numerically showing that the long-term state of the system is sensitive to initial conditions. The bistability state is analogous to the phenomenon of the Allee effect that has been reported in population biology. The model’s results, including the derivation of the threshold parameter of the system, are comprehensively tested via numerical simulations. The output of our model has direct application to mosquito control strategies, for it clearly shows key points in the mosquito’s developmental pathway that can be targeted for control purposes.
{"title":"Mating versus alternative blood sources as determinants to mosquito abundance and population resilience","authors":"Gideon A. Ngwa , Bime M. Ghakanyuy , Miranda I. Teboh-Ewungkem , Jacek Banasiak","doi":"10.1016/j.jtbi.2025.112253","DOIUrl":"10.1016/j.jtbi.2025.112253","url":null,"abstract":"<div><div>A deterministic nonlinear ordinary differential equation model for mosquito dynamics in which the mosquitoes can quest for blood either within a human population or within non-human/vertebrate populations is derived and studied. The model captures both the mosquito’s aquatic and terrestrial forms and includes a mechanism to investigate the impact of mating on mosquito dynamics. The model uses a restricted form of homogeneous mixing based on the idea that the mosquito has a blood-feeding habit determined by its blood-feeding preferences and its gonotrophic cycle. This characterisation allows us to compartmentalise the total mosquito population into distinct compartments according to the spatial location of the mosquito (breeding site, resting places and questing places) as well as blood-fed status. Issues of overcrowding and intraspecific competition both within the aquatic and the terrestrial stages of the mosquito’s life forms are addressed and considered in the model. Results show that the inclusion of mating induces bistability, a phenomenon whereby locally stable trivial and non-trivial equilibria co-exist with an unstable non-zero equilibrium. The local nature of the stable equilibria is demonstrated by numerically showing that the long-term state of the system is sensitive to initial conditions. The bistability state is analogous to the phenomenon of the Allee effect that has been reported in population biology. The model’s results, including the derivation of the threshold parameter of the system, are comprehensively tested via numerical simulations. The output of our model has direct application to mosquito control strategies, for it clearly shows key points in the mosquito’s developmental pathway that can be targeted for control purposes.</div></div>","PeriodicalId":54763,"journal":{"name":"Journal of Theoretical Biology","volume":"615 ","pages":"Article 112253"},"PeriodicalIF":2.0,"publicationDate":"2025-08-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144979364","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-08-24DOI: 10.1016/j.jtbi.2025.112248
Punya Alahakoon , Peter G. Taylor , James M. McCaw
COVID-19, caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has affected millions of people around the globe. We studied the spread of SARS-CoV-2 across six rural counties in North and South Dakota in the United States. The study period was from early March 2020 to mid-June 2021, during which non-pharmaceutical interventions (NPIs) were in place. The end of the study period coincided with the emergence of the Delta variant in the United States. We modelled the transmission dynamics in each county using a stochastic compartmental model and analysed the data within a Bayesian hierarchical statistical framework. We estimated key epidemiological and surveillance parameters including the reproduction number and reporting probability. We conducted a series of counterfactual analyses in which NPIs were lifted earlier and by varying degrees, modelled as an increase in the transmission rate. Under this range of plausible alternative responses, increases in case counts varied from negligible to substantial, underscoring the importance of timely public health measures and compliance with them. From a methodological perspective, our study demonstrates that despite the inherent high variability in epidemic behaviour in small rural communities, the combination of stochastic modelling and application of Bayesian hierarchical analyses enables the estimation of key epidemiological and surveillance parameters and consideration of the potential impact of alternative public health measures in small low population density communities.
{"title":"Stochastic modelling of early-stage COVID-19 epidemic dynamics in rural communities in the United States","authors":"Punya Alahakoon , Peter G. Taylor , James M. McCaw","doi":"10.1016/j.jtbi.2025.112248","DOIUrl":"10.1016/j.jtbi.2025.112248","url":null,"abstract":"<div><div>COVID-19, caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has affected millions of people around the globe. We studied the spread of SARS-CoV-2 across six rural counties in North and South Dakota in the United States. The study period was from early March 2020 to mid-June 2021, during which non-pharmaceutical interventions (NPIs) were in place. The end of the study period coincided with the emergence of the Delta variant in the United States. We modelled the transmission dynamics in each county using a stochastic compartmental model and analysed the data within a Bayesian hierarchical statistical framework. We estimated key epidemiological and surveillance parameters including the reproduction number and reporting probability. We conducted a series of counterfactual analyses in which NPIs were lifted earlier and by varying degrees, modelled as an increase in the transmission rate. Under this range of plausible alternative responses, increases in case counts varied from negligible to substantial, underscoring the importance of timely public health measures and compliance with them. From a methodological perspective, our study demonstrates that despite the inherent high variability in epidemic behaviour in small rural communities, the combination of stochastic modelling and application of Bayesian hierarchical analyses enables the estimation of key epidemiological and surveillance parameters and consideration of the potential impact of alternative public health measures in small low population density communities.</div></div>","PeriodicalId":54763,"journal":{"name":"Journal of Theoretical Biology","volume":"615 ","pages":"Article 112248"},"PeriodicalIF":2.0,"publicationDate":"2025-08-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144979427","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}