RNA secondary structures in the proximal 3'UTR of Indonesian Dengue 1 virus strains.

The characteristics of DENV-1 viruses, isolated during the 2001-2002 outbreak in Indonesia were studied. The secondary structure of the 3'UTR of different DENV-1 strains derived from Indonesian patients was compared with the 3'UTR of previously described DENV-1 sequences. The complete 3'UTR of DENV-1 was sequenced from 13 patients suffering from the severe form of dengue virus infection (dengue hemorrhagic fever). Prediction of RNA secondary structure of the 3'UTR revealed some previously unidentified conserved structures in the proximal region of the 3'UTR, the role of which in viral replication is still unknown. In addition our data suggest that some structural elements previously described in the distal part of the 3'UTR are partly dependent on the proximal part of the UTR. Our data support the existence of previously unidentified conserved secondary structures in the proximal part of the 3'UTR and their roles need to be further investigated.

PseudoBase++: an extension of PseudoBase for easy searching, formatting and visualization of pseudoknots.

Pseudoknots have been recognized to be an important type of RNA secondary structures responsible for many biological functions. PseudoBase, a widely used database of pseudoknot secondary structures developed at Leiden University, contains over 250 records of pseudoknots obtained in the past 25 years through crystallography, NMR, mutational experiments and sequence comparisons. To promptly address the growing analysis requests of the researchers on RNA structures and bring together information from multiple sources across the Internet to a single platform, we designed and implemented PseudoBase++, an extension of PseudoBase for easy searching, formatting and visualization of pseudoknots. PseudoBase++ (http://pseudobaseplusplus.utep.edu) maps the PseudoBase dataset into a searchable relational database including additional functionalities such as pseudoknot type. PseudoBase++ links each pseudoknot in PseudoBase to the GenBank record of the corresponding nucleotide sequence and allows scientists to automatically visualize RNA secondary structures with PseudoViewer. It also includes the capabilities of fine-grained reference searching and collecting new pseudoknot information.

Selective pressures on RNA hairpins in vivo and in vitro.

Comparison of the most stable potential hairpins in the sequences of natural ribozymes with those in the randomized sequences has revealed that the hairpin loop energies are lower than expected by chance. Although these hairpins are not necessarily parts of functional structures, there is a selective pressure to diminish the destabilizing free energies of the hairpin loops. In contrast, no significant bias is observed in the stacking values of the most stable stems. In the ribozymes isolated in vitro the loops of potential hairpins are closer to random values, which can result in less efficient folding rates. Furthermore, the effects of kinetic traps seem to be more significant in the folding pathways of the in vitro isolates due to a potential to form stable stacks incompatible with the functional folds. Similarly to natural ribozyme sequences, the untranslated regions of viral RNAs also form hairpins with relatively low loop free energies. These evolutionary trends suggest ways for efficient engineering of improved RNA constructs on the basis of analysis of in vitro isolates and approaches for the search of regions coding for functional RNA structures in large genome sequences.

Selection of antisense oligodeoxynucleotides against glutathione S-transferase Mu.

The aim of the present study was to identify functional antisense oligodeoxynucleotides (ODNs) against the rat glutathione S-transferase Mu (GSTM) isoforms, GSTM1 and GSTM2. These antisense ODNs would enable the study of the physiological consequences of GSTM deficiency. Because it has been suggested that the effectiveness of antisense ODNs is dependent on the secondary mRNA structures of their target sites, we made mRNA secondary structure predictions with two software packages, Mfold and STAR. The two programs produced only marginally similar structures, which can probably be attributed to differences in the algorithms used. The effectiveness of a set of 18 antisense ODNs was evaluated with a cell-free transcription/translation assay, and their activity was correlated with the predicted secondary RNA structures. Four phosphodiester ODNs specific for GSTM1, two ODNs specific for GSTM2, and four ODNs targeted at both GSTM isoforms were found to be potent, sequence-specific, and RNase H-dependent inhibitors of protein expression. The IC50 value of the most potent ODN was approximately 100 nM. Antisense ODNs targeted against regions that were predicted by STAR to be predominantly single stranded were more potent than antisense ODNs against double-stranded regions. Such a correlation was not found for the Mfold prediction. Our data suggest that simulation of the local folding of RNA facilitates the discovery of potent antisense sequences. In conclusion, we selected several promising antisense sequences, which, when synthesized as biologically stable oligonucleotides, can be applied for study of the physiological impact of reduced GSTM expression.

PseudoBase: structural information on RNA pseudoknots.

PseudoBase is a database containing structural, functional and sequence data related to RNA pseudo-knots. It can be reached at http://wwwbio.LeidenUniv.nl/ approximately Batenburg/PKB.html. For each pseudoknot, thirteen items are stored, for example the relevant sequence, the stem positions of the pseudoknot, the EMBL accession number of the sequence and the support that can be given regarding the reliability of the pseudo-knot. Since the last publication, information on sizes of the stems and the loops in the pseudoknots has been added. Also added are alternative entries that produce surveys of where the pseudoknots are, sorted according to stem size or loop size.

PseudoBase: a database with RNA pseudoknots.

PseudoBase is a database containing structural, functional and sequence data related to RNA pseudo-knots. It can be reached at http://wwwbio. Leiden Univ.nl/ approximately Batenburg/PKB.html. This page will direct the user to a retrieval page from where a particular pseudoknot can be chosen, or to a submission page which enables the user to add pseudoknot information to the database or to an informative page that elaborates on the various aspects of the database. For each pseudoknot, 12 items are stored, e.g. the nucleotides of the region that contains the pseudoknot, the stem positions of the pseudoknot, the EMBL accession number of the sequence that contains this pseudoknot and the support that can be given regarding the reliability of the pseudoknot. Access is via a small number of steps, using 16 different categories. The development process was done by applying the evolutionary methodology for software development rather than by applying the methodology of the classical waterfall model or the more modern spiral model.

An approximation of loop free energy values of RNA H-pseudoknots.

A set of free energy values is suggested for RNA H-pseudoknot loops. The parameters are adjusted to be consistent with the theory of polymer thermodynamics and known data on pseudoknots. The values can be used for estimates of pseudoknot stabilities and computer predictions of RNA structures.

Dynamic competition between alternative structures in viroid RNAs simulated by an RNA folding algorithm.

The folding pathways of viroid RNAs were studied using computer simulations by the genetic algorithm for RNA folding. The folding simulations were performed for PSTVd RNAs of both polarities, using the wild-type sequence and some previously known mutants with suggested changes in the stable or metastable structures. It is shown that metastable multihairpin foldings in the minus strand replicative intermediates are established due to the specific folding pathway that ensures the absence of the most stable rod-like structure. Simulations of the PSTVd minus strand folding during transcription reveal a metastable hairpin, formed in the left terminal domain region of the PSTVd. Despite high sequence variability, this hairpin is conserved in all known large viroids of both subgroups of PSTVd type, and is presumably necessary to guide the folding of the HPII hairpin which is functional in the minus strand. The folding simulations are able to demonstrate the changes in the balance between metastable and stable structures in mutant PSTVd RNAs. The stable rod-like structure of the circular viroid (+) RNA is also folded via a dynamic folding pathway. Furthermore, the simulations show that intermediate steps in the forced evolution of a shortened PSTVd replicon may be reconstructed by a mechanistic model of different folding pathway requirements in plus- and minus-strand RNAs. Thus the formation of viroid RNA structure strongly depends on dynamics of competition between alternative RNA structures. This also suggests that the replication efficiency of viroid sequences may be estimated by a simulation of the folding process.

The influence of a metastable structure in plasmid primer RNA on antisense RNA binding kinetics.

Replication of the ColE1 group plasmids is kinetically regulated by the interaction between plasmid-encoded primer RNA II and antisense RNA I. The binding is dependent on alternative RNA II conformations, formed during the transcription, and effectively inhibits the primer function within some time interval. In this paper, the folding pathways for the wild type and copy number mutants of ColE1 RNA II are studied using simulations by a genetic algorithm. The simulated pathways reveal a transient formation of a metastable structure, which is stabilized by copy number mutations. The folding kinetics of the proposed conformational transitions is calculated using a model of a multistep refolding process with elementary steps of double-helical stem formation or disruption. The approximation shows that the lifetime of the metastable structure is relatively long and is considerably increased in the mutants, resulting in a delay of the formation of the stable RNA II structure, which is the most sensitive to the inhibition by the antisense RNA I. Thus the effect of copy number mutations can be interpreted as a compression of the time window of effective inhibition due to an increased time spent by the RNA II in the metastable state. The implications of metastable foldings in RNA functioning are discussed.

The computer simulation of RNA folding pathways using a genetic algorithm.

A procedure for simulating the RNA folding process using the principles of genetic algorithm is proposed. The method allows one to simulate a folding pathway of RNA, including such processes as disruption of temporarily formed structures, the folding of a molecule during its synthesis and pseudoknot formation. The simulations are able to predict functional metastable foldings and kinetically driven transitions to more stable structures. The analysis of free energies for intermediate foldings allows estimation of the ranges of kinetic refolding barriers and suggests that in some RNAs the selective evolutionary pressure suppresses the possibilities for alternative structures that could form in the course of transcription. It is shown that the folding pathway simulation can result in structure predictions that are more consistent with phylogenetically proven structures than minimum energy solutions. This suggest that RNA folding kinetics is very important for the formation of functional RNA structures. Therefore, apart form its value for predictions of RNA structures, the proposed computer simulations can be a powerful tool in the studies of RNA folding features.

An APL-programmed genetic algorithm for the prediction of RNA secondary structure.

The possibilities of using a genetic algorithm for the prediction of RNA secondary structure were investigated. The algorithm, using the procedure of stepwise selection of the most fit structures (similarly to natural evolution), allows different models of fitness or driving forces determining RNA structure to be easily introduced. This can be used for simulation of the RNA folding process and for the investigation of possible folding pathways. Such an algorithm needs several modifications before it can predict RNA secondary structures. After modification, a fair number of correct stems are predicted, even when using computationally quick, but very crude, fitness criteria such as stem length and stacking energy, including elements of tertiary structure (pseudoknots). The fact that genetic algorithm simulation includes both stem formations and stem disruption allows one to observe intermediate structures that may be used in combination with phylogenetic or experimental research.

The function of host discrimination and superparasitization in parasitoids.

Host discrimination, i.e. the ability to distinguish unparasitized hosts from parasitized ones, and to reject the latter for egg laying is present in many parasitic wasp species. This property is classically considered as an example of contest competition, and is supposed to have a number of functions. However, different species do not react to each other's marks and lay eggs in hosts parasitized by the other species. Apparently the marks used for recognition are specific.Multiparasitization is the best strategy when hosts are scarce and the egg supplies of the parasitoids are not limited. Interspecific host discrimination is not an ESS.Superparasitization within one species would have selective advantage if the number of unparasitized hosts is small and the wasp has a reasonable chance to lay her egg in a host that is not parasitized by herself, and if the chance for her offspring to survive the competitive battle with the first parasitoid larva is not too small. This is shown to be the case.However, marks are not individual and wasps cannot distinguish hosts parasitized by themselves from those parasitized by others. The hypothesis is tested that the egg laying strategy (i.e. the decision to superparasitize) of wasps is dependent on the number of conspecifics that is searching simultaneously for hosts, since this determines the chance that a parasitized host encountered by a wasp is parasitized by herself.It is shown that host discrimination cannot be regarded as a case of contest competition. Other aspects of superparasitization, related to interference and population regulation, sex allocation and encapsulation are briefly discussed.