Evaluating the Substrate Selectivity of Alkyladenine DNA Glycosylase: The Synergistic Interplay of Active Site Flexibility and Water Reorganization

Investigation of the Flipping Dynamics of 1, N6-Ethenoadenine in Alkyladenine DNA Glycosylase

Alkyladenine DNA glycosylase (AAG) is an essential enzyme responsible for maintaining genome integrity by repairing several DNA lesions damaged by alkylation or deamination. Understanding how it can recognize and excise the lesions thus lays the foundation for therapeutic treatment against lesion-associated diseases or cancers. However, the molecular details of how the lesion can be distinguished from the matched base by AAG and how it enters the cleavage site, ready for excision, are not fully elucidated. In this study, we have revealed the molecular details of the flipping dynamics of 1, N6-ethenoadenine (εA) not only in the form of free double-stranded DNA (dsDNA) but also in the form of the AAG-dsDNA complex. Our MD simulations and PMF calculations have shown that the flipping of εA and dA is thermodynamically disfavored in the free dsDNA, even though εA has a lower flipping energy barrier than dA. By sharp contrast, the flipping of εA is thermodynamically favored in AAG with an obvious free energy drop, while dA is equally stabilized before and after the flipping. Moreover, a comparison of the PMFs in the forms of free dsDNA and the AAG-dsDNA complex has pinpointed the role of AAG in discriminating εA against dA and facilitating the flipping of εA. Besides, the flipping process is simulated along the major and minor grooves, and our results have additionally demonstrated that the flipping is not directional in the free dsDNA while flipping along the major groove is kinetically more favorable than the minor groove in the AAG-dsDNA complex. Overall, our study has offered molecular insights into the flipping dynamics of εA and revealed its discrimination mechanism by AAG, which is expected to guide further enzyme engineering for therapeutic applications.

Development of a QM/MM(ABEEM) method combined with a polarizable force field to investigate the excision reaction mechanism of damaged thymine

This paper focuses on the development of a quantum mechanics/molecular mechanics method using the ABEEM polarizable force field (QM/MM(ABEEM) method) to investigate the excision reaction mechanism of damaged thymine. This method does not simply combine the QM method with the polarizable force field. A valence electronegativity piecewise function with the distance between atoms as a variable is introduced to describe the atomic partial charges, and changes greatly during the reaction process. At the same time, the charge transfer effect is treated using the condition of local charge conservation. Compared with the traditional QM/MM method, the QM/MM(ABEEM) method can more accurately simulate the polarization effect and charge transfer effect in the reaction process. Focusing on the controversial problems of the excision of damaged bases, six reaction pathways were designed for monofunctional and difunctional deglycosylation of neutral bases and protonated bases. The results show that the QM/MM(ABEEM) method accurately simulates the polarization effect, charge transfer effect, activation energy and other properties of the reaction process. The process in which the active residue Asp activates the nucleophile H₂O to attack the protonated base is the preferred path. The average activation energy and free activation energy of the protonated base are 7.00-14.00 kcal mol^-1 lower than that of the neutral base. The study in this paper is helpful to understand the mechanism of repair enzymes in repairing bases.

Development of QM/MM (ABEEM polarizable force field) method to simulate the excision reaction mechanism of damaged cytosine

DNA damages are regarded as having harmful effects on cell. The base excision repair mechanism combats these effects by removing damaged bases. The deglycosylation mechanism of excising damaged bases by DNA glycosylase and the state of the leaving base have been controversial. The enzymatic reaction of DNA glycosylase to remove the damaged bases involves not only the formation and breaking of chemical bonds, but also complex polarization effect and charge transfer, which cannot be accurately simulated by the QM/MM method combined with the fixed charge force field. This work has developed the ABEEM fluctuating polarizable force field combining with the QM method, that is (QM/MM[ABEEM]), to accurately simulate the proton transfer, charge transfer and the charge distribution. The piecewise function is used as the valence-state electronegativity in the QM/MM (ABEEM) to realize the accurate fitting of the charge distribution in reaction. And the charge transfer is accurately simulated by the local charge conservation conditions. Four deglycosylation mechanisms including the monofunctional and difunctional mechanisms of four neutral and protonated cytosine derivatives are explored. It is confirmed that the monofunctional mechanism of Asp-activated nucleophile water is a better deglycosylation mechanism and the base is protonated before the reaction occurs. Protonization of the base reduced the activation energy by 10.00-17.00 kcal/mol. Asp provides the necessary charge for the reaction, and DNA glycosylase preferentially cleaves ɛC. This work provides a theoretical basis for the research of excising damaged bases by DNA glycosylase.

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.

Water-Mediated Oxidation of Guanine by a Repair Enzyme: Simulation Using the ABEEM Polarizable Force Field

The recognition mechanism of oxidative damage in organisms has long been a research hotspot. Water is an important medium in the recognition process, but its specific role remains unknown. There is a need to develop a suitable force field that can adequately describe the electrostatic, hydrogen bond, and other interactions among the molecules in the complex system of the repair enzyme and oxidized base. The developing ABEEM polarizable force field (PFF) has been used to simulate the repaired enzyme hOGG1 and oxidized DNA (PDB ID: 1EBM) in a biological environment, and the corresponding results are better than those of the fixed-charge force fields OPLS/AA and AMBER OL15. 8-Oxo-G is recognized by Gln315 of hOGG1 mainly through hydrogen bonds mediated by continuous exchange of 2 water molecules. Phe319 and Cys253 are stacked on both sides of the π planes of bases to form sandwich structures. The charge polarization effect gives an important signal to drive the exchange of water molecules and maintains the recognition of oxidation bases by enzymes. The mediated main water molecule A and mediated auxiliary water molecule B together pull Gln315 to recognize 8-oxo-G by hydrogen bond interactions. Then, the charge polarization signal of solvent water molecule C with a large absolute charge causes the absolute charge of O atoms in water molecule A or B to increase by approximately 0.2 e, and water molecule A or B leaves Gln315 and 8-oxo-G. The other water molecule and water molecule C synergistically recognize 8-oxo-G with Gln315. Even though the water molecules between Gln315 and 8-oxo-G are removed, the MD simulation results show that water molecules appear between Gln315 and 8-oxo-G in a very short time (<2 ps). The dwell time of each water molecule is approximately 60 ps. The radial distribution function and dwell time support the correctness of the above mechanism. These polarization effects and hydrogen bonding interactions cannot be simulated by a fixed-charge force field.

Altered APE1 activity on abasic ribonucleotides is mediated by changes in the nucleoside sugar pucker

Ribonucleotides (rNTPs) are predicted to be incorporated into the genome at a rate of up to 3 million times per cell division, making rNTPs the most common non-standard nucleotide in the human genome. Typically, misinserted ribonucleotides are repaired by the ribonucleotide excision repair (RER) pathway, which is initiated by RNase H2 cleavage. However, rNTPs are susceptible to spontaneous depurination generating abasic ribonucleotides (rAPs), which are unable to be processed by RNase H2. Additionally, rAPs have been found in nascent RNA and coupled to R-loops. Recent work identified that base excision repair (BER) protein AP-Endonuclease 1 (APE1) is responsible for the initial processing of rAPs embedded in DNA and in R-loops. APE1 is a well characterized AP endonuclease that cleaves 5' of abasic sites, but its ability to cleave at rAPs remains poorly understood. Here, we utilize enzyme kinetics, X-ray crystallography, and molecular dynamics simulations to provide insight into rAP processing by APE1. Enzyme kinetics were used to determine pre-steady-state rates of APE1 cleavage on DNA substrates containing rAP, revealing a decrease in activity compared to cleavage at a canonical deoxy-AP substrate. Using X-ray crystallography, we identified novel contacts between the rAP and the APE1 active site. We demonstrate that the rAP sugar pucker is accommodated in the active site in a C3'-endo conformation, influencing its position and contributing to a decrease in activity compared to the deoxy-AP site. Together, this work provides molecular level insights into rAP processing by APE1 and advances our understanding of ribonucleotide processing within genomic DNA.

Insights into the Direct Oxidative Repair of Etheno Lesions: MD and QM/MM Study on the Substrate Scope of ALKBH2 and AlkB

E. coli AlkB and human ALKBH2 belong to the AlkB family enzymes, which contain several α-ketoglutarate (α-KG)/Fe(II)-dependent dioxygenases that repair alkylated DNA. Specifically, the AlkB enzymes catalyze decarboxylation of α-KG to generate a high-valent Fe(IV)-oxo species that oxidizes alkyl groups on DNA adducts. AlkB and ALKBH2 have been reported to differentially repair select etheno adducts, with preferences for 1,N⁶-ethenoadenine (1,N⁶-εA) and 3,N⁴-ethenocytosine (3,N⁴-εC) over 1,N²-ethenoguanine (1,N²-εG). However, N²,3-ethenoguanine (N²,3-εG), the most common etheno adduct, is not repaired by the AlkB enzymes. Unfortunately, a structural understanding of the differential activity of E. coli AlkB and human ALKBH2 is lacking due to challenges acquiring atomistic details for a range of substrates using experiments. This study uses both molecular dynamics (MD) simulations and ONIOM(QM:MM) calculations to determine how the active site changes upon binding each etheno adduct and characterizes the corresponding catalytic impacts. Our data reveal that the preferred etheno substrates (1,N⁶-εA and 3,N⁴-εC) form favorable interactions with catalytic residues that situate the lesion near the Fe(IV)-oxo species and permit efficient oxidation. In contrast, although the damage remains correctly aligned with respect to the Fe(IV)-oxo moiety, repair of 1,N²-εG is mitigated by increased solvation of the active site and a larger distance between Fe(IV)-oxo and the aberrant carbons. Binding of non-substrate N²,3-εG in the active site disrupts key DNA-enzyme interactions, and positions the aberrant carbon atoms even further from the Fe(IV)-oxo species, leading to prohibitively high barriers for oxidative catalysis. Overall, our calculations provide the first structural insight required to rationalize the experimentally-reported substrate specificities of AlkB and ALKBH2 and thereby highlight the roles of several active site residues in the repair of etheno adducts that directly correlates with available experimental data. These proposed catalytic strategies can likely be generalized to other α-KG/Fe(II)-dependent dioxygenases that play similar critical biological roles, including epigenetic and post-translational regulation.

RNA abasic sites in yeast and human cells

RNA abasic sites and the mechanisms involved in their regulation are mostly unknown; in contrast, DNA abasic sites are well-studied. We found surprisingly that, in yeast and human cells, RNA abasic sites are prevalent. When a base is lost from RNA, the remaining ribose is found as a closed-ring or an open-ring sugar with a reactive C1' aldehyde group. Using primary amine-based reagents that react with the aldehyde group, we uncovered evidence for abasic sites in nascent RNA, messenger RNA, and ribosomal RNA from yeast and human cells. Mass spectroscopic analysis confirmed the presence of RNA abasic sites. The RNA abasic sites were found to be coupled to R-loops. We show that human methylpurine DNA glycosylase cleaves N-glycosidic bonds on RNA and that human apurinic/apyrimidinic endonuclease 1 incises RNA abasic sites in RNA-DNA hybrids. Our results reveal that, in yeast and human cells, there are RNA abasic sites, and we identify a glycosylase that generates these sites and an AP endonuclease that processes them.

Heterogeneous Dynamical Environment at the Interface of a Protein-DNA Complex

Binding between protein and DNA is an essential process to regulate different biological activities. Two puzzling questions in protein-DNA recognition are (i) how the protein's binding domain identifies the DNA sequence in an aqueous solution and (ii) how the formation of the complex alters the dynamical environment around it. In this work, we present results obtained from molecular dynamics simulations of the N-terminal α-helical domain of the λ-repressor protein (in dimeric form) bound to the corresponding operator DNA. Effects of formation of the complex in modifying the microscopic dynamics of water as well as the kinetics of hydrogen bonds at the interface have been explored. Locally heterogeneous restricted water motions at the complex interface have been observed, the extent of restriction being more significant around the directly bound residues of the protein and the DNA. In particular, the calculation revealed the existence of significantly constrained motionally restricted water layer that can form either bridges around the directly bound residues of the protein and DNA or are engaged in forming water-mediated contacts between a fraction of the unbound residues. More importantly, it is observed that the restricted water motion around the complex is correlated with the hydrogen bond relaxation time scale at the interface. It is further demonstrated that the kinetics of water-water hydrogen bonds involving the bridged water are influenced more due to complex formation.

Structural explanation for the tunable substrate specificity of an E. coli nucleoside hydrolase: insights from molecular dynamics simulations

Parasitic protozoa rely on nucleoside hydrolases that play key roles in the purine salvage pathway by catalyzing the hydrolytic cleavage of the N-glycosidic bond that connects nucleobases to ribose sugars. Cytidine-uridine nucleoside hydrolase (CU-NH) is generally specific toward pyrimidine nucleosides; however, previous work has shown that replacing two active site residues with Tyr, specifically the Thr223Tyr and Gln227Tyr mutations, allows CU-NH to process inosine. The current study uses molecular dynamics (MD) simulations to gain atomic-level insight into the activity of wild-type and mutant E. coli CU-NH toward inosine. By examining systems that differ in the identity and protonation states of active site catalytic residues, key enzyme-substrate interactions that dictate the substrate specificity of CU-NH are identified. Regardless of the wild-type or mutant CU-NH considered, our calculations suggest that inosine binding is facilitated by interactions of the ribose moiety with active site residues and Ca²⁺, and π-interactions between two His residues (His82 and His239) and the nucleobase. However, the lack of observed activity toward inosine for wild-type CU-NH is explained by no residue being correctly aligned to stabilize the departing nucleobase. In contrast, a hydrogen-bonding network between hypoxanthine and a newly identified general acid (Asp15) is present when the two Tyr mutations are engineered into the active site. Investigation of the single CU-NH mutants reveals that this hydrogen-bonding network is only maintained when both Tyr mutations are present due to a π-interaction between the residues. These results rationalize previous experiments that show the single Tyr mutants are unable to efficiently hydrolyze inosine and explain how the Tyr residues work synergistically in the double mutant to stabilize the nucleobase leaving group during hydrolysis. Overall, our simulations provide a structural explanation for the substrate specificity of nucleoside hydrolases, which may be used to rationally develop new treatments for kinetoplastid diseases.

QM/MM Study of the Reaction Catalyzed by Alkyladenine DNA Glycosylase: Examination of the Substrate Specificity of a DNA Repair Enzyme

Human alkyladenine DNA glycosylase (AAG) functions as part of the base excision repair pathway to excise structurally diverse oxidized and alkylated DNA purines. Specifically, AAG uses a water molecule activated by a general base and a nonspecific active site lined with aromatic residues to cleave the N-glycosidic bond. Despite broad substrate specificity, AAG does not target the natural purines (adenine (A) and guanine (G)). Using the ONIOM(QM:MM) methodology, we provide fundamental atomic level details of AAG bound to DNA-containing a neutral substrate (hypoxanthine (Hx)), a nonsubstrate (G), or a cationic substrate (7-methylguanine (7MeG)) and probe changes in the reaction pathway that occur when AAG targets different nucleotides. We reveal that subtle differences in protein-DNA contacts upon binding different substrates within the flexible AAG active site can significantly affect the deglycosylation reaction. Notably, we predict that AAG excises Hx in a concerted mechanism that is facilitated through correct alignment of the (E125) general base due to hydrogen bonding with a neighboring aromatic amino acid (Y127). Hx departure is further stabilized by π-π interactions with aromatic amino acids and hydrogen bonds with active site water. Despite possessing a similar structure to Hx, G is not excised since the additional exocyclic amino group leads to misalignment of the general base due to disruption of the key E125-Y127 hydrogen bond, the catalytically unfavorable placement of water within the active site, and weakened π-contacts between aromatic amino acids and the nucleobase. In contrast, cationic 7MeG does not occupy the same position within the AAG active site as G due to steric clashes with the additional N7 methyl group, which results in the correct alignment of the general base and permits nucleobase excision as observed for neutral Hx. Overall, our structural data rationalizes the observed substrate specificity of AAG and contributes to our fundamental understanding of enzymes with flexible active sites and broad substrate specificities.

The formation of catalytically competent enzyme-substrate complex is not a bottleneck in lesion excision by human alkyladenine DNA glycosylase

Human alkyladenine DNA glycosylase (AAG) protects DNA from alkylated and deaminated purine lesions. AAG flips out the damaged nucleotide from the double helix of DNA and catalyzes the hydrolysis of the N-glycosidic bond to release the damaged base. To understand better, how the step of nucleotide eversion influences the overall catalytic process, we performed a pre-steady-state kinetic analysis of AAG interaction with specific DNA-substrates, 13-base pair duplexes containing in the 7th position 1-N6-ethenoadenine (εA), hypoxanthine (Hx), and the stable product analogue tetrahydrofuran (F). The combination of the fluorescence of tryptophan, 2-aminopurine, and 1-N6-ethenoadenine was used to record conformational changes of the enzyme and DNA during the processes of DNA lesion recognition, damaged base eversion, excision of the N-glycosidic bond, and product release. The thermal stability of the duplexes characterized by the temperature of melting, T_m, and the rates of spontaneous opening of individual nucleotide base pairs were determined by NMR spectroscopy. The data show that the relative thermal stability of duplexes containing a particular base pair in position 7, (T_m(F/T) < T_m(εA/T) < T_m(Hx/T) < T_m(A/T)) correlates with the rate of reversible spontaneous opening of the base pair. However, in contrast to that, the catalytic lesion excision rate is two orders of magnitude higher for Hx-containing substrates than for substrates containing εA, proving that catalytic activity is not correlated with the stability of the damaged base pair. Our study reveals that the formation of the catalytically competent enzyme-substrate complex is not the bottleneck controlling the catalytic activity of AAG.