Molecular dynamics simulations and binding free energy analysis of DNA minor groove complexes of curcumin

Curcumin is a natural phytochemical that exhibits a wide range of pharmacological properties, including antitumor and anticancer activities. The similarity in the shape of curcumin to DNA minor groove binding drugs is the motivation for exploring its binding affinity in the minor grooves of DNA sequences. Interactions of curcumin with DNA have not been extensively examined, while its pharmacological activities have been studied and documented in depth. Curcumin was docked with two DNA duplexes, d(GTATATAC)₂ and d(CGCGATATCGCG)₂, and molecular dynamics simulations of the complexes were performed in explicit solvent to determine the stability of the binding. In all systems, the curcumin is positioned in the minor groove in the A·T region, and was stably bound throughout the simulation, causing only minor modifications to the structural parameters of DNA. Water molecules were found to contribute to the stability of the binding of the ligand. Free energy analyses of the complexes were performed with MM-PBSA, and the binding affinities that were calculated are comparable to the values reported for other similar nucleic acid–ligand systems, indicating that curcumin is a suitable natural molecule for the development of minor groove binding drugs.     

DNA binding of the nonintercalative ligands SN-6132, SN-6131 and SN-6113: minor variations of the ligand structure may cause changes in the base pair preference.

The DNA binding selectivity of three ligands of a series of antitumor agents of bisquaternary ammonium heterocycles has been investigated by means of CD spectroscopy and melting measurements. From the spectroscopic results and binding data it is concluded that the agents SN-6132, SN-6131 and SN-6113 have relatively high affinity to AT base pair sequences whereas the binding to GC pairs is very low. The binding selectivity to AT base pair sequences decreases in the order netropsin > SN-6132 > SN-6113 > SN-6131. Poly(dA).poly(dT) has the highest binding preference for SN-6132 relative to that of SN-6131. The different binding behavior of the ligands is related to their distinct changes in the chemical structure and to the DNA minor groove properties which determines the adaptability of the ligands in the groove.     

Molecular docking study investigating the possible mode of binding of C.I. Acid Red 73 with DNA

C.I. Acid Red 73 is a reactive azo dye with a variable potential carcinogenicity. The mechanism mediating interactions that occur between the dye and DNA have not been completely understood thus far. In this study, molecular docking techniques were applied to describe the most probable mode of DNA binding as well as the sequence selectivity of the C.I. Acid Red 73 dye. These docking experiments revealed that the dye is capable of interacting with the minor groove of the DNA on the basis of its curved shape, which fits well with the topology of double-stranded DNA. In addition, the dye can bind selectively to the minor groove of the DNA by applying CGT sequence selectivity. Further, the minor groove can be recognized although DNA targets present intercalation gaps. However, intercalative binding can also occur when the DNA target possesses an appropriate intercalation gap. Compared with the other eight DNA sequences that were studied, the DNA dodecamer d(CGCGATATCGCG)₂ (PDB ID: 1DNE) presents a very favorable target for the binding of C.I. Acid Red 73 to the minor groove, with the lowest binding free energy −9.19 kcal/mol. Results reported from this study are expected to provide useful information for research involving further simulations of molecular dynamics and toxicology investigations of the dye.     

Molecular dynamics simulations and free energy calculations of netropsin and distamycin binding to an AAAAA DNA binding site

Molecular dynamics simulations have been performed on netropsin in two different charge states and on distamycin binding to the minor groove of the DNA duplex d(CGCGAAAAACGCG)·d(CGCGTTTTTCGCG). The relative free energy of binding of the two non-covalently interacting ligands was calculated using the thermodynamic integration method and reflects the experimental result. From 2 ns simulations of the ligands free in solution and when bound to DNA, the mobility and the hydrogen-bonding patterns of the ligands were studied, as well as their hydration. It is shown that even though distamycin is less hydrated than netropsin, the loss of ligand–solvent interactions is very similar for both ligands. The relative mobilities of the ligands in their bound and free forms indicate a larger entropic penalty for distamycin when binding to the minor groove compared with netropsin, partially explaining the lower binding affinity of the distamycin molecule. The detailed structural and energetic insights obtained from the molecular dynamics simulations allow for a better understanding of the factors determining ligand–DNA binding.     

Spectroscopic and calorimetric studies on the DNA recognition of pyrrolo[2,1-c][1,4]benzodiazepine hybrids

DNA binding of two hybrid ligands composed of an alkylating pyrrolo[2,1-c][1,4]benzodiazepine (PBD) moiety tethered to either a naphthalimide or a phenyl benzimidazole chromophore was studied by DNA melting experiments, UV and fluorescence titrations, CD spectroscopy and isothermal titration calorimetry (ITC). Binding of both hybrids results in a remarkable thermal stabilization with an increase of DNA melting temperatures by up to 40 °C for duplexes that allow for a covalent attachment of the PBD moiety to guanine bases in their minor groove. CD spectroscopic measurements suggest that the naphthalimide moiety of the drug interacts through intercalation. In contrast, the PBD-benzimidazole hybrid binds in the DNA minor groove with a preference for (A,T)₄G sequences. Whereas the binding of both ligands is enthalpy-driven and associated with a negative entropy, the benzimidazole hybrid exhibits a less favourable binding enthalpy that is counterbalanced by a more favourable entropic term when compared to the naphthalimide hybrid.     

Drug Recognition of DNA. Proposal for GC Minor Groove Specific Ligands: Vinylexins

Abstract

In a previous publication in this journal (ref. 14) we have proposed an isolexin-like prototype of a GC minor groove specific ligand. The present paper is devoted to refinements of this prototype (increase in specificity and in DNA binding energy). It is shown that only a very limited improvement can be obtained by increasing the proton accepting capabilities of the heteroaromatic ring systems of the prototype, although these rings interact directly with the proton donating NH₂ group of guanine. On the other hand a significant increase both in GC specificity and in DNAbinding energy is obtained by replacing the NH linkers of the isolexin by C=C double bonds (yielding what we term “vinylexins”). Specificity is still largely conserved and the DNA binding energy is significantly increased in monocationic vinylexins, which should thus be efficient GC minor groove specific ligands. The outstanding importance for the GC specificity of the C=C linkers is evidenced by the disappearance of this specificity when these linkers are replaced by peptide bonds (peptilexins). On the other hand vinylexins with proton donating heteroaromatic rings are, as expected, AT specific. The vinylexin family may thus represent universal minor groove binding agents susceptible to bind to any given base pair sequence of DNA following the positioning of their proton donor and proton acceptor rings. This study confirms the insufficiency of purely geometrical and/or hydrogen bonding considerations for the correct estimation of GC versus AT specificity of groove binding ligands. These can only be accounted for by taking into consideration the overall electronic properties of the interacting species and explicitly calculating the energies of complex formation including all the relevant contributions.     

Probing the interaction of distamycin A with S100β: the “unexpected” ability of S100β to bind to DNA‐binding ligands

DNA‐minor‐groove‐binding ligands are potent antineoplastic molecules. The antibiotic distamycin A is the prototype of one class of these DNA‐interfering molecules that have been largely used in vitro. The affinity of distamycin A for DNA is well known, and the structural details of the complexes with some B‐DNA and G‐quadruplex‐forming DNA sequences have been already elucidated. Here, we show that distamycin A binds S100β, a protein involved in the regulation of several cellular processes. The reported affinity of distamycin A for the calcium(II)‐loaded S100β reinforces the idea that some biological activities of the DNA‐minor‐groove‐binding ligands arise from the binding to cellular proteins. Copyright © 2015 John Wiley & Sons, Ltd.     

Out-of-shape DNA minor groove binders: induced fit interactions of heterocyclic dications with the DNA minor groove

DB921 and DB911 are benzimidazole-biphenyl isomers with terminal charged amidines. DB911 has a central meta-substituted phenyl that gives it a shape similar to those of known minor groove binding compounds. DB921 has a central para-substituted phenyl with a linear conformation that lacks the appropriate radius of curvature to match the groove shape. It is thus expected that DB911, but not DB921, should be an effective minor groove binder, but we find that DB921 not only binds in the groove but also has an unusually high binding constant in SPR experiments (2.9 x 10(8) M(-)(1), vs 2.1 x 10(7) M(-)(1) for DB911). ITC thermodynamic analysis with an AATT sequence shows that the stronger binding of DB921 is due to a more favorable binding enthalpy relative to that of DB911. CD results support minor groove binding for both compounds but do not provide an explanation for the binding of DB921. X-ray crystallographic analysis of DB921 bound to AATT shows that an induced fit structural change in DB921 reduces the twist of the biphenyl to complement the groove, and places the functional groups in position to interact with bases at the floor of the groove. The phenylamidine of DB921 forms indirect contacts with the bases through a bound water. The DB921-water pair forms a curved binding module that matches the shape of the minor groove and provides a number of strong interactions that are not possible with DB911. This result suggests that traditional views of compound curvature required for minor groove complex formation should be reevaluated.     

Flexible docking of DNA fragments and actinocin derivatives

We have applied molecular docking methods to systems containing nucleic acids as targets and biologically active substances as ligands. The complexes of DNA fragments and actinocin derivatives with different lengths of aminoalkyl side chains were obtained by molecular docking. It was observed that actinocin derivatives could form energetically favourable complexes with DNA both as intercalators and minor groove binders. It was shown that small changes in the binding energy (～1?kcal/mol) could result in complexes with substantially different structure. The complexes of actinocin derivatives and DNA fragments were stabilized by hydrogen bonding upon intercalation and minor groove binding. It was found that the change of solvent-accessible surface area upon binding of the actinocin derivative to DNA linear increased with the growth of methylene groups' number in ligand side chains. The solvation energy change upon binding of actinocin derivatives to DNA calculated by the WSAS method was favourable in the case of small uncharged ligands and unfavourable for positively charged ligands.     

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.     

Recognition of RNA duplex by a neomycin–Hoechst 33258 conjugate

DNA minor groove binding drugs such as Hoechst 33258 have been shown to bind to a number of RNA structures. Similarly, RNA binding ligands such as neomycin have been shown by us to bind to a number of A-form DNA structures. A neomycin–Hoechst 33258 conjugate was recently shown to bind B-DNA, where Hoechst exhibits high affinity for the minor groove of A/T tract DNA and neomycin docks into the major groove. Further studies now indicate that the Hoechst moiety of the conjugate can be driven to bind RNA duplex as a consequence of neomycin binding in the RNA major groove. This is the first example of Hoechst 33258 binding to RNA duplex not containing bulges or loop motifs.     

Manipulative interplay of two adozelesin molecules with d(ATTAAT)₂achieving ligand-stacked Watson-Crick and Hoogsteen base-paired duplex adducts

Previous structural studies of the cyclopropapyrroloindole (CPI) antitumor antibiotics have shown that these ligands bind covalently edge-on into the minor groove of double-stranded DNA. Reversible covalent modification of the DNA via N3 of adenine occurs in a sequence-specific fashion. Early nuclear magnetic resonance and molecular modeling studies with both mono- and bis-alkylating ligands indicated that the ligands fit tightly within the minor groove, causing little distortion of the helix. In this study, we propose a new binding model for several of the CPI-based analogues, in which the aromatic secondary rings form π-stacked complexes within the minor groove. One of the adducts, formed with adozelesin and the d(ATTAAT)(2) sequence, also demonstrates the ability of these ligands to manipulate the DNA of the binding site, resulting in a Hoogsteen base-paired adduct. Although this type of base pairing has been previously observed with the bisfunctional CPI analogue bizelesin, this is the first time that such an observation has been made with a monoalkylating nondimeric analogue. Together, these results provide a new model for the design of CPI-based antitumor antibiotics, which also has a significant bearing on other structurally related and structurally unrelated minor groove-binding ligands. They indicate the dynamic nature of ligand-DNA interactions, demonstrating both DNA conformational flexibility and the ability of two DNA-bound ligands to interact to form stable covalent modified complexes.     

Control of DNA minor groove width and Fis protein binding by the purine 2-amino group

The width of the DNA minor groove varies with sequence and can be a major determinant of DNA shape recognition by proteins. For example, the minor groove within the center of the Fis–DNA complex narrows to about half the mean minor groove width of canonical B-form DNA to fit onto the protein surface. G/C base pairs within this segment, which is not contacted by the Fis protein, reduce binding affinities up to 2000-fold over A/T-rich sequences. We show here through multiple X-ray structures and binding properties of Fis–DNA complexes containing base analogs that the 2-amino group on guanine is the primary molecular determinant controlling minor groove widths. Molecular dynamics simulations of free-DNA targets with canonical and modified bases further demonstrate that sequence-dependent narrowing of minor groove widths is modulated almost entirely by the presence of purine 2-amino groups. We also provide evidence that protein-mediated phosphate neutralization facilitates minor groove compression and is particularly important for binding to non-optimally shaped DNA duplexes.