1
|
Kariyawasam NL, Ploetz EA, Swint-Kruse L, Smith PE. Simulated pressure changes in LacI suggest a link between hydration and functional conformational changes. Biophys Chem 2024; 304:107126. [PMID: 37924711 PMCID: PMC10842697 DOI: 10.1016/j.bpc.2023.107126] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/28/2023] [Revised: 10/18/2023] [Accepted: 10/25/2023] [Indexed: 11/06/2023]
Abstract
The functions of many proteins are associated with interconversions among conformational substates. However, these substates can be difficult to measure experimentally, and determining contributions from hydration changes can be especially difficult. Here, we assessed the use of pressure perturbations to sample the substates accessible to the Escherichia coli lactose repressor protein (LacI) in various liganded forms. In the presence of DNA, the regulatory domain of LacI adopts an Open conformation that, in the absence of DNA, changes to a Closed conformation. Increasing the simulation pressure prevented the transition from an Open to a Closed conformation, in a similar manner to the binding of DNA and anti-inducer, ONPF. The results suggest the hydration of specific residues play a significant role in determining the population of different LacI substates and that simulating pressure perturbation could be useful for assessing the role of hydration changes that accompany functionally-relevant amino acid substitutions.
Collapse
Affiliation(s)
- Nilusha L Kariyawasam
- Department of Chemistry, 213 CBC Building, 1212 Mid-Campus Dr. North, Kansas State University, Manhattan, KS 66506, USA
| | - Elizabeth A Ploetz
- Department of Chemistry, 213 CBC Building, 1212 Mid-Campus Dr. North, Kansas State University, Manhattan, KS 66506, USA
| | - Liskin Swint-Kruse
- Department of Biochemistry and Molecular Biology, MSN 3030, 3901 Rainbow Blvd, The University of Kansas Medical Center, Kansas City, KS 66160, USA
| | - Paul E Smith
- Department of Chemistry, 213 CBC Building, 1212 Mid-Campus Dr. North, Kansas State University, Manhattan, KS 66506, USA.
| |
Collapse
|
2
|
Lüking M, Elf J, Levy Y. Conformational Change of Transcription Factors from Search to Specific Binding: A lac Repressor Case Study. J Phys Chem B 2022; 126:9971-9984. [PMID: 36416228 PMCID: PMC9743208 DOI: 10.1021/acs.jpcb.2c05006] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
Abstract
In a process known as facilitated diffusion, DNA-binding proteins find their target sites by combining three-dimensional diffusion and one-dimensional scanning of the DNA. Following the trade-off between speed and stability, agile exploration of DNA requires loose binding, whereas, at the DNA target site, the searching protein needs to establish tight interactions with the DNA. To enable both efficient search and stable binding, DNA-binding proteins and DNA often switch conformations upon recognition. Here, we study the one-dimensional diffusion and DNA binding of the dimeric lac repressor (LacI), which was reported to adopt two different conformations when binding different conformations of DNA. Using coarse-grained molecular dynamic simulations, we studied the diffusion and the sequence-specific binding of these conformations of LacI, as well as their truncated or monomeric variants, with two DNA conformations: straight and bent. The simulations were compared to experimental observables. This study supports that linear diffusion along DNA combines tight rotation-coupled groove tracking and rotation-decoupled hopping, where the protein briefly dissociates and reassociates just a few base pairs away. Tight groove tracking is crucial for target-site recognition, while hopping speeds up the overall search process. We investigated the diffusion of different LacI conformations on DNA and show how the flexibility of LacI's hinge regions ensures agility on DNA as well as faithful groove tracking. If the hinge regions instead form α-helices at the protein-DNA interface, tight groove tracking is not possible. On the contrary, the helical hinge region is essential for tight binding to bent, specific DNA, for the formation of the specific complex. Based on our study of different encounter complexes, we argue that the conformational change in LacI and DNA bending are somewhat coupled. Our findings underline the importance of two distinct protein conformations for facilitated diffusion and specific binding, respectively.
Collapse
Affiliation(s)
- Malin Lüking
- Department
of Cell- and Molecular Biology-ICM, Uppsala
University, Uppsala, Uppsala County751 24, Sweden
| | - Johan Elf
- Department
of Cell- and Molecular Biology-ICM, Uppsala
University, Uppsala, Uppsala County751 24, Sweden
| | - Yaakov Levy
- Department
of Chemical and Structural Biology, Weizmann
Institute of Science, Rehovot, Central District76100, Israel,. Tel.: 972-8-9344587
| |
Collapse
|
3
|
Revealing atomic-scale molecular diffusion of a plant-transcription factor WRKY domain protein along DNA. Proc Natl Acad Sci U S A 2021; 118:2102621118. [PMID: 34074787 PMCID: PMC8201915 DOI: 10.1073/pnas.2102621118] [Citation(s) in RCA: 17] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022] Open
Abstract
In transcription factors’ search for target genes, one-dimensional diffusion of the protein along DNA is essential. Experimentally, it remains challenging to resolve the individual diffusional steps of protein on DNA. Here, we report mainly all-atom equilibrium simulations of a WRKY domain protein in association with and diffusion along DNA. We demonstrate a complete stepping cycle of the protein for one base pair on DNA within microseconds, along with stochastic motions. Processive protein diffusions on DNA have been further sampled in a coarse-grained model. We have also found preferential DNA-strand association of the domain protein, which becomes most prominent at specific DNA binding, and it can be common for small-domain proteins to balance movements on the DNA with the sequence recognition. Transcription factor (TF) target search on genome is highly essential for gene expression and regulation. High-resolution determination of TF diffusion along DNA remains technically challenging. Here, we constructed a TF model system using the plant WRKY domain protein in complex with DNA from crystallography and demonstrated microsecond diffusion dynamics of WRKY on DNA by employing all-atom molecular-dynamics (MD) simulations. Notably, we found that WRKY preferentially binds to one strand of DNA with significant energetic bias compared with the other, or nonpreferred strand. The preferential DNA-strand binding becomes most prominent in the static process, from nonspecific to specific DNA binding, but less distinct during diffusive movements of the domain protein on the DNA. Remarkably, without employing acceleration forces or bias, we captured a complete one-base-pair stepping cycle of the protein tracking along major groove of DNA with a homogeneous poly-adenosine sequence, as individual hydrogen bonds break and reform at the protein–DNA binding interface. Further DNA-groove tracking motions of the protein forward or backward, with occasional sliding as well as strand crossing to minor groove of DNA, were also captured. The processive diffusion of WRKY along DNA has been further sampled via coarse-grained MD simulations. The study thus provides structural dynamics details on diffusion of a small TF domain protein, suggests how the protein approaches a specific recognition site on DNA, and supports further high-precision experimental detection. The stochastic movements revealed in the TF diffusion also provide general clues about how other protein walkers step and slide along DNA.
Collapse
|
4
|
Liao Q, Lüking M, Krüger DM, Deindl S, Elf J, Kasson PM, Lynn Kamerlin SC. Long Time-Scale Atomistic Simulations of the Structure and Dynamics of Transcription Factor-DNA Recognition. J Phys Chem B 2019; 123:3576-3590. [PMID: 30952192 DOI: 10.1021/acs.jpcb.8b12363] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
Recent years have witnessed an explosion of interest in computational studies of DNA binding proteins, including both coarse-grained and atomistic simulations of transcription factor-DNA recognition, to understand how these transcription factors recognize their binding sites on the DNA with such exquisite specificity. The present study performs microsecond time scale all-atom simulations of the dimeric form of the lactose repressor (LacI), both in the absence of any DNA and in the presence of both specific and nonspecific complexes, considering three different DNA sequences. We examine, specifically, the conformational differences between specific and nonspecific protein-DNA interactions, as well as the behavior of the helix-turn-helix motif of LacI when interacting with the DNA. Our simulations suggest that stable LacI binding occurs primarily to bent A-form DNA, with a loss of LacI conformational entropy and optimization of correlated conformational equilibria across the protein. In addition, binding to the specific operator sequence involves a slightly larger number of stabilizing DNA-protein hydrogen bonds (in comparison to nonspecific complexes), which may account for the experimentally observed specificity for this operator. In doing so, our simulations provide a detailed atomistic description of potential structural drivers for LacI selectivity.
Collapse
Affiliation(s)
- Qinghua Liao
- Science for Life Laboratory, Department of Chemistry-BMC , Uppsala University , BMC Box 576, S-751 24 Uppsala , Sweden
| | - Malin Lüking
- Science for Life Laboratory, Department of Chemistry-BMC , Uppsala University , BMC Box 576, S-751 24 Uppsala , Sweden
| | - Dennis M Krüger
- Science for Life Laboratory, Department of Cell and Molecular Biology , Uppsala University , BMC Box 596, S-751 23 Uppsala , Sweden.,Department for Epigenetics and Systems Medicine in Neurodegenerative Diseases, Bioinformatics Unit , German Center for Neurodegenerative Diseases, Göttingen , von Siebold Strasse 3A , 37075 Göttingen , Germany
| | - Sebastian Deindl
- Science for Life Laboratory, Department of Cell and Molecular Biology , Uppsala University , BMC Box 596, S-751 23 Uppsala , Sweden
| | - Johan Elf
- Science for Life Laboratory, Department of Cell and Molecular Biology , Uppsala University , BMC Box 596, S-751 23 Uppsala , Sweden
| | - Peter M Kasson
- Science for Life Laboratory, Department of Cell and Molecular Biology , Uppsala University , BMC Box 596, S-751 23 Uppsala , Sweden
| | - Shina Caroline Lynn Kamerlin
- Science for Life Laboratory, Department of Chemistry-BMC , Uppsala University , BMC Box 576, S-751 24 Uppsala , Sweden
| |
Collapse
|
5
|
Seckfort D, Montgomery Pettitt B. Price of disorder in the lac repressor hinge helix. Biopolymers 2019; 110:e23239. [PMID: 30485404 PMCID: PMC6335174 DOI: 10.1002/bip.23239] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/11/2018] [Revised: 09/12/2018] [Accepted: 10/04/2018] [Indexed: 12/26/2022]
Abstract
The Lac system of genes has been pivotal in understanding gene regulation. When the lac repressor protein binds to the correct DNA sequence, the hinge region of the protein goes through a disorder to order transition. The structure of this region of the protein is well understood when it is in this bound conformation, but less so when it is not. Structural studies show that this region is flexible. Our simulations show this region is extremely flexible in solution; however, a high concentration of salt can help kinetically trap the hinge helix. Thermodynamically, disorder is more favorable without the DNA present.
Collapse
Affiliation(s)
- Danielle Seckfort
- Structural and Computational Biology and Molecular Biophysics, Baylor College of Medicine, Houston, Texas
| | - B Montgomery Pettitt
- Structural and Computational Biology and Molecular Biophysics, Baylor College of Medicine, Houston, Texas
- Department of Biochemistry and Molecular Biology, Sealy Center for Structural Biology, University of Texas Medical Branch, Galveston, Texas
| |
Collapse
|
6
|
Wieczór M, Czub J. How proteins bind to DNA: target discrimination and dynamic sequence search by the telomeric protein TRF1. Nucleic Acids Res 2017. [PMID: 28633355 PMCID: PMC5737604 DOI: 10.1093/nar/gkx534] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023] Open
Abstract
Target search as performed by DNA-binding proteins is a complex process, in which multiple factors contribute to both thermodynamic discrimination of the target sequence from overwhelmingly abundant off-target sites and kinetic acceleration of dynamic sequence interrogation. TRF1, the protein that binds to telomeric tandem repeats, faces an intriguing variant of the search problem where target sites are clustered within short fragments of chromosomal DNA. In this study, we use extensive (>0.5 ms in total) MD simulations to study the dynamical aspects of sequence-specific binding of TRF1 at both telomeric and non-cognate DNA. For the first time, we describe the spontaneous formation of a sequence-specific native protein-DNA complex in atomistic detail, and study the mechanism by which proteins avoid off-target binding while retaining high affinity for target sites. Our calculated free energy landscapes reproduce the thermodynamics of sequence-specific binding, while statistical approaches allow for a comprehensive description of intermediate stages of complex formation.
Collapse
Affiliation(s)
- Milosz Wieczór
- Department of Physical Chemistry, Gdansk University of Technology, ul. Narutowicza 11/12, 80-233 Gdansk, Poland
| | - Jacek Czub
- Department of Physical Chemistry, Gdansk University of Technology, ul. Narutowicza 11/12, 80-233 Gdansk, Poland
| |
Collapse
|
7
|
Sun L, Tabaka M, Hou S, Li L, Burdzy K, Aksimentiev A, Maffeo C, Zhang X, Holyst R. The Hinge Region Strengthens the Nonspecific Interaction between Lac-Repressor and DNA: A Computer Simulation Study. PLoS One 2016; 11:e0152002. [PMID: 27008630 PMCID: PMC4805274 DOI: 10.1371/journal.pone.0152002] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/13/2015] [Accepted: 03/06/2016] [Indexed: 11/30/2022] Open
Abstract
LacI is commonly used as a model to study the protein-DNA interaction and gene regulation. The headpiece of the lac-repressor (LacI) protein is an ideal system for investigation of nonspecific binding of the whole LacI protein to DNA. The hinge region of the headpiece has been known to play a key role in the specific binding of LacI to DNA, whereas its role in nonspecific binding process has not been elucidated. Here, we report the results of explicit solvent molecular dynamics simulation and continuum electrostatic calculations suggesting that the hinge region strengthens the nonspecific interaction, accounting for up to 50% of the micro-dissociation free energy of LacI from DNA. Consequently, the rate of microscopic dissociation of LacI from DNA is reduced by 2~3 orders of magnitude in the absence of the hinge region. We find the hinge region makes an important contribution to the electrostatic energy, the salt dependence of electrostatic energy, and the number of salt ions excluded from binding of the LacI-DNA complex.
Collapse
Affiliation(s)
- Lili Sun
- Institute of Physical Chemistry PAS, Kasprzaka 44/52, 01–224, Warsaw, Poland
| | - Marcin Tabaka
- Institute of Physical Chemistry PAS, Kasprzaka 44/52, 01–224, Warsaw, Poland
| | - Sen Hou
- Institute of Physical Chemistry PAS, Kasprzaka 44/52, 01–224, Warsaw, Poland
| | - Lin Li
- Computational Biophysics and Bioinformatics, Department of Physics, Clemson University, Clemson, South Carolina, 29634, United States of America
| | - Krzysztof Burdzy
- Department of Mathematics, University of Washington, Seattle, Washington, 98195–4350, United States of America
| | - Aleksei Aksimentiev
- Department of Physics, University of Illinois, Urbana, Illinois, 61801, United States of America
| | - Christopher Maffeo
- Department of Physics, University of Illinois, Urbana, Illinois, 61801, United States of America
| | - Xuzhu Zhang
- Institute of Physical Chemistry PAS, Kasprzaka 44/52, 01–224, Warsaw, Poland
| | - Robert Holyst
- Institute of Physical Chemistry PAS, Kasprzaka 44/52, 01–224, Warsaw, Poland
- * E-mail:
| |
Collapse
|
8
|
Sun L, Tabaka M, Hou S, Li L, Burdzy K, Aksimentiev A, Maffeo C, Zhang X, Holyst R. The Hinge Region Strengthens the Nonspecific Interaction between Lac-Repressor and DNA: A Computer Simulation Study. PLoS One 2016; 11:e0152002. [PMID: 27008630 DOI: 10.1371/joumal.pone.0152002] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/13/2015] [Accepted: 03/06/2016] [Indexed: 05/27/2023] Open
Abstract
LacI is commonly used as a model to study the protein-DNA interaction and gene regulation. The headpiece of the lac-repressor (LacI) protein is an ideal system for investigation of nonspecific binding of the whole LacI protein to DNA. The hinge region of the headpiece has been known to play a key role in the specific binding of LacI to DNA, whereas its role in nonspecific binding process has not been elucidated. Here, we report the results of explicit solvent molecular dynamics simulation and continuum electrostatic calculations suggesting that the hinge region strengthens the nonspecific interaction, accounting for up to 50% of the micro-dissociation free energy of LacI from DNA. Consequently, the rate of microscopic dissociation of LacI from DNA is reduced by 2~3 orders of magnitude in the absence of the hinge region. We find the hinge region makes an important contribution to the electrostatic energy, the salt dependence of electrostatic energy, and the number of salt ions excluded from binding of the LacI-DNA complex.
Collapse
Affiliation(s)
- Lili Sun
- Institute of Physical Chemistry PAS, Kasprzaka 44/52, 01-224, Warsaw, Poland
| | - Marcin Tabaka
- Institute of Physical Chemistry PAS, Kasprzaka 44/52, 01-224, Warsaw, Poland
| | - Sen Hou
- Institute of Physical Chemistry PAS, Kasprzaka 44/52, 01-224, Warsaw, Poland
| | - Lin Li
- Computational Biophysics and Bioinformatics, Department of Physics, Clemson University, Clemson, South Carolina, 29634, United States of America
| | - Krzysztof Burdzy
- Department of Mathematics, University of Washington, Seattle, Washington, 98195-4350, United States of America
| | - Aleksei Aksimentiev
- Department of Physics, University of Illinois, Urbana, Illinois, 61801, United States of America
| | - Christopher Maffeo
- Department of Physics, University of Illinois, Urbana, Illinois, 61801, United States of America
| | - Xuzhu Zhang
- Institute of Physical Chemistry PAS, Kasprzaka 44/52, 01-224, Warsaw, Poland
| | - Robert Holyst
- Institute of Physical Chemistry PAS, Kasprzaka 44/52, 01-224, Warsaw, Poland
| |
Collapse
|
9
|
Anwar MA, Yesudhas D, Shah M, Choi S. Structural and conformational insights into SOX2/OCT4-bound enhancer DNA: a computational perspective. RSC Adv 2016. [DOI: 10.1039/c6ra15176k] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
Abstract
The roles of SOX2 and OCT4 are critical in stem cell maintenance either in the context of iPSCs generation or cancer stem cell growth; therefore, it is imperative to study their cooperative binding and SOX2/OCT4-induced DNA conformational switching.
Collapse
Affiliation(s)
- Muhammad Ayaz Anwar
- Department of Molecular Science and Technology
- Ajou University
- Suwon 443-749
- Korea
| | - Dhanusha Yesudhas
- Department of Molecular Science and Technology
- Ajou University
- Suwon 443-749
- Korea
| | - Masaud Shah
- Department of Molecular Science and Technology
- Ajou University
- Suwon 443-749
- Korea
| | - Sangdun Choi
- Department of Molecular Science and Technology
- Ajou University
- Suwon 443-749
- Korea
| |
Collapse
|
10
|
Etheve L, Martin J, Lavery R. Dynamics and recognition within a protein-DNA complex: a molecular dynamics study of the SKN-1/DNA interaction. Nucleic Acids Res 2015; 44:1440-8. [PMID: 26721385 PMCID: PMC4756839 DOI: 10.1093/nar/gkv1511] [Citation(s) in RCA: 46] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/23/2015] [Accepted: 12/15/2015] [Indexed: 11/13/2022] Open
Abstract
Molecular dynamics simulations of the Caenorhabditis elegans transcription factor SKN-1 bound to its cognate DNA site show that the protein–DNA interface undergoes significant dynamics on the microsecond timescale. A detailed analysis of the simulation shows that movements of two key arginine side chains between the major groove and the backbone of DNA generate distinct conformational sub-states that each recognize only part of the consensus binding sequence of SKN-1, while the experimentally observed binding specificity results from a time-averaged view of the dynamic recognition occurring within this complex.
Collapse
Affiliation(s)
- Loïc Etheve
- BMSSI UMR 5086 CNRS/Univ. Lyon I, Institut de Biologie et Chimie des Protéines, 7 passage du Vercors, Lyon 69367, France
| | - Juliette Martin
- BMSSI UMR 5086 CNRS/Univ. Lyon I, Institut de Biologie et Chimie des Protéines, 7 passage du Vercors, Lyon 69367, France
| | - Richard Lavery
- BMSSI UMR 5086 CNRS/Univ. Lyon I, Institut de Biologie et Chimie des Protéines, 7 passage du Vercors, Lyon 69367, France
| |
Collapse
|
11
|
Tabaka M, Burdzy K, Hołyst R. Method for the analysis of contribution of sliding and hopping to a facilitated diffusion of DNA-binding protein: Application to in vivo data. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2015; 92:022721. [PMID: 26382446 DOI: 10.1103/physreve.92.022721] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/11/2015] [Indexed: 06/05/2023]
Abstract
DNA-binding protein searches for its target, a specific site on DNA, by means of diffusion. The search process consists of many recurrent steps of one-dimensional diffusion (sliding) along the DNA chain and three-dimensional diffusion (hopping) after dissociation of a protein from the DNA chain. Here we propose a computational method that allows extracting the contribution of sliding and hopping to the search process in vivo from the measurements of the kinetics of the target search by the lac repressor in Escherichia coli [P. Hammar et al., Science 336, 1595 (2012)]. The method combines lattice Monte Carlo simulations with the Brownian excursion theory and includes explicitly steric constraints for hopping due to the helical structure of DNA. The simulation results including all experimental data reveal that the in vivo target search is dominated by sliding. The short-range hopping to the same base pair interrupts one-dimensional sliding while long-range hopping does not contribute significantly to the kinetics of the search of the target in vivo.
Collapse
Affiliation(s)
- Marcin Tabaka
- Institute of Physical Chemistry, Polish Academy of Sciences, Kasprzaka 44/52, 01-224 Warsaw, Poland
| | - Krzysztof Burdzy
- Department of Mathematics, University of Washington, Box 354350, Seattle, Washington 98195, USA
| | - Robert Hołyst
- Institute of Physical Chemistry, Polish Academy of Sciences, Kasprzaka 44/52, 01-224 Warsaw, Poland
| |
Collapse
|
12
|
Mahmutovic A, Berg OG, Elf J. What matters for lac repressor search in vivo--sliding, hopping, intersegment transfer, crowding on DNA or recognition? Nucleic Acids Res 2015; 43:3454-64. [PMID: 25779051 PMCID: PMC4402528 DOI: 10.1093/nar/gkv207] [Citation(s) in RCA: 40] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2015] [Accepted: 02/28/2015] [Indexed: 01/04/2023] Open
Abstract
We have investigated which aspects of transcription factor DNA interactions are most important to account for the recent in vivo search time measurements for the dimeric lac repressor. We find the best agreement for a sliding model where non-specific binding to DNA is improbable at first contact and the sliding LacI protein binds at high probability when reaching the specific Osym operator. We also find that the contribution of hopping to the overall search speed is negligible although physically unavoidable. The parameters that give the best fit reveal sliding distances, including hopping, close to what has been proposed in the past, i.e. ∼40 bp, but with an unexpectedly high 1D diffusion constant on non-specific DNA sequences. Including a mechanism of inter-segment transfer between distant DNA segments does not bring down the 1D diffusion to the expected fraction of the in vitro value. This suggests a mechanism where transcription factors can slide less hindered in vivo than what is given by a simple viscosity scaling argument or that a modification of the model is needed. For example, the estimated diffusion rate constant would be consistent with the expectation if parts of the chromosome, away from the operator site, were inaccessible for searching.
Collapse
Affiliation(s)
- Anel Mahmutovic
- Department of Cell and Molecular Biology, Science for Life Laboratory, Uppsala University, 75124 Uppsala, Sweden
| | - Otto G Berg
- Department of Cell and Molecular Biology, Science for Life Laboratory, Uppsala University, 75124 Uppsala, Sweden
| | - Johan Elf
- Department of Cell and Molecular Biology, Science for Life Laboratory, Uppsala University, 75124 Uppsala, Sweden
| |
Collapse
|
13
|
Furini S, Domene C. DNA recognition process of the lactose repressor protein studied via metadynamics and umbrella sampling simulations. J Phys Chem B 2014; 118:13059-65. [PMID: 25341013 DOI: 10.1021/jp505885j] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The lactose repressor, LacI, finds its DNA target sites via a process that is faster than what it is expected from a diffusion-driven mechanism. This is possible thanks to nonspecific binding of LacI to DNA, followed by diffusion along the DNA molecule. The diffusion of the protein along DNA might lead to a fast-searching mechanism only if LacI binds with comparable strength to different nonspecific sequences and if, in addition, the value of the binding energy remarkably decreases in the presence of a binding site. The first condition would be favored by loose interactions with the base edges, while the second would take advantage from the opposite situation. In order to understand how the protein satisfies these two opposing requirements, the DNA recognition process was studied by a combination of umbrella sampling and metadynamics simulations. The simulations revealed that when aligned with a specific sequence, LacI establishes polar interactions with the base edges that require ∼4 kcal/mol to be disrupted. In contrast, these interactions are not stable when the protein is aligned with nonspecific sequences. These results confirm that LacI is able to efficiently recognize a specific sequence while sliding along DNA before any structural change of the protein-DNA complex occurs.
Collapse
Affiliation(s)
- Simone Furini
- Department of Medical Biotechnologies, University of Siena , viale Mario Bracci 16, I-53100, Siena, Italy
| | | |
Collapse
|
14
|
van der Vaart A. Coupled binding-bending-folding: The complex conformational dynamics of protein-DNA binding studied by atomistic molecular dynamics simulations. Biochim Biophys Acta Gen Subj 2014; 1850:1091-1098. [PMID: 25161164 DOI: 10.1016/j.bbagen.2014.08.009] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/25/2014] [Revised: 08/14/2014] [Accepted: 08/18/2014] [Indexed: 12/21/2022]
Abstract
BACKGROUND Protein-DNA binding often involves dramatic conformational changes such as protein folding and DNA bending. While thermodynamic aspects of this behavior are understood, and its biological function is often known, the mechanism by which the conformational changes occur is generally unclear. By providing detailed structural and energetic data, molecular dynamics simulations have been helpful in elucidating and rationalizing protein-DNA binding. SCOPE OF REVIEW This review will summarize recent atomistic molecular dynamics simulations of the conformational dynamics of DNA and protein-DNA binding. A brief overview of recent developments in DNA force fields is given as well. MAJOR CONCLUSIONS Simulations have been crucial in rationalizing the intrinsic flexibility of DNA, and have been instrumental in identifying the sequence of binding events, the triggers for the conformational motion, and the mechanism of binding for a number of important DNA-binding proteins. GENERAL SIGNIFICANCE Molecular dynamics simulations are an important tool for understanding the complex binding behavior of DNA-binding proteins. With recent advances in force fields and rapid increases in simulation time scales, simulations will become even more important for future studies. This article is part of a Special Issue entitled Recent developments of molecular dynamics.
Collapse
Affiliation(s)
- Arjan van der Vaart
- Department of Chemistry, University of South Florida, 4202 East Fowler Avenue CHE 205, Tampa, FL 33620, USA.
| |
Collapse
|
15
|
Marklund EG, Mahmutovic A, Berg OG, Hammar P, van der Spoel D, Fange D, Elf J. Transcription-factor binding and sliding on DNA studied using micro- and macroscopic models. Proc Natl Acad Sci U S A 2013; 110:19796-801. [PMID: 24222688 PMCID: PMC3856812 DOI: 10.1073/pnas.1307905110] [Citation(s) in RCA: 63] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023] Open
Abstract
Transcription factors search for specific operator sequences by alternating rounds of 3D diffusion with rounds of 1D diffusion (sliding) along the DNA. The details of such sliding have largely been beyond direct experimental observation. For this purpose we devised an analytical formulation of umbrella sampling along a helical coordinate, and from extensive and fully atomistic simulations we quantified the free-energy landscapes that underlie the sliding dynamics and dissociation kinetics for the LacI dimer. The resulting potential of mean force distributions show a fine structure with an amplitude of 1 k(B)T for sliding and 12 k(B)T for dissociation. Based on the free-energy calculations the repressor slides in close contact with DNA for 8 bp on average before making a microscopic dissociation. By combining the microscopic molecular-dynamics calculations with Brownian simulation including rotational diffusion from the microscopically dissociated state we estimate a macroscopic residence time of 48 ms at the same DNA segment and an in vitro sliding distance of 240 bp. The sliding distance is in agreement with previous in vitro sliding-length estimates. The in vitro prediction for the macroscopic residence time also compares favorably to what we measure by single-molecule imaging of nonspecifically bound fluorescently labeled LacI in living cells. The investigation adds to our understanding of transcription-factor search kinetics and connects the macro-/mesoscopic rate constants to the microscopic dynamics.
Collapse
Affiliation(s)
- Erik G. Marklund
- Department of Cell and Molecular Biology, Science for Life Laboratory, Uppsala University, 75124 Uppsala, Sweden
| | - Anel Mahmutovic
- Department of Cell and Molecular Biology, Science for Life Laboratory, Uppsala University, 75124 Uppsala, Sweden
| | - Otto G. Berg
- Department of Cell and Molecular Biology, Science for Life Laboratory, Uppsala University, 75124 Uppsala, Sweden
| | - Petter Hammar
- Department of Cell and Molecular Biology, Science for Life Laboratory, Uppsala University, 75124 Uppsala, Sweden
| | - David van der Spoel
- Department of Cell and Molecular Biology, Science for Life Laboratory, Uppsala University, 75124 Uppsala, Sweden
| | - David Fange
- Department of Cell and Molecular Biology, Science for Life Laboratory, Uppsala University, 75124 Uppsala, Sweden
| | - Johan Elf
- Department of Cell and Molecular Biology, Science for Life Laboratory, Uppsala University, 75124 Uppsala, Sweden
| |
Collapse
|
16
|
Yonetani Y, Kono H. Dissociation Free-Energy Profiles of Specific and Nonspecific DNA–Protein Complexes. J Phys Chem B 2013; 117:7535-45. [DOI: 10.1021/jp402664w] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Affiliation(s)
- Yoshiteru Yonetani
- Molecular Modeling and Simulation Group, Quantum Beam
Science Directorate, Japan Atomic Energy Agency, 8-1-7 Umemidai, Kizugawa, Kyoto 619-0215, Japan
| | - Hidetoshi Kono
- Molecular Modeling and Simulation Group, Quantum Beam
Science Directorate, Japan Atomic Energy Agency, 8-1-7 Umemidai, Kizugawa, Kyoto 619-0215, Japan
| |
Collapse
|
17
|
Furini S, Barbini P, Domene C. DNA-recognition process described by MD simulations of the lactose repressor protein on a specific and a non-specific DNA sequence. Nucleic Acids Res 2013; 41:3963-72. [PMID: 23430151 PMCID: PMC3627591 DOI: 10.1093/nar/gkt099] [Citation(s) in RCA: 43] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023] Open
Abstract
The lactose repressor protein may bind DNA in two possible configurations: a specific one, if the DNA sequence corresponds to a binding site, and a non-specific one otherwise. To find its target sequences, the lactose repressor first binds non-specifically to DNA, and subsequently, it rapidly searches for a binding site. Atomic structures of non-specific and specific complexes are available from crystallographic and nuclear magnetic resonance experiments. However, what remains unknown is a detailed description of the steps that transform the non-specific complex into the specific one. Here, how the protein first recognizes its binding site has been studied using molecular dynamics simulations. The picture that emerges is that of a protein that is as mobile when interacting with non-specific DNA sequences as when free in solution. This high degree of mobility allows the protein to rapidly sample different DNA sequences. In contrast, when the protein encounters a binding site, the configuration ensemble collapses, and the protein sliding movements along the DNA sequence become scarce. The binding energies in the specific and non-specific complexes were analysed using the Molecular Mechanics Poisson Boltzmann Surface Area approach. These results represent a first step towards a throughout characterization of the DNA-recognition process.
Collapse
Affiliation(s)
- Simone Furini
- Department of Medical Biotechnology, University of Siena, viale Mario Bracci 12, I-53100 Siena, Italy.
| | | | | |
Collapse
|
18
|
Optical Methods to Study Protein-DNA Interactions in Vitro and in Living Cells at the Single-Molecule Level. Int J Mol Sci 2013; 14:3961-92. [PMID: 23429188 PMCID: PMC3588080 DOI: 10.3390/ijms14023961] [Citation(s) in RCA: 39] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/15/2012] [Revised: 01/13/2013] [Accepted: 02/04/2013] [Indexed: 12/13/2022] Open
Abstract
The maintenance of intact genetic information, as well as the deployment of transcription for specific sets of genes, critically rely on a family of proteins interacting with DNA and recognizing specific sequences or features. The mechanisms by which these proteins search for target DNA are the subject of intense investigations employing a variety of methods in biology. A large interest in these processes stems from the faster-than-diffusion association rates, explained in current models by a combination of 3D and 1D diffusion. Here, we present a review of the single-molecule approaches at the forefront of the study of protein-DNA interaction dynamics and target search in vitro and in vivo. Flow stretch, optical and magnetic manipulation, single fluorophore detection and localization as well as combinations of different methods are described and the results obtained with these techniques are discussed in the framework of the current facilitated diffusion model.
Collapse
|
19
|
Haeusler AR, Goodson KA, Lillian TD, Wang X, Goyal S, Perkins NC, Kahn JD. FRET studies of a landscape of Lac repressor-mediated DNA loops. Nucleic Acids Res 2012; 40:4432-45. [PMID: 22307389 PMCID: PMC3378866 DOI: 10.1093/nar/gks019] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022] Open
Abstract
DNA looping mediated by the Lac repressor is an archetypal test case for modeling protein and DNA flexibility. Understanding looping is fundamental to quantitative descriptions of gene expression. Systematic analysis of LacI•DNA looping was carried out using a landscape of DNA constructs with lac operators bracketing an A-tract bend, produced by varying helical phasings between operators and the bend. Fluorophores positioned on either side of both operators allowed direct Förster resonance energy transfer (FRET) detection of parallel (P1) and antiparallel (A1, A2) DNA looping topologies anchored by V-shaped LacI. Combining fluorophore position variant landscapes allows calculation of the P1, A1 and A2 populations from FRET efficiencies and also reveals extended low-FRET loops proposed to form via LacI opening. The addition of isopropyl-β-d-thio-galactoside (IPTG) destabilizes but does not eliminate the loops, and IPTG does not redistribute loops among high-FRET topologies. In some cases, subsequent addition of excess LacI does not reduce FRET further, suggesting that IPTG stabilizes extended or other low-FRET loops. The data align well with rod mechanics models for the energetics of DNA looping topologies. At the peaks of the predicted energy landscape for V-shaped loops, the proposed extended loops are more stable and are observed instead, showing that future models must consider protein flexibility.
Collapse
Affiliation(s)
- Aaron R Haeusler
- Department of Chemistry and Biochemistry, University of Maryland, College Park, MD 20742-2021, USA
| | | | | | | | | | | | | |
Collapse
|
20
|
Foffano G, Marenduzzo D, Orlandini E. Facilitated diffusion on confined DNA. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2012; 85:021919. [PMID: 22463256 DOI: 10.1103/physreve.85.021919] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/13/2011] [Revised: 12/12/2011] [Indexed: 05/31/2023]
Abstract
In living cells, proteins combine three-dimensional bulk diffusion and one-dimensional sliding along the DNA to reach a target faster. This process is known as facilitated diffusion and we investigate its dynamics in the physiologically relevant case of confined DNA. The confining geometry and DNA elasticity are key parameters: We find that facilitated diffusion is most efficient inside an isotropic volume and on a flexible polymer. By considering the typical copy numbers of proteins in vivo, we show that the speedup due to sliding becomes insensitive to fine tuning of parameters, rendering facilitated diffusion a robust mechanism to speed up intracellular diffusion-limited reactions. The parameter range we focus on is relevant for in vitro systems and for facilitated diffusion on yeast chromatin.
Collapse
Affiliation(s)
- G Foffano
- School of Physics and Astronomy, University of Edinburgh, Mayfield Road, Edinburgh EH9 3JZ, United Kingdom
| | | | | |
Collapse
|
21
|
Towards a molecular view of transcriptional control. Curr Opin Struct Biol 2012; 22:160-7. [PMID: 22296921 DOI: 10.1016/j.sbi.2012.01.004] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/03/2011] [Revised: 01/06/2012] [Accepted: 01/09/2012] [Indexed: 11/21/2022]
Abstract
The accumulation of experimental data over recent years has fueled theoretical work on how transcription factors (TFs) search for and recognise their DNA target sites, how they interact with one another, or with other DNA-binding proteins, and how they cope with the compaction of DNA within bacterial nucleoids or within eukaryotic chromatin. Many models have been built to study the kinetic, thermodynamic and mechanistic aspects of these questions. In some cases they have resulted in a relatively clear consensus view, but a number of questions remain controversial. We present an overview of recent work, with an emphasis on models that provide, or can inspire, a better understanding of transcriptional control at a detailed molecular level.
Collapse
|