Molecular identification of pseudouridine-metabolizing enzymes

RudS: bacterial desulfidase responsible for tRNA 4-thiouridine de-modification

In this study, we present an extensive analysis of a widespread group of bacterial tRNA de-modifying enzymes, dubbed RudS, which consist of a TudS desulfidase fused to a Domain of Unknown Function 1722 (DUF1722). RudS enzymes exhibit specific de-modification activity towards the 4-thiouridine modification (s4U) in tRNA molecules, as indicated by our experimental findings. The heterologous overexpression of RudS genes in Escherichia coli significantly reduces the tRNA 4-thiouridine content and diminishes UVA-induced growth delay, indicating the enzyme's role in regulating photosensitive tRNA s4U modification. Through a combination of protein modeling, docking studies, and molecular dynamics simulations, we have identified amino acid residues involved in catalysis and tRNA binding. Experimental validation through targeted mutagenesis confirms the TudS domain as the catalytic core of RudS, with the DUF1722 domain facilitating tRNA binding in the anticodon region. Our results suggest that RudS tRNA modification eraser proteins may play a role in regulating tRNA during prokaryotic stress responses.

Elucidation of the Gemcitabine Transporters of Escherichia coli K-12 and Gamma-Proteobacteria Linked to Gemcitabine-Related Chemoresistance

Gemcitabine (2',2'-difluoro-2'-deoxycytidine), a widely used anticancer drug, is considered a gold standard in treating aggressive pancreatic cancers. Gamma-proteobacteria that colonize the pancreatic tumors contribute to chemoresistance against gemcitabine by metabolizing the drug to a less active and deaminated form. The gemcitabine transporters of these bacteria are unknown to date. Furthermore, there is no complete knowledge of the gemcitabine transporters in Escherichia coli or any other related proteobacteria. In this study, we investigate the complement of gemcitabine transporters in E. coli K-12 and two common chemoresistance-related bacteria (Klebsiella pneumoniae and Citrobacter freundii). We found that E. coli K-12 has two high-affinity gemcitabine transporters with distinct specificity properties, namely, NupC and NupG, whereas the gemcitabine transporters of C. freundii and K. pneumoniae include the NupC and NupG orthologs, functionally indistinguishable from their counterparts, and, in K. pneumoniae, one additional NupC variant, designated KpNupC2. All these bacterial transporters have a higher affinity for gemcitabine than their human counterparts. The highest affinity (KM 2.5-3.0 μΜ) is exhibited by NupGs of the bacteria-specific nucleoside-H+ symporter (NHS) family followed by NupCs (KM 10-13 μΜ) of the concentrative nucleoside transporter (CNT) family, 15-100 times higher than the affinities reported for the human gemcitabine transporter hENT1/SLC29A1, which is primarily associated with gemcitabine uptake in the pancreatic adenocarcinoma cells. Our results offer a basis for further insight into the role of specific bacteria in drug availability within tumors and for understanding the structure-function differences of bacterial and human drug transporters.

Efficient Synthetic Access to Stable Isotope Labelled Pseudouridine Phosphoramidites for RNA NMR Spectroscopy

Here we report the efficient synthetic access to 13C/15N-labelled pseudouridine phosphoramidites, which were incorporated into a binary H/ACA box guide RNA/product complex comprising 77 nucleotides (nts) in total and into a 75 nt E. coli tRNAGly. The stable isotope (SI) labelled pseudouridines were produced via a highly efficient chemo-enzymatic synthesis. 13C/15N labelled uracils were produced via chemical synthesis and enzymatically converted to pseudouridine 5'-monophosphate (ΨMP) by using YeiN, a Ψ-5'-monophosphate C-glycosidase. Removal of the 5'-phosphate group yielded the desired pseudouridine nucleoside (Ψ), which was transformed into a phosphoramidite building suitable for RNA solid phase synthesis. A Ψ -building block carrying both a 13C and a 15N label was incorporated into a product RNA and the complex formation with a 63 nt H/ACA box RNA could be observed via NMR. Furthermore, the SI labelled pseudouridine building block was used to determine imino proton bulk water exchange rates of a 75 nt E. coli tRNAGly CCmnm5U, identifying the TΨC-loop 5-methyluridine as a modifier of the exchange rates. The efficient synthetic access to SI-labelled Ψ building blocks will allow the solution and solid-state NMR spectroscopic studies of Ψ containing RNAs and will facilitate the mass spectrometric analysis of Ψ-modified nucleic acids.

Impact of suppression of the SOS response on protein expression in clinical isolates of Escherichia coli under antimicrobial pressure of ciprofloxacin

Introduction/objective

Suppression of the SOS response in combination with drugs damaging DNA has been proposed as a potential target to tackle antimicrobial resistance. The SOS response is the pathway used to repair bacterial DNA damage induced by antimicrobials such as quinolones. The extent of lexA-regulated protein expression and other associated systems under pressure of agents that damage bacterial DNA in clinical isolates remains unclear. The aim of this study was to assess the impact of this strategy consisting on suppression of the SOS response in combination with quinolones on the proteome profile of Escherichia coli clinical strains.

Materials and methods

Five clinical isolates of E. coli carrying different chromosomally- and/or plasmid-mediated quinolone resistance mechanisms with different phenotypes were selected, with E. coli ATCC 25922 as control strain. In addition, from each clinical isolate and control, a second strain was created, in which the SOS response was suppressed by deletion of the recA gene. Bacterial inocula from all 12 strains were then exposed to 1xMIC ciprofloxacin treatment (relative to the wild-type phenotype for each isogenic pair) for 1 h. Cell pellets were collected, and proteins were digested into peptides using trypsin. Protein identification and label-free quantification were done by liquid chromatography-mass spectrometry (LC-MS) in order to identify proteins that were differentially expressed upon deletion of recA in each strain. Data analysis and statistical analysis were performed using the MaxQuant and Perseus software.

Results

The proteins with the lowest expression levels were: RecA (as control), AphA, CysP, DinG, DinI, GarL, PriS, PsuG, PsuK, RpsQ, UgpB and YebG; those with the highest expression levels were: Hpf, IbpB, TufB and RpmH. Most of these expression alterations were strain-dependent and involved DNA repair processes and nucleotide, protein and carbohydrate metabolism, and transport. In isolates with suppressed SOS response, the number of underexpressed proteins was higher than overexpressed proteins.

Conclusion

High genomic and proteomic variability was observed among clinical isolates and was not associated with a specific resistant phenotype. This study provides an interesting approach to identify new potential targets to combat antimicrobial resistance.

Biosynthesis and function of 7-deazaguanine derivatives in bacteria and phages

SUMMARYDeazaguanine modifications play multifaceted roles in the molecular biology of DNA and tRNA, shaping diverse yet essential biological processes, including the nuanced fine-tuning of translation efficiency and the intricate modulation of codon-anticodon interactions. Beyond their roles in translation, deazaguanine modifications contribute to cellular stress resistance, self-nonself discrimination mechanisms, and host evasion defenses, directly modulating the adaptability of living organisms. Deazaguanine moieties extend beyond nucleic acid modifications, manifesting in the structural diversity of biologically active natural products. Their roles in fundamental cellular processes and their presence in biologically active natural products underscore their versatility and pivotal contributions to the intricate web of molecular interactions within living organisms. Here, we discuss the current understanding of the biosynthesis and multifaceted functions of deazaguanines, shedding light on their diverse and dynamic roles in the molecular landscape of life.

"Boomerang" Strategy in Carbohydrate Chemistry: Diastereoselective Synthesis of C-Glycosylated Benzothiazoles from ortho-Isocyanophenyl Thioglycosides

This paper reveals a novel "boomerang" strategy in the expedient and diastereoselective synthesis of C-nucleoside analogues. Bench-stable ortho-isocyanophenyl thioglycosides can be converted to glycosyl radicals through rapid and efficient C-S bond homolysis when they are irradiated by visible light. The glycosyl radicals are subsequently trapped by the corresponding leaving group or other radical acceptors to provide diverse C-nucleoside analogues under mild conditions.

Structure and function of the pseudouridine 5'-monophosphate glycosylase PUMY from Arabidopsis thaliana

Pseudouridine is a noncanonical C-nucleoside containing a C-C glycosidic linkage between uracil and ribose. In the two-step degradation of pseudouridine, pseudouridine 5'-monophosphate glycosylase (PUMY) is responsible for the second step and catalyses the cleavage of the C-C glycosidic bond in pseudouridine 5'-monophosphate (ΨMP) into uridine and ribose 5'-phosphate, which are recycled via other metabolic pathways. Structural features of Escherichia coli PUMY have been reported, but the details of the substrate specificity of ΨMP were unknown. Here, we present three crystal structures of Arabidopsis thaliana PUMY in different ligation states and a kinetic analysis of ΨMP degradation. The results indicate that Thr149 and Asn308, which are conserved in the PUMY family, are structural determinants for recognizing the nucleobase of ΨMP. The distinct binding modes of ΨMP and ribose 5'-phosphate also suggest that the nucleobase, rather than the phosphate group, of ΨMP dictates the substrate-binding mode. An open-to-close transition of the active site is essential for catalysis, which is mediated by two α-helices, α11 and α12, near the active site. Mutational analysis validates the proposed roles of the active site residues in catalysis. Our structural and functional analyses provide further insight into the enzymatic features of PUMY towards ΨMP.

Engineering a synthetic energy-efficient formaldehyde assimilation cycle in Escherichia coli

One-carbon (C1) substrates, such as methanol or formate, are attractive feedstocks for circular bioeconomy. These substrates are typically converted into formaldehyde, serving as the entry point into metabolism. Here, we design an erythrulose monophosphate (EuMP) cycle for formaldehyde assimilation, leveraging a promiscuous dihydroxyacetone phosphate dependent aldolase as key enzyme. In silico modeling reveals that the cycle is highly energy-efficient, holding the potential for high bioproduct yields. Dissecting the EuMP into four modules, we use a stepwise strategy to demonstrate in vivo feasibility of the modules in E. coli sensor strains with sarcosine as formaldehyde source. From adaptive laboratory evolution for module integration, we identify key mutations enabling the accommodation of the EuMP reactions with endogenous metabolism. Overall, our study demonstrates the proof-of-concept for a highly efficient, new-to-nature formaldehyde assimilation pathway, opening a way for the development of a methylotrophic platform for a C1-fueled bioeconomy in the future.

RNA pseudouridine modification in plants

Pseudouridine is one of the well-known chemical modifications in various RNA species. Current advances to detect pseudouridine show that the pseudouridine landscape is dynamic and affects multiple cellular processes. Although our understanding of this post-transcriptional modification mainly depends on yeast and human models, the recent findings provide strong evidence for the critical role of pseudouridine in plants. Here, we review the current knowledge of pseudouridine in plant RNAs, including its synthesis, degradation, regulatory mechanisms, and functions. Moreover, we propose future areas of research on pseudouridine modification in plants.

Metabolic Engineering of Escherichia coli for Efficient Production of Pseudouridine

Pseudouridine-incorporated mRNA vaccines can enhance protein expression and reduce immunogenicity, leading to a high demand for pseudouridine to be used in mRNA drug production. To achieve the low-cost production of pseudouridine, Escherichia coli was systematically modified to utilize inexpensive raw materials to efficiently produce pseudouridine. First, in the pyrimidine biosynthesis pathway, genes related to the precursor competing pathway and the negative regulator were deleted, which increased pseudouridine production. Second, two critical genes, pseudouridine-5'-phosphate glycosidase (psuG) and phosphatase genes from different bacteria, were screened and employed in various genetic constructs, and the pseudouridine yield of the optical strain increased to 599 mg/L. The accumulation of pseudouridine was further increased by the deletion of pseudouridine catabolism-related genes. Ultimately, the pseudouridine titer in a 5 L bioreactor reached 7.9 g/L, and the yield of pseudouridine on glucose was 0.15 g/g. Overall, a cell factory producing pseudouridine was successfully constructed and showed potential for industrial production.

RNA modification: mechanisms and therapeutic targets

RNA modifications are dynamic and reversible chemical modifications on substrate RNA that are regulated by specific modifying enzymes. They play important roles in the regulation of many biological processes in various diseases, such as the development of cancer and other diseases. With the help of advanced sequencing technologies, the role of RNA modifications has caught increasing attention in human diseases in scientific research. In this review, we briefly summarized the basic mechanisms of several common RNA modifications, including m6A, m5C, m1A, m7G, Ψ, A-to-I editing and ac4C. Importantly, we discussed their potential functions in human diseases, including cancer, neurological disorders, cardiovascular diseases, metabolic diseases, genetic and developmental diseases, as well as immune disorders. Through the "writing-erasing-reading" mechanisms, RNA modifications regulate the stability, translation, and localization of pivotal disease-related mRNAs to manipulate disease development. Moreover, we also highlighted in this review all currently available RNA-modifier-targeting small molecular inhibitors or activators, most of which are designed against m6A-related enzymes, such as METTL3, FTO and ALKBH5. This review provides clues for potential clinical therapy as well as future study directions in the RNA modification field. More in-depth studies on RNA modifications, their roles in human diseases and further development of their inhibitors or activators are needed for a thorough understanding of epitranscriptomics as well as diagnosis, treatment, and prognosis of human diseases.

Multiplatform Metabolomics Characterization Reveals Novel Metabolites and Phospholipid Compositional Rules of Haemophilus influenzae Rd KW20

Haemophilus influenzae is a gram-negative bacterium of relevant clinical interest. H. influenzae Rd KW20 was the first organism to be sequenced and for which a genome-scale metabolic model (GEM) was developed. However, current H. influenzae GEMs are unable to capture several aspects of metabolome nature related to metabolite pools. To directly and comprehensively characterize the endometabolome of H. influenzae Rd KW20, we performed a multiplatform MS-based metabolomics approach combining LC-MS, GC-MS and CE-MS. We obtained direct evidence of 15-20% of the endometabolome present in current H. influenzae GEMs and showed that polar metabolite pools are interconnected through correlating metabolite islands. Notably, we obtained high-quality evidence of 18 metabolites not previously included in H. influenzae GEMs, including the antimicrobial metabolite cyclo(Leu-Pro). Additionally, we comprehensively characterized and evaluated the quantitative composition of the phospholipidome of H. influenzae, revealing that the fatty acyl chain composition is largely independent of the lipid class, as well as that the probability distribution of phospholipids is mostly related to the conditional probability distribution of individual acyl chains. This finding enabled us to provide a rationale for the observed phospholipid profiles and estimate the abundance of low-level species, permitting the expansion of the phospholipidome characterization through predictive probabilistic modelling.

A selective and atom-economic rearrangement of uridine by cascade biocatalysis for production of pseudouridine

As a crucial factor of their therapeutic efficacy, the currently marketed mRNA vaccines feature uniform substitution of uridine (U) by the corresponding C-nucleoside, pseudouridine (Ψ), in 1-N-methylated form. Synthetic supply of the mRNA building block (1-N-Me-Ψ-5'-triphosphate) involves expedient access to Ψ as the principal challenge. Here, we show selective and atom-economic 1N-5C rearrangement of β-D-ribosyl on uracil to obtain Ψ from unprotected U in quantitative yield. One-pot cascade transformation of U in four enzyme-catalyzed steps, via D-ribose (Rib)-1-phosphate, Rib-5-phosphate (Rib5P) and Ψ-5'-phosphate (ΨMP), gives Ψ. Coordinated function of the coupled enzymes in the overall rearrangement necessitates specific release of phosphate from the ΨMP, but not from the intermediary ribose phosphates. Discovery of Yjjg as ΨMP-specific phosphatase enables internally controlled regeneration of phosphate as catalytic reagent. With driving force provided from the net N-C rearrangement, the optimized U reaction yields a supersaturated product solution (∼250 g/L) from which the pure Ψ crystallizes (90% recovery). Scale up to 25 g isolated product at enzyme turnovers of ∼10⁵mol/mol demonstrates a robust process technology, promising for Ψ production. Our study identifies a multistep rearrangement reaction, realized by cascade biocatalysis, for C-nucleoside synthesis in high efficiency.

Biocatalytic cascade transformations for the synthesis of C-nucleosides and N-nucleoside analogs

Nucleosides and their analogs, including those that feature substitution of the canonical N-glycosidic by a C-glycosidic linkage, provide access to potent antiviral, antibacterial, and antitumor drugs. Furthermore, they are key building blocks of m-RNA vaccines and play a crucial role for vaccine therapeutic effectiveness. As the medicinal applications of nucleosides increase in number and importance, there is a growing need for efficiency-enhanced routes of nucleoside synthesis. Cascade biocatalysis, that is, the application of natural or evolved enzymes promoting complex transformations in multiple steps in one pot and without the need of intermediate purification, emerges as a powerful tool to obtain nucleosides from readily available starting materials. Recent efforts in enzyme discovery and protein engineering expand the toolbox of catalysts active toward nucleosides or nucleotides. In this review, we highlight recent applications, and discuss challenges, of cascade biocatalysis for nucleoside synthesis. We focus on C-nucleosides and important analogs of the canonical N-nucleosides.

Translating the Arabidopsis thaliana Peroxisome Proteome Insights to Solanum lycopersicum: Consensus Versus Diversity

Peroxisomes are small, single-membrane specialized organelles present in all eukaryotic organisms. The peroxisome is one of the nodal centers of reactive oxygen species homeostasis in plants, which are generated in a high amount due to various stress conditions. Over the past decade, there has been extensive study on peroxisomal proteins and their signaling pathways in the model plant Arabidopsis thaliana, and a lot has been deciphered. However, not much impetus has been given to studying the peroxisome proteome of economically important crops. Owing to the significance of peroxisomes in the physiology of plants during normal and stress conditions, understating its proteome is of much importance. Hence, in this paper, we have made a snapshot of putative peroxisomal matrix proteins in the economically important vegetable crop tomato (Solanum lycopersicum, (L.) family Solanaceae). First, a reference peroxisomal matrix proteome map was generated for Arabidopsis thaliana using the available proteomic and localization studies, and proteins were categorized into various groups as per their annotations. This was used to create the putative peroxisomal matrix proteome map for S. lycopersicum. The putative peroxisome proteome in S. lycopersicum retains the basic framework: the bulk of proteins had peroxisomal targeting signal (PTS) type 1, a minor group had PTS2, and the catalase family retained its characteristic internal PTS. Apart from these, a considerable number of S. lycopersicum orthologs did not contain any "obvious" PTS. The number of PTS2 isoforms was found to be reduced in S. lycopersicum. We further investigated the PTS1s in the case of both the plant species and generated a pattern for canonical and non-canonical PTS1s. The number of canonical PTS1 proteins was comparatively lesser in S. lycopersicum. The non-canonical PTS1s were found to be comparable in both the plant species; however, S. lycopersicum showed greater diversity in the composition of the signal tripeptide. Finally, we have tried to address the lacunas and probable strategies to fill those gaps.

De novo assembly provides new insights into the evolution of Elaeagnus angustifolia L

BACKGROUND

Elaeagnus angustifolia L. is a deciduous tree in the family Elaeagnaceae. It is widely used to study abiotic stress tolerance in plants and to improve desertification-affected land because of its ability to withstand diverse types of environmental stress, such as drought, salt, cold, and wind. However, no studies have examined the mechanisms underlying the resistance of E. angustifolia to environmental stress and its adaptive evolution.

METHODS

Here, we used PacBio, Hi-C, resequencing, and RNA-seq to construct the genome and transcriptome of E. angustifolia and explore its adaptive evolution.

RESULTS

The reconstructed genome of E. angustifolia was 526.80 Mb, with a contig N50 of 12.60 Mb and estimated divergence time of 84.24 Mya. Gene family expansion and resequencing analyses showed that the evolution of E. angustifolia was closely related to environmental conditions. After exposure to salt stress, GO pathway analysis showed that new genes identified from the transcriptome were related to ATP-binding, metal ion binding, and nucleic acid binding.

CONCLUSION

The genome sequence of E. angustifolia could be used for comparative genomic analyses of Elaeagnaceae family members and could help elucidate the mechanisms underlying the response of E. angustifolia to drought, salt, cold, and wind stress. Generally, these results provide new insights that could be used to improve desertification-affected land.

Structure Characterization of Escherichia coli Pseudouridine Kinase PsuK

Pseudouridine (Ψ) is one of the most abundant RNA modifications in cellular RNAs that post-transcriptionally impact many aspects of RNA. However, the metabolic fate of modified RNA nucleotides has long been a question. A pseudouridine kinase (PsuK) and a pseudouridine monophosphate glycosylase (PsuG) in Escherichia coli were first characterized as involved in pseudouridine degradation by catalyzing the phosphorylation of pseudouridine to pseudouridine 5′-phosphate (ΨMP) and further hydrolyzing 5′-ΨMP to produce uracil and ribose 5′-phosphate. Recently, their homolog proteins in eukaryotes were also identified, which were named PUKI and PUMY in Arabidopsis. Here, we solved the crystal structures of apo-EcPsuK and its binary complex with Ψ or N¹-methyl-pseudouridine (m1Ψ). The structure of EcPsuK showed a homodimer conformation assembled by its β-thumb region. EcPsuK has an appropriate binding site with a series of hydrophilic and hydrophobic interactions for Ψ. Moreover, our complex structure of EcPsuK-m1Ψ suggested the binding pocket has an appropriate capacity for m1Ψ. We also identified the monovalent ion-binding site and potential ATP-binding site. Our studies improved the understanding of the mechanism of Ψ turnover.

Proteomes of Uropathogenic Escherichia coli Growing in Human Urine and in J82 Urinary Bladder Cells

Uropathogenic Escherichia coli (UPEC) are the most common cause of urinary tract infection (UTI). UPEC normally reside in the intestine, and during establishment of UTI, they undergo metabolic adaptations, first to urine and then upon tissue invasion to the bladder cell interior. To understand these adaptations, we used quantitative proteomic profiling to characterize protein expression of the UPEC strain UTI89 growing in human urine and when inside J82 bladder cells. In order to facilitate detection of UPEC proteins over the excess amount of eukaryotic proteins in bladder cells, we developed a method where proteins from UTI89 grown in MOPS and urine was spiked-in to enhance detection of bacterial proteins. More than 2000 E. coli proteins were detected. During growth in urine, proteins associated with iron acquisition and several amino acid uptake and biosynthesis systems, most prominently arginine metabolism, were significantly upregulated. During growth in J82 cells, proteins related to iron uptake and arginine metabolisms were likewise upregulated together with proteins involved in sulfur compound turnover. Ribosomal proteins were downregulated relative to growth in MOPS in this environment. There was no direct correlation between upregulated proteins and proteins reported to be essential for infections, showing that upregulation during growth does not signify that the proteins are essential for growth under a condition.

Substrate-binding loop interactions with pseudouridine trigger conformational changes that promote catalytic efficiency of pseudouridine kinase PUKI

Pseudouridine, one major RNA modification, is catabolized into uracil and ribose-5′-phosphate by two sequential enzymatic reactions. In the first step, pseudouridine kinase (PUKI) phosphorylates pseudouridine to pseudouridine 5′-monophosphate. High-fidelity catalysis of pseudouridine by PUKI prevents possible disturbance of in vivo pyrimidine homeostasis. However, the molecular basis of how PUKI selectively phosphorylates pseudouridine over uridine with >100-fold greater efficiency despite minor differences in their K_m values has not been elucidated. To investigate this selectivity, in this study we determined the structures of PUKI from Escherichia coli strain B (EcPUKI) in various ligation states. The structure of EcPUKI was determined to be similar to PUKI from Arabidopsis thaliana, including an α/β core domain and β-stranded small domain, with dimerization occurring via the β-stranded small domain. In a binary complex, we show that Ser30 in the substrate-binding loop of the small domain mediates interactions with the hallmark N1 atom of pseudouridine nucleobase, causing conformational changes in its quaternary structure. Kinetic and fluorescence spectroscopic analyses also showed that the Ser30-mediated interaction is a prerequisite for conformational changes and subsequent catalysis by EcPUKI. Furthermore, S30A mutation or EcPUKI complexed with other nucleosides homologous to pseudouridine but lacking the pseudouridine-specific N1 atom did not induce such conformational changes, demonstrating the catalytic significance of the proposed Ser30-mediated interaction. These analyses provide structural and functional evidence for a pseudouridine-dependent conformational change of EcPUKI and its functional linkage to catalysis.

Semi-enzymatic synthesis of pseudouridine

Modifications of RNA molecules have a significant effect on their structure and function. One of the most common modifications is the isomerization from uridine to pseudouridine. Despite its prevalence in natural RNA sequences, organic synthesis of pseudouridine has been challenging because of the stereochemistry requirement and the sensitivity of reaction steps to moisture. Herein, a semi-enzymatic synthetic route is developed for the synthesis of pseudouridine using adenosine 5'-monophosphate and uracil as the starting materials and a reverse reaction catalyzed by the pseudouridine monophosphate glycosidase. This synthetic route has only three steps and the overall yield of β-pseudouridine production was 68.4%.

Structural basis for the substrate specificity and catalytic features of pseudouridine kinase from Arabidopsis thaliana

RNA modifications can regulate the stability of RNAs, mRNA-protein interactions, and translation efficiency. Pseudouridine is a prevalent RNA modification, and its metabolic fate after RNA turnover was recently characterized in eukaryotes, in the plant Arabidopsis thaliana. Here, we present structural and biochemical analyses of PSEUDOURIDINE KINASE from Arabidopsis (AtPUKI), the enzyme catalyzing the first step in pseudouridine degradation. AtPUKI, a member of the PfkB family of carbohydrate kinases, is a homodimeric α/β protein with a protruding small β-strand domain, which serves simultaneously as dimerization interface and dynamic substrate specificity determinant. AtPUKI has a unique nucleoside binding site specifying the binding of pseudourine, in particular at the nucleobase, by multiple hydrophilic interactions, of which one is mediated by a loop from the small β-strand domain of the adjacent monomer. Conformational transition of the dimerized small β-strand domains containing active site residues is required for substrate specificity. These dynamic features explain the higher catalytic efficiency for pseudouridine over uridine. Both substrates bind well (similar Km), but only pseudouridine is turned over efficiently. Our studies provide an example for structural and functional divergence in the PfkB family and highlight how AtPUKI avoids futile uridine phosphorylation which in vivo would disturb pyrimidine homeostasis.

Reverse C-glycosidase reaction provides C-nucleotide building blocks of xenobiotic nucleic acids

C-Analogues of the canonical N-nucleosides have considerable importance in medicinal chemistry and are promising building blocks of xenobiotic nucleic acids (XNA) in synthetic biology. Although well established for synthesis of N-nucleosides, biocatalytic methods are lacking in C-nucleoside synthetic chemistry. Here, we identify pseudouridine monophosphate C-glycosidase for selective 5-β-C-glycosylation of uracil and derivatives thereof from pentose 5-phosphate (D-ribose, 2-deoxy-D-ribose, D-arabinose, D-xylose) substrates. Substrate requirements of the enzymatic reaction are consistent with a Mannich-like addition between the pyrimidine nucleobase and the iminium intermediate of enzyme (Lys166) and open-chain pentose 5-phosphate. β-Elimination of the lysine and stereoselective ring closure give the product. We demonstrate phosphorylation-glycosylation cascade reactions for efficient, one-pot synthesis of C-nucleoside phosphates (yield: 33 - 94%) from unprotected sugar and nucleobase. We show incorporation of the enzymatically synthesized C-nucleotide triphosphates into nucleic acids by RNA polymerase. Collectively, these findings implement biocatalytic methodology for C-nucleotide synthesis which can facilitate XNA engineering for synthetic biology applications.

Mutation accumulation and horizontal gene transfer in Escherichia coli colonizing the gut of old mice

The ecology and environment of the microbes that inhabit the mammalian intestine undergoes several changes as the host ages. Here, we ask if the selection pressure experienced by a new strain colonizing the aging gut differs from that in the gut of young adults. Using experimental evolution in mice after a short antibiotic treatment, as a model for a common clinical situation, we show that a new colonizing E. coli strain rapidly adapts to the aging gut via both mutation accumulation and bacteriophage-mediated horizontal gene transfer (HGT). The pattern of evolution of E. coli in aging mice is characterized by a larger number of transposable element insertions and intergenic mutations compared to that in young mice, which is consistent with the gut of aging hosts harboring a stressful and iron limiting environment.

Divergent Biosynthesis of C-Nucleoside Minimycin and Indigoidine in Bacteria

Minimycin (MIN) is a C-nucleoside antibiotic structurally related to pseudouridine, and indigoidine is a naturally occurring blue pigment produced by diverse bacteria. Although MIN and indigoidine have been known for decades, the logic underlying the divergent biosynthesis of these interesting molecules has been obscure. Here, we report the identification of a minimal 5-gene cluster (min) essential for MIN biosynthesis. We demonstrated that a non-ribosomal peptide synthetase (MinA) governs “the switch” for the divergent biosynthesis of MIN and the cryptic indigoidine. We also demonstrated that MinC_N (the N-terminal phosphatase domain of MinC), MinD (uracil phosphoribosyltransferase), and MinT (transporter) function together as the safeguard enzymes, which collaboratively constitute an unusual self-resistance system. Finally, we provided evidence that MinD, utilizing an unprecedented substrate-competition strategy for self-resistance of the producer cell, maintains competition advantage over the active molecule MIN-5′-monophosphate by increasing the UMP pool in vivo. These findings greatly expand our knowledge regarding natural product biosynthesis.

•

A minimal 5-gene cluster (min) is essential for minimycin biosynthesis

•

Divergent biosynthesis of minimycin and indigoidine is mediated by an NRPS enzyme

•

A cascade of three safeguard enzymes constitutes the unusual self-resistance system

•

MinD functions as the key safeguard enzyme by increasing the UMP pool in vivo

Identification of a 2'-O-Methyluridine Nucleoside Hydrolase Using the Metagenomic Libraries

Ribose methylation is among the most ubiquitous modifications found in RNA. 2'-O-methyluridine is found in rRNA, snRNA, snoRNA and tRNA of Archaea, Bacteria, and Eukaryota. Moreover, 2'-O-methylribonucleosides are promising starting materials for the production of nucleic acid-based drugs. Despite the countless possibilities of practical use for the metabolic enzymes associated with methylated nucleosides, there are very few reports regarding the metabolic fate and enzymes involved in the metabolism of 2'-O-alkyl nucleosides. The presented work focuses on the cellular degradation of 2'-O-methyluridine. A novel enzyme was found using a screening strategy that employs Escherichia coli uracil auxotroph and the metagenomic libraries. A 2'-O-methyluridine hydrolase (RK9NH) has been identified together with an aldolase (RK9DPA)-forming a part of a probable gene cluster that is involved in the degradation of 2'-O-methylated nucleosides. The RK9NH is functional in E. coli uracil auxotroph and in vitro. The RK9NH nucleoside hydrolase could be engineered to enzymatically produce 2'-O-methylated nucleosides that are of great demand as raw materials for production of nucleic acid-based drugs. Moreover, RK9NH nucleoside hydrolase converts 5-fluorouridine, 5-fluoro-2'-deoxyuridine and 5-fluoro-2'-O-methyluridine into 5-fluorouracil, which suggests it could be employed in cancer therapy.

Measurement of Kinetics and Active Site Distances in Metalloenzymes Using Paramagnetic NMR with 13C Hyperpolarization

Paramagnetic relaxation enhancement (PRE) conjoint with hyperpolarized NMR reveals structural information on the enzyme-product complex in an ongoing metalloenzyme-catalyzed reaction. Substrates of pseudouridine monophosphate glycosidase are hyperpolarized using the dynamic nuclear polarization (DNP) method. Time series of ¹³C NMR spectra are subsequently measured with the enzyme containing diamagnetic Mg²⁺ or paramagnetic Mn²⁺ ions in the active site. The differences of the signal evolution and line widths in the Mg²⁺ vs Mn²⁺ reactions are explained through PRE in the enzyme-bound product, which is in fast exchange with its free form. Here, a strong distance dependence of the paramagnetically enhanced relaxation rates enables the calculation of distances from product atoms to the metal center in the complexed structure. The same method can be used to add structural information to real-time characterizations of chemical processes involving compounds with naturally present or artificially introduced paramagnetic sites.

A gene encoding a DUF523 domain protein is involved in the conversion of 2-thiouracil into uracil

Modified nucleotides are present in many RNA species in all Domains of Life. While the biosynthetic pathways of such nucleotides are well studied, much less is known about the degradation of RNAs and the return to the metabolism of modified nucleotides, their respective nucleosides or heterocyclic bases. Using an E. coli uracil auxotroph, we screened the metagenomic libraries for genes, which would allow the conversion of 2-thiouracil to uracil and thereby lead to the growth on a defined synthetic medium. We show that a gene encoding a protein consisting of previously uncharacterized Domain of Unknown Function 523 (DUF523) is responsible for such phenotype. We have purified this recombinant protein and demonstrated that it contains a FeS cluster. The substitution of cysteines, which have been predicted to form such clusters, with alanines abolished the growth phenotype. We conclude that DUF523 is involved in the conversion of 2-thiouracil into uracil in vivo.

Proteome of Plant Peroxisomes

Plant peroxisomes are required for a number of fundamental physiological processes, such as primary and secondary metabolism, development and stress response. Indexing the dynamic peroxisome proteome is prerequisite to fully understanding the importance of these organelles. Mass Spectrometry (MS)-based proteome analysis has allowed the identification of novel peroxisomal proteins and pathways in a relatively high-throughput fashion and significantly expanded the list of proteins and biochemical reactions in plant peroxisomes. In this chapter, we summarize the experimental proteomic studies performed in plants, compile a list of ~200 confirmed Arabidopsis peroxisomal proteins, and discuss the diverse plant peroxisome functions with an emphasis on the role of Arabidopsis MS-based proteomics in discovering new peroxisome functions. Many plant peroxisome proteins and biochemical pathways are specific to plants, substantiating the complexity, plasticity and uniqueness of plant peroxisomes. Mapping the full plant peroxisome proteome will provide a knowledge base for the improvement of crop production, quality and stress tolerance.

Molecular modeling studies of pseudouridine isoxazolidinyl nucleoside analogues as potential inhibitors of the pseudouridine 5'-monophosphate glycosidase

In this paper, we investigated the hypothesis that pseudouridine isoxazolidinyl nucleoside analogues could act as potential inhibitors of the pseudouridine 5'-monophosphate glycosidase. This purpose was pursued using molecular modeling and in silico ADME-Tox profiling. From these studies emerged that the isoxazolidinyl derivative 1 5'-monophosphate can be effectively accommodated within the active site of the enzyme with a ligand efficiency higher than that of the natural substrate. In this context, the poor nucleofugality of the N-protonated isoxazolidine prevents or slows down, the first mechanistic step proposed for the degradation of the pseudouridine 5'-monophosphate glycosidase, leading to the enzyme inhibition. Finally, the results of the physicochemical and ADME-Tox informative analysis pointed out that compound 1 is weakly bounded to plasma protein, only moderately permeate the blood-brain barrier, and is non-carcinogen in rat and mouse. To the best of our knowledge, this is the first paper that introduces the possibility of inhibition of pseudouridine 5'-monophosphate glycosidase by a molecule that competing with the natural substrate hinders the glycosidic C-C bond cleavage.

Evolving insights into RNA modifications and their functional diversity in the brain

In this Perspective, we expand the notion of temporal regulation of RNA in the brain and propose that the qualitative nature of RNA and its metabolism, together with RNA abundance, are essential for the molecular mechanisms underlying experience-dependent plasticity. We discuss emerging concepts in the newly burgeoning field of epitranscriptomics, which are predicted to be heavily involved in cognitive function. These include activity-induced RNA modifications, RNA editing, dynamic changes in the secondary structure of RNA, and RNA localization. Each is described with an emphasis on its role in regulating the function of both protein-coding genes, as well as various noncoding regulatory RNAs, and how each might influence learning and memory.

Discovery of the Showdomycin Gene Cluster from Streptomyces showdoensis ATCC 15227 Yields Insight into the Biosynthetic Logic of C-Nucleoside Antibiotics

Nucleoside antibiotics are a large class of pharmaceutically relevant chemical entities, which exhibit a broad spectrum of biological activities. Most nucleosides belong to the canonical N-nucleoside family, where the heterocyclic unit is connected to the carbohydrate through a carbon-nitrogen bond. However, atypical C-nucleosides were isolated from Streptomyces bacteria over 50 years ago, but the molecular basis for formation of these metabolites has been unknown. Here, we have sequenced the genome of S. showdoensis ATCC 15227 and identified the gene cluster responsible for showdomycin production. Key to the detection was the presence of sdmA, encoding an enzyme of the pseudouridine monophosphate glycosidase family, which could catalyze formation of the C-glycosidic bond. Sequence analysis revealed an unusual combination of biosynthetic genes, while inactivation and subsequent complementation of sdmA confirmed the involvement of the locus in showdomycin formation. The study provides the first steps toward generation of novel C-nucleosides by pathway engineering.

C-Nucleosides To Be Revisited

Two new C-nucleoside analogues, BCX4430, an imino-C-nucleoside, and GS-6620, a phosphoramidate derivative of 1'-cyano-2'-C-methyl-4-aza-7,9-dideazaadenosine C-nucleoside, have been recently described as effective against filovirus infections (Marburg) and hepatitis C virus (HCV), respectively. The first C-nucleoside analogues were described about half a century ago. The C-nucleoside pseudouridine is a natural component of RNA, and various other C-nucleoside analogues have been reported previously for their antiviral and/or anticancer potential, the most prominent being pyrazofurin, tiazofurin, and selenazofurin. In the meantime, showdomycin, formycin, and various triazole, pyrazine, pyridine, dihydroxyphenyl, thienopyrimidine, pyrazolotriazine, and porphyrin C-nucleoside analogues have been described. It would be worth revisiting these C-nucleosides and derivatives thereof, including their phosphoramidates, for their therapeutic potential in the treatment of virus infections and, where appropriate, cancer as well.

Functional Diversity of Haloacid Dehalogenase Superfamily Phosphatases from Saccharomyces cerevisiae: BIOCHEMICAL, STRUCTURAL, AND EVOLUTIONARY INSIGHTS

The haloacid dehalogenase (HAD)-like enzymes comprise a large superfamily of phosphohydrolases present in all organisms. The Saccharomyces cerevisiae genome encodes at least 19 soluble HADs, including 10 uncharacterized proteins. Here, we biochemically characterized 13 yeast phosphatases from the HAD superfamily, which includes both specific and promiscuous enzymes active against various phosphorylated metabolites and peptides with several HADs implicated in detoxification of phosphorylated compounds and pseudouridine. The crystal structures of four yeast HADs provided insight into their active sites, whereas the structure of the YKR070W dimer in complex with substrate revealed a composite substrate-binding site. Although the S. cerevisiae and Escherichia coli HADs share low sequence similarities, the comparison of their substrate profiles revealed seven phosphatases with common preferred substrates. The cluster of secondary substrates supporting significant activity of both S. cerevisiae and E. coli HADs includes 28 common metabolites that appear to represent the pool of potential activities for the evolution of novel HAD phosphatases. Evolution of novel substrate specificities of HAD phosphatases shows no strict correlation with sequence divergence. Thus, evolution of the HAD superfamily combines the conservation of the overall substrate pool and the substrate profiles of some enzymes with remarkable biochemical and structural flexibility of other superfamily members.

In vitro and in vivo models of Staphylococcus aureus endophthalmitis implicate specific nutrients in ocular infection

PURPOSE

To define global transcriptional responses of Staphylococcus aureus and its codY mutant (CodY is a transcription regulator of virulence and metabolic genes in response to branched-chain amino acids) when growing in bovine aqueous (AH) and vitreous humor (VH) in vitro, and to investigate the impact of codY deletion on S. aureus virulence in a novel murine anterior chamber (AC) infection model.

METHODS

For the in vitro model, differential transcriptomic gene expression of S. aureus and its codY mutant grown in chemically defined medium (CDM), AH, and VH was analyzed. Furthermore, the strains were inoculated into the AC of mice. Changes in bacterial growth, electroretinography and inflammation scores were monitored.

RESULTS

Bovine AH and VH provide sufficient nutrition for S. aureus growth in vitro. Transcriptome analysis identified 72 unique open reading frames differentially regulated ≥10-fold between CDM, AH, and VH. In the AC model, we found comparable growth of the codY mutant and wild type strains in vivo. Average inflammation scores and retinal function were significantly worse for codY mutant-infected eyes at 24 h post-infection.

CONCLUSION

Our in vitro bovine AH and VH models identified likely nutrient sources for S. aureus in the ocular milieu. The in vivo model suggests that control of branched-chain amino acid availability has therapeutic potential in limiting S. aureus endophthalmitis severity.

Genome-wide metabolic re-annotation of Ashbya gossypii: new insights into its metabolism through a comparative analysis with Saccharomyces cerevisiae and Kluyveromyces lactis

BACKGROUND

Ashbya gossypii is an industrially relevant microorganism traditionally used for riboflavin production. Despite the high gene homology and gene order conservation comparatively with Saccharomyces cerevisiae, it presents a lower level of genomic complexity. Its type of growth, placing it among filamentous fungi, questions how close it really is from the budding yeast, namely in terms of metabolism, therefore raising the need for an extensive and thorough study of its entire metabolism. This work reports the first manual enzymatic genome-wide re-annotation of A. gossypii as well as the first annotation of membrane transport proteins.

RESULTS

After applying a developed enzymatic re-annotation pipeline, 847 genes were assigned with metabolic functions. Comparatively to KEGG's annotation, these data corrected the function for 14% of the common genes and increased the information for 52 genes, either completing existing partial EC numbers or adding new ones. Furthermore, 22 unreported enzymatic functions were found, corresponding to a significant increase in the knowledge of the metabolism of this organism. The information retrieved from the metabolic re-annotation and transport annotation was used for a comprehensive analysis of A. gossypii's metabolism in comparison to the one of S. cerevisiae (post-WGD - whole genome duplication) and Kluyveromyces lactis (pre-WGD), suggesting some relevant differences in several parts of their metabolism, with the majority being found for the metabolism of purines, pyrimidines, nitrogen and lipids. A considerable number of enzymes were found exclusively in A. gossypii comparatively with K. lactis (90) and S. cerevisiae (13). In a similar way, 176 and 123 enzymatic functions were absent on A. gossypii comparatively to K. lactis and S. cerevisiae, respectively, confirming some of the well-known phenotypes of this organism.

CONCLUSIONS

This high quality metabolic re-annotation, together with the first membrane transporters annotation and the metabolic comparative analysis, represents a new important tool for the study and better understanding of A. gossypii's metabolism.

Pseudouridine: still mysterious, but never a fake (uridine)!

Pseudouridine (Ψ) is the most abundant of >150 nucleoside modifications in RNA. Although Ψ was discovered as the first modified nucleoside more than half a century ago, neither the enzymatic mechanism of its formation, nor the function of this modification are fully elucidated. We present the consistent picture of Ψ synthases, their substrates and their substrate positions in model organisms of all domains of life as it has emerged to date and point out the challenges that remain concerning higher eukaryotes and the elucidation of the enzymatic mechanism.

Chemoenzymatic synthesis of novel C-ribosylated naphthoquinones

The biological activity of many natural products is dependent on the presence of carbohydrate units, which are usually attached via an O-glycosidic linkage by glycosyltransferases. Recently, an exceptional C-ribosylation event was discovered in the biosynthesis of the polyketide antibiotic alnumycin A. The two-step process involves initial attachment of d-ribose-5-phosphate to the polyaromatic aglycone by the C-glycosynthase AlnA and subsequent dephosphorylation by AlnB, an enzyme of the haloacid dehalogenase family. Here, we tested 23 unnatural substrates to probe the C-ribosylation reaction. The chemoenzymatic synthesis of C-ribosylated juglone, 7-methyl juglone, monomethyl naphthazarin, 8-chloro-7-methyl juglone, and 9-hydroxy-1,4-anthraquinone revealed the importance of a 1,4-quinoid system with an adjacent phenolic ring in order for reaction to occur. To further rationalize the molecular basis for reactivity, factors governing substrate recognition were investigated by NMR binding experiments. Additionally, the suitability of substrates for nucleophilic substitution was assessed by molecular modeling using density functional theory (DFT) calculations.

Structural basis for C-ribosylation in the alnumycin A biosynthetic pathway

Alnumycin A is an exceptional aromatic polyketide that contains a carbohydrate-like 4'-hydroxy-5'-hydroxymethyl-2',7'-dioxane moiety attached to the aglycone via a carbon-carbon bond. Recently, we have identified the D-ribose-5-phosphate origin of the dioxane unit and demonstrated that AlnA and AlnB are responsible for the overall C-ribosylation reaction. Here, we provide direct evidence that AlnA is a natural C-glycosynthase, which catalyzes the attachment of D-ribose-5-phosphate to prealnumycin by formation of the C(8)-C(1') bond as demonstrated by the structure of the intermediate alnumycin P. This compound is subsequently dephosphorylated by AlnB, an enzyme of the haloacid dehalogenase superfamily. Structure determination of the native trimeric AlnA to 2.1-Å resolution revealed a highly globular fold encompassing an α/β/α sandwich. The crystal structure of the complex with D-ribose-5-phosphate indicated that the phosphosugar is bound in the open-chain configuration. Identification of residues E29, K86, and K159 near the C-1 carbonyl of the ligand led us to propose that the carbon-carbon bond formation proceeds through a Michael-type addition. Determination of the crystal structure of the monomeric AlnB in the open conformation to 1.25-Å resolution showed that the protein consists of core and cap domains. Modeling of alnumycin P inside the cap domain positioned the phosphate group next to a Mg(2+) ion present at the junction of the domains. Mutagenesis data were consistent with the canonical reaction mechanism for this enzyme family revealing the importance of residues D15 and D17 for catalysis. The characterization of the prealnumycin C-ribosylation illustrates an alternative means for attachment of carbohydrates to natural products.

Pseudouridine monophosphate glycosidase: a new glycosidase mechanism

Pseudouridine (Ψ), the most abundant modification in RNA, is synthesized in situ using Ψ synthase. Recently, a pathway for the degradation of Ψ was described [Preumont, A., Snoussi, K., Stroobant, V., Collet, J. F., and Van Schaftingen, E. (2008) J. Biol. Chem. 283, 25238-25246]. In this pathway, Ψ is first converted to Ψ 5'-monophosphate (ΨMP) by Ψ kinase and then ΨMP is degraded by ΨMP glycosidase to uracil and ribose 5-phosphate. ΨMP glycosidase is the first example of a mechanistically characterized enzyme that cleaves a C-C glycosidic bond. Here we report X-ray crystal structures of Escherichia coli ΨMP glycosidase and a complex of the K166A mutant with ΨMP. We also report the structures of a ring-opened ribose 5-phosphate adduct and a ring-opened ribose ΨMP adduct. These structures provide four snapshots along the reaction coordinate. The structural studies suggested that the reaction utilizes a Lys166 adduct during catalysis. Biochemical and mass spectrometry data further confirmed the existence of a lysine adduct. We used site-directed mutagenesis combined with kinetic analysis to identify roles for specific active site residues. Together, these data suggest that ΨMP glycosidase catalyzes the cleavage of the C-C glycosidic bond through a novel ribose ring-opening mechanism.

Biosynthetic pathway toward carbohydrate-like moieties of alnumycins contains unusual steps for C-C bond formation and cleavage

Carbohydrate moieties are important components of natural products, which are often imperative for the solubility and biological activity of the compounds. The aromatic polyketide alnumycin A contains an extraordinary sugar-like 4'-hydroxy-5'-hydroxymethyl-2',7'-dioxane moiety attached via a carbon-carbon bond to the aglycone. Here we have extensively investigated the biosynthesis of the dioxane unit through (13)C labeling studies, gene inactivation experiments and enzymatic synthesis. We show that AlnA and AlnB, members of the pseudouridine glycosidase and haloacid dehalogenase enzyme families, respectively, catalyze C-ribosylation conceivably through Michael-type addition of d-ribose-5-phosphate and dephosphorylation. The ribose moiety may be attached both in furanose (alnumycin C) and pyranose (alnumycin D) forms. The C(1')-C(2') bond of alnumycin C is subsequently cleaved and the ribose unit is rearranged into an unprecedented dioxolane (cis-bicyclo[3.3.0]-2',4',6'-trioxaoctan-3'β-ol) structure present in alnumycin B. The reaction is catalyzed by Aln6, which belongs to a previously uncharacterized enzyme family. The conversion was accompanied with consumption of O(2) and formation of H(2)O(2), which allowed us to propose that the reaction may proceed via hydroxylation of C1' followed by retro-aldol cleavage and acetal formation. Interestingly, no cofactors could be detected and the reaction was also conducted in the presence of metal chelating agents. The last step is the conversion of alnumycin B into the final end-product alnumycin A catalyzed by Aln4, an NADPH-dependent aldo-keto reductase. This characterization of the dioxane biosynthetic pathway sets the basis for the utilization of C-C bound ribose, dioxolane and dioxane moieties in the generation of improved biologically active compounds.

Mining high-throughput experimental data to link gene and function

Nearly 2200 genomes that encode around 6 million proteins have now been sequenced. Around 40% of these proteins are of unknown function, even when function is loosely and minimally defined as 'belonging to a superfamily'. In addition to in silico methods, the swelling stream of high-throughput experimental data can give valuable clues for linking these unknowns with precise biological roles. The goal is to develop integrative data-mining platforms that allow the scientific community at large to access and utilize this rich source of experimental knowledge. To this end, we review recent advances in generating whole-genome experimental datasets, where this data can be accessed, and how it can be used to drive prediction of gene function.

Synergistic use of plant-prokaryote comparative genomics for functional annotations

BACKGROUND

Identifying functions for all gene products in all sequenced organisms is a central challenge of the post-genomic era. However, at least 30-50% of the proteins encoded by any given genome are of unknown or vaguely known function, and a large number are wrongly annotated. Many of these 'unknown' proteins are common to prokaryotes and plants. We set out to predict and experimentally test the functions of such proteins. Our approach to functional prediction integrates comparative genomics based mainly on microbial genomes with functional genomic data from model microorganisms and post-genomic data from plants. This approach bridges the gap between automated homology-based annotations and the classical gene discovery efforts of experimentalists, and is more powerful than purely computational approaches to identifying gene-function associations.

RESULTS

Among Arabidopsis genes, we focused on those (2,325 in total) that (i) are unique or belong to families with no more than three members, (ii) occur in prokaryotes, and (iii) have unknown or poorly known functions. Computer-assisted selection of promising targets for deeper analysis was based on homology-independent characteristics associated in the SEED database with the prokaryotic members of each family. In-depth comparative genomic analysis was performed for 360 top candidate families. From this pool, 78 families were connected to general areas of metabolism and, of these families, specific functional predictions were made for 41. Twenty-one predicted functions have been experimentally tested or are currently under investigation by our group in at least one prokaryotic organism (nine of them have been validated, four invalidated, and eight are in progress). Ten additional predictions have been independently validated by other groups. Discovering the function of very widespread but hitherto enigmatic proteins such as the YrdC or YgfZ families illustrates the power of our approach.

CONCLUSIONS

Our approach correctly predicted functions for 19 uncharacterized protein families from plants and prokaryotes; none of these functions had previously been correctly predicted by computational methods. The resulting annotations could be propagated with confidence to over six thousand homologous proteins encoded in over 900 bacterial, archaeal, and eukaryotic genomes currently available in public databases.

Toward a definition of the complete proteome of plant peroxisomes: Where experimental proteomics must be complemented by bioinformatics

In the past few years, proteome analysis of Arabidopsis peroxisomes has been established by the complementary efforts of four research groups and has emerged as the major unbiased approach to identify new peroxisomal proteins on a large scale. Collectively, more than 100 new candidate proteins from plant peroxisomes have been identified, including long-awaited low-abundance proteins. More than 50 proteins have been validated as peroxisome targeted, nearly doubling the number of established plant peroxisomal proteins. Sequence homologies of the new proteins predict unexpected enzyme activities, novel metabolic pathways and unknown non-metabolic peroxisome functions. Despite this remarkable success, proteome analyses of plant peroxisomes remain highly material intensive and require major preparative efforts. Characterization of the membrane proteome or post-translational protein modifications poses major technical challenges. New strategies, including quantitative mass spectrometry methods, need to be applied to allow further identifications of plant peroxisomal proteins, such as of stress-inducible proteins. In the long process of defining the complete proteome of plant peroxisomes, the prediction of peroxisome-targeted proteins from plant genome sequences emerges as an essential complementary approach to identify additional peroxisomal proteins that are, for instance, specific to peroxisome variants from minor tissues and organs or to abiotically stressed model and crop plants.

A new approach to assess and predict the functional roles of proteins across all known structures

The three dimensional atomic structures of proteins provide information regarding their function; and codified relationships between structure and function enable the assessment of function from structure. In the current study, a new data mining tool was implemented that checks current gene ontology (GO) annotations and predicts new ones across all the protein structures available in the Protein Data Bank (PDB). The tool overcomes some of the challenges of utilizing large amounts of protein annotation and measurement information to form correspondences between protein structure and function. Protein attributes were extracted from the Structural Biology Knowledgebase and open source biological databases. Based on the presence or absence of a given set of attributes, a given protein's functional annotations were inferred. The results show that attributes derived from the three dimensional structures of proteins enhanced predictions over that using attributes only derived from primary amino acid sequence. Some predictions reflected known but not completely documented GO annotations. For example, predictions for the GO term for copper ion binding reflected used information a copper ion was known to interact with the protein based on information in a ligand interaction database. Other predictions were novel and require further experimental validation. These include predictions for proteins labeled as unknown function in the PDB. Two examples are a role in the regulation of transcription for the protein AF1396 from Archaeoglobus fulgidus and a role in RNA metabolism for the protein psuG from Thermotoga maritima.

HDHD1, which is often deleted in X-linked ichthyosis, encodes a pseudouridine-5'-phosphatase

Pseudouridine, the fifth-most abundant nucleoside in RNA, is not metabolized in mammals, but is excreted intact in urine. The purpose of the present work was to search for an enzyme that would dephosphorylate pseudouridine 5'-phosphate, a potential intermediate in RNA degradation. We show that human erythrocytes contain a pseudouridine-5'-phosphatase displaying a Km ≤ 1 μM for its substrate. The activity of the partially purified enzyme was dependent on Mg2+, and was inhibited by Ca2+ and vanadate, suggesting that it belonged to the 'haloacid dehalogenase' family of phosphatases. Its low molecular mass (26 kDa) suggested that this phosphatase could correspond to the protein encoded by the HDHD1 (haloacid dehalogenase-like hydrolase domain-containing 1) gene, present next to the STS (steroid sulfatase) gene on human chromosome Xp22. Purified human recombinant HDHD1 dephosphorylated pseudouridine 5'-phosphate with a kcat of 1.6 s-1, a Km of 0.3 μM and a catalytic efficiency at least 1000-fold higher than that on which it acted on other phosphate esters, including 5'-UMP. The molecular identity of pseudouridine-5'-phosphatase was confirmed by the finding that its activity was negligible (<10% of controls) in extracts of B-cell lymphoblasts or erythrocytes from X-linked ichthyosis patients harbouring a combined deletion of the STS gene (the X-linked ichthyosis gene) and the HDHD1 gene. Furthermore, pseudouridine-5'-phosphatase activity was 1.5-fold higher in erythrocytes from women compared with men, in agreement with the HDHD1 gene undergoing only partial inactivation in females. In conclusion, HDHD1 is a phosphatase specifically involved in dephosphorylation of a modified nucleotide present in RNA.

Host imprints on bacterial genomes--rapid, divergent evolution in individual patients

Bacteria lose or gain genetic material and through selection, new variants become fixed in the population. Here we provide the first, genome-wide example of a single bacterial strain's evolution in different deliberately colonized patients and the surprising insight that hosts appear to personalize their microflora. By first obtaining the complete genome sequence of the prototype asymptomatic bacteriuria strain E. coli 83972 and then resequencing its descendants after therapeutic bladder colonization of different patients, we identified 34 mutations, which affected metabolic and virulence-related genes. Further transcriptome and proteome analysis proved that these genome changes altered bacterial gene expression resulting in unique adaptation patterns in each patient. Our results provide evidence that, in addition to stochastic events, adaptive bacterial evolution is driven by individual host environments. Ongoing loss of gene function supports the hypothesis that evolution towards commensalism rather than virulence is favored during asymptomatic bladder colonization.