Karthik Vishwanath

Quantitative fluorescence lifetime spectroscopy in turbid media: comparison of theoretical, experimental and computational methods

A Monte Carlo model developed to simulate time-resolved fluorescence propagation in a semi-infinite turbid medium was validated against previously reported theoretical and computational results. Model simulations were compared to experimental measurements of fluorescence spectra and lifetimes on tissue-simulating phantoms for single and dual fibre-optic probe geometries. Experiments and simulations using a single probe revealed that scattering-induced artefacts appeared in fluorescence emission spectra, while fluorescence lifetimes were unchanged. Although fluorescence lifetime measurements are generally more robust to scattering artefacts than are measurements of fluorescence spectra, in the dual-probe geometry scattering-induced changes in apparent lifetime were predicted both from diffusion theory and via Monte Carlo simulation, as well as measured experimentally. In all cases, the recovered apparent lifetime increased with increasing scattering and increasing source-detector separation. Diffusion theory consistently underestimated the magnitude of these increases in apparent lifetime (predicting a maximum increase of approximately 15%), while Monte Carlo simulations and experiment were closely matched (showing increases as large as 30%). These results indicate that quantitative simulations of time-resolved fluorescence propagation in turbid media will be important for accurate recovery of fluorophore lifetimes in biological spectroscopy and imaging applications.

Advances in Quantitative UV-Visible Spectroscopy for Clinical and Pre-clinical Application in Cancer

Methods of optical spectroscopy that provide quantitative, physically or physiologically meaningful measures of tissue properties are an attractive tool for the study, diagnosis, prognosis, and treatment of various cancers. Recent development of methodologies to convert measured reflectance and fluorescence spectra from tissue to cancer-relevant parameters such as vascular volume, oxygenation, extracellular matrix extent, metabolic redox states, and cellular proliferation have significantly advanced the field of tissue optical spectroscopy. The number of publications reporting quantitative tissue spectroscopy results in the UV-visible wavelength range has increased sharply in the past three years, and includes new and emerging studies that correlate optically measured parameters with independent measures such as immunohistochemistry, which should aid in increased clinical acceptance of these technologies.

The effects of internal refractive index variation in near-infrared optical tomography: a finite element modelling approach

Near-infrared (NIR) tomography is a technique used to measure light propagation through tissue and generate images of internal optical property distributions from boundary measurements. Most popular applications have concentrated on female breast imaging, neonatal and adult head imaging, as well as muscle and small animal studies. In most instances a highly scattering medium with a homogeneous refractive index is assumed throughout the imaging domain. Using these assumptions, it is possible to simplify the model to the diffusion approximation. However, biological tissue contains regions of varying optical absorption and scatter, as well as varying refractive index. In this work, we introduce an internal boundary constraint in the finite element method approach to modelling light propagation through tissue that accounts for regions of different refractive indices. We have compared the results to data from a Monte Carlo simulation and show that for a simple two-layered slab model of varying refractive index, the phase of the measured reflectance data is significantly altered by the variation in internal refractive index, whereas the amplitude data are affected only slightly.

A Robust Monte Carlo Model for the Extraction of Biological Absorption and Scattering In Vivo

We have a toolbox to quantify tissue optical properties that is composed of specialized fiberoptic probes for UV-visible diffuse reflectance spectroscopy and a fast, scalable inverse Monte Carlo (MC) model. In this paper, we assess the robustness of the toolbox for quantifying physiologically relevant parameters from turbid tissue-like media. In particular, we consider the effects of using different instruments, fiberoptic probes, and instrument-specific settings for a wide range of optical properties. Additionally, we test the quantitative accuracy of the inverse MC model for extracting the biologically relevant parameters of hemoglobin saturation and total hemoglobin concentration. We also test the effect of double-absorber phantoms (hemoglobin and crocin to model the absorption of hemoglobin and beta carotene, respectively, in the breast) for a range of absorption and scattering properties. We include an assessment on which reference phantom serves as the best calibration standard to enable accurate extraction of the absorption and scattering properties of the target sample. We found the best reference-target phantom combinations to be ones with similar scattering levels. The results from these phantom studies provide a set of guidelines for extracting optical parameters from clinical studies.

Light transport in biological tissue using three-dimensional frequency-domain simplified spherical harmonics equations.

The accuracy of the commonly used diffusion approximation as used in diffuse optical tomography is known to be limited in cases involving strong absorption and in these situations a higher ordered approximation is necessary. In this study, a light transport model has been developed based upon the three-dimensional frequency-domain simplified spherical harmonics (SP(N)) approximation for orders up to N = 7. The SP(N) data are tested against a semi-infinite multi-layered Monte Carlo model. It has been shown that the SP(N) approximation for higher orders (N >1) provides an increase in accuracy over the diffusion equation specifically near sources and at boundaries of regions with increased optical absorption. It is demonstrated that the error of fluence calculated near the sources between the diffusion approximation and the SP(N) model (N = 7) can be as large as 60%, therefore limiting the use of the diffusion approximation for small animal imaging and in situations where optical changes near sources are critical for tomographic reconstructions.

Time-resolved photon migration in bi-layered tissue models

In this article, we describe a novel Monte Carlo code for time-integrated and time-resolved photon migration simulations of excitation and fluorescent light propagation (with reabsorption) in bi-layered models of biological tissues. The code was experimentally validated using bi-layered, tissue-simulating phantoms and the agreement between simulations and experiment was better than 3%. We demonstrate the utility of the code for quantitative clinical optical diagnostics in epithelial tissues by examining design characteristics for clinically compatible waveguides with arbitrarily complex source-detector configurations. Results for human colonic tissues included a quantitative comparison of simulation predictions with time-resolved fluorescence data measured in vivo and spatio-temporal visualizations of photon migration characteristics in tissue models in both two- and three-dimensions for source-detector configurations, including variable waveguide spacing, numerical aperture, and diameter. These results were then extended from surface point spectroscopy to imaging modalities for both time-gated (fluorescence lifetime) and steady-state (fluorescence intensity) experimental conditions. To illustrate the flexibility of this computational approach, time-domain results were extended to simulate predictions for frequency-domain instrumentation. This work is the first demonstration and validation of a time-domain, multi-wavelength photon transport model with these capabilities in layered turbid-media.

Quantitative optical spectroscopy can identify long-term local tumor control in irradiated murine head and neck xenografts

Noninvasive and longitudinal monitoring of tumor oxygenation status using quantitative diffuse reflectance spectroscopy is used to test whether a final treatment outcome could be estimated from early optical signatures in a murine model of head and neck cancer when treated with radiation. Implanted tumors in the flank of 23 nude mice are exposed to 39 Gy of radiation, while 11 animals exposed to sham irradiation serve as controls. Diffuse optical reflectance is measured from the tumors at baseline (prior to irradiation) and then serially until 17 days posttreatment. The fastest and greatest increase in baseline-corrected blood oxygen saturation levels are observed from the animals that show complete tumor regression with no recurrence 90 days postirradiation, relative to both untreated and treated animals with local recurrences. These increases in saturation are observed starting 5 days posttreatment and last up to 17 days posttreatment. This preclinical study demonstrates that diffuse reflectance spectroscopy could provide a practical method far more effective than the growth delay assay to prognosticate treatment outcome in solid tumors and may hold significant translational promise.

Portable, Fiber-Based, Diffuse Reflection Spectroscopy (DRS) Systems for Estimating Tissue Optical Properties

Steady-state diffuse reflection spectroscopy is a well-studied optical technique that can provide a noninvasive and quantitative method for characterizing the absorption and scattering properties of biological tissues. Here, we compare three fiber-based diffuse reflection spectroscopy systems that were assembled to create a light-weight, portable, and robust optical spectrometer that could be easily translated for repeated and reliable use in mobile settings. The three systems were built using a broadband light source and a compact, commercially available spectrograph. We tested two different light sources and two spectrographs (manufactured by two different vendors). The assembled systems were characterized by their signal-to-noise ratios, the source-intensity drifts, and detector linearity. We quantified the performance of these instruments in extracting optical properties from diffuse reflectance spectra in tissue-mimicking liquid phantoms with well-controlled optical absorption and scattering coefficients. We show that all assembled systems were able to extract the optical absorption and scattering properties with errors less than 10%, while providing greater than ten-fold decrease in footprint and cost (relative to a previously well-characterized and widely used commercial system). Finally, we demonstrate the use of these small systems to measure optical biomarkers in vivo in a small-animal model cancer therapy study. We show that optical measurements from the simple portable system provide estimates of tumor oxygen saturation similar to those detected using the commercial system in murine tumor models of head and neck cancer.

A three-dimensional finite element model and image reconstruction algorithm for time-domain fluorescence imaging in highly scattering media

In this work, development and evaluation of a three-dimensional (3D) finite element model (FEM) based on the diffusion approximation of time-domain (TD) near-infrared fluorescence light transport in biological tissue is presented. This model allows both excitation and fluorescence temporal point-spread function (TPSF) data to be generated for heterogeneous scattering and absorbing media of arbitrary geometry. The TD FEM is evaluated via comparisons with analytical and Monte Carlo (MC) calculations and is shown to provide a quantitative accuracy which has less than 0.72% error in intensity and less than 37 ps error for mean time. The use of the Born-Ratio normalized data is demonstrated to reduce data mismatch between MC and FEM to less than 0.22% for intensity and less than 22 ps in mean time. An image reconstruction framework, based on a 3D FEM formulation, is outlined and simulation results based on a heterogeneous mouse model with a source of fluorescence in the pancreas is presented. It is shown that using early photons (i.e. the photons detected within the first 200 ps of the TPSF) improves the spatial resolution compared to using continuous-wave signals. It is also demonstrated, as expected, that the utilization of two time gates (early and latest photons) can improve the accuracy both in terms of spatial resolution and recovered contrast.

Experimental validation of an inverse fluorescence Monte Carlo model to extract concentrations of metabolically relevant fluorophores from turbid phantoms and a murine tumor model

An inverse Monte Carlo based model has been developed to extract intrinsic fluorescence from turbid media. The goal of this work was to experimentally validate the model to extract intrinsic fluorescence of three biologically meaningful fluorophores related to metabolism from turbid media containing absorbers and scatterers. Experimental studies were first carried out on tissue-mimicking phantoms that contained individual fluorophores and their combinations, across multiple absorption, scattering, and fluorophore concentrations. The model was then tested in a murine tumor model to determine both the kinetics of fluorophore uptake as well as overall tissue fluorophore concentration through extraction of the intrinsic fluorescence of an exogenous contrast agent that reports on glucose uptake. Results show the model can be used to recover the true intrinsic fluorescence spectrum with high accuracy (R(2)=0.988) as well as accurately compute fluorophore concentration in both single and multiple fluorophores phantoms when appropriate calibration standards are available. In the murine tumor, the model-corrected intrinsic fluorescence could be used to differentiate drug dose injections between different groups. A strong linear correlation was observed between the extracted intrinsic fluorescence intensity and injected drug dose, compared with the distorted turbid tissue fluorescence.