Computational Design of Auxotrophy-Dependent Microbial Biosensors for Combinatorial Metabolic Engineering Experiments

Combinatorial approaches in metabolic engineering work by generating genetic diversity in a microbial population followed by screening for strains with improved phenotypes. One of the most common goals in this field is the generation of a high rate chemical producing strain. A major hurdle with this approach is that many chemicals do not have easy to recognize attributes, making their screening expensive and time consuming. To address this problem, it was previously suggested to use microbial biosensors to facilitate the detection and quantification of chemicals of interest. Here, we present novel computational methods to: (i) rationally design microbial biosensors for chemicals of interest based on substrate auxotrophy that would enable their high-throughput screening; (ii) predict engineering strategies for coupling the synthesis of a chemical of interest with the production of a proxy metabolite for which high-throughput screening is possible via a designed bio-sensor. The biosensor design method is validated based on known genetic modifications in an array of E. coli strains auxotrophic to various amino-acids. Predicted chemical production rates achievable via the biosensor-based approach are shown to potentially improve upon those predicted by current rational strain design approaches. (A Matlab implementation of the biosensor design method is available via http://www.cs.technion.ac.il/~tomersh/tools)

Flux-sum analysis: a metabolite-centric approach for understanding the metabolic network

<p>Abstract</p> <p>Background</p> <p>Constraint-based flux analysis of metabolic network model quantifies the reaction flux distribution to characterize the state of cellular metabolism. However, metabolites are key players in the metabolic network and the current reaction-centric approach may not account for the effect of metabolite perturbation on the cellular physiology due to the inherent limitation in model formulation. Thus, it would be practical to incorporate the metabolite states into the model for the analysis of the network.</p> <p>Results</p> <p>Presented herein is a metabolite-centric approach of analyzing the metabolic network by including the turnover rate of metabolite, known as flux-sum, as key descriptive variable within the model formulation. By doing so, the effect of varying metabolite flux-sum on physiological change can be simulated by resorting to mixed integer linear programming. From the results, we could classify various metabolite types based on the flux-sum profile. Using the <it>i</it>AF1260 <it>in silico </it>metabolic model of <it>Escherichia coli</it>, we demonstrated that this novel concept complements the conventional reaction-centric analysis.</p> <p>Conclusions</p> <p>Metabolite flux-sum analysis elucidates the roles of metabolites in the network. In addition, this metabolite perturbation analysis identifies the key metabolites, implicating practical application which is achievable through metabolite flux-sum manipulation in the areas of biotechnology and biomedical research.</p

Consistency analysis of metabolic correlation networks

<p>Abstract</p> <p>Background</p> <p>Metabolic correlation networks are derived from the covariance of metabolites in replicates of metabolomics experiments. They constitute an interesting intermediate between topology (i.e. the system's architecture defined by the set of reactions between metabolites) and dynamics (i.e. the metabolic concentrations observed as fluctuations around steady-state values in the metabolic network).</p> <p>Results</p> <p>Here we analyze, how such a correlation network changes over time, and compare the relative positions of metabolites in the correlation networks with those in established metabolic networks derived from genome databases. We find that network similarity indeed decreases with an increasing time difference between these networks during a day/night course and, counter intuitively, that proximity of metabolites in the correlation network is no indicator of proximity of the metabolites in the metabolic network.</p> <p>Conclusion</p> <p>The organizing principles of correlation networks are distinct from those of metabolic reaction maps. Time courses of correlation networks may in the future prove an important data source for understanding these organizing principles.</p

Search for new phenomena in final states with an energetic jet and large missing transverse momentum in pp collisions at √ s = 8 TeV with the ATLAS detector

Results of a search for new phenomena in final states with an energetic jet and large missing transverse momentum are reported. The search uses 20.3 fb−1 of √ s = 8 TeV data collected in 2012 with the ATLAS detector at the LHC. Events are required to have at least one jet with pT > 120 GeV and no leptons. Nine signal regions are considered with increasing missing transverse momentum requirements between Emiss T > 150 GeV and Emiss T > 700 GeV. Good agreement is observed between the number of events in data and Standard Model expectations. The results are translated into exclusion limits on models with either large extra spatial dimensions, pair production of weakly interacting dark matter candidates, or production of very light gravitinos in a gauge-mediated supersymmetric model. In addition, limits on the production of an invisibly decaying Higgs-like boson leading to similar topologies in the final state are presente

The Carbon Assimilation Network in Escherichia coli Is Densely Connected and Largely Sign-Determined by Directions of Metabolic Fluxes

Gene regulatory networks consist of direct interactions but also include indirect interactions mediated by metabolites and signaling molecules. We describe how these indirect interactions can be derived from a model of the underlying biochemical reaction network, using weak time-scale assumptions in combination with sensitivity criteria from metabolic control analysis. We apply this approach to a model of the carbon assimilation network in Escherichia coli. Our results show that the derived gene regulatory network is densely connected, contrary to what is usually assumed. Moreover, the network is largely sign-determined, meaning that the signs of the indirect interactions are fixed by the flux directions of biochemical reactions, independently of specific parameter values and rate laws. An inversion of the fluxes following a change in growth conditions may affect the signs of the indirect interactions though. This leads to a feedback structure that is at the same time robust to changes in the kinetic properties of enzymes and that has the flexibility to accommodate radical changes in the environment

The logic of kinetic regulation in the thioredoxin system

<p>Abstract</p> <p>Background</p> <p>The thioredoxin system consisting of NADP(H), thioredoxin reductase and thioredoxin provides reducing equivalents to a large and diverse array of cellular processes. Despite a great deal of information on the kinetics of individual thioredoxin-dependent reactions, the kinetic regulation of this system as an integrated whole is not known. We address this by using kinetic modeling to identify and describe kinetic behavioral motifs found within the system.</p> <p>Results</p> <p>Analysis of a realistic computational model of the <it>Escherichia coli </it>thioredoxin system revealed several modes of kinetic regulation in the system. In keeping with published findings, the model showed that thioredoxin-dependent reactions were adaptable (i.e. changes to the thioredoxin system affected the kinetic profiles of these reactions). Further and in contrast to other systems-level descriptions, analysis of the model showed that apparently unrelated thioredoxin oxidation reactions can affect each other via their combined effects on the thioredoxin redox cycle. However, the scale of these effects depended on the kinetics of the individual thioredoxin oxidation reactions with some reactions more sensitive to changes in the thioredoxin cycle and others, such as the Tpx-dependent reduction of hydrogen peroxide, less sensitive to these changes. The coupling of the thioredoxin and Tpx redox cycles also allowed for ultrasensitive changes in the thioredoxin concentration in response to changes in the thioredoxin reductase concentration. We were able to describe the kinetic mechanisms underlying these behaviors precisely with analytical solutions and core models.</p> <p>Conclusions</p> <p>Using kinetic modeling we have revealed the logic that underlies the functional organization and kinetic behavior of the thioredoxin system. The thioredoxin redox cycle and associated reactions allows for a system that is adaptable, interconnected and able to display differential sensitivities to changes in this redox cycle. This work provides a theoretical, systems-biological basis for an experimental analysis of the thioredoxin system and its associated reactions.</p

In vivo pharmacological evaluations of novel olanzapine analogues in rats: a potential new avenue for the treatment of schizophrenia

Olanzapine (Olz) is one of the most effective antipsychotic drugs commonly used for treating schizophrenia. Unfortunately, Olz administration is associated with severe weight gain and metabolic disturbances. Both patients and clinicians are highly interested in the development of new antipsychotics which are as effective as atypical antipsychotics but which have a lower propensity to induce metabolic side effects. In the present study, we examined two new derivatives of Olz; OlzEt (2-ethyl-4-(4′-methylpiperazin-1′-yl)-10Hbenzo[b]thieno[2,3-e][1,4]diazepine), and OlzHomo (2-ethyl-4-(4′-methyl-1′,4′-diazepan-1′-yl)-10H-benzo[b]thieno[2,3-e] [1,4]diazepine), for their tendency to induce weight gain in rats. Weight gain and metabolic changes were measured in female Sprague Dawley rats. Animals were treated orally with Olz, OlzEt, OlzHomo (3 or 6 mg/kg/day), or vehicle (n = 8), three times daily at eight-hour intervals for 5 weeks. Furthermore, a phencyclidine (PCP)-treated rat model was used to examine the prevention of PCP-induced hyperlocomotor activity relevant for schizophrenia therapy. Male Sprague Dawley rats were pre-treated with a single dose (3 mg/kg/day) of Olz, OlzEt, OlzHomo, or vehicle (n = 12), for 2 weeks. Locomotor activity was recorded following a subcutaneous injection with either saline or PCP (10 mg/kg). Olz was found to induce weight gain, hyperphagia, visceral fat accumulation, and metabolic changes associated with reduced histamatergic H1 receptor density in the hypothalamus of treated rats. In contrast, OlzEt and OlzHomo presented promising antipsychotic effects, which did not induce weight gain or fat deposition in the treated animals. Behavioural analysis showed OlzEt to attenuate PCP-induced hyperactivity to a level similar to that of Olz; however, OlzHomo showed a lower propensity to inhibit these stereotyped behaviours. Our data suggest that the therapeutic effectiveness of OlzHomo may be delivered at a higher dose than that of Olz and OlzEt. Overall, OlzEt and OlzHomo may offer a better pharmacological profile than Olz for treating patients with schizophrenia. Clinical trials are needed to test this hypothesis

Mathematical modeling of intracellular signaling pathways

Dynamic modeling and simulation of signal transduction pathways is an important topic in systems biology and is obtaining growing attention from researchers with experimental or theoretical background. Here we review attempts to analyze and model specific signaling systems. We review the structure of recurrent building blocks of signaling pathways and their integration into more comprehensive models, which enables the understanding of complex cellular processes. The variety of mechanisms found and modeling techniques used are illustrated with models of different signaling pathways. Focusing on the close interplay between experimental investigation of pathways and the mathematical representations of cellular dynamics, we discuss challenges and perspectives that emerge in studies of signaling systems

Mathematical modelling of clostridial acetone-butanol-ethanol fermentation

Clostridial acetone-butanol-ethanol (ABE) fermentation features a remarkable shift in the cellular metabolic activity from acid formation, acidogenesis, to the production of industrial-relevant solvents, solventogensis. In recent decades, mathematical models have been employed to elucidate the complex interlinked regulation and conditions that determine these two distinct metabolic states and govern the transition between them. In this review, we discuss these models with a focus on the mechanisms controlling intra- and extracellular changes between acidogenesis and solventogenesis. In particular, we critically evaluate underlying model assumptions and predictions in the light of current experimental knowledge. Towards this end, we briefly introduce key ideas and assumptions applied in the discussed modelling approaches, but waive a comprehensive mathematical presentation. We distinguish between structural and dynamical models, which will be discussed in their chronological order to illustrate how new biological information facilitates the ‘evolution’ of mathematical models. Mathematical models and their analysis have significantly contributed to our knowledge of ABE fermentation and the underlying regulatory network which spans all levels of biological organization. However, the ties between the different levels of cellular regulation are not well understood. Furthermore, contradictory experimental and theoretical results challenge our current notion of ABE metabolic network structure. Thus, clostridial ABE fermentation still poses theoretical as well as experimental challenges which are best approached in close collaboration between modellers and experimentalists