Patent Application
20230172501
The invention relates to methods and devices for monitoring the concentration of a substance, preferably, oxygen, in a cell or tissue, e.g., in cells of the human skin. In particular, it provides a method for determining the concentration of a quencher, such as oxygen and/or the concentration of probe, e.g., a heme precursor such as protoporphyrin IX (PpIX), wherein the probe is capable of exhibiting luminescence (delayed fluorescence or phosphorescence) and or transient triplet absorption, preferably, delayed fluorescence in a living cell. The method comprises steps of exciting the probe, measuring the lifetime of the luminescence, e.g., delayed fluorescence, exhibited by said probe, wherein, in the presence of the quencher, the lifetime of an excited triplet state is shortened as compared to the lifetime of an excited triplet state in the absence of the quencher, and correlating said lifetime with said concentration. The invention optimizes the method and leads to more precise results than conventional methods, because of adaptations based on the understanding of the influence of the concentration of the probe, of its triplet state, and its excitation fluence rate (intensity) and fluence on the analysis. Taking into account second order triplet interactions also permits the interpretation of non-exponential decays and further improvement of the quencher and probe concentration estimation. Methods of the invention rely, e.g., on measurement at different emission wavelengths and application of an Adaptive Stern-Volmer relationship, the decay central fitting method and/or a mixed orders approach. Said method can be applied, e.g., for assessment or imaging of tissue oxygenation, viability and/or oxygen utilization, diagnosis of conditions, assessment of treatment effects or titration of treatments, bedside monitoring of patients, for example, for intensive care patients, patients having sepsis or patients undergoing surgery, organ transplants or a tumor therapy such as radiotherapy or photodynamic therapy. The invention also relates to the use of the PpIX precursor 5-aminolevulinic acid (5-ALA), or derivatives thereof, in this method. A device suitable for the method of the invention is also provided.
Health control, diagnosis of disease and/or monitoring of treatment of disease often involves measurement of various parameters. One parameter is the concentration of a certain substance, such as oxygen, within at least part of an organism. Local tissue oxygenation is an important parameter in the diagnosis and treatment of a wide range of diseases. Measurements of the amount of oxygen present in a specific part of a subject are for instance carried out during peri-operative monitoring in the operating room and intensive care and for diagnosis of a wide range of clinical disorders in which tissue oxygenation is central to the development and cure of disease. Examples include diagnosis of cardiovascular disease, monitoring healing of decubitus and diabetic wounds, monitoring hyperbaric correction of radiation wounds and assessment of success of bypass surgery. Monitoring of tissue oxygen pressure (pO2) during critical illness is considered a major need in the adequate treatment of intensive care patients (Siegemund et al., 1999). Assessment of tumor oxygenation is also of interest, as oxygen is an important determinant for success of radiotherapy and photodynamic therapy. Hence, the concentration of oxygen in a tumor is preferably assessed to determine whether radiotherapy, or photodynamic therapy, is recommended and/or to adapt the radiation or light dose. For example, if the concentration of oxygen in parts of the tumor is low (e.g., the partial oxygen pressure being below 10 mm Hg as demonstrated by Epel el al., 2019, Oxygen-Guided Radiation Therapy. Int J Radiat Oncol Biol Phys. 103(4):977-984), an enhanced local radiation can be used to increase treatment success.
Dioxygen is a molecule of utmost biological importance because of its role as the primary biological oxidant. Therefore, oxygen plays a key role in the oxidation/reduction reactions that are coupled to cellular respiration and energy supply. For instance approximatively 22 mol of dioxygen per day (0,254 μmol.s−1) is consumed for a 70 kg 20-year-old male, where, at least 90% is reduced to water into the mitochondria for ATP production (Wag-ner et al., 2011, Free Radical Biology and Medicine 51(3):700-712). Adequate measurement of oxygen concentrations in biological samples like cells, tissues and whole organs is important to gain insight in the determinants of oxygen supply and utilisation under normal and pathological conditions. The clinical interest in methods providing information about blood-flow and oxygen delivery at the (sub-)organ level (e.g. microcirculatory) is growing. This is amongst other things because of increasing insight into the role of the microcirculation in pathogenesis, and the importance of adequate tissue perfusion as end-point of treatment (Siegemund et al., 1999. Intens Care Med 25:1044-1060).
Various techniques have been developed for direct or indirect oxygen measurements in tissue, each having its specific advantages and disadvantages (J.M. Vanderkooi et al., 1991.Am J Physiol 260:C1131-50). Conventionally, measurements of tissue oxygenation have been made by use of oxygen electrodes and spectrophotometry of the hemoglobin or myoglobin molecule. Reflection spectrophotometry records the difference in absorption and scattering between a standard reference sample and a tissue sample. The method is based on the illumination of a tissue sample by light with a known spectral content and detection of the diffusely reflected light from the tissue at several different wavelengths. Drawbacks of these conventional techniques are that they are either me-chanically disruptive (insertion of oxygen electrodes) or qualitative (spectrophotometry). These constrains have led to the development of alternative methods.
One of the most promising techniques in this respect has been the use of oxygen dependent quenching of phosphorescent dyes for measurements in the microcirculation (Vanderkooi et al., 1987. J Biol Chem 262:5476-5482; Vanderkooi et al., 1987. Science 236:568-569; Sinaasappel & Ince, 1996 J Appl Physiol:2297-2303; Sinaasappel et. al., 1999. J Physiol 514: 245-253). The technique is based upon the principle that a molecule, most frequently metallo-porphyrins, that has been excited by light can either re-lease this absorbed energy as light (phosphorescence) or transfer the absorbed energy to oxygen (without light emission). This results in an oxygen-dependent phosphorescence intensity and lifetime. The relationship between the lifetime and the oxygen concentration is usually given by the Stern-Volmer relationship (Vanderkooi, 1989. Biochim Biophys Acta 976:1-27). Calibration constants associated with the Stern-Volmer relationship allow oxygen concentrations to be calculated from the measured lifetimes. The measurement of lifetimes instead of intensities presents many advantages, in particular because it allows quantitative measurements without the influence of tissue optical properties.
Improved methods for determining the concentration of oxygen in living cells, tissue or organs, which, e.g., do not necessitate the injection of metallo-porphyrin, were developed by Mik et al. These methods are based on excitation of heme precursors such as Protoporphyrin IX (PpIX). PpIX is a natural aromatic molecule comprising four pyrrole rings produced inside the mitochondria as an intermediate in the biosynthesis of heme. PpIX is naturally present in mitochondria inside cells at relatively low concentrations. However, for measuring pO2 in cells or tissues, it is advantageous to increase the endogenous syn-thesis of PpIX to raise its concentration to levels facilitating the analysis by measurement of delayed fluorescence lifetime of PpIX. PpIX is generated in tissue in increased amounts after administration of a precursor, e.g., 5-aminolevulinic acid or derivatives thereof, e.g., with a transdermal patch for skin applications, or systemically (e.g. with a solution administered orally, or intravenously) for various solid tumours. Formulations or derivatives can include pharmaceutically acceptable salts (e.g., hydrochloride or phos-phate) or esters (e.g., hexaminolevulinate or methyl aminolevulinate) of aminolevulinic acid.
PpIX presents a strong spin-orbit coupling between the first excited singlet S1 and triplet T1 states (Balzani et al., 2014, Photochemistry and Photophysics. Concepts, Research, Applications. Angew. Chemie Int. Ed. 53:8817-8817, Section 3.6), and, thus a strong delayed fluorescence. Delayed fluorescence is a radiative molecular desexcitation which has the same spectrum as prompt fluorescence, but radiates at much longer time scales. It is the result of four consecutive steps, illustrated by a simplified Jablonski diagram as shown e.g., in
According to the recited state of the art, the DF intensity (IDF) following excitation of PpIX molecules over time (t) follows:
I
DF(t)=I0e−t/τ (1)
With I0 the initial DF intensity, T the first excited triplet state (T1) lifetime, or, equivalently, the DF lifetime (since the fluorescence lifetime is much shorter than the T1 lifetime). While in the T1 state, PpIX interacts with surrounding oxygen molecules. Ground-state triplet oxygen 3O2 efficiently captures the energy of T1, reaching the excited and reactive singlet oxygen state 1O2. This reaction, called quenching by oxygen, can be written as
and reduces the triplet lifetime. kq0
The calculation of the partial pressure of oxygen, pO2 also requires the knowledge of T0, the T1 lifetime in absence of oxygen.
As explained above, phosphorescence also is a triplet-state-based luminescence. Accordingly, instead of DF lifetime, phosphorescence lifetime can be measured, and pO2 calculated on this basis. However, PpIX phosphorescence is more difficult to detect in cells in tissue. Formula (1) can be used in an analogous way.
A further option is the detection of transient triplet-state absorption (also called triplet-triplet absorption). In this case, it is not emission of light that is monitored, rather the time-resolved light absorption of excited triplet states, measured at the same time scale as DF or phosphorescence. Similarly, equation 2 can be used to calculate the pO2, τ being the first excited triplet state T1 lifetime in that case.
Thus, the lifetime of the luminescence and/or transient absorption exhibited by said PpIX, in particular, delayed fluorescence and/or triplet-triplet absorption, can be measured. In the presence of oxygen, the lifetime of an excited triplet state is shortened as compared to the lifetime of an excited triplet state in the absence of oxygen. Accordingly, said lifetime is correlated with said concentration, and consequently, the concentration of oxygen can be determined (EP1742038 and WO 2007/00487, Mik et al., Nature Methods 3: 939-945, 2006; Mik et al. Biophys. J. 95(8):3977-90, 2008 (incorporated herein by reference).
While, in organic solutions, second order mechanisms have been found to play a role in delayed fluorescence lifetime of PpIX (Vinklarek et al., 2017. Photochem. Photobiol. 16:507-518), it is believed in the art that the advantage of lifetime measurements in cells is the independence from the concentration of the chromophore, making quantitative measurements possible in vivo, where the precise concentration of said chromophore cannot be predicted (WO 2007/004873, Vinklarek et al., 2017. Photochem. Photobiol. 16:507-518).
Based on this, methods for assessment of mitochondrial function in a sample including a tissue or organ were developed, wherein the measurement of the mitochondrial oxygen consumption rate (OCR or VO2), alternatively expressed as oxygen disappearance rate (ODR), in particular upon restricting or ceasing the supply of oxygen to the tissue by tem-porary occlusion of microvessels therein, can be used as an indicator of mitochondrial dysfunction, i.e., for assessment of health or disease (WO 2010/019041; Harms et al., 2013. Mitochondrion 13:507-514).
A device designated COMET (Cellular Oxygen METabolism) suitable for assessment of the partial pressure of oxygen in skin, e.g., for clinical measurements in patients, is commercially available from Photonics Healthcare (Utrecht, NL), and employs methods developed by Mik et al.
Thus, the correlation between the delayed fluorescence, phosphorescence or transient triplet-state absorption of probes and suitable quenchers, e.g., between the optical properties of heme precursors such as PpIX and oxygen acting as a quencher thereof, is well known.
However, the inventors noted there is a relative imprecision of measurements using the methods developed by Mik et al. under certain conditions, in particular following long ALA application, high PpIX excitation intensities or a slight change in the measurement spot. For example, the inventors observed that there were high variations with regard to the value of the DF lifetime of intracellular PpIX in absence of oxygen T0 observed in different studies (see eq. (2)). From study to study, values from 200 to 800 μs have been reported. When measuring the DF lifetime on a volunteer at different times during the day, the DF lifetime obtained was not constant, and showed a systematic, relatively re-producible behaviour: starting at a high value after 3 h of ALA application, it decreased down to a minimum around 4 to 5 h, before increasing again. A typical curve is presented on
In light of these imprecisions, the inventors addressed the problem of providing still more precise methods for determining the concentration of oxygen in living cells, tissue or organs, e.g., of human patients.
This problem is solved by the present invention, e.g., by the claimed subject matter.
The invention provides a method for determining the concentration of a quencher (preferably oxygen) and/or the concentration of a probe capable of exhibiting a triplet-state based luminescence or transient triplet absorption (preferably a heme precursor such as PpIX) in a cell, wherein the cell optionally is in a tissue or organ, comprising steps of
In the context of the present invention, different pairs of probes and quenchers may be used. After excitation, some probes (probe substances) are capable of exhibiting a tri-plet-state (T1) based luminescence (i.e., luminescence generated when a molecule in a first excited triplet state T1 is returned to the ground state S1, in a way that involves the emission of a photon, also called delayed fluorescence DF if via intermediate states or triplet phosphorescence or simply phosphorescence (PH) if direct) or transient triplet absorption (TTA, in which a photon of suitable energy is absorbed by a molecule of a probe substance in an excited triplet state (e.g. T1) so that the molecule transitions to a higher excitation triplet state (e.g., T2 or T0). A triplet state is a quantum state of a system such as a molecule of a probe substance with a multiplicity (spin multiplicity) of 3 i.e., with two unpaired electrons. An excited triplet state is triplet state in a system where at least one electron is not in the lowest possible state (i.e., not the ground state). For most molecules all triplet states that exist are excited triplet states, and for the probe we assume a substance where the ground state is not a triplet state.
A quencher is a substance whose molecules are capable of dynamic quenching. In this case that means capturing the energy of molecules of a probe substance that are in an excited triplet state T1. The presence of a quencher reduces the mean time molecules of the probe remain in triplet state (called triplet state lifetime or lifetime of the triplet state of the probe). Thus, in the presence of a suitable quencher for the probe, the triplet-state based luminescence (e.g., delayed fluorescence or phosphorescence) or transient triplet absorption decays more quickly compared to the decay in the absence of the quencher.
For example, when the quencher is oxygen (the ground state of oxygen, O2, is a triplet state at commonly encountered temperatures), the probe preferably is a heme precursor such as protoporphyrin IX (PpIX). Other heme precursors that are alternative to derivatives of 5-aminolevuinic acid that are precursors of PpIX, are molecules such as Proto-porphyrinogen IX. The probe can also be a mixture of heme precursors such as Proto-porphyrinogen IX and PpIX. As known in the art and described herein, oxygen also influences the luminescence emitted by the tetra-pyrols heme precursor as well as the triplet-triplet-state absorption thereof. Preferably, the temporal evolution of the delayed fluorescence is determined.
When the quencher is oxygen, alternative probes may be phosphorescent Pt(II)- and Pd(II)-porphyrins (Papkovsky DB et al., J. of Fluo, 15(4):569-584, 2005), phosphorescent Ru (II) complexes (Geddes CD, et al., 2005. Reviews in fluorescence, vol 2. Ed. Geddes C and Lakowicz JR (New York: Springer):125-151) or cyclometallated complexes of Ir(III) aI (Zitova A, et al., 2010. Analytical Biochemistry, 397(2):144-151). The characteristics of exemplary probes which may be used in the present invention are given in Table 1 of Dmitriev R, Papkovsky D 2012. Cell. Mol. Life Sci., 69:2025-2039.
Other probes for measuring parameters like temperature, calcium or pH, based on spectral variations such as absorption or emission peak shifts and/or fluorescence/phosphorescence lifetime variations are known in the art (and commercially available, e.g., https://www.thermofisher.com/de/de/home/brands/molecular-probes.html). Probes and quenchers are also disclosed, e.g., by Sarder et al. (2015. Bioconjug Chem. 26(6):963-974).
The method of the invention comprises a) exciting the probe by irradiation with light having an excitation wavelength and a fluence rate. These are chosen to correspond to the spectral properties of the probe. Excitation of the probe leads to formation of the excited triplet state of at least a fraction of the probe molecules in the cell that allows for detection of the triplet-state luminescence or transient triplet absorption. Suitable excitation wavelengths and fluence rates for the probes are known in the art. For example, the absorption spectrum of PpIX is provided in
Typically, heme precursors such as Protoporphyrin IX are excited at a wavelength rang-ing between 390 and 430 nm, and/or in the range of 490-550 nm, e.g., 510-530 nm or about 515 nm. In the context of wavelengths, “about” means+/−10 nm. To assess the temporal evolution of the detected signals, the excitation light can be pulsed. The pulses are typically spaced with a timing that allows for full desexcitation of the system, e.g., by about 1 ms or more. However, other modes to assess the temporal evolution of the signals are possible, including high-repetition pulsing and frequency domain approaches.
The excitation fluence rate, i.e. the power per surface unit for excitation is also suitable for formation of the excited triplet state of at least a fraction of the probe molecules in the cell that allows for detection of the triplet-state based luminescence or transient triplet absorption. Typical excitation fluence, which is the time integral of the fluence rate, are well below 1 J/cm2 at 405 nm, which is equivalent to approximately 10 J/cm2 at 515 nm, to avoid forming PpIX photoproducts (V. Huntosova et al., 2016. Journal of Photochemistry and Photobiology B: Biology 164:49-56.).
For example, a probe in skin cells may be excited by a series of 3 to 10 27 ns pulses at 515 nm. The fluence is about 300 μJ/cm2, and the energy per pulse about 15 μJ. In that case, the fluence rate is typically in the order of 10 KW/cm2. For example, a Laser Innolas Mod3 Yb-Yag (passive cooling, frequency doubling, passive Q-switch, fiber coupling (SMA)) may be used. The excitation light may be transmitted via a light diffuser to illuminate the skin, such as a frontal diffuser producing a circular spot (e.g., with an illuminated surface of about 5 mm2), with 1 ms separating each pulse. About generally means+/−one order of magnitude.
The method further comprises b) measuring the temporal evolution of the first excited triplet state of the probe, wherein, in the presence of the quencher, the triplet-state decays more quickly compared to the decay in the absence of the quencher. Thus, the triplet state lifetime is measured. The triplet state lifetime is typically assessed by measuring the triplet-state based luminescence or transient triplet absorption of the probe at at least one emission or absorption wavelength. In the presence of the quencher, the triplet-state based luminescence or transient triplet absorption decays more quickly compared to the decay in the absence of the quencher.
To assess the temporal evolution, the luminescence intensity is measured over time, i.e., at a plurality of timepoints, typically every 20-50 ns, in the relevant time range. For example, for a full resolution of faster decays, notably, during the first part of the decay at a high pO2 level, every 1-10 ns. Repetitive measurements of such pluralities of timepoints, e.g., every ms for DF measurements, may be performed. Alternatively, a “continuous” measurement can be carried out with a detector presenting a response that is modulated at the same frequency (homodyne detection) as the excitation intensity modulation, or at a slightly different (heterodyne detection) frequency, while exciting the probe with light that is modulated in intensity at the same frequency, which usually corresponds to about the inverse of the luminescence lifetimes to measure (frequency do-main-based approach). The time-domain based approach is preferred. Either way, the temporal evolution of the luminescence intensity can be assessed.
In the method of the invention, transient triplet-state absorption, or, preferably, the temporal evolution of triplet-state based luminescence selected from the group consisting of delayed fluorescence or phosphorescence may be measured. Preferably, the temporal evolution of delayed fluorescence (DF) is measured if the probe belongs to the family of PpIX. DF has the advantage of being easily detectable in comparison to phosphorescence for this family of molecules, and is easier to implement than TTA, the latter requir-ing for example an additional light source.
If TTA is measured, two light sources are used. The pump beam excites the probe to the T1 state by emitting light at one of its absorption peaks detailed above. The probe beam, usually white light, is used to measure the evolution of the concentration of T1. Its intensity is continuously measured at one of the excited triplet states T1molecules absorption peaks. These wavelengths depend on the probe. For example, for PpIX, T1 absorption wavelength may be at 450 nm (Chantrell et al., 1976. Journal of Luminescence 12-13:887-891).
If triplet-state based luminescence is measured, the emission of said luminescence, e.g., delayed fluorescence in the case of the PpIX family, is measured. The emission is measured at emission wavelengths dependent on the probe. If the probe is such that, in certain condition of excitation fluence rate, probe and quencher concentration, it forms dimers out of two excited triplet states T1, so-called excimers, then both the luminescence of the probe monomer and excimer may be measured. For example, for delayed fluorescence of PpIX, emission may be measured at emission wavelengths provided in
In step c) of the method of the invention, said temporal evolution is correlated with said concentration of the probe and/or the quencher. The correlation is typically based on the Stern-Volmer relationship known in the art, however, in contrast to prior art methods, second order reactions (or processes) involving two excited triplet states of the probe are considered.
Several processes on which the lifetime of the excited triplet states, and consequently, e.g., delayed fluorescence, depends generally have been described in literature (but not yet in cells in tissues or organs) (Scholz and Dedic, 2016. Singlet-Oxygen-sensitized Delayed Fluorescence, 63-81; Parker and Hatchard, 1963, Transactions of the Faraday S0-ciety 59, 284; Parker, 1968. Photoluminescence of Solutions: With Applications to Photochemistry and Analytical Chemistry, Elsevier Publishing Group University of California, 2nd edn, p. 544). The necessary energy to perform RISC can come from different origins:
Surprisingly, in contrast to the belief held in the state of the art (e.g., EP1742038 and WO 2007/00487, Mik et al., Nature Methods 3: 939-945, 2006; Mik et al. Biophys. J. 95(8):3977-90, 2008; Mik et al., 2013. Anesthesia and Analegesia 117: 834-846; Vinklarek et al., 2017. Photochem. Photobiol. 16:507-518 with an explicit discussion of the situation in cells), the inventors could show that in cells, in particular, in skin, the temporal evolution of the triplet state, measurable e.g., by delayed fluorescence, involves a significant part of second order reactions.
In human skin, the assumption in the state of the art was that the main Reverse InterSystem Crossing (RISC) mechanism is E-type, thermal DF. According to this model, and using the notations introduced in
Thus, the DF lifetime is
and the DF lifetime in absence of oxygen
Note that, for PpIX, the contribution of phosphorescence (kPh) is negligible at body or room temperature. PpIX phosphorescence can be measured in non-physiological environments, according to Gouterman et al., 1974. Journal of Molecular Spectroscopy 53(1):88-100.
Consequently, the unexpected results mentioned above and described in more detail in the examples cannot be explained by the state of the art model of PpIX DF:
In light of these surprising observations, and in light of advances in the understanding of PpIX DF in vitro, the inventors provide a solution to these incoherencies by taking second order reactions involving two excited triplet states of the probe into account, in particular, taking second order desexcitation into account. As shown below, the inventors surprisingly found that formation of excimers and/or triplet-triplet annihilation (e.g., P-type DF) are second order processes that, if taken into account, render the method for determining the concentration of a probe and/or a quencher in a cell more precise.
The inventors have described new biological, chemical and physical principles of the intracellular pO2 measurement method based on the delayed fluorescence of Protoporphyrin IX, further developing the Triplet State Lifetime Technique (TSLT).
In particular, the discovery of the existence of PpIX excimers in vivo, in particular in cells in tissues or organs such as human skin, changes the understanding of DF dynamics. Excimers lead to a detectable DF signal at 670 nm. Furthermore, the inventors derived a new model of triplet population evolution and associated DF signals, when second-order process play a significant role. The framework introduced, namely the concepts of depopulation equation, and measurement/depopulation pathways, predicted non-monoexponential and non-multiexponential decays. These predictions were confirmed in vivo.
Based on these new elements, the inventors proposed several methods to improve the precision and reliability of pO2 measurements: In different aspects, new fitting procedures are provided, either in the time range considered (Decay Central Fitting method), or in the function fitted (Mixed-orders decay fitting), or in the parameters used (Adaptive Stern-Volmer relationship). The method based on the adaptive Stern-Volmer relationship is based on the specific detection of the DF signal at 630 and 670 nm. These methods are further explained below.
While the inventors' experimental observations rely on delayed fluorescence measurements, the same principles apply for methods of the invention based on phosphorescence and transient triplet absorption. Indeed, the improved pO2 measurement methods are ultimately based on the excited triplet state concentration evolution analysis: any experimental method able to access it is fit for this invention.
In a preferred embodiment, the invention provides a method for determining the concentration of oxygen and/or the concentration of a heme precursor, such as PpIX, in a cell, wherein the cell optionally is in a tissue or organ, comprising steps of
wherein second order reactions involving two excited triplet states of the heme precursor are considered.
In most embodiments of the invention, preferable, the concentration of the quencher is determined, in particular, the concentration of oxygen, in particular, if the status of a tissue, an organ or a patient is to be assessed.
However, in case photodynamic treatment, e.g., of cancer cells is desired, there also is a high interest in also determining the concentration of heme precursors such as PpIX in cells, so in that case, the method is preferably carried out to determine both the concentration of the probe, in particular, of heme precursors such as PpIX, and of the quencher, in particular oxygen. The method allows the user to determine whether the probe and quencher concentrations are in the correct range to perform photodynamic therapy. Additionally, it allows to choose the optimal excitation intensity given the probe and quencher concentration to reach the best photodynamic effect, for example the optimal level of singlet oxygen creation. If the PpIX concentration is high, a low light dose should be used. If the PpIX concentration is low, a high light does should be used. The administration of probe precursor, for example ALA, or derivatives thereof, in the case of PpIX, can also be adjusted depending on the probe concentration measured. Appropriate average PpIX concentration for photodynamic therapy ranges from 10 μM to 5 mM.
There are also forms of porphyria which result in a pathologic accumulation of PpIX in many organs, including the skin. In these cases, e.g., for diagnostic applications, the method of the invention may also be used to determine the concentration of heme precursors such as PpIX in cells.
The invention provides a method for determining the concentration of oxygen in a cell, wherein the cell optionally is in a tissue or organ, comprising steps of
wherein second order reactions involving two excited triplet states of the heme precursor are considered. Preferably, the temporal evolution of delayed fluorescence is determined.
Alternatively, the invention provides a method for determining the concentration of a heme precursor, such as those involved in the heme biosynthetic pathway leading to the production of PpIX, in a cell, wherein the cell optionally is in a tissue or organ, comprising steps of
wherein second order reactions involving two excited triplet states of the heme precursor are considered. Preferably, the temporal evolution of delayed fluorescence is determined.
In preferred embodiments of the invention, the concentration of the probe, e.g., a heme precursor such as PpIX is determined, and, optionally, then, the concentration of the quencher, e.g., oxygen, is determined. Thus, the method can serve to determine both the concentration of the probe and the quencher in the cell. In the context of the invention, concentration can also be measured as partial pressure.
As described herein, inventors have surprisingly found that second order reactions of the probe have a significant influence on the decay of the excited triplet state of the probe after excitation. Second order reactions are influenced by the concentration of excited triplet states of the probe. The concentration of excited triplet states of the probe also depends on the excitation intensity and the concentration of the probe.
Thus, in one aspect, the concentration of excited triplet states of the probe is considered for determining the concentration of the probe or the concentration of the quencher.
In a preferred method of the invention, the emission of delayed fluorescence is measured at at least two wavelengths (including two wavelength ranges, as described above), wherein these two wavelengths are suitable for distinguishing delayed fluorescence of probe monomers and of probe excimers. In a preferred embodiment, the probe is a heme precursor such as PpIX, and the first wavelength is in the range of 615-645 nm, optionally, about 630 nm and the second wavelength is in the range of 646 nm-700 nm, optionally, about 670 nm. The emission at the two wavelengths is preferably measured essentially simultaneously, i.e., without any undue delay, and preferably in the same sample volume. Preferably, it is thus measured by two detectors, however, one detector can also be used if the control unit is configured to measure at both wavelengths in close sequence.
Based on this, the adaptive Stern-Volmer relationship may be used for determining the concentration of the quencher. Keeping in mind the central results discovered by the inventors: first, that PpIX lifetime, particularly in absence of oxygen, strongly depends on the local PpIX concentration, in a Stern-Volmer like quenching relationship, and second, that the ratio of DF energy at 670 and 630 nm is linearly linked to the local PpIX concentration.
The Adaptive Stern-Volmer Relationship
wherein, preferably, the probe is a heme precursor such as PpIX and the ratio E670/E630 is used in step i).
In more details the method steps comprise:
calculating the DF energy at the monomer and excimer wavelengths: E630=∫0+∞IDF630(t)dt and E670=∫0+∞IDF670(t)dt
deducing the concentration of PpIX ([PpIX]) in the cells,
wherein the ratio E670/E630 is a linear function of the PpIX concentration: e.g., E670/E630 =a+b[PpIX] with a and b known parameters (see
determining the DF lifetime at 630 nm in absence of oxygen τ0630 from the concentration of PpIX, according to
fitting the best mono-exponential function to the decay
at 630 nm: by determining the most adapted values for I0, τ630 and S0
deducing the pO2 using the adaptive Stern-Volmer relationship (kqO
As an alternative to calculating the DF energy, it is possible to determine the initial DF intensity (l0=IDF(t=0)). Then, the ratio of initial intensities I0670/I0630 is a linear function of the PpIX concentration (see
Note: the “adaptive” nature of this relationship lies in the fact that it does not use a single value for τ0630, but it is adapted depending on the local concentration of PpIX. Determining (ratios of) initial intensities I0670/I0630 or intensities a short time (short compared to the lifetime of the luminescence) after the excitation pulse an be advantageous when compared with the integral intensity E because in first approximation they are not altered by the possible presence of triplet state quenchers, like oxygen. In practical applications, where oxygen is present, the DF integral intensity at 630 and 670 nm may be significantly reduced. Because of stronger prompt fluorescence or limited dynamic range of detectors it may not be possible to measure initial intensities of the delayed fluorescence directly. Initial intensities can be derived from observable parameters such as intensities and/or their decay a short term, e.g., between 50 ns and 50 μs, preferably, between 200 ns ad 20 μs, more preferably, between 1 μs and 10 μs, after the excitation pulse, or fitted functions for the entire decay, preferably correcting for the contributions of quenchers.
In one aspect of the invention, which can be combined with other approaches, e.g., with the approach underlying the adaptive Stern-Volmer relationship, the decay central fitting method or the mixed orders decay as explained in more detail below, the fluence rate (intensity) of the light for excitation is considered for determining the relative contribution of first and second order triplet deactivation and/or luminescence processes, or for determining the concentration of the quencher. If a low fluence rate is used, the excited triplet state is reached for a low fraction of the probe molecules. If a high fluence rate is used, the excited triplet state is reached for a high fraction of the probe molecules. The fluence rate thus directly affects the concentration of excited triplet state in the cells. First order triplet deactivation refers to interactions that involve one molecule of a probe substance in triplet state is involved that—usually via an intermediate step—ends up in its non-triplet ground state. Second order triplet deactivation refers to deactivations that involve two molecules in the first excited triplet state T1.
Thus, in one embodiment of the method of the invention, the probe is excited by irradiation with light having at least two different fluence rates (intensities) or light doses (if the pulse durations are shorter than the deexcitation times). Variable excitation allows for more precise determination of the contribution of first vs second-order measurement pathways e.g., as explained in detail in the examples below. Indeed, for a given excitation light dose, first order processes are proportional to the exciting fluence rate (intensities), while second-order processes are proportional to the square of the fluence rate (intensities). For instance, the relative contribution of E-type and P-type DF can be determined with this method.
Alternatively or additionally, the probe can also be excited by irradiation with light having at least two excitation wavelengths. Like the excitation fluence rate and dose, the wavelength also directly influences the fraction of probe molecules in an excited triplet state, and, thus, the concentration of the excited triplet state. For example, excitation at a wavelength corresponding to a maximum of the probe absorption peaks leads to more energy being transferred to the probe.
However, the excitation wavelength also has effects on the location of the volume of cells for which the determination is carried out, because different wavelengths have different penetration depths. Thus, the excitation wavelength can also be varied in order to analyse the concentration of the probe and/or the concentration of quencher at different penetration depths. The penetration depths of light of different wavelengths are known in the art (S. Jacques, “Optical Properties of Biological Tissues: A Review”, Physics in Medicine and Biology 58(11):R37-R61, 2013.). If this is of interest, the intensity of excitation can be adapted to make up for lower excitation.
In another aspect, which is optionally combined with excitation at variable excitation intensities or doses, the decay central fitting method (DCF) is used for correlating said temporal evolution with said concentration of quencher. In particular, the decay central fitting method rejects the initial non-exponential part of the decay. Preferably, the method further rejects the tail of the decay dominated by noise, thus selecting the central part of the decay of delayed fluorescence intensity over time. The central part of the decay of delayed fluorescence intensity over time comprises the time after t1, when the monoexponentiality indicator first gets below a threshold T, and, optionally, before t2, when the signal to noise ratio gets below a threshold T.
In details, the algorithm determines threshold values and time boundaries as follows:
Assuming that the DF intensity I (t) can be written as the sum of noise ΔI(t), centered around 0, and an average intensity Ī(t), which is typically the function obtained after fitting of the experimental signal, such that I(t)=Ī(t)+ΔI(t).
Defining the monoexponentiality indicator as
With ′ denoting the derivation with respect to time. The thresholds are values such that: M(t1)=0.7=T.
Introducing the signal to noise ration SNR(t), the noise dominated section starts, for example, at the time value such that
In detail, the algorithm comprises:
determining t1 and t2 as detailed herein
between t1 and t2, fitting a mono-exponential to
by determining the most adapted values for I0,τ630 and y0,
deducing pO2 using the conventional Stern-Volmer relationship (kqO
For example, under the conditions used in the COMET setup, the mono-exponential linear section comprises the time between 275 and 600 μs.
The DCF method may be used in combination with a variation of excitation intensities to quantify linear vs. quadratic behaviour, in particular, the contribution of P-type DF and excimer DF.
In another aspect of the invention, called “mixed-orders fitting”, the concentration of the quencher is determined on the basis of the assumption that the excited triplet state of the probe deactivates through a mix of first and second order reactions. The deactivation of the excited triplet state of the probe can be determined by measurements of delayed fluorescence, phosphorescence or TTA.
Preferably, the concentration of the quencher is determined on the basis of the assumption that the delayed fluorescence is caused by a mix of first and second order reactions.
Thus, the concentration of the quencher may be determined based on the formula:
I
630(t)=θS
With
wherein kE, k1C, kqO
Throughout the invention, if the probe is a heme precursor such as PpIX and the quencher is oxygen, before step a) of the method of the invention, a precursor of PpIX, preferably, 5-aminolevulinic acid (ALA), is administered to the cell(s) or has been administered to the cell(s). Instead of ALA, derivatives, in particular, prodrugs thereof may also be used, e.g., esters such as 5-ALA-OMe or 5-ALA-OHex. ALA prodrugs and derivatives are known in the art (e.g., Tewari et al., 2018. Photochemical & Photobiological Sciences 17:1553-1572; Wu et al., 2017. ACS Appl Mater Interfaces 9(17):14596-14605, Li et al., 2019. European Journal of Medicinal Chemistry 181:111582). This increases the level of PpIX in the cells, even though PpIX is also naturally present at low concentrations, as heme precursor. PpIX is typically produced in the mitochondria of cells. Alternatively, the level of PpIX can be increased by administering PpIX. Different methods for administering PpIX, or preferably, ALA and its derivatives are known in the art. The selection depends on the cells of interest, in particular, on the tissue or organ. For example, ALA can be administered to skin in a transdermal patch (e.g., available from Photonamics, Wedel, DE).
The invention thus also provides a precursor of protoporphyrin IX, preferably, 5-aminolevulinic acid or, alternatively, prodrugs thereof, for use of the invention in an in vivo method, wherein the probe is a heme precursor, such as PpIX, and the quencher is oxygen and the cell is in a human or animal subject, preferably, a human subject.
It is clinically very useful that the method of the invention can be used for assessment of mitochondrial function in cells. In particular, the invention provides a method for assessment of mitochondrial function in a sample comprising cells in a tissue or organ, comprising restricting or ceasing the supply of oxygen to said sample and carrying out the method of the invention with said sample, wherein the probe is PpIX and the quencher is oxygen. The method can be carried out as disclosed in WO 2010/01041 with the adaptations described herein, e.g., measuring at two emission wavelengths and using the adaptive Stern-Volmer relationship.
Throughout the invention, cell encompasses at least one cell. Typically, the concentration of probe and/or quencher is determined in a plurality of cells. The cells of interest are determined by the sample for which the emission of luminescence and/or TTA is assessed.
The cell typically is a living cell. It may also be a dying cell.
Preferably, the cell is part of a tissue or organ. The tissue may be part of an organ, or it may be an isolated tissue. Due to the easy accessibility of skin, the cell preferably is part of skin. However, e.g., under perioperative conditions, other tissues or organs are also accessible, e.g., a tumor. The organ may for example also be liver, kidney, brain, bladder or heart. The tissue may be tumor tissue.
Preferably, the cell is in situ in a subject, preferably, a patient, e.g., a human subject or patient. The cell may also be in a tissue or organ intended for transplantation, wherein preferably, the method provides information that allows for assessment of the viability of the tissue or organ. For example, a tissue or organ derived from a subject, e.g., a human subject may be transported, and before transplantation to the recipient, the status, e.g., the concentration or partial pressure of oxygen in the tissue or organ may be assessed. This allows for conclusions on the viability, and thus, depending on the assessment, the tissue or organ is transplanted or rejected.
The invention also provides a method for assessment of the status of a patient, e.g., a sepsis patient, a critically ill patient, a patient undergoing a tumor treatment (e.g., phototherapy or photodynamic therapy), a patient undergoing surgery, a patient suffering from a neurodegenerative condition, such as Alzheimer, Parkinson or Huntington diseases, or a decubitus patient, comprising carrying out the method of the invention, wherein the cell is a cell of said patient. Preferably, the patient is human.
The invention also provides a method of treating a patient, optionally, a human patient, e.g., a sepsis patient, a critically ill patient, a patient undergoing a tumor treatment (e.g., phototherapy or photodynamic therapy), a patient undergoing surgery, a patient suffering from a neurodegenerative condition, such as Alzheimer, Parkinson or Huntington diseases, or a decubitus patient, or for selecting an organ suitable for transplantation to a patient, comprising carrying out the method of any of the invention, wherein the cell is a cell of said patient (or, for selecting an organ, a cell of said organ).
For example, in the context of photodynamic therapy, ALA is applied, and converted to PpIX in a tumor. This can be advantageously combined with the method of the invention, e.g., to find out if a certain kind of treatment, e.g., irradiation or photodynamic treatment is suitable for treatment of the tumor. For instance, based on the oxygen concentration determined using the invention, the light dose to perform photodynamic therapy can be adapted to reach maximal efficiency. For instance if little oxygen is present in (parts of) the tumor, additional local irradiation is likely to improve the treatment effect. Oxygen concentration in a tumor and/or around a tumor site can also be monitored with a method of the invention, to monitor the progress of disease and/or therapy. The treatment can then be adapted to the progress.
In a further aspect, the invention provides a method for determining the temperature at a location of PpIX in a cell, preferably, in the mitochondria of a cell, comprising steps of
a) minimizing the concentration of a quencher, preferably, by restricting or ceasing the supply of oxygen to said cell,
b) exciting a probe capable of exhibiting delayed fluorescence by irradiation with light having a wavelength and a fluence rate suitable for excitation of the probe, wherein the probe preferably is a heme precursor such as PpIX,
c) determining the lifetime of the monoexponential component of the decay according to a method of any of embodiments 1-36, preferably, by the decay central fitting method, or by the mixed-orders fitting method, and
d) calculating the temperature through the Arrhenius relationship.
A method of the invention for determining the temperature in the mitochondria of a cell, may comprise steps of
wherein:
Ea is the energy gap between the first excited singlet and triplet states of PpIX (S1 and T1 respectively, typically Ea=35 kJ. mol−1; Suisalu, et al. Phosphorescence-detected triplet state microwave resonance of a protoporphyrin molecule in n-octane, Chem. Phys. Lett. 101 (1983) 182-186; Gouterman, et al. Porphyrin free base phosphorescence, J. Mol. Spectrosc. 53 (1974) 88-100.)
R=8.314 J.K−1.mol−1 the ideal gas constant
A: [s−1] the E-type reaction rate constant pre-exponential factor:
A temperature variation can be determined even more precisely through the following relationship, which does not require the knowledge of A:
In one aspect, the invention provides a device suitable for carrying out the method of the invention, comprising
The excitation light is provided by an excitation light source, which can be any light source capable of providing light of a, or several, suitable wavelength(s) in continuous, modulated or pulsed fashion, wherein, continuous light is combined with an accessory presenting variable transmission. Examples include a xenon light source with bandpass filter or monochromator, light emitting diodes (LEDs) and several types of laser systems (e.g., diode lasers or other solid state or gas lasers, tunable or not). Typically, the light source is a flash lamp, a LED or pulsed laser system.
The light detector may be any photodetector, such as for instance, a photodiode, ava-lanche photodiode, photomultiplier tube, charge coupled device (e.g., CCD camera). The light detector may comprise an image intensifier or may not comprise an image intensifier. Detection systems may comprise phase-locked detection techniques in order to improve the signal-to-noise ratio. The light detector may be arranged to detect light emission, e.g., luminescence, as it is functionally linked to an optical fiber capturing said light emission, e.g., luminescence. A system of lenses, filters and/or mirrors may also be used to capture the light, e.g., luminescence from the cell. If an optical fiber is used, it is preferably positioned with one end facing the measurement sample comprising the cell or cells of interest, in contact or not, and the other end towards the detector.
For detection of delayed fluorescence, the photodetector is typically gated with a tunable delay to exclude the strong prompt fluorescence and only capture the delayed fluorescence.
The control unit may be configured to obtain repetitive or continuous measurements of the luminescence or triplet triplet absorption. As stated above, the detected parameter preferably is delayed fluorescence, so the device is preferably suitable for measuring delayed fluorescence. Devices comprising these features are known in the state of the art, e.g., the COMET device (Photonics Healthcare).
Based on the modifications to prior art methods for determining more precise concentrations of probe and/or quencher, the device also needs to be modified by at least one feature. A device of the invention thus further comprises
For example, bandpass filters suitable for selecting delayed fluorescence at a wavelength (or wavelength range) suitable for distinguishing the delayed fluorescence of probe monomers and excimers, e.g., at about 630 nm and at about 670 nm for PpIX as the probe, may be employed. Alternatively, or in addition, the excitation light source, e.g., the excitation laser may have a variable excitation power and/or pulse duration. The pO2 calculation algorithm may also be changing to a calculation algorithm of the invention. A device with preferred features is presented in
A preferred device of the invention comprises means for determining the delayed fluorescence lifetimes and/or intensities of a probe simultaneously at two emission wavelengths suitable for distinguishing delayed fluorescence of probe monomers and of probe excimers, wherein the same sample is analysed. If the probe is a heme precursor such as PpIX, and delayed fluorescence is determined, the first detected wavelength is in the range of 615-645 nm, optionally, about 630 nm, and the second wavelength is in the range of 646 nm-700 nm, optionally, about 670 nm. For example, the means for determining the delayed fluorescence of PpIX at two emission wavelengths may be spectral filters, grating or spectrometers, preferably, spectral filters such as band-pass filters.
The device of the invention may comprise means for sequential or simultaneous light detection at at least two different wavelengths by a light detector, e.g., a photomultiplier, and a photomultiplier array. The device may, e.g., comprise means for simultaneous light detection.
The device of the invention may comprise a gated photomultiplier that is configured to reject the initial prompt fluorescence, and then opens to collect the delayed fluorescence, or two such gated photomultipliers.
Alternatively, or additionally, the device of the invention may comprise means for exciting the probe with at least two, preferably, five different excitation fluence rates and/or fluence, optionally, at least 10. The means for exciting PpIX with at least two, preferably, five different excitation fluence rates and/or fluence may be, e.g. filters suitable for placement in the excitation beam, a laser with variable power and energy per pulse, and/or a plurality of lasers with different excitation fluence rates, preferably, a laser with variable excitation power.
The device of the invention is useful for carrying out the method of the invention, e.g., it is suitable for assessing the concentration of a probe and/or quencher in a cell, typically in a tissue or organ, preferably, in situ in a human patient.
Preferably, the device of the invention further comprises a pressure pad for applying local pressure on tissue containing arterioles, veins and/or capillary bed that supply oxygen to the sample volume. This allows for assessment of the status of the mitochondria in said tissue or sample volume, as described herein.
The device of the invention may further comprise
iii) means for exciting the probe with at least two different excitation wavelengths. Such means may comprise spectral filters, grating, a light source, e.g., a laser, with variable excitation wavelength, or two or more lasers with different excitation wavelengths, or spectrometers, preferably, several lasers.
The device of the invention preferably further comprises a processing unit capable of processing the obtained measurements and configured to apply the adaptive Stern-Volmer relationship and/or the decay central fitting method and/or the mixed orders fitting method, preferably, at least the adaptive Stern-Volmer relationship. Alternatively, the device may comprise communication means adapted to communicate with an external processing unit capable of processing the obtained measurements and configured to apply the adaptive Stern-Volmer relationship and/or the decay central fitting method and/or the mixed orders fitting method, preferably, at least the adaptive Stern-Volmer relationship, e.g., a computer.
In one aspect, the invention provides the use of the device of the invention for determining the concentration of a quencher such as oxygen and/or of a probe such as PpIX in a cell according to the method of the invention.
In embodiment 1, the invention provides a method for determining the concentration of a quencher (preferably oxygen) and/or the concentration of a probe capable of exhibiting a triplet-state based luminescence or transient triplet absorption (preferably a heme precursor such as PpIX) in a cell, wherein the cell optionally is in a tissue or organ, comprising steps of
wherein second order reactions involving two excited triplet states of the probe are considered.
In embodiment 2, in the method of embodiment 1, the quencher is oxygen, and the probe is a heme precursor such as protoporphyrin IX (PpIX).
In embodiment 3, in the method of any of embodiments 1 or 2, the temporal evolution of delayed fluorescence of the probe is measured at at least one emission wavelength.
In embodiment 4, in the method of any of embodiments 1 or 2, the temporal evolution of phosphorescence of the probe is measured at at least one emission wavelength.
In embodiment 5, in the method of any of embodiments 1 or 2, the temporal evolution of transient triplet absorption of the probe is measured at at least one absorption wavelength.
In embodiment 6, in the method of any of embodiments 1-5, the concentration of the quencher is determined, preferably, the concentration of oxygen.
In embodiment 7, in the method of any of embodiments 1-6, the concentration of the probe is determined, preferably, the concentration of heme precursors such as PpIX.
In embodiment 8, in the method of any of embodiments 1-7, the concentration of the probe is determined, and, then, the concentration of the quencher is determined, preferably, the concentration of heme precursors such as PpIX and of oxygen.
In embodiment 9, in the method of any of embodiments 1-8, wherein the concentration of excited triplet states of the probe is considered for determining the concentration of the probe and/or the concentration of the quencher.
In embodiment 10, in the method of any of embodiments 1-9, the emission of the delayed fluorescence is measured at at least two emission wavelengths, wherein these two emission wavelengths are suitable for distinguishing delayed fluorescence of probe monomers and of probe excimers.
In embodiment 11, in the method of any of embodiments 1-10, the probe is a heme precursor such as PpIX and the first emission wavelength is in the range of 615-645 nm, optionally, about 630 nm and the second emission wavelength is in the range of 646 nm-700 nm, optionally, about 670 nm.
In embodiment 12, in the method of any of embodiments 1-11, the adaptive Stern-Volmer relationship is used for determining the concentration of the quencher.
In embodiment 13, in the method of any of embodiments 1-12, the method
wherein, preferably, the probe is a heme precursor such as PpIX and the ratio E670/E630 is used in step i).
In embodiment 14, in the method of any of embodiments 1-13, the fluence rate (intensity) of the light for excitation is considered for determining the concentration of the quencher. In embodiment 15, in the method of embodiment 14, the probe is excited by irradiation with light having at least two different fluence rates.
In embodiment 16, in the method of any of embodiments 1-7, the decay central fitting method is used for correlating said temporal evolution with said concentration of quencher, wherein the decay central fitting method rejects the initial non-exponential part of the decay.
In embodiment 17, in the method of embodiment 16, the method further rejects the tail of the decay dominated by noise, thus selecting the central part of the decay of delayed fluorescence intensity over time.
In embodiment 18, in the method of embodiment 17, the central part of the decay of delayed fluorescence intensity over time comprises the time after t1, when the mono-expo-nentiality indicator first gets below a threshold T, and, optionally, before t2, when the signal to noise ratio gets below a threshold T.
In embodiment 19, the method of any of embodiments 16-18 is used in combination with a variation of excitation intensities to quantify linear vs. quadratic behaviour, in particular, the contribution of P-type DF and excimer DF.
In embodiment 20, in the method of any of embodiments 1-7, the concentration of the quencher is determined on the basis of the assumption that the excited triplet state of the probe deactivates through a mix of first and second order reactions, wherein, preferably, said deactivation is determined by delayed fluorescence or phosphorescence measurements or measurements of TTA.
In embodiment 21, in the method of any of embodiments 1-7 or 20, the concentration of the quencher is determined on the basis of the assumption that the delayed fluorescence is caused by a mix of first and second order reactions.
In embodiment 22, in the method of any of embodiments 1-7 or 20-21, the concentration of the quencher is determined based on the formula:
I
630(t)=θS
wherein
In embodiment 23, in the method of any of embodiments 1-22, the probe is a heme precursor such as PpIX and it is excited by irradiation with light having a wavelength of about 515 nm.
In embodiment 24, in the method of any of embodiments 1-23, the probe is a heme precursor such as PpIX and it is excited by irradiation with light having a wavelength of about 405 nm, wherein, preferably, the signal to noise ratio is improved compared to a measurement according to embodiment 23.
In embodiment 25, in the method of any of embodiments 1-24, the probe is excited by irradiation with light having at least two different excitation wavelengths.
In embodiment 26, in the method of any of embodiments 1-25, the selection of the excitation wavelength allows for selection of the penetration depth of the excitation light and, consequently, for selection of the analysed cells or tissue layers.
In embodiment 27, in the method of any of embodiments 1-26, the probe is a heme precursor such as PpIX and the quencher is oxygen, and, before step a), a precursor of PpIX, preferably, 5-aminolevulinic acid or derivatives thereof, is administered to the cell or has been administered to the cell.
In embodiment 28, the invention provides a method for assessment of mitochondrial function in a sample comprising a tissue or organ, comprising restricting or ceasing the supply of oxygen to said sample and carrying out the method of any of embodiments 1-27 with said sample, wherein the probe is PpIX and the quencher is oxygen.
In embodiment 29, in the method of any of embodiments 1-28, the cell is a living cell. In embodiment 30, in the method of any of embodiments 1-29, the cell is part of a tissue.
In embodiment 31, in the method of any of embodiments 1-30, the cell is part of an organ.
In embodiment 32, in the method of any of embodiments 1-31, the cell is part of skin.
In embodiment 33, in the method of any of embodiments 1-32, the cell is in situ in a patient.
In embodiment 34, in the method of embodiment 1-33, the patient is human.
In embodiment 35, in the method of any of embodiments 1-34, the cell is in a tissue or organ intended for transplantation, wherein preferably, the method provides information that allows for assessment of the viability of the tissue or organ.
In embodiment 36, the invention provides a method for assessment of the status of a patient selected from the group comprising a sepsis patient, a critically ill patient, a patient undergoing a tumor treatment (e.g., phototherapy or photodynamic therapy), a patient undergoing surgery, a patient suffering from a neurodegenerative condition, such as Alzheimer, Parkinson or Huntington diseases, or a decubitus patient, or for assessment of the status of an organ potentially suitable for transplantation to a patient, comprising carrying out the method of any of embodiments 1-35, wherein the cell is a cell of said patient or organ.
In embodiment 37, the invention provides a method of treating a patient, optionally, a human patient, selected from the group comprising a sepsis patient, a critically ill patient, a patient undergoing a tumor treatment (e.g., phototherapy or photodynamic therapy), a patient undergoing surgery, a patient suffering from a neurodegenerative condition, such as Alzheimer, Parkinson or Huntington diseases, or a decubitus patient, or for selecting an organ suitable for transplantation to a patient, comprising carrying out the method of any of any of embodiments 1-36, wherein the cell is a cell of said patient or organ.
In embodiment 38, the patient in the method of any of embodiments 36 or 37 is a human patient.
In embodiment 39, the invention provides a method for determining the temperature at a location of PpIX in a cell, preferably, in the mitochondria in a cell, comprising steps of
a) minimizing the concentration of a quencher, preferably, by restricting or ceasing the supply of oxygen to said cell,
b) exciting a probe capable of exhibiting delayed fluorescence by irradiation with light having an excitation wavelength and a fluence rate suitable for excitation of the probe, wherein the probe preferably is a heme precursor such as PpIX,
c) determining the lifetime of the monoexponential component of the decay according to a method of any of embodiments 1-36, preferably, by the decay central fitting method, and
d) calculating the temperature through the Arrhenius relationship.
In embodiment 40, the invention provides a precursor of protoporphyrin IX, preferably, 5-aminolevulinic acid, for use in an in vivo method of any of embodiments 1-39, wherein the probe is a heme precursor such as PpIX and the quencher is oxygen and the cell is in a human or animal subject.
In embodiment 41, the invention provides a device suitable for carrying out the method of any of embodiments 1-39, comprising
In embodiment 42, the device of embodiment 41 comprises means for determining the delayed fluorescence of a probe at two emission wavelengths suitable for distinguishing delayed fluorescence of probe monomers and of probe excimers, wherein preferably the device is suitable for simultaneous time-resolved detection of monomer and excimer delayed fluorescence.
In embodiment 43, in the device of embodiments 42, the probe is a heme precursor such as PpIX, delayed fluorescence is determined, and the first wavelength is in the range of 615-645 nm, optionally, about 630 nm and the second wavelength is in the range of 646 nm-700 nm, optionally, about 670 nm.
In embodiment 44, in the device of any of embodiments 41-43, the means for determining the delayed fluorescence lifetimes and intensities of PpIX at two wavelengths are selected from the group comprising spectral filters, grating and spectrometers, preferably, spectral filters such as band-pass filters.
In embodiment 45, the device of any of embodiments 41-44 further comprises means for sequential or simultaneous light detection at at least two different wavelengths by a light detector, e.g., photomultipliers.
In embodiment 46, the device of embodiment 45 comprises means for simultaneous light detection.
In embodiment 47, in the device of any of embodiments 41-46 comprises a gated photomultiplier that is configured to reject the initial prompt fluorescence, and then opened to collect the delayed fluorescence.
In embodiment 48, the device of any of embodiments 41-47 comprises means for exciting the probe with at least two, preferably five different excitation fluence rates and/or fluences.
In embodiment 49, in the device of embodiment 48, the means for exciting PpIX with at least two, e.g., five different excitation fluence rates and/or fluences are selected from the group consisting of filters suitable for placement in the excitation beam, a laser with variable excitation power and energy per pulse, and a plurality of lasers with different excitation fluence rates and/or fluences.
In embodiment 50, the device of any of embodiments 41-49 further comprises a pressure pad for applying local pressure on tissue containing arterioles, veins and/or capillary bed that supply oxygen to the sample volume.
In embodiment 51, the device of any of embodiments 41-50 further comprises
In embodiment 52, in the device of embodiment 51, the means for exciting PpIX with at least two different excitation wavelengths comprise spectral filters, grating, several light sources, e.g., a laser, with variable excitation wavelengths, or two or more lasers with different excitation wavelengths, and spectrometers, preferably, spectral filters such as band-pass filters.
In embodiment 53, the device of any of embodiments 41-52 further comprises a processing unit capable of processing the obtained measurements and configured to apply the adaptive Stern-Volmer relationship and/or the decay central fitting method and/or the mixed orders fitting method, preferably, at least the adaptive Stern-Volmer relationship.
In embodiment 54, the invention provides the use of the device of any of embodiments 41-53 for determining the concentration of a quencher such as oxygen and/or of a probe such as PpIX in a cell according to the method of any of embodiments 1-35.
Legends
τ0 in s and [PpIX] in μM. We have identified a Stern-Volmer like self-quenching of PpIX of constant kqPPIX=36.8 μM−1s−1.
and T630 over time, on one specific patient. The anticorrelation is clear, strongly suggesting the presence of excimers in vivo. This shows that the local PpIX concentration induces a strong noise-like variation of the lifetime which is not related to physiology. Averaging such graphs for different volunteers is not relevant here: as the variation of the energy ratio/τ630 depends on the local PpIX concentration, it strongly depends on the probe positioning, and is relatively variable between patients. The variation of
reflects more local rpm inhomogeneities than average PpIX concentration trends, contrary to
on 4 other volunteers. The anticorrelation is visible in every graph. The data presented on the bottom left hand corner was obtained while pressing the skin (pO2=0). Continuous line; DF lifetime 625 nm, Dashed line
and the PpIX concentration, as measured in vitro. The curve is linear:
[PpIX] in μM, r2=0.94. This linearity is characteristic of excimer DF.
(defined in
(defined in
In Vitro Discovery of PpIX Excimers
Measurement & Methods
The inventors worked with a solution of Protoporphyrin IX (PpIX, Fluka, >99.5%) diluted in Dimethylformamide (DMF, Sigma Aldrich D4551, >99%). Concentration ranged from 1 to 400 μM. The solution was placed in a 1 cm wide quartz cuvette. Delayed fluorescence curves were acquired on a calibrated Horiba Fluorolog-3 spectrofluorometer, operated in front face mode. The cuvette was illuminated by a pulsed Xenon lamp (pulse duration FWHM: 1 μs, 10000 pulses per curve), filtered using a monochromator. The sample was excited at 405 nm, with a FWHM of 5 nm. The emission at different wavelengths was collected with a FWHM of 3 nm. The temperature of the setup was maintained at 32° C. (or piloted when temperature was the variable) using a Pelletier Carry PCB-150 thermoregulator. The pO2 was brought to zero by bubbling N2 in the cuvette at a flow rate of 3 L/min, using a Brick™ gas mixer. The zone of bubbling was not illuminated by the excitation beam, and the samples were constantly stirred during measurements. The excitation intensity was modified by placing neutral density filters in the excitation beam. Exponential decays were fitted using Origin V8 (OriginLab, Northampton, Mass.): in the expression of a mono-exponential decay
the initial tensity l0, the DF lifetime τ, and the offset y0 were determined using the built-in Levenberg-Mar-quardt algorithm. The delayed fluorescence signals were collected and fitted from 0,13 to 11 ms.
Results
The first and simplest experiment was to quantify the DF lifetime in absence of oxygen as a function of the concentration of PpIX. In theory,
should be independent of the concentration of PpIX ([PpIX]). The curve is presented on
The rate of PpIX catabolization in human cells:
was estimated by inventors as on average, −226 μM/h.
Conversely, by looking at the decrease of DF lifetime between 2 and 4-5 hours of PpIX application, the inventors could extract the average DF lifetime variation:
which gives a PpIX concentration rate of:
These rates are surprisingly high, suggesting that the concentration in cells could easily be above mM. Interestingly, other methods of oxygen sensing using porphyrin nanoparti-cles have reported unexplained quenching effects, which the inventors postulate to be based on the same mechanism (Kubat et al., 2017. ACS Applied Materials and Interfaces 9: 36229-36238).
A second question is the molecular origin of the self-quenching. By measuring the DF spectrum of different concentrations of PpIX solutions, it appeared that a peak at 670 nm was gradually developing as concentration increased, as presented on
It appears that, similarly to other aromatic molecules, PpIX excimers can be formed either by the encounter of two molecules in the first excited triplet state T1, or by the reaction of an excited singlet S1 with a triplet T1. The excited triplet state concentration can influence the relative contribution of different modes of its decay and hence affect the lifetime of the excited triplet state. For high excited triplet state concentrations, such as those expected after active cells have generated PpIX after administration of 5-ALA, it can be important to consider the excited triplet state concentration or parameters that serve as a proxy thereof as described in this application in order to accurately determine quencher (oxygen) concentrations from observed triplet state lifetimes. The initial triplet state concentration is determined by the excitation radiant exposure per pulse H times the molar extinction coefficient at excitation wavelength ε. At low excited triplet state concentration (HE[PpIX]<6,0 μJ. cm−3, E-type, thermal DF is dominant, and as the triplet concentration increases (HE[PpIX]>12 μJ. cm−3), P-type DF and excimer DF contribute predominantly.
The formation and desexcitation of excimers is presented on an updated Jablonski diagram (
Assessement of the Existence of PpIX Excimers In Vivo and Contribution of P-Type Delayed Fluorescence
Materials & Methods
“Dumbo” Setup
Plasters of 8 mg of ALA (Alacare 8 mg medicated plaster, Photonamic GmbH, Wedel) were applied on the skin of 20 healthy volunteers, in a mediothoracic position, on the sternum. The 2×2 cm2 plasters contain a suspension of 5-ALA HCl crystals suspended in a polymer matrix. The delayed fluorescence of PpIX was readily detectable from 3 to >15 h of ALA application, with a maximum intensity around 8 h. The outer layer of dead skin cells (stratum corneum) was removed using an abrasion pad, following the standard microdermabrasion method (Mik et al., 2013). Every hour, from 3 to 8 h of ALA application, the DF in the skin was quantified using a dedicated optical setup. The skin was excited by a series of 3 to 10 27 ns pulses at 515 nm (light doses 300 μJ/cm2. Energy per pulse: 15 μJ according to Innolas Service Test sheet Inno P1166, 18/02/2016. Laser Innolas Mod3 Yb-Yag, passive cooling, frequency doubling, passive Q-switch, fiber coupling (SMA)) via a frontal light diffuser providing a circular spot directly in contact with the skin (illuminated surface: 5 mm2). 1 ms separated each pulse, so that the system would have time to desexcite fully between each pulse.
The fluorescence 10 ns) and delayed fluorescence (5 to 500 μs) were collected by a fiber next to the light diffuser, and subsequently sent to a Photomultiplier (PMT) (Hamamatsu H11526-20-NF) and amplifier (Hamamatsu C6438-01). The PMT was gated to reject the prompt fluorescence. The decays were recorded by a data acquisition board (DAQ NI USB 6259, National Instruments) and processing unit. The emission light was filtered using a band-pass filter in the [600; 700] nm region, to reject excitation light and keep only the DF of PpIX. For certain measurements, the DF signal was additionally filtered to extract specific spectral bands, using bandpass filters at 625±12 nm, 670±12 nm or 700±12 nm. OD outside of bandpass: >4. When needed, the excitation intensity was varied by placing neutral density filters (OD E [0.15; 0.9]) in the excitation beam. In this work, DF energy EDF is understood as the total energy of the DF: if I(t) is the DF intensity as a function of time as measured by the data acquisition, then the intensity of the DF is EDF=∫0+∞IDF(t) dt. When changing the excitation intensity, the initial concentration of T1 states is linearly changed, as all processes leading to its formation are described by linear Einstein equation (absorption and Intersystem Crossing), and the illumination zone is constant in size. In short, [T1]0∝lex. As a consequence, when [T1]0 changes, IDF is likely to vary linearly if the DF mechanism involves one T1 molecule, or quadratically, if it involves two T1 molecules. Unless otherwise mentioned, the DF intensity was manually fitted to a decreasing mono-exponential IDF(t)=0e−t/τ+y0 using the built-in fitting algorithm of the Origin V8 software.
Results
Contrary to the in vitro situation, PpIX concentration is not known in vivo, which makes the analysis more complex: we are unable to compare the DF at different concentrations, we have to track other proofs of the existence of excimers. The inventors chose to inves-tigate the competition process that occurs between the 630 and 670 nm luminescence pathways.
Evolution of E670/E630 compared to τ630
To assess the existence of excimers, the inventors looked at the DF energy at 670 nm. But the monomer form also has a slight DF signal at this wavelength (
The inventors could have chosen to study the ratio of initial intensities, because it would not depend on the quenching by oxygen, which can be slightly varying between the measurement at 670 and 630 nm. But, in reality, the initial intensity is highly dependent on the second-order mechanisms described rigorously below, which do not follow the Stern Volmer equation. The inventors thus concluded that the DF energy (time integral of the intensity) is more robust for such an analysis.
If excimers did not exist, the ratio of energy at 670 and 630 nm would be constant over time. Showing its variations could be enough a proof. But to make it more compelling, the inventors compared the possible variations of this ratio with the DF lifetime at 630 nm τ630: if excimers exist, these two variables should be anti-correlated: When the PpIX concentration is high, τ630 small. At the same time, excimers form and produce DF at 670 nm, which causes E670/E630 to be high. Conversely, when the PpIX concentration is reduced, τ630 is longer, and no excimers form so that E670/E630 is small.
The result is shown in
Nevertheless, to prove that this phenomenon was common in our measurements, we add the corresponding graphs of 4 other volunteers on
The result is presented in
These measurements all speak in favour of the existence of excimers in vivo. The inconsistencies detailed above could potentially have found an answer. However, the fact that excimer formation, a second-order process, plays a significant role in vivo in the tested conditions, opens an unforeseen question: is it sure that P-type DF, a second order DF process (2T1→S1+S0), does not occur? Could one bimolecular process appear and another one, starting with the same reactants, be negligible?
Contribution of P-type DF
The easiest method to discriminate between E and P-type DF is to study the dependence of DF energy on excitation intensity: in the case of a monomolecular process (E-type), this dependency should be linear, while in the case of a bimolecular one (P-type), it should be quadratic. Indeed, the initial concentration of triplet state [T1]0 is proportional to the excitation intensity, as the illuminated volume is constant. By fitting a second order polynome
E
DF
=aI
ex
2
bI
ex
+c
one can compare the linear contribution bIex to the quadratic one aI2ex. An example of such curve is presented in
Several conclusions can be drawn from these experiments:
Influence of Second Order Reactions on DF Processes: The Kinetics of Non-MonoexPonential Decays
Evolution of the Triplet State: The Depopulation Equation
Reminder on Kinetics
Considering the reaction
we define the reaction rate constant:
under the hypothesis of the law of mass action (reaction rate proportional to the product of the concentration of all reactants at the power of their stoichiometric coefficients), the reaction rate is
r=k[A] (6)
with k the reaction rate constant, which, as a first approximation, depends minimally on concentrations of reactants. In our case of excited states dynamics, the situation is iden-tical: reaction rates are also proportional to the reactants concentration, following the linear Einstein equations. When a reactant appears with a coefficient 1, in the equation, its concentration can be computed by solving the linear differential equation:
The result is an exponential function:
[A](t)=[A]0e−kt=[A]0e−t/τ (8)
which allows the introduction of a characteristic time τ=1/k, which does not depend on [A]. It is called the lifetime of A. To be completely clear, the population evolves as a mono-exponential if and only if it is only depopulated by first-order reactions. A mono-exponential kinetics is equivalent to a first-order depopulation. But as soon as two or more coefficients appears in front of A, we are in the non-linear case:
This case follows the kinetics:
where the concentrations evolution is harder to determine, and is not mono-exponential. When a reactant appears with a coefficient one, we will indifferently designate the process as linear, monomolecular or first-order. Conversely, when the stoichiometric coefficient is 2, we use the terms non linear, bimolecular or second-order.
The Traditional Stern-Volmer Model
In the traditional model, the assumption is that only E-type DF, Internal Conversion and oxygen quenching contribute to the triplet state desexcitation:
Taking Second-Order Desexcitation in Account
The inventors have shown above that the contribution of P-type DF was not negligible. They have also identified a potentially strong effect of excimer formation, which is also bimolecular. Based on this changed technical understanding of the processes involved, they include the quadratic desexcitation terms in the depopulation equation, which thus loses its linear nature (see definition of reaction rates on
This solves into a more complex function. As it involves both first and second-order reactions, it is designated the mixed-orders depopulation function:
In this expression, one can spot the exponential term
perturbed by the second-order constants kp+kd. We can study approximate forms of the function at short and long times:
Mixed-Orders Depopulation at Short Times
which yields, [T1](0)=[T1]0, but also
And if we assume that [T1]0 is big enough, this even reduces to
Thus, the beginning of the depopulation is completely controlled by the second-order processes, and its slope reflects the reaction rate constants kp and kd and the initially excited population [T1]0.
Mixed-Orders Depopulation at Long Times
When t→+∞
a mono-exponential decay of lifetime follows. exactly like the Stern-Volmer theory:
This shows that the later part of the desexcitation follows a Stern-Volmer quenching by oxygen.
However, the DF signal IDF(t) is only an indirect image of the triplet population [T1](t): to draw conclusions on how to analyse the decays, the inventors investigated this link.
From Triplet Concentration Evolution to DF Signals
The DF intensity produced by E-type DF is worth:
I
DF
E(t)=θS
with
the fluorescence quantum yield of the S1 state. Conversely the
DF intensity produced by P-type DF is:
I
DF
P(t)=θS
and finally the excimer DF intensity:
I
DF
excimer(t)=θD·(kTd[T1]2(t)+[S1](t)[S0]) (21)
with
with
the quantum yield of the excimer state. The expected DF intensity can be deduced by injecting the expression of [T1] of equation (15) into equations (19) to (21).
Pure Second-Order Desexcitation and DF
The inventors have further explored the result of a process where second-order reactions are so much stronger than first order ones that the depopulation equation becomes:
In that case, the solution is
and finally, it is logical to consider that only P-type will produce DF:
To assess the existence of P-type or excimer DF in vivo, the inventors compared the quality of a mono-exponential fit versus a mixed-orders decay fit on the DF signals collected. This is done in the next section.
Results
The inventors assessed the nature of the DF process in vivo, at standard excitation intensity. The central hypothesis is that the depopulation equation is of mixed orders, and the DF is caused by a mix of E-type and P-type. Thus, the inventors' proposition for the most exact fit is:
I
DF
630(t)=αIDFE+βIDFp=αθS
with [T1](t) from equation (15). Unfortunately, this function requires sophisticated fitting methods.
To allow for fitting with more standard type fitting methods, e.g., as implemented in Origin, the inventors compared the fitting quality of 4 distinct functions: Mono-exponential (1st order depopulation and E-type DF), Mixed-orders+E-type DF, mixed-orders+P-type DF, and 2nd order+P-type DF (equation (24)). In essence, the inventors are assessing the nature of two reactions:
The combination of these two reactions defines the expected kinetics of the DF. The different possibilities and associated kinetics are summarized in Table 2.
To determine which reactions are the measurement and depopulation pathways, the inventors consider standard 630 nm DF signals obtained on volunteers' sternum after 3 to 7 hours of ALA application. The inventors compared the statistical quality of different fits in Table 3.
These measurements further confirm that P-type and excimer DF are so strong that they significantly change the shape of the decays. Moreover, the great success of second-order+P-type fitting function (introduced above) suggests that most triplet state molecules desexcite via P-type RISC and DF. This does not mean that oxygen quenching is negligible: as described in equation 18, it controls the slope of the decay at longer times. Even if it has a moderate influence on the whole decay (and thus the global shape mainly depends on second-order processes), it has a strong influence on the lifetime that one can extract from the signal. It should also be considered that multi-parameters non-linear fittings are challenging: it is possible that, mixed-orders+P-type DF should the best fit.
Conclusion
In this section, the inventors have shown the necessity to take second order processes into account in order to correctly interpret the DF kinetics, and have exhibited the appropriate decay function. This allows for the development of methods that take these effects
8.46 × 10−10
into account and measure the pO2 more precisely. The correction can be done at different levels of accuracy and practical feasibility, described below.
Practical Consequences on Skin Mitochondrial Oxygen Measurement
Beside the new understanding of PpIX DF, the objective of the invention remains to improve the quality of cellular oxygen measurements, in particular, in patients. More specifically, the objective is to solve the inconsistencies described above. The inventors describe several methods in this section that solve this problem.
1. The Decay Central Fitting (DCF) Method to Capture the Linear Quenching Processes
As detailed above, a mixed-order decay (equation (15)) can be temporally subdivided in three distinctive parts:
This subdivision is illustrated in
The determination of t1, the time of transition between the initial non-exponential phase and the monoexponential phase, and t2, the time between the monoexponential phase and the noise-dominated section, works as follows: assuming that the DF intensity I(t) can be written as the sum of noise ΔI(t), centered around 0, and an average intensity Ī(t), which is typically the function obtained after fitting of the experimental signal, such that I(t)=Ī(t)+ΔI(t).
Defining the monoexponentiality indicator as
With ′ denoting the derivation with respect to time. The thresholds are values such that: M(t1)=0.7=T.
Introducing the signal to noise ration SNR(t), the noise dominated section starts, for example, at the time value such that
The inventors set the optimal values for the decay central fitting method to [t1; t2]=[275; 600]μs in the experimental conditions used (i.e., in the Dumbo setup) based on the criteria and considerations explained above.
A last point should be discussed: the lifetime indicated by the plateau (long 225 μs in this experiment on non-compressed skin),
is in reality not the triplet lifetime
When produced by E-type DF, the
DF intensity in the central range is indeed proportional to
but when produced by P-type DF, which represents 40 to 70% of the signal, the decay lifetime is worth half of the triplet lifetime due to the quadratic dependence. Consequently, the observed lifetime is an average between τT and τT/2. The weights of this average cannot be simply measured by this method, and depends on the local PpIX concentration and fluence rate. The variable excitation fluence rate can be used to determine the relative contributions of E-type and P-type DF, and thus calculate TT.
To assess the efficiency of this method, the data of
The lifetime calculated is less erratic than before application of the DCF method, but some variability remains as τ varies from 100 to 350 μs. This could be due to the variable contribution of E and P-type giving a random fluctuation between τT and τT/2. It suggests that the DCF method can still be rendered more precise with a quantitative discrimination between mono and bimolecular processes.
Conclusion
The decay central fitting method constitutes a step-forward towards a stable, precise measurement of the pO2. It is easy to apply, as it only requires to place a 630 nm bandpass filter in the emission beam and to change the mono-exponential fitting range to the time between [t1; t2], e.g., [275; 600] μs.
As its precision is limited by the unknown E-type vs P-type DF ratio, it is preferred to simultaneously vary the excitation intensity, quantify the fraction of linear vs quadratic behaviour (as done in the section “contribution of P-type DF” above), and subsequently deduce a more precise Ty and pO2.
2. Correction of PpIX Self-Quenching: The Adaptive Stern-Volmer Relationship
Another option, which would be conceptually less different from the standard model, but practically more complex, is to conserve the Stern-Volmer model, but to adapt the DF lifetime in absence of oxygen at 630 nm τ0630, in order to compensate for the variable self-quenching of PpIX. Keeping in mind the central results discovered by the inventors: first, that PpIX lifetime in absence of oxygen strongly depends on the local PpIX concentration, in a Stern-Volmer like quenching relationship, and second, that the ratio of DF energy at 670 and 630 nm is linearly linked to the local PpIX concentration. The method comprises:
measuring the DF signal at 630 nm and 670 nm simultaneously and in the same sample volume,
using the ratio of intensities I0670/I0630 and/or related parameters such as energies
or the inverse of such ratios or the ratio's of intensities at a certain time after the excitation pulse, or energies in a certain time period with respect to an excitation pulse (for example 1 μs or between 1 μs and 100 μs after the excitation pulse had ended or fallen back to below 1/10.000th of its peak) and/or with intensities or energies determined in a part of the spectrum including one of the wavelengths 670 nm or 630 nm, but not including wavelengths within more than 1 nm, preferably more than 5 nm, more preferably more than 15 nm, from the other and setting it in a relationship to an intensity or the energy observed in another part of the spectrum or in the entire spectrum, or from determining the phase shift with respect to a modulated excitation source of luminescence in such spectra, or the phase shift of relationships of spectra or the phase shift of the transient absorption of a probing wavelength to deduce the local concentration of PpIX,
deduce τ0630 out of [PpIX] according to the Stern-Volmer quenching relationship established earlier.
injecting τ0630 and τ630 into the adaptative Stern-Volmer relationship, and finally calculating the corrected pO2:
As shown in
An application of this method is provided in
Conclusion
This method has shown very good results: on volunteers whose pO2 had no reason to vary abruptly, the adaptive Stern-Volmer relationship flattened the pO2 around a stable value, which is physiologically more realistic. The method was very robust and yielded good results even on noisy data.
In terms of practical modifications, the method requires the selection of 2 emission wavelengths or spectral bands at about 630 and about 670 nm, thus, e.g., bandpass filters may be installed. The algorithmic part remains simple, as the fitting function is a mono-exponential. The inventors present the adaptative Stern-Volmer relationship as a most satisfactory method, in terms of precision and ease of implementation.
3. The Complete Method: Fitting Mixed-Orders Decays
As detailed above, when taking into account all different reactions, the DF intensity at 630 nm follows:
I
630(t)=θS
with
using notations introduced earlier.
The most precise method would be to fit this function on each decay. Once the different fixed parameters have been determined (OS, kE, kp,k1C, kqO
To reinforce the result, it is also possible to fit the decay at 670 nm, which is worth
I
DF
670(t)=θD·(kTd[T1]2(t)+kSd[S1](t)[S0]) (29)
with
Also, once typical values of [T1]0 and pO2 are known under the experimental conditions used, it is possible to start the optimization algorithm close to these values. A tailor-made algorithm can then be provided that reaches optimal parameters in a short time.
The application of this method in practice is complex: fitting of such a complex nonlinear function is a numerical challenge.
Exclusion of an In Vivo Contribution of Singlet Oxygen Feedback Delayed Fluorescence (SOFDF)
An additional mechanism of PpIX-DF has been recently described in artificial systems, e.g., by Vinklarek et al., 2017. It involves a double action of oxygen: a ground state triplet molecule collects the PpIX triplet energy, as in conventional quenching, but then this singlet oxygen excites another T1 molecule to S1, which subsequently produces SOFDF. The process can be summed up by the following reaction:
This process has been observed in vitro (Vinklarek et al., 2017; Scholz et al., 2013. Photochemical & Photobiological Sciences 12: 1873) and in isolated fibroblasts (Scholz et al., 2017. Photochemical & Photobiological Sciences 16: 1643-1653; Scholz et al., 2015. Photochemical & Photobiological Sciences 14: 700-713), with different photosensitizers, including PpIX. This type of DF is characterized by a rise-decay kinetics, described by
I
SOFDF(t)=C(e−t/τ1−e−t/τ2) (31)
with
the decay time and τ2=τT/2 the rise time (τT: triplet state lifetime and τΔ: singlet oxygen lifetime).
Despite careful observations, the inventors were not able to identify the rise decay kinetics in vivo, suggesting that SOFDF did not significantly occur in human skin. This can be explained by two reasons:
Compatibility with Previous Measurements (Oxygen Consumption)
A natural question after the discovery of PpIX excimers is: how measurements and conclusions in the past have been performed, using an incomplete model of triplet-state based optical properties, such as delayed fluorescence of probes, e.g., PpIX? The case of Oxygen Consumption Measurements (OCM) is typical: as described in the literature described above, this method, based on the compression of the microcirculation to assess the rate of oxygen consumption by skin cells, has shown excellent results in recent years.
How can the results have been so coherent, if the fitting method was incorrect? To evaluate it, the inventors compared the pO2 consumption curve obtained 1) in the physiologic case, according to the models of oxygen consumption as provided by the present invention, and 2) if this physiologic pO2 was deduced using the standard Stern-Volmer method, while in reality the PpIX decay follows the mixed-orders+P-type function (equation (15)):
According to simple models of oxygen consumption by tissue, in case of microcirculation blockage by compression, the pO2 evolution can be simplified to:
pO
2
real(t)=1+e−t (33)
where we arbitrarily set the time constant to 1. Assuming that the pO2 varies as pO2real, the DF signal that would be obtained if the DF process was a mixed-orders+P-type DF (equation (15)) was deduce, as suggested by table 3. This gives us a series of signals DFreal(t) for each value of pO2real(t). Then mono-exponential decays were incorrectly fitted to the different DFreal(t), thus obtaining an incorrect pO2incorrect(t) thanks to the standard Stern-Volmer relationship. Finally, pO2real was compared to pO2incorrect(t) to see how different the pO2 becomes when one uses an incomplete model of PpIX DF.
The result is presented in
As a conclusion, it is perfectly possible that a large number of studies were conducted based on an improper model of PpIX model, giving sometimes satisfactory results: it is only when looking in detail into the DF shape and lifetime variations that the inconsistencies can be detected, and then overcome by the method of the invention.
| Number | Date | Country | Kind |
|---|---|---|---|
| 20176067.5 | May 2020 | EP | regional |
| Filing Document | Filing Date | Country | Kind |
|---|---|---|---|
| PCT/EP2021/063642 | 5/21/2021 | WO |
