Ab Initio Study of Stacking Interactions in A- and B-DNA

2'-O,4'-C-Methylene-Bridged Nucleic Acids Stabilize Metal-Mediated Base Pairing in a DNA Duplex

The 2'-O,4'-C-methylene-bridged or locked nucleic acid (2',4'-BNA/LNA), with an N-type sugar conformation, effectively improves duplex-forming ability. 2',4'-BNA/LNA is widely used to improve gene knockdown in nucleic acid based therapies and is used in gene diagnosis. Metal-mediated base pairs (MMBPs), such as thymine (T)-Hg^II -T and cytosine (C)-Ag^I -C have been developed and used as attractive tools in DNA nanotechnology studies. This study aimed to investigate the application of 2',4'-BNA/LNA in the field of MMBPs. 2',4'-BNA/LNA with 5-methylcytosine stabilized the MMBP of C with Ag^I ions. Moreover, the 2',4'-BNA/LNA sugar significantly improved the duplex-forming ability of the DNA/DNA complex, relative to that by the unmodified sugar. These results suggest that the sugar conformation is important for improving the stability of duplex-containing MMBPs. The results indicate that 2',4'-BNA/LNA can be applied not only to nucleic acid based therapies, but also to MMBP technologies.

Optical and Structural Characterization of a Chronic Myeloid Leukemia DNA Biosensor

Selective base pairing is the foundation of DNA recognition. Here, we elucidate the molecular and structural details of a FRET-based two-component molecular beacon relying on steady-state fluorescence spectroscopy, small-angle X-ray scattering (SAXS), microscale thermophoresis (MST), and differential electrophoretic mobility. This molecular beacon was designed to detect the most common fusion sequences causing chronic myeloid leukemia, e14a2 and e13a2. The emission spectra indicate that the self-assembly of the different components of the biosensor occurs sequentially, triggered by the fully complementary target. We further assessed the structural alterations leading to the specific fluorescence FRET signature by SAXS, MST, and the differential electrophoretic mobility, where the size range observed is consistent with hybridization and formation of a 1:1:1 complex for the probe in the presence of the complementary target and revelator. These results highlight the importance of different techniques to explore conformational DNA changes in solution and its potential to design and characterize molecular biosensors for genetic disease diagnosis.

Geometric Patterns for Neighboring Bases Near the Stacked State in Nucleic Acid Strands

Structural variation in base stacking has been analyzed frequently in isolated double helical contexts for nucleic acids, but not as often in nonhelical geometries or in complex biomolecular environments. In this study, conformations of two neighboring bases near their stacked state in any environment are comprehensively characterized for single-strand dinucleotide (SSD) nucleic acid crystal structure conformations. An ensemble clustering method is used to identify a reduced set of representative stacking geometries based on pairwise distances between select atoms in consecutive bases, with multiple separable conformational clusters obtained for categories divided by nucleic acid type (DNA/RNA), SSD sequence, stacking face orientation, and the presence or absence of a protein environment. For both DNA and RNA, SSD conformations are observed that are either close to the A-form, or close to the B-form, or intermediate between the two forms, or further away from either form, illustrating the local structural heterogeneity near the stacked state. Among this large variety of distinct conformations, several common stacking patterns are observed between DNA and RNA, and between nucleic acids in isolation or in complex with proteins, suggesting that these might be stable stacking orientations. Noncanonical face/face orientations of the two bases are also observed for neighboring bases in the same strand, but their frequency is much lower, with multiple SSD sequences across categories showing no occurrences of such unusual stacked conformations. The resulting reduced set of stacking geometries is directly useful for stacking-energy comparisons between empirical force fields, prediction of plausible localized variations in single-strand structures near their canonical states, and identification of analogous stacking patterns in newly solved nucleic acid containing structures.

Dynamical treatment of charge transfer through duplex nucleic acids containing modified adenines

We address the issue of whether chemical alterations of nucleobases are an effective tool to modulate charge transfer through DNA molecules. Our investigation uses a multilevel computational approach based on classical molecular dynamics and quantum chemistry. We find yet another piece of evidence that structural fluctuations are a key factor to determine the electronic structure of double-stranded DNA. We argue that the electronic structure and charge transfer ability of flexible polymers is the result of a complex intertwining of various structural, dynamical and chemical factors. Chemical intuition may be used to design molecular wires, but this is not the sole component in the complex charge transfer mechanism through DNA.

Theoretical investigation of the coupling between hydrogen-atom transfer and stacking interaction in adenine-thymine dimers

Three different dimers of the adenine-thymine (A-T) base pair are studied to point out the changes of important properties (structure, atomic charge, energy and so on) induced by coupling between the movement of the atoms in the hydrogen bonds and the stacking interaction. The comparison of these results with those for the A-T monomer system explains the role of the stacking interaction in the hydrogen-atom transfer in this biologically important base pair. The results support the idea that this coupling depends on the exact dimer considered and is different for the N-N and N-O hydrogen bonds. In particular, the correlation between the hydrogen transfer and the stacking interaction is more relevant for the N-N bridge than for the N-O one. Also, the two different mechanisms of two-hydrogen transfer (step by step and concerted) can be modified by the stacking interaction between the base pairs.

Impact of geometry optimization on base-base stacking interaction energies in the canonical A- and B-forms of DNA

Base stacking is known to make an important contribution to the stability of DNA and RNA, and accordingly, significant efforts are ongoing to calculate stacking energies using ab initio quantum mechanical methods. To date, impressive improvements have been made in the model chemistries used to perform stacking energy calculations, including extensions that include robust treatments of electron correlation with extended basis sets, as required to treat interactions where dispersion makes a significant contribution. However, those efforts typically use rigid monomer geometries when calculating the interaction energies. To overcome this, in the present work, we describe a novel internal coordinate definition that allows the relative, intermolecular orientation of stacked base monomers to be constrained during geometry optimizations while allowing full optimization of the intramolecular degrees of freedom. Use of the novel reference frame to calculate the impact of full geometry optimization versus constraining the bases to be planar on base monomer stacking energies, combined with density-fitted, spin-component scaling MP2 treatment of electron correlation, shows that full optimization makes the average stacking energy more favorable by -3.4 and -1.5 kcal/mol for the canonical A and B conformations of the 16 5' to 3' base stacked monomers. Thus, treatment of geometry optimization impacts the stacking energies to an extent similar to or greater than the impact of current state of the art increases in the rigor of the model chemistry itself used to treat base stacking. Results also indicate that stacking favors the B-form of DNA, though the average difference versus the A-form decreases from -2.6 to -0.6 kcal/mol when the intramolecular geometry is allowed to fully relax. However, stacking involving cytosine is shown to favor the A-form of DNA, with that contribution generally larger in the fully optimized bases. The present results show the importance of allowing geometry optimization, as well as properly treating the appropriate model chemistry, in studies of nucleic acid base stacking.

yDNA versus yyDNA pyrimidines: computational analysis of the effects of unidirectional ring expansion on the preferred sugar-base orientation, hydrogen-bonding interactions and stacking abilities

The properties of natural, y- and yy-pyrimidines are compared using computational (B3LYP, MP2) methods. Ring expansion upon incorporation of benzene or naphthalene into the natural pyrimidines affects the preferred orientation of the base about the glycosidic bond in the corresponding nucleoside to a similar extent. Specifically, although the natural pyrimidines preferentially adopt the anti orientation with respect to the 2'-deoxyribose moiety, the expanded analogues will likely display (anti/syn) conformational heterogeneity, which may lead to alternate hydrogen-bonding modes in double-stranded duplexes. Nevertheless, the A:T Watson-Crick hydrogen-bond strengths do not significantly change upon base expansion, while the G:C interaction energy is slightly strengthened upon incorporation of either expanded pyrimidine. The largest effect of base expansion occurs in the stacking energies. Specifically, the maximum (most negative) stacking energies in isolated dimers formed by aligning the nucleobase centers of mass can be increased up to 45% by inclusion of a single y-pyrimidine and up to 55% by consideration of a yy-pyrimidine. Similar increases in the stacking interactions are found when a simplified duplex model composed of two stacked (hydrogen-bonded) base pairs is considered, where both the intrastrand and interstrand stacking interactions can be increased and the effects are more pronounced for the yy-pyrimidines. Moreover, the total stability (sum of all hydrogen-bonding and stacking interactions) is greater for duplexes containing expanded yy-pyrimidines compared to y-pyrimidines, which is mainly due to enhanced stacking interactions. Thus, our calculations suggest that multiple unidirectional increases in the size of the nucleobase spacer can continuously enhance the stability of expanded duplexes.

A computational proposal for the experimentally observed discriminatory behavior of hypoxanthine, a weak universal nucleobase

A computational model composed of six nucleobases was used to investigate why hypoxanthine does not yield duplexes of equal stability when paired opposite each of the natural DNA nucleobases. The magnitudes of all nearest-neighbor interactions in a DNA helix were calculated, including hydrogen-bonding, intra- and interstrand stacking interactions, as well as 1-3 intrastrand stacking interactions. Although the stacking interactions in DNA relevant arrangements are significant and account for at least one third of the total stabilization energy in our nucleobase complexes, the trends in the magnitude of the stacking interactions cannot explain the relative experimental melting temperatures previously reported in the literature. Furthermore, although the total hydrogen-bonding interactions explain why hypoxanthine preferentially pairs with cytosine, the experimental trend for the remaining nucleobases (A, T, G) is not explained. In fact, the calculated pairing preference of hypoxanthine matches that determined experimentally only when the sum of all types of nearest-neighbor interactions is considered. This finding highlights a strong correlation between the relative magnitude of the total nucleobase-nucleobase interactions and measured melting temperatures for DNA strands containing hypoxanthine despite the potential role of other factors (including hydration, temperature, sugar-phosphate backbone). By considering a large range of sequence combinations, we reveal that the binding preference of hypoxanthine is strongly dependent on the nucleobase sequence, which may explain the varied ability of hypoxanthine to universally bind to the natural bases. As a result, we propose that future work should closely examine the interplay between the dominant nucleobase-nucleobase interactions and the overall strand stability to fully understand how sequence context affects the universal binding properties of modified bases and to aid the design of new molecules with ambiguous pairing properties.

COMPUTATION OF LARGE SYSTEMS WITH AN ECONOMIC BASIS SET: AB INITIO CALCULATIONS OF BIOLOGICAL NUCLEIC ACID BASE PAIRS

We show that an economic basis set, in which the polarization functions are considered only for oxygen and nitrogen atoms of strong electronegativity, can be used to determine reliable structures of nucleic acid base pairs. Mulliken charge analysis and the HOMO-LUMO gap in single-point energy calculations using standard basis sets on the geometric structures optimized with the economic basis set found reasonable agreements with the ones of standard calculations. This study is expected to provide a general guideline for basis set selection in the computation of large biological systems being performed with considerable high accuracy, using a low cost computation resource.

INTERACTION OF EMODIN WITH DNA BASES: A DENSITY FUNCTIONAL THEORY

In this study, we present work on the physicochemical interaction between the anticancer drug molecule Emodin (ED) and DNA. Comprehending the physicochemical properties of this drug besides the mechanism by which it interacts with DNA should eventually permit the rational design of novel anticancer or antiviral drugs. The final purpose is the clarification of this novel class of drugs as potential pharmaceutical agents. The properties of the isolated intercalator ED and its stacking interactions with adenine⋯thymine (AT) and guanine⋯cytosine (GC) (nucleic acid base pairs) in face-to-face and face-to-back models were studied by means of the density functional tightbinding (DFTB) method. This method was an approximate version of the density functional theory (DFT) method and it includes London dispersion energy. The molecular modeling of the complex formed between ED and DNA indicated that this complex was capable of contributing to the formation of a constant intercalation site. The results exhibit that ED changes affect DNA structure with reference to bond lengths, bond angles, torsion angles, and charges.

Computational model for predicting experimental RNA and DNA nearest-neighbor free energy rankings

Hydrogen-bonding, intrastrand base-stacking, and interstrand base-stacking energies were calculated for RNA and DNA dimers at the MP2(full)/6-311G** level of theory. Standard A-form RNA and B-form DNA geometries from average fiber diffraction data were employed for all base monomer and dimer geometries, and all dimer binding energies were obtained via single-point calculations. The effects of water solvation were considered using the PCM model. The resulting dimer binding energies were used to calculate the 10 unique RNA and 10 unique DNA computational nearest-neighbor energies, and the ranking of these computational nearest neighbor energies are in excellent agreement with the ranking of the experimental nearest-neighbor free energies. These results dispel the notion that average fiber diffraction geometries are insufficient for calculating RNA and DNA stacking energies.

Comparison of intrinsic stacking energies of ten unique dinucleotide steps in A-RNA and B-DNA duplexes. Can we determine correct order of stability by quantum-chemical calculations?

High level ab initio methods have been used to study stacking interactions in ten unique base pair steps both in A-RNA and in B-DNA duplexes. The protocol for selection of geometries based on molecular dynamics (MD) simulations is proposed, and its suitability is demonstrated by comparison with stacking in steps at fiber diffraction geometries. It is shown that fiber diffraction geometries are not sufficiently accurate for interaction energy calculations. In addition, the protocol for selection of geometries based on MD simulations allows for the evaluation of the variability of the intrinsic stacking energies along the MD trajectories. The uncertainty in stacking energies (difference between the most and least stable geometry) due to the dynamical nature of systems can be, in some cases, as large as 3.0 kcal x mol(-1), which is almost 50% of the actual sequence dependence of base stacking energies (the energy difference between the most and least stable sequences). Thus, assessing the relative magnitude of the gas phase stacking energy using a single geometry for each sequence is insufficient to obtain an unambiguous order of gas phase stacking energies in canonical double helices. Though the ordering of ten unique dinucleotide steps cannot be definitive, some general conclusions were drawn. The stacking energies of base pair steps in A-RNA are more evenly separated compared to B-DNA, and their ordering is less sensitive to the dynamics of the system compared to be B-DNA. The most stable step both in B-DNA and A-RNA is the GC/GC [corrected] step that is well separated from the second most stable step CG/CG. [corrected] Also the least stable step (the CC/GG step) is well separated from the rest of the structures. The calculations further show that B-DNA stacking is favorable only marginally (on average by 1.14 kcal x mol(-1) per base pair step) over A-RNA stacking, and this difference vanishes after subtracting the stabilizing van der Waals effect of the thymine 5-methyl group that is absent in RNA. Basically, no correlation between the sequence dependence of gas phase stacking energies and the sequence dependence of DeltaG degrees(37) free energies used in nearest-neighbor models was found either for B-DNA or for A-RNA. This reflects the complexity of the balance of forces that are responsible for the sequence dependence of thermodynamics stability of nucleic acids, which masks the effect of the intrinsic interactions between the stacked base pairs.

Force field parameters for rotation around chi torsion axis in nucleic acids

To raise the accuracy of the force field for nucleic acids, several parameters were elaborated, focusing on the rotation around chi torsion axis. The reliability of molecular dynamics (MD) simulation was significantly increased by improving the torsion parameters at C8--N9--C1'--X (X = H1', C2', O4') in A, G and those at C6--N1--C1'--X in C, T, and U. In this work, we constructed small models representing the chemical structure of A, G, C, T, and U, and estimated energy profile for chi-axis rotation by executing numerous quantum mechanical (QM) calculations. The parameters were derived by discrete Fourier transformation of the calculated QM data. A comparison in energy profile between molecular mechanical (MM) calculation and QM one shows that our presently derived parameters well reproduce the energy surface of QM calculation for all the above torsion terms. Furthermore, our parameters show a good performance in MD simulations of some nucleic acids. Hence, the present refinement of parameters will enable us to perform more accurate simulations for various types of nucleic acids.

A theoretical study on interactions between mitoxantrone as an anticancer drug and DNA: application in drug design

This research is an effort to further understand the physicochemical interaction between the novel drug, mitoxantrone (MTX) and its biologic receptor, DNA. The ultimate goal is to design drugs that interact more with DNA. Understanding the physicochemical properties of the drug as well as the mechanism by which it interacts with DNA, it should ultimately allow the rational design of novel anti-cancer or anti-viral drugs. Molecular modelling on the complex formed between MTX and DNA presented that this complex was indeed fully capable of participating in the formation of a stable intercalation site. Furthermore, the molecular geometries of MTX and the DNA bases (adenine, guanine, cytosine and thymine) were optimized with the aid of the B3LYP/6-31G* method. The properties of the isolated intercalator and its stacking interactions with the adenine...thymine (AT) and guanine...cytosine (GC) nucleic acid base pairs were studied with the DFTB method (density functional tight-binding), an approximate version of the DFT method, that was extended to cover the London dispersion energy. The B3LYP/6-31G* stabilization energies of the intercalator...base pair complexes were found 10.06 kcal/mol and 21.64 kcal/mol for AT...MTX and GC...MTX, respectively. It was concluded that the dispersion energy and the electrostatic interaction contributed to the stability of the intercalator.DNA base pair complexes. The results concluded from the comparison of the DFTB method and the Hartree-fock method point out that these methods show close results and support each other.

Refinement of the AMBER force field for nucleic acids: improving the description of alpha/gamma conformers

We present here the parmbsc0 force field, a refinement of the AMBER parm99 force field, where emphasis has been made on the correct representation of the alpha/gamma concerted rotation in nucleic acids (NAs). The modified force field corrects overpopulations of the alpha/gamma = (g+,t) backbone that were seen in long (more than 10 ns) simulations with previous AMBER parameter sets (parm94-99). The force field has been derived by fitting to high-level quantum mechanical data and verified by comparison with very high-level quantum mechanical calculations and by a very extensive comparison between simulations and experimental data. The set of validation simulations includes two of the longest trajectories published to date for the DNA duplex (200 ns each) and the largest variety of NA structures studied to date (15 different NA families and 97 individual structures). The total simulation time used to validate the force field includes near 1 mus of state-of-the-art molecular dynamics simulations in aqueous solution.

Polarizable model potential function for nucleic acid bases

A polarizable model potential (PMP) function for adenine (A), cytosine (C), guanine (G), thymine (T), and uracil (U) is developed on the basis of ab initio molecular orbital calculations at the MP2/6-31+G* level. The PMP function consists of Coulomb, van der Waals, and polarization terms. The permanent atomic charges of the Coulomb term are determined by using electrostatic potential (ESP) optimization. The multicenter polarizabilities of the polarization term are determined by using polarized one-electron potential (POP) optimization in which the electron density changes induced by a test charge are target. Isotropic and anisotropic polarizabilities are adopted as the multicenter polarizabilities. In the PMP calculations using the optimized parameters, the interaction energies of Watson-Crick type A-T and C-G base pairs were -15.6 and -29.4 kcal/mol, respectively. The interaction energy of Hoogsteen type A-T base pair was -17.8 kcal/mol. These results reproduce well the quantum chemistry calculations at the MP2/6-311++G(3df,2pd) level within the differences of 0.6 kcal/mol. The stacking energies of A-T and C-G were -9.7 and -10.9 kcal/mol. These reproduce well the calculation results at the MP2/6-311++G (2d,2p) level within the differences of 1.3 kcal/mol. The potential energy surfaces of the system in which a sodium ion or a chloride ion is adjacent to the nucleic acid base are calculated. The interaction energies of the PMP function reproduced well the calculation results at the MP2/6-31+G* or MP2/6-311++G(2d,2p) level. The reason why the PMP function reproduces well the high-level quantum mechanical interaction energies is addressed from the viewpoint of each energy terms.

How Do Size-Expanded DNA Nucleobases Enhance Duplex Stability? Computational Analysis of the Hydrogen-Bonding and Stacking Ability of xDNA Bases

Computational chemistry (B3LYP, MP2) is used to study the properties of size-expanded DNA nucleobases generated by inserting a benzene spacer into the natural nucleobases. Although the addition of the spacer does not significantly affect the hydrogen-bonding properties of natural nucleobases, the orientation of the base about the glycosidic bond necessary for Watson-Crick binding is destabilized, which could have implications for the selectivity of expanded bases, as well as the stability of expanded duplexes. Consideration of the (stacked) binding energies in the preferred relative orientation of natural and expanded nucleobases aligned according to their centers of mass reveals that the stacking within natural dimers can be increased by up to 50% upon expansion of one nucleobase and up to 90% upon expansion of two nucleobases. The implications of these findings to the stability of expanded duplexes were revealed by considering simplified models of natural and mixed duplexes composed of four nucleobases. Although intra- and interstrand interactions within double helices are typically less than those predicted when nucleobases are stacked according to their centers of mass, some nucleobases utilize their full stacking potential within double helices, where both intra- and interstrand interactions can be significant. Most importantly, increasing the size of nucleobases within the duplex significantly increases both intra- and interstrand stacking interactions. Specifically, some interactions are double the magnitude of the corresponding intrastrand interactions in natural helices, and even greater increases in interstrand interactions are sometimes found. Thus, our work suggests that mixed duplexes composed of natural bases hydrogen bound to expanded bases may exploit the increase in the inherent stacking ability of the expanded bases in more than one way and thereby afford duplexes with greater stability than natural DNA.

Effects of Hydrogen-Bonding and Stacking Interactions with Amino Acids on the Acidity of Uracil

The effects of hydrogen-bonding interactions with amino acids on the (N1) acidity of uracil are evaluated using (B3LYP) density functional theory. Many different binding arrangements of each amino acid to three uracil binding sites are considered. The effects on the uracil acidity are found to significantly depend upon the nature of the amino acid and the binding orientation, but weakly depend on the binding site. Our results reveal that in some instances small models for the amino acids can be used, while for other amino acids larger models are required to properly describe the binding to uracil. The gas-phase acidity of uracil is found to increase by up to approximately 60 kJ mol(-1) due to discrete hydrogen-bonding interactions. Although (MP2) stacking interactions with aromatic amino acids decrease the acidity of uracil, unexpected increases in the acidity are found when any of the aromatic amino acids, or the backbone, hydrogen bond to uracil. Consideration of enzymatic and aqueous environments leads to decreases in the effects of the amino acids on the acidity of uracil. However, we find that the magnitude of the decrease varies with the nature of the molecule bound, as well as the (gas-phase) binding orientations and strengths, and therefore solvation effects should be considered on a case-by-case basis in future work. Nevertheless, the effects of amino acid interactions within enzymatic environments are as much as approximately 35 kJ mol(-1). The present study has general implications for understanding the nature of active site amino acids in enzymes, such as DNA repair enzymes, that catalyze reactions involving anionic nucleobase intermediates.

A computational characterization of the hydrogen-bonding and stacking interactions of hypoxanthine

The hydrogen-bonding and stacking interactions of hypoxanthine, a potential universal nucleobase, were calculated using a variety of methodologies (CCSD(T), MP2, B3LYP, PWB6K, AMBER). All methods predict that the hydrogen-bonding interaction in the hypoxanthine-cytosine pair is approximately 25 kJ mol(-1) stronger than that in the other dimers. Although the calculations support suggestions from experiments that hypoxanthine preferentially binds with cytosine, the trend in the calculated hydrogen-bond strengths for the remaining natural nucleobases do not show a strong correlation with the experimentally predicted binding preferences. However, our calculations suggest that the stacking interactions of hypoxanthine are similar in magnitude to the hydrogen-bonding interactions at all levels of theory (with the exception of B3LYP, which incorrectly predicts stacked dimers to be unstable). Therefore, stacking interactions should also be considered when analyzing the stability of DNA helices containing hypoxanthine and the use of larger models that account for both hydrogen-bonding and stacking within DNA duplexes will likely result in better agreement with experimental observations. For the majority of the dimers, PWB6K and AMBER provide reasonable binding strengths at reduced computational costs, and therefore will be useful techniques for considering larger models.

Computation of large systems with an economic basis set: Structures and reactivity indices of nucleic acid base pairs from density functional theory

We show here that an economic basis set can describe nucleic acid base pairs involving the hydrogen bond interactions in density functional calculations. The economic basis set in which the polarization function is added only to oxygen and nitrogen atoms of strong electronegativity can predict reliable geometric structures and dipole moment of nucleic acid base pairs, comparable to those obtained from the basis set of 6-31G* in B3LYP calculations. Combining single point calculations with the standard basis set on the geometric structures optimized by the economic basis set, the present approach has predicted accurate natural bond orbital charge, binding energy, electronegativity, hardness, softness, and electrophilicity index. The principle for basis selection presented in this study can be regarded as a general guideline in the computation of large biological systems with considerably high accuracy and low computational expense.

Syntheses, pi-stacking interactions and base-pairings of uracil pyridinium salts and uracilyl betaines with nucleobases

Reaction of 6-chlorouracil with 4-(dimethylamino)pyridine, 4-methylpyridine, and pyridin-4-yl-morpholine yielded pyridinium-substituted uracils as chlorides which were converted into pyridinium uracilates by deprotonation. These heterocyclic mesomeric betaines are cross-conjugated and thus possess separate cationic (pyridinium) and anionic (uracilate) moieties. Calculations and X-ray single crystal analyses were performed in order to characterize these systems and to compare the salts with the betaines. (1)H NMR experiments in D(2)O proved pi-interactions between the uracilyl betaines and adenine, adenosine, as well as adeninium. No pi-stacking interactions were detected between the betaines and guanosine. The acidic N8-H group of the uracil pyridinium salts caused acid-base reactions which were observed in parallel to pi-stacking interactions. Self-complementarity of the modified uracils was detected by (1)H NMR experiments in DMSO-d(6) and electrospray ionisation mass spectrometry (ESIMS). Ab initio calculations predicted base-pairings of the modified uracils with adeninium, cytosine, and guanine. Several geometries of hydrogen-bonded associates were calculated. Hoogsteen pairings between the uracil-4-(dimethylamino)pyridinium salt and adeninium, as well as associates between the corresponding betaine plus cytosine, and the betaine plus guanine were calculated, and the most stable conformations were determined. In the ESI mass spectra, prominent peaks of associates between the modified uracils and adeninium, cytosine, cytidine, guanosine and d(CpGp) were detected.

Intra- and intermolecular interactions between cyclic-AMP receptor protein and DNA:Ab initio fragment molecular orbital study

The ab initio fragment molecular orbital (FMO) calculations were performed for the cAMP receptor protein (CRP) complexed with a cAMP and DNA duplex to elucidate their sequence-specific binding and the stability of the DNA duplex, as determined by analysis of their inter- and intramolecular interactions. Calculations were performed with the AMBER94 force field and at the HF and MP2 levels with several basis sets. The interfragment interaction energies (IFIEs) were analyzed for interactions of CRP-cAMP with each base pair, DNA duplex with each amino acid residue, and each base pair with each residue. In addition, base-base interactions were analyzed including hydrogen bonding and stacking of DNA. In the interaction between DNA and CRP-cAMP, there was a significant charge transfer (CT) from the DNA to CRP, and this CT interaction played an important role as well as the electrostatic interactions. It is necessary to apply a quantum mechanical approach beyond the "classical" force-field approach to describe the sequence specificity. In the DNA intramolecular interaction, the dispersion interactions dominated the stabilization of the base-pair stacking interactions. Strong, attractive 1,2-stacking interactions and weak, repulsive 1,3-stacking interactions were observed. Comparison of the intramolecular interactions of free and complexed DNA revealed that the base-pairing interactions were stronger, and the stacking interactions were weaker, in the complexed structure. Therefore, the DNA duplex stability appears to change due to both the electrostatic and the CT interactions that take place under conditions of DNA-CRP binding.

The binding of guanine–guanine mismatched DNA to naphthyridine dimer immobilized sensor surfaces: kinetic aspects

Naphthyridine dimer composed of two naphthyridine chromophores and a linker connecting them strongly, and selectively, binds to the guanine-guanine mismatch in duplex DNA. The kinetics for the binding of the G-G mismatch to the naphthyridine dimer was investigated by surface plasmon resonance assay. The sensor surface was prepared by immobilizing naphthyridine dimer through a long poly(ethylene oxide) linker with the ligand density of 9.1 x 10(-12) fmolnm(-2). The kinetic analyses revealed that the binding of the G-G mismatch was sequence dependent on the flanking base pairs, and the G-G mismatches flanking at least one G-C base pair bound to the surface via a two-step process with a 1:1 DNA-ligand stoichiometry. The first association rate constant for the binding of the G-G mismatch in the 5'-CGG-3'/3'-GGC-5' sequence to the naphthyridine dimer-immobilized sensor surface was 3.2 x 10(3)M(-1)s(-1) and the first dissociation rate constant was 1.4 x 10(-2)s(-1). The association and dissociation rate constants for the second step were insensitive to the flanking sequences, and were almost of the same order of magnitude as the first dissociation rate constant. This indicates that the second step had only a small energetic contribution to the binding. The association constant calculated from kinetic parameters was 2.7 x 10(5)M(-1), which is significantly smaller than the apparent association constants obtained from experiments in solution. Electrospray ionization time-of-flight (ESI-TOF) mass spectrometry on the complex produced from the G-G mismatch and naphthyridine dimer showed the formation of the 1:1 complex and a 1:2 DNA-ligand complex in solution. The latter complex became the dominant complex when a six-fold excess of naphthyridine dimer was added to DNA.

The SPR sensor detecting cytosine[bond]cytosine mismatches

We have synthesized the first surface plasmon resonance (SPR) sensor that detects cytosine-cytosine (C[bond]C) mismatches in duplex DNA by immobilizing aminonaphthyridine dimer on the gold surface. The ligand consisting of two 2-aminonaphthyridine chromophores and an alkyl linker connecting them strongly stabilized the C[bond]C mismatches regardless of the flanking sequences. The fully matched duplexes were not stabilized at all under the same conditions. The C[bond]T, C[bond]A, and T[bond]T mismatches were also stabilized with a reduced efficiency. SPR analyses of mismatch-containing 27-mer duplexes were performed with the sensor surface on which the aminonaphthyridine dimer was immobilized. The response for the C[bond]C mismatch in 5'-GCC-3'/3'-CCG-5' was about 83 times stronger than that obtained for the fully matched duplex. The sensor successfully detects the C[bond]C mismatch at the concentration of 10 nM. SPR responses are proportional to the concentration of the C[bond]C mismatch in a range up to 200 nM. Aminonaphthyridine dimer could bind strongly to the C[bond]C mismatches having 10 possible flanking sequences with association constants in the order of 10(6) M(-1). The facile protonation of 2-aminonaphthyridine chromophore at pH 7 producing the hydrogen-bonding surface complementary to that of cytosine was most likely due to the remarkably high selectivity of 1 to the C[bond]C mismatch.

Molecular Interactions of Nucleic Acid Bases. A Review of Quantum-Chemical Studies

Ab initio quantum-chemical calculations with inclusion of electron correlation significantly contributed to our understanding of molecular interactions of DNA and RNA bases. Some of the most important findings are introduced in the present overview: structures and energies of hydrogen bonded base pairs, nature of base stacking, interactions between metal cations and nucleobases, nonplanarity of isolated nucleobases and other monomer properties, tautomeric equilibria of nucleobases, out-of-plane hydrogen bonds and amino acceptor interactions. The role of selected molecular interactions in nucleic acids is discussed and representative examples where these interactions occur are given. Also, accuracy of density functional theory, semiempirical methods, distributed multipole analysis and empirical potentials is commented on. Special attention is given to our very recent reference calculations on base stacking and H-bonding. Finally, we briefly comment on the relationship between advanced ab initio quantum-chemical methods and large-scale explicit solvent molecular dynamics simulations of nucleic acids.

The physical basis of nucleic acid base stacking in water

It has been argued that the stacking of adenyl groups in water must be driven primarily by electrostatic interactions, based upon NMR data showing stacking for two adenyl groups joined by a 3-atom linker but not for two naphthyl groups joined by the same linker. In contrast, theoretical work has suggested that adenine stacking is driven primarily by nonelectrostatic forces, and that electrostatic interactions actually produce a net repulsion between adenines stacking in water. The present study provides evidence that the experimental data for the 3-atom-linked bis-adenyl and bis-naphthyl compounds are consistent with the theory indicating that nonelectrostatic interactions drive adenine stacking. First, a theoretical conformational analysis is found to reproduce the observed ranking of the stacking tendencies of the compounds studied experimentally. A geometric analysis identifies two possible reasons, other than stronger electrostatic interactions, why the 3-atom-linked bis-adenyl compounds should stack more than the bis-naphthyl compounds. First, stacked naphthyl groups tend to lie further apart than stacked adenyl groups, based upon both quantum calculations and crystal structures. This may prevent the bis-naphthyl compound from stacking as extensively as the bis-adenyl compound. Second, geometric analysis shows that more stacked conformations are sterically accessible to the bis-adenyl compound than to the bis-naphthyl compound because the linker is attached to the sides of the adenyl groups, but to the ends of the naphthyl groups. Finally, ab initio quantum mechanics calculations and energy decompositions for relevant conformations of adenine and naphthalene dimers support the view that stacking in these compounds is driven primarily by nonelectrostatic interactions. The present analysis illustrates the importance of considering all aspects of a molecular system when interpreting experimental data, and the value of computer models as an adjunct to chemical intuition.