Kinetics of binding of Hoechst dyes to DNA studied by stopped-flow fluorescence techniques

Dissociation kinetics of small-molecule inhibitors in Escherichia coli is coupled to physiological state of cells

Bioactive small-molecule inhibitors represent a treasure chest for future drugs. In vitro high-throughput screening is a common approach to identify the small-molecule inhibitors that bind tightly to purified targets. Here, we investigate the inhibitor-target binding/unbinding kinetics in E. coli cells using a benzimidazole-derivative DNA inhibitor as a model system. We find that its unbinding rate is not constant but depends on cell growth rate. This dependence is mediated by the cellular activity, forming a feedback loop with the inhibitor's activity. In accordance with this feedback, we find cell-to-cell heterogeneity in inhibitor-target interaction, leading to co-existence of two distinct subpopulations: actively growing cells that dissociate the inhibitors from the targets and non-growing cells that do not. We find similar heterogeneity for other clinical DNA inhibitors. Our studies reveal a mechanism that couples inhibitor-target kinetics to cell physiology and demonstrate the significant effect of this coupling on drug efficacy.

DNA minor-groove binder Hoechst 33258 destabilizes base-pairing adjacent to its binding site

Understanding the dynamic interactions of ligands to DNA is important in DNA-based nanotechnologies. By structurally tracking the dissociation of Hoechst 33258-bound DNA (d(CGCAAATTTGCG)₂) complex (H-DNA) with T-jump 2D-IR spectroscopy, the ligand is found to strongly disturb the stability of the three C:G base pairs adjacent to A:T the binding site, with the broken base pairs being more than triple at 100 ns. The strong stabilization effect of the ligand on DNA duplex makes this observation quite striking, which dramatically increases the melting temperature and dissociation time. MD simulations demonstrate an important role of hydration water and counter cations in maintaining the separation of terminal base pairs. The hydrogen bonds between the ligand and thymine carbonyls are crucial in stabilizing H-DNA, whose breaking signal appearing prior to the complete dissociation. Thermodynamic analysis informs us that H-DNA association is a concerted process, where H cooperates with DNA single strands in forming H-DNA.

Molecular recognition of DNA by ligands: roughness and complexity of the free energy profile

Understanding the molecular mechanism by which probes and chemotherapeutic agents bind to nucleic acids is a fundamental issue in modern drug design. From a computational perspective, valuable insights are gained by the estimation of free energy landscapes as a function of some collective variables (CVs), which are associated with the molecular recognition event. Unfortunately the choice of CVs is highly non-trivial because of DNA's high flexibility and the presence of multiple association-dissociation events at different locations and/or sliding within the grooves. Here we have applied a modified version of Locally-Scaled Diffusion Map (LSDMap), a nonlinear dimensionality reduction technique for decoupling multiple-timescale dynamics in macromolecular systems, to a metadynamics-based free energy landscape calculated using a set of intuitive CVs. We investigated the binding of the organic drug anthramycin to a DNA 14-mer duplex. By performing an extensive set of metadynamics simulations, we observed sliding of anthramycin along the full-length DNA minor groove, as well as several detachments from multiple sites, including the one identified by X-ray crystallography. As in the case of equilibrium processes, the LSDMap analysis is able to extract the most relevant collective motions, which are associated with the slow processes within the system, i.e., ligand diffusion along the minor groove and dissociation from it. Thus, LSDMap in combination with metadynamics (and possibly every equivalent method) emerges as a powerful method to describe the energetics of ligand binding to DNA without resorting to intuitive ad hoc reaction coordinates.

Assessment of CA1 injury after global ischemia using supervised 2D analyses of nuclear pyknosis

Selective neuronal vulnerability is a common theme in both acute and chronic diseases affecting the nervous system. This phenomenon is particularly conspicuous after global cerebral ischemia wherein CA1 pyramidal neurons undergo delayed death while surrounding hippocampal regions are relatively spared. While injury in this model can be easily demonstrated using either histological or immunological stains, current methods used to assess the cellular injury present in these biological images lack the precision required to adequately compare treatment effects. To address this shortcoming, we devised a supervised work-flow that can be used to quantify ischemia-induced nuclear condensation using microscopic images. And while we demonstrate the utility of this technique using models of ischemic brain injury, the approach can be readily applied to other paradigms in which programmed cell death is a major component.

DNA interstrand cross-links of an antitumor trinuclear platinum(II) complex: thermodynamic analysis and chemical probing

The trinuclear platinum compound [{trans-PtCl(NH(3))(2)}(2)(μ-trans-Pt(NH(3))(2){NH(2)(CH(2))(6)NH(2)}(2))](4+) (BBR3464) belongs to the polynuclear class of platinum-based anticancer agents. These agents form in DNA long-range (Pt,Pt) interstrand cross-links, whose role in the antitumor effects of BBR3464 predominates. Our results show for the first time that the interstrand cross-links formed by BBR3464 between two guanine bases in opposite strands separated by two base pairs (1,4-interstrand cross-links) exist as two distinct conformers, which are not interconvertible, not only if these cross-links are formed in the 5'-5', but also in the less-usual 3'-3' direction. Analysis of the conformers by differential scanning calorimetry, chemical probes of DNA conformation, and minor groove binder Hoechst 33258 demonstrate that each of the four conformers affects DNA in a distinctly different way and adopts a different conformation. The results also support the thesis that the molecule of antitumor BBR3464 when forming DNA interstrand cross-links may adopt different global structures, including different configurations of the linker chain of BBR3464 in the minor groove of DNA. Our findings suggest that the multiple DNA interstrand cross-links available to BBR3464 may all contribute substantially to its cytotoxicity.

Molecular motions in drug design: the coming age of the metadynamics method

Metadynamics is emerging as a useful free energy method in physics, chemistry and biology. Recently, it has been applied also to investigate ligand binding to biomolecules of pharmacological interest. Here, after introducing the basic idea of the method, we review applications to challenging targets for pharmaceutical intervention. We show that this methodology, especially when combined with a variety of other computational approaches such as molecular docking and/or molecular dynamics simulation, may be useful to predict structure and energetics of ligand/target complexes even when the targets lack a deep binding cavity, such as DNA and proteins undergoing fibrillation in neurodegenerative diseases. Furthermore, the method allows investigating the routes of molecular recognition and the associated binding energy profiles, providing a molecular interpretation to experimental data.

Triple recognition of B-DNA by a neomycin-Hoechst 33258-pyrene conjugate

Recent developments have indicated that aminoglycoside binding is not limited to RNA, but to nucleic acids that, like RNA, adopt conformations similar to its A-form. We further sought to expand the utility of aminoglycoside binding to B-DNA structures by conjugating neomycin, an aminoglycoside antibiotic, with the B-DNA minor groove binding ligand Hoechst 33258. Envisioning a dual groove binding mode, we have extended the potential recognition process to include a third, intercalative moiety. Similar conjugates, which vary in the number of binding moieties but maintain identical linkages to allow direct comparisons to be made, have also been prepared. We report herein novel neomycin- and Hoechst 33258-based conjugates developed in our laboratories for exploring the recognition potential with B-DNA. Spectroscopic studies such as UV melting, differential scanning calorimetry, isothermal fluorescence titrations, and circular dichroism together illustrate the triple recognition of the novel conjugate containing neomycin, Hoechst 33258, and pyrene. This study represents the first example of DNA molecular recognition capable of minor versus major groove recognition in conjunction with intercalation.

Dissociation of minor groove binders from DNA: insights from metadynamics simulations

We have used metadynamics to investigate the mechanism of noncovalent dissociation from DNA by two representatives of alkylating and noncovalent minor groove (MG) binders. The compounds are anthramycin in its anhydrous form (IMI) and distamycin A (DST), which differ in mode of binding, size, flexibility and net charge. This choice enables to evaluate the influence of such factors on the mechanism of dissociation. Dissociation of IMI requires an activation free energy of approximately 12 kcal/mol and occurs via local widening of the MG and loss of contacts between the drug and one DNA strand, along with the insertion of waters in between. The detachment of DST occurs at a larger free energy cost, approximately 16.5 or approximately 18 kcal/mol depending on the binding mode. These values compare well with that of 16.6 kcal/mol extracted from stopped-flow experiments. In contrast to IMI, an intermediate is found in which the ligand is anchored to the DNA through its amidinium tail. From this conformation, binding and unbinding occur almost at the same rate. Comparison between DST and with kinetic models for the dissociation of Hoechst 33258 from DNA uncovers common characteristics across different classes of noncovalent MG ligands.

Hydration changes in the association of Hoechst 33258 with DNA

The role of water in the interaction of Hoechst 33258 with the minor groove binding site of the (AATT)2 sequence was investigated using calorimetric and equilibrium constant measurements. Using isothermal titration calorimetry measurements, the heat capacity change for the reaction is -256 +/- 10 cal/(K mol of Hoechst). Comparison with the heat capacity changes based on area models supports the expulsion of water from the interface of the Hoechst-DNA complex. To further consider the role of water, the osmotic stress method was used to determine if the Hoechst association with DNA was coupled with hydration changes. Using four osmolytes with varying molecular weights and chemical properties, the Hoechst affinity for DNA decreases with increasing osmolyte concentration. From the dependence of the equilibrium constant on the solution osmolality, 60 +/- 13 waters are acquired in the complex relative to the reactants. It is proposed that the osmotic stress technique is measuring weakly bound waters that are not measured via the heat capacity changes.

Base-sequence specificity of Hoechst 33258 and DAPI binding to five (A/T)4 DNA sites with kinetic evidence for more than one high-affinity Hoechst 33258-AATT complex

The binding of Hoechst 33258 and DAPI to five different (A/T)4 sequences in a stable DNA hairpin was studied exploiting the substantial increase in dye fluorescence upon binding. The two dyes have comparable affinities for the AATT site (e.g. association constant K(a)=5.5 x 10(8) M(-1) for DAPI), and their affinities decrease in the series AATT >> TAAT approximately equal to ATAT > TATA approximately equal to TTAA. The extreme values of K(a) differ by a factor of 200 for Hoechst 33258 but only 30 for DAPI. The binding kinetics of Hoechst 33258 were measured by stopped-flow under pseudo-first order conditions with an (A/T)4 site in excess. The lower-resolution experiments can be well represented by single exponential processes, corresponding to a single-step binding mechanism. The calculated association-rate parameters for the five (A/T)4 sites are similar (2.46 x 10(8) M(-1) s(-1) to 0.86 x 10(8) M(-1) s(-1)) and nearly diffusion-controlled, while the dissociation-rate parameters vary from 0.42 s(-1) to 96 s(-1). Thus the association constants are kinetically controlled and are close to their equilibrium-determined values. However, when obtained with increased signal-to-noise ratio, the kinetic traces for Hoechst 33258 binding at the AATT site reveal two components. The concentration dependencies of the two time constants and amplitudes are consistent with two different kinetically equivalent two-step models. In the first model, fast bimolecular binding is followed by an isomerization of the initial complex. In the second model, two single-step associations form two complexes that mutually exclude each other. For both models the four reaction-rate parameters are calculated. Finally, specific dissociation kinetics, using poly[d(A-5BrU)], show that the kinetics are even more complex than either two-step model. We correlate our results with the different binding orientations and locations of Hoechst 33258 in the DNA minor groove found in several structural studies in the literature.