BRENDA, the enzyme information system in 2011

The BRENDA (BRaunschweig ENzyme Database, http://www.brenda-enzymes.org) enzyme information system is the main collection of enzyme functional and property data for the scientific community. The majority of the data are manually extracted from the primary literature. The content covers information on function, structure, occurrence, preparation and application of enzymes as well as properties of mutants and engineered variants. The number of manually annotated references increased by 30% to more than 100 000, the number of ligand structures by 45% to almost 100 000. New query, analysis and data management tools were implemented to improve data processing, data presentation, data input and data access. BRENDA now provides new viewing options such as the display of the statistics of functional parameters and the 3D view of protein sequence and structure features. Furthermore a ligand summary shows comprehensive information on the BRENDA ligands. The enzymes are linked to their respective pathways and can be viewed in pathway maps. The disease text mining part is strongly enhanced. It is possible to submit new, not yet classified enzymes to BRENDA, which then are reviewed and classified by the International Union of Biochemistry and Molecular Biology. A new SBML output format of BRENDA kinetic data allows the construction of organism-specific metabolic models

Genome-scale architecture of small molecule regulatory networks and the fundamental trade-off between regulation and enzymatic activity

Metabolic flux is in part regulated by endogenous small molecules that modulate the catalytic activity of an enzyme, e.g., allosteric inhibition. In contrast to transcriptional regulation of enzymes, technical limitations have hindered the production of a genome-scale atlas of small molecule-enzyme regulatory interactions. Here, we develop a framework leveraging the vast, but fragmented, biochemical literature to reconstruct and analyze the small molecule regulatory network (SMRN) of the model organism Escherichia coli, including the primary metabolite regulators and enzyme targets. Using metabolic control analysis, we prove a fundamental trade-off between regulation and enzymatic activity, and we combine it with metabolomic measurements and the SMRN to make inferences on the sensitivity of enzymes to their regulators. Generalizing the analysis to other organisms, we identify highly conserved regulatory interactions across evolutionarily divergent species, further emphasizing a critical role for small molecule interactions in the maintenance of metabolic homeostasis.P30 CA008748 - NCI NIH HHS; R01 GM121950 - NIGMS NIH HH

Prediction of enzyme kinetic parameters based on statistical learning

Values of enzyme kinetic parameters are a key requisite for the kinetic modelling of biochemical systems. For most kinetic parameters, however, not even an order of magnitude is known, so the estimation of model parameters from experimental data remains a major task in systems biology. We propose a statistical approach to infer values for kinetic parameters across species and enzymes making use of parameter values that have been measured under various conditions and that are nowadays stored in databases. We fit the data by a statistical regression model in which the substrate, the combination enzyme-substrate and the combination organism-substrate have a linear effect on the logarithmic parameter value. As a result, we obtain predictions and error ranges for unknown enzyme parameters. We apply our method to decadic logarithmic Michaelis-Menten constants from the BRENDA database and confirm the results with leave-one-out crossvalidation, in which we mask one value at a time and predict it from the remaining data. For a set of 8 metabolites we obtain a standard prediction error of 1.01 for the deviation of the predicted values from the true values, while the standard deviation of the experimental values is 1.16. The method is applicable to other types of kinetic parameters for which many experimental data are available

Serine biosynthesis with one carbon catabolism represents a novel pathway for ATP generation in cells using alternative glycolysis with zero net ATP production

Recent experimental evidence indicates that some cancer cells have an alternative glycolysis pathway with net zero ATP production, implying that upregulation of glycolysis in these cells may not be related to the generation of ATP. Here we use a genome-scale model of human cell metabolism to investigate the potential metabolic alterations in cells using net zero ATP glycolysis. We uncover a novel pathway for ATP generation that involves reactions from the serine biosynthesis and one-carbon metabolism pathways. This pathway has a predicted two-fold higher flux rate in cells using net zero ATP glycolysis than those using standard glycolysis and generates twice as much ATP with significantly lower rate of lactate- but higher rate of alanine secretion. Thus, in cells using the standard- or the net zero ATP glycolysis pathways a significant portion of the glycolysis flux is always associated with ATP generation, and the ratio between the flux rates of the two pathways determines the rate of ATP generation and lactate and alanine secretion during glycolysis

Serine biosynthesis with one carbon catabolism represents a novel pathway for ATP generation in cells using alternative glycolysis with zero net ATP production

Recent experimental evidence indicates that some cancer cells have an alternative glycolysis pathway with net zero ATP production, implying that upregulation of glycolysis in these cells may not be related to the generation of ATP. Here we use a genome-scale model of human cell metabolism to investigate the potential metabolic alterations in cells using net zero ATP glycolysis. We uncover a novel pathway for ATP generation that involves reactions from the serine biosynthesis and one-carbon metabolism pathways. This pathway has a predicted two-fold higher flux rate in cells using net zero ATP glycolysis than those using standard glycolysis and generates twice as much ATP with significantly lower rate of lactate- but higher rate of alanine secretion. Thus, in cells using the standard- or the net zero ATP glycolysis pathways a significant portion of the glycolysis flux is always associated with ATP generation, and the ratio between the flux rates of the two pathways determines the rate of ATP generation and lactate and alanine secretion during glycolysis

A synthetic biochemistry platform for cell free production of monoterpenes from glucose.

Cell-free systems designed to perform complex chemical conversions of biomass to biofuels or commodity chemicals are emerging as promising alternatives to the metabolic engineering of living cells. Here we design a system comprises 27 enzymes for the conversion of glucose into monoterpenes that generates both NAD(P)H and ATP in a modified glucose breakdown module and utilizes both cofactors for building terpenes. Different monoterpenes are produced in our system by changing the terpene synthase enzyme. The system is stable for the production of limonene, pinene and sabinene, and can operate continuously for at least 5 days from a single addition of glucose. We obtain conversion yields >95% and titres >15 g l-1. The titres are an order of magnitude over cellular toxicity limits and thus difficult to achieve using cell-based systems. Overall, these results highlight the potential of synthetic biochemistry approaches for producing bio-based chemicals

Machine learning applied to enzyme turnover numbers reveals protein structural correlates and improves metabolic models.

Knowing the catalytic turnover numbers of enzymes is essential for understanding the growth rate, proteome composition, and physiology of organisms, but experimental data on enzyme turnover numbers is sparse and noisy. Here, we demonstrate that machine learning can successfully predict catalytic turnover numbers in Escherichia coli based on integrated data on enzyme biochemistry, protein structure, and network context. We identify a diverse set of features that are consistently predictive for both in vivo and in vitro enzyme turnover rates, revealing novel protein structural correlates of catalytic turnover. We use our predictions to parameterize two mechanistic genome-scale modelling frameworks for proteome-limited metabolism, leading to significantly higher accuracy in the prediction of quantitative proteome data than previous approaches. The presented machine learning models thus provide a valuable tool for understanding metabolism and the proteome at the genome scale, and elucidate structural, biochemical, and network properties that underlie enzyme kinetics

Metabolite concentrations, fluxes and free energies imply efficient enzyme usage.

In metabolism, available free energy is limited and must be divided across pathway steps to maintain a negative ΔG throughout. For each reaction, ΔG is log proportional both to a concentration ratio (reaction quotient to equilibrium constant) and to a flux ratio (backward to forward flux). Here we use isotope labeling to measure absolute metabolite concentrations and fluxes in Escherichia coli, yeast and a mammalian cell line. We then integrate this information to obtain a unified set of concentrations and ΔG for each organism. In glycolysis, we find that free energy is partitioned so as to mitigate unproductive backward fluxes associated with ΔG near zero. Across metabolism, we observe that absolute metabolite concentrations and ΔG are substantially conserved and that most substrate (but not inhibitor) concentrations exceed the associated enzyme binding site dissociation constant (Km or Ki). The observed conservation of metabolite concentrations is consistent with an evolutionary drive to utilize enzymes efficiently given thermodynamic and osmotic constraints

Automation on the generation of genome scale metabolic models

Background: Nowadays, the reconstruction of genome scale metabolic models is a non-automatized and interactive process based on decision taking. This lengthy process usually requires a full year of one person's work in order to satisfactory collect, analyze and validate the list of all metabolic reactions present in a specific organism. In order to write this list, one manually has to go through a huge amount of genomic, metabolomic and physiological information. Currently, there is no optimal algorithm that allows one to automatically go through all this information and generate the models taking into account probabilistic criteria of unicity and completeness that a biologist would consider. Results: This work presents the automation of a methodology for the reconstruction of genome scale metabolic models for any organism. The methodology that follows is the automatized version of the steps implemented manually for the reconstruction of the genome scale metabolic model of a photosynthetic organism, {\it Synechocystis sp. PCC6803}. The steps for the reconstruction are implemented in a computational platform (COPABI) that generates the models from the probabilistic algorithms that have been developed. Conclusions: For validation of the developed algorithm robustness, the metabolic models of several organisms generated by the platform have been studied together with published models that have been manually curated. Network properties of the models like connectivity and average shortest mean path of the different models have been compared and analyzed.Comment: 24 pages, 2 figures, 2 table

metaSHARK: software for automated metabolic network prediction from DNA sequence and its application to the genomes of Plasmodium falciparum and Eimeria tenella

The metabolic SearcH And Reconstruction Kit (metaSHARK) is a new fully automated software package for the detection of enzyme-encoding genes within unannotated genome data and their visualization in the context of the surrounding metabolic network. The gene detection package (SHARKhunt) runs on a Linux systemand requires only a set of raw DNA sequences (genomic, expressed sequence tag and/ or genome survey sequence) as input. Its output may be uploaded to our web-based visualization tool (SHARKview) for exploring and comparing data from different organisms. We first demonstrate the utility of the software by comparing its results for the raw Plasmodium falciparum genome with the manual annotations available at the PlasmoDB and PlasmoCyc websites. We then apply SHARKhunt to the unannotated genome sequences of the coccidian parasite Eimeria tenella and observe that, at an E-value cut-off of 10(-20), our software makes 142 additional assertions of enzymatic function compared with a recent annotation package working with translated open reading frame sequences. The ability of the software to cope with low levels of sequence coverage is investigated by analyzing assemblies of the E.tenella genome at estimated coverages from 0.5x to 7.5x. Lastly, as an example of how metaSHARK can be used to evaluate the genomic evidence for specific metabolic pathways, we present a study of coenzyme A biosynthesis in P.falciparum and E.tenella