Comprehensive survey and geometric classification of base triples in RNA structures

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

Three-dimensional motifs from the SCOR, structural classification of RNA database: extruded strands, base triples, tetraloops and U-turns

Release 2.0.1 of the Structural Classification of RNA (SCOR) database, http://scor.lbl.gov, contains a classification of the internal and hairpin loops in a comprehensive collection of 497 NMR and X-ray RNA structures. This report discusses findings of the classification that have not been reported previously. The SCOR database contains multiple examples of a newly described RNA motif, the extruded helical single strand. Internal loop base triples are classified in SCOR according to their three-dimensional context. These internal loop triples contain several examples of a frequently found motif, the minor groove AGC triple. SCOR also presents the predominant and alternate conformations of hairpin loops, as shown in the most well represented tetraloops, with consensus sequences GNRA, UNCG and ANYA. The ubiquity of the GNRA hairpin turn motif is illustrated by its presence in complex internal loops.     

Conformational polymorphism in G-tetraplex structures: strand reversal by base flipover or sugar flipover.

Guanine rich sequences adopt a variety of four stranded structures, which differ in strand orientation and conformation about the glycosidic bond even though they are all stabilised by Hoogsteen hydrogen bonded guanine tetrads. Detailed model building and molecular mechanics calculations have been carried out to investigate various possible conformations of guanines along a strand and different possible orientations of guanine strands in a G-tetraplex structure. It is found that for an oligo G stretch per se, a parallel four stranded structure with all guanines in anti conformation is favoured over other possible tetraplex structures. Hence an alternating syn-anti arrangement of guanines along a strand is likely to occur only in folded back tetraplex structures with antiparallel G strands. Our study provides a theoretical rationale for the observed alternation of glycosidic conformation and the inverted stacking arrangement arising from base flipover, in antiparallel G-tetraplex structures and also highlights the various structural features arising due to different types of strand orientations. The molecular mechanics calculations help in elucidating the various interactions which stabilize different G-tetraplex structures and indicate that screening of phosphate charge by counterions could have a dramatic effect on groove width in these four stranded structures.     

The relation between hairpin formation by mitochondrial WANCY tRNAs and the occurrence of the light strand replication origin in Lepidosauria

Mitochondrial light strand DNA replication is initiated at light strand replication origins (OLs), short stem-loop hairpins formed by the heavy strand DNA. OL-like secondary structures are also formed by heavy strand DNA templating for the five tRNAs adjacent to OLs, the WANCY tRNA cluster. We tested whether natural OL absence associates with greater capacities for formation of OL-like structures by WANCY tRNA genes. Using lepidosaurian taxa (Sphenodon, lizards and amphisbaenids), we compared WANCY tRNA capacities to form OL-like structures between 248 taxa possessing an OL with 131 taxa without OL (from different families). On average, WANCY tRNA genes form more OL-like structures in the absence of a regular OL than in its presence. Formation of OL-like structures by WANCY tRNAs follows hierarchical patterns that may reduce competition between the tRNA's translational function and its secondary OL function: the rarer the tRNA's cognate amino acid, the greater the capacity to form OL-like structures. High OL-forming capacities for neighboring tRNAs are avoided. Because OL absence usually occurs in taxa with reduced genomes, increased formation of OL-like structures by WANCY tRNAs might result from selection for greater metabolic efficiency. Further analyses suggest that OL loss is one of the latest steps in genome reduction, and promotes the increase in formation of OL-like structures by WANCY tRNA genes in Lepidosauria.     

Homing endonucleases: from microbial genetic invaders to reagents for targeted DNA modification

Homing endonucleases are microbial DNA-cleaving enzymes that mobilize their own reading frames by generating double strand breaks at specific genomic invasion sites. These proteins display an economy of size, and yet recognize long DNA sequences (typically 20 to 30 base pairs). They exhibit a wide range of fidelity at individual nucleotide positions in a manner that is strongly influenced by host constraints on the coding sequence of the targeted gene. The activity of these proteins leads to site-specific recombination events that can result in the insertion, deletion, mutation, or correction of DNA sequences. Over the past fifteen years, the crystal structures of representatives from several homing endonuclease families have been solved, and methods have been described to create variants of these enzymes that cleave novel DNA targets. Engineered homing endonucleases proteins are now being used to generate targeted genomic modifications for a variety of biotech and medical applications.     

Molecular Mechanisms of Radiation Induced Dna Damage: H-Addition to Bases, Direct Ionization and Double Strand Break

Structures and properties of C.(-H) and of TH. were obtained from quantum mechanical calculations. New AMBER parameters for these radicals were obtained to fit their structures and charge distributions. Molecular mechanics simulations of the conformational changes induced in a 12-mer of DNA, d(CGCGAATTCGCG), by these radicals show that the distances between the base and the C₂' of the sugar becomes shorter. Such changes suggest that the base radical can abstract the H,' and transfer the radical from the base to the sugar. Once the radical becomes centered on the sugar a strand break can follow. A simultaneous formation of guanine and thymine radicals on opposite strands may lead to a double strand break.     

Analysis of sequences and predicted structures required for viral satellite RNA accumulation by in vivo genetic selection.

In vivo genetic selection was used to study the sequences and structures required for accumulation of subviral sat-RNA C associated with turnip crinkle virus (TCV). This technique is advantageous over site-specific mutagenesis by allowing side-by-side selection from numerous sequence possibilities as well as sequence evolution. A 22 base hairpin and 6 base single-stranded tail located at the 3'-terminus of sat-RNA C were previously identified as the promoter for minus strand synthesis. Approximately 50% of plants co-inoculated with TCV and sat-RNA C containing randomized sequence in place of the 22 base hairpin accumulated sat-RNA in uninoculated leaves. The 22 base region differed in sat-RNA accumulating in all infected plants, but nearly all were predicted to fold into a hairpin structure that maintained the 6 base tail as a single-stranded sequence. Two additional rounds of sat-RNA amplification led to four sequence family 'winners', with three families containing multiple variants, indicating that evolution of these sequences was occurring in plants. Three of the four sequence family winners had the same 3 bp at the base of the stem as wild-type sat-RNA C. Two of the winners shared 15 of 22 identical bases, including the entire stem region and extending two bases into the loop. These results demonstrate the utility of the in vivo selection approach by showing that both sequence and structure contribute to a more active 3'-end region for accumulation of sat-RNA C.     

Formation of novel hairpin structures by telomeric C-strand oligonucleotides.

Telomeres are specialized structures at the ends of chromosomes that are required for long term chromosome stability and replication of the chromosomal terminus. Telomeric DNA consists of simple repetitive sequences with one strand G-rich relative to the other, C-rich, strand. Evolutionary conservation of this feature of telomeric repeat sequences suggests that they have specific structural characteristics involved in telomere function. Absorbance thermal denaturation, chemical modification and non-denaturing gel electrophoretic analyses showed that telomeric C-strand oligonucleotides form stable non-Watson-Crick hairpin structures containing C.C+ base pairs. Formation of such hairpins may facilitate previously reported G-strand exclusive interactions.     

Molecular Mechanisms of Radiation Induced Dna Damage: H-Addition to Bases,Direct Ionization and Double Strand Break

Structures and properties of C.(-H) and of TH. were obtained from quantum mechanical calculations. New AMBER parameters for these radicals were obtained to fit their structures and charge distributions. Molecular mechanics simulations of the conformational changes induced in a 12-mer of DNA, d(CGCGAATTCGCG), by these radicals show that the distances between the base and the C₂' of the sugar becomes shorter. Such changes suggest that the base radical can abstract the H,' and transfer the radical from the base to the sugar. Once the radical becomes centered on the sugar a strand break can follow. A simultaneous formation of guanine and thymine radicals on opposite strands may lead to a double strand break.     

Patterns,structures, and amino acid frequencies in structural building blocks,a protein secondary structure classification scheme

To study local structures in proteins, we previously developed an autoassociative artificial neural network (autoANN) and clustering tool to discover intrinsic features of macromolecular structures. The hidden unit activations computed by the trained autoANN are a convenient low-dimensional encoding of the local protein backbone structure. Clustering these activation vectors results in a unique classification of protein local structural features called Structural Building Blocks (SBBs). Here we describe application of this method to a larger database of proteins, verification of the applicability of this method to structure classification, and subsequent analysis of amino acid frequencies and several commonly occurring patterns of SBBs. The SBB classification method has several interesting properties: 1) it identifies the regular secondary structures, α helix and β strand; 2) it consistently identifies other local structure features (e.g., helix caps and strand caps); 3) strong amino acid preferences are revealed at some positions in some SBBs; and 4) distinct patterns of SBBs occur in the “random coil” regions of proteins. Analysis of these patterns identifies interesting structural motifs in the protein backbone structure, indicating that SBBs can be used as “building blocks” in the analysis of protein structure. This type of pattern analysis should increase our understanding of the relationship between protein sequence and local structure, especially in the prediction of protein structures. © 1997 Wiley-Liss, Inc.     

Base-flipping dynamics in a DNA hairpin processing reaction

Many enzymes that repair or modify bases in double-stranded DNA gain access to their substrates by base flipping. Although crystal structures provide stunning snap shots, biochemical approaches addressing the dynamics have proven difficult, particularly in complicated multi-step reactions. Here, we use protein–DNA crosslinking and potassium permanganate reactivity to explore the base-flipping step in Tn5 transposition. We present a model to suggest that base flipping is driven by a combination of factors including DNA bending and the intrusion of a probe residue. The forces are postulated to act early in the reaction to create a state of tension, relieved by base flipping after cleavage of the first strand of DNA at the transposon end. Elimination of the probe residue retards the kinetics of nicking and reduces base flipping by 50%. Unexpectedly, the probe residue is even more important during the hairpin resolution step. Overall, base flipping is pivotal to the hairpin processing reaction because it performs two opposite but closely related functions. On one hand it disrupts the double helix, providing the necessary strand separation and steric freedom. While on the other, transposase appears to position the second DNA strand in the active site for cleavage using the flipped base as a handle.     

The Effect of Attachment Site Mutations on Strand Exchange in Bacteriophage λ Site-Specific Recombination

Recombination of phage λ attachment sites occurs by sequential exchange of the DNA strands at two specific locations. The first exchange produces a Holliday structure, and the second resolves it to recombinant products. Heterology for base substitution mutations in the region between the two strand exchange points (the overlap region) reduces recombination; some mutations inhibit the accumulation of Holliday structures, others inhibit their resolution to recombinant products. To see if heterology also alters the location of the strand exchange points, we determined the segregation pattern of three single and one multiple base pair substitution mutations of the overlap region in crosses with wild type sites. The mutations are known to differ in the severity of their recombination defect and in the stage of strand exchange they affect. The three single mutations behaved similarly: each segregated into both products of recombination,, and the two products of a single crossover were frequently nonreciprocal in the overlap region. In contrast, the multiple mutation preferentially segregated into one of the two recombinant products, and the two products of a single crossover appeared to be fully reciprocal. The simplest explanation of the segregation pattern of the single mutations is that strand exchanges occur at the normal locations to produce recombinants with mismatched base pairs that are frequently repaired. The segregation pattern of the multiple mutation is consistent with the view that both strand exchanges usually occur to one side of the mutant site. We suggest that the segregation pattern of a particular mutation is determined by which stage of strand exchange it inhibits and by the severity of the inhibition.     

Formation of an intramolecular triple-stranded DNA structure monitored by fluorescence of 2-aminopurine or 6-methylisoxanthopterin

The parallel (recombination) ‘R-triplex’ can accommodate any nucleotide sequence with the two identical DNA strands in parallel orientation. We have studied oligonucleotides able to fold back into such a recombination-like structure. We show that the fluorescent base analogs 2-aminopurine (2AP) and 6-methylisoxanthopterin (6MI) can be used as structural probes for monitoring the integrity of the triple-stranded conformation and for deriving the thermodynamic characteristics of these structures. A single adenine or guanine base in the third strand of the triplex-forming and the control oligonucleotides, as well as in the double-stranded (ds) and single-stranded (ss) reference molecules, was substituted with 2AP or 6MI. The 2AP*(T·A) and 6MI*(C·G) triplets were monitored by their fluorescence emission and the thermal denaturation curves were analyzed with a quasi-two-state model. The fluorescence of 2AP introduced into an oligonucleotide sequence unable to form a triplex served as a negative control. We observed a remarkable similarity between the thermodynamic parameters derived from melting of the secondary structures monitored through absorption of all bases at 260 nm or from fluorescence of the single base analog. The similarity suggests that fluorescence of the 2AP and 6MI base analogs may be used to monitor the structural disposition of the third strand. We consider the data in the light of alternative ‘branch migration’ and ‘strand exchange’ structures and discuss why these are less likely than the R-type triplex.     

The effect of attachment site mutations on strand exchange in bacteriophage lambda site-specific recombination

Recombination of phage lambda attachment sites occurs by sequential exchange of the DNA strands at two specific locations. The first exchange produces a Holliday structure, and the second resolves it to recombinant products. Heterology for base substitution mutations in the region between the two strand exchange points (the overlap region) reduces recombination; some mutations inhibit the accumulation of Holliday structures, others inhibit their resolution to recombinant products. To see if heterology also alters the location of the strand exchange points, we determined the segregation pattern of three single and one multiple base pair substitution mutations of the overlap region in crosses with wild type sites. The mutations are known to differ in the severity of their recombination defect and in the stage of strand exchange they affect. The three single mutations behaved similarly: each segregated into both products of recombination, and the two products of a single crossover were frequently nonreciprocal in the overlap region. In contrast, the multiple mutation preferentially segregated into one of the two recombinant products, and the two products of a single crossover appeared to be fully reciprocal. The simplest explanation of the segregation pattern of the single mutations is that strand exchanges occur at the normal locations to produce recombinants with mismatched base pairs that are frequently repaired. The segregation pattern of the multiple mutation is consistent with the view that both strand exchanges usually occur to one side of the mutant site. We suggest that the segregation pattern of a particular mutation is determined by which stage of strand exchange it inhibits and by the severity of the inhibition.     

Prediction of the structure of the Y+.R-.R(+)-type DNA triple helix by molecular modelling.

Molecular mechanics has been used to predict the structure of the Y+.R-.R(+)-type DNA triple helix, in which a second polypurine strand binds antiparallel to the homopurine strand of a homopurine/homopyrimidine stretch of duplex DNA. From calculations on the sequence d(C)10.d(G)10.d(G)10, two likely structures emerge. One has the glycosidic torsions of the third strand bases in the anti-conformation and Hoogsteen hydrogen-bonds to the purine strand of the duplex, the other has the third strand purines in the syn orientation and uses a reverse-Hoogsteen hydrogen-bonding pattern. Despite the large structural differences between these two types of triplex, calculations performed in vacuo with a distance-dependent dielectric constant to mimic the shielding effect of solvent show them to be energetically very similar, with the latter (syn) slightly preferred. However, if explicit solvent molecules are included in the calculation, the anti conformation is found to be much preferred. This difference in the results seems to stem from an underestimation of short-range electrostatic interactions in the in vacuo simulations. When TAA or TAT base triples are substituted for the sixth CGG triple in the sequence, it is found that, for the solvated model, the third strand base of the TAA triple prefers the syn orientation while that in the TAT triple retains a preference, though reduced, for the anti conformation.     

Reversible dimer formation and stability of the anti-tumour single-chain Fv antibody MFE-23 by neutron scattering,analytical ultracentrifugation,and NMR and FT-IR spectroscopy

Cre recombinase uses two pairs of sequential cleavage and religation reactions to exchange homologous DNA strands between 34 base-pair (bp) LoxP recognition sequences. In the oligomeric recombination complex, a switch between "cleaving" and "non-cleaving" subunit conformations regulates the number, order, and regio-specificity of the strand exchanges. However, the particular sequence of events has been in question. From analysis of strand composition of the Holliday junction (HJ) intermediate, we determined that Cre initiates recombination of LoxP by cleaving the upper strand on the left arm. Cre preferred to react with the left arm of a LoxP suicide substrate, but at a similar rate to the right arm, indicating that the first strand to be exchanged is selected prior to cleavage. We propose that during complex assembly the cleaving subunit preferentially associates with the LoxP left arm, directing the first strand exchange to that side. In addition, this biased assembly would enforce productive orientation of LoxP sites in the recombination synapses. A novel Cre-HJ complex structure in which LoxP was oriented with the left arm bound by the cleaving Cre subunit suggested a physical basis for the strand exchange order. Lys86 and Lys201 interact with the left arm scissile adenine base differently than in structures that have a scissile guanine. These interactions are associated with positioning the 198-208 loop, a structural component of the conformational switch, in a configuration that is specific to the cleaving conformation. Our results suggest that strand exchange order and site alignment are regulated by an "induced fit" mechanism in which the cleaving conformation is selectively stabilized through protein-DNA interactions with the scissile base on the strand that is cleaved first.     

Protein loops on structurally similar scaffolds: database and conformational analysis

A general problem in comparative modeling and protein design is the conformational evaluation of loops with a certain sequence in specific environmental protein frameworks. Loops of different sequences and structures on similar scaffolds are common in the Protein Data Bank (PDB). In order to explore both structural and sequential diversity of them, a data base of loops connecting similar secondary structure fragments is constructed by searching the data base of families of structurally similar proteins and PDB. A total of 84 loop families having 2-13 residues are found among the well-determined structures of resolution better than 2.5 A. Eight alpha-alpha, 20 alpha-beta, 19 beta-alpha, and 37 beta-beta families are identified. Every family contains more than 5 loop motifs. In each family, no loops share same sequence and all the frameworks are well superimposed. Forty-three new loop classes are distinguished in the data base. The structural variability of loops in homologous proteins are examined and shown in 44 families. Motif families are characterized with geometric parameters and sequence patterns. The conformations of loops in each family are clustered into subfamilies using average linkage cluster analysis method. Information such as geometric properties, sequence profile, sequential and structural variability in loop, structural alignment parameters, sequence similarities, and clustering results are provided. Correlations between the conformation of loops and loop sequence, motif sequence, and global sequence of PDB chain are examined in order to find how loop structures depend on their sequences and how they are affected by the local and global environment. Strong correlations (R > 0.75) are only found in 24 families. The best R value is 0.98. The data base is available through the Internet.