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

In vitro Research Concerning Effect of Clopidogrel Alone and on Combination with Aspirin and Dypiridamoleon Sedimentation of Erythrocytes

Based on extended theory of Derjaguin, Landau and Overbeeck (xDLVO) concerning aggregation of colloids and biological cells it was hypothesized that platelet antiaggregant agents have to reduce the aggregation of erythrocytes also. Applying Einstein-Stokes theory of sedimentation of spheres in viscous media it was concluded that sedimentation of erythrocytes is in fact sedimentation of aggregates of the approximately same size. Consequently, an expected outcome was that addition of antiaggregants in vitro to blood samples from patients with rheumatic or cardiovascular diseases will be the decrease of erythrocytes sedimentation. Starting from usual practice of dual antiaggregant therapy (aspirin and clopidogrel) effects of clopidogrel were compared with effects of clopidogrel plus small concentrations of aspirine and dipyridamole (smaller that their concentrations in plasma after in vivo administration) in order to put in evidence a possible synergic effect at platelet membrane level. Whole blood (0.8ml) was collected on 0.1ml 1% EDTA and then was added 25 or 50µl normal saline solution of clopidogrel or of the combination acetylosalycilic acid, clopidogrel and dypyridamole. The final concentrations were 1, 2 and µg/ml, of the same order as cumulated concentration of clopidogrel and its metabolites in clinical pharmacokinetics. Experiments were performed on a number of 40 human blood samples obtained from 2 groups of 20 patients. Sedimentation of erythrocytes was recorded using a camera and captured data were stored on a computer. Global analysis evidenced that in presence of antiaggregants the clusters of the sedimentation curves shifted down and into right, indicating a decrease and delay of sedimentation. Initial slopes and extent of sedimentation decreased linearly on clopidogrel concentration within the 1-3µg/ml range. For comparison of mean curves corresponding to different clopidogrel concentrations it was applied a metric from biopharmacy: areas under plasma concentrations curves (AUC) of drugs. The areas under average sedimentation curves decreased linearly at clopidogrel concentration within the 1-3µg/ml range. The same experiments were performed and similar results were obtained with the triple antiaggregant combination (clopidogrel, acetylosalycilic acid and dipyridamole). Apparently, a synergism between the tested antiaggregants appeared at studied concentration but the number of data was not sufficient to prove the statistical significance of the difference between clopidogrel alone and in triple combination.

Multi-scale biological and physical modelling of the tumour micro-environment

Paced by advances in high performance computing, and algorithms for multi-physics and multi-scale simulation, a number of groups have recently established numerical models of flowing blood systems, where cell-scale interactions are explicitly resolved. To be biologically representative, these models account for some or all of: (1) fluid dynamics of the carrier flow, (2) structural dynamics of the cells and vessel walls, (3) interaction and transport biochemistry, and, (4) methods for scaling to physiologically representative numbers of cells. In this article, our interest is the modelling of the tumour micro-environment. We review the broader area of cell-scale resolving blood flow modelling, while focusing on the particular interactions of tumour cells and white blood cells, known to play an important role in metastasis.

Systems biology of coagulation

Accurate computer simulation of blood function can inform drug target selection, patient-specific dosing, clinical trial design, biomedical device design, as well as the scoring of patient-specific disease risk and severity. These large-scale simulations rely on hundreds of independently measured physical parameters and kinetic rate constants. However, the models can be validated against large-scale, patient-specific laboratory measurements. By validation with high-dimensional data, modeling becomes a powerful tool to predict clinically complex scenarios. Currently, it is possible to accurately predict the clotting rate of plasma or blood in a tube as it is activated with a dose of tissue factor, even as numerous coagulation factors are altered by exogenous attenuation or potentiation. Similarly, the dynamics of platelet activation, as indicated by calcium mobilization or inside-out signaling, can now be numerically simulated with accuracy in cases where platelets are exposed to combinations of agonists. Multiscale models have emerged to combine platelet function and coagulation kinetics into complete physics-based descriptions of thrombosis under flow. Blood flow controls platelet fluxes, delivery and removal of coagulation factors, adhesive bonding, and von Willebrand factor conformation. The field of blood systems biology has now reached a stage that anticipates the inclusion of contact, complement, and fibrinolytic pathways along with models of neutrophil and endothelial activation. Along with '-omics' data sets, such advanced models seek to predict the multifactorial range of healthy responses and diverse bleeding and clotting scenarios, ultimately to understand and improve patient outcomes.

Microfluidics and coagulation biology

The study of blood ex vivo can occur in closed or open systems, with or without flow. Microfluidic devices, which constrain fluids to a small (typically submillimeter) scale, facilitate analysis of platelet function, coagulation biology, cellular biorheology, adhesion dynamics, and pharmacology and, as a result, can be an invaluable tool for clinical diagnostics. An experimental session can accommodate hundreds to thousands of unique clotting, or thrombotic, events. Using microfluidics, thrombotic events can be studied on defined surfaces of biopolymers, matrix proteins, and tissue factor, under constant flow rate or constant pressure drop conditions. Distinct shear rates can be generated on a device using a single perfusion pump. Microfluidics facilitated both the determination of intraluminal thrombus permeability and the discovery that platelet contractility can be activated by a sudden decrease in flow. Microfluidic devices are ideal for multicolor imaging of platelets, fibrin, and phosphatidylserine and provide a human blood analog to mouse injury models. Overall, microfluidic advances offer many opportunities for research, drug testing under relevant hemodynamic conditions, and clinical diagnostics.

Computational approaches to studying thrombus development

In addition to descriptive biological models, many computational models have been developed for hemostasis/thrombosis that provide quantitative characterization of thrombus development. Simulations using computational models that have been developed for coagulation reactions, platelet activation, and fibrinogen assembly have been shown to be in close agreement with experimental data. Models of processes involved in hemostasis/thrombosis are being integrated to simulate the development of the thrombus simultaneously in time and space. Further development of computational approaches can provide quantitative insights leading to predictions that are not obvious from qualitative biological models.

Scalable stirred-suspension bioreactor culture of human pluripotent stem cells

Advances in stem cell biology have afforded promising results for the generation of various cell types for therapies against devastating diseases. However, a prerequisite for realizing the therapeutic potential of stem cells is the development of bioprocesses for the production of stem cell progeny in quantities that satisfy clinical demands. Recent reports on the expansion and directed differentiation of human embryonic stem cells (hESCs) in scalable stirred-suspension bioreactors (SSBs) demonstrated that large-scale production of therapeutically useful hESC progeny is feasible with current state-of-the-art culture technologies. Stem cells have been cultured in SSBs as aggregates, in microcarrier suspension and after encapsulation. The various modes in which SSBs can be employed for the cultivation of hESCs and human induced pluripotent stem cells (hiPSCs) are described. To that end, this is the first account of hiPSC cultivation in a microcarrier stirred-suspension system. Given that cultured stem cells and their differentiated progeny are the actual products used in tissue engineering and cell therapies, the impact of bioreactor's operating conditions on stem cell self-renewal and commitment should be considered. The effects of variables specific to SSB operation on stem cell physiology are discussed. Finally, major challenges are presented which remain to be addressed before the mainstream use of SSBs for the large-scale culture of hESCs and hiPSCs.

Application of Population Dynamics to Study Heterotypic Cell Aggregations in the Near-Wall Region of a Shear Flow

Our research focused on the polymorphonuclear neutrophils (PMNs) tethering to the vascular endothelial cells (EC) and the subsequent melanoma cell emboli formation in a shear flow, an important process of tumor cell extravasation from the circulation during metastasis. We applied population balance model based on Smoluchowski coagulation equation to study the heterotypic aggregation between PMNs and melanoma cells in the near-wall region of an in vitro parallel-plate flow chamber, which simulates in vivo cell-substrate adhesion from the vasculatures by combining mathematical modeling and numerical simulations with experimental observations. To the best of our knowledge, a multiscale near-wall aggregation model was developed, for the first time, which incorporated the effects of both cell deformation and general ratios of heterotypic cells on the cell aggregation process. Quantitative agreement was found between numerical predictions and in vitro experiments. The effects of factors, including: intrinsic binding molecule properties, near-wall heterotypic cell concentrations, and cell deformations on the coagulation process, are discussed. Several parameter identification approaches are proposed and validated which, in turn, demonstrate the importance of the reaction coefficient and the critical bond number on the aggregation process.

Systems biology to predict blood function

Systems biology seeks to provide a quantitative framework to understand blood as a reactive biological fluid whose function is dictated by prevailing haemodynamics, vessel wall characteristics, platelet metabolism, numerous coagulation factors in plasma, and small molecules released during thrombosis. The hierarchical nature of thrombosis requires analysis of adhesive bond dynamics of activated platelets captured from a flow field to a growing thrombus boundary along with the simultaneous assembly of the coagulation pathway. Several kinetic models of protease cascades have been developed. A full bottom-up model of platelet intracellular metabolism is now available to simulate the metabolism of resting platelets and platelets exposed to activators. Monte Carlo algorithms can finally accommodate platelet reaction, dispersion, and convection for full simulation of platelet deposition and clotting under flow. For clinical applications, the systems biology prediction of patient-specific pharmacological response requires the final assembly of platelet intracellular metabolism models with coagulation protease network models.

Similarities between heterophilic and homophilic cadherin adhesion

The mechanism that drives the segregation of cells into tissue-specific subpopulations during development is largely attributed to differences in intercellular adhesion. This process requires the cadherin family of calcium-dependent glycoproteins. A widely held view is that protein-level discrimination between different cadherins on cell surfaces drives this sorting process. Despite this postulated molecular selectivity, adhesion selectivity has not been quantitatively verified at the protein level. In this work, molecular force measurements and bead aggregation assays tested whether differences in cadherin bond strengths could account for cell sorting in vivo and in vitro. Studies were conducted with chicken N-cadherin, canine E-cadherin, and Xenopus C-cadherin. Both qualitative bead aggregation and quantitative force measurements show that the cadherins cross-react. Furthermore, heterophilic adhesion is not substantially weaker than homophilic adhesion, and the measured differences in adhesion do not correlate with cell sorting behavior. These results suggest that the basis for cell segregation during morphogenesis does not map exclusively to protein-level differences in cadherin adhesion.

Mechanism and dynamics of cadherin adhesion

Cadherins are essential cell adhesion molecules involved in tissue morphogenesis and the maintenance of tissue architecture in adults. The adhesion and selectivity functions of cadherins are located in their extracellular regions. Biophysical studies show that the adhesive activity is not confined to a single interface. Instead, multiple cadherin domains contribute to binding. By contrast, the specificity-determining site maps to the N-terminal domains, which adhere by the reciprocal binding of Trp2 residues from opposing proteins. Structural cooperativity can transmit the effects of subtle structural changes or ligand binding over large distances in the protein. Increasingly, studies show that differential cadherin-mediated adhesion, rather than exclusive homophilic binding between identical cadherins, direct cell segregation and the organization of tissue interfaces during morphogenesis. Force measurements quantified both kinetic and strength differences between different classical cadherins that may underlie cell sorting behavior. Despite the complex adhesion mechanisms and differences in binding properties, cadherin-mediated cell adhesion is also regulated by many other biochemical processes. Elucidating the mechanisms by which cadherins organize cell junctions and tissue architecture requires not only quantitative, mechanistic investigations of cadherin function but also investigations of the biochemical and cellular processes that can modulate those functions.

Stochastic Modeling of Blood Coagulation Initiation

A kinetic Monte Carlo simulation was developed using the deterministic reaction network developed by the Mann laboratory for tissue-factor (TF)-initiated blood coagulation. The model predicted thrombin dynamics in recalcified whole blood (3-fold diluted) pretreated with convulxin (platelet GPVI activator) and picomolar levels of TF (0-14 pM). The model did not accurately predict coagulation times at low TF (0-0.7 pM). The simulation revealed that approximately 0.2 pM TF was the critical concentration to cause 50% of reactions containing 3-fold diluted whole blood to reach a clotting threshold of 0.05 U/ml thrombin by 1 h. Simulations of 1 nl of blood (5 pM TF) revealed small stochastic variations in thrombin initiation time, while 16.6 pl simulations were highly stochastic at this level of TF (50 molecules/16.6 pl). Further experiment and simulation will require evaluation of mechanisms of coagulation kinetics at subpicomolar levels of TF.

Preferential binding of platelets to monocytes over neutrophils under flow

This study was undertaken to systematically investigate the binding kinetics of platelet recruitment by monocytes relative to neutrophils in bulk suspensions subjected to shear as well as the molecular requirements of leukocyte-platelet binding. Hydrodynamic shear-induced collisions augment the proportion of monocytes with adherent platelets more drastically than that of neutrophils with bound platelets. These heterotypic interactions are further potentiated by platelet activation with thrombin or to a lesser extent by monocyte stimulation with N-formyl-methionyl-leucyl-phenylalanine (fMLP). Monocyte-platelet heteroaggregation increases with increasing shear rate and shear exposure time. Platelet P-selectin binding to monocyte P-selectin-glycoprotein-ligand-1 is solely responsible for maximal platelet adhesion to unstimulated monocytes in shear flow. However, the enhanced platelet binding to fMLP-treated monocytes involves a sequential two-step process, wherein P-selectin-PSGL-1 interactions are stabilized by CD18-integrin involvement. Blocking platelet alpha(IIb)beta(3) or monocyte beta(1)-integrin function had no effect. This study underscores the preferential recruitment of platelets by monocytes relative to neutrophils in shear flow, and demonstrates that the shear environment of the vasculature coupled to the state of cell activation modulates the dynamics and molecular constituents mediating monocyte-platelet adhesion.

Fluid Shear Regulates the Kinetics and Receptor Specificity ofStaphylococcus aureusBinding to Activated Platelets

The interaction between surface components on the invading pathogen and host cells such as platelets plays a key role in the regulation of endovascular infections. However, the mechanisms mediating Staphylococcus aureus binding to platelets under shear remain largely unknown. This study was designed to investigate the kinetics and molecular requirements of platelet-S. aureus interactions in bulk suspensions subjected to a uniform shear field. Hydrodynamic shear-induced collisions augment platelet-S. aureus binding, which is further potentiated by platelet activation with stromal derived factor-1beta. Peak adhesion efficiency occurs at low shear (100 s(-1)) and decreases with increasing shear. The molecular interaction of platelet alpha(IIb)beta(3) with bacterial clumping factor A through fibrinogen bridging is necessary for stable bacterial binding to activated platelets under shear. Although this pathway is sufficient at low shear (</=400 s(-1)), the involvement of platelet gpIb and staphylococcal protein A through von Willebrand factor bridging is essential for optimal recruitment of S. aureus cells by platelets in the high shear regime. IgG plays an inhibitory role in the adhesion process, presumably by interfering with the binding of von Willebrand factor to staphylococcal protein A. This study demonstrates that platelet activation and a fluid-mechanical environment representative of the vasculature affect platelet-S. aureus cell-adhesive interactions pertinent to the process of S. aureus-induced bloodstream infections.

Kinetics of random aggregation-fragmentation processes with multiple components

A computationally efficient algorithm is presented for exact simulation of the stochastic time evolution of spatially homogeneous aggregation-fragmentation processes featuring multiple components or conservation laws. The algorithm can predict the average size and composition distributions of aggregating particles as well as their fluctuations, regardless of the functional form (e.g., composition dependence) of the aggregation or fragmentation kernels. Furthermore, it accurately predicts the complete time evolutions of all moments of the size and composition distributions, even for systems that exhibit gel transitions. We demonstrate the robustness and utility of the algorithm in case studies of linear and branched polymerization processes, the last of which is a two-component process. These simulation results provide the stochastic description of these processes and give new insights into their gel transitions, fluctuations, and long-time behavior when deterministic approaches to aggregation kinetics may not be reliable.

Fluid shear- and time-dependent modulation of molecular interactions between PMNs and colon carcinomas

This study compares the effects of fluid shear on the kinetics, adhesion efficiency, stability, and molecular requirements of polymorphonuclear leukocyte (PMN) binding to two colon adenocarcinoma cell-lines, the CD54-negative/sLe(x)-bearing LS174T cells and the CD54-expressing/sLe(x)-low HCT-8 cells. The efficiency of PMN-colon carcinoma heteroaggregation decreases with increasing shear, with PMNs binding HCT-8 more efficiently than LS174T cells at low shear (50-200 s(-1)). In the low shear regime, CD11b is sufficient to mediate PMN binding to LS174T cells. In contrast, both CD11a and CD11b contribute to PMN-HCT-8 heteroaggregation, with CD54 on HCT-8 cells acting as a CD11a ligand at early time points. At high shear, only PMN-LS174T heteroaggregation occurs, which is initiated by PMN L-selectin binding to a sialylated, O-linked, protease-sensitive ligand on LS174T cells. PMN-LS174T heteroaggregation is primarily dependent on the intercellular contact duration (or shear rate), whereas PMN-HCT-8 binding is a function of both the intercellular contact duration and the applied force (or shear stress). Cumulatively, these studies suggest that fluid shear modulates the kinetics and molecular mechanisms of PMN-colon carcinoma cell aggregation.

Fluid shear regulates the kinetics and molecular mechanisms of activation-dependent platelet binding to colon carcinoma cells

This study was undertaken to investigate the kinetics and molecular requirements of platelet binding to tumor cells in bulk suspensions subjected to a uniform linear shear field, using a human colon adenocarcinoma cell line (LS174T) as a model. The effects of shear rate (20-1000 s(-1)), shear exposure time (30-300 s), shear stress (at constant shear rate by adjusting the viscosity of the medium from 1.3-2.6 cP), cell concentration, and platelet activation on platelet-LS174T heteroaggregation were assessed. The results indicate that hydrodynamic shear-induced collisions augment platelet-LS174T binding, which is further potentiated by thrombin/GPRP-NH(2). Peak adhesion efficiency occurs at low shear and decreases with increasing shear. Intercellular contact duration is the predominant factor limiting heteroaggregation at shear rates up to 200 s(-1), whereas these interactions become shear stress-sensitive at > or = 400 s(-1). Heteroaggregation increases with platelet concentration due to an elevation of the intercellular collision frequency, whereas adhesion efficiency remains nearly constant. Moreover, hydrodynamic shear affects the receptor specificity of activation-dependent platelet binding to LS174T cells, as evidenced by the transition from a P-selectin-independent/Arg-Gly-Asp (RGD)-dependent process at 100 s(-1) to a P-selectin/alpha(IIb)beta(3)-dependent interaction at 800 s(-1). This study demonstrates that platelet activation and a fluid-mechanical environment representative of the vasculature affect platelet-tumor cell adhesive interactions pertinent to the process of blood-borne metastasis.

Neutrophil enhancement of fibrin deposition under flow through platelet-dependent and -independent mechanisms

We examined the effect of adherent neutrophils on fibrin deposition under laminar flow conditions. Perfusion of recalcified citrated platelet-free plasma (PFP) over neutrophils adherent to fibrinogen-coated glass at a venous wall shear rate of 62.5 s(-1) for 15 minutes resulted in dense deposition of fibrin around each neutrophil, whereas fibrin deposition on glass alone was sparse. Fibrin deposition on neutrophils was markedly reduced by anti-CD18 or anti-CD11b or a higher shear rate (250 s(-1)). Significantly less fibrin was deposited around adherent fibrinogen-coated beads, indicating that nonspecific "cross-sectional capture" effects were not responsible for the massive fibrin deposition on neutrophils. Direct visualization of fibrin capture by neutrophils and elimination of fibrin deposition at 15 minutes by a factor XIIa inhibitor (50 microg/mL corn trypsin inhibitor [CTI]) or elastase/cathepsin G inhibitors (Methoxysuccinyl-Ala-Ala-Pro-Ala-Chloromethyl-Ketone/Z-Gly-Leu-Phe-CMK, 100 micromol/L) indicated that neutrophils can capture short fibrin strands flowing in recalcified PFP lacking CTI and can also promote thrombin generation through pathways attenuated by inhibitors of factor XIIa, elastase, and cathepsin G. When neutrophils were allowed to interact with platelets on a fibrinogen surface before perfusion of recalcified CTI-treated PFP, the fibrin deposition was observed to be dramatic compared with that over surfaces coated with platelets alone or neutrophils alone and compared with that formed on platelets adherent to collagen. This neutrophil promotion of platelet-mediated fibrin formation was attenuated by inhibitors of elastase or cathepsin G but not anti-tissue factor antibody. Neutrophils can interact with platelets via released proteases to increase platelet procoagulant activity and fibrin formation in CTI-treated plasma under the low-flow conditions expected in venous thrombosis or inflammation.