Monte Carlo simulation of NIR diffuse reflectance in the normal and diseased human breast tissues

The spectral reflectance measurements in tissue reveal physiological meaning. Normally, functional changes like, increase in total hemoglobin concentration, decrease in oxygen saturation, etc., are observed when there is an abnormality creeping in the normal tissue. These functional changes can act together to reveal disease by non-invasive near-infrared (NIR) spectroscopy, as it influence its optical properties. In the present study, a simple two dimensional, four layer model of breast is proposed. The four layers are (i) skin (ii) adipose layer (iii) glandular tissue and (iv) muscle. Each layer is modeled with appropriate biological chromophores like hemoglobin, water, lipid and melanin. From the literature, the concentrations and molar extinction coefficients of the chromophores in various layers of the model are obtained. These values are used to calculate the wavelength dependent absorption characteristics of a particular layer. Monte Carlo simulation of diffuse reflectance (percentage of back reflected photons after multiple scattering with the broad variety of angles) are simulated for the modeled breast tissue with and without diseased condition. Near-infrared wavelengths are chosen, as the depth of penetration in tissue is more compared to UV and visible region. Simulations are carried out on the modeled breast tissue for different races (skin colors) at different NIR wavelengths. Results show significant changes in diffuse reflectance and relative absorbance for normal and diseased breast tissues for differently pigmented model. This model can be used to study the photo dynamical therapy, drug delivery and prognosis of cancer.     

Myocardial oxygenation in isolated hearts predicted by an anatomically realistic microvascular transport model

An anatomically realistic model for oxygen transport in cardiac tissue is introduced for analyzing data measured from isolated perfused guinea pig hearts. The model is constructed to match the microvascular anatomy of cardiac tissue based on available morphometric data. Transport in the three-dimensional system (divided into distinct microvascular, interstitial, and parenchymal spaces) is simulated. The model is used to interpret experimental data on mean cardiac tissue myoglobin saturation and to reveal differences in tissue oxygenation between buffer-perfused and red blood cell-perfused isolated hearts. Interpretation of measured mean myoglobin saturation is strongly dependent on the oxygen content of the perfusate (e.g., red blood cell-containing vs. cell-free perfusate). Model calculations match experimental values of mean tissue myoglobin saturation, measured mean myoglobin, and venous oxygen tension and can be used to predict distributions of intracellular oxygen tension. Calculations reveal that approximately 20% of the tissue is hypoxic with an oxygen tension of <0.5 mmHg when the buffer is equilibrated with 95% oxygen to give an arterial oxygen tension of over 600 mmHg. The addition of red blood cells to give a hematocrit of only 5% prevents tissue hypoxia. It is incorrect to assume that the usual buffer-perfused Langendorff heart preparation is adequately oxygenated for flows in the range of < or =10 ml. min-1. ml tissue-1.     

Rapid Determination of Oxygen Saturation and Vascularity for Cancer Detection

A rapid heuristic ratiometric analysis for estimating tissue hemoglobin concentration and oxygen saturation from measured tissue diffuse reflectance spectra is presented. The analysis was validated in tissue-mimicking phantoms and applied to clinical measurements in head and neck, cervical and breast tissues. The analysis works in two steps. First, a linear equation that translates the ratio of the diffuse reflectance at 584 nm and 545 nm to estimate the tissue hemoglobin concentration using a Monte Carlo-based lookup table was developed. This equation is independent of tissue scattering and oxygen saturation. Second, the oxygen saturation was estimated using non-linear logistic equations that translate the ratio of the diffuse reflectance spectra at 539 nm to 545 nm into the tissue oxygen saturation. Correlations coefficients of 0.89 (0.86), 0.77 (0.71) and 0.69 (0.43) were obtained for the tissue hemoglobin concentration (oxygen saturation) values extracted using the full spectral Monte Carlo and the ratiometric analysis, for clinical measurements in head and neck, breast and cervical tissues, respectively. The ratiometric analysis was more than 4000 times faster than the inverse Monte Carlo analysis for estimating tissue hemoglobin concentration and oxygen saturation in simulated phantom experiments. In addition, the discriminatory power of the two analyses was similar. These results show the potential of such empirical tools to rapidly estimate tissue hemoglobin in real-time spectral imaging applications.     

Mapping the myoglobin concentration, oxygenation, and optical pathlength in heart ex vivo using near-infrared imaging

A method that provides maps of absolute concentrations of oxygenated and deoxygenated myoglobin (Mb), its oxygenation, and its near-infrared (NIR) optical pathlength in cardiac tissue was developed. These parameters are available simultaneously. The method is based on NIR diffuse reflectance spectroscopic imaging and specific processing of the NIR images, which included a first derivative of the diffuse reflectance spectrum. Mb oxygenation, total Mb concentration, and NIR light pathlength were found to be in the range of 92%, 0.3 mM, and 12.5 mm, respectively, in beating isolated buffer-perfused and arrested pig hearts. The charge-coupled device camera enables sub-millimeter spatial resolution and spectroscopic imaging in 1.5 to 2.0 min. The technique is noninvasive and nondestructive. The equipment has no mechanical contact with the tissue of interest, leaving it undisturbed.     

Hemoglobin plus myoglobin concentrations and near infrared light pathlength in phantom and pig hearts determined by diffuse reflectance spectroscopy

To noninvasively determine absolute concentrations of hemoglobin (Hb) plus myoglobin (Mb) in cardiac tissue by means of regular near infrared (NIR) light diffuse reflectance measurements, a first derivative approach was applied. The method was developed to separately calculate oxygenated and deoxygenated [Hb + Mb] as well as an effective pathlength, which NIR light passes through in the tissue between optodes. Applying a cotton wool-based phantom, which mimics muscle tissue, it was shown that the intensity of the pseudo-optical density first derivative depends linearly on both oxygenated and deoxygenated Hb concentration, thereby validating the Lambert-Beer law in the range of 0 to 0.25 mM tetrameric Hb. A high correlation (R = 0.995) was found between concentrations of Hb loaded onto the phantom and those determined spectrophotometrically, thereby verifying the first derivative method validity. The efficiency of the method was tested using in vivo pig hearts prior to and after ischemia initiated experimentally by left anterior descending artery branches occlusion. The results showed that the total [Hb + Mb] was 0.9-1.2 mM heme, the average tissue oxygen saturation was approximately 70% (which reduced to nearly 0% after occlusion), and the NIR (700-965 nm) light pathlength was 2.3 mm (differential pathlength factor [DPF] = 2.7-2.8) in a living heart tissue.     

Photosensitization of aqueous model systems by hypericin.

Absorption and fluorescence measurements of purified hypericin (HY) were made in various media. Photosensitization of two aqueous systems was investigated: resealed red blood cell membranes (ghosts) and hen lysozyme (Lys). Solubilization of HY by ghost membranes was shown by means of diffuse reflectance spectroscopy. Visible light irradiation of the ghosts incorporating HY led to lipid peroxidation with evidence of singlet oxygen involvement. A binding model applicable for insoluble ligands is indicative of strong HY binding to HSA. The HY-HSA complex photosensitized inactivation of Lys. The pseudo-first-order reaction kinetics with protection by azide ion are consistent with a Type II mechanism mediated by singlet oxygen. The results are discussed in the context of the HY photodynamic and antiretroviral activities.     

Human head dynamic response to side impact by finite element modeling

The dynamic response of the human head to side impact was studied by 2-dimensional finite element modeling. Three models were formulated in this study. Model I is an axisymmetric model. It simulated closed shell impact of the human head, and consisted of a single-layered spherical shell filled wiht an inviscid fluid. The other two models (Model II and III) are plane strain models of a coronal section of the human head. Model II approximated a 50th percentile male head by an outer layer to simulate cranial bone and an inviscid interior core to simulate the intracranial contents. The configuration of Model III is the same as Model II but more detailed anatomical features of the head interior were added, such as, cerebral spinal fluid (CSF); falx cerebri, dura, and tentorium. Linear elastic material properties were assigned to all three models. All three models were loaded by a triangular pulse with a peak pressure of 40 kPa, effectively producing a peak force of 1954 N (440 lb). The purpose of this study was to determine the effects of the membranes and that of the mechanical properties of the skull, brain, and membrane on the dynamic response of the brain during side impact, and to compare the pressure distributions from the plane strain model with the axisymmetric model. A parametric study was conducted on Model II to characterize fully its response to impact under various conditions.(ABSTRACT TRUNCATED AT 250 WORDS)     

生物组织漫射光模拟计算的改进杂合模型

结合近源处Ｍｏｎｔｅ Ｃａｒｌｏ法的精确性和漫射理论计算的快速性,在用漫射理论计算漫反射系数时,采用了辐射对背半球的积分形式,即考虑光强度和通量两项对它的贡献,建立了基于Ｍｏｎｔｅ Ｃａｒｌｏ法－漫理论的更为精确、有效的改进的杂合模型,模拟计算了板状生物组织漫反射光的分布。结果表明,在保证高计算速度的前提下,进一步提高了计算精度;同时,此模型也适用于强吸收介质漫反射光分布的模拟计算。     

Quantitation of time- and frequency-resolved optical spectra for the determination of tissue oxygenation.

The recent development of near-infrared time- and frequency-resolved tissue spectroscopy techniques to probe tissue oxygenation and tissue oxygenation kinetics has led to the need for further quantitation of spectroscopic signals. In this paper, we briefly review the theory of light transport in strongly scattering media as monitored in the time and frequency domains, and use this theory to develop algorithms for quantitation of hemoglobin saturation from the photon decay rate (delta log R/delta t) obtained using time-resolved spectroscopy, and from the phase-shift (theta) obtained from frequency-resolved, phase-modulated spectroscopy. To test the relationship of these optical parameters, we studied the behavior of delta log R/delta t and theta as a function of oxygenation in model systems which mimicked the optical properties of tissue. Our results show that deoxygenation at varying hemoglobin concentrations can be monitored with the change in the photon decay kinetics, delta delta log R/delta t in the time-resolved measurements, and with the change in phase-shift, delta theta, in the frequency-resolved technique. Optical spectra of the adult human brain obtained with these two techniques show similar characteristics identified from the model systems.     

Binding of cyanide,cyanate, and thiocyanate to human carbonic anhydrase II

Computer simulation techniques are used to address the question of how cyanide and related ions interact with human carbonic anhydrase II (HCAII). Spectroscopic results have suggested that cyanide is coordinated with the zinc ion, while recent X-ray results suggest that the cyanide ion is noncovalently associated with the zinc–water or zinchydroxide form of the enzyme. We have carried out simulations on three models in an attempt to shed light on why the spectroscopic and X-ray results differ. The first model we studied (Model I) has cyanide directly coordinated to the zinc ion, the second has it noncovalently interacting with the zinc–hydroxide (high pH) form of the enzyme (Model II), and the third has cyanide noncovalently interacting with the zinc–water (low pH) form of the enzyme (Model III). None of these models is satisfactory in explaining the available structural data obtained from X-ray crystallography. This leads us to propose an alternative model, in which HCAII hydrates HCN to form an OH^?/HCN complex coordinated to the Zn ion. Ab initio calculations are consistent with this model. Based on these results we are able to explain the observed crystallographic behavior of cyanate and, by inference, thiocyanate. © 1993 Wiley-Liss, Inc.     

用蒙特卡罗方法模拟光在多层组织中的吸收特性

在讨论目前新颖的组织功能成像打骂能性（例如光声成像）时,光子在组织中的吸收和散射特性是一个很重要的问题,鉴于这一点,本文利用一个多层模型研究了光子在皮肤,脂肪和肌肉组织中的吸收和散射特性,得到了在组织中某一深度处光子在一个平面上的吸收分布,以及在不同吸收系数和散射系数的情况下,光子的反射,吸收和透射几率,结果表明在经过多次散射后,大部分的光子被吸收,在本文的模型中只有7.3%的光子从表面反射（包括镜面反射和漫反射）,还讨论了不同光学参灵敏对参流分布的影响。     

Incorporating rhizosphere processes into field-scale (agro)ecosystem models

Three different strategies for incorporating rhizosphere processes within field-scale models are compared, taking triple-cropped irrigated rice production as a common system and CH₄ emission as a common focus of interest. The strategies may be characterised as homogeneous (model I; root C deposition is added to the bulk soil compartment), areal (model II; roots contribute via aerenchymatous exchange to an increased soil–atmosphere interfacial surface area), and volumetric (model III; roots create around themselves a specific rhizosphere compartment). Model I is simpler than model II, which is simpler than model III. With identical parameters all models lead to similar seasonally integrated CH₄ emissions, but when the pattern of emission and the simulated CH₄ concentration in the soil is brought into the reckoning, the following order of precedence (greater is better) becomes clear: model IIImodel II>model I. Current field-scale models of soil organic matter (SOM) transformation, especially in rice soils, could be improved by taking explicit account of the rhizosphere and the processes which occur within it.     

用于胃癌检测的漫反射光谱

胃癌正严重地威胁人类的健康,为了开发一种可以诊断早期胃癌的新型光学检测技术,建立了一套简便的组织体漫反射光谱检测装置.首先,根据组织体的光学特性,介绍光谱检测方法的基本原理.然后,利用组织癌变导致的组织体光学特性变化,建立漫反射光谱检测的实验装置.最后,利用该实验装置分别获取正常胃组织和胃癌组织的漫反射光谱.实验结果表明:癌变和正常胃组织的漫反射光谱在可见光区域有明显差别,特别是在波长为500 nm和630 nm处,但是漫反射光谱检测无法分辨来自组织不同层的光谱信息.这些结果说明漫反射光谱检测装置可用于胃癌的辅助检测.     

Investigation of the source‐detector separation in near infrared spectroscopy for healthy and clinical applications

Understanding near infrared light propagation in tissue is vital for designing next generation optical brain imaging devices. Monte Carlo (MC) simulations provide a controlled mechanism to characterize and evaluate contributions of diverse near infrared spectroscopy (NIRS) sensor configurations and parameters. In this study, we developed a multilayer adult digital head model under both healthy and clinical settings and assessed light‐tissue interaction through MC simulations in terms of partial differential pathlength, mean total optical pathlength, diffuse reflectance, detector light intensity and spatial sensitivity profile of optical measurements. The model incorporated four layers: scalp, skull, cerebrospinal‐fluid and cerebral cortex with and without a customizable lesion for modeling hematoma of different sizes and depths. The effect of source‐detector separation (SDS) on optical measurements' sensitivity to brain tissue was investigated. Results from 1330 separate simulations [(4 lesion volumes × 4 lesion depths for clinical +3 healthy settings) × 7 SDS × 10 simulation = 1330)] each with 100 million photons indicated that selection of SDS is critical to acquire optimal measurements from the brain and recommended SDS to be 25 to 35 mm depending on the wavelengths to obtain optical monitoring of the adult brain function. The findings here can guide the design of future NIRS probes for functional neuroimaging and clinical diagnostic systems.      

Heterogeneous organ models

A theoretical study is made of three organ flow models with heterogeneity of capillary transit times. A new parametrization of Rose and Goresky's Model III facilitates in many cases a reduction to Goresky's Model II, accomplished by a special time shift. The shift parameter defined here is critical in this analysis of Model III. A new expression of the series for outflow concentration in Model III is given and proves useful in examining the model as an operator and in relating it to Models I and II. A result on parameter optimization is given: if then Model III cannot fit better than Model II. This is applied to some data from Rose and Goresky [Circulation Res. 39, 541–544 (1976)] and raises a new question about their model. A heart model of Levin and Bassingthwaighte based on regional flow measurement is shown to be a discretized generalization of Model II. This work supported in part by PHS Grant Nos. HL-19153 (SCOR for Pulmonary Vascular Disease) and HL-19370 at Vanderbilt University.     

Quinohemoprotein alcohol dehydrogenases: structure, function, and physiology

Quino(hemo)protein alcohol dehydrogenases (ADH) that have pyrroloquinoline quinone (PQQ) as the prosthetic group are classified into 3 groups, types I, II, and III. Type I ADH is a simple quinoprotein having PQQ as the only prosthetic group, while type II and type III ADHs are quinohemoprotein having heme c as well as PQQ in the catalytic polypeptide. Type II ADH is a soluble periplasmic enzyme and is widely distributed in Proteobacteria such as Pseudomonas, Ralstonia, Comamonas, etc. In contrast, type III ADH is a membrane-bound enzyme working on the periplasmic surface solely in acetic acid bacteria. It consists of three subunits that comprise a quinohemoprotein catalytic subunit, a triheme cytochrome c subunit, and a third subunit of unknown function. The catalytic subunits of all the quino(hemo)protein ADHs have a common structural motif, a quinoprotein-specific superbarrel domain, where PQQ is deeply embedded in the center. In addition, in the type II and type III ADHs this subunit contains a unique heme c domain. Various type II ADHs each have a unique substrate specificity, accepting a wide variety of alcohols, as is discussed on the basis of recent X-ray crystallographic analyses. Electron transfer within both type II and III ADHs is discussed in terms of the intramolecular reaction from PQQ to heme c and also from heme to heme, and in terms of the intermolecular reaction with azurin and ubiquinone, respectively. Unique physiological functions of both types of quinohemoprotein ADHs are also discussed.     

Single-crystal electron paramagnetic resonance studies of photolyzed oxy- and nitric oxide-cobalt myoglobins

Low-temperature photodissociation of oxygen from oxy-cobalt myoglobin was studied by single-crystal electron paramagnetic resonance (EPR) spectroscopy at 5 K. The photolyzed oxy-cobalt myoglobin exhibited an EPR spectrum consisting of two nonequivalent sets (species I and II) of the principal values and eigenvectors of the g tensors: g1I = 3.55, g2I = 3.47, and g3I = 2.26 for species I, and g1II = 2.04, g2II = 1.93, and g3II = 1.86 for species II, which resembled neither the deoxy nor the oxy form. Possible models of the photodissociated state of oxy-cobalt myoglobin are proposed by comparison with cobalt porphyrin complexes. The photolyzed product of nitric oxide-cobalt myoglobin exhibited new EPR signals at g = 4.3 and a very broad signal at around g = 2. The principal g values have been determined from the single-crystal EPR measurements: g1 = 4.39, g2 = 4.27, and g3 = 4.00. Analysis of another EPR signal around g = 2 was difficult due to its broadness. Magnetic interactions were observed. An isotropic EPR signal at g = 4.3 suggested a weakly spin-coupled system between cobaltous spin (S = 1/2 or 3/2) and nitric oxide spin (S = 1/2).     

Light absorption, electron paramagnetic resonance and resonance Raman characteristics of nitridochromium(V) protoporphyrin-IX and its reconstituted hemoproteins.

A surprisingly stable complex of the photolyzed product of azidochromium(III)protoporphyrin-IX was prepared and examined by light absorption, electron paramagnetic resonance (EPR) and resonance Raman spectroscopies. The characteristic EPR spectrum for this complex was consistent with a nitridochromium(V)-porphyrin complex which was two oxidation equivalents above the resting Cr(III) complex. The Cr(V)-N stretching mode was observed at 1010 cm-1 by resonance Raman spectroscopy. A simple diatomic harmonic oscillation model gave a force constant of 6.7 mdyn/A for the Cr(V)-N bond, in the region characteristic for the metal-nitrogen triple bond. Nitridochromium(V) protoporphyrin-IX reconstituted myoglobin and cytochrome c peroxidase were prepared for the first time. The nitridochromium(V)-porphyrins in these apo-proteins were unstable in contrast with the protein-free chromium(V)porphyrin. Upon irradiation of the azide complexes of the chromium(III) protoporphyrin-IX reconstituted myoglobin and cytochrome c peroxidase with ultraviolet light aerobically at room temperature, the characteristic optical and EPR spectra for nitridochromium(V) derivatives were observed. The optical spectra of these photo-induced products were different from those of the nitridochromium(V) protoporphyrin-IX reconstituted hemoproteins. The electrochemical structures of the unusual metalloporphyrin seemed to be modulated by the heme surrounding amino acid residues.     

Analysis of Fungal Pellets by UV-Visible Spectrum Diffuse Reflectance Spectroscopy

The application of the UV-visible spectrum diffuse reflectance spectroscopy for the determination of intracellular pH in vivo, for determination of cytochrome content, and for the noninvasive in vivo detection of the redox state of fungal mitochondrial cytochromes in filamentous fungi is introduced. The time course of the intracellular pH values, mitochondrial cytochromes, and CO-binding pigments content and the correlations between the actual redox state of cytochrome aa(3) and saturation of growth medium with oxygen in pellets of the basidiomycete Phanerochaete chrysosporium were determined. As the test microorganism, the yeast Saccharomyces cerevisiae was used. UV-visible spectrum diffuse reflectance spectroscopy proved to be a promising method for the quick and simple analysis of light-impermeable biological structures for which the classical transmittance spectrophotometric methods are difficult to implement.     

On selection of functional response models: Holling’s models and more

Model selection is a common and established research method. Statistically rigorous model selection methods are used in a variety of research fields. In contrast, studies that characterize functional response models commonly use a model selection method that is specific to functional response studies. The specific method aims to distinguish between Holling’s type II and type III functional response models. This paper discusses problems associated with the specific method and suggests that it would be better to use general model selection methods that allow to consider a variety of models.