A universal trade-off between growth and lag in fluctuating environments.

The rate of cell growth is crucial for bacterial fitness and drives the allocation of bacterial resources, affecting, for example, the expression levels of proteins dedicated to metabolism and biosynthesis1,2. It is unclear, however, what ultimately determines growth rates in different environmental conditions. Moreover, increasing evidence suggests that other objectives are also important3-7, such as the rate of physiological adaptation to changing environments8,9. A common challenge for cells is that these objectives cannot be independently optimized, and maximizing one often reduces another. Many such trade-offs have indeed been hypothesized on the basis of qualitative correlative studies8-11. Here we report a trade-off between steady-state growth rate and physiological adaptability in Escherichia coli, observed when a growing culture is abruptly shifted from a preferred carbon source such as glucose to fermentation products such as acetate. These metabolic transitions, common for enteric bacteria, are often accompanied by multi-hour lags before growth resumes. Metabolomic analysis reveals that long lags result from the depletion of key metabolites that follows the sudden reversal in the central carbon flux owing to the imposed nutrient shifts. A model of sequential flux limitation not only explains the observed trade-off between growth and adaptability, but also allows quantitative predictions regarding the universal occurrence of such tradeoffs, based on the opposing enzyme requirements of glycolysis versus gluconeogenesis. We validate these predictions experimentally for many different nutrient shifts in E. coli, as well as for other respiro-fermentative microorganisms, including Bacillus subtilis and Saccharomyces cerevisiae.

Global coordination of metabolic pathways in Escherichia coli by active and passive regulation.

Microorganisms adjust metabolic activity to cope with diverse environments. While many studies have provided insights into how individual pathways are regulated, the mechanisms that give rise to coordinated metabolic responses are poorly understood. Here, we identify the regulatory mechanisms that coordinate catabolism and anabolism in Escherichia coli. Integrating protein, metabolite, and flux changes in genetically implemented catabolic or anabolic limitations, we show that combined global and local mechanisms coordinate the response to metabolic limitations. To allocate proteomic resources between catabolism and anabolism, E. coli uses a simple global gene regulatory program. Surprisingly, this program is largely implemented by a single transcription factor, Crp, which directly activates the expression of catabolic enzymes and indirectly reduces the expression of anabolic enzymes by passively sequestering cellular resources needed for their synthesis. However, metabolic fluxes are not controlled by this regulatory program alone; instead, fluxes are adjusted mostly through passive changes in the local metabolite concentrations. These mechanisms constitute a simple but effective global regulatory program that coarsely partitions resources between different parts of metabolism while ensuring robust coordination of individual metabolic reactions.

From coarse to fine: the absolute Escherichia coli proteome under diverse growth conditions.

Accurate measurements of cellular protein concentrations are invaluable to quantitative studies of gene expression and physiology in living cells. Here, we developed a versatile mass spectrometric workflow based on data-independent acquisition proteomics (DIA/SWATH) together with a novel protein inference algorithm (xTop). We used this workflow to accurately quantify absolute protein abundances in Escherichia coli for > 2,000 proteins over > 60 growth conditions, including nutrient limitations, non-metabolic stresses, and non-planktonic states. The resulting high-quality dataset of protein mass fractions allowed us to characterize proteome responses from a coarse (groups of related proteins) to a fine (individual) protein level. Hereby, a plethora of novel biological findings could be elucidated, including the generic upregulation of low-abundant proteins under various metabolic limitations, the non-specificity of catabolic enzymes upregulated under carbon limitation, the lack of large-scale proteome reallocation under stress compared to nutrient limitations, as well as surprising strain-dependent effects important for biofilm formation. These results present valuable resources for the systems biology community and can be used for future multi-omics studies of gene regulation and metabolic control in E. coli.

Overflow metabolism in Escherichia coli results from efficient proteome allocation.

Overflow metabolism refers to the seemingly wasteful strategy in which cells use fermentation instead of the more efficient respiration to generate energy, despite the availability of oxygen. Known as the Warburg effect in the context of cancer growth, this phenomenon occurs ubiquitously for fast-growing cells, including bacteria, fungi and mammalian cells, but its origin has remained unclear despite decades of research. Here we study metabolic overflow in Escherichia coli, and show that it is a global physiological response used to cope with changing proteomic demands of energy biogenesis and biomass synthesis under different growth conditions. A simple model of proteomic resource allocation can quantitatively account for all of the observed behaviours, and accurately predict responses to new perturbations. The key hypothesis of the model, that the proteome cost of energy biogenesis by respiration exceeds that by fermentation, is quantitatively confirmed by direct measurement of protein abundances via quantitative mass spectrometry.

Coordination of bacterial proteome with metabolism by cyclic AMP signalling.

The cyclic AMP (cAMP)-dependent catabolite repression effect in Escherichia coli is among the most intensely studied regulatory processes in biology. However, the physiological function(s) of cAMP signalling and its molecular triggers remain elusive. Here we use a quantitative physiological approach to show that cAMP signalling tightly coordinates the expression of catabolic proteins with biosynthetic and ribosomal proteins, in accordance with the cellular metabolic needs during exponential growth. The expression of carbon catabolic genes increased linearly with decreasing growth rates upon limitation of carbon influx, but decreased linearly with decreasing growth rate upon limitation of nitrogen or sulphur influx. In contrast, the expression of biosynthetic genes showed the opposite linear growth-rate dependence as the catabolic genes. A coarse-grained mathematical model provides a quantitative framework for understanding and predicting gene expression responses to catabolic and anabolic limitations. A scheme of integral feedback control featuring the inhibition of cAMP signalling by metabolic precursors is proposed and validated. These results reveal a key physiological role of cAMP-dependent catabolite repression: to ensure that proteomic resources are spent on distinct metabolic sectors as needed in different nutrient environments. Our findings underscore the power of quantitative physiology in unravelling the underlying functions of complex molecular signalling networks.

Structural Basis for the Serratia marcescens Lipase Secretion System: Crystal Structures of the Membrane Fusion Protein and Nucleotide-Binding Domain.

Serratia marcescens secretes a lipase, LipA, through a type I secretion system (T1SS). The T1SS for LipA, the Lip system, is composed of an inner membrane ABC transporter with its nucleotide-binding domains (NBD), LipB, a membrane fusion protein, LipC, and an outer membrane channel protein, LipD. Passenger protein secreted by this system has been functionally and structurally characterized well, but relatively little information about the transporter complex is available. Here, we report the crystallographic studies of LipC without the membrane anchor region, LipC-, and the NBD of LipB (LipB-NBD). LipC- crystallographic analysis has led to the determination of the structure of the long α-helical and lipoyl domains, but not the area where it interacts with LipB, suggesting that the region is flexible without LipB. The long α-helical domain has three α-helices, which interacts with LipD in the periplasm. LipB-NBD has the common overall architecture and ATP hydrolysis activity of ABC transporter NBDs. Using the predicted models of full-length LipB and LipD, the overall structural insight into the Lip system is discussed.

Next-generation sequencing-based analysis of reverse transcriptase fidelity.

In this study, we devised a simple and rapid method to analyze fidelity of reverse transcriptase (RT) using next-generation sequencing (NGS). The method comprises a cDNA synthesis reaction from standard RNA with a primer containing a tag of 14 randomized bases and the RT to be tested, PCR using high-fidelity DNA polymerase, and NGS. By comparing the sequence of each read with the reference sequence, mutations were identified. The mutation can be identified to be due to an error introduced by either cDNA synthesis, PCR, or NGS based on whether the sequence reads with the same tag contain the same mutation or not. The error rates in cDNA synthesis with Moloney murine leukemia virus (MMLV) RT thermostable variant MM4 or the recently developed 16-tuple variant of family B DNA polymerase with RT activity, RTX, from Thermococcus kodakarensis, were 0.75-1.0 × 10-4 errors/base, while that in the reaction with the wild-type human immunodeficiency virus type 1 (HIV-1) RT was 2.6 × 10-4 errors/base. Overall, our method could precisely evaluate the fidelity of various RTs with different reaction conditions in a high-throughput manner without the use of expensive optics and troublesome adaptor ligation.

A growth-rate composition formula for the growth of E.coli on co-utilized carbon substrates.

When bacteria are cultured in medium with multiple carbon substrates, they frequently consume these substrates simultaneously. Building on recent advances in the understanding of metabolic coordination exhibited by Escherichia coli cells through cAMP-Crp signaling, we show that this signaling system responds to the total carbon-uptake flux when substrates are co-utilized and derive a mathematical formula that accurately predicts the resulting growth rate, based only on the growth rates on individual substrates.

Expression of moloney murine leukemia virus reverse transcriptase in a cell-free protein expression system.

OBJECTIVE: To characterize Moloney murine leukemia virus (MMLV) reverse transcriptases (RTs) expressed in a cell-free system and in Escherichia coli. RESULTS: We previously expressed MMLV RT using an E. coli expression system and generated a highly thermostable quadruple variant MM4 (E286R/E302K/L435R/D524A) by site-directed mutagenesis. In this study, we expressed the wild-type MMLV RT (WT) and MM4 using a cell-free protein expression system from insect cells. WT exhibited DNA polymerase and RNase H activities, while MM4, in which the catalytic residue for RNase H activity, Asp524 is changed into Ala, exhibited only DNA polymerase activity. MM4, when held at 60 °C for 10 min, retained DNA polymerase activity, while WT, held at 54 °C for 10 min, lost this activity. In the cDNA synthesis reaction (0.5 µl) in which WT or MM4 were exposed to various temperatures and amounts of target RNA in a microarray chip, MM4 exhibited higher thermostability than WT. CONCLUSION: MMLV RT expressed in the cell-free system is indistinguishable from that expressed in E. coli.

Need-based activation of ammonium uptake in Escherichia coli.

The efficient sequestration of nutrients is vital for the growth and survival of microorganisms. Some nutrients, such as CO2 and NH3, are readily diffusible across the cell membrane. The large membrane permeability of these nutrients obviates the need of transporters when the ambient level is high. When the ambient level is low, however, maintaining a high intracellular nutrient level against passive back diffusion is both challenging and costly. Here, we study the delicate management of ammonium (NH4+/NH3) sequestration by E. coli cells using microfluidic chemostats. We find that as the ambient ammonium concentration is reduced, E. coli cells first maximize their ability to assimilate the gaseous NH3 diffusing into the cytoplasm and then abruptly activate ammonium transport. The onset of transport varies under different growth conditions, but always occurring just as needed to maintain growth. Quantitative modeling of known interactions reveals an integral feedback mechanism by which this need-based uptake strategy is implemented. This novel strategy ensures that the expensive cost of upholding the internal ammonium concentration against back diffusion is kept at a minimum.

Functional decomposition of metabolism allows a system-level quantification of fluxes and protein allocation towards specific metabolic functions.

Quantifying the contribution of individual molecular components to complex cellular processes is a grand challenge in systems biology. Here we establish a general theoretical framework (Functional Decomposition of Metabolism, FDM) to quantify the contribution of every metabolic reaction to metabolic functions, e.g. the synthesis of biomass building blocks. FDM allowed for a detailed quantification of the energy and biosynthesis budget for growing Escherichia coli cells. Surprisingly, the ATP generated during the biosynthesis of building blocks from glucose almost balances the demand from protein synthesis, the largest energy expenditure known for growing cells. This leaves the bulk of the energy generated by fermentation and respiration unaccounted for, thus challenging the common notion that energy is a key growth-limiting resource. Moreover, FDM together with proteomics enables the quantification of enzymes contributing towards each metabolic function, allowing for a first-principle formulation of a coarse-grained model of global protein allocation based on the structure of the metabolic network.

Efficient integrated production of bioethanol and antiviral glycerolysis lignin from sugarcane trash.

BACKGROUND: Sugarcane trash (SCT) represents up to 18% of the aboveground biomass of sugarcane, surpassing 28 million tons globally per year. The majority of SCT is burning in the fields. Hence, efficient use of SCT is necessary to reduce carbon dioxide emissions and global warming and establish agro-industrial biorefineries. Apart from its low costs, conversion of whole biomass with high production efficiency and titer yield is mandatory for effective biorefinery systems. Therefore, in this study, we developed a simple, integrated method involving a single step of glycerolysis pretreatment to produce antiviral glycerolysis lignin (AGL). Subsequently, we co-fermented glycerol with hydrolyzed glucose and xylose to yield high titers of bioethanol. RESULTS: SCT was subjected to pretreatment with microwave acidic glycerolysis with 50% aqueous (aq.) glycerol (MAG50); this pretreatment was optimized across different temperature ranges, acid concentrations, and reaction times. The optimized MAG50 (opMAG50) of SCT at 1:15 (w/v) in 1% H2SO4, 360 µM AlK(SO4)2 at 140 °C for 30 min (opMAG50) recovered the highest amount of total sugars and the lowest amount of furfural byproducts. Following opMAG50, the soluble fraction, i.e., glycerol xylose-rich solution (GXRS), was separated by filtration. A residual pulp was then washed with acetone, recovering 7.9% of the dry weight (27% of lignin) as an AGL. AGL strongly inhibited the replication of encephalomyocarditis virus (EMCV) in L929 cells without cytotoxicity. The pulp was then saccharified in yeast peptone medium by cellulase to produce a glucose concentration similar to the theoretical yield. The total xylose and arabinose recoveries were 69% and 93%, respectively. GXRS and saccharified sugars were combined and co-fermented through mixed cultures of two metabolically engineered Saccharomyces cerevisiae strains: glycerol-fermenting yeast (SK-FGG4) and xylose-fermenting yeast (SK-N2). By co-fermenting glycerol and xylose with glucose, the ethanol titer yield increased to 78.7 g/L (10% v/v ethanol), with a 96% conversion efficiency. CONCLUSION: The integration of AGL production with the co-fermentation of glycerol, hydrolyzed glucose, and xylose to produce a high titer of bioethanol paves an avenue for the use of surplus glycerol from the biodiesel industry for the efficient utilization of SCT and other lignocellulosic biomasses.

A Genome-Scale Metabolic Model of Marine Heterotroph Vibrio splendidus Strain 1A01.

While Vibrio splendidus is best known as an opportunistic pathogen in oysters, Vibrio splendidus strain 1A01 was first identified as an early colonizer of synthetic chitin particles incubated in seawater. To gain a better understanding of its metabolism, a genome-scale metabolic model (GSMM) of V. splendidus 1A01 was reconstructed. GSMMs enable us to simulate all metabolic reactions in a bacterial cell using flux balance analysis. A draft model was built using an automated pipeline from BioCyc. Manual curation was then performed based on experimental data, in part by gap-filling metabolic pathways and tailoring the model's biomass reaction to V. splendidus 1A01. The challenges of building a metabolic model for a marine microorganism like V. splendidus 1A01 are described. IMPORTANCE A genome-scale metabolic model of V. splendidus 1A01 was reconstructed in this work. We offer solutions to the technical problems associated with model reconstruction for a marine bacterial strain like V. splendidus 1A01, which arise largely from the high salt concentration found in both seawater and culture media that simulate seawater.

Enzyme expression kinetics by Escherichia coli during transition from rich to minimal media depends on proteome reserves.

Bacterial fitness depends on adaptability to changing environments. In rich growth medium, which is replete with amino acids, Escherichia coli primarily expresses protein synthesis machineries, which comprise ~40% of cellular proteins and are required for rapid growth. Upon transition to minimal medium, which lacks amino acids, biosynthetic enzymes are synthesized, eventually reaching ~15% of cellular proteins when growth fully resumes. We applied quantitative proteomics to analyse the timing of enzyme expression during such transitions, and established a simple positive relation between the onset time of enzyme synthesis and the fractional enzyme 'reserve' maintained by E. coli while growing in rich media. We devised and validated a coarse-grained kinetic model that quantitatively captures the enzyme recovery kinetics in different pathways, solely on the basis of proteomes immediately preceding the transition and well after its completion. Our model enables us to infer regulatory strategies underlying the 'as-needed' gene expression programme adopted by E. coli.

Hierarchical and simultaneous utilization of carbon substrates: mechanistic insights, physiological roles, and ecological consequences.

Bacteria grown on a mixture of carbon substrates exhibit two utilization patterns: hierarchical utilization (HU) and simultaneous utilization (SU). How and why cells adopt these different behaviors remains poorly understood despite decades of research. Recent studies address various open questions from multiple viewpoints. From a mechanistic perspective, it was found that flux sensors play a central role in the regulation of substrate utilization, accounting for the known dependences on single-substrate growth rates, substrate concentrations, and the point where the substrate enters central metabolism. From a physiological perspective, several recent studies suggested HU or SU as growth-optimizing strategies through efficient allocation of essential proteome resources. However, other studies demonstrate that a significant fraction of the proteome is dedicated to functions apparently unnecessary for growth, casting doubt on explanations based on slight efficiency gains. From an ecological perspective, recent theoretical studies suggest that HU can help increase species diversity in bacterial communities.

Regulation underlying hierarchical and simultaneous utilization of carbon substrates by flux sensors in Escherichia coli.

Many microorganisms exhibit nutrient preferences, exemplified by the 'hierarchical' consumption of certain carbon substrates. Here, we systematically investigate under which physiological conditions hierarchical substrate utilization occurs and its mechanisms of implementation. We show utilization hierarchy of Escherichia coli to be ordered by the carbon-uptake flux rather than the identity of the substrates. A detailed study of glycerol uptake finds that it is fully suppressed if the uptake flux of another glycolytic substrate exceeds a threshold, which is set to the influx obtained when grown on glycerol alone. Below this threshold, limited glycerol uptake is 'supplemented' such that the total carbon uptake is maintained at the threshold. This behaviour results from total-flux feedback mediated by cAMP-Crp signalling but also requires inhibition by the regulator fructose 1,6-bisphosphate, which senses the upper-glycolytic flux and ensures that glycerol uptake defers to other glycolytic substrates but not to gluconeogenic ones. A quantitative model reproduces all of the observed utilization patterns, including those of key mutants. The proposed mechanism relies on the differential regulation of uptake enzymes and requires a specific operon organization. This organization is found to be conserved across related species for several uptake systems, suggesting the deployment of similar mechanisms for hierarchical substrate utilization by a spectrum of microorganisms.

Predictors of gastric myoelectrical activity in type 2 diabetes mellitus.

BACKGROUND: Previous studies have clearly demonstrated the delayed gastric emptying of solid meals in diabetics, whereas their gastric myoelectrical activity, which primarily determines gastric motility, has not yet been fully confirmed. GOALS: This study aimed to clarify the characteristics and potential predictors of gastric myoelectrical activity in type 2 diabetics. STUDY: Twenty-eight diabetics and 18 healthy controls participated. Duodenal biopsy sample was used for reverse transcription-polymerase chain reaction to evaluate cholecystokinin and motilin mRNA contents. Electrogastrography was performed before and after the test meal, and was assessed in terms of dominant frequency; dominant frequency instability coefficient; and the percentage of bradygastria, normogastria, and tachygastria. RESULTS: Over the entire recording period, dominant frequency was significantly lower, and dominant frequency instability coefficient and the percentage of bradygastria were significantly higher in diabetics than in controls. In diabetics, the multiple regression analysis demonstrated that dominant frequency instability coefficient and the percentage of tachygastria in the fasting period were dependent on fasting plasma glucose level and HbA1c, respectively. Moreover, dominant frequency over the entire period and the postprandial percentage of bradygastria were significantly associated with body mass index; the fasting percentage of bradygastria and postprandial dominant frequency instability coefficient were associated with fasting serum leptin level; the postprandial percentage of bradygastria was also associated with cholecystokinin mRNA content. CONCLUSIONS: Gastric myoelectrical activity in type 2 diabetics is impaired on dominant frequency, dominant frequency instability coefficient, and the percentage of bradygastria and predicted by body mass index, fasting serum leptin level, and cholecystokinin mRNA content besides the glycemic status.

Estimation of the source-by-source effect of autorepression on genetic noise.

Transcriptional autorepression has been thought to be one of the simplest control circuits to attenuate fluctuations in gene expression. Here, we explored the effect of autorepression on fluctuations from different noise sources. We theoretically represent the fluctuations in the copy number of proteins as the sum of several terms, each of which is related to a specific noise source and expressed as the product of the source-specific fluctuations under no autorepression (path gain) and the effect of autorepression on them (loop gain). Inspection of each term demonstrates the source-independent noise-attenuating effect of autorepression as well as its source-dependent efficiency. Our experiments using a synthetic autorepression module revealed that autorepression attenuates fluctuations of various noise compositions. These findings indicate that the noise-attenuating effect of autorepression is robust against variation in noise compositions. We also experimentally estimated the loop gain for mRNA noise, demonstrating that loop gains are measurable parameters. Decomposition of fluctuations followed by experimental estimation of path and loop gains would help us to understand the noise-related feature of design principles underlying loop-containing biological networks.

Reconstructing the single-cell-level behavior of a toggle switch from population-level measurements.

Single-cell-level behaviors of cells are typically inferred from ensemble measurements. However, such inferences implicitly assume a biological version of ergodicity: the percentage of cells in a state is identical to the probability to find a cell in that state. While the ergodicity does not always hold, it has been rarely tested. Here, we reveal that the ergodicity does not necessarily hold even for simple toggle switches and that apparent stabilities of the switches are due to a balance between single-cell-level biased stabilities and growth rates differences. Therefore, verification of the ergodicity and reconstructing single-cell-level behaviors are crucial for understanding intracellular systems.