Mathematical Biosciences

We propose new single and two-strain epidemic models represented by systems of delay differential equations and based on the number of newly exposed individuals. Transitions between exposed, infectious, recovered, and back to susceptible compartments are determined by the corresponding time delays. Existence and positiveness of solutions are proved. Reduction of delay differential equations to integral equations allows the analysis of stationary solutions and their stability. In the case of two strains, they compete with each other, and the strain with a larger individual basic reproduction number dominates the other one...

In this paper, we present a virus infection model that incorporates eclipse-stage and Beddington-DeAngelis function, along with perturbation in infection rate using logarithmic Ornstein-Uhlenbeck process. Rigorous analysis demonstrates that the stochastic model has a unique global solution. Through construction of appropriate Lyapunov functions and a compact set, combined with the strong law of numbers and Fatou's lemma, we obtain the existence of the stationary distribution under a critical condition, which indicates the long-term persistence of T-cells and virions...

Plants in arid environments have evolved many strategies to resist drought. Among them, the developed water storage tissue is an essential characteristic of xerophytes. To clarify the role of water storage capacity in plant performance, we originally formulate a stoichiometric model to describe the interaction between plants and water with explicit water storage. Via an ecological reproductive index, we explore the effects of precipitation and water storage capacity on plant dynamics. The model possesses saddle-node bifurcation and forward or backward bifurcation, and the latter may lead to the emergence of alternative stable states between a stable survival state and a stable extinction state...

Oncologists and applied mathematicians are interested in understanding the dynamics of cancer-immune interactions, mainly due to the unpredictable nature of tumour cell proliferation. In this regard, mathematical modelling offers a promising approach to comprehend this potentially harmful aspect of cancer biology. This paper presents a novel dynamical model that incorporates the interactions between tumour cells, healthy tissue cells, and immune-stimulated cells when subjected to simultaneous chemotherapy and radiotherapy for treatment...

We consider a three-dimensional mathematical model that describes the interaction between the effector cells, tumor cells, and the cytokine (IL-2) of a patient. This is called the Kirschner-Panetta model. Our objective is to explain the tumor oscillations in tumor sizes as well as long-term tumor relapse. We then explore the effects of adoptive cellular immunotherapy on the model and describe under what circumstances the tumor can be eliminated or can remain over time but in a controlled manner. Nonlinear dynamics of immunogenic tumors are given, for example: we prove that the trajectories of the associated system are bounded and defined for all positive time; there are some invariant subsets; there are open subsets of parameters, such that the system in the first octant has at most five equilibrium solutions, one of them is tumor-free and the others are of co-existence...

SARS-CoV-2 has the ability to form large multi-nucleated cells known as syncytia. Little is known about how syncytia affect the dynamics of the infection or severity of the disease. In this manuscript, we extend a mathematical model of cell-cell fusion assays to estimate both the syncytia formation rate and the average duration of the fusion phase for five strains of SARS-CoV-2. We find that the original Wuhan strain has the slowest rate of syncytia formation (6.4×10-4 /h), but takes only 4.0 h to complete the fusion process, while the Alpha strain has the fastest rate of syncytia formation (0...

This article addresses reaction networks in which spatial and stochastic effects are of crucial importance. For such systems, particle-based models allow us to describe all microscopic details with high accuracy. However, they suffer from computational inefficiency if particle numbers and density get too large. Alternative coarse-grained-resolution models reduce computational effort tremendously, e.g., by replacing the particle distribution by a continuous concentration field governed by reaction-diffusion PDEs...

This paper addresses the problem of designing a safe and optimal control strategy for typical cancer using the Control Barrier Function (CBF) technique. Cancer is a complex and highly dynamic disease characterized by uncontrolled cell growth and proliferation. By formulating the cancer dynamics as a control system, this study introduces a CBF-based controller that guides the cancerous tissue toward safe and controlled behaviors. The controller is designed to simultaneously optimize treatment efficacy and patient safety...

Based on a deterministic and stochastic process hybrid model, we use white noises to account for patient variabilities in treatment outcomes, use a hyperparameter to represent patient heterogeneity in a cohort, and construct a stochastic model in terms of Ito stochastic differential equations for tesing the efficacy of three different treatment protocols in CAR T cell therapy. The stochastic model has three ergodic invariant measures which correspond to three unstable equilibrium solutions of the deterministic system, while the ergodic invariant measures are attractors under some conditions for tumor growth...

A dynamic model is proposed to describe a mycoloop in aquatic food webs. The model consists of phytoplankton, chytrids and zooplankton. It characterizes that zooplankton consume both phytoplankton and free-living chytrid spores, and that chytrids infect phytoplankton. The dynamics of the model are investigated containing the dissipativity, existence and stability of equilibria, and persistence. The ecological reproductive indexes for phytoplankton or zooplankton invasion and basic reproduction numbers for chytrid transmission are derived...

The unprecedented scale and rapidity of dissemination of re-emerging and emerging infectious diseases impose new challenges for regulators and health authorities. To curb the dispersal of such diseases, proper management of healthcare facilities and vaccines are core drivers. In the present work, we assess the unified impact of healthcare facilities and vaccination on the control of an infectious disease by formulating a mathematical model. To formulate the model for any region, we consider four classes of human population; namely, susceptible, infected, hospitalized, and vaccinated...

The emergence and maintenance of tree species diversity in tropical forests is commonly attributed to the Janzen-Connell (JC) hypothesis, which states that growth of seedlings is suppressed in the proximity of conspecific adult trees. As a result, a JC distribution due to a density-dependent negative feedback emerges in the form of a (transient) pattern where conspecific seedling density is highest at intermediate distances away from parent trees. Several studies suggest that the required density-dependent feedbacks behind this pattern could result from interactions between trees and soil-borne pathogens...

Engaging in smoking not only leads to substantial health risks but also imposes considerable financial burdens. To deepen our understanding of the mechanisms behind smoking transmission and to address the tobacco epidemic, we examined a five-dimensional smoking epidemic model that accounts for different degrees of smoking under both deterministic and stochastic conditions. In the deterministic case, we determine the basic reproduction number, analyze the stability of equilibria with and without smoking, and investigate the existence of saddle-node bifurcation...

Research into the processes governing species richness has often assumed that the environment is fixed, whereas realistic environments are often characterised by random fluctuations over time. This temporal environmental stochasticity (TES) changes the demographic rates of species populations, with cascading effects on community dynamics and species richness. Theoretical and applied studies have used process-based mathematical models to determine how TES affects species richness, but under a variety of frameworks...

In this paper, a stoichiometric aquatic tri-trophic level model is proposed and analyzed, which incorporates the effect of light and phosphorus, as well as the fear effect in predator-prey interactions. The analysis of the model includes the dissipativity and the existence and stability of equilibria. The influence of environmental factors and fear effect on the dynamics of the system is particularly investigated. The key findings reveal that the coexistence of populations is positively influenced by an appropriate level of light intensity and/or the dissolved phosphorus input concentration; however, excessive levels of phosphorus input can disrupt the system, leading to chaotic behaviors...

Influenza virus infection can result in a wide range of clinical outcomes from asymptomatic infection to severe disease and death. While there are undoubtedly many factors that contribute to the severity of disease, one possible contributing factor that needs more investigation is the route of infection. In this study, we use previously published data from cynomolgus macaques infected with A/Vietnam/1203/04 (H5N1) via either aerosol (with and without bronchoalveolar lavages (BAL)) or a combined intrabronchial, oral, and intranasal route...

A common question in the aquatic sciences is that of how zooplankter movement can be modeled. It is well-established in the literature that there exists a randomness to this movement, but the question is how to characterize this randomness. The most common methods for doing this involve the random walk and correlated random walk (CRW) models. Here, we present a time series model that allows a better description the randomness in Daphnia motion when the amount of time that elapses between observations of their position is small...

Cell migration is one of the most studied phenomena in biology since it plays a fundamental role in many physiological and pathological processes such as morphogenesis, wound healing and tumorigenesis. In recent years, researchers have performed experiments showing that cells can migrate in response to mechanical stimuli of the substrate they adhere to. Motion toward regions of the substrate with higher stiffness is called durotaxis, while motion guided by the stress or the deformation of the substrate itself is called tensotaxis...

Avian influenza scare is a human psychological factor that asserts both positive and negative effects on the transmission of zoonotic avian influenza. In order to study the dichotomous effect of avian influenza scare on disease transmission, taking H7N9 avian influenza as a typical case, a two-patch epidemic model is proposed. The global dynamics and the threshold criteria are established by LaSalle invariant principle and the theory of asymptotic autonomous system. To mitigate the negative effects and curb illegal poultry trade, a game-theoretic model is adopted to explore the optimal policy of culling subsidies to reasonably compensate stakeholders for their economic losses resulting from the scare...

Microorganisms regulate the expression of energetically expensive phenotypes via a collective decision-making mechanism known as quorum sensing (QS). This study investigates the intricate dynamics of biofilm growth and QS-controlled biofilm dispersal in heterogeneous porous media, employing a pore-scale reactive transport modeling approach. Model simulations carried out under various fluid flow conditions and biofilm growth scenarios reveal that QS processes are influenced not only by the biomass density of biofilm colonies but also by a complex interplay between pore architecture, flow velocity, and the rates of biofilm growth and dispersal...