Non-Watson-Crick base pairing in RNA. quantum chemical analysis of the cis Watson-Crick/sugar edge base pair family

Large RNA molecules exhibit an astonishing variability of base-pairing patterns, while many of the RNA base-pairing families have no counterparts in DNA. The cis Watson-Crick/sugar edge (cis WC/SE) RNA base pairing is investigated by ab initio quantum chemical calculations. A detailed structural and energetic characterization of all 13 crystallographically detected members of this family is provided by means of B3LYP/6-31G and RIMP2/aug-cc-pVDZ calculations. Further, a prediction is made for the remaining 3 cis WC/SE base pairs which are yet to be seen in the experiments. The interaction energy calculations point at the key role of the 2'-OH group in stabilizing the sugar-base contact and predict all 16 cis WC/SE base-pairing patterns to be nearly isoenergetic. The perfect correlation of the main geometrical parameters in the gas-phase optimized and X-ray structures shows that the principle of isosteric substitutions in RNA is rooted from the intrinsic structural similarity of the isolated base pairs. The present quantum chemical calculations for the first time analyze base pairs involving the ribose 2'-OH group and unambiguously correlate the structural information known from experiments with the energetics of interactions. The calculations further show that the relative importance and absolute value of the dispersion energy in the cis WC/SE base pairs are enhanced compared to the standard base pairs. This may by an important factor contributing to the strength of such interactions when RNA folds in its polar environment. The calculations further demonstrate that the Cornell et al. force field commonly used in molecular modeling and simulations provides satisfactory performance for this type of RNA interactions.     

Sugar edge/sugar edge base pairs in RNA: stabilities and structures from quantum chemical calculations

Cis and trans sugar edge/sugar edge (SE/SE) binding patterns are essential building units of RNAs. For example, SE/SE interactions form the A-minor motifs, the most important tertiary interaction type in functional RNAs. This study provides an in-depth structure and stability analysis for these two base pair families. Gas-phase-optimized geometries are reported for 12 cis and 7 trans SE/SE base pairs and contrasted to their X-ray counterparts. Interaction energies are computed at the RIMP2 level of theory using the density-functional-theory-optimized geometries. There is a good overall agreement between the optimized and X-ray geometries of the cis SE/SE base pairs. In contrast, only three of the seven trans SE/SE binding patterns could be optimized without a significant distortion of the X-ray geometry. Note, however, that many SE/SE base pairs participate in broader networks of interactions; thus it is not surprising to see some of them to deviate from the X-ray geometry in a complete isolation. Computed interaction energies reveal that all 12 known cis SE/SE binding patterns are very stable. Among the trans SE/SE binding patterns, only the rG/rG, rG/rC, and rA/rG base pairs are sufficiently stable in the crystal geometry. Prediction has been made for some structures not yet detected by crystallography, namely, cis rC/rC, rG/rC, rG/rU, and rU/rU and trans rG/rA base pairs. Interestingly, the new cis SE/SE binding patterns are not necessarily isosteric with the remaining 12 members of this family. The trans rG/rA base pair represents a viable option for base pairing in RNA to be identified by future X-ray studies. In a complete lack of structural information, prediction of other unknown members of the trans SE/SE family was not attempted. Analysis of the interaction energies shows a very large electron correlation component of the interaction energy, pointing at the elevated role of dispersion energy as compared to other types of base pairs. This likely is profitable for stabilization of SE/SE binding patterns in polar environments and could be one of the reasons why the A-minor motif is the leading type of tertiary interactions in RNAs.     

Conformational specificity of non-canonical base pairs and higher order structures in nucleic acids: crystal structure database analysis

Non-canonical base pairs contribute immensely to the structural and functional variability of RNA, which calls for a detailed characterization of their spatial conformation. Intra-base pair parameters, namely propeller, buckle, open-angle, stagger, shear and stretch describe structure of base pairs indicating planarity and proximity of association between the two bases. In order to study the conformational specificities of non-canonical base pairs occurring in RNA crystal structures, we have upgraded NUPARM software to calculate these intra-base pair parameters using a new base pairing edge specific axis system. Analysis of base pairs and base triples with the new edge specific axis system indicate the presence of specific structural signatures for different classes of non-canonical pairs and triples. Differentiating features could be identified for pairs in cis or trans orientation, as well as those involving sugar edges or C-H-mediated hydrogen bonds. It was seen that propeller for all types of base pairs in cis orientation are generally negative, while those for trans base pairs do not have any preference. Formation of a base triple is seen to reduce propeller of the associated base pair along with reduction of overall flexibility of the pairs. We noticed that base pairs involving sugar edge are generally more non-planar, with large propeller or buckle values, presumably to avoid steric clash between the bulky sugar moieties. These specific conformational signatures often provide an insight into their role in the structural and functional context of RNA.     

Theoretical study of the scalar coupling constants across the noncovalent contacts in RNA base pairs: the cis- and trans-watson-crick/sugar edge base pair family

The structure and function of RNA molecules are substantially affected by non-Watson-Crick base pairs actively utilizing the 2'-hydroxyl group of ribose. Here we correlate scalar coupling constants across the noncovalent contacts calculated for the cis- and trans-WC/SE (Watson-Crick/sugar edge) RNA base pairs with the geometry of base to base and sugar to base hydrogen bond(s). 23 RNA base pairs from the 32 investigated were found in RNA crystal structures, and the calculated scalar couplings are therefore experimentally relevant with regard to the binding patterns occurring in this class of RNA base pairs. The intermolecular scalar couplings 1hJ(N,H), 2hJ(N,N), 2hJ(C,H), and 3hJ(C,N) were calculated for the N-H...N and N-H...O=C base to base contacts and various noncovalent links between the sugar hydroxyl and RNA base. Also, the intramolecular 1J(N,H) and 2J(C,H) couplings were calculated for the amino or imino group of RNA base and the ribose 2'-hydroxyl group involved in the noncovalent interactions. The calculated scalar couplings have implications for validation of local geometry, show specificity for the amino and imino groups of RNA base involved in the linkage, and can be used for discrimination between the cis- and trans-WC/SE base pairs. The RNA base pairs within an isosteric subclass of the WC/SE binding patterns can be further sorted according to the scalar couplings calculated across different local noncovalent contacts. The effect of explicit water inserted in the RNA base pairs on the magnitude of the scalar couplings was calculated, and the data for discrimination between the water-inserted and direct RNA base pairs are presented. The calculated NMR data are significant for structural interpretation of the scalar couplings in the noncanonical RNA base pairs.     

Cis Watson-Crick/Watson-Crick Base Pairs in DNA: A Density Functional Theory(DFT) Study

The four nucleic acid DNA bases(adenine, thymine, guanine, cytosine) and ten cis Watson-Crick/Watson-Crick(cis WC/WC) DNA base pairs were investigated by density functional theory(DFT) quantum chemical calculations. Geometry optimizations were carried out on the four bases and ten base pairs at the B3LYP level with 6-31G~(**) basis set. All the optimizations were performed within Cs symmetry. The optimum structures for the four bases and seven cis WC/WC base pairs were obtained, and Natural Bond Orbital analysis(NBO) was based on these structures. The possibilities of matches between any two of the four bases through their Watson-Crick(WC) edges were discussed. The structures of seven cis WC/WC base pairs change to a certain extent relative to these of the four bases due to the formation of hydrogen bonds. These base pairs existing in DNA have an important influence on the structural stability of the double helix. The analysis of the electronic structures and molecular orbitals for seven cis WC/WC base pairs can provide significant information about the relationship between charge transfer along the hydrogen bond and the Frontier orbitals of these base pairs.     

Leading RNA tertiary interactions: structures, energies, and water insertion of A-minor and P-interactions. A quantum chemical view

Complex molecular shapes of ribosomal RNA molecules are stabilized by recurrent types of tertiary interactions involving highly specific and conserved non-Watson-Crick base pairs, triplets, and quartets. We analyzed the intrinsic structure and stability of the P-motif and the four basic types of A-minor interactions (types I, II, III, and 0), which represent the most prominent RNA tertiary interaction patterns refined in the course of evolution. In the studied interactions, the electron correlation component of the stabilization usually exceeds the Hartree-Fock (HF) term, leading to a strikingly different balance of forces as compared to standard base pairing stabilized primarily by the HF term. In other words, the A-minor and P-interactions are considerably more influenced by the dispersion energy as compared to canonical base pairs, which makes them particularly suitable to zip the folded RNA structures that are substantially hydrated even in their interior. Continuum solvent COSMO calculations confirm that the stability of the canonical GC base pair is affected (reduced) by the continuum solvent screening considerably more than the stability of the A-minor interaction. Among the studied systems, the strong A-minor II and weak A-minor III interactions require water molecules to stabilize the experimental geometry. Gas-phase optimization of the canonical A-minor II A/CG triplet without water results in a geometry that is clearly inconsistent with the RNA structure. The gas-phase structure of the P-interaction and the most stable A-minor I interaction nicely agrees with the geometries occurring in the ribosome. A-minor I can also adopt an alternative water-mediated substate rather often observed in X-ray and molecular dynamics studies. The A-minor I water bridge, however, does not appear to stabilize the tertiary contact, and its role is to provide structural flexibility to this binding pattern within the context of the RNA structure. Interestingly, the insertion of a polar water molecule in the A-minor I A/CG tertiary contact occurring in the A/C tertiary pair is stabilized primarily by the HF (electrostatic) interaction energy, while the dispersion-controlled A/G contact remains firmly bound. Thus, the intrinsic balance of forces as revealed by quantum mechanics (QM) calculations nicely correlates with many behavioral aspects of the studied interactions inside RNA. The comparison of interaction energies computed using quantum chemistry and an AMBER force field reveals that common molecular mechanics calculations perform rather well, except that the strength of the P-interaction is modestly overestimated. We also briefly discuss the non-negligible methodological differences when evaluating simple base-base nucleic acids base pairs and the complex RNA tertiary base pairing patterns using QM procedures.     

Hydrogen-bonded nucleic acid base pairs containing unusual base tautomers: complete basis set calculations at the MP2 and CCSD(T) levels

The total interaction energies of altogether 15 hydrogen-bonded nucleic acid base pairs containing unusual base tautomers were calculated. The geometry properties of all selected adenine-thymine and guanine-cytosine hydrogen-bonded base pairs enable their incorporation into DNA. Unusual base pairing patterns were compared with Watson-Crick H-bonded structures of the adenine-thymine and guanine-cytosine pairs. The complete basis set (CBS) limit of the MP2 interaction energy and the CCSD(T) correction term, determined as the difference between the CCSD(T) and MP2 interaction energies, was evaluated. Extrapolation to the MP2 CBS limit was done using the aug-cc-pVDZ and aug-cc-pVTZ results, and the CCSD(T) correction term was determined with the 6-31G*(0.25) basis set. Final interaction energies were corrected while taking into account both tautomeric penalization determined at the CBS level and solvation/desolvation free energies. The situation for the adenine-thymine pairs is straightforward, and tautomeric pairs are significantly less stable than the Watson-Crick pair consisting of the canonical forms. In the case of the guanine-cytosine pair, the Watson-Crick structure made by canonical forms is again the most stable. The other two structures are, however, energetically rather similar (by 5 and 6 kcal/mol), which provides a very small but non-negligible chance of detecting these structures in the DNA double helix (1:5000). Due to the fact that DNA bases and base pairs incorporated into DNA are solvated less favorably than in isolated systems, this probability represents the very upper limit. The results clearly show how precisely the canonical building blocks of DNA molecules were chosen and how well their stability is maintained.     

Base pairing in RNA structures: A computational analysis of structural aspects and interaction energies

The base pairing patterns in RNA structures are more versatile and completely different as compared to DNA. We present here results of ab-initio studies of structures and interaction energies of eight selected RNA base pairs reported in literature. Interaction energies, including BSSE correction, of hydrogen added crystal geometries of base pairs have been calculated at the HF/6-31G** level. The structures and interaction energies of the base pairs in the crystal geometry are compared with those obtained after optimization of the base pairs. We find that the base pairs become more planar on full optimization. No change in the hydrogen bonding pattern is seen. It is expected that the inclusion of appropriate considerations of many of these aspects of RNA base pairing would significantly improve the accuracy of RNA secondary structure prediction.     

An ab initio MO study on the structures and electronic states of hydrogen-bonded O3- HF and SO2-HF complexes

Ab initio molecular orbital (MO) calculations have been carried out for base-hydrogen fluoride (HF) complexes (base = O3 and SO2) in order to elucidate the structures and energetics of the complexes. The ab initio calculations were performed up to the QCISD(T)/6-311++G(d,p) level of theory. In both complexes, hydrogen-bonded structures where the hydrogen of HF orients toward one of the oxygen atoms of bases were obtained as stable forms. The calculations showed that cis and trans isomers exist in both complexes. All calculations for the SO2-HF complex indicated that the cis form is more stable in energy than the trans form. On the other hand, in O3-HF complexes, the stable structures are changed by the ab initio levels of theory used, and the energies of the cis and trans forms are close to each other. From the most sophisticated calculations (QCISD(T)/6-311++G(d,p)//QCISD/6-311+G(d) level), it was predicted that the complex formation energies for cis SO2-HF, trans SO2-HF, cis O3-HF, and trans O3-HF are 6.1, 5.7, 3.4, and 3.6 kcal/mol, respectively, indicating that the binding energy of HF to SO2 is larger than that of O3. The harmonic vibrational frequencies calculated for cis O3-HF and cis SO2-HF complexes were in good agreement with the experimental values measured by Andrews et al. Also, the calculated rotation constants for cis SO2-HF agreed with the experiment.     

Ferrocene-linked thymine/uracil conjugates: base pairing directed self-assembly and supramolecular packing

Ferrocene-linked bis(nucleobase) (1a-c) and chimeric nucleobase (1d) conjugates have been synthesized from mono- and bis(hydroxybutyl)ferrocene 6 via Mitsunobu reaction as the key step. X-ray crystallographic studies of ferrocene bis(nucleobase) conjugates reveal two-dimensional supramolecular organizations of backbones through self-assembled Watson-Crick and reverse Watson-Crick type pairs. Ferrocene-bis(thymine) conjugate self-assembles by reverse Watson-Crick pairing, while the corresponding bis(uracil) conjugate self-assembles by alternating WC and reverse WC type pairing. Such continuous assemblies are not seen in monosubstituted ferrocene nucleobase conjugates which form only planar sheets. The results are interesting from the point of understanding and engineering supramolecular assemblies through rational design of base pairing patterns.     

A theoretical study of the rotational isomers of nitrosomethanol by semiempirical (AM1) and ab initio methods

The conformational potential energy surface as a function of the two internal torsion angles in C-nitrosomethanol has been obtained using the semiempirical AM1 method. Optimized geometries are reported for the local minima on this surface and also for the corresponding points on the HF/6-31G, 6-31G*, and 6-31G** surfaces. All methods predict cis and trans minima which occur in degenerate pairs, each pair being connected by a transition state of C_s symmetry. The AM1 structures are found to compare well with the corresponding ab initio structures. Ab initio HF/6-31G and HF/6-31G* harmonic vibrational frequencies are reported for the cis and trans forms of nitrosomethanol. When scaled appropriately the calculated frequencies are found to compare well with experimental frequencies. The ab initio calculations predict the energy barrier for cis → trans isomerization to be between 5.8 and 6.5 kcal/mol with the trans → cis isomerization barrier lying between 2.3 and 6.5 kcal/mol. The corresponding AM1 energy barriers are around 1 kcal/mol lower in energy. The ab initio calculations predict the barrier to conversion between the two cis rotamers to be very small with the AM1 value being around 1 kcal/mol. Both AM1 and ab initio calculations predict interconversion between trans rotamers to require between 1.2 and 1.4 kcal/mol.     

Structure, stability, and dynamics of canonical and noncanonical base pairs: quantum chemical studies

The importance of non-Watson-Crick base pairs in the three-dimensional structure of RNA is now well established. The structure and stability of these noncanonical base pairs are, however, poorly understood. We have attempted to understand structural features of 33 frequently occurring base pairs using density functional theory. These are of three types, namely (i) those stabilized by two or more polar hydrogen bonds between the bases, (ii) those having one polar and another C-H...O/N type interactions, and (iii) those having one H-bond between the bases and another involving one of the sugars linked to the bases. We found that the base pairs having two polar H-bonds are very stable as compared to those having one C-H...O/N interaction. Our quantitatively analysis of structures of these optimized base pairs indicates that they possess a different amount of nonplanarity with large propeller or buckle values as also observed in the crystal structures. We further found that geometry optimization does not modify the hydrogen-bonding pattern, as values of shear and open angle of the base pairs remain conserved. The structures of initial crystal geometry and final optimized geometry of some base pairs having only one polar H-bond and a C-H...O/N interaction, however, are significantly different, indicating the weak nature of the nonpolar interaction. The base pair flexibility, as measured from normal-mode analysis, in terms of the intrinsic standard deviations of the base pair structural parameters are in conformity with those calculated from RNA crystal structures. We also noticed that deformation of a base pair along the stretch direction is impossible for all of the base pairs, and movements of the base pairs along shear and open are also quite restricted. The base pair opening mode through alteration of propeller or buckle is considerably less restricted for most of the base pairs.     

Cyclohexyl "base pairs" stabilize duplexes and intensify pyrene fluorescence by shielding it from natural base pairs

In this study, we investigated the stability and structure of artificial base pairs that contain cyclohexyl rings. The introduction of a single pair of isopropylcyclohexanes into the middle of DNA slightly destabilized the duplex. Interestingly, as the number of the "base pairs" increased, the duplex was remarkably stabilized. A duplex with six base pairs was even more stable than one containing six A-T pairs. Thermodynamic analysis revealed that changes in entropy and not enthalpy contributed to duplex stability, demonstrating that hydrophobic interactions between isopropyl groups facilitated the base pairing, and thus stabilized the duplex. NOESY of a duplex containing an isopropylcyclohexane-methylcyclohexane pair unambiguously demonstrated its "pairing" in the duplex because distinct NOEs between the protons of cyclohexyl moieties and imino protons of both of the neighboring natural base pairs were observed. CD spectra of duplexes tethering cyclohexyl moieties also showed a positive-negative couplet that is characteristic of the B-form DNA duplex. Taken together, these results showed that cyclohexyl moieties formed base pairs in the DNA duplex without severely disturbing the helical structure of natural DNA. Next, we introduced cyclohexyl base pairs between pyrene and nucleobases as an "insulator" that suppresses electron transfer between them. We found a massive increase in the quantum yield of pyrene due to the efficient shielding of pyrene from nucleobases. The cyclohexyl base pairs reported here have the potential to prepare highly fluorescent labeling agents by multiplying fluorophores and insulators alternately into DNA duplexes.     

Rational design of hetero-ring-expanded guanine analogs with enhanced properties for modified DNA building blocks

The properties and modes of recognition of physiological DNAs associated with the four natural nucleobases might be extended, in principle, by the design of non-natural nucleobase derivatives. The goal is an expansion of the genetic alphabet, with the possible outcome of producing new DNAs with improved physical or biological properties. In this work, a new series of hetero-ring-expanded guanine analogs are proposed, and their relevant structural characteristics and electronic properties are determined by density functional theory. The stabilities of the decamer DNA duplexes (dn.dC)10 (where n represents the corresponding expanded guanine analog designed here) are also examined, using molecular dynamics. The simulations show that the designed motifs can form stable DNA-like structures. We determined the pairing energies for the Watson-Crick (WC) hydrogen-bonded dimers between the expanded G-analogs and the natural C, and found that the pairing energies are close to those of the natural GC pair. The calculated adiabatic ionization potentials (IPs) of the size-expanded guanine analogs and their base pairs, and the corresponding vertical ionization potentials, show that some are distinctly smaller than the corresponding natural versions. The HOMO-LUMO energy gaps for most of the size-expanded guanine analogs and their WC base pairs are considerably lower than those of the corresponding natural base and base pairs. Thus, the expanded G bases may be considered as DNA genetic motifs, and they may serve as building blocks for potential biological applications and the development of molecular electronic devices.     

Parallel DNA double helices incorporating isoG or m5isoC bases studied by FTIR, CD and molecular modeling

FTIR spectroscopy has been used to follow the formation of parallel stranded DNA duplexes incorporating isoG or m5isoC bases and determine their base pairing scheme. The results are discussed in comparison with data concerning anti-parallel duplexes with comparable base composition and sequence. In duplexes containing A-T and isoG-C or m5isoC-G base pairs shifts of the thymine C2=O2 and C4=O4 carbonyl stretching vibrations (to lower and higher wavenumbers, respectively, when compared to their positions in classical cis Watson-Crick (WC) base pairs) reflect the formation of trans Watson-Crick A-T base pairs. All carbonyl groups of cytosines, m5isocytosines, guanines and isoguanines are found to be involved in hydrogen bonds, indicative of the formation of isoG-C and m5isoC-G base pairs with three hydrogen bonds. Molecular modeling shows that both structures form regular right handed helices with C2'endo sugar puckers. The role of the water content on the helical conformation of the parallel duplexes has been studied by FTIR and CD. It is found that a conformational transition similar to the B --> A transition observed for anti-parallel duplexes induced by a decrease of the water content of the samples can occur for these parallel duplexes. Their helical flexibility has been evidenced by FTIR studies on hydrated films by the emergence of absorption bands characteristic of A type geometry, in particular by an S-type --> N-type repuckering of the deoxyribose. All sugars in the parallel duplex with alternating d(isoG-A)/d(C-T) sequence can adopt an N-type geometry in low water content conditions. The conformational transition of the parallel hairpin duplex with alternating d(isoG-A)/d(C-T) sequence was followed by circular dichroism in water/trifluoroethanol solutions and its free energy at 0 degrees C was estimated to be 6.6 +/- 0.3 kcal mol(-1).     

The identification of vicinally substituted cyclohexane isomers in their mixtures by 1H and 13C NMR spectroscopy

The radical addition reactions of organobromine compounds, XBr (X = CH2COOMe, PhCH2, CHBr2 and CCl3) with cyclohexene afforded mixtures of cis/trans isomer pairs of 1-X-2-Br-cyclohexanes. In addition to benzyl benzoyloxy derivatives are formed also, when benzoyl peroxide is used as an initiator. Owing to the great difficulties in separating these cis/trans isomer pairs, they are identified directly in their mixtures by NMR spectroscopy. In addition to one-dimensional (ID) 1H, proton decoupled 13C and DEPT-135, also two-dimensional (2D) 13C-13C INADEQUATE as well as 1H-13C HMQC experiments have been used in assigning the signals of each compound in their mixtures. The identification of each isomer was based on comparison of experimental 3J(H,H) coupling constants with theoretical ones based on the well-known Karplus type relationship. The more stable conformation for each isomer was estimated using the semiempirical AM1 molecular orbital method. The calculations support the isomer pair elucidations.     

Accurate interaction energies of hydrogen-bonded nucleic acid base pairs

Hydrogen-bonded nucleic acids base pairs substantially contribute to the structure and stability of nucleic acids. The study presents reference ab initio structures and interaction energies of selected base pairs with binding energies ranging from -5 to -47 kcal/mol. The molecular structures are obtained using the RI-MP2 (resolution of identity MP2) method with extended cc-pVTZ basis set of atomic orbitals. The RI-MP2 method provides results essentially identical with the standard MP2 method. The interaction energies are calculated using the Complete Basis Set (CBS) extrapolation at the RI-MP2 level. For some base pairs, Coupled-Cluster corrections with inclusion of noniterative triple contributions (CCSD(T)) are given. The calculations are compared with selected medium quality methods. The PW91 DFT functional with the 6-31G basis set matches well the RI-MP2/CBS absolute interaction energies and reproduces the relative values of base pairing energies with a maximum relative error of 2.6 kcal/mol when applied with Becke3LYP-optimized geometries. The Becke3LYP DFT functional underestimates the interaction energies by few kcal/mol with relative error of 2.2 kcal/mol. Very good performance of nonpolarizable Cornell et al. force field is confirmed and this indirectly supports the view that H-bonded base pairs are primarily stabilized by electrostatic interactions.     

Nucleic acid base pair dynamics: the impact of sequence and structure using free-energy calculations

Molecular dynamics free-energy calculations of base pair opening within double helical DNA and RNA are used to explain why A-tracts (oligo-adenine repeats) greatly increase the lifetimes of AT base pairs, whereas the structural and the chemical changes involved in passing from B-DNA to A-RNA have comparatively small effects.     

Hydration and stability of nucleic acid bases and base pairs

Empirical, quantum chemical calculations and molecular dynamics simulations of the role of a solvent on tautomerism of nucleic acid bases and structure and properties of nucleic acid base pairs are summarized. Attention was paid to microhydrated (by one and two water molecules) complexes, for which structures found by scanning of empirical potential surfaces were recalculated at a correlated ab initio level. Additionally, isolated as well as mono- and dihydrated H-bonded, T-shaped and stacked structures of all possible nucleic acid base pairs were studied at the same theoretical levels. We demonstrate the strong influence of a solvent on the tautomeric equilibrium between the tautomers of bases and on the spatial arrangement of the bases in a base pair. The results provide clear evidence that the prevalence of either the stacked or hydrogen-bonded structures of the base pairs in the solvent is not determined only by its bulk properties, but rather by specific hydrophilic interactions of the base pair with a small number of solvent molecules.