FR 2998056 A1 20140516 FR 1260892 A 20121115 PROCÉDÉ D'ESTIMATION AUTOMATIQUE D'UNE VALEUR D'UN TAUX DE SURVIE DE CELLULES BIOLOGIQUES IRRADIÉES PAR UN RAYONNEMENT IONISANT [on] L'invention concerne un procédé d'estimation automatique d'une valeur d'un taux de survie de cellules biologiques irradiées par un rayonnement ionisant. L'invention concerne également un procédé de calcul d'une dose d'irradiation pour irradier ces cellules biologiques. L'invention concerne en outre un support d'enregistrement d'informations contenant des instructions pour l'exécution de ces procédés. L'invention concerne enfin un calculateur électronique programmé pour mettre en oeuvre ce procédé d'estimation. [2] On connaît des procédés pour estimer un taux de survie de cellules biologiques irradiées par un rayonnement ionisant. Ces procédés sont typiquement utilisés dans des systèmes de planification de traitement (« treatment planning system » ou « radiation treatment planning » en langue anglaise) en radiothérapie et, plus particulièrement, en hadronthérapie. Ces systèmes permettent notamment de calculer des paramètres physiques d'un rayonnement ionisant, en vue d'irradier des cellules biologiques, telles que des tumeurs. Un exemple d'un tel procédé est décrit dans l'article, noté Al, de M. Krâmer et M. Scholz, « Rapid calculation of biological effects in ion rad iotherapy », Physics in Medecine and Biology, vol. 51, p. 1959-1970, 2006. [3] Ces procédés comportent typiquement l'application d'un coefficient de correction CrPelf.P2 sur le taux de survie estimé, pour tenir compte de la concentration en oxygène dans les cellules (typiquement mesurée par la pression partielle d'oxygène de ces cellules). En effet, des cellules biologiques possédant une concentration réduite en oxygène peuvent présenter une résistance accrue à des rayonnements ionisants. Cette résistance accrue, connue sous le nom d « effet oxygène », est décrite, par exemple, dans l'article de L. H. Gray et al., « The concentration of oxygen dissolved in tissues at the time of irradiation as a factor in radiotherapy », British Journal of Radiology, vol. 26, p. 638-648,1953. Cependant, ce coefficient de correction est typiquement fourni uniquement pour un rayonnement ionisant dit de référence, et ne dépend pas à la fois de la concentration en oxygène de ces cellules et du type et de l'énergie du rayonnement ionisant. Ce coefficient de correction doit ainsi être à nouveau déterminé si des paramètres expérimentaux (tels que l'énergie, ou le type de particules du rayonnement ionisant) sont modifiés. Une telle détermination comprend généralement des étapes expérimentales, longues et coûteuses à mettre en oeuvre. En outre, une estimation erronée de l'amplitude de l'effet oxygène sur ce taux de survie peut compromettre l'efficacité du traitement. [4] II existe donc un besoin pour un procédé d'estimation d'un taux de survie de cellules biologiques irradiées par un rayonnement ionisant, qui prenne en compte de manière fiable la concentration réelle en oxygène de ces cellules biologiques pour le calcul du taux de survie, ce procédé pouvant être mis en oeuvre de manière simplifiée. [5] L'invention concerne donc un procédé d'estimation automatique d'une valeur STP,ECD d'un taux de survie de cellules biologiques irradiées par une dose D d'un premier rayonnement ionisant comportant des particules de type T et d'énergie E, ces cellules présentant une pression partielle P d'oxygène, caractérisé en ce que ce 10 procédé comporte : -l'acquisition d'un taux de survie SE des cellules biologiques lorsqu'elles sont irradiées par une dose D du premier rayonnement ionisant et qu'elles présentent une pression partielle P1 d'oxygène, chaque valeur STP,1E,CD de ce taux de survie STP,iE étant fonction de la valeur de la dose D; 15 -l'acquisition d'un coefficient de correction CrPelfP2 , ce coefficient étant une fonction de la valeur de la dose D, chaque valeur CrPelfP2CD de ce coefficient CrPelf.P2 pour une valeur particulière de la dose D exprimant la dépendance du taux de survie des cellules à la pression partielle d'oxygène dans les cellules, lorsque ces cellules sont irradiées par cette dose D d'un second rayonnement ionisant, ce second 20 rayonnement ionisant comportant des particules d'un type différent du premier rayonnement ionisant, chaque valeur CrPelfP2CD permettant de convertir : - une valeur SrPelf(D d'un taux de survie SrPelf des cellules biologiques, irradiées par la dose D du second rayonnement ionisant, lorsqu'elles présentent une pression partielle P1; 25 - en une valeur Sf (D d'un taux de survie Sf de ces mêmes cellules biologiques, irradiées par la même dose D du second rayonnement ionisant lorsqu'elles présentent une pression partielle d'oxygène P2, distincte de la pression partielle P1; -la détermination d'une quantité Qi représentative d'une quantité d'une espèce 30 chimique donnée R, produite par radiolyse du milieu intracellulaire des cellules biologiques sous l'effet de la dose D du premier rayonnement ionisant, lorsque ces cellules présentent la pression partielle P d'oxygène ; -la détermination d'une dose D du second rayonnement ionisant qu'il faudrait appliquer sur les mêmes cellules biologiques pour obtenir, par radiolyse dudit milieu 35 cellulaire, la même quantité Qi représentative de la quantité de l'espèce chimique R, lorsque ces cellules présentent la pression partielle P2 d'oxygène ; -le calcul d'une valeur CrPeliP2(D7) du coefficient de correction CrPelf.P2 pour la dose D' ; -le calcul de la valeur STP,ECD pour la dose D, selon la relation : STPE(D =CP'iiP2(13' xSPTiE,(13 [6] Dans ce procédé, pour obtenir la valeur sTP,E(D), au lieu d'utiliser directement la valeur cTP,1,'P (D) d'un coefficient de correction pour corriger la valeur sTP,1E(D) acquise du taux de survie sTP,I, , on utilise le coefficient de correction CrPelfP2 défini pour le second rayonnement. Ainsi, ce procédé permet de calculer la valeur du coefficient de correction cTP,1,'P (D) pour le premier rayonnement, en utilisant des données déjà disponibles pour le second rayonnement. Cela permet de réduire le nombre de mesures expérimentales nécessaires pour déterminer la valeur cTP,ii(D) du coefficient de correction pour un nouveau rayonnement, car on peut utiliser des données existantes. La mise en oeuvre du procédé d'estimation du taux de survie des cellules biologiques est donc simplifiée. [7] Les modes de réalisation de l'invention peuvent comporter une ou plusieurs des caractéristiques suivantes : -l'espèce chimique R est un radical libre produit par une ou une suite de réactions chimiques dont un réactif est une molécule d'oxygène et un autre réactif est un électron aqueux eaq ou un radical hydrogène H', cet électron aqueux eaq ou ce radical hydrogène H étant un produit de la radiolyse du milieu intracellulaire des cellules par la dose D du premier rayonnement ; -l'élément R est compris dans le groupe composé de l'anion superoxyde 02-, du radical hydroperoxyde H02* et des radicaux G02* et GS02' issus, respectivement, de réactions de l'oxygène avec des radicaux G' et GS*, ces radicaux G' et GS* étant issus eux-mêmes des réactions du glutathion GSH avec l'électron aqueux ou avec le radical hydrogène ; -le coefficient de correction CP1P2 est défini, pour toute dose D, par la formule : ref C rPeliP 2CD =SPre2f(D S rPeif (D ; -le second rayonnement présente un transfert d'énergie linéique inférieur ou égal à 1 keV/pm sur la majorité de la longueur de sa trajectoire dans le milieu contenant les cellules biologiques ; -le premier rayonnement présente un transfert d'énergie linéique supérieur ou égal à 10keV/pm au niveau des cellules ; -les particules de type T sont des hadrons ; -les particules de type T sont des ions carbone 1-2C6+. posj Ces caractéristiques présentent en outre les avantages suivants : -le choix de l'espèce chimique R comme étant l'une de celles issues d'une ou une suite de réactions chimiques dont un réactif est une molécule d'oxygène et un autre réactif est un électron aqueux eaq ou un radical hydrogène H' permet d'accroître la précision de l'évaluation du taux de survie ; -l'utilisation de hadrons, et plus particulièrement d'ions carbone, permet d'irradier efficacement des cellules biologiques tout en réduisant l'endommagement du tissu biologique voisin de ces cellules biologiques. [009] L'invention concerne également un procédé de calcul automatique d'une valeur d'une dose D d'irradiation de cellules biologiques avec un rayonnement ionisant, pour que ces cellules présentent, après irradiation, une valeur s prédéfinie d'un taux de survie, ce procédé comportant : a] l'acquisition de la valeur s prédéfinie et d'une cartographie des valeurs de la pression partielle P d'oxygène des cellules biologiques à irradier ; b] le calcul de paramètres du rayonnement ionisant nécessaires pour que les cellules biologiques, une fois irradiées avec ce rayonnement ionisant, présentent la valeur s prédéfinie de taux de survie, ce rayonnement ionisant comprenant, au niveau des cellules biologiques , des particules de type T et d'énergie E et ce calcul comprenant notamment la détermination de la valeur de la dose D de ce rayonnement à appliquer sur les cellules biologiques ; c] l'estimation d'une valeur STP,EiD d'un taux de survie STP,E, des cellules biologiques irradiées par la dose D du rayonnement ionisant déterminée lors de l'étape b] avant application de cette dose D du rayonnement sur les cellules biologiques : d] la comparaison de la valeur STP,EiD estimée lors de l'étape c] avec la valeur S 20 prédéfinie, puis, - la fourniture de la valeur de la dose D calculée lors de l'étape b] si la valeur estimée sTP,EiD correspond à la valeur S prédéfinie, et, - dans le cas contraire, la réitération de l'étape b] pour déterminer une nouvelle valeur de la dose D; 25 caractérisé en ce que la valeur SÇ(D du taux de survie STP,E est estimé selon un procédé conforme à l'invention. [0olo] Selon un autre aspect, l'invention concerne un support d'enregistrement d'informations, caractérisé en ce qu'il comprend des instructions pour l'exécution d'un procédé conforme à l'une quelconque des revendications précédentes lorsque ces 30 instructions sont exécutées par un calculateur électronique. [0ou] Selon un autre aspect, l'invention concerne un calculateur électronique pour la mise en oeuvre d'un procédé d'estimation conforme à l'invention, ce calculateur étant programmé pour : -acquérir d'un taux de survie STP,1E, des cellules biologiques lorsqu'elles sont irradiées 35 par une dose D du premier rayonnement ionisant et qu'elles présentent une pression partielle P1 d'oxygène, chaque valeur STP,1ElD de ce taux de survie STP,1E, étant fonction de la valeur de la dose D ; -acquérir un coefficient de correction CrPelfP2 , ce coefficient étant une fonction de la valeur de la dose D, chaque valeur CrPelfP2CD de ce coefficient CrPelfP2 pour une valeur particulière de la dose D exprimant la dépendance du taux de survie des cellules à la pression partielle d'oxygène dans les cellules, lorsque ces cellules sont irradiées par cette dose D d'un second rayonnement ionisant, ce second rayonnement ionisant comportant des particules d'un type différent du premier rayonnement ionisant, chaque valeur CrPelfP2CD permettant de convertir : - une valeur SrPelf(D d'un taux de survie SrPelf. des cellules biologiques, irradiées par la dose D du second rayonnement ionisant, lorsqu'elles présentent une pression partielle P1; - en une valeur SrPe2f. (D d'un taux de survie SrPe2f. de ces mêmes cellules biologiques, irradiées par la même dose D du second rayonnement ionisant lorsqu'elles présentent une pression partielle d'oxygène P2, distincte de la pression partielle P1; -déterminer une quantité Qi représentative d'une quantité d'une espèce chimique donnée R, produite par radiolyse du milieu intracellulaire des cellules biologiques sous l'effet de la dose D du premier rayonnement ionisant, lorsque ces cellules présentent la pression partielle P d'oxygène ; -déterminer une dose D du second rayonnement ionisant qu'il faudrait appliquer sur 20 les mêmes cellules biologiques pour obtenir, par radiolyse dudit milieu intracellulaire, la même quantité Qi représentative de la quantité de l'espèce chimique R, lorsque ces cellules présentent la pression partielle P2 d'oxygène ; -calculer une valeur CrPeliP2CD' du coefficient de correction CrPelfP2 pour la dose D' ; -calculer la valeur sTP,ECD pour la dose D, selon la relation : 25 STPECD =Cl;P2(13' xSTP1E,CD . [0012] L'invention sera mieux comprise à la lecture de la description qui va suivre, donnée uniquement à titre d'exemple non limitatif et faite en se référant aux dessins sur lesquels : -la figure 1 est une illustration schématique d'un dispositif pour irradier des cellules 30 biologiques avec un rayonnement ionisant ; -la figure 2 illustre deux taux de survie de cellules biologiques de la figure 1 irradiées par un rayonnement ionisant, ces taux de survie correspondant chacun à une valeur différente de la pression partielle d'oxygène présente dans les cellules biologiques de la figure 1 ; 35 -la figure 3 illustre une courbe d'un rendement de production d'une espèce chimique par radiolyse d'un milieu aqueux représentatif du milieu intracellulaire des cellules biologiques de la figure 1, en fonction de la pression partielle d'oxygène contenue dans ces cellules ; -la figure 4 illustre une courbe du « rapport d'amplification par oxygène » des cellules biologiques de la figure 1 en fonction de la pression partielle en oxygène que présentent ces cellules biologiques ; -la figure 5 est une illustration schématique d'un dispositif de calcul pour calculer une 5 dose d'irradiation à fournir aux cellules biologiques de la figure 1; -la figure 6 est un organigramme d'un procédé d'utilisation du dispositif de la figure 5 ; -la figure 7 est un organigramme détaillant une étape du procédé de la figure 6. [0013] Dans ces figures, les mêmes références sont utilisées pour désigner les mêmes éléments. 10 [0014] Dans la suite de cette description, les caractéristiques et fonctions bien connues de l'homme du métier ne sont pas décrites en détail. ponj La figure 1 représente un exemple d'irradiation de cellules biologiques par un rayonnement ionisant 2. Les cellules sont ici comprises dans un tissu biologique 4. Pour simplifier la figure, les cellules biologiques 6 sont représentées de façon 15 schématique sous la forme d'un amas de cellules. Le rayonnement 2 est ici émis par une source de rayonnement 8. [0016] Le rayonnement 2 émis par la source 8 est un faisceau de particules d'un type T et d'une énergie Eo. Pour simplifier cette description, on considère, dans un premier temps, que les cellules 6 reçoivent uniquement des particules de type T et 20 dont l'énergie E, est identique à Eo. Dans cet exemple, les particules T du rayonnement 2 sont des hadrons. Plus précisément, ces particules T sont des ions carbone 1206+. Ce rayonnement 2 présente ici, à la sortie de la source 8, une énergie Eo supérieure ou égale à 1 MeV par nucléon ou à 10 MeV par nucléon. [0017] Le rayonnement 2 présente ici un transfert d'énergie linéique (TEL) (« linear 25 energy transfer » en langue anglaise) moyen, en entrée dans les cellules 6, supérieur ou égal à 1keV/pm, ou à 10keV/pm ou à 200keV/pm. Cette valeur de TEL est ici définie pour le pic de Bragg du rayonnement 2 dans la partie du tissu 4 contenant les cellules 6. ponj Dans cet exemple, les cellules 6 sont des cellules cancéreuses de la lignée 30 cellulaire SQ20B issue d'une tumeur de la région tête et cou d'un être humain. [0019] Les cellules 6 contiennent de l'oxygène (02). Dans cette description, la concentration en oxygène dans les cellules 6 est quantifiée au moyen de la pression partielle d'oxygène qui règne dans ces cellules 6. Ces valeurs de pressions partielles d'oxygène sont, par exemple, obtenues à partir d'un cliché d'imagerie à résonance 35 magnétique de la zone 4, ce cliché étant réalisé après avoir injecté un marqueur sensible au niveau d'oxygène dans cette zone 4. Ces cellules 6 présentent ici une pression partielle d'oxygène P uniforme. Dans le cas où différentes cellules appartenant à la population totale de cellules à irradier de la zone 4 présentent des valeurs de pressions partielles différentes, alors la concentration en oxygène peut se 40 présenter sous la forme d'une cartographie de ces différentes valeurs de pressions partielles. Dans cet exemple, pour simplifier, l'amas de cellules 6 est un sous-ensemble de cette population totale des cellules, présentant la même pression partielle d'oxygène P. [0020] Pour simplifier, dans la suite de cette description, sauf précision contraire, le terme de « pression P » désignera cette pression partielle d'oxygène P. [0021] La source 8 est configurée pour fournir une dose prédéfinie du rayonnement 2 sous forme d'irradiation en particules de type T et d'énergie E sur les cellules 6, en vue de détruire une partie de ces cellules 6, tout en limitant la destruction de portions du tissu 4 ne contenant pas de cellules 6. Cette source 8 est, par exemple, un accélérateur de particules. [0022] Dans cette description, en l'absence de précision contraire, le terme de « dose d'un rayonnement » fait référence à la dose absorbée de ce rayonnement, telle que définie par la Commission internationale des unités et mesures radiologiques (« International Commission on Radiation Units and Measurements » en langue anglaise). Cette dose est exprimée en Gray. [0023] Typiquement, l'efficacité de l'irradiation des cellules 6 est quantifiée à l'aide du taux de survie que présentent les cellules 6 après irradiation. Dans cette description, le taux de survie d'une population de cellules 6 est défini comme le pourcentage de ces cellules 6 ayant survécu après que ces cellules 6 aient été irradiées avec une dose D. [0024] Par exemple, on définit le taux de survie sTP,E pour les cellules 6 présentant une pression partielle d'oxygène P et irradiées par le rayonnement 2. Ce taux de survie sTP,E est une fonction qui associe, à toute valeur d'une dose D du rayonnement 2, la valeur sTP,EiD du taux de survie des cellules 6, présentant la pression P et recevant cette valeur de dose D du rayonnement 2. Ce taux sTP,E dépend également du type T et de l'énergie E du rayonnement appliqué. C'est pourquoi les lettres T et E sont notées en indice de cette fonction SE . [0025] Pour simplifier, dans la suite de cette description, la convention de notation suivante est adoptée : -on note par la lettre S la fonction taux de survie des cellules 6 ; -la pression P pour laquelle est calculé ce taux de survie est notée en exposant de la lettre S; -l'identifiant du rayonnement ionisant pour lequel est calculé ce taux de survie est noté en indice de la lettre S. [0026] Dans cette description, un rayonnement de référence Ref est identifié par l'indice ref. Un rayonnement ionisant quelconque, de particules de type T et d'énergie E, est noté par les lettres T,E. Ainsi, la notation STP,E désigne le taux de survie des cellules 6 présentant une pression P lorsqu'elles sont irradiées par le rayonnement 2. Ce taux est une fonction mathématique de la dose D. La valeur de ce taux sTP,E , pour une valeur de la dose D, est noté sTP,EiD . [0027] Par exemple, dans le modèle dit linéaire-quadratique, le taux de survie sTP,E, s'écrit, en fonction de la dose D, selon la formule : -1n(sTP,E(D))=0(Puxp+pTP,E,> A+B. L'espèce I est, par exemple, une protéine ou un antioxydant. Cette réaction conduit à la formation d'une espèce chimique radicalaire A, et d'un produit B. Cette espèce radicalaire A est potentiellement apte à réagir chimiquement à son tour 40 avec de l'oxygène 02 présent dans les cellules 6, pour former, directement ou indirectement, les éléments toxiques R. Par exemple, cette espèce chimique A réagit avec de l'oxygène 02 selon la réaction suivante : A-+02-> X02+C . Cette réaction conduit à la formation d'un produit C et, surtout, de l'espèce X02 , qui présente une toxicité importante pour les cellules 6. De telles réactions sont décrites en détail, par exemple, dans l'article de A. Colliaux et al. « 02 and glutathione effects on water rad iolysis : a simulation study », Journal of Physics: Conference Series, vol. 261, p. 012007, 2011. [0032] La production de ces espèces chimiques par radiolyse est typiquement quantifiée au moyen de rendements de production. Le rendement de production d'une espèce chimique dans un milieu donné, sous l'effet indirect de la radiolyse induite par un rayonnement, est ici défini comme étant la quantité de molécules de cette espèce produites par 100eV d'énergie transmise au milieu par ce rayonnement. Ce rendement est, par exemple, obtenu pour un milieu aqueux, tel que de l'eau, à l'aide de codes de calcul numériques décrits dans l'article, noté A2, de B. Gervais et al., « Numerical simulation of multiple ionization and high LET effects in liquid water radiolysis », Radiation Physics and Chemistry, vol. 75, n° 4, p. 493-513, 2006. Par la suite, on admet que le rendement de production obtenu pour ce milieu aqueux est représentatif du rendement de production qui serait obtenu pour le milieu intracellulaire. En effet, ce rendement qui serait obtenu pour le milieu intracellulaire est fortement corrélé au rendement obtenu pour un milieu aqueux. Ces rendements sont, par exemple, liés entre eux par un coefficient de proportionnalité. [0033] La figure 3 illustre plus en détail un exemple de courbe 140 du rendement de production, noté Gt, de l'espèce R, en fonction de la pression partielle d'oxygène P02 dans le milieu aqueux. Cette courbe 140 présente une région de croissance 141 et une région de saturation 142. Dans la région de saturation, le rendement n'augmente quasiment plus, même lorsque la valeur de la pression P02 augmente. Cette région 142 permet de définir une pression de saturation PsGat . Par exemple, cette pression PsGat est définie comme étant la valeur de la pression P02 correspondant à l'intersection entre la droite tangente à la courbe 140 à l'origine d'une part, et la droite tangente à l'infini de à cette courbe 140 d'autre part. Ici, chaque valeur C(P02) du rendement Gt de cette courbe est enregistrée au bout d'une durée t suivant l'instant d'impact d'une particule T (ou d'un paquet de particules) du rayonnement 2 sur le milieu aqueux. Ce rendement Gt varie en fonction de la valeur choisie pour la durée t. Dans cet exemple, cette durée t est comprise entre 5ps et 15ps et, de préférence, égale à 10ps. Une manière plus générale de déterminer cette durée t sera décrite plus en détail en référence à la figure 4. [0034] Il est souhaitable de connaître le taux de survie STP,E dans le cas général, c'est-à-dire pour des conditions de pression P quelconques, et/ou pour un rayonnement ionisant comportant des particules de type T et d'énergie E 40 quelconques. Un tel calcul peut, par exemple, être nécessaire pour prédire avec précision la dose D à appliquer sur les cellules 6 dans le cadre d'un traitement de radiothérapie. Cependant, ce taux de survie SE n'est pas connu dans le cas général. Cette expression pourrait être déterminée expérimentalement, mais il faudrait pour cela mesurer le taux de survie d'une population de cellules biologiques analogues aux cellules 6, et ce pour chacune des valeurs possibles de pression P, de type T et d'énergie E. De telles expériences seraient particulièrement fastidieuses et coûteuses à mettre en place. [0035] En revanche, il est possible de calculer le taux de survie des cellules 6 pour une pression P quelconque lorsqu'elles sont irradiées par un rayonnement ionisant dit de référence, noté ici Ref. On définit ainsi le taux de survie srPef qui, à toute valeur d'une dose D du rayonnement Ref, associe la valeur srPefID du taux de survie si.Pef des cellules 6, lorsqu'elles présentent la pression P et qu'elles sont irradiées par la dose D du rayonnement Ref. Ce taux SrPef peut s'exprimer en fonction de la dose D, selon le modèle linéaire-quadratique précédemment décrit. Typiquement, ce taux SrPef est connu pour des conditions normalisées et présentées, sous forme de figures de courbes de survie rassemblant des données expérimentales, disponibles dans la littérature scientifique (comme par exemple, pour la lignée cellulaire V79, dans l'article de I. J. Spiro et al., « Oxygen radiosensitisation at 10w-dose rate », British Journal of Radiology, vol. 58, n°688 p. 357-363, avril 1985, dont la figure 1.2 présente des taux de survie pour différentes pressions en oxygène), ou bien est déterminé expérimentalement. Par exemple, on peut obtenir ce taux pour la valeur P1 de pression partielle d'oxygène. On le note alors SrPeif i D . [0036] Pour quantifier la contribution de l'effet oxygène sur ce taux de survie, on définit le rapport d'amplification par oxygène, dit OER (pour « Oxygen Enhancement Ratio » en langue anglaise). Par exemple, pour les deux valeurs P1 et P2 de la pression P, et pour un taux de survie de référence (typiquement 10%), la valeur OER(P2) de ce rapport OER est définie comme le rapport : -de la valeur de dose pour laquelle le taux de survie associé à la valeur P1 de pression est égal au taux de survie de référence, -sur la valeur de la dose pour laquelle le taux de survie associé à la valeur P2 de pression P est égal au taux de survie de référence. [0037] Par exemple, dans le cas illustré sur la figure 2, le taux de survie de référence est choisi égal à 10 % et l'OER pour la valeur P2 de pression partielle est défini comme étant le quotient Dp1/Dp2. En répétant le calcul de l'OER pour différentes valeurs P2 de la pression P, il est possible de disposer d'une fonction exprimant l'OER en fonction de la pression P. [0038] La figure 4 représente un exemple d'une courbe 15 de cette fonction OER. Cette courbe 15 présente une région de croissance 17 et une région de saturation 19. Dans la région 19, l'OER n'augmente quasiment plus, même si la pression P 40 augmente. Cette région 19 permet de définir une pression de saturation 1)50aE,R . Cette pression de saturation CR est ici définie de manière analogue à la pression de saturation 1),Ga, , à partir de droites tangentes à la courbe 15 en l'origine et à l'infini. [0039] La durée t peut ainsi être déterminée à partir des courbes de rendements et d'OER en fonction de la pression P. Par exemple, plusieurs courbes de rendement Gt 5 sont enregistrées pour chaque valeur de pression P, pour différentes valeurs de la durée t. Puis, parmi toutes les courbes de rendement Gt de production d'une espèce par une dose d'un rayonnement donné, dont on dispose ainsi pour des valeurs différentes de la durée t, on choisit la courbe de rendement associée à la valeur de la durée t qui est telle que la pression de saturation PsGa, soit égale à, ou le plus proche 10 possible de la pression de saturation P500R,R . [0040] Par ailleurs, pour une valeur P2 de pression P, on peut définir un coefficient de correction CrPelfP2 qui, à toute valeur d'une dose D du rayonnement Ref, associe une valeur notée CrPelfP2CD du coefficient de correction du taux de survie srPelf , pour prendre en compte la différence d'impact de l'effet oxygène entre des valeurs de 15 pression d'oxygène P2 et P1, pour une dose D du rayonnement Ref. Ce coefficient PiP2 P P Cref permet notamment, pour chaque valeur Cref 113 associée à une dose D du rayonnement Ref, de déterminer : -une valeur SrP:f(D) du taux de survie srPe2f des cellules 6, irradiées par une dose D du rayonnement Ref et présentant une pression partielle d'oxygène P2, distincte de la 20 pression partielle P1, -à partir d'une valeur connue SrPelf(D) du taux de survie srPelf pour la même dose D. [0041] Par exemple, le coefficient CrPelfP2 est calculé pour quelques valeurs de la dose D à l'aide de la relation suivante : CrPefP2CD =SP'2f(D /SrPeif (D , où srPelf(D) et S'2f(D) sont des valeurs connues ou mesurées expérimentalement. Ensuite, pour 25 toute valeur de la dose D, il est possible de construire la valeur du coefficient CrPelfP2 par interpolation numérique, ici entre deux valeurs connues. Pour simplifier, la même convention de notation que celle définie pour le taux SE est adoptée, sauf qu'en exposant de la lettre C sont notées les pressions partielles d'oxygène entre lesquelles ce coefficient est calculé. 30 [0042] Le rayonnement de référence Ref est notamment choisi comme étant un rayonnement pour lequel le coefficient CrPelfP2 est connu, ou facilement déterminable. Par exemple, le rayonnement Ref est un rayonnement ionisant à faible TEL. Ce rayonnement Ref présente un TEL inférieur ou égal à 1,5keV/pm ou à 1keV/pm. Cette valeur de TEL est ici définie pour la majorité de la longueur de la trajectoire du 35 rayonnement Ref dans le tissu 4. Ici, ce rayonnement Ref est un faisceau de rayons X, dont l'essentiel du spectre en énergie a des composantes dont l'énergie est supérieure à 100keV. [0043] La figure 5 représente un dispositif de calcul 12, pour calculer une dose D de rayonnement 2 à fournir aux cellules 6. Ce dispositif 12 comprend : 40 -un calculateur électronique 14 programmable ; -un support d'enregistrement 16 d'informations ; -une interface 18 d'échange de données. [0044] L'interface 18 permet notamment l'acquisition et la fourniture de données numériques, tels que des représentations numériques de fonctions mathématiques. Le calculateur 14 est notamment programmé pour exécuter des instructions enregistrées dans le support d'enregistrement 16. Le support d'enregistrement 16 contient notamment des instructions pour l'exécution du procédé des figures 6 et 7. [0045] Le fonctionnement du dispositif 12 va maintenant être décrit en référence au procédé de la figure 6 et à l'aide de la figure 1. [0046] Lors d'une étape 100, des paramètres d'initialisation du procédé sont acquis par le dispositif 12. Ces paramètres comprennent notamment des caractéristiques expérimentales des cellules 6. Cette étape 100 comprend ainsi : -une opération 102 d'acquisition de la pression P des cellules 6 ; -une opération 104 d'acquisition d'informations sur le type des cellules 6, telles que 15 l'identification de la lignée cellulaire auquel appartiennent les cellules 6 ; -une opération 106 d'acquisition d'une représentation spatiale du tissu 4, et notamment de la représentation spatiale d'une zone précise du tissu 4 sur laquelle le rayonnement 2 doit être appliqué ; -une opération 108 d'acquisition du taux sTP,IE et du coefficient CrPelfP2 ; 20 -une opération 110 d'acquisition de valeurs du(P) , et Ge( P2) de rendement de production de l'espèce R dans le milieu aqueux représentant le milieu intracellulaire des cellules 6, respectivement, présentant la pression P pour le rayonnement 2 et présentant la pression P2 pour le rayonnement Ref. [0047] Le taux STP,'E, est, par exemple, obtenu expérimentalement. Le coefficient P 25 CPreit 2 est obtenu expérimentalement, ou à partir de données bibliographiques. Ce taux et ce coefficient sont ici acquis sous la forme d'une table de valeurs numériques. [0048] Dans cet exemple, l'espèce R correspond au couple formé par l'anion superoxyde 02* - et le radical hydroperoxyde H02*. Ce couple est produit par les réactions suivantes : eaq +02-> 02- et 11. + 02-> HO; , à partir des radicaux eaq et H' 30 produits par radiolyse du milieu aqueux. Les courbes de rendement GtTE et Coef correspondent donc ici à la production de ce couple superoxyde 02* - et hydroperoxyde H02*. Par exemple, ces courbes GtT,E et doef sont obtenues à partir des simulations numériques décrites dans l'article précité, noté A2, de B. Gervais et al. Ces courbes de rendement sont construites pour la durée t choisie comme 35 précédemment décrit. Les valeurs GT,t,(P) et Ge( P2) sont ensuite obtenues à partir de ces courbes. [0049] Les fonctions donnant les taux, rendements et coefficients en fonction de la dose D sont ici, lors des opérations 108 et 110, acquises sous la forme de tables de valeurs numériques, ou bien de toute autre forme qui permet d'associer, à chaque 40 valeur de la dose D, la valeur correspondante de la fonction. posol Lors d'une étape 120, une valeur prédéfinie S d'un taux de survie est acquis par le dispositif 12. Cette valeur S est par la suite simplement appelée « taux de survie S ». Ce taux de survie S correspond au taux de survie que doivent présenter les cellules 6 après irradiation par le rayonnement 2. [0051] Par exemple, ce taux de survie S est déterminé à partir d'une dose prescrite fournie par un praticien dans le cadre d'un traitement par radiothérapie, sans prendre en compte l'effet oxygène et, plus précisément, le fait que la pression d'oxygène dans les cellules 6 peut-être inférieure à la pression normale. Le taux de survie S est ensuite déterminé automatiquement à partir de cette dose prescrite, selon un modèle classique, ne prenant pas en compte l'effet oxygène et, plus précisément, la radiorésistance augmentée par le manque d'oxygène. [0052] Lors d'une étape 122, des paramètres physiques du rayonnement (tels que des propriétés énergétiques ou de focalisation) sont déterminés, pour que les cellules 6 présentent le taux de survie S après irradiation. Cette détermination est réalisée au moyen de procédés d'estimation connus. Ces paramètres comprennent notamment une valeur de la dose D du rayonnement 2 à appliquer. [0053] A ce stade du procédé, cette valeur de la dose D déterminée lors de l'étape 122 n'est cependant qu'une approximation et doit faire l'objet d'une optimisation et d'une validation avant de pouvoir être fournie à la source 8. En effet cette dose D a été déterminée lors de l'étape 122 sans prendre en compte l'effet oxygène. Aussi, la valeur de la dose D doit être optimisée, notamment, pour tenir compte de cet effet oxygène. [0054] Cette optimisation est ici réalisée en utilisant, comme contrainte d'optimisation, la valeur STP,EiD du taux de survie des cellules 6 si elles recevaient 25 cette valeur de la dose D du rayonnement 2 en tenant compte de l'effet oxygène. [0055] A cet effet, lors d'une étape 124, la valeur sTP,EiD du taux de survie sTP,E est calculée pour la dose D du rayonnement 2. Cette valeur sTP,EiD peut se décomposer sous la forme suivante : STPECD =SPT1ECD xCTP1E'PCD , où : - sTP,1E(D) est la valeur, pour la dose D du rayonnement 2, du taux sTP,IE acquis lors de 30 l'opération 108; - CTP1E'P (D est une valeur d'un coefficient de correction CTP1E'P . [0056] Ce coefficient CTP,E1'P associe, à toute valeur d'une dose D du rayonnement 2, un coefficient de correction du taux de survie sTP,1E(D) , pour prendre en compte l'impact de l'effet oxygène, dans le cas où la pression d'oxygène P des cellules 35 acquise est différente de la valeur P1 pour une dose D du rayonnement 2. [0057] Cette décomposition facilite le calcul du taux sTP,EiD . En effet, l'expression exacte des coefficients aPT,E et PPTE n'est généralement pas connue. Le calcul de la valeur sTP,EiD se ramène donc simplement au calcul de STP,I,CD et de CTP,E1'PCD . Le calcul de la valeur CTP1E'PCD est décrit en référence à la figure 7. [0058] A l'issue de l'étape 124, lors d'une étape 130, la valeur sTP,EiD calculée est comparée à la valeur du taux de survie S prédéfini, acquis lors de l'étape 120. Si la valeur sTP,EiD ne correspond pas à la valeur de S, alors la valeur de D déterminée lors de l'étape 122 est rejetée. On retourne ensuite à l'étape 122 pour déterminer une 5 nouvelle valeur de D. La valeur sTP,EiD est dite correspondre à la valeur S si elle est égale à S à moins de 10 % près, ou à moins de 5 % près ou, de préférence, à moins de 2 % près. Dans le cas contraire, si la valeur sTP,EiD correspond à la valeur de S, alors la valeur de D déterminée lors de l'étape 122 est validée. Cette valeur de D validée est ici fournie sur l'interface 18 lors d'une étape 132 et est utilisée pour 10 commander la source 8, de manière à appliquer la dose D du rayonnement 2 sur les cellules 6. Par exemple, l'interface 18 transmet automatiquement la dose D à la source 8, qui applique ensuite cette dose D du rayonnement 2 sur les cellules 6. [0059] La figure 7 décrit plus en détail l'étape 124. [0060] Lors d'une opération 150, la quantité Qi représentative de la quantité de 15 radicaux R produits dans les cellules 6 est calculée, lorsque ces cellules 6 sont irradiées par la dose D du rayonnement 2, ces cellules 6 présentant la pression P d'oxygène. Cette quantité Qi est ici calculée, à partir de la valeur du(P) acquise lors de l'opération 110, par la formule suivante : =G;-,E (P )x D - 20 [0061] En radiothérapie, et plus spécifiquement en hadronthérapie, l'ordre de grandeur des valeurs des doses de rayonnement est typiquement inférieur ou égal à 200Gy et, de préférence, inférieur ou égal à 50Gy. Dans cette gamme de valeurs de doses, il est ainsi permis d'approximer, avec une précision suffisante, la quantité Qi comme étant égale au produit de la valeur du(P) par la dose D. 25 [0062] Lors d'une opération 152, une dose équivalente D est déterminée. Cette dose D' est définie comme étant la valeur de la dose du rayonnement Ref qu'il faudrait appliquer sur ces cellules 6 pour produire la même quantité Qi représentative de la quantité de radicaux R dans les cellules 6 présentant une valeur P2 de pression P. 30 [0063] Cette dose D' est ici déterminée, à partir de la valeur Ge( P2) acquise lors de l'opération 110, par la formule suivante : D'- Gt'f ,P2) 35 [0064] Cette dose D' permet, en passant par les données radiochimiques liées à la radiolyse du milieu aqueux représentant le milieu intracellulaire des cellules 6, de ramener le calcul de la valeur CTPEPP (D) du coefficient CTP.EPP au calcul de la valeur CrijelfP2CD' du coefficient CrPeifP2 Selon la formule suivante : CTP'EPI(D)=CP'liP2(D1 . [0065] En effet, on considère ici que l'effet oxygène dépend uniquement de la 40 quantité de radicaux produits dans les cellules 6. Par contre, que cette quantité de radicaux ait été produite par la dose D de rayonnement 2 ou par la dose D du rayonnement Ref est considéré comme sans conséquence sur l'effet oxygène. Dès lors, la valeur Cul' (D du coefficient CT1' peut être prise égale à la valeur P P( P P - - - - P ,P P ,P Crelf ID' du coefficient Crplf. . Ainsi, ici, la valeur Cul ( du coefficient Cul est égale à la valeur CrijelfP2CD' . Le passage des valeurs P1 et P2 à une valeur quelconque de la pression P est rendu possible du fait que les valeurs G6(P) et C'f(1)2) prennent en compte la pression P. [0066] Le procédé décrit ci-dessus permet de calculer la valeur CTP,ErPCD du coefficient CTPI'P pour un rayonnement T, E quelconque, sans qu'il ne soit nécessaire de réaliser de nombreuses expériences pour construire des valeurs du le taux de survie STP,E, pour ce rayonnement T,E quelconque. En effet, ce coefficient CrPeifP2 présente l'avantage de pouvoir être calculé à partir de données de référence, disponibles dans la littérature et qui ont déjà été déterminées, au prix d'expériences lourdes et longues à mettre en oeuvre. Calculer littéralement directement ce coefficient CTP1E'P impliquerait notamment de devoir refaire toutes ces expériences. [0067] Ainsi, lors d'une opération 154, une valeur CrPelfP2CD' du coefficient C1/j' est calculée pour la dose D. [0068] Enfin, lors de l'opération 156, la valeur STP,EiD du taux de survie STP,E est calculée pour la valeur de la dose D, suivant la formule suivante : 20 T SPECD =C.1 1j2E r 'P(D xS.PiE (D =CrP ef T 1P2CD' xSP1ECD . [0069] De nombreux autres modes de réalisation sont possibles. [0070] En variante, le rayonnement 2 comporte des particules de type et/ou d'énergie différents de ceux décrits. Par exemple, les particules de type T sont des ions formés à partir d'un élément atomique quelconque. Par exemple, les particules 25 de type T sont des ions hydrogène, hélium ou lithium. [0071] Les cellules 6 peuvent être d'un type différent de celui décrit. Par exemple, les cellules 6 sont des cellules tumorales du sein, des poumons, du système nerveux central, du cerveau, ou encore de la région tête et cou du corps humain. [0072] Il existe d'autres procédés de calcul dans lesquels l'étape 124 peut être 30 utilisée. Par exemple, l'étape 124 est utilisée dans un procédé ne comprenant pas d'étapes d'optimisation. Le taux de survie est alors estimé directement. [0073] L'espèce chimique R peut être différente. En fait, n'importe quelle espèce chimique R dont la quantité dans la cellule est corrélée avec l'effet oxygène peut être utilisée. De préférence, l'espèce chimique R est choisie comme étant un radical dont 35 la formation implique, d'une part, l'électron aqueux eaq ou le radical hydrogène H' et, d'autre part, l'oxygène 02. Par exemple, l'espèce chimique R est un radical G02* ou un radical GS02' issu d'une réaction chimique de l'oxygène avec, respectivement, le radical G' ou le radical GS*. Ces radicaux G' et GS* sont eux-mêmes issus d'une réaction chimique du glutathion GSH avec l'électron aqueux ou avec le radical 40 hydrogène. [0074] Les taux de survie peuvent être exprimés littéralement suivant un modèle différent du modèle linéaire-quadratique. [0075] Typiquement, le rayonnement 2 peut également comprendre au niveau des cellules 6, en plus des particules de type T et d'énergie E, une pluralité de particules secondaires. Ces particules secondaires sont, par exemple, issues de réactions d'interaction non souhaitées du rayonnement 2 avec son environnement, lors de la propagation de ce rayonnement 2 vers les cellules 6. Un tel phénomène est connu sous le terme de « fragmentation ». Par exemple, dans le cas où les particules de type T sont des ions carbone 12L,,-,6+, ces particules secondaires sont des noyaux atomiques dont le numéro atomique est inférieur ou égal à celui du carbone, tels que le lithium, le béryllium ou le bore. En toute rigueur, la contribution de ces particules secondaires devrait être prise en compte dans le calcul du taux de survie. Toutefois, cette contribution présente généralement une contribution en dose inférieure ou égale à 10 % ou à 5 % de la dose totale déposée par le rayonnement 2. Aussi, pour simplifier, seul le calcul du taux de survie STP,E, pour les particules de type T a été décrit en détail ici. Cependant, en variante, la contribution des particules secondaires peut être prise en compte pour calculer le taux de survie des particules 6. Pour ce faire, l'étape 124 est appliquée successivement pour chacun des types de particules secondaires, de la même manière que pour les particules de type T. Les taux de survie ainsi obtenus peuvent ensuite être combinés entre eux et avec le taux SE selon une méthode dite des champs mixtes (« mixed radiation fields » en langue anglaise). Un exemple de méthode des champs mixtes est décrit à la section 2.3 de l'article précité, noté Al, de M. Krâmer et M. Scholz. [0076] En variante, le rayonnement 2 peut comprendre, au niveau des cellules 6, des particules présentant une pluralité d'énergies E distinctes de l'énergie Eo. Typiquement, des particules d'énergie Eo du rayonnement 2 émises par la source 8 subissent un certain nombre d'interactions aléatoires au cours de leur déplacement depuis la source 8 vers les cellules 6, ce qui modifie leur énergie E. Ainsi, au niveau des cellules 6, le rayonnement 2 comporte des particules de type T mais dont les énergies E se répartissent suivant une distribution statistique, typiquement centrée sur Eo. En toute rigueur, la contribution de ces particules dont l'énergie est différente de Eo doit être prise en compte lors de l'étape d'optimisation du calcul de la dose, et donc dans le calcul du taux de survie. Ainsi, en variante, une méthode de champs mixtes est mise en place pour calculer un taux de survie global, en combinant les valeurs STP,EiD du taux de survie sTP,E , calculées pour chacune des énergies E. Cette méthode des champs mixtes est par exemple celle décrite dans l'article noté Al. Cette variante s'applique également au cas où la source 8 est configurée pour émettre un rayonnement comportant des particules avec plusieurs énergies Eo distinctes. [0077] Le rayonnement 2 peut présenter un TEL quelconque. [0078] Les pressions P1 et/ou P2 peuvent prendre des valeurs différentes, selon le type des cellules 6. [0079] Par exemple, les fonctions donnant les taux, rendements et coefficients en fonction de la dose D sont acquises sous la forme de valeurs d'interpolations 5 numériques de ces fonctions. posol En variante, la durée t présente une valeur différente. Cette durée t peut aussi être déterminée différemment. Par exemple une valeur de la durée t est systématiquement choisie égale à une valeur fixe, par exemple égale à 10ùs. [0081] Les courbes de rendement GtTE et dref peuvent être déterminées de façon 10 différente. Par exemple, ces rendements sont mesurés expérimentalement. [0082] La quantité Qi peut être déterminée différemment. Par exemple, cette quantité est obtenue à partir d'abaques définis expérimentalement. [0083] D'autres formules pour calculer la valeur STP,ECD sont possibles. Par exemple, le taux sTP,1E peut ne pas être connu. Dans ce cas, lors de l'opération 108, 15 un taux connu sTP,°E est acquis à la place pour une valeur Po de pression, distincte de P1. De plus, lors de l'étape 124, une valeur CTP:EP°CD d'un coefficient de correction additionnel est introduite, de la même manière que pour la valeur CTP,El'PCD , sauf que la valeur de la pression P est remplacée par la valeur Po. Dans ce cas, la valeur sTP,EiD est calculée au moyen de la formule suivante : 20 ( SP (D =SP° (D X DT,E = S T1JD X CP P2CD T,E T,E ref CTE0(D où CrPe'f.P2CD ' est une valeur d'un coefficient de référence C rPe'fij , qui peut être calculé à partir du coefficient CrPplfP2 . 25 [0084] La valeur STP,EiD du taux de survie peut ainsi être calculée de manière analogue aux opérations de l'étape 124, mais en substituant la valeur Po à la valeur P1 de pression et en définissant une dose équivalente D" de façon similaire à la dose équivalente D mais par rapport à la valeur Po de pression. [0085] En effet, la valeur STP,°E(D) peut s'écrire de la façon suivante : 30 STP1CD =SPT1ECD xCT1E' 113)=STE1 (D xCrPeliP2(13 . [0086] En outre, la valeur STPEiD peut s'écrire de la façon suivante : STP.ECD =SPTIECD xCPT.IE'PCD =SPTIECD xCPPP2(D'). CP1;P2(13 [0087] On en déduit donc l'expression STP,E,CD =SPT°E(D x ref p 2(D et, de là, Crei ") l'expression de la valeur CrP:f.P2CD ' comme suit : 35 P P ( CP P2 (D') = ref C )3 Cref DI I) ref Pl)P2( REVENDICATIONS1. Procédé d'estimation automatique d'une valeur STP,ECD d'un taux de survie de cellules biologiques irradiées par une dose D d'un premier rayonnement ionisant 5 comportant des particules de type T et d'énergie E, ces cellules présentant une pression partielle P d'oxygène, caractérisé en ce que ce procédé comporte : -l'acquisition (108) d'un taux de survie STP,1E, des cellules biologiques lorsqu'elles sont irradiées par une dose D du premier rayonnement ionisant et qu'elles présentent une pression partielle P1 d'oxygène, chaque valeur STP,1E,CD de ce taux de survie STP,iE 10 étant fonction de la valeur de la dose D ; -l'acquisition (108) d'un coefficient de correction CrPelfP2 , ce coefficient étant une fonction de la valeur de la dose D, chaque valeur CrPelf.P2CD de ce coefficient CrPelfP2 pour une valeur particulière de la dose D exprimant la dépendance du taux de survie des cellules à la pression partielle d'oxygène dans les cellules, lorsque ces cellules 15 sont irradiées par cette dose D d'un second rayonnement ionisant, ce second rayonnement ionisant comportant des particules d'un type différent du premier rayonnement ionisant, chaque valeur CrPelf.P2CD permettant de convertir : - une valeur SP1 (D d'un taux de survie SP1 des cellules biologiques, irradiées ref ref par la dose D du second rayonnement ionisant, lorsqu'elles présentent une 20 pression partielle P1; - en une valeur SrPe2f(D d'un taux de survie Sre2f de ces mêmes cellules biologiques, irradiées par la même dose D du second rayonnement ionisant lorsqu'elles présentent une pression partielle d'oxygène P2, distincte de la pression partielle P1; 25 -la détermination (150) d'une quantité Qi représentative d'une quantité d'une espèce chimique donnée R, produite par radiolyse du milieu intracellulaire des cellules biologiques sous l'effet de la dose D du premier rayonnement ionisant, lorsque ces cellules présentent la pression partielle P d'oxygène ; -la détermination (152) d'une dose D du second rayonnement ionisant qu'il faudrait 30 appliquer sur les mêmes cellules biologiques pour obtenir, par radiolyse dudit milieu cellulaire, la même quantité Qi représentative de la quantité de l'espèce chimique R, lorsque ces cellules présentent la pression partielle P2d'oxygène ; -le calcul (154) d'une valeur CPf PP2(Dr) du coefficient de correction CP1P pour la re re dose D' ; 35 -le calcul (156) de la valeur STP,ECD pour la dose D, selon la relation : STPECD =CrPeif2CD' xSTPiECD . 2. Procédé selon la revendication 1, dans lequel l'espèce chimique R est un radical libre produit par une ou une suite de réactions chimiques dont un réactif est une molécule d'oxygène et un autre réactif est un électron aqueux eaq ou un radical hydrogène H', cet électron aqueux eaq ou ce radical hydrogène N'étant un produit de la radiolyse du milieu intracellulaire des cellules par la dose D du premier rayonnement. 3. Procédé selon la revendication 2, dans lequel l'élément R est compris dans le groupe composé de l'anion superoxyde 02-, du radical hydroperoxyde H02 et des radicaux G02* et GS02' issus,respectivement, de réactions de l'oxygène avec des radicaux G' et GS, ces radicaux G' et GS* étant issus eux-mêmes des réactions du glutathion GSH avec l'électron aqueux ou avec le radical hydrogène. 4. Procédé selon l'une quelconque des revendications précédentes, dans lequel le 15 coefficient de correction CrPelf.P2 est défini, pour toute dose D, par la formule : CrPeif:P2( D =SP'2f(D ISrP elf(D 5. Procédé selon l'une quelconque des revendications précédentes, dans lequel le second rayonnement présente un transfert d'énergie linéique inférieur ou égal à 1 20 keV/pm sur la majorité de la longueur de sa trajectoire dans le milieu contenant les cellules biologiques. 6. Procédé selon l'une quelconque des revendications précédentes, dans lequel le premier rayonnement présente un transfert d'énergie linéique supérieur ou égal à 25 10keV/pm au niveau des cellules. 7. Procédé selon l'une quelconque des revendications précédentes, dans lequel les particules de type T sont des hadrons. 30 8. Procédé selon la revendication 7, dans lequel les particules de type T sont des ions carbone 1206+ . 9. Procédé de calcul automatique d'une valeur d'une dose D d'irradiation de cellules biologiques avec un rayonnement ionisant, pour que ces cellules présentent, après 35 irradiation, une valeur S prédéfinie d'un taux de survie, ce procédé comportant : a] l'acquisition (120) de la valeur S prédéfinie et d'une cartographie des valeurs de la pression partielle P d'oxygène des cellules biologiques à irradier ;b] le calcul (122) de paramètres du rayonnement ionisant nécessaires pour que les cellules biologiques, une fois irradiées avec ce rayonnement ionisant, présentent la valeur S prédéfinie de taux de survie, ce rayonnement ionisant comprenant, au niveau des cellules biologiques , des particules de type T et d'énergie E, ce calcul comprenant notamment la détermination de la valeur de la dose D de ce rayonnement à appliquer sur les cellules biologiques ; c] l'estimation (124) d'une valeur STP,ECD d'un taux de survie SPT,E des cellules biologiques irradiées par la dose D du rayonnement ionisant déterminée lors de l'étape b] avant application de cette dose D du rayonnement sur les cellules 10 biologiques : d] la comparaison (130) de la valeur STP,ECD estimée lors de l'étape c] avec la valeur S prédéfinie, puis, - la fourniture (132) de la valeur de la dose D calculée lors de l'étape b] si la valeur estimée STP.ECD correspond à la valeur S prédéfinie, et, 15 - dans le cas contraire, la réitération de l'étape b] pour déterminer une nouvelle valeur de la dose D; caractérisé en ce que la valeur SPT,ECD du taux de survie STP,E est estimée selon un procédé conforme à l'une quelconque des revendications précédentes. 20 10. Support d'enregistrement d'informations (16), caractérisé en ce qu'il comprend des instructions pour l'exécution d'un procédé conforme à l'une quelconque des revendications précédentes lorsque ces instructions sont exécutées par un calculateur électronique. 25 11. Calculateur électronique (12) pour la mise en oeuvre d'un procédé d'estimation conforme à l'une quelconque des revendications 1 à 8, ce calculateur étant programmé pour : -acquérir d'un taux de survie STP,1E, des cellules biologiques lorsqu'elles sont irradiées par une dose D du premier rayonnement ionisant et qu'elles présentent une pression 30 partielle P1 d'oxygène, chaque valeur STP,1E,CD de ce taux de survie STP,1E, étant fonction de la valeur de la dose D ; P1 -acquérir un coefficient de correction CP2 ce coefficient étant une fonction de la ref valeur de la dose D, chaque valeur CrPeifP2CD de ce coefficient CrPeifP2 pour une valeur particulière de la dose D exprimant la dépendance du taux de survie des cellules à la 35 pression partielle d'oxygène dans les cellules, lorsque ces cellules sont irradiées par cette dose D d'un second rayonnement ionisant, ce second rayonnement ionisantcomportant des particules d'un type différent du premier rayonnement ionisant, chaque valeur CrPelf.PID permettant de convertir : - une valeur SrefCD d'un taux de survie Sreit des cellules biologiques, irradiées par la dose D du second rayonnement ionisant, lorsqu'elles présentent une pression partielle P1; - en une valeur SrPe2f(D d'un taux de survie Sre2f de ces mêmes cellules biologiques, irradiées par la même dose D du second rayonnement ionisant lorsqu'elles présentent une pression partielle d'oxygène P2, distincte de la pression partielle P1; -déterminer une quantité Qi représentative d'une quantité d'une espèce chimique donnée R, produite par radiolyse du milieu intracellulaire des cellules biologiques sous l'effet de la dose D du premier rayonnement ionisant, lorsque ces cellules présentent la pression partielle P d'oxygène ; -déterminer une dose D du second rayonnement ionisant qu'il faudrait appliquer sur 15 les mêmes cellules biologiques pour obtenir, par radiolyse dudit milieu intracellulaire, la même quantité Qi représentative de la quantité de l'espèce chimique R, lorsque ces cellules présentent la pression partielle P2 d'oxygène ; -calculer une valeur CrPeliP2CD' du coefficient de correction C rPelfP pour la dose D' ; -calculer la valeur STP,ECD pour la dose D, selon la relation : ,P 20 ST-E,,13 =C'f xSTP1ECD .