Next Article in Journal
Multi-Scale X-ray Imaging of the Pigment Discoloration Processes Triggered by Chlorine Compounds in the Upper Basilica of Saint Francis of Assisi
Previous Article in Journal
Ecotoxicity of Diazinon and Atrazine Mixtures after Ozonation Catalyzed by Na+ and Fe2+ Exchanged Montmorillonites on Lemna minor
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Theoretical Assessment of the Ligand/Metal/Quadruplex Recognition in the Non-Canonical Nucleic Acids Structures

1
Institute of Optical Materials and Technologies “Acad. J. Malinowski”, Bulgarian Academy of Sciences, 1113 Sofia, Bulgaria
2
Institute of Organic Chemistry with Centre of Phytochemistry, Bulgarian Academy of Sciences, 1113 Sofia, Bulgaria
*
Author to whom correspondence should be addressed.
Molecules 2023, 28(16), 6109; https://doi.org/10.3390/molecules28166109
Submission received: 5 July 2023 / Revised: 11 August 2023 / Accepted: 14 August 2023 / Published: 17 August 2023
(This article belongs to the Section Physical Chemistry)

Abstract

:
Quadruplexes (GQs), peculiar DNA/RNA motifs concentrated in specific genomic regions, play a vital role in biological processes including telomere stability and, hence, represent promising targets for anticancer therapy. GQs are formed by folding guanine-rich sequences into square planar G-tetrads which stack onto one another. Metal cations, most often potassium, further stabilize the architecture by coordinating the lone electron pairs of the O atoms. The presence of additional nucleic acid bases, however, has been recently observed experimentally and contributes substantially to the structural heterogeneity of quadruplexes. Therefore, it is of paramount significance to understand the factors governing the underlying complex processes in these structures. The current study employs DFT calculations to model the interactions between metal cations (K+, Na+, Sr2+) and diverse tetrads composed of a guanine layer in combination with a guanine (G)-, adenine (A)-, cytosine (C)-, thymine (T)-, or uracil (U)-based tetrad layer. Moreover, the addition of 4-(3,4-dihydroisoquinolin-2-yl)-2-(quinolin-2-yl)quinazoline to the modeled quadruplexes as a possible mechanism of its well-exerted antitumor effect is assessed. The calculations imply that the metal cation competition and ligand complexation are influenced by the balance between electronic and implicit/explicit solvation effects, the composition of the tetrad layers, as well as by the solvent exposure to the surrounding environment expressed in terms of different dielectric constant values. The provided results significantly enhance our understanding of quadruplex diversity, ligand recognition, and the underlying mechanisms of stabilization at an atomic level.

1. Introduction

Ever since the fundamental suggestion of the DNA double-stranded right-handed helix structure by Watson and Crick [1], the quest for unraveling its enigmatic and fascinating composition has tempted many researchers from diverse scientific fields. Since then, however, numerous peculiar DNA/RNA motifs have been discovered, with one such example posing the so-called quadruplexes. Notably, these motifs are enriched in certain regions rather than being randomly distributed throughout the genome. This strictly corresponds to the key part they play in the regulation of various biological processes, e.g., mitochondrial transcription, telomere biogenesis and homeostasis, mRNA translocation, transcription, replication, and chromatin remodeling, just to name a few [2,3,4,5,6,7]. In the context of their diverse functions, especially in telomere regulation and homeostasis, quadruplexes have been recognized as a plausible drug target in anticancer therapy against most tumor formations [8,9,10,11]. Therefore, it is of utmost interest and importance to understand the underlying factors that govern the intra- and intermolecular interactions observed in their complex architecture.
Quadruplexes are most commonly formed by folding guanine-rich nucleic acid sequences in both RNA and DNA in the arrangement of square planar compositions due to the presence of complementing hydrogen bonds (N1-O6)(N2-N7) on their Hoogsteen (HG) and Watson–Crick (WC) faces (Figure 1).
The arising G-quartets (GQs) or G-tetrads stack onto one another. It is further stabilized by the addition of a metal cation, usually potassium (K+) or, much less frequently, sodium (Na+), although other (mostly alien) metal species have additionally been observed to have an effect [12,13,14,15,16]. The metal cofactor is a highly required component, which stabilizes the intricate architecture through coordination with the oxygen lone pairs in the central channel of the quadruplex. Depending on the chemical characteristics of the metal cation (ionic radius, charge, coordination preferences) and the structure of the GQs, the Mn+ is positioned either within the approximate tetrad plane or between two planes. In the nucleic sequence, the intervening bases between the G-tetrads build the so-called loops whose function is to link the stacked G-quartets [17,18,19]. Hence, the typical composition of a quadruplex consists of three main constituents: (1) layers of rigid coplanar G-tetrads, (2) flexible linking loops, and (3) a central alkali metal cation, which is present in almost all structures. This structural polymorphism results in great diversity in loop and tetrad arrangements, capping structures, and strand orientation. Furthermore, their various conformations (parallel, antiparallel, hybrid-I and II) are not constant but could undergo conversion when the environmental conditions change [20]. With the development of experimental techniques like nuclear magnetic resonance (NMR), circular dichroism (CD), and gel electrophoresis, some interesting nontypical tetrads have been recently observed, namely A-/C-/U- and T-tetrads [21,22,23,24,25,26,27,28,29,30]. These are commonly stacked over the adjacent G-quartets and add significantly to the structural diversity and recognition especially in telomeres.
The unceasing search for novel anticancer therapeutics that exert not only selective effect but possess a different mechanism of action that could circumvent the increasing tumor resistance against conventional drugs has led to the recognition of the crescent-shaped planar heterocycles targeting the quadruplex architecture [31,32,33,34,35]. Thus, the resulting stabilization and/or induction of quadruplexes due to interaction with the ligand may lead to (1) maintaining telomere problems (decreasing telomerase activity and increasing telomere stability); (2) reducing oncogene expression by inhibiting transcription and translation; (3) increasing genome instability; hence, inducing apoptosis (cell death). In this regard, immense effort has been made, but still, no applicable drug has reached the market up to date. Herewith, we continue our previous work [36,37] on studying the possible complexation between heterocyclic molecules (quinazoline derivatives) and G-quadruplexes. Our main goal in the current study is, however, to decipher the principal factors that govern the recognition between three metal cations (K+, Na+, and Sr2+) and the intriguing non-canonical quadruplexes by modeling their interaction with combinations of tetrads. Furthermore, we evaluate the possible stacking of the most active (lowest IC50) 2-quinolinyl-quinazoline derivative reported previously in [37] to the modeled tetrads. The results are obtained by utilizing a well-studied and validated DFT methodology that has proven extensively reliable when the system under study contains competing metal cofactors which bind to biological or biomimetic molecules through mostly electrostatic interactions [38,39,40,41,42]. Thus, we aim to contribute to the understanding of the underlying processes in quadruplex stabilization and its further complexation to drug molecules.

2. Results

2.1. Modeled Reactions

The main objective of the current study is to outline the specific effects that control the metal recognition in nucleic tetrads, as well as the influence that these conditions may have upon the addition of an outer ligand with potential anticancer activity. In order to achieve this goal, we modeled the following reactions of competition and addition:
Na(H2O)6Ⴈ+ + K@GQႨ+ → Na@GQႨ+ + K(H2O)6Ⴈ+
Sr(H2O)8Ⴈ2+ + K@GQႨ+ → Sr@GQႨ2+ + K(H2O)6Ⴈ+ + 2H2O
L + K@GQ → L@K@GQႨ+
where Na(H2O)6Ⴈ+, K(H2O)6Ⴈ+, and Sr(H2O)8Ⴈ2+ represent the corresponding hydrated species of the studied metal cations; M@GQႨn+ stands for the metal-containing structure of two stacked tetrads; the considered ligand L is 4-(3,4-dihydroisoquinolin-2-yl)-2-(quinolin-2-yl)quinazoline in its “concave” form; and the L@K@GQႨ+ complex depicts the resulting adduct between the ligand and the initial metal-containing quadruplex structure yielded in Reaction (R3). This particular ligand was chosen as it was previously demonstrated to express the most potent cytotoxicity effect on A375 melanoma cells (81% and 91% reduction at the low and higher concentrations, respectively) according to Ref. [37]. Four different rotamers of the ligand were furthermore studied, and the most stable one was considered as the structure of choice in the current calculations (Figure S1).
The number of water molecules surrounding the metal species under study varies greatly (between 5 and 7 for Na+, from 6 to 8 for K+, and from 5 to 9 for Sr2+) depending on the environment and experimental/theoretical approach used [43,44,45,46]. For the purposes of the current study and based on our previous work with these metal species, the Na(H2O)6Ⴈ+, K(H2O)6Ⴈ+, and Sr(H2O)8Ⴈ2+ constructs were taken into account when modeling the envisioned reactions [41,47]. These substantially correspond to experimental and theoretical data provided and reproduce structural characteristics such as bond lengths and valence angles: the 4 + 2 composition of the sodium and potassium-hydrated species as seen in [43,48,49] along with the 6 + 2 preference of the strontium cation [50]. These particular compositions are consistent with structures in which the first coordination sphere consists of four/six water molecules for the mono-and divalent cations, respectively, whereas the remaining two water molecules migrate to the second one, distanced at more than 4 Å away from the metal cation but bound to the structure by hydrogen bonds. These were found to be the most energetically favorable geometries for the studied ions among others taken into consideration, as presented in Figure S1. Furthermore, the quadruplex architectures are composed by two layers, one of which is a G-quartet, whereas the other one is built by four guanines, adenines, cytosines, thymines, or uracils, henceforth referred to as GG, AG, CG, TG, and UG in the text. Their initial geometries are taken from the data deposited in the Protein Data Bank structures—PDB Entries 3IBK [51], 1EVN [22], 1EVO [24], 1EMQ [23], and 6GE1 [52], correspondingly, followed by full geometry optimization and vibrational frequency analysis (see Materials and Methods). Note that all chosen initial structures accept parallel geometry. The architectures of the stacked tetrads are considered as two layers of four building nucleic bases each, held by the arising Hoogsteen interactions and the complementing metal cation complexation due to the computational restriction of the DFT methodology concerning the size of the studied systems. Thus, the linking loops as well as the phosphate moieties at the side chain of the quadruplexes are omitted. Note, however, that this simplification does not significantly alter the overall architecture, as seen in Figure 2, where the PDB deposited and the optimized K@GGႨ+ are juxtaposed.

2.2. Replacement of K+ with Na+

The first step in the current study is modeling the substitution of K+ with Na+ in the GQ structures. The optimized geometries of the M@GQႨ+ complexes, along with the obtained results at the wB97XD/6-31+G(d,p)//wB97XD/6-31G(d,p) level of theory in accordance to Reaction (R1), are presented in Figure 3.
The calculations reveal that the most common and widely distributed G-quadruplex built by G-stacked tetrads is protected against alien attack: the ∆Gε values stay on positive ground ranging from 5.3 to 6.4 kcal mol−1. This result corresponds well to experimental data claiming the inability of the sodium cation to compete with the native potassium one in these structures [13,14,16,19,35]. Moreover, the obtained positive changes in the ∆Gε values for the substitution of K+ with Na+ in the CG, TG, and UG compositions add further evidence to this concept. Interestingly, the AG structure appears vulnerable to a sodium attack, as the ∆G2/4/30 values are −1.1/−0.4/0.6. Calculations at the lower wB97XD/6-31G(d,p)//wB97XD/6-31G(d,p) level show a similar trend, although with lower ∆G numbers in absolute value (compare results in Table S1). Note that Na+ takes the position in the plane of the G-tetrad in the Na@AGႨ+ and Na@CGႨ+ complexes (at about 0.2 Å closer to the guanine quadruplex), whereas it stands between two layers in all other Na@GQႨ+. The substitution of the native potassium with the alien sodium ion in these constructs leads either to energy gain (in the case of an AG tetrad) or yields the lowest positive ∆Gε values. Hence, the former octahedral configuration appears more stable for the smaller sodium cation as opposed to its latter bipyramidal geometry, as it provides greater interaction with the bases of the layer. Both of these compositions are common, in addition to the unsymmetrical geometry of Na+ [12,19]; however, the Na+/K+ rivalry is most probably dominated by the intracellular concentration of the cations (5–15 mM Na+ and 140 mM K+ [53]) followed by the position of the sodium ion.

2.3. Replacement of K+ with Sr2+

The next step in the current study is modeling the substitution of K+ with Sr2+ in the GQ structures. The optimized geometries of the M@GQႨ+/2+ complexes, along with the obtained results at the wB97XD/6-31+G(d,p)//wB97XD/6-31G(d,p) level of theory in accordance to Reaction (R2), are presented in Figure 4.
The performed calculations at the two different levels (Figure 4 and Table S1) provide interesting results regarding the importance of diverse factors controlling the native/abiogenic ion competition. The alien strontium cation effectively replaces the native K+ in all studied constructs in a nonpolar environment characterized by a dielectric constant ε = 2, as all of the obtained ∆G2 values are strongly negative: ranging from −0.9 kcal mol−1 for the substitution in the M@AGႨ+/2+ complex to −15.5 kcal mol−1 in the M@CGႨ+/2+ structure. This outcome should be attributed to the greater complexation ability of Sr2+ due to its higher charge in combination with the plausible ionic radius resembling that of the biogenic potassium (1.32 Å vs. 1.38 Å for Sr2+/K+, respectively [54]). Nevertheless, with increasing the dielectric constant of the medium, corresponding to greater solvent exposure, the substitution of K+ with Sr2+ becomes thermodynamically unfavorable in all modeled reactions (positive ∆G4/30 values), with the M@CGႨ+/2+ structure posing the only exception where the reaction yields a positive ∆G4. The explanation lies in the enhanced desolvation penalty for the strontium cation in the more polar environment. Noteworthy, the polarity of the medium plays a significant role, as with increasing the dielectric constant the energy gain decreases substantially. Interestingly, these results correlate to contradicting experimental data, as some authors provide clear evidence of the ability of strontium to substitute the native potassium [12], whereas others claim the opposite trend [13]. Our computational approach delineates the importance of the medium surrounding these DNA/RNA motifs, as well as the structure of the building tetrad that undoubtedly affects the outcome of the Sr2+/K+ rivalry.

2.4. Ligand Complexation by the Studied K@GQႨ+ Structures

As a continuation of our previous work [36,37], we further modeled the reaction of complexation between the most potent heterocyclic molecule exhibiting the greatest cytotoxic effect and all five types of K@GQႨ+ complexes in order to assess its ability to bind them. The most stable rotamer of ligand L was used in modeling the complexes with K@GQႨ+ (the relative stabilities of possible rotamers of L are shown in (Figure S1, Supplementary Materials). The optimized structures of the ternary complexes denoted as L@K@GQႨ+ are depicted in Figure 5. The obtained ∆Gε values regarding Reaction (R3) are provided as well.
The calculated Gibbs energies in the nonpolar environment stay on negative ground or close to zero, indicating a thermodynamically favorable or probable reaction of formation. The ∆Gε values vary between −2.8 and 1.6 kcal mol−1 at ε = 2, between 0.2 and 4.0 kcal mol−1 at ε = 4, and between 2.6 and 5.8 kcal mol−1 at ε = 30. This fluctuation further reveals the significant role of the surrounding medium (the reactions of formation of the L@K@GGႨ+ and L@K@AGႨ+ complexes appear thermodynamically favorable (exergonic) in the deeply buried nonpolar environment (∆G2 = −0.1/−2.8 kcal mol−1, respectively) whereas the ∆G30 values become positive (3.4/2.6 kcal mol−1, correspondingly) but still stay close to zero. Moreover, the obtained result through the currently applied computational protocol of ∆G2 = −2.4/−0.1 kcal mol−1 (higher level vs. lower level of theory) for the formation of L@K@GGႨ+ provides evidence of the expressed hypothesis that the ligand’s mechanism of action is due to stacking to the quadruplex architecture, hence the previously reported low IC50 value of 1.43 × 10−3 M. Notably, the rest of the GQ-stacked tetrads further appear prone to positive interaction with the ligand under study, evidenced by the negative ∆G2 values at the wB97XD/6-31G(d,p)//wB97XD/6-31G(d,p), and, all the more, the K@AGႨ+ motif with negative values even at the higher wB97XD/6-31+G(d,p)//wB97XD/6-31G(d,p) level of theory. This outcome unambiguously proves that 4-(3,4-dihydroisoquinolin-2-yl)-2-(quinolin-2-yl)quinazoline is indeed a potent anticancer drug exerting its effect (at least partially) by successfully binding the quadruplex architectures in telomeres.

3. Discussion

3.1. Effect of the Tetrad Nature

G-stacked tetrads constitute most commonly quadruplexes in telomeres playing a role of high significance in cell apoptosis. Recent studies, however, have revealed even more peculiar motifs, where in addition to the G-layers, tetrads composed of other nucleic bases, both of a purine and pyrimidine nature, are observed. The current study delineates the importance of the present bases for the outcome in metal competition as well as ligand addition to the structures. The absence of carbonyl groups in the adenine tetrad, as well as the presence of only one per base in the cytidine quadruplex, predispose the smaller sodium cation to migration, thus occupying the position in the plane of the G-quartet in an octahedral geometry. This results in a positive outcome (negative ∆Gε values) of the competition with K+ in the M@AGႨ+ construct, which is the most commonly present metal cation stabilizing the quadruplex architecture. Notably, the presence of more carbonyl groups in the remaining studied bases inclines the metal cations to accommodate in the position between two layers, which negatively affects the outcome in the Na+ → K+ substitution (positive ∆Gε values). This is especially underlined in the case of the Na@UGႨ+ complex, where the attacking ion is unable to substitute the native potassium, evidenced by the highest yielded ∆Gε values. The same conclusion can be drawn regarding the Sr2+/K+ competition for binding the AG -composition, where ∆G2 equals −0.9 kcal mol−1. Although the result stays below zero, this is the smallest in absolute value yielded for the substitution of K+ with Sr2+ at ε = 2. Therefore, the optimization of the structures under study points out that the tetrad nature affects mostly the position of the smaller Na+, while the bulkier K+ and Sr2+ experience only loss of interaction but remain between two layers due to their ionic radii. Furthermore, the arising interaction with the NH2 groups favors the addition of 4-(3,4-dihydroisoquinolin-2-yl)-2-(quinolin-2-yl)quinazoline to the K@AGႨ+ complex to the greatest extent, indicated by the lowest ∆Gε values. Still, it should be noted that the studied ligand successfully binds to the most abundant quadruplex construct, the GG tetrad, displaying a potential for application in antitumor therapy aimed at interaction with quadruplexes in telomeres. Overall, the K@AGႨ+ compositions appear the most vulnerable towards an outer attack, either by a competing metal cation or by an external ligand binding to the complex, evidenced by the performed calculations.

3.2. Effect of the Metal Cation Properties

The conducted study reveals the significance of the metal cation properties for elucidating the outcome of the Na+/K+ and Sr2+/K+ rivalry for binding diverse G-quadruplex architectures by utilizing a well-defined and previously applied DFT methodology. This particular approach outlines the most important factors: the ionic radius, charge, and de/hydration energy of the metal cation. Notably, the obtained results fall strictly in line with previously observed experimental trends but draw a clearer picture at an atomic level. The data demonstrate that metal cations characterized by an ionic radius of around 1.35 Å bind most effectively the lone electron pairs of the O/N atoms of the tetrad bases when positioned between two layers in a bipyramidal geometry. Furthermore, the charge of the cation should also be taken into account as it corresponds to an increased ion–dipole interaction demonstrated by a Hirshfeld analysis [55]: 0.89 e/0.77 e vs. 1.98 e for the monocationic K+/Na+, respectively, vs. the dicationic Sr2+ in the M@GG complexes. The analysis draws an analogous picture in all other compositions, as the charge transferred from the ligands to the metal cation varies between 0.81 and 0.86 e for K+, 0.71 and 0.75 e for Na+, and 1.93 and 1.98 e for Sr2+. The data clearly show the increased interaction between the potassium cation in comparison with the sodium one as well as the greater ion–dipole interaction with strontium, attributed to its higher charge. Nonetheless, the desolvation penalty of the Mn+ attenuates the free energy gain, especially in the polar medium, thus settling the outcome of the metal cation rivalry. Consequently, the native potassium cation appears suitable for stabilizing the ion channel in the quadruplex architecture, but its extremely high cellular concentration should further be acknowledged. An even clearer picture regarding the significance of the nature of the metal cation in quadruplexes could be drawn in a future study focusing on the competition between K+ and a whole series of metal cations including trivalent ones, when the ionic polarity index presented in Ref. [56] could be applied.

3.3. Effect of the Dielectric Constant of the Medium

The results regarding ligand/metal/quadruplex recognition presented in the current study underline the considerable importance of the surrounding environment when modeling interactions in biological systems (inter/intracellular space, cytoplasm, nucleus, proteins). Firstly, the dielectric constant of the medium strongly influences the outcome of the metal competition at both levels of theory (Table S1). The interplay between the energy loss due to the desolvation of the incoming metal species and the energy gain from the solvation of the Mn+ products results in effective competition of the alien strontium cation against the native potassium in an environment resembling the gas phase (negative ∆G2 in all cases and negative ∆G4 values in most of the studied tetrads for the substitution of K+ with Sr2+). Yet, this trend is reversed in the more polar surroundings, which substantially amplify the ion–dipole interaction in the Sr(H2O)8Ⴈ2+ construct. An analogous conclusion can be drawn regarding the Na+/K+ rivalry where the de/solvation energy of the much smaller sodium strongly surpasses that of the bulkier potassium. Moreover, the addition of the ligand to the K@GQႨ+ complexes is affected by the change in the dielectric constant of the medium, especially in the cases of the GG and AG architectures, where the change of the polarity of the environment reverses the sign of the Gibbs energy, thus indicating a thermodynamically unfavorable reaction under these conditions. Still, the calculations imply a plausible recognition between some K@GQႨ+ structures and the incoming ligand, especially in the less polar medium characterized by a dielectric constant of two.

3.4. Effect of the Ligand Addition

The previously reported cytotoxic analysis [37] provided evidence for the anticancer activity of the 4-(3,4-dihydroisoquinolin-2-yl)-2-(quinolin-2-yl)quinazoline in addition to theoretical assessment of its affinity to the quadruplex guanine tetrads. The computational protocol utilized herewith further assessed the change in the Gibbs energy in different media, thus providing negative ∆G2 value for the formation of L@K@GGႨ+ that corresponds well to the experimentally observed low IC50 values, suggesting a probable mechanism of antitumor activity. Moreover, we demonstrate the ability of the ligand to bind effectively not only the GG-quartets, but also the AG-type non-canonical motif observed in telomeres, as the calculations stay firmly on negative ground: −2.8 kcal mol−1 in the nonpolar environment. The computed results differ with the change of the dielectric constant of the medium; hence, the values fluctuate between 0.2 and 4.0 kcal mol−1 when ε equals 4, and between 2.6 and 5.8 kcal mol−1 in the solvent-exposed environment. The calculations further imply that the tetrad most prone to ligand addition is the adenine-based quadruplex with ∆G2/4/30 values = −2.8/0.2/2.6 kcal mol−1, correspondingly, due to the arising interaction between the NH2 groups from the A bases and the N atoms in the heterocyclic rings in the ligand. This outcome provides strong evidence of the susceptibility of the non-canonical motifs to complexations to the potential drug molecule under study.

4. Materials and Methods

The DFT calculations were performed using the Gaussian 09 suite of programs [57]. The hybrid long-range and dispersion-corrected wB97XD functional [58,59] basis set, in conjunction with the 6-31G(d,p) basis set and the pseudopotential SDD [60] for the heavier Sr2+ in all cases, was employed for the full geometry optimization of the modeled structures, followed by vibrational frequency analysis where no negative values were observed. This outcome indicates a local minimum of the potential energy surface necessary for obtaining the electronic energy, Eel, the thermal energy, Eth, including the non-scaled zero-point energy, and entropy, S. They were further utilized in the calculation of the Gibbs energy value, ∆G1, in the gas phase at room temperature T = 298 K and atmospheric pressure 1 atm according to Equation (1):
∆G1 = ∆Eel + ∆Eth − T∆S
where ΔEel, ΔEth, and ΔS represent the corresponding differences between the products and the reactants in consistency with Reactions R1 to R3. The change in the number of moles ∆n during the reaction is also accounted for in ∆G1 (∆nRT ≈ P∆V). Further assessment of the effect of the surrounding medium was obtained through additional SMD calculations [61]. Three different media characterized by a dielectric constant ε = 2, 4, and 30 were chosen, corresponding to a deeply buried, partially solvent-exposed, and solvent-accessible environment, respectively. “Generic” solvent option was used with Eps (2/4/30) and EpsInf (2) parameters specified. Single-point calculations for each of the optimized structures were performed and used for the evaluation of the solvation energies, ΔGsolvε, defined as the difference between the condensed-phase and the gas-phase energies of the respective constructs. These were employed for obtaining the energies for the substitution of K+ with Na+/Sr2+ and ligand addition to the already formed K@GQႨ+ structures in accordance with the equation:
∆Gε = ∆G1 + ∆∆Gsolvε
where
∆∆Gsolvε = ∆Gsolvε (products) − ∆Gsolvε (reactants)
A positive ∆Gsolvε value characterizes a thermodynamically unfavorable reaction, whereas a negative ∆Gsolvε suggests a favorable one. In terms of selectivity, this outcome corresponds to a potassium-selective construct, while the opposite distinguishes a preference toward the foreign intruder and the addition of the outer ligand. In order to more adequately assess the underlying factors, additional single-point calculations with the triple-zeta 6-31+G(d,p) basis set for all optimized structures in the gas phase, as well as in different media, were employed. The obtained results are denoted as calculated at the wB97XD/6-31+G(d,p)//wB97XD/6-31G(d,p) level of theory. Further comparison between the two utilized schemes of computational methodologies is provided in Table S1 (Supplementary Materials), where additional presentation of the terms ∆H (accounting for ∆Eel, ∆Eth, and the change of the number of molecules at the two sides of the chemical equation) and T∆S is provided. This particular combination was chosen as it yielded positive ∆G2 value for the substitution of K+ with Na+ in the GG structure, indicating the inability of sodium to substitute the native potassium, falling in line with various experimental results observed in the literature. Note, however, that the aim of the performed calculations is not to reproduce specific values of ∆Gsolvε but, rather, to provide reliable trends that outline thermodynamic descriptors affecting the metal and ligand selectivity in nucleic quadruplexes. The PyMol molecular graphics system [62] is implemented for the visualization of the obtained results.

5. Conclusions

Guanine-rich motifs in telomeres have recently attracted the attention of scientists, not only due to their peculiarity, but also by their role in cells’ apoptosis, turning them into a plausible target in antitumor therapy. By utilizing the powerful tools of DFT computations, the current study adds to our understanding of the diversity of effects that play a role in defining the complex structure in quadruplexes. Overall, tetrad nature is of crucial importance for the modeled ligand/metal/quadruplex recognition. Apparently, the main function has been bestowed upon guanine as a leading base forming a great percentage of the quadruplex architectures by no chance, since it provides solid structural integrity evidenced by the performed computations: the layered composition is not altered either by the addition of metal cations or by the complexation of a ligand. Hence, potassium is protected by an attack from the fellow metal cation sodium in the K@GGႨ+ composition, which falls in line with experimental data. On the other hand, the presence of additional nucleic bases adds to the heterogeneity of structure and function of the quadruplexes with arising novel interactions with the constituents. Furthermore, the metal cation rivalry between the native potassium and the less observed sodium in addition to the alien strontium strongly correlates with the balance between electronic and solvation effects. Hence, the replacement of K+ with Sr2+ occurs spontaneously in the nonpolar environment characterized by ε = 2, whereas this reaction becomes thermodynamically unfavorable with the enhanced polarity of the medium, where the desolvation penalty for the doubly charged strontium increases substantially. A thought-provoking trend is observed in the case of the Na+/K+ competition in the M@AGႨ+ complex: the outer sodium effectively substitutes its contender only when its interaction with the G-layer is increased due to a migration to the plane of the tetrad. Consequently, the arising electronic interaction with the lone pairs of the O-atom from the carbonyl groups compensates for the high desolvation penalty of the small cation. This outcome is attenuated, nonetheless, by enhancing the polarity of the surrounding medium. The computations reveal that the dielectric constant of the environment correlating with exposure to the aqueous environment in the cellular interior has a crucial role, not only on the outcome in metal cation rivalry, but also on the ligand addition to these structures. In combination with previously reported cytotoxic analysis, the presented results provide further strong evidence at an atomic level of the complexation ability of 4-(3,4-dihydroisoquinolin-2-yl)-2-(quinolin-2-yl)quinazoline to the K@GGႨ+ and K@AGႨ+ compositions, thus providing a solid premise for its application as an anticancer drug. Yet, the presented DFT calculations, although reliable, shed light only on the static aspect of the intricate biochemical processes under study. Hence, as a future perspective, the utilization of MD simulations as an additional approach to the problem and/or the assessment of the effect of trivalent metal cations would undoubtedly add to our understanding of ligand/metal/quadruplex recognition in non-canonical nucleic acids structures.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules28166109/s1, Figure S1: Structures and relative stabilities of the hydrated Na+, K+, and Sr2+ ion clusters: Na(H2O)6Ⴈ+ (A), K(H2O)6Ⴈ+ (B), and Sr(H2O)8Ⴈ2+ (C); Structures and relative stabilities of rotamers of Ligand L. All structures are optimized at wb97xd/6-31g(d,p) level of theory; Table S1: ∆H, T∆S, ∆G1, ∆G2, ∆G4, and ∆G30 values in kcal mol−1 calculated at the wB97XD/6-31G(d,p)//wB97XD/6-31G(d,p) and wB97XD/6-31+G(d,p)//wB97XD/6-31G(d,p) levels of theory for reaction R1: substitution of K+ by Na+; reaction R2: substitution of K+ by Sr2+; reaction R3: addition of the ligand 4-(3,4-dihydroisoquinolin-2-yl)-2-(quinolin-2-yl)quinazoline L to the K@GQႨ+ constructs; methodology adopted for implicit solvent calculations; starting geometries for K@GGႨ+, K@AGႨ+, K@CGႨ+, K@TGႨ+, and K@UGႨ+ optimization. References [39,42,63] are cited in Supplementary Materials.

Author Contributions

Conceptualization, S.B., S.A., N.K. and J.K.; investigation, N.K., S.D., V.P. and S.A.; writing—original draft preparation, N.K.; writing—review and editing, S.A., S.B. and J.K.; visualization, S.D. and V.P.; supervision, S.A.; project administration, S.B.; funding acquisition, S.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Bulgarian National Science Fund, Grant number KP-06-N59/1 of 15 November 2021.

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

The authors gratefully acknowledge the access provided to the e-infrastructure of the NCHDC, part of the Bulgarian National Roadmap for RIs, with the financial support by the Grant No D01−168/28.07.2022.

Conflicts of Interest

The authors declare no conflict of interest.

Sample Availability

Not applicable.

References

  1. Watson, J.D.; Crick, F.H.C. Molecular Structure of Nucleic Acids. Nature 1953, 171, 737–738. [Google Scholar] [CrossRef] [PubMed]
  2. Huppert, J.L. Four-Stranded Nucleic Acids: Structure, Function and Targeting of G-Quadruplexes. Chem. Soc. Rev. 2008, 37, 1375–1384. [Google Scholar] [CrossRef] [PubMed]
  3. Bugaut, A.; Balasubramanian, S. 5′-UTR RNA G-Quadruplexes: Translation Regulation and Targeting. Nucleic Acids Res. 2012, 40, 4727–4741. [Google Scholar] [CrossRef] [PubMed]
  4. Blackburn, E.H. Structure and Function of Telomeres. Nature 1991, 350, 569–573. [Google Scholar] [CrossRef]
  5. Huppert, J.L.; Balasubramanian, S. Prevalence of Quadruplexes in the Human Genome. Nucleic Acids Res. 2005, 33, 2908–2916. [Google Scholar] [CrossRef]
  6. Fernando, H.; Reszka, A.P.; Huppert, J.; Ladame, S.; Rankin, S.; Venkitaraman, A.R.; Neidle, S.; Balasubramanian, S. A Conserved Quadruplex Motif Located in a Transcription Activation Site of the Human C-Kit Oncogene. Biochemistry 2006, 45, 7854–7860. [Google Scholar] [CrossRef]
  7. Cogoi, S.; Xodo, L.E. G-Quadruplex Formation within the Promoter of the KRAS Proto-Oncogene and Its Effect on Transcription. Nucleic Acids Res. 2006, 34, 2536–2549. [Google Scholar] [CrossRef]
  8. Han, H.; Hurley, L.H. G-Quadruplex DNA: A Potential Target for Anti-Cancer Drug Design. Trends Pharmacol. Sci. 2000, 21, 136–142. [Google Scholar] [CrossRef]
  9. Neidle, S. Human Telomeric G-Quadruplex: The Current Status of Telomeric G-Quadruplexes as Therapeutic Targets in Human Cancer. FEBS J. 2010, 277, 1118–1125. [Google Scholar] [CrossRef]
  10. Che, T.; Wang, Y.Q.; Huang, Z.L.; Tan, J.H.; Huang, Z.S.; Chen, S. Bin Natural Alkaloids and Heterocycles as G-Quadruplex Ligands and Potential Anticancer Agents. Molecules 2018, 23, 493. [Google Scholar] [CrossRef]
  11. Hu, M.H.; Chen, S.B.; Wang, B.; Ou, T.M.; Gu, L.Q.; Tan, J.H.; Huang, Z.S. Specific Targeting of Telomeric Multimeric G-Quadruplexes by a New Triaryl-Substituted Imidazole. Nucleic Acids Res. 2017, 45, 1606–1618. [Google Scholar] [CrossRef] [PubMed]
  12. Guiset Miserachs, H.; Donghi, D.; Börner, R.; Johannsen, S.; Sigel, R.K.O. Distinct Differences in Metal Ion Specificity of RNA and DNA G-Quadruplexes. J. Biol. Inorg. Chem. 2016, 21, 975–986. [Google Scholar] [CrossRef]
  13. Campbell, N.H.; Neidle, S. G-Quadruplexes and Metal Ions. In Metal Ions in Life Sciences, Interplay between Metal Ions and Nucleic Acids; Springer: Berlin/Heidelberg, Germany, 2012; pp. 119–134. ISBN 9789400721722. [Google Scholar]
  14. Smargiasso, N.; Hsia, W.; Colson, P.; Baker, E.S.; Bowers, M.T.; Pauw, E. De G-Quadruplex DNA Assemblies: Loop Length, Cation Identity, and Multimer Formation. J. Am. Chem. Soc. 2008, 130, 10208–10216. [Google Scholar] [CrossRef]
  15. Wang, Z.; Liu, J.P. Effects of the Central Potassium Ions on the G-Quadruplex and Stabilizer Binding. J. Mol. Graph. Model. 2017, 72, 168–177. [Google Scholar] [CrossRef] [PubMed]
  16. Bhattacharyya, D.; Arachchilage, G.M.; Basu, S. Metal Cations in G-Quadruplex Folding and Stability. Front. Chem. 2016, 4, 38. [Google Scholar] [CrossRef] [PubMed]
  17. Neidle, S. Quadruplex Nucleic Acids as Novel Therapeutic Targets. J. Med. Chem. 2016, 59, 5987–6011. [Google Scholar] [CrossRef] [PubMed]
  18. Neidle, S. The Structures of Quadruplex Nucleic Acids and Their Drug Complexes. Curr. Opin. Struct. Biol. 2009, 19, 239–250. [Google Scholar] [CrossRef]
  19. 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]
  20. Phan, A.T. Human Telomeric G-Quadruplex: Structures of DNA and RNA Sequences. FEBS J. 2010, 277, 1107–1117. [Google Scholar] [CrossRef]
  21. Cheong, C.; Moore, P.B. Solution Structure of an Unusually Stable RNA Tetraplex Containing G- and U-Quartet Structures. Biochemistry 1992, 31, 8406–8414. [Google Scholar] [CrossRef]
  22. Patel, P.K.; Koti, A.S.R.; Hosur, R. V NMR Studies on Truncated Sequences of Human Telomeric DNA: Observation of a Novel A-Tetrad. Nucleic Acids Res. 1999, 27, 3836–3843. [Google Scholar] [CrossRef] [PubMed]
  23. Patel, P.K.; Hosur, R. V NMR Observation of T-Tetrads in a Parallel Stranded DNA Quadruplex Formed by Saccharomyces cerevisiae Telomere Repeats. Nucleic Acids Res. 1999, 27, 2457–2464. [Google Scholar] [CrossRef]
  24. Patel, P.K.; Bhavesh, N.S.; Hosur, R.V. NMR Observation of a Novel C-Tetrad in the Structure of the SV40 Repeat Sequence GGGCGG. Biochem. Biophys. Res. Commun. 2000, 270, 967–971. [Google Scholar] [CrossRef] [PubMed]
  25. Pan, B.; Xiong, Y.; Shi, K.; Deng, J.; Sundaralingam, M.; Avenue, W.; Haven, N.; The, X. Crystal Structure of an RNA Purine-Rich Tetraplex Containing Adenine Tetrads: Implications for Specific Binding in RNA Tetraplexes. Structure 2003, 11, 815–823. [Google Scholar] [CrossRef] [PubMed]
  26. Pan, B.; Shi, K.; Sundaralingam, M.; Hall, J. Crystal Structure of an RNA Quadruplex Containing Inosine Tetrad: Implications for the Roles of NH2 Group in Purine Tetrads. J. Mol. Biol. 2006, 363, 451–459. [Google Scholar] [CrossRef]
  27. Kimura, T.; Xu, Y.; Komiyama, M. Human Telomeric RNA r(UAGGGU) Sequence Forms Parallel Tetraplex Structure with U-Qurtet. Nucleic Acids Symp. Ser. 2009, 53, 239–240. [Google Scholar] [CrossRef]
  28. Virgilio, A.; Esposito, V.; Citarella, G.; Mayol, L.; Galeone, A. Structural Investigations on the Anti-HIV G-Quadruplex-Forming Oligonucleotide TGGGAG and Its Analogues: Evidence for the Presence of an A-Tetrad. ChemBioChem 2012, 13, 2219–2224. [Google Scholar] [CrossRef]
  29. Xiao, C.-D.; Ishizuka, T.; Zhu, X.-Q.; Li, Y.; Sugiyama, H.; Xu, Y. An Unusual Topological RNA Architecture with an Eight—An Unusual Topological RNA Architecture with an Eight-Stranded Helical Fragment Containing A-, G- and U-Tetrads. J. Am. Chem. Soc. 2017, 139, 2565–2568. [Google Scholar] [CrossRef]
  30. Liu, H.; Wang, R.; Yu, X.; Shen, F.; Lan, W.; Haruehanroengra, P.; Yao, Q.; Zhang, J.; Chen, Y.; Li, S.; et al. High-Resolution DNA Quadruplex Structure Containing All the A-, G-, C-, T-Tetrads. Nucleic Acids Res. 2018, 46, 11627–11638. [Google Scholar] [CrossRef]
  31. Frasson, I.; Pirota, V.; Richter, S.N.; Doria, F. Multimeric G-Quadruplexes: A Review on Their Biological Roles and Targeting. Int. J. Biol. Macromol. 2022, 204, 89–102. [Google Scholar] [CrossRef]
  32. Alessandrini, I.; Recagni, M.; Zaffaroni, N.; Folini, M. On the Road to Fight Cancer: The Potential of G-quadruplex Ligands as Novel Therapeutic Agents. Int. J. Mol. Sci. 2021, 22, 5947. [Google Scholar] [CrossRef]
  33. Kosiol, N.; Juranek, S.; Brossart, P.; Heine, A.; Paeschke, K. G-Quadruplexes: A Promising Target for Cancer Therapy. Mol. Cancer 2021, 20, 1–18. [Google Scholar] [CrossRef]
  34. Liang, T.; Sun, X.; Li, W.; Hou, G.; Gao, F. 1,2,3-Triazole-Containing Compounds as Anti–Lung Cancer Agents: Current Developments, Mechanisms of Action, and Structure–Activity Relationship. Front. Pharmacol. 2021, 12, 661173. [Google Scholar] [CrossRef] [PubMed]
  35. Spiegel, J.; Adhikari, S.; Balasubramanian, S. The Structure and Function of DNA G-Quadruplexes. Trends Chem. 2020, 2, 123–136. [Google Scholar] [CrossRef] [PubMed]
  36. Kaneti, J.; Georgieva, M.; Rangelov, M.; Philipova, I.; Vasileva, B.; Angelov, I.; Staneva, D.; Miloshev, G.; Bakalova, S. Biological Activity of Quinazoline Analogues and Molecular Modeling of Their Interactions with G-Quadruplexes. Biochim. Biophys. Acta Gen. Subj. 2021, 1865, 129773–129787. [Google Scholar] [CrossRef] [PubMed]
  37. Kaneti, J.; Kurteva, V.; Georgieva, M.; Krasteva, N.; Miloshev, G.; Tabakova, N.; Petkova, Z.; Bakalova, S.M. Small Heterocyclic Ligands as Anticancer Agents: QSAR with a Model G-Quadruplex. Molecules 2022, 27, 7577. [Google Scholar] [CrossRef]
  38. Kircheva, N.; Dobrev, S.; Dasheva, L.; Nikolova, V.; Angelova, S.; Dudev, T. Metal-Assisted Complexation of Fluorogenic Dyes by Cucurbit[7]Uril and Cucurbit[8]Uril: A DFT Evaluation of the Key Factors Governing the Host–Guest Recognition. Molecules 2023, 28, 1540. [Google Scholar] [CrossRef]
  39. Kircheva, N.; Toshev, N.; Dudev, T. Holo-Chromodulin: Competition between the Native Cr3+ and Other Biogenic Cations (Fe3+, Fe2+, Mg2+, and Zn2+) for the Binding Sites. Metallomics 2022, 14, mfac082. [Google Scholar] [CrossRef]
  40. Nikolova, V.; Kircheva, N.; Dobrev, S.; Angelova, S. Lanthanides as Calcium Mimetic Species in Calcium-Signaling/Buffering Proteins: The Effect of Lanthanide Type on the Ca2+/Ln3+ Competition. Int. J. Mol. Sci. 2023, 24, 6297. [Google Scholar] [CrossRef]
  41. Vologzhannikova, A.A.; Shevelyova, M.P.; Kazakov, A.S.; Sokolov, A.S.; Borisova, N.I.; Permyakov, E.A.; Kircheva, N.; Nikolova, V.; Dudev, T.; Permyakov, S.E. Strontium Binding to α-Parvalbumin, a Canonical Calcium-Binding Protein of the “EF-Hand” Family. Biomolecules 2021, 11, 1158. [Google Scholar] [CrossRef]
  42. Kircheva, N.; Dobrev, S.; Nikolova, V.; Angelova, S.; Dudev, T. Zinc and Its Critical Role in Retinitis Pigmentosa: Insights from DFT/SMD Calculations. Inorg. Chem. 2020, 59, 17347–17355. [Google Scholar] [CrossRef] [PubMed]
  43. Wang, P.; Shi, R.; Su, Y.; Tang, L.; Huang, X.; Zhao, J. Hydrated Sodium Ion Clusters [Na+(H2O)n (n = 1–6)]: An Ab Initio Study on Structures and Non-Covalent Interaction. Front. Chem. 2019, 7, 624. [Google Scholar] [CrossRef] [PubMed]
  44. Rodriguez-Cruz, S.E.; Jockusch, R.A.; Williams, E.R. Hydration Energies and Structures of Alkaline Earth Metal Ions, M2+(H2O)n, n = 5–7, M = Mg, Ca, Sr, and Ba. J. Am. Chem. Soc. 1999, 121, 8898–8906. [Google Scholar] [CrossRef] [PubMed]
  45. Carrillo-Tripp, M.; Saint-Martin, H.; Ortega-Blake, I. A Comparative Study of the Hydration of Na+ and K+ with Refined Polarizable Model Potentials. J. Chem. Phys. 2003, 118, 7062–7073. [Google Scholar] [CrossRef]
  46. Bush, M.F.; Saykally, R.J.; Williams, E.R. Hydration of the Calcium Dication: Direct Evidence for Second Shell Formation from Infrared Spectroscopy. ChemPhysChem 2007, 8, 2245–2253. [Google Scholar] [CrossRef]
  47. Kircheva, N.; Dobrev, S.; Dasheva, L.; Koleva, I.; Nikolova, V.; Angelova, S.; Dudev, T. Complexation of Biologically Essential (Mono- And Divalent) Metal Cations to Cucurbiturils: A DFT/SMD Evaluation of the Key Factors Governing the Host-Guest Recognition. RSC Adv. 2020, 10, 28139–28147. [Google Scholar] [CrossRef]
  48. Tahoon, M.A.; Gomaa, E.A.; Suleiman, M.H.A. Aqueous Micro-Hydration of Na+(H2O)n = 1–7 Clusters: DFT Study. Open Chem. 2019, 17, 260–269. [Google Scholar]
  49. Jia, N.; Wang, P.; Su, Y. Structures and Spectroscopic Properties of K+(H2O)n with n = 1–10 Clusters Based on Density Functional Theory. Chem. Phys. Lett. 2022, 801, 139735. [Google Scholar] [CrossRef]
  50. Boda, A.; De, S.; Ali, S.M.; Tulishetti, S.; Khan, S.; Singh, J.K. From Microhydration to Bulk Hydration of Sr2+ Metal Ion: DFT, MP2 and Molecular Dynamics Study. J. Mol. Liq. 2012, 172, 110–118. [Google Scholar] [CrossRef]
  51. Collie, G.W.; Haider, S.M.; Neidle, S.; Parkinson, G.N. A Crystallographic and Modelling Study of a Human Telomeric RNA (TERRA) Quadruplex. Nucleic Acids Res. 2010, 38, 5569–5580. [Google Scholar] [CrossRef]
  52. 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] [PubMed]
  53. Alberts, B.; Johnson, A.; Lewis, J.; Raff, M.; Roberts, K.; Walter, P. Molecular Biology of the Cell, 4th ed.; Garland Science: New York, NY, USA, 2008. [Google Scholar]
  54. Shannon, R.D. Revised Effective Ionic Radii and Systematic Studies of Interatomic Distances in Halides and Chalcogenides. Acta Crystallogr. Sect. A 1976, 32, 751–767. [Google Scholar] [CrossRef]
  55. Hirshfeld, F.L. Bonded-Atom Fragments for Describing Molecular Charge Densities. Theor. Chim. Acta 1977, 44, 129–138. [Google Scholar] [CrossRef]
  56. Liu, X.; O’Harra, K.E.; Bara, J.E.; Turner, C.H. Screening Ionic Liquids Based on Ionic Volume and Electrostatic Potential Analyses. J. Phys. Chem. B 2021, 125, 3653–3664. [Google Scholar] [CrossRef]
  57. Frisch, M.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R. Gaussian 09, Revision D.01, Gaussian 2013; Gaussian, Inc.: Wallingford, CT, USA, 2009. [Google Scholar]
  58. Grimme, S. Improved Second-Order Møller-Plesset Perturbation Theory by Separate Scaling of Parallel- and Antiparallel-Spin Pair Correlation Energies. J. Chem. Phys. 2003, 118, 9095–9102. [Google Scholar] [CrossRef]
  59. Chai, J.D.; Head-Gordon, M. Long-Range Corrected Hybrid Density Functionals with Damped Atom-Atom Dispersion Corrections. Phys. Chem. Chem. Phys. 2008, 10, 6615–6620. [Google Scholar] [CrossRef]
  60. Andrae, D.; Häußermann, U.; Dolg, M.; Stoll, H.; Preuß, H. Energy-Adjusted Ab Initio Pseudopotentials for the Second and Third Row Transition Elements. Theor. Chim. Acta 1990, 77, 123–141. [Google Scholar] [CrossRef]
  61. Marenich, A.V.; Cramer, C.J.; Truhlar, D.G. Universal Solvation Model Based on Solute Electron Density and on a Continuum Model of the Solvent Defined by the Bulk Dielectric Constant and Atomic Surface Tensions. J. Phys.Chem. B 2009, 113, 6378–6396. [Google Scholar] [CrossRef]
  62. Schrödinger, L.; DeLano, W. The PyMol Graphics System, Version 2.2.3; Schrödinger, Inc.: New York, NY, USA, 2018. Available online: http://www.pymol.org/pymol(accessed on 10 August 2023).
  63. Kircheva, N.; Angelova, S.; Dobrev, S.; Petkova, V.; Nikolova, V.; Dudev, T. Cu+/Ag+ Competition in Type I Copper Proteins (T1Cu). Biomolecules 2023, 13, 681. [Google Scholar] [CrossRef]
Figure 1. Guanine structures constructed via HG base pairs from guanine: 2D model of a G-quartet (A) and schematic structure of a quadruplex (B); function of quadruplexes (C).
Figure 1. Guanine structures constructed via HG base pairs from guanine: 2D model of a G-quartet (A) and schematic structure of a quadruplex (B); function of quadruplexes (C).
Molecules 28 06109 g001
Figure 2. Structures of the fully optimized Na(H2O)6Ⴈ+, K(H2O)6Ⴈ+, and Sr(H2O)8Ⴈ2+ constructs at the wB97XD/6-31G(d,p) level used in the current study (A); juxtaposed PDB deposited (3IBK) and GG quadruplex models, fully optimized in the gas phase, colored in light blue and purple, correspondingly (B).
Figure 2. Structures of the fully optimized Na(H2O)6Ⴈ+, K(H2O)6Ⴈ+, and Sr(H2O)8Ⴈ2+ constructs at the wB97XD/6-31G(d,p) level used in the current study (A); juxtaposed PDB deposited (3IBK) and GG quadruplex models, fully optimized in the gas phase, colored in light blue and purple, correspondingly (B).
Molecules 28 06109 g002
Figure 3. Gibbs energies of substitution of K+ with Na+ (in kcal mol−1) in three environments of different polarity for wB97XD/6-31+G(d,p)//wB97XD/6-31G(d,p) optimized M@GQႨ+ complexes: M@GGႨ+ (A); M@AGႨ+ (B); M@CGႨ+ (C); M@TGႨ+ (D); M@UGႨ+ (E). The color scheme is further applied in all presented figures.
Figure 3. Gibbs energies of substitution of K+ with Na+ (in kcal mol−1) in three environments of different polarity for wB97XD/6-31+G(d,p)//wB97XD/6-31G(d,p) optimized M@GQႨ+ complexes: M@GGႨ+ (A); M@AGႨ+ (B); M@CGႨ+ (C); M@TGႨ+ (D); M@UGႨ+ (E). The color scheme is further applied in all presented figures.
Molecules 28 06109 g003
Figure 4. Gibbs energies of substitution of K+ with Sr2+ (in kcal mol−1) in three environments of different polarity for wB97XD/6-31+G(d,p)//wB97XD/6-31G(d,p) optimized M@GQႨ+/2+ complexes: M@GGႨ+/2+ (A); M@AGႨ+/2+ (B); M@CGႨ+/2+ (C); M@TGႨ+/2+ (D); M@UGႨ+/2+ (E).The applied color scheme has been previously presented in Figure 3.
Figure 4. Gibbs energies of substitution of K+ with Sr2+ (in kcal mol−1) in three environments of different polarity for wB97XD/6-31+G(d,p)//wB97XD/6-31G(d,p) optimized M@GQႨ+/2+ complexes: M@GGႨ+/2+ (A); M@AGႨ+/2+ (B); M@CGႨ+/2+ (C); M@TGႨ+/2+ (D); M@UGႨ+/2+ (E).The applied color scheme has been previously presented in Figure 3.
Molecules 28 06109 g004
Figure 5. Gibbs energies of formation (in kcal mol−1) in three environments of different polarity for wB97XD/6-31+G(d,p)//wB97XD/6-31G(d,p)optimized L@K@GQ Ⴈ+ ternary complexes: L@M@GGႨ+ (A); L@M@AGႨ+ (B); L@M@CGႨ+ (C); L@M@TGႨ+ (D); L@M@UGႨ+ (E). The applied color scheme has been previously presented in Figure 3.
Figure 5. Gibbs energies of formation (in kcal mol−1) in three environments of different polarity for wB97XD/6-31+G(d,p)//wB97XD/6-31G(d,p)optimized L@K@GQ Ⴈ+ ternary complexes: L@M@GGႨ+ (A); L@M@AGႨ+ (B); L@M@CGႨ+ (C); L@M@TGႨ+ (D); L@M@UGႨ+ (E). The applied color scheme has been previously presented in Figure 3.
Molecules 28 06109 g005
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kircheva, N.; Dobrev, S.; Petkova, V.; Bakalova, S.; Kaneti, J.; Angelova, S. Theoretical Assessment of the Ligand/Metal/Quadruplex Recognition in the Non-Canonical Nucleic Acids Structures. Molecules 2023, 28, 6109. https://doi.org/10.3390/molecules28166109

AMA Style

Kircheva N, Dobrev S, Petkova V, Bakalova S, Kaneti J, Angelova S. Theoretical Assessment of the Ligand/Metal/Quadruplex Recognition in the Non-Canonical Nucleic Acids Structures. Molecules. 2023; 28(16):6109. https://doi.org/10.3390/molecules28166109

Chicago/Turabian Style

Kircheva, Nikoleta, Stefan Dobrev, Vladislava Petkova, Snezhana Bakalova, Jose Kaneti, and Silvia Angelova. 2023. "Theoretical Assessment of the Ligand/Metal/Quadruplex Recognition in the Non-Canonical Nucleic Acids Structures" Molecules 28, no. 16: 6109. https://doi.org/10.3390/molecules28166109

APA Style

Kircheva, N., Dobrev, S., Petkova, V., Bakalova, S., Kaneti, J., & Angelova, S. (2023). Theoretical Assessment of the Ligand/Metal/Quadruplex Recognition in the Non-Canonical Nucleic Acids Structures. Molecules, 28(16), 6109. https://doi.org/10.3390/molecules28166109

Article Metrics

Back to TopTop