Rapid reaction kinetics of proline dehydrogenase in the multifunctional proline utilization A protein

Metabolic collaboration between cells in the tumor microenvironment has a negligible effect on tumor growth

The tumor microenvironment is composed of a complex mixture of different cell types interacting under conditions of nutrient deprivation, but the metabolism therein is not fully understood due to difficulties in measuring metabolic fluxes and exchange of metabolites between different cell types in vivo. Genome-scale metabolic modeling enables estimation of such exchange fluxes as well as an opportunity to gain insight into the metabolic behavior of individual cell types. Here, we estimated the availability of nutrients and oxygen within the tumor microenvironment using concentration measurements from blood together with a metabolite diffusion model. In addition, we developed an approach to efficiently apply enzyme usage constraints in a comprehensive metabolic model of human cells. The combined modeling reproduced severe hypoxic conditions and the Warburg effect, and we found that limitations in enzymatic capacity contribute to cancer cells' preferential use of glutamine as a substrate to the citric acid cycle. Furthermore, we investigated the common hypothesis that some stromal cells are exploited by cancer cells to produce metabolites useful for the cancer cells. We identified over 200 potential metabolites that could support collaboration between cancer cells and cancer-associated fibroblasts, but when limiting to metabolites previously identified to participate in such collaboration, no growth advantage was observed. Our work highlights the importance of enzymatic capacity limitations for cell behaviors and exemplifies the utility of enzyme-constrained models for accurate prediction of metabolism in cells and tumor microenvironments.

Computational insights on the hydride and proton transfer mechanisms of L-proline dehydrogenase

L-Proline dehydrogenase (ProDH) is a flavin-dependent oxidoreductase, which catalyzes the oxidation of L-proline to (S)-1-pyrroline-5-carboxylate. Based on the experimental studies, a stepwise proton and hydride transfer mechanism is supported. According to this mechanism, the amino group of L-proline is deprotonated by a nearby Lys residue, which is followed by the hydride transfer process from C5 position of L-proline to N5 position of isoalloxazine ring of FAD. It was concluded that the hydride transfer step is rate limiting in the reductive half-reaction, however, in the overall reaction, the oxidation of FAD is the rate limiting step. In this study, we performed a computational mechanistic investigation based on ONIOM method to elucidate the mechanism of the reductive half-reaction corresponding to the oxidation of L-proline into iminoproline. Our calculations support the stepwise mechanism in which the deprotonation occurs initially as a fast step as result of a proton transfer from L-proline to the Lys residue. Subsequently, a hydride ion transfers from L-proline to FAD with a higher activation barrier. The enzyme-product complex showed a strong interaction between reduced FAD and iminoproline, which might help to explain why a step in the oxidative half-reaction is rate-limiting.

Proline Dehydrogenase and Pyrroline 5 Carboxylate Dehydrogenase from Mycobacterium tuberculosis: Evidence for Substrate Channeling

In Mycobacterium tuberculosis, proline dehydrogenase (PruB) and ∆1-pyrroline-5-carboxylate (P5C) dehydrogenase (PruA) are monofunctional enzymes that catalyze proline oxidation to glutamate via the intermediates P5C and L-glutamate-γ-semialdehyde. Both enzymes are essential for the replication of pathogenic M. tuberculosis. Highly active enzymes were expressed and purified using a Mycobacterium smegmatis expression system. The purified enzymes were characterized using natural substrates and chemically synthesized analogs. The structural requirements of the quinone electron acceptor were examined. PruB displayed activity with all tested lipoquinone analogs (naphthoquinone or benzoquinone). In PruB assays utilizing analogs of the native naphthoquinone [MK-9 (II-H2)] specificity constants Kcat/Km were an order of magnitude greater for the menaquinone analogs than the benzoquinone analogs. In addition, mycobacterial PruA was enzymatically characterized for the first time using exogenous chemically synthesized P5C. A Km value of 120 ± 0.015 µM was determined for P5C, while the Km value for NAD+ was shown to be 33 ± 4.3 µM. Furthermore, proline competitively inhibited PruA activity and coupled enzyme assays, suggesting that the recombinant purified monofunctional PruB and PruA enzymes of M. tuberculosis channel substrate likely increase metabolic flux and protect the bacterium from methylglyoxal toxicity.

CRISPR-assisted rational flux-tuning and arrayed CRISPRi screening of an L-proline exporter for L-proline hyperproduction

Development of hyperproducing strains is important for biomanufacturing of biochemicals and biofuels but requires extensive efforts to engineer cellular metabolism and discover functional components. Herein, we optimize and use the CRISPR-assisted editing and CRISPRi screening methods to convert a wild-type Corynebacterium glutamicum to a hyperproducer of L-proline, an amino acid with medicine, feed, and food applications. To facilitate L-proline production, feedback-deregulated variants of key biosynthetic enzyme γ-glutamyl kinase are screened using CRISPR-assisted single-stranded DNA recombineering. To increase the carbon flux towards L-proline biosynthesis, flux-control genes predicted by in silico analysis are fine-tuned using tailored promoter libraries. Finally, an arrayed CRISPRi library targeting all 397 transporters is constructed to discover an L-proline exporter Cgl2622. The final plasmid-, antibiotic-, and inducer-free strain produces L-proline at the level of 142.4 g/L, 2.90 g/L/h, and 0.31 g/g. The CRISPR-assisted strain development strategy can be used for engineering industrial-strength strains for efficient biomanufacturing.

Evidence for Proline Catabolic Enzymes in the Metabolism of Thiazolidine Carboxylates

Thiazolidine carboxylates such as thiazolidine-4-carboxylate (T4C) and thiazolidine-2-carboxylate (T2C) are naturally occurring sulfur analogues of proline. These compounds have been observed to have both beneficial and toxic effects in cells. Given that proline dehydrogenase has been proposed to be a key enzyme in the oxidative metabolism of thioprolines, we characterized T4C and T2C as substrates of proline catabolic enzymes using proline utilization A (PutA), which is a bifunctional enzyme with proline dehydrogenase (PRODH) and l-glutamate-γ-semialdehyde dehydrogenase (GSALDH) activities. PutA is shown here to catalyze the FAD-dependent PRODH oxidation of both T4C and T2C with catalytic efficiencies significantly higher than with proline. Stopped-flow experiments also demonstrate that l-T4C and l-T2C reduce PutA-bound FAD at rates faster than proline. Unlike proline, however, oxidation of T4C and T2C does not generate a substrate for NAD⁺-dependent GSALDH. Instead, PutA/PRODH oxidation of T4C leads to cysteine formation, whereas oxidation of T2C generates an apparently stable Δ⁴-thiazoline-2-carboxylate species. Our results provide new insights into the metabolism of T2C and T4C.

Probing the function of a ligand-modulated dynamic tunnel in bifunctional proline utilization A (PutA)

In many bacteria, the reactions of proline catabolism are catalyzed by the bifunctional enzyme known as proline utilization A (PutA). PutA catalyzes the two-step oxidation of l-proline to l-glutamate using distinct proline dehydrogenase (PRODH) and l-glutamate-γ-semialdehyde dehydrogenase (GSALDH) active sites, which are separated by over 40 Å and connected by a complex tunnel system. The tunnel system consists of a main tunnel that connects the two active sites and functions in substrate channeling, plus six ancillary tunnels whose functions are unknown. Here we used tunnel-blocking mutagenesis to probe the role of a dynamic ancillary tunnel (tunnel 2a) whose shape is modulated by ligand binding to the PRODH active site. The 1.90 Å resolution crystal structure of Geobacter sulfurreducens PutA variant A206W verified that the side chain of Trp206 cleanly blocks tunnel 2a without perturbing the surrounding structure. Steady-state kinetic measurements indicate the mutation impaired PRODH activity without affecting the GSALDH activity. Single-turnover experiments corroborated a severe impairment of PRODH activity with flavin reduction decreased by nearly 600-fold in A206W relative to wild-type. Substrate channeling is also significantly impacted as A206W exhibited a 3000-fold lower catalytic efficiency in coupled PRODH-GSALDH activity assays, which measure NADH formation as a function of proline. The structure suggests that Trp206 inhibits binding of the substrate l-proline by preventing the formation of a conserved glutamate-arginine ion pair and closure of the PRODH active site. Our data are consistent with tunnel 2a serving as an open space through which the glutamate of the ion pair travels during the opening and closing of the active site in response to binding l-proline. These results confirm the essentiality of the conserved ion pair in binding l-proline and support the hypothesis that the ion pair functions as a gate that controls access to the PRODH active site.

Production of l-glutamate family amino acids in Corynebacterium glutamicum: Physiological mechanism, genetic modulation, and prospects

l-glutamate family amino acids (GFAAs), consisting of l-glutamate, l-arginine, l-citrulline, l-ornithine, l-proline, l-hydroxyproline, γ-aminobutyric acid, and 5-aminolevulinic acid, are widely applied in the food, pharmaceutical, cosmetic, and animal feed industries, accounting for billions of dollars of market activity. These GFAAs have many functions, including being protein constituents, maintaining the urea cycle, and providing precursors for the biosynthesis of pharmaceuticals. Currently, the production of GFAAs mainly depends on microbial fermentation using Corynebacterium glutamicum (including its related subspecies Corynebacterium crenatum), which is substantially engineered through multistep metabolic engineering strategies. This review systematically summarizes recent advances in the metabolic pathways, regulatory mechanisms, and metabolic engineering strategies for GFAA accumulation in C. glutamicum and C. crenatum, which provides insights into the recent progress in l-glutamate-derived chemical production.

Detailed evaluation of pyruvate dehydrogenase complex inhibition in simulated exercise conditions

The mammalian pyruvate dehydrogenase complex (PDC) is a mitochondrial multienzyme complex that connects glycolysis to the tricarboxylic acid cycle by catalyzing pyruvate oxidation to produce acetyl-CoA, NADH, and CO₂. This reaction is required to aerobically utilize glucose, a preferred metabolic fuel, and is composed of three core enzymes: pyruvate dehydrogenase (E1), dihydrolipoyl transacetylase (E2), and dihydrolipoyl dehydrogenase (E3). The pyruvate-dehydrogenase-specific kinase (PDK) and pyruvate-dehydrogenase-specific phosphatase (PDP) are considered the main control mechanism of mammalian PDC activity. However, PDK and PDP activity are allosterically regulated by several effectors fully overlapping PDC substrates and products. This collection of positive and negative feedback mechanisms confounds simple predictions of relative PDC flux, especially when all effectors are dynamically modulated during metabolic states that exist in physiologically realistic conditions, such as exercise. Here, we provide, to our knowledge, the first globally fitted, pH-dependent kinetic model of the PDC accounting for the PDC core reaction because it is regulated by PDK, PDP, metal binding equilibria, and numerous allosteric effectors. The model was used to compute PDH regulatory complex flux as a function of previously determined metabolic conditions used to simulate exercise and demonstrates increased flux with exercise. Our model reveals that PDC flux in physiological conditions is primarily inhibited by product inhibition (∼60%), mostly NADH inhibition (∼30-50%), rather than phosphorylation cycle inhibition (∼40%), but the degree to which depends on the metabolic state and PDC tissue source.

Mechanisms of ligand binding

Many processes in chemistry and biology involve interactions of a ligand with its molecular target. Interest in the mechanism governing such interactions has dominated theoretical and experimental analysis for over a century. The interpretation of molecular recognition has evolved from a simple rigid body association of the ligand with its target to appreciation of the key role played by conformational transitions. Two conceptually distinct descriptions have had a profound impact on our understanding of mechanisms of ligand binding. The first description, referred to as induced fit, assumes that conformational changes follow the initial binding step to optimize the complex between the ligand and its target. The second description, referred to as conformational selection, assumes that the free target exists in multiple conformations in equilibrium and that the ligand selects the optimal one for binding. Both descriptions can be merged into more complex reaction schemes that better describe the functional repertoire of macromolecular systems. This review deals with basic mechanisms of ligand binding, with special emphasis on induced fit, conformational selection, and their mathematical foundations to provide rigorous context for the analysis and interpretation of experimental data. We show that conformational selection is a surprisingly versatile mechanism that includes induced fit as a mathematical special case and even captures kinetic properties of more complex reaction schemes. These features make conformational selection a dominant mechanism of molecular recognition in biology, consistent with the rich conformational landscape accessible to biological macromolecules being unraveled by structural biology.

De Novo Engineering of Corynebacterium glutamicum for l-Proline Production

l-Proline is an important amino acid that has various industrial applications. Industrial l-proline-producing strains are obtained by the mutagenesis of Corynebacterium glutamicum. In this study, the optimized C. glutamicum genome-editing tools were further applied in the de novo construction of a hyper-l-proline-producing strain. Overexpression of a feedback inhibition-resistant γ-glutamic kinase mutant ProBG149K, deletion of a proline dehydrogenase to block l-proline degradation, overexpression of glutamate dehydrogenase to increase glutamate synthesis flux, the mutation of 6-phosphate gluconate dehydrogenase and glucose-6-phosphate-dehydrogenase in the pentose phosphate pathway to enhance NADPH supply, the deletion of pyruvate aminotransferase to decrease the byproduct l-alanine synthesis, and weakening of α-ketoglutarate dehydrogenase to regulate the TCA cycle were combined to obtain ZQJY-9. ZQJY-9 produced 19.68 ± 0.22 g/L of l-proline in flask fermentation and was also demonstrated at the 3 L bioreactor level by fed-batch fermentation producing 120.18 g/L of l-proline at 76 h with the highest productivity of 1.581 g/L/h.

Covalent Modification of the Flavin in Proline Dehydrogenase by Thiazolidine-2-Carboxylate

Proline dehydrogenase (PRODH) catalyzes the first step of proline catabolism, the FAD-dependent 2-electron oxidation of l-proline to Δ¹-pyrroline-5-carboxylate. PRODH has emerged as a possible cancer therapy target, and thus the inhibition of PRODH is of interest. Here we show that the proline analogue thiazolidine-2-carboxylate (T2C) is a mechanism-based inactivator of PRODH. Structures of the bifunctional proline catabolic enzyme proline utilization A (PutA) determined from crystals grown in the presence of T2C feature strong electron density for a 5-membered ring species resembling l-T2C covalently bound to the N5 of the FAD in the PRODH domain. The modified FAD exhibits a large butterfly bend angle, indicating that the FAD is locked into the 2-electron reduced state. Reduction of the FAD is consistent with the crystals lacking the distinctive yellow color of the oxidized enzyme and stopped-flow kinetic data showing that T2C is a substrate for the PRODH domain of PutA. A mechanism is proposed in which PRODH catalyzes the oxidation of T2C at the C atom adjacent to the S atom of the thiazolidine ring (C5). Then, the N5 atom of the reduced FAD attacks the C5 of the oxidized T2C species, resulting in the covalent adduct observed in the crystal structure. To our knowledge, this is the first report of T2C inactivating (or inhibiting) PRODH or any other flavoenzyme. These results may inform the design of new mechanism-based inactivators of PRODH for use as chemical probes to study the roles of proline metabolism in cancer.

Role of Proline in Pathogen and Host Interactions

SIGNIFICANCE

Proline metabolism has complex roles in a variety of biological processes, including cell signaling, stress protection, and energy production. Proline also contributes to the pathogenesis of various disease-causing organisms. Understanding the mechanisms of how pathogens utilize proline is important for developing new strategies against infectious diseases. Recent Advances: The ability of pathogens to acquire amino acids is critical during infection. Besides protein biosynthesis, some amino acids, such as proline, serve as a carbon, nitrogen, or energy source in bacterial and protozoa pathogens. The role of proline during infection depends on the physiology of the host/pathogen interactions. Some pathogens rely on proline as a critical respiratory substrate, whereas others exploit proline for stress protection.

CRITICAL ISSUES

Disruption of proline metabolism and uptake has been shown to significantly attenuate virulence of certain pathogens, whereas in other pathogens the importance of proline during infection is not known. Inhibiting proline metabolism and transport may be a useful therapeutic strategy against some pathogens. Developing specific inhibitors to avoid off-target effects in the host, however, will be challenging. Also, potential treatments that target proline metabolism should consider the impact on intracellular levels of Δ1-pyrroline-5-carboxylate, a metabolite intermediate that can have opposing effects on pathogenesis.

FUTURE DIRECTIONS

Further characterization of how proline metabolism is regulated during infection would provide new insights into the role of proline in pathogenesis. Biochemical and structural characterization of proline metabolic enzymes from different pathogens could lead to new tools for exploring proline metabolism during infection and possibly new therapeutic compounds.

Structural Biology of Proline Catabolic Enzymes

SIGNIFICANCE

Proline catabolism refers to the 4-electron oxidation of proline to glutamate catalyzed by the enzymes proline dehydrogenase (PRODH) and l-glutamate γ-semialdehyde dehydrogenase (GSALDH, or ALDH4A1). These enzymes and the intermediate metabolites of the pathway have been implicated in tumor growth and suppression, metastasis, hyperprolinemia metabolic disorders, schizophrenia susceptibility, life span extension, and pathogen virulence and survival. In some bacteria, PRODH and GSALDH are combined into a bifunctional enzyme known as proline utilization A (PutA). PutAs are not only virulence factors in some pathogenic bacteria but also fascinating systems for studying the coordination of metabolic enzymes via substrate channeling. Recent Advances: The past decade has seen an explosion of structural data for proline catabolic enzymes. This review surveys these structures, emphasizing protein folds, substrate recognition, oligomerization, kinetic mechanisms, and substrate channeling in PutA.

CRITICAL ISSUES

Major unsolved structural targets include eukaryotic PRODH, the complex between monofunctional PRODH and monofunctional GSALDH, and the largest of all PutAs, trifunctional PutA. The structural basis of PutA-membrane association is poorly understood. Fundamental aspects of substrate channeling in PutA remain unknown, such as the identity of the channeled intermediate, how the tunnel system is activated, and the roles of ancillary tunnels.

FUTURE DIRECTIONS

New approaches are needed to study the molecular and in vivo mechanisms of substrate channeling. With the discovery of the proline cycle driving tumor growth and metastasis, the development of inhibitors of proline metabolic enzymes has emerged as an exciting new direction. Structural biology will be important in these endeavors.

The Proline Cycle As a Potential Cancer Therapy Target

Interest in how proline contributes to cancer biology is expanding because of the emerging role of a novel proline metabolic cycle in cancer cell survival, proliferation, and metastasis. Proline biosynthesis and degradation involve the shared intermediate Δ1-pyrroline-5-carboxylate (P5C), which forms l-glutamate-γ-semialdehyde (GSAL) in a reversible non-enzymatic reaction. Proline is synthesized from glutamate or ornithine through GSAL/P5C, which is reduced to proline by P5C reductase (PYCR) in a NAD(P)H-dependent reaction. The degradation of proline occurs in the mitochondrion and involves two oxidative steps catalyzed by proline dehydrogenase (PRODH) and GSAL dehydrogenase (GSALDH). PRODH is a flavin-dependent enzyme that couples proline oxidation with reduction of membrane-bound quinone, while GSALDH catalyzes the NAD+-dependent oxidation of GSAL to glutamate. PRODH and PYCR form a metabolic relationship known as the proline-P5C cycle, a novel pathway that impacts cellular growth and death pathways. The proline-P5C cycle has been implicated in supporting ATP production, protein and nucleotide synthesis, anaplerosis, and redox homeostasis in cancer cells. This Perspective details the structures and reaction mechanisms of PRODH and PYCR and the role of the proline-P5C cycle in cancer metabolism. A major challenge in the field is to discover inhibitors that specifically target PRODH and PYCR isoforms for use as tools for studying proline metabolism and the functions of the proline-P5C cycle in cancer. These molecular probes could also serve as lead compounds in cancer drug discovery targeting the proline-P5C cycle.

Discovery of the Membrane Binding Domain in Trifunctional Proline Utilization A

Escherichia coli proline utilization A (EcPutA) is the archetype of trifunctional PutA flavoproteins, which function both as regulators of the proline utilization operon and bifunctional enzymes that catalyze the four-electron oxidation of proline to glutamate. EcPutA shifts from a self-regulating transcriptional repressor to a bifunctional enzyme in a process known as functional switching. The flavin redox state dictates the function of EcPutA. Upon proline oxidation, the flavin becomes reduced, triggering a conformational change that causes EcPutA to dissociate from the put regulon and bind to the cellular membrane. Major structure/function domains of EcPutA have been characterized, including the DNA-binding domain, proline dehydrogenase (PRODH) and l-glutamate-γ-semialdehyde dehydrogenase catalytic domains, and an aldehyde dehydrogenase superfamily fold domain. Still lacking is an understanding of the membrane-binding domain, which is essential for EcPutA catalytic turnover and functional switching. Here, we provide evidence for a conserved C-terminal motif (CCM) in EcPutA having a critical role in membrane binding. Deletion of the CCM or replacement of hydrophobic residues with negatively charged residues within the CCM impairs EcPutA functional and physical membrane association. Furthermore, cell-based transcription assays and limited proteolysis indicate that the CCM is essential for functional switching. Using fluorescence resonance energy transfer involving dansyl-labeled liposomes, residues in the α-domain are also implicated in membrane binding. Taken together, these experiments suggest that the CCM and α-domain converge to form a membrane-binding interface near the PRODH domain. The discovery of the membrane-binding region will assist efforts to define flavin redox signaling pathways responsible for EcPutA functional switching.

Structure, function, and mechanism of proline utilization A (PutA)

Proline has important roles in multiple biological processes such as cellular bioenergetics, cell growth, oxidative and osmotic stress response, protein folding and stability, and redox signaling. The proline catabolic pathway, which forms glutamate, enables organisms to utilize proline as a carbon, nitrogen, and energy source. FAD-dependent proline dehydrogenase (PRODH) and NAD⁺-dependent glutamate semialdehyde dehydrogenase (GSALDH) convert proline to glutamate in two sequential oxidative steps. Depletion of PRODH and GSALDH in humans leads to hyperprolinemia, which is associated with mental disorders such as schizophrenia. Also, some pathogens require proline catabolism for virulence. A unique aspect of proline catabolism is the multifunctional proline utilization A (PutA) enzyme found in Gram-negative bacteria. PutA is a large (>1000 residues) bifunctional enzyme that combines PRODH and GSALDH activities into one polypeptide chain. In addition, some PutAs function as a DNA-binding transcriptional repressor of proline utilization genes. This review describes several attributes of PutA that make it a remarkable flavoenzyme: (1) diversity of oligomeric state and quaternary structure; (2) substrate channeling and enzyme hysteresis; (3) DNA-binding activity and transcriptional repressor function; and (4) flavin redox dependent changes in subcellular location and function in response to proline (functional switching).

Identification of a Conserved Histidine As Being Critical for the Catalytic Mechanism and Functional Switching of the Multifunctional Proline Utilization A Protein

Proline utilization A from Escherichia coli (EcPutA) is a multifunctional flavoenzyme that oxidizes proline to glutamate through proline dehydrogenase (PRODH) and Δ¹-pyrroline-5-carboxylate dehydrogenase (P5CDH) activities, while also switching roles as a DNA-bound transcriptional repressor and a membrane-bound catabolic enzyme. This phenomenon, termed functional switching, occurs through a redox-mediated mechanism in which flavin reduction triggers a conformational change that increases EcPutA membrane binding affinity. Structural studies have shown that reduction of the FAD cofactor causes the ribityl moiety to undergo a crankshaft motion, indicating that the orientation of the ribityl chain is a key element of PutA functional switching. Here, we test the role of a conserved histidine that bridges from the FAD pyrophosphate to the backbone amide of a conserved leucine residue in the PRODH active site. An EcPutA mutant (H487A) was characterized by steady-state and rapid-reaction kinetics, and cell-based reporter gene experiments. The catalytic activity of H487A is severely diminished (>50-fold) with membrane vesicles as the electron acceptor, and H487A exhibits impaired lipid binding and in vivo transcriptional repressor activity. Rapid-reaction kinetic experiments demonstrate that H487A is 3-fold slower than wild-type EcPutA in a conformational change step following reduction of the FAD cofactor. Furthermore, the reduction potential (E_m) of H487A is ∼40 mV more positive than that of wild-type EcPutA, and H487A has an attenuated ability to catalyze the reverse PRODH chemical step of reoxidation by P5C. In this process, significant red semiquinone forms in contrast to the same reaction with wild-type EcPutA, in which facile two-electron reoxidation occurs without the formation of a measurable amount of semiquinone. These results indicate that His487 is critically important for the proline/P5C chemical step, conformational change kinetics, and functional switching in EcPutA.

A new genome-scale metabolic model of Corynebacterium glutamicum and its application

BACKGROUND

Corynebacterium glutamicum is an important platform organism for industrial biotechnology to produce amino acids, organic acids, bioplastic monomers, and biofuels. The metabolic flexibility, broad substrate spectrum, and fermentative robustness of C. glutamicum make this organism an ideal cell factory to manufacture desired products. With increases in gene function, transport system, and metabolic profile information under certain conditions, developing a comprehensive genome-scale metabolic model (GEM) of C. glutamicum ATCC13032 is desired to improve prediction accuracy, elucidate cellular metabolism, and guide metabolic engineering.

RESULTS

Here, we constructed a new GEM for ATCC13032, iCW773, consisting of 773 genes, 950 metabolites, and 1207 reactions. Compared to the previous model, iCW773 supplemented 496 gene-protein-reaction associations, refined five lumped reactions, balanced the mass and charge, and constrained the directionality of reactions. The simulated growth rates of C. glutamicum cultivated on seven different carbon sources using iCW773 were consistent with experimental values. Pearson's correlation coefficient between the iCW773-simulated and experimental fluxes was 0.99, suggesting that iCW773 provided an accurate intracellular flux distribution of the wild-type strain growing on glucose. Furthermore, genetic interventions for overproducing l-lysine, 1,2-propanediol and isobutanol simulated using OptForce_MUST were in accordance with reported experimental results, indicating the practicability of iCW773 for the design of metabolic networks to overproduce desired products. In vivo genetic modifications of iCW773-predicted targets resulted in the de novo generation of an l-proline-overproducing strain. In fed-batch culture, the engineered C. glutamicum strain produced 66.43 g/L l-proline in 60 h with a yield of 0.26 g/g (l-proline/glucose) and a productivity of 1.11 g/L/h. To our knowledge, this is the highest titer and productivity reported for l-proline production using glucose as the carbon resource in a minimal medium.

CONCLUSIONS

Our developed iCW773 serves as a high-quality platform for model-guided strain design to produce industrial bioproducts of interest. This new GEM will be a successful multidisciplinary tool and will make valuable contributions to metabolic engineering in academia and industry.

Visualization of Metabolic Interaction Networks in Microbial Communities Using VisANT 5.0

The complexity of metabolic networks in microbial communities poses an unresolved visualization and interpretation challenge. We address this challenge in the newly expanded version of a software tool for the analysis of biological networks, VisANT 5.0. We focus in particular on facilitating the visual exploration of metabolic interaction between microbes in a community, e.g. as predicted by COMETS (Computation of Microbial Ecosystems in Time and Space), a dynamic stoichiometric modeling framework. Using VisANT's unique metagraph implementation, we show how one can use VisANT 5.0 to explore different time-dependent ecosystem-level metabolic networks. In particular, we analyze the metabolic interaction network between two bacteria previously shown to display an obligate cross-feeding interdependency. In addition, we illustrate how a putative minimal gut microbiome community could be represented in our framework, making it possible to highlight interactions across multiple coexisting species. We envisage that the "symbiotic layout" of VisANT can be employed as a general tool for the analysis of metabolism in complex microbial communities as well as heterogeneous human tissues. VisANT is freely available at: http://visant.bu.edu and COMETS at http://comets.bu.edu.

Analytical methods for kinetic studies of biological interactions: A review

The rates at which biological interactions occur can provide important information concerning the mechanism and behavior of these processes in living systems. This review discusses several analytical methods that can be used to examine the kinetics of biological interactions. These techniques include common or traditional methods such as stopped-flow analysis and surface plasmon resonance spectroscopy, as well as alternative methods based on affinity chromatography and capillary electrophoresis. The general principles and theory behind these approaches are examined, and it is shown how each technique can be utilized to provide information on the kinetics of biological interactions. Examples of applications are also given for each method. In addition, a discussion is provided on the relative advantages or potential limitations of each technique regarding its use in kinetic studies.

Proline metabolism increases katG expression and oxidative stress resistance in Escherichia coli

The oxidation of l-proline to glutamate in Gram-negative bacteria is catalyzed by the proline utilization A (PutA) flavoenzyme, which contains proline dehydrogenase (PRODH) and Δ(1)-pyrroline-5-carboxylate (P5C) dehydrogenase domains in a single polypeptide. Previous studies have suggested that aside from providing energy, proline metabolism influences oxidative stress resistance in different organisms. To explore this potential role and the mechanism, we characterized the oxidative stress resistance of wild-type and putA mutant strains of Escherichia coli. Initial stress assays revealed that the putA mutant strain was significantly more sensitive to oxidative stress than the parental wild-type strain. Expression of PutA in the putA mutant strain restored oxidative stress resistance, confirming that depletion of PutA was responsible for the oxidative stress phenotype. Treatment of wild-type cells with proline significantly increased hydroperoxidase I (encoded by katG) expression and activity. Furthermore, the ΔkatG strain failed to respond to proline, indicating a critical role for hydroperoxidase I in the mechanism of proline protection. The global regulator OxyR activates the expression of katG along with several other genes involved in oxidative stress defense. In addition to katG, proline increased the expression of grxA (glutaredoxin 1) and trxC (thioredoxin 2) of the OxyR regulon, implicating OxyR in proline protection. Proline oxidative metabolism was shown to generate hydrogen peroxide, indicating that proline increases oxidative stress tolerance in E. coli via a preadaptive effect involving endogenous hydrogen peroxide production and enhanced catalase-peroxidase activity.

Characterization of the kinetics and electron paramagnetic resonance spectroscopic properties of Acidithiobacillus ferrooxidans sulfide:quinone oxidoreductase (SQR)

Acidithiobacillus ferrooxidans sulfide:quinone oxidoreductase (SQR) catalyzes the oxidation of sulfide to polysulfide chains or elemental sulfur coupled to quinone reduction via a non-covalent FAD cofactor. We investigated the role of the FAD using kinetics and EPR spectroscopy. The properties of the enzyme were compared with alanine and/or serine variants of conserved cysteine residues (Cys128, Cys160, Cys356) structurally close to the FAD cofactor and histidine residues (His132, His198) implicated in function. When the pre-steady state reduction of FAD was monitored, variants of Cys128 and His132 had similar rates to wild-type enzyme confirming they do not participate in the reductive half reaction whereas variants of Cys160, Cys356 and His198 had greatly reduced activity. Using steady state kinetics of Na2S-dependent decylubiquinone (DUQ) reduction we measured a kcat of 6.5s(-1) and a Km (Na2S) of 3.0μM and a Km (DUQ) of 3.4μM. Variants of Cys160, Cys356 and His198 had greatly diminished DUQ reduction activity whereas variants of Cys128 and His132 were less affected. A neutral flavin semiquinone was observed in the EPR spectrum of SQR reduced with Na2S which was enhanced in the Cys160Ala variant suggesting the presence of a Cys356-S(γ)-S-C(4A)-FAD adduct. Potentiometric titrations of the FAD semiquinone revealed an Em of -139±4mV at pH 7.0.

Kinetic and structural characterization of tunnel-perturbing mutants in Bradyrhizobium japonicum proline utilization A

Proline utilization A from Bradyrhizobium japonicum (BjPutA) is a bifunctional flavoenzyme that catalyzes the oxidation of proline to glutamate using fused proline dehydrogenase (PRODH) and Δ(1)-pyrroline-5-carboxylate dehydrogenase (P5CDH) domains. Recent crystal structures and kinetic data suggest an intramolecular channel connects the two active sites, promoting substrate channeling of the intermediate Δ(1)-pyrroline-5-carboxylate/glutamate-γ-semialdehyde (P5C/GSA). In this work, the structure of the channel was explored by inserting large side chain residues at four positions along the channel in BjPutA. Kinetic analysis of the different mutants revealed replacement of D779 with Tyr (D779Y) or Trp (D779W) significantly decreased the overall rate of the PRODH-P5CDH channeling reaction. X-ray crystal structures of D779Y and D779W revealed that the large side chains caused a constriction in the central section of the tunnel, thus likely impeding the travel of P5C/GSA in the channel. The D779Y and D779W mutants have PRODH activity similar to that of wild-type BjPutA but exhibit significantly lower P5CDH activity, suggesting that exogenous P5C/GSA enters the channel upstream of Asp779. Replacement of nearby Asp778 with Tyr (D778Y) did not impact BjPutA channeling activity. Consistent with the kinetic results, the X-ray crystal structure of D778Y shows that the main channel pathway is not impacted; however, an off-cavity pathway is closed off from the channel. These findings provide evidence that the off-cavity pathway is not essential for substrate channeling in BjPutA.

Structures of the PutA peripheral membrane flavoenzyme reveal a dynamic substrate-channeling tunnel and the quinone-binding site

Proline utilization A (PutA) proteins are bifunctional peripheral membrane flavoenzymes that catalyze the oxidation of L-proline to L-glutamate by the sequential activities of proline dehydrogenase and aldehyde dehydrogenase domains. Located at the inner membrane of Gram-negative bacteria, PutAs play a major role in energy metabolism by coupling the oxidation of proline imported from the environment to the reduction of membrane-associated quinones. Here, we report seven crystal structures of the 1,004-residue PutA from Geobacter sulfurreducens, along with determination of the protein oligomeric state by small-angle X-ray scattering and kinetic characterization of substrate channeling and quinone reduction. The structures reveal an elaborate and dynamic tunnel system featuring a 75-Å-long tunnel that links the two active sites and six smaller tunnels that connect the main tunnel to the bulk medium. The locations of these tunnels and their responses to ligand binding and flavin reduction suggest hypotheses about how proline, water, and quinones enter the tunnel system and where L-glutamate exits. Kinetic measurements show that glutamate production from proline occurs without a lag phase, consistent with substrate channeling and implying that the observed tunnel is functionally relevant. Furthermore, the structure of reduced PutA complexed with menadione bisulfite reveals the elusive quinone-binding site. The benzoquinone binds within 4.0 Å of the flavin si face, consistent with direct electron transfer. The location of the quinone site implies that the concave surface of the PutA dimer approaches the membrane. Altogether, these results provide insight into how PutAs couple proline oxidation to quinone reduction.

Evidence for hysteretic substrate channeling in the proline dehydrogenase and Δ1-pyrroline-5-carboxylate dehydrogenase coupled reaction of proline utilization A (PutA)

PutA (proline utilization A) is a large bifunctional flavoenzyme with proline dehydrogenase (PRODH) and Δ(1)-pyrroline-5-carboxylate dehydrogenase (P5CDH) domains that catalyze the oxidation of l-proline to l-glutamate in two successive reactions. In the PRODH active site, proline undergoes a two-electron oxidation to Δ(1)-pyrroline-5-carboxlylate, and the FAD cofactor is reduced. In the P5CDH active site, l-glutamate-γ-semialdehyde (the hydrolyzed form of Δ(1)-pyrroline-5-carboxylate) undergoes a two-electron oxidation in which a hydride is transferred to NAD(+)-producing NADH and glutamate. Here we report the first kinetic model for the overall PRODH-P5CDH reaction of a PutA enzyme. Global analysis of steady-state and transient kinetic data for the PRODH, P5CDH, and coupled PRODH-P5CDH reactions was used to test various models describing the conversion of proline to glutamate by Escherichia coli PutA. The coupled PRODH-P5CDH activity of PutA is best described by a mechanism in which the intermediate is not released into the bulk medium, i.e., substrate channeling. Unexpectedly, single-turnover kinetic experiments of the coupled PRODH-P5CDH reaction revealed that the rate of NADH formation is 20-fold slower than the steady-state turnover number for the overall reaction, implying that catalytic cycling speeds up throughput. We show that the limiting rate constant observed for NADH formation in the first turnover increases by almost 40-fold after multiple turnovers, achieving half of the steady-state value after 15 turnovers. These results suggest that EcPutA achieves an activated channeling state during the approach to steady state and is thus a new example of a hysteretic enzyme. Potential underlying causes of activation of channeling are discussed.

Essential role of conformational selection in ligand binding

Two competing and mutually exclusive mechanisms of ligand recognition - conformational selection and induced fit - have dominated our interpretation of ligand binding in biological macromolecules for almost six decades. Conformational selection posits the pre-existence of multiple conformations of the macromolecule from which the ligand selects the optimal one. Induced fit, on the other hand, postulates the existence of conformational rearrangements of the original conformation into an optimal one that are induced by binding of the ligand. In the former case, conformational transitions precede the binding event; in the latter, conformational changes follow the binding step. Kineticists have used a facile criterion to distinguish between the two mechanisms based on the dependence of the rate of relaxation to equilibrium, kobs, on the ligand concentration, [L]. A value of kobs decreasing hyperbolically with [L] has been seen as diagnostic of conformational selection, while a value of kobs increasing hyperbolically with [L] has been considered diagnostic of induced fit. However, this simple conclusion is only valid under the rather unrealistic assumption of conformational transitions being much slower than binding and dissociation events. In general, induced fit only produces values of kobs that increase with [L] but conformational selection is more versatile and is associated with values of kobs that increase with, decrease with or are independent of [L]. The richer repertoire of kinetic properties of conformational selection applies to kinetic mechanisms with single or multiple saturable relaxations and explains the behavior of nearly all experimental systems reported in the literature thus far. Conformational selection is always sufficient and often necessary to account for the relaxation kinetics of ligand binding to a biological macromolecule and is therefore an essential component of any binding mechanism. On the other hand, induced fit is never necessary and only sufficient in a few cases. Therefore, the long assumed importance and preponderance of induced fit as a mechanism of ligand binding should be reconsidered.

Conformational selection is a dominant mechanism of ligand binding

Molecular recognition in biological macromolecules is achieved by binding interactions coupled to conformational transitions that precede or follow the binding step, two limiting mechanisms known as conformational selection and induced fit, respectively. Sorting out the contribution of these mechanisms to any binding interaction remains a challenging task of general interest in biochemistry. Here we show that conformational selection is associated with a vast repertoire of kinetic behaviors, can never be disproved a priori as a mechanism of ligand binding, and is sufficient to explain the relaxation kinetics documented experimentally for a large number of systems. On the other hand, induced fit features a narrow spectrum of kinetic behaviors and can be disproved in many cases in which conformational selection offers the only possible explanation. This conclusion offers a paradigm shift in the analysis of relaxation kinetics, with conformational selection acquiring preeminence as a mechanism of ligand binding. The dominant role of conformational selection supports the emerging structural view of the macromolecule as a conformational ensemble from which the ligand selects the initial optimal fit to produce a biological response.

Kinetic and isotopic characterization of L-proline dehydrogenase from Mycobacterium tuberculosis

The monofunctional proline dehydrogenase (ProDH) from Mycobacterium tuberculosis performs the flavin-dependent oxidation of l-proline to Δ(1)-pyrroline-5-carboxylate in the proline catabolic pathway. The ProDH gene, prub, was cloned into the pYUB1062 vector, and the C-terminal His-tagged 37 kDa protein was expressed and purified by nickel affinity chromatography. A steady-state kinetic analysis revealed a ping-pong mechanism with an overall kcat of 33 ± 2 s(-1) and Km values of 5.7 ± 0.8 mM and 3.4 ± 0.3 μM for l-proline and 2,6-dichlorophenolindophenol (DCPIP), respectively. The pH dependence of kcat revealed that one enzyme group exhibiting a pK value of 6.8 must be deprotonated for optimal catalytic activity. Site-directed mutagenesis suggests that this group is Lys110. The primary kinetic isotope effects on V/KPro and V of 5.5 and 1.1, respectively, suggest that the transfer of hydride from l-proline to FAD is rate-limiting for the reductive half-reaction, but that FAD reoxidation is the rate-limiting step in the overall reaction. Solvent and multiple kinetic isotope effects suggest that l-proline oxidation occurs in a stepwise rather than concerted mechanism. Pre-steady-state kinetics reveal an overall kred of 88.5 ± 0.7 s(-1), and this rate is subject to a primary kinetic isotope effect of 5.2. These data confirm that the overall reaction is limited by reduced flavin reoxidation in the second half-reaction.

Ornithine cyclodeaminase-based proline production by Corynebacterium glutamicum

BACKGROUND

The soil bacterium Corynebacterium glutamicum, best known for its glutamate producing ability, is suitable as a producer of a variety of bioproducts. Glutamate is the precursor of the amino acid proline. Proline biosynthesis typically involves three enzymes and a spontaneous cyclisation reaction. Alternatively, proline can be synthesised from ornithine, an intermediate of arginine biosynthesis. The direct conversion of ornithine to proline is catalysed by ornithine cyclodeaminase. An ornithine overproducing platform strain with deletions of argR and argF (ORN1) has been employed for production of derived compounds such as putrescine. By heterologous expression of ocd this platform strain can be engineered further for proline production.

RESULTS

Plasmid-based expression of ocd encoding the putative ornithine cyclodeaminase of C. glutamicum did not result in detectable proline accumulation in the culture medium. However, plasmid-based expression of ocd from Pseudomonas putida resulted in proline production with yields up to 0.31 ± 0.01 g proline/g glucose. Overexpression of the gene encoding a feedback-alleviated N-acetylglutamate kinase further increased proline production to 0.36 ± 0.01 g/g. In addition, feedback-alleviation of N-acetylglutamate kinase entailed growth-coupled production of proline and reduced the accumulation of by-products in the culture medium.

CONCLUSIONS

The product spectrum of the platform strain C. glutamicum ORN1 was expanded to include the amino acid L-proline. Upon further development of the ornithine overproducing platform strain, industrial production of amino acids of the glutamate family and derived bioproducts such as diamines might become within reach.

Crystal structures and kinetics of monofunctional proline dehydrogenase provide insight into substrate recognition and conformational changes associated with flavin reduction and product release

Proline dehydrogenase (PRODH) catalyzes the FAD-dependent oxidation of proline to Δ(1)-pyrroline-5-carboxylate, which is the first step of proline catabolism. Here, we report the structures of proline dehydrogenase from Deinococcus radiodurans in the oxidized state complexed with the proline analogue L-tetrahydrofuroic acid and in the reduced state with the proline site vacant. The analogue binds against the si face of the FAD isoalloxazine and is protected from bulk solvent by helix α8 and the β1-α1 loop. The FAD ribityl chain adopts two conformations in the E-S complex, which is unprecedented for flavoenzymes. One of the conformations is novel for the PRODH superfamily and may contribute to the low substrate affinity of Deinococcus PRODH. Reduction of the crystalline enzyme-inhibitor complex causes profound structural changes, including 20° butterfly bending of the isoalloxazine, crankshaft rotation of the ribityl, shifting of α8 by 1.7 Å, reconfiguration of the β1-α1 loop, and rupture of the Arg291-Glu64 ion pair. These changes dramatically open the active site to facilitate product release and allow electron acceptors access to the reduced flavin. The structures suggest that the ion pair, which is conserved in the PRODH superfamily, functions as the active site gate. Mutagenesis of Glu64 to Ala decreases the catalytic efficiency 27-fold, which demonstrates the importance of the gate. Mutation of Gly63 decreases the efficiency 140-fold, which suggests that flexibility of the β1-α1 loop is essential for optimal catalysis. The large conformational changes that are required to form the E-S complex suggest that conformational selection plays a role in substrate recognition.