Metal complex–DNA binding: Insights from molecular dynamics and DFT/MM calculations

In-vitro and in-silico analysis and antitumor studies of novel Cu(II) and V(V) complexes of N-p-Tolylbenzohydroxamic acid

Copper(L2Cu) and vanadium(L2VOCl) complexes of N-p-tolylbenzohydroxamic acid (LH) ligand have been investigated for DNA binding efficacy by multiple analytical, spectral, and computational techniques. The results revealed that complexes as groove binders as evidenced by UV absorption. Fluorescence studies including displacement assay using classical intercalator ethidium bromide as fluorescent probe also confirmed as groove binders. The viscometric analysis too supports the inferences as strong groove binders for both the complexes. Molecular docking too exposed DNA as a target to the complexes which precisely binds L2Cu, in the minor groove region while L2VOCl in major groove region. Molecular dynamic simulation performed on L2Cu complex revealing the interaction of complex with DNA within 20 ns time. The complex stacked into the nitrogen bases of oligonucleotides and the bonding features were intrinsically preserved for longer simulation times. In-vitro cytotoxicity study was undertaken employing MTT assay against the breast cancer cell line (MCF-7). Potential cytotoxic activities were observed for L2Cu and L2VOCl complexes with IC50 values of showing 71 % and 74 % of inhibition respectively.

An insight into the mixed quantum mechanical-molecular dynamic simulation of a ZnII-Curcumin complex with a chosen DNA sequence that supports experimental DNA binding investigations

An important aspect of research pertaining to Curcumin (HCur) is the need to arrest its degradation in aqueous solution and in biological milieu. This may be achieved through complex formation with metal ions. For this reason, a complex of HCur was prepared with Zn^II, that is not likely to be active in redox pathways, minimizing further complications. The complex is monomeric, tetrahedral, with one HCur, an acetate and a molecule of water bound to Zn^II. It arrests degradation of HCur to a considerable extent that was realized by taking it in phosphate buffer and in biological milieu. The structure was obtained by DFT calculations. Stable adduct formation was identified between optimized structures of HCur and [Zn(Cur)] with DNA (PDB ID: 1BNA) through experiments validated with multiscale modeling approach. Molecular docking studies provide 2D and 3D representations of binding of HCur and [Zn(Cur)] through different non-covalent interactions with the nucleotides of the chosen DNA. Through molecular dynamics simulation, a detailed understanding of binding pattern and key structural characteristics of the generated DNA-complex was obtained following analysis by RMSD, RMSF, radius of gyration, SASA and aspects like formation of hydrogen bonds. Experimental studies provide binding constants for [Zn(Cur)] with calf thymus DNA at 25 °C that effectively helps one to realize its high affinity towards DNA. In the absence of an experimental binding study of HCur with DNA, owing to its tendency to degrade in solution, a theoretical analysis of the binding of HCur to DNA is extremely helpful. Besides, both experimental and simulated binding of [Zn(Cur)] to DNA may be considered as a case of pseudo-binding of HCur to DNA. In a way, such studies on interaction with DNA helps one to identify HCur's affinity for cellular target DNA, not realized through experiments. The entire investigation is an understanding of experimental and theoretical approaches that has been compared continuously, being particularly useful when a molecule's interaction with a biological target cannot realized experimentally.

Intercalation Ability of Novel Monofunctional Platinum Anticancer Drugs: A Key Step in Their Biological Action

Phenanthriplatin (PtPPH) is a monovalent platinum(II)-based complex with a large cytotoxicity against cancer cells. Although the aqua-activated drug has been assumed to be the precursor for DNA damage, it is still under debate whether the way in which that metallodrug attacks to DNA is dominated by a direct binding to a guanine base or rather by an intercalated intermediate product. Aiming to capture the mechanism of action of PtPPH, the present contribution used theoretical tools to systematically assess the sequence of all possible mechanisms on drug activation and reactivity, for example, hydrolysis, intercalation, and covalent damage to DNA. Ab initio quantum mechanical (QM) methods, hybrid QM/QM′ schemes, and independent gradient model approaches are implemented in an unbiased protocol. The performed simulations show that the cascade of reactions is articulated in three well-defined stages: (i) an early and fast intercalation of the complex between the DNA bases, (ii) a subsequent hydrolysis reaction that leads to the aqua-activated form, and (iii) a final formation of the covalent bond between PtPPH and DNA at a guanine site. The permanent damage to DNA is consequently driven by that latter bond to DNA but with a simultaneous π–π intercalation of the phenanthridine into nucleobases. The impact of the DNA sequence and the lateral backbone was also discussed to provide a more complete picture of the forces that anchor the drug into the double helix.

Organoplatinum(II) Complexes Self-Assemble and Recognize AT-Rich Duplex DNA Sequences

The specific recognition of AT-rich DNA sequences opens up the door to promising diagnostic and/or therapeutic strategies against gene-related diseases. Here, we demonstrate that amphiphilic Pt^II complexes of the type [Pt(dmba)(N∧N)]NO₃ (dmba = N,N-dimethylbenzylamine-κN, κC; N∧N = dpq (3), dppz (4), and dppn (5)) recognize AT-rich oligonucleotides over other types of DNA, RNA, and model proteins. The crystal structure of 4 shows the presence of significant π-stacking interactions and a distorted coordination sphere of the d⁸ Pt^II atom. Complex 5, containing the largest π-conjugated ligand, forms supramolecular assemblies at high concentrations under aqueous environment. However, its aggregation can be promoted in the presence of DNA at concentrations as low as 10 μM in a process that “turns on” its excimer emission around 600 nm. Viscometry, gel electrophoresis, and theoretical calculations demonstrate that 5 binds to minor groove when self-assembled, while the monomers of 3 and 4 intercalate into the DNA. The complexes also inhibit cancer cell growth with low-micromolar IC₅₀ values in 2D tissue culture and suppress tumor growth in 3D tumor spheroids with a multicellular resistance (MCR) index comparable to that of cisplatin.

Cyclometalated Pt^II complexes containing π-conjugated ligands form supramolecular assemblies under aqueous environment, and DNA-induced aggregation occurs for the one containing the highest conjugated N,N-diimine ligand. The complexes recognize AT-rich DNA sequences over others in DNA, RNA, and proteins. Their DNA binding mode switches from intercalation to minor groove binding when self-assembled. The complexes suppress tumor growth in 3D tumor spheroids.

Human DNA Telomeres in Presence of Oxidative Lesions: The Crucial Role of Electrostatic Interactions on the Stability of Guanine Quadruplexes

By using all atom molecular dynamics simulations, we studied the behavior of human DNA telomere sequences in guanine quadruplex (G4) conformation and in the presence of oxidative lesions, namely abasic sites. In particular, we evidenced that while removing one guanine base induces a significant alteration and destabilization of the involved leaflet, human telomere oligomers tend, in most cases, to maintain at least a partial quadruplex structure, eventually by replacing the empty site with undamaged guanines of different leaflets. This study shows that (i) the disruption of the quadruplex leaflets induces the release of at least one of the potassium cations embedded in the quadruplex channel and that (ii) the electrostatic interactions of the DNA sequence with the aforementioned cations are fundamental to the maintenance of the global quadruplex structure.

Interaction of osmium(ii) redox probes with DNA: insights from theory

In the course of developing ultrasensitive and quantitative electrochemical point-of-care analytical tools for genetic detection of infectious diseases, osmium(ii) metallointercalators were revealed to be suitable and efficient redox probes to monitor the in vitro DNA amplification [Defever etal, Anal. Chem., 2011, 83, 1815-1821]. In this work, we thus propose a complete computational protocol in order to evaluate the affinity between Os(ii) complexes with double-stranded DNA. This protocol is based on molecular dynamics, with the parametrization of the GAFF force field for the Os(ii) complexes presenting an octahedral environment with polypyridine ligands, and QM/QM' calculations to evaluate the binding energy. For three Os(ii) probes and different binding sites, molecular dynamics simulations and interaction energies calculated at the QM/QM' level are successively discussed and compared to experimental data in order to identify the most stable binding sites. The computational protocol we propose should then be used to design more efficient Os(ii) metallointercalators.

Interaction of Cd(ii) and Ni(ii) terpyridine complexes with model polynucleotides: a multidisciplinary approach

The study of the intercalation of both complexes, evidenced by CD and fluorescence spectroscopy and supported by QM/MM calculations, broadens the experimental and theoretical background on drugs/DNA interactions.

Theoretical studies on binding modes of copper-based nucleases with DNA

In the present work, molecular simulations were performed for the purpose of predicting the binding modes of four types of copper nucleases (a total 33 compounds) with DNA. Our docking results accurately predicted the groove binding and electrostatic interaction for some copper nucleases with B-DNA. The intercalation modes were also reproduced by "gap DNA". The obtained results demonstrated that the ligand size, length, functional groups and chelate ring size bound to the copper center could influence the binding affinities of copper nucleases. The binding affinities obtained from the docking calculations herein also replicated results found using MM-PBSA approach. The predicted DNA binding modes of copper nucleases with DNA will ultimately help us to better understand the interaction of copper compounds with DNA.

Intercalation processes of copper complexes in DNA

The family of anticancer complexes that include the transition metal copper known as Casiopeínas® shows promising results. Two of these complexes are currently in clinical trials. The interaction of these compounds with DNA has been observed experimentally and several hypotheses regarding the mechanism of action have been developed, and these include the generation of reactive oxygen species, phosphate hydrolysis and/or base-pair intercalation. To advance in the understanding on how these ligands interact with DNA, we present a molecular dynamics study of 21 Casiopeínas with a DNA dodecamer using 10 μs of simulation time for each compound. All the complexes were manually inserted into the minor groove as the starting point of the simulations. The binding energy of each complex and the observed representative type of interaction between the ligand and the DNA is reported. With this extended sampling time, we found that four of the compounds spontaneously flipped open a base pair and moved inside the resulting cavity and four compounds formed stacking interactions with the terminal base pairs. The complexes that formed the intercalation pocket led to more stable interactions.

Photophysics of acetophenone interacting with DNA: why the road to photosensitization is open

Deoxyribonucleic acid photosensitization, i.e. the photoinduced electron- or energy-transfer of chromophores interacting with DNA, is a crucial phenomenon that triggers important DNA lesions such as pyrimidine dimerization, even upon absorption of relatively low-energy radiation. Oxidative lesions may also be produced via the photoinduced production of reactive oxygen species. Aromatic ketones, and acetophenone in particular, are well known for their sensitization effects. In this contribution we model the structural and dynamical properties of the acetophenone/DNA aggregates as well as their spectroscopic and photophysical properties using high-level hybrid quantum mechanics/molecular mechanics methods. We show that the key steps of the photochemistry of acetophenone in gas phase are conserved in the macromolecular environment and thus an ultrafast singlet-triplet conversion of acetophenone is expected prior to the transfer to DNA.

Interaction of doxorubicin with polynucleotides. A spectroscopic study

The interaction of doxorubicin (DX) with model polynucleotides poly(dG-dC)·poly(dG-dC) (polyGC), poly(dA-dT)·poly(dA-dT) (polyAT), and calf thymus DNA has been studied by several spectroscopic techniques in phosphate buffer aqueous solutions. UV-vis, circular dichroism, and fluorescence spectroscopic data confirm that intercalation is the prevailing mode of interaction, and also reveal that the interaction with AT-rich regions leads to the transfer of excitation energy to DX not previously documented in the literature. Moreover, the DX affinity for AT sites has been found to be on the same order of magnitude as that reported for GC sites.

Toward a rationale for the PTC124 (Ataluren) promoted readthrough of premature stop codons: a computational approach and GFP-reporter cell-based assay

The presence in the mRNA of premature stop codons (PTCs) results in protein truncation responsible for several inherited (genetic) diseases. A well-known example of these diseases is cystic fibrosis (CF), where approximately 10% (worldwide) of patients have nonsense mutations in the CF transmembrane regulator (CFTR) gene. PTC124 (3-(5-(2-fluorophenyl)-1,2,4-oxadiazol-3-yl)-benzoic acid), also known as Ataluren, is a small molecule that has been suggested to allow PTC readthrough even though its target has yet to be identified. In the lack of a general consensus about its mechanism of action, we experimentally tested the ability of PTC124 to promote the readthrough of premature termination codons by using a new reporter. The reporter vector was based on a plasmid harboring the H2B histone coding sequence fused in frame with the green fluorescent protein (GFP) cDNA, and a TGA stop codon was introduced in the H2B-GFP gene by site-directed mutagenesis. Additionally, an unprecedented computational study on the putative supramolecular interaction between PTC124 and an 11-codon (33-nucleotides) sequence corresponding to a CFTR mRNA fragment containing a central UGA nonsense mutation showed a specific interaction between PTC124 and the UGA codon. Altogether, the H2B-GFP-opal based assay and the molecular dynamics (MD) simulation support the hypothesis that PTC124 is able to promote the specific readthrough of internal TGA premature stop codons.

Selective G-quadruplex stabilizers: Schiff-base metal complexes with anticancer activity

Molecular dynamics simulations and quantum mechanics/molecular mechanics calculations provided a mechanism for G-quadruplex binding of three transition metal complexes.

B-DNA Structure and Stability as Function of Nucleic Acid Composition: Dispersion-Corrected DFT Study of Dinucleoside Monophosphate Single and Double Strands

We have computationally investigated the structure and stability of all 16 combinations of two out of the four natural DNA bases A, T, G and C in a di-2′-deoxyribonucleoside-monophosphate model DNA strand as well as in 10 double-strand model complexes thereof, using dispersion-corrected density functional theory (DFT-D). Optimized geometries with B-DNA conformation were obtained through the inclusion of implicit water solvent and, in the DNA models, of sodium counterions, to neutralize the negative charge of the phosphate groups. The results obtained allowed us to compare the relative stability of isomeric single and double strands. Moreover, the energy of the Watson–Crick pairing of complementary single strands to form double-helical structures was calculated. The latter furnished the following increasing stability trend of the double-helix formation energy: d(TpA)₂ <d(CpA)₂ <d(ApT)₂ <d(ApA)₂ <d(GpT)₂ <d(GpA)₂ <d(ApG)₂ <d(CpG)₂ <d(GpG)₂ <d(GpC)₂, where the energy differences between the last four dimers, d(ApG)₂, d(CpG)₂, d(GpG)₂ and d(GpC)₂, is within 4.0 kcal mol⁻¹, and the energy between the most and the least stable isomers is 13.4 kcal mol⁻¹. This trend shows that the formation energy essentially increases with the number of hydrogen bonds per base pair, that is two between A and T and three between G and C. Superimposed on this main trend are more subtle effects that depend on the order in which bases occur within a strand from the 5’- to the 3’-end.