Density functional study of guanine and uracil quartets and of guanine quartet/metal ion complexes

Interaction of 5-fluorouracil anticancer drug with nucleobases: insight from DFT, TD-DFT, and AIM calculations

In this study, the interaction of 5-fluorouracil (5FU) drug with adenine (A), guanine(G), cytosine(C), uracil (U), and thymine (T) nucleobases of DNA and RNA are surveyed at the ωB97XD/LANL2DZ, M06-2X/6-31G (d, p), MPW1PWQ1/6-31G(d, p), PBEPBE/6-31(d, p) and ωB97XD/6-31G(d, p) levels of density functional theory (DFT). The considered complexes of 5FU drug with nucleobases are optimized at the above level of theories. Max Force and RMS of optimization criteria are 0.00035 (Ha), and 0.0003 respectively. From optimized structures, the adsorption energy, thermodynamic parameters in gas and solvent media, quantum theory atom in molecule (QTAIM), electron localized function (ELF), and reduced density gradient (RDG) are calculated at ωB97XD/LANL2DZ and M06-2X/6-31G (d, p) level of DFT theory. The QTAIM, ELF, and RDG results confirm that the nature of bonding between 5FU drug with A, C, G, U, and T nucleobases is electrostatic or hydrogen bond type. The adsorption and thermodynamic energy results demonstrate that the interaction of the 5FU drug with C and G nucleobases is stronger than other nucleobases. The results of this study can be suggested the mechanism of interaction of the 5FU drug with nucleobases of DNA and RNA.Communicated by Ramaswamy H. Sarma.

Ion-Dependent Conformational Plasticity of Telomeric G-Hairpins and G-Quadruplexes

Telomeric DNA is guanine-rich and can adopt structures such as G-quadruplexes (GQs) and G-hairpins. Telomeric GQs influence genome stability and telomerase activity, making understanding of enzyme-GQ interactions and dynamics important for potential drug design. GQs have a characteristic tetrad core, which is connected by loop regions. Within this architecture are G-hairpins, fold-back motifs that are thought to represent the first intermediate in GQ folding. To better understand the relationship between G-hairpin motifs and GQs, we performed polarizable simulations of a two-tetrad telomeric GQ and an isolated SC11 telomeric G-hairpin. The telomeric GQ contains a G-triad, which functions as part of the tetrad core or linker regions, depending on local conformational change. This triad and another motif below the tetrad core frequently bound ions and may represent druggable sites. Further, we observed the unbiased formation of a G-triad and a G-tetrad in simulations of the SC11 G-hairpin and found that cations can be partially hydrated while facilitating the formation of these motifs. Finally, we demonstrated that K⁺ ions form specific interactions with guanine bases, while Na⁺ ions interact nonspecifically with bases in the structure. Together, these simulations provide new insights into the influence of ions on GQs, G-hairpins, and G-triad motifs.

Proton-Transfer Reactions in One-Electron-Oxidized G-Quadruplexes: A Density Functional Theory Study

Recently, G-quadruplexes (Gq) formed in B-DNA as secondary structures are found to be important therapeutic targets and material for developing nanodevices. Gq are guanine-rich and thus susceptible to oxidative damage by producing short-lived intermediate radicals via proton-transfer reactions. Understanding the mechanisms of radical formation in Gq is of fundamental interest to understand the early stages of DNA damage. Herein, we used density functional theory including aqueous phase (ωB97XD-PCM/6-31++G**) and considered single layer of Gq [G-quartets (G4): association of four guanines in a cyclic Hoogsteen hydrogen-bonded arrangement (Scheme 1)] to unravel the mechanisms of formation of intermediates by calculating the relative Gibbs free energies and spin density distributions of one-electron-oxidized G4 and its various proton-transfer states: G^•+, G(N₁-H)^•, G(N₂-H')^•, G(N₂-H″)^•, G(N₁-H)^•-(H⁺O₆)G, and G(N₂-H)^•-(H⁺N₇)G. The present calculation predicts the formation of G(N₂-H)^•-(H⁺N₇)G, which is only ca. 0.8 kcal/mol higher in energy than the initially formed G^•+. The formation of G(N₂-H)^•-(H⁺N₇)G plays a key role in explaining the formation of 8-OG along with G(N₁-H)^• formation via tautomerization from G(N₂-H)^•, as proposed recently.

Alkaline earth cations binding mode tailors excited-state charge transfer properties of guanine quadruplex: A TDDFT study

Quadruplexes formed by nucleic acids and their derivates tend to chelate different monovalent and bivalent cations, which simultaneously affect their excited electronic states properties. Cation binding to every and every other cavity of the central ion channel could be exploited for tuning exited-state charge transfer properties. In this work we utilize set of descriptors constructed on the basis of the one-electron transition density matrix obtained using linear-response TDDFT to study excited states properties of four crystallized tetramolecular quadruplexes that chelate alkaline earth cations (Ca²⁺, Sr²⁺ and Ba²⁺). Here, we show that alkaline earth cations situated at adjacent vacancies promote existence of the nucleobase-metal charge separation (CS) states, contrary to the structures with cations that occupy every second available vacancy. We argued that stabilization of these CS states is due to the strong electric field that stabilizes d orbitals of the cations which accept an excited-electron. Moreover, CS content is increased and redshifted below the first bright transition when number of the chelated cations is increased. Hydration effects stabilized CS states and increased their relative content. We also identified electron detachment states in the broad energy range for the Ca²⁺ containing system. These findings are valuable for understanding and development of the novel nanostructures based on the quadruplex scaffold with adjustable optical properties.

Cation competition and recruitment around the c-kit1 G-quadruplex using polarizable simulations

Nucleic acid-ion interactions are fundamentally important to the physical, energetic, and conformational properties of DNA and RNA. These interactions help fold and stabilize highly ordered secondary and tertiary structures, such as G-quadruplexes (GQs), which are functionally relevant in telomeres, replication initiation sites, and promoter sequences. The c-kit proto-oncogene encodes for a receptor tyrosine kinase and is linked to gastrointestinal stromal tumors, mast cell disease, and leukemia. This gene contains three unique GQ-forming sequences that have proposed antagonistic effects on gene expression. The dominant GQ, denoted c-kit1, has been shown to decrease expression of c-kit transcripts, making the c-kit1 GQ a promising drug target. Toward disease intervention, more information is needed regarding its conformational dynamics and ion binding properties. Therefore, we performed molecular dynamics simulations of the c-kit1 GQ with K⁺, Na⁺, Li⁺, and mixed salt solutions using the Drude-2017 polarizable force field. We evaluated GQ structure, ion sampling, core energetics, ion dehydration and binding, and ion competition and found that each analysis supported the known GQ-ion specificity trend (K⁺ > Na⁺ > Li⁺). We also found that K⁺ ions coordinate in the tetrad core antiprismatically, whereas Na⁺ and Li⁺ align coplanar to guanine tetrads, partially because of their attraction to surrounding water. Further, we showed that K⁺ occupancy is higher around the c-kit1 GQ and its nucleobases than Na⁺ and Li⁺, which tend to interact with backbone and sugar moieties. Finally, we showed that K⁺ binding to the c-kit1 GQ is faster and more frequent than Na⁺ and Li⁺. Such descriptions of GQ-ion dynamics suggest the rate of dehydration as the dominant factor for preference of K⁺ by DNA GQs and provide insight into noncanonical nucleic acids for which little experimental data exist.

The significance of the metal cation in guanine-quartet – metalloporphyrin complexes

The distinct positions of the divalent metal ions with respect to the porphyrin ring are responsible for different interaction energies between metalloporphyrins and the guanine quartet.

Halogen-Bonded Guanine Base Pairs, Quartets and Ribbons

Halogen bonding is studied in different structures consisting of halogenated guanine DNA bases, including the Hoogsteen guanine–guanine base pair, two different types of guanine ribbons (R-I and R-II) consisting of two or three monomers, and guanine quartets. In the halogenated base pairs (except the Cl-base pair, which has a very non-planar structure with no halogen bonds) and R-I ribbons (except the At trimer), the potential N-X•••O interaction is sacrificed to optimise the N-X•••N halogen bond. In the At trimer, the astatines originally bonded to N1 in the halogen bond donating guanines have moved to the adjacent O6 atom, enabling O-At•••N, N-At•••O, and N-At•••At halogen bonds. The brominated and chlorinated R-II trimers contain two N-X•••N and two N-X•••O halogen bonds, whereas in the iodinated and astatinated trimers, one of the N-X•••N halogen bonds is lost. The corresponding R-II dimers keep the same halogen bond patterns. The G-quartets display a rich diversity of symmetries and halogen bond patterns, including N-X•••N, N-X•••O, N-X•••X, O-X•••X, and O-X•••O halogen bonds (the latter two facilitated by the transfer of halogens from N1 to O6). In general, halogenation decreases the stability of the structures. However, the stability increases with the increasing atomic number of the halogen, and the At-doped R-I trimer and the three most stable At-doped quartets are more stable than their hydrogenated counterparts. Significant deviations from linearity are found for some of the halogen bonds (with halogen bond angles around 150°).

A vibrational spectroscopic and computational study of gaseous protonated and alkali metal cationized G-C base pairs

The structures and properties of metal cationized complexes of 9-ethylguanine (9eG) and 1-methylcytosine (1mC), (9eG:1mC)M⁺, where M⁺ = Li⁺, Na⁺, K⁺, Rb⁺, Cs⁺ as well as the protonated complex, (9eG:1mC)H⁺, have been studied using a combination of IRMPD spectroscopy and computational methods. For (9eG:1mC)H⁺, the dominant structure is a Hoogsteen type complex with the proton covalently bound to N3 of 1mC despite this being the third best protonation site of the two bases; based on proton affinities N7 of 9eG should be protonated. However, this structural oddity can be explained considering both the number of hydrogen bonds that can be formed when N3 of 1mC is protonated as well as the strong ion-induced dipole interaction that exists between an N3 protonated 1mC and 9eG due to the higher polarizability of 9eG. The anomalous dissociation of (9eG:1mC)H⁺, forming much more (1mC)H⁺ than would be predicted based on the computed thermochemistry, can be explained as being due to the structural oddity of the protonation site and that the barrier to proton transfer from N3 of 1mC to N7 of 9eG grows dramatically as the base pair begins to dissociate. For the (9eG:1mC)M⁺; M = Li⁺, Na⁺, K⁺, Rb⁺, Cs⁺ complexes, single unique structures could not be assigned. However, the experimental spectra were consistent with the computed spectra. For (9eG:1mC)Li⁺, the lowest energy structure is one in which Li⁺ is bound to O6 of 9eG and both O2 and N3 of 1mC; there is also an interbase hydrogen bond from the amine of 1mC to N7 of 9eG. For Na⁺, K⁺, and Rb⁺, similar binding of the metal cation to 1mC is calculated but, unlike Li⁺, the lowest energy structure is one in which the metal cation is bound to N7 of 9eG; there is also an interbase hydrogen bond between the amine of 1mC and the carbonyl of 9eG. The lowest energy structure for the Cs complex is the Watson-Crick type base pairing with Cs⁺ binding only to 9eG through O6 and N7 and with three hydrogen bonds between 9eG and 1mC. It also interesting to note that the Watson-Crick base pairing structure gets lower in Gibbs energy relative to the lowest energy complexes as the metal gets larger. This indicates that the smaller, more densely charged cations have a greater propensity to interfere with Watson-Crick base pairing than do the larger, less densely charged metal cations.

Integration of Fuzzy Matter-Element Method and 3D-QSAR Model for Generation of Environmentally Friendly Quinolone Derivatives

The environmental pollution of quinolone antibiotics (QAs) has caused rising public concern due to their widespread usage. In this study, Gaussian 09 software was used to obtain the infrared spectral intensity (IRI) and ultraviolet spectral intensity (UVI) of 24 QAs based on the Density Functional Theory (DFT). Rather than using two single-factor inputs, a fuzzy matter-element method was selected to calculate the combined effects of infrared and ultraviolet spectra (CI). The Comparative Molecular Field Analysis (CoMFA) was then used to construct a three-dimensional quantitative structure–activity relationship (3D-QSAR) with QAs’ molecular structure as the independent variable and CI as the dependent variable. Using marbofloxacin and levofloxacin as target molecules, the molecular design of 87 QA derivatives was carried out. The developed models were further used to determine the stability, functionality (genetic toxicity), and the environmental effects (bioaccumulation, biodegradability) of these designed QA derivatives. Results indicated that all QA derivatives are stable in the environment with their IRI, UVI, and CI enhanced. Meanwhile, the genetic toxicity of the 87 QA derivatives increased by varying degrees (0.24%–29.01%), among which the bioaccumulation and biodegradability of 43 QA derivatives were within the acceptable range. Through integration of fuzzy matter-element method and 3D-QSAR, this study advanced the QAs research with the enhanced CI and helped to generate the proposed environmentally friendly quinolone derivatives so as to aid the management of this class of antibiotics.

Polarizable Molecular Dynamics Simulations of Two c-kit Oncogene Promoter G-Quadruplexes: Effect of Primary and Secondary Structure on Loop and Ion Sampling

G-quadruplexes (GQs) are highly ordered nucleic acid structures that play fundamental roles in regulating gene expression and maintaining genomic stability. GQs are topologically diverse and enriched in promoter sequences of growth regulatory genes and proto-oncogenes, suggesting that they may serve as attractive targets for drug design at the level of transcription rather than inhibiting the activity of the protein products of these genes. The c-kit promoter contains three adjacent GQ-forming sequences that have proposed antagonistic effects on gene expression and thus are promising drug targets for diseases such as gastrointestinal stromal tumors, mast cell disease, and leukemia. Because GQ stability is influenced by primary structure, secondary structure, and ion interactions, a greater understanding of GQ structure, dynamics, and ion binding properties is needed to develop novel, GQ-targeting therapeutics. Here, we performed molecular dynamics simulations to systematically study the c-kit2 and c-kit* GQs, evaluating nonpolarizable and polarizable force fields (FFs) and examining the effects of base substitutions and cation type (K⁺, Na⁺, and Li⁺) on the dynamics of their isolated and linked structures. We found that the Drude polarizable FF outperformed the additive CHARMM36 FF in two- and three-tetrad GQs and solutions of KCl, NaCl, and LiCl. Drude simulations with different cations agreed with the known GQ stabilization preference (K⁺ > Na⁺ > Li⁺) and illustrated that tetrad core-ion coordination differs as a function of cation type. Finally, we showed that differences in primary and secondary structure influence loop sampling, ion binding, and core-ion energetics of GQs.

Self-assembling purine and pteridine quartets: how do π-conjugation patterns affect resonance-assisted hydrogen bonding?

Computed association strengths for 43 purine and pteridine quartets (38 to 100 kcal mol-1) show excellent linear correlation with π-conjugation gain in the assembled monomers (r2 = 0.965). Even quartets having the same secondary electrostatic interactions can display very different association strengths depending on the π-conjugation patterns of the monomeric units.

An IRMPD spectroscopic and computational study of protonated guanine-containing mismatched base pairs in the gas phase

Infrared multiple photon dissociation spectroscopy has been used to probe the structures of the three protonated base-pair mismatches containing 9-ethylguanine (9eG) in the gas phase. Some of these protonated base-pairs have been identified in RNA.

Quantum Mechanical Investigation of the G-Quadruplex Systems of Human Telomere

The three G-quadruplexes involved in the human telomere have been studied with an accurate quantum mechanical approach, and the possibility of reducing them to a simpler model has been tested. The similarities and the differences of these three systems are shown and discussed. Each system has been analyzed through different properties and compared to the others. In particular, we have considered: (1) the shape of the cavity and the atomic charges around it; (2) the electric field in and out of the cavity; (3) the stabilization energy due to the stacking of G-tetrads, to the H-bonds and to the ion interactions; and, finally, (4) to study the mechanism of the process of the ion inclusion in the cavity, the curves of potential energy due to the movement of the Na+ and K+ ions toward the cavity. The results suggest that a detailed study is essential in order to obtain the quantitative properties of these complex systems, but also that some qualitative behaviors can be schematized. Our study makes it clear that the entry of an ion in the cavity of these systems is a complex process, where it is possible to find stable structures with the ion out and in the cavity. Moreover, it is possible that more than one diabatic state is involved in this process.

Umbrella sampling molecular dynamics simulations reveal concerted ion movement through G-quadruplex DNA channels

We have applied the umbrella sampling (US) method in all-atom molecular dynamics (MD) simulations to obtain potential of mean force (PMF) profiles for ion transport through three representative G-quadruplex DNA channels: [d(TG₄T)]₄, [d(G₃T₄G₄)]₂, and d[G₄(T₄G₄)₃]. The US MD results are in excellent agreement with those obtained previously with the adaptive biasing force (ABF) method. We then utilized the unique features in the US MD method to investigate multi-ion effects in [d(G₃T₄G₄)]₂ and discovered that the concerted ion movement is crucial for fully explaining the unusual experimental results on ion movement in this particular G-quadruplex system. We anticipate that these modern free-energy methods will be useful tools in evaluating ion transport properties of other G-quadruplex DNA channels.

Theoretical study of gas and solvent phase stability and molecular adsorption of noncanonical guanine bases on graphene

The gas and solvent phase stability of noncanonical (Gua)_n nucleobases is investigated in the framework of dispersion-corrected density functional theory (DFT). The calculated results strongly support the high tendency for the dimerization of (Gua)_n bases in both gas and solvent phases. An interplay between intermolecular and bifurcated H-bonds is suggested to govern the stability of (Gua)_n bases which bears a correlation with the description of dispersion correction terms employed in the DFT calculations. For example, a higher polarity is predicted for (Gua)_n bases by the dispersion-corrected DFT in contrast to the non-polar nature of (Gua)₃ and (Gua)₄ predicted by the hybrid meta-GGA calculations. This distinct variation becomes significant under physiological conditions as polar (Gua)_n is likely to exhibit greater stabilization in the gas phase compared to solvated (Gua)_n. Graphene acting as a substrate induces modification in base configurations via maximization of π-orbital overlap between the base and substrate. In solvent, the substrate-induced effects are further heightened with lowering of the dipole moments of (Gua)_n as also displayed by the corresponding isosurface of the electrostatic potential. The graphene-induced stability in both gas and solvent phases appears to fulfill one of the prerequisite criteria for molecular self-assembly. The DFT results therefore provide atomistic insights into the stability and molecular assembly of free-standing noncanonical (Gua)_n nucleobases which can be extended to understanding the self-assembly process of functional biomolecules on 2D materials for potential biosensing applications.

Implication of the solvent effect, metal ions and topology in the electronic structure and hydrogen bonding of human telomeric G-quadruplex DNA

We present a first-principles density functional study elucidating the effects of solvent, metal ions and topology on the electronic structure and hydrogen bonding of 12 well-designed three dimensional G-quadruplex (G4-DNA) models in different environments. Our study shows that the parallel strand structures are more stable in dry environments and aqueous solutions containing K(+) ions within the tetrad of guanine but conversely, that the anti-parallel structure is more stable in solutions containing the Na(+) ions within the tetrad of guanine. The presence of metal ions within the tetrad of the guanine channel always enhances the stability of the G4-DNA models. The parallel strand structures have larger HOMO-LUMO gaps than antiparallel structures, which are in the range of 0.98 eV to 3.11 eV. Partial charge calculations show that sugar and alkali ions are positively charged whereas nucleobases, PO4 groups and water molecules are all negatively charged. Partial charges on each functional group with different signs and magnitudes contribute differently to the electrostatic interactions involving G4-DNA and favor the parallel structure. A comparative study between specific pairs of different G4-DNA models shows that the Hoogsteen OH and NH hydrogen bonds in the guanine tetrad are significantly influenced by the presence of metal ions and water molecules, collectively affecting the structure and the stability of G4-DNA.

Temperature and pressure limits of guanosine monophosphate self-assemblies

Guanosine monophosphate, among the nucleotides, has the unique property to self-associate and form nanoscale cylinders consisting of hydrogen-bonded G-quartet disks, which are stacked on top of one another. Such self-assemblies describe not only the basic structural motif of G-quadruplexes formed by, e.g., telomeric DNA sequences, but are also interesting targets for supramolecular chemistry and nanotechnology. The G-quartet stacks serve as an excellent model to understand the fundamentals of their molecular self-association and to unveil their application spectrum. However, the thermodynamic stability of such self-assemblies over an extended temperature and pressure range is largely unexplored. Here, we report a combined FTIR and NMR study on the temperature and pressure stability of G-quartet stacks formed by disodium guanosine 5′-monophosphate (Na₂5′-GMP). We found that under abyssal conditions, where temperatures as low as 5 °C and pressures up to 1 kbar are reached, the self-association of Na₂5′-GMP is most favoured. Beyond those conditions, the G-quartet stacks dissociate laterally into monomer stacks without significantly changing the longitudinal dimension. Among the tested alkali cations, K⁺ is the most efficient one to elevate the temperature as well as the pressure limits of GMP self-assembly.

Effect of Monovalent Ion Parameters on Molecular Dynamics Simulations of G-Quadruplexes

G-quadruplexes (GQs) are key noncanonical DNA and RNA architectures stabilized by desolvated monovalent cations present in their central channels. We analyze extended atomistic molecular dynamics simulations (∼580 μs in total) of GQs with 11 monovalent cation parametrizations, assessing GQ overall structural stability, dynamics of internal cations, and distortions of the G-tetrad geometries. Majority of simulations were executed with the SPC/E water model; however, test simulations with TIP3P and OPC water models are also reported. The identity and parametrization of ions strongly affect behavior of a tetramolecular d[GGG]₄ GQ, which is unstable with several ion parametrizations. The remaining studied RNA and DNA GQs are structurally stable, though the G-tetrad geometries are always deformed by bifurcated H-bonding in a parametrization-specific manner. Thus, basic 10-μs-scale simulations of fully folded GQs can be safely done with a number of cation parametrizations. However, there are parametrization-specific differences and basic force-field errors affecting the quantitative description of ion-tetrad interactions, which may significantly affect studies of the ion-binding processes and description of the GQ folding landscape. Our d[GGG]₄ simulations indirectly suggest that such studies will also be sensitive to the water models. During exchanges with bulk water, the Na⁺ ions move inside the GQs in a concerted manner, while larger relocations of the K⁺ ions are typically separated. We suggest that the Joung-Cheatham SPC/E K⁺ parameters represent a safe choice in simulation studies of GQs, though variation of ion parameters can be used for specific simulation goals.

Cooperative vibrational properties of hydrogen bonds in Watson–Crick DNA base pairs

For the AT pair, Symst and Strech peaks further shift toward the red, giving the H-bonds an amplified effect (orange arrows).

Quantum mechanical investigation of G-quartet systems of DNA

Minima of the electric field and positions of K⁺ and Na⁺ (zero of the x-coordinate is the center of the cavity).

A comparison of four different conformations adopted by human telomeric G-quadruplex using computer simulations

The telomeric G-quadruplexes for their unique structural features are considered as potential anticancer drug targets. These, however, exhibit structural polymorphism as different topology types for the intra-molecular G-quadruplexes from human telomeric G-rich sequences have been reported based on NMR spectroscopy and X-ray crystallography. These techniques provide detailed atomic-level information about the molecule but relative conformational stability of the different topologies remains unsolved. Therefore, to understand the conformational preference, we have carried out quantum chemical calculations on G-quartets; used all-atom molecular dynamics (MD) simulations and steered molecular dynamics (SMD) simulations to characterize the four human telomeric G-quadruplex topologies based on its G-tetrad core-types, viz., parallel, anti-parallel, mixed-(3 + 1)-form1 and mixed-(3 + 1)-form2. We have also studied a non-telomeric sequence along with these telomeric forms giving a comparison between the two G-rich forms. The structural properties such as base pairing, stacking geometry and backbone conformations have been analyzed. The quantum calculations indicate that presence of a sodium ion inside the G-tetrad plane or two potassium ions on both sides of the plane give it an overall planarity which is much needed for good stacking to form a helix. MD simulations indicate that capping of the G-tetrad core by the TTA loops keep the terminal guanine bases away from water. The SMD simulations along with equilibrium MD studies indicate that the parallel and non-telomeric forms are comparatively less stable. We could come to the conclusion that the anti-parallel form and also the mixed-(3 + 1)-form1 topology are most likely to represent the major conformation., 2016. © 2015 Wiley Periodicals, Inc. Biopolymers 105: 83-99, 2016.

MoD-QM/MM Structural Refinement Method: Characterization of Hydrogen Bonding in the Oxytricha nova G-Quadruplex

A generalization of the Moving-Domain Quantum Mechanics/Molecular Mechanics (MoD-QM/MM) hybrid method [Gascon, J. A.; Leung, S. S. F.; Batista, E. R.; Batista, V. S. J. Chem. Theory Comput. 2006, 2, 175-186] is introduced to provide a self-consistent computational protocol for structural refinement of extended systems. The method partitions the system into molecular domains that are iteratively optimized as quantum mechanical (QM) layers embedded in their surrounding molecular environment to obtain an ab initio quality description of the geometry and the molecular electrostatic potential of the extended system composed of those constituent fragments. The resulting methodology is benchmarked as applied to model systems that allow for full QM optimization as well as through refinement of the hydrogen bonding geometry in Oxytricha nova guanine quadruplex for which several studies have been reported, including the X-ray structure and NMR data. Calculations of (1)H NMR chemical shifts based on the gauge independent atomic orbital (GIAO) method and direct comparisons with experiments show that solvated MoD-QM/MM structures, sampled from explicit solvent molecular dynamics simulations, allow for NMR simulations in much improved agreement with experimental data than models based on the X-ray structure or those optimized using classical molecular mechanics force fields.

Spectroscopic evidence for a gas-phase librating G-quartet–Na+ complex

The Na⁺–G₄ adduct exists as a mixture of two metastable populations, rapidly interconverting at room temperature.

Structural dynamics of human telomeric G-quadruplex loops studied by molecular dynamics simulations

Loops which are linkers connecting G-strands and supporting the G-tetrad core in G-quadruplex are important for biological roles of G-quadruplexes. TTA loop is a common sequence which mainly resides in human telomeric DNA (hTel) G-quadruplex. A series of molecular dynamics (MD) simulations were carried out to investigate the structural dynamics of TTA loops. We found that (1) the TA base pair formed in TTA loops are very stable, the occupied of all hydrogen bonds are more than 0.95. (2) The TA base pair makes the adjacent G-quartet more stable than others. (3) For the edgewise loop and the diagonal loop, most loop bases are stacking with others, only few bases have considerable freedom. (4) The stabilities of these stacking structures are distinct. Part of the loops, especially TA base pairs, and bases stacking with the G-quartet, maintain certain stable conformations in the simulation, but other parts, like TT and TA stacking structures, are not stable enough. For the first time, spontaneous conformational switches of TTA edgewise loops were observed in our long time MD simulations. (5) For double chain reversal loop, it is really hard to maintain a stable conformation in the long time simulation under present force fields (parm99 and parmbsc0), as it has multiple conformations with similar free energies.

Relative stability of different DNA guanine quadruplex stem topologies derived using large-scale quantum-chemical computations

We provide theoretical predictions of the intrinsic stability of different arrangements of guanine quadruplex (G-DNA) stems. Most computational studies of nucleic acids have applied Molecular Mechanics (MM) approaches using simple pairwise-additive force fields. The principle limitation of such calculations is the highly approximate nature of the force fields. In this study, we for the first time apply accurate QM computations (DFT-D3 with large atomic orbital basis sets) to essentially complete DNA building blocks, seven different folds of the cation-stabilized two-quartet G-DNA stem, each having more than 250 atoms. The solvent effects are approximated by COSMO continuum solvent. We reveal sizable differences between MM and QM descriptions of relative energies of different G-DNA stems, which apparently reflect approximations of the DNA force field. Using the QM energy data, we propose correction to earlier free energy estimates of relative stabilities of different parallel, hybrid, and antiparallel G-stem folds based on classical simulations. The new energy ranking visibly improves the agreement between theory and experiment. We predict the 5'-anti-anti-3' GpG dinucleotide step to be the most stable one, closely followed by the 5'-syn-anti-3' step. The results are in good agreement with known experimental structures of 2-, 3-, and 4-quartet G-DNA stems. Besides providing specific results for G-DNA, our study highlights basic limitations of force field modeling of nucleic acids. Although QM computations have their own limitations, mainly the lack of conformational sampling and the approximate description of the solvent, they can substantially improve the quality of calculations currently relying exclusively on force fields.

How to understand quantum chemical computations on DNA and RNA systems? A practical guide for non-specialists

In this review primarily written for non-experts we explain basic methodological aspects and interpretation of modern quantum chemical (QM) computations applied to nucleic acids. We introduce current reference QM computations on small model systems consisting of dozens of atoms. Then we comment on recent advance of fast and accurate dispersion-corrected density functional theory methods, which will allow computations of small but complete nucleic acids building blocks in the near future. The qualitative difference between QM and molecular mechanics (MM, force field) computations is discussed. We also explain relation of QM and molecular simulation computations to experiments.

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.

Thymine quintets and their higher order assemblies studied by electrospray ionization mass spectrometry and theoretical calculation

We previously reported that thymine molecules can specifically form a pentameric magic number cluster named as thymine quintet in the presence of K(+) , Rb(+) and Cs(+) . Actually, thymine decamer and doubly charged thymine 15-mer metaclusters can be observed along with thymine quintet in the ESI mass spectra of thymine with the addition of K(+) , Rb(+) and Cs(+) . The product ion spectra of these metaclusters, especially the 15-mer with hetero central ions, indicate that they are higher order assemblies of thymine quintets. The collision-induced dissociation experiments show that the gas-phase stabilities of these metaclusters depend on the size of the central ions, following the order Cs(+) > Rb(+) > K(+) , while K(+) leads to the highest dissociation energy of a thymine quintet. The optimized structures of thymine quintet and decamer were provided by density functional theory calculations, which showed that thymine quintet is bowl-shaped and its tilting angle increases with the size of the central ion. Furthermore, the chirality of thymine quintet was defined for the first time and the resulting different diastereoisomers of thymine decamers were also revealed by the calculation study. Copyright © 2011 John Wiley & Sons, Ltd.

Decameric uracil complexes around Li+

Electrospray ionization (ESI) in combination with mass spectrometry (MS) experiments were carried out to study decameric uracil complexes cationized with Li(+) ion. A previous study has shown that, under specific experimental conditions, a particularly intense peak of the decamer U(10)Li(+) is formed, which was referred to as an indication for so-called 'magic number' cluster. In order to gain more insight on the structure of this decameric complex, here, we report experimental studies concerning the kinetics of the fragmentation. In accordance with the new experimental data, structural models were constructed and fully optimized using ab initio and density functional theory quantum chemistry calculations. The theoretical study allowed us to propose a stable gas-phase structure which is compatible with all experimental findings.

Tailoring the properties of guanosine-based supramolecular hydrogels

We demonstrate that very stable hydrogels can be formed in aqueous potassium chloride solution by mixing a well-known gelator (guanosine, G) with a nongelator of similar structure (2',3',5'-tri-O-acetylguanosine, TAcG), and through a variety of characterization methods including rheology, small-angle neutron scattering, differential scanning calorimetry, and atomic force microscopy, we report substantial progress toward elucidating the factors that control the structure and stability of this fibrous gel system. The results suggest that the tailorability, long lifetime stability, and thermomechanical behavior of these gels derives from a reduction in the driving force toward crystallization with increased hydrophobic (TAcG) content, accompanied by a simultaneous decrease in fiber length and an increase in fiber width.

Formation of complexes between uracil and calcium ions: an ESI/MS/MS study in combination with theoretical calculations

Cationized uracil clusters around calcium metal ions were generated in the gas phase by electrospray ionization (ESI). A previous study showed that with particular experimental conditions, hexamer, octamer, decamer, dodecamer and tetradecamer uracil clusters are present in high quantities. New experiments were carried out to understand the reasons for the particular stability of these complexes. MS/MS experiments suggested that these uracil clusters belong to the same family. Based on ab initio and DFT quantum chemistry calculations, structures in agreement with experimental results are proposed for these clusters.

Theoretical and experimental studies of cationized uracil complexes in the gas phase

Cationized uracil clusters were generated in the gas phase by electrospray ionization (ESI). Mass spectrometry experiments showed that with particular experimental conditions, decameric uracil clusters are magic number clusters. MS/MS experiments demonstrated that the structure of these decameric uracil clusters depends substantially on the size and the charge of the cation. On the basis of the ab initio and density functional theory (DFT) quantum chemistry calculations, structures for these decameric clusters were proposed. These structures are in agreement with the experimental mass spectra of modified nucleobases. Theoretical calculations showed that complexes experimentally observed using ESI-MS techniques, are not naturally the most stable in the gas phase.

Molecular structure of guanine-quartet supramolecular assemblies in a gel-state based on a DFT calculation of infrared and vibrational circular dichroism spectra

The infrared (IR) and vibrational circular dichroism (VCD) spectra of guanosine-5'-hydrazide ( G-1), a powerful hydrogelator, have been measured and analyzed on the basis of ab initio modeling. B3LYP/6-31G** DFT calculations predict that G-1, forming a clear solution in deuterated DMSO, is present in monomeric form in this solvent, whereas strong gelation in a phosphate buffer is due to the formation of a guanine-quartet structure, ( G-1)4, in which the four G-1 are linked by hydrogen-bonded guanine moieties and stabilized by an alkali metal cation. The B3LYP/6-31G** IR and VCD spectra of the nearly planar G-quartet, whose structure is slightly distorted from the C4h symmetry, in which the G-bases interact via four Hoogsteen-type hydrogen bonds and a sodium cation is positioned in the middle of the G-quartet, are in very good agreement with the experimental spectra, indicating that this structure is the predominant structure in the gel state. The geometric parameters are discussed. This study is the first to use IR and VCD spectroscopies coupled with DFT calculations to elucidate the structure of a supramolecular species in a gel state and shows the VCD spectroscopy as a powerful method for investigating the structure of complex supramolecular self-assemblies where the use of other structural methods is limited.