High base pair opening rates in tracts of GC base pairs

Sequence-dependent structural features of the DNA double helix have a strong influence on the base pair opening dynamics. Here we report a detailed study of the kinetics of base pair breathing in tracts of GC base pairs in DNA duplexes derived from 1H NMR measurements of the imino proton exchange rates upon titration with the exchange catalyst ammonia. In the limit of infinite exchange catalyst concentration, the exchange times of the guanine imino protons of the GC tracts extrapolate to much shorter base pair lifetimes than commonly observed for isolated GC base pairs. The base pair lifetimes in the GC tracts are below 5 ms for almost all of the base pairs. The unusually rapid base pair opening dynamics of GC tracts are in striking contrast to the behavior of AT tracts, where very long base pair lifetimes are observed. The implication of these findings for the structural principles governing spontaneous helix opening as well as the DNA-binding specificity of the cytosine-5-methyltransferases, where flipping of the cytosine base has been observed, are discussed.     

Base pair opening in three DNA-unwinding elements

DNA-unwinding elements are specific base sequences that are located in the origin of DNA replication where they provide the start point for strand separation and unwinding of the DNA double helix. In the present work we have obtained the first characterization of the opening of individual base pairs in DNA-unwinding elements. The three DNA molecules investigated reproduce the 13-mer DNA-unwinding elements present in the Escherichia coli chromosome. The base sequences of the three 13-mers are conserved in the origins of replication of enteric bacterial chromosomes. The exchange of imino protons with solvent protons was measured for each DNA as a function of the concentration of exchange catalyst using nuclear magnetic resonance spectroscopy. The exchange rates provided the rates and the equilibrium constants for opening of individual base pairs in each DNA at 20 degrees C. The results reveal that the kinetics and energetics of the opening reactions for AT/TA base pairs are different in the three DNA-unwinding elements due to long range effects of the base sequence. These differences encompass the AT/TA base pairs that are conserved in various bacterial genomes. Furthermore, a qualitative correlation is observed between the kinetics and energetics of opening of AT/TA base pairs and the location of the corresponding DNA-unwinding element in the origin of DNA replication.     

Enhanced base-pair opening in the adenine tract of a RNA double helix

Proton exchange and nuclear magnetic resonance spectroscopy are being used to characterize the kinetics and energetics of base-pair opening in two nucleic acid double helices. One is the RNA duplex 5'-r(GCGAUAAAAAGGCC)-3'/5'-r(GGCCUUUUUAUCGC)-3', which contains a central tract of five AU base pairs. The other is the homologous DNA duplex with a central tract of five AT base pairs. The rates and the equilibrium constants of the opening reaction of each base pair are measured from the dependence of the exchange rates of imino protons on ammonia concentration, at 10 °C. The results reveal that the tract of AU base pairs in the RNA duplex differs from the homologous tract of AT base pairs in DNA in several ways. The rates of opening of AU base pairs in RNA are high and increase progressively along the tract, reaching their largest values at the 3'-end of the tract. In contrast, the opening rates of AT base pairs in DNA are much lower than those of AU base pairs. Within the tract, the largest opening rate is observed for the AT base pair at the 5'-end of the tract. These differences in opening kinetics are paralleled by differences in the stabilities of individual base pairs. All AU base pairs in the RNA are less stable than the AT base pairs in the DNA. The presence of the tract enhances these differences by increasing the stability of AT base pairs in DNA while decreasing the stability of AU base pairs in RNA. Due to these divergent trends, along the tracts, the AU base pairs become progressively less stable than AT base pairs. These findings demonstrate that tracts of AU base pairs in RNA have specific dynamic and energetic signatures that distinguish them from similar tracts of AT base pairs in DNA.     

RNA structure and dynamics: A base pairing perspective

RNA is now known to possess various structural, regulatory and enzymatic functions for survival of cellular organisms. Functional RNA structures are generally created by three-dimensional organization of small structural motifs, formed by base pairing between self-complementary sequences from different parts of the RNA chain. In addition to the canonical Watson–Crick or wobble base pairs, several non-canonical base pairs are found to be crucial to the structural organization of RNA molecules. They appear within different structural motifs and are found to stabilize the molecule through long-range intra-molecular interactions between basic structural motifs like double helices and loops. These base pairs also impart functional variation to the minor groove of A-form RNA helices, thus forming anchoring site for metabolites and ligands. Non-canonical base pairs are formed by edge-to-edge hydrogen bonding interactions between the bases. A large number of theoretical studies have been done to detect and analyze these non-canonical base pairs within crystal or NMR derived structures of different functional RNA. Theoretical studies of these isolated base pairs using ab initio quantum chemical methods as well as molecular dynamics simulations of larger fragments have also established that many of these non-canonical base pairs are as stable as the canonical Watson–Crick base pairs. This review focuses on the various structural aspects of non-canonical base pairs in the organization of RNA molecules and the possible applications of these base pairs in predicting RNA structures with more accuracy.     

A review of recent experiments on step-to-step "hand-off" of the DNA intermediates in mammalian base excision repair pathways

The current "working model" for mammalian base excision repair involves two sub-pathways termed single-nucleotide base excision repair and long patch base excision repair that are distinguished by their repair patch sizes and the enzymes/co-factors involved. These base excision repair sub-pathways are designed to sequester the various DNA intermediates, passing them along from one step to the next without allowing these toxic molecules to trigger cell cycle arrest, necrotic cell death, or apoptosis. Although a variety of DNA-protein and protein-protein interactions are known for the base excision repair intermediates and enzymes/co-factors, the molecular mechanisms accounting for step-to-step coordination are not well understood. In this review, we explore the question of whether there is an actual step-to-step "hand-off" of the DNA intermediates during base excision repair in vitro. The results show that when base excision repair enzymes are pre-bound to the initial single-nucleotide base excision repair intermediate, the DNA is channeled from apurinic/apyrimidinic endonuclease 1 to DNA polymerase beta and then to DNA ligase. In the long patch base excision repair sub-pathway, where the 5'-end of the incised strand is blocked, the intermediate after polymerase beta gap filling is not channeled from polymerase beta to the subsequent enzyme, flap endonuclease 1. Instead, flap endonuclease 1 must recognize and bind to the intermediate in competition with other molecules.     

Energetics of base pair opening in a DNA dodecamer containing an A3T3 tract.

Nuclear magnetic resonance spectroscopy has been used to characterize the kinetics and energetics of opening of base pairs in the DNA dodecamer [d(CGCAAATTTGCG)]2. The dodecamer contains an A3T3 tract that induces intrinsic curvature of the helix axis. Previous studies from this and other laboratories have shown that the kinetics of base pair opening in AnTn tracts is unique: the opening rates of the A.T base pairs in the interior of the tract are much lower than that of the A.T base pair at the 5'-end of the tract. In the present work, we have investigated the energetics of the pathways for opening of the A.T base pairs in the A3T3 tract. The energetic parameters of the activated state(s) are obtained from the temperature dependence of the opening rate constants. The lower opening rates for the A.T base pairs situated in the interior of the tract are shown to originate from higher activation enthalpies which are compensated, in part, by increases in the activation entropies. We have also obtained an energetic characterization of the open state(s) of the A.T base pairs in the dodecamer by measuring the equilibrium constants for base pair opening and their temperature dependence. The results suggest that the transitions from closed to open state(s) in the A.T base pairs of the A3T3 tract are energetically similar.     

Mechanistic details of the DNA recognition by the Dnmt1 DNA methyltransferase

A recently solved Dnmt1-DNA crystal structure revealed several enzyme-DNA contacts and large structural rearrangements of the DNA at the target site, including the flipping of the non-target strand base of the base pair flanking the CpG site and formation of a non-canonical base pair between the non-target strand Gua and the flanking base pair. Here, we show that the contacts of the enzyme to the target base and the Gua:5mC base pair that are observed in the structure are very important for catalytic activity. The contacts to the non-target strand Gua are not important since its exchange by Ade stimulated activity. Except target base flipping, we could not find evidence that the DNA rearrangements have a functional role.     

New relations for steady-state enzyme kinetics and their application to rapid equilibrium assumption

With the use of graph theory, new relations for steady-state enzyme kinetics are derived and strictly proved for an arbitrary mechanism of enzyme-catalyzed reaction containing a reversible segment. Using these relations, a general principle for rapid equilibrium assumption is formulated and proved: the reversible bound segment can be considered as an equilibrium segment only when the values of the base trees that are not proper to this segment can be neglected (within prescribed accuracy) in relation to the values of the base trees that belong to this segment. In contrast with the foreign base trees, the base trees that are proper to the segment have the following properties: the tree that is directed to the base within this segment does not contain edges leaving this segment; and the tree that is directed to the base outside the segment contains only one edge leaving this segment. Equilibrium variations are assessed for steady-state concentrations of intermediates in the equilibrium segment, numerical expressions are obtained for the accuracy of determination of intermediate concentrations as well as for the accuracy of determination of the rate of enzyme-catalyzed reaction under rapid equilibrium assumption.     

Using an RNA secondary structure partition function to determine confidence in base pairs predicted by free energy minimization

A partition function calculation for RNA secondary structure is presented that uses a current set of nearest neighbor parameters for conformational free energy at 37 degrees C, including coaxial stacking. For a diverse database of RNA sequences, base pairs in the predicted minimum free energy structure that are predicted by the partition function to have high base pairing probability have a significantly higher positive predictive value for known base pairs. For example, the average positive predictive value, 65.8%, is increased to 91.0% when only base pairs with probability of 0.99 or above are considered. The quality of base pair predictions can also be increased by the addition of experimentally determined constraints, including enzymatic cleavage, flavin mono-nucleotide cleavage, and chemical modification. Predicted secondary structures can be color annotated to demonstrate pairs with high probability that are therefore well determined as compared to base pairs with lower probability of pairing.     

Mutagenesis and comparative sequence analysis of a base triple joining the two domains of group I ribozymes.

Tertiary interactions are important in the higher-order folding of catalytic RNAs. Recently, a base triple, joining the two major domains of the catalytic core, was determined in group I introns from the cyanobacterium Anabaena PCC7120 and the eukaryote Tetrahymena thermophila. This base triple involves the fifth base pair of P4 and the fifth base of the single-stranded region J8/7. We made base pair and single-nucleotide substitutions in the fifth base pair of P4, a G-C in the wild-type Anabaena intron, and tested them for self-splicing activity. The results suggest a hydrogen bonding model in which only the C of the base pair interacts directly with the fifth base of J8/7. Comparative sequence analysis was used to determine the different combinations of base triples that occur in approximately 450 natural group I introns identified to date. About 94% of the base triples analyzed are compatible with the proposed hydrogen bonding model. Disrupting this base triple in the Tetrahymena intron resulted in the disappearance of splicing intermediates (intron 3' exon and 5' exon), even though the first step of splicing was not affected. Restoration of the base triple by a compensatory mutation reverted the intermediates to wild-type levels. These results suggest that disruption of the base triple increases the rate of the second step of splicing or of a conformational change preceding the second step. Repositioning of the base triple to form a new set of interactions may be required for the second step of splicing.     

An investigation of genomic base distribution

While veritable oceans of ink have been spilled over the base distributions within genes, the literature is virtually silent on large scale intra genomic base distribution. To address this issue, we have examined approximately 3400 chromosomal sequences from approximately 2000 entire genomes-including DNA and RNA, single- and double-stranded, coding and non-coding genomes. For each sequence the mean, variance, skewness, and kurtosis for each base were computed along with the genome base composition. The main findings are: (1) there is no simple relationship between these statistics and the base composition of the genome, (2) in non-viral genomes, base distribution is non-uniform, (3) base distribution in non-eukaryotic genomes obeys a number of simple rules, (4) these rules are not dependent on the presence of coding sequences, (5) bacterial genomes in particular are unusually compliant with these rules, and (6) eukaryotes have a unique pattern of base distribution.     

Synthesis and characterization of oligodeoxynucleotides containing the mutagenic base analogue 4-O-ethylthymine.

A method for the preparation of oligonucleotides containing the mutagenic base 4-O-ethylthymine is described for the first time. Use of p-nitrophenylethyl type base protecting groups together with phosphitetriester solid-phase methodology makes possible the rapid and efficient preparation of oligonucleotides bearing 4-O-ethylthymine, while standard base protecting groups are not compatible with the presence of this base. Possible applications of this methodology are discussed.     

A review of recent experiments on step-to-step “hand-off” of the DNA intermediates in mammalian base excision repair pathways

The current “working model” for mammalian base excision repair involves two sub-pathways termed single-nucleotide base excision repair and long patch base excision repair that are distinguished by their repair patch sizes and the enzymes/co-factors involved. These base excision repair sub-pathways are designed to sequester the various DNA intermediates, passing them along from one step to the next without allowing these toxic molecules to trigger cell cycle arrest, necrotic cell death, or apoptosis. Although a variety of DNA-protein and protein-protein interactions are known for the base excision repair intermediates and enzymes/co-factors, the molecular mechanisms accounting for step-to-step coordination are not well understood. In this review, we explore the question of whether there is an actual step-to-step “hand-off” of the DNA intermediates during base excision repair in vitro. The results show that when base excision repair enzymes are pre-bound to the initial single-nucleotide base excision repair intermediate, the DNA is channeled from apurinic/apyrimidinic endonuclease 1 to DNA polymerase β and then to DNA ligase. In the long patch base excision repair sub-pathway, where the 5′-end of the incised strand is blocked, the intermediate after polymerase β gap filling is not channeled from polymerase β to the subsequent enzyme, flap endonuclease 1. Instead, flap endonuclease 1 must recognize and bind to the intermediate in competition with other molecules.     

Thermodynamics and kinetics for base-pair opening in the P1 duplex of the Tetrahymena group I ribozyme

The thermodynamics and kinetics for base-pair opening of the P1 duplex of the Tetrahymena group I ribozyme were studied by NMR hydrogen exchange experiments. The apparent equilibrium constants for base pair opening were measured for most of the imino protons in the P1 duplex using the base catalysts NH3, HPO4(2-) or TRIS. These equilibrium constants were also measured for several modified P1 duplexes, and the C-2.G23 base pair was the most stable base pair in all the duplexes. The conserved U-1*G22 base pair is required for activity of the ribozyme and the data here show that this wobble base pair destabilizes neighboring base pairs on only one side of the wobble. A 2'-OMe modification on the U-3 residue stabilized its own base pair but had little effect on the neighboring base pairs. Three base pairs, U-1*G22, C-2*G23 and A2*U21 showed unusual equilibrium constants for opening and possible implications of the opening thermodynamics of these base pairs on the undocking rates of the P1 helix with catalytic core are discussed.     

Structure of transfer RNA molecules containing the long variable loop

A structure is proposed for the type II tRNA molecules containing the long variable loop and the tertiary base interactions here are compared with type I tRNAs having the short variable loop. The type II tRNAs are similar to the type I tRNAs in their tertiary base pairing interactions but differ from them generally by not having the tertiary base triples. The long variable loop, which is comprised of a helical stem and a loop at the end of it, emerges from the deep groove side of the dihydrouridine helix, and is tilted roughly 30° to the plane formed by the amino acid-pseudouridine and anticodon-dihydrouridine helices found in yeast tRNA^Phe. The fact that many of the type I tRNAs also lack the full compliment of base triples suggests that the tertiary base pairs may alone suffice to sustain the tRNA fold required for its biological function. The base triples and the variable loop appear to have little functional significance. The base type at position 9 is correlated with the number of base triples and G-C base pairs in the dihydrouridine stem.     

Published data on mutagenesis by ionizing radiation of plasmids in solution probably reflect in part the specificity of adventitious transition metal ions complexed to the DNA

A number of recent papers show that single base changes induced by mutagenesis with ionizing radiation of genes on plasmids in solution, followed by transfection into mammalian or bacterial cells for assay, are mostly at G:C base pairs, with mutagenic hot spots. Genes irradiated in mammalian or bacterial cells, on the other hand, have comparable numbers of base changes at all sites, with no evidence for hot spots. The differences are ascribed to induction of many base change mutations in vitro by reactions catalyzed by adventitious transition metal ions complexed to the DNA. Reasons are given why this process should play a much smaller role in vivo.     

Influence of a fluorobenzene nucleobase analogue on the conformational flexibility of RNA studied by molecular dynamics simulations

Chemically modified bases are frequently used to stabilize nucleic acids, to study the driving forces for nucleic acid structure formation and to tune DNA and RNA hybridization conditions. In particular, fluorobenzene and fluorobenzimidazole base analogues can act as universal bases able to pair with any natural base and to stabilize RNA duplex formation. Although these base analogues are compatible with an A-form RNA geometry, little is known about the influence on the fine structure and conformational dynamics of RNA. In the present study, nano-second molecular dynamics (MD) simulations have been performed to characterize the dynamics of RNA duplexes containing a central 1'-deoxy-1'-(2,4-difluorophenyl)-beta-D-ribofuranose base pair or opposite to an adenine base. For comparison, RNA with a central uridine:adenine pair and a 1'-deoxy-1'-(phenyl)-beta-D-ribofuranose opposite to an adenine was also investigated. The MD simulations indicate a stable overall A-form geometry for the RNAs with base analogues. However, the presence of the base analogues caused a locally enhanced mobility of the central bases inducing mainly base pair shear and opening motions. No stable 'base-paired' geometry was found for the base analogue pair or the base analogue:adenine pairs, which explains in part the universal base character of these analogues. Instead, the conformational fluctuations of the base analogues lead to an enhanced accessibility of the bases in the major and minor grooves of the helix compared with a regular base pair.     

On the nature of stability of the nucleotide base associates in water solution

Expression for the long-range intermolecular interaction energy obtained by the perturbation theory method in atomic dipole approximation are used for the study of the nature of base interaction in stacked dimers formed of two neighbouring DNA base pairs. Base wave functions are computed by the CNDO-CI method. The inplane interactions are shown to give the dominant contribution into the DNA stabilization energy in vacuum. The estimations performed for the solvent effect on intermolecular interaction energy allowed us to draw conclusions about the decisive role of hydrophobic interactions in a base stacking.     

Published data on mutagenesis by ionizing radiation of plasmids in solution probably reflect in part the specificity of adventitious transition metal ions complexed to the DNA.

A number of recent papers show that single base changes induced by mutagenesis with ionizing radiation of genes on plasmids in solution, followed by transfection into mammalian or bacterial cells for assay, are mostly at G:C base pairs, with mutagenic hot spots. Genes irradiated in mammalian or bacterial cells, on the other hand, have comparable numbers of base changes at all sites, with no evidence for hot spots. The differences are ascribed to induction of many base change mutations in vitro by reactions catalyzed by adventitious transition metal ions complexed to the DNA. Reasons are given why this process should play a much smaller role in vivo.