Nature of Nucleic Acid−Base Stacking: Nonempirical ab Initio and Empirical Potential Characterization of 10 Stacked Base Dimers. Comparison of Stacked and H-Bonded Base Pairs

Predicting noncovalent interactions between aromatic biomolecules with London-dispersion-corrected DFT

Within the framework of Kohn-Sham density functional theory, interaction energies of hydrogen bonded and pi-pi stacked supramolecular complexes of aromatic heterocycles, nucleobase pairs, and complexes of nucleobases with the anti-cancer agent ellipticine as well as its derivatives are evaluated. Dispersion-corrected atom-centered potentials (DCACPs) are employed together with a generalized gradient approximation to the exchange correlation functional. For all systems presented, the DCACP calculations are in very good agreement with available post Hartree-Fock quantum chemical results. Estimates of 3-body contributions (<15% of the respective interaction energy) and deformation energies (5-15% of the interaction energy) are given. Based on our results, we predict a strongly bound interaction energy profile for the ellipticine intercalation process with a stabilization of nearly 40 kcal/mol (deformation energy not taken into account) when fully intercalated. The frontier orbitals of the intercalator-nucleobase complex and the corresponding non-intercalated nucleobases are investigated and show significant changes upon intercalation. The results not only offer some insights into the systems investigated but also suggest that DCACPs can serve as an effective way to achieve higher accuracy in density functional theory without incurring an unaffordable computational overhead, paving ways for more realistic studies on biomolecular complexes in the condensed phase.

Pi-pi stacking tackled with density functional theory

Through comparison with ab initio reference data, we have evaluated the performance of various density functionals for describing pi-pi interactions as a function of the geometry between two stacked benzenes or benzene analogs, between two stacked DNA bases, and between two stacked Watson-Crick pairs. Our main purpose is to find a robust and computationally efficient density functional to be used specifically and only for describing pi-pi stacking interactions in DNA and other biological molecules in the framework of our recently developed QM/QM approach "QUILD". In line with previous studies, most standard density functionals recover, at best, only part of the favorable stacking interactions. An exception is the new KT1 functional, which correctly yields bound pi-stacked structures. Surprisingly, a similarly good performance is achieved with the computationally very robust and efficient local density approximation (LDA). Furthermore, we show that classical electrostatic interactions determine the shape and depth of the pi-pi stacking potential energy surface.

Vibrational dynamics of DNA. II. Deuterium exchange effects and simulated IR absorption spectra

In Paper I, we studied vibrational properties of normal bases, base derivatives, Watson-Crick base pairs, and multiple layer base pair stacks in the frequency range of 1400-1800 cm(-1). However, typical IR absorption spectra of single- and double-stranded DNA have been measured in D(2)O solution. Consequently, the more relevant bases and base pairs are those with deuterium atoms in replacement with labile amino hydrogen atoms. Thus, we have carried out density functional theory vibrational analyses of properly deuterated bases, base pairs, and stacked base pair systems. In the frequency range of interest, both aromatic ring deformation modes and carbonyl stretching modes appear to be strongly IR active. Basis mode frequencies and vibrational coupling constants are newly determined and used to numerically simulate IR absorption spectra. It turns out that the hydration effects on vibrational spectra are important. The numerically simulated vibrational spectra are directly compared with experiments. Also, the (18)O-isotope exchange effect on the poly(dG):poly(dC) spectrum is quantitatively described. The present calculation results will be used to further simulate two-dimensional IR photon echo spectra of DNA oligomers in the companion Paper III.

An analysis of the electrostatic interaction between nucleic acid bases

Results from several commonly used approximate methods of evaluating electrostatic interactions have been compared to the rigorous, nonexpanded electrostatic energies at both uncorrelated and correlated levels of theory. We examined a number of energy profiles for both hydrogen bonded and stacked configurations of the nucleic acid base pairs. We found that the penetration effects play an extremely important role and the expanded electrostatic energies are significantly underestimated with respect to the ab initio values. Apart from the inability to reproduce the magnitudes of the ab initio electrostatic energy, there are other problems with the available approximate electrostatic models. For example, the distributed multipole analysis, one of the most advanced methods, is extremely sensitive to the basis set and level of theory used to evaluate the multipole moments. Detailed ab initio results are provided that other researchers could use to test their approximate models.

Stabilization energies of the hydrogen-bonded and stacked structures of nucleic acid base pairs in the crystal geometries of CG, AT, and AC DNA steps and in the NMR geometry of the 5'-d(GCGAAGC)-3' hairpin: Complete basis set calculations at the MP2 and CCSD(T) levels

Stabilization energies of the H-bonded and stacked structures of a DNA base pair were studied in the crystal structures of adenine-thymine, cytosine-guanine, and adenine-cytosine steps as well as in the 5'-d(GCGAAGC)-3' hairpin (utilizing the NMR geometry). Stabilization energies were determined as the sum of the complete basis set (CBS) limit of MP2 stabilization energies and the Delta E(CCSD(T)) - Delta E(MP2) correction term evaluated with the 6-31G*(0.25) basis set. The CBS limit was determined by a two-point extrapolation using the aug-cc-pVXZ basis sets for X = D and T. While the H-bonding energies are comparable to those of base pairs in a crystal and a vacuum, the stacking energies are considerably smaller in a crystal. Despite this, the stacking is still important and accounts for a significant part of the overall stabilization. It contributes equally to the stability of DNA as does H-bonding for AT-rich DNAs, while in the case of GC-rich DNAs it forms about one-third of the total stabilization. Interstrand stacking reaches surprisingly large values, well comparable to the intrastrand ones, and thus contributes significantly to the overall stabilization. The hairpin structure is characterized by significant stacking, and both guanine...cytosine pairs possess stacking energies larger than 11.5 kcal/mol. A high portion of stabilization in the studied hairpin comes from stacking (similar to that found for AT-rich DNAs) despite the fact that it contains two GC Watson-Crick pairs having very large H-bonding stabilization. The DFT/B3LYP/6-31G** method yields satisfactory values of interaction energies for H-bonded structures, while it fails completely for stacking.

Vibrational dynamics of DNA: IV. Vibrational spectroscopic characteristics of A-, B-, and Z-form DNA’s

Linear and nonlinear IR spectroscopic studies of nucleic acids can provide crucial information on solution conformations of DNA double helix and its complex with other molecules. Carrying out density functional theory calculations of A-, B-, and Z-form DNA's, the authors obtained vibrational spectroscopic properties as well as coupling constants between different basis modes. The vibrational couplings that determine the extent of exciton delocalization are strongly dependent on DNA conformation mainly because the interlayer distance between two neighboring base pairs changes with respect to the DNA conformation. The Z-DNA has comparatively small interlayer vibrational coupling constants so that its vibrational spectrum depends little on the number of base pairs, whereas the A-DNA shows a notable dependency on the size. Furthermore, it is shown that a few distinctively different line shape changes in both IR and two-dimensional IR spectra as the DNA conformation changes from B to A or from B to Z can be used as marker bands and characteristic features distinguishing different DNA conformations.

Hydrogen-Atom Abstraction from the Adenine−Uracil Base Pair†

The hydrogen-abstracted radicals from the adenine-uracil (AU) base pair have been studied at the B3LYP/DZP++ level of theory. The A(N9)-U and A-U(N1) radicals, which correspond to hydrogen-atom abstraction at the adenine N9 and uracil N1 atoms, respectively, were predicted to be the two lowest-lying among the nine (AU-H) radicals studied in this study. The removal of the amino hydrogen of the adenine moiety that forms a hydrogen bond with the uracil O4 atom in the AU pair resulted in radical A(N6a)-U, which has the smallest base-pair dissociation energy, 5.9 kcal mol(-1). This radical is more likely to dissociate into the two isolated bases than to recover the hydrogen bond with the O4 atom through N6-H bond rotation along the C6-N6 bond. In general, the radicals generated by C-H bond breaking were higher in energy than those arising from N-H bond cleavage, because the unpaired electrons in the carbon-centered radicals were mainly localized on the carbon atom from which the hydrogen atom was removed. However, the highest-lying radical was found to arise from removal of the N3 hydrogen of uracil. The most remarkable structural feature of this radical is a very short C-H...O distance of 2.094 A, consistent with a substantial hydrogen bond. Although this radical lost the N1...H-N3 hydrogen bond between the two bases, its dissociation energy was predicted to be 12.9 kcal mol(-1), similar to that of the intact AU base pair. This is due to the transfer of electron density from the adenine N1 atom to the uracil N3 atom.

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.

Charge Localization in Stacked Radical Cation DNA Base Pairs and the Benzene Dimer Studied by Self-Interaction Corrected Density-Functional Theory

The incomplete cancellation of the electron self-interaction can be a serious shortcoming of density-functional theory especially when treating odd-electron systems. In this work, several popular and potentially viable correction schemes are applied in order to characterize the electronic structure of stacked molecular pairs, consisting of a neutral molecule and adjacent radical cation, as a function of separation distance. The unphysical sharing of the positive charge between adjacent molecules separated by 6-7 A is corrected for by applying a new empirical scheme proposed by VandeVondele and Sprik [Phys. Chem. Chem. Phys. 2005, 7, 1363] with a unique choice of parameters. This method is subsequently applied to characterize the electronic structure of two neighboring guanines excised from a canonical Arnott B-DNA structure and will be used in future investigations of certain model DNA fibers.

Synthesis of 2'-O-methyl-RNAs incorporating a 3-deazaguanine, and UV melting and computational studies on its hybridization properties

2′-O-Methyl-RNAs incorporating 3-deazaguanine (c³G) were synthesized by use of N,N-diphenylcarbamoyl and N,N-dimethylaminomethylene as its base protecting groups to suppress sheared-type 5′-GA-3′/5′-GA-3′ tandem mismatched base pairing which requires the N³ atom. These modified RNAs hybridized more weakly with the complementary and single mismatch-containing RNAs than the unmodified RNAs. The T_m experiments were performed to clarify the effects of replacement of the fifth G with c³G on stabilization of 2′-O-methyl-(5′-CGGCGAGGAG-3′)/5′-CUCCGAGCCG-3′ and 2′-O-methyl-(5′-CGGGGACGAG-3′)/5′-CUCGGACCCG-3′duplexes, which form sheared-type and face-to-face type 5′-GA-3′/5′-GA-3′ tandem mismatched base pairs, respectively. Consequently, this replacement led to more pronounced destabilization of the former duplex that needs the N³ atom for the sheared-type base pair than the latter that does not need it for the face-to-face type base pair. A similar tendency was observed for 2′-O-methyl-RNA/DNA duplexes. These results suggest that the N³ atom of G plays an important role in stabilization of the canonical G/C base pair as well as the base discrimination and its loss suppressed formation of the undesired sheared-type mismatched base pair. Computational studies based on ab initio calculations suggest that the weaker hydrogen bonding ability and larger dipole moment of c³G can be the origin of the lower T_m.

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.

Electrical conductance of DNA molecules with varied density of itinerant π electrons

The electrical transport of DNA is closely related to the density of itinerant pi electrons because of the strong electron-lattice interaction. The resistivities of two typical DNA molecules [poly(dG)-poly(dC) and lambda-DNA] with varied densities of itinerant pi electrons are calculated. It is found that the dependence of the resistivity on the density of itinerant pi electrons is symmetrical about the half-filling state of itinerant pi electrons in poly(dG)-poly(dC). At the half-filling state, the Peierls phase transition takes place and poly(dG)-poly(dC) has a large resistivity. When the density of itinerant pi electrons departs far from the half-filling state, the resistivity of poly(dG)-poly(dC) becomes small. For lambda-DNA, there is no Peierls phase transition due to the aperiodicity of its base pair arrangement. The resistivity of poly(dG)-poly(dC) decreases with increasing length of the molecular chain, but the resistivity of lambda-DNA increases with increasing length. The conducting mechanisms for poly(dG)-poly(dC) and a few lambda-DNA molecules with varied densities of itinerant pi electrons are analyzed.

Electron density distribution in stacked benzene dimers: a new approach towards the estimation of stacking interaction energies

The potential energy surface for the benzene dimer in stacked conformations (84 points calculated) was computed at the MP2(FC)6-31+G(2d,2p) level of theory. Electron density (ED) distribution computed using the MP2(FC), B3LYP, and Hartree-Fock methods with the same basis set is studied in the frame of topological analysis. It is found that ED topology does not depend on the method of calculation. The values of the ED and its Laplacian in the cage critical point calculated using different methods are determined to be linearly dependent with the slope depending on basis set. Correlation equations based on these properties allow the interaction energy between benzene rings to be predicted with 8% mean relative error in the energy for the given region of the potential energy surface. This provides a new method for the estimation of stacking interaction energy using ED properties calculated with low level quantum-chemical methods.

Nature of base stacking: reference quantum-chemical stacking energies in ten unique B-DNA base-pair steps

Base-stacking energies in ten unique B-DNA base-pair steps and some other arrangements were evaluated by the second-order Møller-Plesset (MP2) method, complete basis set (CBS) extrapolation, and correction for triple (T) electron-correlation contributions. The CBS(T) calculations were compared with decade-old MP2/6-31G*(0.25) reference data and AMBER force field. The new calculations show modest increases in stacking stabilization compared to the MP2/6-31G*(0.25) data and surprisingly large sequence-dependent variation of stacking energies. The absolute force-field values are in better agreement with the new reference data, while relative discrepancies between quantum-chemical (QM) and force-field values increase modestly. Nevertheless, the force field provides good qualitative description of stacking, and there is no need to introduce additional pair-additive electrostatic terms, such as distributed multipoles or out-of-plane charges. There is a rather surprising difference of about 0.1 A between the vertical separation of base pairs predicted by quantum chemistry and derived from crystal structures. Evaluations of different local arrangements of the 5'-CG-3' step indicate a sensitivity of the relative stacking energies to the level of calculation. Thus, describing quantitative relations between local DNA geometrical variations and stacking may be more complicated than usually assumed. The reference calculations are complemented by continuum-solvent assessment of solvent-screening effects.

Quantum mechanical description of the interactions between DNA and water

In recent years, a lot of attention has been focused on the electronic properties of DNA. With recent advances in linear scaling quantum mechanics there are now new tools available to enhance our understanding of the electronic properties of DNA among other biomolecules. Using both explicit solvent models and implicit (continuum) solvent models, the electronic characteristics of a dodecamer duplex DNA have been fully studied using both divide and conquer (D&C), semi-empirical quantum mechanics and non-D&C semi-empirical quantum mechanics. According to the AM1 Hamiltonian, approximately 3.5 electrons (approximately 0.3 electron/base pair) are transferred from the duplex to the solvent. According to the density of state (DOS) analysis, in vacuo DNA has a band gap of approximately 1 eV showing that in the absence of solvent, the DNA may exhibit similar properties to those of a semiconductor. Upon increasing solvation (2.5-5.5 A), the band gap ranges from approximately 3 eV to approximately 6 eV. For the implicit solvent model, the band gap continues this widening trend to approximately 7 eV. Therefore, upon solvation and in the absence of dopants, the DNA should begin to loose its conductive properties. Finally, when one considers the energy and localization of the frontier orbitals (HOMO and LUMO), solvent has a stabilizing effect on the DNA system. The energy of the HOMO drops from approximately 15 eV in vacuo to approximately 2 eV for 5.5 A of water to approximately -8 eV for the implicit solvent model. Similarly, the LUMO drops from approximately 16 eV for in vacuo to approximately 9 eV for 5.5 A of water to approximately -1 eV for the implicit model. Beyond the importance of the computed results on the materials properties of DNA, the present work also shows that the behavior of intercalators will be affected by the electronic properties of DNA. This could have an impact on our understanding of how DNA based drugs interact with DNA and on the design of new DNA based small molecule drugs.

Bioactive principles in the bark of Pilidiostigma tropicum

The crude dichloromethane bark extract of Pilidiostigma tropicum (Myrtaceae) from north Queensland, Australia, shows antibacterial and cytotoxic activity. Bioactivity-directed separation led to the isolation of rhodomyrtoxin B and ursolic acid-3-p-coumarate as the biologically active materials. The structures of these compounds were elucidated on the basis of spectral analysis. The intercalation interaction of rhodomyrtoxin B with DNA was investigated using molecular mechanics and ab initio molecular-orbital techniques. A favorable pi-pi interaction between rhodomyrtoxin B and the cytosine-guanine base pair is predicted, but the orientation of the interaction cannot be predicted based on frontier molecular orbitals.

Influence of Thioketo Substitution on the Properties of Uracil and Its Noncovalent Interactions with Alkali Metal Ions: Threshold Collision-Induced Dissociation and Theoretical Studies†

Experimental and theoretical studies are carried out to determine the influence of thioketo substitution on the properties of uracil and its noncovalent interactions with alkali metal ions. Bond dissociation energies of alkali metal ion-thiouracil complexes, M(+)(SU), are determined using threshold collision-induced dissociation techniques in a guided ion beam mass spectrometer, where M(+) = Li(+), Na(+), and K(+) and SU = 2-thiouracil, 4-thiouracil, 2,4-dithiouracil, 5-methyl-2-thiouracil, and 6-methyl-2-thiouracil. Ab initio electronic structure calculations are performed to determine the structures and theoretical bond dissociation energies of these complexes and provide molecular constants necessary for thermodynamic analysis of the experimental data. Theoretical calculations are also performed to examine the influence of thioketo substitution on the acidities, proton affinities, and A::SU Watson-Crick base pairing energies. In general, thioketo substitution leads to an increase in both the proton affinity and the acidity of uracil. 2-Thio substitution generally results in an increase in the alkali metal ion binding affinities but has almost no affect on the stability of the A::SU base pair. In contrast, 4-thio substitution results in a decrease in the alkali metal ion binding affinities and a significant decrease in the stability of the A::SU base pair. In addition, alkali metal ion binding is expected to lead to an increase in the stability of both single-stranded and double-stranded nucleic acids by reducing the charge on the nucleic acid in a zwitterion effect as well as through additional noncovalent interactions between the alkali metal ion and the nucleobases.

Large-amplitude fast motions in double-stranded DNA driven by solvent thermal fluctuations

The nature of the internal dynamics of double-stranded DNA in aqueous environment remains to be established. We consider the motions to stem from thermal fluctuations/dissipations of the harmonic modes of beads (bases and sugars) in a cylindrical geometry that are tracked through the stochastic Langevin trajectories; these are characterized by parameters obtained from published data. The present approach has allowed a comparative study of the dynamics for DNA lengths in the range of 20-600 base pairs. For this range, we find that rotational motions about directions parallel to the helix axis (opening, twist) and perpendicular to it (propeller-twist, roll) contribute significantly to the dynamics. For a 20-mer at a solvent viscosity of 1 cP, the calculated fluorescence anisotropy profile exhibits a fast decay in the subnanosecond range due to large-amplitude fluctuations at the mesoscopic level. This feature reproduces the experimental behavior well, and suggests a possible way for the initiation of biological processes: they may be suddenly triggered on this scale through the occurrence of favorable thermal fluctuations. This analysis also reveals that, as is the case for a 20-mer, the dynamics of longer N-mers are dominated by internal motions, and are modulated by the viscosity of the solvent, in agreement with our previous experimental observations. Moreover, the model indicates that occurrence of partially concerted rotations of the bases due to thermal fluctuations can possibly be sustained over a DNA length of the order of 100 A at 1 ns, suggesting a possible mechanism for action-at-a-distance in transcription.

The World of Non-Covalent Interactions: 2006

The review focusses on the fundamental importance of non-covalent interactions in nature by illustrating specific examples from chemistry, physics and the biosciences. Laser spectroscopic methods and both ab initio and molecular modelling procedures used for the study of non-covalent interactions in molecular clusters are briefly outlined. The role of structure and geometry, stabilization energy, potential and free energy surfaces for molecular clusters is extensively discussed in the light of the most advanced ab initio computational results for the CCSD(T) method, extrapolated to the CBS limit. The most important types of non-covalent complexes are classified and several small and medium size non-covalent systems, including H-bonded and improper H-bonded complexes, nucleic acid base pairs, and peptides and proteins are discussed with some detail. Finally, we evaluate the interpretation of experimental results in comparison with state of the art theoretical models: this is illustrated for phenol...Ar, the benzene dimer and nucleic acid base pairs. A review with 270 references.

Electrospray ionization of nucleic acid aptamer/small molecule complexes for screening aptamer selectivity

Molecular recognition of small molecule ligands by the nucleic acid aptamers for tobramycin, ATP, and FMN has been examined using electrospray ionization mass spectrometry (ESI-MS). Mass spectrometric data for binding stoichiometry and relative binding affinity correlated well with solution data for tobramycin aptamer complexes, in which aptamer/ligand interactions are mediated by hydrogen bonds. For the ATP and FMN aptamers, where ligand interactions involve both hydrogen bonding and significant pi-stacking, the relative binding affinities determined by MS did not fully correlate with results obtained from solution experiments. Some high-affinity aptamer/ligand complexes appeared to be destabilized in the gas phase by internal Coulombic repulsion. In CAD experiments, complexes with a greater number of intermolecular hydrogen bonds exhibited greater gas-phase stability even in cases when solution binding affinities were equivalent. These results indicate that in at least some cases, mass spectrometric data on aptamer/ligand binding affinities should be used in conjunction with complementary techniques to fully assess aptamer molecular recognition properties.

Potential Energy Surfaces of an Adenine−Thymine Base Pair and Its Methylated Analogue in the Presence of One and Two Water Molecules: Molecular Mechanics and Correlated Ab Initio Study

Potential energy surfaces of monohydrated and dihydrated adenine-thymine and 9-methyladenine-1-methylthymine base pairs were examined by the molecular dynamics/quenching technique using the Cornell et al. force field (J. Am. Chem. Soc. 1995, 117, 5179). Long runs of molecular dynamics/quenching calculations allowed us to evaluate the free energy surface. The most stable and populated structures found were fully reoptimized at the correlated ab initio level employing the resolution of identity Møller-Plesset method. A systematic study of the base pairs' microhydration using both the empirical and the high-level correlated ab initio approaches is presented for the first time. We show that the occurrence of water molecules and their gradually increasing number as well as the methylation of the bases favor stacked structures over the planar hydrogen-bonded ones. These results based on the correlated ab initio calculations are in the excellent agreement with data obtained from our previous empirical potential molecular dynamics study (Kabelác et al. Chem.-Eur. J. 2001, 7, 2067).

Influence of halogenation on the properties of uracil and its noncovalent interactions with alkali metal ions. Threshold collision-induced dissociation and theoretical studies

The influence of halogenation on the properties of uracil and its noncovalent interactions with alkali metal ions is investigated both experimentally and theoretically. Bond dissociation energies of alkali metal ion-halouracil complexes, M+(XU), are determined using threshold collision-induced dissociation techniques in a guided ion beam mass spectrometer, where M+ = Li+, Na+, and K+ and XU = 5-fluorouracil, 5-chlorouracil, 6-chlorouracil, 5-bromouracil, and 5-iodouracil. The structures and theoretical bond dissociation energies of these complexes are determined from ab initio calculations. Theoretical calculations are also performed to examine the influence of halogenation on the acidities, proton affinities, and Watson-Crick base pairing energies. Halogenation of uracil is found to produce a decrease in the proton affinity, an increase in the alkali metal ion binding affinities, an increase in the acidity, and stabilization of the A::U base pair. In addition, alkali metal ion binding is expected to lead to an increase in the stability of nucleic acids by reducing the charge on the nucleic acid in a zwitterion effect as well as through additional noncovalent interactions between the alkali metal ion and the nucleobases.

Insight into the structural role of carotenoids in the photosystem I: a quantum chemical analysis

The structural stabilization role of carotenoids in the formation of photosynthetic pigment-protein complexes is investigated theoretically. The pi-pi stacking and CH-pi interactions between beta-carotenes and their surrounding chlorophylls (and/or aromatic residues) in Photosystem I (PS1) from the cyanobacterium Synechococcus elongatus were studied by means of the supermolecular approach at the level of the second-order Møller-Plesset perturbation method. PS1 features a core integral antenna system consisting of 22 beta-carotenes intertwined with 90 chlorophyll molecules. The binding environments of all 22 beta-carotenes were systematically analyzed. For 21 out of the 22 cases, one or more chlorophyll molecules exist within van der Waals' contacts of the beta-carotene molecule. The calculated strengths of pi-pi stacking interactions between the conjugated core of beta-carotene and the aromatic tetrapyrrole rings of chlorophyll are substantial, ranging from -3.54 kcal/mol for the perpendicular-positioned BCR4004...CHL1217 pair to -16.01 kcal/mol for the parallel-oriented BCR4007...CHL1122 pair. A strong dependence of the pi-pi stacking interaction energies on the intermolecular configurations of the two interacting pi-planes is observed. The parallel-oriented beta-carotene and chlorophyll pair is energetically much more stable than the perpendicular-positioned pair. The larger the extent of pi-pi overlapping, the stronger the interaction strength. In many cases, the beta-ring ends of beta-carotene molecules are found to interact with the tetrapyrrole rings of chlorophyll via CH-pi interactions. For the latter interactions, the calculated interaction strengths vary from -7.03 to -11.03 kcal/mol, depending on the intermolecular configuration. This work leads to the conclusion that pi-pi stacking and CH-pi interactions between beta-carotene and their surrounding chlorophylls and aromatic residues play an essential role in binding beta-carotenes in PS1 from S. elongatus. Consequently, the molecular basis of the structural stabilization function of carotenoids in formation of the photosynthetic pigment-protein complexes is established.

Stair motifs at protein-DNA interfaces: nonadditivity of H-bond, stacking, and cation-pi interactions

At the interface between protein and double-stranded DNA, stair motifs simultaneously involve three different types of pairwise interactions: aromatic base stacking, hydrogen bonding, and cation-pi. The relative importance of these interactions is studied in the stair motif occurring in the 1TC3 crystal structure, which involves an arginine and two stacked guanines, by means of Hartree-Fock (HF) and Møller-Plesset energy and free energy calculations, including vibrational, rotational, translational contributions, both in a vacuum and various solvents. The results obtained show an anti-cooperative tendency of the HF energy and vibrational free energy terms, and the cooperativity of the rotational, translational, and solvation free energies. Hence, the cooperativity of the stair motif interactions, in the context of protein-DNA recognition, can be viewed as arising from the environment.

A local second-order Møller–Plesset method with localized orbitals: A parallelized efficient electron correlation method

Using orthogonal localized occupied orbitals we have developed and implemented a parallelized local second-order Møller-Plesset (MP2) method based on the idea developed by Head-Gordon and co-workers. A subset of nonorthogonal correlation functions (the orbital domain) was assigned to each of the localized occupied orbitals using a distance criterion and excitations from localized occupied orbitals that were arranged into subsets. The correlation energy was estimated using a partial diagonalization and an iterative efficient method for solving large-scale linear equations. Some illustrative calculations are provided for molecules with up to 1484 Cartesian basis sets. The orbital domain sizes were found to be independent of the molecular size, and the present local MP2 method covered about 98%-99% of the correlation energy of the conventional canonical MP2 method.

True stabilization energies for the optimal planar hydrogen-bonded and stacked structures of guanine...cytosine, adenine...thymine, and their 9- and 1-methyl derivatives: complete basis set calculations at the MP2 and CCSD(T) levels and comparison with experiment

Planar H-bonded and stacked structures of guanine...cytosine (G.C), adenine...thymine (A...T), 9-methylguanine...1-methylcytosine (mG...mC), and 9-methyladenine...1-methylthymine (mA...mT) were optimized at the RI-MP2 level using the TZVPP ([5s3p2d1f/3s2p1d]) basis set. Planar H-bonded structures of G...C, mG...mC, and A...T correspond to the Watson-Crick (WC) arrangement, in contrast to mA...mT for which the Hoogsteen (H) structure is found. Stabilization energies for all structures were determined as the sum of the complete basis set limit of MP2 energies and a (DeltaE(CCSD(T)) - DeltaE(MP2)) correction term evaluated with the cc-pVDZ(0.25,0.15) basis set. The complete basis set limit of MP2 energies was determined by two-point extrapolation using the aug-cc-pVXZ basis sets for X = D and T and X = T and Q. This procedure is required since the convergency of the MP2 interaction energy for the present complexes is rather slow, and it is thus important to include the extrapolation to the complete basis set limit. For the MP2/aug-cc-pVQZ level of theory, stabilization energies for all complexes studied are already very close to the complete basis set limit. The much cheaper D-->T extrapolation provided a complete basis set limit close (by less than 0.7 kcal/mol) to the more accurate T-->Q term, and the D-->T extrapolation can be recommended for evaluation of complete basis set limits of more extended complexes (e.g. larger motifs of DNA). The convergency of the (DeltaE(CCSD(T)) - DeltaE(MP2)) term is known to be faster than that of the MP2 or CCSD(T) correlation energy itself, and the cc-pVDZ(0.25,0.15) basis set provides reasonable values for planar H-bonded as well as stacked structures. Inclusion of the CCSD(T) correction is essential for obtaining reliable relative values for planar H-bonding and stacking interactions; neglecting the CCSD(T) correction results in very considerable errors between 2.5 and 3.4 kcal/mol. Final stabilization energies (kcal/mol) for the base pairs studied are very substantial (A...T WC, 15.4; mA...mT H, 16.3; A...T stacked, 11.6; mA...mT stacked, 13.1; G...C WC, 28.8; mG...mC WC, 28.5; G...C stacked, 16.9; mG...mC stacked, 18.0), much larger than published previously. On the basis of comparison with experimental data, we conclude that our values represent the lower boundary of the true stabilization energies. On the basis of error analysis, we expect the present H-bonding energies to be fairly close to the true values, while stacked energies are still expected to be about 10% too low. The stacking energy for the mG...mC pair is considerably lower than the respective H-bonding energy, but it is larger than the mA...mT H-bonding energy. This conclusion could significantly change the present view on the importance of specific H-bonding interactions and nonspecific stacking interactions in nature, for instance, in DNA. Present stabilization energies for H-bonding and stacking energies represent the most accurate and reliable values and can be considered as new reference data.

Electronic transitions of guanine tautomers, their stacked dimers, trimers and sodium complexes

Planar and nonplanar geometries of the keto-N9H and keto-N7H tautomers of the guanine base of DNA as well as the hydrogen bonded complexes of these species with three water molecules each were optimized using the density functional theory at the B3LYP/6-31G** level. Geometries of the isolated bases were also optimized using the ab initio approach at the MP2/6-31G** level. The isolated keto-N9H and keto-N7H tautomers as well as their hydrogen bonded complexes with three water molecules each were solvated in bulk water employing the polarized continuum model (PCM) of the self-consistent reaction field theory (SCRF). Stacked dimers and trimers of both the tautomers of guanine were generated by placing the planar forms of the species at a fixed distance of 3.5 A from the neighboring one and rotating one molecule with respect to the other by 110 degrees for the keto-N9H form and 90 degrees for the keto-N7H form which corresponded to total energy minima at the B3LYP/6-31G** level. Geometry optimization for the cation of the monomer of guanine was performed at the same level of theory, and its solvation in bulk water was treated using the PCM model of the SCRF theory. The geometries of complexes of the two tautomers of guanine with a Na+ ion each were optimized at the B3LYP/6-31G** level, and the Na+ ion is predicted to bind with the keto-N9H tautomer preferentially. While the complex of the keto-N7H form of guanine with three water molecules in gas phase is slightly more stable than the corresponding complex of the keto-N9H form of guanine, the reverse is true in bulk water. Stacking interactions enhance the relative stability of the keto-N9H tautomer over that of the keto-N7H tautomer, suggesting that in bulk solutions, the former would be dominant. Electronic spectra of the isolated tautomers of guanine, those of their complexes with three water molecules each, the (keto-n9h and keto-n7h) cation of guanine, the complexes of the tautomers with a Na+ ion each, the stacked dimers and trimers of the two tautomers were calculated using configuration interaction involving single electron excitations (CIS). The relative absorption intensities of the two tautomers of guanine near 275 and 248 nm in the monomer, dimer, and trimer are predicated to be in the opposite order. Thus the absorption intensity oscillation observed using a guanine aqueous solution can be explained in terms of oscillation of relative populations of the two tautomers of the molecule. The 248 nm absorption peak would be appreciably red-shifted on formation of the cation of guanine. Binding of the Na+ ion with the two tautomers of guanine reduces intensities of their transitions appreciably and also it causes large red-shifts in the same.