N-H...O, O-H...O, and C-H...O hydrogen bonds in protein-ligand complexes: strong and weak interactions in molecular recognition

The characteristics of N-H...O, O-H...O, and C-H...O hydrogen bonds are examined in a group of 28 high-resolution crystal structures of protein-ligand complexes from the Protein Data Bank and compared with interactions found in small-molecule crystal structures from the Cambridge Structural Database. It is found that both strong and weak hydrogen bonds are involved in ligand binding. Because of the prevalence of multifurcation, the restrictive geometrical criteria set up for hydrogen bonds in small-molecule crystal structures may need to be relaxed in macromolecular structures. For example, there are definite deviations from linearity for the hydrogen bonds in protein-ligand complexes. The formation of C-H...O hydrogen bonds is influenced by the activation of the C(alpha)-H atoms and by the flexibility of the side-chain atoms. In contrast to small-molecule structures, anticooperative geometries are common in the macromolecular structures studied here, and there is a gradual lengthening as the extent of furcation increases. C-H...O bonds formed by Gly, Phe, and Tyr residues are noteworthy. The numbers of hydrogen bond donors and acceptors agree with Lipinski's "rule of five" that predicts drug-like properties. Hydrogen bonds formed by water are also seen to be relevant in ligand binding. Ligand C-H...O(w) interactions are abundant when compared to N-H...O(w) and O-H...O(w). This suggests that ligands prefer to use their stronger hydrogen bond capabilities for use with the protein residues, leaving the weaker interactions to bind with water. In summary, the interplay between strong and weak interactions in ligand binding possibly leads to a satisfactory enthalpy-entropy balance. The implications of these results to crystallographic refinement and molecular dynamics software are discussed.     

Molecular recognition motifs in cytidinium and 2'-deoxycytidinium salts with composite anions.

In the crystal structures of N3-protonated cytidinium and 2'-deoxycytidinium salts with composite XYn anions capable of accepting hydrogen bonds through their Y atoms, the dominating motif of cytosinium...anion interactions consists of a pair of hydrogen bonds donated from the N3+ -H protonation site and from the exoamino N4-H41 group cis to N3, and accepted by two Y centers of one anion. This multipoint recognition pattern is stable and robust and thus can be classified as a supramolecular synthon. In a broader group of N3-protonated, N1-substituted cytosinium salts with composite anions it occurs with 70% frequency. The C5 side of the cytosine ring mimics the N3+ -H type synthon and shows a propensity to form an analogous motif in which a C5-H5...Y hydrogen bond replaces the strong N3+ -H...Y interaction. Since the C-H...Y bond is much weaker, the secondary motif shows higher deformability and is less frequent (44%).     

Conformational flexibility in a triazole nucleoside derivative: 4-cyano-5-cyanomethyl-1-(2,3,5-tri-O-acetyl-beta-D-ribofuranosyl)-1,2,3-triazole

The crystal-structure determination of the molecular structure of the hydrophobic compound, 4-cyano-5-cyanomethyl-1-(2,3,5-tri-O-acetyl-beta-D-ribofuranosyl)-1,2,3-triazole, C16H17N5O7, provides us with two different conformations of ribofuranosyl moieties [(C2'-exo, C3'-endo) and C2'-exo] with two markedly different N-glycosidic angles. There are two molecules in the asymmetric unit of the crystal and the overall stereochemistry of the molecules are influenced predominantly by weak intramolecular bifurcated and trifurcated hydrogen bonds of the type C-H...O and C-H...N, where endo-H atoms attached to C2' and C3' are involved. The molecules in the crystal are interconnected with longer intermolecular bonds of the same type. There are empty channels (occupying 14.0% of the whole volume of the unit cell), which are extended along b-axis of the entire crystal.     

The propeller DNA conformation of poly(dA).poly(dT).

Physical properties of the DNA duplex, poly(dA).poly(dT) differ considerably from the alternating copolymer poly(dAT). A number of molecular models have been used to describe these structures obtained from fiber X-ray diffraction data. The recent solutions of single crystal DNA dodecamer structures with segments of oligo-A.oligo-T have revealed the presence of a high propeller twist in the AT regions which is stabilized by the formation of bifurcated (three-center) hydrogen bonds on the floor of the major groove, involving the N6 amino group of adenine hydrogen bonding to two O4 atoms of adjacent thymine residues on the opposite strand. Here we show that it is possible to incorporate the features of the single crystal analysis, specifically high propeller twist, bifurcated hydrogen bonds, and a narrow minor groove, as well as the close interstrand NMR signal between adenine HC2 and ribose HC1' of the opposite strand, into a model that is fully compatible with the diffraction data obtained from poly(dA).poly(dT).     

Synthesis and X-ray crystallographic investigation of N-(3-deoxy-3-acetamido-��-D-glycopyranosyl)alkanamides as analogs of N-glycoprotein linkage region

As part of our ongoing program aimed at understanding the structural significance of GlcNAcβAsn linkage conserved in all eukaryotic N-glycoproteins, the present study reports on the synthesis and X-ray crystal structures of N-(3-deoxy-3-acetamido-β-D-glycopyranosyl)acetamide (Glc3NAcβNHAc) and the corresponding propionamide (Glc3NAcβNHPr). Comparative analysis of these structures with those of the corresponding GlcNAc (C2 acetamido) compounds showed that the bifurcated anti-parallel pattern involving N-H…O and C-H…O hydrogen bonds, the hallmark feature of the N-glycoprotein models, GlcNAcβNHAc and GlcNAcβAsn, is absent in both the C3 acetamido analogs. The extended (anti) conformation of the amido aglycon moiety as defined by χ(2) seen in the case of C2 acetamido derivative, GlcNAcβNHPr, turns into gauche for the C3 acetamido analog (Glc3NAcβNHPr). This observation is consistent with the earlier work on the critical role of the C2-NHAc group of GlcNAcβAsn in controlling χ(2) at the linkage region of N-glycoproteins.     

Self base pairing in a complementary deoxydinucleoside monophosphate duplex: crystal and molecular structure of deoxycytidylyl-(3'-5')-deoxyguanosine

The molecular structure of ammonium deoxycytidylyl-(3'-5')-deoxyguanosine, crystallized from aqueous acetone near pH 4, was determined for X-ray diffraction data. The crystals were tetragonal, space group P43212 with a = b = 11.078 (1) A and c = 45.826 (4) A. The structure was solved by tangent expansion of phases based on a derived phosphorus position and refined to R = 0.060 by full matrix least squares. Molecules related by a 2-fold symmetry axis are connected by hydrogen bonds between the bases and form parallel right-handed duplexes. Pairs of cytosines share a proton at N(3) and are joined by three hydrogen bonds: N(4)-H...O(2)...H-N(4), and N(3)-H...N(3). Guanines are joined by two hydrogen bonds: N(2)-H...N(3) and N(3)...H-N(2). Base-stacking interactions within the duplex are weak with the cytosine and guanine ring planes inclined at 24 degrees to each other in each monomer. Despite the unusual arrangement of the molecules, the sugar phosphate backbone has the g-g- conformation normally associated with right-handed double helical structures. Conformational parameters of the nucleosides are also typical with both sugars C(2')-endo and glycosidic torsion angles 55 degrees for cytidine and 94 degrees for guanosine. The bonding geometry of the bases is influenced by hydrogen bonding and charge-transfer networks in the crystal lattice. The solvent molecules interact with the dimer in three fused circular hydrogen bonding domains with a single disordered ammonium cation per d(CpG) dimer. Parallels with the formation of self base pairs and their implications in molecular biology are discussed.     

Evaluation of C-H cdots,three dots,centered O hydrogen bonds in native and misfolded proteins

Non-traditional C-H cdots, three dots, centered Y hydrogen bonds, in which a carbon atom acts as the hydrogen donor and an electronegative atom Y (Y=N, O or S) acts as the acceptor, have been reported in proteins, but their importance in protein structures is not well established. Here, we present the results of three computational tests that examine the significance of C-H cdots, three dots, centered Y bonds involving the C(alpha) in proteins. First, we compared the number of C(alpha)-H cdots, three dots, centered Y bonds in native structures with two sets of compact, energy-minimized decoy structures. The decoy structures contain about as many C(alpha)-H cdots, three dots, centered Y bonds as the native structures, indicating that the constraints of chain connectivity and compactness can lead to incidental formation of C(alpha)-H cdots, three dots, centered Y bonds. Secondly, we examined whether short C(alpha)-H cdots, three dots, centered Y bonds have a tendency to be linear, as is expected for a cohesive hydrogen-bonding interaction. The native structures do show this trend, but so does one of the decoy sets, suggesting that this criterion is also not sufficient to indicate a stabilizing interaction. Finally, we examined the preference for C(alpha)-H cdots, three dots, centered Y bond donors to be near to strong hydrogen bond acceptors. In the native proteins, the alpha protons attract strong acceptors like oxygen atoms more than weak acceptors. In contrast, hydrogen bond donors in the decoy structures do not distinguish between strong and weak acceptors. Thus, any individual C(alpha)-H cdots, three dots, centered Y bond may be fortuitous and occur due to the polypeptide connectivity and compactness. Taken collectively, however, C(alpha)-H cdots, three dots, centered Y bonds provide a weakly cohesive force that stabilizes proteins.     

C-H...O hydrogen bonds in the nuclear receptor RARgamma--a potential tool for drug selectivity

Hydrogen bonds between polarized atoms play a crucial role in protein interactions and are often used in drug design, which usually neglects the potential of C-H...O hydrogen bonds. The 1.4 A resolution crystal structure of the ligand binding domain of the retinoic acid receptor RARgamma complexed with the retinoid SR11254 reveals several types of C-H...O hydrogen bonds. A striking example is the hydroxyl group of the ligand that acts as an H bond donor and acceptor, leading to a synergy between classical and C-H...O hydrogen bonds. This interaction introduces both specificity and affinity within the hydrophobic ligand pocket. The similarity of intraprotein and protein-ligand C-H...O interactions suggests that such bonds should be considered in rational drug design approaches.     

Crystal structure of N-(1-deoxy-beta-D-fructopyranos-1-yl)-L-proline-an Amadori compound

Here we report the crystal structure data on N-(1-deoxy-beta-D-fructopyranos-1-yl)-L-proline (Fru-Pro)-an Amadori compound. X-ray crystal and molecular structures of its two isomorphous crystalline forms, (Fru-Pro)xMeOH, C(11)H(19)NO(7)xCH(4)O (1a) and (Fru-Pro)x2H(2)O, C(11)H(19)NO(7)x2H(2)O (1b) were determined. In 1a and 1b the compound crystallizes as the beta-anomer with the overall geometry of Fru-Pro zwitterions being very similar. Fructose ring adopts the chair (2)C(5) conformation with the proline moiety bonded to equatorial C-1 atom and remaining in a trans-gauche (tg) orientation with respect to the sugar ring. The five-membered pyrrolidine ring adopts an envelope conformation, with the Cbeta atom puckered. Fructosyl and carboxylate groups are in bisectional and axial positions of pyrrolidine ring, respectively. The overall molecular geometry of Fru-Pro zwitterions, especially the relative orientation of sugar and amino acid moieties, is stabilized by intramolecular, three-centred N-H...O(Fru)/O(Pro) hydrogen bonds (with bifurcated acceptor) formed between aminium and hydroxyl/carboxylate groups. The packing diagrams are very similar in both 1a and 1b with the adjacent zwitterions linked to each other by the extensive network of O-H...O and C-H...O hydrogen bonds to form channels along the a-axis, filled up with solvent molecules.     

Strong and weak hydrogen bonds in protein-ligand complexes of kinases: a comparative study

Strong and weak hydrogen bonds between protein and ligand are analyzed in a group of 233 X-ray crystal structures of the kinase family. These kinases are from both eukaryotic and prokaryotic organisms. The dataset comprises of 44 sub-families, out of which 35 are of human origin and the rest belong to other organisms. Interaction analysis was carried out in the active sites, defined here as a sphere of 10 A radius around the ligand. A majority of the interactions are observed between the main chain of the protein and the ligand atoms. As a donor, the ligand frequently interacts with amino acid residues like Leu, Glu and His. As an acceptor, the ligand interacts often with Gly, and Leu. Strong hydrogen bonds N-H...O, O-H...O, N-H...N and weak bonds C-H...O, C-H...N are common between the protein and ligand. The hydrogen bond donor capacity of Gly in N-H...O and C-H...O interactions is noteworthy. Similarly, the acceptor capacity of main chain Glu is ubiquitous in several kinase sub-families. Hydrogen bonds between protein and ligand form characteristic hydrogen bond patterns (supramolecular synthons). These synthon patterns are unique to each sub-family. The synthon locations are conserved across sub-families due to a higher percentage of conserved sequences in the active sites. The nature of active site water molecules was studied through a novel classification scheme, based on the extent of exposure of water molecules. Water which is least exposed usually participates in hydrogen bond formation with the ligand. These findings will help structural biologists, crystallographers and medicinal chemists to design better kinase inhibitors.     

Functional role of C(alpha)-H...O hydrogen bonds between transmembrane alpha-helices in photosystem I

The crystal structure at 2.5A resolution of the membrane-intrinsic, homotrimeric photosystem I (PSI) isolated from the thermophilic cyanobacterium Synechococcus elongatus shows that each monomer is composed of 12 protein subunits of which nine are embedded in the membrane and feature a total of 34 transmembrane alpha-helices (TMH). Hence, PSI provides an ideal case to study "conventional" and C(alpha)-H...O hydrogen bonds between TMH engaged in intra- and intersubunit interactions. Of the total of 75 C(alpha)-H...O hydrogen bonds between TMHs, 72 are intrasubunit and only three are intersubunit. The two largest subunits PsaA and PsaB are each folded into 11 TMHs showing 29 and 24 intrasubunit C(alpha)-H...O hydrogen bonds, respectively, that are not distributed randomly but many of them flank chlorophyll a (Chl a) co-ordinating amino acids, suggesting stabilisation of these structural segments. As major constituent of the trimerisation domain, subunit PsaL is located next to the 3-fold axis relating the three monomers of PSI. PsaL features a unique number of 19 intrasubunit C(alpha)-H...O hydrogen bonds that connect two of its three TMHs but there are no intersubunit C(alpha)-H...O hydrogen bonds between the three PsaL. Of the three intersubunit C(alpha)-H...O hydrogen bonds, two are formed between PsaA and PsaB and one between PsaB and PsaM. The large number of 75 C(alpha)-H...O hydrogen bonds contrasts the 49 conventional hydrogen bonds, indicating that the former and van der Waals contacts determine association and orientation of TMHs in PSI.     

Inclusion complexes of V-amylose with undecanoic acid and dodecanol at atomic resolution: X-ray structures with cycloamylose containing 26 D-glucoses (cyclohexaicosaose) as host

Crystal structures are reported of cycloamylose containing 26 D-glucose residues (CA26, cyclohexaicosaose, C156H260O130) in complexes with undecanoic acid (CA26 x 2C10H21COOH x 34.95 H2O, orthorhombic P2(1)2(1)2(1), one CA26 and two bound undecanoic acids F1 and F2 in the asymmetric unit, resolution 0.95 angstroms) and with dodecanol ((CA26)(0.5) x C12H25OH x 32.0H2O, monoclinic C2, half a CA26 binding one dodecanol, A, in the asymmetric unit, resolution 1.0 angstroms). The macrocycle of CA26 is folded like the figure '8' into two 10 D-glucoses long left-handed V-amylose helices forming approximately 5A wide V-channels that are occupied by undecanoic acid (F1, F2) or dodecanol (A) as guest molecules. The functional head groups of the guests near the O(6) ends of the V-channels are hydrogen bonded with d-glucose O(6)n-H; the aliphatic termini beyond C(9) protrude from the O(2), O(3) ends. Parts of the aliphatic chains enclosed in the V-channels are all-trans except for one torsion angle each (approximately 130 degrees ) in undecanoic acid molecules F1 and F2. There are several (guest)C-H...O hydrogen bonds to O(4) and O(6) of CA26 in both complexes, and H...H van der Waals interactions with d-glucose C(3)-H and C(5)-H dominate. C(5)-H determine the position of the aliphatic chains of undecanoic acid F1 and of dodecanol A in contrast to F2 where both C(3)-H and C(5)-H contribute equally, probably because the V-channel is narrower than in F1 and in dodecanol. Complexes of polymeric V-amylose with fatty acids and alcohols studied by X-ray fiber diffraction could not provide the here described high resolution.     

Interaction geometry involving planar groups in protein-protein interfaces

The geometry of interactions of planar residues is nonrandom in protein tertiary structures and gives rise to conventional, as well as nonconventional (X--H...pi, X--H...O, where X = C, N, or O) hydrogen bonds. Whether a similar geometry is maintained when the interaction is across the protein-protein interface is addressed here. The relative geometries of interactions involving planar residues, and the percentage of contacts giving rise to different types of hydrogen bonds are quite similar in protein structures and the biological interfaces formed by protein chains in homodimers and protein-protein heterocomplexes--thus pointing to the similarity of chemical interactions that occurs during protein folding and binding. However, the percentage is considerably smaller in the nonspecific and nonphysiological interfaces that are formed in crystal lattices of monomeric proteins. The C--H...O interaction linking the aromatic and the peptide groups is quite common in protein structures as well as the three types of interfaces. However, as the interfaces formed by crystal contacts are depleted in aromatic residues, the weaker hydrogen bond interactions would contribute less toward their stability.     

C-H ... O hydrogen bonding in solutions of methylated nucleic acid base analogs as revealed by NMR

Formation and thermodynamic characteristics of C-H ... O hydrogen bonding of methylated uracils and caffeine have been studied by nmr along two lines. 1. The concentration and temperature dependencies of the PMR spectra of 1,3-dimethyluracil (m2 1,3Ura), 1,3-dimethylthymine (m2 1,3Thy), and 1,3,6-trimethyluracil (m3 1,3,6Ura) in chloroform at high concentrations of base analogs indicated the self-association of m2 1,3Ura and m2 1,3Thy via C(6)H ... O hydrogen bonding and the competitive formation of C-H ... O bonds between carbonyl oxygens and chloroform. The intermolecular interaction energy and the arrangement of molecules in the local minima of various m2 1,3Ura dimers were calculated by the method of atom-atom potentials. The deepest minimum for the m2 1,3Ura coplanar dimer corresponds to a C(6)-H ... O hydrogen-bond formation. 2. At low concentration of m2 1,3Ura and caffeine in CCl4, C(6)-H ... O bonding for m2 1,3Ura and C(8)-H ... O bonding for caffeine with oxygens of dimethyl sulfoxide (DMSO) and acetone were observed. The association constants of these complexes were obtained at different temperatures. The enthalpies delta H, of the m2 1,3Ura-DMSO, m2 1,3Ura-accetone, caffeine-DMSO, and caffeine-acetone complexes were -2 +/- 0.1 kcal/mol. The calculations showed that the deepest minimum of the caffeine-acetone coplanar complex corresponds to C(8)-H ... O bonding with energy of -3.5 kcal/mol and that of the m2 1,3Ura-acetone complexes corresponds to C(6)-H ... O bonding with energy of -3.4 kcal/mol. The approximate correction for the solvent effect provides good agreement of the experimental data with the calculations.     

C-H.O hydrogen bonds in minor groove of A-tracts in DNA double helices

Analysis of available B-DNA type oligomeric crystal structures as well as protein-bound DNA fragments (solved using data with resolution <2.6 A) indicates that in both data sets, a majority of the (3'-Ade) H2..O2(3'-Thy/Cyt) distances in AA.TT and GA.TC dinucleotide steps, are considerably shorter than their values in a uniform fibre model, and are smaller than their optimum separation distance. Since the electropositive C2-H2 group of adenine is in close proximity of the electronegative keto oxygen atoms of both pyrimidine bases in the antiparallel strand of the double-helical DNA structures, it suggests the possibility of intra-base-pair as well as cross-strand C-H..O hydrogen bonds in the minor groove. The C2-H2..O2 hydrogen bonds within the A.T base-pairs could be a natural consequence of Watson-Crick pairing. However, the close cross-strand interactions between the bases at the 3'-ends of the AA.TT and GA.TC steps arise due to the local sequence-dependent geometry of these steps. While the base-pair propeller twist in these steps is comparable to the fibre model, some of the other local parameters such as base-pair opening angle and inter-base-pair slide show coordinated changes, leading to these shorter C2-H2..O2 distances. Hence, in addition to the well-known minor groove hydration, it appears that favourable C2-H2..O2 cross-strand interactions may play a role in imparting a characteristic geometry to AA.TT and GA.TC steps, as well as An.Tn and GAn.TnC tracts, which leads to a narrow minor groove in these regions.     

Synthesis and characterization of 4,6-O-butylidene-N-(2-hydroxybenzylidene)-beta-D-glucopyranosylamine: crystal structures of 4,6-O-butylidene-alpha-D-glucopyranose, 4,6-O-butylidene-beta-D-glucopyranosylamine and 4,6-O-butylidene-N-(2-hydroxybenzylidene)-beta-D-glucopyranosylamine

4,6-O-Butylidene-N-(2-hydroxybenzylidene)-beta-D-glucopyranosylamine was synthesized and characterized using analytical, spectral and single-crystal X-ray diffraction methods. 1H and 13C NMR studies showed the presence of the beta-anomer, which has also been confirmed by the crystal structure. The molecular structure of this compound showed the presence of the tridentate ONO ligation-core. Both precursors, 4,6-O-butylidene-alpha-D-glucopyranose and 4,6-O-butylidene-beta-D-glucopyranosylamine were characterized using single crystal X-ray diffraction. The alpha-anomeric nature of the former and beta-anomeric nature of the latter were proposed based on 1H NMR studies and were confirmed by determining the crystal structures. In addition, the crystal structure of 4,6-O-butylidene-beta-D-glucopyranosylamine revealed the C-1-N-glycosylation. In all the three molecules, the saccharide unit exhibits a 4C(1) chair conformation. In the lattice, the molecules are connected by hydrogen-bond interactions. The conformation of 4,6-O-butylidene-N-(2-hydroxybenzylidene)-beta-D-glucopyranosylamine is stabilized via an O-H...N intramolecular interaction, and each molecule in the lattice interacts with three neighboring molecules through hydrogen bonds of the type O-H...O and C-H...O.     

Comparative structural analysis of cytidine, ethenocytidine and their protonated salts. II. IR spectral studies.

The IR spectra of crystalline cytidine (Cyd), ethenocytidine (epsilon Cyd), and their hydrochlorides (Cyd-Hcl and epsilon CyD-HCl) have been analyzed to determine the spectroscopic manifestations of the structural differences that were previously established for these nucleosides from X-ray studies. O,N-Deuteration of the samples turned out to be a successful approach to obtaining interpretable spectra. The analysis was carried out in three frequency ranges: (i) The 2600-1900 cm-1 range originating from the vO-D and VN-D vibrations. All intermolecular hydrogen bonds could be recognized here. The positions of the individual vO-D (vN-D) bands were correlated with the geometrical delta HB parameters presenting the strengths of hydrogen bonds in which these groups act as donors (ii) The 1750-1500 cm-1 region originating from the stretching vibrations of double bonds. All absorption bands in this region were interpreted in terms of electronic structures of the base fragments. (iii) The region of the C-H stretching vibrations of the base fragments (3200-3000 cm-1) and sugar moieties (3000-2800 cm-1). The Csp2-H vibrations also reflect the electronic structures of the base fragments, whereas the vCsp-H frequencies seem to be sensitive to etheno-bridging and to the presence of an intramolecular C6-H...05' hydrogen bond.     

Occurrence of bifurcated three-center hydrogen bonds in proteins

Analysis of 13 high-resolution protein X-ray crystal structures shows that 1204 (24%) of all the 4974 hydrogen bonds are of the bifurcated three-center type with the donor X-H opposing two acceptors A1, A2. They occur systematically in alpha-helices where 90% of the hydrogen bonds are of this type; the major component is (n + 4)N-H ... O = C(n) as expected for a 3.6(13) alpha-helix, and the minor component is (n + 4)N-H ... O = C(n + 1), as observed in 3(10) helices; distortions at the C-termini of alpha-helices are stabilized by three-center bonds. In beta-sheets 40% of the hydrogen bonds are three-centered. The frequent occurrence of three-center hydrogen bonds suggests that they should not be neglected in protein structural studies.     

Estimation of strength in different extra Watson-Crick hydrogen bonds in DNA double helices through quantum chemical studies

It was shown earlier, from database analysis, model building studies, and molecular dynamics simulations that formation of cross-strand bifurcated or Extra Watson-Crick hydrogen (EWC) bonds between successive base pairs may lead to extra rigidity to DNA double helices of certain sequences. The strengths of these hydrogen bonds are debatable, however, as they do not have standard linear geometry criterion. We have therefore carried out detailed ab initio quantum chemical studies using RHF/6-31G(2d,2p) and B3LYP/6-31G(2p,2d) basis sets to determine strengths of several bent hydrogen bonds with different donor and acceptors. Interaction energy calculations, corrected for the basis set superposition errors, suggest that N-H...O type bent EWC hydrogen bonds are possible along same strands or across the strands between successive base pairs, leading to significant stability (ca. 4-9 kcal/mol). The N-H...N and C-H...O type interactions, however, are not so stabilizing. Hence, consideration of EWC N-H...O H-bonds can lead to a better understanding of DNA sequence directed structural features.