Sequence-Specific Positions of Water Molecules at the Interface between DNA and Minor Groove Binders

Solvent effects on the NMR shieldings of stacked DNA base pairs

Stacking effects are among the most important effects in DNA. We have recently studied their influence in fragments of DNA through the analysis of NMR magnetic shieldings, firstly in vacuo. As a continuation of this line of research we show here the influence of solvent effects on the shieldings through the application of both explicit and implicit models. We found that the explicit solvent model is more appropriate for consideration due to the results matching better in general with experiments, as well as providing clear knowledge of the electronic origin of the value of the shieldings. Our study is grounded on a recently developed theoretical model of our own, by which we are able to learn about the magnetic effects of given fragments of DNA molecules on selected base pairs. We use the shieldings of the atoms of a central base pair (guanine-cytosine) of a selected fragment of DNA molecules as descriptors of physical effects, like π-stacking and solvent effects. They can be taken separately and altogether. The effect of π-stacking is introduced through the addition of some pairs above and below of the central base pair, and now, the solvent effect is considered including a network of water molecules that consist of two solvation layers, which were fixed in the calculations performed in all fragments. We show that the solvent effects enhance the stacking effects on the magnetic shieldings of atoms that belong to the external N-H bonds. The net effect is of deshielding on both atoms. There is also a deshielding effect on the carbon atoms that belong to CO bonds, for which the oxygen atom has an explicit hydrogen bond (HB) with a solvent water molecule. Solvent effects are found to be no higher than a few percent of the total value of the shieldings (between 1% and 5%) for most atoms, although there are few for which such an effect can be higher. There is one nitrogen atom, the acceptor of the HB between guanine and cytosine, that is more highly shielded (around 15 ppm or 10%) when the explicit solvent is considered. In a similar manner, the most external nitrogen atom of cytosine and the hydrogen atom that is bonded to it are highly deshielded (around 10 ppm for nitrogen and around 3 ppm for hydrogen).

Explicit solvation thermodynamics in ionic solution: extending grid inhomogeneous solvation theory to solvation free energy of salt-water mixtures

Hydration thermodynamics play a fundamental role in fields ranging from the pharmaceutical industry to environmental research. Numerous methods exist to predict solvation thermodynamics of compounds ranging from small molecules to large biomolecules. Arguably the most precise methods are those based on molecular dynamics (MD) simulations in explicit solvent. One theory that has seen increased use is inhomogeneous solvation theory (IST). However, while many applications require accurate description of salt-water mixtures, no implementation of IST is currently able to estimate solvation properties involving more than one solvent species. Here, we present an extension to grid inhomogeneous solvation theory (GIST) that can take salt contributions into account. At the example of carbazole in 1 M NaCl solution, we compute the solvation energy as well as first and second order entropies. While the effect of the first order ion entropy is small, both the water-water and water-ion entropies contribute strongly. We show that the water-ion entropies are efficiently approximated using the Kirkwood superposition approximation. However, this approach cannot be applied to the water-water entropy. Furthermore, we test the quantitative validity of our method by computing salting-out coefficients and comparing them to experimental data. We find a good correlation to experimental salting-out constants, while the absolute values are overpredicted due to the approximate second order entropy. Since ions are frequently used in MD, either to neutralize the system or as a part of the investigated process, our method greatly extends the applicability of GIST. The use-cases range from biopharmaceuticals, where many assays require high salt concentrations, to environmental research, where solubility in sea water is important to model the fate of organic substances.

Calcium Ions Affect Water Molecular Structures Surrounding an Oligonucleotide Duplex as Revealed by Sum Frequency Generation Vibrational Spectroscopy

The solvation of DNA in water facilitates the formation of a hydration layer surrounding it, thus stabilizing the DNA duplex in the biological aqueous environment. In this study, via using the lipid bilayer as a soft substrate to accommodate the duplex oligonucleotide, the structure of the water layer surrounding the oligonucleotide was detected under the perturbation of the calcium ions (Ca²⁺) with chiral and achiral sum frequency generation (SFG) vibrational spectroscopy. With increasing Ca²⁺ concentration, both the chiral and achiral water vibrational signals had similar concentration-dependent changes, i.e., an initial decreasing phase followed by an increasing phase. However, when the Ca²⁺ concentrations were adjusted to within the range comparable to those in the human serum, the chiral water vibrational signals remained nearly unchanged, whereas the achiral water vibrational signals still changed as a function of the Ca²⁺ concentration. Therefore, the current experimental result supports the possible protection function of the chiral hydration layer against the Ca²⁺ ions, which generally exist in the cell sap and play important roles in many biological functions.

Bound Compound, Interfacial Water, and Phenyl Ring Rotation Dynamics of a Compound in the DNA Minor Groove

DB2277, a heterocyclic diamidine, is a successful design for mixed base pair (bp) DNA sequence recognition. The compound has a central aza-benzimidazole group that forms two H-bonds with a GC bp that has flanking AT bps. The nuclear magnetic resonance structure of the DB2277-DNA complex with an AAGATA recognition site sequence was determined, and here we report extended molecular dynamics (MD) simulations of the structure. DB2277 has two terminal phenyl-amidine groups, one of which is directly linked to the DB2277 heterocyclic core and the other through a flexible -OCH₂- group. The flexibly linked phenyl is too far from the minor groove floor to make direct H-bonds but is linked to an AT bp through water-mediated H-bonds. The flexibly linked phenyl-amidine with water-mediated H-bonds to the bases at the floor of the minor groove suggested that it might rotate in time spans accessible in MD. To test this idea, we conducted multimicrosecond MD simulations to determine if these phenyl rotations could be observed for a bound compound. In a 3 μs simulation, highly dynamic torsional motions were observed for the -OCH₂-linked phenyl but not for the other phenyl. The dynamics periodically reached a level to allow 180° rotation of the phenyl while it was still bound in the minor groove. This is the first observation of rotation of a phenyl bound to DNA, and the results provide mechanistic details about how a rotation can occur as well as how mixed bp recognition can occur for monomer compounds bound to the minor groove.

Conformational rearrangement of 1,2-d(GG) intrastrand cis-diammineplatinum crosslinked DNA is driven by counter-ion penetration within the minor groove of the modified site

The major structural aberrations of DNA induced by a cis-diammineplatinum (II) 1,2-d(GG) intrastrand cross-link (CPT) have been known for decades. To gain deeper insights into the structural dynamics of the sequence-dependent DNA distortions adjacent to the CPT adduct, we employed molecular modeling and molecular dynamics (MD) simulations. The structural dynamics of native (N-DNA) and cisPt 1,2-d(GG) crosslinked (CPT-DNA) in the form of symmetric 36 nt d(G₂T₁₅G*G*T₁₅G₂)●C₂A₁₅CCA₁₅C₂) oligonucleotide duplexes is compared. The selected sequence context enabled tracking of the origin of the DNA axis curvature at the YpR flexible points (N-DNA), the enhancement of axis bending, and further distortions due to steric/electrostatic perturbations arising from the CPT-crosslink. In addition to the known structural distortions of CPT-DNA: helix bend towards the major groove; local helix unwinding; high roll angle between cross-linked guanine bases; and adoption of A-form DNA on the 5'-side of the CPT-crosslink (TpG junction); our results show the existence of a singular irreversible and reproducible conformational rearrangement, not previously observed, resulting in two stable CPT-DNA1 and CPT-DNA2 conformers. The CPT-DNA2 conformation presents an enhanced DNA axis bend and a wider and shallower minor grove with increased solvent accessibility within the modified site. It is concluded that the polymorphous (unstable) DNA environment near the cisPt 1,2-d(GG) unit in synergy with specific dynamic events, such as prolonged minor groove retention of particular Na⁺ ions and water redistribution within the d(TG*G*T) site, together with the formation of extra and more stable H-bonds between Pt(NH₃)₂ amines and neighboring nucleotides, are cooperatively responsible for the initiation of the conformational rearrangement leading to the CPT-DNA2 conformer, which, surprisingly, closely resembles the HMGB1-bound CPT-DNA structure. Graphical abstract Superimposed averaged structures of normal (N-DNA, green) and cisplatin intrastrand cross-linked (CPT-DNA, orange). Global DNA axes: N-DNA (blue beads); CPT-DNA (red beads); PT (yellow sphere); crosslinked dGs viewed from the minor groove (bold).

DNA's Chiral Spine of Hydration

The iconic helical structure of DNA is stabilized by the solvation environment, where a change in the hydration state can lead to dramatic changes to the DNA structure. X-ray diffraction experiments at cryogenic temperatures have shown crystallographic water molecules in the minor groove of DNA, which has led to the notion of a spine of hydration of DNA. Here, chiral nonlinear vibrational spectroscopy of two DNA sequences shows that not only do such structural water molecules exist in solution at ambient conditions but that they form a chiral superstructure: a chiral spine of hydration. This is the first observation of a chiral water superstructure templated by a biomolecule. While the biological relevance of a chiral spine of hydration is unknown, the method provides a direct way to interrogate the properties of the hydration environment of DNA and water in biological settings without the use of labels.

Stimuli-responsive colorimetric and NIR fluorescence combination probe for selective reporting of cellular hydrogen peroxide

Hydrogen peroxide (H2O2) is a key reactive oxygen species and a messenger in cellular signal transduction apart from playing a vital role in many biological processes in living organisms. In this article, we present phenyl boronic acid-functionalized quinone-cyanine (QCy-BA) in combination with AT-rich DNA (exogenous or endogenous cellular DNA), i.e., QCy-BA⊂DNA as a stimuli-responsive NIR fluorescence probe for measuring in vitro levels of H2O2. In response to cellular H2O2 stimulus, QCy-BA converts into QCy-DT, a one-donor-two-acceptor (D2A) system that exhibits switch-on NIR fluorescence upon binding to the DNA minor groove. Fluorescence studies on the combination probe QCy-BA⊂DNA showed strong NIR fluorescence selectively in the presence of H2O2. Furthermore, glucose oxidase (GOx) assay confirmed the high efficiency of the combination probe QCy-BA⊂DNA for probing H2O2 generated in situ through GOx-mediated glucose oxidation. Quantitative analysis through fluorescence plate reader, flow cytometry and live imaging approaches showed that QCy-BA is a promising probe to detect the normal as well as elevated levels of H2O2 produced by EGF/Nox pathways and post-genotoxic stress in both primary and senescent cells. Overall, QCy-BA, in combination with exogenous or cellular DNA, is a versatile probe to quantify and image H2O2 in normal and disease-associated cells.

Sequence-specific recognition of DNA minor groove by an NIR-fluorescence switch-on probe and its potential applications

In molecular biology, understanding the functional and structural aspects of DNA requires sequence-specific DNA binding probes. Especially, sequence-specific fluorescence probes offer the advantage of real-time monitoring of the conformational and structural reorganization of DNA in living cells. Herein, we designed a new class of D2A (one-donor-two-acceptor) near-infrared (NIR) fluorescence switch-on probe named quinone cyanine–dithiazole (QCy–DT) based on the distinctive internal charge transfer (ICT) process for minor groove recognition of AT-rich DNA. Interestingly, QCy–DT exhibited strong NIR-fluorescence enhancement in the presence of AT-rich DNA compared to GC-rich and single-stranded DNAs. We show sequence-specific minor groove recognition of QCy–DT for DNA containing 5′-AATT-3′ sequence over other variable (A/T)4 sequences and local nucleobase variation study around the 5′-X(AATT)Y-3′ recognition sequence revealed that X = A and Y = T are the most preferable nucleobases. The live cell imaging studies confirmed mammalian cell permeability, low-toxicity and selective staining capacity of nuclear DNA without requiring RNase treatment. Further, Plasmodium falciparum with an AT-rich genome showed specific uptake with a reasonably low IC₅₀ value (<4 µM). The ease of synthesis, large Stokes shift, sequence-specific DNA minor groove recognition with switch-on NIR-fluorescence, photostability and parasite staining with low IC₅₀ make QCy–DT a potential and commercially viable DNA probe.

Polyamide-scorpion cyclam lexitropsins selectively bind AT-rich DNA independently of the nature of the coordinated metal

Cyclam was attached to 1-, 2- and 3-pyrrole lexitropsins for the first time through a synthetically facile copper-catalyzed "click" reaction. The corresponding copper and zinc complexes were synthesized and characterized. The ligand and its complexes bound AT-rich DNA selectively over GC-rich DNA, and the thermodynamic profile of the binding was evaluated by isothermal titration calorimetry. The metal, encapsulated in a scorpion azamacrocyclic complex, did not affect the binding, which was dominated by the organic tail.

Hydrogen-bonding patterns of minor groove-binder-DNA complexes reveal criteria for discovery of new scaffolds

Minor groove-binding ligands are able to control gene expression and are of great interest for therapeutic applications. We extracted hydrogen-bonding geometries from all available structures of minor groove-binder-DNA complexes of two noncovalent binding modes, namely 1:1 (including hairpin and cyclic ligands) and 2:1 ligand/DNA binding. Positions of the ligand atoms involved in hydrogen bonding deviate from idealized hydrogen bond geometries and do not exploit the possibilities indicated by water molecules. Therefore, we suggest the inclusion of shape-based descriptors rather than hydrogen-bond patterns in virtual screening protocols for the identification of innovative minor groove-binding scaffolds.