EP3698124A1

EP3698124A1 - Method of stimulated emission depletion microscophy having high spatial resolution

Info

Publication number: EP3698124A1
Application number: EP18815011.4A
Authority: EP
Inventors: Luca LANZANO'; Paolo Bianchini; Giuseppe Vicidomini; Alberto DIASPRO
Original assignee: Fondazione Istituto Italiano di Tecnologia
Current assignee: Fondazione Istituto Italiano di Tecnologia
Priority date: 2017-10-19
Filing date: 2018-10-19
Publication date: 2020-08-26
Anticipated expiration: 2038-10-19
Also published as: WO2019077556A1; EP3698124B1; IT201700118432A1; US11372224B2; US20200333573A1

Abstract

Method for increasing the optical resolution of a stimulated emission depletion microscope, or STED microscope (Stimulated Emission Depletion), based on the modulation of the intensity of a STED beam on an arbitrary time scale during the acquisition of an image and the analysis of the induced dynamics, without increasing the intensity of the STED beam and in a simple and economic manner.

Description

METHOD OF STI MULATED EMISSION DEPLETION MICROSCOPHY HAVING H IGH SPATIAL

RESOLUTION

* * *

The present invention relates to a method for increasing the optical resolution of a stimulated emission depletion microscope, or STED microscope (Stimulated Emission Depletion) based on the modulation of the intensity of a STED beam on an arbitrary time scale during the acquisition of an image and the analysis of the induced dynamics, without increasing the intensity of the STED beam and in a simple a nd economic manner. Such method is based on the principle that the fluorescence at the centre of a n observation volume is not modulated, while the fluorescence at the periphery of said volume is modulated. The exploitation of this difference in the modulation of the fluorescence signal entails an increase in the spatial resolution.

STED microscopy is one of the many different types of super-resolution microscopy techniques which have been recently developed to increase spatial resolution through exceeding the limit of diffraction of conventional light microscopy.

The main application of STED microscopy is directed to the study of sub-cellular architectures and dynamics and takes advantage of the non-linear response of fluorophores commonly used for marking biological samples. One sample is excited with a first radiation beam, the excitation beam, and the fluorescence is measured in a conventional manner. One donut-shaped second radiation beam, the STED beam, de-energizes by stimulated emission the fluorophores in an annular region around the axis of the first beam, allowing receiving the signal only from the central region of the first beam. The mechanism responsible for the increase in resolution is the saturation of the fluorescence reduction by stimulated emission. With STED microscopy high spatial resolutions were reached.

The spatial resolution w of a STED microscope depends on the intensity ISTED of the STED beam according to the equation :

1 _ 1

W I sat Eq. l

wherein wo is the spatial resolution of the confocal microscope used in STED microscopy and Isat is a constant representing the intensity value of the STED beam, required to increase the resolution of a factor Λ/2 that depends on the sa mple. In theory, a n unlimited spatial resolution could be obtained by increasing the intensity value ISTED of the STED beam. I n most practical cases, the maximum resolution of a STED microscope is limited by the maximum power "tolerated" by the sample due to detrimental photo-damage effects induced by illumination with the STED beam, such as for example phototoxicity and photobleaching. For this reason strategies were developed to increase the spatial resolution without increasing the intensity of the STEAD beam. For example, document WO2015/022635 Al discloses a method that uses a pulsed excitation beam (in the picoseconds time scale) and a continuous wave, by varying the detection time window. Since the average life of the fluorescent decay after a pulsed excitation is always in the time scale of nanoseconds, such method involves the use of an ultra fast detection electronics in the time scale of nanoseconds. In document L.V. Doronina-Amitonova et a I, 'Ultrahigh- contrast imaging by temporally modulated stimulated emission depletion', Optics Letters 40,725-728 (2015), STED light beams are used with a Gaussian time-modulated spatial pattern, whereby the whole fluorescence signal is modulated in order to distinguish it from a background signal, thereby increasing the contrast of fluorescence images, but without increasing the resolution thereof. Also document Chaoyang Fan et al, 'All-optical fluorescence image recovery using modulated stimulated emission Depletion', Chemical Science 2, 1080- 1085 (2011), uses STED light beams with a Gaussian time-modulated spatial pattern, only modestly increasing the resolution. In the document by Lanzano et al, "Encoding and decoding spatio-temporal information for super-resolution microscopy". Nat Commun 6, 6701 (2015), a method of analysis called SPLIT (Separation of photons by Lifetime Tuning) is described, which improves the spatial resolution by recording the temporal decay of the fluorescence signal and using the additional information encoded in the fluorescence decay dynamics. This is obtained by separating one super-resolved component, characterized by a specific temporal footprint. In fact, the key idea on which the SPLIT analysis method is based is that for each pixel the fluorescence intensity decay can be described as the linear superposition of one slow decay component (due to fluorophores in the centre of the actual observation volume) and one fast decay component (due to fluorophores at the periphery of the actual observation volume). By using an algorithm based on Fourier transform of decays, the fraction associated with the slow component can be extracted effectively without using minimization algorithms, with consequent increase in the spatial resolution. Even the SPLIT analysis method requires the use of a pulsed excitation and ultra fast detection electronics, in order to record the fluorescence dynamics in the time scale of nanoseconds.

The object of the present invention is to overcome the disadvantages described heretofore, allowing to increase the optical resolution of a STED microscope in a reliable, simple and economical way, without increasing the intensity of the STED beam and without work in the time scale of nanoseconds.

It is therefore one specific object of the present invention a method for increasing the optical resolution of a microscope STED based on the modulation of an intensity of a STED beam during the acquisition of an image as defined by independent claim 1.

Further embodiments of the method according to the invention are defined in the dependent claims.

The advantages offered by the method according to the invention with respect to the solutions of the prior art are numerous and important.

First of all, the method uses one STED light beam having a donut-shaped spatial pattern (in English doughnut or donut) which has the advantage of acting only on the periphery of the observation volume of the microscope. Thanks to the "donut" shape, since only the signal at the periphery is modulated, the latter can be subtracted from the whole signal and a smaller observation volume can be obtained. In other words, the method advantageously allows increasing the spatial resolution without increasing the intensity of the STED beam, thus avoiding damaging the sample. Moreover, the method of the present invention based on the slow modulation of the intensity of the STED beam and the analysis of the dynamics induced entails a general simplification of procedure and can be integrated into any STED microscope. The fluorescence modulation at the periphery of the observation volume is performed on an arbitrary time scale, whereby it is not necessary to use a pulsed excitation beam. This entails the advantage of being able to use light sources such as for example laser, which are inexpensive compared to those required for a pulsed excitation of other techniques. Another important advantage is that it is not necessary to detect the fluorescence dynamics in the time scale of nanoseconds and therefore a simple and economical electronic detection can be used.

The present invention will be now described, for illustrative but not limiting purposes, according to its preferred embodiments, with particular reference to the Figures of the accompanying drawings, wherein:

Figure 1 is a diagram of the principles on which standard STED microscopy (la) and STED microscopy integrated with the present method (lb) are based;

Figure 2 is a flow diagram of a first embodiment of the method for increasing the optical resolution of a STED microscope, according to the invention; Figure 3 schematically shows some of the results of a simulated first experiment, obtained by performing the first embodiment of the method of Figure 2.

Figure 4 is a flow diagram of a second embodiment of the method for increasing the optical resolution of a STED microscope, according to the invention;

Figure 5 schematically shows some of the results of a simulated second experiment, obtained by performing the second embodiment of the method of Figure 4;

Figure 6 schematically shows some of the experimental results, obtained by performing the first embodiment of the method of Figure 2;

Figure 7 illustrates the results of processing performed on simulated data in the first experiment of Figure 3.

In the Figures the same reference numerals will be used for similar elements.

With reference to Figure la, the principle on which STED microscopy is based is that excited fluorophores at the periphery of a confocal observation volume are "turned off" in a selective manner by means of a (donut-shaped) second beam, the STED beam. Since the fluorophores at the centre are not turned off, a smaller actual STED observation volume is obtained and an increase in the spatial resolution. The maximum resolution of a STED microscope is determined by the ability to turn off the molecules at the periphery of the observation volume, and this is generally determined by the intensity ISTED of the STED beam. The method of the present invention provides for a significant increase of the resolution by modulation of the STED beam intensity but without increasing the average intensity of the STED beam, as shown schematically in Figure lb. The STED beam modulation generates spatial variations in the temporal dynamic within the actua l STED observation volume. The analysis of the signal temporal evolution allows to separate smaller actual STED-modulated observation volume characterized by a specific temporal footprint.

Fig. 2 is a flow diagram of a first embodiment of the method for increasing the optical resolution of a STED microscope. Advantageously, the method is performed automatically. I n an initial step 100 a sam ple is illuminated with either a pulsed or continuous excitation light beam with constant excitation intensity, in order to excite fluorophores in the sample. I n a step 200 the sa mple is illuminated with a depletion light beam, either pulsed or continuous, STED beam, modulated to de-energize fluorophores by stim ulated emission. I .e. the intensity ISTED of the STED beam is modulated over time, ISTED = IsTEo(t), by generating spatia l variations in the temporal dynamic inside an actual observation volume of the sample. On an image plane xy of pixels, the image is given by an image function F=F(x,y) which is expressed as the convolution of a distribution of a sample object p (χ', y') a nd a point spread function, (PSF, Point Spread Function) of the microscope. It should be noted that in general the distribution of a sample object is three-dimensional, but since the preferred embodiments of the method described later have an effect only on the xy plane, the distribution is approximated on the xy plane. It should be noted that it is possible to extend the method also to the third spatial axis to increase also the lateral resolution of the image. By assuming a PSF of Gaussian shape, the convolution equation is: herein dS is an infinitesimal surface element a nd w is the spatia l resolution expressed by Eq whereby the Eq. 2 can be rewritten as:

^FSTED

with FSTED an image function obtained with STED microscopy. It should be noted that the Eq. 3 and 4 represent the integral of the contributions of all fluorophores to the fluorescence signal. The intensity modulation of the STED beam performed at step 200 also causes the image function of the Eq. 3 to be a function of time:

t- Eq. 4, wherein the subscript "t-STED" indicates an image resolved in time, obtained through modulation of the intensity of the STED beam. The modulation of the image function F is then obtained by the modulation of IsTEo(t) for each pixel.

In other words, if in a given pixel (x, y), the intensity ISTED = lsTEo(x, y, t) of the STED beam is modulated over time, a modulation M(x,y) of the value of the image function in that pixel, F = F (x, y, t) will result. Modulation M(x,y) of an image is defined as the quantity:

I ')e hat J m

The term m(x,y) contains sub-diffractive information on the distribution of the object p(x, y) in the vicinity of the position (x,y). The multiplication factor 5(lsTEo(t)/lsat) is a function of the modulation of the intensity of IsTEo(t) of the STED beam.

In the first preferred embodiment of the method of figure 2, the modulation M(x,y) is calculated from the signal F_t-sTED(x,y,t) through an algorithm based on Fourier transforms. For each pixel, the amounts g(x,y) and s(x,y) according to the following formulas are first calculated:

N,

∑^FtsTED(^x> y, t)∞i2 (j - l)/ N_t)

N,

∑^Ft-STED(^X> y> t) ^27T (^j - ^Nt )

wherein N_t is the total number of points in time wherein the signal F_t-sTED(x,y,t) is sampled. The modulation M(x,y) is then calculated in each pixel, as:

M{x, y) = ^(g(x, y))² + (s(x, y) Eq.7

By the term M(x,y), one fraction _Cenier(x,y) of intensity of the fluorescence signal, corresponding to the contribution to the fluorescence signal coming from the fluorophores at the centre of the observation volume, is calculated. In fact the modulation M(x,y) in a given pixel can be expressed as a combination of the modulation value at the centre of the observation volume M_center and the modulation value in the periphery of the observation

Volume Mperiph, according tO the formula: (1 - fcenterfcy)) Mperiph, wherein M_center and Mperiph are two constant having non negative value. From this formula it is obtained:

(x, y) = - y) - M_CENTER center Eq.8

The values of the consta nts Mcente and Mpe ph that allow to obtain a significant increase in resolution without the introduction of artefacts in the image, must be determined in the specific embodiment of the method, since they depend on some experimental features such as the modulation amplitude of the of the intensity of the STED beam and the signal to noise ratio of the acquired signals.

In the first preferred embodiment of the method of figure 2, the values of constants Mcenter and Mpe ph are determined as follows. The value of Mcenter corresponds to the minimum modulation va lue measured in the sample to the given experimental conditions. I n other words, the va lue of Mcenter corresponds to the minimum modulation value for which the number of pixels is greater than or equal to a percentage δ, optionally equal to 10%, of the peak value of the histogram of of M(x,y) values. The value of Mpe ph is calculated according to the equation :

wherein σι and σ₂ are the values of the standard deviation of the bi-dimensional gaussian distribution approximating a bi-dimensional histogram of g(x,y) and s(x,y) values.

I n other em bodiments, the values of constants Mcenter and Mpehph can be determined by using a known sa mple of fluorescent beads a nd evaluating the quality of the images obtained at different values of Mcenter and Mpehph.

Figure 7 illustrates some processing of data obtained in a first simulated experiment, corresponding to a first experiment ca rried out with the first preferred em bodiment described above of the method according to the invention. Said first experiment will be described in detail below with reference to figure 3. Figure 7a shows the histogram of the measured M(x,y) values from which Mcenter is obtained, Figure 7b shows the g(x,y) and s(x,y) values that are approximated with a bi-dimensional Gaussian distribution shown in figure 7c, whose standard deviation values are used for the calculation of Mpe ph.

By multiplying the fraction f_Center(x,y) for a STED image F'sTED(x,y) a final image of higher resolution is obtained : wherein the subscript "m-STED" indicates the STED technique with modulation of the intensity of the STED beam. The STED image F'sTED(x,y) used to generate the image FM-sTED(x,y) can be chosen between an image averaged in time <F-r-sTED(x,y,t)> and any image of the temporal series F_T-sTED(x,y,t).

All the above equations 2 to 9 are valid for slow modulations, i.e. for temporal variations slow with respect to r that is the fluorescence lifetime of the singlet excited state. Sincer is in the range of nanoseconds, it is in theory possible to modulate the STED intensity on any time scale T higher than nanoseconds, T»lns. The speed of the actual modulation is limited by the specific hardware device used to modulate the STED beam. In general, the response of the modulation hardware device represents the lower limit at the time scale that can be used by the method. In fact, a too slow modulation could be disadvantageous in applications wherein a data acquisition high speed is required. In the first preferred embodiment of the method according to the invention, acoustic-optical modulators (AOM, acoustic-optical modulator) commonly employed in confocal microscopes are used, which are much faster than any mechanical device. The response of the AOM is approximately limited to the transit time of the sound wave through the beam, typically < 50 ns. In further embodiments, opto-electronic modulators (EOM, Electro-optic modulator) are used, wherein the time limit of the response is equal to tens of nanoseconds.

In the first preferred embodiment, the modulation in step 200 has a temporal pattern corresponding to a linear ramp whereby from a minimum value \_mm- , optionally equal to zero at time to up to a maximum value l_max at time ti=to + T, i.e. in a time window T, wherein T is the modulation period. For which the intensity of the STED beam, set for the sake of simplicity also to=0, is given by:

In the first preferred embodiment of the method, the time window T coincides with the typical time of integration of intensity in a pixel of a confocal or STED microscope (pixel dwell time), whose values are equal to about 10 μ≤ up to about 100 μ≤.

The linear ramp temporal pattern has the advantage of making the analysis of the signal in the time domain more direct. In other embodiments of the present invention the modulation 200 of the STED beam is obtained with different temporal pattern, such as sinusoidal patterns. Returning to the flow chart of Figure 2, in a step 300 of the method, a series of images of fluorescence are acquired, i.e. are scanned, corresponding to the image function f_t-sTED(x,y,t) of Eq. 4, having increasing spatial resolution whereby the image at time T has the higher spatial resolution. The acquisition, i.e. the scanning, 300 of the fluorescence signal, is solved in time and synchronised with the modulation of the STED bea m. Hence, the modulation period equa l to T and the number of acquired time instants equal to n being set, the required temporal resolution is T/n. In a scanning STED microscope, moreover, the timed acquisition must be synchronized with the scanning parameters. I n the first preferred embodiment, T coincides with the dwell time per total pixel. I n other embodiments, the synchronization ca n be performed even with other parameters, for example with the line scan time. I n this latter case, the scanning of a horizontal line in a given vertical position is performed consecutively for n times, using a different intensity level at each scan, a nd then it goes to a line on a different vertical position. I n a ny case, at the end of acquisition 300, n fluorescence intensity values are associated with each pixel, which values will be processed in a step 400 to obtain the modulation M (x, y) of the image on each pixel given by Eq. 7. Optionally, at step 400 the intensity fluorescence values are processed to a lso calculate a n averaged in time image <FT- sTED(x,y,t)>. Once the modulation M of the image is obtained, in a step 500 the contribution cenier(x,y) to fluorescence, given by fluorophores positioned at the center of the observation volume, is calculated, whereby the fraction is obtained. I n fact, the spatial distribution a_IsTED (x', y') of the intensity of the STED beam on the sample is given by:

^(JI_STED ^ ^ ^{= I}STED(^t)^hSTED (^x y) _{Eq 12} wherein hsTEoiX' ) is a donut-shaped PSF of the STED microscope.

The fluorescence signal f(X', Y') from each single fluorophore inside the observation volume is approximated with an ex onential decay:

/(*',/,f) oc exp

with a time constant of depletion Xde_P(x', y') which depends accordingly on the donut-shape of the PSF, on the saturation va lue hat, on the maxim um intensity ISTED of the STED beam a nd the size of the time window T according to the following equation: Eq. 14 It should be noted that, in this case, (x',y') corresponds to the distance from the center of a given pixel (x,y), whereby one more fluorophores can contribute to the fluorescence signal measured in a single pixel F(x,y) .

The Equation 14 also shows that a spatial distribution of the depletion decay rate Tdep (^x'_> y') '^s provided determined by the spatial distribution of hsTED (χ', y').

Since hsTED (χ', y') has a donut shape, the signal from the fluorophores located at the periphery of the observation volume (maximum values of hsTEo) will decay more quickly while the signal from the center of the observation volume (values close to zero of hsTEo) will remain approximately constant, as schematically shown in the graph of Figure lb. In other words, the signal of the fluorophores located on the periphery of the observation volume will be modulated while the signal from the center of the observation volume will not be modulated. This modulation gradient between the center and the periphery of the volume is used to separate from the total signal, the signal coming from the center. Since the modulation of the linear ramp of the STED beam induces an exponential decay of the fluorophores during the time window T, the separation at step 500 is performed with the same formalism used in the SPLIT technique, with the difference that the period T is not determined by the fluorescence decay of the fluorophores in the order of nanoseconds, but rather it is arbitrarily selected by setting the period T of the ramp (T> 100ns). In order to obtain significant spatial variations of Tapp (χ', y'), from Eq. 14 It is clear that /_mox must be of the magnitude order of hat. By the fluorescence image given by fluorophores positioned at the center of the observation volume, in a step 600 the final image FM-sTED(x,y) given by Eq. 10 is obtained, with increased resolution with respect to the image that would be obtained without integrating the present method to the STED microscope, as can be seen schematically by comparing the actual observation volume obtained with standard STED microscopy of Figure la and that obtained with STED microscopy integrated with the method of the present invention of Figure lb.

Figure 3 shows the results of a first experiment simulated by the inventors who have applied a first preferred embodiment of the method of the present invention. The first simulated experiment consisted in having randomly located punctiform fluorophores in a sample and simulating an acquisition according to the first embodiment of the method with a modulation of the intensity of the STED beam in the time domain. The simulated image comprises 32x32 pixels, the pixel size is 20nm, the value of wo = 160nm, the value of Imax/isat = 4, n = 8. Close to zero, the donut shape hsTED(x',y') is simulated with a parabolic approximation hsTED(x',y') ~2r²/wo². A sample is illuminated in step 100 with constant excitation intensity and is illuminated in step 200 with a modulated STED intensity with a linear ramp changing the STED intensity from value 0 to value /_m0x in the time window T, shown in Figure 3a. Figure 3b shows a series of STED images having an increasing spatial resolution acquired at step 300. The image at t = 0 has the lower spatial resolution, while the image at t = T has the highest spatial resolution. The shown variations of the fluorescence intensity for each pixel were then used to calculate at step 400 a modulation M(x,y) of the image shown in pseudo-coloured scale in Figure 3c, which can be used to improve the STED image resolution that would occur with a continuous STED beam, with standard STED technique, equal to I = Lax and whose fluorescence image is shown in Figure 3d. By processing the modulation M of the image at step 500 through equation 8, which allows to calculate fcenter (x,y) described above, the contribution of the fluorophores located at the center of the observation volume, shown in the image of Figure 3a, is separated from the contribution of those situated at the periphery, whose fluorescence image is shown in figure 3f. The final image called m-STED, which is that of Figure 3and given by the contribution of the fluorescence of the fluorophores at the center, has a higher resolution than the STED image with the maximum resolution of Figure 3d, as shown by a comparison of the line profiles shown in Figure 3g (wherein it is also reported that obtained with confocal microscopy without STED technique).

In a second preferred embodiment of the method of figure 4 of the present invention, the time_ pattern of the modulation performed in a step 200' has a sinusoidal pattern with frequency ω :

where is a value of average intensity and Alsted is the intensity of modulation amplitude of the STED beam. The fluorescence signal f(x',Yy) coming from each single fluorophore inside the observation volume is also modulated according to the equation:

The fluorescence signal f (χ', y', t) is phase shifted with respect to a signal of the STED beam and will contain the main frequency ω and the higher harmonics due to the non-linearity of the exponential factor. The modulation M of the image can be obtained by analysing the fluorescence signal at a frequency ω or at higher frequencies. This can be done via a suitable hardware setting, wherein an analysis instrument acquires in a step 300' fluorescence signals from a sample and outputs, in a step 400', the values of the average intensity <FTSTED (X , y, t)> and modulation Μ_Ω (x,y) of the fluorescence signal directly. A suitable hardware setting can be for example a standard lock-in detecting apparatus at frequency ω, wherein a lock-in amplifier directly provides, in addition to the average value of intensity <FTSTED (X , y, t)>, the modulation values Μ_Ω (x,y) and phase φ_ω (x,y) during an acquisition 300' of fluorescence signals. Therefore, given an image F_t-srfD(x,y,t) obtained by modulating the intensity of the STED beam at a frequency ω, the increase of the resolution can be obtained in a step 600' by using the average intensity <FTSTED (X , y, t)> in pixel (x,y) and the signal modulation Μ_Ω (x,y) on the same pixel:

^m-SIED (^■*·^{; =} f center (^X? )(^i-S7"£D (^■*·^; 0) £q ]_Q'

where the fraction fcenter is calculated at a step 500' according to Equation 8:

f center ^ = ^{1 "} [^M _ω (^Χ )^{~ M} center ]/^ ^' periph ^{~ M} center \ _Eq g'

The values of the constants M_center and Mperiph can be determined as described above, taking into account that the values of g(x,y) and s(x,y) can be obtained from Μ_Ω (x,y) and φ_ω (x,y) according to the following formulas:

g(x, y) = Μ_ω (x, γ)∞^φ_ω (x, y))

s(x,y) =M_a(x,y)sm((l)_a(x,y)) _{Eq M}

Similarly to the case of the modulation of the intensity of the STED beam in the time domain, it is in theory possible to modulate the STED intensity on any frequency scale lower to Gigahertz, G>«2n/lns i.e. ω/2π«1 GHz. The frequency of the actual modulation is limited by the specific hardware device used to modulate the STED beam.

It should be noted that, in contrast to the first preferred embodiment of the method, in the case of the second preferred embodiment of the method step 400 for calculating the modulation M(x,y) of the image to obtain the final image F/w-srfD(x,y) at higher resolution is not necessary.

In a further embodiment of the method wherein the modulation is sinusoidal with frequency ω, M_Center = 0 is set, for which it is

J center ^x> y) — 1 ~ Eq. 8" where K corresponds to a constant having a positive value Mpe ph above defined between 0 and 1. Thus the final image that is obtained from Eq. 10' becomes:

Fm-STED (^X> y) ^Ft-STED (^X, y ))

K Eq. 10'

By changing the value of the constant K, it is possible to adjust the increase in resolution,

Figure 5 shows the results of a second experiment simulated by the inventors who have applied a second preferred embodiment of the method. The second simulated experiment consisted in using the same sample of the first experiment, i.e. with randomly located punctiform fluorophores, and simulating an acquisition according to the second embodiment of the method with a modulation of the intensity of the STED beam in the frequency domain. Similarly to the first simulated experiment the results of which are shown schematically in Figure 3, the simulated image is 32x32 pixels, the pixel size is 20nm, the value of wo = 160nm, the value of Imax/isat = 4, a modulation frequency f= /2 π is equal to three times the inverse of the pixel dwell time. In this specific case the authors have used (so that the STED intensity varies from 0 to a value l_max= IO+AISTED) . Close to zero, the donut shape hsTED(x',y') is simulated with a parabolic approximation hsTED(x',y')

A sample is illuminated at step 100 with constant excitation intensity and is illuminated at step 200' with a STED intensity modulated with a sine wave at a frequency co, shown in Figure 5a. At step 300' STED images are acquired, with detection techniques in lock-in at frequency co, whereby the fluorescence signal is modulated and phase shifted, with the highest modulation at the periphery of the observation volume as schematically shown in Figure 5b. The variations of the fluorescence intensity at frequency co for each pixel generate at step 400' a modulation M(x,y) of the image shown in pseudo-coloured scale in Figure 5c, which can be used to improve the STED image resolution that would occur with a continuous STED beam, with standard STED technique, equal to I = Lax and whose fluorescence image is shown in Figure 5d. By using Eq. 10', the contribution of the fluorophores located at the center of the observation volume, shown in the image of Figure 5e, is separated from the contribution of those situated at the periphery, whose fluorescence image is shown in figure 5f. The final image called m-STED, which is that of Figure 5e given by the contribution of the fluorescence of the fluorophores at the center, has a higher resolution than the STED image with the maximum resolution of Figure 5d, as shown by a comparison of the line profiles shown in Figure 5g.

The authors have also performed experimental tests on fluorescent beads analysed with a STED microscopy apparatus integrated with the first preferred embodiment of the method. Some experimental results are reported in Figure 6. The experiment was performed with a standard STED apparatus on a sample of fluorescent beads of about 40nm in size. The intensity of the STED beam was modulated (step 200) with an AOM driven in such a way as to generate a linear ramp, from a value 0 to a va lue l_max, shown in Figure 6a. The modulation of the intensity of the STED bea m was synchronized with the line scanning time in such a way as to acquire n = 16 consecutive lines at different power (see always Figure 6a). A series of fluorescence STED images at increasing spatial resolution were then acquired (step 300), of which the first for ISTED=0 and last for lsTED=lmax a re shown in Figure 6b. These series of images were subsequently processed (step 400) to obtain an image modulation M shown in Figure 6c. Figure 6d shows the STED image of fluorescence at maximum resolution. By applying, then, the SPLIT algorithm (step 500) to the image M of modulation the contribution by fluorophores positioned in the center of the observation volume, shown in Fig. 6e, is separated from that of the fluorophores positioned at the periphery of the observation volume, shown in Figure 6f. As for the results of the simulations reported in Figures 3 and 5, the final image called m-STED, which is that of Figure 6e given by the contribution of the fluorescence of the fluorophores at the center, has a higher resolution than the STED image with the maximum resolution of Figure 6d, as shown by a comparison of the line profiles shown in Figure 6g.

I n the foregoing the preferred embodiments were described a nd some modifications of this invention have been suggested, but it should be understood that those skilled in the art can make modifications and changes without departing from the relative scope of protection, as defined by the appended claims.

Claims

1. Method of stimulated emission depletion microscopy modulated at high spatial resolution (STED microscopy) comprising the following steps:

A. illuminating (100) a sample including one or more fluorophores in an observation volume with an excitation light beam having a constant excitation intensity, to energize said one or more fluorophores in the sample;

B. illuminating (200, 200') the sample with a donut-shaped light beam, to de-energize by stimulated emission said one or more excited fluorophores, the depletion light beam having an intensity ISTED modulated over time ISTED = IsTEo(t) with a temporal pattern having a modulation period T;

C. acquiring (300, 300') one or more fluorescence signals emitted from the sample by an image plane (x, y) of pixels;

D. obtaining (400, 400') an image modulation M(x,y) of the sample on the image plane (x, y) of pixel on the basis of said one or more fluorescence signals acquired at step C;

E. processing (500, 500') the image modulation M(x,y) to obtain a fraction /_cefiier(x,y) of fluorescence intensity on the image plane (x, y) of pixel, corresponding to a contribution to said one or more fluorescence signals, acquired at step C, originating from fluorophores present in the center of an observation volume of the sample;

F. calculating (600, 600') a final image FM-STED(X,Y) of the sample on the image plane (x, y) of pixel, from said one or more fluorescence signals acquired at step C and from the image modulation function M(x,y) obtained at step D, as follows:

Fm-STED(^X> y) ⁼ fcenter(^x> J⁷)^' ' STED (^x> J⁷)

wherein F'sTED(x,y) is selected between a fluorescence image <F-r-sTED(x,y,t)> averaged in time and any image of a temporal series F-r sTED(x,y,t).

2. Method according to claim 1, characterized in that step C includes acquiring a time series of fluorescence images F-r sTED(x,y,t), whereby step D includes calculating the image modulation M(x,y) from variations |6F-r-sTED(x,y,t)/f_t-sTED(x,y,t)| of intensity of the acquired fluorescence images F-r sTED(x,y,t).

3. The method according to claim 1 or 2, characterized in that step C is carried out by scanning an xy area of the sample in a synchronized manner with the modulation period T, whereby, a number of time points acquired being set equal to n, the scan temporal resolution required is T/n.

4. The method according to claim 3, characterized in that T is equal to a dwell time for the total pixel.

5. Method according to any one of claims 3 or 4, characterized in that T/n is equal to a line scan time and n consecutive lines at different STED power are acquired.

6. Method according to any one of the preceding claims 3 to 5, characterized in that the temporal pattern of ISTED has a linear pattern with a period T » Ins, optionally comprised between 10 μ≤ and 100μ≤, still more optionally equal to 50 μ≤.

7. Method according to claim 1, characterized in that the intensity ISTED is modulated at a frequency ω and step C includes a frequency acquisition of fluorescence signals with a lock-in acquisition apparatus and analysis directly following step D, providing as output a modulation M_ro(x,y) and also providing an output fluorescence image <FT-STED (x,y,t)> averaged in time, whereby the fi step F as:

where K is a positive constant ranging between 0 and 1.

8. Method according to claim 7, characterized in that the temporal pattern of ISTED has a sinusoidal pattern at frequency ω/2π «1 GHz, optionally between 10kHz to IMHz, still more optionally equal to 100kHz.