Docking studies on DNA-ligand interactions: building and application of a protocol to identify the binding mode

In-silico modeling studies of G-quadruplex with soy isoflavones having anticancerous activity

Telomere forms t-loop and G-quadruplex as the protective structure and the formation of these structures hinder the telomerase enzyme action. The binding affinities of ligand which stabilize the G-quadruplex represent good correlation with telomerase inhibition depicted in the anti-cancerous action. Most of the potent G-quadruplex stabilizing compounds suffer from the poor drug like properties. Herein, natural dietary compounds isoflavones were taken for the theoretical study to examine their stabilizing effect on G-quadruplex structure. The experimental G-quadruplex complexes were reproduced to obtain and validate the theoretical parameters. The obtained theoretical binding energies are in significant correlation with the experimental data. Analysis of binding shows isoflavones to be groove binders, and differential nature of quadruplex grooves might be beneficial in the selectivity aspects. Among all, derrubone was found to have better selectivity as well as affinity for the G-quadruplex comparable to well known ligand TMPyP4. The GBSA rescoring result enlightens the various interaction terms involved in the binding process. Cumulative stabilizing effects coming from VDW, ES, and GB energy terms attest to optimal binding of derrubone molecule which can be considered as a lead for the higher phases of drug designing. These findings are of great value in terms of unexplored groove binding modes and the studied natural compounds might be helpful to direct the focus of synthetic chemists in designing of new generation of antitumor agents.

Novel insights into the apoptosis mechanism of DNA topoisomerase I inhibitor isoliquiritigenin on HCC tumor cell

The inhibitory effect of DNA topoisomerase (Top I) by isoliquiritigenin(ISO) were investigated and their interaction mechanism was evaluated using methods including UV-vis absorption, fluorescence, coupled with molecular simulation, and using the MTT method of inhibition rate of HCC tumor cell SNU475 proliferation assay, finally, the interaction of ISO with calf thymus DNA was investigated by melting measurements and molecular docking studies. It was found that isoliquiritigenin reversibly inhibited DNA Top I in a competitive manner with the concentrations of ISO resulting in 50% activity lost (IC50) were estimated to be 0.178 ± 0.12 mM. Isoliquiritigenin exhibited a strong ability to quench the intrinsic fluorescence of Top I through a static quenching procedure. The positive values of enthalpy change and entropy change suggested that the binding of isoliquiritigenin to Top I was driven mainly by hydrophobic interactions. The molecular docking results revealed isoliquiritigenin actually interacted with the primary amino acid residues on the active site of Top I, and the detection results of fluorescence staining and the inhibitory effect on the growth of HCC SUN475 showed that isoliquiritigenin induced the apoptosis cells increased gradually. The interaction of ISO with DNA can cause the denaturation temperature to be increased, which indicated that the stabilization of the DNA helix was increased in the presence of ISO, which indicated that the results provide strong evidence for intercalative binding of ISO with DNA.

Water-soluble DNA minor groove binders as potential chemotherapeutic agents: synthesis, characterization, DNA binding and cleavage, antioxidation, cytotoxicity and HSA interactions

Two new water-soluble copper(ii)-dipeptide complexes: [Cu(glygly)(PyTA)]ClO4·1.5H2O (1) and [Cu(glygly)(PzTA)]ClO4·1.5H2O (2) (glygly = glycylglycine anion, PyTA = 2,4-diamino-6-(2'-pyridyl)-1,3,5-triazine and PzTA = 2,4-diamino-6-(2'-pyrazino)-1,3,5-triazine), utilizing two interrelated DNA base-like ligands (PyTA and PzTA), have been synthesized and characterized. The structure elucidation for 1 performed by single crystal X-ray diffraction showed a one dimensional chain conformation in which the central copper ions arrange in a five-coordinate distorted square-pyramidal geometry. Spectroscopic titration, viscosity and electrophoresis measurements revealed that the complexes bound to DNA via an outside groove binding mode, and cleaved pBR322 DNA efficiently in the presence of ascorbate, probably via an oxidative mechanism with the involvement of ˙OH and ˙O2(-). Notably, the complexes exhibited considerable in vitro cytotoxicity against four human carcinoma cell lines (HepG2, HeLa, A549 and U87) with IC50 values ranging from 41.68 to 159.17 μM, in addition to their excellent SOD mimics (IC50 ~ 0.091 and 0.114 μM). Besides, multispectroscopic evidence suggested their HSA-binding at the cavity containing Trp-214 in subdomain IIA with moderate affinity, mainly via hydrophobic interaction. Further, the molecular docking technique utilized for ascertaining the mechanism and mode of action towards DNA and HSA theoretically verified the experimental results.

Computational docking simulations of a DNA-aptamer for argininamide and related ligands

The binding properties of sequence-specific nucleic acids (aptamers) to low-molecular-weight ligands, macromolecules and even cells attract substantial scientific interest. These ligand-DNA complexes found different applications for sensing, nanomedicine, and DNA nanotechnology. Structural information on the aptamer-ligand complexes is, however, scarce, even though it would open-up the possibilities to design novel features in the complexes. In the present study we apply molecular docking simulations to probe the features of an experimentally documented L-argininamide aptamer complex. The docking simulations were performed using AutoDock 4.0 and YASARA Structure software, a well-suited program for following intermolecular interactions and structures of biomolecules, including DNA. We explored the binding features of a DNA aptamer to L-argininamide and to a series of arginine derivatives or arginine-like ligands. We find that the best docking results are obtained after an energy-minimization of the parent ligand-aptamer complexes. The calculated binding energies of all mono-substituted guanidine-containing ligands show a good correlation with the experimentally determined binding constants. The results provide valuable guidelines for the application of docking simulations for the prediction of aptamer-ligand structures, and for the design of novel features of ligand-aptamer complexes.

Expanding the family of heteroleptic oxidovanadium(IV) compounds with salicylaldehyde semicarbazones and polypyridyl ligands showing anti-Trypanosoma cruzi activity

Searching for prospective vanadium-based drugs for the treatment of Chagas disease, a new series of heteroleptic [V(IV)O(L-2H)(NN)] compounds was developed by including the lipophilic 3,4,7,8-tetramethyl-1,10-phenanthroline (tmp) NN ligand and seven tridentate salicylaldehyde semicarbazone derivatives (L1-L7). The compounds were characterized in the solid state and in solution. EPR spectroscopy suggests that the NN ligand is bidentate bound through both nitrogen donor atoms in an axial-equatorial mode. The EPR and (51)V-NMR spectra of aerated solutions at room temperature indicate that the compounds are stable to hydrolysis and that no significant oxidation of V(IV) to V(V) takes place at least in 24h. The complexes are more active in vitro against Trypanosoma cruzi, the parasite responsible for Chagas disease, than the reference drug Nifurtimox and most of them are more active than previously reported [V(IV)O(L-2H)(NN)] complexes of other NN co-ligands. Selectivity towards the parasite was analyzed using J-774 murine macrophages as mammalian cell model. Due to both, high activity and high selectivity, L2, L4, L5 and L7 complexes could be considered new hits for further drug development. Lipophilicity probably plays a relevant role in the bioactivity of the new compounds. The [V(IV)O(L-2H)(NN)] compounds were designed aiming DNA as potential molecular target. Therefore, the novel L1-L7 tmp complexes were screened by computational modeling, comparing their DNA-binding features with those of previously reported [V(IV)O(L-2H)(NN)] compounds with different NN co-ligands. Whereas all the complexes interact well with DNA, with binding modes and strength tuned in different extents by the NN and semicarbazone co-ligands, molecular docking suggests that the observed anti-T. cruzi activity cannot be explained upon DNA intercalation as the sole mechanism of action.

Validation of a computational docking methodology to identify the non-covalent binding site of ligands to DNA

Despite the biomedical consequences of carcinogen-DNA interactions and the potential of DNA as a drug target in medicinal chemistry, only a small number of studies have validated or used docking methods for the prediction of the physical binding of small molecules to DNA. Knowledge of the DNA-physically-bound ligand geometry can lead to the elucidation of the molecular-level mechanism of drugs as well as predicting the subsequent chemical interactions that lead to DNA damage from carcinogens. We sought to validate AutoDock 4.2, a docking method that includes a physics-based free energy function and a Lamarckian Genetic Algorithm, for the prediction of ligand geometries upon physical binding to DNA. We performed simulations by systematically changing the length of the search process for a comprehensive set of 32 ligand-DNA molecular systems with different physico-chemical properties, and we used a free-energy-based convergence criterion to terminate our simulations. For 11 out of 28 molecular systems for which convergence was achieved, the lowest binding free energy geometries were within 2 Å of the experimentally determined geometry. Considering all predicted sites with free energy changes within 20% of the lowest binding free energy site, we found a site within 2 Å of the experimentally determined geometry for 24 out of the 28 systems. However, the predicted hydrogen bonding interactions were different for most molecular systems compared to the same interactions in the experimentally determined geometry. We discuss reasons for the successes and failures, implications, and the importance of ensuring an adequate search in docking calculations. Overall, we concluded that AutoDock 4.2 can be used to predict the non-covalent binding geometry of a small molecule to DNA with some limitations.

In search of a docking protocol to distinguish between DNA intercalators and groove binders: genetic algorithm vs. shape-complementarity based docking methods

A comparison of docking of DNA groove binders and intercalators using genetic algorithm and shape-complementarity based methods is presented.

Metal complexes of tosyl sulfonamides: design, X-ray structure, biological activities and molecular docking studies

CompoundNA3bound inside the active site of the enzyme.

New palladium(ii) and platinum(ii) 5,5-diethylbarbiturate complexes with 2-phenylpyridine, 2,2′-bipyridine and 2,2′-dipyridylamine: synthesis, structures, DNA binding, molecular docking, cellular uptake, antioxidant activity and cytotoxicity

DNA interaction, antioxidant and cytotoxic activities of Pd(ii)/Pt(ii) 5,5-diethylbarbiturate complexes were evaluated.

Exploration of the binding interaction of a potential nervous system stimulant with calf-thymus DNA and dissociation of the drug–DNA complex by detergent sequestration

The binding interaction of a potential nervous system stimulant with calf-thymus DNA has been divulged and dissociation of the drug–DNA complex has been achieved by the detergent sequestration method.

Evidence for the contribution of non-covalent steroid interactions between DNA and topoisomerase in the genotoxicity of steroids

Fifty two steroids and 9 Vitamin D analogs were docked into ten crystallographically-defined DNA dinucleotide sites and two human topoisomerase II ATP binding sites using two computational programs, Autodock and Surflex. It is shown that both steroids and Vitamin D analogs exhibit a propensity for non-covalent intercalative binding to DNA. A higher predicted binding affinity was found, however, for steroids and the ATP binding site of topoisomerase; in fact these drugs exhibited among the highest topo II binding observed in over 1370 docked drugs. These findings along with genotoxicity data from 26 additional steroids not subjected to docking analysis, support a mechanism wherein the long known, but poorly understood, clastogenicity of steroids may be attributable to inhibition of topoisomerase. A "proof of principle" experiment with dexamethasone demonstrated this to be the likely mechanism of clastogenicity of, at least, this steroid. The generality of this proposed mechanism of genotoxicity across the steroids and vitamin-D analogs is discussed.

Recent advances in small organic molecules as DNA intercalating agents: synthesis, activity, and modeling

The interaction of small molecules with DNA plays an essential role in many biological processes. As DNA is often the target for majority of anticancer and antibiotic drugs, study about the interaction of drug and DNA has a key role in pharmacology. Moreover, understanding the interactions of small molecules with DNA is of prime significance in the rational design of more powerful and selective anticancer agents. Two of the most important and promising targets in cancer chemotherapy include DNA alkylating agents and DNA intercalators. For these last the DNA recognition is a critical step in their anti-tumor action and the intercalation is not only one kind of the interactions in DNA recognition but also a pivotal step of several clinically used anti-tumor drugs such as anthracyclines, acridines and anthraquinones. To push clinical cancer therapy, the discovery of new DNA intercalators has been considered a practical approach and a number of intercalators have been recently reported. The intercalative binding properties of such molecules can also be harnessed as diagnostic probes for DNA structure in addition to DNA-directed therapeutics. Moreover, the problem of intercalation site formation in the undistorted B-DNA of different length and sequence is matter of tremendous importance in molecular modeling studies and, nowadays, three models of DNA intercalation targets have been proposed that account for the binding features of intercalators. Finally, despite DNA being an important target for several drugs, most of the docking programs are validated only for proteins and their ligands. Therefore, a default protocol to identify DNA binding modes which uses a modified canonical DNA as receptor is needed.

Binding of an anionic fluorescent probe with calf thymus DNA and effect of salt on the probe–DNA binding: a spectroscopic and molecular docking investigation

Binding mode of biologically relevant anionic probe, ANS, with ctDNA is divulged from spectroscopic and molecular docking studies.

Probing the binding of insecticide permethrin to calf thymus DNA by spectroscopic techniques merging with chemometrics method

The binding of permethrin (PE) with calf thymus DNA (ctDNA) in physiological buffer (pH 7.4) was investigated by ultraviolet-visible (UV-vis) absorption, fluorescence, circular dichroism (CD), and Fourier transform infrared (FT-IR) spectroscopy merging with multivariate curve resolution-alternating least-squares (MCR-ALS) chemometrics approach. The MCR-ALS was applied to resolve the combined spectroscopic data matrix, which was obtained by UV-vis and fluorescence methods. The concentration profiles of PE, ctDNA, and PE-ctDNA complex and their pure spectra were then successfully obtained. The PE molecular was found to be able to intercalate into the base pairs of ctDNA as evidenced by decreases in resonance light-scattering signal and iodide-quenching effect and increase in ctDNA viscosity. The results of FT-IR spectra indicated that PE was prone to bind to G-C base pairs of ctDNA, and the molecular docking studies were used to validate and clarify the specific binding. The observed changes in CD signals revealed that the DNA turned into a more highly wound form of B-conformation. The calculated thermodynamic parameters, enthalpy change (ΔH°) and entropy change (ΔS°), suggested that hydrogen bonds and van der Waals forces played a predominant role in the binding of PE to ctDNA.

Antimicrobial activity of PVP from an Antarctic bacterium, Janthinobacterium sp. Ant5-2, on multi-drug and methicillin resistant Staphylococcus aureus

Multiple drug resistant (MDR) and methicillin-resistant Staphylococcus aureus (MRSA) have become increasingly prevalent as a community acquired infection. As a result limited treatment options are available with conventional synthetic antibiotics. Bioprospecting natural products with potent antimicrobial activity show promise for developing new drugs against this pathogen. In this study, we have investigated the antimicrobial activity of a purple violet pigment (PVP) from an Antarctic bacterium, Janthinobacterium sp. Ant5-2 on 15 clinical MDR and MRSA strains. The colorimetric resazurin assay was employed to determine the minimum inhibitory concentration (MIC₉₀) of PVP against MDR and MRSA. The MIC₉₀ ranged between 1.57 µg/mL and 3.13 µg/mL, which are significantly lower than many antimicrobials tested from natural sources against this pathogen. The spectrophotometrically determined growth analysis and total microscopic counts using Live/dead® BacLight™ fluorescent stain exhibited a steady decrease in viability of both MDR and MRSA cultures following treatment with PVP at the MIC levels. In silico predictive molecular docking study revealed that PVP could be a DNA-targeting minor groove binding antimicrobial compound. The continued development of novel antimicrobials derived from natural sources with the combination of a suite of conventional antibiotics could stem the rising pandemic of MDR and MRSA along with other deadly microbial pathogens.

Existence of a new emitting singlet state of proflavine: femtosecond dynamics of the excited state processes and quantum chemical studies in different solvents

Proflavine (3,6-diaminoacridine) shows fluorescence emission with lifetime, 4.6 ± 0.2 ns, in all the solvents irrespective of the solvent polarity. To understand this unusual photophysical property, investigations were carried out using steady state and time-resolved fluorescence spectroscopy in the pico- and femtosecond time domain. Molecular geometries in the ground and low-lying excited states of proflavine were examined by complete structural optimization using ab initio quantum chemical computations at HF/6-311++G** and CIS/6-311++G** levels. Time dependent density functional theory (TDDFT) calculations were performed to study the excitation energies in the low-lying excited states. The steady state absorption and emission spectral details of proflavine are found to be influenced by solvents. The femtosecond fluorescence decay of the proflavine in all the solvents follows triexponential function with two ultrafast decay components (τ(1) and τ(2)) in addition to the nanosecond component. The ultrafast decay component, τ(1), is attributed to the solvation dynamics of the particular solvent used. The second ultrafast decay component, τ(2), is found to vary from 50 to 215 ps depending upon the solvent. The amplitudes of the ultrafast decay components vary with the wavelength and show time dependent spectral shift in the emission maximum. The observation is interpreted that the time dependent spectral shift is not only due to solvation dynamics but also due to the existence of more than one emitting state of proflavine in the solvent used. Time resolved area normalized emission spectral (TRANES) analysis shows an isoemissive point, indicating the presence of two emitting states in homogeneous solution. Detailed femtosecond fluorescence decay analysis allows us to isolate the two independent emitting components of the close lying singlet states. The CIS and TDDFT calculations also support the existence of the close lying emitting states. The near constant lifetime observed for proflavine in different solvents is suggested to be due to the similar dipole moments of the ground and the evolved emitting singlet state of the dye from the Franck-Condon excited state.

Noncovalent interactions between hydroxylated polycyclic aromatic hydrocarbon and DNA: molecular docking and QSAR study

Polycyclic aromatic hydrocarbons (PAHs) can be hydroxylated by CYP450-oxidases (1A1 and 1B1 mainly) and may cause DNA damage and cancer. However, the mechanism of such interactions has not been fully understood. In this study, an integrated molecular docking and QSAR approach was employed to further investigate the binding interactions between hydroxylated PAHs (HO-PAHs) and calf thymus DNA (CT-DNA). Molecular docking, hydrogen-bonding, hydrophobic and π-π interactions were observed to be characteristic interactions between HO-PAHs and DNA. An optimum QSAR model with good robustness and predictability was developed based on the molecular structural parameters calculated by the density function theory and partial least squares. Additionally, the developed QSAR model indicated that the molecular size, polarizability and electrostatic potential of HO-PAHs were related to the binding affinities to DNA.