Anass Bouchnita

Hybrid approach to model the spatial regulation of T cell responses

BackgroundMoving from the molecular and cellular level to a multi-scale systems understanding of immune responses requires the development of novel approaches to integrate knowledge and data from different biological levels into mechanism-based integrative mathematical models. The aim of our study is to present a methodology for a hybrid modelling of immunological processes in their spatial context.MethodsA two-level hybrid mathematical model of immune cell migration and interaction integrating cellular and organ levels of regulation for a 2D spatial consideration of idealized secondary lymphoid organs is developed. It considers the population dynamics of antigen-presenting cells, CD4 + and CD8 + T lymphocytes in naive-, proliferation- and differentiated states. Cell division is assumed to be asymmetric and regulated by the extracellular concentration of interleukin-2 (IL-2) and type I interferon (IFN), together controlling the balance between proliferation and differentiation. The cytokine dynamics is described by reaction-diffusion PDEs whereas the intracellular regulation is modelled with a system of ODEs.ResultsThe mathematical model has been developed, calibrated and numerically implemented to study various scenarios in the regulation of T cell immune responses to infection, in particular the change in the diffusion coefficient of type I IFN as compared to IL-2. We have shown that a hybrid modelling approach provides an efficient tool to describe and analyze the interplay between spatio-temporal processes in the emergence of abnormal immune response dynamics.DiscussionVirus persistence in humans is often associated with an exhaustion of T lymphocytes. Many factors can contribute to the development of exhaustion. One of them is associated with a shift from a normal clonal expansion pathway to an altered one characterized by an early terminal differentiation of T cells. We propose that an altered T cell differentiation and proliferation sequence can naturally result from a spatial separation of the signaling events delivered via TCR, IL-2 and type I IFN receptors. Indeed, the spatial overlap of the concentration fields of extracellular IL-2 and IFN in lymph nodes changes dynamically due to different migration patterns of APCs and CD4 + T cells secreting them.ConclusionsThe proposed hybrid mathematical model of the immune response represents a novel analytical tool to examine challenging issues in the spatio-temporal regulation of cell growth and differentiation, in particular the effect of timing and location of activation signals.

On the regimes of blood coagulation

Abstract Reaction–diffusion system of equations describing blood clotting is studied. Different regimes of clot growth are identified in a quiescent plasma and in blood flow depending on the relative strength of initiation, propagation and inhibition of the thrombin production.

Conditions of microvessel occlusion for blood coagulation in flow.

Vessel occlusion is a perturbation of blood flow inside a blood vessel because of the fibrin clot formation. As a result, blood circulation in the vessel can be slowed down or even stopped. This can provoke the risk of cardiovascular events. In order to explore this phenomenon, we used a previously developed mathematical model of blood clotting to describe the concentrations of blood factors with a reaction-diffusion system of equations. The Navier-Stokes equations were used to model blood flow, and we treated the clot as a porous medium. We identify the conditions of partial or complete occlusion in a small vessel depending on various physical and physiological parameters. In particular, we were interested in the conditions on blood flow and diameter of the wounded area. The existence of a critical flow velocity separating the regimes of partial and complete occlusion was demonstrated through the mathematical investigation of a simplified model of thrombin wave propagation in Poiseuille flow. We observed different regimes of vessel occlusion depending on the model parameters both for the numerical simulations and in the theoretical study. Then, we compared the rate of clot growth in flow obtained in the simulations with experimental data. Both of them showed the existence of different regimes of clot growth depending on the velocity of blood flow.

Reaction-diffusion waves of blood coagulation ☆

One of the main characteristics of blood coagulation is the speed of clot growth. In the current work we consider a mathematical model of the coagulation cascade and study existence, stability and speed of propagation of the reaction-diffusion waves of blood coagulation. We also develop a simplified one-equation model that reflects the main features of the thrombin wave propagation. For this equation we estimate the wave speed analytically. The resulting formulas provide a good approximation for the speed of wave propagation in a more complex model as well as for the experimental data.

Bone marrow infiltration by multiple myeloma causes anemia by reversible disruption of erythropoiesis

Multiple myeloma (MM) infiltrates bone marrow and causes anemia by disrupting erythropoiesis, but the effects of marrow infiltration on anemia are difficult to quantify. Marrow biopsies of newly diagnosed MM patients were analyzed before and after four 28‐day cycles of nonerythrotoxic remission induction chemotherapy. Complete blood cell counts and serum paraprotein concentrations were measured at diagnosis and before each chemotherapy cycle. At diagnosis, marrow area infiltrated by myeloma correlated negatively with hemoglobin, erythrocytes, and marrow erythroid cells. After successful chemotherapy, patients with less than 30% myeloma infiltration at diagnosis had no change in these parameters, whereas patients with more than 30% myeloma infiltration at diagnosis increased all three parameters. Clinical data were used to develop mathematical models of the effects of myeloma infiltration on the marrow niches of terminal erythropoiesis, the erythroblastic islands (EBIs). A hybrid discrete‐continuous model of erythropoiesis based on EBI structure/function was extended to sections of marrow containing multiple EBIs. In the model, myeloma cells can kill erythroid cells by physically destroying EBIs and by producing proapoptotic cytokines. Following chemotherapy, changes in serum paraproteins as measures of myeloma cells and changes in erythrocyte numbers as measures of marrow erythroid cells allowed modeling of myeloma cell death and erythroid cell recovery, respectively. Simulations of marrow infiltration by myeloma and treatment with nonerythrotoxic chemotherapy demonstrate that myeloma‐mediated destruction and subsequent reestablishment of EBIs and expansion of erythroid cell populations in EBIs following chemotherapy provide explanations for anemia development and its therapy‐mediated recovery in MM patients. Am. J. Hematol. 91:371–378, 2016.

A Hybrid Computation Model to Describe the Progression of Multiple Myeloma and Its Intra-Clonal Heterogeneity

Multiple myeloma (MM) is a genetically complex hematological cancer that is characterized by proliferation of malignant plasma cells in the bone marrow. MM evolves from the clonal premalignant disorder monoclonal gammopathy of unknown significance (MGUS) by sequential genetic changes involving many different genes, resulting in dysregulated growth of multiple clones of plasma cells. The migration, survival, and proliferation of these clones require the direct and indirect interactions with the non-hematopoietic cells of the bone marrow. We develop a hybrid discrete-continuous model of MM development from the MGUS stage. The discrete aspect of the modelisobservedatthecellularlevel: cellsarerepresentedasindividualobjectswhichmove,interact, divide, and die by apoptosis. Each of these actions is regulated by intracellular and extracellular processes as described by continuous models. The hybrid model consists of the following submodels that have been simplified from the much more complex state of evolving MM: cell motion due to chemotaxis, intracellular regulation of plasma cells, extracellular regulation in the bone marrow, and acquisition of mutations upon cell division. By extending a previous, simpler model in which the extracellular matrix was considered to be uniformly distributed, the new hybrid model provides a more accurate description in which cytokines are produced by the marrow microenvironment and consumed by the myeloma cells. The complex multiple genetic changes in MM cells and the numerous cell-cell and cytokine-mediated interactions between myeloma cells and their marrow microenviroment are simplified in the model such that four related but evolving MM clones can be studied as they compete for dominance in the setting of intraclonal heterogeneity.