1
|
Eymard N, Bessonov N, Volpert V, Kurbatova P, Gueyffier F, Nony P. Pharmacokinetic/pharmacodynamic model of a methionine starvation based anti-cancer drug. Med Biol Eng Comput 2023:10.1007/s11517-023-02786-2. [PMID: 36882575 DOI: 10.1007/s11517-023-02786-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/28/2020] [Accepted: 01/17/2023] [Indexed: 03/09/2023]
Abstract
A new therapeutic approach against cancer is developed by the firm Erytech. This approach is based on starved cancer cells of an amino acid essential to their growth (the L-methionine). The depletion of plasma methionine level can be induced by an enzyme, the methionine-γ-lyase. The new therapeutic formulation is a suspension of erythrocytes encapsulating the activated enzyme. Our work reproduces a preclinical trial of a new anti-cancer drug with a mathematical model and numerical simulations in order to replace animal experiments and to have a deeper insight on the underlying processes. With a combination of a pharmacokinetic/pharmacodynamic model for the enzyme, substrate, and co-factor with a hybrid model for tumor, we develop a "global model" that can be calibrated to simulate different human cancer cell lines. The hybrid model includes a system of ordinary differential equations for the intracellular concentrations, partial differential equations for the concentrations of nutrients and drugs in the extracellular matrix, and individual based model for cancer cells. This model describes cell motion, division, differentiation, and death determined by the intracellular concentrations. The models are developed on the basis of experiments in mice carried out by Erytech. Parameters of the pharmacokinetics model were determined by fitting a part of experimental data on the concentration of methionine in blood. Remaining experimental protocols effectuated by Erytech were used to validate the model. The validated PK model allowed the investigation of pharmacodynamics of cell populations. Numerical simulations with the global model show cell synchronization and proliferation arrest due to treatment similar to the available experiments. Thus, computer modeling confirms a possible effect of treatment based on the decrease of methionine concentration. The main goal of the study is the development of an integrated pharmacokinetic/pharmacodynamic model for encapsulated methioninase and of a mathematical model of tumor growth/regression in order to determine the kinetics of L-methionine depletion after co-administration of Erymet product and Pyridoxine.
Collapse
Affiliation(s)
- N Eymard
- Institut Camille Jordan, UMR 5208 CNRS, University Lyon 1, 69622, Villeurbanne, France.
| | - N Bessonov
- Institute of Mechanical Engineering Problems, 199178, Saint Petersburg, Russia
| | - V Volpert
- Institut Camille Jordan, UMR 5208 CNRS, University Lyon 1, 69622, Villeurbanne, France.,Peoples Friendship University of Russia (RUDN University), 6 Miklukho-Maklaya St, Moscow, 117198, Russia
| | - P Kurbatova
- Institut Camille Jordan, UMR 5208 CNRS, University Lyon 1, 69622, Villeurbanne, France
| | - F Gueyffier
- CNRS, Laboratoire de Biométrie et Biologie Evolutive UMR 5558, Université Lyon 1, F-69622, Villeurbanne, France
| | - P Nony
- Institut Camille Jordan, UMR 5208 CNRS, University Lyon 1, 69622, Villeurbanne, France.,CNRS, Laboratoire de Biométrie et Biologie Evolutive UMR 5558, Université Lyon 1, F-69622, Villeurbanne, France
| |
Collapse
|
2
|
Etheve L, Courcelles E, Lefaudeux D, Kurbatova P, Gomez G, Kahoul R, Boissel JP, Kulesza A, Dwivedi G, Schmidely N, Bertocchio JP, Gittoes N, Siebert U, Bechet E. 4Essais cliniques in silico, une approche innovante visant à compléter les essais cliniques dans le domaine des maladies rares : validation d'un modèle computationnel chez les patients atteints d'hypoparathyroïdie. Rev Epidemiol Sante Publique 2022. [DOI: 10.1016/j.respe.2022.09.063] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022] Open
|
3
|
Eymard N, Volpert V, Kurbatova P, Volpert V, Bessonov N, Ogungbenro K, Aarons L, Janiaud P, Nony P, Bajard A, Chabaud S, Bertrand Y, Kassaï B, Cornu C, Nony P. Mathematical model of T-cell lymphoblastic lymphoma: disease, treatment, cure or relapse of a virtual cohort of patients. Math Med Biol 2018; 35:25-47. [PMID: 28082512 DOI: 10.1093/imammb/dqw019] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/27/2015] [Accepted: 10/12/2016] [Indexed: 12/19/2022]
Abstract
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.
Collapse
Affiliation(s)
- N Eymard
- Institut Camille Jordan, UMR, CNRS, University Lyon 1, Villeurbanne, France
| | - V Volpert
- Institut Camille Jordan, UMR, CNRS, University Lyon 1, Villeurbanne, France
| | - P Kurbatova
- Institut Camille Jordan, UMR, CNRS, University Lyon 1, Villeurbanne, France
| | - V Volpert
- INRIA Team Dracula, INRIA Antenne Lyon la Doua 69603 Villeurbanne, France
| | - N Bessonov
- Institute of Mechanical Engineering Problems, Saint Petersburg, Russia
| | - K Ogungbenro
- Centre for Applied Pharmacokinetic Research, Manchester Pharmacy School The University of Manchester, Manchester, UK
| | - L Aarons
- Centre for Applied Pharmacokinetic Research, Manchester Pharmacy School The University of Manchester, Manchester, UK
| | - P Janiaud
- University of Lyon 1, UMR, CNRS, Lyon, France
| | - P Nony
- University of Lyon 1, UMR, CNRS, Lyon, France
| | - A Bajard
- Unité de Biostatistique et d'Evaluation des Thérapeutiques Centre Léon Bérard, Lyon, France
| | - S Chabaud
- Unité de Biostatistique et d'Evaluation des Thérapeutiques Centre Léon Bérard, Lyon, France
| | - Y Bertrand
- Institute of Hematology and Oncology Paediatrics, Hospices Civils de Lyon, University Claude Bernard Lyon I, Lyon, France
| | - B Kassaï
- Hospices Civils de Lyon, Centre d'Investigation Clinique, INSERM CIC1407, Lyon, France
| | - C Cornu
- Hospices Civils de Lyon, Centre d'Investigation Clinique, INSERM CIC1407, Lyon, France
| | - P Nony
- CHU Lyon, Service de Pharmacologie Clinique et Essais Thérapeutiques, Lyon, France
| | | |
Collapse
|
4
|
Bouchnita A, Galochkina T, Kurbatova P, Nony P, Volpert V. Conditions of microvessel occlusion for blood coagulation in flow. Int J Numer Method Biomed Eng 2017; 33:e2850. [PMID: 27863131 DOI: 10.1002/cnm.2850] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/13/2016] [Revised: 09/26/2016] [Accepted: 10/25/2016] [Indexed: 06/06/2023]
Abstract
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.
Collapse
Affiliation(s)
- A Bouchnita
- Institut Camille Jordan, UMR 5208 CNRS, University Lyon 1, Villeurbanne, 69622, France
- Laboratoire de Biométrie et Biologie Evolutive, UMR 5558 CNRS, University Lyon 1, Lyon, 69376, France
- Laboratory of Study and Research in Applied Mathematics, Mohammadia School of Engineers, Mohamed V University, Rabat, Morocco
| | - T Galochkina
- Institut Camille Jordan, UMR 5208 CNRS, University Lyon 1, Villeurbanne, 69622, France
- Department of Biophysics, Faculty of Biology, M.V. Lomonosov Moscow State University, Leninskie gory 1, Moscow, Russia
- Federal Research Clinical Center of Federal Medical & Biological Agency of Russia, Orekhovy boulevard 28, Moscow, Russia
| | - P Kurbatova
- Institut Camille Jordan, UMR 5208 CNRS, University Lyon 1, Villeurbanne, 69622, France
- Laboratoire de Biométrie et Biologie Evolutive, UMR 5558 CNRS, University Lyon 1, Lyon, 69376, France
| | - P Nony
- Laboratoire de Biométrie et Biologie Evolutive, UMR 5558 CNRS, University Lyon 1, Lyon, 69376, France
- Service de Pharmacologie Clinique, Hospices Civils de Lyon, Lyon, France
| | - V Volpert
- Institut Camille Jordan, UMR 5208 CNRS, University Lyon 1, Villeurbanne, 69622, France
- INRIA Team Dracula, INRIA Lyon La Doua, 69603 Villeurbanne, France
- Laboratoire Poncelet, UMI 2615 CNRS, Bolshoy Vlasyevskiy Pereulok 11, 119002 Moscow, Russia
| |
Collapse
|
5
|
Kurbatova P, Wendling F, Kaminska A, Rosati A, Nabbout R, Guerrini R, Dulac O, Pons G, Cornu C, Nony P, Chiron C, Benquet P. Dynamic changes of depolarizing GABA in a computational model of epileptogenic brain: Insight for Dravet syndrome. Exp Neurol 2016; 283:57-72. [PMID: 27246997 DOI: 10.1016/j.expneurol.2016.05.037] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/18/2015] [Revised: 05/12/2016] [Accepted: 05/26/2016] [Indexed: 10/21/2022]
Abstract
Abnormal reemergence of depolarizing GABAA current during postnatal brain maturation may play a major role in paediatric epilepsies, Dravet syndrome (DS) being among the most severe. To study the impact of depolarizing GABA onto distinct patterns of EEG activity, we extended a neural mass model as follows: one sub-population of pyramidal cells was added as well as two sub-populations of interacting interneurons, perisomatic-projecting interneurons (basket-like) with fast synaptic kinetics GABAA (fast, I1) and dendritic-projecting interneurons with slow synaptic kinetics GABAA (slow, I2). Basket-like cells were interconnected to reproduce mutual inhibition mechanisms (I1➔I1). The firing rate of interneurons was adapted to mimic the genetic alteration of voltage gated sodium channels found in DS patients, SCN1A(+/-). We implemented the "dynamic depolarizing GABAA" mediated post-synaptic potential in the model, as some studies reported that the chloride reversal potential can switch from negative to more positive value depending on interneuron activity. The "shunting inhibition" promoted by GABAA receptor activation was also implemented. We found that increasing the proportion of depolarizing GABAA mediated IPSP (I1➔I1 and I1➔P) only (i.e., other parameters left unchanged) was sufficient to sequentially switch the EEG activity from background to (1) interictal isolated polymorphic epileptic spikes, (2) fast onset activity, (3) seizure like activity and (4) seizure termination. The interictal and ictal EEG patterns observed in 4 DS patients were reproduced by the model via tuning the amount of depolarizing GABAA postsynaptic potential. Finally, we implemented the modes of action of benzodiazepines and stiripentol, two drugs recommended in DS. Both drugs blocked seizure-like activity, partially and dose-dependently when applied separately, completely and with a synergic effect when combined, as has been observed in DS patients. This computational modeling study constitutes an innovative approach to better define the role of depolarizing GABA in infantile onset epilepsy and opens the way for new therapeutic hypotheses, especially in Dravet syndrome.
Collapse
Affiliation(s)
- P Kurbatova
- University Lyon 1, UMR 5558, CRNS, Lyon, France
| | - F Wendling
- UMR 1099, Inserm-University Rennes1, LTSI, Rennes, France
| | - A Kaminska
- UMR 1129, Inserm-Paris Descartes University-CEA, Paris, France
| | - A Rosati
- Pediatric Neurology Unit and Laboratories, Children's Hospital A. Meyer-University of Florence, Firenze, Italy
| | - R Nabbout
- UMR 1129, Inserm-Paris Descartes University-CEA, Paris, France
| | - R Guerrini
- Pediatric Neurology Unit and Laboratories, Children's Hospital A. Meyer-University of Florence, Firenze, Italy; IRCCS Fondazione Stella Maris, Pisa, Italy
| | - O Dulac
- UMR 1129, Inserm-Paris Descartes University-CEA, Paris, France
| | - G Pons
- UMR 1129, Inserm-Paris Descartes University-CEA, Paris, France
| | - C Cornu
- Hôpital Louis Pradel, Centre d'Investigation Clinique, INSERM CIC201/UMR5558, Bron, France; CHU Lyon, Service de Pharmacologie Clinique, Lyon, France; University Lyon 1, UMR 5558, CRNS, Lyon, France
| | - P Nony
- CHU Lyon, Service de Pharmacologie Clinique, Lyon, France; University Lyon 1, UMR 5558, CRNS, Lyon, France
| | | | - C Chiron
- UMR 1129, Inserm-Paris Descartes University-CEA, Paris, France
| | - P Benquet
- UMR 1099, Inserm-University Rennes1, LTSI, Rennes, France
| |
Collapse
|
6
|
Kurbatova P, Bessonov N, Volpert V, Tiddens HAWM, Cornu C, Nony P, Caudri D. Model of mucociliary clearance in cystic fibrosis lungs. J Theor Biol 2015; 372:81-8. [PMID: 25746843 DOI: 10.1016/j.jtbi.2015.02.023] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2014] [Revised: 02/06/2015] [Accepted: 02/25/2015] [Indexed: 11/19/2022]
Abstract
Mucus clearance is a primary innate defense mechanism in the human airways. Cystic fibrosis (CF) is a genetic disease caused by mutations in the gene encoding the cystic fibrosis transmembrane conductance regulator (CFTR) protein. CF is characterized by dehydration of airway surface liquid and impaired mucociliary clearance. As a result, microorganisms are not efficiently removed from the airways, and patients experience chronic pulmonary infections and inflammation. We propose a new physiologically based mathematical model of muco-ciliary transport consisting of the two major components of the mucociliary clearance system: (i) periciliary liquid layer (PCL) and (ii) mucus layer. We study mucus clearance under normal conditions and in CF patients. Restoring impaired clearance of airway secretions in one of the major goals of therapy in patients with CF. We consider the action of the aerosolized and inhaled medication dornase alfa, which reduces the viscosity of cystic fibrosis mucus, by selectively cleaving the long DNA strands it contains. The results of the model simulations stress the potential relevance of the location of the drug deposition in the central or peripheral airways. Mucus clearance was increased in case the drug was primarily deposited peripherally, i.e. in the small airways.
Collapse
Affiliation(s)
- P Kurbatova
- University of Lyon 1, UMR 5558, CRNS Lyon, 8 rue Guillaume Paradin, BP8071, 69376 cedex 08, Lyon, France; UMR 663, Inserm-University Paris Descartes-CEA, Necker Hospital, Paris, France.
| | - N Bessonov
- Institute of Mechanical Engineering Problems, 199178 Saint Petersburg, Russia
| | - V Volpert
- University of Lyon 1, CNRS UMR 5208, Institut Camille Jordan 43 blvd du 11 novembre 1918, F-69622 Villeurbanne-Cedex, France
| | - H A W M Tiddens
- Erasmus University Medical Centre-Sophia Children׳s Hospital, PO Box 2060, 3000 CB Rotterdam, Netherlands
| | - C Cornu
- University of Lyon 1, UMR 5558, CRNS Lyon, 8 rue Guillaume Paradin, BP8071, 69376 cedex 08, Lyon, France; CHU Lyon, Service de Pharmacologie Clinique, 8 rue Guillaume Paradin, BP8071, 69376 cedex 08, Lyon, France; Hôpital Louis Pradel, Centre d׳Investigation Clinique, INSERM CIC201/UMR5558, Bron, France
| | - P Nony
- University of Lyon 1, UMR 5558, CRNS Lyon, 8 rue Guillaume Paradin, BP8071, 69376 cedex 08, Lyon, France; CHU Lyon, Service de Pharmacologie Clinique, 8 rue Guillaume Paradin, BP8071, 69376 cedex 08, Lyon, France
| | - D Caudri
- Erasmus University Medical Centre-Sophia Children׳s Hospital, PO Box 2060, 3000 CB Rotterdam, Netherlands
| |
Collapse
|
7
|
Kurbatova P, Bajard A, Tiddens H, Volpert V, Cornu C, Bessonov N, Kassai B, Chabaud S, Nony P, Caudri D. SFP PC-29 - Modélisation de plans expérimentaux d’essais cliniques randomisés dans la mucoviscidose. (Projet Européen CRESIM). Arch Pediatr 2014. [DOI: 10.1016/s0929-693x(14)72179-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
8
|
Abstract
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.
Collapse
|
9
|
Kurbatova P, Benquet P, Wendling F, Cornu C, Guerrini R, Kaminska A, Dulac O, Pons G, Nabbout R, Chiron C, Nony P. PP270—Computational modeling of dravet syndrome. Clin Ther 2013. [DOI: 10.1016/j.clinthera.2013.07.298] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
10
|
Fischer S, Kurbatova P, Bessonov N, Gandrillon O, Volpert V, Crauste F. Modeling erythroblastic islands: Using a hybrid model to assess the function of central macrophage. J Theor Biol 2012; 298:92-106. [DOI: 10.1016/j.jtbi.2012.01.002] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/14/2011] [Revised: 11/10/2011] [Accepted: 01/03/2012] [Indexed: 11/30/2022]
|