The width of the minor groove affects the binding of the bisquaternary heterocycle SN-6999 to duplex DNA.

A complex between d(GGGAAAAACGG).d(CCGTTTTTCCC) and the minor groove binding drug SN-6999 has been studied by 1H nuclear magnetic resonance spectroscopy. The drug is found to bind in the d(A)5 tract, but with interactions extending one residue in the 3'-direction along each strand. Doubling of resonances in the complex indicates slow to intermediate exchange between two binding modes. An orientational preference (7:3) is found, the first such example in an SN-6999 complex. Furthermore, the upper limit of the lifetime for the major species is longer than was found for SN-6999 with other DNA duplexes. The preferred orientation of SN-6999 has the pyridinium ring near the 5'-end of the (+) strand; the minor binding mode has the reverse orientation. The orientational preference and slower exchange rate relative to other SN-6999 complexes is attributed to increased stabilization from van der Waals interactions due to better shape complementarity between the DNA duplex and ligand. The comparison of these results with studies of SN-6999 complexed to other DNA duplexes reveals the sensitivity of the binding properties to the delicate interplay between the molecular structure of the ligand and the specific characteristics of the DNA minor groove.     

Anthracycline antibiotic arugomycin binds in both grooves of the DNA helix simultaneously: an NMR and molecular modelling study.

Perturbations to the 1H and 31P chemical shifts of DNA resonances together with twenty-four intermolecular nuclear Overhauser effects show that the anthracycline antibiotic arugomycin intercalates between the basepairs of the hexamer duplex d(5'-GCATGC)2 at the 5'-CpA and 5'-TpG binding sites. In the complex two drug molecules are bound per duplex with full retention of the dyad symmetry. Arugomycin adopts a threaded binding orientation with chains of sugars positioned in both the major and minor groove of the helix simultaneously. The complex is stabilized by hydrogen bonding, electrostatic and van der Waals interactions principally in the major groove and involving substituents on the rigidly oriented bicycloamino-glucose sugar of the antibiotic. A specific hydrogen bond is identified between the C2'-hydroxyl and the guanine N7 at the intercalation site. Together, interactions in the major groove appear to account for the intercalation specificity of arugomycin that requires both a guanine and thymine at the intercalation site. We are unable to identify any sequence specific interactions between the minor groove and the arugarose sugar (S1) which binds only weakly, through van der Walls contacts, over the d(GCA).d(TGC) trinucleotide sequence. The data indicate that the sugar chains of arugomycin are flexible and play little part in the interaction of the antibiotic with DNA. The intensity of sequential internucleotide NOEs identifies the intercalation site as being assymmetric. A family of conformers computed using restrained energy minimisation and molecular dynamics indicate that basepair buckling is a feature of the anthracycline intercalation site that may serve to maximise intermolecular van der Waals interactions by wrapping the basepairs around the antibiotic chromophore.     

Multi-spectroscopic methods combined with molecular modeling dissect the interaction mechanisms of ractopamine and calf thymus DNA

The toxic interaction of ractopamine (RAC) with calf thymus DNA (ct DNA) was studied in vitro using multi-spectroscopic methods and molecular modeling methods. The hypochromic effect without a noticeable shift in UV-vis absorption indicated that the minor groove binding mode existed in the interaction between RAC and DNA. The fluorescence quenching of RAC was observed with the increasing addition of DNA and was proved to be the static quenching. The binding constant and the binding site sizes were 4.13 × 10(3) and 0.97, respectively. The thermodynamic calculation demonstrated that the hydrogen bond and van der Waals were main acting forces. This result further confirmed the existence of groove binding mode. Afterwards, we found another interaction mode, electrostatic binding mode through the fluorescence polarization, ionic effects and denatured DNA experiments. Circular dichroism spectroscopy (CD) was then employed to monitor the conformation changes of DNA. Molecular modeling studies illustrated the visual display of the binding mode and the detailed information of the H-bond.     

Induced fit DNA recognition by a minor groove binding analogue of Hoechst 33258: fluctuations in DNA A tract structure investigated by NMR and molecular dynamics simulations

NMR analysis and molecular dynamics simulations of d(GGTAATTACC)₂ and its complex with a tetrahydropyrimidinium analogue of Hoechst 33258 suggest that DNA minor groove recognition in solution involves a combination of conformational selection and induced fit, rather than binding to a preorganised site. Analysis of structural fluctuations in the bound and unbound states suggests that the degree of induced fit observed is primarily a consequence of optimising van der Waals contacts with the walls of the minor groove resulting in groove narrowing through: (i) changes in base step parameters, including increased helical twist and propeller twist; (ii) changes to the sugar–phosphate backbone conformation to engulf the bound ligand; (iii) suppression of bending modes at the TpA steps. In contrast, the geometrical arrangement of hydrogen bond acceptors on the groove floor appears to be relatively insensitive to DNA conformation (helical twist and propeller twist). We suggest that effective recognition of DNA sequences (in this case an A tract structure) appears to depend to a significant extent on the sequence being flexible enough to be able to adopt the geometrically optimal conformation compatible with the various binding interactions, rather than involving ‘lock and key’ recognition.     

Variability in DNA minor groove width recognised by ligand binding: the crystal structure of a bis-benzimidazole compound bound to the DNA duplex d(CGCGAATTCGCG)2.

An analogue of the DNA-binding compound Hoechst 33258, in which the piperazine ring has been replaced by an imidazoline group, has been cocrystallized with the dodecanucleotide sequence d(CGCGAATTCGCG)2. The structure has been solved by X-ray diffraction analysis and has been refined to an R-factor of 19.7% at a resolution of 2.0 A. The ligand is found to bind in the minor groove, at the central four AATT base pairs of the B-DNA double helix, with the involvement of a number of van der Waals contacts and hydrogen bonds. There are significant differences in minor groove width for the two compounds, along much of the AATT region. In particular this structure shows a narrower groove at the 3' end of the binding site consistent with the narrower cross-section of the imidazole group compared with the piperazine ring of Hoechst 33258 and therefore a smaller perturbation in groove width. The higher binding affinity to DNA shown by this analogue compared with Hoechst 33258 itself, has been rationalised in terms of these differences.     

The Mu repressor-DNA complex contains an immobilized 'wing' within the minor groove

We have determined the solution structure of the complex between the 'winged-helix' enhancer binding domain of the Mu repressor protein and its cognate DNA site. The structure reveals an unusual use for the 'wing' which becomes immobilized upon DNA binding where it makes intermolecular hydrogen bond contacts deep within the minor groove. Although the wing is mobile in the absence of DNA, it partially negates the large entropic penalty associated with its burial by maintaining a small degree of structural order in the DNA-free state. Extensive contacts are also formed between the recognition helix and the DNA, which reads the major groove of a highly conserved region of the binding site through a single base-specific hydrogen bond and van der Waals contacts.     

A demonstration of the intrinsic importance of stabilizing hydrophobic binding and non-covalent van der Waals contacts dominant in the non-covalent CC-1065/B-DNA binding

The comparative DNA binding properties and cytotoxic activity of CDPIn methyl esters (n = 1-5) vs. PDE-In methyl esters (n = 1-3) are detailed in studies which provide experimental evidence for the intrinsic importance of stabilizing hydrophobic binding and non-covalent van der Waals contacts dominant in the CC-1065/B-DNA minor groove binding. High affinity minor groove binding to DNA was established through: (1) the observation of CDPI3 binding (UV) but not unwinding of supercoiled DNA (phi 174 RFI DNA) thus excluding intercalative binding; (2) the observation of CDPI3 binding to T4 phage DNA (UV, delta Tm) in which the major groove is occluded by glycosylation thus excluding major groove binding; (3) the observation of salt (Na+) concentration independent high affinity CDPI3 binding to poly(dA . poly(dT) thus excluding simple electrostatic binding to the DNA phosphate backbone; and further inferred through (4) the observation of an intense induced dichroism (ICD, poly(dA) . poly(dT) and poly(dG) . poly(dC) [phi]23(358) = 24,000 and 23,500). This high affinity minor groove binding is sufficient to produce a potent cytotoxic effect.(ABSTRACT TRUNCATED AT 400 WORDS)     

Molecular modelling of the interaction of carbocyclic analogues of netropsin and distamycin with d(CGCGAATTCGCG)2

A molecular mechanics and molecular dynamics approach was used to examine the structure of complexes formed between the d(CGCGAATTCGCG)2 duplex and netropsin, distamycin, and four carbocyclic analogues of netropsin and distamycin (1-4). The resulting structures of the ligand-DNA model complexes and their energetics were examined. It is predicted that the compounds 1-4 should have a decreased affinity for the minor groove of AT-rich regions in comparison to netropsin and distamycin. From the energetic analysis it appears that van der Waals and electrostatic interactions are more important than specific hydrogen bonds in stabilizing the ligand-duplex complexes. We predict that compounds 1 and 2 are effectively isohelical with the DNA minor groove. The superior DNA-binding afforded by 1 and 2 in comparison to 3 and 4 results from their more effective penetration into the minor groove and smaller perturbation of molecular structure upon complex formation.     

Molecular modeling study of the norfloxacin-DNA complex

Molecular modeling and molecular dynamics were performed to investigate the interaction of norfloxacin with the DNA oligonucleotide 5'-d(ATACGTAT)(2). Eight quinolone-DNA binding structures were built by molecular modeling on the basis of experimental results. A 100ps molecular dynamics calculation was carried out on two groove binding models and six partially intercalating models. The resulting average structures were compared with each other and to free DNA structure as a reference. The favorable binding mode of norfloxacin to a DNA substrate was pursued by structural assess including steric hindrance, presence of hydrogen-bonding, non-bonding energies of the complex and presence of abnormal structural distortion. Although two of the intercalative models showed the highest binding energy and the lowest non-bonding interaction energy, they presented structural features which contrast with experimental results. On the other hand, one groove binding model demonstrated the most acceptable structure when the experimental observation was accounted. In this model, hydrogen bonding of the carbonyl and carboxyl group of the norfloxacin rings with the DNA bases was present, and norfloxacin binds to the amine group of the guanine base which protrudes toward the minor groove of B-DNA.     

Targeting Telomeres: Molecular Dynamics and Free Energy Simulation of Gold-Carbene Binding to DNA

DNA sequences in regulatory regions and in telomers at the ends of chromosomes frequently contain tandem repeats of guanine nucleotides that can form stacked structures stabilized by Hoogsten pairing and centrally bound monovalent cations. The replication and elongation of telomeres requires the disruption of these G-quadruplex structures. Hence, drug molecules such as gold (Au)-carbene that stabilize G-quadruplexes may also interfere with the elongation of telomeres and, in turn, could be used to control cell replication and growth. To better understand the molecular mechanism of Au-carbene binding to G-quadruplexes, we employed molecular dynamics simulations and free energy simulations. Whereas very restricted mobility of two Au-carbene ligands was found upon binding as a doublet to one side of the G-quadruplex, much larger translational and orientational mobility was observed for a single Au-carbene binding at the second G-quadruplex surface. Comparative simulations on duplex DNA in the presence of Au-carbene ligands indicates a preference for the minor groove and weaker unspecific and more salt-dependent binding than to the G-quadruplex surface. Analysis of energetic contributions reveals a dominance of nonpolar and van der Waals interactions to drive binding. The simulations can also be helpful for proposing possible modifications that could improve Au-carbene affinity and specificity for G-quadruplex binding.     

Designer DNA-binding drugs: the crystal structure of a meta-hydroxy analogue of Hoechst 33258 bound to d(CGCGAATTCGCG)2.

An analogue of the DNA binding compound Hoechst 33258, which has the para hydroxyl group altered to be at the meta position, together with the replacement of one benzimidazole group by pyridylimidazole, has been cocrystallized with the dodecanucleotide sequence d(CGCGAATTCGCG)2. The X-ray structure has been determined at 2.2 A resolution and refined to an R factor of 20.1%. The ligand binds in the minor groove at the sequence 5'-AATTC with the bulky piperazine group extending over the CxG base pair. This binding is stabilised by hydrogen bonding and numerous close van der Waals contacts to the surface of the groove walls. The meta-hydroxyl group was found in two distinct orientations, neither of which participates in direct hydrogen bonds to the exocyclic amino group of a guanine base. The conformation of the drug differs from that found previously in other X-ray structures of Hoechst 33258-DNA complexes. There is significant variation between the minor groove widths in the complexes of Hoechst 33258 and the meta-hydroxyl derivative as a result of these conformational differences. Reasons are discussed for the inability of this derivative to actively recognise guanine.     

Towards an understanding of DNA recognition by the methyl-CpG binding domain 1

CpG methylation determines a variety of biological functions of DNA. The methylation signal is interpreted by proteins containing a methyl-CpG binding domain (MBDs). Based on the NMR structure of MBD1 complexed with methylated DNA we analysed the recognition mode by means of molecular dynamics simulations. As the protein is monomeric and recognizes a symmetrically methylated CpG step, the recognition mode is an asymmetric one. We find that the two methyl groups do not contribute equally to the binding energy. One methyl group is associated with the major part of the binding energy and the other one nearly does not contribute at all. The contribution of the two cytosine methyl groups to binding energy is calculated to be -3.6 kcal/mol. This implies a contribution of greater than two orders of magnitude to the binding constant. The conserved amino acid Asp32 is known to be essential for DNA binding by MBD1, but so far no direct contact with DNA has been observed. We detected a direct DNA base contact to Asp32. This could be the main reason for the importance of this amino acid. MBD contacts DNA exclusively in the major groove, the minor groove is reserved for histone contacts. We found a deformation of the minor groove shape due to complexation by MBD1, which indicates an information transfer between the major and the minor groove.     

Towards an Understanding of DNA Recognition by the Methyl-CpG Binding Domain 1

Abstract

CpG methylation determines a variety of biological functions of DNA. The methylation signal is interpreted by proteins containing a methyl-CpG binding domain (MBDs). Based on the NMR structure of MBD1 complexed with methylated DNA we analysed the recognition mode by means of molecular dynamics simulations.

As the protein is monomeric and recognizes a symmetrically methylated CpG step, the recognition mode is an asymmetric one. We find that the two methyl groups do not contribute equally to the binding energy. One methyl group is associated with the major part of the binding energy and the other one nearly does not contribute at all. The contribution of the two cytosine methyl groups to binding energy is calculated to be ?3.6 kcal/mol. This implies a contribution of greater than two orders of magnitude to the binding constant. The conserved amino acid Asp32 is known to be essential for DNA binding by MBD1, but so far no direct contact with DNA has been observed. We detected a direct DNA base contact to Asp32. This could be the main reason for the importance of this amino acid. MBD contacts DNA exclusively in the major groove, the minor groove is reserved for histone contacts. We found a deformation of the minor groove shape due to complexation by MBD1, which indicates an information transfer between the major and the minor groove.     

Molecular Modeling Study of the Norfloxacin-DNA Complex

Abstract

Molecular modeling and molecular dynamics were performed to investigate the interaction of norfloxacin with the DNA oligonucleotide 5′-d(ATACGTAT)₂. Eight quinolone-DNA binding structures were built by molecular modeling on the basis of experimental results. A 100ps molecular dynamics calculation was carried out on two groove binding models and six partially intercalating models. The resulting average structures were compared with each other and to free DNA structure as a reference. The favorable binding mode of norfloxacin to a DNA substrate was pursued by structural assess including steric hindrance, presence of hydrogen-bonding, non-bonding energies of the complex and presence of abnormal structural distortion. Although two of the intercalative models showed the highest binding energy and the lowest non-bonding interaction energy, they presented structural features which contrast with experimental results. On the other hand, one groove binding model demonstrated the most acceptable structure when the experimental observation was accounted. In this model, hydrogen bonding of the carbonyl and carboxyl group of the norfloxacin rings with the DNA bases was present, and norfloxacin binds to the amine group of the guanine base which protrudes toward the minor groove of B-DNA.     

Probing the interaction mechanism of small molecule inhibitors with matriptase based on molecular dynamics simulation and free energy calculations

Matriptase is a serine protease associated with a wide variety of human tumors and carcinoma progression. Up to now, many promising anti-cancer drugs have been developed. However, the detailed structure–function relationship between inhibitors and matriptase remains elusive. In this work, molecular dynamics simulation and binding free energy studies were performed to investigate the biochemistry behaviors of two class inhibitors binding to matriptase. The binding free energies predicted by MM/GBSA methods are in good agreement with the experimental bioactivities, and the analysis of the individual energy terms suggests that the van der Waals interaction is the major driving force for ligand binding. The key residues responsible for achieving strong binding have been identified by the MM/GBSA free energy decomposition analysis. Especially, Trp215 and Phe99 had an important impact on active site architecture and ligand binding. The results clearly identify the two class inhibitors exist different binding modes. Through summarizing the two different modes, we have mastered some important and favorable interaction patterns between matriptase and inhibitors. Our findings would be helpful for understanding the interaction mechanism between the inhibitor and matriptase and afford important guidance for the rational design of potent matriptase inhibitors.     

A molecular thermodynamic view of DNA-drug interactions: a case study of 25 minor-groove binders

Developing a molecular view of the thermodynamics of DNA recognition is essential to the design of ligands for regulating gene expression. In a first comprehensive attempt at sketching an atlas of DNA-drug energetics, we present here a detailed thermodynamic view of minor-groove recognition by small molecules via a computational study on 25 DNA-drug complexes. The studies are configured in the MMGBSA (Molecular Mechanics-Generalized Born-Solvent Accessibility) framework at the current state of the art and facilitate a structure-energy component correlation. Analyses were conducted on both energy minimized structures of DNA-drug complexes and molecular dynamics trajectories developed for the purpose of this study. While highlighting the favorable role of packing, shape complementarity, and van der Waals and hydrophobic interactions of the drugs in the minor groove in conformity with experiment, the studies reveal an interesting annihilation of favorable electrostatics by desolvation. Structural modifications attempted on the ligands point to the requisite physico-chemical factors for obtaining improved binding energies. Hydrogen bonds predicted to be important for specificity based on structural considerations do not always turn out to be significant to binding in post facto analyses of molecular dynamics trajectories, which treat thermal averaging, solvent, and counterion effects rigorously. The strength of the hydrogen bonds retained between the DNA and drug during the molecular dynamics simulations is approximately 1kcal/mol. Overall, the study reveals the compensatory nature of the diverse binding free energy components, possible threshold limits for some of these properties, and the availability of a computationally viable free energy methodology which could be of value in drug-design endeavors.     

Binding properties of palmatine to DNA: spectroscopic and molecular modeling investigations

Palmatine, an isoquinoline alkaloid, is an important medicinal herbal extract with diverse pharmacological and biological properties. In this work, spectroscopic and molecular modeling approaches were employed to reveal the interaction between palmatine and DNA isolated from herring sperm. The absorption spectra and iodide quenching results indicated that groove binding was the main binding mode of palmatine to DNA. Fluorescence studies indicated that the binding constant (K) of palmatine and DNA was ~ 10⁴ L·mol^?1. The associated thermodynamic parameters, ΔG, ΔH, and ΔS, indicated that hydrogen bonds and van der Waals forces played major roles in the interaction. The effects of chemical denaturant, thermal denaturation and pH on the interaction were investigated and provided further support for the groove binding mode. In addition to experimental approaches, molecular modeling was conducted to verify binding pattern of palmatine–DNA. Copyright © 2015 John Wiley & Sons, Ltd.     

Interaction of berenil with the tyrT DNA sequence studied by footprinting and molecular modelling. Implications for the design of sequence-specific DNA recognition agents.

We have developed a technique of partially-restrained molecular mechanics enthalpy minimisation which enables the sequence-dependence of the DNA binding of a non-intercalating ligand to be studied for arbitrary sequences of considerable length (greater than = 60 base-pairs). The technique has been applied to analyse the binding of berenil to the minor groove of a 60 base-pair sequence derived from the tyrT promoter; the results are compared with those obtained by DNAse I and hydroxyl radical footprinting on the same sequence. The calculated and experimentally observed patterns of binding are in good agreement. Analysis of the modelling data highlights the importance of DNA flexibility in ligand binding. Further, the electrostatic component of the interaction tends to favour binding to AT-rich regions, whilst the van der Waals interaction energy term favours GC-rich ones. The results also suggest that an important contribution to the observed preference for binding in AT-rich regions arises from lower DNA perturbation energies and is not accompanied by reduced DNA structural perturbations in such sequences. It is therefore concluded that those modes of DNA distortion favourable to binding are probably more flexible in AT-rich regions. The structure of the modelled DNA sequence has also been analysed in terms of helical parameters. For the DNA energy-minimised in the absence of berenil, certain helical parameters show marked sequence-dependence. For example, purine-pyrimidine (R-Y) base pairs show a consistent positive buckle whereas this feature is consistently negative for Y-R pairs. Further, CG steps show lower than average values of slide while GC steps show lower than average values of rise. Similar analysis of the modelling data from the calculations including berenil highlights the importance of DNA flexibility in ligand binding. We observe that the binding of berenil induces characteristic responses in different helical parameters for the base-pairs around the binding site. For example, buckle and tilt tend to become more negative to the 5'-side of the binding site and more positive to the 3'-side, while the base steps at either side of the centre of the site show increased twist and decreased roll.     

Spectroscopic studies of 9-hydroxyellipticine binding to DNA

The binding of 9-hydroxyellipticine to calf thymus DNA, poly[d(A-T)]₂, and poly-[d(G-C)]₂ has been studied in detail by means of CD, linear dichroism, resonance light scattering, and molecular dynamics. The transition moment polarizations of 9-hydroxyelliptiycine were determined in polyvinyl alcohol stretched film. Spectroscopic solution studies of the DNA/drug complex are combined with theoretical CD calculations using the final 50 ps of a series of molecular dynamics simulations as input. The spectroscopic data shows 9-hydroxyellipticine to adopt two main binding modes, one intercalative and the other a stacked binding mode involving the formation of drug oligomers in the DNA major groove. Analysis of the intercalated binding mode in poly[d(A-T)]₂ suggests the 9-hydroxyellipticine hydroxyl group lies in the minor groove and hydrogen bonds to water with the pyridine ring protruding into the major groove. The stacked binding mode was examined using resonance light scattering and it was concluded that the drug was forming small oligomer stacks rather than extended aggregates. Reduced linear dichroism measurements suggested a binding geometry that precluded a minor groove binding mode where the plane of the drug makes a 45° angle with the plane of the bases. Thus it was concluded that the drug stacks in the major groove. No obvious differences in the mode of binding of 9-hydroxyellipticine were observed between different DNA sequences; however, the stacked binding mode appeared to be more favorable for calf thymus DNA and poly[d(G-C)]₂ than for poly[d(A-T)]₂, an observation that could be explained by the slightly greater steric hindrance of the poly[d(A-T)]₂ major groove. A strong concentration dependence was observed for the two binding modes where intercalation is favored at very low drug load, with stacking interactions becoming more prominent as the drug concentration is increased. Even at DNA : drug mixing ratios of 70:1 the stacked binding mode was still important for GC-rich DNAs. © 1998 John Wiley & Sons, Inc. Biopoly 46: 127–143, 1998