Mathematical Biosciences

Computational models of brain regions are crucial for understanding neuronal network dynamics and the emergence of cognitive functions. However, current supercomputing limitations hinder the implementation of large networks with millions of morphological and biophysical accurate neurons. Consequently, research has focused on simplified spiking neuron models, ranging from the computationally fast Leaky Integrate and Fire (LIF) linear models to more sophisticated non-linear implementations like Adaptive Exponential (AdEX) and Izhikevic models, through Generalized Leaky Integrate and Fire (GLIF) approaches...

We investigated a system of ordinary differential equations that describes the dynamics of prey and predator populations, taking into account the Allee effect affecting the reproduction of the predator population, and mutual interference amongst predators, which is modeled with the Bazykin-Crowley-Martin (BCM) trophic function. Bifurcation analysis revealed a rich spectrum of bifurcations occurring in the system. In particular, analytical conditions for the saddle-node, Hopf, cusp, and Bogdanov-Takens bifurcations were derived for the model parameters, quantifying the strength of the predator interference, the Allee effect, and the predation efficiency...

This paper proposes a bidimensional modeling framework for Wolbachia invasion, assuming imperfect maternal transmission, incomplete cytoplasmic incompatibility, and direct infection loss due to thermal stress. Our model adapts to various Wolbachia strains and retains all properties of higher-dimensional models. The conditions for the durable coexistence of Wolbachia-carrying and wild mosquitoes are expressed using the model's parameters in a compact closed form. When the Wolbachia bacterium is locally established, the size of the remanent wild population can be assessed by a direct formula derived from the model...

Antibiotics Time Machine is an important problem to understand antibiotic resistance and how it can be reversed. Mathematically, it can be modelled as follows: Consider a set of genotypes, each of which contain a set of mutated and unmutated genes. Suppose that a set of growth rate measurements of each genotype under a set of antibiotics are given. The transition probabilities of a 'realization' of a Markov chain associated with each arc under each antibiotic are computable via a predefined function given the growth rate realizations...

More than 60% of individuals recovering from substance use disorder relapse within one year. Some will resume drug consumption even after decades of abstinence. The cognitive and psychological mechanisms that lead to relapse are not completely understood, but stressful life experiences and external stimuli that are associated with past drug-taking are known to play a primary role. Stressors and cues elicit memories of drug-induced euphoria and the expectation of relief from current anxiety, igniting an intense craving to use again; positive experiences and supportive environments may mitigate relapse...

The mosquito-borne disease (malaria) imposes significant challenges on human health, healthcare systems, and economic growth/productivity in many countries. This study develops and analyzes a model to understand the interplay between malaria dynamics, economic growth, and transient events. It uncovers varied effects of malaria and economic parameters on model outcomes, highlighting the interdependence of the reproduction number (R0 ) on both malaria and economic factors, and a reciprocal relationship where malaria diminishes economic productivity, while higher economic output is associated with reduced malaria prevalence...

Metronomic chemotherapy refers to the frequent administration of chemotherapeutic agents at a lower dose and presents an attractive alternative to conventional chemotherapy with encouraging response rates. However, the schedule of the therapy, including the dosage of the drug, is usually based on empiricism. The confounding effects of tumor-endothelial-immune interactions during metronomic administration of drugs have not yet been explored in detail, resulting in an incomplete assessment of drug dose and frequency evaluations...

A basic mathematical model for IL-2-based cancer immunotherapy is proposed and studied. Our analysis shows that the outcome of therapy is mainly determined by three parameters, the relative death rate of CD4+ T cells, the relative death rate of CD8+ T cells, and the dose of IL-2 treatment. Minimal equilibrium tumor size can be reached with a large dose of IL-2 in the case that CD4+ T cells die out. However, in cases where CD4+ and CD8+ T cells persist, the final tumor size is independent of the IL-2 dose and is given by the relative death rate of CD4+ T cells...

Understanding the conditions for maintaining cooperation in groups of unrelated individuals despite the presence of non-cooperative members is a major research topic in contemporary biological, sociological, and economic theory. The N-person snowdrift game models the type of social dilemma where cooperative actions are costly, but there is a reward for performing them. We study this game in a scenario where players move between play groups following the casual group dynamics, where groups grow by recruiting isolates and shrink by losing individuals who then become isolates...

We have designed a stochastic model of embryonic neurogenesis in the mouse cerebral cortex, using the formalism of compound Poisson processes. The model accounts for the dynamics of different progenitor cell types and neurons. The expectation and variance of the cell number of each type are derived analytically and illustrated through numerical simulations. The effects of stochastic transition rates between cell types, and stochastic duration of the cell division cycle have been investigated sequentially. The model does not only predict the number of neurons, but also their spatial distribution into deeper and upper cortical layers...

We propose a continuum model for pattern formation, based on the multiphase model framework, to explore in vitro cell patterning within an extracellular matrix. We demonstrate that, within this framework, chemotaxis-driven cell migration can lead to the formation of cell clusters and vascular-like structures in 1D and 2D respectively. The influence on pattern formation of additional mechanisms commonly included in multiphase tissue models, including cell-matrix traction, contact inhibition, and cell-cell aggregation, are also investigated...

We use a compartmental model with a time-varying transmission parameter to describe county level COVID-19 transmission in the greater St. Louis area of Missouri and investigate the challenges in fitting such a model to time-varying processes. We fit this model to synthetic and real confirmed case and hospital discharge data from May to December 2020 and calculate uncertainties in the resulting parameter estimates. We also explore non-identifiability within the estimated parameter set. We determine that that death rate of infectious non-hospitalized individuals, the testing parameter and the initial number of exposed individuals are not identifiable based on an investigation of correlation coefficients between pairs of parameter estimates...

The collective foraging behavior of ant colonies is a central focus in behavioral ecology. This paper enhances the classical model of foraging dynamics in harvester ant colonies by introducing a nonlinear recruitment rate and considering environmental variability. Initially, we analyze the existence and stability of steady states in the deterministic model. The results suggest that an increase in mean recruitment time can reduce the foraging threshold, leading to both forward and backward bifurcations. Furthermore, both average recruitment time and the interference intensity of recruiters impact the number of workers in each subgroup...

Efficient and accurate large-scale networks are a fundamental tool in modelling brain areas, to advance our understanding of neuronal dynamics. However, their implementation faces two key issues: computational efficiency and heterogeneity. Computational efficiency is achieved using simplified neurons, whereas there are no practical solutions available to solve the problem of reproducing in a large-scale network the experimentally observed heterogeneity of the intrinsic properties of neurons. This is important, because the use of identical nodes in a network can generate artifacts which can hinder an adequate representation of the properties of a real network...

The development or remission of diet-induced overweight type 2 diabetes involves many biological changes which occur over very different timescales. Remission, defined by HbA1c <6.5%, or fasting plasma glucose concentration G<126 mg/dl, may be achieved rapidly by following weight loss guidelines. However, remission is often short-term, followed by relapse. Mathematical modelling provides a way of investigating a typical situation, in which patients are advised to lose weight and then maintain fat mass, a slow variable...

Interactions between SARS-CoV-2 and the immune system during infection are complex. However, understanding the within-host SARS-CoV-2 dynamics is of enormous importance for clinical and public health outcomes. Current mathematical models focus on describing the within-host SARS-CoV-2 dynamics during the acute infection phase. Thereby they ignore important long-term post-acute infection effects. We present a mathematical model, which not only describes the SARS-CoV-2 infection dynamics during the acute infection phase, but extends current approaches by also recapitulating clinically observed long-term post-acute infection effects, such as the recovery of the number of susceptible epithelial cells to an initial pre-infection homeostatic level, a permanent and full clearance of the infection within the individual, immune waning, and the formation of long-term immune capacity levels after infection...

Drug resistance is one of the most intractable issues to the targeted therapy for cancer diseases. To explore effective combination therapy schemes, we propose a mathematical model to study the effects of different treatment schemes on the dynamics of cancer cells. Then we characterize the dynamical behavior of the model by finding the equilibrium points and exploring their local stability. Lyapunov functions are constructed to investigate the global asymptotic stability of the model equilibria. Numerical simulations are carried out to verify the stability of equilibria and treatment outcomes using a set of collected model parameters and experimental data on murine colon carcinoma...

Substance use disorder (SUD) is a complex disease involving nontrivial biological, psychological, environmental, and social factors. While many mathematical studies have proposed compartmental models for SUD, almost all of these exclusively model new cases as the result of an infectious process, neglecting any SUD that was primarily developed in social isolation. While these decisions were likely made to facilitate mathematical analysis, isolated SUD development is critical for the most common substances of abuse today, including opioid use disorder developed through prescription use and alcoholism developed primarily due to genetic factors or stress, depression, and other psychological factors...

We here propose a hybrid computational framework to reproduce and analyze aspects of the avascular progression of a generic solid tumor. Our method first employs an individual-based approach to represent the population of tumor cells, which are distinguished in viable and necrotic agents. The active part of the disease is in turn differentiated according to a set of metabolic states. We then describe the spatio-temporal evolution of the concentration of oxygen and of tumor-secreted proteolytic enzymes using partial differential equations (PDEs)...

The sterile insect technique (SIT) is a technique to control pests and vectors of diseases by releasing mainly sterile males. Several challenges need to be solved before large-scale field application in order to guarantee its success. In this paper we intend to focus on two important issues: residual fertility in released (sterile) males and contamination of each release by sterile females. Indeed, sterile males are never 100% sterile, that is there is always a small proportion, ɛ, of fertile males (sperm of) within the sterile males population...