Four tables available at the above site.

The kinetic parameters of enzymes are key to understanding the rate and specificity of most biological processes. Although specific trends are frequently studied for individual enzymes, global trends are rarely addressed. We performed an analysis of k(cat) and K(M) values of several thousand enzymes collected from the literature. We found that the "average enzyme" exhibits a k(cat) of ~0 s(-1) and a k(cat)/K(M) of ~10(5) s(-1) M(-1), much below the diffusion limit and the characteristic textbook portrayal of kinetically superior enzymes. Why do most enzymes exhibit moderate catalytic efficiencies? Maximal rates may not evolve in cases where weaker selection pressures are expected. We find, for example, that enzymes operating in secondary metabolism are, on average, ~30-fold slower than those of central metabolism. We also find indications that the physicochemical properties of substrates affect the kinetic parameters. Specifically, low molecular mass and hydrophobicity appear to limit K(M) optimization. In accordance, substitution with phosphate, CoA, or other large modifiers considerably lowers the K(M) values of enzymes utilizing the substituted substrates. It therefore appears that both evolutionary selection pressures and physicochemical constraints shape the kinetic parameters of enzymes. It also seems likely that the catalytic efficiency of some enzymes toward their natural substrates could be increased in many cases by natural or laboratory evolution.

Noise in gene expression is generated at multiple levels, such as transcription and translation, chromatin remodeling and pathway-specific regulation. Studies of individual promoters have suggested different dominating noise sources, raising the question of whether a general trend exists across a large number of genes and conditions. We examined the variation in the expression levels of 43 Saccharomyces cerevisiae proteins, in cells grown under 11 experimental conditions. For all classes of genes and under all conditions, the expression variance was approximately proportional to the mean; the same scaling was observed at steady state and during the transient responses to the perturbations. Theoretical analysis suggests that this scaling behavior reflects variability in mRNA copy number, resulting from random 'birth and death' of mRNA molecules or from promoter fluctuations. Deviation of coexpressed genes from this general trend, including high noise in stress-related genes and low noise in proteasomal genes, may indicate fluctuations in pathway-specific regulators or a differential activation pattern of the underlying gene promoters.

The laws of thermodynamics constrain the action of biochemical systems. However, thermodynamic data on biochemical compounds can be difficult to find and is cumbersome to perform calculations with manually. Even simple thermodynamic questions like ‘how much Gibbs energy is released by ATP hydrolysis at pH 5?’ are complicated excessively by the search for accurate data. To address this problem, eQuilibrator couples a comprehensive and accurate database of thermodynamic properties of biochemical compounds and reactions with a simple and powerful online search and calculation interface. The web interface to eQuilibrator (http://equilibrator.weizmann.ac.il) enables easy calculation of Gibbs energies of compounds and reactions given arbitrary pH, ionic strength and metabolite concentrations. The eQuilibrator code is open-source and all thermodynamic source data are freely downloadable in standard formats. Here we describe the database characteristics and implementation and demonstrate its use.

Contrary to the textbook portrayal of glycolysis as a single pathway conserved across all domains of life, not all sugar-consuming organisms use the canonical Embden-Meyerhoff-Parnass (EMP) glycolytic pathway. Prokaryotic glucose metabolism is particularly diverse, including several alternative glycolytic pathways, the most common of which is the Entner-Doudoroff (ED) pathway. The prevalence of the ED pathway is puzzling as it produces only one ATP per glucose-half as much as the EMP pathway. We argue that the diversity of prokaryotic glucose metabolism may reflect a tradeoff between a pathway's energy (ATP) yield and the amount of enzymatic protein required to catalyze pathway flux. We introduce methods for analyzing pathways in terms of thermodynamics and kinetics and show that the ED pathway is expected to require several-fold less enzymatic protein to achieve the same glucose conversion rate as the EMP pathway. Through genomic analysis, we further show that prokaryotes use different glycolytic pathways depending on their energy supply. Specifically, energy-deprived anaerobes overwhelmingly rely upon the higher ATP yield of the EMP pathway, whereas the ED pathway is common among facultative anaerobes and even more common among aerobes. In addition to demonstrating how protein costs can explain the use of alternative metabolic strategies, this study illustrates a direct connection between an organism's environment and the thermodynamic and biochemical properties of the metabolic pathways it employs.evolution | enzyme cost G lycolysis is the process by which glucose is broken down anaerobically into incompletely oxidized compounds like pyruvate, a process which is usually coupled to the synthesis of ATP. Although the Embden-Meyerhof-Parnas pathway (EMP, often simply "glycolysis") is the nearly ubiquitous glycolytic route among eukaryotes (1, 2), it is not the only game in town. Prokaryotes display impressive diversity in glucose metabolism (2, 3) and natural glycolytic alternatives like the Entner-Doudoroff (ED), and phosphoketolase pathways attest to the fact that there are multiple biologically feasible routes for glucose metabolism (2, 4-8). Natural glycolytic pathways vary in their reaction sequence and in how much ATP they produce per glucose metabolized, ranging from zero to three ATP molecules in most cases (7).The EMP and ED pathways ( Fig. 1 A and B and Fig. S1) are the most common bacterial glycolytic pathways (2, 4, 9), and their general schemes are quite similar: glucose is phosphorylated and then cleaved into two three-carbon units, which are further metabolized to produce ATP (4). In some organisms, these pathways differ slightly in the specific redox cofactors they use (e.g., NAD + vs. NADP + ; Fig. 1B, Fig. S2, and Table S1), but here we focus on the prominent difference in ATP yield. If we take lactate as a representative final product, these two pathways have the same net reaction:and differ primarily in n, the number of ATP produced, and the specific intermediate reaction steps (Fig. 1B, SI T...

Carbon fixation is the process by which CO 2 is incorporated into organic compounds. In modern agriculture in which water, light, and nutrients can be abundant, carbon fixation could become a significant growth-limiting factor. Hence, increasing the fixation rate is of major importance in the road toward sustainability in food and energy production. There have been recent attempts to improve the rate and specificity of Rubisco, the carboxylating enzyme operating in the Calvin-Benson cycle; however, they have achieved only limited success. Nature employs several alternative carbon fixation pathways, which prompted us to ask whether more efficient novel synthetic cycles could be devised. Using the entire repertoire of approximately 5,000 metabolic enzymes known to occur in nature, we computationally identified alternative carbon fixation pathways that combine existing metabolic building blocks from various organisms. We compared the natural and synthetic pathways based on physicochemical criteria that include kinetics, energetics, and topology. Our study suggests that some of the proposed synthetic pathways could have significant quantitative advantages over their natural counterparts, such as the overall kinetic rate. One such cycle, which is predicted to be two to three times faster than the Calvin-Benson cycle, employs the most effective carboxylating enzyme, phosphoenolpyruvate carboxylase, using the core of the naturally evolved C4 cycle. Although implementing such alternative cycles presents daunting challenges related to expression levels, activity, stability, localization, and regulation, we believe our findings suggest exciting avenues of exploration in the grand challenge of enhancing food and renewable fuel production via metabolic engineering and synthetic biology. metabolic engineering | synthetic biology | photosynthesis | carboxylation | biological optimization I n the process of transforming sunlight into stored chemical energy, plants absorb approximately 10 times more CO 2 from the atmosphere than the total amount emitted by human activities globally (1). Moreover, agriculture, which is effectively a massive carbon fixation industry, makes use of the majority of cultivatable land on earth and accounts for most of the fresh water used by humanity (2). These figures point to the central role that carbon fixation by plants plays in our global ecological footprint. In nature, the factors limiting the growth of photosynthetic organisms vary among habitats and often include the availability of water, light, fixed nitrogen, iron, and phosphorus (e.g., ref.3). However, in agriculture today, the use of fertilizers and irrigation can make carbon fixation a rate-limiting factor. For example, many C3 plants have shown a significant increase in biomass when exposed to twice the atmospheric CO 2 concentration (4). Impressive growth enhancements have been demonstrated by addressing several biochemical limiting factors, related both to the light-dependent and light-independent reactions (5-9). For instance, tran...

Turnover numbers, also known as k cat values, are fundamental properties of enzymes. However, k cat data are scarce and measured in vitro, thus may not faithfully represent the in vivo situation. A basic question that awaits elucidation is: how representative are k cat values for the maximal catalytic rates of enzymes in vivo? Here, we harness omics data to calculate k vivo max , the observed maximal catalytic rate of an enzyme inside cells. Comparison with k cat values from Escherichia coli, yields a correlation of r 2 = 0.62 in log scale (p < 10 −10 ), with a root mean square difference of 0.54 (3.5-fold in linear scale), indicating that in vivo and in vitro maximal rates generally concur. By accounting for the degree of saturation of enzymes and the backward flux dictated by thermodynamics, we further refine the correspondence between k vivo max and k cat values. The approach we present here characterizes the quantitative relationship between enzymatic catalysis in vitro and in vivo and offers a highthroughput method for extracting enzyme kinetic constants from omics data. (1-6). Many models of cellular metabolism include k cat values, the maximal turnover rates of enzymes, as key inputs to predict the behavior of metabolic pathways and networks (7-9). However, most values have never been measured experimentally. Escherichia coli is the most intensely biochemically characterized organism, but k cat values are available for only about 10% of its ≈ 2, 000 enzyme-reaction pairs (Dataset S1). Indeed, k cat values are missing for several central metabolic enzymes. The scarcity of kinetic data limits the scope of models and necessitates generic parameter assignments that significantly reduce the predictive power of cellular models.Even if a larger collection of k cat values was made available, their current use poses a major difficulty: k cat values are measured through in vitro enzyme assays, representing the initial rate of the reaction, i.e., full substrates saturation and negligible levels of products. Such assays may underrepresent factors like cellular metabolite concentrations, thermodynamic constraints, posttranslational modifications, chaperones, cellular crowding, and activating and inhibiting molecules, which can substantially alter enzyme kinetics in vivo. These omissions call into question the relevance of k cat measurements in vivo (10-12). Furthermore, an effort to measure a large number of k cat values under in vivo-like conditions presents a daunting challenge, given how many unknown biochemical factors might be involved.Several studies grapple with missing k cat values by sampling from the distribution of k cat values measured in vitro or by using measurements of the same enzyme from related species (13-16). These approximations systematically ignore any errors resulting from the differences between in vitro and in vivo environments. Approximations of this sort may also introduce significant errors, as k cat values can deviate by orders of magnitude between isozymes in the same organism as well a...

In metabolism research, thermodynamics is usually used to determine the directionality of a reaction or the feasibility of a pathway. However, the relationship between thermodynamic potentials and fluxes is not limited to questions of directionality: thermodynamics also affects the kinetics of reactions through the flux-force relationship, which states that the logarithm of the ratio between the forward and reverse fluxes is directly proportional to the change in Gibbs energy due to a reaction (ΔrG′). Accordingly, if an enzyme catalyzes a reaction with a ΔrG′ of -5.7 kJ/mol then the forward flux will be roughly ten times the reverse flux. As ΔrG′ approaches equilibrium (ΔrG′ = 0 kJ/mol), exponentially more enzyme counterproductively catalyzes the reverse reaction, reducing the net rate at which the reaction proceeds. Thus, the enzyme level required to achieve a given flux increases dramatically near equilibrium. Here, we develop a framework for quantifying the degree to which pathways suffer these thermodynamic limitations on flux. For each pathway, we calculate a single thermodynamically-derived metric (the Max-min Driving Force, MDF), which enables objective ranking of pathways by the degree to which their flux is constrained by low thermodynamic driving force. Our framework accounts for the effect of pH, ionic strength and metabolite concentration ranges and allows us to quantify how alterations to the pathway structure affect the pathway's thermodynamics. Applying this methodology to pathways of central metabolism sheds light on some of their features, including metabolic bypasses (e.g., fermentation pathways bypassing substrate-level phosphorylation), substrate channeling (e.g., of oxaloacetate from malate dehydrogenase to citrate synthase), and use of alternative cofactors (e.g., quinone as an electron acceptor instead of NAD). The methods presented here place another arrow in metabolic engineers' quiver, providing a simple means of evaluating the thermodynamic and kinetic quality of different pathway chemistries that produce the same molecules.