Tetraloop Hairpin Research Articles

DeGenPrime provides robust primer design and optimization unlocking the biosphere.

Polymerase chain reaction (PCR) is the world's most important molecular diagnostic with applications ranging from medicine to ecology. PCR can fail because of poor primer design. The nearest-neighbor thermodynamic properties, picking conserved regions, and filtration via penalty of oligonucleotides form the basis for good primer design. DeGenPrime is a console-based high-quality PCR primer design tool that can utilize MSA formats and degenerate bases expanding the target range for a single primer set. Our software utilizes thermodynamic properties, filtration metrics, penalty scoring, and conserved region finding of any proposed primer. It has degeneracy, repeated k-mers, relative GC content, and temperature range filters. Minimal penalty scoring is included according to secondary structure self-dimerization metrics, GC clamping, tri- and tetra-loop hairpins, and internal repetition. We compared PrimerDesign-M, DegePrime, ConsensusPrimer, and DeGenPrime on acceptable primer yield. PrimerDesign-M, DegePrime, and ConsensusPrimer provided 0%, 11%, and 17% yield, respectively, for the alternative iron nitrogenase (anfD) gene target. DeGenPrime successfully identified quality primers within the conserved regions of the T4-like phage major capsid protein (g23), conserved regions of molybdenum-based nitrogenase (nif), and its alternatives vanadium (vnf) and iron (anf) nitrogenase. DeGenPrime provides a universal and scalable primer design tool for the entire tree of life. DeGenPrime is written in C++ and distributed under a BSD-3-Clause license. The source code for DeGenPrime is freely available on www.github.com/raw-lab/degenprime.

Multiscale Modeling of Phosphate···π Contacts in RNA U-Turns Exposes Differences between Quantum-Chemical and AMBER Force Field Descriptions.

Phosphate···π, also called anion···π, contacts occur between nucleobases and anionic phosphate oxygens (OP2) in r(GNRA) and r(UNNN) U-turn motifs (N = A,G,C,U; R = A,G). These contacts were investigated using state-of-the-art quantum-chemical methods (QM) to characterize their physicochemical properties and to serve as a reference to evaluate AMBER force field (AFF) performance. We found that phosphate···π interaction energies calculated with the AFF for dimethyl phosphate···nucleobase model systems are less stabilizing in comparison with double-hybrid DFT and that minimum contact distances are larger for all nucleobases. These distance stretches are also observed in large-scale AFF vs QM/MM computations and classical molecular dynamics (MD) simulations on several r(gcGNRAgc) tetraloop hairpins when compared to experimental data extracted from X-ray/cryo-EM structures (res. ≤ 2.5 Å) using the WebFR3D bioinformatic tool. MD simulations further revealed shifted OP2/nucleobase positions. We propose that discrepancies between the QM and AFF result from a combination of missing polarization in the AFF combined with too large AFF Lennard-Jones (LJ) radii of nucleobase carbon atoms in addition to an exaggerated short-range repulsion of the r-12 LJ repulsive term. We compared these results with earlier data gathered on lone pair···π contacts in CpG Z-steps occurring in r(UNCG) tetraloops. In both instances, charge transfer calculations do not support any significant n → π* donation effects. We also investigated thiophosphate···π contacts that showed reduced stabilizing interaction energies when compared to phosphate···π contacts. Thus, we challenge suggestions that the experimentally observed enhanced thermodynamic stability of phosphorothioated r(GNRA) tetraloops can be explained by larger London dispersion.

Structural effects of m6A modification of the Xist A-repeat AUCG tetraloop and its recognition by YTHDC1.

The A-repeat region of the lncRNA Xist is critical for X inactivation and harbors several N6-methyladenosine (m6A) modifications. How the m6A modification affects the conformation of the conserved AUCG tetraloop hairpin of the A-repeats and how it can be recognized by the YTHDC1 reader protein is unknown. Here, we report the NMR solution structure of the (m6A)UCG hairpin, which reveals that the m6A base extends 5′ stacking of the A-form helical stem, resembling the unmethylated AUCG tetraloop. A crystal structure of YTHDC1 bound to the (m6A)UCG tetraloop shows that the (m6A)UC nucleotides are recognized by the YTH domain of YTHDC1 in a single-stranded conformation. The m6A base inserts into the aromatic cage and the U and C bases interact with a flanking charged surface region, resembling the recognition of single-stranded m6A RNA ligands. Notably, NMR and fluorescence quenching experiments show that the binding requires local unfolding of the upper stem region of the (m6A)UCG hairpin. Our data show that m6A can be readily accommodated in hairpin loop regions, but recognition by YTH readers requires local unfolding of flanking stem regions. This suggests how m6A modifications may regulate lncRNA function by modulating RNA structure.

High-throughput thermal stability assessment of DNA hairpins based on high resolution melting

On the basis of high-resolution melting, a high-throughput approach to measure melting temperatures (Tms) of short DNA hairpins was developed. With this method, Tms of thousands of triloop, tetraloop, and pentaloop hairpins involving various loop sequences and various closing base pairs (cbp) were obtained in hours. The stability of triloop hairpins decreased with the change of cbp (5′–3′) in the order of c-g > g-c > t-a ≥ a-t, showing that the cbp of 5′-Pyr-Pur-3′ (Pyr = pyrimidine, Pur = purine) contributed more stability than 5′-Pur-Pyr-3′. For tetraloop hairpins, GNNA, GNAB, and CNNG (N = A, G, C, or T; B = G, C, or T) were found to be highly stable irrespective of the cbp type. TNNA was also stable in both g-c and a-t families, while CGNA only in the c-g family. Pentaloop hairpins of cTGNAGg, cGNYNAg (Y = T or C) and cCGNNAg were exceptionally stable motifs. In most cases, pyrimidine-rich loops were more favorable to stabilize the whole structure than purine-rich ones. The present approach showed a good performance in assessing the thermal stability of large amounts of DNA hairpins comprehensively. These data are useful to understand the sequence dependence of the stability of DNA secondary structures and promising to improve the structure simulation by consummating basic databases.

Molecular simulation of cyclohexanyl nucleic acid (CNA) duplexes with CNA, DNA and RNA and CNA triloop and tetraloop hairpin structures

As part of a selection strategy for artificial nucleic acids (XNA) (to be considered as potential new information systems in vivo), we have carried out a modelling study on cyclohexanyl nucleic acids (CNA) duplexes and hairpins. CNA may form a duplex as well as hairpin structures, having the carbocyclic nucleosides in the 4C1 conformation (with equatorial basis). The geometry of ds CNA is close to that of a HNA:RNA duplex. We demonstrated that CNA triphosphates function as a substrate for polymerases. Modelling experiments indicate that the monomers are probably presented to the polymerase in the 1C4 conformation.

Energy landscapes, folding mechanisms, and kinetics of RNA tetraloop hairpins.

RNA hairpins play a pivotal role in a diverse range of cellular functions, and are integral components of ribozymes, mRNA, and riboswitches. However, the mechanistic and kinetic details of RNA hairpin folding, which are key determinants of most of its biological functions, are poorly understood. In this work, we use the discrete path sampling (DPS) approach to explore the energy landscapes of two RNA tetraloop hairpins, and provide insights into their folding mechanisms and kinetics in atomistic detail. Our results show that the potential energy landscapes have a distinct funnel-like bias toward the folded hairpin state, consistent with efficient structure-seeking properties. Mechanistic and kinetic information is analyzed in terms of kinetic transition networks. We find microsecond folding times, consistent with temperature jump experiments, for hairpin folding initiated from relatively compact unfolded states. This process is essentially driven by an initial collapse, followed by rapid zippering of the helix stem in the final phase. Much lower folding rates are predicted when the folding is initiated from extended chains, which undergo longer excursions on the energy landscape before nucleation events can occur. Our work therefore explains recent experiments and coarse-grained simulations, where the folding kinetics exhibit precisely this dependency on the initial conditions.

Tospovirus ambisense genomic RNA segments use almost complete repertoire of stable tetraloops in the intergenic region

The intergenic regions of the ambisense RNA segments of viruses from the Tospovirus genus form large extended RNA structures that regulate virus replication. Using comparative structure analysis, we show the presence of conserved alternative conformations at the apical parts of these structures. In one conformation, a branched Y-shape, the 5'-proximal hairpin arms are mostly capped by exceptionally stable tetraloop motifs. The tetraloop hairpins are folded in both virus and virus-complementary sense RNAs, and different tetraloops can functionally replace each other. Folding simulations show that the branched Y-shape structures can undergo a conformational transition to alternative extended rod-like conformations. Functional importance of both alternatives is supported by nucleotide covariations. The balanced equilibrium between alternative structures is evidenced by native gel electrophoresis of mutant RNA transcripts with shifted equilibria. The tetraloops play a role in the stability and dynamics of structures but may also be recognized by proteins involved in translation and/or replication.

Combining temperature and force to study folding of an RNA hairpin

RNA folding in cells typically occurs at mesophilic temperatures. However, in vitro, RNA can be unfolded either by increasing temperature to values that are much higher than physiological, or by mechanically pulling structures apart at ambient temperature. To directly study RNA folding at physiological temperatures and to unify thermodynamics measured by melting and pulling, we developed temperature-controlled optical tweezers (thermal tweezers) that can be used to mechanically unfold single RNA molecules at mesophilic temperatures. Folding of a 20-base-pair tetraloop hairpin was studied under different ionic conditions and at temperatures ranging from 22 °C to 42 °C. At each temperature, single hairpin molecules were held at constant force, and their two-state folding equilibria were monitored. The change in free energy derived from these measurements was used to construct a phase diagram of RNA structure using force and temperature as variables. Furthermore, we derived ΔG(0pN,T), the folding free energy at zero force and temperature T, by subtracting the stretching energy of unfolded RNA from the reversible mechanical work done to unfold the hairpin. ΔG(0pN,T) and its salt dependence agree reasonably well with the predictions by the nearest neighbor model. Under each ionic condition, ΔG(0pN,T) depended linearly on temperature, yielding ΔH(exp) and ΔS(exp) that also matched the predictions. The combination of force and temperature to study RNA folding is a step toward unifying thermodynamics measured by thermal melting and mechanical unfolding, and opens a new path for directly monitoring temperature induced RNA structural changes, as it occurs often in biology.

Investigating the Folding Dynamics of RNA Pseudoknot Structural Motif via Massively Parallel Molecular Dynamics

RNA pseudoknots compose a three-dimensional structural motif that is present in the catalytic cores of some ribozymes, and are also capable of stimulating ribosomal frameshifts. Their complex topology and non-canonical hairpin-loop composition make pseudoknots an ideal structural motif with which to study the RNA folding process. Here we report our analysis of nearly 20,000 independent all-atom molecular dynamics simulations of the ribosomal frame-shifting pseudoknot of Luteovirus and the tmRNA pseudoknot from Aquifex aeolicus, which share global topology but have only ∼50% sequence similarity. Using the Folding@Home distributed computing network and a novel Pathway Enumeration sampling method, a cumulative sampling time of over 115 μs was achieved for each of these pseudoknots. K-means clustering identified 27 conformational microstates for each pseudoknot, which reached conformational equilibrium after ∼6 ns of ensemble sampling. Multiple folding metrics were used to identify 9 macrostates participating in the folding process, including previously undescribed misfolded and intermediate states. In agreement with our previous studies of tetraloop hairpins and tRNA, the similar folding behavior of these pseudoknots suggests that native state topology is a predominant factor in the RNA folding mechanism.

Counterion and Polythymidine Loop-Length-Dependent Folding and Thermodynamic Stability of DNA Hairpins Reveal the Unusual Counterion-Dependent Stability of Tetraloop Hairpins

Stem-loop DNA hairpins containing a 5-base-pair (bp) stem and single-stranded polythymidine loop were investigated using thermodynamic melting analysis and stopped-flow kinetics. These studies revealed the thermodynamic stability and folding kinetics as a function of loop length and counterion concentration. Our results show the unusually high thermodynamic stability for tetraloop or 4 poly(dT) loop hairpin as compared with longer loop length hairpins. Furthermore, this exceptional stability is highly counterion-dependent. For example, in the higher counterion concentration regime of 50 mM NaCl and above, the tetraloop hairpin displays enhanced stability as compared with longer loop length hairpins. However, at lower counterion concentration of 25 mM NaCl and below, the thermal stability of tetraloop hairpin is consistent with the longer loop hairpins. The enhanced stability of tetraloop hairpins at higher counterion concentration can be explained on the basis of the combined entropic effect of loop closure as well as base stacking in the loop regions. The stability of longer loop length hairpins at all counterion concentrations as well as tetraloop hairpin at lower counterion concentration can be explained on the basis of entropic effect of loop closure alone. The thermodynamic parameters at lower and higher counterion concentrations were determined to quantify the enhanced stability of base-stacking effects occurring at higher counterion concentrations. For example, for 100 mM NaCl, excess Gibbs energy and enthalpy due to base stacking within the tetraloops were measured to be -1.2 ± 0.14 and -3.28 ± 0.32 kcal/mol, respectively, whereas, no excess of Gibbs energy and enthalpy was observed for 0, 5, 10, and 25 mM NaCl. These findings suggest significant base-stacking interactions occurring in the loop region of the tetraloop hairpins at higher counterion concentration and less significant base-stacking interactions in the lower counterion concentration regime. We suggest that at higher counterion concentrations, hydrophobic collapse of the nucleotides in the loop may be enhanced due to the increased polarity of the solvent, thereby enhancing base-stacking interactions that contribute to unusually high stability.

The effects of loop size on Sac7d-hairpin DNA interactions

Hairpin structure is a common feature of DNA molecules. They are located near functional loci, such as regulation and promotion sites, as well as in cruciform structures, and they provide potential binding sites for endogenous proteins. The effects of different hairpin loops that are composed of one to five thymidines, designated as L1–L5, and have a common self-complementary stem, CTATATAG, on the interactions with Sac7d were studied. In thermostability studies, Sac7d stabilized a tetra-loop hairpin DNA and hairpin DNA with GTTC tetra-loop regions better than it stabilized tri- and penta-loops. Circular dichroism (CD) spectra showed that hairpins retained primarily a B-type conformation upon Sac7d binding. Intermolecular interactions between hairpins were likely decreased, due to the Sac7d-induced kinks, as shown by an increase at 220nm in the CD spectra. Surface plasmon resonance (SPR) observations suggested that the rates of Sac7d binding to hairpin DNA depend on the loop size of the hairpin duplexes. At a fixed stem length, Sac7d binds to tetra-loop hairpin DNA duplexes with a higher association rate and lower dissociation rate, compared with their tri- and penta-loop counterparts. In addition, the tri-loop and GTC tri-loop hairpin DNA had lower affinity for Sac7d because of the smaller and tighter loop size. Our study indicates that Sac7d binding affinity to hairpin DNA is primarily determined by loop size and stem integrity, and the results presented here provide a model for studies concerning other minor groove DNA-binding proteins that kink hairpin DNA.

The Conformationally Constrained N-Methanocarba-dT Analogue Adopts an Unexpected C4′-exo Sugar Pucker in the Structure of a DNA Hairpin

Incorporation of a bicyclo[3.1.0]hexane scaffold into the nucleoside sugar was devised to lock the embedded cyclopentane ring in conformations that mimic the furanose North and South sugar puckers. To analyze the effects of North-methanocarba-2'-deoxythymidine (N-MCdT) on the B-form DNA, we crystallized d(CGCGAA[mcTmcT]CGCG) with two N-MCdTs. Instead of a duplex, the 12mer forms a tetraloop hairpin, whereby loop N-MCdTs adopt the C4'-exo pucker (NE; P = 50°). Thus, the bicyclic framework does not limit the pucker to the anticipated C2'-exo range (NNW; P = -18°).

Structure of a human Xist RNA A-repeat AUCG tetraloop hairpin essential for X-inactivation

Structure of a human Xist RNA A-repeat AUCG tetraloop hairpin essential for X-inactivation

Folding Kinetics for the Conformational Switch Between Alternative RNA Structures

The conformational switching between different conformational states is intrinsic to RNA catalytic and regulatory functions, which oftenly occurs on time-scales of several seconds. In combination with the recent real-time NMR experiments (Wenter et al. Angew. Chem. Int. Ed. (2005). 44, 2600; Wenter et al. ChemBioChem. (2006). 7, 417) for the transitions between bistable RNA conformations, we combine the master equation method with the kinetic cluster method to investigate the detailed kinetic mechanism and the factors that govern the folding kinetics. Based on the computational studies, we propose that heat capacity change upon RNA folding may be important for RNA folding kinetics. In addition, we find that noncanonical (tertiary) intraloop interactions in tetraloop hairpins are important to determine the folding kinetics. Furthermore, through theory-experiment comparisons, we find that the different rate models for the fundamental steps (i.e., formation/disruption of a base pair or stack) can cause contrasting results in the theoretical predictions.

Folding kinetics for the conformational switch between alternative RNA structures.

Transitions between different conformational states, so-called conformational switching, are intrinsic to RNA catalytic and regulatory functions. Often, conformational switching occurs on time scales of several seconds. In combination with the recent real-time NMR experiments (Wenter et al. Angew. Chem. Int. Ed. 2005, 44, 2600; Wenter et al. ChemBioChem 2006, 7, 417) for the transitions between bistable RNA conformations, we combine the master equation method with the kinetic cluster method to investigate the detailed kinetic mechanism and the factors that govern the folding kinetics. We propose that heat capacity change (ΔC(p)) upon RNA folding may be important for RNA folding kinetics. In addition, we find that, for tetraloop hairpins, noncanonical (tertiary) intraloop interactions are important to determine the folding kinetics. Furthermore, through theory-experiment comparisons, we find that the different rate models for the fundamental steps (i.e., formation/disruption of a base pair or stack) can cause contrasting results in the theoretical predictions.

Folding of a Small RNA Hairpin Based on Simulation with Replica Exchange Molecular Dynamics

The folding of a small RNA tetraloop hairpin is studied based on intensive molecular dynamics simulation, aiming to understand the folding mechanism of this small and fast RNA folder. Our results showed that this RNA hairpin has very complicated folding behavior in spite of its small size. It is found that the folding transition has low cooperativity. Instead of a two-state folding, four major states are observed, including the native state, the intermediate, the unfolded state, and the misfolded state. The misfolded state is mainly stabilized by the non-native hydrogen bonds, and is more compact. Two potential folding pathways, in which two basepairs formed with different order, are observed, and the pathway with the inboard basepair formed before the terminal one is much more favorable, and dominates the folding of the RNA hairpin.

High-resolution NMR structure of an RNA model system: the 14-mer cUUCGg tetraloop hairpin RNA

We present a high-resolution nuclear magnetic resonance (NMR) solution structure of a 14-mer RNA hairpin capped by cUUCGg tetraloop. This short and very stable RNA presents an important model system for the study of RNA structure and dynamics using NMR spectroscopy, molecular dynamics (MD) simulations and RNA force-field development. The extraordinary high precision of the structure (root mean square deviation of 0.3 Å) could be achieved by measuring and incorporating all currently accessible NMR parameters, including distances derived from nuclear Overhauser effect (NOE) intensities, torsion-angle dependent homonuclear and heteronuclear scalar coupling constants, projection-angle-dependent cross-correlated relaxation rates and residual dipolar couplings. The structure calculations were performed with the program CNS using the ARIA setup and protocols. The structure quality was further improved by a final refinement in explicit water using OPLS force field parameters for non-bonded interactions and charges. In addition, the 2′-hydroxyl groups have been assigned and their conformation has been analyzed based on NOE contacts. The structure currently defines a benchmark for the precision and accuracy amenable to RNA structure determination by NMR spectroscopy. Here, we discuss the impact of various NMR restraints on structure quality and discuss in detail the dynamics of this system as previously determined.

Bending and Internal Tension of an Intramolecular Kissing Complex

RNA tertiary base pairs are often bent and twisted by relative orientation of the interacting domains and by structures connecting these domains. Kissing complex, a basic RNA tertiary interaction, is formed by base pairing between distal hairpin loops. It remains unclear how much structure distortion a kissing complex can withstand. Crystal and NMR structures of several RNA kissing dimers show either a head-to-head coaxial helix or a slightly bent structure. An adenine riboswitch forms a well-bent two-base-pair kissing complex around its ligand. We designed two RNAs both consisting of a pair of hairpins connected by a single-strand. The two tetraloop hairpins in each RNA can form two GC kissing base pairs. The two RNAs, KC15 and KC30, have 15- and 30-nucleotide single-stranded linkers, respectively. Using optical tweezers, we examined mechanical folding of the two structures. With applied force, both RNAs can form an intramolecular kissing complex but show distinct mechanical unfolding patterns. Based on measured size and stability of the structures, we propose a mechanical model that such intramolecular kissing complexes, particularly the KC15, are bent to a blunt angle by internal tension of the single strand. Formation of such RNA triangle requires externally applied force by optical tweezers to offset the internal tension. Measurement of internal tension and structure distortion may facilitate understanding and prediction of tertiary interaction.

Structural Insight into RNA Hairpin Folding Intermediates

Hairpins are a ubiquitous secondary structure motif in RNA molecules. Despite their simple structure, there is some debate over whether they fold in a two-state or multi-state manner. We have studied the folding of a small tetraloop hairpin using a serial version of replica exchange molecular dynamics on a distributed computing environment. On the basis of these simulations, we have identified a number of intermediates that are consistent with experimental results. We also find that folding is not simply the reverse of high-temperature unfolding and suggest that this may be a general feature of biomolecular folding.

Recognition of a Flipped Base in a Hairpinloop DNA by a Small Peptide

Two tiny hairpin DNAs, CORE (dAGGCTTCGGCCT) and AP2 (dAGGCTXCGGCCT; X: abasic nucleotide), fold into almost the same tetraloop hairpin structure with one exception, that is, the sixth thymine (T6) of CORE is exposed to the solvent water (Kawakami, J. et al., Chem. Lett. 2001, 258–259). In the present study, we selected small peptides that bind to CORE or AP2 from a combinatorial pentapeptide library with 2.5 × 106 variants. On the basis of the structural information, the selected peptide sequences should indicate the essential qualifications for recognition of the hairpin loop DNA with and without a flipped base. In the selected DNA binding peptides, aromatic amino acids such as histidine for CORE and glutamine/aspartic acid for AP2 were found to be abundant amino acids. This amino acid preference suggests that CORE-binding peptides use π–π stacking to recognize the target while hydrogen bonding is dominant for AP2-binding peptides. To investigate the binding properties of the selected peptide to the target, surface plasmon resonance was used. The binding constant of the interaction between CORE and a CORE-binding peptide (HWHHE) was about 1.1 × 106 M−1 at 25°C and the resulting binding free energy change at 25°C (ΔG°25) was −8.2 kcal mol−1. The binding of the peptide to AP2 was also analyzed and the resulting binding constant and ΔG°25 were about 4.2 × 104 M−1 and −6.3 kcal mol−1, respectively. The difference in the binding free energy changes (ΔΔG°25) of 1.9 kcal mol−1 was comparable to the values reported in other systems and was considered a consequence of the loss of π–π stacking. Moreover, the stabilization effect by stacking affected the dissociation step as well as the association step. Our results suggest that the existence of an aromatic ring (T6 base) produces new dominant interactions between peptides and nucleic acids, although hydrogen bonding is the preferable mode of interaction in the absence of the flipping base. These findings regarding CORE and AP2 recognition are expected to give useful information in the design of novel artificial DNA binding peptides.