The effect of spermine binding on the reactivity of DNA towards carcinogenic alkylating agents

The effect of spermine binding on the electrostatic potential of DNA is evaluated. The calculations are performed for the essential reactive sites, atoms N7 and O6 of guanine, N3 and N7 of adenine, of the nucleic acid and for its surface envelope. An important weakening of the potential is found affecting all the important reactive sites in both grooves and spreading moreover along the polynucleotide chain far away from the site of binding of the ligand. These results are discussed in connection with the experimentally observed inhibitory effect of spermine binding on DNA methylation by carcinogenic agents.     

Substitution for Asn460 cripples β-galactosidase (Escherichia coli) by increasing substrate affinity and decreasing transition state stability

Substrate initially binds to β-galactosidase (Escherichia coli) at a 'shallow' site. It then moves ～3? to a 'deep' site and the transition state forms. Asn460 interacts in both sites, forming a water bridge interaction with the O3 hydroxyl of the galactosyl moiety in the shallow site and a direct H-bond with the O2 hydroxyl of the transition state in the deep site. Structural and kinetic studies were done with β-galactosidases with substitutions for Asn460. The substituted enzymes have enhanced substrate affinity in the shallow site indicating lower E·substrate complex energy levels. They have poor transition state stabilization in the deep site that is manifested by increased energy levels of the E·transition state complexes. These changes in stability result in increased activation energies and lower k(cat) values. Substrate affinity to N460D-β-galactosidase was enhanced through greater binding enthalpy (stronger H-bonds through the bridging water) while better affinity to N460T-β-galactosidase occurred because of greater binding entropy. The transition states are less stable with N460S- and N460T-β-galactosidase because of the weakening or loss of the important bond to the O2 hydroxyl of the transition state. For N460D-β-galactosidase, the transition state is less stable due to an increased entropy penalty.     

Stable isotope labeling-mass spectrometry analysis of methyl- and pyridyloxobutyl-guanine adducts of 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone in p53-derived DNA sequences

The p53 tumor suppressor gene is a primary target in smoking-induced lung cancer. Interestingly, p53 mutations observed in lung tumors of smokers are concentrated at guanine bases within endogenously methylated (Me)CG dinucleotides, e.g., codons 157, 158, 245, 248, and 273 ((Me)C = 5-methylcytosine). One possible mechanism for the increased mutagenesis at these sites involves targeted binding of metabolically activated tobacco carcinogens to (Me)CG sequences. In the present work, a stable isotope labeling HPLC-ESI(+)-MS/MS approach was employed to analyze the formation of guanine lesions induced by the tobacco-specific lung carcinogen 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone (NNK) within DNA duplexes representing p53 mutational "hot spots" and surrounding sequences. Synthetic DNA duplexes containing p53 codons 153-159, 243-250, and 269-275 were prepared, where (Me)C was incorporated at all physiologically methylated CG sites. In each duplex, one of the guanine bases was replaced with [1,7,NH(2)-(15)N(3)-2-(13)C]-guanine, which served as an isotope "tag" to enable specific quantification of guanine lesions originating from that position. After incubation with NNK diazohydroxides, HPLC-ESI(+)-MS/MS analysis was used to determine the yields of NNK adducts at the isotopically labeled guanine and at unlabeled guanine bases elsewhere in the sequence. We found that N7-methyl-2'-deoxyguanosine and N7-[4-oxo-4-(3-pyridyl)but-1-yl]guanine lesions were overproduced at the 3'-guanine bases within polypurine runs, while the formation of O(6)-methyl-2'-deoxyguanosine and O(6)-[4-oxo-4-(3-pyridyl)but-1-yl]-2'-deoxyguanosine adducts was specifically preferred at the 3'-guanine base of 5'-GG and 5'-GGG sequences. In contrast, the presence of 5'-neighboring (Me)C inhibited O(6)-guanine adduct formation. These results indicate that the N7- and O(6)-guanine adducts of NNK are not overproduced at the endogenously methylated CG dinucleotides within the p53 tumor suppressor gene, suggesting that factors other than NNK adduct formation are responsible for mutagenesis at these sites.     

Informational parameters of an exact DNA base sequence

A series of hierarchical chemical reactivity calculations have been performed to elucidate the alkylation properties of a methyldiazonium ion toward DNA base sites. Both MINDO/3 and CNDO/2 approximate methods have been employed. For the isolated bases the O6 of guanine is predicted to be the most reactive site. This prediction may also be relevant to single-stranded DNA chains containing guanine. For base-pairs, the N7 and O6 sites on guanine are about equally favored for alkylation. The previous study of aziridinium ion alkylation gave about the same results with N7 guanine modestly favored as the preferred site of alkylation for base-pairs. In composite we conclude that N7 guanine and/or O6 guanine will be the preferred sites for alkylation by a methyldiazonium ion but cannot distinguish between these two in terms of chemical specificity.     

A theoretical evaluation of the effect of netropsin binding on the reactivity of DNA towards alkylating agents.

The effect of netropsin binding on the electrostatic potential of DNA reactive sites is presented. Calculations are performed for atoms N7 and O6 of guanine, N3 and N7 of adenine of model, 25 base pair long, DNA-netropsin complexes. An important weakening of the potential is found spreading along all the oligonucleotide chain studied. The results are discussed in connection with the inhibitory effect of a related ligand, distamycin A, on DNA methylation.     

The Effect of Spermine Binding on the Reactivity of DNA Towards Carcinogenic Alkylating Agents

Abstract

The effect of spermine binding on the electrostatic potential of DNA is evaluated. The calculations are performed for the essential reactive sites, atoms N7 and 06 of guanine, N3 and N7 of adenine, of the nucleic acid and for its surface envelope. An important weakening of the potential is found affecting all the important reactive sites in both grooves and spreading moreover along the polynucleotide chain far away from the site of binding of the ligand. These results are discussed in connection with the experimentally observed inhibitory effect of spermine binding on DNA methylation by carcinogenic agents.     

Binding of transition metal complexes to guanine and guanine–cytosine: hydrogen bonding and covalent effects

Density functional calculations and Atoms in Molecules analysis are used to investigate the role of covalent and hydrogen bondings in determining the binding of transition metal complexes to guanine, and the subsequent effect on pairing with cytosine. Hydrogen bonding is ubiquitous, and typically contributes ca. 10% to overall binding, a value that varies with the coordination site on guanine, as well as metal and ligands. Early transition metals show a clear preference for the O6 position, while later ones prefer N7, the crossover point coming at the vanadium group. Metallation at N7 causes a redistribution of hydrogen bonding strength between guanine and cytosine, but does not greatly affect the overall pairing energy. In contrast, metallation at O6 strongly reduces the pairing energy, as may be expected given the role of O6 in pairing guanine with cytosine. This effect can be quantified using electron density properties, and seems to be due to both electrostatic repulsion from the positive metal centre and a redistribution of electron density within guanine itself. Qualitative agreement with experimental mass spectroscopic results is obtained.     

Insight into the phosphodiesterase mechanism from combined QM/MM free energy simulations

Molecular dynamics simulations employing a combined quantum mechanical and molecular mechanical potential have been carried out to elucidate the reaction mechanism of the hydrolysis of a cyclic nucleotide cAMP substrate by phosphodiesterase 4B (PDE4B). PDE4B is a member of the PDE superfamily of enzymes that play crucial roles in cellular signal transduction. We have determined a two-dimensional potential of mean force (PMF) for the coupled phosphoryl bond cleavage and proton transfer through a general acid catalysis mechanism in PDE4B. The results indicate that the ring-opening process takes place through an S(N)2 reaction mechanism, followed by a proton transfer to stabilize the leaving group. The computed free energy of activation for the PDE4B-catalyzed cAMP hydrolysis is about 13 kcal·mol(-1) and an overall reaction free energy is about -17 kcal·mol(-1), both in accord with experimental results. In comparison with the uncatalyzed reaction in water, the enzyme PDE4B provides a strong stabilization of the transition state, lowering the free energy barrier by 14 kcal·mol(-1). We found that the proton transfer from the general acid residue His234 to the O3' oxyanion of the ribosyl leaving group lags behind the nucleophilic attack, resulting in a shallow minimum on the free energy surface. A key contributing factor to transition state stabilization is the elongation of the distance between the divalent metal ions Zn(2+) and Mg(2+) in the active site as the reaction proceeds from the Michaelis complex to the transition state.     

A selective recognition mode of a nucleic acid base by an aromatic amino acid: L-phenylalanine-7-methylguanosine 5''-monophosphate stacking interaction.

The conformation of 7-methylguanosine 5'-monophosphate (m7GMP) and its interaction with L-phenylalanine (Phe) have been investigated by X-ray crystallographic, 1H-nuclear magnetic resonance, and energy calculation methods. The N(7) methylation of the guanine base shifts m7GMP toward an anti--gauche, gauche conformation about the glycosyl and exocyclic C(4')-C(5') bonds, respectively. The prominent stacking observed between the benzene ring of Phe and guanine base of m7GMP is primarily due to the N(7) guarternization of the guanine base. The formation of a hydrogen bonding pair between the anionic carboxyl group and the guanine base further stabilizes this stacking interaction. The present results imply the importance of aromatic amino acids as a hallmark for the selective recognition of a nucleic acid base.     

Chromomycin dimer-DNA oligomer complexes. Sequence selectivity and divalent cation specificity

This paper reports on a solution NMR characterization of the sequence selectivity and metal ion specificity in chromomycin-DNA oligomer complexes in the presence of divalent cations. The sequence selectivity studies have focused on chromomycin complexes with the self-complementary d(A1-A2-G3-G4-C5-C6-T7-T8) duplex containing a pair of adjacent (G3-G4).(C5-C6) steps and the self-complementary d(A1-G2-G3-A4-T5-C6-C7-T8) duplex containing a pair of separated (G2-G3).(C6-C7) steps in aqueous solution. The antitumor agent (chromomycin) and nucleic acid protons have been assigned following analysis of distance connectivities in NOESY spectra and coupling connectivities in DQF-COSY spectra for both complexes in H2O and D2O solution. The observed intermolecular NOEs establish that chromomycin binds as a Mg(II)-coordinated dimer [1 Mg(II) per complex] and contacts the minor-groove edge with retention of 2-fold symmetry centered about the (G3-G4-C5-C6).(G3-G4-C5-C6) segment of the d(A2G2C2T2) duplex. By contrast, complex formation is centered about the (G2-G3-A4-T5).(A4-T5-C6-C7) segment and results in removal of the two fold symmetry of the d(AG2ATC2T) duplex. Thus, the binding of one subunit of the chromomycin dimer at its preferred (G-G).(C-C) site assists in the binding of the second subunit to the less preferred adjacent (A-T).(A-T) site. These observations suggest a hierarchy of chromomycin binding sites, with a strong site detected at the (G-G) step due to the hydrogen-bonding potential of acceptor N3 and donor NH2 groups of guanosine that line the minor groove. The divalent cation specificity has been investigated by studies on the symmetric chromomycin-d(A2G2C2T2) complex in the presence of diamagnetic Mg(II), Zn(II), and Cd(II) cations and paramagnetic Ni(II) and Co(II) cations. A comparative NOESY study of the Mg(II) and Ni(II) symmetric complexes suggests that a single tightly bound divalent cation aligns the two chromomycins in the dimer through coordination to the C1 carbonyl and C9 enolate ions on the hydrophilic edge of each aglycon ring. Secondary divalent cation binding sites involve coordination to the major-groove N7 atoms on adjacent guanosines in G-G steps. This coordination is perturbed on lowering the pH below 6.0, presumably due to protonation of the N7 atoms. The midpoint of the thermal dissociation of the symmetric complex is dependent on the divalent cation with the stability for reversible transitions decreasing in the order Mg(II) greater than Zn(II) greater than Cd(II) complexes.(ABSTRACT TRUNCATED AT 400 WORDS)     

Simulation of DNA bases in water: Comparison of the Monte Carlo algorithm with molecular mechanics force fields

The interaction between the nucleic acid bases and solvent molecules has an important effect in various biochemical processes. We have calculated total energy and free energy of the solvation of DNA bases in water by Monte Carlo simulation. Adenine, guanine, cytosine, and thymine were first optimized in the gas phase and then placed in a cubic box of water. We have used the TIP3 model for water and OPLS for the nucleic acid bases. The canonical (T, V, N) ensemble at 25°C and Metropolis sampling technique have been used. Good agreement with other available computational data was obtained. Radial distribution functions of water around each site of adenine, guanine, cytosine, and thymine have been computed and the results have shown the ability of the sites for hydrogen bonding and other interactions. The computations have shown that guanine has the highest value of solvation free energy and N7 and N6 in adenine and guanine, N3 in cytosine, and N3 and O4 in thymine have the largest radial distribution function. Monte Carlo simulation has also been performed using the CHARMM program under the same conditions, and the results of two procedures are compared.     

O6-ethylguanine carcinogenic lesions in DNA: an NMR study of O6etG.T pairing in dodecanucleotide duplexes

High-resolution two-dimensional NMR studies are reported on the self-complementary d-(C1-G2-C3-O6etG4-A5-G6-C7-T8-T9-G10-C11-G12) duplex (designated O6etG.T 12-mer) containing two symmetrically related O6etG.T lesion sites located four base pairs in from either end of the duplex. Parallel studies were undertaken on a related sequence containing O6meG.T lesion sites (designated O6meG.T 12-mer) in order to evaluate the influence of the size of the alkyl substituent on the structure of the duplex and were undertaken on a related sequence containing G.T mismatch sites (designated G.T 12-mer duplex), which served as the control duplex. The exchangeable and nonexchangeable proton and the phosphorus nuclei have been assigned from an analysis of two-dimensional nuclear Overhauser enhancement (NOE) and correlated spectra of the O6etG.T 12-mer, O6meG.T 12-mer, and G.T 12-mer duplexes in H2O and D2O solutions. The distance connectivities observed in the NOESY spectra of the O6alkG.T 12-mer duplexes establish that the helix is right-handed and all of the bases adopt an anti conformation of the glycosidic torsion angle including the O6alkG4 and T9 bases at the lesion site. The imino proton of T9 at the O6alkG.T lesion sites resonates at 8.85 ppm in the O6etG.T 12-mer duplex and at 9.47 ppm in the O6meG.T 12-mer duplex. The large upfield shift of the T9 imino proton resonance at the O6alkG4.T9 lesion site relative to that of the same proton in the G4.T9 wobble pair (11.99 ppm) and the A4.T9 Watson-Crick pair (13.95 ppm) in related sequences establishes that the hydrogen bonding of the imino proton of T9 to O6alkG4 is either very weak or absent. The imino proton of T9 develops NOEs to the CH3 protons of the O6etG and O6meG alkyl groups across the base pair, as well as to the imino and H5 protons of the flanking C3.G10 base pair and the imino and CH3 protons of the flanking A5.T8 base pair in the O6alkG.T 12-mer duplexes. These observations establish that the O6alkG4 and T9 residues are stacked into the duplex and that the O6CH3 and O6CH2CH3 groups of O6alkG4 adopt a syn orientation with respect to the N1 of the alkylated guanine.(ABSTRACT TRUNCATED AT 400 WORDS)     

Simulation of DNA bases in water: Comparison of the Monte Carlo algorithm with molecular mechanics force fields

The interaction between the nucleic acid bases and solvent molecules has an important effect in various biochemical processes. We have calculated total energy and free energy of the solvation of DNA bases in water by Monte Carlo simulation. Adenine, guanine, cytosine, and thymine were first optimized in the gas phase and then placed in a cubic box of water. We have used the TIP3 model for water and OPLS for the nucleic acid bases. The canonical (T, V, N) ensemble at 25 degrees C and Metropolis sampling technique have been used. Good agreement with other available computational data was obtained. Radial distribution functions of water around each site of adenine, guanine, cytosine, and thymine have been computed and the results have shown the ability of the sites for hydrogen bonding and other interactions. The computations have shown that guanine has the highest value of solvation free energy and N7 and N6 in adenine and guanine, N3 in cytosine, and N3 and O4 in thymine have the largest radial distribution function. Monte Carlo simulation has also been performed using the CHARMM program under the same conditions, and the results of two procedures are compared.     

O6-ethylguanine carcinogenic lesions in DNA: an NMR study of O6etG.C pairing in dodecanucleotide duplexes

The pairing of O6etG with C located four base pairs in from either end of the self-complementary d(C1-G2-C3-O6etG4-A5-G6-C7-T8-C9-G10-C11-G12) duplex (designated O6etG.C 12-mer) has been investigated from an analysis of proton and phosphorus two-dimensional NMR experiments. The structural consequences of increasing the alkyl group size were elucidated from a comparative study of the pairing of O6meG4 with C9 in a related sequence (designated O6meG.C 12-mer). The NMR parameters for both O6alkG-containing dodecanucleotides are also compared with those of the control sequence containing G4.C9 base pairs (designated G.C 12-mer). The NOE cross-peaks detected in the two-dimensional NOESY spectra of the O6alkG.C 12-mer duplexes in H2O solution establish that the O6etG4/O6meG4 and C9 bases at the lesion site stack into the helix between the flanking C3.G10 and A5.T8 Watson-Crick base pairs. The amino protons of C9 at the O6alkG4-C9 lesion site resonate as an average resonance at 7.78 and 7.63 ppm in the O6etG.C 12-mer and O6meG.C 12-mer duplexes, respectively. The observed NOEs between the amino protons of C9 and the CH3 protons of O6alkG4 establish a syn orientation of the O6-alkyl group with respect to the N1 of alkylated guanine. A wobble alignment of the O6alkG4.C9 base pair stablized by two hydrogen bonds, one between the amino group of C9 and N1 of O6alkG and the other between the amino group of O6alkG and N3 of C9, is tentatively proposed on the basis of the NOEs between the amino protons of C9 at the lesion site and the imino protons of flanking Watson-Crick base pairs. The proton and phosphorus chemical shift differences between the O6etG.C 12-mer and O6meG.C 12-mer duplexes are small compared to the differences between these O6alkG-containing duplexes and the control G.C 12-mer duplex.     

Mechanism of CpG DNA methyltransferases M.SssI and Dnmt3a studied by DNA containing 2-aminopurine

Murine DNA methyltransferases Dnmt3a-CD and M.SssI from Spiroplasma methylate cytosines at CpG sites. The role of 6-oxo groups of guanines in DNA methylation by these enzymes has been studied using DNA substrates, which contained 2-aminopurine at different positions. Removal of the 6-oxo group of the guanine located adjacent to the target cytosine in the CpG site dramatically reduces the stability of the methyltransferase-DNA complexes and leads to a significant decrease in the methylation. Apparently, O6 of this guanine is involved in the recognition of CpG sites by the enzymes. Cooperative binding of Dnmt3a-CD to 2-aminopurine-containing DNA and the formation of nonproductive enzyme-substrate complexes were observed.     

Guanine-06 methylation reduces the reactivity of d(GpG) towards platinum complexes.

6-methylated guanine dinucleotides were used to study the influence of hydrogen bonding on the specific binding of the antitumor drug cDDP, cis-PtCl2(NH3)2, to DNA. In this interaction, the guanine-06 site appears to be important in explaining the preference for a pGpG-N7(1),N7(2) chelate, which results from H-bridge formation with the ammine ligand of cDDP. Guanine-06 methylated dinucleotides and the nonmodified dinucleotides were reacted with [Pt(dien)Cl]+, cis-PtCl2(NH3)2, and cis-[Pt(NH3)2(H2O)2]2+ and the reaction products were characterized by 1H NMR using pH titrations. Methylation at guanine-06 clearly reduces the preference for the guanine. In competition experiments monitored by NMR and experiments using UV spectrophotometry a decreasing reactivity towards [Pt(dien)(H2O)]2+ and cis-[Pt(NH3)2(H2O)2]2+ was found, in the order of d(GpG) greater than d(GomepG) greater than d(GpGome) greater than d(GomepGome). The difference in reactivity between 5' guanine methylation and 3' guanine methylation is ascribed to differences in the H-bond formation with the backbone phosphate. The resulting reduced stacking of the bases in both modified dinucleotides, compared to the bases in d(GpG), results in a preference for the 3' guanine over 5'.     

The interactions of ruthenium hexaammine with Z-DNA: crystal structure of a Ru(NH3)6+3 salt of d(CGCGCG) at 1.2 A resolution

A crystal of d(CGCGCG) in the Z-DNA lattice was soaked with ruthenium(III) hexaammine and its structure refined at 1.2 A resolution. Three unique metal complexes were found absorbed to each hexamer duplex. In addition, two symmetry-related binding sites were located, yielding a total of five ruthenium complexes bound to each d(CGCGCG) duplex. One unique site and its symmetry related site are nearly identical to the binding site of cobalt(III) hexaammine on Z-DNA. At that position, the metal complex bridges the convex surfaces of two adjacent Z-DNA strands by hydrogen bonds to the N7 and O6 functional groups of the guanine bases. The remaining three ruthenium three ruthenium(III) hexaammine binding sites are not present in the cobalt(III) hexaammine Z-DNA structure. Of these, two are related by symmetry and span the gap between the convex outer surface of one Z-DNA strand and the helical groove crevice of a neighboring strand. The third ruthenium site has no symmetry mate and involves interactions with only the deep groove. In this interaction, the metal complex hydrogen bonds to both the phosphate backbone and to a set of primary shell water molecules that extend the hydrogen bonding potential of the deep groove crevice out to the surface of the molecule. Solution studies comparing the circular dichroism spectra of low salt poly(dG-dC).poly(dG-dC) samples in the presence of ruthenium(III) and cobalt(III) hexammine show that the ruthenium complex does stabilize Z-DNA in solution, but not as effectively as the cobalt analogue. This suggests that some of the interactions available for the larger ruthenium complex may not be important for stabilization of the left-handed DNA conformation.     

Transition state structure of purine nucleoside phosphorylase and principles of atomic motion in enzymatic catalysis

Immucillin-H [ImmH; (1S)-1-(9-deazahypoxanthin-9-yl)-1,4-dideoxy-1,4-imino-D-ribitol] is a 23 pM inhibitor of bovine purine nucleoside phosphorylase (PNP) specifically designed as a transition state mimic [Miles, R. W., Tyler, P. C., Furneaux, R. H., Bagdassarian, C. K., and Schramm, V. L. (1998) Biochemistry 37, 8615-8621]. Cocrystals of PNP and the inhibitor are used to provide structural information for each step through the reaction coordinate of PNP. The X-ray crystal structure of free ImmH was solved at 0.9 A resolution, and a complex of PNP.ImmH.PO(4) was solved at 1.5 A resolution. These structures are compared to previously reported complexes of PNP with substrate and product analogues in the catalytic sites and with the experimentally determined transition state structure. Upon binding, ImmH is distorted to a conformation favoring ribosyl oxocarbenium ion formation. Ribosyl destabilization and transition state stabilization of the ribosyl oxocarbenium ion occur from neighboring group interactions with the phosphate anion and the 5'-hydroxyl of the ribosyl group. Leaving group activation of hypoxanthine involves hydrogen bonds to O6, N1, and N7 of the purine ring. Ordered water molecules provide a proton transfer bridge to O6 and N7 and permit reversible formation of these hydrogen bonds. Contacts between PNP and catalytic site ligands are shorter in the transition state analogue complex of PNP.ImmH.PO(4) than in the Michaelis complexes of PNP.inosine.SO(4) or PNP.hypoxanthine.ribose 1-PO(4). Reaction coordinate motion is dominated by translation of the carbon 1' of ribose between relatively fixed phosphate and purine groups. Purine and pyrimidine phosphoribosyltransferases and nucleoside N-ribosyl hydrolases appear to operate by a similar mechanism.     

Mechanism of CpG DNA Methyltransferases M.SssI and Dnmt3a Studied by DNA Containing 2-Aminopurine

Murine DNA methyltransferases Dnmt3a-CD and M.SssI from Spiroplasma methylate cytosines at CpG sites. The role of 6-oxo groups of guanines in DNA methylation by these enzymes has been studied using DNA substrates, which contained 2-aminopurine at different positions. Removal of the 6-oxo group of the guanine located adjacent to the target cytosine in the CpG site dramatically reduces the stability of the methyltransferase–DNA complexes and leads to a significant decrease in the methylation. Apparently, O6 of this guanine is involved in the recognition of CpG sites by the enzymes. Cooperative binding of Dnmt3a-CD to 2-aminopurine-containing DNA and the formation of nonproductive enzyme–substrate complexes were observed.