From Bugs to Bioplastics: Total (+)-Dihydrocarvide Biosynthesis by Engineered Escherichia coli.

The monoterpenoid lactone derivative (+)-dihydrocarvide ((+)-DHCD) can be polymerised to form shape-memory polymers. Synthetic biology routes from simple, inexpensive carbon sources are an attractive, alternative route over chemical synthesis from (R)-carvone. We have demonstrated a proof-of-principle in vivo approach for the complete biosynthesis of (+)-DHCD from glucose in Escherichia coli (6.6 mg L-1 ). The pathway is based on the Mentha spicata route to (R)-carvone, with the addition of an 'ene'-reductase and Baeyer-Villiger cyclohexanone monooxygenase. Co-expression with a limonene synthesis pathway enzyme enables complete biocatalytic production within one microbial chassis. (+)-DHCD was successfully produced by screening multiple homologues of the pathway genes, combined with expression optimisation by selective promoter and/or ribosomal binding-site screening. This study demonstrates the potential application of synthetic biology approaches in the development of truly sustainable and renewable bioplastic monomers.

Biocatalytic Routes to Lactone Monomers for Polymer Production.

Monoterpenoids offer potential as biocatalytically derived monomer feedstocks for high-performance renewable polymers. We describe a biocatalytic route to lactone monomers menthide and dihydrocarvide employing Baeyer-Villiger monooxygenases (BVMOs) from Pseudomonas sp. HI-70 (CPDMO) and Rhodococcus sp. Phi1 (CHMOPhi1) as an alternative to organic synthesis. The regioselectivity of dihydrocarvide isomer formation was controlled by site-directed mutagenesis of three key active site residues in CHMOPhi1. A combination of crystal structure determination, molecular dynamics simulations, and mechanistic modeling using density functional theory on a range of models provides insight into the origins of the discrimination of the wild type and a variant CHMOPhi1 for producing different regioisomers of the lactone product. Ring-opening polymerizations of the resultant lactones using mild metal-organic catalysts demonstrate their utility in polymer production. This semisynthetic approach utilizing a biocatalytic step, non-petroleum feedstocks, and mild polymerization catalysts allows access to known and also to previously unreported and potentially novel lactone monomers and polymers.

Regulation of the activity of lactate dehydrogenases from four lactic acid bacteria.

Despite high similarity in sequence and catalytic properties, the l-lactate dehydrogenases (LDHs) in lactic acid bacteria (LAB) display differences in their regulation that may arise from their adaptation to different habitats. We combined experimental and computational approaches to investigate the effects of fructose 1,6-bisphosphate (FBP), phosphate (Pi), and ionic strength (NaCl concentration) on six LDHs from four LABs studied at pH 6 and pH 7. We found that 1) the extent of activation by FBP (Kact) differs. Lactobacillus plantarum LDH is not regulated by FBP, but the other LDHs are activated with increasing sensitivity in the following order: Enterococcus faecalis LDH2 ≤ Lactococcus lactis LDH2 < E. faecalis LDH1 < L. lactis LDH1 ≤ Streptococcus pyogenes LDH. This trend reflects the electrostatic properties in the allosteric binding site of the LDH enzymes. 2) For L. plantarum, S. pyogenes, and E. faecalis, the effects of Pi are distinguishable from the effect of changing ionic strength by adding NaCl. 3) Addition of Pi inhibits E. faecalis LDH2, whereas in the absence of FBP, Pi is an activator of S. pyogenes LDH, E. faecalis LDH1, and L. lactis LDH1 and LDH2 at pH 6. These effects can be interpreted by considering the computed binding affinities of Pi to the catalytic and allosteric binding sites of the enzymes modeled in protonation states corresponding to pH 6 and pH 7. Overall, the results show a subtle interplay among the effects of Pi, FBP, and pH that results in different regulatory effects on the LDHs of different LABs.

A Biological Route to Conjugated Alkenes: Microbial Production of Hepta-1,3,5-triene.

Conjugated alkenes such as dienes and polyenes have a range of applications as pharmaceutical agents and valuable building blocks in the polymer industry. Development of a renewable route to these compounds provides an alternative to fossil fuel derived production. The enzyme family of the UbiD decarboxylases offers substantial scope for alkene production, readily converting poly unsaturated acids. However, biochemical pathways producing the required substrates are poorly characterized, and UbiD-application has hitherto been limited to biological styrene production. Herein, we present a proof-of-principle study for microbial production of polyenes using a bioinspired strategy employing a polyketide synthase (PKS) in combination with a UbiD-enzyme. Deconstructing a bacterial iterative type II PKS enabled repurposing the broad-spectrum antibiotic andrimid biosynthesis pathway to access the metabolic intermediate 2,4,6-octatrienoic acid, a valuable chemical for material and pharmaceutical industry. Combination with the fungal ferulic acid decarboxylase (Fdc1) led to a biocatalytic cascade-type reaction for the production of hepta-1,3,5-triene in vivo. Our approach provides a novel route to generate unsaturated hydrocarbons and related chemicals and provides a blue-print for future development and application.

Systematic integration of experimental data and models in systems biology.

BACKGROUND: The behaviour of biological systems can be deduced from their mathematical models. However, multiple sources of data in diverse forms are required in the construction of a model in order to define its components and their biochemical reactions, and corresponding parameters. Automating the assembly and use of systems biology models is dependent upon data integration processes involving the interoperation of data and analytical resources. RESULTS: Taverna workflows have been developed for the automated assembly of quantitative parameterised metabolic networks in the Systems Biology Markup Language (SBML). A SBML model is built in a systematic fashion by the workflows which starts with the construction of a qualitative network using data from a MIRIAM-compliant genome-scale model of yeast metabolism. This is followed by parameterisation of the SBML model with experimental data from two repositories, the SABIO-RK enzyme kinetics database and a database of quantitative experimental results. The models are then calibrated and simulated in workflows that call out to COPASIWS, the web service interface to the COPASI software application for analysing biochemical networks. These systems biology workflows were evaluated for their ability to construct a parameterised model of yeast glycolysis. CONCLUSIONS: Distributed information about metabolic reactions that have been described to MIRIAM standards enables the automated assembly of quantitative systems biology models of metabolic networks based on user-defined criteria. Such data integration processes can be implemented as Taverna workflows to provide a rapid overview of the components and their relationships within a biochemical system.

KiPar, a tool for systematic information retrieval regarding parameters for kinetic modelling of yeast metabolic pathways.

MOTIVATION: Most experimental evidence on kinetic parameters is buried in the literature, whose manual searching is complex, time consuming and partial. These shortcomings become particularly acute in systems biology, where these parameters need to be integrated into detailed, genome-scale, metabolic models. These problems are addressed by KiPar, a dedicated information retrieval system designed to facilitate access to the literature relevant for kinetic modelling of a given metabolic pathway in yeast. Searching for kinetic data in the context of an individual pathway offers modularity as a way of tackling the complexity of developing a full metabolic model. It is also suitable for large-scale mining, since multiple reactions and their kinetic parameters can be specified in a single search request, rather than one reaction at a time, which is unsuitable given the size of genome-scale models. RESULTS: We developed an integrative approach, combining public data and software resources for the rapid development of large-scale text mining tools targeting complex biological information. The user supplies input in the form of identifiers used in relevant data resources to refer to the concepts of interest, e.g. EC numbers, GO and SBO identifiers. By doing so, the user is freed from providing any other knowledge or terminology concerned with these concepts and their relations, since they are retrieved from these and cross-referenced resources automatically. The terminology acquired is used to index the literature by mapping concepts to their synonyms, and then to textual documents mentioning them. The indexing results and the previously acquired knowledge about relations between concepts are used to formulate complex search queries aiming at documents relevant to the user's information needs. The conceptual approach is demonstrated in the implementation of KiPar. Evaluation reveals that KiPar performs better than a Boolean search. The precision achieved for abstracts (60%) and full-text articles (48%) is considerably better than the baseline precision (44% and 24%, respectively). The baseline recall is improved by 36% for abstracts and by 100% for full text. It appears that full-text articles are a much richer source of information on kinetic data than are their abstracts. Finally, the combined results for abstracts and full text compared with the curated literature provide high values for relative recall (88%) and novelty ratio (92%), suggesting that the system is able to retrieve a high proportion of new documents. AVAILABILITY: Source code and documentation are available at: (http://www.mcisb.org/resources/kipar/).

Magnetic field effects as a result of the radical pair mechanism are unlikely in redox enzymes.

Environmental exposure to electromagnetic fields is potentially carcinogenic. The radical pair mechanism is considered the most feasible mechanism of interaction between weak magnetic fields encountered in our environment and biochemical systems. Radicals are abundant in biology, both as free radicals and reaction intermediates in enzyme mechanisms. The catalytic cycles of some flavin-dependent enzymes are either known or potentially involve radical pairs. Here, we have investigated the magnetic field sensitivity of a number of flavoenzymes with important cellular roles. We also investigated the magnetic field sensitivity of a model system involving stepwise reduction of a flavin analogue by a nicotinamide analogue-a reaction known to proceed via a radical pair. Under the experimental conditions used, magnetic field sensitivity was not observed in the reaction kinetics from stopped-flow measurements in any of the systems studied. Although widely implicated in radical pair chemistry, we conclude that thermally driven, flavoenzyme-catalysed reactions are unlikely to be influenced by exposure to external magnetic fields.

The photochemical mechanism of a B12-dependent photoreceptor protein.

The coenzyme B12-dependent photoreceptor protein, CarH, is a bacterial transcriptional regulator that controls the biosynthesis of carotenoids in response to light. On binding of coenzyme B12 the monomeric apoprotein forms tetramers in the dark, which bind operator DNA thus blocking transcription. Under illumination the CarH tetramer dissociates, weakening its affinity for DNA and allowing transcription. The mechanism by which this occurs is unknown. Here we describe the photochemistry in CarH that ultimately triggers tetramer dissociation; it proceeds via a cob(III)alamin intermediate, which then forms a stable adduct with the protein. This pathway is without precedent and our data suggest it is independent of the radical chemistry common to both coenzyme B12 enzymology and its known photochemistry. It provides a mechanistic foundation for the emerging field of B12 photobiology and will serve to inform the development of a new class of optogenetic tool for the control of gene expression.

A new regulatory principle for in vivo biochemistry: pleiotropic low affinity regulation by the adenine nucleotides--illustrated for the glycolytic enzymes of Saccharomyces cerevisiae.

Enzymology tends to focus on highly specific effects of substrates, allosteric modifiers, and products occurring at low concentrations, because these are most informative about the enzyme's catalytic mechanism. We hypothesized that at relatively high in vivo concentrations, important molecular monitors of the state of living cells, such as ATP, affect multiple enzymes of the former and that these interactions have gone unnoticed in enzymology. We test this hypothesis in terms of the effect that ATP, ADP, and AMP might have on the major free-energy delivering pathway of the yeast Saccharomyces cerevisiae. Assaying cell-free extracts, we collected a comprehensive set of quantitative kinetic data concerning the enzymes of the glycolytic and the ethanol fermentation pathways. We determined systematically the extent to which the enzyme activities depend on the concentrations of the adenine nucleotides. We found that the effects of the adenine nucleotides on enzymes catalysing reactions in which they are not directly involved as substrate or product, are substantial. This includes effects on the Michaelis-Menten constants, adding new perspective on these, 100 years after their introduction.

A model of yeast glycolysis based on a consistent kinetic characterisation of all its enzymes.

We present an experimental and computational pipeline for the generation of kinetic models of metabolism, and demonstrate its application to glycolysis in Saccharomyces cerevisiae. Starting from an approximate mathematical model, we employ a "cycle of knowledge" strategy, identifying the steps with most control over flux. Kinetic parameters of the individual isoenzymes within these steps are measured experimentally under a standardised set of conditions. Experimental strategies are applied to establish a set of in vivo concentrations for isoenzymes and metabolites. The data are integrated into a mathematical model that is used to predict a new set of metabolite concentrations and reevaluate the control properties of the system. This bottom-up modelling study reveals that control over the metabolic network most directly involved in yeast glycolysis is more widely distributed than previously thought.

Towards a full quantitative description of yeast metabolism a systematic approach for estimating the kinetic parameters of isoenzymes under in vivo like conditions.

In order to produce a full quantitative description of yeast metabolism, a number of kinetic parameters of enzymes that are important for energy metabolism must be determined experimentally. We aim to determine the prospective in vivo kinetic properties of a range of yeast-purified isoenzymes that are important in energy metabolism, with respect to the concentration of their substrates and products. This endeavor forms part of our systems biology pipeline to facilitate the production of bottom-up models of metabolism. Within this workflow, we implement an infrastructure for medium- to high-throughput determination of the kinetic properties of purified isoenzymes in in vivo like conditions. This includes the use of the KineticsWizard software for data capture and analysis. The captured experimental data are analyzed by the software and subsequently stored in appropriate repositories (MeMo-RK and SABIO-RK). While we focus initially on glycolysis in Saccharomyces cerevisiae, our methodology is generic and can be widely applied to the study of other enzymes and pathways in yeast and other organisms.

Building a kinetic model of trehalose biosynthesis in Saccharomyces cerevisiae.

In this chapter, we describe the steps needed to create a kinetic model of a metabolic pathway based on kinetic data from experimental measurements and literature review. Our methodology is presented by utilizing the example of trehalose metabolism in yeast. The biology of the trehalose cycle is briefly reviewed and discussed.

Enzyme kinetics informatics: from instrument to browser.

A limited number of publicly available resources provide access to enzyme kinetic parameters. These have been compiled through manual data mining of published papers, not from the original, raw experimental data from which the parameters were calculated. This is largely due to the lack of software or standards to support the capture, analysis, storage and dissemination of such experimental data. Introduced here is an integrative system to manage experimental enzyme kinetics data from instrument to browser. The approach is based on two interrelated databases: the existing SABIO-RK database, containing kinetic data and corresponding metadata, and the newly introduced experimental raw data repository, MeMo-RK. Both systems are publicly available by web browser and web service interfaces and are configurable to ensure privacy of unpublished data. Users of this system are provided with the ability to view both kinetic parameters and the experimental raw data from which they are calculated, providing increased confidence in the data. A data analysis and submission tool, the kineticswizard, has been developed to allow the experimentalist to perform data collection, analysis and submission to both data resources. The system is designed to be extensible, allowing integration with other manufacturer instruments covering a range of analytical techniques.

Role of active site residues and solvent in proton transfer and the modulation of flavin reduction potential in bacterial morphinone reductase.

The reactions of several active site mutant forms of bacterial morphinone reductase (MR) with NADH and 2-cyclohexen-1-one as substrates have been studied by stopped-flow and steady-state kinetic methods and redox potentiometry. The enzymes were designed to (i) probe a role for potential proton donors (Tyr-72 and Tyr-356) in the oxidative half-reaction of MR; (ii) assess the function of a highly conserved tryptophan residue (Trp-106) in catalysis; (iii) investigate the role of Thr-32 in modulating the FMN reduction potential and catalysis. The Y72F and Y356F enzymes retained activity in both steady-state and stopped-flow kinetic studies, indicating they do not serve as key proton donors in the oxidative reaction of MR. Taken together with our recently published data (Messiha, H. L., Munro, A. W., Bruce, N. C., Barsukov, I., and Scrutton, N. S. (2005) J. Biol. Chem. 280, 4627-4631) that rule out roles for Cys-191 (corresponding with the proton donor, Tyr-196, in the structurally related OYE1 enzyme) and His-186 as proton donors, we infer solvent is the source of the proton in the oxidative half-reaction of MR. We demonstrate a key role for Thr-32 in modulating the reduction potential of the FMN, which is decreased approximately 50 mV in the T32A mutant MR. This effects a change in rate-limiting step in the catalytic cycle of the T32A enzyme with the oxidizing substrate 2-cyclohexenone. Despite the conservation of Trp-106 throughout the OYE family, we show this residue does not play a major role in catalysis, although affects on substrate and coenzyme binding are observed in a W106F enzyme. Our studies show some similarities, but also major differences, in the catalytic mechanism of MR and OYE1, and emphasize the need for caution in inferring mechanism by structural comparison of highly related enzymes in the absence of solution studies.