Michiel Bertsch;Bruno Franchi;Norina Marcello;Maria Carla Tesi;Andrea Tosin
{"title":"Alzheimer's disease: a mathematical model for onset and progression","authors":"Michiel Bertsch;Bruno Franchi;Norina Marcello;Maria Carla Tesi;Andrea Tosin","doi":"10.1093/imammb/dqw003","DOIUrl":"10.1093/imammb/dqw003","url":null,"abstract":"In this article we propose a mathematical model for the onset and progression of Alzheimer's disease based on transport and diffusion equations. We regard brain neurons as a continuous medium and structure them by their degree of malfunctioning. Two different mechanisms are assumed to be relevant for the temporal evolution of the disease: i) diffusion and agglomeration of soluble polymers of amyloid, produced by damaged neurons and ii) neuron-to-neuron prion-like transmission. We model these two processes by a system of Smoluchowski equations for the amyloid concentration, coupled to a kinetic-type transport equation for the distribution function of the degree of malfunctioning of neurons. The second equation contains an integral term describing the random onset of the disease as a jump process localized in particularly sensitive areas of the brain. Our numerical simulations are in good qualitative agreement with clinical images of the disease distribution in the brain which vary from early to advanced stages.","PeriodicalId":94130,"journal":{"name":"Mathematical medicine and biology : a journal of the IMA","volume":"34 2","pages":"193-214"},"PeriodicalIF":0.0,"publicationDate":"2017-06-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://sci-hub-pdf.com/10.1093/imammb/dqw003","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"34461501","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Marc Artzrouni;Vasiliy N. Leonenko;Thierry A. Mara
{"title":"A syringe-sharing model for the spread of HIV: application to Omsk, Western Siberia","authors":"Marc Artzrouni;Vasiliy N. Leonenko;Thierry A. Mara","doi":"10.1093/imammb/dqv036","DOIUrl":"10.1093/imammb/dqv036","url":null,"abstract":"A system of two differential equations is used to model the transmission dynamics of human immunodeficiency virus between ‘persons who inject drugs’ (PWIDs) and their syringes. Our vector-borne disease model hinges on a metaphorical urn from which PWIDs draw syringes at random which may or may not be infected and may or may not result in one of the two agents becoming infected. The model's parameters are estimated with data mostly from the city of Omsk in Western Siberia. A linear trend in PWID prevalence in Omsk could only be fitted by considering a time-dependent version of the model captured through a secular decrease in the probability that PWIDs decide to share a syringe. A global sensitivity analysis is performed with 14 parameters considered random variables in order to assess their impact on average numbers infected over a 50-year projection. With obvious intervention implications the drug injection rate and the probability of syringe-cleansing are the only parameters whose coefficients of correlations with numbers of infected PWIDs and infected syringes have an absolute value close to or larger than 0.40.","PeriodicalId":94130,"journal":{"name":"Mathematical medicine and biology : a journal of the IMA","volume":"34 1","pages":"15-37"},"PeriodicalIF":0.0,"publicationDate":"2017-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://sci-hub-pdf.com/10.1093/imammb/dqv036","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"34107572","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Insights from mathematical modelling for T cell migration into the central nervous system","authors":"T. Ruck;S. Bittner;S. G. Meuth;M. Herty","doi":"10.1093/imammb/dqv038","DOIUrl":"10.1093/imammb/dqv038","url":null,"abstract":"The migration of immune cells from peripheral immune organs into the central nervous system (CNS) through the blood–brain barrier (BBB) is a tightly regulated process. The complex interplay between cells of the BBB and immune cells coordinates cell migration as a part of normal immune surveillance while its dysregulation is critically involved in the pathogenesis of various CNS diseases. To develop tools for a deeper understanding of distribution and migratory pattern of immune cells regulated by the BBB, we made use of a mathematical modelling approach derived from Markov chain theory. We present a data-driven model using a derivation of kinetic differential equations from a particle game. According to the theory of gases, these equations allow one to predict the mean behaviour of a large class of cells by modelling cell–cell interactions. We used this model to assess the distribution of naive, central memory and effector memory T lymphocytes in the peripheral blood and cerebrospinal fluid. Our model allows us to evaluate the impact of activation status, migratory capacity and cell death for cell distribution in the peripheral blood and the CNS.","PeriodicalId":94130,"journal":{"name":"Mathematical medicine and biology : a journal of the IMA","volume":"34 1","pages":"39-58"},"PeriodicalIF":0.0,"publicationDate":"2017-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://sci-hub-pdf.com/10.1093/imammb/dqv038","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"34200032","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"What mechanisms of tau protein transport could be responsible for the inverted tau concentration gradient in degenerating axons?","authors":"I. A. Kuznetsov;A. V. Kuznetsov","doi":"10.1093/imammb/dqv041","DOIUrl":"https://doi.org/10.1093/imammb/dqv041","url":null,"abstract":"In tauopathies, such as Alzheimer's disease (AD), microtubule (MT)-associated protein tau detaches from MTs and aggregates, eventually forming insoluble neurofibrillary tangles. In a healthy axon, the tau concentration increases toward the axon terminal, but in a degenerating axon, the tau concentration gradient is inverted and the highest tau concentration is in the soma. In this article, we developed a mathematical model of tau transport in axons. We calibrated and tested the model by using published distributions of tau concentration and tau average velocity in a healthy axon. According to published research, the inverted tau concentration gradient may be one of the reasons leading to AD. We therefore used the model to investigate what modifications in tau transport can lead to the inverted tau concentration gradient. We investigated whether tau detachment from MTs due to tau hyperphosphorylation can cause the inverted tau concentration gradient. We found that the assumption that most tau molecules are detached from MTs does not consistently predict the inverted tau concentration gradient; the predicted tau distribution becomes more uniform if the axon length is increased. We then hypothesized that in degenerating axons some tau remains bound to MTs and participates in the component ‘a’ of slow axonal transport but that the rate of tau reversals from anterograde to retrograde motion increases. We demonstrated that this hypothesis results in a tau distribution where the tau concentration has its maximum value at the axon hillock and its minimum value at the axon terminal, in agreement with what is observed in AD. Our results thus suggest that defects in active transport of tau may be a contributing factor to the onset of neural degeneration.","PeriodicalId":94130,"journal":{"name":"Mathematical medicine and biology : a journal of the IMA","volume":"34 1","pages":"125-150"},"PeriodicalIF":0.0,"publicationDate":"2017-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://sci-hub-pdf.com/10.1093/imammb/dqv041","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"49992393","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Sensitivity analysis for dose deposition in radiotherapy via a Fokker–Planck model","authors":"Richard C. Barnard;Martin Frank;Kai Krycki","doi":"10.1093/imammb/dqv039","DOIUrl":"https://doi.org/10.1093/imammb/dqv039","url":null,"abstract":"In this paper, we study the sensitivities of electron dose calculations with respect to stopping power and transport coefficients. We focus on the application to radiotherapy simulations. We use a Fokker–Planck approximation to the Boltzmann transport equation. Equations for the sensitivities are derived by the adjoint method. The Fokker–Planck equation and its adjoint are solved numerically in slab geometry using the spherical harmonics expansion (\u0000<tex>$P_N$</tex>\u0000) and an Harten-Lax-van Leer finite volume method. Our method is verified by comparison to finite difference approximations of the sensitivities. Finally, we present numerical results of the sensitivities for the normalized average dose deposition depth with respect to the stopping power and the transport coefficients, demonstrating the increase in relative sensitivities as beam energy decreases. This in turn gives estimates on the uncertainty in the normalized average deposition depth, which we present.","PeriodicalId":94130,"journal":{"name":"Mathematical medicine and biology : a journal of the IMA","volume":"34 1","pages":"109-123"},"PeriodicalIF":0.0,"publicationDate":"2017-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://sci-hub-pdf.com/10.1093/imammb/dqv039","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"49992394","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A biochemical and mechanical model of injury-induced intimal thickening","authors":"Pak-Wing Fok;Rebecca Sanft","doi":"10.1093/imammb/dqv040","DOIUrl":"https://doi.org/10.1093/imammb/dqv040","url":null,"abstract":"In this paper, we investigate an axisymmetric model of intimal thickening using hyperelasticity theory. Our model describes the growth of the arterial intima due to cell proliferation which, in turn, is driven by the release of a cytokine such as platelet-derived growth factor (PDGF). With the growth rate tied to both local stress and the local concentration of PDGF, we derive a quadruple free boundary problem with different regions of the vessel wall characterized by different homeostatic stress. We compare our model predictions to rabbit and rodent models of atherosclerosis and find that in order to achieve the growth rates reported in the experiments, growth must be mainly cytokine induced rather than stress induced. Our model is also able to reproduce Glagov remodelling where, as a vessel becomes more diseased, the lumen expands before rapidly contracting.","PeriodicalId":94130,"journal":{"name":"Mathematical medicine and biology : a journal of the IMA","volume":"34 1","pages":"77-108"},"PeriodicalIF":0.0,"publicationDate":"2017-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://sci-hub-pdf.com/10.1093/imammb/dqv040","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"49992395","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Lydia Hill;Mark A. J. Chaplain;Roland Wolf;Yury Kapelyukh
{"title":"The usage of a three-compartment model to investigate the metabolic differences between hepatic reductase null and wild-type mice","authors":"Lydia Hill;Mark A. J. Chaplain;Roland Wolf;Yury Kapelyukh","doi":"10.1093/imammb/dqv029","DOIUrl":"10.1093/imammb/dqv029","url":null,"abstract":"The Cytochrome P450 (CYP) system is involved in 90% of the human body's interactions with xenobiotics and due to this, it has become an area of avid research including the creation of transgenic mice. This paper proposes a three-compartment model which is used to explain the drug metabolism in the Hepatic Reductase Null (HRN) mouse developed by the University of Dundee (Henderson, C. J., Otto, D. M. E., Carrie, D., Magnuson, M. A., McLaren, A. W., Rosewell, I. and Wolf, C. R. (2003) Inactivation of the hepatic cytochrome p450 system by conditional deletion of hepatic cytochrome p450 reductase. J. Biol. Chem. 278, 13480–13486). The model is compared with a two-compartment model using experimental data from studies using wild-type and HRN mice. This comparison allowed for metabolic differences between the two types of mice to be isolated. The three sets of drug data (Gefitinib, Midazolam and Thalidomide) showed that the transgenic mouse has a decreased rate of metabolism.","PeriodicalId":94130,"journal":{"name":"Mathematical medicine and biology : a journal of the IMA","volume":"34 1","pages":"1-13"},"PeriodicalIF":0.0,"publicationDate":"2017-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://sci-hub-pdf.com/10.1093/imammb/dqv029","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"34236577","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A non-local model for cancer stem cells and the tumour growth paradox","authors":"I. Borsi;A. Fasano;M. Primicerio;T. Hillen","doi":"10.1093/imammb/dqv037","DOIUrl":"https://doi.org/10.1093/imammb/dqv037","url":null,"abstract":"The tumour growth paradox refers to the observation that incomplete treatment of cancers can enhance their growth. As shown here and elsewhere, the existence of cancer stem cells (CSCs) can explain this effect. CSC are less sensitive to treatments, hence any stress applied to the tumour selects for CSC, thereby increasing the fitness of the tumour. In this paper, we use a mathematical model to understand the role of CSC in the progression of cancer. Our model is a rather general system of integro-differential equations for tumour growth and tumour spread. Such a model has never been analysed, and we prove results on local and global existence of solutions, their uniqueness and their boundedness. We show numerically that this model exhibits the tumour growth paradox for all parameters tested. This effect becomes more relevant for small renewal rate of the CSC.","PeriodicalId":94130,"journal":{"name":"Mathematical medicine and biology : a journal of the IMA","volume":"34 1","pages":"59-75"},"PeriodicalIF":0.0,"publicationDate":"2017-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://sci-hub-pdf.com/10.1093/imammb/dqv037","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"49992396","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Uncovering the natural history of cancer from post-mortem cross-sectional diameters of hepatic metastases","authors":"Leonid Hanin;Jason Rose","doi":"10.1093/imammb/dqv026","DOIUrl":"10.1093/imammb/dqv026","url":null,"abstract":"We develop a mathematical and statistical methodology for estimation of important unobservable characteristics of the individual natural history of cancer from a sample of cross-sectional diameters of liver metastases measured at autopsy. Estimation of the natural history of cancer is based on a previously proposed stochastic model of cancer progression tailored to this type of observations. The model accounts for primary tumour growth, shedding of metastases, their selection, latency and growth in a given secondary site. The model was applied to the aforementioned data on 428 liver metastases detected in one untreated small cell lung cancer patient. Identifiable model parameters were estimated by the method of maximum likelihood and through minimizing the \u0000<tex>$L^{2}$</tex>\u0000 distance between theoretical and empirical cumulative distribution functions. The model with optimal parameters provided an excellent fit to the data. Results of data analysis support, if only indirectly, the hypothesis of the existence of stem-like cancer cells in the case of small cell lung carcinoma and point to the possibility of suppression of metastatic growth by a large primary tumour. They also lead to determination of the lower and upper bounds for the age of cancer onset and expected duration of metastatic latency. Finally, model-based inference on the patient's natural history of cancer allowed us to conclude that resection of the primary tumour would most likely not have had a curative effect.","PeriodicalId":94130,"journal":{"name":"Mathematical medicine and biology : a journal of the IMA","volume":"33 4","pages":"397-416"},"PeriodicalIF":0.0,"publicationDate":"2016-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://sci-hub-pdf.com/10.1093/imammb/dqv026","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"33892862","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A fluid dynamics multidimensional model of biofilm growth: stability, influence of environment and sensitivity","authors":"F. Clarelli;C. Di Russo;R. Natalini;M. Ribot","doi":"10.1093/imammb/dqv024","DOIUrl":"10.1093/imammb/dqv024","url":null,"abstract":"In this article, we study in detail the fluid dynamics system proposed in Clarelli et al. (2013, J. Math. Biol., 66, 1387–1408) to model the formation of cyanobacteria biofilms. After analysing the linear stability of the unique non-trivial equilibrium of the system, we introduce in the model the influence of light and temperature, which are two important factors for the development of a cyanobacteria biofilm. Since the values of the coefficients we use for our simulations are estimated through information found in the literature, some sensitivity and robustness analyses on these parameters are performed. All these elements enable us to control and to validate the model we have already derived and to present some numerical simulations in the 2D and the 3D cases.","PeriodicalId":94130,"journal":{"name":"Mathematical medicine and biology : a journal of the IMA","volume":"33 4","pages":"371-395"},"PeriodicalIF":0.0,"publicationDate":"2016-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://sci-hub-pdf.com/10.1093/imammb/dqv024","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"34300927","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}