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.     

Studies on the cross stand hydrogen bonds in DNA double helices.

A systematic analysis has been carried out to examine all the stereochemically possible bifurcated hydrogen bonds including those of cross strand type between propeller twisted base pairs in DNA double helices by stereochemical considerations involving base pairs alone and by molecular mechanics studies on dimer and trimer duplexes. The results show that there are limited number of combinations of adjacent base pairs that would facilitate bifurcated cross-strand hydrogen bond (CSH). B-type helices concomitant with negative propeller twist seem to be more favored for the occurrence of CSH than canonical A-type helices because of slide in the latter. The results also demonstrate that helices with appropriate sequences may possess continuous run of these propeller twist driven cross strand hydrogen bonds indicating that they may in fact be considered as yet another general structural feature of DNA helices.     

Sequence directed flexibility of DNA and the role of cross-strand hydrogen bonds

Persistence length and torsional rigidity for different B-DNA sequences have been calculated by analysing crystal structure database. The values of these parameters for mixed sequence DNA are in good agreement with those estimated by others. Persistence lengths for the homopolymeric sequences, namely poly(dA).poly(dT) and poly(dG).poly(dC), are significantly large compared to those of others as expected from the inability of these sequences to form nucleosome under normal conditions. The heteropolymeric sequences poly(dA-dC).poly(dG-dT) and poly(dG-dC).poly(dG-dC), on the other hand, have smaller persistence lengths. This implies larger flexibility of the d(AC).d(GT), d(CA).d(TG), d(GC).d(GC) and d(CG).d(CG) doublets, some of which constitute the genetic disease forming triplet repeats d(CTG).d(CAG) and d(CGG).d(CCG). Thus it is expected that these triplet repeat sequences are also flexible and wrap around the histone octamer efficiently. Persistence length calculations also indicate larger flexibility for these triplet repeat sequences. Furthermore, our computations reveal that the rigidity of a given DNA sequence is controlled by its ability to form cross-strand bifurcated hydrogen bonds between the successive base pairs. Molecular orbital calculations suggest that these hydrogen bonds are generally extended with bond lengths around 3A.     

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.     

A two-stranded template-based approach to G.(C-A) triad formation: designing novel structural elements into an existing DNA framework

We have designed a DNA sequence, d(G-G-G-T-T-C-A-G-G), which dimerizes to form a 2-fold symmetric G-quadruplex in which G(syn). G(anti).G(syn).G(anti) tetrads are sandwiched between all trans G. (C-A) triads. The NMR-based solution structural analysis was greatly aided by monitoring hydrogen bond alignments across N-H...N and N-H...O==C hydrogen bonds within the triad and tetrad, in a uniformly ((13)C,(15)N)-labeled sample of the d(G-G-G-T-T-C-A-G-G) sequence. The solution structure establishes that the guanine base-pairs with the cytosine through Watson-Crick G.C pair formation and with adenine through sheared G.A mismatch formation within the G.(C-A) triad. A model of triad DNA was constructed that contains the experimentally determined G.(C-A) triad alignment as the repeating stacked unit.     

Effect of flanking residues on CA and AA dinucleotides: some rationale

Effect of flanking base pairs on CA and AA dinluceotide step-geometry has been studied using molecular dynamics method. Sixteen dodecameric sequences are constructed for each doublet with all possible bases at their 5' and 3' positions along with their complementary sequences. Structural parameter, formation of Extra Watson Crick bifurcated hydrogen bonds along or across the strands and effect of sodium ions are studied for these sequences. It is found that geometry of CA doublet step is perturbed by the neighboring base pairs, which might be due to inherent flexibility of the step. Flexible character of CA step is reflected in its low bifurcated hydrogen bond formation capability and lower preference of sodium ions to enter in minor or major grooves. AA step on the other hand is quite rigid according to different structural parameters and respond much less to environmental changes due to formation of strong Extra Watson-Crick hydrogen bonds.     

Studies On The Cross Strand Hydrogen Bonds In DNA Double Helices

Abstract

A systematic analysis has been carried out to examine all the stereochemically possible bifurcated hydrogen bonds including those of cross strand type between propeller twisted base pairs in DNA double helices by stereochemical considerations involving base pairs alone and by molecular mechanics studies on dimer and trimer duplexes. The results show that there are limited number of combinations of adjacent base pairs that would facilitate bifurcated cross- strand hydrogen bond (CSH). B-type helices concomitant with negative propeller twist seem to be more favored for the occurrence of CSH than canonical A-type helices because of slide in the latter. The results also demonstrate that helices with appropriate sequences may possess continuous run of these propeller twist driven cross strand hydrogen bonds indicating that they may infact be considered as yet another general structural feature of DNA helices.     

Effect of Flanking Residues on CA and AA Dinucleotides: Some Rationale

Abstract

Effect of flanking base pairs on CA and AA dinluceotide step-geometry has been studied using molecular dynamics method. Sixteen dodecameric sequences are constructed for each doublet with all possible bases at their 5′ and 3′ positions along with their complementary sequences. Structural parameter, formation of Extra Watson Crick bifurcated hydrogen bonds along or across the strands and effect of sodium ions are studied for these sequences. It is found that geometry of CA doublet step is perturbed by the neighboring base pairs, which might be due to inherent flexibility of the step. Flexible character of CA step is reflected in its low bifurcated hydrogen bond formation capability and lower preference of sodium ions to enter in minor or major grooves. AA step on the other hand is quite rigid according to different structural parameters and respond much less to environmental changes due to formation of strong Extra Watson-Crick hydrogen bonds.     

Hydrogen bonds in crystal structures of amino acids, peptides and related molecules

The results of a survey of 439 hydrogen bonds in 95 recently determined crystal structures of amino acids, peptides and related molecules suggest that the following generalizations hold true for linear (angle X-H---Y greater than 150 degrees) hydrogen bonds. (1) The charge on the acceptor group does not influence the length of a hydrogen bond. (2) For a given acceptor group, the hydrogen bond lengths increase in the order imidazolium N--H less than ammonium N-H less than guanidinium N-H; this order holds true for oxygen anion acceptor groups. Cl-ions and the uncharged oxygen of water molecules. (3) The uncharged imidazole N-H group forms shorter hydrogen than the amide N-H GROUP. (4) The carboxyl O-H groups form shorter hydrogen bonds than other hydroxyl groups. (5) The hydrogen bonds involving a halogen ion are longer than hydrogen bonds with other acceptors when corrected for their longer van der Walls radii. The observed differences between the lengths of hydrogen bonds formed by different donor and acceptor groups in amino acids and peptides, imply differences in the energetics of their formation.     

Simulation of interactions between nucleic acid bases by refined atom-atom potential functions

Energy of interaction between nitrogen bases of nucleic acid has been calculated as a function of parameters determining the mutual position of two bases. Refined atom-atom potential functions are suggested. These functions contain terms proportional to the first (electrostatics), sixth (or tenth for the atoms forming a hydrogen bond) and twelfth (repulsion of all atoms) powers of interatomic distance. Calculations have shown that there are two groups of minima of the base interaction energy. The minima of the first group correspond to coplanar arrangement of the base pairs and hydrogen bond formation. The minima of the second group correspond to the position of bases one above the other in almost parallel planes. There are 28 energy minima corresponding to the formation of coplanar pairs with two (three for the G:C pair) almost linear N-H . . . O and (or) N-H . . . N hydrogen bonds. The position of nitrogen bases paired by two such H-bonds in any crystal of nucleic acid component in polynucleotide complexes and in tRNA is close to the position in one of these minima. Besides, for each pair there are energy minima corresponding to the formation of a single N-H . . . O or N-H . . . N and one C-H . . . O or C-H . . . N hydrogen bond. The form of potential surface in the vicinity of minima has been characterized. The results of calculations agree with the experimental data and with more rigorous calculations based on quantum-mechanical approach.     

基于密度泛函理论研究Watson–Crick碱基对中的氢键特征(英文)

基于密度泛函理论(DFT)研究腺嘌呤、胸腺嘧啶、鸟嘌呤、胞嘧啶以及腺嘌呤胸腺嘧啶碱基对、鸟嘌呤胞嘧啶碱基对。在DFT-B3LYP/6-31G**水平上利用自然键轨道理论分析研究结果显示,互补碱基对的结构和电子特征有利于氢键的形成。本文中讨论几何结构、电子结构、分子轨道和能量对于氢键形成的影响。此研究结果将有助于更好的理解AT和GC碱基对中氢键与它们的结构特性之间的关系。     

The effect of anesthetics on hydrogen bonds An infrared study at low anesthetic concentrations

It is shown that a striking parallelism exists between the anesthetic potency of general halocarbon anesthetics and their influence on the hydrogen bond association constants in N-H…OC type hydrogen bonds, important for shaping the ion channels. It is further shown that the effect of potent anesthetics (which contain an acidic hydrogen) on the free/associated ratio in such hydrogen bonds is still significant at clinical anesthetic concentrations. It is argued that the results are in keeping with a pluralistic theory of anesthesia based on both hydrophobic and polar interactions.     

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.     

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.     

Theoretical analysis of noncanonical base pairing interactions in RNA molecules

Noncanonical base pairs in RNA have strong structural and functional implications but are currently not considered for secondary structure predictions. We present results of comparative ab initio studies of stabilities and interaction energies for the three standard and 24 selected unusual RNA base pairs reported in the literature. Hydrogen added models of isolated base pairs, with heavy atoms frozen in their 'away from equilibrium' geometries, built from coordinates extracted from NDB, were geometry optimized using HF/6-31G** basis set, both before and after unfreezing the heavy atoms. Interaction energies, including BSSE and deformation energy corrections, were calculated, compared with respective single point MP2 energies, and correlated with occurrence frequencies and with types and geometries of hydrogen bonding interactions. Systems having two or more N-H...O/N hydrogen bonds had reasonable interaction energies which correlated well with respective occurrence frequencies and highlighted the possibility of some of them playing important roles in improved secondary structure prediction methods. Several of the remaining base pairs with one N-H...O/N and/or one C-H...O/N interactions respectively, had poor interaction energies and negligible occurrences. High geometry variations on optimization of some of these were suggestive of their conformational switch like characteristics.     

Structure and assembly of designed beta-hairpin peptides in crystals as models for beta-sheet aggregation

De novo designed beta-hairpin peptides have generally been recalcitrant to crystallization. The crystal structures of four synthetic peptide beta-hairpins, Boc-Leu-Val-Val-DPro-Gly-Leu-Phe-Val-OMe (1), Boc-Leu-Phe-Val-DPro-Ala-Leu-Phe-Val-OMe (2), Boc-Leu-Val-Val-DPro-Aib-Leu-Val-Val-OMe (3), and Boc-Met-Leu-Phe-Val-DPro-Ala-Leu-Val-Val-Phe-OMe (4), are described. The centrally positioned DPro-Xxx segment promotes prime beta-turn formation, thereby nucleating beta-hairpin structures. In all four peptides well-defined beta-hairpins nucleated by central type II' DPro-Xxx beta-turns have been characterized by X-ray diffraction, providing a view of eight crystallographically independent hairpins. In peptides 1-3 three intramolecular cross-strand hydrogen bonds stabilized the observed beta-hairpin, with some fraying of the structures at the termini. In peptide 4, four intramolecular cross-strand hydrogen bonds stabilized the hairpin. Peptides 1-4 reveal common features of packing of beta-hairpins into crystals. Two-dimensional sheet formation mediated by intermolecular hydrogen bonds formed between antiparallel strands of adjacent molecule is a recurrent theme. The packing of two-dimensional sheets into the crystals is mediated in the third dimension by bridging solvents and interactions of projecting side chains, which are oriented on either face of the sheet. In all cases, solvation of the central DPro-Xxx peptide unit beta-turn is observed. The hairpins formed in the octapeptides are significantly buckled as compared to the larger hairpin in peptide 4, which is much flatter. The crystal structures provide insights into the possible modes of beta-sheet packing in regular crystalline arrays, which may provide a starting point for understanding beta-sandwich and cross-beta-structures in amyloid fibrils.     

Simulation of Interactions between Nucleic Acid Bases by Refined Atom-Atom Potential Functions

Abstract

Energy of interaction between nitrogen bases of nucleic acids has been calculated as a function of parameters determining the mutual position of two bases. Refined atom-atom potential functions are suggested. These functions contain terms proportional to the first (electrostatics), sixth (or tenth for the atoms forming a hydrogen bond) and twelfth (repulsion of all atoms) powers of interatomic distance. Calculations have shown that there are two groups of minima of the base interaction energy. The minima of the first group correspond to coplanar arrangement of the base pairs and hydrogen bond formation. The minima of the second group correspond to the position of bases one above the other in almost parallel planes. There are 28 energy minima corresponding to the formation of coplanar pairs with two (three for the G:C pair) almost linear N-H … O and (or) N-H … N hydrogen bonds. The position of nitrogen bases paired by two such H-bonds in any crystal of nucleic acid component, in polynucleotide complexes and in tRNA is close to the position in one of these minima. Besides, for each pair there are energy minima corresponding to the formation of a single N-H … O or N-H … N and one C-H … O or C-H … N hydrogen bond. The form of potential surface in the vicinity of minima has been characterized. The results of calculations agree with the experimental data and with more rigorous calculations based on quantum- mechanical approach.     

Theoretical analysis of noncanonical base pairing interactions in RNA molecules

Noncanonical base pairs in RNA have strong structural and functional implications but are currently not considered for secondary structure predictions. We present results of comparative ab initio studies of stabilities and interaction energies for the three standard and 24 selected unusual RNA base pairs reported in the literature. Hydrogen added models of isolated base pairs, with heavy atoms frozen in their ‘away from equilibrium’ geometries, built from coordinates extracted from NDB, were geometry optimized using HF/6-31G** basis set, both before and after unfreezing the heavy atoms. Interaction energies, including BSSE and deformation energy corrections, were calculated, compared with respective single point MP2 energies, and correlated with occurrence frequencies and with types and geometries of hydrogen bonding interactions. Systems having two or more N-H...O/N hydrogen bonds had reasonable interaction energies which correlated well with respective occurrence frequencies and highlighted the possibility of some of them playing important roles in improved secondary structure prediction methods. Several of the remaining base pairs with one N-H...O/N and/or one C-H...O/N interactions respectively, had poor interaction energies and negligible occurrences. High geometry variations on optimization of some of these were suggestive of their conformational switch like characteristics.     

Amino acid–base interactions: a three-dimensional analysis of protein–DNA interactions at an atomic level

To assess whether there are universal rules that govern amino acid–base recognition, we investigate hydrogen bonds, van der Waals contacts and water-mediated bonds in 129 protein–DNA complex structures. DNA–backbone interactions are the most numerous, providing stability rather than specificity. For base interactions, there are significant base–amino acid type correlations, which can be rationalised by considering the stereochemistry of protein side chains and the base edges exposed in the DNA structure. Nearly two-thirds of the direct read-out of DNA sequences involves complex networks of hydrogen bonds, which enhance specificity. Two-thirds of all protein–DNA interactions comprise van der Waals contacts, compared to about one-sixth each of hydrogen and water-mediated bonds. This highlights the central importance of these contacts for complex formation, which have previously been relegated to a secondary role. Although common, water-mediated bonds are usually non-specific, acting as space-fillers at the protein–DNA interface. In conclusion, the majority of amino acid–base interactions observed follow general principles that apply across all protein–DNA complexes, although there are individual exceptions. Therefore, we distinguish between interactions whose specificities are ‘universal’ and ‘context-dependent’. An interactive Web-based atlas of side chain–base contacts provides access to the collected data, including analyses and visualisation of the three-dimensional geometry of the interactions.     

Probing the role of hydrogen bonds in the stability of base pairs in double-helical DNA

Aromatic stacking and hydrogen bonding between nucleobases are two of the key interactions responsible for stabilization of DNA double-helical structures. The present work aims at defining the specific contributions of these interactions to the stability of individual base pairs in DNA. The two DNA double helices investigated are formed, respectively, by the palindromic base sequences 5'-dCCAACGTTGG-3' and 5'-dCGCAGATCTGCG-3'. The strength of the N==H...N inter-base hydrogen bond in each base pair is characterized from the measurement of the protium-deuterium fractionation factor of the corresponding imino proton using NMR spectroscopy. The structural stability of each base pair is evaluated from the exchange rate of the imino proton, measured by NMR. The results reveal that the fractionation factors of the imino protons in the two DNA double helices investigated fall within a narrow range of values, between 0.92 and 1.0. In contrast, the free energies of structural stabilization for individual base pairs span 3.5 kcal/mol, from 5.2 to 8.7 kcal/mol (at 15 degrees C). These findings indicate that, in the two DNA double helices investigated, the strength of N==H...N inter-base hydrogen bonds does not change significantly depending on the nature or the sequence context of the base pair. Hence, the variations in structural stability detected by proton exchange do not involve changes in the strength of inter-base hydrogen bonds. Instead, the results suggest that the energetic identity of a base pair is determined by the number of inter-base hydrogen bonds, and by the stacking interactions with neighboring base pairs.