Atomistic account of structural and dynamical changes induced by small binders in the double helix of a short DNA

Predicting the effect of binding molecules on the shape and mechanical properties of structured DNA assemblies

Chemo-mechanical deformation of structured DNA assemblies driven by DNA-binding ligands has offered promising avenues for biological and therapeutic applications. However, it remains elusive how to effectively model and predict their effects on the deformation and mechanical properties of DNA structures. Here, we present a computational framework for simulating chemo-mechanical change of structured DNA assemblies. We particularly quantify the effects of ethidium bromide (EtBr) intercalation on the geometry and mechanical properties of DNA base-pairs through molecular dynamics simulations and integrated them into finite-element-based structural analysis to predict the shape and properties of DNA objects. The proposed model captures various structural changes induced by EtBr-binding such as shape variation, flexibility modulation, and supercoiling instability. It enables a rational design of structured DNA assemblies with tunable shapes and mechanical properties by binding molecules.

Combining steady state and temperature jump IR spectroscopy to investigate the allosteric effects of ligand binding to dsDNA

Changes in the structural dynamics of double stranded (ds)DNA upon ligand binding have been linked to the mechanism of allostery without conformational change, but direct experimental evidence remains elusive. To address this, a combination of steady state infrared (IR) absorption spectroscopy and ultrafast temperature jump IR absorption measurements has been used to quantify the extent of fast (∼100 ns) fluctuations in (ds)DNA·Hoechst 33258 complexes at a range of temperatures. Exploiting the direct link between vibrational band intensities and base stacking shows that the absolute magnitude of the change in absorbance caused by fast structural fluctuations following the temperature jump is only weakly dependent on the starting temperature of the sample. The observed fast dynamics are some two orders of magnitude faster than strand separation and associated with all points along the 10-base pair duplex d(GCATATATCC). Binding the Hoechst 33258 ligand causes a small but consistent reduction in the extent of these fast fluctuations of base pairs located outside of the ligand binding region. These observations point to a ligand-induced reduction in the flexibility of the dsDNA near the binding site, consistent with an estimated allosteric propagation length of 15 Å, about 5 base pairs, which agrees well with both molecular simulation and coarse-grained statistical mechanics models of allostery leading to cooperative ligand binding.

Ethidium bromide interactions with DNA: an exploration of a classic DNA-ligand complex with unbiased molecular dynamics simulations

Visualization of double stranded DNA in gels with the binding of the fluorescent dye ethidium bromide has been a basic experimental technique in any molecular biology laboratory for >40 years. The interaction between ethidium and double stranded DNA has been observed to be an intercalation between base pairs with strong experimental evidence. This presents a unique opportunity for computational chemistry and biomolecular simulation techniques to benchmark and assess their models in order to see if the theory can reproduce experiments and ultimately provide new insights. We present molecular dynamics simulations of the interaction of ethidium with two different double stranded DNA models. The first model system is the classic sequence d(CGCGAATTCGCG)₂ also known as the Drew–Dickerson dodecamer. We found that the ethidium ligand binds mainly stacked on, or intercalated between, the terminal base pairs of the DNA with little to no interaction with the inner base pairs. As the intercalation at the terminal CpG steps is relatively rapid, the resultant DNA unwinding, rigidification, and increased stability of the internal base pair steps inhibits further intercalation. In order to reduce these interactions and to provide a larger groove space, a second 18-mer DNA duplex system with the sequence d(GCATGAACGAACGAACGC) was tested. We computed molecular dynamics simulations for 20 independent replicas with this sequence, each with ∼27 μs of sampling time. Results show several spontaneous intercalation and base-pair eversion events that are consistent with experimental observations. The present work suggests that extended MD simulations with modern DNA force fields and optimized simulation codes are allowing the ability to reproduce unbiased intercalation events that we were not able to previously reach due to limits in computing power and the lack of extensively tested force fields and analysis tools.

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.

Chirality of a rhodamine heterodimer linked to a DNA scaffold: an experimental and computational study

The chiroptical properties of multi-chromophoric systems are governed by the intermolecular arrangement of the monomeric units. We report on a computational and experimental study of the linear optical properties and supramolecular structure of a rhodamine heterodimer assembled on a DNA scaffold. The experimental absorption and circular dichroism (CD) profiles confirm the dimer formation. Computationally, starting from low-cost DFT/TDDFT simulations of the bare dimer we attribute the measured -/+ CD sign sequence of the S1/S2 bands to a specific chiral conformation of the heterodimer. In the monomers, as typical for rhodamine dyes, the electric transition dipole of the lowest π-π* transition is parallel to the long axis of the xanthene planes. We show that in the heterodimer the sign sequence of the two CD bands is related to the orientation of these long axes. To account explicitly for environment effects, we use molecular dynamics (MD) simulations for characterizing the supramolecular structure of the two optical isomers tethered on DNA. Average absorption and CD-profiles were modeled using ab initio TDDFT calculations at the geometries sampled along a few nanosecond MD run. The absorption profiles computed for both optical isomers are in good agreement with the experimental absorption spectrum and do not allow one to discriminate between them. The computed averaged CD profiles provide the orientation of monomers in the enantiomer that is dominant under the experimental conditions.

Understanding enhanced mechanical stability of DNA in the presence of intercalated anticancer drug: Implications for DNA associated processes

Most of the anticancer drugs bind to double-stranded DNA (dsDNA) by intercalative-binding mode. Although experimental studies have become available recently, a molecular-level understanding of the interactions between the drug and dsDNA that lead to the stability of the intercalated drug is lacking. Of particular interest are the modifications of the mechanical properties of dsDNA observed in experiments. The latter could affect many biological functions, such as DNA transcription and replication. Here, we probe, via all-atom molecular dynamics (MD) simulations, the change in the mechanical properties of intercalated drug-DNA complexes for two intercalators, daunomycin and ethidium. We find that, upon drug intercalation, the stretch modulus of DNA increases significantly, whereas its persistence length and bending modulus decrease. Steered MD simulations reveal that it requires higher forces to stretch the intercalated dsDNA complexes than the normal dsDNA. Adopting various pulling protocols to study force-induced DNA melting, we find that the dissociation of dsDNA becomes difficult in the presence of intercalators. The results obtained here provide a plausible mechanism of function of the anticancer drugs, i.e., via altering the mechanical properties of DNA. We also discuss long-time consequences of using these drugs, which require further in vivo investigations.

Coherent electronic and nuclear dynamics in a rhodamine heterodimer-DNA supramolecular complex

Elucidating the role of quantum coherences in energy migration within biological and artificial multichromophoric antenna systems is the subject of an intense debate. It is also a practical matter because of the decisive implications for understanding the biological processes and engineering artificial materials for solar energy harvesting. A supramolecular rhodamine heterodimer on a DNA scaffold was suitably engineered to mimic the basic donor-acceptor unit of light-harvesting antennas. Ultrafast 2D electronic spectroscopic measurements allowed identifying clear features attributable to a coherent superposition of dimer electronic and vibrational states contributing to the coherent electronic charge beating between the donor and the acceptor. The frequency of electronic charge beating is found to be 970 cm^-1 (34 fs) and can be observed for 150 fs. Through the support of high level ab initio TD-DFT computations of the entire dimer, we established that the vibrational modes preferentially optically accessed do not drive subsequent coupling between the electronic states on the 600 fs of the experiment. It was thereby possible to characterize the time scales of the early time femtosecond dynamics of the electronic coherence built by the optical excitation in a large rigid supramolecular system at a room temperature in solution.

2D-IR Spectroscopy Shows that Optimized DNA Minor Groove Binding of Hoechst33258 Follows an Induced Fit Model

The induced fit binding model describes a conformational change occurring when a small molecule binds to its biomacromolecular target. The result is enhanced noncovalent interactions between the ligand and biomolecule. Induced fit is well-established for small molecule-protein interactions, but its relevance to small molecule-DNA binding is less clear. We investigate the molecular determinants of Hoechst33258 binding to its preferred A-tract sequence relative to a suboptimal alternating A-T sequence. Results from two-dimensional infrared spectroscopy, which is sensitive to H-bonding and molecular structure changes, show that Hoechst33258 binding results in loss of the minor groove spine of hydration in both sequences, but an additional perturbation of the base propeller twists occurs in the A-tract binding region. This induced fit maximizes favorable ligand-DNA enthalpic contributions in the optimal binding case and demonstrates that controlling the molecular details that induce subtle changes in DNA structure may hold the key to designing next-generation DNA-binding molecules.

A Cu(ii) complex targeting the translocator protein: in vitro and in vivo antitumor potential and mechanistic insights

A Cu complex targeting the translocator protein induces a 98% reduction of tumor mass in a murine tumor model.

Modeling and Analysis of Intercalant Effects on Circular DNA Conformation

The large-scale conformation of DNA molecules plays a critical role in many basic elements of cellular functionality and viability. By targeting the structural properties of DNA, many cancer drugs, such as anthracyclines, effectively inhibit tumor growth but can also produce dangerous side effects. To enhance the development of innovative medications, rapid screening of structural changes to DNA can provide important insight into their mechanism of interaction. In this study, we report changes to circular DNA conformation from intercalation with ethidium bromide using all-atom molecular dynamics simulations and characterized experimentally by translocation through a silicon nitride solid-state nanopore. Our measurements reveal three distinct current blockade levels and a 6-fold increase in translocation times for ethidium bromide-treated circular DNA as compared to untreated circular DNA. We attribute these increases to changes in the supercoiled configuration hypothesized to be branched or looped structures formed in the circular DNA molecule. Further evidence of the conformational changes is demonstrated by qualitative atomic force microscopy analysis. These results expand the current methodology for predicting and characterizing DNA tertiary structure and advance nanopore technology as a platform for deciphering structural changes of other important biomolecules.

Force measurements reveal how small binders perturb the dissociation mechanisms of DNA duplex sequences

The force-driven separation of double-stranded DNA is crucial to the accomplishment of cellular processes like genome transactions. Ligands binding to short DNA sequences can have a local stabilizing or destabilizing effect and thus severely affect these processes. Although the design of ligands that bind to specific sequences is a field of intense research with promising biomedical applications, so far, their effect on the force-induced strand separation has remained elusive. Here, by means of AFM-based single molecule force spectroscopy, we show the co-existence of two different mechanisms for the separation of a short DNA duplex and demonstrate how they are perturbed by small binders. With the support of Molecular Dynamics simulations, we evidence that above a critical pulling rate one of the dissociation pathways becomes dominant, with a dramatic effect on the rupture forces. Around the critical threshold, we observe a drop of the most probable rupture forces for ligand-stabilized duplexes. Our results offer a deep understanding of how a stable DNA-ligand complex behaves under force-driven strand separation.

Structure and dynamics of proflavine association around DNA

Proflavine is a small molecule that intercalates into DNA and, thereby, acts as an anticancer agent. Intercalation of proflavine is shown to be a two-step process in which the first step is believed to be the formation of a pre-intercalative outside bound state. Experimental studies so far have been unable to capture the nature of the outside bound state. However, the sub-millisecond timescale observed in fluorescence kinetic experiments is often attributed to the binding of proflavine outside of DNA. Here, we have performed molecular dynamics simulations with multiple proflavine molecules to study the structure and dynamics of the formation of the outside bound state of DNA at different ion concentrations. We observed that the timescale of the outside bound state formation is, at least, five orders of magnitude faster (in nanoseconds) than the experimentally reported timescale (sub-milliseconds) attributed to binding outside DNA. Moreover, we also observed the stacked arrangement of proflavine all around DNA, which is different from the experimentally predicted stacking arrangement perpendicular to the helical axis of DNA in the close vicinity of the phosphate groups. This study, therefore, provides insight into the molecular structure and dynamics of the pre-intercalative outside bound state and will help in understanding the overall intercalation mechanism.

Understanding the kinetic mechanism of RNA single base pair formation

RNA functions are intrinsically tied to folding kinetics. The most elementary step in RNA folding is the closing and opening of a base pair. Understanding this elementary rate process is the basis for RNA folding kinetics studies. Previous studies mostly focused on the unfolding of base pairs. Here, based on a hybrid approach, we investigate the folding process at level of single base pairing/stacking. The study, which integrates molecular dynamics simulation, kinetic Monte Carlo simulation, and master equation methods, uncovers two alternative dominant pathways: Starting from the unfolded state, the nucleotide backbone first folds to the native conformation, followed by subsequent adjustment of the base conformation. During the base conformational rearrangement, the backbone either retains the native conformation or switches to nonnative conformations in order to lower the kinetic barrier for base rearrangement. The method enables quantification of kinetic partitioning among the different pathways. Moreover, the simulation reveals several intriguing ion binding/dissociation signatures for the conformational changes. Our approach may be useful for developing a base pair opening/closing rate model.

DNA Electronic Circular Dichroism on the Inter-Base Pair Scale: An Experimental-Theoretical Case Study of the AT Homo-Oligonucleotide

A successful elucidation of the near-ultraviolet electronic circular dichroism spectrum of a short double-stranded DNA is reported. Time-dependent density functional theory methods are shown to accurately predict spectra and assign bands on the microscopic base-pair scale, a finding that opens the field for using circular dichroism spectroscopy as a sensitive nanoscale probe of DNA to reveal its complex interactions with the environment.