The perinucleolar compartment (PNC) is a unique nuclear substructure, forming predominantly in cancer cells both in vitro and in vivo. PNC prevalence (percentage of cells containing at least one PNC) has been found to positively correlate with disease progression in several cancers (breast, ovarian, and colon). While there is a clear association between PNCs and cancer, the molecular function of the PNC remains unclear. Here we summarize the current understanding of the association of PNCs with cancer and its possible functions in cancer cells.
PNC; cancer; nuclear substructure; gene expression regulation; structure and function
The developmental dynamics of multicellular organisms is a process that takes place in a multi-stable system in which each attractor state represents a cell type, and attractor transitions correspond to cell differentiation paths. This new understanding has revived the idea of a quasi-potential landscape, first proposed by Waddington as a metaphor. To describe development, one is interested in the ‘relative stabilities’ of N attractors (N > 2). Existing theories of state transition between local minima on some potential landscape deal with the exit part in the transition between two attractors in pair-attractor systems but do not offer the notion of a global potential function that relates more than two attractors to each other. Several ad hoc methods have been used in systems biology to compute a landscape in non-gradient systems, such as gene regulatory networks. Here we present an overview of currently available methods, discuss their limitations and propose a new decomposition of vector fields that permits the computation of a quasi-potential function that is equivalent to the Freidlin–Wentzell potential but is not limited to two attractors. Several examples of decomposition are given, and the significance of such a quasi-potential function is discussed.
multi-stable dynamical system; non-equilibrium dynamics; quasi-potential; state transition; epigenetic landscape; Freidlin–Wentzell theory
Changes in extracellular matrix (ECM) structure or mechanics can actively drive cancer progression; however, the underlying mechanism remains unknown. Here we explore whether this process could be mediated by changes in cell shape that lead to increases in genetic noise, given that both factors have been independently shown to alter gene expression and induce cell fate switching. We do this using a computer simulation model that explores the impact of physical changes in the tissue microenvironment under conditions in which physical deformation of cells increases gene expression variability among genetically identical cells. The model reveals that cancerous tissue growth can be driven by physical changes in the microenvironment: when increases in cell shape variability due to growth-dependent increases in cell packing density enhance gene expression variation, heterogeneous autonomous growth and further structural disorganization can result, thereby driving cancer progression via positive feedback. The model parameters that led to this prediction are consistent with experimental measurements of mammary tissues that spontaneously undergo cancer progression in transgenic C3(1)-SV40Tag female mice, which exhibit enhanced stiffness of mammary ducts, as well as progressive increases in variability of cell-cell relations and associated cell shape changes. These results demonstrate the potential for physical changes in the tissue microenvironment (e.g., altered ECM mechanics) to induce a cancerous phenotype or accelerate cancer progression in a clonal population through local changes in cell geometry and increased phenotypic variability, even in the absence of gene mutation.
Polarization of cell phenotypes, a common strategy to achieve cell type diversity in metazoa, results from binary cell-fate decisions in the branching pedigree of development. Such “either-or” fate decisions are controlled by two opposing cell fate-determining transcription factors. Each of the two distinct “master regulators” promotes differentiation of its respective sister lineage. But they also suppress one other, leading to their mutually exclusive expression in the two ensuing lineages. Thus, promiscuous coexistence of the antagonist regulators in the same cell, the hallmark of the common “undecided” progenitor of two sister lineages, is considered unstable. This antagonism ensures robust polarization into two discretely distinct cell types. But now the immune system's T-helper (Th) cells and their two canonical subtypes, Th1 and Th2 cells, tell a different story, as revealed in three papers recently published in PLOS Biology. The intermediate state that co-expresses the two opposing master regulators of the Th1 and Th2 subtypes, T-bet and Gata3, is highly stable and is not necessarily an undecided precursor. Instead, the Th1/Th2 hybrid cell is a robust new type with properties of both Th1 and Th2 cells. These hybrid cells are functionally active and possess the benefit of moderation: self-limitation of effector T cell function to prevent excessive inflammation, a permanent risk in host defense that can cause tissue damage or autoimmunity. Gene regulatory network analysis suggests that stabilization of the intermediate center in a polarizing system can be achieved by minor tweaking of the architecture of the mutual suppression gene circuit, and thus is a design option readily available to evolution.
Alu repeats within human genes may potentially alter gene expression. Here, we show that 3′-UTR-located inverted Alu repeats significantly reduce expression of an AcGFP reporter gene. Mutational analysis demonstrates that the secondary structure, but not the primary nucleotide sequence, of the inverted Alu repeats is critical for repression. The expression levels and nucleocytoplasmic distribution of reporter mRNAs with or without 3′-UTR inverted Alu repeats are similar; suggesting that reporter gene repression is not due to changes in mRNA levels or mRNA nuclear sequestration. Instead, reporter gene mRNAs harboring 3′-UTR inverted Alu repeats accumulate in cytoplasmic stress granules. These findings may suggest a novel mechanism whereby 3′-UTR-located inverted Alu repeats regulate human gene expression through sequestration of mRNAs within stress granules.
3′ untranslated region; Alu; inverted repeats; mRNA; stress granules
Multifactorial injuries, such as ischemia, trauma, etc., have proven stubbornly elusive to clinical therapeutics, in spite of the binary outcome of recovery or death. This may be due, in part, to the lack of formal approaches to cell injury. We present a minimal system of nonlinear ordinary differential equations describing a theory of cell injury dynamics. A mutual antagonism between injury-driven total damage and total induced stress responses gives rise to attractors representing recovery or death. Solving across a range of injury magnitudes defines an ‘injury course' containing a well-defined tipping point between recovery and death. Via the model, therapeutics is the diverting of a system on a pro-death trajectory to a pro-survival trajectory on bistable phase planes. The model plausibly explains why laboratory-based therapies have tended to fail clinically. A survival outcome is easy to achieve when lethal injury is close to the tipping point, but becomes progressively difficult as injury magnitudes increase, and there is an upper limit to salvageable injuries. The model offers novel insights into cell injury that may assist in overcoming barriers that have prevented development of clinically effective therapies for multifactorial conditions, as exemplified by brain ischemia.
bistability; brain ischemia; cell injury; nonlinear dynamics; therapeutics
It has been shown that nerve growth factor-β (NGF-β) promoted the initiation and progression of many tumors, and we have previously demonstrated that the expression of NGF-β was associated with tumor stage, nerve infiltration and lymph node metastasis in human hilar cholangiocarcinoma. However, whether NGF-β promotes tumor progression in human cholangiocarcinoma requires further investigation. Therefore, we aimed to determine the effects of NGF-β on the progression of human cholangiocarcinoma.
Human cholangiocarcinoma QBC939 stable cell lines with over-expressed or silenced NGF-β genes were generated with pEGFP-N1-NGF-β and pGPU6/GFP/Neo-NGF-β-shRNA recombinant plasmids. Cell proliferation assay, colony formation assay, cell cycle analysis, apoptosis assay and tumorigenicity assay were performed to evaluate the role of NGF-β in the progression of human cholangiocarcinoma. In addition, human lymphatic endothelial cells were co-cultured with QBC939 culture supernatants, and the cell proliferation and migration abilities of the lymphatic endothelial cells were evaluated.
Forced expression of NGF-β in QBC939 cell lines promoted proliferation, colony formation and tumorigenicity in these cells and inhibited the apoptosis. However, down-regulation of NGF-β inhibited proliferation, colony formation and tumorigenicity, and increased the apoptotic rate of QBC939 cells. In addition, the NGF-β gain-of-function induced a high expression of vascular endothelial growth factor C and enhanced the proliferation and migration of lymphatic endothelial cells, while NGF-β loss-of-function showed opposite effects.
We concluded that NGF-β promoted tumor progression in human cholangiocarcinoma QBC939 cells. Our results provided a new concept to understand the role of NGF-β in cholangiocarcinoma progression, and might provide important information for the development of new targeted therapies in human cholangiocarcinoma.
Reverse Transcription - quantitative Polymerase Chain Reaction (RT-qPCR) is a standard technique in most laboratories. The selection of reference genes is essential for data normalization and the selection of suitable reference genes remains critical. Our aim was to 1) review the literature since implementation of the MIQE guidelines in order to identify the degree of acceptance; 2) compare various algorithms in their expression stability; 3) identify a set of suitable and most reliable reference genes for a variety of human cancer cell lines. A PubMed database review was performed and publications since 2009 were selected. Twelve putative reference genes were profiled in normal and various cancer cell lines (n = 25) using 2-step RT-qPCR. Investigated reference genes were ranked according to their expression stability by five algorithms (geNorm, Normfinder, BestKeeper, comparative ΔCt, and RefFinder). Our review revealed 37 publications, with two thirds patient samples and one third cell lines. qPCR efficiency was given in 68.4% of all publications, but only 28.9% of all studies provided RNA/cDNA amount and standard curves. GeNorm and Normfinder algorithms were used in 60.5% in combination. In our selection of 25 cancer cell lines, we identified HSPCB, RRN18S, and RPS13 as the most stable expressed reference genes. In the subset of ovarian cancer cell lines, the reference genes were PPIA, RPS13 and SDHA, clearly demonstrating the necessity to select genes depending on the research focus. Moreover, a cohort of at least three suitable reference genes needs to be established in advance to the experiments, according to the guidelines. For establishing a set of reference genes for gene normalization we recommend the use of ideally three reference genes selected by at least three stability algorithms. The unfortunate lack of compliance to the MIQE guidelines reflects that these need to be further established in the research community.
Sensory systems often detect multiple types of inputs. For example, a receptor in a cell-signaling system often binds multiple kinds of ligands, and sensory neurons can respond to different types of stimuli. How do sensory systems compare these different kinds of signals? Here, we consider this question in a class of sensory systems – including bacterial chemotaxis- which have a property known as fold-change detection: their output dynamics, including amplitude and response time, depends only on the relative changes in signal, rather than absolute changes, over a range of several decades of signal. We analyze how fold-change detection systems respond to multiple signals, using mathematical models. Suppose that a step of fold F1 is made in input 1, together with a step of F2 in input 2. What total response does the system provide? We show that when both input signals impact the same receptor with equal number of binding sites, the integrated response is multiplicative: the response dynamics depend only on the product of the two fold changes, F1F2. When the inputs bind the same receptor with different number of sites n1 and n2, the dynamics depend on a product of power laws, . Thus, two input signals which vary over time in an inverse way can lead to no response. When the two inputs affect two different receptors, other types of integration may be found and generally the system is not constrained to respond according to the product of the fold-change of each signal. These predictions can be readily tested experimentally, by providing cells with two simultaneously varying input signals. The present study suggests how cells can compare apples and oranges, namely by comparing each to its own background level, and then multiplying these two fold-changes.
The increasing availability and maturity of DNA microarray technology has led to an explosion of cancer profiling studies for identifying cancer biomarkers, and predicting treatment response. Uncovering complex relationships, however, remains the most challenging task as it requires compiling and efficiently querying data from various sources. Here, we describe the Stress Response Array Profiler (StRAP), an open-source, web-based resource for storage, profiling, visualization, and sharing of cancer genomic data. StRAP houses multi-cancer microarray data with major emphasis on radiotherapy studies, and takes a systems biology approach towards the integration, comparison, and cross-validation of multiple cancer profiling studies. The database is a comprehensive platform for comparative analysis of gene expression data. For effective use of arrays, we provide user-friendly and interactive visualization tools that can display the data and query results. StRAP is web-based, platform-independent, and freely accessible at http://strap.nci.nih.gov/.
Widespread unexplained variations in clinical practices and patient outcomes suggest major opportunities for improving the quality and safety of medical care. However, there is little consensus regarding how to best identify and disseminate healthcare improvements and a dearth of theory to guide the debate. Many consider multicenter randomized controlled trials to be the gold standard of evidence-based medicine, although results are often inconclusive or may not be generally applicable due to differences in the contexts within which care is provided. Increasingly, others advocate the use “quality improvement collaboratives”, in which multi-institutional teams share information to identify potentially better practices that are subsequently evaluated in the local contexts of specific institutions, but there is concern that such collaborative learning approaches lack the statistical rigor of randomized trials. Using an agent-based model, we show how and why a collaborative learning approach almost invariably leads to greater improvements in expected patient outcomes than more traditional approaches in searching simulated clinical fitness landscapes. This is due to a combination of greater statistical power and more context-dependent evaluation of treatments, especially in complex terrains where some combinations of practices may interact in affecting outcomes. The results of our simulations are consistent with observed limitations of randomized controlled trials and provide important insights into probable reasons for effectiveness of quality improvement collaboratives in the complex socio-technical environments of healthcare institutions. Our approach illustrates how modeling the evolution of medical practice as search on a clinical fitness landscape can aid in identifying and understanding strategies for improving the quality and safety of medical care.
AIM: To investigate the effect and the possible mechanism of ginsenoside Rb1 on small intestinal smooth muscle motility in mice.
METHODS: Intestinal smooth muscle strips were isolated from male ICR mice (5 wk old), and the effect of ginsenoside Rb1 on spontaneous contraction was recorded with an electrophysiolograph. The effect of ginsenoside Rb1 on ion channel currents, including the voltage-gated K+ channel current (IKV), calcium-activated potassium channel currents (IKCa), spontaneous transient outward currents and ATP-sensitive potassium channel current (IKATP), was recorded on freshly isolated single cells using the whole-cell patch clamp technique.
RESULTS: Ginsenoside Rb1 dose-dependently inhibited the spontaneous contraction of intestinal smooth muscle by 21.15% ± 3.31%, 42.03% ± 8.23% and 67.23% ± 5.63% at concentrations of 25 μmol/L, 50 μmol/L and 100 μmol/L, respectively (n = 5, P < 0.05). The inhibitory effect of ginsenoside Rb1 on spontaneous contraction was significantly but incompletely blocked by 10 mmol/L tetraethylammonium or 0.5 mmol/L 4-aminopyridine, respectively (n = 5, P < 0.05). However, the inhibitory effect of ginsenoside Rb1 on spontaneous contraction was not affected by 10 μmol/L glibenclamide or 0.4 μmol/L tetrodotoxin. At the cell level, ginsenoside Rb1 increased outward potassium currents, and IKV was enhanced from 1137.71 ± 171.62 pA to 1449.73 ± 162.39 pA by 50 μmol/L Rb1 at +60 mV (n = 6, P < 0.05). Ginsenoside Rb1 increased IKCa and enhanced the amplitudes of spontaneous transient outward currents from 582.77 ± 179.09 mV to 788.12 ± 278.34 mV (n = 5, P < 0.05). However, ginsenoside Rb1 (50 μmol/L) had no significant effect on IKATP (n = 3, P < 0.05).
CONCLUSION: These results suggest that ginsenoside Rb1 has an inhibitory effect on the spontaneous contraction of mouse intestinal smooth muscle mediated by the activation of IKV and IKCa, but the KATP channel was not involved in this effect.
Ginsenoside Rb1; Intestinal smooth muscle; Intestinal smooth muscle cell; Potassium channel; Spontaneous contraction; Whole-cell patch clamp technique
Accumulating experimental evidence suggests that the gene regulatory networks of living organisms operate in the critical phase, namely, at the transition between ordered and chaotic dynamics. Such critical dynamics of the network permits the coexistence of robustness and flexibility which are necessary to ensure homeostatic stability (of a given phenotype) while allowing for switching between multiple phenotypes (network states) as occurs in development and in response to environmental change. However, the mechanisms through which genetic networks evolve such critical behavior have remained elusive. Here we present an evolutionary model in which criticality naturally emerges from the need to balance between the two essential components of evolvability: phenotype conservation and phenotype innovation under mutations. We simulated the Darwinian evolution of random Boolean networks that mutate gene regulatory interactions and grow by gene duplication. The mutating networks were subjected to selection for networks that both (i) preserve all the already acquired phenotypes (dynamical attractor states) and (ii) generate new ones. Our results show that this interplay between extending the phenotypic landscape (innovation) while conserving the existing phenotypes (conservation) suffices to cause the evolution of all the networks in a population towards criticality. Furthermore, the networks produced by this evolutionary process exhibit structures with hubs (global regulators) similar to the observed topology of real gene regulatory networks. Thus, dynamical criticality and certain elementary topological properties of gene regulatory networks can emerge as a byproduct of the evolvability of the phenotypic landscape.
Dynamically critical systems are those which operate at the border of a phase transition between two behavioral regimes often present in complex systems: order and disorder. Critical systems exhibit remarkable properties such as fast information processing, collective response to perturbations or the ability to integrate a wide range of external stimuli without saturation. Recent evidence indicates that the genetic networks of living cells are dynamically critical. This has far reaching consequences, for it is at criticality that living organisms can tolerate a wide range of external fluctuations without changing the functionality of their phenotypes. Therefore, it is necessary to know how genetic criticality emerged through evolution. Here we show that dynamical criticality naturally emerges from the delicate balance between two fundamental forces of natural selection that make organisms evolve: (i) the existing phenotypes must be resilient to random mutations, and (ii) new phenotypes must emerge for the organisms to adapt to new environmental challenges. The joint effect of these two forces, which are essential for evolvability, is sufficient in our computational models to generate populations of genetic networks operating at criticality. Thus, natural selection acting as a tinkerer of evolvable systems naturally generates critical dynamics.
Stem cell behaviours, such as stabilization of the undecided state of pluripotency or multipotency, the priming towards a prospective fate, binary fate decisions and irreversible commitment, must all somehow emerge from a genome-wide gene-regulatory network. Its unfathomable complexity defies the standard mode of explanation that is deeply rooted in molecular biology thinking: the reduction of observables to linear deterministic molecular pathways that are tacitly taken as chains of causation. Such culture of proximate explanation that uses qualitative arguments, simple arrow–arrow schemes or metaphors persists despite the ceaseless accumulation of ‘omics’ data and the rise of systems biology that now offers precise conceptual tools to explain emergent cell behaviours from gene networks. To facilitate the embrace of the principles of physics and mathematics that underlie such systems and help to bridge the gap between the formal description of theorists and the intuition of experimental biologists, we discuss in qualitative terms three perspectives outside the realm of their familiar linear-deterministic view: (i) state space (ii), high-dimensionality and (iii) heterogeneity. These concepts jointly offer a new vista on stem cell regulation that naturally explains many novel, counterintuitive observations and their inherent inevitability, obviating the need for ad hoc explanations of their existence based on natural selection. Hopefully, this expanded view will stimulate novel experimental designs.
gene-regulatory networks; dynamics; high-dimensional state space; heterogeneity; cell fate
Gradual or sudden transitions among different states as exhibited by cell populations in a biological sample under particular conditions or stimuli can be detected and profiled by flow cytometric time course data. Often such temporal profiles contain features due to transient states that present unique modeling challenges. These could range from asymmetric non-Gaussian distributions to outliers and tail subpopulations, which need to be modeled with precision and rigor.
To ensure precision and rigor, we propose a parametric modeling framework StateProfiler based on finite mixtures of skew t-Normal distributions that are robust against non-Gaussian features caused by asymmetry and outliers in data. Further, we present in StateProfiler a new greedy EM algorithm for fast and optimal model selection. The parsimonious approach of our greedy algorithm allows us to detect the genuine dynamic variation in the key features as and when they appear in time course data. We also present a procedure to construct a well-fitted profile by merging any redundant model components in a way that minimizes change in entropy of the resulting model. This allows precise profiling of unusually shaped distributions and less well-separated features that may appear due to cellular heterogeneity even within clonal populations.
By modeling flow cytometric data measured over time course and marker space with StateProfiler, specific parametric characteristics of cellular states can be identified. The parameters are then tested statistically for learning global and local patterns of spatio-temporal change. We applied StateProfiler to identify the temporal features of yeast cell cycle progression based on knockout of S-phase triggering cyclins Clb5 and Clb6, and then compared the S-phase delay phenotypes due to differential regulation of the two cyclins. We also used StateProfiler to construct the temporal profile of clonal divergence underlying lineage selection in mammalian hematopoietic progenitor cells.
Demand for high quality gene expression data has driven the development of revolutionary microarray technologies. The quality of the data is affected by the performance of the microarray platform as well as how the nucleic acid targets are prepared. The most common method for target nucleic acid preparation includes in vitro transcription amplification of the sample RNA. Although this method requires a small amount of starting material and is reported to have high reproducibility, there are also technical disadvantages such as amplification bias and the long, laborious protocol. Using RNA derived from human brain, breast and colon, we demonstrate that a non-amplification method, which was previously shown to be inferior, could be transformed to a highly quantitative method with a dynamic range of five orders of magnitude. Furthermore, the correlation coefficient calculated by comparing microarray assays using non-amplified samples with qRT-PCR assays was approximately 0.9, a value much higher than when samples were prepared using amplification methods. Our results were also compared with data from various microarray platforms studied in the MicroArray Quality Control (MAQC) project. In combination with micro-columnar 3D-Gene™ microarray, this non-amplification method is applicable to a variety of genetic analyses, including biomarker screening and diagnostic tests for cancer.
Head and neck squamous carcinoma (HNSCC) tumors carry dismal long-term prognosis and the role of tumor initiating cells (TICs) in this cancer is unclear. We investigated in HNSCC xenografts whether specific tumor subpopulations contributed to tumor growth. We used a CFSE-based label retentions assay, CD49f (α6-integrin) surface levels and aldehyde dehydrogenase (ALDH) activity to profile HNSCC subpopulations. The tumorigenic potential of marker-positive and -negative subpopulations was tested in nude (Balb/c nu/nu) and NSG (NOD.Cg-Prkdcscid Il2rgtm1Wjl/SzJ) mice and chicken embryo chorioallantoic membrane (CAM) assays. Here we identified in HEp3, SQ20b and FaDu HNSCC xenografts a subpopulation of G0/G1-arrested slow-cycling CD49fhigh/ALDH1A1high/H3K4/K27me3low subpopulation (CD49f+) of tumor cells. A strikingly similar CD49fhigh/H3K27me3low subpopulation is also present in primary human HNSCC tumors and metastases. While only sorted CD49fhigh/ALDHhigh, label retaining cells (LRC) proliferated immediately in vivo, with time the CD49flow/ALDHlow, non-LRC (NLRC) tumor cell subpopulations were also able to regain tumorigenic capacity; this was linked to restoration of CD49fhigh/ALDHhigh, label retaining cells. In addition, CD49f is required for HEp3 cell tumorigenicity and to maintain low levels of H3K4/K27me3. CD49f+ cells also displayed reduced expression of the histone-lysine N-methyltransferase EZH2 and ERK1/2phosphorylation. This suggests that although transiently quiescent, their unique chromatin structure is poised for rapid transcriptional activation. CD49f− cells can “reprogram” and also achieve this state eventually. We propose that in HNSCC tumors, epigenetic mechanisms likely driven by CD49f signaling dynamically regulate HNSCC xenograft phenotypic heterogeneity. This allows multiple tumor cell subpopulations to drive tumor growth suggesting that their dynamic nature renders them a “moving target” and their eradication might require more persistent strategies.
Epoxyeicosatrienoic acids (EETs) are small molecules produced by cytochrome P450 epoxygenases. They are lipid mediators that act as autocrine or paracrine factors to regulate inflammation and vascular tone. As a result, drugs that raise EET levels are in clinical trials for the treatment of hypertension and many other diseases. However, despite their pleiotropic effects on cells, little is known about the role of these epoxyeicosanoids in cancer. Here, using genetic and pharmacological manipulation of endogenous EET levels, we demonstrate that EETs are critical for primary tumor growth and metastasis in a variety of mouse models of cancer. Remarkably, we found that EETs stimulated extensive multiorgan metastasis and escape from tumor dormancy in several tumor models. This systemic metastasis was not caused by excessive primary tumor growth but depended on endothelium-derived EETs at the site of metastasis. Administration of synthetic EETs recapitulated these results, while EET antagonists suppressed tumor growth and metastasis, demonstrating in vivo that pharmacological modulation of EETs can affect cancer growth. Furthermore, inhibitors of soluble epoxide hydrolase (sEH), the enzyme that metabolizes EETs, elevated endogenous EET levels and promoted primary tumor growth and metastasis. Thus, our data indicate a central role for EETs in tumorigenesis, offering a mechanistic link between lipid signaling and cancer and emphasizing the critical importance of considering possible effects of EET-modulating drugs on cancer.
Acquiring abundant and high-purity cancer stem cells (CSCs) is an important prerequisite for CSC research. At present, researchers usually gain high-purity CSCs through flow cytometry sorting and expand them by short-term sphere culture. However, it is still uncertain whether we can amplify high-purity CSCs through long-term sphere culture. We have proposed a mathematical model using ordinary differential equations to derive the continuous variation of the CSC ratio in long-term sphere culture and estimated the model parameters based on a long-term sphere culture of MCF-7 stem cells. We found that the CSC ratio in long-term sphere culture presented as gradually decreased drift and might be stable at a lower level. Furthermore, we found that fitted model parameters could explain the main growth pattern of CSCs and differentiated cancer cells in long-term sphere culture.
The PNC is a nuclear body that forms in malignant cells. We characterize a newly identified complex in the PNC; determine the dynamics of PNC proteins; and describe the recruitment, localization, and sedimentation properties of PNC components.
The perinucleolar compartment (PNC) forms in cancer cells and is highly enriched with a subset of polymerase III RNAs and RNA-binding proteins. Here we report that PNC components mitochondrial RNA–processing (MRP) RNA, pyrimidine tract–binding protein (PTB), and CUG-binding protein (CUGBP) interact in vivo, as demonstrated by coimmunoprecipitation and RNA pull-down experiments. Glycerol gradient analyses show that this complex is large and sediments at a different fraction from known MRP RNA–containing complexes, the MRP ribonucleoprotein ribozyme and human telomerase reverse transcriptase. Tethering PNC components to a LacO locus recruits other PNC components, further confirming the in vivo interactions. These interactions are present both in PNC-containing and -lacking cells. High-resolution localization analyses demonstrate that MRP RNA, CUGBP, and PTB colocalize at the PNC as a reticulated network, intertwining with newly synthesized RNA. Furthermore, green fluorescent protein (GFP)–PTB and GFP-CUGBP show a slower rate of fluorescence recovery after photobleaching at the PNC than in the nucleoplasm, illustrating the different molecular interaction of the complexes associated with the PNC. These findings support a working model in which the MRP RNA–protein complex becomes nucleated at the PNC in cancer cells and may play a role in gene expression regulation at the DNA locus that associates with the PNC.
The gene regulatory circuit motif in which two opposing fate-determining transcription factors inhibit each other but activate themselves has been used in mathematical models of binary cell fate decisions in multipotent stem or progenitor cells. This simple circuit can generate multistability and explains the symmetric “poised” precursor state in which both factors are present in the cell at equal amounts as well as the resolution of this indeterminate state as the cell commits to either cell fate characterized by an asymmetric expression pattern of the two factors. This establishes the two alternative stable attractors that represent the two fate options. It has been debated whether cooperativity of molecular interactions is necessary to produce such multistability.
Here we take a general modeling approach and argue that this question is not relevant. We show that non-linearity can arise in two distinct models in which no explicit interaction between the two factors is assumed and that distinct chemical reaction kinetic formalisms can lead to the same (generic) dynamical system form. Moreover, we describe a novel type of bifurcation that produces a degenerate steady state that can explain the metastable state of indeterminacy prior to cell fate decision-making and is consistent with biological observations.
The general model presented here thus offers a novel principle for linking regulatory circuits with the state of indeterminacy characteristic of multipotent (stem) cells.
Amonafide is a DNA intercalator in clinical development for the treatment of cancer. The drug has a 5-position amine that is variably acetylated to form a toxic metabolite in humans, increasing adverse effects and complicating the dosing of amonafide. Numonafides, 6-amino derivatives of amonafide that avoid the toxic acetylation, also show in vitro anticancer activity, as we have previously described. Here, we report the in vitro and in vivo activities of two numonafides, 6-methoxyethylamino-numonafide (MEAN) and 6-amino-numonafide (AN) with comparisons to amonafide. The in vitro potencies and cellular anticancer mechanisms are similar for the two numonafides and amonafide. Results from several mouse models of human cancer demonstrate that AN and MEAN require slightly higher doses than amonafide for equal efficacy in short-term dosing models, but the same dose of all three compounds in long-term dosing models are equally efficacious. MEAN is tolerated much better than amonafide and AN at equally efficacious doses based on weight change, activity, stool consistency, and dose tolerance with survival as the end point. The studies presented here demonstrate that MEAN is much less toxic than amonafide or AN in mouse models of human liver and gastric cancers while being equally efficacious in vivo and inhibiting cancer cells through similar mechanisms. These findings demonstrate that numonafides can be less toxic than amonafide and support further preclinical development and novel anticancer agents or as replacements or amonafide.
Predator-prey system, as an essential element of ecological dynamics, has been recently studied experimentally with synthetic biology. We developed a global probabilistic landscape and flux framework to explore a synthetic predator-prey network constructed with two Escherichia coli populations. We developed a self consistent mean field method to solve multidimensional problem and uncovered the potential landscape with Mexican hat ring valley shape for predator-prey oscillations. The landscape attracts the system down to the closed oscillation ring. The probability flux drives the coherent oscillations on the ring. Both the landscape and flux are essential for the stable and coherent oscillations. The landscape topography characterized by the barrier height from the top of Mexican hat to the closed ring valley provides a quantitative measure of global stability of system. The entropy production rate for the energy dissipation is less for smaller environmental fluctuations or perturbations. The global sensitivity analysis based on the landscape topography gives specific predictions for the effects of parameters on the stability and function of the system. This may provide some clues for the global stability, robustness, function and synthetic network design.
Cell fate reprogramming, such as the generation of insulin-producing β cells from other pancreas cells, can be achieved by external modulation of key transcription factors. However, the known gene regulatory interactions that form a complex network with multiple feedback loops make it increasingly difficult to design the cell reprogramming scheme because the linear regulatory pathways as schemes of causal influences upon cell lineages are inadequate for predicting the effect of transcriptional perturbation. However, sufficient information on regulatory networks is usually not available for detailed formal models. Here we demonstrate that by using the qualitatively described regulatory interactions as the basis for a coarse-grained dynamical ODE (ordinary differential equation) based model, it is possible to recapitulate the observed attractors of the exocrine and β, δ, α endocrine cells and to predict which gene perturbation can result in desired lineage reprogramming. Our model indicates that the constraints imposed by the incompletely elucidated regulatory network architecture suffice to build a predictive model for making informed decisions in choosing the set of transcription factors that need to be modulated for fate reprogramming.