Conformational Analysis of Canonical 2-Deoxyribonucleotides. 1. Pyrimidine Nucleotides

Abstract

The molecular structure and relative stability of north and south conformers of 2′-deoxyribonucleotides containing pyrimidine nucleic acid bases (2′-deoxythymidilic (pdT), 2′- deoxycytidilic (pdC) acids and their mono- and dianions) have been obtained and analyzed at the DFT/B3LYP level using the standard 6–31G(d) basis set. We have revealed that, when the nucleobase moiety is incorporated into the nucleotides, it maintains a nonplanar and nonrigid conformation due to out-of-plane deformation of the amino group and pyrimidine ring. It has been demonstrated that an increase of negative charge of the phosphate group results in increase of amino group pyramidalization, discrimination between conformers with syn and anti orientation of base with respect to sugar, strengthening of intramolecular C-H…O hydrogen bonds leading to deformation and fixation of geometry of nucleotides, and weakening of phosphodiester bond. These results allow to make suggestions about sources of twist and buckle deformations of base pairs, mechanisms of repaire of DNA via change of base orientation, and conditions for breakage of the P-O bonds during hydrolysis.     

Conformational analysis of canonical 2-deoxyribonucleotides. 2. Purine nucleotides

The molecular structure of different conformers of isolated canonical purine 2'-deoxyribonucleotides 2-deoxyadenosine-5'-phosphate (pdA) and 2'-deoxyguanosine-5'-phosphate (pdG) was optimized using the B3LYP/6-31G(d) method. The results of the calculations reveal that the geometrical parameters and relative stability of the conformers significantly depend on the nature of the nucleobase, its orientation, the conformation of the furanose ring, the charge of the phosphate group and the character of the intramolecular hydrogen bonds. Analysis of the electron density distribution in purine nucleotides reveals the existence of a number of intramolecular hydrogen bonds. In general, the south conformer has a lower energy at the anti orientation of the base, and both conformers occur as the most stable for the syn orientation of the nucleobases. The results of the calculations reveal that the geometry and relative energy of the conformers of purine DNTs may be easily tuned by the charge of the phosphate group.     

Dependence of Deformability of Geometries and Characteristics of Intramolecular Hydrogen Bonds in Canonical 2′-Deoxyribonucleotides on DNA Conformations

Abstract

The molecular structure and deformability (with respect to average geometry) of methyl ethers of canonical 2′-deoxyribonucleotides thymidine-5′-phosphate (mTMP), 2-deoxycytidine-5′- phosphate (mCMP), 2-deoxyadenosine-5′-phosphate (mAMP) and 2′-deoxyguanosine-5′- phosphate (mGMP) in different types of DNA have been calculated using B3LYP/cc-pvdz method. Comparison of energy at equilibrium conformations of nucleotides and conformations with torsion angles of backbone fixed to average values for different types of DNA reveals that incorporation of nucleotides to A-DNA macromolecules requires the minimum amount of deformation energy. Therefore, this type of DNA should be the least strained from viewpoint of intramolecular deformations of monomers. Modeling of environmental effects within the PCM approach reveals that the immersion of nucleotides in polar medium results in significant decrease of energy differences between anti conformers of all DNTs and syn conformers of mGMP This also leads to reduction by almost a half nucleotides' deformation energy facilitating formation of DNA macromolecule. Change of DNTs conformation causes switch between different types of intramolecular H bonds. Every type of DNA possesses unique set of intramolar hydrogen bonds in nucleotides.     

Configurational effects on conformational properties of cyclic nucleotides. I. Theoretical calculations of conformer preferences in α-nucleoside 3′,5′ cyclic monophosphates

A comparative study has been made of the configurational effects on the conformational properties of α- and β-anomers of purine and pyrimidine nucleoside 3′,5′,-cyclic monophosphates and their 2′-arabino epimers. Correlation between orientation of the base and the 2′-hydroxyl group have been studied theoretically using the PCILO (Perturbative Configuration Interaction using Localized Orbitals) method. The effect of change in ribose puckering on the base-hydroxyl interaction has also been studied. The result show that steric repulsions and stabilizing effects of intramolecular hydrogen bonding between the base and the 2′-hydroxyl (OH) group are of major importance in determining configurations of α-anomers and 2′-arabino-β-epimers. For example, hydrogen bonding between the 2′-hydroxyl group and polar centers on the base ring is clearly implicated as a determinant of syn-anti preferences of the purine (adenine) or pyrimidine (uracil) bases in α-nucleoside 3′,5′-cyclic monophosphates. Moreover, barrier heights for interconversion between conformers are sensitive to ribose pucker and 2′-OH orientations. The result clearly show that a change in ribose-ring pucker plays an essential role in relieving repulsive interaction between the base and the 2′-hydroxyl group. Thus a C2′-exo-C3′-endo (₂T³) pucker is favored for α-anomers in contrast with the C4′-exo-C3′-endo (₄T³) from found in β-compounds.     

A proton-relaxation-time study of the conformation of some purine and pyrimidine 5'-nucleotides in aqueous solution.

The measurement of proton relaxation time T1 in a series of purine and pyrimidine 5'-nucleotides has been carried out to investigate their conformation in dilute neutral aqueous solutions. The interpretation of relaxation data has been performed with the help of computer calculations taking into account the different conformers of the ribose ring and of the exocyclic group. It has been found that all nucleotides under study show nearly the same preferential orientations of the base defined by gamma0 = 70 degrees +/- 10 degrees in the syn range. A more elaborate treatment, using an angular distribution derived from calculated potentials on 5" -GMP gives theoretical relaxation times in close agreement with experimental ones.     

Crystal structures of aspartate carbamoyltransferase ligated with phosphonoacetamide, malonate, and CTP or ATP at 2.8-A resolution and neutral pH

The R-state structures of the ATP and CTP complexes of aspartate carbamoyltransferase ligated with phosphonoacetamide and malonate have been determined at 2.8-A resolution and neutral pH. These structures were solved by the method of molecular replacement and were refined to crystallographic residuals between 0.167 and 0.182. The triphosphate, the ribose, and the purine and pyrimidine moieties of ATP and CTP interact with similar regions of the allosteric domain of the regulatory dimer. ATP and CTP relatively increase and decrease the size of the allosteric site in the vicinity of the base, respectively. For both CTP and ATP at pH 7, the gamma-phosphates are bound to His20 and are also near Lys94, while the alpha-phosphates interact exclusively with Lys94. The 2'-hydroxyls of both CTP and ATP are near the amino group of Lys60. The pyrimidine ring of CTP makes specific hydrogen bonds at the allosteric site: the NH2 group donates hydrogen bonds to the main-chain carbonyls of Ile12 and Tyr89 and the pyrimidine ring carbonyl oxygen accepts a hydrogen bond from the amino group of Lys60; the nitrogen at position 3 in the pyrimidine ring is hydrogen bonded to a main-chain NH group of Ile12. The purine ring of ATP also makes numerous interactions with residues at the allosteric site: the purine NH2 (analogous to the amino group of CTP) donates a hydrogen bond to the main-chain carbonyl oxygen of Ile12, the N3 nitrogen interacts with the amino group of Lys60, and the N1 nitrogen hydrogen bonds to the NH group of Ile12. The binding of CTP and ATP to the allosteric site in the presence of phosphonoacetamide and malonate does not dramatically alter the structure of the allosteric binding site or of the allosteric domain. Nonetheless, in the CTP-ligated structure, the average separation between the catalytic trimers decreases by approximately 0.5 A, indicating a small shift of the quaternary structure toward the T state. In the CTP- and ATP-ligated R-state structures, the binding and occupancy of phosphonoacetamide and malonate are similar and the structures of the active sites are similar at the current resolution of 2.8 A.     

Effects of the minor groove pyrimidine nucleobase functional groups on the stability of duplex DNA: the impact of uncompensated minor groove amino groups

DNA sequences containing four types of analog nucleosides are described. All four are pyridine derivatives constructed as C-nucleosides so that they mimic the pyrimidine derivatives 2'-deoxyuridine, thymidine or 2'-deoxycytidine, but in all cases the analogs lack the corresponding O2-carbonyls that in duplex DNA are located in the minor groove. In place of the O2-carbonyl is a hydrogen atom, a polar fluorine atom, or a nonpolar methyl group. The described C-nucleosides have native-like bidentate Watson-Crick hydrogen-bonding faces and can form essentially normal W-C base pairs of varying stability with A or G. In each modified base pair, two inter-residue hydrogen bonds should be present. In spite of a common number of interstrand hydrogen bonds, the thermodynamic stabilities of the prepared duplexes, each containing two analog base pairs, vary dramatically. Most notably, base pairs containing uncompensated purine amino groups (those lacking a hydrogen-bonding partner) in the minor groove exhibit the most dramatic reductions in thermodynamic stability. Removal of such uncompensated amino groups results in increased duplex stability. Base pairs containing fluorine in the minor groove positioned adjacent to an amino group seem to enhance duplex stability marginally (relative to --H or --CH(3)), but there is little evidence to suggest that fluorine is an effective hydrogen-bonding partner in these systems. The presence of minor groove methyl groups results in the least stable duplexes in each series of sequences.     

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.     

Protein and drug interactions in the minor groove of DNA

Interactions between proteins, drugs, water and B-DNA minor groove have been analyzed in crystal structures of 60 protein–DNA and 14 drug–DNA complexes. It was found that only purine N3, pyrimidine O2, guanine N2 and deoxyribose O4′ are involved in the interactions, and that contacts to N3 and O2 are most frequent and more polar than contacts to O4′. Many protein contacts are mediated by water, possibly to increase the DNA effective surface. Fewer water-mediated contacts are observed in drug complexes. The distributions of ligands around N3 are significantly more compact than around O2, and distributions of water molecules are the most compact. Distributions around O4′ are more diffuse than for the base atoms but most distributions still have just one binding site. Ligands bind to N3 and O2 atoms in analogous positions, and simultaneous binding to N3 and N2 in guanines is extremely rare. Contacts with two consecutive nucleotides are much more frequent than base–sugar contacts within one nucleotide. The probable reason for this is the large energy of deformation of hydrogen bonds for the one nucleotide motif. Contacts of Arg, the most frequent amino acid ligand, are stereochemically indistinguishable from the binding of the remaining amino acids except asparagine (Asn) and phenylalanine (Phe). Asn and Phe bind in distinct ways, mostly to a deformed DNA, as in the complexes of TATA-box binding proteins. DNA deformation concentrates on dinucleotide regions with a distinct deformation of the δ and backbone torsion angles for the Asn and δ, , ζ and χ for the Phe-contacted regions.     

Turnover of carnitine by rat tissues.

A small release of Pi from a diphenylamine-formic acid digest of DNA was detected after elimination of interpurine phosphodiester bonds was complete. Minor components in the DNA digest were identified as pyrimidine oligonucleotides which had lost one terminal phosphate. Isolated pyrimidine tracts released Pi on redigestion with the formic acid-diphenylamine reagent in amounts that increased with the number of nucleotides in the oligonucleotide taken. The oligonucleotides were also partially degraded by the formic acid-diphenylamine reagent and the degradation (2-3% of phosphodiester bonds between consecutive nucleotides) was almost independent of chain length. The cleavage was random with no preference for a phosphodiester bond flanked by particular nucleosides. This minor lack of specificity in the formic acid-diphenylamine-catalysed degradation of DNA can, however, account for the low recoveries of long pyrimidine tracts previously reported. Any analysis of pyrimidine tracts in a DNA molecule should make some correction for this small degree of degradation if exact assignments of the numbers of pyrimidine tracts are to be made.     

The secondary structure of nucleotidyl-(5' bound to N)-amino acids containing different heterocyclic bases and amino acids]

Nucleotidyl-(5' leads to N)-amino acids containing different heterocycle bases: adenine, guanine, hypoxanthine, cytosine, uracyl, and aromatic amino acids: phenylalanine, tyrosine and tryptophan, have been investigated by proton magnetic resonance and circular dichroism. For all the compounds studied folded conformation have been shown stabilized by hydrophobic interaction in aqueous solution. The comparison of the results of the studied nucleotidyl-(5' leads to N)-amino acids unable us to build four secondary structure types in these very compounds. Phenylalanine and tyrosine derivatives of purine nucleotides can be regarded as the first type, tryptophan derivatives of purine nucleotides as the second type, phenylalanine and tyrosine derivatives of pyrimidine nucleotides as the third type and tryptophan derivatives of pyrimidine nucleotides as the fourth type. For each group of these compounds conformational models have been built. In all these compounds the anti-conformation has been proved to exist.     

Molecular mechanics calculations on a triple stranded DNA involving C+.G-T and T.A+-C mismatched base triples

We have carried out molecular modeling of a triple stranded pyrimidine(Y). purine(R): pyrimidine(Y) (where ':' refers to Watson-Crick and '.' to Hoogsteen bonding) DNA, formed by a homopurine (d-TGAGGAAAGAAGGT) and homo-pyrimidine (d-CTCCTTTCTTCC). Molecular mechanics calculations using NMR constraints have provided a detailed three dimensional structure of the triplex. The entire stretches of purine and the pyrimidine nucleotides have a conformation close to B-DNA. The three strands are held by the canonical C+.G:C and T.A:T hydrogen bonds. The structure also contains two mismatch C+.G-T and T.A+-C base triples which have been characterized for the first time. In the A+-C base-pair of the T.A+-C triple, both hydrogen donors are situated on the purine (A+(1N) and A+(6N)). We observe a unique hydrogen bonding interaction scheme in case of C+.G-T where one acceptor, G(60), is bonded to three donors (C+(3NH), C+(4NH2) and T(3NH)). Though the C+.G-T base triple is less stable than C+.G:C, it is significantly more stable than T.A:T. On the other hand, T.A+-C is as stable as the T.A:T base triad.     

The dimensions and shapes of the furanose rings in nucleic acids

A survey was made of the geometry of furanose rings in beta-nucleotides and beta-nucleosides (as monomers related to nucleic acids) for which structures have been determined by X-ray crystallography. Mean values, and estimated standard deviations from them, were calculated for bond-lengths, bond-angles and conformation-angles. For parameters with values dependent on ring-puckering, separate calculations were made for each ring type. (The rings are puckered in one of three conformations: C-2- or C-3-endo or C-3-exo; C-2-exo has not been observed.) The results were used to compute standard furanose rings with C-2-endo, C-3-endo and C-3-exo conformations for use in nucleic acid molecular model-building. The survey also showed that the only other conformation-angle in nucleotides dependent on the furanose ring conformation corresponds to the relative orientation of the purine (or pyrimidine) base and the ring.     

Binding of the antitumor drug nogalamycin and its derivatives to DNA: structural comparison

The three-dimensional molecular structures of the complexes between a novel antitumor drug nogalamycin and its derivative U-58872 with a modified DNA hexamer d[m5CGT(pS)Am5CG] have been determined at 1.7- and 1.8-A resolution, respectively, by X-ray diffraction analyses. Both structures (in space group P6(1)) have been refined with constrained refinement procedure to final R factors of 0.208 (3386 reflections) and 0.196 (2143 reflections). In both complexes, two nogalamycins bind to the DNA hexamer double helix in a 2:1 ratio with the elongated aglycon chromophore intercalated between the CpG steps at both ends of the helix. The aglycon chromophore spans across the GC Watson-Crick base pairs with its nogalose lying in the minor groove and the aminoglucose lying in the major groove of the distorted B-DNA double helix. Most of the sugars remain in the C2'-endo pucker family, except three deoxycytidine residues (terminal C1, C7, and internal C5). All nucleotides are in the anti conformation. Specific hydrogen bonds are found in the complex between the drug and guanine-cytosine bases in both grooves of the helix. One hydroxyl group of the aminoglucose donates a hydrogen bond to the N7 of guanine, while the other receives a hydrogen bond from the N4 amino group of cytosine. The orientation of these two hydrogen bonds suggests that nogalamycin prefers a GC base pair with its aglycon chromophore intercalating at the 5'-side of a guanine (between NpG), or at the 3'-side of a cytosine (between CpN) with the sugars pointing toward the GC base pair. The binding of nogalamycin to DNA requires that the base pairs in DNA open up transiently to allow the bulky sugars to go through, suggesting that nogalamycin prefers GC sequences embedded in a stretch of AT sequences.     

Structure of a beta-alanine-linked polyamide bound to a full helical turn of purine tract DNA in the 1:1 motif

Polyamides composed of N-methylpyrrole (Py), N-methylimidazole (Im) and N-methylhydroxypyrrole (Hp) amino acids linked by beta-alanine (beta) bind the minor groove of DNA in 1:1 and 2:1 ligand to DNA stoichiometries. Although the energetics and structure of the 2:1 complex has been explored extensively, there is remarkably less understood about 1:1 recognition beyond the initial studies on netropsin and distamycin. We present here the 1:1 solution structure of ImPy-beta-Im-beta-ImPy-beta-Dp bound in a single orientation to its match site within the DNA duplex 5'-CCAAAGAGAAGCG-3'.5'-CGCTTCTCTTTGG-3' (match site in bold), as determined by 2D (1)H NMR methods. The representative ensemble of 12 conformers has no distance constraint violations greater than 0.13 A and a pairwise RMSD over the binding site of 0.80 A. Intermolecular NOEs place the polyamide deep inside the minor groove, and oriented N-C with the 3'-5' direction of the purine-rich strand. Analysis of the high-resolution structure reveals the ligand bound 1:1 completely within the minor groove for a full turn of the DNA helix. The DNA is B-form (average rise=3.3 A, twist=38 degrees ) with a narrow minor groove closing down to 3.0-4.5 A in the binding site. The ligand and DNA are aligned in register, with each polyamide NH group forming bifurcated hydrogen bonds of similar length to purine N3 and pyrimidine O2 atoms on the floor of the minor groove. Each imidazole group is hydrogen bonded via its N3 atom to its proximal guanine's exocyclic amino group. The important roles of beta-alanine and imidazole for 1:1 binding are discussed.     

The reactivity of lactyl-oxythiamin implies the role of the amino-pyrimidine in thiamin catalyzed decarboxylation

It has previously been established that the deprotonated amino substituent of the pyrimidine of thiamin diphosphate (ThDP) acts as an internal base to accept the C2H of the thiazolium in ThDP-dependent enzymes. The amino group has also been implicated in assisting the departure of the aldehydic product formed after loss of CO₂ from ketoacid substrates. However, the potential role for the pyrimidine amino group in the key decarboxylation step has not been assessed. Oxythiamin contains a hydroxyl group in place of the pyrimidine amino group in thiamin, providing a basis for comparison of reactivity. Lactyl-oxythiamin (LOTh), the conjugate of pyruvic acid and oxythiamin was prepared by condensation of ethyl pyruvate and hydroxyl-protected oxythiamin followed by deprotection and acidic hydrolysis of the ethyl ester. The rate constants observed for the decarboxylation of LOTh in neutral and acidic solutions are about four times smaller than those for the corresponding compound that contains the amino group, lactylthiamin. The difference in reactivity is consistent with the amino group’s participation in facilitating the decarboxylation step by allowing a competitive addition pathway that produces bicarbonate and has implications for the corresponding enzymic reaction.     

Thermodynamics of 2'-ribose substitutions in UUCG tetraloops

The ribose 2'-hydroxyl group confers upon RNA many unique molecular properties. To better appreciate its contribution to structure and stability and to monitor how substitutions of the 2' hydroxyl can alter an RNA molecule, each loop pyrimidine ribonucleotide in the UUCG tetraloop was substituted with a nucleotide containing either a fluorine (2'-F), hydrogen (2'-H), amino (2'-NH2), or methoxy (2'-OCH3) group, in the context of both the C:G and G:C loop-closing base pair. The thermodynamic parameters of these tetraloop variants have been determined and NMR experiments used to monitor the structural changes resulting from the substitutions. The modified riboses are better tolerated in the G[UUCG]C tetraloop, which may be due to its increased loop flexibility relative to the C[UUCG]G loop. Even for these simple substitutions, the free-energy change reflects a complex interplay of hydrogen bonding, solvation effects, and intrinsic pucker preferences of the nucleotides.     

Structure of a G.T.A triplet in an intramolecular DNA triplex.

A 32-base DNA oligonucleotide has been studied by one- and two-dimensional 1H NMR spectroscopy and is shown to form a stable, pyr.pur.pyr, intramolecular triple helical structure, with a four C loop and a TATA loop connecting the Watson-Crick- and Hoogsteen-paired strands, respectively. This triplex contains five T.A.T base triplets, two C+.G.C base triplets, and an unusual G.T.A base triplet which disrupts the pyr.pur.pyr motif. The G.T.A triplet consists of a Watson-Crick T.A base pair, with the T situated in the "purine strand" and the A situated in the "pyrimidine strand" and a G situated in the Hoogsteen-base-paired "pyrimidine strand" hydrogen bonded to the T. The base-pairing structure of the G.T.A triplet has been investigated and has been found to involve a single hydrogen bond from the guanine amino group to the O4 carbonyl of the thymine, leaving the guanine imino proton free. The specific amino proton involved in the hydrogen bond is the H2(2) proton. This orients the guanine such that its sugar is near the thymine methyl group. The guanine sugar adopts an N-type (C3'-endo) sugar pucker in this triplet. The stability of the G.T.A triplet within pyr.pur.pyr triplexes is discussed.     

Pyrimidine-specific cleavage by an endoribonuclease of Saccharomyces cerevisiae.

An endoribonuclease with pyrimidine cleavage site specificity was isolated from Saccharomyces cerevisiae. The enzyme had a pH optimum of 6 to 7 and did not require a divalent cation. It was inhibited by 5 X 10(-5) M ethidium bromide, although it appeared to be single strand specific. The enzyme gave a limited cleavage of yeast mRNA and rRNA, yielding products that were terminated with pyrimidine nucleoside 2',3'-cyclic phosphate. The bonds between pyrimidine and A residues constituted more than 90% of the scission sites when the average product size was 50 nucleotides. Homopolyribonucleotides were cleaved poorly. Poly(A,U) was cleaved rapidly, and analysis of the products of poly(A,U) hydrolysis showed a very stringent cleavage of U-A bonds.     

Conformational sub-states in B-DNA.

Theoretical studies of the sequence-dependent conformation of B-DNA have been carried out using Jumna, a helicoidal co-ordinate minimization algorithm. The results obtained for a series of six oligomers with repetitive sequences show that, with the exception of the homopolymers (dA)n.(dT)n and (dG)n.(dC)n, all sequences can adopt a variety of conformations characterized by considerable changes in helicoidal parameters and also in sugar puckers which adopt C(2')-endo (falling into 2 classes) or, in the case of pyrimidine nucleotides, O(1')-endo forms. These studies lead to an improved understanding of the role of base sequence on DNA conformation and point to a number of interesting correlations between the various structural parameters describing the double helix.