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

Large-scale molecular dynamics (MD) simulations have been utilized to study G-DNA quadruplex molecules containing mixed GCGC and all-guanine GGGG quartet layers. Incorporation of mixed GCGC quartets into G-DNA stems substantially enhances their sequence variability. The mixed quadruplexes form rigid assemblies that require integral monovalent cations for their stabilization. The interaction of cations with the all-guanine quartets is the leading contribution for the stability of the four-stranded assemblies, while the mixed quartets are rather tolerated within the structure. The simulations predict that two cations are preferred to stabilize a four-layer quadruplex stem composed of two GCGC and two all-guanine quartets. The distribution of cations in the structure is influenced by the position of the GCGC quartets within the quadruplex, the presence and arrangement of thymidine loops connecting the guanine/cytosine stretches forming the stems, and the cation type present (Na(+) or K(+)). The simulations identify multiple nanosecond-scale stable arrangements of the thymidine loops present in the molecules investigated. In these thymidine loops, several structured pockets are identified capable of temporarily coordinating cations. However, no stable association of cations to a loop has been observed. The simulations reveal several paths through the thymidine loop regions that can be followed by the cations when exchanging between the central ion channel in the quadruplex stem and the surrounding solvent. We have carried out 20 independent simulations while the length of simulations reaches a total of 90 ns, rendering this study one of the most extensive MD investigations carried out on nucleic acids so far. The trajectories provide a largely converged characterization of the structural dynamics of these four-stranded G-DNA molecules.     

Fast gas-phase hydrogen/deuterium exchange observed for a DNA G-quadruplex

The gas-phase hydrogen/deuterium (H/D) exchange kinetics of DNA G-quadruplexes has been investigated using Fourier transform ion cyclotron resonance mass spectrometry (FTICRMS). The quadruplex [(TGGGGT)4 . 3NH4+] undergoes very fast H/D exchange, in both the positive and in the negative ion modes, compared to DNA duplexes and other quadruplexes tested, and compared to the corresponding single-stranded TGGGGT. Substitution of NH4+ for K+ did not alter this fast H/D exchange, indicating that the hydrogens of the ammonium ions are not those exchanged. However, stripping of the interior cations of the quadruplex by source collision-induced dissociation (CID) in the positive ion mode showed that the presence of the inner cations is essential for the fast exchange to be possible. Molecular dynamics simulations show that the G-quadruplex is very rigid in the gas phase with NH4+ ions inside the tetrads. We suggest that the fast H/D exchange is favored by this rigid quadruplex conformation. This example illustrates that the concept that compact DNA structures exchange H for D slower than unfolded ones is a misconception.     

Spectroscopic detection of DNA quadruplexes by vibrational circular dichroism

The four-stranded G-quadruplex motif is a conformation frequently adopted by guanine-rich nucleic acids that plays an important role in biology, medicine, and nanotechnology. Although vibrational spectroscopy has been widely used to investigate nucleic acid structure, association of particular spectral features with the quadruplex structure has to date been ambiguous. In this work, experimental IR absorption and vibrational circular dichroism (VCD) spectra of the model quadruplex systems d(G)(8) and deoxyguanosine-5'-monophosphate (5'-dGMP) were analyzed using molecular dynamics (MD) and quantum-chemical modeling. The experimental spectra were unambiguously assigned to the quadruplex DNA arrangement, and several IR and VCD bands related to this structural motif were determined. Involvement of MD in the modeling was essential for realistic simulation of the spectra. The VCD signal was found to be more sensitive to dynamical structural variations than the IR signal. The combination of the spectroscopic techniques with multiscale simulations provides extended information about nucleic acid conformations and their dynamics.     

Nature of minor-groove binders-DNA complexes in the gas phase

The structure of noncovalent complexes of DNA duplex with minor groove binders (mG-binders) has been analyzed by state of the art molecular dynamics (MD) simulations. More than 3.3 micros of MD trajectories (including 4 x 0.5 micros trajectories) were collected for the Dickerson's dodecamer bound to DAPI, Hoechst 33258, and Netropsin. Comparison of these trajectories with control simulations in water allowed us to determine that the extreme dehydration and partial neutralization occurring during electrospray experiments does not produce the disruption of the DNA:mG-binder complexes or the dissociation of the two strands of the duplex. Irrespective of the drug and the simulation conditions the mG-binders remains bound to the DNA near the preferential binding position in aqueous conditions. Large distortions appear in the two DNA strands, which maintain however a memory of the original DNA duplex structure in water, and a general helical-like conformation.     

Gas-phase stability of G-quadruplex DNA determined by electrospray ionization tandem mass spectrometry and molecular dynamics simulations

The relative gas-phase stabilities of seven quadruplex DNA structures, [d(TG(4)T)](4), [d(T(2)G(3)T)](4), [d(G(4)T(4)G(4))](2), [d(T(2)AG(3))(2)](2), d(T(2)AG(3))(4), d(T(2)G(4))(4), and d(G(2)T(4))(4), were investigated using molecular dynamics simulations and electrospray ionization mass spectrometry (ESI-MS). MD simulations revealed that the G-quadruplexes maintained their structures in the gas phase although the G-quartets were distorted to some degree and ammonium ions, retained by [d(TG(4)T)](4) and [d(T(2)G(3)T)](4), played a key role in stabilizing the tetrad structure. Energy-variable collisional activated dissociation was used to assess the relative stabilities of each quadruplex based on E(1/2) values, and the resulting order of relative stabilities was found to be [d(TG(4)T)](4) > d(T(2)AG(3))(4) approximately d(T(2)G(4))(4) > [d(T(2)G(3)T)](4) > [d(T(2)AG(3))(2)](2) approximately d(G(2)T(4))(4) approximately [d(G(4)T(4)G(4))](2.) The stabilities from the E(1/2) values generally paralleled the RMSD and relative free energies of the quadruplexes based on the MD energy analysis. One exception to the general agreement is [d(G(4)T(4)G(4))](2), which had the lowest E(1/2) value, but was determined to be the most stable quadruplex according to the free-energy analysis and ranked fourth based on the RMSD comparison. This discrepancy is attributed to differences in the fragmentation pathway of the quadruplex.     

Car-Parrinello MD simulations for the Na+ -phenylalanine complex in aqueous solution

Car-Parrinello molecular dynamics (CPMD) calculations are presented for a Na (+)(Phe) complex in aqueous solution and for various stable Na (+)(Phe) complexes and Na (+)(H 2O) n clusters in the gas phase (with up to six water molecules). The CPMD results are compared to available experimental and ab initio reference data, to DFT results obtained with various combinations of density functionals and basis sets, and to previous classical mechanics MD simulations. The agreement with the reference data in the gas phase validates the CPMD method, showing that it is a valid approach for studying these systems and that it describes correctly the competing Na (+)-Phe and Na (+)-H 2O interactions. Analysis of MD trajectories reveals that the Na (+)(Phe) complex in aqueous solution maintains a stable configuration in which the Na (+) cation hovers above the phenyl ring, at an average distance of 3.85 A from the ring center, while remaining strongly bound to one of the carboxylic oxygens of Phe. Constrained MD simulations indicate that the free energy barrier opposing dissociation of the complex exceeds 5.5 kcal/mol. We thus confirm that "cation- pi" interactions between alcali cations and the pi ring, combined with other kinds of interactions, may allow aromatic amino acids to overcome the competition with water in binding a cation.     

Implementations of Nosé-Hoover and Nosé-Poincaré thermostats in mesoscopic dynamic simulations with the united-residue model of a polypeptide chain

Molecular dynamics (MD) simulations generate a canonical ensemble only when integration of the equations of motion is coupled to a thermostat. Three extended phase space thermostats, one version of Nose-Hoover and two versions of Nose-Poincare, are compared with each other and with the Berendsen thermostat and Langevin stochastic dynamics. Implementation of extended phase space thermostats was first tested on a model Lennard-Jones fluid system; subsequently, they were implemented with our physics-based protein united-residue (UNRES) force field MD. The thermostats were also implemented and tested for the multiple-time-step reversible reference system propagator (RESPA). The velocity and temperature distributions were analyzed to confirm that the proper canonical distribution is generated by each simulation. The value of the artificial mass constant, Q, of the thermostat has a large influence on the distribution of the temperatures sampled during UNRES simulations (the velocity distributions were affected only slightly). The numerical stabilities of all three algorithms were compared with each other and with that of microcanonical MD. Both Nose-Poincare thermostats, which are symplectic, were not very stable for both the Lennard-Jones fluid and UNRES MD simulations started from nonequilibrated structures which implies major changes of the potential energy throughout a trajectory. Even though the Nose-Hoover thermostat does not have a canonical symplectic structure, it is the most stable algorithm for UNRES MD simulations. For UNRES with RESPA, the "extended system inside-reference system propagator algorithm" of the RESPA implementation of the Nose-Hoover thermostat was the only stable algorithm, and enabled us to increase the integration time step.     

Structure of triplex DNA in the gas phase

Extensive (more than 90 microseconds) molecular dynamics simulations complemented with ion-mobility mass spectrometry experiments have been used to characterize the conformational ensemble of DNA triplexes in the gas phase. Our results suggest that the ensemble of DNA triplex structures in the gas phase is well-defined over the experimental time scale, with the three strands tightly bound, and for the most abundant charge states it samples conformations only slightly more compact than the solution structure. The degree of structural alteration is however very significant, mimicking that found in duplex and much larger than that suggested for G-quadruplexes. Our data strongly supports that the gas phase triplex maintains an excellent memory of the solution structure, well-preserved helicity, and a significant number of native contacts. Once again, a linear, flexible, and charged polymer as DNA surprises us for its ability to retain three-dimensional structure in the absence of solvent. Results argue against the generally assumed roles of the different physical interactions (solvent screening of phosphate repulsion, hydrophobic effect, and solvation of accessible polar groups) in modulating the stability of DNA structures.     

Macroscopic evidences for non-Rice-Ramsperger-Kassel effects in the reaction between H3O+ and D2O: the occurrence of nonstatistical isotopic branching ratio

The dynamics of the isotopic scrambling in the energized and metastable complex D2O-H3O+ has been studied using classical molecular dynamics (MD) trajectories starting from regions of phase space corresponding to an already formed collisional complex. The simulations cover the range of internal energies spanned by gas phase collision experiments. Rate constants for the isotopic exchange and the complex dissociation have been computed; the isotopic branching ratio R=[HD2O+]/[H2DO+] has also been obtained from MD simulations and found to deviate substantially from an equivalent prediction based on a previously proposed kinetic scheme. This finding suggests the possibility that details of the reaction dynamics play a role in defining the isotopic branching ratio. The analysis of trajectory results indicated a relatively long lifetime for the collisional complex and the presence of multiple time scales for the exchange process, with a large fraction of the exchange events being separated only by a single oxygen-oxygen vibration or half of it. The occurrence of these fast consecutive jumps and their different probabilities as a function of the relative direction between first and second jumps suggest the presence of ballistic motion in the complex following each reactive event. This can be explained on the basis of overlapping regions in phase space and it is used to provide an explanation of the difference between kinetic and MD branching ratios.     

Simulations of vibrational spectra from classical trajectories: calibration with ab initio force fields

An algorithm allowing simulating vibrational spectra from classical time-dependent trajectories was applied for infrared absorption, vibrational circular dichroism, Raman, and Raman optical activity of model harmonic systems. The implementation of the theory within the TINKER molecular dynamics (MD) program package was tested with ab initio harmonic force fields in order to determine the feasibility for more extended MD simulations. The results suggest that sufficiently accurate frequencies can be simulated with integration time steps shorter than about 0.5 fs. For a given integration time step, lower vibrational frequencies ( approximately 0-2000 cm(-1)) could be reproduced with a higher accuracy than higher-frequency vibrational modes (e.g., O-H and C-H stretching). In principle, the algorithm also provides correct intensities for ideal systems. In applied simulations, however, the intensity profiles are affected by an unrealistic energy distribution between normal modes and a slow energy relaxation. Additionally, the energy fluctuations may cause weakening of the intensities on average. For ab initio force fields, these obstacles could be overcome by an arbitrary normal mode energy correction. For general MD simulations, averaging of many shorter MD trajectories started with randomly distributed atomic velocities provided the best spectral shapes. alpha-pinene, D-gluconic acid, formaldehyde dimer, and the acetylprolineamide molecule were used in the tests.     

Stability and cations coordination of DNA and RNA 14-mer G-quadruplexes: a multiscale computational approach

Molecular dynamics simulations have been used to study the differences between two DNA and RNA 14-mer quadruplexes of analogous sequences. Their structures present a completely different fold: DNA forms a bimolecular quadruplex containing antiparallel strands and diagonal loops; RNA forms an intrastrand parallel quadruplex containing a G-tetrad and an hexad, which dimerizes by hexad stacking. We used a multiscale computational approach combining classical Molecular dynamics simulations and density functional theory calculations to elucidate the difference in stability of the 2-folds and their ability in coordinating cations. The presence of 2'-OH groups in the RNA promotes the formation of a large number of intramolecular hydrogen bonds that account for the difference in fold and stability of the two 14-mers. We observe that the adenines in the RNA quadruplex play a key role in conserving the geometry of the hexad. We predict the cation coordination mode of the two quadruplexes, not yet observed experimentally, and we offer a rationale for the corresponding binding energies involved.     

Formation of Human Telomeric G-quadruplex Structures Induced by the Quaternary Benzophenanthridine Alkaloids:Sanguinarine,Nitidine,and Chelerythrine

The ligands which can facilitate the formation and stabilize G‐quadruplex structures have attracted enormous attention due to their potential ability of inhibiting the telomerase activity and halting tumor cell proliferation. It is noteworthy that the abilities of the quaternary benzophenanthridine alkaloids (QBAs), the very important G‐quadruplex binders, in inducing the formation of human telomeric DNA G‐quadruplex structures, have not been reported. Herein, the interaction between single‐strand human telomeric DNA and three QBAs: Sanguinarine (San), Nitidine (Nit) and Chelerythrine (Che), has been investigated. Although these molecules are very similar in structure, they exhibit significantly different abilities in inducing oligonucleotide d(TTAGGG)₄ (HT4) to specific G‐quadruplex structures. Our experimental results indicated that the best ligand San could convert HT4 into antiparallel G‐quadruplex structure completely, followed by Nit, which could transform to mixed‐type or hybrid G‐quadruplex structure partially, whereas Che could only transform to antiparallel G‐quadruplex structure in small quantities. The relative QBAs' inducing abilities as indicated by the CD data are in the order of San>Nit>Che. Further investigation revealed that the G‐quadruplex structures from HT4 induced by QBAs are of intramolecular motif. And only sequences with certain length could be induced by QBAs because of their positive charges which could not attract short chain DNA molecules to close to each other and form intermolecular G‐quadruplex. In addition, the factors that affect the interaction between HT4 and QBAs were discussed. It is proposed that the thickness of the molecular frame and the steric hindrance are the primary reasons why the subtle differences in QBAs' structure lead to their remarkable differences in inducing the formation of the G‐quadruplex structures.     

Theoretical study of a new DNA structure: the antiparallel Hoogsteen duplex

The structure of a new form of duplex DNA, the antiparallel Hoogsteen duplex, is studied in polyd(AT) sequences by means of state-of-the-art molecular dynamics simulations in aqueous solution. The structure, which was found to be stable in all of the simulations, has many similarities with the standard Watson-Crick duplex in terms of general structure, flexibility, and molecular recognition patterns. Accurate MM-PB/SA (and MM-GB/SA) analysis shows that the new structure has an effective energy similar to that of the B-type duplex, while it is slightly disfavored by intramolecular entropic considerations. Overall, MD simulations strongly suggest that the antiparallel Hoogsteen duplex is an accessible structure for a polyd(AT) sequence, which might compete under proper experimental conditions with normal B-DNA. MD simulations also suggest that chimeras containing Watson-Crick duplex and Hoogsteen antiparallel helices might coexist in a common structure, but with the differential characteristics of both type of structures preserved.     

A molecular dynamics study of the pentacyclo-undecane cage amino acid tripeptide

In this paper, we report on the conformational profile of the pentacyclo-undecane (PCU) cage tripeptide carried out by molecular dynamics (MD) simulation using water as an explicit solvent. The MD solution phase studies carried on the model peptide analogues (A)=Ac–Ala–Ala–Ala–NHMe; (B)=Ac–Cage–Cage–Cage–NHMe; (C)=Ac–Ala–Cage–Ala–NHMe and (D)=Ac–Ala–Pro–Ala–NHMe, are used as a complimentary technique to the corresponding gas phase simulated annealing (SA) study previously carried out in our laboratory. No significant structural changes were observed over the MD trajectories. However, the results reported here provide further evidence that the (PCU) cage amino acid exhibits C_7eq, C_7aq, _R and _L conformations, and the theoretical results suggest that the PCU cage amino acid is a strong β-turn inducer. These results support the prediction that when the PCU cage residues are in the (i) and (i+2) positions, the β-turn can be extended in either direction to form anti-parallel β-pleated sheets, thereby forming the basis of the mechanism for the folding back of the chain in a cross-β-turn structure.     

Rectification of the current in alpha-hemolysin pore depends on the cation type: the alkali series probed by MD simulations and experiments

A striking feature of the alpha-hemolysin channel-a prime candidate for biotechnological applications-is the dependence of its ionic conductance on the magnitude and direction of the applied bias. Through a combination of lipid bilayer single-channel recording and molecular dynamics (MD) simulations, we characterized the current-voltage relationship of alpha-hemolysin for all alkali chloride salts at neutral pH. The rectification of the ionic current was found to depend on the type of cations and increase from Li(+) to Cs(+). Analysis of the MD trajectories yielded a simple quantitative model that related the ionic current to the electrostatic potential, the concentration and effective mobility of ions in the channel. MD simulations reveal that the major contribution to the current asymmetry and rectification properties originates from the cationic contribution to the current that is significantly reduced in a cationic dependent way when the membrane polarity is reversed. The variation of chloride current was found to be less important. We report that the differential affinity of cations for the charged residues positioned at the channel's end modulates the number of ions inside the channel stem thus affecting the current properties. Through direct comparison of simulation and experiment, this study evaluates the accuracy of the MD method for prediction of the asymmetric, voltage dependent conductances of a membrane channel.     

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.