Molecular mechanics parameters for the FapydG DNA lesion

FapydG is a common oxidative DNA lesion involving opening of the imidazole ring. It shares the same precursor as 8-oxodG and can be excised by the same enzymes as 8-oxodG. However, the loss of the aromatic imidazole in FapydG results in a reduction of the double bond character between C5 and N7, with an accompanying increase in conformational flexibility. Experimental characterization of FapydG is hampered by high reactivity, and thus it is desirable to investigate structural details through computer simulation. We show that the existing Amber force field parameters for FapydG do not reproduce X-ray structural data. We employed quantum mechanics energy profile calculations to derive new molecular mechanics parameters for the rotation of the dihedral angles in the eximidazole moiety. Using these parameters, all-atom simulations in explicit water reproduce the nonplanar conformation of cFapydG in the crystal structure of the complex with L. lactis glycosylase Fpg. We note that the nonplanar structure is stabilized by an acidic residue that is not present in most Fpg sequences. Simulations of the E-->S mutant, as present in E. coli, resulted in a more planar conformation, suggesting that the highly nonplanar form observed in the crystal structure may not have direct biological relevance for FapydG.     

Bacterial base excision repair enzyme Fpg recognizes bulky N7-substituted-FapydG lesion via unproductive binding mode

Fpg is a bacterial base excision repair enzyme that removes oxidized purines from DNA. This work shows that Fpg and its eukaryote homolog Ogg1 recognize with high affinity FapydG and bulky N7-benzyl-FapydG (Bz-FapydG). The comparative crystal structure analysis of stable complexes between Fpg and carbocyclic cFapydG or Bz-cFapydG nucleoside-containing DNA provides the molecular basis of the ability of Fpg to bind both lesions with the same affinity and to differently process them. To accommodate the steric hindrance of the benzyl group, Fpg selects the adequate rotamer of the extrahelical Bz-cFapydG formamido group, forcing the bulky group to go outside the binding pocket. Contrary to the binding mode of cFapydG, the particular recognition of Bz-cFapydG leads the BER enzymes to unproductive complexes which would hide the lesion and slow down its repair by the NER machinery.     

Probing the requirements for recognition and catalysis in Fpg and MutY with nonpolar adenine isosteres

The Escherichia coli DNA repair enzymes Fpg and MutY are involved in the prevention of mutations resulting from 7,8-dihydro-8-oxo-2'-deoxyguanosine (OG) in DNA. The nonpolar isosteres of 2'-deoxyadenosine, 4-methylbenzimidazole beta-deoxynucleoside (B), and 9-methyl-1H-imidazo[4,5-b]pyridine beta-deoxynucleoside (Q), were used to examine the importance of hydrogen bonding within the context of DNA repair. Specifically, the rate of base removal under single-turnover conditions by the MutY and Fpg glycosylases from duplexes containing OG:B and OG:Q mismatches, relative to OG:A mismatches, was evalulated. The reaction of Fpg revealed a 5- and 10-fold increase in rate of removal of OG from duplexes containing OG:B and OG:Q base pairs, respectively, relative to an OG:A mispair. These results suggest that the lack of the ability to hydrogen bond to the opposite base facilitates removal of OG. In contrast, adenine removal catalyzed by MutY was much more efficient from an OG:A mispair-containing duplex (k2 = 12 +/- 2 min(-1)) compared to the removal of B from an OG:B duplex (k(obs) < 0.002 min(-1)). Surprisingly, MutY was able to catalyze base removal from the OG:Q-containing substrate (k2 = 1.2 +/- 0.2 min(-1)). Importantly, the B and Q analogues are not deleterious to high-affinity DNA binding by MutY. In addition, the B and Q analogues are more susceptible to acid-catalyzed depurination illustrating that the enzyme-catalyzed mechanism is distinct from the nonenzymatic mechanism. Taken together, these results point to the importance of both N7 and N3 in the mechanism of adenine excision catalyzed by MutY.     

Dynamic behavior of DNA base pairs containing 8-oxoguanine

The process by which DNA repair enzymes recognize and selectively excise damaged bases in duplex DNA is fundamental to our mechanistic understanding of these critical biological reactions. 8-Oxoguanine (8-oxoG) is the most common form of oxidative DNA damage; unrepaired, this lesion generates a G:C-->T:A mutation. Central to the recognition and repair of DNA damage is base extrusion, a process in which the damaged base lesion or, in some cases, its partner disengages from the helix and is bound to the enzyme's active site where base excision takes place. The conformation adopted by 8-oxoG in duplex DNA is affected by the base positioned opposite this lesion; conformational changes may also take place when the damaged base binds to its cognate repair enzyme. We performed unrestrained molecular dynamics simulations for several 13-mer DNA duplexes. Oligomers containing G:C and 8oxoG:C pairs adopted Watson-Crick geometries in stable B-form duplexes; 8oxoG showed increased local and global flexibility and a reduced barrier to base extrusion. Duplexes containing the G:A mismatch showed much larger structural fluctuations and failed to adopt a well-defined structure. For the 8oxoG:A mismatch that is recognized by the DNA glycosylase MutY, the damaged nucleoside underwent spontaneous and reproducible anti-->syn transitions. The syn conformation is thermodynamically preferred. Steric hindrance and unfavorable electrostatics associated with the 8oxoG O8 atom in the anti conformation were the major driving forces for this transition. Transition events follow two qualitatively different pathways. The overall anti-->syn transition rate and relative probability of the two transition paths were dependent on local sequence context. These simulations indicate that both the dynamic and equilibrium behavior of the duplex change as a result of oxidation; these differences may provide valuable new insight into the selective action of enzymes on damaged DNA.     

Surprising repair activities of nonpolar analogs of 8-oxoG expose features of recognition and catalysis by base excision repair glycosylases

Repair glycosylases locate and excise damaged bases from DNA, playing central roles in preservation of the genome and prevention of disease. Two key glycosylases, Fpg and hOGG1, function to remove the mutagenic oxidized base 8-oxoG (OG) from DNA. To investigate the relative contributions of conformational preferences, leaving group ability, enzyme-base hydrogen bonding, and nucleobase shape on damage recognition by these glycosylases, a series of four substituted indole nucleosides, based on the parent OG nonpolar isostere 2Cl-4F-indole, were tested as possible direct substrates of these enzymes in the context of 30 base pair duplexes paired with C. Surprisingly, single-turnover experiments revealed that Fpg-catalyzed base removal activity of two of the nonpolar analogs was superior to the native OG substrate. The hOGG1 glycosylase was also found to catalyze removal of three of the nonpolar analogs, albeit considerably less efficiently than removal of OG. Of note, the analog that was completely resistant to hOGG1-catalyzed excision has a chloro-substituent at the position of NH7 of OG, implicating the importance of recognition of this position in catalysis. Both hOGG1 and Fpg retained high affinity for the duplexes containing the nonpolar isosteres. These studies show that hydrogen bonds between base and enzyme are not needed for efficient damage recognition and repair by Fpg and underscore the importance of facile extrusion from the helix in its damaged base selection. In contrast, damage removal by hOGG1 is sensitive to both hydrogen bonding groups and nucleobase shape. The relative rates of excision of the analogs with the two glycosylases highlight key differences in their mechanisms of damaged base recognition and removal.     

Molecular dynamics simulations and thermodynamics analysis of DNA-drug complexes. Minor groove binding between 4',6-diamidino-2-phenylindole and DNA duplexes in solution

An extended set of nanosecond-scale molecular dynamics simulations of DNA duplex sequences in explicit solvent interacting with the minor groove binding drug 4',6-diamidino-2-phenylindole (DAPI) are investigated for four different and sequence specific binding modes. Force fields for DAPI have been parametrized to properly reflect its internal nonplanarity. Sequences investigated include the binding modes observed experimentally, that is, AATT in d(CGCGAATTCGCG)(2) and ATTG in d(GGCCAATTGG)(2) and alternative shifted binding modes ATTC and AATT, respectively. In each case, stable MD simulations are obtained, well reproducing specific hydration patterns seen in the experiments. In contrast to the 2.4 A d(CGCGAATTCGCG)(2) crystal structure, the DAPI is nonplanar, consistent with its gas-phase geometry and the higher resolution crystal structure. The simulations also suggest that the DAPI molecule is able to adopt different conformational substates accompanied by specific hydration patterns that include long-residing waters. The MM_PBSA technology for estimating relative free energies was utilized. The most consistent free energy results were obtained with an approach that uses a single trajectory of the DNA-DAPI complex to estimate all free energy terms. It is demonstrated that explicit inclusion of a subset of bound water molecules shifts the calculated relative binding free energies in favor of both crystallographically observed binding modes, underlining the importance of structured hydration.     

Oxidation of guanine bases mediated by a manganese-porphyrin-oligonucleotide conjugate onto a G-rich DNA target: 8-oxo-7,8-dihydroguanine is not the major lesion

The oxidative alkali-labile G lesions mediated by manganese porphyrin oligonucleotide conjugate on a DNA target could not be attibuted to the formation of 8-oxo-7,8-dihydroguanine since they were not substrate of the Fpg protein. In order to identify the nature of these lesions, an analysis of the oxidized derivatives of 2′-deoxyguanosine, imidazolone and oxazolone, generated by photooxidation, was efficiently performed by using the positive electrospray ionisation-mass spectrometry method.     

A prebiotic role for 8-oxoguanosine as a flavin mimic in pyrimidine dimer photorepair

Redox-active enzyme cofactors derived from ribonucleotides have been called "fossils of the RNA world," suggesting that early catalysts employed modified nucleobases to facilitate redox chemistry in primitive metabolism. Here, we show that the common oxidative damage product 8-oxo-7,8-dihydroguanine (OG), when incorporated into a DNA or RNA strand in proximity to a cyclobutane pyrimidine dimer, can mimic the function of a flavin in photorepair. The OG nucleotide acts catalytically in a mechanism consistent with that of photolyase in which the photoexcited state of the purine donates an electron to a pyrimidine dimer to initiate bond cleavage; subsequent back electron transfer regenerates OG. This unusual example of one form of DNA damage, oxidation, functioning to repair another, photodimerization, may provide insight into the origins of prebiotic redox processes.     

Identification of Escherichia coli 8-oxoguanine endonuclease

7,8-Dihydro-8-oxoguanine (oh8Gua) endonuclease is a DNA repair enzyme in Escherichia coli to remove oh8Gua, a promutagenic DNA adduct. Due to the unique mode of enzyme action and substrate specificity, this DNA repair enzyme has been suggested to be identical to 2,6-diamino-4-hydroxyformamidopyrimidine (Fapy)-DNA glycosylase (Fpg). However, oh8Gua endonuclease had not been definitely identified because it had not been homogeneously purified. In this study, we attempted to purify and identify the enzyme. Through several purification procedures, we obtained two proteins (32 kD and 29 kD). The larger protein co-migrated with Fpg in 12% SDS-PAGE gel. Sequences of N-terminal amino acids of these two proteins were identical to that of Fpg; the smaller one is a degraded product of oh8Gua endonuclease during purification steps. These results indicate that oh8Gua endonuclease is identical to Fpg, implying that oh8Gua in oxidatively damaged DNA rather than Fapy is more physiologically relevant substrate for Fpg.     

Formation of tricyclic [4.3.3.0] adducts between 8-oxoguanosine and tyrosine under conditions of oxidative DNA-protein cross-linking

8-Oxo-7,8-dihydro-2'-deoxyguanosine (OG), a prevalent product of oxidative stress on cellular DNA, is readily further oxidized forming adducts with nucleophiles. In the presence of tyrosine or p-cresol, an unusual tricyclo[4.3.3.0] adduct has been characterized in both nucleoside and oligodeoxynucleotide studies. The adduct is more stable in oligomers than nucleosides and undergoes slow reversion and hydration to spiroiminodihydantoin.     

Unzipping kinetics of duplex DNA containing oxidized lesions in an α-hemolysin nanopore

The unzipping kinetics for lesion-containing DNA duplexes was studied in an α-hemolysin (α-HL) nanopore. The lesion of focus was the guanine two-electron oxidation product, 8-oxo-7,8-dihydroguanine (OG), and its further oxidation products, the hydantoins guanidinohydantoin (Gh) and spiroiminodihydantoin (Sp). The voltage-driven unzipping of individual duplex DNA molecules with symmetrical overhangs was carried out by pulling one strand of the duplex through the α-HL channel using an electrical field. Entry from the 3' or 5' end produced distinct current blockages, allowing directional effects on unzipping kinetics to be investigated. We find that the strand dissociation of complementary duplexes or duplexes containing the slightly destabilizing lesion OG follows a first-order kinetic model, while opening of duplexes that contain the highly destabilizing lesions Gh or Sp is described by two sequential first-order reactions, in which the intermediate state is proposed to correspond to the duplex unzipped to the lesion site within the channel. The rate constants for strand separation of the duplexes containing single lesions were obtained from kinetic model fits to histograms of unzipping duration. For all duplexes, the rate constants for strand separation displayed a significant dependence on the direction of entry into the nanopore. For duplexes containing Gh, truncated duplexes were used to assign the measured rate constants for the first and second unzipping steps of symmetrically designed duplexes.     

Nanopore detection of 8-oxo-7,8-dihydro-2'-deoxyguanosine in immobilized single-stranded DNA via adduct formation to the DNA damage site

The ability to detect DNA damage within the context of the surrounding sequence is an important goal in medical diagnosis and therapies, but there are no satisfactory methods available to detect a damaged base while providing sequence information. One of the most common base lesions is 8-oxo-7,8-dihydroguanine, which occurs during oxidation of guanine. In the work presented here, we demonstrate the detection of a single oxidative damage site using ion channel nanopore methods employing α-hemolysin. Hydantoin lesions produced from further oxidation of 8-oxo-7,8-dihydroguanine, as well as spirocyclic adducts produced from covalently attaching a primary amine to the spiroiminodihydantoin lesion, were detected by tethering the damaged DNA to streptavidin via a biotin linkage and capturing the DNA inside an α-hemolysin ion channel. Spirocyclic adducts, in both homo- and heteropolymer background single-stranded DNA sequences, produced current blockage levels differing by almost 10% from those of native base current blockage levels. These preliminary studies show the applicability of ion channel recordings not only for DNA sequencing, which has recently received much attention, but also for detecting DNA damage, which will be an important component to any sequencing efforts.     

Substrate binding to mammalian 15-lipoxygenase

Lipoxygenases (LOs) are implicated in the regulation of metabolic processes and in several human diseases. Revealing their exact role is hindered by an incomplete understanding of their activity, including substrate specificity and substrate alignment in the active site. Recently, it has been proposed that the change in substrate specificity for arachidonic acid (AA) or linoleic acid (LA) could be part of an auto-regulatory mechanism related to cancer grow. Kinetic differences between reactions of 15-hLO with AA and LA have also led to the suggestion that the two substrates could present mechanistic differences. In the absence of a crystal structure for the substrate:15-LO complex, here we present an atomic-level study of catalytically competent binding modes for LA to rabbit 15-LO (15-rLO-1) and compare the results to our previous work on AA. Docking calculations, molecular dynamics simulations, re-docking and cross-docking calculations are all used to analyze the differences and similarities between the binding modes of the two substrates. Interestingly, LA seems to adapt more easily to the enzyme structure and differs from AA on some dynamical aspects that could introduce kinetic differences, as observed experimentally. Still, our study concludes that, despite the different chain lengths and number of insaturations between these two physiological substrates of 15-rLO-1, the enzyme seems to catalyze their hydroperoxidation by binding them with a common binding mode that leads to similar catalytically competent complexes.     

Synthesis, stability, and conformation of the formamidopyrimidine G DNA lesion

The formamidopyrimidine (FapydGua) lesion, derived from the nucleobase guanine, is a major DNA lesion involved in mutagenesis and carcinogenesis. To date, the chemical information available about this main lesion is very limited. Herein, we describe a synthesis and a detailed characterization of the acetyl-protected monomer of the FapydGua lesion. Stability studies in DMSO and in water/acetonitrile show that the N-glycosidic bond, previously thought to be highly labile, is much more stable than anticipated. Decomposition of the FapydGua lesion proceeds with half-life times of 37.8 h for the beta-anomer and 65.2 h for the alpha-anomer in water/acetonitrile. The relaxation time for the anomerization reaction was determined to tau = 6.5 h at room temperature. Most important, it was found that the formamido group, which is critical for the lesion recognition process by repair enzymes, is fixed in the cis-conformation in apolar solvents such as chloroform. This conformation enables the formation of a hydrogen bond between the carbonyl oxygen of the formamide and the NH of the N-glycosidic bond within the framework of a seven-membered ring system. This has consequences for the recognition of the lesion by repair enzymes (hOGG1 and Fpg protein). These enzymes were so far believed to recognize the carbonyl group of the FapydGua lesion. Our investigations show that this carbonyl group is not readily accessible because it is almost buried in the dominating cis-conformation. In agreement with the recent X-ray structure of hOGG1 in complex with 8-oxo-7,8-dihydroguanine-containing DNA, we can conclude that repair enzymes can contact both lesions only via the N(7)-H group, which is a hydrogen-bond acceptor in guanine.     

Insights into intrastrand cross-link lesions of DNA from QM/MM molecular dynamics simulations

DNA damages induced by oxidative intrastrand cross-links have been the subject of intense research during the past decade. Yet, the currently available experimental protocols used to isolate such lesions only allow to get structural information about linked dinucleotides. The detailed structure of the damaged DNA macromolecule has remained elusive. In this study we generated in silico the most frequent oxidative intrastrand cross-link adduct, G[8,5-Me]T, embedded in a solvated DNA dodecamer by means of quantum mechanics/molecular mechanics (QM/MM) Car-Parrinello simulations. The free energy of activation required to bring the reactant close together and to form the C-C covalent-bond is estimated to be ~10 kcal/mol. We observe that the G[8,5-Me]T tandem lesion is accommodated with almost no perturbation of the Watson-Crick hydrogen-bond network and induces bend and unwinding angles of ~20° and 8°, respectively. This rather small structural distortion of the DNA macromolecule compared to other well characterized intrastrand cross-links, such as cyclobutane pyrimidines dimers or cisplatin-DNA complex adduct, is a probable rationale for the known lack of efficient repair of oxidative damages.     

Oxidation of guanine in double-stranded DNA by [Ru(bpy)2dppz]Cl2 in cationic reverse micelles

DNA oxidation has been investigated in the medium of cationic reverse micelles (RMs). The oxidative chemistry is photochemically initiated using the DNA intercalator bis(bipyridine)dipyridophenazine ruthenium(II) chloride ([Ru(bpy)2dppz]Cl2) bound to duplex DNA in the RMs. High-resolution polyacrylamide gel electrophoresis (PAGE) is used to reveal and quantify guanine (G) oxidation products, including 8-oxo-7,8-dihydroguanine (8OG). In buffer solution, the addition of the oxidative quenchers potassium ferricyanide or pentaamminechlorocobalt(III) dichloride leads to an increase in the amount of piperidine-labile G oxidation products generated via one-electron oxidation. In RMs, however, the yield of oxidatively generated damage is attenuated. With or without ferricyanide quencher in the RMs, the yield of oxidatively generated products is approximately the same. Inclusion of the cationic quencher [CoCl(NH3)5]2+ in the RMs increases the amount of oxidation products generated but not to the extent that it does in buffer solution. Under anaerobic conditions, all of the samples in RMs, with or without added oxidative quenchers, show decreased levels of piperidine-labile oxidation products, suggesting that the primary oxidant in RMs is singlet oxygen. G oxidation is enhanced in D2O and deuterated heptane and is diminished in the presence of sodium azide in RMs, also supporting 1O2 as the main G oxidant in RMs. Isotopic labeling experiments show that the oxygen atom in 8OG produced in RMs is not from water. The observed change in the G oxidation mechanism from a one-electron process in buffer to mostly 1O2 in RMs illustrates the importance of both DNA structure and DNA environment on the chemistry of G oxidation.     

Base excision repair synthesis of DNA containing 8-oxoguanine in Escherichia coli

8-oxo-7,8-dihydroguanine (8-oxo-G) in DNA is a mutagenic adduct formed by reactive oxygen species. In Escherichia coli, 2,6-dihydroxy-5N-formamidopyrimidine (Fapy)-DNA glycosylase (Fpg) removes this mutagenic adduct from DNA. In this report, we demonstrate base excision repair (BER) synthesis of DNA containing 8-oxo-G with Fpg in vitro. Fpg cut the oligonucleotide at the site of 8-oxo-G, producing one nucleotide gap with 3' and 5' phosphate termini. Next, 3' phosphatase(s) in the supernatant obtained by precipitating a crude extract of E. coli with 40% ammonium sulfate, removed the 3' phosphate group at the gap, thus exposing the 3' hydroxyl group to prime DNA synthesis. DNA polymerase and DNA ligase then completed the repair. These results indicate the biological significance of the glycosylase and apurinic/ apyrimidinic (AP) lyase activities of Fpg, in concert with 3' phosphatase(s) to create an appropriately gapped substrate for efficient BER synthesis of DNA containing 8-oxo-G.     

Mass spectrometric studies of alkali metal ion binding on thrombin-binding aptamer DNA

The binding sites and consecutive binding constants of alkali metal ions, (M⁺ = Na⁺, K⁺, Rb⁺, and Cs⁺), to thrombin-binding aptamer (TBA) DNA were studied by Fourier-transform ion cyclotron resonance spectrometry. TBA-metal complexes were produced by electrospray ionization (ESI) and the ions of interest were mass-selected for further characterization. The structural motif of TBA in an ESI solution was checked by circular dichroism. The metal-binding constants and sites were determined by the titration method and infrared multiphoton dissociation (IRMPD), respectively. The binding constant of potassium is 5–8 times greater than those of other alkali metal ions, and the potassium binding site is different from other metal binding sites. In the 1:1 TBA-metal complex, potassium is coordinated between the bottom G-quartet and two adjacent TT loops of TBA. In the 1:2 TBA—metal complex, the second potassium ion binds at the TGT loop of TBA, which is in line with the antiparallel G-quadruplex structure of TBA. On the other hand, other alkali metal ions bind at the lateral TGT loop in both 1:1 and 1:2 complexes, presumably due to the formation of ion-pair adducts. IRMPD studies of the binding sites in combination with measurements of the consecutive binding constants help elucidate the binding modes of alkali metal ions on DNA aptamer at the molecular level.     

Forever Young: Structural Stability of Telomeric Guanine Quadruplexes in the Presence of Oxidative DNA Lesions**

Human telomeric DNA, in G-quadruplex (G4) conformation, is characterized by a remarkable structural stability that confers it the capacity to resist to oxidative stress producing one or even clustered 8-oxoguanine (8oxoG) lesions. We present a combined experimental/computational investigation, by using circular dichroism in aqueous solutions, cellular immunofluorescence assays and molecular dynamics simulations, that identifies the crucial role of the stability of G4s to oxidative lesions, related also to their biological role as inhibitors of telomerase, an enzyme overexpressed in most cancers associated to oxidative stress.