Kinetic and mechanistic analysis of the Escherichia coli ribD-encoded bifunctional deaminase-reductase involved in riboflavin biosynthesis

Improving pathway prediction accuracy of constraints-based metabolic network models by treating enzymes as microcompartments

Metabolic network models have become increasingly precise and accurate as the most widespread and practical digital representations of living cells. The prediction functions were significantly expanded by integrating cellular resources and abiotic constraints in recent years. However, if unreasonable modeling methods were adopted due to a lack of consideration of biological knowledge, the conflicts between stoichiometric and other constraints, such as thermodynamic feasibility and enzyme resource availability, would lead to distorted predictions. In this work, we investigated a prediction anomaly of EcoETM, a constraints-based metabolic network model, and introduced the idea of enzyme compartmentalization into the analysis process. Through rational combination of reactions, we avoid the false prediction of pathway feasibility caused by the unrealistic assumption of free intermediate metabolites. This allowed us to correct the pathway structures of l-serine and l-tryptophan. A specific analysis explains the application method of the EcoETM-like model and demonstrates its potential and value in correcting the prediction results in pathway structure by resolving the conflict between different constraints and incorporating the evolved roles of enzymes as reaction compartments. Notably, this work also reveals the trade-off between product yield and thermodynamic feasibility. Our work is of great value for the structural improvement of constraints-based models.

Inhibitors of riboflavin biosynthetic pathway enzymes as potential antibacterial drugs

Multiple drug resistance is the main obstacle in the treatment of bacterial diseases. Resistance against antibiotics demands the exploration of new antimicrobial drug targets. A variety of in silico and genetic approaches show that the enzymes of the riboflavin biosynthetic pathway are crucial for the survival of bacteria. This pathway is absent in humans thus enzymes of the riboflavin biosynthetic pathway are emerging drug targets for resistant pathogenic bacterial strains. Exploring the structural details, their mechanism of action, intermediate elucidation, and interaction analysis would help in designing suitable inhibitors of these enzymes. The riboflavin biosynthetic pathway consists of seven distinct enzymes, namely, 3,4-dihydroxy-2-butanone 4-phosphate synthase, GTP cyclohydrolase II, pyrimidine deaminase/reductase, phosphatase, lumazine synthase, and riboflavin synthase. The present review summarizes the research work that has been carried out on these enzymes in terms of their structures, active site architectures, and molecular mechanism of catalysis. This review also walks through small molecule inhibitors that have been developed against several of these enzymes.

An interbacterial DNA deaminase toxin directly mutagenizes surviving target populations

When bacterial cells come in contact, antagonism mediated by the delivery of toxins frequently ensues. The potential for such encounters to have long-term beneficial consequences in recipient cells has not been investigated. Here, we examined the effects of intoxication by DddA, a cytosine deaminase delivered via the type VI secretion system (T6SS) of Burkholderia cenocepacia. Despite its killing potential, we observed that several bacterial species resist DddA and instead accumulate mutations. These mutations can lead to the acquisition of antibiotic resistance, indicating that even in the absence of killing, interbacterial antagonism can have profound consequences on target populations. Investigation of additional toxins from the deaminase superfamily revealed that mutagenic activity is a common feature of these proteins, including a representative we show targets single-stranded DNA and displays a markedly divergent structure. Our findings suggest that a surprising consequence of antagonistic interactions between bacteria could be the promotion of adaptation via the action of directly mutagenic toxins.

Maize Dek33 encodes a pyrimidine reductase in riboflavin biosynthesis that is essential for oil-body formation and ABA biosynthesis during seed development

The maize (Zea mays) defective kernel 33 (dek33) mutant produces defective and occasionally viviparous kernel phenotypes. In this study, we cloned Dek33 by positional cloning and found that it encodes a pyrimidine reductase in riboflavin biosynthesis. In dek33, a single-base mutation (G to A) in the C-terminal COG3236 domain caused a premature stop codon (TGA), producing a weak mutant allele with only a truncated form of the DEK33 protein that occurred at much lower levels that the completed WT form, and with a reduced riboflavin content. The dek33 mutation significantly affected oil-body formation and suppressed endoreduplication. It also disrupted ABA biosynthesis, resulting in lower ABA content that might be responsible for the viviparous embryo. In addition, our results indicated that the COG3236 domain is important for the protein stability of DEK33. Yeast two-hybrid experiments identified several proteins that interacted with DEK33, including RGLG2 and SnRK1, suggesting possible post-translational regulation of DEK33 stability. The interaction between DEK33 and these proteins was further confirmed by luciferase complementation image assays. This study provides a weak mutant allele that can be utilized to explore cellular responses to impaired riboflavin biosynthesis during seed development. Our findings indicate that the COG3236 domain might be an essential regulatory structure for DEK33 stability in maize.

Identification and characterization of the missing phosphatase on the riboflavin biosynthesis pathway in Arabidopsis thaliana

Despite the importance of riboflavin as the direct precursor of the cofactors flavin adenine dinucleotide (FAD) and flavin mononucleotide (FMN), the physiologically relevant catalyst dephosphorylating the riboflavin biosynthesis pathway intermediate 5-amino-6-ribitylamino-2,4(1H,3H) pyrimidinedione 5'-phosphate (ARPP) has not been characterized from any organism. By using as the query sequence a previously identified plastidial FMN hydrolase AtcpFHy1 (At1g79790), belonging to the haloacid dehalogenase (HAD) superfamily, seven candidates for the missing ARPP phosphatase were found, cloned, recombinantly expressed, and purified. Activity screening showed that the enzymes encoded by AtcpFHy1, At4g11570, and At4g25840 catalyze dephosphorylation of ARPP. AtcpFHy1 was renamed AtcpFHy/PyrP1, At4g11570 and At4g25840 were named AtPyrP2 and AtGpp1/PyrP3, respectively. Subcellular localization in planta indicated that AtPyrP2 was localized in plastids and AtGpp1/PyrP3 in mitochondria. Biochemical characterization of AtcpFHy/PyrP1 and AtPyrP2 showed that they have similar K_m values for the substrate ARPP, with AtcpFHy/PyrP1 having higher catalytic efficiency. Screening of 21 phosphorylated substrates showed that AtPyrP2 is specific for ARPP. Molecular weights of AtcpFHy/PyrP1 and AtPyrP2 were estimated at 46 and 72 kDa, suggesting dimers. pH and temperature optima for AtcpFHy/PyrP1 and AtPyrP2 were ~7.0-8.5 and 40-50°C. T-DNA knockout of AtcpFHy/PyrP1 did not affect the flavin profile of the transgenic plants, whereas silencing of AtPyrP2 decreased accumulation of riboflavin, FMN, and FAD. Our results strongly support AtPyrP2 as the missing phosphatase on the riboflavin biosynthesis pathway in Arabidopsis thaliana. The identification of this enzyme closes a long-standing gap in understanding of the riboflavin biosynthesis in plants.

‘Nothing of chemistry disappears in biology’: the Top 30 damage-prone endogenous metabolites

Many common metabolites are intrinsically unstable and reactive, and hence prone to chemical (i.e. non-enzymatic) damage in vivo. Although this fact is widely recognized, the purely chemical side-reactions of metabolic intermediates can be surprisingly hard to track down in the literature and are often treated in an unprioritized case-by-case way. Moreover, spontaneous chemical side-reactions tend to be overshadowed today by side-reactions mediated by promiscuous (‘sloppy’) enzymes even though chemical damage to metabolites may be even more prevalent than damage from enzyme sloppiness, has similar outcomes, and is held in check by similar biochemical repair or pre-emption mechanisms. To address these limitations and imbalances, here we draw together and systematically integrate information from the (bio)chemical literature, from cheminformatics, and from genome-scale metabolic models to objectively define a ‘Top 30’ list of damage-prone metabolites. A foundational part of this process was to derive general reaction rules for the damage chemistries involved. The criteria for a ‘Top 30’ metabolite included predicted chemical reactivity, essentiality, and occurrence in diverse organisms. We also explain how the damage chemistry reaction rules (‘operators’) are implemented in the Chemical-Damage-MINE (CD-MINE) database (minedatabase.mcs.anl.gov/#/top30) to provide a predictive tool for many additional potential metabolite damage products. Lastly, we illustrate how defining a ‘Top 30’ list can drive genomics-enabled discovery of the enzymes of previously unrecognized damage-control systems, and how applying chemical damage reaction rules can help identify previously unknown peaks in metabolomics profiles.

Metabolite Damage and Metabolite Damage Control in Plants

It is increasingly clear that (a) many metabolites undergo spontaneous or enzyme-catalyzed side reactions in vivo, (b) the damaged metabolites formed by these reactions can be harmful, and (c) organisms have biochemical systems that limit the buildup of damaged metabolites. These damage-control systems either return a damaged molecule to its pristine state (metabolite repair) or convert harmful molecules to harmless ones (damage preemption). Because all organisms share a core set of metabolites that suffer the same chemical and enzymatic damage reactions, certain damage-control systems are widely conserved across the kingdoms of life. Relatively few damage reactions and damage-control systems are well known. Uncovering new damage reactions and identifying the corresponding damaged metabolites, damage-control genes, and enzymes demands a coordinated mix of chemistry, metabolomics, cheminformatics, biochemistry, and comparative genomics. This review illustrates the above points using examples from plants, which are at least as prone to metabolite damage as other organisms.

Structural Insights into Mycobacterium tuberculosis Rv2671 Protein as a Dihydrofolate Reductase Functional Analogue Contributing to para-Aminosalicylic Acid Resistance

Mycobacterium tuberculosis (Mtb) Rv2671 is annotated as a 5-amino-6-ribitylamino-2,4(1H,3H)-pyrimidinedione 5'-phosphate (AROPP) reductase (RibD) in the riboflavin biosynthetic pathway. Recently, a strain of Mtb with a mutation in the 5' untranslated region of Rv2671, which resulted in its overexpression, was found to be resistant to dihydrofolate reductase (DHFR) inhibitors including the anti-Mtb drug para-aminosalicylic acid (PAS). In this study, a biochemical analysis of Rv2671 showed that it was able to catalyze the reduction of dihydrofolate (DHF) to tetrahydrofolate (THF), which explained why the overexpression of Rv2671 was sufficient to confer PAS resistance. We solved the structure of Rv2671 in complex with the NADP(+) and tetrahydrofolate (THF), which revealed the structural basis for the DHFR activity. The structures of Rv2671 complexed with two DHFR inhibitors, trimethoprim and trimetrexate, provided additional details of the substrate binding pocket and elucidated the differences between their inhibitory activities. Finally, Rv2671 was unable to catalyze the reduction of AROPP, which indicated that Rv2671 and its closely related orthologues are not involved in riboflavin biosynthesis.

A directed-overflow and damage-control N-glycosidase in riboflavin biosynthesis

Plants and bacteria synthesize the essential human micronutrient riboflavin (vitamin B2) via the same multi-step pathway. The early intermediates of this pathway are notoriously reactive and may be overproduced in vivo because riboflavin biosynthesis enzymes lack feedback controls. In the present paper, we demonstrate disposal of riboflavin intermediates by COG3236 (DUF1768), a protein of previously unknown function that is fused to two different riboflavin pathway enzymes in plants and bacteria (RIBR and RibA respectively). We present cheminformatic, biochemical, genetic and genomic evidence to show that: (i) plant and bacterial COG3236 proteins cleave the N-glycosidic bond of the first two intermediates of riboflavin biosynthesis, yielding relatively innocuous products; (ii) certain COG3236 proteins are in a multi-enzyme riboflavin biosynthesis complex that gives them privileged access to riboflavin intermediates; and (iii) COG3236 action in Arabidopsis thaliana and Escherichia coli helps maintain flavin levels. COG3236 proteins thus illustrate two emerging principles in chemical biology: directed overflow metabolism, in which excess flux is diverted out of a pathway, and the pre-emption of damage from reactive metabolites.

Covalent docking predicts substrates for haloalkanoate dehalogenase superfamily phosphatases

Enzyme function prediction remains an important open problem. Though structure-based modeling, such as metabolite docking, can identify substrates of some enzymes, it is ill-suited to reactions that progress through a covalent intermediate. Here we investigated the ability of covalent docking to identify substrates that pass through such a covalent intermediate, focusing particularly on the haloalkanoate dehalogenase superfamily. In retrospective assessments, covalent docking recapitulated substrate binding modes of known cocrystal structures and identified experimental substrates from a set of putative phosphorylated metabolites. In comparison, noncovalent docking of high-energy intermediates yielded nonproductive poses. In prospective predictions against seven enzymes, a substrate was identified for five. For one of those cases, a covalent docking prediction, confirmed by empirical screening, and combined with genomic context analysis, suggested the identity of the enzyme that catalyzes the orphan phosphatase reaction in the riboflavin biosynthetic pathway of Bacteroides.

Kinetic modeling of riboflavin biosynthesis in Bacillus subtilis under production conditions

To study the network dynamics of the riboflavin biosynthesis pathway and to identify potential bottlenecks in the system, an ordinary differential equation-based model was constructed using available literature data for production strains. The results confirmed that the RibA protein is rate limiting in the pathway. Under the conditions investigated, we determined a potential limiting order of the remaining enzymes under increased RibA concentration (>0.102 mM) and therefore higher riboflavin production (>0.045 mmol g(CDW)(-1) h(-1) and 0.0035 mM s(-1), respectively). The reductase activity of RibG and lumazine synthase (RibH) might be the next most limiting steps. The computational minimization of the enzyme concentrations of the pathway suggested the need for a greater RibH concentration (0.251 mM) compared with the other enzymes (RibG: 0.188 mM, RibB: 0.023 mM).

Structure of diaminohydroxyphosphoribosylaminopyrimidine deaminase/5-amino-6-(5-phosphoribosylamino)uracil reductase from Acinetobacter baumannii

The bifunctional diaminohydroxyphosphoribosylaminopyrimidine deaminase/5-amino-6-(5-phosphoribosylamino)uracil reductase (RibD) represents a potential antibacterial drug target. The structure of recombinant Acinetobacter baumannii RibD is reported in orthorhombic and tetragonal crystal forms at 2.2 and 2.0 Å resolution, respectively. Comparisons with orthologous structures in the Protein Data Bank indicated close similarities. The tetragonal crystal form was obtained in the presence of guanosine monophosphate, which surprisingly was observed to occupy the adenine-binding site of the reductase domain.

Structural and functional insights into Saccharomyces cerevisiae riboflavin biosynthesis reductase RIB7

Saccharomyces cerevisiae RIB7 (ScRIB7) is a potent target for anti-fungal agents because of its involvement in the riboflavin biosynthesis pathway as a NADPH-dependent reductase. However, the catalytic mechanism of riboflavin biosynthesis reductase (RBSRs) is controversial, and enzyme structure information is still lacking in eukaryotes. Here we report the crystal structure of Saccharomyces cerevisiae RIB7 at 2.10 Å resolution and its complex with NADPH at 2.35 Å resolution. ScRIB7 exists as a stable homodimer, and each subunit consists of nine central β-sheets flanked by five helices, resembling the structure of RIB7 homologues. A conserved G(76)-X-G(78)-Xn-G(181)-G(182) motif is present at the NADPH pyrophosphate group binding site. Activity assays confirmed the necessity of Thr79, Asp83, Glu180 and Gly182 for the activity of ScRIB7. Substrate preference of ScRIB7 was altered by mutating one residue (Thr35) to a Lysine, implying that ScRIB7 Thr35 and its corresponding residue, a lysine in bacteria, are important in substrate-specific recognition.

Identification and characterization of the missing pyrimidine reductase in the plant riboflavin biosynthesis pathway

Riboflavin (vitamin B₂) is the precursor of the flavin coenzymes flavin mononucleotide and flavin adenine dinucleotide. In Escherichia coli and other bacteria, sequential deamination and reduction steps in riboflavin biosynthesis are catalyzed by RibD, a bifunctional protein with distinct pyrimidine deaminase and reductase domains. Plants have two diverged RibD homologs, PyrD and PyrR; PyrR proteins have an extra carboxyl-terminal domain (COG3236) of unknown function. Arabidopsis (Arabidopsis thaliana) PyrD (encoded by At4g20960) is known to be a monofunctional pyrimidine deaminase, but no pyrimidine reductase has been identified. Bioinformatic analyses indicated that plant PyrR proteins have a catalytically competent reductase domain but lack essential zinc-binding residues in the deaminase domain, and that the Arabidopsis PyrR gene (At3g47390) is coexpressed with riboflavin synthesis genes. These observations imply that PyrR is a pyrimidine reductase without deaminase activity. Consistent with this inference, Arabidopsis or maize (Zea mays) PyrR (At3g47390 or GRMZM2G090068) restored riboflavin prototrophy to an E. coli ribD deletant strain when coexpressed with the corresponding PyrD protein (At4g20960 or GRMZM2G320099) but not when expressed alone; the COG3236 domain was unnecessary for complementing activity. Furthermore, recombinant maize PyrR mediated NAD(P)H-dependent pyrimidine reduction in vitro. Import assays with pea (Pisum sativum) chloroplasts showed that PyrR and PyrD are taken up and proteolytically processed. Ablation of the maize PyrR gene caused early seed lethality. These data argue that PyrR is the missing plant pyrimidine reductase, that it is plastid localized, and that it is essential. The role of the COG3236 domain remains mysterious; no evidence was obtained for the possibility that it catalyzes the dephosphorylation that follows pyrimidine reduction.

Genetic control of biosynthesis and transport of riboflavin and flavin nucleotides and construction of robust biotechnological producers

Riboflavin [7,8-dimethyl-10-(1'-d-ribityl)isoalloxazine, vitamin B₂] is an obligatory component of human and animal diets, as it serves as the precursor of flavin coenzymes, flavin mononucleotide, and flavin adenine dinucleotide, which are involved in oxidative metabolism and other processes. Commercially produced riboflavin is used in agriculture, medicine, and the food industry. Riboflavin synthesis starts from GTP and ribulose-5-phosphate and proceeds through pyrimidine and pteridine intermediates. Flavin nucleotides are synthesized in two consecutive reactions from riboflavin. Some microorganisms and all animal cells are capable of riboflavin uptake, whereas many microorganisms have distinct systems for riboflavin excretion to the medium. Regulation of riboflavin synthesis in bacteria occurs by repression at the transcriptional level by flavin mononucleotide, which binds to nascent noncoding mRNA and blocks further transcription (named the riboswitch). In flavinogenic molds, riboflavin overproduction starts at the stationary phase and is accompanied by derepression of enzymes involved in riboflavin synthesis, sporulation, and mycelial lysis. In flavinogenic yeasts, transcriptional repression of riboflavin synthesis is exerted by iron ions and not by flavins. The putative transcription factor encoded by SEF1 is somehow involved in this regulation. Most commercial riboflavin is currently produced or was produced earlier by microbial synthesis using special selected strains of Bacillus subtilis, Ashbya gossypii, and Candida famata. Whereas earlier RF overproducers were isolated by classical selection, current producers of riboflavin and flavin nucleotides have been developed using modern approaches of metabolic engineering that involve overexpression of structural and regulatory genes of the RF biosynthetic pathway as well as genes involved in the overproduction of the purine precursor of riboflavin, GTP.