1
|
van Heesch T, van de Lagemaat EM, Vreede J. Deciphering Sequence-Specific DNA Binding by H-NS Using Molecular Simulation. Methods Mol Biol 2024; 2819:585-609. [PMID: 39028525 DOI: 10.1007/978-1-0716-3930-6_27] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/20/2024]
Abstract
H-NS is a DNA organizing protein that occurs in Gram-negative bacteria. It can form long filaments between two DNA duplexes by first binding to a high-affinity AT-rich nucleotide sequence and extending from there. Using molecular dynamics simulations and steered molecular dynamics, we are able to determine the free energy of formation and dissociation of a protein-DNA complex comprising an H-NS DNA-binding domain and a specific nucleotide sequence. The molecular dynamics simulations allow detailed characterization of the interactions between the protein and a specific nucleotide sequence. To quantify the strength of the interaction, we employ an additional potential based on protein-DNA contacts to speed up dissociation of the protein-DNA complex. The work required for the dissociation results in an estimate of the free energy of dissociation/complex formation. Our protocol can provide quantitative prediction of protein-DNA complex stability, while also providing high-resolution insights into recognition mechanisms. In this chapter, we have used this approach to quantify the sequence specificity of H-NS DNA-binding domains to various nucleotide sequences, thus elucidating the mechanism with which H-NS can specifically bind to AT-rich DNA.
Collapse
Affiliation(s)
- Thor van Heesch
- van 't Hoff Institute for Molecular Sciences, University of Amsterdam, Amsterdam, The Netherlands
| | - Eline M van de Lagemaat
- van 't Hoff Institute for Molecular Sciences, University of Amsterdam, Amsterdam, The Netherlands
| | - Jocelyne Vreede
- van 't Hoff Institute for Molecular Sciences, University of Amsterdam, Amsterdam, The Netherlands.
| |
Collapse
|
2
|
van Heesch T, Bolhuis PG, Vreede J. Decoding dissociation of sequence-specific protein-DNA complexes with non-equilibrium simulations. Nucleic Acids Res 2023; 51:12150-12160. [PMID: 37953329 PMCID: PMC10711434 DOI: 10.1093/nar/gkad1014] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/27/2023] [Revised: 10/13/2023] [Accepted: 10/19/2023] [Indexed: 11/14/2023] Open
Abstract
Sequence-specific protein-DNA interactions are crucial in processes such as DNA organization, gene regulation and DNA replication. Obtaining detailed insights into the recognition mechanisms of protein-DNA complexes through experiments is hampered by a lack of resolution in both space and time. Here, we present a molecular simulation approach to quantify the sequence specificity of protein-DNA complexes, that yields results fast, and is generally applicable to any protein-DNA complex. The approach is based on molecular dynamics simulations in combination with a sophisticated steering potential and results in an estimate of the free energy difference of dissociation. We provide predictions of the nucleotide specific binding affinity of the minor groove binding Histone-like Nucleoid Structuring (H-NS) protein, that are in agreement with experimental data. Furthermore, our approach offers mechanistic insight into the process of dissociation. Applying our approach to the major groove binding ETS domain in complex with three different nucleotide sequences identified the high affinity consensus sequence, quantitatively in agreement with experiments. Our protocol facilitates quantitative prediction of protein-DNA complex stability, while also providing high resolution insights into recognition mechanisms. As such, our simulation approach has the potential to yield detailed and quantitative insights into biological processes involving sequence-specific protein-DNA interactions.
Collapse
Affiliation(s)
- Thor van Heesch
- Van ’t Hoff Institute for Molecular Sciences, University of Amsterdam, Netherlands
| | - Peter G Bolhuis
- Van ’t Hoff Institute for Molecular Sciences, University of Amsterdam, Netherlands
| | - Jocelyne Vreede
- Van ’t Hoff Institute for Molecular Sciences, University of Amsterdam, Netherlands
| |
Collapse
|
3
|
Al Masri C, Wan B, Yu J. Nonspecific vs. specific DNA binding free energetics of a transcription factor domain protein. Biophys J 2023; 122:4476-4487. [PMID: 37897044 PMCID: PMC10722393 DOI: 10.1016/j.bpj.2023.10.025] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2023] [Revised: 07/06/2023] [Accepted: 10/20/2023] [Indexed: 10/29/2023] Open
Abstract
Transcription factor (TF) proteins regulate gene expression by binding to specific sites on the genome. In the facilitated diffusion model, an optimized search process is achieved by the TF alternating between 3D diffusion in the bulk and 1D diffusion along DNA. While undergoing 1D diffusion, the protein can switch from a search mode for fast diffusion along nonspecific DNA to a recognition mode for stable binding to specific DNA. It was recently noticed that, for a small TF domain protein, reorientations on DNA happen between the nonspecific and specific DNA binding. We here conducted all-atom molecular dynamics simulations with steering forces to reveal the protein-DNA binding free energetics, confirming that the search and recognition modes are distinguished primarily by protein orientations on the DNA. As the binding free energy difference between the specific and nonspecific DNA system slightly deviates from that being estimated directly from dissociation constants on 15-bp DNA constructs, we hypothesize that the discrepancy can come from DNA sequences flanking the 6-bp central binding sites that impact on the dissociation kinetics measurements. The hypothesis is supported by a simplified spherical protein-DNA model along with stochastic simulations and kinetic modeling.
Collapse
Affiliation(s)
- Carmen Al Masri
- Department of Physics and Astronomy, University of California, Irvine, California
| | - Biao Wan
- Wenzhou Institute, University of Chinese Academy of Sciences, Wenzhou, China
| | - Jin Yu
- Department of Physics and Astronomy, University of California, Irvine, California; Department of Physics and Astronomy, Department of Chemistry, NSF-Simons Center for Multiscale Cell Fate Research, University of California, Irvine, California.
| |
Collapse
|
4
|
Singh RK, Mukherjee A. Molecular Mechanism of the Intercalation of the SOX-4 Protein into DNA Inducing Bends and Kinks. J Phys Chem B 2021; 125:3752-3762. [PMID: 33848164 DOI: 10.1021/acs.jpcb.0c11496] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Abstract
DNA-protein interactions regulate several biophysical functions, yet the mechanism of only a few is investigated in molecular detail. An important example is the intercalation of transcription factor proteins into DNA that produce bent and kinked DNA. Here, we have studied the molecular mechanism of the intercalation of a transcription factor SOX4 into DNA with a goal to understand the sequence of molecular events that precede the bending and kinking of the DNA. Our long well-tempered metadynamics and molecular dynamics (MD) simulations show that the protein primarily binds to the backbone of DNA and rotates around it to form an intercalative native state. We show that although there are multiple pathways for intercalation, the deintercalation pathway matches with the most probable intercalation pathway. In both cases, bending and kinking happen simultaneously, driven by the onset of the intercalation of the amino acid.
Collapse
Affiliation(s)
- Reman Kumar Singh
- Department of Chemistry, Indian Institute of Science Education and Research, Pune 411008, India
| | - Arnab Mukherjee
- Department of Chemistry, Indian Institute of Science Education and Research, Pune 411008, India
| |
Collapse
|
5
|
Chen J, Qi Y, Duan Y, Duan M, Yang M. C1188D mutation abolishes specific recognition between MLL1-CXXC domain and CpG site by inducing conformational switch of flexible N-terminal. Proteins 2020; 88:1401-1412. [PMID: 32519403 DOI: 10.1002/prot.25960] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/05/2019] [Revised: 05/22/2020] [Accepted: 06/06/2020] [Indexed: 01/19/2023]
Abstract
Mixed lineage leukemia protein (MLL1 protein) recognizes the CpG site via its CXXC domain and is frequently associated with leukemia. The specific recognition is abolished by C1188D mutation, which also prevents MLL-related leukemia. In this paper, multiple molecular dynamic (MD) simulations were performed to investigate the mechanism of recognition and influences of C1188D mutation. Started from fully dissociated DNA and MLL1-CXXC domain, remarkably, the center of mass (COM) of MLL1-CXXC domain quickly concentrates on the vicinity of the CpG site in all 53 short MD simulations. Extended simulations of the wild type showed that the native complex formed in 500 ns among 4 of 53 simulations. In contrast, the C1188D mutant COM distributed broadly around the DNA and the native complex was not observed in any of the extended simulations. Simulations on the apo MLL1-CXXC domain further suggest that the wild type protein remained predominantly in an open form that closely resembles its structure in the native complex whereas C1188D mutant formed predominantly compact structures in which the N- terminal bends to D1188. This conformational switch hinders the formation of encounter complex, thus abolishes the recognition. Our study also provides clues to the study mechanism of recognition, by the CXXC domain from proteins like DNA methyltransferase and ten-eleven translocation enzymes.
Collapse
Affiliation(s)
- Jiawen Chen
- Key Laboratory of Magnetic Resonance in Biological Systems, National Center for Magnetic Resonances in Wuhan, State Key laboratory of Magnetic Resonance and Atomic and Molecular Physics, Wuhan Institute of Physics and Mathematics, Innovation Academy for Precision Measurement Science and Technology, Chinese Academy of Sciences, Wuhan, China
| | - Yanping Qi
- Key Laboratory of Magnetic Resonance in Biological Systems, National Center for Magnetic Resonances in Wuhan, State Key laboratory of Magnetic Resonance and Atomic and Molecular Physics, Wuhan Institute of Physics and Mathematics, Innovation Academy for Precision Measurement Science and Technology, Chinese Academy of Sciences, Wuhan, China.,College of Physical Science and Technology, Central China Normal University, Wuhan, China
| | - Yong Duan
- Department of Biomedical Engineering and UC Davis Genome Center, University of California at Davis, Davis, California, USA
| | - Mojie Duan
- Key Laboratory of Magnetic Resonance in Biological Systems, National Center for Magnetic Resonances in Wuhan, State Key laboratory of Magnetic Resonance and Atomic and Molecular Physics, Wuhan Institute of Physics and Mathematics, Innovation Academy for Precision Measurement Science and Technology, Chinese Academy of Sciences, Wuhan, China
| | - Minghui Yang
- Key Laboratory of Magnetic Resonance in Biological Systems, National Center for Magnetic Resonances in Wuhan, State Key laboratory of Magnetic Resonance and Atomic and Molecular Physics, Wuhan Institute of Physics and Mathematics, Innovation Academy for Precision Measurement Science and Technology, Chinese Academy of Sciences, Wuhan, China.,Wuhan National Laboratory for Optoelectronics, Huazhong University of Science and Technology, Wuhan, China
| |
Collapse
|
6
|
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
|
7
|
Jakubec D, Vondrášek J. Can All-Atom Molecular Dynamics Simulations Quantitatively Describe Homeodomain-DNA Binding Equilibria? J Chem Theory Comput 2019; 15:2635-2648. [PMID: 30807142 DOI: 10.1021/acs.jctc.8b01144] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
We systematically investigate the applicability of a molecular dynamics-based setup for the calculations of standard binding free energies of biologically relevant protein-DNA complexes. The free energies are extracted from a potential of mean force calculated using umbrella sampling simulations. Two protein-DNA systems derived from a homeodomain transcription factor complex are studied in order to investigate the binding of both disordered and globular proteins. Free energies and trajectories obtained using two modern molecular mechanical force fields are compared to each other and to experimental data. The temperature dependence of the calculated standard binding free energies is investigated by performing all simulations over a range of temperatures. We show that the values of standard binding free energies obtained from these simulations are overestimated compared to experimental results. Significant differences are observed between the two protein-DNA systems and between the two force fields, which are explained by different propensities to form inter- and intramolecular contacts. The number of protein-DNA contacts increases with increasing temperature, in agreement with the experimentally known temperature dependence of enthalpies of binding. However, conclusions about the temperature dependence of the standard binding free energies cannot be made with confidence, as the differences among the values are on the order of statistical uncertainty.
Collapse
Affiliation(s)
- David Jakubec
- Bioinformatics Group, Institute of Organic Chemistry and Biochemistry , Czech Academy of Sciences , 166 10 Praha 6, Czech Republic.,Department of Physical and Macromolecular Chemistry, Faculty of Science , Charles University , 128 43 Praha 2, Czech Republic
| | - Jiří Vondrášek
- Bioinformatics Group, Institute of Organic Chemistry and Biochemistry , Czech Academy of Sciences , 166 10 Praha 6, Czech Republic
| |
Collapse
|
8
|
Yesudhas D, Anwar MA, Choi S. Structural mechanism of DNA-mediated Nanog–Sox2 cooperative interaction. RSC Adv 2019; 9:8121-8130. [PMID: 35521171 PMCID: PMC9061787 DOI: 10.1039/c8ra10085c] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/08/2018] [Accepted: 03/04/2019] [Indexed: 01/06/2023] Open
Abstract
The efficiency of stem cell transcriptional regulation always depends on the cooperative association and expression of transcription factors (TFs). Among these, Oct4, Sox2, and Nanog play major roles. Their cooperativity is facilitated via direct protein–protein interactions or DNA-mediated interactions, yet the mechanism is not clear. Most biochemical studies have examined Oct4/Sox2 cooperativity, whereas few studies have evaluated how Nanog competes in the connection between these TFs. In this study, using computational models and molecular dynamics simulations, we built a framework representing the DNA-mediated cooperative interaction between Nanog and Sox2 and analyzed the plausible interaction factors experienced by Nanog because of Sox2, its cooperative binding partner. Comparison of a wild-type and mutant Nanog/Sox2 model with the Nanog crystal structure revealed the regulatory structural mechanism between Nanog/Sox2–DNA-mediated cooperative bindings. Along with the transactivation domains interaction, the DNA-mediated allosteric interactions are also necessary for Nanog cooperative binding. DNA-mediated Nanog–Sox2 cooperativity influences the protein conformational changes and a stronger interaction profile was observed for Nanog-Mut (L103E) in comparison with the Nanog-WT complex. The efficiency of stem cell transcriptional regulation always depends on the cooperative association and expression of transcription factors (TFs).![]()
Collapse
Affiliation(s)
- Dhanusha Yesudhas
- Department of Molecular Science and Technology
- Ajou University
- Suwon
- Korea
| | | | - Sangdun Choi
- Department of Molecular Science and Technology
- Ajou University
- Suwon
- Korea
| |
Collapse
|
9
|
Yonetani Y. Water access and ligand dissociation at the binding site of proteins. J Chem Phys 2018; 149:175102. [PMID: 30408972 DOI: 10.1063/1.5042491] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023] Open
Abstract
Although water is undoubtedly an essential mediator of protein-ligand interactions, whether or not such water molecules are critical for the progress of ligand dissociation remains unclear. To gain a more complete understanding, molecular dynamics simulations are performed with two molecular systems, rigid model binding sites and trypsin-benzamidine. Free-energy landscapes are calculated with a suitably chosen solvent coordinate, which well describes water access to the ligand binding site. The results of free energy provided clear description of water-ligand exchange process, where two different mechanisms appear depending on whether the binding site is buried or not. As the site is more buried, water access is more difficult. When water does not access the site, ligand dissociation produces a large energy barrier, i.e., slow dissociation kinetics. This indicates that control of ligand dissociation kinetics becomes possible with burying the binding site. However, the results also showed that appropriate burying is important because burying reduces not only water access but also ligand binding. The role of the protein structural change is also discussed; it likely plays a similar role to water access because during ligand dissociation, it can make new coordination with the ligand binding site like water. These results contribute to the future pharmaceutical drug design and will be useful for fundamental exploration of various molecular events.
Collapse
Affiliation(s)
- Yoshiteru Yonetani
- Quantum Beam Science Research Directorate, National Institutes for Quantum and Radiological Science and Technology (QST), Tokai-mura, Ibaraki 319-1195, Japan
| |
Collapse
|
10
|
In-depth study of DNA binding of Cys2His2 finger domains in testis zinc-finger protein. PLoS One 2017; 12:e0175051. [PMID: 28384299 PMCID: PMC5383199 DOI: 10.1371/journal.pone.0175051] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2016] [Accepted: 03/20/2017] [Indexed: 11/26/2022] Open
Abstract
Previously, we identified that both fingers 1 and 2 in the three Cys2His2 zinc-finger domains (TZD) of testis zinc-finger protein specifically bind to its cognate DNA; however, finger 3 is non-sequence–specific. To gain insights into the interaction mechanism, here we further investigated the DNA-binding characteristics of TZD bound to non-specific DNAs and its finger segments bound to cognate DNA. TZD in non-specific DNA binding showed smaller chemical shift perturbations, as expected. However, the direction of shift perturbation, change of DNA imino-proton NMR signal, and dynamics on the 15N backbone atom significantly differed between specific and non-specific binding. Using these unique characteristics, we confirmed that the three single-finger segments (TZD1, TZD2 and TZD3) and the two-finger segment (TZD23) non-specifically bind to the cognate DNA. In comparison, the other two-finger segment (TZD12) binding to the cognate DNA features simultaneous non-specific and semi-specific binding, both slowly exchanged in terms of NMR timescale. The process of TZD binding to the cognate DNA is likely stepwise: initially TZD non-specifically binds to DNA, then fingers 1 and 2 insert cooperatively into the major groove of DNA by semi-specific binding, and finally finger 3 non-specifically binds to DNA, which promotes the specific binding on fingers 1 and 2 and stabilizes the formation of a specific TZD–DNA complex.
Collapse
|
11
|
Mechanism of pyrophosphate ion release in T7 RNA polymerase revealed by free energy simulations. COMPUT THEOR CHEM 2016. [DOI: 10.1016/j.comptc.2016.09.005] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
12
|
Dans PD, Walther J, Gómez H, Orozco M. Multiscale simulation of DNA. Curr Opin Struct Biol 2016; 37:29-45. [DOI: 10.1016/j.sbi.2015.11.011] [Citation(s) in RCA: 99] [Impact Index Per Article: 12.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/01/2015] [Revised: 11/23/2015] [Accepted: 11/25/2015] [Indexed: 01/05/2023]
|
13
|
Yonetani Y. Distinct dissociation kinetics between ion pairs: Solvent-coordinate free-energy landscape analysis. J Chem Phys 2015; 143:044506. [DOI: 10.1063/1.4927093] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Yoshiteru Yonetani
- Quantum Beam Science Center, Japan Atomic Energy Agency, 8-1-7 Umemidai, Kizugawa, Kyoto 619-0215, Japan
| |
Collapse
|
14
|
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
|
15
|
Terakawa T, Takada S. RESPAC: Method to Determine Partial Charges in Coarse-Grained Protein Model and Its Application to DNA-Binding Proteins. J Chem Theory Comput 2014; 10:711-21. [DOI: 10.1021/ct4007162] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Tsuyoshi Terakawa
- Department of Biophysics,
Graduate School of Science, Kyoto University, Kitashirakawa Owake-cho, Sakyo, Kyoto, 606-8501, Japan
| | - Shoji Takada
- Department of Biophysics,
Graduate School of Science, Kyoto University, Kitashirakawa Owake-cho, Sakyo, Kyoto, 606-8501, Japan
| |
Collapse
|