Co-motif discovery identifies an Esrrb-Sox2-DNA ternary complex as a mediator of transcriptional differences between mouse embryonic and epiblast stem cells.

Transcription factors (TF) often bind in heterodimeric complexes with each TF recognizing a specific neighboring cis element in the regulatory region of the genome. Comprehension of this DNA motif grammar is opaque, yet recent developments have allowed the interrogation of genome-wide TF binding sites. We reasoned that within this data novel motif grammars could be identified that controlled distinct biological programs. For this purpose, we developed a novel motif-discovery tool termed fexcom that systematically interrogates ChIP-seq data to discover spatially constrained TF-TF composite motifs occurring over short DNA distances. We applied this to the extensive ChIP-seq data available from mouse embryonic stem cells (ESCs). In addition to the well-known and most prevalent sox-oct motif, we also discovered a novel constrained spacer motif for Esrrb and Sox2 with a gap of between 2 and 8 bps that Essrb and Sox2 cobind in a selective fashion. Through the use of knockdown experiments, we argue that the Esrrb-Sox2 complex is an arbiter of gene expression differences between ESCs and epiblast stem cells (EpiSC). A number of genes downregulated upon dual Esrrb/Sox2 knockdown (e.g., Klf4, Klf5, Jam2, Pecam1) are similarly downregulated in the ESC to EpiSC transition and contain the esrrb-sox motif. The prototypical Esrrb-Sox2 target gene, containing an esrrb-sox element conserved throughout eutherian and metatherian mammals, is Nr0b1. Through positive regulation of this transcriptional repressor, we argue the Esrrb-Sox2 complex promotes the ESC state through inhibition of the EpiSC transcriptional program and the same trio may also function to maintain trophoblast stem cells.

Sox transcription factors require selective interactions with Oct4 and specific transactivation functions to mediate reprogramming.

The unique ability of Sox2 to cooperate with Oct4 at selective binding sites in the genome is critical for reprogramming somatic cells into induced pluripotent stem cells (iPSCs). We have recently demonstrated that Sox17 can be converted into a reprogramming factor by alteration of a single amino acid (Sox17EK) within its DNA binding HMG domain. Here we expanded this study by introducing analogous mutations to 10 other Sox proteins and interrogated the role of N-and C-termini on the reprogramming efficiency. We found that point-mutated Sox7 and Sox17 can convert human and mouse fibroblasts into iPSCs, but Sox4, Sox5, Sox6, Sox8, Sox9, Sox11, Sox12, Sox13, and Sox18 cannot. Next we studied regions outside the HMG domain and found that the C-terminal transactivation domain of Sox17 and Sox7 enhances the potency of Sox2 in iPSC assays and confers weak reprogramming potential to the otherwise inactive Sox4EK and Sox18EK proteins. These results suggest that the glutamate (E) to lysine (K) mutation in the HMG domain is necessary but insufficient to swap the function of Sox factors. Moreover, the HMG domain alone fused to the VP16 transactivation domain is able to induce reprogramming, albeit at low efficiency. By molecular dissection of the C-terminus of Sox17, we found that the ß-catenin interaction region contributes to the enhanced reprogramming efficiency of Sox17EK. To mechanistically understand the enhanced reprogramming potential of Sox17EK, we analyzed ChIP-sequencing and expression data and identified a subset of candidate genes specifically regulated by Sox17EK and not by Sox2.

Conversion of Sox17 into a pluripotency reprogramming factor by reengineering its association with Oct4 on DNA.

Very few proteins are capable to induce pluripotent stem (iPS) cells and their biochemical uniqueness remains unexplained. For example, Sox2 cooperates with other transcription factors to generate iPS cells, but Sox17, despite binding to similar DNA sequences, cannot. Here, we show that Sox2 and Sox17 exhibit inverse heterodimerization preferences with Oct4 on the canonical versus a newly identified compressed sox/oct motif. We can swap the cooperativity profiles of Sox2 and Sox17 by exchanging single amino acids at the Oct4 interaction interface resulting in Sox2KE and Sox17EK proteins. The reengineered Sox17EK now promotes reprogramming of somatic cells to iPS, whereas Sox2KE has lost this potential. Consistently, when Sox2KE is overexpressed in embryonic stem cells it forces endoderm differentiation similar to wild-type Sox17. Together, we demonstrate that strategic point mutations that facilitate Sox/Oct4 dimer formation on variant DNA motifs lead to a dramatic swap of the bioactivities of Sox2 and Sox17.

Evolution of the vertebrate gene regulatory network controlled by the transcriptional repressor REST.

Specific wiring of gene-regulatory networks is likely to underlie much of the phenotypic difference between species, but the extent of lineage-specific regulatory architecture remains poorly understood. The essential vertebrate transcriptional repressor REST (RE1-Silencing Transcription Factor) targets many neural genes during development of the preimplantation embryo and the central nervous system, through its cognate DNA motif, the RE1 (Repressor Element 1). Here we present a comparative genomic analysis of REST recruitment in multiple species by integrating both sequence and experimental data. We use an accurate, experimentally validated Position-Specific Scoring Matrix method to identify REST binding sites in multiply aligned vertebrate genomes, allowing us to infer the evolutionary origin of each of 1,298 human RE1 elements. We validate these findings using experimental data of REST binding across the whole genomes of human and mouse. We show that one-third of human RE1s are unique to primates: These sites recruit REST in vivo, target neural genes, and are under purifying evolutionary selection. We observe a consistent and significant trend for more ancient RE1s to have higher affinity for REST than lineage-specific sites and to be more proximal to target genes. Our results lead us to propose a model where new transcription factor binding sites are constantly generated throughout the genome; thereafter, refinement of their sequence and location consolidates this remodeling of networks governing neural gene regulation.

Crystal optimization and preliminary diffraction data analysis of the Smad1 MH1 domain bound to a palindromic SBE DNA element.

The bone morphogenetic protein (BMP) signalling pathway regulates diverse processes such as cell differentiation, anterior/posterior axis specification, cell growth and the formation of extra-embryonic tissues. The transcription factor Smad1 relays the BMP signal from the cytoplasm to the nucleus, where it binds short DNA-sequence motifs and regulates gene expression. However, how Smad1 selectively targets particular genomic regions is poorly understood. In order to understand the physical basis of the specific interaction of Smad1 with DNA and to contrast it with the highly homologous but functionally distinct Smad3 protein, the DNA-binding Mad-homology 1 (MH1) domain of Smad1 was cocrystallized with a 17-mer palindromic Smad-binding element (SBE). The extensive optimizations of the length, binding-site spacing and terminal sequences of the DNA element in combination with the other crystallization parameters necessary for obtaining diffraction-quality crystals are described here. A 2.7 angstrom resolution native data set was collected at the National Synchrotron Radiation Research Centre, Taiwan, from crystals grown in a solution containing 0.2 M ammonium tartrate dibasic, 20% PEG 3350, 3% 2-propanol and 10% glycerol. The data set was indexed and merged in space group P222, with unit-cell parameters a = 73.94, b = 77.49, c = 83.78 angstrom, alpha = beta = gamma = 90 degrees. The solvent content in the unit cell is consistent with the presence of two Smad1 MH1 molecules bound to the duplex DNA in the asymmetric unit.

Purification, crystallization and preliminary X-ray diffraction analysis of the HMG domain of Sox17 in complex with DNA.

Sox17 is a member of the SRY-related high-mobility group (HMG) of transcription factors that have been shown to direct endodermal differentiation in early mammalian development. The LAMA1 gene encoding the alpha-chain of laminin-1 has been reported to be directly bound and regulated by Sox17. This paper describes the details of initial crystallization attempts with the HMG domain of mouse Sox17 (mSox17-HMG) with a 16-mer DNA element derived from the LAMA1 enhancer and optimization strategies to obtain a better diffracting crystal. The best diffracting crystal was obtained in a condition containing 0.1 M Tris-HCl pH 7.4, 0.2 M MgCl(2), 30% PEG 3350 using the hanging-drop vapour-diffusion method. A highly redundant in-house data set was collected to 2.75 A resolution with 99% completeness. The presence of the mSox17-HMG-DNA complex within the crystals was confirmed and Matthews analysis indicated the presence of one complex per asymmetric unit.

Structural basis for the SOX-dependent genomic redistribution of OCT4 in stem cell differentiation.

In pluripotent cells, OCT4 associates with SOX2 to maintain pluripotency or with SOX17 to induce primitive endoderm commitment. The OCT4-SOX2 and OCT4-SOX17 combinations bind mutually exclusive to two distinct composite DNA elements, known as the "canonical" and "compressed" motifs, respectively. The structural basis for the OCT4-SOX17 cooperativity is unknown. Whereas SOX17 has been engineered to replace SOX2 in the pluripotency circuitry, all generated SOX2 mutants have failed to act like SOX17. From molecular simulations, we revealed the OCT4-SOX17 interaction interface and elucidated the SOX-dependent motif preference of OCT4. Moreover, we designed a SOX2 mutant that we predicted and confirmed experimentally to bind cooperatively with OCT4 to the compressed motif. Ultimately, we found a strong correlation between the experimental and calculated relative cooperative-binding free energies of 12 OCT4-SOX-DNA complexes. Therefore, we validated the OCT4-SOX interfaces and demonstrated that in silico design of DNA-binding cooperativity is suitable for altering transcriptional circuitries.

The structure of Sox17 bound to DNA reveals a conserved bending topology but selective protein interaction platforms.

Sox17 regulates endodermal lineage commitment and is thought to function antagonistically to the pluripotency determinant Sox2. To investigate the biochemical basis for the distinct functions of Sox2 and Sox17, we solved the crystal structure of the high mobility group domain of Sox17 bound to a DNA element derived from the Lama1 enhancer using crystals diffracting to 2.7 A resolution. Sox17 targets the minor groove and bends the DNA by approximately 80 degrees . The DNA architecture closely resembles the one seen for Sox2/DNA structures, suggesting that the degree of bending is conserved between both proteins and nucleotide substitutions have only marginal effects on the bending topology. Accordingly, affinities of Sox2 and Sox17 for the Lama1 element were found to be identical. However, when the Oct1 contact interface of Sox2 is compared with the corresponding region of Sox17, a significantly altered charge distribution is observed, suggesting differential co-factor recruitment that may explain their biological distinctiveness.

Crystal structure and DNA binding of the homeodomain of the stem cell transcription factor Nanog.

The transcription factor Nanog is an upstream regulator in early mammalian development and a key determinant of pluripotency in embryonic stem cells. Nanog binds to promoter elements of hundreds of target genes and regulates their expression by an as yet unknown mechanism. Here, we report the crystal structure of the murine Nanog homeodomain (HD) and analysis of its interaction with a DNA element derived from the Tcf3 promoter. Two Nanog amino acid pairs, unique among HD sequences, appear to affect the mechanism of nonspecific DNA recognition as well as maintain the integrity of the structural scaffold. To assess selective DNA recognition by Nanog, we performed electrophoretic mobility shift assays using a panel of modified DNA binding sites and found that Nanog HD preferentially binds the TAAT(G/T)(G/T) motif. A series of rational mutagenesis experiments probing the role of six variant residues of Nanog on its DNA binding function establish their role in affecting binding affinity but not binding specificity. Together, the structural and functional evidence establish Nanog as a distant member of a Q50-type HD despite having considerable variation at the sequence level.