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.     

Characterization of nucleobase-amino acid stacking interactions utilized by a DNA repair enzyme

The present work characterizes the gas-phase stacking interactions between four aromatic amino acid residues (histidine, phenylalanine, tyrosine, and tryptophan) and adenine or 3-methyladenine due to the proposed utilization of these interactions by enzymes that repair DNA alkylation damage. The MP2 potential energy surfaces of the stacked dimers are considered as a function of four variables (vertical displacement, angle of rotation, horizontal displacement, and tilt angle) using a variety of basis sets. It is found that the maximum stacking interaction energy decreases with the amino acid according to TRP > TYR approximately HIS > PHE for both nucleobases. However, the magnitude of the stacking interaction significantly increases upon alkylation (by 50-115%). Comparison of the stacking energies calculated using our surface scans to those estimated from experimental crystal structures indicates that the stacking interactions within the active site of 3-methyladenine DNA glycosylase can account for 65-75% of the maximum possible stacking interaction between the relevant molecules. The decrease in stacking in the crystal structure arises due to significant differences in the relative orientations of the nucleobase and amino acid. Nevertheless, alkylation is found to significantly increase the stacking energy when the crystal structure geometries are considered. Our calculations provide computational support for suggestions that alkylation enhances the stacking interactions within the active site of DNA repair enzymes, and they give a measure of the magnitude of this enhancement. Our results suggest that alkylation likely plays a more important role in substrate identification and removal than the nature of the aromatic amino acid that interacts with the substrate via stacking interactions.     

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.     

Insights into DNA binding of ruthenium arene complexes: role of hydrogen bonding and pi stacking

Density functional theory (DFT) methods are used to investigate the binding of ruthenium arene complexes, proposed as promising anticancer drugs, to isolated nucleobases. This shows a clear preference for binding at guanine over any other base and an approximately 100 kJ mol (-1) difference in binding between guanine and adenine in the gas phase, while binding to cytosine and inosine are intermediate in energy between these extremes. Solvation reduces binding energies and the discrimination between bases but maintains the overall pattern of binding. DFT and ab initio data on arene-base interactions in the absence of ruthenium show that stacking and hydrogen-bonding interactions play a significant role but cannot account for all of the energy difference between bases observed. Atoms-in-molecules analysis allows further decomposition of binding energies into contributions from covalent-binding, hydrogen-bonding, and pi-stacking interactions. Larger arenes undergo stabilizing stacking interactions, whereas N-H...X hydrogen bonding is independent of arene. Pairing of guanine to cytosine is affected by ruthenium complexation, with individual hydrogen-bonding energies being altered but the overall pairing energy remaining almost constant.     

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.     

Scope and limitations of the SCS-MP2 method for stacking and hydrogen bonding interactions

Fluorobenzenes are pi-acceptor synthons that form pi-stacked structures in molecular crystals as well as in artificial DNAs. We investigate the competition between hydrogen bonding and pi-stacking in dimers consisting of the nucleobase mimic 2-pyridone (2PY) and all fluorobenzenes from 1-fluorobenzene to hexafluorobenzene (n-FB, with n = 1-6). We contrast the results of high level ab initio calculations with those obtained using ultraviolet (UV) and infrared (IR) laser spectroscopy of isolated and supersonically cooled dimers. The 2PY.n-FB complexes with n = 1-5 prefer double hydrogen bonding over pi-stacking, as diagnosed from the UV absorption and IR laser depletion spectra, which both show features characteristic of doubly H-bonded complexes. The 2-pyridone.hexafluorobenzene dimer is the only pi-stacked dimer, exhibiting a homogeneously broadened UV spectrum and no IR bands characteristic for H-bonded species. MP2 (second-order M?ller-Plesset perturbation theory) calculations overestimate the pi-stacked dimer binding energies by about 10 kJ/mol and disagree with the experimental observations. In contrast, the MP2 treatment of the H-bonded dimers appears to be quite accurate. Grimme's spin-component-scaled MP2 approach (SCS-MP2) is an improvement over MP2 for the pi-stacked dimers, reducing the binding energy by approximately 10 kJ/mol. When applied to explicitly correlated MP2 theory (SCS-MP2-R12 approach), agreement with the corresponding coupled-cluster binding energies [at the CCSD(T) level] is very good for the pi-stacked dimers, within +/- 1 kJ/mol for the 2PY complexes with 1-fluorobenzene, 1,2-difluorobenzene, 1,2,4,5-tetrafluorobenzene, pentafluorobenzene and hexafluorobenzene. Unfortunately, the SCS-MP2 approach also reduces the binding energy of the H-bonded species, leading to disagreement with both coupled-cluster theory and experiment. The SCS-MP2-R12 binding energies follow the SCS-MP2 binding energies closely, being about 0.5 and 0.7 kJ/mol larger for the H-bonded and pi-stacked forms, respectively, in an augmented correlation-consistent polarized valence quadruple-zeta basis. It seems that the SCS-MP2 and SCS-MP2-R12 methods cannot provide sufficient accuracy to replace the CCSD(T) method for intermolecular interactions where H-bonding and pi-stacking are competitive.     

The solvation, partitioning, hydrogen bonding, and dimerization of nucleotide bases: a multifaceted challenge for quantum chemistry

We present M06-2X density functional calculations of the chloroform/water partition coefficients of cytosine, thymine, uracil, adenine, and guanine and calculations of the free energies of association of selected unsubstituted and alkylated nucleotide base pairs in chloroform and water. Both hydrogen bonding and π-π stacking interactions are considered. Solvation effects are treated using the continuum solvent models SM8, SM8AD, and SMD, including geometry optimization in solution. Comparison of theoretical results with available experimental data indicates that all three of these solvation models predict the chloroform-water partition coefficients for the studied nucleobases qualitatively well, with mean unsigned errors in the range of 0.4-1.3 log units. All three models correctly predict the preference for hydrogen bonding over stacking for nucleobase pairs solvated in chloroform, and SM8, SM8AD, and SMD show similar accuracy in predicting the corresponding free energies of association. The agreement between theory and experiment for the association free energies of the dimers in water is more difficult to assess, as the relevant experimental data are indirect. Theory predicts that the stacking interaction of nucleobases in water is more favorable than hydrogen bonding for only two out of three tested hetero-dimers.     

Cross-modulation of the pKa of nucleobases in a single-stranded hexameric-RNA due to tandem electrostatic nearest-neighbor interactions

The pH titration studies (pH 6.7-12.1) in a series of dimeric, trimeric, tetrameric, pentameric, and hexameric oligo-RNA molecules [GpA (2a), GpC (3a), GpApC (5), GpA(1)pA(2)pC (6), GpA(1)pA(2)pA(3)pC (7), GpA(1)pA(2)pA(3)pA(4)pC (8)] have shown that the pK(a) of N(1)-H of 9-guaninyl could be measured not only from its own deltaH8G, but also from the aromatic marker protons of other constituent nucleobases. The relative chemical shift differences [Deltadelta((N)(-)(D))] between the protons in various nucleotide residues in the oligo-RNAs at the neutral (N) and deprotonated (D) states of the guanine moiety show that the generation of the 5'-(9-guanylate ion) in oligo-RNAs 2-8 reduces the stability of the stacked helical RNA conformation owing to the destabilizing anion(G(-))-pi/dipole(Im(delta)(-)) interaction. This destabilizing effect in the deprotonated RNA is, however, opposed by the electrostatically attractive atom-pisigma (major) as well as the anion(G(-))-pi/dipole(Py(delta)(+)) (minor) interactions. Our studies have demonstrated that the electrostatically repulsive anion(G(-))-pi/dipole(Im(delta)(-)) interaction propagates from the first to the third nucleobase quite strongly in the oligo-RNAs 6-8, causing destacking of the helix, and then its effect is gradually reduced, although it is clearly NMR detectable along the RNA chain. Thus, such specific generation of a charge at a single nucleobase moiety allows us to explore the relative strength of stacking within a single-stranded helix. The pK(a) of 5'-Gp residue from its own deltaH8G in the hexameric RNA 8 is found to be 9.76 +/- 0.01; it, however, varies from 9.65 +/- 0.01 to 10.5 +/- 0.07 along the RNA chain as measured from the other marker protons (H2, H8, H5, and H6) of 9-adeninyl and 1-cytosinyl residues. This nucleobase-dependent modulation of pK(a)s (DeltapK(a) +/- 0.9) of 9-guaninyl obtained from other nucleobases in the hexameric RNA 8 represents a difference of ca. 5.1 kJ mol(-)(1), which has been attributed to the variable strength of electrostatic interactions between the electron densities of the involved atoms in the offset stacked nucleobases as well as with that of the phosphates. The chemical implication of this variable pK(a) for guanin-9-yl deprotonation as obtained from all other marker protons of each nucleotide residue within a ssRNA molecule is that it enables us to experimentally understand the variation of the electronic microenvironment around each constituent nucleobase along the RNA chain in a stepwise manner with very high accuracy without having to make any assumption. This means that the pseudoaromaticity of neighboring 9-adeninyl and next-neighbor nucleobases within a polyanionic sugar-phosphate backbone of a ssRNA can vary from one case to another due to cross-modulation of an electronically coupled pi system by a neighboring nucleobase. This modulation may depend on the sequence context, spatial proximity of the negatively charged phosphates, as well as whether the offset stacking is ON or OFF. The net outcome of this electrostatic interaction between the neighbors is creation of new sequence-dependent hybrid nucleobases in an oligo- or polynucleotide whose properties are unlike the monomeric counterpart, which may have considerable biological implications.     

yDNA versus xDNA pyrimidine nucleobases: computational evidence for dependence of duplex stability on spacer location

The structural and binding properties of the natural and x- and y-pyrimidines were compared using computational methods. Our calculations show that although the x-pyrimidines favor different orientations about the glycosidic bond compared to the natural pyrimidines, which could have implications for the formation and resulting stability of xDNA duplexes and jeopardize the selectivity of expanded nucleobases, y-pyrimidines have rotational profiles more similar to the natural bases. Increasing the pyrimidine size using a benzene spacer leads to relatively minor changes in the hydrogen-bond strength of isolated Watson-Crick base pairs. However, differences in the anomeric carbon distances in pairs composed of x- or y-pyrimidines suggest yDNA may yield a more optimal expanded structure. By stacking two monomers via their centers of mass, we find that the expanded nucleobases stack much stronger than the natural bases. Additionally, although replacing xT by yT changes the stacking energy by less than 5 kJ mol (-1), replacing xC by yC significantly strengthens complexes with the natural nucleobases (by up to 30%). Calculations on larger duplex models composed of four nucleobases reveal that x- and y-pyrimidines can increase duplex stability of natural helices by strengthening both the intra and interstrand stacking interactions. Furthermore, when the total stability (sum of all hydrogen-bonding and (intrastrand and interstrand) stacking interactions) of the larger models is considered, y-pyrimidines lead to more stable complexes than x-pyrimidines for all but three duplex sequences. Thus, through analysis of a variety of properties, our calculations suggest that the location of the benzene spacer affects the properties of expanded nucleobases and the stability of expanded duplexes, and therefore should be carefully considered when designing future expanded analogues.     

Toward true DNA base-stacking energies: MP2, CCSD(T), and complete basis set calculations

Stacking energies in low-energy geometries of pyrimidine, uracil, cytosine, and guanine homodimers were determined by the MP2 and CCSD(T) calculations utilizing a wide range of split-valence, correlation-consistent, and bond-functions basis sets. Complete basis set MP2 (CBS MP2) stacking energies extrapolated using aug-cc-pVXZ (X = D, T, and for pyrimidine dimer Q) basis sets equal to -5.3, -12.3, and -11.2 kcal/mol for the first three dimers, respectively. Higher-order correlation corrections estimated as the difference between MP2 and CCSD(T) stacking energies amount to 2.0, 0.7, and 0.9 kcal/mol and lead to final estimates of the genuine stacking energies for the three dimers of -3.4, -11.6, and -10.4 kcal/mol. The CBS MP2 stacking-energy estimate for guanine dimer (-14.8 kcal/mol) was based on the 6-31G(0.25) and aug-cc-pVDZ calculations. This simplified extrapolation can be routinely used with a meaningful accuracy around 1 kcal/mol for large aromatic stacking clusters. The final estimate of the guanine stacking energy after the CCSD(T) correction amounts to -12.9 kcal/mol. The MP2/6-31G(0.25) method previously used as the standard level to calculate aromatic stacking in hundreds of geometries of nucleobase dimers systematically underestimates the base stacking by ca. 1.0-2.5 kcal/mol per stacked dimer, covering 75-90% of the intermolecular correlation stabilization. We suggest that this correction is to be considered in calibration of force fields and other cheaper computational methods. The quality of the MP2/6-31G(0.25) predictions is nevertheless considerably better than suggested on the basis of monomer polarizability calculations. Fast and very accurate estimates of the MP2 aromatic stacking energies can be achieved using the RI-MP2 method. The CBS MP2 calculations and the CCSD(T) correction, when taken together, bring only marginal changes to the relative stability of H-bonded and stacked base pairs, with a slight shift of ca. 1 kcal/mol in favor of H-bonding. We suggest that the present values are very close to ultimate predictions of the strength of aromatic base stacking of DNA and RNA bases.     

X-ray diffraction and theoretical studies for the quantitative assessment of intermolecular arene-perfluoroarene stacking interactions

The arene-perfluoroarene stacking interaction was studied by experimental and theoretical methods. A series of compounds with different possibilities for formation of this recognition motif in the solid state were synthesized, and their crystal structures determined by single-crystal X-ray diffraction. The crystal packing of these compounds, as well as the packing of related compounds retrieved from crystallographic databases, were analyzed with quantitative crystal potentials: total lattice energies and the cohesive energies of closest molecular pairs in the crystals were calculated. The arene-perfluoroarene recognition motif emerges as a dominant interaction in the non-hydrogen-bonding compounds studied here, to the point that asymmetric dimers formed over the stacking motif carry over to asymmetric units made of two molecules in the crystal both for pure compounds and for molecular complexes; however, inter-ring distances and angles range from 3.70 to 4.85 A and from 5 to 21 degrees , respectively. Pixel energy partitioning reveals that whenever aromatic rings stack, the largest cohesive energy contribution comes from dispersion, which roughly amounts to 20 kJ mol(-1) per phenyl ring, while the coulombic term is minor but significant enough to make a difference between the arene-arene or perfluoroarene-perfluoroarene interactions on the one hand, and arene-perfluoroarene interactions on the other, whereby the latter are favored by about 10 kJ mol(-1) per phenyl ring. No evidence of special interaction which can be attributed to HF confrontation was recognizable.     

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.     

Isoelectronic homologues and isomers: tropolone, 5-azatropolone, 1-H-azepine-4,5-dione, saddle points, and ions

Computational studies of 12 64-electron homologues and isomers of tropolone in the S(0) electronic ground state are reported. Three minimum-energy structures, tropolone (Tp), 5-azatropolone (5Azt), and 5-H-5-azatropolonium (5AztH(+)), have an internal H-bond and planar C(s)) geometry, and three, tropolonate (TpO(-)), 5-azatropolonate (5AzO(-)), and 1-H-azepine-4,5-dione (45Di), lack the H-bond and have twisted C(2) geometry. All 6 substances have an equal double-minimum potential energy surface and a saddle point with planar C(2)(v) geometry. The energy for the gas-phase isomerization reaction 45Di --> 5Azt is near +4 kJ mol(-1) at the MP4(SDQ)/6-311++G(df,pd)//MP2/6-311++G(df,pd) (energy//geometry) theoretical level and around -20 kJ mol(-1) at lower theoretical levels. The dipole moments computed for 45Di and 5Azt are 9.6 and 2.1 D, respectively, and this large difference contributes to MO-computed free energies of solvation that strongly favor--as experimentally observed--45Di over 5Azt in chloroform solvent. The MO-computed energy for the gas-phase protonation reaction 45Di + H(+) --> 5AztH(+) is -956.4 kJ mol(-1), leading to 926.8 kJ mol(-1) as the estimated proton affinity for 45Di at 298 K and 1 atm. The intramolecular dynamical properties predicted for 5Azt and 5AztH(+) parallel those observed for tropolone. They are therefore expected to exhibit spectral tunneling doublets. Once they are synthesized, they should contribute importantly to the understanding of multidimensional intramolecular H transfer and dynamical coupling processes.     

Environmental effects on the enhancement in natural and damaged DNA nucleobase acidity because of discrete hydrogen-bonding interactions

The present study uses density functional theory to carefully consider the effects of the environment on the enhancement in (natural and damaged) DNA nucleobase acidities because of multiple hydrogen-bonding interactions. Although interactions with one small molecule can increase the acidity of the nucleobases by up to 60 kJ mol-1 in the gas phase, the maximum increase in enzymatic-like environments is expected to be approximately 40 kJ mol-1, which reduces to approximately 30 kJ mol-1 in water. Furthermore, the calculated (simultaneous) effects of two, three, or four molecules are increasingly less than the sum of the individual (additive) effects with an increase in the number and acidity of the small molecules bound or the dielectric constant of the solvent. Regardless of these trends, our calculations reveal that additional hydrogen-bonding interactions will have a significant effect on nucleobase acidity in a variety of environments, where the exact magnitude of the effect depends on the properties of the small molecule bound, the nucleobase binding site, and the solvent. The maximum increase in nucleobase acidity because of interactions with up to four small molecules is approximately 80 kJ mol-1 in enzymatic-like environments (or 65 kJ mol-1 in water). These results suggest that hydrogen-bonding interactions likely play an important role in many biological processes by changing the physical and chemical properties of the nucleobases.     

Theoretical and vibrational study of the conformation of 2-methoxy- 1,2-diphenylethanone.

The conformation and vibrational properties of 2-methoxy-1,2-diphenylethanone (MDPE) are investigated in the gas phase and in organic solvents. Ab initio calculations carried out at the B3LYP/6-31G(d) level demonstrate that three stable conformers having cisoid, skewed and transoid structures are present in the gas phase. In the gas phase, the conformers are separated by a low energy barrier and their relative energies do not differ by more than 7.2 kJ mol (-1) Like in crystalline MDPE 'Acta Crystallogr. Sect. C 44 (1988) 894', weak CH...O hydrogen bonds are present in the cisoid conformation. The IR and Raman spectra of solid MDPE are discussed. Several vibrational modes are split in organic solvents. A comparison between the theoretical data and the experimental dipole moments indicates that two conformers are present in solution, the population of the cisoid form increasing with the permittivity of the medium.     

How flexible are fleximer nucleobases? A computational study

Density functional theory was used to study the potential energy surface for rotation about the carbon-carbon bonds in a variety of guanosine, adenosine, and inosine fleximers, which are modified purines with the imidazole and pyrimidine rings separated by a single carbon-carbon bond. Various connectivities between C4 or C5 of the imidazole ring and C5' or C6' of the pyrimidine ring were considered. Calculations on fleximer nucleobases in the absence of the ribose moiety suggest that a planar relative arrangement of the imidazole and pyrimidine rings is favored, and that all fleximers are indeed very flexible with regards to rotation about the carbon-carbon bond, where calculated barriers are generally less than 40 kJ mol(-1). Furthermore, calculated binding energies of fleximer-pyrimidine pairs indicate that the hydrogen-bonding properties of these modified nucleobases mimic those of the corresponding natural purine. Inclusion of the sugar moiety often leads to a favored nonplanar orientation of the two rings, and either a reduction in the rotational barrier height or small changes in the rotational surface depending on the connectivity and nucleobase considered. It is concluded that several connectivities may have favorable properties for biochemical applications where flexible nucleobases would be beneficial.     

Local Structures and Chemical Properties of Deprotonated Arginine

The potential energy surface of gaseous deprotonated arginine has been systematically in- vestigated by first principles calculations. At the B3LYP/6-31G(d) level, apart from the identification of several stable local structures, a new global minimum is located which is about 6.56 kJ/mol more stable than what has been reported. The deprotonated arginine molecule has two distinct forms with the deprotonation at the carboxylate group (COO^-). These two forms are bridged by a very high energy barrier and possess very different IR spectral profiles. Our calculated proton dissociation energy and gas-phase acidity of argi-nine molecule are found to be in good agreement with the corresponding experimental results. The predicted geometries, dipole moments, rotational constants, vertical ionization energies and IR spectra of low energy conformers will be useful for future experimental measurements.     

How Can Rotaxanes Be Modified by Varying Functional Groups at the Axle?—A Combined Theoretical and Experimental Analysis of Thermochemistry and Electronic Effects

We present theoretically as well as experimentally determined thermochemical data of the non-covalent interactions in different axle-substituted pseudorotaxanes. The overall interaction energy lies in the region of 35 kJ mol(-1), independent of the substitution pattern at the axle. Because rearrangement energies of 7 and 3 kJ mol(-1) are required for wheel and axle, respectively, the sum of the net interactions of individual non-covalent bonds must exceed 10 kJ mol(-1) to achieve a successful host-guest interaction. The geometrical analysis shows three hydrogen bonds, and the close inspection of the individual dipole moments as well as the individual hydrogen bonds reveals trends according to the different functional groups at the axle. The individual trends for the different hydrogen bonds almost lead to a cancellation of the substitution effects. From solvent-effect considerations it can be predicted that the pseudorotaxane is stable in CHCl(3) and CH(2)Cl(2), whereas it would dethread in water. Comparing experimentally and theoretically calculated Gibbs free enthalpies, we find reasonable agreement if an exchange reaction of one solvent molecule instead of the direct formation reaction is considered.     

Cooperative and diminutive unusual weak bonding in F3CX···HMgH···Y and F3CX···Y···HMgH trimers (X = Cl, Br; Y = HCN, and HNC)

MP2 calculations with cc-pVTZ basis set were used to analyze intermolecular interactions in F(3)CX···HMgH···Y and F(3)CX···Y···HMgH triads (X = Cl, Br; Y = HCN, and HNC) which are connecting with three kinds of unusual weak interactions, namely halogen-hydride, dihydrogen, and σ-hole. To understand the properties of the systems better, the corresponding dyads are also studied. Molecular geometries, binding energies, and infrared spectra of monomers, dyads, and triads were investigated at the MP2/cc-pVTZ computational level. Particular attention is given to parameters such as cooperative energies, cooperative dipole moments, and many-body interaction energies. Those complexes with simultaneous presence of a σ-hole bond and a dihydrogen bond show cooperativity energy ranging between -1.02 and -2.31 kJ mol(-1), whereas those with a halogen-hydride bond and a dihydrogen bond are diminutive, with this energetic effect between 0.1 and 0.63 kJ mol(-1). The electronic properties of the complexes have been analyzed using the molecular electrostatic potential (MEP), the electron density shift maps, and the parameters derived from the atoms in molecules (AIM) methodology.