Simulation of ultrasound backscattering by red cell aggregates: effect of shear rate and anisotropy

Tissue characterization using ultrasound (US) scattering allows extraction of relevant cellular biophysical information noninvasively. Characterization of the level of red blood cell (RBC) aggregation is one of the proposed application. In the current paper, it is hypothesized that the microstructure of the RBCs is a main determinant of the US backscattered power. A simulation model was developed to study the effect of various RBC configurations on the backscattered power. It is an iterative dynamical model that considers the effect of the adhesive and repulsive forces between RBCs, and the effect of the flow. The method is shown to be efficient to model polydispersity in size, shape, and orientation of the aggregates due to the flow, and to relate these variations to the US backscattering properties. Three levels of aggregability at shear rates varying between 0.05 and 10 s(-1) were modeled at 40% hematocrit. The simulated backscattered power increased with a decrease in the shear rate or an increase in the RBC aggregability. Angular dependence of the backscattered power was observed. It is the first attempt to model the US power backscattered by RBC aggregates polydisperse in size and shape due to the shearing of the flow.     

生物组织散射元平均间距估计的一种新方法

生物组织散射元平均间中划描述生物组织微观结构特性和生物组织超微散射特性的重要参数。本文在对生物组织超声背向散射随机模的基础上,提出了基于生物组织超声背向散射信号突变点检测的工用射元平均间距估计的新方法。该方法是生物组织超声散射分析的有效方法。     

Dynamic evaluation of blood flow microcirculation by combined use of the laser Doppler flowmetry and high‐speed videocapillaroscopy methods

The dynamic light scattering methods are widely used in biomedical diagnostics involving evaluation of blood flow. However, there exist some difficulties in quantitative interpretation of backscattered light signals from the viewpoint of diagnostic information. This study considers the application of the high‐speed videocapillaroscopy (VCS) method that provides the direct measurement of the red blood cells (RBCs) velocity into a capillary. The VCS signal presents true oscillation nature of backscattered light caused by moving RBCs. Thus, the VCS signal can be assigned as a reference one with respect to more complicated signals like in laser Doppler flowmetry (LDF). An essential correlation between blood flow velocity oscillations in a separate human capillary and the integral perfusion estimate obtained by the LDF method has been found. The observation of blood flow by the VCS method during upper arm occlusion has shown emergence of the reverse blood flow effect in capillaries that corresponds to the biological zero signal in the LDF. The reverse blood flow effect has to be taken into account in interpretation of LDF signals.      

Spatial and spectral resolution of the plasma Doppler reflectometry

An approach to numerical simulations of microwave scattering in Doppler reflectometry is developed that is based on calculations of the two-dimensional spatial weighting function in the Born approximation. The knowledge of the spatial weighting (or instrumental) function, which contains complete information on both forward and backward scattering processes, allows one to calculate the output signal of a quadrature diagnostic detector for any given distribution of the electron plasma density fluctuations. The weighting functions were computed for axisymmetric distributions of the background plasma density and for arbitrarily specified fields in the mouths of the emitting and receiving antennas. Simulations carried out for different model representations of the scattering fluctuations yielded quantitative estimates of both the spatial and wavenumber resolutions of Doppler reflectometry. The simulations showed that the resolution of the method proposed worsens with increasing curvature of the cut-off surface. Conditions are determined under which the resolution of the diagnostics in small-and medium-size tokamaks can be substantially improved by using converging microwave beams. The possibility of utilizing focused probing microwave beams in Doppler reflectometry on ITER-scale devices is discussed.     

Size determination of red blood cell aggregates induced by dextran using ultrasound backscattering phenomenon.

Different methods are commonly used to study the red blood cell aggregation phenomenon. The major interest of the ultrasonic method presently discussed is to assess the mean size of red blood cell (RBC) aggregates by measuring ultrasonic intensity backscattered by blood. Applying Rayleigh theory of sound to blood medium, one can show that the scattered ultrasonic intensity is proportional to the 6th power of the size of the RBC aggregates. The ultrasonic method is used to evaluate the mean size of RBC aggregates induced by dextrans. RBCs are suspended at various hematocrits H, in solution of dextrans of different molecular weights M and at different weight concentrations Cw. Results are presented by using the ultrasonic backscattering coefficient chi which is a relevant quantity in a scattering experiment. For suspensions of RBCs aggregated with dextran of molecular weight 70,000 dalton (dextran 70) at concentration Cw = 40 g/l, variations of chi as a function of H are similar to those obtained for normal blood. At a fixed hematocrit, variation of chi versus Cw for dextran 70 exhibits a maximum at 40 g/l. In the case of RBCs suspended at hematocrit 20% and aggregated with dextrans of molecular weight M, 70,000 less than or equal to M less than or equal to 2,000,000, the variations of chi versus molar concentration Cm are similar to those of the microscopic aggregation index defined by Chien (1). Finally, a statistical model of the blood structure previously described (2) is applied to evaluate the mean size of the aggregates. According to this model, the mean size of aggregates is independent of hematocrit for H less than or equal to 40% and independent of the molecular weight of dextran for M greater than or equal to 150,000 dalton.     

超声彩色血流成像的计算机快速仿真方法

研究超声彩色血流成像的快速仿真方法,克服原先仿真方法非常耗时的缺点。方法超声彩色血流成像计算机仿真中,血流信号是对成像区间内所有点散射体的回波信号累加而得到的。通过引入新的等效散射体模型,可以大大降低散射体的密度,从而减少计算回波信号所需时间。在计算机上用Matlab编程来进行仿真实验,对以往仿真方法和基于等效散射体模型方法的性能进行比较。结果实验表明:基于等效散射体模型的仿真,在保证相同流速精度的前提下,仿真速度比传统方法提高了10倍以上。结论基于等效散射体模型的仿真方法能极大地提高超声彩色血流成像的仿真速度,可以为超声彩色血流成像的方法研究提供便利。     

Quasi-elastic light scattering studies of membrane motion in single red blood cells.

Studies of red blood cells (RBCs) and RBC ghosts, using a quasi-elastic light scattering (QELS) microscope spectrometer, have identified the membrane as the primary source of the light scattering signal. This is the first report in which motion of the cell membrane has been demonstrated to be the primary source of the QELS signal from a cell. Cytoplasmic changes induced in the RBC by varying the osmotic strength of the medium were also detected using this technique. Comparison of the data from white blood cells (WBCs) with the RBC data demonstrated significant differences between different types of cells.     

Size Distribution of Air Bubbles Entering the Brain during Cardiac Surgery

Background

Thousands of air bubbles enter the cerebral circulation during cardiac surgery, but whether high numbers of bubbles explain post-operative cognitive decline is currently controversial. This study estimates the size distribution of air bubbles and volume of air entering the cerebral arteries intra-operatively based on analysis of transcranial Doppler ultrasound data.

Methods

Transcranial Doppler ultrasound recordings from ten patients undergoing heart surgery were analysed for the presence of embolic signals. The backscattered intensity of each embolic signal was modelled based on ultrasound scattering theory to provide an estimate of bubble diameter. The impact of showers of bubbles on cerebral blood-flow was then investigated using patient-specific Monte-Carlo simulations to model the accumulation and clearance of bubbles within a model vasculature.

Results

Analysis of Doppler ultrasound recordings revealed a minimum of 371 and maximum of 6476 bubbles entering the middle cerebral artery territories during surgery. This was estimated to correspond to a total volume of air ranging between 0.003 and 0.12 mL. Based on analysis of a total of 18667 embolic signals, the median diameter of bubbles entering the cerebral arteries was 33 μm (IQR: 18 to 69 μm). Although bubble diameters ranged from ~5 μm to 3.5 mm, the majority (85%) were less than 100 μm. Numerous small bubbles detected during cardiopulmonary bypass were estimated by Monte-Carlo simulation to be benign. However, during weaning from bypass, showers containing large macro-bubbles were observed, which were estimated to transiently affect up to 2.2% of arterioles.

Conclusions

Detailed analysis of Doppler ultrasound data can be used to provide an estimate of bubble diameter, total volume of air, and the likely impact of embolic showers on cerebral blood flow. Although bubbles are alarmingly numerous during surgery, our simulations suggest that the majority of bubbles are too small to be harmful.     

Simulations of NMR-detected diffusion in suspensions of red cells: the "signatures" in q-space plots of various lattice arrangements

Coherence effects from pulsed field-gradient spin-echo (PGSE) nuclear magnetic resonance diffusion experiments have been observed and characterized for diffusants in many heterogeneous systems, ranging from porous materials to cell suspensions. The resulting coherence patterns appear in plots of the normalized PGSE signal intensities as a function of the spatial wave vector Q in a so-called q-space plot. The origin of these phenomena and their mathematical and physical underpinnings are now well established. We have conducted a number of studies of diffusion-coherence phenomena in suspensions of red blood cells and have made extensive use of computer simulations of molecular diffusion in virtual lattices of cells to aid in the interpretation and analysis of experimental data. In the current work we extended the canonical model used in these studies to investigate the effect that varying the packing arrangement of cells in the suspension has on the coherence patterns, as seen in q-space plots. We show that changes in the packing arrangement of cells are reflected in the q-space plots and in the results of diffusion tensor analysis and thus we speculate upon the possible clinical importance of these findings.     

Advanced spectral estimators for detailed blood flow studies.

Recent publications have emphasized the relationship between the spectrum of the backscattered acoustic signal, beam geometry, and flow patterns in the measurement of blood flow by Doppler ultrasound. On this basis, we believe that in the future more importance will be placed on analyzing various characteristics of the spectral shape rather than absolute parameters of measurement, such as the mean frequency. The potential of this approach for extracting more information from the raw Doppler signal is introduced by considering the Spectral Broadening Index (SBI). We explain the use of the SBI parameter for measuring flow angle under restricted flow conditions. This is done by using an analytic/computational model for prediction of the spectral broadening effect. By simulation study, the performance of various spectral estimators for determining the SBI from finite Doppler signal segments is evaluated.     

A two-compartment model of osmotic lysis in Plasmodium falciparum-infected erythrocytes

We recently identified a voltage-dependent anion channel on the surface of human red blood cells (RBCs) infected with the malaria parasite, Plasmodium falciparum. This channel, the plasmodial erythrocyte surface anion channel (PESAC), likely accounts for the increased permeability of infected RBCs to various small solutes, as assessed quantitatively with radioisotope flux and patch-clamp studies. Whereas this increased permeability has also been studied by following osmotic lysis of infected cells in permeant solutes, these experiments have been limited to qualitative comparisons of lysis rates. To permit more quantitative examination of lysis rates, we have developed a mathematical model for osmotic fragility of infected cells based on diffusional uptake via PESAC and the two-compartment geometry of infected RBCs. This model, combined with a simple light scattering assay designed to track osmotic lysis precisely, produced permeability coefficients that match both previous isotope flux and patch-clamp estimates. Our model and light scattering assay also revealed Michaelian kinetics for inhibition of PESAC by furosemide, suggesting a 1:1 stoichiometry for their interaction.     

Inflow/Outflow Boundary Conditions for Particle-Based Blood Flow Simulations: Application to Arterial Bifurcations and Trees

When blood flows through a bifurcation, red blood cells (RBCs) travel into side branches at different hematocrit levels, and it is even possible that all RBCs enter into one branch only, leading to a complete separation of plasma and RBCs. To quantify this phenomenon via particle-based mesoscopic simulations, we developed a general framework for open boundary conditions in multiphase flows that is effective even for high hematocrit levels. The inflow at the inlet is duplicated from a fully developed flow generated in a pilot simulation with periodic boundary conditions. The outflow is controlled by adaptive forces to maintain the flow rate and velocity gradient at fixed values, while the particles leaving the arteriole at the outlet are removed from the system. Upon validation of this approach, we performed systematic 3D simulations to study plasma skimming in arterioles of diameters 20 to 32 microns. For a flow rate ratio 6:1 at the branches, we observed the “all-or-nothing” phenomenon with plasma only entering the low flow rate branch. We then simulated blood-plasma separation in arteriolar bifurcations with different bifurcation angles and same diameter of the daughter branches. Our simulations predict a significant increase in RBC flux through the main daughter branch as the bifurcation angle is increased. Finally, we demonstrated the effectiveness of the new methodology in simulations of blood flow in vessels with multiple inlets and outlets, constructed using an angiogenesis model.     

A methodology for estimating the shape of biconcave red blood cells using multicolor scattering images

In this paper, a novel methodology for estimating the shape of human biconcave red blood cells (RBCs), using color scattering images, is presented. The information retrieval process includes, image normalization, features extraction using two-dimensional discrete transforms, such as angular radial transform (ART), Zernike moments and Gabor filters bank and features dimension reduction using both independent component analysis (ICA) and principal component analysis (PCA). A radial basis neural network (RBF-NN) estimates the RBC geometrical properties. The proposed method is evaluated in both regression and identification tasks by processing images of a simulated device used to acquire scattering phenomena of moving RBCs. The simulated device consists of a tricolor light source (light emitting diode – LED) and moving RBCs in a thin glass. The evaluation database includes 23,625 scattering images, obtained by means of the boundary element method. The regression and identification accuracy of the actual RBC shape is estimated using three feature sets in the presence of additive white Gaussian noise from 60 to 10 dB SNR and systematic distortion, giving a mean error rate less than 1% of the actual RBC shape, and more than 99% mean identification rate.     

Spatial arrangement of DNA-dependent RNA polymerase of Escherichia coli and DNA in the specific complex. A neutron small angle scattering study

In this paper we demonstrate that neutron small angle scattering is a suitable method to study the spatial arrangement of large specific protein-DNA complexes. We studied the complex of DNA-dependent RNA polymerase of Escherichia coli and a 130 base-pair DNA fragment containing the strong promoter A1 of bacteriophage T7. Contrast variation of the complex with deuterium allowed us to "visualize" either RNA polymerase, or DNA, or both components in situ. From the corresponding scattering curves information was derived about: (1) Conformational changes of RNA polymerase and DNA by complex formation: comparison of the scattering profiles of the isolated and complexed components showed that by specific complex formation the cross-section of RNA polymerase decreases, while the DNA fragment does not undergo a gross conformational change. (2) The spatial arrangement of RNA polymerase and DNA in the specific complex from the cross-sectional radii of gyration of the complex the normal distance dn between the centre of gravity of the RNA polymerase and the axis of the DNA fragment was derived as 5.0 (+/- 0.3) nm. On the basis of these and footprinting data a low resolution model of the RNA polymerase-promoter complex is proposed. The main feature of this model is the positioning of RNA polymerase to only one side of the DNA.     

Spectra of enhanced scattering by spontaneous density fluctuations in a tokamak

The wavenumber-resolved radar backscattering diagnostics in the upper hybrid resonance (UHR) region was used to study low-frequency short-scale turbulence in the FT-1 tokamak. The scattered spectra were measured for different delay times of the scattered signals. It is shown that the width of the spectrum of enhanced scattering by spontaneous fluctuations is proportional to the delay time. Possible mechanisms for the formation of the scattered spectra are proposed and discussed. The results of simulations and additional experiments were used to determine the dominant mechanisms governing the formation of the scattered spectra in the FT-1 tokamak. These mechanisms are related to the effect of multiple small-angle scattering of both the probing wave and the waves backscattered in the UHR region by long-scale density fluctuations and to the Doppler effect caused by the entrainment of short-scale fluctuations by the long-scale turbulent flow.     

Regulation of phosphatidylserine exposure in red blood cells

The exposure of phosphatidylserine (PS) on the outer membrane leaflet of red blood cells (RBCs) serves as a signal for eryptosis, a mechanism for the RBC clearance from blood circulation. The process of PS exposure was investigated as function of the intracellular Ca(2+) content and the activation of PKCα in human and sheep RBCs. Cells were treated with lysophosphatidic acid (LPA), 4-bromo-A23187, or phorbol-12 myristate-13 acetate (PMA) and analysed by flow cytometry, single cell fluorescence video imaging, or confocal microscopy. For human RBCs, no clear correlation existed between the number of cells with an elevated Ca(2+) content and PS exposure. Results are explained by three different mechanisms responsible for the PS exposure in human RBCs: (i) Ca(2+)-stimulated scramblase activation (and flippase inhibition) by LPA, 4-bromo-A23187, and PMA; (ii) PKC activation by LPA and PMA; and (iii) enhanced lipid flop caused by LPA. In sheep RBCs, only the latter mechanism occurs suggesting absence of scramblase activity.     

生物组织光散射等效颗粒模型及Mie相函数计算

为研究生物组织的散射特性,将其从散射效果上等效为离散的球形散射体的集合,结合经典的Mie散射理论对生物组织散射相函数进行数值计算。计算结果表明:Mie散射相函数能够描述生物组织后向(大角度)散射光强振荡特性与等效粒径的对应关系,可为基于后向散射光的无创伤或微创伤诊疗提供理论依据;Mie散射相函数能够解释生物组织散射光空间分布与波长的相关性,为医学诊疗上入射光波长选择提供参考;合理选择集群散射体粒度分布参数,可实现对复杂生物组织散射相函数的精确描述。     

Temporal adaptability and the inverse relationship to sensitivity: a parameter identification model

Following a prolonged period of visual adaptation to a temporally modulated sinusoidal luminance pattern, the threshold contrast of a similar visual pattern is elevated. The adaptive elevation in threshold contrast is selective for spatial frequency, may saturate at low adaptor contrast, and increases as a function of the spatio-temporal frequency of the adapting signal. A model for signal extraction that is capable of explaining these threshold contrast effects of adaptation is proposed. Contrast adaptation in the model is explained by the identification of the parameters of an environmental model: the autocorrelation function of the visualized signal. The proposed model predicts that the adaptability of threshold contrast is governed by unpredicted signal variations present in the visual signal, and thus represents an internal adjustment by the visual system that takes into account these unpredicted signal variations given the additional possibility for signal corruption by additive noise.     

Red blood cells under flow show maximal ATP release for specific hematocrit

ATP release by red blood cells (RBCs) under shear stress (SS) plays a pivotal role in endothelial biochemical signaling cascades. The aim of this study is to investigate through numerical simulation how RBC spatiotemporal organization depends on flow and geometrical conditions to generate ATP patterns. Numerical simulations were conducted in a straight channel by considering both plasma and explicit presence of RBCs, their shape deformation and cell-cell interaction, and ATP release by RBCs. Two ATP release pathways through cell membrane are taken into account: pannexin 1 channel, sensitive to SS, and cystic fibrosis transmembrane conductance regulator, which responds to cell deformation. Several flow and hematocrit conditions are explored. The problem is solved by the lattice Boltzmann method. Application of SS to the RBC suspension triggers a nontrivial spatial RBC organization and ATP patterns. ATP localizes preferentially in the vicinity of the cell-free layer close to channel wall. Conditions for maximal ATP release per cell are identified, which depend on vessel size and hematocrit Ht. Increasing further Ht beyond optimum enhances the total ATP release but should degrade oxygen transport capacity, a compromise between an efficient ATP release and minimal blood dissipation. Moreover, ATP is boosted in capillaries, suggesting a vasomotor activity coordination throughout the resistance network.