Duplex-tetraplex equilibria in guanine- and cytosine-rich DNA

In-cell NMR suggests that DNA i-motif levels are strongly depleted in living human cells

I-Motifs (iM) are non-canonical DNA structures potentially forming in the accessible, single-stranded, cytosine-rich genomic regions with regulatory roles. Chromatin, protein interactions, and intracellular properties seem to govern iM formation at sites with i-motif formation propensity (iMFPS) in human cells, yet their specific contributions remain unclear. Using in-cell NMR with oligonucleotide iMFPS models, we monitor iM-associated structural equilibria in asynchronous and cell cycle-synchronized HeLa cells at 37 °C. Our findings show that iMFPS displaying pHT < 7 under reference in vitro conditions occur predominantly in unfolded states in cells, while those with pHT > 7 appear as a mix of folded and unfolded states depending on the cell cycle phase. Comparing these results with previous data obtained using an iM-specific antibody (iMab) reveals that cell cycle-dependent iM formation has a dual origin, and iM formation concerns only a tiny fraction (possibly 1%) of genomic sites with iM formation propensity. We propose a comprehensive model aligning observations from iMab and in-cell NMR and enabling the identification of iMFPS capable of adopting iM structures under physiological conditions in living human cells. Our results suggest that many iMFPS may have biological roles linked to their unfolded states.

Structural Transitions in Complementary G-Rich and C-Rich Strands and Their Mixture at Various pH Conditions

We used circular dichroism spectroscopy, UV spectrophotometry, and differential scanning calorimetry to investigate pH-dependent structural transitions in an equimolar mixture of complementary G-rich d[5'-A(GGGTTA)3GGG-3'] (TelG) and C-rich d[3'-T(CCCAAT)3CCC-5'] (TelC) human telomeric DNA strands. Our studies were conducted at neutral (pH 7.0) and slightly acidic (pH 5.5 and 6.5) pH. We analyzed the melting thermodynamics of TelG and TelC and their equimolar mixture. Our analysis revealed that the preferred conformation of an equimolar mixture of TelG and TelC is the duplex. At pH 5.5, however, in addition to the duplex state, we observed a significant population of the i-motif state formed by TelC. Our results are consistent with the picture in which an increase in pH from 5.5 to 7.0 has little effect on the melting enthalpy of an isolated G-quadruplex while causing a strong reduction in the melting enthalpy of an isolated i-motif (the latter diminishes to 0 at pH 7.0). These effects summarily lead to a decrease in the contribution of the i-motif to the melting enthalpy of the mixture and, hence, an increase in the apparent melting enthalpy and overall stability of the duplex state.

Transition Mutations in the hTERT Promoter Are Unrelated to Potential i-motif Formation in the C-Rich Strand

Increased expression of the human telomere reverse transcriptase (hTERT) in tumors promotes tumor cell survival and diminishes the survival of patients. Cytosine-to-thymine (C-to-T) transition mutations (C250T or C228T) in the hTERT promoter create binding sites for transcription factors, which enhance transcription. The G-rich strand of the hTERT promoter can form G-quadruplex structures, whereas the C-rich strand can form an i-motif in which multiple cytosine residues are protonated. We considered the possibility that i-motif formation might promote cytosine deamination to uracil and C-to-T mutations. We computationally probed the accessibility of cytosine residues in an i-motif to attack by water. We experimentally examined regions of the C-rich strand to form i-motifs using pH-dependent UV and CD spectra. We then incubated the C-rich strand with and without the G-rich complementary strand DNA under various conditions, followed by deep sequencing. Surprisingly, deamination rates did not vary substantially across the 46 cytosines examined, and the two mutation hotspots were not deamination hotspots. The appearance of mutational hotspots in tumors is more likely the result of the selection of sequences with increased promoter binding affinity and hTERT expression.

Kinetic and Thermodynamic Control of G-Quadruplex Polymorphism by Na+ and K+ Cations

G-Quadruplexes (G4s) are ubiquitous nucleic acid folding motifs that exhibit structural diversity that is dependent on cationic conditions. In this work, we exploit temperature-controlled single-molecule fluorescence resonance energy transfer (smFRET) to elucidate the kinetic and thermodynamic mechanisms by which monovalent cations (K+ and Na+) impact folding topologies for a simple G-quadruplex sequence (5'-GGG-(TAAGGG)3-3') with a three-state folding equilibrium. Kinetic measurements indicate that Na+ and K+ influence G4 formation in two distinctly different ways: the presence of Na+ modestly enhances an antiparallel G4 topology through an induced fit (IF) mechanism with a low affinity (Kd = 228 ± 26 mM), while K+ drives G4 into a parallel/hybrid topology via a conformational selection (CS) mechanism with much higher affinity (Kd = 1.9 ± 0.2 mM). Additionally, temperature-dependent studies of folding rate constants and equilibrium ratios reveal distinctly different thermodynamic driving forces behind G4 binding to K+ (ΔH°bind > 0, ΔS°bind > 0) versus Na+ (ΔH°bind < 0, ΔS°bind < 0), which further illuminates the diversity of the possible pathways for monovalent facilitation of G-quadruplex folding.

Assessing the Potential for DNA Quadruplex Formation in the Predatory Bacterium Bdellovibrio bacteriovorus

During its life cycle, the predatory bacterium Bdellovibrio bacteriovorus switches between an attack and a growth phase, each of which is characterized by a distinct pattern of gene expression. Twenty-one potential G-quadruplex-forming sequences (PQFS) have been identified in the Bdellovibrio genome. These G-rich sequences are prevalent within open reading frames and nearly evenly distributed between the template and the coding strand, suggesting that they could play a role in gene expression and life cycle switching. Published transcriptomic data show that the genes nearest these sequences are not (de)activated together during the same phases of the life cycle. We explored the biophysical properties of three identified PQFS using circular dichroism (CD) spectroscopy and gel electrophoresis and demonstrated that all three sequences fold into stable unimolecular quadruplexes with distinct topologies. In the presence of their complementary strands, each forms an equilibrium mixture of duplex and quadruplex in which quadruplex formation is favored at higher temperatures. Once the quadruplexes are folded, they are slow to form a duplex when the complementary strand is added, with one sequence requiring the equivalent of many Bdellovibrio lifetimes to do so. Using a variety of cosolutes, we showed that molecular crowding mimicking cellular conditions stabilizes the quadruplex structures and induces structural transitions to the parallel topology regardless of the original topology. Taken together, these experiments suggest that Bdellovibrio PQFS are capable of forming quadruplexes in vivo and thereby playing a role in gene expression.

Distribution of Conformational States Adopted by DNA from the Promoter Regions of the VEGF and Bcl-2 Oncogenes

We employed a previously described procedure, based on circular dichroism (CD) spectroscopy, to quantify the distribution of conformational states adopted by equimolar mixtures of complementary G-rich and C-rich DNA strands from the promoter regions of the VEGF and Bcl-2 oncogenes. Spectra were recorded at different pHs, concentrations of KCl, and temperatures. The temperature dependences of the fractional populations of the duplex, G-quadruplex, i-motif, and coiled conformations of each promoter were then analyzed within the framework of a thermodynamic model to obtain the enthalpy and melting temperature of each folded-to-unfolded transition involved in the equilibrium. A comparison of the conformational data on the VEGF and Bcl-2 DNA with similar results on the c-MYC DNA, which we reported previously, reveals that the distribution of conformational states depends on the specific DNA sequence and is modulated by environmental factors. Under the physiological conditions of room temperature, neutral pH, and elevated concentrations of potassium ions, the duplex conformation coexists with the G-quadruplex conformation in proportions that depend on the sequence. This observed conformational diversity has biological implications, and it further supports our previously proposed thermodynamic hypothesis of gene regulation. In that hypothesis, a specific distribution of duplex and tetraplex conformations in a promoter region is fine-tuned to maintain the healthy level of gene expression. Any deviation from a healthy distribution of conformational states may result in pathology stemming from up- or downregulation of the gene.

Modulation of Aptamer-Ligand-Binding by Complementary Oligonucleotides: A G-Quadruplex Anti-Ochratoxin A Aptamer Case Study

Short oligonucleotides are widely used for the construction of aptamer-based sensors and logical bioelements to modulate aptamer-ligand binding. However, relationships between the parameters (length, location of the complementary region) of oligonucleotides and their influence on aptamer-ligand interactions remain unclear. Here, we addressed this task by comparing the effects of short complementary oligonucleotides (ssDNAs) on the structure and ligand-binding ability of an aptamer and identifying ssDNAs' features that determine these effects. Within this, the interactions between the OTA-specific G-quadruplex aptamer 1.12.2 (5'-GATCGGGTGTGGGTGGCGTAAAGGGA GCATCGGACA-3') and 21 single-stranded DNA (ssDNA) oligonucleotides complementary to different regions of the aptamer were studied. Two sets of aptamer-ssDNA dissociation constants were obtained in the absence and in the presence of OTA by isothermal calorimetry and fluorescence anisotropy, respectively. In both sets, the binding constants depend on the number of hydrogen bonds formed in the aptamer-ssDNA complex. The ssDNAs' having more than 23 hydrogen bonds with the aptamer have a lower aptamer dissociation constant than for aptamer-OTA interactions. The ssDNAs' having less than 18 hydrogen bonds did not affect the aptamer-OTA affinity. The location of ssDNA's complementary site in the aptamer affeced the kinetics of the interaction and retention of OTA-binding in aptamer-ssDNA complexes. The location of the ssDNA site in the aptamer G-quadruplex led to its unfolding. In the presence of OTA, the unfolding process was longer and takes from 20 to 70 min. The refolding in the presence of OTA was possible and depends on the length and location of the ssDNA's complementary site. The location of the ssDNA site in the tail region led to its rapid displacement and wasn't affecting the G-qaudruplex's integrity. It makes the tail region more perspective for the development of ssDNA-based tools using this aptamer.

Characterization of a G-quadruplex from hepatitis B virus and its stabilization by binding TMPyP4, BRACO19 and PhenDC3

Specific guanine rich nucleic acid sequences can form non-canonical structures, like the four stranded G-quadruplex (GQ). We studied the GQ-forming sequence (named HepB) found in the genome of the hepatitis B virus. Fluorescence-, infrared- and CD-spectroscopy were used. HepB shows a hybrid form in presence of K⁺, but Na⁺, Li⁺, and Rb⁺ induce parallel structure. Higher concentrations of metal ions increase the unfolding temperature, which was explained by a short thermodynamic calculation. Temperature stability of the GQ structure was determined for all these ions. Na⁺ has stronger stabilizing effect on HepB than K⁺, which is highly unusual. The transition temperatures were 56.6, 53.8, 58.5 and 54.4 °C for Na⁺, K⁺, Li⁺, and Rb⁺ respectively. Binding constants for Na⁺ and K⁺ were 10.2 mM and 7.1 mM respectively. Study of three ligands designed in cancer research for GQ targeting (TMPyP4, BRACO19 and PhenDC3) showed unequivocally their binding to HepB. Binding was proven by the increased stability of the bound form. The stabilization was higher than 20 °C for TMPyP4 and PhenDC3, while it was considerably lower for BRACO19. These results might have medical importance in the fight against the hepatitis B virus.

Effect of DNA modifications on the transition between canonical and non-canonical DNA structures in CpG islands during senescence

Patterns and levels of DNA modifications play important roles in senescence. Two major epigenetic modifications of DNA, 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC), target CpG sites. Importantly, CpG concentrated regions, known as CpG islands, contain GC-rich sequences, which have the potential to fold into non-canonical DNA structures such as i-motifs and G-quadruplexes. In this study, we investigated the effect of 5mC and 5hmC modifications on the transition between a duplex, and i-motif and G-quadruplexes. To examine the transition, we firstly investigated the stability and structure of the i-motif and G-quadruplexes, considering the molecular environment in senescent cells. Analyses of their stability showed that the modifications did not drastically affect the stability. However, noteworthily, the modification can weaken the (de)stabilisation effect on G-quadruplexes caused by cosolute(s) and cations. Circular dichroism analyses indicated that the surrounding environments, including the molecular crowding and the type of cations such as K⁺ and Na⁺, regulate the topology of G-quadruplexes, while neither 5mC nor 5hmC had a drastic effect. On the other hand, the modifications changed the transition between duplexes and quadruplexes. Unmodified DNA preferred to fold into quadruplexes, whereas DNA with 5mC and 5hmC preferred to fold into duplexes in the absence of PEG200; on the other hand, DNA with or without modifications tended to fold into i-motifs under crowded conditions. Furthermore, an investigation of quadruplexes forming sequences in CpG islands, which are hyper- or hypomethylated during senescence, followed by gene ontology enrichment analysis for each gene group classified by the presence of quadruplexes, showed a difference in function between genes with and without quadruplexes in the CpG region. These results indicate that it is important to consider the effects of patterns and levels of DNA modifications on the transition between canonical and non-canonical DNA structures to understand gene regulation by epigenetic modification during senescence.

The modification of DNA can regulate the transition between a duplex and quadruplexes during senescence responding to surrounding environments.

Oxygen vacancy assisted condensation of DNA molecule observed on ZnO thin film

An exotic condensation of DNA molecules is observed on the nanostructured ZnO surface. The ZnO nanostructures (NS) fabricated by thermal vapor deposition technique were associated with a large number of oxygen vacancies on the surface. These oxygen vacancies induced changes in the DNA conformation which further reflected through changes in the persistence length of the DNA molecules. This indicates a reinforcement of the bonds and binding in both the phosphate and the base regions of the DNA molecules with the positively charged core vacancy sites on the ZnO nanostructured surface through strong interaction mediated via long-range electrostatic forces which effectively reduced the end-to-end distance of the λ-DNA molecule. This strongly suggests a transition of the λ-DNA molecule through structural modification into a more compact higher-order fractal dimension from its native state.

Volumetric Properties of Four-Stranded DNA Structures

Four-stranded non-canonical DNA structures including G-quadruplexes and i-motifs have been found in the genome and are thought to be involved in regulation of biological function. These structures have been implicated in telomere biology, genomic instability, and regulation of transcription and translation events. To gain an understanding of the molecular determinants underlying the biological role of four-stranded DNA structures, their biophysical properties have been extensively studied. The limited libraries on volume, expansibility, and compressibility accumulated to date have begun to provide insights into the molecular origins of helix-to-coil and helix-to-helix conformational transitions involving four-stranded DNA structures. In this article, we review the recent progress in volumetric investigations of G-quadruplexes and i-motifs, emphasizing how such data can be used to characterize intra-and intermolecular interactions, including solvation. We describe how volumetric data can be interpreted at the molecular level to yield a better understanding of the role that solute-solvent interactions play in modulating the stability and recognition events of nucleic acids. Taken together, volumetric studies facilitate unveiling the molecular determinants of biological events involving biopolymers, including G-quadruplexes and i-motifs, by providing one more piece to the thermodynamic puzzle describing the energetics of cellular processes in vitro and, by extension, in vivo.

Untangling the interaction of α-synuclein with DNA i-motifs and hairpins by volume-sensitive single-molecule FRET spectroscopy

The intrinsically disordered protein α-synuclein causes Parkinson's disease by forming toxic oligomeric aggregates inside neurons. Single-molecule FRET experiments revealed conformational changes of noncanonical DNA structures, such as i-motifs and hairpins, in the presence of α-synuclein. Volumetric analyses revealed differences in binding mode, which is also affected by cellular osmolytes.

The conformational landscape of noncanonical DNA structures is markedly affected by monomeric and aggregated α-synuclein, by osmolytes such as TMAO, and by pressure.

Volumetric Interplay between the Conformational States Adopted by Guanine-Rich DNA from the c-MYC Promoter

The kinetic and thermodynamic stabilities of G-quadruplex structures have been extensively studied. In contrast, systematic investigations of the volumetric properties of G-quadruplexes determining their pressure stability are still relatively scarce. The G-rich strand from the promoter region of the c-MYC oncogene (G-strand) is known to adopt a range of conformational states including the duplex, G-quadruplex, and coil states depending on the presence of the complementary C-rich strand (C-strand) and solution conditions. In this work, we report changes in volume, ΔV, and adiabatic compressibility, ΔK_S, accompanying interconversions of G-strand between the G-quadruplex, duplex, and coil conformations in the presence and absence of C-strand. We rationalize these volumetric characteristics in terms of the hydration and intrinsic properties of the DNA in each of the sampled conformational states. We further use our volumetric results in conjunction with the reported data on changes in expansibility, ΔE, and heat capacity, ΔC_P, associated with G-quadruplex-to-coil transitions to construct the pressure-temperature phase diagram describing the stability of the G-quadruplex. The phase diagram is elliptic in shape, resembling the classical elliptic phase diagram of a globular protein, and is distinct from the phase diagram for duplex DNA. The observed similarity of the pressure-temperature phase diagrams of G-quadruplexes and globular proteins stems from their shared structural and hydration features that, in turn, result in the similarity of their volumetric properties. To the best of our knowledge, this is the first pressure-temperature stability diagram reported for a G-quadruplex.

Folding and persistence times of intramolecular G-quadruplexes transiently embedded in a DNA duplex

G-quadruplex (G4) DNA structures have emerged as important regulatory elements during DNA metabolic transactions. While many in vitro studies have focused on the kinetics of G4 formation within DNA single-strands, G4 are found in vivo in double-stranded DNA regions, where their formation is challenged by the complementary strand. Since the energy of hybridization of Watson-Crick structures dominates the energy of G4 folding, this competition should play a critical role on G4 persistence. To address this, we designed a single-molecule assay allowing to measure G4 folding and persistence times in the presence of the complementary strand. We quantified both folding and unfolding rates of biologically relevant G4 sequences, such as the cMYC and cKIT oncogene promoters, human telomeres and an avian replication origin. We confirmed that G4s are found much more stable in tested replication origin and promoters than in human telomere repeats. In addition, we characterized how G4 dynamics was affected by G4 ligands and showed that both folding rate and persistence time increased. Our assay opens new perspectives for the measurement of G4 dynamics in double-stranded DNA mimicking a replication fork, which is important to understand their role in DNA replication and gene regulation at a mechanistic level.

Unraveling the binding characteristics of small ligands to telomeric DNA by pressure modulation

Recently, non-canonical DNA structures, such as G-quadruplexes (GQs), were found to be highly pressure sensitive, suggesting that pressure modulation studies can provide additional mechanistic details of such biomolecular systems. Using FRET and CD spectroscopy as well as binding equilibrium measurements, we investigated the effect of pressure on the binding reaction of the ligand ThT to the quadruplex 22AG in solutions containing different ionic species and a crowding agent mimicking the intracellular milieu. Pressure modulation helped us to identify the different conformational substates adopted by the quadruplex at the different solution conditions and to determine the volumetric changes during complex formation and the conformational transitions involved. The magnitudes of the binding volumes are a hallmark of packing defects and hydrational changes upon ligand binding. The conformational substates of the GQ as well as the binding strength and the stoichiometry of complex formation depend strongly on the solution conditions as well as on pressure. High hydrostatic pressure can also impact GQs inside living cells and thus affect expression of genetic information in deep sea organisms. We show that sub-kbar pressures do not only affect the conformational dynamics and structures of GQs, but also their ligand binding reactions.

Watson-Crick versus Hoogsteen Base Pairs: Chemical Strategy to Encode and Express Genetic Information in Life

Nucleic acids typically form a double helix structure through Watson-Crick base-pairing. In contrast, non-Watson-Crick base pairs can form other three-dimensional structures. Although it is well-known that Watson-Crick base pairs may be more unstable than non-Watson-Crick base pairs under some conditions, the importance of non-Watson-Crick base pairs has not been widely examined. Hoogsteen base pairs, the non-Watson-Crick base pairs, contain important hydrogen-bond patterns that form the helices of nucleic acids, such as in Watson-Crick base pairs, and can form non-double helix structures such as triplexes and quadruplexes. In recent years, non-double helix structures have been discovered in cells and were reported to considerably influence gene expression. The complex behavior of these nucleic acids in cells is gradually being revealed, but the underlying mechanisms remain almost unknown.Quantitatively analyzing the structural stability of nucleic acids is important for understanding their behavior. A nucleic acid is an anionic biopolymer composed of a sugar, base, and phosphoric acid. The physicochemical factors that determine the stability of nucleic acid structures include those derived from the interactions of nucleic acid structures and those derived from the environments surrounding nucleic acids. The Gibbs free energy change (ΔG) of structure formation is the most commonly used physicochemical parameter for analyzing quantitative stability. Quantitatively understanding the intracellular behavior of nucleic acids involves describing the formation of nucleic acid structures and related reactions as ΔG. Based on this concept, we quantitatively analyzed the stability of double helix and non-double helix structures and found that decreased water activity, an important factor in crowded cellular conditions, significantly destabilize the formation of Watson-Crick base pairs but stabilizes Hoogsteen base pairs.Here, we describe a physicochemical approach to understand the regulation of gene expressions based on the stability of nucleic acid structures. We developed new methods for predicting the stability of double and non-double helices in various molecular environments by mimicking intracellular environments. Furthermore, the physicochemical approach used for analyzing gene expression regulated by non-double helix structures is useful for not only determining how gene expression is controlled by cellular environments but also for developing new technologies to chemically regulate gene expression by targeting non-double helix structures. We discuss the roles of Watson-Crick and Hoogsteen base pairs in cells based on our results and why both types of base pairing are required for life. Finally, a new concept in nucleic acid science beyond that of Watson and Crick base pairing is introduced.

Marshall's nucleic acid: From double-helical structure to a potent intercalator

Deoxyribonucleic acid (DNA) not only stores genetic information but also emerged as a popular drug target. Modified nucleotides/nucleosides have been extensively studied in recent years wherein the sugar/nucleobase/phosphate-backbone has been altered. Several such molecules are FDA approved, capable of targeting nucleic acids and proteins. In this article, we modified negatively charged phosphate backbone to marshall's acid-based neutral backbone and analyzed the resultant structures by utilizing Gaussian accelerated molecular dynamics simulations (1 μs) in aqueous media at 150 mM salt concentration. We noted that the double-helical marshall's nucleic acid structure was partially denatured during the course of simulations, however, after using conformationally locked sugar, the marshall's nucleic acid (hereby called MNA) maintained the double-helical structure throughout the simulations. Despite the fact that MNA has a more extended backbone than the regular DNA, surprisingly, both showed similar helical rise (~3.4 Å) along with a comparable Watson-Crick hydrogen bond profile. The backbone difference was majorly compensated in terms of helical twist (~56° (MNA) and ~ 35° (control DNA)). Further, we examined a few MNA based ss-dinucleotides as intercalating ligands for a regular B-DNA. Quite strikingly, the ligands unwinded the DNA and showed intercalating properties with high DNA binding affinities. Hence, the use of small fragments of MNA based molecules in DNA targeted drug discovery is foreseen.