Scale invariance in biology: coincidence or footprint of a universal mechanism?

In this article, we present a self-contained review of recent work on complex biological systems which exhibit no characteristic scale. This property can manifest itself with fractals (spatial scale invariance), flicker noise or 1/f-noise where f denotes the frequency of a signal (temporal scale invariance) and power laws (scale invariance in the size and duration of events in the dynamics of the system). A hypothesis recently put forward to explain these scale-free phenomomena is criticality, a notion introduced by physicists while studying phase transitions in materials, where systems spontaneously arrange themselves in an unstable manner similar, for instance, to a row of dominoes. Here, we review in a critical manner work which investigates to what extent this idea can be generalized to biology. More precisely, we start with a brief introduction to the concepts of absence of characteristic scale (power-law distributions, fractals and 1/f-noise) and of critical phenomena. We then review typical mathematical models exhibiting such properties: edge of chaos, cellular automata and self-organized critical models. These notions are then brought together to see to what extent they can account for the scale invariance observed in ecology, evolution of species, type III epidemics and some aspects of the central nervous system. This article also discusses how the notion of scale invariance can give important insights into the workings of biological systems.     

Cellular Sensory Mechanisms for Detecting Specific Fold-Changes in Extracellular Cues

Cellular sensory systems often respond not to the absolute levels of inputs but to the fold-changes in inputs. Such a property is called fold-change detection (FCD) and is important for accurately sensing dynamic changes in environmental signals in the presence of fluctuations in their absolute levels. Previous studies defined FCD as input-scale invariance and proposed several biochemical models that achieve such a condition. Here, we prove that the previous FCD models can be approximated by a log-differentiator. Although the log-differentiator satisfies the input-scale invariance requirement, its response amplitude and response duration strongly depend on the input timescale. This creates limitations in the specificity and repeatability of detecting fold-changes in inputs. Nevertheless, FCD with specificity and repeatability by cells has been reported in the context of Drosophila wing development. Motivated by this fact and by extending previous FCD models, we here propose two possible mechanisms to achieve FCD with specificity and repeatability. One is the integrate-and-fire type: a system integrates the rate of temporal change in input and makes a response when the integrated value reaches a constant threshold, and this is followed by the reset of the integrated value. The other is the dynamic threshold type: a system response occurs when the input level reaches a threshold, whose value is multiplied by a certain constant after each response. These two mechanisms can be implemented biochemically by appropriately combining feed-forward and feedback loops. The main difference between the two models is their memory of input history; we discuss possible ways to distinguish between the two models experimentally.     

Sloppiness, robustness, and evolvability in systems biology

The functioning of many biochemical networks is often robust-remarkably stable under changes in external conditions and internal reaction parameters. Much recent work on robustness and evolvability has focused on the structure of neutral spaces, in which system behavior remains invariant to mutations. Recently we have shown that the collective behavior of multiparameter models is most often sloppy: insensitive to changes except along a few 'stiff' combinations of parameters, with an enormous sloppy neutral subspace. Robustness is often assumed to be an emergent evolved property, but the sloppiness natural to biochemical networks offers an alternative nonadaptive explanation. Conversely, ideas developed to study evolvability in robust systems can be usefully extended to characterize sloppy systems.     

Identification of a Topological Characteristic Responsible for the Biological Robustness of Regulatory Networks

Attribution of biological robustness to the specific structural properties of a regulatory network is an important yet unsolved problem in systems biology. It is widely believed that the topological characteristics of a biological control network largely determine its dynamic behavior, yet the actual mechanism is still poorly understood. Here, we define a novel structural feature of biological networks, termed ‘regulation entropy’, to quantitatively assess the influence of network topology on the robustness of the systems. Using the cell-cycle control networks of the budding yeast (Saccharomyces cerevisiae) and the fission yeast (Schizosaccharomyces pombe) as examples, we first demonstrate the correlation of this quantity with the dynamic stability of biological control networks, and then we establish a significant association between this quantity and the structural stability of the networks. And we further substantiate the generality of this approach with a broad spectrum of biological and random networks. We conclude that the regulation entropy is an effective order parameter in evaluating the robustness of biological control networks. Our work suggests a novel connection between the topological feature and the dynamic property of biological regulatory networks.     

Cellular Sensory Mechanisms for Detecting Specific Fold-Changes in Extracellular Cues

Cellular sensory systems often respond not to the absolute levels of inputs but to the fold-changes in inputs. Such a property is called fold-change detection (FCD) and is important for accurately sensing dynamic changes in environmental signals in the presence of fluctuations in their absolute levels. Previous studies defined FCD as input-scale invariance and proposed several biochemical models that achieve such a condition. Here, we prove that the previous FCD models can be approximated by a log-differentiator. Although the log-differentiator satisfies the input-scale invariance requirement, its response amplitude and response duration strongly depend on the input timescale. This creates limitations in the specificity and repeatability of detecting fold-changes in inputs. Nevertheless, FCD with specificity and repeatability by cells has been reported in the context of Drosophila wing development. Motivated by this fact and by extending previous FCD models, we here propose two possible mechanisms to achieve FCD with specificity and repeatability. One is the integrate-and-fire type: a system integrates the rate of temporal change in input and makes a response when the integrated value reaches a constant threshold, and this is followed by the reset of the integrated value. The other is the dynamic threshold type: a system response occurs when the input level reaches a threshold, whose value is multiplied by a certain constant after each response. These two mechanisms can be implemented biochemically by appropriately combining feed-forward and feedback loops. The main difference between the two models is their memory of input history; we discuss possible ways to distinguish between the two models experimentally.     

Reduced models of networks of coupled enzymatic reactions

The Michaelis-Menten equation has played a central role in our understanding of biochemical processes. It has long been understood how this equation approximates the dynamics of irreversible enzymatic reactions. However, a similar approximation in the case of networks, where the product of one reaction can act as an enzyme in another, has not been fully developed. Here we rigorously derive such an approximation in a class of coupled enzymatic networks where the individual interactions are of Michaelis-Menten type. We show that the sufficient conditions for the validity of the total quasi-steady state assumption (tQSSA), obtained in a single protein case by Borghans, de Boer and Segel can be extended to sufficient conditions for the validity of the tQSSA in a large class of enzymatic networks. Secondly, we derive reduced equations that approximate the network's dynamics and involve only protein concentrations. This significantly reduces the number of equations necessary to model such systems. We prove the validity of this approximation using geometric singular perturbation theory and results about matrix differentiation. The ideas used in deriving the approximating equations are quite general, and can be used to systematize other model reductions.     

Pleiotropic and Epistatic Behavior of a Ring-Hydroxylating Oxygenase System in the Polycyclic Aromatic Hydrocarbon Metabolic Network from Mycobacterium vanbaalenii PYR-1

Despite the considerable knowledge of bacterial high-molecular-weight (HMW) polycyclic aromatic hydrocarbon (PAH) metabolism, the key enzyme(s) and its pleiotropic and epistatic behavior(s) responsible for low-molecular-weight (LMW) PAHs in HMW PAH-metabolic networks remain poorly understood. In this study, a phenotype-based strategy, coupled with a spray plate method, selected a Mycobacterium vanbaalenii PYR-1 mutant (6G11) that degrades HMW PAHs but not LMW PAHs. Sequence analysis determined that the mutant was defective in pdoA2, encoding an aromatic ring-hydroxylating oxygenase (RHO). A series of metabolic comparisons using high-performance liquid chromatography (HPLC) analysis revealed that the mutant had a lower rate of degradation of fluorene, anthracene, and pyrene. Unlike the wild type, the mutant did not produce a color change in culture media containing fluorene, phenanthrene, and fluoranthene. An Escherichia coli expression experiment confirmed the ability of the Pdo system to oxidize biphenyl, the LMW PAHs naphthalene, phenanthrene, anthracene, and fluorene, and the HMW PAHs pyrene, fluoranthene, and benzo[a]pyrene, with the highest enzymatic activity directed toward three-ring PAHs. Structure analysis and PAH substrate docking simulations of the Pdo substrate-binding pocket rationalized the experimentally observed metabolic versatility on a molecular scale. Using information obtained in this study and from previous work, we constructed an RHO-centric functional map, allowing pleiotropic and epistatic enzymatic explanation of PAH metabolism. Taking the findings together, the Pdo system is an RHO system with the pleiotropic responsibility of LMW PAH-centric hydroxylation, and its epistatic functional contribution is also crucial for the metabolic quality and quantity of the PAH-MN.     

Stability of Metabolic Correlations under Changing Environmental Conditions in Escherichia coli – A Systems Approach

Background

Biological systems adapt to changing environments by reorganizing their cellular and physiological program with metabolites representing one important response level. Different stresses lead to both conserved and specific responses on the metabolite level which should be reflected in the underlying metabolic network.

Methodology/Principal Findings

Starting from experimental data obtained by a GC-MS based high-throughput metabolic profiling technology we here develop an approach that: (1) extracts network representations from metabolic condition-dependent data by using pairwise correlations, (2) determines the sets of stable and condition-dependent correlations based on a combination of statistical significance and homogeneity tests, and (3) can identify metabolites related to the stress response, which goes beyond simple observations about the changes of metabolic concentrations. The approach was tested with Escherichia coli as a model organism observed under four different environmental stress conditions (cold stress, heat stress, oxidative stress, lactose diauxie) and control unperturbed conditions. By constructing the stable network component, which displays a scale free topology and small-world characteristics, we demonstrated that: (1) metabolite hubs in this reconstructed correlation networks are significantly enriched for those contained in biochemical networks such as EcoCyc, (2) particular components of the stable network are enriched for functionally related biochemical pathways, and (3) independently of the response scale, based on their importance in the reorganization of the correlation network a set of metabolites can be identified which represent hypothetical candidates for adjusting to a stress-specific response.

Conclusions/Significance

Network-based tools allowed the identification of stress-dependent and general metabolic correlation networks. This correlation-network-based approach does not rely on major changes in concentration to identify metabolites important for stress adaptation, but rather on the changes in network properties with respect to metabolites. This should represent a useful complementary technique in addition to more classical approaches.     

Dimerization and Bifunctionality Confer Robustness to the Isocitrate Dehydrogenase Regulatory System in Escherichia coli

An important goal of systems biology is to develop quantitative models that explain how specific molecular features give rise to systems-level properties. Metabolic and regulatory pathways that contain multifunctional proteins are especially interesting to study from this perspective because they have frequently been observed to exhibit robustness: the ability for a system to perform its proper function even as levels of its components change. In this study, we use extensive biochemical data and algebraic modeling to develop and analyze a model that shows how robust behavior arises in the isocitrate dehydrogenase (IDH) regulatory system of Escherichia coli, which was shown in 1985 to experimentally exhibit robustness. E. coli IDH is regulated by reversible phosphorylation catalyzed by the bifunctional isocitrate dehydrogenase kinase/phosphatase (IDHKP), and the level of IDH activity determines whether carbon flux is directed through the glyoxylate bypass (for growth on two-carbon substrates) or the full tricarboxylic acid cycle. Our model, which incorporates recent structural data on IDHKP, identifies several specific biochemical features of the system (including homodimerization of IDH and bifunctionality of IDHKP) that provide a potential explanation for robustness. Using algebraic techniques, we derive an invariant that summarizes the steady-state relationship between the phospho-forms of IDH. We use the invariant in combination with kinetic data on IDHKP to calculate IDH activity at a range of total IDH levels and find that our model predicts robustness. Our work unifies much of the known biochemistry of the IDH regulatory system into a single quantitative framework and highlights the importance of constructing biochemically realistic models in systems biology.     

Using multilayer network analysis to explore the temporal dynamics of collective behavior

Social organisms often show collective behaviors such as group foraging or movement.Collective behaviors can emerge from interactions between group members and may depend on the behavior of key individuals.When social interactions change over time,collective behaviors may change because these behaviors emerge from interactions among individuals.Despite the importance of,and growing interest in,the temporal dynamics of social interactions,it is not clear how to quantify changes in interactions over time or measure their stability.Furthermore,the temporal scale at which we should observe changes in social networks to detect biologically meaningful changes is not always apparent.Here we use multilayer network analysis to quantify temporal dynamics of social networks of the social spider Stegodyphus dumicola and determine how these dynamics relate to individual and group behaviors.We found that social interactions changed over time at a constant rate.Variation in both network structure and the identity of a keystone individual was not related to the mean or variance of the collective prey attack speed.Individuals that maintained a large and stable number of connections,despite changes in network structure,were the boldest individuals in the group.Therefore,social interactions and boldness are linked across time,but group collective behavior is not influenced by the stability of the social network.Our work demonstrates that dynamic social networks can be modeled in a multilayer framework.This approach may reveal biologically important temporal changes to social structure in other systems.     

Biochemical properties and expression profile of human prolyl dipeptidase DPP9

Dipetidyl peptidase 9 (DPP9) is a prolyl dipeptidase preferentially cleaving the peptide bond after the penultimate proline residue. The biological function of DPP9 is unknown. In this study, we have significantly improved the yield using Strep·Tactin^® purification system and characterized the biochemical property of DPP9. Moreover, the dimer interaction mode was investigated by introducing a mutation (F842A) at the dimer interface, which abolished the enzymatic activity without disrupting its quaternary structure. Furthermore, DPP9 was found ubiquitously expressed in fibroblasts, epithelial, and blood cells. Surprisingly, contrary to previous report, we found that the expression levels of DPP8 and DPP9 did not change upon the activation of the PBMC or Jurkat cells. These results indicate that the biochemical property of DPP9 is very similar to that of DPP8, its homologous protease. DPP9 and DPP8 are likely redundant proteins carrying out overlapping functions in vivo.     

Queueing up for enzymatic processing: correlated signaling through coupled degradation

High‐throughput technologies have led to the generation of complex wiring diagrams as a post‐sequencing paradigm for depicting the interactions between vast and diverse cellular species. While these diagrams are useful for analyzing biological systems on a large scale, a detailed understanding of the molecular mechanisms that underlie the observed network connections is critical for the further development of systems and synthetic biology. Here, we use queueing theory to investigate how ‘waiting lines’ can lead to correlations between protein ‘customers’ that are coupled solely through a downstream set of enzymatic ‘servers’. Using the E. coli ClpXP degradation machine as a model processing system, we observe significant cross‐talk between two networks that are indirectly coupled through a common set of processors. We further illustrate the implications of enzymatic queueing using a synthetic biology application, in which two independent synthetic networks demonstrate synchronized behavior when common ClpXP machinery is overburdened. Our results demonstrate that such post‐translational processes can lead to dynamic connections in cellular networks and may provide a mechanistic understanding of existing but currently inexplicable links.     

Exploring multiplicity conditions in enzymatic reaction networks

In this work, a novel algorithmic approach to detect multiplicity of steady states in enzymatic reaction networks is presented. The method exploits the structural properties of networks derived from the Chemical Reaction Network Theory. In first instance, the space of parameters is divided in different regions according to the qualitative behavior induced by the parameters in the long term dynamics of the network. Once the regions are identified, a condition for the appearance of multiplicities is checked in the different regions by solving a given optimization problem. In this way, the method allows the characterization of the whole parameter space of biochemical networks in terms of the appearance or not of multistability. The approach is illustrated through a well‐known case of enzymatic catalysis with substrate inhibition. © 2009 American Institute of Chemical Engineers Biotechnol. Prog., 2009     

From intracellular signaling to population oscillations: bridging size- and time-scales in collective behavior

Collective behavior in cellular populations is coordinated by biochemical signaling networks within individual cells. Connecting the dynamics of these intracellular networks to the population phenomena they control poses a considerable challenge because of network complexity and our limited knowledge of kinetic parameters. However, from physical systems, we know that behavioral changes in the individual constituents of a collectively behaving system occur in a limited number of well-defined classes, and these can be described using simple models. Here, we apply such an approach to the emergence of collective oscillations in cellular populations of the social amoeba Dictyostelium discoideum. Through direct tests of our model with quantitative in vivo measurements of single-cell and population signaling dynamics, we show how a simple model can effectively describe a complex molecular signaling network at multiple size and temporal scales. The model predicts novel noise-driven single-cell and population-level signaling phenomena that we then experimentally observe. Our results suggest that like physical systems, collective behavior in biology may be universal and described using simple mathematical models.     

Stability from Structure: Metabolic Networks Are Unlike Other Biological Networks

In recent work, attempts have been made to link the structure of biochemical networks to their complex dynamics. It was shown that structurally stable network motifs are enriched in such networks. In this work, we investigate to what extent these findings apply to metabolic networks. To this end, we extend a previously proposed method by changing the null model for determining motif enrichment, by using interaction types directly obtained from structural interaction matrices, by generating a distribution of partial derivatives of reaction rates and by simulating enzymatic regulation on metabolic networks. Our findings suggest that the conclusions drawn in previous work cannot be extended to metabolic networks, that is, structurally stable network motifs are not enriched in metabolic networks.     

Learning Contrast-Invariant Cancellation of Redundant Signals in Neural Systems

Cancellation of redundant information is a highly desirable feature of sensory systems, since it would potentially lead to a more efficient detection of novel information. However, biologically plausible mechanisms responsible for such selective cancellation, and especially those robust to realistic variations in the intensity of the redundant signals, are mostly unknown. In this work, we study, via in vivo experimental recordings and computational models, the behavior of a cerebellar-like circuit in the weakly electric fish which is known to perform cancellation of redundant stimuli. We experimentally observe contrast invariance in the cancellation of spatially and temporally redundant stimuli in such a system. Our model, which incorporates heterogeneously-delayed feedback, bursting dynamics and burst-induced STDP, is in agreement with our in vivo observations. In addition, the model gives insight on the activity of granule cells and parallel fibers involved in the feedback pathway, and provides a strong prediction on the parallel fiber potentiation time scale. Finally, our model predicts the existence of an optimal learning contrast around 15% contrast levels, which are commonly experienced by interacting fish.     

Characterizing multistationarity regimes in biochemical reaction networks

Switch like responses appear as common strategies in the regulation of cellular systems. Here we present a method to characterize bistable regimes in biochemical reaction networks that can be of use to both direct and reverse engineering of biological switches. In the design of a synthetic biological switch, it is important to study the capability for bistability of the underlying biochemical network structure. Chemical Reaction Network Theory (CRNT) may help at this level to decide whether a given network has the capacity for multiple positive equilibria, based on their structural properties. However, in order to build a working switch, we also need to ensure that the bistability property is robust, by studying the conditions leading to the existence of two different steady states. In the reverse engineering of biological switches, knowledge collected about the bistable regimes of the underlying potential model structures can contribute at the model identification stage to a drastic reduction of the feasible region in the parameter space of search. In this work, we make use and extend previous results of the CRNT, aiming not only to discriminate whether a biochemical reaction network can exhibit multiple steady states, but also to determine the regions within the whole space of parameters capable of producing multistationarity. To that purpose we present and justify a condition on the parameters of biochemical networks for the appearance of multistationarity, and propose an efficient and reliable computational method to check its satisfaction through the parameter space.