Unusual evolutionary history of the tRNA splicing endonuclease EndA: relationship to the LAGLIDADG and PD-(D/E)XK deoxyribonucleases

Integrated evolution of ribosomal RNAs, introns, and intron nurseries

The initial components of ribosomes first appeared more than 3.8 billion years ago during a time when many types of RNAs were evolving. While modern ribosomes are complex molecular machines consisting of rRNAs and proteins, they were assembled during early evolution by the association and joining of small functional RNA units. Introns may have provided the means to ligate many of these pieces together. All four classes of introns (group I, group II, spliceosomal, and archaeal) are present in many rRNA gene loci over a broad phylogenetic range. A survey of rRNA intron sequences across the three major life domains suggests that some of the classes of introns may have diverged from one another within rRNA gene loci. Analyses of rRNA sequences revealed self-splicing group I and group II introns are present in ancestral regions of the SSU (small subunit) and LSU (large subunit), whereas spliceosomal and archaeal introns appeared in sections of the rRNA that evolved later. Most classes of introns increased in number for approximately 1 billion years. However, their frequencies are low in the most recently evolved regions added to the SSU and LSU rRNAs. Furthermore, many of the introns appear to have been in the same locations for billions of years, suggesting an ancient origin for these sequences. In this Perspectives paper, I reviewed and analyzed rRNA intron sequences, locations, structural characteristics, and splicing mechanisms; and suggest that rRNA gene loci may have served as evolutionary nurseries for intron formation and diversification.

The CDI toxin of Yersinia kristensenii is a novel bacterial member of the RNase A superfamily

Contact-dependent growth inhibition (CDI) is an important mechanism of inter-bacterial competition found in many Gram-negative pathogens. CDI⁺ cells express cell-surface CdiA proteins that bind neighboring bacteria and deliver C-terminal toxin domains (CdiA-CT) to inhibit target-cell growth. CDI⁺ bacteria also produce CdiI immunity proteins, which specifically neutralize cognate CdiA-CT toxins to prevent self-inhibition. Here, we present the crystal structure of the CdiA-CT/CdiI^Ykris complex from Yersinia kristensenii ATCC 33638. CdiA-CT^Ykris adopts the same fold as angiogenin and other RNase A paralogs, but the toxin does not share sequence similarity with these nucleases and lacks the characteristic disulfide bonds of the superfamily. Consistent with the structural homology, CdiA-CT^Ykris has potent RNase activity in vitro and in vivo. Structure-guided mutagenesis reveals that His175, Arg186, Thr276 and Tyr278 contribute to CdiA-CT^Ykris activity, suggesting that these residues participate in substrate binding and/or catalysis. CdiI^Ykris binds directly over the putative active site and likely neutralizes toxicity by blocking access to RNA substrates. Significantly, CdiA-CT^Ykris is the first non-vertebrate protein found to possess the RNase A superfamily fold, and homologs of this toxin are associated with secretion systems in many Gram-negative and Gram-positive bacteria. These observations suggest that RNase A-like toxins are commonly deployed in inter-bacterial competition.

Functional Diversity of Cytotoxic tRNase/Immunity Protein Complexes from Burkholderia pseudomallei

Contact-dependent growth inhibition (CDI) is a widespread mechanism of inter-bacterial competition. CDI(+) bacteria deploy large CdiA effector proteins, which carry variable C-terminal toxin domains (CdiA-CT). CDI(+) cells also produce CdiI immunity proteins that specifically neutralize cognate CdiA-CT toxins to prevent auto-inhibition. Here, we present the crystal structure of the CdiA-CT/CdiI(E479) toxin/immunity protein complex from Burkholderia pseudomallei isolate E479. The CdiA-CT(E479) tRNase domain contains a core α/β-fold that is characteristic of PD(D/E)XK superfamily nucleases. Unexpectedly, the closest structural homolog of CdiA-CT(E479) is another CDI toxin domain from B. pseudomallei 1026b. Although unrelated in sequence, the two B. pseudomallei nuclease domains share similar folds and active-site architectures. By contrast, the CdiI(E479) and CdiI(1026b) immunity proteins share no significant sequence or structural homology. CdiA-CT(E479) and CdiA-CT(1026b) are both tRNases; however, each nuclease cleaves tRNA at a distinct position. We used a molecular docking approach to model each toxin bound to tRNA substrate. The resulting models fit into electron density envelopes generated by small-angle x-ray scattering analysis of catalytically inactive toxin domains bound stably to tRNA. CdiA-CT(E479) is the third CDI toxin found to have structural homology to the PD(D/E)XK superfamily. We propose that CDI systems exploit the inherent sequence variability and active-site plasticity of PD(D/E)XK nucleases to generate toxin diversity. These findings raise the possibility that many other uncharacterized CDI toxins may belong to the PD(D/E)XK superfamily.

Crystal structure of the Varkud satellite ribozyme

Varkud Satellite (VS) ribozyme mediates rolling circle replication of a plasmid found in the Neurospora mitochondria. We report crystal structures of this ribozyme at 3.1Å resolution, revealing an intertwined dimer formed by an exchange of substrate helices. Within each protomer, an arrangement of three-way helical junctions organizes seven helices into a global fold that creates a docking site for the substrate helix of the other protomer, resulting in the formation of two active sites in trans. This mode of RNA-RNA association resembles the process of domain swapping in proteins and has implications for RNA regulation and evolution. Within each active site, adenine and guanine nucleobases abut the scissile phosphate, poised to serve direct roles in catalysis. Similarities to the active sites of the hairpin and hammerhead ribozymes highlight the functional significance of active site features, underscore the ability of RNA to access functional architectures from distant regions of sequence space, and suggest convergent evolution.

Novel mutations in TSEN54 in pontocerebellar hypoplasia type 2

Pontocerebellar hypoplasias represent a group of neurodegenerative autosomal recessive disorders characterized by hypoplasia/atrophy of the cerebellum, hypoplastic ventral pons, and microcephaly and associated with various clinical features. Pontocerebellar hypolasia type 2 is the most common form, and different mutations in genes encoding subunits of the transfer ribonucleic acid (RNA)-splicing endonuclease (TSEN) complex were identified in patients. The authors report clinical, imaging, and molecular studies in 2 unrelated patients with different clinical pictures of the pontocerebellar hypoplasia type 2 spectrum and novel mutations in TSEN54, aiming to further define the clinical spectrum of the disease and possible indicators of more favorable progression. They identified a novel missense mutation c.355T>G/p.Y119D in compound heterozygosity with the "common" c.919G>T/p.A307S (patient 1) and a novel homozygous c.7ins6(CCGGAG)/p.E2-P3insPE variant (patient 2). An expanded array of mutations might contribute in defining possible differences in severity and phenotype-genotype correlations.

Phylogenomic analysis of the GIY-YIG nuclease superfamily

Background

The GIY-YIG domain was initially identified in homing endonucleases and later in other selfish mobile genetic elements (including restriction enzymes and non-LTR retrotransposons) and in enzymes involved in DNA repair and recombination. However, to date no systematic search for novel members of the GIY-YIG superfamily or comparative analysis of these enzymes has been reported.

Results

We carried out database searches to identify all members of known GIY-YIG nuclease families. Multiple sequence alignments together with predicted secondary structures of identified families were represented as Hidden Markov Models (HMM) and compared by the HHsearch method to the uncharacterized protein families gathered in the COG, KOG, and PFAM databases. This analysis allowed for extending the GIY-YIG superfamily to include members of COG3680 and a number of proteins not classified in COGs and to predict that these proteins may function as nucleases, potentially involved in DNA recombination and/or repair. Finally, all old and new members of the GIY-YIG superfamily were compared and analyzed to infer the phylogenetic tree.

Conclusion

An evolutionary classification of the GIY-YIG superfamily is presented for the very first time, along with the structural annotation of all (sub)families. It provides a comprehensive picture of sequence-structure-function relationships in this superfamily of nucleases, which will help to design experiments to study the mechanism of action of known members (especially the uncharacterized ones) and will facilitate the prediction of function for the newly discovered ones.

An inactivated nuclease-like domain in RecC with novel function: implications for evolution

BACKGROUND

The PD-(D/E)xK superfamily, containing a wide variety of other exo- and endonucleases, is a notable example of general function conservation in the face of extreme sequence and structural variation. Almost all members employ a small number of shared conserved residues to bind catalytically essential metal ions and thereby effect DNA cleavage. The crystal structure of the RecBCD prokaryotic DNA repair machinery shows that RecB contains such a nuclease domain at its C-terminus. The RecC C-terminal region was reported as having a novel fold.

RESULTS

The RecC C-terminal region can be divided into an alpha/beta domain and a smaller alpha-helical bundle domain. Here we show that the alpha/beta domain is homologous to the RecB nuclease domain but lacks the features necessary for catalysis. Instead, the domain has a novel function within the nuclease superfamily--providing a hoop through which single-stranded DNA passes. Comparison with other structures of nuclease domains bound to DNA reveals strikingly different modes of ligand binding. The alpha-helical bundle domain contributes the pin which splits the DNA duplex.

CONCLUSION

The demonstrated homology of RecB and RecC shows how evolution acted to produce the present RecBCD complex through aggregation of new domains as well as functional divergence and structural redeployment of existing domains. Distantly homologous nuclease(-like) domains bind DNA in highly diverse manners.

Identification of a new family of putative PD-(D/E)XK nucleases with unusual phylogenomic distribution and a new type of the active site

Background

Prediction of structure and function for uncharacterized protein families by identification of evolutionary links to characterized families and known structures is one of the cornerstones of genomics. Theoretical assignment of three-dimensional folds and prediction of protein function even at a very general level can facilitate the experimental determination of the molecular mechanism of action and the role that members of a given protein family fulfill in the cell. Here, we predict the three-dimensional fold and study the phylogenomic distribution of members of a large family of uncharacterized proteins classified in the Clusters of Orthologous Groups database as COG4636.

Results

Using protein fold-recognition we found that members of COG4636 are remotely related to Holliday junction resolvases and other nucleases from the PD-(D/E)XK superfamily. Structure modeling and sequence analyses suggest that most members of COG4636 exhibit a new, unusual variant of the putative active site, in which the catalytic Lys residue migrated in the sequence, but retained similar spatial position with respect to other functionally important residues. Sequence analyses revealed that members of COG4636 and their homologs are found mainly in Cyanobacteria, but also in other bacterial phyla. They undergo horizontal transfer and extensive proliferation in the colonized genomes; for instance in Gloeobacter violaceus PCC 7421 they comprise over 2% of all protein-encoding genes. Thus, members of COG4636 appear to be a new type of selfish genetic elements, which may fulfill an important role in the genome dynamics of Cyanobacteria and other species they invaded. Our analyses provide a platform for experimental determination of the molecular and cellular function of members of this large protein family.

Conclusion

After submission of this manuscript, a crystal structure of one of the COG4636 members was released in the Protein Data Bank (code 1wdj; Idaka, M., Wada, T., Murayama, K., Terada, T., Kuramitsu, S., Shirouzu, M., Yokoyama, S.: Crystal structure of Tt1808 from Thermus thermophilus Hb8, to be published). Our analysis of the Tt1808 structure reveals that we correctly predicted all functionally important features of the COG4636 family, including the membership in the PD-(D/E)xK superfamily of nucleases, the three-dimensional fold, the putative catalytic residues, and the unusual configuration of the active site.

Extensive domain shuffling in transcription regulators of DNA viruses and implications for the origin of fungal APSES transcription factors

BACKGROUND

Viral DNA-binding proteins have served as good models to study the biochemistry of transcription regulation and chromatin dynamics. Computational analysis of viral DNA-binding regulatory proteins and identification of their previously undetected homologs encoded by cellular genomes might lead to a better understanding of their function and evolution in both viral and cellular systems.

RESULTS

The phyletic range and the conserved DNA-binding domains of the viral regulatory proteins of the poxvirus D6R/N1R and baculoviral Bro protein families have not been previously defined. Using computational analysis, we show that the amino-terminal module of the D6R/N1R proteins defines a novel, conserved DNA-binding domain (the KilA-N domain) that is found in a wide range of proteins of large bacterial and eukaryotic DNA viruses. The KilA-N domain is suggested to be homologous to the fungal DNA-binding APSES domain. We provide evidence for the KilA-N and APSES domains sharing a common fold with the nucleic acid-binding modules of the LAGLIDADG nucleases and the amino-terminal domains of the tRNA endonuclease. The amino-terminal module of the Bro proteins is another, distinct DNA-binding domain (the Bro-N domain) that is present in proteins whose domain architectures parallel those of the KilA-N domain-containing proteins. A detailed analysis of the KilA-N and Bro-N domains and the associated domains points to extensive domain shuffling and lineage-specific gene family expansion within DNA virus genomes.

CONCLUSIONS

We define a large class of novel viral DNA-binding proteins and their cellular homologs and identify their domain architectures. On the basis of phyletic pattern analysis we present evidence for a probable viral origin of the fungus-specific cell-cycle regulatory transcription factors containing the APSES DNA-binding domain. We also demonstrate the extensive role of lineage-specific gene expansion and domain shuffling, within a limited set of approximately 24 domains, in the generation of the diversity of virus-specific regulatory proteins.

Three-dimensional modeling of the I-TevI homing endonuclease catalytic domain, a GIY-YIG superfamily member, using NMR restraints and Monte Carlo dynamics

Using a recent version of the SICHO algorithm for in silico protein folding, we made a blind prediction of the tertiary structure of the N-terminal, independently folded, catalytic domain (CD) of the I-TevI homing endonuclease, a representative of the GIY-YIG superfamily of homing endonucleases. The secondary structure of the I-TevI CD has been determined using NMR spectroscopy, but computational sequence analysis failed to detect any protein of known tertiary structure related to the GIY-YIG nucleases (Kowalski et al., Nucleic Acids Res., 1999, 27, 2115-2125). To provide further insight into the structure-function relationships of all GIY-YIG superfamily members, including the recently described subfamily of type II restriction enzymes (Bujnicki et al., Trends Biochem. Sci., 2000, 26, 9-11), we incorporated the experimentally determined and predicted secondary and tertiary restraints in a reduced (side chain only) protein model, which was minimized by Monte Carlo dynamics and simulated annealing. The subsequently elaborated full atomic model of the I-TevI CD allows the available experimental data to be put into a structural context and suggests that the GIY-YIG domain may dimerize in order to bring together the conserved residues of the active site.