Structural studies of DNA fragments: the G.T wobble base pair in A, B and Z DNA; the G.A base pair in B-DNA

The crystal structures of five double helical DNA fragments containing non-Watson-Crick complementary base pairs are reviewed. They comprise four fragments containing G.T base pairs: two deoxyoctamers d(GGGGCTCC) and d(GGGGTCCC) which crystallise as A type helices; a deoxydodecamer d(CGCGAATTTGCG) which crystallises in the B-DNA conformation; and the deoxyhexamer d(TGCGCG), which crystallises as a Z-DNA helix. In all four duplexes the G and T bases form wobble base pairs, with bases in the major tautomer forms and hydrogen bonds linking N1 of G with O2 of T and O6 of G with N3 of T. The X-ray analyses establish that the G.T wobble base pair can be accommodated in the A, B or Z double helix with minimal distortion of the global conformation. There are, however, changes in base stacking in the neighbourhood of the mismatched bases. The fifth structure, d(CGCGAATTAGCG), contains the purine purine mismatch G.A where G is in the anti and A in the syn conformation. The results represent the first direct structure determinations of base pair mismatches in DNA fragments and are discussed in relation to the fidelity of replication and mismatch recognition.     

Structural Studies of DNA Fragments: The G·T Wobble Base Pair in A,B and Z DNA; The G·A Base Pair in B-DNA

Abstract

The crystal structures of five double helical DNA fragments containing non-Watson-Crick complementary base pairs are reviewed. They comprise four fragments containing G·T base pairs: two deoxyoctamers d(GGGGCTCC) and d(GGGGTCCC) which crystallise as A type helices; a deoxydodecamer d(CGCGAATTTGCG) which crystallises in the B-DNA conformation; and the deoxyhexamer d(TGCGCG), which crystallises as a Z-DNA helix. In all four duplexes the G and T bases form wobble base pairs, with bases in the major tautomer forms and hydrogen bonds linking N1 of G with 02 of T and 06 of G with N3 of T. The X-ray analyses establish that the G·T wobble base pair can be accommodated in the A, B or Z double helix with minimal distortion of the global conformation. There are, however, changes in base stacking in the neighbourhood of the mismatched bases. The fifth structure, d(CGCGAATTAGCG), contains the purine purine mismatch G·A where G is in the anti and A in the syn conformation. The results represent the first direct structure determinations of base pair mismatches in DNA fragments and are discussed in relation to the fidelity of replication and mismatch recognition.     

Noninteger pitch and nuclease sensitivity of chromatin DNA.

Assuming that variation of nuclease sensitivity along nucleosomal DNA can basically be attributed to orientations of sugar--phosphate bonds relative to histone core, the pitch of chromatin DNA is estimated to be 10.33--10.40 base pairs. This is in accordance both with the known measured average distance between cleavage sites (10.3--10.4 base pairs) and with published data on variation of relative sensitivities of these sites to nuclease attack. The variation can be explained solely as a result of the systematic change of orientation of sugar--phosphate bonds of sensitive sites without additional suggestions about local steric hindrances by histone molecules. According to the analysis locations of sites least sensitive to nuclease attack should not depend on kind of endonuclease though the stagger could differ. We conclude that the nucleosome core particle is axially symmetrical. The results strongly support the suggestion that DNA is wrapped around the histone octamer smoothly, without interruption of base-stacking interactions.     

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

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

Mesoscopic model parametrization of hydrogen bonds and stacking interactions of RNA from melting temperatures

Information about molecular interactions in DNA can be obtained from experimental melting temperature data by using mesoscopic statistical physics models. Here, we extend the technique to RNA and show that the new parameters correctly reproduce known properties such as the stronger hydrogen bonds of AU base pairs. We also were able to calculate a complete set of elastic constants for all 10 irreducible combinations of nearest neighbours (NNs). We believe that this is particularly useful as experimentally derived information about RNA elasticity is relatively scarce. The melting temperature prediction using the present model improves over those from traditional NN model, providing thus an alternative way to calculate these temperatures for RNA. Additionally, we calculated the site-dependent base pair oscillation to explain why RNA shows larger oscillation amplitudes despite having stronger AU hydrogen bonds.     

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.     

Evaluation of Elastic Rod Models with Long Range Interactions for Predicting Nucleosome Stability

Abstract

The ability of a dinucleotide-step based elastic-rod model of DNA to predict nucleosome binding free energies is investigated using four available sets of elastic parameters. We compare the predicted free energies to experimental values derived from nucleosome reconstitution experiments for 84 DNA sequences. Elastic parameters (conformation and stiffnessess) obtained from MD simulations are shown to be the most reliable predictors, as compared to those obtained from analysis of base-pair step melting temperatures, or from analysis of x-ray structures. We have also studied the effect of varying the folded conformation of nucleosomal DNA by means of our Fourier filtering knock-out and knock-in procedure. This study confirmed the above ranking of elastic parameters, and helped to reveal problems inherent in models using only a local elastic energy function. Long-range interactions were added to the elastic-rod model in an effort to improve its predictive ability. For this purpose a Debye-Huckel energy term with a single, homogenous point charge per base- pair was introduced. This term contains only three parameters,—its weight relative to the elastic energy, the Debye screening length, and a minimum sequence distance for including pairwise interactions between charges. After optimization of these parameters, our Debye-Huckel term is attractive, and yields the same level of correlation with experiment (R = 0.75) as was achieved merely by varying the nucleosomal shape in the elastic-rod model. We suggest this result indicates a linker DNA—histone attraction or, possibly, entropic effects, that lead to a stabilization of a nucleosome away from the ends of DNA segments longer than 147 bp. Such effects are not accounted for by a localized elastic energy model.     

Sequence-dependent Nucleosome Positioning

Eukaryotic DNA is organized into a macromolecular structure called chromatin. The basic repeating unit of chromatin is the nucleosome, which consists of two copies of each of the four core histones and DNA. The nucleosomal organization and the positions of nucleosomes have profound effects on all DNA-dependent processes. Understanding the factors that influence nucleosome positioning is therefore of general interest. Among the many determinants of nucleosome positioning, the DNA sequence has been proposed to have a major role. Here, we analyzed more than 860,000 nucleosomal DNA sequences to identify sequence features that guide the formation of nucleosomes in vivo. We found that both a periodic enrichment of AT base pairs and an out-of-phase oscillating enrichment of GC base pairs as well as the overall preference for GC base pairs are determinants of nucleosome positioning. The preference for GC pairs can be related to a lower energetic cost required for deformation of the DNA to wrap around the histones. In line with this idea, we found that only incorporation of both signal components into a sequence model for nucleosome formation results in maximal predictive performance on a genome-wide scale. In this manner, one achieves greater predictive power than published approaches. Our results confirm the hypothesis that the DNA sequence has a major role in nucleosome positioning in vivo.     

The effect of intrinsic curvature on conformational properties of circular DNA.

Both thermal fluctuations and the intrinsic curvature of DNA contribute to conformations of the DNA axis. We looked for a way to estimate the relative contributions of these two components of the double-helix curvature for DNA with a typical sequence. We developed a model and Monte Carlo procedure to simulate the Boltzmann distribution of DNA conformations with a specific intrinsic curvature. Two steps were used to construct the equilibrium conformation of the model chain. We first specified the equilibrium DNA conformation at the base pair level of resolution, using a set of the equilibrium dinucleotide angles and DNA sequence. This conformation was then approximated by the conformation of the model chain consisting of a reduced number of longer, straight cylindrical segments. Each segment of the chain corresponded to a certain number of DNA base pairs. We simulated conformational properties of nicked circular DNA for different sets of equilibrium dinucleotide angles, different random DNA sequences, and lengths. Only random sequences of DNA generated with equal probability of appearance for all types of bases at any site of the sequence were used. The results showed that for a broad range of intrinsic curvature parameters, the radius of gyration of DNA circles should be nearly independent of DNA sequence for all DNA lengths studied. We found, however, a DNA properly that should strongly depend on DNA sequence if the double helix has essential intrinsic curvature. This property is the equilibrium distribution of the linking number for DNA circles that are 300-1000 bp in length. We found that a large fraction of the distributions corresponding to random DNA sequences should have two separate maxima. The physical nature of this unexpected effect is discussed. This finding opens new opportunities for joined experimental and theoretical studies of DNA intrinsic curvature.     

Low temperature solution structures and base pair stacking of double helical d(CGTACG)(2)

Solution structures and base pair stacking of a self- complementary DNA hexamer d(CGTACG)(2) have been studied at 5, 10 and 15 degrees C, respectively. The stacking interactions among the center base pair steps of the DNA duplex are found to improve when the terminal base pairs became less stable due to end fraying. A new structural quantity, the stacking sum (Sigma(s)), is introduced to indicate small changes in the stacking overlaps between base pairs. The improvements in the stacking overlaps to maintain the double helical conformation are probably the cause for the observed temperature dependent structural changes in double helical DNA molecule. A detailed analysis of the helical parameters, backbone torsion angles, base orientations and sugar conformations of these structures has been performed.     

Versatile method employing basic techniques of genetic engineering to study the ability of low-molecular-weight compounds to bind covalently with DNA in cell-free systems

Numerous antitumor and carcinogenic compounds and free radicals are able to modify DNA by forming covalent bonds, mainly with nucleophilic centers in nucleobases. Such a binding is usually of utmost importance for the biological outcome. The level of DNA adducts formed by a given agent is in most cases extremely low; hence their detection is very difficult. Here we propose a simple approach, exploiting techniques widely used in genetic engineering, to demonstrate and characterize the covalent modification of a DNA fragment by any low-molecular-weight compound of interest in a cell-free system. The specifically designed, several-hundred-base-pairs-long double-stranded deoxyoligonucleotide (PCR amplified)--subject to modification--includes two restriction sites: one containing only GC base pairs recognized by restriction endonuclease MspI and the other including only AT base pairs recognized by restriction endonuclease Tru1I. The covalent modification of the restriction sites abolishes their recognition and thus cleavage by the endonucleases applied. The formation of DNA adducts is induced by incubating the oligonucleotide with increasing concentrations of a studied compound, in the appropriate activating system if required. Then, the modified oligonucleotide is submitted to digestion by the above-mentioned restriction endonucleases and the DNA fragments are separated by polyacrylamide gel electrophoresis. The inhibition of cleavage indicates the occurrence of covalent modification of the restriction site(s) while simultaneously pointing at the kind of base pairs involved in DNA adduct formation. The validation of the method was performed for two DNA binding antitumor compounds, cisplatin and CC-1065, which form adducts preferentially with guanine and adenine, respectively.     

Database of non-canonical base pairs found in known RNA structures

Atomic resolution RNA structures are being published at an increasing rate. It is common to find a modest number of non-canonical base pairs in these structures in addition to the usual Watson-Crick pairs. This database summarizes the occurrence of these rare base pairs in accordance with standard nomenclature. The database, http://prion.bchs.uh.edu/, contains information such as sequence context, sugar pucker conformation, anti / syn base conformations, chemical shift, p K (a)values, melting temperature and free energy. Of the 29 anticipated pairs with two or more hydrogen bonds, 20 have been encountered to date. In addition, four unexpected pairs with two hydrogen bonds have been reported bringing the total to 24. Single hydrogen bond versions of five of the expected geometries have been encountered among the single hydrogen bond interactions. In addition, 18 different types of base triplets have been encountered, each of which involves three to six hydrogen bonds. The vast majority of the rare base pairs are antiparallel with the bases in the anti configuration relative to the ribose. The most common are the GU wobble, the Sheared GA pair, the Reverse Hoogsteen pair and the GA imino pair.     

1H NMR investigation of the conformational states of DNA in nucleosome core particles.

In this study 1H NMR has been used to investigate the conformational state of DNA in nucleosome core particles. The nucleosome core particles exhibit partially resolved low field (10-15 ppm) spectra due to imino protons in Watson-Crick base pairs (one resonance per GC or AT base pair). To a first approximation, the spectrum is virtually identical with that of protein-free 140 base pair DNA, and from this observation we draw two important conclusions: (i) Since the low field spectra of DNA are known to be sensitive to conformation, the conformation of DNA in the core particles is essentially the same as that of free DNA (presumably B-form), (ii) since kinks occurring at a frequency at 1 in 10 or 1 in 20 base pairs would result in a core particle spectrum different from that of free DNA we find no NMR evidence supporting either the Crick-Klug or the Sobell models for kinking DNA around the core histones. Linewidth considerations indicate that the rotational correlation time for the core particles is approximately 1.5 X 10(-7) sec, whereas the end-over-end tumbling time of the free 140 base pair DNA is 3 X 10(-7) sec.     

Mismatch-induced conformational distortions in polymerase beta support an induced-fit mechanism for fidelity

Molecular dynamics simulations of DNA polymerase (pol) beta complexed with different incorrect incoming nucleotides (G x G, G x T, and T x T template base x incoming nucleotide combinations) at the template-primer terminus are analyzed to delineate structure-function relationships for aberrant base pairs in a polymerase active site. Comparisons, made to pol beta structure and motions in the presence of a correct base pair, are designed to gain atomically detailed insights into the process of nucleotide selection and discrimination. In the presence of an incorrect incoming nucleotide, alpha-helix N of the thumb subdomain believed to be required for pol beta's catalytic cycling moves toward the open conformation rather than the closed conformation as observed for the correct base pair (G x C) before the chemical reaction. Correspondingly, active-site residues in the microenvironment of the incoming base are in intermediate conformations for non-Watson-Crick pairs. The incorrect incoming nucleotide and the corresponding template residue assume distorted conformations and do not form Watson-Crick bonds. Furthermore, the coordination number and the arrangement of ligands observed around the catalytic and nucleotide binding magnesium ions are mismatch specific. Significantly, the crucial nucleotidyl transferase reaction distance (P(alpha)-O3') for the mismatches between the incoming nucleotide and the primer terminus is not ideally compatible with the chemical reaction of primer extension that follows these conformational changes. Moreover, the extent of active-site distortion can be related to experimentally determined rates of nucleotide misincorporation and to the overall energy barrier associated with polymerase activity. Together, our studies provide structure-function insights into the DNA polymerase-induced constraints (i.e., alpha-helix N conformation, DNA base pair bonding, conformation of protein residues in the vicinity of dNTP, and magnesium ions coordination) during nucleotide discrimination and pol beta-nucleotide interactions specific to each mispair and how they may regulate fidelity. They also lend further support to our recent hypothesis that additional conformational energy barriers are involved following nucleotide binding but prior to the chemical reaction.     

Vibrational fluctuations of hydrogen bonds in a DNA double helix with nonuniform base pairs.

The Green's function technique is applied to a study of breathing modes in a DNA double helix which contains a region of different base pairs from the rest of the double helix. The calculation is performed on a G-C helix in the B conformation with four consecutive base pairs replaced by A-T. The average stretch in hydrogen bonds is found amplified around the A-T base pair region compared with that of poly(dG)-poly(dC). This is likely related to the A-T regions lower stability against hydrogen bond melting. The A-T region may be considered to be the initiation site for melting in such a helix.     

Formation of purine-purine mispairs by Sulfolobus solfataricus DNA polymerase IV

DNA damage that stalls replicative polymerases can be bypassed with the Y-family polymerases. These polymerases have more open active sites that can accommodate modified nucleotides. The lack of protein-DNA interactions that select for Watson-Crick base pairs correlate with the lowered fidelity of replication. Interstrand hydrogen bonds appear to play a larger role in dNTP selectivity. The mechanism by which purine-purine mispairs are formed and extended was examined with Solfolobus solfataricus DNA polymerase IV, a member of the RAD30A subfamily of the Y-family polymerases, as is pol eta. The structures of the purine-purine mispairs were examined by comparing the kinetics of mispair formation with adenine versus 1-deaza- and 7-deazaadenine and guanine versus 7-deazaguanine at four positions in the DNA, the incoming dNTP, the template base, and both positions of the terminal base pair. The time course of insertion of a single dNTP was examined with a polymerase concentration of 50 nM and a DNA concentration of 25 nM with various concentrations of dNTP. The time courses were fitted to a first-order equation, and the first-order rate constants were plotted against the dNTP concentration to produce k pol and K d (dNTP) values. A decrease in k pol/ K d (dNTP) associated with the deazapurine substitution would indicate that the position is involved in a crucial hydrogen bond. During correct base pair formation, the adenine to 1-deazaadenine substitution in both the incoming dNTP and template base resulted in a >1000-fold decrease in k pol/ K d (dNTP), indicating that interstrand hydrogen bonds are important in correcting base pair formation. During formation of purine-purine mispairs, the k pol/ K d (dNTP) values for the insertion of dATP and dGTP opposite 7-deazaadenine and 7-deazaguanine were decreased >10-fold with respect to those of the unmodified nucleotides. In addition, the rate of incorporation of 1-deaza-dATP opposite guanine was decreased 5-fold. These results suggest that during mispair formation the newly forming base pair is in a Hoogsteen geometry with the incoming dNTP in the anti conformation and the template base in the syn conformation. These results indicate that Dpo4 holds the incoming dNTP in the normal anti conformation while allowing the template nucleotide to change conformations to allow reaction to occur. This result may be functionally relevant in the replication of damaged DNA in that the polymerase may allow the template to adopt multiple configurations.     

The role of the minor groove substituents in indirect readout of DNA sequence by 434 repressor

The sequence of non-contacted bases at the center of the 434 repressor binding site affects the strength of the repressor-DNA complex by influencing the structure and flexibility of DNA (Koudelka, G. B., and Carlson, P. (1992) Nature 355, 89-91). We synthesized 434 repressor binding sites that differ in their central sequence base composition to test the importance of minor groove substituents and/or the number of base pair hydrogen bonds between these base pairs on DNA structure and strength of the repressor-DNA complex. We show here that the number of base pair H-bonds between the central bases apparently has no role in determining the relative affinity of a DNA site for repressor. Instead we find that the affinity of DNA for repressor depends on the absence or presence the N2-NH(2) group on the purine bases at the binding site center. The N2-NH(2) group on bases at the center of the 434 binding site appears to destabilize 434 repressor-DNA complexes by decreasing the intimacy of the specific repressor-DNA contacts, while increasing the reliance on protein contacts to the DNA phosphate backbone. Thus, the presence of an N2-NH(2) group on the purines at the center of a binding site globally alters the precise conformation of the protein-DNA interface.