Dracula is a joint research team between Inria, Université Claude Bernard Lyon 1 (UCBL) and CNRS (Institut Camille-Jordan (ICJ, UMR 5208) and Laboratoire de Biologie et Modélisation de la Cellule (LBMC, UMR 5239)).

The Dracula project is devoted to multi-scale modeling in biology and medicine, and more specifically to the development of tools and methods to describe multi-scale processes in biology and medicine. Applications include normal and pathological hematopoiesis (for example leukemia), immune response, and other biological processes, like: tissue renewal, morphogenesis, atherosclerosis, prion disease, hormonal regulation of food intake, and so on. Multi-scale modeling implies simultaneous modeling of several levels of descriptions of biological processes: intra-cellular networks (molecular level), cell behavior (cellular level), dynamics of cell populations (organ or tissue) with the control by other organs (organism) (see Figure ). Such modeling represents one of the major challenges in modern science due to its importance and because of the complexity of biological phenomena and of the presence of very different interconnected scales.

Scheme of multi-scale models of cell dynamics

Although multi-scale modeling holds a great potential for biology and medicine, and despite the fact that a variety of techniques exists to deal with such problems, the complexity of the systems poses new challenges and needs the development of new tools. Moreover, different biological questions usually require different types of multi-scale modeling. The expected results of these studies are numerous. On one hand, they will shed new light on the understanding of specific biological and medical questions (for instance, what is the behavior of hematopoietic stem cells under pathological conditions? Or how to efficiently stimulate an immune response in order to design new vaccines?). On the other hand, the modeling methods developed here for specific processes are relevant to study other complex biological systems. We pay a special attention on developing methods that are not restricted to one or two applications.

An important part of our researches is performed in close collaboration with biologists and physicians in order to stay in contact with the biological and medical goals. The presence, within the project, of a biologist (Olivier Gandrillon) who has acquired over the years the know-how required for interacting with mathematicians is probably one of the main assets of the project. He participates actively in many tasks of our program, stimulates interactions between members of the project and biologists, and everyone benefits from his expertise in molecular and cell biology.


Multi-scale modeling; Hybrid modeling; Mathematical Biology; Computational Biology; Immune response modeling; Normal and pathological hematopoiesis; Multi-scale cancer modeling; Regulatory networks; Reaction-diffusion equation; Structured partial differential equations; Delay differential equations; Agent-based modeling; Dynamical systems.

Research axis 1: Mathematical modeling for cell population dynamics

Executive summary

Stem cells are essential for development and keep the maintenance of many tissues homeostasis. They are characterized by their ability to self-renew as well as to produce differentiated cells. They vary enormously, for each organ, in their proliferation capacity, their potency to produce different cell lineage and their response to various environmental cues. How a cell will react to a given external signal does not depend only on its current state but also on its environment. Understanding the effect of cell-to-cell heterogeneity and the spatial organization of cell populations is therefore necessary to help keeping the normal function of an organ.

We develop mathematical tools and methods to study cell population dynamics and other biological processes: stability of steady sates, existence of bifurcations, kinetic properties, spatial organization, in finely detailed cell populations. The main tools we use are hybrid discrete-continuous models, reaction-diffusion equations, structured models (in which the population is endowed with relevant structures or traits), delay differential systems, agent-based models. Our team has acquired an international expertise in the fields of analysis of reaction-diffusion and structured equations, particularly integro-differential and delay differential equations.

The mathematical methods we develop are not restricted to hematopoietic system (Research axis 2), and immune response (Research axis 3), rather we apply them in many other biological phenomena, for example: tissue renewal, morphogenesis, prion disease, atherosclerosis, hormonal regulation of food intake, cancer, and others.

Project-team positioning

The focus of this objective is the development, analysis and application of hybrid discrete-continuous, reaction-diffusion and structured partial differential models. The structured equations allow a fine description of a population as some structures (age, maturity, intracellular content) change with time. In many cases, structured equations can be partially integrated to yield integro-differential equations (ordinary or partial differential equations involving non-local integral terms), time-delay differential or time-delay partial differential, or coupled differential-difference models. Analysis of integro-differential and time-delay systems deals with existence of solutions and their stability. Applications are found in the study of normal and pathological hematopoietic system (Research axis 2), immune response (Research axis 3), morphogenesis, prion disease, cancer development and treatment, and generally in tissue renewal problems. Models based on structured equations are especially useful to take into account the effect of finite time cells take to divide, die or become mature. Reaction-diffusion equations are used in order to describe spatial distribution of cell populations. It is a well developed area of research in our team which includes qualitative properties of travelling waves for reaction-diffusion systems with or without delay, and complex nonlinear dynamics.

Our team has developed a solid expertise in mathematical analysis of reaction-diffusion with or without delay and structured equations (in particular, delay differential equations) and one of the most prolific. Other major groups are the teams of Benoit Perthame (Pierre et Marie CURIE University and Mamba, Paris,, Emmanuel Grenier (Ecole normale supérieure de Lyon and NUMED,, Odo Diekmann (Utrecht University, The Netherlands,, Avner Friedman (The Ohio State University, USA,, Jianhong Wu (York University, Canada,, Glenn Webb (Vanderbilt University, Nashville, USA,, Philip K. Maini (University of Oxford, England,, Mark Chaplain (University of St Andrews, Scotland,, Nicola Bellomo (University of Turin, Italy, Most of the members of all these groups and of our team belong to the same mathematical community working on partial differential equations and dynamical systems with applications to biology and medicine.


  •  University of Toronto, Canada; Mathematical analysis and applications of reaction-diffusion equations (more than 30 joint papers).

  •  Institute of Problems of Mechanical Engineering, St.Petersburg, Russia; Dynamics of cell renewal (more than 10 joint papers).

  •  Department of Cell and Molecular Biology and Department of Forensic Medicine, Stockholm, Sweden; Dynamics of cell generation and turnover (3 joint papers).

  •  Universities of Tlemcen (Algeria) and Marrakech (Morocco); Delay differential equations (7 joint papers)

Research axis 2: Multi-scale modeling of hematopoiesis and leukemia

Executive summary

Hematopoiesis is a complex process that begins with hematopoietic stem cells (HSCs) and results in formation of mature cells: red blood cells, white cells and platelets. Blood cells are produced in the bone marrow, from where mature cells are released into the blood stream. Hematopoiesis is based on a balance between cell proliferation (including self-renewal), differentiation and apoptosis. The choice between these three possibilities is determined by intra-cellular regulatory networks and by numerous control mechanisms in the bone marrow or carried out by other organs. Intra-cellular regulatory networks are complex biochemical reactions involving proteins, enzymes and signalling molecules. The deregulation of hematopoiesis can result in numerous blood diseases including leukemia (a cancer of blood cells). One important type of leukemia is Chronic Myeloid Leukemia (CML). The strong tyrosine kinase activity of the BCR-ABL protein is the basis for the main cell effects that are observed in CML: significant proliferation, anti-apoptotic effect, disruption of stroma adhesion properties. This explains the presence in CML blood of a very important number of cells belonging to the myeloid lineage, at all stages of maturation.

Multi-scale modeling in hematopoiesis holds a great potential. A variety of techniques exists to deal with this problem. However, the complexity of the system poses new difficulties and leads to the development of new tools. The expected results of this study are numerous. On one hand, it will shed new light on the different physiological mechanisms that converge toward the continuous regeneration of blood cells, for example: the understanding of deregulation of erythropoiesis (the process of red blood cell production) under drug treatments (this can lead to lack of red blood cells (anemia), or a surplus of red blood cells), the dynamic of leukemic cells under the action of drugs and the control of their resistance to these treatments.

Project team positioning

Multi-scale modeling of hematopoiesis is one of the key points of the project that has started in the early stage of the Dracula team. Investigated by all the team members, it took many years of close discussion with biologists to get the best understanding of the key role played by the most important molecules, hormones, kinase cascade, cell communication up to the latest knowledge. One of the important questions here is to identify particular biological mechanisms (intracellular regulation, control mechanisms) and to integrate them in the different models. Our main work consisted in the development of a hybrid (continuous/discrete) model for red blood cell progenitor proliferation, survival/death, differentiation, and migration. Cells are modeled as discrete objects, and the extracellular medium is described by continuous equations for extracellular concentrations. This is to our knowledge the most complete model for erythropoiesis to date, and the only one using a multi-scale formalism. Other models published by our group and others for hematopoiesis are population-based models, mostly population structured equations (transport partial differential equations or delay differential equations). The interest in modeling hematopoiesis dates back to the 70’s and two groups have been responsible for most of development in the past 40 years: Markus Loeffer’s team in Leipzig, Germany (Wichmann et al. 1976, in Mathematical Models in Medicine) and Michael Mackey’s team at McGill University, Montreal, Canada (Mackey 1978, Blood). Our model differs from population based models in that the regulation is directly modeled at the molecular level (See Figure ) rather than acting on rates at the population level. Thus we can take into account non-predictable effects of interactions between different molecular pathways and between cells that would otherwise be lost in the global population rates.

Regarding modeling leukemia, we concentrated on Chronic Myeloid Leukemia (CML) and its treatment. We considered models based on ordinary differential equations for the action of the main proteins involved in CML (as BCR-ABL protein), and of transport equations (with or without delay, physiologically structured or not) to represent healthy and leukemic cell populations, take into account many interactions between proteins (especially BCR-ABL), cells (anti-apoptotic effect, etc.). The development of models for CML allowed us to interact with Franck Nicolini in Lyon (Centre Hospitalier de Lyon) and Doron Levy (Maryland University, Different schools developed models for CML and its treatment. The three leading groups are the ones of Franziska Michor (Harvard School of public health,, Ingo Roeder (Institute for Medical Informatics and Biometry, Dresden, and Michael Mackey (McGill University,


Members of the team have worked for several years in collaboration with biologists (François Morlé, University Lyon 1) and hematologists (Charles Dumontet, Lyon and Mark Koury, Nashville, on the Modelling of normal and pathological hematopoiesis .

The work on modeling Leukemia is based on two major collaborations: firstly, an ongoing (since 2011) mathematical collaboration with the University of Maryland through the program Associate Teams Inria project, “Modelling Leukemia” ( Secondly, an ongoing (since 2012) collaboration with a clinician from Hospices Civils de Lyon (Dr. F.E. Nicolini). In this framework, we shall have soon access to the data of the clinical trial PETALs ( 2×100 patients).

Research axis 3: Multi-scale modeling of the immune response

Executive summary

Vaccination represents a worldwide health, social and economical challenge as it has allowed the eradication or the strong containment of several devastating diseases over the past century. However to date, most of the effective vaccines rely on the generation of neutralizing antibody responses and such vaccines have proven largely unsuccessful in the prevention against some pathogens, such as HIV or malaria. In such cases, vaccines geared towards the generation of CD8 T cell immunity may provide a better protection. The generation of memory CD8 T cells following antigenic immunization is a long process (lasting up to month in murine preclinical models), therefore strongly slowing the process of vaccine monitoring in preclinical studies. Thus, the dynamical modeling of the CD8 T cell immune response both at the cellular and molecular levels should provide an important tool to better understand the dynamics of the response and to speed-up the process and reduce costs of vaccine development.

However, currently published cellular models of the immune response are either over-simplified, not predicting important parameters of this response, or too complicated for most of their parameters to be accessible for experimental measurements, thus impeding their biological validation. Dynamical models of the CD8 T cell response at the molecular level are very scarce and there is no multi-scale model of the immune response giving insights into both the regulation at the molecular scale and the consequences on cell population dynamics.

The objective of this research axis is therefore to develop a predictive multi-scale model of the CD8 T cell response, by confronting the model at different stages to in vivo-acquired experimental data, in order to be able to investigate the influence of early molecular events on cell population dynamics few days or weeks later.

Project-team positioning

We are aiming at building and analyzing a multi-scale model of the CD8 T cell immune response, from the molecular to the cellular and potentially organismal scale. This consists in describing the dynamics at each scale with relevant formalisms as well as the careful description of the couplings between scales.

Only few research groups are actually working on the CD8 T cell immune response around the world, and none of them deals with multi-scale modeling of this response. A network developed around Alan Perelson’s work in theoretical immunology in the last decades, at Los Alamos National Laboratory, and involves mainly people in various US universities or institutes. In Europe, Rob De Boer’s group ( of theoretical immunology in Utrecht, Netherlands, is the historical leader in the CD8 T cell dynamics modeling. We considered the models developed in these groups when we started our project, and we contributed to improve them by using nonlinearities accounting for cell population interactions to regulate the response. Also, our initial focus was on the generation of memory cells associated with vaccine development so we modeled CD8 T cell responses against influenza and vaccinia viruses, whereas other groups usually consider LCMV in its chronic form.

Ron Germain’s group at the NIH, and Grégoire Altan-Bonnet in subsequent works, focused on the molecular regulation of the CD4 and CD8 T cell immune responses. In particular, they built the Simmune software, which allows the modeling and simulation of molecular interactions ( This software is not really devoted to multi-scale modeling yet it provides an interesting tool to describe molecular interactions. Since our aim is to couple molecular and cellular scales at the tissue level, and we do not want to consider large networks but rather small-simplified informative interaction networks, we are confident that our approach is complementary of these works.

Within Inria project-teams, NUMED develops multi-scale approaches for biological problems, and MAMBA and MONC ( mention models of cancer progression and treatment including immune responses. In the first case the methodology is similar, and collaborations between NUMED and DRACULA already exist (both teams are located in Lyon), but applications differ. In the second case, MAMBA and MONC are mainly focused on cancer modeling and up to now are motivated by including an action of the immune system in the fight against cancer, which is very different from what we are developing. However, both modeling approaches are complementary and could lead to interactions, in particular in the light of recent advances in medical research pointing towards an important role – and high expectations – of the immune reaction in fighting cancers. Finally, SISTM ( also focuses on the modeling of the immune response, mainly against HIV, but the motivation is very similar to ours: the objective is to provide tools and methods in order to efficiently develop vaccines. They consider the CD4 T cell response instead of the CD8 T cell response, and biostatistics to achieve their goals instead of multi-scale models, yet even though there is no interaction between SISTM and DRACULA at this moment our methods and objectives are close enough to foreshadow future collaborations.


On this topic our main collaborators are members of Jacqueline Marvel’s team in Lyon in the CIRI (Centre International de Recherche en Infectiologie INSERM U1111): Dr. Jacqueline Marvel, head of the team, Dr. Christophe Arpin (CR CNRS), and other technicians and engineers of the team. They are all immunologists, specialists of the CD8 T cell response and of the generation of memory CD8 T cells.

We also interact with private companies: AltraBio (, that provides tools for data analysis, and CosmoTech, that develops a modeling and simulating platform that should allow transferring our model on an easy-to-use platform devoted to commercial uses.

Evolution of research direction during the last evaluation

Reminder of the objectives given for the last evaluation

The aim of this project is the development of modern tools for multi-scale modeling in biological phenomena. During the period 2014-2017, the objectives we had fixed were to develop modern tools for multi-scale modeling of biological phenomena, as detailed hereafter:

  1. Multi-scale modeling of erythropoiesis, the process of red blood cell production, in order to describe normal, stress, and pathological erythropoiesis, using mathematical and computational models. This led to:

  2. The modeling of hemoglobin instability in dialysis patients: Thomas Lepoutre has been progressively taking part in this theme through a collaboration with P. Kim (University of Sydney, Australia);

  3. Multi-scale modeling of the CD8 T cell immune response, in order to develop a predictive model of the CD8 T cell response, by confronting the model at different stages to in vivo-acquired experimental data;

  4. Population dynamics modeling, with the aim to develop general mathematical tools to study them. The main tools we were using were structured equations, in which the cell population is endowed with relevant structures, or traits. We identified limitations in using these formalisms, this is why we started developing multi-scale approaches;

  5. Modeling of Chronic Myeloid Leukemia (CML) treatment, using ordinary differential equations models. Our team had already developed a first model of mutant leukemic cells being resistant to chemotherapy. A next step would be to identify the parameters using experimental data;

  6. Multi-scale modeling carried out on the basis of hybrid discrete-continuous models, where dissipative particle dynamics (DPD) are used in order to describe individual cells and relatively small cell populations, partial differential equations (PDE) are used to describe concentrations of bio-chemical substances in the extracellular matrix, and ordinary differential equations for intracellular regulatory networks (Figure ). An emphasis would be made on developing codes that are both flexible and powerful enough to implement variants of the model, perform simulations, produce desired outputs, and provide tools for analysis; to do so:

  7. We planned to contribute to a recent project named chronos, whose code (written in C++) represents heterogeneous populations of individual cells evolving in time and interacting physically and biochemically, and the objective is to make the code flexible enough to implement different formalisms within the same model, so that different components of the model can be represented in the most appropriate way;

  8. Partial differential equations (PDE) analysis, with a focus on reaction-diffusion equations, transport equations (hyperbolic PDEs) in which the structure can be age, maturity, protein concentration, etc., with particular cases where transport equations are reduced to delay differential equations (DDE).

Comments on these objectives over the evaluation period

We have had strong contributions to objectives , , , , and consequently to objective , as well as to objective , as mentioned in previous sections. These contributions represented the core of the team’s research activity over the evaluation period, as stressed by our publications. It is however noticeable that multi-scale modeling of the immune response and of pathological hematopoiesis (leukemia) has come to represent a proportionally more important part of our activity.

Objective has been cancelled few months after the previous evaluation, following meetings with clinicians who did not show any particular interest in our approaches. The modeling of chronic myeloid leukemia instead took a bigger part of the team’s research activity, both project being at the time coordinated by Thomas Lepoutre.

Objective has been pursued, the project chronos evolved to a better defined project SiMuScale that is currently being developed and aims at structuring the team’s activity and providing a simulation platform that could be adapted to various biological questions necessitating multi-scale modeling.

Objectives for the next four years

The main objectives for the next four years are to continue to improve the 3 previous points: 1) Mathematical and computational modeling for cell population dynamics; 2) Multi-scale modeling of hematopoiesis and leukemia; 3) Multi-scale modeling of the immune response. In addition, we will pursue our effort to develop a simulation platform for multi-scale models (SiMuScale) and we intend to develop the use of mixed effect models and other statistical approaches to deal with the challenges offered by modern biology, in particular the generation of single cell data.

Last activity report : 2019

Les résultats

New Results

Mathematical models describing the interaction between cancer and immune cells in the lymph node

To study the interplay between tumor progression and the immune response, we develop in [article] two new models describing the interaction between cancer and immune cells in the lymph node. The first model consists of partial differential equations (PDEs) describing the populations of the different types of cells. The second one is a hybrid discrete-continuous model integrating the mechanical and biochemical mechanisms that define the tumor-immune interplay in the lymph node. We use the continuous model to determine the conditions of the regimes of tumor-immune interaction in the lymph node. While we use the hybrid model to elucidate the mechanisms that contribute to the development of each regime at the cellular and tissue levels. We study the dynamics of tumor growth in the absence of immune cells. Then, we consider the immune response and we quantify the effects of immunosuppression and local EGF concentration on the fate of the tumor.

WASABI: a dynamic iterative framework for gene regulatory network inference

Background Inference of gene regulatory networks from gene expression data has been a long-standing and notoriously difficult task in systems biology. Recently, single-cell transcriptomic data have been massively used for gene regulatory network inference, with both successes and limitations. In the work [article], we propose an iterative algorithm called WASABI, dedicated to inferring a causal dynamical network from time-stamped single-cell data, which tackles some of the limitations associated with current approaches. We first introduce the concept of waves, which posits that the information provided by an external stimulus will affect genes one-by-one through a cascade, like waves spreading through a network. This concept allows us to infer the network one gene at a time, after genes have been ordered regarding their time of regulation. We then demonstrate the ability of WASABI to correctly infer small networks, which have been simulated in silico using a mechanistic model consisting of coupled piecewise-deterministic Markov processes for the proper description of gene expression at the single-cell level. We finally apply WASABI on in vitro generated data on an avian model of erythroid differentiation. The structure of the resulting gene regulatory network sheds a new light on the molecular mechanisms controlling this process. In particular, we find no evidence for hub genes and a much more distributed network structure than expected. Interestingly, we find that a majority of genes are under the direct control of the differentiation-inducing stimulus. Conclusions Together, these results demonstrate WASABI versatility and ability to tackle some general gene regulatory networks inference issues. It is our hope that WASABI will prove useful in helping biologists to fully exploit the power of time-stamped single-cell data.

A multiscale model of platelet-fibrin thrombus growth in the flow

Thrombosis is a life-threatening clinical condition characterized by the obstruction of blood flow in a vessel due to the formation of a large thrombus. The pathogenesis of thrombosis is complex because the type of formed clots depends on the location and function of the corresponding blood vessel. To explore this phenomenon, we develop in [article] a novel multiscale model of platelet-fibrin thrombus growth in the flow. In this model, the regulatory network of the coagulation cascade is described by partial differential equations. Blood flow is introduced using the Navier–Stokes equations and the clot is treated as a porous medium. Platelets are represented as discrete spheres that migrate with the flow. Each platelet can attach to the thrombus, aggregate, become activated, express proteins on its surface, detach, and/or become non-adhesive. The interaction of platelets with blood flow is captured using the Immersed Boundary Method (IBM). We use the model to investigate the role of flow conditions in shaping the dynamics of venous and arterial thrombi. We describe the formation of red and white thrombi under venous and arterial flow respectively and highlight the main characteristics of each type. We identify the different regimes of normal and pathological thrombus formation depending on flow conditions.

Mathematical modeling of platelet production

  • In [article], we analyze the existence of oscillating solutions and the asymptotic convergence for a nonlinear delay differential equation arising from the modeling of platelet production. We consider four different cell compartments corresponding to different cell maturity levels: stem cells, megakaryocytic progenitors, megakaryocytes, and platelets compartments, and the quantity of circulating thrombopoietin (TPO), a platelet regulation cytokine.

  • In [article], we analyze the stability of a differential equation with two delays originating from a model for a population divided into two subpopulations, immature and mature, and we apply this analysis to a model for platelet production. The dynamics of mature individuals is described by the following nonlinear differential equation with two delays: x(t)=λx(t)+g(x(tτ1))g(x(tτ1τ2))eλτ2 . The method of D -decomposition is used to compute the stability regions for a given equilibrium. The center manifold theory is used to investigate the steady-state bifurcation and the Hopf bifurcation. Similarly, analysis of the center manifold associated with a double bifurcation is used to identify a set of parameters such that the solution is a torus in the pseudo- phase space. Finally, the results of the local stability analysis are used to study the impact of an increase of the death rate γ or of a decrease of the survival time τ2 of platelets on the onset of oscillations. We show that the stability is lost through a small decrease of survival time (from 8.4 to 7 days), or through an important increase of the death rate (from 0.05 to 0.625 days 1 ).

  • In [article], we analyze the stability of a system of differential equations with a threshold-defined delay arising from a model for platelet production. We consider a maturity-structured population of megakaryocyte progenitors and an age-structured population of platelets, where the cytokine thrombopoietin (TPO) increases the maturation rate of progenitors. Using the quasi-steady-state approximation for TPO dynamics and the method of characteristics, partial differential equations are reduced to a system of two differential equations with a state-dependent delay accounting for the variable maturation rate. We start by introducing the model and proving the positivity and boundedness of the solutions. Then we use a change of variables to obtain an equivalent system of two differential equations with a constant delay, from which we prove existence and uniqueness of the solution. As linearization around the unique positive steady state yields a transcendental characteristic equation of third degree, we introduce the main result, a new framework for stability analysis on models with fixed delays. This framework is then used to describe the stability of the megakaryopoiesis with respect to its parameters. Finally, with parameters being obtained and estimated from data, we give an example in which oscillations appear when the death rate of progenitors is increased 10-fold.

Nonlinear analysis of a model for yeast cell communication

In [article], we study the non-linear stability of a coupled system of two non-linear transport-diffusion equations set in two opposite half-lines. This system describes some aspects of yeast pairwise cellular communication, through the concentration of some protein in the cell bulk and at the cell boundary. We show that it is of bistable type, provided that the intensity of active molecular transport is large enough. We prove the non-linear stability of the most concentrated steady state, for large initial data, by entropy and comparison techniques. For small initial data we prove the self-similar decay of the molecular concentration towards zero. Informally speaking, the rise of a dialog between yeast cells requires enough active molecular transport in this model. Besides, if the cells do not invest enough in the communication with their partner, they do not respond to each other; but a sufficient initial input from each cell in the dialog leads to the establishment of a stable activated state in both cells.

Alzheimer’s disease and prion: An in vitro mathematical model

Alzheimer’s disease ( AD ) is a fatal incurable disease leading to progressive neuron destruction. AD is caused in part by the accumulation in the brain of Aβ monomers aggregating into oligomers and fibrils. Oligomers are amongst the most toxic structures as they can interact with neurons via membrane receptors, including PrPc proteins. This interaction leads to the misconformation of PrPc into pathogenic oligomeric prions, PrPol . In [article], we develop a model describing in vitro Aβ polymerization process. We include interactions between oligomers and PrPc , causing the misconformation of PrPc into PrPol . The model consists of nine equations, including size structured transport equations, ordinary differential equations and delayed differential equations. We analyse the well-posedness of the model and prove the existence and uniqueness of solutions of our model using Schauder fixed point theorem and Cauchy-Lipschitz theorem. Numerical simulations are also provided to give an illustration of the profiles that can be obtained with this model.

Calibration, Selection and Identifiability Analysis of a Mathematical Model of the in vitro Erythropoiesis in Normal and Perturbed Contexts

The in vivo erythropoiesis, which is the generation of mature red blood cells in the bone marrow of whole organisms, has been described by a variety of mathematical models in the past decades. However, the in vitro erythropoiesis, which produces red blood cells in cultures, has received much less attention from the modelling community. In the paper [article], we propose the first mathematical model of in vitro erythropoiesis. We start by formulating different models and select the best one at fitting experimental data of in vitro erythropoietic differentiation obtained from chicken erythroid progenitor cells. It is based on a set of linear ODE, describing 3 hypothetical populations of cells at different stages of differentiation. We then compute confidence intervals for all of its parameters estimates, and conclude that our model is fully identifiable. Finally, we use this model to compute the effect of a chemical drug called Rapamycin, which affects all states of differentiation in the culture, and relate these effects to specific parameter variations. We provide the first model for the kinetics of in vitro cellular differentiation which is proven to be identifiable. It will serve as a basis for a model which will better account for the variability which is inherent to the experimental protocol used for the model calibration.

Model-based assessment of the role of uneven partitioning of molecular content on heterogeneity and regulation of differentiation in CD8 T-cell immune responses

Activation of naive CD8 T-cells can lead to the generation of multiple effector and memory subsets. Multiple parameters associated with activation conditions are involved in generating this diversity that is associated with heterogeneous molecular contents of activated cells. Although naive cell polarisation upon antigenic stimulation and the resulting asymmetric division are known to be a major source of heterogeneity and cell fate regulation, the consequences of stochastic uneven partitioning of molecular content upon subsequent divisions remain unclear yet. In [article], we aim at studying the impact of uneven partitioning on molecular-content heterogeneity and then on the immune response dynamics at the cellular level. To do so, we introduce a multiscale mathematical model of the CD8 T-cell immune response in the lymph node. In the model, cells are described as agents evolving and interacting in a 2D environment while a set of differential equations, embedded in each cell, models the regulation of intra and extracellular proteins involved in cell differentiation. Based on the analysis of in silico data at the single cell level, we 1 show that immune response dynamics can be explained by the molecular-content heterogeneity generated by uneven partitioning at cell division. In particular, uneven partitioning acts as a regulator of cell differentiation and induces the emergence of two coexisting sub-populations of cells exhibiting antagonistic fates. We show that the degree of unevenness of molecular partitioning, along all cell divisions, affects the outcome of the immune response and can promote the generation of memory cells.

Spatial lymphocyte dynamics in lymph nodes predicts the cytotoxic T-Cell frequency needed for HIV infection control

The surveillance of host body tissues by immune cells is central for mediating their defense function. In vivo imaging technologies have been used to quantitatively characterize target cell scanning and migration of lymphocytes within lymph nodes (LNs). The translation of these quantitative insights into a predictive understanding of immune system functioning in response to various perturbations critically depends on computational tools linking the individual immune cell properties with the emergent behavior of the immune system. By choosing the Newtonian second law for the governing equations, we developed in [article] a broadly applicable mathematical model linking individual and coordinated T-cell behaviors. The spatial cell dynamics is described by a superposition of autonomous locomotion, intercellular interaction, and viscous damping processes. The model is calibrated using in vivo data on T-cell motility metrics in LNs such as the translational speeds, turning angle speeds, andmeandering indices. Themodel is applied to predict the impact of T-cell motility on protection against HIV infection, i.e., to estimate the threshold frequency of HIV-specific cytotoxic T cells (CTLs) that is required to detect productively infected cells before the release of viral particles starts. With this, it provides guidance for HIV vaccine studies allowing for the migration of cells in fibrotic LNs.

Drugs modulating stochastic gene expression affect the erythroid differentiation process

To better understand the mechanisms behind cells decision-making to differentiate, we assessed in [article] the influence of stochastic gene expression (SGE) modulation on the erythroid dif- ferentiation process. It has been suggested that stochastic gene expression has a role in cell fate decision-making which is revealed by single-cell analyses but studies dedicated to demonstrate the consistency of this link are still lacking. Recent observations showed that SGE significantly increased during differentiation and a few showed that an increase of the level of SGE is accompanied by an increase in the differentiation process. However, a con- sistent relation in both increasing and decreasing directions has never been shown in the same cellular system. Such demonstration would require to be able to experimentally manipulate simultaneously the level of SGE and cell differentiation in order to observe if cell behavior matches with the current theory. We identified three drugs that modulate SGE in primary erythroid progenitor cells. Both Artemisinin and Indomethacin decreased SGE and reduced the amount of differentiated cells. On the contrary, a third component called MB-3 simultaneously increased the level of SGE and the amount of differentiated cells. We then used a dynamical modelling approach which confirmed that differentiation rates were indeed affected by the drug treatment. Using single-cell analysis and modeling tools, we provide experimental evidence that, in a physiologically relevant cellular system, SGE is linked to differentiation.

Stochastic gene expression with a multistate promoter: breaking down exact distributions

We consider in [article] a stochastic model of gene expression in which transcription depends on a multistate promoter, including the famous two-state model and refractory promoters as special cases, and focus on deriving the exact stationary distribution. Building upon several successful approaches, we present a more unified viewpoint that enables us to simplify and generalize existing results. In particular, the original jump process is deeply related to a multivariate piecewise-deterministic Markov process that may also be of interest beyond the biological field. In a very particular case of promoter configuration, this underlying process is shown to have a simple Dirichlet stationary distribution. In the general case, the corresponding marginal distributions extend the well-known class of Beta products, involving complex parameters that directly relate to spectral properties of the promoter transition matrix. Finally, we illustrate these results with biologically plausible examples.

Cell generation dynamics underlying naive T-cell homeostasis in adult humans

Thymic involution and proliferation of naive T-cells both contribute to shaping the naive T-cell repertoire as humans age, but a clear understanding of the roles of each throughout a human life span has been difficult to determine. By measuring nuclear bomb test-derived 14 C in genomic DNA, we determined in [article] the turnover rates of CD4 + and CD8 + naive T-cell populations and defined their dynamics in healthy individuals ranging from 20 to 65 years of age. We demonstrate that naive T-cell generation decreases with age because of a combination of declining peripheral division and thymic production during adulthood. Concomitant decline in T-cell loss compensates for decreased generation rates. We investigated putative mechanisms underlying age-related changes in homeostatic regulation of CD4+ naive T-cell turnover, using mass cytometry to profile candidate signaling pathways involved in T-cell activation and proliferation relative to CD31 expression, a marker of thymic proximity for the CD4+ naive T-cell population. We show that basal nuclear factor kB (NF-kB) phosphory-lation positively correlated with CD31 expression and thus is decreased in peripherally expanded naive T-cell clones. Functionally, we found that NF-kB signaling was essential for naive T-cell proliferation to the homeostatic growth factor interleukin (IL)-7, and reduced NF-kB phosphorylation in CD4 + CD31 − naive T cells is linked to reduced homeostatic proliferation potential. Our results reveal an age-related decline in naive T-cell turnover as a putative regulator of naive T-cell diversity and identify a molecular pathway that restricts proliferation of peripherally expanded naive T-cell clones that accumulate with age.

Erythroid differentiation displays a peak of energy consumption concomitant with glycolytic metabolism rearrangements

Our previous single-cell based gene expression analysis pointed out significant variations of LDHA level during erythroid differentiation. Deeper investigations highlighted that a metabolic switch occurred along differentiation of erythroid cells. More precisely we showed in [article] that self-renewing progenitors relied mostly upon lactate-productive glycolysis, and required LDHA activity, whereas differentiating cells, mainly involved mitochondrial oxidative phosphorylation (OXPHOS). These metabolic rearrangements were coming along with a particular temporary event, occurring within the first 24h of erythroid differentiation. The activity of glycolytic metabolism and OXPHOS rose jointly with oxgene consumption dedicated to ATP production at 12-24h of the differentiation process before lactate-productive glycolysis sharply fall down and energy needs decline. Finally, we demonstrated that the metabolic switch mediated through LDHA drop and OXPHOS upkeep might be necessary for erythroid differentiation. We also discuss the possibility that metabolism, gene expression and epigenetics could act together in a circular manner as a driving force for differentiation.

Les commentaires sont clos.