Structural bioinformatics analysis of enzymes involved in the biosynthesis pathway of the hypermodified nucleoside ms(2)io(6)A37 in tRNA

Landscape of Post-Transcriptional tRNA Modifications in Streptomyces albidoflavus J1074 as Portrayed by Mass Spectrometry and Genomic Data Mining

Actinobacterial genus Streptomyces (streptomycetes) represents one of the largest cultivable group of bacteria famous for their ability to produce valuable specialized (secondary) metabolites. Regulation of secondary metabolic pathways inextricably couples the latter to essential cellular processes that determine levels of amino acids, carbohydrates, phosphate, etc. Post-transcriptional tRNA modifications remain one of the least studied aspects of streptomycete physiology, albeit a few of them were recently shown to impact antibiotic production. In this study, we describe the diversity of post-transcriptional tRNA modifications in model strain Streptomyces albus (albidoflavus) J1074 by combining mass spectrometry and genomic data. Our results show that J1074 can produce more chemically distinct tRNA modifications than previously thought. An in silico approach identified orthologs for enzymes governing most of the identified tRNA modifications. Yet, genetic control of certain modifications remained elusive, suggesting early divergence of tRNA modification pathways in Streptomyces from the better studied model bacteria, such as Escherichia coli and Bacillus subtilis. As a first point in case, our data point to the presence of a non-canonical MiaE enzyme performing hydroxylation of prenylated adenosines. A further finding concerns the methylthiotransferase MiaB, which requires previous modification of adenosines by MiaA to i6A for thiomethylation to ms²i6A. We show here that the J1074 ortholog, when overexpressed, yields ms²A in a ΔmiaA background. Our results set the working ground for and justify a more detailed studies of biological significance of tRNA modification pathways in streptomycetes. IMPORTANCE Post-transcriptional tRNA modifications (PTTMs) play an important role in maturation and functionality of tRNAs. Little is known about tRNA modifications in the antibiotic-producing actinobacterial genus Streptomyces, even though peculiar tRNA-based regulatory mechanisms operate in this taxon. We provide a first detailed description of the chemical diversity of PTTMs in the model species, S. albidoflavus J1074, and identify most plausible genes for these PTTMs. Some of the PTTMs are described for the first time for Streptomyces. Production of certain PTTMs in J1074 appears to depend on enzymes that show no sequence similarity to known PTTM enzymes from model species. Our findings are of relevance for interrogation of genetic basis of PTTMs in pathogenic actinobacteria, such as M. tuberculosis.

A tRNA modifying enzyme as a tunable regulatory nexus for bacterial stress responses and virulence

Post-transcriptional modifications can impact the stability and functionality of many different classes of RNA molecules and are an especially important aspect of tRNA regulation. It is hypothesized that cells can orchestrate rapid responses to changing environmental conditions by adjusting the specific types and levels of tRNA modifications. We uncovered strong evidence in support of this tRNA global regulation hypothesis by examining effects of the well-conserved tRNA modifying enzyme MiaA in extraintestinal pathogenic Escherichia coli (ExPEC), a major cause of urinary tract and bloodstream infections. MiaA mediates the prenylation of adenosine-37 within tRNAs that decode UNN codons, and we found it to be crucial to the fitness and virulence of ExPEC. MiaA levels shifted in response to stress via a post-transcriptional mechanism, resulting in marked changes in the amounts of fully modified MiaA substrates. Both ablation and forced overproduction of MiaA stimulated translational frameshifting and profoundly altered the ExPEC proteome, with variable effects attributable to UNN content, changes in the catalytic activity of MiaA, or availability of metabolic precursors. Cumulatively, these data indicate that balanced input from MiaA is critical for optimizing cellular responses, with MiaA acting much like a rheostat that can be used to realign global protein expression patterns.

Biochemical and Structural Aspects of Cytokinin Biosynthesis and Degradation in Bacteria

It has been known for quite some time that cytokinins, hormones typical of plants, are also produced and metabolized in bacteria. Most bacteria can only form the tRNA-bound cytokinins, but there are examples of plant-associated bacteria, both pathogenic and beneficial, that actively synthesize cytokinins to interact with their host. Similar to plants, bacteria produce diverse cytokinin metabolites, employing corresponding metabolic pathways. The identification of genes encoding the enzymes involved in cytokinin biosynthesis and metabolism facilitated their detailed characterization based on both classical enzyme assays and structural approaches. This review summarizes the present knowledge on key enzymes involved in cytokinin biosynthesis, modifications, and degradation in bacteria, and discusses their catalytic properties in relation to the presence of specific amino acid residues and protein structure.

The Hulks and the Deadpools of the Cytokinin Universe: A Dual Strategy for Cytokinin Production, Translocation, and Signal Transduction

Cytokinins are plant hormones, derivatives of adenine with a side chain at the N⁶-position. They are involved in many physiological processes. While the metabolism of trans-zeatin and isopentenyladenine, which are considered to be highly active cytokinins, has been extensively studied, there are others with less obvious functions, such as cis-zeatin, dihydrozeatin, and aromatic cytokinins, which have been comparatively neglected. To help explain this duality, we present a novel hypothesis metaphorically comparing various cytokinin forms, enzymes of CK metabolism, and their signalling and transporter functions to the comics superheroes Hulk and Deadpool. Hulk is a powerful but short-lived creation, whilst Deadpool presents a more subtle and enduring force. With this dual framework in mind, this review compares different cytokinin metabolites, and their biosynthesis, translocation, and sensing to illustrate the different mechanisms behind the two CK strategies. This is put together and applied to a plant developmental scale and, beyond plants, to interactions with organisms of other kingdoms, to highlight where future study can benefit the understanding of plant fitness and productivity.

The Origins and Roles of Methylthiolated Cytokinins: Evidence From Among Life Kingdoms

Cytokinins (CKs) are a group of adenine-derived, small signaling molecules of crucial importance for growth and multiple developmental processes in plants. Biological roles of classical CKs: isopentenyladenine (iP), trans and cis isomers of zeatin (tZ, cZ), and dihydrozeatin, have been studied extensively and their functions are well defined in many aspects of plant physiology. In parallel, extensive knowledge exists for genes involved in tRNA modifications that lead to the production of tRNA-bound methylthiolated CKs, especially in bacterial and mammalian systems. However, not much is known about the origins, fates, and possible functions of the unbound methylthiolated CKs (2MeS-CKs) in biological systems. 2MeS-CKs are the free base or riboside derivatives of iP or Z-type CKs, modified by the addition of a thiol group (–SH) at position 2 of the adenine ring that is subsequently methylated. Based on the evidence to date, these distinctive CK conjugates are derived exclusively via the tRNA degradation pathway. This review summarizes the knowledge on the probable steps involved in the biosynthesis of unbound 2MeS-CKs across diverse kingdoms of life. Furthermore, it provides examples of CK profiles of organisms from which the presence of 2MeS-CKs have been detected and confirms a close association and balance between the production of classical CKs and 2MeS-CKs. Finally, it discusses available reports regarding the possible physiological functions of 2MeS-CKs in different biological systems.

Structural, biochemical and functional analyses of tRNA-monooxygenase enzyme MiaE from Pseudomonas putida provide insights into tRNA/MiaE interaction

MiaE (2-methylthio-N6-isopentenyl-adenosine37-tRNA monooxygenase) is a unique non-heme diiron enzyme that catalyzes the O2-dependent post-transcriptional allylic hydroxylation of a hypermodified nucleotide 2-methylthio-N6-isopentenyl-adenosine (ms2i6A37) at position 37 of selected tRNA molecules to produce 2-methylthio-N6-4-hydroxyisopentenyl-adenosine (ms2io6A37). Here, we report the in vivo activity, biochemical, spectroscopic characterization and X-ray crystal structure of MiaE from Pseudomonas putida. The investigation demonstrates that the putative pp-2188 gene encodes a MiaE enzyme. The structure shows that Pp-MiaE consists of a catalytic diiron(III) domain with a four alpha-helix bundle fold. A docking model of Pp-MiaE in complex with tRNA, combined with site directed mutagenesis and in vivo activity shed light on the importance of an additional linker region for substrate tRNA recognition. Finally, krypton-pressurized Pp-MiaE experiments, revealed the presence of defined O2 site along a conserved hydrophobic tunnel leading to the diiron active center.

Genetic analysis of Streptomyces albus J1074 mia mutants suggests complex relationships between post-transcriptional tRNAXXA modifications and physiological traits

Proteins MiaA and MiaB catalyze sequential isopentenylation and methylthiolation, respectively, of adenosine residue in 37th position of tRNA_XXA. The mia mutations were recently shown by us to affect secondary metabolism and morphology of Streptomyces. However, it remained unknown as to whether both or one of the aforementioned modifications is critical for colony development and antibiotic production. Here, we addressed this issue through analysis of Streptomyces albus J1074 strains carrying double miaAmiaB knockout or extra copy of miaB gene. The double mutant differed from wild-type and miaA-minus strains in severity of morphological defects, growth dynamics, and secondary metabolism. Introduction of extra copy of miaB gene into miaA mutant restored aerial mycelium formation to the latter on certain solid media. Hence, miaB gene might be involved in tRNA thiomethylation in the absence of miaA; either MiaA- or MiaB-mediated modification appears to be enough to support normal metabolic and morphological processes in Streptomyces.

Gene ssfg_01967 (miaB) for tRNA modification influences morphogenesis and moenomycin biosynthesis in Streptomyces ghanaensis ATCC14672

Streptomyces ghanaensis ATCC14672 is remarkable for its production of phosphoglycolipid compounds, moenomycins, which serve as a blueprint for the development of a novel class of antibiotics based on inhibition of peptidoglycan glycosyltransferases. Here we employed mariner transposon (Tn) mutagenesis to find new regulatory genes essential for moenomycin production. We generated a library of 3000 mutants which were screened for altered antibiotic activity. Our focus centred on a single mutant, HIM5, which accumulated lower amounts of moenomycin and was impaired in morphogenesis as compared to the parental strain. HIM5 carried the Tn insertion within gene ssfg_01967 for putative tRNA (N6-isopentenyl adenosine(37)-C2)-methylthiotransferase, or MiaB, and led to a reduced level of thiomethylation at position 37 in the anticodon of S. ghanaensis transfer ribonucleic acid (tRNA). It is likely that the mutant phenotype of HIM5 stems from the way in which ssfg_01967::Tn influences translation of the rare leucine codon UUA in several genes for moenomycin production and life cycle progression in S. ghanaensis. This is the first report showing that quantitative changes in tRNA modification status in Streptomyces have physiological consequences.

Diverse Mechanisms of Sulfur Decoration in Bacterial tRNA and Their Cellular Functions

Sulfur-containing transfer ribonucleic acids (tRNAs) are ubiquitous biomolecules found in all organisms that possess a variety of functions. For decades, their roles in processes such as translation, structural stability, and cellular protection have been elucidated and appreciated. These thionucleosides are found in all types of bacteria; however, their biosynthetic pathways are distinct among different groups of bacteria. Considering that many of the thio-tRNA biosynthetic enzymes are absent in Gram-positive bacteria, recent studies have addressed how sulfur trafficking is regulated in these prokaryotic species. Interestingly, a novel proposal has been given for interplay among thionucleosides and the biosynthesis of other thiocofactors, through participation of shared-enzyme intermediates, the functions of which are impacted by the availability of substrate as well as metabolic demand of thiocofactors. This review describes the occurrence of thio-modifications in bacterial tRNA and current methods for detection of these modifications that have enabled studies on the biosynthesis and functions of S-containing tRNA across bacteria. It provides insight into potential modes of regulation and potential evolutionary events responsible for divergence in sulfur metabolism among prokaryotes.

Biosynthesis of Sulfur-Containing tRNA Modifications: A Comparison of Bacterial, Archaeal, and Eukaryotic Pathways

Post-translational tRNA modifications have very broad diversity and are present in all domains of life. They are important for proper tRNA functions. In this review, we emphasize the recent advances on the biosynthesis of sulfur-containing tRNA nucleosides including the 2-thiouridine (s²U) derivatives, 4-thiouridine (s⁴U), 2-thiocytidine (s²C), and 2-methylthioadenosine (ms²A). Their biosynthetic pathways have two major types depending on the requirement of iron–sulfur (Fe–S) clusters. In all cases, the first step in bacteria and eukaryotes is to activate the sulfur atom of free l-cysteine by cysteine desulfurases, generating a persulfide (R-S-SH) group. In some archaea, a cysteine desulfurase is missing. The following steps of the bacterial s²U and s⁴U formation are Fe–S cluster independent, and the activated sulfur is transferred by persulfide-carrier proteins. By contrast, the biosynthesis of bacterial s²C and ms²A require Fe–S cluster dependent enzymes. A recent study shows that the archaeal s⁴U synthetase (ThiI) and the eukaryotic cytosolic 2-thiouridine synthetase (Ncs6) are Fe–S enzymes; this expands the role of Fe–S enzymes in tRNA thiolation to the Archaea and Eukarya domains. The detailed reaction mechanisms of Fe–S cluster depend s²U and s⁴U formation await further investigations.

Identification of 2-methylthio cyclic N6-threonylcarbamoyladenosine (ms2ct6A) as a novel RNA modification at position 37 of tRNAs

Transfer RNA modifications play pivotal roles in protein synthesis. N6-threonylcarbamoyladenosine (t6A) and its derivatives are modifications found at position 37, 3΄-adjacent to the anticodon, in tRNAs responsible for ANN codons. These modifications are universally conserved in all domains of life. t6A and its derivatives have pleiotropic functions in protein synthesis including aminoacylation, decoding and translocation. We previously discovered a cyclic form of t6A (ct6A) as a chemically labile derivative of t6A in tRNAs from bacteria, fungi, plants and protists. Here, we report 2-methylthio cyclic t6A (ms2ct6A), a novel derivative of ct6A found in tRNAs from Bacillus subtilis, plants and Trypanosoma brucei. In B. subtilis and T. brucei, ms2ct6A disappeared and remained to be ms2t6A and ct6A by depletion of tcdA and mtaB homologs, respectively, demonstrating that TcdA and MtaB are responsible for biogenesis of ms2ct6A.

CDK5RAP1 deficiency induces cell cycle arrest and apoptosis in human breast cancer cell line by the ROS/JNK signaling pathway

Cyclin-dependent kinase 5 regulatory subunit associated protein 1 (CDK5RAP1) is an enzyme which post-synthetically converts the RNA modification N6-iso-pentenyladenosine (i6A) into 2-methylthio-N6-isopentenyladenosine (ms2i6A). However, the interaction between CDK5RAP1 deficiency and cell apoptosis has not been studied. Breast cancer has long been a leading cause of mortality in the world. Therefore, in the present study, CDK5RAP1 deficiency in a human breast cancer cell line was investigated. CDK5RAP1 small interfering RNA (siRNA) and negative control siRNA were transfected into MCF-7 cells, and the cells were further incubated for 48 h. CDK5RAP1 deficiency suppressed tumor growth in MCF-7 cells and arrested the cells at G2/M phase. CDK5RAP1 deficiency also induced cell apoptosis and reactive oxygen species (ROS) generation. Furthermore, western blot analysis showed that the expression of phospho-c-Jun N-terminal kinase (p-JNK), p53, caspase-9 and caspase-3 were upregulated in CDK5RAP1-deficient MCF-7 cells. Pretreatment with N-acetyl-cysteine (NAC), the inhibitor of ROS, or with SP600125, the inhibitor of JNK, prevented the apoptosis and the high expression of p-JNK, p53, caspase-9 and caspase-3 in CDK5RAP1-deficient MCF-7 cells. Taken together, these data indicated that CDK5RAP1 deficiency induced cell cycle arrest and apoptosis in human breast cancer MCF-7 cells by the ROS/JNK signaling pathway. Our findings indicated a novel therapeutic strategy for cancer.

Steady-state kinetics and spectroscopic characterization of enzyme-tRNA interactions for the non-heme diiron tRNA-monooxygenase, MiaE

MiaE [2-methylthio-N(6)-isopentenyl-adenosine(37)-tRNA monooxygenase] isolated from Salmonella typhimurium is a unique non-heme diiron enzyme that catalyzes the O2-dependent post-transcriptional allylic hydroxylation of a hypermodified nucleotide (ms(2)i(6)A37) at position 37 of selected tRNA molecules to produce 2-methylthio-N(6)-(4-hydroxyisopentenyl)-adenosine(37). In this work, isopentenylated tRNA substrates for MiaE were produced from small RNA oligomers corresponding to the anticodon stem loop (ACSL) region of tRNA(Trp) using recombinant MiaA and dimethylallyl pyrophosphate. Steady-state rates for MiaE-catalyzed substrate hydroxylation were determined using recombinant ferredoxin (Fd) and ferredoxin reductase (FdR) to provide a catalytic electron transport chain (ETC) using NADPH as the sole electron source. As with previously reported peroxide-shunt assays, steady-state product formation retains nearly stoichiometric (>98%) E stereoselectivity. MiaE-catalyzed i(6)A-ACSL(Trp) hydroxylation follows Michaelis-Menten saturation kinetics with kcat, KM, and V/K determined to be 0.10 ± 0.01 s(-1), 9.1 ± 1.5 μM, and ∼11000 M(-1) s(-1), respectively. While vastly slower, MiaE-catalyzed hydroxylation of free i(6)A nucleoside could also be observed using the (Fd/FdR)-ETC assay. By comparison to the V/K determined for i(6)A-ACSL substrates, an ∼6000-fold increase in enzymatic efficiency is imparted by ACSL(Trp)-MiaE interactions. The impact of substrate tRNA-MiaE interactions on protein secondary structure and active site electronic configuration was investigated using circular dichroism, dual-mode X-band electron paramagnetic resonance, and Mössbauer spectroscopies. These studies demonstrate that binding of tRNA to MiaE induces a protein conformational change that influences the electronic structure of the diiron site analogous to what has been observed for various bacterial multicomponent diiron monooxygenases upon titration with their corresponding effector proteins. These observations suggest that substrate-enzyme interactions may play a pivotal role in modulating the reactivity of the MiaE diiron active site. Moreover, the simplified monomeric (α) protein configuration exhibited by MiaE provide an unparalleled opportunity to study the impact of protein-effector interactions on non-heme diiron site geometry and reactivity.

Transfer RNA Modification: Presence, Synthesis, and Function

Transfer RNA (tRNA) from all organisms on this planet contains modified nucleosides, which are derivatives of the four major nucleosides. tRNA from Escherichia coli/Salmonella enterica serovar Typhimurium contains 33 different modified nucleosides, which are all, except one (Queuosine [Q]), synthesized on an oligonucleotide precursor, which by specific enzymes later matures into tRNA. The structural genes for these enzymes are found in mono- and polycistronic operons, the latter of which have a complex transcription and translation pattern. The synthesis of the tRNA-modifying enzymes is not regulated similarly, and it is not coordinated to that of their substrate, the tRNA. The synthesis of some of them (e.g., several methylated derivatives) is catalyzed by one enzyme, which is position and base specific, whereas synthesis of some has a very complex biosynthetic pathway involving several enzymes (e.g., 2-thiouridines, N  6-cyclicthreonyladenosine [ct6A], and Q). Several of the modified nucleosides are essential for viability (e.g., lysidin, ct6A, 1-methylguanosine), whereas the deficiency of others induces severe growth defects. However, some have no or only a small effect on growth at laboratory conditions. Modified nucleosides that are present in the anticodon loop or stem have a fundamental influence on the efficiency of charging the tRNA, reading cognate codons, and preventing missense and frameshift errors. Those that are present in the body of the tRNA primarily have a stabilizing effect on the tRNA. Thus, the ubiquitous presence of these modified nucleosides plays a pivotal role in the function of the tRNA by their influence on the stability and activity of the tRNA.

Peroxide-shunt substrate-specificity for the Salmonella typhimurium O2-dependent tRNA modifying monooxygenase (MiaE)

Post-transcriptional modifications of tRNA are made to structurally diversify tRNA. These modifications alter noncovalent interactions within the ribosomal machinery, resulting in phenotypic changes related to cell metabolism, growth, and virulence. MiaE is a carboxylate bridged, nonheme diiron monooxygenase, which catalyzes the O2-dependent hydroxylation of a hypermodified-tRNA nucleoside at position 37 (2-methylthio-N(6)-isopentenyl-adenosine(37)-tRNA) [designated ms(2)i(6)A37]. In this work, recombinant MiaE was cloned from Salmonella typhimurium , purified to homogeneity, and characterized by UV-visible and dual-mode X-band EPR spectroscopy for comparison to other nonheme diiron enzymes. Additionally, three nucleoside substrate-surrogates (i(6)A, Cl(2)i(6)A, and ms(2)i(6)A) and their corresponding hydroxylated products (io(6)A, Cl(2)io(6)A, and ms(2)io(6)A) were synthesized to investigate the chemo- and stereospecificity of this enzyme. In the absence of the native electron transport chain, the peroxide-shunt was utilized to monitor the rate of substrate hydroxylation. Remarkably, regardless of the substrate (i(6)A, Cl(2)i(6)A, and ms(2)i(6)A) used in peroxide-shunt assays, hydroxylation of the terminal isopentenyl-C4-position was observed with >97% E-stereoselectivity. No other nonspecific hydroxylation products were observed in enzymatic assays. Steady-state kinetic experiments also demonstrate that the initial rate of MiaE hydroxylation is highly influenced by the substituent at the C2-position of the nucleoside base (v0/[E] for ms(2)i(6)A > i(6)A > Cl(2)i(6)A). Indeed, the >3-fold rate enhancement exhibited by MiaE for the hydroxylation of the free ms(2)i(6)A nucleoside relative to i(6)A is consistent with previous whole cell assays reporting the ms(2)io(6)A and io(6)A product distribution within native tRNA-substrates. This observation suggests that the nucleoside C2-substituent is a key point of interaction regulating MiaE substrate specificity.

Structure and function of noncanonical nucleobases

DNA and RNA contain, next to the four canonical nucleobases, a number of modified nucleosides that extend their chemical information content. RNA is particularly rich in modifications, which is obviously an adaptation to their highly complex and variable functions. In fact, the modified nucleosides and their chemical structures establish a second layer of information which is of central importance to the function of the RNA molecules. Also the chemical diversity of DNA is greater than originally thought. Next to the four canonical bases, the DNA of higher organisms contains a total of four epigenetic bases: m(5) dC, hm(5) dC, f(5) dC und ca(5) dC. While all cells of an organism contain the same genetic material, their vastly different function and properties inside complex higher organisms require the controlled silencing and activation of cell-type specific genes. The regulation of the underlying silencing and activation process requires an additional layer of epigenetic information, which is clearly linked to increased chemical diversity. This diversity is provided by the modified non-canonical nucleosides in both DNA and RNA.

The CDK5 repressor CDK5RAP1 is a methylthiotransferase acting on nuclear and mitochondrial RNA

The unusual cyclin-dependent protein kinase 5 (CDK5) was discovered based on its sequence homology to cell cycle regulating CDKs. CDK5 was found to be active in brain tissues, where it is not involved in cell cycle regulation but in the regulation of neuronal cell differentiation and neurocytoskeleton dynamics. An aberrant regulation of CDK5 leads to the development of various neurodegenerative diseases including Alzheimer’s disease. Although CDK5 is not regulated by cyclins, its activity does depend on the association with a protein activator and the presence or absence of further inhibitory factors. Recently, CDK5RAP1 was discovered to inhibit the active CDK5 kinase. Here, we show that CDK5RAP1 is a radical SAM enzyme, which postsynthetically converts the RNA modification N6-isopentenyladenosine (i⁶A) into 2-methylthio-N6-isopentenyladenosine (ms²i⁶A). This conversion is surprisingly not limited to mitochondrial tRNA, where the modification was known to exist. Instead, CDK5RAP1 introduces the modification also into nuclear RNA species establishing a link between postsynthetic kinase-based protein modification and postsynthetic RNA modification.

Protein structural domains: definition and prediction

Recognition and prediction of structural domains in proteins is an important part of structure and function prediction. This unit lists the range of tools available for domain prediction, and describes sequence and structural analysis tools that complement domain prediction methods. Also detailed are the basic domain prediction steps, along with suggested strategies for different protein sequences and potential pitfalls in domain boundary prediction. The difficult problem of domain orientation prediction is also discussed. All the resources necessary for domain boundary prediction are accessible via publicly available Web servers and databases and do not require computational expertise.

Evolution in an oncogenic bacterial species with extreme genome plasticity: Helicobacter pylori East Asian genomes

BACKGROUND

The genome of Helicobacter pylori, an oncogenic bacterium in the human stomach, rapidly evolves and shows wide geographical divergence. The high incidence of stomach cancer in East Asia might be related to bacterial genotype. We used newly developed comparative methods to follow the evolution of East Asian H. pylori genomes using 20 complete genome sequences from Japanese, Korean, Amerind, European, and West African strains.

RESULTS

A phylogenetic tree of concatenated well-defined core genes supported divergence of the East Asian lineage (hspEAsia; Japanese and Korean) from the European lineage ancestor, and then from the Amerind lineage ancestor. Phylogenetic profiling revealed a large difference in the repertoire of outer membrane proteins (including oipA, hopMN, babABC, sabAB and vacA-2) through gene loss, gain, and mutation. All known functions associated with molybdenum, a rare element essential to nearly all organisms that catalyzes two-electron-transfer oxidation-reduction reactions, appeared to be inactivated. Two pathways linking acetyl~CoA and acetate appeared intact in some Japanese strains. Phylogenetic analysis revealed greater divergence between the East Asian (hspEAsia) and the European (hpEurope) genomes in proteins in host interaction, specifically virulence factors (tipα), outer membrane proteins, and lipopolysaccharide synthesis (human Lewis antigen mimicry) enzymes. Divergence was also seen in proteins in electron transfer and translation fidelity (miaA, tilS), a DNA recombinase/exonuclease that recognizes genome identity (addA), and DNA/RNA hybrid nucleases (rnhAB). Positively selected amino acid changes between hspEAsia and hpEurope were mapped to products of cagA, vacA, homC (outer membrane protein), sotB (sugar transport), and a translation fidelity factor (miaA). Large divergence was seen in genes related to antibiotics: frxA (metronidazole resistance), def (peptide deformylase, drug target), and ftsA (actin-like, drug target).

CONCLUSIONS

These results demonstrate dramatic genome evolution within a species, especially in likely host interaction genes. The East Asian strains appear to differ greatly from the European strains in electron transfer and redox reactions. These findings also suggest a model of adaptive evolution through proteome diversification and selection through modulation of translational fidelity. The results define H. pylori East Asian lineages and provide essential information for understanding their pathogenesis and designing drugs and therapies that target them.

Functional characterization of the YmcB and YqeV tRNA methylthiotransferases of Bacillus subtilis

Methylthiotransferases (MTTases) are a closely related family of proteins that perform both radical-S-adenosylmethionine (SAM) mediated sulfur insertion and SAM-dependent methylation to modify nucleic acid or protein targets with a methyl thioether group (-SCH(3)). Members of two of the four known subgroups of MTTases have been characterized, typified by MiaB, which modifies N(6)-isopentenyladenosine (i(6)A) to 2-methylthio-N(6)-isopentenyladenosine (ms(2)i(6)A) in tRNA, and RimO, which modifies a specific aspartate residue in ribosomal protein S12. In this work, we have characterized the two MTTases encoded by Bacillus subtilis 168 and find that, consistent with bioinformatic predictions, ymcB is required for ms(2)i(6)A formation (MiaB activity), and yqeV is required for modification of N(6)-threonylcarbamoyladenosine (t(6)A) to 2-methylthio-N(6)-threonylcarbamoyladenosine (ms(2)t(6)A) in tRNA. The enzyme responsible for the latter activity belongs to a third MTTase subgroup, no member of which has previously been characterized. We performed domain-swapping experiments between YmcB and YqeV to narrow down the protein domain(s) responsible for distinguishing i(6)A from t(6)A and found that the C-terminal TRAM domain, putatively involved with RNA binding, is likely not involved with this discrimination. Finally, we performed a computational analysis to identify candidate residues outside the TRAM domain that may be involved with substrate recognition. These residues represent interesting targets for further analysis.

Origin and evolution of peptide-modifying dioxygenases and identification of the wybutosine hydroxylase/hydroperoxidase

Unlike classical 2-oxoglutarate and iron-dependent dioxygenases, which include several nucleic acid modifiers, the structurally similar jumonji-related dioxygenase superfamily was only known to catalyze peptide modifications. Using comparative genomics methods, we predict that a family of jumonji-related enzymes catalyzes wybutosine hydroxylation/peroxidation at position 37 of eukaryotic tRNAPhe. Identification of this enzyme raised questions regarding the emergence of protein- and nucleic acid-modifying activities among jumonji-related domains. We addressed these with a natural classification of DSBH domains and reconstructed the precursor of the dioxygenases as a sugar-binding domain. This precursor gave rise to sugar epimerases and metal-binding sugar isomerases. The sugar isomerase active site was exapted for catalysis of oxygenation, with a radiation of these enzymes in bacteria, probably due to impetus from the primary oxygenation event in Earth’s history. 2-Oxoglutarate-dependent versions appear to have further expanded with rise of the tricarboxylic acid cycle. We identify previously under-appreciated aspects of their active site and multiple independent innovations of 2-oxoacid-binding basic residues among these superfamilies. We show that double-stranded β-helix dioxygenases diversified extensively in biosynthesis and modification of halogenated siderophores, antibiotics, peptide secondary metabolites and glycine-rich collagen-like proteins in bacteria. Jumonji-related domains diversified into three distinct lineages in bacterial secondary metabolism systems and these were precursors of the three major clades of eukaryotic enzymes. The specificity of wybutosine hydroxylase/peroxidase probably relates to the structural similarity of the modified moiety to the ancestral amino acid substrate of this superfamily.

Adenosyl radical: reagent and catalyst in enzyme reactions

Adenosine is undoubtedly an ancient biological molecule that is a component of many enzyme cofactors: ATP, FADH, NAD(P)H, and coenzyme A, to name but a few, and, of course, of RNA. Here we present an overview of the role of adenosine in its most reactive form: as an organic radical formed either by homolytic cleavage of adenosylcobalamin (coenzyme B(12), AdoCbl) or by single-electron reduction of S-adenosylmethionine (AdoMet) complexed to an iron-sulfur cluster. Although many of the enzymes we discuss are newly discovered, adenosine's role as a radical cofactor most likely arose very early in evolution, before the advent of photosynthesis and the production of molecular oxygen, which rapidly inactivates many radical enzymes. AdoCbl-dependent enzymes appear to be confined to a rather narrow repertoire of rearrangement reactions involving 1,2-hydrogen atom migrations; nevertheless, mechanistic insights gained from studying these enzymes have proved extremely valuable in understanding how enzymes generate and control highly reactive free radical intermediates. In contrast, there has been a recent explosion in the number of radical-AdoMet enzymes discovered that catalyze a remarkably wide range of chemically challenging reactions; here there is much still to learn about their mechanisms. Although all the radical-AdoMet enzymes so far characterized come from anaerobically growing microbes and are very oxygen sensitive, there is tantalizing evidence that some of these enzymes might be active in aerobic organisms including humans.

Deciphering synonymous codons in the three domains of life: co-evolution with specific tRNA modification enzymes

The strategies organisms use to decode synonymous codons in cytosolic protein synthesis are not uniform. The complete isoacceptor tRNA repertoire and the type of modified nucleoside found at the wobble position 34 of their anticodons were analyzed in all kingdoms of life. This led to the identification of four main decoding strategies that are diversely used in Bacteria, Archaea and Eukarya. Many of the modern tRNA modification enzymes acting at position 34 of tRNAs are present only in specific domains and obviously have arisen late during evolution. In an evolutionary fine-tuning process, these enzymes must have played an essential role in the progressive introduction of new amino acids, and in the refinement and standardization of the canonical nuclear genetic code observed in all extant organisms (functional convergent evolutionary hypothesis).

Assessment of domain boundary predictions and the prediction of intramolecular contacts in CASP8

This article details the assessment process and evaluation results for two categories in the 8th Critical Assessment of Protein Structure Prediction experiment (CASP8). The domain prediction category was evaluated with a range of scores including the Normalized Domain Overlap score and a domain boundary distance measure. Residue-residue contact predictions were evaluated with standard CASP measures, prediction accuracy, and Xd. In the domain boundary prediction category, prediction methods still make reliable predictions for targets that have structural templates, but continue to struggle to make good predictions for the few ab initio targets in CASP. There was little indication of improvement in the domain prediction category. The contact prediction category demonstrated that there was renewed interest among predictors and despite the small sample size the results suggested that there had been an increase in prediction accuracy. In contrast to CASP7 contact specialists predicted contacts more accurately than the majority of tertiary structure predictors. Despite this small success, the lack of free modeling targets makes it unlikely that either category will be included in their present form in CASP9.

Characterization of RimO, a new member of the methylthiotransferase subclass of the radical SAM superfamily

RimO, encoded by the yliG gene in Escherichia coli, has been recently identified in vivo as the enzyme responsible for the attachment of a methylthio group on the beta-carbon of Asp88 of the small ribosomal protein S12 [Anton, B. P., Saleh, L., Benner, J. S., Raleigh, E. A., Kasif, S., and Roberts, R. J. (2008) Proc. Natl. Acad. Sci. U.S.A. 105, 1826-1831]. To date, it is the only enzyme known to catalyze methylthiolation of a protein substrate; the four other naturally occurring methylthio modifications have been observed on tRNA. All members of the methylthiotransferase (MTTase) family, to which RimO belongs, have been shown to contain the canonical CxxxCxxC motif in their primary structures that is typical of the radical S-adenosylmethionine (SAM) family of proteins. MiaB, the only characterized MTTase, and the enzyme experimentally shown to be responsible for methylthiolation of N(6)-isopentenyladenosine of tRNA in E. coli and Thermotoga maritima, has been demonstrated to harbor two distinct [4Fe-4S] clusters. Herein, we report in vitro biochemical and spectroscopic characterization of RimO. We show by analytical and spectroscopic methods that RimO, overproduced in E. coli in the presence of iron-sulfur cluster biosynthesis proteins from Azotobacter vinelandii, contains one [4Fe-4S](2+) cluster. Reconstitution of this form of RimO (RimO(rcn)) with (57)Fe and sodium sulfide results in a protein that contains two [4Fe-4S](2+) clusters, similar to MiaB. We also show by mass spectrometry that RimO(rcn) catalyzes the attachment of a methylthio group to a peptide substrate analogue that mimics the loop structure bearing aspartyl 88 of the S12 ribosomal protein from E. coli. Kinetic analysis of this reaction shows that the activity of RimO(rcn) in the presence of the substrate analogue does not support a complete turnover. We discuss the possible requirement for an assembled ribosome for fully active RimO in vitro. Our findings are consistent with those of other enzymes that catalyze sulfur insertion, such as biotin synthase, lipoyl synthase, and MiaB.

Glutathione production by efficient ATP-regenerating Escherichia coli mutants

There is an ongoing demand to improve the ATP-regenerating system for industrial ATP-driven bioprocesses because of the low efficiency of ATP regeneration. To address this issue, we investigated the efficiency of ATP regeneration in Escherichia coli using the Permeable Cell Assay. This assay identified 40 single-gene deletion strains that had over 150% higher total cellular ATP synthetic activity relative to the parental strain. Most of them also showed higher ATP-driven glutathione synthesis. The deleted genes of the identified strains that showed increased efficiency of ATP regeneration for glutathione production could be divided into the following four groups: (1) glycolytic pathway-related genes, (2) genes related to degradation of ATP or adenosine, (3) global regulatory genes, and (4) genes whose contribution to the ATP regeneration is unknown. Furthermore, the high glutathione productivity of DeltanlpD, the highest glutathione-producing mutant strain, was due to its reduced sensitivity to the externally added ATP for ATP regeneration. This study showed that the Permeable Cell Assay was useful for improving the ATP-regenerating activity of E. coli for practical applications in various ATP-driven bioprocesses, much as that of glutathione production.

Type II restriction endonuclease R.Hpy188I belongs to the GIY-YIG nuclease superfamily, but exhibits an unusual active site

BACKGROUND

Catalytic domains of Type II restriction endonucleases (REases) belong to a few unrelated three-dimensional folds. While the PD-(D/E)XK fold is most common among these enzymes, crystal structures have been also determined for single representatives of two other folds: PLD (R.BfiI) and half-pipe (R.PabI). Bioinformatics analyses supported by mutagenesis experiments suggested that some REases belong to the HNH fold (e.g. R.KpnI), and that a small group represented by R.Eco29kI belongs to the GIY-YIG fold. However, for a large fraction of REases with known sequences, the three-dimensional fold and the architecture of the active site remain unknown, mostly due to extreme sequence divergence that hampers detection of homology to enzymes with known folds.

RESULTS

R.Hpy188I is a Type II REase with unknown structure. PSI-BLAST searches of the non-redundant protein sequence database reveal only 1 homolog (R.HpyF17I, with nearly identical amino acid sequence and the same DNA sequence specificity). Standard application of state-of-the-art protein fold-recognition methods failed to predict the relationship of R.Hpy188I to proteins with known structure or to other protein families. In order to increase the amount of evolutionary information in the multiple sequence alignment, we have expanded our sequence database searches to include sequences from metagenomics projects. This search resulted in identification of 23 further members of R.Hpy188I family, both from metagenomics and the non-redundant database. Moreover, fold-recognition analysis of the extended R.Hpy188I family revealed its relationship to the GIY-YIG domain and allowed for computational modeling of the R.Hpy188I structure. Analysis of the R.Hpy188I model in the light of sequence conservation among its homologs revealed an unusual variant of the active site, in which the typical Tyr residue of the YIG half-motif had been substituted by a Lys residue. Moreover, some of its homologs have the otherwise invariant Arg residue in a non-homologous position in sequence that nonetheless allows for spatial conservation of the guanidino group potentially involved in phosphate binding.

CONCLUSION

The present study eliminates a significant "white spot" on the structural map of REases. It also provides important insight into sequence-structure-function relationships in the GIY-YIG nuclease superfamily. Our results reveal that in the case of proteins with no or few detectable homologs in the standard "non-redundant" database, it is useful to expand this database by adding the metagenomic sequences, which may provide evolutionary linkage to detect more remote homologs.

MetaMQAP: a meta-server for the quality assessment of protein models

Background

Computational models of protein structure are usually inaccurate and exhibit significant deviations from the true structure. The utility of models depends on the degree of these deviations. A number of predictive methods have been developed to discriminate between the globally incorrect and approximately correct models. However, only a few methods predict correctness of different parts of computational models. Several Model Quality Assessment Programs (MQAPs) have been developed to detect local inaccuracies in unrefined crystallographic models, but it is not known if they are useful for computational models, which usually exhibit different and much more severe errors.

Results

The ability to identify local errors in models was tested for eight MQAPs: VERIFY3D, PROSA, BALA, ANOLEA, PROVE, TUNE, REFINER, PROQRES on 8251 models from the CASP-5 and CASP-6 experiments, by calculating the Spearman's rank correlation coefficients between per-residue scores of these methods and local deviations between C-alpha atoms in the models vs. experimental structures. As a reference, we calculated the value of correlation between the local deviations and trivial features that can be calculated for each residue directly from the models, i.e. solvent accessibility, depth in the structure, and the number of local and non-local neighbours. We found that absolute correlations of scores returned by the MQAPs and local deviations were poor for all methods. In addition, scores of PROQRES and several other MQAPs strongly correlate with 'trivial' features. Therefore, we developed MetaMQAP, a meta-predictor based on a multivariate regression model, which uses scores of the above-mentioned methods, but in which trivial parameters are controlled. MetaMQAP predicts the absolute deviation (in Ångströms) of individual C-alpha atoms between the model and the unknown true structure as well as global deviations (expressed as root mean square deviation and GDT_TS scores). Local model accuracy predicted by MetaMQAP shows an impressive correlation coefficient of 0.7 with true deviations from native structures, a significant improvement over all constituent primary MQAP scores. The global MetaMQAP score is correlated with model GDT_TS on the level of 0.89.

Conclusion

Finally, we compared our method with the MQAPs that scored best in the 7th edition of CASP, using CASP7 server models (not included in the MetaMQAP training set) as the test data. In our benchmark, MetaMQAP is outperformed only by PCONS6 and method QA_556 – methods that require comparison of multiple alternative models and score each of them depending on its similarity to other models. MetaMQAP is however the best among methods capable of evaluating just single models.

We implemented the MetaMQAP as a web server available for free use by all academic users at the URL

RimO, a MiaB-like enzyme, methylthiolates the universally conserved Asp88 residue of ribosomal protein S12 in Escherichia coli

Ribosomal protein S12 undergoes a unique posttranslational modification, methylthiolation of residue D88, in Escherichia coli and several other bacteria. Using mass spectrometry, we have identified the enzyme responsible for this modification in E. coli, the yliG gene product. This enzyme, which we propose be called RimO, is a radical-S-adenosylmethionine protein that bears strong sequence similarity to MiaB, which methylthiolates tRNA. We show that RimO and MiaB represent two of four subgroups of a larger, ancient family of likely methylthiotransferases, the other two of which are typified by Bacillus subtilis YqeV and Methanococcus jannaschii Mj0867, and we predict that RimO is unique among these subgroups in its modification of protein as opposed to tRNA. Despite this, RimO has not significantly diverged from the other three subgroups at the sequence level even within the C-terminal TRAM domain, which in the methyltransferase RumA is known to bind the RNA substrate and which we presume to be responsible for substrate binding and recognition in all four subgroups of methylthiotransferases. To our knowledge, RimO and MiaB represent the most extreme known case of resemblance between enzymes modifying protein and nucleic acid. The initial results presented here constitute a bioinformatics-driven prediction with preliminary experimental validation that should serve as the starting point for several interesting lines of further inquiry.