Spontaneous reaction silencing in metabolic optimization

Targeting cancer metabolism through synthetic lethality-based combinatorial treatment strategies

PURPOSE OF REVIEW

Targeting cancer metabolism for therapy has received much attention over the last decade with various small molecule inhibitors entering clinical trials. The present review highlights the latest strategies to target glucose and glutamine metabolism for cancer therapy with a particular emphasis on novel combinatorial treatment approaches.

RECENT FINDINGS

Inhibitors of glucose, lactate, and glutamine transport and the ensuing metabolism are in preclinical to clinical trial stages of investigation. Recent advances in our understanding of cell-intrinsic and cell-extrinsic factors that dictate dependence on these targets have informed the development of rational, synthetic lethality-based strategies to exploit these metabolic vulnerabilities.

SUMMARY

Cancer cells exhibit a number of metabolic alterations with functional consequences beyond that of sustaining cellular energetics and biosynthesis. Elucidating context-specific metabolic dependencies and their connections to oncogenic signaling and epigenetic programs in tumor cells represents a promising approach to identify new metabolic drug targets for cancer therapy.

Flux tope analysis: studying the coordination of reaction directions in metabolic networks

Motivation

Elementary flux mode (EFM) analysis allows an unbiased description of metabolic networks in terms of minimal pathways (involving a minimal set of reactions). To date, the enumeration of EFMs is impracticable in genome-scale metabolic models. In a complementary approach, we introduce the concept of a flux tope (FT), involving a maximal set of reactions (with fixed directions), which allows one to study the coordination of reaction directions in metabolic networks and opens a new way for EFM enumeration.

Results

A FT is a (nontrivial) subset of the flux cone specified by fixing the directions of all reversible reactions. In a consistent metabolic network (without unused reactions), every FT contains a 'maximal pathway', carrying flux in all reactions. This decomposition of the flux cone into FTs allows the enumeration of EFMs (of individual FTs) without increasing the problem dimension by reaction splitting. To develop a mathematical framework for FT analysis, we build on the concepts of sign vectors and hyperplane arrangements. Thereby, we observe that FT analysis can be applied also to flux optimization problems involving additional (inhomogeneous) linear constraints. For the enumeration of FTs, we adapt the reverse search algorithm and provide an efficient implementation. We demonstrate that (biomass-optimal) FTs can be enumerated in genome-scale metabolic models of B.cuenoti and E.coli, and we use FTs to enumerate EFMs in models of M.genitalium and B.cuenoti.

Availability and implementation

The source code is freely available at https://github.com/mpgerstl/FTA.

Supplementary information

Supplementary data are available at Bioinformatics online.

Experimental evolution of diverse Escherichia coli metabolic mutants identifies genetic loci for convergent adaptation of growth rate

Cell growth is determined by substrate availability and the cell’s metabolic capacity to assimilate substrates into building blocks. Metabolic genes that determine growth rate may interact synergistically or antagonistically, and can accelerate or slow growth, depending on genetic background and environmental conditions. We evolved a diverse set of Escherichia coli single-gene deletion mutants with a spectrum of growth rates and identified mutations that generally increase growth rate. Despite the metabolic differences between parent strains, mutations that enhanced growth largely mapped to core transcription machinery, including the β and β’ subunits of RNA polymerase (RNAP) and the transcription elongation factor, NusA. The structural segments of RNAP that determine enhanced growth have been previously implicated in antibiotic resistance and in the control of transcription elongation and pausing. We further developed a computational framework to characterize how the transcriptional changes that occur upon acquisition of these mutations affect growth rate across strains. Our experimental and computational results provide evidence for cases in which RNAP mutations shift the competitive balance between active transcription and gene silencing. This study demonstrates that mutations in specific regions of RNAP are a convergent adaptive solution that can enhance the growth rate of cells from distinct metabolic states.

The loss of a metabolic function caused by gene deletion can be compensated, in certain cases, by the concurrent mutation of a second gene. Whether such gene pairs share a local chemical or regulatory relationship or interact via a non-local mechanism has implications for the co-evolution of genetic changes, development of alternatives to gene therapy, and the design of combination antimicrobial therapies that select against resistance. Yet, we lack a comprehensive knowledge of adaptive responses to metabolic mutations, and our understanding of the mechanisms underlying genetic rescue remains limited. We present results of a laboratory evolution approach that has the potential to address both challenges, showing that mutations in specific regions of RNA polymerase enhance growth rates of distinct mutant strains of Escherichia coli with a spectrum of growth defects. Several of these adaptive mutations are deleterious when engineered directly into the original wild-type strain under alternative cultivation conditions, and thus have epistatic rescue properties when paired with the corresponding primary metabolic gene deletions. Our combination of adaptive evolution, directed genetic engineering, and mathematical analysis of transcription and growth rate distinguishes between rescue interactions that are specific or non-specific to a particular deletion. Our study further supports a model for RNA polymerase as a locus of convergent adaptive evolution from different sub-optimal metabolic starting points.

Metabolic Model Reconstruction and Analysis of an Artificial Microbial Ecosystem

Microbial communities are widespread in the environment, and to isolate and identify species or to determine relations among microorganisms, some 'omics methods like metagenomics, proteomics, and metabolomics have been used. When combined with various 'omics data, models known as artificial microbial ecosystems (AME) are powerful methods that can make functional predictions about microbial communities. Reconstruction of an AME model is the first step for model analysis. Many techniques have been applied to the construction of AME models, e.g., the compartmentalization approach, community objectives method, and dynamic analysis approach. Of these approaches, species compartmentalization is the most relevant to genetics. Besides, some algorithms have been developed for the analysis of AME models. In this chapter, we present a general protocol for the use of the species compartmentalization method to reconstruct a model of microbial communities. Then, the analysis of an AME is discussed.

A simple strategy guides the complex metabolic regulation in Escherichia coli

A way to decipher the complexity of the cellular metabolism is to study the effect of different external perturbations. Through an analysis over a sufficiently large set of gene knockouts and growing conditions, one aims to find a unifying principle that governs the metabolic regulation. For instance, it is known that the cessation of the microorganism proliferation after a gene deletion is only transient. However, we do not know the guiding principle that determines the partial or complete recovery of the growth rate, the corresponding redistribution of the metabolic fluxes and the possible different phenotypes. In spite of this large variety in the observed metabolic adjustments, we show that responses of E. coli to several different perturbations can always be derived from a sequence of greedy and myopic resilencings. This simple mechanism provides a detailed explanation for the experimental dynamics both at cellular (proliferation rate) and molecular level (¹³C-determined fluxes), also in case of appearance of multiple phenotypes. As additional support, we identified an example of a simple network motif that is capable of implementing this myopic greediness in the regulation of the metabolism.

Cellular Growth Arrest and Persistence from Enzyme Saturation

Metabolic efficiency depends on the balance between supply and demand of metabolites, which is sensitive to environmental and physiological fluctuations, or noise, causing shortages or surpluses in the metabolic pipeline. How cells can reliably optimize biomass production in the presence of metabolic fluctuations is a fundamental question that has not been fully answered. Here we use mathematical models to predict that enzyme saturation creates distinct regimes of cellular growth, including a phase of growth arrest resulting from toxicity of the metabolic process. Noise can drive entry of single cells into growth arrest while a fast-growing majority sustains the population. We confirmed these predictions by measuring the growth dynamics of Escherichia coli utilizing lactose as a sole carbon source. The predicted heterogeneous growth emerged at high lactose concentrations, and was associated with cell death and production of antibiotic-tolerant persister cells. These results suggest how metabolic networks may balance costs and benefits, with important implications for drug tolerance.

In bacteria, changes in gene expression, with resulting changes in protein concentration, can drastically change how fast cells and cellular populations grow. This fact has big implications for how we treat infectious disease, which types of organisms make up our microbiomes, and what patterns of gene regulation have undergone evolutionary selection. Here, we show how, in principle, the expression level of a single enzyme can affect bacterial population growth by creating a threshold where cells grow optimally fast just below it, but rapidly reach a state of no growth just above it because metabolic byproducts build up and halt growth. The narrow margin between these two states makes entering either of them possible for the same bacterium because of intrinsic uncertainty, or "noise", in gene expression. The predicted result is a variety of growth rates in a single population of genetically identical cells, manifested as a mix of fast- and slow-growing cells. We created laboratory conditions that reproduce the effect in the model organism E. coli, and showed that there may be a benefit to having slower growing cells, because they can survive antibiotic exposure for longer.

Metabolic Adaptation Processes That Converge to Optimal Biomass Flux Distributions

In simple organisms like E.coli, the metabolic response to an external perturbation passes through a transient phase in which the activation of a number of latent pathways can guarantee survival at the expenses of growth. Growth is gradually recovered as the organism adapts to the new condition. This adaptation can be modeled as a process of repeated metabolic adjustments obtained through the resilencings of the non-essential metabolic reactions, using growth rate as selection probability for the phenotypes obtained. The resulting metabolic adaptation process tends naturally to steer the metabolic fluxes towards high growth phenotypes. Quite remarkably, when applied to the central carbon metabolism of E.coli, it follows that nearly all flux distributions converge to the flux vector representing optimal growth, i.e., the solution of the biomass optimization problem turns out to be the dominant attractor of the metabolic adaptation process.

In modeling metabolic networks, concepts like biomass optimization are often used to determine flux distributions of simple organisms such as E.coli. Although they often give good results in practice, they normally rely on heuristic considerations like “evolution has tuned metabolic fluxes to optimize growth, hence optimizing growth gives reasonable fluxes”. The main result of this paper is to show that metabolic adaptation naturally leads to optimal growth, in the sense that the flux distribution associated to optimal growth is the dominant attractor of the fitness landscape of the metabolic adaptation process.

Metabolic flux analysis of Escherichia coli knockouts: lessons from the Keio collection and future outlook

Cellular metabolic and regulatory systems are of fundamental interest to biologists and engineers. Incomplete understanding of these complex systems remains an obstacle to progress in biotechnology and metabolic engineering. An established method for obtaining new information on network structure, regulation and dynamics is to study the cellular system following a perturbation such as a genetic knockout. The Keio collection of all viable Escherichia coli single-gene knockouts is facilitating a systematic investigation of the regulation and metabolism of E. coli. Of all omics measurements available, the metabolic flux profile (the fluxome) provides the most direct and relevant representation of the cellular phenotype. Recent advances in (13)C-metabolic flux analysis are now permitting highly precise and accurate flux measurements for investigating cellular systems and guiding metabolic engineering efforts.

Parallel exploitation of diverse host nutrients enhances Salmonella virulence

Pathogen access to host nutrients in infected tissues is fundamental for pathogen growth and virulence, disease progression, and infection control. However, our understanding of this crucial process is still rather limited because of experimental and conceptual challenges. Here, we used proteomics, microbial genetics, competitive infections, and computational approaches to obtain a comprehensive overview of Salmonella nutrition and growth in a mouse typhoid fever model. The data revealed that Salmonella accessed an unexpectedly diverse set of at least 31 different host nutrients in infected tissues but the individual nutrients were available in only scarce amounts. Salmonella adapted to this situation by expressing versatile catabolic pathways to simultaneously exploit multiple host nutrients. A genome-scale computational model of Salmonella in vivo metabolism based on these data was fully consistent with independent large-scale experimental data on Salmonella enzyme quantities, and correctly predicted 92% of 738 reported experimental mutant virulence phenotypes, suggesting that our analysis provided a comprehensive overview of host nutrient supply, Salmonella metabolism, and Salmonella growth during infection. Comparison of metabolic networks of other pathogens suggested that complex host/pathogen nutritional interfaces are a common feature underlying many infectious diseases.

Basic and applied uses of genome-scale metabolic network reconstructions of Escherichia coli

The genome-scale model (GEM) of metabolism in the bacterium Escherichia coli K-12 has been in development for over a decade and is now in wide use. GEM-enabled studies of E. coli have been primarily focused on six applications: (1) metabolic engineering, (2) model-driven discovery, (3) prediction of cellular phenotypes, (4) analysis of biological network properties, (5) studies of evolutionary processes, and (6) models of interspecies interactions. In this review, we provide an overview of these applications along with a critical assessment of their successes and limitations, and a perspective on likely future developments in the field. Taken together, the studies performed over the past decade have established a genome-scale mechanistic understanding of genotype-phenotype relationships in E. coli metabolism that forms the basis for similar efforts for other microbial species. Future challenges include the expansion of GEMs by integrating additional cellular processes beyond metabolism, the identification of key constraints based on emerging data types, and the development of computational methods able to handle such large-scale network models with sufficient accuracy.

Design constraints on a synthetic metabolism

A metabolism is a complex network of chemical reactions that converts sources of energy and chemical elements into biomass and other molecules. To design a metabolism from scratch and to implement it in a synthetic genome is almost within technological reach. Ideally, a synthetic metabolism should be able to synthesize a desired spectrum of molecules at a high rate, from multiple different nutrients, while using few chemical reactions, and producing little or no waste. Not all of these properties are achievable simultaneously. We here use a recently developed technique to create random metabolic networks with pre-specified properties to quantify trade-offs between these and other properties. We find that for every additional molecule to be synthesized a network needs on average three additional reactions. For every additional carbon source to be utilized, it needs on average two additional reactions. Networks able to synthesize 20 biomass molecules from each of 20 alternative sole carbon sources need to have at least 260 reactions. This number increases to 518 reactions for networks that can synthesize more than 60 molecules from each of 80 carbon sources. The maximally achievable rate of biosynthesis decreases by approximately 5 percent for every additional molecule to be synthesized. Biochemically related molecules can be synthesized at higher rates, because their synthesis produces less waste. Overall, the variables we study can explain 87 percent of variation in network size and 84 percent of the variation in synthesis rate. The constraints we identify prescribe broad boundary conditions that can help to guide synthetic metabolism design.

iAK692: a genome-scale metabolic model of Spirulina platensis C1

Background

Spirulina (Arthrospira) platensis is a well-known filamentous cyanobacterium used in the production of many industrial products, including high value compounds, healthy food supplements, animal feeds, pharmaceuticals and cosmetics, for example. It has been increasingly studied around the world for scientific purposes, especially for its genome, biology, physiology, and also for the analysis of its small-scale metabolic network. However, the overall description of the metabolic and biotechnological capabilities of S. platensis requires the development of a whole cellular metabolism model. Recently, the S. platensis C1 (Arthrospira sp. PCC9438) genome sequence has become available, allowing systems-level studies of this commercial cyanobacterium.

Results

In this work, we present the genome-scale metabolic network analysis of S. platensis C1, iAK692, its topological properties, and its metabolic capabilities and functions. The network was reconstructed from the S. platensis C1 annotated genomic sequence using Pathway Tools software to generate a preliminary network. Then, manual curation was performed based on a collective knowledge base and a combination of genomic, biochemical, and physiological information. The genome-scale metabolic model consists of 692 genes, 837 metabolites, and 875 reactions. We validated iAK692 by conducting fermentation experiments and simulating the model under autotrophic, heterotrophic, and mixotrophic growth conditions using COBRA toolbox. The model predictions under these growth conditions were consistent with the experimental results. The iAK692 model was further used to predict the unique active reactions and essential genes for each growth condition. Additionally, the metabolic states of iAK692 during autotrophic and mixotrophic growths were described by phenotypic phase plane (PhPP) analysis.

Conclusions

This study proposes the first genome-scale model of S. platensis C1, iAK692, which is a predictive metabolic platform for a global understanding of physiological behaviors and metabolic engineering. This platform could accelerate the integrative analysis of various “-omics” data, leading to strain improvement towards a diverse range of desired industrial products from Spirulina.

Load-induced modulation of signal transduction networks

Biological signal transduction networks are commonly viewed as circuits that pass along information--in the process amplifying signals, enhancing sensitivity, or performing other signal-processing tasks--to transcriptional and other components. Here, we report on a "reverse-causality" phenomenon, which we call load-induced modulation. Through a combination of analytical and experimental tools, we discovered that signaling was modulated, in a surprising way, by downstream targets that receive the signal and, in doing so, apply what in physics is called a load. Specifically, we found that non-intuitive changes in response dynamics occurred for a covalent modification cycle when load was present. Loading altered the response time of a system, depending on whether the activity of one of the enzymes was maximal and the other was operating at its minimal rate or whether both enzymes were operating at submaximal rates. These two conditions, which we call "limit regime" and "intermediate regime," were associated with increased or decreased response times, respectively. The bandwidth, the range of frequency in which the system can process information, decreased in the presence of load, suggesting that downstream targets participate in establishing a balance between noise-filtering capabilities and a circuit's ability to process high-frequency stimulation. Nodes in a signaling network are not independent relay devices, but rather are modulated by their downstream targets.

The molecular origins of evolutionary innovations

The history of life is a history of evolutionary innovations, qualitatively new phenotypic traits that endow their bearers with new, often game-changing abilities. We know many individual examples of innovations and their natural history, but we know little about the fundamental principles of phenotypic variability that permit new phenotypes to arise. Most phenotypic innovations result from changes in three classes of systems: metabolic networks, regulatory circuits, and macromolecules. I here highlight two important features that these classes of systems share. The first is the ubiquity of vast genotype networks - connected sets of genotypes with the same phenotype. The second is the great phenotypic diversity of small neighborhoods around different genotypes in genotype space. I here explain that both features are essential for the phenotypic variability that can bring forth qualitatively new phenotypes. Both features emerge from a common cause, the robustness of phenotypes to perturbations, whose origins are linked to life in changing environments.

Metabolic network reconstruction and genome-scale model of butanol-producing strain Clostridium beijerinckii NCIMB 8052

BACKGROUND

Solventogenic clostridia offer a sustainable alternative to petroleum-based production of butanol--an important chemical feedstock and potential fuel additive or replacement. C. beijerinckii is an attractive microorganism for strain design to improve butanol production because it (i) naturally produces the highest recorded butanol concentrations as a byproduct of fermentation; and (ii) can co-ferment pentose and hexose sugars (the primary products from lignocellulosic hydrolysis). Interrogating C. beijerinckii metabolism from a systems viewpoint using constraint-based modeling allows for simulation of the global effect of genetic modifications.

RESULTS

We present the first genome-scale metabolic model (iCM925) for C. beijerinckii, containing 925 genes, 938 reactions, and 881 metabolites. To build the model we employed a semi-automated procedure that integrated genome annotation information from KEGG, BioCyc, and The SEED, and utilized computational algorithms with manual curation to improve model completeness. Interestingly, we found only a 34% overlap in reactions collected from the three databases--highlighting the importance of evaluating the predictive accuracy of the resulting genome-scale model. To validate iCM925, we conducted fermentation experiments using the NCIMB 8052 strain, and evaluated the ability of the model to simulate measured substrate uptake and product production rates. Experimentally observed fermentation profiles were found to lie within the solution space of the model; however, under an optimal growth objective, additional constraints were needed to reproduce the observed profiles--suggesting the existence of selective pressures other than optimal growth. Notably, a significantly enriched fraction of actively utilized reactions in simulations--constrained to reflect experimental rates--originated from the set of reactions that overlapped between all three databases (P = 3.52 × 10-9, Fisher's exact test). Inhibition of the hydrogenase reaction was found to have a strong effect on butanol formation--as experimentally observed.

CONCLUSIONS

Microbial production of butanol by C. beijerinckii offers a promising, sustainable, method for generation of this important chemical and potential biofuel. iCM925 is a predictive model that can accurately reproduce physiological behavior and provide insight into the underlying mechanisms of microbial butanol production. As such, the model will be instrumental in efforts to better understand, and metabolically engineer, this microorganism for improved butanol production.

Dispensability of Escherichia coli's latent pathways

Gene-knockout experiments on single-cell organisms have established that expression of a substantial fraction of genes is not needed for optimal growth. This problem acquired a new dimension with the recent discovery that environmental and genetic perturbations of the bacterium Escherichia coli are followed by the temporary activation of a large number of latent metabolic pathways, which suggests the hypothesis that temporarily activated reactions impact growth and hence facilitate adaptation in the presence of perturbations. Here, we test this hypothesis computationally and find, surprisingly, that the availability of latent pathways consistently offers no growth advantage and tends, in fact, to inhibit growth after genetic perturbations. This is shown to be true even for latent pathways with a known function in alternate conditions, thus extending the significance of this adverse effect beyond apparently nonessential genes. These findings raise the possibility that latent pathway activation is in fact derivative of another, potentially suboptimal, adaptive response.

Use of genome-scale metabolic models in evolutionary systems biology

One of the major aims of the nascent field of evolutionary systems biology is to test evolutionary hypotheses that are not only realistic from a population genetic point of view but also detailed in terms of molecular biology mechanisms. By providing a mapping between genotype and phenotype for hundreds of genes, genome-scale systems biology models of metabolic networks have already provided valuable insights into the evolution of metabolic gene contents and phenotypes of yeast and other microbial species. Here we review the recent use of these computational models to predict the fitness effect of mutations, genetic interactions, evolutionary outcomes, and to decipher the mechanisms of mutational robustness. While these studies have demonstrated that even simplified models of biochemical reaction networks can be highly informative for evolutionary analyses, they have also revealed the weakness of this modeling framework to quantitatively predict mutational effects, a challenge that needs to be addressed for future progress in evolutionary systems biology.

Improving the iMM904 S. cerevisiae metabolic model using essentiality and synthetic lethality data

BACKGROUND

Saccharomyces cerevisiae is the first eukaryotic organism for which a multi-compartment genome-scale metabolic model was constructed. Since then a sequence of improved metabolic reconstructions for yeast has been introduced. These metabolic models have been extensively used to elucidate the organizational principles of yeast metabolism and drive yeast strain engineering strategies for targeted overproductions. They have also served as a starting point and a benchmark for the reconstruction of genome-scale metabolic models for other eukaryotic organisms. In spite of the successive improvements in the details of the described metabolic processes, even the recent yeast model (i.e., iMM904) remains significantly less predictive than the latest E. coli model (i.e., iAF1260). This is manifested by its significantly lower specificity in predicting the outcome of grow/no grow experiments in comparison to the E. coli model.

RESULTS

In this paper we make use of the automated GrowMatch procedure for restoring consistency with single gene deletion experiments in yeast and extend the procedure to make use of synthetic lethality data using the genome-scale model iMM904 as a basis. We identified and vetted using literature sources 120 distinct model modifications including various regulatory constraints for minimal and YP media. The incorporation of the suggested modifications led to a substantial increase in the fraction of correctly predicted lethal knockouts (i.e., specificity) from 38.84% (87 out of 224) to 53.57% (120 out of 224) for the minimal medium and from 24.73% (45 out of 182) to 40.11% (73 out of 182) for the YP medium. Synthetic lethality predictions improved from 12.03% (16 out of 133) to 23.31% (31 out of 133) for the minimal medium and from 6.96% (8 out of 115) to 13.04% (15 out of 115) for the YP medium.

CONCLUSIONS

Overall, this study provides a roadmap for the computationally driven correction of multi-compartment genome-scale metabolic models and demonstrates the value of synthetic lethals as curation agents.

Thermal adaptation of viruses and bacteria

A previously established multiscale population genetics model posits that fitness can be inferred from the physical properties of proteins under the physiological assumption that a loss of stability by any protein confers the lethal phenotype to an organism. Here, we develop this model further by positing that replication rate (fitness) of a bacterial or viral strain directly depends on the copy number of folded proteins, which determine its replication rate. Using this model, and both numerical and analytical approaches, we studied the adaptation process of bacteria and viruses at varied environmental temperatures. We found that a broad distribution of protein stabilities observed in the model and in experiment is the key determinant of thermal response for viruses and bacteria. Our results explain most of the earlier experimental observations: the striking asymmetry of thermal response curves; the absence of evolutionary tradeoff, which was expected but not found in experiments; correlation between denaturation temperature for several protein families and the optimal growth temperature of their carrier organisms; and proximity of bacterial or viral optimal growth temperatures to their evolutionary temperatures. Our theory quantitatively and with high accuracy described thermal response curves for 35 bacterial species using, for each species, only two adjustable parameters-the number of rate-determining genes and the energy barrier for metabolic reactions.

The renaissance of continuous culture in the post-genomics age

The development of continuous culture techniques 60 years ago and the subsequent formulation of theory and the diversification of experimental systems revolutionised microbiology and heralded a unique period of innovative research. Then, progressively, molecular biology and thence genomics and related high-information-density omics technologies took centre stage and microbial growth physiology in general faded from educational programmes and research funding priorities alike. However, there has been a gathering appreciation over the past decade that if the claims of systems biology are going to be realised, they will have to be based on rigorously controlled and reproducible microbial and cell growth platforms. This revival of continuous culture will be long lasting because its recognition as the growth system of choice is firmly established. The purpose of this review, therefore, is to remind microbiologists, particularly those new to continuous culture approaches, of the legacy of what I call the first age of continuous culture, and to explore a selection of researches that are using these techniques in this post-genomics age. The review looks at the impact of continuous culture across a comprehensive range of microbiological research and development. The ability to establish (quasi-) steady state conditions is a frequently stated advantage of continuous cultures thereby allowing environmental parameters to be manipulated without causing concomitant changes in the specific growth rate. However, the use of continuous cultures also enables the critical study of specified transition states and chemical, physical or biological perturbations. Such dynamic analyses enhance our understanding of microbial ecology and microbial pathology for example, and offer a wider scope for innovative drug discovery; they also can inform the optimization of batch and fed-batch operations that are characterized by sequential transitions states.

Network synchronization landscape reveals compensatory structures, quantization, and the positive effect of negative interactions

Synchronization, in which individual dynamical units keep in pace with each other in a decentralized fashion, depends both on the dynamical units and on the properties of the interaction network. Yet, the role played by the network has resisted comprehensive characterization within the prevailing paradigm that interactions facilitating pairwise synchronization also facilitate collective synchronization. Here we challenge this paradigm and show that networks with best complete synchronization, least coupling cost, and maximum dynamical robustness, have arbitrary complexity but quantized total interaction strength, which constrains the allowed number of connections. It stems from this characterization that negative interactions as well as link removals can be used to systematically improve and optimize synchronization properties in both directed and undirected networks. These results extend the recently discovered compensatory perturbations in metabolic networks to the realm of oscillator networks and demonstrate why "less can be more" in network synchronization.

Improved network performance via antagonism: From synthetic rescues to multi-drug combinations

Recent research shows that a faulty or sub-optimally operating metabolic network can often be rescued by the targeted removal of enzyme-coding genes – the exact opposite of what traditional gene therapy would suggest. Predictions go as far as to assert that certain gene knockouts can restore the growth of otherwise nonviable gene-deficient cells. Many questions follow from this discovery: What are the underlying mechanisms? How generalizable is this effect? What are the potential applications? Here, I approach these questions from the perspective of compensatory perturbations on networks. Relations are drawn between such synthetic rescues and naturally occurring cascades of reaction inactivation, as well as their analogs in physical and other biological networks. I specially discuss how rescue interactions can lead to the rational design of antagonistic drug combinations that select against resistance and how they can illuminate medical research on cancer, antibiotics, and metabolic diseases.

Applications of genome-scale metabolic reconstructions

The availability and utility of genome-scale metabolic reconstructions have exploded since the first genome-scale reconstruction was published a decade ago. Reconstructions have now been built for a wide variety of organisms, and have been used toward five major ends: (1) contextualization of high-throughput data, (2) guidance of metabolic engineering, (3) directing hypothesis-driven discovery, (4) interrogation of multi-species relationships, and (5) network property discovery. In this review, we examine the many uses and future directions of genome-scale metabolic reconstructions, and we highlight trends and opportunities in the field that will make the greatest impact on many fields of biology.

Lethal mutagenesis in viruses and bacteria

In this work we study how mutations that change physical properties of cell proteins (stability) affect population survival and growth. We present a model in which the genotype is presented as a set folding free energies of cell proteins. Mutations occur upon replication, so stabilities of some proteins in daughter cells differ from those in the parent cell by amounts deduced from the distribution of mutational effects on protein stability. The genotype-phenotype relationship posits that the cell's fitness (replication rate) is proportional to the concentration of its folded proteins and that unstable essential proteins result in lethality. Simulations reveal that lethal mutagenesis occurs at a mutation rate close to seven mutations in each replication of the genome for RNA viruses and at about half that rate for DNA-based organisms, in accord with earlier predictions from analytical theory and experimental results. This number appears somewhat dependent on the number of genes in the organisms and the organism's natural death rate. Further, our model reproduces the distribution of stabilities of natural proteins, in excellent agreement with experiments. We find that species with high mutation rates tend to have less stable proteins compared to species with low mutation rates.

Conservation of high-flux backbone in alternate optimal and near-optimal flux distributions of metabolic networks

Constraint-based flux balance analysis (FBA) has proven successful in predicting the flux distribution of metabolic networks in diverse environmental conditions. FBA finds one of the alternate optimal solutions that maximizes the biomass production rate. Almaas et al. have shown that the flux distribution follows a power law, and it is possible to associate with most metabolites two reactions which maximally produce and consume a given metabolite, respectively. This observation led to the concept of high-flux backbone (HFB) in metabolic networks. In previous work, the HFB has been computed using a particular optima obtained using FBA. In this paper, we investigate the conservation of HFB of a particular solution for a given medium across different alternate optima and near-optima in metabolic networks of E. coli and S. cerevisiae. Using flux variability analysis (FVA), we propose a method to determine reactions that are guaranteed to be in HFB regardless of alternate solutions. We find that the HFB of a particular optima is largely conserved across alternate optima in E. coli, while it is only moderately conserved in S. cerevisiae. However, the HFB of a particular near-optima shows a large variation across alternate near-optima in both organisms. We show that the conserved set of reactions in HFB across alternate near-optima has a large overlap with essential reactions and reactions which are both uniquely consuming (UC) and uniquely producing (UP). Our findings suggest that the structure of the metabolic network admits a high degree of redundancy and plasticity in near-optimal flow patterns enhancing system robustness for a given environmental condition.