Richárd Kicsiny, Tamás Bódai, László Székely, Zoltán Varga
{"title":"Extended Discrete-Time Population Model to Describe the Competition of Nutrient-Producing Protocells.","authors":"Richárd Kicsiny, Tamás Bódai, László Székely, Zoltán Varga","doi":"10.1007/s11538-025-01488-0","DOIUrl":"10.1007/s11538-025-01488-0","url":null,"abstract":"<p><p>Modeling the behavior of simple communities of protocells (as basic life-like organisms) is of vital importance since their better understanding may help to describe more complex (artificial and real) ecological systems. In this paper, we extend a recently developed discrete-time dynamic population model (called preliminary model) to a more general, completely reformulated version for describing the competition in a community of three protocell species (one generalist and two specialists). The advantage for the generalist is that it produces more kinds of nutrients than the specialists. In contrast to the preliminary model, the reproduction times and the times of (first) appearance of the three species can be all different in the extended model. The aim is to achieve the most basic/fundamental model that already displays complex population phenomena, like competitive exclusion, keystone species and an interesting \"anomaly\" regarding the connection between the survival of certain species and the decreasing rates of certain nutrients in the environment. Although, we could achieve this aim with a three-species model, at the simplest level, the model can be easily extended for more species in the future. The mentioned \"anomaly\" is a new discovery as it was not observed in the preliminary model. A particular equilibrium, when only the generalist survives, is exactly analyzed, where, interestingly, the golden ratio arises regarding the densities of the protocells of different ages. In future works, the extended model may serve as a useful tool for studying further phenomena in ecosystems, in their pure/abstract form.</p>","PeriodicalId":9372,"journal":{"name":"Bulletin of Mathematical Biology","volume":"87 8","pages":"111"},"PeriodicalIF":2.2,"publicationDate":"2025-07-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC12274253/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144658471","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}
Daniel J Netherwood, Alexander K Y Tam, Campbell W Gourlay, Tea Knežević, Jennifer M Gardner, Vladimir Jiranek, Benjamin J Binder, J Edward F Green
{"title":"Accidental and Regulated Cell Death in Yeast Colony Biofilms.","authors":"Daniel J Netherwood, Alexander K Y Tam, Campbell W Gourlay, Tea Knežević, Jennifer M Gardner, Vladimir Jiranek, Benjamin J Binder, J Edward F Green","doi":"10.1007/s11538-025-01470-w","DOIUrl":"10.1007/s11538-025-01470-w","url":null,"abstract":"<p><p>The yeast species Saccharomyces cerevisiae is one of the most intensively studied organisms on the planet due to it being an excellent eukaryotic model organism in molecular and cell biology. In this work, we investigate the growth and morphology of yeast colony biofilms, where proliferating yeast cells reside within a self-produced extracellular matrix. This research area has garnered significant scientific interest due to its applicability in the biological and biomedical sectors. A central feature of yeast colony biofilm expansion is cellular demise, which is onset by one of two independent mechanisms: either accidental cell death (ACD) or regulated cell death (RCD). In this article, we generalise a continuum model for the nutrient-limited growth of a yeast colony biofilm to include the effects of ACD and RCD. This new model involves a system of four coupled nonlinear reaction-diffusion equations for the yeast-cell density, the nutrient concentration, and two species of dead cells. Numerical solutions of the spatially one and two-dimensional governing equations reveal the impact that ACD and RCD have on expansion speed, morphology and cell distribution within the colony biofilm. Our results are in good qualitative agreement with our own experiments.</p>","PeriodicalId":9372,"journal":{"name":"Bulletin of Mathematical Biology","volume":"87 8","pages":"110"},"PeriodicalIF":2.2,"publicationDate":"2025-07-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC12271256/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144648657","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}
Ariane Ernst, Anastasia Bankowski, Meida Jusyte, Toluwani Okunola, Tino Petrov, Alexander M Walter, Stefanie Winkelmann
{"title":"Parameter Optimization for a Neurotransmission Recovery Model.","authors":"Ariane Ernst, Anastasia Bankowski, Meida Jusyte, Toluwani Okunola, Tino Petrov, Alexander M Walter, Stefanie Winkelmann","doi":"10.1007/s11538-025-01486-2","DOIUrl":"10.1007/s11538-025-01486-2","url":null,"abstract":"<p><p>We assess the empirical applicability of a simplified model for neurotransmitter release that incorporates maturation, fusion, and recovery of both release sites and vesicles. Model parameters are optimized by fitting the model to experimental data obtained from neuromuscular junction synapses of 3rd-instar Drosophila melanogaster larvae. In particular, the mean-squared error between the local extrema of the simulated total junction current and its experimental counterpart is minimized. We compare three estimation approaches, differing in the choice of optimized parameters and the fusion rate function. Despite the model's minimalistic structure, it demonstrates a compelling ability to replicate experimental data, yielding plausible parameter estimates for five different animals. An additional identifiability analysis based on the profile likelihood reveals practical non-identifiabilities for several parameters, highlighting the need for additional constraints or data to improve estimation accuracy.</p>","PeriodicalId":9372,"journal":{"name":"Bulletin of Mathematical Biology","volume":"87 8","pages":"109"},"PeriodicalIF":2.2,"publicationDate":"2025-07-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144574872","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}
{"title":"The Detailed Balance Property and Chemical Systems out of Equilibrium.","authors":"E Franco, J J L Velázquez","doi":"10.1007/s11538-025-01487-1","DOIUrl":"10.1007/s11538-025-01487-1","url":null,"abstract":"<p><p>The detailed balance property is a fundamental property that must be satisfied in all the macroscopic systems with a well defined temperature at each point. On the other hand, many biochemical networks work in non-equilibrium conditions and they can be effectively modelled using sets of equations in which the detailed balance condition fails. In this paper we study a class of \"out of equilibrium\" chemical networks that can be obtained freezing the concentration of some substances in chemical networks for which the detailed balance property holds. In particular, we prove that any chemical system with bidirectional chemical reactions can be extended to a system having additional substances and for which the detailed balance property holds.</p>","PeriodicalId":9372,"journal":{"name":"Bulletin of Mathematical Biology","volume":"87 8","pages":"108"},"PeriodicalIF":2.2,"publicationDate":"2025-07-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC12213986/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144539055","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}
Molly Brennan, Andrew L Krause, Edgardo Villar-Sepúlveda, Christopher B Prior
{"title":"Pattern Formation as a Resilience Mechanism in Cancer Immunotherapy.","authors":"Molly Brennan, Andrew L Krause, Edgardo Villar-Sepúlveda, Christopher B Prior","doi":"10.1007/s11538-025-01485-3","DOIUrl":"10.1007/s11538-025-01485-3","url":null,"abstract":"<p><p>Mathematical and computational modelling in oncology has played an increasingly important role in not only understanding the impact of various approaches to treatment on tumour growth, but in optimizing dosing regimens and aiding the development of treatment strategies. However, as with all modelling, only an approximation is made in the description of the biological and physical system. Here we show that tissue-scale spatial structure can have a profound impact on the resilience of tumours to immunotherapy using a classical model incorporating IL-2 compounds and effector cells as treatment parameters. Using linear stability analysis, numerical continuation, and direct simulations, we show that diffusing cancer cell populations can undergo pattern-forming (Turing) instabilities, leading to spatially-structured states that persist far into treatment regimes where the corresponding spatially homogeneous systems would uniformly predict a cancer-free state. These spatially-patterned states persist in a wide range of parameters, as well as under time-dependent treatment regimes. Incorporating treatment via domain boundaries can increase this resistance to treatment in the interior of the domain, further highlighting the importance of spatial modelling when designing treatment protocols informed by mathematical models. Counter-intuitively, this mechanism shows that increased effector cell mobility can increase the resilience of tumours to treatment. We conclude by discussing practical and theoretical considerations for understanding this kind of spatial resilience in other models of cancer treatment, in particular those incorporating more realistic spatial transport. This paper belongs to the special collection: Problems, Progress and Perspectives in Mathematical and Computational Biology.</p>","PeriodicalId":9372,"journal":{"name":"Bulletin of Mathematical Biology","volume":"87 8","pages":"106"},"PeriodicalIF":2.2,"publicationDate":"2025-07-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC12214011/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144539054","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}
{"title":"Investigating Tumour Responses to Combinations of Radiotherapy and Hyperthermia.","authors":"Chloé Colson, Philip K Maini, Helen M Byrne","doi":"10.1007/s11538-025-01449-7","DOIUrl":"10.1007/s11538-025-01449-7","url":null,"abstract":"<p><p>Hyperthermia (HT) is a promising candidate for enhancing the efficacy of radiotherapy (RT), but its use in the clinic has been limited by incomplete understanding of its interactions with RT. In this work, we investigate tumour responses to high temperature HT alone and combined with RT, focussing on how two different mechanisms for growth control may impact tumour sensitivity to these treatments. We extend an existing ordinary differential equation model of tumour growth and RT response to include high HT. In the absence of treatment, this model distinguishes between growth arrest due to nutrient insufficiency and competition for space, and exhibits three growth regimes: nutrient limited (NL), space limited (SL) and bistable (BS), where both mechanisms for growth arrest coexist. We construct three virtual tumour populations corresponding to the NL, SL and BS regimes and, for each population, we identify the treatment (RT, HT or RT + HT) and dosing regimen that maximise the reduction in tumour burden at the treatment end-point. We thus generate experimentally testable predictions that may explain highly variable experimental and clinical responses to RT and HT and assist patient-specific treatment design.</p>","PeriodicalId":9372,"journal":{"name":"Bulletin of Mathematical Biology","volume":"87 8","pages":"107"},"PeriodicalIF":2.2,"publicationDate":"2025-07-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC12213896/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144539053","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}
{"title":"Inferring Birth Versus Death Dynamics for Ecological Interactions in Stochastic Heterogeneous Populations.","authors":"Erin Beckman, Heyrim Cho, Linh Huynh","doi":"10.1007/s11538-025-01477-3","DOIUrl":"10.1007/s11538-025-01477-3","url":null,"abstract":"<p><p>In this paper, we study the significance of ecological interactions and separation of birth and death dynamics in stochastic heterogeneous populations via general birth-death processes. Interactions can manifest through the birth dynamics, the death dynamics, or some combination of the two. The underlying microscopic mechanisms are important but often implicit in population-level data. We propose an inference method for disambiguating the types of interaction and the birth and death processes from population size time series data of a stochastic n-type heterogeneous population. The interspecies interactions considered can be competitive, antagonistic, or mutualistic. We show that different pairs of birth and death rates with the same net growth rate result in different time series statistics. Then, the inference method is validated in the example of a birth-death process inspired by the two-type Lotka-Volterra interaction dynamics. Utilizing stochastic fluctuations enables us to estimate additional parameters in this stochastic Lotka-Volterra model, which are not identifiable in a deterministic model.</p>","PeriodicalId":9372,"journal":{"name":"Bulletin of Mathematical Biology","volume":"87 8","pages":"105"},"PeriodicalIF":2.2,"publicationDate":"2025-06-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144526511","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}
{"title":"On Modeling Ovarian Aging and Menopause Timing.","authors":"Sean D Lawley, Nanette Santoro, Joshua Johnson","doi":"10.1007/s11538-025-01481-7","DOIUrl":"10.1007/s11538-025-01481-7","url":null,"abstract":"<p><p>Mathematical modeling of ovarian aging and menopause timing has a long history, dating back a half-century to the models of Nobel Prize winner Robert G. Edwards. More recently, such models have been used to investigate clinical interventions for women, which underscores the importance of scientific rigor in model development and analysis. In this paper, we analyze a recent model published in the biophysics literature. We first correct an error which invalidates claims about menopause age in different populations. We then use stochastic analysis to show how this model is a reparameterization of a prior model and put it in the framework of several prior models, which enables the application of extreme value theory. We prove some general extreme value theory results and use them to obtain detailed estimates of menopause age in this model. In particular, we derive a new expected menopause age formula which is orders of magnitude more accurate than the previous heuristic estimate. We further obtain rigorous analytical estimates of the full menopause age distribution and all its moments. We conclude by using these mathematical results to elucidate the physiological sources of menopause age variability.</p>","PeriodicalId":9372,"journal":{"name":"Bulletin of Mathematical Biology","volume":"87 8","pages":"104"},"PeriodicalIF":2.2,"publicationDate":"2025-06-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144526512","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}
Hiroshi Toki, Yoshiharu Yonekura, Yuichi Tsunoyama, Masako Bando
{"title":"Theoretical Study of Retinoblastoma in the Hereditary and Non-hereditary Processes Including the Cancer Growth.","authors":"Hiroshi Toki, Yoshiharu Yonekura, Yuichi Tsunoyama, Masako Bando","doi":"10.1007/s11538-025-01483-5","DOIUrl":"10.1007/s11538-025-01483-5","url":null,"abstract":"<p><p>The two-hit model proposed by Knudson for retinoblastoma has been widely recognized as a standard model for cancer incidence. It successfully predicted the existence of the tumor suppressor gene known as \"Rb1\" by effectively demonstrating the overall patterns observed in clinical data covering both bilateral and unilateral retinoblastoma cases. However, it is important to note that the model's prediction currently deviates significantly from clinical data, both qualitatively and quantitatively. Regrettably, this disparity has remained unresolved. In light of this, we conducted a thorough re-evaluation of Knudson's two-hit model and arrived at a plausible solution that an additional somatic mutation mechanism is required to accurately replicate the magnitude and age dependence observed in both bilateral and unilateral retinoblastoma cases. This revelation offers a fresh and valuable perspective on the development of cancer, highlighting the significance of mutations not only during the cell growth period but also after the retina organ has reached maturity. We refer to this phase as the \"mature period,\" during which the mutation rate has been observed to surpass that of the growth period. With this enhanced understanding of retinoblastoma (Rb), we believe we have shed light on the intricate relationship between somatic and germline mutations. Moreover, this insight provides a promising clue for further exploration into the broader context of cancer incidence resulting from genetic mutations.</p>","PeriodicalId":9372,"journal":{"name":"Bulletin of Mathematical Biology","volume":"87 8","pages":"103"},"PeriodicalIF":2.2,"publicationDate":"2025-06-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC12202645/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144494687","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}
{"title":"Nutrient-Driven Adaptive Evolution of Foraging Traits Impacts Producer-Grazer Dynamics.","authors":"Oluwagbemisola Oladepo, Angela Peace","doi":"10.1007/s11538-025-01482-6","DOIUrl":"10.1007/s11538-025-01482-6","url":null,"abstract":"<p><p>This study investigates the nutrient-driven adaptability of foraging efforts in producer-grazer dynamics. We develop two stoichiometric producer-grazer models: a base model incorporating a fixed energetic cost of feeding and an adaptive model where feeding costs vary over time in response to environmental conditions. By comparing these models, we examine the effects of adaptive foraging strategies on population dynamics. Our adaptive model suggests a potential mechanism for evolutionary rescue, where the population dynamically adjusts to environmental changes, such as fluctuations in food quality, by modifying its feeding strategies. However, when population densities oscillate in predator-prey limit cycles, fast adaptation can lead to very wide amplitude cycles, where populations are in danger of stochastic extinction. Overall, this increases our understanding of the conditions under which nutrient-driven adaptive foraging strategies can yield benefits to grazers.</p>","PeriodicalId":9372,"journal":{"name":"Bulletin of Mathematical Biology","volume":"87 8","pages":"102"},"PeriodicalIF":2.2,"publicationDate":"2025-06-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC12198312/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144483217","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}