Perturbation of Hydrogen Bonds in the Adenine…Thymine Base Pair by Na+, Mg2+, Ca2+ and NH4 + Cations

Abstract

Perturbation of the hydrogen bonds in the adenine…thymine base pair by Na⁺, Mg²⁺, Ca²⁺ and NH₄ ⁺ cations has been investigated by means of ab initio SCF calculations with the STO-3G basis set. The geometry of adenine…thymine, as well as those of the perturbed pairs were optimized. Approach of any cation to thymine at 06 leads to destabilization of the adenine…thy mine pair; divalent cations (Mg²⁺, Ca²⁺) have a profound effect on the structure of the base pair. The approach of a cation to other available sites (thymine: O2, adenine N1 and N3) leads, on the other hand, to stabilization of the base pair. If a water molecule is placed between the cation and the base pair, the structure and stability of the base pair are changed only negligibly.     

Molecular mechanical studies of base-pair opening in d(CGCGC):d(GCGCG), dG5·dC5, d(TATAT):d(ATATA), and dA5·dT5 in the B and Z forms of DNA

Molecular mechanical energy refinement of double-helical pentanucleotide tetra-phosphates, d(CGCGC):d(GCGCG), dG₅·dC₅, d(TATAT):d(ATATA), and dA₅ ·dT₅ geometries, are presented in order to examine the energy required to open the Nl(purine) …? N3(pyrimidine) distance (base-pair opening) of a Watson-Crick base pair from its normal value of 3 Å to a value of 6 Å. The structural consequences of forcing base-pair opening is sequence dependent. For both dA₅ ·dT₅ and d(TATAT):d(ATATA), forcing the Nl (AdeKN3 (Thy) distance of the central base pair to a value of 6 Å slides the bases perpendicular to the helix axis forming a low-energy non-Watson-Crick base pair having an adenine amine hydrogen …? thymine carbonyl oxygen hydrogen bond. The two GC sequences behave differently from both AT sequences and differently from each other. Forcing the Nl(Gua) …? N3(Cyt) distance to 6 Å leads to unconventional structures in which hydrogen bonds are formed between the separated bases and the bases above or below them. These structures appear to be trapped in true local minima 6–10 kcal/mol higher in energy than the Watson-Crick structures. Preliminary simulations on d(CGCGC):d(GCGCG) in the Z geometry suggest the reason the Z form may be more refractory to proton exchange than the B form, consistent with experimental observations.     

A Theoretical Model for the Binding of Cis-Pt(NH3)2 +2 to DNA

Abstract

The binding of cis-Pt(NH₃)₂B₁B2 to the bases B₁ and B₂, i.e., guanine (G), cytosine (C), adenine (A), and thymine (T), of DNA is studied theoretically. The components of the binding are analyzed and a model structure is proposed for the intrastrand binding to the dB,pdB₂ sequence of a kinked double helical DNA. Quantum mechanical calculations of the ligand binding energy indicates that cw-Pt(NH₃)₂ ⁺² (cis-PDA) binds to N7(G), N3(C), 02(C), 06(G), N3(A), N7(A), 04(T) and 02(T) in order of decreasing binding energy. Conformational analysis provides structures of kinked DNA in which adjacent bases chelate to cis-PDA. Only bending toward the major groove allows the construction of acceptable square planar complexes. Examples are presented for kinks of ?70° and ?40° at the receptor site to orient the base pairs for ligand binding to B, and B₂ to form a nearly square planar complex. The energies for complex formation of cis-PDA to the various intra-strand base sites in double stranded DNA are estimated. At least 32 kcal/mole separates the energetically favorable dGpdG·cis-PDA chelate from the dCpdG·cis-PDA chelate. All other possible chelate structures are much higher in energy which correlates with their lack of observation in competition with the preferred dGpdG chelate.

The second most favorable ligand energy occurs with N3(C). A novel binding site involving dC(N3)pdG(N7) is examined. Denaturation can result in an anti ? syn rotation of C about its glycosidic bond to place N3(C) in the major groove for intrastrand binding in duplex DNA. This novel intrastrand dCpdG complex and the most favored dGpdG structure are illustrated with stereographic projections.     

Investigation of the base-pairing structure of the anticodon hairpin from E. coli initiator tRNA by high-resolution nmr

The high-resolution nmr spectrum of the anticodon hairpin from E. coli tRNA^fMet has been obtained at a number of different temperatures. The positions of the resonances from interior Watson-Crick base pairs are well accounted for (within 0.1 ppm) by a semi-empirical ring current shift theory, but the terminal base pairs are susceptible to the exact orientation of adjacent bases in single-stranded regions. From a careful examination of the exact way in which resonances disappear at elevated temperatures, we conclude that melting in the nmr experiments occurs when the lifetime of a base pair is reduced to several milliseconds. On the basis of these experiments we are able to assign an nmr T_m to each individual base pair and these should be useful in interpreting the melting behavior of the intact molecule. An “extra” resonance is observed at ～11.3 ppm and, on the basis of its position and temperature sensitivity, it is tentatively assigned to the ring nitrogen proton of a “protected” U residue in the anticodon loop. A strong preference for stacking of a nonbase-paired A residue on an adjacent GC base pair is observed even at temperatures in excess of 52°C.     

Persistence analysis of the static and dynamical helix deformations of DNA oligonucleotides: Application to the crystal structure and molecular dynamics simulation of d(CGCGAATTCGCG)2

A theory and graphical presentation for the analysis of helix structure and deformations in oligonucleotides is presented. The parameters “persistence” and “flexibility” as defined in the configurational statistics of polymers of infinite length are reformulated at the oligonucleotide level in an extension of J. A. Schellman's method [(1974) Biopolymers, Vol. 17, pp. 217–226], and used as a basis for a systematic “Persistence Analysis” of the helix deformation properties for all possible subsequences in the structure. The basis for the analysis is a set of link vectors referenced to individual base pairs, and is limited to sequences exhibiting only perturbed rod-like behavior, i.e., below the threshold for supercoiling. The present application of the method is concerned with a physical model for the angular component of bending, so the link vectors are defined as the unit components of a global helix axis obtained by the procedure “Curves” of R. Lavery and H. Sklenar [(1988) J. Biomol. Struct. Dynam., Vol. 6, pp. 63–91; (1989) J. Biomol. Struct. Dynam., Vol. 6, pp. 655–667]. A discussion, of the relationship between global bending and relative orientation of base pairs is provided. Our approach is illustrated by analysis of some model oligonucleotide structures with intrinsic kinks, the crystal structure of the dodecamer d (CGCGAATTCGCG)₂, and the results of two molecular dynamics simulations on this dodecamer using two variations of the GROMOS force field. The results indicate that essentially all aspects of curvature in short oligonucleotides can be determined, such as the position and orientation of each bend, the sharpness or smoothness, and the location and linearity of subsequences. In the case of molecular dynamics simulations, where a Boltzmann ensemble of structures is analyzed, the spatial extent of the deformations (flexibility) is also considered. © 1993 John Wiley & Sons, Inc.     

Interactions between Al12X (X = Al, C, N and P) nanoparticles and DNA nucleobases/base pairs: implications for nanotoxicity

The interactions between neutral Al₁₂X(I _h ) (X = Al, C, N and P) nanoparticles and DNA nucleobases, namely adenine (A), thymine (T), guanine (G) and cytosine (C), as well as the Watson−Crick base pairs (BPs) AT and GC, were investigated by means of density functional theory computations. The Al₁₂X clusters can tightly bind to DNA bases and BPs to form stable complexes with negative binding Gibbs free energies at room temperature, and considerable charge transfers occur between the bases/BPs and the Al₁₂X clusters. These strong interactions, which are also expected for larger Al nanoparticles, may have potentially adverse impacts on the structure and stability of DNA and thus cause its dysfunction.     

AC-motif: a DNA motif containing adenine and cytosine repeat plays a role in gene regulation

I-motif or C4 is a four-stranded DNA structure with a protonated cytosine:cytosine base pair (C+:C) found in cytosine-rich sequences. We have found that oligodeoxynucleotides containing adenine and cytosine repeats form a stable secondary structure at a physiological pH with magnesium ion, which is similar to i-motif structure, and have named this structure ‘adenine:cytosine-motif (AC-motif)’. AC-motif contains C⁺:C base pairs intercalated with putative A⁺:C base pairs between protonated adenine and cytosine. By investigation of the AC-motif present in the CDKL3 promoter (AC-motif_CDKL3), one of AC-motifs found in the genome, we confirmed that AC-motif_CDKL3 has a key role in regulating CDKL3 gene expression in response to magnesium. This is further supported by confirming that genome-edited mutant cell lines, lacking the AC-motif formation, lost this regulation effect. Our results verify that adenine-cytosine repeats commonly present in the genome can form a stable non-canonical secondary structure with a non-Watson–Crick base pair and have regulatory roles in cells, which expand non-canonical DNA repertoires.     

Propeller-twisted adenine.thymine and guanine.cytosine base pairs tend to buckle and stagger in opposite directions.

Base pairs are propeller-twisted, buckled and staggered in DNA fragment crystals. These deformations were analyzed with isolated Watson-Crick base pairs using empirical potentials and buckle was found to almost linearly correlate with propeller. Interestingly, the thymine.adenine pair favours negative buckling for propellers mostly observed in DNA crystals while positive buckling is preferred by the cytosine.guanine pair. The propeller also induces opposite staggers in the adenine.thymine and guanine.cytosine base pairs.     

Perturbation of hydrogen bonds in the adenine ... thymine base pair by Na+, Mg2+, Ca2+ and NH4+ cations

Perturbation of the hydrogen bonds in the adenine ... thymine base pair by Na+, Mg2+, Ca2+ and NH4+ cations has been investigated by means of ab initio SCF calculations with the STO-3G basis set. The geometry of adenine...thymine, as well as those of the perturbed pairs were optimized. Approach of any cation to thymine at O6 leads to destabilization of the adenine...thymine pair; divalent cations (Mg2+, Ca2+) have a profound effect on the structure of the base pair. The approach of a cation to other available sites (thymine: O2, adenine N1 and N3) leads, on the other hand, to stabilization of the base pair. If a water molecule is placed between the cation and the base pair, the structure and stability of the base pair are changed only negligibly.     

NMR studies on an oligodeoxynucleotide containing 2-aminopurine opposite adenine

A heteroduplex containing the mismatch 2-aminopurine (AP)-adenine has been synthesized and studied by proton NMR. The mismatch was incorporated into the sequence d[CGG(AP)GGC].d-(GCCACCG). One-dimensional nuclear Overhauser effect measurements in H2O and two-dimensional nuclear Overhauser effect spectra in D2O show AP.A base pairs in a wobble structure in which both bases are in the anti conformation. The adenine is stacked well in the helix, but the helix twist between the adenine and neighboring cytosine in the 3' direction is unusually small. As a result, the aminopurine on the opposite strand is somewhat pushed out of the helix. From the measurements of the imino proton line widths, the two adjacent G.C base pairs are not found to be significantly destabilized by the presence of the purine-purine wobble pair.     

Structure of the deoxytetranucleotide d-pApTpApT and a sequence-dependent model for poly(dA-dT)

The x-ray structure of the hydrated ammonium salt of the deoxytetranucleotide d-pApTpApT was determined by Patterson and direct methods at a resolution of 1 Å. The crystal structure contains right-handed double-helical segments formed by complementary Watson-Crick-type hydrogen bonding between the adenine and thymine bases of neighboring molecules. The minihelix contains two base pairs. The chains are antiparallel. The A-T and T-A sequences have different phosphodiester conformations. The deoxyribose-pucker and the sugar–base orientation alternate along the chain depending on the nature of the base (3′-endo for purine and 2′-endo for pyrimidine). The extended structure is stabilized by base–base, base–sugar, and hydrogen-bond interactions. The minihelix of two base pairs provides starting coordinates for model-building studies of the dA-dT polymer. A B-DNA-type polymer structure is described, which has sequence-dependent alternations of both the deoxyribose pucker and the phosphate diester bridge conformation. Such sequence-dependent DNA structures, if present locally in regions such as operator sequences, could facilitate sequence-specific interactions. The crystal study also suggests possible geometrical parameters for the replication fork.     

Ab initio molecular orbital study of the mispairing ability of a nucleotide base analogue, N4-aminocytosine

The intrinsic properties of N4-aminocytosine, a base analogue of cytosine, are analyzed by an ab initio molecular orbital method. Relative stabilities of four possible isomeric structures of N4-aminocytosine are shown. The more stable isomer has the smaller dipole moment, so the relative stabilities of the isomers in solutions are subject to solvent polarity. The mutagenicity of this base analogue must arise because it can behave like either cytosine or thymine. It can form a guanine-cytosine-like base pair more easily than cytosine, and an adenine-thymine-like base pair less easily than thymine.     

Synthesis and properties of oligonucleotides containing 2'-deoxynebularine and 2'-deoxyxanthosine.

The preparation of synthetic oligonucleotides containing 2'-deoxynebularine (dN) and 2'-deoxyxanthosine (dX) is described. The thermal stabilities of duplexes containing dX, dN, and 2'-deoxyinosine (dI) base-paired with the four natural bases have been measured. Xanthine base pairs have stabilities at pH 5.5 that are similar to those of dI-containing duplexes at neutral pH. When xanthine is paired with adenine or cytosine an unusual stabilization of the duplex structure is observed at acid pH. Incorporation of base mispairs opposite template xanthine sites were measured using Drosophila DNA polymerase alpha. The relative nucleoside incorporation rates are in the order: T greater than C much greater than A approximately equal to G. These rates do not correlate with relative thermodynamic stabilities of base mispairs with xanthine obtained from Tm measurements: T greater than G greater than A approximately equal to C. We suggest that DNA polymerase misinsertion rates are greatest when the base mispair can be formed in accordance with Watson-Crick as opposed to other base pairing geometries even though other geometries, e.g. wobble, may result in a more stable final DNA product.     

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

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

Structure of the 2-aminopurine-cytosine base pair formed in the polymerase active site of the RB69 Y567A-DNA polymerase

The adenine base analogue 2-aminopurine (2AP) is a potent base substitution mutagen in prokaryotes because of its enhanceed ability to form a mutagenic base pair with an incoming dCTP. Despite more than 50 years of research, the structure of the 2AP-C base pair remains unclear. We report the structure of the 2AP-dCTP base pair formed within the polymerase active site of the RB69 Y567A-DNA polymerase. A modified wobble 2AP-C base pair was detected with one H-bond between N1 of 2AP and a proton from the C4 amino group of cytosine and an apparent bifurcated H-bond between a proton on the 2-amino group of 2-aminopurine and the ring N3 and O2 atoms of cytosine. Interestingly, a primer-terminal region rich in AT base pairs, compared to GC base pairs, facilitated dCTP binding opposite template 2AP. We propose that the increased flexibility of the nucleotide binding pocket formed in the Y567A-DNA polymerase and increased "breathing" at the primer-terminal junction of A+T-rich DNA facilitate dCTP binding opposite template 2AP. Thus, interactions between DNA polymerase residues with a dynamic primer-terminal junction play a role in determining base selectivity within the polymerase active site of RB69 DNA polymerase.     

Complex of cis-(NH3)2PtCl2 with guanylyl(3′-5′)-cytidine

The predominant complex formed by the reaction of cis-(NH₃)₂PtCl₂ and guanylyl(3′-5′)cytidine has been isolated. The molar ratio of the binding of cis-(NH₃)₂PtCl₂ to guanylyl(3′-5′)-cytidine is 1:2. The values of proton dissociation constant due to guanine and cytosine bases provide useful information for determining the binding site of the isolated complex. In addition, nmr and ir spectral data were used to determine the binding site. cis-(NH₃)₂PtCl₂ coordinates to guanylyl(3′-5′)cytidine through N(7) position of the guanine base, but cytosine base does not participate in the binding to cis-(NH₃)₂Pt²⁺. Interbase crosslink has not been detected. The binding specificity of cis-(NH₃)₂PtCl₂ to guanine base is discussed.     

The Effect of Methylation of the 6 Oxygen of Guanine on the Structure and Stability of Double Helical DNA

Abstract

The effect of methylation of the 0–6 position of guanine in short segments of double helical DNA has been investigated by molecular mechanical simulations on the sequences d(CGCGCG)₂, d(CGC+AFs-OMG+AF0-CG)₂, d(CGT+AFs-OMG+AF0-CG)₂, d(CGC+AFs-OMC+AF0-CG/(CGCGCG), d(CGC+AFs-OMG+AF0-CG/d(CGTGCG), d(CGCGAATTCGCG)₂ and d(CGCGAATTC+AFs-OMG+AF0-CG)₂. Guanines methylated at the 0–6 position are found to form hydrogen bonds of roughly equal strength to cytosine and thymine. The optimum structure of these modified base pairs are not dramatically different from normal GC pairs, but both involve some bifurcation of the proton donors of cytosine (4NH₂) or thymine (3NH) between the guanine N3 and O6 groups.     

Assignment of the proton Nmr chemical shifts of the T-N3H and G-N1H proton resonances in isolated AT and GC Watson-Crick base pairs in double-stranded deoxy oligonucleotides in aqueous solution

The 300-MHz proton nmr spectra (between 11 and 14 ppm) of a series of double-stranded deoxy oligonucleotides of known sequence have been recorded in H₂O solution. These resonances have been assigned to the G? N₁H and T? N₃H protons of specific base pairs from an evaluation of the temperature dependence of the ring NH linewidths and from the selective ring NH chemical shift changes on actinomycin-D binding. The deoxy oligonucleotides exist predominantly in the DNA-B conformation as evaluated from antibiotic binding studies. Ring-current calculations have been utilized to evaluate the up-field shifts of the G? N₁H and T? N₃H protons in Watson-Crick base pairs due to the ring currents from the pyrimidine and purine rings of nearest neighbor base pairs in regular DNA-B- and RNA-A-type helices. The perturbations on these up-field ring-current contributions that arise from twisting and tilting a base pair adjacent to the ring NH under study have been evaluated and found to change the calculated chemical shift by ±0.6 ppm for twist and tilt distortions of <30°C in a single adjacent base pair. A knowledge of the experimentally assigned ring NH chemical shifts of specific base pairs in known sequences of double-stranded deoxy oligonucleotides coupled with the ring-current tables for the DNA-B helical structure permit the assignment of 13.6 ± 0.1 ppm and 14.6 ± 0.2 ppm for the G? N₁H proton of an isolated GC base pair and the T? N₃H proton of an isolated AT base pair, respectively.     

Insights into the DNA stabilizing contributions of a bicyclic cytosine analogue: crystal structures of DNA duplexes containing 7,8-dihydropyrido [2,3-d]pyrimidin-2-one

The incorporation of the bicyclic cytosine analogue 7,8-dihydropyrido[2,3-d]pyrimidin-2-one (X) into DNA duplexes results in a significant enhancement of their stability (3–4 K per modification). To establish the effects of X on the local hydrogen-bonding and base stacking interactions and the overall DNA conformation, and to obtain insights into the correlation between the structure and stability of X-containing DNA duplexes, the crystal structures of [d(CGCGAATT-X-GCG)]₂ and [d(CGCGAAT-X-CGCG)]₂ have been determined at 1.9–2.9 Å resolutions. In all of the structures, the analogue X base pairs with the purine bases on the opposite strands through Watson–Crick and/or wobble type hydrogen bonds. The additional ring of the X base is stacked on the thymine bases at the 5′-side and overall exhibits greatly enhanced stacking interactions suggesting that this is a major contribution to duplex stabilization.