Conformational studies of 13 trinucleoside bisphosphates by 360-MHz 1H-NMR spectroscopy. 1. Ribose protons

The ribose protons of 13 trinucleoside bisphosphates (trimers) were studied, using 360-MHz proton nuclear magnetic resonance spectroscopy. Complete assignments and analyses of the NMR signals of these protons were carried out by the methods of homonuclear decoupling and computer line-shape simulations. It was shown that the trinucleotides preferred the anti, 3' endo, gamma +, beta t and epsilon t/epsilon- conformations for the glycosidic torsions, the ribose rings, the C4'-C5' bonds, the C5'-O5' bonds, and the C3'-O3' bonds, respectively. It was also found that the trimers, especially those which had noticeable population of 'bulged' structures, did not necessarily have a higher population of these preferred local conformations than their component dimers. The overall conformations of the trinucleotides are classified into two categories. The conformations in the first category involve the nearest-neighbor interactions. Each dinucleotide moiety can assume one of the four stable conformations (I, I', II and III) or the open forms of dinucleoside monophosphates. However, due to steric hindrance, there are only four cases in which both dinucleotide moieties can assume one of the four stable conformations at the same time. These four combinations of conformations are I-I, I'-I', I-II and III-I', where the first Roman numeral represents the conformation of the NpN'p-moiety and the second one, that of the -pN'pN' moiety of the trimers. Among them, I-I and I'-I' are helical structures, capable of forming a double helix. The second category contains conformations with bulged structures which have the two dinucleotide moieties in open forms (i.e. no nearest-neighbor interactions) and the bases of the two terminal residues stacking on each other while the middle residue is bulged out. These bulged conformations may serve as structural models for frame-shift mutations.     

Single-stranded adenine-rich DNA and RNA retain structural characteristics of their respective double-stranded conformations and show directional differences in stacking pattern

The structural preorganization of isosequential ssDNA and ssRNA hexamers d/r(GAAAAC)(1) [J. Am. Chem. Soc. 2003, 125, 9948] have been investigated by NMR and molecular dynamics simulations. Analysis of the nuclear Overhauser effect spectrometry (NOESY) footprints in the aqueous solution has shown that there is a substantial population of ordered right-handed helical structure in both hexameric single-stranded DNA and RNA, which are reminiscent of their respective right-handed helical duplex form, despite the fact these single-stranded molecules are devoid of any intermolecular hydrogen bonds. The NMR-constrained molecular dynamics (1.5 ns) derived geometries of the adenine-adenine overlaps at each dinucleotide step of the hexameric ssDNA (1a) and ssRNA (1b) show that the relatively electron-rich imidazole stacks above the electron-deficient pyrimidine in 5' to 3' direction in ssDNA (1a) while, in contradistinction, the pyrimidine stacks above the imidazole in the 5' to 3' direction in ssRNA (1b). This also means that the pi-frame of the 5'-pyrimidine can interact with the relatively positively charged imino and amino protons in the 3' direction in ssRNA and in the 5' direction in ssDNA, thereby stabilizing the twist and slide observed in the stacked oligonucleotides. The differently preferred stacking geometries in ssDNA and ssRNA have direct physicochemical implications for self-assembly and pK(a) modulation by the nearest-neighbor interactions, as well as for the dangling-end stabilization effects and imino-proton reactivity.     

Stability of 3' double nucleotide overhangs that model the 3' ends of siRNA

Thermodynamic parameters are reported for duplex formation in 1 M NaCl for 16 RNA sequences, each containing a core tetramer duplex, GGCC, and a 3' overhang consisting of two bases. The results indicate additional double-helical stability is conferred by the double 3' terminal overhang relative to the single 3' terminal overhang. A nearest-neighbor analysis of the data indicates that the free energy contribution at 37 degrees C of the second base in the double 3' terminal overhang varies from 0 to 0.7 kcal/mol. The second base in the 3' double overhang can contribute nearly the same stability to a duplex as a base pair or a 3' dangling overhang. Stability contribution of a dangling base, two nucleotides removed from the 3' end of a duplex, is dependent upon both the identity of the base as well as that of the dangling base that it neighbors. A second dangling base only increases the stability of the duplex when it is neighboring a 3' purine dangling nucleotide. Furthermore, a second dangling pyrimidine provides a greater contribution to duplex stability than a purine. A nearest-neighbor model was developed to predict the influence of 3' double overhang on the stability of duplex formation. The model improves the prediction of free energy and melting temperature when tested against six sequences with different core duplexes.     

Exchange mechanisms for hydrogen bonding protons of cytidylic and guanylic acids.

The pH dependence of buffer catalysis of exchange of the C-4 amino protons of cyclic cytosine 2',3'-monophosphate (cCMP) and the N-1 proton of cyclic guanosine 2',3'-monophosphate (cGMP) conforms to an exchange mechanism, in which protonation of the nucleobases at C(N-3) AND G(N-7) establishes the important intermediates at neutral to acidic pH. Rate constants for transfer of the G(N-1) proton to H2O, OH-, phosphate, acetate, chloracetate, lactate, and cytosine (N-3) were obtained from 1H nuclear magnetic resonance line width measurements at 360 MHz and were used to estimate the pK or acidity of the exchange site in both the protonated and unprotonated nucleobase. These estimates reveal an increase in acidity of the G(N-1) site corresponding to 2 to 3 pK units as the G(N-7) site is protonated: At neutral pH the G(N-1) site of the protonated purine would be ionized (pK = 6.3). Determinations of phosphate, imidazole, and methylimidazole rate constants for transfer of the amino protons of cCMP provide a more approximate estimate of pK = 7 to 9 for the amino of the protonated pyrimidine. A comparison of the intrinsic amino acidity in the neutral and protonated cytosine is vitiated by the observation that OH- catalyzed exchange in the neutral base is not diffusion limited. This leads to the conclusion that protonation of the nucleobase effects a qualitative increase in the ability of the amino protons to form hydrogen bonds: from very poor in the neutral base to "normal" in the conjugate acid.     

A structural model for sequence-specific proflavin-DNA interactions during in vitro frameshift mutagenesis.

Molecular models describing intermediates that may lead to proflavin-induced 1 bp deletions during in vitro polymerization by E. coli DNA polymerase I Klenow fragment are proposed. The models provide structural explanations for the fact that the induced frameshifts always occur opposite template bases that are adjacent to 5' pyrimidines and are based on the underlying hypothesis that the deletions arise because the polymerase passes by a template base without copying it. Because the most frequent mutations are opposite Pu in the template sequence 5' Py Pu 3', a single-strand loop-out model was constructed for this sequence and proflavin was added, using structures found in crystalline oligonucleotides and their complexes with proflavin. The model seeks to rationalize the roles of the 5' pyrimidine and proflavin in facilitating the bypass. Four potential roles for proflavin in mutagenesis are described: 1) stacking on the looped-out base; 2) stacking on the base pair immediately preceding the site of mutation; 3) hydrogen bonding with the 5' pyrimidine; 4) hydrogen bonding with the phosphate backbone. These models point to the possibility that a number of proflavin-DNA interactions may be involved. In contrast, modeling does not suggest a role for classically intercalated proflavin in frameshift mutagenesis arising during in vitro DNA polymerization.     

Potential of mean force calculations of the stacking-unstacking process in single-stranded deoxyribodinucleoside monophosphates.

The free energy of the stacking-unstacking process of deoxyribodinucleoside monophosphates in aqueous solution has been investigated by potential of mean force calculations along a reaction coordinate, defined by the distance between the glycosidic nitrogen atoms of the bases. The stacking-unstacking process of a ribodinucleoside monophosphate was observed to be well characterized by this coordinate, which has the advantage that it allows for a dynamical backbone and flexible bases. All 16 naturally occurring DNA dimers composed of the adenine, cytosine, guanine, or thymine bases in both the 5' and the 3' positions were studied. From the free-energy profiles we observed the deepest minima for the stacked states of the purine-purine dimers, but good stacking was also observed for the purine-pyrimidine and pyrimidine-purine dimers. Substantial stacking ability was found for the dimers composed of a thymine base and a purine base and also for the deoxythymidylyl-3',5'-deoxythymidine dimer. Very poor stacking was observed for the dCpdC dimer. Conformational properties and solvent accessibility are discussed for the stacked and unstacked dimers. The potential of mean force profiles of the stacking-unstacking process for the DNA dimers are compared with the RNA dimers.     

Influence of uracil on the conformational behaviour of RNA oligonucleotides in solution

NMR studies were carried out on some alternating pyrimidine-purine sequences: the single-stranded tetramers CACA and UGUG and the self-complementary octamer CACAUGUG. Assignments, based upon COSY, homonuclear Hartmann-Hahn, and NOESY experiments, are given for the resonances of all base protons and of several sugar protons. Chemical shift vs temperature profiles were used to obtain thermodynamic parameters for the single-stranded stack in equilibrium with random coil and the duplex in equilibrium with random coil equilibria. The populations of N-type conformer of the ribose rings were estimated from the observed J1'2'. Comparisons with another alternating pyrimidine-purine sequence Um2(6)AUm2(6)A and with the deoxyribose counterparts d(CACA), d(TGTG) and d(CACATGTG) are given. Previous 1H-NMR investigations of Um2(6)AUm2(6)A revealed that the population of bulge-out structure diminishes compared to m2(6)AUm2(6)A due to the U(1)-m2(6)A(2) stacking interaction. In CACA a strong stacking proclivity (Tm = 310 K) together with a clear preference for N-type ribose is observed. However, the stacking interactions in UGUG are relatively less stable (Tm = 288 K) and a bias towards S-type sugar is present. Besides a small amount of stack, a significant contribution of bulge out structure is proposed for UGUG. We conclude that the nature of the pyrimidine base mainly determines the formation of bulge-out structures. The poor stacking properties of uracil now appear to be mainly responsible for this phenomenon. Comparison with the deoxyribose counterparts shows a reasonable agreement between the Tm values of CACA and d(CACA), whereas the Tm of UGUG (288 K) is much lower than the Tm of d(TGTG) (315 K). It is suggested that the absence of bulge-out structures in DNA purine-pyrimidine-purine sequences is related to the relatively strong stacking proclivity of dT residues compared to that of U residues. The Tm values (average 341 K) for the duplex in equilibrium with random coil transition obtained for each residue of CACAUGUG appear very similar. All ribose rings, except the G(8), adopt a pure N conformer in the duplex. This is taken to mean that the differences in conformational behaviour of the constituent tetramers disappear upon duplex formation.     

Stabilities of nearest-neighbor doublets in double-helical DNA determined by fitting calculated melting profiles to observed profiles

Melting profiles were calculated for restriction fragments of ?X174 and fd phage DNAs and compared with experimental profiles. The algorithm of Fixman and Freire was slightly modified so that a stability parameter was assigned not to a base pair but to each nearest-neighbor doublet. Stabilities of the 10 kinds of nearest-neighbor doublets were estimated by fitting the calculated profiles to the observed ones. Agreement of the calculated and observed profiles was much improved by this modification. The most interesting finding was that purine (3′-5′) pyrimidine stackings are much more stable than their respective reverses. The order of nearest-neighbor stabilities is in excellent agreement with that of negative stacking energies calculated by Rein and coworkers by a quantum-chemical method.     

Resonance assignments of non-exchangeable protons in B type DNA oligomers, an overview.

The chemical shifts of 1H resonances of non exchangeable protons (except H5', H5" and adenine H2) of over six hundred nucleotides have been collected. The influence which the base of the nucleotide itself as well as the bases on its 5' and 3' side exert on the chemical shifts of the various resonances has been investigated. Most of the resonances appear to be predominantly influenced by only one base. For H2', H2", H3', H4' and H6/H8 this is the base of the central nucleotide, for H5(C) and CH3(T) it is the one on the 5' side and for H1' it is the one on the 3' side. Chemical shift distribution profiles are presented which allow an estimation of the probability of finding a particular resonance at a particular position in the spectrum.     

Conformational determinants of tandem GU mismatches in RNA: insights from molecular dynamics simulations and quantum mechanical calculations

Structure and energetic properties of base pair mismatches in duplex RNA have been the focus of numerous investigations due to their role in many important biological functions. Such efforts have contributed to the development of models for secondary structure prediction of RNA, including the nearest-neighbor model. In RNA duplexes containing GU mismatches, 5'-GU-3' tandem mismatches have a different thermodynamic stability than 5'-UG-3' mismatches. In addition, 5'-GU-3' mismatches in some sequence contexts do not follow the nearest-neighbor model for stability. To characterize the underlying atomic forces that determine the structural and thermodynamic properties of GU tandem mismatches, molecular dynamics (MD) simulations were performed on a series of 5'-GU-3' and 5'-UG-3' duplexes in different sequence contexts. Overall, the MD-derived structural models agree well with experimental data, including local deviations in base step helicoidal parameters in the region of the GU mismatches and the model where duplex stability is associated with the pattern of GU hydrogen bonding. Further analysis of the simulations, validated by data from quantum mechanical calculations, suggests that the experimentally observed differences in thermodynamic stability are dominated by GG interstrand followed by GU intrastrand base stacking interactions that dictate the one versus two hydrogen bonding scenarios for the GU pairs. In addition, the inability of 5'-GU-3' mismatches in different sequence contexts to all fit into the nearest-neighbor model is indicated to be associated with interactions of the central four base pairs with the surrounding base pairs. The results emphasize the role of GG and GU stacking interactions on the structure and thermodynamics of GU mismatches in RNA.     

Sequence-specific recognition of deoxyribonucleic acid. Chemical synthesis and nuclear magnetic resonance assignment of the imino protons of lambda OR3 operator deoxyribonucleic acid

Using solid-phase phosphite triester methods, we have synthesized both strands of the phage lambda OR3 DNA sequence, reannealed them, and studied the native operator duplex by high-resolution NMR at 500 MHz. At 7 degrees C the imino protons of the two terminal base pairs at each end have disappeared from the spectrum by exchange broadening. The 13 detectable imino resonances have been assigned to their respective base pairs in the duplex by using sequential nearest-neighbor NOE connectivity methods described previously. In cases where two imino protons overlap in the spectrum, spin diffusion was used to drive the cross-saturation further afield in order to produce second-order next-nearest-neighbor effects. The results show that the imino connectivity method can be used to unambiguously assign the imino proton spectrum of operator DNAs containing one to two full turns of the helix.     

Conformation and interaction of short nucleic acid double-stranded helices. I. Proton magnetic resonance studies on the nonexchangeable protons of ribosyl ApApGpCpUpU.

1H nuclear magnetic resonance (NMR) spectra of a self-complementary ribosyl hexanucleotide, A2GCU2, are investigated as a function of temperature and ionic strength in D2O. Seventeen nonexchangeable base and ribose-H1' resonances are resolved, and unequivocally assigned by a systematic comparison with the spectra of a series of oligonucleotide fragments of the A2GCU2 sequence varying in chain length from 2 to 5. Changes in the chemical shifts of the 17 protons from the hexamer as well as the six H1'-H2' coupling constants are followed throughout a thermally induced helix-coil transition. These sigma vs. T and J vs. T (degrees C) profiles indicate that the transition is not totally cooperative and that substantial populations of partially bonded structures must exist at intermediate temperatures, with the central G-C region being most stable. Transitions in chemical shift for protons in the same base pair exhibit considerable differences in their Tm values as the data reflect both thermodynamic and local magnetic field effects in the structural transition, which are not readily separable. However, an average of the Tm values agrees well with the value predicted from studies of the thermally induced transition made by optical methods. The values of J1'-2' for all six residues become very small (less than 1.5 Hz) at low temperatures indicating that C3'-endo is the most heavily populated furanose conformation in the helix. The sigma values of protons in the duplex were compared with those calculated from the ring current magnetic anisotropies of nearest and next-nearest neighboring bases using the geometrical parameters of the A'-RNA and B-DNA models. The sigma values of the base protons in the duplex calculated assuming the A'-RNA geometry agree (+/- approximately 0.1 ppm) with the observed values much more accurately than those calculated on the basis of B-DNA geometry. The measured sigma values of the H1' are not accurately predicted from either model. The synthesis of 35 mg of A2GCU2 using primer-dependent polynucleotide phosphorylase is described in detail with extensive discussion in the microfilm edition.     

Structure-related properties of the mutagenic lesion 6-O-methylguanine in DNA.

Chemical probing of the structures of a few very similar 30 base-pair duplexes containing a 6-O-methylguanine (meG) residue at the 16th position reveals that the modified base simultaneously perturbs the helical structure in two ways; it preferentially unstacks the 3' neighbouring base residue (thymine in this study) on the same strand and it unstacks the pyrimidine to which it is base-paired. Depending on its neighbouring 5' base residue and the base-pairing pyrimidine, this perturbation can extend to a few base-pairs in both 3' and 5' directions from the abnormal base-pair. These perturbations can be detected by cleavage at the site for the restriction enzyme MaeII. The unstaking of the C in the meG.C and A.C base-pairs may explain the de novo methylation of these helices by the human DNA-(cytosine-5-)methyltransferase. Interestingly, the kinetics of repair of the 6-O-methylguanine-containing dinucleotides by the cloned human methylguanine methyltransferase appears to be largely determined by the strength of the stacking interaction between the 6-O-methylguanine and the 5' neighbouring base.     

Utilization of DNA photolyase, pyrimidine dimer endonucleases, and alkali hydrolysis in the analysis of aberrant ABC excinuclease incisions adjacent to UV-induced DNA photoproducts.

ABC excinuclease of Escherichia coli removes 6-4 photoproducts and pyrimidine dimers from DNA by making two single strand incisions, one 8 phosphodiester bonds 5' and another 4 or 5 phosphodiester bonds 3' to the lesion. We describe in this communication a method, which utilizes DNA photolyase from E. coli, pyrimidine dimer endonucleases from M. luteus and bacteriophage T4, and alkali hydrolysis, for analyzing the ABC excinuclease incision pattern corresponding to each of these photoproducts in a DNA fragment. On occasion, ABC excinuclease does not incise DNA exclusively 8 phosphodiester bonds 5' or 4 or 5 phosphodiester bonds 3' to the photoproduct. Both the nature of the adduct (6-4 photoproduct or pyrimidine dimer) and the sequence of neighboring nucleotides influence the incision pattern of ABC excinuclease. We show directly that photolyase stimulates the removal of pyrimidine dimers (but not 6-4 photoproducts) by the excinuclease. Also, photolyase does not repair CC pyrimidine dimers efficiently while it does repair TT or TC pyrimidine dimers.     

Enzymatic analysis of isomeric trithymidylates containing ultraviolet light-induced cyclobutane pyrimidine dimers. II. Phosphorylation by phage T4 polynucleotide kinase

Phage T4 polynucleotide kinase (EC 2.7.1.78) proved incapable of catalyzing the phosphorylation of thymidylyl-(3'----5')-thymidine containing either a cis-syn-cyclobutane pyrimidine dimer (d-T less than p greater than T) or a 6-4'-[pyrimidin-2'-one]pyrimidine photoproduct (d-T[p]-T), and similarly the UV-modified compounds of (dT)3 bearing either photoproduct at their 5'-end (d-T less than p greater than TpT and d-T[p]TpT). In contrast, the 3'-structural isomers of these trinucleotides (d-TpT less than p greater than T and d-TpT[p]T) were phosphorylated at the same rate as the parent compound. These phosphorylatable lesion-containing oligonucleotides are quantitatively released from UV-irradiated poly(dA):poly(dT) by enzymatic hydrolysis with snake venom phosphodiesterase and alkaline phosphatase (Liuzzi, M., Weinfeld, M., and Paterson, M. C. (1989) J. Biol. Chem. 264, 6355-6363). By combining this digestion regimen with phosphorylation by polynucleotide kinase and [gamma-32P]ATP, pyrimidine dimers were quantitated at the fmol level following exposure of poly(dA):poly(dT) and herring sperm DNA to biologically relevant UV fluences. The rate of dimer induction in the synthetic polymer, approximately 10 dimers/10(6) nucleotides/Jm-2, was in close agreement with that obtained by conventional methods. Dimers were induced at one-fourth of this rate in the natural DNA. Further treatment of the phosphorylated oligonucleotides derived from irradiated herring sperm DNA with nuclease P1 released the labeled 5'-nucleotide, thus permitting analysis of the nearest-neighbor bases 5' to the lesions. We observed a ratio for pyrimidine-to-purine bases of almost 6:1, implicating tripyrimidine stretches as hotspots for UV-induced DNA damage.     

Thermodynamics of unpaired terminal nucleotides on short RNA helixes correlates with stacking at helix termini in larger RNAs.

Free energies for stacking of unpaired nucleotides (dangling ends) at the termini of oligoribonucleotide Watson-Crick helixes (DeltaG(0)37,stack) depend on sequence for 3' ends but are always small for 5' ends. Here, these free energies are correlated with stacking at helix termini in a database of 34 RNA structures determined by X-ray crystallography and NMR spectroscopy. Stacking involving GA pairs is considered separately. A base is categorized as stacked by its distance from (相似文献   

Micrococcus luteus correndonucleases. II. Mechanism of action of two endonucleases specific for DNA containing pyrimidine dimers.

Py pyrimidine dimers Py correndonucleases I and II from Micrococcus luteus act exclusively on thymine-thymine, cytosine-cytosine, and thymine-cytosine cyclobutyl dimers in DNA, catalyzing incision 5' to the damage and generating 3'-hydroxyl and 5'-phosphoryl termini. Both enzymes initiate excision of pyrimidine dimers in vitro by correxonucleases and DNA polymerase I. The respective incised DNAs, however, differ in their ability to act as substrate for phage T4 polynucleotide ligase or bacterial alkaline phosphatase, suggesting that each endonuclease is specific for a conformationally unique site. The possibility that their respective action generates termini which represent different degrees of single strandedness is suggested by the unequal protection by Escherichia coli binding protein from the hydrolytic action of exonuclease VII.     

Different intrastrand and interstrand contributions to stacking account for roll variations at the alternating purine-pyrimidine sequences in A-DNA and A-RNA.

An explanation is suggested for the roll alternation between low and high values in A-type nucleic acid duplexes containing alternating sequences of purine and pyrimidine residues. The explanation combines two points. (1) Roll inevitably occurs in A-type duplexes due to geometrical reasons. (2) Intrastrand base stacking is much more impaired by roll than interstrand base stacking in A-type duplexes. Therefore purine-pyrimidine steps, whose bases mainly exhibit an intrastrand stacking, resist roll and decrease it. By contrast, bases at pyrimidine-purine steps exhibit a significant interstrand stacking that is tolerant to roll in A-type nucleic acid duplexes. In consequence, it is favourable if the purine-pyrimidine and pyrimidine-purine steps adopt low and high rolls, respectively in A-conformations of DNA and RNA molecules containing alternating purine-pyrimidine sequences. This is actually observed in the relevant molecular crystal structures.     

alpha-DNA-III. Characterization by high field 1H-NMR, anti-parallel self-recognition and conformation of the unnatural hexadeoxyribonucleotides alpha-[d(CpApTpGpCpG)] and alpha-[d(CpGpCpApTpG)]. Alpha-oligodeoxynucleotides as potential cellular probes for gene control.

High field 2-D-1H-NMR techniques permitted the assignment of all non-exchangeable protons of the unnatural deoxyribonucleotides alpha-[d(CpApTpGpCpG)] and alpha-[d(CpGpCpApTpG)]. 1-D and 2-D NOESY experiments show strong H6H8-H4' dipolar interactions for all nucleotides in both sequences. These data, together with COSY and J-resolved spectra, indicate that these two alpha-oligomers adopt 3'-exo conformations of the sugar moieties in solution with anti conformations of the glycosyl linkages. Both 1H-NMR data, and hypochromocity comparison of alpha-CATGCG and beta-CATGCG, demonstrate a higher degree of base stacking in the case of the alpha-sequence. The UV hyperchromicity at 260 nm, and symmetry considerations in the imino proton NMR experiments reveal antiparallel self-recognition and duplex annealing at positions 1-4 for alpha-[d(CATGCG)] and positions 3-6 for alpha-[d(CGCATG)]. The temperature variation of the imino proton NMR signals suggests that the hydrogen bonding in self-recognition is comparable in strength with that in a beta-DNA duplex, and NOE data are in accord with Watson-Crick rather than Hoogsteen base pairing.