Steady-state fluorescence and molecular-modeling studies of tomaymycin-DNA adducts

The interaction of tomaymycin and 8-O-methyltomaymycin with calf thymus DNA was studied by steady-state fluorescence techniques. The 8-phenolic proton of tomaymycin has a pK = 8.0, and the phenolate anion is essentially nonfluorescent. However, the fluorescence of the DNA adduct does not decrease until pH greater than 10.5, when the DNA double helix denatures. Acrylamide quenches the fluorescence of the free antibiotic with a quenching rate constant kq = 7 x 10(9) M-1 s-1. In DNA adducts, the quenching rate constant is reduced about 50-fold, indicating that the aromatic ring of the drug is shielded from the solvent. The four possible binding modes of the antibiotics were modeled on a 6-mer duplex by molecular mechanics calculations in the absence and presence of water and counterions. The modeling studies show that the antibiotic is buried in the minor groove in all binding modes, with the 8-substituent pointing away from the DNA core. Three or five waters are displaced from the minor groove, depending on the orientation of the drug on the DNA.     

Structure of a beta-alanine-linked polyamide bound to a full helical turn of purine tract DNA in the 1:1 motif

Polyamides composed of N-methylpyrrole (Py), N-methylimidazole (Im) and N-methylhydroxypyrrole (Hp) amino acids linked by beta-alanine (beta) bind the minor groove of DNA in 1:1 and 2:1 ligand to DNA stoichiometries. Although the energetics and structure of the 2:1 complex has been explored extensively, there is remarkably less understood about 1:1 recognition beyond the initial studies on netropsin and distamycin. We present here the 1:1 solution structure of ImPy-beta-Im-beta-ImPy-beta-Dp bound in a single orientation to its match site within the DNA duplex 5'-CCAAAGAGAAGCG-3'.5'-CGCTTCTCTTTGG-3' (match site in bold), as determined by 2D (1)H NMR methods. The representative ensemble of 12 conformers has no distance constraint violations greater than 0.13 A and a pairwise RMSD over the binding site of 0.80 A. Intermolecular NOEs place the polyamide deep inside the minor groove, and oriented N-C with the 3'-5' direction of the purine-rich strand. Analysis of the high-resolution structure reveals the ligand bound 1:1 completely within the minor groove for a full turn of the DNA helix. The DNA is B-form (average rise=3.3 A, twist=38 degrees ) with a narrow minor groove closing down to 3.0-4.5 A in the binding site. The ligand and DNA are aligned in register, with each polyamide NH group forming bifurcated hydrogen bonds of similar length to purine N3 and pyrimidine O2 atoms on the floor of the minor groove. Each imidazole group is hydrogen bonded via its N3 atom to its proximal guanine's exocyclic amino group. The important roles of beta-alanine and imidazole for 1:1 binding are discussed.     

The determinants of carboxyl pKa values in turkey ovomucoid third domain

A computational methodology for protein pK(a) predictions, based on ab initio quantum mechanical treatment of part of the protein and linear Poisson-Boltzmann equation treatment of the bulk solvent, is presented. The method is used to predict and interpret the pK(a) values of the five carboxyl residues (Asp7, Glu10, Glu19, Asp27, and Glu43) in the serine protease inhibitor turkey ovomucoid third domain. All the predicted pK(a) values are within 0.5 pH units of experiment, with a root-mean-square deviation of 0.31 pH units. We show that the decreased pK(a) values observed for some of the residues are primarily due to hydrogen bonds to the carboxyl oxygens. Hydrogen bonds involving amide protons are shown to be particularly important, and the effect of hydrogen bonding is shown to be nonadditive. Hydrophobic effects are also shown to be important in raising the pK(a). Interactions with charged residues are shown to have relatively little effect on the carboxyl pK(a) values in this protein, in general agreement with experiment.     

Investigations of Sso7d catalytic residues by NMR titration shifts and electrostatic calculations

Sso7d is a small basic protein consisting of 62 amino acids isolated from the thermoacidophilic archeobacterium Sulfolobus solfataricus. The protein is endowed with DNA binding properties, RNase activity, and the capability of rescuing aggregated proteins in the presence of ATP. In this study, the electrostatic properties of Sso7d are investigated by using the Poisson-Boltzmann calculation of the surface potential distribution and following by NMR spectroscopy the proton chemical shift pH titration of acidic residues. Although the details of the catalytic mechanism still have to be defined, the results from NMR experiments confirm the possible involvement of Glu35 as the proton acceptor in the catalytic reaction, as seen by its abnormally high pK(a) value. Poisson-Boltzmann calculations and NMR titration shifts suggest the presence of a possible hydrogen bond between Glu35 and Tyr33, with a consequent rather rigid arrangement at these positions. Comparison with RNase T1 suggests that Tyr7 may be a good candidate for acting as a proton donor in the active site of Sso7d as shown by its low phenolic pK(a) of approximately 9.3. Titration experiments performed with the UpA, a RNA dinucleotide model, showed that the protein residues affected by the interaction are mainly located in a different region with respect to the surface affected by DNA recognition, in good agreement with the surface potential distribution found with electrostatic calculations.     

A buried lysine that titrates with a normal pKa: role of conformational flexibility at the protein-water interface as a determinant of pKa values

Previously we reported that Lys, Asp, and Glu residues at positions 66 and 92 in staphylococcal nuclease (SNase) titrate with pK(a) values shifted by up to 5 pK(a) units in the direction that promotes the neutral state. In contrast, the internal Lys-38 in SNase titrates with a normal pK(a). The crystal structure of the L38K variant shows that the side chain of Lys-38 is buried. The ionizable moiety is approximately 7 A from solvent and ion paired with Glu-122. This suggests that the pK(a) value of Lys-38 is normal because the energetic penalty for dehydration is offset by a favorable Coulomb interaction. However, the pK(a) of Lys-38 was also normal when Glu-122 was replaced with Gln or with Ala. Continuum electrostatics calculations were unable to reproduce the pK(a) of Lys-38 unless the protein was treated with an artificially high dielectric constant, consistent with structural reorganization being responsible for the normal pK(a) value of Lys-38. This reorganization must be local because circular dichroism and NMR spectroscopy indicate that the L38K protein is native-like under all conditions studied. In molecular dynamics simulations, the ion pair between Lys-38 and Glu-122 is unstable. The simulations show that a minor rearrangement of a loop is sufficient to allow penetration of water to the amino moiety of Lys-38. This illustrates both the important roles of local flexibility and water penetration as determinants of pK(a) values of ionizable groups buried near the protein-water interface, and the challenges faced by structure-based pK(a) calculations in reproducing these effects.     

Ionization properties of mixed lipid membranes: a Gouy-Chapman model of the electrostatic-hydrogen bond switch

The dissociation state of phosphatidic acid (PA) in a lipid bilayer is governed by the competition of proton binding and formation of a hydrogen bond through a mechanism termed the electrostatic-hydrogen bond switch. This mechanism has been suggested to provide the basis for the specific recognition of PA by proteins. Even in bare lipid bilayers the electrostatic-hydrogen bond switch is present if the membrane contains lipids like phosphatidylethanolamine that act as hydrogen bond donors. In this case, the dissociation state (pK(a)) of PA depends strongly on membrane composition. In the present work we incorporate the electrostatic-hydrogen bond switch mechanism into the Gouy-Chapman model for a membrane that is composed of PA and a hydrogen bond-donating zwitterionic lipid. To this end, our model integrates into the Gouy-Chapman approach a recently suggested electrostatic model for zwitterionic lipids. Hydrogen bond formation is incorporated phenomenologically as an additional non-electrostatic interaction between the phosphomonoester headgroup of PA and the zwitterionic lipid headgroup. We express the energetics of the composite membrane in terms of a free energy functional whose minimization leads to a modified non-linear Poisson-Boltzmann equation that we have solved numerically. Our calculations focus on the influence of the membrane environment on the dissociation state of PA. This influence can be expressed as a shift of the second pK(a) of PA, which we calculate as function of membrane composition, following experimental observation. The shift is large and negative if PA is the minor component in the membrane, and it changes over four pH units as function of the mole fraction of PA in the membrane. In contrast, the shift of the second pK(a) of PA remains small and is always positive if the zwitterionic lipid is unable to act as hydrogen bond donor. Hence, we find that the electrostatic-hydrogen bond switch mechanism regulates the dissociation state of PA with much greater sensitivity than would be possible based on a pure electrostatic regulation through the membrane potential.     

Liquid phase method of synthesizing polyamides binding DNA using oligocarboxamide blocks

A method was developed for the synthesis of sequence-specific polyamides on the basis of 4-amino-1-methylpyrrole-2-carboxylic acid, 4-amino-1-methylimidazole-2-carboxylic acid, beta-alanine, and gamma-aminobutyric acid. Dimeric and trimeric oligocarboxamides were used as building blocks. Our synthetic scheme was applied for the synthesis of DNA minor groove binders containing up to twelve carboxamide units. The English version of the paper: Russian Journal of Bioorganic Chemistry, 2004, vol. 30, no. 1; see also http://www.maik.ru.     

Analysing DNA complexes by circular and linear dichroism

The application of linear and circular dichroism (LD and CD) in nucleic acid research id illustrated by recent results aimed at answering specific structural problem in the interaction of DNA with molecules of biological importance. We first consider the circumstances under which ligands, such as DAPI (4′, 6-diamidino-2-phenylindole), change their preferred binding mode in the minor groove to major groove binding or intercalation. As an extension of this problem we refer to the switch between groove binding and intercalation of structurally similar ligands such as ellipticines and trigonal ruthenium complexes. We also explore the use of LD and CD in the determination of the structure of the complex formed between the polynucleotide poly(dA) and the novel ‘peptide nucleic acid’, consisting of nucleic acid bases joined by a polyamide homomorphous with the deoxyribose-phosphate backbone of DNA. Finally, the structure and interaction of the recombination enzyme RecA with DNA is discussed, in particular the influence of the presence of the intercalators, groove binders or covalent DNA adducts.     

Kinetics of DNA alkylation, depurination and hydrolysis of anti diol epoxide of benzo(a)pyrene and the effect of cadmium on DNA alkylation

Anti benzo[a]pyrene diol epoxide (BPDE) alkylates guanines of DNA at N7 in the major groove and at the exocyclic amino group in the minor groove. In this report we investigated the rates of BPDE hydrolysis, DNA alkylation and subsequent depurination of BPDE-adducted pBR322 DNA fragment using polyacrylamide gel electrophoresis. Preincubation studies showed that it hydrolyzed completely in triethanolamine buffer in <2 min. The depurination kinetics showed that a fraction of the N7 alkylated guanine depurinated rapidly; however a significant amount of N7 guanine alkylation remained stable to spontaneous depurination over a 4-h period. Similar results were obtained for the hydrolysis and alkylation rates of syn isomer but it required nearly 500 times more concentration to induce similar levels of N7 guanine alkylation. Cadmium ion strongly inhibited the N7 guanine alkylation of both isomers. But the minor groove alkylation was not affected as demonstrated by postlabeling assay which confirmed the presence of heat-and cadmium-stable minor groove adducts in BPDE-treated calf thymus DNA. Based on these and our earlier findings, we propose a mechanism for the synergistic effect of cadmium in chemically induced carcinogenesis.     

GABA transaminase inhibitors enhance the release of endogenous GABA but decrease the release of beta-alanine evoked by electrical stimulation of slices of the rat medulla oblongata

Slices of the rat medulla oblongata were superfused and electrically stimulated. The amount of endogenous GABA, beta-alanine and glutamate release from the slices was determined by high performance liquid chromatography with fluorometric detection. Inhibitors of GABA-transaminase (GABA-T), aminooxyacetic acid (10(-5) M), gamma-acetylenic GABA (10(-4) and 10(-3) M) and gabaculine (10(-5) M), enhanced the stimulus-evoked release of GABA and reduced that of beta-alanine, while no change was observed in the release of glutamate. These changes in amino acid release from the slices were accompanied by an increase in the content of GABA and a decrease in that of beta-alanine. The stimulus-evoked release of these amino acids was abolished by Ca2+-deprivation, in either the presence or absence of GABA-T inhibitors. These results suggest a modulatory role of GABA-T for synaptically releasable GABA and involvement of this enzyme in the synthesis of releasable beta-alanine.     

Oligonucleotide--minor groove binder 1:2 conjugates: side by side parallel minor groove binder motif in stabilization of DNA duplex

Synthetic polycarboxamides consisting of N-methylpyrrole (Py), N-methylimidazole (Im), N-methyl-3-hydroxypyrrole (Hp) and beta-alanine (beta) show strong and sequence-specific interaction with the DNA minor groove when they form hairpin structures with side-by-side antiparallel motifs. In the present paper, new conjugates containing two ligands linked to the same terminal phosphate of DNA strand were constructed. The paper describes optimized synthesis and properties of oligonucleotide-linked polyamide strands that insert into the minor groove of a duplex in a parallel or antiparallel orientation. Strong stabilization of DNA duplexes by two attached minor groove ligands is demonstrated by the thermal denaturation method. The unmodified duplex 5'-CGTTTATTp-3'/5'-AATAAACG-3' melts at 20 degrees C. When one tetra(Py) residue was attached to the first strand of this duplex, denaturation temperature was increased to 46 degrees C; attachment of the second tetra(Py) in a parallel orientation resulted in denaturation temperature of 60 degrees C. It is even higher than in case of "classic" octapyrrole hairpin ligand (Tm = 58 degrees C). Sequence-specific character of stabilization by two conjugated ligands was demonstrated for G:C-containing oligonucleotides attached to tetracarboxamide and octacarboxamide ligands constructed from Py, Im and beta units according to established recognition rules (deltaTm = 20 degrees C). The two-strand parallel minor groove binder constructions attached to addressing oligonucleotides could be considered as site-specific ligands recognizing single- and double-stranded DNA similarly to already described hairpin MGB structures with antiparallel orientation of carboxamide units.     

A new protein domain for binding to DNA through the minor groove.

Protein p6 of the Bacillus subtilis phage phi 29 binds with low sequence specificity to DNA through the minor groove, forming a multimeric nucleoprotein complex that activates the initiation of phi 29 DNA replication. Deletion analysis suggested that the N-terminal part of protein p6, predicted to form an amphipathic alpha-helix, is involved in DNA binding. We have constructed site-directed mutants at the polar side of the putative alpha-helix. DNA binding and activation of initiation of phi 29 DNA replication were impaired in most of the mutant proteins obtained. A 19 amino acid peptide comprising the N-terminus of protein p6 interacted with a DNA fragment containing high-affinity signals for protein p6 binding with approximately 50-fold higher affinity than the peptide corresponding to an inactive mutant. Both wild-type peptide and protein p6 recognized the same sequences in this DNA fragment. This result, together with distamycin competition experiments, suggested that the wild-type peptide also binds to DNA through the minor groove. In addition, CD spectra of the wild-type peptide showed an increase in the alpha-helical content when bound to DNA. All these results indicate that an alpha-helical structure located in the N-terminal region of protein p6 is involved in DNA binding through the minor groove.     

Different effects on human topoisomerase I by minor groove and intercalated deoxyguanosine adducts derived from two polycyclic aromatic hydrocarbon diol epoxides at or near a normal cleavage site

Topoisomerase I (top1) relieves supercoiling in DNA by forming transient covalent cleavage complexes. These cleavage complexes can accumulate in the presence of damaged DNA or anticancer drugs that either intercalate or lie in the minor groove. Recently we reported that covalent diol epoxide (DE) adducts of benzo[a]pyrene (BaP) at the exocyclic amino group of G(+1) block cleavage at a preferred cleavage site ( approximately CTT-G(+1)G(+2)A approximately ) and cause accumulation of cleavage products at remote sites. In the present study, we have found that the 10S G(+2) adduct of BaP DE, which lies toward the scissile bond in the minor groove, blocks normal cleavage, whereas the 10R isomer, which orients away from this bond, allows normal cleavage but blocks religation. In contrast to BaP, the pair of benzo[c] phenanthrene (BcPh) DE adducts at G(+2), which intercalate from the minor groove either between G(+1)/G(+2) or between G(+2)/A, allow normal cleavage but block religation. Both intercalated BcPh DE adducts at G(+1) suppress normal cleavage, as do both groove bound BaP DE adducts at this position. These studies demonstrate that these DE adducts provide a novel set of tools to study DNA topoisomerases and emphasize the importance of contacts between the minor groove and top1's catalytic site.     

Interaction of synthetic analogs of distamycin with DNA. Role of the conjugated N-methylpyrrole system in specificity of binding

Interaction of two synthetic analogs of distamycin (Dst), PPA and PAP, containing a saturated beta-alanine moiety replacing one N-methylpyrrole ring, with different polynucleotides and natural DNAs were studied using UV and CD spectroscopy. The results indicate that, similar to Dst, these analogs bind to DNA via the minor groove with a specificity towards AT-base pairs. It may be proposed that pyrrole chromophores in Dst probably do not play a role in the AT-base selectivity exhibited by Dst.     

Responses to GABA, glycine and beta-alanine induced in Xenopus oocytes by messenger RNA from chick and rat brain

Poly (A)+ messenger RNA (mRNA) was extracted from rat and chick brains, and injected into oocytes of Xenopus laevis. This led to the expression of receptors that evoked membrane currents in response to gamma-aminobutyric acid (GABA), glycine and beta-alanine. These currents all inverted at about the chloride equilibrium potential in the oocyte, and showed a marked rectification at negative potentials. Oocytes injected with mRNA from chick optic lobe gave large responses to GABA and beta-alanine, but small responses to glycine. In contrast, one fraction of mRNA from rat cerebral cortex (obtained by sucrose density gradient centrifugation) caused oocytes to develop sensitivity to GABA, glycine and beta-alanine, but very little to GABA. The pharmacological properties of the three amino acid responses also differed. Barbiturate and benzodiazepines potentiated the responses to GABA and beta-alanine, but not to glycine. Strychnine reduced the responses to glycine and beta-alanine, but not to GABA, whereas bicuculline reduced the responses to GABA and beta-alanine, but not to glycine. We conclude that different species of mRNA code for receptors to GABA and glycine, and possibly also for separate beta-alanine receptors.     

Property of the M-DNA probed by a minor groove binding dye 4',6-diamidino-2-phenylindole

Spectral properties including circular and linear dichroism (CD and LD) of M-DNA, a molecular electric wire, formed at a high Zn(2+) concentration have been studied using a minor groove binding drug 4',6-diamidino-2-phenylindole (DAPI) as a probe. As the Zn(2+) concentration increased, the magnitude of LD in the DNA absorption region decreased at pH 8.5, implying the aggregation of DNA, which is in contrast with the retained LD magnitude at pH 7.0. As the M-DNA formed, change in the secondary structure of DNA was observed by CD spectrum, which resembles that of the C-form DNA, although overall structure seemed to remain as a right handed double helix. The DAPI-DNA complex in the presence of high concentration of Zn(2+) ions at pH 7.0 exhibited the similar CD spectrum with that in the absence of Zn(2+) ion, consisting of type I, II and III. In contrast, at pH 8.5 at a high Zn(2+) concentration in which DNA is in its M-form, DNA bound DAPI produced only the type III CD, suggesting that DAPI binds at the surface of M-DNA: the presence of Zn(2+) ions prevents the minor groove binding of DAPI.     

Salt effects on ionization equilibria of histidines in myoglobin

The salt dependence of histidine pK(a) values in sperm whale and horse myoglobin and in histidine-containing peptides was measured by (1)H-NMR spectroscopy. Structure-based pK(a) calculations were performed with continuum methods to test their ability to capture the effects of solution conditions on pK(a) values. The measured pK(a) of most histidines, whether in the protein or in model compounds, increased by 0.3 pH units or more between 0.02 M and 1.5 M NaCl. In myoglobin two histidines (His(48) and His(36)) exhibited a shallower dependence than the average, and one (His(113)) showed a steeper dependence. The (1)H-NMR data suggested that the salt dependence of histidine pK(a) values in the protein was determined primarily by the preferential stabilization of the charged form of histidine with increasing salt concentrations rather than by screening of electrostatic interactions. The magnitude and salt dependence of interactions between ionizable groups were exaggerated in pK(a) calculations with the finite-difference Poisson-Boltzmann method applied to a static structure, even when the protein interior was treated with arbitrarily high dielectric constants. Improvements in continuum methods for calculating salt effects on pK(a) values will require explicit consideration of the salt dependence of model compound pK(a) values used for reference in the calculations.     

The electrostatic potential of B-DNA

Electrostatic potentials around DNA are obtained by solving the nonlinear Poisson-Boltzmann (PB) equation. The detailed charge distribution of the DNA and the different polarizabilities of the macromolecule and solvent are included explicitly in the calculations. The PB equation is solved using extensions of a finite difference approach applied previously to proteins. Electrical potentials and ion concentrations are compared to those obtained with simpler models. It is found that the shape of the dielectric boundary between the macromolecule and solvent has significant effects on the calculated potentials near the surface, particularly in the grooves. Sequence-specific patterns are found, the most surprising result being the existence of positive regions of potential near the bases in both the major and minor grooves. The effect of solvent and ionic atmosphere screening of phosphate-phosphate repulsions is studied, and an effective dielectric function, appropriate for molecular mechanics simulations, is derived.     

DAPI binding to the DNA minor groove: a continuum solvent analysis

A continuum solvent model based on the generalized Born (GB) or finite-difference Poisson-Boltzmann (FDPB) approaches has been employed to compare the binding of 4'-6-diamidine-2-phenyl indole (DAPI) to the minor groove of various DNA sequences. Qualitative agreement between the results of GB and FDPB approaches as well as between calculated and experimentally observed trends regarding the sequence specificity of DAPI binding to B-DNA was obtained. Calculated binding energies were decomposed into various contributions to solvation and DNA-ligand interaction. DNA conformational adaptation was found to make a favorable contribution to the calculated total interaction energy but did not change the DAPI binding affinity ranking of different DNA sequences. The calculations indicate that closed complex formation is mainly driven by nonpolar contributions and was found to be disfavored electrostatically due to a desolvation penalty that outbalances the attractive Coulomb interaction. The calculated penalty was larger for DAPI binding to GC-rich sequences compared with AT-rich target sequences and generally larger for the FDPB vs the GB continuum model. A radial interaction profile for DAPI at different distances from the DNA minor groove revealed an electrostatic energy minimum a few Angstroms farther away from the closed binding geometry. The calculated electrostatic interaction up to this distance is attractive and it may stabilize a nonspecific binding arrangement.