Dynamic behavior of the telomerase RNA hairpin structure and its relationship to dyskeratosis congenita

Effects of flanking regions on HDV cotranscriptional folding kinetics

Hepatitis delta virus (HDV) ribozyme performs the self-cleavage activity through folding to a double pseudoknot structure. The folding of functional RNA structures is often coupled with the transcription process. In this work, we developed a new approach for predicting the cotranscriptional folding kinetics of RNA secondary structures with pseudoknots. We theoretically studied the cotranscriptional folding behavior of the 99-nucleotide (nt) HDV sequence, two upstream flanking sequences, and one downstream flanking sequence. During transcription, the 99-nt HDV can effectively avoid the trap intermediates and quickly fold to the cleavage-active state. It is different from its refolding kinetics, which folds into an intermediate trap state. For all the sequences, the ribozyme regions (from 1 to 73) all fold to the same structure during transcription. However, the existence of the 30-nt upstream flanking sequence can inhibit the ribozyme region folding into the active native state through forming an alternative helix Alt1 with the segments 70-90. The longer upstream flanking sequence of 54 nt itself forms a stable hairpin structure, which sequesters the formation of the Alt1 helix and leads to rapid formation of the cleavage-active structure. Although the 55-nt downstream flanking sequence could invade the already folded active structure during transcription by forming a more stable helix with the ribozyme region, the slow transition rate could keep the structure in the cleavage-active structure to perform the activity.

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview

With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.

Quantitative sampling of conformational heterogeneity of a DNA hairpin using molecular dynamics simulations and ultrafast fluorescence spectroscopy

Molecular dynamics (MD) simulations and time resolved fluorescence (TRF) spectroscopy were combined to quantitatively describe the conformational landscape of the DNA primary binding sequence (PBS) of the HIV-1 genome, a short hairpin targeted by retroviral nucleocapsid proteins implicated in the viral reverse transcription. Three 2-aminopurine (2AP) labeled PBS constructs were studied. For each variant, the complete distribution of fluorescence lifetimes covering 5 orders of magnitude in timescale was measured and the populations of conformers experimentally observed to undergo static quenching were quantified. A binary quantification permitted the comparison of populations from experimental lifetime amplitudes to populations of aromatically stacked 2AP conformers obtained from simulation. Both populations agreed well, supporting the general assumption that quenching of 2AP fluorescence results from pi-stacking interactions with neighboring nucleobases and demonstrating the success of the proposed methodology for the combined analysis of TRF and MD data. Cluster analysis of the latter further identified predominant conformations that were consistent with the fluorescence decay times and amplitudes, providing a structure-based rationalization for the wide range of fluorescence lifetimes. Finally, the simulations provided evidence of local structural perturbations induced by 2AP. The approach presented is a general tool to investigate fine structural heterogeneity in nucleic acid and nucleoprotein assemblies.

Structure Prediction: New Insights into Decrypting Long Noncoding RNAs

Long noncoding RNAs (lncRNAs), which form a diverse class of RNAs, remain the least understood type of noncoding RNAs in terms of their nature and identification. Emerging evidence has revealed that a small number of newly discovered lncRNAs perform important and complex biological functions such as dosage compensation, chromatin regulation, genomic imprinting, and nuclear organization. However, understanding the wide range of functions of lncRNAs related to various processes of cellular networks remains a great experimental challenge. Structural versatility is critical for RNAs to perform various functions and provides new insights into probing the functions of lncRNAs. In recent years, the computational method of RNA structure prediction has been developed to analyze the structure of lncRNAs. This novel methodology has provided basic but indispensable information for the rapid, large-scale and in-depth research of lncRNAs. This review focuses on mainstream RNA structure prediction methods at the secondary and tertiary levels to offer an additional approach to investigating the functions of lncRNAs.

Understanding the kinetic mechanism of RNA single base pair formation

RNA functions are intrinsically tied to folding kinetics. The most elementary step in RNA folding is the closing and opening of a base pair. Understanding this elementary rate process is the basis for RNA folding kinetics studies. Previous studies mostly focused on the unfolding of base pairs. Here, based on a hybrid approach, we investigate the folding process at level of single base pairing/stacking. The study, which integrates molecular dynamics simulation, kinetic Monte Carlo simulation, and master equation methods, uncovers two alternative dominant pathways: Starting from the unfolded state, the nucleotide backbone first folds to the native conformation, followed by subsequent adjustment of the base conformation. During the base conformational rearrangement, the backbone either retains the native conformation or switches to nonnative conformations in order to lower the kinetic barrier for base rearrangement. The method enables quantification of kinetic partitioning among the different pathways. Moreover, the simulation reveals several intriguing ion binding/dissociation signatures for the conformational changes. Our approach may be useful for developing a base pair opening/closing rate model.

A U⋅U Pair-to-U⋅C Pair Mutation-Induced RNA Native Structure Destabilisation and Stretching-Force-Induced RNA Misfolding

Little is known about how a non-Watson-Crick pair affects the RNA folding dynamics. We studied the effects of a U⋅U-to-U⋅C pair mutation on the folding of a hairpin in human telomerase RNA. The ensemble thermal melting of the hairpins shows an on-pathway intermediate with the disruption of the internal loop structure containing the U⋅U/U⋅C pairs. By using optical tweezers, we applied a stretching force on the terminal ends of the hairpins to probe directly the non-nearest-neighbour effects upon the mutations. The single U⋅U to U⋅C mutations are observed to 1) lower the mechanical unfolding force by approximately 1 picoNewton (pN) per mutation without affecting the unfolding reaction transition-state position (thus suggesting that removing a single hydrogen bond affects the structural dynamics at least two base pairs away), 2) result in more frequent misfolding into a small hairpin at approximately 10 pN and 3) shift the folding reaction transition-state position towards the native hairpin structure and slightly increase the mechanical folding kinetics (thus suggesting that untrapping from the misfolded state is not the rate-limiting step).

Large-scale study of long non-coding RNA functions based on structure and expression features

Mammals and other complex organisms can transcribe an abundance of long non-coding RNAs (lncRNAs) that fulfill a wide variety of regulatory roles in many biological processes. These roles, including as scaffolds and as guides for protein-coding genes, mainly depend on the structure and expression level of lncRNAs. In this review, we focus on the current methods for analyzing lncRNA structure and expression, which is basic but necessary information for in-depth, large-scale analysis of lncRNA functions.

Interactions of cations with RNA loop-loop complexes

RNA loop-loop interactions are essential in many biological processes, including initiation of RNA folding into complex tertiary shapes, promotion of dimerization, and viral replication. In this article, we examine interactions of metal ions with five RNA loop-loop complexes of unique biological significance using explicit-solvent molecular-dynamics simulations. These simulations revealed the presence of solvent-accessible tunnels through the major groove of loop-loop interactions that attract and retain cations. Ion dynamics inside these loop-loop complexes were distinctly different from the dynamics of the counterion cloud surrounding RNA and depend on the number of basepairs between loops, purine sequence symmetry, and presence of unpaired nucleotides. The cationic uptake by kissing loops depends on the number of basepairs between loops. It is interesting that loop-loop complexes with similar functionality showed similarities in cation dynamics despite differences in sequence and loop size.

RNA stability under different combinations of amber force fields and solvation models

The proper matching of force field and solvent is critical to obtain correct result in molecular dynamics simulation of bio-molecules. This problem has been intensively investigated for protein but not for RNA yet. In this paper, we use standard molecular dynamics and replica exchange molecular dynamics to take a series of tests on the RNA stability under different combinations of Amber force field parameters (ff98, ff99 and ff99bsc0) and the general Born implicit solvent models (igb1, igb2 and igb5). It is found that only ff98 and ff99bsc0 with igb1 can keep the native conformations of RNA hairpin and duplex. Our results suggest that ff98 plus igb1 may be reasonable choice for molecular dynamics simulation of RNA dynamics.

Structural dynamics of the box C/D RNA kink-turn and its complex with proteins: the role of the A-minor 0 interaction, long-residency water bridges, and structural ion-binding sites revealed by molecular simulations

Kink-turns (K-turns) are recurrent elbow-like RNA motifs that participate in protein-assisted RNA folding and contribute to RNA dynamics. We carried out a set of molecular dynamics (MD) simulations using parm99 and parmbsc0 force fields to investigate structural dynamics of the box C/D RNA and its complexes with two proteins: native archaeal L7ae protein and human 15.5 kDa protein, originally bound to very similar structure of U4 snRNA. The box C/D RNA forms K-turn with A-minor 0 tertiary interaction between its canonical (C) and noncanonical (NC) stems. The local K-turn architecture is thus different from the previously studied ribosomal K-turns 38 and 42 having A-minor I interaction. The simulations reveal visible structural dynamics of this tertiary interaction involving altogether six substates which substantially contribute to the elbow-like flexibility of the K-turn. The interaction can even temporarily shift to the A-minor I type pattern; however, this is associated with distortion of the G/A base pair in the NC-stem of the K-turn. The simulations show reduction of the K-turn flexibility upon protein binding. The protein interacts with the apex of the K-turn and with the NC-stem. The protein-RNA interface includes long-residency hydration sites. We have also found long-residency hydration sites and major ion-binding sites associated with the K-turn itself. The overall topology of the K-turn remains stable in all simulations. However, in simulations of free K-turn, we observed instability of the key C16(O2')-A7(N1) H-bond, which is a signature interaction of K-turns and which was visibly more stable in simulations of K-turns possessing A-minor I interaction. It may reflect either some imbalance of the force field or it may be a correct indication of early stages of unfolding since this K-turn requires protein binding for its stabilization. Interestingly, the 16(O2')-7(N1) H- bond is usually not fully lost since it is replaced by a water bridge with a tightly bound water, which is adenine-specific similarly as the original interaction. The 16(O2')-7(N1) H-bond is stabilized by protein binding. The stabilizing effect is more visible with the human 15.5 kDa protein, which is attributed to valine to arginine substitution in the binding site. The behavior of the A-minor interaction is force-field-dependent because the parmbsc0 force field attenuates the A-minor fluctuations compared to parm99 simulations. Behavior of other regions of the box C/D RNA is not sensitive to the force field choice. Simulation with net-neutralizing Na(+) and 0.2 M excess salt conditions appear in all aspects equivalent. The simulations show loss of a hairpin tetraloop, which is not part of the K-turn. This was attributed to force field limitations.

Effect of oligonucleotide length on the assembly of DNA materials: molecular dynamics simulations of layer-by-layer DNA films

DNA strand length has been found to be an important factor in many DNA-based nanoscale systems. Here, we apply molecular dynamics simulations in a synergistic effort with layer-by-layer experimental data to understand the effect of DNA strand length on the assembly of DNA films. The results indicate that short (less than 10 bases) and long (more than 30 bases) single-stranded DNAs do not exhibit optimal film growth, and this can be associated with the limited accessibility of the bases on the surface due to formation of self-protected interactions that prevent efficient hybridization. Interestingly, the presence of a duplex attached to a single strand significantly alters the persistence length of the polyT strands. Our study suggests that restrained polyT, compared to labile suspensions of free polyT, are more capable of hybridization and hence DNA-based assembly.

A CA(+) pair adjacent to a sheared GA or AA pair stabilizes size-symmetric RNA internal loops

RNA internal loops are often important sites for folding and function. Residues in internal loops can have pK_a values shifted close to neutral pH because of the local structural environment. A series of RNA internal loops were studied at different pH by UV absorbance versus temperature melting experiments and imino proton nuclear magnetic resonance (NMR). A stabilizing CA pair forms at pH 7 in the and nearest neighbors when the CA pair is the first noncanonical pair (loop-terminal pair) in 3 × 3 nucleotide and larger size-symmetric internal loops. These and nearest neighbors, with CA adjacent to a closing Watson−Crick pair, are further stabilized when the pH is lowered from 7 to 5.5. The results are consistent with a significantly larger fraction (from ∼20% at pH 7 to ∼90% at pH 5.5) of adenines being protonated at the N1 position to form stabilizing wobble CA⁺ pairs adjacent to a sheared GA or AA pair. The noncanonical pair adjacent to the GA pair in can either stabilize or destabilize the loop, consistent with the sequence-dependent thermodynamics of GA pairs. No significant pH-dependent stabilization is found for most of the other nearest neighbor combinations involving CA pairs (e.g., and ), which is consistent with the formation of various nonwobble pairs observed in different local sequence contexts in crystal and NMR structures. A revised free-energy model, including stabilization by wobble CA⁺ pairs, is derived for predicting stabilities of medium-size RNA internal loops.

Triplex structures in an RNA pseudoknot enhance mechanical stability and increase efficiency of -1 ribosomal frameshifting

Many viruses use programmed -1 ribosomal frameshifting to express defined ratios of structural and enzymatic proteins. Pseudoknot structures in messenger RNAs stimulate frameshifting in upstream slippery sequences. The detailed molecular determinants of pseudoknot mechanical stability and frameshifting efficiency are not well understood. Here we use single-molecule unfolding studies by optical tweezers, and frameshifting assays to elucidate how mechanical stability of a pseudoknot and its frameshifting efficiency are regulated by tertiary stem-loop interactions. Mechanical unfolding of a model pseudoknot and mutants designed to dissect specific interactions reveals that mechanical stability depends strongly on triplex structures formed by stem-loop interactions. Combining single-molecule and mutational studies facilitates the identification of pseudoknot folding intermediates. Average unfolding forces of the pseudoknot and mutants ranging from 50 to 22 picoNewtons correlated with frameshifting efficiencies ranging from 53% to 0%. Formation of major-groove and minor-groove triplex structures enhances pseudoknot stem stability and torsional resistance, and may thereby stimulate frameshifting. Better understanding of the molecular determinants of frameshifting efficiency may facilitate the development of anti-virus therapeutics targeting frameshifting.

Discovery of ligands for a novel target, the human telomerase RNA, based on flexible-target virtual screening and NMR

The human ribonucleoprotein telomerase is a validated anticancer drug target, and hTR-P2b is a part of the human telomerase RNA (hTR) essential for its activity. Interesting ligands that bind hTR-P2b were identified by iteratively using a tandem structure-based approach: docking of potential ligands from small databases to hTR-P2b via the program MORDOR, which permits flexibility in both ligand and target, with subsequent NMR screening of high-ranking compounds. A high percentage of the compounds tested experimentally were found via NMR to bind to the U-rich region of hTR-P2b; most have MW < 500 Da and are from different compound classes, and several possess a charge of 0 or +1. Of the 48 ligands identified, 24 exhibit a decided preference to bind hTR-P2b RNA rather than A-site rRNA and 10 do not bind A-site rRNA at all. Binding affinity was measured by monitoring RNA imino proton resonances for some of the compounds that showed hTR binding preference.

Theoretical studies on the interaction of modified pyrimidines and purines with purine riboswitch

Recent experimental study [S.D. Gilbert, S.J. Mediatore, R.T. Batey, Modified pyrimidine specifically bind the purine riboswitch, J. Am. Chem. Soc. 128 (2006) 14214-14215] demonstrated that the purine riboswitch could specifically bind some ligands other than purines such as amino-pyrimidines, and the authors proposed that the five-membered ring of purine was not required for recognition. To get insight into the interaction details, we used molecular docking method to investigate the interactions of a mutant form of guanine riboswitch with a series of amino-purines, amino-pyrimidines and imidazole derivatives, and employed molecular simulation method to study the dynamic behavior of the selected complexes. The calculation results reveal that (1) all the amino-purines and amino-pyrimidines bind in a same cavity composed of four nucleobases including U22, U47, U51 and U74, which is consistent with the experimental results, while the two imidazole derivatives adopt other binding modes; (2) the purines are engulfed within three-way junction motifs, but most pyrimidines only form two-way junctions with the riboswitch; (3) the number and position of amino substituents could seriously affect the binding of pyrimidines. As riboswitches are potentially excellent candidates for antibiotic therapeutics, these findings may be useful for understanding the range of compounds that riboswitch can specifically recognize.

NMR-assisted prediction of RNA secondary structure: identification of a probable pseudoknot in the coding region of an R2 retrotransposon

As the rate of functional RNA sequence discovery escalates, high-throughput techniques for reliable structural determination are becoming crucial for revealing the essential features of these RNAs in a timely fashion. Computational predictions of RNA secondary structure quickly generate reasonable models but suffer from several approximations, including overly simplified models and incomplete knowledge of significant interactions. Similar problems limit the accuracy of predictions for other self-folding polymers, including DNA and peptide nucleic acid (PNA). The work presented here demonstrates that incorporating unassigned data from simple nuclear magnetic resonance (NMR) experiments into a dynamic folding algorithm greatly reduces the potential folding space of a given RNA and therefore increases the confidence and accuracy of modeling. This procedure has been packaged into an NMR-assisted prediction of secondary structure (NAPSS) algorithm that can produce pseudoknotted as well as non-pseudoknotted secondary structures. The method reveals a probable pseudoknot in the part of the coding region of the R2 retrotransposon from Bombyx mori that orchestrates second-strand DNA cleavage during insertion into the genome.

Protocols for the in silico design of RNA nanostructures

Recent developments in the field of nanobiology have significantly expanded the possibilities for new modalities in the treatment of many diseases, including cancer. Ribonucleic acid (RNA) represents a relatively new molecular material for the development of these biologically oriented nanodevices. In addition, RNA nanobiology presents a relatively new approach for the development of RNA-based nanoparticles that can be used as crystallization substrates and scaffolds for RNA-based nanoarrays. Presented in this chapter are some methodological shaped-based protocols for the design of such RNA nanostructures. Included are descriptions and background materials describing protocols that use a database of three-dimensional RNA structure motifs; designed RNA secondary structure motifs; and a combination of the two approaches. An example is also given illustrating one of the protocols.

Single-molecule mechanical unfolding and folding of a pseudoknot in human telomerase RNA

RNA unfolding and folding reactions in physiological conditions can be facilitated by mechanical force one molecule at a time. By using force-measuring optical tweezers, we studied the mechanical unfolding and folding of a hairpin-type pseudoknot in human telomerase RNA in a near-physiological solution, and at room temperature. Discrete two-state folding transitions of the pseudoknot are seen at approximately 10 and approximately 5 piconewtons (pN), with ensemble rate constants of approximately 0.1 sec(-1), by stepwise force-drop experiments. Folding studies of the isolated 5'-hairpin construct suggested that the 5'-hairpin within the pseudoknot forms first, followed by formation of the 3'-stem. Stepwise formation of the pseudoknot structure at low forces are in contrast with the one-step unfolding at high forces of approximately 46 pN, at an average rate of approximately 0.05 sec(-1). In the constant-force folding trajectories at approximately 10 pN and approximately 5 pN, transient formation of nonnative structures were observed, which is direct experimental evidence that folding of both the hairpin and pseudoknot takes complex pathways. Possible nonnative structures and folding pathways are discussed.

The impact of dyskeratosis congenita mutations on the structure and dynamics of the human telomerase RNA pseudoknot domain

The pseudoknot domain is a functionally crucial part of telomerase RNA and influences the activity and stability of the ribonucleoprotein complex. Autosomal dominant dyskeratosis congenita (DKC) is an inherited disease that is linked to mutations in telomerase RNA and impairs telomerase function. In this paper, we present a computational prediction of the influence of two base DKC mutations on the structure, dynamics, and stability of the pseudoknot domain. We use molecular dynamics simulations, MM-GBSA free energy calculations, static analysis, and melting simulations analysis. Our results show that the DKC mutations stabilize the hairpin form and destabilize the pseudoknot form of telomerase RNA. Moreover, the P3 region of the predicted DKC-mutated pseudoknot structure is unstable and fails to form as a defined helical stem. We directly compare our predictions with experimental observations by calculating the enthalpy of folding and melting profiles for each structure. The enthalpy values are in very good agreement with values determined by thermal denaturation experiments. The melting simulations and simulations at elevated temperatures show the existence of an intermediate structure, which involves the formation of two UU base pairs observed in the hairpin form of the pseudoknot domain.

Biphasic folding kinetics of RNA pseudoknots and telomerase RNA activity

Using a combined master equation and kinetic cluster approach, we investigate RNA pseudoknot folding and unfolding kinetics. The energetic parameters are computed from a recently developed Vfold model for RNA secondary structure and pseudoknot folding thermodynamics. The folding kinetics theory is based on the complete conformational ensemble, including all the native-like and non-native states. The predicted folding and unfolding pathways, activation barriers, Arrhenius plots, and rate-limiting steps lead to several findings. First, for the PK5 pseudoknot, a misfolded 5' hairpin emerges as a stable kinetic trap in the folding process, and the detrapping from this misfolded state is the rate-limiting step for the overall folding process. The calculated rate constant and activation barrier agree well with the experimental data. Second, as an application of the model, we investigate the kinetic folding pathways for human telomerase RNA (hTR) pseudoknot. The predicted folding and unfolding pathways not only support the proposed role of conformational switch between hairpin and pseudoknot in hTR activity, but also reveal molecular mechanism for the conformational switch. Furthermore, for an experimentally studied hTR mutation, whose hairpin intermediate is destabilized, the model predicts a long-lived transient hairpin structure, and the switch between the transient hairpin intermediate and the native pseudoknot may be responsible for the observed hTR activity. Such finding would help resolve the apparent contradiction between the observed hTR activity and the absence of a stable hairpin.

Molecular dynamics simulations of RNA: an in silico single molecule approach

RNA molecules are now known to be involved in the processing of genetic information at all levels, taking on a wide variety of central roles in the cell. Understanding how RNA molecules carry out their biological functions will require an understanding of structure and dynamics at the atomistic level, which can be significantly improved by combining computational simulation with experiment. This review provides a critical survey of the state of molecular dynamics (MD) simulations of RNA, including a discussion of important current limitations of the technique and examples of its successful application. Several types of simulations are discussed in detail, including those of structured RNA molecules and their interactions with the surrounding solvent and ions, catalytic RNAs, and RNA-small molecule and RNA-protein complexes. Increased cooperation between theorists and experimentalists will allow expanded judicious use of MD simulations to complement conceptually related single molecule experiments. Such cooperation will open the door to a fundamental understanding of the structure-function relationships in diverse and complex RNA molecules. .

Conformational analysis of the telomerase RNA pseudoknot hairpin by Raman spectroscopy

We have measured the temperature-dependent Raman spectra of two 30-mer ribonucleotides that represent the wild-type (WT) and dyskeratosis congenita (DKC) mutant (MT) GC (107-108) --> AG structures of the pseudoknot hairpin region of human telomerase RNA. We have used these structures, previously characterized by UV-melting and NMR, as a model system for our Raman investigation. We observe that Raman hypochromism of vibrational bands, previously assigned to specific bases or conformational RNA markers, reflect temperature-dependent alterations in the pentaloop and stem structures of these two oligonucleotides. We also observe that the intense nu(s)(O-P-O) band at 812 cm(-1) indicates the presence of A-form backbone structure at relatively low temperatures in both the WT and MT RNA sequences. The mutation induces a decrease in the intensity of the uridine (rU) band at 1244 cm(-1) associated with C2'-endo/anti ribose conformation in the pentaloop. Two transition temperatures (T(m) ) were determined from the analysis of Raman difference intensity-temperature profiles of the 1256 cm(-1) band, which is associated with vibrations of cytidine (rC) residues, in particular, the C2'-endo/anti ribose conformation (T(m) 1 = 23.6 +/- 1.6 degrees C for WT and 19.7 +/- 2.8 degrees C for MT; T(m) 2 = 68.9 +/- 1.8 degrees C for WT and 70.9 +/- 1.1 degrees C for MT). From these results we can conclude that the DKC mutant 30-mer exhibits a lower stability in the pentaloop region and a slightly higher stability in the stem region than the WT 30-mer. This demonstrates that Raman bands, previously assigned to specific bases or conformational RNA markers, can be used to probe local structural features of the telomerase pseudoknot hairpin sequence.

High sensitivity RNA pseudoknot prediction

Most ab initio pseudoknot predicting methods provide very few folding scenarios for a given RNA sequence and have low sensitivities. RNA researchers, in many cases, would rather sacrifice the specificity for a much higher sensitivity for pseudoknot detection. In this study, we introduce the Pseudoknot Local Motif Model and Dynamic Partner Sequence Stacking (PLMM_DPSS) algorithm which predicts all PLM model pseudoknots within an RNA sequence in a neighboring-region-interference-free fashion. The PLM model is derived from the existing Pseudobase entries. The innovative DPSS approach calculates the optimally lowest stacking energy between two partner sequences. Combined with the Mfold, PLMM_DPSS can also be used in predicting complicated pseudoknots. The test results of PLMM_DPSS, PKNOTS, iterated loop matching, pknotsRG and HotKnots with Pseudobase sequences have shown that PLMM_DPSS is the most sensitive among the five methods. PLMM_DPSS also provides manageable pseudoknot folding scenarios for further structure determination.

The prediction of the wild-type telomerase RNA pseudoknot structure and the pivotal role of the bulge in its formation

In this study, the three-dimensional structure of the wild-type human telomerase RNA pseudoknot was predicted via molecular modeling. The wild-type pseudoknot structure is then compared to the recent NMR solution structure of the telomerase pseudoknot, which does not contain the U177 bulge. The removal of the bulge from the pseudoknot structure results in higher stability and significant reduction of activity of telomerase. We show that the effect of the bulge on the structure results in a significant transformation of the pseudoknot junction region where the starting base pairs are disrupted and unique triple base pairs are formed. We found that the formation of the junction region is greatly influenced by interactions of the U177 bulge with loop residues and rotation of residue A174. Moreover, this is the first study to our knowledge where a structure as complex as the pseudoknot has been solved by purely theoretical methods.

Dyskeratosis Congenita: A Disorder of Defective Telomere Maintenance?

Dyskeratosis congenita (DC) is a rare multisystem bone marrow failure syndrome that displays marked clinical and genetic heterogeneity. X-linked recessive, autosomal dominant and autosomal recessive forms of the disease are recognized. The gene that is mutated in the X-linked form of the disease is DKC1. The DKC1-encoded protein, dyskerin, is a component of small nucleolar ribonucleoprotein particles, which are important in ribosomal RNA processing, and of the telomerase complex. The autosomal dominant form of DC is due to mutations in the gene for the RNA component of telomerase (TERC). Because both dyskerin and TERC are components of the telomerase complex and all patients with DC have short telomeres, the principal pathology of DC appears to relate to telomerase dysfunction, although defects in ribosomal processing via dyskerin's involvement in pseudouridylation cannot be completely ruled out. The gene or genes involved in autosomal recessive DC remain elusive, although genes whose products are required for telomere maintenance remain strong candidates. The study of DC highlights the importance of telomerase in humans and how its deficiency results in multiple abnormalities, including premature aging, bone marrow failure, and cancer.