The guanosine binding site of the Tetrahymena ribozyme

Programmable trans-splicing riboregulators for complex cellular logic computation

Synthetic genetic circuits program the cellular input-output relationships to execute customized functions. However, efforts to scale up these circuits have been hampered by the limited number of reliable regulatory mechanisms with high programmability, performance, predictability and orthogonality. Here we report a class of split-intron-enabled trans-splicing riboregulators (SENTRs) based on de novo designed external guide sequences. SENTR libraries provide low leakage expression, wide dynamic range, high predictability with machine learning and low crosstalk at multiple component levels. SENTRs can sense RNA targets, process signals by logic computation and transduce them into various outputs, either mRNAs or noncoding RNAs. We subsequently demonstrate that digital logic operation with up to six inputs can be implemented using multiple orthogonal SENTRs to regulate a single gene simultaneously and coupling SENTRs with split intein-mediated protein trans-splicing. SENTR represents a powerful and versatile regulatory tool at the post-transcriptional level in Escherichia coli, suggesting broad biotechnological applications.

Efficient circularization of protein-encoding RNAs via a novel cis-splicing system

Circular RNAs (circRNAs) have emerged as a promising alternative to linear mRNA, owing to their unique properties and potential therapeutic applications, driving the development of novel approaches for their production. This study introduces a cis-splicing system that efficiently produces circRNAs by incorporating a ribozyme core at one end of the precursor, thereby eliminating the need for additional spacer elements between the ribozyme and the gene of interest (GOI). In this cis-splicing system, sequences resembling homologous arms at both ends of the precursor are crucial for forming the P9.0 duplex, which in turn facilitates effective self-splicing and circularization. We demonstrate that the precise recognition of the second transesterification site depends more on the structural characteristics of P9.0 adjacent to the ωG position than on the nucleotide composition of the P9.0-ωG itself. Further optimization of structural elements, like P10 and P1-ex, significantly improves circularization efficiency. The circRNAs generated through the cis-splicing system exhibit prolonged protein expression and minimal activation of the innate immune response. This study provides a comprehensive exploration of circRNA generation via a novel strategy and offers valuable insights into the structural engineering of RNA, paving the way for future advancements in circRNA-based applications.

Mechanistic Insights on Metformin and Arginine Implementation as Repurposed Drugs in Glioblastoma Treatment

As the most common and aggressive primary malignant brain tumor, glioblastoma is still lacking a satisfactory curative approach. The standard management consisting of gross total resection followed by radiotherapy and chemotherapy with temozolomide only prolongs patients' life moderately. In recent years, many therapeutics have failed to give a breakthrough in GBM treatment. In the search for new treatment solutions, we became interested in the repurposing of existing medicines, which have established safety profiles. We focused on the possible implementation of well-known drugs, metformin, and arginine. Metformin is widely used in diabetes treatment, but arginine is mainly a cardiovascular protective drug. We evaluated the effects of metformin and arginine on total DNA methylation, as well as the oxidative stress evoked by treatment with those agents. In glioblastoma cell lines, a decrease in 5-methylcytosine contents was observed with increasing drug concentration. When combined with temozolomide, both guanidines parallelly increased DNA methylation and decreased 8-oxo-deoxyguanosine contents. These effects can be explained by specific interactions of the guanidine group with m5CpG dinucleotide. We showed that metformin and arginine act on the epigenetic level, influencing the foreground and potent DNA regulatory mechanisms. Therefore, they can be used separately or in combination with temozolomide, in various stages of disease, depending on desired treatment effects.

Highly Reactive Group I Introns Ubiquitous in Pathogenic Fungi

Systemic fungal infections are a growing public health threat, and yet viable antifungal drug targets are limited as fungi share a similar proteome with humans. However, features of RNA metabolism and the noncoding transcriptomes in fungi are distinctive. For example, fungi harbor highly structured RNA elements that humans lack, such as self-splicing introns within key housekeeping genes in the mitochondria. However, the location and function of these mitochondrial riboregulatory elements has largely eluded characterization. Here we used an RNA-structure-based bioinformatics pipeline to identify the group I introns interrupting key mitochondrial genes in medically relevant fungi, revealing their fixation within a handful of genetic hotspots and their ubiquitous presence across divergent phylogenies of fungi, including all highest priority pathogens such as Candida albicans, Candida auris, Aspergillus fumigatus and Cryptococcus neoformans. We then biochemically characterized two representative introns from C. albicans and C. auris, demonstrating their exceptionally efficient splicing catalysis relative to previously-characterized group I introns. Indeed, the C. albicans mitochondrial intron displays extremely rapid catalytic turnover, even at ambient temperatures and physiological magnesium ion concentrations. Our results unmask a significant new set of players in the RNA metabolism of pathogenic fungi, suggesting a promising new type of antifungal drug target.

Snapshots of the second-step self-splicing of Tetrahymena ribozyme revealed by cryo-EM

Group I introns are catalytic RNAs that coordinate two consecutive transesterification reactions for self-splicing. To understand how the group I intron promotes catalysis and coordinates self-splicing reactions, we determine the structures of L-16 Tetrahymena ribozyme in complex with a 5'-splice site analog product and a 3'-splice site analog substrate using cryo-EM. We solve six conformations from a single specimen, corresponding to different splicing intermediates after the first ester-transfer reaction. The structures reveal dynamics during self-splicing, including large conformational changes of the internal guide sequence and the J5/4 junction as well as subtle rearrangements of active-site metals and the hydrogen bond formed between the 2'-OH group of A261 and the N2 group of guanosine substrate. These results help complete a detailed structural and mechanistic view of this paradigmatic group I intron undergoing the second step of self-splicing.

A DFT study of vibrational spectra of 5-chlorouracil with molecular structure, HOMO-LUMO, MEPs/ESPs and thermodynamic properties

The density functional theory calculation has been carried out for the analysis of 5-chlorouracil using DFT/Gaussian 09 with GAR2PED. Recorded experimental spectra for Raman and IR of 5-chlorouracil have been analyzed all fundamental vibrational modes using the outcome results of DFT at 6-311++G** of Gaussian 09 calculations and the GaussView 5.09. To help the analysis of vibrational modes, GAR2PED program has been used in the calculation of PEDs. The charge transfer properties of 5-chlorouracil have been analyzed using HOMO and LUMO level energy analysis. HOMO and LUMO energy gap study supports the charge transfer possibility in molecule. These have been made to study for reactivity and stability of heterocyclic molecules for the analysis of antiviral drugs against the new corona virus: COVID-19. Here, the smaller energy gap of 5-chlorouracil is more responsible for charge transfer interaction in the heterocyclic drug molecules and a reason of more bioactivity. The electron density mapping within molecular electrostatic potential plot and electrostatic potential plotting within iso-surface plot have been evaluated the charge distribution concept in the molecule as the nucleophilic reactions and electrophilic sites. These computations have been used to produce the molecular charges, structure and thermodynamic functions of biomolecule. This study has been made to all internal modes of chloro group substituent at pyrimidine ring of C₅ atom. The splitting of frequencies has arisen in the two species for the normal distribution modes.

Graphical abstract

Refining Mitochondrial Intron Classification With ERPIN: Identification Based on Conservation of Sequence Plus Secondary Structure Motifs

Mitochondrial genomes—in particular those of fungi—often encode genes with a large number of Group I and Group II introns that are conserved at both the sequence and the RNA structure level. They provide a rich resource for the investigation of intron and gene structure, self- and protein-guided splicing mechanisms, and intron evolution. Yet, the degree of sequence conservation of introns is limited, and the primary sequence differs considerably among the distinct intron sub-groups. It makes intron identification, classification, structural modeling, and the inference of gene models a most challenging and error-prone task—frequently passed on to an “expert” for manual intervention. To reduce the need for manual curation of intron structures and mitochondrial gene models, computational methods using ERPIN sequence profiles were initially developed in 2007. Here we present a refinement of search models and alignments using the now abundant publicly available fungal mtDNA sequences. In addition, we have tested in how far members of the originally proposed sub-groups are clearly distinguished and validated by our computational approach. We confirm clearly distinct mitochondrial Group I sub-groups IA1, IA3, IB3, IC1, IC2, and ID. Yet, IB1, IB2, and IB4 ERPIN models are overlapping substantially in predictions, and are therefore combined and reported as IB. We have further explored the conversion of our ERPIN profiles into covariance models (CM). Current limitations and prospects of the CM approach will be discussed.

The Biochemical Landscape of Riboswitch Ligands

More than 55 distinct classes of riboswitches that respond to small metabolites or elemental ions have been experimentally validated to date. The ligands sensed by these riboswitches are biased in favor of fundamental compounds or ions that are likely to have been relevant to ancient forms of life, including those that might have populated the "RNA World", which is a proposed biochemical era that predates the evolutionary emergence of DNA and proteins. In the following text, I discuss the various types of ligands sensed by some of the most common riboswitches present in modern bacterial cells and consider implications for ancient biological processes centered on the proven capabilities of these RNA-based sensors. Although most major biochemical aspects of metabolism are represented by known riboswitch classes, there are striking sensory gaps in some key areas. These gaps could reveal weaknesses in the performance capabilities of RNA that might have hampered RNA World evolution, or these could highlight opportunities to discover additional riboswitch classes that sense essential metabolites.

Optimized Loopable Translation as a Platform for the Synthesis of Repetitive Proteins

The expression of long proteins with repetitive amino acid sequences often presents a challenge in recombinant systems. To overcome this obstacle, we report a genetic construct that circularizes mRNA in vivo by rearranging the topology of a group I self-splicing intron from T4 bacteriophage, thereby enabling "loopable" translation. Using a fluorescence-based assay to probe the translational efficiency of circularized mRNAs, we identify several conditions that optimize protein expression from this system. Our data suggested that translation of circularized mRNAs could be limited primarily by the rate of ribosomal initiation; therefore, using a modified error-prone PCR method, we generated a library that concentrated mutations into the initiation region of circularized mRNA and discovered mutants that generated markedly higher expression levels. Combining our rational improvements with those discovered through directed evolution, we report a loopable translator that achieves protein expression levels within 1.5-fold of the levels of standard vectorial translation. In summary, our work demonstrates loopable translation as a promising platform for the creation of large peptide chains, with potential utility in the development of novel protein materials.

Fitness landscape of a dynamic RNA structure

RNA structures are dynamic. As a consequence, mutational effects can be hard to rationalize with reference to a single static native structure. We reasoned that deep mutational scanning experiments, which couple molecular function to fitness, should capture mutational effects across multiple conformational states simultaneously. Here, we provide a proof-of-principle that this is indeed the case, using the self-splicing group I intron from Tetrahymena thermophila as a model system. We comprehensively mutagenized two 4-bp segments of the intron. These segments first come together to form the P1 extension (P1ex) helix at the 5' splice site. Following cleavage at the 5' splice site, the two halves of the helix dissociate to allow formation of an alternative helix (P10) at the 3' splice site. Using an in vivo reporter system that couples splicing activity to fitness in E. coli, we demonstrate that fitness is driven jointly by constraints on P1ex and P10 formation. We further show that patterns of epistasis can be used to infer the presence of intramolecular pleiotropy. Using a machine learning approach that allows quantification of mutational effects in a genotype-specific manner, we demonstrate that the fitness landscape can be deconvoluted to implicate P1ex or P10 as the effective genetic background in which molecular fitness is compromised or enhanced. Our results highlight deep mutational scanning as a tool to study alternative conformational states, with the capacity to provide critical insights into the structure, evolution and evolvability of RNAs as dynamic ensembles. Our findings also suggest that, in the future, deep mutational scanning approaches might help reverse-engineer multiple alternative or successive conformations from a single fitness landscape.

Thermodynamic functions with other properties and vibrational spectra of pyrimidine ring of uracil for RNA and bio-molecule 5-aminouracil

Vibrational spectra (Infra Red and Raman) of uracil and 5-aminouracil have been recorded as well as reproduced in solid phase for the pyrimidinal region (200-2000 cm⁻¹) and the internal modes were analyzed especially for pyrimidine ring of bio-molecules theoretically by using the density functional theory (DFT) calculations on the different basis sets for Restricted Hartree–Fock (RHF) by employing the Gaussian-03 program. These quantum chemical calculations have been employed to yield the Mulliken atomic charges, atomic polarizability tensor (APT) charges and thermodynamic functional properties of these bio-molecules. These DFT calculations have been computed for the some important internal ring modes especially to pyrimidine ring of uracil with a substituent amino-group (NH₂) at site of C₅ atom on pyrimidine ring.

Enhancement of RNA/Ligand Association Kinetics via an Electrostatic Anchor

The diverse biological processes mediated by RNA rest upon its recognition of various ligands, including small molecules and nucleic acids. Nevertheless, a recent literature survey suggests that RNA molecular recognition of these ligands is slow, with association rate constants orders of magnitude below the diffusional limit. Thus, we were prompted to consider strategies for increasing RNA association kinetics. Proteins can accelerate ligand association via electrostatic forces, and here, using the Tetrahymena group I ribozyme, we provide evidence that electrostatic forces can accelerate RNA/ligand association. This RNA enzyme (E) catalyzes cleavage of an oligonucleotide substrate (S) by an exogenous guanosine (G) cofactor. The G 2′- and 3′-OH groups interact with an active site metal ion, termed M_C, within E·S·G, and we perturbed each of these contacts via −NH₃⁺ substitution. New and prior data indicate that G(2′NH₃⁺) and G(3′NH₃⁺) bind as strongly as G, suggesting that the −NH₃⁺ substituents of these analogues avoid repulsive interactions with M_C and make alternative interactions. Unexpectedly, removal of the adjacent −OH via −H substitution to give G(2′H,3′NH₃⁺) and G(2′NH₃⁺_,3′H) enhanced binding, in stark contrast to the deleterious effect of these substitutions on G binding. Pulse–chase experiments indicate that the −NH₃⁺ moiety of G(2′H,3′NH₃⁺) increases the rate of G association. These results suggest that the positively charged −NH₃⁺ group can act as a molecular “anchor” to increase the residence time of the encounter complex and thereby enhance productive binding. Electrostatic anchors may provide a broadly applicable strategy for the development of fast binding RNA ligands and RNA-targeted therapeutics.

Multi-metal-dependent nucleic acid enzymes

Nucleic acid enzymes require metal ions for activity, and many recently discovered enzymes can use multiple metals, either binding to the scissile phosphate or also playing an allosteric role.

Structure of the yeast spliceosomal postcatalytic P complex

The spliceosome undergoes dramatic changes in a splicing cycle. Structures of B, Bact, C, C*, and intron lariat spliceosome complexes revealed mechanisms of 5'-splice site (ss) recognition, branching, and intron release, but lacked information on 3'-ss recognition, exon ligation, and exon release. Here we report a cryo-electron microscopy structure of the postcatalytic P complex at 3.3-angstrom resolution, revealing that the 3' ss is mainly recognized through non-Watson-Crick base pairing with the 5' ss and branch point. Furthermore, one or more unidentified proteins become stably associated with the P complex, securing the 3' exon and potentially regulating activity of the helicase Prp22. Prp22 binds nucleotides 15 to 21 in the 3' exon, enabling it to pull the intron-exon or ligated exons in a 3' to 5' direction to achieve 3'-ss proofreading or exon release, respectively.

The Genetic Code and RNA-Amino Acid Affinities

A significant part of the genetic code likely originated via a chemical interaction, which should be experimentally verifiable. One possible verification relates bound amino acids (or perhaps their activated congeners) and ribonucleotide sequences within cognate RNA binding sites. To introduce this interaction, I first summarize how amino acids function as targets for RNA binding. Then the experimental method for selecting relevant RNA binding sites is characterized. The selection method’s characteristics are related to the investigation of the RNA binding site model treated at the outset. Finally, real binding sites from selection and also from extant natural RNAs (for example, the Sulfobacillus guanidinium riboswitch) are connected to the genetic code, and by extension, to the evolutionary progression that produced the code. During this process, peptides may have been produced directly on an instructive amino acid binding RNA (a DRT; Direct RNA Template). Combination of observed stereochemical selectivity with adaptation and co-evolutionary refinement is logically required, and also potentially sufficient, to create the striking order conserved throughout the present coding table.

Use of a Fluorescent Aptamer RNA as an Exonic Sequence to Analyze Self-Splicing Ability of aGroup I Intron from Structured RNAs

Group I self-splicing intron constitutes an important class of functional RNA molecules that can promote chemical transformation. Although the fundamental mechanism of the auto-excision from its precursor RNA has been established, convenient assay systems for its splicing activity are still useful for a further understanding of its detailed mechanism and of its application. Because some host RNA sequences, to which group I introns inserted form stable three-dimensional (3D) structures, the effects of the 3D structures of exonic elements on the splicing efficiency of group I introns are important but not a fully investigated issue. We developed an assay system for group I intron self-splicing by employing a fluorescent aptamer RNA (spinach RNA) as a model exonic sequence inserted by the Tetrahymena group I intron. We investigated self-splicing of the intron from spinach RNA, serving as a model exonic sequence with a 3D structure.

An active site rearrangement within the Tetrahymena group I ribozyme releases nonproductive interactions and allows formation of catalytic interactions

Biological catalysis hinges on the precise structural integrity of an active site that binds and transforms its substrates and meeting this requirement presents a unique challenge for RNA enzymes. Functional RNAs, including ribozymes, fold into their active conformations within rugged energy landscapes that often contain misfolded conformers. Here we uncover and characterize one such "off-pathway" species within an active site after overall folding of the ribozyme is complete. The Tetrahymena group I ribozyme (E) catalyzes cleavage of an oligonucleotide substrate (S) by an exogenous guanosine (G) cofactor. We tested whether specific catalytic interactions with G are present in the preceding E•S•G and E•G ground-state complexes. We monitored interactions with G via the effects of 2'- and 3'-deoxy (-H) and -amino (-NH(2)) substitutions on G binding. These and prior results reveal that G is bound in an inactive configuration within E•G, with the nucleophilic 3'-OH making a nonproductive interaction with an active site metal ion termed MA and with the adjacent 2'-OH making no interaction. Upon S binding, a rearrangement occurs that allows both -OH groups to contact a different active site metal ion, termed M(C), to make what are likely to be their catalytic interactions. The reactive phosphoryl group on S promotes this change, presumably by repositioning the metal ions with respect to G. This conformational transition demonstrates local rearrangements within an otherwise folded RNA, underscoring RNA's difficulty in specifying a unique conformation and highlighting Nature's potential to use local transitions of RNA in complex function.

A single mutation in the 15S rRNA gene confers non sense suppressor activity and interacts with mRF1 the release factor in yeast mitochondria

We have determined the nucleotide sequence of the mim3-1 mitochondrial ribosomal suppressor, acting on ochre mitochondrial mutations and one frameshift mutation in Saccharomyces cerevisiae. The 15s rRNA suppressor gene contains a G633 to C transversion. Yeast mitochondrial G633 corresponds to G517 of the E.coli 15S rRNA, which is occupied by an invariant G in all known small rRNA sequences. Interestingly, this mutation has occurred at the same position as the known MSU1 mitochondrial suppressor which changes G633 to A. The suppressor mutation lies in a highly conserved region of the rRNA, known in E.coli as the 530-loop, interacting with the S4, S5 and S12 ribosomal proteins. We also show an interesting interaction between the mitochondrial mim3-1 and the nuclear nam3-1 suppressors, both of which have the same action spectrum on mitochondrial mutations: nam3-1 abolishes the suppressor effect when present with mim3-1 in the same haploid cell. We discuss these results in the light of the nature of Nam3, identified by 1 as the yeast mitochondrial translation release factor. A hypothetical mechanism of suppression by "ribosome shifting" is also discussed in view of the nature of mutations suppressed and not suppressed.

FT-IR and Raman spectra, ab initio and density functional computations of the vibrational spectra, molecular geometries and atomic charges of uracil and 5-methyluracil (thymine)

FT-IR (400-4000 cm(-1)) and Raman spectra (200-4000 cm(-1)) of uracil and 5-methyluracil (thymine) have been recorded and analyzed. The optimized molecular geometries, atomic polar tensor (APT) charges and vibrational characteristics have been studied theoretically using restricted Hartree-Fock (RHF) and density functional theory (DFT) methods. Using the Becke's exchange in conjunction with Lee-Yang-Parr's correlation functional and Becke's three-parameter hybrid method (B3LYP), the ab initio and DFT calculations were carried out to study the optimized molecular fundamental vibrational frequencies for uracil and 5-methyluracil (thymine) by employing Gaussian-03 program. The fundamental vibrational frequencies along with their corresponding intensities in IR and Raman activities and depolarization ratios of the Raman lines have also been calculated using the RHF and DFT methods employing different basis sets. In quantum chemical calculations, most of the B3LYP/6-311++G(∗∗) vibrational frequencies are in excellent agreement with the available experimental assignments and helped to propose in the reassignments of some missing frequencies in experimental study. Assuming under the Cs point group for both molecules, the distribution of normal mode of vibrations between the two species as planar (a') and non-planar (a″) for all 39 normal vibrational modes of 5-methyluracil are given by 26a'+13a″, of which 30 modes (21a'+9a″) correspond to the uracil moiety and 9 modes (5a'+4a″) to the CH3 group. Consistent assignments have been made for the internal modes of CH3 group, especially for the anti-symmetric CH3 stretching and bending modes. A possible explanation could be the planarity of pyrimidine ring and non-planarity at carbon site of methyl group which might cause the splitting of frequencies including three components due to the substitution of CH3 group at the site of C5 atom on pyrimidine ring of uracil. The three non-equivalent CH bonds of CH3 group are distinctly separated from the CH/NH ring stretching frequencies. Kekule ring stretching mode is found to be comparatively higher frequency magnitude than the mode of uracil due to the involvement of hydrogen bonding of methyl group. But, the ring breathing is found to be lower frequency magnitude compared to those for uracil which could be due to mass effect of the CH3 group in place of the hydrogen atom at the site of C5 atom on pyrimidine ring of uracil.

IR and Raman spectra, ab initio and density functional computations of the vibrational spectra, molecular geometries and atomic charges of uracil and 5-aminouracil

Infrared (IR) and Raman spectra of uracil and 5-aminouracil have been recorded and analyzed between the region 200-4000 cm(-1). The optimized molecular geometries, atomic polar tensor (APT) charges and vibrational characteristics have been studied theoretically using restricted Hartree-Fock (RHF) and density functional theory (DFT) methods. Using the Becke's exchange in conjunction with Lee-Yang-Parr's correlation functional and Becke's three-parameter hybrid method (B3LYP), the ab initio and DFT calculations were carried out to study the optimized molecular fundamental vibrational frequencies for uracil and 5-aminouracil by employing Gaussian-03 program. The fundamental vibrational frequencies along with their corresponding intensities in IR and Raman activities and depolarization ratios of the Raman lines have also been calculated using the RHF and DFT methods employing different basis sets. In quantum chemical calculations, the most of B3LYP/6-311++G** vibrational frequencies are in the excellent agreement with available experimental assignments and helped in the reassignments of some fundamental vibrational modes. On the basis of calculated results, the assignments of some missing frequencies in the experimental study are proposed. Assuming under the Cs point group for both molecules, the distribution of normal mode of vibrations between the two species as planar (a') and non-planar (a″) are given by 25a'+11a″, of which 30 modes (21a'+9a″) correspond to the uracil moiety and 6 modes (4a'+2a″) to the NH2 group. Kekule ring stretching mode is found to be comparatively higher frequency magnitude than the mode of uracil due to the involvement of hydrogen bonding of amino group. But, the ring breathing is found to be lower frequency magnitude compared to those for uracil which could be due to mass effect of the NH2 group in place of the hydrogen atom. All other bands have also been assigned different fundamentals/overtones/combinations.

Spliceozymes: ribozymes that remove introns from pre-mRNAs in trans

Group I introns are pre-mRNA introns that do not require the spliceosome for their removal. Instead, they fold into complex three-dimensional structures and catalyze two transesterification reactions, thereby excising themselves and joining the flanking exons. These catalytic RNAs (ribozymes) have been modified previously to work in trans, whereby the ribozymes can recognize a splice site on a substrate RNA and replace the 5′- or 3′-portion of the substrate. Here we describe a new variant of the group I intron ribozyme from Tetrahymena that recognizes two splice sites on a substrate RNA, removes the intron sequences between the splice sites, and joins the flanking exons, analogous to the action of the spliceosome. This ‘group I spliceozyme’ functions in vitro and in vivo, and it is able to mediate a growth phenotype in E. coli cells. The intron sequences of the target pre-mRNAs are constrained near the splice sites but can carry a wide range of sequences in their interior. Because the splice site recognition sequences can be adjusted to different splice sites, the spliceozyme may have the potential for wide applications as tool in research and therapy.

FT-IR and Raman spectra, ab initio and density functional computations of the vibrational spectra, molecular geometries and atomic charges of uracil and 5-halogenated uracils (5-X-uracils; X=F, Cl, Br, I)

Raman (200-4000 cm(-1)) and FT-IR (400-4000 cm(-1)) spectra of uracil and 5-halogenated uracils (5-X-uracils; X=F, Cl, Br, I) have been recorded and analyzed in the range 200-4000 cm(-1). The optimized molecular geometries, atomic polar tensor (APT) charges and vibrational characteristics have been studied theoretically using restricted Hartree-Fock (RHF) and density functional theory (DFT) methods. Ab initio and DFT calculations [using Becke's exchange in conjunction with Lee-Yang-Parr's correlation functional and Becke's three-parameter hybrid method (B3LYP)] were carried out to study the optimized molecular fundamental vibrational frequencies for uracil and 5-halogenated uracils by employing Gaussian-03 program. Gauss View software was used to make the vibrational analysis. Raman and IR spectra have been computed theoretically for the uracil and 5-halogenated molecules. The fundamental vibrational frequencies along with their corresponding intensities in IR and Raman activities and depolarization ratios of the Raman lines have also been calculated using the RHF and DFT methods employing different basis sets. Quantum chemical calculations helped in the reassignments of some fundamental vibrational modes. Most of the B3LYP/6-311++G(**) vibrational frequencies are in excellent agreement with available experimental assignments. The ring breathing and kekule stretching modes are found to lower magnitudes compared to those for uracil which could be due to mass effect of halogen atom in place of the hydrogen atom. The C-X (X=F, Cl, Br, I) stretching frequency is distinctly separated from the CH/NH ring stretching frequencies on the pyrimidine ring. All other bands have also been assigned different fundamentals/overtones/combinations.

An RNA folding motif: GNRA tetraloop-receptor interactions

Nearly two decades after Westhof and Michel first proposed that RNA tetraloops may interact with distal helices, tetraloop–receptor interactions have been recognized as ubiquitous elements of RNA tertiary structure. The unique architecture of GNRA tetraloops (N=any nucleotide, R=purine) enables interaction with a variety of receptors, e.g., helical minor grooves and asymmetric internal loops. The most common example of the latter is the GAAA tetraloop–11 nt tetraloop receptor motif. Biophysical characterization of this motif provided evidence for the modularity of RNA structure, with applications spanning improved crystallization methods to RNA tectonics. In this review, we identify and compare types of GNRA tetraloop–receptor interactions. Then we explore the abundance of structural, kinetic, and thermodynamic information on the frequently occurring and most widely studied GAAA tetraloop–11 nt receptor motif. Studies of this interaction have revealed powerful paradigms for structural assembly of RNA, as well as providing new insights into the roles of cations, transition states and protein chaperones in RNA folding pathways. However, further research will clearly be necessary to characterize other tetraloop–receptor and long-range tertiary binding interactions in detail – an important milestone in the quantitative prediction of free energy landscapes for RNA folding.

Multiple group I introns in the small-subunit rDNA of Botryosphaeria dothidea: implication for intraspecific genetic diversity

Botryosphaeria dothidea is a widespread and economically important pathogen on various fruit trees, and it often causes die-back and canker on limbs and fruit rot. In characterizing intraspecies genetic variation within this fungus, group I introns, rich in rDNA of fungi, may provide a productive region for exploration. In this research, we analysed complete small subunit (SSU) ribosomal DNA (rDNA) sequences of 37 B. dothidea strains, and found four insertions, designated Bdo.S943, Bdo.S1199-A, Bdo.S1199-B and Bdo.S1506, at three positions. Sequence analysis and structure prediction revealed that both Bdo.S943 and Bdo.S1506 belonged to subgroup IC1 of group I introns, whereas Bdo.S1199-A and Bdo.S1199-B corresponded to group IE introns. Moreover, Bdo.S1199-A was found to host an open reading frame (ORF) for encoding the homing endonuclease (HE), whereas Bdo.S1199-B, an evolutionary descendant of Bdo.S1199-A, included a degenerate HE. The above four introns were novel, and were the first group I introns observed and characterized in this species. Differential distribution of these introns revealed that all strains could be separated into four genotypes. Genotype III (no intron) and genotype IV (Bdo.S1199-B) were each found in only one strain, whereas genotype I (Bdo.S1199-A) and genotype II (Bdo.S943 and Bdo.S1506) occurred in 95% of the strains. There is a correlation between B. dothidea genotypes and hosts or geographic locations. Thus, these newly discovered group I introns can help to advance understanding of genetic differentiation within B. dothidea.

The GA-minor submotif as a case study of RNA modularity, prediction, and design

Complex natural RNAs such as the ribosome, group I and group II introns, and RNase P exemplify the fact that three-dimensional (3D) RNA structures are highly modular and hierarchical in nature. Tertiary RNA folding typically takes advantage of a rather limited set of recurrent structural motifs that are responsible for controlling bends or stacks between adjacent helices. Herein, the GA minor and related structural motifs are presented as a case study to highlight several structural and folding principles, to gain further insight into the structural evolution of naturally occurring RNAs, as well as to assist the rational design of artificial RNAs.

Cooperative tertiary interaction network guides RNA folding

Noncoding RNAs form unique 3D structures, which perform many regulatory functions. To understand how RNAs fold uniquely despite a small number of tertiary interaction motifs, we mutated the major tertiary interactions in a group I ribozyme by single-base substitutions. The resulting perturbations to the folding energy landscape were measured using SAXS, ribozyme activity, hydroxyl radical footprinting, and native PAGE. Double- and triple-mutant cycles show that most tertiary interactions have a small effect on the stability of the native state. Instead, the formation of core and peripheral structural motifs is cooperatively linked in near-native folding intermediates, and this cooperativity depends on the native helix orientation. The emergence of a cooperative interaction network at an early stage of folding suppresses nonnative structures and guides the search for the native state. We suggest that cooperativity in noncoding RNAs arose from natural selection of architectures conducive to forming a unique, stable fold.

Exploring purine N7 interactions via atomic mutagenesis: the group I ribozyme as a case study

Atomic mutagenesis has emerged as a powerful tool to unravel specific interactions in complex RNA molecules. An early extensive study of analogs of the exogenous guanosine nucleophile in group I intron self-splicing by Bass and Cech demonstrated structure-function relationships analogous to those seen for protein ligands and provided strong evidence for a well-formed substrate binding site made of RNA. Subsequent functional and structural studies have confirmed these interacting sites and extended our understanding of them, with one notable exception. Whereas 7-methyl guanosine did not affect reactivity in the original study, a subsequent study revealed a deleterious effect of the seemingly more conservative 7-deaza substitution. Here we investigate this paradox, studying these and other analogs with the more thoroughly characterized ribozyme derived from the Tetrahymena group I intron. We found that the 7-deaza substitution lowers binding by ~20-fold, relative to the cognate exogenous guanosine nucleophile, whereas binding and reaction with 7-methyl and 8-aza-7-deaza substitutions have no effect. These and additional results suggest that there is no functionally important contact between the N7 atom of the exogenous guanosine and the ribozyme. Rather, they are consistent with indirect effects introduced by the N7 substitution on stacking interactions and/or solvation that are important for binding. The set of analogs used herein should be valuable in deciphering nucleic acid interactions and how they change through reaction cycles for other RNAs and RNA/protein complexes.

Site-specific isotope labeling of long RNA for structural and mechanistic studies

A site-specific isotope labeling technique of long RNA molecules was established. This technique is comprised of two simple enzymatic reactions, namely a guanosine transfer reaction of group I self-splicing introns and a ligation with T4 DNA ligase. The trans-acting group I self-splicing intron with its external cofactor, ‘isotopically labeled guanosine 5′-monophosphate’ (5′-GMP), steadily gave a 5′-residue-labeled RNA fragment. This key reaction, in combination with a ligation of 5′-remainder non-labeled sequence, allowed us to prepare a site-specifically labeled RNA molecule in a high yield, and its production was confirmed with ¹⁵N NMR spectroscopy. Such a site-specifically labeled RNA molecule can be used to detect a molecular interaction and to probe chemical features of catalytically/structurally important residues with NMR spectroscopy and possibly Raman spectroscopy and mass spectrometry.

Comprehensive survey and geometric classification of base triples in RNA structures

Base triples are recurrent clusters of three RNA nucleobases interacting edge-to-edge by hydrogen bonding. We find that the central base in almost all triples forms base pairs with the other two bases of the triple, providing a natural way to geometrically classify base triples. Given 12 geometric base pair families defined by the Leontis-Westhof nomenclature, combinatoric enumeration predicts 108 potential geometric base triple families. We searched representative atomic-resolution RNA 3D structures and found instances of 68 of the 108 predicted base triple families. Model building suggests that some of the remaining 40 families may be unlikely to form for steric reasons. We developed an on-line resource that provides exemplars of all base triples observed in the structure database and models for unobserved, predicted triples, grouped by triple family, as well as by three-base combination (http://rna.bgsu.edu/Triples). The classification helps to identify recurrent triple motifs that can substitute for each other while conserving RNA 3D structure, with applications in RNA 3D structure prediction and analysis of RNA sequence evolution.

Tightening of active site interactions en route to the transition state revealed by single-atom substitution in the guanosine-binding site of the Tetrahymena group I ribozyme

Protein enzymes establish intricate networks of interactions to bind and position substrates and catalytic groups within active sites, enabling stabilization of the chemical transition state. Crystal structures of several RNA enzymes also suggest extensive interaction networks, despite RNA's structural limitations, but there is little information on the functional and the energetic properties of these inferred networks. We used double mutant cycles and presteady-state kinetic analyses to probe the putative interaction between the exocyclic amino group of the guanosine nucleophile and the N7 atom of residue G264 of the Tetrahymena group I ribozyme. As expected, the results supported the presence of this interaction, but remarkably, the energetic penalty for introducing a CH group at the 7-position of residue G264 accumulates as the reaction proceeds toward the chemical transition state to a total of 6.2 kcal/mol. Functional tests of neighboring interactions revealed that the presence of the CH group compromises multiple contacts within the interaction network that encompass the reactive elements, apparently forcing the nucleophile to bind and attack from an altered, suboptimal orientation. The energetic consequences of this indirect disruption of neighboring interactions as the reaction proceeds demonstrate that linkage between binding interactions and catalysis hinges critically on the precise structural integrity of a network of interacting groups.

Occurrence and characteristics of group 1 introns found at three different positions within the 28S ribosomal RNA gene of the dematiaceous Phialophora verrucosa: phylogenetic and secondary structural implications

BACKGROUND

Group 1 introns (ribozymes) are among the most ancient and have the broadest phylogenetic distribution among the known self-splicing ribozymes. Fungi are known to be rich in rDNA group 1 introns. In the present study, five sequences of the 28S ribosomal RNA gene (rDNA) regions of pathogenic dematiaceous Phialophora verrucosa were analyzed using PCR by site-specific primers and were found to have three insertions, termed intron-F, G and H, at three positions of the gene. We investigated the distribution of group 1 introns in this fungus by surveying 34 strains of P. verrucosa and seven strains of Phialophora americana as the allied species.

RESULTS

Intron-F's (inserted at L798 position) were found in 88% of P. verrucosa strains, while intron-G's (inserted at L1921) at 12% and intron-H's (inserted at L2563) at 18%. There was some correlation between intron distribution and geographic location. In addition, we confirmed that the three kinds of introns are group 1 introns from results of BLAST search, alignment analysis and Reverse Transcriptase-Polymerase Chain Reaction (RT-PCR). Prediction of secondary structures and phylogenetic analysis of intron sequences identified introns-F and G as belonging to subgroup IC1. In addition, intron-H was identified as IE.

CONCLUSION

The three intron insertions and their insertion position in the 28S rDNA allowed the characterization of the clinical and environmental isolates of P. verrucosa and P. americana into five genotypes. All subgroups of introns-F and G and intron-H were characterized and observed for the first time in both species.

Predicting and modeling RNA architecture

A general approach for modeling the architecture of large and structured RNA molecules is described. The method exploits the modularity and the hierarchical folding of RNA architecture that is viewed as the assembly of preformed double-stranded helices defined by Watson-Crick base pairs and RNA modules maintained by non-Watson-Crick base pairs. Despite the extensive molecular neutrality observed in RNA structures, specificity in RNA folding is achieved through global constraints like lengths of helices, coaxiality of helical stacks, and structures adopted at the junctions of helices. The Assemble integrated suite of computer tools allows for sequence and structure analysis as well as interactive modeling by homology or ab initio assembly with possibilities for fitting within electronic density maps. The local key role of non-Watson-Crick pairs guides RNA architecture formation and offers metrics for assessing the accuracy of three-dimensional models in a more useful way than usual root mean square deviation (RMSD) values.

A folding "framework structure" of Tetrahymena group I intron

We have published the dynamic extended folding (DEF) method, which is a RNA secondary structure prediction approach-to simulate the in vivo RNA co-transcriptional folding process. In order to verify the reliability of the method, we selected the X-ray-determined Tetrahymena group I intron as a sample to construct the framework of its folding secondary structure. Our prediction coincides well with the secondary structure predicted by T.R. Cech and the X-ray diffraction crystal structure determined by Lehnert V. Our results show that the DEF framework structure of Tetrahymena group I intron reflects its function sites in a concise and straightforward manner, and the scope of the simulation was expanded.

A rearrangement of the guanosine-binding site establishes an extended network of functional interactions in the Tetrahymena group I ribozyme active site

Protein enzymes appear to use extensive packing and hydrogen bonding interactions to precisely position catalytic groups within active sites. Because of their inherent backbone flexibility and limited side chain repertoire, RNA enzymes face additional challenges relative to proteins in precisely positioning substrates and catalytic groups. Here, we use the group I ribozyme to probe the existence, establishment, and functional consequences of an extended network of interactions in an RNA active site. The group I ribozyme catalyzes a site-specific attack of guanosine on an oligonucleotide substrate. We previously determined that the hydrogen bond between the exocyclic amino group of guanosine and the 2'-hydroxyl group at position A261 of the Tetrahymena group I ribozyme contributes to overall catalysis. We now use functional data, aided by double mutant cycles, to probe this hydrogen bond in the individual reaction steps of the catalytic cycle. Our results indicate that this hydrogen bond is not formed upon guanosine binding to the ribozyme but instead forms at a later stage of the catalytic cycle. Formation of this hydrogen bond is correlated with other structural rearrangements in the ribozyme's active site that are promoted by docking of the oligonucleotide substrate into the ribozyme's active site, and disruption of this interaction has deleterious consequences for the chemical transformation within the ternary complex. These results, combined with earlier results, provide insight into the nature of the multiple conformational steps used by the Tetrahymena group I ribozyme to achieve its active structure and reveal an intricate, extended network of interactions that is used to establish catalytic interactions within this RNA's active site.

Exploring RNA structure by integrative molecular modelling

RNA molecular modelling is adequate to rapidly tackle the structure of RNA molecules. With new structured RNAs constituting a central class of cellular regulators discovered every year, the need for swift and reliable modelling methods is more crucial than ever. The pragmatic method based on interactive all-atom molecular modelling relies on the observation that specific structural motifs are recurrently found in RNA sequences. Once identified by a combination of comparative sequence analysis and biochemical data, the motifs composing the secondary structure of a given RNA can be extruded in three dimensions (3D) and used as building blocks assembled manually during a bioinformatic interactive process. Comparing the models to the corresponding crystal structures has validated the method as being powerful to predict the RNA topology and architecture while being less accurate regarding the prediction of base-base interactions. These aspects as well as the necessary steps towards automation will be discussed.

RNA reprogramming and repair based on trans-splicing group I ribozymes

While many traditional gene therapy strategies attempt to deliver new copies of wild-type genes back to cells harboring the defective genes, RNA-directed strategies offer a range of novel therapeutic applications. Revision or reprogramming of mRNA is a form of gene therapy that modifies mRNA without directly changing the transcriptional regulation or the genomic gene sequence. Group I ribozymes can be engineered to act in trans by recognizing a separate RNA molecule in a sequence-specific manner, and to covalently link a new RNA sequence to this separate RNA molecule. Group I ribozymes have been shown to repair defective transcripts that cause human genetic or malignant diseases, as well as to replace transcript sequences by foreign RNA resulting in new cellular functions. This review provides an overview of current strategies using trans-splicing group I ribozymes in RNA repair and reprogramming.

One-electron reduction of 2-aminopurine in the aqueous phase. A DFT and pulse radiolysis study

The electron affinity and the subsequent proton affinity of the electron adducts of 2-aminopurine (abbreviated 2AP) and adenine (A) are calculated with density functional theory (DFT). According to these calculations 2AP*- and A*- have similar thermochemical parameters leading to the conclusion that their reaction pathways should be close to analogous. Using the pulse radiolysis technique 2AP*- is formed by reaction with the hydrated electron (e(-)aq) and the resulting transient absorption spectrum is assigned to 2AP(NH)*. Additionally, it was found, employing the laser flash photolysis technique, that the excited singlet state of 2AP is incapable of oxidizing guanine in the aqueous phase. However, the one-electron oxidized 2AP (2AP*+) has sufficient energy to ionize guanine according to the DFT results in agreement with experimental data from the literature.

A unique group I intron in Coxiella burnetii is a natural splice mutant

Cbu.L1917, a group I intron present in the 23S rRNA gene of Coxiella burnetii, possesses a unique 3'-terminal adenine in place of a conserved guanine. Here, we show that, unlike all other group I introns, Cbu.L1917 utilizes a different cofactor for each splicing step and has a decreased self-splicing rate in vitro.

Coordination of two sequential ester-transfer reactions: exogenous guanosine binding promotes the subsequent omegaG binding to a group I intron

Self-splicing of group I introns is accomplished by two sequential ester-transfer reactions mediated by sequential binding of two different guanosine ligands, but it is yet unclear how the binding is coordinated at a single G-binding site. Using a three-piece trans-splicing system derived from the Candida intron, we studied the effect of the prior GTP binding on the later ωG binding by assaying the ribozyme activity in the second reaction. We showed that adding GTP simultaneously with and prior to the esterified ωG in a substrate strongly accelerated the second reaction, suggesting that the early binding of GTP facilitates the subsequent binding of ωG. GTP-mediated facilitation requires C2 amino and C6 carbonyl groups on the Watson–Crick edge of the base but not the phosphate or sugar groups, suggesting that the base triple interactions between GTP and the binding site are important for the subsequent ωG binding. Strikingly, GTP binding loosens a few local structures of the ribozyme including that adjacent to the base triple, providing structural basis for a rapid exchange of ωG for bound GTP.

Splicing and evolution of an unusually small group I intron

Introns are common in the rRNA gene loci of fungal genomes, but biochemical studies to investigate splicing are rare. Here, self-splicing of a very small (67 nucleotide) group I intron is demonstrated. The PaSSU intron (located within the rRNA small subunit gene of Phialophora americana) splices in vitro under group I intron conditions. Most group I ribozymes contain pairing regions P1-P10, with a conserved G.U pair at the 5' splice site, and a G at the 3' intron border. The PaSSU intron contains only P1, P7, and P10. While it contains the G.U pair at the 5' splice, a U is found at the 3' end of the intron instead of a G. Phylogenetic analysis places it within subgroup IC1, whose members are found in the nuclear rRNA genes of fungi. The structural elements are similar to those in the centermost regions of other group I introns. Its size can be explained by a single large deletion that removed P2 through much of P9. Part of the original P9 region has assumed the function of P7. Its small size and genealogy makes it an excellent model to study RNA catalysis and evolution.

The 2'-hydroxyl group of the guanosine nucleophile donates a functionally important hydrogen bond in the tetrahymena ribozyme reaction

In the first step of self-splicing, group I introns utilize an exogenous guanosine nucleophile to attack the 5'-splice site. Removal of the 2'-hydroxyl of this guanosine results in a 10 (6)-fold loss in activity, indicating that this functional group plays a critical role in catalysis. Biochemical and structural data have shown that this hydroxyl group provides a ligand for one of the catalytic metal ions at the active site. However, whether this hydroxyl group also engages in hydrogen-bonding interactions remains unclear, as attempts to elaborate its function further usually disrupt the interactions with the catalytic metal ion. To address the possibility that this 2'-hydroxyl contributes to catalysis by donating a hydrogen bond, we have used an atomic mutation cycle to probe the functional importance of the guanosine 2'-hydroxyl hydrogen atom. This analysis indicates that, beyond its role as a ligand for a catalytic metal ion, the guanosine 2'-hydroxyl group donates a hydrogen bond in both the ground state and the transition state, thereby contributing to cofactor recognition and catalysis by the intron. Our findings continue an emerging theme in group I intron catalysis: the oxygen atoms at the reaction center form multidentate interactions that function as a cooperative network. The ability to delineate such networks represents a key step in dissecting the complex relationship between RNA structure and catalysis.

Kinetic characterization of the first step of the ribozyme-catalyzed trans excision-splicing reaction

Group I introns catalyze the self-splicing reaction, and their derived ribozymes are frequently used as model systems for the study of RNA folding and catalysis, as well as for the development of non-native catalytic reactions. Utilizing a group I intron-derived ribozyme from Pneumocystis carinii, we previously reported a non-native reaction termed trans excision-splicing (TES). In this reaction, an internal segment of RNA is excised from an RNA substrate, resulting in the covalent reattachment of the flanking regions. TES proceeds through two consecutive phosphotransesterification reactions, which are similar to the reaction steps of self-splicing. One key difference is that TES utilizes the 3'-terminal guanosine of the ribozyme as the first-step nucleophile, whereas self-splicing utilizes an exogenous guanosine. To further aid in our understanding of ribozyme reactions, a kinetic framework for the first reaction step (substrate cleavage) was established. The results demonstrate that the substrate binds to the ribozyme at a rate expected for simple helix formation. In addition, the rate constant for the first step of the TES reaction is more than one order of magnitude lower than the analogous step in self-splicing. Results also suggest that a conformational change, likely similar to that in self-splicing, exists between the two reaction steps of TES. Finally, multiple turnover is curtailed because dissociation of the cleavage product is slower than the rate of chemistry.

A Pneumocystis carinii group I intron-derived ribozyme utilizes an endogenous guanosine as the first reaction step nucleophile in the trans excision-splicing reaction

In the trans excision-splicing reaction, a Pneumocystis carinii group I intron-derived ribozyme binds an RNA substrate, excises a specific internal segment, and ligates the flanking regions back together. This reaction can occur both in vitro and in vivo. In this report, the first of the two reaction steps was analyzed to distinguish between two reaction mechanisms: ribozyme-mediated hydrolysis and nucleotide-dependent intramolecular transesterification. We found that the 3'-terminal nucleotide of the ribozyme is the first-reaction step nucleophile. In addition, the 3'-half of the RNA substrate becomes covalently attached to the 3'-terminal nucleotide of the ribozyme during the reaction, both in vitro and in vivo. Results also show that the identity of the 3'-terminal nucleotide influences the rate of the intramolecular transesterification reaction, with guanosine being more effective than adenosine. Finally, expected products of the hydrolysis mechanism do not form during the reaction. These results are consistent with only the intramolecular transesterification mechanism. Unexpectedly, we also found that ribozyme constructs become truncated in vivo, probably through intramolecular 3'-hydrolysis (self-activation), to create functional 3'-terminal nucleotides.

Structural inference of native and partially folded RNA by high-throughput contact mapping

The biological behaviors of ribozymes, riboswitches, and numerous other functional RNA molecules are critically dependent on their tertiary folding and their ability to sample multiple functional states. The conformational heterogeneity and partially folded nature of most of these states has rendered their characterization by high-resolution structural approaches difficult or even intractable. Here we introduce a method to rapidly infer the tertiary helical arrangements of large RNA molecules in their native and non-native solution states. Multiplexed hydroxyl radical (.OH) cleavage analysis (MOHCA) enables the high-throughput detection of numerous pairs of contacting residues via random incorporation of radical cleavage agents followed by two-dimensional gel electrophoresis. We validated this technology by recapitulating the unfolded and native states of a well studied model RNA, the P4-P6 domain of the Tetrahymena ribozyme, at subhelical resolution. We then applied MOHCA to a recently discovered third state of the P4-P6 RNA that is stabilized by high concentrations of monovalent salt and whose partial order precludes conventional techniques for structure determination. The three-dimensional portrait of a compact, non-native RNA state reveals a well ordered subset of native tertiary contacts, in contrast to the dynamic but otherwise similar molten globule states of proteins. With its applicability to nearly any solution state, we expect MOHCA to be a powerful tool for illuminating the many functional structures of large RNA molecules and RNA/protein complexes.