High-resolution profiling of homing endonuclease binding and catalytic specificity using yeast surface display

Experimental analysis and manipulation of protein–DNA interactions pose unique biophysical challenges arising from the structural and chemical homogeneity of DNA polymers. We report the use of yeast surface display for analytical and selection-based applications for the interaction between a LAGLIDADG homing endonuclease and its DNA target. Quantitative flow cytometry using oligonucleotide substrates facilitated a complete profiling of specificity, both for DNA-binding and catalysis, with single base pair resolution. These analyses revealed a comprehensive segregation of binding specificity and affinity to one half of the pseudo-dimeric interaction, while the entire interface contributed specificity at the level of catalysis. A single round of targeted mutagenesis with tandem affinity and catalytic selection steps provided mechanistic insights to the origins of binding and catalytic specificity. These methods represent a dynamic new approach for interrogating specificity in protein–DNA interactions.     

Function and ribosomal localization of aIF6, a translational regulator shared by archaea and eukarya

The translation factor IF6 is shared by the Archaea and the Eukarya, but is not found in Bacteria. The properties of eukaryal IF6 (eIF6) have been extensively studied, but remain somewhat elusive. eIF6 behaves as a ribosome-anti-association factor and is involved in miRNA-mediated gene silencing; however, it also seems to participate in ribosome synthesis and export. Here we have determined the function and ribosomal localization of the archaeal (Sulfolobus solfataricus) IF6 homologue (aIF6). We find that aIF6 binds specifically to the 50S ribosomal subunits, hindering the formation of 70S ribosomes and strongly inhibiting translation. aIF6 is uniformly expressed along the cell cycle, but it is upregulated following both cold- and heat shock. The aIF6 ribosomal binding site lies in the middle of the 30-S interacting surface of the 50S subunit, including a number of critical RNA and protein determinants involved in subunit association. The data suggest that the IF6 protein evolved in the archaeal–eukaryal lineage to modulate translational efficiency under unfavourable environmental conditions, perhaps acquiring additional functions during eukaryotic evolution.     

Evolutionary Dynamics of the mS952 Intron: A Novel Mitochondrial Group II Intron Encoding a LAGLIDADG Homing Endonuclease Gene

Examination of the mitochondrial small subunit ribosomal RNA (rns) gene of five species of the fungal genus Leptographium revealed that the gene has been invaded at least once at position 952 by a group II intron encoding a LAGLIDADG homing endonuclease gene. Phylogenetic analyses of the intron and homing endonuclease sequences indicated that each element in Leptographium species forms a single clade and is closely related to the group II intron/homing endonuclease gene composite element previously reported at position 952 of the mitochondrial rns gene of Cordyceps species and of Cryphonectria parasitica. The results of an intron survey of the mt rns gene of Leptographium species superimposed onto the phylogenetic analysis of the host organisms suggest that the composite element was transmitted vertically in Leptographium lundbergii. However, its stochastic distribution among strains of L. wingfieldii, L. terebrantis, and L. truncatum suggests that it has been horizontally transmitted by lateral gene transfer among these species, although the random presence of the intron may reflect multiple random loss events. A model is proposed describing the initial invasion of the group II intron in the rns gene of L. lundbergii by a LAGLIDADG homing endonuclease gene and subsequent evolution of this gene to recognize a novel DNA target site, which may now promote the mobility of the intron and homing endonuclease gene as a composite element.     

Simplification of Ribosomes in Bacteria with Tiny Genomes

The ribosome is an essential cellular machine performing protein biosynthesis. Its structure and composition are highly conserved in all species. However, some bacteria have been reported to have an incomplete set of ribosomal proteins. We have analyzed ribosomal protein composition in 214 small bacterial genomes (<1 Mb) and found that although the ribosome composition is fairly stable, some ribosomal proteins may be absent, especially in bacteria with dramatically reduced genomes. The protein composition of the large subunit is less conserved than that of the small subunit. We have identified the set of frequently lost ribosomal proteins and demonstrated that they tend to be positioned on the ribosome surface and have fewer contacts to other ribosome components. Moreover, some proteins are lost in an evolutionary correlated manner. The reduction of ribosomal RNA is also common, with deletions mostly occurring in free loops. Finally, the loss of the anti-Shine–Dalgarno sequence is associated with the loss of a higher number of ribosomal proteins.     

Rapid evolution of the DNA-binding site in LAGLIDADG homing endonucleases

Sequence analysis of chloroplast and mitochondrial large subunit rRNA genes from over 75 green algae disclosed 28 new group I intron-encoded proteins carrying a single LAGLIDADG motif. These putative homing endonucleases form four subfamilies of homologous enzymes, with the members of each subfamily being encoded by introns sharing the same insertion site. We showed that four divergent endonucleases from the I-CreI subfamily cleave the same DNA substrates. Mapping of the 66 amino acids that are conserved among the members of this subfamily on the 3-dimensional structure of I-CreI bound to its recognition sequence revealed that these residues participate in protein folding, homodimerization, DNA recognition and catalysis. Surprisingly, only seven of the 21 I-CreI amino acids interacting with DNA are conserved, suggesting that I-CreI and its homologs use different subsets of residues to recognize the same DNA sequence. Our sequence comparison of all 45 single-LAGLIDADG proteins identified so far suggests that these proteins share related structures and that there is a weak pressure in each subfamily to maintain identical protein–DNA contacts. The high sequence variability we observed in the DNA-binding site of homologous LAGLIDADG endonucleases provides insight into how these proteins evolve new DNA specificity.     

Mycobacterium tuberculosis RecA intein,a LAGLIDADG homing endonuclease,displays Mn(2+) and DNA-dependent ATPase activity

Mycobacterium tuberculosis RecA intein (PI-MtuI), a LAGLIDADG homing endonuclease, displays dual target specificity in response to alternative cofactors. While both ATP and Mn²⁺ were required for optimal cleavage of an inteinless recA allele (hereafter referred to as cognate DNA), Mg²⁺ alone was sufficient for cleavage of ectopic DNA sites. In this study, we have explored the ability of PI-MtuI to catalyze ATP hydrolysis in the presence of alternative metal ion cofactors and DNA substrates. Our results indicate that PI-MtuI displays maximum ATPase activity in the presence of cognate but not ectopic DNA. Kinetic analysis revealed that Mn²⁺ was able to stimulate PI-MtuI catalyzed ATP hydrolysis, whereas Mg²⁺ failed to do so. Using UV crosslinking, limited proteolysis and amino acid sequence analysis, we show that ³²P-labeled ATP was bound to a 14 kDa peptide containing the putative Walker A motif. Furthermore, the limited proteolysis approach disclosed that cognate DNA was able to induce structural changes in PI-MtuI. Mutation of the presumptive metal ion-binding ligands (Asp122 and Asp222) in the LAGLIDADG motifs of PI-MtuI impaired its affinity for ATP, thus resulting in a reduction in or loss of its endonuclease activity. Together, these results suggest that PI-MtuI is a (cognate) DNA- and Mn²⁺-dependent ATPase, unique from the LAGLIDADG family of homing endonucleases, and implies a possible role for ATP hydrolysis in the recognition and/or cleavage of homing site DNA sequence.     

Coevolution of a homing endonuclease and its host target sequence

We have determined the specificity profile of the homing endonuclease I-AniI and compared it to the conservation of its host gene. Homing endonucleases are encoded within intervening sequences such as group I introns. They initiate the transfer of such elements by cleaving cognate alleles lacking the intron, leading to their transfer via homologous recombination. Each structural homing endonuclease family has arrived at an appropriate balance of specificity and fidelity that avoids toxicity while maximizing target recognition and invasiveness. I-AniI recognizes a strongly conserved target sequence in a host gene encoding apocytochrome B and has fine-tuned its specificity to correlate with wobble versus nonwobble positions across that sequence and to the amount of degeneracy inherent in individual codons. The physiological target site in the host gene is not the optimal substrate for recognition and cleavage: at least one target variant identified during a screen is bound more tightly and cleaved more rapidly. This is a result of the periodic cycle of intron homing, which at any time can present nonoptimal combinations of endonuclease specificity and insertion site sequences in a biological host.     

RNomics in Archaea reveals a further link between splicing of archaeal introns and rRNA processing

The bulge–helix–bulge (BHB) motif recognised by the archaeal splicing endonuclease is also found in the long processing stems of archaeal rRNA precursors in which it is cleaved to generate pre-16S and pre-23S rRNAs. We show that in two species, Archaeoglobus fulgidus and Sulfolobus solfataricus, representatives from the two major archaeal kingdoms Euryarchaeota and Crenarchaeota, respectively, the pre-rRNA spacers cleaved at the BHB motifs surrounding pre-16S and pre-23S rRNAs subsequently become ligated. In addition, we present evidence that this is accompanied by circularisation of ribosomal pre-16S and pre-23S rRNAs in both species. These data reveal a further link between intron splicing and pre-rRNA processing in Archaea, which might reflect a common evolutionary origin of the two processes. One spliced RNA species designated 16S-D RNA, resulting from religation at the BHB motif of 16S pre-rRNA, is a highly abundant and stable RNA which folds into a three-stem structure interrupted by two single-stranded regions as assessed by chemical probing. It spans a region of the pre-rRNA 5′ external transcribed spacer exhibiting a highly conserved folding pattern in Archaea. Surprisingly, 16S-D RNA contains structural motifs found in archaeal C/D box small RNAs and binds to the L7Ae protein, a core component of archaeal C/D box RNPs. This supports the notion that it might have an important but still unknown role in pre-rRNA biogenesis or might even target RNA molecules other than rRNA.     

Archaeal ribosomal protein L7 is a functional homolog of the eukaryotic 15.5kD/Snu13p snoRNP core protein

Recent investigations have identified homologs of eukaryotic box C/D small nucleolar RNAs (snoRNAs) in Archaea termed sRNAs. Archaeal homologs of the box C/D snoRNP core proteins fibrillarin and Nop56/58 have also been identified but a homolog for the eukaryotic 15.5kD snoRNP protein has not been described. Our sequence analysis of archaeal genomes reveals that the highly conserved ribosomal protein L7 exhibits extensive homology with the eukaryotic 15.5kD protein. Protein binding studies demonstrate that recombinant Methanoccocus jannaschii L7 protein binds the box C/D snoRNA core motif with the same specificity and affinity as the eukaryotic 15.5kD protein. Identical to the eukaryotic 15.5kD core protein, archaeal L7 requires a correctly folded box C/D core motif and intact boxes C and D. Mutational analysis demonstrates that critical features of the box C/D core motif essential for 15.5kD binding are also required for L7 interaction. These include stem I which juxtaposes boxes C and D, as well as the sheared G:A pairs and protruded pyrimidine nucleotide of the asymmetric bulge region. The demonstrated presence of L7Ae in the Haloarcula marismortui 50S ribosomal subunit, taken with our demonstration of the ability of L7 to bind to the box C/D snoRNA core motif, indicates that this protein serves a dual role in Archaea. L7 functioning as both an sRNP core protein and a ribosomal protein could potentially regulate and coordinate sRNP assembly with ribosome biogenesis.     

Cryo-EM structure of the archaeal 50S ribosomal subunit in complex with initiation factor 6 and implications for ribosome evolution

Translation of mRNA into proteins by the ribosome is universally conserved in all cellular life. The composition and complexity of the translation machinery differ markedly between the three domains of life. Organisms from the domain Archaea show an intermediate level of complexity, sharing several additional components of the translation machinery with eukaryotes that are absent in bacteria. One of these translation factors is initiation factor 6 (IF6), which associates with the large ribosomal subunit. We have reconstructed the 50S ribosomal subunit from the archaeon Methanothermobacter thermautotrophicus in complex with archaeal IF6 at 6.6?? resolution using cryo-electron microscopy (EM). The structure provides detailed architectural insights into the 50S ribosomal subunit from a methanogenic archaeon through identification of the rRNA expansion segments and ribosomal proteins that are shared between this archaeal ribosome and eukaryotic ribosomes but are mostly absent in bacteria and in some archaeal lineages. Furthermore, the structure reveals that, in spite of highly divergent evolutionary trajectories of the ribosomal particle and the acquisition of novel functions of IF6 in eukaryotes, the molecular binding of IF6 on the ribosome is conserved between eukaryotes and archaea. The structure also provides a snapshot of the reductive evolution of the archaeal ribosome and offers new insights into the evolution of the translation system in archaea.     

Structure-function relationships in the ribosomal protein L12 in the archaeon Sulfolobus acidocaldarius.

A series of mutant L12 ribosomal proteins was prepared by site-directed mutations in the L12 protein gene of the archaeon Sulfolobus acidocaldarius. The mutant protein genes were overexpressed in Escherichia coli, and the products purified and incorporated into ribosomal cores which had been ethanol extracted to remove wild-type L12 protein. Measurements were made to determine if the mutation affected the binding of the L12 protein to the ribosome core or affected the translational activity of the resulting ribosome. Changing tyrosine [3] or tyrosine [5], conserved in all archaea and present in all eukarya in positions [3] and [7], to phenylalanine had no effect on binding or translational activity while changes to glycine significantly reduced binding and translational activity. Changing the single arginine [37] residue, conserved in almost all archaeal and eukaryal L12 proteins, to lysine, glutamic acid, glutamine, or glycine had no effect on binding to the core and had little or no significant effect on translational activity. The same was true when lysine [39], conserved in all archaeal L12 proteins, was changed to arginine, glutamic acid, glutamine, or glycine. Changing phenylalanine [104], the penultimate amino acid at the C-terminal end, which is conserved in all archaeal and eukaryal L12 proteins, to tyrosine or glycine had no effect on binding but lowered the translational activity by 60 and 75%, respectively, suggesting that this amino acid plays an important role in translation. Deletion of the highly charged region in the C-terminal domain, which is present in all archaeal and eukaryal L12 proteins, decreased transitional activity by 50%, suggesting this region is also involved in factor interactions.     

Conserved functional domains and a novel tertiary interaction near the pseudoknot drive translational activity of hepatitis C virus and hepatitis C virus-like internal ribosome entry sites

The translational activity of the hepatitis C virus (HCV) internal ribosome entry site (IRES) and other HCV-like IRES RNAs depends on structured RNA elements in domains II and III, which serve to recruit the ribosomal 40S subunit, eukaryotic initiation factor (eIF) 3 and the ternary eIF2/Met-tRNA_i^Met/GTP complex and subsequently domain II assists subunit joining. Porcine teschovirus-1 talfan (PTV-1) is a member of the Picornaviridae family, with a predicted HCV-like secondary structure, but only stem-loops IIId and IIIe in the 40S-binding domain display significant sequence conservation with the HCV IRES. Here, we use chemical probing to show that interaction sites with the 40S subunit and eIF3 are conserved between HCV and HCV-like IRESs. In addition, we reveal the functional role of a strictly conserved co-variation between a purine–purine mismatch near the pseudoknot (A–A/G) and the loop sequence of domain IIIe (GAU/CA). These nucleotides are involved in a tertiary interaction, which serves to stabilize the pseudoknot structure and correlates with translational efficiency in both the PTV-1 and HCV IRES. Our data demonstrate conservation of functional domains in HCV and HCV-like IRESs including a more complex structure surrounding the pseudoknot than previously assumed.     

An evolutionary conserved pattern of 18S rRNA sequence complementarity to mRNA 5′ UTRs and its implications for eukaryotic gene translation regulation

There are several key mechanisms regulating eukaryotic gene expression at the level of protein synthesis. Interestingly, the least explored mechanisms of translational control are those that involve the translating ribosome per se, mediated for example via predicted interactions between the ribosomal RNAs (rRNAs) and mRNAs. Here, we took advantage of robustly growing large-scale data sets of mRNA sequences for numerous organisms, solved ribosomal structures and computational power to computationally explore the mRNA–rRNA complementarity that is statistically significant across the species. Our predictions reveal highly specific sequence complementarity of 18S rRNA sequences with mRNA 5′ untranslated regions (UTRs) forming a well-defined 3D pattern on the rRNA sequence of the 40S subunit. Broader evolutionary conservation of this pattern may imply that 5′ UTRs of eukaryotic mRNAs, which have already emerged from the mRNA-binding channel, may contact several complementary spots on 18S rRNA situated near the exit of the mRNA binding channel and on the middle-to-lower body of the solvent-exposed 40S ribosome including its left foot. We discuss physiological significance of this structurally conserved pattern and, in the context of previously published experimental results, propose that it modulates scanning of the 40S subunit through 5′ UTRs of mRNAs.     

Statistical modeling and analysis of the LAGLIDADG family of site-specific endonucleases and identification of an intein that encodes a site-specific endonuclease of the HNH family.

The LAGLIDADG and HNH families of site-specific DNA endonucleases encoded by viruses, bacteriophages as well as archaeal, eucaryotic nuclear and organellar genomes are characterized by the sequence motifs 'LAGLIDADG' and 'HNH', respectively. These endonucleases have been shown to occur in different environments: LAGLIDADG endonucleases are found in inteins, archaeal and group I introns and as free standing open reading frames (ORFs); HNH endonucleases occur in group I and group II introns and as ORFs. Here, statistical models (hidden Markov models, HMMs) that encompass both the conserved motifs and more variable regions of these families have been created and employed to characterize known and potential new family members. A number of new, putative LAGLIDADG and HNH endonucleases have been identified including an intein-encoded HNH sequence. Analysis of an HMM-generated multiple alignment of 130 LAGLIDADG family members and the three-dimensional structure of the I- Cre I endonuclease has enabled definition of the core elements of the repeated domain (approximately 90 residues) that is present in this family of proteins. A conserved negatively charged residue is proposed to be involved in catalysis. Phylogenetic analysis of the two families indicates a lack of exchange of endonucleases between different mobile elements (environments) and between hosts from different phylogenetic kingdoms. However, there does appear to have been considerable exchange of endonuclease domains amongst elements of the same type. Such events are suggested to be important for the formation of elements of new specficity.     

Homology modeling and mutational analysis of Ho endonuclease of yeast

Ho endonuclease is a LAGLIDADG homing endonuclease that initiates mating-type interconversion in yeast. Ho is encoded by a free-standing gene but shows 50% primary sequence similarity to the intein (protein-intron encoded) PI-SceI. Ho is unique among LAGLIDADG endonucleases in having a 120-residue C-terminal putative zinc finger domain. The crystal structure of PI-SceI revealed a bipartite enzyme with a protein-splicing domain (Hint) and intervening endonuclease domain. We made a homology model for Ho on the basis of the PI-SceI structure and performed mutational analysis of putative critical residues, using a mating-type switch as a bioassay for activity and GFP-fusion proteins to detect nuclear localization. We found that residues of the N-terminal sequence of the Hint domain are important for Ho activity, in particular the DNA recognition region. C-terminal residues of the Hint domain are dispensable for Ho activity; however, the C-terminal putative zinc finger domain is essential. Mutational analysis indicated that residues in Ho that are conserved relative to catalytic, active-site residues in PI-SceI and other related homing endonucleases are essential for Ho activity. Our results indicate that in addition to the conserved catalytic residues, Hint domain residues and the zinc finger domain have evolved a critical role in Ho activity.     

Protein footprinting approach to mapping DNA binding sites of two archaeal homing enzymes: evidence for a two-domain protein structure.

The archaeal intron-encoded homing enzymes I-PorI and I-DmoI belong to a family of endonucleases that contain two copies of a characteristic LAGLIDADG motif. These endonucleases cleave their intron- or intein- alleles site-specifically, and thereby facilitate homing of the introns or inteins which encode them. The protein structure and the mechanism of DNA recognition of these homing enzymes is largely unknown. Therefore, we examined these properties of I-PorI and I-DmoI by protein footprinting. Both proteins were susceptible to proteolytic cleavage within regions that are equidistant from each of the two LAGLIDADG motifs. When complexed with their DNA substrates, a characteristic subset of the exposed sites, located in regions immediately after and 40-60 amino acids after each of the LAGLIDADG motifs, were protected. Our data suggest that the enzymes are structured into two, tandemly repeated, domains, each containing both the LAGLIDADG motif and two putative DNA binding regions. The latter contains a potentially novel DNA binding motif conserved in archaeal homing enzymes. The results are consistent with a model where the LAGLIDADG endonucleases bind to their non-palindromic substrates as monomeric enzymes, with each of the two domains recognizing one half of the DNA substrate.     

From monomeric to homodimeric endonucleases and back: engineering novel specificity of LAGLIDADG enzymes

Monomeric homing endonucleases of the LAGLIDADG family recognize DNA in a bipartite manner, reflecting the underlying structural assembly of two protein domains (A and B) related by pseudo 2-fold symmetry. This architecture allows for changes in DNA specificity via the distinct combination of these half-site domains. The key to engineering such hybrid proteins lies in the LAGLIDADG two-helix bundle that forms both the domain interface and the endonuclease active site. In this study, we utilize domain A of the monomeric I-DmoI to demonstrate the feasibility of generating functional homodimeric endonucleases that recognize palindromic DNA sequences derived from the original, non-palindromic target. Wild-type I-DmoI domain A is capable of forming a homodimer (H-DmoA) that binds tightly to, but does not cleave efficiently, its anticipated DNA target. Partial restoration of DNA cleavage ability was obtained by re-engineering the LAGLIDADG dimerization interface (H-DmoC). Upon fusing two copies of H-DmoC via a short peptide linker, a novel, site-specific DNA endonuclease was created (H-DmoC2). Like I-DmoI, H-DmoC2 is thermostable and cleaves the new target DNA to generate the predicted 4 nt 3'-OH overhangs but, unlike I-DmoI, H-DmoC2 retains stringent cleavage specificity when substituting Mn2+ for Mg2+ as co-factor. This novel endonuclease allows speculation regarding specificity of monomeric LAGLIDADG proteins, while it supports the evolutionary genesis of these proteins by a gene duplication event.     

The archaeal eIF2 homologue: functional properties of an ancient translation initiation factor

The eukaryotic translation initiation factor 2 (eIF2) is pivotal for delivery of the initiator tRNA (tRNAi) to the ribosome. Here, we report the functional characterization of the archaeal homologue, a/eIF2. We have cloned the genes encoding the three subunits of a/eIF2 from the thermophilic archaeon Sulfolobus solfataricus, and have assayed the activities of the purified recombinant proteins in vitro. We demonstrate that the trimeric factor reconstituted from the recombinant polypeptides has properties similar to those of its eukaryal homologue: it interacts with GTP and Met-tRNAi, and stimulates binding of the latter to the small ribosomal subunit. However, the archaeal protein differs in some functional aspects from its eukaryal counterpart. In contrast to eIF2, a/eIF2 has similar affinities for GDP and GTP, and the β-subunit does not contribute to tRNAi binding. The detailed analysis of the complete trimer and of its isolated subunits is discussed in light of the evolutionary history of the eIF2-like proteins.     

Purification and characterization of the DNA cleavage and recognition site of I-ScaI mitochondrial group I intron encoded endonuclease produced in Escherichia coli

The second intron in the mitochondrial cytb gene of Saccharomyces capensis, belonging to group I, encodes a 280 amino acid protein containing two LAGLIDADG motifs. Genetic and molecular studies have previously shown that this protein has a dual function in the wild-type strain. It acts as a specific homing endonuclease I-ScaI promoting intron mobility and as a maturase promoting intron splicing. Here we describe the synthesis of a universal code equivalent to the mitochondrial sequence coding for this protein and the in vitro characterization of I-ScaI endonuclease activity, using a truncated mutant form of the protein p28bi2 produced in Escherichia coli. We have also determined the cleavage pattern as well as the recognition site of p28bi2. It was found that p28bi2 generates a double-strand cleavage downstream from the intron insertion site with 4 nt long 3′-overhangs. Mutational analysis of the DNA target site shows that p28bi2 recognizes a 16–19 bp sequence from positions –11 to +8 with respect to the intron insertion site.     

I-OmiI and I-OmiII: Two intron-encoded homing endonucleases within the Ophiostoma minus rns gene

The mitochondrial small subunit ribosomal RNA (rns) gene of the ascomycetous fungus Ophiostoma minus [strain WIN(M)371] was found to contain a group IC2 and a group IIB1 intron at positions mS569 and mS952 respectively. Both introns have open reading frames (ORFs) embedded that encode double motif LAGLIDADG homing endonucleases (I-OmiI and I-OmiII respectively). Codon-optimized versions of I-OmiI and I-OmiII were synthesized for overexpression in Escherichia coli. The in vitro characterization of I-OmiII showed that it is a functional homing endonuclease that cleaves the rns target site two nucleotides upstream (sense strand) of the intron insertion site generating 4 nucleotide 3′ overhangs. The endonuclease activity of I-OmiII was tested using linear and circular substrates and cleavage activity was evaluated at various temperatures. The I-OmiI protein was expressed in E. coli, but purification was difficult, thus the endonuclease activity of this protein was tested via in vivo assays. Overall this study showed that there are many native forms of functional homing endonucleases yet to be discovered among fungal mtDNA genomes.