Evolutionary analysis of the TPP-dependent enzyme family

Biosynthesis of the benzylpyrrolidine precursor in anisomycin by a unique ThDP-dependent enzyme

Anisomycin (compound 1), a multifunctional pyrrolidine antibiotic, primarily inhibits protein biosynthesis by binding to the ribosome. Upon binding to the ribosome, the para-phenol moiety of anisomycin inserts completely into the hydrophobic crevice of the A-site and blocks the access of the incoming aminoacyl-tRNAs, disrupting peptide bond formation. Hence, the para-methoxyphenyl group serves as a starting point for developing novel anisomycin analogs with potent antifungal and insecticidal properties. However, the activation and condensation mechanism of phenylpyruvic acid has not yet been elucidated. In this study, genetic manipulations of aniP and its homologue siAniP confirmed their indispensable role in 1 biosynthesis. Bioinformatics analysis suggested that AniP and siAniP function as transketolase. siAniP was found to catalyzed condensation between 4-hydroxyphenylpyruvic acid (3) and glyceraldehyde (GA), initiating pyrrolidine synthesis. siAniP was specific for aromatic keto acids and tolerant of aliphatic and aromatic aldehydes, and was able to catalyze the asymmetric intermolecular condensation of two keto acids, leading to the formation of 24 α-hydroxy ketone. To the best of our knowledge, siAniP is the first TK that catalyzes the transfer of a C2 ketol and symmetrical intermolecular coupling using aromatic keto acids as donor substrates. Structural analysis, docking model construction, and site-directed mutagenesis identified that I220, H275, R322 and W391 were crucial for substrate binding. Moreover, sequence similarity network (SSN)-based genome neighborhood network (GNN) analyses of AniP suggested the widespread occurrence of the AniP-like-mediated reaction in the biosynthesis of 1 and its analogs, particularly in the assembly of benzylpyrrolidine. These findings not only expand the repertoire of TKs but also provide a potent biocatalyst that could be used for the structural innovation of 1 and its derivatives.

A Novel Pathway for Biosynthesis of the Herbicidal Phosphonate Natural Product Phosphonothrixin Is Widespread in Actinobacteria

Phosphonothrixin is an herbicidal phosphonate natural product with an unusual, branched carbon skeleton. Bioinformatic analyses of the ftx gene cluster, which is responsible for synthesis of the compound, suggest that early steps of the biosynthetic pathway, up to production of the intermediate 2,3-dihydroxypropylphosphonic acid (DHPPA) are identical to those of the unrelated phosphonate natural product valinophos. This conclusion was strongly supported by the observation of biosynthetic intermediates from the shared pathway in spent media from two phosphonothrixin producing strains. Biochemical characterization of ftx-encoded proteins confirmed these early steps, as well as subsequent steps involving the oxidation of DHPPA to 3-hydroxy-2-oxopropylphosphonate and its conversion to phosphonothrixin by the combined action of an unusual heterodimeric, thiamine-pyrophosphate (TPP)-dependent ketotransferase and a TPP-dependent acetolactate synthase. The frequent observation of ftx-like gene clusters within actinobacteria suggests that production of compounds related to phosphonothrixin is common within these bacteria. IMPORTANCE Phosphonic acid natural products, such as phosphonothrixin, have great potential for biomedical and agricultural applications; however, discovery and development of these compounds requires detailed knowledge of the metabolism involved in their biosynthesis. The studies reported here reveal the biochemical pathway phosphonothrixin production, which enhances our ability to design strains that overproduce this potentially useful herbicide. This knowledge also improves our ability to predict the products of related biosynthetic gene clusters and the functions of homologous enzymes.

Two distinct gut microbial pathways contribute to meta-organismal production of phenylacetylglutamine with links to cardiovascular disease

Recent studies show gut microbiota-dependent metabolism of dietary phenylalanine into phenylacetic acid (PAA) is critical in phenylacetylglutamine (PAGln) production, a metabolite linked to atherosclerotic cardiovascular disease (ASCVD). Accordingly, microbial enzymes involved in this transformation are of interest. Using genetic manipulation in selected microbes and monocolonization experiments in gnotobiotic mice, we identify two distinct gut microbial pathways for PAA formation; one is catalyzed by phenylpyruvate:ferredoxin oxidoreductase (PPFOR) and the other by phenylpyruvate decarboxylase (PPDC). PPFOR and PPDC play key roles in gut bacterial PAA production via oxidative and non-oxidative phenylpyruvate decarboxylation, respectively. Metagenomic analyses revealed a significantly higher abundance of both pathways in gut microbiomes of ASCVD patients compared with controls. The present studies show a role for these two divergent microbial catalytic strategies in the meta-organismal production of PAGln. Given the numerous links between PAGln and ASCVD, these findings will assist future efforts to therapeutically target PAGln formation in vivo.

Characterisation of plasmodial transketolases and identification of potential inhibitors: an in silico study

BACKGROUND

Plasmodial transketolase (PTKT) enzyme is one of the novel pharmacological targets being explored as potential anti-malarial drug target due to its functional role and low sequence identity to the human enzyme. Despite this, features contributing to such have not been exploited for anti-malarial drug design. Additionally, there are no anti-malarial drugs targeting PTKTs whereas the broad activity of these inhibitors against PTKTs from other Plasmodium spp. is yet to be reported. This study characterises different PTKTs [Plasmodium falciparum (PfTKT), Plasmodium vivax (PvTKT), Plasmodium ovale (PoTKT), Plasmodium malariae (PmTKT) and Plasmodium knowlesi (PkTKT) and the human homolog (HsTKT)] to identify key sequence and structural based differences as well as the identification of selective potential inhibitors against PTKTs.

METHODS

A sequence-based study was carried out using multiple sequence alignment, phylogenetic tree calculations and motif discovery analysis. Additionally, TKT models of PfTKT, PmTKT, PoTKT, PmTKT and PkTKT were modelled using the Saccharomyces cerevisiae TKT structure as template. Based on the modelled structures, molecular docking using 623 South African natural compounds was done. The stability, conformational changes and detailed interactions of selected compounds were accessed viz all-atom molecular dynamics (MD) simulations and binding free energy (BFE) calculations.

RESULTS

Sequence alignment, evolutionary and motif analyses revealed key differences between plasmodial and the human TKTs. High quality homodimeric three-dimensional PTKTs structures were constructed. Molecular docking results identified three compounds (SANC00107, SANC00411 and SANC00620) which selectively bind in the active site of all PTKTs with the lowest (better) binding affinity ≤ - 8.5 kcal/mol. MD simulations of ligand-bound systems showed stable fluctuations upon ligand binding. In all systems, ligands bind stably throughout the simulation and form crucial interactions with key active site residues. Simulations of selected compounds in complex with human TKT showed that ligands exited their binding sites at different time steps. BFE of protein-ligand complexes showed key residues involved in binding.

CONCLUSIONS

This study highlights significant differences between plasmodial and human TKTs and may provide valuable information for the development of novel anti-malarial inhibitors. Identified compounds may provide a starting point in the rational design of PTKT inhibitors and analogues based on these scaffolds.

A 'Split-Gene' Transketolase From the Hyper-Thermophilic Bacterium Carboxydothermus hydrogenoformans: Structure and Biochemical Characterization

A novel transketolase has been reconstituted from two separate polypeptide chains encoded by a ‘split-gene’ identified in the genome of the hyperthermophilic bacterium, Carboxydothermus hydrogenoformans. The reconstituted active α₂β₂ tetrameric enzyme has been biochemically characterized and its activity has been determined using a range of aldehydes including glycolaldehyde, phenylacetaldehyde and cyclohexanecarboxaldehyde as the ketol acceptor and hydroxypyruvate as the donor. This reaction proceeds to near 100% completion due to the release of the product carbon dioxide and can be used for the synthesis of a range of sugars of interest to the pharmaceutical industry. This novel reconstituted transketolase is thermally stable with no loss of activity after incubation for 1 h at 70°C and is stable after 1 h incubation with 50% of the organic solvents methanol, ethanol, isopropanol, DMSO, acetonitrile and acetone. The X-ray structure of the holo reconstituted α₂β₂ tetrameric transketolase has been determined to 1.4 Å resolution. In addition, the structure of an inactive tetrameric β₄ protein has been determined to 1.9 Å resolution. The structure of the active reconstituted α₂β₂ enzyme has been compared to the structures of related enzymes; the E1 component of the pyruvate dehydrogenase complex and D-xylulose-5-phosphate synthase, in an attempt to rationalize differences in structure and substrate specificity between these enzymes. This is the first example of a reconstituted ‘split-gene’ transketolase to be biochemically and structurally characterized allowing its potential for industrial biocatalysis to be evaluated.

X-ray crystallography-based structural elucidation of enzyme-bound intermediates along the 1-deoxy-d-xylulose 5-phosphate synthase reaction coordinate

1-Deoxy-d-xylulose 5-phosphate synthase (DXPS) uses thiamine diphosphate (ThDP) to convert pyruvate and d-glyceraldehyde 3-phosphate (d-GAP) into 1-deoxy-d-xylulose 5-phosphate (DXP), an essential bacterial metabolite. DXP is not utilized by humans; hence, DXPS has been an attractive antibacterial target. Here, we investigate DXPS from Deinococcus radiodurans (DrDXPS), showing that it has similar kinetic parameters K_m^d-GAP and K_m^pyruvate (54 ± 3 and 11 ± 1 μm, respectively) and comparable catalytic activity (k_cat = 45 ± 2 min^-1) with previously studied bacterial DXPS enzymes and employing it to obtain missing structural data on this enzyme family. In particular, we have determined crystallographic snapshots of DrDXPS in two states along the reaction coordinate: a structure of DrDXPS bound to C2α-phosphonolactylThDP (PLThDP), mimicking the native pre-decarboxylation intermediate C2α-lactylThDP (LThDP), and a native post-decarboxylation state with a bound enamine intermediate. The 1.94-Å-resolution structure of PLThDP-bound DrDXPS delineates how two active-site histidine residues stabilize the LThDP intermediate. Meanwhile, the 2.40-Å-resolution structure of an enamine intermediate-bound DrDXPS reveals how a previously unknown 17-Å conformational change removes one of the two histidine residues from the active site, likely triggering LThDP decarboxylation to form the enamine intermediate. These results provide insight into how the bi-substrate enzyme DXPS limits side reactions by arresting the reaction on the less reactive LThDP intermediate when its cosubstrate is absent. They also offer a molecular basis for previous low-resolution experimental observations that correlate decarboxylation of LThDP with protein conformational changes.

Biosynthesis of Diverse Antimicrobial and Antiproliferative Acyloins in Anaerobic Bacteria

Metabolic profiling and genome mining revealed that anaerobic bacteria have the potential to produce acyloin natural products. In addition to sattazolin A and B, three new sattazolin congeners and a novel acyloin named clostrocyloin were isolated from three strains of Clostridium beijerinckii, a bacterium used for industrial solvent production. Bioactivity profiling showed that the sattazolin derivatives possess antimicrobial activities against mycobacteria and pseudomonads with only low cytotoxicity. Clostrocyloin was found to be mainly active against fungi. The thiamine diphosphate (ThDP)-dependent sattazolin-producing synthase was identified in silico and characterized both in vivo and in in vitro enzyme assays. A related acyloin synthase from the clostrocyloin producer was shown to be responsible for the production of the acyloin core of clostrocyloin. The biotransformation experiments provided first insights into the substrate scope of the clostrocyloin synthase and revealed biosynthetic intermediates.

Navigating within thiamine diphosphate-dependent decarboxylases: Sequences, structures, functional positions, and binding sites

Thiamine diphosphate-dependent decarboxylases catalyze both cleavage and formation of CC bonds in various reactions, which have been assigned to different homologous sequence families. This work compares 53 ThDP-dependent decarboxylases with known crystal structures. Both sequence and structural information were analyzed synergistically and data were analyzed for global and local properties by means of statistical approaches (principle component analysis and principal coordinate analysis) enabling complexity reduction. The different results obtained both locally and globally, that is, individual positions compared with the overall protein sequence or structure, revealed challenges in the assignment of separated homologous families. The methods applied herein support the comparison of enzyme families and the identification of functionally relevant positions. The findings for the family of ThDP-dependent decarboxylases underline that global sequence identity alone is not sufficient to distinguish enzyme function. Instead, local sequence similarity, defined by comparisons of structurally equivalent positions, allows for a better navigation within several groups of homologous enzymes. The differentiation between homologous sequences is further enhanced by taking structural information into account, such as BioGPS analysis of the active site properties or pairwise structural superimpositions. The methods applied herein are expected to be transferrable to other enzyme families, to facilitate family assignments for homologous protein sequences.

Exploiting correlated molecular-dynamics networks to counteract enzyme activity-stability trade-off

Rigidifying flexible sites is a powerful method to improve enzyme stability. However, if the highly flexible regions form the active site, modifying them risks losing activity due to the activity–stability trade-off. We hypothesized here that regions outside the active site whose dynamics were highly correlated to flexible active sites, would provide good targets for stabilizing mutations. To test this hypothesis, six variants were constructed in the 3M variant of Escherichia coli transketolase. The best variant had a 10.8-fold improved half-life at 55 °C, and increased the T_m and T_agg by 3 °C and 4.3 °C, respectively. The variants even increased the activity, by up to threefold. This study highlights how protein engineering strategies could be potentially improved by considering long-range dynamics.

The directed evolution of enzymes for improved activity or substrate specificity commonly leads to a trade-off in stability. We have identified an activity–stability trade-off and a loss in unfolding cooperativity for a variant (3M) of Escherichia coli transketolase (TK) engineered to accept aromatic substrates. Molecular dynamics simulations of 3M revealed increased flexibility in several interconnected active-site regions that also form part of the dimer interface. Mutating the newly flexible active-site residues to regain stability risked losing the new activity. We hypothesized that stabilizing mutations could be targeted to residues outside of the active site, whose dynamics were correlated with the newly flexible active-site residues. We previously stabilized WT TK by targeting mutations to highly flexible regions. These regions were much less flexible in 3M and would not have been selected a priori as targets using the same strategy based on flexibility alone. However, their dynamics were highly correlated with the newly flexible active-site regions of 3M. Introducing the previous mutations into 3M reestablished the WT level of stability and unfolding cooperativity, giving a 10.8-fold improved half-life at 55 °C, and increased midpoint and aggregation onset temperatures by 3 °C and 4.3 °C, respectively. Even the activity toward aromatic aldehydes increased up to threefold. Molecular dynamics simulations confirmed that the mutations rigidified the active-site via the correlated network. This work provides insights into the impact of rigidifying mutations within highly correlated dynamic networks that could also be useful for developing improved computational protein engineering strategies.

Optimisation of enzyme cascades for chiral amino alcohol synthesis in aid of host cell integration using a statistical experimental design approach

Chiral amino alcohols are compounds of pharmaceutical interest as they are building blocks of sphingolipids, antibiotics, and antiviral glycosidase inhibitors. Due to the challenges of chemical synthesis we recently developed two TK-TAm reaction cascades using natural and low cost feedstocks as substrates: a recycling cascade comprising of 2 enzymes and a sequential 3-step enzyme cascade yielding 30% and 1% conversion, respectively. In order to improve the conversion yield and aid the future host strain engineering for whole cell biocatalysis, we used a combination of microscale experiments and statistical experimental design. For this we implemented a full factorial design to optimise pH, temperature and buffer type, followed by the application of Response Surface Methodology for the optimisation of substrates and enzymes concentrations. Using purified enzymes we achieved 60% conversion for the recycling cascade and 3-fold improvement using the sequential pathway. Based on the results, limiting steps and individual requirements for host cell metabolic integration were identified expanding the understanding of the cascades without implementing extensive optimisation modelling. Therefore, the approach described here is well suited for optimising reaction conditions as well as defining the relative enzyme expression levels required for construction of microbial cell factories.

Binding site for coenzyme A revealed in the structure of pyruvate:ferredoxin oxidoreductase from Moorella thermoacetica

Pyruvate:ferredoxin oxidoreductase (PFOR) is a microbial enzyme that uses thiamine pyrophosphate (TPP), three [4Fe-4S] clusters, and coenzyme A (CoA) in the reversible oxidation of pyruvate to generate acetyl-CoA and carbon dioxide. The two electrons that are generated as a result of pyruvate decarboxylation are used in the reduction of low potential ferredoxins, which provide reducing equivalents for central metabolism, including the Wood-Ljungdahl pathway. PFOR is a member of the 2-oxoacid:ferredoxin oxidoreductase (OFOR) superfamily, which plays major roles in both microbial redox reactions and carbon dioxide fixation. Here, we present a set of crystallographic snapshots of the best-studied member of this superfamily, the PFOR from Moorella thermoacetica (MtPFOR). These snapshots include the native structure, those of lactyl-TPP and acetyl-TPP reaction intermediates, and the first of an OFOR with CoA bound. These structural data reveal the binding site of CoA as domain III, the function of which in OFORs was previously unknown, and establish sequence motifs for CoA binding in the OFOR superfamily. MtPFOR structures further show that domain III undergoes a conformational change upon CoA binding that seals off the active site and positions the thiolate of CoA directly adjacent to the TPP cofactor. These structural findings provide a molecular basis for the experimental observation that CoA binding accelerates catalysis by 10⁵-fold.

A ThDP-dependent enzymatic carboligation reaction involved in Neocarazostatin A tricyclic carbazole formation

Although the biosynthetic pathway of Neocarazostatin A (1) has been identified, the detailed enzymatic reactions underlying the assembly of the carbazole ring still remain largely unknown. We demonstrate here that NzsH, a putative thiamine diphosphate dependent enzyme, can catalyze an acyloin coupling reaction between indole-3-pyruvate and pyruvate to generate a β-ketoacid intermediate. Our findings thus shed light on further characterization of the unusual biosynthetic pathway of the bacterial tricyclic carbazole alkaloids.

Crosstalk between the tricarboxylic acid cycle and peptidoglycan synthesis in Caulobacter crescentus through the homeostatic control of α-ketoglutarate

To achieve robust replication, bacteria must integrate cellular metabolism and cell wall growth. While these two processes have been well characterized, the nature and extent of cross-regulation between them is not well understood. Here, using classical genetics, CRISPRi, metabolomics, transcriptomics and chemical complementation approaches, we show that a loss of the master regulator Hfq in Caulobacter crescentus alters central metabolism and results in cell shape defects in a nutrient-dependent manner. We demonstrate that the cell morphology phenotype in the hfq deletion mutant is attributable to a disruption of α-ketoglutarate (KG) homeostasis. In addition to serving as a key intermediate of the tricarboxylic acid (TCA) cycle, KG is a by-product of an enzymatic reaction required for the synthesis of peptidoglycan, a major component of the bacterial cell wall. Accumulation of KG in the hfq deletion mutant interferes with peptidoglycan synthesis, resulting in cell morphology defects and increased susceptibility to peptidoglycan-targeting antibiotics. This work thus reveals a direct crosstalk between the TCA cycle and cell wall morphogenesis. This crosstalk highlights the importance of metabolic homeostasis in not only ensuring adequate availability of biosynthetic precursors, but also in preventing interference with cellular processes in which these intermediates arise as by-products.

Streptomyces virginiae PPDC Is a New Type of Phenylpyruvate Decarboxylase Composed of Two Subunits

Streptomyces virginiae phenylpyruvate decarboxylase (PPDC) has not been identified before. Two putative branched-chain α-keto acid dehydrogenase subunit genes bkdC and bkdD from S. virginiae are similar to halves of other PPDC coding sequences. We cloned and characterized them biochemically in this work. The two proteins formed a stable complex attested by pull-down assay, consistent with the finding that their soluble expression was obtained only when they were coexpressed in Escherichia coli. The subunits were redesignated as SvPPDCα and SvPPDCβ, because the SvPPDCα/β complex catalyzed the conversion of phenylpyruvate to phenylacetaldehyde, reflecting the nature of the enzyme. Moreover, mutations of conserved residues in either of the two subunits led to inactivation or decreased specific activity of the enzymatic reaction. All previously identified PPDCs are encoded by a single gene. Here, we identified a new type of PPDC that contains two subunits, which gives new insights into the PPDC family.

Phosphonodifluoropyruvate is a mechanism-based inhibitor of phosphonopyruvate decarboxylase from Bacteroides fragilis

Bacteroides fragilis, a human pathogen, helps in the formation of intra-abdominal abscesses and is involved in 90% of anaerobic peritoneal infections. Phosphonopyruvate decarboxylase (PnPDC), a thiamin diphosphate (ThDP)-dependent enzyme, plays a key role in the formation of 2-aminoethylphosphonate, a component of the cell wall of B. fragilis. As such PnPDC is a possible target for therapeutic intervention in this, and other phosphonate producing organisms. However, the enzyme is of more general interest as it appears to be an evolutionary forerunner to the decarboxylase family of ThDP-dependent enzymes. To date, PnPDC has proved difficult to crystallize and no X-ray structures are available. In the past we have shown that ThDP-dependent enzymes will often crystallize if the cofactor has been irreversibly inactivated. To explore this possibility, and the utility of inhibitors of phosphonate biosynthesis as potential antibiotics, we synthesized phosphonodifluoropyruvate (PnDFP) as a prospective mechanism-based inhibitor of PnPDC. Here we provide evidence that PnDFP indeed inactivates the enzyme, that the inactivation is irreversible, and is accompanied by release of fluoride ion, i.e., PnDFP bears all the hallmarks of a mechanism-based inhibitor. Unfortunately, the enzyme remains refractive to crystallization.

Molecular evolution of acetohydroxyacid synthase in bacteria

Acetohydroxyacid synthase (AHAS) is the key enzyme in the biosynthetic pathways of branched chain amino acids in bacteria. Since it does not exist in animal and plant cells, AHAS is an attractive target for developing antimicrobials and herbicides. In some bacteria, there is a single copy of AHAS, while in others there are multiple copies. Therefore, it is necessary to investigate the origin and evolutionary pathway of various AHASs in bacteria. In this study, all the available protein sequences of AHAS in bacteria were investigated, and an evolutionary model of AHAS in bacteria is proposed, according to gene structure, organization and phylogeny. Multiple copies of AHAS in some bacteria might be evolved from the single copy of AHAS, the ancestor. Gene duplication, domain deletion and horizontal gene transfer might occur during the evolution of this enzyme. The results show the biological significance of AHAS, help to understand the functions of various AHASs in bacteria, and would be useful for developing industrial production strains of branched chain amino acids or novel antimicrobials.

Two strategies to engineer flexible loops for improved enzyme thermostability

Flexible sites are potential targets for engineering the stability of enzymes. Nevertheless, the success rate of the rigidifying flexible sites (RFS) strategy is still low due to a limited understanding of how to determine the best mutation candidates. In this study, two parallel strategies were applied to identify mutation candidates within the flexible loops of Escherichia coli transketolase (TK). The first was a “back to consensus mutations” approach, and the second was computational design based on ΔΔG calculations in Rosetta. Forty-nine single variants were generated and characterised experimentally. From these, three single-variants I189H, A282P, D143K were found to be more thermostable than wild-type TK. The combination of A282P with H192P, a variant constructed previously, resulted in the best all-round variant with a 3-fold improved half-life at 60 °C, 5-fold increased specific activity at 65 °C, 1.3-fold improved k_cat and a T_m increased by 5 °C above that of wild type. Based on a statistical analysis of the stability changes for all variants, the qualitative prediction accuracy of the Rosetta program reached 65.3%. Both of the two strategies investigated were useful in guiding mutation candidates to flexible loops, and had the potential to be used for other enzymes.

Lyme disease spirochaete Borrelia burgdorferi does not require thiamin

Thiamin pyrophosphate (ThDP), the active form of thiamin (vitamin B₁), is believed to be an essential cofactor for all living organisms^1,2. Here, we report the unprecedented result that thiamin is dispensable for the growth of the Lyme disease pathogen Borrelia burgdorferi (Bb)³. Bb lacks genes for thiamin biosynthesis and transport as well as known ThDP-dependent enzymes⁴, and we were unable to detect thiamin or its derivatives in Bb cells. We showed that eliminating thiamin in vitro and in vivo using BcmE, an enzyme that degrades thiamin, has no impact on Bb growth and survival during its enzootic infectious cycle. Finally, high-performance liquid chromatography analysis reveals that the level of thiamin and its derivatives in Ixodes scapularis ticks, the enzootic vector of Bb, is extremely low. These results suggest that by dispensing with use of thiamin, Borrelia, and perhaps other tick-transmitted bacterial pathogens, are uniquely adapted to survive in tick vectors before transmitting to mammalian hosts. To our knowledge, such a mechanism has not been reported previously in any living organisms.

Thiamin Diphosphate Activation in 1-Deoxy-d-xylulose 5-Phosphate Synthase: Insights into the Mechanism and Underlying Intermolecular Interactions

1-Deoxy-d-xylulose 5-phosphate synthase (DXS) is a thiamin diphosphate (TDP) dependent enzyme that marks the beginning of the methylerythritol 4-phosphate isoprenoid biosynthesis pathway. The mechanism of action for DXS is still poorly understood and begins with the formation of a thiazolium ylide. This TDP activation step is thought to proceed through an intramolecular deprotonation by the 4'-aminopyrimidine ring of TDP; however, this step would occur only after an initial deprotonation of its own 4'-amino group. The mechanism of the initial deprotonation has been hypothesized, by analogy to transketolases, to occur via a histidine or an active site water molecule. Results from hybrid quantum mechanical/molecular mechanical (QM/MM) reaction path calculations reveal an ∼10 kcal/mol difference in transition state energies, favoring a water mediated mechanism over direct deprotonation by histidine. This difference was determined to be largely governed by electrostatic changes induced by conformational variations in the active site. Additionally, mutagenesis studies reveal DXS to be an evolutionarily resilient enzyme. Particularly, we hypothesize that residues H82 and H304 may act in a compensatory fashion if the other is lost due to mutation. Further, nucleus-independent chemical shifts (NICSs) and aromatic stabilization energy (ASE) calculations suggest that reduction in TDP aromaticity also serves as a factor for regulating ylide formation and controlling reactivity.

Chopping and Changing: the Evolution of the Flavin-dependent Monooxygenases

Flavin-dependent monooxygenases play a variety of key physiological roles and are also very powerful biotechnological tools. These enzymes have been classified into eight different classes (A–H) based on their sequences and biochemical features. By combining structural and sequence analysis, and phylogenetic inference, we have explored the evolutionary history of classes A, B, E, F, and G and demonstrate that their multidomain architectures reflect their phylogenetic relationships, suggesting that the main evolutionary steps in their divergence are likely to have arisen from the recruitment of different domains. Additionally, the functional divergence within in each class appears to have been the result of other mechanisms such as a complex set of single-point mutations. Our results reinforce the idea that a main constraint on the evolution of cofactor-dependent enzymes is the functional binding of the cofactor. Additionally, a remarkable feature of this family is that the sequence of the key flavin adenine dinucleotide-binding domain is split into at least two parts in all classes studied here. We propose a complex set of evolutionary events that gave rise to the origin of the different classes within this family.

•

Changes in domain architectures reflect the phylogeny of flavin monooxygenases.

•

Recruitment of different domains has been the main force driving its evolution.

•

A notable feature of flavin monooxygenases is that the flavin adenine dinucleotide-binding domain is split.

•

Classes of monooxygenases emerged from an ancestral domain by structural changes.

Helper component-proteinase enhances the activity of 1-deoxy-D-xylulose-5-phosphate synthase and promotes the biosynthesis of plastidic isoprenoids in Potato virus Y-infected tobacco

Virus-infected plants show strong morphological and physiological alterations. Many physiological processes in chloroplast are affected, including the plastidic isoprenoid biosynthetic pathway [the 2C-methyl-D-erythritol-4-phosphate (MEP) pathway]; indeed, isoprenoid contents have been demonstrated to be altered in virus-infected plants. In this study, we found that the levels of photosynthetic pigments and abscisic acid (ABA) were altered in Potato virus Y (PVY)-infected tobacco. Using yeast two-hybrid assays, we demonstrated an interaction between virus protein PVY helper component-proteinase (HC-Pro) and tobacco chloroplast protein 1-deoxy-D-xylulose-5-phosphate synthase (NtDXS). This interaction was confirmed using bimolecular fluorescence complementation (BiFC) assays and pull-down assays. The Transket_pyr domain (residues 394-561) of NtDXS was required for interaction with HC-Pro, while the N-terminal region of HC-Pro (residues 1-97) was necessary for interaction with NtDXS. Using in vitro enzyme activity assays, PVY HC-Pro was found to promote the synthase activity of NtDXS. We observed increases in photosynthetic pigment contents and ABA levels in transgenic plants with HC-Pro accumulating in the chloroplasts. During virus infection, the enhancement of plastidic isoprenoid biosynthesis was attributed to the enhancement of DXS activity by HC-Pro. Our study reveals a new role of HC-Pro in the host plant metabolic system and will contribute to the study of host-virus relationships.

The Structure of an Oxalate Oxidoreductase Provides Insight into Microbial 2-Oxoacid Metabolism

Thiamine pyrophosphate (TPP), a derivative of vitamin B1, is a versatile and ubiquitous cofactor. When coupled with [4Fe-4S] clusters in microbial 2-oxoacid:ferredoxin oxidoreductases (OFORs), TPP is involved in catalyzing low-potential redox reactions that are important for the synthesis of key metabolites and the reduction of N2, H(+), and CO2. We have determined the high-resolution (2.27 Å) crystal structure of the TPP-dependent oxalate oxidoreductase (OOR), an enzyme that allows microbes to grow on oxalate, a widely occurring dicarboxylic acid that is found in soil and freshwater and is responsible for kidney stone disease in humans. OOR catalyzes the anaerobic oxidation of oxalate, harvesting the low-potential electrons for use in anaerobic reduction and fixation of CO2. We compare the OOR structure to that of the only other structurally characterized OFOR family member, pyruvate:ferredoxin oxidoreductase. This side-by-side structural analysis highlights the key similarities and differences that are relevant for the chemistry of this entire class of TPP-utilizing enzymes.

A metabolic strategy to enhance long-term survival by Phx1 through stationary phase-specific pyruvate decarboxylases in fission yeast

In the fission yeast Schizosaccharomyces pombe, the stationary phase-specific transcription factor Phx1 contributes to long-term survival, stress tolerance, and meiosis. We identified Phx1-dependent genes through transcriptome analysis, and further analyzed those related with carbohydrate and thiamine metabolism, whose expression decreased in Δphx1. Consistent with mRNA changes, the level of thiamine pyrophosphate (TPP) and TPP-utilizing pyruvate decarboxylase activity that converts pyruvate to acetaldehyde were also reduced in the mutant. Therefore, Phx1 appears to shift metabolic flux by diverting pyruvate from the TCA cycle and respiration to ethanol fermentation. Among the four predicted genes for pyruvate decarboxylase, only the Phx1-dependent genes (pdc201⁺ and pdc202⁺) contributed to long-term survival as judged by mutation and overexpression studies. These findings indicate that the Phx1-mediated long-term survival is achieved primarily through increasing the synthesis and activity of pyruvate decarboxylase. Consistent with this hypothesis, we observed that Phx1 curtailed respiration when cells entered stationary phase. Introduction of Δphx1 mutation compromised the long-lived phenotypes of Δpka1 and Δsck2 mutants that are devoid of pro-aging kinases of nutrient-signalling pathways, and of the Δpyp1 mutant with constitutively activated stress-responsive kinase Sty1. Therefore, achievement of long-term viability through both nutrient limitation and anti-stress response appears to be dependent on Phx1.

Perturbation of the monomer-monomer interfaces of the benzoylformate decarboxylase tetramer

The X-ray structure of benzoylformate decarboxylase (BFDC) from Pseudomonas putida ATCC 12633 shows it to be a tetramer. This was believed to be typical of all thiamin diphosphate-dependent decarboxylases until recently when the structure of KdcA, a branched-chain 2-keto acid decarboxylase from Lactococcus lactis, showed it to be a homodimer. This lent credence to earlier unfolding experiments on pyruvate decarboxylase from Saccharomyces cerevisiae that indicated that it might be active as a dimer. To investigate this possibility in BFDC, we sought to shift the equilibrium toward dimer formation. Point mutations were made in the noncatalytic monomer–monomer interfaces, but these had a minimal effect on both tetramer formation and catalytic activity. Subsequently, the R141E/Y288A/A306F variant was shown by analytical ultracentrifugation to be partially dimeric. It was also found to be catalytically inactive. Further experiments revealed that just two mutations, R141E and A306F, were sufficient to markedly alter the dimer–tetramer equilibrium and to provide an ∼450-fold decrease in k_cat. Equilibrium denaturation studies suggested that the residual activity was possibly due to the presence of residual tetramer. The structures of the R141E and A306F variants, determined to <1.5 Å resolution, hinted that disruption of the monomer interfaces will be accompanied by movement of a loop containing Leu109 and Leu110. As these residues contribute to the hydrophobicity of the active site and the correct positioning of the substrate, it seems that tetramer formation may well be critical to the catalytic activity of BFDC.

Structure and functioning mechanism of transketolase

Studies of thiamine diphosphate-dependent enzymes appear to have commenced in 1937, with the isolation of the coenzyme of yeast pyruvate decarboxylase, which was demonstrated to be a diphosphoric ester of thiamine. For quite a long time, these studies were largely focused on enzymes decarboxylating α-keto acids, such as pyruvate decarboxylase and pyruvate dehydrogenase complexes. Transketolase, discovered independently by Racker and Horecker in 1953 (and named by Racker) [1], did not receive much attention until 1992, when crystal X-ray structure analysis of the enzyme from Saccharomyces cerevisiae was performed [2]. These data, together with the results of site-directed mutagenesis, made it possible to understand in detail the mechanism of thiamine diphosphate-dependent catalysis. Some progress was also made in studies of the functional properties of transketolase. The last review on transketolase, which was fairly complete, appeared in 1998 [3]. Therefore, the publication of this paper should not seem premature.

Decarboxylation of pyruvate to acetaldehyde for ethanol production by hyperthermophiles

Pyruvate decarboxylase (PDC encoded by pdc) is a thiamine pyrophosphate (TPP)-containing enzyme responsible for the conversion of pyruvate to acetaldehyde in many mesophilic organisms. However, no pdc/PDC homolog has yet been found in fully sequenced genomes and proteomes of hyper/thermophiles. The only PDC activity reported in hyperthermophiles was a bifunctional, TPP- and CoA-dependent pyruvate ferredoxin oxidoreductase (POR)/PDC enzyme from the hyperthermophilic archaeon Pyrococcus furiosus. Another enzyme known to be involved in catalysis of acetaldehyde production from pyruvate is CoA-acetylating acetaldehyde dehydrogenase (AcDH encoded by mhpF and adhE). Pyruvate is oxidized into acetyl-CoA by either POR or pyruvate formate lyase (PFL), and AcDH catalyzes the reduction of acetyl-CoA to acetaldehyde in mesophilic organisms. AcDH is present in some mesophilic (such as clostridia) and thermophilic bacteria (e.g., Geobacillus and Thermoanaerobacter). However, no AcDH gene or protein homologs could be found in the released genomes and proteomes of hyperthermophiles. Moreover, no such activity was detectable from the cell-free extracts of different hyperthermophiles under different assay conditions. In conclusion, no commonly-known PDCs was found in hyperthermophiles. Instead of the commonly-known PDC, it appears that at least one multifunctional enzyme is responsible for catalyzing the non-oxidative decarboxylation of pyruvate to acetaldehyde in hyperthermophiles.

The effect of MEP pathway and other inhibitors on the intracellular localization of a plasma membrane-targeted, isoprenylable GFP reporter protein in tobacco BY-2 cells

We have established anin vivovisualization system for the geranylgeranylation of proteins in a stably transformed tobacco BY-2 cell line, based on the expression of a dexamethasone-inducible GFP fused to the carboxy-terminal basic domain of the rice calmodulin CaM61, which naturally bears a CaaL geranylgeranylation motif (GFP-BD-CVIL). By using pathway-specific inhibitors it was demonstrated that inhibition of the methylerythritol phosphate (MEP) pathway with known inhibitors like oxoclomazone and fosmidomycin, as well as inhibition of the protein geranylgeranyltransferase type 1 (PGGT-1), shifted the localization of the GFP-BD-CVIL protein from the membrane to the nucleus. In contrast, the inhibition of the mevalonate (MVA) pathway with mevinolin did not affect the localization. During the present work, this test system has been used to examine the effect of newly designed inhibitors of the MEP pathway and inhibitors of sterol biosynthesis such as squalestatin, terbinafine and Ro48-8071. In addition, we also studied the impact of different post-prenylation inhibitors or those suspected to affect the transport of proteins to the plasma membrane on the localization of the geranylgeranylable fusion protein GFP-BD-CVIL.

The effect of MEP pathway and other inhibitors on the intracellular localization of a plasma membrane-targeted, isoprenylable GFP reporter protein in tobacco BY-2 cells

We have established an in vivo visualization system for the geranylgeranylation of proteins in a stably transformed tobacco BY-2 cell line, based on the expression of a dexamethasone-inducible GFP fused to the carboxy-terminal basic domain of the rice calmodulin CaM61, which naturally bears a CaaL geranylgeranylation motif (GFP-BD-CVIL). By using pathway-specific inhibitors it was demonstrated that inhibition of the methylerythritol phosphate (MEP) pathway with known inhibitors like oxoclomazone and fosmidomycin, as well as inhibition of the protein geranylgeranyltransferase type 1 (PGGT-1), shifted the localization of the GFP-BD-CVIL protein from the membrane to the nucleus. In contrast, the inhibition of the mevalonate (MVA) pathway with mevinolin did not affect the localization. During the present work, this test system has been used to examine the effect of newly designed inhibitors of the MEP pathway and inhibitors of sterol biosynthesis such as squalestatin, terbinafine and Ro48-8071. In addition, we also studied the impact of different post-prenylation inhibitors or those suspected to affect the transport of proteins to the plasma membrane on the localization of the geranylgeranylable fusion protein GFP-BD-CVIL.

Structure and function of the transketolase from Mycobacterium tuberculosis and comparison with the human enzyme

The transketolase (TKT) enzyme in Mycobacterium tuberculosis represents a novel drug target for tuberculosis treatment and has low homology with the orthologous human enzyme. Here, we report on the structural and kinetic characterization of the transketolase from M. tuberculosis (TBTKT), a homodimer whose monomers each comprise 700 amino acids. We show that TBTKT catalyses the oxidation of donor sugars xylulose-5-phosphate and fructose-6-phosphate as well as the reduction of the acceptor sugar ribose-5-phosphate. An invariant residue of the TKT consensus sequence required for thiamine cofactor binding is mutated in TBTKT; yet its catalytic activities are unaffected, and the 2.5 Å resolution structure of full-length TBTKT provides an explanation for this. Key structural differences between the human and mycobacterial TKT enzymes that impact both substrate and cofactor recognition and binding were uncovered. These changes explain the kinetic differences between TBTKT and its human counterpart, and their differential inhibition by small molecules. The availability of a detailed structural model of TBTKT will enable differences between human and M. tuberculosis TKT structures to be exploited to design selective inhibitors with potential antitubercular activity.

Probing a coral genome for components of the photoprotective scytonemin biosynthetic pathway and the 2-aminoethylphosphonate pathway

Genome sequences of the reef-building coral, Acropora digitifera, have been decoded. Acropora inhabits an environment with intense ultraviolet exposure and hosts the photosynthetic endosymbiont, Symbiodinium. Acropora homologs of all four genes necessary for biosynthesis of the photoprotective cyanobacterial compound, shinorine, are present. Among metazoans, these genes are found only in anthozoans. To gain further evolutionary insights into biosynthesis of photoprotective compounds and associated coral proteins, we surveyed the Acropora genome for 18 clustered genes involved in cyanobacterial synthesis of the anti-UV compound, scytonemin, even though it had not previously been detected in corals. We identified candidates for only 6 of the 18 genes, including tyrP, scyA, and scyB. Therefore, it does not appear that Acropora digitifera can synthesize scytonemin independently. On the other hand, molecular phylogenetic analysis showed that one tyrosinase gene is an ortholog of vertebrate tyrosinase genes and that the coral homologs, scyA and scyB, are similar to bacterial metabolic genes, phosphonopyruvate (ppyr) decarboxylase and glutamate dehydrogenase (GDH), respectively. Further genomic searches for ppyr gene-related biosynthetic components indicate that the coral possesses a metabolic pathway similar to the bacterial 2-aminoethylphosphonate (AEP) biosynthetic pathway. The results suggest that de novo synthesis of carbon-phosphorus compounds is performed in corals.

A standard numbering scheme for thiamine diphosphate-dependent decarboxylases

Background

Standard numbering schemes for families of homologous proteins allow for the unambiguous identification of functionally and structurally relevant residues, to communicate results on mutations, and to systematically analyse sequence-function relationships in protein families. Standard numbering schemes have been successfully implemented for several protein families, including lactamases and antibodies, whereas a numbering scheme for the structural family of thiamine-diphosphate (ThDP) -dependent decarboxylases, a large subfamily of the class of ThDP-dependent enzymes encompassing pyruvate-, benzoylformate-, 2-oxo acid-, indolpyruvate- and phenylpyruvate decarboxylases, benzaldehyde lyase, acetohydroxyacid synthases and 2-succinyl-5-enolpyruvyl-6-hydroxy-3-cyclohexadiene-1-carboxylate synthase (MenD) is still missing.

Despite a high structural similarity between the members of the ThDP-dependent decarboxylases, their sequences are diverse and make a pairwise sequence comparison of protein family members difficult.

Results

We developed and validated a standard numbering scheme for the family of ThDP-dependent decarboxylases. A profile hidden Markov model (HMM) was created using a set of representative sequences from the family of ThDP-dependent decarboxylases. The pyruvate decarboxylase from S. cerevisiae (PDB: 2VK8) was chosen as a reference because it is a well characterized enzyme. The crystal structure with the PDB identifier 2VK8 encompasses the structure of the ScPDC mutant E477Q, the cofactors ThDP and Mg²⁺ as well as the substrate analogue (2S)-2-hydroxypropanoic acid. The absolute numbering of this reference sequence was transferred to all members of the ThDP-dependent decarboxylase protein family. Subsequently, the numbering scheme was integrated into the already established Thiamine-diphosphate dependent Enzyme Engineering Database (TEED) and was used to systematically analyze functionally and structurally relevant positions in the superfamily of ThDP-dependent decarboxylases.

Conclusions

The numbering scheme serves as a tool for the reliable sequence alignment of ThDP-dependent decarboxylases and the unambiguous identification and communication of corresponding positions. Thus, it is the basis for the systematic and automated analysis of sequence-encoded properties such as structural and functional relevance of amino acid positions, because the analysis of conserved positions, the identification of correlated mutations and the determination of subfamily specific amino acid distributions depend on reliable multisequence alignments and the unambiguous identification of the alignment columns. The method is reliable and robust and can easily be adapted to further protein families.

Directed evolution to re-adapt a co-evolved network within an enzyme

We have previously used targeted active-site saturation mutagenesis to identify a number of transketolase single mutants that improved activity towards either glycolaldehyde (GA), or the non-natural substrate propionaldehyde (PA). Here, all attempts to recombine the singles into double mutants led to unexpected losses of specific activity towards both substrates. A typical trade-off occurred between soluble expression levels and specific activity for all single mutants, but many double mutants decreased both properties more severely suggesting a critical loss of protein stability or native folding. Statistical coupling analysis (SCA) of a large multiple sequence alignment revealed a network of nine co-evolved residues that affected all but one double mutant. Such networks maintain important functional properties such as activity, specificity, folding, stability, and solubility and may be rapidly disrupted by introducing one or more non-naturally occurring mutations. To identify variants of this network that would accept and improve upon our best D469 mutants for activity towards PA, we created a library of random single, double and triple mutants across seven of the co-evolved residues, combining our D469 variants with only naturally occurring mutations at the remaining sites. A triple mutant cluster at D469, E498 and R520 was found to behave synergistically for the specific activity towards PA. Protein expression was severely reduced by E498D and improved by R520Q, yet variants containing both mutations led to improved specific activity and enzyme expression, but with loss of solubility and the formation of inclusion bodies. D469S and R520Q combined synergistically to improve k_cat 20-fold for PA, more than for any previous transketolase mutant. R520Q also doubled the specific activity of the previously identified D469T to create our most active transketolase mutant to date. Our results show that recombining active-site mutants obtained by saturation mutagenesis can rapidly destabilise critical networks of co-evolved residues, whereas beneficial single mutants can be retained and improved upon by randomly recombining them with natural variants at other positions in the network.

Expressed sequence tag analysis of khat (Catha edulis) provides a putative molecular biochemical basis for the biosynthesis of phenylpropylamino alkaloids

Khat (Catha edulis Forsk.) is a flowering perennial shrub cultivated for its neurostimulant properties resulting mainly from the occurrence of (S)-cathinone in young leaves. The biosynthesis of (S)-cathinone and the related phenylpropylamino alkaloids (1S,2S)-cathine and (1R,2S)-norephedrine is not well characterized in plants. We prepared a cDNA library from young khat leaves and sequenced 4,896 random clones, generating an expressed sequence tag (EST) library of 3,293 unigenes. Putative functions were assigned to > 98% of the ESTs, providing a key resource for gene discovery. Candidates potentially involved at various stages of phenylpropylamino alkaloid biosynthesis from L-phenylalanine to (1S,2S)-cathine were identified.

Structural stability of E. coli transketolase to temperature and pH denaturation

We have previously shown that the denaturation of TK with urea follows a non-aggregating though irreversible denaturation pathway in which the cofactor binding appears to become altered but without dissociating, then followed at higher urea by partial denaturation of the homodimer prior to any further unfolding or dissociation of the two monomers. Urea is not typically present during biocatalysis, whereas access to TK enzymes that retain activity at increased temperature and extreme pH would be useful for operation under conditions that increase substrate and product stability or solubility. To provide further insight into the underlying causes of its deactivation in process conditions, we have characterised the effects of temperature and pH on the structure, stability, aggregation and activity of Escherichia coli transketolase. The activity of TK was initially found to progressively improve after pre-incubation at increasing temperatures. Loss of activity at higher temperature and low pH resulted primarily from protein denaturation and subsequent irreversible aggregation. By contrast, high pH resulted in the formation of a native-like state that was only partially inactive. The apo-TK enzyme structure content also increased at pH 9 to converge on that of the holo-TK. While cofactor dissociation was previously proposed for high pH deactivation, the observed structural changes in apo-TK but not holo-TK indicate a more complex mechanism.

Identification and characterization of oxalate oxidoreductase, a novel thiamine pyrophosphate-dependent 2-oxoacid oxidoreductase that enables anaerobic growth on oxalate

Moorella thermoacetica is an anaerobic acetogen, a class of bacteria that is found in the soil, the animal gastrointestinal tract, and the rumen. This organism engages the Wood-Ljungdahl pathway of anaerobic CO(2) fixation for heterotrophic or autotrophic growth. This paper describes a novel enzyme, oxalate oxidoreductase (OOR), that enables M. thermoacetica to grow on oxalate, which is produced in soil and is a common component of kidney stones. Exposure to oxalate leads to the induction of three proteins that are subunits of OOR, which oxidizes oxalate coupled to the production of two electrons and CO(2) or bicarbonate. Like other members of the 2-oxoacid:ferredoxin oxidoreductase family, OOR contains thiamine pyrophosphate and three [Fe(4)S(4)] clusters. However, unlike previously characterized members of this family, OOR does not use coenzyme A as a substrate. Oxalate is oxidized with a k(cat) of 0.09 s(-1) and a K(m) of 58 μM at pH 8. OOR also oxidizes a few other 2-oxoacids (which do not induce OOR) also without any requirement for CoA. The enzyme transfers its reducing equivalents to a broad range of electron acceptors, including ferredoxin and the nickel-dependent carbon monoxide dehydrogenase. In conjunction with the well characterized Wood-Ljungdahl pathway, OOR should be sufficient for oxalate metabolism by M. thermoacetica, and it constitutes a novel pathway for oxalate metabolism.

The Thiamine diphosphate dependent Enzyme Engineering Database: a tool for the systematic analysis of sequence and structure relations

BACKGROUND

Thiamine diphosphate (ThDP)-dependent enzymes form a vast and diverse class of proteins, catalyzing a wide variety of enzymatic reactions including the formation or cleavage of carbon-sulfur, carbon-oxygen, carbon-nitrogen, and especially carbon-carbon bonds. Although very diverse in sequence and domain organisation, they share two common protein domains, the pyrophosphate (PP) and the pyrimidine (PYR) domain. For the comprehensive and systematic comparison of protein sequences and structures the Thiamine diphosphate (ThDP)-dependent Enzyme Engineering Database (TEED) was established.

DESCRIPTION

The TEED http://www.teed.uni-stuttgart.de contains 12048 sequence entries which were assigned to 9443 different proteins and 379 structure entries. Proteins were assigned to 8 different superfamilies and 63 homologous protein families. For each family, the TEED offers multisequence alignments, phylogenetic trees, and family-specific HMM profiles. The conserved pyrophosphate (PP) and pyrimidine (PYR) domains have been annotated, which allows the analysis of sequence similarities for a broad variety of proteins. Human ThDP-dependent enzymes are known to be involved in many diseases. 20 different proteins and over 40 single nucleotide polymorphisms (SNPs) of human ThDP-dependent enzymes were identified in the TEED.

CONCLUSIONS

The online accessible version of the TEED has been designed to serve as a navigation and analysis tool for the large and diverse family of ThDP-dependent enzymes.