Molecular modeling of immersion optical clearing of biological tissues

The interaction of six low-molecular tissue-clearing agents (1,2 and 1,3-propanediol, ethylene glycol, glycerol, xylitol, sorbitol) with the collagen mimetic peptide (GPH)₃ was studied by applying the methods of classical molecular dynamics (GROMACS), molecular docking (AutoDock Vina) and quantum chemistry (PM6 and B3LYP). The spatial configurations of intermolecular complexes were determined and interaction energies calculated. The dependence of the volume occupied by the collagen peptide on the clearing agent concentration in an aqueous solution was calculated. This dependence is not linear, and has a maximum for almost all the agents in the study. The correlations between the optical clearing potential and intermolecular interactions parameters, such as the time of an agent being in a hydrogen-bonded state, and the relative probability of formation of double hydrogen bonds and interaction energies, were determined. Using the correlations determined, we predicted the numeric value of the optical clearing potential of dextrose molecules in rat skin, which correlates with experimental data. A molecular mechanism of tissue optical clearing within the post-diffusion stage is suggested.

Open image in new window

Graphical abstract The molecular modeling of the interaction between clearing agents and collagen

     

Dependence of optical attenuation coefficient and mechanical tension of irradiated human cartilage measured by optical coherence tomography

As banked human tissues are not widely available, the development of new non-destructive and contactless techniques to evaluate the quality of allografts before distribution for transplantation is very important. Also, tissues will be processed accordingly to standard procedures and to minimize disease transmission most tissue banks will include a decontamination or sterilization step such as ionizing radiation. In this work, we present a new method to evaluate the internal structure of frozen or glycerol processed human cartilages, submitted to various dosis of irradiation, using the total optical attenuation coefficient retrieved from optical coherence tomography (OCT) images. Our results show a close relationship between tensile properties and the total optical attenuation coefficient of cartilages. Therefore, OCT associated with the total optical attenuation coefficient open a new window to evaluate quantitatively biological changes in processed tissues.     

Lymph vessels visualization from optical coherence tomography data using depth-resolved attenuation coefficient calculation

Multimodal optical coherent tomography grows popularity with researchers and clinicians over the past decade. One of the modalities is lymphangiography, which allows visualization of the lymphatic vessel networks within optical coherence tomography (OCT) imaging volume. In the present study, it is shown that lymphatic vessel visualization obtained from the depth-resolved attenuation coefficient distributions, corrected for the noise, shows improved contrast and detail in comparison with previously proposed approaches. We also argue that the two most popular approaches for lymphatic vessel visualization, namely simple intensity thresholding and vesselness calculation based on local Hessian matrix eigenvalues, imply different definitions of the lymphatic vessel's appearance in the OCT volume and lead to the different networks.     

Characterizing thrombus with multiple red blood cell compositions by optical coherence tomography attenuation coefficient

Embolectomy is one of the emergency procedures performed to remove emboli. Assessing the composition of human blood clots is an important diagnostic factor and could provide guidance for an appropriate treatment strategy for interventional physicians. Immunostaining has been used to identity compositions of clots as a gold-standard procedure, but it is time-consuming and cannot be performed in situ. Here, we proposed that the optical attenuation coefficient of optical coherence tomography (OCT) can be a reliable indicator as a new imaging modality to differentiate clot compositions. Fifteen human blood clots with multiple red blood cell (RBC) compositions from 21% to 95% were prepared using healthy human whole blood. A homogeneous gelatin phantom experiment and numerical simulation based on the Lambert–Beer's law were examined to verify the validity of the attenuation coefficient estimation. The results displayed that optical attenuation coefficients were strongly correlated with RBC compositions. We reported that attenuation coefficients could be a promising biomarker to guide the choice of an appropriate interventional device in a clinical setting and assist in characterizing blood clots.     

Revisiting the variation of clustering coefficient of biological networks suggests new modular structure

ABSTRACT: A central idea in biology is the hierarchical organization of cellular processes. A commonly used method to identify the hierarchical modular organization of network relies on detecting a global signature known as variation of clustering coefficient (so-called modularity scaling). Although several studies have suggested other possible origins of this signature, it is still widely used nowadays to identify hierarchical modularity, especially in the analysis of biological networks. Therefore, a further and systematical investigation of this signature for different types of biological networks is necessary. RESULTS: We analyzed a variety of biological networks and found that the commonly used signature of hierarchical modularity is actually the reflection of spoke-like topology, suggesting a different view of network architecture. We proved that the existence of super-hubs is the origin that the clustering coefficient of a node follows a particular scaling law with degree k in metabolic networks. To study the modularity of biological networks, we systematically investigated the relationship between repulsion of hubs and variation of clustering coefficient. We provided direct evidences for repulsion between hubs being the underlying origin of the variation of clustering coefficient, and found that for biological networks having no anti-correlation between hubs, such as gene co-expression network, the clustering coefficient doesn't show dependence of degree. CONCLUSIONS: Here we have shown that the variation of clustering coefficient is neither sufficient nor exclusive for a network to be hierarchical. Our results suggest the existence of spoke-like modules as opposed to "deterministic model" of hierarchical modularity, and suggest the need to reconsider the organizational principle of biological hierarchy.     

The effect of particle characteristics on the beam attenuation coefficient and output from an optical backscatter sensor

Experiments were conducted to measure optical backscatter and beam transmission of suspensions of 180, 150 and 90 μm sand, and 40 μm clay, in a recirculation tank designed to house an optical backscatterance sensor (O.B.S.) and a beam transmissometer. Particle size was determined using gravimetric techniques and Coulter counter. By contriving known sediment distributions from the fractionated sediment samples, it was found that both the O.B.S. and beam transmissometer responded approximately linearly to narrow band and broad band particle suspensions. The beam transmissometer showed greater sensitivity to the fine-grain fraction of a poly-disperse suspension than the O.B.S.     

Three‐dimensional mapping of the attenuation coefficient in optical coherence tomography to enhance breast tissue microarchitecture contrast

Effective intraoperative tumor margin assessment is needed to reduce re‐excision rates in breast‐conserving surgery (BCS). Mapping the attenuation coefficient in optical coherence tomography (OCT) throughout a sample to create an image (attenuation imaging) is one promising approach. For the first time, three‐dimensional OCT attenuation imaging of human breast tissue microarchitecture using a wide‐field (up to ~45 × 45 × 3.5 mm) imaging system is demonstrated. Representative results from three mastectomy and one BCS specimen (from 31 specimens) are presented with co‐registered postoperative histology. Attenuation imaging is shown to provide substantially improved contrast over OCT, delineating nuanced features within tumors (including necrosis and variations in tumor cell density and growth patterns) and benign features (such as sclerosing adenosis). Additionally, quantitative micro‐elastography (QME) images presented alongside OCT and attenuation images show that these techniques provide complementary contrast, suggesting that multimodal imaging could increase tissue identification accuracy and potentially improve tumor margin assessment.     

Predicting the light attenuation coefficient through Secchi disk depth and beam attenuation coefficient in a large, shallow, freshwater lake

The diffuse attenuation coefficient of photosynthetically active radiation (PAR) (400–700 nm) (K _d(PAR)) is one of the most important optical properties of water. Our purpose was to create K _d(PAR) prediction models from the Secchi disk depth (SDD) and beam attenuation coefficient of particulate and dissolved organic matter (C _t−w(PAR), excluding pure water) in the PAR range. We compare their performance and prediction precision by using the determination coefficient (r ²), relative root mean square error (RRMSE), and mean relative error (MRE). Our dataset comprised 1,067 measurements, including K _d(PAR), SDD, and C _t−w(PAR) taken in shallow, eutrophic, Lake Taihu, China, from 2005 to 2010. The prediction models of K _d(PAR) were based on the linear model with an intercept of zero, using the inverse SDD, and the nonlinear model using SDD. The linear model generated a slope of 1.369, which was not significantly different from 1.7, the index used worldwide, but significantly lower than the value of 2.26. The nonlinear model gave a slightly more reliable prediction of K _d(PAR) with a r ² of 0.804. Compared to the SDD, C _t−w(PAR) was more significantly correlated to K _d(PAR) based on the linear model, with a significantly higher r ² and lower RMSE and RE. Considering the measurement simplicity of C _t−w(PAR) and data acquisition feasibility from high-frequency autonomous buoys and satellites, our results demonstrated that this prediction model reliably estimates K _d(PAR), and could be used to significantly expand optical observations in an environment where the conditions for underwater PAR measurement are limited.     

K‐distribution three‐dimensional mapping of biological tissues in optical coherence tomography

Probability density function (PDF) analysis with K‐distribution model of optical coherence tomography (OCT) intensity signals has previously yielded a good representation of the average number of scatterers in a coherence volume for microspheres‐in‐water systems, and has shown initial promise for biological tissue characterization. In this work, we extend these previous findings, based on single point M‐mode or two‐dimenstional slice analysis, to full three‐dimensional (3D) imaging maps of the shape parameter α of the K‐distribution PDF. After selecting a suitably sized 3D evaluation window, and verifying methodology in phantoms, the resultant parametric α images obtained in different animal tissues (rat liver and brain) show new contrasting ability not seen in conventional OCT intensity images.      

Image reconstruction of optical computed tomography by using the algebraic reconstruction technique for dose readouts of polymer gel dosimeters

Optical computed tomography (optical CT) has been proven to be a useful tool for dose readouts of polymer gel dosimeters. In this study, the algebraic reconstruction technique (ART) for image reconstruction of gel dosimeters was used to improve the image quality of optical CT. Cylindrical phantoms filled with N-isopropyl-acrylamide polymer gels were irradiated using a medical linear accelerator. A circular dose distribution and a hexagonal dose distribution were produced by applying the VMAT technique and the six-field dose delivery, respectively. The phantoms were scanned using optical CT, and the images were reconstructed using the filtered back-projection (FBP) algorithm and the ART. For the circular dose distribution, the ART successfully reduced the ring artifacts and noise in the reconstructed image. For the hexagonal dose distribution, the ART reduced the hot spots at the entrances of the beams and increased the dose uniformity in the central region. Within 50% isodose line, the gamma pass rates for the 2 mm/3% criteria for the ART and FBP were 99.2% and 88.1%, respectively. The ART could be used for the reconstruction of optical CT images to improve image quality and provide accurate dose conversion for polymer gel dosimeters.     

Determination of linamarin in biological tissues

A new paper chromatography method for determination of intact cyanogenic glucoside linamarin is based on a reaction with p-anisaldehyde at 85°C; the reaction produces a pink color which is brightly fluorescent under long-wave uv light.     

Accuracy and repeatability of an optical motion analysis system for measuring small deformations of biological tissues

Optical motion analysis techniques have been widely used in biomechanics for measuring large-scale motions such as gait, but have not yet been significantly explored for measuring smaller movements such as the tooth displacements under load. In principle, very accurate measurements could be possible and this could provide a valuable tool in many engineering applications. The aim of this study was to evaluate accuracy and repeatability of the Qualisys ProReflex-MCU120 system when measuring small displacements, as a step towards measuring tooth displacements to characterise the properties of the periodontal ligament. Accuracy and repeatability of the system was evaluated using a wedge comparator with a resolution of 0.25 microm to provide measured marker displacements in three orthogonal directions. The marker was moved in ten steps in each direction, for each of seven step sizes (0.5, 1, 2, 3, 5, 10, and 20 microm), repeated five times. Spherical and diamond markers were tested. The system accuracy (i.e. percentage of maximum absolute error in range/measurement range), in the 20-200 microm ranges, was +/-1.17%, +/-1.67% and +/-1.31% for the diamond marker in x, y and z directions, while the system accuracy for the spherical marker was +/-1.81%, +/-2.37% and +/-1.39%. The system repeatability (i.e. maximum standard deviation in the measurement range) measured under the different days, light intensity and temperatures for five times, carried out step up and then step down measurements for the same step size, was +/-1.7, +/-2.3 and +/-1.9 microm for the diamond marker, and +/-2.6, +/-3.9 and +/-1.9 microm for the spherical marker in x, y and z directions, respectively. These results demonstrate that the system suffices accuracy for measuring tooth displacements and could potentially be useful in many other applications.     

Overcoming technical variation and biological variation in quantitative proteomics

Quantitative proteomics investigates physiology at the molecular level by measuring relative differences in protein expression between samples under different experimental conditions. A major obstacle to reliably determining quantitative changes in protein expression is to overcome error imposed by technical variation and biological variation. In drug discovery and development the issue of biological variation often rises in concordance with the developmental stage of research, spanning from in vitro assays to clinical trials. In this paper we present case studies to raise awareness to the issues of technical variation and biological variation and the impact this places on applying quantitative proteomics. We defined the degree of technical variation from the process of two-dimensional electrophoresis as 20-30% coefficient of variation. On the other hand, biological variation observed experiment-to-experiment showed a broader degree of variation depending upon the sample type. This was demonstrated with case studies where variation was monitored across experiments with bacteria, established cell lines, primary cultures, and with drug treated human subjects. We discuss technical variation and biological variation as key factors to consider during experimental design, and offer insight into preparing experiments that overcome this challenge to provide statistically significant outcomes for conducting quantitative proteomic research.     

Ultrasonic characterization of biological tissues

Ultrasonic imaging has become increasingly important as a diagnostic tool in medicine because it is noninvasive and it can provide valuable information otherwise unattainable. However, at present, clinical interpretation of an ultrasonic image still mostly relies on recognition of boundaries and positional relationship of anatomical structures and a subjective analysis of the distribution or texture of echo amplitudes. Other potentially useful information carried back by the echoes is completely discarded. The aim of ultrasonic tissue characterization research is to develop methods to extract additional information from the returned echoes so that tissue pathology or abnormality can be reliably identifed and severity of the pathology objectively assessed with quantitative criteria. A number of ultrasonic parameters including acoustic velocity, impedance, attentuation and scattering, have been utilized in attempting to achieve this goal. In this paper, recent progress in this research will be discussed and relevant results presented.     

Lipid peroxidation products in biological tissues

Over the past twenty-years of lipid peroxidation research in this laboratory, considerable effort has gone into development of new methods, with emphasis on measurement of lipid-soluble fluorophores and the volatile hydrocarbons ethane and pentane. Application of these and other methods has been made to biological materials and living animals. Although the various methodologies used in lipid peroxidation research do not necessarily measure the same class of products, and although agreement of results is not always 100%, there is substantial documentation of good correlations between measurements; for example, of trace volatile hydrocarbons with thiobarbituric acid-reacting substances, of pentane production with dietary and/or tissue vitamin E content, and of pentane production with lipid-soluble fluorophores accumulated in spleen as a function of oxidant stress. Individual methodologies do have their inherent limitations. However, measurements of multiple products and their correlations have added significantly to the base of information on biological damage and protection by dietary antioxidants against nutritional and toxicological insults to tissues, cells, and macromolecules as a result of peroxidative and oxidative reactions.     

Comparison of skewness coefficient,coefficient of variation,and Gini coefficient as inequality measures within populations

Summary The moment skewness coefficient, coefficient of variation and Gini coefficient are contrasted as statistical measures of inequality among members of plant populations. Constructed examples, real data examples, and distributional considerations are used to illustrate pertinent properties of these statistics to assess inequality. All three statistics possess some undesirable properties but these properties are shown to be often unimportant with real data. If the underlying distribution of the variable follows the often assumed two-parameter lognormal model, it is shown that all three statistics are likely to be highly and positively correlated. In contrast, for distributions which are not two-parameter lognormally distributed, and when the distribution is not concentrated near zero, the coefficient of variation and Gini coefficient, which are sensitive to small shifts in the mean, are often of little practical use in ordering the equality of populations. The coefficent of variation is more sensitive to individuals in the right-hand tail of a distribution than is the Gini coefficient. Therefore, the coefficient of variation may often be recommended over the Gini coefficient if a measure of relative precision is selected to assess inequality. The skewness coeficient is suggested when the distribution is either three-parameter lognormally distributed (or close to such), or when a measure of relative precision is not indicated.Scientific Paper no 7830. College of Agriculture and Home Economics Research Center, Washington State University     

Antigenic variation in its biological context

The biology of antigenic variation is discussed, and the problems that must be solved to provide a full understanding of antigenic variation are considered. These are (i) the induction of v.s.g. synthesis in the salivary glands of the tsetse fly; (ii) the nature of the restriction on v.s.g. genes that allows only some of them to be expressed in the salivary glands; (iii) the nature of 'predominance' in v.s.g. expression in the mammalian host, and the mechanism by which it operates; (iv) the repression of v.s.g. synthesis in the insect midgut; (v) the anamnestic response that produces expression of the ingested variant in the first patent parasitaemia in the mammalian host; (vi) the mechanism by which only one v.s.g. gene at a time is expressed; (vii) the relationship if any of v.s.g. structure to v.s.g.-associated differences in growth rate and host range; (viii) the role of v.s.g. release within the life cycle and to pathogenesis.