FR 2998480 A1 20140530 FR 1261339 A 20121128 « Procédé pour contrôler l'activation d'une substance photosensibilisante dans un tissu biologique, et dispositif mettant en oeuvre le procédé » Domaine technique La présente invention concerne un procédé pour contrôler l'activation d'une substance photo-sensibilisante dans un tissu biologique. Elle concerne également un dispositif permettant de mettre en oeuvre le procédé. Le domaine de l'invention est plus particulièrement mais de manière non 10 limitative celui des thérapies photodynamiques, en particulier pour le traitement des cellules cancéreuses. Etat de la technique antérieure Les thérapies photodynamiques (PDT) sont des techniques utilisées pour traiter des tissus biologiques et éliminer sélectivement des cellules. Elles sont 15 utilisées notamment pour traiter des lésions cancéreuses et précancéreuses, en particulier de la peau. Elles sont également utilisées pour traiter d'autres affections telles que l'acné kystique, les kératoses actiniques, les tâches pigmentaires, la dégénération maculaire liées à l'âge, les rugosités, les vaisseaux dilatés, le psoriasis, ou pour des traitements anti-infectieux, des 20 petites rides, ou l'épilation. Le principe des thérapies photodynamiques est le suivant : - dans un premier temps, un produit photosensibilisant (PS) est introduit dans les tissus biologiques ; - la vitesse d'élimination du produit photosensibilisant étant variable en 25 fonction de la nature des cellules et de leur environnement (vascularisation, ...), après un certain temps le produit restant dans le tissu biologique est essentiellement concentré autour des cellules cibles (par exemple dans la tumeur); - le tissus biologique et le produit photosensibilisant sont alors illuminés 30 par une lumière à une longueur d'onde appropriée ; - le produit est conçu de telle sorte que sous l'effet de l'illumination il est « activé » et génère des espèces réactives de l'oxygène, tels que de l'oxygène singulet, à haute dose ; - ces espèces réactives de l'oxygène ont pour effet de détruire ou de 35 dégrader les molécules voisines. Ainsi, si le produit est concentré plus -2- spécifiquement dans certaines zones du tissu biologique, ces zones sont sélectivement endommagées ou détruites. Les produits photosensibilisants utilisés peuvent par exemple être des composés de la famille des porphyrines, ou des chlorines. Les applications de la thérapie photodynamique concernent en particulier des affections et des cancers de la peau, ou des traitements de la dégénérescence maculaire due à l'âge. Dans ce cas les tissus sont facilement accessibles pour l'illumination. La technique peut également être utilisée pour le traitement de tumeurs 10 interstitielles, c'est-à-dire dans les tissus. Dans ce cas l'illumination peut être effectuée au moyen de fibres optiques introduites dans les tissus. Un inconvénient majeur des thérapies photodynamiques, qui en limite le développement pour des applications cliniques, est leur faible reproductibilité in-vivo. En effet, le comportement et la réponse des tissus au traitement 15 varient fortement d'un sujet à l'autre, ce qui entraine une grande incertitude et une grande variabilité sur les réponses thérapeutiques. De manière habituelle, les traitements consistent à illuminer la zone à traiter pendant une durée prédéterminée, suivant une ou plusieurs séquences. Il n'y a donc pas de contrôle de l'efficacité en cours de traitement. 20 Les applications cliniques actuelles de la thérapie photodynamique reposent sur un modèle macroscopique statique simplifié du phénomène photophysique inhérent pour déterminer la concentration de produit photosensibilisant nécessaire et d'éclairement énergétique ou irradiance à appliquer : 25 [R] = IcsbapAT[So]d)f , (Eq- 1) où [R] est une concentration d'espèces réactives de l'oxygène qui, si elles sont situées dans une région sensible des cellules cancéreuses, déclenchent une cascade de réactions conduisant à la mort de la cellule. (pA.A désigne l'irradiance ou la puissance lumineuse par unité de surface reçue par la 30 surface du tissu à traiter et ÀA correspond à une longueur d'onde d'absorption du PS. T est le temps d'exposition à la lumière et [S0] est la concentration intracellulaire d'agent photosensibilisant. Ics est le coefficient de rétro-diffusion du tissu, b est un rendement de conversion, E est le coefficient d'extinction de l'agent photosensibilisant et t est un rendement quantique de conversion du 35 médicament en radicaux oxydants. f désigne la fraction d'espèces réactives -3- de l'oxygène générées attaquant les cibles sensibles des cellules, tandis que l'autre fraction (1- f) ne cause que des effets mineurs. En dépit de son utilisation actuelle dans les applications cliniques, ce modèle est sujet à plusieurs controverses : - l'Eq. 1 montre une simple réciprocité entre la concentration de l'agent photosensibilisant et la lumière, alors que plusieurs expériences ont montré des résultats contradictoires ; - le rendement quantique t est en réalité une fonction du taux d'oxygénation de la tumeur. Ce dernier est souvent variable dans le temps et 10 reste difficile à mesurer durant le traitement ; - les sites intracellulaires endommagés par la PDT dépendent beaucoup de la localisation du PS dans les cellules. Les sites d'action de la PDT sont principalement les mitochondries, le reticulum endoplasmique, l'appareil de Golgi, les lysosomes et les lipides membranaires. Une re-localisation du PS 15 peut se produire pendant la PDT. Certains sites sont plus critiques que d'autres mais ce facteur n'est pas pris en compte dans l'Eq. 1; - plusieurs quantités impliquées dans cette équation sont en réalité variables dans le temps. En particulier, il a été montré qu'une réduction de la concentration du PS durant le traitement était un facteur non négligeable, 20 mais absent de cette formule. Plus généralement, l'Eq. 1 n'est qu'une description statique du processus dynamique de la PDT. Il est donc nécessaire de mettre au point des méthodes de dosimétrie pour améliorer la reproductibilité des réponses thérapeutiques de la PDT. Une telle méthode de dosimétrie doit permettre de déterminer un dosage des 25 facteurs incidents (tels que la concentration locale en produit photosensibilisant, la densité des photons envoyés, le temps d'illumination, l'intervalle entre l'administration du produit et les phases d'illumination), pour obtenir un effet souhaité (par exemple la destruction sélective de certaines cellules). 30 La détermination d'une méthode de dosimétrie précise constitue une voie de recherche primordiale. Il existe plusieurs classes de dosimétrie potentiellement utilisables en PDT, et toutes reposent sur un modèle : La dosimétrie explicite Ce type de technique implique la mesure des principaux composants de 35 la thérapie, c'est-à-dire les doses de lumière, d'oxygène et la concentration -4- intracellulaire de l'agent photosensibilisant. Ces trois constituants sont regroupés dans un modèle de dose des espèces réactives de l'oxygène produites, défini par l'Eq. 1. Cependant, aucun de ces éléments n'est mesurable facilement. En outre, d'autres inconvénients déjà évoqués précédemment, comme la variation des doses des trois éléments au cours du traitement, constituent autant de limites au développement de cette approche. La dosimétrie implicite Ce type d'approche propose d'utiliser une seule grandeur dont l'évolution dépend de l'effet conjugué des trois variables de base de la PDT (doses de lumière, de produit administré et d'oxygène). Cette grandeur doit être en mesure de renseigner sur les dommages biologiques prévisibles et donc sur les issues thérapeutiques du traitement. Une variable intéressante pour jouer ce rôle est la concentration intratumorale de l'agent photosensibilisant (PS). Une technique de mesure connue consiste à utiliser une méthode invasive, nécessitant l'exérèse de tumeurs xénogreffées chez des souris « nude » (avec un système immunitaire déficient) ayant reçu l'agent photosensibilisant, puis d'estimer les concentrations du PS avant et après illumination par chromatographie liquide haute performance. Cette approche fournit des résultats précis sur la concentration du PS et des photoproduits. Mais elle ne permet pas d'observer en temps réel l'évolution du PS durant le traitement. La mesure des réponses biologiques/biophysiques des tissus Il s'agit d'une approche alternative pour suivre l'évolution du traitement 25 et prédire les dommages sur les tissus comme les dommages vasculaires et les nécroses induites par le traitement. Les méthodes utilisées dans ce but incluent la spectroscopie par impédance électrique, la mesure du flux sanguin par laser Doppler, l'analyse de la diffusion par spectroscopie de corrélation de fluorescence, la 30 tomographie Doppler par cohérence optique, aussi bien que les méthodes d'imagerie moléculaire ou radiologique telles que la tomographie à émission de positron (PET). Il est encore difficile de savoir si ces techniques pourront prédire efficacement les réponses thérapeutiques de la PDT car elles nécessitent en particulier des changements rapides et significatifs des 35 caractéristiques biologiques et spectrales des tissus. -5- On connaît par exemple le document EP 2 431 072 Al qui divulgue un procédé de contrôle de l'illumination en thérapie photodynamique de tumeur interstitielle basé sur une mesure de propriétés optiques des tissus et une méthode de calcul pour déterminer l'état du tissu pendant le traitement. La dosimétrie directe Cette technique implique la mesure directe des concentrations des espèces réactives de l'oxygène et en particulier l'oxygène singulet. Comme pour l'approche de dosimétrie implicite, cette stratégie ramène le problème de dosimétrie à une seule dimension, celle de l'oxygène singulet mais directement cette fois. Plusieurs travaux ont montré des premiers résultats in vivo en laboratoire mais l'instrumentation actuelle requise pour ces travaux ne permet pas d'envisager des applications cliniques dans l'immédiat. Un objet de la présente invention est de proposer une méthode d'ajustement automatique du rayonnement électromagnétique d'activation de 15 l'agent photosensibilisant pendant une session de traitement en thérapie photodynamique in-vivo, et qui résolve les inconvénients de l'art antérieur. Un autre objet de la présente invention est de proposer une méthode d'ajustement automatique du rayonnement électromagnétique d'activation de l'agent photosensibilisant qui permette de contrôler de manière reproductible 20 l'application d'un traitement en thérapie photodynamique. Un autre objet de la présente invention est de proposer une méthode d'ajustement automatique du rayonnement électromagnétique d'activation de l'agent photosensibilisant en thérapie photodynamique qui puisse être mise en oeuvre avec une instrumentation simple et adaptée aux applications in-vivo. 25 Exposé de l'invention Cet objectif est atteint avec un procédé pour contrôler l'activation d'une substance photosensibilisante préalablement introduite dans un tissu biologique, comprenant une étape d'exposition de ladite substance photosensibilisante à un rayonnement électromagnétique d'activation de telle 30 sorte à l'activer, caractérisé en ce qu'il comprend en outre des étapes : - d'obtention d'une information de luminescence issue de ladite substance photosensibilisante, et - d'ajustement de l'exposition de la substance photosensibilisante au 35 rayonnement électromagnétique d'activation de telle sorte que ladite -6- information de luminescence suive une évolution temporelle correspondant à une consigne temporelle prédéterminée. Il est entendu que le terme de luminescence regroupe tous modes d'émission de lumière dits « froids », tels que par exemple la 5 chimiluminescence, la phosphorescence ou la fluorescence. La substance photosensibilisante peut être introduite dans le tissu biologique notamment par injection, ou par simple contact dans le cas d'une surface accessible telle que la peau. La consigne temporelle prédéterminée peut comprendre toutes sortes de 10 profils de courbes, décrites par exemple par des équations (droites, courbes linéaires par morceaux, polynômes ...) ou des tableaux de points. Cette consigne temporelle peut par exemple être issue de protocoles de mesure. Le rayonnement électromagnétique d'activation peut comprendre un 15 rayonnement ionisant (rayons X), et/ou un rayonnement non ionisant. Suivant des modes de réalisation, le rayonnement électromagnétique d'activation peut comprendre un rayonnement lumineux d'activation avec des longueurs d'ondes optiques. Suivant des modes de réalisation, le procédé selon l'invention peut 20 comprendre en outre une étape de mesure d'une information de luminescence issue de l'oxygène singulet généré par l'activation de la substance photosensibilisante. Cette information de luminescence peut être notamment issue d'une phosphorescence à une longueur d'onde de 1,27 dam de l'oxygène singulet qui 25 est généré par des substances photosensibilisantes telles que les porphyrines ou les chlorines. Suivant des modes de réalisation, le procédé selon l'invention peut comprendre en outre des étapes : - d'exposition de la substance photosensibilisante à un rayonnement 30 lumineux d'activation appliqué sous forme d'impulsions lumineuses d'activation, et - de mesure de luminescence entre lesdites impulsions. Cette exposition de la substance photosensibilisante à un rayonnement lumineux d'activation appliqué sous forme d'impulsions lumineuses séparées 35 par des périodes d'obscurité permet ainsi d'effectuer des mesures de -7- luminescence en temps réel, au cours du traitement, et y compris dans les cas où de telles mesures ne sont pas possibles en présence du rayonnement lumineux d'activation. La mesure de luminescence peut comprendre des étapes : - d'illumination de la substance photosensibilisante avec une lumière d'excitation à une longueur d'onde d'excitation, - de mesure d'une intensité de fluorescence. Dans ce cas, la mesure de luminescence comprend une mesure de fluorescence effectuée sur la substance photosensibilisante. L'évolution temporelle de cette fluorescence permet d'évaluer le phénomène de photoblanchiment de la substance photosensibilisante qui résulte de sa dégradation. Le procédé selon l'invention peut comprendre en outre une étape de calcul d'un observateur de luminescence utilisant en entrée au moins une 15 mesure de luminescence et une commande des impulsions lumineuses d'activation, et produisant en sortie une estimation de la luminescence. De préférence, ces mesures de luminescence comprennent des mesures de fluorescence effectuées sur la substance photosensibilisante. L'observateur de luminescence peut être modélisé sous la forme d'un 20 observateur d'état comprenant une variable d'état correspondant à une estimation de la luminescence (ou de la fluorescence) associée à l'effet seul du traitement. Cet observateur d'état peut être dérivé d'un modèle autorégressif avec entrées externes comme la structure de modèle ARX (« Auto Regressive 25 model with eXternal inputs » en Anglais) ou un modèle de type représentation d'état. L'observateur de luminescence peut comprendre des paramètres qui sont déterminés lors d'une étape préalable de calibration à partir de mesures expérimentales. 30 L'observateur de luminescence permet de réaliser une interpolation intelligente des mesures de luminescence ou de fluorescence. En effet, en pratique ces mesures ne peuvent pas être réalisées pendant les phases d'impulsions du rayon électromagnétique car ce dernier aveugle de par sa grande intensité le système de mesure de luminescence. 35 L'observateur de luminescence permet aussi de réduire très -8- sensiblement l'intensité du bruit de mesure de luminescence. L'observateur de luminescence permet ainsi d'améliorer considérablement le fonctionnement d'une boucle de régulation de la lumière d'activation. En outre, l'observateur de luminescence est déterminé à partir de paramètres expérimentaux relativement faciles à obtenir. Une grande précision n'est pas nécessairement requise puisque l'observateur de luminescence peut être recalé périodiquement par les mesures de luminescence. Ainsi, ses paramètres peuvent être déterminés par exemple pour un type de tissu biologique donné, et applicables directement à différents sujets traités. Le procédé selon l'invention peut comprendre en outre une étape de comparaison de l'estimation de la luminescence issue de l'observateur de luminescence avec la consigne temporelle de telle sorte à générer un signal d'erreur, lequel étant utilisé pour l'ajustement de l'exposition au rayonnement électromagnétique d'activation. L'ajustement de l'exposition au rayonnement électromagnétique d'activation peut comprendre un ajustement de l'amplitude et/ou de la durée des impulsions lumineuses. Ainsi, le procédé selon l'invention permet un asservissement en boucle fermée de l'exposition au rayonnement électromagnétique d'activation de telle sorte qu'une variable de mesure (la luminescence mesurée ou estimée par l'observateur d'état) suive une trajectoire prédéterminée (la consigne temporelle). Le procédé selon l'invention peut être est mis en oeuvre avec une substance photosensibilisante dont l'activation produit des substances aptes à détruire des cellules. Il peut notamment être est mis en oeuvre avec des substances photosensibilisantes telles que les porphyrines ou les chlorines, qui libèrent de l'oxygène singulet. Suivant un autre aspect, il est proposé un dispositif pour contrôler l'activation d'une substance photosensibilisante préalablement introduite dans un tissu biologique, comprenant : - des moyens d'exposition aptes à permettre une exposition de ladite substance photosensibilisante à un rayonnement électromagnétique 35 d'activation de telle sorte à l'activer, -9- - des moyens d'acquisition aptes à produire une information de luminescence issue de ladite substance photosensibilisante, et - des moyens électroniques et de calcul aptes à contrôler les moyens d'exposition de telle sorte à permettre un ajustement de l'exposition de la substance photosensibilisante au rayonnement électromagnétique d'activation de telle sorte que ladite information de luminescence suive une évolution temporelle correspondant à une consigne temporelle prédéterminée. Les moyens d'exposition peuvent comprendre une source lumineuse d'activation apte à produire des impulsions de lumière en fonction d'un signal 10 de commande. Les moyens d'acquisition peuvent comprendre un spectromètre ou un spectrofluorimètre. Les moyens électroniques et de calcul peuvent implémenter une boucle de contrôle comprenant un observateur de luminescence, 15 - lequel observateur de luminescence utilisant en entrée au moins une mesure de luminescence et un signal de commande de la source lumineuse d'activation, et produisant en sortie une estimation de la luminescence, - laquelle estimation de la luminescence étant utilisée pour générer un signal de commande de la source lumineuse d'activation. 20 Description des figures et modes de réalisation D'autres avantages et particularités de l'invention apparaîtront à la lecture de la description détaillée des mises en oeuvre et des modes de réalisation nullement limitatifs, et des dessins annexés suivants : - la Fig. 1 présente des résultats de simulation d'un processus de 25 photoréaction, avec, Fig. 1(a) la concentration d'agent photosensibilisant PS à l'état fondamental (variable de sortie), Fig. 1(b) le signal d'irradiance (variable d'entrée), Fig. 1(c) la concentration d'oxygène singulet (variable d'état), et Fig. 1(d) la quantité cumulée d'oxygène singulet produite, - la Fig. 2 présente un schéma synoptique de l'invention, 30 - la Fig. 3 présente un organigramme du procédé selon l'invention, - La Fig. 4 illustre une comparaison entre des mesures de photoblanchiment in-vitro et une estimation du photoblanchiment obtenue avec le modèle d'état selon l'invention, - la Fig. 5 illustre les performances de l'observateur de 35 photoblanchiment, au travers de mesures simulées, -10- - la Fig. 6 présente des résultats de l'asservissement du photoblanchiment, en conditions in-vitro. Modèle cytotoxique Le processus de réaction à la base des thérapies photodynamiques est le 5 suivant : - un photosensibiliseur (PS) est introduit dans les tissus, par exemple par voie intraveineuse ou même par application directe, pour des applications sur la peau. Dans le cadre des modes de réalisation présentés, on considère que ce PS peut être par exemple une porphyrine ou une hématoprophyrine ; 10 - le PS se concentre plus ou moins sélectivement dans les tissus à détruire (par exemple les tissus tumoraux), ou il s'en élimine plus lentement ce qui produit le même résultat ; - la zone à traiter est ensuite illuminée avec une lumière d'activation qui contient les longueurs d'onde nécessaires au déclanchement des réactions. 15 Idéalement, ces longueurs d'onde correspondent à des bandes d'absorption du PS, ou du moins à un compromis entre les longueurs d'onde aptes à pénétrer les tissus biologiques (dans le rouge) et ces bandes d'absorption. Par exemple, pour des hématoporphyrines, l'illumination est souvent effectuée à des longueurs d'onde de l'ordre de 630-650 nm ; 20 - l'action de la lumière sur le PS entraîne la formation d'oxygène singulet 102. Cette molécule est une forme d'oxygène de plus haute énergie que l'oxygène classique, qui est formée par excitation de ce dernier. C'est une molécule à courte durée de vie qui est très réactive vis-à-vis des composantes cellulaires et donc très toxique. 25 L'ensemble de ces processus entraîne la dégradation d'un grand nombre de constituants cellulaires, ce qui conduit à la mort des cellules par nécrose ou par apoptose. Le processus photochimique à l'origine de la formation d'oxygène singulet est le suivant : 30 - le PS absorbe des photons de la lumière d'activation ; - il passe alors d'un état singulet fondamental So à un état singulet excité S1 ; - le retour à l'état initial So peut s'opérer notamment en passant à un état triplet T1 par conversion intersystème (ISC), donc une transition non 35 radiative ; - la durée de vie de cet état triplet T1 étant relativement longue, la molécule de PS peut alors interagir avec des substrats situés dans son environnement proche, comme par exemple l'oxygène présent dans les tissus. Deux types de réactions, connues sous les noms de photoréactions de type I et II, sont possibles. Les photoréactions de type II sont prédominantes en PDT avec des PS de type porphyrines ; - dans une photoréaction de type II, l'état triplet T1 du PS réagit directement avec l'oxygène présent en lui transférant son excès d'énergie, ce qui fait passer l'oxygène de son état initial 302 à son état singulet 102. Il est à noter que l'oxygène singulet formé peut également réagir avec le PS et le détruire. Cette dégradation, appelée photoblanchiment, est une conséquence des réactions de type I et II. Ce photoblanchiment conduit à une diminution de l'absorbance du PS, et donc à une diminution de son efficacité. Le photoblanchiment met en jeu deux mécanismes : - la photodégradation, qui correspond à une modification profonde de la structure du PS conduisant à la formation de photoproduits n'absorbant peu ou pas la lumière d'activation; - la phototransformation qui correspond à une perte de l'absorbance et de la fluorescence du photosensibilisateur à certaines longueurs d'onde mais 20 le chromophore est conservé. Par ailleurs, si l'on irradie un PS à une longueur d'onde d'excitation adéquate, il est possible d'observer une fluorescence. Cette fluorescence correspond à un retour du PS à un état singulet excité S1 vers l'état fondamental So par émission de photons de fluorescence. 25 Cette fluorescence peut notamment être utilisée pour déterminer la localisation du photosensibilisateur dans les cellules, et pour vérifier sa complète élimination après traitement. L'intensité de fluorescence est également affectée par le photoblanchiment. Ainsi, le phénomène de photoblanchiment peut être évalué 30 par la mesure de la décroissance de l'intensité de fluorescence photoinduite du PS. Dans le cadre de l'invention, un modèle non linéaire des processus de photoréaction à la base des thérapies photodynamiques (PDT) a été développé et décrit par une représentation d'état. Ce modèle a été obtenu à 35 partir d'équations cinétiques, connues dans la littérature, qui décrivent les (Eq. 2) dt [M] = [l] 0 Y = Y [S o] + V .7V'(0, s2) Les variables du modèle sont les suivantes: [Se] : concentration de PS à l'état singulet fondamental (millimolaire, mM) ; : concentration de PS à l'état singulet excité (mM) ; [Tl] : concentration de PS à l'état triplet excité (mM) ; [ 302] : concentration d'oxygène à l'état triplet initial (mM) ; [1-02] : concentration d'oxygène à l'état singulé excité (mM) ; [M] : concentration de cellules organiques (mM) ; Up : vitesse d'absorption des molécules de PS dans les tissus (mM.s-1) , uL : éclairement énergétique de l'illumination (W-cm-2 ); uo : vitesse d'absorption des molécules d'oxygène (mM.s ); y : intensité de fluorescence (unités arbitraires, u.a.). x(t) = [[So], [Tl], [ 302], [ 102], [M]]T E R6 sont les variables d'état du modèle. La pénétration des molécules photosensibilisantes (PS) dans les cellules est un processus avec une dynamique très lente, comparativement aux photoréactions. Aussi, la vitesse d'absorption des molécules photosensibilisantes Up peut être considérée comme nulle (Up = 0) pendant la durée du traitement, qui correspond à une période de temps relativement courte, inférieure à 1 heure. La vitesse d'absorption de molécules d'oxygène uo(t) peut être supposée constante durant le traitement. Le bruit de mesure v(ti) peut être considéré comme étant une séquence de réalisation de variables gaussiennes indépendantes et uniformément distribuées aux instants de mesure ftil, selon une loi normale .7v-(0,62). -12- réactions de type II spécifiques à la PDT : d[so] Up + kf - kpb[Se][ 102] + k p [Ti] + k [ 302] - kAUL [S0] dt dt dt 3 d[Ti] k sC [S11 k p [Ti] - C [Ti] 02] - kTM Ti [M] d[ 3 2 0] 110 kT [T1] [ 302] + kd 102] dt d[S1] kAUL[So] kISC)[S1] kSM[S1] [M] d[ 102] kT [Ti] [ 302] - k1[ 1-02] -kox [M] [ 102] kPb [S0 2 ][ 1°] -13- Le modèle de l'Eq. 2 met également en oeuvre un ensemble de paramètres qui sont listés ci-après, avec les valeurs et les conditions initiales utilisées pour effectuer les calculs : Symbole Définition Valeur kA Constante d'absorption du photosensibiliseur pour un éclairement uL = 76mW-cm-2 1 kpb Constante de photoblanchiment 1,6.106 mM-is-I. kT Constante de vitesse bimoléculaire pour la réaction de 302 avec T1 105 mM-is-1 kf Constante de fluorescence pour la réaction 2.107 s-1 Sl So kP Constante de phosphorescence pour la 1250 s-1 réaction T1 So k15 Constante de vitesse pour la réaction Sl T1 8.107 s-1 k1 Constante de luminescence pour la réaction 105 s-1 102 302 k0 Constante bimoléculaire de la réaction d'oxydation entre 102 et M 105 mM-is-1 Y Coefficient de mesure 43,15 u.a. [Solo Condition initiale pour [So] 8,5.10-3 mM [S1]0 Condition initiale pour [Si] 0 mM [T1]0 Condition initiale pour [T1] 0 mM [3010 Condition initiale pour [302] 83.10-3 mM [102]0 Condition initiale pour [102] 0 mM [M]o Condition initiale pour [M] 830.10-3 mM La Fig. 1 présente les résultats d'une simulation numérique d'un processus de photoréaction avec un agent photosensibilisant présent dans un tissu biologique. Cette simulation est calculée en implémentant le modèle de l'Eq. 2. Les axes des abscisses correspondent au temps t(s). Le signal lumineux d'activation utilisé pour illuminer les tissus biologiques et le PS est présenté à la Fig. 1(b). Il s'agit d'un signal carré de période 60s, de rapport cyclique 1/2 et d'irradiance ou d'éclairement -14- énergétique uL=95 mW-cm-2. La Fig. 1(a) illustre l'évolution temporelle de la concentration [S0] de l'agent photosensibilisant à l'état singulet fondamental. Cette courbe correspond également à la courbe de photoblanchiment du PS. La Fig. 1(c) illustre l'évolution de la concentration [1021 d'oxygène singulet produit par l'activation du PS par le signal lumineux d'activation. La Fig. 1(d) illustre l'évolution de la quantité cumulée d'oxygène singulet produite par l'activation du PS par le signal lumineux d'activation. Dans la mesure où l'oxygène singulet est l'espèce photocytotoxique véritablement à l'origine des dommages biologiques, une mesure de sa concentration au cours du traitement peut permettre de contrôler avec une grande précision l'application du traitement et ainsi de développer des techniques de dosimétrie capables d'amener une meilleure reproductibilité des résultats. Une mesure directe de la concentration d'oxygène singulet peut être effectuée notamment en mesurant la luminescence produite à une longueur d'onde de 1,27 pm par cet oxygène singulet. Cependant, en pratique, cette concentration d'oxygène singulet est très difficilement mesurable dans des conditions in-vivo. Les simulations de la Fig. 1 montrent également que la courbe de photoblanchiment du PS (Fig. 1(a)) est un bon indicateur de la production d'oxygène singulet (Fig. 1(d)). Ainsi suivant un mode de réalisation préférentiel de l'invention, c'est le signal de photoblanchiment qui est utilisé comme marqueur temporel de l'efficacité du traitement. Comme expliqué précédemment, l'invention a pour de proposer une méthode de contrôle en boucle fermée qui permette, à partir de la surveillance en cours de traitement d'une variable témoin de l'efficacité du traitement (comme le photoblanchiment), d'ajuster les modalités d'application de l'irradiation pour que ce paramètre objectif suive une évolution temporelle prédéterminée. Ainsi, le traitement peut être adapté en temps réel de telle sorte que le signal de photoblanchiment suive une évolution prédéterminée en étant peu ou pas affecté par les variabilités biologiques entre les sujets, de telle sorte à améliorer la reproductibilité du traitement. La détermination de l'évolution temporelle ou de la courbe souhaitée dépend de considérations qui sortent du champ de la présente invention. -15- Comme expliqué précédemment, le photoblanchiment est estimé en mesurant une intensité de fluorescence du PS. Cette mesure ne peut pas être effectuée en même temps que l'illumination car le signal de fluorescence ne serait pas discernable. Ainsi, dans le cadre de l'invention, l'illumination est appliquée sous la forme d'impulsions lumineuses telles que présentées à la Fig. 2(b), de telle sorte à pouvoir effectuer des mesures de photoblanchiment en temps réel entre ces impulsions. Dosimétrie Le fonctionnement du dispositif et du procédé selon l'invention est illustré à la Fig. 2. Le dispositif selon l'invention comprend une source de lumière d'activation 2 apte à générer une lumière d'activation 10. Cette lumière d'activation 10 est dirigée par exemple au moyen de fibres optiques vers les tissus biologiques 3 de la zone à traiter, dans lesquels un agent photosensibilisant PS a été préalablement introduit. Comme expliqué précédemment, ce PS excité par la lumière d'activation provoque la génération d'espèces phototoxiques qui entraînent la destruction des cellules. La source de lumière d'activation 2 est pilotée par un signal de commande d'activation 9 de forme temporelle u(t), issu de la boucle d'asservissement 1. Ce signal de commande d'activation 9 se présente sous la forme d'une séquence d'impulsions telle qu'illustrée à la fig. 1(b). Ainsi, la lumière d'activation issue de la source 2 pilotée par le signal d'activation 9 est appliquée sous la forme d'une suite temporelle L(t) d'impulsions lumineuses séparées par des moments d'extinction. La durée (et/ou l'amplitude) des impulsions lumineuses est contrôlée par le signal de commande 9. Le dispositif selon l'invention comprend également un spectrofluorimètre 4 pour effectuer des mesures de photoblanchiment du PS. Ce spectrofluorimètre 4 comprend une source de lumière d'excitation qui émet une lumière à des longueurs d'onde aptes à exciter la fluorescence du PS. La source lumineuse peut par exemple être une diode laser à 410 nm. Cette lumière d'excitation est dirigée vers les tissus biologiques 3 de la zone à traiter au moyen d'au moins une fibre optique d'illumination. La fluorescence 11 générée par le PS est collectée par des fibres optiques regroupées en faisceau autour de la ou des fibre(s) optique(s) d'illumination, pour être -16- amenée au spectromètre du spectrofluorimètre 4. Un spectre de fluorescence de forme /(t,.1), où À est la longueur d'onde d'émission de fluorescence, est collecté pour plusieurs instants t avec une période d'échantillonnage définie par l'utilisateur. L'intensité de fluorescence y(t) émise par l'agent photosensibilisant peut alors être déterminée par la formule suivante : y(t) = , .12) - .11) (Eq. 3) où J(t,À) est une intégrale du spectre : J(t, .1) = fco I (t, 1) dl , (Eq. 4) et l'intervalle [.11;.12] correspond à la largeur de bande d'émission de 10 fluorescence de l'agent photosensibilisant. Cette intensité de fluorescence y(t) est ensuite utilisée comme signal de mesure du photoblanchiment 12. Le dispositif selon l'invention peut être mis en oeuvre pour traiter des affections à la surface de la peau ou d'une paroi biologique (oesophage, ...). Il peut également être mise en oeuvre pour traiter des affections ou des 15 tumeurs interstitielles. Dans ce cas, la lumière d'activation 10 peut être amenée dans la zone à traiter 3 et la fluorescence 11 peut être collectée avec des fibres optiques insérées dans les tissus biologiques. Comme expliqué précédemment, la mesure du photoblanchiment 12 ne peut être effectuée que lorsque la source d'activation 2 est éteinte, car dans 20 le cas contraire la fluorescence 11 n'est pas discernable. Le dispositif selon l'invention comprend enfin une boucle d'asservissement 1 implémentée dans un calculateur numérique. La fonction de cette boucle d'asservissement 1 est d'ajuster le signal de commande d'activation 9 de forme temporelle u(t) de telle sorte que le signal de 25 photoblanchiment 12 de forme temporelle y(t) suive une consigne 5 de forme temporelle R(t) prédéterminée. Bien entendu, dans le cadre de l'invention cette consigne peut être de toute forme temporelle R(t) . Dans la boucle d'asservissement 1, une grandeur d'état 15 de forme 30 temporelle x(t) qui dépend du signal de photoblanchiment 12 est comparée à la consigne 5 dans un comparateur 6, de telle sorte à générer un signal d'erreur de poursuite 7 de forme temporelle e(t). Le signal d'erreur de poursuite 7 est introduit dans un correcteur 8 qui génère un signal de commande 9 de la source d'activation 2 visant à -17- minimiser cette erreur aux itérations suivantes. Ce signal de commande d'activation 9 a une forme temporelle u(t) qui correspond à une suite d'impulsions avec une période constante mais dont l'amplitude et/ou le rapport cyclique est ajusté en fonction du signal d'erreur 5 de poursuite 7. Dans le cas d'un ajustement du rapport cyclique, qui est la solution mise en oeuvre dans le mode de réalisation présenté, ces impulsions sont obtenues à partir d'un signal sinusoïdal auquel est ajouté un offset dépendant de l'amplitude de l'erreur de poursuite e (t) . Ce signal sinusoïdal décalé en amplitude est ensuite transformé en impulsions de période 10 constante et de rapport cyclique dépendant de l'offset calculé au moyen d'un comparateur, pour générer le signal de commande u(t). Suivant un premier mode de réalisation de l'invention, la boucle d'asservissement 1 comprend un interpolateur 13 apte à générer la grandeur d'état 15 à partir d'une interpolation, par une méthode connue de l'homme du 15 métier, des mesures de photoblanchiment 12. Cette interpolation consiste à calculer des valeurs de la grandeur d'état x(t) à partir d'un historique de mesures de photoblanchiment y (t) , éventuellement avec une fréquence d'échantillonnage plus élevée, et en effectuant un moyennage du bruit de mesure. Des interpolations de type polynômiales ou splines peuvent 20 notamment être utilisées. Suivant un second mode de réalisation préférentiel de l'invention, la boucle d'asservissement comprend un observateur d'état 14 dont le fonctionnement sera expliqué ci-après. Bien entendu, l'interpolateur 13 et l'observateur d'état 14 peuvent être 25 tous deux implémentés dans une boucle d'asservissement 1. Le choix de l'utilisation de l'un ou de l'autre peut alors être déterminé par un paramètre accessible à l'utilisateur. En pratique, les mesures de fluorescence dans des conditions in-vivo sont des mesures difficiles. L'intensité de fluorescence 11 est très faible. Ces 30 mesures nécessitent des temps d'intégration de l'ordre de la seconde, et demeurent très bruitées. En outre, les périodes d'illumination et d'obscurité du signal lumineux d'activation 10 sont de l'ordre de 1s à quelques dizaines de secondes. Sachant que les mesures de fluorescence 11 ne peuvent être effectuées que pendant les périodes d'obscurité du signal lumineux 35 d'activation 10, il s'ensuit que le signal de photoblanchiment 12 qui peut être -18- obtenu est en général très bruité et comprend peu d'échantillons. Dans ces conditions, un interpolateur classique 13 risque de donner de mauvais résultats, et de générer une grandeur d'état 15 également très bruitée qui conduit à un mauvais fonctionnement global de la boucle de régulation 1. L'estimation de la grandeur d'état 15 à partir des mesures de photoblanchiment 12 peut être considérablement améliorée en mettant en oeuvre un observateur d'état 14. Il s'agit d'un capteur logiciel qui est basé sur un modèle dynamique 16 du phénomène de photoblanchiment. Il accepte en 10 entrée le signal de commande d'activation 9 et les mesures de photoblanchiment 12, et fournit en sortie une estimation de la grandeur d'état 15. Le modèle dynamique 16 a pour fonction de simuler la chaine d'activation du PS et de mesure du photoblanchiment, qui comprend la 15 génération du signal lumineux d'activation 10 par la source 2, les interactions avec les tissus biologiques 3 et la mesure du photoblanchiment 12 avec le spectrofluorimètre 4. Le modèle biologique de l'Eq. 2 est difficilement applicable, car il dépend d'un grand nombre de paramètres difficiles à déterminer, et variables entre 20 les sujets et les types de tissus biologiques. Modèle dynamique simplifié Il a donc été développé dans le cadre de l'invention un modèle dynamique simplifié, obtenu directement à partir de données expérimentales, et non pas à partir des équations théoriques des réactions photochimiques 25 comme le modèle de l'Eq. 2. Ce modèle dynamique simplifié permet à l'observateur d'état 14 de reconstruire en temps réel un signal de photoblanchiment "propre" à partir des données partielles et bruitées 12 mesurées pendant le traitement. Dans le cas de la thérapie photodynamique, obtenir un modèle simplifié 30 du comportement dynamique du phénomène de photoblanchiment est complexe. Il y a deux raisons à cela: - le comportement non linéaire du phénomène de photoblanchiment visible sur les simulations de la Fig. 1(a) et 1(b). En effet, on observe clairement un comportement asymétrique de la réponse lorsque le signal 35 d'excitation lumineuse est activé ou pas. Or, ce changement d'allure ne peut -19- être décrit par un modèle linéaire classique; - de plus, comme expliqué précédemment, on ne peut mesurer les spectres de fluorescence de l'agent photosensibilisant que pendant les phases où la lumière d'activation 10 est éteinte sinon elle sature le signal en raison 5 de sa forte puissance. On est donc "aveugle" une grande partie du temps. Ce modèle dynamique a été développé sur la base d'une structure ARX (« Auto Regressive Model with eXternal inputs » en Anglais ou modèle autorégressif avec des entrées externes). Cette structure ARX est définie par l'équation aux différences : y(t) + a1y(t - 1) + --- + anay(t - na) 10 = biu(t - nk) + --- + bnbu(t - nk - nb + 1) + v(t). (Eq. 5) y(t) désigne la variable de sortie qui est ici le signal de photoblanchiment 12). u(t) désigne la variable d'entrée qui est ici le signal de commande 9 de la source lumineuse d'activation 2. L'erreur v(t) est une variable aléatoire dont la distribution suit une loi normale de moyenne nulle et de variance a2 , 15 N(0,0-2). t est la variable du temps qui représente les instants discrets échantillonnés avec une période d'échantillonnage T. na,nb,nk sont des indices de structure, comprenant respectivement les ordres des polynômes des parties auto-régressives et externes ainsi que le retard de la sortie par rapport à l'entrée. 20 Pour décrire le comportement asymétrique du phénomène de photoblanchiment lorsque le signal d'excitation lumineuse est activé ou pas, on introduit un modèle ARX par morceaux, défini par les équations suivantes, y(t) + ctiay(t - 1) + ... + aimuy(t - na) = biau(t - nk) + --- + bi,nbu(t - nk - nb + 1) + v1(t) (Eq- 6) si 8u(t) 0, et y(t) + a2,07(t - 1) + ... + a2 ,,,y(t - na) 25 = b2,1u(t - nk) + --- + b2,,bu(t - nk - nb + 1) + v2(t) (Eq. 7) si 8u(t) 0. 8u(t) = u(t) - u(t - 1) correspond aux variations du signal de commande 9. v1(t) et v2(t) sont deux processus Gaussiens centrés indépendants et identiquement distribués. 30 Dans le cas présent, on peut se limiter à na = nb = nk = 1. Les Eq. 6 et Eq. 7 peuvent être regroupées au sein d'une même relation, -20- y(t) + ctiay(t - 1) + (a2,1 - ctia)y(t - 1)b (t) = biau(t - 1) + (b2,1 - 1)u(t - 1)b (t) + v1(t) + (v2 (t) - v1(t))b(t). (Eq. 8) La variable de séquencement b(t) est une variable binaire nécessaire à la sélection des deux modèles linéaires. Elle est définie de la manière suivante : {b(t) = 0 si su(t) > 0 b(t) = 1 si 8u(t) Finalement, pour décrire complètement le modèle, il reste les paramètres a1,1 et b1,1 à déterminer. Ces paramètres peuvent être estimés à partir de mesures expérimentales réalisées dans des conditions proches des conditions d'application du traitement, par exemple lors d'expérimentations in-vivo sur des animaux. La détermination de ces paramètres peut être effectuée avec le montage de la Fig. 2, dans une phase d'identification effectuée en boucle ouverte, en exploitant des mesures de photoblanchiment y(t) obtenues en réponse à un signal de commande u(t) connu. Les paramètres a1,1 et b1,1 sont ajustés pour que les Eq. 10 et Eq. 11 décrivent au mieux les mesures. La Fig. 4 illustre la qualité d'ajustement de ce type de modèle sur des données de photoblanchiment obtenues en condition in-vitro. La courbe 35 représente des mesures de photoblanchiment y(t) obtenues en réponse à un signal de commande u(t), et la courbe 36 représente la réponse du modèle des Eq. 10 et Eq. 11 dans ces conditions. La correspondance entre la réponse du modèle et les mesures (coefficient de détermination) est de l'ordre de 95%. Les Eq. 10 et Eq. 11 constituent ainsi le modèle dynamique 16 du -21- phénomène de photoblanchiment, donc on peut voir qu'il peut être obtenu de manière simple en faisant un minimum d'hypothèses sur des paramètres physiologiques. On peut également en déduire directement un modèle pour l'observateur 5 d'état 14 comme suit : x(t + 1) = -ctiax(t) B u(t) +KE(t) (Eq. 12) 3î(t) = x(t) (Eq. 13) E(t) = y(t) - y(t), (Eq. 14) Avec 3-(0) = y(0). y(t) est le photoblanchiment estimé par l'observateur 10 d'état 14, et K un facteur multiplicatif (gain d'innovation). Il est à noter que ce modèle est également en mesure de prendre en compte le photoblanchiment du PS dû à la lumière d'excitation émise par le spectrocolorimètre 4 pour générer la fluorescence. Le déroulement des étapes du procédé d'application de la PDT selon 15 l'invention est illustré à la figure 3 : - étape 20 : dans un premier temps, une trajectoire de photoblanchiment R(t) prédéterminée est programmée ; - étape 21 : une mesure de photoblanchiment est effectuée au moyen du spectrofluorimètre 4 pour obtenir une valeur de photoblanchiment mesurée 20 Y(t) ; - étape 22 : la variable d'état x(t) (qui correspond à une valeur estimée de photoblanchiment 3-(t)) est calculée par l'observateur d'état 14, en prenant en compte la dernière valeur du signal de commande de la source lumineuse d'activation u(t) ; 25 - étape 23 : l'erreur de poursuite e(t) est calculée ; - étape 24 : le signal de commande de la source lumineuse d'activation u(t) est calculé par le correcteur 8 ; - étape 25 : le signal lumineux d'activation L(t) est appliqué par la source d'activation 2 à la zone à traiter 3 sous la forme d'une impulsion 30 l'illumination dont la durée est déterminée par le signal de commande u(t) ; - étape 26 : le système attend la fin de l'impulsion lumineuse L(t) ; - étape 27 : si la trajectoire de photoblanchiment R(t) est complètement parcourue, le traitement est arrêté (étape 28). Dans le cas contraire l'exécution de l'algorithme est poursuivie de manière itérative à l'étape 21 par -22- une nouvelle mesure de photoblanchiment (retour à l'étape 21). Il est à noter que cet algorithme comprend deux boucles imbriquées : - une première boucle de mesure 29 entre les étapes 27 et 21 qui est cadencée par les mesures de photoblanchiment ; - une deuxième boucle de mesure 30 entre les étapes 24 et 22 qui permet de prendre en compte les variations du signal de commande u(t) dans les estimations du photoblanchiment ÿ(t)) fournies par l'observateur d'état 14. La boucle de contre réaction 30 peut ainsi être exécutée à une cadence très supérieure à la boucle de mesure 29 (par exemple de l'ordre de 10 fois) 10 en se basant sur des estimations produites par l'observateur à cadence régulière, même lorsque les mesures de photoblanchiment sont impossibles (durant l'illumination). Ainsi, le bruit dans la boucle de contre réaction peut être fortement diminué. La période de la boucle de mesure 29 détermine la période du signal de 15 commande u(t). La durée des impulsions de ce signal de commande u(t) ou son rapport cyclique est ajustée par la boucle de contre réaction 30. Enfin, l'observateur d'état 14 peut être « recalé » à chaque itération de la boucle de mesure, ce qui permet de relâcher considérablement les contraintes sur la précision nécessaire pour le modèle dynamique 16. 20 La Fig. 5 illustre les performances de l'observateur de photoblanchiment. La courbe 40 est une courbe de photoblanchiment théorique calculée avec le modèle cytotoxique de l'Eq. 2. Cette courbe simule, de manière idéale et sans bruit de mesure, un processus d'asservissement du photoblanchiment effectué selon l'invention, qui suit une trajectoire de photoblanchiment R(t) 25 qui est une droite. A partir de cette courbe, des mesures bruitées 41 ont été générées en y rajoutant du bruit blanc, pour simuler des mesures de photoblanchiment y(t) telles qu'elles auraient été obtenues avec le spectrofluorimètre 4 dans des conditions réelles. La courbe 42 présente les valeurs estimées de photoblanchiment 3-(t) produites par l'observateur d'état 30 14 à partir de la commande u(t) et des mesures bruitées y(t) dans ces conditions. On observe bien la capacité de l'observateur d'état 14 à estimer précisément la courbe du photoblanchiment théorique 40 malgré le bruit de mesure. La Fig. 6 présente des mesures de l'asservissement du photoblanchiment 35 obtenu avec le procédé et le dispositif selon l'invention, dans des conditions -23- in-vitro. La courbe 52 correspond à la trajectoire de photoblanchiment souhaitée R(t), la courbe 50 correspond aux mesures de photoblanchiment y (t) , et la courbe 51 correspond au signal de commande de l'activation u(t). On peut observer que malgré un écart important au départ entre la consigne 5 R(t) et les mesures y (t) , le système est capable de faire converger rapidement puis d'asservir le photoblanchiment au plus près de la trajectoire de consigne. Suivant des variantes de modes de réalisation, la source de lumière d'activation 2 et le correcteur 8 peuvent être conçus de telle sorte à générer : - des impulsions de lumière 10 d'amplitude variable en fonction de 10 l'erreur de poursuite 7. Ces impulsions peuvent être de rapport cyclique constant (par exemple de 50%) et de période constante ; - des impulsions de lumière 10 de durée et de période variable en fonction de l'erreur de poursuite 7. Bien sûr, l'invention n'est pas limitée aux exemples qui viennent d'être 15 décrits et de nombreux aménagements peuvent être apportés à ces exemples sans sortir du cadre de l'invention. REVENDICATIONS1. Procédé pour contrôler un rayonnement électromagnétique d'activation (10) d'une substance photosensibilisante, caractérisé en ce qu'il comprend des étapes : - d'obtention d'une information de luminescence (11, 12, 15) issue de ladite substance photosensibilisante, et - d'ajustement de l'exposition de la substance photosensibilisante au rayonnement électromagnétique d'activation (10) de telle sorte que ladite information de luminescence (11, 12, 15) suive une évolution temporelle correspondant à une consigne temporelle prédéterminée (5). 2. Le procédé de la revendication 1, dans lequel le rayonnement électromagnétique d'activation (10) comprend un rayonnement lumineux d'activation avec des longueurs d'onde optiques. 3. Le procédé de l'une des revendications 1 ou 2, qui comprend en outre une étape de mesure d'une information de luminescence issue de l'oxygène singulet généré par l'activation de la substance photosensibilisante. 4. Le procédé de la revendication 2, qui comprend en outre des étapes : - d'exposition de la substance photosensibilisante à un rayonnement lumineux d'activation (10) appliqué sous forme d'impulsions lumineuses d'activation, et - de mesure de luminescence (11) entre lesdites impulsions. 5. Le procédé de la revendication 4, dans lequel la mesure de luminescence (11) comprend des étapes : - d'illumination de la substance photosensibilisante avec une lumière d'excitation à une longueur d'onde d'excitation, - de mesure d'une intensité de fluorescence (11). 6. Le procédé de la revendication 5, qui comprend en outre une étape de calcul d'un observateur de luminescence (14) utilisant en entrée au moins une-25- mesure de luminescence (12) et une commande des impulsions lumineuses d'activation (9), et produisant en sortie une estimation de la luminescence (15). 7. Le procédé de la revendication 6, dans lequel l'observateur de luminescence (14) est modélisé sous la forme d'un observateur d'état (14) comprenant une variable d'état (15) correspondant à une estimation de la luminescence. 8. Le procédé de la revendication 7, dans lequel l'observateur d'état (14) est dérivé d'un modèle autorégressif (16) avec entrées externes ARX, ou d'un modèle de type représentation d'état. 9. Le procédé de l'une des revendications 6 à 8, dans lequel l'observateur de luminescence (14) comprend des paramètres qui sont déterminés lors d'une étape préalable de calibration à partir de mesures expérimentales. 10. Le procédé de l'une des revendications 6 à 9, qui comprend en outre une étape de comparaison (6) de l'estimation de la luminescence (15) issue de l'observateur de luminescence (14) avec la consigne temporelle (5) de telle sorte à générer un signal d'erreur (7), lequel étant utilisé pour l'ajustement de l'exposition au rayonnement électromagnétique d'activation (10). 11. Le procédé de l'une des revendications 4 à 10, dans lequel l'ajustement de l'exposition au rayonnement électromagnétique d'activation (10) comprend un ajustement de l'amplitude et/ou de la durée des impulsions lumineuses. 12. Le procédé de l'une des revendications précédentes, qui est mis en oeuvre avec une substance photosensibilisante dont l'activation produit des substances aptes à détruire des cellules. 13. Dispositif pour contrôler l'activation d'une substance-26- photosensibilisante préalablement introduite dans un tissu biologique (3), comprenant des moyens d'exposition (2) aptes à permettre une exposition de ladite substance photosensibilisante à un rayonnement électromagnétique d'activation (10) de telle sorte à l'activer, caractérisé en ce qu'il comprend en outre : - des moyens d'acquisition (4) aptes à produire une information de luminescence (11, 12) issue de ladite substance photosensibilisante, et - des moyens électroniques et de calcul (1) aptes à contrôler les moyens d'exposition (2) de telle sorte à permettre un ajustement de l'exposition de la substance photosensibilisante au rayonnement électromagnétique d'activation (10) de telle sorte que ladite information de luminescence (11, 12) suive une évolution temporelle correspondant à une consigne temporelle prédéterminée (5). 14. Le dispositif de la revendication 13, dans lequel les moyens d'exposition comprennent une source lumineuse d'activation (2) apte à produire des impulsions de lumière (10) en fonction d'un signal de commande (9). 15. Le dispositif de l'une des revendications 13 ou 14, dans lequel les moyens d'acquisition (4) comprennent un spectromètre ou un spectrofluorimètre. 16. Le dispositif de l'une des revendications 14 ou 15, dans lequel les moyens électroniques et de calcul (1) implémentent une boucle de contrôle comprenant un observateur de luminescence (14), - lequel observateur de luminescence utilisant en entrée au moins une mesure de luminescence (12) et un signal de commande (9) de la source lumineuse d'activation (2), et produisant en sortie une estimation de la luminescence (15), - laquelle estimation de la luminescence (15) étant utilisée pour générer un signal de commande (9) de la source lumineuse d'activation (2).