Three essential and conserved regions of the group II intron are proximal to the 5'-splice site

DNASCANNER v2: A Web-Based Tool to Analyze the Characteristic Properties of Nucleotide Sequences

Throughout the process of evolution, DNA undergoes the accumulation of distinct mutations, which can often result in highly organized patterns that serve various essential biological functions. These patterns encompass various genomic elements and provide valuable insights into the regulatory and functional aspects of DNA. The physicochemical, mechanical, thermodynamic, and structural properties of DNA sequences play a crucial role in the formation of specific patterns. These properties contribute to the three-dimensional structure of DNA and influence their interactions with proteins, regulatory elements, and other molecules. In this study, we introduce DNASCANNER v2, an advanced version of our previously published algorithm DNASCANNER for analyzing DNA properties. The current tool is built using the FLASK framework in Python language. Featuring a user-friendly interface tailored for nonspecialized researchers, it offers an extensive analysis of 158 DNA properties, including mono/di/trinucleotide frequencies, structural, physicochemical, thermodynamics, and mechanical properties of DNA sequences. The tool provides downloadable results and offers interactive plots for easy interpretation and comparison between different features. We also demonstrate the utility of DNASCANNER v2 in analyzing splice-site junctions, casposon insertion sequences, and transposon insertion sites (TIS) within the bacterial and human genomes, respectively. We also developed a deep learning module for the prediction of potential TIS in a given nucleotide sequence. In the future, we aim to optimize the performance of this prediction model through extensive training on larger data sets.

KnotAli: informed energy minimization through the use of evolutionary information

BACKGROUND

Improving the prediction of structures, especially those containing pseudoknots (structures with crossing base pairs) is an ongoing challenge. Homology-based methods utilize structural similarities within a family to predict the structure. However, their prediction is limited to the consensus structure, and by the quality of the alignment. Minimum free energy (MFE) based methods, on the other hand, do not rely on familial information and can predict structures of novel RNA molecules. Their prediction normally suffers from inaccuracies due to their underlying energy parameters.

RESULTS

We present a new method for prediction of RNA pseudoknotted secondary structures that combines the strengths of MFE prediction and alignment-based methods. KnotAli takes a multiple RNA sequence alignment as input and uses covariation and thermodynamic energy minimization to predict possibly pseudoknotted secondary structures for each individual sequence in the alignment. We compared KnotAli's performance to that of three other alignment-based programs, two that can handle pseudoknotted structures and one control, on a large data set of 3034 RNA sequences with varying lengths and levels of sequence conservation from 10 families with pseudoknotted and pseudoknot-free reference structures. We produced sequence alignments for each family using two well-known sequence aligners (MUSCLE and MAFFT).

CONCLUSIONS

We found KnotAli's performance to be superior in 6 of the 10 families for MUSCLE and 7 of the 10 for MAFFT. While both KnotAli and Cacofold use background noise correction strategies, we found KnotAli's predictions to be less dependent on the alignment quality. KnotAli can be found online at the Zenodo image: https://doi.org/10.5281/zenodo.5794719.

Titration of Adenine in a GA Mismatch with Grand Canonical Simulations

Computational prediction of the pKa of ionizable groups remains a central challenge in biomolecular modeling. Although all-atom fixed-charge force fields could be accurate to describe the interaction network within the biomolecules, proper sampling techniques are required to obtain the thermodynamic information in the (de)protonation event. Sufficient sampling requires an ensemble of structures from simulations with proper treatments of the acid-base equilibria, and the grand canonical simulation technique could be used to model the growth/annihilation of hydrogen atoms by merging Hamiltonians of different protonation states into one simulation ensemble. The electrostatic feature of nucleotide systems is especially difficult to model, and the situation becomes more challenging when the ionizable site is highly perturbed. Although there are many successful predictions obtained from the grand canonical constant pH simulations, few reports focus on highly perturbed nucleotide systems with unconventional base-pair features. In this work, with the discrete constant pH method, we investigate the titration thermodynamics of an adenine in the catalytic triad in a 35-nucleotide single-stranded RNA hairpin, featuring an unconventional GA mismatch and a substantially shifted pKa value. Validation tests are performed with two system setups, both of which provide pKa predictions in good agreement with the experimental value. A single-configuration-based technique is used to calculate the pKa for comparison. The current success indicates the predictive power of the current nucleotide modeling framework.

Visualizing group II intron dynamics between the first and second steps of splicing

Group II introns are ubiquitous self-splicing ribozymes and retrotransposable elements evolutionarily and chemically related to the eukaryotic spliceosome, with potential applications as gene-editing tools. Recent biochemical and structural data have captured the intron in multiple conformations at different stages of catalysis. Here, we employ enzymatic assays, X-ray crystallography, and molecular simulations to resolve the spatiotemporal location and function of conformational changes occurring between the first and the second step of splicing. We show that the first residue of the highly-conserved catalytic triad is protonated upon 5’-splice-site scission, promoting a reversible structural rearrangement of the active site (toggling). Protonation and active site dynamics induced by the first step of splicing facilitate the progression to the second step. Our insights into the mechanism of group II intron splicing parallels functional data on the spliceosome, thus reinforcing the notion that these evolutionarily-related molecular machines share the same enzymatic strategy.

Group II introns are self-splicing ribozymes. Here, the authors employ enzymatic assay, X-ray crystallography and molecular dynamics simulations to show that protonation of the group II intron catalytic triad plays an important role for the transition from the first to the second step of splicing.

The mechanism of splicing as told by group II introns: Ancestors of the spliceosome

Spliceosomal introns and self-splicing group II introns share a common mechanism of intron splicing where two sequential transesterification reactions remove intron lariats and ligate exons. The recent revolution in cryo-electron microscopy (cryo-EM) has allowed visualization of the spliceosome's ribozyme core. Comparison of these cryo-EM structures to recent group II intron crystal structures presents an opportunity to draw parallels between the RNA active site, substrate positioning, and product formation in these two model systems of intron splicing. In addition to shared RNA architectural features, structural similarity between group II intron encoded proteins (IEPs) and the integral spliceosomal protein Prp8 further support a shared catalytic core. These mechanistic and structural similarities support the long-held assertion that group II introns and the eukaryotic spliceosome have a common evolutionary origin. In this review, we discuss how recent structural insights into group II introns and the spliceosome facilitate the chemistry of splicing, highlight similarities between the two systems, and discuss their likely evolutionary connections. This article is part of a Special Issue entitled: RNA structure and splicing regulation edited by Francisco Baralle, Ravindra Singh and Stefan Stamm.

Group II Intron Self-Splicing

Group II introns are large, autocatalytic ribozymes that catalyze RNA splicing and retrotransposition. Splicing by group II introns plays a major role in the metabolism of plants, fungi, and yeast and contributes to genetic variation in many bacteria. Group II introns have played a major role in genome evolution, as they are likely progenitors of spliceosomal introns, retroelements, and other machinery that controls genetic variation and stability. The structure and catalytic mechanism of group II introns have recently been elucidated through a combination of genetics, chemical biology, solution biochemistry, and crystallography. These studies reveal a dynamic machine that cycles progressively through multiple conformations as it stimulates the various stages of splicing. A central active site, containing a reactive metal ion cluster, catalyzes both steps of self-splicing. These studies provide insights into RNA structure, folding, and catalysis, as they raise new questions about the behavior of RNA machines.

Protonation-Dependent Base Flipping at Neutral pH in the Catalytic Triad of a Self-Splicing Bacterial Group II Intron

NMR spectroscopy has revealed pH-dependent structural changes in the highly conserved catalytic domain 5 of a bacterial group II intron. Two adenines with pK(a) values close to neutral pH were identified in the catalytic triad and the bulge. Protonation of the adenine opposite to the catalytic triad is stabilized within a G(syn)-AH(+) (anti) base pair. The pH-dependent anti-to-syn flipping of this G in the catalytic triad modulates the known interaction with the linker region between domains 2 and 3 (J23) and simultaneously the binding of the catalytic Mg(2+) ion to its backbone. Hence, this here identified shifted pK(a) value controls the conformational change between the two steps of splicing.

Evidence for a group II intron-like catalytic triplex in the spliceosome

To catalyze pre-mRNA splicing, U6 small nuclear RNA positions two metals that interact directly with the scissile phosphates. U6 metal ligands correspond stereospecifically to metal ligands within the catalytic domain V of a group II self-splicing intron. Domain V ligands are organized by base-triple interactions, which also juxtapose the 3' splice site with the catalytic metals. However, in the spliceosome, the mechanism for organizing catalytic metals and recruiting the substrate has remained unclear. Here we show by genetics, cross-linking and biochemistry in yeast that analogous triples form in U6 and promote catalytic-metal binding and both chemical steps of splicing. Because the triples include an element that defines the 5' splice site, they also provide a mechanism for juxtaposing the pre-mRNA substrate with the catalytic metals. Our data indicate that U6 adopts a group II intron-like tertiary conformation to catalyze splicing.

Principles of ion recognition in RNA: insights from the group II intron structures

Metal ions promote both RNA folding and catalysis, thus being essential in stabilizing the structure and determining the function of large RNA molecules, including group II introns. The latter are self-splicing metalloribozymes, containing a heteronuclear four-metal-ion center within the active site. In addition to these catalytic ions, group II introns bind many other structural ions, including delocalized ions that bind the RNA diffusively and well-ordered ions that bind the RNA tightly with high occupancy. The latter ions, which can be studied by biophysical methods, have not yet been analyzed systematically. Here, we compare crystal structures of the group IIC intron from Oceanobacillus iheyensis and classify numerous site-bound ions, which are primarily localized in the intron core and near long-range tertiary contacts. Certain ion-binding sites resemble motifs observed in known RNA structures, while others are idiosyncratic to the group II intron. Particularly interesting are (1) ions proximal to the active site, which may participate in splicing together with the catalytic four-metal-ion center, (2) organic ions that bind regions predicted to interact with intron-encoded proteins, and (3) unusual monovalent ions bound to GU wobble pairs, GA mismatches, the S-turn, the tetraloop-receptor, and the T-loop. Our analysis extends the general principles by which ions participate in RNA structural organization and it will aid in the determination and interpretation of future RNA structures.

Visualizing group II intron catalysis through the stages of splicing

Group II introns are self-splicing ribozymes that share a reaction mechanism and a common ancestor with the eukaryotic spliceosome, thereby providing a model system for understanding the chemistry of pre-mRNA splicing. Here we report 14 crystal structures of a group II intron at different stages of catalysis. We provide a detailed mechanism for the first step of splicing, we describe a reversible conformational change between the first and the second steps of splicing, and we present the ligand-free intron structure after splicing in an active state that corresponds to the retrotransposable form of the intron. During each reaction, the reactants are aligned and activated by a heteronuclear four-metal-ion center that contains a metal cluster and obligate monovalent cations, and they adopt a structural arrangement similar to that of protein endonucleases. Based on our data, we propose a model for the splicing cycle and show that it is applicable to the eukaryotic spliceosome.

Emerging views about the molecular structure of the spliceosomal catalytic center

Pre-mRNA splicing occurs in two chemical steps that are catalyzed by a large, dynamic RNA-protein complex called the spliceosome. Initially assembled in a catalytically inactive form, the spliceosome undergoes massive compositional and conformational remodeling, through which disparate RNA elements are re-configured and juxtaposed into a functional catalytic center. The intricate construction of the catalytic center requires the assistance of spliceosomal proteins. Recent structure-function analyses have demonstrated that the yeast-splicing factor Cwc2 is a main player that contacts and shapes the catalytic center of the spliceosome into a functional conformation. With this advance, corroborated by the atomic structure of the evolutionarily related group IIC introns, our understanding of the organization and formation of the spliceosomal catalytic center has progressed to a new level.

Crystal structure of a group II intron in the pre-catalytic state

Group II introns are self-splicing catalytic RNAs that are thought to be ancestral to the spliceosome. Here we report the 3.65-Å crystal structure of the group II intron from Oceanobacillus iheyensis in the pre-catalytic state. The structure reveals the conformation of the 5' splice site in the catalytic core and represents the first structure of an intron prior to the first step of splicing.

DEAD-box protein facilitated RNA folding in vivo

In yeast mitochondria the DEAD-box helicase Mss116p is essential for respiratory growth by acting as group I and group II intron splicing factor. Here we provide the first structure-based insights into how Mss116p assists RNA folding in vivo. Employing an in vivo chemical probing technique, we mapped the structure of the ai5γ group II intron in different genetic backgrounds to characterize its intracellular fold. While the intron adopts the native conformation in the wt yeast strain, we found that the intron is able to form most of its secondary structure, but lacks its tertiary fold in the absence of Mss116p. This suggests that ai5γ is largely unfolded in the mss116-knockout strain and requires the protein at an early step of folding. Notably, in this unfolded state misfolded substructures have not been observed. As most of the protein-induced conformational changes are located within domain D1, Mss116p appears to facilitate the formation of this largest domain, which is the scaffold for docking of other intron domains. These findings suggest that Mss116p assists the ordered assembly of the ai5γ intron in vivo.

RNA trans-splicing: identification of components of a putative chloroplast spliceosome

Group II introns with highly complex RNA structures have been discovered in both prokaryotes and eukaryotic organelles. Usually, excision of non-coding group II intron sequences occurs by cis-splicing, the intramolecular ligation of exons in the same precursor RNA, but some group II introns are excised by intermolecular ligation. This process is called trans-splicing, and genome sequencing predicted that this type of RNA processing occurs in more than 180 organelle genomes from eukaryotes. A well characterised trans-spliced intron RNA is represented by the chloroplast psaA gene of the model alga Chlamydomonas reinhardtii. The psaA gene is split into three exons, which are widely distributed over the plastome and transcribed independently. PsaA exons are flanked by sequences typical for group II introns and joined by trans-splicing via two transesterification reactions. Although it is known that some group II introns are able to splice autocatalytically, trans-splicing of the psaA RNA depends on several nucleus and chloroplast encoded factors. The phylogenetic relationship between group II introns and nuclear spliceosomal RNA led to the hypothesis that these factors are part of large multiprotein and ribonucleoprotein complexes akin to the nuclear spliceosome. Here, we give a concise overview of experimental strategies to identify novel factors involved in trans-splicing of psaA RNA and review recent results that have elucidated the composition and function of a putative chloroplast spliceosome involved in processing of chloroplast precursor RNAs.

Structural insights into RNA splicing

Intron splicing is a fundamental biological process whereby noncoding sequences are removed from precursor RNAs. Recent work has provided new insights into the structural features and reaction mechanisms of two introns that catalyze their own splicing from precursor RNA: the group I and II introns. In addition, there is an increasing amount of structural information on the spliceosome, which is a ribonucleoprotein machine that catalyzes nuclear pre-mRNA splicing in eukaryotes. Here, we compare structures and catalytic mechanisms of self-splicing RNAs and we discuss the possible implications for spliceosomal reaction mechanisms.

A conserved 3' extension in unusual group II introns is important for efficient second-step splicing

The B.c.I4 group II intron from Bacillus cereus ATCC 10987 harbors an unusual 3′ extension. Here, we report the discovery of four additional group II introns with a similar 3′ extension in Bacillus thuringiensis kurstaki 4D1 that splice at analogous positions 53/56 nt downstream of domain VI in vivo. Phylogenetic analyses revealed that the introns are only 47–61% identical to each other. Strikingly, they do not form a single evolutionary lineage even though they belong to the same Bacterial B class. The extension of these introns is predicted to form a conserved two-stem–loop structure. Mutational analysis in vitro showed that the smaller stem S1 is not critical for self-splicing, whereas the larger stem S2 is important for efficient exon ligation and lariat release in presence of the extension. This study clearly demonstrates that previously reported B.c.I4 is not a single example of a specialized intron, but forms a new functional class with an unusual mode that ensures proper positioning of the 3′ splice site.

A glimpse into the active site of a group II intron and maybe the spliceosome, too

The X-ray crystal structure of an excised group II self-splicing intron was recently solved by the Pyle group. Here we review some of the notable features of this structure and what they may tell us about the catalytic active site of the group II ribozyme and potentially the spliceosome. The new structure validates the central role of domain V in both the structure and catalytic function of the ribozyme and resolves several outstanding puzzles raised by previous biochemical, genetic and structural studies. While lacking both exons as well as the cleavage sites and nucleophiles, the structure reveals how a network of tertiary interactions can position two divalent metal ions in a configuration that is ideal for catalysis.

A conserved element that stabilizes the group II intron active site

The internal loop at the base of domain 3 (D3) is one of the most conserved and catalytically important elements of a group II intron. However, the location and molecular nature of its tertiary interaction partners has remained unknown. By employing a combination of site-directed photo-cross-linking and nucleotide analog interference suppression (NAIS), we show that the domain 3 internal loop (D3IL) interacts with the epsilon-epsilon' duplex, which is an active-site element located near the 5'-splice site in D1. Our data also suggest that the D3IL may interact with the bulge of D5, which is a critical active site component. The results of this and other recent studies indicate that the D3IL participates in a complex network of tertiary interactions involving epsilon-epsilon', the bulge of D5 and J23, and that it helps to optimize active site architecture by supporting interactions among these catalytic motifs. Our results are consistent with the role of D3 as a catalytic effector that enhances intron reactivity through active site stabilization.

Crystal structure of a self-spliced group II intron

Group II introns are self-splicing ribozymes that catalyze their own excision from precursor transcripts and insertion into new genetic locations. Here we report the crystal structure of an intact, self-spliced group II intron from Oceanobacillus iheyensis at 3.1 angstrom resolution. An extensive network of tertiary interactions facilitates the ordered packing of intron subdomains around a ribozyme core that includes catalytic domain V. The bulge of domain V adopts an unusual helical structure that is located adjacent to a major groove triple helix (catalytic triplex). The bulge and catalytic triplex jointly coordinate two divalent metal ions in a configuration that is consistent with a two-metal ion mechanism for catalysis. Structural and functional analogies support the hypothesis that group II introns and the spliceosome share a common ancestor.