Structural Mechanism behind Distinct Efficiency of Oct4/Sox2 Proteins in Differentially Spaced DNA Complexes

The octamer-binding transcription factor 4 (Oct4) and sex-determining region Y (SRY)-box 2 (Sox2) proteins induce various transcriptional regulators to maintain cellular pluripotency. Most Oct4/Sox2 complexes have either 0 base pairs (Oct4/Sox20bp) or 3 base pairs (Oct4/Sox23bp) separation between their DNA-binding sites. Results from previous biochemical studies have shown that the complexes separated by 0 base pairs are associated with a higher pluripotency rate than those separated by 3 base pairs. Here, we performed molecular dynamics (MD) simulations and calculations to determine the binding free energy and per-residue free energy for the Oct4/Sox20bp and Oct4/Sox23bp complexes to identify structural differences that contribute to differences in induction rate. Our MD simulation results showed substantial differences in Oct4/Sox2 domain movements, as well as secondary-structure changes in the Oct4 linker region, suggesting a potential reason underlying the distinct efficiencies of these complexes during reprogramming. Moreover, we identified key residues and hydrogen bonds that potentially facilitate protein-protein and protein-DNA interactions, in agreement with previous experimental findings. Consequently, our results confess that differential spacing of the Oct4/Sox2 DNA binding sites can determine the magnitude of transcription of the targeted genes during reprogramming.

[1]  P Stothard,et al.  The sequence manipulation suite: JavaScript programs for analyzing and formatting protein and DNA sequences. , 2000, BioTechniques.

[2]  H. Kondoh,et al.  Pairing SOX off: with partners in the regulation of embryonic development. , 2000, Trends in genetics : TIG.

[3]  Ray Luo,et al.  Molecular dynamics simulations of p53 DNA-binding domain. , 2007, The journal of physical chemistry. B.

[4]  Sangdun Choi,et al.  Structure-Function Relationship of Cytoplasmic and Nuclear IκB Proteins: An In Silico Analysis , 2010, PloS one.

[5]  David J. Arenillas,et al.  JASPAR 2014: an extensively expanded and updated open-access database of transcription factor binding profiles , 2013, Nucleic Acids Res..

[6]  Christophe Blanchet,et al.  CURVES+ web server for analyzing and visualizing the helical, backbone and groove parameters of nucleic acid structures , 2011, Nucleic Acids Res..

[7]  Conrad C. Huang,et al.  UCSF Chimera—A visualization system for exploratory research and analysis , 2004, J. Comput. Chem..

[8]  D. Tantin,et al.  Oct transcription factors in development and stem cells: insights and mechanisms , 2013, Development.

[9]  Debra L. Fulton,et al.  The Transcription Factor Encyclopedia , 2012, Genome Biology.

[10]  Hiroshi Handa,et al.  A general mechanism for transcription regulation by Oct1 and Oct4 in response to genotoxic and oxidative stress. , 2009, Genes & development.

[11]  P. Robson,et al.  Transcriptional Regulation of Nanog by OCT4 and SOX2* , 2005, Journal of Biological Chemistry.

[12]  Olivier Michielin,et al.  Comparison between computational alanine scanning and per‐residue binding free energy decomposition for protein–protein association using MM‐GBSA: Application to the TCR‐p‐MHC complex , 2007, Proteins.

[13]  Stefan Fischer,et al.  Analyzing large‐scale structural change in proteins: Comparison of principal component projection and sammon mapping , 2006, Proteins.

[14]  Giovanni Chillemi,et al.  Thr729 in human topoisomerase I modulates anti-cancer drug resistance by altering protein domain communications as suggested by molecular dynamics simulations , 2008, Nucleic acids research.

[15]  K Schulten,et al.  VMD: visual molecular dynamics. , 1996, Journal of molecular graphics.

[16]  A. Wolffe,et al.  The structure of DNA in a nucleosome. , 1990, Proceedings of the National Academy of Sciences of the United States of America.

[17]  N. D. Clarke,et al.  Integration of External Signaling Pathways with the Core Transcriptional Network in Embryonic Stem Cells , 2008, Cell.

[18]  Junmei Wang,et al.  Development and testing of a general amber force field , 2004, J. Comput. Chem..

[19]  R. Mann,et al.  The role of DNA shape in protein-DNA recognition , 2009, Nature.

[20]  L. Nilsson,et al.  Conformational states of the glucocorticoid receptor DNA‐binding domain from molecular dynamics simulations , 2002, Proteins.

[21]  J. H. Maddocks,et al.  Conformational analysis of nucleic acids revisited: Curves+ , 2009, Nucleic acids research.

[22]  Megan F. Cole,et al.  Core Transcriptional Regulatory Circuitry in Human Embryonic Stem Cells , 2005, Cell.

[23]  Marcos J. Araúzo-Bravo,et al.  A unique Oct4 interface is crucial for reprogramming to pluripotency , 2013, Nature Cell Biology.

[24]  G. Marius Clore,et al.  Molecular Basis for Synergistic Transcriptional Activation by Oct1 and Sox2 Revealed from the Solution Structure of the 42-kDa Oct1·Sox2·Hoxb1-DNA Ternary Transcription Factor Complex* , 2004, Journal of Biological Chemistry.

[25]  Peter M. Kasson,et al.  GROMACS 4.5: a high-throughput and highly parallel open source molecular simulation toolkit , 2013, Bioinform..

[26]  Jiaxuan Chen,et al.  Conversion of Sox17 into a Pluripotency Reprogramming Factor by Reengineering Its Association with Oct4 on DNA , 2011, Stem cells.

[27]  Felipe Merino,et al.  Structural basis for the SOX-dependent genomic redistribution of OCT4 in stem cell differentiation. , 2014, Structure.

[28]  R. Dickerson,et al.  Structure of the B-DNA decamer C-C-A-A-C-G-T-T-G-G and comparison with isomorphous decamers C-C-A-A-G-A-T-T-G-G and C-C-A-G-G-C-C-T-G-G. , 1991, Journal of molecular biology.

[29]  P Shing Ho,et al.  How sequence defines structure: a crystallographic map of DNA structure and conformation. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[30]  Felipe Merino,et al.  Cooperative DNA Recognition Modulated by an Interplay between Protein-Protein Interactions and DNA-Mediated Allostery , 2015, PLoS Comput. Biol..

[31]  Tingjun Hou,et al.  Assessing the Performance of the MM/PBSA and MM/GBSA Methods. 1. The Accuracy of Binding Free Energy Calculations Based on Molecular Dynamics Simulations , 2011, J. Chem. Inf. Model..

[32]  T Darden,et al.  New tricks for modelers from the crystallography toolkit: the particle mesh Ewald algorithm and its use in nucleic acid simulations. , 1999, Structure.

[33]  Sangdun Choi,et al.  Molecular Modeling-Based Evaluation of hTLR10 and Identification of Potential Ligands in Toll-Like Receptor Signaling , 2010, PloS one.

[34]  Balachandran Manavalan,et al.  Random Forest-Based Protein Model Quality Assessment (RFMQA) Using Structural Features and Potential Energy Terms , 2014, PloS one.

[35]  A. Travers,et al.  DNA conformation and protein binding. , 1989, Annual review of biochemistry.

[36]  T. Richmond,et al.  The structure of DNA in the nucleosome core , 2003, Nature.

[37]  Prasanna R Kolatkar,et al.  The structure of Sox17 bound to DNA reveals a conserved bending topology but selective protein interaction platforms. , 2009, Journal of molecular biology.

[38]  Leo S. D. Caves,et al.  Bio3d: An R Package , 2022 .

[39]  K. Plath,et al.  The roles of the reprogramming factors Oct4, Sox2 and Klf4 in resetting the somatic cell epigenome during induced pluripotent stem cell generation , 2012, Genome Biology.

[40]  Qingchuan Zheng,et al.  Insights into the effects of mutations on Cren7-DNA binding using molecular dynamics simulations and free energy calculations. , 2015, Physical chemistry chemical physics : PCCP.

[41]  S. Upadhyay Dynamics of Gal80p in the Gal80p-Gal3p complex differ significantly from the dynamics in the Gal80p-Gal1p complex: implications for the higher specificity of Gal3p. , 2014, Molecular bioSystems.

[42]  H. Schöler,et al.  OCT4: dynamic DNA binding pioneers stem cell pluripotency. , 2014, Biochimica et biophysica acta.

[43]  Matthias Wilmanns,et al.  Crystal structure of a POU/HMG/DNA ternary complex suggests differential assembly of Oct4 and Sox2 on two enhancers. , 2003, Genes & development.

[44]  E. Papaleo,et al.  Intramolecular interactions stabilizing compact conformations of the intrinsically disordered kinase-inhibitor domain of Sic1: a molecular dynamics investigation , 2012, Front. Physio..

[45]  R. Dror,et al.  Improved side-chain torsion potentials for the Amber ff99SB protein force field , 2010, Proteins.

[46]  Holger Gohlke,et al.  MMPBSA.py: An Efficient Program for End-State Free Energy Calculations. , 2012, Journal of chemical theory and computation.

[47]  H. Kondoh,et al.  Sox proteins: regulators of cell fate specification and differentiation , 2013, Development.

[48]  V. Hilser,et al.  Structural Dynamics, Intrinsic Disorder, and Allostery in Nuclear Receptors as Transcription Factors* , 2011, The Journal of Biological Chemistry.

[49]  Berk Hess,et al.  LINCS: A linear constraint solver for molecular simulations , 1997, J. Comput. Chem..

[50]  J. Berg,et al.  Molecular dynamics simulations of biomolecules , 2002, Nature Structural Biology.

[51]  V. Hornak,et al.  Comparison of multiple Amber force fields and development of improved protein backbone parameters , 2006, Proteins.

[52]  M. Bianchi,et al.  Spatially Precise DNA Bending Is an Essential Activity of the Sox2 Transcription Factor* , 2001, The Journal of Biological Chemistry.

[53]  Sándor Pongor,et al.  DNA analysis servers: plot.it, bend.it, model.it and IS , 2003, Nucleic Acids Res..

[54]  R. Nussinov,et al.  p53-Induced DNA bending: the interplay between p53-DNA and p53-p53 interactions. , 2008, The journal of physical chemistry. B.

[55]  A. Liwo,et al.  Principal component analysis for protein folding dynamics. , 2009, Journal of molecular biology.

[56]  H. Berendsen,et al.  Essential dynamics of proteins , 1993, Proteins.

[57]  T. Ichisaka,et al.  Induction of Pluripotent Stem Cells from Adult Human Fibroblasts by Defined Factors , 2007, Cell.

[58]  H. Schöler,et al.  Stem cell pluripotency and transcription factor Oct4 , 2002, Cell Research.

[59]  Rajendra Kumar,et al.  g_mmpbsa - A GROMACS Tool for High-Throughput MM-PBSA Calculations , 2014, J. Chem. Inf. Model..

[60]  Debbie L C van den Berg,et al.  An Oct4-Centered Protein Interaction Network in Embryonic Stem Cells , 2010, Cell stem cell.