Atomistic Picture for the Folding Pathway of a Hybrid-1 Type Human Telomeric DNA G-quadruplex

In this work we studied the folding process of the hybrid-1 type human telomeric DNA G-quadruplex with solvent and ions explicitly modeled. Enabled by the powerful bias-exchange metadynamics and large-scale conventional molecular dynamic simulations, the free energy landscape of this G-DNA was obtained for the first time and four folding intermediates were identified, including a triplex and a basically formed quadruplex. The simulations also provided atomistic pictures for the structures and cation binding patterns of the intermediates. The results showed that the structure formation and cation binding are cooperative and mutually supporting each other. The syn/anti reorientation dynamics of the intermediates was also investigated. It was found that the nucleotides usually take correct syn/anti configurations when they form native and stable hydrogen bonds with the others, while fluctuating between two configurations when they do not. Misfolded intermediates with wrong syn/anti configurations were observed in the early intermediates but not in the later ones. Based on the simulations, we also discussed the roles of the non-native interactions. Besides, the formation process of the parallel conformation in the first two G-repeats and the associated reversal loop were studied. Based on the above results, we proposed a folding pathway for the hybrid-1 type G-quadruplex with atomistic details, which is new and more complete compared with previous ones. The knowledge gained for this type of G-DNA may provide a general insight for the folding of the other G-quadruplexes.     

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.     

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.     

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₇) 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.     

Formation pathways of a guanine-quadruplex DNA revealed by molecular dynamics and thermodynamic analysis of the substates

The formation of a cation-stabilized guanine quadruplex (G-DNA) stem is an exceptionally slow process involving complex kinetics that has not yet been characterized at atomic resolution. Here, we investigate the formation of a parallel stranded G-DNA stem consisting of four strands of d(GGGG) using molecular dynamics simulations with explicit inclusion of counterions and solvent. Due to the limitations imposed by the nanosecond timescale of the simulations, rather than watching for the spontaneous formation of G-DNA, our approach probes the stability of possible supramolecular intermediates (including two-, three-, and four-stranded assemblies with out-of-register base pairing between guanines) on the formation pathway. The simulations suggest that "cross-like" two-stranded assemblies may serve as nucleation centers in the initial formation of parallel stranded G-DNA quadruplexes, proceeding through a series of rearrangements involving trapping of cations, association of additional strands, and progressive slippage of strands toward the full stem. To supplement the analysis, approximate free energies of the models are obtained with explicit consideration of the integral cations. The approach applied here serves as a prototype for qualitatively investigating other G-DNA molecules using molecular dynamics simulation and free-energy analysis.     

The ionic atmosphere around A-RNA: Poisson-Boltzmann and molecular dynamics simulations

The distributions of different cations around A-RNA are computed by Poisson-Boltzmann (PB) equation and replica exchange molecular dynamics (MD). Both the nonlinear PB and size-modified PB theories are considered. The number of ions bound to A-RNA, which can be measured experimentally, is well reproduced in all methods. On the other hand, the radial ion distribution profiles show differences between MD and PB. We showed that PB results are sensitive to ion size and functional form of the solvent dielectric region but not the solvent dielectric boundary definition. Size-modified PB agrees with replica exchange molecular dynamics much better than nonlinear PB when the ion sizes are chosen from atomistic simulations. The distribution of ions 14 Å away from the RNA central axis are reasonably well reproduced by size-modified PB for all ion types with a uniform solvent dielectric model and a sharp dielectric boundary between solvent and RNA. However, this model does not agree with MD for shorter distances from the A-RNA. A distance-dependent solvent dielectric function proposed by another research group improves the agreement for sodium and strontium ions, even for shorter distances from the A-RNA. However, Mg²⁺ distributions are still at significant variances for shorter distances.     

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.     

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

Multinuclear NMR study has demonstrated that G-quadruplex adopted by d(G3T4G4) exhibits two cation binding sites between three of its G-quartets. Titration of tighter binding K+ ions into the solution of d(G3T4G4)2 folded in the presence of 15NH4+ ions uncovered a mixed mono-K+-mono-15NH4+ form that represents intermediate in the conversion of di-15NH4+ into di-K+ form. Analogously, 15NH4+ ions were found to replace Na+ ions inside d(G3T4G4)2 quadruplex. The preference of 15NH4+ over Na+ ions for the two binding sites is considerably smaller than the preference of K+ over 15NH4+ ions. The two cation binding sites within the G-quadruplex core differ to such a degree that 15NH4+ 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(G3T4G4)2 is more favorable with respect to 15NH4+ ions by Gibbs free energies of approximately -24 and -18 kJ mol(-1) which includes differences in cation dehydration energies, respectively.     

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.     

NMR structure of the thrombin-binding DNA aptamer stabilized by Sr2+

The structure of thrombin-binding DNA aptamer complexed with a single Sr2+ ion (Sr2+:TBA complex) has been determined using NMR spectroscopy and restrained molecular dynamics simulations. The quadruplex structure for the Sr2+:TBA complex is similar in topology, but distinct in structure, from that previously reported for the K+:TBA complex. The inter-tetrad distance of the Sr2+:TBA complex is 3.8 angstroms, or 0.7 angstroms larger than in the K+:TBA complex. This substantial difference can be attributed to a different binding site for Sr2+ in the Sr2+:TBA complex than for K+ in the K+:TBA complex. The Sr2+:TBA complex assumes a 1:1 stoichiometry, and it is very likely that the Sr2+ ion simultaneously interacts with the eight O6 atoms of the two G-tetrads. The results indicate that quadruplex DNA structures are highly sensitive to the presence of specific metal ions. The binding of specific metal ions may modulate the biological activity of quadruplex DNA structures in vivo.     

Protein grabs a ligand by extending anchor residues: molecular simulation for Ca2+ binding to calmodulin loop

The structural difference in proteins between unbound and bound forms directly suggests the importance of the conformational plasticity of proteins. However, pathways that connect two-end structures and how they are coupled to the binding reaction are not well understood at atomic resolution. Here, we analyzed the free-energy landscape, explicitly taking into account coupling between binding and conformational change by performing atomistic molecular dynamics simulations for Ca2+ binding to a calmodulin loop. Using the AMBER force field with explicit water solvent, we conducted umbrella sampling for the free-energy surface and steered molecular dynamics for the pathway search. We found that, at an early stage of binding, some key residue side chains extend their "arms" to catch Ca2+ and, after catching, they carry the Ca2+ to the center of the binding pocket. This grabbing motion resulted in smooth and stepwise exchange in coordination partners of Ca2+ from water oxygen to atoms in the calmodulin loop. The key residue that first caught the ion was one of the two acidic residues, which are highly conserved. In the pathway simulations, different pathways were observed between binding and dissociation reactions: The former was more diverse than the latter.     

The NorM MATE Transporter from N. gonorrhoeae: Insights into Drug and Ion Binding from Atomistic Molecular Dynamics Simulations

The multidrug and toxic compound extrusion transporters extrude a wide variety of substrates out of both mammalian and bacterial cells via the electrochemical gradient of protons and cations across the membrane. The substrates transported by these proteins include toxic metabolites and antimicrobial drugs. These proteins contribute to multidrug resistance in both mammalian and bacterial cells and are therefore extremely important from a biomedical perspective. Although specific residues of the protein are known to be responsible for the extrusion of solutes, mechanistic details and indeed structures of all the conformational states remain elusive. Here, we report the first, to our knowledge, simulation study of the recently resolved x-ray structure of the multidrug and toxic compound extrusion transporter, NorM from Neisseria gonorrhoeae (NorM_NG). Multiple, atomistic simulations of the unbound and bound forms of NorM in a phospholipid lipid bilayer allow us to identify the nature of the drug-protein/ion-protein interactions, and secondly determine how these interactions contribute to the conformational rearrangements of the protein. In particular, we identify the molecular rearrangements that occur to enable the Na⁺ ion to enter the cation-binding cavity even in the presence of a bound drug molecule. These include side chain flipping of a key residue, GLU-261 from pointing toward the central cavity to pointing toward the cation binding side when bound to a Na⁺ ion. Our simulations also provide support for cation binding in the drug-bound and apo states of NorM_NG.     

The NorM MATE Transporter from N.?gonorrhoeae: Insights into Drug and Ion Binding from Atomistic Molecular Dynamics Simulations

The multidrug and toxic compound extrusion transporters extrude a wide variety of substrates out of both mammalian and bacterial cells via the electrochemical gradient of protons and cations across the membrane. The substrates transported by these proteins include toxic metabolites and antimicrobial drugs. These proteins contribute to multidrug resistance in both mammalian and bacterial cells and are therefore extremely important from a biomedical perspective. Although specific residues of the protein are known to be responsible for the extrusion of solutes, mechanistic details and indeed structures of all the conformational states remain elusive. Here, we report the first, to our knowledge, simulation study of the recently resolved x-ray structure of the multidrug and toxic compound extrusion transporter, NorM from Neisseria gonorrhoeae (NorM_NG). Multiple, atomistic simulations of the unbound and bound forms of NorM in a phospholipid lipid bilayer allow us to identify the nature of the drug-protein/ion-protein interactions, and secondly determine how these interactions contribute to the conformational rearrangements of the protein. In particular, we identify the molecular rearrangements that occur to enable the Na⁺ ion to enter the cation-binding cavity even in the presence of a bound drug molecule. These include side chain flipping of a key residue, GLU-261 from pointing toward the central cavity to pointing toward the cation binding side when bound to a Na⁺ ion. Our simulations also provide support for cation binding in the drug-bound and apo states of NorM_NG.     

Molecular dynamics simulation of the induced‐fit binding process of DNA aptamer and L‐argininamide

Aptamers are rare functional nucleic acids with binding affinity to and specificity for target ligands. Recent experiments have lead to the proposal of an induced‐fit binding mechanism for L ‐argininamide (Arm) and its binding aptamer. However, at the molecular level, this mechanism between the aptamer and its coupled ligand is still poorly understood. The present study used explicit solvent molecular dynamics (MD) simulations to examine the critical bases involved in aptamer‐Arm binding and the induced‐fit binding process at atomic resolution. The simulation results revealed that the Watson‐Crick pair (G10‐C16), C9, A12, and C17 bases play important roles in aptamer‐Arm binding, and that binding of Arm results in an aptamer conformation optimized through a general induced‐fit process. In an aqueous solution, the mechanism has the following characteristic stages: (a) adsorption stage, the Arm anchors to the binding site of aptamer with strong electrostatic interaction; (b) binding stage, the Arm fits into the binding site of aptamer by hydrogen‐bond formation; and (c) complex stabilization stage, the hydrogen bonding and electrostatic interactions cooperatively stabilize the complex structure. This study provides dynamics information on the aptamer‐ligand induced‐fit binding mechanism. The critical bases in aptamer‐ligand binding may provide a guideline in aptamer design for molecular recognition engineering.     

Determination of the number and location of the manganese binding sites of DNA quadruplexes in solution by EPR and NMR.

The paramagnetic metal ion Mn2+ has been used to probe the electrostatic potentials of a DNA quadruplex that has two quartets with an overall fold of the chair type. A quadruplex with a basket type structure has also been examined. The binding of the paramagnetic ion manganese to these quadruplex DNAs has been investigated by solution state electron paramagnetic resonance (EPR) and nuclear magnetic resonance (NMR) spectroscopies. The EPR results indicate that the DNA aptamer, d(GGTTGGTGTGGTTGG), binds two manganese ions and that the binding constants for each of these sites is approximately 10(5) M-1. The NMR results indicate that the binding sites of the manganese are in the narrow grooves of this quadruplex DNA. The binding sites of the DNA quadruplex formed by dimers of d(GGGGTTTTGGGG) which forms a basket structure are also in the narrow groove. These results indicate that the close approach of phosphates in the narrow minor grooves of the quadruplex structures provide strong binding sites for the manganese ions and that EPR and NMR monitoring of manganese binding can be used to distinguish between the different types of quadruplex structures.     

Molecular dynamics simulations of RNA: an in silico single molecule approach

RNA molecules are now known to be involved in the processing of genetic information at all levels, taking on a wide variety of central roles in the cell. Understanding how RNA molecules carry out their biological functions will require an understanding of structure and dynamics at the atomistic level, which can be significantly improved by combining computational simulation with experiment. This review provides a critical survey of the state of molecular dynamics (MD) simulations of RNA, including a discussion of important current limitations of the technique and examples of its successful application. Several types of simulations are discussed in detail, including those of structured RNA molecules and their interactions with the surrounding solvent and ions, catalytic RNAs, and RNA-small molecule and RNA-protein complexes. Increased cooperation between theorists and experimentalists will allow expanded judicious use of MD simulations to complement conceptually related single molecule experiments. Such cooperation will open the door to a fundamental understanding of the structure-function relationships in diverse and complex RNA molecules. .     

Dynamics, hydration, and motional averaging of a loop-gated artificial protein cavity: the W191G mutant of cytochrome c peroxidase in water as revealed by molecular dynamics simulations

Five molecular dynamics simulations of the W191G cavity mutant of cytochrome c peroxidase in explicit water reveal distinct dynamic and hydration behavior depending on the closed or open state of the flexible loop gating the cavity, the binding of (K+ or small molecule) cations, and the system temperature. The conformational spaces sampled by the loop region and by the cavity significantly reduce upon binding. The largest ordering factor on water dynamics is the presence of the K+ ion occupying the gated cavity. Considerable water exchange occurs for the open-gate cavity when no ligand or cation is bound. In all cases, good correspondence is found between the calculated (ensemble-averaged) location of water molecules and the water sites determined by X-ray crystallography experiments. However, our simulations suggest that these sites do not necessarily correspond to the presence of bound water molecules. In fact, individual water molecules may repeatedly exchange within the cavity volume yet occupy on average these water sites. Four major conclusions emerge. First, it seems misleading to interpret the conformation of protein loop regions in terms of single dominant structures. Second, our simulations support the general picture of Pro 190 cis-trans isomerization as a determinant of the loop-opening mechanism. Third, receptor flexibility is fundamental for ligand binding and molecular recognition, and our results suggest its importance for the docking of small compounds to the artificial cavity. Fourth, after validation against the available experimental data, molecular dynamics simulations can be used to characterize the dynamics and exchange of water molecules and ions, providing atomic level and time-dependent information otherwise inaccessible to experiments.     

Molecular dynamics of DNA quadruplex molecules containing inosine, 6-thioguanine and 6-thiopurine

The ability of the four-stranded guanine (G)-DNA motif to incorporate nonstandard guanine analogue bases 6-oxopurine (inosine, I), 6-thioguanine (tG), and 6-thiopurine (tI) has been investigated using large-scale molecular dynamics simulations. The simulations suggest that a G-DNA stem can incorporate inosines without any marked effect on its structure and dynamics. The all-inosine quadruplex stem d(IIII)(4) shows identical dynamical properties as d(GGGG)(4) on the nanosecond time scale, with both molecular assemblies being stabilized by monovalent cations residing in the channel of the stem. However, simulations carried out in the absence of these cations show dramatic differences in the behavior of d(GGGG)(4) and d(IIII)(4). Whereas vacant d(GGGG)(4) shows large fluctuations but does not disintegrate, vacant d(IIII)(4) is completely disrupted within the first nanosecond. This is a consequence of the lack of the H-bonds involving the N2 amino group that is not present in inosine. This indicates that formation of the inosine quadruplex could involve entirely different intermediate structures than formation of the guanosine quadruplex, and early association of cations in this process appears to be inevitable. In the simulations, the incorporation of 6-thioguanine and 6-thiopurine sharply destabilizes four-stranded G-DNA structures, in close agreement with experimental data. The main reason is the size of the thiogroup leading to considerable steric conflicts and expelling the cations out of the channel of the quadruplex stem. The G-DNA stem can accommodate a single thioguanine base with minor perturbations. Incorporation of a thioguanine quartet layer is associated with a large destabilization of the G-DNA stem whereas the all-thioguanine quadruplex immediately collapses.