Folding of Intramolecular DNA Hairpin Loops: Enthalpy−Entropy Compensations and Hydration Contributions

Ab initio predictions for 3D structure and stability of single- and double-stranded DNAs in ion solutions

The three-dimensional (3D) structure and stability of DNA are essential to understand/control their biological functions and aid the development of novel materials. In this work, we present a coarse-grained (CG) model for DNA based on the RNA CG model proposed by us, to predict 3D structures and stability for both dsDNA and ssDNA from the sequence. Combined with a Monte Carlo simulated annealing algorithm and CG force fields involving the sequence-dependent base-pairing/stacking interactions and an implicit electrostatic potential, the present model successfully folds 20 dsDNAs (≤52nt) and 20 ssDNAs (≤74nt) into the corresponding native-like structures just from their sequences, with an overall mean RMSD of 3.4Å from the experimental structures. For DNAs with various lengths and sequences, the present model can make reliable predictions on stability, e.g., for 27 dsDNAs with/without bulge/internal loops and 24 ssDNAs including pseudoknot, the mean deviation of predicted melting temperatures from the corresponding experimental data is only ~2.0°C. Furthermore, the model also quantificationally predicts the effects of monovalent or divalent ions on the structure stability of ssDNAs/dsDNAs.

To determine 3D structures and quantify stability of single- (ss) and double-stranded (ds) DNAs is essential to unveil the mechanisms of their functions and to further guide the production and development of novel materials. Although many DNA models have been proposed to reproduce the basic structural, mechanical, or thermodynamic properties of dsDNAs based on the secondary structure information or preset constraints, there are very few models can be used to investigate the ssDNA folding or dsDNA assembly from the sequence. Furthermore, due to the polyanionic nature of DNAs, metal ions (e.g., Na⁺ and Mg²⁺) in solutions can play an essential role in DNA folding and dynamics. Nevertheless, ab initio predictions for DNA folding in ion solutions are still an unresolved problem. In this work, we developed a novel coarse-grained model to predict 3D structures and thermodynamic stabilities for both ssDNAs and dsDNAs in monovalent/divalent ion solutions from their sequences. As compared with the extensive experimental data and available existing models, we showed that the present model can successfully fold simple DNAs into their native-like structures, and can also accurately reproduce the effects of sequence and monovalent/divalent ions on structure stability for ssDNAs including pseudoknot and dsDNAs with/without bulge/internal loops.

Amino Acid Stabilization of Nucleic Acid Secondary Structure: Kinetic Insights from Single-Molecule Studies

Amino acid and nucleic acid interactions are central in biology and may have played a role in the evolutionary development of protein-based life from an early "RNA Universe." To explore the possible role of single amino acids in promoting nucleic acid folding, single-molecule Förster resonance energy transfer experiments have been implemented with a DNA hairpin construct (7 nucleotide double strand with a 40A loop) as a simple model for secondary structure formation. Exposure to positively charged amino acids (arginine and lysine) is found to clearly stabilize the secondary structure. Kinetically, each amino acid promotes folding by generating a large increase in the folding rate with little change in the unfolding rate. From analysis as a function of temperature, arginine and lysine are found to significantly increase the overall exothermicity of folding while imposing only a small entropic penalty on the folding process. Detailed investigations into the kinetics and thermodynamics of this amino acid-induced folding stability reveal arginine and lysine to interact with nucleic acids in a manner reminiscent of monovalent cations. Specifically, these observations are interpreted in the context of an ion atmosphere surrounding the nucleic acid, in which amino acid salts stabilize folding qualitatively like small monovalent cations but also exhibit differences because of the composition of their side chains.

Effect of GCAA stabilizing loops on three- and four-way intramolecular junctions

Tetraloops are a common way of changing the melting behavior of a DNA or RNA structure without changing the sequence of the stem. Because of the ubiquitous nature of tetraloops, our goal is to understand the effect a GCAA tetraloop, which belongs to the GNRA family of tetraloops, has on the unfolding thermodynamics of intramolecular junctions. Specifically, we have described the melting behavior of intramolecular three-way and four-way junctions where a T₅ loop has been replaced with a GCAA tetraloops in different positions. Their thermodynamic profiles, including Δn_Na⁺ and Δn_W, were analyzed based on the position of the tetraloop. We obtained between -16.7 and -27.5 kcal mol^-1 for all junctions studied. The experimental data indicates the influence of the GCAA tetraloop is primarily dictated by the native unfolding of the junction; if the tetraloop is placed on a stem that unfolds as a single domain when the tetraloop is not present, it will unfold as a single domain when the tetraloop is present but with a higher thermal stability. Conversely, if the tetraloop is placed on a stem which unfolds cooperatively with other stems when the tetraloop is not present, the tetraloop will increase the thermal stability of all the stems in the melting domain. The oligonucleotide structure and not the tetraloop itself affects ion uptake; three-way junctions do not gain an increase in ion uptake, but four-way junctions do. This is not the case for water immobilization, where the position of the tetraloop dictates the amount of water immobilized.

Effect of Loop Length and Sequence on the Stability of DNA Pyrimidine Triplexes with TAT Base Triplets

We report the thermodynamic contributions of loop length and loop sequence to the overall stability of DNA intramolecular pyrimidine triplexes. Two sets of triplexes were designed: in the first set, the C₅ loop closing the triplex stem was replaced with ^5'-CT_nC loops (n = 1-5), whereas in the second set, both the duplex and triplex loops were replaced with a ^5'-GCAA or ^5'-AACG tetraloop. For the triplexes with a ^5'-CT_nC loop, the triplex with five bases in the loop has the highest stability relative to the control. A loop length lower than five compromises the strength of the base-pair stacks without decreasing the thermal stability, leading to a decreased enthalpy, whereas an increase in the loop length leads to a decreased enthalpy and a higher entropic penalty. The incorporation of the GCAA loop yielded more stable triplexes, whereas the incorporation of AACG in the triplex loop yielded a less stable triplex due to an unfavorable enthalpy term. Thus, addition of the GCAA tetraloop can cause an increase in the thermodynamics of the triplex without affecting the sequence or melting behavior and may result in an additional layer of genetic regulation.

Investigation of the Melting Behavior of DNA Three-Way Junctions in the Closed and Open States

Intramolecular three-way junctions are commonly found in both DNA and RNA. These structures are functionally relevant in ribozymes, riboswitches, rRNA, and during replication. In this work, we present a thermodynamic description of the unfolding of DNA intramolecular three-way junctions. We used a combination of spectroscopic and calorimetric techniques to investigate the folding/unfolding thermodynamics of two three-way junctions with a closed (Closed-J) or open (Open-J) junction and their appropriate control stem-loop motifs (GAAATT-Hp, CTATC-Hp, and Dumbbell). The overall results show that both junctions are stable over a wide range of salt concentrations. However, Open-J is more stable due to a higher enthalpy contribution from the formation of a higher number of basepair stacks whereas Closed-J has a defined structure and retains the basepair stacking of all three stems. The comparison of the experimental results of Closed-J and Open-J with those of their component stem-loop motifs allowed us to be more specific about their cooperative unfolding. For instance, Closed-J sacrifices thermal stability of the Dumbbell structure to maintain an overall folded state. At higher salt concentration, the simultaneous unfolding of the above domains is lost, resulting in the unfolding of the three separate stems. In contrast, the junction of Open-J in low salt retains the thermal and enthalpic stability of the Dumbbell structure although sacrificing stability of the CTATC stem. The relative stability of Dumbbell is the primary reason for the higher ΔG°₍₅₎, or free energy, value seen for Open-J at low salt. Higher salt not only maintains thermal stability of the Dumbbell structure in Open-J but causes the CTATC stem to fully fold.

Impact of thymine glycol damage on DNA duplex energetics: Correlations with lesion-induced biochemical and structural consequences

The magnitude and nature of lesion-induced energetic perturbations empirically correlate with mutagenicity/cytotoxicity profiles and can be predictive of lesion outcomes during polymerase-mediated replication in vitro. In this study, we assess the sequence and counterbase-dependent energetic impact of the Thymine glycol (Tg) lesion on a family of deoxyoligonucleotide duplexes. Tg damage arises from thymine and methyl-cytosine exposure to oxidizing agents or radiation-generated free-radicals. The Tg lesion blocks polymerase-mediated DNA replication in vitro and the unrepaired site elicits cytotoxic lethal consequences in vivo. Our combined calorimetric and spectroscopic characterization correlates Tg -induced energetic perturbations with biological and structural properties. Specifically, we incorporate a 5R-Tg isomer centered within the tridecanucleotide sequence 5'-GCGTACXCATGCG-3' (X = Tg or T) which is hybridized with the corresponding complementary sequence 5'-CGCATGNGTACGC-3' (N = A, G, T, C) to generate families of Tg -damaged (Tg ·N) and lesion-free (T·N) duplexes. We demonstrate that the magnitude and nature of the Tg destabilizing impact is dependent on counterbase identity (i.e., A ∼ G < T < C). The observation that a Tg lesion is less destabilizing when positioned opposite purines suggests that favorable counterbase stacking interactions may partially compensate lesion-induced perturbations. Moreover, the destabilizing energies of Tg ·N duplexes parallel their respective lesion-free T·N mismatch counterparts (i.e., G < T < C). Elucidation of Tg-induced destabilization relative to the corresponding undamaged mismatch energetics allows resolution of lesion-specific and sequence-dependent impacts. The Tg-induced energetic perturbations are consistent with its replication blocking properties and may serve as differential recognition elements for discrimination by the cellular repair machinery.

The size of the internal loop in DNA hairpins influences their targeting with partially complementary strands

Targeting of noncanonical DNA structures, such as hairpin loops, may have significant diagnostic and therapeutic potential. Oligonucleotides can be used for binding to mRNA, forming a DNA/RNA hybrid duplex that inhibits translation. This kind of modulation of gene expression is called the antisense approach. In order to determine the best strategy to target a common structural motif in mRNA, we have designed a set of stem-loop DNA molecules with sequence: d(GCGCT_nGTAAT₅GTTACT_nGCGC), where n = 1, 3, or 5, “T₅” is an end loop of five thymines. We used a combination of calorimetric and spectroscopy techniques to determine the thermodynamics for the reaction of a set of hairpins containing internal loops with their respective partially complementary strands. Our aim was to determine if internal- and end-loops are promising regions for targeting with their corresponding complementary strands. Indeed, all targeting reactions were accompanied by negative changes in free energy, indicating that reactions proceed spontaneously. Further investigation showed that these negative free energy terms result from a net balance of unfavorable entropy and favorable enthalpy contributions. In particular, unfolding of hairpins and duplexes is accompanied by positive changes in heat capacity, which may be a result of exposure of hydrophobic groups to the solvent. This study provides a new method for the targeting of mRNA in order to control gene expression.

Design and evaluation of an i-motif-based allosteric control mechanism in DNA-hairpin molecular devices

Molecular devices designed to assess and manipulate biologically relevant conditions with required accuracy and precision play an essential role in life sciences research. Incorporating allosteric regulation mechanism is an attractive strategy toward more efficient artificial sensing and switching systems. Herein, we report on a new principle of regulating switching parameters of a DNA-based molecular device based on allosteric interaction between spatially separated hairpin stem and a tetraplexed fragment (i.e., i-motif). We characterized thermodynamic and kinetic effects arising from interaction between functional domains of the device and demonstrated the potential of applying the allosteric control principle for rational design of sensors and switches with precisely defined operational characteristics.

Stable conformations of a single stranded deprotonated DNA i-motif

We present molecular dynamics simulations of a single stranded deprotonated DNA i-motif in explicit solvent. Our results indicate that hairpin structures are stable equilibrium conformations at 300 K. The entropic preference of these configurations is explained by strong water ordering effects due to the present number of hydrogen bonds. We observe a full unfolding at higher temperatures in good agreement with experimental results.

Recognition of thymine in DNA bulges by a Zn(II) macrocyclic complex

A Zn(II) macrocyclic complex with appended quinoline is a bifunctional recognition agent that uses both the Zn(II) center and the pendent aromatic group to bind to thymine in bulges with good selectivity over DNA containing G, C or A bulges. Spectroscopic studies show that the stem containing the bulge stays largely intact in a DNA hairpin with the Zn(II) complex bound to the thymine bulge.

Calorimetric unfolding of intramolecular triplexes: length dependence and incorporation of single AT --> TA substitutions in the duplex domain

DNA triplexes have been the subject of great interest due to their ability to interfere with gene expression. The inhibition of gene expression involves the design of stable triplexes under physiological conditions; therefore, it is important to have a clear understanding of the energetic contributions controlling their stability. We have used a combination of UV spectroscopy and differential scanning calorimetric (DSC) techniques to investigate the unfolding of intramolecular triplexes, d(A(n)C5T(n)C5T(n)), where n is 5-7, 9, and 11, and related triplexes with a single AT --> TA substitution in their duplex stem. Specifically, we obtain standard thermodynamic profiles for the unfolding of each triplex in buffer solutions containing 0.1 M or 1 M NaCl. The triplexes unfold in monophasic or biphasic transitions (triplex --> duplex --> coil) depending on the concentration of salt used and position of the substitution, and their transition temperatures are independent of strand concentration. The DSC curves of the unsubstituted triplexes yielded an unfolding heat of 13.9 kcal/mol for a TAT/TAT base-triplet stack and a heat capacity of 505 cal/ degrees C.mol. The incorporation of a single substitution destabilizes triplex formation (association of the third strand) to a larger extent in 0.1 M NaCl, and the magnitude of the effects also depends on the position of the substitution. The combined results show that a single AT --> TA substitution in a homopurine/homopyrimidine duplex does not allow triplex formation of the neighboring five TAT base triplets, indicating that the in vivo formation of triplexes, such as H-DNA, is exclusive to homopurine/homopyrimidine sequences.

Identification of high-stringency DNA hairpin probes by partial gene folding

Hairpin DNA sequences are widely used as probes for oligonucleotides in a broad range of assays, often as "molecular beacons". A potential disadvantage of the standard methodology for molecular beacon design is the need to add several self-complementary bases to each end of the probe, since these do not correspond to the target sequence. We describe a conceptually new method of hairpin DNA probe identification, in which a secondary structure prediction algorithm is employed to identify oligonucleotide sequences within an expressed gene having the requisite hairpin structure. Intuitively, such probes should have significantly improved performance over "traditional" hairpin probes, because they are fully complementary with the target. We present experimental evidence verifying this hypothesis for a series of hairpin probes targeting the pag gene of Bacillus anthracis.