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.     

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.     

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 neighboring bases on base-pair stacking orientation: a molecular dynamics study

It is generally believed that base-pair stacking interaction in DNA double helix is one of the strongest interactions that governs sequence directed structural variability. However, X-ray crystal structures of some base-paired doublet sequences have been seen to adopt different structures when flanked by different base-pairs. DNA crystal database, however, is still too small to make good statistical inference about effect of such flanking residues. Influence of neighboring residue on the local helical geometry of a base-paired doublet in B-DNA has been investigated here using molecular dynamics simulation. We have generated ensembles of structures for d(CA).d(TG) and d(AA).d(TT) base-paired doublets located at the centers of d(CGCGCAAAGCG).d(CGCTTTGCGCG) and d(CGCGAAAACGCG).d(CGCGTTTTCGCG) sequences along with their analogs by varying the bases either at 5'- or 3'- position to the central doublet. Comparison of base paired doublet parameters for the ensembles of structures show that stacking geometry of d(CA).d(TG) doublet depends on some of the flanking base-pairs. On the other hand d(AA).d(TT) doublet remains nearly unperturbed when the flanking residues are altered.     

Sequence dependence of DNA conformational flexibility

By using conformational free energy calculations, we have studied the sequence dependence of flexibility and its anisotropy along various conformational variables of DNA base pairs. The results show the AT base step to be very flexible along the twist coordinate. On the other hand, homonucleotide steps, GG(CC) and AA(TT), are among the most rigid sequences. For the roll motion that would correspond to a bend, the TA step is most flexible, while the GG(CC) step is least flexible. The flexibility of roll is quite anisotropic; the ratio of fluctuations toward the major and minor grooves is the largest for the GC step and the smallest for the AA(TT) and CG steps. Propeller twisting of base pairs is quite flexible, especially of A.T base pairs; propeller twist can reach 19 degrees by thermal fluctuation. We discuss the effect of electrostatic parameters, comparison with available experimental results, and biological relevance of these results.     

Structural basis of DNA flexibility

It has been recently shown by us, on the basis of crystal structure database that the flexibility of B-DNA double helices depends significantly on their base sequence. Our model building studies further indicated that the existence of bifurcated cross-strand hydrogen bonds between successive base pairs is possibly the main factor behind the sequence directed DNA flexibility. These cross-strand hydrogen bonds are, of course, weaker than the usual Watson-Crick hydrogen bonds and their bond geometry is characterized by relatively larger bond lengths and smaller bond angles. We have tried to improve our model structures by incorporating non-planarity of the amino groups in DNA bases due to the presence of lone pair electrons at the nitrogen atoms. Energy minimization studies have been carried out by using different quantum chemical methods, whereby it is found that in all cases of N-H....O type cross-strand hydrogen bonds, the bond geometry improves significantly. In the cases of N-H....N type hydrogen bonds, however, no such consistent improvements can be noticed. Perhaps the true picture would emerge only if all the other interactions present in the DNA macromolecule could be appropriately taken into account.     

Geometric criteria of hydrogen bonds in proteins and identification of "bifurcated" hydrogen bonds

Empirical criteria for identification of hydrogen bonds were analyzed to produce a set of geometrically consistent criteria. For a data set of 30 structures, application of a set of purely geometrical criteria, along with exclusion of abnormal backbone conformations, also excluded a common interaction of Ser/Thr side chains with Asp/Glu side chains ([ST]/[DE] pairs). These interactions were termed "bifurcated hydrogen bonds", which implies delocalization of a positively charged hydrogen of hydroxyl between the two acceptor atoms of the carboxylic group. These "bifurcated" interactions are among the most common packing patterns for [ST]/[DE] pairs of side chains. Therefore, the identification of hydrogen bonds cannot be based on geometrical criteria only and requires introduction of some physico-chemical criteria.     

Solution structure of an A-tract DNA bend

The solution structure of a DNA dodecamer d(GGCAAAAAACGG)/d(CCGTTTTTTGCC) containing an A-tract has been determined by NMR spectroscopy with residual dipolar couplings. The structure shows an overall helix axis bend of 19 degrees in a geometry consistent with solution and gel electrophoresis experiments. Fourteen degrees of the bending occurs in the GC regions flanking the A-tract. The remaining 5 degrees is spread evenly over its six AT base-pairs. The A-tract is characterized by decreasing minor groove width from the 5' to the 3' direction along the A strand. This is a result of propeller twist in the AT pairs and the increasing negative inclination of the adenine bases at the 3' side of the run of adenine bases. The four central thymine bases all have negative inclination throughout the A-tract with an average value of -6.1 degrees. Although this negative inclination makes the geometry of the A-tract different from all X-ray structures, the proton on N6 of adenine and the O4 of thymine one step down the helix are within distance to form bifurcated hydrogen bonds. The 5' bend of 4 degrees occurs at the junction between the GC flank and the A-tract through a combination of tilt and roll. The larger 3' bend, 10 degrees, occurs in two base steps: the first composed of tilt, -4.1 degrees, and the second a combination of tilt, -4.2 degrees, and roll, 6.0 degrees. This second step is a direct consequence of the change in inclination between an adjacent cytosine base, which has an inclination of -12 degrees, and the next base, a guanine, which has 3 degrees inclination. This bend is a combination of tilt and roll. The large change in inclination allows the formation of a hydrogen bond between the protons of N4 of the 3' cytosine and the O6 of the next 3' base, a guanine, stabilizing the roll component in the bend. These structural features differ from existing models for A-tract bends.For comparison, we also determined the structure of the control sequence, d(GGCAAGAAACGG)/d(CCGTTTCTTGCC), with an AT to GC transition in the center of the A-tract. This structure has no negative inclination in most of the bases within the A-tract, resulting in a bend of only 9 degrees. When ligated in phase, the control sequence has nearly normal mobility in gel electrophoresis experiments.     

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.     

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.     

Crystallographic study of one turn of G/C-rich B-DNA

The DNA decamer d(CCAGGCCTGG) has been studied by X-ray crystallography. At a nominal resolution of 1.6 A, the structure was refined to R = 16.9% using stereochemical restraints. The oligodeoxyribonucleotide forms a straight B-DNA double helix with crystallographic dyad symmetry and ten base-pairs per turn. In the crystal lattice, DNA fragments stack end-to-end along the c-axis to form continuous double helices. The overall helical structure and, notably, the groove dimensions of the decamer are more similar to standard, fiber diffraction-determined B-DNA than A-tract DNA. A unique stacking geometry is observed at the CA/TG base-pair step, where an increased rotation about the helix axis and a sliding motion of the base-pairs along their long axes leads to a superposition of the base rings with neighboring carbonyl and amino functions. Three-center (bifurcated) hydrogen bonds are possible at the CC/GG base-pair steps of the decamer. In their common sequence elements, d(CCAGGCCTGG) and the related G.A mismatch decamer d(CCAAGATTGG) show very similar three-dimensional structures, except that d(CCAGGCCTGG) appears to have a less regularly hydrated minor groove. The paucity of minor groove hydration in the center of the decamer may be a general feature of G/C-rich DNA and explain its relative instability in the B-form of DNA.     

Recognition of hairpin-containing single-stranded DNA by oligonucleotides containing internal acridine derivatives.

Oligodeoxynucleotides with an internal intercalating agent have been targeted to single-stranded sequences containing hairpin structures. The oligonucleotide binds to nonadjacent single-stranded sequences on both sides of the hairpin structure in such a way as to form a three-way junction. The acridine derivative is inserted at a position that allows it to interact with the three-way junction. The melting temperature (Tm) of complexes formed between the hairpin-containing target and oligonucleotides containing one internal acridine derivative was higher than that obtained with the same target and an unmodified oligonucleotide (DeltaTm = +13 degrees C). The internal acridine provided the oligonucleotide with a higher affinity than covalent attachment to the 5' end. Oligonucleotides could also be designed to recognize a hairpin-containing single-stranded nucleic acid by formation of Watson-Crick hydrogen bonds with a single-stranded part and Hoogsteen hydrogen bonds with the stem of the hairpin. An internal acridine derivative was introduced at the junction between the two domains, the double helix domain with Watson-Crick base pairs and the triple helix domain involving Hoogsteen base triplets in the major groove of the hairpin stem. Oligonucleotides with an internal acridine or an acridine at their 5' end have similar binding affinities for the stem-loop-containing target. The bis-modified oligonucleotide containing two acridines, one at the 5' end and one at an internal site, did not exhibit a higher affinity than the oligonucleotides with only one intercalating agent. The design of oligonucleotides with an internal intercalating agent might be of interest to control gene expression through recognition of secondary structures in single-stranded targets.     

Thermodynamics of DNA duplexes with adjacent G.A mismatches.

The sequence 5'-d(ATGAGCGAAT) forms a very stable self-complementary duplex with four G.A mismatch base pairs (underlined) out of ten total base pairs [Li et al. (1991) Proc. Natl. Acad. Sci. U.S.A. 88, 26-30]. The conformation is in the general B-family and is stabilized by base-pair hydrogen bonding of an unusual type, by favorable base dipole orientations, and by extensive purine-purine stacking at the mismatched sites. We have synthesized 13 decamers with systematic variations in the sequence above to determine how the flanking sequences, the number of G.A mismatches, and the mismatch sequence order (5'-GA-3' or 5'-AG-3') affect the duplex stability. Changing A.T to G.C base pairs in sequences flanking the mismatches stabilizes the duplexes, but only to the extent observed with B-form DNA. The sequence 5'-pyrimidine-GA-purine-3', however, is considerably more stable than 5'-purine-GA-pyrimidine-3'. The most stable sequences with two pairs of adjacent G.A mismatches have thermodynamic parameters for duplex formation that are comparable to those for fully Watson-Crick base-paired duplexes. Similar sequences with single G.A pairs are much less stable than sequences with adjacent G.A mismatches. Reversing the mismatch order from 5'-GA-3' to 5'-AG-3' results in an oligomer that does not form a duplex. These results agree with predictions from the model derived from NMR and molecular mechanics and indicate that the sequence 5'-pyrimidine-GA-purine-3' forms a stable conformational unit that fits quite well into a B-form double helix.(ABSTRACT TRUNCATED AT 250 WORDS)     

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

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

Silver(I)-mediated Hoogsteen-type base pairs

Metal-mediated Hoogsteen-type base pairs are useful for the construction of DNA duplexes containing contiguous stretches of metal ions along the helical axis. To fine-tune the stability of such base pairs and the selectivity toward different metal ions, the availability of a selection of artificial nucleobases is highly desirable. In this study, we follow a theoretical approach utilizing dispersion-corrected density functional methods to evaluate a variety of artificial nucleobases as candidates for metal-mediated Hoogsteen-type base pairs. We focus on silver(I)-mediated Hoogsteen- and reverse Hoogsteen-type base pairs formed between 1-deaza- and 1,3-dideazapurine-derived nucleobases, respectively, and cytosine. Apart from two coordinative bonds, these base pairs are stabilized by a hydrogen bond. We elucidate the impact of different substituents at the C6 position and the presence or absence of an endocyclic N3 nitrogen atom on the overall stability of a base pair and concomitantly on the strength of the hydrogen and coordinative bonds. All artificial base pairs investigated in this study are less stable than the experimentally established benchmark base pair C–Ag⁺–G. The base pair formed from 1,3-dideaza-6-methoxypurine is isoenergetic to the experimentally observed C–Ag⁺–C base pair. This makes 1,3-dideaza-6-methoxypurine a promising candidate for the use as an artificial nucleobase in DNA.     

Crystal structure of a B-DNA dodecamer containing inosine, d(CGCIAATTCGCG), at 2.4 A resolution and its comparison with other B-DNA dodecamers.

The crystal structure of the dodecamer, d(CGCIAATTCGCG), has been determined at 2.4 A resolution by molecular replacement, and refined to an R-factor of 0.174. The structure is isomorphous with that of the B-DNA dodecamer, d(CGCGAATTCGCG), in space group P2(1)2(1)2(1) with cell dimensions of a = 24.9, b = 40.4, and c = 66.4 A. The initial difference Fourier maps clearly indicated the presence of inosine instead of guanine. The structure was refined with 44 water molecules, and compared to the parent dodecamer. Overall the two structures are very similar, and the I:C forms Watson-Crick base pairs with similar hydrogen bond geometry to the G:C base pairs. The propeller twist angle is low for I4:C21 and relatively high for the I16:C9 base pair (-3.2 degrees compared to -23.0 degrees), and the buckle angles alter, probably due to differences in the contacts with symmetry related molecules in the crystal lattice. The central base pairs of d(CGCIAATTCGCG) show the large propeller twist angles, and the narrow minor groove that characterize A-tract DNA, although I:C base pairs cannot form the major groove bifurcated hydrogen bonds that are possible for A:T base pairs.     

Observation of a unique pattern of bifurcated hydrogen bonds in the crystal structures of the N-glycoprotein linkage region models

Elucidation of the intra- and intermolecular carbohydrate-protein interactions would greatly contribute toward obtaining a better understanding of the structure-function correlations of the protein-linked glycans. The weak interactions involving C-H...O have recently been attracting immense attention in the domain of biomolecular recognition. However, there has been no report so far on the occurrence of C-H...O hydrogen bonds in the crystal structures of models and analogs of N-glycoproteins. We present herein an analysis of C-H...O interactions in the crystal structures of all N-glycoprotein linkage region models and analogs. The study reveals a cooperative network of bifurcated hydrogen bonds consisting of N-H...O and C-H...O interactions seen uniquely for the models. The cooperative network consists of two antiparallel chains of bifurcated hydrogen bonds, one involving N1-H, C2'-H and O1' of the aglycon moiety and the other involving N2-H, C1-H and O1' of the sugar. Such bifurcated hydrogen bonds between the core glycan and protein are likely to play an important role in the folding and stabilization of proteins.     

DNA curvature does not require bifurcated hydrogen bonds or pyrimidine methyl groups.

Short tracts of the homopolymer dA.dT confer intrinsic curvature on the axis of the DNA double helix. This phenomenon is assumed to be a consequence of such tracts adopting a stable B'-DNA conformation that is distinct from B-form structure normally assumed by other DNA sequences. The more stable B' structure of dA.dT tracts has been attributed to several possible stabilizing factors: (1) optimal base stacking interactions consequent upon the high propeller twist, (2) bifurcated hydrogen bonds between adjacent dA.dT base-pairs, (3) stacking interactions involving the dT methyl groups, and finally (4) a putative spine of ordered water molecules in the minor groove. DNA oligodeoxynucleotides have been synthesized that enable these hypotheses to be tested; of particular interest is the combination of effects due to bifurcation (2) and methylation of the pyrimidines nucleotides (3). The data indicate that neither bifurcated hydrogen bonds nor pyrimidine methyl groups nor both are essential for DNA curvature. The data further suggest that the influence of the minor groove spine of hydration on the B'-formation is small. The experiments favor the hypothesis that base stacking interactions are the dominant force in stabilizing the B'-form structure.     

Atomic resolution analysis of a 2:1 complex of CpG and acridine orange.

Cytidylyl-3', 5'-guanosine and acridine orange crystallize in a highly-ordered triclinic lattice which diffracts X-rays to 0.85 angstrom resolution. The crystal structure has been solved and refined to a residual factor of 9.5%. The two dinucleoside phosphate molecules form an antiparallel double helix with the acridine orange intercalated between them. The two base pairs of the double helical fragment have a twist angle of 10 degrees and it is found to have a C3' endo-(3', 5')-C2' endo mixed sugar puckering along the nucleotide backbone as has been observed for other simple intercalator complexes. Twenty-five water molecules have been located in the lattice together with a sodium ion. The intercalator double helical fragments form sheets which are held together by van der Waals interactions in one direction and hydrogen bonding interactions in the other. The crystal lattice contains aqueous channels in which sixteen water molecules are hydrogen bonded to the nucleotide, none to the intercalator, five water molecules are coordinated about the sodium ion and four water molecules bind solely to other water molecules. The bases in the base pairs have a dihedral angle of 7 to 8 degrees between them.