Electrostatic free energy landscapes for nucleic acid helix assembly

Metal ions are crucial for nucleic acid folding. From the free energy landscapes, we investigate the detailed mechanism for ion-induced collapse for a paradigm system: loop-tethered short DNA helices. We find that Na⁺ and Mg²⁺ play distinctive roles in helix–helix assembly. High [Na⁺] (>0.3 M) causes a reduced helix–helix electrostatic repulsion and a subsequent disordered packing of helices. In contrast, Mg²⁺ of concentration >1 mM is predicted to induce helix–helix attraction and results in a more compact and ordered helix–helix packing. Mg²⁺ is much more efficient in causing nucleic acid compaction. In addition, the free energy landscape shows that the tethering loops between the helices also play a significant role. A flexible loop, such as a neutral loop or a polynucleotide loop in high salt concentration, enhances the close approach of the helices in order to gain the loop entropy. On the other hand, a rigid loop, such as a polynucleotide loop in low salt concentration, tends to de-compact the helices. Therefore, a polynucleotide loop significantly enhances the sharpness of the ion-induced compaction transition. Moreover, we find that a larger number of helices in the system or a smaller radius of the divalent ions can cause a more abrupt compaction transition and a more compact state at high ion concentration, and the ion size effect becomes more pronounced as the number of helices is increased.     

Conformations of an RNA Helix-Junction-Helix Construct Revealed by SAXS Refinement of MD Simulations

RNA is involved in a broad range of biological processes that extend far beyond translation. Many of RNA’s recently discovered functions rely on folding to a specific conformation or transitioning between conformations. The RNA structure contains rigid, short basepaired regions connected by more flexible linkers. Studies of model constructs such as small helix-junction-helix (HJH) motifs are useful in understanding how these elements work together to determine RNA conformation. Here, we reveal the full ensemble of solution structures assumed by a model RNA HJH. We apply small-angle x-ray scattering and an ensemble optimization method to selectively refine models generated by all-atom molecular dynamics simulations. The expectation of a broad distribution of helix orientations, at and above physiological ionic strength, is not met. Instead, this analysis shows that the HJH structures are dominated by two distinct conformations at moderate to high ionic strength. Atomic structures, selected from the molecular dynamics simulations, reveal strong base-base interactions in the junction that critically constrain the conformational space available to the HJH molecule and lead to a surprising re-extension at high salt. These results are corroborated by comparison with previous single-molecule fluorescence resonance energy transfer experiments on the same constructs.     

Three-way RNA junctions with remote tertiary contacts: A recurrent and highly versatile fold

Three-way junction RNAs adopt a recurrent Y shape when two of the helices form a coaxial stack and the third helix establishes one or more tertiary contacts several base pairs away from the junction. In this review, the structure, distribution, and functional relevance of these motifs are examined. Structurally, the folds exhibit conserved junction topologies, and the distal tertiary interactions play a crucial role in determining the final shape of the structures. The junctions and remote tertiary contacts behave as flexible hinge motifs that respond to changes in the other region, providing these folds with switching mechanisms that have been shown to be functionally useful in a variety of contexts. In addition, the juxtaposition of RNA domains at the junction and at the distal tertiary complexes enables the RNA helices to adopt unusual conformations that are frequently used by proteins, RNA molecules, and antibiotics as platforms for specific binding. As a consequence of these properties, Y-shaped junctions are widely distributed in all kingdoms of life, having been observed in small naked RNAs such as riboswitches and ribozymes or embedded in complex ribonucleoprotein systems like ribosomal RNAs, RNase P, or the signal recognition particle. In all cases, the folds were found to play an essential role for the functioning or assembly of the RNA or ribonucleoprotein systems that contain them.     

New insights into the fundamental role of topological constraints as a determinant of two-way junction conformation

Recent studies have shown that topological constraints encoded at the RNA secondary structure level involving basic steric and stereochemical forces can significantly restrict the orientations sampled by helices across two-way RNA junctions. Here, we formulate these topological constraints in greater quantitative detail and use this topological framework to rationalize long-standing but poorly understood observations regarding the basic behavior of RNA two-way junctions. Notably, we show that the asymmetric nature of the A-form helix and the finite length of a bulge provide a physical basis for the experimentally observed directionality and bulge-length amplitude dependence of bulge induced inter-helical bends. We also find that the topologically allowed space can be modulated by variations in sequence, particularly with the addition of non-canonical GU base pairs at the junction, and, surprisingly, by the length of the 5' and 3' helices. A survey of two-way RNA junctions in the protein data bank confirms that junction residues have a strong preference to adopt looped-in, non-canonically base-paired conformations, providing a route for extending our bulge-directed framework to internal loop motifs and implying a simplified link between secondary and tertiary structure. Finally, our results uncover a new simple mechanism for coupling junction-induced topological constraints with tertiary interactions.     

The Snakelike Chain Character of Unstructured RNA

In the absence of base-pairing and tertiary structure, ribonucleic acid (RNA) assumes a random-walk conformation, modulated by the electrostatic self-repulsion of the charged, flexible backbone. This behavior is often modeled as a Kratky-Porod “wormlike chain” (WLC) with a Barrat-Joanny scale-dependent persistence length. In this study we report measurements of the end-to-end extension of poly(U) RNA under 0.1 to 10 pN applied force and observe two distinct elastic-response regimes: a low-force, power-law regime characteristic of a chain of swollen blobs on long length scales and a high-force, salt-valence-dependent regime consistent with ion-stabilized crumpling on short length scales. This short-scale structure is additionally supported by force- and salt-dependent quantification of the RNA ion atmosphere composition, which shows that ions are liberated under stretching; the number of ions liberated increases with increasing bulk salt concentration. Both this result and the observation of two elastic-response regimes directly contradict the WLC model, which predicts a single elastic regime across all forces and, when accounting for scale-dependent persistence length, the opposite trend in ion release with salt concentration. We conclude that RNA is better described as a “snakelike chain,” characterized by smooth bending on long length scales and ion-stabilized crumpling on short length scales. In monovalent salt, these two regimes are separated by a characteristic length that scales with the Debye screening length, highlighting the determining importance of electrostatics in RNA conformation.     

3D maps of RNA interhelical junctions

More than 50% of RNA secondary structure is estimated to be A-form helices, which are linked together by various junctions. Here we describe a protocol for computing three interhelical Euler angles describing the relative orientation of helices across RNA junctions. 5' and 3' helices, H1 and H2, respectively, are assigned based on the junction topology. A reference canonical helix is constructed using an appropriate molecular builder software consisting of two continuous idealized A-form helices (iH1 and iH2) with helix axis oriented along the molecular Z-direction running toward the positive direction from iH1 to iH2. The phosphate groups and the carbon and oxygen atoms of the sugars are used to superimpose helix H1 of a target interhelical junction onto the corresponding iH1 of the reference helix. A copy of iH2 is then superimposed onto the resulting H2 helix to generate iH2'. A rotation matrix R is computed, which rotates iH2' into iH2 and expresses the rotation parameters in terms of three Euler angles α(h), β(h) and γ(h). The angles are processed to resolve a twofold degeneracy and to select an overall rotation around the axis of the reference helix. The three interhelical Euler angles define clockwise rotations around the 5' (-γ(h)) and 3' (α(h)) helices and an interhelical bend angle (β(h)). The angles can be depicted graphically to provide a 'Ramachandran'-type view of RNA global structure that can be used to identify unusual conformations as well as to understand variations due to changes in sequence, junction topology and other parameters.     

The Lipid-binding Domain of Wild Type and Mutant α-Synuclein: COMPACTNESS AND INTERCONVERSION BETWEEN THE BROKEN AND EXTENDED HELIX FORMS*

α-Synuclein (αS) is linked to Parkinson disease through its deposition in an amyloid fibril form within Lewy Body deposits, and by the existence of three αS point mutations that lead to early onset autosomal dominant Parkinsonism. The normal function of αS is thought to be linked to the ability of the protein to bind to the surface of synaptic vesicles. Upon binding to vesicles, αS undergoes a structural reorganization from a dynamic and disordered ensemble to a conformation consisting of a long extended helix. In the presence of small spheroidal detergent micelles, however, this extended helix conformation can convert into a broken helix state, in which a region near the middle of the helix unwinds to form a linker between the two resulting separated helices. Membrane-bound conformations of αS likely mediate the function of the protein, but may also play a role in the aggregation and toxicity of the protein. Here we have undertaken a study of the effects of the three known PD-linked mutations on the detergent- and membrane-bound conformations of αS, as well as factors that govern the transition of the protein between the extended helix and broken helix states. Using pulsed dipolar ESR measurements of distances up to 8.7 nm, we show that all three PD-linked αS mutants retain the ability to transition from the broken helix to the extended helix conformation. In addition, we find that the ratio of protein to detergent, rather than just the absolute detergent concentration, determines whether the protein adopts the broken or extended helix conformation.     

Triple helical DNA in a duplex context and base pair opening

It is fundamental to explore in atomic detail the behavior of DNA triple helices as a means to understand the role they might play in vivo and to better engineer their use in genetic technologies, such as antigene therapy. To this aim we have performed atomistic simulations of a purine-rich antiparallel triple helix stretch of 10 base triplets flanked by canonical Watson–Crick double helices. At the same time we have explored the thermodynamic behavior of a flipping Watson–Crick base pair in the context of the triple and double helix. The third strand can be accommodated in a B-like duplex conformation. Upon binding, the double helix changes shape, and becomes more rigid. The triple-helical region increases its major groove width mainly by oversliding in the negative direction. The resulting conformations are somewhere between the A and B conformations with base pairs remaining almost perpendicular to the helical axis. The neighboring duplex regions maintain a B DNA conformation. Base pair opening in the duplex regions is more probable than in the triplex and binding of the Hoogsteen strand does not influence base pair breathing in the neighboring duplex region.     

Do conformational biases of simple helical junctions influence RNA folding stability and specificity?

Structured RNAs must fold into their native structures and discriminate against a large number of alternative ones, an especially difficult task given the limited information content of RNA''s nucleotide alphabet. The simplest motifs within structured RNAs are two helices joined by nonhelical junctions. To uncover the fundamental behavior of these motifs and to elucidate the underlying physical forces and challenges faced by structured RNAs, we computationally and experimentally studied a tethered duplex model system composed of two helices joined by flexible single- or double-stranded polyethylene glycol tethers, whose lengths correspond to those typically observed in junctions from structured RNAs. To dissect the thermodynamic properties of these simple motifs, we computationally probed how junction topology, electrostatics, and tertiary contact location influenced folding stability. Small-angle X-ray scattering was used to assess our predictions. Single- or double-stranded junctions, independent of sequence, greatly reduce the space of allowed helical conformations and influencing the preferred location and orientation of their adjoining helices. A double-stranded junction guides the helices along a hinge-like pathway. In contrast, a single-stranded junction samples a broader set of conformations and has different preferences than the double-stranded junction. In turn, these preferences determine the stability and distinct specificities of tertiary structure formation. These sequence-independent effects suggest that properties as simple as a junction''s topology can generally define the accessible conformational space, thereby stabilizing desired structures and assisting in discriminating against misfolded structures. Thus, junction topology provides a fundamental strategy for transcending the limitations imposed by the low information content of RNA primary sequence.     

Helix capping.

Helix-capping motifs are specific patterns of hydrogen bonding and hydrophobic interactions found at or near the ends of helices in both proteins and peptides. In an alpha-helix, the first four >N-H groups and last four >C=O groups necessarily lack intrahelical hydrogen bonds. Instead, such groups are often capped by alternative hydrogen bond partners. This review enlarges our earlier hypothesis (Presta LG, Rose GD. 1988. Helix signals in proteins. Science 240:1632-1641) to include hydrophobic capping. A hydrophobic interaction that straddles the helix terminus is always associated with hydrogen-bonded capping. From a global survey among proteins of known structure, seven distinct capping motifs are identified-three at the helix N-terminus and four at the C-terminus. The consensus sequence patterns of these seven motifs, together with results from simple molecular modeling, are used to formulate useful rules of thumb for helix termination. Finally, we examine the role of helix capping as a bridge linking the conformation of secondary structure to supersecondary structure.     

RNA helix stability in mixed Na+/Mg2+ solution

A recently developed tightly bound ion model can account for the correlation and fluctuation (i.e., different binding modes) of bound ions. However, the model cannot treat mixed ion solutions, which are physiologically relevant and biologically significant, and the model was based on B-DNA helices and thus cannot directly treat RNA helices. In the present study, we investigate the effects of ion correlation and fluctuation on the thermodynamic stability of finite length RNA helices immersed in a mixed solution of monovalent and divalent ions. Experimental comparisons demonstrate that the model gives improved predictions over the Poisson-Boltzmann theory, which has been found to underestimate the roles of multivalent ions such as Mg2+ in stabilizing DNA and RNA helices. The tightly bound ion model makes quantitative predictions on how the Na+-Mg2+ competition determines helix stability and its helix length-dependence. In addition, the model gives empirical formulas for the thermodynamic parameters as functions of Na+/Mg2+ concentrations and helix length. Such formulas can be quite useful for practical applications.     

Solution structural studies on human erythrocyte alpha-spectrin tetramerization site

We have determined the solution NMR structure of a recombinant peptide that consists of the first 156 residues of erythroid alpha-spectrin. The first 20 residues preceding the first helix (helix C') are in a disordered conformation. The subsequent three helices (helices A1, B1, and C1) form a triple helical bundle structural domain that is similar, but not identical, to previously published structures for spectrin from Drosophila and chicken brain. Paramagnetic spin label-induced NMR resonance broadening shows that helix C', the partial domain involved in alpha- and beta-spectrin association, exhibits little interaction with the structural domain. Surprisingly, helix C' is connected to helix A1 of the structural domain by a segment of 7 residues (the junction region) that exhibits a flexible disordered conformation, in contrast to the predicted rigid helical structure. We suggest that the flexibility of this particular junction region may play an important role in modulating the association affinity of alpha- and beta-spectrin at the tetramerization site of different isoforms, such as erythroid spectrin and brain spectrin. These findings may provide insight for explaining various physiological and pathological conditions that are a consequence of varying alpha- and beta-subunit self-association affinities in their formation of the various spectrin tetramers.     

Dynamic Motions of the HIV-1 Frameshift Site RNA

The HIV-1 frameshift site (FS) plays a critical role in viral replication. During translation, the HIV-1 FS transitions from a 3-helix to a 2-helix junction RNA secondary structure. The 2-helix junction structure contains a GGA bulge, and purine-rich bulges are common motifs in RNA secondary structure. Here, we investigate the dynamics of the HIV-1 FS 2-helix junction RNA. Interhelical motions were studied under different ionic conditions using NMR order tensor analysis of residual dipolar couplings. In 150 mM potassium, the RNA adopts a 43°(±4°) interhelical bend angle (β) and displays large amplitude, anisotropic interhelical motions characterized by a 0.52(±0.04) internal generalized degree of order (GDO_int) and distinct order tensor asymmetries for its two helices (η = 0.26(±0.04) and 0.5(±0.1)). These motions are effectively quenched by addition of 2 mM magnesium (GDO_int = 0.87(±0.06)), which promotes a near-coaxial conformation (β = 15°(±6°)) of the two helices. Base stacking in the bulge was investigated using the fluorescent purine analog 2-aminopurine. These results indicate that magnesium stabilizes extrahelical conformations of the bulge nucleotides, thereby promoting coaxial stacking of helices. These results are highly similar to previous studies of the HIV transactivation response RNA, despite a complete lack of sequence similarity between the two RNAs. Thus, the conformational space of these RNAs is largely determined by the topology of their interhelical junctions.     

Conformation of human apolipoprotein C-I in a lipid-mimetic environment determined by CD and NMR spectroscopy.

The high-resolution conformation of human apoC-I in complexes with sodium dodecyl sulfate (SDS) is presented. As estimated from CD data, apoC-I adopts 54% helical secondary structure when bound to SDS, which is similar to the helical content previously found with phospholipids. The NMR-derived conformation of apoC-I is composed of two amphipathic helices, residues 7-29 and 38-52, separated by a flexible linker. The N-terminal helix contains a mobile hinge involving residues 12-15. The hydrophobic side chains cluster on the nonpolar face of both helices, thus forming two discrete lipid-binding sites in the N-terminal helix and one in the C-terminal helix. As suggested by amide proton resonance line widths and deuterium exchange rates, the N-terminal helix is more flexible and may bind less tightly to the detergent than the C-terminal helix. The different mobility of both helices appears to be related to side-chain composition, rather than length of the amphipathic helix, and may play a role in the function of apoC-I as an activator of lecithin:cholesterol acyltransferase (LCAT). A model is suggested in which the C-terminal helix serves as a lipid anchor while the N-terminal helix may hinge off the lipid surface to make specific contacts with LCAT.     

Effects of single-base bulges on intercalator binding to small RNA and DNA hairpins and a ribosomal RNA fragment

The way in which a single-base bulge might affect the structure of an RNA helix has been examined by preparing a series of six RNA hairpins, all with seven base pairs and a four-nucleotide loop. Five of the hairpins have single-base bulges at different positions. The intercalating cleavage reagent (methidiumpropyl)-EDTA-Fe(II) [MPE-Fe(II)] binds preferentially at a CpG sequence in the helix lacking a bulge and in four of the five hairpins with bulges. Hairpins with a bulge one or two bases to the 3' side of the CpG sequence bind ethidium 4-5-fold more strongly than the others. V1 RNase, which is sensitive to RNA backbone conformation in helices, detects a conformational change in all of the helices when ethidium binds; the most dramatic changes, involving the entire hairpin stem, are in one of the two hairpins with enhanced ethidium affinity. Only a slight conformational change is detected in the hairpin lacking a bulge. A bulge adjacent to a CpG sequence in a 100-nucleotide ribosomal RNA fragment enhances MPE-Fe(II) binding by an order of magnitude. These results extend our previous observations of bulges at a single position in an RNA hairpin [White, S. A., & Draper, D.E. (1987) Nucleic Acids Res. 15, 4049] and show that (1) a structural change in an RNA helix may be propagated for several base pairs, (2) bulges tend to increase the number of conformations available to a helix, and (3) the effects observed in small RNA hairpins are relevant to larger RNAs with more extensive structure. A bulge in a DNA hairpin identical in sequence with the RNA hairpins does not enhance MPE-Fe(II) binding affinity, relative to a control DNA hairpin. The effects of bulges on ethidium intercalation are evidently modulated by helix structure.     

Determination of thermodynamic parameters for HIV DIS type loop-loop kissing complexes

The HIV-1 type dimerization initiation signal (DIS) loop was used as a starting point for the analysis of the stability of Watson–Crick (WC) base pairs in a tertiary structure context. We used ultraviolet melting to determine thermodynamic parameters for loop–loop tertiary interactions and compared them with regular secondary structure RNA helices of the same sequences. In 1 M Na⁺ the loop–loop interaction of a HIV-1 DIS type pairing is 4 kcal/mol more stable than its sequence in an equivalent regular and isolated RNA helix. This difference is constant and sequence independent, suggesting that the rules governing the stability of WC base pairs in the secondary structure context are also valid for WC base pairs in the tertiary structure context. Moreover, the effect of ion concentration on the stability of loop–loop tertiary interactions differs considerably from that of regular RNA helices. The stabilization by Na⁺ and Mg²⁺ is significantly greater if the base pairing occurs within the context of a loop–loop interaction. The dependence of the structural stability on salt concentration was defined via the slope of a T_m/log [ion] plot. The short base-paired helices are stabilized by 8°C/log [Mg²⁺] or 11°C/log [Na⁺], whereas base-paired helices forming tertiary loop–loop interactions are stabilized by 16°C/log [Mg²⁺] and 26°C/log [Na⁺]. The different dependence on ionic strength that is observed might reflect the contribution of specific divalent ion binding to the preformation of the hairpin loops poised for the tertiary kissing loop–loop contacts.     

The Snakelike Chain Character of Unstructured RNA

In the absence of base-pairing and tertiary structure, ribonucleic acid (RNA) assumes a random-walk conformation, modulated by the electrostatic self-repulsion of the charged, flexible backbone. This behavior is often modeled as a Kratky-Porod “wormlike chain” (WLC) with a Barrat-Joanny scale-dependent persistence length. In this study we report measurements of the end-to-end extension of poly(U) RNA under 0.1 to 10 pN applied force and observe two distinct elastic-response regimes: a low-force, power-law regime characteristic of a chain of swollen blobs on long length scales and a high-force, salt-valence-dependent regime consistent with ion-stabilized crumpling on short length scales. This short-scale structure is additionally supported by force- and salt-dependent quantification of the RNA ion atmosphere composition, which shows that ions are liberated under stretching; the number of ions liberated increases with increasing bulk salt concentration. Both this result and the observation of two elastic-response regimes directly contradict the WLC model, which predicts a single elastic regime across all forces and, when accounting for scale-dependent persistence length, the opposite trend in ion release with salt concentration. We conclude that RNA is better described as a “snakelike chain,” characterized by smooth bending on long length scales and ion-stabilized crumpling on short length scales. In monovalent salt, these two regimes are separated by a characteristic length that scales with the Debye screening length, highlighting the determining importance of electrostatics in RNA conformation.