key: cord-0734583-xdmfpljp authors: Molnár, Orsolya Réka; Végh, András; Somkuti, Judit; Smeller, László title: Characterization of a G-quadruplex from hepatitis B virus and its stabilization by binding TMPyP4, BRACO19 and PhenDC3 date: 2021-12-01 journal: Sci Rep DOI: 10.1038/s41598-021-02689-y sha: 280d6b157753698acc19bfb80022c95e4bd94295 doc_id: 734583 cord_uid: xdmfpljp 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. contacted the virus and currently about 350-400 million people are infected 35 . In 2010, HBV infection was the tenth leading cause of death globally 36 . HBV is the only member of the Hepadnaviridae family that can infect human cells; it contains partially double-stranded DNA (dsDNA) and possesses reverse transcriptase (RT) enzyme activity. It has a small circular DNA genome with the length of approximately 3.2 kilobase pairs 35, 37 . HBV has four overlapping reading frames that encode seven proteins (including S for surface or envelope gene, C for the core gene and P for the polymerase gene) 37, 38 . It is classified into ten genotypes (A-J) based on nucleotide variation 36, 37 . Lavezzo et al. 31 analysed genomes of several viruses. For the HBV the highest G-score was associated to the GGC TGG GGC TTG GTC ATG GGC CAT CAG (NC_003977.2:1204..1230 (+ strand)) sequence. This sequence can be found in the coding region of the polymerase protein. This paper focuses on the stability of this GQ and on the targeting of this GQ with ligands developed for cancer therapy: TMPyP4, BRACO19 and PhenDC3. These ligands were developed with the aim to stabilize the telomere and oncogene promoter GQs in order to reduce the cancerous development. We investigated whether these ligands also influence the viral HepB GQ. Materials. The oligonucleotide GGC TGG GGC TTG GTC ATG GGC CAT CAG was named HepB and purchased from IDT (NY, USA) and Sigma-Aldrich Kft (Hungary). The oligo labeled with the FRET pair FAM and TAMRA (HepB_FRET) was also purchased from the same sources. The oligos were obtained from the manufacturers in lyophilized form. For the fluorescence experiments the HepB_FRET samples were first dissolved in MilliQ water in a concentration of 100 μM, according to the suggestion of IDT. This stock solution was kept frozen and diluted with an appropriate buffer during sample preparation. The final concentration of the oligos in the samples was 1 μM, unless stated otherwise. For the heating experiments 100 mM K-phosphate and Na-phosphate buffers (pH 7.4) were used, because of their insensitivity to temperature. For experiments with Li + and Rb + Tris buffer was used. In this case the pH of the solution was 7.4 at the transition temperature, taking into account the − 0.024 pH unit/°C drift value (calculated from the Sigma product information page). The ion titration experiments were performed in 1 mM Tris-HCl buffer (pH 7.4). TMPyP4 was purchased EMD (USA) and ChemCruz (Dallas, USA), BRACO19 and PhenDC3 were purchased from Sigma-Aldrich Kft. (Hungary). In the infrared experiments much higher concentration of HepB was required. The oligo was dissolved in D 2 O based Bis-Tris buffer (100 mM, pD 7.4) in a concentration of 20 mg/ml. The metal ion concentrations were 100 mM. For the CD experiments the same buffers were used as for the fluorescence experiments. The concentration of the oligo was 12 µM. All chemicals not specified above differently were purchased from Sigma-Aldrich. Fluorescence experiments were performed as described earlier in detail 39, 40 . Two spectrom- www.nature.com/scientificreports/ In the absorption spectroscopic measurements the TMPyP4 concentration was determined by the absorption of the Soret band using the extinction coefficient of 2.26 × 10 5 M −1 cm −141 . FTIR spectra were measured with a Bruker Vertex 80v spectrometer. 256 spectra were averaged at 2 cm −1 resolution. D 2 O buffer was used as solvent, in order to avoid the large absorption band of water around 1640 cm −1 . The samples were measured in a temperature controlled diamond anvil cell (Diacell, UK) to reduce the sample volume. CD spectra were measured using a JASCO 1500 CD/LD spectrometer in CD mode. The sample was injected into a cuvette with 1 mm path length. Three spectra recorded in 1 nm steps were averaged for each sample. The average spectra were smoothed by a boxcar function. These spectra were recorded at room temperature. Absorption spectra were measured by a Cary4 UV-Vis spectrometer. Spectra were smoothed with 9-point Savitzky-Golay filter 42 . Spectral parameters, like peak position and amplitude were evaluated with the PIW program using the Savitzky-Golay peak finding algorithm 42, 43 . Here c L and c G denote the total concentrations of the ligand and the GQ. We used this equation for the determination of K d , since c L and c G both change during the titration experiment. Similarly: In the metal ion binding experiments, the labeled oligo gave the measured fluorescent signal. In this case the donor fluorescence was used to measure the free (unfolded) oligo concentration. [G] was obtained from the fluorescence signal in the following way: Similarly, if the ligand (TMPyP4, BRACO19 and PhenDC3) did have fluorescence signal, we used: where F means the fluorescence intensity at the given concentration, F bound is the intensity in case of complete binding and a is a fitted parameter. In the metal ion titration experiments the ion concentration was increased, while the oligo concentration decreased only due to dilution of the sample. In the ligand binding experiments the GQ concentration was increased, while the total ligand concentration of the solution changed only due to the dilution. It has to be mentioned that the fitting is mathematically ill-conditioned, if the resulting K d is smaller than c L . In this case the K d value was accepted as the highest estimate of the dissociation constant. Temperature stability experiments. The unfolding temperature of the GQ was determined from the fit of the following equation: Here, y is the physical parameter to be fitted (e.g. fluorescence intensity, or ratio of fluorescence intensities at two different wavelengths), a and b are the parameters describing the linear dependence of y(T) below the transition, T is the thermodynamic temperature, Δa and Δb are the changes of a and b during the transition, ΔH is the enthalpy change, R is the universal gas constant, and T m is the transition midpoint. Characterization of HepB structure and stability by FRET and FTIR measurements. FRET has been proven to be suitable for studying GQ conformational changes 6, 44 . Figure 1a and Supplementary Fig. S1 shows the fluorescence spectra of the HepB_FRET at selected temperatures in presence of 140 mM K + ion. Increase of the donor intensity relative to the acceptor intensity shows the loss of the fluorescence energy transfer due to increased distance between the fluorophores. The increase of the distance indicates unfolding of the GQ structure. The donor intensity vs. temperature is plotted in Fig. 1b . Fitting the Eq. (6) results in a transition temperature of T m = 53.1 °C. The measured temperature is comparable with the typical values measured for twoquartet GQs. One of the most known two-quartet GQ is the thrombin binding aptamer (TBA). Sugimoto's lab studied this GQ extensively using optical spectroscopic techniques. Their measured T m value in KCl is 51.2 °C, which is very close to our results. Although the stability of GQs is also influenced by the composition and length of the loops 45,46 typical unfolding temperature of two staged GQs lies around 50°. In contrary, 3-quartet GQs like Htel, c-MYC or KIT1 have an unfolding point higher than 60 °C in presence of K + ions 39, 40 . www.nature.com/scientificreports/ Observing the sequence of HepB we can hypothesize a two-quartet structure, since the shortest guanine repeats contain only two guanine bases. Although there is a GGGG repeat present and the last guanine in the sequence can fold back, the formation of a complete third G-quartet is not possible. This GQ has loops with a various length from two to six bases. Effect of loop length on the stability of GQs was systematically investigated by Guedin et al. 45 They investigated the G 3 T x G 3 T y G 3 T z G 3 sequence, where x,y,z, are integer numbers indicating the loop lengths. They found a clear destabilizing effect of loops larger than 3-4 bases. In presence of Na + the T m values were lower compared to those in presence of K + , and the destabilizing effect is also higher in case of K + . As mentioned previously, GQs form only in presence of cations. HepB adopts the GQ structure in presence of several stabilizing metal ions. We compare stabilizing effect of four ions and ask the question: to what extent other monovalent cations stabilize the folded form of the oligo. It is a common belief that K + or Na + ions are able to stabilize the folded form, and K + leads to a more stable structure. This was measured in case of several inter and intramolecular GQs 47-50 . Li + is believed to be too small to stabilize the structure 51 , while stabilizing effect of Rb + ion is between that of Na + and K +48 . Our results however show an interesting deviation from this common trend. Table 1 . shows the unfolding temperature (T m ) of HepB_FRET in presence of 140 mM of monovalent ions: Na + , K + Li + and Rb + . As it can be seen all the ions are able to stabilize the GQ structure at 140 mM concentration. Higher concentration increases the unfolding temperature both in the case of Na + and K + , which www.nature.com/scientificreports/ meets our expectations. The surprising result is the reversal of the effects of K + and Na + compared to the expectations. Contrary to other oligos, in case of HepB_FRET Na + seems to be a more effective stabilizer than K + . The transition temperatures of Na-stabilized GQs are higher by 2-3 °C. This difference can be observed in the whole concentration range we studied. Such a behavior is unique as far as we know. As mentioned above all the literature data report higher T m for K + stabilized GQs than for Na + stabilized ones. In our earlier works we also investigated a series of GQs 39,40 and K + was the more potent stabilizer in all cases. We obtained more than 15 °C higher T m in KCl for our above mentioned oligos. McGregor's lab also reported a 10 °C stabilization of human telomere GQ in K + compared to Na +52 . Several other papers also reported higher stability in presence of K +45, 53 . So according to the literature a higher T m value in KCl is a general phenomenon and HepB is the first example that violates this rule. Li + is known to be too small to stabilize the GQ structure, but in this case the stabilizing effect of Li + is even higher than that of Na + and K + . We hypothesize that this might be the result of the relatively large loops, which allow a more relaxed core structure. Rb + has the biggest ionic radius of 152 pm. Although its dimension is considerably larger than that of Na + (95 pm) and K + (133 pm), Rb + was also proven to be capable of stabilizing some of the GQs 54 . According to the literature its stabilizing effect is comparable with that of Na + . In case of HepB it builds a GQ with similar stability as K + , but weaker than that of Na + . Figure 2a shows the normalized fluorescence spectra of HepB_FRET at 30 °C in case of different cations used in this study. The spectra are normalized to the donor emission intensity. In this case the intensity of the acceptor emission shows the efficiency of the energy transfer. The higher the acceptor intensity is, the smaller is the donor-acceptor distance. Since the FRET efficiency depends on the distance very strongly, the observed 40% variation in the acceptor intensity indicates a very slight difference in the distance, i.e. very small distortion of the structure. The general trend of higher acceptor intensity in case of higher temperature stability can however be observed for K + , Rb + and Li + . This indicates a more compact structure in case of higher temperature stability. Na + does not fit in this trend. The energy transfer is the most efficient in the presence of Na + that suggests that the structure is the most tightly packed when it is stabilized by Na + .One can argue of course that different ions induce different structures, which can explain the alterations in the stability. It was reported that K + stabilizes the parallel, while Na + induces the antiparallel structures 49, 55 . On the other hand, in case of Htel K + stabilizes a hybrid structure, while Na + prefers the antiparallel one 56 . In order to check if there are different structures formed in case of different cations, we performed CD experiments. As the CD results (Fig. 2b) show, Na + , Li + and Rb + www.nature.com/scientificreports/ induce the same structure, while the CD spectra indicate a different conformation in case of K + . The positive band around 260 nm accompanied with a negative band at 240 nm indicates a parallel structure in case of all cations except K +57, 58 . The CD spectrum of the sample with K + shows a positive band at 290 nm, a shoulder at 270 nm and a negative band at 240 nm. This is indicative for the hybrid mixed parallel and antiparallel structure. Such hybrid structure was found in case of human telomere in presence of K +59 . From our CD experiments we can conclude, that Na + Li + and Rb + induces a parallel structure, while HepB shows a hybrid form in presence of K + . Infrared spectra of HepB in presence of different metal ions can be seen in Fig. 3 . The band at 1672 cm −1 belongs to the C 6 =O 6 vibration. The position of this band implies the presence of GQ structure 40 . The most prominent difference between the spectra is the appearance of the band at 1563 cm −1 in case of Na + , Li + and Rb + . This vibration is absent in the K + induced GQ form. The band can be assigned to C=N and C-N stretching not including the N 7 atom 40,60 . This band was present in our previous experiments on Htel but in much lower amount 40 . This intense band seems therefore to be characteristic for the presence of the parallel structure. Effect of Na + and K + concentrations on the thermal stability of the GQ was studied in the range of 100 to 250 mM. Figure 4 shows the thermal stability vs. Na + and K + concentrations. The unfolding curves can be seen in Supplementary Figs. S2 and S3 . In order to understand how and why the transition midpoint depends on the concentration of the metal ions, a short thermodynamic calculation is necessary. The concentration dependence of T m can be explained by considering the Gibbs free energy change of the unfolding: ΔG u = G ss − G GQ , where G ss and G GQ are the Gibbs free energies of the unfolded single stranded oligo, and the folded GQ respectively. ΔG u can be written as: where µ and ν are the chemical potential and the amount of substance. M denotes the free metal ion Na + or K + in the solution. The Na + and K + concentration dependence of the transition temperature was fitted with the above function, and we obtained a quite good fit. This means that the concentration dependent stabilization of the GQs can be explained by a simple thermodynamic model described above. Similar stabilization of GQs was observed by Risitano and Fox for oligos similar to Htel, but they did not fit any theoretical curve to their data 61 . The binding constants of Na + and K + were also measured. Figure 5 shows the folding of HepB as function of the concentration of Na + and K + at 30 °C (The spectra can be seen in Supplementary Fig. S4 ). The donor intensity decreases at around the dissociation constant (K d ). The lines show the fit of Eqs. (2)-(5). The dissociation constants determined from the fit are 10.2 mM and 7.1 mM for K + and Na + respectively. The smaller dissociation constant of Na + indicates a tighter binding to the GQ resulting in a higher temperature stability e.g. in T m value. This means the titration results of HepB with K + and Na + ions are consistent with the temperature stability data. Binding of ligands to HepB. GQs started attracting attention of researchers, when their presence in the telomere region was proven. Several small molecules have been developed and investigated, which can stabilize the human telomere GQ. We have chosen three of them to see whether these can be used in case of our viral GQ. We investigated the three most important representative ligands: TMPyP4, BRACO19, PhenDC3. Influence of ligands on the temperature stability. All of the ligands increased the temperature stability of HepB. This was measured by fluorescence experiments. The ligands were added in twofold excess. TmpyP4 and PhenDC3 increased the unfolding temperature of HepB_FRET by 23 °C and 24 °C respectively. On the contrary, BRACO19 had a slight stabilization effect of 5 °C. This means all the studied ligands can bind the HepB oligo, and their binding stabilizes the folded structure. Our 5 °C stabilizing effect in case of BRACO19 is in agreement with the similar data in the literature. Majee et al. found stabilization effects between 3.6 and 13.1 °C for BRACO19 and on several GQs found in the genome of the Zika virus 62 . performed at different TMPyP4 concentrations. Both absorption and fluorescence spectroscopies were used. Figure 6 shows the absorption spectra of TMPyP4 at different HepB concentrations. Binding of HepB causes a red shift of 21 nm of the Soret band (from 422 to 443 nm) and 40% hipochromicity. Similar bathochromic shift has been observed when TMPyP4 bound to other GQs. Nagesh et al. 41 observed 18 nm red shift and 60% hypochromicity for the Bcl-2 promoter GQ. The appearance of the isobestic point in the plot of absorption spectra shows the two state character of the binding. These results suggest that the stoichiometry of the binding of TMPyP4 to HepB is 1:1. The same result can be obtained from the Job plot (Fig. S5) . Unfortunately, similar experiments could not be performed for BRACO19 and PhenDC3, since their spectra overlap with that of the DNA oligo. �G u = �v(µ 00 + RT × ln((c ss × c M )/(c GQ × c 0 ))). 63 . These results compared with our values suggest, that TMPyP4 binds very strongly to HepB GQ, its affinity is similar or even higher than that of measured in case of human telomere GQ. Binding of ligand to unlabeled HepB. The above results show clearly that all the three ligands bind to the fluorescent labeled HepB. In order to prove that the fluorescence labeling does not considerably influence the binding, we used the competition 64 assay described by Luo et al. This assay can prove that the binding is not restricted to the fluorescently labeled oligo, but the non-labeled oligos can also bind the three ligands we investigated. The method was slightly modified. Instead of investigating the binding of PhenDC3 to different oligos, we measured the binding of different ligands to the same oligo. The main point of the method is the following. Thermal unfolding curves are measured for three solutions: 1. oligo labeled by a FRET pair. 2. The same solution with the ligand. 3. The previous solution together with excess of unlabeled oligo. The binding of a ligand is expected to increase the stability of the oligo, which is measured as an increased unfolding temperature (T m ). T m returns to its original value (or close to it) if the excess unlabeled oligo will bind the ligand. If only the labeled oligo binds the ligand, and it does not bind the unlabeled, the third sample shows the same T m as the second one. The S-factor defined by Luo et al. as S = (T m3 − T m1 )/(T m2 − T m1 ) is close to zero if the unlabeled oligo binds the ligand, while it is close to 1 if only the FRET labeled oligo can bind the ligand (The indices of T m in the definition of S refer to the solutions described above). Figure 7 shows the three unfolding curves in case of BRACO19 (The spectra can be seen in Supplementary Fig. S9 ). As it can be seen the unlabeled oligo captured all the ligands, and the FRET labeled oligo showed the same fluorescence intensity profile as without BRACO19. Same experiments were performed with PhenDC3 and TMPyP4 too. The unfolding curves for PhenDC3 and TMPyP4 are shown in Figs. S10 and S11. Table 3 shows the temperature increases compared to the sole HepB_FRET solution, and also the calculated S values. It can be seen, that HepB binds all the three oligos we studied. In case of TMPyP4 and PhenDC3 an interesting effect has been observed: they bind to both the labeled and unlabeled GQ with high affinity. Normally, during the preparation www.nature.com/scientificreports/ of the third solution we added the unlabeled oligo in the last step. In this case we obtained a high S value, which indicated, that the ligands stay bound to the FRET labeled GQ, and do not switch for the unlabeled one. The experiments were repeated adding the labeled HepB in the last step; we obtained no increase of T m compared to the first sample (when only the labeled HepB was present). This indicates, that both the labeled and unlabeled oligo bind TMPyP4 strongly. In case of BRACO19, the unlabeled oligo binds it better than the labeled. All the studied metal ions (Na + , K + , Li + , and Rb + ) induce GQ structure in the HepB sequence in the genome of the Hepatitis B virus. HepB shows a hybrid form in presence of K + , but all the other ions induce parallel structure. Higher concentrations of metal ions increase the unfolding temperature. Study of three ligands designed for GQ targeting (TMPyP4, BRACO19 and PhenDC3) showed clearly their binding to HepB and to its fluorescently labeled variant. The binding to the unlabeled HepB was proven by the competitive assay. Additionally, the above results show an increased stability of the ligand bound GQs. The stabilization was higher than 20 °C for TMPyP4 and PhenDC3, while it was considerably lower for BRACO19. Binding of the TMPyP4 and PhendDC3 to this viral GQ might have important medical relevance. These ligands, which were developed for cancer treatment, could have a potential role in the fight against the HepB virus. This hypothesis should however be confirmed by several further studies, but we believe this might be a promising new perspective. 64 ). T m2 − T m1 = increase of HepB_FRET when adding the ligand in twofold excess. T m3 − T m1 = increase of the transition midpoint when both the ligand and the unlabeled HepB were added. S = (T m3 − T m1 )/(T m2 − T m1 ). First the ligand then the unlabeled HepB was added except of cases marked by *, where the unlabeled oligo was added first. T m2 − T m1 (°C) T m3 − T m1 (°C) S Tetraplex DNA and its interacting proteins DNA: Sequence, topology and structure Characterization of parallel and antiparallel G-tetraplex structures by vibrational spectroscopy Stability prediction of canonical and non-canonical structures of nucleic acids in various molecular environments and cells Effect of molecular crowding on the stability of RNA G-quadruplexes with various numbers of quartets and lengths of loops Exploring the effects of cosolutes and crowding on the volumetric and kinetic profile of the conformational dynamics of a poly dA loop DNA hairpin: A single-molecule FRET study Effect of pressure on RNA G-quadruplex structures The deep sea osmolyte trimethylamine N-oxide and macromolecular crowders rescue the antiparallel conformation of the human telomeric G-quadruplex from urea and pressure stress Duplex-tetraplex equilibria in guanine-and cytosine-rich DNA Thermodynamics of the G-quadruplex formation of modified human telomeric sequences Unraveling the binding characteristics of small ligands to telomeric DNA by pressure modulation Peroxidase-mimicking DNAzymes for biosensing applications: A review The structure of telomeric DNA Effect of pressure on the stability of G-quadruplex DNA: Thermodynamics under crowding conditions Targeting G-quadruplexes in gene promoters: A novel anticancer strategy? Direct evidence for a G-quadruplex in a promoter region and its targeting with a small molecule to repress c-MYC transcription Structure of an unprecedented G-quadruplex scaffold in the human c-kit promoter Facilitation of a structural transition in the polypurine/polypyrimidine tract within the proximal promoter region of the human VEGF gene by the presence of potassium and G-quadruplex-interactive agents Replication of G quadruplex DNA Platinum(II) and palladium(II) complexes of tridentate hydrazone-based ligands as selective guanine quadruplex binders Dual functional dinuclear platinum complex with selective reactivity towards c-myc G-quadruplex Quadruplex nucleic acids as targets for anticancer therapeutics Targeting glioma stem cells in vivo by a G-quadruplex-stabilizing synthetic macrocyclic hexaoxazole The different biological effects of telomestatin and TMPyP4 can be attributed to their selectivity for interaction with intramolecular or intermolecular G-quadruplex structures The G-quadruplex-interactive molecule BRACO-19 inhibits tumor growth, consistent with telomere targeting and interference with telomerase function Genetic instability triggered by G-quadruplex interacting Phen-DC compounds in Saccharomyces cerevisiae More than just a kink in microbial genomes G-quadruplexes in pathogens: A common route to virulence control? G-quadruplex DNA motifs in the malaria parasite Plasmodium falciparum and their potential as novel antimalarial drug targets Multiple G-quartet structures in pre-edited mRNAs suggest evolutionary driving force for RNA editing in trypanosomes G-quadruplex forming sequences in the genome of all known human viruses: A comprehensive guide G-quadruplexes in viruses: Function and potential therapeutic applications Targets and tools in antiviral therapy RNA synthesis is modulated by G-quadruplex formation in Hepatitis C virus negative RNA strand Hepatitis B virus infection Hepatitis B virus infection A G-quadruplex motif in an envelope gene promoter regulates transcription and virion secretion in HBV genotype B Negative volume changes of human G-quadruplexes at unfolding Self-crowding influences the temperature-pressure stability of the human telomere G-quadruplex Studies on the site and mode of TMPyP4 interactions with Bcl-2 promoter sequence G-Quadruplexes Smoothing and differentiation of data by simplified least squares procedures How precise are the positions of computer-determined peaks? The effect of osmolytes and small molecule on Quadruplex-WC duplex equilibrium: A fluorescence resonance energy transfer study How long is too long? Effects of loop size on G-quadruplex stability A Sequence-independent analysis of the loop length dependence of intramolecular RNA G-quadruplex stability and topology Cation-dependent transition between the quadruplex and Watson-Crick Hairpin forms of D(Cgcg3gcg) Parallel and antiparallel G-DNA structures from a complex telomeric sequence The diverse structural landscape of quadruplexes Metal cations in G-quadruplex folding and stability Quadruplex Nucleic Acids Preface The role of loops and cation on the volume of unfolding of G-quadruplexes related to HTel Effect of loops and G-quartets on the stability of RNA G-quadruplexes Effect of rubidium and cesium ions on the dimeric quaduplex formed by the Oxytricha nova telomeric repeat oligonucleotide d(GGG GTT TTG GGG ) Effects of metal ions and cosolutes on G-quadruplex topology Assessment of human telomeric G-quadruplex structures using surface-enhanced Raman spectroscopy Role of alkali metal ions in G-quadruplex nucleic acid structure and stability Topological characterization of nucleic acid g-quadruplexes by UV absorption and circular dichroism Human telomeric sequence forms a hybrid-type intramolecular G-quadruplex structure with mixed parallel/ antiparallel strands in potassium solution Tetraplex structure formation in the thrombin-binding DNA aptamer by metal cations measured by vibrational spectroscopy Stability of intramolecular DNA quadruplexes: Comparison with DNA duplexes Inhibition of Zika virus replication by G-quadruplex-binding ligands Evidence for the binding mode of porphyrins to G-quadruplex DNA A fluorescence melting competition assay for studying G4 structures in vitro The authors thank Tünde Juhász for her help in the CD experiments, Levente Herényi for discussions of the thermodynamic model and Anna Grád for carefully reading the manuscript. This work was supported by the National Research, Development and Innovation Office of Hungary [NKFI K-124697]. Open Access funding provided by Semmelweis University. The authors declare no competing interests. The online version contains supplementary material available at https:// doi. org/ 10. 1038/ s41598-021-02689-y.Correspondence and requests for materials should be addressed to L.S. Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.