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

Back-Up Base Excision DNA Repair in Human Cells Deficient in the Major AP Endonuclease, APE1

Apurinic/apyrimidinic (AP) sites are abundant DNA lesions generated both by spontaneous base loss and as intermediates of base excision DNA repair. In human cells, they are normally repaired by an essential AP endonuclease, APE1, encoded by the APEX1 gene. Other enzymes can cleave AP sites by either hydrolysis or β-elimination in vitro, but it is not clear whether they provide the second line of defense in living cells. Here, we studied AP site repairs in APEX1 knockout derivatives of HEK293FT cells using a reporter system based on transcriptional mutagenesis in the enhanced green fluorescent protein gene. Despite an apparent lack of AP site-processing activity in vitro, the cells efficiently repaired the tetrahydrofuran AP site analog resistant to β-elimination. This ability persisted even when the second AP endonuclease homolog, APE2, was also knocked out. Moreover, APEX1 null cells were able to repair uracil, a DNA lesion that is removed via the formation of an AP site. If AP site hydrolysis was chemically blocked, the uracil repair required the presence of NTHL1, an enzyme that catalyzes β-elimination. Our results suggest that human cells possess at least two back-up AP site repair pathways, one of which is NTHL1-dependent.

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.

Fluorescent adenine analogues with ESPT characteristic utilized for real-time detecting DNA adduct

The 8-oxo-7,8-dihydro-2'-deoxyguanine (8-oxoG) is the representative damaged nucleoside that may increase the risk of developing diseases. Accordingly, the selective detection of 8-oxoG in DNA with minimal disturbance to the native structure is important to have an in-depth understanding of the formation mechanism and becomes an attractive tool for genomic research. To identify the DNA adduct in real-time efficiently, a series of quasi-intrinsic optical probes are performed based on the natural adenine, which has preference to form a stable base pair with 8-oxoG in the syn conformation. The calculations revealed that the A-analogues in solution could bring red-shifted absorption spectra and bright photoluminescence arisen from the additional π-conjugation by means of fluorophore modification and the ring expansion. Especially, A1 possesses large Stokes shifts and the highest fluorescence intensity in emission, which is proposed as the biosensor to monitor the optical changes in the presence and absence of the considered 8-oxoG. It is found that the fluorescence is insensitive to base pairing with thymine, while the excited state intermolecular proton transfer (ESPT) induced efficient fluorescence quenching is observed upon pairing with the 8-oxoG. To evaluate the direct usefulness of the bright adenine analogues in biological environment, we further examined the influences of linking deoxyribose on the absorption and emission, which are consistent with the experimental data.

Markov state models elucidate the stability of DNA influenced by the chiral 5S-Tg base

The static and dynamic structures of DNA duplexes affected by 5S-Tg (Tg, Thymine glycol) epimers were studied using MD simulations and Markov State Models (MSMs) analysis. The results show that the 5S,6S-Tg base caused little perturbation to the helix, and the base-flipping barrier was determined to be 4.4 kcal mol-1 through the use of enhanced sampling meta-eABF calculations, comparable to 5.4 kcal mol-1 of the corresponding thymine flipping. Two conformations with the different hydrogen bond structures between 5S,6R-Tg and A19 were identified in several independent MD trajectories. The 5S,6R-Tg:O6HO6•••N1:A19 hydrogen bond is present in the high-energy conformation displaying a clear helical distortion, and near barrier-free Tg base flipping. The low-energy conformation always maintains Watson-Crick base pairing between 5S,6R-Tg and A19, and 5S-Tg base flipping is accompanied by a small barrier of ca. 2.0 KBT (T = 298 K). The same conformations are observed in the MSMs analysis. Moreover, the transition path and metastable structures of the damaged base flipping are for the first time verified through MSMs analysis. The data clearly show that the epimers have completely different influence on the stability of the DNA duplex, thus implying different enzymatic mechanisms for DNA repair.

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.

Dynamics of 5R-Tg Base Flipping in DNA Duplexes Based on Simulations─Agreement with Experiments and Beyond

Damaged or mismatched DNA bases are normally thought to be able to flip out of the helical stack, providing enzymes with access to the faulty genetic information otherwise hidden inside the helix. Thymine glycol (Tg) is one of the most common products of nucleic acid damage. However, the static and dynamic structures of DNA duplexes affected by 5R-Tg epimers are still not clearly understood, including the ability of these to undergo spontaneous base flipping. Structural effects of the 5R-Tg epimers on the duplex DNA are herein studied using molecular dynamics together with reliable DFT based calculations. In comparison with the corresponding intact DNA, the cis-5R,6S-Tg epimer base causes little perturbation to the duplex DNA, and a barrier of 4.9 kcal mol^–1 is obtained by meta-eABF for cis-5R,6S-Tg base flipping out of the duplex DNA, comparable to the 5.4 kcal mol^–1 obtained for the corresponding thymine flipping in intact DNA. For the trans-5R,6R-Tg epimer, three stable local structures were identified, of which the most stable disrupts the Watson–Crick hydrogen-bonded G5/C20 base pair, leading to conformational distortion of the duplex. Interestingly, the relative barrier height of the 5R-Tg flipping is only 1.0 kcal mol^–1 for one of these trans-5R,6R-Tg epimers. Water bridge interactions were identified to be essential for 5R-Tg flipping. The study clearly demonstrates the occurrence of partial trans-5R,6R-Tg epimer flipping in solution.

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.

A perspective on the molecular simulation of DNA from structural and functional aspects

As genetic material, DNA not only carries genetic information by sequence, but also affects biological functions ranging from base modification to replication, transcription and gene regulation through its structural and dynamic properties and variations. The motion and structural properties of DNA involved in related biological processes are also multi-scale, ranging from single base flipping to local DNA deformation, TF binding, G-quadruplex and i-motif formation, TAD establishment, compartmentalization and even chromosome territory formation, just to name a few. The sequence-dependent physical properties of DNA play vital role in all these events, and thus it is interesting to examine how simple sequence information affects DNA and the formation of the chromatin structure in these different hierarchical orders. Accordingly, molecular simulations can provide atomistic details of interactions and conformational dynamics involved in different biological processes of DNA, including those inaccessible by current experimental methods. In this perspective, which is mainly based on our recent studies, we provide a brief overview of the atomistic simulations on how the hierarchical structure and dynamics of DNA can be influenced by its sequences, base modifications, environmental factors and protein binding in the context of the protein-DNA interactions, gene regulation and structural organization of chromatin. We try to connect the DNA sequence, the hierarchical structures of DNA and gene regulation.

Reaction Coordinate and Thermodynamics of Base Flipping in RNA

Base flipping is a key biophysical event involved in recognition of various ligands by ribonucleic acid (RNA) molecules. However, the mechanism of base flipping in RNA remains poorly understood, in part due to the lack of atomistic details on complex rearrangements in neighboring bases. In this work, we applied transition path sampling (TPS) methods to study base flipping in a double-stranded RNA (dsRNA) molecule that is known to interact with RNA-editing enzymes through this mechanism. We obtained an ensemble of 1000 transition trajectories to describe the base-flipping process. We used the likelihood maximization method to determine the refined reaction coordinate (RC) consisting of two collective variables (CVs), a distance and a dihedral angle between nucleotides that form stacking interactions with the flipping base. The free energy profile projected along the refined RC revealed three minima, two corresponding to the initial and final states and one for a metastable state. We suggest that the metastable state likely represents a wobbled conformation of nucleobases observed in NMR studies that is often characterized as the flipped state. The analyses of reactive trajectories further revealed that the base flipping is coupled to a global conformational change in a stem-loop of dsRNA.

Molecular dynamics approach to identification of new OGG1 cancer-associated somatic variants with impaired activity

DNA of living cells is always exposed to damaging factors. To counteract the consequences of DNA lesions, cells have evolved several DNA repair systems, among which base excision repair is one of the most important systems. Many currently used antitumor drugs act by damaging DNA, and DNA repair often interferes with chemotherapy and radiotherapy in cancer cells. Tumors are usually extremely genetically heterogeneous, often bearing mutations in DNA repair genes. Thus, knowledge of the functionality of cancer-related variants of proteins involved in DNA damage response and repair is of great interest for personalization of cancer therapy. Although computational methods to predict the variant functionality have attracted much attention, at present, they are mostly based on sequence conservation and make little use of modern capabilities in computational analysis of 3D protein structures. We have used molecular dynamics (MD) to model the structures of 20 clinically observed variants of a DNA repair enzyme, 8-oxoguanine DNA glycosylase. In parallel, we have experimentally characterized the activity, thermostability, and DNA binding in a subset of these mutant proteins. Among the analyzed variants of 8-oxoguanine DNA glycosylase, three (I145M, G202C, and V267M) were significantly functionally impaired and were successfully predicted by MD. Alone or in combination with sequence-based methods, MD may be an important functional prediction tool for cancer-related protein variants of unknown significance.

Lost in the Crowd: How Does Human 8-Oxoguanine DNA Glycosylase 1 (OGG1) Find 8-Oxoguanine in the Genome?

The most frequent DNA lesion resulting from an oxidative stress is 7,8-dihydro-8-oxoguanine (8-oxoG). 8-oxoG is a premutagenic base modification due to its capacity to pair with adenine. Thus, the repair of 8-oxoG is critical for the preservation of the genetic information. Nowadays, 8-oxoG is also considered as an oxidative stress-sensor with a putative role in transcription regulation. In mammalian cells, the modified base is excised by the 8-oxoguanine DNA glycosylase (OGG1), initiating the base excision repair (BER) pathway. OGG1 confronts the massive challenge that is finding rare occurrences of 8-oxoG among a million-fold excess of normal guanines. Here, we review the current knowledge on the search and discrimination mechanisms employed by OGG1 to find its substrate in the genome. While there is considerable data from in vitro experiments, much less is known on how OGG1 is recruited to chromatin and scans the genome within the cellular nucleus. Based on what is known of the strategies used by proteins searching for rare genomic targets, we discuss the possible scenarios allowing the efficient detection of 8-oxoG by OGG1.

Inhibitors of DNA Glycosylases as Prospective Drugs

DNA glycosylases are enzymes that initiate the base excision repair pathway, a major biochemical process that protects the genomes of all living organisms from intrinsically and environmentally inflicted damage. Recently, base excision repair inhibition proved to be a viable strategy for the therapy of tumors that have lost alternative repair pathways, such as BRCA-deficient cancers sensitive to poly(ADP-ribose)polymerase inhibition. However, drugs targeting DNA glycosylases are still in development and so far have not advanced to clinical trials. In this review, we cover the attempts to validate DNA glycosylases as suitable targets for inhibition in the pharmacological treatment of cancer, neurodegenerative diseases, chronic inflammation, bacterial and viral infections. We discuss the glycosylase inhibitors described so far and survey the advances in the assays for DNA glycosylase reactions that may be used to screen pharmacological libraries for new active compounds.

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.

Quantifying the stability of oxidatively damaged DNA by single-molecule DNA stretching

One of the most common DNA lesions is created when reactive oxygen alters guanine. 8-oxo-guanine may bind in the anti-conformation with an opposing cytosine or in the syn-conformation with an opposing adenine paired by transversion, and both conformations may alter DNA stability. Here we use optical tweezers to measure the stability of DNA hairpins containing 8-oxoguanine (8oxoG) lesions, comparing the results to predictive models of base-pair energies in the absence of the lesion. Contrasted with either a canonical guanine-cytosine or adenine-thymine pair, an 8oxoG-cytosine base pair shows significant destabilization of several k_BT. The magnitude of destabilization is comparable to guanine-thymine ‘wobble’ and cytosine-thymine mismatches. Furthermore, the measured energy of 8oxoG-adenine corresponds to theoretical predictions for guanine-adenine pairs, indicating that oxidative damage does not further destabilize this mismatch in our experiments, in contrast to some previous observations. These results support the hypothesis that oxidative damage to guanine subtly alters the direction of the guanine dipole, base stacking interactions, the local backbone conformation, and the hydration of the modified base. This localized destabilization under stress provides additional support for proposed mechanisms of enzyme repair.

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.

Electrophilic oligodeoxynucleotide synthesis using dM-Dmoc for amino protection

Solid-phase synthesis of electrophilic oligodeoxynucleotides (ODNs) was achieved using dimethyl-Dmoc (dM-Dmoc) as amino protecting group. Due to the high steric hindrance of the 2-(propan-2-ylidene)-1,3-dithiane side product from deprotection, the use of excess nucleophilic scavengers such as aniline to prevent Michael addition of the side product to the deprotected ODN during ODN cleavage and deprotection was no longer needed. The improved technology was demonstrated by the synthesis and characterization of five ODNs including three modified ones. The modified ODNs contained the electrophilic groups ethyl ester, α-chloroamide, and thioester. Using the technology, the sensitive groups can be installed at any location within the ODN sequences without using any sequence- or functionality-specific conditions and procedures.

Sensitive detection of formamidopyrimidine-DNA glycosylase activity based on target-induced self-primed rolling circle amplification and magnetic nanoprobes

We developed a novel approach to determine formamidopyrimidine DNA glycosylase (FPG) activity by taking advantage of target-induced self-primed rolling circle amplification (RCA) and magnetic nanoprobes. Herein, a unique nick (8-oxoguanine, 8-oxoG) was positioned in duplex DNA containing P-circle and P1, which together serve as a FPG substrate, RCA template, and RCA primer probe. The presence of FPG specifically binds 8-oxoG and cleaves the P-circle into two parts, producing 5'-phosphoryl termini. A phosphodiester bond between the 5'-phosphoryl and 3'-hydroxyl termini was formed with the addition of T4 DNA ligase, producing an unnicked circular strand. Using the unnicked strand as the RCA template, the P1 hybridized with the circle probe as a primer will trigger the RCA process. The RCA reaction produces amounts of long tandem-repeat DNA tiles with multiple recognizing regions for the FAM modified DNA probes (FP) and biotin-modified DNA probes (BP). With the streptavidin-biotin interaction, the BPs and FPs can be easily immobilized on the surface of streptavidin-modified magnetic microbeads (MBs). Due to the RCA enhanced and highly-concentrated fluorescence accumulation on the MBs, an ultralow detection limit of 1.033 U mL-1 for FPG was obtained. Combined with the high tolerance capability of human blood serum owing to magnetic isolation, the FPG assays in human blood serum were also obtained using fluorescence and confocal laser scanning microscopy. These results indicate that this robust self-primed RCA combined with magnetic nanoprobes is an excellent candidate for quantitatively monitoring the FPG activity responsible for DNA oxidative damage-related clinical diagnosis and therapy.

Structural interpretation of the 31P NMR chemical shifts in thiophosphate and phosphate: key effects due to spin–orbit and explicit solvent

Structural interpretation of the ³¹P NMR shifts measured in different molecules including thiophosphate or phosphate group was obtained by means of theoretical calculations including the effects of geometry, molecular dynamics, solvent, relativistic effects and the effect of NMR reference.

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.

Role of Conformational Fluctuations of Protein toward Methylation in DNA by Cytosine-5-methyltransferase

Methylation of cytosine is the common epigenetic modification in genomes ranging from bacteria to mammals, and aberrant methylation leads to human diseases including cancer. Recognition of a cognate DNA sequence by DNA methyltransferases and flipping of a target base into the enzyme active site pocket are the key steps in DNA methylation. Using molecular dynamics simulations and enhanced sampling techniques here we elucidate the role of conformational fluctuations of protein and active or passive involvement of protein elements that mediate base flipping and formation of the closed catalytic complex. The free energy profiles for the flipping of target cytosine into the enzyme active site support the major groove base eversion pathway; and the results show that the closed state of enzyme increases the free energy barrier, whereas the open state reduces it. We found that the interactions of the key loop residues of protein with cognate DNA altered the protein motions, and modulation of protein fluctuations relates to the closed catalytic complex formation. Methylation of cytosine in the active site of the closed complex destabilizes the interactions of catalytic loop residues with cognate DNA and reduces the stability of the closed state. Our study provides microscopic insights on the base flipping mechanism coupled with enzyme's loop motions and provides evidence for the role of conformational fluctuations of protein in the enzyme-catalyzed DNA processing mechanism.

Structural Insight into the Discrimination between 8-Oxoguanine Glycosidic Conformers by DNA Repair Enzymes: A Molecular Dynamics Study of Human Oxoguanine Glycosylase 1 and Formamidopyrimidine-DNA Glycosylase

hOgg1 and FPG are the primary DNA repair enzymes responsible for removing the major guanine (G) oxidative product, namely, 7,8-dihydro-8-oxoguanine (OG), in humans and bacteria, respectively. While natural G adopts the anti conformation and forms a Watson-Crick pair with cytosine (C), OG can also adopt the syn conformation and form a Hoogsteen pair with adenine (A). hOgg1 removes OG paired with C but is inactive toward the OG:A pair. In contrast, FPG removes OG from OG:C pairs and also exhibits appreciable (although diminished) activity toward OG:A pairs. As a first step toward understanding this difference in activity, we have employed molecular dynamics simulations to examine how the anti and syn conformers of OG are accommodated in the hOgg1 and FPG active sites. When anti-OG is bound, hOgg1 active site residues are properly aligned to initiate catalytic base departure, while geometrical parameters required for the catalytic reaction are not conserved for syn-OG. On the other hand, the FPG catalytic residues are suitably aligned for both OG conformers, with anti-OG being more favorably bound. Thus, our data suggests that the differential ability of hOgg1 and FPG to accommodate the anti- and syn-OG glycosidic conformations is an important factor that contributes to the relative experimental excision rates. Nevertheless, the positions of the nucleophiles with respect to the lesion in the active sites suggest that the reactant complex is poised to initiate catalysis through a similar mechanism for both repair enzymes and supports a recently proposed mechanism in which sugar-ring opening precedes nucleoside deglycosylation.