Structural dynamics and cation interactions of DNA quadruplex molecules containing mixed guanine/cytosine quartets revealed by large-scale MD simulations

Computational and Experimental Characterization of rDNA and rRNA G-Quadruplexes

DNA G-quadruplexes in human telomeres and gene promoters are being extensively studied for their role in controlling the growth of cancer cells. G-quadruplexes have been unambiguously shown to exist both in vitro and in vivo, including in the guanine (G)-rich DNA genes encoding pre-ribosomal RNA (pre-rRNA), which is transcribed in the cell's nucleolus. Recent studies strongly suggest that these DNA sequences ("rDNA"), and the transcribed rRNA, are a potential anticancer target through the inhibition of RNA polymerase I (Pol I) in ribosome biogenesis, but the structures of ribosomal G-quadruplexes at atomic resolution are unknown and very little biophysical characterization has been performed on them to date. In the present study, circular dichroism (CD) spectroscopy is used to show that two putative rDNA G-quadruplex sequences, NUC 19P and NUC 23P and their counterpart rRNAs, predominantly adopt parallel topologies, reminiscent of the analogous telomeric quadruplex structures. Based on this information, we modeled parallel topology atomistic structures of the putative ribosomal G-quadruplexes. We then validated and refined the modeled ribosomal G-quadruplex structures using all-atom molecular dynamics (MD) simulations with the CHARMM36 force field in the presence and absence of stabilizing K⁺. Motivated by preliminary MD simulations of the telomeric parallel G-quadruplex (TEL 24P) in which the K⁺ ion is expelled, we used updated CHARMM36 force field K⁺ parameters that were optimized, targeting the data from quantum mechanical calculations and the polarizable Drude model force field. In subsequent MD simulations with optimized CHARMM36 parameters, the K⁺ ions are predominantly in the G-quadruplex channel and the rDNA G-quadruplexes have more well-defined, predominantly parallel-topology structures as compared to rRNA. In addition, NUC 19P is more structured than NUC 23P, which contains extended loops. Results from this study set the structural foundation for understanding G-quadruplex functions and the design of novel chemotherapeutics against these nucleolar targets and can be readily extended to other DNA and RNA G-quadruplexes.

An atomistic investigation on the interaction of distamycin A and its derivative with the telomeric G-Quadruplex as anticancer agents by molecular dynamics simulation

Human telomerase that activates within cancer cells has a telomeric sequence at the 3' end. Each factor that stabilizes the G-quadruplex in guanine-rich telomeric sequences can inhibit the regular telomerase activity. Therefore, the telomeric G-quadruplex is known as a promising target in cancer treatment. In this work, we studied the binding of positively charged distamycin A and its uncharged derivative to the G-quadruplex in a solution environment by Molecular Dynamics (MD) simulation. The binding mechanism and subtle conformational changes were investigated as a result of the ligand attachment. Moreover, binding free energy and clustering analysis describe the stability and flexibility of G-quadruplexes upon ligand binding. Structural analyses displayed that the favorable binding of both ligands imposes significant stability and rigidity in G-quadruplex conformation compared to free G-quadruplex, especially charged distamycin. Hydration pattern and ion distribution were different for free G-quadruplex and both of the ligand complexes. Energy decomposition reveals the electrostatic effect on the stability of G-quadruplex. The radial distribution function displayed the solvent shell and ion moving away from the groove. The hydrogen bond played an essential role in the binding of both ligands, especially for the charged derivative. van der Waals interaction is the only factor that is more important in binding uncharged distamycin into G-quadruplex than the charged one. The calculated ΔG_bind showed the stability of both ligands within grooves and good agreement with the experimental binding free energy data. Finally, the results suggest that ligand modification improves the binding mode toward stabilizing G-quadruplexes.

Understanding The Role of Reline, a Natural DES, on Temperature-Induced Conformational Changes of C-Kit G-Quadruplex DNA: A Molecular Dynamics Study

The noncanonical guanine-rich DNAs have drawn particular attention to the scientific world due to their controllable diverse and polymorphic structures. Apart from biological and medical significance, G-quadruplex DNAs are widely used in various fields such as nanotechnology, nanomachine, biosensors, and biocatalyst. So far, the applications of the G-quadruplex DNA are mainly limited in the water medium. Recently, a new generation of solvent named deep eutectic solvent (DES) has become very popular and has been widely used as a reaction medium of biocatalytic reactions and long-term storage medium for nucleic acids, even at high temperature. Hence, it is essential to understand the role of DES on temperature-induced conformational changes of a G-quadruplex DNA. In this research work, we have explored the temperature-mediated conformational dynamics of c-kit oncogene promoter G-quadruplex DNA in reline medium in the temperature range of 300-500 K, using a total of 10 μs unbiased all-atom molecular dynamics simulation. Here, from RMSD, RMSF, R_g and principal component analyses, we notice that the c-kit G-quadruplex DNA is stable up to 450 K in reline medium. However, it unfolds in water medium at 450 K. It is found that the hydrogen bonding interactions between c-kit G-quadruplex DNA and reline play a key role in the stabilization of the G-quadruplex DNA even at high temperature. Furthermore, in this work we have observed a very interesting and distinctive phenomenon of the central cation of the G-quadruplex DNA. Its position was seen to fluctuate between the two tetrad cores, that is, the region between tetrad-1 and tetrad-2 and that between tetrad-2 and tetrad-3 and vice versa at 450 and 500 K in reline medium which is absent in water medium at 450 K. Moreover, the rate of its oscillation is increased when temperature is increased.

Charge-dipole interactions in G-quadruplex thrombin-binding aptamer

DNAs form various structures through hydrogen-bonding, base-stacking and electrostatic interactions. Although these noncovalent interactions are known to be cooperative in stabilizing a G-quadruplex (G4) structure of DNA, we find from all-atom molecular dynamics simulations that the electrostatic charge-dipole interaction is competitive with both hydrogen-bonding and base-stacking interactions. For the thrombin-binding aptamer (TBA) forming a chair-type antiparallel G4 structure, we have examined effects of an intercalating metal ion [K+, Sr2+, Mn+: an ion having a charge of n+ (n = 1-4) with the ionic radius of K+] on structural properties and noncovalent interactions. When K+ in the TBA·K+ complex is replaced with Sr2+, guanine dipoles in the two G-tetrads are realigned toward the central metal ion, thereby distorting the planar G4 geometry. Replacing K+ with Sr2+ significantly enhances the charge-dipole interaction but substantially reduces the number of hydrogen bonds in the G-tetrads. In the case of TBA·Mn+ complexes, as the charge n increases, the charge-dipole interaction increases but both of the hydrogen-bonding and base-stacking interactions decrease. These results suggest that the charge-dipole interaction realigning guanine dipoles in the G-tetrads is not cooperative but competitive with both hydrogen-bonding and base-stacking interactions favoring the planar G-tetrad geometry. Obviously, the charge state of an intercalating metal ion is as important as the ionic radius in forming a stable G4 structure. Thus, a delicate balance between these competing noncovalent interactions makes the chair-type antiparallel G4 structure of TBA selective for intercalating metal ions.

Probing the Limits of Supramolecular G-Quadruplexes Using Atomistic Molecular Dynamics Simulations

Guanosine and related derivatives self-assemble in the presence of cations like potassium into supramolecular G-quadruplexes (SGQs), where four guanine moieties form planar tetrads (T) that coaxially stack into columnar aggregates with broad size distributions. However, SGQs made from 8-aryl-2'-deoxyguanosine derivatives (8ArGs), form mostly octamers, or two-tetrad (2T)-SGQs, while some form dodecamers (3T-SGQs), or hexadecamers (4T-SGQs), and none reported to date form higher assemblies. A theoretical model that addresses the configurational space available for the multiple pathways available for 8ArGs to self-assemble into SGQs is used to frame a series of molecular dynamics simulations (MDS) with selected SGQs. Some key insights from this work include: (a) The predicted entropic costs are not significantly higher for SGQs with more subunits due to their hierarchical assembly pathways; (b) The multiple isomeric SGQs vary in the interfacial contacts between consecutive tetrads, due to their two distinct sides (head, h; tail, t), with the MDS supporting the predicted order of stability of hh > ht > tt for octamers. (c) Such order also applies to dodecamers and hexadecamers, but with context-dependent exceptions due to strong allosteric effects. (d) The main factor disfavoring the tt interface is the repulsive dipolar interactions between the O4' from ribose moieties on adjacent tetrads. (e) SGQs with 5 or more tetrads are disfavored because the attractive interactions are not large or strong enough to overcome the many repulsive forces resulting from the addition of further tetrads. We expect these findings provide some guidelines to enable the further development of SGQs into functional materials.

Origin of Ion Specificity of Telomeric DNA G-Quadruplexes Investigated by Free-Energy Simulations

Telomeric DNA consists of tandem repeats of the sequence d(TTAGGG) that form G-quadruplex structures made of stacked guanines with monovalent cations bound at a central cavity. Although different ions can stabilize a G-quadruplex structure, the preferred bound ions are typically K⁺ or Na⁺. Several different strand-folding topologies have been reported for Q-quadruplexes formed from telomeric repeats depending on DNA length and ion solution condition. This suggests a possible dependence of the ion selectivity of the central pore on the folding topology of the quadruplex. Molecular dynamics free energy perturbation has been employed to systematically study the relative affinity of the central quadruplex pore for different cation types and the associated energetic and solvation contributions to ion selectivity. The calculations have been performed on two different common quadruplex folding topologies. For both topologies, the same ion selectivity was found with a preference for K⁺ followed by Rb⁺ and Na⁺, which agrees with the experimentally determined preference for most investigated quadruplexes. The selectivity is determined by a balance between attractive Coulomb interactions and loss of hydration but also modulated by van der Waals contributions. Specificity is mediated by the central guanines and no significant contribution of the nucleic acid backbone. The simulations indicate that different topologies might be stabilized by ions bound at the surface or alternative sites of the quadruplex because the ion specificity of the central pore does not depend on the strand folding topology.

Exploring ion induced folding of a single-stranded DNA oligomer from molecular simulation studies

One crucial issue in DNA hydration is the effect of salts on its conformational features. This has relevance in biology as cations present in the cellular environment shield the negative charges on the DNA backbone, thereby reducing the repulsive force between them. By screening the negative charges along the backbone, cations stabilize the folded structure of DNA. To study the effect of the added salt on single-stranded DNA (ss-DNA) conformations, we have performed room temperature molecular dynamics simulations of an aqueous solution containing the ss-DNA dodecamer with the 5'-CGCGAATTCGCG-3' sequence in the presence of 0.2, 0.5, and 0.8 M NaCl. Our calculations reveal that in the presence of the salt, the DNA molecule forms more collapsed coil-like conformations due to the screening of negative charges along the backbone. Additionally, we demonstrated that the formation of an octahedral inner-sphere complex by the strongly bound ion plays an important role in the stabilization of such folded conformation of DNA. Importantly, it is found that ion-DNA interactions can also explain the formation of non-sequential base stackings with longer lifetimes. Such non-sequential base stackings further stabilize the collapsed coil-like folded form of the DNA oligomer.

Structural characterization of a dimer of RNA duplexes composed of 8-bromoguanosine modified CGG trinucleotide repeats: a novel architecture of RNA quadruplexes

Fragile X syndrome and fragile X-associated tremor/ataxia syndrome (FXTAS) are neurodegenerative disorders caused by the pathogenic expansion of CGG triplet repeats in the FMR1 gene. FXTAS is likely to be caused by a 'toxic' gain-of-function of the FMR1 mRNA. We provide evidence for the existence of a novel quadruplex architecture comprising CGG repeats. The 8-bromoguanosine ((Br)G)-modified molecule GC(Br)GGCGGC forms a duplex in solution and self-associates via the major groove to form a four-stranded, antiparallel (GC(Br)GGCGGC)4 RNA quadruplex with (Br)G3:G6:(Br)G3:G6 tetrads sandwiched between mixed G:C:G:C tetrads. Self-association of Watson-Crick duplexes to form a four-stranded structure has previously been predicted; however, no experimental evidence was provided. This novel four-stranded RNA structure was characterized using a variety of experimental methods, such as native gel electrophoresis, NMR spectroscopy, small-angle X-ray scattering and electrospray ionization mass spectrometry.

Sodium and Potassium Interactions with Nucleic Acids

Metal ions are essential cofactors for the structure and functions of nucleic acids. Yet, the early discovery in the 70s of the crucial role of Mg(2+) in stabilizing tRNA structures has occulted for a long time the importance of monovalent cations. Renewed interest in these ions was brought in the late 90s by the discovery of specific potassium metal ions in the core of a group I intron. Their importance in nucleic acid folding and catalytic activity is now well established. However, detection of K(+) and Na(+) ions is notoriously problematic and the question about their specificity is recurrent. Here we review the different methods that can be used to detect K(+) and Na(+) ions in nucleic acid structures such as X-ray crystallography, nuclear magnetic resonance or molecular dynamics simulations. We also discuss specific versus non-specific binding to different structures through various examples.

Structural Dynamics of Thrombin-Binding DNA Aptamer d(GGTTGGTGTGGTTGG) Quadruplex DNA Studied by Large-Scale Explicit Solvent Simulations

The thrombin-binding aptamer (15-TBA) is a 15-mer DNA oligonucleotide with sequence d(GGTTGGTGTGGTTGG). 15-TBA folds into a quadruplex DNA (G-DNA) structure with two planar G-quartets connected by three single-stranded loops. The arrangement of the 15-TBA-thrombin complex is unclear, particularly with respect to the precise 15-TBA residues that interact with the thrombin structure. Our present understanding suggests either the 15-TBA single stranded loops containing sequential thymidines (TT) or alternatively a single-stranded loop, containing a guanine flanked by 2 thymidines (TGT), physically associates with thrombin protein. In the present study, the explicit solvent molecular dynamics (MD) simulation method was utilized to further analyze the 15-TBA-thrombin three-dimensional structure. Functional annotation of the loop residues was made with long simulations in the parmbsc0 force field. In total, the elapsed time of simulations carried out in this study exceeds 12 microseconds, substantially surpassing previous G-DNA simulation reports. Our simulations suggest that the TGT-loop function is to stabilize the structure of the aptamer, while the TT-loops participate in direct binding to thrombin. The findings of the present report advance our understanding of the molecular structure of the 15-TBA-thrombin structure further enabling the construction of biosensors for aptamer bases and the development of anticoagulant agents.

Extended molecular dynamics of a c-kit promoter quadruplex

The 22-mer c-kit promoter sequence folds into a parallel-stranded quadruplex with a unique structure, which has been elucidated by crystallographic and NMR methods and shows a high degree of structural conservation. We have carried out a series of extended (up to 10 μs long, ∼50 μs in total) molecular dynamics simulations to explore conformational stability and loop dynamics of this quadruplex. Unfolding no-salt simulations are consistent with a multi-pathway model of quadruplex folding and identify the single-nucleotide propeller loops as the most fragile part of the quadruplex. Thus, formation of propeller loops represents a peculiar atomistic aspect of quadruplex folding. Unbiased simulations reveal μs-scale transitions in the loops, which emphasizes the need for extended simulations in studies of quadruplex loops. We identify ion binding in the loops which may contribute to quadruplex stability. The long lateral-propeller loop is internally very stable but extensively fluctuates as a rigid entity. It creates a size-adaptable cleft between the loop and the stem, which can facilitate ligand binding. The stability gain by forming the internal network of GA base pairs and stacks of this loop may be dictating which of the many possible quadruplex topologies is observed in the ground state by this promoter quadruplex.

Role of the central cations in the mechanical unfolding of DNA and RNA G-quadruplexes

Cations are known to mediate diverse interactions in nucleic acids duplexes but they are critical in the arrangement of four-stranded structures. Here, we use all-atom molecular dynamics simulations with explicit solvent to analyse the mechanical unfolding of representative intramolecular G-quadruplex structures: a parallel, a hybrid and an antiparallel DNA and a parallel RNA, in the presence of stabilising cations. We confirm the stability of these conformations in the presence of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$\rm {K}^+$\end{document} central ions and observe distortions from the tetrad topology in their absence. Force-induced unfolding dynamics is then investigated. We show that the unfolding events in the force-extension curves are concomitant to the loss of coordination between the central ions and the guanines of the G-quadruplex. We found lower ruptures forces for the parallel configuration with respect to the antiparallel one, while the behaviour of the force pattern of the parallel RNA appears similar to the parallel DNA. We anticipate that our results will be essential to interpret the fine structure rupture profiles in stretching assays at high resolution and will shed light on the mechanochemical activity of G-quadruplex-binding machinery.

Molecular dynamics and principal components of potassium binding with human telomeric intra-molecular G-quadruplex

Telomere assumes intra-molecular G-quadruplex that is a significant drug target for inhibiting telomerase maintenance of telomeres in cancer. Metal cations have been recognized as playing important roles in stabilizing G-quadruplex, but their binding processes to human telomeric G-quadruplex remain uncharacterized. To investigate the detailed binding procedures, molecular dynamics simulations were conducted on the hybrid [3 + 1] form-one human telomeric intra-molecular G-quadruplex. We show here that the binding of a potassium ion to a G-tetrad core is mediated by two alternative pathways. Principal component analysis illustrated the dominant concerted motions of G-quadruplex occurred at the loop domains. MM-PBSA calculations revealed that binding was energetically favorable and driven by the electrostatic interactions. The lower binding site was found more constructive favorable for binding. Our data provide useful information on a potassium-mediated stable structure of human telomeric intra-molecular G-quadruplex, implicating in ion disorder associated conformational changes and targeted drug design.

Structural dynamics of possible late-stage intermediates in folding of quadruplex DNA studied by molecular simulations

Explicit solvent molecular dynamics simulations have been used to complement preceding experimental and computational studies of folding of guanine quadruplexes (G-DNA). We initiate early stages of unfolding of several G-DNAs by simulating them under no-salt conditions and then try to fold them back using standard excess salt simulations. There is a significant difference between G-DNAs with all-anti parallel stranded stems and those with stems containing mixtures of syn and anti guanosines. The most natural rearrangement for all-anti stems is a vertical mutual slippage of the strands. This leads to stems with reduced numbers of tetrads during unfolding and a reduction of strand slippage during refolding. The presence of syn nucleotides prevents mutual strand slippage; therefore, the antiparallel and hybrid quadruplexes initiate unfolding via separation of the individual strands. The simulations confirm the capability of G-DNA molecules to adopt numerous stable locally and globally misfolded structures. The key point for a proper individual folding attempt appears to be correct prior distribution of syn and anti nucleotides in all four G-strands. The results suggest that at the level of individual molecules, G-DNA folding is an extremely multi-pathway process that is slowed by numerous misfolding arrangements stabilized on highly variable timescales.

New insights from molecular dynamic simulation studies of the multiple binding modes of a ligand with G-quadruplex DNA

G-quadruplexes are higher-order DNA and RNA structures formed from guanine-rich sequences. These structures have recently emerged as a new class of potential molecular targets for anticancer drugs. An understanding of the three-dimensional interactions between small molecular ligands and their G-quadruplex targets in solution is crucial for rational drug design and the effective optimization of G-quadruplex ligands. Thus far, rational ligand design has been focused mainly on the G-quartet platform. It should be noted that small molecules can also bind to loop nucleotides, as observed in crystallography studies. Hence, it would be interesting to elucidate the mechanism underlying how ligands in distinct binding modes influence the flexibility of G-quadruplex. In the present study, based on a crystal structure analysis, the models of a tetra-substituted naphthalene diimide ligand bound to a telomeric G-quadruplex with different modes were built and simulated with a molecular dynamics simulation method. Based on a series of computational analyses, the structures, dynamics, and interactions of ligand-quadruplex complexes were studied. Our results suggest that the binding of the ligand to the loop is viable in aqueous solutions but dependent on the particular arrangement of the loop. The binding of the ligand to the loop enhances the flexibility of the G-quadruplex, while the binding of the ligand simultaneously to both the quartet and the loop diminishes its flexibility. These results add to our understanding of the effect of a ligand with different binding modes on G-quadruplex flexibility. Such an understanding will aid in the rational design of more selective and effective G-quadruplex binding ligands.

Molecular dynamics simulations to provide new insights into the asymmetrical ammonium ion movement inside of the [d(G3T4G4)]2 G-quadruplex DNA structure

We have used both adaptive biasing force (ABF) and regular molecular dynamics (MD) simulations to investigate the asymmetrical NH(4)(+) ion movement inside of a bimolecular G-quadruplex DNA structure [d(G(3)T(4)G(4))](2). The free-energy landscapes obtained from ABF MD simulations suggest that the NH(4)(+) ion exiting the [d(G(3)T(4)G(4))](2) G-quadruplex stem in the direction toward the edge-type loop (denoted as the upper direction) experiences a lower free-energy barrier than that toward the diagonal loop (denoted as the lower direction) by approximately 3-4 kcal mol(-1). This result is in qualitative agreement with the previous discovery made by Šket and Plavec on the same G-quadruplex structure from (15)N NMR experiments (J. Am. Chem. Soc. 2007, 129, 8794). In the Na(+) form of the same G-quadruplex, Na(+) ion movement was found to be symmetrical, with a free-energy barrier of only 5-7 kcal mol(-1) to cross all three G-quartets, that is, [d(G(3)T(4)G(4))](2) still exhibits ion-channel-like behaviors for Na(+) ions. On the basis of the new computational results, we hypothesize that the stiffness of a G-quartet is primarily determined by the base stacking interactions within the G-quadruplex stem. Therefore, the structural origin for the asymmetrical NH(4)(+) ion movement in [d(G(3)T(4)G(4))](2) is the presence of two different modes of base stacking around the NH(4)(+) binding sites, a more stable 5'-syn-anti mode between lower and central G-quartets and a less stable 5'-anti-anti mode between upper and central G-quartets. Simulations also suggest that loop topology at the end of a G-quadruplex stem only controls the direction at which an exiting NH(4)(+) ion reaches bulk solution but does not impose significant free-energy barriers.

Molecular dynamics simulations of G-DNA and perspectives on the simulation of nucleic acid structures

The article reviews the application of biomolecular simulation methods to understand the structure, dynamics and interactions of nucleic acids with a focus on explicit solvent molecular dynamics simulations of guanine quadruplex (G-DNA and G-RNA) molecules. While primarily dealing with these exciting and highly relevant four-stranded systems, where recent and past simulations have provided several interesting results and novel insight into G-DNA structure, the review provides some general perspectives on the applicability of the simulation techniques to nucleic acids.

Cation binding to 15-TBA quadruplex DNA is a multiple-pathway cation-dependent process

A combination of explicit solvent molecular dynamics simulation (30 simulations reaching 4 µs in total), hybrid quantum mechanics/molecular mechanics approach and isothermal titration calorimetry was used to investigate the atomistic picture of ion binding to 15-mer thrombin-binding quadruplex DNA (G-DNA) aptamer. Binding of ions to G-DNA is complex multiple pathway process, which is strongly affected by the type of the cation. The individual ion-binding events are substantially modulated by the connecting loops of the aptamer, which play several roles. They stabilize the molecule during time periods when the bound ions are not present, they modulate the route of the ion into the stem and they also stabilize the internal ions by closing the gates through which the ions enter the quadruplex. Using our extensive simulations, we for the first time observed full spontaneous exchange of internal cation between quadruplex molecule and bulk solvent at atomistic resolution. The simulation suggests that expulsion of the internally bound ion is correlated with initial binding of the incoming ion. The incoming ion then readily replaces the bound ion while minimizing any destabilization of the solute molecule during the exchange.

Explaining the varied glycosidic conformational, G-tract length and sequence preferences for anti-parallel G-quadruplexes

Guanine-rich DNA sequences tend to form four-stranded G-quadruplex structures. Characteristic glycosidic conformational patterns along the G-strands, such as the 5'-syn-anti-syn-anti pattern observed with the Oxytricha nova telomeric G-quadruplexes, have been well documented. However, an explanation for these featured glycosidic patterns has not emerged. This work presents MD simulation and free energetic analyses for simplified two-quartet [d(GG)](4) models and suggests that the four base pair step patterns show quite different relative stabilities: syn-anti > anti-anti > anti-syn > syn-syn. This suggests the following rule: when folding, anti-parallel G-quadruplexes tend to maximize the number of syn-anti steps and avoid the unfavorable anti-syn and syn-syn steps. This rule is consistent with most of the anti-parallel G-quadruplex structures in the Protein Databank (PDB). Structural polymorphisms of G-quadruplexes relate to these glycosidic conformational patterns and the lengths of the G-tracts. The folding topologies of G2- and G4-tracts are not very polymorphic because each strand tends to populate the stable syn-anti repeat. G3-tracts, on the other hand, cannot present this repeating pattern on each G-tract. This leads to smaller energy differences between different geometries and helps explain the extreme structural polymorphism of the human telomeric G-quadruplexes.

Enantioselective recognition between polydiacetylene nucleolipid monolayers and complementary oligonucleotides

A two-dimensional bio/synthetic hybrid system at the air-solution interface made of a polymerized diacetylene Langmuir film with nucleobase modified headgroups is presented. The polymerized film presents a crystalline array of nucleobases, capable of specific binding of complementary mononucleoside or oligonucleotide sequences. Mixed monolayers of the linear polyconjugated polydiacetylene (PDA) films derivatized with cytosine (10,12-pentacosadiyne-cytidyl, PDC) monomers and alcohol-terminated diacetylene lipid (10,12-pentacosadiynol, PDOH) at a 3:1 ratio (PDC 75%) were compressed and polymerized at the air-water interface with circular polarized light (CPL) or nonpolarized UV light. Here we report a grazing incidence X-ray diffraction (GIXD) investigation of PDC films polymerized to different chirality and hybridized with complementary ssDNA strands. We have demonstrated enantioselective interactions on synthetic structured interfaces produced by Langmuir surface compression followed by polymerization with circular polarized UV light (CPL). The left- and right-CPL polymerized light exhibit the same well-defined crystalline structure. The observed difference between left- and right-CPL polymerized PDC 75% Langmuir films compressed over the complementary mononucleotide guanosine or hybridized with fully complementary ssG(12)T(5) oligonucleotide in the subphase suggests that they are indeed enantiomeric structures, capable of enantioselective binding of their natural ligand, guanosine, solely as a result of surface induced asymmetry in "left" but not in "right" form. This observation may also be related to the intriguing question of chiral selection during the early period of "Origin of Life". We show that achiral compounds, as a result of irradiation with circular polarized light, can organize in chiral surface structures capable of amplification of biopolymer binding of particular handedness.

Molecular modeling and simulation of G-quadruplexes and quadruplex-ligand complexes

Methods for the molecular modeling and simulation of G-quadruplex structures and their drug/ligand complexes are discussed, and a range of protocols is presented for undertaking a variety of tasks including model-building, ligand docking, dynamics simulation, continuum solvent modeling, energetic calculations, principal component analysis, and quantum chemical computations. The scope and limitations of these approaches are discussed.

The impact of monovalent ion force field model in nucleic acids simulations

Different classical models for monovalent ions (typically used to neutralize proteins or nucleic acids) are available in the literature and are widely used in molecular dynamics simulations without a great knowledge of their quality, consistency with the macromolecular force field and impact on the global simulation results. In this paper the ability of several of the most popular ion models to reproduce both quantum mechanics and experimental results is examined. Artefacts due to the use of incorrect ion models in molecular dynamics simulations of concentrated solutions of NaCl and KCl in water and of a short DNA duplex in 500 mM aqueous solutions of NaCl and KCl have been analyzed. Our results allow us to discuss the robustness and reliability of different ion models and to highlight the source of potential errors arising from non-optimal models. However, it is also found that the structural and dynamic characteristics of DNA (as an example of a heavily charged macromolecule) in near-physiological conditions are quite independent of the ion model used, providing support to most already-published simulations of macromolecules.

Determination of alkali and halide monovalent ion parameters for use in explicitly solvated biomolecular simulations

Alkali (Li(+), Na(+), K(+), Rb(+), and Cs(+)) and halide (F(-), Cl(-), Br(-), and I(-)) ions play an important role in many biological phenomena, roles that range from stabilization of biomolecular structure, to influence on biomolecular dynamics, to key physiological influence on homeostasis and signaling. To properly model ionic interaction and stability in atomistic simulations of biomolecular structure, dynamics, folding, catalysis, and function, an accurate model or representation of the monovalent ions is critically necessary. A good model needs to simultaneously reproduce many properties of ions, including their structure, dynamics, solvation, and moreover both the interactions of these ions with each other in the crystal and in solution and the interactions of ions with other molecules. At present, the best force fields for biomolecules employ a simple additive, nonpolarizable, and pairwise potential for atomic interaction. In this work, we describe our efforts to build better models of the monovalent ions within the pairwise Coulombic and 6-12 Lennard-Jones framework, where the models are tuned to balance crystal and solution properties in Ewald simulations with specific choices of well-known water models. Although it has been clearly demonstrated that truly accurate treatments of ions will require inclusion of nonadditivity and polarizability (particularly with the anions) and ultimately even a quantum mechanical treatment, our goal was to simply push the limits of the additive treatments to see if a balanced model could be created. The applied methodology is general and can be extended to other ions and to polarizable force-field models. Our starting point centered on observations from long simulations of biomolecules in salt solution with the AMBER force fields where salt crystals formed well below their solubility limit. The likely cause of the artifact in the AMBER parameters relates to the naive mixing of the Smith and Dang chloride parameters with AMBER-adapted Aqvist cation parameters. To provide a more appropriate balance, we reoptimized the parameters of the Lennard-Jones potential for the ions and specific choices of water models. To validate and optimize the parameters, we calculated hydration free energies of the solvated ions and also lattice energies (LE) and lattice constants (LC) of alkali halide salt crystals. This is the first effort that systematically scans across the Lennard-Jones space (well depth and radius) while balancing ion properties like LE and LC across all pair combinations of the alkali ions and halide ions. The optimization across the entire monovalent series avoids systematic deviations. The ion parameters developed, optimized, and characterized were targeted for use with some of the most commonly used rigid and nonpolarizable water models, specifically TIP3P, TIP4P EW, and SPC/E. In addition to well reproducing the solution and crystal properties, the new ion parameters well reproduce binding energies of the ions to water and the radii of the first hydration shells.

Geometrical and electronic structure variability of the sugar-phosphate backbone in nucleic acids

The anionic sugar-phosphate backbone of nucleic acids substantially contributes to their structural flexibility. To model nucleic acid structure and dynamics correctly, the potentially sampled substates of the sugar-phosphate backbone must be properly described. However, because of the complexity of the electronic distribution in the nucleic acid backbone, its representation by classical force fields is very challenging. In this work, the three-dimensional potential energy surfaces with two independent variables corresponding to rotations around the alpha and gamma backbone torsions are studied by means of high-level ab initio methods (B3LYP/6-31+G*, MP2/6-31+G*, and MP2 complete basis set limit levels). The ability of the AMBER ff99 [Wang, J. M.; Cieplak, P.; Kollman, P. A. J. Comput. Chem. 2000, 21, 1049-1074] and parmbsc0 [Perez, A.; Marchan, I.; Svozil, D.; Sponer, J.; Cheatham, T. E.; Laughten, C. A.; Orozco, M. Biophys. J. 2007, 92, 3817-3829] force fields to describe the various alpha/gamma conformations of the DNA backbone accurately is assessed by comparing the results with those of ab initio quantum chemical calculations. Two model systems differing in structural complexity were used to describe the alpha/gamma energetics. The simpler one, SPM, consisting of a sugar and methyl group linked through a phosphodiester bond was used to determine higher-order correlation effects covered by the CCSD(T) method. The second, more complex model system, SPSOM, includes two deoxyribose residues (without the bases) connected via a phosphodiester bond. It has been shown by means of a natural bond orbital analysis that the SPSOM model provides a more realistic representation of the hyperconjugation network along the C5'-O5'-P-O3'-C3' linkage. However, we have also shown that quantum mechanical investigations of this model system are nontrivial because of the complexity of the SPSOM conformational space. A comparison of the ab initio data with the ff99 potential energy surface clearly reveals an incorrect ff99 force-field description in the regions where the gamma torsion is in the trans conformation. An explanation is proposed for why the alpha/gamma flips are eliminated so successfully when the parmbsc0 force-field modification is used.

Computer simulation study of probe-target hybridization in model DNA microarrays: effect of probe surface density and target concentration

We use lattice Monte Carlo simulations to study the thermodynamics of hybridization of single-stranded "target" genes in solution with complementary "probe" DNA molecules immobilized on a microarray surface. The target molecules in our system contain 48 segments and the probes tethered on a hard surface contain 8-24 segments. The segments on the probe and target are distinct, with each segment representing a sequence of nucleotides that interacts exclusively with its unique complementary target segment with a single hybridization energy; all other interactions are zero. We examine how surface density (number of probes per unit surface area) and concentration of target molecules affect the extent of hybridization. For short probe lengths, as the surface density increases, the probability of binding long stretches of target segments increases at low surface density, reaches a maximum at an intermediate surface density, and then decreases at high surface density. Furthermore, as the surface density increases, the target is less likely to bind completely to one probe; instead, it binds simultaneously to multiple probes. At short probe lengths, as the target concentration increases, the fraction of targets binding completely to the probes (specificity) decreases. At long probe lengths, varying the target concentration does not affect the specificity. At all target concentrations as the probe length increases, the fraction of target molecules bound to the probes by at least one segment (sensitivity) increases while the fraction of target molecules completely bound to the probes (specificity) decreases. This work provides general guidelines to maximizing microarray sensitivity and specificity. Our results suggest that the sensitivity and specificity can be maximized by using probes 130-180 nucleotides long at a surface density in the range of 7 x 10(-5)- 3 x 10(-4) probe molecules per nm(2).

Molecular dynamics simulations and their application to four-stranded DNA

This review provides a critical assessment of the advantages and limitations of modeling methods available for guanine quadruplex (G-DNA) molecules. We characterize the relations of simulations to the experimental techniques and explain the actual meaning and significance of the results. The following aspects are discussed: pair-additive approximation of the empirical force fields, sampling limitations stemming from the simulation time and accuracy of description of base stacking, H-bonding, sugar-phosphate backbone and ions by force fields. Several methodological approaches complementing the classical explicit solvent molecular dynamics simulations are commented on, including enhanced sampling methods, continuum solvent methods, free energy calculations and gas phase simulations. The successes and pitfalls of recent simulation studies of G-DNA are demonstrated on selected results, including studies of cation interactions and dynamics of G-DNA stems, studies of base substitutions (inosine, thioguanine and mixed tetrads), analysis of possible kinetic intermediates in folding pathway of a G-DNA stem and analysis of loop regions of G-DNA molecules.

Stability and migration of metal ions in G4-wires by molecular dynamics simulations

We present a molecular dynamics investigation of guanine quadruple helices based on classical force fields. We analyze the dependence of the helical conformation on various compositional factors, such as the length of the G4-wire, as well as the incorporation into the helix channel of alkali ions of different species and in different amounts. In compliance with previous indications, our results suggest that monovalent alkali cations assist the stability of the quadruplex arrangement against disruption on the few nanoseconds time scale in the order of increasing van der Waals radius. Whereas very short G4-wire fragments immediately unfold in the absence of coordinating metal ions or in the presence of tiny ions (e.g., Li+) in agreement with the experimental evidence that empty short guanine quadruplexes are not formed in any synthetic conditions, our simulations show that longer empty helices do not discompose. This finding supports the possibility of producing long G4-wires with different guanine-cation stoichiometries than those routinely known. The classical trajectories allow us to identify different stationary axial sites for the different metal species, which are confirmed by complementary quantum calculations.

Density functional study of isoguanine tetrad and pentad sandwich complexes with alkali metal ions

Isoguanine tetraplexes and pentaplexes contain two or more stacked polyads with intercalating metal ions. We report here the results of a density functional study of sandwiched isoguanine tetrad and pentad complexes consisting of two polyads with Na(+), K(+) and Rb(+) ions at the B3LYP level. In comparison to single polyad metal ion complexes, there is a trend towards increased non-planarity of the polyads in the sandwich complexes. In general, the pentad sandwiches have relatively planar polyad structures, whereas the tetrad complexes contain highly non-planar polyad building blocks. As in other sandwich complexes and in metal ion complexes with single polyads, the metal ion-base interaction energy plays an essential role. In iG sandwich structures, this interaction energy is slightly larger than in the corresponding guanine sandwich complexes. Because the base-base interaction energy is even more increased in passing from guanine to isoguanine, the isoguanine sandwiches are thus far the only examples where the base-base interaction energy is larger than the base-metal ion interaction energy. Stacking interactions have been studied in smaller models consisting of two bases, retaining the geometry from the complete complex structures. From the data obtained at the B3LYP and BH&H levels and with Møller-Plesset perturbation theory, one can conclude that the B3LYP method overestimates the repulsion in stacked base dimers. For the complexes studied in this work, this is only of minor importance because the direct inter-tetrad or inter-pentad interaction is supplemented by a strong metal ion-base interaction. Using a microsolvation model, the metal ion preference K(+) approximately Rb(+) > Na(+) is found for tetrad complexes. On the other hand, for pentads the ordering is Rb(+) > K(+) > Na(+). In the latter case experimental data are available that agree with this prediction.

Quadruplex DNA: sequence, topology and structure

G-quadruplexes are higher-order DNA and RNA structures formed from G-rich sequences that are built around tetrads of hydrogen-bonded guanine bases. Potential quadruplex sequences have been identified in G-rich eukaryotic telomeres, and more recently in non-telomeric genomic DNA, e.g. in nuclease-hypersensitive promoter regions. The natural role and biological validation of these structures is starting to be explored, and there is particular interest in them as targets for therapeutic intervention. This survey focuses on the folding and structural features on quadruplexes formed from telomeric and non-telomeric DNA sequences, and examines fundamental aspects of topology and the emerging relationships with sequence. Emphasis is placed on information from the high-resolution methods of X-ray crystallography and NMR, and their scope and current limitations are discussed. Such information, together with biological insights, will be important for the discovery of drugs targeting quadruplexes from particular genes.

Computer simulation study of molecular recognition in model DNA microarrays

DNA microarrays have been widely adopted by the scientific community for a variety of applications. To improve the performance of microarrays there is a need for a fundamental understanding of the interplay between the various factors that affect microarray sensitivity and specificity. We use lattice Monte Carlo simulations to study the thermodynamics and kinetics of hybridization of single-stranded target genes in solution with complementary probe DNA molecules immobilized on a microarray surface. The target molecules in our system contain 48 segments and the probes tethered on a hard surface contain 8-24 segments. The segments on the probe and target are distinct and each segment represents a sequence of nucleotides ( approximately 11 nucleotides). Each probe segment interacts exclusively with its unique complementary target segment with a single hybridization energy; all other interactions are zero. We examine how the probe length, temperature, or hybridization energy, and the stretch along the target that the probe segments complement, affect the extent of hybridization. For systems containing single probe and single target molecules, we observe that as the probe length increases, the probability of binding all probe segments to the target decreases, implying that the specificity decreases. We observe that probes 12-16 segments ( approximately 132-176 nucleotides) long gave the highest specificity and sensitivity. This agrees with the experimental results obtained by another research group, who found an optimal probe length of 150 nucleotides. As the hybridization energy increases, the longer probes are able to bind all their segments to the target, thus improving their specificity. The hybridization kinetics reveals that the segments at the ends of the probe are most likely to start the hybridization. The segments toward the center of the probe remain bound to the target for a longer time than the segments at the ends of the probe.

Predictive modelling of topology and loop variations in dimeric DNA quadruplex structures

We have used a combination of simulated annealing (SA), molecular dynamics (MD) and locally enhanced sampling (LES) methods in order to predict the favourable topologies and loop conformations of dimeric DNA quadruplexes with T2 or T3 loops. This follows on from our previous MD simulation studies on the influence of loop lengths on the topology of intramolecular quadruplex structures [P. Hazel et al. (2004) J. Am. Chem. Soc., 126, 16 405–16 415], which provided results consistent with biophysical data. The recent crystal structures of d(G4T3G4)2 and d(G4BrUT2G4) (P. Hazel et al. (2006) J. Am. Chem. Soc., in press) and the NMR-determined topology of d(TG4T2G4T)2 [A.T. Phan et al. (2004) J. Mol. Biol., 338, 93–102] have been used in the present study for comparison with simulation results. These together with MM-PBSA free-energy calculations indicate that lateral T3 loops are favoured over diagonal loops, in accordance with the experimental structures; however, distinct loop conformations have been predicted to be favoured compared to those found experimentally. Several lateral and diagonal loop conformations have been found to be similar in energy. The simulations suggest an explanation for the distinct patterns of observed dimer topology for sequences with T3 and T2 loops, which depend on the loop lengths, rather than only on G-quartet stability.

G-Quadruplexes Can Maintain Their Structure in the Gas Phase

Several very extended (0.5-1 micros) molecular dynamics (MD) simulations of parallel and antiparallel G-quadruplex DNA strongly suggest that in the presence of suitable cations the quadruplex not only remains stable in the gas phase, but also displays a structure that closely resembles that found in extended (25-ns long) trajectories in aqueous solution. In the absence of the crucial cations, the trajectories become unstable and in general the quadruplex structure is lost. To our knowledge, this is the first physiologically relevant structure of DNA for which very large MD simulations suggest that the structure in water and in the gas phase are indistinguishable.

Trisubstituted Acridines as G-quadruplex Telomere Targeting Agents. Effects of Extensions of the 3,6- and 9-Side Chains on Quadruplex Binding, Telomerase Activity, and Cell Proliferation

The synthesis is reported of a group of 3,6,9-trisubstituted acridine compounds as telomeric quadruplex-stabilizing ligands with systematic variations at the 3-, 6-, and 9-positions. A new microwave-assisted methodology has been developed for trisubstituted acridine synthesis. Structure-activity relationships are reported using surface plasmon resonance and a fluorescence melting assay to examine quadruplex binding, together with a telomerase inhibition assay. These reveal relationships between G-quadruplex stabilization and telomerase inhibition and optimal 3,6- and 9-substituent side-chain lengths for maximal activity. Qualitative molecular modeling using molecular dynamics simulations has been undertaken on four quadruplex-DNA complexes. Long-term exposure of MCF7 cancer cells to a subset of the most active compounds, at doses lower than the IC(50) values, showed that one compound produced a marked decrease in population growth, accompanied by senescence, which is consistent with telomere targeting by this agent.

Structural dynamics of precursor and product of the RNA enzyme from the hepatitis delta virus as revealed by molecular dynamics simulations

The hepatitis delta virus (HDV) ribozyme is a self-cleaving RNA enzyme involved in the replication of a human pathogen, the hepatitis delta virus. Recent crystal structures of the precursor and product of self-cleavage, together with detailed kinetic analyses, have led to hypotheses on the catalytic strategies employed by the HDV ribozyme. We report molecular dynamics (MD) simulations (approximately 120 ns total simulation time) to test the plausibility that specific conformational rearrangements are involved in catalysis. Site-specific self-cleavage requires cytidine in position 75 (C75). A precursor simulation with unprotonated C75 reveals a rather weak dynamic binding of C75 in the catalytic pocket with spontaneous, transient formation of a H-bond between U-1(O2') and C75(N3). This H-bond would be required for C75 to act as the general base. Upon protonation in the precursor, C75H+ has a tendency to move towards its product location and establish a firm H-bonding network within the catalytic pocket. However, a C75H+(N3)-G1(O5') H-bond, which would be expected if C75 acted as a general acid catalyst, is not observed on the present simulation timescale. The adjacent loop L3 is relatively dynamic and may serve as a flexible structural element, possibly gated by the closing U20.G25 base-pair, to facilitate a conformational switch induced by a protonated C75H+. L3 also controls the electrostatic environment of the catalytic core, which in turn may modulate C75 base strength and metal ion binding. We find that a distant RNA tertiary interaction involving a protonated cytidine (C41) becomes unstable when left unprotonated, leading to disruptive conformational rearrangements adjacent to the catalytic core. A Na ion temporarily compensates for the loss of the protonated hydrogen bond, which is strikingly consistent with the experimentally observed synergy between low pH and high Na+ concentrations in mediating residual self-cleavage of the HDV ribozyme in the absence of divalents.

Recognition of RNA by amide modified backbone nucleic acids: molecular dynamics simulations of DNA-RNA hybrids in aqueous solution

Thermodynamic and structural properties of a chemically modified DNA-RNA hybrid in which a phosphodiester linkage is replaced by a neutral amide-3 linkage (3'-CH(2)-CONH-5') were investigated using UV melting experiments, molecular dynamics simulations in explicit water, and continuum solvent models. van't Hoff analysis of the experimental UV melting curves suggests that the significant increase of the thermodynamic stability of a 15-mer DNA-RNA with seven alternated amide-3 modifications (+11 degrees C) is mainly due to an increased binding enthalpy. To further evaluate the origin in the observed affinities differences, the electrostatic contribution to the binding free energy was calculated by solving the Poisson-Boltzmann equation numerically. The nonelectrostatic contribution was estimated as the product of a hydrophobic surface tension coefficient and the surface area that is buried upon double strand formation. Structures were taken from 10 ns molecular dynamics simulations computed in a consistent fashion using explicit solvent, counterions, and the particle-mesh Ewald procedure. The present preliminary thermodynamic study suggests that the favorable binding free energy of the amide-3 DNA single strand to the complementary RNA is equally driven by electrostatic and nonpolar contributions to the binding compared to their natural analogues. In addition, molecular dynamics simulations in explicit water were performed on an amide-3 DNA single strand and the corresponding natural DNA. Results from the conformations cluster analysis of the simulated amide-3 DNA single strand ensembles suggest that the 25% of the population sampled within 10 ns has a pre-organized conformation where the sugar C3' endo pucker is favored at the 3'-flanking nucleotides. These structural and thermodynamic features contribute to the understanding of the observed increased affinities of the amide-3 DNA-RNA hybrids at the microscopic level.

Identification of mixed di-cation forms of G-quadruplex in solution

Multinuclear NMR study has demonstrated that G-quadruplex adopted by d(G₃T₄G₄) exhibits two cation binding sites between three of its G-quartets. Titration of tighter binding K⁺ ions into the solution of d(G₃T₄G₄)₂ folded in the presence of NH4+15 ions uncovered a mixed mono-K⁺-mono-NH4+15 form that represents intermediate in the conversion of di-NH4+15 into di-K⁺ form. Analogously, NH4+15 ions were found to replace Na⁺ ions inside d(G₃T₄G₄)₂ quadruplex. The preference of NH4+15 over Na⁺ ions for the two binding sites is considerably smaller than the preference of K⁺ over NH4+15 ions. The two cation binding sites within the G-quadruplex core differ to such a degree that NH4+15 ions bound to the site, which is closer to the edge-type loop, are always replaced first during titration by K⁺ ions. The second binding site is not taken up by K⁺ ion until K⁺ ion already resides at the first binding site. Quantitative analysis of concentrations of the three di-cation forms, which are in slow exchange on the NMR time scale, at 12 K⁺ ion concentrations afforded equilibrium binding constants. K⁺ ion binding to sites U and L within d(G₃T₄G₄)₂ is more favorable with respect to NH4+15 ions by Gibbs free energies of approximately −24 and −18 kJ mol⁻¹ which includes differences in cation dehydration energies, respectively.