Role of the closing base pair for d(GCA) hairpin stability: free energy analysis and folding simulations

Quantitative Effects of the Loop Region on Topology, Thermodynamics, and Cation Binding of DNA G-quadruplexes

The thermal stability of G-quadruplexes is important for their biological roles. G-quadruplexes are stable in the presence of cations such as K+ and Na+ because these cations coordinate in the G-quartet of four guanine bases. It is well known that the number of G-quartets and the configuration of the guanine bases affect the binding affinity of the cation. Recently, structures formed in the loop regions connecting the guanine stretches have attracted significant attention, because the loop region affects G-quadruplex properties, such as topology, thermal stability, and interactions with proteins and small molecules. Considering these effects, the loop region can also affect the binding affinity of the cations. Here, we designed a series of G-quadruplex-forming DNA sequences that contain a hairpin in a loop region and investigated the effects of the sequence and structure of the loop region on the cation binding affinity as well as the thermal stability of the G-quadruplex as a whole. First, structural analysis of the DNA sequences showed that the hairpin at the loop plays a key role in determining G4 topology (strand orientation). Second, in the case of the G-quadruplexes with the hairpin-forming loop region, it was found that a longer loop length led to a higher thermodynamic stability of the G-quadruplex as well as higher cation binding affinity. In contrast, an unstructured loop region did not lead to such effects. Interestingly, the cation binding affinity was correlated to the thermodynamic stability of the hairpin structure at the loop region. It was quantitatively demonstrated that the stable loop region stabilized the whole G-quadruplex structure, which induced higher cation binding affinity. These systematic and quantitative results showed that the loop region is one of the determinants of cation binding and expanded the possibilities of drug development targeting G4s by stabilizing the loop region.

Force-induced melting of DNA hairpin: Unfolding pathways and phase diagrams

Using the exact enumeration technique, we have studied the force-induced melting of a DNA hairpin on the face centered cubic lattice for two different sequences which differ in terms of loop closing base pairs. The melting profiles obtained from the exact enumeration technique is consistent with the Gaussian network model and Langevin dynamics simulations. Probability distribution analysis based on the exact density of states revealed the microscopic details of the opening of the hairpin. We showed the existence of intermediate states near the melting temperature. We further showed that different ensembles used to model single-molecule force spectroscopy setups may give different force-temperature diagrams. We delineate the possible reasons for the observed discrepancies.

Influence of Single Deuterium Replacement on Frequency of Hydrogen Bond Dissociation in IFNA17 under the Highest Critical Energy Range

The effect of single substitutions of protium for deuterium in hydrogen bonds between pairs of nitrogenous bases on the open states occurrence probability at high critical breaking energies of these bonds has been studied. The study was carried out using numerical methods based on the angular mathematical model of DNA. The IFNA17 gene was divided into three approximately equal parts. A comparison of the open states occurrence probability in these parts of the gene was done. To improve the accuracy of the results, a special data processing algorithm was developed. The developed methods have shown their suitability for taking into account the occurrence of open states in the entire range of high critical energies. It has been established that single ²H/¹H substitutions in certain nitrogenous bases can be a mechanism for maintaining the vital activity of IFNA17 under critical conditions. In general, the developed method of the mathematical modeling provide unprecedented insight into the DNA behavior under the highest critical energy range, which greatly expands scientific understanding of nucleobases interaction.

The effect of hairpin loop on the structure and gene expression activity of the long-loop G-quadruplex

G-quadruplex (G4), a four-stranded DNA or RNA structure containing stacks of guanine tetrads, plays regulatory roles in many cellular functions. So far, conventional G4s containing loops of 1–7 nucleotides have been widely studied. Increasing experimental evidence suggests that unconventional G4s, such as G4s containing long loops (long-loop G4s), play a regulatory role in the genome by forming a stable structure. Other secondary structures such as hairpins in the loop might thus contribute to the stability of long-loop G4s. Therefore, investigation of the effect of the hairpin-loops on the structure and function of G4s is required. In this study, we performed a systematic biochemical investigation of model G4s containing long loops with various sizes and structures. We found that the long-loop G4s are less stable than conventional G4s, but their stability increased when the loop forms a hairpin (hairpin-G4). We also verified the biological significance of hairpin-G4s by showing that hairpin-G4s present in the genome also form stable G4s and regulate gene expression as confirmed by in cellulo reporter assays. This study contributes to expanding the scope and diversity of G4s, thus facilitating future studies on the role of G4s in the human genome.

Effect of loop sequence on unzipping of short DNA hairpins

The dependence of stability on the sequence of a DNA hairpin has been investigated through atomistic simulations. For this, a sequence of 16 bases of a hairpin, which consists of a loop of four bases and a stem of six base pairs, has been considered. We have taken eight different sequences, where the first five base pairs were kept fixed in all sequences, whereas the loop sequence and the identity of the duplex base pair closing the loop have been varied. For these hairpin structures, force-induced melting (unzipping) studies were carried out to investigate the effect of the variables on the stability of hairpin. The temperature at which half of the base pairs are open is termed the melting temperature. We defined the unzipping force F_{h} (half of the base pairs are open) and showed that it may not provide the effect of closing the base pair or loop sequence on the stability of the DNA hairpin. In order to have a better understanding of the stability of a DNA hairpin, the closing base pair or hairpin loop must be open. This requires complete opening of the stem. We defined a force F_{c} at which all base pairs of the stem are open, and we showed that the F_{c} gives better understanding of DNA hairpin stability.

Dramatic changes in the electrochemical and spectroscopic behavior of DNA fragments induced by the loss of a single nucleotide

In this paper, the d(GCGAAGC) heptamer and the closely related d(GCGAGC) hexamer are examined via electrochemical (cyclic voltammetry) and spectroscopic (circular dichroism) methods. Dramatic changes in the CD spectroscopic and CV electrochemical properties, induced by the loss of only one single nucleotide (A), are detected. The CD spectra and native polyacrylamide gel electrophoresis (PAGE) confirmed structural changes taking place in the relevant chain-like oligodeoxynucleotide assemblies. Dedicated studies suggest that the heptamer (Hp) possesses a hairpin structure, whereas the hexamer (Hx) appears to be rather a duplex. Both of the structures exhibited completely different adsorption behavior at the hanging mercury drop electrode, and this factor was readily confirmed by means of elimination voltammetry with linear scan (EVLS). We established that the Hp hairpin (~-1300 mV), compared to the Hx duplex (~-1360 mV), is the thermodynamically favored electron acceptor. The adsorption isotherms were constructed based on the voltammetric peak height values, reflecting the reduction of the adenine (A) and cytosine (C) moieties as well as the oxidation of the 7,8-dihydroguanine (7,8-DHG) moieties. Finally, as revealed by the spectroscopic and electrochemical results, Hx forms a bimolecular antiparallel homo-duplex carrying both Watson-Crick base pairs (CG or GC) and mismatched edge-to-edge base pairs (GA or AG).

Refinement of Generalized Born Implicit Solvation Parameters for Nucleic Acids and Their Complexes with Proteins

The Generalized Born (GB) implicit solvent model has undergone significant improvements in accuracy for modeling of proteins and small molecules. However, GB still remains a less widely explored option for nucleic acid simulations, in part because fast GB models are often unable to maintain stable nucleic acid structures or they introduce structural bias in proteins, leading to difficulty in application of GB models in simulations of protein-nucleic acid complexes. Recently, GB-neck2 was developed to improve the behavior of protein simulations. In an effort to create a more accurate model for nucleic acids, a similar procedure to the development of GB-neck2 is described here for nucleic acids. The resulting parameter set significantly reduces absolute and relative energy error relative to Poisson-Boltzmann for both nucleic acids and nucleic acid-protein complexes, when compared to its predecessor GB-neck model. This improvement in solvation energy calculation translates to increased structural stability for simulations of DNA and RNA duplexes, quadruplexes, and protein-nucleic acid complexes. The GB-neck2 model also enables successful folding of small DNA and RNA hairpins to near native structures as determined from comparison with experiment. The functional form and all required parameters are provided here and also implemented in the AMBER software.

Unfolding and conformational variations of thrombin-binding DNA aptamers: synthesis, circular dichroism and molecular dynamics simulations

Thrombin-binding DNA aptamer (TBA), with a consensus 15-base sequence: d(GGTTGGTGTGGTTGG), can fold into an antiparallel unimolecular G-quadruplex structure that is necessary for its interaction with thrombin. For the first time, using steered molecular dynamics (SMD) simulations, we have successfully simulated the unfolding process of native TBA G-quadruplex. The unfolding pathway proposed is in agreement with previously reported experimental NMR data. Moreover, the critical intermediate structure in the unfolding pathway, predicted by the NMR results, was identified. The structural characteristics of several TBA oligonucleotides modified with locked nucleoside (LNA) or 2'-O-methyl-nucleoside (MNA) at different positions and number were also investigated by CD spectroscopy. An oligonucleotide substituted with either LNA or MNA at position 2 folds into a native-like G-quadruplex, while doubly substituted derivatives of TBA where LNA or MNA is incorporated at positions 11 and 14 are no longer able to form a G-quadruplex. Starting from the same initial intermediate structure, we successfully overcame sampling limitations, and simulated the large conformational variations in structures of native TBA and modified TBAs by classic MD. Analysis of the models showed that inversion of the glycosyl orientation at position 14 contributes significantly to the disruption of G-quadruplex formation in both of the di-substituted modified TBA systems. Our calculations provide a simple and reliable theoretical model that can be used to investigate and predict the effects of the modifications of an oligonucleotide on the resultant G-quadruplex structure. In addition, the computational protocol described can also be applied for other G-quadruplex systems.

The effect of loop residues in four-stranded dimeric structures stabilized by minor groove tetrads

Some DNA oligonucleotides can fold back and self-associate forming dimeric structures stabilized by intermolecular base pairs. The resulting antiparallel dimer is a tightly packed four-stranded structure formed by a core of minor groove tetrads connected by short loops of unpaired nucleotides. We have explored the sequential requirements for the loop residues and have found that this family of structures is only stable with one- and two-residue loops, with the stability of the former ones being only marginal. Two-residue loops with purines in the first position give rise to the most stable structures due to their enhanced stacking interaction with the adjacent minor groove tetrad. On the other hand, pyrimidines confer more stability than purines in the second position of the loop.

Reliable oligonucleotide conformational ensemble generation in explicit solvent for force field assessment using reservoir replica exchange molecular dynamics simulations

Molecular dynamics force field development and assessment requires a reliable means for obtaining a well-converged conformational ensemble of a molecule in both a time-efficient and cost-effective manner. This remains a challenge for RNA because its rugged energy landscape results in slow conformational sampling and accurate results typically require explicit solvent which increases computational cost. To address this, we performed both traditional and modified replica exchange molecular dynamics simulations on a test system (alanine dipeptide) and an RNA tetramer known to populate A-form-like conformations in solution (single-stranded rGACC). A key focus is on providing the means to demonstrate that convergence is obtained, for example, by investigating replica RMSD profiles and/or detailed ensemble analysis through clustering. We found that traditional replica exchange simulations still require prohibitive time and resource expenditures, even when using GPU accelerated hardware, and our results are not well converged even at 2 μs of simulation time per replica. In contrast, a modified version of replica exchange, reservoir replica exchange in explicit solvent, showed much better convergence and proved to be both a cost-effective and reliable alternative to the traditional approach. We expect this method will be attractive for future research that requires quantitative conformational analysis from explicitly solvated simulations.

Conformational changes of non-B DNA

In contrast to B-DNA that has a right-handed double helical structure with Watson-Crick base pairing under the ordinary physiological conditions, repetitive DNA sequences under certain conditions have the potential to fold into non-B DNA structures such as hairpin, triplex, cruciform, left-handed Z-form, tetraplex, A-motif, etc. Since the non-B DNA-forming sequences induce the genetic instability and consequently can cause human diseases, the molecular mechanism for their genetic instability has been extensively investigated. On the contrary, non-B DNA can be widely used for application in biotechnology because many DNA breakage hotspots are mapped in or near the sequences that have the potential to adopt non-B DNA structures. In addition, they are regarded as a fascinating material for the nanotechnology using non-B DNAs because they do not produce any toxic byproducts and are robust enough for the repetitive working cycle. This being the case, an understanding on the mechanism and dynamics of their structural changes is important. In this critical review, we describe the latest studies on the conformational dynamics of non-B DNAs, with a focus on G-quadruplex, i-motif, Z-DNA, A-motif, hairpin and triplex (189 references).