Interpreting scratch assays using pair density dynamics and approximate Bayesian computation

Quantifying the impact of biochemical compounds on collective cell spreading is an essential element of drug design, with various applications including developing treatments for chronic wounds and cancer. Scratch assays are a technically simple and inexpensive method used to study collective cell spreading; however, most previous interpretations of scratch assays are qualitative and do not provide estimates of the cell diffusivity, D, or the cell proliferation rate,l. Estimating D and l is important for investigating the efficacy of a potential treatment and provides insight into the mechanism through which the potential treatment acts. While a few methods for estimating D and l have been proposed, these previous methods lead to point estimates of D and l, and provide no insight into the uncertainty in these estimates. Here, we compare various types of information that can be extracted from images of a scratch assay, and quantify D and l using discrete computational simulations and approximate Bayesian computation. We show that it is possible to robustly recover estimates of D and l from synthetic data, as well as a new set of experimental data. For the first time, our approach also provides a method to estimate the uncertainty in our estimates of D and l. We anticipate that our approach can be generalized to deal with more realistic experimental scenarios in which we are interested in estimating D and l, as well as additional relevant parameters such as the strength of cell-to-cell adhesion or the strength of cell-to-substrate adhesion.

Regulatory interplay between pap operons in uropathogenic Escherichia coli

Pathogenic bacteria have a large repertoire of surface organelles involved in adherence, motility and protein export, but how individual bacteria co-ordinate surface organelle expression to prevent interference and excessive immune stimulation is unclear. Phase variation is a mechanism by which expression of surface factors is limited to a fraction of the bacterial population; however, the presence of multiple homologous surface structures controlled by related mechanisms and regulators antagonizes the independent expression achieved by phase variation. To investigate whether other mechanisms have evolved to sort out the bacterial cell surface, we examined regulatory cross-talk between multiple phase-variable pyelonephritis-associated pili (pap) operons in Escherichia coli isolates associated with urinary tract infections. Allelic variation identified in the regulatory regions and regulators acts synergistically to limit coexpression of homologous fimbrial operons. In particular, there is evidence that papI is under positive selection and PapI variants displayed differences in their capacity to activate related pap operons. Alleles of the high-affinity binding site for PapB were shown to contain a variable number of (T/A)3 repeats occurring every 9 bp that altered the sensitivity of pap operon activation. Taken together with other examples of surface organelle cross-talk, we illustrate how this regulation could promote sequential expression.

Genotoxicity of inorganic lead salts and disturbance of microtubule function

Lead compounds are known genotoxicants, principally affecting the integrity of chromosomes. Lead chloride and lead acetate induced concentration-dependent increases in micronucleus frequency in V79 cells, starting at 1.1 μM lead chloride and 0.05 μM lead acetate. The difference between the lead salts, which was expected based on their relative abilities to form complex acetato-cations, was confirmed in an independent experiment. CREST analyses of the micronuclei verified that lead chloride and acetate were predominantly aneugenic (CREST-positive response), which was consistent with the morphology of the micronuclei (larger micronuclei, compared with micronuclei induced by a clastogenic mechanism). The effects of high concentrations of lead salts on the microtubule network of V79 cells were also examined using immunofluorescence staining. The dose effects of these responses were consistent with the cytotoxicity of lead(II), as visualized in the neutral-red uptake assay. In a cell-free system, 20-60 μM lead salts inhibited tubulin assembly dose-dependently. The no-observed-effect concentration of lead(II) in this assay was 10 μM. This inhibitory effect was interpreted as a shift of the assembly/disassembly steady-state toward disassembly, e.g., by reducing the concentration of assembly-competent tubulin dimers. The effects of lead salts on microtubule-associated motor-protein functions were studied using a kinesin-gliding assay that mimics intracellular transport processes in vitro by quantifying the movement of paclitaxel-stabilized microtubules across a kinesin-coated glass surface. There was a dose-dependent effect of lead nitrate on microtubule motility. Lead nitrate affected the gliding velocities of microtubules starting at concentrations above 10 μM and reached half-maximal inhibition of motility at about 50 μM. The processes reported here point to relevant interactions of lead with tubulin and kinesin at low dose levels.

Chromosomal genotoxicity of nitrobenzene and benzonitrile

In order to investigate the chromosomal genotoxicity of nitrobenzene and benzonitrile, we studied the induction of micronuclei (MN) by these test compounds in V79 cells, as well as effects on the formation and stability of microtubules and on motor protein functions. No cytotoxicity was seen in V79 cell cultures in terms of Neutral red uptake after 18 h treatment with up to 1 mM nitrobenzene or 1 mM benzonitrile. Subsequently, a concentration range up to 100 μM was used in the experiments on induction of MN. Both test compounds exhibit a weak, but definitely positive test result compared to the solvent (DMSO) control. Minimal effect concentrations of nitrobenzene and benzonitrile appeared as low as 0.01 μM, and no-effect-concentrations were between 0.001 and 0.005 μM. Clearly enhanced MN rates were found at 0.1 μM and higher. Both, nitrobenzene and benzonitrile, induced mostly kinetochor (CREST)-positive micronuclei, thus characterising the chromosomal effects as aneugenic. In cell-free assays, a slight effect on tubulin assembly was observed at 1 mM nitrobenzene without addition of DMSO. Higher concentrations (5 mM) led to secondary effects. In presence of 1% DMSO, nitrobenzene exerted no detectable effect on tubulin assembly up to the solubility limit in water of about 15 mM. For benzonitrile in presence of DMSO, a clear dose-response of inhibition of tubulin assembly at 37°C was seen above the no-effect-concentration of 2 mM, with an IC50 of 13 mM and protein denaturation starting above a level of about 20 mM. The nature of the effects of nitrobenzene and benzonitrile on the association of tubulin to form microtubules was confirmed by electron microscopy. Treatment by either 5 mM nitrobenzene or 13 mM benzonitrile plus 1% DMSO left the microtubular structure intact whereas 5 mM nitrobenzene, in absence of DMSO, led to irregular cluster formations. The experiments demonstrate that both nitrobenzene and benzonitrile, in millimolar concentration ranges, may lead to interference with tubulin assembly in a cell-free system. The functionality of the tubulin-kinesin motor protein system was assessed using the microtubule gliding assay. Nitrobenzene affected the gliding velocity in a concentration-dependent manner, starting at about 7.5 μM and reaching complete inhibition of motility at 30 μM, whereas benzonitrile up to 200 μM did not affect the kinesin-driven gliding velocity. The micronucleus assay data demonstrate a chromosomal endpoint of genotoxicity of nitrobenzene and benzonitrile. Aneugenic effects of both compounds occur at remarkably low concentrations, with lowest-effect-concentrations being 0.1 μM. This points to the relevance of interactions with the cellular spindle apparatus.

Genotoxicity of inorganic mercury salts based on disturbed microtubule function

This study investigated the hypothesis that the chromosomal genotoxicity of inorganic mercury results from interaction(s) with cytoskeletal proteins. Effects of Hg2+ salts on functional activities of tubulin and kinesin were investigated by determining tubulin assembly and kinesin-driven motility in cell-free systems. Hg2+ inhibits microtubule assembly at concentrations above 1 μM, and inhibition is complete at about 10 μM. In this range, the tubulin assembly is fully (up to 6 μM) or partially (∼6-10 μM) reversible. The inhibition of tubulin assembly by mercury is independent of the anion, chloride or nitrate. The no-observed-effect- concentration for inhibition of microtubule assembly in vitro was 1 μM Hg2+, the IC50 5.8 μM. Mercury(II) salts at the IC 50 concentrations partly inhibiting tubulin assembly did not cause the formation of aberrant microtubule structures. Effects of mercury salts on the functionality of the microtubule motility apparatus were studied with the motor protein kinesin. By using a "gliding assay" mimicking intracellular movement and transport processes in vitro, HgCl2 affected the gliding velocity of paclitaxel-stabilised microtubules in a clear dose-dependent manner. An apparent effect is detected at a concentration of 0.1 μM and a complete inhibition is reached at 1 μM. Cytotoxicity of mercury chloride was studied in V79 cells using neutral red uptake, showing an influence above 17 μM HgCl2. Between 15 and 20 μM HgCl2 there was a steep increase in cell toxicity. Both mercury chloride and mercury nitrate induced micronuclei concentration-dependently, starting at concentrations above 0.01 μM. CREST analyses on micronuclei formation in V79 cells demonstrated both clastogenic (CREST-negative) and aneugenic effects of Hg2+, with some preponderance of aneugenicity. A morphological effect of high Hg2+ concentrations (100 μM HgCl2) on the microtubule cytoskeleton was verified in V79 cells by immuno-fluorescence staining. The overall data are consistent with the concept that the chromosomal genotoxicity could be due to interaction of Hg2+ with the motor protein kinesin mediating cellular transport processes. Interactions of Hg 2+ with the tubulin shown by in vitro investigations could also partly influence intracellular microtubule functions leading, together with the effects on the kinesin, to an impaired chromosome distribution as shown by the micronucleus test.

Disturbed microtubule function and induction of micronuclei by chelate complexes of mercury(II)

Interactions of mercury(II) with the microtubule network of cells may lead to genotoxicity. Complexation of mercury(II) with EDTA is currently being discussed for its employment in detoxification processes of polluted sites. This prompted us to re-evaluate the effects of such complexing agents on certain aspects of mercury toxicity, by examining the influences of mercury(II) complexes on tubulin assembly and kinesin-driven motility of microtubules. The genotoxic effects were studied using the micronucleus assay in V79 Chinese hamster fibroblasts. Mercury(II) complexes with EDTA and related chelators interfered dose-dependently with tubulin assembly and microtubule motility in vitro. The no-effect-concentration for assembly inhibition was 1 μM of complexed Hg(II), and for inhibition of motility it was 0.05 μM, respectively. These findings are supported on the genotoxicity level by the results of the micronucleus assay, with micronuclei being induced dose-dependently starting at concentrations of about 0.05 μM of complexed Hg(II). Generally, the no-effect-concentrations for complexed mercury(II) found in the cell-free systems and in cellular assays (including the micronucleus test) were identical with or similar to results for mercury tested in the absence of chelators. This indicates that mercury(II) has a much higher affinity to sulfhydryls of cytoskeletal proteins than to this type of complexing agents. Therefore, the suitability of EDTA and related compounds for remediation of environmental mercury contamination or for other detoxification purposes involving mercury has to be questioned.

Interaction of metal salts with cytoskeletal motor protein systems

Interactions of chemicals with the microtubular network of cells may lead to genotoxicity. Micronuclei (MN) might be caused by interaction of metals with tubulin and/or kinesin. The genotoxic effects of inorganic lead and mercury salts were studied using the MN assay and the CREST analysis in V79 Chinese hamster fibroblasts. Effects on the functional activity of motor protein systems were examined by measurement of tubulin assembly and kinesin-driven motility. Lead and mercury salts induced MN dose-dependently. The no-effect-concentration for MN induction was 1.1 μM PbCl2, 0.05 μM Pb(OAc)2 and 0.01 μM HgCl2. The in vitro results obtained for PbCl2 correspond to reported MN induction in workers occupationally exposed to lead, starting at 1.2 μM Hg(II) (Vaglenov et al., 2001, Environ. Health Perspect. 109, 295-298). The CREST Analysis indicate aneugenic effects of Pb(II) and aneugenic and additionally clastogenic effects of Hg(II). Lead (chloride, acetate, and nitrate) and mercury (chloride and nitrate) interfered dose-dependently with tubulin assembly in vitro. The no-effect-concentration for lead salts in this assay was 10 μM. Inhibition of tubulin assembly by mercury started at 2 μM. The gliding velocity of microtubules along immobilised kinesin molecules was affected by 25 μM Pb(NO3)2 and 0.1 μM HgCl2 in a dose-dependent manner. Our data support the hypothesis that lead and mercury genotoxicity may result, at least in part, via disturbance of chromosome segregation via interaction with cytoskeletal proteins.

Tetraspanins as regulators of the tumour microenvironment : implications for metastasis and therapeutic strategies

One of the hallmarks of cancer is the ability to activate invasion and metastasis (Hanahan et al., 2011). Cancer morbidity and mortality are largely related to the spread of the primary, localised tumour to adjacent and distant sites (Pantel et al., 2004). Appropriate management and treatment decisions of predicting metastatic disease at the time of diagnosis is thus crucial, which supports better understanding of the metastatic process. There are common events that occur during metastasis: dissociation from the primary tumour mass, reorganisation/remodelling of extracellular matrix, cell migration, recognition and transversal of endothelial cells and the vascular circulation and lodgement and proliferation within ectopic stroma (Wells, 2006). One of the key and initial events is the increased capability of cancer cells to move, escaping the regulation of normal physiological control. The cellular cytoskeleton plays an important role in cancer cell motility and active cytoskeletal rearrangement can result in metastatic disease. This active change in cytoskeletal dynamics results in manipulation of plasma membrane and cellular balance between cellular adhesion and motility which in turn determines cancer cell movement. Members of the tetraspanins play important roles in regulation of cancer migration and cancer-endothelial cell interactions, which are critical for cancer invasion and metastasis. Their involvements in active cytoskeletal dynamics, cancer metastasis and potential clinical application will be discussed in this review. In particular, tetraspanin member, CD151, is highlighted for its major role in cancer invasion and metastasis

Menstrual effluent induces epithelial-mesenchymal transitions in mesothelial cells

Background Menstrual effluent affects mesothelial cell (MC) morphology. We evaluated whether these changes were consistent with epithelial-mesenchymal transitions (EMT). Methods Monolayer cultures of MC were incubated overnight in conditioned media, prepared from cells isolated form menstrual effluent, with or without kinase and ATP inhibitors. Changes in cell morphology were monitored using time-lapse video microscopy and immunohistochemistry. Effects on the expression of EMT-associated molecules were evaluated using real-time RT-PCR and/or Western blot analysis. Results Incubation in conditioned media disrupted cell-cell contacts, and increased MC motility. The changes were reversible. During the changes the distribution of cytokeratins, fibrillar actin and α-tubulin changed. Sodium azide, an inhibitor of ATP production, and Genistein, a general tyrosine kinase inhibitor, antagonized these effects. Wortmannin, a phosphatidylinositol 3-kinase inhibitor, and SU6656, an Src tyrosine kinase inhibitor, only partially antagonized the effect. The expression of Snail and vimentin was markedly up-regulated, whereas the expression of E-cadherin was decreased and cytokeratins were altered. Conclusions In MC, menstrual effluent initiates a reversible, energy-dependent transition process from an epithelial to a mesenchymal phenotype. Involvement of the (Src) tyrosine kinase signalling pathway and the changes in the expression of cytokeratins, Snail, vimentin and E-cadherin demonstrate that the morphological changes are EMT.

Touch, grasp, deliver and control : functional cross-talk between microtubules and cell adhesions

Cross-talk between microtubule networks and sites of cell-matrix and cell-cell adhesion has profound impact on these structures and is essential for proper cell organization, polarization and motility. Components of adhesion sites can interact directly with microtubules or with proteins that specifically associate with microtubule plus ends and minus ends and in this way capture, stabilize or destabilize microtubules. In their turn, microtubules can serve as routes for delivery of structural and regulatory factors that control adhesion site turnover. In addition, the microtubule lattice or growing microtubule plus ends can serve as diffusional sinks that accumulate and scaffold regulatory molecules, thereby affecting their activity in the vicinity of adhesions. Combination of these mechanisms underlies the functional co-operation between microtubules and adhesion sites and defines their dynamic behavior.

Choosing an appropriate modelling framework for analysing multispecies co-culture cell biology experiments

In vitro cell biology assays play a crucial role in informing our understanding of the migratory, proliferative and invasive properties of many cell types in different biological contexts. While mono-culture assays involve the study of a population of cells composed of a single cell type, co-culture assays study a population of cells composed of multiple cell types (or subpopulations of cells). Such co-culture assays can provide more realistic insights into many biological processes including tissue repair, tissue regeneration and malignant spreading. Typically, system parameters, such as motility and proliferation rates, are estimated by calibrating a mathematical or computational model to the observed experimental data. However, parameter estimates can be highly sensitive to the choice of model and modelling framework. This observation motivates us to consider the fundamental question of how we can best choose a model to facilitate accurate parameter estimation for a particular assay. In this work we describe three mathematical models of mono-culture and co-culture assays that include different levels of spatial detail. We study various spatial summary statistics to explore if they can be used to distinguish between the suitability of each model over a range of parameter space. Our results for mono-culture experiments are promising, in that we suggest two spatial statistics that can be used to direct model choice. However, co-culture experiments are far more challenging: we show that these same spatial statistics which provide useful insight into mono-culture systems are insuffcient for co-culture systems. Therefore, we conclude that great care ought to be exercised when estimating the parameters of co-culture assays.

Association of heparan sulfate proteoglycans SDC1 and SDC4 polymorphisms with breast cancer in an Australian Caucasian population

Breast cancer is a common disease in both developing and developed countries with early identification and treatment improving prognosis and survival. Heparan sulfate proteoglycans (HSPGs) are key components of the extracellular matrix (ECM) that mediate cell adhesion, motility, proliferation, invasion and cell signalling. Members of the syndecan family of HSPGs have been identified to be involved in breast cancer progression through their varied interactions with a number of growth factors, ligands and receptors. Specifically, high expression levels of syndecan-1 (SDC1) have been demonstrated in more invasive breast tumours while elevated syndecan-4 (SDC4) levels have been identified to correspond with improved prognosis. With genetic changes in the syndecans and their association with breast cancers plausible, we examined two single nucleotide polymorphisms in SDC1 (rs1131351) and SDC4 (rs67068737) within an Australian Caucasian breast cancer case/control population. No association was found with SDC4 and breast cancer in our population. However, a significant association between SDC1 and breast cancer was identified in both our case/control population and in a replication cohort. When both populations were combined for analysis, this association became more significant (genotype, p = 0.0003; allele, p = 0.0001). This data suggests an increased risk of developing breast cancer associated with the presence of the C allele of the SDC1 rs1131351 single nucleotide polymorphism (SNP) and may provide a marker toward early breast cancer detection.

Kallikrein 4 (hK4) and prostate-specific antigen (PSA) are associated with the loss of E-cadherin and an epithelial-mesenchymal transition (EMT)-like effect in prostate cancer cells

Prostate-specific antigen (PSA) and the related kallikrein family of serine proteases are current or emerging biomarkers for prostate cancer detection and progression. Kallikrein 4 (KLK4/hK4) is of particular interest, as KLK4 mRNA has been shown to be elevated in prostate cancer. In this study, we now show that the comparative expression of hK4 protein in prostate cancer tissues, compared with benign glands, is greater than that of PSA and kallikrein 2 (KLK2/hK2), suggesting that hK4 may play an important functional role in prostate cancer progression in addition to its biomarker potential. To examine the roles that hK4, as well as PSA and hK2, play in processes associated with progression, these kallikreins were separately transfected into the PC-3 prostate cancer cell line, and the consequence of their stable transfection was investigated. PC-3 cells expressing hK4 had a decreased growth rate, but no changes in cell proliferation were observed in the cells expressing PSA or hK2. hK4 and PSA, but not hK2, induced a 2.4-fold and 1.7-fold respective increase, in cellular migration, but not invasion, through Matrigel, a synthetic extracellular matrix. We hypothesised that this increase in motility displayed by the hK4 and PSA-expressing PC-3 cells may be related to the observed change in structure in these cells from a typical rounded epithelial-like cell to a spindle-shaped, more mesenchymal-like cell, with compromised adhesion to the culture surface. Thus, the expression of E-cadherin and vimentin, both associated with an epithelial-mesenchymal transition (EMT), was investigated. E-cadherin protein was lost and mRNA levels were significantly decreased in PC-3 cells expressing hK4 and PSA (10-fold and 7-fold respectively), suggesting transcriptional repression of E-cadherin, while the expression of vimentin was increased in these cells. The loss of E-cadherin and associated increase in vimentin are indicative of EMT and provides compelling evidence that hK4, in particular, and PSA have a functional role in the progression of prostate cancer through their promotion of tumour cell migration.

Quantifying uncertainty in parameter estimates for stochastic models of collective cell spreading using approximate Bayesian computation

Wound healing and tumour growth involve collective cell spreading, which is driven by individual motility and proliferation events within a population of cells. Mathematical models are often used to interpret experimental data and to estimate the parameters so that predictions can be made. Existing methods for parameter estimation typically assume that these parameters are constants and often ignore any uncertainty in the estimated values. We use approximate Bayesian computation (ABC) to estimate the cell diffusivity, D, and the cell proliferation rate, λ, from a discrete model of collective cell spreading, and we quantify the uncertainty associated with these estimates using Bayesian inference. We use a detailed experimental data set describing the collective cell spreading of 3T3 fibroblast cells. The ABC analysis is conducted for different combinations of initial cell densities and experimental times in two separate scenarios: (i) where collective cell spreading is driven by cell motility alone, and (ii) where collective cell spreading is driven by combined cell motility and cell proliferation. We find that D can be estimated precisely, with a small coefficient of variation (CV) of 2–6%. Our results indicate that D appears to depend on the experimental time, which is a feature that has been previously overlooked. Assuming that the values of D are the same in both experimental scenarios, we use the information about D from the first experimental scenario to obtain reasonably precise estimates of λ, with a CV between 4 and 12%. Our estimates of D and λ are consistent with previously reported values; however, our method is based on a straightforward measurement of the position of the leading edge whereas previous approaches have involved expensive cell counting techniques. Additional insights gained using a fully Bayesian approach justify the computational cost, especially since it allows us to accommodate information from different experiments in a principled way.

The tumour-promoting receptor tyrosine kinase, EphB4, regulates expression of Integrin-β8 in prostate cancer cells

Background The EphB4 receptor tyrosine kinase is overexpressed in many cancers including prostate cancer. The molecular mechanisms by which this ephrin receptor influences cancer progression are complex as there are tumor-promoting ligand-independent mechanisms in place as well as ligand-dependent tumor suppressive pathways. Methods We employed transient knockdown of EPHB4 in prostate cancer cells, coupled with gene microarray analysis, to identify genes that were regulated by EPHB4 and may represent linked tumor-promoting factors. We validated target genes using qRT-PCR and employed functional assays to determine their role in prostate cancer migration and invasion. Results We discovered that over 500 genes were deregulated upon EPHB4 siRNA knockdown, with integrin β8 (ITGB8) being the top hit (29-fold down-regulated compared to negative non-silencing siRNA). Gene ontology analysis found that the process of cell adhesion was highly deregulated and two other integrin genes, ITGA3 and ITGA10, were also differentially expressed. In parallel, we also discovered that over-expression of EPHB4 led to a concomitant increase in ITGB8 expression. In silico analysis of a prostate cancer progression microarray publically available in the Oncomine database showed that both EPHB4 and ITGB8 are highly expressed in prostatic intraepithelial neoplasia, the precursor to prostate cancer. Knockdown of ITGB8 in PC-3 and 22Rv1 prostate cancer cells in vitro resulted in significant reduction of cell migration and invasion. Conclusions These results reveal that EphB4 regulates integrin β8 expression and that integrin β8 plays a hitherto unrecognized role in the motility of prostate cancer cells and thus targeting integrin β8 may be a new treatment strategy for prostate cancer.