Nucleic acid-mutagen interactions: crystal structure of adenylyl-3',5'-uridine plus 9-aminoacridine

Ned Seeman and the prediction of amino acid-basepair motifs mediating protein-nucleic acid recognition

Fifty years ago, the first atomic-resolution structure of a nucleic acid double helix, the mini-duplex (ApU)2, revealed details of basepair geometry, stacking, sugar conformation, and backbone torsion angles, thereby superseding earlier models based on x-ray fiber diffraction, including the original DNA double helix proposed by Watson and Crick. Just 3 years later, in 1976, Ned Seeman, John Rosenberg, and Alex Rich leapt from their structures of mini-duplexes and H-bonding motifs between bases in small-molecule structures and transfer RNA to predicting how proteins could sequence specifically recognize double helix nucleic acids. They proposed interactions between amino acid side chains and nucleobases mediated by two hydrogen bonds in the major or minor grooves. One of these, the arginine-guanine pair, emerged as the most favored amino acid-base interaction in experimental structures of protein-nucleic acid complexes determined since 1986. In this brief review we revisit the pioneering work by Seeman et al. and discuss the importance of the arginine-guanine pairing motif.

RNA–Intercalating Agent Interactions: in vitro Antiviral Activity Studies

Twenty intercalating agents were tested to examine the effects of intercalating dye-induced perturbations upon the antiviral activity of poly (adenylate–uridylate) [poly (A-U)]. Neither poly (A-U) alone nor each intercalative dye was an efficacious antiviral agent. When poly (A-U) was combined with major groove intercalating dyes (acridine orange or proflavine), no synergism was observed. When poly (A-U) was combined with minor groove intercalating dyes [ethidium (EB), propidium (PI), adriamycin (ADR) or daunomycin (DMN)] or minor/major groove intercalating dyes [9-aminoacridine (9-AA), N2-methyl-9-hydroxy-ellipticine (NMHE) or N2,N6-dimethyl-9-hydroxy-ellipticine (DMHE)] the 50% effective doses (ED50) of the poly (A-U), 9-AA, ADR, DMHE, DMN, EB, NMHE and PI decreased 18-, 22-, 60-, 274-, 61-, 154-, 113- and 299-fold, respectively. When poly (A-U) was combined individually with 11 dyes whose mode of intercalation was not known, the ED50 of ametantrone (HAQ), chloroquine (CHL), mitoxantrone (DHAQ) and quinine (QUI) decreased 125-, 65-, 251- and 32-fold, respectively. These results suggest that the four dyes may intercalate into poly (A-U) from the minor groove. Ten (ADR, CHL, DMN, DHAQ, DMHE, EB, HAQ, NMHE, PI, QUI) of the 20 dyes evaluated exhibited significant synergism with poly (A-U), as quantified by the fractional inhibitory concentration index. Interferon (IFN) neutralization assays demonstrated that the IFN-inducing capability of the dye/poly (A-U) combinations approximated the sum of the capabilities of the poly (A-U) and the dyes employed. These results suggest that the majority of the dyes tested potentiate the antiviral activity of poly (A-U) without affecting the amount of IFN induced.

Conformational analysis of GpA and GpAp in aqueous solution by molecular dynamics and statistical methods

Barnase, an extracellular endoribonuclease from Bacillus amyloliquefaciens, hydrolyses single-stranded RNA. Its very low catalytic activity toward GpN dinucleotides, where N stands for any nucleoside, is markedly increased when a phosphate is added to the 3'-end, as in GpNp. Here we investigate the conformational properties of GpA and GpAp in solution, in order to determine whether differences in these properties may be related to the changes in enzymatic activity. Two independent 1.3 ns molecular dynamics trajectories are generated for each dinucleotide in the presence of explicit water molecules and counter ions. These trajectories are analysed by monitoring molecular properties, such as the solvent accessible surface area, the distance and orientation between the bases, the behaviour of torsion angles and formation of intramolecular H-bonds. To identify relevant correlations between these parameters, statistical techniques, comprising multiple regression, clustering and discriminant analysis are used. Results show that GpA has a significant propensity to form folded conformations (approximately 50%), fostered by a small number of intramolecular H-bonds, whereas GpAp remains essentially extended. The latter behaviour seems to be due to an H-bond between the terminal phosphate and adenosine ribose group, which restricts rotation about the adenine Agamma angle. We also find that GpA folding is induced by a concerted motion of specific torsion angles, which is closely coupled to the formation of a network of flexible hydrogen bonds. Finally, on the basis of an expression for barnase KM, which incorporates the folded/extended conformational equilibria of the dinucleotide substrates, it is argued that our findings on the differences between these equilibria, can qualitatively rationalize the experimentally measured differences in enzymatic properties.

Photochemical and photophysical studies of 3-amino-6-iodoacridine and the inactivation of lambda phage

The photochemistry and photophysics of 3-amino-6-iodoacridine (Acr-I) was studied. Photolysis (350 nm) of Acr-I (free base) generates products consistent with a free radical intermediate in methanol, benzene and carbon tetrachloride. The Acr-I hydrochloride is shown to bind to calf thymus DNA and to the self-complementary dinucleotide cytidylyl-(3'-5')-guanosine (CpG) miniduplex in a manner similar to that of proflavine (Acr-NH2), a known DNA intercalator. The Acr-I is shown to more efficiently nick supercoiled plasmid DNA pBR322 upon 350 nm or 420 nm photolysis than Acr-NH2. The efficiency of Acr-I-sensitized DNA nicking is not oxygen dependent. Photolysis of the Acr-I/(CpG)2 complex leads to cleavage of the dinucleotide and to cytidine base release by selective damage to a specific ribose moiety. Dinucleotide cleavage occurs equally well in the presence or absence of oxygen, thereby eliminating a singlet oxygen- or peroxyl radical-mediated process. Photolysis of Acr-I in the presence of a mononucleotide (GMP) or a non-self-complementary dinucleotide (uridylyl-[3'-5']-cytidine-UpC) does not lead to fragmentation and base release. Similarly, photolysis of the Acr-NH2/(CpG)2 complex does not lead to fragmentation and base release. The data indicate that photolysis of an iodinated intercalator bound to CpG or plasmid DNA generates an intercalated aryl radical and that the reactive intermediate initiates a sequence of reactions that efficiently nick nucleic acids. The inactivation of lambda phage sensitized by Acr-I with UV (350 nm) light is oxygen independent but with visible (420 nm) light is strongly oxygen dependent. The Acr-I fluoresces more intensely when excited at 446 than at 376 nm. Thus, UV photolysis may lead to C-I bond homolysis and free radical formation, a process that is not energetically feasible with visible light. The results demonstrate the difficulty of extrapolating model studies involving simple molecules and DNA to understanding the mechanism of viral inactivation with a particular sensitizer.

System for the pH-dependent release of a dye in model dental restorations

We are developing a system for detecting recurrent caries under dental restorations. The controlled release of dyes under conditions of likely demineralization will alert the dentist to potential secondary caries. Production of acidic species is a characteristic of caries activity; hence, the system uses pH-sensitive polymers to release markers when the pH at the cavity wall of the restored tooth is below 6.5. The objectives of this investigation were to test the hypotheses that (1) the proposed system can be designed to release detectable marker continuously for at least six months in a simulated carious environment, and (2) the transient pH changes in the oral cavity caused by simulated dietary intake will not induce premature marker release from the pH-sensitive polymer placed beneath restorations. Two types of dye-loaded microspheres based on styrene, vinylpyridine, and divinylbenzene were prepared and placed on the floor of model cavity preparations made from an acrylic rod. Each model cavity was restored with a hybrid dental composite, placed in a vial with 5 mL of sodium-lactate/lactic-acid base buffer solution, and stored at 37 degrees C. Solutions of three different pH values were used: 2.86, 4.73, and 6.39. The dye release into storage media was monitored periodically with a UV/VIS spectrophotometer. Results showed that the duration could extend beyond six months for pH > 4.73, and that transient oral pH changes are not likely to result in premature dye release. The data indicate that it would take approximately 21 days for the acidic agent external to the restoration to initiate dye release from restored sites.

Ligand effects on platinum binding to DNA. A comparison of DNA binding properties for cis- and trans-[PtCl2(amine)2] (amine = NH3, pyridine)

The DNA binding properties of cis- and trans-[PtCl2(pyridine)2] have been examined and compared with their NH3 analogs, cis- and trans-DDP. The presence of a planar ligand reduces the rates of DNA binding but does not greatly affect the overall conformation of CT DNA, as measured by circular dichroism spectroscopy. The sequence specificity of trans-[PtCl2(py)2] includes alternating purine-pyrimidine sequences. The sequence specificity is further different between the two pyridine isomers, and the steric effects of two cis-pyridine groups are demonstrated by the appearance of relatively few binding sites in the 49-bp duplex. The effects of the pyridine ligand are further manifested by a greatly enhanced DNA-DNA interstrand cross-linking efficiency for the trans isomer, with a cross-link per adduct frequency of between 0.14 and 0.23, depending on the rb of the sample. The unwinding of closed circular pUC19 DNA by trans-[PtCl2(pyridine)2] is also more efficient than that by either DDP isomer, with an unwinding angle calculated at phi = 17 degrees (compare cis-DDP with phi = 13 degrees and trans-DDP with phi = 9-10 degrees). In contrast, little unwinding is induced by cis-[PtCl2(pyridine)2], with phi = 4 degrees. These results in particular invert the standard cis/trans structure-activity relationships observed previously for [PtCl2(NH3)2]. The results are discussed with respect to the previously demonstrated effect of activation of the trans-platinum geometry using sterically hindered ligands.

Correlation of frameshift mutagenicity with DNA intercalation by CGS 20928A using an in vitro DNA unwinding assay

A compound's mutagenicity in different Salmonella tester strains can suggest its mechanism of reaction with DNA. Clear confirmation of such a mechanism, however, requires a direct test of the compound's reaction with DNA, often relying on specific in vitro studies. We report the use of a rapid in vitro test designed to measure DNA unwinding, a characteristic of DNA intercalators and many frameshift mutagens. CGS 20928A, an adenosine antagonist, produced a significant (> 2-fold) increase in revertants only for Salmonella tester strain TA1537, and only without metabolic activation. These data indicated that the compound was a direct acting frameshift mutagen and possibly intercalated into DNA. Our DNA unwinding assay indicated that at concentrations of > 0.1 mM CGS 20928A behaved like known intercalating compounds in that it unwound DNA. These concentrations of compound are comparable to those found mutagenic to TA1537. By comparison, the frameshift mutagen and known intercalating compound 9-aminoacridine unwound DNA in this assay in a concentration dependent fashion between 6-12 microM. ICR-191, another acridine frameshift mutagen, also unwound DNA. A compound structurally related to CGS 20928A, which was not mutagenic in Salmonella tester strains, did not produce any DNA unwinding even at 10 mM. Because the assay uses microgram quantities of material, it should be ideal for screening small amounts of congeneric series suspected of frameshift mutagenicity.

A structural model for sequence-specific proflavin-DNA interactions during in vitro frameshift mutagenesis

Molecular models describing intermediates that may lead to proflavin-induced 1 bp deletions during in vitro polymerization by E. coli DNA polymerase I Klenow fragment are proposed. The models provide structural explanations for the fact that the induced frameshifts always occur opposite template bases that are adjacent to 5' pyrimidines and are based on the underlying hypothesis that the deletions arise because the polymerase passes by a template base without copying it. Because the most frequent mutations are opposite Pu in the template sequence 5' Py Pu 3', a single-strand loop-out model was constructed for this sequence and proflavin was added, using structures found in crystalline oligonucleotides and their complexes with proflavin. The model seeks to rationalize the roles of the 5' pyrimidine and proflavin in facilitating the bypass. Four potential roles for proflavin in mutagenesis are described: 1) stacking on the looped-out base; 2) stacking on the base pair immediately preceding the site of mutation; 3) hydrogen bonding with the 5' pyrimidine; 4) hydrogen bonding with the phosphate backbone. These models point to the possibility that a number of proflavin-DNA interactions may be involved. In contrast, modeling does not suggest a role for classically intercalated proflavin in frameshift mutagenesis arising during in vitro DNA polymerization.

Conformational characteristics of mixed sugar puckered deoxydinucleoside triphosphate units d-pCpGp and d-pGpCp from energy minimization studies

The deoxydinucleoside triphosphate units d-pCpGp and d-pGpCp were subjected to a rigorous theoretical investigation with a view to describing their distinctive conformational characteristics. For each unit 216 probable three-dimensional forms defined by the backbone-base dihedral angles and sugar pucker modes were considered for conformational energy minimization process and scrutinized with reference to properties, such as base-stacking, hydrogen-bonding, internal flexibility and base sequence-phosphate influence. The P-O bond torsions and the phosphate groups were treated with special attention. The results reveal a number of preferred conformational states other than the known helical forms, such as, A-, B-, C-, Z-, and Watson-Crick conformation. Many interesting one-step (change in only one of the dihedral angles or sugar puckers) conformational transitions which involve just about a kcal/mol of energy came to light. The two base sequences CG and GC were noted to differ strikingly in many of their conformational characteristics.

Intercalation model for DNA-cross linking in a 1-nitro-9-aminoacridine derivative, an analog of the antitumor agent "ledakrin" (nitracrine)

Ledakrin (nitracrine), C-283, is a 1-nitro-9-aminoacridine derivative that is used in Poland as an antitumor agent. In order to investigate the basis of the activity of this compound the structure of another analog, [9-(3-dimethyl-1-methylpropylimino)-1-nitro-9,10-dihydroacridin e], C-829, that has similar activity, was determined by X-ray crystallographic techniques and was compared with that of ledakrin, already reported in the literature. In both molecules the proximity of the 1-nitro to the substituted 9-aminoacridine group causes extensive distortions. These compounds are believed to act, after metabolic "activation", by cross-linking DNA. Such cross-linking does not occur in the absence of the 1-nitro group or if the nitro group is moved to the 2-, 3- or 4-position. Computer-assisted model-building has been used to test possible intercalative models. It has shown that functional groups on C-829 or C-283 are, when the acridine portion of the molecule is intercalated as in a proflavine dinucleoside phosphate complex, in positions suitable for DNA cross-linking by activated 1-nitro-9-aminoacridine derivatives.

Structure of a novel drug-nucleic acid crystalline complex: 1, 10-phenanthroline-platinum (II) ethylenediamine--5'-phosphoryl-thymidylyl(3'-5') deoxyadenosine

1,10-Phenanthroline-platinum (II) ethylenediamine (PEPt) forms a 1:2 crystalline complex with 5'-phosphorylthymidylyl (3'-5') deoxyadeno sine (d-pTpA). Crystals are monoclinic, P2, with a = 10.204 A, b = 24.743 A, c = 21.064 A, Beta = 94.6 degrees. The structure has been determined by Patterson and Fourier methods, and refined by least squares to a residual of 0.128 on 2,367 observed reflections. PEPt molecules form sandwich-like stacks with adenine-thymine hydrogen-bonded pairs along the alpha axis. Intercalation in the classic sense is not observed in this structure. Instead, d-pTpA molecules form an open chain structure in which adenine-thymine residues hydrogenbond together with the reversed Hoogsteen type base-pairing configuration. Deoxyadenosine residues exist in the syn conformation and are C3' endo and C1' exo. Thymidine residues are in the high anti conformation with C2' endo puckers. The structure is heavily hydrated, forming a channel-like water network along the alpha axis. Other features of the structure are described.

Structural and sequence-dependent aspects of drug intercalation into nucleic acids

Information gained from X-ray crystallographic studies on drug-nucleic acid complexes is described, with emphasis on the intercalation process. Relevant data from NMR experiments are examined in order to highlight similarities and differences between solution and solid-state structures. Theoretical analyses of intercalation complexes are also discussed and evaluated, with respect to the structural methods, with special reference being made to nucleic acid conformation and positions of drug molecules in the binding sites.

Correlation between glycosyl torsion angle and sugar ring pucker does not always exist

A survey of the conformational parameters of deoxypyrimidine nucleosides and nucleotides shows that the correlation between glycosyl torsion angle and sugar pucker, which has often been considered to be well-established, does not always exist. This may be of significance when interactions between DNA and other molecules are considered.

The molecular structure of d(ICpCpGpG), a fragment of right-handed double helical A-DNA

The DNA tetramer d(ICpCpGpG) or ICCGG crystallizes as a double-stranded 4-base pair (bp) segment of an A helix. Two such tetramer helices are packed together in the crystal with local helix axes nearly coincident, simulating an 8-bp helix, and four such octamers make up the tetragonal unit cell. Restrained energy and reciprocal space refinement has led to an R factor of 20.5% at 2.1 A resolution. The ICCGG helix has a twist corresponding to 10.7 bp per turn, a 19 degree base tilt and a 2.3 A rise per base pair along the helix axis. The mean propeller twist of 18 degree is comparable with, and has the same rotational sense as that observed in the B-DNA dodecamer CGCGAATTCGCG at similarly high alcohol concentration. Backbone phosphate groups in A-DNA are extensively hydrated, including a network across the opening of the major groove, whereas base edge N and O groups in major and minor grooves are less hydrated than in B-DNA. The minor groove spine of hydration observed in B-DNA is totally absent. These observations of relative hydration confirm and extend the model for the B- to A- helix transition proposed earlier on the basis of the B helix structure.

Conformational characteristics of dimeric subunits of RNA from energy minimization studies. Mixed sugar-puckered ApG, ApU, CpG, and CpU

Following the procedure described in the preceding article, the low energy conformations located for the four dimeric subunits of RNA, ApG, ApU, CpG, and CpU are presented. The A-RNA type and Watson-Crick type helical conformations and a number of different kinds of loop promoting ones were identified as low energy in all the units. The 3E-3E and 3E-2E pucker sequences are found to be more or less equally preferred; the 2E-2E sequence is occasionally preferred, while the 2E-3E is highly prohibited in all the units. A conformation similar to the one observed in the drug-dinucleoside monophosphate complex crystals becomes a low energy case only for the CpG unit. The low energy conformations obtained for the four model units were used to assess the stability of the conformational states of the dinucleotide segments in the four crystal models of the tRNAPhe molecule. Information on the occurrence of the less preferred sugar-pucker sequences in the various loop regions in the tRNAPhe molecule has been obtained. A detailed comparison of the conformational characteristics of DNA and RNA subunits at the dimeric level is presented on the basis of the results.

Theoretical calculations of conformational preferences in dinucleoside monophosphates Up-U and Ap-A. Significance of intramolecular base-backbone interaction

Conformational preferences of dinucleoside monophosphates Up-U and Ap-A have been studied theoretically using the Quantum Chemical PCILO (Perturbative Configuration Interaction Using Localized Orbitals) method. The importance of intramolecular interactions between the bases and ribose phosphate backbone has been indicated. Additionally, a close similarity is predicted between overall Up-U and Ap-A conformational structures. Specific prefrences about the key conformation bonds are C4'-C5' (psi = 60 degrees), C5'-O5' (phi = 180 degrees), O5'-P (omega = 110 degrees), P-O3' (omega' = 310 degrees), O3'-C3' (phi' = 210 degrees), C1'-N (chi = 20 degrees), C2'-O2' (theta = 300 degrees) for nucleotidyl units with C3' -endo ribose pucker. Except for the O5'-P torsion angles (omega = 110 degrees), the remaining values fall in ranges observed for monomers and polynucleotides. Stabilization for this unique omega = 110 degrees preference is provided by intramolecular hydrogen bonding between the 3'-OH (hydroxyl group) and O2 (in uracil) or N3 (in adenine) of the second base from the 3'-terminal. The possibility of such interaction mechanisms at 3'-terminals of ribonucleic acids and polynucleotides is also discussed.

Spectroscopic properties of ethidium monoazide: a fluorescent photoaffinity label for nucleic acids

The non-covalent binding of ethidium monoazide to nucleic acids is entirely analogous to that of ethidium (binding constant approximately 2-3 X 10(5) M). The ethidium monoazide can be photochemically covalently linked to nucleic acids in high yield, up to 75%, by long wavelength light. The fluorescence of ethidium monoazide and ethidium crosslinked to nucleic acids show the same environmental sensitivity as does the fluorescence of ethidium. These properties of ethidium monoazide indicate its use as a fluorescent photoaffinity label for nucleic acids. Ethidium diazide can be photochemically linked to nucleic acids but appears to have properties substantially different from those of ethidium.

A 1:2 crystalline complex of ApA:proflavine: a model for binding to single-stranded regions in RNA

The structure of a 1"2 complex of adenylyl-(3',5')-adenosine phosphate and proflavine hemisulfate has been determined using the methods of x-ray crystallography. Since the ApA does not form a mini double helix, it may serve as a model for the interaction of planar molecules with single stranded nucleic acids. The dinucleotide adopts an extended conformation with the adenines in adjacent molecules forming base pairs. A most unusual feature of the molecule is that it does not obey the "rigid nucleotide" concept although none of the torsion angles occur in energetically unfavourable regions. This is most probably due to the strong interactions between the proflavine and the oligonucleotide.

Spatial configuration of deoxyribotrinucleoside diphosphates in aqueous solution

The detailed conformational features and dynamics of the naturally occurring deoxyribotrinucleoside diphosphates d-TpTpT and d-TpTpC have been investigated at 20 degrees C and 80 degrees C in aqueous solution by nuclear magnetic resonance spectroscopy. The observed NMR parameters indicate that the conformational properties of the trimers are very similar to those of the constituent dimers and monomers, i.e., the monomers and dimers conserve their intrinsic conformational features when they become incorporated into oligomers. Model building indicate that the distant shieldings can originate from spatial configurations in which the central nucleotidyl unit is bulged out and the w'1w1, w'2w2 occupy /g+g+, g+g+/ domains.

Structure of mutagen nucleic acid complexes in solution. Proton chemical shifts in 9-aminoacridine complexes with dG-dC, dC-dG, and dA-dT-dG-dC-dA-dT

The influence of self-complementary oligodeoxynucleotides on the chemical shifts of protons of the mutagenic acridine dye 9-aminoacridine has been measured. Upfield shifts indicative of intercalative binding are found in the cases of dG-dC, dC-dG, and dA-dT-dG-dC-dA-dT but not in dA-dT. Geometries for the complexes that are compatibile with the chemical-shift data and the X-ray structure of the complex between ri5C-rG and 9-aminoacridine determined by Sakore et al. [Sakore, T.D., Jain, S.C., Tsai, C., and Sobell, H.M. (1977), Proc. Natl. Acad. Sci. U.S.A. 74, 188--192] can be identified using recent theoretical estimates of shifts induced by nucleotide bases.

Origins of the specificity in the intercalation of ethidium into nucleic acids. A theoretical analysis

The sequence preferences observed in the intercalative binding of ethidium to dinucleoside phosphates have been theoretically examined. This specificity is for pyrimidine (3'--5') purine sequences as compared to their purine (3'--5') pyrimidine sequence isomers. It is shown that the stacking energies between the ethidium cation and the base pairs are fairly constant for all combinations of bases at the intercalation site. In contrast, the energy of unwinding the double helix to assume the geometry of the intercalation complex shows substantial sequence differences. Thus, the specificity observed is more readily explained in terms of these conformation energy changes than by preferential stacking interactions. These results imply that there may be a large class of intercalating molecules which exhibit similar pyrimidine (3'--5') purine sequence specificity.

Stereochemical model for proflavin intercalation in A-DNA

Linked-atom molecular modelling was employed to determine the steric and torsional requirements for intercalation of proflavine into a double-stranded region of DNA compatible with adjacent regions of cohelical A-DNA. The optimum intercalation conformation is characterized by the dihedral angles xi and psi becoming trans, with all sugars retaining the characteristics C3'-endo pucker. This extended conformation results in virtually no helical unwinding, suggesting it may be an appropriate model for an intercalative intermediary in mutagenesis by virtue of its similarity to standard helical DNA.

RNA structure

This review is concerned primarily with the physical structure and changes in the structure of RNA molecules. It will be evident that we have not attempted comprehensive coverage of what amounts to a vast literature. We have tried to stay away from particular areas that have been recently reviewed elsewhere. Citations to and information from them are included, however, so that access to the literature is available. Much of what we treat in depth deals with the crystal structures and solution behaviour of model RNA compounds, including synthetic polymers and molecular fragments such as dinucleoside phosphates. Sequence data on natural RNA are cited, but not in detail. Similarly, apart from tRNA, natural RNAs the structural determinations of which are presently not so far advanced, are not dwelt upon. We have tried to present in detail the available structural data with scaled drawings that permit facile comparisons of molecular geometries.