No abstract

During the last 10 years, systems biology has matured from a fuzzy concept combining omics, mathematical modeling and computers into a scientific field on its own right. In spite of its incredible potential, the multilevel complexity of its objects of study makes it very difficult to establish a reliable connection between data and models. The great number of degrees of freedom often results in situations, where many different models can explain/fit all available datasets. This has resulted in a shift of paradigm from the initially dominant, maybe naive, idea of inferring the system out of a number of datasets to the application of different techniques that reduce the degrees of freedom before any data set is analyzed. There is a wide variety of techniques available, each of them can contribute a piece of the puzzle and include different kinds of experimental information. But the challenge that remains is their meaningful integration. Here we show some theoretical results that enable some of the main modeling approaches to be applied sequentially in a complementary manner, and how this workflow can benefit from evolutionary reasoning to keep the complexity of the problem in check. As a proof of concept, we show how the synergies between these modeling techniques can provide insight into some well studied problems: Ammonia assimilation in bacteria and an unbranched linear pathway with end-product inhibition.

Salt tolerance in the γ-proteobacterium Halomonas elongata is linked to its ability to produce the compatible solute ectoine. The metabolism of ectoine production is of great interest since it can shed light on the biochemical basis of halotolerance as well as pave the way for the improvement of the biotechnological production of such compatible solute. The ectoine production pathway uses oxaloacetate as a precursor, thereby connecting ectoine production to the anaplerotic reactions that refill carbon into the TCA cycle. This places a high demand on these reactions and creates the need to regulate them not only in response to growth but also in response to extracellular salt concentration. In this work we combine modeling and experiments to analyze how these different needs shape the anaplerotic reactions in H. elongata. First, the stoichiometric and thermodynamic factors that condition the flux distributions are analyzed, then the optimal patterns of operation for oxaloacetate production are calculated. Finally, the phenotype of two deletion mutants lacking potentially relevant anaplerotic enzymes: Phosphoenolpyruvate carboxylase (Ppc) and Oxaloacetate decarboxylase (Oad) is experimentally characterized. The results show that the anaplerotic reactions in H. elongata are indeed subject to different evolutionary pressures than those of other gram-negative bacteria. Ectoine producing halophiles must meet a higher metabolic demand for oxaloacetate and the reliance of many marine bacteria on the Entner-Doudoroff pathway compromises the anaplerotic efficiency of Ppc, which is usually one of the main enzymes fulfilling this role. The anaplerotic flux in H. elongata is contributed not only by Ppc but also by Oad, an enzyme that has not yet been shown to play this role in vivo. Ppc is necessary for H. elongata to grow normally at low salt concentrations but it is not required to achieve near maximal growth rates as long as there is a steep sodium gradient. On the other hand, the lack of Oad presents serious difficulties to grow at high salt concentrations. This points to a shared role of these two enzymes in guaranteeing the supply of OAA for biosynthetic reactions.