- Open Access
Investigation of the essential role of platelet-tumor cell interactions in metastasis progression using an agent-based model
Theoretical Biology and Medical Modelling volume 11, Article number: 17 (2014)
Metastatic tumors are a major source of morbidity and mortality for most cancers. Interaction of circulating tumor cells with endothelium, platelets and neutrophils play an important role in the early stages of metastasis formation. These complex dynamics have proven difficult to study in experimental models. Prior computational models of metastases have focused on tumor cell growth in a host environment, or prediction of metastasis formation from clinical data. We used agent-based modeling (ABM) to dynamically represent hypotheses of essential steps involved in circulating tumor cell adhesion and interaction with other circulating cells, examine their functional constraints, and predict effects of inhibiting specific mechanisms.
We developed an ABM of Early Metastasis (ABMEM), a descriptive semi-mechanistic model that replicates experimentally observed behaviors of populations of circulating tumor cells, neutrophils, platelets and endothelial cells while incorporating representations of known surface receptor, autocrine and paracrine interactions. Essential downstream cellular processes were incorporated to simulate activation in response to stimuli, and calibrated with experimental data. The ABMEM was used to identify potential points of interdiction through examination of dynamic outcomes such as rate of tumor cell binding after inhibition of specific platelet or tumor receptors.
The ABMEM reproduced experimental data concerning neutrophil rolling over endothelial cells, inflammation-induced binding between neutrophils and platelets, and tumor cell interactions with these cells. Simulated platelet inhibition with anti-platelet drugs produced unstable aggregates with frequent detachment and re-binding. The ABMEM replicates findings from experimental models of circulating tumor cell adhesion, and suggests platelets play a critical role in this pre-requisite for metastasis formation. Similar effects were observed with inhibition of tumor integrin αV/β3. These findings suggest that anti-platelet or anti-integrin therapies may decrease metastasis by preventing stable circulating tumor cell adhesion.
Circulating tumor cell adhesion is a complex, dynamic process involving multiple cell-cell interactions. The ABMEM successfully captures the essential interactions necessary for this process, and allows for in-silico iterative characterization and invalidation of proposed hypotheses regarding this process in conjunction with in-vitro and in-vivo models. Our results suggest that anti-platelet therapies and anti-integrin therapies may play a promising role in inhibiting metastasis formation.
Metastatic disease presents a significant burden in cancer treatment and could be considered perhaps the most critical aspect of oncology in terms of control of established cancers. Currently, strategies for treatment of metastatic cancer focus on limiting the growth and further extension of metastatic disease, either by removal surgically or using chemotherapy or radiation therapy to eradicate or limit existing sites. Given the clinical challenges of dealing with established metastases, there is an increasing focus in limiting the generation of metastatic foci through better understanding of the processes leading from a primary tumor to metastasis. It is well recognized that the process of metastasis is a complex, multi-step process involving tumor cell interactions with the host’s vascular and immune systems at multiple levels: a non-comprehensive list of the steps required to develop a metastasis include tumor shedding from the primary cancer, the lifecycle/dynamics of circulating tumor cells and the process of metastasis colonization. (Reviewed in[2–9]). While each of these processes are the subject of intense study, it clear that there is a defined sequence of events with dependencies and consequences between levels of processes. This suggests that being able to effectively modulate the generation of metastases requires an integrated view of the multiple events involved. Challenges to an integrated view of the formation of a metastasis include the rarity of these events, the variability of tumor metastatic and growth potential, and the heterogeneity of host organs for metastasis[1, 2, 4, 10–14].
We focus on the vital step between the conditions leading to the generation of circulating tumor cells from the primary tumor and the process of invasion of those circulating tumor cells: the initial adhesion of shed and circulating tumor cells to their potential target tissue beds. Though millions of tumor cells are detectable in the circulation of patients with primary tumors, very few have the ability to successfully adhere, invade and proliferate[4, 13]. Coagulation events are known to play a significant role, with one important recognized interaction being the role of platelet activation in circulating tumor cells adhesion to the endothelium[15–18]. Given the relationship between coagulation and inflammation, we posit that there are potential down-stream effects from the initial tumor-coagulation-endothelial interactions on endothelial permeability, tumor cell killing and the host tissue inflammatory/immune milieu, with consequent influences on the establishment of metastases. Additionally, there exists some evidence that interdiction at the early events of tumor cell adhesion has an effect on the establishment of metastases. Notably, aspirin, a potent inhibitor of platelet activation, has been found to reduce the incidence of metastasis in patients[10, 14]. In addition, inhibition of tumor cell expression of integrin αV/β3 has been shown to decrease melanoma, colorectal and breast cancer cell adhesion to endothelium and platelets, resulting in decreased number of metastases[16, 19–21]. The characterization of the establishment of metastases as a cascading series of events prompts our focus on the most upstream interactions associated with the establishment of metastases in order to provide a functional context for complementary investigations in the downstream processes of metastasis formation.
However, there are significant challenges to the experimental study of the establishment of metastases from circulating tumor cells. In particular, longitudinal observation of single cells in circulation to metastasis formation is technically difficult in animal models, and impossible in humans. This results in a host of potential variables extremely difficult to account for in experimental models. Identifying which of myriad genetic variations is responsible for a phenotype is resource intensive and potentially intractable without use of new methodologies.
These challenges can potentially be met by dynamic computational modeling to represent and interrogate relevant variables of existing mechanistic knowledge within a functional context of complex cellular interactions[22–25]. Dynamic modeling allows for variables to be adjusted in silico and resulting behaviors observed with more ease and at a higher degree of spatial and temporal resolution than can be achieved with standard biological models. This allows for more rapid consideration of the plausibility of potential mechanisms, discarding those clearly not correct and allowing experimental resources to be focused on the most plausible hypotheses[23, 26–29].
One method used for computational dynamic knowledge representation is agent-based modeling[30–35]. Agent-based models (ABMs) can be used to simulate complex interactions as they are made of populations of computational agents, mimicking cells, that follow programmed rules, in parallel, that regulate their interaction with the environment and one another. Variability in response to certain inputs and production of outputs simulates the diversity of cellular behavior in a complex environment. The effect of altering specific variables on the complex dynamics generated can be examined in simulation runs. The outputs of experiments are provided continuously in a visual format that can be compared to biological experiments.
We have developed a descriptive, first-generation agent-based computational model that incorporates observed cellular behaviors and phenomenon in order to simulate the basic dynamics of circulating tumor cell adhesion in the context of endothelial, neutrophil and platelet interactions: the Agent-Based Model of early metastasis (ABMEM). Circulating tumor cell adhesion involves recruitment of neutrophils and platelets, multiple cell-cell interactions, initiation of cellular processes by cytokines, and activation of the coagulation cascade. These processes culminate in the stable binding of tumor cells to endothelial cells, a necessary precursor for subsequent tumor cell invasion into the host organ. Though not a predictive model, the ABMEM allows us to propose which mechanisms are essential for stable tumor cell adhesion and thus may represent potential therapeutic targets for anti-metastasis therapy.
Overview of the Agent-Based Model of Early Metastasis (ABMEM)
The ABMEM integrates currently known mechanistic knowledge observed in published biological models of tumor, neutrophil, platelet and endothelial interactions (see Table 1 and the Materials and Methods for a list of components of the model). Development of the ABMEM was performed in an iterative manner, with successive layers of validation in regards to known behaviors, a procedure referred to as the Iterative Refinement Protocol[19, 28, 36–39]. Initial iterations of the ABMEM focused on producing face validity, the ability of the model to behave in a biologically plausible fashion[40, 41]. Subsequent iterations emphasized experimental validity through addition of mechanistic details if the existing model is unable to reproduce the behaviors of interest observed in experimental systems[42, 43].
The computational agents represent the four cell types experimentally shown to be involved in tumor cell adhesion: neutrophils, platelets, endothelial cells and the tumor cells themselves. Figure 1 displays the interactions between cell types included in the ABMEM. Agents have internal state variables representing molecular components of cells: receptors, cell-surface proteins, genes, and signaling molecules. The pathways modeled were abstractly represented with logic-based, algebraic rules: these are demonstrated in Figure 2 for platelet agents and Figure 3 for neutrophil agents. The level of available active surface receptors and ligands were represented as state variables, and a threshold necessary for activation was fit for each variable. More details of the ABMEM can be seen in Materials and Methods. Following initial development and calibration we used the ABMEM to predict the effects of inhibiting the different adhesion mechanisms on tumor cell adhesion in the terms of dependencies on neutrophil, platelet and endothelial interactions.
Calibration and initial validation of neutrophil, platelet and endothelial cell dynamics
The ABMEM was calibrated by confirming that the dynamics between neutrophils, platelets and endothelial cells behaved in a plausible fashion by binding and activating in expected manners in response to endothelial inflammation, thrombin generation and cytokine production. This step establishes face validity of the ABMEM, i.e. establishing that the model performs in an intuitively plausible fashion as compared to existing real world reference systems[41, 44]. Rates for signal molecule production and diffusion, and adhesion molecule binding levels for cell-cell interactions were adjusted to fit neutrophil rolling, endothelial binding, platelet binding and thrombin generation in appropriate timescales as observed in vitro. The desired plausible behavioral criteria included: 1) the maintenance of a dynamic, steady state of neutrophil rolling and platelet inactivation on non-inflamed endothelium[23, 27, 45]; 2) neutrophil and platelet activation in presence of inflamed endothelium[33, 39, 40, 46, 47]; 3) inhibition of platelet activation with inhibition of thromboxane or adenosine diphosphate (ADP) receptors[44, 48].
The baseline ABMEM calibration of platelet and neutrophil behavior was followed by the addition of tumor cell agents to the model. Further heuristic tuning of tumor cell adhesion thresholds and signaling molecule production was performed to achieve similarity to interaction steps observed in vitro and in vivo[40, 49–51]. This included circulating tumor cell adhesion to activated neutrophils and platelets, initial adhesion of tumor cells to the endothelium directly or through these mechanisms, and stable binding over time (see Additional file1 for similarity assessments). Figure 4 demonstrates the interactions between cell types over the course of a simulation run.
Simulation of platelet inhibition
The next simulations were performed to examine the effect of platelet inhibition on the observed dynamics of tumor cell adhesion. Specifically, we hypothesized that platelet inhibition would destabilize tumor cell binding to the endothelium, a pre-requisite step for metastatic colonization of a host organ. Two different types of platelet inhibition were examined, each corresponding to a currently clinically available anti-platelet drug: inhibition of thromboxane by aspirin and inhibition of platelet ADP receptors by clopidorel. These experiments were run for 1,000 iterations, representing a time period of approximately one hour, corresponding with the time-course of in-vitro adhesion and invasion observed in multiple tumor types[21, 28, 50, 52–54]. Simulations (n = 20 per condition) were run to simulate the baseline platelet activity, inhibition of thromboxane by aspirin, inhibition of ADP receptors by clopidogrel and inhibition of both pathways. Outcome measures were the percentage of tumor cells stably adhered to the endothelium at the end of each run (Figure 5). Simulation experiments with baseline platelets resulted in 95 ± 5% tumor adhesion vs. 5.5 ± 7% (p < 0.001) with thromboxane inhibition, 32 ± 26% (p < 0.05) with ADP receptor inhibition and 4 ± 5% (p < 0.001) with combined inhibition. Though both ADP receptor inhibition and thromboxane inhibition prevented platelet adhesion in the presence of only neutrophils and inflamed endothelium, thromboxane inhibition alone had significantly more effect on tumor cell adhesion. The inability to maintain integrin GpIIb activation through thromboxane resulted in failure of platelets to stably adhere to tumor cells and endothelium. In contrast, inhibition of the less stable P-selectin bond was partially compensated by binding via activated integrin GpIIb. These results suggest that aspirin may have a strong anti-metastatic effect through inhibition of stable tumor cell binding to activated platelets, while clopidogrel may have a smaller effect.
Evaluating the role of integrins: simulation of αV/β3 inhibition
Another mechanism of tumor metastasis inhibition identified in animal models is targeting of the integrin αV/β3, which is able to bind surface receptors of both platelets and neutrophils[21, 32, 55, 56]. Monoclonal antibodies to integrin αV/β3 have been used in clinical trials to improve progression-free survival of glioblastoma purportedly through inhibition of angiogenesis[1, 13, 40, 57–60]. Its efficacy against established metastatic disease, however, has been limited[61, 62]. We hypothesize another potential mechanism of integrin αV/β3 blockade in anti-metastastis therapy: preventing circulating tumor cell adhesion to neutrophils and platelets. As the interactions between tumor cells and these cells promote stable endothelium adhesion, we hypothesized that inhibiting these bonds would attenuate the subsequent adherence of the tumor cell complex to the endothelium. Simulations were performed where integrin αV/β3 was fully inhibited to mimic the effects of the potent monoclonal antibodies used in animal and clinical studies[63, 64]. These experiments were run for 1,000 steps to simulate a time course of 1 hour. Outcomes measured were stable tumor cell adhesion to endothelium and number of platelets bound to individual tumor cells. Varying degrees of inhibition of αν/β3 were examined, and 20 runs were performed for each condition. The stable binding of tumor cells was almost entirely inhibited by αV/β3 integrin suppression, though tumor cells were able to intermittently bind to PMNs through selectins (Figure 6). These simulations indicate that inhibition of αV/β3 may play an earlier role in prevention of metastasis formation than currently attempted in clinical trials by inhibiting stable circulating tumor cell adhesion.
Overview of tumor cell adhesion in the context of metastasis
The formation of tumor metastasis is a complex process involving multiple interactions between tumor cells and the potential target tissue environment. Each of these steps presents an opportunity for interdiction of tumor cells progressing towards metastasis. Tumor cell adhesion to the vasculature of distant organs is an essential step in this process that exemplifies the co-opting of local tissue factors in a fashion that promotes metastasis. It involves multiple cell types and mechanisms, all of which may play a role in the ability of a given tumor cell adhering to the endothelium for sufficient time to begin the next step of intravasation. This presents a rich field of potential therapeutic targets that may curtail metastatic dissemination before tumor cells are able to colonize distant organs.
Tumor cell adhesion involves binding to endothelial cells directly and indirectly via platelets and neutrophils. Tumor cells of various primary histologies are able to activate and bind directly to endothelial cells via integrins[40, 57, 59, 60]. Further studies showed that circulating tumor cells are able to activate and adhere to neutrophils through integrin interactions[54, 65–67]. Inhibition of neutrophil activation and integrin expression results in decreased tumor cell binding to neutrophils and subsequent endothelial adhesion within the microvasculature[40, 49, 50]. In addition to neutrophil binding, circulating tumor cell adhesion and subsequent metastasis formation is significantly enhanced by platelet binding through integrins[16, 18, 68]. Inhibition of this interaction results in decreased stable adhesion to endothelium[2–9]. In addition to direct binding, tumor cells can activate platelets or neutrophils through paracrine mechanisms such as interleukin-8 (IL-8) secretion, matrix metalloproteinase secretion and thrombin induction[50, 59, 69–71]. These pathways represent the major known mechanisms through which tumor cells are able to stably adhere to vasculature, an essential step in the formation of hematogenous metastases. A more detailed understanding of the interplay between these mechanisms is essential for the development of therapies to limit metastatic progression. However, the complexity of this system and its multiple redundancies limits the ability of traditional reductionist biological experiments to predict downstream effects of a perturbation of a single mechanism.
Complex, dynamic systems such as circulating tumor cell adhesion can be explored through the use of computer models. Incorporating these models with traditional experimentation allows for hypotheses to be discarded or advanced for further examination in a less resource-intensive manner. Significant efforts have been made in computational modeling the multiple steps of tumor growth, with the majority of efforts examining the local invasion and growth seen in primary tumors[3, 6–9]. Computational modeling of metastases has examined cell population dynamics based on a combination of proliferation, mutation and metastasis rates derived from clinical data to predict the number and time distribution of metastases from primary tumors,[72–76]. Others use mathematical models of cell growth and specific biologic processes such as matrix modeling and migration to predict metastasis growth over time[77–80]. Agent-based modeling of metastases has been employed to characterize the selective forces involved in the generation of circulating, potentially metastatic tumor cells[81, 82], implicate the role of host immunity in the generation of satellite metastases[1, 2, 4, 10–14], and examine the interactions of metastatic tumor cells with the host immune system for optimizing tumor vaccine delivery[15–18]. Other models of individual tumor cell interactions with host environment, employing knowledge of biomechanics and enzyme kinetics in a system of differential equations that describe integrin interactions and tumor cell growth patterns[83, 84]. In addition, several models of single circulating cell adhesion to vasculature have been developed incorporating biomechanic principles, fluid dynamics and knowledge of integrin activation (Reviewed in[22–25]). However, these models were designed to examine single cell adhesion, not the complex processes of tumor cell interaction with other cell types. Given the multiple interactions shown to be necessary for tumor cell adhesion, a computational model incorporating these various interactions can lead to a more detailed understanding of these early adhesion events.
This ABMEM is the first computational model to examine the multiple cell-cell interactions specifically related to adhesion of circulating tumor cells to vascular endothelium. This allows for an examination of population dynamics not possible with current single-cell models biomechanical models. Incorporation of intra-cellular pathways allows for a higher-resolution understanding of the functional signaling and receptor-binding events leading to stable tumor cell adhesion. The ABMEM, as an ABM, has a modular structure, both with respect to the agents/cells included, as well as the agents’ rules. This modular design allows for the addition of further cell types (as was done with the addition of tumor cells to the ABMEM) or molecular mechanisms as necessary to validate the model against prior experimental findings. This property of ABMs is consistent with the general iterative refinement protocol described in the development of dynamic computational models[23, 26–29]. This approach may promote a better understanding of the essential underlying biology by both limiting the initial set of components represented by the model and identifying when a particular model becomes insufficient to represent desired features and phenotypes of the real world system under investigation. Used in this adjunctive fashion, as a means of dynamic knowledge representation, models like the ABMEM can assist assessing the sufficiency of existing mechanistic hypotheses and aid in the generation of new hypotheses that may lead to the development of novel therapies to mitigate metastatic spread.
The ABMEM successfully represents the key mechanisms involved in tumor cell adhesion within the context of neutrophil and platelet interactions. Simulations of the inhibition of key platelet activation pathways support the concept that they may be necessary for stable tumor adhesion[16, 68]. The inhibition of thromboxane by aspirin resulted in significant inhibition of tumor cell adhesion by preventing activation of GpIIb/IIIa and P-selectin. This suggests that aspirin may reduce the formation of metastases in distant organs. These findings are supported by in vivo models of hepatocellular carcinoma and colon adenocarcinoma[85, 86]. Though thromboxane has many effects on platelets, integrin GpIIb/IIIa inhibition as modeled in the ABMEM is supported as a mechanism of reducing platelet adhesion to tumor cells, and thus decreased metastasis formation[31, 33–35]. The anti-metastastic effect of aspirin is also observed in retrospective reviews of randomized controlled trials, which revealed that aspirin was a predictor of decreased metastastic burden and subsequent death in patients with primary adenocarcinomas[10, 12, 14]. This suggests that therapies directed towards inhibition of platelet GpIIb/IIIa activation via thromboxane may potentially reduce metastatic spread through inhibition of adhesion.
The ABMEM also revealed clopidogrel as another potential inhibitor of tumor adhesion, though its mechanism of GpIIb/IIIa and P-selectin activation is dependent on calcium production and thus is upstream of thromboxane inhibition. This novel finding has not been described by in vivo models or clinical data, but these results from the ABMEM simulations suggest that a retrospective review of metastatic disease in the relatively large number of patients receiving clopidogrel post-cardiac and vascular procedures might be fruitful. Additionally, the combined effect of aspirin and clopidigrel did not significantly further reduce tumor cell adhesion in the ABMEM compared to either agent alone, as a portion of tumor cells were able to adhere to the endothelium without platelet binding. As with the effect of clopidogrel alone, this is a novel finding that may be subjected to experimental testing in the future.
The ABMEM also identified direct inhibition of integrin αν/β3 expressed by tumor cells themselves as potential mechanism of tumor cell adhesion, even in the presence of platelets. This resulted in intermittent binding of tumor cells to other cells through selectins, but failure of stable integrin-based bonds to form. As noted with combined anti-platelet inhibition, a small fraction of tumor cells were still able to stably adhere to the endothelium. This is supported by in vivo models of colon metastasis to liver, wherein antibodies to αV/β3 were able to significantly reduce, but not completely suppress, tumor cell adhesion[19, 28, 37, 39]. In addition, αV/β3 integrin expression alone was shown to enable non-tumor CHO cells to adhere and invade into to the liver[30, 32, 40, 41]. These results suggest that inhibitors of integrin αV/β3 may be potent therapies for prevention of metastases. Though clinical trials of the monoclonal antibody to integrin αV, intetumumab, have shown limited success with metastatic disease, trials of this agent in patients with early-stage tumors for prevention has not been carried out[36, 38, 41, 44]. The ABMEM demonstrated that integrin αV/β3 is a potential target for tumor metastasis prevention, and suggests a novel use of a therapeutic agent undergoing clinical trials.
The finding of an apparent saturation point in the effect of inhibition of platelet activation pathways on tumor cell adhesion is an example of the multi-component and interaction-dependent phenomena that can be examined using agent-based modeling. This finding suggests a mechanism by which tumor cells are able to adhere directly to endothelial cells via integrins, thus producing metastases despite potent inhibition of platelet adhesion within in vivo models[68, 87].
The ABMEM, as with all computational models, has limitations in its scope and the parameters included within its design. This model does not incorporate mechanisms not specifically related to adhesion or activation of the cell types. Mechanisms such as cell cycle and replication are not included as these occur in time scales well beyond the model timeframe. Energy consumption is also not incorporated as it is assumed to be un-constrained by the environment of oxygenated blood. In addition, the range of cell adhesion molecules expressed on various cell types is quite diverse and previously unknown interactions are rapidly being discovered. Many of these are currently only described in certain contexts, while others appear to play a lesser role in the dynamics in question. These mechanisms may be incorporated in future versions of the ABMEM if proven to be essential and within its scope. Finally, some of the biophysical properties of the cells are abstracted, such as flow rate and resulting shear stress incorporated as cell movement and degeneration of bonds over time. Biomechanical properties such as cell membrane deformity and adhesion forces between cells are not included as the model focuses on the functional effects of cell pathways rather than the specific properties of these parameters. Despite these limitations, the model accurately reflects the known mechanistic knowledge of circulating tumor cell adhesion, correlates with findings from biological models, and presents evidence to support anti-platelet therapies as anti-metastatic.
The ABMEM is a novel proof of concept model that provides insights into the dynamic processes resulting in metastasis. A potential extension of this model includes incorporation of information regarding tumor-specific cell adhesion molecules expression patterns. In addition, the model can be expanded to include the effects of early inflammation induced by neutrophils and adherent tumor cells on vascular endothelial disruption and subsequent extravasation. These effects may modulate the pre-metastatic niche and early metastatic tumor survival.
Materials and methods
The Materials & Methods section consists of a description of the ABM’s development and a description of the simulation experiments carried out with the model. The development followed the general process described in the Overview, Design Concepts, Details (ODD) protocol of Grimm, et al., modified for agent-based modeling of biological systems. Additionally, the iterative nature of model development emphasizes the successive addition of features to match more detailed observed variables to improve the model’s representation of the biological system. The ABM was implemented using NetLogo 5.0, which can be obtained online at http://ccl.northwestern.edu/netlogo/. Initial simulation experiments include testing the validity of the ABM to reproduce normal circulating neutrophil, platelet and vascular endothelial interactions. Subsequent experiments introduced tumor cell interactions with these populations to validate known behaviors, and finally the model was used to examine the effects of defects in platelet, neutrophil or tumor cell pathways.
Overview of ABM architecture
The ABM represents a two-dimensional patch of vascular endothelium that has been split open and laid flat into a toroidal square grid. Immobile endothelial cell agents along with circulating neutrophils, platelets and tumor cells populate the grid. Grid spaces within the model provide locations for agents, and possess environmental variables representing concentration of cytokines, reactive oxygen species and soluble proteins. The ABM incorporates spatial effect to examine the relationships between agents rather than specific histologic details. This model can be viewed as a cellular interaction network. Initial cell population ratios are based on physiologic ratios of platelets to neutrophils. The agent types and their internal variables are listed in Table 1. These variables are abstractions of molecular pathways manipulated by agents using logic-based and simple algebraic rules. They are updated and the resulting commands executed by each agent with each time step of the model. This ABM utilizes 1 time step to represent 1 second in order to reflect the timescale of the cell-cell interactions observed in vitro and in vivo. The complete code of the model can be found online upon publication.
Description of secreted mediators implemented as environmental variables
Much of cell-to-cell interactions occur through secretion and binding of mediators secreted into circulation. The ABM simulates this paracrine behavior through the use of environmental (“Patch”) variables representing thrombin, fibrin, ADP, MMP2, interleukin-8, and reactive oxygen species. These variables are listed in Table 1. They represent the amount of a given signaling molecule in a given area, affecting those cells within that grid space. These variables are “produced” by cellular agents acting on their specific rules, and “bound” by other agents acting on their set of rules. This results in increasing or decreasing levels of a given variable on a grid space with each time step. To simulate the spread of these mediators through circulation, the ABM utilizes NetLogo’s “Diffuse” function, which reduces the value of a given environmental variable by a set percentage and distributes this amount to neighboring spaces. These parameters are calibrated in relation to each other and the model’s timescale, and against the behavior of agents during calibration simulations. Their values are listed in the Additional file1.
Description of cell types implemented as agents
The ABM emphasizes the interactions between circulating cells and interactions with the endothelium in a time scale representative of cellular adhesion. These interactions are described below, and globally depicted in Figure 1.
Description of each cell type
Vascular endothelial cells
Vascular endothelial cell functions in the ABM are heavily abstracted, focusing on their role in circulating cell adhesion and response to inflammatory mediators. They are non-motile, but interact through paracrine means via environmental variables.
Their rules include:
Degradation of circulating thrombin by anti-thrombin 3. This is an abstraction of the degradation of thrombin by endothelial cells through multiple mechanisms, including activated protein C and tissue factor pathway inhibitor [91, 92].
Platelets play an important role in adhesion of circulating tumor cells to endothelium. They can bind to tumor cells during circulation, enhancing their binding to endothelial cells. In addition, platelets bound to activated endothelium can also bind circulating tumor cells. Activated platelets also produce thrombin, which trigger inactive platelets in a paracrine manner. At initialization, platelets are randomly distributed throughout the field in an inactive form.
Rules of circulation:
Platelets circulate quasi-randomly through the model at set rate.
Platelets circulating off the model-space are destroyed and replaced.
Platelets bound to a neutrophil or tumor cell will move with the target cell.
Rules of activation:
Activation of PAR-1, PAR-4 or P2Y12 results in increased intra-cellular calcium levels, triggering production of thromboxane, which further activates integrin GpIIb. GpIIb is deactivated at a set rate over time [30, 37, 95–98].
Rules of interaction with endothelium:
Rules of interaction with neutrophils:
Platelets are assigned, upon creation, a given binding threshold representing the amount of activated P-selectin necessary for binding to neutrophil integrins [112–114]. These thresholds are normally distributed across the population and calibrated within the model relative to other cell types. Platelets with activated P-selectin above their set threshold will bind to circulating neutrophils activated above their given threshold. They will disengage if levels of activation of either surface protein decreases below given thresholds.
Rules of interaction with tumor cells:
Platelets are assigned, as above, a threshold for GpIIb/IIIa activation and P-selectin binding necessary for tumor cell binding. These molecules bind with tumor cells expressing sufficient levels of αV/β3 in a manner similar to platelet-neutrophil interactions [5, 16, 115].
Neutrophils also play an important role in adhesion of circulating tumor cells to endothelium. They can bind to tumor cells during circulation, enhancing their binding to endothelial cells. In addition, neutrophils roll across and bind to inflamed endothelium. These bound neutrophils primed by exposure to inflammatory cytokines will produce reactive oxygen species, a processed termed “respiratory burst”. These neutrophils also release inflammatory cytokines and pro-coagulants, promoting further neutrophil and platelet activation. At initialization of the model, neutrophils are randomly distributed throughout the field in an inactive form.
Rules of circulation:
Neutrophils circulate quasi- randomly through the model at a set rate.
Rules of interaction with endothelium:
Neutrophils exposed to the inflammatory cytokine IL-8 will be “primed”, activating L-selectin and N-cadherin, rolling along endothelial cells at a slow speed [40, 99, 101–103]. L-selectin is decreased as the neutrophil rolls to represent shedding of this molecule, which diminishes further neutrophil adhesion by competitive inhibition [116, 117].
Primed neutrophils will also activate Macrophage-1 antigen (Mac-1 or integrin αM/β2) in response to IL-8 [23, 47, 118–120]. Rolling neutrophils with sufficient activation will bind to an endothelial cell in their patch if the endothelial cell is also sufficiently activated [121, 122].
Rules for interaction with tumor cells:
Neutrophils are assigned a threshold for binding to tumor cells representing activation of Mac-1 and Lymphocyte-Function Associated Antigen 1 (LFA-1 or integrin αL/β2). Sufficiently activate neutrophils bind tumor cells through αV/β3 integrin and ICAM-1, also represented by an activation threshold [49, 119, 123]. This is calibrated in relation to other cell-cell interactions within the model.
Rules for interaction with platelets:
Neutrophils will bind to platelets in a reciprocal relation described above.
Rules for activation of oxidative burst:
Neutrophils exposed to IL-8 or thrombin and bound to other cell types through β2 integrin will become progressively more activated, thus increasing both integrin expression and reactive oxygen species [99, 120, 127–129].
This process is partially inhibited by endothelial cells’ release of nitric oxide .
Once sufficiently activated above a set threshold, the neutrophil will produce reactive oxygen species until deactivated or killed from oxidative damage at a set threshold.
ROS production is capable of killing tumor cells and endothelial cells if the local concentration exceeds set thresholds for each cell type. These thresholds were set to a level that resulted in death of cell types over the timescale of in vitro and in vivo experiments (30 minutes to 4 hours) [116, 130–133].
Tumor cells in this simulation possess abstract representations of the major cell surface molecules experimentally associated with endothelial, platelet and neutrophil adhesion. Though modeled on melanoma, breast and colorectal tumor cells, they are generalizable to most tumor types. They also constitutively produce inflammatory cytokines, represented by IL-8, and induce thrombin generation. The specific internal details of the tumor cells are not modeled. At initialization of the model, tumor cells are randomly distributed throughout the field in an inactive form.
Rules of circulation:
Tumor cells circulate quasi- randomly through the model at a set rate.
Rules of interaction with endothelium:
Rules for interaction with neutrophils and platelets:
Tumor cells bind with neutrophils and platelets in reciprocal interactions as described above.
Process of Fitting and Calibration:
Each cell adhesion molecule is controlled by two sets of rules governing the activation of the molecules and the thresholds necessary for effective binding. To determine the optimal parameters for replicating experimentally observed behaviors between cell types, a range of values for each threshold was inputted (parameter sweeping), and effects on cell-cell adhesion in the time frame used in the reference experiments were quantified. These parameters were then adjusted to qualitatively replicate referenced results. Criteria for fitting included 1) replicating the observed adhesion interactions in the time frame observed experimentally and 2) replicating the effects of inhibition of each molecular component of the cell-cell interaction. Refer to Additional file1 for examples of model calibration. Adjustment of parameters over a broad range of values resulted in minor changes in quantitative effects (example in Additional file1: Figure S4).
The Agent-Based Model of Early Metastasis (ABMEM) provides an effective platform for hypothesis generation concerning the early events of metastasis formation, namely the initial adhesion of circulating tumor cell complexes in the target organs, and examination of potential anti-metastatic therapies. It allows for dynamic representation of current mechanistic knowledge of complex tumor cell interactions with the host environment. The ABMEM suggests that inhibition of platelet mechanisms of activation, mediated by thromboxane and ADP receptors, may strongly inhibit stable circulating tumor cell adhesion. In addition, the ABMEM suggests a novel use of inhibiting circulating tumor cell Integrin αV/β3 to also prevent stable adhesion to the endothelium. Inhibition of stable tumor cell adhesion, a critical step in the metastatic cascade, may lead to novel anti-metastasis therapies. These findings demonstrate the utility of Agent-Based Modeling in understanding phenomena that are difficult to examine with traditional experimental methods.
Agent Based Model
Agent Based Model of Early Metastasis
Talmadge JE, Fidler IJ: AACR centennial series: the biology of cancer metastasis: historical perspective. Cancer Res. 2010, 70: 5649-5669. 10.1158/0008-5472.CAN-10-1040.
Langley RR, Fidler IJ: Tumor cell-organ microenvironment interactions in the pathogenesis of cancer metastasis. J Biol Chem. 2007, 28: 297-321.
Kam Y, Rejniak KA, Anderson ARA: Cellular modeling of cancer invasion: integration of in silico and in vitro approaches. J Cell Physiol. 2011, 227: 431-438.
Maheswaran S, Haber DA: Circulating tumor cells: a window into cancer biology and metastasis. Curr Opin Genet Dev. 2010, 20: 96-99. 10.1016/j.gde.2009.12.002.
Zhang C, Liu Y, Gao Y, Shen J, Zheng S, Wei M, Zeng X: Modified heparins inhibit integrin αIIb/β3 mediated adhesion of melanoma cells to platelets in vitroand in vivo. J Pathol. 2009, 125: 2058-2065.
Wynn ML, Kulesa PM, Schnell S: Computational modelling of cell chain migration reveals mechanisms that sustain follow-the-leader behaviour. J R Soc Interface. 2012, 9: 1576-10.1098/rsif.2011.0726.
Frieboes HB: An Integrated Computational/Experimental Model of Tumor Invasion. Cancer Res. 2006, 66: 1597-1604. 10.1158/0008-5472.CAN-05-3166.
Brodland GW, Veldhuis JH: The Mechanics of Metastasis: Insights from a Computational Model. PLoS One. 2012, 7: e44281-10.1371/journal.pone.0044281.
Lefor AT: Computational Oncology. Jpn J Clin Oncol. 2011, 41: 937-947. 10.1093/jjco/hyr082.
Rothwell PM, Wilson M, Price JF, Belch JF, Meade TW, Mehta Z: Effect of daily aspirin on risk of cancer metastasis: a study of incident cancers during randomised controlled trials. Lancet. 2012, 379: 1591-1601. 10.1016/S0140-6736(12)60209-8.
Heiko Enderling LHPH: Immunoediting: evidence of the multifaceted role of the immune system in self-metastatic tumor growth. Theor Biol Med Model. 2012, 9: 31-10.1186/1742-4682-9-31.
Cole BF, Logan RF, Halabi S, Benamouzig R, Sandler RS, Grainge MJ, Chaussade S, Baron JA: Aspirin for the Chemoprevention of Colorectal Adenomas: meta-analysis of the Randomized Trials. J Natl Cancer Inst. 2009, 101: 256-266. 10.1093/jnci/djn485.
Allard WJ, Matera J, Miller MC, Repollet M, Connelly MC, Rao C, Tibbe AGJ, Uhr JW, Terstappen LWMM: Tumor Cells Circulate in the Peripheral Blood of All Major Carcinomas but not in Healthy Subjects or Patients With Nonmalignant Diseases. Clin Cancer Res. 2004, 10: 6897-6904. 10.1158/1078-0432.CCR-04-0378.
Jonsson F, Yin L, Lundholm C, Smedby KE, Czene K, Pawitan Y: Low-dose aspirin use and cancer characteristics: a population-based cohort study. Br J Cancer. 2013, 109: 1921-1925. 10.1038/bjc.2013.411.
Makó V, Czúcz J, Weiszhár Z, Herczenik E, Matkó J, Prohászka Z, Cervenak L: Proinflammatory activation pattern of human umbilical vein endothelial cells induced by IL-1β, TNF-α, and LPS. J Pathol. 2010, 77A: 962-970.
Lonsdorf AS, Krämer BF, Fahrleitner M, Schönberger T, Gnerlich S, Ring S, Gehring S, Schneider SW, Kruhlak MJ, Meuth SG, Nieswandt B, Gawaz M, Enk AH, Langer HF: Engagement of αIIbβ3 (GPIIb/IIIa) with ανβ3 integrin mediates interaction of melanoma cells with platelets: a connection to hematogenous metastasis. J Biol Chem. 2012, 287: 2168-2178. 10.1074/jbc.M111.269811.
Pennisi M, Pappalardo F, Palladini A, Nicoletti G, Nanni P, Lollini PL, Motta S: Modeling the competition between lung metastases and the immune system using agents. BMC Bioinforma. 2010, 11: S13-
Gay LJ, Felding-Habermann B: Contribution of platelets to tumour metastasis. Nat Rev Cancer. 2011, 11: 123-134. 10.1038/nrc3004.
Enns A, Korb T, Schlüter K, Gassmann P, Spiegel H-U, Senninger N, Mitjans F, Haier J: ανβ5-integrins mediate early steps of metastasis formation. Eur J Cancer. 2005, 41: 1065-1072. 10.1016/j.ejca.2004.12.031.
Gomes N, Vassy J, Lebos C, Arbeille B, Legrand C, Fauvel-Lafeve F: Breast adenocarcinoma cell adhesion to the vascular subendothelium in whole blood and under flow conditions: effects of ανβ3 and αIIbβ3 antagonists. Clin Exp Metastasis. 2004, 21: 553-561. 10.1007/s10585-004-3756-4.
Felding-Habermann B, O’Toole TE, Smith JW, Fransvea E, Ruggeri ZM, Ginsberg MH, Hughes PE, Pampori N, Shattil SJ, Saven A, Mueller BM: Integrin activation controls metastasis in human breast cancer. Proc Natl Acad Sci U S A. 2001, 98: 1853-1858. 10.1073/pnas.98.4.1853.
Li A, Varney ML, Valasek J, Godfrey M, Dave BJ, Singh RK: Autocrine role of interleukin-8 in induction of endothelial cell proliferation, survival, migration and MMP-2 production and angiogenesis. Angiogenesis. 2005, 8: 63-71. 10.1007/s10456-005-5208-4.
Zarbock A, Ley K: Neutrophil adhesion and activation under flow. Microcirculation. 2009, 16: 31-42. 10.1080/10739680802350104.
Michael J, Mitchell MRK: Computational and Experimental Models of Cancer Cell Response to Fluid Shear Stress. Frontiers in Oncology. 2013, 3: 44-
Vodovotz Y, Csete M, Bartels J, Chang S, An G: Translational Systems Biology of Inflammation. PLoS Comput Biol. 2008, 4: e1000014-10.1371/journal.pcbi.1000014.
Yeh M, Leitinger N, de Martin R, Onai N, Matsushima K, Vora DK, Berliner JA, Reddy ST: Increased transcription of IL-8 in endothelial cells is differentially regulated by TNF- α and oxidized phospholipids. Arterioscler Thromb Vasc Biol. 2001, 21: 1585-1591. 10.1161/hq1001.097027.
Lee D, Schultz JB, Knauf PA, King MR: Mechanical Shedding of L-selectin from the Neutrophil Surface during Rolling on Sialyl Lewis x under Flow. J Biol Chem. 2006, 282: 4812-4820.
Sheikh-Bahaei S, Hunt CA: Enabling Clearance Predictions to Emerge from In Silico Actions of Quasi-Autonomous Hepatocyte Components. Drug Metab Dispos. 2011, 39: 1910-1920. 10.1124/dmd.111.038703.
An G: Closing the Scientific Loop: bridging Correlation and Causality in the Petaflop Age. Sci Transl Med. 2010, 2: 34-41.
Wu C-C, Hwang T-L, Liao C-H, Kuo S-C, Lee F-Y, Teng C-M: The role of PAR4 in thrombin-induced thromboxane production in human platelets. Thromb Haemost. 2003, 90: 299-308.
Koshio O, Nagao T, Mabuchi A, Ono Y, Suzuki K: Apoptotic signaling in endothelial cells with neutrophil activation. Mol Cell Biochem. 2011, 363: 269-280.
Kikkawa H, Kaihou M, Horaguchi N, Uchida T, Imafuku H, Takiguchi A, Yamazaki Y, Koike C, Kuruto R, Kakiuchi T, Tsukada H, Takada Y, Matsuura N, Oku N: Role of integrin ανβ3 in the early phase of liver metastasis: PET and IVM analyses. Clin Exp Metastasis. 2002, 19: 717-725. 10.1023/A:1021356019563.
Jung U, Norman KE, Scharffetter-Kochanek K, Beaudet AL, Ley K: Transit time of leukocytes rolling through venules controls cytokine-induced inflammatory cell recruitment in vivo. J Clin Invest. 1998, 102: 1526-1533. 10.1172/JCI119893.
Nierodzik ML, Klepfish A, Karpatkin S: Role of platelets, thrombin, integrin IIb-IIIa, fibronectin and von Willebrand factor on tumor adhesion in vitro and metastasis in vivo. Thromb Haemost. 1995, 74: 282-290.
An G: Introduction of an agent-based multi-scale modular architecture for dynamic knowledge representation of acute inflammation. Theor Biol Med Model. 2008, 5: 11-10.1186/1742-4682-5-11.
Shapiro MJ, Weiss EJ, Faruqi TR, Coughlin SR: Protease-activated receptors 1 and 4 are shut off with distinct kinetics after activation by thrombin. J Biol Chem. 2000, 275: 25216-25221. 10.1074/jbc.M004589200.
Covic L, Gresser AL, Kuliopulos A: Biphasic kinetics of activation and signaling for PAR1 and PAR4 thrombin receptors in platelets. Biochemistry. 2000, 39: 5458-5467. 10.1021/bi9927078.
Heidenreich A, Rawal SK, Szkarlat K, Bogdanova N, Dirix L, Stenzl A, Welslau M, Wang G, Dawkins F, de Boer CJ, Schrijvers D: A randomized, double-blind, multicenter, phase 2 study of a human monoclonal antibody to human αν integrins (intetumumab) in combination with docetaxel and prednisone for the first-line treatment of patients with metastatic castration-resistant prostate cancer. Ann Oncol. 2012, 24: 329-336.
Darbousset R, Thomas GM, Mezouar S, Frere C, Bonier R, Mackman N, Renne T, Dignat-George F, Dubois C, Panicot-Dubois L: Tissue factor-positive neutrophils bind to injured endothelial wall and initiate thrombus formation. Blood. 2012, 120: 2133-2143. 10.1182/blood-2012-06-437772.
Strell C, Lang K, Niggemann B, Zaenker KS, Entschladen F: Surface molecules regulating rolling and adhesion to endothelium of neutrophil granulocytes and MDA-MB-468 breast carcinoma cells and their interaction. Cell Mol Life Sci. 2007, 64: 3306-3316. 10.1007/s00018-007-7402-6.
Klügl F: A validation methodology for agent-based simulations. 2008, New York, NY: Association for Computing Machinery, 39-43.
Balci O: Verification, Validation and Testing. 1998, New York: John Wiley & Sons, 335-396.
Balci O: A methodology for certification of modeling and simulation applications. ACM Trans Model Comput Simul. 2001, 11: 352-377. 10.1145/508366.508369.
Jin J, Daniel JL, Kunapuli SP: Molecular Basis for ADP-induced Platelet Activation. II. The p2y1 receptor mediates adp-induced intracellular calcium mobilization and shape change in platelets. J Biol Chem. 1998, 273: 2030-2034. 10.1074/jbc.273.4.2030.
Goerge T, Barg A, Schnaeker E-M, Poppelmann B, Shpacovitch V, Rattenholl A, Maaser C, Luger TA, Steinhoff M, Schneider SW: Tumor-derived matrix metalloproteinase-1 targets endothelial proteinase-activated receptor 1 promoting endothelial cell activation. Cancer Res. 2006, 66: 7766-7774. 10.1158/0008-5472.CAN-05-3897.
Frenette PS, Moyna C, Hartwell DW, Lowe JB, Hynes RO, Wagner DD: Platelet-endothelial interactions in inflamed mesenteric venules. Blood. 1998, 91: 1318-1324.
Ding Z-M, Babensee JE, Simon SI, Lu H, Perrard JL, Bullard DC, Dai XY, Bromley SK, Dustin ML, Entman ML, Smith CW, Ballantyne CM: Relative Contribution of LFA-1 and Mac-1 to Neutrophil Adhesion and Migration. J Immunol. 1999, 163: 5029-5038.
Nitelius E, Brantmark B, Fredholm B, Hedner U, Forshell GP, Wåhlin-Boll E, Melander A: Actions and interactions of acetylsalicylic acid, salicylic acid and diflunisal on platelet aggregation. Eur J Clin Pharmacol. 1984, 27: 165-168. 10.1007/BF00544040.
Zhang P, Ozdemir T, Chung C-Y, Robertson GP, Dong C: Sequential binding of αVβ3 and ICAM-1 determines fibrin-mediated melanoma capture and stable adhesion to CD11b/CD18 on neutrophils. J Immunol. 2011, 186: 242-254. 10.4049/jimmunol.1000494.
Huh SJ, Liang S, Sharma A, Dong C, Robertson GP: Transiently entrapped circulating tumor cells interact with neutrophils to facilitate lung metastasis development. Cancer Res. 2010, 70: 6071-6082. 10.1158/0008-5472.CAN-09-4442.
Liang S, Hoskins M, Khanna P, Kunz RF, Dong C: Effects of the Tumor-Leukocyte Microenvironment on Melanoma–Neutrophil Adhesion to the Endothelium in a Shear Flow. Cel Mol Bioeng. 2008, 1: 189-200. 10.1007/s12195-008-0016-8.
Enns A, Gassmann P, Schlüter K, Korb T, Spiegel H-U, Senninger N, Haier J: Integrins can directly mediate metastatic tumor cell adhesion within the liver sinusoids. J Gastrointest Surg. 2004, 8: 1049-1059. 10.1016/j.gassur.2004.08.016. discussion 1060
Gassmann P, Kang M-L, Mees ST, Haier J: In vivo tumor cell adhesion in the pulmonary microvasculature is exclusively mediated by tumor cell–endothelial cell interaction. BMC Cancer. 2010, 10: 177-10.1186/1471-2407-10-177.
ten Kate M, Aalbers AGJ, Sluiter W, Hofland LJ, Sonneveld P, van Jeekel J, Eijck CHJ: Polymorphonuclear leukocytes increase the adhesion of circulating tumor cells to microvascular endothelium. Anticancer Res. 2007, 27: 17-22.
Reinmuth N, Liu W, Ahmad SA, Fan F, Stoeltzing O, Parikh AA, Bucana CD, Gallick GE, Nickols MA, Westlin WF, Ellis LM: αvβ3 Integrin Antagonist S247 Decreases Colon Cancer Metastasis and Angiogenesis and Improves Survival in Mice. Cancer Res. 2003, 63: 2079-2087.
Wu YJ, Muldoon LL, Gahramanov S, Kraemer DF, Marshall DJ, Neuwelt EA: Targeting αV-integrins decreased metastasis and increased survival in a nude rat breast cancer brain metastasis model. J Neurooncol. 2012, 110: 27-36. 10.1007/s11060-012-0942-0.
Wang H: Tumor cell α3β1 integrin and vascular laminin-5 mediate pulmonary arrest and metastasis. J Cell Biol. 2004, 164: 935-941. 10.1083/jcb.200309112.
Desgrosellier JS, Cheresh DA: Integrins in cancer: biological implications and therapeutic opportunities. Nat Rev Cancer. 2010, 10: 9-22. 10.1038/nrc2748.
Kerk N, Strozyk EA, Poppelmann B, Schneider SW: The mechanism of melanoma-associated thrombin activity and von Willebrand factor release from endothelial cells. J Invest Dermatol. 2010, 130: 2259-2268. 10.1038/jid.2010.136.
Giavazzi R, Foppolo M, Dossi R, Remuzzi A: Rolling and adhesion of human tumor cells on vascular endothelium under physiological flow conditions. J Clin Invest. 1993, 92: 3038-3044. 10.1172/JCI116928.
Manegold C, Vansteenkiste J, Cardenal F, Schuette W, Woll PJ, Ulsperger E, Kerber A, Eckmayr J, von Pawel J: Randomized phase II study of three doses of the integrin inhibitor cilengitide versus docetaxel as second-line treatment for patients with advanced non-small-cell lung cancer. Investig New Drugs. 2013, 31: 175-10.1007/s10637-012-9842-6.
Bradley DA, Daignault S, Ryan CJ, DiPaola RS, Smith DC, Small E, Gross ME, Stein MN, Chen A, Hussain M: Cilengitide (EMD 121974, NSC 707544) in asymptomatic metastatic castration resistant prostate cancer patients: a randomized phase II trial by the prostate cancer clinical trials consortium. Investig New Drugs. 2010, 29: 1432-1440.
Trikha M, Zhou Z, Nemeth JA, Chen Q, Sharp C, Emmell E, Giles-Komar J, Nakada MT: CNTO 95, a fully human monoclonal antibody that inhibits? v integrins, has antitumor and antiangiogenic activityin vivo. Int J Cancer. 2004, 110: 326-335. 10.1002/ijc.20116.
DeNardo SJ, Burke PA, Leigh BR, O’Donnell RT, Miers LA, Kroger LA, Goodman SL, Matzku S, Jonczyk A, Lamborn KR, DeNardo GL: Neovascular targeting with cyclic RGD peptide (cRGDf-ACHA) to enhance delivery of radioimmunotherapy. Cancer Biother Radiopharm. 2000, 15: 71-79. 10.1089/cbr.2000.15.71.
Liang S, Slattery MJ, Dong C: Shear stress and shear rate differentially affect the multi-step process of leukocyte-facilitated melanoma adhesion. Exp Cell Res. 2005, 310: 282-292. 10.1016/j.yexcr.2005.07.028.
Wu QD, Wang JH, Condron C, Bouchier-Hayes D, Redmond HP: Human neutrophils facilitate tumor cell transendothelial migration. Am J Physiol Cell Physiol. 2001, 280: C814-C822.
Ozdemir T, Zhang P, Fu C, Dong C: Fibrin serves as a divalent ligand that regulates neutrophil-mediated melanoma cells adhesion to endothelium under shear conditions. Am J Physiol, Cell Physiol. 2012, 302: C1189-1201. 10.1152/ajpcell.00346.2011.
Coupland LA, Chong BH, Parish CR: Platelets and P-Selectin Control Tumor Cell Metastasis in an Organ-Specific Manner and Independently of NK Cells. Cancer Res. 2012, 72: 4662-4671. 10.1158/0008-5472.CAN-11-4010.
Goerge T, Kleinerüschkamp F, Barg A, Schnaeker E-M, Huck V, Schneider MF, Steinhoff M, Schneider SW: Microfluidic reveals generation of platelet-strings on tumor-activated endothelium. Thromb Haemost. 2007, 98: 283-286.
Gerotziafas GT, Galea V, Mbemba E, Khaterchi A, Sassi M, Baccouche H, Prengel C, van Dreden P, Hatmi M, Bernaudin JF, Elalamy I: Tissue factor over-expression by human pancreatic cancer cells BXPC3 is related to higher prothrombotic potential as compared to breast cancer cells MCF7. Thromb Res. 2012, 129: 779-786. 10.1016/j.thromres.2011.07.049.
Gil-Bernabé AM, Lucotti S, Muschel RJ: Coagulation and metastasis: what does the experimental literature tell us?. Br J Haematol. 2013, 162: 433-441. 10.1111/bjh.12381.
Haeno H, Gonen M, Davis MB, Herman JM, Iacobuzio-Donahue CA, Michor F: Computational Modeling of Pancreatic Cancer Reveals Kinetics of Metastasis Suggesting Optimum Treatment Strategies. Cell. 2012, 148: 362-375. 10.1016/j.cell.2011.11.060.
Newton PK, Mason J, Bethel K, Bazhenova LA, Nieva J, Kuhn P: A Stochastic Markov Chain Model to Describe Lung Cancer Growth and Metastasis. PLoS One. 2012, 7: e34637-10.1371/journal.pone.0034637.
Hanin L, Rose J, Zaider M: A stochastic model for the sizes of detectable metastases. J Theor Biol. 2006, 243: 407-417. 10.1016/j.jtbi.2006.07.005.
Hanin L, Korosteleva O: Does extirpation of the primary breast tumor give boost to growth of metastases? Evidence revealed by mathematical modeling. Math Biosci. 2010, 223: 133-141. 10.1016/j.mbs.2009.11.006.
Michor F, Nowak MA, Iwasa Y: Stochastic dynamics of metastasis formation. J Theor Biol. 2006, 240: 521-530. 10.1016/j.jtbi.2005.10.021.
Volker Haustein US: A dynamic model for tumour growth and metastasis formation. J Clin Bioinform. 2012, 2: 11-
Leon Chen L, Beck C: A superstatistical model of metastasis and cancer survival. Physica A: Stat Mechan Appl. 2008, 387: 3162-3172. 10.1016/j.physa.2008.01.116.
Withers HR, Lee SP: Modeling Growth Kinetics and Statistical Distribution of Oligometastases. Semin Radiat Oncol. 2006, 16: 111-119. 10.1016/j.semradonc.2005.12.006.
Anderson ARA, Chaplain MAJ, Newman EL, Steele RJC, Thompson AM: Mathematical Modelling of Tumour Invasion and Metastasis. J Theor Med. 2000, 2: 129-154. 10.1080/10273660008833042.
Chen J, Sprouffske K, Huang Q, Maley CC: Solving the Puzzle of Metastasis: The Evolution of Cell Migration in Neoplasms. PLoS One. 2011, 6: e17933-10.1371/journal.pone.0017933.
Athena C, Aktipis CCMJWP: Dispersal evolution in neoplasms: The role of disregulated metabolism in the evolution of cell motility. Cancer Prev Res (Phila). 2012, 5: 266-10.1158/1940-6207.CAPR-11-0004.
Geier F, Fengos G, Iber D: A computational analysis of the dynamic roles of talin, Dok1, and PIPKI for integrin activation. PLoS One. 2011, 6: e24808-10.1371/journal.pone.0024808.
Singh J, Hussain F, Decuzzi P: Role of differential adhesion in cell cluster evolution: from vasculogenesis to cancer metastasis. Computer Methods In Biomechanics and Biomedical Engineering. 2013
Futakuchi M, Ogawa K, Sano M, Tamano S, Takeshita F, Shirai T: Suppression of lung metastasis by aspirin but not indomethacin in an in vivo model of chemically induced hepatocellular carcinoma. Jpn J Cancer Res. 2002, 93: 1175-1181. 10.1111/j.1349-7006.2002.tb01220.x.
Yokoyama I, Hayashi S, Kobayashi T, Negita M, Yasutomi M, Uchida K, Takagi H: Prevention of experimental hepatic metastasis with thromboxane synthase inhibitor. Res Exp Med (Berl). 1995, 195: 209-215. 10.1007/BF02576790.
De Candia E, Hall SW, Rutella S, Landolfi R, Andrews RK, De Cristofaro R: Binding of thrombin to glycoprotein Ib accelerates the hydrolysis of Par-1 on intact platelets. J Biol Chem. 2001, 276: 4692-4698. 10.1074/jbc.M008160200.
Sellak H, Franzini E, Hakim J, Pasquier C: Reactive oxygen species rapidly increase endothelial ICAM-1 ability to bind neutrophils without detectable upregulation. Blood. 1994, 83: 2669-2677.
Holland JA, Meyer JW, Chang M-M, O’Donnell RW, Johnson DK, Ziegler LM: Thrombin Stimulated Reactive Oxygen Species Production in Cultured Human Endothelial Cells. Endothelium. 1998, 6: 113-121. 10.3109/10623329809072198.
Niu XF, Smith CW, Kubes P: Intracellular oxidative stress induced by nitric oxide synthesis inhibition increases endothelial cell adhesion to neutrophils. Circ Res. 1994, 74: 1133-1140. 10.1161/01.RES.74.6.1133.
Chan V, Chan TK: Antithrombin III in fresh and cultured human endothelial cells: A natural anticoagulant from the vascular endothelium. Thromb Res. 1979, 15: 209-213. 10.1016/0049-3848(79)90066-5.
van Hinsbergh VWM: Endothelium—role in regulation of coagulation and inflammation. Semin Immunopathol. 2012, 34: 93-10.1007/s00281-011-0285-5.
Irani K: Oxidant signaling in vascular cell growth, death, and survival: a review of the roles of reactive oxygen species in smooth muscle and endothelial cell mitogenic and apoptotic signaling. Circ Res. 2000, 87: 179-183. 10.1161/01.RES.87.3.179.
Lefer AM, Lefer DJ: The role of nitric oxide and cell adhesion molecules on the microcirculation in ischaemia-reperfusion. Cardiovasc Res. 1996, 32: 743-751. 10.1016/0008-6363(96)00073-9.
Jin J, Quinton TM, Zhang J, Rittenhouse SE, Kunapuli SP: Adenosine diphosphate (ADP)–induced thromboxane A2generation in human platelets requires coordinated signaling through integrin αIIbβ3 and ADP receptors. Blood. 2002, 99: 193-198. 10.1182/blood.V99.1.193.
Henriksen RA, Hanks VK: PAR-4 agonist AYPGKF stimulates thromboxane production by human platelets. Arterioscler Thromb Vasc Biol. 2002, 22: 861-866. 10.1161/01.ATV.0000014742.56572.25.
Falker K, Haglund L, Gunnarsson P, Nylander M, Lindahl TL, Grenegard M: Protease-activated receptor 1 (PAR1) signalling desensitization is counteracted via PAR4 signalling in human platelets. Biochem J. 2011, 436: 469-480. 10.1042/BJ20101360.
Jamieson GA, Tandon N, Greco N: Human Platelet Thrombin Receptors and the Two Receptor Model for Platelet Activation. Handb Platelet Physiol Pharmacol. 1999, 1: 21-37.
Guichard C, Pedruzzi E, Dewas C, Fay M, Pouzet C, Bens M, Vandewalle A, Ogier-Denis E, Gougerot-Pocidalo M-A, Elbim C: Interleukin-8-induced priming of neutrophil oxidative burst requires sequential recruitment of NADPH oxidase components into lipid rafts. J Biol Chem. 2005, 280: 37021-37032. 10.1074/jbc.M506594200.
Holinstat M, Voss B, Bilodeau ML, McLaughlin JN, Cleator J, Hamm HE: PAR4, but not PAR1, signals human platelet aggregation via Ca2+ mobilization and synergistic P2Y12 receptor activation. J Biol Chem. 2006, 281: 26665-26674. 10.1074/jbc.M602174200.
El Habbal MH, Smith L, Elliott MJ, Strobel S: Effect of heparin anticoagulation on neutrophil adhesion molecules and release of IL8: C3 is not essential. Cardiovasc Res. 1995, 30: 676-681. 10.1016/S0008-6363(95)00069-0.
Tedder TF: L-selectin-deficient mice have impaired leukocyte recruitment into inflammatory sites. J Exp Med. 1995, 181: 2259-2264. 10.1084/jem.181.6.2259.
Pillay J, Kamp VM, Pennings M, Oudijk E-J, Leenen LP, Ulfman LH, Koenderman L: Acute-phase concentrations of soluble fibrinogen inhibit neutrophil adhesion under flow conditions in vitro through interactions with ICAM-1 and MAC-1 (CD11b/CD18). J Thromb Haemost. 2013, 11: 1172-1182. 10.1111/jth.12250.
Tarandovskiy ID, Artemenko EO, Panteleev MA, Sinauridze EI, Ataullakhanov FI: Antiplatelet agents can promote two-peaked thrombin generation in platelet rich plasma: mechanism and possible applications. PLoS One. 2013, 8: e55688-10.1371/journal.pone.0055688.
Keuren JFW, Wielders SJH, Ulrichts H, Hackeng T, Heemskerk JWM, Deckmyn H, Bevers EM, Lindhout T: Synergistic effect of thrombin on collagen-induced platelet procoagulant activity is mediated through protease-activated receptor-1. Arterioscler Thromb Vasc Biol. 2005, 25: 1499-1505. 10.1161/01.ATV.0000167526.31611.f6.
Lentz BR: Exposure of platelet membrane phosphatidylserine regulates blood coagulation. Prog Lipid Res. 2003, 42: 423-438. 10.1016/S0163-7827(03)00025-0.
Vanschoonbeek K, Feijge MAH, Van Kampen RJW, Kenis H, Hemker HC, Giesen PLA, Heemskerk JWM: Initiating and potentiating role of platelets in tissue factor-induced thrombin generation in the presence of plasma: subject-dependent variation in thrombogram characteristics. J Thromb Haemost. 2004, 2: 476-484. 10.1111/j.1538-7933.2004.00618.x.
Bombeli T, Schwartz BR, Harlan JM: Adhesion of activated platelets to endothelial cells: evidence for a GPIIbIIIa-dependent bridging mechanism and novel roles for endothelial intercellular adhesion molecule 1 (ICAM-1), ανβ3 integrin, and GPIbα. J Exp Med. 1998, 187: 329-339. 10.1084/jem.187.3.329.
Kubes P, Suzuki M, Granger DN: Nitric oxide: an endogenous modulator of leukocyte adhesion. Proc Natl Acad Sci U S A. 1991, 88: 4651-4655. 10.1073/pnas.88.11.4651.
Alon R, Hammer DA, Springer TA: Lifetime of the P-selectin-carbohydrate bond and its response to tensile force in hydrodynamic flow. Nature. 1995, 374: 539-542. 10.1038/374539a0.
Lefer DJ, Jones SP, Girod WG, Baines A, Grisham MB, Cockrell AS, Huang PL, Scalia R: Leukocyte-endothelial cell interactions in nitric oxide synthase-deficient mice. Am J Physiol Heart Circ Physiol. 1999, 276: H1943-H1950.
Brown KK, Henson PM, Maclouf J, Moyle M, Ely JA, Worthen GS: Neutrophil-platelet adhesion: relative roles of platelet P-selectin and neutrophil β2 (DC18) integrins. Am J Respir Cell Mol Biol. 1998, 18: 100-110. 10.1165/ajrcmb.18.1.2314.
Piccardoni P, Sideri R, Manarini S, Piccoli A, Martelli N, de Gaetano G, Cerletti C, Evangelista V: Platelet/polymorphonuclear leukocyte adhesion: a new role for SRC kinases in Mac-1 adhesive function triggered by P-selectin: presented in part at the 42nd Annual Meeting of the American Society of Hematology, December 1-5, 2000, San Francisco, CA. Blood. 2001, 98: 108-116. 10.1182/blood.V98.1.108.
Patko Z, Csaszar A, Acsady G, Peter K, Schwarz M: Roles of Mac-1 and glycoprotein IIb/IIIa integrins in leukocyte–platelet aggregate formation: stabilization by Mac-1 and inhibition by GpIIb/IIIa blockers. Platelets. 2012, 23: 368-375. 10.3109/09537104.2011.625098.
Liu Y, Zhao F, Gu W, Yang H, Meng Q, Zhang Y, Yang H, Duan Q: The Roles of Platelet GPIIb/IIIa and αvβ3 Integrins during HeLa Cells Adhesion, Migration, and Invasion to Monolayer Endothelium under Static and Dynamic Shear Flow. J Biomed Biotech. 2009, 2009: 1-10.
Cao D, Boxer LA, Petty HR: Deposition of reactive oxygen metabolites onto and within living tumor cells during neutrophil-mediated antibody-dependent cellular cytotoxicity. J Cell Physiol. 1993, 156: 428-436. 10.1002/jcp.1041560227.
Kuijper PH, Gallardo Torres HI, van der Linden JA, Lammers JW, Sixma JJ, Zwaginga JJ, Koenderman L: Neutrophil adhesion to fibrinogen and fibrin under flow conditions is diminished by activation and L-selectin shedding. Blood. 1997, 89: 2131-2138.
DiVietro JA, Smith MJ, Smith BRE, Petruzzelli L, Larson RS, Lawrence MB: Immobilized IL-8 Triggers Progressive Activation of Neutrophils Rolling In Vitro on P-Selectin and Intercellular Adhesion Molecule-1. J Immunol. 2001, 167: 4017-4025. 10.4049/jimmunol.167.7.4017.
Strell C, Lang K, Niggemann B, Zaenker KS, Entschladen F: Neutrophil granulocytes promote the migratory activity of MDA-MB-468 human breast carcinoma cells via ICAM-1. Exp Cell Res. 2010, 316: 138-148. 10.1016/j.yexcr.2009.09.003.
A Wozniak WHBGAMMR: Interleukin-8 primes human neutrophils for enhanced superoxide anion production. Immunology. 1993, 79: 608-
Lee YS, Choi I, Ning Y, Kim NY, Khatchadourian V, Yang D, Chung HK, Choi D, LaBonte MJ, Ladner RD, Nagulapalli Venkata KC, Rosenberg DO, Petasis NA, Lenz H-J, Hong Y-K: Interleukin-8 and its receptor CXCR2 in the tumour microenvironment promote colon cancer growth, progression and metastasis. Br J Cancer. 2012, 106: 1833-1841. 10.1038/bjc.2012.177.
Fu C, Tong C, Wang M, Gao Y, Zhang Y, Lü S, Liang S, Dong C, Long M: Determining β2-integrin and intercellular adhesion molecule 1 binding kinetics in tumor cell adhesion to leukocytes and endothelial cells by a gas-driven micropipette assay. J Biol Chem. 2011, 286: 34777-34787. 10.1074/jbc.M111.281642.
Spicer JD, McDonald B, Cools-Lartigue JJ, Chow SC, Giannias B, Kubes P, Ferri LE: Neutrophils Promote Liver Metastasis via Mac-1-Mediated Interactions with Circulating Tumor Cells. Cancer Res. 2012, 72: 3919-3927. 10.1158/0008-5472.CAN-11-2393.
Zago AC, Simon DI, Wang Y, Sakuma M, Chen Z, Croce K, Ustinov V, Shi C, Martinez Filho EE: The importance of the interaction between leukocyte integrin Mac-1 and platelet glycoprotein Ib-a for leukocyte recruitment by platelets and for the inflammatory response to vascular injury. Arq Bras Cardiol. 2008, 90: 54-63.
Naum CC, Kaplan SS, Basford RE: Platelets and ATP prime neutrophils for enhanced O2- generation at low concentrations but inhibit O2- generation at high concentrations. J Leukoc Biol. 1991, 49: 83-89.
Detmers PA, Zhou D, Polizzi E, Thieringer R, Hanlon WA, Vaidya S, Bansal V: Role of stress-activated mitogen-activated protein kinase (p38) in β2-integrin-dependent neutrophil adhesion and the adhesion-dependent oxidative burst. J Immunol. 1998, 161: 1921-1929.
Babior BM, Lambeth JD, Nauseef W: The Neutrophil NADPH Oxidase. Arch Biochem Biophys. 2002, 397: 342-344. 10.1006/abbi.2001.2642.
Tsen A, Kirschenbaum LA, LaRow C, Khan R, Kurtz S, Bansal S, Astiz ME: The effect of anticoagulants and the role of thrombin on neutrophil-endothelial cell interactions in septic shock. Shock. 2009, 31: 120-124. 10.1097/SHK.0b013e318180ffe9.
Jenkins AL, Howells GL, Scott E, Le Bonniec BF, Curtis MA, Stone SR: The response to thrombin of human neutrophils: evidence for two novel receptors. J Cell Sci. 1995, 108 (Pt 9): 3059-3066.
Clancy RM, Leszczynska-Piziak J, Abramson SB: Nitric oxide, an endothelial cell relaxation factor, inhibits neutrophil superoxide anion production via a direct action on the NADPH oxidase. J Clin Invest. 1992, 90: 1116-1121. 10.1172/JCI115929.
Lichtenstein A: Neutrophil-mediated nonoxidative tumor lysis stimulated by high concentrations of phorbol myristate acetate. Clin Immunol Immunopathol. 1988, 47: 296-309. 10.1016/S0090-1229(88)80008-4.
Lichtenstein A, Kahle J: Anti-tumor effect of inelammatory neutrophils: Characteristics ofin vivo generation andin vitro tumor cell lysis. J Pathol. 1985, 35: 121-127.
Bruchelt G, Handgretinger R, Kimmig A, Goeke B, Siedner R, Reisfeld RA, Niethammer D, Treuner J: Effects of granulocytes on human neuroblastoma cells measured by chemiluminescence and chromium-51 release assay. J Biolumin Chemilumin. 1989, 3: 93-96. 10.1002/bio.1170030212.
De Larco JE, Wuertz BRK, Rosner KA, Erickson SA, Gamache DE, Manivel JC, Furcht LT: A Potential Role for Interleukin-8 in the Metastatic Phenotype of Breast Carcinoma Cells. Am J Pathol. 2001, 158: 639-10.1016/S0002-9440(10)64005-9.
Ning Y, Manegold PC, Hong YK, Zhang W, Pohl A, Lurje G, Winder T, Yang D, LaBonte MJ, Wilson PM, Ladner RD, Lenz H-J: Interleukin-8 is associated with proliferation, migration, angiogenesis and chemosensitivity in vitro and in vivo in colon cancer cell line models. Int J Cancer. 2010, 128: 2038-2049.
Waugh DJJ, Wilson C: The interleukin-8 pathway in cancer. Clin Cancer Res. 2008, 14: 6735-6741. 10.1158/1078-0432.CCR-07-4843.
Horowitz NA, Blevins EA, Miller WM, Perry AR, Talmage KE, Mullins ES, Flick MJ, Queiroz KCS, Shi K, Spek CA, Conway EM, Monia BP, Weiler H, Degen JL, Palumbo JS: Thrombomodulin is a determinant of metastasis through a mechanism linked to the thrombin binding domain but not the lectin-like domain. Blood. 2011, 118: 2889-2895. 10.1182/blood-2011-03-341222.
Zigler M, Kamiya T, Brantley EC, Villares GJ, Bar-Eli M: PAR-1 and thrombin: the ties that bind the microenvironment to melanoma metastasis. Cancer Res. 2011, 71: 6561-6566. 10.1158/0008-5472.CAN-11-1432.
Abecassis J, Millon-Collard R, Klein-Soyer C, Nicora F, Fricker JP, Beretz A, Eber M, Muller D, Cazenave JP: Adhesion of human breast cancer cell line MCF-7 to human vascular endothelial cells in culture. Enhancement by activated platelets. Int J Cancer. 1987, 40: 525-531. 10.1002/ijc.2910400416.
Al-Mondhiry H, McGarvey V, Leitzel K: Interaction of human tumor cells with human platelets and the coagulation system. Thromb Haemost. 1983, 50: 726-730.
Schulze EB, Hedley BD, Goodale D, Postenka CO, Al-Katib W, Tuck AB, Chambers AF, Allan AL: The thrombin inhibitor Argatroban reduces breast cancer malignancy and metastasis via osteopontin-dependent and osteopontin-independent mechanisms. Breast Cancer Res Treat. 2008, 112: 243-254. 10.1007/s10549-007-9865-4.
Marchiò S, Soster M, Cardaci S, Muratore A, Bartolini A, Barone V, Ribero D, Monti M, Bovino P, Sun J, Giavazzi R, Asioli S, Cassoni P, Capussotti L, Pucci P, Bugatti A, Rusnati M, Pasqualini R, Arap W, Bussolino F: A complex of α (6) integrin and E-cadherin drives liver metastasis of colorectal cancer cells through hepatic angiopoietin-like 6. EMBO Mol Med. 2012, 4: 1156-1175. 10.1002/emmm.201101164.
All contributors were involved sufficiently involved in the design of the experiment to be listed as an author. This work was supported in part by the University of Chicago Digestive Disease Research Core Center (P30 DK42086), which provided funding to GA. It had no role in the design, in the collection, analysis, and interpretation of data; in the writing of the manuscript; or in the decision to submit the manuscript for publication.
The authors declare that they have no competing interests.
AU and GA conceived of and participated in the design of the model, and helped draft the manuscript. SCW, SG and RRW participated in the design of the model, and helped revise the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Figure S1: Calibration of Neutrophils binding to endothelium. Calibration of Platelet Binding to Neutrophils in circulation. Figure S3. Calibration of tumor cell binding to platelets in circulation. Figure S4. Representative parameter sweeps. (PDF 181 KB)
Authors’ original submitted files for images
About this article
Cite this article
Uppal, A., Wightman, S.C., Ganai, S. et al. Investigation of the essential role of platelet-tumor cell interactions in metastasis progression using an agent-based model. Theor Biol Med Model 11, 17 (2014). https://doi.org/10.1186/1742-4682-11-17
- Agent-based modeling
- Computational modeling