Systematic search for recipes to generate induced pluripotent stem cells

Modulating the dynamics of NFκB and PI3K enhances the ensemble-level TNFR1 signaling mediated apoptotic response

Cell-to-cell variability during TNFα stimulated Tumor Necrosis Factor Receptor 1 (TNFR1) signaling can lead to single-cell level pro-survival and apoptotic responses. This variability stems from the heterogeneity in signal flow through intracellular signaling entities that regulate the balance between these two phenotypes. Using systematic Boolean dynamic modeling of a TNFR1 signaling network, we demonstrate that the signal flow path variability can be modulated to enable cells favour apoptosis. We developed a computationally efficient approach "Boolean Modeling based Prediction of Steady-state probability of Phenotype Reachability (BM-ProSPR)" to accurately predict the network's ability to settle into different phenotypes. Model analysis juxtaposed with the experimental observations revealed that NFκB and PI3K transient responses guide the XIAP behaviour to coordinate the crucial dynamic cross-talk between the pro-survival and apoptotic arms at the single-cell level. Model predicted the experimental observations that ~31% apoptosis increase can be achieved by arresting Comp1 - IKK* activity which regulates the NFκB and PI3K dynamics. Arresting Comp1 - IKK* activity causes signal flow path re-wiring towards apoptosis without significantly compromising NFκB levels, which govern adequate cell survival. Priming an ensemble of cancerous cells with inhibitors targeting the specific interaction involving Comp1 and IKK* prior to TNFα exposure could enable driving them towards apoptosis.

Minimal frustration underlies the usefulness of incomplete regulatory network models in biology

Regulatory networks as large and complex as those implicated in cell-fate choice are expected to exhibit intricate, very high-dimensional dynamics. Cell-fate choice, however, is a macroscopically simple process. Additionally, regulatory network models are almost always incomplete and/or inexact, and do not incorporate all the regulators and interactions that may be involved in cell-fate regulation. In spite of these issues, regulatory network models have proven to be incredibly effective tools for understanding cell-fate choice across contexts and for making useful predictions. Here, we show that minimal frustration-a feature of biological networks across contexts but not of random networks-can compel simple, low-dimensional steady-state behavior even in large and complex networks. Moreover, the steady-state behavior of minimally frustrated networks can be recapitulated by simpler networks such as those lacking many of the nodes and edges and those that treat multiple regulators as one. The present study provides a theoretical explanation for the success of network models in biology and for the challenges in network inference.

Covariant influences for finite discrete dynamical systems

We develop a rigorous theory of external influences on finite discrete dynamical systems, going beyond the perturbation paradigm, in that the external influence need not be a small contribution. Indeed, the covariance condition can be stated as follows: If we evolve the dynamical system for n time steps and then disturb it, then it is the same as first disturbing the system with the same influence and then letting the system evolve for n time steps. Applying the powerful machinery of resource theories, we develop a theory of covariant influences both when there is a purely deterministic evolution and when randomness is involved. Subsequently, we provide necessary and sufficient conditions for the transition between states under deterministic covariant influences and necessary conditions in the presence of stochastic covariant influences, predicting which transitions between states are forbidden. Our approach, for the first time, employs the framework of resource theories, borrowed from quantum information theory, to the study of finite discrete dynamical systems. The laws we articulate unify the behavior of different types of finite discrete dynamical systems, and their mathematical flavor makes them rigorous and checkable.

Distance covariance entropy reveals primed states and bifurcation dynamics in single-cell RNA-Seq data

Cell-fate transitions are fundamental to development and differentiation. Studying them with single-cell omic data is important to advance our understanding of the cell-fate commitment process, yet this remains challenging. Here we present a computational method called DICE, which analyzes the entropy of expression covariation patterns and which is applicable to static and dynamically changing cell populations. Using only single-cell RNA-Seq data, DICE is able to predict multipotent primed states and their regulatory factors, which we subsequently validate with single-cell epigenomic data. DICE reveals that primed states are often defined by epigenetic regulators or pioneer factors alongside lineage-specific transcription factors. In developmental time course single-cell RNA-Seq datasets, DICE can pinpoint the timing of bifurcations more precisely than lineage-trajectory inference algorithms or competing variance-based methods. In summary, by studying the dynamic changes of expression covariation entropy, DICE can help elucidate primed states and bifurcation dynamics without the need for single-cell epigenomic data.

Finite-Time Set Reachability of Probabilistic Boolean Multiplex Control Networks

This study focuses on the finite-time set reachability of probabilistic Boolean multiplex control networks (PBMCNs). Firstly, based on the state transfer graph (STG) reconstruction technique, the PBMCNs are extended to random logic dynamical systems. Then, a necessary and sufficient condition for the finite-time set reachability of PBMCNs is obtained. Finally, the obtained results are effectively illustrated by an example.

Effective connectivity determines the critical dynamics of biochemical networks

Living systems comprise interacting biochemical components in very large networks. Given their high connectivity, biochemical dynamics are surprisingly not chaotic but quite robust to perturbations-a feature C.H. Waddington named canalization. Because organisms are also flexible enough to evolve, they arguably operate in a critical dynamical regime between order and chaos. The established theory of criticality is based on networks of interacting automata where Boolean truth values model presence/absence of biochemical molecules. The dynamical regime is predicted using network connectivity and node bias (to be on/off) as tuning parameters. Revising this to account for canalization leads to a significant improvement in dynamical regime prediction. The revision is based on effective connectivity, a measure of dynamical redundancy that buffers automata response to some inputs. In both random and experimentally validated systems biology networks, reducing effective connectivity makes living systems operate in stable or critical regimes even though the structure of their biochemical interaction networks predicts them to be chaotic. This suggests that dynamical redundancy may be naturally selected to maintain living systems near critical dynamics, providing both robustness and evolvability. By identifying how dynamics propagates preferably via effective pathways, our approach helps to identify precise ways to design and control network models of biochemical regulation and signalling.

Predictive network modeling in human induced pluripotent stem cells identifies key driver genes for insulin responsiveness

Insulin resistance (IR) precedes the development of type 2 diabetes (T2D) and increases cardiovascular disease risk. Although genome wide association studies (GWAS) have uncovered new loci associated with T2D, their contribution to explain the mechanisms leading to decreased insulin sensitivity has been very limited. Thus, new approaches are necessary to explore the genetic architecture of insulin resistance. To that end, we generated an iPSC library across the spectrum of insulin sensitivity in humans. RNA-seq based analysis of 310 induced pluripotent stem cell (iPSC) clones derived from 100 individuals allowed us to identify differentially expressed genes between insulin resistant and sensitive iPSC lines. Analysis of the co-expression architecture uncovered several insulin sensitivity-relevant gene sub-networks, and predictive network modeling identified a set of key driver genes that regulate these co-expression modules. Functional validation in human adipocytes and skeletal muscle cells (SKMCs) confirmed the relevance of the key driver candidate genes for insulin responsiveness.

Insulin resistance is characterized by a defective response (“resistance”) to normal insulin concentrations to uptake the glucose present in the blood, and is the underlying condition that leads to type 2 diabetes (T2D) and increases the risk of cardiovascular disease. It is estimated that 25–33% of the US population are insulin resistant enough to be at risk of serious clinical consequences. For more than a decade, large population studies have tried to discover the genes that participate in the development of insulin resistance, but without much success. It is now increasingly clear that the complex genetic nature of insulin resistance requires novel approaches centered in patient specific cellular models. To fill this gap, we have generated an induced pluripotent stem cell (iPSC) library from individuals with accurate measurements of insulin sensitivity, and performed gene expression and key driver analyses. Our work demonstrates that iPSCs can be used as a revolutionary technology to model insulin resistance and to discover key genetic drivers. Moreover, they can develop our basic knowledge of the disease, and are ultimately expected to increase the therapeutic targets to treat insulin resistance and type 2 diabetes.

Biological Networks Regulating Cell Fate Choice Are Minimally Frustrated

Characterization of the differences between biological and random networks can reveal the design principles that enable the robust realization of crucial biological functions including the establishment of different cell types. Previous studies, focusing on identifying topological features that are present in biological networks but not in random networks, have, however, provided few functional insights. We use a Boolean modeling framework and ideas from the spin glass literature to identify functional differences between five real biological networks and random networks with similar topological features. We show that minimal frustration is a fundamental property that allows biological networks to robustly establish cell types and regulate cell fate choice, and that this property can emerge in complex networks via Darwinian evolution. The study also provides clues regarding how the regulation of cell fate choice can go awry in a disease like cancer and lead to the emergence of aberrant cell types.

Computational advances of tumor marker selection and sample classification in cancer proteomics

Cancer proteomics has become a powerful technique for characterizing the protein markers driving transformation of malignancy, tracing proteome variation triggered by therapeutics, and discovering the novel targets and drugs for the treatment of oncologic diseases. To facilitate cancer diagnosis/prognosis and accelerate drug target discovery, a variety of methods for tumor marker identification and sample classification have been developed and successfully applied to cancer proteomic studies. This review article describes the most recent advances in those various approaches together with their current applications in cancer-related studies. Firstly, a number of popular feature selection methods are overviewed with objective evaluation on their advantages and disadvantages. Secondly, these methods are grouped into three major classes based on their underlying algorithms. Finally, a variety of sample separation algorithms are discussed. This review provides a comprehensive overview of the advances on tumor maker identification and patients/samples/tissues separations, which could be guidance to the researches in cancer proteomics.

Algorithms for the Sequential Reprogramming of Boolean Networks

Cellular reprogramming, a technique that opens huge opportunities in modern and regenerative medicine, heavily relies on identifying key genes to perturb. Most of the existing computational methods for controlling which attractor (steady state) the cell will reach focus on finding mutations to apply to the initial state. However, it has been shown, and is proved in this article, that waiting between perturbations so that the update dynamics of the system prepares the ground, allows for new reprogramming strategies. To identify such sequential perturbations, we consider a qualitative model of regulatory networks, and rely on Binary Decision Diagrams to model their dynamics and the putative perturbations. Our method establishes a set identification of sequential perturbations, whether permanent (mutations) or only temporary, to achieve the existential or inevitable reachability of an arbitrary state of the system. We apply an implementation for temporary perturbations on models from the literature, illustrating that we are able to derive sequential perturbations to achieve trans-differentiation.

Exploring induced pluripotency in human fibroblasts via construction, validation, and application of a gene regulatory network

Reprogramming of somatic cells to induced pluripotent stem cells, by overexpressing certain factors referred to as the reprogramming factors, can revolutionize regenerative medicine. To provide a coherent description of induced pluripotency from the gene regulation perspective, we use 35 microarray datasets to construct a reprogramming gene regulatory network. Comprising 276 nodes and 4471 links, the resulting network is, to the best of our knowledge, the largest gene regulatory network constructed for human fibroblast reprogramming and it is the only one built using a large number of experimental datasets. To build the network, a model that relates the expression profiles of the initial (fibroblast) and final (induced pluripotent stem cell) states is proposed and the model parameters (link strengths) are fitted using the experimental data. Twenty nine additional experimental datasets are collectively used to test the model/network, and good agreement between experimental and predicted gene expression profiles is found. We show that the model in conjunction with the constructed network can make useful predictions. For example, we demonstrate that our approach can incorporate the effect of reprogramming factor stoichiometry and that its predictions are consistent with the experimentally observed trends in reprogramming efficiency when the stoichiometric ratios vary. Using our model/network, we also suggest new (not used in training of the model) candidate sets of reprogramming factors, many of which have already been experimentally verified. These results suggest our model/network can potentially be used in devising new recipes for induced pluripotency with higher efficiencies. Additionally, we classify the links of the network into three classes of different importance, prioritizing them for experimental verification. We show that many of the links in the top ranked class are experimentally known to be important in reprogramming. Finally, comparing with other methods, we show that using our model is advantageous.

Ensembles, dynamics, and cell types: Revisiting the statistical mechanics perspective on cellular regulation

Genetic regulatory networks control ontogeny. For fifty years Boolean networks have served as models of such systems, ranging from ensembles of random Boolean networks as models for generic properties of gene regulation to working dynamical models of a growing number of sub-networks of real cells. At the same time, their statistical mechanics has been thoroughly studied. Here we recapitulate their original motivation in the context of current theoretical and empirical research. We discuss ensembles of random Boolean networks whose dynamical attractors model cell types. A sub-ensemble is the critical ensemble. There is now strong evidence that genetic regulatory networks are dynamically critical, and that evolution is exploring the critical sub-ensemble. The generic properties of this sub-ensemble predict essential features of cell differentiation. In particular, the number of attractors in such networks scales as the DNA content raised to the 0.63 power. Data on the number of cell types as a function of the DNA content per cell shows a scaling relationship of 0.88. Thus, the theory correctly predicts a power law relationship between the number of cell types and the DNA contents per cell, and a comparable slope. We discuss these new scaling values and show prospects for new research lines for Boolean networks as a base model for systems biology.

Rare-event sampling of epigenetic landscapes and phenotype transitions

Stochastic simulation has been a powerful tool for studying the dynamics of gene regulatory networks, particularly in terms of understanding how cell-phenotype stability and fate-transitions are impacted by noisy gene expression. However, gene networks often have dynamics characterized by multiple attractors. Stochastic simulation is often inefficient for such systems, because most of the simulation time is spent waiting for rare, barrier-crossing events to occur. We present a rare-event simulation-based method for computing epigenetic landscapes and phenotype-transitions in metastable gene networks. Our computational pipeline was inspired by studies of metastability and barrier-crossing in protein folding, and provides an automated means of computing and visualizing essential stationary and dynamic information that is generally inaccessible to conventional simulation. Applied to a network model of pluripotency in Embryonic Stem Cells, our simulations revealed rare phenotypes and approximately Markovian transitions among phenotype-states, occurring with a broad range of timescales. The relative probabilities of phenotypes and the transition paths linking pluripotency and differentiation are sensitive to global kinetic parameters governing transcription factor-DNA binding kinetics. Our approach significantly expands the capability of stochastic simulation to investigate gene regulatory network dynamics, which may help guide rational cell reprogramming strategies. Our approach is also generalizable to other types of molecular networks and stochastic dynamics frameworks.

Bayesian Networks Predict Neuronal Transdifferentiation

We employ the language of Bayesian networks to systematically construct gene-regulation topologies from deep-sequencing single-nucleus RNA-Seq data for human neurons. From the perspective of the cell-state potential landscape, we identify attractors that correspond closely to different neuron subtypes. Attractors are also recovered for cell states from an independent data set confirming our models accurate description of global genetic regulations across differing cell types of the neocortex (not included in the training data). Our model recovers experimentally confirmed genetic regulations and community analysis reveals genetic associations in common pathways. Via a comprehensive scan of all theoretical three-gene perturbations of gene knockout and overexpression, we discover novel neuronal trans-differrentiation recipes (including perturbations of SATB2, GAD1, POU6F2 and ADARB2) for excitatory projection neuron and inhibitory interneuron subtypes.

Nanog induced intermediate state in regulating stem cell differentiation and reprogramming

Background

Heterogeneous gene expressions of cells are widely observed in self-renewing pluripotent stem cells, suggesting possible coexistence of multiple cellular states with distinct characteristics. Though the elements regulating cellular states have been identified, the underlying dynamic mechanisms and the significance of such cellular heterogeneity remain elusive.

Results

We present a gene regulatory network model to investigate the bimodal Nanog distribution in stem cells. Our model reveals a novel role of dynamic conversion between the cellular states of high and low Nanog levels. Model simulations demonstrate that the low-Nanog state benefits cell differentiation through serving as an intermediate state to reduce the barrier of transition. Interestingly, the existence of low-Nanog state dynamically slows down the reprogramming process, and additional Nanog activation is found to be essential to quickly attaining the fully reprogrammed cell state.

Conclusions

Nanog has been recognized as a critical pluripotency gene in stem cell regulation. Our modeling results quantitatively show a dual role of Nanog during stem cell differentiation and reprogramming, and the importance of the intermediate state during cell state transitions. Our approach offers a general method for analyzing key regulatory factors controlling cell differentiation and reprogramming.

Electronic supplementary material

The online version of this article (10.1186/s12918-018-0552-3) contains supplementary material, which is available to authorized users.

Relative stability of network states in Boolean network models of gene regulation in development

Progress in cell type reprogramming has revived the interest in Waddington's concept of the epigenetic landscape. Recently researchers developed the quasi-potential theory to represent the Waddington's landscape. The Quasi-potential U(x), derived from interactions in the gene regulatory network (GRN) of a cell, quantifies the relative stability of network states, which determine the effort required for state transitions in a multi-stable dynamical system. However, quasi-potential landscapes, originally developed for continuous systems, are not suitable for discrete-valued networks which are important tools to study complex systems. In this paper, we provide a framework to quantify the landscape for discrete Boolean networks (BNs). We apply our framework to study pancreas cell differentiation where an ensemble of BN models is considered based on the structure of a minimal GRN for pancreas development. We impose biologically motivated structural constraints (corresponding to specific type of Boolean functions) and dynamical constraints (corresponding to stable attractor states) to limit the space of BN models for pancreas development. In addition, we enforce a novel functional constraint corresponding to the relative ordering of attractor states in BN models to restrict the space of BN models to the biological relevant class. We find that BNs with canalyzing/sign-compatible Boolean functions best capture the dynamics of pancreas cell differentiation. This framework can also determine the genes' influence on cell state transitions, and thus can facilitate the rational design of cell reprogramming protocols.

Effects of Collective Histone State Dynamics on Epigenetic Landscape and Kinetics of Cell Reprogramming

Cell reprogramming is a process of transitions from differentiated to pluripotent cell states via transient intermediate states. Within the epigenetic landscape framework, such a process is regarded as a sequence of transitions among basins on the landscape; therefore, theoretical construction of a model landscape which exhibits experimentally consistent dynamics can provide clues to understanding epigenetic mechanism of reprogramming. We propose a minimal gene-network model of the landscape, in which each gene is regulated by an integrated mechanism of transcription-factor binding/unbinding and the collective chemical modification of histones. We show that the slow collective variation of many histones around each gene locus alters topology of the landscape and significantly affects transition dynamics between basins. Differentiation and reprogramming follow different transition pathways on the calculated landscape, which should be verified experimentally via single-cell pursuit of the reprogramming process. Effects of modulation in collective histone state kinetics on transition dynamics and pathway are examined in search for an efficient protocol of reprogramming.

Discrete Logic Modelling Optimization to Contextualize Prior Knowledge Networks Using PRUNET

High-throughput technologies have led to the generation of an increasing amount of data in different areas of biology. Datasets capturing the cell's response to its intra- and extra-cellular microenvironment allows such data to be incorporated as signed and directed graphs or influence networks. These prior knowledge networks (PKNs) represent our current knowledge of the causality of cellular signal transduction. New signalling data is often examined and interpreted in conjunction with PKNs. However, different biological contexts, such as cell type or disease states, may have distinct variants of signalling pathways, resulting in the misinterpretation of new data. The identification of inconsistencies between measured data and signalling topologies, as well as the training of PKNs using context specific datasets (PKN contextualization), are necessary conditions to construct reliable, predictive models, which are current challenges in the systems biology of cell signalling. Here we present PRUNET, a user-friendly software tool designed to address the contextualization of a PKNs to specific experimental conditions. As the input, the algorithm takes a PKN and the expression profile of two given stable steady states or cellular phenotypes. The PKN is iteratively pruned using an evolutionary algorithm to perform an optimization process. This optimization rests in a match between predicted attractors in a discrete logic model (Boolean) and a Booleanized representation of the phenotypes, within a population of alternative subnetworks that evolves iteratively. We validated the algorithm applying PRUNET to four biological examples and using the resulting contextualized networks to predict missing expression values and to simulate well-characterized perturbations. PRUNET constitutes a tool for the automatic curation of a PKN to make it suitable for describing biological processes under particular experimental conditions. The general applicability of the implemented algorithm makes PRUNET suitable for a variety of biological processes, for instance cellular reprogramming or transitions between healthy and disease states.

Summary of the DREAM8 Parameter Estimation Challenge: Toward Parameter Identification for Whole-Cell Models

Whole-cell models that explicitly represent all cellular components at the molecular level have the potential to predict phenotype from genotype. However, even for simple bacteria, whole-cell models will contain thousands of parameters, many of which are poorly characterized or unknown. New algorithms are needed to estimate these parameters and enable researchers to build increasingly comprehensive models. We organized the Dialogue for Reverse Engineering Assessments and Methods (DREAM) 8 Whole-Cell Parameter Estimation Challenge to develop new parameter estimation algorithms for whole-cell models. We asked participants to identify a subset of parameters of a whole-cell model given the model’s structure and in silico “experimental” data. Here we describe the challenge, the best performing methods, and new insights into the identifiability of whole-cell models. We also describe several valuable lessons we learned toward improving future challenges. Going forward, we believe that collaborative efforts supported by inexpensive cloud computing have the potential to solve whole-cell model parameter estimation.

Whole-cell models promise to enable rational bioengineering by predicting how cells behave. Even for simple bacteria, whole-cell models require thousands of parameters, many of which are poorly characterized or unknown. New approaches are needed to estimate these parameters. We organized the Dialogue for Reverse Engineering Assessments and Methods (DREAM) 8 Whole-Cell Parameter Estimation Challenge to develop new approaches for whole-cell model parameter identification. Here we describe the challenge, the best performing methods, new insights into the identifiability of whole-cell models, and several lessons we learned for improving future challenges. Going forward, we believe that collaborative efforts have the potential to produce powerful tools for identifying whole-cell models.

Quantifying the underlying landscape and paths of cancer

Cancer is a disease regulated by the underlying gene networks. The emergence of normal and cancer states as well as the transformation between them can be thought of as a result of the gene network interactions and associated changes. We developed a global potential landscape and path framework to quantify cancer and associated processes. We constructed a cancer gene regulatory network based on the experimental evidences and uncovered the underlying landscape. The resulting tristable landscape characterizes important biological states: normal, cancer and apoptosis. The landscape topography in terms of barrier heights between stable state attractors quantifies the global stability of the cancer network system. We propose two mechanisms of cancerization: one is by the changes of landscape topography through the changes in regulation strengths of the gene networks. The other is by the fluctuations that help the system to go over the critical barrier at fixed landscape topography. The kinetic paths from least action principle quantify the transition processes among normal state, cancer state and apoptosis state. The kinetic rates provide the quantification of transition speeds among normal, cancer and apoptosis attractors. By the global sensitivity analysis of the gene network parameters on the landscape topography, we uncovered some key gene regulations determining the transitions between cancer and normal states. This can be used to guide the design of new anti-cancer tactics, through cocktail strategy of targeting multiple key regulation links simultaneously, for preventing cancer occurrence or transforming the early cancer state back to normal state.

Modeling the epigenetic attractors landscape: toward a post-genomic mechanistic understanding of development

Robust temporal and spatial patterns of cell types emerge in the course of normal development in multicellular organisms. The onset of degenerative diseases may result from altered cell fate decisions that give rise to pathological phenotypes. Complex networks of genetic and non-genetic components underlie such normal and altered morphogenetic patterns. Here we focus on the networks of regulatory interactions involved in cell-fate decisions. Such networks modeled as dynamical non-linear systems attain particular stable configurations on gene activity that have been interpreted as cell-fate states. The network structure also restricts the most probable transition patterns among such states. The so-called Epigenetic Landscape (EL), originally proposed by C. H. Waddington, was an early attempt to conceptually explain the emergence of developmental choices as the result of intrinsic constraints (regulatory interactions) shaped during evolution. Thanks to the wealth of molecular genetic and genomic studies, we are now able to postulate gene regulatory networks (GRN) grounded on experimental data, and to derive EL models for specific cases. This, in turn, has motivated several mathematical and computational modeling approaches inspired by the EL concept, that may be useful tools to understand and predict cell-fate decisions and emerging patterns. In order to distinguish between the classical metaphorical EL proposal of Waddington, we refer to the Epigenetic Attractors Landscape (EAL), a proposal that is formally framed in the context of GRNs and dynamical systems theory. In this review we discuss recent EAL modeling strategies, their conceptual basis and their application in studying the emergence of both normal and pathological developmental processes. In addition, we discuss how model predictions can shed light into rational strategies for cell fate regulation, and we point to challenges ahead.

Implementation of proteomics for cancer research: past, present, and future

Cancer is the leading cause of the death, accounts for about 13% of all annual deaths worldwide. Many different fields of science are collaborating together studying cancer to improve our knowledge of this lethal disease, and find better solutions for diagnosis and treatment. Proteomics is one of the most recent and rapidly growing areas in molecular biology that helps understanding cancer from an omics data analysis point of view. The human proteome project was officially initiated in 2008. Proteomics enables the scientists to interrogate a variety of biospecimens for their protein contents and measure the concentrations of these proteins. Current necessary equipment and technologies for cancer proteomics are mass spectrometry, protein microarrays, nanotechnology and bioinformatics. In this paper, we provide a brief review on proteomics and its application in cancer research. After a brief introduction including its definition, we summarize the history of major previous work conducted by researchers, followed by an overview on the role of proteomics in cancer studies. We also provide a list of different utilities in cancer proteomics and investigate their advantages and shortcomings from theoretical and practical angles. Finally, we explore some of the main challenges and conclude the paper with future directions in this field.

Network-Based Drug-Target Interaction Prediction with Probabilistic Soft Logic

Drug-target interaction studies are important because they can predict drugs' unexpected therapeutic or adverse side effects. In silico predictions of potential interactions are valuable and can focus effort on in vitro experiments. We propose a prediction framework that represents the problem using a bipartite graph of drug-target interactions augmented with drug-drug and target-target similarity measures and makes predictions using probabilistic soft logic (PSL). Using probabilistic rules in PSL, we predict interactions with models based on triad and tetrad structures. We apply (blocking) techniques that make link prediction in PSL more efficient for drug-target interaction prediction. We then perform extensive experimental studies to highlight different aspects of the model and the domain, first comparing the models with different structures and then measuring the effect of the proposed blocking on the prediction performance and efficiency. We demonstrate the importance of rule weight learning in the proposed PSL model and then show that PSL can effectively make use of a variety of similarity measures. We perform an experiment to validate the importance of collective inference and using multiple similarity measures for accurate predictions in contrast to non-collective and single similarity assumptions. Finally, we illustrate that our PSL model achieves state-of-the-art performance with simple, interpretable rules and evaluate our novel predictions using online data sets.

Stem cell differentiation as a many-body problem

Stem cell differentiation has been viewed as coming from transitions between attractors on an epigenetic landscape that governs the dynamics of a regulatory network involving many genes. Rigorous definition of such a landscape is made possible by the realization that gene regulation is stochastic, owing to the small copy number of the transcription factors that regulate gene expression and because of the single-molecule nature of the gene itself. We develop an approximation that allows the quantitative construction of the epigenetic landscape for large realistic model networks. Applying this approach to the network for embryonic stem cell development explains many experimental observations, including the heterogeneous distribution of the transcription factor Nanog and its role in safeguarding the stem cell pluripotency, which can be understood by finding stable steady-state attractors and the most probable transition paths between those attractors. We also demonstrate that the switching rate between attractors can be significantly influenced by the gene expression noise arising from the fluctuations of DNA occupancy when binding to a specific DNA site is slow.

A general strategy for cellular reprogramming: the importance of transcription factor cross-repression

Transcription factor cross-repression is an important concept in cellular differentiation. A bistable toggle switch constitutes a molecular mechanism that determines cellular commitment and provides stability to transcriptional programs of binary cell fate choices. Experiments support that perturbations of these toggle switches can interconvert these binary cell fate choices, suggesting potential reprogramming strategies. However, more complex types of cellular transitions could involve perturbations of combinations of different types of multistable motifs. Here, we introduce a method that generalizes the concept of transcription factor cross-repression to systematically predict sets of genes, whose perturbations induce cellular transitions between any given pair of cell types. Furthermore, to our knowledge, this is the first method that systematically makes these predictions without prior knowledge of potential candidate genes and pathways involved, providing guidance on systems where little is known. Given the increasing interest of cellular reprogramming in medicine and basic research, our method represents a useful computational methodology to assist researchers in the field in designing experimental strategies.

Epigenetic state network approach for describing cell phenotypic transitions

Recent breakthroughs of cell phenotype reprogramming impose theoretical challenges on unravelling the complexity of large circuits maintaining cell phenotypes coupled at many different epigenetic and gene regulation levels, and quantitatively describing the phenotypic transition dynamics. A popular picture proposed by Waddington views cell differentiation as a ball sliding down a landscape with valleys corresponding to different cell types separated by ridges. Based on theories of dynamical systems, we establish a novel 'epigenetic state network' framework that captures the global architecture of cell phenotypes, which allows us to translate the metaphorical low-dimensional Waddington epigenetic landscape concept into a simple-yet-predictive rigorous mathematical framework of cell phenotypic transitions. Specifically, we simplify a high-dimensional epigenetic landscape into a collection of discrete states corresponding to stable cell phenotypes connected by optimal transition pathways among them. We then apply the approach to the phenotypic transition processes among fibroblasts (FBs), pluripotent stem cells (PSCs) and cardiomyocytes (CMs). The epigenetic state network for this case predicts three major transition pathways connecting FBs and CMs. One goes by way of PSCs. The other two pathways involve transdifferentiation either indirectly through cardiac progenitor cells or directly from FB to CM. The predicted pathways and multiple intermediate states are supported by existing microarray data and other experiments. Our approach provides a theoretical framework for studying cell phenotypic transitions. Future studies at single-cell levels can directly test the model predictions.

Detecting cellular reprogramming determinants by differential stability analysis of gene regulatory networks

BACKGROUND

Cellular differentiation and reprogramming are processes that are carefully orchestrated by the activation and repression of specific sets of genes. An increasing amount of experimental results show that despite the large number of genes participating in transcriptional programs of cellular phenotypes, only few key genes, which are coined here as reprogramming determinants, are required to be directly perturbed in order to induce cellular reprogramming. However, identification of reprogramming determinants still remains a combinatorial problem, and the state-of-art methods addressing this issue rests on exhaustive experimentation or prior knowledge to narrow down the list of candidates.

RESULTS

Here we present a computational method, without any preliminary selection of candidate genes, to identify reduced subsets of genes, which when perturbed can induce transitions between cellular phenotypes. The method relies on the expression profiles of two stable cellular phenotypes along with a topological analysis stability elements in the gene regulatory network that are necessary to cause this multi-stability. Since stable cellular phenotypes can be considered as attractors of gene regulatory networks, cell fate and cellular reprogramming involves transition between these attractors, and therefore current method searches for combinations of genes that are able to destabilize a specific initial attractor and stabilize the final one in response to the appropriate perturbations.

CONCLUSIONS

The method presented here represents a useful framework to assist researchers in the field of cellular reprogramming to design experimental strategies with potential applications in the regenerative medicine and disease modelling.

Quantifying Waddington landscapes and paths of non-adiabatic cell fate decisions for differentiation, reprogramming and transdifferentiation

Cellular differentiation, reprogramming and transdifferentiation are determined by underlying gene regulatory networks. Non-adiabatic regulation via slow binding/unbinding to the gene can be important in these cell fate decision-making processes. Based on a stem cell core gene network, we uncovered the stem cell developmental landscape. As the binding/unbinding speed decreases, the landscape topography changes from bistable attractors of stem and differentiated states to more attractors of stem and other different cell states as well as substates. Non-adiabaticity leads to more differentiated cell types and provides a natural explanation for the heterogeneity observed in the experiments. We quantified Waddington landscapes with two possible cell fate decision mechanisms by changing the regulation strength or regulation timescale (non-adiabaticity). Transition rates correlate with landscape topography through barrier heights between different states and quantitatively determine global stability. We found the optimal speeds of these cell fate decision-making processes. We quantified biological paths and predict that differentiation and reprogramming go through an intermediate state (IM1), whereas transdifferentiation goes through another intermediate state (IM2). Some predictions are confirmed by recent experimental studies.

Quantifying cell fate decisions for differentiation and reprogramming of a human stem cell network: landscape and biological paths

Cellular reprogramming has been recently intensively studied experimentally. We developed a global potential landscape and kinetic path framework to explore a human stem cell developmental network composed of 52 genes. We uncovered the underlying landscape for the stem cell network with two basins of attractions representing stem and differentiated cell states, quantified and exhibited the high dimensional biological paths for the differentiation and reprogramming process, connecting the stem cell state and differentiated cell state. Both the landscape and non-equilibrium curl flux determine the dynamics of cell differentiation jointly. Flux leads the kinetic paths to be deviated from the steepest descent gradient path, and the corresponding differentiation and reprogramming paths are irreversible. Quantification of paths allows us to find out how the differentiation and reprogramming occur and which important states they go through. We show the developmental process proceeds as moving from the stem cell basin of attraction to the differentiation basin of attraction. The landscape topography characterized by the barrier heights and transition rates quantitatively determine the global stability and kinetic speed of cell fate decision process for development. Through the global sensitivity analysis, we provided some specific predictions for the effects of key genes and regulation connections on the cellular differentiation or reprogramming process. Key links from sensitivity analysis and biological paths can be used to guide the differentiation designs or reprogramming tactics.

Cellular differentiation and reprogramming have been extensively studied using experimental methods. We developed a landscape and kinetic path approach to explore the global stability of a stem cell developmental network. The cell fates are quantified by the basins of attractions of the underlying landscape. The developmental process can be quantitatively described and uncovered by the biological paths on the landscape from the progenitor state to the differentiation state. This allows us to trace the underlying detailed kinetic process and obtain the recipe for engineering differentiation and reprogramming. By quantifying the landscape topography by the barrier heights and dynamic transition speed, we can evaluate the stability and kinetics of cell fate decision making process of the development and reprogramming. The global sensitivity analysis provided predictions about the effects of the key genes and regulation links of the network on the stability of differentiation and reprogramming process. This can be tested in the experiments. Results from sensitivity analysis and biological paths acquired can be used to guide the differentiation designs or reprogramming tactics.

Modeling to optimize terminal stem cell differentiation

Embryonic stem cell (ESC), iPCs, and adult stem cells (ASCs) all are among the most promising potential treatments for heart failure, spinal cord injury, neurodegenerative diseases, and diabetes. However, considerable uncertainty in the production of ESC-derived terminally differentiated cell types has limited the efficiency of their development. To address this uncertainty, we and other investigators have begun to employ a comprehensive statistical model of ESC differentiation for determining the role of intracellular pathways (e.g., STAT3) in ESC differentiation and determination of germ layer fate. The approach discussed here applies the Baysian statistical model to cell/developmental biology combining traditional flow cytometry methodology and specific morphological observations with advanced statistical and probabilistic modeling and experimental design. The final result of this study is a unique tool and model that enhances the understanding of how and when specific cell fates are determined during differentiation. This model provides a guideline for increasing the production efficiency of therapeutically viable ESCs/iPSCs/ASC derived neurons or any other cell type and will eventually lead to advances in stem cell therapy.

Therapeutic hints from analyzing the attractor landscape of the p53 regulatory circuit

Genes are interconnected in the cell to form a genetic network that regulates cell fate. Targeting multiple genes is expected to be more effective in developing therapeutics than targeting single genes. A recent study demonstrated the possibility of systematically searching for such combinatorial treatments by characterizing the attractor landscape of the p53 regulatory circuit.

Expression of pluripotency factors in larval epithelia of the frog Xenopus: evidence for the presence of cornea epithelial stem cells

Understanding the biology of somatic stem cells in self renewing tissues represents an exciting field of study, especially given the potential to harness these cells for tissue regeneration and repair in treating injury and disease. The mammalian cornea contains a population of basal epithelial stem cells involved in cornea homeostasis and repair. Research has been restricted to mammalian systems and little is known about the presence or function of these stem cells in other vertebrates. Therefore, we carried out studies to characterize frog cornea epithelium. Careful examination shows that the Xenopus larval cornea epithelium consists of three distinct layers that include an outer epithelial layer and underlying basal epithelium, in addition to a deeper fibrous layer that contains the main sensory nerve trunks that give rise to numerous branches that extend into these epithelia. These nerves convey sensory and presumably also autonomic innervation to those tissues. The sensory nerves are all derived as branches of the trigeminal nerve/ganglion similar to the situation encountered in mammals, though there appear to be some potentially interesting differences, which are detailed in this paper. We show further that numerous pluripotency genes are expressed by cells in the cornea epithelium, including: sox2, p63, various oct4 homologs, c-myc, klf4 and many others. Antibody localization revealed that p63, a well known mammalian epithelial stem cell marker, was localized strictly to all cells in the basal cornea epithelium. c-myc, was visualized in a smaller subset of basal epithelial cells and adjacent stromal tissue predominately at the periphery of the cornea (limbal zone). Finally, sox2 protein was found to be present throughout all cells of both the outer and basal epithelia, but was much more intensely expressed in a distinct subset of cells that appeared to be either multinucleate or possessed multi-lobed nuclei that are normally located at the periphery of the cornea. Using a thymidine analog (EdU), we were able to label mitotically active cells, which revealed that cell proliferation takes place throughout the cornea epithelium, predominantly in the basal epithelial layer. Species of Xenopus and one other amphibian are unique in their ability to replace a missing lens from cells derived from the basal cornea epithelium. Using EdU we show, as others have previously, that proliferating cells within the cornea epithelium do contribute to the formation of these regenerated lenses. Furthermore, using qPCR we determined that representatives of various pluripotency genes (i.e., sox2, p63 and oct60) are upregulated early during the process of lens regeneration. Antibody labeling showed that the number of sox2 expressing cells increased dramatically within 4 h following lens removal and these cells were scattered throughout the basal layer of the cornea epithelium. Historically, the process of lens regeneration in Xenopus had been described as one involving transdifferentiation of cornea epithelial cells (i.e., one involving cellular dedifferentiation followed by redifferentiation). Our combined observations provide evidence that a population of stem cells exists within the Xenopus cornea. We hypothesize that the basal epithelium contains oligopotent epithelial stem cells that also represent the source of regenerated lenses in the frog. Future studies will be required to clearly identify the source of these lenses.

Feeder-free derivation of human induced pluripotent stem cells with messenger RNA

The therapeutic promise of induced pluripotent stem cells (iPSCs) has spurred efforts to circumvent genome alteration when reprogramming somatic cells to pluripotency. Approaches based on episomal DNA, Sendai virus, and messenger RNA (mRNA) can generate "footprint-free" iPSCs with efficiencies equaling or surpassing those attained with integrating viral vectors. The mRNA method uniquely affords unprecedented control over reprogramming factor (RF) expression while obviating a cleanup phase to purge residual traces of vector. Currently, mRNA-based reprogramming is relatively laborious due to the need to transfect daily for ~2 weeks to induce pluripotency, and requires the use of feeder cells that add complexity and variability to the procedure while introducing a route for contamination with non-human-derived biological material. We accelerated the mRNA reprogramming process through stepwise optimization of the RF cocktail and leveraged these kinetic gains to establish a feeder-free, xeno-free protocol which slashes the time, cost and effort involved in iPSC derivation.

A stochastic model of epigenetic dynamics in somatic cell reprogramming

Somatic cell reprogramming has dramatically changed stem cell research in recent years. The high pace of new findings in the field and an ever increasing amount of data from new high throughput techniques make it challenging to isolate core principles of the process. In order to analyze such mechanisms, we developed an abstract mechanistic model of a subset of the known regulatory processes during cell differentiation and production of induced pluripotent stem cells. This probabilistic Boolean network describes the interplay between gene expression, chromatin modifications, and DNA methylation. The model incorporates recent findings in epigenetics and partially reproduces experimentally observed reprogramming efficiencies and changes in methylation and chromatin remodeling. It enables us to investigate, how the temporal progression of the process is regulated. It also explicitly includes the transduction of factors using viral vectors and their silencing in reprogrammed cells, since this is still a standard procedure in somatic cell reprogramming. Based on the model we calculate an epigenetic landscape for probabilities of cell states. Simulation results show good reproduction of experimental observations during reprogramming, despite the simple structure of the model. An extensive analysis and introduced variations hint toward possible optimizations of the process that could push the technique closer to clinical applications. Faster changes in DNA methylation increase the speed of reprogramming at the expense of efficiency, while accelerated chromatin modifications moderately improve efficiency.