DNA binding redistributes activation domain ensemble and accessibility in pioneer factor Sox2

More than 1600 human transcription factors orchestrate the transcriptional machinery to control gene expression and cell fate. Their function is conveyed through intrinsically disordered regions (IDRs) containing activation or repression domains but lacking quantitative structural ensemble models prevents their mechanistic decoding. Here we integrate single-molecule FRET and NMR spectroscopy with molecular simulations showing that DNA binding can lead to complex changes in the IDR ensemble and accessibility. The C-terminal IDR of pioneer factor Sox2 is highly disordered but its conformational dynamics are guided by weak and dynamic charge interactions with the folded DNA binding domain. Both DNA and nucleosome binding induce major rearrangements in the IDR ensemble without affecting DNA binding affinity. Remarkably, interdomain interactions are redistributed in complex with DNA leading to variable exposure of two activation domains critical for transcription. Charged intramolecular interactions allowing for dynamic redistributions may be common in transcription factors and necessary for sensitive tuning of structural ensembles.

Structural Plasticity of Pioneer Factor Sox2 and DNA Bendability Modulate Nucleosome Engagement and Sox2-Oct4 Synergism

Pioneer transcription factors (pTFs) can bind directly to silent chromatin and promote vital transcriptional programs. Here, by integrating high-resolution nuclear magnetic resonance (NMR) spectroscopy with biochemistry, we reveal new structural and mechanistic insights into the interaction of pluripotency pTFs and functional partners Sox2 and Oct4 with nucleosomes. We find that the affinity and conformation of Sox2 for solvent-exposed nucleosome sites depend strongly on their position and DNA sequence. Sox2, which is partially disordered but becomes structured upon DNA binding and bending, forms a super-stable nucleosome complex at superhelical location +5 (SHL+5) with similar affinity and conformation to that with naked DNA. However, at suboptimal internal and end-positioned sites where DNA may be harder to deform, Sox2 favors partially unfolded and more dynamic states that are encoded in its intrinsic flexibility. Importantly, Sox2 structure and DNA bending can be stabilized by synergistic Oct4 binding, but only on adjacent motifs near the nucleosome edge and with the full Oct4 DNA-binding domain. Further mutational studies reveal that strategically impaired Sox2 folding is coupled to reduced DNA bending and inhibits nucleosome binding and Sox2-Oct4 cooperation, while increased nucleosomal DNA flexibility enhances Sox2 association. Together, our findings fit a model where the site-specific DNA bending propensity and structural plasticity of Sox2 govern distinct modes of nucleosome engagement and modulate Sox2-Oct4 synergism. The principles outlined here can potentially guide pTF site selection in the genome and facilitate interaction with other chromatin factors or chromatin opening in vivo.

Reentrant liquid condensate phase of proteins is stabilized by hydrophobic and non-ionic interactions

Liquid-liquid phase separation of proteins underpins the formation of membraneless compartments in living cells. Elucidating the molecular driving forces underlying protein phase transitions is therefore a key objective for understanding biological function and malfunction. Here we show that cellular proteins, which form condensates at low salt concentrations, including FUS, TDP-43, Brd4, Sox2, and Annexin A11, can reenter a phase-separated regime at high salt concentrations. By bringing together experiments and simulations, we demonstrate that this reentrant phase transition in the high-salt regime is driven by hydrophobic and non-ionic interactions, and is mechanistically distinct from the low-salt regime, where condensates are additionally stabilized by electrostatic forces. Our work thus sheds light on the cooperation of hydrophobic and non-ionic interactions as general driving forces in the condensation process, with important implications for aberrant function, druggability, and material properties of biomolecular condensates.

In Silico Identification of SOX1 Post-Translational Modifications Highlights a Shared Protein Motif

The transcription factor SOX1 is a key regulator of neural stem cell development, acting to keep neural stem cells (NSCs) in an undifferentiated state. Postnatal expression of Sox1 is typically confined to the central nervous system (CNS), however, its expression in non-neural tissues has recently been implicated in tumorigenesis. The mechanism through which SOX1 may exert its function is not fully understood, and studies have mainly focused on changes in SOX1 expression at a transcriptional level, while its post-translational regulation remains undetermined. To investigate this, data were extracted from different publicly available databases and analysed to search for putative SOX1 post-translational modifications (PTMs). Results were compared to PTMs associated with SOX2 in order to identify potentially key PTM motifs common to these SOXB1 proteins, and mapped on SOX1 domain structural models. This approach identified several putative acetylation, phosphorylation, glycosylation and sumoylation sites within known functional domains of SOX1. In particular, a novel SOXB1 motif (xKSExSxxP) was identified within the SOX1 protein, which was also found in other unrelated proteins, most of which were transcription factors. These results also highlighted potential phospho-sumoyl switches within this SOXB1 motif identified in SOX1, which could regulate its transcriptional activity. This analysis indicates different types of PTMs within SOX1, which may influence its regulatory role as a transcription factor, by bringing changes to its DNA binding capacities and its interactions with partner proteins. These results provide new research avenues for future investigations on the mechanisms regulating SOX1 activity, which could inform its roles in the contexts of neural stem cell development and cancer.

Nucleosome-bound SOX2 and SOX11 structures elucidate pioneer factor function

'Pioneer' transcription factors are required for stem-cell pluripotency, cell differentiation and cell reprogramming1,2. Pioneer factors can bind nucleosomal DNA to enable gene expression from regions of the genome with closed chromatin. SOX2 is a prominent pioneer factor that is essential for pluripotency and self-renewal of embryonic stem cells3. Here we report cryo-electron microscopy structures of the DNA-binding domains of SOX2 and its close homologue SOX11 bound to nucleosomes. The structures show that SOX factors can bind and locally distort DNA at superhelical location 2. The factors also facilitate detachment of terminal nucleosomal DNA from the histone octamer, which increases DNA accessibility. SOX-factor binding to the nucleosome can also lead to a repositioning of the N-terminal tail of histone H4 that includes residue lysine 16. We speculate that this repositioning is incompatible with higher-order nucleosome stacking, which involves contacts of the H4 tail with a neighbouring nucleosome. Our results indicate that pioneer transcription factors can use binding energy to initiate chromatin opening, and thereby facilitate nucleosome remodelling and subsequent transcription.

Regulation of the Sox3 Gene in an X0/X0 Mammal without Sry, the Amami Spiny Rat, Tokudaia osimensis

Two species of spiny rats, Tokudaia osimensis and Tokudaia tokunoshimensis, show an X0/X0 sex chromosome constitution due to the lack of a Y chromosome. The Sry gene has been completely lost from the genome of these species. We hypothesized that Sox3, which is thought to be originally a homologue of Sry, could function in sex determination in these animals in the absence of Sry. Sox3 was localized in a region of the X chromosome in T. osimensis homologous to mouse. A similar testis- and ovary-specific pattern of expression was observed in mouse and T. osimensis. Although the sequence of the Sox3 gene and its promoter are highly conserved, a 13-bp deletion was specifically found in the promoter region of the 2 spiny rat species. Reporter gene assays were performed to examine the effect of the 13-bp deletion in the promoter region on Sox3 regulation. Although an approximately 60% decrease in activity was observed using the Tokudaia promoters with the 13-bp deletion, the activity was recovered using a mutated promoter in which the deletion was filled with mouse sequence. To evaluate whether SOX3 could regulate Sox9 expression, a reporter gene assay was carried out using testis-specific enhancer of Sox9 core (TESCO). Co-transfection with a combination of mouse SF1 and mouse SOX3 or T. osimensis SOX3 resulted in a greater than 2-fold increase in activity of mouse and T. osimensis TESCO. These results support the idea that the function of SOX3 as a transcription factor, as has been reported in mice and humans, is conserved in T. osimensis. Therefore, we conclude that the Sox3 gene has no function in sex determination in Sry-lacking Tokudaia species.

Increased transactivation and impaired repression of β-catenin-mediated transcription associated with a novel SOX3 missense mutation in an X-linked hypopituitarism pedigree with modest growth failure

SOX3, a transcription factor of the SRY-related high mobility group box family, has been implicated in the etiology of X-linked hypopituitarism. Here, we report a Chinese pedigree of X-linked hypopituitarism with variable phenotypes. Despite the complete growth hormone deficiency, the growth failure of the patients was relatively modest. A rare point variant of SOX3 (c.424C > A; p. P142T) was identified in the pedigree via target panel sequencing. An in vitro study showed that both the expression and nuclear targeting of SOX3 remained unaffected by the variant. However, increased transcriptional activation and impaired repression of β-catenin-mediated transcription were noticed as a result of the SOX3 variant. This is the first study to report that the rare SOX3 missense variant associated with hypopituitarism possibly due to increased activation of SOX3 target genes and disregulation of β-catenin target genes. In addition, we have expanded the phenotypic spectrum associated with SOX3 mutations.

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/Sox2^0bp) or 3 base pairs (Oct4/Sox2^3bp) 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/Sox2^0bp and Oct4/Sox2^3bp 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.

Molecular Cloning, Promoter Analysis and Expression Profiles of the sox3 Gene in Japanese Flounder, Paralichthys olivaceus

Sox3, which belongs to the SoxB1 subgroup, plays major roles in neural and gonadal development. In the present study, Japanese flounder Paralichthys olivaceus sox3 gene (Posox3) and its promoter sequence were isolated and characterized. The deduced PoSox3 protein contained 298 amino acids with a characteristic HMG-box domain. Alignment and phylogenetic analyses indicated that PoSox3 shares highly identical sequence with Sox3 homologues from different species. The promoter region of Posox3 has many potential transcription factor (TF) binding sites. The expression profiles of Posox3 in different developmental stages and diverse adult tissues were analyzed by quantitative real-time RT-PCR (qRT-PCR). Posox3 mRNA was maternally inherited, and maintained at a considerably high expression level between the blastula stage and the hatching stage during embryonic development. Posox3 was abundantly expressed in the adult brain and showed sexually dimorphic expression pattern. In situ hybridization (ISH) was carried out to investigate the cellular distribution of Posox3 in the ovary, and results showed the uniform distribution of Posox3 throughout the cytoplasm of oogonia and stage I–III oocytes. These results indicate that Posox3 has potentially vital roles in embryonic and neural development and may be involved in the oogenesis process. Our work provides a fundamental understanding of the structure and potential functions of Sox3 in Paralichthys olivaceus.

Novel Function of Lysine Methyltransferase G9a in the Regulation of Sox2 Protein Stability

G9a is a lysine methyltransferase (KMTase) for histone H3 lysine 9 that plays critical roles in a number of biological processes. Emerging evidence suggests that aberrant expression of G9a contributes to tumor metastasis and maintenance of a malignant phenotype in cancer by inducing epigenetic silencing of tumor suppressor genes. Here, we show that G9a regulates Sox2 protein stability in breast cancer cells. When G9a lysine methyltransferase activity was chemically inhibited in the ER(+) breast cancer cell line MCF7, Sox2 protein levels were decreased. In addition, ectopic overexpression of G9a induced accumulation of Sox2. Changes in cell migration, invasion, and mammosphere formation by MCF7 cells were correlated with the activity or expression level of G9a. Ectopic expression of G9a also increased Sox2 protein levels in another ER(+) breast cancer cell line, ZR-75-1, whereas it did not affect Sox2 expression in MDA-MB-231 cells, an ER(-) breast cancer cell line, or in glioblastoma cell lines. Furthermore, treatment of mouse embryonic stem cells with a KMT inhibitor, BIX-01294, resulted in a rapid reduction in Sox2 protein expression despite increased Sox2 transcript levels. This finding suggests that G9a has a novel function in the regulation of Sox2 protein stability in a cell type-dependent manner.

Cooperative DNA Recognition Modulated by an Interplay between Protein-Protein Interactions and DNA-Mediated Allostery

Highly specific transcriptional regulation depends on the cooperative association of transcription factors into enhanceosomes. Usually, their DNA-binding cooperativity originates from either direct interactions or DNA-mediated allostery. Here, we performed unbiased molecular simulations followed by simulations of protein-DNA unbinding and free energy profiling to study the cooperative DNA recognition by OCT4 and SOX2, key components of enhanceosomes in pluripotent cells. We found that SOX2 influences the orientation and dynamics of the DNA-bound configuration of OCT4. In addition SOX2 modifies the unbinding free energy profiles of both DNA-binding domains of OCT4, the POU specific and POU homeodomain, despite interacting directly only with the first. Thus, we demonstrate that the OCT4-SOX2 cooperativity is modulated by an interplay between protein-protein interactions and DNA-mediated allostery. Further, we estimated the change in OCT4-DNA binding free energy due to the cooperativity with SOX2, observed a good agreement with experimental measurements, and found that SOX2 affects the relative DNA-binding strength of the two OCT4 domains. Based on these findings, we propose that available interaction partners in different biological contexts modulate the DNA exploration routes of multi-domain transcription factors such as OCT4. We consider the OCT4-SOX2 cooperativity as a paradigm of how specificity of transcriptional regulation is achieved through concerted modulation of protein-DNA recognition by different types of interactions.

Pluripotent stem cells can give rise to all somatic lineages. When taken out of the context of the embryo they can be maintained and for this a core transcriptional regulatory circuitry is crucial. OCT4 and SOX2, two factors of this network, are also critical for the induction of pluripotency in somatic cells. In pluripotent cells, OCT4 and SOX2 associate on DNA regulatory regions, enhancing or modifying each other's sequence specificity. In contrast, in the early stages during induction of pluripotency, it was proposed that OCT4 explores the genome independent of SOX2. Here we report the mechanism by which SOX2 influences the orientation, dynamics, and unbinding free energy profile of OCT4. This involves an interplay of protein-protein interactions and DNA-mediated allostery. We consider that this mechanism enables OCT4 to use its DNA binding domains and the interaction partners available in a certain biological context to access alternative genome exploration routes. This study enhances the understanding of the context specific function of OCT4 and provides a general perspective on how DNA-binding cooperativity is modulated by different types of interactions.

Differential Nucleosome Occupancies across Oct4-Sox2 Binding Sites in Murine Embryonic Stem Cells

The binding sequence for any transcription factor can be found millions of times within a genome, yet only a small fraction of these sequences encode functional transcription factor binding sites. One of the reasons for this dichotomy is that many other factors, such as nucleosomes, compete for binding. To study how the competition between nucleosomes and transcription factors helps determine a functional transcription factor site from a predicted transcription factor site, we compared experimentally-generated in vitro nucleosome occupancy with in vivo nucleosome occupancy and transcription factor binding in murine embryonic stem cells. Using a solution hybridization enrichment technique, we generated a high-resolution nucleosome map from targeted regions of the genome containing predicted sites and functional sites of Oct4/Sox2 regulation. We found that at Pax6 and Nes, which are bivalently poised in stem cells, functional Oct4 and Sox2 sites show high amounts of in vivo nucleosome displacement compared to in vitro. Oct4 and Sox2, which are active, show no significant displacement of in vivo nucleosomes at functional sites, similar to nonfunctional Oct4/Sox2 binding. This study highlights a complex interplay between Oct4 and Sox2 transcription factors and nucleosomes among different target genes, which may result in distinct patterns of stem cell gene regulation.

Identification of highly conserved putative developmental enhancers bound by SOX3 in neural progenitors using ChIP-Seq

The transcription factor SOX3 is expressed within most neural progenitor (NP) cells of the vertebrate central nervous system (CNS) and is essential for normal brain development in mice and humans. However, despite the widespread expression of Sox3, CNS defects in null mice are relatively mild due to functional redundancy with the other SOXB1 sub-group members Sox1 and Sox2. To further understand the molecular function of SOX3, we investigated the genome-wide binding profile of endogenous SOX3 in NP cells using ChIP-seq. SOX3 binding was identified at over 8,000 sites, most of which were intronic or intergeneic and were significantly associated with neurodevelopmental genes. The majority of binding sites were moderately or highly conserved (phastCons scores >0.1 and 0.5, respectively) and included the previously characterised, SOXB1-binding Nestin NP cell enhancer. Comparison of SOX3 and published ChIP-Seq data for the co-activator P300 in embryonic brain identified hundreds of highly conserved putative enhancer elements. In addition, we identified a subset of highly conserved putative enhancers for CNS development genes common to SOXB1 members in NP cells, all of which contained the SOX consensus motif (ACAAWR). Together these data implicate SOX3 in the direct regulation of hundreds of NP genes and provide molecular insight into the overlapping roles of SOXB1 proteins in CNS development.

Single molecule visualization and characterization of Sox2-Pax6 complex formation on a regulatory DNA element using a DNA origami frame

We report the use of atomic force microscopy (AFM) to study Sox2-Pax6 complex formation on the regulatory DNA element at a single molecule level. Using an origami DNA scaffold containing two DNA strands with different levels of tensile force, we confirmed that DNA bending is necessary for Sox2 binding. We also demonstrated that two transcription factors bind cooperatively by observing the increased occupancy of Sox2-Pax6 on the DNA element compared to that of Sox2 alone.

Nano-analysis of DNA conformation changes induced by transcription factor complex binding using plasmonic nanodimers

The plasmon resonant wavelength for a pair of gold nanoparticles, or gold nanodimer, increases inversely with the gap distance between the two nanoparticles. Taking advantage of this property, we performed nanoscale measurements of DNA conformation changes induced by transcription factor binding. Gold nanoparticles were bridged by double-stranded DC5 DNA that included binding sequences for the transcription factors SOX2 and PAX6, which interact on the DC5 enhancer sequence and activate transcription. The gold nanodimers bound by SOX2 shifted the plasmon resonant wavelength from 586.8 to 604.1 nm, indicating that SOX2 binding induces DNA bending. When the SOX2 formed a ternary complex with PAX6 on DC5, the plasmon resonant wavelength showed a further shift to 611.6 nm, indicating additional bending in the DC5 sequence. Furthermore, we investigated DNA conformation changes induced by SOX2 and PAX6 on the DC5-con sequence, which is a consensus sequence of DC5 for the PAX6 binding region that strengthens the PAX6 binding but at the same time disrupts SOX2-PAX6-dependent transcriptional activation. When the PAX6 binding sequence in DC5 was altered to DC5-con, the plasmon resonant wavelength shifted much less to 606.5 nm, which is more comparable to the 603.9 nm by SOX2 alone. These results demonstrate that SOX2-PAX6 cobinding induces a large conformation change in DC5 DNA.

Complex architecture and regulated expression of the Sox2ot locus during vertebrate development

The Sox2 gene is a key regulator of pluripotency embedded within an intron of a long noncoding RNA (ncRNA), termed Sox2 overlapping transcript (Sox2ot), which is transcribed in the same orientation. However, this ncRNA remains uncharacterized. Here we show that Sox2ot has multiple transcription start sites associated with genomic features that indicate regulated expression, including highly conserved elements (HCEs) and chromatin marks characteristic of gene promoters. To identify biological processes in which Sox2ot may be involved, we analyzed its expression in several developmental systems, compared to expression of Sox2. We show that Sox2ot is a stable transcript expressed in mouse embryonic stem cells, which, like Sox2, is down-regulated upon induction of embryoid body (EB) differentiation. However, in contrast to Sox2, Sox2ot is up-regulated during EB mesoderm-lineage differentiation. In adult mouse, Sox2ot isoforms were detected in tissues where Sox2 is expressed, as well as in different tissues, supporting independent regulation of expression of the ncRNA. Sox2dot, an isoform of Sox2ot transcribed from a distal HCE located >500 kb upstream of Sox2, was detected exclusively in the mouse brain, with enrichment in regions of adult neurogenesis. In addition, Sox2ot isoforms are transcribed from HCEs upstream of Sox2 in other vertebrates, including in several regions of the human brain. We also show that Sox2ot is dynamically regulated during chicken and zebrafish embryogenesis, consistently associated with central nervous system structures. These observations provide insight into the structure and regulation of the Sox2ot gene, and suggest conserved roles for Sox2ot orthologs during vertebrate development.

[Construction and expression of vector encoding Sox2 with mutated SUMO accepter site]

AIM

To construct eukaryotic expression plasmids encoding Sox2 and Sox2 K247R and identify their expression and SUMOylation.

METHODS

With gift plasmid encoding Sox2 gene as a template, Sox2 K247R was obtained by overlapping extension PCR, followed by construction of pCMV-HA-Sox2 and pCMV-HA-Sox2 K247R. After enzyme digestion analysis and DNA sequencing, these two constructs were transfected or co-transfected with pCMV-Myc-SUMO1 into 293FT cells by lipofectin method. Western blot was employed to analyze expression and SUMO ylation of Sox2.

RESULTS

It was revealed that eukaryotic expression vectors were constructed with correct sequence, where in mutant Sox2, the AAG codon was switched to CGG codon. Western blot results showed that good expression of both wt and mut Sox2, of which the latter could not be modified by SUMO1.

CONCLUSION

Successful construction and expression of Sox2 and Sox2 K247R. Sox2 could be SUMO lyated in vitro but Sox2 K247R not.

Regulated fluctuations in nanog expression mediate cell fate decisions in embryonic stem cells

The notion that the differentiated state of a cell population is determined simply by expression of specific marker genes is changing. In this work, the authors reveal that a pluripotent cell population comprises cells with temporal fluctuations in the expression of Nanog.

There is evidence that pluripotency of mouse embryonic stem (ES) cells is associated with the activity of a network of transcription factors with Sox2, Oct4, and Nanog at the core. Using fluorescent reporters for the expression of Nanog, we observed that a population of ES cells is best described by a dynamic distribution of Nanog expression characterized by two peaks defined by high (HN) and low (LN) Nanog expression. Typically, the LN state is 5%–20% of the total population, depending on the culture conditions. Modelling of the activity of Nanog reveals that a simple network of Oct4/Sox2 and Nanog activity can account for the observed distribution and its properties as long as the transcriptional activity is tuned by transcriptional noise. The model also predicts that the LN state is unstable, something that is born out experimentally. While in this state, cells can differentiate. We suggest that transcriptional fluctuations in Nanog expression are an essential element of the pluripotent state and that the function of Sox2, Oct4, and Nanog is to act as a network that promotes and maintains transcriptional noise to interfere with the differentiation signals.

Embryonic stem (ES) cells are a pluripotent cell population derived from early mammalian embryos. An intrinsic feature of ES cells is their phenotypic heterogeneity: they display promiscuous activation of lineage-specific genes and exhibit a fluctuating flow of differentiating cells. A gene regulatory network (GRN) centred around the transcription factors Sox2, Oct4, and Nanog is essential for the establishment and the maintenance of the pluripotent state. Previous studies had suggested that ES cells can reversibly change their state of Nanog expression without losing pluripotency. Here, we extend these studies by quantifying and monitoring the expression of Nanog in a Nanog-GFP reporter cell line. We show that Nanog levels undergo slow, random fluctuations in ES cells, giving rise to heterogeneous cell populations. We identify two states, one stable, characterized by high levels of expression (HN), and another with low levels of Nanog expression (LN), which is highly unstable. While in the LN state, cells are more likely to differentiate depending on the culture medium. Mathematical modelling shows that a simple excitable system driven by transcriptional noise can account for the observed distributions and behaviours in gene expression. Our study suggests that rather than a discrete state dependent on the fixed expression of a small set of genes, pluripotency is best represented by a state of dynamic heterogeneity of a population driven by transcriptional noise, and that the function of the gene regulatory network centred around Nanog might be to generate dynamic heterogeneities at the population level.