A dynamic checkpoint in oxidative lesion discrimination by formamidopyrimidine-DNA glycosylase

Dynamics of 8-Oxoguanine in DNA: Decisive Effects of Base Pairing and Nucleotide Context

8-Oxo-7,8-dihydroguanine (oxoG), an abundant DNA lesion, can mispair with adenine and induce mutations. To prevent this, cells possess DNA repair glycosylases that excise either oxoG from oxoG:C pairs (bacterial Fpg, human OGG1) or A from oxoG:A mispairs (bacterial MutY, human MUTYH). Early lesion recognition steps remain murky and may include enforced base pair opening or capture of a spontaneously opened pair. We adapted the CLEANEX-PM NMR protocol to detect DNA imino proton exchange and analyzed the dynamics of oxoG:C, oxoG:A, and their undamaged counterparts in nucleotide contexts with different stacking energy. Even in a poorly stacking context, the oxoG:C pair did not open easier than G:C, arguing against extrahelical base capture by Fpg/OGG1. On the contrary, oxoG opposite A significantly populated the extrahelical state, which may assist recognition by MutY/MUTYH.

Distinct Mechanisms of Target Search by Endonuclease VIII-like DNA Glycosylases

Proteins that recognize specific DNA sequences or structural elements often find their cognate DNA lesions in a processive mode, in which an enzyme binds DNA non-specifically and then slides along the DNA contour by one-dimensional diffusion. Opposite to the processive mechanism is distributive search, when an enzyme binds, samples and releases DNA without significant lateral movement. Many DNA glycosylases, the repair enzymes that excise damaged bases from DNA, use processive search to find their cognate lesions. Here, using a method based on correlated cleavage of multiply damaged oligonucleotide substrates we investigate the mechanism of lesion search by three structurally related DNA glycosylases—bacterial endonuclease VIII (Nei) and its mammalian homologs NEIL1 and NEIL2. Similarly to another homologous enzyme, bacterial formamidopyrimidine–DNA glycosylase, NEIL1 seems to use a processive mode to locate its targets. However, the processivity of Nei was notably lower, and NEIL2 exhibited almost fully distributive action on all types of substrates. Although one-dimensional diffusion is often regarded as a universal search mechanism, our results indicate that even proteins sharing a common fold may be quite different in the ways they locate their targets in DNA.

Noncatalytic Domains in DNA Glycosylases

Many proteins consist of two or more structural domains: separate parts that have a defined structure and function. For example, in enzymes, the catalytic activity is often localized in a core fragment, while other domains or disordered parts of the same protein participate in a number of regulatory processes. This situation is often observed in many DNA glycosylases, the proteins that remove damaged nucleobases thus initiating base excision DNA repair. This review covers the present knowledge about the functions and evolution of such noncatalytic parts in DNA glycosylases, mostly concerned with the human enzymes but also considering some unique members of this group coming from plants and prokaryotes.

DNA Deformation Exerted by Regulatory DNA-Binding Motifs in Human Alkyladenine DNA Glycosylase Promotes Base Flipping

Human alkyladenine DNA glycosylase (AAG) is a key enzyme that corrects a broad range of alkylated and deaminated nucleobases to maintain genomic integrity. When encountering the lesions, AAG adopts a base-flipping strategy to extrude the target base from the DNA duplex to its active site, thereby cleaving the glycosidic bond. Despite its functional importance, the detailed mechanism of such base extrusion and how AAG distinguishes the lesions from an excess of normal bases both remain elusive. Here, through the Markov state model constructed on extensive all-atom molecular dynamics simulations, we find that the alkylated nucleobase (N3-methyladenine, 3MeA) everts through the DNA major groove. Two key AAG motifs, the intercalation and E131-N146 motifs, play active roles in bending/pressing the DNA backbone and widening the DNA minor groove during 3MeA eversion. In particular, the intercalated residue Y162 is involved in buckling the target site at the early stage of 3MeA eversion. Our traveling-salesman based automated path searching algorithm further revealed that a non-target normal adenine tends to be trapped in an exo site near the active site, which however barely exists for a target base 3MeA. Collectively, these results suggest that the Markov state model combined with traveling-salesman based automated path searching acts as a promising approach for studying complex conformational changes of biomolecules and dissecting the elaborate mechanism of target recognition by this unique enzyme.

Free Energy Landscapes from SARS-CoV-2 Spike Glycoprotein Simulations Suggest that RBD Opening Can Be Modulated via Interactions in an Allosteric Pocket

The SARS-CoV-2 coronavirus is an enveloped, positive-sense single-stranded RNA virus that is responsible for the COVID-19 pandemic. The spike is a class I viral fusion glycoprotein that extends from the viral surface and is responsible for viral entry into the host cell and is the primary target of neutralizing antibodies. The receptor binding domain (RBD) of the spike samples multiple conformations in a compromise between evading immune recognition and searching for the host-cell surface receptor. Using atomistic simulations of the glycosylated wild-type spike in the closed and 1-up RBD conformations, we map the free energy landscape for RBD opening and identify interactions in an allosteric pocket that influence RBD dynamics. The results provide an explanation for experimental observation of increased antibody binding for a clinical variant with a substitution in this pocket. Our results also suggest the possibility of allosteric targeting of the RBD equilibrium to favor open states via binding of small molecules to the hinge pocket. In addition to potential value as experimental probes to quantify RBD conformational heterogeneity, small molecules that modulate the RBD equilibrium could help explore the relationship between RBD opening and S1 shedding.

Free Energy Landscapes from SARS-CoV-2 Spike Glycoprotein Simulations Suggest that RBD Opening Can Be Modulated via Interactions in an Allosteric Pocket

The SARS-CoV-2 coronavirus is an enveloped, positive-sense single-stranded RNA virus that is responsible for the COVID-19 pandemic. The spike is a class I viral fusion glycoprotein that extends from the viral surface and is responsible for viral entry into the host cell and is the primary target of neutralizing antibodies. The receptor binding domain (RBD) of the spike samples multiple conformations in a compromise between evading immune recognition and searching for the host-cell surface receptor. Using atomistic simulations of the glycosylated wild-type spike in the closed and 1-up RBD conformations, we map the free energy landscape for RBD opening and identify interactions in an allosteric pocket that influence RBD dynamics. The results provide an explanation for experimental observation of increased antibody binding for a clinical variant with a substitution in this pocket. Our results also suggest the possibility of allosteric targeting of the RBD equilibrium to favor open states via binding of small molecules to the hinge pocket. In addition to potential value as experimental probes to quantify RBD conformational heterogeneity, small molecules that modulate the RBD equilibrium could help explore the relationship between RBD opening and S1 shedding.

Recognition of a tandem lesion by DNA bacterial formamidopyrimidine glycosylases explored combining molecular dynamics and machine learning

The combination of several closely spaced DNA lesions, which can be induced by a single radical hit, constitutes a hallmark in the DNA damage landscape and radiation chemistry. The occurrence of such a tandem base lesion gives rise to a strong coupling with the double helix degrees of freedom and induces important structural deformations, in contrast to DNA strands containing a single oxidized nucleobase. Although such complex lesions are known to be refractory to repair by DNA glycosylases, there is still a lack of structural evidence to rationalize these phenomena. In this contribution, we explore, by numerical modeling and molecular simulations, the behavior of the bacterial glycosylase responsible for base excision repair (MutM), specialized in excising oxidatively-damaged defects such as 7,8-dihydro-8-oxoguanine (8-oxoG). The difference in lesion recognition between a simple damage and a tandem lesion featuring an additional abasic site is assessed at atomistic resolution owing to microsecond molecular dynamics simulations and machine learning postprocessing, allowing to extensively pinpoint crucial differences in the interaction patterns of the damaged bases. Our results reveal substantial changes in the interaction network surrounding the 8-oxoG upon addition of an adjacent abasic site, leading to the perturbation of the intercalation triad which is crucial for lesion recognition and processing. The recognition process might also be impacted by a more constrained MutM-DNA binding upon tandem damage, as shown by the machine learning post-processing. This work advocates for the use of such high throughput numerical simulations for exploring the complex combinatorial chemistry of tandem DNA lesions repair and more generally local multiple damaged sites of the utmost significance in radiation chemistry.

Catalytically Competent Conformation of the Active Site of Human 8-Oxoguanine-DNA Glycosylase

8-Oxoguanine-DNA N-glycosylase (OGG1) is a eukaryotic DNA repair enzyme responsible for the removal of 8-oxoguanine (oxoG), one of the most abundant oxidative DNA lesions. OGG1 catalyzes two successive reactions - N-glycosidic bond hydrolysis (glycosylase activity) and DNA strand cleavage on the 3'-side of the lesion by β-elimination (lyase activity). The enzyme also exhibits lyase activity with substrates containing apurinic/apyrimidinic (AP) sites (deoxyribose moieties lacking the nucleobase). OGG1 is highly specific for the base opposite the lesion, efficiently excising oxoG and cleaving AP sites located opposite to C, but not opposite to A. The activity is also profoundly decreased by amino acid changes that sterically interfere with oxoG binding in the active site of the enzyme after the lesion is everted from the DNA duplex. Earlier, the molecular dynamics approach was used to study the conformational dynamics of such human OGG1 mutants in complexes with the oxoG:C-containing substrate DNA, and the population density of certain conformers of two OGG1 catalytic residues, Lys249 and Asp268, was suggested to determine the enzyme activity. Here, we report the study of molecular dynamics of human OGG1 bound to the oxoG:A-containing DNA and OGG1 mutants bound to the AP:C-containing DNA. We showed that the enzyme low activity is associated with a decrease in the populations of Lys249 and Asp268 properly configured for catalysis. The experimentally measured rate constants for the OGG1 mutants show a good agreement with the models. We conclude that the enzymatic activity of OGG1 is determined majorly by the population density of the catalytically competent conformations of the active site residues Lys249 and Asp268.

Fast Implementation of the Nudged Elastic Band Method in AMBER

We present a fast implementation of the nudged elastic band (NEB) method into the particle mesh Ewald molecular dynamics module of the Amber software package for both central processing units (CPU) and graphics processing units (GPU). The accuracy of the new implementation has been validated for three cases: a conformational change of alanine dipeptide, the α-helix to β-sheet transition in polyalanine, and a large conformational transition in the human 8-oxoguanine-DNA glycosylase with DNA complex (OGG1-DNA). Timing benchmark tests were performed on the explicitly solvated OGG1-DNA system containing ∼50 000 atoms. The GPU-optimized implementation of NEB achieves a more than two orders of magnitude speedup compared with the previous CPU implementation performed with a two-core CPU processor. The speed and scalable features of this implementation will enable NEB applications on larger and more complex systems.

Critical Sites of DNA Backbone Integrity for Damaged Base Removal by Formamidopyrimidine-DNA Glycosylase

DNA glycosylases, the enzymes that initiate base excision DNA repair, recognize damaged bases through a series of precisely orchestrated movements. Most glycosylases sharply kink the DNA axis at the lesion site and extrude the target base from the DNA double helix into the enzyme's active site. Little attention has been paid so far to the role of the physical continuity of the DNA backbone in allowing the required conformational distortion. Here, we analyze base excision by formamidopyrimidine-DNA glycosylase (Fpg) from substrates keeping all phosphates but containing a nick within three nucleotides of the lesion in either DNA strand. Four phosphoester linkages at the damaged nucleotide and two nucleotides 3' to it were essential for Fpg activity, while the breakage of the others, even at the same critical phosphates, had no effect or even stimulated the reaction. Reduction of the likelihood of hydrogen bonding at the nicks by using dideoxynucleotides as their 3'-terminal groups was more detrimental for the activity. All phosphoester bonds in the complementary strand were dispensable for base excision, but nicks close to the orphaned nucleotide caused early termination of damaged strand cleavage. Elastic network analysis of Fpg-DNA structures showed that the vibrational motions of the critical phosphates are strongly correlated, in part due to the presence of the protein. Overall, our results suggest that mechanical forces propagating along the DNA backbone play a critical role in the correct conformational distortion of DNA by Fpg and possibly by other target base-everting DNA glycosylases.

Reading and Misreading 8-oxoguanine, a Paradigmatic Ambiguous Nucleobase

7,8-Dihydro-8-oxoguanine (oxoG) is the most abundant oxidative DNA lesion with dual coding properties. It forms both Watson–Crick (anti)oxoG:(anti)C and Hoogsteen (syn)oxoG:(anti)A base pairs without a significant distortion of a B-DNA helix. DNA polymerases bypass oxoG but the accuracy of nucleotide incorporation opposite the lesion varies depending on the polymerase-specific interactions with the templating oxoG and incoming nucleotides. High-fidelity replicative DNA polymerases read oxoG as a cognate base for A while treating oxoG:C as a mismatch. The mutagenic effects of oxoG in the cell are alleviated by specific systems for DNA repair and nucleotide pool sanitization, preventing mutagenesis from both direct DNA oxidation and oxodGMP incorporation. DNA translesion synthesis could provide an additional protective mechanism against oxoG mutagenesis in cells. Several human DNA polymerases of the X- and Y-families efficiently and accurately incorporate nucleotides opposite oxoG. In this review, we address the mutagenic potential of oxoG in cells and discuss the structural basis for oxoG bypass by different DNA polymerases and the mechanisms of the recognition of oxoG by DNA glycosylases and dNTP hydrolases.

Characterization of the Search Complex and Recognition Mechanism of the AlkD-DNA Glycosylase

DNA damage is a routine problem for cells, and pathways such as base excision repair have evolved to protect the genome by using DNA glycosylases to first recognize and excise lesions. The search mechanism of these enzymes is of particular interest due to the seemingly intractable problem of probing the billions of base pairs in the genome for potential damage. It has been hypothesized that glycosylases form multiple protein-DNA conformational states to efficiently search and recognize DNA lesions, ultimately only flipping out the damaged substrate into the active site. A unique DNA glycosylase, the Bacillus cereus AlkD enzyme, has been shown to excise damaged DNA without flipping the nucleobase into a protein binding pocket following lesion recognition. Here, we use microsecond-scale all-atom molecular dynamics simulations to characterize the AlkD recognition mechanism, putting it in perspective with other DNA glycosylases. We first identify and describe two distinct enzyme-DNA conformations of AlkD: the search complex (SC) and excision complex (EC). The SC is distinguished by the linearity of DNA, changes in four helical parameters in the vicinity of the lesion, and changes in distance between active site residues and the DNA. Free DNA simulations are used to demonstrate that the DNA structural deviations and increased active site interactions present in the EC are initiated by the recognition of a methylation-induced signal in the rises both 5' to the methylation and opposing this base. Our results support the hypothesis that subtle geometric distortions in DNA are recognized by AlkD and are consequently probed to initiate concerted protein and DNA conformational changes which prime excise without additional intermediate states. This mechanism is shown to be consistent among the three methylated DNA sequences that have been crystallized bound to AlkD.

Base-Independent DNA Base-Excision Repair of 8-Oxoguanine

Living organisms protect their genome from gene mutation by excising damaged DNA bases. Here, 8-oxoguanine (8OG) is one of the most abundant DNA lesions. In bacteria the base excision is catalyzed by the enzyme formamidopyrimidine-DNA- glycosylase (Fpg), for which two different orientations of 8OG binding into the active site of Fpg have been proposed: syn- and anti-conformation. Here, we present a new ribose-protonated repair mechanism for 8OG that is base-independent and can excise 8OG in both conformations. Using high-level QM/MM calculations with up to 588/573 atoms in the QM sphere, the activation barrier is computed in excellent agreement with the experimentally measured value. Since the excised base itself is not directly involved in the mechanism, this implies that lesion discrimination does not occur within the active site of the enzyme.

Aberrant base excision repair pathway of oxidatively damaged DNA: Implications for degenerative diseases

In cellular organisms composition of DNA is constrained to only four nucleobases A, G, T and C, except for minor DNA base modifications such as methylation which serves for defence against foreign DNA or gene expression regulation. Interestingly, this severe evolutionary constraint among other things demands DNA repair systems to discriminate between regular and modified bases. DNA glycosylases specifically recognize and excise damaged bases among vast majority of regular bases in the base excision repair (BER) pathway. However, the mismatched base pairs in DNA can occur from a spontaneous conversion of 5-methylcytosine to thymine and DNA polymerase errors during replication. To counteract these mutagenic threats to genome stability, cells evolved special DNA repair systems that target the non-damaged DNA strand in a duplex to remove mismatched regular DNA bases. Mismatch-specific adenine- and thymine-DNA glycosylases (MutY/MUTYH and TDG/MBD4, respectively) initiated BER and mismatch repair (MMR) pathways can recognize and remove normal DNA bases in mismatched DNA duplexes. Importantly, in DNA repair deficient cells bacterial MutY, human TDG and mammalian MMR can act in the aberrant manner: MutY and TDG removes adenine and thymine opposite misincorporated 8-oxoguanine and damaged adenine, respectively, whereas MMR removes thymine opposite to O⁶-methylguanine. These unusual activities lead either to mutations or futile DNA repair, thus indicating that the DNA repair pathways which target non-damaged DNA strand can act in aberrant manner and introduce genome instability in the presence of unrepaired DNA lesions. Evidences accumulated showing that in addition to the accumulation of oxidatively damaged DNA in cells, the aberrant DNA repair can also contribute to cancer, brain disorders and premature senescence. For example, the aberrant BER and MMR pathways for oxidized guanine residues can lead to trinucleotide expansion that underlies Huntington's disease, a severe hereditary neurodegenerative syndrome. This review summarises the present knowledge about the aberrant DNA repair pathways for oxidized base modifications and their possible role in age-related diseases.

Repair of 8-oxo-7,8-dihydroguanine in prokaryotic and eukaryotic cells: Properties and biological roles of the Fpg and OGG1 DNA N-glycosylases

Oxidatively damaged DNA results from the attack of sugar and base moieties by reactive oxygen species (ROS), which are formed as byproducts of normal cell metabolism and during exposure to endogenous or exogenous chemical or physical agents. Guanine, having the lowest redox potential, is the DNA base the most susceptible to oxidation, yielding products such as 8-oxo-7,8-dihydroguanine (8-oxoG) and 2-6-diamino-4-hydroxy-5-formamidopyrimidine (FapyG). In DNA, 8-oxoG was shown to be mutagenic yielding GC to TA transversions upon incorporation of dAMP opposite this lesion by replicative DNA polymerases. In prokaryotic and eukaryotic cells, 8-oxoG is primarily repaired by the base excision repair pathway (BER) initiated by a DNA N-glycosylase, Fpg and OGG1, respectively. In Escherichia coli, Fpg cooperates with MutY and MutT to prevent 8-oxoG-induced mutations, the "GO-repair system". In Saccharomyces cerevisiae, OGG1 cooperates with nucleotide excision repair (NER), mismatch repair (MMR), post-replication repair (PRR) and DNA polymerase η to prevent mutagenesis. Human and mouse cells mobilize all these pathways using OGG1, MUTYH (MutY-homolog also known as MYH), MTH1 (MutT-homolog also known as NUDT1), NER, MMR, NEILs and DNA polymerases η and λ, to prevent 8-oxoG-induced mutations. In fact, mice deficient in both OGG1 and MUTYH develop cancer in different organs at adult age, which points to the critical impact of 8-oxoG repair on genetic stability in mammals. In this review, we will focus on Fpg and OGG1 proteins, their biochemical and structural properties as well as their biological roles. Other DNA N-glycosylases able to release 8-oxoG from damaged DNA in various organisms will be discussed. Finally, we will report on the role of OGG1 in human disease and the possible use of 8-oxoG DNA N-glycosylases as therapeutic targets.

DNA Deformation-Coupled Recognition of 8-Oxoguanine: Conformational Kinetic Gating in Human DNA Glycosylase

8-Oxoguanine (8-oxoG), a mutagenic DNA lesion generated under oxidative stress, differs from its precursor guanine by only two substitutions (O⁸ and H⁷). Human 8-oxoguanine glycosylase 1 (OGG1) can locate and remove 8-oxoG through extrusion and excision. To date, it remains unclear how OGG1 efficiently distinguishes 8-oxoG from a large excess of undamaged DNA bases. We recently showed that formamidopyrimidine-DNA glycosylase (Fpg), a bacterial functional analog of OGG1, can selectively facilitate eversion of oxoG by stabilizing several intermediate states, and it is intriguing whether OGG1 also employs a similar mechanism in lesion recognition. Here, we use molecular dynamics simulations to explore the mechanism by which OGG1 discriminates between 8-oxoG and guanine along the base-eversion pathway. The MD results suggest an important role for kinking of the DNA by the glycosylase, which positions DNA phosphates in a way that assists lesion recognition during base eversion. The computational predictions were validated through experimental enzyme assays on phosphorothioate substrate analogs. Our simulations suggest that OGG1 distinguishes between 8-oxoG and G using their chemical dissimilarities not only at the active site but also at earlier stages during base eversion, and this mechanism is at least partially conserved in Fpg despite a lack of structural homology. The similarity also suggests that lesion recognition through multiple gating steps may be a common theme in DNA repair. Our results provide new insight into how enzymes can exploit kinetics and DNA conformational changes to probe the chemical modifications present in DNA lesions.