Shu-Yun Le

A common structural core in the internal ribosome entry sites of picornavirus, hepatitis C virus, and pestivirus

Cap-independent translations of viral RNAs of enteroviruses and rhinoviruses, cardioviruses and aphthoviruses, hepatitis A and C viruses (HAV and HCV), and pestivirus are initiated by the direct binding of 40S ribosomal subunits to acis-acting genetic element termed theinternal ribosome entry site (IRES) orribosome landing pad (RLP) in the 5′ noncoding region (5′NCR). RNA higher ordered structure models for these IRES elements were derived by a combined approach using thermodynamic RNA folding, Monte Carlo simulation, and phylogenetic comparative analysis. The structural differences among the three groups of picornaviruses arise not only from point mutations, but also from the addition or deletion of structural domains. However, a common core can be identified in the proposed structural models of these IRES elements from enteroviruses and rhinoviruses, cardioviruses and aphthoviruses, and HAV. The common structural core identified within the picornavirus IRES is also conserved in the 5′NCR of the divergent viruses, HCV, and pestiviruses. Furthermore, the proposed structural motif shares a structural feature similar to that observed in the catalytic core of the group I intron. The conserved structural motif from these divergent sequences that looks like the common core region of group I introns is probably a crucial element involved in the IRES-dependent translation.

A computational procedure for assessing the significance of RNA secondary structure

In our recent series of papers, we have used the structures of statistical significance from Monte Carlo simulations to improve the predictions of secondary structure of RNA and to analyze the possible role of locally significant structures in the life cycle of human immunodeficiency virus. Because of intensive computational requirements for Monte Carlo simulation, it becomes impractical even using a supercomputer to assess the significance of a structure with a window size greater than 200 along an RNA sequence of 1000 bases or more. In this paper, we have developed a new procedure that drastically reduces the time needed to assess the significance of structures. In fact, the efficiency of this new method allows us to assess structures on the VAX as well as the CRAY.

Conserved tertiary structure elements in the 5′ untranslated region of human enteroviruses and rhinoviruses

Abstract A combination of comparative sequence analysis and thermodynamic methods reveals the conservation of tertiary structure elements in the 5′ untranslated region (UTR) of human enteroviruses and rhinoviruses. The predicted common structural elements occur in the 3′ end of a segment that is critical for internal ribosome binding, termed “ribosome landing pad” (RLP), of polioviruses. Base pairings between highly conserved 17-nucleotide (nt) and 21-nt sequences in the 5′ UTR of human enteroviruses and rhinoviruses constitute a predicted pseudoknot that is significantly more stable than those that can be formed from a large set of randomly shuffled sequences. A conserved single-stranded polypyrimidine tract is located between two conserved tertiary elements. R. Nicholson, J. Pelletier, S.-Y. Le, and N. Sonenberg (1991, J. Virol. 65, 5886–5894) demonstrated that the point mutations of 3-nt UUU out of an essential 4-nt pyrimidine stretch sequence UUUC abolished translation. Structural analysis of the mutant sequence indicates that small point mutations within the short polypyrimidine sequence would destroy the tertiary interaction in the predicted, highly ordered structure. The proposed common tertiary structure can offer experimentalists a model upon which to extend the interpretations for currently available data. Based on these structural features possible base-pairing models between human enteroviruses and 18 S rRNA and between human rhinoviruses and 18 S rRNA are proposed. The proposed common structure implicates a biological function for these sequences in translational initiation.

Discovering well-ordered folding patterns in nucleotide sequences

MOTIVATION Growing evidence demonstrates that local well-ordered structures are closely correlated with cis-acting elements in the post-transcriptional regulation of gene expression. The prediction of a well-ordered folding sequence (WFS) in genomic sequences is very helpful in the determination of local RNA elements with structure-dependent functions in mRNAs. RESULTS In this study, the quality of local WFS is assessed by the energy difference (E(diff)) between the free energies of the global minimal structure folded in the segment and its corresponding optimal restrained structure (ORS). The ORS is an optimal structure under the condition in which none of the base-pairs in the global minimal structure is allowed to form. Those WFSs in HIV-1 mRNA, various ferritin mRNAs and genomic sequences containing let-7 RNA gene were searched by a novel method, ed_scan. Our results indicate that the detected WFSs are coincident with known Rev response element in HIV-1 mRNA, iron-responsive elements in ferritin mRNAs and small let-7 RNAs in Caenorhabditis elegans, Caenorhabditis briggsae and Drosophila melanogaster genomic sequences. Statistical significance of the WFS is addressed by a quantitative measure Zscr(e) that is a z-score of E(diff) and extensive random simulations. We suggest that WFSs with high statistical significance have structural roles involving their sequence information. AVAILABILITY The source code of ed_scan is available via anonymous ftp as ftp://ftp.ncifcrf.gov/pub/users/shuyun/scan/ed_scan.tar.

The human endogenous retrovirus K Rev response element coincides with a predicted RNA folding region.

Human endogenous retrovirus K (HERV-K) is the name given to an approximately 30-million-year-old family of endogenous retroviruses present at >50 copies per haploid human genome. Previously, the HERV-K were shown to encode a nuclear RNA export factor, termed K-Rev, that is the functional equivalent of the H-Rev protein encoded by human immunodeficiency virus type 1. HERV-K was also shown to contain a cis-acting target element, the HERV-K Rev response element (K-RRE), that allowed the nuclear export of linked RNA transcripts in the presence of either K-Rev or H-Rev. Here, we demonstrate that the functionally defined K-RRE coincides with a statistically highly significant unusual RNA folding region and present a potential RNA secondary structure for the approximately 416-nt K-RRE. Both in vitro and in vivo assays of sequence specific RNA binding were used to map two primary binding sites for K-Rev, and one primary binding site for H-Rev, within the K-RRE. Of note, all three binding sites map to discrete predicted RNA stem-loop subdomains within the larger K-RRE structure. Although almost the entire 416-nt K-RRE was required for the activation of nuclear RNA export in cells expressing K-Rev, mutational inactivation of the binding sites for K-Rev resulted in the selective loss of the K-RRE response to K-Rev but not to H-Rev. Together, these data strongly suggest that the K-RRE, like the H-RRE, coincides with an extensive RNA secondary structure and identify specific sites within the K-RRE that can recruit either K-Rev or H-Rev to HERV-K RNA transcripts.

Ion-RNA Interactions in the RNA Pseudoknot of a Ribosomal Frameshifting Site: Molecular Modeling Studies

The three-dimensional (3-D) structure of a RNA pseudoknot that causes the efficient ribosomal frameshifting in the gag-pro region of mouse mammary tumor virus (MMTV) has been determined recently by nuclear magnetic resonance (NMR) studies. But since the structure refinement in the studies did not use metal ions and waters, it is not clear how metal ions participate in the stabilization of the pseudoknot, and what kind of ion-RNA interactions dominate in the tertiary contacts for the RNA pseudoknotting. Based on the reported structure data of the pseudoknot VPK of MMTV, we gradually refined the structure by restrained molecular dynamics (MD) using NMR distance restraints. Restrained MD simulation of the RNA pseudoknot was performed with sodium ions and water molecules. Our results are in good agreement with known NMR data and delineate the importance of the metal ion coordination in the stability of the pseudoknot. In the non-coaxially stacking pseudoknot, stem 1 (S1), stem 2 (S2), and the intervening A14 involves unconventional stacking of base pairs coordinated by Na+ and/or bridging water molecules. A6 and G7 of loop L1 make a perfect base stacking in the major groove and are further stabilized by coordinated Na+ ions and water molecules. The first 4-nucleotide (nt) ACUC of loop L2 form a sharp turn and the following 4-nt AAAA cross the minor groove of S1 and are steadied by interactions with the nucleotides of S , bridging water molecules and coordinated Na+ ions. Our studies suggest that the metal ion plays a crucial role in the RNA pseudoknotting of VPK. In the stacking interior of S1 and S2, the Na+ ion is positioned in the major groove and interacts directly with the carbonyl group O6 of G28 and carbonyl group O4 of U13 in the wobble base pair U13:G28. The ion-RNA interactions in MMTV VPK not only stabilize the RNA pseudoknot but also modify the electrostatic properties of the nucleotides at the critical parts of the pseudoknot VPK.

A data mining approach to discover unusual folding regions in genome sequences

Numerous experiments and analyses of RNA structures have revealed that the local distinct structure closely correlates with the biological function. In this study, we present a data mining approach to discover such unusual folding regions (UFRs) in genome sequences. Our approach is a three-step procedure. During the first step, the quality of a local structure different from a random folding in a genomic sequence is evaluated by two z-scores, significance score (SIGSCR) and stability score (STBSCR) of the local segment. The two scores are computed by sliding a fixed window stepped a base along the sequence from the start to end position. Next, based on the non-central Students t distribution theory we derive a linearly transformed non-central Students t distribution (LTNSTD) to describe the distribution of SIGSCR and STBSCR computed in the sequence. In the third step, we extract these significant UFRs from the sequence whose SIGSCR and/or STBSCR are greater or less than a given threshold calculated from the derived LTNSTD. Our data mining approach is successfully applied to the complete genome of Mycoplasma genitalium (M. gen) and discovers these statistical extremes in the genome. By comparisons with the two scores computed from randomly shuffled sequences of the entire M. gen genome, our results demonstrate that the UFRs in the M. gen sequence are not selected by chance. These UFRs may imply an important structure role involved in their sequence information.

Statistical inference for well-ordered structures in nucleotide sequences

Distinct, local structures are frequently correlated with functional RNA elements involved in post-transcriptional regulation of gene expression. Discovery of microRNAs (miRNAs) suggests that there are a large class of small noncoding RNAs in eukaryotic genomes. These miRNAs have the potential to form distinct fold-back stem-loop structures. The prediction of those well-ordered folding sequences (WFS) in genomic sequences is very helpful for our understanding of RNA-based gene regulation and the determination of local RNA elements with structure-dependent functions. In this study, we describe a novel method for discovering the local WFS in a nucleotide sequence by Monte Carlo simulation and RNA folding. In the approach the quality of a local WFS is assessed by the energy difference (E/sub diff/) between the optimal structure folded in the local segment and its corresponding optimal, restrained structure where all the previous base pairings formed in the optimal structure are prohibited. Distinct WFS can be discovered by scanning successive segments along a sequence for evaluating the difference between Ediff of the natural sequence and those computed from randomly shuffled sequences. Our results indicate that the statistically significant WFS detected in the genomic sequences of Caenorhabditis elegans (C.elegans) F49E12, T07C5, T07D1, T10H9, Y56A3A and Y71G12B are coincident with known fold-back stem-loops found in miRNA precursors. The potential and implications of our method in searching for miRNAs in genomes is discussed.

An improved secondary structure computation method and its application to intervening sequences in the human alpha-like globin mRNA precursors

Current secondary structure prediction computations have a serious drawback. The calculated thermodynamically most stable structure often differs from that observed in solution or in crystal form. In this paper we suggest a way to partially over-come some of these limitations by simulating the RNA folding process and calculating the frequencies of occurrence of the various substructures obtained. The frequently recurring substructures are then selected to construct the secondary structure of the whole RNA. 142 tRNA molecules and an E. coli 16S rRNA molecule have been examined by this method. The percentage of successful prediction of the correct helices are significantly higher than those calculated previously. The secondary structures of intervening sequences (IVSs) excised from human alpha-like globin pre-mRNAs are also computed. Thus, in this method the secondary structures obtained are composed of the statistically more significant substructures. This has also been demonstrated by using randomly shuffled sequences. The secondary structures of each of the randomized sequences are computed and their mean and standard deviations are used in evaluating the significance of the substructures obtained in the folding of the biological sequence. Some potentially appealing structural features aligning adjacent exons for ligation have been found.

DETECTION OF UNUSUAL RNA FOLDING REGIONS IN HIV AND SIV SEQUENCES

We have developed a method for detecting more stable and significant folding regions relative to others in the sequence. The algorithm is based on the calculation of the lowest free energy of RNA secondary structures and Monte Carlo simulation. For any given RNA segment, the stability and statistical significance of RNA folding are assessed by two measures: the stability score and the significance score. The stability score measures the degree of thermodynamic stability of the segment between all possible biological segments in the RNA sequence. The significance score characterizes the specific arrangement of the nucleotides in the segment that could imply a structural role for the sequence information. Using these two measures, we are able to detect a series of distinct folding regions where highly stable and statistically significant secondary structures occur in human immunodeficiency virus (HIV) and simian immunodeficiency virus (SIV) sequences.