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.

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.

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.

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.

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.

Accurate fidelity analysis of the reverse transcriptase by a modified next-generation sequencing.

We evaluated fidelity of various reverse transcriptases (RTs) by a novel method with modified next-generation sequencing (NGS). In the optimized condition, one NGS run could handle cDNA products from multiple cDNA synthesis reactions performed at different conditions. This was achieved using a primer containing not only the tag of 14 randomized bases to label each cDNA molecule but also a tag of five bases to label each reaction condition. With this method, we quantitated the error rates of 44 cDNA synthesis reactions by retroviral RTs or genetically engineered DNA polymerases with RT activity under different conditions. The results indicated that high concentrations of MgCl2, Mn(OCOCH3)2, and dNTP decrease the fidelity and that these effects are more pronounced in reactions using RT from human immunodeficiency virus type 1. This is the first report about a precise fidelity monitoring of various RTs by a direct sequence determination.

Slowdown of Translational Elongation in Escherichia coli under Hyperosmotic Stress.

In nature, bacteria frequently experience many adverse conditions, including heat, oxidation, acidity, and hyperosmolarity, which all tend to slow down if not outright stop cell growth. Previous work on bacterial stress mainly focused on understanding gene regulatory responses. Much less is known about how stresses compromise protein synthesis, which is the major driver of cell growth. Here, we quantitatively characterize the translational capacity of Escherichia coli cells growing exponentially under hyperosmotic stress. We found that hyperosmotic stress affects bacterial protein synthesis through reduction of the translational elongation rate, which is largely compensated for by an increase in the cellular ribosome content compared with nutrient limitation at a similar growth rate. The slowdown of translational elongation is attributed to a reduction in the rate of binding of tRNA ternary complexes to the ribosomes.IMPORTANCE Hyperosmotic stress is a common stress condition confronted by E. coli during infection of the urinary tract. It can significantly compromise the bacterial growth rate. Protein translation capacity is a critical component of bacterial growth. In this study, we find for the first time that hyperosmotic stress causes substantial slowdown in bacterial ribosome translation elongation. The slowdown of translation elongation originates from a reduced binding rate of tRNA ternary complex to the ribosomes.

High sensitive RNA detection by one-step RT-PCR using the genetically engineered variant of DNA polymerase with reverse transcriptase activity from hyperthermophilies.

One-step RT-PCR has not been widely used even though some thermostable DNA polymerases with reverse transcriptase (RT) activity were developed from bacterial and archaeal polymerases, which is owing to low cDNA synthesis activity from RNA. In the present study, we developed highly-sensitive one-step RT-PCR using the single variant of family A DNA polymerase with RT activity, K4polL329A (L329A), from the hyperthermophilic bacterium Thermotoga petrophila K4 or the 16-tuple variant of family B DNA polymerase with RT activity, RTX, from the hyperthermophilic archaeon Thermococcus kodakarensis. Optimization of reaction condition revealed that the activities for cDNA synthesis and PCR of K4polL329A and RTX were highly affected by the concentrations of MgCl2 and Mn(OCOCH3)2 as well as those of K4polL329A or RTX. Under the optimized condition, 300 copies/µl of target RNA in 10 µl reaction volumes were successfully detected by the one-step RT-PCR with K4polL329A or RTX, which was almost equally sensitive enough compared with the current RT-PCR condition using retroviral RT and thermostable DNA polymerase. Considering that K4polL329A and RTX are stable even at 90-100°C, our results suggest that the one-step RT-PCR with K4polL329A or RTX is more advantageous than the current one.

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.

Enhanced detection of RNA by MMLV reverse transcriptase coupled with thermostable DNA polymerase and DNA/RNA helicase.

Detection of mRNA is a valuable method for monitoring the specific gene expression. In this study, we devised a novel cDNA synthesis method using three enzymes, the genetically engineered thermostable variant of reverse transcriptase (RT), MM4 (E286R/E302K/L435R/D524A) from Moloney murine leukemia virus (MMLV), the genetically engineered variant of family A DNA polymerase with RT activity, K4polL329A from thermophilic Thermotoga petrophila K4, and the DNA/RNA helicase Tk-EshA from a hyperthermophilic archaeon Thermococcus kodakarensis. By optimizing assay conditions for three enzymes using Taguchi's method, 100 to 1000-fold higher sensitivity was achieved for cDNA synthesis than conventional assay condition using only RT. Our results suggest that DNA polymerase with RT activity and DNA/RNA helicase are useful to increase the sensitivity of cDNA synthesis.

Reduction of translating ribosomes enables Escherichia coli to maintain elongation rates during slow growth.

Bacteria growing under different conditions experience a broad range of demand on the rate of protein synthesis, which profoundly affects cellular resource allocation. During fast growth, protein synthesis has long been known to be modulated by adjusting the ribosome content, with the vast majority of ribosomes engaged at a near-maximal rate of elongation. Here, we systematically characterize protein synthesis by Escherichia coli, focusing on slow-growth conditions. We establish that the translational elongation rate decreases as growth slows, exhibiting a Michaelis-Menten dependence on the abundance of the cellular translational apparatus. However, an appreciable elongation rate is maintained even towards zero growth, including the stationary phase. This maintenance, critical for timely protein synthesis in harsh environments, is accompanied by a drastic reduction in the fraction of active ribosomes. Interestingly, well-known antibiotics such as chloramphenicol also cause a substantial reduction in the pool of active ribosomes, instead of slowing down translational elongation as commonly thought.

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.

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.

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.