4π light scattering flow cytometry: enhancing the identification and characterization of individual cells.

The analysis of individual particles with complex morphologies from light scattering is crucial in disperse systems studies, such as blood cells. Characterization, which assumes determining particle characteristics, has a higher likelihood of succeeding in solving the inverse light-scattering problem if an instrument provides enough light-scattering data. In this study, we demonstrate how we extend the operating angular interval for the 4π Scanning Flow Cytometer (4πSFC), which measures angle-resolved light-scattering profiles (LSPs) of individual particles. The angular interval is extended by additionally measuring light scattering for the backward hemisphere. Currently, the 4πSFC setup uses three lasers, a single optical cell, and three photomultipliers. It enables the measurement of the LSP of individual particles within the angular interval of 10 to 170° for polar angles with integration over azimuth angles, which covers the spatial angle of 98.5% of the 4π angle. We demonstrate the 4πSFC's performance in measuring LSPs from the analysis of polymer beads, mature and spherized erythrocytes, and platelets. The 4πSFC has the potential to be very useful in identifying platelet dimers and granulocytes without labels, characterizing lymphocytes, monocytes, and abnormal erythrocytes.

Kinetic titration method in flow cytometry for quantification of cell receptors.

For the quantitative determination of cell receptors by fluorescence flow cytometry, we proposed a new method, which takes into account the reaction kinetics. The binding reaction of the ligand with receptors begins after placing the cells in the ligand solution. In the proposed method, there are several samples with the same concentration of cells and different initial concentrations of fluorescently labeled ligand, and each sample is measured by a flow cytometer once at the time when the following condition is met: the product of the incubation time (cells with ligand) and the initial concentration of ligand is the same for all samples. The proposed approach eliminates disadvantages and combines advantages of both kinetic and titration methods for quantification of receptors on single cells without the use of traditional calibration fluorescent beads. Practical application of the method was demonstrated in quantification of CD8 and CD14 on peripheral blood human leukocytes. Particularly, we found decreased (by a factor of two) mean number of CD14 on monocytes and granulocytes in patients with atherosclerosis (treated in the hospital) compared to conditionally healthy donors, whereas no difference was found in the mean CD8 expression on leukocytes between the same patient and donor groups.

A light scatter based model relating erythrocyte vesiculation to lifetime in circulation.

Methods for measuring erythrocyte age distribution are not available as a simple analytical tool. Most of them utilize the fluorescence or radioactive isotopes labeling to construct the age distribution and support physicians with aging indices of donor's erythrocytes. The age distribution of erythrocyte may be a useful snapshot of patient state over 120-days period of life. Previously, we introduced the enhanced assay of erythrocytes with measurement of 48 indices in four categories: concentration/content, morphology, aging and function (10.1002/cyto.a.24554). The aging category was formed by the indices based on the evaluation of the derived age of individual cells. The derived age does not exactly mean the real age of erythrocytes and its evaluation utilizes changes of cellular morphology during a lifespan. In this study, we are introducing the improved methodological approach that allows us to retrieve the derived age of individual erythrocytes, to construct the aging distribution, and to reform the aging category consisting of eight indices. The approach is based on the analysis of the erythrocyte vesiculation. The erythrocyte morphology is analyzed by scanning flow cytometry that measures the primary characteristics (diameter, thickness, and waist) of individual cells. The surface area (S) and sphericity index (SI) are calculated from the primary characteristics and the scattering diagram SI versus S is used in the evaluation of the derived age of each erythrocyte in a sample. We developed the algorithm to evaluate the derived age that provides eight indices in the aging category based on a model using light scatter features. The novel erythrocyte indices were measured for simulated cells and blood samples of 50 donors. We determined the first-ever reference intervals for these indices.

Proposed Dynamics of CDB3 Activation in Human Erythrocytes by Nifedipine Studied with Scanning Flow Cytometry.

Nifedipine is calcium channels and pumps blocker widely used in medicine. However, mechanisms of nifedipine action in blood are not clear. In particular, the influence of nifedipine on erythrocytes is far from completely understood. In this work, applying scanning flow cytometry, we observed experimentally for the first time the dynamics behind a significant increase of HCO3- /Cl- transmembrane exchange rate of CDB3 (main anion exchanger, AE1, Band 3, SLC4A1) of human erythrocytes in the presence of nifedipine in blood. It was found that the rate of CDB3 activation is not limited by the rate of nifedipine binding and/or Ca2+ transport. In order to explain the experimental data, we suggested a kinetic model assuming that the rate of CDB3 activation is limited by the dynamics of the balance between two intracellular processes (1) the activation of CDB3 limited by its interaction with intracellular Ca2+ , and (2) the spontaneous deactivation of CDB3. Thus the use of scanning flow cytometry allowed to clarify quantitatively the molecular kinetic mechanism of nifedipine action on human erythrocytes. In particular, the efficiency (~30) and rates of activation (~0.3 min-1 ) and deactivation (~10-3 min-1 ) of CDB3 in human erythrocytes was evaluated for two donors. © 2019 International Society for Advancement of Cytometry.

Calibration-free quantitative immunoassay by flow cytometry: Theoretical consideration and practical implementation for IgG antibody binding to CD14 receptors on human leukocytes.

We propose a calibration-free method to determine the number of receptors per cell, as well as the direct and the reverse reaction rate constants for a single receptor. The method is based on the analysis of the temporal evolution of the cells mean fluorescent intensity measured by a flow cytometer during the ligand-receptor (antigen-antibody) binding under the conditions of their comparable concentrations. We developed the kinetic approach accounting both for the delay between the dilution and the measurement and for the practical duration of the measurement itself. The method was applied to determine thenumber of CD14 receptors on human blood mononuclear (granulocytes, monocytes, lymphocytes) cells of several donors. We also obtained the direct ( k+= (5.6 ± 0.2) × 107 M-1 min-1 ) and reverse ( k-= (1.3 ± 0.2) × 10-2 min-1 ) rate constants of ligand-receptor interaction, and estimated the size of the binding site as b = 0.5 nm. The latter allows one to recalculate the rate constants for a different ligand, fluorescent label, medium viscosity, and/or temperature. The knowledge of the rate constants is essential for the calibration-free determination of the number of receptors per cell from a single kinetic curve of the cells mean fluorescence intensity.

Nuclear apoptotic volume decrease in individual cells: Confocal microscopy imaging and kinetic modeling.

The dynamics of nuclear morphology changes during apoptosis remains poorly investigated and understood. Using 3D time-lapse confocal microscopy we performed a study of early-stage apoptotic nuclear morphological changes induced by etoposide in single living HepG2 cells. These observations provide a definitive evidence that nuclear apoptotic volume decrease (AVD) is occurring simultaneously with peripheral chromatin condensation (so called "apoptotic ring"). In order to describe quantitatively the dynamics of nuclear morphological changes in the early stage of apoptosis we suggest a general molecular kinetic model, which fits well the obtained experimental data in our study. Results of this work may clarify molecular mechanisms of nuclear morphology changes during apoptosis.

Chylomicrons against light scattering: The battle for characterization.

Chylomicrons (CMs) are lipoprotein particles circulating in blood and transporting dietary lipids. Optically speaking, CMs are small compared to the wavelength of visible light and widely distributed by the size and refractive index (RI). Consequently, intensity of light scattered by the CMs scales with up to the sixth power of their size, hampering simultaneous analysis of 60 and 600 nm CMs. We present an accurate method for quantitative characterization of large-size CM subpopulation by the distributions over size and RI. For the first time the CM characteristics have been determined at a single particle level based on angle-resolved light-scattering measurements. We applied the developed method to 2 key processes relating to CM metabolism, namely in vivo dynamics of CMs in blood plasma after a meal and in vitro lipolysis of CMs by the lipoprotein lipase in postheparin plasma. We have observed the substantial variations in CM concentration, size and RI distributions. This opens the way for a multitude of medical applications involving screening of CM metabolism, which we exemplified by revealing large differences in CM characteristics after a 12-hour fast between a healthy volunteer and a patient with atherosclerosis.

Method for the simulation of blood platelet shape and its evolution during activation.

We present a simple physically based quantitative model of blood platelet shape and its evolution during agonist-induced activation. The model is based on the consideration of two major cytoskeletal elements: the marginal band of microtubules and the submembrane cortex. Mathematically, we consider the problem of minimization of surface area constrained to confine the marginal band and a certain cellular volume. For resting platelets, the marginal band appears as a peripheral ring, allowing for the analytical solution of the minimization problem. Upon activation, the marginal band coils out of plane and forms 3D convoluted structure. We show that its shape is well approximated by an overcurved circle, a mathematical concept of closed curve with constant excessive curvature. Possible mechanisms leading to such marginal band coiling are discussed, resulting in simple parametric expression for the marginal band shape during platelet activation. The excessive curvature of marginal band is a convenient state variable which tracks the progress of activation. The cell surface is determined using numerical optimization. The shapes are strictly mathematically defined by only three parameters and show good agreement with literature data. They can be utilized in simulation of platelets interaction with different physical fields, e.g. for the description of hydrodynamic and mechanical properties of platelets, leading to better understanding of platelets margination and adhesion and thrombus formation in blood flow. It would also facilitate precise characterization of platelets in clinical diagnosis, where a novel optical model is needed for the correct solution of inverse light-scattering problem.

Light-scattering gating and characterization of plasma microparticles.

Flow cytometry method (FCM) is widely used for analysis of cell-derived microparticles (MPs). Numerous efforts are currently aimed to standardize these measurements among different instruments. We push the FCM characterization of MPs to the limit based on rigorous simulation of measured signals. We measured forward- and side-scatter (FSC/SSC) signals and angle-resolved light-scattering profiles (LSPs) of polystyrene microspheres and MPs, including their aggregates, using a scanning flow cytometer (SFC). We used the Mie theory to (1) accurately evaluate instrument detection limits; (2) construct FSC/SSC gates for MPs in absolute scales of size and refractive index (RI); and (3) determine size and RI of individual spherical MPs. LSPs were used for advanced characterization, including differentiation of spherical and nonspherical particles. The proposed absolute FSC/SSC gating is naturally standardized for any FCM instrument, given the knowledge of its optical system and leads to instrument-independent analysis of MPs. The inverse Mie problem has a unique solution only for some regions of size and RI and uncertainties rapidly increase with decreasing size and RI. The developed methods are applicable to any flow cytometer, but are limited by assumption of particle sphericity. The latter can be relaxed only if additional signals, such as LSP, are measured.

Mature red blood cells: from optical model to inverse light-scattering problem.

We propose a method for characterization of mature red blood cells (RBCs) morphology, based on measurement of light-scattering patterns (LSPs) of individual RBCs with the scanning flow cytometer and on solution of the inverse light-scattering (ILS) problem for each LSP. We considered a RBC shape model, corresponding to the minimal bending energy of the membrane with isotropic elasticity, and constructed an analytical approximation, which allows rapid simulation of the shape, given the diameter and minimal and maximal thicknesses. The ILS problem was solved by the nearest-neighbor interpolation using a preliminary calculated database of 250,000 theoretical LSPs. For each RBC in blood sample we determined three abovementioned shape characteristics and refractive index, which also allows us to calculate volume, surface area, sphericity index, spontaneous curvature, hemoglobin concentration and content.

Influence of magnesium sulfate on HCO3/Cl transmembrane exchange rate in human erythrocytes.

Magnesium sulfate (MgSO4) is widely used in medicine but molecular mechanisms of its protection through influence on erythrocytes are not fully understood and are considerably controversial. Using scanning flow cytometry, in this work for the first time we observed experimentally (both in situ and in vitro) a significant increase of HCO3(-)/Cl(-) transmembrane exchange rate of human erythrocytes in the presence of MgSO4 in blood. For a quantitative analysis of the obtained experimental data, we introduced and verified a molecular kinetic model, which describes activation of major anion exchanger Band 3 (or AE1) by its complexation with free intracellular Mg(2+) (taking into account Mg(2+) membrane transport and intracellular buffering). Fitting the model to our in vitro experimental data, we observed a good correspondence between theoretical and experimental kinetic curves that allowed us to evaluate the model parameters and to estimate for the first time the association constant of Mg(2+) with Band 3 as KB~0.07mM, which is in agreement with known values of the apparent Mg(2+) dissociation constant (from 0.01 to 0.1mM) that reflects experiments on enrichment of Mg(2+) at the inner erythrocyte membrane (Gunther, 2007). Results of this work partly clarify the molecular mechanisms of MgSO4 action in human erythrocytes. The method developed allows one to estimate quantitatively a perspective of MgSO4 treatment for a patient. It should be particularly helpful in prenatal medicine for early detection of pathologies associated with the risk of fetal hypoxia.

Super-resolved calibration-free flow cytometric characterization of platelets and cell-derived microparticles in platelet-rich plasma.

Importance of microparticles (MPs), also regarded as extracellular vesicles, in many physiological processes and clinical conditions motivates one to use the most informative and precise methods for their characterization. Methods based on individual particle analysis provide statistically reliable distributions of MP population over characteristics. Although flow cytometry is one of the most powerful technologies of this type, the standard forward-versus-side-scattering plots of MPs and platelets (PLTs) overlap considerably because of similarity of their morphological characteristics. Moreover, ordinary flow cytometry is not capable of measurement of size and refractive index (RI) of MPs. In this study, we 1) employed the potential of the scanning flow cytometer (SFC) for identification and characterization of MPs from light scattering; 2) suggested the reference method to characterize MP morphology (size and RI) with high precision; and 3) determined the lowest size of a MP that can be characterized from light scattering with the SFC. We equipped the SFC with 405 and 488 nm lasers to measure the light-scattering profiles and side scattering from MPs, respectively. The developed two-stage method allowed accurate separation of PLTs and MPs in platelet-rich plasma. We used two optical models for MPs, a sphere and a bisphere, in the solution of the inverse light-scattering problem. This solution provides unprecedented precision in determination of size and RI of individual spherical MPs-median uncertainties (standard deviations) were 6 nm and 0.003, respectively. The developed method provides instrument-independent quantitative information on MPs, which can be used in studies of various factors affecting MP population.

Additivity of light-scattering patterns of aggregated biological particles.

The paper is focused on light scattering by aggregates of optically soft particles with a size larger than the wavelength, in particular, blood platelets. We conducted a systematic simulation of light scattering by dimers and larger aggregates of blood platelets, each modeled as oblate spheroids, using the discrete dipole approximation. Two-dimensional (2-D) light scattering patterns (LSPs) and internal fields showed that the multiple scattering between constituent particles can be neglected. Additionally, we derived conditions of the scattering angle and orientation of the dimer, under which the averaging of the 2-D LSPs over the azimuthal scattering angle washes out interference in the far field, resulting in averaged LSPs of the aggregate being equal to the sum of that for its constituents. We verified theoretical conclusions using the averaged LSPs of blood platelets measured with the scanning flow cytometer (SFC). Moreover, we obtained similar results for a model system of aggregates of polystyrene beads, studied both experimentally and theoretically. Finally, we discussed the potential of discriminating platelet aggregates from monomers using the SFC.

Brownian aggregation rate of colloid particles with several active sites.

We theoretically analyze the aggregation kinetics of colloid particles with several active sites. Such particles (so-called "patchy particles") are well known as chemically anisotropic reactants, but the corresponding rate constant of their aggregation has not yet been established in a convenient analytical form. Using kinematic approximation for the diffusion problem, we derived an analytical formula for the diffusion-controlled reaction rate constant between two colloid particles (or clusters) with several small active sites under the following assumptions: the relative translational motion is Brownian diffusion, and the isotropic stochastic reorientation of each particle is Markovian and arbitrarily correlated. This formula was shown to produce accurate results in comparison with more sophisticated approaches. Also, to account for the case of a low number of active sites per particle we used Monte Carlo stochastic algorithm based on Gillespie method. Simulations showed that such discrete model is required when this number is less than 10. Finally, we applied the developed approach to the simulation of immunoagglutination, assuming that the formed clusters have fractal structure.

Accurate measurement of volume and shape of resting and activated blood platelets from light scattering.

We introduce a novel approach for determination of volume and shape of individual blood platelets modeled as an oblate spheroid from angle-resolved light scattering with flow-cytometric technique. The light-scattering profiles (LSPs) of individual platelets were measured with the scanning flow cytometer and the platelet characteristics were determined from the solution of the inverse light-scattering problem using the precomputed database of theoretical LSPs. We revealed a phenomenon of parameter compensation, which is partly explained in the framework of anomalous diffraction approximation. To overcome this problem, additional a priori information on the platelet refractive index was used. It allowed us to determine the size of each platelet with subdiffraction precision and independent of the particular value of the platelet aspect ratio. The shape (spheroidal aspect ratio) distributions of platelets showed substantial differences between native and activated by 10 µM adenosine diphosphate samples. We expect that the new approach may find use in hematological analyzers for accurate measurement of platelet volume distribution and for determination of the platelet activation efficiency.

Arrangement of nuclear structures is not transmitted through mitosis but is identical in sister cells.

Although it is well known that chromosomes are non-randomly organized during interphase, it is not completely clear whether higher-order chromatin structure is transmitted from mother to daughter cells. Therefore, we addressed the question of how chromatin is rearranged during interphase and whether heterochromatin pattern is transmitted after mitosis. We additionally tested the similarity of chromatin arrangement in sister interphase nuclei. We noticed a very active cell rotation during interphase, especially when histone hyperacetylation was induced or transcription was inhibited. This natural phenomenon can influence the analysis of nuclear arrangement. Using photoconversion of Dendra2-tagged core histone H4 we showed that the distribution of chromatin in daughter interphase nuclei differed from that in mother cells. Similarly, the nuclear distribution of heterochromatin protein 1ß (HP1ß) was not completely identical in mother and daughter cells. However, identity between mother and daughter cells was in many cases evidenced by nucleolar composition. Moreover, morphology of nucleoli, HP1ß protein, Cajal bodies, chromosome territories, and gene transcripts were identical in sister cell nuclei. We conclude that the arrangement of interphase chromatin is not transmitted through mitosis, but the nuclear pattern is identical in naturally synchronized sister cells. It is also necessary to take into account the possibility that cell rotation and the degree of chromatin condensation during functionally specific cell cycle phases might influence our view of nuclear architecture.

Light-scattering flow cytometry for identification and characterization of blood microparticles.

We describe a novel approach to study blood microparticles using the scanning flow cytometer, which measures light scattering patterns (LSPs) of individual particles. Starting from platelet-rich plasma, we separated spherical microparticles from non-spherical plasma constituents, such as platelets and cell debris, based on similarity of their LSP to that of sphere. This provides a label-free method for identification (detection) of microparticles, including those larger than 1 µm. Next, we rigorously characterized each measured particle, determining its size and refractive index including errors of these estimates. Finally, we employed a deconvolution algorithm to determine size and refractive index distributions of the whole population of microparticles, accounting for largely different reliability of individual measurements. Developed methods were tested on a blood sample of a healthy donor, resulting in good agreement with literature data. The only limitation of this approach is size detection limit, which is currently about 0.5 µm due to used laser wavelength of 0.66 µm.

Distribution function approach to the study of the kinetics of IgM antibody binding to FcγRIIIb (CD16b) receptors on neutrophils by flow cytometry.

Though flow cytometry provides the entire distribution of cellular fluorescence (i.e., "fluorescence profile"), only mean fluorescence data are usually considered in studies of ligand-receptor binding. In this study, we presented a method of the treatment of the temporal evolution of the whole fluorescence profile with a comprehensive statistical approach extended to the reversible binding case. The method was demonstrated in the study of the 1D3 IgM monoclonal antibodies binding to FcγRIIIb receptors (CD16b) on neutrophils. Kinetic experiments were carried out using a FACSCalibur (Becton Dickinson, USA) flow cytometer. For each of four donors, we obtained the distribution of the number of FcγRIIIb surface receptors for neutrophils and the rate constants per receptor: the association rate constant of (2.7±0.4)×10(7) M(-1) min(-1), and the dissociation rate constant of (1.3±0.4)×10(-1) min(-1). Based on the obtained values, the size of the receptor reaction site was estimated at approximately 1 nm. It was found, that cell receptors distributions differed sufficiently between donors in mean and the skewness values, whereas the coefficient of variation (i.e., the ratio of the standard deviation to the mean) did not vary significantly.

A nonfitting method using a spatial sine window transform for inhomogeneous effective-diffusion measurements by FRAP.

Determining averaged effective diffusion constants from experimental measurements of fluorescent proteins in an inhomogeneous medium in the presence of ligand-receptor interactions poses problems of analytical tractability. Here, we introduced a nonfitting method to evaluate the averaged effective diffusion coefficient of a region of interest (which may include a whole nucleus) by mathematical processing of the entire cellular two-dimensional spatial pattern of recovered fluorescence. Spatially and temporally resolved measurements of protein transport inside cells were obtained using the fluorescence recovery after photobleaching technique. Two-dimensional images of fluorescence patterns were collected by laser-scanning confocal microscopy. The method was demonstrated by applying it to an estimation of the mobility of green fluorescent protein-tagged heterochromatin protein 1 in the nuclei of living mouse embryonic fibroblasts. This approach does not require the mathematical solution of a corresponding system of diffusion-reaction equations that is typical of conventional fluorescence recovery after photobleaching data processing, and is most useful for investigating highly inhomogeneous areas, such as cell nuclei, which contain many protein foci and chromatin domains.

In situ optical monitoring of RDX nanoparticles formation during rapid expansion of supercritical CO2 solutions.

Nanoparticles of RDX (cyclotrimethylenetrinitramine) generated by RESS (rapid expansion of supercritical solutions) using supercritical CO2 were characterized in situ by a pulsed laser light scattering imaging technique using a gated ICCD (intensified CCD) camera. The absolute sensitivity calibration was performed using Rayleigh light scattering from air as well as light scattering from standard polystyrene spheres. The size distribution functions of the particles formed in the RESS jet were determined using the calibrated sensitivity. The diameter of RDX particles formed at the pre-expansion pressure of 180 bar was 73 nm at the maximum of the size distribution function. Assuming that the particles near the nozzle consisted mainly of CO2 and the size distribution was log-normal, the diameter of the particles near the nozzle (7.5 mm from the nozzle) at the distribution maximum was 3.3 microm at the pre-expansion pressure of 180 bar. The number densities of the particles in the RESS jet were determined by counting individual particles in the light scattering images. Based on the measured particle size distributions and the number density of particles along the RESS jet, the mechanism of particle formation in RESS is discussed. The homogeneous nucleation mechanism is rejected as it fails to explain the large particle size experimentally observed. Instead, a modified "spray-drying" mechanism is suggested.