Experimental determination of photon propagation in highly absorbing and scattering media

Optical parameters estimation in inhomogeneous turbid media using backscattered light: for transcutaneous scattering measurement of intravascular blood

In our earlier research, a technique was developed to estimate the effective attenuation coefficient of subcutaneous blood vessels from the skin surface using the spatial distribution of backscattered near-infrared (NIR) light. The scattering effect in surrounding tissues was suppressed through the application of a differential principle, provided that the in vivo structure is known. In this study, a new method is proposed enabling the separate estimation of both scattering and absorption coefficients using NIR light of different wavelengths. The differential technique is newly innovated to make it applicable to the subcutaneous structure without requiring explicit geometrical information. Suppression of the scattering effect from surrounding tissue can be incorporated into the process of estimating the scattering and absorption coefficients. The validity of the proposed technique can be demonstrated through Monte Carlo simulations using both homogeneous and inhomogeneous tissue-simulating models. The estimated results exhibit good coherence with theoretical values (r2 = 0.988-0.999). Moreover, the vulnerability and robustness of the proposed technique against different measurement errors are verified. Optimal conditions for practical measurement are specified under various light-detection conditions. Separate estimation of scattering and absorption coefficients improves the accuracy of turbidity measurements and spectroscopy in biomedical applications considerably, particularly for noninvasive measurements and analysis of blood, lipids, and other components in subcutaneous blood vessels.

Optogenetic Modulation of Arrhythmia Triggers: Proof-of-Concept from Computational Modeling

Introduction

Early afterdepolarizations (EADs) are secondary voltage depolarizations associated with reduced repolarization reserve (RRR) that can trigger lethal arrhythmias. Relating EADs to triggered activity is difficult to study, so the ability to suppress or provoke EADs would be experimentally useful. Here, we use computational simulations to assess the feasibility of subthreshold optogenetic stimulation modulating the propensity for EADs (cell-scale) and EAD-associated ectopic beats (organ-scale).

Methods

We modified a ventricular ionic model by reducing rapid delayed rectifier potassium (0.25-0.1 × baseline) and increasing L-type calcium (1.0-3.5 × baseline) currents to create RRR conditions with varying severity. We ran simulations in models of single cardiomyocytes and left ventricles from post-myocardial infarction patient MRI scans. Optogenetic stimulation was simulated using either ChR2 (depolarizing) or GtACR1 (repolarizing) opsins.

Results

In cell-scale simulations without illumination, EADs were seen for 164 of 416 RRR conditions. Subthreshold stimulation of GtACR1 reduced EAD incidence by up to 84.8% (25/416 RRR conditions; 0.1 μW/mm2); in contrast, subthreshold ChR2 excitation increased EAD incidence by up to 136.6% (388/416 RRR conditions; 50 μW/mm2). At the organ scale, we assumed simultaneous, uniform illumination of the epicardial and endocardial surfaces. GtACR1-mediated suppression (10-50 μW/mm2) and ChR2-mediated unmasking (50-100 μW/mm2) of EAD-associated ectopic beats were feasible in three distinct ventricular models.

Conclusions

Our findings suggest that optogenetics could be used to silence or provoke both EADs and EAD-associated ectopic beats. Validation in animal models could lead to exciting new experimental regimes and potentially to novel anti-arrhythmia treatments.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12195-023-00781-z.

Optogenetic Stimulation Using Anion Channelrhodopsin (GtACR1) Facilitates Termination of Reentrant Arrhythmias With Low Light Energy Requirements: A Computational Study

Optogenetic defibrillation of hearts expressing light-sensitive cation channels (e.g., ChR2) has been proposed as an alternative to conventional electrotherapy. Past modeling work has shown that ChR2 stimulation can depolarize enough myocardium to interrupt arrhythmia, but its efficacy is limited by light attenuation and high energy needs. These shortcomings may be mitigated by using new optogenetic proteins like Guillardia theta Anion Channelrhodopsin (GtACR1), which produces a repolarizing outward current upon illumination. Accordingly, we designed a study to assess the feasibility of GtACR1-based optogenetic arrhythmia termination in human hearts. We conducted electrophysiological simulations in MRI-based atrial or ventricular models (n = 3 each), with pathological remodeling from atrial fibrillation or ischemic cardiomyopathy, respectively. We simulated light sensitization via viral gene delivery of three different opsins (ChR2, red-shifted ChR2, GtACR1) and uniform endocardial illumination at the appropriate wavelengths (blue, red, or green light, respectively). To analyze consistency of arrhythmia termination, we varied pulse timing (three evenly spaced intervals spanning the reentrant cycle) and intensity (atrial: 0.001–1 mW/mm²; ventricular: 0.001–10 mW/mm²). In atrial models, GtACR1 stimulation with 0.005 mW/mm² green light consistently terminated reentry; this was 10–100x weaker than the threshold levels for ChR2-mediated defibrillation. In ventricular models, defibrillation was observed in 2/3 models for GtACR1 stimulation at 0.005 mW/mm² (100–200x weaker than ChR2 cases). In the third ventricular model, defibrillation failed in nearly all cases, suggesting that attenuation issues and patient-specific organ/scar geometry may thwart termination in some cases. Across all models, the mechanism of GtACR1-mediated defibrillation was voltage forcing of illuminated tissue toward the modeled channel reversal potential of −40 mV, which made propagation through affected regions impossible. Thus, our findings suggest GtACR1-based optogenetic defibrillation of the human heart may be feasible with ≈2–3 orders of magnitude less energy than ChR2.

Axially slanted laser illumination scheme for direct and accurate Raman spectroscopic determination of gemcitabine concentration in freeze-dried gemcitabine injection powder housed in a glass container

When Raman spectroscopy is employed for a direct in situ determination of ingredient concentration for a product stored in a glass container, minimization of the interfering glass background in the collected spectrum is demanding to secure a more accurate analysis. To meet this request, an axially slanted illumination (ASI) scheme slantingly irradiating laser on the headspace side of a glass container and positioning a detector beneath the container was demonstrated in this study. This ASI scheme was basically designed to increase the distance between the laser illumination spot and detector location to minimize the number of glass photons reaching the detector. The analytical utility of the scheme was evaluated for the determination of gemcitabine concentration (42.9-58.2 wt%) in the gemcitabine injection powder housed in a glass container. Using the ASI scheme, the spectral features of the gemcitabine powder became distinct with only a weak underlying glass background signal. For comparative purpose, when an axially perpendicular offset (APO) scheme perpendicularly irradiating laser on the side wall where the sample was filled was used, the magnitude of glass background was higher, and the most intense gemcitabine peak was largely buried in the glass peak. The accuracy for determination of gemcitabine concentration using the ASI scheme was superior with an error of 0.20 wt%, while 0.33 wt% with employing the APO scheme. Overall, this study demonstrates that the ASI scheme is a potentially versatile Raman spectroscopic tool for fast non-sampling analysis of other products stored in a glass container.

Impact of the emission wavelengths on in vivo multiphoton imaging of mouse brains

Tissue scattering and absorption impact the excitation and emission light in different ways for multiphoton imaging. The collected fluorescence includes both ballistic photons and scattered photons whereas multiphoton excited signal within the focal volume is mostly generated by ballistic photons. The impact of excitation wavelengths on multiphoton imaging has been extensively investigated before; however, experimental data is lacking to evaluate the impact of emission wavelengths on fluorescence attenuation in deep imaging. Here we perform three-photon imaging of mouse brain vasculature in vivo using green, red, and near-infrared emission fluorophores, and compare quantitatively the attenuation of the fluorescence signal in the mouse brain at the emission wavelengths of 520 nm, 615 nm and 711 nm. Our results show that the emission wavelengths do not significantly influence the fluorescence collection efficiency. For the green, red and near-infrared fluorophores investigated, the difference in fluorescence collection efficiency is less than a factor of 2 at imaging depths between 0.6 and 1 mm. The advantage of long wavelength dyes for multiphoton deep imaging is almost entirely due to the long excitation wavelengths.

Signal-to-noise ratio analysis and improvement for fluorescence tomography imaging

CCD-based fluorescence tomography is widely used for small animal whole-body imaging. In this report, systematic signal-to-noise ratio (SNR) analyses of a fluorescence tomography imaging (FTI) system were performed, resulting in an easy-to-follow strategy to optimize hardware configurations and operational conditions for acquiring high-quality imaging data and for improving the overall system performance. Phantom experiments were conducted to demonstrate the performance improvement by these optimizations. The improved performance was further verified by imaging a tumor-bearing mouse in vivo. This report provides general and practical guidelines for setting up a high-performance electron multiplying charge coupled device based FTI system to achieve an optimized SNR, which can be useful for future FTI technology development.

High sensitivity optical measurement of skin gloss

We demonstrate a low-cost optical method for measuring the gloss properties with improved sensitivity in the low gloss regime, relevant for skin gloss properties. The gloss estimation method is based on, on the one hand, the slope of the intensity gradient in the transition regime between specular and diffuse reflection and on the other on the sum over the intensities of pixels above threshold, derived from a camera image obtained using unpolarized white light illumination. We demonstrate the improved sensitivity of the two proposed methods using Monte Carlo simulations and experiments performed on ISO gloss calibration standards with an optical prototype. The performance and linearity of the method was compared with different professional gloss measurement devices based on the ratio of specular to diffuse intensity. We demonstrate the feasibility for in-vivo skin gloss measurements by quantifying the temporal evolution of skin gloss after application of standard paraffin cream bases on skin. The presented method opens new possibilities in the fields of cosmetology and dermatopharmacology for measuring the skin gloss and resorption kinetics and the pharmacodynamics of various external agents.

Enhanced adaptive focusing through semi-transparent media

Adaptive optics can focus light through opaque media by compensating the random phase delay acquired while crossing a scattering curtain. The technique is commonly exploited in many fields, including astrophysics, microscopy, biomedicine and biology. A turbid lens has the capability of producing foci with a resolution higher than conventional optics, however it has a fundamental limit: to obtain a sharp focus one has to introduce a strongly scattering medium in the optical path. Indeed a tight focusing needs strong scattering and, as a consequence, high resolution focusing is obtained only for weakly transmitting samples. Here we describe a novel method allowing to obtain highly concentrated optical spots even by introducing a minimum amount of scattering in the beam path with semi-transparent materials. By filtering the pseudo-ballistic components of the transmitted beam we are able to experimentally overcome the limits of the adaptive focus resolution, gathering light on a spot with a diameter which is one third of the original speckle correlation function.

Non-invasive assessment of dairy products using spatially resolved diffuse reflectance spectroscopy

The quality of a dairy product is largely determined by its microstructure which also affects its optical properties. Consequently, an assessment of the optical properties during production may be part of a feedback system for ensuring the quality of the production process. This paper presents a novel camera-based measurement technique that enables robust quantification of a wide range of reduced scattering coefficients and absorption coefficients. Measurements are based on hyperspectral images of diffuse reflectance in the wavelength range of 470 to 1020 nm. The optical properties of commercially available milk and yogurt products with three different levels of fat content are measured. These constitute a relevant range of products at a dairy plant. The measured reduced scattering properties of the samples are presented and show a clear discrimination between levels of fat contents as well as fermentation. The presented measurement technique and method of analysis is thus suitable for a rapid, non-contact, and non-invasive inspection that can deduce physically interpretable properties.

Steady-state total diffuse reflectance with an exponential decaying source

The increasing preclinical and clinical utilization of digital cameras for photographic measurements of tissue conditions motivates the study of reflectance measurements obtained with planar illumination. We examine herein a formula that models the total diffuse reflectance measured from a semi-infinite medium using an exponentially decaying source, assuming continuous plane wave epi-illumination. The model is validated with experimental reflectance measurements from tissue mimicking phantoms. The need for adjusting the blood absorption spectrum due to pigment packaging is discussed along with the potential applications of the proposed formulation.

Modeling and experimental verification for a broad beam light transport in optical tomography

This paper describes a general theoretical model for computing a broad beam excitation light transport in a 3D diffusion medium. The model is based on the diffusion approximation of the radiative transport equation. An analytical approach for the light propagation is presented by deriving a corresponding Green's function. A finite cylindrical domain with a rectangular cross section was considered as a 3D homogeneous phantom model. The results of the model are compared with corresponding experimental data. The measurements are done on solid and liquid phantoms replicating tissue-like optical properties.

In vivo quantitative reconstruction studies of bioluminescence tomography: effects of peak-wavelength shift and model deviation

Bioluminescence tomography is a novel optical molecular imaging technology. The corresponding system, theory, and algorithmic frames have been set up. In the present study, we concentrated on the analysis of quantitative reconstruction deviation from peak-wavelength shift of luminescent source and the deviation of heterogeneous mouse model. The findings suggest that the reconstruction results are significantly affected by the peak-wavelength shift and deviation of anatomical structure animal models. Furthermore, the model deviations exhibit much more influence than the wavelength shift on the reconstruction results.

Truncated total least squares method with a practical truncation parameter choice scheme for bioluminescence tomography inverse problem

In bioluminescence tomography (BLT), reconstruction of internal bioluminescent source distribution from the surface optical signals is an ill-posed inverse problem. In real BLT experiment, apart from the measurement noise, the system errors caused by geometry mismatch, numerical discretization, and optical modeling approximations are also inevitable, which may lead to large errors in the reconstruction results. Most regularization techniques such as Tikhonov method only consider measurement noise, whereas the influences of system errors have not been investigated. In this paper, the truncated total least squares method (TTLS) is introduced into BLT reconstruction, in which both system errors and measurement noise are taken into account. Based on the modified generalized cross validation (MGCV) criterion and residual error minimization, a practical parameter-choice scheme referred to as improved GCV (IGCV) is proposed for TTLS. Numerical simulations with different noise levels and physical experiments demonstrate the effectiveness and potential of TTLS combined with IGCV for solving the BLT inverse problem.

A fiberoptic reflectance probe with multiple source-collector separations to increase the dynamic range of derived tissue optical absorption and scattering coefficients

Measurement of tissue optical absorption and (transport) reduced scattering coefficients (mu(a) and mu(s)', respectively) is fundamental to many applications of light in medicine and biology. We report a handheld fiberoptic probe to determine these coefficients by measuring the diffuse reflectance at multiple source-collector distances, which allows for a larger dynamic range than a single source-collector separation. Diffusion theory and a priori knowledge of the spectral shape of mu(a) and mu(s)' are used in a forward model of the diffuse reflectance. The dynamic range and accuracy of this method were evaluated using Monte Carlo simulations, phantom experiments and tissues in vivo.

Real-time intraoperative fluorescence imaging system using light-absorption correction

We present a novel fluorescence imaging system developed for real-time interventional imaging applications. The system implements a correction scheme that improves the accuracy of epi-illumination fluorescence images for light intensity variation in tissues. The implementation is based on the use of three cameras operating in parallel, utilizing a common lens, which allows for the concurrent collection of color, fluorescence, and light attenuation images at the excitation wavelength from the same field of view. The correction is based on a ratio approach of fluorescence over light attenuation images. Color images and video is used for surgical guidance and for registration with the corrected fluorescence images. We showcase the performance metrics of this system on phantoms and animals, and discuss the advantages over conventional epi-illumination systems developed for real-time applications and the limits of validity of corrected epi-illumination fluorescence imaging.

Drug development in oncology assisted by noninvasive optical imaging

Early and accurate detection of tumors, like the development of targeted treatments, is a major field of research in oncology. The generation of specific vectors, capable of transporting a drug or a contrast agent to the primary tumor site as well as to the remote (micro-) metastasis would be an asset for early diagnosis and cancer therapy. Our goal was to develop new treatments based on the use of tumor-targeted delivery of large biomolecules (DNA, siRNA, peptides, or nanoparticles), able to induce apoptosis while dodging the specific mechanisms developed by tumor cells to resist this programmed cell death. Nonetheless, the insufficient effectiveness of the vectorization systems is still a crucial issue. In this context, we generated new targeting vectors for drug and biomolecules delivery and developed several optical imaging systems for the follow-up and evaluation of these vectorization systems in live mice. Based on our recent work, we present a brief overview of how noninvasive optical imaging in small animals can accelerate the development of targeted therapeutics in oncology.

Integral equations of the photon fluence rate and flux based on a generalized Delta-Eddington phase function

We present a generalized Delta-Eddington phase function to simplify the radiative transfer equation to integral equations with respect to both photon fluence rate and flux vector. The photon fluence rate and flux can be solved from the system of integral equations. By comparing to the Monte Carlo simulation results, the solutions of the system of integral equations accurately model the photon propagation in biological tissue over a wide range of optical parameters.

An Adaptive Multigrid method for modeling photon transport through biological tissues in bioluminescence tomography

In this paper, we develop an Adaptive Multigrid method (AMGM) to model photon transport through biological tissues in bioluminescence tomography. In our method, the smoothing operation on fine levels and residual correction on coarse levels in V-Cycle offer fast convergence rate for this forward problem. Using a heterogeneous phantom, the methodology is validated by Monte Carlo simulations, and the computation speed is much higher than conventional smoothing iteration methods on a single grid. In actual biomedical imaging applications, especially when there are many sources in small animal body, AMGM is potential to accurately simulate the forward problem at very low computation cost.

Fluorescence molecular imaging

There is a wealth of new fluorescent reporter technologies for tagging of many cellular and subcellular processes in vivo. This imposed contrast is now captured with an increasing number of available imaging methods that offer new ways to visualize and quantify fluorescent markers distributed in tissues. This is an evolving field of imaging sciences that has already achieved major advances but is also facing important challenges. It is nevertheless well poised to significantly impact the ways of biological research, drug discovery, and clinical practice in the years to come. Herein, the most pertinent technologies associated with in vivo noninvasive or minimally invasive fluorescence imaging of tissues are summarized. Focus is given to small-animal imaging. However, while a broad spectrum of fluorescence reporter technologies and imaging methods are outlined, as necessary for biomedical research, and clinical translation as well.

Image analysis methods for diffuse optical tomography

Three major analytical tools in imaging science are summarized and demonstrated relative to optical imaging in vivo. Standard resolution testing is optimal when infinite contrast is used and hardware evaluation is the goal. However, deep tissue imaging of absorption or fluorescent contrast agents in vivo often presents a different problem, which requires contrast-detail analysis. This analysis shows that the minimum detectable sizes are in the range of 1/10 the outer diameter, whereas minimum detectable contrast values are in the range of 10 to 20% relative to the continuous background values. This is estimated for objects being in the center of the domain being imaged, and as the heterogeneous region becomes closer to the surface, the lower limit on size and contrast can become arbitrarily low and more dictated by hardware specifications. Finally, if human observer detection of abnormalities in the images is the goal, as is standard in most radiological practice, receiver operating characteristic (ROC) curve and location receiver operating characteristic curve (LROC) are used. Each of these three major areas of image interpretation and analysis are reviewed in the context of medical imaging as well as how they are used to quantify the performance of diffuse optical imaging of tissue.

Photon diffusion coefficient in scattering and absorbing media

We present a unified derivation of the photon diffusion coefficient for both steady-state and time-dependent transport in disordered absorbing media. The derivation is based on a modal analysis of the time-dependent radiative transfer equation. This approach confirms that the dynamic diffusion coefficient is given by the random-walk result D = cl(*)/3, where l(*) is the transport mean free path and c is the energy velocity, independent of the level of absorption. It also shows that the diffusion coefficient for steady-state transport, often used in biomedical optics, depends on absorption, in agreement with recent theoretical and experimental works. These two results resolve a recurrent controversy in light propagation and imaging in scattering media.

Boundary integral method for bioluminescence tomography

Bioluminescence tomography (BLT) allows in vivo localization and quantification of bioluminescent sources inside a small animal to reveal various molecular and cellular activities. We develop a reconstruction method to identify such a bioluminescent source distribution using the boundary integral method. Based on the diffusion model of the photon propagation in the biological tissue, this method incorporates a priori knowledge to define the permissible source region, and establish a direct linear relationship between measured body surface data and an unknown bioluminescent source distribution to enhance numerical stability and efficiency. The feasibility of the proposed BLT algorithm is demonstrated in heterogeneous mouse chest phantom studies.

Volumetric tomography of fluorescent proteins through small animals in vivo

Volumetric detection and accurate quantification of fluorescent proteins in entire animals would greatly enhance our ability to monitor biological processes in vivo. Here we present a quantitative tomographic technique for visualization of superficial and deep-seated (>2-3 mm) fluorescent protein activity in vivo. We demonstrate noninvasive imaging of lung tumor progression in a murine model, as well as imaging of gene delivery using a herpes virus vector. This technology can significantly improve imaging capacity over the current state of the art and should find wide in vivo imaging applications in drug discovery, immunology, and cancer research.