Application of the stretched exponential function to fluorescence lifetime imaging

Photochemistry and Photophysics in Silica-Based Materials: Ultrafast and Single Molecule Spectroscopy Observation

Silica-based materials (SBMs) are widely used in catalysis, photonics, and drug delivery. Their pores and cavities act as hosts of diverse guests ranging from classical dyes to drugs and quantum dots, allowing changes in the photochemical behavior of the confined guests. The heterogeneity of the guest populations as well as the confinement provided by these hosts affect the behavior of the formed hybrid materials. As a consequence, the observed reaction dynamics becomes significantly different and complex. Studying their photobehavior requires advanced laser-based spectroscopy and microscopy techniques as well as computational methods. Thanks to the development of ultrafast (spectroscopy and imaging) tools, we are witnessing an increasing interest of the scientific community to explore the intimate photobehavior of these composites. Here, we review the recent theoretical and ultrafast experimental studies of their photodynamics and discuss the results in comparison to those in homogeneous media. The discussion of the confined dynamics includes solvation and intra- and intermolecular proton-, electron-, and energy transfer events of the guest within the SBMs. Several examples of applications in photocatalysis, (photo)sensors, photonics, photovoltaics, and drug delivery demonstrate the vast potential of the SBMs in modern science and technology.

Unveiling the Self-Assembling Behavior of 5-Fluorouracil and its N,N'-Dimethyl Derivative: A Spectroscopic and Microscopic Approach

Under physiological conditions, 5-fluorouracil (5-FU), an anticancer drug, self-assembles into fibrils by strong hydrogen-bonding network, whereas its N,N'-dimethyl derivative, 5-fluoro-1,3-dimethyluracil (5-FDMU), does not make fibrils due to lack of strong hydrogen-bonding motif. In vitro, 5-FU self-assembly is sensitive to physicochemical conditions like the pH and ionic strength of the solution, which tune the strength of the noncovalent driving forces. Here we report a surprising finding that the buffer, which is necessary to control the pH and is typically considered to be inert, also significantly influences 5-FU self-assembly, which indicates an important role of counterions in the fibril formation. We have also monitored concentration- and time-dependent fibrillar growth of 5-FU. Again, fibril growth process is probed under dynamic conditions using microfluidic platform. The self-assembly of 5-FU compared with its N,N'-dimethyl derivative shows lower cytotoxicity to the cultured human erythroleukemic cells (K562 cells), which plausibly states the reason behind the greater effectiveness of 5-FU derivative drugs than 5-FU itself.

Fluorescence lifetime imaging ophthalmoscopy

Imaging techniques based on retinal autofluorescence have found broad applications in ophthalmology because they are extremely sensitive and noninvasive. Conventional fundus autofluorescence imaging measures fluorescence intensity of endogenous retinal fluorophores. It mainly derives its signal from lipofuscin at the level of the retinal pigment epithelium. Fundus autofluorescence, however, can not only be characterized by the spatial distribution of the fluorescence intensity or emission spectrum, but also by a characteristic fluorescence lifetime function. The fluorescence lifetime is the average amount of time a fluorophore remains in the excited state following excitation. Fluorescence lifetime imaging ophthalmoscopy (FLIO) is an emerging imaging modality for in vivo measurement of lifetimes of endogenous retinal fluorophores. Recent reports in this field have contributed to our understanding of the pathophysiology of various macular and retinal diseases. Within this review, the basic concept of fluorescence lifetime imaging is provided. It includes technical background information and correlation with in vitro measurements of individual retinal metabolites. In a second part, clinical applications of fluorescence lifetime imaging and fluorescence lifetime features of selected retinal diseases such as Stargardt disease, age-related macular degeneration, choroideremia, central serous chorioretinopathy, macular holes, diabetic retinopathy, and retinal artery occlusion are discussed. Potential areas of use for fluorescence lifetime imaging ophthalmoscopy will be outlined at the end of this review.

Development of Low-Cost Instrumentation for Single Point Autofluorescence Lifetime Measurements

Autofluorescence lifetime measurements, which can provide label-free readouts in biological tissues, contrasting e.g. different types and states of tissue matrix components and different cellular metabolites, may have significant clinical potential for diagnosis and to provide surgical guidance. However, the cost of the instrumentation typically used currently presents a barrier to wider implementation. We describe a low-cost single point time-resolved autofluorescence instrument, exploiting modulated laser diodes for excitation and FPGA-based circuitry for detection, together with a custom constant fraction discriminator. Its temporal accuracy is compared against a "gold-standard" instrument incorporating commercial TCSPC circuitry by resolving the fluorescence decays of reference fluorophores presenting single and double exponential decay profiles. To illustrate the potential to read out intrinsic contrast in tissue, we present preliminary measurements of autofluorescence lifetime measurements of biological tissues ex vivo. We believe that the lower cost of this instrument could enhance the potential of autofluorescence lifetime metrology for clinical deployment and commercial development.

Interaction of amino acid-functionalized silver nanoparticles and Candida albicans polymorphs: A deep-UV fluorescence imaging study

The interaction of the tryptophan functionalized Ag nanoparticles and live Candida albicans cells was studied by synchrotron excitation deep-ultraviolet (DUV) fluorescence imaging at the DISCO beamline of Synchrotron SOLEIL. DUV imaging showed that incubation of the fungus with functionalized nanoparticles results in significant increase in the fluorescence signal. The analysis of the images revealed that the interaction of the nanoparticles with (pseudo)hyphae polymorphs of the diploid fungus was less pronounced than in the case of yeast cells or budding spores. The changes in the intensity of the fluorescence signals of the cells after incubation were followed in [327-353nm] and [370-410nm] spectral ranges that correspond to the fluorescence of tryptophan in non-polar and polar environment, respectively. As a consequence of the environmental sensitivity of the silver-tryptophan fluorescent nanoprobe, we were able to determine the possible accumulation sites of the nanoparticles. The analysis of the intensity decay kinetics showed that the photobleaching effects were more pronounced in the case of the functionalized nanoparticle treated cells. The results of time-integrated emission in the mentioned spectral ranges suggested that the nanoparticles penetrate the cells, but that the majority of the nanoparticles attach to the cells' surfaces.

Microscopic understanding of the conformational features of a protein–DNA complex

Protein–DNA interactions play crucial roles in different stages of genetic activities, such as replication of genome, initiation of transcription,etc.

Quantitative Viscosity Mapping Using Fluorescence Lifetime Measurements

Lubricant viscosity is a key driver in both the tribological performance and energy efficiency of a lubricated contact. Elastohydrodynamic (EHD) lubrication produces very high pressures and shear rates, conditions hard to replicate using conventional rheometry. In situ rheological measurements within a typical contact are therefore important to investigate how a fluid behaves under such conditions. Molecular rotors provide such an opportunity to extract the local viscosity of a fluid under EHD lubrication. The validity of such an application is shown by comparing local viscosity measurements obtained using molecular rotors and fluorescence lifetime measurements, in a model EHD lubricant, with reference measurements using conventional rheometry techniques. The appropriateness of standard methods used in tribology for high-pressure rheometry (combining friction and film thickness measurements) has been verified when the flow of EHD lubricant is homogeneous and linear. A simple procedure for calibrating the fluorescence lifetime of molecular rotors at elevated pressure for viscosity measurements is proposed.

Towards unsupervised fluorescence lifetime imaging using low dimensional variable projection

Analyzing large fluorescence lifetime imaging (FLIM) data is becoming overwhelming; the latest FLIM systems easily produce massive amounts of data, making an efficient analysis more challenging than ever. In this paper we propose the combination of a custom-fit variable projection method, with a Laguerre expansion based deconvolution, to analyze bi-exponential data obtained from time-domain FLIM systems. Unlike nonlinear least squares methods, which require a suitable initial guess from an experienced researcher, the new method is free from manual interventions and hence can support automated analysis. Monte Carlo simulations are carried out on synthesized FLIM data to demonstrate the performance compared to other approaches. The performance is also illustrated on real-life FLIM data obtained from the study of autofluorescence of daisy pollen and the endocytosis of gold nanorods (GNRs) in living cells. In the latter, the fluorescence lifetimes of the GNRs are much shorter than the full width at half maximum of the instrument response function. Overall, our proposed method contains simple steps and shows great promise in realising automated FLIM analysis of large data sets.

Imaging studies of temperature dependent photodegradation and self-healing in disperse orange 11 dye-doped polymers

Using confocal transmission imaging microscopy, we measure the temperature dependence of photodegradation and self-healing in disperse orange 11 (DO11) dye-doped (poly)methyl-methacrylate (PMMA) and polystyrene (PS). In both dye-doped polymers, an increase in sample temperature results in a greater photodegradation rate and degree of degradation, while also resulting in a slower recovery rate and larger recovery fraction. These results confirm the temperature dependence predictions of the modified correlated chromophore domain model (mCCDM) [B. R. Anderson and M. G. Kuzyk, Phys. Rev. E 89, 032601 (2014)]. Additionally, using quantitative fitting of the imaging data for DO11/PMMA, we determine the domain density parameter to be ρ = 1.19 (±0.25) × 10(-2) and the domain free energy advantage to be λ = 0.282 ± 0.015 eV, which are within the uncertainty of the values previously determined using amplified spontaneous emission as the probe method [S. K. Ramini et al., Polym. Chem. 4, 4948 (2013)]. Finally, while we find photodegradation and self-healing of DO11/PS to be qualitatively consistent with the mCCDM, we find that it is quantitatively incompatible with the mCCDM as recovery in DO11/PS is found to behave as a stretched (or double) exponential as a function of time.

Optimizing Laguerre expansion based deconvolution methods for analysing bi-exponential fluorescence lifetime images

Fast deconvolution is an essential step to calibrate instrument responses in big fluorescence lifetime imaging microscopy (FLIM) image analysis. This paper examined a computationally effective least squares deconvolution method based on Laguerre expansion (LSD-LE), recently developed for clinical diagnosis applications, and proposed new criteria for selecting Laguerre basis functions (LBFs) without considering the mutual orthonormalities between LBFs. Compared with the previously reported LSD-LE, the improved LSD-LE allows to use a higher laser repetition rate, reducing the acquisition time per measurement. Moreover, we extended it, for the first time, to analyze bi-exponential fluorescence decays for more general FLIM-FRET applications. The proposed method was tested on both synthesized bi-exponential and realistic FLIM data for studying the endocytosis of gold nanorods in Hek293 cells. Compared with the previously reported constrained LSD-LE, it shows promising results.

Thermo-Mechanical Compatibility of Viscoelastic Mortars for Stone Repair

The magnitude of the thermal stresses that originate in an acrylic-based repair material used for the reprofiling of natural sandstone is analyzed. This kind of artificial stone was developed in the late 1970s for its peculiar property of reversibility in an organic solvent. However, it displays a high thermal expansion coefficient, which can be a matter of concern for the durability either of the repair or of the underlying original stone. To evaluate this risk we propose an analytical solution that considers the viscoelasticity of the repair layer. The temperature profile used in the numerical evaluation has been measured in a church where artificial stone has been used in a recent restoration campaign. The viscoelasticity of the artificial stone has been characterized by stress relaxation experiments. The numerical analysis shows that the relaxation time of the repair mortar, originating from a low T g , allows relief of most of the thermal stresses. It explains the good durability of this particular repair material, as observed by the practitioners, and provides a solid scientific basis for considering that the problem of thermal expansion mismatch is not an issue for this type of stone under any possible conditions of natural exposure.

Lifetime of fluorescent dye molecules in dense aqueous suspensions of polystyrene nanoparticles

We study the lifetime of two common fluorescent dye molecules from the Alexa Fluor NHS Ester family dissolved in water in an opaque aqueous dispersion of dielectric polystyrene nanoparticles. We investigate the role of the dispersion composition by varying the particle concentration and adding SDS (sodium dodecyl sulfate) surfactant molecules. The observed strong changes in lifetime of Alexa 430 can be attributed to the relative contribution of radiative and non-radiative decay channels while the lifetime of the Alexa 488 dye depends only weakly on the sample composition. For Alexa 430, a dye with a rather low quantum yield in aqueous solution, the addition of polystyrene nanoparticles leads to a significant enhancement in quantum yield and an associated increase of the fluorescent lifetime by up to 55 %. We speculate that the increased quantum yield can be attributed to the hydrophobic effect on the structure of water in the boundary layer around the polystyrene particles in suspension. Adding SDS acts as a quencher. Over a range of particle concentrations the particle induced increase of the lifetime can be completely compensated by adding SDS.

Deconvolution of fluorescence lifetime imaging microscopy by a library of exponentials

Fluorescence lifetime microscopy imaging (FLIM) is an optic technique that allows a quantitative characterization of the fluorescent components of a sample. However, for an accurate interpretation of FLIM, an initial processing step is required to deconvolve the instrument response of the system from the measured fluorescence decays. In this paper, we present a novel strategy for the deconvolution of FLIM data based on a library of exponentials. Our approach searches for the scaling coefficients of the library by non-negative least squares approximations plus Thikonov/l₂ or l₁ regularization terms. The parameters of the library are given by the lower and upper bounds in the characteristic lifetimes of the exponential functions and the size of the library, where we observe that this last variable is not a limiting factor in the resulting fitting accuracy. We compare our proposal to nonlinear least squares and global non-linear least squares estimations with a multi-exponential model, and also to constrained Laguerre-base expansions, where we visualize an advantage of our proposal based on Thikonov/l₂ regularization in terms of estimation accuracy, computational time, and tuning strategy. Our validation strategy considers synthetic datasets subject to both shot and Gaussian noise and samples with different lifetime maps, and experimental FLIM data of ex-vivo atherosclerotic plaques and human breast cancer cells.

Thermal Recovery of Colloidal Quantum Dot Ensembles Following Photoinduced Dimming

Colloidal CdSe quantum dot (QD) core ensembles were photodimmed and allowed to recover in the dark using ambient thermal energy at a range of temperatures. Nonlinear thermal recovery is well described by a stretched exponential function, and further analysis yields an underlying probability distribution of rate constants. Casting the rate constants as a collection of first-order activated processes provides an activation barrier probability distribution with significant density at room-temperature thermal energy that peaks at 200 meV before decaying to zero. This treatment for the recovery transition intuitively describes the distributed kinetics observed and complements commonly proposed blinking mechanisms.

eSIP: A Novel Solution-Based Sectioned Image Property Approach for Microscope Calibration

Fluorescence confocal microscopy represents one of the central tools in modern sciences. Correspondingly, a growing amount of research relies on the development of novel microscopic methods. During the last decade numerous microscopic approaches were developed for the investigation of various scientific questions. Thereby, the former qualitative imaging methods became replaced by advanced quantitative methods to gain more and more information from a given sample. However, modern microscope systems being as complex as they are, require very precise and appropriate calibration routines, in particular when quantitative measurements should be compared over longer time scales or between different setups. Multispectral beads with sub-resolution size are often used to describe the point spread function and thus the optical properties of the microscope. More recently, a fluorescent layer was utilized to describe the axial profile for each pixel, which allows a spatially resolved characterization. However, fabrication of a thin fluorescent layer with matching refractive index is technically not solved yet. Therefore, we propose a novel type of calibration concept for sectioned image property (SIP) measurements which is based on fluorescent solution and makes the calibration concept available for a broader number of users. Compared to the previous approach, additional information can be obtained by application of this extended SIP chart approach, including penetration depth, detected number of photons, and illumination profile shape. Furthermore, due to the fit of the complete profile, our method is less susceptible to noise. Generally, the extended SIP approach represents a simple and highly reproducible method, allowing setup independent calibration and alignment procedures, which is mandatory for advanced quantitative microscopy.

FLIMX: A Software Package to Determine and Analyze the Fluorescence Lifetime in Time-Resolved Fluorescence Data from the Human Eye

Fluorescence lifetime imaging ophthalmoscopy (FLIO) is a new technique for measuring the in vivo autofluorescence intensity decays generated by endogenous fluorophores in the ocular fundus. Here, we present a software package called FLIM eXplorer (FLIMX) for analyzing FLIO data. Specifically, we introduce a new adaptive binning approach as an optimal tradeoff between the spatial resolution and the number of photons required per pixel. We also expand existing decay models (multi-exponential, stretched exponential, spectral global analysis, incomplete decay) to account for the layered structure of the eye and present a method to correct for the influence of the crystalline lens fluorescence on the retina fluorescence. Subsequently, the Holm-Bonferroni method is applied to FLIO measurements to allow for group comparisons between patients and controls on the basis of fluorescence lifetime parameters. The performance of the new approaches was evaluated in five experiments. Specifically, we evaluated static and adaptive binning in a diabetes mellitus patient, we compared the different decay models in a healthy volunteer and performed a group comparison between diabetes patients and controls. An overview of the visualization capabilities and a comparison of static and adaptive binning is shown for a patient with macular hole. FLIMX’s applicability to fluorescence lifetime imaging microscopy is shown in the ganglion cell layer of a porcine retina sample, obtained by a laser scanning microscope using two-photon excitation.

Blind deconvolution estimation of fluorescence measurements through quadratic programming

Time-deconvolution of the instrument response from fluorescence lifetime imaging microscopy (FLIM) data is usually necessary for accurate fluorescence lifetime estimation. In many applications, however, the instrument response is not available. In such cases, a blind deconvolution approach is required. An iterative methodology is proposed to address the blind deconvolution problem departing from a dataset of FLIM measurements. A linear combination of a base conformed by Laguerre functions models the fluorescence impulse response of the sample at each spatial point in our formulation. Our blind deconvolution estimation (BDE) algorithm is formulated as a quadratic approximation problem, where the decision variables are the samples of the instrument response and the scaling coefficients of the basis functions. In the approximation cost function, there is a bilinear dependence on the decision variables. Hence, due to the nonlinear nature of the estimation process, an alternating least-squares scheme iteratively solves the approximation problem. Our proposal searches for the samples of the instrument response with a global perspective, and the scaling coefficients of the basis functions locally at each spatial point. First, the iterative methodology relies on a least-squares solution for the instrument response, and quadratic programming for the scaling coefficients applied just to a subset of the measured fluorescence decays to initially estimate the instrument response to speed up the convergence. After convergence, the final stage computes the fluorescence impulse response at all spatial points. A comprehensive validation stage considers synthetic and experimental FLIM datasets of ex vivo atherosclerotic plaques and human breast cancer cell samples that highlight the advantages of the proposed BDE algorithm under different noise and initial conditions in the iterative scheme and parameters of the proposal.

Systematically constructing kinetic transition network in polypeptide from top to down: trajectory mapping

Molecular dynamics (MD) simulation is an important tool for understanding bio-molecules in microscopic temporal/spatial scales. Besides the demand in improving simulation techniques to approach experimental scales, it becomes more and more crucial to develop robust methodology for precisely and objectively interpreting massive MD simulation data. In our previous work [J Phys Chem B 114, 10266 (2010)], the trajectory mapping (TM) method was presented to analyze simulation trajectories then to construct a kinetic transition network of metastable states. In this work, we further present a top-down implementation of TM to systematically detect complicate features of conformational space. We first look at longer MD trajectory pieces to get a coarse picture of transition network at larger time scale, and then we gradually cut the trajectory pieces in shorter for more details. A robust clustering algorithm is designed to more effectively identify the metastable states and transition events. We applied this TM method to detect the hierarchical structure in the conformational space of alanine-dodeca-peptide from microsecond to nanosecond time scales. The results show a downhill folding process of the peptide through multiple pathways. Even in this simple system, we found that single common-used order parameter is not sufficient either in distinguishing the metastable states or predicting the transition kinetics among these states.

Applications of second-harmonic generation imaging microscopy in ovarian and breast cancer

In this perspective, we discuss how the nonlinear optical technique of second-harmonic generation (SHG) microscopy has been used to greatly enhance our understanding of the tumor microenvironment (TME) of breast and ovarian cancer. Striking changes in collagen architecture are associated with these epithelial cancers, and SHG can image these changes with great sensitivity and specificity with submicrometer resolution. This information has not historically been exploited by pathologists but has the potential to enhance diagnostic and prognostic capabilities. We summarize the utility of image processing tools that analyze fiber morphology in SHG images of breast and ovarian cancer in human tissues and animal models. We also describe methods that exploit the SHG physical underpinnings that are effective in delineating normal and malignant tissues. First we describe the use of polarization-resolved SHG that yields metrics related to macromolecular and supramolecular structures. The coherence and corresponding phase-matching process of SHG results in emission directionality (forward to backward), which is related to sub-resolution fibrillar assembly. These analyses are more general and more broadly applicable than purely morphology-based analyses; however, they are more computationally intensive. Intravital imaging techniques are also emerging that incorporate all of these quantitative analyses. Now, all these techniques can be coupled with rapidly advancing miniaturization of imaging systems to afford their use in clinical situations including enhancing pathology analysis and also in assisting in real-time surgical determination of tumor margins.

256 × 2 SPAD line sensor for time resolved fluorescence spectroscopy

We present a CMOS chip 256 × 2 single photon avalanche diode (SPAD) line sensor, 23.78 µm pitch, 43.7% fill factor, custom designed for time resolved emission spectroscopy (TRES). Integrating time-to-digital converters (TDCs) implement on-chip mono-exponential fluorescence lifetime pre-calculation allowing timing of 65k photons/pixel at 200 Hz line rate at 40 ps resolution using centre-of-mass method (CMM). Per pixel time-correlated single-photon counting (TCSPC) histograms can also be generated with 320 ps bin resolution. We characterize performance in terms of dark count rate, instrument response function and lifetime uniformity for a set of fluorophores with lifetimes ranging from 4 ns to 6 ns. Lastly, we present fluorescence lifetime spectra of multicolor microspheres and skin autofluorescence acquired using a custom built spectrometer. In TCSPC mode, time-resolved spectra are acquired within 5 minutes whilst in CMM mode spectral lifetime signatures are acquired within 2 ms for fluorophore in cuvette and 200 ms for skin autofluorescence. We demonstrate CMOS line sensors to be a versatile tool for time-resolved fluorescence spectroscopy by providing parallelized and flexible spectral detection of fluorescence decay.

A solution to debates over the behavior of mortality at old ages

As humans live longer, the precise modeling of mortality curves in very old age is becoming more important in aging research and public health. Here, we address a methodology that utilizes a modified stretched exponential survival function where a stretched exponent is relevant to heterogeneity in human populations. This function allows better estimation of the maximum human lifespan by providing a good description of the mortality curves in very old age. Demographic analysis of Swedish females over three recent decades revealed an important trend: the maximum human lifespan (existing around 125 years) gradually decreased at a constant rate of ~1.6 years per decade, while the characteristic life gradually increased at a constant rate of ~1.2 years per decade. This trend indicates that the number of aging people is increasingly concentrated at very old age, which is consistent with the definition of population aging. Importantly analyzing the stretched exponents would help in evaluating the heterogeneity trends in human populations.

In situ viscosity measurement of confined liquids

Viscosity heterogeneity in a thin glycerol film (170 nm) at high pressure was observed with fluorescence lifetime measurements.

Regulation of Sec16 levels and dynamics links proliferation and secretion

We currently lack a broader mechanistic understanding of the integration of the early secretory pathway with other homeostatic processes such as cell growth. Here, we explore the possibility that Sec16A, a major constituent of endoplasmic reticulum exit sites (ERES), acts as an integrator of growth factor signaling. Surprisingly, we find that Sec16A is a short-lived protein that is regulated by growth factors in a manner dependent on Egr family transcription factors. We hypothesize that Sec16A acts as a central node in a coherent feed-forward loop that detects persistent growth factor stimuli to increase ERES number. Consistent with this notion, Sec16A is also regulated by short-term growth factor treatment that leads to increased turnover of Sec16A at ERES. Finally, we demonstrate that Sec16A depletion reduces proliferation, whereas its overexpression increases proliferation. Together with our finding that growth factors regulate Sec16A levels and its dynamics on ERES, we propose that this protein acts as an integrator linking growth factor signaling and secretion. This provides a mechanistic basis for the previously proposed link between secretion and proliferation.

Blind end-member and abundance extraction for multispectral fluorescence lifetime imaging microscopy data

This paper proposes a new blind end-member and abundance extraction (BEAE) method for multispectral fluorescence lifetime imaging microscopy (m-FLIM) data. The chemometrical analysis relies on an iterative estimation of the fluorescence decay end-members and their abundances. The proposed method is based on a linear mixture model with positivity and sum-to-one restrictions on the abundances and end-members to compensate for signature variability. The synthesis procedure depends on a quadratic optimization problem, which is solved by an alternating least-squares structure over convex sets. The BEAE strategy only assumes that the number of components in the analyzed sample is known a spriori. The proposed method is first validated by using synthetic m-FLIM datasets at 15, 20, and 25 dB signal-to-noise ratios. The samples simulate the mixed response of tissue containing multiple fluorescent intensity decays. Furthermore, the results were also validated with six m-FLIM datasets from fresh postmortem human coronary atherosclerotic plaques. A quantitative evaluation of the BEAE was made against two popular techniques: minimum volume constrained nonnegative matrix factorization (MVC-NMF) and multivariate curve resolution-alternating least-squares (MCR-ALS). Our proposed method (BEAE) was able to provide more accurate estimations of the end-members: 0.32% minimum relative error and 13.82% worst-case scenario, despite different initial conditions in the iterative optimization procedure and noise effect. Meanwhile, MVC-NMF and MCR-ALS presented more variability in estimating the end-members: 0.35% and 0.34% for minimum errors and 15.31% and 13.25% in the worst-case scenarios, respectively. This tendency was also maintained for the abundances, where BEAE obtained 0.05 as the minimum absolute error and 0.12 in the worst-case scenario; MCR-ALS and MVC-NMF achieved 0.04 and 0.06 for the minimum absolute errors, and 0.15 and 0.17 under the worst-case conditions, respectively. In addition, the average computation time was evaluated for the synthetic datasets, where MVC-NMF achieved the fastest time, followed by BEAE and finally MCR-ALS. Consequently, BEAE improved MVC-NMF in convergence to a local optimal solution and robustness against signal variability, and it is roughly 3.6 time faster than MCR-ALS.

Estimation of the number of fluorescent end-members for quantitative analysis of multispectral FLIM data

Multispectral fluorescence lifetime imaging (m-FLIM) can potentially allow identifying the endogenous fluorophores present in biological tissue. Quantitative description of such data requires estimating the number of components in the sample, their characteristic fluorescent decays, and their relative contributions or abundances. Unfortunately, this inverse problem usually requires prior knowledge about the data, which is seldom available in biomedical applications. This work presents a new methodology to estimate the number of potential endogenous fluorophores present in biological tissue samples from time-domain m-FLIM data. Furthermore, a completely blind linear unmixing algorithm is proposed. The method was validated using both synthetic and experimental m-FLIM data. The experimental m-FLIM data include in-vivo measurements from healthy and cancerous hamster cheek-pouch epithelial tissue, and ex-vivo measurements from human coronary atherosclerotic plaques. The analysis of m-FLIM data from in-vivo hamster oral mucosa identified healthy from precancerous lesions, based on the relative concentration of their characteristic fluorophores. The algorithm also provided a better description of atherosclerotic plaques in term of their endogenous fluorophores. These results demonstrate the potential of this methodology to provide quantitative description of tissue biochemical composition.

Engineering of weak helper interactions for high-efficiency FRET probes

Fluorescence resonance energy transfer (FRET)-based detection of protein interactions is limited by the very narrow range of FRET-permitting distances. We show two different strategies for the rational design of weak helper interactions that co-recruit donor and acceptor fluorophores for a more robust detection of bimolecular FRET: (i) in silico design of electrostatically driven encounter complexes and (ii) fusion of tunable domain-peptide interaction modules based on WW or SH3 domains. We tested each strategy for optimization of FRET between (m)Citrine and mCherry, which do not natively interact. Both approaches yielded comparable and large increases in FRET efficiencies with little or no background. Helper-interaction modules can be fused to any pair of fluorescent proteins and could, we found, enhance FRET between mTFP1 and mCherry as well as between mTurquoise2 and mCitrine. We applied enhanced helper-interaction FRET (hiFRET) probes to study the binding between full-length H-Ras and Raf1 as well as the drug-induced interaction between Raf1 and B-Raf.

Dynamics of geminate rebinding of CO to cytochrome c in guanidine HCl probed by femtosecond vibrational spectroscopy

Femtosecond vibrational spectroscopy was used to probe the rebinding dynamics of CO to cytochrome c (Cytc) in 1.8 and 7 M guanidine HCl (GdnHCl) after photodeligation of the corresponding CO-bound protein in D2O buffer (pD = 7.4) at 283 K. Geminate rebinding (GR) dynamics of CO to the folded Cytc in 1.8 M GdnHCl (nCytc) is similar to that to chemically modified cytochrome c (cCytc), suggesting that the overall conformations of nCytcCO and cCytcCO are similar. About 86% of the dissociated CO molecules were geminately rebound to nCytc nonexponentially within 1 ns. The efficient GR of CO to the folded Cytc can be attributed to the organized protein matrix near the active site of nCytc that provides an efficient trap for the diffusing CO ligand after photodissociation. Although the concentration of nCytc did not affect its GR yield of CO, GR yield of CO to the unfolded Cytc in 7 M GdnHCl (uCytc) increased from 5 to 30% as the protein concentration increased from 0.3 to 9 mM. Time-resolved spectra of the (13)CO dissociated from both 9 mM nCytc(13)CO and 9 mM uCytc(13)CO showed a growing band with a peak at 2090 cm(-1) on the picosecond time scale, which was assigned to (13)CO in D2O solvent. At 1 ns, the fraction of the CO band in the solvent was about 10% of the nascent photodeligated protein in nCytc and more than 50% in the concentrated uCytc. Whereas a small opening in the active site of nCytc is responsible for the ultrafast escape of CO to solution in the folded protein, a large fraction of the CO escape to the solvent in uCytc results from the denatured structure of the active site in the unfolded protein. The spectrum of the CO dissociated from the concentrated uCytcCO contained a band that decayed as efficiently as that for the folded protein, suggesting that some fraction of uCytcCO might form aggregates even in 7 M denaturant, such that the aggregate acts as an efficient trap for the diffusing CO after deligation. No hint of precipitate in the concentrated uCytcCO and protein refolding upon dilution of the GdnHCl indicate that the aggregate does not grow continuously but remains as a soluble oligomer. The delayed appearance of the solvated CO and the inefficient GR of CO in uCytcCO suggest that the monomeric unfolded CytcCO so loosely arranged that the protein matrix cannot trap CO efficiently but the bound CO is still buried within hydrophobic residues even under the harsh denaturation condition.

Parallelized TCSPC for dynamic intravital fluorescence lifetime imaging: quantifying neuronal dysfunction in neuroinflammation

Two-photon laser-scanning microscopy has revolutionized our view on vital processes by revealing motility and interaction patterns of various cell subsets in hardly accessible organs (e.g. brain) in living animals. However, current technology is still insufficient to elucidate the mechanisms of organ dysfunction as a prerequisite for developing new therapeutic strategies, since it renders only sparse information about the molecular basis of cellular response within tissues in health and disease. In the context of imaging, Förster resonant energy transfer (FRET) is one of the most adequate tools to probe molecular mechanisms of cell function. As a calibration-free technique, fluorescence lifetime imaging (FLIM) is superior for quantifying FRET in vivo. Currently, its main limitation is the acquisition speed in the context of deep-tissue 3D and 4D imaging. Here we present a parallelized time-correlated single-photon counting point detector (p-TCSPC) (i) for dynamic single-beam scanning FLIM of large 3D areas on the range of hundreds of milliseconds relevant in the context of immune-induced pathologies as well as (ii) for ultrafast 2D FLIM in the range of tens of milliseconds, a scale relevant for cell physiology. We demonstrate its power in dynamic deep-tissue intravital imaging, as compared to multi-beam scanning time-gated FLIM suitable for fast data acquisition and compared to highly sensitive single-channel TCSPC adequate to detect low fluorescence signals. Using p-TCSPC, 256×256 pixel FLIM maps (300×300 µm(2)) are acquired within 468 ms while 131×131 pixel FLIM maps (75×75 µm(2)) can be acquired every 82 ms in 115 µm depth in the spinal cord of CerTN L15 mice. The CerTN L15 mice express a FRET-based Ca-biosensor in certain neuronal subsets. Our new technology allows us to perform time-lapse 3D intravital FLIM (4D FLIM) in the brain stem of CerTN L15 mice affected by experimental autoimmune encephalomyelitis and, thereby, to truly quantify neuronal dysfunction in neuroinflammation.

A new method to estimate abundances of multiple components using multi-spectral fluorescence lifetime imaging microscopy

Multi-Spectral Fluorescent Lifetime Imaging Microscopy (m-FLIM) is a technique that aims to perform noninvasive in situ clinical diagnosis of several diseases. It measures the endogenous fluorescence of molecules, recording their lifetime decay in different wavelength bands. This signal is a mixed response of multiple fluorescent components present in a tissue sample. The goal is to decompose the mixture and estimate the proportional contributions of its constituents. Estimation of such quantitative description will help to characterize the molecular constitution of a given sample. This paper presents a new method to estimate the abundances of multiple components present in a mixture measured using m-FLIM data. It provides a closed-form solution under the fully constrained linear unmixing model and assuming the number of components as well as their ideal lifetime decays are known. Its performance is tested using synthetic samples with three components, where performance can be measured accurately and the percentage error is around 6%. The algorithm was also validated performing unmixing of ex vivo data samples from atherosclerotic human tissue containing collagen, elastin and low-density lipoproteins. These experiments were validated against ground-truth maps, which only give a quantitative description, and the estimated accuracy was around 88%.

A fully constrained optimization method for time-resolved multispectral fluorescence lifetime imaging microscopy data unmixing

This paper presents a new unmixing methodology of multispectral fluorescence lifetime imaging microscopy (m-FLIM) data, in which the spectrum is defined as the combination of time-domain fluorescence decays at multiple emission wavelengths. The method is based on a quadratic constrained optimization (CO) algorithm that provides a closed-form solution under equality and inequality restrictions. In this paper, it is assumed that the time-resolved fluorescence spectrum profiles of the constituent components are linearly independent and known a priori. For comparison purposes, the standard least squares (LS) solution and two constrained versions nonnegativity constrained least squares (NCLS) and fully constrained least squares (FCLS) (Heinz and Chang, 2001) are also tested. Their performance was evaluated by using synthetic simulations, as well as imaged samples from fluorescent dyes and ex vivo tissue. In all the synthetic evaluations, the CO obtained the best accuracy in the estimations of the proportional contributions. CO could achieve an improvement ranging between 41% and 59% in the relative error compared to LS, NCLS, and FCLS at different signal-to-noise ratios. A liquid mixture of fluorescent dyes was also prepared and imaged in order to provide a controlled scenario with real data, where CO and FCLS obtained the best performance. The CO and FCLS were also tested with 20 ex vivo samples of human coronary arteries, where the expected concentrations are qualitatively known. A certainty measure was employed to assess the confidence in the estimations made by each algorithm. The experiments confirmed a better performance of CO, since this method is optimal with respect to equality and inequality restrictions in the linear unmixing formulation. Thus, the evaluation showed that CO achieves an accurate characterization of the samples. Furthermore, CO is a computational efficient alternative to estimate the abundance of components in m-FLIM data, since a global optimal solution is always guaranteed in a closed form.

Multimodal in vivo imaging of oral cancer using fluorescence lifetime, photoacoustic and ultrasound techniques

This work reports a multimodal system for label-free tissue diagnosis combining fluorescence lifetime imaging (FLIm), ultrasound backscatter microscopy (UBM), and photoacoustic imaging (PAI). This system provides complementary biochemical, structural and functional features allowing for enhanced in vivo detection of oral carcinoma. Results from a hamster oral carcinoma model (normal, precancer and carcinoma) are presented demonstrating the ability of FLIm to delineate biochemical composition at the tissue surface, UBM and related radiofrequency parameters to identify disruptions in the tissue microarchitecture and PAI to map optical absorption associated with specific tissue morphology and physiology.

Fluorescence lifetime imaging microscopy for quantitative biological imaging

Fluorescence lifetime imaging microscopy (FLIM) is a method for measuring fluorophore lifetimes with microscopic spatial resolution, providing a useful tool for cell biologists to detect, visualize, and investigate structure and function of biological systems. In this chapter, we begin by introducing the basic theory of fluorescence lifetime, including the characteristics of fluorophore decay, followed by a discussion of factors affecting fluorescence lifetimes and the potential advantages of fluorescence lifetime as a source of image contrast. Experimental methods for creating lifetime maps, including both time- and frequency-domain experimental approaches, are then introduced. Then, FLIM data analysis methods are discussed, including rapid lifetime determination, multiexponential fitting, Laguerre polynomial fitting, and phasor plot analysis. After, data analysis methods are introduced that improve lifetime precision of FLIM maps based upon optimal virtual gating and total variation denoising. The chapter concludes by highlighting several recent FLIM applications for quantitative biological imaging, including Förster resonance energy transfer-FLIM, fluorescence correlation spectroscopy-FLIM, multispectral-FLIM, and multiphoton-FLIM.

Quantitative fluorescence loss in photobleaching for analysis of protein transport and aggregation

Background

Fluorescence loss in photobleaching (FLIP) is a widely used imaging technique, which provides information about protein dynamics in various cellular regions. In FLIP, a small cellular region is repeatedly illuminated by an intense laser pulse, while images are taken with reduced laser power with a time lag between the bleaches. Despite its popularity, tools are lacking for quantitative analysis of FLIP experiments. Typically, the user defines regions of interest (ROIs) for further analysis which is subjective and does not allow for comparing different cells and experimental settings.

Results

We present two complementary methods to detect and quantify protein transport and aggregation in living cells from FLIP image series. In the first approach, a stretched exponential (StrExp) function is fitted to fluorescence loss (FL) inside and outside the bleached region. We show by reaction–diffusion simulations, that the StrExp function can describe both, binding/barrier–limited and diffusion-limited FL kinetics. By pixel-wise regression of that function to FL kinetics of enhanced green fluorescent protein (eGFP), we determined in a user-unbiased manner from which cellular regions eGFP can be replenished in the bleached area. Spatial variation in the parameters calculated from the StrExp function allow for detecting diffusion barriers for eGFP in the nucleus and cytoplasm of living cells. Polyglutamine (polyQ) disease proteins like mutant huntingtin (mtHtt) can form large aggregates called inclusion bodies (IB’s). The second method combines single particle tracking with multi-compartment modelling of FL kinetics in moving IB’s to determine exchange rates of eGFP-tagged mtHtt protein (eGFP-mtHtt) between aggregates and the cytoplasm. This method is self-calibrating since it relates the FL inside and outside the bleached regions. It makes it therefore possible to compare release kinetics of eGFP-mtHtt between different cells and experiments.

Conclusions

We present two complementary methods for quantitative analysis of FLIP experiments in living cells. They provide spatial maps of exchange dynamics and absolute binding parameters of fluorescent molecules to moving intracellular entities, respectively. Our methods should be of great value for quantitative studies of intracellular transport.

Application of non-negative matrix factorization to multispectral FLIM data analysis

Existing methods of interpreting fluorescence lifetime imaging microscopy (FLIM) images are based on comparing the intensity and lifetime values at each pixel with those of known fluorophores. This method becomes unwieldy and subjective in many practical applications where there are several fluorescing species contributing to the bulk fluorescence signal, and even more so in the case of multispectral FLIM. Non-negative matrix factorization (NMF) is a multivariate data analysis technique aimed at extracting non-negative signatures of pure components and their non-negative abundances from an additive mixture of those components. In this paper, we present the application of NMF to multispectral time-domain FLIM data to obtain a new set of FLIM features (relative abundance of constituent fluorophores). These features are more intuitive and easier to interpret than the standard fluorescence intensity and lifetime values. The proposed approach, unlike several FLIM data analysis methods, is not limited by the number of constituent fluorescing species or their possibly complex decay dynamics. Moreover, the new set of FLIM features can be obtained by processing raw multispectral FLIM intensity data, thereby rendering time deconvolution unnecessary and resulting in lesser computational time and relaxed SNR requirements. The performance of the NMF method was validated on simulated and experimental multispectral time-domain FLIM data. The NMF features were also compared against the standard intensity and lifetime features, in terms of their ability to discriminate between different types of atherosclerotic plaques.

A fluorescence approach to investigate repartitioning of coalescing agents in acrylic polymer emulsions

Repartitioning of co-solvents between particles of latex emulsions was investigated by means of a fluorescence method based on the detection of the amount of co-solvent via the solvatochromic shift of the emission maximum of a fluorescent probe, copolymerized at a low concentration. Complete repartitioning of co-solvents between particles of latex materials with a low T_g (ca. 25 °C) occurred within minutes. For a hydrophilic latex with a T_g of 68 °C, equilibration was achieved within an hour. Repartitioning was faster for more hydrophobic co-solvents. For a hydrophobic latex of similar T_g, co-solvent repartitioning took place on the same time scale, but complete equilibration was not reached. Possibly, there is an additional slow component in the repartitioning, or the prolonged presence of co-solvent causes a structural change in the latex particles that affects the outcome of the experiment.

Two-step anomalous diffusion tensor imaging

We extend the formalism of anomalous diffusion imaging to include directional anisotropy of fitted parameters. The resulting technique is termed anomalous diffusion tensor imaging (aDTI), and allows the directional properties of the distributed diffusion coefficient (α) and the anomalous diffusion exponent, (γ) to be analysed using the same analytical techniques as regular diffusion tensor imaging (DTI). Together, these parameters quantify the rate of diffusion (α) and the complexity of the diffusion environment (γ). We generated tensor images for the anomalous exponent tensor (Γ) and distributed diffusivity tensor (A) from in vivo human brain data and present images of eigenvalues, eigenvectors, Trace/3 (Tr), fractional anisotropy (FA) and tensor shape measures. In white matter, A is found to have a median Tr = 0.56 × 10(- 3) mm(2) s(- 1), FA = 0.58 and Γ Tr = 0.69, FA = 0.13. We observed that white matter shows a similar anisotropic geometry for the distributed diffusion tensor as for the regular diffusion tensor, whereas the anomalous exponent tensor exhibits a different shape characteristic which may be informative of tissue microstructure.

A novel method for fast and robust estimation of fluorescence decay dynamics using constrained least-squares deconvolution with Laguerre expansion

We report a novel method for estimating fluorescence impulse response function (fIRF) from noise-corrupted time-domain fluorescence measurements of biological tissue. This method is based on the use of high-order Laguerre basis functions and a constrained least-squares approach that addresses the problem of overfitting due to increased model complexity. The new method was extensively evaluated on fluorescence data from simulation, fluorescent standard dyes, ex vivo tissue samples of atherosclerotic plaques and in vivo oral carcinoma. Current results demonstrate that this method allows for rapid and accurate deconvolution of multiple channel fluorescence decays without adaptively adjusting the Laguerre scale parameter. The appropriate choice of the scale parameter is essential for accurate estimation of the fIRF. The method described here is anticipated to play an important role in the development of computational techniques for real-time analysis of time-resolved fluorescence data from biological tissues and to support the advancement of fluorescence lifetime instrumentation for biomedical diagnostics by providing a means for on-line robust analysis of fluorescence decay.

Phasor imaging with a widefield photon-counting detector

Fluorescence lifetime can be used as a contrast mechanism to distinguish fluorophores for localization or tracking, for studying molecular interactions, binding, assembly, and aggregation, or for observing conformational changes via Förster resonance energy transfer (FRET) between donor and acceptor molecules. Fluorescence lifetime imaging microscopy (FLIM) is thus a powerful technique but its widespread use has been hampered by demanding hardware and software requirements. FLIM data is often analyzed in terms of multicomponent fluorescence lifetime decays, which requires large signals for a good signal-to-noise ratio. This confines the approach to very low frame rates and limits the number of frames which can be acquired before bleaching the sample. Recently, a computationally efficient and intuitive graphical representation, the phasor approach, has been proposed as an alternative method for FLIM data analysis at the ensemble and single-molecule level. In this article, we illustrate the advantages of combining phasor analysis with a widefield time-resolved single photon-counting detector (the H33D detector) for FLIM applications. In particular we show that phasor analysis allows real-time subsecond identification of species by their lifetimes and rapid representation of their spatial distribution, thanks to the parallel acquisition of FLIM information over a wide field of view by the H33D detector. We also discuss possible improvements of the H33D detector's performance made possible by the simplicity of phasor analysis and its relaxed timing accuracy requirements compared to standard time-correlated single-photon counting (TCSPC) methods.