Proliferation of cells with aggregation and communication.

Cell proliferation is often considered to occur via front propagation with constant velocity. This scenario proposed by Fisher, Kolmogorov, Petrovsky, and Piskunov is based on the solution of the corresponding mean-field reaction-diffusion equations and does not take into account that due to adhesion the cells have tendency to aggregate and that the rate of cell division may depend on the cell-cell communication. Herein, the author presents extensive Monte Carlo simulations taking both these factors into account and illustrating that the former factor can dramatically modify the spatio-temporal kinetics of cell proliferation. In particular, the conventional relation between the front velocity and diffusion coefficient may fail, the front velocity may appreciably increase with increasing time, and/or the front may be partly or fully smeared on the realistic length scales.

Modeling tensional homeostasis in multicellular clusters.

Homeostasis of mechanical stress in cells, or tensional homeostasis, is essential for normal physiological function of tissues and organs and is protective against disease progression, including atherosclerosis and cancer. Recent experimental studies have shown that isolated cells are not capable of maintaining tensional homeostasis, whereas multicellular clusters are, with stability increasing with the size of the clusters. Here, we proposed simple mathematical models to interpret experimental results and to obtain insight into factors that determine homeostasis. Multicellular clusters were modeled as one-dimensional arrays of linearly elastic blocks that were either jointed or disjointed. Fluctuating forces that mimicked experimentally measured cell-substrate tractions were obtained from Monte Carlo simulations. These forces were applied to the cluster models, and the corresponding stress field in the cluster was calculated by solving the equilibrium equation. It was found that temporal fluctuations of the cluster stress field became attenuated with increasing cluster size, indicating that the cluster approached tensional homeostasis. These results were consistent with previously reported experimental data. Furthermore, the models revealed that key determinants of tensional homeostasis in multicellular clusters included the cluster size, the distribution of traction forces, and mechanical coupling between adjacent cells. Based on these findings, we concluded that tensional homeostasis was a multicellular phenomenon. Copyright © 2016 John Wiley & Sons, Ltd.

Kinetic Monte Carlo and cellular particle dynamics simulations of multicellular systems.

Computer modeling of multicellular systems has been a valuable tool for interpreting and guiding in vitro experiments relevant to embryonic morphogenesis, tumor growth, angiogenesis and, lately, structure formation following the printing of cell aggregates as bioink particles. Here we formulate two computer simulation methods: (1) a kinetic Monte Carlo (KMC) and (2) a cellular particle dynamics (CPD) method, which are capable of describing and predicting the shape evolution in time of three-dimensional multicellular systems during their biomechanical relaxation. Our work is motivated by the need of developing quantitative methods for optimizing postprinting structure formation in bioprinting-assisted tissue engineering. The KMC and CPD model parameters are determined and calibrated by using an original computational-theoretical-experimental framework applied to the fusion of two spherical cell aggregates. The two methods are used to predict the (1) formation of a toroidal structure through fusion of spherical aggregates and (2) cell sorting within an aggregate formed by two types of cells with different adhesivities.

Morphological assessment of basic multicellular unit resorption parameters in dogs shows additional mechanisms of bisphosphonate effects on bone.

Bisphosphonates (BPs) slow bone loss by reducing initiation of new basic multicellular units (BMUs). Whether or not BPs simply prevent osteoclasts from initiating new BMUs that resorb bone or also reduce the amount of bone they resorb at the BMU level is not clear. The goal of this study was to determine the effects of BPs on three morphological parameters of individual BMUs, resorption depth (Rs.De), area (Rs.Ar), and width (Rs.Wi). After 1 year of treatment with vehicle (VEH), alendronate (ALN; 0.10, 0.20, or 1.00 mg/kg/day), or risedronate (RIS; 0.05, 0.10, or 0.50 mg/kg/day), resorption cavity morphology was assessed in vertebral trabecular bone of beagle dogs by histology. Animals treated with ALN or RIS at the doses representing those used to treat postmenopausal osteoporosis (0.20 and 0.10 mg/kg/day, respectively) had significantly lower Rs.Ar (-27%) and Rs.Wi (-17%), with no difference in Rs.De, compared to VEH-treated controls. Low doses of ALN and RIS did not affect any parameters, whereas higher doses resulted in similar changes to those of the clinical dose. There were no significant differences in the resorption cavity measures between RIS and ALN at any of the dose equivalents. These results highlight the importance of examining parameters beyond erosion depth for assessment of resorption parameters. Furthermore, these results suggest that in addition to the well-known effects of BPs on reducing the number of active BMUs, these drugs also reduce the activity of osteoclasts at the individual BMU level at doses at and above those used clinically for the treatment of postmenopausal osteoporosis.

Development of a technique for quantification of reticulocytes and assessment of erythrocyte regenerative capacity in birds.

OBJECTIVE: To develop a reticulocyte classification scheme, optimize an avian reticulocyte staining protocol, and compare the percentages of reticulocyte types with polychromatophil percentage in blood samples from birds. SAMPLE POPULATION: Blood samples from a red-tailed hawk and 31 ill birds. PROCEDURES: A single blood sample obtained from a red-tailed hawk (Buteo jamaicensis) was used to optimize the staining protocol. For optimization of the staining protocol, 4 dilutions of whole blood with new methylene blue stain and 4 incubation times were evaluated. From samples submitted for avian CBCs, EDTA-anticoagulated whole blood samples from 31 ill birds were randomly selected and examined to compare polychromatophil and reticulocyte percentages. Reticulocyte staining was performed in all samples by use of a 1:3 (whole blood to new methylene blue) dilution with incubation for 10 minutes at room temperature (approx 22 degrees C); reticulocytes were assessed as a percentage of 1,000 RBCs by 2 independent observers. In Wright-Giemsa-stained blood smears, a polychromatophil percentage was similarly determined. RESULTS: 4 avian reticulocyte types were defined: ring-form reticulocytes, aggregate reticulocytes, and 2 subcategories of punctate reticulocytes. A reticulocyte-staining protocol was optimized. Interobserver and intraobserver variations in assessment of reticulocyte and polychromatophil percentages were not significant. A strong positive correlation (Spearman coefficient of rank correlation [rho] = 0.978) was identified between the percentage of polychromatophils and the percentage of ring-form reticulocytes. CONCLUSIONS AND CLINICAL RELEVANCE: Results indicated that quantification of ring-form reticulocytes provides an accurate assessment of erythrocyte regenerative capacity in birds.

Assessment of pluripotency and multilineage differentiation potential of NTERA-2 cells as a model for studying human embryonic stem cells.

Embryonal carcinoma cells are pluripotent stem cells derived from teratocarcinomas and are considered to be the malignant counterparts of human embryonic stem cells. As there are few reliable experimental systems available to study the molecular mechanisms governing normal embryogenesis, well-characterized human embryonal carcinoma stem cell lines may provide a robust and simple model to study certain aspects of pluripotency and cellular differentiation. Here, we have analysed NTERA-2 cL.D1 cells at molecular and cellular levels during expansion and differentiation, via formation of cell aggregates similar to embryoid bodies in embryonic stem cells. Thus, human embryonal carcinoma cells may provide a valuable insight into cell fate determination, into the embryonic ectoderm, mesoderm and endoderm and their downstream derivatives.

Effect of red cell clustering and anisotropy on ultrasound blood backscatter: a Monte Carlo study.

When flowing at a low shear rate, blood appears hyperechogenic on ultrasound B-scans. The formation of red blood cell (RBC) aggregates that also alters blood viscosity is the microscopic mechanism explaining this acoustical phenomenon. In this study, Monte Carlo simulations were performed to predict how RBC clustering increases ultrasound scattering by blood. A bidimensional Gibbs-Markov random point process parameterized by the adhesion energy epsilon and an anisotropy index nu was used to describe RBC positions for a hematocrit H = 40%. The frequency dependence of the backscattering coefficient chi(f) was computed using Born approximation. The backscattering coefficient chi0 at 5 MHz and the spectral slopes n(x) and n(y) (chi alpha f(nx) or f(ny)) measured, respectively, when the insonification is parallel and perpendicular with the RBC cluster axis were then extracted. Under isotropic conditions, chi0 increased up to 7 dB with epsilon and n(x) = n(y) decreased from 4.2 to 3.4. Under anisotropic conditions, the backscattering was stronger perpendicularly to aggregate axis, resulting in n(x) < n(y). The anisotropy in scattering appeared more pronounced when epsilon or nu increased. These two dimensional results generally predict that low-frequency blood backscatter is related to cluster dimension, and higher-frequency properties are affected by finer morphological features as anisotropy. This numerically establishes that ultrasound backscatter spectroscopy on a large frequency range is pertinent to characterize in situ hemorheology.

The influence of psychological stress and socioeconomic status on platelet activation in men.

OBJECTIVE: Circulating monocyte- and neutrophil-platelet aggregates are sensitive markers of in vivo platelet activation. Socioeconomic status is inversely associated with risk of coronary heart disease. We assessed the impact of psychological stress on leukocyte-platelet aggregates in men from higher and lower socioeconomic status groups. METHODS: Participants were 37 healthy non-smoking men aged 30-59 years, divided by occupation into higher and lower social status groups. Blood was drawn at baseline, immediately following stressful behavioural tasks, and at 30 and 75 min post-stress, and aggregates were analysed using flow cytometry. Cardiovascular and subjective stress responses were also monitored. RESULTS: There were significant increases following stress in monocyte-, neutrophil-, lymphocyte- and total leukocyte-platelet aggregates (all P<0.05). The largest responses were in monocyte-platelet (21% increase) and neutrophil-platelet (16.7% increase) aggregates. Lower socioeconomic status men had greater numbers of leukocyte-platelet aggregates throughout, but the magnitude of stress responses did not vary with social status. The increase in monocyte- and leukocyte-platelet aggregates was associated with systolic blood pressure stress responsivity. CONCLUSIONS: Psychological stress induces platelet activation as indexed by leukocyte-platelet aggregates, and correlations with cardiovascular stress reactions suggest that sympathoadrenal responses may be responsible. Platelet activation may be a mechanism through which social position influences cardiovascular disease risk.

Monte Carlo simulation of the heterotypic aggregation kinetics of platelets and neutrophils.

The heterotypic aggregation of cell mixtures or colloidal particles such as proteins occurs in a variety of settings such as thrombosis, immunology, cell separations, and diagnostics. Using the set of population balance equations (PBEs) to predict dynamic aggregate size and composition distributions is not feasible. The stochastic algorithm of Gillespie for chemical reactions (. J. Comput. Phys. 22:403-434) was reformulated to simulate the kinetic behavior of aggregating systems. The resulting Monte Carlo (MC) algorithm permits exact calculation of the decay rates of monomers and the temporally evolving distribution of sizes and compositions of the aggregates. Moreover, it permits calculation of all moments of these distributions. Using this method, we explored the heterotypic aggregation of fully activated platelets and neutrophils in a linear shear flow of shear rate G = 335 s(-1). At plasma concentrations, the half-lives of homotypically aggregating platelet and neutrophil singlets were 8.5 and 2.4 s, respectively. However, for heterotypic aggregation, the half-lives for platelets and neutrophils decreased to 2.0 and 0.11 s, respectively, demonstrating that flowing neutrophils accelerate capture of platelets and growth of aggregates. The required number of calculations per time step of the MC algorithm was typically a small fraction of Omega(1/2), where Omega is the initial number of particles in the system, making this the fastest MC method available. The speed of the algorithm makes feasible the deconvolution of kernels for general biological heterotypic aggregation processes.

Assessment of neutrophil aggregation by Coulter STKR and STKS haematological analysers.

We have studied an alternative method to aggregometry for the assessment of human polymorphonuclear (PMN) leucocyte aggregation. This simple, rapid and reliable procedure counts unaggregated cells on both Coulter STKS and STKR haematological analysers by the impedance principle. Aggregation of PMN was induced by 15 min incubation with fresh autologous serum (FAS) after a 10 min phorbol myristate acetate (PMA) activation of neutrophils in small aliquots (0.25 ml) of suspension containing about 4.0 x 10(9) PMN/1. Differences (x 100) between count of resting and PMA+FAS treated neutrophils/count of resting PMN reflect percent aggregation. By this procedure, PMN aggregation did not occur in autologous plasma from EDTA anticoagulated whole blood; it was partially inhibited by hydrocortisone, whereas inactivated or Zymosan activated sera gave values similar to those from FAS induced aggregation. PMA aggregation was dependent on Ca2+ + Mg2+ concentration. Intra-assay analytical variability did not exceed 4% on either instrument. Reference values (n = 20) of percent PMN aggregation were 50.7 +/- 4.7 on STKS and 47.1 +/- 4.8 on STKR. Most probably, the interindividual variance was due to the physiological variability of Mg2+ and/or Ca2+ concentrations in FAS. Thus, this procedure reflects the true PMN aggregability status in a given subject, and in a given electrolyte environment.