Nathalie Eymard

The role of spatial organization of cells in erythropoiesis

Erythropoiesis, the process of red blood cell production, occurs mainly in the bone marrow. The functional unit of mammalian erythropoiesis, the erythroblastic island, consists of a central macrophage surrounded by adherent erythroid progenitor cells (CFU-E/Pro-EBs) and their differentiating progeny, the erythroblasts. Central macrophages display on their surface or secrete various growth or inhibitory factors that influence the fate of the surrounding erythroid cells. CFU-E/Pro-EBs have three possible fates: (a) expansion of their numbers without differentiation, (b) differentiation into reticulocytes that are released into the blood, (c) death by apoptosis. CFU-E/Pro-EB fate is under the control of a complex molecular network, that is highly dependent upon environmental conditions in the erythroblastic island. In order to assess the functional role of space coupled with the complex network behavior in erythroblastic islands, we developed hybrid discrete-continuous models of erythropoiesis. A model was developed in which cells are considered as individual physical objects, intracellular regulatory networks are modeled with ordinary differential equations and extracellular concentrations by partial differential equations. We used the model to investigate the impact of an important difference between humans and mice in which mature late-stage erythroblasts produce the most Fas-ligand in humans, whereas early-stage erythroblasts produce the most Fas-ligand in mice. Although the global behaviors of the erythroblastic islands in both species were similar, differences were found, including a relatively slower response time to acute anemia in humans. Also, our modeling approach was very consistent with in vitro culture data, where the central macrophage in reconstituted erythroblastic islands has a strong impact on the dynamics of red blood cell production. The specific spatial organization of erythroblastic islands is key to the normal, stable functioning of mammalian erythropoiesis, both in vitro and in vivo. Our model of a simplified molecular network controlling cell decision provides a realistic functional unit of mammalian erythropoiesis that integrates multiple microenvironmental influences within the erythroblastic island with those of circulating regulators of erythropoiesis, such as EPO and glucocorticosteroids, that are produced at remote sites.

Hybrid Model of Erythropoiesis

A hybrid model of cell dynamics is presented. It is illustrated by model examples and applied to study erythropoiesis (red blood cell production). In this approach, cells are considered as discrete objects while intra-cellular proteins and extra-cellular biochemical substances are described with continuous models. Spatial organization of erythropoiesis occurring in specific structures of the bone marrow, called erythroblastic island, is investigated.

A methodological framework for drug development in rare diseases.

IntroductionDeveloping orphan drugs is challenging because of their severity and the requisite for effective drugs. The small number of patients does not allow conducting adequately powered randomized controlled trials (RCTs). There is a need to develop high quality, ethically investigated, and appropriately authorized medicines, without subjecting patients to unnecessary trials.Aims and ObjectivesThe main aim is to develop generalizable framework for choosing the best-performing drug/endpoint/design combinations in orphan drug development using an in silico modeling and trial simulation approach. The two main objectives were (i) to provide a global strategy for each disease to identify the most relevant drugs to be evaluated in specific patients during phase III RCTs, (ii) and select the best design for each drug to be used in future RCTs.Methodological approachIn silico phase III RCT simulation will be used to find the optimal trial design and was carried out in two steps: (i) statistical analysis of available clinical databases and (ii) integrative modeling that combines mathematical models for diseases with pharmacokinetic-pharmacodynamics models for the selected drug candidates.ConclusionThere is a need to speed up the process of orphan drug development, develop new methods for translational research and personalized medicine, and contribute to European Medicines Agency guidelines. The approach presented here offers many perspectives in clinical trial conception.

An in silico approach helped to identify the best experimental design, population, and outcome for future randomized clinical trials.

OBJECTIVES The main objective of our work was to compare different randomized clinical trial (RCT) experimental designs in terms of power, accuracy of the estimation of treatment effect, and number of patients receiving active treatment using in silico simulations. STUDY DESIGN AND SETTING A virtual population of patients was simulated and randomized in potential clinical trials. Treatment effect was modeled using a dose-effect relation for quantitative or qualitative outcomes. Different experimental designs were considered, and performances between designs were compared. One thousand clinical trials were simulated for each design based on an example of modeled disease. RESULTS According to simulation results, the number of patients needed to reach 80% power was 50 for crossover, 60 for parallel or randomized withdrawal, 65 for drop the loser (DL), and 70 for early escape or play the winner (PW). For a given sample size, each design had its own advantage: low duration (parallel, early escape), high statistical power and precision (crossover), and higher number of patients receiving the active treatment (PW and DL). CONCLUSION Our approach can help to identify the best experimental design, population, and outcome for future RCTs. This may be particularly useful for drug development in rare diseases, theragnostic approaches, or personalized medicine.

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.

Mathematical model of T-cell lymphoblastic lymphoma: disease, treatment, cure or relapse of a virtual cohort of patients

T lymphoblastic lymphoma (T-LBL) is a rare type of lymphoma with a good prognosis with a remission rate of 85%. Patients can be completely cured or can relapse during or after a 2-year treatment. Relapses usually occur early after the remission of the acute phase. The median time of relapse is equal to 1 year, after the occurrence of complete remission (range 0.2-5.9 years) (Uyttebroeck et al., 2008). It can be assumed that patients may be treated longer than necessary with undue toxicity.The aim of our model was to investigate whether the duration of the maintenance therapy could be reduced without increasing the risk of relapses and to determine the minimum treatment duration that could be tested in a future clinical trial.We developed a mathematical model of virtual patients with T-LBL in order to obtain a proportion of virtual relapses close to the one observed in the real population of patients from the EuroLB database. Our simulations reproduced a 2-year follow-up required to study the onset of the disease, the treatment of the acute phase and the maintenance treatment phase.