{"title":"Extending Mathematical Frameworks to Investigate Neuronal Dynamics in the Presence of Microglial Ensheathment.","authors":"Nellie Garcia, Silvie Reitz, Gregory Handy","doi":"10.1007/s11538-025-01438-w","DOIUrl":"https://doi.org/10.1007/s11538-025-01438-w","url":null,"abstract":"<p><p>Recent experimental evidence has shown that glial cells, including microglia and astrocytes, can ensheathe specific synapses, positioning them to disrupt neurotransmitter flow between pre- and post-synaptic terminals. This study, as part of the special issue \"Problems, Progress and Perspectives in Mathematical and Computational Biology,\" expands micro- and network-scale theoretical frameworks to incorporate these new experimental observations that introduce substantial heterogeneities into the system. Specifically, we aim to explore how varying degrees of synaptic ensheathment affect synaptic communication and network dynamics. Consistent with previous studies, our microscale model shows that ensheathment accelerates synaptic transmission while reducing its strength and reliability, with the potential to effectively switch off synaptic connections. Building on these findings, we integrate an \"effective\" glial cell model into a large-scale neuronal network. Specifically, we analyze a network with highly heterogeneous synaptic strengths and time constants, where glial proximity parametrizes synaptic properties. This parametrization results in a multimodal distribution of synaptic parameters across the network, introducing significantly greater variability compared to previous modeling efforts that assumed a normal distribution. This framework is applied to large networks of exponential integrate-and-fire neurons, extending linear response theory to analyze not only firing rate distributions but also noise correlations across the network. Despite the significant heterogeneity in the system, a mean-field approximation accurately captures network statistics. We demonstrate the utility of our model by reproducing experimental findings, showing that microglial ensheathment leads to post-anesthesia hyperactivity in excitatory neurons of mice. Furthermore, we explore how glial ensheathment may be used in the visual cortex to target specific neuronal subclasses, tuning higher-order network statistics.</p>","PeriodicalId":9372,"journal":{"name":"Bulletin of Mathematical Biology","volume":"87 5","pages":"63"},"PeriodicalIF":2.0,"publicationDate":"2025-04-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143779216","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":"Impact of ADE and Dengue Vaccination with Screening on Cost and Disease Burden for Homoserotypic Dengue and Zika.","authors":"Christopher M Kribs, Parker Mays","doi":"10.1007/s11538-025-01440-2","DOIUrl":"https://doi.org/10.1007/s11538-025-01440-2","url":null,"abstract":"<p><p>The tetravalent dengue vaccine Dengvaxia<sup>®</sup> may prime dengue-seronegative vaccinees for antibody-dependent enhancement (ADE) of any subsequent dengue (in case of vaccine failure) or Zika infections. Many researchers associate ADE of such cases with more severe outcomes including death. This study uses a mathematical model of transmission dynamics that distinguishes ADE and non-ADE cases for each virus, to identify the potential impact of a dengue screening and vaccination campaign on the economic cost and disease burden of a dual dengue-Zika outbreak, under the hypothesis that severe outcomes are associated with ADE. Results indicate that when all dengue exposure is to a single serotype, in most cases vaccination increases both cost and burden because they are dominated by the high costs associated with complications from ADE Zika cases. However, if per-case ADE Zika costs are lower than estimated (a real possibility given the limited data available), by a factor ranging from 1 to 6 (for cost, except in Vietnam) or 8 (for burden), sufficiently high vaccination coverage can reduce total cost and burden substantially over a year. Analysis also identifies variations across countries, dengue serotypes, and timeframes of evaluation.</p>","PeriodicalId":9372,"journal":{"name":"Bulletin of Mathematical Biology","volume":"87 5","pages":"62"},"PeriodicalIF":2.0,"publicationDate":"2025-04-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143779217","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 Impact of T-cell Exhaustion Dynamics on Tumour-Immune Interactions and Tumour Growth.","authors":"Nicholas Lai, Alexis Farman, Helen M Byrne","doi":"10.1007/s11538-025-01433-1","DOIUrl":"10.1007/s11538-025-01433-1","url":null,"abstract":"<p><p>Tumours evade immune surveillance through a number of different immunosuppressive mechanisms. One such mechanism causes cytotoxic T-cells, a major driving force of the immune system, to differentiate to a state of 'exhaustion', rendering them less effective at killing tumour cells. We present a structured mathematical model that focuses on T-cell exhaustion and its effect on tumour growth. We compartmentalise cytotoxic T-cells into discrete subgroups based on their exhaustion level, which affects their ability to kill tumour cells. We show that the model reduces to a simpler system of ordinary differential equations (ODEs) that describes the time evolution of the total number of T-cells, their mean exhaustion level and the total number of tumour cells. Numerical simulations of the model equations reveal how the exhaustion distribution of T-cells changes over time and how it influences the tumour's growth dynamics. Complementary bifurcation analysis shows how altering key parameters significantly reduces the tumour burden, highlighting exhaustion as a promising target for immunotherapy. Finally, we derive a continuum approximation of the discrete ODE model, which admits analytical solutions that provide complementary insight into T-cell exhaustion dynamics and their effect on tumour growth.</p>","PeriodicalId":9372,"journal":{"name":"Bulletin of Mathematical Biology","volume":"87 5","pages":"61"},"PeriodicalIF":2.0,"publicationDate":"2025-04-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC11965189/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143762953","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":"A Framework for Parameter Estimation and Uncertainty Quantification in Systems Biology Using Quantile Regression and Physics-Informed Neural Networks.","authors":"Haoran Hu, Qianru Cheng, Shuli Guo, Huifang Wen, Jing Zhang, Yongqi Song, Kaiqun Wang, Di Huang, Hui Zhang, Chaofeng Zhang, Yanhu Shan","doi":"10.1007/s11538-025-01439-9","DOIUrl":"https://doi.org/10.1007/s11538-025-01439-9","url":null,"abstract":"<p><p>A framework for parameter estimation and uncertainty quantification is crucial for understanding the mechanisms of biological interactions within complex systems and exploring their dynamic behaviors beyond what can be experimentally observed. Despite recent advances, challenges remain in achieving the high accuracy of parameter estimation and uncertainty quantification at moderate computational costs. To tackle these challenges, we developed a novel approach that integrates the quantile method with Physics-Informed Neural Networks (PINNs). This method utilizes a network architecture with multiple parallel outputs, each corresponding to a distinct quantile, facilitating a comprehensive characterization of parameter estimation and its associated uncertainty. The effectiveness of the proposed approach was validated across three study cases, where it was compared to the Monte Carlo dropout (MCD) and the Bayesian methods. Furthermore, a larger-scale model was employed to further demonstrate the excellent performance of the proposed approach. Our approach exhibited significantly superior efficacy in parameter estimation and uncertainty quantification. This highlights its great promise to broaden the scope of applications in system biology modeling.</p>","PeriodicalId":9372,"journal":{"name":"Bulletin of Mathematical Biology","volume":"87 5","pages":"60"},"PeriodicalIF":2.0,"publicationDate":"2025-03-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143735553","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":"How Missed Doses of Antibiotics Affect Bacteria Growth Dynamics.","authors":"Hwai-Ray Tung, Sean D Lawley","doi":"10.1007/s11538-025-01430-4","DOIUrl":"https://doi.org/10.1007/s11538-025-01430-4","url":null,"abstract":"<p><p>What should you do if you miss a dose of antibiotics? Despite the prevalence of missed antibiotic doses, there is vague or little guidance on what to do when a dose is forgotten. In this paper, we consider the effects of different patient responses after missing a dose using a mathematical model that links antibiotic concentration with bacteria dynamics. We show using simulations that, in some circumstances, (a) missing just a few doses can cause treatment failure, and (b) this failure can be remedied by simply taking a double dose after a missed dose. We then develop an approximate model that is analytically tractable and use it to understand when it might be advisable to take a double dose after a missed dose.</p>","PeriodicalId":9372,"journal":{"name":"Bulletin of Mathematical Biology","volume":"87 5","pages":"58"},"PeriodicalIF":2.0,"publicationDate":"2025-03-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143735555","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}
Daria Stepanova, Meritxell Brunet Guasch, Helen M Byrne, Tomás Alarcón
{"title":"Understanding How Chromatin Folding and Enzyme Competition Affect Rugged Epigenetic Landscapes.","authors":"Daria Stepanova, Meritxell Brunet Guasch, Helen M Byrne, Tomás Alarcón","doi":"10.1007/s11538-025-01434-0","DOIUrl":"https://doi.org/10.1007/s11538-025-01434-0","url":null,"abstract":"<p><p>Epigenetics plays a key role in cellular differentiation and maintaining cell identity, enabling cells to regulate their genetic activity without altering the DNA sequence. Epigenetic regulation occurs within the context of hierarchically folded chromatin, yet the interplay between the dynamics of epigenetic modifications and chromatin architecture remains poorly understood. In addition, it remains unclear what mechanisms drive the formation of rugged epigenetic patterns, characterised by alternating genomic regions enriched in activating and repressive marks. In this study, we focus on post-translational modifications of histone H3 tails, particularly H3K27me3, H3K4me3, and H3K27ac. We introduce a mesoscopic stochastic model that incorporates chromatin architecture and competition of histone-modifying enzymes into the dynamics of epigenetic modifications in small genomic loci comprising several nucleosomes. Our approach enables us to investigate the mechanisms by which epigenetic patterns form on larger scales of chromatin organisation, such as loops and domains. Through bifurcation analysis and stochastic simulations, we demonstrate that the model can reproduce uniform chromatin states (open, closed, and bivalent) and generate previously unexplored rugged profiles. Our results suggest that enzyme competition and chromatin conformations with high-frequency interactions between distant genomic loci can drive the emergence of rugged epigenetic landscapes. Additionally, we hypothesise that bivalent chromatin can act as an intermediate state, facilitating transitions between uniform and rugged landscapes. This work offers a powerful mathematical framework for understanding the dynamic interactions between chromatin architecture and epigenetic regulation, providing new insights into the formation of complex epigenetic patterns.</p>","PeriodicalId":9372,"journal":{"name":"Bulletin of Mathematical Biology","volume":"87 5","pages":"59"},"PeriodicalIF":2.0,"publicationDate":"2025-03-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143735557","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}
Arianna Ceccarelli, Alexander P Browning, Ruth E Baker
{"title":"Approximate Solutions of a General Stochastic Velocity-Jump Model Subject to Discrete-Time Noisy Observations.","authors":"Arianna Ceccarelli, Alexander P Browning, Ruth E Baker","doi":"10.1007/s11538-025-01437-x","DOIUrl":"10.1007/s11538-025-01437-x","url":null,"abstract":"<p><p>Advances in experimental techniques allow the collection of high-resolution spatio-temporal data that track individual motile entities over time. These tracking data motivate the use of mathematical models to characterise the motion observed. In this paper, we aim to describe the solutions of velocity-jump models for single-agent motion in one spatial dimension, characterised by successive Markovian transitions within a finite network of n states, each with a specified velocity and a fixed rate of switching to every other state. In particular, we focus on obtaining the solutions of the model subject to noisy, discrete-time, observations, with no direct access to the agent state. The lack of direct observation of the hidden state makes the problem of finding the exact distributions generally intractable. Therefore, we derive a series of approximations for the data distributions. We verify the accuracy of these approximations by comparing them to the empirical distributions generated through simulations of four example model structures. These comparisons confirm that the approximations are accurate given sufficiently infrequent state switching relative to the imaging frequency. The approximate distributions computed can be used to obtain fast forwards predictions, to give guidelines on experimental design, and as likelihoods for inference and model selection.</p>","PeriodicalId":9372,"journal":{"name":"Bulletin of Mathematical Biology","volume":"87 5","pages":"57"},"PeriodicalIF":2.0,"publicationDate":"2025-03-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC11937228/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143708630","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}
Rowan L Hassman, Iona M H McCabe, Kaia M Smith, Linda J S Allen
{"title":"Correction: Stochastic Models of Zoonotic Avian Influenza with Multiple Hosts, Environmental Transmission, and Migration in the Natural Reservoir.","authors":"Rowan L Hassman, Iona M H McCabe, Kaia M Smith, Linda J S Allen","doi":"10.1007/s11538-025-01428-y","DOIUrl":"10.1007/s11538-025-01428-y","url":null,"abstract":"","PeriodicalId":9372,"journal":{"name":"Bulletin of Mathematical Biology","volume":"87 5","pages":"56"},"PeriodicalIF":2.0,"publicationDate":"2025-03-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143691318","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":"Vaccination and Collective Action Under Social Norms.","authors":"Bryce Morsky","doi":"10.1007/s11538-025-01436-y","DOIUrl":"https://doi.org/10.1007/s11538-025-01436-y","url":null,"abstract":"<p><p>Social dynamics are an integral part of the spread of disease affecting contact rates as well as the adoption of pharmaceutical and non-pharmaceutical interventions. When vaccines provide waning immunity, efficient and timely uptake of boosters is required to maintain protection and reduce infections. How then do social dynamics affect the timely uptake of vaccines and thereby the course of an epidemic? This paper explores this scenario through a behavioural-epidemiological model. It features a tipping-point dynamic for the uptake of vaccines that combines the risk of infection, perceived morbidity risk of the vaccine, and social payoffs for deviating from the vaccination decision-making of others. The social payoffs are derived from a social norm of conformity, and they create a collective action problem. A key finding driven by this dilemma is that waves of vaccine uptake and infections can occur due to inefficient and delayed uptake of boosters. This results in a nonlinear response of the infection load to the transmission rate: an intermediate transmission rate can result in greater prevalence of disease relative to more or less transmissible diseases. Further, global information about the prevalence of the disease and vaccine uptake can increase the infection load and peak relative to information restricted to individuals' contact networks. Thus, decisions driven by local information can mitigate the collective action problem across the population. Finally, the optimal public policy program to promote boosters is shown to be one that focuses on overcoming the social inertia to vaccinate at the start of an epidemic.</p>","PeriodicalId":9372,"journal":{"name":"Bulletin of Mathematical Biology","volume":"87 5","pages":"55"},"PeriodicalIF":2.0,"publicationDate":"2025-03-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143699768","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":"Mathematical Modelling of the First HIV/ZIKV Co-infection Cases in Colombia and Brazil.","authors":"Jhoana P Romero-Leiton, Idriss Sekkak, Julien Arino, Iain Moyles, Bouchra Nasri","doi":"10.1007/s11538-025-01429-x","DOIUrl":"https://doi.org/10.1007/s11538-025-01429-x","url":null,"abstract":"<p><p>This paper presents a mathematical model to investigate the co-infection with human immunodeficiency virus (HIV) and Zika virus (ZIKV) in Colombia and Brazil, where the first cases were reported in 2015. The model considers the sexual transmission dynamics of both viruses and vector-host interactions. We begin by exploring the qualitative behaviour of each model separately. We then analyze the dynamics of the co-infection model using the thresholds and results defined separately for each model. The model also considers the impact of intervention strategies, such as personal protection, antiretroviral therapy (ART), and sexual protection (condom use). Using available and assumed parameter values for Colombia and Brazil, the model is calibrated to investigate the long-term co-infection dynamics, the influence of specific parameters, and the potential effect of implementing these intervention strategies on co-infection spread. The study's results revealed that the duration of Zika infection is a critical factor influencing the burden of co-infection cases. Additionally, bed nets and use of condoms are essential for disease control, while ART is less emphasized due to the cost-effectiveness of condom use.</p>","PeriodicalId":9372,"journal":{"name":"Bulletin of Mathematical Biology","volume":"87 5","pages":"54"},"PeriodicalIF":2.0,"publicationDate":"2025-03-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143673386","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}