Next Article in Journal
Neuropathic and Psychogenic Components of Burning Mouth Syndrome: A Systematic Review
Next Article in Special Issue
Impact of G-Quadruplexes on the Regulation of Genome Integrity, DNA Damage and Repair
Previous Article in Journal
Perineuronal Nets and Metal Cation Concentrations in the Microenvironments of Fast-Spiking, Parvalbumin-Expressing GABAergic Interneurons: Relevance to Neurodevelopment and Neurodevelopmental Disorders
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Impact of a Single Nucleotide Change or Non-Nucleoside Modifications in G-Rich Region on the Quadruplex–Duplex Hybrid Formation

Institute of Bioorganic Chemistry, Polish Academy of Sciences, 61-704 Poznan, Poland
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Biomolecules 2021, 11(8), 1236; https://doi.org/10.3390/biom11081236
Submission received: 19 July 2021 / Revised: 13 August 2021 / Accepted: 15 August 2021 / Published: 18 August 2021

Abstract

:
In this paper, a method to discriminate between two target RNA sequences that differ by one nucleotide only is presented. The method relies on the formation of alternative structures, i.e., quadruplex–duplex hybrid (QDH) and duplex with dangling ends (Dss), after hybridization of DNA or RNA G-rich oligonucleotides with target sequences containing 5′–GGGCUGG–3′ or 5′–GGGCGGG–3′ fragments. Using biophysical methods, we studied the effect of oligonucleotide types (DNA, RNA), non-nucleotide modifications (aliphatic linkers or abasic), and covalently attached G4 ligand on the ability of G-rich oligonucleotides to assemble a G-quadruplex motif. We demonstrated that all examined non-nucleotide modifications could mimic the external loops in the G-quadruplex domain of QDH structures without affecting their stability. Additionally, some modifications, in particular the presence of two abasic residues in the G-rich oligonucleotide, can induce the formation of non-canonical QDH instead of the Dss structure upon hybridization to a target sequence containing the GGGCUGG motif. Our results offer new insight into the sequential requirements for the formation of G-quadruplexes and provide important data on the effects of non-nucleotide modifications on G-quadruplex formation.

1. Introduction

In addition to the well-known double helix, nucleic acids can form a variety of different conformations. G-quadruplexes (G4) are structures that have received significant attention in recent years. The canonical, unimolecular G4 structures are formed from DNA or RNA G-quadruplex forming sequences (G4FS) comprising four G-tracts, and they fold into stacks of at least two guanine tetrads stabilized by Hoogsteen hydrogen bonds. Four G-tracts are necessary to form a G-quadruplex; however, they can also be part of two (bimolecular G4) or four separate strands (tetramolecular G4) [1]. The folding of G-quadruplexes into parallel, antiparallel, or hybrid-type topologies and their stabilities depends on the intramolecular factors, including sequence (number of G residues in G-tracts, the sequence of the loops, 5′/3′-flanking sequences) [2], environmental conditions (type of cations, pH, co-solvent) [3], and presence of small molecules (G4-stabilizing or G4-destabilizing ligands) [4] or proteins [5,6]. Polymorphism is inherent in G4 formation for most DNA G4-forming sequences. In contrast, RNA G-quadruplexes have long been considered as structurally monomorphic. However, RNA G4s with unusual motifs and arrangements of G-tetrads have been reported recently [7,8]. For example, in G4 cores of Spinach and Mango aptamers, some guanosine residues assume the syn conformation, consequently inducing the unusual antiparallel topology of G4 RNA [9,10,11]. Many modified nucleotides or non-nucleotide linkers have been used to manipulate G-quadruplex structures [12,13,14,15,16,17,18]. Rational incorporation of these modifications into G4FS can be used to stabilize particular conformers. Although the presence of G-tetrads is required to form G-quadruplexes, non-canonical tetrads have also been observed in high-resolution structures or models of G4s [8,19,20,21,22,23,24,25,26,27]. G4s are more tolerant to mutations than was previously thought. For example, it was found that G4 structures with G-tetrads replaced by non-guanosine tetrads can retain their biological functions [28]. Environmental conditions are one of the factors that can trigger the formation of unusual tetrads. For example, molecular crowding can destabilize the hairpin structure, thus promoting the formation of a G-quadruplex containing a GGUU tetrad [29]. It was also observed that intracellular mRNA can be entrapped in the formation of kinetically favored metastable hairpin-like structures that disturb the formation of the thermodynamically favored G-quadruplex [30]. At present, it is increasingly clear that the principles of predicting G4FS are more complicated than originally thought. Recently, the G4-seq technique, developed to detect G4 structures (next-generation sequencing, NGS), was used to identify over 700,000 DNA G4-forming sequences in the human genome [31]. This G4 mapping doubled the number of DNA G-quadruplexes predicted in silico by standard algorithms [32]. Using a similar rG4-seq approach, more than 13,000 G4FS were identified in 3000 human mRNAs [33]. In addition to mRNA, RNA G4 may also be present in mitochondrial RNA, tiRNA, piRNA, lncRNA, and miRNA [34,35,36,37,38]. Taking into account that these newly identified G4-forming sequences have the potential to fold into an unconventional class of G4 structures containing two-layer G-tetrads, bulges, or longer loops, the known repertoire of G4 structures may still be incomplete [31,33]. The most recent research has revealed that G-tracts involved in G-quadruplex formation can be located as far as 7–20 nucleotides away in the human genome [39]. The long loop sequences may adopt stem-loop secondary structures, and their formation was reported to accelerate DNA G4 folding [40]. G-quadruplexes containing stem-loop motifs are called quadruplex–duplex hybrid structures (QDHs). Structural studies highlighted a large diversity across their conformations [41,42,43,44,45,46,47]. The same G-rich sequence can adopt more than one hybrid form. For example, the coexistence of two different structures of DNA QDHs was observed for fragments of PIM1 or EGFR oncogenes [48,49]. Two separate oligonucleotides can assemble into bimolecular QDHs involving DNA-DNA, RNA-RNA, or DNA-RNA strands [50,51,52,53,54]. The endogenous QDH structures have been proposed as an attractive target for the regulation of oncogene expression in cells. When designing new small-molecule drugs, differences in QDH structure generating different potential sites for interaction with G4 stabilizing ligands should be taken into account [54,55,56]. For example, for quadruplex–duplex hybrid structures, an interesting approach based on the simultaneous binding of two ligands connected by a linker, PIP and cIKP, with different affinities for the duplex and G-quadruplex, has been proposed [55].
Previous reports suggested that most RNA G-quadruplexes exist in living cells in an unfolded state; however, recent studies proposed a model of a dynamic G4 RNA folding equilibrium controlled mainly by ions and G-quadruplex binding proteins [5,57,58,59]. The visualization of G4 RNA using the QUMA-1 ligand that binds only to the existing G4 structure confirmed a dynamic equilibrium between G4 RNA folded and unfolded states [60]. In general, the number of G4 RNA structures detected in cells depends on the technique used [33,60]. Nevertheless, detection of G4 RNA with G4-specific antibodies [61], G4 ligands [57,60], and a G-quadruplex-triggered fluorogenic hybridization probe [62,63] unequivocally confirmed their presence, both in vitro and in cells. Currently, G4 RNAs are important objects of research in biology, with a particular emphasis on the role of G4 in the flow of genetic information in cells [30] and G4-associated diseases [36,64,65,66,67,68,69]. Recently, we have demonstrated the sequence-specific targeting of a 56 nt long EGFR mRNA fragment comprising two distant GGGG tracts by RNA oligonucleotides composed of a chemically modified G-rich segment and a flanking 16 nt fragment complementary to mRNA EGFR [52]. As a consequence, the formation of a bimolecular RNA quadruplex–duplex hybrid structure containing a 28 nt long external loop adopting two duplex-stem structures was observed in vitro and in cells. We also noted the possibility of using RNA G-rich oligonucleotides conjugated to a fluorescent probe in a visualization of the density of the endogenous EGFR mRNA in MCF-7, HeLa, and A431 cancer cells.
Here, we demonstrate that target G-rich sequences that differ by a single nucleotide (5′–GGGCUGG–3′ vs. 5′–GGGCGGG–3′) can be targeted by G-rich oligonucleotides (QD) in a structure-specific manner. The change in a single nucleotide from uridine to guanosine results in the formation of a second 3 nt GGG tract in a target sequence and triggers the formation of a G-quadruplex motif (Figure 1). Using biophysical methods, we demonstrated the formation of two alternative structures, quadruplex–duplex hybrid (QDH) and duplex with dangling ends (Dss), depending on the target sequence (Figure 1). We undertook systematic studies on the effects of the type of oligonucleotide (DNA, RNA), non-nucleotide modifications (aliphatic linkers or abasic), and the covalently attached G4 ligand on the ability of G-rich oligonucleotides to recognize target molecules and to fold into QDH and Dss structures. Additionally, we proved that the replacement of a single nucleotide loop or junction with a non-nucleotide modification led to the formation of the QDH structures without affecting their stabilities. We also found that the presence of chemical modifications can induce the formation of the G4 motif on the UT target sequence despite the lack of one of the guanosine residues in the G-tract. The presented results may have implications for the structure-based design of G-rich antisense oligonucleotides and enable a more rational design of G-rich oligonucleotides in anticancer therapy.

2. Materials and Methods

All RNA and DNA sequences used in the study are shown in Table 1.

2.1. Synthesis and Purification of DNA and RNA Oligonucleotides

Synthesis of DMT-ON DNA and RNA oligonucleotides at the 1.0 µmole scale was performed according to routine procedures on a MerMade12 (BioAutomation, Irving, TX, USA) or a Gene World DNA synthesizer (K&A, Schaafhein, Germany) under the conditions recommended by the manufacturer. Samples were cleaved from the solid support, deprotected using standard procedures (Glen Research, Glen-Pak RNA, or DNA Cartridge Purification (DMT-ON), Sterling, VA, USA). The only exception was MMTr-C6-amino oligoribonucleotide (NH2-CCC), which was removed from the solid support and deprotected by overnight treatment with concentrated ammonia and ethanol (3:1) at 55 °C. Next, the solvent was evaporated to dryness, and the 2′-silyl protection was removed by treatment with 1.0 M triethylammonium fluoride at 65 °C for 2.5 h. In the next step, the precipitation was made by the addition of 5 mL of 1-butanol, and the samples were stored at −20 °C for 1 h. The precipitate was separated from the solution by spinning at 5000 rpm, 4 °C, for 10 min. The MMTr group was removed using 80% acetic acid/water for 30 min, and the solvent was evaporated to dryness. The NH2-CCC oligonucleotide was precipitated in the presence of 2% NaClO4/acetone. All oligonucleotides (Table 1) were then purified by RP-HPLC on a 1260 Infinity HPLC system (Agilent Technologies, Santa Clara, CA, USA) using XTERRA 5 µm, C18, 150 × 4.6 mm column with buffer A (0.1 M NH4HCO3/H2O) and buffer B (100% CH3CN) at a 1.5 mL/min flow rate, 70 °C. The buffer gradient was as follows: (1) 0–2 min 0% B; (2) 2–10 min 0–10% B; (3) 10–12 min 10–50% B; (4) 12–13 min 50–0% B; and (5) 13–30 min 0% B. UV detection was performed at λmax 268 nm.
All oligonucleotides were desalted using Amicon® Ultra 3K centrifugal filters (Merck, Millipore, Darmstadt, Germany) by loading on the filter, washing several times with 4000 µL MilliQ water, and then washing successively against ∼150 mM LiCl and against water.
To prevent dimerization of G-quadruplexes, experiments were performed in the presence of 50 mM KCl instead of 150 mM KCl [22,52].

2.2. NH2-CCC Labeling with o-BMVC-C3

o-BMVC-C3-NHS was coupled to the NH2-CCC according to a published procedure [52]. The purity and homogeneity of the oligo-o-BMVC-C3 conjugates were verified by 15% denaturing gel electrophoresis (Figure S1).

2.3. Denaturing Electrophoresis of RNA Oligo-o-BMVC-C3

First, 350 pmole RNA oligo-o-BMVC-C3 was suspended in 4 µL water and 4 µL of 8 M urea (in 8 µL final volume). Next, the sample was heated at 95 °C for 4 min, cooled to room temperature, and loaded on 15% TBE-Urea gels (Invitrogen, Thermo Fisher Scientific, Pittsburgh, PA, USA). Denaturing gel electrophoresis experiments were performed in 1.5 × TBE buffer. The electrophoresis experiment was run at 180 V for 2 h, and the gel was visualized by UV shadowing (Figure S1).

2.4. Non-Denaturing Electrophoresis of Oligonucleotides

First, 350 pmole of oligonucleotides were suspended in 50 mM KCl, 10 mM phosphate potassium buffer, pH 6.8 or 150 mM NaCl, 10 mM phosphate sodium buffer, pH 6.8 (in 6 µL final volume). Next, samples were heated at 95 °C for 3 min and gradually cooled to room temperature. Subsequently, samples were mixed with 2 µL of 50% glycerol and loaded on 20% TBE gels (Invitrogen, Thermo Fisher Scientific, Pittsburgh, PA, USA). O’RangeRuler 5 bp and DNA ladder containing 50 and 100 bp as brighter bands (Thermo Fisher Scientific, Pittsburgh, PA, USA) were used as molecular markers. Gel electrophoresis experiments performed in a 0.5 × TBE buffer were run at 160 V for ~2 h at 4 °C (an ice bath), and the gels were first visualized by UV shadow (image) and then stained by N-methyl mesoporphyrin IX (NMM; Frontier Scientific, Newark, DE, USA). After post-staining the gels with NMM, the gels were scanned with a Fuji FLA-5100 imaging system (Fujifilm Life Sciences, Cambridge, MA, USA).
The last gel was viewed by UV shadow to visualize all RNA (image) and exposed to 473 nm light to visualize fluorescently labeled G-quadruplex (o-BMVC-C3-CCC/GT hybrid).

2.5. NMR Experiments

NMR experiments were performed on a 700 MHz Bruker AVANCE III spectrometer (Bruker Corporation, Billerica, MA, USA) equipped with a QCI CryoProbe. The oligonucleotides were annealed by heating to 90 °C and then slowly cooled to room temperature. The 1H NMR spectra of oligonucleotides were recorded at 25 °C or in the 25–75 °C range in 3 mm thin wall tubes with a sample volume of 200 μL. The final concentration of QD/GT, QD/UT was ~0.1 mM in a buffer containing 50 mM KCl, 10 mM phosphate potassium, or 150 mM NaCl, 10 mM phosphate sodium at pH 6.8. A mixture of 90% H2O and 10% D2O was used for experiments undertaken to study imino protons. The water signal was suppressed by excitation sculpting with a gradient pulse. Spectra were processed and prepared with TopSpin 3.2 Bruker Software.

2.6. UV Thermal Denaturation Curves

First, 590 pmole of QD/GT and QD/UT oligonucleotides were dissolved in a buffer containing 50 mM KCl, 10 mM phosphate potassium buffer, or 150 mM NaCl, 10 mM phosphate sodium buffer at pH 6.8. Thermal denaturation curves were obtained by monitoring at 260 and 295 nm with a JASCO V-650 spectrophotometer (JASCO International Co., Ltd., Tokyo, Japan) using quartz optical cuvettes of 0.5 path length with the sample volume of 150 µL. Samples were protected against evaporation by silicone oil. Before measurements, the cuvettes filled with samples were spun at 5000 revolutions per minute for 3 min to avoid the formation of air bubbles during measurements. The temperature range was 20–90 °C, using a scan rate of 0.5 °C min−1. Spectra were processed and prepared using Origin 8 Software (OriginLab Corporation, Northampton, MA, USA). The melting temperature was determined by the local maximum of the first derivatives of the absorbance vs. temperature curve.

2.7. Fluorescence Measurements

One equivalent of NMM was mixed with QD/GT or QD/UT (25 µmol/L) in a buffer containing 10 mM potassium phosphate, 50 mM KCl, 6.8 pH at 25 °C. Fluorescence spectra were measured from 550 to 750 nm with a 2 nm step, 800 V detector sensitivity using 400 nm excitation. Fluorescence spectra were recorded using a JASCO J-815 CD/fluorescence spectropolarimeter (JASCO International Co., Ltd., Tokyo, Japan).

2.8. CD Measurements

CD spectra of oligonucleotides were recorded using a JASCO J-815 spectropolarimeter (JASCO International Co., Ltd., Tokyo, Japan) equipped with a temperature controller. For each sample, 3 spectral scans were accumulated at 25 °C over wavelengths in the range from 220 to 320 nm. Samples of each oligonucleotide were prepared at a concentration of 8.0 μM in 50 mM KCl, 10 mM potassium phosphate buffer, or 150 mM NaCl, 10 mM sodium phosphate buffer at 6.8 pH using a 0.5 cm path length quartz cuvette with a volume of 1750 µL. Next, samples were prepared by heating the oligonucleotides at 90 °C for 5 min and gradually cooling to room temperature. CD spectra were expressed in the units of molar ellipticity Δɛ (cm2 mmol−1), without normalization by the number of residues in the oligonucleotide. Spectra were processed and prepared using the Origin 8 Software.

3. Results

3.1. Formation of Alternative Structures by Hybridization of G-Rich Oligonucleotide to Target Sequences That Differ by a Single Nucleotide Residue

The proposed method to discriminate between two target RNA sequences that differ by one nucleotide only is schematically presented in Figure 1. Generally, it relies on the formation of alternative structures, i.e., a quadruplex–duplex hybrid (QDH) or duplex with long dangling ends (Dss) (Figure 1D,E), after the binding of G-rich oligonucleotides (Figure 1A) to target sequences that differ by a single U/G change (Figure 1B,C). We designed two model target sequences (Table 1), GT and UT, containing 5′–GGGCGGG–3′ or 5′–GGGCUGG–3′ fragments, respectively. The CCC oligonucleotide (Figure 1A, Table 1) is composed of two segments, duplex forming (D) and G-quadruplex forming (Q). The duplex forming segment consists of 14 nt and is complementary to 14 nt fragments of both targets (GT and UT, Figure 1B,C) being responsible for the formation of a duplex domain (Figure 1D,E). The G-rich segment of the CCC strand was designed to assemble into a G-quadruplex motif only with the GT target (CCC/GT, Figure 1D) but not with UT (CCC/UT, Figure 1E).

3.2. Evaluation of Secondary Structures of CCC/GT and CCC/UT RNA:RNA Complexes

Both target sequences, GT and UT, were hybridized to CCC RNA G-rich oligonucleotides (Table 1) and characterized by several experimental methods. To determine the type of structure formed in the solution, we used 1H NMR spectroscopy. In general, imino proton signals associated with Watson–Crick base pairs typically appear at 12–15 ppm, and those between 10.5 and 12 ppm are characteristic of guanosine imino protons involved in G-tetrad formation [70]. The imino region of 1H NMR spectra of CCC/GT and CCC/UT are shown in Figure 2A and compared with 1H NMR spectra of individual strands of the CCC/GT and CCC/UT molecules (UT, GT, and CCC), the CGGGCGGGC G-quadruplex (2Q), and the 14 bp duplex (DX). In the presence of K+ ions, the NMR spectrum of CCC/GT displayed imino proton signals in both duplex and G-quadruplex regions, and the peak distribution patterns resembled that of the individual components, duplex, and G-quadruplex. This was consistent with the formation of the quadruplex–duplex hybrid structure (QDH) by CCC/GT. When CCC RNA oligonucleotide was hybridized to the UT target, only signals characteristic of the imino protons involved in the formation of Watson–Crick base pairs were observed. This suggests that the replacement of one of the guanosine residues by uridine, such as the UT target, can prevent G-quadruplex formation (Figure 1C,E and Figure 2A). For the CCC/UT complex, only signals corresponding to the duplex domain were observed while G-rich fragments remained unstructured (Figure 1E and Figure 2A), indicating the formation of a duplex structure with long dangling ends (Dss).
Next, we performed native polyacrylamide gel electrophoresis (PAGE) in the presence of 50 mM KCl to assess the molecularity of CCC/GT and CCC/UT and to verify structural differences between these two complexes (Figure 3). The bands were first visualized by UV shadowing at 254 nm and immediately after electrophoresis were stained with NMM, a dye that exhibits significantly increased fluorescence only upon binding to parallel G-quadruplexes [71,72]. The first lane in Figure 3 corresponds to DNA markers (10–100 bp), and lanes 2–6 correspond to CCC, GT, UT, CCC/GT, and CCC/UT, respectively. The mobilities of the bands corresponding to oligonucleotides CCC, GT, and UT (lanes 2–4) were in agreement with the formation of the hairpin structures as predicted by the RNAstructure software (Figure 2B). As expected, the main bands in lanes 5 and 6 migrated similarly to that of the 50 bp DNA marker, indicating the formation of bimolecular structures. In lane 6, bands of non-associated strands of CCC and UT were observed. After staining with NMM, only the band corresponding to CCC/GT (Figure 3, lane 5′) was clearly visible, thus supporting the formation of the G-quadruplex domain.
The binding of NMM to the G-quadruplex domain of CCC/GT suggested its parallel topology [71,72], which was further confirmed by the presence of a characteristic positive band at ~265 nm and a negative band at 240 nm in the CD spectrum (Figure S2). The CD spectrum of CCC/UT differed from that of CCC/GT only in the intensity of the band at ~265 nm. Unfortunately, CD spectra of parallel G-quadruplexes are very similar to those of A-form duplex structures. As a result, they have little use in structural analysis of RNA quadruplex–duplex hybrids [73].
The thermal stabilities of CCC/GT and CCC/UT structures were determined by analyzing UV-melting profiles at 260 and 295 nm in the presence of potassium cations (Table 2, Figure S3). The melting temperature of CCC/GT determined at 260, and 295 nm was 69.4 °C and 68.8 °C, respectively (Table 2, Figure S3). The profile of the melting curve at 295 nm did not reflect the typical reverse sigmoid. This resulted from the overlapping of two opposing effects, hypochromic and hyperchromic, for the melting of G-quadruplex and duplex motifs at this wavelength (Figure S3). A similar profile of the melting curve at 295 nm was observed previously for the RNA hairpin structure in equilibrium with the G-quadruplex, but this phenomenon was not discussed [74]. The melting temperature of the CCC/UT duplex determined at 260 nm was 61.0 °C. The higher stability of the CCC/GT hybrid relative to the CCC/UT duplex (ΔTm = 8.4 °C) confirmed the G4-dependent stabilization effect.

3.3. Evaluation of Secondary Structures of DNA-CCC/GT and DNA-CCC/UT Complexes

Next, we compared the tendency of G-rich RNA and DNA molecules to form QDH and Dss structures. We hybridized the DNA equivalent of the CCC sequence, DNA-CCC (Table 1), with GT and UT RNA targets in the presence of potassium cations.
Figure 4A shows the 1H NMR spectra of two DNA:RNA hetero-complexes: DNA-CCC/UT and DNA-CCC/GT. The spectrum of the DNA-CCC/GT displayed signals both in duplex (11.5–15 ppm) and G-quadruplex (10.5–12 ppm) regions, indicating the formation of a hybrid quadruplex–duplex structure. The 1H NMR spectrum of DNA-CCC/UT showed that signals in the region of 10.5–12 ppm characteristic of G-quadruplex formation were missing, thus suggesting the formation of a Dss structure. However, comparison of the 1H NMR spectrum of DNA-CCC/UT with DNA-CCC and UT clearly showed that it was a sum of signals characteristic of two separate strands, i.e., DNA-CCC and UT, thus suggesting that the structure of a DNA: RNA heteroduplex with dangling ends was not formed (Figure 4A,B).
Figure 5 compares the electrophoretic mobility of hetero-complexes DNA-CCC/GT, DNA-CCC/UT (lanes 4–5), and single-strand DNA-CCC (lane 6) relative to RNA homo-complexes CCC/GT and CCC/UT (lanes 2–3). Bands corresponding to bimolecular structures DNA-CCC/GT and DNA-CCC/UT migrated faster than their RNA:RNA counterparts. The smaller intensity of the band corresponding to DNA-CCC/UT (lane 5) than that of CCC/UT (lane 3), and the presence of unhybridized single-strands DNA-CCC and UT in lane 5, were indicative of a poor hybridization of DNA to RNA strands. When the gel was post-stained with NMM, fluorescence bands were observed only for CCC/GT and DNA-CCC/GT and indicated the formation of the G-quadruplex domain in both cases (Figure 5, bands 2ʹ, 4ʹ). Additionally, the presence of a minor band characterized by slow electrophoretic mobility for CCC/GT indicated the formation of higher-order structures (HOS) (Figure 5, lane 2 and lane 2′). An increase in intensity with time was observed for those retarded bands for CCC/GT (Figure S4). This structure could be a dimer of two G-quadruplex domains of QDH structures. The parallel topology of the G4 domain of the DNA-CCC/GT hybrid was additionally confirmed by a characteristic CD spectrum (Figure S5). The DNA-CCC/GT structure was less stable than that of CCC/GT, as indicated by the comparison of melting temperatures, which were 61.4 and 69.4 °C, respectively (Figure S6, Table 2). Data obtained from the NMR spectra showed that the DNA-CCC/UT complex was not formed at room temperature (Figure 4). Therefore, for further studies, we decided to use RNA G-rich oligonucleotides because they are definitely more stable when hybridized to RNA targets, thus enabling the determination of secondary structures. In addition, the sequential requirements for the formation of RNA quadruplex–duplex hybrids are less known.

3.4. Impact of the Chemical Modifications on the Ability to Form Alternate QDH and Dss Structures

The modified G-rich oligoribonucleotides used in this study, and their corresponding names and notations used to describe the bimolecular complexes, are summarized in Table 1. Bolded cytidine residues in the CCC sequence point to the modification sites (Table 1, Figure 6). Non-nucleotide modifications, such as abasic (a), aliphatic 1,2-ethanediol (L2), 1,3-propanediol (L3), and 1,4-butanediol (L4) linkers, were chosen to study their effect on the formation of the G-quadruplex domain. It was previously reported that the presence of abasic residues and long aliphatic linkers in the loop region promoted and stabilized the formation of parallel G-quadruplex topologies [16,17,75]. In turn, the presence of a phosphate group (p) or additional nucleotide residues at the 5′-end should prevent dimerization of G-quadruplexes [76,77]. A 2′O-Me modification in a duplex segment was used to estimate its potential to stabilize QDH and Dss structures.
To determine the impact of the chemical modifications on the formation and stability of G-quadruplex motifs, we compared the 1H NMR spectra of unmodified and modified G-rich oligoribonucleotides hybridized to GT and UT targets in the presence of K+ and Na+ cations. The corresponding 1H NMR spectra of QD/GT complexes recorded in buffers containing K+ or Na+ cations are presented in Figure 7 and Figure S6, respectively. All modified G-rich oligoribonucleotides after hybridization to GT target display the imino peak patterns typical of the structure containing Watson–Crick and Hoogsteen-type hydrogen bonds. In general, all these spectra were similar to that of the model CCC/GT hybrid, regardless of the type of modification, their modification site, and the type of ions present in the solution, K+ or Na+. For two complexes in a solution containing potassium ions, AaC/GT and Aaa/GT, a slightly different pattern was observed in the Hoogsteen region, which could be due to the presence of adenosine residues at the 5′-end in AaC and Aaa strands (Figure 7).
When modified G-rich oligoribonucleotides were hybridized to the UT target (QD/UT complexes) in the presence of K+, in addition to the imino protons from Watson–Crick base pairs, broad resonances were observed in the 10.5–12.0 ppm region for several complexes containing modified residues (for instance: paa/UT, Aaa/UT, pL4C/UT, AaC/UT, or pCa/UT) (Figure 7). In sodium phosphate buffer, a broad resonance suggesting the presence of a G4 motif was clearly visible only for molecules containing two abasic residues (paa/UT, Aaa/UT), whereas for AaC/UT these G4 imino signals were very broad (Figure S7). These results suggest that the incorporation of chemical modifications into G-rich strands can promote the formation of structures other than Dss after hybridization to the UT strand. These structures may contain G-quadruplex motifs instead of dangling ends.
The effect of the chemical modification on the formation of QDH and Dss structures was further analyzed by native PAGE (Figure 8, UV and NMM visualization). In the presence of K+ cations, the migration rate of bands corresponding to modified QD/GT complexes (Figure 8A, lines 3–20) was in agreement with the formation of the CCC/GT quadruplex–duplex hybrid (Figure 8A, lines 1) and these bands were efficiently stained with NMM (Figure 8B, lines 1′). Slower migration of pL2C/GT, pL3C/GT, and pL4C/GT bands probably reflected the reduced negative charge of these hybrids resulting from the introduction of aliphatic linkers. Moreover, for QD/GT complexes comprising two abasic residues or L3 and L4 aliphatic linkers and a phosphate group at the 5′-end (paa/GT, pL3C/GT, pL4C/GT, Figure 8B, lines: 7′, 13′, 15′ respectively), the main bands were split into two bands with similar fluorescence intensity after staining with NMM. The presence of two bands suggested the formation of two different conformers. For this parallel G-quadruplex motif, there are two forms possible that differ in the directionality of the hydrogen bonds within the G-tetrads (Figure S8). Additionally, similarly to CCC/GT molecule, the presence of minor bands characterized by slow electrophoretic mobility for most modified QD/GT complexes indicated the formation of higher-order structures (HOS) (Figure 8B).
A comparison of the migration of QD/GT hybrids containing a phosphate group or adenosine at the 5′-end of QD (Figure 8B, band 3′ vs. 5′ and 7′ vs. 9′) showed that the presence of adenosine prevented the dimerization of two G-quadruplex units of the QDH structure more than the phosphate group. The QD/UT complexes containing non-nucleotide linkers migrated with a mobility similar to that of the model CCC/UT (Figure 8). Most of the bands corresponding to bimolecular structures did not stain with NMM (even lines), confirming the formation of duplexes with dangling ends as the dominant form. However, for paa/UT (lane 8′), Aaa/UT (lane 10′) low-intensity fluorescence bands were observed after staining with NMM (Figure 8B). To confirm the formation of a G-quadruplex motif in paa/UT and Aaa/UT, we performed an experiment in which we compared fluorescence emissions after adding equimolar ratios of representative QD/GT and QD/UT to a solution of NMM (Figure S9). The intensity of the fluorescence signal of NMM and NMM after the addition of CCC/UT (1:1) was similarly low, indicating a weak affinity of NMM towards CCC/UT. For NMM complexes with Aaa/GT and paa/GT, we observed a moderate 7- to 9-fold enhancement of NMM fluorescence emission, whereas, for Aaa/UT and paa/UT, a 4–6.5-fold enhancement of NMM fluorescence emission was determined. The intensities of the fluorescence signals of the NMM in complex with Aaa/UT and paa/UT were significantly higher compared to those obtained for NMM in complex with CCC/UT duplex. These results confirmed the selective binding of NMM to the non-canonical G-quadruplex domain of Aaa/UT and paa/UT.
Additionally, we recorded 1H NMR spectra after the addition of one equivalent of NMM to NMR samples containing paa/GT, paa/UT, or Aaa/UT (Figure S10). Ligands binding to paa/GT resulted in a significant broadening of the imino signals corresponding to G-tetrads and the appearance of signals from the NMM ligand in the 10.3–9.0 ppm region (Figure S10a). NMM-dependent broadening of the imino signals was also observed in the 1H NMR spectra of Aaa/UT and paa/UT (Figure S10b,c). These results strongly supported data obtained from native PAGE and fluorescence experiments, indicating that a noncanonical G-quadruplex domain could form in the Aaa/UT and paa/UT molecules. The formation of G-quadruplex motifs in Aaa/UT and paa/UT structures, which was observed in the presence of K+ ions, was also observed in native PAGE performed in the presence of Na+ ions (Figure S11, lines 5′ vs. 6′ and 7′ vs. 8′).
To further characterize QD/GT and QD/UT complexes, we analyzed their UV-melting profiles in solutions containing K+ or Na+ ions. The apparent melting temperatures determined based on measurements at 260 nm are collected in Table 2 and Table S1 and presented graphically in Figures S12a,b and S13a,b (melting curves at 260 and 295 nm). For several quadruplex–duplex hybrid structures, accurate determination of Tm at 295 nm was not possible (Figures S12 and S13). However, analysis of melting profiles at 295 nm was used to estimate the stability of the G-quadruplex motif and to confirm the formation of QDH or Dss structures in both K+ and Na+ solutions. Profiles of melting curves at 295 nm indicated the presence of G-quadruplex motifs for all QD/GT hybrid molecules and for paa/UT, Aaa/UT, and AaC/UT in potassium phosphate buffer (Figures S12a and S13a). All QD/GT structures had comparable Tm values in the range of 67.4–69.0 °C, regardless of the type and site of modification, and were very close to that of CCC/GT (69.4 °C) (Table 2). The only exception was observed for 2’OMe modification, which is well known to stabilize RNA duplex structures (Table 2). In addition, we confirmed that the introduction of non-nucleotide modifications does not affect the kinetics of QDH structure formation (no hysteresis was observed, K+, data not shown). In general, in potassium phosphate buffer, all QD/GT structures were more stable than the corresponding QD/UT complexes; their Tm values were higher at least by 6.0 °C. For paa/UT, Aaa/UT, and AaC/UT, molecules for which the presence of noncanonical G-quadruplex motifs was suggested, analysis of the melting curves at 295 nm indicated that these motifs were less stable than their QD/GT counterparts (in average by 7.25 °C). Apparent Tm values determined for the quadruplex–duplex hybrid structures in the presence of sodium cations (Table S1) were comparable to the corresponding melting temperatures of duplexes with dangling ends. This was due to the increased stability of the duplex motif ascribed to the different ionic strengths of the buffers (60 mM K+ versus 160 mM Na+). By comparison, the exchange of K+ with Na+ caused a significant decrease in the stability of the G-quadruplex domain of QD/GT, which could be estimated from the melting curve recorded at 295 nm and in 1H NMR spectra (Tables S1, S12b, and S14b).
Different thermal stabilities of duplex and G-quadruplex motifs depending on the solution conditions were best illustrated by 1H NMR spectra recorded as a function of temperature. Figure S14 shows the 1H NMR spectra of CCC/GT, paa/GT, and pL3C/GT obtained in the range of 25–75 °C in buffers containing 60 mM K+ and 160 mM Na+ (Figure S14a,b). In the potassium environment, with increasing temperature, Watson–Crick imino signals disappeared faster than those from G-tetrads. Moreover, the stability of the G-quadruplex motifs appeared to be higher for molecules containing modified residues. By comparison, in 1H NMR spectra recorded in the presence of sodium cations, the G4 signals tend to disappear at a faster rate compared to Watson–Crick base pairs.
Although CD spectra of RNA G-quadruplexes were very similar to those of A-form duplexes, the presence of a clear positive peak at 260 nm and a negative peak at 240 nm for all QD/GT complexes studied (Figure S15, Na+, K+) verified the parallel topology of their G-quadruplex motifs. For QDH structures containing chemically modified G-rich strands, an increase in the intensities of the Cotton effects at 265 nm was observed compared to CCC/GT. As expected, the shape of the CD spectrum of the CCC/UT duplex with dangling ends was similar to that of the quadruplex–duplex hybrids; however, the bands at 265 nm were of lower intensity than the spectra of the corresponding QDH structures (Figures S2 and S15). For paa/UT and Aaa/UT molecules, the observation of increased intensity of the ellipticity values in CD spectra was consistent with the presence of the G-quadruplex motif in their structures (Figure S16).

3.5. Influence of CCC Oligonucleotide with Covalently Attached G4 Ligand on the Formation and Stability of QDH and Dss Structures

Several small-molecule ligands are known to specifically bind and stabilize G-quadruplexes and are considered to be promising therapeutic targets. However, few ligands are known that are specific to a given G-quadruplex topology [78]. Therefore, we decided to covalently attach a well-known o-BMVC ligand to the CCC oligoribonucleotide and investigate its impact on the secondary structure formation after hybridization to the target sequence GT or UT. We chose the o-BMVC ligand due to its properties, which are ideal for the purpose of our study, i.e., a large difference in fluorescent intensity after binding to G4 structures or duplexes [79,80]. Derivative of o-BMVC, propionic acid NHS ester (o-BMVC-C3-NHS) was made in-house [52,81] and then attached to the 5′-end of the CCC molecule via a C6-aminolinker (o-BMVC-C3-CCC) (Figure S17).
The structure and stability of complexes of o-BMVC-C3-CCC (L-CCC, Table 1) annealed to GT or UT (Figure 9) were investigated using 1H NMR, UV, CD, and fluorescence spectroscopy. 1H NMR spectra of CCC/GT, o-BMVC-C3-CCC/GT, and o-BMVC-C3-CCC/UT are shown in Figure 10A. The presence of imino signals in the o-BMVC-C3-CCC/GT spectrum in the region typical of Watson–Crick base pairs and of G-tetrads confirmed the formation of both duplex and G-quadruplex domains (Figure 10A). Furthermore, a significant broadening of signals only in the 10.5–12 ppm region indicated that ligand bounds selectively to the G-quadruplex motif. In the 1H NMR spectrum of the o-BMVC-C3-CCC/UT, signals corresponding to the duplex domain remained sharp, but a broad signal of low intensity additionally appeared in the 10.5–12 ppm region. This suggested that o-BMVC-C3 could possibly interact nonspecifically with G-rich dangling ends.
Structures of o-BMVC-C3-CCC/GT and o-BMVC-C3-CCC/UT were further investigated using native PAGE gel electrophoresis. The migration of o-BMVC-C3-CCC/GT (Figure 10B, lane 2) and o-BMVC-C3-CCC/UT (Figure 10B, lane 3) visualized in UV light appeared to be undoubtedly slower than that of the o-BMVC-C3-CCC single-strand (Figure 10B, lane 1), clearly showing the presence of bimolecular complexes. After exposure of the gel to 473 nm light, strong fluorescence was observed only for o-BMVC-C3-CCC/GT when visualized at 532 nm. The large difference in fluorescence intensity observed between o-BMVC-C3-CCC/GT and o-BMVC-C3-CCC/UT molecules confirmed the formation of a G-quadruplex domain only in the o-BMVC-C3-CCC/GT hybrid (Figure 10B). The presence of a weak fluorescent band corresponding to the o-BMVC-C3-CCC/UT conjugate supported the earlier observation from analysis of the 1H NMR spectra, suggesting that the covalently attached ligand o-BMVC-C3 could interact nonspecifically with G-rich dangling ends. Additionally, the formation of higher-order species was observed for both complexes. To determine the impact of the covalently attached ligand on the thermal stabilities of the studied molecules, we performed UV-melting experiments in the presence of potassium cations. The melting temperature obtained for o-BMVC-C3-CCC/GT (Tm = 72.1 °C) was 2.7 °C higher than that for CCC/GT (Figure S18). In turn, attachment of the G4 ligand to CCC/UT did not affect its thermal stability; Tm values of o-BMVC-C3-CCC/UT and CCC/UT were the same within experimental error, 61.2 °C and 61.0 °C, respectively.
The CD spectra of o-BMVC-C3-CCC/GT and o-BMVC-C3-CCC/UT are compared to the CD spectrum of CCC/GT in Figure S19. The presence of a strong maximum at 265 nm indicated that the G-quadruplex domain of o-BMVC-C3-CCC/GT maintained a parallel topology. As can be seen, the intensity of this band was significantly higher than that of o-BMVC-C3-CCC/UT, which further supported the formation of a stable G-quadruplex domain. By comparison, the lower intensity of this band observed for CCC/GT relative to o-BMVC-C3-CCC/GT suggested that in a molecule without an attached ligand, the G-quadruplex domain could be less structured.

4. Discussion

4.1. Quadruplex–Duplex Hybrid Structures (RNA vs. DNA)

Thousands of RNA and DNA sequences that are able to form stem-loop-containing G-quadruplexes have been identified in the human genome and transcriptome using bioinformatic methods [39]. Systematic studies of the unimolecular DNA QDH structures comprising a hairpin or multiple stems in place of at least one of the loops suggest that the junction between the G-quadruplex and duplex motifs is characteristic of a particular G4 topology: parallel, antiparallel, or hybrid [41,44]. Recently, QDH structures with a duplex stem incorporated into a bulge of a G-quadruplex or structures with a duplex stem build from the G4 flanking regions have also been demonstrated [45,46,47,82]. QDH structures identified in RNA aptamers (Spinach, Mango, and FMRP-binding) were found to have unique G-core architectures [82]. By comparison, all known bimolecular RNA:RNA or DNA:RNA quadruplex–duplex structures tend to be structurally conservative and are constructed of a parallel G-quadruplex domain and a duplex stem located in the place of an external loop [50,51,52]. This is in agreement with our results showing that all RNA-RNA homohybrids studied in this work and DNA-CCC/GT heterohybrids adopt this architecture (Figure 1D, Figure 2B and Figure 4B). This structural conservatism of bimolecular RNA:RNA or DNA:RNA QDHs allows stable hybrid structures to be predicted, and this could help in the future to design drugs that selectively stabilize G-quadruplex motifs.

4.2. Stability of QDH and Dss Structures

Although a type of G-rich oligonucleotide (DNA or RNA) did not change the structural preferences of QD/GT complexes, the use of DNA strands instead of RNA reduced the stability of DNA-CCC/GT heterohybrid by 8 °C. This destabilizing effect was even more pronounced in the case of the DNA-CCC/UT complex. We showed that this complex was unstable and did not form at room temperature. In the presence of potassium cations, the thermal stability of the examined quadruplex–duplex hybrid structures increased by ~8 °C compared to the corresponding duplexes with dangling ends. Due to the significant destabilization of the G-quadruplex domain (Figure S14b) in a solution containing Na+ ions, this duplex domain was responsible for the stability of the QDH. The other important factor that could influence the stability and/or structure of G-quadruplexes was the incorporation of non-nucleoside modifications, such as an aliphatic linker or abasic residues, into a G-rich region. It was previously shown that the replacement of TTA loops with aliphatic linkers of different lengths in the human telomeric repeat sequences (GGGTTA)3G3 resulted in a conversion of G-quadruplex topology from antiparallel to parallel [16]. The stability of the modified G-quadruplexes was strongly correlated with the length of the linker in the order 1,8-octanediol > hexaethylene glycol > 1,3-propanediol. Interestingly, the G-quadruplex bearing the 1,8-octanediol loop was found to adopt several stable conformations, including the formation of the two-layer G-quadruplex structure as proposed by the authors. Replacement of one nucleoside loop with abasic residue resulted in the most stable DNA G-quadruplexes. Furthermore, it was shown that RNA G-quadruplexes containing 1 nt external loops are the most stable [83]. In the present study, the replacement of a single nucleoside loop with a non-nucleotide modification led to the formation of stable QDH structures irrespective of the modification type or the length of the linker. When the external loop was replaced by L3 or L4 aliphatic linkers, the formation of two bands corresponding to two conformers of pL3C/GT and pL4C/GT quadruplex–duplex hybrids were observed in the native PAGE experiment. The formation of two conformers was also detected for the paa/GT complex containing abasic residues located in loop and junction positions. It has been previously reported that base composition proximal to the junction between duplex and G-quadruplex motifs can also play an important role in both the structure and stability of the hybrids. For example, when the terminal C:G base pair was substituted with a less stable T:A base pair in a junction, two bands appeared in gel electrophoresis, probably due to two different conformations present in solution [84]. In general, non-nucleotide linkers can contribute to the increased flexibility of the loop or influence the interaction between the loop and G-tetrads. Our results indicate that the replacement of the loop or junction residues with a flexible linker did not change the structure and stability of QDHs. Furthermore, flexible linkers appeared to have a high propensity to form external loops. When the G4-selective ligand was conjugated to the G-rich oligonucleotide, a moderate increase in stability was observed for o-BMVC-C3-CCC/GT relative to CCC/GT. In contrast, the presence of the non-nucleotide modifications or G4 ligand did not affect the stability of duplexes with dangling ends.

4.3. Formation of Non-Canonical Hybrid Structures

Here, we demonstrated that one nucleotide change in a G-rich oligonucleotide (change from 5′–GGGCGGG–3′ to 5′–GGGCUGG–3′) prevented the formation of the G-quadruplex domain in the CCC/UT complex. For the CCC/UT molecule, we considered the possibility of the formation of structures other than Dss. For this purpose, we used the RNAstructure software, which suggested several stable secondary structures with additional canonical Watson-Crick or non-canonical base pairs (Figure S20). However, analysis of the imino region of the 1H NMR spectrum (Figure 2) indicated that the duplex with the dangling ends was a dominant form of CCC/UT. The number of imino signals observed in the 1H NMR spectrum precluded the formation of additional stable base pairs in the structure of CCC/UT. However, in the 1H NMR spectra of paa/UT and Aaa/UT, we surprisingly noticed an appearance of signals indicating the formation of G-tetrads, even in the solution containing Na+ cations. The formation of this G-quadruplex motif was also confirmed by the observation of low-intensity fluorescence bands in native PAGE after staining with NMM. UV-melting experiments showed that thermal stabilities of paa/UT and Aaa/UT were lower than those determined for their counterparts paa/GT and Aaa/GT. These results suggest that G-quadruplex domains in paa/UT and Aaa/UT could be composed of two G-tetrads, or two G-tetrads and one mixed GGGU tetrad (Figure S21). The binding of NMM to this atypical G-quadruplex motif was additionally confirmed by 1H NMR experiments and an increased emission in fluorescence spectra after the addition of paa/UT or Aaa/UT to a solution containing NMM. It was previously reported that the hairpin in dilute solution could undergo the transition to a G-quadruplex containing one mixed GGUU tetrad. This type of G-quadruplex exhibited reduced thermal stability and fluorescence intensity after NMM binding [29]. In summary, our data suggest that the formation of G-quadruplexes is possible even when the sequences lack regular G-tracts, which were not predicted by algorithms developed to search for putative G-quadruplex forming sequences. Moreover, the obtained results advance our knowledge regarding the influence of non-nucleotide modifications on the formation of non-canonical G-quadruplexes. In this respect, a recent report should be mentioned that demonstrated that the presence of non-guanosine residues or bulges in the G-core generally destabilized the G-quadruplex structure but nevertheless did not prevent G-quadruplex formation [33,85,86]. Another study has shown that G-quadruplexes with mixed central tetrads were able to compensate for the lack of a G-tetrad in the context of both GTP-binding and peroxidase activity [28]. These data shed light on the possible biological functionality of non-canonical G-quadruplexes.

4.4. Application Potential of G-Rich Oligonucleotides

We have recently shown that the formation of a stable G4 motif on an mRNA template was an effective steric hindrance for ribosomes and results in stalling and inhibition of protein synthesis. The silencing activity of oligoribonucleotides modified with 1,6-hexanediol and an abasic linker toward the EGFR mRNA target appeared to be at least 20% higher than for classical 16 nt antisense oligonucleotides [52]. Several studies indicated that a correlation exists between G-quadruplex stability and efficacy of gene expression silencing [50,87,88]. Our data revealed that RNA:RNA QDH was more stable, by 8 °C, than the corresponding DNA:RNA heterohybrid. Additionally, we demonstrated that abasic residues and L2-L4 aliphatic linkers could perfectly mimic external loops. A previous study indicated that the incorporation of the methylene unit (L1) into the TBA sequence could enhance the lifetime of the G-quadruplex in blood [89]. Another study has shown that the presence of a long aliphatic linker between G-tracts can prevent the hybridization of a G-rich oligonucleotide to a C-rich sequence [16]. We strongly believe that chemically modified RNA G-rich antisense oligonucleotides forming stable quadruplex–duplex hybrid structures with mRNA have significant application potential [52]. The ability of these G-rich oligonucleotides to recognize a single nucleotide change in the target sequences in a secondary structure-dependent manner (Figure 1) increases their attractiveness. For therapeutic purposes, the advantages of using G-rich oligonucleotides with a duplex forming segment that is complementary to the target sequence rely on the ability to guide these oligonucleotides to a specific mRNA site. Reviewing the group of representative human mRNAs, we found that over 90% of mRNA annotated in the human genome (RefSeq) contain at least one potential target sequence composed of two neighboring blocks of consecutive guanosine residues separated by one or two non-guanosine residues: HGGGH(H)GGGH, HGGGGH(H)GGGGH (H denotes A, C or T). With the use of mRNA sequences retrieved from the RefSeq database (www.ncbi.nlm.nih.gov/refseq) and the IGV tool (www.broadinstitute.org/igv (accessed on 31th October 2019)), we identified potential target sequences in protooncogenes EGFR, KRAS, HER2, and genes implicated in human Mendelian dominant disorders, such as PTPN11 and FGF23. The G-rich oligonucleotides may also be considered for allele-specific targeting. For example, the well-recognized cancer-driving EGFR mutation, L858R (c.2573T>G), introduces a 5′–GGGCUGG–3′ to 5′–GGGCGGG–3′ change in the sequence of the EGFR mRNA, creating a potential target for G-rich oligonucleotides [90]. L858R next to exon 19 in-frame deletions are the most frequent EGFR mutations (cumulatively accounting for ~90% of all EGFR mutations), well recognized as biomarkers of targeted therapy of non-small cell lung cancer and other cancers with the use of EGFR-specific tyrosine kinase inhibitors.

4.5. G4 Ligand Conjugated to Oligonucleotide

Ligands interacting with G4s have attracted significant attention as potential anticancer therapeutics [91,92]. In view of the high number of G-quadruplex forming sequences identified in a human genome, the possibility of non-specific interactions with G4 ligands cannot be neglected [93]. The attractiveness of G4 ligands may be enhanced by determining a means to increase their ability to specifically bind and stabilize G4 structures. In this paper, we presented the concept of covalent attachment of G4 ligands to G-rich oligoribonucleotides. The use of such conjugates creates the possibility of inducing and stabilizing the bimolecular G-quadruplexes on the G-rich mRNA template in a sequence-specific manner. We showed recently that oligonucleotides with a covalently attached fluorescent carbazole derivative (o-BMVC-Q-ASO) recognized the target site on EGFR mRNA comprising two G-tracts separated by 28 residues by forming a ligand–quadruplex–duplex hybrid structure. The observed level of green fluorescence of an o-BMVC moiety in three different cancer cells correlated well with the amount of EGFR mRNA [52]. We believe that the formation of QDH or Dss structures is dependent on a single nucleotide change in the target sequence (Figure 1), and the possibility to selectively stabilize the G-quadruplex domain by attaching the G4 ligand (Figure 9) may become an attractive alternative therapy for patients with an EGFR-L858R mutation. Furthermore, o-BMVC-C3-oligonucleotides can be used as fluorescent hybridization probes to visualize the single nucleotide EGFR-L858R mutation of mRNA. The use of sequence-guided G4 ligands acting as a G4 stabilizer can minimize the problem of the dynamic nature of the RNA G-quadruplex by slowing the action of G4-helicases [60,94,95]. In our opinion, the results obtained in this work provide new perspectives to solve the problem of multi-target binding of G4 ligands.

5. Conclusions

In this article, we proposed a method to discriminate between two similar target RNA sequences, GT and UT, that differ in one nucleotide only, based on the formation of alternative structures, i.e., quadruplex–duplex hybrids or duplexes with dangling ends, respectively. We designed DNA and RNA G-rich oligonucleotides with the ability to trigger the formation of a G-quadruplex motif only when hybridized to the GT target and leaving an unstructured G-rich fragment when hybridized to the UT sequence. We also showed that the replacement of RNA G-rich oligonucleotides with DNA decreased the thermal stability of the quadruplex–duplex hybrid structure and the duplex with dangling ends. Subsequently, we demonstrated that abasic residues and aliphatic linkers can mimic the external loops of bimolecular RNA G-quadruplexes. The quadruplex–duplex hybrid structures containing these non-nucleotide modifications exhibited similar stability as their unmodified counterparts. Unexpectedly, we noticed that the presence of two abasic modifications in G-rich strands induced the formation of non-canonical G-quadruplexes after hybridization to the UT target. Finally, RNA G-rich oligonucleotides with a covalently attached carbazole derivative, o-BMVC-C3, were shown to selectively bind and stabilize a G-quadruplex domain of QDH. The obtained results advance our knowledge and ability to predict structures adopted by G-rich sequences and can be used as a starting point to design anti-EGFR G-rich antisense oligonucleotides. Moreover, o-BMVC-C3 covalently attached to an oligonucleotide allows fluorescent probe visualization of sequences containing two GGG tracts.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/biom11081236/s1, Figure S1: The purity and homogeneity of o-BMVC-C3-CCC oligonucleotide verified by 15% denaturing gel electrophoresis, Figure S2: CD spectra of CCC/GT and CCC/UT in the presence of 50 mM KCl, 10 mM potassium phosphate and 0.1mM EDTA, pH 6.8; 25 °C, Figure S3: UV-melting profiles of CCC/GT and CCC/UT in the presence of 50 mM KCl, 10 mM potassium phosphate and 0.1mM EDTA, pH 6.8, Figure S4: Monitoring the increasing level of higher-order structure as a function of time. The migration of GT, CCC, and CCC/GT molecules in the presence of 50 mM KCl by non-denaturing PAGE 20%. Gel was post-stained with NMM solution to detect G4 motifs, Figure S5: CD spectra of DNA-CCC/GT and DNA-CCC/UT in the presence of 50 mM KCl, 10 mM potassium phosphate and 0.1mM EDTA, pH 6.8; 25 °C, Figure S6: UV-melting curves of DNA-CCC/GT hybrid in the presence of 50 mM KCl, 10 mM potassium phosphate and 0.1mM EDTA, pH 6.8, Figure S7: Comparison of the imino region of the 1H NMR spectra of CCC/GT, CCC-OMe/GT, paa/GT, Aaa/GT, paC/GT, AaC/GT, pCa/GT, pL2C/GT, pL3C/GT, pL4C/GT (A) and CCC/UT, CCC-OMe/UT, paa/UT, Aaa/UT, paC/UT, AaC/UT, pCa/UT, pL2C/UT, pL3C/UT, pL4C/UT complexes (B) at 25 °C in the presence of 150 mM NaCl, 10 mM sodium phosphate, 0.1 mM EDTA, pH 6.8, Figure S8: Schematic representation of the two possible models of the quadruplex–duplex hybrid, Figure S9: Fluorescence spectra of NMM and NMM mixed with QD/GT (solid lines) and QD/UT (dotted lines) molecules in a buffer containing 50 mM KCl, 10 mM potassium phosphate, pH 6.8 at 25 °C, Figure S10: Comparison of the imino region of 1H NMR spectra of paa/GT (A), Aaa/UT (B), and paa/UT (C) without and with a 1:1 molar ratio of NMM in a buffer containing 50 mM KCl, 10 mM potassium phosphate, pH 6.8 at 25 °C, Figure S11: Migration of CCC/GT, CCC/UT and chemically modified QD/GT and QD/UT complexes in the presence of 150 mM NaCl in non-denaturing PAGE 20%; gels were visualized by UV light (254 nm) to detect all RNAs and post-stained with NMM solution to detect G4 motifs, Table S1: Thermal stability for QDH and Dss structures in sodium phosphate buffer, Figure S12: The UV-melting profiles of QD/GT at 260 and 295 nm in the presence of 50 mM KCl, 10 mM potassium phosphate, 0.1 mM EDTA, pH 6.8 (A) and 150 mM NaCl, 10 mM sodium phosphate, 0.1 mM EDTA, pH 6.8 (B), Figure S13: The UV-melting profiles of QD/UT at 260 and 295 nm in the presence of 50 mM KCl, 10 mM potassium phosphate, 0.1 mM EDTA, pH 6.8 (A) and 150 mM NaCl, 10 mM sodium phosphate, 0.1 mM EDTA, pH 6.8 (B), Figure S14: Temperature dependence of the imino region of 1H NMR spectra for CCC/GT, paa/GT, and pL3C/GT in the presence of 50 mM KCl (A) and 150 mM NaCl (B), Figure S15: Circular dichroism (CD) spectra of QDH in the presence of potassium and sodium cations, Figure S16. Circular dichroism (CD) spectra of molecules indicated in the legend in the presence of 50 mM KCl, 10 mM potassium phosphate, 0.1 mM EDTA, pH 6.8, Figure S17: o-BMVC-C3-CCC oligoribonucleotide synthesis scheme, Figure S18: The UV-melting profiles of CCC/GT, o-BMVC-C3-CCC/GT, and o-BMVC-C3-CCC/UT at 260 nm. The spectra were recorded in the presence of 50 mM KCl, 10 mM potassium phosphate and 0.1 mM EDTA, pH 6.8, Figure S19: Circular dichroism (CD) spectra of o-BMVC-C3-CCC/GT, o-BMVC-C3-CCC/UT, and CCC/GT at 25 °C in the presence of 50 mM KCl, 10 mM potassium phosphate and 0.1 mM EDTA, pH 6.8, Figure S20. Secondary structures of CCC/UT predicted by RNAstructure, Figure S21. Putative models of bimolecular paa/UT and Aaa/UT hybrid.

Author Contributions

Conceptualization and supervision, D.G. and Z.G.; investigation, analyzed the data, methodology, visualization, data curation, D.G. and K.Z.; synthesis of the o-BMVC-C3-NHS ligand, D.B.; selecting the target sequence, P.K.; synthesis of the oligonucleotides, R.K.; writing—original draft preparation, D.G., Z.G. and D.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Science Centre, grant number 2014/13/B/ST5/04144 to Z.G., 2017/01/X/ST5/00589 and 2020/39/D/ST4/03177 to D.G.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Acknowledgments

We thank Karol Pasternak for the synthesis of the DNA oligonucleotide.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Burge, S.; Parkinson, G.N.; Hazel, P.; Todd, A.K.; Neidle, S. Quadruplex DNA: Sequence, topology and structure. Nucleic Acids Res. 2006, 34, 5402–5415. [Google Scholar] [CrossRef] [Green Version]
  2. Parkinson, G.N. Fundamentals of Quadruplex Structures. In Quadruplex Nucleic Acids; Neidle, S., Balasubramanian, S., Eds.; Royal Society of Chemistry: London, UK, 2007; pp. 1–30. [Google Scholar]
  3. Han, J.H.; Kim, J.H.; Kim, S.K.; Jang, Y.J. Conformational change of a G-quadruplex under molecular crowding conditions. J. Biomol. Struct. Dyn. 2020, 38, 2575–2581. [Google Scholar] [CrossRef]
  4. Chilka, P.; Desai, N.; Datta, B. Small Molecule Fluorescent Probes for G- Quadruplex Visualization as Potential Cancer Theranostic Agents. Molecules 2019, 24, 752. [Google Scholar] [CrossRef] [Green Version]
  5. Kwok, C.K.; Merrick, C.J. G-Quadruplexes: Prediction, Characterization, and Biological Application. Trends Biotechnol. 2017, 35, 997–1013. [Google Scholar] [CrossRef] [PubMed]
  6. Dumas, L.; Herviou, P.; Dassi, E.; Cammas, A.; Millevoi, S. G-Quadruplexes in RNA Biology: Recent Advances and Future Directions. Trends Biochem. Sci. 2021, 46, 270–283. [Google Scholar] [CrossRef] [PubMed]
  7. Xiao, C.-D.; Shibata, T.; Yamamoto, Y.; Xu, Y. An intramolecular antiparallel G-quadruplex formed by human telomere RNA. Chem. Commun. 2018, 54, 3944–3946. [Google Scholar] [CrossRef] [PubMed]
  8. Banco, M.T.; Ferré-D’Amaré, A.R. The emerging structural complexity of G-quadruplex RNAs. RNA 2021, 27, 390–402. [Google Scholar] [CrossRef] [PubMed]
  9. Trachman, R.J.; Autour, A.; Jeng, S.C.Y.; Abdolahzadeh, A.; Andreoni, A.; Cojocaru, R.; Garipov, R.; Dolgosheina, E.V.; Knutson, J.R.; Ryckelynck, M.; et al. Structure and functional reselection of the Mango-III fluorogenic RNA aptamer. Nat. Chem. Biol. 2019, 15, 472–479. [Google Scholar] [CrossRef]
  10. Trachman, R.J.; Demeshkina, N.A.; Lau, M.W.L.; Panchapakesan, S.S.S.; Jeng, S.C.Y.; Unrau, P.J.; Ferré-D’Amaré, A.R. Structural basis for high-affinity fluorophore binding and activation by RNA Mango. Nat. Chem. Biol. 2017, 13, 807–813. [Google Scholar] [CrossRef] [Green Version]
  11. Huang, H.; Suslov, N.B.; Li, N.-S.; Shelke, S.A.; Evans, M.E.; Koldobskaya, Y.; Rice, P.A.; Piccirilli, J.A. A G-quadruplex–containing RNA activates fluorescence in a GFP-like fluorophore. Nat. Chem. Biol. 2014, 10, 686–691. [Google Scholar] [CrossRef]
  12. Haase, L.; Karg, B.; Weisz, K. Manipulating DNA G-Quadruplex Structures by Using Guanosine Analogues. ChemBioChem 2019, 20, 985–993. [Google Scholar] [CrossRef]
  13. Sagi, J. G-quadruplexes incorporating modified constituents: A review. J. Biomol. Struct. Dyn. 2014, 32, 477–511. [Google Scholar] [CrossRef] [PubMed]
  14. Lech, C.J.; Li, Z.; Heddi, B.; Phan, A.T. 2′-F-ANA-guanosine and 2′-F-guanosine as powerful tools for structural manipulation of G-quadruplexes. Chem. Commun. 2012, 48, 11425. [Google Scholar] [CrossRef] [PubMed]
  15. Li, Z.; Lech, C.J.; Phan, A.T. Sugar-modified G-quadruplexes: Effects of LNA-, 2′F-RNA– and 2′F-ANA-guanosine chemistries on G-quadruplex structure and stability. Nucleic Acids Res. 2014, 42, 4068–4079. [Google Scholar] [CrossRef] [PubMed]
  16. Risitano, A. Influence of loop size on the stability of intramolecular DNA quadruplexes. Nucleic Acids Res. 2004, 32, 2598–2606. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Rachwal, P.A.; Brown, T.; Fox, K.R. Sequence effects of single base loops in intramolecular quadruplex DNA. FEBS Lett. 2007, 581, 1657–1660. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Gudanis, D.; Popenda, L.; Szpotkowski, K.; Kierzek, R.; Gdaniec, Z. Structural characterization of a dimer of RNA duplexes composed of 8-bromoguanosine modified CGG trinucleotide repeats: A novel architecture of RNA quadruplexes. Nucleic Acids Res. 2016, 44, 2409–2416. [Google Scholar] [CrossRef] [Green Version]
  19. Ji, X.; Sun, H.; Zhou, H.; Xiang, J.; Tang, Y.; Zhao, C. Research progress of RNA quadruplex. Nucleic Acid Ther. 2011, 21, 185–200. [Google Scholar] [CrossRef] [PubMed]
  20. Bhavesh, N.S.; Patel, P.K.; Karthikeyan, S.; Hosur, R.V. Distinctive features in the structure and dynamics of the DNA repeat sequence GGCGGG. Biochem. Biophys. Res. Commun. 2004, 317, 625–633. [Google Scholar] [CrossRef] [PubMed]
  21. Mergny, J.L.; De Cian, A.; Amrane, S.; da Silva, M.W. Kinetics of double-chain reversals bridging contiguous quartets in tetramolecular quadruplexes. Nucleic Acids Res. 2006, 34, 2386–2397. [Google Scholar] [CrossRef] [Green Version]
  22. Malgowska, M.; Gudanis, D.; Kierzek, R.; Wyszko, E.; Gabelica, V.; Gdaniec, Z. Distinctive structural motifs of RNA G-quadruplexes composed of AGG, CGG and UGG trinucleotide repeats. Nucleic Acids Res. 2014, 42, 10196–10207. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Lightfoot, H.L.; Hagen, T.; Tatum, N.J.; Hall, J. The diverse structural landscape of quadruplexes. FEBS Lett. 2019, 593, 2083–2102. [Google Scholar] [CrossRef] [PubMed]
  24. Malgowska, M.; Czajczynska, K.; Gudanis, D.; Tworak, A.; Gdaniec, Z. Overview of the RNA G-quadruplex structures. Acta Biochim. Pol. 2016, 63, 609–621. [Google Scholar] [CrossRef]
  25. Andrałojć, W.; Małgowska, M.; Sarzyńska, J.; Pasternak, K.; Szpotkowski, K.; Kierzek, R.; Gdaniec, Z. Unraveling the structural basis for the exceptional stability of RNA G-quadruplexes capped by a uridine tetrad at the 3′ terminus. RNA 2019, 25, 121–134. [Google Scholar] [CrossRef]
  26. Andrałojć, W.; Pasternak, K.; Sarzyńska, J.; Zielińska, K.; Kierzek, R.; Gdaniec, Z. The origin of the high stability of 3′-terminal uridine tetrads: Contributions of hydrogen bonding, stacking interactions, and steric factors evaluated using modified oligonucleotide analogs. RNA 2020, 26, 2000–2016. [Google Scholar] [CrossRef]
  27. Esposito, V.; Martino, L.; Citarella, G.; Virgilio, A.; Mayol, L.; Giancola, C.; Galeone, A. Effects of abasic sites on structural, thermodynamic and kinetic properties of quadruplex structures. Nucleic Acids Res. 2009, 38, 2069–2080. [Google Scholar] [CrossRef] [PubMed]
  28. Švehlová, K.; Lawrence, M.S.; Bednárová, L.; Curtis, E.A. Altered biochemical specificity of G-quadruplexes with mutated tetrads. Nucleic Acids Res. 2016, 44, 10789–10803. [Google Scholar] [CrossRef]
  29. Endoh, T.; Kawasaki, Y.; Sugimoto, N. Stability of RNA quadruplex in open reading frame determines proteolysis of human estrogen receptor α. Nucleic Acids Res. 2013, 41, 6222–6231. [Google Scholar] [CrossRef] [Green Version]
  30. Endoh, T.; Sugimoto, N. Conformational Dynamics of the RNA G-Quadruplex and its Effect on Translation Efficiency. Molecules 2019, 24, 1613. [Google Scholar] [CrossRef] [Green Version]
  31. Chambers, V.S.; Marsico, G.; Boutell, J.M.; Di Antonio, M.; Smith, G.P.; Balasubramanian, S. High-throughput sequencing of DNA G-quadruplex structures in the human genome. Nat. Biotechnol. 2015, 33, 877–881. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Huppert, J.L. Prevalence of quadruplexes in the human genome. Nucleic Acids Res. 2005, 33, 2908–2916. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Kwok, C.K.; Marsico, G.; Sahakyan, A.B.; Chambers, V.S.; Balasubramanian, S. rG4-seq reveals widespread formation of G-quadruplex structures in the human transcriptome. Nat. Methods 2016, 13, 841–844. [Google Scholar] [CrossRef] [PubMed]
  34. Kharel, P.; Becker, G.; Tsvetkov, V.; Ivanov, P. Properties and biological impact of RNA G-quadruplexes: From order to turmoil and back. Nucleic Acids Res. 2020, 1–22. [Google Scholar] [CrossRef]
  35. Tassinari, M.; Richter, S.N.; Gandellini, P. Biological relevance and therapeutic potential of G-quadruplex structures in the human noncoding transcriptome. Nucleic Acids Res. 2021, 49, 3617–3633. [Google Scholar] [CrossRef]
  36. Kharel, P.; Balaratnam, S.; Beals, N.; Basu, S. The role of RNA G-quadruplexes in human diseases and therapeutic strategies. Wiley Interdiscip. Rev. RNA 2020, 11, 1–20. [Google Scholar] [CrossRef] [PubMed]
  37. Lyons, S.M.; Gudanis, D.; Coyne, S.M.; Gdaniec, Z.; Ivanov, P. Identification of functional tetramolecular RNA G-quadruplexes derived from transfer RNAs. Nat. Commun. 2017, 8, 1127. [Google Scholar] [CrossRef] [Green Version]
  38. Jackowiak, P.; Hojka-Osinska, A.; Gasiorek, K.; Stelmaszczuk, M.; Gudanis, D.; Gdaniec, Z.; Figlerowicz, M. Effects of G-quadruplex topology on translational inhibition by tRNA fragments in mammalian and plant systems in vitro. Int. J. Biochem. Cell Biol. 2017, 92, 148–154. [Google Scholar] [CrossRef]
  39. Lim, K.W.; Jenjaroenpun, P.; Low, Z.J.; Khong, Z.J.; Ng, Y.S.; Kuznetsov, V.A.; Phan, A.T. Duplex stem-loop-containing quadruplex motifs in the human genome: A combined genomic and structural study. Nucleic Acids Res. 2015, 43, 5630–5646. [Google Scholar] [CrossRef] [Green Version]
  40. Nguyen, T.Q.N.; Lim, K.W.; Phan, A.T. Folding Kinetics of G-Quadruplexes: Duplex Stem Loops Drive and Accelerate G-Quadruplex Folding. J. Phys. Chem. B 2020, 124, 5122–5130. [Google Scholar] [CrossRef]
  41. Lim, K.W.; Phan, A.T. Structural basis of DNA quadruplex-duplex junction formation. Angew. Chemie Int. Ed. 2013, 52, 8566–8569. [Google Scholar] [CrossRef]
  42. Lim, K.W.; Khong, Z.J.; Phan, A.T. Thermal stability of DNA quadruplex-duplex hybrids. Biochemistry 2014, 53, 247–257. [Google Scholar] [CrossRef]
  43. Butovskaya, E.; Heddi, B.; Bakalar, B.; Richter, S.N.; Phan, A.T. Major G-Quadruplex Form of HIV-1 LTR Reveals a (3 + 1) Folding Topology Containing a Stem-Loop. J. Am. Chem. Soc. 2018, 140, 13654–13662. [Google Scholar] [CrossRef] [Green Version]
  44. Lim, K.W.; Nguyen, T.Q.N.; Phan, A.T. Joining of Multiple Duplex Stems at a Single Quadruplex Loop. J. Am. Chem. Soc. 2014, 136, 17969–17973. [Google Scholar] [CrossRef] [PubMed]
  45. Ngoc Nguyen, T.Q.; Lim, K.W.; Phan, A.T. Duplex formation in a G-quadruplex bulge. Nucleic Acids Res. 2020, 1–9. [Google Scholar] [CrossRef]
  46. Vianney, Y.M.; Preckwinkel, P.; Mohr, S.; Weisz, K. Quadruplex–Duplex Junction: A High-Affinity Binding Site for Indoloquinoline Ligands. Chem. A Eur. J. 2020, 26, 16910–16922. [Google Scholar] [CrossRef]
  47. Karg, B.; Mohr, S.; Weisz, K. Duplex-Guided Refolding into Novel G-Quadruplex (3+1) Hybrid Conformations. Angew. Chemie Int. Ed. 2019, 58, 11068–11071. [Google Scholar] [CrossRef]
  48. Greco, M.L.; Kotar, A.; Rigo, R.; Cristofari, C.; Plavec, J.; Sissi, C. Coexistence of two main folded G-quadruplexes within a single G-rich domain in the EGFR promoter. Nucleic Acids Res. 2017, 45, 10132–10142. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Tan, D.J.Y.; Winnerdy, F.R.; Lim, K.W.; Phan, A.T. Coexistence of two quadruplex–duplex hybrids in the PIM1 gene. Nucleic Acids Res. 2020, 48, 11162–11171. [Google Scholar] [CrossRef] [PubMed]
  50. Hagihara, M.; Yamauchi, L.; Seo, A.; Yoneda, K.; Senda, M.; Nakatani, K. Antisense-induced guanine quadruplexes inhibit reverse transcription by HIV-1 reverse transcriptase. J. Am. Chem. Soc. 2010, 132, 11171–11178. [Google Scholar] [CrossRef]
  51. Bhattacharyya, D.; Nguyen, K.; Basu, S. Rationally Induced RNA:DNA G-Quadruplex Structures Elicit an Anticancer Effect by Inhibiting Endogenous eIF-4E Expression. Biochemistry 2014, 53, 5461–5470. [Google Scholar] [CrossRef]
  52. Gudanis, D.; Kaniowski, D.; Kulik, K.; Baranowski, D.; Gdaniec, Z.; Nawrot, B. Formation of an RNA Quadruplex-Duplex Hybrid in Living Cells between mRNA of the Epidermal Growth Factor Receptor (EGFR) and a G-Rich Antisense Oligoribonucleotide. Cells 2020, 9, 2375. [Google Scholar] [CrossRef]
  53. Tan, D.J.Y.; Das, P.; Winnerdy, F.R.; Lim, K.W.; Phan, A.T. Guanine anchoring: A strategy for specific targeting of a G-quadruplex using short PNA, LNA and DNA molecules. Chem. Commun. 2020, 56, 5897–5900. [Google Scholar] [CrossRef] [PubMed]
  54. Russo Krauss, I.; Ramaswamy, S.; Neidle, S.; Haider, S.; Parkinson, G.N. Structural Insights into the Quadruplex-Duplex 3′ Interface Formed from a Telomeric Repeat: A Potential Molecular Target. J. Am. Chem. Soc. 2016, 138, 1226–1233. [Google Scholar] [CrossRef] [Green Version]
  55. Asamitsu, S.; Obata, S.; Phan, A.T.; Hashiya, K.; Bando, T.; Sugiyama, H. Simultaneous Binding of Hybrid Molecules Constructed with Dual DNA-Binding Components to a G-Quadruplex and Its Proximal Duplex. Chem. A Eur. J. 2018, 24, 4428–4435. [Google Scholar] [CrossRef]
  56. Nguyen, T.Q.N.; Lim, K.W.; Phan, A.T. A Dual-Specific Targeting Approach Based on the Simultaneous Recognition of Duplex and Quadruplex Motifs. Sci. Rep. 2017, 7, 11969. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Xu, S.; Li, Q.; Xiang, J.; Yang, Q.; Sun, H.; Guan, A.; Wang, L.; Liu, Y.; Yu, L.; Shi, Y.; et al. Directly lighting up RNA G-quadruplexes from test tubes to living human cells. Nucleic Acids Res. 2015, 43, 9575–9586. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Murat, P.; Marsico, G.; Herdy, B.; Ghanbarian, A.; Portella, G.; Balasubramanian, S. RNA G-quadruplexes at upstream open reading frames cause DHX36- and DHX9-dependent translation of human mRNAs. Genome Biol. 2018, 19, 229. [Google Scholar] [CrossRef]
  59. Benhalevy, D.; Gupta, S.K.; Danan, C.H.; Ghosal, S.; Sun, H.-W.; Kazemier, H.G.; Paeschke, K.; Hafner, M.; Juranek, S.A. The Human CCHC-type Zinc Finger Nucleic Acid-Binding Protein Binds G-Rich Elements in Target mRNA Coding Sequences and Promotes Translation. Cell Rep. 2017, 18, 2979–2990. [Google Scholar] [CrossRef] [Green Version]
  60. Chen, X.-C.; Chen, S.-B.; Dai, J.; Yuan, J.-H.; Ou, T.-M.; Huang, Z.-S.; Tan, J.-H. Tracking the Dynamic Folding and Unfolding of RNA G-Quadruplexes in Live Cells. Angew. Chemie Int. Ed. 2018, 57, 4702–4706. [Google Scholar] [CrossRef]
  61. Biffi, G.; Di Antonio, M.; Tannahill, D.; Balasubramanian, S. Visualization and selective chemical targeting of RNA G-quadruplex structures in the cytoplasm of human cells. Nat. Chem. 2013, 6, 75–80. [Google Scholar] [CrossRef]
  62. Lyu, K.; Chen, S.B.; Chan, C.Y.; Tan, J.H.; Kwok, C.K. Structural analysis and cellular visualization of APP RNA G-quadruplex. Chem. Sci. 2019, 10, 11095–11102. [Google Scholar] [CrossRef] [Green Version]
  63. Chen, S.-B.; Hu, M.-H.; Liu, G.-C.; Wang, J.; Ou, T.-M.; Gu, L.-Q.; Huang, Z.-S.; Tan, J.-H. Visualization of NRAS RNA G-Quadruplex Structures in Cells with an Engineered Fluorogenic Hybridization Probe. J. Am. Chem. Soc. 2016, 138, 10382–10385. [Google Scholar] [CrossRef] [PubMed]
  64. Collie, G.W.; Parkinson, G.N.; Moyzis, R.K.; Buckingham, J.M.; Cram, L.S.; Dani, M.; Deaven, L.L.; Jones, M.D.; Meyne, J.; Ratliff, R.L.; et al. The application of DNA and RNA G-quadruplexes to therapeutic medicines. Chem. Soc. Rev. 2011, 40, 5867. [Google Scholar] [CrossRef] [PubMed]
  65. Millevoi, S.; Moine, H.; Vagner, S. G-quadruplexes in RNA biology. Wiley Interdiscip. Rev. RNA 2012, 3, 495–507. [Google Scholar] [CrossRef] [PubMed]
  66. Agarwala, P.; Pandey, S.; Maiti, S. The tale of RNA G-quadruplex. Org. Biomol. Chem. 2015, 13, 5570–5585. [Google Scholar] [CrossRef] [PubMed]
  67. Kumari, S.; Bugaut, A.; Huppert, J.L.; Balasubramanian, S. An RNA G-quadruplex in the 5′ UTR of the NRAS proto-oncogene modulates translation. Nat. Chem. Biol. 2007, 3, 218–221. [Google Scholar] [CrossRef] [Green Version]
  68. Fay, M.M.; Lyons, S.M.; Ivanov, P. RNA G-Quadruplexes in Biology: Principles and Molecular Mechanisms. J. Mol. Biol. 2017, 429, 2127–2147. [Google Scholar] [CrossRef]
  69. Tian, T.; Chen, Y.Q.; Wang, S.R.; Zhou, X. G-Quadruplex: A Regulator of Gene Expression and Its Chemical Targeting. Chem 2018, 4, 1314–1344. [Google Scholar] [CrossRef] [Green Version]
  70. Adrian, M.; Heddi, B.; Phan, A.T. NMR spectroscopy of G-quadruplexes. Methods 2012, 57, 11–24. [Google Scholar] [CrossRef]
  71. Kreig, A.; Calvert, J.; Sanoica, J.; Cullum, E.; Tipanna, R.; Myong, S.; Gqs, G. G-quadruplex formation in double strand DNA probed by NMM and CV fluorescence. Nucleic Acids Res. 2015, 43, 7961–7970. [Google Scholar] [CrossRef]
  72. Sabharwal, N.C.; Savikhin, V.; Turek-Herman, J.R.; Nicoludis, J.M.; Szalai, V.A.; Yatsunyk, L.A. N -methylmesoporphyrin IX fluorescence as a reporter of strand orientation in guanine quadruplexes. FEBS J. 2014, 281, 1726–1737. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Kypr, J.; Kejnovska, I.; Renciuk, D.; Vorlickova, M. Circular dichroism and conformational polymorphism of DNA. Nucleic Acids Res. 2009, 37, 1713–1725. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. Granqvist, L.; Virta, P. Characterization of G-Quadruplex/Hairpin Transitions of RNAs by 19 F NMR Spectroscopy. Chem. A Eur. J. 2016, 22, 15360–15372. [Google Scholar] [CrossRef]
  75. Kejnovská, I.; Bednářová, K.; Renčiuk, D.; Dvořáková, Z.; Školáková, P.; Trantírek, L.; Fiala, R.; Vorlíčková, M.; Sagi, J. Clustered abasic lesions profoundly change the structure and stability of human telomeric G-quadruplexes. Nucleic Acids Res. 2017, 45, 4294–4305. [Google Scholar] [CrossRef] [PubMed]
  76. Kejnovská, I.; Stadlbauer, P.; Trantírek, L.; Renčiuk, D.; Gajarský, M.; Krafčík, D.; Palacký, J.; Bednářová, K.; Šponer, J.; Mergny, J.-L.; et al. G-Quadruplex formation by DNA sequences deficient in guanines: Two tetrad parallel quadruplexes do not fold intramolecularly. Chem. A Eur. J. 2021. [Google Scholar] [CrossRef]
  77. Uddin, M.K. Phosphorylation at 5’ end of guanosine stretches inhibits dimerization of G-quadruplexes and formation of a G-quadruplex interferes with the enzymatic activities of DNA enzymes. Nucleic Acids Res. 2004, 32, 4618–4629. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  78. Asamitsu, S.; Bando, T.; Sugiyama, H. Ligand Design to Acquire Specificity to Intended G-Quadruplex Structures. Chem. A Eur. J. 2019, 25, 417–430. [Google Scholar] [CrossRef]
  79. Tseng, T.-Y.; Chien, C.-H.; Chu, J.-F.; Huang, W.-C.; Lin, M.-Y.; Chang, C.-C.; Chang, T.-C. Fluorescent probe for visualizing guanine-quadruplex DNA by fluorescence lifetime imaging microscopy. J. Biomed. Opt. 2013, 18, 101309. [Google Scholar] [CrossRef] [Green Version]
  80. Tseng, T.-Y.; Chu, I.-T.; Lin, S.-J.; Li, J.; Chang, T.-C. Binding of Small Molecules to G-quadruplex DNA in Cells Revealed by Fluorescence Lifetime Imaging Microscopy of o-BMVC Foci. Molecules 2018, 24, 35. [Google Scholar] [CrossRef] [Green Version]
  81. Koralewska, N.; Szczepanska, A.; Ciechanowska, K.; Wojnicka, M.; Pokornowska, M.; Milewski, M.; Gudanis, D.; Baranowski, D.; Nithin, C.; Bujnicki, J.; et al. RNA and DNA G-quadruplexes bind to human Dicer and inhibit its activity. Cell. Mol. Life Sci. 2020, 1–24. [Google Scholar] [CrossRef]
  82. Winnerdy, F.R.; Phan, A.T. Quadruplex structure and diversity. Annu. Rep. Med. Chem. 2020, 54, 45–73. [Google Scholar] [CrossRef]
  83. Pandey, S.; Agarwala, P.; Maiti, S. Effect of Loops and G-Quartets on the Stability of RNA G-Quadruplexes. J. Phys. Chem. B 2013, 117, 6896–6905. [Google Scholar] [CrossRef] [PubMed]
  84. Fotticchia, I.; Amato, J.; Pagano, B.; Novellino, E.; Petraccone, L.; Giancola, C. How are thermodynamically stable G-quadruplex-duplex hybrids? J. Therm. Anal. Calorim. 2015, 121, 1121–1127. [Google Scholar] [CrossRef]
  85. Agarwala, P.; Kumar, S.; Pandey, S.; Maiti, S. Human Telomeric RNA G-Quadruplex Response to Point Mutation in the G-Quartets. J. Phys. Chem. B 2015, 119, 4617–4627. [Google Scholar] [CrossRef] [PubMed]
  86. Gros, J.; Rosu, F.; Amrane, S.; De Cian, A.; Gabelica, V.; Lacroix, L.; Mergny, J.-L. Guanines are a quartet’s best friend: Impact of base substitutions on the kinetics and stability of tetramolecular quadruplexes. Nucleic Acids Res. 2007, 35, 3064–3075. [Google Scholar] [CrossRef] [PubMed]
  87. Wieland, M.; Hartig, J.S. RNA Quadruplex-Based Modulation of Gene Expression. Chem. Biol. 2007, 14, 757–763. [Google Scholar] [CrossRef] [Green Version]
  88. Halder, K.; Wieland, M.; Hartig, J.S. Predictable suppression of gene expression by 5′-UTR-based RNA quadruplexes. Nucleic Acids Res. 2009, 37, 6811–6817. [Google Scholar] [CrossRef]
  89. He, G.-X.; Williams, J.P.; Postich, M.J.; Swaminathan, S.; Shea, R.G.; Terhorst, T.; Law, V.S.; Mao, C.T.; Sueoka, C.; Coutré, S.; et al. In Vitro and in Vivo Activities of Oligodeoxynucleotide-Based Thrombin Inhibitors Containing Neutral Formacetal Linkages. J. Med. Chem. 1998, 41, 4224–4231. [Google Scholar] [CrossRef]
  90. Zhang, T.; Wan, B.; Zhao, Y.; Li, C.; Liu, H.; Lv, T.; Zhan, P.; Song, Y. Treatment of uncommon EGFR mutations in non-small cell lung cancer: New evidence and treatment. Transl. Lung Cancer Res. 2019, 8, 302–316. [Google Scholar] [CrossRef]
  91. Asamitsu, S.; Obata, S.; Yu, Z.; Bando, T.; Sugiyama, H. Recent Progress of Targeted G-Quadruplex-Preferred Ligands Toward Cancer Therapy. Molecules 2019, 24, 429. [Google Scholar] [CrossRef] [Green Version]
  92. Tao, Y.; Zheng, Y.; Zhai, Q.; Wei, D. Recent advances in the development of small molecules targeting RNA G-quadruplexes for drug discovery. Bioorg. Chem. 2021, 110, 104804. [Google Scholar] [CrossRef] [PubMed]
  93. Armas, P.; David, A.; Calcaterra, N.B. Transcriptional control by G-quadruplexes: In vivo roles and perspectives for specific intervention. Transcription 2017, 8, 21–25. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  94. Guo, J.U.; Bartel, D.P. RNA G-quadruplexes are globally unfolded in eukaryotic cells and depleted in bacteria. Science 2016, 353, aaf5371. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  95. Chen, M.C.; Murat, P.; Abecassis, K.; Ferre-D’Amare, A.R.; Balasubramanian, S. Insights into the mechanism of a G-quadruplex-unwinding DEAH-box helicase. Nucleic Acids Res. 2015, 43, 2223–2231. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Schematic presentation of the recognition of GT and UT targets by the CCC oligonucleotide. CCC oligonucleotide (A) and target sequence with two GGG-tracts, GT (B) forms a quadruplex–duplex structural hybrid (QDH) (D). CCC oligonucleotide hybridized to the UT target (C) forms a duplex with long dangling ends (Dss) (E).
Figure 1. Schematic presentation of the recognition of GT and UT targets by the CCC oligonucleotide. CCC oligonucleotide (A) and target sequence with two GGG-tracts, GT (B) forms a quadruplex–duplex structural hybrid (QDH) (D). CCC oligonucleotide hybridized to the UT target (C) forms a duplex with long dangling ends (Dss) (E).
Biomolecules 11 01236 g001
Figure 2. The imino region of the 1H NMR spectra of UT, GT, CCC, 2Q, DX, CCC/GT, and CCC/UT, recorded in 90% H2O/10% D2O (v/v) in the presence of 50 mM KCl, 10 mM potassium phosphate, and 0.1 mM EDTA, pH 6.8; 25 °C (A). Proposed secondary structures (2Q, QDH, Dss) and predicted by RNAstructure (hairpins, DX) (B).
Figure 2. The imino region of the 1H NMR spectra of UT, GT, CCC, 2Q, DX, CCC/GT, and CCC/UT, recorded in 90% H2O/10% D2O (v/v) in the presence of 50 mM KCl, 10 mM potassium phosphate, and 0.1 mM EDTA, pH 6.8; 25 °C (A). Proposed secondary structures (2Q, QDH, Dss) and predicted by RNAstructure (hairpins, DX) (B).
Biomolecules 11 01236 g002
Figure 3. Analysis of the migration of CCC, GT, UT, CCC/GT, and CCC/UT by non-denaturing PAGE 20%; gel was visualized by UV light (254 nm) to detect all RNAs (A) and post-stained with NMM solution to detect G4 motifs (B).
Figure 3. Analysis of the migration of CCC, GT, UT, CCC/GT, and CCC/UT by non-denaturing PAGE 20%; gel was visualized by UV light (254 nm) to detect all RNAs (A) and post-stained with NMM solution to detect G4 motifs (B).
Biomolecules 11 01236 g003
Figure 4. The imino region of the 1H NMR spectra of DNA-CCC/UT, DNA-CCC/GT, DNA-CCC, and UT, recorded in 90% H2O/10% D2O (v/v) in the presence of 50 mM KCl, 10 mM potassium phosphate, and 0.1 mM EDTA, pH 6.8; 25 °C (A). Red stars—imino signals of DNA-CCC, blue dots—imino signals of UT. Proposed secondary structures (2Q, QDH, Dss) and predicted by RNAstructure (hairpins, DX) (B).
Figure 4. The imino region of the 1H NMR spectra of DNA-CCC/UT, DNA-CCC/GT, DNA-CCC, and UT, recorded in 90% H2O/10% D2O (v/v) in the presence of 50 mM KCl, 10 mM potassium phosphate, and 0.1 mM EDTA, pH 6.8; 25 °C (A). Red stars—imino signals of DNA-CCC, blue dots—imino signals of UT. Proposed secondary structures (2Q, QDH, Dss) and predicted by RNAstructure (hairpins, DX) (B).
Biomolecules 11 01236 g004
Figure 5. Analysis of the migration of CCC/GT, CCC/UT, DNA-CCC/GT, DNA-CCC/UT, and DNA-CCC by non-denaturing PAGE 20%; gel was visualized by UV light (254 nm) to detect all RNAs and post-stained with NMM solution to detect G4 motifs. Red star—DNA-CCC, blue star—UT, HOS—higher-order structure, QDH—quadruplex-duplex hybrid. The first lane corresponds to the DNA ladder marker in the range of 10 to 100 bp (W).
Figure 5. Analysis of the migration of CCC/GT, CCC/UT, DNA-CCC/GT, DNA-CCC/UT, and DNA-CCC by non-denaturing PAGE 20%; gel was visualized by UV light (254 nm) to detect all RNAs and post-stained with NMM solution to detect G4 motifs. Red star—DNA-CCC, blue star—UT, HOS—higher-order structure, QDH—quadruplex-duplex hybrid. The first lane corresponds to the DNA ladder marker in the range of 10 to 100 bp (W).
Biomolecules 11 01236 g005
Figure 6. The chemical formulas and the sites of chemical modifications in the QD strand.
Figure 6. The chemical formulas and the sites of chemical modifications in the QD strand.
Biomolecules 11 01236 g006
Figure 7. Comparison of imino region of the 1H NMR spectra of CCC/GT, CCC-OMe/GT, paa/GT, Aaa/GT, paC/GT, AaC/GT, pCa/GT, pL2C/GT, pL3C/GT, pL4C/GT (A) and CCC/UT, CCC-OMe/UT, paa/UT, Aaa/UT, paC/UT, AaC/UT, pCa/UT, pL2C/UT, pL3C/UT, pL4C/U complexes (B) at 25 °C in the presence of 50 mM KCl, 10 mM potassium phosphate, 0.1 mM EDTA, pH 6.8.
Figure 7. Comparison of imino region of the 1H NMR spectra of CCC/GT, CCC-OMe/GT, paa/GT, Aaa/GT, paC/GT, AaC/GT, pCa/GT, pL2C/GT, pL3C/GT, pL4C/GT (A) and CCC/UT, CCC-OMe/UT, paa/UT, Aaa/UT, paC/UT, AaC/UT, pCa/UT, pL2C/UT, pL3C/UT, pL4C/U complexes (B) at 25 °C in the presence of 50 mM KCl, 10 mM potassium phosphate, 0.1 mM EDTA, pH 6.8.
Biomolecules 11 01236 g007
Figure 8. Analysis of the migration of CCC/GT, CCC/UT, and chemically modified QD/GT and QD/UT complexes in potassium phosphate buffer by non-denaturing PAGE 20%; gels were visualized by UV light (254 nm) to detect all RNAs (A) and post-stained with NMM solution to detect G4 motifs (B).
Figure 8. Analysis of the migration of CCC/GT, CCC/UT, and chemically modified QD/GT and QD/UT complexes in potassium phosphate buffer by non-denaturing PAGE 20%; gels were visualized by UV light (254 nm) to detect all RNAs (A) and post-stained with NMM solution to detect G4 motifs (B).
Biomolecules 11 01236 g008
Figure 9. The recognition of GT and UT targets by a model L–CCC oligonucleotide resulting in the formation of a ligand–quadruplex–duplex structural hybrid (L–CCC/GT) or duplex with long dangling ends (L–CCC/UT).
Figure 9. The recognition of GT and UT targets by a model L–CCC oligonucleotide resulting in the formation of a ligand–quadruplex–duplex structural hybrid (L–CCC/GT) or duplex with long dangling ends (L–CCC/UT).
Biomolecules 11 01236 g009
Figure 10. The imino region of the 1H NMR spectra of o-BMVC-C3-CCC/GT, o-BMVC-C3-CCC/UT, and CCC/GT recorded in 90% H2O/10% D2O (v/v); 25 °C (A). Migration of o-BMVC-C3-CCC, o-BMVC-C3-CCC/GT, and o-BMVC-C3-CCC/GT in non-denaturing PAGE 20%. Gel was visualized by UV light (254 nm) to detect all RNAs (B, left), and the same gel was exposed to 473 nm light to visualize the RNA structure adopting the ligand–quadruplex domain (B, right) in the presence of 50 mM KCl, 10 mM potassium phosphate and 0.1mM EDTA, pH 6.8.
Figure 10. The imino region of the 1H NMR spectra of o-BMVC-C3-CCC/GT, o-BMVC-C3-CCC/UT, and CCC/GT recorded in 90% H2O/10% D2O (v/v); 25 °C (A). Migration of o-BMVC-C3-CCC, o-BMVC-C3-CCC/GT, and o-BMVC-C3-CCC/GT in non-denaturing PAGE 20%. Gel was visualized by UV light (254 nm) to detect all RNAs (B, left), and the same gel was exposed to 473 nm light to visualize the RNA structure adopting the ligand–quadruplex domain (B, right) in the presence of 50 mM KCl, 10 mM potassium phosphate and 0.1mM EDTA, pH 6.8.
Biomolecules 11 01236 g010
Table 1. List of RNA and DNA Sequences (5′3′) Analyzed in This Work.
Table 1. List of RNA and DNA Sequences (5′3′) Analyzed in This Work.
NameRNA G-Rich Oligonucleotide Sequences (QD) 5′-3′Bimolecular Complexes
CCCC-GGG-C-GGG-C-CUUCAAGUCCGGCACCC/GT, CCC/UT
CCC-OMeC-GGG-C-GGG-C-(CUUCAAGUCCGGCA)2′OMeCCC-OMe/GT, CCC-OMe/UT
paap-GGG-a-GGG-a-CUUCAAGUCCGGCApaa/GT, paa/UT
AaaA-GGG-a-GGG-a-CUUCAAGUCCGGCAAaa/GT, Aaa/UT
paCp-GGG-a-GGG-C-CUUCAAGUCCGGCApaC/GT, paC/UT
AaCA-GGG-a-GGG-C-CUUCAAGUCCGGCAAaC/GT, AaC/UT
pCap-GGG-C-GGG-a-CUUCAAGUCCGGCApCa/GT, pCa/UT
pL2Cp-GGG-L2-GGG-a-CUUCAAGUCCGGCApL2C/GT, pL2C/UT
pL3Cp-GGG-L3-GGG-a-CUUCAAGUCCGGCApL3C/GT, pL3C/UT
pL4Cp-GGG-L4-GGG-a-CUUCAAGUCCGGCApL4C/GT, pL4C/UT
DNA G-rich Oligonucleotide Sequence
DNA-CCCC-GGG-C-GGG-C-CTTCAAGTCCGGCADNA-CCC/GT, DNA-CCC/UT
G4 Ligand Bearing RNA G-rich Oligonucleotide
L-CCCo-BMVC-C3-aminolinker-C-GGG-C-GGG-C-CUUCAAGUCCGGCAL-CCC/GT, L-CCC/UT
Target Oligonucleotide Sequences
GTUGCCGGACUUGAAG-UA-GGGCGGGC
UTUGCCGGACUUGAAG-UA-GGGCUGGC
Control Oligonucleotide Sequences
DXUGCCGGACUUGAAG/CUUCAAGUCCGGCA
2QCGGGCGGGC
p: phosphate group; a: abasic; L2: 1,2-ethanediol; L3: 1,3-propanediol; L4: 1,4-butanediol; o-BMVC: G-quadruplex binder.
Table 2. Thermal stability for QDH and Dss structures in potassium phosphate buffer.
Table 2. Thermal stability for QDH and Dss structures in potassium phosphate buffer.
QD/GT60 mM K+QD/UT60 mM K+
Tm (°C) Tm (°C)
CCC/GT69.4CCC/UT61.0
CCC-OMe/GT71.7CCC-OMe/UT63.5
paa/GT68.0paa/UT61.5
Aaa/GT69.0Aaa/UT61.0
paC/GT67.5paC/GT61.0
AaC/GT68.5AaC/GT62.0
pCa/GT67.4pCa/GT60.5
pL2C/GT68.0pL2C/GT61.5
pL3C/GT68.5pL3C/GT61.5
pL4C/GT68.5pL4C/GT60.5
DNA-CCC/GT61.4DNA-CCC/UT59.5
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Gudanis, D.; Zielińska, K.; Baranowski, D.; Kierzek, R.; Kozłowski, P.; Gdaniec, Z. Impact of a Single Nucleotide Change or Non-Nucleoside Modifications in G-Rich Region on the Quadruplex–Duplex Hybrid Formation. Biomolecules 2021, 11, 1236. https://doi.org/10.3390/biom11081236

AMA Style

Gudanis D, Zielińska K, Baranowski D, Kierzek R, Kozłowski P, Gdaniec Z. Impact of a Single Nucleotide Change or Non-Nucleoside Modifications in G-Rich Region on the Quadruplex–Duplex Hybrid Formation. Biomolecules. 2021; 11(8):1236. https://doi.org/10.3390/biom11081236

Chicago/Turabian Style

Gudanis, Dorota, Karolina Zielińska, Daniel Baranowski, Ryszard Kierzek, Piotr Kozłowski, and Zofia Gdaniec. 2021. "Impact of a Single Nucleotide Change or Non-Nucleoside Modifications in G-Rich Region on the Quadruplex–Duplex Hybrid Formation" Biomolecules 11, no. 8: 1236. https://doi.org/10.3390/biom11081236

APA Style

Gudanis, D., Zielińska, K., Baranowski, D., Kierzek, R., Kozłowski, P., & Gdaniec, Z. (2021). Impact of a Single Nucleotide Change or Non-Nucleoside Modifications in G-Rich Region on the Quadruplex–Duplex Hybrid Formation. Biomolecules, 11(8), 1236. https://doi.org/10.3390/biom11081236

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop