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.     

Binding free energy of selected anticancer compounds to DNA--theoretical calculations

Many studies have elucidated structures and thermodynamics of complexes formed by different ligands with DNA. However, in most cases structural and free energy binding studies were not correlated with each other because of the problem of identifying which experimental free energy of binding corresponds to which experimental DNA-ligand structure. In the present work, Poisson-Boltzmann and solvent-accessible surface area methods were used to predict unknown modes of interaction between DNA and three different ligands: mitoxantrone and two pyrimidoacridine derivatives. In parallel, experimental measurements of binding free energy for the studied complexes were performed to compare experimental and calculated values. Our studies showed that the calculated values of free energy are only close to experimental data for some models of interaction between ligands and DNA. Based on this correlation, the most likely models of DNA-ligand complexes were postulated: (i) mitoxantrone and one derivative of pyrimidoacridine, both with two charged side chains, intercalate from the minor groove of DNA and bind with both chains in this groove; (ii) pyrimidoacridine, with only one side chain, very likely does not intercalate into DNA at all. Additionally, the non-electrostatic and electrostatic parts of the calculated binding free energy for the DNA-ligands studied are discussed.     

Interaction of DNA with benzocrown derivatives of actinocin

Complexes of DNA with actinocin derivatives containing benzocrown groups at the 1 and/or 9 positions of the chromophore were studied by spectrophotometric titration and circular dichroism. The actinocin chromophore and the crown fragments are the binding sites of the ligands with DNA. The mode of ligand-DNA binding is shown to depend on the size of the crown group, its distance to the actinocin chromophore, and the ionic strength of the medium. Selective toward Na+ ion benzocrown fragments combine with DNA phosphate groups. The simultaneous interaction of the actinocin chromophore with the DNA bases is possible only at optimal distance between both the binding sites of ligand molecule.     

Interaction of DNA with Benzocrown Derivatives of Actinocin

Complexes of DNA with actinocin derivatives containing benzocrown groups at positions 1 and/or 9 of the chromophore were studied by spectrophotometric titration and circular dichroism. The actinocin chromophore and the crown fragments are the DNA-binding sites of the ligands. The mode of ligand–DNA binding is shown to depend on the size of the crown group, its distance to the actinocin chromophore, and the ionic strength of the medium. Na⁺-selective benzocrown fragments combine with DNA phosphate groups. The simultaneous interaction of the actinocin chromophore with the DNA bases is possible only at an optimal distance between the two binding sites of ligand molecule.     

Molecular dynamics of anthraquinone DNA intercalators with polyethylene glycol side chains

The interactions of a homologous series of four anthraquinone (AQ) intercalators with increasing lengths of polyethylene glycol (PEG) side chains with DNA have been studied via molecular dynamics (MD) simulations. The geometry, conformation, interactions, and hydration of the complexes were examined. The geometries of the four ligands were similar with parallel stacking to the long axis of the adjacent DNA base pairs. Hydrogen bonding between the AQ amide and DNA led to a preference for the trans-syn conformer. As the side chain lengthened, binding to DNA reduced the conformational space, resulting in an increase in unfavorable entropy. Increased localization of the PEG side chain in the DNA groove, indicating some interaction of the side chain with DNA, also contributed unfavorably to the entropy. The changes in free energy of binding due to entropic considerations (-3.9 to -6.3 kcal/mol) of AQ I-IV were significant. The hydration of the PEG side chain decreased upon binding to DNA. Understanding of side chain conformations, interactions, and hydration changes that accompany the formation of a ligand-DNA complex may be important in the development of new applications of pegylated small molecules that target biological macromolecules.     

Structural basis for the binding affinity of a homologous series of synthetic phenoxazone drugs with DNA: NMR and molecular mechanics analysis

The molecular basis of the marked structure-activity relationship for a homologous series of DNA-binding phenoxazone drugs (ActII-ActIV) has been investigated by NMR spectroscopy and molecular mechanics. The spatial structures of the complexes between the drugs and a model deoxytetranucleotide, 5'-d(TpGpCpA), have been determined by molecular mechanics methods using homonuclear (1)H-(1)H 2D-NOESY and heteronuclear (1)H-(31)P (HMBC) NMR spectroscopic data. Observed intermolecular NOE contacts and equilibrium binding studies confirm that the binding affinity of the synthetic phenoxazone derivatives with d(TGCA) decreases with an increase in the number of CH(2) groups in the dimethylaminoalkyl side chains, i.e., ActII > ActIII > ActIV, in agreement with the observed biological activity of these compounds. Molecular mechanics calculations of the spatial structures of the intercalated complexes of ActII-ActIV with d(TGCA) indicate that the different binding constants of the phenoxazone derivatives with the DNA oligomer are due to the different degrees of intercalation of the chromophore and the different steric arrangements of aminoalkyl side chains in the minor groove of the tetramer duplex; this results in different distances between the negatively-charged phosphates of the DNA duplex and the terminal positively-charged N(CH(3))(2) groups of the side chains.     

Molecular Dynamics of Anthraquinone DNA Intercalators with Polyethylene Glycol Side Chains

Abstract

The interactions of a homologous series of four anthraquinone (AQ) intercalators with increasing lengths of polyethylene glycol (PEG) side chains with DNA have been studied via molecular dynamics (MD) simulations. The geometry, conformation, interactions, and hydration of the complexes were examined. The geometries of the four ligands were similar with parallel stacking to the long axis of the adjacent DNA base pairs. Hydrogen bonding between the AQ amide and DNA led to a preference for the trans-syn conformer. As the side chain lengthened, binding to DNA reduced the conformational space, resulting in an increase in unfavorable entropy. Increased localization of the PEG side chain in the DNA groove, indicating some interaction of the side chain with DNA, also contributed unfavorably to the entropy. The changes in free energy of binding due to entropic considerations (— 3.9 to—6.3 kcal/mol) of AQ I-IV were significant. The hydration of the PEG side chain decreased upon binding to DNA. Understanding of side chain conformations, interactions, and hydration changes that accompany the formation of a ligand-DNA complex may be important in the development of new applications of pegylated small molecules that target biological macromolecules.     

Is There Any Interaction Between Telomeric DNA Structures,G-Quadruplex and I-Motif,with Saffron Active Metabolites?

Telomeric DNA contains some unique secondary structures, such as G-quadruplex and I-motif. These structures may be stabilized or changed by binding to specific proteins or small molecules. Herein, we report the in vitro effect of crocin, crocetin, picrocrocin, and safranal on these structures. Circular dichroism (CD) data indicate that crocetin has higher affinity for these structures. Safranal and crocin induce little change in the I-motif and G-quadruplex, respectively. The molecular docking confirms the experimental data and indicates the minor groove binding of ligands with G-quadruplex. The possibility for application of these ligands as sequence-specific drugs should be further investigated.     

The interaction of ellipticine derivatives with nucleic acids studied by optical and 1H-nmr spectroscopy: Effect of size of the heterocyclic ring system

The DNA interaction of derivatives of ellipticine with heterocyclic ring systems with three, four, or five rings and a dimethylaminoethyl side chain was studied. Optical spectroscopy of drug complexes with calf thymus DNA, poly [(dA-dT) · (dA-dT)], or poly [(dG-dC) · (dG-dC)] showed a 10 nm bathochromic shift of the light absorption bands of the pentacyclic and tetracyclic compounds upon binding to the nucleic acids, which indicates binding by intercalation. For the tricyclic compound a smaller shift of 1–3 nm was observed upon binding to the nucleic acids. Flow linear dichroism studies show that the geometry of all complexes is consistent with intercalation of the ring system, except for the DNA and poly [(dG-dC) · (dG-dC)] complexes of the tricyclic compound, where the average angle between the drug molecular plane and the DNA helix axis was found to be 65°. One-dimensional ¹H-nmr spectroscopy was used to study complexes between d(CGCGATCGCG)₂ and the tricyclic and pentacyclic compounds. The results on the pentacyclic compound show nonselective broadening due to intermediate chemical exchange of most oligonucleotide resonances upon drug binding. The imino proton resonances are in slow chemical exchange, and new resonances with upfield shifts approaching 1 ppm appear upon drug binding, which supports intercalative binding of the pentacyclic compound. The results on the tricyclic compound show more rapid binding kinetics and very selective broadening of resonances. The data suggest that the tricyclic compound is in an equilibrium between intercalation and minor groove binding, with a preference to bind close to the AT base pairs with the side chain residing in the minor groove. © 1994 John Wiley & Sons, Inc.     

Surfactant–copper(II) Schiff base complexes: synthesis,structural investigation,DNA interaction,docking studies,and cytotoxic activity

A series of surfactant–copper(II) Schiff base complexes (1–6) of the general formula, [Cu(sal-R₂)₂] and [Cu(5-OMe-sal-R₂)₂], {where, sal?=?salicylaldehyde, 5-OMe-sal?=?5-methoxy- salicylaldehyde, and R₂?=?dodecylamine (DA), tetradecylamine (TA), or cetylamine (CA)} have been synthesized and characterized by spectroscopic, ESI-MS, and elemental analysis methods. For a special reason, the structure of one of the complexes (2) was resolved by single crystal X-ray diffraction analysis and it indicates the presence of a distorted square-planar geometry in the complex. Analysis of the binding of these complexes with DNA has been carried out adapting UV-visible-, fluorescence-, as well as circular dichroism spectroscopic methods and viscosity experiments. The results indicate that the complexes bind via minor groove mode involving the hydrophobic surfactant chain. Increase in the length of the aliphatic chain of the ligands facilitates the binding. Further, molecular docking calculations have been performed to understand the nature as well as order of binding of these complexes with DNA. This docking analysis also suggested that the complexes interact with DNA through the alkyl chain present in the Schiff base ligands via the minor groove. In addition, the cytotoxic property of the surfactant–copper(II) Schiff base complexes have been studied against a breast cancer cell line. All six complexes reduced the visibility of the cells but complexes 2, 3, 5, and 6 brought about this effect at fairly low concentrations. Analyzed further, but a small percentage of cells succumbed to necrosis. Of these complexes (6) proved to be the most efficient aptotoxic agent.     

Binding of netropsin to several DNA constructs: evidence for at least two different 1:1 complexes formed from an -AATT-containing ds-DNA construct and a single minor groove binding ligand

Isothermal titration calorimetry, ITC, has been used to determine the thermodynamics (DeltaG, DeltaH, and -TDeltaS) for binding netropsin to a number of DNA constructs. The DNA constructs included: six different 20-22mer hairpin forming sequences and an 8-mer DNA forming a duplex dimer. All DNA constructs had a single -AT-rich netropsin binding with one of the following sequences, (A(2)T(2))(2), (ATAT)(2), or (AAAA/TTTT). Binding energetics are less dependent on site sequence than on changes in the neighboring single stranded DNA (hairpin loop size and tail length). All of the 1:1 complexes exhibit an enthalpy change that is dependent on the fractional saturation of the binding site. Later binding ligands interact with a significantly more favorable enthalpy change (partial differential DeltaH(1-2) from 2 to 6 kcal/mol) and a significantly less favorable entropy change (partial differential (-TDeltaS(1-2))) from -4 to -9 kcal/mol). The ITC data could only be fit within expected experimental error by use of a thermodynamic model that includes two independent binding processes with a combined stoichiometry of 1 mol of ligand per 1 mol of oligonucleotide. Based on the biophysical evidence reported here, including theoretical calculations for the energetics of "trapping" or structuring of a single water molecule and molecular docking computations, it is proposed that there are two modes by which flexible ligands can bind in the minor groove of duplex DNA. The higher affinity binding mode is for netropsin to lay along the floor of the minor groove in a bent conformation and exclude all water from the groove. The slightly weaker binding mode is for the netropsin molecule to have a slightly more linear conformation and for the required curvature to be the result of a water molecule that bridges between the floor of the minor groove and two of the amidino nitrogens located at one end of the bound netropsin molecule.     

Distamycin-induced inhibition of homeodomain-DNA complexes.

The mobility shift assay was used to study the competition of the minor groove binder distamycin A with either an Antennapedia homeodomain (Antp HD) peptide or derivatives of a fushi tarazu homeodomain (ftz HD) peptide for their AT-rich DNA binding site. The results show that distamycin and the homeodomain peptides compete under the conditions: (i) preincubation of DNA with distamycin and subsequent addition of HD peptide; (ii) simultaneous incubation of DNA with distamycin and HD peptide; and (iii) preincubation of DNA with HD peptide and subsequent addition of distamycin. There is also competition when using a peptide which lacks the N-terminal arm of ftz HD that is involved in contacts in the minor groove. It is proposed that the protein's binding affinity is diminished by distamycin-induced conformational changes of the DNA. The feasibility of the propagation of conformational changes upon binding in the minor groove is also shown for the inhibition of restriction endonucleases differing in the AT content of their recognition site and of their flanking DNA sequences. Thus, it is demonstrated that minor groove binders can compete with the binding of proteins in the major groove, providing an experimental indication for the influence of biological activities exerted by DNA ligands binding in the minor groove.     

Side-chain flexibility in protein-ligand binding: the minimal rotation hypothesis

The goal of this work is to learn from nature about the magnitudes of side-chain motions that occur when proteins bind small organic molecules, and model these motions to improve the prediction of protein-ligand complexes. Following analysis of protein side-chain motions upon ligand binding in 63 complexes, we tested the ability of the docking tool SLIDE to model these motions without being restricted to rotameric transitions or deciding which side chains should be considered as flexible. The model tested is that side-chain conformational changes involving more atoms or larger rotations are likely to be more costly and less prevalent than small motions due to energy barriers between rotamers and the potential of large motions to cause new steric clashes. Accordingly, SLIDE adjusts the protein and ligand side groups as little as necessary to achieve steric complementarity. We tested the hypothesis that small motions are sufficient to achieve good dockings using 63 ligands and the apo structures of 20 different proteins and compared SLIDE side-chain rotations to those experimentally observed. None of these proteins undergoes major main-chain conformational change upon ligand binding, ensuring that side-chain flexibility modeling is not required to compensate for main-chain motions. Although more frugal in the number of side-chain rotations performed, this model substantially mimics the experimentally observed motions. Most side chains do not shift to a new rotamer, and small motions are both necessary and sufficient to predict the correct binding orientation and most protein-ligand interactions for the 20 proteins analyzed.     

Ligand‐based pharmacophore modelling and screening of DNA minor groove binders targeting Staphylococcus aureus

The recognition of DNA by small molecules is of special importance in the design of new drugs. Many natural and synthetic compounds have the ability to interact with the minor groove of DNA. In the present study, identification of minor groove binding compounds was attained by the combined approach of pharmacophore modelling, virtual screening and molecular dynamics approach. Experimentally reported 32 minor groove binding compounds were used to develop the pharmacophore model. Based on the fitness score, best three pharmacophore hypotheses were selected and used as template for screening the compounds from drug bank database. This pharmacophore‐based screening provides many compounds with the same pharmacological properties. All these compounds were subjected to four phases of docking protocols with combined Glide‐quantum‐polarized ligand docking approach. Molecular dynamics results indicated that selected compounds are more active and showed good interaction in the binding site of DNA. Based on the scoring parameters and energy values, the best compounds were selected, and antibacterial activity of these compounds was identified using in vitro antimicrobial techniques. Copyright © 2014 John Wiley & Sons, Ltd.     

Structural Basis for the Binding Affinity of a Homologous Series of Synthetic Phenoxazone Drugs with DNA: NMR and Molecular Mechanics Analysis

Abstract

The molecular basis of the marked structure-activity relationship for a homologous series of DNA-binding phenoxazone drugs (ActII-ActIV) has been investigated by NMR spectroscopy and molecular mechanics. The spatial structures of the complexes between the drugs and a model deoxytetranucleotide, 5′-d(TpGpCpA), have been determined by molecular mechanics methods using homonuclear ¹H-¹H 2D-NOESY and heteronuclear ¹H-³¹P (HMBC) NMR spectroscopic data. Observed intermolecular NOE contacts and equilibrium binding studies confirm that the binding affinity of the synthetic phenoxazone derivatives with d(TGCA) decreases with an increase in the number of CH₂ groups in the dimethylami- noalkyl side chains, i.e., ActII > ActIII > ActIV, in agreement with the observed biological activity of these compounds. Molecular mechanics calculations of the spatial structures of the intercalated complexes of ActII-ActIV with d(TGCA) indicate that the different binding constants of the phenoxazone derivatives with the DNA oligomer are due to the different degrees of intercalation of the chromophore and the different steric arrangements of aminoalkyl side chains in the minor groove of the tetramer duplex; this results in different distances between the negatively-charged phosphates of the DNA duplex and the terminal positively-charged N(CH₃)₂ groups of the side chains.     

Interaction of DNA with actinocin amides containing cationoid centers in the amide groups

Complexes of DNA with actinocin derivatives containing ω-dialkylaminoalkyl groups in the1 and/or 9 positions of the chromophore were studied by spectrophotometric titration, circular dichroism, and viscometry. Induced circular dichroism (ICD) spectra of the DNA-ligand complexes were compared for the cases of the complexes of known structure established by other methods. It was shown that the presence of an isoelliptic point in the long-wavelength absorption band of the ICD spectra of the ligand under monomeric binding conditions could indicate intercalation of the actinocin chromophore into DNA. The separation of the cationoid center and the chromosphore upon elongation of the methylene chain increases the aggregability of the ligand pn the surface of the DNA double helix, which prevents the intercalation of the chromophore.     

Structural studies on ligand–DNA systems: A robust approach in drug design

Molecular docking, molecular mechanics, molecular dynamics and relaxation matrix simulation protocols have been extensively used to generate the structural details of ligand-receptor complexes in order to understand the binding interactions between the two entities. Experimental methods like NMR spectroscopy and X-ray crystallography are known to provide structural information about ligand-receptor complexes. In addition, fluorescence spectroscopy, circular dichroism (CD) spectroscopy and molecular docking have also been utilized to decode the phenomenon of the ligand-DNA interactions, with good correlation between experimental and computational results. The DNA binding affinity was demonstrated by analysing fluorescence spectral data. Structural rigidity of DNA upon ligand binding was identified by CD spectroscopy. Docking is carried out using the DNA-Dock program which results in the binding affinity data along with structural information like interatomic distances and H-bonding, etc. The complete structural analyses of various drug-DNA complexes have afforded results that indicate a specific DNA binding pattern of these ligands. It also exhibited that certain structural features of ligands can make a ligand to be AT- or GC-specific. It was also demonstrated that changing specificity from AT base pairs to GC base pairs further improved the DNA topoisomerase inhibiting activity in certain ligands. Thus, a specific molecular recognition signature encrypted in the structure of ligand can be decoded and can be effectively employed in designing more potent antiviral and antitumour agents.     

The Influence of Sequence Microvariation Among HLA-DR3 Alleles on their Structure and Complexes Formation with Presentation Peptides

The three-dimensional structure of human leukocyte antigens HLA-DR*0301 and HLA-DR*0302 have been calculated using the homology modeling approach. General structural features of our models are similar to those of related HLA molecules. The typical layout of segments of the secondary structure is well preserved. However polypeptide chains are less tightly bound, which causes slightly broader opening of the binding groove. It also results in the modified layout of pockets in the binding groove. Amino acids defining the restricted sequence diversity of studied proteins are easily available for interactions with ligands.A set of docking simulations was performed using modeled structures of both HLA molecules and various specific peptide ligands. The control docking of influenza hemagglutinin peptide into HLA-DR*0101 molecule gives the complex structure which is in good agreement with that from crystallographic studies. The extensive analysis of the structure of modeled complexes of HLA-DR*0301 and HLA-DR*0302 with various ligands indicates that sequence microvariation of both alleles is not directly controlling the binding specificity. Preferences for binding of specific ligands, as evaluated from interactions in modeled complexes, agree qualitatively with experimental observations. Thus the computer aided docking simulations can be successfully used to calculate the three-dimensional structure of HLA-ligand complexes. However detailed explanation of binding specificity can not be achieved using presently available modeling procedures.Electronic Supplementary Material available.     

Design of Sequence-Specific DNA Binding Ligands That Use a Two-Stranded Peptide Motif for DNA Sequence Recognition

Abstract

The design and DNA binding activity of β-structure-forming peptides and netropsin-peptide conjugates are reported. It is found that a pair of peptides - S,S'-bis(Lys-Gly-Val-Cys-Val- NH-NH-Dns) - bridged by an S-S bond binds at least 10 times more strongly to poly(dG)?poly(dC) than to poly(dA)?poly(dT). This peptide can also discriminate between 5′-GpG-3′ and 5′-GpC-3′ steps in the DNA minor groove. Based on these observations, new synthetic ligands, bis-netropsins, were constructed in which two netropsin-like fragments were attached by means of short linkers to a pair of peptides - Gly-Cys-Gly- or Val-Cys-Val - bridged by S-S bonds. These compounds possess a composite binding specificity: the peptide chains recognize 5′-GpG-3′ steps on DNA, whereas the netropsin-like fragments bind preferentially to tuns of 4 AT base pairs. Our data indicate that combining the AT-base-pair specific properties of the netropsin-type structure with the 5′-GpG-3′-specific properties of certain oligopeptides offers a new approach to the synthesis of ligands capable of recognizing mixed sequences of AT- and GC-base pairs in the DNA minor groove. These compounds are potential models for DNA-binding domains in proteins which specifically recognize base pair sequences in the minor groove of DNA.