About this Journal Submit a Manuscript Table of Contents
Journal of Theoretical Chemistry
Volume 2013 (2013), Article ID 526569, 8 pages
http://dx.doi.org/10.1155/2013/526569
Research Article

DFT Description of Intermolecular Forces between 9-Aminoacridines and DNA Base Pairs

1Departamento de Química y Biología, Universidad del Norte, Km 5 Antigua vía a Puerto Colombia, Barranquilla, Colombia
2Grupo de Investigación Max Planck, Universidad del Atlántico, Km 7 Antigua vía a Puerto Colombia, Barranquilla, Colombia
3Universidad del Atlántico, Km 7 Antigua vía a Puerto Colombia, Barranquilla, Colombia

Received 15 March 2013; Revised 29 May 2013; Accepted 16 July 2013

Academic Editors: A. M. Lamsabhi, G. Monard, and H.-Y. Zhang

Copyright © 2013 Sandra Cotes Oyaga et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

The B3LYP method with 6-31G* basis set was used to predict the geometries of five 9-aminoacridines (9-AA 1(a–e)), DNA base pairs, and respective complexes. Polarizabilities, charge distribution, frontier molecular orbital (FMO), and dipole moments were used to analyze the nature of interactions that allow reasonable drug diffusion levels. The results showed that charge delocalization, high polarizabilities, and high dipole moments play an important role in intermolecular interactions with DNA. The interactions of 9-AA 1(a–e) with GC are the strongest. 9-AA 1(d) displayed the strongest interaction and 9-AA 1(b) the weakest.

1. Introduction

Intermolecular recognition is a key process for the interaction with biological systems. The strength of intermolecular weak interactions such as hydrogen bonding, Van der Waals forces, aromatic stacking, binding, and intercalation is closely related to biological activity. Acridines are known as antitumor [15], antiviral [68], antiprion [9, 10], antimicrobial [11], anti-inflammatory, and analgesic [12].

9-Aminoacridines (9-AA) have also been considered for the treatment of protozoal infections [13], and anticancer activity was first considered in the 1940s [14]. Since then, a large number of 9-aminoacridine drugs, natural alkaloids, or synthetic molecules have been tested as antitumoral agents. 9-AA activity is due to the intercalation of the tricyclic aromatic ring between adjacent base-pairs [15]. Thereby high affinity for DNA is the result of stacking between base-pairs and enhanced by electrostatic and hydrophobic interactions with other substituents. Strong binding correlates cytotoxicity; however, extravascular distribution is limited by low levels of circulating free drug. Thus biological activity is the result of a compromise between binding and circulating drug levels.

Sebestik et al. [16, 17] reported the synthesis of acridine conjugates 9-AA 1(a–e) possessing amino acids bonded to an ethylenediamine linker that showed low affinity to nucleic acids. Sebestik et al. reduced the DNA-binding constant of 9-aminoacridines by introduction of DNA binders with different base-pair affinity. 9-AA as a GC-DNA binder and peptide residues as AT-DNA binders were connected to the ethylenediamine linker. The peptide conjugates 9-AA 1(a–e) showed weak interactions with DNA due to the competitive role of substituents.

DFT methods are accurate methods for computational studies [18, 19], and electrostatic interactions are described by available force fields [2022]. However, the quest for DFT methods which accurately describe long-range electron correlation continues to be a challenge in computational chemistry. Modeling of DNA interactions with 9-AA offers an insight about the physical and chemical phenomena that are not measurable through experimentation. In this work, the structures and electrostatic interactions of 9-AA 1(a–e) with AT and GC (Figure 1) are calculated and analyzed. To our knowledge, this is the first report showing the predicted electrostatic interactions of these molecules with DNA base-pairs that might be responsible for the observed weak interactions. The effect of polarizability, dipole moment, FMO energies of the base-pair with 9-AA 1(a–e) was also studied by DFT. Further improvements in the quantum chemistry calculations of base-pair interactions with 9-AA 1(a–e) including intermolecular electron correlation effects remain to be done as a continuing effort to estimate base-pairing energies sufficiently accurate for most applications.

526569.fig.001
Figure 1: Substituted 9-aminoacridines 1(a–e).

2. Computational Methods

Calculations on the isolated molecules and molecular complexes were performed using Spartan ’08 [18] molecular mechanic MMFF method. Geometries optimizations were performed at the B3LYP level using 6-31G* basis set. All geometries correspond to potential energy minima as shown by harmonic normal-mode analysis at the same level of theory. Molecular descriptors such as dipole moment, polarizability, and FMO calculations have been carried out using the B3LYP/6-31G** optimized geometries within Gaussian 09W [23].

3. Results and Discussion

Numbering of the DNA base-pairs is shown in Figure 2. Optimized structures of the molecular complexes between 9-AA 1(a–e) and base-pairs are shown in Figures 3 and 4. Equilibrium geometries energetics are shown in Table 1. From equilibrium geometries shown in Figures 3 and 4, it is seen that coplanarity is not retained in the complex of 9-AA 1(c) with AT and complexes with GC.

tab1
Table 1: Energies of optimized 9-AA 1(a–e)—base-pair complexes calculated at the B3LYP/6-31G* basis set.
526569.fig.002
Figure 2: Base-pair numbering AT and GC.
526569.fig.003
Figure 3: Optimized structures of complexes 9-AA 1(a–e) AT calculated at the B3LYP/6-31G* basis set.
526569.fig.004
Figure 4: Optimized structures of complexes 9-AA 1(a–e) GC calculated at the B3LYP/6-31G* basis set.

Atomic charges, polarizability, dipole moment, and FMO are shown in Tables 24. Among the atoms of the heterocyclic ring in the isolated 9-AA 1(a–e), C13 and C6 have the maximum positive charges which are caused by the bonding to a negatively charged electronegative N14. Instead 9-AA 1(c) exhibits a maximum positive charge in C11, which is bonded to a negatively charged N15. See Table 3. In general, acridine atomic charges in 9-AA 1(a–e) are delocalized by resonance effects.

tab2
Table 2: AT, GC atomic charges* and complex with 9-AA 1(a–e).
tab3
Table 3: 9-AA 1(a–e)atomic charges* and respective complexes with AT and GC.
tab4
Table 4: HOMO and LUMO energies (ev), polarizabilities (au), and dipole moments (D) of 9-AA 1(a–e), base pairs and isolated bases, calculated at the B3LYP/6-31G** basis set.

Significant changes in atomic charges are observed in the carbon atoms adjacent to the acridine N14 in the complexes of 9-AA 1(a–e) with DNA base-pairs. See Table 3. However the highest changes are displayed in the complexes with GC. Note, for instance, N14 in 9-AA 1(d) changes from −0.611 to 0.101 by interacting with GC while with AT it changes to −0.127. See Table 3.

The study of the charges in the complexes of 9-AA 1(a–e) with base-pairs showed that the strongest interaction with GC is 9-AA 1(d) while the weakest is with 9-AA 1(b). Accordingly, the study of atom charges involved in hydrogen bonding in GC and respective complexes with 9-AA 1(a–e) denoted the highest charge changes in N5. Note also N5 in GC changes from −0.618 to −0.084 in the complex with 9-AA 1(d). See Table 2.

On the other hand, peptide chains in 9-AA 1(a–e) exhibit intramolecular hydrogen bonding. Significant atomic charges and torsion angles changes are also observed in the peptide chains after complex formation with base-pairs. Torsion angle changes in 9-AA 1(a) and (b) are observed from C25, causing a hindrance to intramolecular hydrogen-bonding formation. Oxygen atoms in the peptide chains of 9-AA 1(d) and (e) are oriented towards hydrogen in GC causing a hindrance in the formation of intramolecular hydrogen bonding as well.

9-AA 1(a–e) and DNA base-pairs are electron-deficient aromatics, and the quadrupole moment is reversed due to the electronegativity of N14 and N15 in 9-AA and electron-withdrawing functional groups in base-pairs; see Figure 5. Face-centered stacking is disfavored in the complexes of 9-AA with base-pairs by repulsive quadrupole/quadrupole interactions. By offsetting one of the rings, a displaced configuration reduces the repulsive interactions and stabilizes complexes configuration enhancing electrostatic interactions and hydrogen bonding.

526569.fig.005
Figure 5: Electrostatic potential maps of 9-AA and base-pairs AT and GC, respectively.

In particular 9-AA 1(a, b, and e) exhibit and electrostatic interaction between positively charged hydrogen and N17 in their complexes with AT, giving rise to in plane complexes. In the complex of 9-AA 1(c) and AT, repulsive interaction is overcome by hydrogen bonding between AT H30 and the carbonyl of the first Ile residue of the 9-AA. Also a hydrogen bonding between H32 of 9-AA 1(c) with second Ile carbonyl oxygen and another hydrogen bonding between oxygen of first Ile residue with H30 in AT and with H16 in 9-AA 1(c) allows the stabilization of the calculated configuration.

The complex of 9-AA 1(d) with AT has shows an inverted position of the base-pair, allowing direct interaction between hydrogen joined to the most negatively charged N21 of AT and N14 of 9-AA 1(d). Also both Ile carbonyl residues are directed towards the positively charged 9-AA aromatic rings.

The series of complexes with GC has shown orientation of 9-AA carbonyl oxygen of residues towards base-pair core rings as can be seen in Figure 4. Hydrogen bonding between N14 of 9-AA 1(b) and H14 in GC stabilizes the configuration. The complex between 9-AA 1(c) and GC shows 3 hydrogen bonding: intramolecular hydrogen bonding between a carbonyl oxygen of the Ile residue and H16 of 9-AA 1(c), hydrogen bonding between a carbonyl oxygen of the Ile residue and H28 in GC, and hydrogen bonding between 9-AA 1(c) amino terminal group and H14 in GC.

Hydrogen bonding between carbonyl oxygen of Gly residue in 9-AA 1(d) and H14 in GC stabilizes the observed configuration in Figure 4. 9-AA 1(e) shows orientation of the three carbonyl oxygen atoms of the 9-AA residues towards the base-pair ring core as stated before. This complex also has a hydrogen bonding N14 of 9-AA 1(e) and H2 of GC.

Interactions of 9-AA 1(a–e) with base-pairs are not based upon inherent attraction between pi cloud electron density, so pi stacking is not appropriate to describe them.

C–C bonds in 9-AA 1(a–e) rings showed characteristic sp2 bond distances of polycyclic aromatic systems. 9-AA 1(a–e) are rigid rings so bond distances changes are not significant by interacting with bases pairs. Base-pair bond lengths were kept almost invariable after complex interactions with 9-AA 1(a–e).

9-AA 1(a–e) have shown an asymmetric distribution of atomic charges in electrostatic potential maps; this is supported by the high dipole moments calculated; see Table 4. Symmetric distribution of atomic charges of the polycyclic aromatic rings in 9-AA 1(a–e) was observed due to resonance effects.

The sequence of the intermolecular interaction energies is the same on all base-pairs: 9-AA 1(a) < 9-AA 1(b) < 9-AA 1(d) < 9-AA 1(c) < 9-AA 1(e) in Table 5. From the interaction energies it is seen that 9-AA 1(c) and 9-AA 1(e) exhibit the strongest interactions with base-pairs. These results are supported on the number of hydrogen bonding allowed for the adopted conformation in complexes of 9-AA 1(c) and 9-AA 1(e).

tab5
Table 5: Interaction energies of optimized 9-AA 1(a–e)—base-pair complexes calculated at the B3LYP/6-31G* basis set.

The lowest energy difference between FMO was shown by 9-AA 1(d), (0,127 eV), in Table 4. 9-AA 1(a–e) showed negative energy values for LUMO and HOMO, indicating that they are good electron acceptors. LUMO energies of isolated bases and base-pairs are more positive than LUMO energies in 9-AA 1(a–e), so isolated bases and base-pairs possess electron-donating abilities, which are supported by the high ionization potential [24, 25]. Among the isolated bases Guanine is the best nucleophile having the highest LUMO energy value. Thus the electron donor-electron acceptor contribution could play a role in the interacting systems.

The highest dipole moment in the 9-AA 1(a–e) series was found in 9-AA 1(c), favoring electrostatic attractions with base-pairs. As long as polarizability is in isolated bases, Guanine also exhibits the highest value (120,114 au) and GC is the most polarizable base-pair (228,859 au); see Table 4. 9-AA 1(a–e) have polarizabilities from 350.554 to 202.861 au, 9-AA being 1(c) the largest.

Finally, high polarizability favors interactions between the systems and also high dipole moment favors electrostatic interactions. From Table 4 it can be noted that 9-AA 1(e) has a large polarizability (298,449 au) but also has a lower dipole moment when compared to 9-AA 1(d) and 1(c). The strongest interactions of the complexes are then a result of mainly polarizability, hydrogen bonding, and dipole moments. The weakest interaction is found with low dipole moments and low polarizability. Charge transfer, on the other hand, could play a less significant role as there is no significant difference between FMO energies of the systems.

4. Conclusion

9-AA 1(a–e) have been shown to have high polarizabilities, high dipole moments, and electron acceptors abilities; on the opposite side, base-pairs are electron donors, resulting in favorable interactions between the systems under examination in this paper. From the 9-AA under consideration, the strongest interaction was exhibited by 9-AA 1(c) and 9-AA 1(e). The weakest interaction was observed with 9-AA 1(b).

Complex interactions are stabilized by electrostatic and hydrogen-bonding contributions while charge transfer is not significant for the type of molecules under study. Theoretical procedures that cover the type of interactions described before are therefore a first good approximation towards the study of 9-aminoacridines.

Acknowledgments

The authors greatly appreciate the financial support of the Dirección de Investigación, Desarrollo e Innovación, and DIDI at the Universidad del Norte.

References

  1. P. Belmont, J. Bosson, T. Godet, and M. Tiano, “Acridine and acridone derivatives, anticancer properties and synthetic methods: where are we now?” Anti-Cancer Agents in Medicinal Chemistry, vol. 7, no. 2, pp. 139–169, 2007. View at Publisher · View at Google Scholar · View at Scopus
  2. F. Charmantray, M. Demeunynck, D. Carrez et al., “4-Hydroxymethyl-3-aminoacridine derivatives as a new family of anticancer agents,” Journal of Medicinal Chemistry, vol. 46, no. 6, pp. 967–977, 2003. View at Publisher · View at Google Scholar · View at Scopus
  3. R. A. Heald and M. F. G. Stevens, “Antitumour polycyclic acridines. Palladium(0) mediated syntheses of quino[4,3,2-kl]acridines bearing peripheral substituents as potential telomere maintenance inhibitors,” Organic & Biomolecular Chemistry, vol. 19, pp. 3377–3389, 2003. View at Publisher · View at Google Scholar
  4. R. Hegde, P. Thimmaiah, M. C. Yerigeri, G. Krishnegowda, K. N. Thimmaiah, and P. J. Houghton, “Anti-calmodulin acridone derivatives modulate vinblastine resistance in multidrug resistant (MDR) cancer cells,” European Journal of Medicinal Chemistry, vol. 39, no. 2, pp. 161–177, 2004. View at Publisher · View at Google Scholar · View at Scopus
  5. M. Kimura, A. Kato, and I. Okabayashi, “Acridine derivatives. IV. Synthesis, molecular structure, and antitumor activity of the novel 9-anilino-2,3-methylenedioxyacridines,” Journal of Heterocyclic Chemistry, vol. 29, no. 1, pp. 73–80, 1992. View at Scopus
  6. M. Fujiwara, M. Okamoto, M. Okamoto et al., “Acridone derivatives are selective inhibitors of HIV-1 replication in chronically infected cells,” Antiviral Research, vol. 43, no. 3, pp. 189–199, 1999. View at Scopus
  7. N. W. Luedtke, Q. Liu, and Y. Tor, “RNA—ligand interactions: affinity and specificity of aminoglycoside dimers and acridine conjugates to the HIV-1 Rev response element,” Biochemistry, vol. 42, no. 39, pp. 11391–11403, 2003. View at Publisher · View at Google Scholar · View at Scopus
  8. I. B. Taraporewala, “Thiazolo[5,4-b]acridines and thiazolo[4,5-b]acridines: probable pharmacophores of antiviral and anti-tumor marine alkaloids,” Tetrahedron Letters, vol. 32, no. 1, pp. 39–42, 1991. View at Publisher · View at Google Scholar
  9. C. Korth, B. C. May, F. E. Cohen, and S. B. Prusiner, “Acridine and phenothiazine derivatives as pharmacotherapeutics for prion disease,” Proceedings of the National Academy of Sciences of the United States of America, vol. 98, no. 17, pp. 9836–9841, 2001.
  10. B. C. May, A. T. Fafarman, S. B. Hong et al., “Potent inhibition of scrapie prion replication in cultured cells by bis-acridines,” Proceedings of the National Academy of Sciences of the United States of America, vol. 100, no. 6, pp. 3416–3421, 2003.
  11. H. P. Kavitha, “Synthesis and antimicrobial activity of 1-(9′-acridinyl)-5-(4-substituted phenyl) tetrazoles,” Asian Journal of Chemistry, vol. 16, no. 2, pp. 1191–1193, 2004. View at Scopus
  12. S. M. Sondhi, M. Johar, N. Singh, R. Shukla, R. Raghubir, and S. G. Dastidar, “Synthesis of sulpha drug acridine derivatives and their evaluation for anti-inflammatory, analgesic and anticancer activity,” Indian Journal of Chemistry B, vol. 41, no. 12, pp. 2659–2666, 2002. View at Scopus
  13. A. Albert, Ricard Clay and Company, Ltda, 1951.
  14. I. Lasnitzki and J. H. Wilkinson, “The effect of acridine derivatives on growth and mitoses of cells in vitro,” British Journal of Cancer, vol. 2, pp. 369–375, 1948. View at Publisher · View at Google Scholar
  15. J. B. Le Pecq, M. Le Bret, J. Barbet, and B. Roques, “DNA polyintercalating drugs: DNA binding of diacridine derivatives,” Proceedings of the National Academy of Sciences of the United States of America, vol. 72, no. 8, pp. 2915–2919, 1975. View at Scopus
  16. J. Sebestik, I. Stibor, and J. Hlavcek, “New peptide conjugates with 9-aminoacridine: synthesis and binding to DNA,” Journal of Peptide Science, vol. 12, no. 7, pp. 472–480, 2006. View at Publisher · View at Google Scholar
  17. C. B. Carlson and P. A. Beal, “Point of sttachment and sequence of immobilized peptide-acridine conjugates control affinity for nucleic acids,” Journal of the American Chemical Society, vol. 124, no. 29, pp. 8510–8511, 2002. View at Publisher · View at Google Scholar
  18. J. Sponer, J. Leszczynski, and O. Hobza, “Nature of nucleic acid-base stacking: nonempirical ab initio and empirical potential characterization of 10 stacked base dimers. Comparison of stacked and H-bonded base pairs,” The Journal of Physical Chemistry A, vol. 100, no. 13, pp. 5590–5596, 1996. View at Publisher · View at Google Scholar
  19. P. Hobza, M. Kabelac, J. Sponer, P. Mejzlik, and J. Vondrasek, “Performance of empirical potentials (AMBER, CFF95, CVFF, CHARMM, OPLS, POLTEV), semiempirical quantum chemical methods (AM1, MNDO/M, PM3), and ab initio Hartree–Fock method for interaction of DNA bases: comparison with nonempirical beyond Hartree–Fock results,” Journal of Computational Chemistry, vol. 18, no. 9, pp. 1136–1150, 1997.
  20. S. Riahi, M. R. Ganjali, A. B. Moghaddam, P. Norouzi, and M. Niasari, “Determination of the electrode potentials for substituted 1,2-dihydroxybenzenes in aqueous solution: theory and experiment,” Journal of Molecular Structure, vol. 774, no. 1–3, pp. 107–111, 2006. View at Publisher · View at Google Scholar · View at Scopus
  21. S. Riahi, M. R. Ganjali, A. B. Moghaddam, and P. Norouzi, “Theoretical and experimental study of electrical and electrochemical properties of (E)-3-(4,5-dihydroxy-2-(phenylsulphonyl) phenyl) acrylic acid as a new caffeic acid derivative,” Journal of Theoretical and Computational Chemistry, vol. 6, no. 2, pp. 255–268, 2007. View at Publisher · View at Google Scholar · View at Scopus
  22. S. Riahi, A. B. Moghaddam, M. R. Ganjali, and P. Norouzi, “Determination of the oxidation potentials of pyrogallol and some of its derivatives: theory and experiment,” Journal of Theoretical and Computational Chemistry, vol. 6, no. 2, pp. 331–340, 2007. View at Publisher · View at Google Scholar · View at Scopus
  23. M. J. Frisch, G. W. Trucks, H. B. Schlegel, et al., “Gaussian 09, Revision A02,” Gaussian, Inc, Wallingford, UK, 2009.
  24. F. Prat, K. N. Houk, and C. S. Foote, “Effect of Guanine stacking on the oxidation of 8-oxoguanine in B-DNA,” Journal of the American Chemical Society, vol. 120, no. 4, pp. 845–846, 1998. View at Publisher · View at Google Scholar
  25. A. O. Colson and M. D. Sevilla, “Elucidation of primary radiation damage in DNA through application of Ab initio molecular orbital theory,” International Journal of Radiation Biology, vol. 67, no. 6, pp. 627–645, 1995. View at Publisher · View at Google Scholar