The origins of high-affinity enzyme binding to an extrahelical DNA base

Base flipping is a highly conserved strategy used by enzymes to gain catalytic access to DNA bases that would otherwise be sequestered in the duplex structure. A classic example is the DNA repair enzyme uracil DNA glycosylase (UDG) which recognizes and excises unwanted uracil bases from DNA using a flipping mechanism. Previous work has suggested that enzymatic base flipping begins with dynamic breathing motions of the enzyme-bound DNA substrate, and then, only very late during the reaction trajectory do strong specific interactions with the extrahelical uracil occur. Here we report that UDG kinetically and thermodynamically prefers substrate sites where the uracil is paired with an unnatural adenine analogue that lacks any Watson-Crick hydrogen-bonding groups. The magnitude of the preference is a striking 43000-fold as compared to an adenine analogue that forms three H-bonds. Transient kinetic and fluorescence measurements suggest that preferential recognition of uracil in the context of a series of incrementally destabilized base pairs arises from two distinct effects: weak or absent hydrogen bonding, which thermodynamically assists extrusion, and, most importantly, increased flexibility of the site which facilitates DNA bending during base flipping. A coupled, stepwise reaction coordinate is implicated in which DNA bending precedes base pair rupture and flipping.     

Mutational analysis of the base-flipping mechanism of uracil DNA glycosylase

The DNA repair enzyme uracil DNA glycosylase (UDG) locates unwanted uracil bases in genomic DNA using a remarkable base-flipping mechanism in which the entire deoxyuridine nucleotide is rotated from the DNA base stack into the enzyme active site. Enzymatic base flipping has been described as a three-step process involving phosphodiester backbone pinching, base extrusion through active pushing and plugging by a leucine side chain that inserts in the DNA minor groove, and, finally, pulling by hydrogen-bonding groups that interact with the extrahelical base. Here we employ mutagenesis in combination with transient kinetic approaches to assess the functional roles of six conserved enzymatic groups of UDG that have been implicated in the "pinch, push, plug, and pull" base-flipping mechanism. Our results show that these mutant enzymes are capable of flipping the uracil base from the duplex, but that many of these mutations prevent a subsequent induced fit conformational step in which catalytic groups of UDG dock with the flipped-out base. These studies support our previous model for base flipping in which a conformational gating step closely follows base extrusion from the DNA duplex [Stivers, J. T., et al. (1999) Biochemistry 38, 952-963]. A model that accounts for the temporal and functional roles of these side chain interactions along the reaction pathway for base flipping is presented.     

The P1 phage replication protein RepA contacts an otherwise inaccessible thymine N3 proton by DNA distortion or base flipping

The RepA protein from bacteriophage P1 binds DNA to initiate replication. RepA covers one face of the DNA and the binding site has a completely conserved T that directly faces RepA from the minor groove at position +7. Although all four bases can be distinguished through contacts in the major groove of B-form DNA, contacts in the minor groove cannot easily distinguish between A and T bases. Therefore the 100% conservation at this position cannot be accounted for by direct contacts approaching into the minor groove of B-form DNA. RepA binding sites with modified base pairs at position +7 were used to investigate contacts with RepA. The data show that RepA contacts the N3 proton of T at position +7 and that the T=A hydrogen bonds are already broken in the DNA before RepA binds. To accommodate the N3 proton contact the T₊₇ /A_+7^′ base pair must be distorted. One possibility is that T₊₇ is flipped out of the helix. The energetics of the contact allows RepA to distinguish between all four bases, accounting for the observed high sequence conservation. After protein binding, base pair distortion or base flipping could initiate DNA melting as the second step in DNA replication.     

Solution structure of a nitrous acid induced DNA interstrand cross-link

Nitrous acid is a mutagenic agent. It can induce interstrand cross-links in duplex DNA, preferentially at d(CpG) steps: two guanines on opposite strands are linked via a single shared exocyclic imino group. Recent synthetic advances have led to the production of large quantities of such structurally homogenous cross-linked duplex DNA. Here we present the high resolution solution structure of the cross-linked dodecamer [d(GCATCCGGATGC)]2 (the cross-linked guanines are underlined), determined by 2D NMR spectroscopy, distance geometry, restrained molecular dynamics and iterative NOE refinement. The cross-linked guanines form a nearly planar covalently linked 'G:G base pair' with only minor propeller twisting, while the cytidine bases of their normal base pairing partners have been flipped out of the helix and adopt well defined extrahelical positions in the minor groove. On the 5'-side of the cross-link, the minor groove is widened to accommodate these extrahelical bases, and the major groove becomes quite narrow at the cross-link. The cross-linked 'G:G base pair' is well stacked on the spatially adjacent C:G base pairs, particularly on the 3'-side guanines. In addition to providing the first structure of a nitrous acid cross-link in DNA, these studies could be of major importance to the understanding of the mechanisms of nitrous acid cross-linking and mutagenicity, as well as the mechanisms responsible for its repair in intracellular environments. It is also the shortest DNA cross-link structure to be described.     

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.     

Mechanical Properties of Base-Modified DNA Are Not Strictly Determined by Base Stacking or Electrostatic Interactions

This work probes the mystery of what balance of forces creates the extraordinary mechanical stiffness of DNA to bending and twisting. Here we explore the relationship between base stacking, functional group occupancy of the DNA minor and major grooves, and DNA mechanical properties. We study double-helical DNA molecules substituting either inosine for guanosine or 2,6-diaminopurine for adenine. These DNA variants, respectively, remove or add an amino group from the DNA minor groove, with corresponding changes in hydrogen-bonding and base stacking energy. Using the techniques of ligase-catalyzed cyclization kinetics, atomic force microscopy, and force spectroscopy with optical tweezers, we show that these DNA variants have bending persistence lengths within the range of values reported for sequence-dependent variation of the natural DNA bases. Comparison with seven additional DNA variants that modify the DNA major groove reveals that DNA bending stiffness is not correlated with base stacking energy or groove occupancy. Data from circular dichroism spectroscopy indicate that base analog substitution can alter DNA helical geometry, suggesting a complex relationship among base stacking, groove occupancy, helical structure, and DNA bend stiffness.     

Mechanical Properties of Base-Modified DNA Are Not Strictly Determined by Base Stacking or Electrostatic Interactions

This work probes the mystery of what balance of forces creates the extraordinary mechanical stiffness of DNA to bending and twisting. Here we explore the relationship between base stacking, functional group occupancy of the DNA minor and major grooves, and DNA mechanical properties. We study double-helical DNA molecules substituting either inosine for guanosine or 2,6-diaminopurine for adenine. These DNA variants, respectively, remove or add an amino group from the DNA minor groove, with corresponding changes in hydrogen-bonding and base stacking energy. Using the techniques of ligase-catalyzed cyclization kinetics, atomic force microscopy, and force spectroscopy with optical tweezers, we show that these DNA variants have bending persistence lengths within the range of values reported for sequence-dependent variation of the natural DNA bases. Comparison with seven additional DNA variants that modify the DNA major groove reveals that DNA bending stiffness is not correlated with base stacking energy or groove occupancy. Data from circular dichroism spectroscopy indicate that base analog substitution can alter DNA helical geometry, suggesting a complex relationship among base stacking, groove occupancy, helical structure, and DNA bend stiffness.     

Solution structure of an O6-[4-oxo-4-(3-pyridyl)butyl]guanine adduct in an 11 mer DNA duplex: evidence for formation of a base triplex

The pyridyloxobutylating agents derived from metabolically activated tobacco-specific nitrosamines can covalently modify guanine bases in DNA at the O(6) position. The adduct formed, O(6)-[4-oxo-4-(3-pyridyl)butyl]guanine ([POB]dG), results in mutations that can lead to tumor formation, posing a significant cancer risk to humans exposed to tobacco smoke. A combined NMR-molecular mechanics computational approach was used to determine the solution structure of the [POB]dG adduct within an 11mer duplex sequence d(CCATAT-[POB]G-GCCC).d(GGGCCATATGG). In agreement with the NMR results, the POB ligand is located in the major groove, centered between the flanking 5'-side dT.dA and the 3'-side dG.dC base pairs and thus in the plane of the modified [POB]dG.dC base pair, which is displaced slightly into the minor groove. The modified base pair in the structure adopts wobble base pairing (hydrogen bonds between [POB]dG(N1) and dC(NH4) amino proton and between [POB]dG(NH2) amino proton and dC(N3)). A hydrogen bond appears to occur between the POB carbonyl oxygen and the partner dC's second amino proton. The modified guanine purine base, partner cytosine pyrimidine base, and POB pyridyl ring form a triplex via this unusual hydrogen-bonding pattern. The phosphodiester backbone twists at the lesion site, accounting for the unusual phosphorus chemical shift differences relative to those for the control DNA duplex. The helical distortions and wobble base pairing induced by the covalent binding of POB to the O(6)-position of dG help explain the significant decrease of 17.6 degrees C in melting temperature of the modified duplex relative to the unmodified control.     

Conformational features of benzo‐homologated yDNA duplexes by molecular dynamics simulation

yDNA is a base‐modified nucleic acid duplex containing size‐expanded nucleobases. Base‐modified nucleic acids could expand the genetic alphabet and thereby enhance the functional potential of DNA. Unrestrained 100 ns MD simulations were performed in explicit solvent on the yDNA NMR sequence [5′(yA T yA yA T yA T T yA T)₂] and two modeled yDNA duplexes, [5′(yC yC G yC yC G G yC G G)₂] and [(yT5′ G yT A yC yG C yA yG T3′)?(yA5′ C T C yG C G yT A yC A3′)]. The force field parameters for the yDNA bases were derived in consistent with the well‐established AMBER force field. Our results show that DNA backbone can withstand the stretched size of the bases retaining the Watson‐Crick base pairing in the duplexes. The duplexes retained their double helical structure throughout the simulations accommodating the strain due to expanded bases in the backbone torsion angles, sugar pucker and helical parameters. The effect of the benzo‐expansion is clearly reflected in the extended C1′‐C1′ distances and enlarged groove widths. The size expanded base modification leads to reduction in base pair twist resulting in larger overlapping area between the stacked bases, enhancing inter and intra strand stacking interactions in yDNA in comparison with BDNA. This geometry could favour enhanced interactions with the groove binders and DNA binding proteins., 2016. © 2015 Wiley Periodicals, Inc. Biopolymers 105: 55–64, 2016     

Base pair opening within B-DNA: free energy pathways for GC and AT pairs from umbrella sampling simulations

The conformational pathways and the free energy variations for base opening into the major and minor grooves of a B-DNA duplex are studied using umbrella sampling molecular dynamics simulations. We compare both GC and AT base pair opening within a double-stranded d(GAGAGAGAGAGAG)· d(CTCTCTCTCTCTC) oligomer, and we are also able to study the impact of opening on the conformational and dynamic properties of DNA and on the surrounding solvent. The results indicate a two-stage opening process with an initial coupling of the movements of the bases within the perturbed base pair. Major and minor groove pathways are energetically comparable in the case of the pyrimidine bases, but the major groove pathway is favored for the larger purine bases. Base opening is coupled to changes in specific backbone dihedrals and certain helical distortions, including untwisting and bending, although all these effects are dependent on the particular base involved. Partial opening also leads to well defined water bridging sites, which may play a role in stabilizing the perturbed base pairs.     

Flexible 5-guanidino-4-nitroimidazole DNA lesions: structures and thermodynamics

5-Guanidino-4-nitroimidazole (NI), derived from guanine oxidation by reactive oxygen and nitrogen species, contains an unusual flexible ring-opened structure, with nitro and guanidino groups which possess multiple hydrogen bonding capabilities. In vitro primer extension experiments with bacterial and mammalian polymerases show that NI incorporates C as well as A and G opposite the lesion, depending on the polymerase. To elucidate structural and thermodynamic properties of the mutagenic NI lesion, we have investigated the structure of the modified base itself and the NI-containing nucleoside with high-level quantum mechanical calculations and have employed molecular modeling and molecular dynamics simulations in solution for the lesion in B-DNA duplexes, with four partner bases opposite the NI. Our results show that NI adopts a planar structure at the damaged base level. However, in the nucleoside and in DNA duplexes, steric hindrance between the guanidino group and its linked sugar causes NI to be nonplanar. The NI lesion can adopt both syn and anti conformations on the DNA duplex level, with the guanidino group positioned in the DNA major and minor grooves, respectively; the specific preference depends on the partner base. On the basis of hydrogen bonding and stacking interactions, groove dimensions, and bending, we find that the least distorted NI-modified duplex contains partner C, consistent with observed incorporation of C opposite NI. However, hydrogen bonding interactions between NI and partner G or A are also found, which would be compatible with the observed mismatches.     

Dynamics in the U6 RNA intramolecular stem-loop: a base flipping conformational change

The U6 RNA intramolecular stem-loop (ISL) structure is an essential component of the spliceosome and binds a metal ion required for pre-messenger RNA splicing. The metal binding internal loop region of the stem contains a partially protonated C67-(+)A79 base pair (pK(a) = 6.5) and an unpaired U80 nucleotide that is stacked within the helix at pH 7.0. Here, we determine that protonation occurs with an exchange lifetime of approximately 20 micros and report the solution structures of the U6 ISL at pH 5.7. The differences between pH 5.7 and 7.0 structures reveal that the pH change significantly alters the RNA conformation. At lower pH, U80 is flipped out into the major groove. Base flipping involves a purine stacking interaction of flanking nucleotides, inversion of the sugar pucker 5' to the flipped base, and phosphodiester backbone rearrangement. Analysis of residual dipolar couplings as a function of pH indicates that base flipping is not restricted to a local conformational change. Rather, base flipping alters the alignment of the upper and lower helices. The alternative conformations of the U6 ISL reveal striking structural similarities with both the NMR and crystal structures of domain 5 of self-splicing group II introns. These structures suggest that base flipping at an essential metal binding site is a conserved feature of the splicing machinery for both the spliceosome and group II self-splicing introns.     

Turning On uracil-DNA glycosylase using a pyrene nucleotide switch

Base flipping is a highly conserved process by which enzymes swivel an entire nucleotide from the DNA base stack into their active site pockets. Uracil DNA glycosylase (UDG) is a paradigm enzyme that uses a base flipping mechanism to catalyze the hydrolysis of the N-glycosidic bond of 2'-deoxyuridine (2'-dUrd) in DNA as the first step in uracil base excision repair. Flipping of 2'-dUrd by UDG has been proposed to follow a "pushing" mechanism in which a completely conserved leucine side chain (Leu-191) is inserted into the DNA minor groove to expel the uracil. Here we report a novel implementation of the "chemical rescue" approach to show that the weak binding affinity and low catalytic activity of L191A or L191G can be completely or partially restored by substitution of a pyrene (Y) nucleotide wedge on the DNA strand opposite to the uracil base (U/A to U/Y). These results indicate that pyrene acts both as a wedge to push the uracil from the base stack in the free DNA and as a "plug" to hinder its reinsertion after base flipping. Pyrene rescue should serve as a useful and novel tool to diagnose the functional roles of other amino acid side chains involved in base flipping.     

三链DNA稳定性的影响因素

介绍了碱基组成、碱基修饰或替代、DNA骨架的修饰、DNA配体的结合及反应体系中的盐离子和pH值等因素对三链DNA稳定性的影响。对三链DNA稳定性研究中应注意的几个问题也作了讨论。