Mathematical Biosciences and Engineering最新文献

{"title":"Evolutionary dynamics of acid resistance in tumors: A mathematical model.","authors":"Prithvi Anickode, Fabio Augusto Milner","doi":"10.3934/mbe.2026023","DOIUrl":"10.3934/mbe.2026023","url":null,"abstract":"<p><p>Acidosis in tumors arises from reprogrammed metabolism and compromised vasculature, creating a harsh, acidic microenvironment that drives the evolutionary selection of acid-resistant cell phenotypes. A mathematical model was proposed to integrate phenotypic evolution, microenvironmental acidification, and tumor density dynamics. Three key mechanisms were incorporated in it: frequency-dependent selection favoring acid-resistant cells below a critical pH, stress-induced phenotypic switching, and a positive feedback loop where resistant cells produce excess acid that intensifies selection pressure. The well-posedness of the model was established. Through numerical simulations across biologically relevant parameter regimes, we identified two therapeutically targetable parameters: the baseline acid clearance rate (a proxy for vascular perfusion) and a protection factor (representing acid-resistance \"machinery\" effectiveness) as critical bifurcation parameters for resistance evolution. The model exhibits qualitatively distinct dynamics depending on phenotypic plasticity levels. In low-plasticity tumors, both parameters exhibit sharp bifurcations with strong parameter interactions: clearance and protection effects are context-dependent, with therapeutic interventions effective only within specific parameter ranges. In high-plasticity tumors, both parameters produce continuous, monotonic responses with independent, additive effects. These regime-dependent dynamics suggest that treatment strategies should adapt to tumor plasticity: in the former, targeting perfusion alone is typically sufficient, though sequential therapy may be required if the perfusion rate approaches or exceeds the bifurcation threshold, whereas in the latter, treatment might benefit from combination therapies addressing both parameters simultaneously. These findings suggest that a low-dimensional model can identify therapeutically targetable parameters governing resistance evolution, suggesting interventions that may prevent or reverse the harmful effect of acid-resistant phenotypes, which are associated with chemotherapy failure, immune evasion, and metastatic progression.</p>","PeriodicalId":49870,"journal":{"name":"Mathematical Biosciences and Engineering","volume":"23 3","pages":"594-618"},"PeriodicalIF":2.6,"publicationDate":"2026-01-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"147366980","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":"Effects of multi-component small permanent charges on the dynamics of Poisson-Nernst-Planck models.","authors":"Guojian Lin","doi":"10.3934/mbe.2026025","DOIUrl":"https://doi.org/10.3934/mbe.2026025","url":null,"abstract":"<p><p>In this paper, the classical Poisson-Nernst-Planck (PNP) model describing ion transport through a membrane channel is used to study the effects of small permanent charges and the structures of ion channels on ionic flows. The model under study includes two oppositely charged ion species, and the permanent charge in this model is a piecewise constant function with two nonzero regions. By rescaling, the classical PNP model can be viewed as a singularly perturbed differential equation system. Therefore, the geometric singular perturbation theory is employed to get a singular orbit. Assuming that the permanent charge density is small, a regular perturbation expansion is used to obtain the first-order approximation of the individual flux, which acts as a basis for our analysis. Then, the effects of small permanent charges on the fluxes and the current-voltage relation, which not only depend on the boundary conditions, but also depend on the structures of ion channels and the ratio between two nonzero permanent charge densities, are analyzed in this paper. Particularly, our results indicate that the geometric structures of three-dimensional ion channels have a short and narrow cross-section, which is explained in [1]. Also, our results indicate that the ratio between two nonzero permanent charge densities can change the position of a short and narrow cross-section in ion channels.</p>","PeriodicalId":49870,"journal":{"name":"Mathematical Biosciences and Engineering","volume":"23 3","pages":"636-677"},"PeriodicalIF":2.6,"publicationDate":"2026-01-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"147366971","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":"Bayesian chain graph models to characterize microbe-environment dynamics.","authors":"Yunyi Shen, Claudia Solís-Lemus","doi":"10.3934/mbe.2026020","DOIUrl":"10.3934/mbe.2026020","url":null,"abstract":"<p><p>Microbiome data require statistical models that can simultaneously decode microbes' reaction to the environment and interactions among microbes. While a multiresponse linear regression model seems like a straight-forward solution, we argue that treating it as a graphical model is problematic given that the regression coefficient matrix does not encode the conditional dependence structure between response and predictor nodes. This observation is especially important in biological settings when we have prior knowledge on the edges from specific experimental interventions that can only be properly encoded under a conditional dependence model. Here, we propose a chain graph model with two sets of nodes (predictors and responses) whose solution yields a graph with edges that indeed represent conditional dependence, thus agreeing with the experimenter's intuition on the average behavior of nodes under treatment. The solution to our model is sparse via the Bayesian linear regression (LASSO). In addition, we propose an adaptive extension so that different shrinkages can be applied to different edges to incorporate edge-specific prior knowledge. Our model is computationally inexpensive through an efficient Gibbs sampling algorithm and can account for binary, counting, and compositional responses via an appropriate hierarchical structure. We test the performance of our model in a variety of simulated datasets, thereby showing superior performance to state-of-the-art approaches. We further apply our model to human gut and soil microbial compositional datasets, and we highlight that CG-LASSO can estimate biologically meaningful network structures in the data. Our software is available as an R package at https://github.com/YunyiShen/CAR-LASSO.</p>","PeriodicalId":49870,"journal":{"name":"Mathematical Biosciences and Engineering","volume":"23 2","pages":"499-546"},"PeriodicalIF":2.6,"publicationDate":"2026-01-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146087876","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":"Augmented data and neural networks for robust epidemic forecasting: Application to COVID-19 in Italy.","authors":"Giacomo Dimarco, Federica Ferrarese, Lorenzo Pareschi","doi":"10.3934/mbe.2026019","DOIUrl":"10.3934/mbe.2026019","url":null,"abstract":"<p><p>In this work, we propose a data augmentation strategy aimed at improving the training phase of neural networks and, consequently, the accuracy of their predictions. Our approach relies on generating synthetic data through a suitable compartmental model combined with the incorporation of uncertainty. Available data are used to calibrate the model, which is further integrated with deep learning techniques to produce additional synthetic data for training. The results show that neural networks trained on these augmented datasets exhibit significantly improved predictive performances. In particular, we focus on two different neural network architectures: Physics-Informed Neural Networks (PINNs) and Nonlinear Autoregressive (NAR) models. The NAR approach proves especially effective for short-term forecasting, thereby providing accurate quantitative estimates by directly learning the dynamics from data and avoiding the additional computational cost of embedding physical constraints into the training. In contrast, PINNs yield less accurate quantitative predictions but capture the qualitative long-term behavior of the system, thus making them more suitable to explore broader dynamical trends. Numerical simulations of the second phase of the COVID-19 pandemic in the Lombardy region (Italy) validate the effectiveness of the proposed approach.</p>","PeriodicalId":49870,"journal":{"name":"Mathematical Biosciences and Engineering","volume":"23 2","pages":"474-498"},"PeriodicalIF":2.6,"publicationDate":"2026-01-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146087802","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":"Eye2Heart: A reduced mathematical model bridging cardiovascular and ocular hemodynamics.","authors":"Lorenzo Sala, Mohamed Zaid, Faith Hughes, Marcela Szopos, Virginia H Huxley, Alon Harris, Giovanna Guidoboni, Sergey Lapin","doi":"10.3934/mbe.2026017","DOIUrl":"10.3934/mbe.2026017","url":null,"abstract":"<p><p>The cardiovascular and ocular systems are intricately connected, with hemodynamic interactions playing a crucial role in both physiological regulation and pathological conditions. However, existing models often treat these systems separately, thus limiting the understanding of their interdependence. In this study, we present the Eye2Heart model, which is a novel closed-loop mathematical framework that integrates cardiovascular and ocular dynamics. Using an electrical-hydraulic analogy, the model describes the interactions between the heart and retinal circulation through a nonlinear system of ordinary differential equations. The model is tested against clinical and experimental data, thus demonstrating its ability to reproduce key cardiovascular parameters (e.g., stroke volume, cardiac output) and ocular hemodynamics (e.g., retinal blood flow). Additionally, we explore in silico the effects of intraocular pressure and left ventricular compliance on both local ocular and global systemic circulation, thus revealing critical dependencies between cardiovascular and ocular health. The results highlight the model's potential for studying cardiovascular diseases with ocular manifestations and support emerging research in oculomics by providing a mechanistic basis to interpret ocular biomarkers within a systemic context. This paves the way for patient-specific data integration and broader applications in personalized medicine.</p>","PeriodicalId":49870,"journal":{"name":"Mathematical Biosciences and Engineering","volume":"23 2","pages":"421-448"},"PeriodicalIF":2.6,"publicationDate":"2026-01-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146087846","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":"Qualitative optimization of oncolytic virotherapy and immune therapy combination treatments.","authors":"Negar Mohammadnejad, Thomas Hillen","doi":"10.3934/mbe.2026016","DOIUrl":"10.3934/mbe.2026016","url":null,"abstract":"<p><p>Oncolytic viruses (OVs) are designed to selectively target and destroy cancer cells while sparing normal, healthy tissue. Several viruses for oncolytic virotherapy are currently developed. In this paper, we will use mathematical modeling to consider key strategies that can improve the efficacy of oncolytic virotherapy. These include the integration of immunotherapy approaches with virotherapy to amplify anti-tumor immune responses, as well as optimizing the timing, dosage, and sequencing of viral administrations. Specifically, we consider strategies that increase the burst size of the virus, immunostimulation and immunosuppression, we optimize for different weekly virus injection schedules, and we consider the combination of OV therapy with chimeric antigen receptor (CAR) T-cell therapy. A limiting factor is the availability of data. We parametrize the model using several different data sets. These, however, correspond to different cancers and experimental setups. Hence our model cannot be considered to be validated. Consequently, our results are qualitative. Our results highlight the critical importance of timing for virotherapy's efficacy and overall success. They outline strong evidence for promising treatment scenarios that needs to be further tested experimentally in the future.</p>","PeriodicalId":49870,"journal":{"name":"Mathematical Biosciences and Engineering","volume":"23 2","pages":"388-420"},"PeriodicalIF":2.6,"publicationDate":"2026-01-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146087835","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":"Stability and bifurcation of difference equations from stochastic logistic models.","authors":"Haiyan Wang, Emily Wang","doi":"10.3934/mbe.2026018","DOIUrl":"https://doi.org/10.3934/mbe.2026018","url":null,"abstract":"<p><p>This paper investigates the stochastic logistic difference equation, $ X_{n+1} = r X_n (1 - X_n)varepsilon_n $, where $ X_n $ is a random variable of population size, and $ {varepsilon_n} $ represents independent random perturbations with $ E[varepsilon_n] = 1 $ and $ E[varepsilon_n^2] = v > 1 $. Under the Gaussian moment-closure approximation, we derived a closed system of difference equations for the mean and variance of $ X_n $. The analysis of the system of difference equations identified two classes of equilibria: a trivial equilibrium $ (0, 0) $ representing extinction, and nontrivial equilibria corresponding to positive steady population levels. Explicit conditions for the existence and local stability of these equilibria were obtained, showing that the extinction state is stable when $ r^2v < 1 $, whereas nontrivial equilibria arise for $ r > 1 $ with stability dependent on the stochastic intensity $ v $. The saddle-node (fold) bifurcation induced by variations in the stochastic intensity $ v $ was explicitly formulated. Monte Carlo simulations confirmed the analytical analysis.</p>","PeriodicalId":49870,"journal":{"name":"Mathematical Biosciences and Engineering","volume":"23 2","pages":"449-473"},"PeriodicalIF":2.6,"publicationDate":"2026-01-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146087805","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":"Incentives for self-isolation based on incidence rather than prevalence could help to flatten the curve: A modeling study.","authors":"Giulia de Meijere, Hugo Martin","doi":"10.3934/mbe.2026014","DOIUrl":"10.3934/mbe.2026014","url":null,"abstract":"<p><p>In recent years, numerous advances have been made in understanding how epidemic dynamics are affected by changes in individual behaviors. We propose a Susceptible-Infected-Susceptible (SIS) based compartmental model to tackle the simultaneous and coupled evolution of an outbreak and of the adoption by individuals of the isolation measure. The compliance with self-isolation is described with the help of the imitation dynamics framework. Individuals are incentivized to isolate based on the prevalence and the incidence rate of the outbreak, and are tempted to defy isolation recommendations depending on the duration of the isolation and on the cost of putting social interactions on hold. We are able to derive analytical results on the equilibria of the model under the homogeneous mean-field approximation. Simulating the compartmental model on empirical networks, we also perform a preliminary check of the impact of a network structure on our analytical predictions. We find that the dynamics collapse to surprisingly simple regimes where either the imitation dynamics no longer plays a role or the equilibrium prevalence depends on only two parameters of the model, namely the cost and the relative time spent in isolation. Whether individuals prioritize disease prevalence or incidence as an indicator of the state of the outbreak appears to play no role on the equilibria of the dynamics. However, it turns out that favoring incidence may help to flatten the curve in the transient phase of the dynamics. We also find a fair agreement between our analytical predictions and simulations run on an empirical multiplex network.</p>","PeriodicalId":49870,"journal":{"name":"Mathematical Biosciences and Engineering","volume":"23 2","pages":"333-365"},"PeriodicalIF":2.6,"publicationDate":"2026-01-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146087818","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":"An enhanced triclustering $ delta $-Trimax method with fuzzy cuckoo search based on Lévy flight and Gaussian distribution for gene expression data.","authors":"Titin Siswantining, Muhamad Ido Raskapati, Nisa Nurul Hidayah, Gianinna Ardaneswari, Saskya Mary Soemartojo, Siti Nurrohmah, Devvi Sarwinda, Setia Pramana","doi":"10.3934/mbe.2026015","DOIUrl":"10.3934/mbe.2026015","url":null,"abstract":"<p><p>The triclustering method employed in this study integrates the $ delta $-Trimax approach with the fuzzy cuckoo search (FCS), thereby leveraging the Lévy flight and Gaussian distribution to analyze gene expression data in three dimensions. In this framework, the initial triclusters produced by $ delta $-Trimax are further optimized using FCS, where the Lévy flight enhances global exploration and the Gaussian distribution intensifies local exploitation, thus achieving a balanced search for optimal solutions. Each tricluster set is evaluated using the tricluster quality index (TQI) to ensure coherence across genes, conditions, and time points. The method was applied to gene expression datasets from primary fibroblast cells and heart disease samples. In the fibroblast dataset, the best tricluster set was obtained with $ delta = 0.015 $ and yielded the lowest average TQI value. For the heart disease dataset, the most optimal solution was achieved with $ delta = 0.026 $, which yielded the lowest average TQI, and the best tricluster showed large gene coverage across multiple time points. A functional analysis of the selected triclusters using gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways uncovered significant enrichment in pathways such as the NF-$ kappa $B signaling pathway (hsa04064), TGF-$ beta $ signaling pathway (hsa04350), and calcium signaling pathway (hsa04020), all of which are mechanistically relevant to immune modulation, extracellular matrix organization, and cardiac muscle function. These findings highlight the utility of the proposed hybrid framework in uncovering biologically meaningful gene modules and provide valuable insights into the molecular mechanisms underlying fibrotic and cardiovascular diseases.</p>","PeriodicalId":49870,"journal":{"name":"Mathematical Biosciences and Engineering","volume":"23 2","pages":"366-387"},"PeriodicalIF":2.6,"publicationDate":"2026-01-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146087878","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 the statistical description of large populations.","authors":"Yuri Kozitsky, Krzysztof Pilorz","doi":"10.3934/mbe.2026013","DOIUrl":"https://doi.org/10.3934/mbe.2026013","url":null,"abstract":"<p><p>A wide variety of works exists on the dynamics of large populations, ranging from simple heuristic modeling to those based on advanced computer-supported methods. However, their interconnections remain mostly vague, which, in particular, limits the effectiveness of computer methods in this domain. In this work, we propose and justify the following concept. Typically, the description of the population dynamics is based on the sole use of low-order correlations. As we demonstrate here, in important cases, where the local population structure is shaped by strong interactions, higher-order correlations become essential. To verify when one can or cannot rely on studying low-order correlations only, we suggest to explicitly use probability measures as micro-states. Among such states may be those whose adequate characterization is based on their low-order correlation functions. In particular, this is the case for sub-Poissonian states where the large $ n $ asymptotics of the probability of finding $ n $ particles in a given vessel is similar to that for non-interacting entities, which can completely be described by the density of the particles. To illustrate this concept, a general individual-based model of an infinite population of interacting entities is analyzed. The evolution of this model preserves the sub-Poissonian states, which allows one to describe it through the correlation functions of such states for which a chain of evolution equations is obtained. The corresponding kinetic equation is derived, numerically solved, and analyzed.</p>","PeriodicalId":49870,"journal":{"name":"Mathematical Biosciences and Engineering","volume":"23 2","pages":"312-332"},"PeriodicalIF":2.6,"publicationDate":"2026-01-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146087851","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}