On a three-dimensional and two four-dimensional oncolytic viro-therapy models

We revisit here and carry out further works on tumor-virotherapy compartmental models of [Tian, 2011, Wang et al., 2013, Phan and Tian, 2017, Guo et al., 2019]. The results of these papers are only slightly pushed further. However, what is new is the fact that we make public our electronic notebooks, since we believe that easy electronic reproducibility is crucial in an era in which the role of the software becomes very important.Comment: 41 pages, 15 figure

Dynamical Models of Biology and Medicine

Mathematical and computational modeling approaches in biological and medical research are experiencing rapid growth globally. This Special Issue Book intends to scratch the surface of this exciting phenomenon. The subject areas covered involve general mathematical methods and their applications in biology and medicine, with an emphasis on work related to mathematical and computational modeling of the complex dynamics observed in biological and medical research. Fourteen rigorously reviewed papers were included in this Special Issue. These papers cover several timely topics relating to classical population biology, fundamental biology, and modern medicine. While the authors of these papers dealt with very different modeling questions, they were all motivated by specific applications in biology and medicine and employed innovative mathematical and computational methods to study the complex dynamics of their models. We hope that these papers detail case studies that will inspire many additional mathematical modeling efforts in biology and medicin

Mathematical Models of Androgen Resistance in Prostate Cancer Patients under Intermittent Androgen Suppression Therapy

abstract: Predicting resistant prostate cancer is critical for lowering medical costs and improving the quality of life of advanced prostate cancer patients. I formulate, compare, and analyze two mathematical models that aim to forecast future levels of prostate-specific antigen (PSA). I accomplish these tasks by employing clinical data of locally advanced prostate cancer patients undergoing androgen deprivation therapy (ADT). I demonstrate that the inverse problem of parameter estimation might be too complicated and simply relying on data fitting can give incorrect conclusions, since there is a large error in parameter values estimated and parameters might be unidentifiable. I provide confidence intervals to give estimate forecasts using data assimilation via an ensemble Kalman Filter. Using the ensemble Kalman Filter, I perform dual estimation of parameters and state variables to test the prediction accuracy of the models. Finally, I present a novel model with time delay and a delay-dependent parameter. I provide a geometric stability result to study the behavior of this model and show that the inclusion of time delay may improve the accuracy of predictions. Also, I demonstrate with clinical data that the inclusion of the delay-dependent parameter facilitates the identification and estimation of parameters.Dissertation/ThesisDoctoral Dissertation Applied Mathematics 201

Dynamics and precursor signs for phase transitions in neural systems

This thesis investigates neural state transitions associated with sleep, seizure and anaesthesia. The aim is to address the question: How does a brain traverse the critical threshold between distinct cortical states, both healthy and pathological? Specifically we are interested in sub-threshold neural behaviour immediately prior to state transition. We use theoretical neural modelling (single spiking neurons, a network of these, and a mean-field continuum limit) and in vitro experiments to address this question. Dynamically realistic equations of motion for thalamic relay neuron, reticular nuclei, cortical pyramidal and cortical interneuron in different vigilance states are developed, based on the Izhikevich spiking neuron model. A network of cortical neurons is assembled to examine the behaviour of the gamma-producing cortical network and its transition to lower frequencies due to effect of anaesthesia. Then a three-neuron model for the thalamocortical loop for sleep spindles is presented. Numerical simulations of these networks confirms spiking consistent with reported in vivo measurement results, and provides supporting evidence for precursor indicators of imminent phase transition due to occurrence of individual spindles. To complement the spiking neuron networks, we study the Wilson–Cowan neural mass equations describing homogeneous cortical columns and a 1D spatial cluster of such columns. The abstract representation of cortical tissue by a pair of coupled integro-differential equations permits thorough linear stability, phase plane and bifurcation analyses. This model shows a rich set of spatial and temporal bifurcations marking the boundary to state transitions: saddle-node, Hopf, Turing, and mixed Hopf–Turing. Close to state transition, white-noise-induced subthreshold fluctuations show clear signs of critical slowing down with prolongation and strengthening of autocorrelations, both in time and space, irrespective of bifurcation type. Attempts at in vitro capture of these predicted leading indicators form the last part of the thesis. We recorded local field potentials (LFPs) from cortical and hippocampal slices of mouse brain. State transition is marked by the emergence and cessation of spontaneous seizure-like events (SLEs) induced by bathing the slices in an artificial cerebral spinal fluid containing no magnesium ions. Phase-plane analysis of the LFP time-series suggests that distinct bifurcation classes can be responsible for state change to seizure. Increased variance and growth of spectral power at low frequencies (f < 15 Hz) was observed in LFP recordings prior to initiation of some SLEs. In addition we demonstrated prolongation of electrically evoked potentials in cortical tissue, while forwarding the slice to a seizing regime. The results offer the possibility of capturing leading temporal indicators prior to seizure generation, with potential consequences for understanding epileptogenesis. Guided by dynamical systems theory this thesis captures evidence for precursor signs of phase transitions in neural systems using mathematical and computer-based modelling as well as in vitro experiments

2020 Conference Abstracts: Annual Undergraduate Research Conference at the Interface of Biology and Mathematics

Schedule and abstract book for the Twelfth Annual Undergraduate Research Conference at the Interface of Biology and Mathematics Date: October 31 - November 1, 2020Location: The 2020 conference was conducted remotely due to COVID-19 concerns, utilizing the sococo platform that allows personal avatars to move between rooms and sessions, interact in small groups and also participate in zoom sessions.Keynote Speaker: Gerardo Chowell, Population Health Sciences, Georgia State Univ. School of Public Health, AtlantaFeatured Speaker: Olivia Prosper, Mathematics, Univ. of Tennessee, Knoxvill

Methods in Computational Biology

Modern biology is rapidly becoming a study of large sets of data. Understanding these data sets is a major challenge for most life sciences, including the medical, environmental, and bioprocess fields. Computational biology approaches are essential for leveraging this ongoing revolution in omics data. A primary goal of this Special Issue, entitled “Methods in Computational Biology”, is the communication of computational biology methods, which can extract biological design principles from complex data sets, described in enough detail to permit the reproduction of the results. This issue integrates interdisciplinary researchers such as biologists, computer scientists, engineers, and mathematicians to advance biological systems analysis. The Special Issue contains the following sections:•Reviews of Computational Methods•Computational Analysis of Biological Dynamics: From Molecular to Cellular to Tissue/Consortia Levels•The Interface of Biotic and Abiotic Processes•Processing of Large Data Sets for Enhanced Analysis•Parameter Optimization and Measuremen

Identifying stochastic oscillations in single-cell live imaging time series using Gaussian processes

Multiple biological processes are driven by oscillatory gene expression at different time scales. Pulsatile dynamics are thought to be widespread, and single-cell live imaging of gene expression has lead to a surge of dynamic, possibly oscillatory, data for different gene networks. However, the regulation of gene expression at the level of an individual cell involves reactions between finite numbers of molecules, and this can result in inherent randomness in expression dynamics, which blurs the boundaries between aperiodic fluctuations and noisy oscillators. Thus, there is an acute need for an objective statistical method for classifying whether an experimentally derived noisy time series is periodic. Here we present a new data analysis method that combines mechanistic stochastic modelling with the powerful methods of non-parametric regression with Gaussian processes. Our method can distinguish oscillatory gene expression from random fluctuations of non-oscillatory expression in single-cell time series, despite peak-to-peak variability in period and amplitude of single-cell oscillations. We show that our method outperforms the Lomb-Scargle periodogram in successfully classifying cells as oscillatory or non-oscillatory in data simulated from a simple genetic oscillator model and in experimental data. Analysis of bioluminescent live cell imaging shows a significantly greater number of oscillatory cells when luciferase is driven by a {\it Hes1} promoter (10/19), which has previously been reported to oscillate, than the constitutive MoMuLV 5' LTR (MMLV) promoter (0/25). The method can be applied to data from any gene network to both quantify the proportion of oscillating cells within a population and to measure the period and quality of oscillations. It is publicly available as a MATLAB package.Comment: 36 pages, 17 figure

Mathematical models for immunology:current state of the art and future research directions

The advances in genetics and biochemistry that have taken place over the last 10 years led to significant advances in experimental and clinical immunology. In turn, this has led to the development of new mathematical models to investigate qualitatively and quantitatively various open questions in immunology. In this study we present a review of some research areas in mathematical immunology that evolved over the last 10 years. To this end, we take a step-by-step approach in discussing a range of models derived to study the dynamics of both the innate and immune responses at the molecular, cellular and tissue scales. To emphasise the use of mathematics in modelling in this area, we also review some of the mathematical tools used to investigate these models. Finally, we discuss some future trends in both experimental immunology and mathematical immunology for the upcoming years