Molecular dynamics simulations of A. T-rich oligomers: sequence-specific binding of Na+ in the minor groove of B-DNA

Molecular dynamics simulations have been employed to probe the sequence-specific binding of sodium ions to the minor groove of B-DNA of three A. T-rich oligomers having identical compositions but different orders of the base pairs: C(AT)(4)G, CA(4)T(4)G, and CT(4)A(4)G. Recent experimental investigations, either in crystals or in solution, have shown that monovalent cations bind to DNA in a sequence-specific mode, preferentially in the narrow minor groove regions of uninterrupted sequences of four or more adenines (A-tracts), replacing a water molecule of the ordered hydration structure, the hydration spine. Following this evidence, it has been hypothesized that in A-tracts these events may be responsible for structural peculiarities such as a narrow minor groove and a curvature of the helix axis. The present simulations confirm a sequence specificity of the binding of sodium ions: Na(+) intrusions in the first layer of hydration of the minor groove, with long residence times, up to approximately 3 ns, are observed only in the minor groove of A-tracts but not in the alternating sequence. The effects of these intrusions on the structure of DNA depend on the ion coordination: when the ion replaces a water molecule of the spine, the minor groove becomes narrower. Ion intrusions may also disrupt the hydration spine modifying the oligomer structure to a large extent. However, in no case intrusions were observed to locally bend the axis toward the minor groove. The simulations also show that ions may reside for long time periods in the second layer of hydration, particularly in the wider regions of the groove, often leading to an opening of the groove.     

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.     

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.     

Role of Minor Groove Width and Hydration Pattern on Amsacrine Interaction with DNA

Amsacrine is an anilinoacridine derivative anticancer drug, used to treat a wide variety of malignancies. In cells, amsacrine poisons topoisomerase 2 by stabilizing DNA-drug-enzyme ternary complex. Presence of amsacrine increases the steady-state concentration of these ternary complexes which in turn hampers DNA replication and results in subsequent cell death. Due to reversible binding and rapid slip-out of amsacrine from DNA duplex, structural data is not available on amsacrine-DNA complexes. In the present work, we designed five oligonucleotide duplexes, differing in their minor groove widths and hydration pattern, and examined their binding with amsacrine using attenuated total reflection Fourier transform infrared (ATR-FTIR) spectroscopy. Complexes of amsacrine with calf thymus DNA were also evaluated for a comparison. Our results demonstrate for the first time that amsacrine is not a simple intercalator; rather mixed type of DNA binding (intercalation and minor groove) takes place between amsacrine and DNA. Further, this binding is highly sensitive towards the geometries and hydration patterns of different minor grooves present in the DNA. This study shows that ligand binding to DNA could be very sensitive to DNA base composition and DNA groove structures. Results demonstrated here could have implication for understanding cytotoxic mechanism of aminoacridine based anticancer drugs and provide directions to modify these drugs for better efficacy and few side-effects.     

Interaction of procarbazine with calf thymus DNA—a biophysical and molecular docking study

Interaction of procarbazine (PCZ) with calf thymus DNA was studied using biophysical and molecular docking studies. Procarbazine was to interact with DNA with a binding constant of 6.52 × 10³ M⁻¹ as calculated using ultraviolet‐visible spectroscopy. To find out the binding mode, molecular docking was performed that predicted PCZ to interact with DNA through groove binding mode with binding affinity of −6.7 kcal/mole. To confirm the groove binding nature, different experiments were performed. Dye displacement assays confirmed the non‐intercalative binding mode. Procarbazine displaced Hoechst dye from the minor groove of DNA while it was unable to displace intercalating dyes. There was no increase in the viscosity of DNA solution in presence of PCZ. Also, negligible change in the secondary structure of DNA was observed in presence of PCZ as evident by circular dichroism spectra. Procarbazine caused decrease in the melting temperature of DNA possibly because of decrease in the stability of DNA caused by groove binding interaction of PCZ with DNA.     

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.     

A study of the interaction of Cr(III) complexes and their selective binding with B-DNA: a molecular modeling approach

Molecular modeling and energy minimisation calculations have been used to investigate the interaction of chromium(III) complexes in different ligand environments with various sequences of B-DNA. The complexes are [Cr(salen)(H(2)O)(2)](+); salen denotes 1, 2 bis-salicylideneaminoethane, [Cr(salprn)(H(2)O)(2)](+); salprn denotes 1, 3 bis- salicylideneaminopropane, [Cr(phen)(3)](3+); phen denotes 1, 10 phenanthroline and [Cr(en)(3)](3+); en denotes ethylenediamine. All the chromium(III) complexes are interacted with the minor groove and major groove of d(AT)(12), d(CGCGAATTCGCG)(2) and d(GC)(12) sequences of DNA. The binding energy and hydrogen bond parameters of DNA-Cr complex adduct in both the groove have been determined using molecular mechanics approach. The binding energy and formation of hydrogen bonds between chromium(III) complex and DNA has shown that all complexes of chromium(III) prefer minor groove interaction as the favourable binding mode.     

Methylene blue binding to DNA with alternating GC base sequence: continuum treatment of salt effects

Methylene blue (MB), an efficient singlet oxygen generating photoactive dye, binds to DNA and allows photosensitized reactions to be used for sequence-specific cleavage of the DNA backbone. Intercalation and groove binding are possible binding modes of the dye, depending on base sequences and environmental conditions. In a recent modeling study of methylene blue binding to a double stranded DNA decamer with an alternating GC sequence, six structural models for intercalation structures and for minor and major groove binding have been obtained. By estimating the binding energies (including electrostatic reaction field contributions of a salt-free aqueous solvent), symmetric intercalation at the 5'-CpG-3' and 5'-GpC-3' steps was found as the predominant binding mode, followed by a slightly weaker binding of the dye in the minor groove. In this study, the stability of the modeled structures has been analysed as a function of salt concentration. The results of finite difference numerical solutions of the non-linear Poisson-Boltzmann equation show that the stabilizing effect of salt is larger for free DNA than for the modeled MB-DNA complexes. Accordingly, the estimated binding energies decrease with increasing ionic strength. A slightly higher stabilization of the groove binding complexes results in comparable binding energies for symmetric intercalation and minor groove binding at high salt concentration. Both results are in qualitative agreement with experimental data.     

Binding properties and DNA sequence-specific recognition of two bithiazole-linked netropsin hybrid molecules.

We report the DNA binding properties of two hybrid molecules which result from the combination of the DNA sequence-specific minor groove ligand netropsin with the bithiazole moiety of the antitumor drug bleomycin. The drug-DNA interaction has been investigated by means of electric linear dichroism (ELD) spectroscopy and DNase I footprinting. In compound 1 the two moieties are linked by a flexible aliphatic tether while in compound 2 the two aromatic ring systems are directly coupled by a rigid peptide bond. The results are consistent with a model in which the netropsin moiety of compound 1 resides in the minor groove of DNA and where the appended bithiazole moiety is projected away from the DNA groove. This monocationic hybrid compound has a weak affinity for DNA and shows a strict preference for A and T stretches. ELD measurements indicate that in the presence of DNA compound 2 has an orientation typical of a minor groove binder. Similar orientation angles were measured for netropsin and compound 2. This ligand which has a biscationic nature tightly binds to DNA (Ka = 6.3 x 10(5) M-1) and is mainly an AT-specific groove binder. But, depending on the nature of the sequence flanking the AT site first targeted by its netropsin moiety, the bithiazole moiety of 2 can accommodate various types of nucleotide motifs with the exception of homooligomeric sequences. As evidenced by footprinting data, the bithiazole group of bleomycin acts as a DNA recognition element, offering opportunities to recognize GC bp-containing DNA sequences with apparently a preference (although not absolute) for a pyrimidine-G-pyrimidine motif. Thus, the bithiazole unit of bleomycin provides an additional anchor for DNA binding and is also capable of specifically recognizing particular DNA sequences when it is appended to a strongly sequence selective groove binding entity. Finally, a model which schematizes the binding of compound 2 to the sequence 5'-TATGC is proposed. This model readily explains the experimentally observed specificity of this netropsin-bithiazole conjugate.     

Synthesis and DNA‐binding study of imidazole linked thiazolidinone derivatives

A novel series of imidazole‐linked thiazolidinone hybrid molecules were designed and synthesized through a feasible synthetic protocol. The molecules were characterized with Fourier transform infrared (FT‐IR), ¹H nuclear magnetic resonance (NMR), ¹³C NMR and high‐resolution mass spectrometry (HRMS) techniques. In vitro susceptibility tests against Gram‐positive (S. aureus and B. subtilis ) and Gram‐negative bacteria (E. coli and P. aeruginosa ) gave highly promising results. The most active molecule (3e) gave a minimal inhibitory concentration (MIC) value of 3.125 μg/mL which is on par with the reference drug streptomycin. Structure–activity relationships revealed activity enhancement by nitro and chloro groups when they occupied meta position of the arylidene ring in 2‐((3‐(imidazol‐1‐yl)propyl)amino)‐5‐benzylidenethiazolidin‐4‐ones. DNA‐binding study of the most potent molecule 3e with salmon milt DNA (sm‐DNA) under simulated physiological pH was probed with UV–visible absorption, fluorescence quenching, gel electrophoresis and molecular docking techniques. These studies established that compound 3e has a strong affinity towards DNA and binds at DNA minor groove with a binding constant (K_b) 0.18 × 10² L mol^?1. Molecular docking simulations predicted strong affinity of 3e towards DNA with a binding affinity (ΔG) ‐8.5 kcal/mol. Van der Waals forces, hydrogen bonding and hydrophobic interactions were predicted as the main forces of interaction. The molecule 3e exhibited specific affinity towards adenine–thiamine base pairs. Copyright © 2016 John Wiley & Sons, Ltd.     

Study on the inhibitory mechanism and binding mode of the hydroxycoumarin compound NSC158393 to HIV‐1 integrase by molecular modeling

Human immunodeficiency virus type 1 integrase (IN) is an essential enzyme in the life cycle of this virus and also an important target for the study of anti‐HIV drugs. In this work, the binding modes of the wild type IN core domain and the two mutants, that is, W132G and C130S, with the 4‐hydroxycoumarin compound NSC158393 were evaluated by using the “relaxed complex” molecular docking approach combined with molecular dynamics (MD) simulations. Based on the monomer MD simulations, both of the two substitutions affect not only the stability of the 128–136 peptides, but also the flexibility of the functional 140s loop. In principle, NSC158393 binds the 128–136 peptides of IN; however, the specific binding modes for the three systems are various. According to the binding mode of NSC158393 with WT, NSC158393 can effectively interfere with the stability of the IN dimer by causing a steric hindrance around the monomer interface. Additionally, through the comparative analysis of the MD trajectories of the wild type IN and the IN‐NSC158393 complex, we found that NSC15893 may also exert its inhibitory function by diminishing the mobility of the function loop of IN. Three key binding residues, that is, W131, K136, and G134, were discovered by energy decomposition calculated with the Molecular Mechanics Generalized Born Surface Area method. Characterized by the largest binding affinity, W131 is likely to be indispensable for the ligand binding. All the above results are consistent with experiment data, providing us some helpful information for understanding the mechanism of the coumarin‐based inhibitors. © 2009 Wiley Periodicals, Inc. Biopolymers 91: 700–709, 2009. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com