Improving the use of G-CSF during chemotherapy using physiological mathematical modelling : a quantitative systems pharmacology approach

La diminution des doses administrées ou même la cessation complète d'un traitement chimiothérapeutique est souvent la conséquence de la réduction du nombre de neutrophiles, qui sont les globules blancs les plus fréquents dans le sang. Cette réduction dans le nombre absolu des neutrophiles, aussi connue sous le nom de myélosuppression, est précipitée par les effets létaux non spécifiques des médicaments anti-cancéreux, qui, parallèlement à leur effet thérapeutique, produisent aussi des effets toxiques sur les cellules saines. Dans le but d'atténuer cet impact myélosuppresseur, on administre aux patients un facteur de stimulation des colonies de granulocytes recombinant humain (rhG-CSF), une forme exogène du G-CSF, l'hormone responsable de la stimulation de la production des neutrophiles et de leurs libération dans la circulation sanguine. Bien que les bienfaits d'un traitement prophylactique avec le G-CSF pendant la chimiothérapie soient bien établis, les protocoles d'administration demeurent mal définis et sont fréquemment déterminés ad libitum par les cliniciens. Avec l'optique d'améliorer le dosage thérapeutique et rationaliser l'utilisation du rhG-CSF pendant le traitement chimiothérapeutique, nous avons développé un modèle physiologique du processus de granulopoïèse, qui incorpore les connaissances actuelles de pointe relatives à la production des neutrophiles des cellules souches hématopoïétiques dans la moelle osseuse. À ce modèle physiologique, nous avons intégré des modèles pharmacocinétiques/pharmacodynamiques (PK/PD) de deux médicaments: le PM00104 (Zalypsis®), un médicament anti-cancéreux, et le rhG-CSF (filgrastim). En se servant des principes fondamentaux sous-jacents à la physiologie, nous avons estimé les paramètres de manière exhaustive sans devoir recourir à l'ajustement des données, ce qui nous a permis de prédire des données cliniques provenant de 172 patients soumis au protocol CHOP14 (6 cycles de chimiothérapie avec une période de 14 jours où l'administration du rhG-CSF se fait du jour 4 au jour 13 post-chimiothérapie). En utilisant ce modèle physio-PK/PD, nous avons démontré que le nombre d'administrations du rhG-CSF pourrait être réduit de dix (pratique actuelle) à quatre ou même trois administrations, à condition de retarder le début du traitement prophylactique par le rhG-CSF. Dans un souci d'applicabilité clinique de notre approche de modélisation, nous avons investigué l'impact de la variabilité PK présente dans une population de patients, sur les prédictions du modèle, en intégrant des modèles PK de population (Pop-PK) des deux médicaments. En considérant des cohortes de 500 patients in silico pour chacun des cinq scénarios de variabilité plausibles et en utilisant trois marqueurs cliniques, soient le temps au nadir des neutrophiles, la valeur du nadir, ainsi que l'aire sous la courbe concentration-effet, nous avons établi qu'il n'y avait aucune différence significative dans les prédictions du modèle entre le patient-type et la population. Ceci démontre la robustesse de l'approche que nous avons développée et qui s'apparente à une approche de pharmacologie quantitative des systèmes (QSP). Motivés par l'utilisation du rhG-CSF dans le traitement d'autres maladies, comme des pathologies périodiques telles que la neutropénie cyclique, nous avons ensuite soumis l'étude du modèle au contexte des maladies dynamiques. En mettant en évidence la non validité du paradigme de la rétroaction des cytokines pour l'administration exogène des mimétiques du G-CSF, nous avons développé un modèle physiologique PK/PD novateur comprenant les concentrations libres et liées du G-CSF. Ce nouveau modèle PK a aussi nécessité des changements dans le modèle PD puisqu’il nous a permis de retracer les concentrations du G-CSF lié aux neutrophiles. Nous avons démontré que l'hypothèse sous-jacente de l'équilibre entre la concentration libre et liée, selon la loi d'action de masse, n'est plus valide pour le G-CSF aux concentrations endogènes et mènerait en fait à la surestimation de la clairance rénale du médicament. En procédant ainsi, nous avons réussi à reproduire des données cliniques obtenues dans diverses conditions (l'administration exogène du G-CSF, l'administration du PM00104, CHOP14). Nous avons aussi fourni une explication logique des mécanismes responsables de la réponse physiologique aux deux médicaments. Finalement, afin de mettre en exergue l’approche intégrative en pharmacologie adoptée dans cette thèse, nous avons démontré sa valeur inestimable pour la mise en lumière et la reconstruction des systèmes vivants complexes, en faisant le parallèle avec d’autres disciplines scientifiques telles que la paléontologie et la forensique, où une approche semblable a largement fait ses preuves. Nous avons aussi discuté du potentiel de la pharmacologie quantitative des systèmes appliquées au développement du médicament et à la médecine translationnelle, en se servant du modèle physio-PK/PD que nous avons mis au point.

Détermination de la spéciation du samarium dans l'environnement

Les éléments de terres rares (REEs) sont de plus en plus utilisés dans une multitude d’applications, notamment la fabrication d’aimants, de batteries rechargeables et les écrans de téléviseurs. Ils sont pour la plupart des métaux trivalents peu solubles dans les eaux naturelles. Comme pour les métaux divalents, le risque écologique des REEs est très probablement étroitement lié à leurs spéciations chimiques. Cependant, le comportement du samarium (Sm) dans les matrices environnementales est très peu connu et il n'existe actuellement aucune technique pour évaluer sa spéciation chimique. Dans cette optique, la technique d'échange d'ions (IET) sur la résine Dowex a été optimisée pour mesurer le samarium libre en solution. Les temps d'équilibre ont d'abord été déterminés pour des solutions tamponnées de samarium (Sm 6,7x10-8 M ; MES 1,0 mM M ; pH 6,0) en présence du nitrate de sodium (de 0,01M à 0,5 M). Pour ces diverses forces ioniques, l’équilibre thermodynamique n’est atteint que pour NaNO3 0,5M. Un autre mode d’utilisation de la résine (mode dynamique) a donc été développé pour tenir compte des conditions environnementales et évaluer efficacement le samarium libre. Les impacts des ligands organiques tels le NTA, l’EDTA, le citrate, l’acide malique et l’acide fulvique Suwannee River Standard I (SRFA) ont été étudiés par l’IET en mode dynamique. Une grande corrélation a été trouvée plus entre les taux d’accumulation de samarium sur la résine d’échange pour différents rapports NTA : Sm, EDTA : Sm, SRFA : Sm et le samarium libre. Par contre, aucune corrélation significative n’a été observée pour les ligands citrate et acide malique compte tenu des complexes qu’ils forment avec le samarium et qui s’adsorbent aussi sur la résine Dowex. Les concentrations Sm3+ mesurées par la technique IET ont été fortement corrélées avec celles prédites par le modèle WHAM 7.0 en utilisant la constante de stabilité obtenue par titration de SRFA par extinction de la fluorescence. Par ailleurs, la formation de colloïdes de samarium en fonction du pH influe grandement sur la détermination du samarium libre et doit être prise en compte dans la spéciation du samarium. L'IET assisté par des techniques auxiliaires comme le dosage par extinction de la fluorescence et le SP-ICPMS pourrait être une technique utile pour évaluer les concentrations de Sm biodisponible dans les eaux naturelles.