Next Article in Journal
The Incidence and Associated Risk Factors of Contrast-Induced Nephropathy after Contrast-Enhanced Computed Tomography in the Emergency Setting: A Systematic Review
Next Article in Special Issue
Characterization of Aurintricarboxylic Acid (ATA) Interactions with Plasma Transporter Protein and SARS-CoV-2 Viral Targets: Correlation of Functional Activity and Binding Energetics
Previous Article in Journal
Ecdysteroids as Potent Enzyme Inhibitors and Verification of Their Activity Using In Vitro and In Silico Docking Studies
Previous Article in Special Issue
The Pressure Dependence of the Stability of the G-quadruplex Formed by d(TGGGGT)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Thermally Induced Transitions of d(G4T4G3) Quadruplexes Can Be Described as Kinetically Driven Processes

1
Biotechnical Faculty, University of Ljubljana, Jamnikarjeva 101, 1000 Ljubljana, Slovenia
2
Faculty of Chemistry and Chemical Technology, University of Ljubljana, Vecna Pot 113, 1000 Ljubljana, Slovenia
*
Authors to whom correspondence should be addressed.
Life 2022, 12(6), 825; https://doi.org/10.3390/life12060825
Submission received: 6 April 2022 / Revised: 24 May 2022 / Accepted: 27 May 2022 / Published: 1 June 2022

Abstract

:
DNA sequences that are rich in guanines and can form four-stranded structures are called G-quadruplexes. Due to the growing evidence that they may play an important role in several key biological processes, the G-quadruplexes have captured the interest of several researchers. G-quadruplexes may form in the presence of different metal cations as polymorphic structures formed in kinetically governed processes. Here we investigate a complex polymorphism of d(G4T4G3) quadruplexes at different K+ concentrations. We show that population size of different d(G4T4G3) quadruplex conformations can be manipulated by cooling rate and/or K+ concentration. We use a kinetic model to describe data obtained from DSC, CD and UV spectroscopy and PAGE experiments. Our model is able to describe the observed thermally induced conformational transitions of d(G4T4G3) quadruplexes at different K+ concentrations.

1. Introduction

Understanding the relationship between the structure of biological macromolecules and the thermodynamic or kinetic properties that dictate stability and binding with other molecules remains one of the most important problems in biochemistry and biotechnology [1]. Such knowledge is crucial to understanding biological processes and to design more efficient pharmaceutical ligands. The laws of thermodynamics and kinetics do not tell us directly about molecular structures or mechanisms, although we can use the results of thermodynamic property measurements to help us interpret molecular mechanisms of a system [2].
Over the past couple of decades, DNA sequences that are both rich in guanines and can form four-stranded structures, called G-quadruplexes, have captured the interest of several major laboratories all over the world. This is mainly due to the growing evidence that G-quadruplexes may play an important role in several key biological processes [3,4]. Even though much progress has been made in understanding kinetic and thermodynamic factors responsible for the structural interconversion of several quadruplex structures, there are still many factors that govern the formation of G-quadruplexes and their physico-chemical properties that are poorly understood.
Four guanine bases linked by Hoogsteen type hydrogen bonds form a cyclic coplanar G-tetrads which can stack on top of another tetrad to form G-quadruplex structures [5,6]. These structures are additionally stabilized by metal cations that are selectively bound in the central cavity between the G-quartets. Cation coordination is essential for the stabilization of G-quadruplexes, and monovalent and divalent ions have been shown to influence the structure and stability of G-quadruplexes. By compiling a number of studies, one can estimate that G-quadruplex stabilization follows the general trend: Sr2+ > K+ > Ca2+ > NH4+, Na+, Rb+ > Mg2+ > Li+ ≥ Cs+ [7].
Changing the type or concentration of cations can induce the formation of new long-lived G-quadruplex conformations with identical nucleotide sequence. For instance, oligonucleotide d(G2AG2AG) forms a bimolecular G-quadruplex at low Na+ concentration and a tetramolecular G-quadruplex at higher Na+ concentration [8]. Furthermore, d(G3T4G3) forms only hairpin dimer G-quadruplex structures in the presence of Na+ ions, while in the presence of K+ ions both hairpin dimer and linear, four-stranded G-quadruplex structures were observed [9]. This phenomenon is known as G4 polymorphism [10,11,12] and can even lead to concurrent folding isomers in heterogenous ensembles [13,14,15]. The polymorphism of GQ structures is determined by strand orientation (i.e., parallel vs. antiparallel), the conformation of the glycosidic bonds, and the loop topology (e.g., lateral loops vs. diagonal loops). These variations in the structure might be due to several reasons, such as the cationic coordination, π-stacking interactions, hydrogen bonding and hydrophobic effects [16]. Different mechanisms have been proposed to describe folding/unfolding pathways of G-quadruplexes. Some of them are sequential and do not consider side-reactions [17] while others consist of several parallel pathways [18,19].
The behavior of G-quadruplexes in the presence of potassium ions is of particular interest due to the high intracellular concentrations. Chaires et al. have suggested, based on NMR, fluorescence, CD experiments, and molecular modeling, that model telomeric human DNA, 5′-AGGGTTAGGGTTAGGGTTAGGG-3′ (Tel22) in K+ solutions appears as an equilibrium mixture of two (3 + 1) hybrid-type G-quadruplex structures, Hybrid-1 and Hybrid-2 [20]. This was later confirmed by Bončina et al., who successfully described the unfolding/folding mechanism of Tel22 and showed that each conformational transition depends on K+ concentration [11]. Grün et al. were studying the cMYC gene-promoter region and used a kinetic model to describe a novel, noncanonical folding isomerism of G-quadruplexes with more than four G-rich tracts [21]. Grey et al. utilized a variety of spectroscopies to follow the kinetics of human telomeric DNA sequences folding and unfolding, and identified previously unreported slow kinetic steps [17]. Kinetics of potassium binding to telomeric (hybrid) and c-myc (parallel) G-quadruplexes have been investigated with electrospray mass spectrometry where authors have reported that misfolded states can linger for several minutes [18]. Different folding timescales of DNA and RNA G-quadruplexes have been reported by Zhang and Balasubramanian, who discussed their findings in relation to a biological timescale [22]. The folding rate of quadruplexes and folding timescale were also investigated by Nguyen et al., who showed that a stem loop accelerates the folding of quadruplexes significantly as compared to a non-structured loop [23]. One can conclude that the formation of G-quadruplexes takes place through different mechanisms on various timescales, making the folding rate of a G-quadruplex and the effect of the environment and loop properties very important parameters to investigate.
Our laboratory has used kinetic models to describe the polymorphism of d(G4T4G3) quadruplexes in Na+ and K+ solutions [24,25]. To test whether this approach can be used under different environmental conditions, we followed in this work the thermally induced folding/unfolding transitions d(G4T4G3) quadruplexes at different K+ concentrations, and tried to interpret them in terms of a kinetically governed coexistence of several folded structures and their unfolded forms. We used a variety of experimental techniques (UV and CD spectroscopy, DSC, gel electrophoresis) in order to improve the significance of the proposed model and calculated parameters. We demonstrate that the d(G4T4G3) sequence can fold into different quadruplex structures and that the mechanism of thermally induced folding/unfolding depends on the cooling/heating rate and potassium concentration.

2. Materials and Methods

2.1. Sample Preparation and UV Spectroscopy

The d(G4T4G3) oligonucleotide was obtained HPLC pure from Invitrogen Co., (Carlsbad, CA, USA) and Midland Co. (Midland, CO, USA)., and used without further purification. Its concentrations in buffer solutions were determined at 25 °C spectrophotometrically using for the extinction coefficient of its single stranded form at 25 °C the value of ε260 = 105,100 M−1cm−1 estimated from the nearest-neighbor data of Cantor et al. [26]. The quadruplex concentration used in this study were 0.82, 0.60 and 0.52 mM in single strands in 25, 35 and 50 mM K+ solutions, respectively. The buffer used in all experiments consisted of 10 mM K-cacodylic buffer, 1 mM EDTA and different concentrations of KCl (15, 25 and 45 mM) at pH = 6.9. The starting samples were first unfolded into single-stranded form by heating in an outer thermostat at 95 °C for 5 min, cooled down to 4 °C at the cooling rates of 0.05 or 1.0 °C/min to form quadruplex structures and then used in the UV, CD, DSC and PAGE experiments.

2.2. CD Spectroscopy Melting Experiments

CD spectra of d(G4T4G3) quadruplexes in K+ solutions were measured as a function of temperature in an Jasco J-1500 CD Spectropolarimeter. Ellipticity, Θ, was measured between 5 and 95 °C in the temperature intervals of 3 °C at the average heating rate of 0.5 °C/min. CD spectra of samples (0.82, 0.60 and 0.52 mM in single strands) prepared at cooling rates of 0.05 or 1.0 °C/min, corrected for the corresponding buffer contribution, were collected between 215 and 320 nm in a 0.25 mm cuvette at 60 nm/min, with a signal averaging time of 2 s and a 5 nm bandwidth. Melting curves were obtained by plotting ellipticity at λ = 290 nm vs. temperature.

2.3. Gel Electrophoresis

G-quadruplex structures formed upon the cooling of single-stranded DNA at different rates of 0.05 or 1.0 °C/min and at different potassium concentrations were studied by non-denaturing PAGE performed on 20% polyacrylamide gels supplemented with 25, 35 and 55 mM KCl. G-quadruplex samples were loaded on gels and the electrophoreses were run at 5 °C (5 h), 20 °C (3.5 h) and 35 °C (1.7 h), all at 10 V/cm (I = 300 mA). Bands in the gels were followed by UV shadowing at λ = 254 nm. To facilitate comparisons between the bands observed with different samples, the single-stranded d(5′-AGAAGAAAAGA-3′) and d(5′-TCTTTTCTTCT-3′) and double-stranded (5′-AGAAGAAAAGA-3 mixed with 5′-TCTTTTCTTCT-3′) 11-mere control oligonucleotides were used.
After UV shadowing, experiment gels were submerged in ethidium bromide solution and bands were recorded in fluorescent mode at an excitation wavelength of 290 nm.

2.4. Differential Scanning Calorimetry (DSC)

DSC experiments were performed using the Nano DSC III instrument (Calorimetry Sciences Corp., Lindon, UT, USA) and the Nano DSC instrument (TA Instruments, New Castle, DE, USA) on samples prepared by cooling in an outer thermostat. Samples were first heated to 95 °C and cooled down to four °C at the selected rate (0.05 or 1.0 °C/min). Samples were then loaded to DSC and the first DSC melting scans were recorded at the selected heating rate (0.5, 1.0 and 2.0 °C/min). Next, an annealing scan was recorded at a cooling rate of 1.0 °C, followed by another melting scan at the selected heating rate. The corresponding baseline (buffer-buffer) scans were subtracted from the unfolding/folding scans prior to their normalization and analysis. The total enthalpy of unfolding or folding, Δ H tot , was obtained from the measured DSC thermograms as the area under the Δ c p = c ¯ p , 2 c ¯ p , S versus T curve, where Δ c p is the measured heat capacity c ¯ p , 2 corrected for the baseline and normalized to 1 mole of quadruplex in single strands, and c ¯ p , S is the corresponding partial molar heat capacity of the unfolded single stranded state extrapolated from high temperatures over the whole measured temperature interval. For samples prepared either at the cooling rate of 0.05 or 1.0 °C/min, the measured DSC heating and cooling curves were highly reproducible at all measured heating and cooling rates above ~10 °C.

2.5. Model Analysis of Structural Transitions

We attempted to interpret the observed DSC data in terms of the already described complex model (Figure 1) [24], which considers the thermally induced folding/unfolding transitions of d(G4T4G3) quadruplexes in the presence of K+ ions as a global kinetic transition process that involves one parallel tetramolecular quadruplex structure (T4), three bimolecular quadruplexes that possibly exhibit characteristics of parallel, anti-parallel and/or hybrid structures (A2, B2 and C2), and the corresponding single strands (S).
In order to describe the DSC data, we used the previously reported model function [24]:
Δ c p = c ¯ p , 2 c ¯ p , S = d α A d T Δ H A d α B d T Δ H B d α C d T Δ H C d α T d T Δ H T
in which at any T in the measured temperature interval c ¯ p , 2 is the measured partial molar heat capacity of the solute (DNA) expressed per moles of single strands, c ¯ p , S is the corresponding measured heat capacity of the unfolded quadruplex form occurring at high temperatures extrapolated to T, and α A , α B , α C and α T are the corresponding fractions of quadruplex species present in the solution. The d α d T terms needed to calculate Δ c p are obtained by taking into account the rates of reactions predicted by the model and the heating or cooling rate, r = d T d t , at which the DSC experiment is performed. The rates of reactions can be expressed using the Arrhenius relation k i j = A i j e ( E i j / R T ) = e A i j e ( E i j / R T ) = e ( A i j E i j / R T ) , where k i j represents rate constant, A i j = e A i j frequency factor and E i j activation energy for corresponding i j step. Unit of e A i j is s−1 for unfolding and s−1 M−1 for folding transitions of A2, B2 and C2 (s−1 M−3 for folding transition of T4 from S, i.e. e A ST ) and unit of E i j is cal/mol.
For a given total concentration of DNA in the single-stranded form,   c tot (in unit mol/L), expressed as c tot = c S + 2 c A + 2 c B + 2 c C + 4 c T , the d α / d T terms needed to calculate Δ c p are obtained by taking into account the rates of reactions predicted by the model (Figure 1):
d α A d T = 1 r [ ( k AB + k AS + k AC ) α A + C tot k AT α A 2 k BA α B k TA α T k CA α C 2 k SA C tot ( 1 α A α B α C α T ) 2 ]
d α B d T = 1 r [ ( k BS + k BC + k BA ) α B + C tot k BT α B 2 k AB α A k TB α T k CB α C 2 k SB C tot ( 1 α A α B α C α T ) 2 ]
d α C d T = 1 r [ ( k CB + k CS + k CA ) α C + C tot k CT α C 2 k BC α B k TC α T k AC α A 2 k SC C tot ( 1 α A α B α C α T ) 2 ]
d α T d T = 1 r [ ( k TA + k TB + k TC + k TS ) α T C tot k AT α A 2 C tot k BT α B 2 C tot k CT α C 2 4 k ST C tot 3 ( 1 α A α B α C α T ) 4 ]
When fitting the kinetic model (Figure 1) to our experimental data, all parameters were allowed to be used. We found out that not all of them were necessary to obtain an acceptable fit of model functions to experimental data. We neglected these parameters (set k ij to 0), thus simplifying the system of equations and increasing the significance of the remaining parameters:
d α A d T = 1 r [ ( k AB + k AS ) α A k BA α B ]
d α B d T = 1 r [ ( k BS + k BC + k BA ) α B + C tot k BT α B 2 k AB α A 2 k SB C tot ( 1 α A α B α C α T ) 2 ]
d α C d T = 1 r [ k CS α C k BC α B ]
d α T d T = 1 r [ k TS α T C tot k BT α B 2 ]
To solve this system of differential equations for a given set of adjustable parameters at each measured heating and cooling rate, the Cash-Karp adaptive step-size controlled Runge-Kutta method was employed [27], and the obtained solutions ( α i and d α i d T ) were used to calculate the corresponding model function. The calculation started at a high temperature and a set cooling rate. The initial composition was α S = 1 . When the low temperature was reached, the obtained composition was a starting point for heating. The “best fit” adjustable parameters were calculated from global fitting model functions to the experimental DSC curves using the non-linear minimization of the corresponding χ 2 function. The values of adjustable parameters at the global minimum of χ 2 are considered to be the best descriptors of the experimental Δ c p vs. T curve, and therefore are used to characterize the kinetics and thermodynamics for all steps in the suggested model mechanism.

3. Results and Discussion

3.1. Observing Polymorphism of d(G4T4G3) Quadruplexes in K+ Solution

We measured the CD spectra of d(G4T4G3) samples prepared in 25, 35 and 50 mM K+ solutions at the cooling rate of 0.05 °C/min (Figure 2a), and at the cooling rate of 1.0 °C/min (Figure 2b). The CD spectra is often used to interpret the structure of G-quadruplexes. CD-spectra with positive peak at ~290 nm and negative peak at ~263 nm are characteristic of the anti-parallel type quadruplex structures. On the other hand, the parallel type structures exhibit CD spectra with a positive peak at ~263 nm and negative peak at ~240 nm [28,29]. Hybrid quadruplex structures may also be involved in thermally induced folding/unfolding transitions [11,30,31]. The comparison of three spectra in Figure 2a shows that the peak at 263 nm increases, and the peak at 290 nm decreases with increasing K+ concentration. The same behavior can be observed when samples are prepared at the cooling rate of 1.0 °C/min (Figure 2b). Thus, we can conclude that increasing K+ concentration favors the formation of a structure with spectral properties similar to parallel type quadruplexes. Interestingly, when comparing the CD spectra of samples prepared at different cooling rates, we can observe that a slower cooling rate results in increased peak at 263 nm and a decreased peak at 290 nm. These results show that besides changing K+ concentration, different cooling rates can be used to manipulate the structure of d(G4T4G3) quadruplexes during thermally induced folding. The cooling rate dependence of structural properties of d(G4T4G3) quadruplexes point to the conclusion that thermally induced folding/unfolding transitions of d(G4T4G3) quadruplexes are a kinetically driven process. Furthermore, because of the increased peak at 263 nm after a slow cooling rate, we can assume that a slower cooling rate also favors the formation of a structure with spectral properties similar to parallel-type quadruplexes.
The CD spectra suggest that at least two structures with antiparallel and parallel spectral properties (for reasons of clarity we will call them “antiparallel” and “parallel” structures) coexist in solution after cooling, and that the ratio between antiparallel and parallel structures depends on the potassium concentration and a cooling rate. To confirm these results, we performed PAGE experiments where slow and fast cooled d(G4T4G3) samples in 25 mM and 50 mM K+ solutions were run on the gels at 5 °C, the same temperature at which CD-spectra were recorded (Figure 3). Control oligonucleotide markers (single-stranded d(5′-TCTTTTCTTCT-3′), single-stranded d(5′-AGAAGAAAAGA-3′) and double-stranded (5′-AGAAGAAAAGA-3′ mixed with 5′-TCTTTTCTTCT-3′)) were used to ascribe the size of the molecules to the positions of the bands, and to compare different gels. The gel experiments show that all samples exhibit at least three bands where two fast migrating bands are very close to each other, indicating the presence of at least two structures similar in size and charge. Furthermore, the samples prepared at the cooling rate of 0.05 °C/min (2nd lane) exhibit a more pronounced band, positioned just behind the fastest migrating band, which in accordance with CD-spectra probably represents a parallel structure. This observation can be further confirmed by comparing gels at different K+ concentrations since the slower of the two fast migrating bands is more pronounced at 50 mM K+ concentration. Besides the two fast migrating bands, one slow migrating band can be observed at both K+ concentrations and cooling rates, but the intensity is less pronounced in samples prepared at cooling rate of 1.0 °C/min (1st lane). Interestingly, staining the gels with ethidium bromide results in high intensity fluorescence signals for the slowest migration band and double stranded marker, whereas very weak signals can be observed for the fastest migrating bands. Based on the fluorescence intensity, we can conclude that ethidium bromide interacts with the structure in the slowest migrating band in a similar way to the double stranded marker, and that interaction with quadruplexes is weaker compared to duplex DNA as already observed before [32]. The position and difference in fluorescence intensity between fast and slow migrating bands indicates the significant differences of complexes, and we believe that the slower migrating band corresponds to some highly ordered structure, for example the parallel tetramolecular d(G4T4G3)4 quadruplex.

3.2. Thermally Induced Structural Transitions and Model Analysis

To observe the thermal behavior of d(G4T4G3) starting samples, we have conducted a series of DSC and CD experiments. DSC heating thermograms of d(G4T4G3) samples measured in 25, 35, and 50 mM K+ solutions and prepared at the cooling rate of either 0.05 or 1.0 °C/min exhibit at least three peaks suggesting that the melting processes consist of several conformational transitions (Figure 4). Samples prepared at the cooling rate of 0.05 °C/min exhibit peaks at ~20 °C, ~45 °C and ~65 °C, suggesting the coexistence of at least three quadruplex structures, which is in line with PAGE experiments. Increasing the K+ concentration decreases the intensity of the second peak and increases the intensity of the third. Together with the data from CD experiments, this suggests that the third peak corresponds to the unfolding of a parallel structure and that the parallel structure is thermally more stable than the antiparallel. This can be further proved by PAGE experiments of samples prepared at a slow cooling rate at 50 mM K+ concentration (Figure 3b) where the increased intensity of the slower of the two fast migrating bands can be observed. Samples prepared at the cooling rate of 1.0 °C/min also exhibit three peaks, but the third peak is not as expressed as in the case of a slower cooling rate. These results, together with the CD (Figure 2) and PAGE (Figure 3) experiments, strongly suggest that the thermal unfolding transition of d(G4T4G3) in the presence of K+ ions may be considered as a combination of several kinetically governed steps.
The total enthalpy of unfolding, Δ H tot , was calculated from the experimentally obtained DSC thermograms (total area under the Δ C P versus T curve). Analysis of the total area under measured heating DSC thermograms resulted in very similar overall enthalpies of unfolding, ΔHtot ≈ 56 kcal/mol of double stranded quadruplex, regardless of the cooling rate or K+ concentration used. From these ΔHtot values, the enthalpy of quadruplex unfolding is estimated to be about 19 kcal/mol of G-quartets, which agrees well with the literature values reported for G-quartet formation in the presence of K+ ions [33,34,35]. After the first melting scan was recorded at the selected heating rate (0.5, 1.0 and 2.0 °C/min), the sample was left in DSC and a cooling scan was recorded at the rate of 1.0 °C/min, followed by another melting scan at the selected heating rate. Almost no differences were observed when comparing the second melting scans to the melting scan at corresponding heating rates obtained from samples prepared in an outer thermostat at the cooling rate of 1 °C/min. This confirms the kinetic nature of the folding/unfolding transition of d(G4T4G3) in the presence of K+ ions.
Interestingly, DSC thermograms show transition at low temperatures. Even though the reproducibility and reliability of the measured DSC peaks between 4 °C and ~20 °C is rather poor, PAGE experiments confirmed that there is an additional structure present at low temperatures that disappears at temperatures over 20 °C (Figure S2). The position and difference in fluorescence intensity between fast and slow migrating bands (Figure 3) indicates significant differences in structures of complexes, and we believe that a slower migrating band corresponds to some highly ordered structure, for example, the parallel tetramolecular d(G4T4G3)4 quadruplex.
The results of DSC, PAGE and CD experiments performed on d(G4T4G3) quadruplexes in K+ solutions reveal that their thermally induced folding/unfolding transitions are kinetically governed and include the participation of at least three structures. A model function (1) was used to describe thermally induced folding/unfolding transitions of d(G4T4G3) quadruplexes in the presence of K+ ions at different concentrations. Figure 4 shows the best fit of model function to the DSC data. To further test the validity of our model, we have collected and fitted data at several heating rates (0.5, 1.0 and 2.0 °C/min) for DSC experiments. By increasing the heating rate, we can observe shifting of the DSC peaks to higher temperatures (Figure 5), which is another indication that conformational transitions of d(G4T4G3) in the presence of K+ ions are kinetically driven processes. The model function (1) with the same set of parameters (Table 1) was used to describe experimental data at all heating rates at selected K+ concentration.
By comparing the parameters in Table 1 at different potassium concentrations, we can observe an increase of parameter ABC, which is connected to the frequency factor for the transition from structure B2 to structure C2 by equation A BC = e A BC . Also, we can observe a decrease of parameter ACS, which is connected to frequency factor for transition from structure B2 to structure C2 by equation A BC = e A BC . Frequency factors are proportional to the rate constant, kij, which means that increase in potassium concentration will increase the population of structure C2. This behavior of model-based parameters correlates nicely with the experimental data which show higher populations of parallel quadruplex at CK+ = 50 mM after slow and fast rates of cooling. Calculated population distributions of bimolecular structures A2 and C2 at slower and at higher cooling rates at different potassium concentrations are shown in Figure 6. By comparing the data from Figure 4, Figure 5 and Figure 6 we can see that increase of the third peak in DSC thermograms correlates to the increase of the ratio between the calculated fractions of bimolecular quadruplex species C2 and A2. This means that the position of the third peak in DSC heating thermograms, which increases in intensity when slow cooling or high concentrations of potassium ions are used, corresponds to the unfolding of structure C2.
The correlation between fitting parameters was calculated for all potassium concentrations (Tables S1–S3). Most of the parameters show weak correlation but stronger correlation is observed between parameters Aij and Eij. Since the choice of baseline is the largest source of errors accompanying the DSC measurements, the errors were estimated by changing the baseline and plotting “top and bottom” DSC curves which served as confidence limits for our experiment (Figure S1). Applying the calculated errors simultaneously to highly correlated parameters resulted in poorer agreement between the model function and the experiment (data not shown). This suggests that within the estimated error margins the observed higher correlation between parameters Aij and Eij has no significant effect on their reliability.
Fitting the model function to experimental DSC data yields temperature dependence of population of predicted quadruplex structures involved in the measured folding/unfolding process. In order to successfully interpret the observed DSC data, we had to include one additional bimolecular quadruplex structure B2 in our model (Figure 1). The calculated population of quadruplex structure B2 is negligibly small at all temperatures, suggesting that B2 acts as an unstable intermediate. Even though the population of B2 is negligibly small, its presence in our model is necessary to fit the entire cycle of DSC curves at selected potassium concentration (two different cooling rates and three different heating rates) with a single set of parameters. Efforts were made to fit the experimental data without quadruplex structure B2, but this less complex model was only able to satisfactorily describe either the melting or cooling curves.
The calculated populations of other predicted structures can be used to construct a CD-melting curve and compare it to the experimental CD-curve. Figure 7 shows normalized CD melting curves of d(G4T4G3) quadruplexes obtained from experiments and the CD-melting curves calculated from the concentrations of predicted structures in the solution. Good agreement between the calculated and experimental CD-melting curves can be observed at 25 and 35 mM K+ concentrations (Figure 7a,b). However, the agreement of model melting curve with the experimental melting curve at 50 mM K+ concentration is not so good but still acceptable (Figure 7c). We attribute this discrepancy to the error in estimated heating rate during CD-measurements.
Even though our proposed model is able to satisfactory describe experimental data, it is worth mentioning that the rate of quadruplex transitions can be slower than the DSC timescale [21,36]. Our DSC experiments were performed at different scanning rates, but other than that, we did not consider this issue in experimental setup and kinetic modeling. For example, if the transition between structures A2 and C2 is very slow, DSC would be unable to record it. Although the proposed kinetic model (Figure 1) allowed multiple pathways, only a minimum number of parameters was used to obtain a good agreement between model prediction and experimental results. The elimination of parameters might be interpreted as the nonexistence of corresponding transitions whereas in reality, the transitions were too slow to be recorded by DSC in dynamic mode or the change in enthalpy was close to zero. Additional experiments should be performed to quantify the real number of species in mixture and to define which (if any) are kinetically trapped intermediates. Fully matched quadruplex assemblies are the most stable structures but mismatched species can form at a faster rate and have long lifetimes [37,38,39]. For example, if we use data from Table 1 at C K + = 25   mM   and T = 281 K, we can calculate the formation and dissociation rates of structures A2 ( k BA = 19.6   s 1 , k AS = 2.6 · Δ 10 4   s 1 ) and C2 ( k BC = 0.05   s 1 , k CS = 5.0 · Δ 10 6   s 1 ). Since k BA > k BC and k AS < k CS , we may conclude that although folding of structure A2 is faster, C2 is thermodynamically more stable. When the concentration of potassium is increased ( C K + = 50   mM ), the folding rate of A2 is decreased ( k BA = 10.7   s 1 ) while the folding rate of C2 is increased ( k BC = 0.66   s 1 ). At the same time, the increased potassium concentration further stabilizes C2 by decreasing the dissociation rate ( k CS = 9.2 × 10 8   s 1 ). One might jump to a conclusion that A2 is kinetically trapped and C2 is a thermodynamically stable structure, but further experiments on a longer timescale are needed to confirm this hypothesis.

4. Conclusions

The results of DSC, PAGE and CD experiments performed on d(G4T4G3) quadruplexes at different K+ concentrations suggest a complex thermally induced folding/unfolding mechanism that includes the participation of at least three quadruplex structures. To describe the effects of cooling/heating rates and potassium concentration on thermally induced formation/disruption of d(G4T4G3) quadruplex structures, we developed a kinetic model that involves one parallel tetramolecular quadruplex structure, three bimolecular quadruplexes and the corresponding single strands. The model was able to describe an interesting behavior of the d(G4T4G3) quadruplex where the distribution of structures with different spectral and gel mobility properties could be manipulated by cooling rate and the concentration of potassium ions. We would like to emphasize that following thermally induced folding/unfolding transitions by DSC is a dynamic measurement. There could be additional transitions present in our system but if their rate is slow relative to the DSC timescale they remained undetected.
We believe that our research provides new insights into how conformational changes of highly ordered G-quadruplexes can be manipulated by the cooling rate of thermally induced folding and concentration of potassium ions. Also, we have shown how to use a kinetic model to describe the thermally induced folding/unfolding processes of a d(G4T4G3)4 quadruplex and obtain kinetic and thermodynamic parameters for individual elementary steps. It would be interesting to include the concentration of potassium ions as the fitting parameter in our model, although additional experiments would have to be performed at several other K+ concentrations. Furthermore, additional experiments on a longer timescale should be performed to assess whether the system is at equilibrium and not kinetically trapped. However, this remains the focus of our future work.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/life12060825/s1, Figure S1. In order to estimate experimental error, two baselines were used to subtract DSC data. First baseline yielded largest overall area under thermogram (empty circles) and the second baseline yielded smallest overall area under thermogram (full circles). The Δcp vs. T curves were measured at the heating rate of 1.0 °C/min for samples prepared at the cooling rate of 0.05 °C/min; (a) in 25 mM K+ solution (b) in 35 mM K+ solution; (c) in 50 mM K+ solution, Figure S2. Non denaturing PAGE (20%) of d(G4T4G3) quadruplexes prepared at the cooling rate of either 1.0 °C/min (1st lane) or 0.05 °C/min (2nd lane) in 50 mM K+ solution and followed by UV shadowing (λ = 254 nm) and Fluorescence (λ = 590 nm): (a) d(G4T4G3) quadruplexes measured at constant temperature of 5 °C; (b) d(G4T4G3) quadruplexes prepared at constant temperature of 20 °C. (c) d(G4T4G3) quadruplexes prepared at constant temperature of 35 °C. As control oligonucleotide double stranded d(5’-AGAAGAAAAGA-3’; 5’-TCTTTTCTTCT-3’) (3rd lane), single stranded d(5’-TCTTTTCTTCT-3’) (4th lane) and single stranded d(5’-AGAAGAAAAGA-3’) (5th lane) markers were used. Table S1: Correlation matrix corresponding to the model analysis of the DSC melting and cooling thermograms for concentration K+ ions equal to 25 mM, Table S2: Correlation matrix corresponding to the model analysis of the DSC melting and cooling thermograms for concentration K+ ions equal to 35 mM, Table S3: Correlation matrix corresponding to the model analysis of the DSC melting and cooling thermograms for concentration K+ ions equal to 50 mM.

Author Contributions

Data curation, I.P.; Formal analysis, I.P. and T.U.; Funding acquisition, N.P.U.; Methodology, I.P. and T.U.; Software, I.P. and T.U.; Writing—original draft, I.P.; Writing—review & editing, T.U. and N.P.U. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Slovenian Research Agency, grant numbers P4-0121 and P1-0201.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

Nataša Poklar Ulrih would like to thank Kenneth J. Breslauer for giving me the opportunity to learn from the best.

Conflicts of Interest

The authors declare that they have no conflict of interest.

References

  1. Bruylants, G.; Wouters, J.; Michaux, C. Differential Scanning Calorimetry in Life Science: Thermodynamics, Stability, Molecular Recognition and Application in Drug Design. Curr. Med. Chem. 2005, 12, 2011–2020. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Ouldridge, T.E. The Importance of Thermodynamics for Molecular Systems, and the Importance of Molecular Systems for Thermodynamics. Nat. Comput. 2018, 17, 3–29. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Rhodes, D.; Lipps, H.J. Survey and Summary G-Quadruplexes and Their Regulatory Roles in Biology. Nucleic Acids Res. 2015, 43, 8627–8637. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Abiri, A.; Lavigne, M.; Rezaei, M.; Nikzad, S.; Zare, P.; Mergny, J.L.; Rahimi, H.R. Unlocking G-Quadruplexes as Antiviral Targets. Pharmacol. Rev. 2021, 73, 897–923. [Google Scholar] [CrossRef]
  5. Sasisekharan, V.; Zimmerman, S.; Davies, D.R. The Structure of Helical 5′-Guanosine Monophosphate. J. Mol. Biol. 1975, 92, 171–179. [Google Scholar] [CrossRef]
  6. Gellert, M.; Lipsett, M.N.; Davies, D.R. Helix Formation BY Guanylic Acid. Proc. Natl. Acad. Sci. USA 1962, 48, 2013–2018. [Google Scholar] [CrossRef] [Green Version]
  7. Largy, E.; Mergny, J.L.; Gabelica, V. Role of Alkali Metal Ions in G-Quadruplex Nucleic Acid Structure and Stability. In Metal Ions in Life Sciences; Springer: Berlin/Heidelberg, Germany, 2016; Volume 16, pp. 203–258. [Google Scholar]
  8. Kettani, A.; Gorin, A.; Majumdar, A.; Hermann, T.; Skripkin, E.; Zhao, H.; Jones, R.; Patel, D.J. A Dimeric DNA Interface Stabilized by Stacked A·(G·G·G·G)·A Hexads and Coordinated Monovalent Cations. J. Mol. Biol. 2000, 297, 627–644. [Google Scholar] [CrossRef] [Green Version]
  9. Strahan, G.D.; Keniry, M.A.; Shafer, R.H. NMR Structure Refinement and Dynamics of the K+-[d(G3T4G3)]2 Quadruplex via Particle Mesh Ewald Molecular Dynamics Simulations. Biophys. J. 1998, 75, 968–981. [Google Scholar] [CrossRef] [Green Version]
  10. Grün, J.T.; Schwalbe, H. Folding Dynamics of Polymorphic G-Quadruplex Structures. Biopolymers 2022, 113, e23477. [Google Scholar] [CrossRef]
  11. Bončina, M.; Lah, J.; Prislan, I.; Vesnaver, G. Energetic Basis of Human Telomeric DNA Folding into G-Quadruplex Structures. J. Am. Chem. Soc. 2012, 134, 9657–9663. [Google Scholar] [CrossRef]
  12. Čeru, S.; Šket, P.; Prislan, I.; Lah, J.; Plavec, J. A New Pathway of DNA G-Quadruplex Formation. Angew. Chem. 2014, 126, 4981–4984. [Google Scholar] [CrossRef]
  13. Marquevielle, J.; Robert, C.; Lagrabette, O.; Wahid, M.; Bourdoncle, A.; Xodo, L.E.; Mergny, J.L.; Salgado, G.F. Structure of Two G-Quadruplexes in Equilibrium in the KRAS Promoter. Nucleic Acids Res. 2020, 48, 9336–9345. [Google Scholar] [CrossRef] [PubMed]
  14. 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 EGFRpromoter. Nucleic Acids Res. 2017, 45, 10132–10142. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Lim, K.W.; Lacroix, L.; Yue, D.J.E.; Lim, J.K.C.; Lim, J.M.W.; Phan, A.T. Coexistence of Two Distinct G-Quadruplex Conformations in the HTERT Promoter. J. Am. Chem. Soc. 2010, 132, 12331–12342. [Google Scholar] [CrossRef]
  16. Bhattacharyya, D.; Arachchilage, G.M.; Basu, S. Metal Cations in G-Quadruplex Folding and Stability. Front. Chem. 2016, 4, 38. [Google Scholar] [CrossRef] [Green Version]
  17. Gray, R.D.; Trent, J.O.; Chaires, J.B. Folding and Unfolding Pathways of the Human Telomeric G-Quadruplex. J. Mol. Biol. 2014, 426, 1629–1650. [Google Scholar] [CrossRef] [Green Version]
  18. Marchand, A.; Gabelica, V. Folding and Misfolding Pathways of G-Quadruplex DNA. Nucleic Acids Res. 2016, 44, 10999–11012. [Google Scholar] [CrossRef]
  19. Stadlbauer, P.; Kührová, P.; Banáš, P.; Koča, J.; Bussi, G.; Trantírek, L.; Otyepka, M.; Šponer, J. Hairpins Participating in Folding of Human Telomeric Sequence Quadruplexes Studied by Standard and T-REMD Simulations. Nucleic Acids Res. 2015, 43, 9626–9644. [Google Scholar] [CrossRef] [Green Version]
  20. Miller, M.C.; Le, H.T.; Dean, W.L.; Holt, P.A.; Chaires, J.B.; Trent, J.O. Polymorphism and Resolution of Oncogene Promoter Quadruplex-Forming Sequences. Org. Biomol. Chem. 2011, 9, 7633–7637. [Google Scholar] [CrossRef] [Green Version]
  21. Grün, J.T.; Blümler, A.; Burkhart, I.; Wirmer-Bartoschek, J.; Heckel, A.; Schwalbe, H. Unraveling the Kinetics of Spare-Tire DNA G-Quadruplex Folding. J. Am. Chem. Soc. 2021, 143, 6185–6193. [Google Scholar] [CrossRef]
  22. Zhang, A.Y.Q.; Balasubramanian, S. The Kinetics and Folding Pathways of Intramolecular G-Quadruplex Nucleic Acids. J. Am. Chem. Soc. 2012, 134, 19297–19308. [Google Scholar] [CrossRef] [PubMed]
  23. 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] [PubMed]
  24. Prislan, I.; Lah, J.; Milanic, M.; Vesnaver, G. Kinetically Governed Polymorphism of d(G4T4G 3) Quadruplexes in K+ Solutions. Nucleic Acids Res. 2011, 39, 1933–1942. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Prislan, I.; Lah, J.; Vesnaver, G. Diverse Polymorphism of G-Quadruplexes as a Kinetic Phenomenon. J. Am. Chem. Soc. 2008, 130, 14161–14169. [Google Scholar] [CrossRef]
  26. Dalgleish, D. Biophysical Chemistry: Part III ’The Behaviour of Biological Macromolecules. Biochem. Educ. 1981, 9, 157. [Google Scholar] [CrossRef]
  27. Seiler, M.C.; Seiler, F.A. Numerical Recipes in C: The Art of Scientific Computing. Risk Anal. 1989, 9, 415–416. [Google Scholar] [CrossRef]
  28. Rujan, I.N.; Meleney, J.C.; Bolton, P.H. Vertebrate Telomere Repeat DNAs Favor External Loop Propeller Quadruplex Structures in the Presence of High Concentrations of Potassium. Nucleic Acids Res. 2005, 33, 2022–2031. [Google Scholar] [CrossRef] [Green Version]
  29. Kypr, J.; Kejnovská, I.; Renčiuk, D.; Vorlíčková, M. Circular Dichroism and Conformational Polymorphism of DNA. Nucleic Acids Res. 2009, 37, 1713–1725. [Google Scholar] [CrossRef] [Green Version]
  30. Chang, C.C.; Chien, C.W.; Lin, Y.H.; Kang, C.C.; Chang, T.C. Investigation of Spectral Conversion of d(TTAGGG)4 and d(TTAGGG)13 upon Potassium Titration by a G-Quadruplex Recognizer BMVC Molecule. Nucleic Acids Res. 2007, 35, 2846–2860. [Google Scholar] [CrossRef] [Green Version]
  31. Luu, K.N.; Phan, A.T.; Kuryavyi, V.; Lacroix, L.; Patel, D.J. Structure of the Human Telomere in K+ Solution: An Intramolecular (3 + 1) G-Quadruplex Scaffold. J. Am. Chem. Soc. 2006, 128, 9963–9970. [Google Scholar] [CrossRef] [Green Version]
  32. Sun, X.; Cao, E.; He, Y.; Qin, J. Fluorescence Studies on the Interaction of Ethidium Bromide with Duplex, Triplex and Quadruplex DNA Structures. Sci. China Ser. B Chem. 1999, 42, 62–69. [Google Scholar] [CrossRef]
  33. Petraccone, L.; Pagano, B.; Esposito, V.; Randazzo, A.; Piccialli, G.; Barone, G.; Mattia, C.A.; Giancola, C. Thermodynamics and Kinetics of PNA-DNA Quadruplex-Forming Chimeras. J. Am. Chem. Soc. 2005, 127, 16215–16223. [Google Scholar] [CrossRef] [PubMed]
  34. Ren, J.; Qu, X.; Trent, J.O.; Chaires, J.B. Tiny Telomere DNA. Nucleic Acids Res. 2002, 30, 2307–2315. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Jin, R.; Gaffney, B.L.; Wang, C.; Jones, R.A.; Breslauer, K.J. Thermodynamics and Structure of a DNA Tetraplex: A Spectroscopic and Calorimetric Study of the Tetramolecular Complexes of d(TG3T) and d(TG3T2G3T). Proc. Natl. Acad. Sci. USA 1992, 89, 8832–8836. [Google Scholar] [CrossRef] [Green Version]
  36. Carrino, S.; Hennecker, C.D.; Murrieta, A.C.; Mittermaier, A. Frustrated Folding of Guanine Quadruplexes in Telomeric DNA. Nucleic Acids Res. 2021, 49, 3063–3076. [Google Scholar] [CrossRef]
  37. Long, X.; Stone, M.D. Kinetic Partitioning Modulates Human Telomere DNA G-Quadruplex Structural Polymorphism. PLoS ONE 2013, 8, e83420. [Google Scholar] [CrossRef] [Green Version]
  38. Müller, D.; Bessi, I.; Richter, C.; Schwalbe, H. The Folding Landscapes of Human Telomeric RNA and DNA G-Quadruplexes Are Markedly Different. Angew. Chem. Int. Ed. 2021, 60, 10895–10901. [Google Scholar] [CrossRef]
  39. Bardin, C.; Leroy, J.L. The Formation Pathway of Tetramolecular G-Quadruplexes. Nucleic Acids Res. 2008, 36, 477–488. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Proposed kinetic model for the thermally induced folding/unfolding transitions of d(G4T4G3) in K+ solutions which assumes the coexistence and interconversions of one tetramolecular quadruplex T4, three bimolecular quadruplexes A2, B2 and C2 and the unfolded single strands S. kij are rate constants for each i→j transition step.
Figure 1. Proposed kinetic model for the thermally induced folding/unfolding transitions of d(G4T4G3) in K+ solutions which assumes the coexistence and interconversions of one tetramolecular quadruplex T4, three bimolecular quadruplexes A2, B2 and C2 and the unfolded single strands S. kij are rate constants for each i→j transition step.
Life 12 00825 g001
Figure 2. Temperature dependence of CD spectra of d(G4T4G3) quadruplexes in K+ solutions (25, 35 and 50 mM): (a) d(G4T4G3) quadruplexes prepared at the cooling rate of 0.05 °C/min; (b) d(G4T4G3) quadruplexes prepared at the cooling rate of 1 °C/min. The measured ellipticity was normalized to 1 M single strand concentration and 1 cm light path length.
Figure 2. Temperature dependence of CD spectra of d(G4T4G3) quadruplexes in K+ solutions (25, 35 and 50 mM): (a) d(G4T4G3) quadruplexes prepared at the cooling rate of 0.05 °C/min; (b) d(G4T4G3) quadruplexes prepared at the cooling rate of 1 °C/min. The measured ellipticity was normalized to 1 M single strand concentration and 1 cm light path length.
Life 12 00825 g002
Figure 3. Non denaturing PAGE (20%) of d(G4T4G3) quadruplexes prepared at the cooling rate of either 1 °C/min (1st lane) or 0.05 °C/min (2nd lane) and followed by UV shadowing (λ = 254 nm) and Fluorescence (λ = 590 nm): (a) d(G4T4G3) quadruplexes prepared in 25 mM K+ solution and measured at a constant temperature of 5 °C. As a control oligonucleotide double stranded d(5′-AGAAGAAAAGA-3′; 5′-TCTTTTCTTCT-3′) (3rd lane) and 11-mere single stranded d(5′-TCTTTTCTTCT-3′) (4th lane) markers were used;): (b) d(G4T4G3) quadruplexes prepared in 50 mM K+ solution and measured at a constant temperature of 5 °C. As control oligonucleotide, double stranded d(5′-AGAAGAAAAGA-3′; 5′-TCTTTTCTTCT-3′) (3rd lane), single stranded d(5′-TCTTTTCTTCT-3′) (4th lane) and single stranded d(5′-AGAAGAAAAGA-3′) (5th lane) markers were used.
Figure 3. Non denaturing PAGE (20%) of d(G4T4G3) quadruplexes prepared at the cooling rate of either 1 °C/min (1st lane) or 0.05 °C/min (2nd lane) and followed by UV shadowing (λ = 254 nm) and Fluorescence (λ = 590 nm): (a) d(G4T4G3) quadruplexes prepared in 25 mM K+ solution and measured at a constant temperature of 5 °C. As a control oligonucleotide double stranded d(5′-AGAAGAAAAGA-3′; 5′-TCTTTTCTTCT-3′) (3rd lane) and 11-mere single stranded d(5′-TCTTTTCTTCT-3′) (4th lane) markers were used;): (b) d(G4T4G3) quadruplexes prepared in 50 mM K+ solution and measured at a constant temperature of 5 °C. As control oligonucleotide, double stranded d(5′-AGAAGAAAAGA-3′; 5′-TCTTTTCTTCT-3′) (3rd lane), single stranded d(5′-TCTTTTCTTCT-3′) (4th lane) and single stranded d(5′-AGAAGAAAAGA-3′) (5th lane) markers were used.
Life 12 00825 g003
Figure 4. Fitting the model function to the heating and cooling DSC thermograms of d(G4T4G3) quadruplexes in the presence of 25, 35 and 50 mM K+ ions: (a) The Δcp vs. T curves measured at the heating rate of 1 °C/min in the presence of 25 mM K+ ions for samples prepared by controlled cooling in an outer thermostat at the rate of 0.05 °C/min (red line) or 1.0 °C/min (blue line); (b) The Δcp vs. T curves measured at the heating rate of 1 °C/min in the presence of 35 mM K+ ions for samples prepared by controlled cooling in an outer thermostat at the rate of 0.05 °C/min (red line) or 1.0 °C/min (blue line); (c) The Δcp vs. T curves measured at the heating rate of 1 °C/min in the presence of 50 mM K+ ions for samples prepared by controlled cooling in an outer thermostat at the rate of 0.05 °C/min (red line) or 1.0 °C/min (blue line); (d) The Δcp vs. T curves measured at the cooling rate of 1 °C/min in the presence of 25 (blue line), 35 (red line) and 50 (green line) mM K+ ions. In all panels full lines represent model–based Δcp vs. T curves calculated from (1).
Figure 4. Fitting the model function to the heating and cooling DSC thermograms of d(G4T4G3) quadruplexes in the presence of 25, 35 and 50 mM K+ ions: (a) The Δcp vs. T curves measured at the heating rate of 1 °C/min in the presence of 25 mM K+ ions for samples prepared by controlled cooling in an outer thermostat at the rate of 0.05 °C/min (red line) or 1.0 °C/min (blue line); (b) The Δcp vs. T curves measured at the heating rate of 1 °C/min in the presence of 35 mM K+ ions for samples prepared by controlled cooling in an outer thermostat at the rate of 0.05 °C/min (red line) or 1.0 °C/min (blue line); (c) The Δcp vs. T curves measured at the heating rate of 1 °C/min in the presence of 50 mM K+ ions for samples prepared by controlled cooling in an outer thermostat at the rate of 0.05 °C/min (red line) or 1.0 °C/min (blue line); (d) The Δcp vs. T curves measured at the cooling rate of 1 °C/min in the presence of 25 (blue line), 35 (red line) and 50 (green line) mM K+ ions. In all panels full lines represent model–based Δcp vs. T curves calculated from (1).
Life 12 00825 g004
Figure 5. Fitting the model function to heating DSC thermograms of d(G4T4G3) quadruplexes at different rates in the presence of 25 (blue color), 35 (red color) and 50 (green color) mM K+ ions: (a) at the heating rate of 0.5 °C/min for samples prepared at the cooling rate of 0.05 °C/min; (b) at the heating rate of 0.5 °C/min for samples prepared at the cooling rate of 1.0 °C/min; (c) at the heating rate of 2.0 °C/min for samples prepared at the cooling rate of 0.05 °C/min; (d) at the heating rate of 2.0 °C/min for samples prepared at the cooling rate of 1.0 °C/min. In all panels full lines represent the model—based curves calculated from (1).
Figure 5. Fitting the model function to heating DSC thermograms of d(G4T4G3) quadruplexes at different rates in the presence of 25 (blue color), 35 (red color) and 50 (green color) mM K+ ions: (a) at the heating rate of 0.5 °C/min for samples prepared at the cooling rate of 0.05 °C/min; (b) at the heating rate of 0.5 °C/min for samples prepared at the cooling rate of 1.0 °C/min; (c) at the heating rate of 2.0 °C/min for samples prepared at the cooling rate of 0.05 °C/min; (d) at the heating rate of 2.0 °C/min for samples prepared at the cooling rate of 1.0 °C/min. In all panels full lines represent the model—based curves calculated from (1).
Life 12 00825 g005
Figure 6. Ratio between calculated fractions of bimolecular quadruplex species C2 and A2 at different K+ concentrations prepared at cooling rates of 0.05 and 1.0 °C/min on a logarithm scale.
Figure 6. Ratio between calculated fractions of bimolecular quadruplex species C2 and A2 at different K+ concentrations prepared at cooling rates of 0.05 and 1.0 °C/min on a logarithm scale.
Life 12 00825 g006
Figure 7. Normalized CD melting curves at λ = 290 nm of d(G4T4G3) quadruplexes obtained from experiments (dots) and calculated from fractions of quadruplex species present in the solution ( α A , α B , α C and α T ) (full line) after slow (blue color) and fast (red color) cooling: (a) in the presence of 25 mM K+ ions; (b) in the presence of 35 mM K+ ions; (c) in the presence of 50 mM K+ ions.
Figure 7. Normalized CD melting curves at λ = 290 nm of d(G4T4G3) quadruplexes obtained from experiments (dots) and calculated from fractions of quadruplex species present in the solution ( α A , α B , α C and α T ) (full line) after slow (blue color) and fast (red color) cooling: (a) in the presence of 25 mM K+ ions; (b) in the presence of 35 mM K+ ions; (c) in the presence of 50 mM K+ ions.
Life 12 00825 g007
Table 1. Calculated adjustable parameters for the kinetic model in Figure 1 which were used to describe thermally induced folding/unfolding transitions. The unit of e A i j is s−1 for unfolding and s−1 M−1 for folding transitions of A2, B2 and C2 and unit of E i j is cal/mol. Experimentally determined enthalpies of transition, ΔH (cal/mol single strands) were as follows: ΔHA = ΔHB = ΔHC = (3.1 ± 0.2)⋅104, ΔHT = (3.6 ± 0.2)⋅104.
Table 1. Calculated adjustable parameters for the kinetic model in Figure 1 which were used to describe thermally induced folding/unfolding transitions. The unit of e A i j is s−1 for unfolding and s−1 M−1 for folding transitions of A2, B2 and C2 and unit of E i j is cal/mol. Experimentally determined enthalpies of transition, ΔH (cal/mol single strands) were as follows: ΔHA = ΔHB = ΔHC = (3.1 ± 0.2)⋅104, ΔHT = (3.6 ± 0.2)⋅104.
ParameterCK+ = 25 mMCK+ = 35 mMCK+ = 50 mM
ABS74 ± 579 ± 581 ± 5
ASB19 ± 319 ± 319 ± 3
AAS49 ± 249 ± 249 ± 2
ATS100 ± 10100 ± 10100 ± 10
ACS54 ± 253 ± 250 ± 2
ABA15.5 ± 0.815.0 ± 0.814.9 ± 0.8
AAB7.2 ± 0.47.2 ± 0.47.2 ± 0.4
ABT66 ± 466 ± 466 ± 4
ABC24 ± 225 ± 230 ± 2
EBS(3.9 ± 0.2) × 104(4.1 ± 0.2) × 104(4.0 ± 0.2) × 104
EAS(3.2 ± 0.1) × 104(3.2 ± 0.1) × 104(3.2 ± 0.1) × 104
ETS(6.0 ± 0.2) × 104(6.0 ± 0.2) × 104(6.0 ± 0.2) × 104
ECS(3.7 ± 0.1) × 104(3.6 ± 0.1) × 104(3.7 ± 0.1) × 104
EBA(7.0 ± 0.5) × 103(7.0 ± 0.5) × 103(7.0 ± 0.5) × 103
EBT(3.2 ± 0.1) × 104(3.2 ± 0.1) × 104(3.2 ± 0.1) × 104
EBC(1.5 ± 0.1) × 104(1.4 ± 0.1) × 104(1.7 ± 0.1) × 104
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Prislan, I.; Urbic, T.; Poklar Ulrih, N. Thermally Induced Transitions of d(G4T4G3) Quadruplexes Can Be Described as Kinetically Driven Processes. Life 2022, 12, 825. https://doi.org/10.3390/life12060825

AMA Style

Prislan I, Urbic T, Poklar Ulrih N. Thermally Induced Transitions of d(G4T4G3) Quadruplexes Can Be Described as Kinetically Driven Processes. Life. 2022; 12(6):825. https://doi.org/10.3390/life12060825

Chicago/Turabian Style

Prislan, Iztok, Tomaz Urbic, and Natasa Poklar Ulrih. 2022. "Thermally Induced Transitions of d(G4T4G3) Quadruplexes Can Be Described as Kinetically Driven Processes" Life 12, no. 6: 825. https://doi.org/10.3390/life12060825

APA Style

Prislan, I., Urbic, T., & Poklar Ulrih, N. (2022). Thermally Induced Transitions of d(G4T4G3) Quadruplexes Can Be Described as Kinetically Driven Processes. Life, 12(6), 825. https://doi.org/10.3390/life12060825

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