Hidden entropic contribution in the thermodynamics of molecular complexation

It has become an axiom that the thermodynamic analysis of non-covalent molecular complexation is intrinsically model-dependent, i.e. the set of implicitly or explicitly introduced assumptions may strongly affect the thermodynamic parameters.

Dynamic hydrogen bonding and DNA flexibility in minor groove binders: molecular dynamics simulation of the polyamide f-ImPyIm bound to the Mlu1 (MCB) sequence 5'-ACGCGT-3' in 2:1 motif

Molecular dynamics simulations of the DNA 10-mer 5'-CCACGCGTGG-3' alone and complexed with the formamido-imidazole-pyrrole-imidazole (f-ImPyIm) polyamide minor groove binder in a 2:1 fashion were conducted for 50 ns using the pbsc0 parameters within the AMBER 12 software package. The change in DNA structure upon binding of f-ImPyIm was evaluated via minor groove width and depth, base pair parameters of Slide, Twist, Roll, Stretch, Stagger, Opening, Propeller, and x-displacement, dihedral angle distributions of ζ, ε, α, and γ determined using the Curves+ software program, and hydrogen bond formation. The dynamic hydrogen bonding between the f-ImPyIm and its cognate DNA sequence was compared to the static image used to predict sequence recognition by polyamide minor groove binders. Many of the predicted hydrogen bonds were present in less than 50% of the simulation; however, persistent hydrogen bonds between G5/15 and the formamido group of f-ImPyIm were observed. It was determined that the DNA is wider in the Complex than without the polyamide binder; however, there is flexibility in this particular sequence, even in the presence of the f-ImPyIm as evidenced by the range of minor groove widths the DNA exhibits and the dynamics of the hydrogen bonding that binds the two f-ImPyIm ions to the minor groove. The Complex consisting of the DNA and the 2 f-ImPyIm binders shows slight fraying of the 5' end of the 10-mer at the end of the simulation, but the portion of the oligomer responsible for recognition and binding is stable throughout the simulation. Several structural changes in the Complex indicate that minor groove binders may have a more active role in inhibiting transcription than just preventing binding of important transcription factors.

Recognition of the DNA minor groove by thiazotropsin analogues

Solution-phase self-association characteristics and DNA molecular-recognition properties are reported for three close analogues of minor-groove-binding ligands from the thiazotropsin class of lexitropsin molecules; they incorporate isopropyl thiazole as a lipophilic building block. Thiazotropsin B (AcImPy(iPr) ThDp) shows similar self-assembly characteristics to thiazotropsin A (FoPyPy(iPr) ThDp), although it is engineered, by incorporation of imidazole in place of N-methyl pyrrole, to swap its DNA recognition target from 5'-ACTAGT-3' to 5'-ACGCGT-3'. Replacement of the formamide head group in thiazotropsin A by nicotinamide in AIK-18/51 results in a measureable difference in solution-phase self-assembly character and substantially enhanced DNA association characteristics. The structures and associated thermodynamic parameters of self-assembled ligand aggregates and their complexes with their respective DNA targets are considered in the context of cluster targeting of DNA by minor-groove complexes.

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.

Complexation of DNA with ruthenium organometallic compounds: the high complexation ratio limit

Organometallic compounds possess two modes of interaction with DNA: intercalation at low complexation ratios and electrostatic adsorption at high ratios.

Thiazotropsin aggregation and its relationship to molecular recognition in the DNA minor groove

Aggregated states have been alluded to for many DNA minor groove binders but details of the molecule-on-molecule relationship have either been under-reported or ignored. Here we report our findings from ITC and NMR measurements carried out with AIK-18/51, a compound representative of the thiazotropsin class of DNA minor groove binders. The free aqueous form of AIK-18/51 is compared with that found in its complex with cognate DNA duplex d(CGACTAGTCG)2. Molecular self-association of AIK-18/51 is consistent with anti-parallel, face-to-face dimer formation, the building block on which the molecule aggregates. This underlying structure is closely allied to the form found in the ligand's DNA complex. NMR chemical shift and diffusion measurements yield a self-association constant Kass=(61±19)×10(3)M(-1) for AIK-18/51 that fits with a stepwise self-assembly model and is consistent with ITC data. The deconstructed energetics of this assembly process are reported with respect to a design strategy for ligand/DNA recognition.

Complexation of C60 fullerene with aromatic drugs

The contributions of various physical factors to the energetics of complexation of aromatic drug molecules with C(60) fullerene are investigated in terms of the calculated magnitudes of equilibrium complexation constants and the components of the net Gibbs free energy. Models of complexation are developed taking into account the polydisperse nature of fullerene solutions in terms of the continuous or discrete (fractal) aggregation of C(60) molecules. Analysis of the energetics has shown that stabilization of the ligand-fullerene complexes in aqueous solution is mainly determined by intermolecular van der Waals interactions and, to lesser extent, by hydrophobic interactions. The results provide a physicochemical basis for a potentially new biotechnological application of fullerenes as modulators of biological activity of aromatic drugs.