Interaction of unfused tricyclic aromatic cations with DNA: a new class of intercalators

Unfused tricyclic aromatic ring systems 1-6 with one or two cationic side chains have been synthesized and their interactions with DNA and synthetic polymers probed with a variety of techniques. Molecular mechanics calculations indicate that the torsional angle between ring planes in the minimum energy conformation of the tricyclic molecules can range from 0 degree to as high as 50 degrees depending on the type of rings and substituents. Viscometric titrations with linear and supercoiled DNA, linear dichroism, and NMR studies indicated that all compounds with torsional angles of approximately 20 degrees or less bind to DNA by intercalation. The more highly twisted intercalators caused significant perturbation of DNA structure. Unfused intercalators with twist angles of approximately 20 degrees have reduced binding constants, suggesting that they could not form an optimum interaction with the DNA base pairs. Unfused intercalators with twist less than 20 degrees formed strong complexes with DNA. The structures of these unfused intercalators are more analogous to typical groove-binding molecules, and an analysis of their interaction with DNA provides a better understanding of the subtle differences between intercalation and groove-binding modes for aromatic cations. The results indicate that intercalation and groove-binding modes should be viewed as two potential wells on a continuous energy surface. The results also suggest design strategies for intercalators that can optimally complement DNA base pair propeller twist or that can induce bends in DNA at the intercalation site.     

Imino 1H- and 31P-NMR analysis of the interaction of propidium and ethidium with DNA

At low temperature and low salt concentration, both imino proton and ³¹p-nmr spectra of DNA complexes with the intercalators ethidium and propidium are in the slow-exchange region. Increasing temperature and/or increasing salt concentration results in an increase in the site exchange rate. Ring-current effects from the intercalated phenanthridinium ring of ethidium and propidium cause upfield shifts of the imino protons of A · T and G · C base pairs, which are quite similar for the two intercalators. The limiting induced chemical shifts for propidium and ethidium at saturation of DNA binding sites are approximately 0.9 ppm for A · T and 1.1 ppm for G · C base pairs. The similarity of the shifts for ethidium and propidium, in both the slow- and fast-exchange regions over the entire titration of DNA, shows that a binding model for propidium with neighbor-exclusion binding and negative ligand cooperativity is correct. The fact that a unique chemical shift is obtained for imino protons at intercalated sites over the entire titration and that no unshifted imino proton peaks remain at saturation binding of ethidium and propidium supports a neighbor-exclusion binding model with intercalators bound at alternating sites rather than in clusters on the double helix. Addition of ethidium and propidium to DNA results in downfield shifts in ³¹P-nmr spectra. At saturation ratios of intercalator to DNA base pairs in the titration, a downfield shoulder (approximately ?2.7 ppm) is apparent, which accounts for approximately 15% of the spectral area. The main peak is at ?3.9 to ?4.0 ppm relative to ?4.35 in uncomplexed DNA. The simplest neighbor-binding model predicts a downfield peak with approximately 50% of the spectral area and an upfield peak, near the chemical shift for uncomplexed DNA, with 50% of the area. This is definitely not the case with these intercalators. The observed chemical shifts and areas for the DNA complexes can be explained by models, for example, that involve spreading the intercalation-induced unwinding of the double helix over several base pairs and/or a DNA sequence- and conformation-dependent heterogeneity in intercalation-induced chemical shifts and resulting exchange rates.     

Hybrid simulation approach incorporating microscopic interaction along with rigid body degrees of freedom for stacking between base pairs

Stacking interaction between the aromatic heterocyclic bases plays an important role in the double helical structures of nucleic acids. Considering the base as rigid body, there are total of 18 degrees of freedom of a dinucleotide step. Some of these parameters show sequence preferences, indicating that the detailed atomic interactions are important in the stacking. Large variants of non‐canonical base pairs have been seen in the crystallographic structures of RNA. However, their stacking preferences are not thoroughly deciphered yet from experimental results. The current theoretical approaches use either the rigid body degrees of freedom where the atomic information are lost or computationally expensive all atom simulations. We have used a hybrid simulation approach incorporating Monte‐Carlo Metropolis sampling in the hyperspace of 18 stacking parameters where the interaction energies using AMBER‐parm99bsc0 and CHARMM‐36 force‐fields were calculated from atomic positions. We have also performed stacking energy calculations for structures from Monte‐Carlo ensemble by Dispersion corrected density functional theory. The available experimental data with Watson–Crick base pairs are compared to establish the validity of the method. Stacking interaction involving A:U and G:C base pairs with non‐canonical G:U base pairs also were calculated and showed that these structures were also sequence dependent. This approach could be useful to generate multiscale modeling of nucleic acids in terms of coarse‐grained parameters where the atomic interactions are preserved. This method would also be useful to predict structure and dynamics of different base pair steps containing non Watson–Crick base pairs, as found often in the non‐coding RNA structures. © 2015 Wiley Periodicals, Inc. Biopolymers 105: 212–226, 2016.     

Specificity in the genetic code. The role of nucleotide base-amino acid interaction

The mechanism of the unique and specific association of a given amino acid to its t-RNA is investigated by theoretical methods. Several possible schemes are proposed to explain specificity. The physical forces which act within these mechanisms are illustrated by the computer simulation of probable interactions between glycine and nucleotide bases and base pairs. It is demonstrated that glycine has direct and selective affinities for the nucleotide bases and that these interactions are principally determined by the polar groups. Energies have been calculated for the interaction of glycine with several base pairs. From these, the possibility that specificity arises through direct complexing of an amino acid with its anticodon is evaluated.     

Kinetic and equilibrium analysis of a threading intercalation mode: DNA sequence and ion effects

The interaction of a symmetric naphthalene diimide with alkylamino substituents at each imide position was investigated with the alternating sequence polymers, poly[d(A-T)]2 and poly[d(G-C)]2. Spectrophotometric binding studies indicate strong binding of the diimide to both sequences although the GC binding constant is 20-25 times larger than the AT binding constant. Analysis of the effects of salt concentration on the binding equilibria shows that the diimide forms two ion pairs in its complex with both polymers as expected for a simple dication. Stopped-flow kinetics experiments demonstrate that the diimide both associates and dissociates from DNA more slowly than classical intercalators with similar binding constants. Analysis of salt concentration effects on dissociation kinetics rate constants (kd) reveals that slopes in log kd versus log [Na+] plots are only approximately half the value obtained for classical dicationic intercalators that have both charged groups in the same groove. These kinetics results support a threading intercalation model, with one charged diimide substituent in each of the DNA grooves rather than with both side chains in the same groove, for the diimide complex with DNA. In the rate-determining step of the mechanism for dissociation of a threading complex only one ion pair is broken; the free side chain can then slide between base pairs to put both diimide side chains in the same groove, and this is followed by rapid full dissociation of the diimide. This sequential release of ion pairs makes the dissociation slope for dicationic threading intercalators more similar to the slope for classical monocationic intercalating ligands.(ABSTRACT TRUNCATED AT 250 WORDS)     

Post Hartree-Fock studies of the canonical Watson-Crick DNA base pairs: molecular structure and the nature of stability

Gas-phase gradient optimization was carried out on the canonical Watson-Crick DNA base pairs using the second-order M?ller-Plesset perturbation method at the 6-31G(d) and 6-31G(d,p) basis sets. It is detected that full geometry optimization at the MP2 level leads to an intrinsically nonplanar propeller-twisted and buckled geometry of G-C and A-T base pairs; while HF and DFT methods predict perfect planar or almost planar geometry of the base pairs. Supposedly the nonplanarity of the pairs is caused by pyramidalization of the amino nitrogen atoms, which is underestimated by the HF and DFT methods. This justifies the importance of geometry optimization at the MP2 level for obtaining reliable prediction of the charge distribution, molecular dipole moments and geometrical structure of the base pairs. The Morokuma-Kitaura and the Reduced Variational Space methods of the decomposition for molecular HF interaction energies were used for investigation of the hydrogen bonding in the Watson-Crick base pairs. It is shown that the HF stability of the hydrogen-bonded DNA base pairs originates mainly from electrostatic interactions. At the same time, the calculated magnitude of the second order intramolecular correlation correction to the Coulomb energy showed that electron correlation reduces the contribution of the electrostatic term to the attractive interaction for the A-T and G-C base pairs. Polarization, charge transfer and dispersion interactions also make considerable contribution to the attraction energy of bases.     

Model for the porphyrin-DNA binding site: ENDOR investigations of Cu-porphyrins binding to DNA.

Proton ENDOR has been observed from frozen solutions (ca. 38K degrees) of copper meso-(4-N-tetra-methylpyridyl)porphyrin (CuTMpyP(4)) complexed with Salmon sperm DNA in water and D2O. Lines from exchangeable protons of the DNA bases have been observed in these ENDOR spectra. Analyses of these ENDOR data show that the separations of these DNA protons from the copper atom are between 3.76 and 3.84 A with angles of 19.5 to 22.5 degrees between the Cu-H vectors and the gz axis. A distant ENDOR response has also been observed from phosphorous nuclei in the DNA backbone. We estimate that the phosphorous atoms producing this ENDOR signal are 7.5-10 A from the copper center of the porphyrin. These ENDOR data combined with results from an earlier NMR investigation have been used to construct a computer simulated model of the binding site in which the porphyrin is partially intercalated and extends into the major groove of DNA. The two GC base pairs at this site are slightly inequivalent. For each, the G imino proton and one of the C amino protons are at appropriate positions to account for the ENDOR signals arising from exchangeable protons. It is unlikely that this inequivalence would persist at room temperature where dynamic processes would give an apparently symmetric interaction. Although the model accounts for all reported experimental data involving tetracationic porphyrin species which have been suggested to be intercalators, it is not a unique solution.     

Binding of a macrocyclic bisacridine and ametantrone to CGTACG involves similar unusual intercalation platforms

The binding of a macrocyclic bisacridine and an antitumor intercalator ametantrone to DNA has been studied. We carried out X-ray diffraction analyses of the complexes between both intercalators and CGTACG. We have determined the crystal structure, by the multiple-wavelength anomalous diffraction (MAD) method, of bisacridine complexed with CGTA[br(5)C]G at 1.8 A resolution. The refined native crystal structure at 1.1 A resolution (space group C222, a = 29.58 A, b = 54.04 A, c = 40.22 A, and R-factor = 0.163) revealed that only one acridine of the bisacridine drug binds at the C5pG6 step of the DNA, with the other acridine plus both linkers completely disordered. Surprisingly, both terminal G.C base pairs are unraveled. The C1 nucleotide is disordered, and the G2 base is bridged to its own phosphate P2 through a hydrated Co(2+) ion. G12 is swung toward the minor groove with its base stacked over the backbone. The C7 nucleotide is flipped away from the duplex part and base paired to a 2-fold symmetry-related G6. The central four base pairs adopt the B-DNA conformation. An unusual intercalator platform is formed by bringing four complexes together (involving the 222 symmetry) such that the intercalator cavity is flanked by two sets of G x C base pairs (i.e., C5 x G8 and G6 x C7) on each side, joined together by G6 x G8 tertiary base pairing interactions. In the bisacridine-CGTACG complex, the intercalation platform is intercalated with two acridines, whereas in the ametantrone-CGTACG complex, only one ametantrone is bound. NMR titration of the bisacridine to AACGATCGTT suggests that the bisacridine prefers to bridge more than one DNA duplex by intercalating each acridine to different duplexes. The results may be relevant in understanding binding of certain intercalators to DNA structure associated with the quadruplet helix and Holliday junction.     

Half‐Intercalation Stabilizes Slipped Mispairing and Explains Genome Vulnerability to Frameshift Mutagenesis by Endogenous “Molecular Bookmarks”

Some 60 years ago chemicals that intercalate between base pairs of duplex DNA were found to amplify frameshift mutagenesis. Surprisingly, the robust induction of frameshifts by intercalators still lacks a mechanistic model, leaving this classic phenomenon annoyingly intractable. A promising idea of asymmetric half‐intercalation‐stabilizing frameshift intermediates during DNA synthesis has never been developed into a model. Instead, researchers of frameshift mutagenesis embraced the powerful slipped‐mispairing concept that unexpectedly struggled with the role of intercalators in frameshifting. It is proposed that the slipped mispairing and the half‐intercalation ideas are two sides of the same coin. Further, existing findings are reviewed to test predictions of the combined “half‐intercalation into the slipped‐mispairing intermediate” model against accumulated knowledge. The existence of potential endogenous intercalators and the phenomenon of “DNA bookmarks” reveal ample possibilities for natural frameshift mutagenisis in the cell. From this alarming perspective, it is discussed how the cell could prevent genome deterioration from frameshift mutagenesis.     

Post Hartree-Fock Studies of the Canonical Watson-Crick DNA Base Pairs: Molecular Structure and the Nature of Stability

Abstract

Gas-phase gradient optimization was carried out on the canonical Watson-Crick DNA base pairs using the second-order Møller-Plesset perturbation method at the 6–31G(d) and 6- 31G(d,p) basis sets. It is detected that full geometry optimization at the MP2 level leads to an intrinsically nonplanar propeller-twisted and buckled geometry of G-C and A-T base pairs; while HF and DFT methods predict perfect planar or almost planar geometry of the base pairs. Supposedly the nonplanarity of the pairs is caused by pyramidalization of the amino nitrogen atoms, which is underestimated by the HF and DFT methods. This justifies the importance of geometry optimization at the MP2 level for obtaining reliable prediction of the charge distribution, molecular dipole moments and geometrical structure of the base pairs. The Morokuma-Kitaura and the Reduced Variational Space methods of the decomposition for molecular HF interaction energies were used for investigation of the hydrogen bonding in the Watson-Crick base pairs. It is shown that the HF stability of the hydrogen-bonded DNA base pairs originates mainly from electrostatic interactions. At the same time, the calculated magnitude of the second order intramolecular correlation correction to the Coulomb energy showed that electron correlation reduces the contribution of the electrostatic term to the attractive interaction for the A-T and G-C base pairs. Polarization, charge transfer and dispersion interactions also make considerable contribution to the attraction energy of bases.     

Theoretical analysis of noncanonical base pairing interactions in RNA molecules

Noncanonical base pairs in RNA have strong structural and functional implications but are currently not considered for secondary structure predictions. We present results of comparative ab initio studies of stabilities and interaction energies for the three standard and 24 selected unusual RNA base pairs reported in the literature. Hydrogen added models of isolated base pairs, with heavy atoms frozen in their 'away from equilibrium' geometries, built from coordinates extracted from NDB, were geometry optimized using HF/6-31G** basis set, both before and after unfreezing the heavy atoms. Interaction energies, including BSSE and deformation energy corrections, were calculated, compared with respective single point MP2 energies, and correlated with occurrence frequencies and with types and geometries of hydrogen bonding interactions. Systems having two or more N-H...O/N hydrogen bonds had reasonable interaction energies which correlated well with respective occurrence frequencies and highlighted the possibility of some of them playing important roles in improved secondary structure prediction methods. Several of the remaining base pairs with one N-H...O/N and/or one C-H...O/N interactions respectively, had poor interaction energies and negligible occurrences. High geometry variations on optimization of some of these were suggestive of their conformational switch like characteristics.     

Model for the Porphyrin-DNA Binding Site: ENDOR Investigations of Cu-Porphyrins Binding to DNA

Abstract

Proton ENDOR has been observed from frozen solutions (ca. 38K°) of copper meso-(4-N-tetra-methylpyridyl)porphyrin (CuTMpyP(4)) complexed with Salmon sperm DNA in water and D₂O. Lines from exchangeable protons of the DNA bases have been observed in these ENDOR spectra. Analyses of these ENDOR data show that the separations of these DNA protons from the copper atom are between 3.76 and 3.84 A with angles of 19.5 to 22.5 degrees between the Cu-H vectors and the g_z axis. A distant ENDOR response has also been observed from phosphorous nuclei in the DNA backbone. We estimate that the phosphorous atoms producing this ENDOR signal are 7.5–10 Å from the copper center of the porphyrin. These ENDOR data combined with results from an earlier NMR investigation (1) have been used to construct a computer simulated model of the binding site in which the porphyrin is partially intercalated and extends into the major groove of DNA. The two GC base pairs at this site are slightly inequivalent. For each, the G imino proton and one of the C amino protons are at appropriate positions to account for the ENDOR signals arising from exchangeable protons. It is unlikely that this inequivalence would persist at room temperature where dynamic processes would give an apparently symmetric interaction. Although the model accounts for all reported experimental data involving tetracationic porphyrin species which have been suggested to be intercalators, it is not a unique solution.     

Theoretical studies on protein-nucleic acid interactions. II. Hydrogen bonding of amino acid side chains with bases and base pairs of nucleic acids

With a view to understanding the role of hydrogen bonds in the recognition of nucleic acids by proteins, hydrogen bonding between the bases and base pairs of nucleic acids and the amino acids (Asn, Gln, Asp and Glu, and charged residues Arg⁺, Glu^?, and Asp^?) has been studied by a second-order perturbation theory. Binding energies have been calculated for all possible configurations involving a pair of hydrogen bonds between the base (or base pair) and the amino acid residue. Our results show that the hydrogen bonding in these cases has a large contribution from electrostatic interaction. In general, the charged amino acids, compared to the uncharged ones, form more stable complexes with bases or base pairs. The hydrogen-bond energies are an order of magnitude smaller than the Coulombic interaction energies between basic amino acids (Lys⁺, Arg⁺, and His⁺) and the phosphate groups of nucleic acids. The stabilities of the complexes of amino acids Asn, Gln, Asp, and Glu with bases are in the order: G–X > C–X > A–X U–X or T–X, and G · C–X > A · T(U)–X, where X is one of these amino acid residues. It has been shown that Glu^? and Asp^? can recognize guanine in single-stranded nucleic acids; Arg⁺ can recognize G · C base pairs from A · T base pairs in double-stranded structures.     

Theoretical analysis of noncanonical base pairing interactions in RNA molecules

Noncanonical base pairs in RNA have strong structural and functional implications but are currently not considered for secondary structure predictions. We present results of comparative ab initio studies of stabilities and interaction energies for the three standard and 24 selected unusual RNA base pairs reported in the literature. Hydrogen added models of isolated base pairs, with heavy atoms frozen in their ‘away from equilibrium’ geometries, built from coordinates extracted from NDB, were geometry optimized using HF/6-31G** basis set, both before and after unfreezing the heavy atoms. Interaction energies, including BSSE and deformation energy corrections, were calculated, compared with respective single point MP2 energies, and correlated with occurrence frequencies and with types and geometries of hydrogen bonding interactions. Systems having two or more N-H...O/N hydrogen bonds had reasonable interaction energies which correlated well with respective occurrence frequencies and highlighted the possibility of some of them playing important roles in improved secondary structure prediction methods. Several of the remaining base pairs with one N-H...O/N and/or one C-H...O/N interactions respectively, had poor interaction energies and negligible occurrences. High geometry variations on optimization of some of these were suggestive of their conformational switch like characteristics.     

Theoretical studies on protein-nucleic acid interactions. III. Stacking of aromatic amino acids with bases and base pairs of nucleic acids

Stacking of aromatic amino acids tryptophan (Trp), tyrosine (Tyr), phenylalanine (Phe), and histidine (His) with bases and base pairs of nucleic acids has been studied. Stacking energies of the amino acid–base (or base pair) complexes have been calculated by second-order perturbation theory. Our results show that, in general, the predominant contribution to the total stacking energy comes from the dispersion terms. In these cases, repulsion energy is greater than the sum of electrostatic and polarization energies. In contrast to this, interaction of histidine with the bases and base pairs is largely Coulombic in nature. The complexes of guanine with aromatic amino acids are more stable than the corresponding complexes of adenine. Among pyrimidines, cytosine forms the most stable complexes with the aromatic amino acids. The G · C base pair has the highest affinity with aromatic amino acids among various sets of base pairs. Optimized geometries of the stacked complexes show that the aromatic moieties overlap only partially. The heteroatom of one residue generally overlaps with the other aromatic moiety. There is a considerable degree of configurational freedom in the stacked geometries. The role of stacking in specific recognition of base sequences by proteins is discussed.     

Electrostatic factors in DNA intercalation

The factors that determine the binding of a chromophore between the base pairs in DNA intercalation complexes are dissected. The electrostatic potential in the intercalation plane is calculated using an accurate ab initio based distributed multipole electrostatic model for a range of intercalation sites, involving different sequences of base pairs and relative twist angles. There will be a significant electrostatic contribution to the binding energy for chromophores with a predominantly positive electrostatic potential, but this varies significantly with sequence, and somewhat with twist angle. The usefulness of these potential maps for understanding the binding of intercalators is explored by calculating the electrostatic binding energy for 9-aminoacridine, ethidium, and daunomycin in a variety of model binding sites. The electrostatic forces play a major role in the positioning of an intercalating 9-aminoacridine and a significant stabilizing role in the binding of ethidium in its sterically constrained position, but the intercalation of daunomycin is determined by the side-chain binding. Sequence preferences are likely to be determined by a complex and subtle mixture of effects, with electrostatics being just one component. The electrostatic binding energy is also unlikely to be a major determinant of the twist angle, as its variation with angle is modest for most intercalation sites. Overall, the electrostatic potential maps give guidance on how positively charged chromophores can be chemically adapted by heteroatomic substitution to optimise their binding.     

NMR study of the drug-base overlap geometry in the dark complex of 8-methoxypsoralen and d(pApT)4

The dark binding of 8-methoxypsoralen (MOP) to d(pApT)4 was investigated by 270-MHz 1H nuclear magnetic resonance (NMR) spectra. The continuous high-field shifts of the MOP resonances by d(pApT)4 at low temperatures indicate fast exchange between free and bound drug. The limiting complexation shifts of the various MOP protons between 0.36 (CH3) and 1.20 ppm (H5) are in the range expected for an intercalation complex. The NMR line widths of the MOP ring protons vary with the square of the observed complexation shifts (maximum at H5), indicating a dominant effect of the fast exchange between free and bound drug. The corresponding kinetic parameters agree with the values previously reported for a variety of other intercalators. The observed exchange broadenings were also used as a criterion to limit the uncertainty connected with fast averaging of the signals of the drug in potential multiple binding modes: A qualitatively different pattern of broadenings (minimum at H5) is expected from fast exchange between the two binding modes related by the short 2-fold quasi-symmetry axis of MOP. The measured complexation shifts were compared to theoretical values calculated on the basis of coplanar intercalation with base pair arrangements derived from typical published intercalation site geometries. The standard deviation between observed and calculated shifts was considerably smaller for asymmetrical intercalation between the bases of the same strand (less than or equal to 0.11 ppm) than for symmetrical intercalation between the base pairs (greater than or equal to 0.28 ppm).(ABSTRACT TRUNCATED AT 250 WORDS)     

Selective catalysis of A.T base pair proton exchange in DNA complexes: imino proton NMR analysis.

Polyaromatic molecules with amino chain substituents, upon binding with DNA, selectively catalyze exchange of the A.T base pair protons with bulk water protons. The amine-catalyzed exchange is mediated by compounds which are A.T and G.C base sequence specific, intercalators, and outside binders. A mechanism for the selective exchange, involving transient opening and closing of individual A.T base pairs in the duplex, is discussed.