Chiara Damiani

Dynamical properties of a Boolean model of gene regulatory network with memory

Classical random Boolean networks (RBN) are not well suited to describe experimental data from time-course microarray, mainly because of the strict assumptions about the synchronicity of the regulatory mechanisms. In order to overcome this setback, a generalization of the RBN model is described and analyzed. Gene products (e.g., regulatory proteins) are introduced, with each one characterized by a specific decay time, thereby introducing a form of memory in the system. The dynamics of these networks is analyzed, and it is shown that the distribution of the decay times has a strong effect that can be adequately described and understood. The implications for the dynamical criticality of the networks are also discussed.

Computational Strategies for a System-Level Understanding of Metabolism

Cell metabolism is the biochemical machinery that provides energy and building blocks to sustain life. Understanding its fine regulation is of pivotal relevance in several fields, from metabolic engineering applications to the treatment of metabolic disorders and cancer. Sophisticated computational approaches are needed to unravel the complexity of metabolism. To this aim, a plethora of methods have been developed, yet it is generally hard to identify which computational strategy is most suited for the investigation of a specific aspect of metabolism. This review provides an up-to-date description of the computational methods available for the analysis of metabolic pathways, discussing their main advantages and drawbacks. In particular, attention is devoted to the identification of the appropriate scale and level of accuracy in the reconstruction of metabolic networks, and to the inference of model structure and parameters, especially when dealing with a shortage of experimental measurements. The choice of the proper computational methods to derive in silico data is then addressed, including topological analyses, constraint-based modeling and simulation of the system dynamics. A description of some computational approaches to gain new biological knowledge or to formulate hypotheses is finally provided.

Cell–cell interaction and diversity of emergent behaviours

Despite myriads of possible gene expression profiles, cells tend to be found in a confined number of expression patterns. The dynamics of Boolean models of gene regulatory networks has proven to be a likely candidate for the description of such self-organisation phenomena. Because cells do not live in isolation, but they constantly shape their functions to adapt to signals from other cells, this raises the question of whether the cooperation among cells entails an expansion or a reduction of their possible steady states. Multi random Boolean networks are introduced here as a model for interaction among cells that might be suitable for the investigation of some generic properties regarding the influence of communication on the diversity of cell behaviours. In spite of its simplicity, the model exhibits a non-obvious phenomenon according to which a moderate exchange of products among adjacent cells fosters the variety of their possible behaviours, which on the other hand are more similar to one another. On the contrary, a more invasive coupling would lead cells towards homogeneity.

Information transfer among coupled random boolean networks

Information processing and information flow occur at many levels in the course of an organisms development and throughout its lifespan. Biological networks inside cells transmit information from their inputs (e.g. the concentrations of proteins or other signaling molecules) to their outputs (e.g. the expression levels of various genes). Moreover, cells do not exist in isolation, but they constantly interact with one another. We study the information flow in a model of interacting genetic networks, which are represented as Boolean graphs. It is observed that the information transfer among the networks is not linearly dependent on the amount of nodes that are able to influence the state of genes in surrounding cells.

Zooming-in on cancer metabolic rewiring with tissue specific constraint-based models

The metabolic rearrangements occurring in cancer cells can be effectively investigated with a Systems Biology approach supported by metabolic network modeling. We here present tissue-specific constraint-based core models for three different types of tumors (liver, breast and lung) that serve this purpose. The core models were extracted and manually curated from the corresponding genome-scale metabolic models in the Human Metabolic Atlas database with a focus on the pathways that are known to play a key role in cancer growth and proliferation. Along similar lines, we also reconstructed a core model from the original general human metabolic network to be used as a reference model. A comparative Flux Balance Analysis between the reference and the cancer models highlighted both a clear distinction between the two conditions and a heterogeneity within the three different cancer types in terms of metabolic flux distribution. These results emphasize the need for modeling approaches able to keep up with this tumoral heterogeneity in order to identify more suitable drug targets and develop effective treatments. According to this perspective, we identified key points able to reverse the tumoral phenotype toward the reference one or vice-versa.

Growth and Division in a Dynamic Protocell Model

In this paper a new model of growing and dividing protocells is described, whose main features are (i) a lipid container that grows according to the composition of the molecular milieu (ii) a set of “genetic memory molecules” (GMMs) that undergo catalytic reactions in the internal aqueous phase and (iii) a set of stochastic kinetic equations for the GMMs. The mass exchange between the external environment and the internal phase is described by simulating a semipermeable membrane and a flow driven by the differences in chemical potentials, thereby avoiding to resort to sometimes misleading simplifications, e.g., that of a flow reactor. Under simple assumptions, it is shown that synchronization takes place between the rate of replication of the GMMs and that of the container, provided that the set of reactions hosts a so-called RAF (Reflexive Autocatalytic, Food-generated) set whose influence on synchronization is hereafter discussed. It is also shown that a slight modification of the basic model that takes into account a rate-limiting term, makes possible the growth of novelties, allowing in such a way suitable evolution: so the model represents an effective basis for understanding the main abstract properties of populations of protocells.

CABeRNET: a Cytoscape app for augmented Boolean models of gene regulatory NETworks

BackgroundDynamical models of gene regulatory networks (GRNs) are highly effective in describing complex biological phenomena and processes, such as cell differentiation and cancer development. Yet, the topological and functional characterization of real GRNs is often still partial and an exhaustive picture of their functioning is missing.ResultsWe here introduce CABeRNET, a Cytoscape app for the generation, simulation and analysis of Boolean models of GRNs, specifically focused on their augmentation when a only partial topological and functional characterization of the network is available. By generating large ensembles of networks in which user-defined entities and relations are added to the original core, CABeRNET allows to formulate hypotheses on the missing portions of real networks, as well to investigate their generic properties, in the spirit of complexity science.ConclusionsCABeRNET offers a series of innovative simulation and modeling functions and tools, including (but not being limited to) the dynamical characterization of the gene activation patterns ruling cell types and differentiation fates, and sophisticated robustness assessments, as in the case of gene knockouts. The integration within the widely used Cytoscape framework for the visualization and analysis of biological networks, makes CABeRNET a new essential instrument for both the bioinformatician and the computational biologist, as well as a computational support for the experimentalist. An example application concerning the analysis of an augmented T-helper cell GRN is provided.

An ensemble evolutionary constraint-based approach to understand the emergence of metabolic phenotypes

Constraint-based modeling is largely used in computational studies of metabolism. We propose here a novel approach that aims to identify ensembles of flux distributions that comply with one or more target phenotype(s). The methodology has been tested on a small-scale model of yeast energy metabolism. The target phenotypes describe the differential pattern of ethanol production and O2 consumption observed in “Crabtree-positive” and “Crabtree-negative” yeasts in changing environment (i.e., when the upper limit of glucose uptake is varied). The ensembles were obtained either by selection among sampled flux distributions or by means of a search heuristic (genetic algorithm). The former approach provided indication about the probability to observe a given phenotype, but the resulting ensembles could not be unambiguously partitioned into “Crabtree-positive” and “Crabtree-negative” clusters. On the contrary well-separated clusters were obtained with the latter method. The cluster analysis further allowed identification of distinct groups within each target phenotype. The method may thus prove useful in characterizing the design principles underlying metabolic plasticity arising from evolving physio-pathological or developmental constraints.

The Diffusion of Perturbations in a Model of Coupled Random Boolean Networks

Deciphering the influence of the interaction among the constituents of a complex system on the overall behaviour is one of the main goals of complex systems science. The model we present in this work is a 2D square cellular automaton whose of each cell is occupied by a complete random Boolean network. Random Boolean networks are a well-known simplified model of genetic regulatory networks and this model of interacting RBNs may be therefore regarded as a simplified model of a tissue or a monoclonal colony. The mechanism of cell-to-cell interaction is here simulated letting some nodes of a particular network being influenced by the state of some nodes belonging to its neighbouring cells. One possible means to investigate the overall dynamics of a complex system is studying its response to perturbations. Our analyses follow this methodological approach. Even though the dynamics of the system is far from trivial we could show in a clear way how the interaction affects the dynamics and the global degree of order.

A stochastic model of catalytic reaction networks in protocells

Protocells are supposed to have played a key role in the self-organizing processes leading to the emergence of life. Existing models either (i) describe protocell architecture and dynamics, given the existence of sets of collectively self-replicating molecules for granted, or (ii) describe the emergence of the aforementioned sets from an ensemble of random molecules in a simple experimental setting (e.g. a closed system or a steady-state flow reactor) that does not properly describe a protocell. In this paper we present a model that goes beyond these limitations by describing the dynamics of sets of replicating molecules within a lipid vesicle. We adopt the simplest possible protocell architecture, by considering a semi-permeable membrane that selects the molecular types that are allowed to enter or exit the protocell and by assuming that the reactions take place in the aqueous phase in the internal compartment. As a first approximation, we ignore the protocell growth and division dynamics. The behavior of catalytic reaction networks is then simulated by means of a stochastic model that accounts for the creation and the extinction of species and reactions. While this is not yet an exhaustive protocell model, it already provides clues regarding some processes that are relevant for understanding the conditions that can enable a population of protocells to undergo evolution and selection.