ABSTRACT
In 2004, we used NMR to solve the structure of the minor groove binder thiazotropsin A bound in a 2:1 complex to the DNA duplex, d(CGACTAGTCG)2. In this current work, we have combined theory and experiment to confirm the binding thermodynamics of this system. Molecular dynamics simulations that use polarizable or non-polarizable force fields with single and separate trajectory approaches have been used to explore complexation at the molecular level. We have shown that the binding process invokes large conformational changes in both the receptor and ligand, which is reflected by large adaptation energies. This is compensated for by the net binding free energy, which is enthalpy driven and entropically opposed. Such a conformational change upon binding directly impacts on how the process must be simulated in order to yield accurate results. Our MM-PBSA binding calculations from snapshots obtained from MD simulations of the polarizable force field using separate trajectories yield an absolute binding free energy (-15.4 kcal mol(-1)) very close to that determined by isothermal titration calorimetry (-10.2 kcal mol(-1)). Analysis of the major energy components reveals that favorable non-bonded van der Waals and electrostatic interactions contribute predominantly to the enthalpy term, whilst the unfavorable entropy appears to be driven by stabilization of the complex and the associated loss of conformational freedom. Our results have led to a deeper understanding of the nature of side-by-side minor groove ligand binding, which has significant implications for structure-based ligand development.
Subject(s)
DNA/chemistry , Thiazoles/chemistry , Calorimetry , Hydrogen Bonding , Ligands , Molecular Dynamics Simulation , ThermodynamicsABSTRACT
The development of new drugs to treat disease by binding directly to DNA offers much promise but is reliant on methods to determine the relative affinity of the putative drug for different DNA sequences. Such methods should ideally be rapid and inexpensive as well as reliable. Use of capillary electrophoresis in simple silica columns offers such a method. The development of systems in which the solvent carries a soluble polymer allows the reliable separation of DNA oligomers, of 12-20 bp in length, which can then be titrated with the ligand in competition experiments. The results obtained are comparable with those obtained by footprinting and give direct graphical output, easily analysed for relative binding affinity.
Subject(s)
DNA/chemistry , DNA/metabolism , Electrophoresis, Capillary/methods , Pharmaceutical Preparations/chemistry , Pharmaceutical Preparations/metabolism , Base Sequence , DNA Footprinting , Drug Interactions , Electrophoresis, Capillary/instrumentation , Ligands , Mass Spectrometry , Molecular Sequence DataABSTRACT
The idea that liquid water is not a uniform and random arrangement of molecules has been taken very seriously by the scientific community. Many experimental and computational investigations show that clathrate- or ice-like structures probably exist at a short time scale in solution. We have designed a new program to simulate water structure around solutes. Our model is based on the geometrical constraints of hydrogen bonding in order to be capable of producing clathrate-like structures. Simulations with small molecules and bio-molecules, using the new software, produce networks of water with specific patterns made of small water rings. The water structures built are consistent with the classification of molecules in terms of structure breaking and making. This approach may give insight into, and a more accurate description of, drug-receptor interactions. The results also suggest that water structure may impart sufficient energy to modify the conformational space of organic molecules through hydrogen bonding.
Subject(s)
Computer Simulation , DNA/chemistry , Hexanones/chemistry , Propanolamines/chemistry , Water/chemistry , Crystallography, X-Ray , Hydrogen Bonding , Methane/chemistry , Models, Molecular , Molecular Structure , SoftwareABSTRACT
The sequence selectivity of small molecules binding to the minor groove of DNA can be predicted by "in silico footprinting". Any potential ligand can be docked in the minor groove and then moved along it using simple simulation techniques. By applying a simple scoring function to the trajectory after energy minimization, the preferred binding site can be identified. We show application to all known noncovalent binding modes, namely 1:1 ligand:DNA binding (including hairpin ligands) and 2:1 side-by-side binding, with various DNA base pair sequences and show excellent agreement with experimental results from X-ray crystallography, NMR, and gel-based footprinting.
