Monte Carlo simulation of near infrared autofluorescence measurements of in vivo skin

The autofluorescence properties of normal human skin in the near-infrared (NIR) spectral range were studied using Monte Carlo simulation. The light-tissue interactions including scattering, absorption and anisotropy propagation of the regenerated autofluorescence photons in the skin tissue were taken into account in the theoretical modeling. Skin was represented as a turbid seven-layered medium. To facilitate the simulation, ex vivo NIR autofluorescence spectra and images from different skin layers were measured from frozen skin vertical sections to define the intrinsic fluorescence properties. Monte Carlo simulation was then used to study how the intrinsic fluorescence spectra were distorted by the tissue reabsorption and scattering during in vivo measurements. We found that the reconstructed model skin spectra were in good agreement with the measured in vivo skin spectra from the same anatomical site as the ex vivo tissue sections, demonstrating the usefulness of this modeling. We also found that difference exists over the melanin fluorescent wavelength range (880-910 nm) between the simulated spectrum and the measured in vivo skin spectrum from a different anatomical site. This difference suggests that melanin contents may affect in vivo skin autofluorescence properties, which deserves further investigation.     

In Vivo Optical Properties of Melanocytic Skin Lesions: Common Nevi, Dysplastic Nevi and Malignant Melanoma

We present an in vivo study of the optical properties of common nevi, dysplastic nevi and malignant melanoma skin lesions in human subjects. Reflectance spectra were measured on 1379 skin lesions, in the visible and near-infrared spectral regions, using a spectral imaging system, in a clinical setting. Analysis of the data using a reflectance model revealed differences between the optical properties of melanin present in nevi and melanoma lesions. These differences, which are in agreement with our previous observations on average reflectance spectra, may be potentially useful for the noninvasive characterization of pigmented skin lesions and the early diagnosis of melanoma.     

A Geant4 Monte Carlo method for synthesizing radioxenon beta-gamma coincidence spectra

This paper describes a Monte Carlo (MC) simulation method for calculating radioxenon beta-gamma coincidence spectral information. These spectral components include detector response simulations by Geant4 modeling, detector energy and resolution calibrations using the histograms of detector response, beta-gamma coincidence efficiency values and spectral interference ratios. The work presented in this paper demonstrates the feasibility of using the spectral information to create beta-gamma coincidence spectra at various radioxenon activity concentrations. The analysis of these synthetic spectra by XECON software shows an excellent correlation between the analysed radioxenon activity concentration and number of MC samplings.     

Status of the Monte Carlo library least-squares (MCLLS) approach for non-linear radiation analyzer problems

The Center for Engineering Applications of Radioisotopes (CEAR) has been working for over a decade on the Monte Carlo library least-squares (MCLLS) approach for treating non-linear radiation analyzer problems including: (1) prompt gamma-ray neutron activation analysis (PGNAA) for bulk analysis, (2) energy-dispersive X-ray fluorescence (EDXRF) analyzers, and (3) carbon/oxygen tool analysis in oil well logging. This approach essentially consists of using Monte Carlo simulation to generate the libraries of all the elements to be analyzed plus any other required background libraries. These libraries are then used in the linear library least-squares (LLS) approach with unknown sample spectra to analyze for all elements in the sample. Iterations of this are used until the LLS values agree with the composition used to generate the libraries. The current status of the methods (and topics) necessary to implement the MCLLS approach is reported. This includes: (1) the Monte Carlo codes such as CEARXRF, CEARCPG, and CEARCO for forward generation of the necessary elemental library spectra for the LLS calculation for X-ray fluorescence, neutron capture prompt gamma-ray analyzers, and carbon/oxygen tools; (2) the correction of spectral pulse pile-up (PPU) distortion by Monte Carlo simulation with the code CEARIPPU; (3) generation of detector response functions (DRF) for detectors with linear and non-linear responses for Monte Carlo simulation of pulse-height spectra; and (4) the use of the differential operator (DO) technique to make the necessary iterations for non-linear responses practical. In addition to commonly analyzed single spectra, coincidence spectra or even two-dimensional (2-D) coincidence spectra can also be used in the MCLLS approach and may provide more accurate results.     

SPECTROSCOPIC AND MICROSCOPIC CHARACTERISTICS OF HUMAN SKIN AUTOFLUORESCENCE EMISSION

To improve the understanding of human skin autofluorescence emission, the spectroscopic and microscopic characteristics of skin autofluorescence were studied using a combined fluorescence and reflectance spectroanalyzer and a fiber optic microspectrophotometer. The autofluorescence spectra of in vivo human skin were measured over a wide excitation wavelength range (350–470 nm). The excitation–emission matrices of in vivo skin were obtained. An excitation–emission maximum pair (380 nm, 470 nm) was identified. It was revealed that the most probable energy of skin autofluorescence emission photons increases monotonically and near linearly with increasing excitation photon energy. It was demonstrated that the diffuse reflectance, R, can be used as a first order approximation of the fluorescence distortion factor f to correct the measured in vivo autofluorescence spectra for the effect of tissue reabsorption and scattering. The microscopic in vitro autofluorescence properties of excised skin tissue sections were examined using 442 nm He–Cd laser light excitation as an example. It was demonstrated that the fluorophore distribution inside the skin tissue is not uniform and the shapes of the autofluorescence spectra of different anatomical skin layers vary. The result of this study confirms that the major skin fluorophores are located in the dermis and provides an excellent foundation for Monte Carlo modeling of in vivo autofluorescence measurements.     

Accuracy of Noninvasive in vivo Measurements of Photosensitizer Uptake Based on a Diffusion Model of Reflectance Spectroscopy

Abstract— This study compares the photosensitizer concentration measured noninvasively in vivo by diffuse reflectance spectroscopy with the results of postmortem tissue solubilization and fluorometric assay. The reflectance spectrometer consists of a fiber optic surface probe, spectrometer and charge-coupled device (CCD) array detector. The surface probe has eight detection fibers separated from the light source fiber by distances ranging from 0.85 to 10 mm. The imaging spectrometer disperses the light from each detector fiber onto the two-dimensional CCD array, while maintaining spatial separation of each individual spectrum. A single exposure of the CCD therefore captures the reflectance spectrum at eight distances and over a range of 300 nm. From the spectra, the tissue's optical scattering and absorption coefficients are determined using a diffusion model of light propagation. Changes in the tissue absorption are used to estimate the photosensitizer concentration. Normal New Zealand White rabbits were injected with aluminum phthalocyanine tetrasulfonate (AlPcS4) and probe measurements made 24 h after injection on the dorsal skin, on muscle after surgically turning the skin back and on liver. For skin, the noninvasive estimate is proportional to the true concentration but low by a factor of 3. Based on Monte Carlo modeling of multilayered systems, this underestimate is attributed to the layered structure of the skin and nonuniform AIPcS4 distribution. A comparison of the noninvasive concentration estimates to the postmortem assay results finds good agreement for liver tissue even though application of the diffusion model is not strictly justified.     

Optical properties of in vitro epidermis and their possible relationship with optical properties of in vivo skin.

Using a spectrophotometer with an internal integrating sphere, the absorption (mu a) and reduced scattering (mu s') coefficients of in vitro epidermis were evaluated from reflectance and transmittance measurements. mu a and mu s' varied from 24 to 0.2 cm-1, and from 32 to 21 cm-1 respectively, on passing from 400 to 800 nm. Moreover, using an external integrating sphere, the reflectance spectrum of in vivo skin was compared with the reflectance spectrum calculated with a Monte Carlo model, in which the mean values of mu a and mu s' and different anisotropy parameters were used as input data. In vivo results show that the principle of similarity is entirely valid for wavelengths greater than 600 nm and may be considered a good approximation in the 400-600 nm band, and suggest that optical characteristics of in vivo skin may be inferred from reflectance measurements.     

IMMEDIATE PIGMENT DARKENING: VISUAL AND REFLECTANCE SPECTROPHOTOMETRIC ANALYSIS OF ACTION SPECTRUM

Immediate pigment darkening (IPD) occurs in human skin upon exposure to ultraviolet-A and visible radiation. The spectral changes that occur during IPD were measured with a rapid scanning reflectance spectrophotometer (RS) which employs optical fiber bundles for delivery and detection of light between 400 and 750 nm. The radiation dose dependence and wavelength dependence (334-549 nm irradiation) of IPD were studied by both the classical visual grading method and by spectrophotometric scoring using the RS system. The spectral changes that occur at long wavelengths with IPD mimic the natural absorption spectrum of melanin. Therefore, the IPD was scored in terms of the apparent change in melanin optical density, using the method Kollias and Baqer [Photochem. Photobiol. 43, 49-54 (1986)], based on reflectance in the 620-720 nm range. The nonlinearity of the visual grading method is demonstrated. The degree of IPD is first-order with respect to delivered dose and saturates after high doses. The maximum amount of IPD attained at saturation is greater for shorter wavelengths. Extrapolation of the reflectance data suggests the longest wavelength capable of eliciting IPD is about 470 nm.     

An ensemble of Monte Carlo uninformative variable elimination for wavelength selection

An improved method based on an ensemble of Monte Carlo uninformative variable elimination (EMCUVE) is presented for wavelength selection in multivariate calibration of spectral data. The proposed algorithm introduces Monte Carlo (MC) strategy to uninformative variable elimination-PLS (UVE-PLS) instead of leave-one-out strategy for estimating the contributions of each wavelength variable in the PLS model. In EMCUVE wavelength variables are evaluated by different Monte Carlo uninformative variable elimination (MCUVE) models. Moreover, a fusion of MCUVE and the vote rule can obtain an improvement over the original uninformative variable elimination method. Results obtained from simulated data and real data sets demonstrate that EMCUVE can properly carry out wavelength selection in the course of data analysis and improve predictive ability for multivariate calibration model.     

Porphyrin bleaching and PDT-induced spectral changes are irradiance dependent in ALA-sensitized normal rat skin in vivo

Photobleaching kinetics of aminolevulinic acid-induced protoporphyrin IX (PpIX) were measured in the normal skin of rats in vivo using a technique in which fluorescence spectra were corrected for the effects of tissue optical properties in the emission spectral window through division by reflectance spectra acquired in the same geometry and wavelength interval and for changes in excitation wavelength optical properties using diffuse reflectance measured at the excitation wavelength. Loss of PpIX fluorescence was monitored during photodynamic therapy (PDT) performed using 514 nm irradiation. Bleaching in response to irradiances of 1, 5 and 100 mW cm-2 was evaluated. The results demonstrate an irradiance dependence to the rate of photobleaching vs irradiation fluence, with the lowest irradiance leading to the most efficient loss of fluorescence. The kinetics for the accumulation of the primary fluorescent photoproduct of PpIX also exhibit an irradiance dependence, with greater peak accumulation at higher irradiance. These findings are consistent with a predominantly oxygen-dependent photobleaching reaction mechanism in vivo, and they provide spectroscopic evidence that PDT delivered at low irradiance deposits greater photodynamic dose for a given irradiation fluence. We also observed an irradiance dependence to the appearance of a fluorescence emission peak near 620 nm, consistent with accumulation of uroporphyrin/coproporphyrin in response to mitochondrial damage.     

The Dynamics of Laser-Induced Changes in Human Skin Autofluorescence—Experimental Measurements and Theoretical Modeling

To study the temporal dynamics of human skin autofluorescence photobleaching, we measured the autofluorescence spectral changes of skin in vivo during continuous exposure to 442 nm (He-Cd) laser light. Integral intensities were calculated for various spectral wavelength bands and plotted as a function of time. Mathematical analysis of the time function revealed a double-exponential photobleaching process: I(t) = a exp(-t/τ¹,) + b exp(-t/τ²) + c, in which t¹, and t² differed by an order of magnitude. A hypothesis for the mechanism of the double-exponential photobleaching dynamics was proposed and evaluated using Monte Carlo modeling of light propagation in the skin and autofluorescence escape from skin. By combining the fluorophore microdistributions, Monte Carlo simulation results and the variation in fluorescence decrease parameters (a, b, c, τ¹τ²) with increasing exposure intensities a biophysical explanation for the double-exponential photobleaching function was elucidated. The fast decrease term corresponds to laser-induced photobleaching in the stratum corneum, while the slow decrease term represents fluorophore changes in the dermis. The measured autofluorescence photobleaching dynamics can be used to determine the fractional contributions of different skin layers to the total autofluorescence signal measured in vivo.     

The effect of UV absorbing sunscreens on the reflectance and the consequent protection of skin

The in vivo reflectance spectra of Caucasian skin, coated with preparations containing sunscreen vehicle, vehicle with olive oil and vehicle with the UVB and UVA absorbers 2-ethylhexyl-4-methoxycinnamate and 4-t-butyl-4'-methoxydibenzoylmethane were determined. All preparations reduced the reflectance of skin throughout the UVA spectral range (320 to 400 nm), with the sunscreen preparations containing the UVB and UVB plus UVA absorbers reducing the reflectance more than the sunscreen vehicle alone. This phenomenon, which facilitates the penetration of UV radiation to the lower epidermis and dermal layers of skin and therefore lessens sunscreen efficacy, is attributed to optical coupling mediated by refractive index matching of the sunscreen to the upper epidermis. The greater reduction in skin diffuse reflectance caused by sunscreens containing methoxycinnamate is associated with this compound's high refractive index. Also, by determining the excitation spectra of the autofluorescence originating from the dermal layer of skin, the transmission spectra of the various components of sunscreen on skin were established, and these were in good general agreement with previously published spectra.     

A variable differential consensus method for improving the quantitative near-infrared spectroscopic analysis

Consensus methods have presented promising tools for improving the reliability of quantitative models in near-infrared(NIR) spectroscopic analysis.A strategy for improving the performance of consensus methods in multivariate calibration of NIR spectra is proposed.In the approach,a subset of non-collinear variables is generated using successive projections algorithm(SPA) for each variable in the reduced spectra by uninformative variables elimination(UVE).Then sub-models are built using the variable subsets and the calibration subsets determined by Monte Carlo(MC) re-sampling,and the sub-model that produces minimal error in cross validation is selected as a member model.With repetition of the MC re-sampling,a series of member models are built and a consensus model is achieved by averaging all the member models.Since member models are built with the best variable subset and the randomly selected calibration subset,both the quality and the diversity of the member models are insured for the consensus model.Two NIR spectral datasets of tobacco lamina are used to investigate the proposed method.The superiority of the method in both accuracy and reliability is demonstrated.     

Distribution of counterions in the double layer around a cylindrical polyion

The modified Poisson-Boltzmann (MPB) equation, and the Monte Carlo (MC) method, were applied to the cell model of a polyelectrolyte solution in order to calculate the distribution of counterions around a cylindrical polyion. Both methods suggest stronger binding of counterions to the polyion than predicted by the ordinary Poisson-Boltzmann (PB) equation. The inclusion of counterion-counterion correlation being neglected in the PB equation, leads to a better agreement of the calculated osmotic coefficients with those measured.     

Three-dimensional modeling of EXAFS spectral mixtures by combining Monte Carlo simulations and target transformation factor analysis

We have developed a new method for the three-dimensional modeling of extended X-ray absorption fine structure (EXAFS) spectra which enables the extraction of the local structure of aqueous metal complexes from spectral mixtures of several components. The new method combines two techniques: Monte Carlo simulation and target transformation factor analysis (TFA). Monte Carlo simulation is used to create random arrangements between the X-ray absorbing metal ion and the ligand atoms, and to calculate the theoretical EXAFS spectrum of each arrangement. The theoretical EXAFS spectrum is then introduced as test spectrum in the TFA procedure, to test whether or not the test spectrum is likely to be a component of the spectral mixtures. This coupled procedure is repeated until the error in the test spectrum is minimized. The new method can thus be used to isolate and refine the structure of complexes from spectral mixtures and to determine their relative concentrations, solely on the basis of an estimate of a ligand structure. The performance of the proposed method is validated using uranium Liii-edge EXAFS spectra of binary mixtures of two uranium(VI) 3,4-dihydroxybenzoic acid complexes.     

Mapping the Distribution of Melanin Concentration in Different Fitzpatrick Skin Types Using Hyperspectral Imaging Technique

Measuring skin melanin concentration in order to assess skin phototype according to Fitzpatrick's classification is a constant research goal. In this study, a new approach for assessing skin melanin concentration based on hyperspectral imaging combined with an appropriate analytical model that exploits specific spectral bands to generate maps of melanin content distribution on different Fitzpatrick skin phototypes is presented. Hyperspectral images from the proximal inner side of the forearms of 51 young volunteers covering the first four classes of Fitzpatrick's phototypes were acquired using a hyperspectral imaging system. The images were analyzed using a modified Beer–Lambert law that segregates the contribution of melanin from the other constituents to the skin absorption spectrum. The performance of the model was evaluated using the coefficient of determination (r-squared). The results revealed that the approach proposed in this study generated accurate melanin concentration distribution maps that allowed a correct classification of skin phototype. In conclusion, the proposed approach for assessing skin melanin concentration proved to be very reliable for classifying skin phototypes, and, as it provides maps that are easily read, it has the advantage of a possible extension of its applications to other research concerning skin pigmentation.     

Multiparticle moves in acceptance rate optimized monte carlo

Molecular Dynamics (MD) and Monte Carlo (MC) based simulation methods are widely used to investigate molecular and nanoscale structures and processes. While the investigation of systems in MD simulations is limited by very small time steps, MC methods are often stifled by low acceptance rates for moves that significantly perturb the system. In many Metropolis MC methods with hard potentials, the acceptance rate drops exponentially with the number of uncorrelated, simultaneously proposed moves. In this work, we discuss a multiparticle Acceptance Rate Optimized Monte Carlo approach (AROMoCa) to construct collective moves with near unit acceptance probability, while preserving detailed balance even for large step sizes. After an illustration of the protocol, we demonstrate that AROMoCa significantly accelerates MC simulations in four model systems in comparison to standard MC methods. AROMoCa can be applied to all MC simulations where a gradient of the potential is available and can help to significantly speed up molecular simulations. © 2015 Wiley Periodicals, Inc.     

Noninvasive Assessment of UV-induced Skin Damage: Comparison of Optical Measurements to Histology and MMP Expression

Acute exposure to UV radiation (UVR) causes visible skin damage such as erythema and results in local and systemic immunosuppression while chronic exposure can result in photocarcinogenesis. These deleterious effects can be quantified by histology and by bioassays of key biological markers, including matrix metalloproteinases (MMPs), or tryptophan moieties. We now report our results in quantifying UV skin damage with noninvasive optical methods based on reflectance and fluorescence spectroscopy and compare these noninvasive measurements to histopathology and MMP-13 expression. A solar simulator with spectral output nearly identical to that of solar radiation was developed and used in our experiments. SKH1 hairless mice were exposed to solar-simulated UVR at a total dose of 21 MED delivered over 10 weeks. Changes in oxygenated and deoxygenated hemoglobin were measured by diffuse reflectance spectroscopy, and tryptophan changes were monitored via a fluorescence monitor. Our results show that there is an increase in erythema, skin fluorescence, sunburn cells and MMP-13 after a series of suberythemal doses of UV irradiation on a hairless mouse animal model. Increased skin fluorescence is observed with increasing UV exposure. The levels of MMP-13 increase as the cumulative UV dose increases but their increase does not correspond to noninvasively measured changes.     

Modeling and verification of melanin concentration on human skin type

Lasers are used in the minimalistic or noninvasive diagnosis and treatment of skin disorders. Less laser light reaches the deeper skin layers in dark skin types, due to its higher epidermal melanin concentration compared with lighter skin. Laser–tissue interaction modeling software can correct for this by adapting the dose applied to the skin. This necessitates an easy and reliable method to determine the skin’s type. Noninvasive measurement of the skin’s melanin content is the best method. However, access to samples of all skin types is often limited and skin‐like phantoms are used instead. This study’s objective is to compare experimentally measured absorption features of liquid skin‐like phantoms representing Skin Types I–VI with a realistic skin computational model component of ASAP^®. Sample UV–VIS transmittance spectra were measured from 370 to 900 nm and compared with simulated results from ASAP^® using the same optical parameters. Results indicated nonmonotonic absorption features towards shorter wavelengths, which may allow for more accurate ways of determining melanin concentration and expected absorption through the epidermal layer. This suggests possible use in representing optical characteristics of real skin. However, a more comprehensive model and phantoms are necessary to account for the effects of sun exposure.