The impact of "search precision" in an agent-based tumor model

Several cell surface receptors are overexpressed in malignant brain tumors and reportedly involved in tumor progression and invasion. It is unclear, however, whether such an improvement of cellular signal reception leads to a monotonic increase in the tumor system's average velocity during invasion or whether there is a maximum threshold beyond which the average velocity starts to decelerate. To gain more insight into the systemic effects of such cellular search precision modulations, this study proposes a two-dimensional agent-based model in which the spatio-temporal expansion of malignant brain tumor cells is guided by environmental heterogeneities in mechanical confinement, toxic metabolites and nutrient sources. Here, the spatial field of action is represented by an adaptive grid lattice, which corresponds to the experimental finding that tumor cells are more likely to follow each other along preformed pathways. Another prominent feature is the dual threshold concept for both nutrient level and toxicity, which determine whether cells proliferate, migrate, remain quiescent or die in the next period. The numerical results from varying the key parameters encoding the capability of tumor cells to invade and their ability to proliferate indicate an emergent behavior. Specifically, increasing invasiveness not only leads to an increase in maximum expansion velocity, but also requires a more precise spatial search process, corresponding to an improved cell signal reception, in order to obtain maximum velocity. To increase cellular invasiveness beyond the maximum that can be achieved by exclusively tuning the motility parameter, it requires an additional reduction in the cells' proliferation rate and prompts an even more biased search process. Most interestingly, however, a prominent phase transition suggests that tumor cells do not employ a 100 percent search precision to attain maximum spatial velocity. These findings argue for a selection advantage conferred by limited randomness in processing spatial search and indicate that our computational platform may prove valuable in investigating emergent, multicellular tumor patterns caused by alterations on the molecular level.     

PATHWAYS OF CELLULAR MORPHOGENESIS : A Diversity in Nitella

Evidence is presented to show that a given change in cell form or size may generally be brought about by a variety of patterns of local surface distortion and expansion. Structural and chemical features of the cell which are important in morphogenesis may thus be expected to relate not to form per se but to the kinetics of surface behavior which establish form. These kinetics evaluate both the rate at which local regions of cell surface expand and the directed character (anisotropy) of this expansion. These variables have been studied in model systems and, through marking experiments, in growing cells of various shapes in Phycomyces, Clypeaster, and particularly Nitella. In the latter plant, prominent "giant internodes" display a well defined longitudinal anisotropic expansion devoid of sizeable gradients in expansion rate. These cells have their origin, however, in apical cells which have a pronounced gradient in area expansion rate (maximal at the tip). The great part of the expansion in the apical cell is apparently isotropic (equal in all directions), but the basal region often shows predominant expansion laterally. This transverse stretching in the apical cell could align cell wall texture and possibly fibrous cytoplasmic constituents, such as microtubules, into configurations significant in later morphogenetic stages, including the elongation of the internodes.     

Reprogramming of human cancer cells to pluripotency for models of cancer progression

The ability to study live cells as they progress through the stages of cancer provides the opportunity to discover dynamic networks underlying pathology, markers of early stages, and ways to assess therapeutics. Genetically engineered animal models of cancer, where it is possible to study the consequences of temporal‐specific induction of oncogenes or deletion of tumor suppressors, have yielded major insights into cancer progression. Yet differences exist between animal and human cancers, such as in markers of progression and response to therapeutics. Thus, there is a need for human cell models of cancer progression. Most human cell models of cancer are based on tumor cell lines and xenografts of primary tumor cells that resemble the advanced tumor state, from which the cells were derived, and thus do not recapitulate disease progression. Yet a subset of cancer types have been reprogrammed to pluripotency or near‐pluripotency by blastocyst injection, by somatic cell nuclear transfer and by induced pluripotent stem cell (iPS) technology. The reprogrammed cancer cells show that pluripotency can transiently dominate over the cancer phenotype. Diverse studies show that reprogrammed cancer cells can, in some cases, exhibit early‐stage phenotypes reflective of only partial expression of the cancer genome. In one case, reprogrammed human pancreatic cancer cells have been shown to recapitulate stages of cancer progression, from early to late stages, thus providing a model for studying pancreatic cancer development in human cells where previously such could only be discerned from mouse models. We discuss these findings, the challenges in developing such models and their current limitations, and ways that iPS reprogramming may be enhanced to develop human cell models of cancer progression.     

Local interactions shape plant cells

Plant cell expansion is usually attributed to the considerable osmotic pressure that develops within and impinges upon the cell boundary. Whereas turgor containment within expandable walls explains global expansion, the scalar nature of turgor does not directly suggest a mechanism for achieving the localized, differential growth that is responsible for the diversity of plant-cell forms. The key to achieving local growth in plant cells appears to lie not in harnessing turgor but in using it to identify weak regions in the cell boundary and thus creating discrete intracellular domains for targeting the growth machinery. Membrane-interacting phospholipases, Rho-like proteins and their interactors, an actin-modulating ARP2/3 complex with its upstream regulators, and actin-microtubule interactions play important roles in the intracellular cooperation to shape plant cells.     

A Mechanistic Collective Cell Model for Epithelial Colony Growth and Contact Inhibition

We present a mechanistic hybrid continuum-discrete model to simulate the dynamics of epithelial cell colonies. Collective cell dynamics are modeled using continuum equations that capture plastic, viscoelastic, and elastic deformations in the clusters while providing single-cell resolution. The continuum equations can be viewed as a coarse-grained version of previously developed discrete models that treat epithelial clusters as a two-dimensional network of vertices or stochastic interacting particles and follow the framework of dynamic density functional theory appropriately modified to account for cell size and shape variability. The discrete component of the model implements cell division and thus influences cell size and shape that couple to the continuum component. The model is validated against recent in vitro studies of epithelial cell colonies using Madin-Darby canine kidney cells. In good agreement with experiments, we find that mechanical interactions and constraints on the local expansion of cell size cause inhibition of cell motion and reductive cell division. This leads to successively smaller cells and a transition from exponential to quadratic growth of the colony that is associated with a constant-thickness rim of growing cells at the cluster edge, as well as the emergence of short-range ordering and solid-like behavior. A detailed analysis of the model reveals a scale invariance of the growth and provides insight into the generation of stresses and their influence on the dynamics of the colonies. Compared to previous models, our approach has several advantages: it is independent of dimension, it can be parameterized using classical elastic properties (Poisson’s ratio and Young’s modulus), and it can easily be extended to incorporate multiple cell types and general substrate geometries.     

Activation of a novel proto-oncogene, Frat1, contributes to progression of mouse T-cell lymphomas.

Acceleration of lymphomagenesis in oncogene-bearing transgenic mice by slow-transforming retroviruses has proven a valuable tool in identifying cooperating oncogenes. We have modified this protocol to search for genes that can collaborate effectively with the transgene in later stages of tumor development. Propagation of tumors induced by Moloney murine leukemia virus (M-MuLV) in E mu-Pim1 or H2-K-myc transgenic mice by transplantation to syngeneic hosts permitted proviral tagging of 'progression' genes. Molecular cloning of common proviral insertion sites that were detected preferentially in transplanted tumors led to the identification of a novel gene, designated Frat1. The initial selection for integrations near Frat1 occurs in primary tumor cells that have already acquired proviruses in other common insertion sites, yielding primary lymphomas that contain only a minor fraction of tumor cells with an activated Frat1 allele. Transplantation of such primary lymphomas allows for a further expansion of tumor cell clones carrying a proviral insertion near Frat1, resulting in detectable Frat1 rearrangements in 17% of the transplanted E mu-Pim1 tumors and 30% of the transplanted H2-K-myc tumors, respectively. We have cloned and sequenced both the mouse Frat1 gene and its human counterpart. The proteins encoded by Frat1 and FRAT1 are highly homologous and their functions are thus far unknown. Tumor cell lines with high expression of Myc and Pim1 acquired an additional selective advantage in vivo upon infection with a Frat1-IRES-lacZ retrovirus, thus underscoring the role of Frat1 in tumor progression, and the ability of Frat1 to collaborate with Pim1 and Myc in lymphomagenesis.     

Cellular search migrations in normal development and carcinogenesis

This review describes the large group of morphogenetic processes designated as search migrations. Search migrations typically include two stages: i) search, when a group of cells or of the cytoplasmic processes migrate over the cell-free spaces, and ii) choice, the stage when migrating cells reach specific loci where they stop and undergo specific differentiations induced by local factors such as cell-cell contacts and humoral agents. Migrating cells that do not meet their targets usually undergo apoptosis. Numerous examples of search migrations range from gastrulation to formation of axon-muscle connections. Critical stages of carcinogenesis such as acquisition of cell ability for invasion may be regarded as the genetic aberration of normal search migration: cancer cells perform an endless search but cannot make final choice.     

Finite element simulations of the active stress in the imaginal disc of the Drosophila Melanogaster

During the larval stages of development, the imaginal disc of Drosphila Melanogaster is composed by a monolayer of epithelial cells, which undergo a strain actively produced by the cells themselves. The well-organized collective contraction produces a stress field that seemingly has a double morphogenetic role: it orchestrates the cellular organization towards the macroscopic shape emergence while simultaneously providing a local information on the organ size. Here we perform numerical simulations of such a mechanical control on morphogenesis at a continuum level, using a three-dimensional finite model that accounts for the active cell contraction. The numerical model is able to reproduce the (few) known qualitative characteristics of the tensional patterns within the imaginal disc of the fruit fly. The computed stress components slightly deviate from planarity, thus confirming the previous theoretical assumptions of a nonlinear elastic analytical model, and enforcing the hypothesis that the spatial variation of the mechanical stress may act as a size regulating signal that locally scales with the global dimension of the domain.     

Activation of CD95L fusion protein prodrugs by tumor-associated proteases

To achieve tumor cell-restricted activation of CD95, we developed a CD95L fusion protein format, in which CD95L activity is only unmasked upon antibody-mediated binding to tumor cells and subsequent processing by tumor-associated proteases, such as matrix metalloproteases (MMPs) and urokinase plasminogen activator (uPA). On target-negative, but MMP- and uPA-expressing HT1080 tumor cells, the CD95L prodrugs were virtually inactive. On target antigen-expressing HT1080 cells, however, the CD95L prodrugs showed an apoptotic activity comparable to soluble CD95L artificially activated by crosslinking. CD95 activation by the CD95L prodrugs was preceded by prodrug processing. Apoptosis was blocked by inhibitors of MMPs or uPA and by neutralizing antibodies recognizing the targeted cell surface antigen or the CD95L moiety of the prodrugs. In a xenotransplantation tumor model, local application of the prodrug reduced the growth of target antigen-expressing, but not antigen-negative tumor cells, verifying targeted CD95L prodrug activation in vivo.     

Nanoscale and geometric influences on the microtubule cytoskeleton in plants: thinking inside and outside the box

The dynamic microtubule (MT) cytoskeleton found in the cell cortex of plants drives cell expansion via cell wall modifications. In the last decade, live cell imaging studies employing green fluorescent protein have helped unravel the mechanisms behind how cells arrange cortical MTs into complex arrays and shape cell expansion. In this review, we explore the reverse scenario: how cell geometry and organelles influence and constrain the organization and behavior of cortical MTs. This newly emerging principle explains how cells perceive local nanoscale structural input from MT-organizing centers, such as the nucleus, endomembranes, and cell edges, and translate this into global cell-wide order via MT self-organization. Studies primarily using the model plant Arabidopsis thaliana and tobacco BY-2 suspension cultures have broadened our understanding of how cells form not only elegant parallel arrays but also more complex MT configurations, including the prominent MT bundles found in preprophase bands, leaf epidermal cells, and developing xylem.     

Nanoscale and geometric influences on the microtubule cytoskeleton in plants: thinking inside and outside the box

The dynamic microtubule (MT) cytoskeleton found in the cell cortex of plants drives cell expansion via cell wall modifications. In the last decade, live cell imaging studies employing green fluorescent protein have helped unravel the mechanisms behind how cells arrange cortical MTs into complex arrays and shape cell expansion. In this review, we explore the reverse scenario: how cell geometry and organelles influence and constrain the organization and behavior of cortical MTs. This newly emerging principle explains how cells perceive local nanoscale structural input from MT-organizing centers, such as the nucleus, endomembranes, and cell edges, and translate this into global cell-wide order via MT self-organization. Studies primarily using the model plant Arabidopsis thaliana and tobacco BY-2 suspension cultures have broadened our understanding of how cells form not only elegant parallel arrays but also more complex MT configurations, including the prominent MT bundles found in preprophase bands, leaf epidermal cells, and developing xylem.     

Criticality in intracellular calcium signaling in cardiac myocytes

Calcium (Ca) is a ubiquitous second messenger that regulates many biological functions. The elementary events of local Ca signaling are Ca sparks, which occur randomly in time and space, and integrate to produce global signaling events such as intra- and intercellular Ca waves and whole-cell Ca oscillations. Despite extensive experimental characterization in many systems, the transition from local random to global synchronous events is still poorly understood. Here we show that criticality, a ubiquitous dynamical phenomenon in nature, is responsible for the transition from local to global Ca signaling. We demonstrate this first in a computational model of Ca signaling in a cardiac myocyte and then experimentally in mouse ventricular myocytes, complemented by a theoretical agent-based model to delineate the underlying dynamics. We show that the interaction between the Ca release units via Ca-induced Ca release causes self-organization of Ca spark clusters. When the coupling between Ca release units is weak, the cluster-size distribution is exponential. As the interactions become strong, the cluster-size distribution changes to a power-law distribution, which is characteristic of criticality in thermodynamic and complex nonlinear systems, and facilitates the formation and propagation of Ca waves and whole-cell Ca oscillations. Our findings illustrate how criticality is harnessed by a biological cell to regulate Ca signaling via self-organization of random subcellular events into cellular-scale oscillations, and provide a general theoretical framework for the transition from local Ca signaling to global Ca signaling in biological cells.     

Induction and large-scale expansion of CD8+ tumor specific cytotoxic T lymphocytes from peripheral blood lymphocytes by in vitro stimulation with CD80-transfected autologous melanoma cells.

Human melanoma cell lines may induce a specific T cell response against tumor cells in vitro. However, after repeated restimulation with autologous tumor cells, expansion of CTL is limited and often apoptosis of the T cells occurs. In order to improve conditions inducing primary T cell responses and thus allowing further expansion of tumor specific T cells for an adoptive transfer, we transfected human melanoma cells with the B7.1 gene (CD80), known to be a potent costimulatory molecule for T cell activation. CD80 expression on melanoma cells resulted in improved primary T cell activation, especially of CD8+ T cells. Furthermore, restimulation with CD80+ tumor cells gave rise to long term proliferating CD8+ T cell lines demonstrating an 100-fold expansion of T cells compared to the 20-30-fold increased numbers obtained with the controls (parental tumor cells +/- anti-CD28). T cells stimulated with CD80+ melanoma cells were found to display a MHC class I-restricted cytotoxic activity against the autologous tumor cells. In conclusion, these studies demonstrate the requirement of costimulation in generating large numbers of tumor specific T cells in vitro that may be used for an adoptive transfer in tumor immunotherapy.     

Nanoscale imaging and quantification of local proteolytic activity

Proteolytic cleavage of extracellular matrix (ECM) is a critical feature of tumor cell invasion, and affects cancer cell growth, differentiation, apoptosis, and migration. Malignant cells secrete most proteases as inactive proenzymes that undergo proteolytic cleavage for activation, and proteolytic activity is elevated in close proximity to these cells. Therefore, local activity rather than protease concentration determines ECM proteolysis. Precise quantification of local proteolytic activity, functional investigation, and high resolution imaging of morphological ECM alterations have proven difficult. In this study, we present a novel approach for measuring proteolytic activity in the microenvironment of cells by using atomic force microscopy (AFM). Amelanotic melanoma cells (A7-clone) were seeded on fluorescent gelatin or collagen-IV coatings. Proteolysis reduced fluorescence of these coatings. Fluorescence microscopy (FM) in combination with AFM was used to maneuver the AFM-tip to tumor cell induced proteolytic spots. AFM enabled nanoscale volume measurement, three-dimensional reconstruction of single proteins and demonstrated that ECM cleavage is restricted to the proteolytic microenvironment of cancer cells. This method detected significant decreases in molecular weight of protein clusters (-76.6%), matrix volume (-46.6%), and height (-38.1%) between intact and proteolyzed gelatin. Similar parameter changes were demonstrated without FM, by AFM-scanning gelatin in close proximity to invasive cells. Furthermore, AFM depicted significantly stronger local degradation of gelatin than collagen-IV by A7-cells. Taken together, AFM allows specific quantification and imaging of local proteolytic processes at a nanometer level, thus providing a unique method for the functional evaluation of invasiveness and metastatic potential of tumor cells in small scale samples.     

A large number of T lymphocytes recognize Moloney-murine leukemia virus-induced antigens, but a few mediate long-lasting tumor immunosurveillance

The CD8(+) T cell response to Moloney-murine leukemia virus (M-MuLV)-induced Ags is almost entirely dominated by the exclusive expansion of lymphocytes that use preferential TCRVbeta chain rearrangements. In mice lacking T cells expressing these TCRVbeta, we demonstrate that alternative TCRVbeta can substitute for the lack of the dominant TCRVbeta in the H-2-restricted M-MuLV Ag recognition. We show that, at least for the H-2(b)-restricted response, the shift of TCR usage is not related to a variation of the immunodominant M-MuLV epitope recognition. After virus immunization, all the potentially M-MuLV-reactive lymphocytes are primed, but only the deletion of dominant Vbeta rescues the alternative Vbeta response. The mechanism of clonal T cell "immunodomination" that guides the preferential Vbeta expansion is likely the result of a proliferative advantage of T cells expressing dominant Vbeta, due to differences in TCR affinity and/or cosignal requirements. In this regard, a CD8 involvement is strictly required for the virus-specific cytotoxic activity of CTL expressing alternative, but not dominant, Vbeta gene rearrangements. The ability of T cells expressing alternative TCRVbeta rearrangements to mediate tumor protection was evaluated by a challenge with M-MuLV tumor cells. Although T cells expressing alternative Vbeta chains were activated and expanded, they were not able to control tumor growth in a long-lasting manner due to their incapacity of conversion and accumulation in the T central memory pool.     

Actin polymerization and intracellular solvent flow in cell surface blebbing

The cortical actin gel of eukaryotic cells is postulated to control cell surface activity. One type of protrusion that may offer clues to this regulation are the spherical aneurysms of the surface membrane known as blebs. Blebs occur normally in cells during spreading and alternate with other protrusions, such as ruffles, suggesting similar protrusive machinery is involved. We recently reported that human melanoma cell lines deficient in the actin filament cross-linking protein, ABP-280, show prolonged blebbing, thus allowing close study of blebs and their dynamics. Blebs expand at different rates of volume increase that directly predict the final size achieved by each bleb. These rates decrease as the F-actin concentration of the cells increase over time after plating on a surface, but do so at lower concentrations in ABP-280 expressing cells. Fluorescently labeled actin and phalloidin injections of blebbing cells indicate that a polymerized actin structure is not present initially, but appears later and is responsible for stopping further bleb expansion. Therefore, it is postulated that blebs occur when the fluid-driven expansion of the cell membrane is sufficiently rapid to initially outpace the local rate of actin polymerization. In this model, the rate of intracellular solvent flow driving this expansion decreases as cortical gelation is achieved, whether by factors such as ABP-280, or by concentrated actin polymers alone, thereby leading to decreased size and occurrence of blebs. Since the forces driving bleb extension would always be present in a cell, this process may influence other cell protrusions as well.     

Effect of platelet-activating factor receptor expression on CHO cell motility

Tumor cell migration may favor local mass expansion and metastasis dissemination. Several tumors were found to express the receptor for platelet-activating factor (PAF), a potent mediator of leukocyte chemotaxis and endothelial cell migration. However, its functional role on tumor cells is largely unexplored. In the present study, we evaluated the motogenic effect of PAF on Chinese hamster ovarian (CHO) cancer cells transfected with the human PAF-receptor cDNA (CHO PAF-R). By using time-lapse recording, we detected a rapid motogenic response to PAF stimulation on CHO PAF-R, whereas no effect was evident on vector-only transfected cells. Such an effect was observed on scattered cell motility, on cells seeded on a fibronectin- or collagen-coated surface, and on migration of confluent monolayer cells. Cell speed increased at 1 h and was maximal 6-8 h after PAF stimulation on CHO PAF-R. Concomitantly, PAF induced marked changes in cytoskeleton actin distribution with cell contraction, assembling of stress fibers, and polar foci of adhesion. In conclusion, the present study demonstrates that PAF is a potent inducer of tumor cell motility, thus suggesting a role for this mediator in tumor growth and dissemination.     

fMRI evidence for a cortical hierarchy of pitch pattern processing

Pitch patterns, such as melodies, consist of two levels of structure: a global level, comprising the pattern of ups and downs, or contour; and a local level, comprising the precise intervals that make up this contour. An influential neuropsychological model suggests that these two levels of processing are hierarchically linked, with processing of the global structure occurring within the right hemisphere in advance of local processing within the left. However, the predictions of this model and its anatomical basis have not been tested in neurologically normal individuals. The present study used fMRI and required participants to listen to consecutive pitch sequences while performing a same/different one-back task. Sequences, when different, either preserved (local) or violated (global) the contour of the sequence preceding them. When the activations for the local and global conditions were contrasted directly, additional activation was seen for local processing in right planum temporale and posterior superior temporal sulcus (pSTS). The presence of additional activation for local over global processing supports the hierarchical view that the global structure of a pitch sequence acts as a "framework" on which the local detail is subsequently hung. However, the lateralisation of activation seen in the present study, with global processing occurring in left pSTS and local processing occurring bilaterally, differed from that predicted by the neuroanatomical model. A re-examination of the individual lesion data on which the neuroanatomical model is based revealed that the lesion data equally well support the laterality scheme suggested by our data. While the present study supports the hierarchical view of local and global processing, there is an evident need for further research, both in patients and neurologically normal individuals, before an understanding of the functional lateralisation of local and global processing can be considered established.     

Multi-scale agent-based brain cancer modeling and prediction of TKI treatment response: Incorporating EGFR signaling pathway and angiogenesis

ABSTRACT: BACKGROUND: The epidermal growth factor receptor (EGFR) signaling pathway and angiogenesis in brain cancer act as an engine for tumor initiation, expansion and response to therapy. Since the existing literature does not have any models that investigate the impact of both angiogenesis and molecular signaling pathways on treatment, we propose a novel multi-scale, agent-based computational model that includes both angiogenesis and EGFR modules to study the response of brain cancer under tyrosine kinase inhibitors (TKIs) treatment. RESULTS: The novel angiogenesis module integrated into the agent-based tumor model is based on a set of reaction--diffusion equations that describe the spatio-temporal evolution of the distributions of micro-environmental factors such as glucose, oxygen, TGFalpha, VEGF and fibronectin. These molecular species regulate tumor growth during angiogenesis. Each tumor cell is equipped with an EGFR signaling pathway linked to a cell-cycle pathway to determine its phenotype. EGFR TKIs are delivered through the blood vessels of tumor microvasculature and the response to treatment is studied. CONCLUSIONS: Our simulations demonstrated that entire tumor growth profile is a collective behaviour of cells regulated by the EGFR signaling pathway and the cell cycle. We also found that angiogenesis has a dual effect under TKI treatment: on one hand, through neo-vasculature TKIs are delivered to decrease tumor invasion; on the other hand, the neo-vasculature can transport glucose and oxygen to tumor cells to maintain their metabolism, which results in an increase of cell survival rate in the late simulation stages.     

The life history of 21 breast cancers

Cancer evolves dynamically as clonal expansions supersede one another driven by shifting selective pressures, mutational processes, and disrupted cancer genes. These processes mark the genome, such that a cancer's life history is encrypted in the somatic mutations present. We developed algorithms to decipher this narrative and applied them to 21 breast cancers. Mutational processes evolve across a cancer's lifespan, with many emerging late but contributing extensive genetic variation. Subclonal diversification is prominent, and most mutations are found in just a fraction of tumor cells. Every tumor has a dominant subclonal lineage, representing more than 50% of tumor cells. Minimal expansion of these subclones occurs until many hundreds to thousands of mutations have accumulated, implying the existence of long-lived, quiescent cell lineages capable of substantial proliferation upon acquisition of enabling genomic changes. Expansion of the dominant subclone to an appreciable mass may therefore represent the final rate-limiting step in a breast cancer's development, triggering diagnosis.