Solvation Dynamics of Hoechst 33258 in Water: An Equilibrium and Nonequilibrium Molecular Dynamics Study

Molecular Mechanism of Ligand Binding to the Minor Groove of DNA

Although DNA-ligand binding is pervasive in biology, little is known about molecular-level binding mechanisms. Using all-atom, explicit-solvent molecular dynamics simulations in conjunction with weighted ensemble (WE)-enhanced sampling, an ensemble of 2562 binding trajectories of Hoechst 33258 (H33258) to d(CGC AAA TTT GCG) was generated from which the binding mechanism was extracted. In particular, the electrostatic interaction between the positively charged H33258 and the negatively charged DNA backbone drives the formation of initial H33258-DNA contacts. After this initial contact, a hinge-like intermediate state is formed in which one end of H33258 inserts into the minor groove of DNA. Following hinge state formation is a concerted motion whereby the second end of H33258 swings into the minor groove and the spine of hydration along the minor groove causing dehydration. This study illustrates how WE-enhanced simulations of biomolecular ligation processes can offer novel mechanistic insights by generating ensembles of binding events.

Explicit Projection of Stokes Shifts onto Solvent Motion in an Aqueous Liquid and Linear Response Theory

We investigate the molecular origin of the fluorescence Stokes shift in an aqueous liquid. By examining the speed of energy change, the solvation response function is explicitly projected onto the translational and rotational motions of water molecules for both nonequilibrium relaxation and equilibrium fluctuations. Molecular dynamics simulations of a tryptophan solution show that these two processes have highly consistent dynamics, not only for the total response function but also for the decomposed components in terms of specific molecular movements. We found that the rotational mode governs the relaxation of the Stokes shift, whereas the translational mode contributes non-negligibly with slower dynamics. This consistency implies the similarity of the underlying translational and rotational movements of water molecules as the system is far away from and at equilibrium, supporting the validity of the linear response theory at the molecular level. The decomposition methodology is also applicable to a rigid solvent.

To unravel the connection between the non-equilibrium and equilibrium solvation dynamics of tryptophan: success and failure of the linear response theory of fluorescence Stokes shift

The connections between the non-equilibrium solvation dynamics upon optical transitions and the system's equilibrium fluctuations are explored in aqueous liquid. Linear response theory correlates time-dependent fluorescence with the equilibrium time correlation functions. In the previous work [T. Li, J. Chem. Theory Comput., 2017, 13, 1867], Stokes shift was explicitly decomposed into the contributions of various order time correlation functions on the excited state surface. Gaussian fluctuations of the solute-solvent interactions validate linear response theory. Correspondingly, the deviation of the Gaussian statistics causes the inefficiency of linear response evaluation. The above mechanism is thoroughly tested in this study. By employing molecular simulations, multiple non-equilibrium processes, not necessarily initiated from the ground state equilibrium minimum, were examined for tryptophan. Both the success and failure of linear response theory are found for this simple system and the mechanism is analyzed. These observations, assisted by the width dynamics, the initial state linear response approach, and the variation of the solvation structures, integrally verify the virtue of the excited state Gaussian statistics on the dynamics of Stokes shift.

Linear Response Theory for Decomposition Energies of Stokes Shift in Proteins

In aqueous solution, fluorescence Stokes shift experiments monitor the relaxation of the solute-solvent interactions upon photon excitation of the solute chromophore. Linear response (LR) theory expects the identical dynamics between the Stokes shift and the system's spontaneous fluctuations. Whether this identity guarantees similar dynamics between the nonequilibrated and equilibrium processes for the decomposition energy of the Stokes shift is the main focus of this study. In our previous work [Li, T. J. Chem. Theory Comput. 2017, 13, 1867-1873], Stokes shift is properly correlated with various order time-correlation functions. As a continuation, its decomposition energy from the subsystem is further represented as the full summation of all of the cross-time correlation functions between the decomposition energy and the total solute-solvent interactions. Gaussian statistics of the total solute-solvent interactions ensure the same decay rates among the odd orders not only for the time-correlation functions but also for the cross-time correlation functions, validating the LR of the Stokes shift and the decompositions, respectively. The above mechanism is verified by molecular dynamics simulations in the protein Staphylococcus nuclease and is robust even as the decomposed energy associated with an individual residue exhibits typical non-Gaussian properties. Further examinations reveal the consistent molecular motions for a specific residue over the nonequilibrium and equilibrium processes, which are responsible for the nonequilibrium dynamics of the associated decomposed energy. Our results show the appropriateness of LR on finer molecular scales.

On the validity of linear response approximations regarding the solvation dynamics of polyatomic solutes

Gaussian statistics and linear response predictions of the nonequilibrium solvation dynamics are tested for numerous solute/solvent combinations.

Dispersed dynamics of solvation in G-quadruplex DNA: comparison of dynamic Stokes shifts of probes in parallel and antiparallel quadruplex structures

G-quadruplex DNA (GqDNA) structures play an important role in many specific cellular functions and are promising anti-tumor targets for small molecules (ligands). Here, we measured the dynamic Stokes shift of a ligand (Hoechst) bound to parallel c-Myc (mPu22) GqDNA over five decades of time from 100 fs to 10 ns, and compared it with the previously reported dynamics of DAPI bound to antiparallel human telomeric (hTelo22) GqDNA (Pal et al 2015 J. Phys. Chem. Lett. 6 1754). Stokes shift data from fluorescence up-conversion and time-correlated single photon counting experiments was combined to cover the broad dynamic range. The results show that the solvation dynamics of Hoechst in parallel mPu22 GqDNA follow a power law relaxation, added to fast 2 ps exponential relaxation, from 100 fs to 10 ns, with only a subtle difference of power law exponents in the two ligand-GqDNA systems (0.06 in Hoechst-mPu22 compared to 0.16 in DAPI-hTelo22). We measured steady-state fluorescence spectra and time-resolved anisotropy decays which confirm the tight binding of Hoechst to parallel mPu22 with a binding constant of ~1  ×  10⁵ M^-1. The molecular docking of Hoechst in parallel GqDNA followed by a 50 ns molecular dynamics (MD) simulation on a Hoechst-GqDNA complex reveals that Hoechst binds to one of the outer G-tetrads by end-stacking near G13 and G4, which is different from the binding site of DAPI inside a groove of antiparallel hTelo22 GqDNA. Reconciling previous experimental and simulation results, we assign the 2 ps component to the hydration dynamics of only weakly perturbed water near mPu22 and the power law relaxation to the coupled motion of water and DNA (i.e. DNA backbone, unpaired bases and loops connecting G-tetrads) which come near the Hoechst inside parallel GqDNA.

Diabatic Population Matrix Formalism for Performing Molecular Mechanics Style Simulations with Multiple Electronic States

An accurate description of nonbonded interactions is important in investigating dynamics of molecular systems. In many situations, fixed point charge models are successfully applied to explaining various chemical phenomena. However, these models with conventional formulations will not be appropriate in elucidating the detailed dynamics during nonadiabatic events. This is mainly because the chemical properties of any molecule, especially its electronic populations, significantly change with respect to molecular distortions in the vicinity of the surface crossing. To overcome this issue in molecular simulations yet within the framework of the fixed point charge model, we define a diabatic electronic population matrix and substitute it for the conventional adiabatic partial charges. We show that this matrix can be readily utilized toward attaining more reliable descriptions of Coulombic interactions, in combination with the interpolation formalism for obtaining the intramolecular interaction potential. We demonstrate how the mixed formalism with the diabatic charges and the interpolation can be applied to molecular simulations by conducting adiabatic and nonadiabatic molecular dynamics trajectory calculations of the green fluorescent protein chromophore anion in aqueous environment.

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

Nucleic acids are flexible molecules and their dynamical properties play a key role in molecular recognition events. Small binders interacting with DNA fragments induce both structural and dynamical changes in the double helix. We study the dynamics of a DNA dodecamer and of its complexes with Hoechst 33258, which is a minor groove binder, and with the ethidium cation, which is an intercalator, by molecular dynamics simulation. The thermodynamics of DNA-drug interaction is evaluated in connection with the structure and the dynamics of the resulting complexes. We identify and characterize the relevant changes in the configurational distribution of the DNA helix and relate them to the corresponding entropic contributions to the binding free energy. The binder Hoechst locks the breathing motion of the minor groove inducing a reduction of the configurational entropy of the helix, which amounts to 20 kcal mol(-1). In contrast, intercalations with the ethidium cation enhance the flexibility of the double helix. We show that the balance between the energy required to deform the helix for the intercalation and the gain in configurational entropy is the origin of cooperativity in the binding of a second ethidium and of anti-cooperativity in the binding of a third one. The results of our study provide an understanding of the relation between structure, dynamics and energetics in the interaction between DNA fragments and small binders, highlighting the role of dynamical changes and consequent variation of the configurational entropy of the DNA double helix for both types of binders.

Molecular dynamics investigation of the vibrational spectroscopy of isolated water in an ionic liquid

Experimental studies examining the structure and dynamics of water in ionic liquids (ILs) have revealed local ion rearrangements that occur an order of magnitude faster than complete randomization of the liquid structure. Simulations of an isolated water molecule embedded in 1-butyl-3-methyl imidazolium hexafluorophosphate, [bmim][PF6], were performed to shed insight into the nature of these coupled water-ion dynamics. The theoretical calculations of the spectral diffusion dynamics and the infrared absorption spectra of the OD stretch of isolated HOD in [bmim][PF6] agree well with experiment. The infrared absorption line shape of the OD stretch is narrower and blue-shifted in the IL compared to those in aqueous solution. Decomposition of the OD frequency time correlation function revealed that translational motions of the anions dominate the spectral diffusion dynamics.

Translational repression of thymidylate synthase by targeting its mRNA

Resistance to drugs targeting human thymidylate synthase (TS) poses a major challenge in the field of anti-cancer therapeutics. Overexpression of the TS protein has been implicated as one of the factors leading to the development of resistance. Therefore, repressing translation by targeting the TS mRNA could help to overcome this problem. In this study, we report that the compound Hoechst 33258 (HT) can reduce cellular TS protein levels without altering TS mRNA levels, suggesting that it modulates TS expression at the translation level. We have combined nuclear magnetic resonance, UV-visible and fluorescence spectroscopy methods with docking and molecular dynamics simulations to study the interaction of HT with a region in the TS mRNA. The interaction predominantly involves intercalation of HT at a CC mismatch in the region near the translational initiation site. Our results support the use of HT-like compounds to guide the design of therapeutic agents targeting TS mRNA.

Microscopic rearrangement of bound minor groove binders detected by NMR

Thermodynamic and structural studies are commonly utilized to optimize small molecules for specific DNA interactions, and, thus, a significant amount of binding data is available. However, the dynamic processes that are involved in minor groove complex formation and maintenance are not fully understood. To help define the processes involved, we have conducted 1D and 2D NMR in conjunction with biosensor-SPR experiments with a variety of compounds and symmetric, as well as asymmetric, AT tract DNA sequences. Surprisingly, the NMR data clearly show exchange between equivalent binding sites for strongly binding compounds like netropsin and DB921 (Ka > 10(8) M(-1)) that does not involve dissociation off the DNA. A quantitative analysis of the data revealed that these bound exchange rates are indeed much faster than the macroscopic dissociation rates which were independently determined by biosensor-SPR. Additionally, we could show the existence of at least two 1:1 compound DNA complexes at the same site for the interaction of these compounds with an asymmetric DNA sequence. To explain this behavior we introduced a model in which the ligand is rapidly flipping between two orientations while in close association with the DNA. The ligand reorientation will contribute favorably to the binding entropy. As the potential of minor groove binders to form more than a single complex with asymmetric, as well as symmetric, duplexes is widely unknown, the consequences for binding thermodynamics and compound design are discussed.

Multiple drug delivery hydrogel system for spinal cord injury repair strategies

The multifactorial pathological progress of spinal cord injury (SCI) is probably the main reason behind the absence of efficient therapeutic approaches. Hence, very recent highlights suggest the use of new multidrug delivery systems capable of local controlled release of therapeutic agents. In this work, a biocompatible hydrogel-based system was developed as multiple drug delivery tool, specifically designed for SCI repair strategies. Multiple release profiles were achieved by loading gel with a combination of low and high steric hindrance molecules. In vitro, in vivo and ex vivo release studies showed an independent combination of fast diffusion-controlled kinetics for smaller molecules together with slow diffusion-controlled kinetics for bigger ones. A preserved functionality of loaded substances was always achieved, confirming the absence of any chemical stable interactions between gel matrix and loaded molecules. Moreover, the relevant effect of the cerebrospinal fluid flux dynamics on the drug diffusion in the spinal cord tissue was here revealed for the first time: an oriented delivery of the released molecules in the spinal cord tract caudally to the gel site is demonstrated, thus suggesting a more efficient gel positioning rostrally to the lesion.

Condensed phase molecular dynamics using interpolated potential energy surfaces with application to the resolvation process of coumarin 153

Interpolated potential energy surfaces (PESs) have been used for performing reliable molecular dynamics (MD) simulations of small molecular reactions. In this article, we extend this method to MD simulations in condensed phase and show that the same scheme can also be feasibly used when it is supplemented with additional terms for describing intermolecular interactions. We then apply the approach for studying the resolvation process of coumarin 153 in a number of polar solvents. We find that the interpolated surface actually reproduces experimentally found features much better than the conventional force field based potential especially in terms of both dynamics Stokes shift in the short time limit and solute vibrational decoherence. This shows that the solute vibrational effect is important to some degree along the resolvation and should be modeled properly for accurate description of the related dynamics. The stability issue of trajectories on the interpolated PESs is also discussed, in regard to the goal of reliably performing long time simulations. Operational limitations of the present scheme are also discussed.

Effects of an unnatural base pair replacement on the structure and dynamics of DNA and neighboring water and ions

Incorporating small molecule probes into biomolecular systems to report on local structure and dynamics is a powerful strategy that underlies a wide variety of experimental techniques, including fluorescence, electron paramagnetic resonance (EPR), and Forster resonance energy transfer (FRET) measurements. When an unnatural probe is inserted into a protein or DNA, the degree to which the presence of the probe has perturbed the local structure and dynamics it was intended to study is always an important concern. Here, molecular dynamics (MD) simulations are used to systematically study the effect of replacing a DNA base pair with a fluorescent probe, coumarin 102 deoxyriboside, at six unique sites along an A-tract DNA dodecamer. While the overall structure of the DNA oligonucleotide remains intact, replacement of A*T base pairs leads to widespread structural and dynamic perturbations up to four base pairs away from the probe site, including widening of the minor groove and increased DNA flexibility. New DNA conformations, not observed in the native sequence, are sometimes found in the vicinity of the probe and its partner abasic site analog. Strong correlations are demonstrated between DNA surface topology and water mobility.