Macromolecular intelligence in microorganisms

On paradoxes between optimal growth, metabolic control analysis, and flux balance analysis

In Microbiology it is often assumed that growth rate is maximal. This may be taken to suggest that the dependence of the growth rate on every enzyme activity is at the top of an inverse-parabolic function, i.e. that all flux control coefficients should equal zero. This might seem to imply that the sum of these flux control coefficients equals zero. According to the summation law of Metabolic Control Analysis (MCA) the sum of flux control coefficients should equal 1 however. And in Flux Balance Analysis (FBA) catabolism is often limited by a hard bound, causing catabolism to fully control the fluxes, again in apparent contrast with a flux control coefficient of zero. Here we resolve these paradoxes (apparent contradictions) in an analysis that uses the 'Edinburgh pathway', the 'Amsterdam pathway', as well as a generic metabolic network providing the building blocks or Gibbs energy for microbial growth. We review and show that (i) optimization depends on so-called enzyme control coefficients rather than the 'catalytic control coefficients' of MCA's summation law, (ii) when optimization occurs at fixed total protein, the former differ from the latter to the extent that they may all become equal to zero in the optimum state, (iii) in more realistic scenarios of optimization where catalytically inert biomass is compensating or maintenance metabolism is taken into consideration, the optimum enzyme concentrations should not be expected to equal those that maximize the specific growth rate, (iv) optimization may be in terms of yield rather than specific growth rate, which resolves the paradox because the sum of catalytic control coefficients on yield equals 0, (v) FBA effectively maximizes growth yield, and for yield the summation law states 0 rather than 1, thereby removing the paradox, (vi) furthermore, FBA then comes more often to a 'hard optimum' defined by a maximum catabolic flux and a catabolic-enzyme control coefficient of 1. The trade-off between maintenance metabolism and growth is highlighted as worthy of further analysis.

Macromolecular networks and intelligence in microorganisms

Living organisms persist by virtue of complex interactions among many components organized into dynamic, environment-responsive networks that span multiple scales and dimensions. Biological networks constitute a type of information and communication technology (ICT): they receive information from the outside and inside of cells, integrate and interpret this information, and then activate a response. Biological networks enable molecules within cells, and even cells themselves, to communicate with each other and their environment. We have become accustomed to associating brain activity - particularly activity of the human brain - with a phenomenon we call "intelligence." Yet, four billion years of evolution could have selected networks with topologies and dynamics that confer traits analogous to this intelligence, even though they were outside the intercellular networks of the brain. Here, we explore how macromolecular networks in microbes confer intelligent characteristics, such as memory, anticipation, adaptation and reflection and we review current understanding of how network organization reflects the type of intelligence required for the environments in which they were selected. We propose that, if we were to leave terms such as "human" and "brain" out of the defining features of "intelligence," all forms of life - from microbes to humans - exhibit some or all characteristics consistent with "intelligence." We then review advances in genome-wide data production and analysis, especially in microbes, that provide a lens into microbial intelligence and propose how the insights derived from quantitatively characterizing biomolecular networks may enable synthetic biologists to create intelligent molecular networks for biotechnology, possibly generating new forms of intelligence, first in silico and then in vivo.

Nitrogen assimilation in Escherichia coli: putting molecular data into a systems perspective

We present a comprehensive overview of the hierarchical network of intracellular processes revolving around central nitrogen metabolism in Escherichia coli. The hierarchy intertwines transport, metabolism, signaling leading to posttranslational modification, and transcription. The protein components of the network include an ammonium transporter (AmtB), a glutamine transporter (GlnHPQ), two ammonium assimilation pathways (glutamine synthetase [GS]-glutamate synthase [glutamine 2-oxoglutarate amidotransferase {GOGAT}] and glutamate dehydrogenase [GDH]), the two bifunctional enzymes adenylyl transferase/adenylyl-removing enzyme (ATase) and uridylyl transferase/uridylyl-removing enzyme (UTase), the two trimeric signal transduction proteins (GlnB and GlnK), the two-component regulatory system composed of the histidine protein kinase nitrogen regulator II (NRII) and the response nitrogen regulator I (NRI), three global transcriptional regulators called nitrogen assimilation control (Nac) protein, leucine-responsive regulatory protein (Lrp), and cyclic AMP (cAMP) receptor protein (Crp), the glutaminases, and the nitrogen-phosphotransferase system. First, the structural and molecular knowledge on these proteins is reviewed. Thereafter, the activities of the components as they engage together in transport, metabolism, signal transduction, and transcription and their regulation are discussed. Next, old and new molecular data and physiological data are put into a common perspective on integral cellular functioning, especially with the aim of resolving counterintuitive or paradoxical processes featured in nitrogen assimilation. Finally, we articulate what still remains to be discovered and what general lessons can be learned from the vast amounts of data that are available now.

Why in vivo may not equal in vitro - new effectors revealed by measurement of enzymatic activities under the same in vivo-like assay conditions

Does the understanding of the dynamics of biochemical networks in vivo, in terms of the properties of their components determined in vitro, require the latter to be determined all under the same conditions? An in vivo-like assay medium for enzyme activity determination was designed based on the concentrations of the major ionic constituents of the Escherichia coli cytosol: K(+), Na(+), Mg(2+), phosphate, glutamate, sulfate and Cl(-). The maximum capacities (V(max)) of the extracted enzymes of two pathways were determined using both this in vivo-like assay medium and the assay medium specific for each enzyme. The enzyme activities differed between the two assay conditions. Most of the differences could be attributed to unsuspected, pleiotropic effects of K(+) and phosphate. K(+) activated some enzymes (aldolase, enolase and glutamate dehydrogenase) and inhibited others (phosphoglucose isomerase, phosphofructokinase, triosephosphate isomerase, glyceraldehyde 3-phosphate dehydrogenase, phosphoglycerate kinase, phosphoglycerate mutase), whereas phosphate inhibited all glycolytic enzymes and glutamine synthetase but only activated glutamine 2-oxoglutarate amidotransferase. Neither a high glutamate concentration, nor macromolecular crowding affected the glycolytic or nitrogen assimilation enzymes, other than through the product inhibition of glutamate dehydrogenase by glutamate. This strategy of assessing all pathway enzymes kinetically under the same conditions may be necessary to avoid inadvertent differences between in vivo and in vitro biochemistry. It may also serve to reveal otherwise unnoticed pleiotropic regulation, such as that demonstrated in the present study by K(+) and phosphate.

Linking the proteins--elucidation of proteome-scale networks using mass spectrometry

Proteomes are intricate. Typically, thousands of proteins interact through physical association and post-translational modifications (PTMs) to give rise to the emergent functions of cells. Understanding these functions requires one to study proteomes as "systems" rather than collections of individual protein molecules. The abstraction of the interacting proteome to "protein networks" has recently gained much attention, as networks are effective representations, that lose specific molecular details, but provide the ability to see the proteome as a whole. Mostly two aspects of the proteome have been represented by network models: proteome-wide physical protein-protein-binding interactions organized into Protein Interaction Networks (PINs), and proteome-wide PTM relations organized into Protein Signaling Networks (PSNs). Mass spectrometry (MS) techniques have been shown to be essential to reveal both of these aspects on a proteome-wide scale. Techniques such as affinity purification followed by MS have been used to elucidate protein-protein interactions, and MS-based quantitative phosphoproteomics is critical to understand the structure and dynamics of signaling through the proteome. We here review the current state-of-the-art MS-based analytical pipelines for the purpose to characterize proteome-scale networks.

Modeling the role of covalent enzyme modification in Escherichia coli nitrogen metabolism

In the bacterium Escherichia coli, the enzyme glutamine synthetase (GS) converts ammonium into the amino acid glutamine. GS is principally active when the cell is experiencing nitrogen limitation, and its activity is regulated by a bicyclic covalent modification cascade. The advantages of this bicyclic-cascade architecture are poorly understood. We analyze a simple model of the GS cascade in comparison to other regulatory schemes and conclude that the bicyclic cascade is suboptimal for maintaining metabolic homeostasis of the free glutamine pool. Instead, we argue that the lag inherent in the covalent modification of GS slows the response to an ammonium shock and thereby allows GS to transiently detoxify the cell, while maintaining homeostasis over longer times.

A network model for the correlation between epistasis and genomic complexity

The study of genetic interactions (epistasis) is central to the understanding of genome organization and evolution. A general correlation between epistasis and genomic complexity has been recently shown, such that in simpler genomes epistasis is antagonistic on average (mutational effects tend to cancel each other out), whereas a transition towards synergistic epistasis occurs in more complex genomes (mutational effects strengthen each other). Here, we use a simple network model to identify basic features explaining this correlation. We show that, in small networks with multifunctional nodes, lack of redundancy, and absence of alternative pathways, epistasis is antagonistic on average. In contrast, lack of multi-functionality, high connectivity, and redundancy favor synergistic epistasis. Moreover, we confirm the previous finding that epistasis is a covariate of mutational robustness: in less robust networks it tends to be antagonistic whereas in more robust networks it tends to be synergistic. We argue that network features associated with antagonistic epistasis are typically found in simple genomes, such as those of viruses and bacteria, whereas the features associated with synergistic epistasis are more extensively exploited by higher eukaryotes.

Taking the example of computer systems engineering for the analysis of biological cell systems

In this paper, we discuss the potential for the use of engineering methods that were originally developed for the design of embedded computer systems, to analyse biological cell systems. For embedded systems as well as for biological cell systems, design is a feature that defines their identity. The assembly of different components in designs of both systems can vary widely. In contrast to the biology domain, the computer engineering domain has the opportunity to quickly evaluate design options and consequences of its systems by methods for computer aided design and in particular design space exploration. We argue that there are enough concrete similarities between the two systems to assume that the engineering methodology from the computer systems domain, and in particular that related to embedded systems, can be applied to the domain of cellular systems. This will help to understand the myriad of different design options cellular systems have. First we compare computer systems with cellular systems. Then, we discuss exactly what features of engineering methods could aid researchers with the analysis of cellular systems, and what benefits could be gained.

Control of MAPK signalling: from complexity to what really matters

Oncogenesis results from changes in kinetics or in abundance of proteins in signal transduction networks. Recently, it was shown that control of signalling cannot reside in a single gene product, and might well be dispersed over many components. Which of the reactions in these complex networks are most important, and how can the existing molecular information be used to understand why particular genes are oncogenes whereas others are not? We implement a new method to help address such questions. We apply control analysis to a detailed kinetic model of the epidermal growth factor-induced mitogen-activated protein kinase network. We determine the control of each reaction with respect to three biologically relevant characteristics of the output of this network: the amplitude, duration and integrated output of the transient phosphorylation of extracellular signal-regulated kinase (ERK). We confirm that control is distributed, but far from randomly: a small proportion of reactions substantially control signalling. In particular, the activity of Raf is in control of all characteristics of the transient profile of ERK phosphorylation, which may clarify why Raf is an oncogene. Most reactions that really matter for one signalling characteristic are also important for the other characteristics. Our analysis also predicts the effects of mutations and changes in gene expression.

A network of net-workers: report of the Euresco conference on ‘Bacterial Neural Networks’ held at San Feliu (Spain) from 8 to 14 May 2004

In May 2004, over 100 bacteriologists from 19 different countries discussed recent progress in identification and understanding of individual signal transfer mechanisms in bacteria and in the mutual interactions between these systems to form a functional living cell. The meeting was held in San Feliu and supported by ESF and EMBO. In part through the extensive sequencing efforts of the past few years, the bulk of the bacterial signal transfer systems have been resolved and their detailed characterization is revealing such characteristics as signal specificity, signalling rate constants, molecular interaction affinities, subcellular localization, etc., which should provide a solid basis to a computational extension of this field of studies. In parallel, the new genomics techniques are providing tools to characterize the way a collection of such systems interact in an individual cell, to give rise to 'life'. Systems theory provides rational and convenient ways to bring order to the wide range of observables thus obtained. Ultimately, the performance of engineered design will have to prove whether or not we know enough about the processes involved.

Environmental conditions and transcriptional regulation inEscherichia coli: a physiological integrative approach

Bacteria develop a number of devices for sensing, responding, and adapting to different environmental conditions. Understanding within a genomic perspective how the transcriptional machinery of bacteria is modulated, as a response for changing conditions, is a major challenge for biologists. Knowledge of which genes are turned on or turned off under specific conditions is essential for our understanding of cell behavior. In this study we describe how the information pertaining to gene expression and associated growth conditions (even with very little knowledge of the associated regulatory mechanisms) is gathered from the literature and incorporated into RegulonDB, a database on transcriptional regulation and operon organization in E. coli. The link between growth conditions, signal transduction, and transcriptional regulation is modeled in the database in a simple format that highlights biological relevant information. As far as we know, there is no other database that explicitly clarifies the effect of environmental conditions on gene transcription. We discuss how this knowledge constitutes a benchmark that will impact future research aimed at integration of regulatory responses in the cell; for instance, analysis of microarrays, predicting culture behavior in biotechnological processes, and comprehension of dynamics of regulatory networks. This integrated knowledge will contribute to the future goal of modeling the behavior of E. coli as an entire cell. The RegulonDB database can be accessed on the web at the URL: http://www.cifn.unam.mx/Computational_Biology/regulondb/.

Selectivity in overlapping MAP kinase cascades

Some protein kinases operate in more than one mitogen-activated protein-kinase (MAPK) cascade. We here address the question whether specificity of the cascades necessitates physical sequestration of these "promiscuous" kinases (e.g. by binding to scaffolds). A model is constructed, in which two MAPK cascades depend on a single MAP-kinase kinase that is not sequestered in two subpopulations. We show that in this model selective signal transduction is possible provided that there is an additional interaction at the MAP-kinase level, there is no simultaneous activation of more than one response by either signal. We discuss a number of additional interactions that can generate the selectivity, as well as some kinetic features by which this mechanism may be recognized experimentally.

Redundancy, antiredundancy, and the robustness of genomes

Genetic mutations that lead to undetectable or minimal changes in phenotypes are said to reveal redundant functions. Redundancy is common among phenotypes of higher organisms that experience low mutation rates and small population sizes. Redundancy is less common among organisms with high mutation rates and large populations, or among the rapidly dividing cells of multicellular organisms. In these cases, one even observes the opposite tendency: a hypersensitivity to mutation, which we refer to as antiredundancy. In this paper we analyze the evolutionary dynamics of redundancy and antiredundancy. Assuming a cost of redundancy, we find that large populations will evolve antiredundant mechanisms for removing mutants and thereby bolster the robustness of wild-type genomes; whereas small populations will evolve redundancy to ensure that all individuals have a high chance of survival. We propose that antiredundancy is as important for developmental robustness as redundancy, and is an essential mechanism for ensuring tissue-level stability in complex multicellular organisms. We suggest that antiredundancy deserves greater attention in relation to cancer, mitochondrial disease, and virus infection.

Product dependence and bifunctionality compromise the ultrasensitivity of signal transduction cascades

Covalent modification cycles are ubiquitous. Theoretical studies have suggested that they serve to increase sensitivity. However, this suggestion has not been corroborated experimentally in vivo. Here, we demonstrate that the assumptions of the theoretical studies, i.e., irreversibility and absence of product inhibition, were not trivial: when the conversion reactions are close to equilibrium or saturated by their product, "zero-order" ultrasensitivity disappears. For high sensitivities to arise, not only substrate saturation (zero-order) but also high equilibrium constants and low product saturation are required. Many covalent modification cycles are catalyzed by one bifunctional 'ambiguous' enzyme rather than by two independent proteins. This makes high substrate concentration and low product concentration for both reactions of the cycle inconsistent; such modification cycles cannot have high responses. Defining signal strength as ratios of modified (e.g., phosphorylated) over unmodified protein, signal-to-signal response sensitivity equals 1: signal strength should remain constant along a cascade of ambiguous modification cycles. We also show that the total concentration of a signalling effector protein cannot affect the signal emanating from a modification cycle catalyzed by an ambiguous enzyme if the ratio of the two forms of the effector protein is not altered. This finding may explain the experimental result that the pivotal signal transduction protein PII plus its paralogue GlnK do not control steady-state N-signal transduction in Escherichia coli. It also rationalizes the absence of strong phenotypes for many signal-transduction proteins. Emphasis on extent of modification of these proteins is perhaps more urgent than transcriptome analysis.