Optimization models and methods for real-time transportation planning in forestry

Lors du transport du bois de la forêt vers les usines, de nombreux événements imprévus peuvent se produire, événements qui perturbent les trajets prévus (par exemple, en raison des conditions météo, des feux de forêt, de la présence de nouveaux chargements, etc.). Lorsque de tels événements ne sont connus que durant un trajet, le camion qui accomplit ce trajet doit être détourné vers un chemin alternatif. En l’absence d’informations sur un tel chemin, le chauffeur du camion est susceptible de choisir un chemin alternatif inutilement long ou pire, qui est lui-même "fermé" suite à un événement imprévu. Il est donc essentiel de fournir aux chauffeurs des informations en temps réel, en particulier des suggestions de chemins alternatifs lorsqu’une route prévue s’avère impraticable. Les possibilités de recours en cas d’imprévus dépendent des caractéristiques de la chaîne logistique étudiée comme la présence de camions auto-chargeurs et la politique de gestion du transport. Nous présentons trois articles traitant de contextes d’application différents ainsi que des modèles et des méthodes de résolution adaptés à chacun des contextes. Dans le premier article, les chauffeurs de camion disposent de l’ensemble du plan hebdomadaire de la semaine en cours. Dans ce contexte, tous les efforts doivent être faits pour minimiser les changements apportés au plan initial. Bien que la flotte de camions soit homogène, il y a un ordre de priorité des chauffeurs. Les plus prioritaires obtiennent les volumes de travail les plus importants. Minimiser les changements dans leurs plans est également une priorité. Étant donné que les conséquences des événements imprévus sur le plan de transport sont essentiellement des annulations et/ou des retards de certains voyages, l’approche proposée traite d’abord l’annulation et le retard d’un seul voyage, puis elle est généralisée pour traiter des événements plus complexes. Dans cette ap- proche, nous essayons de re-planifier les voyages impactés durant la même semaine de telle sorte qu’une chargeuse soit libre au moment de l’arrivée du camion à la fois au site forestier et à l’usine. De cette façon, les voyages des autres camions ne seront pas mo- difiés. Cette approche fournit aux répartiteurs des plans alternatifs en quelques secondes. De meilleures solutions pourraient être obtenues si le répartiteur était autorisé à apporter plus de modifications au plan initial. Dans le second article, nous considérons un contexte où un seul voyage à la fois est communiqué aux chauffeurs. Le répartiteur attend jusqu’à ce que le chauffeur termine son voyage avant de lui révéler le prochain voyage. Ce contexte est plus souple et offre plus de possibilités de recours en cas d’imprévus. En plus, le problème hebdomadaire peut être divisé en des problèmes quotidiens, puisque la demande est quotidienne et les usines sont ouvertes pendant des périodes limitées durant la journée. Nous utilisons un modèle de programmation mathématique basé sur un réseau espace-temps pour réagir aux perturbations. Bien que ces dernières puissent avoir des effets différents sur le plan de transport initial, une caractéristique clé du modèle proposé est qu’il reste valable pour traiter tous les imprévus, quelle que soit leur nature. En effet, l’impact de ces événements est capturé dans le réseau espace-temps et dans les paramètres d’entrée plutôt que dans le modèle lui-même. Le modèle est résolu pour la journée en cours chaque fois qu’un événement imprévu est révélé. Dans le dernier article, la flotte de camions est hétérogène, comprenant des camions avec des chargeuses à bord. La configuration des routes de ces camions est différente de celle des camions réguliers, car ils ne doivent pas être synchronisés avec les chargeuses. Nous utilisons un modèle mathématique où les colonnes peuvent être facilement et naturellement interprétées comme des itinéraires de camions. Nous résolvons ce modèle en utilisant la génération de colonnes. Dans un premier temps, nous relaxons l’intégralité des variables de décision et nous considérons seulement un sous-ensemble des itinéraires réalisables. Les itinéraires avec un potentiel d’amélioration de la solution courante sont ajoutés au modèle de manière itérative. Un réseau espace-temps est utilisé à la fois pour représenter les impacts des événements imprévus et pour générer ces itinéraires. La solution obtenue est généralement fractionnaire et un algorithme de branch-and-price est utilisé pour trouver des solutions entières. Plusieurs scénarios de perturbation ont été développés pour tester l’approche proposée sur des études de cas provenant de l’industrie forestière canadienne et les résultats numériques sont présentés pour les trois contextes.

Novel surface-tethered estrogen polymeric platforms in cardiovascular regenerative medicine

L’estradiol (E2) est une hormone femelle qui joue un rôle essentiel, à la fois dans la régulation et dans la détermination de certaines conditions physiologiques in vivo, telle que la différenciation et la prolifération cellulaire. Lorsque l’E2 est donné en supplément, par exemple dans le cas de thérapie hormonale, deux effets sont observés, un effet génomique et un effet non-génomique, de par son interaction avec les récepteurs à œstrogène du noyau ou de la membrane cellulaire, respectivement. L’effet non-génomique est plus difficile à étudier biologiquement parce que l’effet se produit sur une échelle de temps extrêmement courte et à cause de la nature hydrophobe de l’E2 qui réduit sa biodisponibilité et donc son accessibilité aux cellules cibles. C’est pourquoi il est nécessaire de développer des systèmes d’administration de l’E2 qui permettent de n’étudier que l’effet non-génomique de l’œstrogène. Une des stratégies employée consiste à greffer l’E2 à des macromolécules hydrophiles, comme de l’albumine de sérum bovin (BSA) ou des dendrimères de type poly(amido)amine, permettant de maintenir l’interaction de l’E2 avec les récepteurs d’œstrogène de la membrane cellulaire et d’éviter la pénétration de l’E2 dans le noyau des cellules. Toutefois, ces systèmes macromolécules-E2 sont critiquables car ils sont peu stables et l’E2 peut se retrouver sous forme libre, ce qui affecte sa localisation cellulaire. L’objectif de cette thèse est donc de développer de nouvelles plateformes fonctionnalisées avec de l’E2 en utilisant les approches de synthèses ascendantes et descendantes. Le but de ces plateformes est de permettre d’étudier le mécanisme de l’effet non-génomique de l’E2, ainsi que d’explorer des applications potentielles dans le domaine biomédical. L’approche ascendante est basée sur un ligand d’E2 activé, l’acide 17,α-éthinylestradiol-benzoïque, attaché de façon covalente à un polymère de chitosan avec des substitutions de phosphorylcholine (CH-PC-E2). L’estradiol est sous forme de pro-drogue attachée au polymère qui s’auto-assembler pour former un film. L’effet biologique de la composition chimique du film de chitosan-phosphorylcholine a été étudié sur des cellules endothéliales. Les films de compositions chimiques différentes ont préalablement été caractérisés de façon physicochimique. La topographie de la surface, la charge de surface, ainsi que la rhéologie des différents films contenant 15, 25, ou 40% molaires de phosphorylcholine, ont été étudiés par microscopie à force atomique (AFM), potentiel zêta, résonance plasmonique de surface et par microbalance à cristal de quartz avec dissipation (QCM-D). Les résultats de QCM-D ont montré que plus la part molaire en phosphorylcholine est grande moins il y a de fibrinogène qui s’adsorbe sur le film de CH-PC. Des cellules humaines de veine ombilicale (HUVECs) cultivées sur des films de CH-PC25 et de CH-PC40 forment des amas cellulaire appelés sphéroïdes au bout de 4 jours, alors que ce n’est pas le cas lorsque ces cellules sont cultivées sur des films de CH-PC15. L’attachement de l’estradiol au polymère a été caractérisé par plusieurs techniques, telles que la résonance magnétique nucléaire de proton (1H NMR), la spectroscopie infrarouge avec transformée de Fourier à réfraction totale atténuée (FTIR-ATR) et la spectroscopie UV-visible. La nature hydrogel des films (sa capacité à retenir l’eau) ainsi que l’interaction des films avec des récepteurs à E2, ont été étudiés par la QCM-D. Des études d’imagerie cellulaires utilisant du diacétate de diaminofluoresceine-FM ont révélé que les films hydrogels de CH-PC-E2 stimulent la production d’oxyde nitrique par les cellules endothéliales, qui joue un rôle protecteur pour le système cardiovasculaire. L’ensemble de ces études met en valeur les rôles différents et les applications potentielles qu’ont les films de type CH-PC-E2 et CH-PC dans le cadre de la médecine cardiovasculaire régénérative. L’approche descendante est basée sur l’attachement de façon covalente d’E2 sur des ilots d’or de 2 μm disposés en rangées et espacés par 12 μm sur un substrat en verre. Les ilots ont été préparés par photolithographie. La surface du verre a quant à elle été modifiée à l’aide d’un tripeptide cyclique, le cRGD, favorisant l’adhésion cellulaire. L’attachement d’E2 sur les surfaces d’or a été suivi et confirmé par les techniques de SPR et de QCM-D. Des études d’ELISA ont montré une augmentation significative du niveau de phosphorylation de la kinase ERK (marqueur important de l’effet non-génomique) après 1 heure d’exposition des cellules endothéliales aux motifs alternant l’E2 et le cRGD. Par contre lorsque des cellules cancéreuses sont déposées sur les surfaces présentant des motifs d’E2, ces cellules ne croissent pas, ce qui suggère que l’E2 n’exerce pas d’effet génomique. Les résultats de l’approche descendante montrent le potentiel des surfaces présentant des motifs d’E2 pour l’étude des effets non-génomiques de l’E2 dans un modèle in vitro.

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.