1
|
Jia Z, Li L, Chakravorty A, Alexov E. Treating ion distribution with Gaussian-based smooth dielectric function in DelPhi. J Comput Chem 2017; 38:1974-1979. [PMID: 28602026 PMCID: PMC5495612 DOI: 10.1002/jcc.24831] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/14/2017] [Revised: 03/21/2017] [Accepted: 04/22/2017] [Indexed: 11/06/2022]
Abstract
The standard treatment of ions in the framework of the Poisson-Boltzmann equation relies on molecular surfaces, which are commonly constructed along with the Stern layer. The molecular surface determines where ions can be present. In the Gaussian-based smooth dielectric function in DelPhi, smooth boundaries between the solute and solvent take the place of molecular surface. Therefore, this invokes the question of how to model mobile ions in the water phase without a definite solute-solvent boundary. This article reports a natural extension of the Gaussian-based smooth dielectric function approach that treats mobile ions via Boltzmann distribution with an added desolvation penalty. Thus, ion concentration near macromolecules is governed by the local electrostatic potential and the desolvation penalty (from being partially desolvated). The approach is tested against the experimental salt dependence of binding free energy on 7 protein-protein complexes and 12 DNA-protein complexes, resulting in Pearson correlations of 0.95 and 0.88, respectively. © 2017 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Zhe Jia
- Computational Biophysics and Bioinformatics, Department of Physics and Astronomy, Clemson University, Clemson, South Carolina, United States, 29634
| | - Lin Li
- Computational Biophysics and Bioinformatics, Department of Physics and Astronomy, Clemson University, Clemson, South Carolina, United States, 29634
| | - Arghya Chakravorty
- Computational Biophysics and Bioinformatics, Department of Physics and Astronomy, Clemson University, Clemson, South Carolina, United States, 29634
| | - Emil Alexov
- Computational Biophysics and Bioinformatics, Department of Physics and Astronomy, Clemson University, Clemson, South Carolina, United States, 29634
| |
Collapse
|
2
|
Dragan AI, Read CM, Crane-Robinson C. Enthalpy-entropy compensation: the role of solvation. EUROPEAN BIOPHYSICS JOURNAL: EBJ 2016; 46:301-308. [PMID: 27796417 PMCID: PMC5384952 DOI: 10.1007/s00249-016-1182-6] [Citation(s) in RCA: 72] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 08/22/2016] [Revised: 10/06/2016] [Accepted: 10/17/2016] [Indexed: 11/24/2022]
Abstract
Structural modifications to interacting systems frequently lead to changes in both the enthalpy (heat) and entropy of the process that compensate each other, so that the Gibbs free energy is little changed: a major barrier to the development of lead compounds in drug discovery. The conventional explanation for such enthalpy–entropy compensation (EEC) is that tighter contacts lead to a more negative enthalpy but increased molecular constraints, i.e., a compensating conformational entropy reduction. Changes in solvation can also contribute to EEC but this contribution is infrequently discussed. We review long-established and recent cases of EEC and conclude that the large fluctuations in enthalpy and entropy observed are too great to be a result of only conformational changes and must result, to a considerable degree, from variations in the amounts of water immobilized or released on forming complexes. Two systems exhibiting EEC show a correlation between calorimetric entropies and local mobilities, interpreted to mean conformational control of the binding entropy/free energy. However, a substantial contribution from solvation gives the same effect, as a consequence of a structural link between the amount of bound water and the protein flexibility. Only by assuming substantial changes in solvation—an intrinsically compensatory process—can a more complete understanding of EEC be obtained. Faced with such large, and compensating, changes in the enthalpies and entropies of binding, the best approach to engineering elevated affinities must be through the addition of ionic links, as they generate increased entropy without affecting the enthalpy.
Collapse
Affiliation(s)
- Anatoliy I Dragan
- Institute of High Technologies, Taras Shevchenko National University of Kyiv, 64, Volodymyrs'ka St., Kiev, 01601, Ukraine.,Institute of Molecular Biology and Genetics, NASU, 150, Zabolotnogo St., Kiev, 03680, Ukraine
| | - Christopher M Read
- Biophysics Laboratories, School of Biology, University of Portsmouth, Portsmouth, PO1 2DT, UK
| | - Colyn Crane-Robinson
- Biophysics Laboratories, School of Biology, University of Portsmouth, Portsmouth, PO1 2DT, UK.
| |
Collapse
|
3
|
Shukla H, Vaitiekunas P, Majumdar AK, Dragan AI, Dimitriadis EK, Kotova S, Crane-Robinson C, Privalov PL. The linker of the interferon response factor 3 transcription factor is not unfolded. Biochemistry 2012; 51:6320-7. [PMID: 22812703 DOI: 10.1021/bi300260s] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Interferon response factor 3 (IRF-3) is a transcription factor that plays an essential role in controlling the synthesis of interferon-β (IFN-β) and is a protein consisting of two well-defined domains, the N-terminal DNA-binding and the C-terminal dimerization domains, connected by a 75-residue linker, supposedly unfolded. However, it was not clear whether in intact IRF-3 this linker segment of the chain, which carries the nuclear export signal and includes a region of high helical propensity, remains unfolded. This has been investigated using nuclear magnetic resonance by ligating the (15)N-labeled linker to the unlabeled N-terminal and C-terminal domains. It was found that, while the linker alone is indeed in a completely unfolded state, when ligated to the C-terminal domain it shows some ordering, and this ordering becomes much more pronounced when the linker is also ligated to the N-terminal domain. Thus, in intact IRF-3, the linker represents a folded structural domain; i.e., IRF-3 is a three-domain globular protein. Light scattering studies of wild-type IRF-3 showed that these three domains are tightly packed, and therefore, the dimer of IRF-3, which is formed upon phosphorylation of its C-terminal domains following virus invasion, must be a rather rigid and compact construction. One would then expect that binding of such a dimer to its tandem recognition sites PRDIII and PRDI, which are located on opposing faces of the IFN-β enhancer DNA, should result in deformation of the DNA. Analysis of the characteristics of binding of the monomeric and dimeric IRF-3 to the enhancer DNA indeed showed that formation of this complex requires considerable work for deformation of its components, most likely bending of the DNA. Such bending was confirmed by atomic force microscopy of dimeric IRF-3 bound to the PRDII-PRDI tandem recognition sites placed at the middle of a 300 bp DNA probe. Bending of DNA by IRF-3 must be significant in the assembly and function of the IFN-β enhancer.
Collapse
Affiliation(s)
- Hem Shukla
- Department of Biology, Johns Hopkins University, Baltimore, MD 21218, USA
| | | | | | | | | | | | | | | |
Collapse
|
4
|
Abstract
In the present paper, we outline how to construct the partition function for a protein using empirical heat capacity data. The procedure is based on the calculation of a set of energy moments from the temperature dependence of the heat capacity. Given a set of energy moments, one can then use the maximum-entropy method to calculate an approximate energy distribution for the protein; the more energy moments one has, the better the approximation. The energy distribution can then be used to calculate the probability that the molecule is in a given energy level, which, using standard statistical mechanics, gives the degeneracy of the particular energy level. The degeneracy as a function of energy is the central ingredient in the construction of the partition function. Given the partition function, one can calculate all of the thermodynamic functions of the protein (free energy, energy, entropy, heat capacity, and energy probability distribution) as a function of temperature. The three-dimensional plot of the probability that the protein has a given energy at a given temperature tells one graphically (without imposing the assumption) whether or not it is a good approximation to divide the terms in the partition function into two or more groups, reflecting, for example, the presence of distinct native and denatured species.
Collapse
Affiliation(s)
- Douglas Poland
- Department of Chemistry, The Johns Hopkins University, Baltimore, Maryland 21218, USA
| |
Collapse
|
5
|
De Ioannes P, Escalante CR, Aggarwal AK. Structures of apo IRF-3 and IRF-7 DNA binding domains: effect of loop L1 on DNA binding. Nucleic Acids Res 2011; 39:7300-7. [PMID: 21596780 PMCID: PMC3167601 DOI: 10.1093/nar/gkr325] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/03/2022] Open
Abstract
Interferon regulatory factors IRF-3 and IRF-7 are transcription factors essential in the activation of interferon-β (IFN-β) gene in response to viral infections. Although, both proteins recognize the same consensus IRF binding site AANNGAAA, they have distinct DNA binding preferences for sites in vivo. The X-ray structures of IRF-3 and IRF-7 DNA binding domains (DBDs) bound to IFN-β promoter elements revealed flexibility in the loops (L1–L3) and the residues that make contacts with the target sequence. To characterize the conformational changes that occur on DNA binding and how they differ between IRF family members, we have solved the X-ray structures of IRF-3 and IRF-7 DBDs in the absence of DNA. We found that loop L1, carrying the conserved histidine that interacts with the DNA minor groove, is disordered in apo IRF-3 but is ordered in apo IRF-7. This is reflected in differences in DNA binding affinities when the conserved histidine in loop L1 is mutated to alanine in the two proteins. The stability of loop L1 in IRF-7 derives from a unique combination of hydrophobic residues that pack against the protein core. Together, our data show that differences in flexibility of loop L1 are an important determinant of differential IRF-DNA binding.
Collapse
Affiliation(s)
- Pablo De Ioannes
- Department of Structural and Chemical Biology, Mount Sinai School of Medicine, Box 1677, 1425 Madison Avenue, NY 10029, USA
| | | | | |
Collapse
|
6
|
Privalov PL, Dragan AI, Crane-Robinson C. The cost of DNA bending. Trends Biochem Sci 2009; 34:464-70. [PMID: 19726198 DOI: 10.1016/j.tibs.2009.05.005] [Citation(s) in RCA: 35] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/12/2009] [Revised: 05/13/2009] [Accepted: 05/14/2009] [Indexed: 12/22/2022]
Abstract
Experimental data on protein-DNA interactions highlight a surprising peculiarity of protein binding to the minor groove: in contrast to major groove binding, which proceeds with heat release and does not induce substantial deformation of DNA, minor groove binding takes place at AT-rich sites, proceeds with heat absorption and results in significant DNA bending. By forming a highly ordered and dense spine in the minor groove of AT-rich DNA, water plays an essential role in defining the energetic signature of protein-minor groove binding. Removal of this water requires minimal work and results in significant loss of rigidity in the DNA, which can then easily acquire the conformation imposed by the bound protein. Therefore the introduction of substantial bends into the DNA is not energetically expensive.
Collapse
Affiliation(s)
- Peter L Privalov
- Department of Biology, Johns Hopkins University, Baltimore, MD 21218, USA.
| | | | | |
Collapse
|
7
|
Favicchio R, Dragan AI, Kneale GG, Read CM. Fluorescence spectroscopy and anisotropy in the analysis of DNA-protein interactions. Methods Mol Biol 2009; 543:589-611. [PMID: 19378154 DOI: 10.1007/978-1-60327-015-1_35] [Citation(s) in RCA: 26] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/07/2023]
Abstract
Fluorescence spectroscopy can be used as a sensitive non-destructive technique for the characterisation of protein-DNA interactions. A comparison of the intrinsic emission spectra obtained for a protein-DNA complex and for free protein can be informative about the environment of tryptophan and tyrosine residues in the two states. Often there is quenching of the fluorescence intensity of an intrinsic emission spectrum and/or a shift in the wavelength maximum on protein binding to DNA. A step-by-step protocol describes the determination of a DNA-binding curve by measurement of the quenching of the intrinsic protein fluorescence.Fluorescence anisotropy can also be used to obtain a DNA-binding curve if the molecular size of the protein-DNA complex is sufficiently different from the free fluorescing component. Typically an extrinsic fluorophore attached to one or both 5' ends of single-stranded or duplex DNA is used, for this increases the sensitivity of measurement.Fitting of the binding curves, assuming a model, can often yield the stoichiometry and association constant of the interaction. The approach is illustrated using the interaction of the DNA-binding domains (HMG boxes) of mouse Sox-5 and mammalian HMGB1 with short DNA duplexes.
Collapse
Affiliation(s)
- Rosy Favicchio
- Institute of Molecular Biology and Biotechnology, Foundation for Research and Technology Hellas, Heraklion, Crete, Greece
| | | | | | | |
Collapse
|
8
|
Dragan A, Carrillo R, Gerasimova T, Privalov P. Assembling the Human IFN-β Enhanceosome in Solution. J Mol Biol 2008; 384:335-48. [DOI: 10.1016/j.jmb.2008.09.015] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/13/2008] [Revised: 07/28/2008] [Accepted: 09/02/2008] [Indexed: 10/21/2022]
|
9
|
Hargreaves VV, Schleif RF. The salt dependence of the interferon regulatory factor 1 DNA binding domain binding to DNA reveals ions are localized around protein and DNA. Biochemistry 2008; 47:4119-28. [PMID: 18324782 DOI: 10.1021/bi702082q] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The equilibrium dissociation constant of the DNA binding domain of interferon regulatory factor 1 (IRF1 DBD) for its DNA binding site depends strongly on salt concentration and salt type. These dependencies are consistent with IRF1 DBD binding to DNA, resulting in the release of cations from the DNA and both release of anions from the protein and uptake of a cation by the protein. We demonstrated this by utilizing the fact that the release of fluoride from protein upon complex formation does not contribute to the salt concentration dependence of binding and by studying mutants in which charged residues in IRF1 DBD that form salt bridges with DNA phosphates are changed to alanine. The salt concentration dependencies of the dissociation constants of wild-type IRF1 DBD and the mutants R64A, D73A, K75A, and D73A/K75A were measured in buffer containing NaF, NaCl, or NaBr. The salt concentration and type dependencies of the mutants relative to wild-type IRF1 DBD provide evidence of charge neutralization by solution ions for R64 and by a salt bridge between D73 and K75 in buffer containing chloride or bromide salts. These data also allowed us to determine the number, type, and localization of condensed ions around both IRF1 DBD and its DNA binding site.
Collapse
|
10
|
Dragan AI, Hargreaves VV, Makeyeva EN, Privalov PL. Mechanisms of activation of interferon regulator factor 3: the role of C-terminal domain phosphorylation in IRF-3 dimerization and DNA binding. Nucleic Acids Res 2007; 35:3525-34. [PMID: 17483521 PMCID: PMC1920236 DOI: 10.1093/nar/gkm142] [Citation(s) in RCA: 41] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The interferon regulatory transcription factor (IRF-3) is activated by phosphorylation of Ser/Thr residues clustered in its C-terminal domain. Phosphorylation of these residues, which increases the negative charge of IRF-3, results in its dimerization and association with DNA, despite the increase in repulsive electrostatic interactions. To investigate this surprising effect, the dimerization of IRF-3 and two phosphomimetic mutants, 2D (S396D, S398D) and 5D (S396D, S398D, S402D, T404D and S405D), and their binding to single-site PRDI and double-site PRDIII-PRDI DNA sequences from the IFN-beta enhancer have been studied. It was found that: (a) the mutations in the C-terminal domain do not affect the state of the DNA-binding N-terminal domain or its ability to bind target DNA; (b) in the 5D-mutant, the local increase of negative charge in the C-terminal domain induces restructuring, resulting in the formation of a stable dimer; (c) dimerization of IRF-3 is the basis of its strong binding to PRDIII-PRDI sites since binding of 5D to the single PRDI site is similar to that of inactivated IRF-3. Analysis of the binding characteristics leads to the conclusion that binding of dimeric IRF-3 to the DNA with two tandem-binding sites, which are twisted by approximately 100 degrees relative to each other, requires considerable work to untwist and/or bend the DNA.
Collapse
Affiliation(s)
- Anatoly I Dragan
- Department of Biology, Johns Hopkins University, Baltimore, MD 21218, USA.
| | | | | | | |
Collapse
|
11
|
Privalov PL, Dragan AI. Microcalorimetry of biological macromolecules. Biophys Chem 2007; 126:16-24. [PMID: 16781052 DOI: 10.1016/j.bpc.2006.05.004] [Citation(s) in RCA: 105] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2006] [Accepted: 05/08/2006] [Indexed: 10/24/2022]
Abstract
The capabilities of contemporary differential scanning and isothermal titration microcalorimetry for studying the thermodynamics of protein unfolding/refolding and their association with partners, particularly target DNA duplexes, are considered. It is shown that the predenaturational changes of proteins must not be ignored in studying the thermodynamics of formation of their native structure and their complexes with partners, particularly their cognate DNA duplexes.
Collapse
|
12
|
Privalov PL, Dragan AI, Crane-Robinson C, Breslauer KJ, Remeta DP, Minetti CASA. What drives proteins into the major or minor grooves of DNA? J Mol Biol 2006; 365:1-9. [PMID: 17055530 PMCID: PMC1934558 DOI: 10.1016/j.jmb.2006.09.059] [Citation(s) in RCA: 155] [Impact Index Per Article: 8.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2006] [Revised: 09/21/2006] [Accepted: 09/22/2006] [Indexed: 01/22/2023]
Abstract
The energetic profiles of a significant number of protein-DNA systems at 20 degrees C reveal that, despite comparable Gibbs free energies, association with the major groove is primarily an enthalpy-driven process, whereas binding to the minor groove is characterized by an unfavorable enthalpy that is compensated by favorable entropic contributions. These distinct energetic signatures for major versus minor groove binding are irrespective of the magnitude of DNA bending and/or the extent of binding-induced protein refolding. The primary determinants of their different energetic profiles appear to be the distinct hydration properties of the major and minor grooves; namely, that the water in the A+T-rich minor groove is in a highly ordered state and its removal results in a substantial positive contribution to the binding entropy. Since the entropic forces driving protein binding into the minor groove are a consequence of displacing water ordered by the regular arrangement of polar contacts, they cannot be regarded as hydrophobic.
Collapse
Affiliation(s)
- Peter L Privalov
- Department of Biology, Johns Hopkins University, Baltimore, MD 21218, USA.
| | | | | | | | | | | |
Collapse
|