Macromolecular Crowding, Phase Separation, and Homeostasis in the Orchestration of Bacterial Cellular Functions

Macromolecular crowding affects the activity of proteins and functional macromolecular complexes in all cells, including bacteria. Crowding, together with physicochemical parameters such as pH, ionic strength, and the energy status, influences the structure of the cytoplasm and thereby indirectly macromolecular function. Notably, crowding also promotes the formation of biomolecular condensates by phase separation, initially identified in eukaryotic cells but more recently discovered to play key functions in bacteria. Bacterial cells require a variety of mechanisms to maintain physicochemical homeostasis, in particular in environments with fluctuating conditions, and the formation of biomolecular condensates is emerging as one such mechanism. In this work, we connect physicochemical homeostasis and macromolecular crowding with the formation and function of biomolecular condensates in the bacterial cell and compare the supramolecular structures found in bacteria with those of eukaryotic cells. We focus on the effects of crowding and phase separation on the control of bacterial chromosome replication, segregation, and cell division, and we discuss the contribution of biomolecular condensates to bacterial cell fitness and adaptation to environmental stress.

Potential of native and bioprocessed brewers' spent grains as organic soil amendments

Introduction

The use of novel soil amendments and the exploitation of plant growth-promoting microorganisms are considered promising tools for developing a more sustainable agriculture in times when ensuring high-yield productions with limited resources is essential.

Methods

In this study, the potential of brewers' spent grain (BSG), the major by-product of the brewing industry, as organic soil amendment, was investigated. Bioprocessed BSG, obtained by an enzymatic treatment coupled with fermentation, together with native BSG, were used as amendments in a pot-trial. An integrated analytical approach aimed at assessing the modification of the physicochemical properties of a typical Mediterranean alkaline agricultural soil, and the plant growth-promoting effect on escarole (Cichorium endivia var. Cuartana), was carried out.

Results

The use of biomasses led to soil organic content and total nitrogen content up to 72 and 42% higher, compared to the unamended soils. Moreover, the lower pH and the higher organic acids content doubled phosphorus availability. Although the number of leaves per plant in escaroles from pots amended with native and bioprocessed BSG did not show any difference compared to plants cultivated on unamended pots, the average fresh weight per escarole head, was higher in pots amended with bioprocessed BSG.

Discussion

Hence, the results collected so far encourage BSG application for agricultural purpose, while solving the problem of disposing of such abundant side stream.

Sugar Shock: Probing Streptococcus pyogenes Metabolism Through Bioluminescence Imaging

Streptococcus pyogenes (S. pyogenes) can thrive in its host during an infection, and, as a result, it must be able to respond to external stimuli and available carbon sources. The preclinical use of engineered pathogens capable of constitutive light production may provide real-time information on microbial-specific metabolic processes. In this study, we mapped the central metabolism of a luxABCDE-modified S. pyogenes Xen20 (Strep. Xen20) to its de novo synthesis of luciferase substrates as assessed by the rate of light production in response to different environmental triggers. Previous characterization predicted that the lux operon was under the myo-inositol iolE promotor. In this study, we revealed that supplementation with myo-inositol generated increased Strep. Xen20 luminescence. Surprisingly, when supplemented with infection-relevant carbon sources, such as glucose or glycine, light production was diminished. This was presumably due to the scavenging of pyruvate by L-lactate dehydrogenase (LDH). Inhibition of LDH by its inhibitor, oxamate, partially restored luminescent signal in the presence of glucose, presumably by allowing the resulting pyruvate to proceed to acetyl-coenzyme A (CoA). This phenomenon appeared specific to the lactic acid bacterial metabolism as glucose or glycine did not reduce signal in an analogous luxABCDE-modified Gram-positive pathogen, Staph. Xen29. The Strep. Xen20 cells produced light in a concentration-dependent manner, inversely related to the amount of glucose present. Taken together, our measures of microbial response could provide new information regarding the responsiveness of S. pyogenes metabolism to acute changes in its local environments and cellular health.

MICROBIAL isoprene production: an overview

Isoprene, a volatile C5 hydrocarbon, is a precursor of synthetic rubber and an important building block for a variety of natural products, solely being produced by petrochemical routes. To mitigate the ever-increasing contribution of petrochemical industry to global warming through significant carbon (CO₂) evolution, bio-based process for isoprene production using microbial cell factories have been explored. Highly efficient fermentation-based processes have been studied for little over a decade now with extensive research on the rational strain development for creating robust strains for commercial isoprene production. Most of these studies involved sugars as feedstocks and using naturally occurring isoprene pathways viz., mevalonate and methyl erythritol pathway in E. coli. Recent advances, driven by efforts in reducing environmental pollution, have focused on utilization of inorganic CO₂ by cyanobacteria or syngas from waste gases by acetogens for isoprene production. This review endeavors to capture the latest relevant progress made in rational strain development, metabolic engineering and synthetic biology strategies used, challenges in fermentation process development at lab and commercial scale production of isoprene along with a future perspective pertaining to this area of research.

The Non-phosphorylating Glyceraldehyde-3-Phosphate Dehydrogenase GapN Is a Potential New Drug Target in Streptococcus pyogenes

The strict human pathogen Streptococcus pyogenes causes infections of varying severity, ranging from self-limiting suppurative infections to life-threatening diseases like necrotizing fasciitis or streptococcal toxic shock syndrome. Here, we show that the non-phosphorylating glyceraldehyde-3-phosphate dehydrogenase GapN is an essential enzyme for S. pyogenes. GapN converts glyceraldehyde 3-phosphate into 3-phosphoglycerate coupled to the reduction of NADP to NADPH. The knock-down of gapN by antisense peptide nucleic acids (asPNA) significantly reduces viable bacterial counts of S. pyogenes laboratory and macrolide-resistant clinical strains in vitro. As S. pyogenes lacks the oxidative part of the pentose phosphate pathway, GapN appears to be the major NADPH source for the bacterium. Accordingly, other streptococci that carry a complete pentose phosphate pathway are not prone to asPNA-based gapN knock-down. Determination of the crystal structure of the S. pyogenes GapN apo-enzyme revealed an unusual cis-peptide in proximity to the catalytic binding site. Furthermore, using a structural modeling approach, we correctly predicted competitive inhibition of S. pyogenes GapN by erythrose 4-phosphate, indicating that our structural model can be used for in silico screening of specific GapN inhibitors. In conclusion, the data provided here reveal that GapN is a potential target for antimicrobial substances that selectively kill S. pyogenes and other streptococci that lack the oxidative part of the pentose phosphate pathway.

C/N Ratio and Specific Growth Rate Plays Important Role on Enhancing Isoprene Production in Recombinant Escherichia coli

Effect of fermentation parameters such as C/N ratio, specific growth rate, phosphate limitation, and plasmid instability on enhancing isoprene production is the focus of the current study. Isoprene productivity in the recombinant Escherichia coli K12_MVA strain showed a bell-shaped relationship with specific growth rate in bioreactor studies with isoprene volumetric productivity peaking at 0.35/h. This behavior was depicted by a production inhibition kinetic model which envisaged a serious competition between the cellular growth, acetic acid production, and isoprene biosynthesis. The model equation derived showed a reasonable fit with the experimental values. Judicious control of the growth rates and acetate accumulation by optimizing C/N ratio, phosphate concentration, and intermittent feeding strategy resulted in maximizing the carbon flux towards isoprene. Plasmid instability caused by metabolic burden posed by the presence of dual plasmids on the bacteria was simulated using first-order degradation kinetics. The experimental plasmid loss trend was in accordance with the model simulated trend, where higher plasmid loss correlated with higher specific growth rates. Modulating the growth rate, acetate accumulation, and plasmid instability resulted in achieving maximum isoprene volumetric productivity of 1.125 g/l/h with 46.67% of carbon flux towards isoprene and a isoprene titre of 18 g/l in 16 h fermentation run.

Updates on understanding of probiotic lactic acid bacteria responses to environmental stresses and highlights on proteomic analyses

Probiotics are defined as live microorganisms that improve the health of the host when administered in adequate quantities. Nonetheless, probiotics encounter extreme environmental conditions during food processing or along the gastrointestinal tract. This review discusses different environmental stresses that affect probiotics during food preparation, storage, and along the alimentary canal, including high temperature, low temperature, low and alkaline pH, oxidative stress, high hydrostatic pressure, osmotic pressure, and starvation. The understanding of how probiotics deal with environmental stress and thrive provides useful information to guide the selection of the strains with enhanced performance in specific situations, in food processing or during gastrointestinal transit. In most cases, multiple biological functions are affected upon exposure of the cell to environmental stress. Sensing of sublethal environmental stress can allow for adaptation processes to occur, which can include alterations in the expression of specific proteins.

Delving deeper: Relating the behaviour of a metabolic system to the properties of its components using symbolic metabolic control analysis

High-level behaviour of metabolic systems results from the properties of, and interactions between, numerous molecular components. Reaching a complete understanding of metabolic behaviour based on the system’s components is therefore a difficult task. This problem can be tackled by constructing and subsequently analysing kinetic models of metabolic pathways since such models aim to capture all the relevant properties of the system components and their interactions. Symbolic control analysis is a framework for analysing pathway models in order to reach a mechanistic understanding of their behaviour. By providing algebraic expressions for the sensitivities of system properties, such as metabolic flux or steady-state concentrations, in terms of the properties of individual reactions it allows one to trace the high level behaviour back to these low level components. Here we apply this method to a model of pyruvate branch metabolism in Lactococcus lactis in order to explain a previously observed negative flux response towards an increase in substrate concentration. With this method we are able to show, first, that the sensitivity of flux towards changes in reaction rates (represented by flux control coefficients) is determined by the individual metabolic branches of the pathway, and second, how the sensitivities of individual reaction rates towards their substrates (represented by elasticity coefficients) contribute to this flux control. We also quantify the contributions of enzyme binding and mass-action to enzyme elasticity separately, which allows for an even finer-grained understanding of flux control. These analytical tools allow us to analyse the control properties of a metabolic model and to arrive at a mechanistic understanding of the quantitative contributions of each of the enzymes to this control. Our analysis provides an example of the descriptive power of the general principles of symbolic control analysis.

Production of a highly concentrated probiotic culture of Lactococcus lactis CECT 539 containing high amounts of nisin

In this study, probiotic biomass and nisin productions by Lactococcus lactis CECT 539 were followed in two realkalized fed-batch cultures in diluted whey (DW) supplemented with KH₂PO₄ up to a total phosphorus concentration of 0.459 g/L. Increased biomass (5.12 g/L, 2.18 × 10¹⁰ CFU/mL) and nisin (235.23 BU/mL) concentrations were obtained in the culture fed concentrated whey and concentrated mussel processing waste (CMPW) medium supplemented with glucose up to a concentration of 400 g/L (CMPW + G medium) compared with similar fed-batch fermentations in DW medium. In the second fed-batch fermentation, the feeding medium CMPW + G was supplemented with KH₂PO₄ up to a TP concentration of 3.21 g/L. With this approach, increased production of biomass (5.49 g/L, 2.33 × 10¹⁰ CFU/mL) and nisin (257.59 BU/mL) was obtained. Considering the substantial availability of these wastes at very low prices from local dairy and mussel processing plants in Galicia, their use as culture media could offer an attractive alternative for a low-cost production of probiotic biomass and nisin at a high scale.

Deletion of the L-Lactate Dehydrogenase Gene ldh in Streptococcus pyogenes Leads to a Loss of SpeB Activity and a Hypovirulent Phenotype

Streptococcus pyogenes uses lactic acid fermentation for the generation of ATP. Here, we analyzed the impact of a deletion of the L-lactate dehydrogenase gene ldh on the virulence of S. pyogenes M49. While the ldh deletion does not cause a general growth deficiency in laboratory media, the growth in human blood and plasma is significantly hampered. The ldh deletion strain is furthermore less virulent in a Galleria mellonella infection model. We show that the ldh deletion leads to a decrease in the activity of the cysteine protease SpeB, an important secreted virulence factor of S. pyogenes. The reduced SpeB activity is caused by a hampered autocatalytic activation of the SpeB zymogen into the mature SpeB. The missing SpeB activity furthermore leads to increased plasmin activation and a reduced activation of the contact system on the surface of S. pyogenes. All these effects can be reversed when ldh is reintroduced into the mutant via a plasmid. The results demonstrate a previously unappreciated role for LDH in modulation of SpeB maturation.

13C based proteinogenic amino acid (PAA) and metabolic flux ratio analysis of Lactococcus lactis reveals changes in pentose phosphate (PP) pathway in response to agitation and temperature related stresses

Lactococcus lactis subsp. cremoris MG1363 is an important starter culture for dairy fermentation. During industrial fermentations, L. lactis is constantly exposed to stresses that affect the growth and performance of the bacterium. Although the response of L. lactis to several stresses has been described, the adaptation mechanisms at the level of in vivo fluxes have seldom been described. To gain insights into cellular metabolism, ¹³C metabolic flux analysis and gas chromatography mass spectrometry (GC-MS) were used to measure the flux ratios of active pathways in the central metabolism of L. lactis when subjected to three conditions varying in temperature (30°C, 37°C) and agitation (with and without agitation at 150 rpm). Collectively, the concentrations of proteinogenic amino acids (PAAs) and free fatty acids (FAAs) were compared, and Pearson correlation analysis (r) was calculated to measure the pairwise relationship between PAAs. Branched chain and aromatic amino acids, threonine, serine, lysine and histidine were correlated strongly, suggesting changes in flux regulation in glycolysis, the pentose phosphate (PP) pathway, malic enzyme and anaplerotic reaction catalysed by pyruvate carboxylase (pycA). Flux ratio analysis revealed that glucose was mainly converted by glycolysis, highlighting the stability of L. lactis’ central carbon metabolism despite different conditions. Higher flux ratios through oxaloacetate (OAA) from pyruvate (PYR) reaction in all conditions suggested the activation of pyruvate carboxylate (pycA) in L. lactis, in response to acid stress during exponential phase. Subsequently, more significant flux ratio differences were seen through the oxidative and non-oxidative pentose phosphate (PP) pathways, malic enzyme, and serine and C1 metabolism, suggesting NADPH requirements in response to environmental stimuli. These reactions could play an important role in optimization strategies for metabolic engineering in L. lactis. Overall, the integration of systematic analysis of amino acids and flux ratio analysis provides a systems-level understanding of how L. lactis regulates central metabolism under various conditions.

Transcriptional Reprogramming at Genome-Scale of Lactobacillus plantarum WCFS1 in Response to Olive Oil Challenge

Dietary fats may exert selective pressures on Lactobacillus species, however, knowledge on the mechanisms of adaptation to fat stress in these organisms is still fragmentary. This study was undertaken to gain insight into the mechanisms of adaptation of Lactobacillus plantarum WCFS1 to olive oil challenge by whole genome transcriptional profiling using DNA microarrays. A set of 230 genes were differentially expressed by L. plantarum WCFS1 to respond to this vegetable oil. This response involved elements typical of the stringent response, as indicated by the induction of genes involved in stress-related pathways and downregulation of genes related to processes associated with rapid growth. A set of genes involved in the transport and metabolism of compatible solutes were downregulated, indicating that this organism does not require osmoprotective mechanisms in presence of olive oil. The fatty acid biosynthetic pathway was thoroughly downregulated at the transcriptional level, which coincided with a diminished expression of genes controlled by this pathway in other organisms and that are required for the respiratory function, pyruvate dehydrogenase activity, RNA processing and cell size setting. Finally, a set of genes involved in host-cell signaling by L. plantarum were differentially regulated indicating that olive oil can influence the expression of metabolic traits involved in the crosstalk between this bacterium and the host.

Integrating highly quantitative proteomics and genome-scale metabolic modeling to study pH adaptation in the human pathogen Enterococcus faecalis

Genome-scale metabolic models represent the entirety of metabolic reactions of an organism based on the annotation of the respective genome. These models commonly allow all reactions to proceed concurrently, disregarding the fact that at no point all proteins will be present in a cell. The metabolic reaction space can be constrained to a more physiological state using experimentally obtained information on enzyme abundances. However, high-quality, genome-wide protein measurements have been challenging and typically transcript abundances have been used as a surrogate for protein measurements. With recent developments in mass spectrometry-based proteomics, exemplified by SWATH-MS, the acquisition of highly quantitative proteome-wide data at reasonable throughput has come within reach. Here we present methodology to integrate such proteome-wide data into genome-scale models. We applied this methodology to study cellular changes in Enterococcus faecalis during adaptation to low pH. Our results indicate reduced proton production in the central metabolism and decreased membrane permeability for protons due to different membrane composition. We conclude that proteomic data constrain genome-scale models to a physiological state and, in return, genome-scale models are useful tools to contextualize proteomic data.

Genomic characterization of symbiotic mycoplasmas from the stomach of deep-sea isopod bathynomus sp

Deep-sea isopod scavengers such as Bathynomus sp. are able to live in nutrient-poor environments, which is likely attributable to the presence of symbiotic microbes in their stomach. In this study we recovered two draft genomes of mycoplasmas, Bg1 and Bg2, from the metagenomes of the stomach contents and stomach sac of a Bathynomus sp. sample from the South China Sea (depth of 898 m). Phylogenetic trees revealed a considerable genetic distance to other mycoplasma species for Bg1 and Bg2. Compared with terrestrial symbiotic mycoplasmas, the Bg1 and Bg2 genomes were enriched with genes encoding phosphoenolpyruvate-dependent phosphotransferase systems (PTSs) and sodium-driven symporters responsible for the uptake of sugars, amino acids and other carbohydrates. The genome of mycoplasma Bg1 contained sialic acid lyase and transporter genes, potentially enabling the bacteria to attach to the stomach sac and obtain organic carbons from various cell walls. Both of the mycoplasma genomes contained multiple copies of genes related to proteolysis and oligosaccharide degradation, which may help the host survive in low-nutrient conditions. The discovery of the different types of mycoplasma bacteria in the stomach of this deep-sea isopod affords insights into symbiotic model of deep-sea animals and genomic plasticity of mycoplasma bacteria.

Metabolic traits of pathogenic streptococci

Invasive and noninvasive diseases caused by facultative pathogenic streptococci depend on their equipment with virulence factors and on their ability to sense and adapt to changing nutrients in different host environments. The knowledge of the principal metabolic mechanisms which allow these bacteria to recognize and utilize nutrients in host habitats is a prerequisite for our understanding of streptococcal pathogenicity and the development of novel control strategies. This review aims to summarize and compare the central carbohydrate metabolic and amino acid biosynthetic pathways of a selected group of streptococcal species, all belonging to the naso-oropharyngeal microbiome in humans and/or animals. We also discuss the urgent need of comprehensive metabolomics approaches for a better understanding of the streptococcal metabolism during host-pathogen interaction.

Genome-scale reconstruction of the Streptococcus pyogenes M49 metabolic network reveals growth requirements and indicates potential drug targets

Genome-scale metabolic models comprise stoichiometric relations between metabolites, as well as associations between genes and metabolic reactions and facilitate the analysis of metabolism. We computationally reconstructed the metabolic network of the lactic acid bacterium Streptococcus pyogenes M49. Initially, we based the reconstruction on genome annotations and already existing and curated metabolic networks of Bacillus subtilis, Escherichia coli, Lactobacillus plantarum and Lactococcus lactis. This initial draft was manually curated with the final reconstruction accounting for 480 genes associated with 576 reactions and 558 metabolites. In order to constrain the model further, we performed growth experiments of wild type and arcA deletion strains of S. pyogenes M49 in a chemically defined medium and calculated nutrient uptake and production fluxes. We additionally performed amino acid auxotrophy experiments to test the consistency of the model. The established genome-scale model can be used to understand the growth requirements of the human pathogen S. pyogenes and define optimal and suboptimal conditions, but also to describe differences and similarities between S. pyogenes and related lactic acid bacteria such as L. lactis in order to find strategies to reduce the growth of the pathogen and propose drug targets.

Targeting glycolysis in the malaria parasite Plasmodium falciparum

Glycolysis is the main pathway for ATP production in the malaria parasite Plasmodium falciparum and essential for its survival. Following a sensitivity analysis of a detailed kinetic model for glycolysis in the parasite, the glucose transport reaction was identified as the step whose activity needed to be inhibited to the least extent to result in a 50% reduction in glycolytic flux. In a subsequent inhibitor titration with cytochalasin B, we confirmed the model analysis experimentally and measured a flux control coefficient of 0.3 for the glucose transporter. In addition to the glucose transporter, the glucokinase and phosphofructokinase had high flux control coefficients, while for the ATPase a small negative flux control coefficient was predicted. In a broader comparative analysis of glycolytic models, we identified a weakness in the P. falciparum pathway design with respect to stability towards perturbations in the ATP demand.

DATABASE

The mathematical model described here has been submitted to the JWS Online Cellular Systems Modelling Database and can be accessed at http://jjj.bio.vu.nl/database/vanniekerk1. The SEEK-study including the experimental data set is available at DOI 10.15490/seek.1.

INVESTIGATION

56 (http://dx.doi.org/10.15490/seek.1.

INVESTIGATION

56).

Tracing regulatory routes in metabolism using generalised supply-demand analysis

Background

Generalised supply-demand analysis is a conceptual framework that views metabolism as a molecular economy. Metabolic pathways are partitioned into so-called supply and demand blocks that produce and consume a particular intermediate metabolite. By studying the response of these reaction blocks to perturbations in the concentration of the linking metabolite, different regulatory routes of interaction between the metabolite and its supply and demand blocks can be identified and their contribution quantified. These responses are mediated not only through direct substrate/product interactions, but also through allosteric effects. Here we subject previously published kinetic models of pyruvate metabolism in Lactococcus lactis and aspartate-derived amino acid synthesis in Arabidopsis thaliana to generalised supply-demand analysis.

Results

Multiple routes of regulation are brought about by different mechanisms in each model, leading to behavioural and regulatory patterns that are generally difficult to predict from simple inspection of the reaction networks depicting the models. In the pyruvate model the moiety-conserved cycles of ATP/ADP and NADH/NAD ⁺ allow otherwise independent metabolic branches to communicate. This causes the flux of one ATP-producing reaction block to increase in response to an increasing ATP/ADP ratio, while an NADH-consuming block flux decreases in response to an increasing NADH/NAD ⁺ ratio for certain ratio value ranges.

In the aspartate model, aspartate semialdehyde can inhibit its supply block directly or by increasing the concentration of two amino acids (Lys and Thr) that occur as intermediates in demand blocks and act as allosteric inhibitors of isoenzymes in the supply block. These different routes of interaction from aspartate semialdehyde are each seen to contribute differently to the regulation of the aspartate semialdehyde supply block.

Conclusions

Indirect routes of regulation between a metabolic intermediate and a reaction block that either produces or consumes this intermediate can play a much larger regulatory role than routes mediated through direct interactions. These indirect routes of regulation can also result in counter-intuitive metabolic behaviour. Performing generalised supply-demand analysis on two previously published models demonstrated the utility of this method as an entry point in the analysis of metabolic behaviour and the potential for obtaining novel results from previously analysed models by using new approaches.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-015-0236-1) contains supplementary material, which is available to authorized users.

Rites of passage: requirements and standards for building kinetic models of metabolic phenotypes

The overarching ambition of kinetic metabolic modeling is to capture the dynamic behavior of metabolism to such an extent that systems and synthetic biology strategies can reliably be tested in silico. The lack of kinetic data hampers the development of kinetic models, and most of the current models use ad hoc reduced stoichiometry or oversimplified kinetic rate expressions, which may limit their predictive strength. There is a need to introduce the community-level standards that will organize and accelerate the future developments in this area. We introduce here a set of requirements that will ensure the model quality, we examine the current kinetic models with respect to these requirements, and we propose a general workflow for constructing models that satisfy these requirements.

Microbial production of lactic acid

Lactic acid is an important commodity chemical having a wide range of applications. Microbial production effectively competes with chemical synthesis methods because biochemical synthesis permits the generation of either one of the two enantiomers with high optical purity at high yield and titer, a result which is particularly beneficial for the production of poly(lactic acid) polymers having specific properties. The commercial viability of microbial lactic acid production relies on utilization of inexpensive carbon substrates derived from agricultural or waste resources. Therefore, optimal lactic acid formation requires an understanding and engineering of both the competing pathways involved in carbohydrate metabolism, as well as pathways leading to potential by-products which both affect product yield. Recent research leverages those biochemical pathways, while researchers also continue to seek strains with improved tolerance and ability to perform under desirable industrial conditions, for example, of pH and temperature.

Monte-Carlo modeling of the central carbon metabolism of Lactococcus lactis: insights into metabolic regulation

Metabolic pathways are complex dynamic systems whose response to perturbations and environmental challenges are governed by multiple interdependencies between enzyme properties, reactions rates, and substrate levels. Understanding the dynamics arising from such a network can be greatly enhanced by the construction of a computational model that embodies the properties of the respective system. Such models aim to incorporate mechanistic details of cellular interactions to mimic the temporal behavior of the biochemical reaction system and usually require substantial knowledge of kinetic parameters to allow meaningful conclusions. Several approaches have been suggested to overcome the severe data requirements of kinetic modeling, including the use of approximative kinetics and Monte-Carlo sampling of reaction parameters. In this work, we employ a probabilistic approach to study the response of a complex metabolic system, the central metabolism of the lactic acid bacterium Lactococcus lactis, subject to perturbations and brief periods of starvation. Supplementing existing methodologies, we show that it is possible to acquire a detailed understanding of the control properties of a corresponding metabolic pathway model that is directly based on experimental observations. In particular, we delineate the role of enzymatic regulation to maintain metabolic stability and metabolic recovery after periods of starvation. It is shown that the feedforward activation of the pyruvate kinase by fructose-1,6-bisphosphate qualitatively alters the bifurcation structure of the corresponding pathway model, indicating a crucial role of enzymatic regulation to prevent metabolic collapse for low external concentrations of glucose. We argue that similar probabilistic methodologies will help our understanding of dynamic properties of small-, medium- and large-scale metabolic networks models.

Streptococcus pyogenes polymyxin B-resistant mutants display enhanced ExPortal integrity

The ExPortal protein secretion organelle in Streptococcus pyogenes is an anionic phospholipid-containing membrane microdomain enriched in Sec translocons and postsecretion protein biogenesis factors. Polymyxin B binds to and disrupts ExPortal integrity, resulting in defective secretion of several toxins. To gain insight into factors that influence ExPortal organization, a genetic screen was conducted to select for spontaneous polymyxin B-resistant mutants displaying enhanced ExPortal integrity. Whole-genome resequencing of 25 resistant mutants revealed from one to four mutations per mutant genome clustered primarily within a core set of 10 gene groups. Construction of mutants with individual deletions or insertions demonstrated that 7 core genes confer resistance and enhanced ExPortal integrity through loss of function, while 3 were likely due to gain of function and/or combinatorial effects. Core resistance genes include a transcriptional regulator of lipid biosynthesis, several genes involved in nutrient acquisition, and a variety of genes involved in stress responses. Two members of the latter class also function as novel regulators of the secreted SpeB cysteine protease. Analysis of the most frequently isolated mutation, a single nucleotide deletion in a track of 9 consecutive adenine residues in pstS, encoding a component of a high-affinity Pi transporter, suggests that this sequence functions as a molecular switch to facilitate stress adaptation. Together, these data suggest the existence of a membrane stress response that promotes enhanced ExPortal integrity and resistance to cationic antimicrobial peptides.

Lost in transition: start-up of glycolysis yields subpopulations of nongrowing cells

Cells need to adapt to dynamic environments. Yeast that fail to cope with dynamic changes in the abundance of glucose can undergo growth arrest. We show that this failure is caused by imbalanced reactions in glycolysis, the essential pathway in energy metabolism in most organisms. The imbalance arises largely from the fundamental design of glycolysis, making this state of glycolysis a generic risk. Cells with unbalanced glycolysis coexisted with vital cells. Spontaneous, nongenetic metabolic variability among individual cells determines which state is reached and, consequently, which cells survive. Transient ATP (adenosine 5'-triphosphate) hydrolysis through futile cycling reduces the probability of reaching the imbalanced state. Our results reveal dynamic behavior of glycolysis and indicate that cell fate can be determined by heterogeneity purely at the metabolic level.

An extended dynamic model of Lactococcus lactis metabolism for mannitol and 2,3-butanediol production

Biomedical research and biotechnological production are greatly benefiting from the results provided by the development of dynamic models of microbial metabolism. Although several kinetic models of Lactococcus lactis (a Lactic Acid Bacterium (LAB) commonly used in the dairy industry) have been developed so far, most of them are simplified and focus only on specific metabolic pathways. Therefore, the application of mathematical models in the design of an engineering strategy for the production of industrially important products by L. lactis has been very limited. In this work, we extend the existing kinetic model of L. lactis central metabolism to include industrially relevant production pathways such as mannitol and 2,3-butanediol. In this way, we expect to study the dynamics of metabolite production and make predictive simulations in L. lactis. We used a system of ordinary differential equations (ODEs) with approximate Michaelis-Menten-like kinetics for each reaction, where the parameters were estimated from multivariate time-series metabolite concentrations obtained by our team through in vivo Nuclear Magnetic Resonance (NMR). The results show that the model captures observed transient dynamics when validated under a wide range of experimental conditions. Furthermore, we analyzed the model using global perturbations, which corroborate experimental evidence about metabolic responses upon enzymatic changes. These include that mannitol production is very sensitive to lactate dehydrogenase (LDH) in the wild type (W.T.) strain, and to mannitol phosphoenolpyruvate: a phosphotransferase system (PTS(Mtl)) in a LDH mutant strain. LDH reduction has also a positive control on 2,3-butanediol levels. Furthermore, it was found that overproduction of mannitol-1-phosphate dehydrogenase (MPD) in a LDH/PTS(Mtl) deficient strain can increase the mannitol levels. The results show that this model has prediction capability over new experimental conditions and offers promising possibilities to elucidate the effect of alterations in the main metabolism of L. lactis, with application in strain optimization.

Food environments select microorganisms based on selfish energetic behavior

Nutrient richness, and specifically the abundance of mono- and disaccharides that characterize several food matrixes, such as milk and grape juice, has allowed the speciation of lactic acid bacteria and yeasts with a high fermentation capacity instead of energetically favorable respiratory metabolism. In these environmental contexts, rapid sugar consumption and lactic acid or ethanol production, accumulation, and tolerance, together with the ability to propagate in the absence of oxygen, are several of the "winning" traits that have apparently evolved and become specialized to perfection in these fermenting microorganisms. Here, we summarize and discuss the evolutionary context that has driven energetic metabolism in food-associated microorganisms, using the dairy species Lactococcus lactis and Streptococcus thermophilus among prokaryotes and the bakers' yeast Saccharomyces cerevisiae among eukaryotes as model organisms.

Organism-adapted specificity of the allosteric regulation of pyruvate kinase in lactic acid bacteria

Pyruvate kinase (PYK) is a critical allosterically regulated enzyme that links glycolysis, the primary energy metabolism, to cellular metabolism. Lactic acid bacteria rely almost exclusively on glycolysis for their energy production under anaerobic conditions, which reinforces the key role of PYK in their metabolism. These organisms are closely related, but have adapted to a huge variety of native environments. They include food-fermenting organisms, important symbionts in the human gut, and antibiotic-resistant pathogens. In contrast to the rather conserved inhibition of PYK by inorganic phosphate, the activation of PYK shows high variability in the type of activating compound between different lactic acid bacteria. System-wide comparative studies of the metabolism of lactic acid bacteria are required to understand the reasons for the diversity of these closely related microorganisms. These require knowledge of the identities of the enzyme modifiers. Here, we predict potential allosteric activators of PYKs from three lactic acid bacteria which are adapted to different native environments. We used protein structure-based molecular modeling and enzyme kinetic modeling to predict and validate potential activators of PYK. Specifically, we compared the electrostatic potential and the binding of phosphate moieties at the allosteric binding sites, and predicted potential allosteric activators by docking. We then made a kinetic model of Lactococcus lactis PYK to relate the activator predictions to the intracellular sugar-phosphate conditions in lactic acid bacteria. This strategy enabled us to predict fructose 1,6-bisphosphate as the sole activator of the Enterococcus faecalis PYK, and to predict that the PYKs from Streptococcus pyogenes and Lactobacillus plantarum show weaker specificity for their allosteric activators, while still having fructose 1,6-bisphosphate play the main activator role in vivo. These differences in the specificity of allosteric activation may reflect adaptation to different environments with different concentrations of activating compounds. The combined computational approach employed can readily be applied to other enzymes.

Kinetic characterization of arginine deiminase and carbamate kinase from Streptococcus pyogenes M49

Streptococcus pyogenes (group A Streptococcus, GAS) is an important human pathogen causing mild superficial infections of skin and mucous membranes, but also life-threatening systemic diseases. S. pyogenes and other prokaryotic organisms use the arginine deiminase system (ADS) for survival in acidic environments. In this study, the arginine deiminase (AD), and carbamate kinase (CK) from S. pyogenes M49 strain 591 were heterologously expressed in Escherichia coli DH5α, purified, and kinetically characterized. AD and CK from S. pyogenes M49 share high amino acid sequence similarity with the respective enzymes from Lactococcus lactis subsp. lactis IL1403 (45.6% and 53.5% identical amino acids) and Enterococcus faecalis V583 (66.8% and 66.8% identical amino acids). We found that the arginine deiminase of S. pyogenes is not allosterically regulated by the intermediates and products of the arginine degradation (e.g., ATP, citrulline, carbamoyl phosphate). The Km and Vmax values for arginine were 1.13±0.12mM (mean±SD) and 1.51±0.07μmol/min/mg protein. The carbamate kinase is inhibited by ATP but unaffected by arginine and citrulline. The Km and Vmax values for ADP were 0.72±0.08mM and 1.10±0.10μmol/min/mg protein and the Km for carbamoyl phosphate was 0.65±0.07mM. The optimum pH and temperature for both enzymes were 6.5 and 37°C, respectively.

Regulation of the activity of lactate dehydrogenases from four lactic acid bacteria

Despite high similarity in sequence and catalytic properties, the l-lactate dehydrogenases (LDHs) in lactic acid bacteria (LAB) display differences in their regulation that may arise from their adaptation to different habitats. We combined experimental and computational approaches to investigate the effects of fructose 1,6-bisphosphate (FBP), phosphate (Pi), and ionic strength (NaCl concentration) on six LDHs from four LABs studied at pH 6 and pH 7. We found that 1) the extent of activation by FBP (Kact) differs. Lactobacillus plantarum LDH is not regulated by FBP, but the other LDHs are activated with increasing sensitivity in the following order: Enterococcus faecalis LDH2 ≤ Lactococcus lactis LDH2 < E. faecalis LDH1 < L. lactis LDH1 ≤ Streptococcus pyogenes LDH. This trend reflects the electrostatic properties in the allosteric binding site of the LDH enzymes. 2) For L. plantarum, S. pyogenes, and E. faecalis, the effects of Pi are distinguishable from the effect of changing ionic strength by adding NaCl. 3) Addition of Pi inhibits E. faecalis LDH2, whereas in the absence of FBP, Pi is an activator of S. pyogenes LDH, E. faecalis LDH1, and L. lactis LDH1 and LDH2 at pH 6. These effects can be interpreted by considering the computed binding affinities of Pi to the catalytic and allosteric binding sites of the enzymes modeled in protonation states corresponding to pH 6 and pH 7. Overall, the results show a subtle interplay among the effects of Pi, FBP, and pH that results in different regulatory effects on the LDHs of different LABs.

A Systems Biology Study on NFκB Signaling in Primary Mouse Hepatocytes

The cytokine tumor necrosis factor-alpha (TNFα) is one of the key factors during the priming phase of liver regeneration as well as in hepatocarcinogenesis. TNFα activates the nuclear factor κ-light-chain-enhancer of activated B cells (NFκB) signaling pathway and contributes to the conversion of quiescent hepatocytes to activated hepatocytes that are able to proliferate in response to growth factor stimulation. Different mathematical models have been previously established for TNFα/NFκB signaling in the context of tumor cells. Combining these mathematical models with time-resolved measurements of expression and phosphorylation of TNFα/NFκB pathway constituents in primary mouse hepatocytes revealed that an additional phosphorylation step of the NFκB isoform p65 has to be considered in the mathematical model in order to sufficiently describe the dynamics of pathway activation in the primary cells. Also, we addressed the role of basal protein turnover by experimentally measuring the degradation rate of pivotal players in the absence of TNFα and including this information in the model. To elucidate the impact of variations in the protein degradation rates on TNFα/NFκB signaling on the overall dynamic behavior we used global sensitivity analysis that accounts for parameter uncertainties and showed that degradation and translation of p65 had a major impact on the amplitude and the integral of p65 phosphorylation. Finally, our mathematical model of TNFα/NFκB signaling was able to predict the time-course of the complex formation of p65 and of the inhibitor of NFκB (IκB) in primary mouse hepatocytes, which was experimentally verified. Hence, we here present a mathematical model for TNFα/NFκB signaling in primary mouse hepatocytes that provides an important basis to quantitatively disentangle the complex interplay of multiple factors in liver regeneration and tumorigenesis.

Engineering strategies aimed at control of acidification rate of lactic acid bacteria

The ability of lactic acid bacteria to produce lactic acid from various sugars plays an important role in food fermentations. Lactic acid is derived from pyruvate, the end product of glycolysis and thus a fast lactic acid production rate requires a high glycolytic flux. In addition to lactic acid, alternative end products--ethanol, acetic acid and formic acid--are formed by many species. The central role of glycolysis in lactic acid bacteria has provoked numerous studies aiming at identifying potential bottleneck(s) since knowledge about flux control could be important not only for optimizing food fermentation processes, but also for novel applications of lactic acid bacteria, such as cell factories for the production of green fuels and chemicals. With respect to the control and regulation of the fermentation mode, some progress has been made, but the question of which component(s) control the main glycolytic flux remains unanswered.

Glycolytic oscillations in a model of a lactic acid bacterium metabolism

Glycolytic oscillations in yeast have been extensively studied. It is still unclear, if these oscillations are caused by the allosteric enzyme phosphofructokinase or the stoichiometry of glycolysis which contains an autocatalysis with respect to ATP. Bacterial glycolysis shows a different stoichiometry, however, also containing a stoichiometric autocatalysis. For Escherichia coli, the regulation of the enzyme phosphofructokinase is also assumed to be a major reason for oscillations to occur. We investigated glycolytic oscillations in a quantitative kinetic model for Streptococcus pyogenes set-up on the basis of experimental data. We found oscillations within physiologically feasible parameter ranges. We investigated the origin of these oscillations and conclude that, again, both the stoichiometry of the system, as well as its allosterically regulated enzymes can give rise to these oscillations. For the analysis we employed established and new optimization methods for finding oscillatory regimes and present these in the context of this study.

Towards metagenome-scale models for industrial applications--the case of Lactic Acid Bacteria

We review the uses and limitations of modelling approaches that are in use in the field of Lactic Acid Bacteria (LAB). We describe recent developments in model construction and computational methods, starting from application of such models to monocultures. However, since most applications in food biotechnology involve complex nutrient environments and mixed cultures, we extend the scope to discuss developments in modelling such complex systems. With metagenomics and meta-functional genomics data becoming available, the developments in genome-scale community models are discussed. We conclude that exploratory tools are available and useful, but truly predictive mechanistic models will remain a major challenge in the field.

Metabolic shifts: a fitness perspective for microbial cell factories

Performance of industrial microorganisms as cell factories is limited by the capacity to channel nutrients to desired products, of which optimal production usually requires careful manipulation of process conditions, or strain improvement. The focus in process improvement is often on understanding and manipulating the regulation of metabolism. Nonetheless, one encounters situations where organisms are remarkably resilient to further optimization or their properties become unstable. Therefore it is important to understand the origin of these apparent limitations to find whether and how they can be improved. We argue that by considering fitness effects of regulation, a more generic explanation for certain behaviour can be obtained. In this view, apparent process limitations arise from trade-offs that cells faced as they evolved to improve fitness. A deeper understanding of such trade-offs using a systems biology approach can ultimately enhance performance of cell factories.

Effects of the ERES pathogenicity region regulator Ralp3 on Streptococcus pyogenes serotype M49 virulence factor expression

Streptococcus pyogenes (group A streptococcus [GAS]) is a highly virulent Gram-positive bacterium. For successful infection, GAS expresses many virulence factors, which are clustered together with transcriptional regulators in distinct genomic regions. Ralp3 is a central regulator of the ERES region. In this study, we investigated the role of Ralp3 in GAS M49 pathogenesis. The inactivation of Ralp3 resulted in reduced attachment to and internalization into human keratinocytes. The Δralp3 mutant failed to survive in human blood and serum, and the hyaluronic acid capsule was slightly decreased. In addition, the mutant showed a lower binding capacity to human plasminogen, and the SpeB activity was significantly decreased. Complementation of the Δralp3 mutant restored the wild-type phenotype. The transcriptome and quantitative reverse transcription-PCR analysis of the serotype M49 GAS strain and its isogenic Δralp3 mutant identified 16 genes as upregulated, and 43 genes were found to be downregulated. Among the downregulated genes, there were open reading frames encoding proteins involved in metabolism (e.g., both lac operons and the fru operon), genes encoding lantibiotics (e.g., the putative salivaricin operon), and ORFs encoding virulence factors (such as the whole Mga core regulon and further genes under Mga control). In summary, the ERES region regulator Ralp3 is an important serotype-specific transcriptional regulator for virulence and metabolic control.

Systems solutions by lactic acid bacteria: from paradigms to practice

Lactic acid bacteria are among the powerhouses of the food industry, colonize the surfaces of plants and animals, and contribute to our health and well-being. The genomic characterization of LAB has rocketed and presently over 100 complete or nearly complete genomes are available, many of which serve as scientific paradigms. Moreover, functional and comparative metagenomic studies are taking off and provide a wealth of insight in the activity of lactic acid bacteria used in a variety of applications, ranging from starters in complex fermentations to their marketing as probiotics. In this new era of high throughput analysis, biology has become big science. Hence, there is a need to systematically store the generated information, apply this in an intelligent way, and provide modalities for constructing self-learning systems that can be used for future improvements. This review addresses these systems solutions with a state of the art overview of the present paradigms that relate to the use of lactic acid bacteria in industrial applications. Moreover, an outlook is presented of the future developments that include the transition into practice as well as the use of lactic acid bacteria in synthetic biology and other next generation applications.