New insights into structure-function relationships of oxalyl CoA decarboxylase from Escherichia coli

A conserved oxalyl-coenzyme A decarboxylase in oxalate catabolism

The ability to biosynthesize oxalic acid can provide beneficial functions to plants; however, uncontrolled or prolonged exposure to this strong organic acid results in multiple physiological problems. Such problems include a disruption of membrane integrity, mitochondrial function, metal chelation, and free radical formation. Recent work suggests that a CoA-dependent pathway of oxalate catabolism plays a critical role in regulating tissue oxalate concentrations in plants. Although this CoA-dependent pathway of oxalate catabolism is important, large gaps in our knowledge of the enzymes catalyzing each step remain. Evidence that an oxalyl-CoA decarboxylase (OXC) catalyzes the second step in this pathway, accelerating the conversion of oxalyl-CoA to formyl-CoA, has been reported. Induction studies revealed that OXC gene expression was upregulated in response to an exogenous oxalate supply. Phylogenetic analysis indicates that OXCs are conserved across plant species. Evolutionarily the plant OXCs can be separated into dicot and monocot classes. Multiple sequence alignments and molecular modeling suggest that OXCs have similar functionality with three conserved domains, the N-terminal PYR domain, the middle R domain, and the C-terminal PP domain. Further study of this CoA-dependent pathway of oxalate degradation would benefit efforts to develop new strategies to improve the nutrition quality of crops.

Probiotics in the Prevention of the Calcium Oxalate Urolithiasis

Nephrolithiasis ranks third among urological diseases in terms of prevalence, making up about 15% of cases. The continued increase in the incidence of nephrolithiasis is most probably due to changes in eating habits (high protein, sodium, and sugar diets) and lifestyle (reduced physical activity) in all developed countries. Some 80% of all kidney stones cases are oxalate urolithiasis, which is also characterized by the highest risk of recurrence. Frequent relapses of nephrolithiasis contribute to severe complications and high treatment costs. Unfortunately, there is no known effective way to prevent urolithiasis at present. In cases of diet-related urolithiasis, dietary changes may prevent recurrence. However, in some patients, the condition is unrelated to diet; in such cases, there is evidence to support the use of stone-related medications. Interestingly, a growing body of evidence indicates the potential of the microbiome to reduce the risk of developing renal colic. Previous studies have primarily focused on the use of Oxalobacterformigenes in patients with urolithiasis. Unfortunately, this bacterium is not an ideal probiotic due to its antibiotic sensitivity and low pH. Therefore, subsequent studies sought to find bacteria which are capable of oxalate degradation, focusing on well-known probiotics including Lactobacillus and Bifidobacterium strains, Eubacterium lentum, Enterococcus faecalis, and Escherichia coli.

The identification and characterization of an oxalyl-CoA synthetase from grass pea (Lathyrus sativus L.)

Oxalic acid is a small metabolite that can be found in many plants in which it serves as protection from herbivores, a chelator of metal ions, a regulator of calcium...

Mechanistic details of the actinobacterial lyase-catalyzed degradation reaction of 2-hydroxyisobutyryl-CoA

Actinobacterial 2-hydroxyacyl-CoA lyase reversibly catalyzes the thiamine diphosphate-dependent cleavage of 2-hydroxyisobutyryl-CoA to formyl-CoA and acetone. This enzyme has great potential for use in synthetic one-carbon assimilation pathways for sustainable production of chemicals, but lacks details of substrate binding and reaction mechanism for biochemical reengineering. We determined crystal structures of the tetrameric enzyme in the closed conformation with bound substrate, covalent postcleavage intermediate, and products, shedding light on active site architecture and substrate interactions. Together with molecular dynamics simulations of the covalent precleavage complex, the complete catalytic cycle is structurally portrayed, revealing a proton transfer from the substrate acyl Cβ hydroxyl to residue E493 that returns it subsequently to the postcleavage Cα-carbanion intermediate. Kinetic parameters obtained for mutants E493A, E493Q, and E493K confirm the catalytic role of E493 in the WT enzyme. However, the 10- and 50-fold reduction in lyase activity in the E493A and E493Q mutants, respectively, compared with WT suggests that water molecules may contribute to proton transfer. The putative catalytic glutamate is located on a short α-helix close to the active site. This structural feature appears to be conserved in related lyases, such as human 2-hydroxyacyl-CoA lyase 2. Interestingly, a unique feature of the actinobacterial 2-hydroxyacyl-CoA lyase is a large C-terminal lid domain that, together with active site residues L127 and I492, restricts substrate size to ≤C5 2-hydroxyacyl residues. These details about the catalytic mechanism and determinants of substrate specificity pave the ground for designing tailored catalysts for acyloin condensations for one-carbon and short-chain substrates in biotechnological applications.

Engineering a Highly Efficient Carboligase for Synthetic One-Carbon Metabolism

One of the biggest challenges to realize a circular carbon economy is the synthesis of complex carbon compounds from one-carbon (C1) building blocks. Since the natural solution space of C1–C1 condensations is limited to highly complex enzymes, the development of more simple and robust biocatalysts may facilitate the engineering of C1 assimilation routes. Thiamine diphosphate-dependent enzymes harbor great potential for this task, due to their ability to create C–C bonds. Here, we employed structure-guided iterative saturation mutagenesis to convert oxalyl-CoA decarboxylase (OXC) from Methylobacterium extorquens into a glycolyl-CoA synthase (GCS) that allows for the direct condensation of the two C1 units formyl-CoA and formaldehyde. A quadruple variant MeOXC4 showed a 100 000-fold switch between OXC and GCS activities, a 200-fold increase in the GCS activity compared to the wild type, and formaldehyde affinity that is comparable to natural formaldehyde-converting enzymes. Notably, MeOCX4 outcompetes all other natural and engineered enzymes for C1–C1 condensations by more than 40-fold in catalytic efficiency and is highly soluble in Escherichia coli. In addition to the increased GCS activity, MeOXC4 showed up to 300-fold higher activity than the wild type toward a broad range of carbonyl acceptor substrates. When applied in vivo, MeOXC4 enables the production of glycolate from formaldehyde, overcoming the current bottleneck of C1–C1 condensation in whole-cell bioconversions and paving the way toward synthetic C1 assimilation routes in vivo.

Allosteric regulation of menaquinone (vitamin K2) biosynthesis in the human pathogen Mycobacterium tuberculosis

Menaquinone (vitamin K₂) plays a vital role in energy generation and environmental adaptation in many bacteria, including the human pathogen Mycobacterium tuberculosis (Mtb). Although menaquinone levels are known to be tightly linked to the cellular redox/energy status of the cell, the regulatory mechanisms underpinning this phenomenon are unclear. The first committed step in menaquinone biosynthesis is catalyzed by MenD, a thiamine diphosphate-dependent enzyme comprising three domains. Domains I and III form the MenD active site, but no function has yet been ascribed to domain II. Here, we show that the last cytosolic metabolite in the menaquinone biosynthesis pathway, 1,4-dihydroxy-2-naphthoic acid (DHNA), binds to domain II of Mtb-MenD and inhibits its activity. Using X-ray crystallography of four apo- and cofactor-bound Mtb-MenD structures, along with several spectroscopy assays, we identified three arginine residues (Arg-97, Arg-277, and Arg-303) that are important for both enzyme activity and the feedback inhibition by DHNA. Among these residues, Arg-277 appeared to be particularly important for signal propagation from the allosteric site to the active site. This is the first evidence of feedback regulation of the menaquinone biosynthesis pathway in bacteria, identifying a protein-level regulatory mechanism that controls menaquinone levels within the cell and may therefore represent a good target for disrupting menaquinone biosynthesis in M. tuberculosis.

Machine learning and structural analysis of Mycobacterium tuberculosis pan-genome identifies genetic signatures of antibiotic resistance

Mycobacterium tuberculosis is a serious human pathogen threat exhibiting complex evolution of antimicrobial resistance (AMR). Accordingly, the many publicly available datasets describing its AMR characteristics demand disparate data-type analyses. Here, we develop a reference strain-agnostic computational platform that uses machine learning approaches, complemented by both genetic interaction analysis and 3D structural mutation-mapping, to identify signatures of AMR evolution to 13 antibiotics. This platform is applied to 1595 sequenced strains to yield four key results. First, a pan-genome analysis shows that M. tuberculosis is highly conserved with sequenced variation concentrated in PE/PPE/PGRS genes. Second, the platform corroborates 33 genes known to confer resistance and identifies 24 new genetic signatures of AMR. Third, 97 epistatic interactions across 10 resistance classes are revealed. Fourth, detailed structural analysis of these genes yields mechanistic bases for their selection. The platform can be used to study other human pathogens.

Modeling time-series data from microbial communities

As sequencing technologies have advanced, the amount of information regarding the composition of bacterial communities from various environments (for example, skin or soil) has grown exponentially. To date, most work has focused on cataloging taxa present in samples and determining whether the distribution of taxa shifts with exogenous covariates. However, important questions regarding how taxa interact with each other and their environment remain open thus preventing in-depth ecological understanding of microbiomes. Time-series data from 16S rDNA amplicon sequencing are becoming more common within microbial ecology, but methods to infer ecological interactions from these longitudinal data are limited. We address this gap by presenting a method of analysis using Poisson regression fit with an elastic-net penalty that (1) takes advantage of the fact that the data are time series; (2) constrains estimates to allow for the possibility of many more interactions than data; and (3) is scalable enough to handle data consisting of thousands of taxa. We test the method on gut microbiome data from white-throated woodrats (Neotoma albigula) that were fed varying amounts of the plant secondary compound oxalate over a period of 22 days to estimate interactions between OTUs and their environment.

A functional analysis of the formyl-coenzyme A (frc) gene from Lactobacillus reuteri 100-23C

AIM

To examine the role of the Lactobacillus reuteri 100-23C frc gene product in oxalate metabolism, host colonization and the acid stress response.

METHODS AND RESULTS

Genes encoding putative formyl-CoA transferase (frc) and oxalyl-CoA decarboxylase (oxc) enzymes are present in the genome sequences of Lact. reuteri strains. Two strains isolated from humans harboured an IS200 insertion sequence in the frc ORF and a group 2 intron-associated transposase downstream of the frc gene, both of which were lacking in two strains of animal origin, which contained intact frc and oxc genes. An frc(-) insertional mutant of Lact. reuteri 100-23C was compared with the parent strain with respect to oxalate degradation, colonization of an RLF-mouse host model and growth in the presence of acids. Neither parent nor mutant degraded oxalate in vitro or in vivo. However, the parent outcompeted the frc(-) mutant in the mouse intestine during co-colonization and the frc(-) mutant showed a reduced growth rate in the presence of hydrochloric acid.

CONCLUSIONS

Intact oxc and frc genes do not ensure oxalate degradation under the conditions tested. The frc gene product is important during host colonization and survival of acid stress by Lact. reuteri 100-23C.

SIGNIFICANCE AND IMPACT OF THE STUDY

Oxalate metabolism by oxalate-degrading intestinal bacterial strains may be important in preventing urolithiasis and might lead to the derivation of probiotic products. To produce safe and efficacious probiotics, however, an understanding of the genetic characteristics of potential oxalate degraders must be obtained, together with knowledge of their functional ramifications.

YfdW and YfdU are required for oxalate-induced acid tolerance in Escherichia coli K-12

Escherichia coli has several mechanisms for surviving low-pH stress. We report that oxalic acid, a small-chain organic acid (SCOA), induces a moderate acid tolerance response (ATR) in two ways. Adaptation of E. coli K-12 at pH 5.5 with 50 mM oxalate and inclusion of 25 mM oxalate in pH 3.0 minimal challenge medium separately conferred protection, with 67% ± 7% and 87% ± 17% survival after 2 h, respectively. The combination of oxalate adaptation and oxalate supplementation in the challenge medium resulted in increased survival over adaptation or oxalate in the challenge medium alone. The enzymes YfdW, a formyl coenzyme A (CoA) transferase, and YfdU, an oxalyl-CoA decarboxylase, are required for the adaptation effect but not during challenge. Unlike other SCOAs, this oxalate ATR is not a part of the RpoS regulon but appears to be linked to the signal protein GadE. We theorize that this oxalate ATR could enhance the pathogenesis of virulent E. coli consumed with oxalate-containing foods like spinach.

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.

In silico identification of epitopes in Mycobacterium avium subsp. paratuberculosis proteins that were upregulated under stress conditions

Johne's disease in ruminants is caused by Mycobacterium avium subsp. paratuberculosis. Diagnosis of M. avium subsp. paratuberculosis infection is difficult, especially in the early stages. To date, ideal antigen candidates are not available for efficient immunization or immunodiagnosis. This study reports the in silico selection and subsequent analysis of epitopes of M. avium subsp. paratuberculosis proteins that were found to be upregulated under stress conditions as a means to identify immunogenic candidate proteins. Previous studies have reported differential regulation of proteins when M. avium subsp. paratuberculosis is exposed to stressors which induce a response similar to dormancy. Dormancy may be involved in evading host defense mechanisms, and the host may also mount an immune response against these proteins. Twenty-five M. avium subsp. paratuberculosis proteins that were previously identified as being upregulated under in vitro stress conditions were analyzed for B and T cell epitopes by use of the prediction tools at the Immune Epitope Database and Analysis Resource. Major histocompatibility complex class I T cell epitopes were predicted using an artificial neural network method, and class II T cell epitopes were predicted using the consensus method. Conformational B cell epitopes were predicted from the relevant three-dimensional structure template for each protein. Based on the greatest number of predicted epitopes, eight proteins (MAP2698c [encoded by desA2], MAP2312c [encoded by fadE19], MAP3651c [encoded by fadE3_2], MAP2872c [encoded by fabG5_2], MAP3523c [encoded by oxcA], MAP0187c [encoded by sodA], and the hypothetical proteins MAP3567 and MAP1168c) were identified as potential candidates for study of antibody- and cell-mediated immune responses within infected hosts.

Substrate specificity in thiamin diphosphate-dependent decarboxylases

Thiamin diphosphate (ThDP) is the biologically active form of vitamin B(1), and ThDP-dependent enzymes are found in all forms of life. The catalytic mechanism of this family requires the formation of a common intermediate, the 2α-carbanion-enamine, regardless of whether the enzyme is involved in C-C bond formation or breakdown, or even formation of C-N, C-O and C-S bonds. This demands that the enzymes must screen substrates prior to, and/or after, formation of the common intermediate. This review is focused on the group for which the second step is the protonation of the 2α-carbanion, i.e., the ThDP-dependent decarboxylases. Based on kinetic data, sequence/structure alignments and mutagenesis studies the factors involved in substrate specificity have been identified.