Identification of small molecule synthetic inhibitors of DNA polymerase beta by NMR chemical shift mapping

DNA polymerase beta (beta-pol) plays a central role in repair of damaged DNA bases by base excision repair (BER) pathways. A predominant phenotype of beta-pol null mouse fibroblasts is hypersensitivity to the DNA-methylating agent methyl methanesulfonate. Residues in the 8-kDa domain of beta-pol that seem to interact with a known natural product beta-pol inhibitor, koetjapic acid, were identified by NMR chemical shift mapping. The data implicate the binding pocket as the hydrophobic cleft between helix-2 and helix-4, which provides the DNA binding and deoxyribose phosphate lyase activities of the enzyme. Nine structurally related synthetic compounds, containing aromatic or other hydrophobic groups in combination with two carboxylate groups, were then tested. They were found to bind to the same or a very similar region on the surface of the enzyme. The ability of these compounds to potentiate methyl methanesulfonate cytotoxicity, an indicator of cellular BER capacity, in wild-type and beta-pol null mouse fibroblasts, was next ascertained. The most active and beta-pol-specific of these agents, pamoic acid, was further characterized and found to be an inhibitor of the deoxyribose phosphate lyase and DNA polymerase activities of purified beta-pol on a BER substrate. Our results illustrate that NMR-based mapping techniques can be used in the design of small molecule enzyme inhibitors including those with potential use in a clinical setting.     

Determination of the structure of the DNA binding domain of gamma delta resolvase in solution.

The DNA binding domain (DBD) of gamma delta resolvase (residues 141-183) is responsible for the interaction of this site-specific DNA recombinase with consensus site DNA within the gamma delta transposable element in Escherichia coli. Based on chemical-shift comparisons, the proteolytically isolated DBD displays side-chain interactions within a hydrophobic core that are highly similar to those of this domain when part of the intact enzyme (Liu T, Liu DJ, DeRose EF, Mullen GP, 1993, J Biol Chem 268:16309-16315). The structure of the DBD in solution has been determined using restraints obtained from 2-dimensional proton NMR data and is represented by 17 conformers. Experimental restraints included 458 distances based on analysis of nuclear Overhauser effect connectivities, 17 phi and chi 1 torsion angles based on analysis of couplings, and 17 backbone hydrogen bonds determined from NH exchange data. With respect to the computed average structure, these conformers display an RMS deviation of 0.67 A for the heavy backbone atoms and 1.49 A for all heavy atoms within residues 149-180. The DBD consists of 3 alpha-helices comprising residues D149-Q157, S162-T167, and R172-N183. Helix-2 and helix-3 form a backbone fold, which is similar to the canonical helix-turn-helix motif. The conformation of the NH2-terminal residues, G141-R148, appears flexible in solution. A hydrophobic core is formed by side chains donated by essentially all hydrophobic residues within the helices and turns. Helix-1 and helix-3 cross with a right-handed folding topology. The structure is consistent with a mechanism of DNA binding in which contacts are made by the hydrophilic face of helix-3 in the major groove and the amino-terminal arm in the minor groove. This structure represents an important step toward analysis of the mechanism of DNA interaction by gamma delta resolvase and provides initial structure-function comparisons among the divergent DBDs of related resolvases and invertases.     

Sequence-specific oxidative cleavage of DNA by a designed metalloprotein, Ni(II).GGH(Hin139-190).

A 55-residue protein containing the DNA binding domain of Hin recombinase, residues 139-190, with the tripeptide Gly-Gly-His (GGH) at the NH2 terminus was synthesized by stepwise solid-phase methods. GGH(Hin139-190) binds sequence specifically to DNA at four 13 base pair sites (termed hixL and secondary) and, in the presence of Ni(OAc)2 and monoperoxyphthalic acid, reacts predominantly at a single deoxyribose position on one strand of each binding site [Mack, D.P., & Dervan, P.B. (1990) J. Am. Chem. Soc. 112, 4604]. We find that, upon treatment with n-butylamine, the DNA termini at the cleavage site are 3'- and 5'-phosphate, consistent with oxidative degradation of the deoxyribose backbone. The nickel-mediated oxidation can be activated with peracid, iodosylbenzene, or hydrogen peroxide. The sequence specificity of the reaction is not dependent on oxidant, but the rates of cleavage differ, decreasing in the order peracid greater than iodosylbenzene greater than hydrogen peroxide. Optimal cleavage conditions for a 1 microM concentration of protein are 50 microM peracid, pH 8.0, and 1 equiv of Ni(OAc)2. The preferential cleavage at a single base pair position on one strand of the minor groove indicates a nondiffusible oxidizing species. A change of absolute configuration in the GGH metal binding domain from L-His to D-His [Ni(II).GG-(-D-)H(Hin139-190)] affords cleavage at similar base pair locations but opposite with regard to strand specificity.     

Retinoic acid receptor: a simulation analysis of retinoic acid binding and the resulting conformational changes.

The binding/escape mechanism of all- trans retinoic acid with respect to the ligand-binding domain of the nuclear receptor RARgamma has been studied by molecular dynamic simulations. The entry/exit channel was shown to be on the side of the activation helix by the use of multiple copy dynamics. Three independent minimum energy paths from the liganded structure to a model for the unliganded structure were calculated with the conjugate peak refinement method. Ligand escape takes place in the early steps of the transition during rearrangement of the binding pocket; the latter involves inward motion of the beta-sheet and outward motions of the Omega-loop and helix H6. The correlated rearrangements involved in the escape phase are similar and occur in the same order for the different paths. After the escape phase, the conformational changes affect primarily the C-terminal helices H11-H12 and the Omega-loop. The three paths are significantly different for this reorganization phase and reveal a multiplicity of possibilities, in agreement with the idea that the apo state is structurally less constrained. The present calculations extend the crystallographic results, confirming the "mouse trap" mechanism and stressing the importance of the helix H3 conformation and of the contacts between the Omega-loop and helices H11 and H6. They are in good agreement with known mutants and point to other functionally important residues, especially in helices H3 and H11, suggesting mutations that may affect the ligand-binding function and the associated conformational changes.     

Structural insights into abasic site for Fpg specific binding and catalysis: comparative high-resolution crystallographic studies of Fpg bound to various models of abasic site analogues-containing DNA

Fpg is a DNA glycosylase that recognizes and excises the mutagenic 8-oxoguanine (8-oxoG) and the potentially lethal formamidopyrimidic residues (Fapy). Fpg is also associated with an AP lyase activity which successively cleaves the abasic (AP) site at the 3′ and 5′ sides by βδ-elimination. Here, we present the high-resolution crystal structures of the wild-type and the P1G defective mutant of Fpg from Lactococcus lactis bound to 14mer DNA duplexes containing either a tetrahydrofuran (THF) or 1,3-propanediol (Pr) AP site analogues. Structures show that THF is less extrahelical than Pr and its backbone C5′–C4′–C3′ diverges significantly from those of Pr, rAP, 8-oxodG and FapydG. Clearly, the heterocyclic oxygen of THF is pushed back by the carboxylate of the strictly conserved E2 residue. We can propose that the ring-opened form of the damaged deoxyribose is the structure active form of the sugar for Fpg catalysis process. Both structural and functional data suggest that the first step of catalysis mediated by Fpg involves the expulsion of the O4′ leaving group facilitated by general acid catalysis (involving E2), rather than the immediate cleavage of the N-glycosic bond of the damaged nucleoside.     

Footprinting of Chlorella virus DNA ligase bound at a nick in duplex DNA.

The 298-amino acid ATP-dependent DNA ligase of Chlorella virus PBCV-1 is the smallest eukaryotic DNA ligase known. The enzyme has intrinsic specificity for binding to nicked duplex DNA. To delineate the ligase-DNA interface, we have footprinted the enzyme binding site on DNA and the DNA binding site on ligase. The size of the exonuclease III footprint of ligase bound a single nick in duplex DNA is 19-21 nucleotides. The footprint is asymmetric, extending 8-9 nucleotides on the 3'-OH side of the nick and 11-12 nucleotides on the 5'-phosphate side. The 5'-phosphate moiety is essential for the binding of Chlorella virus ligase to nicked DNA. Here we show that the 3'-OH moiety is not required for nick recognition. The Chlorella virus ligase binds to a nicked ligand containing 2',3'-dideoxy and 5'-phosphate termini, but cannot catalyze adenylation of the 5'-end. Hence, the 3'-OH is important for step 2 chemistry even though it is not itself chemically transformed during DNA-adenylate formation. A 2'-OH cannot substitute for the essential 3'-OH in adenylation at a nick or even in strand closure at a preadenylated nick. The protein side of the ligase-DNA interface was probed by limited proteolysis of ligase with trypsin and chymotrypsin in the presence and absence of nicked DNA. Protease accessible sites are clustered within a short segment from amino acids 210-225 located distal to conserved motif V. The ligase is protected from proteolysis by nicked DNA. Protease cleavage of the native enzyme prior to DNA addition results in loss of DNA binding. These results suggest a bipartite domain structure in which the interdomain segment either comprises part of the DNA binding site or undergoes a conformational change upon DNA binding. The domain structure of Chlorella virus ligase inferred from the solution experiments is consistent with the structure of T7 DNA ligase determined by x-ray crystallography.     

Sulfoquinovosylmonoacylglycerol inhibitory mode analysis of rat DNA polymerase beta

We have previously reported that sulfoquinovosylmonoacylglycerol (SQMG) is a potent inhibitor of mammalian DNA polymerases. DNA polymerase beta (pol beta) is one of the most important enzymes protecting the cell against DNA damage by base excision repair. In this study, we characterized the inhibitory action of SQMG against rat pol beta. SQMG competed with both the substrate and the template-primer for binding to pol beta. A gel mobility shift assay and a polymerase activity assay showed that SQMG competed with DNA for a binding site on the N-terminal 8-kDa domain of pol beta, subsequently inhibiting its catalytic activity. Fragments of SQMG such as sulfoquinovosylglycerol (SQG) and fatty acid (myristoleic acid, MA) weakly inhibited pol beta activity and the inhibitory effect of a mixture of SQG and MA was stronger than that of SQG or MA. To characterize this inhibition more precisely, we attempted to identify the interaction interface between SQMG and the 8-kDa domain by NMR chemical shift mapping. Firstly, we determined the binding site on a fragment of SQMG, the SQG moiety. We observed chemical shift changes primarily at two sites, the residues comprising the C-terminus of helix-1 and the N-terminus of helix-2, and residues in helix-4. Finally, based on our present results and our previously reported study of the interaction interface of fatty acids, we constructed two three-dimensional models of a complex between the 8-kDa domain and SQMG and evaluated them by the mutational analysis. The models show a SQMG interaction interface that is consistent with the data.     

The Leu22Pro tumor-associated variant of DNA polymerase beta is dRP lyase deficient

Approximately 30% of human tumors characterized to date express DNA polymerase beta (pol β) variant proteins. Two of the polymerase beta cancer-associated variants are sequence-specific mutators, and one of them binds to DNA but has no polymerase activity. The Leu22Pro (L22P) DNA polymerase beta variant was identified in a gastric carcinoma. Leu22 resides within the 8 kDa amino terminal domain of DNA polymerase beta, which exhibits dRP lyase activity. This domain catalyzes the removal of deoxyribose phosphate during short patch base excision repair. We show that this cancer-associated variant has very little dRP lyase activity but retains its polymerase activity. Although residue 22 has no direct contact with the DNA, we report here that the L22P variant has reduced DNA-binding affinity. The L22P variant protein is deficient in base excision repair. Molecular dynamics calculations suggest that alteration of Leu22 to Pro changes the local packing, the loop connecting helices 1 and 2 and the overall juxtaposition of the helices within the N-terminal domain. This in turn affects the shape of the binding pocket that is required for efficient dRP lyase catalysis.     

Mechanism of DNA strand nicking at apurinic/apyrimidinic sites by Escherichia coli [formamidopyrimidine]DNA glycosylase.

Escherichia coli [formamidopyrimidine]DNA glycosylase catalyses the nicking of both the phosphodiester bonds 3' and 5' of apurinic or apyrimidinic sites in DNA so that the base-free deoxyribose is replaced by a gap limited by 3'-phosphate and 5'-phosphate ends. The two nickings are not the results of hydrolytic processes; the [formamidopyrimidine]DNA glycosylase rather catalyses a beta-elimination reaction that is immediately followed by a delta-elimination. The enzyme is without action on a 3'-terminal base-free deoxyribose or on a 3'-terminal base-free unsaturated sugar produced by a beta-elimination reaction nicking the DNA strand 3' to an apurinic or apyrimidinic site.     

Product analyses in DNA strand scission by antitumor antibiotic elsamicin A.

Elsamicin A is an antitumor antibiotic with fascinating chemical structure and a good candidate for pharmaceutical development. Molecular mechanism of DNA backbone cleavage mediated by Fe(II)-elsamicin A has been examined. Product analysis using DNA sequencing gels and HPLC reveals the production of damaged DNA fragments bearing 3'-/5'-phosphate and 3'-phosphoglycolate termini associated with formation of free base. In addition, hydrazine-trapping experiments indicate that C-4' hydroxylated abasic sites are formed concomitant with DNA degradation by Fe(II)-elsamicin A. The results lead to the conclusion that the hydroxyl radical formed in Fe(II)-elsamicin A plus dithiothreitol system oxidizes the deoxyribose moiety via hydrogen abstraction predominantly at the C-4' carbon of the deoxyribose backbone and ultimately produces strand breakage of DNA.     

A plant phytotoxin, solanapyrone A, is an inhibitor of DNA polymerase beta and lambda.

Solanapyrone A, a phytotoxin and enzyme inhibitor isolated from a fungus (SUT 01B1-2) selectively inhibits the activities of mammalian DNA polymerase beta and lambda (pol beta and lambda) in vitro. The IC50 values of the compound were 30 microm for pol beta and 37 microm for pol lambda. Because pol beta and lambda are in a family and their three-dimensional structures are thought to be highly similar to each other, we used pol beta to analyze the biochemical relationship with solanapyrone A. On pol beta, solanapyrone A antagonistically competed with both the DNA template and the nucleotide substrate. BIAcore analysis demonstrated that solanapyrone A bound selectively to the N-terminal 8-kDa domain of pol beta. This domain is known to bind single-stranded DNA, provide 5'-phosphate recognition of gapped DNA, and cleave the sugar-phosphate bond 3' to an intact apurinic/apyrimidinic (AP) site (i.e. AP lyase activity) including 5'-deoxyribose phosphate lyase activity. Solanapyrone A inhibited the single-stranded DNA-binding activity but did not influence the activities of the 5'-phosphate recognition in gapped DNA structures and the AP lyase. Based on these results, the inhibitory mechanism of solanapyrone A is discussed.     

Rapid segmental and subdomain motions of DNA polymerase beta

DNA polymerase (pol) beta is a two-domain DNA repair enzyme that undergoes structural transitions upon binding substrates. Crystallographic structures indicate that these transitions include movement of the amino-terminal 8-kDa lyase domain relative to the 31-kDa polymerase domain. Additionally, a polymerase subdomain moves toward the nucleotide-binding pocket after nucleotide binding, resulting in critical contacts between alpha-helix N and the nascent base pair. Kinetic and structural characterization of pol beta has suggested that these conformational changes participate in stabilizing the ternary enzyme-substrate complex facilitating chemistry. To probe the microenvironment and dynamics of both the lyase domain and alpha-helix N in the polymerase domain, the single native tryptophan (Trp-325) of wild-type enzyme was replaced with alanine, and tryptophan was strategically substituted for residues in the lyase domain (F25W/W325A) or near the end of alpha-helix N (L287W/W325A). Influences of substrate on the fluorescence anisotropy decay of these single tryptophan forms of pol beta were determined. The results revealed that the segmental motion of alpha-helix N was rapid ( approximately 1 ns) and far more rapid than the step that limits chemistry. Binding of Mg(2+) and/or gapped DNA did not cause a noticeable change in the rotational correlation time or angular amplitude of tryptophan in alpha-helix N. More important, binding of a correct nucleotide significantly limited the angular range of the nanosecond motion within alpha-helix N. In contrast, the segmental motion of the 8-kDa domain was "frozen" upon DNA binding alone, and this restriction did not increase further upon nucleotide binding. The dynamics of alpha-helix N are discussed from the perspective of the "open" to "closed" conformational change of pol beta deduced from crystallography, and the results are more generally discussed in the context of reaction cycle-regulated flexibility for proteins acting as molecular motors.     

The N-terminal domain of the human Rad51 protein binds DNA: structure and a DNA binding surface as revealed by NMR.

Human Rad51 protein (HsRad51) is a homolog of Escherichia coli RecA protein, and functions in DNA repair and recombination. In higher eukaryotes, Rad51 protein is essential for cell viability. The N-terminal region of HsRad51 is highly conserved among eukaryotic Rad51 proteins but is absent from RecA, suggesting a Rad51-specific function for this region. Here, we have determined the structure of the N-terminal part of HsRad51 by NMR spectroscopy. The N-terminal region forms a compact domain consisting of five short helices, which shares structural similarity with a domain of endonuclease III, a DNA repair enzyme of E. coli. NMR experiments did not support the involvement of the N-terminal domain in HsRad51-HsBrca2 interaction or the self-association of HsRad51 as proposed by previous studies. However, NMR tiration experiments demonstrated a physical interaction of the domain with DNA, and allowed mapping of the DNA binding surface. Mutation analysis showed that the DNA binding surface is essential for double-stranded and single-stranded DNA binding of HsRad51. Our results suggest the presence of a DNA binding site on the outside surface of the HsRad51 filament and provide a possible explanation for the regulation of DNA binding by phosphorylation within the N-terminal domain.     

Structural insight into the DNA polymerase beta deoxyribose phosphate lyase mechanism

A large number of biochemical and genetic studies have demonstrated the involvement of DNA polymerase beta (Pol beta) in mammalian base excision repair (BER). Pol beta participates in BER sub-pathways by contributing gap filling DNA synthesis and lyase removal of the 5'-deoxyribose phosphate (dRP) group from the cleaved abasic site. To better understand the mechanism of the dRP lyase reaction at an atomic level, we determined a crystal structure of Pol beta complexed with 5'-phosphorylated abasic sugar analogs in nicked DNA. This DNA ligand represents a potential BER intermediate. The crystal structure reveals that the dRP group is bound in a non-catalytic binding site. The catalytic nucleophile in the dRP lyase reaction, Lys72, and all other potential secondary nucleophiles, are too far away to participate in nucleophilic attack on the C1' of the sugar. An approximate model of the dRP group in the expected catalytic binding site suggests that a rotation of 120 degrees about the dRP 3'-phosphate is required to position the epsilon-amino Lys72 close to the dRP C1'. This model also suggests that several other side chains are in position to facilitate the beta-elimination reaction. From results of mutational analysis of key residues in the dRP lyase active site, it appears that the substrate dRP can be stabilized in the observed non-catalytic binding conformation, hindering dRP lyase activity.     

Protein dynamics and monomer-monomer interactions in AntR activation by electron paramagnetic resonance and double electron-electron resonance

The Anthracis repressor (AntR) is a Mn(II)-activated DNA binding protein that is involved in the regulation of Mn(II) homeostasis in Bacillus anthracis. AntR is structurally and functionally homologous to Mn(II)-activated repressor from Bacillus subtillis (MntR). Our studies on AntR focus on metal-regulated activation of the protein. Line shape analysis of continuous wave electron paramagnetic resonance (EPR) spectra showed that metal binding resulted in a general reduction of backbone dynamics and that there were no further changes in backbone motion upon DNA binding. Double electron-electron resonance (DEER) pulsed EPR spectroscopy was used to measure distances between nitroxide spin labels strategically placed in dimeric AntR. The DEER data were analyzed assuming Gaussian distributions for discrete populations of spins. A structural model for AntR was built from homology to MntR, and the experimentally measured distances were simulated to distinguish between spin label and backbone motions. Together with the computational analysis, the DEER results for apo-AntR indicated relatively narrow conformational distributions for backbone residues at the dimer interface and near the metal binding site. No significant changes were observed on these sites in the presence of metal or DNA. On the other hand, the distribution of the conformers and the distances between the putative DNA binding helices decreased upon metal binding. These results suggest that the DNA binding region of AntR shows large amplitude backbone motions in the absence of metal, which may preclude sequence-specific binding to promoter sites. Metal binding narrows the range of conformations accessible in this region and shortens the mean distance between the DNA binding helices, probably resulting in alignment that optimizes promoter recognition and binding.     

Backbone dynamics of the c-Myb DNA-binding domain complexed with a specific DNA

The DNA-binding domain of c-Myb consists of three imperfect tandem repeats, R1, R2, and R3. Each repeat contains three helices. The minimal DNA-binding domain is an R2R3 fragment. Here, we have examined the backbone dynamics of R2R3 in its DNA-bound form by NMR. Upon binding to DNA, the N- and C-termini, and the linker between R2 and R3 become less flexible. In the free form the third helix of R2 exhibits slow conformational exchange fluctuations owing to a cavity in the hydrophobic core of R2. Upon binding to DNA, the conformational exchange contributions in R2 are reduced but remain significant in NMR relaxation measurements. Upon binding to DNA, the third helix of R3 comes to exhibit significant chemical exchange contributions. These findings suggest that the orientations of the third helices of both R2 and R3 as to DNA are being chemically exchanged. In the DNA-bound form both R2 and R3 exhibit similar dynamical characters, except for amino acids Trp 95, Thr 96, and Val 103 of R2, which are located around the cavity of the unbound form. Upon binding to DNA, since Trp 95 moves into the cavity to fill it up, the local conformational exchange contributions seem to be still observable around the filled cavity.     

T4 DNA ligase can seal a nick in double-stranded DNA limited by a 5''-phosphorylated base-free deoxyribose residue.

The 5' AP endodeoxyribonucleases hydrolyze the phosphodiester bond 5' to AP (apurinic or apyrimidinic) sites in double-stranded DNA leaving 3'-OH and 5'-phosphate ends. These nicks are sealed by T4 DNA ligase although the 5'-phosphate end belongs to a base-free deoxyribose.