Deoxyxylulose 5-Phosphate Synthase Does Not Play a Major Role in Regulating the Methylerythritol 4-Phosphate Pathway in Poplar.

The plastidic 2-C-methylerythritol 4-phosphate (MEP) pathway supplies the precursors of a large variety of essential plant isoprenoids, but its regulation is still not well understood. Using metabolic control analysis (MCA), we examined the first enzyme of this pathway, 1-deoxyxylulose 5-phosphate synthase (DXS), in multiple grey poplar (Populus × canescens) lines modified in their DXS activity. Single leaves were dynamically labeled with 13CO2 in an illuminated, climate-controlled gas exchange cuvette coupled to a proton transfer reaction mass spectrometer, and the carbon flux through the MEP pathway was calculated. Carbon was rapidly assimilated into MEP pathway intermediates and labeled both the isoprene released and the IDP+DMADP pool by up to 90%. DXS activity was increased by 25% in lines overexpressing the DXS gene and reduced by 50% in RNA interference lines, while the carbon flux in the MEP pathway was 25-35% greater in overexpressing lines and unchanged in RNA interference lines. Isoprene emission was also not altered in these different genetic backgrounds. By correlating absolute flux to DXS activity under different conditions of light and temperature, the flux control coefficient was found to be low. Among isoprenoid end products, isoprene itself was unchanged in DXS transgenic lines, but the levels of the chlorophylls and most carotenoids measured were 20-30% less in RNA interference lines than in overexpression lines. Our data thus demonstrate that DXS in the isoprene-emitting grey poplar plays only a minor part in controlling flux through the MEP pathway.

Quantifying redox transcription factor dynamics as a tool to investigate redox signalling.

A critical feature of the cellular antioxidant response is the induction of gene expression by redox-sensitive transcription factors. In many cells, activating these transcription factors is a dynamic process involving multiple redox steps, but it is unclear how these dynamics should be measured. Here, we show how the dynamic profile of the Schizosaccharomyces pombe Pap1 transcription factor is quantifiable by three parameters: signal amplitude, signal time and signal duration. In response to increasing hydrogen peroxide concentrations, the Pap1 amplitude decreased while the signal time and duration showed saturable increases. In co-response plots, these parameters showed a complex, non-linear relationship to the mRNA levels of four Pap1-regulated genes. We also demonstrate that hydrogen peroxide and tert-butyl hydroperoxide trigger quantifiably distinct Pap1 activation profiles and transcriptional responses. Based on these findings, we propose that different oxidants and oxidant concentrations modulate the Pap1 dynamic profile, leading to specific transcriptional responses. We further show how the effect of combination and pre-exposure stresses on Pap1 activation dynamics can be quantified using this approach. This method is therefore a valuable addition to the redox signalling toolbox that may illuminate the role of dynamics in determining appropriate responses to oxidative stress.

A special issue of Essays in Biochemistry on computational biology.

Computational biology is a diverse research field that has gained increasing importance over the last two decades. Broadly, it aims to apply computational approaches to advance our understanding of biological systems. This can take place on multiple levels, for example, by creating computational models of specific biological systems, by developing algorithms that assist in the analysis of experimental data, or by investigating fundamental biological design principles through modelling. The articles in this special issue highlight and review four such distinct applications of computational biology.

Computational models as catalysts for investigating redoxin systems.

Thioredoxin, glutaredoxin and peroxiredoxin systems play central roles in redox regulation, signaling and metabolism in cells. In these systems, reducing equivalents from NAD(P)H are transferred by coupled thiol-disulfide exchange reactions to redoxins which then reduce a wide array of targets. However, the characterization of redoxin activity has been unclear, with redoxins regarded as enzymes in some studies and redox metabolites in others. Consequently, redoxin activities have been quantified by enzyme kinetic parameters in vitro, and redox potentials or redox ratios within cells. By analyzing all the reactions within these systems, computational models showed that many kinetic properties attributed to redoxins were due to system-level effects. Models of cellular redoxin networks have also been used to estimate intracellular hydrogen peroxide levels, analyze redox signaling and couple omic and kinetic data to understand the regulation of these networks in disease. Computational modeling has emerged as a powerful complementary tool to traditional redoxin enzyme kinetic and cellular assays that integrates data from a number of sources into a single quantitative framework to accelerate the analysis of redoxin systems.

Modelling the Decamerisation Cycle of PRDX1 and the Inhibition-like Effect on Its Peroxidase Activity.

Peroxiredoxins play central roles in the detoxification of reactive oxygen species and have been modelled across multiple organisms using a variety of kinetic methods. However, the peroxiredoxin dimer-to-decamer transition has been underappreciated in these studies despite the 100-fold difference in activity between these forms. This is due to the lack of available kinetics and a theoretical framework for modelling this process. Using published isothermal titration calorimetry data, we obtained association and dissociation rate constants of 0.050 µM-4·s-1 and 0.055 s-1, respectively, for the dimer-decamer transition of human PRDX1. We developed an approach that greatly reduces the number of reactions and species needed to model the peroxiredoxin decamer oxidation cycle. Using these data, we simulated horse radish peroxidase competition and NADPH-oxidation linked assays and found that the dimer-decamer transition had an inhibition-like effect on peroxidase activity. Further, we incorporated this dimer-decamer topology and kinetics into a published and validated in vivo model of PRDX2 in the erythrocyte and found that it almost perfectly reconciled experimental and simulated responses of PRDX2 oxidation state to hydrogen peroxide insult. By accounting for the dimer-decamer transition of peroxiredoxins, we were able to resolve several discrepancies between experimental data and available kinetic models.

Atypical network topologies enhance the reductive capacity of pathogen thiol antioxidant defense networks.

Infectious diseases are a significant health burden for developing countries, particularly with the rise of multidrug resistance. There is an urgent need to elucidate the factors underlying the persistence of pathogens such as Mycobacterium tuberculosis, Plasmodium falciparum and Trypanosoma brucei. In contrast to host cells, these pathogens traverse multiple and varied redox environments during their infectious cycles, including exposure to high levels of host-derived reactive oxygen species. Pathogen antioxidant defenses such as the peroxiredoxin and thioredoxin systems play critical roles in the redox stress tolerance of these cells. However, many of the kinetic rate constants obtained for the pathogen peroxiredoxins are broadly similar to their mammalian homologs and therefore, their contributions to the redox tolerances within these cells are enigmatic. Using graph theoretical analysis, we show that compared to a canonical Escherichia coli redoxin network, pathogen redoxin networks contain unique network connections (motifs) between their thioredoxins and peroxiredoxins. Analysis of these motifs reveals that they increase the hydroperoxide reduction capacity of these networks and, in response to an oxidative insult, can distribute fluxes into specific thioredoxin-dependent pathways. Our results emphasize that the high oxidative stress tolerance of these pathogens depends on both the kinetic parameters for hydroperoxide reduction and the connectivity within their thioredoxin/peroxiredoxin systems.

HDR, the last enzyme in the MEP pathway, differently regulates isoprenoid biosynthesis in two woody plants.

Dimethylallyl diphosphate (DMADP) and isopentenyl diphosphate (IDP) serves as the universal C5 precursors of isoprenoid biosynthesis in plants. These compounds are formed by the last step of the 2-C-methyl-D-erythritol 4-phosphate (MEP) pathway, catalyzed by (E)-4-hydroxy-3-methylbut-2-en-1-yl diphosphate reductase (HDR). In this study, we investigated the major HDR isoforms of two woody plant species, Norway spruce (Picea abies) and gray poplar (Populus × canescens), to determine how they regulate isoprenoid formation. Since each of these species has a distinct profile of isoprenoid compounds, they may require different proportions of DMADP and IDP with proportionally more IDP being needed to make larger isoprenoids. Norway spruce contained two major HDR isoforms differing in their occurrence and biochemical characteristics. PaHDR1 produced relatively more IDP than PaHDR2 and it encoding gene was expressed constitutively in leaves, likely serving to form substrate for production of carotenoids, chlorophylls, and other primary isoprenoids derived from a C20 precursor. On the other hand, Norway spruce PaHDR2 produced relatively more DMADP than PaHDR1 and its encoding gene was expressed in leaves, stems, and roots, both constitutively and after induction with the defense hormone methyl jasmonate. This second HDR enzyme likely forms a substrate for the specialized monoterpene (C10), sesquiterpene (C15), and diterpene (C20) metabolites of spruce oleoresin. Gray poplar contained only one dominant isoform (named PcHDR2) that produced relatively more DMADP and the gene of which was expressed in all organs. In leaves, where the requirement for IDP is high to make the major carotenoid and chlorophyll isoprenoids derived from C20 precursors, excess DMADP may accumulate, which could explain the high rate of isoprene (C5) emission. Our results provide new insights into the biosynthesis of isoprenoids in woody plants under conditions of differentially regulated biosynthesis of the precursors IDP and DMADP.

EnzymeML: seamless data flow and modeling of enzymatic data.

The design of biocatalytic reaction systems is highly complex owing to the dependency of the estimated kinetic parameters on the enzyme, the reaction conditions, and the modeling method. Consequently, reproducibility of enzymatic experiments and reusability of enzymatic data are challenging. We developed the XML-based markup language EnzymeML to enable storage and exchange of enzymatic data such as reaction conditions, the time course of the substrate and the product, kinetic parameters and the kinetic model, thus making enzymatic data findable, accessible, interoperable and reusable (FAIR). The feasibility and usefulness of the EnzymeML toolbox is demonstrated in six scenarios, for which data and metadata of different enzymatic reactions are collected and analyzed. EnzymeML serves as a seamless communication channel between experimental platforms, electronic lab notebooks, tools for modeling of enzyme kinetics, publication platforms and enzymatic reaction databases. EnzymeML is open and transparent, and invites the community to contribute. All documents and codes are freely available at https://enzymeml.org .

Manganese Privation-Induced Transcriptional Upregulation of the Class IIa Bacteriocin Plantaricin 423 in Lactobacillus plantarum Strain 423.

Plantaricin 423 is produced by Lactobacillus plantarum 423 using the pla biosynthetic operon located on the 8,188-bp plasmid pPLA4. As with many class IIa bacteriocin operons, the pla operon carries biosynthetic genes (plaA, precursor peptide; plaB, immunity; plaC, accessory; and plaD, ABC transporter) but does not carry local regulatory genes. Little is known about the regulatory mechanisms involved in the expression of the apparently regulationless class IIa bacteriocins, such as plantaricin 423. In this study, phylogenetic analysis of class IIa immunity proteins indicated that at least three distinct clades exist, which were then used to subgroup the class IIa operons. It became evident that the absence of classical quorum-sensing genes on mobile bacteriocin-encoding elements is a predisposition of the subgroup that includes plantaricin 423, pediocin AcH/PA-1, divercin V41, enterocin A, leucocin-A and -B, mesentericin Y105, and sakacin G. Further analysis of the subgroup suggested that the regulation of these class IIa operons is linked to transition metal homeostasis in the host. By using a fluorescent promoter-reporter system in Lactobacillus plantarum 423, transcriptional regulation of plantaricin 423 was shown to be upregulated in response to manganese privation. IMPORTANCE Lactic acid bacteria hold huge industrial application and economic value, especially bacteriocinogenic strains, which further aids in the exclusion of specific foodborne pathogens. Since bacteriocinogenic strains are sought after, it is equally important to understand the mechanism of bacteriocin regulation. This is currently an understudied aspect of class IIa operons. Our research suggests the existence of a previously undescribed mode of class IIa bacteriocin regulation, whereby bacteriocin expression is linked to management of the producer's transition metal homeostasis. This delocalized metalloregulatory model may fundamentally affect the selection of culture conditions for bacteriocin expression and change our understanding of class IIa bacteriocin gene transfer dynamics in a given microbiome.

Functional Characterisation of Three Glycine N-Acyltransferase Variants and the Effect on Glycine Conjugation to Benzoyl-CoA.

The glycine conjugation pathway in humans is involved in the metabolism of natural substrates and the detoxification of xenobiotics. The interactions between the various substrates in this pathway and their competition for the pathway enzymes are currently unknown. The pathway consists of a mitochondrial xenobiotic/medium-chain fatty acid: coenzyme A (CoA) ligase (ACSM2B) and glycine N-acyltransferase (GLYAT). The catalytic mechanism and substrate specificity of both of these enzymes have not been thoroughly characterised. In this study, the level of evolutionary conservation of GLYAT missense variants and haplotypes were analysed. From these data, haplotype variants were selected (156Asn > Ser, [17Ser > Thr,156Asn > Ser] and [156Asn > Ser,199Arg > Cys]) in order to characterise the kinetic mechanism of the enzyme over a wide range of substrate concentrations. The 156Asn > Ser haplotype has the highest frequency and the highest relative enzyme activity in all populations studied, and hence was used as the reference in this study. Cooperative substrate binding was observed, and the kinetic data were fitted to a two-substrate Hill equation. The coding region of the GLYAT gene was found to be highly conserved and the rare 156Asn > Ser,199Arg > Cys variant negatively affected the relative enzyme activity. Even though the 156Asn > Ser,199Arg > Cys variant had a higher affinity for benzoyl-CoA (s0.5,benz = 61.2 µM), kcat was reduced to 9.8% of the most abundant haplotype 156Asn > Ser (s0.5,benz = 96.6 µM), while the activity of 17Ser > Thr,156Asn > Ser (s0.5,benz = 118 µM) was 73% of 156Asn > Ser. The in vitro kinetic analyses of the effect of the 156Asn > Ser,199Arg > Cys variant on human GLYAT enzyme activity indicated that individuals with this haplotype might have a decreased ability to metabolise benzoate when compared to individuals with the 156Asn > Ser variant. Furthermore, the accumulation of acyl-CoA intermediates can inhibit ACSM2B leading to a reduction in mitochondrial energy production.

Effect of Drought on the Methylerythritol 4-Phosphate (MEP) Pathway in the Isoprene Emitting Conifer Picea glauca.

The methylerythritol 4-phosphate (MEP) pathway of isoprenoid biosynthesis produces chlorophyll side chains and compounds that function in resistance to abiotic stresses, including carotenoids, and isoprene. Thus we investigated the effects of moderate and severe drought on MEP pathway function in the conifer Picea glauca, a boreal species at risk under global warming trends. Although moderate drought treatment reduced the photosynthetic rate by over 70%, metabolic flux through the MEP pathway was reduced by only 37%. The activity of the putative rate-limiting step, 1-deoxy-D-xylulose-5-phosphate synthase (DXS), was also reduced by about 50%, supporting the key role of this enzyme in regulating pathway metabolic flux. However, under severe drought, as flux declined below detectable levels, DXS activity showed no significant decrease, indicating a much-reduced role in controlling flux under these conditions. Both MEP pathway intermediates and the MEP pathway product isoprene incorporate administered 13CO2 to high levels (75-85%) under well-watered control conditions indicating a close connection to photosynthesis. However, this incorporation declined precipitously under drought, demonstrating exploitation of alternative carbon sources. Despite the reductions in MEP pathway flux and intermediate pools, there was no detectable decline in most major MEP pathway products under drought (except for violaxanthin under moderate and severe stress and isoprene under severe stress) suggesting that the pathway is somehow buffered against this stress. The resilience of the MEP pathway under drought may be a consequence of the importance of the metabolites formed under these conditions.

The thioredoxin redox potential and redox charge are surrogate measures for flux in the thioredoxin system.

The thioredoxin system plays a central role in intracellular redox regulation and its dysregulation is associated with a number of pathologies. However, the connectivity within this system poses a significant challenge for quantification and consequently several disparate measures have been used to characterize the system. For in vitro studies, the thioredoxin system flux has been measured by NADPH oxidation while the thioredoxin redox state has been used to estimate the activity of the system in vivo. The connection between these measures has been obscure although substrate saturation in the thioredoxin system results from the saturation of the thioredoxin redox cycle. We used computational modeling and in vitro kinetic assays to clarify the relationship between flux and the current in vivo measures of the thioredoxin system together with a novel measure, the thioredoxin redox charge (reduced thioredoxin/total thioredoxin). Our results revealed that the thioredoxin redox potential and redox charge closely tracked flux perturbations showing that these indices could be used as surrogate measures of the flux in vivo and, provide a mechanistic explanation for the previously observed correlations between thioredoxin oxidation and certain pathologies. While we found no significant difference in the linear correlations obtained for the thioredoxin redox potential and redox charge with the flux, the redox charge may be preferred because it is bounded between zero and one and can be determined over a wider range of conditions allowing for quantitative flux comparisons between cell types and conditions.

Investigation of the methylerythritol 4-phosphate pathway for microbial terpenoid production through metabolic control analysis.

BACKGROUND: Terpenoids are of high interest as chemical building blocks and pharmaceuticals. In microbes, terpenoids can be synthesized via the methylerythritol phosphate (MEP) or mevalonate (MVA) pathways. Although the MEP pathway has a higher theoretical yield, metabolic engineering has met with little success because the regulation of the pathway is poorly understood. RESULTS: We applied metabolic control analysis to the MEP pathway in Escherichia coli expressing a heterologous isoprene synthase gene (ispS). The expression of ispS led to the accumulation of isopentenyl pyrophosphate (IPP)/dimethylallyl pyrophosphate (DMAPP) and severely impaired bacterial growth, but the coexpression of ispS and isopentenyl diphosphate isomerase (idi) restored normal growth and wild-type IPP/DMAPP levels. Targeted proteomics and metabolomics analysis provided a quantitative description of the pathway, which was perturbed by randomizing the ribosome binding site in the gene encoding 1-deoxyxylulose 5-phosphate synthase (Dxs). Dxs has a flux control coefficient of 0.35 (i.e., a 1% increase in Dxs activity resulted in a 0.35% increase in pathway flux) in the isoprene-producing strain and therefore exerted significant control over the flux though the MEP pathway. At higher dxs expression levels, the intracellular concentration of 2-C-methyl-D-erythritol-2,4-cyclopyrophosphate (MEcPP) increased substantially in contrast to the other MEP pathway intermediates, which were linearly dependent on the abundance of Dxs. This indicates that 4-hydroxy-3-methylbut-2-en-1-yl diphosphate synthase (IspG), which consumes MEcPP, became saturated and therefore limited the flux towards isoprene. The higher intracellular concentrations of MEcPP led to the efflux of this intermediate into the growth medium. DISCUSSION: These findings show the importance of Dxs, Idi and IspG and metabolite export for metabolic engineering of the MEP pathway and will facilitate further approaches for the microbial production of valuable isoprenoids.

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.

An empirical analysis of enzyme function reporting for experimental reproducibility: Missing/incomplete information in published papers.

A key component of enzyme function experiments is reporting of considerable metadata, to allow other researchers to replicate, interpret properly or use fully the results. This paper evaluates the completeness of enzyme function data reporting for reproducibility. We present a detailed examination of 11 recent papers (and their supplementary material) from two leading journals. We found that in every paper we were not able to collect some critical information necessary to reproduce the enzyme function findings. Study of 100 papers used by the SABIO-RK database confirmed some of the most common omissions: concentration of enzyme or its substrates, identity of counter-ions in buffers. A computer system should be better at preventing such omissions, helping secure the scientific record. Many of the omissions found would be trapped by the currently available version of STRENDA DB.

STRENDA DB: enabling the validation and sharing of enzyme kinetics data.

Standards for reporting enzymology data (STRENDA) DB is a validation and storage system for enzyme function data that incorporates the STRENDA Guidelines. It provides authors who are preparing a manuscript with a user-friendly, web-based service that checks automatically enzymology data sets entered in the submission form that they are complete and valid before they are submitted as part of a publication to a journal.

PySCeSToolbox: a collection of metabolic pathway analysis tools.

Summary: PySCeSToolbox is an extension to the Python Simulator for Cellular Systems (PySCeS) that includes tools for performing generalized supply-demand analysis, symbolic metabolic control analysis, and a framework for investigating the kinetic and thermodynamic aspects of enzyme-catalyzed reactions. Each tool addresses a different aspect of metabolic behaviour, control, and regulation; the tools complement each other and can be used in conjunction to better understand higher level system behaviour. Availability and implementation: PySCeSToolbox is available on Linux, Mac OS X and Windows. It is licensed under the BSD 3-clause licence. Code, setup instructions and a link to documentation can be found at https://github.com/PySCeS/PyscesToolbox. Contact: jr@sun.ac.za. Supplementary information: Supplementary data are available at Bioinformatics online.

Quantitative measures for redox signaling.

Redox signaling is now recognized as an important regulatory mechanism for a number of cellular processes including the antioxidant response, phosphokinase signal transduction and redox metabolism. While there has been considerable progress in identifying the cellular machinery involved in redox signaling, quantitative measures of redox signals have been lacking, limiting efforts aimed at understanding and comparing redox signaling under normoxic and pathogenic conditions. Here we have outlined some of the accepted principles for redox signaling, including the description of hydrogen peroxide as a signaling molecule and the role of kinetics in conferring specificity to these signaling events. Based on these principles, we then develop a working definition for redox signaling and review a number of quantitative methods that have been employed to describe signaling in other systems. Using computational modeling and published data, we show how time- and concentration- dependent analyses, in particular, could be used to quantitatively describe redox signaling and therefore provide important insights into the functional organization of redox networks. Finally, we consider some of the key challenges with implementing these methods.

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.