Chain length-dependent association of distamycin-type oligopeptides with A X T and G X C pairs in polydeoxynucleotide duplexes

Different binding affinities of various distamycin analogs including the deformylated derivative with poly(dA-dC) X poly(dG-dT) were investigated using CD measurements. The inhibitory effect of distamycins on the DNAase I cleavage activity of DNA duplexes strongly supports the binding data. The base specificity of the ligand interaction with duplex DNA depends on the chain length of distamycin analogs. Netropsin, distamycin-2 and the deformylated distamycin-3 show no binding to dG X dC containing sequences at moderate ionic strength and are classified as highly dA X dT specific. In contrast distamycin having three, four or five methylpyrrolecarboxamide groups also forms more or less stable complexes with dG X dC-containing duplexes. These ligands possess a lower basepair specificity. The correlation between binding behavior and oligopeptide structure shows that presence of the number of hydrogen acceptor and donor sites determines the basepair and sequence specificity. The additional interaction with dG X dC pairs becomes essential when the number of hydrogen acceptor sites exceeds n = 3.     

Chain length-dependent association of distamycin-type oligopeptides with A·T and G·C pairs in polydeoxynucleotide duplexes

Different binding affinities of various distamycin analogs including the deformylated derivative with poly(dA-dC)·poly(dG-dT) were investigated using CD measurements. The inhibitory effect of distamycins on the DNAase I cleavage activity of DNA duplexes strongly supports the binding data. The base specificity of the ligand interaction with duplex DNA depends on the chain length of distamycin analogs. Netropsin, distamycin-2 and the deformylated distamycin-3 show no binding to dG·dC containing sequences at moderate ionic strength and are classified as highly dA·dT specific. In contrast distamycin having three, four or five methylpyrrolecarboxamide groups also forms more or less stable complexes with dG·dC-containing duplexes. These ligands possess a lower basepair specificity. The correlation between binding behavior and oligopeptide structure shows that presence of the number of hydrogen acceptor and donor sites determines the basepair and sequence specificity. The additional interaction with dG·dC pairs becomes essential when the number of hydrogen acceptor sites exceeds n = 3.     

Z-DNA and other non-B-DNA structures are reversed to B-DNA by interaction with netropsin

The interaction between the B-form specific ligands netropsin (Nt) and distamycin-3 (Dst-3) and DNA duplexes has been studied under conditions of salt concentration and low water activity that modify the polymer conformation into a non-B DNA form, putatively a Z-like form. Three polymers with strict alternating purine-pyrimidine sequences and GC content from 100-0% have been tested: poly(dG-dC) . poly(dG-dC), poly(dA-dC) . poly(dG-dT) and poly(dA-dT) . poly(dA-dT). The titrations by Nt and Dst-3 were followed by circular dichroism. Although specific binding of Nt to the Z-form of poly(dG-dC) . poly(dG-dC) does not occur, Nt reverses this Z structure to the B-type conformation; Dst-3 is, however, totally inefficient. The presumed non-B or Z-like structure of poly(dA-dC) . poly(dG-dT) is reversed to the B-form upon interaction with Nt; Dst-3 also induces this reversal but at higher ligand ratios. The modified B-structure of poly(dA-dT) . poly(dA-dT) in low water activity is efficiently reversed to the B-form by interaction with both Nt and Dst-3.     

Binding of Non-Intercalating Antibiotics to B-DNA: A Theoretical Study Taking Into Account Nucleic Acid Flexibility

Abstract

A detailed theoretical study has been made for five antibiotics which all bind selectively to AT sequences in the minor groove of B-DNA: SN-18071, NSCT-101327, distamycin-2, distamycin-3 and netropsin. The optimal complexes were found for systems in which the flexibility of DNA, as well as that of the antibiotics, was taken into account. Explicit, mobile counterions and a dielectric function modelling aqueous solution were also included. The binding geometries of the most strongly interacting antibiotics, distamycin-3 and netropsin, are compared in considerable detail and it is shown that notable differences exist between them. The results for netropsin are also discussed in the light of recent disagreements concerning its exact binding location within DNA.     

Specific interaction of netropsin, distamycin-3 and analogs with LC duplexes: reversion towards the B form of the 2''-deoxy-.2''-deoxy-2''-fluoro-hybrid duplexes upon specific interaction with netropsin, distamycin-3 and analogs.

Binding of the B-form specific ligands netropsin and distamycin-3, -4 and -5 has been used to monitor the presence and/or the inducibility of a B-type structure in various poly-inosinic.poly-cytidilic double stranded polymers with deoxyribose, ribose or 2'-deoxy-2'-fluororibose as sugar on either strand. The efficiency of binding was followed by circular dichroism and further evaluated by the increase in melting temperature of the complexes. The efficient binding of netropsin and distamycins to the hybrid polymer (dIfl)n. (dC)n demonstrated that the fluorine carrying strand may undergo a A to B-type transition reflecting a change of the 2'-deoxy-2'-fluororibose from the 3'-endo to the 1'-exo or 2'-endo pucker. The less efficient binding of the same ligands to the reverse hybrid (dI)n.(dCfl)n showed that the geometry of the pyrimidine strand is the most critical for the specific interaction. Taking into account the recent findings about the regular hydration in the minor groove of the B-type dodecamer dCGCGAATTCGCG in solid-state, the different binding modes observed between the different polymers and antibiotics are explained by differences in their possibilities of hydration. Binding of netropsin to a double stranded deoxypolymer is interpreted as a local replacement of water molecules by netropsin in the minor groove hydration network which is typical of the B-form.     

Binding of non-intercalating antibiotics to B-DNA: a theoretical study taking into account nucleic acid flexibility

A detailed theoretical study has been made for five antibiotics which all bind selectively to AT sequences in the minor groove of B-DNA: SN-18071, NSC-101327, distamycin-2, distamycin-3 and netropsin. The optimal complexes were found for systems in which the flexibility of DNA, as well as that of the antibiotics, was taken into account. Explicit, mobile counterions and a dielectric function modelling aqueous solution were also included. The binding geometries of the most strongly interacting antibiotics, distamycin-3 and netropsin, are compared in considerable detail and it is shown that notable differences exist between them. The results for netropsin are also discussed in the light of recent disagreements concerning its exact binding location within DNA.     

A QM protein–ligand investigation of antipsychotic drugs with the dopamine D2 Receptor (D2R)

The dopamine D2 Receptor (D2R) is a member of the G-Protein-Coupled Receptor family and plays a critical role in neurotransmission activities in the human brain. Dysfunction in dopamine receptor signaling may lead to mental health illnesses such as schizophrenia and Parkinson’s disease. D2R is the target protein of the commonly used antipsychotic drugs such as risperidone, clozapine, aripiprazole, olanzapine, ziprasidone, and quetiapine. Due to their significant side effects and non-selective profiles, the discovery of novel drugs has become a challenge for researchers working in this field. Recently, our group has focused on the interactions of these drug molecules in the active site of the D2R using different in silico approaches. We here compare the performances of different approaches in estimating the drug binding affinities using quantum chemical approaches. Conformations of drug molecules (ligands) at the binding site of the D2R taken from the preliminary docking studies and molecular dynamics simulations were used to generate protein–ligand interaction models. In a first approach, the BSSE-corrected interaction energies of the ligands with the most critical amino acid Asp114 and with the other amino acids closest to ligands in the binding cavity were calculated separately by density functional theory method in implicit water environment at the M06-2X/6-31 g(d,p) level of the theory. In a second approach, ligand binding affinities were calculated by taking into consideration not only the interaction energies but also deformation and desolvation energies of ligands with surrounding amino acid residues, in a radius of 5 Å of the protein-bound ligand. The quantum mechanically obtained results were compared with the experimentally obtained binding affinity values. We concluded that although H-bond interactions of ligands with Asp114 are the most dominant interaction in the binding site, if van der Waals and steric interactions of ligands which have cumulative effect on the ligand binding are not included in the calculations, the interaction energies are overestimated.     

Study of effect of water on the interaction of DNA and actinocin derivatives with various aminoalkyl chain length by infrared spectrophotometry and computer modeling

A comparative study of the effect of water on the interaction of DNA with actinocin derivatives having different numbers of methylene groups in side chains was performed by IR spectroscopy. It was found that, as relative humidity increases, water molecules simultaneously bind to hydrate-active sites of DNA and ligands. The absorption band at v = 1137 cm-1, caused by oscillations of the C-O and P-O groups of atoms in the DNA-ligand complex having two methylene groups, is due to the interactions between the cationic groups of the ligand and the sugar-phosphate backbone of DNA, which may be one of the reasons for the high stability of this complex. Using computer simulation of interaction of DNA fragments and actinocin derivatives in water environment, molecular models of the formation of their complexes for two ways of binding were constructed.     

C60 binds to and deforms nucleotides

Atomistic molecular dynamics simulations are performed for up to 20 ns to monitor the formation and the stability of complexes composed of single- or double-strand DNA molecules and C60 in aqueous solution. Despite the hydrophobic nature of C60, our results show that fullerenes strongly bind to nucleotides. The binding energies are in the range -27 to -42 kcal/mol; by contrast, the binding energy of two fullerenes in aqueous solution is only -7.5 kcal/mol. We observe the displacement of water molecules from the region between the nucleotides and the fullerenes and we attribute the large favorable interaction energies to hydrophobic interactions. The features of the DNA-C60 complexes depend on the nature of the nucleotides: C60 binds to double-strand DNA, either at the hydrophobic ends or at the minor groove of the nucleotide. C60 binds to single-strand DNA and deforms the nucleotides significantly. Unexpectedly, when the double-strand DNA is in the A-form, fullerenes penetrate into the double helix from the end, form stable hybrids, and frustrate the hydrogen bonds between end-group basepairs in the nucleotide. When the DNA molecule is damaged (specifically, a gap was created by removing a piece of the nucleotide from one helix), fullerenes can stably occupy the damaged site. We speculate that this strong association may negatively impact the self-repairing process of the double-strand DNA. Our results clearly indicate that the association between C60 and DNA is stronger and more favorable than that between two C60 molecules in water. Therefore, our simulation results suggest that C60 molecules have potentially negative impact on the structure, stability, and biological functions of DNA molecules.     

Study of DNA base-Li doped SiC nanotubes in aqueous solutions: a computer simulation study

Due to the importance of soluble nanotubes in biological systems, computational research on DNA base functionalized nanotubes is of interest. This study presents the quantitative results of Monte Carlo simulations of Li-doped silicon carbide nanotubes and its nucleic acid base complexes in water. Each species was first modeled by quantum mechanical calculations and then Monte Carlo simulations were applied to study their properties in aqueous solution. Solvation free energies were computed to indicate the solvation behavior of these compounds. The computations show that solvation free energies of the complexes of DNA bases with Li-doped SiC nanotubes are in the order: thymine > cytosine > adenine > guanine. The results of complexation free energies were also used to study the stability of related structures, which indicate that thymine-Li-doped SiC nanotubes produce the most stable compound among the four DNA base complexes.     

Flexible docking using tabu search and an empirical estimate of binding affinity

This article describes the implementation of a new docking approach. The method uses a Tabu search methodology to dock flexibly ligand molecules into rigid receptor structures. It uses an empirical objective function with a small number of physically based terms derived from fitting experimental binding affinities for crystallographic complexes. This means that docking energies produced by the searching algorithm provide direct estimates of the binding affinities of the ligands. The method has been tested on 50 ligand-receptor complexes for which the experimental binding affinity and binding geometry are known. All water molecules are removed from the structures and ligand molecules are minimized in vacuo before docking. The lowest energy geometry produced by the docking protocol is within 1.5 Å root-mean square of the experimental binding mode for 86% of the complexes. The lowest energies produced by the docking are in fair agreement with the known free energies of binding for the ligands. Proteins 33:367–382, 1998. © 1998 Wiley-Liss, Inc.     

Oxidative DNA cleavage by Schiff base tetraazamacrocyclic oxamido nickel(II) complexes

Nickel is considered a weak carcinogen. Some researches have shown that bound proteins or synthetic ligands may increase the toxic effect of nickel ions. A systematic study of ligand effects on the interaction between nickel complexes and DNA is necessary. Here, we compared the interactions between DNA and six closely related Schiff base tetraazamacrocyclic oxamido nickel(II) complexes NiL(1-3a,1-3b). The structure of one of the six complexes, NiL(3b) has been characterized by single crystal X-ray analysis. All of the complexes can cleave plasmid DNA under physiological conditions in the presence of H(2)O(2). NiL(3b) shows the highest DNA cleavage activity. It can convert supercoiled DNA to nicked DNA then linear DNA in a sequential manner as the complex concentration or reaction time is increased. The cleavage reaction is a typical pseudo-first-order consecutive reaction with the rate constants of 3.27+/-0.14h(-1) (k(1)) and 0.0966+/-0.0042h(-1) (k(2)), respectively, when a complex concentration of 0.6mM is used. The cleavage mechanism between the complex and plasmid DNA is likely to involve hydroxyl radicals as reactive oxygen species. Circular dichronism (CD), fluorescence spectroscopy and gel electrophoresis indicate that the complexes bind to DNA by partial intercalative and groove binding modes, but these binding interactions are not the dominant factor in determining the DNA cleavage abilities of the complexes.     

The V(D)J recombination signal sequence and kappa B binding protein Rc binds DNA as dimers and forms multimeric structures with its DNA ligands.

The murine DNA binding protein Rc binds to the heptamer motif of the V(D)J recombination signal sequences and to the kappa B motif of the immunoglobulin enhancer. Bacterial fusion proteins for Rc and DNA ligands of Rc form multiple protein-DNA complexes in electrophoretic mobility shift assays (EMSA). Large complexes formation is favored by an increased Rc concentration. In order to determine the architecture of these complexes, the apparent molecular weights of the protein-DNA complexes were first determined by their gel mobilities. The data suggest that Rc binds to its DNA ligands as dimers, tetramers, and multiples of tetramers. The inference that Rc binds DNA as dimers was substantiated by the formation of chimeric complexes when two electrophoretically distinguishable Rc proteins were employed in EMSA. Methylation interference experiments show that there are no contiguous protein binding sites evident in the DNA of the larger complexes. Apparently, multimerization occurs via protein-protein interactions. Such interaction was demonstrated by the formation of Rc dimers and tetramers in a chemical crosslinking experiment. Significantly, the multimerization of DNA-bound Rc could be involved in bringing the variable region gene segments together for the somatic V(D)J recombination.     

Theoretical study on binding of Hoechst 33258 with oligonucleotides

Computer modelling with an energy minimization procedure is used here to obtain stereochemical and energetic details for complexes of the dye Hoechst 33258 with different oligonucleotide sequences. An optimised model of the dye with d(A)5 X d(T)5 is in conformity with previous proposed models. It has bifurcated hydrogen bonds between N2H and N4H of benzimidazole rings with N3 of adenine and O2 of thymine. Relative binding energies with different oligonucleotides show preference for AT containing sequences, with an intermediate affinity between that for netropsin and distamycin-2. Reduced binding is observed at high ionic concentration. The benzimidazole rings are twisted with respect to the phenol ring in the optimal model. This gives desired curvature to the molecule which is stabilised by intermolecular forces.     

Binding of nonintercalative antitumor drugs to DNA-polymers: structural effects of bisquaternary ammonium heterocycles

The binding of the antitumor agents SN-16814 nd SN-13232 to various DNA's in solution was monitored by CD and UV absorption measurements. In addition comparative studies with dA.dT containing duplex DNA of the related ligands SN-6136 and SN-6324 were included with respect to effects of structural variations. In general all four ligands show a dA.dT preference in their binding affinity to DNA. Differences were observed for the reaction of SN-16814 which contains bicyclic ring system: it has a lower base pair selectivity, shows some affinity to poly(dG-dC).poly(dG-dC), poly(rA).poly(rU) and poly(rU). The binding mechanism of SN-16814 is associated with a significant time dependent binding effect in CD spectra and UV absorption in case of reaction with poly(dA).poly(dT) and poly(dI).poly(dC) indicating a slow kinetics. The preferred binding to dA.dT base pairs in DNA decreases in the order from SN-61367 greater than SN-13232 greater than SN-6324,SN-16814 as judged from CD titration studies, salt dissociation and melting temperature data. Competitive binding experiments with netropsin (Nt) or distamycin-5 revealed that SN-16814 and SN-13232 are displaced from poly(dA.dT).poly(dA-dT) suggesting that both ligands are less strongly bound than Nt and Dst-5 within the minor groove of B-DNA. These studies are consistent with results of the DNAse I cleavage of poly(dA-dT).poly(dA-dT) which show the same relative order of inhibition of the cleavage reaction due to ligand binding. The results suggest that the variability of the DNA binding and dA.dT sequence specificity may reside in the adaptability of benzamide-type ligands in the helical groove which is influenced by distinct structural modifications of the ligand conformation.     

Carboxylate binding modes in zinc proteins: A theoretical study

The relative energies of different coordination modes (bidentate, monodentate, syn, and anti) of a carboxylate group bound to a zinc ion have been studied by the density functional method B3LYP with large basis sets on realistic models of the active site of several zinc proteins. In positively charged four-coordinate complexes, the mono- and bidentate coordination modes have almost the same energy (within 10 kJ/mol). However, if there are negatively charged ligands other than the carboxylate group, the monodentate binding mode is favored. In general, the energy difference between monodentate and bidentate coordination is small, 4-24 kJ/mol, and it is determined more by hydrogen-bond interactions with other ligands or second-sphere groups than by the zinc-carboxylate interaction. Similarly, the activation energy for the conversion between the two coordination modes is small, approximately 6 kJ/mol, indicating a very flat Zn-O potential surface. The energy difference between syn and anti binding modes of the monodentate carboxylate group is larger, 70-100 kJ/mol, but this figure again strongly depends on interactions with second-sphere molecules. Our results also indicate that the pK(a) of the zinc-bound water ligand in carboxypeptidase and thermolysin is 8-9.     

Combined theoretical studies on solvation and hydrogen bond interactions in glycine tripeptide

The theoretical conformational analysis of glycine tripeptide (GT) has been carried out by molecular dynamics (MD) method in order to find minimum energy conformations. The MD studies on GT with water have been carried out for over 10 ns with a time step of 2 fs using fixed charge force field (AMBER ff03). By adding the solvation effect using water as a solvent, the GT conformers identified in this study exhibit α-helical conformation. Compared with the earlier reports, this MD study is able to identify the energetically favourable GT conformations. The obtained geometry of the five most stable GT conformations was optimised using the density functional theory method at B3LYP/6-311G** level of theory. Subsequently, the effects of solvation on the conformational characteristics of five most stable GT conformers with four water molecules (the number of water molecules in the first solvation shell of GT obtained from MD study) were investigated using the same method and the same level of theory. The effect of microsolvation on the fifth GT conformer has been studied with a cluster of 11 water molecules as the first hydration shell which generates folded structure. The interaction energies of all the complexes are calculated by correcting the basis set superposition error. The strong hydrogen bond mainly contributes to the interaction energies. The atoms in molecules theory and natural bond orbital analysis were used to study the origin of H-bonds. A good correlation between the structural parameters and the properties of charge density is found. NMR calculations show that the C = O carbons of the amine groups of the first and middle glycine fragments have maximum chemical shifts.     

Studies on the interaction of extended terpyridyl and triazine metal complexes with DNA

The interaction with duplex DNA of a series of 12 previously synthesised ruthenium (II) complexes incorporating extended terpyridines or analogous ligands is reported. The interaction was investigated using absorption spectroscopy and also, when possible, luminescence spectroscopy. Binding parameters for the interactions were estimated using established models. It was found that binding affinities are relatively large and their magnitude is somewhat dependent on the metal centre attached to the ligand. Viscosity studies reveal that all the complexes interact with DNA through groove binding and not intercalation. The variation in binding affinities of these complexes is rationalised by a consideration of electrostatic and steric factors.     

Hydration effects accompanying the formation of DNA complexes with some ligands

In the range of millimeter wavelengths the dielectric properties of aqueous solutions of some biologically active ligands (potential anticarcinogen chlorophyllin; pharmacological drug caffeine; polyamine putrescine; mutagens proflavine and ethidium bromide; actinocin derivative, an analogue of antitumor antibiotic actinomycin D) and DNA complexes with these substances were studied. It was shown that complex formation is accompanied by the change in dielectric properties of the solution. These changes during interaction of DNA with the first three compounds correspond to a decrease in hydration (compared with the total hydration of free components), and in other cases they cause an increase in hydration. The number of water molecules bound with both the ligand and DNA nucleotide in the complex was estimated. The results were compared with existing models of DNA interaction with the studied substances.