Modelling the three-dimensional structure of ribonucleic acids

An intricate RNA structure with two tRNA-derived motifs directs complex formation between yeast aspartyl-tRNA synthetase and its mRNA

Accurate translation of genetic information necessitates the tuned expression of a large group of genes. Amongst them, controlled expression of the enzymes catalyzing the aminoacylation of tRNAs, the aminoacyl-tRNA synthetases, is essential to insure translational fidelity. In the yeast Saccharomyces cerevisiae, expression of aspartyl-tRNA synthetase (AspRS) is regulated in a process necessitating recognition of the 5' extremity of AspRS messenger RNA (mRNA(AspRS)) by its translation product and adaptation to the cellular tRNA(Asp) concentration. Here, we have established the folding of the approximately 300 nucleotides long 5' end of mRNA(AspRS) and identified the structural signals involved in the regulation process. We show that the regulatory region in mRNA(AspRS) folds in two independent and symmetrically structured domains spaced by two single-stranded connectors. Domain I displays a tRNA(Asp) anticodon-like stem-loop structure with mimics of the aspartate identity determinants, that is restricted in domain II to a short double-stranded helix. The overall mRNA structure, based on enzymatic and chemical probing, supports a three-dimensional model where each monomer of yeast AspRS binds one individual domain and recognizes the mRNA structure as it recognizes its cognate tRNA(Asp). Sequence comparison of yeast genomes shows that the features within the mRNA recognized by AspRS are conserved in different Saccharomyces species. In the recognition process, the N-terminal extension of each AspRS subunit plays a crucial role in anchoring the tRNA-like motifs of the mRNA on the synthetase.

Untying the FIV frameshifting pseudoknot structure by MS3D

The structure of the putative feline immunodeficiency virus (FIV) ribosomal frameshifting pseudoknot (PK) has been investigated by a mass spectrometric three-dimensional (MS3D) approach, which involves the application of established solvent-accessibility probes and chemical crosslinkers with detection by electrospray ionization (ESI) Fourier transform mass spectrometry (FTMS). Regardless of their size, probed substrates can be treated with ribonucleases and analyzed by ESI-FTMS to obtain the correct position of chemically modified nucleotides. Protection maps and distance information can be utilized to generate 3D models using the constraint satisfaction algorithm provided by MC-SYM and the energy minimization modules included in CNS. Control experiments were performed on a mutant of mouse mammary tumor virus pseudoknot (VPK), for which an NMR structure is available. Comparison between the MS3D model and the high-resolution structure provided a approximately 3A root-mean-square deviation calculated from all the atoms present in double-stranded regions. Applied to FIV-PK, the MS3D approach confirmed that the selected sequence could fold into an actual pseudoknot, supporting the sequence alignment predictions. Characteristic features of H-type pseudoknots were recognized immediately, but a putative A13-U30 pair was not observed at the stem junction, making FIV-PK resemble VPK more closely than the initially suggested simian retrovirus type-1 pseudoknot. In our model, the unpaired U30 protrudes into the medium, while the hinging A13 assumes a stacked conformation that enables the stems to form a approximately 60 degrees bend and relieve the strain caused by a short loop 1. The model provided the basis to explain the different alkylation patterns observed in the absence and presence of Mg(2+), suggesting the possible formation of a specific metal-binding site between loop 1 and stem 2. This instance illustrates how the MS3D model of FIV-PK can be utilized effectively to generate hypotheses and support functional observations in the absence of a high-resolution structure.

Architecture and folding mechanism of the Azoarcus Group I Pre-tRNA

Self-splicing RNAs must evolve to function in their specific exon context. The conformation of a group I pre-tRNA(ile) from the bacterium Azoarcus was probed by ribonuclease T(1) and hydroxyl radical cleavage, and by native gel electrophoresis. Biochemical data and three-dimensional models of the pre-tRNA showed that the tRNA is folded, and that the tRNA and intron sequences form separate tertiary domains. Models of the active site before steps 1 and 2 of the splicing reaction predict that exchange of the external G-cofactor and the 3'-terminal G is accomplished by a slight conformational change in P9.0 of the Azoarcus group I intron. Kinetic assays showed that the pre-tRNA folds in minutes, much more slowly than the intron alone. The dependence of the folding kinetics on Mg(2+) and the concentration of urea, and RNase T(1) experiments showed that formation of native pre-tRNA is delayed by misfolding of P3-P9, including mispairing between residues in P9 and the tRNA. Thus, although the intron and tRNA sequences form separate domains in the native pre-tRNA, their folding is coupled via metastable non-native base-pairs. This could help prevent premature processing of the 5' and 3' ends of unspliced pre-tRNA.

The three-dimensional architecture of the class I ligase ribozyme

The class I ligase ribozyme catalyzes a Mg(++)-dependent RNA-ligation reaction that is chemically analogous to a single step of RNA polymerization. Indeed, this ribozyme constitutes the catalytic domain of an accurate and general RNA polymerase ribozyme. The ligation reaction is also very rapid in both single- and multiple-turnover contexts and thus is informative for the study of RNA catalysis as well as RNA self-replication. Here we report the initial characterization of the three-dimensional architecture of the ligase. When the ligase folds, several segments become protected from hydroxyl-radical cleavage, indicating that the RNA adopts a compact tertiary structure. Ribozyme folding was largely, though not completely, Mg(++) dependent, with a K(1/2[Mg]) < 1 mM, and was observed over a broad temperature range (20 degrees C -50 degrees C). The hydroxyl-radical mapping, together with comparative sequence analyses and analogy to a region within 23S ribosomal RNA, were used to generate a three-dimensional model of the ribozyme. The predictive value of the model was tested and supported by a photo-cross-linking experiment.

Assembly of core helices and rapid tertiary folding of a small bacterial group I ribozyme

Compact but non-native intermediates have been implicated in the hierarchical folding of several large RNAs, but there is little information on their structure. In this article, ribonuclease and hydroxyl radical cleavage protection assays showed that base pairing of core helices stabilize a compact state of a small group I ribozyme from Azoarcus pre-tRNA(ile). Base pairing of the ribozyme core requires 10-fold less Mg(2+) than stable tertiary interactions, indicating that assembly of helices in the catalytic core represents a distinct phase that precedes the formation of native tertiary structure. Tertiary folding occurs in <100 ms at 37 degrees C. Such rapid folding is unprecedented among group I ribozymes and illustrates the association between structural complexity and folding time. A 3D model of the Azoarcus ribozyme was constructed by identifying homologous sequence motifs in rRNA. The model reveals distinct structural features, such as a large interface between the P4-P6 and P3-P9 domains, that may explain the unusual stability of the Azoarcus ribozyme and the cooperativity of folding.

The modular structure of Escherichia coli threonyl-tRNA synthetase as both an enzyme and a regulator of gene expression

In addition to its role in tRNA aminoacylation, Escherichia coli threonyl-tRNA synthetase is a regulatory protein which binds a site, called the operator, located in the leader of its own mRNA and inhibits translational initiation by competing with ribosome binding. This work shows that the two essential steps of regulation, operator recognition and inhibition of ribosome binding, are performed by different domains of the protein. The catalytic and the C-terminal domain of the protein are involved in binding the two anticodon arm-like structures in the operator whereas the N-terminal domain of the enzyme is responsible for the competition with the ribosome. This is the first demonstration of a modular structure for a translational repressor and is reminiscent of that of transcriptional regulators. The mimicry between the operator and tRNA, suspected on the basis of previous experiments, is further supported by the fact that identical regions of the synthetase recognize both the operator and the tRNA anticodon arm. Based on these results, and recent structural data, we have constructed a computer-derived molecular model for the operator-threonyl-tRNA synthetase complex, which sheds light on several essential aspects of the regulatory mechanism.

Molecular modeling of the three-dimensional structure of the bacterial RNase P holoenzyme

Bacterial ribonuclease P (RNase P), an enzyme involved in tRNA maturation, consists of a catalytic RNA subunit and a protein cofactor. Comparative phylogenetic analysis and molecular modeling have been employed to derive secondary and tertiary structure models of the RNA subunits from Escherichia coli (type A) and Bacillus subtilis (type B) RNase P. The tertiary structure of the protein subunit of B.subtilis and Staphylococcus aureus RNase P has recently been determined. However, an understanding of the structure of the RNase P holoenzyme (i.e. the ribonucleoprotein complex) is lacking. We have now used an EDTA-Fe-based footprinting approach to generate information about RNA-protein contact sites in E.coli RNase P. The footprinting data, together with results from other biochemical and biophysical studies, have furnished distance constraints, which in turn have enabled us to build three-dimensional models of both type A and B versions of the bacterial RNase P holoenzyme in the absence and presence of its precursor tRNA substrate. These models are consistent with results from previous studies and provide both structural and mechanistic insights into the functioning of this unique catalytic RNP complex.

Monitoring intermediate folding states of the td group I intron in vivo

Group I introns consist of two major structural domains, the P4-P6 and P3-P9 domains, which assemble through interactions with peripheral extensions to fold into an active ribozyme. To assess group I intron folding in vivo, we probed the structure of td wild-type and mutant introns using dimethyl sulfate. The results suggest that the majority of the intron population is in the native state in accordance with the current structural model, which was refined to include two novel tertiary contacts. The importance of the loop E motif in the P7.1-P7.2 extension in assisting ribozyme folding was deduced from modeling and mutational analyses. Destabilization of stem P6 results in a deficiency in tertiary structure formation in both major domains, while weakening of stem P7 only interferes with folding of the P3-P9 domain. The different impact of mutations on the tertiary structure suggests that they interfere with folding at different stages. These results provide a first insight into the structure of folding intermediates and suggest a putative order of events in a hierarchical folding pathway in vivo.

Protein-dependent transition states for ribonucleoprotein assembly

Native folding and splicing by the Saccharomyces cerevisiae mitochondrial bI5 group I intron RNA is facilitated by both the S. cerevisiae CBP2 and Neurospora crassa CYT-18 protein cofactors. Both protein-bI5 RNA complexes splice at similar rates, suggesting that the RNA active site structure is similar in both ribonucleoproteins. In contrast, the two proteins assemble with the bI5 RNA by distinct mechanisms and bind opposing, but partially overlapping, sides of the group I intron catalytic core. Assembly with CBP2 is limited by a slow, unimolecular RNA folding step characterized by a negligible activation enthalpy. We show that assembly with CYT-18 shows four distinctive features. (1) CYT-18 binds stably to the bI5 RNA at the diffusion controlled limit, but assembly to a catalytically active RNA structure is still limited by RNA folding, as visualized directly using time-resolved footprinting. (2) This mechanism of rapid stable protein binding followed by subsequent assembly steps has a distinctive kinetic signature: the apparent ratio of k(off) to k(on), determined in a partitioning experiment, differs from the equilibrium K(d) by a large factor. (3) Assembly with CYT-18 is characterized by a large activation enthalpy, consistent with a rate limiting conformational rearrangement. (4) Because assembly from the kinetically trapped state is faster at elevated temperature, we can identify conditions where CYT-18 accelerates (catalyzes) bI5 RNA folding relative to assembly with CBP2.

Four-way junctions in antisense RNA-mRNA complexes involved in plasmid replication control: a common theme?

In several groups of bacterial plasmids, antisense RNAs regulate copy number through inhibition of replication initiator protein synthesis. In plasmid R1, we have recently shown that the inhibitory complex between the antisense RNA (CopA) and its target mRNA (CopT) is characterized by the formation of two intermolecular helices, resulting in a four-way junction structure and a side-by-side helical alignment. Based on lead-induced cleavage and ribonuclease (RNase) V(1) probing combined with molecular modeling, a strikingly similar topology is supported for the complex formed between the antisense RNA (Inc) and mRNA (RepZ) of plasmid Col1b-P9. In particular, the position of the four-way junction and the location of divalent ion-binding site(s) indicate that the structural features of these two complexes are essentially the same in spite of sequence differences. Comparisons of several target and antisense RNAs in other plasmids further indicate that similar binding pathways are used to form the inhibitory antisense-target RNA complexes. Thus, in all these systems, the structural features of both antisense and target RNAs determine the topologically possible and kinetically favored pathway that is essential for efficient in vivo control.

MANIP: an interactive tool for modelling RNA

Large RNA structures can be viewed as assemblies of smaller units or modules that are usually clearly identified (helices, hairpin loops, other recurrent motifs, etc.). We have developed a program, MANIP, which allows the rapid assembly of separate motifs (each with a specified sequence) into a complex three-dimensional architecture. The already determined modules are present in a database from which they can be extracted with the appropriate sequence. Their assembly is performed in real time on the computer screen with buttons and dials that command rotation and translation of any chosen fragment with respect to the chosen pivot, or that generate all possible variations of any torsion angle within a specified segment either in the 5' or in the 3' direction. The possible in-built manipulations follow the general stereochemical rules of RNA structure. MANIP automatically recognizes and displays the allowed and nonallowed hydrogen bonds between the residues. The program is interfaced with a rapid and automatic online refinement tool of partial or full assemblies, NUCLIN-NUCLSQ. The refinement protocol incorporates canonical as well as noncanonical base pairing constraints together with restraints imposed by covalent geometry, stereochemistry, and van der Waals contacts. The computer package runs on UNIX Silicon Graphics workstations and is written in C with OpenGL and X11/Motif libraries.

Derivation of the three-dimensional architecture of bacterial ribonuclease P RNAs from comparative sequence analysis

The secondary structure of bacterial RNase P RNA, a ribozyme responsible for the maturation of the 5' end of tRNAs, is well established on the basis of sequence comparison analysis. RNase P RNA secondary structures fall into two types, A and B, which share a common core formed by the assembly of two main folding domains, but differ in their peripheral elements.A revised alignment of 137 available sequences reveals new covariations allowing for the refinement of both types of secondary structures. Phylogenetic evidence is thus provided for the extension of stems P11, P14, P19, P10.1 and P15.1 through further canonical base-pairs or GAellipsisGA mismatches. These refinements led in turn to a new organization of the catalytic core, with coaxial stackings of helices P2 and P19 as well as P1 and P4. New inter-domain tertiary interactions involve loop L9 and helix P1 and loop L8 with helix P4. These features were incorporated into atomic-scale 3D models of RNase P RNA for representatives of each structural type, namely Escherichia coli and Bacillus subtilis. In each model, the juxtaposition of the core helices creates a cradle onto which the pre-tRNA substrate binds with most evolutionarily conserved residues converging towards the cleavage site. The inner cores of both types are stabilized similarly, albeit by different peripheral elements, emphasizing the modular and hierarchical organisation of the architecture of RNase P RNAs. Similarities are thus apparent between the type A modules, P16/P17/P6 and P13/P14, and their type B analogs, P5.1/P15.1 and P10. 1/P10.1a, respectively. Other noteworthy features of these models include compactness and good agreement with published crosslinking data.

Modified nucleotides of tRNAPro restrict interactions in the binary primer/template complex of M-MuLV

In all retroviruses, reverse transcription is primed by a cellular tRNA, which is base-paired through its 3'-terminal 18 nucleotides to a complementary sequence on the viral RNA genome termed the primer binding site (PBS). Evidence for specific primer-template interactions in addition to this standard interaction has recently been demonstrated for several retroviruses. Here, we used chemical and enzymatic probing to investigate the interactions between Moloney murine leukemia virus (M-MuLV) RNA and its natural primer tRNAPro. The existence of extended interactions was further tested by comparing the viral RNA/tRNAPro complex with simplified complexes in which viral RNA or tRNA were reduced to the 18 nt of the PBS or to the complementary tRNA sequence. These data, combined with computer modeling provide important clues on the secondary structure and three-dimensional folding of the M-MuLV RNA/tRNAPro complex. In contrast with other retroviruses, we found that the interaction between tRNAPro and the M-MuLV RNA template is restricted to the standard PBS interaction. In this binary complex, the viral RNA is highly constrained and the rest of tRNAPro is rearranged, with the exception of the anticodon arm, leading to a very compact structure. Unexpectedly, when a synthetic tRNAPro lacking the post-transcriptional modifications is substituted for the natural tRNAPro primer, the interactions between the primer and the viral RNA are extended. Hence, our data suggest that the post-transcriptional modifications of natural tRNAPro prevent additional contacts between tRNAPro and the U5 region of M-MuLV RNA.

Inter-domain cross-linking and molecular modelling of the hairpin ribozyme

The hairpin ribozyme is a small catalytic RNA composed of two helical domains containing a small and a large internal loop and, thus, constitutes a valuable paradigm for the study of RNA structure and catalysis. We have carried out molecular modelling of the hairpin ribozyme to learn how the two domains (A and B) might fold and approach each other. To help distinguish alternative inter-domain orientations, we have chemically synthesized hairpin ribozymes containing 2'-2' disulphide linkages of known spacing (12 or 16 A) between defined ribose residues in the internal loop regions of each domain. The abilities of cross-linked ribozymes to carry out RNA cleavage under single turnover conditions were compared to the corresponding disulphide-reduced, untethered ribozymes. Ribozymes were classed in three categories according to whether their cleavage rates were marginally, moderately, or strongly affected by cross-linking. This rank order of activity guided the docking of the two domains in the molecular modelling process. The proposed three-dimensional model of the hairpin ribozyme incorporates three different crystallographically determined structural motifs: in domain A, the 5'-GAR-3'-motif of the hammerhead ribozyme, in domain B, the J4/5 motif of group I ribozymes, and connecting the two domains, a "ribose zipper", another group I ribozyme feature, formed between the hydroxyl groups of residues A10, G11 of domain A and C25, A24 of domain B. This latter feature might be key to the selection and precise orientation of the inter-domain docking necessary for the specific phosphodiester cleavage. The model provides an important basis for further studies of hairpin ribozyme structure and function.

Non-canonical interactions in a kissing loop complex: the dimerization initiation site of HIV-1 genomic RNA

Retroviruses encapsidate two molecules of genomic RNA that are noncovalently linked close to their 5' ends in a region called the dimer linkage structure (DLS). The dimerization initiation site (DIS) of human immunodeficiency virus type 1 (HIV-1) constitutes the essential part of the DLS in vitro and is crucial for efficient HIV-1 replication in cell culture. We previously identified the DIS as a hairpin structure, located upstream of the major splice donor site, that contains in the loop a six-nucleotide self-complementary sequence preceded and followed by two and one purines, respectively. Two RNA monomers form a kissing loop complex via intermolecular interactions of the six nucleotide self-complementary sequence. Here, we introduced compensatory mutations in the self-complementary sequence and/or a mutation in the flanking purines. We determined the kinetics of dimerization, the thermal stabilities and the apparent equilibrium dissociation constants of wild-type and mutant dimers and used chemical probing to obtain structural information. Our results demonstrate the importance of the 5'-flanking purine and of the two central bases of the self-complementary sequence in the dimerization process. The experimental data are rationalized by triple interactions between these residues in the deep groove of the kissing helix and are incorporated into a three-dimensional model of the kissing loop dimer. In addition, chemical probing and molecular modeling favor the existence of a non-canonical interaction between the conserved adenine residues at the first and last positions in the DIS loop. Furthermore, we show that destabilization of the kissing loop complex at the DIS can be compensated by interactions involving sequences located downstream of the splice donor site of the HIV-1 genomic RNA.

A three-dimensional working model for a guide RNA from Trypanosoma brucei

RNA editing in protozoan parasites is a mitochondrial RNA processing reaction in which exclusively uridylate residues are inserted into, and less frequently deleted from, pre-mRNAs. Molecules central to the process are so-called guide RNAs (gRNAs) which function as templates in the reaction. For a detailed molecular understanding of the mechanism of the editing process knowledge of structural features of gRNAs will be essential. Here we report on a computer-assisted molecular modelling approach to construct the first three-dimensional gRNA model for gND7-506, a ND7-specific gRNA from Trypanosoma brucei. The modelling process relied on chemical modification and enzymatic probing data and was validated by in vitro mutagenesis experiments. The model predicts a reasonably compact structure, where two stem/loop secondary structure elements are brought into close proximity by a triple A tertiary interaction, forming a core element within the centre of the molecule. The model further suggests that the surface of the gRNA is primarily made up of the sugar-phoshate backbone. On the basis of the model, footprinting experiments of gND7-506 in a complex with the gRNA binding protein gBP21 could successfully be interpreted and provide a first picture for the assembly of gRNAs within a ribonucleoprotein complex.

Context dependent RNA-RNA recognition in a three-dimensional model of the 16S rRNA core

A 3-D model of the core of the 16S rRNA of Escherichia coli containing 328 residues has been built in the protein map derived from neutron scattering data with the help of all the available phylogenetic, biochemical, and cross-linking data. The three pseudoknots of the 16S-core cluster, through the arrangement of complex three-, four- and five-way junctions, around the neck and at the subunit interface. The roles in assembly, initiation or elongation of the three pseudoknots in ribosomal dynamics are emphasized. The 530-loop, localized on the periphery of the 30S particle, could be built with and without a pseudoknot independently of the state of the particle. The pseudoknot of the central domain controls the dynamics of an helix connected to the subunit interface which could trigger some mechanism during translation. The process of the model construction is compatible with a folding scenario in which the 5'-terminal pseudoknot controls the assembly of the central junction and the subsequent folding of the 3'-major domain. The modelling, together with the phylogenetic analysis and the experimental data, point to several potential RNA-RNA contacts which depend on the structural and sequence context in which they occur.

New loop-loop tertiary interactions in self-splicing introns of subgroup IC and ID: a complete 3D model of the Tetrahymena thermophila ribozyme

BACKGROUND

Group I introns self-splice via two consecutive trans-esterification reactions in the presence of guanosine cofactor and magnesium ions. Comparative sequence analysis has established that a catalytic core of about 120 nucleotides is conserved in all known group I introns. This core is generally not sufficient for activity, however, and most self-splicing group I introns require non-conserved peripheral elements to stabilize the complete three-dimensional (3D) structure. The physico-chemical properties of group I introns make them excellent systems for unraveling the structural basis of the RNA-RNA interactions responsible for promoting the self-assembly of complex RNAs.

RESULTS

We present phylogenetic and experimental evidence for the existence of three additional tertiary base pairings between hairpin loops within peripheral components of subgroup IC1 and ID introns. Each of these new long range interactions, called P13, P14 and P16, involves a terminal loop located in domain 2. Although domains 2 of IC and ID introns share very strong sequence similarity, their terminal loops interact with domains 5 and 9 (subgroup IC1) and domain 6 (subgroup ID). Based on these tertiary contacts, comparative sequence analysis, and published experimental results such as Fe(II)-EDTA protection patterns, we propose 3D models for two entire group I introns, the subgroup IC1 intron in the large ribosomal precursor RNA of Tetrahymena thermophila and the SdCob.1 subgroup ID intron found in the cytochrome b gene of Saccharomyces douglasii.

CONCLUSIONS

Three-dimensional models of group I introns belonging to four different subgroups are now available. They all emphasize the modular and hierarchical organization of the architecture of group I introns and the widespread use of base-pairings between terminal hairpin loops for stabilizing the folded and active structures of large and complex RNA molecules.

A geometrical approach in folding a pseudoknot motif within the E. coli 16S-RNA

Except for tRNA, the tertiary structure of RNA molecules are very little known. The many possibilities in the arrangement of different helices in space and the flexibility in the single-stranded loops that connect the helical regions make the modeling of the tertiary structure of RNA molecule a very complex task. Here, we introduce an approach to fold RNA tertiary structure based only on the information of the secondary structure and the stereochemistry of the molecule. This approach was used to construct an atomic structure of a pseudoknot (bases 500-545) in the E. coli 16S RNA. The resulting structure is a closely packed molecule that is consistent with the predicted secondary structure and stereochemically feasible. This new approach is very general and easily adaptable. Experimental data (e.g., NMR, fluorescence energy transfer, etc.), as they become available, can be incorporated directly into the approach to improve the accuracy of the modeled structure.

Determining the conformation of RNAs in solution. Application to a retroviral system: structure of the HIV-1 primer binding site region and effect of tRNA(3Lys) binding

RNAs play a crucial and central role in a large variety of biological functions obviously linked to the wide variety of structures that they can adopt. Understanding the function of RNAs thus requires the knowledge of their two- and three-dimensional structures. We describe in detail the way to access the secondary structure of RNAs, by combining sequence comparison, secondary structure prediction by computer and, mainly, experimental data obtained by probing with chemicals and ribonucleases. These approaches were used to investigate secondary structure of the region containing the primer binding site of HIV-1 genomic RNA either free or involved in the binary complex with the replication primer tRNA(3Lys).

Calculations of nucleic acid conformations

The present computational power and sophistication of theoretical approaches to nucleic acid structural investigation are sufficient for the realization of static and dynamic models that correlate accurately with current crystallographic, NMR and solution-probing structural data, and consequently are able to provide valuable insights and predictions for a variety of nucleic acid conformational families. In molecular dynamics simulations, the year 1995 was marked by the foray of fast Ewald methods, an accomplishment resulting from several years' work in the search for an adequate treatment of the electrostatic long-range forces so primordial in nucleic acid behavior. In very large systems, and particularly in the RNA-folding field, techniques originating from artificial intelligence research, like constraint satisfaction programming or genetic algorithms, have established their utility and potential.

DRAWNA: a program for drawing schematic views of nucleic acids

A program for drawing automatically exact and schematic views of nucleic acids is described. The program is written in C ANSI and uses the Silicon Graphics GL and Xirisw libraries within the X11/Motif environment. Through menus, the user can choose, specify, and manipulate in real time the three-dimensional views to be displayed. Drawing options include partitioning of structures into differently colored or shaped fragments, representation of backbones as flat or with conic-section ribbons, display of paired or free bases as rods, and display of surfaces as filled or outlined and stereo or depth-cued views.

A three-dimensional model of hepatitis delta virus ribozyme based on biochemical and mutational analyses

BACKGROUND

Hepatitis delta virus (HDV), which has a single-stranded RNA genome about 1700 nucleotides long, is a satellite virus of hepatitis B, and is associated with a high incidence of fulminant hepatitis and death in infected humans. Like certain pathogenic subviral RNAs that infect plants, HDV RNA features a closed-circular conformation, a rolling-circle mechanism of replication and RNA-catalyzed self-cleaving reactions of both genomic and anti-genomic strands in vitro. The catalytic domains cannot be folded into either the hammerhead or hairpin secondary-structure motifs that have been found in other self-cleaving RNAs.

RESULTS

A pseudoknot secondary-structure model has been suggested for the catalytic domain (ribozyme) of HDV RNA. We conducted extensive mutational analyses of regions of the HDV ribozyme predicted in this model to be single stranded, and found that several of them are important for catalytic activity. We used these data, sequence comparisons between different isolates and previously published structural analyses to produce a computer graphic model of the three-dimensional architecture of the HDV ribozyme.

CONCLUSIONS

Our model supports the pseudoknotted structure and rationalizes several observations relating to the lengths of the various stems and the sequence requirements of the single-stranded regions. It also provides insight into the catalytic mechanism of the HDV ribozyme. We specifically propose that residues C75, U20 and C21 form the basis of the catalytic region and are close to the cleavable phosphate.