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.

Melanoma cell colony expansion parameters revealed by approximate Bayesian computation

In vitro studies and mathematical models are now being widely used to study the underlying mechanisms driving the expansion of cell colonies. This can improve our understanding of cancer formation and progression. Although much progress has been made in terms of developing and analysing mathematical models, far less progress has been made in terms of understanding how to estimate model parameters using experimental in vitro image-based data. To address this issue, a new approximate Bayesian computation (ABC) algorithm is proposed to estimate key parameters governing the expansion of melanoma cell (MM127) colonies, including cell diffusivity, D, cell proliferation rate, λ, and cell-to-cell adhesion, q, in two experimental scenarios, namely with and without a chemical treatment to suppress cell proliferation. Even when little prior biological knowledge about the parameters is assumed, all parameters are precisely inferred with a small posterior coefficient of variation, approximately 2–12%. The ABC analyses reveal that the posterior distributions of D and q depend on the experimental elapsed time, whereas the posterior distribution of λ does not. The posterior mean values of D and q are in the ranges 226–268 µm2h−1, 311–351 µm2h−1 and 0.23–0.39, 0.32–0.61 for the experimental periods of 0–24 h and 24–48 h, respectively. Furthermore, we found that the posterior distribution of q also depends on the initial cell density, whereas the posterior distributions of D and λ do not. The ABC approach also enables information from the two experiments to be combined, resulting in greater precision for all estimates of D and λ.

Nanofiber orientation and surface functionalization modulate human mesenchymal stem cell behavior in vitro

Electrospun nanofiber meshes have emerged as a new generation of scaffold membranes possessing a number of features suitable for tissue regeneration. One of these features is the flexibility to modify their structure and composition to orchestrate specific cellular responses. In this study, we investigated the effects of nanofiber orientation and surface functionalization on human mesenchymal stem cell (hMSC) migration and osteogenic differentiation. We used an in vitro model to examine hMSC migration into a cell-free zone on nanofiber meshes and mitomycin C treatment to assess the contribution of proliferation to the observed migration. Poly (ɛ-caprolactone) meshes with oriented topography were created by electrospinning aligned nanofibers on a rotating mandrel, while randomly oriented controls were collected on a stationary collector. Both aligned and random meshes were coated with a triple-helical, type I collagen-mimetic peptide, containing the glycine-phenylalanine-hydroxyproline-glycine-glutamate-arginine (GFOGER) motif. Our results indicate that nanofiber GFOGER peptide functionalization and orientation modulate cellular behavior, individually, and in combination. GFOGER significantly enhanced the migration, proliferation, and osteogenic differentiation of hMSCs on nanofiber meshes. Aligned nanofiber meshes displayed increased cell migration along the direction of fiber orientation compared to random meshes; however, fiber alignment did not influence osteogenic differentiation. Compared to each other, GFOGER coating resulted in a higher proliferation-driven cell migration, whereas fiber orientation appeared to generate a larger direct migratory effect. This study demonstrates that peptide surface modification and topographical cues associated with fiber alignment can be used to direct cellular behavior on nanofiber mesh scaffolds, which may be exploited for tissue regeneration.

Analysis of focal adhesion turnover: A quantitative live-cell imaging example

Recent advances in optical and fluorescent protein technology have rapidly raised expectations in cell biology, allowing quantitative insights into dynamic intracellular processes like never before. However, quantitative live-cell imaging comes with many challenges including how best to translate dynamic microscopy data into numerical outputs that can be used to make meaningful comparisons rather than relying on representative data sets. Here, we use analysis of focal adhesion turnover dynamics as a straightforward specific example on how to image, measure, and analyze intracellular protein dynamics, but we believe this outlines a thought process and can provide guidance on how to understand dynamic microcopy data of other intracellular structures.

Spectral analysis of pair-correlation bandwidth: Application to cell biology images

Images from cell biology experiments often indicate the presence of cell clustering, which can provide insight into the mechanisms driving the collective cell behaviour. Pair-correlation functions provide quantitative information about the presence, or absence, of clustering in a spatial distribution of cells. This is because the pair-correlation function describes the ratio of the abundance of pairs of cells, separated by a particular distance, relative to a randomly distributed reference population. Pair-correlation functions are often presented as a kernel density estimate where the frequency of pairs of objects are grouped using a particular bandwidth (or bin width), Δ>0. The choice of bandwidth has a dramatic impact: choosing Δ too large produces a pair-correlation function that contains insufficient information, whereas choosing Δ too small produces a pair-correlation signal dominated by fluctuations. Presently, there is little guidance available regarding how to make an objective choice of Δ. We present a new technique to choose Δ by analysing the power spectrum of the discrete Fourier transform of the pair-correlation function. Using synthetic simulation data, we confirm that our approach allows us to objectively choose Δ such that the appropriately binned pair-correlation function captures known features in uniform and clustered synthetic images. We also apply our technique to images from two different cell biology assays. The first assay corresponds to an approximately uniform distribution of cells, while the second assay involves a time series of images of a cell population which forms aggregates over time. The appropriately binned pair-correlation function allows us to make quantitative inferences about the average aggregate size, as well as quantifying how the average aggregate size changes with time.

The Neurofibromatosis 2 tumor suppressor merlin in cytoskeleton organization and cell cycle regulation

Neurofibromatosis 2 (NF2) is a dominantly inherited disorder, which predisposes to multiple tumours of the nervous system, typically schwannomas and meningiomas. Biallelic inactivation of the NF2 gene occurs both in sporadic and NF2-related schwannomas and in most meningiomas. The NF2 gene product merlin (or schwannomin) is structurally related to the ERM proteins, ezrin, radixin and moesin, which act as molecular linkers between the actin cytoskeleton and the plasma membrane. Merlin is a tumor suppressor that participates in cell cycle regulation. Merlin s phosphorylation status appears to be associated with its tumour suppressor activity, i.e. non-phosphorylated merlin functions as a tumour suppressor, whereas protein phosphorylation results in loss of functional activity. This thesis study was initiated to investigate merlin s role as a tumor suppressor and growth inhibitor. These studies show, that like many other tumor suppressors, also merlin is targeted to the nucleus at some stages of the cell cycle. Merlin s nuclear localization is regulated by cell cycle phase, contact inhibition and adhesion. In addition, a potential nuclear binding partner for merlin was identified, Human Enhancer of Invasion 10 (HEI10), a cyclin B interacting protein. Many tumor suppressors interact with microtubules and this thesis work shows that also merlin colocalizes with microtubules in mitotic structures. Merlin binds microtubules directly, and increases their polymerization in vitro and in vivo. In addition, primary mouse Schwann cells lacking merlin displays disturbed microtubule cytoskeleton. Fourth part of this thesis work began from the notion that PKA phosphorylates an unidentified site from the merlin N-terminus. Our studies show that serine 10 is a target for PKA and modulation of this residue regulates cytoskeletal organization, lamellipodia formation and cell migration. In summary, this thesis work shows that merlin s role is much more versatile than previously thought. It has a yet unidentified role in the nucleus and it participates in the regulation of both microtubules and the actin cytoskeleton. These studies have led to a better understanding of this enigmatic tumor suppressor, which eventually will aid in the design of specific drugs for the NF2 disease.

Bending the Rules of Cell Protrusions : Molecular Mechanisms and Biological Roles of Inverse-BAR Proteins in Cell Morphogenesis

Plasma membrane adopts myriad of different shapes to carry out essential cellular processes such as nutrient uptake, immunological defence mechanisms and cell migration. Therefore, the details how different plasma membrane structures are made and remodelled are of the upmost importance. Bending of plasma membrane into different shapes requires substantial amount of force, which can be provided by the actin cytoskeleton, however, the molecules that regulate the interplay between the actin cytoskeleton and plasma membrane have remained elusive. Recent findings have placed new types of effectors at sites of plasma membrane remodelling, including BAR proteins, which can directly bind and deform plasma membrane into different shapes. In addition to their membrane-bending abilities, BAR proteins also harbor protein domains that intimately link them to the actin cytoskeleton. The ancient BAR domain fold has evolved into at least three structurally and functionally different sub-groups: the BAR, F-BAR and I-BAR domains. This thesis work describes the discovery and functional characterization of the Inverse-BAR domains (I-BARs). Using synthetic model membranes, we have shown that I-BAR domains bind and deform membranes into tubular structures through a binding-surface composed of positively charged amino acids. Importantly, the membrane-binding surface of I-BAR domains displays an inverse geometry to that of the BAR and F-BAR domains, and these structural differences explain why I-BAR domains induce cell protrusions whereas BAR and most F-BAR domains induce cell invaginations. In addition, our results indicate that the binding of I-BAR domains to membranes can alter the spatial organization of phosphoinositides within membranes. Intriguingly, we also found that some I-BAR domains can insert helical motifs into the membrane bilayer, which has important consequences for their membrane binding/bending functions. In mammals there are five I-BAR domain containing proteins. Cell biological studies on ABBA revealed that it is highly expressed in radial glial cells during the development of the central nervous system and plays an important role in the extension process of radial glia-like C6R cells by regulating lamellipodial dynamics through its I-BAR domain. To reveal the role of these proteins in the context of animals, we analyzed MIM knockout mice and found that MIM is required for proper renal functions in adult mice. MIM deficient mice displayed a severe urine concentration defect due to defective intercellular junctions of the kidney epithelia. Consistently, MIM localized to adherens junctions in cultured kidney epithelial cells, where it promoted actin assembly through its I-BAR andWH2 domains. In summary, this thesis describes the mechanism how I-BAR proteins deform membranes and provides information about the biological role of these proteins, which to our knowledge are the first proteins that have been shown to directly deform plasma membrane to make cell protrusions.

N-syndecan and HB-GAM in neural migration and differentiation : Modulation of growth factor activity in brain

The juvenile sea squirt wanders through the sea searching for a suitable rock or hunk of coral to cling to and make its home for life. For this task it has a rudimentary nervous system. When it finds its spot and takes root, it doesn't need its brain any more so it eats it. It's rather like getting tenure. Daniel C. Dennett (from Consciousness Explained, 1991) The little sea squirt needs its brain for a task that is very simple and short. When the task is completed, the sea squirt starts a new life in a vegetative state, after having a nourishing meal. The little brain is more tightly structured than our massive primate brains. The number of neurons is exact, no leeway in neural proliferation is tolerated. Each neuroblast migrates exactly to the correct position, and only a certain number of connections with the right companions is allowed. In comparison, growth of a mammalian brain is a merry mess. The reason is obvious: Squirt brain needs to perform only a few, predictable functions, before becoming waste. The more mobile and complex mammals engage their brains in tasks requiring quick adaptation and plasticity in a constantly changing environment. Although the regulation of nervous system development varies between species, many regulatory elements remain the same. For example, all multicellular animals possess a collection of proteoglycans (PG); proteins with attached, complex sugar chains called glycosaminoglycans (GAG). In development, PGs participate in the organization of the animal body, like in the construction of parts of the nervous system. The PGs capture water with their GAG chains, forming a biochemically active gel at the surface of the cell, and in the extracellular matrix (ECM). In the nervous system, this gel traps inside it different molecules: growth factors and ECM-associated proteins. They regulate the proliferation of neural stem cells (NSC), guide the migration of neurons, and coordinate the formation of neuronal connections. In this work I have followed the role of two molecules contributing to the complexity of mammalian brain development. N-syndecan is a transmembrane heparan sulfate proteoglycan (HSPG) with cell signaling functions. Heparin-binding growth-associated molecule (HB-GAM) is an ECM-associated protein with high expression in the perinatal nervous system, and high affinity to HS and heparin. N-syndecan is a receptor for several growth factors and for HB-GAM. HB-GAM induces specific signaling via N-syndecan, activating c-Src, calcium/calmodulin-dependent serine protein kinase (CASK) and cortactin. By studying the gene knockouts of HB-GAM and N-syndecan in mice, I have found that HB-GAM and N-syndecan are involved as a receptor-ligand-pair in neural migration and differentiation. HB-GAM competes with the growth factors fibriblast growth factor (FGF)-2 and heparin-binding epidermal growth factor (HB-EGF) in HS-binding, causing NSCs to stop proliferation and to differentiate, and affects HB-EGF-induced EGF receptor (EGFR) signaling in neural cells during migration. N-syndecan signaling affects the motility of young neurons, by boosting EGFR-mediated cell migration. In addition, these two receptors form a complex at the surface of the neurons, probably creating a motility-regulating structure.

Missing-In-Metastasis (MIM) regulates cell morphology by promoting plasma membrane and actin cytoskeleton dynamics

The cells of multicellular organisms have differentiated to carry out specific functions that are often accompanied by distinct cell morphology. The actin cytoskeleton is one of the key regulators of cell shape subsequently controlling multiple cellular events including cell migration, cell division, endo- and exocytosis. A large set of actin regulating proteins has evolved to achieve and tightly coordinate this wide range of functions. Some actin regulator proteins have so-called house keeping roles and are essential for all eukaryotic cells, but some have evolved to meet the requirements of more specialized cell-types found in higher organisms enabling complex functions of differentiated organs, such as liver, kidney and brain. Often processes mediated by the actin cytoskeleton, like formation of cellular protrusions during cell migration, are intimately linked to plasma membrane remodeling. Thus, a close cooperation between these two cellular compartments is necessary, yet not much is known about the underlying molecular mechanisms. This study focused on a vertebrate-specific protein called missing-in-metastasis (MIM), which was originally characterized as a metastasis suppressor of bladder cancer. We demonstrated that MIM regulates the dynamics of actin cytoskeleton via its WH2 domain, and is expressed in a cell-type specific manner. Interestingly, further examination showed that the IM-domain of MIM displays a novel membrane tubulation activity, which induces formation of filopodia in cells. Following studies demonstrated that this membrane deformation activity is crucial for cell protrusions driven by MIM. In mammals, there are five members of IM-domain protein family. Functions and expression patterns of these family members have remained poorly characterized. To understand the physiological functions of MIM, we generated MIM knockout mice. MIM-deficient mice display no apparent developmental defects, but instead suffer from progressive renal disease and increased susceptibility to tumors. This indicates that MIM plays a role in the maintenance of specific physiological functions associated with distinct cell morphologies. Taken together, these studies implicate MIM both in the regulation of the actin cytoskeleton and the plasma membrane. Our results thus suggest that members of MIM/IRSp53 protein family coordinate the actin cytoskeleton:plasma membrane interface to control cell and tissue morphogenesis in multicellular organisms.

Reproducibility of scratch assays is affected by the initial degree of confluence: Experiments, modelling and model selection

Scratch assays are difficult to reproduce. Here we identify a previously overlooked source of variability which could partially explain this difficulty. We analyse a suite of scratch assays in which we vary the initial degree of confluence (initial cell density). Our results indicate that the rate of re-colonisation is very sensitive to the initial density. To quantify the relative roles of cell migration and proliferation, we calibrate the solution of the Fisher–Kolmogorov model to cell density profiles to provide estimates of the cell diffusivity, D, and the cell proliferation rate, λ. This procedure indicates that the estimates of D and λ are very sensitive to the initial density. This dependence suggests that the Fisher–Kolmogorov model does not accurately represent the details of the collective cell spreading process, since this model assumes that D and λ are constants that ought to be independent of the initial density. Since higher initial cell density leads to enhanced spreading, we also calibrate the solution of the Porous–Fisher model to the data as this model assumes that the cell flux is an increasing function of the cell density. Estimates of D and λ associated with the Porous–Fisher model are less sensitive to the initial density, suggesting that the Porous–Fisher model provides a better description of the experiments.

Acto-myosin contractility rotates the cell nucleus

The nucleus of the eukaryotic cell functions amidst active cytoskeletal �laments, but its response to the stresses carried by these �laments is largely unexplored. We report here the results of studies of the translational and rotational dynamics of the nuclei of single �broblast cells, with the e�ects of cell migration suppressed by plating onto �bronectin-coated micro-fabricated patterns. Patterns of the same area but di�erent shapes and/or aspect ratio were used to study the e�ect of cell geometry on the dynamics. On circles, squares and equilateral triangles, the nucleus undergoes persistent rotational motion, while on high-aspect-ratio rectangles of the same area it moves only back and forth. The circle and the triangle showed respectively the largest and the smallest angular speed. We show that our observations can be understood through a hydrodynamic approach in which the nucleus is treated as a highly viscous inclusion residing in a less viscous uid of orientable �laments endowed with active stresses. Lowering actin contractility selectively by introducing blebbistatin at low concentrations drastically reduced the speed and persistence time of the angular motion of the nucleus. Time-lapse imaging of actin revealed a correlated hydrodynamic ow around the nucleus, with pro�le and magnitude consistent with the results of our theoretical approach. Coherent intracellular ows and consequent nuclear rotation thus appear to be a generic property that cells must balance by speci�c mechanisms in order to maintain nuclear homeostasis

Actomyosin contractility rotates the cell nucleus

The cell nucleus functions amidst active cytoskeletal filaments, but its response to their contractile stresses is largely unexplored. We study the dynamics of the nuclei of single fibroblasts, with cell migration suppressed by plating onto micro-fabricated patterns. We find the nucleus undergoes noisy but coherent rotational motion. We account for this observation through a hydrodynamic approach, treating the nucleus as a highly viscous inclusion residing in a less viscous fluid of orientable filaments endowed with active stresses. Lowering actin contractility selectively by introducing blebbistatin at low concentrations drastically reduced the speed and coherence of the angular motion of the nucleus. Time-lapse imaging of actin revealed a correlated hydrodynamic flow around the nucleus, with profile and magnitude consistent with the results of our theoretical approach. Coherent intracellular flows and consequent nuclear rotation thus appear to be an intrinsic property of cells.

Inhibition of p300 lysine acetyltransferase activity by luteolin reduces tumor growth in head and neck squamous cell carcinoma (HNSCC) xenograft mouse model

Chromatin acetylation is attributed with distinct functional relevance with respect to gene expression in normal and diseased conditions thereby leading to a topical interest in the concept of epigenetic modulators and therapy. We report here the identification and characterization of the acetylation inhibitory potential of an important dietary flavonoid, luteolin. Luteolin was found to inhibit p300 acetyltransferase with competitive binding to the acetyl CoA binding site. Luteolin treatment in a xenografted tumor model of head and neck squamous cell carcinoma (HNSCC), led to a dramatic reduction in tumor growth within 4 weeks corresponding to a decrease in histone acetylation. Cells treated with luteolin exhibit cell cycle arrest and decreased cell migration. Luteolin treatment led to an alteration in gene expression and miRNA profile including up-regulation of p53 induced miR-195/215, let7C; potentially translating into a tumor suppressor function. It also led to down regulation of oncomiRNAs such as miR-135a, thereby reflecting global changes in the microRNA network. Furthermore, a direct correlation between the inhibition of histone acetylation and gene expression was established using chromatin immunoprecipitation on promoters of differentially expressed genes. A network of dysregulated genes and miRNAs was mapped along with the gene ontology categories, and the effects of luteolin were observed to be potentially at multiple levels: at the level of gene expression, miRNA expression and miRNA processing.

Advancing Edge Speeds of Epithelial Monolayers Depend on Their Initial Confining Geometry

Collective cell migrations are essential in several physiological processes and are driven by both chemical and mechanical cues. The roles of substrate stiffness and confinement on collective migrations have been investigated in recent years, however few studies have addressed how geometric shapes influence collective cell migrations. Here, we address the hypothesis that the relative position of a cell within the confinement influences its motility. Monolayers of two types of epithelial cells-MCF7, a breast epithelial cancer cell line, and MDCK, a control epithelial cell line-were confined within circular, square, and cross-shaped stencils and their migration velocities were quantified upon release of the constraint using particle image velocimetry. The choice of stencil geometry allowed us to investigate individual cell motility within convex, straight and concave boundaries. Cells located in sharp, convex boundaries migrated at slower rates than those in concave or straight edges in both cell types. The overall cluster migration occurred in three phases: an initial linear increase with time, followed by a plateau region and a subsequent decrease in cluster speeds. An acto-myosin contractile ring, present in the MDCK but absent in MCF7 monolayer, was a prominent feature in the emergence of leader cells from the MDCK clusters which occurred every similar to 125 mu m from the vertex of the cross. Further, coordinated cell movements displayed vorticity patterns in MDCK which were absent in MCF7 clusters. We also used cytoskeletal inhibitors to show the importance of acto-myosin bounding cables in collective migrations through translation of local movements to create long range coordinated movements and the creation of leader cells within ensembles. To our knowledge, this is the first demonstration of how bounding shapes influence long-term migratory behaviours of epithelial cell monolayers. These results are important for tissue engineering and may also enhance our understanding of cell movements during developmental patterning and cancer metastasis.