Sequence-Dependent Configurational Entropy Change of DNA upon Intercalation

On the recognition of Yersinia protein tyrosine phosphatase by carboxylic acid derivatives

Some members of Yersinia (Y), a genus of bacteria in the family Yersiniaceae, are pathogenic in humans, causing a range of health problems, from gastrointestinal syndromes to the plague. The Y protein tyrosine phosphatase (PTP) YopH is a crucial virulence determinant, considering the vital roles of PTPs in the intracellular signal transduction pathways and cell cycle control. The structural understanding of YopH as a cellular target in pathogenic conditions caused by Y infection is a prerequisite for designing potent and selective YopH inhibitors. Thus, by using molecular docking simulations, the open and closed conformations of the so-called 'WPD loop' (352-Gly-Asn-Trp-Pro-Asp-Gln-Thr-Ala-Val-Ser-361), located nearby the active site (403-Cys-Arg-Ala-Gly-Val-Gly-Arg-Thr-410) in YopH structure, are shown to be relevant for recognition by carboxylic acid derivatives, and the closed conformation is a more preferable receptor in terms of the quantitative correlation with experimental data. In both cases, aurintricarboxylic acid (ATA) has the greatest affinity to YopH. Consequently, a quantum mechanics/molecular mechanics (QM/MM) molecular model is derived to see into the extent of the ATA-induced open-closed conformational change. Active site residues and the WPD loop, as well as ATA are treated using SCC-DFTB-D (QM level), while the rest of the complex is treated using AMBER force field (MM level). The active/inactive functional behavior of YopH is explored by observing the interaction mode of ATA with the wild-type (wt)/Cys403Ser receptor and evaluating the competitive inhibition parameters. Implications of the present study for experimental research are discussed. Communicated by Ramaswamy H. Sarma.

Comparative evaluation of 11 in silico models for the prediction of small molecule mutagenicity: role of steric hindrance and electron-withdrawing groups

The goal of this investigation was to perform a comparative analysis on how accurately 11 routinely-used in silico programs correctly predicted the mutagenicity of test compounds that contained either bulky or electron-withdrawing substituents. To our knowledge this is the first study of its kind in the literature. Such substituents are common in many pharmaceutical agents so there is a significant need for reliable in silico programs to predict precisely whether they truly pose a risk for mutagenicity. The predictions from each program were compared to experimental data derived from the Ames II test, a rapid reverse mutagenicity assay with a high degree of agreement with the traditional Ames assay. Eleven in silico programs were evaluated and compared: Derek for Windows, Derek Nexus, Leadscope Model Applier (LSMA), LSMA featuring the in vitro microbial Escherichia coli-Salmonella typhimurium TA102 A-T Suite (LSMA+), TOPKAT, CAESAR, TEST, ChemSilico (±S9 suites), MC4PC and a novel DNA docking model. The presence of bulky or electron-withdrawing functional groups in the vicinity of a mutagenic toxicophore in the test compounds clearly affected the ability of each in silico model to predict non-mutagenicity correctly. This was because of an over reliance on the part of the programs to provide mutagenicity alerts when a particular toxicophore is present irrespective of the structural environment surrounding the toxicophore. From this investigation it can be concluded that these models provide a high degree of specificity (ranging from 71% to 100%) and are generally conservative in their predictions in terms of sensitivity (ranging from 5% t o 78%). These values are in general agreement with most other comparative studies in the literature. Interestingly, the DNA docking model was the most sensitive model evaluated, suggesting a potentially useful new mode of screening for mutagens. Another important finding was that the combination of a quantitative structure-activity relationship and an expert rules system appeared to offer little advantage in terms of sensitivity, despite of the requirement for such a screening paradigm under the ICH M7 regulatory guideline.

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.

Spatial assignment of symmetry adapted perturbation theory interaction energy components: The atomic SAPT partition

We develop a physically-motivated assignment of symmetry adapted perturbation theory for intermolecular interactions (SAPT) into atom-pairwise contributions (the A-SAPT partition). The basic precept of A-SAPT is that the many-body interaction energy components are computed normally under the formalism of SAPT, following which a spatially-localized two-body quasiparticle interaction is extracted from the many-body interaction terms. For electrostatics and induction source terms, the relevant quasiparticles are atoms, which are obtained in this work through the iterative stockholder analysis (ISA) procedure. For the exchange, induction response, and dispersion terms, the relevant quasiparticles are local occupied orbitals, which are obtained in this work through the Pipek-Mezey procedure. The local orbital atomic charges obtained from ISA additionally allow the terms involving local orbitals to be assigned in an atom-pairwise manner. Further summation over the atoms of one or the other monomer allows for a chemically intuitive visualization of the contribution of each atom and interaction component to the overall noncovalent interaction strength. Herein, we present the intuitive development and mathematical form for A-SAPT applied in the SAPT0 approximation (the A-SAPT0 partition). We also provide an efficient series of algorithms for the computation of the A-SAPT0 partition with essentially the same computational cost as the corresponding SAPT0 decomposition. We probe the sensitivity of the A-SAPT0 partition to the ISA grid and convergence parameter, orbital localization metric, and induction coupling treatment, and recommend a set of practical choices which closes the definition of the A-SAPT0 partition. We demonstrate the utility and computational tractability of the A-SAPT0 partition in the context of side-on cation-π interactions and the intercalation of DNA by proflavine. A-SAPT0 clearly shows the key processes in these complicated noncovalent interactions, in systems with up to 220 atoms and 2845 basis functions.

Sequence-dependent binding of flavonoids to duplex DNA

The whole family of structurally distinct flavonoids has been recognized as a valuable source of prospective anticancer agents. There is experimental evidence demonstrating that some flavonoids, like flavopiridol (FLP) and quercetin (QUE), bind to DNA influencing their key physiological function. FLP is involved in the combined mode of interaction (intercalation and minor groove binding), while QUE is viewed as a minor groove binder. From a physical standpoint, experimental and theoretical studies have not so far provided a sufficiently consistent picture of the nature of interaction with DNA. Herein the sequence-dependent binding of FLP and of QUE (two representative examples of the structurally different flavonoids) with duplex DNA, containing a variety of the sequences of eight nucleotides (I: GGGGCCCC, II: GGCCGGCC, III: AAAATTTT, IV: AAGCGCTT, V: GCGCGCGC) in the 5'-strand, is investigated using a sophisticated molecular dynamics (MD) approach. For various parts (helix, backbone, bases) of the DNA structure, the change of asymptotic (in terms of an infinite length of MD simulation) configurational entropy, being the thermodynamic consequence of DNA flexibility change due to ligand binding, is explored. As far as the sequence-dependent extent of DNA flexibility change upon QUE (or FLP) binding is concerned, for the entire double helix, increased flexibility is observed for I (or I ≈ II), while increased rigidity is found to be in the order of V > III > II > IV (or III > V > IV) for the rest of sequences. For the backbone, increased rigidity in the order of V > III > II > IV > I (or III > V > IV > I > II) is generally observed. For the nucleobases, increased flexibility is determined for I and II (I > II for both ligands), while increased rigidity in the order of V ≈ III > IV (or III > V > IV) is reported for the other sequences. Of the overall increased rigidity of the DNA structure upon ligand binding that is observed for the sequences III, IV, and V, about 50-70% comes from the sugar-phosphate backbone. Noteworthy is that the increased flexibility of the entire double helix and of the complete system of nucleobases upon ligand binding is only established for sequence I. The insights are further subtly substantiated by considering the configurational entropy contributions at the level of individual nucleobase pairs and of individual nucleo-base pair steps and by analyzing the sequence dependent estimates of intra-base pair entropy and inter-base pair entropy. The GGC triplet, which is part of the central tetramer (GGCC) of I, is concluded to be critical for binding of flavonoids, while the effect of the presence of ligand to the flexibility of nucleobases is localized through the intra-base pair motion of the intercalation site and its immediate vicinity. G-rich DNA sequences with consecutive Gs going before and/or after the critical GGC code (such as I: GGGGCCCC) are proposed to be uniquely specific for flavonoids. The configurational entropy contribution, as an upper bound of the true entropy contribution to the free energy in noncovalent binding, is demonstrated to influence the fundamental discrimination (intercalation vs groove binding) of DNA-flavonoid recognition modes. Some interesting implications for the structure-based design of optimal DNA binders are discussed.

Role of microscopic flexibility in tightly curved DNA

The genetic material in living cells is organized into complex structures in which DNA is subjected to substantial contortions. Here we investigate the difference in structure, dynamics, and flexibility between two topological states of a short (107 base pair) DNA sequence in a linear form and a covalently closed, tightly curved circular DNA form. By employing a combination of all-atom molecular dynamics (MD) simulations and elastic rod modeling of DNA, which allows capturing microscopic details while monitoring the global dynamics, we demonstrate that in the highly curved regime the microscopic flexibility of the DNA drastically increases due to the local mobility of the duplex. By analyzing vibrational entropy and Lipari–Szabo NMR order parameters from the simulation data, we propose a novel model for the thermodynamic stability of high-curvature DNA states based on vibrational untightening of the duplex. This novel view of DNA bending provides a fundamental explanation that bridges the gap between classical models of DNA and experimental studies on DNA cyclization, which so far have been in substantial disagreement.

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.

The allosteric modulation of thyroxine-binding globulin affinity is entropy driven

BACKGROUND

Thyroxine-binding globulin (TBG) is a non-inhibitory member of the serpin family of proteins whose main structural element is the reactive center loop (RCL), that, upon cleavage by proteases, is inserted into the protein core adopting a β-strand conformation (stressed to relaxed transition, S-to-R). After S-to-R transition thyroxine (T4) affinity decreases. However, crystallographic studies in the presence or absence of the hormone in different states are unable to show significant differences in the structure and interactions of the binding site. Experimental results also suggest the existence of several S states (differing in the number of inserted RCL residues), associated with a differential affinity.

METHODS

To shed light into the molecular basis that regulates T4 affinity according to the degree of RCL insertion in TBG, we performed extended molecular dynamics simulations combined with several thermodynamic analysis of the T4 binding to TBG in three different S states, and in the R state.

RESULTS

Our results show that, despite T4 binding in the protein by similar interactions in all states, a good correlation between the degree of RCL insertion and the binding affinity, driven by a change in TBG conformational entropy, was observed.

CONCLUSION

TBG allosteric regulation is entropy driven. The presence of multiple S states may allow more efficient T4 release due to protease activity.

GENERAL SIGNIFICANCE

The presented results are clear examples of how computer simulation methods can reveal the thermodynamic basis of allosteric effects, and provide a general framework for understanding serpin allosteric affinity regulation.

Drug-DNA intercalation: from discovery to the molecular mechanism

The ability of small molecules to perturb the natural structure and dynamics of nucleic acids is intriguing and has potential applications in cancer therapeutics. Intercalation is a special binding mode where the planar aromatic moiety of a small molecule is inserted between a pair of base pairs, causing structural changes in the DNA and leading to its functional arrest. Enormous progress has been made to understand the nature of the intercalation process since its idealistic conception five decades ago. However, the biological functions were detected even earlier. In this review, we focus mainly on the acridine and anthracycline types of drugs and provide a brief overview of the development in the field through various experimental methods that led to our present understanding of the subject. Subsequently, we discuss the molecular mechanism of the intercalation process, free-energy landscapes, and kinetics that was revealed recently through detailed and rigorous computational studies.

Relation between the change in DNA elasticity on ligand binding and the binding energetics

The widespread use of tweezers for measurement of ligand-DNA binding parameters is based on the McGhee-von Hippel treatment of the DNA contour and persistence length as a function of concentration. The McGhee-von Hippel approach contains the basic assumption that the binding constant K is independent of the number of already bound ligands. However, the change in elasticity of DNA on binding affects the entropic part of the Gibbs free energy and, hence, the K value in a concentration-dependent manner, making the whole approach inconsistent. In the present work we show that the energetic effect of DNA stiffening on noncovalent binding of small ligands is negligible with respect to the net energy of reaction, whereas the DNA stiffening on binding of large ligands must always be considered in each particular case.

Applications of isothermal titration calorimetry in pure and applied research--survey of the literature from 2010

Isothermal titration calorimetry (ITC) is a biophysical technique for measuring the formation and dissociation of molecular complexes and has become an invaluable tool in many branches of science from cell biology to food chemistry. By measuring the heat absorbed or released during bond formation, ITC provides accurate, rapid, and label-free measurement of the thermodynamics of molecular interactions. In this review, we survey the recent literature reporting the use of ITC and have highlighted a number of interesting studies that provide a flavour of the diverse systems to which ITC can be applied. These include measurements of protein-protein and protein-membrane interactions required for macromolecular assembly, analysis of enzyme kinetics, experimental validation of molecular dynamics simulations, and even in manufacturing applications such as food science. Some highlights include studies of the biological complex formed by Staphylococcus aureus enterotoxin C3 and the murine T-cell receptor, the mechanism of membrane association of the Parkinson's disease-associated protein α-synuclein, and the role of non-specific tannin-protein interactions in the quality of different beverages. Recent developments in automation are overcoming limitations on throughput imposed by previous manual procedures and promise to greatly extend usefulness of ITC in the future. We also attempt to impart some practical advice for getting the most out of ITC data for those researchers less familiar with the method.

Large-scale symmetry-adapted perturbation theory computations via density fitting and Laplace transformation techniques: investigating the fundamental forces of DNA-intercalator interactions

Symmetry-adapted perturbation theory (SAPT) provides a means of probing the fundamental nature of intermolecular interactions. Low-orders of SAPT (here, SAPT0) are especially attractive since they provide qualitative (sometimes quantitative) results while remaining tractable for large systems. The application of density fitting and Laplace transformation techniques to SAPT0 can significantly reduce the expense associated with these computations and make even larger systems accessible. We present new factorizations of the SAPT0 equations with density-fitted two-electron integrals and the first application of Laplace transformations of energy denominators to SAPT. The improved scalability of the DF-SAPT0 implementation allows it to be applied to systems with more than 200 atoms and 2800 basis functions. The Laplace-transformed energy denominators are compared to analogous partial Cholesky decompositions of the energy denominator tensor. Application of our new DF-SAPT0 program to the intercalation of DNA by proflavine has allowed us to determine the nature of the proflavine-DNA interaction. Overall, the proflavine-DNA interaction contains important contributions from both electrostatics and dispersion. The energetics of the intercalator interaction are are dominated by the stacking interactions (two-thirds of the total), but contain important contributions from the intercalator-backbone interactions. It is hypothesized that the geometry of the complex will be determined by the interactions of the intercalator with the backbone, because by shifting toward one side of the backbone, the intercalator can form two long hydrogen-bonding type interactions. The long-range interactions between the intercalator and the next-nearest base pairs appear to be negligible, justifying the use of truncated DNA models in computational studies of intercalation interaction energies.