Theoretical studies on interaction of anticancer drugs (dacarbazine,procarbazine and triethylenemelamine) with normal (AT and GC) and mismatch (GG,CC, AA and TT) base pairs

This study is an attempt to gain a better understanding of the physicochemical interaction between novel anticancer drugs and DNA bases. We have employed quantum chemical tools to explore the interaction of a few anticancer drugs [namely procarbazine (PR), dacarbazine (DC) and triethylenemelamine (TR)] with isolated normal (GC and AT) and mismatch (AA, CC, GG and TT) base pairs. The molecular geometries, electronic structural stability, vibrational energies, chemical reactivity and other electronic properties were studied using MP2/6-311+G^**, B3LYP/6-311+G^** and M05-2X/6-311+G^** methods. The optimised geometries of the usual and mismatch base pairs are almost planar whereas the geometries of drug-interacting complexes deviate from planarity. The presence of steric hindrance and π-bond overlaps between C–C bonds in the complexes has distorted the planarity of the four- and five-member rings in the base pairs. Among the three drugs chosen, DC and PR bond well with normal and mismatch base pairs with large interaction energy. The electron density (ED) difference maps of the most stable GG–DC, GG–PR and GG–TR drug-interacting complexes show the information about sharing of ED and gain or loss of ED within the interacting molecules. The stabilisation energy of the charge transfer interaction between the relevant donor–acceptor orbital of GG–DC and GC–DC complexes has been found to be around 16 kcal/mol and GG–PR and GC–PR complexes has been found to be around 12 kcal/mol. But, for the GG–TR and GC–TR complexes, the stabilisation energy is found to be less than 6 kcal/mol. Moreover, the topological analysis of hydrogen bond network of DC and PR drug-interacting complexes have high electron and Laplacian density with structural stability at the bond critical points (BCPs), while compared TR drug-interacting complexes by atoms in molecules and natural bond orbital analysis. Finally, we may conclude that the drugs DC and PR are highly efficient drugs to target normal and mismatch base pair for control and inhibition of DNA replication.     

Alteration of the groove width of DNA induced by the multimodal hydrogen bonding of denaturants with DNA bases in its grooves affects their stability

BackgroundDenaturants, namely, urea and guanidinium chloride (GdmCl) affect the stability as well as structure of DNA. Critical assessment of the role of hydrogen bonding of these denaturants with the different regions of DNA is essential in terms of its stability and structural aspect. However, the understanding of the mechanistic aspects of structural change of DNA induced by the denaturants is not yet well understood.MethodsIn this study, various spectroscopic along with molecular dynamics (MD) simulation techniques were employed to understand the role of hydrogen bonding of these denaturants with DNA bases in their stability and structural change.Results and conclusionIt has been found that both, GdmCl and urea intrude into groove region of DNA by striping surrounding water. The hydrogen bonding pattern of Gdm⁺ and urea with DNA bases in its groove region is multimodal and distinctly different from each other. The interaction of GdmCl with DNA is stabilized by electrostatic interaction whereas electrostatic and Lennard-Jones interactions both contribute for urea. Gdm⁺ forms direct hydrogen bond with the bases in the minor groove of DNA whereas direct and water assisted hydrogen bond takes place with urea. The hydrogen bond formed between Gdm⁺ with bases in the groove region of DNA is stronger than urea due to strong electrostatic interaction along with less self-aggregation of Gdm⁺ than urea. The distinct hydrogen bonding capability of Gdm⁺ and urea with DNA bases in its groove region affects its width differently. The interaction of Gdm⁺ decreases the width of the minor and major groove which probably increases the strength of hydrogen bond between the Watson-Crick base pairs of DNA leading to its stability. In contrast, the interaction of urea does not affect much to the width of the grooves except the marginal increase in the minor groove width which probably decreases the strength of hydrogen bond between Watson Crick base pairs leading to the destabilization of DNA.General significanceOur study clearly depicts the role of hydrogen bonding between DNA bases and denaturants in their stability and structural change which can be used further for designing of the guanidinium based drug molecules.     

Distamycin and Penta-N-Methylpyrrolecarboxamide Binding Sites on Native DNA A Comparison of Methidiumproyl-EDTA-Fe(II) Footprinting and DNA Affinity Cleaving

Abstract

Using two direct methods we have studied the binding locations and site sizes of distamycin and penta-N-methylpyrrolecarboxamide on three DNA restriction fragments from pBR322 plasmid. We find that methidiumpropyl-EDTA·Fe(II) footprinting and DNA affinity cleaving methods report common binding locations and site sizes for the tri- and pentapeptides bound to heterogeneous DNA. The tripeptide distamycin binds 5-base-pair sites with a preference for poly(dA)·poly(dT) regions. The pentapeptide binds 6–7-base-pair sites with a preference for poly(dA)·poly(dT) regions. These results are consistent with distamycin binding as an isogeometric helix to the minor groove of DNA with the four carboxamide N-H's hydrogen bonding five A+T base pairs. The data supports a model where each of the carboxamide N-H's can hydrogen bond to two bases, either O(2) of thymine or N(3) of adenine, located on adjacent base pairs on opposite strands of the helix. In most (but not all) cases the tri- and pentapeptide can adopt two orientations at each A+T rich binding site.     

Theoretical study on the binding mechanism between N6-methyladenine and natural DNA bases

N6-methyladenine (m⁶A) is a rare base naturally occurring in DNA. It is different from the base adenine due to its N-CH₃. Therefore, the base not only pairs with thymine, but also with other DNA bases (cytosine, adenine and guanine). In this work, Møller-Plesset second-order (MP2) method has been used to investigate the binding mechanism between m⁶A and natural DNA bases in gas phase and in aqueous solution. The results show that N-CH₃ changed the way of N6-methyladenine binding to natural DNA bases. The binding style significantly influences the stability of base pairs. The trans-m⁶A:G and trans-m⁶A:C conformers are the most stable among all the base pairs. The existence of solvent can remarkably reduce the stability of the base pairs, and the DNA bases prefer pairing with trans-m⁶A to cis-m⁶A. Besides, the properties of these hydrogen bonds have been analyzed by atom in molecules (AIM) theory, natural bond orbital (NBO) analysis and Wiberg bond indexes (WBI). In addition, pairing with m⁶A decreases the binding energies compared to the normal Watson-Crick base pairs, it may explain the instability of the N6 site methylated DNA in theory.

Oligonucleotide probes containing pyrimidine analogs reveal diminished hydrogen bonding capacity of the DNA adduct O6-methyl-G in DNA duplexes

Oligonucleotide hybridization probes containing nucleoside analogs offer a potential strategy for binding specific DNA sequences that bear pro-mutagenic O⁶-G alkylation adducts. To optimize O⁶-Me-G-targeting probes, an understanding of how base pairs with O⁶-Me-G are stabilized is needed. In this study, we compared the ability of O⁶-Me-G and G to hydrogen bond with three pyrimidine-like nucleobases (Z, 4-thio-U, and 3-deaza-C) bearing varied hydrogen bond donor and acceptor groups. We found that duplexes containing the pyrimidine analog nucleoside:G pairs were more thermodynamically stable than those containing pyrimidine analog nucleoside:O⁶-alkyl-G pairs. Thus, hydrogen bonding alone was not sufficient to impart selectivity to probes that target O⁶-G alkylation adducts in DNA.     

Structural insights into the anti-cancer activity of quercetin on G-tetrad,mixed G-tetrad,and G-quadruplex DNA using quantum chemical and molecular dynamics simulations

Abstract

Human telomerase referred as ‘terminal transferase’ is a nucleoprotein enzyme which inhibits the disintegration of telomere length and act as a drug target for the anticancer therapy. The tandem repeating structure of telomere sequence forms the guanine-rich quadruplex structures that stabilize stacked tetrads. In our present work, we have investigated the interaction of quercetin with DNA tetrads using DFT. Geometrical analysis revealed that the influence of quercetin drug induces the structural changes into the DNA tetrads. Among DNA tetrads, the quercetin stacked with GCGC tetrad has the highest interaction energy of ?88.08?kcal/mol. The binding mode and the structural stability are verified by the absorption spectroscopy method. The longer wavelength was found at 380?nm and it exhibits bathochromic shift. The findings help us to understand the binding nature of quercetin drug with DNA tetrads and it also inhibits the telomerase activity. Further, the quercetin drug interacted with G-quadruplex DNA by using molecular dynamics (MD) simulation studies for 100?ns simulation at different temperatures and different pH levels (T?=?298 K, 320?K and pH = 7.4, 5.4). The structural stability of the quercetin with G-quadruplex structure is confirmed by RMSD. For the acidic condition (pH = 5.4), the binding affinity is higher toward G-quadruplex DNA, this result resembles that the quercetin drug is well interacted with G-quadruplex DNA at acidic condition (pH = 7.4) than the neutral condition. The obtained results show that quercetin drug stabilizes the G-quadruplex DNA, which regulates telomerase enzyme and it potentially acts as a novel anti-cancer agent.

Communicated by Ramaswamy H. Sarma     

Systematic Analysis of Possible Hydrogen Bonds between Amino Acid Side Chains and B-form DNA

Abstract

The ways in which amino acid side chains could make a pair of hydrogen bonds within the major groove of B DNA are systematically analyzed. Hydrogen bond donors within the major groove are characterized by determining the idealized position of the hydrogen bond acceptors that they might bond with. It appears that an amino acid side chain could, at most, contact a pair of base pairs. The ten possible pairs of base pairs are analyzed to determine how they could be recognized by the amino acid side chains.     

Effects of deficient of the Hoogsteen base-pairs on the G-quadruplex stabilization and binding mode of a cationic porphyrin

BackgroundIn stabilization of the G-quadruplex, formation of a Hoogsteen base-pair between the guanine (G) bases is essential. However, the contribution of each Hoogsteen base-pair at different positions to whole stability of the G-quadruplex has not been known. In this study, the effect of a deficiency of the Hoogsteen type hydrogen bond in the G-quadruplex stability was investigated. Spectral properties of meso-tetrakis(1-methylpyridinium-4-yl)porphyrin (TMPyP) associated with various G-quadruplexes were also examined.MethodsThe thermal stability of the thrombin-binding DNA aptamer 5′G₁G₂TTG₅G₆TG₈TG₁₀G₁₁TTG₁₄G₁₅ G-quadruplex, in which the guanine (G) base at 1, 2, 5, 6 and 8th positions was replaced with an inosine (I) base, one at a time, was investigated by circular dichroism (CD). The absorption, CD and fluorescence decay curve for the G-quadruplex associated TMPyP were also measured.ResultsThe transition from the G-quadruplex to a single stranded form was endothermic and induced by an increase in entropy. The order in stability was 0>8>6>2>5>1, where the numbers denote the position of the replacement and 0 represents no replacements of the G base, suggesting the significant contribution of the G₁ base in the stability of the G-quadruplex. Alteration in the spectral property of TMPyP briefly followed the order in thermal stability.ConclusionsReplacement of a G base with an I base resulted in destabilization of the G-quadruplex. The missing hydrogen bond at position 1 destabilized the G-quadruplex most efficiently. TMPyP binds near the I base-replaced location namely, the side of the G-quadruplex.General significanceThe Hoogsteen base-pairing is confirmed to be essential in stabilization of G-quadruplex. When G is replaced with I, the latter base is mobile to interact with cationic porphyrin.     

Silver(I)-mediated Hoogsteen-type base pairs

Metal-mediated Hoogsteen-type base pairs are useful for the construction of DNA duplexes containing contiguous stretches of metal ions along the helical axis. To fine-tune the stability of such base pairs and the selectivity toward different metal ions, the availability of a selection of artificial nucleobases is highly desirable. In this study, we follow a theoretical approach utilizing dispersion-corrected density functional methods to evaluate a variety of artificial nucleobases as candidates for metal-mediated Hoogsteen-type base pairs. We focus on silver(I)-mediated Hoogsteen- and reverse Hoogsteen-type base pairs formed between 1-deaza- and 1,3-dideazapurine-derived nucleobases, respectively, and cytosine. Apart from two coordinative bonds, these base pairs are stabilized by a hydrogen bond. We elucidate the impact of different substituents at the C6 position and the presence or absence of an endocyclic N3 nitrogen atom on the overall stability of a base pair and concomitantly on the strength of the hydrogen and coordinative bonds. All artificial base pairs investigated in this study are less stable than the experimentally established benchmark base pair C–Ag⁺–G. The base pair formed from 1,3-dideaza-6-methoxypurine is isoenergetic to the experimentally observed C–Ag⁺–C base pair. This makes 1,3-dideaza-6-methoxypurine a promising candidate for the use as an artificial nucleobase in DNA.     

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.     

Chemical consequences of incorporation of 5-fluorouracil into DNA as studied by NMR

The cytotoxic analogue of thymine, 5-fluorouracil (Uf), is known to be incorporated into DNA in biological systems. This abnormal base has been synthetically incorporated into short DNA oligomers. The ionization of the N-3 proton of this base within DNA oligomers was measured by observation of the 19F chemical shift at varying pH values. The pKa values for the Uf ring of dTpdUfpdT and dApdUfpdA were determined to be 7.84 and 7.9, respectively. The self-complementary 12-mers d(G-C-G-C-A-A-T-Uf-G-C-G-C) and d(C-G-A-T-Uf-A-T-A-A-T-C-G) were synthesized, and 1H NMR was used to compare the helix dynamics and stability of the interstrand imino proton hydrogen bonds with those of the 12-mers d(G-C-G-C-A-A-T-T-G-C-G-C) and d(C-G-A-T-T-A-T-A-A-T-C-G). The N-3 hydrogen bond of the A-Uf base pair was less stable than the corresponding hydrogen bond in A-T base pairs in the same helix, and the A-Uf base pair was less stable than the A-T base pair in the analogous position of the control helix. The observed temperature-dependent dynamics and NMR melting temperatures of the control and dUf-containing oligomers were similar.     

Studies On The Cross Strand Hydrogen Bonds In DNA Double Helices

Abstract

A systematic analysis has been carried out to examine all the stereochemically possible bifurcated hydrogen bonds including those of cross strand type between propeller twisted base pairs in DNA double helices by stereochemical considerations involving base pairs alone and by molecular mechanics studies on dimer and trimer duplexes. The results show that there are limited number of combinations of adjacent base pairs that would facilitate bifurcated cross- strand hydrogen bond (CSH). B-type helices concomitant with negative propeller twist seem to be more favored for the occurrence of CSH than canonical A-type helices because of slide in the latter. The results also demonstrate that helices with appropriate sequences may possess continuous run of these propeller twist driven cross strand hydrogen bonds indicating that they may infact be considered as yet another general structural feature of DNA helices.     

NMR based structural studies decipher stacking of the alkaloid coralyne to terminal guanines at two different sites in parallel G-quadruplex DNA, [d(TTGGGGT)]4 and [d(TTAGGGT)]4

BackgroundTelomere elongation by telomerase gets inhibited by G-quadruplex DNA found in its guanine rich region. Stabilization of G-quadruplex DNA upon ligand binding has evolved as a promising strategy to target cancer cells in which telomerase is over expressed.MethodsInteraction of anti-leukemic alkaloid, coralyne, to tetrameric parallel [d(TTGGGGT)]₄ (Ttel7), [d(TTAGGGT)]₄ (Htel7) and monomeric anti-parallel [dGGGG(TTGGGG)₃] (Ttel22) G-quadruplex DNA has been studied using Circular Dichroism (CD) spectroscopy. Titrations of coralyne with Ttel7 and Htel7 were monitored by ¹H and ³¹P NMR spectroscopy. Solution structure of coralyne-Ttel7 complex was obtained by restrained Molecular Dynamics (rMD) simulations using distance restraints from 2D NOESY spectra. Thermal stabilization of DNA was determined by absorption, CD and ¹H NMR.Results and conclusionsBinding of coralyne to Ttel7/Htel7 induces negative CD band at 315/300 nm. A significant upfield shift in all GNH, downfield shift in T2/T7 base protons and upfield shift (1.8 ppm) in coralyne protons indicates stacking interactions. ³¹P chemical shifts and NOE contacts of G3, G6, T2, T7 protons with methoxy protons reveal proximity of coralyne to T2pG3 and G6pT7 sites. Solution structure reveals stacking of coralyne at G6pT7 and T2pG3 steps with two methoxy groups of coralyne located in the grooves along with formation of a hydrogen bond. Binding stabilizes Ttel7/Htel7 by ~ 25–35 °C in 2:1 coralyne-Ttel7/Htel7 complex.General significanceThe present study is the first report on solution structure of coralyne-Ttel7 complex showing stacking of coralyne with terminal guanine tetrads leading to significant thermal stabilization, which may be responsible for telomerase inhibition.     

Density functional theory study on the interaction between keto-9H guanine and aspartic acid

A theoretical study was performed using density functional theory (DFT) to investigate hydrogen bonding interactions in signature complexes formed between keto-9H guanine (Gua) and aspartic acid (Asp) at neutral pH. Optimized geometries, binding energies and the theoretical IR spectra of guanine, aspartic acid and their corresponding complexes (Gua-Asp) were calculated using the B3LYP method and the 6-31+G(d) basis set. Stationary points found to be at local minima on the potential energy surface were verified by second derivative harmonic vibrational frequency calculations at the same level of theory. AIM theory was used to analyze the hydrogen bonding characteristics of these DNA base complex systems. Our results show that the binding motif for the most stable complex is strikingly similar to a Watson-Crick motif observed in the guanine-cytosine base pair. We have found a range of hydrogen bonding interactions between guanine and aspartic acid in the six complexes. This was further verified by theoretical IR spectra of ω(C-H---O-H) cm⁻¹ stretches for the Gua-Asp complexes. The electron density plot indicates strong hydrogen bonding as shown by the 2p _z dominant HOMO orbital character.     

Quantum Chemical Studies of Structures and Binding in Noncanonical RNA Base pairs: The Trans Watson-Crick:Watson-Crick Family

Abstract

The trans Watson-Crick/Watson-Crick family of base pairs represent a geometric class that play important structural and possible functional roles in the ribosome, tRNA, and other functional RNA molecules. They nucleate base triplets and quartets, participate as loop closing terminal base pairs in hair pin motifs and are also responsible for several tertiary interactions that enable sequentially distant regions to interact with each other in RNA molecules. Eleven representative examples spanning nine systems belonging to this geometric family of RNA base pairs, having widely different occurrence statistics in the PDB database, were studied at the HF/6–31G (d, p) level using Morokuma decomposition, Atoms in Molecules as well as Natural Bond Orbital methods in the optimized gas phase geometries and in their crystal structure geometries, respectively. The BSSE and deformation energy corrected interaction energy values for the optimized geometries are compared with the corresponding values in the crystal geometries of the base pairs. For non protonated base pairs in their optimized geometry, these values ranged from ?8.19 kcal/mol to ?21.84 kcal/mol and compared favorably with those of canonical base pairs. The interaction energies of these base pairs, in their respective crystal geometries, were, however, lesser to varying extents and in one case, that of A:A W:W trans, it was actually found to be positive. The variation in RMSD between the two geometries was also large and ranged from 0.32–2.19 Å. Our analysis shows that the hydrogen bonding characteristics and interaction energies obtained, correlated with the nature and type of hydrogen bonds between base pairs; but the occurrence frequencies, interaction energies, and geometric variabilities were conspicuous by the absence of any apparent correlation. Instead, the nature of local interaction energy hyperspace of different base pairs as inferred from the degree of their respective geometric variability could be correlated with the identities of free and bound hydrogen bond donor/acceptor groups present in interacting bases in conjunction with their tertiary and neighboring group interaction potentials in the global context. It also suggests that the concept of isostericity alone may not always determine covariation potentials for base pairs, particularly for those which may be important for RNA dynamics. These considerations are more important than the absolute values of the interaction energies in their respective optimized geometries in rationalizing their occurrences in functional RNAs. They highlight the importance of revising some of the existing DNA based structure analysis approaches and may have significant implications for RNA structure and dynamics, especially in the context of structure prediction algorithms.     

Investigation on the interaction of letrozole with herring sperm DNA through spectroscopic and modeling methods

The interaction of letrozole, an efficient and safe aromatase inhibitor, with herring sperm DNA (hsDNA) was investigated in vitro through spectroscopy analysis and molecular modeling to elucidate the binding mechanism of anticancer drugs and DNA. The binding constant and the number of binding sites were 2.13 × 10⁴ M^?1 and 1.09, respectively, at 298 K. Thermodynamic parameters (ΔG, ΔH and ΔS) exhibited negative values, which indicated that binding was spontaneous and Van der Waals forces and hydrogen bond were the main interaction forces. Fourier transform infrared spectroscopy and other spectroscopy analysis methods illustrated that letrozole could intercalate into the phosphate backbone of hsDNA and interact with the nitrogenous bases. Consistent with the experimental findings, molecular modeling results demonstrated that the interaction was dominated by intercalation and hydrogen bonding. Copyright © 2015 John Wiley & Sons, Ltd.     

Intercalation of a Zn(II) complex containing ciprofloxacin drug between DNA base pairs

In this study, an attempt has been made to study the interaction of a Zn(II) complex containing an antibiotic drug, ciprofloxacin, with calf thymus DNA using spectroscopic methods. It was found that Zn(II) complex could bind with DNA via intercalation mode as evidenced by: hyperchromism in UV–Vis spectrum; these spectral characteristics suggest that the Zn(II) complex interacts with DNA most likely through a mode that involves a stacking interaction between the aromatic chromophore and the base pairs of DNA. DNA binding constant (K_b = 1.4 × 10⁴ M^?1) from spectrophotometric studies of the interaction of Zn(II) complex with DNA is comparable to those of some DNA intercalative polypyridyl Ru(II) complexes 1.0 ?4.8 × 10⁴ M^?1. CD study showed stabilization of the right-handed B form of DNA in the presence of Zn(II) complex as observed for the classical intercalator methylene blue. Thermodynamic parameters (ΔH < 0 and ΔS < 0) indicated that hydrogen bond and Van der Waals play main roles in this binding prose. Competitive fluorimetric studies with methylene blue (MB) dye have shown that Zn(II) complex exhibits the ability of this complex to displace with DNA-MB, indicating that it binds to DNA in strong competition with MB for the intercalation.     

Rational Design of Pyrrolo[1,2-a]benzimidazole-Based Antitumor Agents Targeting the DNA Major Groove

Models have been developed for the interaction of the pyrrolo[1,2-a]benzimidazole (PBI) antitumor agents with the two-electron activating enzyme DT-diaphorase and the DNA major groove. The DT-diaphorase model and experimental results indicate that the S-enantiomer of 3-carbamido PBI can enantioselect ovarian cancers. The reduced PBI interacts with the DNA major groove at AT base pairs by forming Hoogsteen-like hydrogen bonds. The reduced 3-amino PBI forms three hydrogen bonds in the major groove with the amino group acting as an H-bond donor to the thymine carbonyl. The DNA-binding model will permit the design of major groove recognition agents.