Effect of Coordinated Ions on Structure and Flexibility of Parallel G-quadruplexes: A Molecular Dynamics Study

Abstract Single tract guanine residues can associate to form stable parallel quadruplex structures in the presence of certain cations. Nanosecond scale molecular dynamics simulations have been performed on fully solvated fibre model of parallel d(G(7)) quadruplex structures with Na(+) or K(+) ions coordinated in the cavity formed by the O6 atoms of the guanine bases. The AMBER 4.1 force field and Particle Mesh Ewald technique for electrostatic interactions have been used in all simulations. These quadruplex structures are stable during the simulation, with the middle four base tetrads showing root mean square deviation values between 0.5 to 0.8 ? from the initial structure as well the high resolution crystal structure. Even in the absence of any coordinated ion in the initial structure, the G-quadruplex structure remains intact throughout the simulation. During the 1.1 ns MD simulation, one Na(+) counter ion from the solvent as well as several water molecules enter the central cavity to occupy the empty coordination sites within the parallel quadruplex and help stabilize the structure. Hydrogen bonding pattern depends on the nature of the coordinated ion, with the G-tetrad undergoing local structural variation to accommodate cations of different sizes. In the absence of any coordinated ion, due to strong mutual repulsion, O6 atoms within G-tetrad are forced farther apart from each other, which leads to a considerably different hydrogen bonding scheme within the G-tetrads and very favourable interaction energy between the guanine bases constituting a G-tetrad. However, a coordinated ion between G-tetrads provides extra stacking energy for the G-tetrads and makes the quadruplex structure more rigid. Na(+) ions, within the quadruplex cavity, are more mobile than coordinated K(+) ions. A number of hydrogen bonded water molecules are observed within the grooves of all quadruplex structures.     

Effect of coordinated ions on structure and flexiblity of parallel G-quandruplexes: a molecular dynamics study

Single tract guanine residues can associate to form stable parallel quadruplex structures in the presence of certain cations. Nanosecond scale molecular dynamics simulations have been performed on fully solvated fibre model of parallel d(G7) quadruplex structures with Na+ or K+ ions coordinated in the cavity formed by the 06 atoms of the guanine bases. The AMBER 4.1 force field and Particle Mesh Ewald technique for electrostatic interactions have been used in all simulations. These quadruplex structures are stable during the simulation, with the middle four base tetrads showing root mean square deviation values between 0.5 to 0.8 A from the initial structure as well the high resolution crystal structure. Even in the absence of any coordinated ion in the initial structure, the G-quadruplex structure remains intact throughout the simulation. During the 1.1 ns MD simulation, one Na+ counter ion from the solvent as well as several water molecules enter the central cavity to occupy the empty coordination sites within the parallel quadruplex and help stabilize the structure. Hydrogen bonding pattern depends on the nature of the coordinated ion, with the G-tetrad undergoing local structural variation to accommodate cations of different sizes. In the absence of any coordinated ion, due to strong mutual repulsion, 06 atoms within G-tetrad are forced farther apart from each other, which leads to a considerably different hydrogen bonding scheme within the G-tetrads and very favourable interaction energy between the guanine bases constituting a G-tetrad. However, a coordinated ion between G-tetrads provides extra stacking energy for the G-tetrads and makes the quadruplex structure more rigid. Na+ ions, within the quadruplex cavity, are more mobile than coordinated K+ ions. A number of hydrogen bonded water molecules are observed within the grooves of all quadruplex structures.     

A Nanosecond Molecular Dynamics Study of Antiparallel d(G)7 Quadruplex Structures: Effect of the Coordinated Cations

Abstract

Nanosecond scale molecular dynamics simulations have been performed on antiparallel Greek key type d(G₇) quadruplex structures with different coordinated ions, namely Na⁺ and K⁺ ion, water and Na⁺ counter ions, using the AMBER force field and Particle Mesh Ewald technique for electrostatic interactions. Antiparallel structures are stable during the simulation, with root mean square deviation values of ～ 1.5 Å from the initial structures. Hydrogen bonding patterns within the G-tetrads depend on the nature of the coordinated ion, with the G-tetrad undergoing local structural variation to accommodate different cations. However, alternating syn-anti arrangement of bases along a chain as well as in a quartet is maintained through out the MD simulation. Coordinated Na+ ions, within the quadruplex cavity are quite mobile within the central channel and can even enter or exit from the quadruplex core, whereas coordinated K⁺ ions are quite immobile. MD studies at 400K indicate that K⁺ ion cannot come out from the quadruplex core without breaking the terminal G-tetrads. Smaller grooves in antiparallel structures are better binding sites for hydrated counter ions, while a string of hydrogen bonded water molecules are observed within both the small and large grooves. The hydration free energy for the K⁺ ion coordinated structure is more favourable than that for the Na⁺ ion coordinated antiparallel quadruplex structure.     

A nanosecond molecular dynamics study of antiparallel d(G)7 quadruplex structures: effect of the coordinated cations

Nanosecond scale molecular dynamics simulations have been performed on antiparallel Greek key type d(G7) quadruplex structures with different coordinated ions, namely Na+ and K+ ion, water and Na+ counter ions, using the AMBER force field and Particle Mesh Ewald technique for electrostatic interactions. Antiparallel structures are stable during the simulation, with root mean square deviation values of approximately 1.5 A from the initial structures. Hydrogen bonding patterns within the G-tetrads depend on the nature of the coordinated ion, with the G-tetrad undergoing local structural variation to accommodate different cations. However, alternating syn-anti arrangement of bases along a chain as well as in a quartet is maintained through out the MD simulation. Coordinated Na+ ions, within the quadruplex cavity are quite mobile within the central channel and can even enter or exit from the quadruplex core, whereas coordinated K+ ions are quite immobile. MD studies at 400K indicate that K+ ion cannot come out from the quadruplex core without breaking the terminal G-tetrads. Smaller grooves in antiparallel structures are better binding sites for hydrated counter ions, while a string of hydrogen bonded water molecules are observed within both the small and large grooves. The hydration free energy for the K+ ion coordinated structure is more favourable than that for the Na+ ion coordinated antiparallel quadruplex structure.     

Interaction of Cyclic Cytosine-, Guanine-, Thymine-, Uracil- and Mixed Guanine-Cytosine Base Tetrads with K+, Na+ and Li+ Ions—A Density Functional Study

Abstract

We have carried out B3LYP hybrid density functional studies of complexes formed by cyclic cytosine-, guanine-, thymine-, uracil- and mixed guanine cytosine-tetrads with Li⁺, Na⁺ and K⁺ ions to determine their structures and interaction energies. The conformations studied have been restricted to a hydrogen bond pattern closely related to the tetrads observed in experimental nucleic acid structures. A comparison of the alkali metal ion/tetrad complexes with the tetrads without cations indicates that alkali metal ions modulate the tetrad structures significantly and that even the hydrogen bond pattern may change. Guanine-tetrad cation complexes show the strongest interaction energy compared to other tetrads that occur less frequently in experimental structures. The most stable G-tetrad/metal ion structure adopts a nearly planar geometry that is especially suitable for tetraplex formation, which requires approximately parallel tetrad planes. In the cytosine-tetrad there is a very large central cavity suitable for cation recognition, but the complexes adopt a non-planar structure unsuitable for stacking, except possibly for ions with very large radii. Uracil and thymine tetrads show a significant different characteristics which may contribute to the differences between DNA and RNA.     

A COMPARISON OF THE PROPERTIES AND THE SOLUTION STRUCTURE FOR RNA AND DNA QUADRUPLEXES WHICH CONTAIN TWO GGAGG SEQUENCES JOINED WITH A TETRANUCLEOTIDE LINKER

Abstract

We have determined solution structure of r(GGAGGUUUUGGAGG) (R14) by NMR; the RNA 14-mer forms an intra-strand parallel quadruplex with a G-tetrad and a hexad, in which a G-tetrad core is augmented by association of two A residues. The quadruplex further forms a dimer through stacking interaction between the hexads. In order to obtain insight into the difference between RNA and DNA quadruplexes, we synthesized the corresponding DNA 14-mer, d(GGAGGTTTTGGAGG) (D14), and examined its properties and structure by CD, gel electrophoresis, and NMR. K⁺ ions increased the thermal stability of both R14 and D14 structures. The binding affinity of K⁺ ions to R14 was much higher than that to D14. The CD and gel electrophoretic studies suggest that D14 forms a quadruplex entirely different from that of R14 in the presence of K⁺ ions; two molecules of D14 form a quadruplex with both antiparallel and parallel strand alignments and with diagonal loops at both ends of the stacked G-tetrads. The NMR study also gave results that are consistent with such structure: alternate glycosidic conformation, 5′G(syn)-G(anti)3′, and characteristic chemical shift data observed for many quadruplexes containing diagonal TTTT loops.     

A parallel stranded G‐quadruplex composed of threose nucleic acid (TNA)

G‐rich sequences can adopt four‐stranded helical structures, called G‐quadruplexes, that self‐assemble around monovalent cations like sodium (Na⁺) and potassium (K⁺). Whether similar structures can be formed from xeno‐nucleic acid (XNA) polymers with a shorter backbone repeat unit is an unanswered question with significant implications on the fold space of functional XNA polymers. Here, we examine the potential for TNA (α‐l ‐threofuranosyl nucleic acid) to adopt a four‐stranded helical structure based on a planar G‐quartet motif. Using native polyacrylamide gel electrophoresis (PAGE), circular dichroism (CD) and solution‐state nuclear magnetic resonance (NMR) spectroscopy, we show that despite a backbone repeat unit that is one atom shorter than the backbone repeat unit found in DNA and RNA, TNA can self‐assemble into stable G‐quadruplex structures that are similar in thermal stability to equivalent DNA structures. However, unlike DNA, TNA does not appear to discriminate between Na⁺ and K⁺ ions, as G‐quadruplex structures form equally well in the presence of either ion. Together, these findings demonstrate that despite a shorter backbone repeat unit, TNA is capable of self‐assembling into stable G‐quadruplex structures.     

A comparison of the properties and the solution structure for RNA and DNA quadruplexes which contain two GGAGG sequences joined with a tetranucleotide linker

We have determined solution structure of r(GGAGGUUUUGGAGG) (R14) by NMR; the RNA 14-mer forms an intra-strand parallel quadruplex with a G-tetrad and a hexad, in which a G-tetrad core is augmented by association of two A residues. The quadruplex further forms a dimer through stacking interaction between the hexads. In order to obtain insight into the difference between RNA and DNA quadruplexes, we synthesized the corresponding DNA 14-mer, d(GGAGGTTTTGGAGG) (D14), and examined its properties and structure by CD, gel electrophoresis, and NMR. K+ ions increased the thermal stability of both R14 and D14 structures. The binding affinity of K+ ions to R14 was much higher than that to D14. The CD and gel electrophoretic studies suggest that D14 forms a quadruplex entirely different from that of R14 in the presence of K+ ions; two molecules of D14 form a quadruplex with both antiparallel and parallel strand alignments and with diagonal loops at both ends of the stacked G-tetrads. The NMR study also gave results that are consistent with such structure: alternate glycosidic conformation, 5'G(syn)-G(anti)3', and characteristic chemical shift data observed for many quadruplexes containing diagonal TTTT loops.     

Conformational switch of insulin-binding aptamer into G-quadruplex induced by K+ and Na+: an experimental and theoretical approach

Guanine-rich sequences can form the G-quadruplex structure in the presence of specific metal ions. Here, circular dichroism, UV–vis absorption, fluorescence, and molecular dynamics simulation studies revealed that insulin-binding aptamer (IBA) could form an intramolecular G-quadruplex structure after binding K⁺. Circular dichroism (CD) spectra demonstrated that IBA could fold into a parallel G-quadruplex with a strong positive peak at 263?nm. Analysis of equilibrium titration data revealed that cation binding was cooperative with the Hill coefficient of 2.01 in K⁺ and 1.90 in Na⁺. Thermal denaturation assays indicated that K⁺-induced G-quadruplex is more stable than Na⁺-induced structure. Folding of IBA into G-quadruplex leading to the contact quenching occurs as a result of the formation of a nonfluorescent complex between donor and acceptor. Based on fluorescence quenching of IBA folding, a potassium-sensing aptasensor in the range of 0–1.4?mM was proposed. Since the quenching process was predominantly static, the binding constant and the number of binding sites were determined. In this research, based on the experimental data, the initial model of IBA G-quadruplex was constructed by molecular modeling method. The modeling structure of IBA is an intramolecular parallel-strand quadruplex conformation with two guanine tetrads. The extended molecular dynamics simulation for the model indicated that the G-quadruplex maintains its structure very well in aqueous solution in presence of K⁺ in the central cavity. In contrast, it was demonstrated that the G-quadruplex structure of model in the water collapses in absence of this cation.     

Interaction Energies and Dynamics of Alkali and Alkaline-Earth Cations in Quadruplex-DNA-Structures

We have investigated the nonbonded interaction energies and dynamical properties of different types of cations in quadruplex DNA structures using the GROMOS force field [1]. Quadruplex structures consist of planar guanine-quartets stacking together and causing the formation of a channel, large enought to enclose several cations (Figure 1). In recent years many experimental studies have indicated a prefered formation of this unusually stabel complexes with K⁺-ions. However, the high selectivity of this cation has not yet been understood [2].To determine the most stable coordination sites and the mobility of cations, we have calculated the pair potential energy of alkali and alkaline-earth cations along the helical axis of a model quadruplex structure (Figure 2). Our force field calculations indicate that small ions like Li⁺, Na⁺, Mg²⁺ and Ca²⁺ are free to move throught the channel. In contrast, for K⁺ and larger ions a high potential barrier appears, located in the plane of the tetramer unit. These findings are in agreement with data from X-ray crystallography, indicating that K⁺ cations are located between two planes while Na⁺ ions also can occupy coordination sites in the G-quartet plane.Considering solvent atoms in our calculations leads to the observation that a cation at the end of the quadruplex strongly interacts with one water molecule located near the entrance of the cage. Snapshots taken at different times of the MD simulation provide configurations which differ mainly in the position of this complexing water molecule. Moving away from the entrance of the cage causes a significant decrease of the potential barrier for K⁺ and smaller cage cations (Figure 3). For the larger ions the potential barrier is much higher than the thermal energy (not shown), preventing the cations from leaving and entering the cage.This conclusion is in agreement with results from our MD simulations. We followed the dynamics of different cations. While K⁺ is able to leave as well as to re-migrate into the channel (movie I), this was not observed for other types of cations. Figure 4 shows the time history of the positional fluctuations of potassium along the helical axis, and in Figure 5 we have monitored the distance between the O6-oxygen atoms of the outer G-quartet. It becomes clearly evident that the cation movement through the planes is correlated with the dynamic behaviour of the tetrameric planes. When the K⁺-ion penetrates the tetrameric unit to enter the quadruplex, the O6-O6-distance - a measure of size of the hole of the plane - increases. After a while the cation has reached the cage position and the G-quartet contracts to the initial value. That means the tetrameric planes perform a kind of breathing motion.For lithium ions we find a much higher mobility of the cation within the quadruplex channel. Two of three ions are leaving the cage instantaniously (not shown).Another indication of the experimentally observed much weaker complexation tendency of quadruplexes with lithium is the change in the distances of planes (a measure for the cage size) with time. While in the case of potassium the distance of planes is nearly the same for all three cages, for lithium the central occupied cage is much smaller than the cage in the starting structure, indicating that the DNA structure has to adjust its conformation to the cation size. On the other side the outer unoccupied cages are much greater and less stable. Due to this cation induced quadruplex deformation we observe an unwinding of the DNA-structure in the presens of lithium ions at longer simulation periods (400 ps).     

Arrangements of human telomere DNA quadruplex in physiologically relevant K+ solutions

The arrangement of the human telomeric quadruplex in physiologically relevant conditions has not yet been unambiguously determined. Our spectroscopic results suggest that the core quadruplex sequence G₃(TTAG₃)₃ forms an antiparallel quadruplex of the same basket type in solution containing either K⁺ or Na⁺ ions. Analogous sequences extended by flanking nucleotides form a mixture of the antiparallel and hybrid (3 + 1) quadruplexes in K⁺-containing solutions. We, however, show that long telomeric DNA behaves in the same way as the basic G₃(TTAG₃)₃ motif. Both G₃(TTAG₃)₃ and long telomeric DNA are also able to adopt the (3 + 1) quadruplex structure: Molecular crowding conditions, simulated here by ethanol, induced a slow transition of the K⁺-stabilized quadruplex into the hybrid quadruplex structure and then into a parallel quadruplex arrangement at increased temperatures. Most importantly, we demonstrate that the same transitions can be induced even in aqueous, K⁺-containing solution by increasing the DNA concentration. This is why distinct quadruplex structures were detected for AG₃(TTAG₃)₃ by X-ray, nuclear magnetic resonance and circular dichrosim spectroscopy: Depending on DNA concentration, the human telomeric DNA can adopt the antiparallel quadruplex, the (3 + 1) structure, or the parallel quadruplex in physiologically relevant concentrations of K⁺ ions.     

A thymine tetrad in d(TGGGGT) quadruplexes stabilized with Tl+/Na+ ions

We report two new structures of the quadruplex d(TGGGGT)₄ obtained by single crystal X-ray diffraction. In one of them a thymine tetrad is found. Thus the yeast telomere sequences d(TG_1–3) might be able to form continuous quadruplex structures, involving both guanine and thymine tetrads. Our study also shows substantial differences in the arrangement of thymines when compared with previous studies. We find five different types of organization: (i) groove binding with hydrogen bonds to guanines from a neighbour quadruplex; (ii) partially ordered groove binding, without any hydrogen bond; (iii) stacked thymine triads, formed at the 3′ends of the quadruplexes; (iv) a thymine tetrad between two guanine tetrads. Thymines are stabilized in pairs by single hydrogen bonds. A central sodium ion interacts with two thymines and contributes to the tetrad structure. (v) Completely disordered thymines which do not show any clear location in the crystal. The tetrads are stabilized by either Na⁺ or Tl⁺ ions. We show that by using MAD methods, Tl⁺ can be unambiguously located and distinguished from Na⁺. We can thus determine the preference for either ion in each ionic site of the structure under the conditions used by us.     

Small-Molecule Selectively Recognizes Human Telomeric G-Quadruplex DNA and Regulates Its Conformational Switch

Structural complexity is an inherent feature of the human telomeric sequence, and it presents a major challenge for developing ligands of pharmaceutical interest. Recent studies have pointed out that the induction of a quadruplex or change of a quadruplex conformation on binding may be the most powerful method to exert the desired biological effect. In this study, we demonstrate a quadruplex ligand that binds selectively to different forms of the human telomeric G-quadruplex structure and regulates its conformational switch. The results show that not only can oxazine750 selectively induce parallel quadruplex formation from a random coil telomeric oligonucleotide in the absence of added cations, it also can easily surpass the energy barrier between two structures and change the G-quadruplex conformation in Na⁺ or K⁺ solution. The combination of its unique properties, including the size and shape of the G-quadruplex and the small molecule, is proposed as the predominant force for regulating the special structural formation and transitions. These results may stimulate the design of new quadruplex binders that would be capable of discriminating different G-quadruplex structures as well as controlling biological phenomena, functional molecules, and nanomaterials.     

Interaction of cyclic cytosine-, guanine-, thymine-, uracil- and mixed guanine-cytosine base tetrads with K+, Na+ and Li+ ions -- a density functional study

We have carried out B3LYP hybrid density functional studies of complexes formed by cyclic cytosine-, guanine-, thymine-, uracil- and mixed guanine cytosine-tetrads with Li+, Na+ and K+ ions to determine their structures and interaction energies. The conformations studied have been restricted to a hydrogen bond pattern closely related to the tetrads observed in experimental nucleic acid structures. A comparison of the alkali metal ion/tetrad complexes with the tetrads without cations indicates that alkali metal ions modulate the tetrad structures significantly and that even the hydrogen bond pattern may change. Guanine-tetrad cation complexes show the strongest interaction energy compared to other tetrads that occur less frequently in experimental structures. The most stable G-tetrad/metal ion structure adopts a nearly planar geometry that is especially suitable for tetraplex formation, which requires approximately parallel tetrad planes. In the cytosine-tetrad there is a very large central cavity suitable for cation recognition, but the complexes adopt a non-planar structure unsuitable for stacking, except possibly for ions with very large radii. Uracil and thymine tetrads show a significant different characteristics which may contribute to the differences between DNA and RNA     

Drug recognition and stabilisation of the parallel-stranded DNA quadruplex d(TTAGGGT)4 containing the human telomeric repeat

The NMR structure of the parallel-stranded DNA quadruplex d(TTAGGGT)(4), containing the human telomeric repeat, has been determined in solution in complex with a fluorinated pentacyclic quino[4,3,2-kl]acridinium cation (RHPS4). RHPS4 has been identified as a potent inhibitor of telomerase at submicromolar levels (IC(50) value of 0.33(+/-0.13)microM), exhibiting a wide differential between telomerase inhibition and acute cellular toxicity. All of the data point to RHPS4 exerting its chemotherapeutic potency through interaction with, and stabilisation of, four-stranded G-quadruplex structures. RHPS4 forms a dynamic interaction with d(TTAGGGT)(4), as evident from 1H and 19F linewidths, with fast exchange between binding sites induced at 318 K. Perturbations to DNA chemical shifts and 24 intermolecular nuclear Overhauser effects (NOEs) identify the 5'-ApG and 5'-GpT steps as the principle intercalation sites; a structural model has been refined using NOE-restrained molecular dynamics. The central G-tetrad core remains intact, with drug molecules stacking at the ends of the G-quadruplex. The partial positive charge on position 13-N of the acridine ring appears to act as a "pseudo" potassium ion and is positioned above the centre of the G-tetrad in the region of high negative charge density. In both ApG and GpT intercalation sites, the drug is seen to converge to the same orientation in which the pi-system of the drug overlaps primarily with two bases of each G-tetrad. The drug is held in place by stacking interactions with the G-tetrads; however, there is some evidence for a more dynamic, weakly stabilised A-tetrad that stacks partially on top of the drug at the 5'-end of the sequence. Together, the interactions of RHPS4 increase the t(m) of the quadruplex by approximately 20 degrees C. There is no evidence for drug intercalation within the G-quadruplex; however, the structural model strongly supports end-stacking interactions with the terminal G-tetrads.     

Distinct differences in metal ion specificity of RNA and DNA G-quadruplexes

RNA G-quadruplexes, as their well-studied DNA analogs, require the presence of cations to fold and remain stable. This is the first comprehensive study on the interaction of RNA quadruplexes with metal ions. We investigated the formation and stability of two highly conserved and biologically relevant RNA quadruplex-forming sequences (24nt-TERRA and 18nt-NRAS) in the presence of several monovalent and divalent metal ions, namely Li⁺, Na⁺, K⁺, Rb⁺, Cs⁺, NH₄ ⁺, Mg²⁺, Ca²⁺, Sr²⁺, and Ba²⁺. Circular dichroism was used to probe the influence of these metal ions on the folded fraction of the parallel G-quadruplexes, and UV thermal melting experiments allowed to assess the relative stability of the structures in each cationic condition. Our results show that the RNA quadruplexes are more stable than their DNA counterparts under the same buffer conditions. We have observed that the addition of mainly Na⁺, K⁺, Rb⁺, NH₄ ⁺, as well as Sr²⁺ and Ba²⁺ in water, shifts the equilibrium to the folded quadruplex form, whereby the NRAS sequence responds stronger than TERRA. However, only K⁺ and Sr²⁺ lead to a significant increase in the stability of the folded structures, which is consistent with their coordination to the O6 atoms from the G-quartet guanosines. Compared to the respective DNA motives, dNRAS and htelo, the RNA sequences are not stabilized by Na⁺ ions. Finally, the difference in response between NRAS and TERRA, as well as to the corresponding DNA sequences with respect to different metal ions, could potentially be exploited for selective targeting purposes.     

The formation pathway of tetramolecular G-quadruplexes

Oligonucleotides containing guanosine stretches associate into tetrameric structures stabilized by monovalent ions. In order to describe the sequence of reactions leading to association of four identical strands, we measured by NMR the formation and dissociation rates of (TGnT)₄ quadruplexes (n = 3–6), their dissociation constants and the reaction orders for quadruplex formation. The quadruplex formation rates increase with the salt concentration but weakly depend on the nature (K⁺, Na⁺ or Li⁺) of the counter ions. The activation energies for quadruplex formation are negative. The quadruplex lifetimes strongly increase with the G-tract length and are much more longer in K⁺ solution than in Na⁺ or Li⁺ solutions. The reaction order for quadruplex formation is 3 in 0.125 M KCl and 4 in LiCl solutions. The kinetics measurements suggest that quadruplex formation proceeds step by step via sequential strand association into duplex and triplex intermediate species. Triplex formation is rate limiting in 0.125 M KCl solution. In LiCl, each step of the association process depends on the strand concentration. Parallel reactions to formation of the fully matched canonical quadruplex may result in kinetically trapped mismatched quadruplexes making the canonical quadruplex practically inaccessible in particular at low temperature in KCl solution.     

Molecular dynamics simulation of hepatitis C virus IRES IIId domain: structural behavior,electrostatic and energetic analysis

The dynamic behavior of the HCV IRES IIId domain is analyzed by means of a 2.6-ns molecular dynamics simulation, starting from an NMR structure. The simulation is carried out in explicit water with Na⁺ counterions, and particle-mesh Ewald summation is used for the electrostatic interactions. In this work, we analyze selected patterns of the helix that are crucial for IRES activity and that could be considered as targets for the intervention of inhibitors, such as the hexanucleotide terminal loop (more particularly its three consecutive guanines) and the loop-E motif. The simulation has allowed us to analyze the dynamics of the loop substructure and has revealed a behavior among the guanine bases that might explain the different role of the third guanine of the GGG triplet upon molecular recognition. The accessibility of the loop-E motif and the loop major and minor groove is also examined, as well as the effect of Na⁺ or Mg²⁺ counterion within the simulation. The electrostatic analysis reveals several ion pockets, not discussed in the experimental structure. The positions of these ions are useful for locating specific electrostatic recognition sites for potential inhibitor binding. Figure Superposition of 14 structures representative of the evolution of IRES IIId RNA along 2.6-ns MD simulation