41 The essential adenosine stacking in a two-base-pair minimal kissing complex

A stable RNA helix requires at least three base pairs. Surprisingly, a tertiary kissing complex formed between two GACG hairpin loops contains only two GC pairs. In the NMR structure of this complex, the two flanking adenosines stack on the kissing GC pair. This observation raised a possibility that the 5’-dangling adenines contribute to the formation and stability of the kissing interaction. To test this hypothesis, we took a two-pronged approach to examine the effects of various mutational and chemical modifications of the flanking adenosines on the folding of the kissing complex. Using mass spectrometry, we studied formation of kissing dimers formed by different hairpins. Using optical tweezers, we monitored mechanical unfolding of intramolecular kissing complex at single-molecule level. In both experiments, replacing adenine with uridine abolished the kissing interaction, suggesting that a minimal kissing complex must contain two GC pairs flanked by inter-strand stacking adenines. The stabilizing effect by the adenines can be explained by the fact that the stacking purine nucleobases shield the hydrogen bonds of the adjacent GC pairs, preventing them from fraying. Unlike in the context of secondary structure, the 5’-unpaired adenines in the tertiary structure are structurally constrained in a way that allows for effective stacking onto the adjacent base pairs.     

Amino Acid Insertion Reveals a Necessary Three-Helical Intermediate in the Folding Pathway of the Colicin E7 Immunity Protein Im7

The small (87-residue) α-helical protein Im7 (an inhibitor protein for colicin E7 that provides immunity to cells producing colicin E7) folds via a three-state mechanism involving an on-pathway intermediate. This kinetic intermediate contains three of four native helices that are oriented in a non-native manner so as to minimise exposed hydrophobic surface area at this point in folding. The short (6-residue) helix III has been shown to be unstructured in the intermediate ensemble and does not dock onto the developing hydrophobic core until after the rate-limiting transition state has been traversed. After helix III has docked, it adopts an α-helical secondary structure, and the side chains of residues within this region provide contacts that are crucial to native-state stability. In order to probe further the role of helix III in the folding mechanism of Im7, we created a variant that contains an eight-amino-acid polyalanine-like helix stabilised by a Glu-Arg salt bridge and an Asn-Pro-Gly capping motif, juxtaposed C-terminal to the natural 6-residue helix III. The effect of this insertion on the structure of the native protein and its folding mechanism were studied using NMR and ?-value analysis, respectively. The results reveal a robust native structure that is not perturbed by the presence of the extended helix III. Mutational analysis performed to probe the folding mechanism of the redesigned protein revealed a conserved mechanism involving the canonical three-helical intermediate. The results suggest that folding via a three-helical species stabilised by both native and non-native interactions is an essential feature of Im7 folding, independent of the helical propensity of helix III.     

Structure of the Tetrahymena thermophila telomerase RNA helix II template boundary element

Telomere addition by telomerase requires an internal templating sequence located in the RNA subunit of telomerase. The correct boundary definition of this template sequence is essential for the proper addition of the nucleotide repeats. Incorporation of incorrect telomeric repeats onto the ends of chromosomes has been shown to induce chromosomal instability in ciliate, yeast and human cells. A 5′ template boundary defining element (TBE) has been identified in human, yeast and ciliate telomerase RNAs. Here, we report the solution structure of the TBE element (helix II) from Tetrahymena thermophila telomerase RNA. Our results indicate that helix II and its capping pentaloop form a well-defined structure including unpaired, stacked adenine nucleotides in the stem and an unusual syn adenine nucleotide in the loop. A comparison of the T.thermophila helix II pentaloop with a pentaloop of the same sequence found in the 23S rRNA of the Haloarcula marismortui ribosome suggests possible RNA and/or protein interactions for the helix II loop within the Tetrahymena telomerase holoenzyme.     

Alpha helix capping in synthetic model peptides by reciprocal side chain–main chain interactions: Evidence for an N terminal “capping box”

A significant fraction of the amino acids in proteins are alpha helical in conformation. Alpha helices in globular proteins are short, with an average length of about twelve residues, so that residues at the ends of helices make up an important fraction of all helical residues. In the middle of a helix, H-bonds connect the NH and CO groups of each residue to partners four residues along the chain. At the ends of a helix, the H-bond potential of the main chain remains unfulfilled, and helix capping interactions involving bonds from polar side chains to the NH or CO of the backbone have been proposed and detected. In a study of synthetic helical peptides, we have found that the sequence Ser-Glu-Asp-Glu stabilizes the alpha helix in a series of helical peptides with consensus sequences. Following the report by Harper and Rose, which identifies SerXaaXaaGlu as a member of a class of common motifs at the N termini of alpha helices in proteins that they refer to as “capping boxes,” we have reexamined the side chain–main chain interactions in a varient sequence using ¹H NMR, and find that the postulated reciprocal side chain-backbone bonding between the first Ser and last Glu side chains and their peptide NH partners can be resolved: Deletion of two residues N terminal to the Ser-Glu-Asp-Glu sequence in these peptides has no effect on the initiation of helical structure, as defined by two-dimensional (2D) NMR experiments on this variant. Thus the capping box sequence Ser-Glu-Asp-Glu inhibits N terminal fraying of the N terminus of alpha helix in these peptides, and shows the side chain–main chain interactions proposed by Harper and Rose. It thus acts as a helix initiating signal. Since normal a helix cannot propagate beyond the N terminus of this structure, the box acts as a termination signal in this direction as well. © 1994 John Wiley & Sons, Inc.     

Crystal Structure of the Bacteriophage Qβ Coat Protein in Complex with the RNA Operator of the Replicase Gene

The coat proteins of single-stranded RNA bacteriophages specifically recognize and bind to a hairpin structure in their genome at the beginning of the replicase gene. The interaction serves to repress the synthesis of the replicase enzyme late in infection and contributes to the specific encapsidation of phage RNA. While this mechanism is conserved throughout the Leviviridae family, the coat protein and operator sequences from different phages show remarkable variation, serving as prime examples for the co-evolution of protein and RNA structure. To better understand the protein–RNA interactions in this virus family, we have determined the three-dimensional structure of the coat protein from bacteriophage Qβ bound to its cognate translational operator. The RNA binding mode of Qβ coat protein shares several features with that of the widely studied phage MS2, but only one nucleotide base in the hairpin loop makes sequence-specific contacts with the protein. Unlike in other RNA phages, the Qβ coat protein does not utilize an adenine-recognition pocket for binding a bulged adenine base in the hairpin stem but instead uses a stacking interaction with a tyrosine side chain to accommodate the base. The extended loop between β strands E and F of Qβ coat protein makes contacts with the lower part of the RNA stem, explaining the greater length dependence of the RNA helix for optimal binding to the protein. Consequently, the complex structure allows the proposal of a mechanism by which the Qβ coat protein recognizes and discriminates in favor of its cognate RNA.     

On loop folding in nucleic acid hairpin-type structures

In a series of studies, combining NMR, optical melting and T-jump experiments, it was found that DNA hairpins display a maximum stability when the loop part of the molecule comprises four or five nucleotide residues. This is in contrast with the current notion based on RNA hairpin studies, from which it had been established that a maximum hairpin stability is obtained for six or seven residues in the loop. Here we present a structural model to rationalize these observations. This model is based on the notion that to a major extent base stacking interactions determine the stability of nucleic acid conformations. The model predicts that loop folding in RNA is characterized by an extension of the base stacking at the 5'-side of the double helix by five or six bases; the remaining gap can then easily be closed by two nucleotides. Conversely, loop folding in DNA is characterized by extending base stacking at the 3'-side of the double helical stem by two or three residues; again bridging of the remaining gap can then be achieved by one or two nucleotides. As an example of loop folding in RNA the anticodon loop of yeast tRNAPhe is discussed. For the DNA hairpin formed by d(ATCCTAT4TAGGAT) it is shown that the loop structure obtained from molecular mechanics calculations obeys the above worded loop folding principles.     

Thermodynamics of unpaired terminal nucleotides on short RNA helixes correlates with stacking at helix termini in larger RNAs.

Free energies for stacking of unpaired nucleotides (dangling ends) at the termini of oligoribonucleotide Watson-Crick helixes (DeltaG(0)37,stack) depend on sequence for 3' ends but are always small for 5' ends. Here, these free energies are correlated with stacking at helix termini in a database of 34 RNA structures determined by X-ray crystallography and NMR spectroscopy. Stacking involving GA pairs is considered separately. A base is categorized as stacked by its distance from (相似文献   

The k-junction motif in RNA structure

The k-junction is a structural motif in RNA comprising a three-way helical junction based upon kink turn (k-turn) architecture. A computer program written to examine relative helical orientation identified the three-way junction of the Arabidopsis TPP riboswitch as an elaborated k-turn. The Escherichia coli TPP riboswitch contains a related k-junction, and analysis of >11 000 sequences shows that the structure is common to these riboswitches. The k-junction exhibits all the key features of an N1-class k-turn, including the standard cross-strand hydrogen bonds. The third helix of the junction is coaxially aligned with the C (canonical) helix, while the k-turn loop forms the turn into the NC (non-canonical) helix. Analysis of ligand binding by ITC and global folding by gel electrophoresis demonstrates the importance of the k-turn nucleotides. Clearly the basic elements of k-turn structure are structurally well suited to generate a three-way helical junction, retaining all the key features and interactions of the k-turn.     

Helix capping in the GCN4 leucine zipper.

Capping interactions associated with specific sequences at or near the ends of alpha-helices are important determinants of the stability of protein secondary and tertiary structure. We investigate here the role of the helix-capping motif Ser-X-X-Glu, a sequence that occurs frequently at the N termini of alpha helices in proteins, on the conformation and stability of the GCN4 leucine zipper. The 1.8 A resolution crystal structure of the capped molecule reveals distinct conformations, packing geometries and hydrogen-bonding networks at the amino terminus of the two helices in the leucine zipper dimer. The free energy of helix stabilization associated with the hydrogen-bonding and hydrophobic interactions in this capping structure is -1.2 kcal/mol, evaluated from thermal unfolding experiments. A single cap thus contributes appreciably to stabilizing the terminated helix and thereby the native state. These results suggest that helix capping plays a further role in protein folding, providing a sensitive connector linking alpha-helix formation to the developing tertiary structure of a protein.     

Effects of alanine substitutions in alpha-helices of sperm whale myoglobin on protein stability.

The peptide backbones in folded native proteins contain distinctive secondary structures, alpha-helices, beta-sheets, and turns, with significant frequency. One question that arises in folding is how the stability of this secondary structure relates to that of the protein as a whole. To address this question, we substituted the alpha-helix-stabilizing alanine side chain at 16 selected sites in the sequence of sperm whale myoglobin, 12 at helical sites on the surface of the protein, and 4 at obviously internal sites. Substitution of alanine for bulky side chains at internal sites destabilizes the protein, as expected if packing interactions are disrupted. Alanine substitutions do not uniformly stabilize the protein, either in capping positions near the ends of helices or at mid-helical sites near the surface of myoglobin. When corrected for the extent of exposure of each side chain replaced by alanine at a mid-helix position, alanine replacement still has no clear effect in stabilizing the native structure. Thus linkage between the stabilization of secondary structure and tertiary structure in myoglobin cannot be demonstrated, probably because of the relatively small free energy differences between side chains in stabilizing isolated helix. By contrast, about 80% of the variance in free energy observed can be accounted for by the loss in buried surface area of the native residue substituted by alanine. The differential free energy of helix stabilization does not account for any additional variation.     

Is there a special function for U.G basepairs in ribosomal RNA?

U.G basepairs are well-established elements of RNA structure. The geometry of this pair is different, however, from classical Watson-Crick basepairs. This leads to an unusual stacking of the basepair: overlap with the basepair at the 5' side of the U (and the 3' side of the G) is strong (stacked) while it is weak with the basepair on the other side (destacked). The closure of an RNA helix by a U.G pair will be energetically unfavourable when the U residue occupies the 5' end. In transfer RNA there is a strong selection against a 'destacked' U.G pair at helix ends. In the 16S rRNA model of Escherichia coli there are 72 U.G pairs of which 36 or 22 occupy a helix end, depending on how such an end is defined. There is a slight preference for 'stacked' U.G's in these positions. It is remarkable, however, that of 13 very conserved U.G pairs in the 16S (-like) rRNA, 7 occur at helix ends and that 5 of these have the 'destacked' configuration. It is suggested that these pairs, if they exist at all in a hydrogen-bounded form, are stabilized by co-axial stacking with other helices or by interaction with a protein.     

Analyzing the evolution of RNA secondary structures in vertebrate introns using Kimura's model of compensatory fitness interactions

Previous studies have shown that splicing efficiency, and thus maturation of pre-mRNA, depends on the correct folding of the RNA molecule into a secondary or higher order structure. When disrupted by a mutation, aberrant folding may result in a lower splicing efficiency. However, the structure can be restored by a second, compensatory mutation. Here, we present a logistic regression approach to analyze the evolutionary dynamics of RNA secondary structures. We apply our approach to a set of computationally predicted RNA secondary structures in vertebrate introns. Our results are consistent with the hypothesis of a negative influence of the physical distance between pairing nucleotides on the occurrence of covariations, as predicted by Kimura's model of compensatory evolution. We also confirm the hypothesis that longer local secondary structure elements (helices) can accommodate a larger number of covariations, wobbles, and mismatches. Furthermore, we find that wobbles and mismatches are more frequent in the middle of a helix, whereas covariations occur preferentially at the helix ends. The GC content is a major determinant of this pattern.     

Analysis of Four-Way Junctions in RNA Structures

RNA secondary structures can be divided into helical regions composed of canonical Watson-Crick and related base pairs, as well as single-stranded regions such as hairpin loops, internal loops, and junctions. These elements function as building blocks in the design of diverse RNA molecules with various fundamental functions in the cell. To better understand the intricate architecture of three-dimensional (3D) RNAs, we analyze existing RNA four-way junctions in terms of base-pair interactions and 3D configurations. Specifically, we identify nine broad junction families according to coaxial stacking patterns and helical configurations. We find that helices within junctions tend to arrange in roughly parallel and perpendicular patterns and stabilize their conformations using common tertiary motifs such as coaxial stacking, loop-helix interaction, and helix packing interaction. Our analysis also reveals a number of highly conserved base-pair interaction patterns and novel tertiary motifs such as A-minor-coaxial stacking combinations and sarcin/ricin motif variants. Such analyses of RNA building blocks can ultimately help in the difficult task of RNA 3D structure prediction.     

A cis-acting replication element in the sequence encoding the NS5B RNA-dependent RNA polymerase is required for hepatitis C virus RNA replication

RNA structures play key roles in the replication of RNA viruses. Sequence alignment software, thermodynamic RNA folding programs, and classical comparative phylogenetic analysis were used to build models of six RNA elements in the coding region of the hepatitis C virus (HCV) RNA-dependent RNA polymerase, NS5B. The importance of five of these elements was evaluated by site-directed mutagenesis of a subgenomic HCV replicon. Mutations disrupting one of the predicted stem-loop structures, designated 5BSL3.2, blocked RNA replication, implicating it as an essential cis-acting replication element (CRE). 5BSL3.2 is about 50 bases in length and is part of a larger predicted cruciform structure (5BSL3). As confirmed by RNA structure probing, 5BSL3.2 consists of an 8-bp lower helix, a 6-bp upper helix, a 12-base terminal loop, and an 8-base internal loop. Mutational analysis and structure probing were used to explore the importance of these features. Primary sequences in the loops were shown to be important for HCV RNA replication, and the upper helix appears to serve as an essential scaffold that helps maintain the overall RNA structure. Unlike certain picornavirus CREs, whose function is position independent, 5BSL3.2 function appears to be context dependent. Understanding the role of 5BSL3.2 and determining how this new CRE functions in the context of previously identified elements at the 5' and 3' ends of the RNA genome should provide new insights into HCV RNA replication.     

All-atom Monte Carlo simulation of GCAA RNA folding

We report a detailed all-atom simulation of the folding of the GCAA RNA tetraloop. The GCAA tetraloop motif is a very common and thermodynamically stable secondary structure in natural RNAs. We use our simulation methods to study the folding behavior of a 12-base GCAA tetraloop structure with a four-base helix adjacent to the tetraloop proper. We implement an all-atom Monte Carlo (MC) simulation of RNA structural dynamics using a Go potential. Molecular dynamics (MD) simulation of RNA and protein has realistic energetics and sterics, but is extremely expensive in terms of computational time. By coarsely treating non-covalent energetics, but retaining all-atom sterics and entropic effects, all-atom MC techniques are a useful method for the study of protein and now RNA. We observe a sharp folding transition for this structure, and in simulations at room temperature the state histogram shows three distinct minima: an unfolded state (U), a more narrow intermediated state (I), and a narrow folded state (F). The intermediate consists primarily of structures with the GCAA loop and some helix hydrogen bonds formed. Repeated kinetic folding simulations reveal that the number of helix base-pairs forms a simple 1D reaction coordinate for the I-->N transition.     

A novel loop-loop recognition motif in the yeast ribosomal protein L30 autoregulatory RNA complex

The yeast Saccharomyces cerevisiae ribosomal protein L30 negatively autoregulates its production by binding to a helix-loop-helix structure formed in its pre-mRNA and its mRNA. A three-dimensional solution structure of the L30 protein in complex with its regulatory RNA has been solved using NMR spectroscopy. In the complex, the helix-loop-helix RNA adopts a sharply bent conformation at the internal loop region. Unusual RNA features include a purine stack, a reverse Hoogsteen base pair (G11anti-G56syn) and highly distorted backbones. The L30 protein is folded in a three-layer alpha/beta/alpha sandwich topology, and three loops at one end of the sandwich make base-specific contacts with the RNA internal loop. The protein-RNA binding interface is divided into two clusters, including hydrophobic and aromatic stacking interactions centering around G56, and base-specific hydrogen-bonding contacts to A57, G58 and G10-U60 wobble base pair. Both the protein and the RNA exhibit a partially induced fit for binding, where loops in the protein and the internal loop in the RNA become more ordered upon complex formation. The specific interactions formed between loops on L30 and the internal loop on the mRNA constitute a novel loop-loop recognition motif where an intimate RNA-protein interface is formed between regions on both molecules that lack regular secondary structure.     

The influence of the local sequence environment on RNA loop structures

RNA folding is assumed to be a hierarchical process. The secondary structure of an RNA molecule, signified by base-pairing and stacking interactions between the paired bases, is formed first. Subsequently, the RNA molecule adopts an energetically favorable three-dimensional conformation in the structural space determined mainly by the rotational degrees of freedom associated with the backbone of regions of unpaired nucleotides (loops). To what extent the backbone conformation of RNA loops also results from interactions within the local sequence context or rather follows global optimization constraints alone has not been addressed yet. Because the majority of base stacking interactions are exerted locally, a critical influence of local sequence on local structure appears plausible. Thus, local loop structure ought to be predictable, at least in part, from the local sequence context alone. To test this hypothesis, we used Random Forests on a nonredundant data set of unpaired nucleotides extracted from 97 X-ray structures from the Protein Data Bank (PDB) to predict discrete backbone angle conformations given by the discretized η/θ-pseudo-torsional space. Predictions on balanced sets with four to six conformational classes using local sequence information yielded average accuracies of up to 55%, thus significantly better than expected by chance (17%-25%). Bases close to the central nucleotide appear to be most tightly linked to its conformation. Our results suggest that RNA loop structure does not only depend on long-range base-pairing interactions; instead, it appears that local sequence context exerts a significant influence on the formation of the local loop structure.