Nucleosomal DNA Dynamics Mediate Oct4 Pioneer Factor Binding

Conformational Dynamics of the Nucleosomal Histone H2B Tails Revealed by Molecular Dynamics Simulations

Epigenetic modifications of histone N-terminal tails play a critical role in regulating the chromatin structure and biological processes such as transcription and DNA repair. One of the key post-translational modifications (PTMs) is the acetylation of lysine residues on histone tails. Epigenetic modifications are ubiquitous in the development of diseases, such as cancer and neurological disorders. Histone H2B tails are critical regulators of nucleosome dynamics, biological processes, and certain diseases. Here, we report all-atomistic molecular dynamics (MD) simulations of the nucleosome to demonstrate that acetylation of the histone tails changes their conformational space and interaction with DNA. We perform simulations of H2B tails, critical regulators of gene regulation, in both the lysine-acetylated (ACK) and unacetylated wild type (WT) states. To explore the effects of salt concentration, we use two different NaCl concentrations to perform simulations at microsecond time scales. Salt can modulate the effects of electrostatic interactions between the DNA phosphate backbone and histone tails. Upon acetylation, H2B tails shift their secondary structure helical propensity. The number of contacts between the DNA and the H2B tail decreases. We characterize the conformational dynamics of the H2B tails by principal component analysis (PCA). The ACK tails become more compact at increased salt concentrations, but conformations from the WT tails display the most contacts with DNA at both salt concentrations. Mainly, H2B acetylation may increase the DNA accessibility for regulatory proteins to bind, which can aid in gene regulation and NCP stability.

Pioneer factors: nature or nurture?

Chromatin is densely packed with nucleosomes, which limits the accessibility of many chromatin-associated proteins. Pioneer factors (PFs) are usually viewed as a special group of sequence-specific transcription factors (TFs) that can recognize nucleosome-embedded motifs, invade compact chromatin, and generate open chromatin regions. Through this process, PFs initiate a cascade of events that play key roles in gene regulation and cell differentiation. A current debate in the field is if PFs belong to a unique subset of TFs with intrinsic "pioneering activity", or if all TFs have the potential to function as PFs within certain cellular contexts. There are also different views regarding the key feature(s) that define pioneering activity. In this review, we present evidence from the literature related to these alternative views and discuss how to potentially reconcile them. It is possible that both intrinsic properties, like tight nucleosome binding and structural compatibility, and cellular conditions, like concentration and co-factor availability, are important for PF function.

Regulation of chromatin architecture by protein binding: insights from molecular modeling

Histone and non-histone proteins play key roles in the activation and repression of genes. In addition to experimental studies of their regulation of gene expression, molecular modeling at the nucleosome, chromatin, and chromosome levels can contribute insights into the molecular mechanisms involved. In this review, we provide an overview for protein-bound chromatin modeling, and describe how our group has integrated protein binding into genome systems across the scales, from all-atom to coarse-grained models, using explicit to implicit descriptions. We describe the associated applications to protein binding effects and biological mechanisms of genome folding and gene regulation. We end by illustrating the application of machine learning tools like AlphaFold2 to proteins relevant to chromatin systems.

Pioneer factors: roles and their regulation in development

Pioneer factors are a subclass of transcription factors that can bind and initiate opening of silent chromatin regions. Pioneer factors subsequently regulate lineage-specific genes and enhancers and, thus, activate the zygotic genome after fertilization, guide cell fate transitions during development, and promote various forms of human cancers. As such, pioneer factors are useful in directed cell reprogramming. In this review, we define the structural and functional characteristics of pioneer factors, how they bind and initiate opening of closed chromatin regions, and the consequences for chromatin dynamics and gene expression during cell differentiation. We also discuss emerging mechanisms that modulate pioneer factors during development.

Regulation of chromatin architecture by transcription factor binding

Transcription factors (TF) bind to chromatin and regulate the expression of genes. The pair Myc:Max binds to E-box regulatory DNA elements throughout the genome to control the transcription of a large group of specific genes. We introduce an implicit modeling protocol for Myc:Max binding to mesoscale chromatin fibers at nucleosome resolution to determine TF effect on chromatin architecture and shed light into its mechanism of gene regulation. We first bind Myc:Max to different chromatin locations and show how it can direct fiber folding and formation of microdomains, and how this depends on the linker DNA length. Second, by simulating increasing concentrations of Myc:Max binding to fibers that differ in the DNA linker length, linker histone density, and acetylation levels, we assess the interplay between Myc:Max and other chromatin internal parameters. Third, we study the mechanism of gene silencing by Myc:Max binding to the Eed gene loci. Overall, our results show how chromatin architecture can be regulated by TF binding. The position of TF binding dictates the formation of microdomains that appear visible only at the ensemble level. At the same time, the level of linker histone and tail acetylation, or different linker DNA lengths, regulates the concentration-dependent effect of TF binding. Furthermore, we show how TF binding can repress gene expression by increasing fiber folding motifs that help compact and occlude the promoter region. Importantly, this effect can be reversed by increasing linker histone density. Overall, these results shed light on the epigenetic control of the genome dictated by TF binding.

Highly cooperative chimeric super-SOX induces naive pluripotency across species

Our understanding of pluripotency remains limited: iPSC generation has only been established for a few model species, pluripotent stem cell lines exhibit inconsistent developmental potential, and germline transmission has only been demonstrated for mice and rats. By swapping structural elements between Sox2 and Sox17, we built a chimeric super-SOX factor, Sox2-17, that enhanced iPSC generation in five tested species: mouse, human, cynomolgus monkey, cow, and pig. A swap of alanine to valine at the interface between Sox2 and Oct4 delivered a gain of function by stabilizing Sox2/Oct4 dimerization on DNA, enabling generation of high-quality OSKM iPSCs capable of supporting the development of healthy all-iPSC mice. Sox2/Oct4 dimerization emerged as the core driver of naive pluripotency with its levels diminished upon priming. Transient overexpression of the SK cocktail (Sox+Klf4) restored the dimerization and boosted the developmental potential of pluripotent stem cells across species, providing a universal method for naive reset in mammals.

Regulation of Chromatin Architecture by Transcription Factor Binding

Transcription factors (TF) bind to chromatin and regulate the expression of genes. The pair Myc:Max binds to E-box regulatory DNA elements throughout the genome, controlling transcription of a large group of specific genes. We introduce an implicit modeling protocol for Myc:Max binding to mesoscale chromatin fibers to determine TF effect on chromatin architecture and shed light on its mechanism of gene regulation. We first bind Myc:Max to different chromatin locations and show how it can direct fiber folding and formation of microdomains, and how this depends on the linker DNA length. Second, by simulating increasing concentrations of Myc:Max binding to fibers that differ in the DNA linker length, linker histone density, and acetylation levels, we assess the interplay between Myc:Max and other chromatin internal parameters. Third, we study the mechanism of gene silencing by Myc:Max binding to the Eed gene loci. Overall, our results show how chromatin architecture can be regulated by TF binding. The position of TF binding dictates the formation of microdomains that appear visible only at the ensemble level. On the other hand, the presence of linker histone, acetylations, or different linker DNA lengths regulates the concentration-dependent effect of TF binding. Furthermore, we show how TF binding can repress gene expression by increasing fiber folding motifs that help compact and occlude the promoter region. Importantly, this effect can be reversed by increasing linker histone density. Overall, these results shed light on the epigenetic control of the genome dictated by TF binding.

Structural and Dynamic Changes of Nucleosome upon GATA3 Binding

Pioneer factors, which can directly bind to nucleosomes, have been considered to change chromatin conformations. However, the binding impact on the nucleosome is little known. Here, we show how the pioneer factor GATA3 binds to nucleosomal DNA and affects the conformation and dynamics of nucleosomes by using a combination of SAXS, molecular modeling, and molecular dynamics simulations. Our structural models, consistent with the SAXS data, indicate that only one of the two DNA binding domains, N- and C-fingers, of GATA3 binds to an end of the DNA in solution. Our MD simulations further showed that the other unbound end of the DNA increases the fluctuation and enhances the DNA dissociation from the histone core when the N-finger binds to a DNA end, a site near the entry or exit of the nucleosome. However, this was not true for the binding of the C-finger that binds to a location about 15 base pairs distant from the DNA end. In this case, DNA dissociation occurred on the bound end. Taken together, we suggest that the N-finger and C-finger bindings of GATA3 commonly enhance DNA dissociation at one of the two DNA ends (the bound end for the C-finger binding and the unbound end for the N-finger binding), leading to triggering a conformational change in the chromatin.

Structural mechanism of LIN28B nucleosome targeting by OCT4

Pioneer transcription factors are essential for cell fate changes by targeting closed chromatin. OCT4 is a crucial pioneer factor that can induce cell reprogramming. However, the structural basis of how pioneer factors recognize the in vivo nucleosomal DNA targets is unknown. Here, we determine the high-resolution structures of the nucleosome containing human LIN28B DNA and its complexes with the OCT4 DNA binding region. Three OCT4s bind the pre-positioned nucleosome by recognizing non-canonical DNA sequences. Two use their POUS domains while the other uses the POUS-loop-POUHD region; POUHD serves as a wedge to unwrap ∼25 base pair DNA. Our analysis of previous genomic data and determination of the ESRRB-nucleosome-OCT4 structure confirmed the generality of these structural features. Moreover, biochemical studies suggest that multiple OCT4s cooperatively open the H1-condensed nucleosome array containing the LIN28B nucleosome. Thus, our study suggests a mechanism of how OCT4 can target the nucleosome and open closed chromatin.

Molecular Mechanism of Nucleosome Recognition by the Pioneer Transcription Factor Sox

Pioneer transcription factors (PTFs) have the remarkable ability to directly bind to chromatin to stimulate vital cellular processes. In this work, we dissect the universal binding mode of Sox PTF by combining extensive molecular simulations and physiochemistry approaches, along with DNA footprinting techniques. As a result, we show that when Sox consensus DNA is located at the solvent-facing DNA strand, Sox binds to the compact nucleosome without imposing any significant conformational changes. We also reveal that the base-specific Sox:DNA interactions (base reading) and Sox-induced DNA changes (shape reading) are concurrently required for sequence-specific nucleosomal DNA recognition. Among three different nucleosome positions located on the positive DNA arm, a sequence-specific reading mechanism is solely satisfied at the superhelical location 2 (SHL2). While SHL2 acts transparently for solvent-facing Sox binding, among the other two positions, SHL4 permits only shape reading. The final position, SHL0 (dyad), on the other hand, allows no reading mechanism. These findings demonstrate that Sox-based nucleosome recognition is essentially guided by intrinsic nucleosome properties, permitting varying degrees of DNA recognition.

Machines on Genes through the Computational Microscope

Macromolecular machines acting on genes are at the core of life's fundamental processes, including DNA replication and repair, gene transcription and regulation, chromatin packaging, RNA splicing, and genome editing. Here, we report the increasing role of computational biophysics in characterizing the mechanisms of "machines on genes", focusing on innovative applications of computational methods and their integration with structural and biophysical experiments. We showcase how state-of-the-art computational methods, including classical and ab initio molecular dynamics to enhanced sampling techniques, and coarse-grained approaches are used for understanding and exploring gene machines for real-world applications. As this review unfolds, advanced computational methods describe the biophysical function that is unseen through experimental techniques, accomplishing the power of the "computational microscope", an expression coined by Klaus Schulten to highlight the extraordinary capability of computer simulations. Pushing the frontiers of computational biophysics toward a pragmatic representation of large multimegadalton biomolecular complexes is instrumental in bridging the gap between experimentally obtained macroscopic observables and the molecular principles playing at the microscopic level. This understanding will help harness molecular machines for medical, pharmaceutical, and biotechnological purposes.

Studies of the Mechanism of Nucleosome Dynamics: A Review on Multifactorial Regulation from Computational and Experimental Cases

The nucleosome, which organizes the long coil of genomic DNA in a highly condensed, polymeric way, is thought to be the basic unit of chromosomal structure. As the most important protein–DNA complex, its structural and dynamic features have been successively revealed in recent years. However, its regulatory mechanism, which is modulated by multiple factors, still requires systemic discussion. This study summarizes the regulatory factors of the nucleosome’s dynamic features from the perspective of histone modification, DNA methylation, and the nucleosome-interacting factors (transcription factors and nucleosome-remodeling proteins and cations) and focuses on the research exploring the molecular mechanism through both computational and experimental approaches. The regulatory factors that affect the dynamic features of nucleosomes are also discussed in detail, such as unwrapping, wrapping, sliding, and stacking. Due to the complexity of the high-order topological structures of nucleosomes and the comprehensive effects of regulatory factors, the research on the functional modulation mechanism of nucleosomes has encountered great challenges. The integration of computational and experimental approaches, the construction of physical modes for nucleosomes, and the application of deep learning techniques will provide promising opportunities for further exploration.

Genome modeling: From chromatin fibers to genes

The intricacies of the 3D hierarchical organization of the genome have been approached by many creative modeling studies. The specific model/simulation technique combination defines and restricts the system and phenomena that can be investigated. We present the latest modeling developments and studies of the genome, involving models ranging from nucleosome systems and small polynucleosome arrays to chromatin fibers in the kb-range, chromosomes, and whole genomes, while emphasizing gene folding from first principles. Clever combinations allow the exploration of many interesting phenomena involved in gene regulation, such as nucleosome structure and dynamics, nucleosome-nucleosome stacking, polynucleosome array folding, protein regulation of chromatin architecture, mechanisms of gene folding, loop formation, compartmentalization, and structural transitions at the chromosome and genome levels. Gene-level modeling with full details on nucleosome positions, epigenetic factors, and protein binding, in particular, can in principle be scaled up to model chromosomes and cells to study fundamental biological regulation.

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.

The homeodomain of Oct4 is a dimeric binder of methylated CpG elements

Oct4 is essential to maintain pluripotency and has a pivotal role in establishing the germline. Its DNA-binding POU domain was recently found to bind motifs with methylated CpG elements normally associated with epigenetic silencing. However, the mode of binding and the consequences of this capability has remained unclear. Here, we show that Oct4 binds to a compact palindromic DNA element with a methylated CpG core (CpGpal) in alternative states of pluripotency and during cellular reprogramming towards induced pluripotent stem cells (iPSCs). During cellular reprogramming, typical Oct4 bound enhancers are uniformly demethylated, with the prominent exception of the CpGpal sites where DNA methylation is often maintained. We demonstrate that Oct4 cooperatively binds the CpGpal element as a homodimer, which contrasts with the ectoderm-expressed POU factor Brn2. Indeed, binding to CpGpal is Oct4-specific as other POU factors expressed in somatic cells avoid this element. Binding assays combined with structural analyses and molecular dynamic simulations show that dimeric Oct4-binding to CpGpal is driven by the POU-homeodomain whilst the POU-specific domain is detached from DNA. Collectively, we report that Oct4 exerts parts of its regulatory function in the context of methylated DNA through a DNA recognition mechanism that solely relies on its homeodomain.

Structural mechanism of LIN28B nucleosome targeting by OCT4 for pluripotency

Pioneer transcription factors are essential for cell fate changes by targeting closed chromatin. OCT4 is a crucial pioneer factor that can induce cell reprogramming. However, the structural basis of how pioneer factors recognize the in vivo nucleosomal DNA targets is unknown. Here, we determine the high-resolution structures of the nucleosome containing human LIN28B DNA and its complexes with the OCT4 DNA binding region. Three OCT4s bind the pre-positioned nucleosome by recognizing non-canonical DNA motifs. Two use their POUS domains by forming extensive hydrogen bonds. The other uses the POUS-loop-POUHD region; POUHD serves as a wedge to unwrap ∼25 base pair DNA. Biochemical studies suggest that multiple OCT4s cooperatively open the H1-condensed nucleosome array containing the LIN28B nucleosome. Our study suggests a mechanism whereby OCT4s target the LIN28B nucleosome by forming multivalent interactions with nucleosomal motifs, unwrapping nucleosomal DNA, evicting H1, and cooperatively open closed chromatin to initiate cell reprogramming.

Molecular dynamics simulations for the study of chromatin biology

The organization of Eukaryotic DNA into chromatin has profound implications for the processing of genetic information. In the past years, molecular dynamics (MD) simulations proved to be a powerful tool to investigate the mechanistic basis of chromatin biology. We review recent all-atom and coarse-grained MD studies revealing how the structure and dynamics of chromatin underlie its biological functions. We describe the latest method developments; the structural fluctuations of nucleosomes and the various factors affecting them; the organization of chromatin fibers, with particular emphasis on its liquid-like character; the interactions and dynamics of transcription factors on chromatin; and how chromatin organization is modulated by molecular motors acting on DNA.

OCT4 interprets and enhances nucleosome flexibility

Pioneer transcription factors are proteins that induce cellular identity transitions by binding to inaccessible regions of DNA in nuclear chromatin. They contribute to chromatin opening and recruit other factors to regulatory DNA elements. The structural features and dynamics modulating their interaction with nucleosomes are still unresolved. From a combination of experiments and molecular simulations, we reveal here how the pioneer factor and master regulator of pluripotency, Oct4, interprets and enhances nucleosome structural flexibility. The magnitude of Oct4’s impact on nucleosome dynamics depends on the binding site position and the mobility of the unstructured tails of nucleosomal histone proteins. Oct4 uses both its DNA binding domains to propagate and stabilize open nucleosome conformations, one for specific sequence recognition and the other for nonspecific interactions with nearby regions of DNA. Our findings provide a structural basis for the versatility of transcription factors in engaging with nucleosomes and have implications for understanding how pioneer factors induce chromatin dynamics.

Insights into the multi-faceted role of Pioneer transcription factors in glioma formation and progression with targeting options

Pioneer transcription factors (TFs) present an important subtype of transcription factors which are vital for cell programming during embryonic development and cellular memory during mitotic growth, as well as cell fate reprogramming. Pioneer TFs can engage specific target binding sites on nucleosomal DNA to attract chromatin remodeling complexes, cofactors, and other transcription factors, ultimately controlling gene expression by shaping locally the epigenome. The priority of binding that they exhibit in contrast to other transcription factors and their involvement in crucial events regarding cell fate, has implicated their aberrant function in the pathogenesis of several disorders including carcinogenesis. Emerging experimental data indicate that certain Pioneer TFs are highly implicated in gliomas development, in neoplastic cell proliferation, angiogenic processes, resistance to therapy, and patient survival. Herein, we describe the main structural characteristics and functional mechanisms of pioneer TFs, focusing on their central role in the pathogenesis and progression of gliomas. We further highlight the current treatment options ranging from natural agents (oleanolic acid) to a variety of chemical compounds (APR-246, COTI-2) and discuss potential delivery systems, including nanoparticles, viral vectors, and intracellular protein delivery techniques.

H2A-H2B Histone Dimer Plasticity and Its Functional Implications

The protein core of the nucleosome is composed of an H3-H4 histone tetramer and two H2A-H2B histone dimers. The tetramer organizes the central 60 DNA bp, while H2A-H2B dimers lock the flanking DNA segments. Being positioned at the sides of the nucleosome, H2A-H2B dimers stabilize the overall structure of the nucleosome and modulate its dynamics, such as DNA unwrapping, sliding, etc. Such modulation at the epigenetic level is achieved through post-translational modifications and the incorporation of histone variants. However, the detailed connection between the sequence of H2A-H2B histones and their structure, dynamics and implications for nucleosome functioning remains elusive. In this work, we present a detailed study of H2A-H2B dimer dynamics in the free form and in the context of nucleosomes via atomistic molecular dynamics simulations (based on X. laevis histones). We supplement simulation results by comparative analysis of information in the structural databases. Particularly, we describe a major dynamical mode corresponding to the bending movement of the longest H2A and H2B α-helices. This overall bending dynamics of the H2A-H2B dimer were found to be modulated by its interactions with DNA, H3-H4 tetramer, the presence of DNA twist-defects with nucleosomal DNA and the amino acid sequence of histones. Taken together, our results shed new light on the dynamical mechanisms of nucleosome functioning, such as nucleosome sliding, DNA-unwrapping and their epigenetic modulation.

Structures and consequences of pioneer factor binding to nucleosomes

Pioneer transcription factors are able to bind a partially exposed motif on the surface of a nucleosome, enabling the proteins to target sites in silent regions of chromatin that have been compacted by linker histone. The targeting of nucleosomal DNA by pioneer factors has been observed in vitro and in vivo, where binding can promote local nucleosome exposure that allows other transcription factors, nucleosome remodelers, and histone modifiers to engage the chromatin and elicit gene activation or further repression. Pioneer factors thereby establish new gene expression programs during cell fate changes that occur during embryonic development, regeneration, and cancer. Here, we review recent biophysical studies that reveal the structural features and strategies used by pioneer factors to accomplish nucleosome binding and the consequential changes to nucleosomes that can lead to DNA accessibility.

Integrated Chromatin Accessibility and Transcriptome Landscapes of 5-Fluorouracil-Resistant Colon Cancer Cells

Background: 5-Fluorouracil (5-FU) is one of the most effective and widely used chemotherapeutic drugs in the treatment of colon cancer, yet chemoresistance is a common feature of colon cancer treatment, resulting in poor prognosis and short survival. Dynamic reprogramming of chromatin accessibility is crucial for proper regulation of gene transcription associated with cancer drug resistance by providing the gene regulatory machinery with rapid access to the open genomic DNA.

Methods: Here, we explored the global chromatin accessibility and transcription changes by the assay for transposase-accessible chromatin using sequencing (ATAC-seq) in combination with transcriptome sequencing of both parental and 5-FU-resistant HCT15 cells, followed by integrative analysis to better understand the regulatory network underlying 5-FU resistance in colon cancer cells.

Results: A total of 3,175 differentially expressed mRNAs (DEGs), lncRNAs (DELs), and miRNAs (DEMs) related to 5-FU resistance were identified, including significantly upregulated IL33, H19, and miR-17-5p; the downregulated AKR1B10, LINC01012, and miR-125b-5p; and chromatin modifiers such as INO80C, HDAC6, and KDM5A. The construction of the ceRNA regulatory network revealed that H19, HOXA11-AS, and NEAT1 might function as ceRNAs associated with 5-FU resistance in HCT15 cells. Moreover, 9,868 differentially accessible regions (DARs) were obtained, which were positively (r = 0.58) correlated with their nearest DEGs and DELs. The upregulated genes related to 4,937 hyper-accessible regions were significantly enriched in signaling pathways of MAPK, FOX, and WNT, while the 4,931 hypo-accessible regions were considered to be involved in declined biosynthesis of amino acids and nucleotide sugars, signaling pathways of Notch, and HIF-1. Analyses of the DAR sequences revealed that besides the AP-1 family, the TF motifs of FOX and KLF family members were highly enriched in hyper- and hypo-accessible regions, respectively. Finally, we obtained several critical TFs and their potential targets associated with DARs and 5-FU resistance, including FOXA1 and KLF3.

Conclusion: These data provided clear insights and valuable resources for an improved understanding of the non-genetic landscape of 5-FU-resistant colon cancer cells based on chromatin accessibility and transcript levels, which allowed for genome-wide detection of TF binding sites, potential cis-regulatory elements and therapeutic targets.

A test of the pioneer factor hypothesis using ectopic liver gene activation

The pioneer factor hypothesis (PFH) states that pioneer factors (PFs) are a subclass of transcription factors (TFs) that bind to and open inaccessible sites and then recruit non-pioneer factors (non-PFs) that activate batteries of silent genes. The PFH predicts that ectopic gene activation requires the sequential activity of qualitatively different TFs. We tested the PFH by expressing the endodermal PF FOXA1 and non-PF HNF4A in K562 lymphoblast cells. While co-expression of FOXA1 and HNF4A activated a burst of endoderm-specific gene expression, we found no evidence for a functional distinction between these two TFs. When expressed independently, both TFs bound and opened inaccessible sites, activated endodermal genes, and ‘pioneered’ for each other, although FOXA1 required fewer copies of its motif for binding. A subset of targets required both TFs, but the predominant mode of action at these targets did not conform to the sequential activity predicted by the PFH. From these results, we hypothesize an alternative to the PFH where ‘pioneer activity’ depends not on categorically different TFs but rather on the affinity of interaction between TF and DNA.

Cells only use a fraction of their genetic information to make the proteins they need. The rest is carefully packaged away and tightly bundled in structures called nucleosomes. This physically shields the DNA from being accessed by transcription factors – the molecular actors that can read genes and kickstart the protein production process. Effectively, the genetic sequences inside nucleosomes are being silenced.

However, during development, transcription factors must overcome this nucleosome barrier and activate silent genes to program cells. The pioneer factor hypothesis describes how this may be possible: first, ‘pioneer’ transcription factors can bind to and ‘open up’ nucleosomes to make target genes accessible. Then, non-pioneer factors can access the genetic sequence and recruit cofactors that begin copying the now-exposed genetic information.

The widely accepted theory is based on studies of two proteins – FOXA1, an archetypal pioneer factor, and HNF4A, a non-pioneer factor – but the predictions of the pioneer factor hypothesis have yet to be explicitly tested. To do so, Hansen et al. expressed FOXA1 and HNF4A, separately and together, in cells which do not usually make these proteins. They then assessed how the proteins could bind to DNA and impact gene accessibility and transcription.

The experiments demonstrate that FOXA1 and HNF4A do not necessarily follow the two-step activation predicted by the pioneer factor hypothesis. When expressed independently, both transcription factors bound and opened inaccessible sites, activated target genes, and ‘pioneered’ for each other. Similar patterns were observed across the genome. The only notable distinction between the two factors was that FOXA1, the archetypal pioneering factor, required fewer copies of its target sequence to bind DNA than HNF4A.

These findings led Hansen et al. to propose an alternative theory to the pioneer factor hypothesis which eliminates the categorical distinction between pioneer and non-pioneer factors. Overall, this work has implications for how biologists understand the way that transcription factors activate silent genes during development.

Changes in chromatin accessibility landscape and histone H3 core acetylation during valproic acid-induced differentiation of embryonic stem cells

Directed differentiation of mouse embryonic stem cells (mESCs) or induced pluripotent stem cells (iPSCs) provides powerful models to dissect the molecular mechanisms leading to the formation of specific cell lineages. Treatment with histone deacetylase inhibitors can significantly enhance the efficiency of directed differentiation. However, the mechanisms are not well understood. Here, we use CUT&RUN in combination with ATAC-seq to determine changes in both histone modifications and genome-wide chromatin accessibility following valproic acid (VPA) exposure. VPA induced a significant increase in global histone H3 acetylation (H3K56ac), a core histone modification affecting nucleosome stability, as well as enrichment at loci associated with cytoskeletal organization and cellular morphogenesis. In addition, VPA altered the levels of linker histone H1 subtypes and the total histone H1/nucleosome ratio indicative of initial differentiation events. Notably, ATAC-seq analysis revealed changes in chromatin accessibility of genes involved in regulation of CDK serine/threonine kinase activity and DNA duplex unwinding. Importantly, changes in chromatin accessibility were evident at several key genomic loci, such as the pluripotency factor Lefty, cardiac muscle troponin Tnnt2, and the homeodomain factor Hopx, which play critical roles in cardiomyocyte differentiation. Massive parallel transcription factor (TF) footprinting also indicates an increased occupancy of TFs involved in differentiation toward mesoderm and endoderm lineages and a loss of footprints of POU5F1/SOX2 pluripotency factors following VPA treatment. Our results provide the first genome-wide analysis of the chromatin landscape following VPA-induced differentiation in mESCs and provide new mechanistic insight into the intricate molecular processes that govern departure from pluripotency and early lineage commitment.

Dissecting OCT4 defines the role of nucleosome binding in pluripotency

Pioneer transcription factors such as OCT4 can target silent genes embedded in nucleosome-dense regions. How nucleosome interaction enables transcription factors to target chromatin and determine cell identity remains elusive. Here, we systematically dissect OCT4 to show that nucleosome binding is encoded within the DNA-binding domain and yet can be uncoupled from free-DNA binding. Furthermore, accelerating the binding kinetics of OCT4 to DNA enhances nucleosome binding. In cells, uncoupling nucleosome binding diminishes the ability of OCT4 to individually access closed chromatin, while more dynamic nucleosome binding results in expansive genome scanning within closed chromatin. However, both uncoupling and enhancing nucleosome binding are detrimental to inducing pluripotency from differentiated cells. Remarkably, stable interactions between OCT4 and nucleosomes are continuously required for maintaining the accessibility of pluripotency enhancers in stem cells. Our findings reveal how the affinity and residence time of OCT4-nucleosome complexes modulate chromatin accessibility during cell fate changes and maintenance.

Integrated Chromatin Accessibility and Transcriptome Landscapes of Doxorubicin-Resistant Breast Cancer Cells

BACKGROUND

Doxorubicin is one of the most effective chemotherapeutic drugs for breast cancer while its common drug resistance leads to poor patient prognosis and survival. Growing evidence indicate dynamically reorganized chromatin allows rapid access of the gene regulatory machinery to open genomic regions facilitating subsequent gene expression through direct transcription factor (TF) activation and regulatory element binding.

METHODS

To better understand the regulatory network underlying doxorubicin resistance in breast cancer cells, we explored the systematic alterations of chromatin accessibility and gene expression by the assay for transposase-accessible chromatin using sequencing (ATAC-seq) in combination with RNA sequencing, followed by integrative analysis to identify potential regulators and their targets associated with differentially accessible regions (DARs) in doxorubicin-resistant MCF7 (MCF7-DR) cells.

RESULTS

A total of 3,963 differentially expressed genes (DEGs) related to doxorubicin resistance were identified, including dramatically up-regulated MT1E, GSTP1, LDHB, significantly down-regulated TFF1, UBB, DSCAM-AS1, and histone-modifying enzyme coding genes HDAC2, EZH2, PRMT5, etc. By integrating with transcriptomic datasets, we identified 18,228 DARs in MCF7-DR cells compared to control, which were positively correlated with their nearest DEGs (r = 0.6). There were 11,686 increased chromatin-accessible regions, which were enriched in up-regulated genes related to diverse KEGG pathways, such as the cell cycle, regulation of actin cytoskeleton, signaling pathways of MAPK, PI3K/Akt and Hippo, which play essential roles in regulating cell apoptosis, proliferation, metabolism, and inflammatory responses. The 6,542 decreased chromatin-accessible regions were identified for the declined doxorubicin-associated biological processes, for instance, endocrine and insulin resistance, central carbon metabolism, signaling pathways of TGF-beta and P53. Combining data from TCGA, analyses of the DAR sequences associated with the DNA-binding motifs of significantly enriched TF families including AP-1, TEAD and FOX, indicated that the loss-function of FOXA1 might play a critical role in doxorubicin-resistant breast cancer cells (DOX-R BCCs).

CONCLUSION

These data exhibit the non-genetic landscape of chromatin accessibility and transcript levels in the DOX-R BCCs, and provide clear insights and resources for the detection of critical TFs and potential cis-regulatory elements-based putative therapeutic targets.

Histone tails cooperate to control the breathing of genomic nucleosomes

Genomic DNA is packaged in chromatin, a dynamic fiber variable in size and compaction. In chromatin, repeating nucleosome units wrap 145–147 DNA basepairs around histone proteins. Genetic and epigenetic regulation of genes relies on structural transitions in chromatin which are driven by intra- and inter-nucleosome dynamics and modulated by chemical modifications of the unstructured terminal tails of histones. Here we demonstrate how the interplay between histone H3 and H2A tails control ample nucleosome breathing motions. We monitored large openings of two genomic nucleosomes, and only moderate breathing of an engineered nucleosome in atomistic molecular simulations amounting to 24 μs. Transitions between open and closed nucleosome conformations were mediated by the displacement and changes in compaction of the two histone tails. These motions involved changes in the DNA interaction profiles of clusters of epigenetic regulatory aminoacids in the tails. Removing the histone tails resulted in a large increase of the amplitude of nucleosome breathing but did not change the sequence dependent pattern of the motions. Histone tail modulated nucleosome breathing is a key mechanism of chromatin dynamics with important implications for epigenetic regulation.

In the cell, the DNA is packed in chromatin. Chromatin is a highly dynamic fiber structure made of arrays of nucleosomes with different degrees of compaction. Each nucleosome has 145–147 basepairs of DNA wrapped around a protein octamer made of four unique histone proteins. Each histone is present twice and has a structured part and one or two disordered terminal tails. The regulation of gene expression in the cell and during cellular transitions depends on dynamic changes in chromatin structure. Chromatin dynamics are modulated by intra and inter nucleosome motions and by posttranslational chemical modifications of the histone tails. Here we reveal how histone tails control the intra nucleosome dynamics at atomic resolution. From extensive sampling of nucleosome dynamics in atomistic molecular simulations, we show that genomic nucleosomes breath more extensively than engineered ones and we describe how two histone tails cooperate to control nucleosome breathing through interactions between clusters of positively charged residues and the DNA. Nucleosome conformations with different degrees of opening are associated with different conformations, positions, and DNA interaction patterns of the tails. With this mechanism, we contribute to the understanding of chromatin dynamics at atomic resolution.

Epigenetic Landscape of Liquid Biopsy in Colorectal Cancer

Colorectal cancer (CRC) is one of the most common malignancies and is a major cause of cancer-related deaths worldwide. Thus, there is a clinical need to improve early detection of CRC and personalize therapy for patients with this disease. In the era of precision oncology, liquid biopsy has emerged as a major approach to characterize the circulating tumor elements present in body fluids, including cell-free DNA and RNA, circulating tumor cells, and extracellular vesicles. This non-invasive tool has allowed the identification of relevant molecular alterations in CRC patients, including some indicating the disruption of epigenetic mechanisms. Epigenetic alterations found in solid and liquid biopsies have shown great utility as biomarkers for early detection, prognosis, monitoring, and evaluation of therapeutic response in CRC patients. Here, we summarize current knowledge of the most relevant epigenetic mechanisms associated with cancer development and progression, and the implications of their deregulation in cancer cells and liquid biopsy of CRC patients. In particular, we describe the methodologies used to analyze these epigenetic alterations in circulating tumor material, and we focus on the clinical utility of epigenetic marks in liquid biopsy as tumor biomarkers for CRC patients. We also discuss the great challenges and emerging opportunities of this field for the diagnosis and personalized management of CRC patients.

Breaths, Twists, and Turns of Atomistic Nucleosomes

Gene regulation programs establish cellular identity and rely on dynamic changes in the structural packaging of genomic DNA. The DNA is packaged in chromatin, which is formed from arrays of nucleosomes displaying different degree of compaction and different lengths of inter-nucleosomal linker DNA. The nucleosome represents the repetitive unit of chromatin and is formed by wrapping 145-147 basepairs of DNA around an octamer of histone proteins. Each of the four histones is present twice and has a structured core and intrinsically disordered terminal tails. Chromatin dynamics are triggered by inter- and intra-nucleosome motions that are controlled by the DNA sequence, the interactions between the histone core and the DNA, and the conformations, positions, and DNA interactions of the histone tails. Understanding chromatin dynamics requires studying all these features at the highest possible resolution. For this, molecular dynamics simulations can be used as a powerful complement or alternative to experimental approaches, from which it is often very challenging to characterize the structural features and atomic interactions controlling nucleosome motions. Molecular dynamics simulations can be performed at different resolutions, by coarse graining the molecular system with varying levels of details. Here we review the successes and the remaining challenges of the application of atomic resolution simulations to study the structure and dynamics of nucleosomes and their complexes with interacting partners.

The Role and Specific Mechanism of OCT4 in Cancer Stem Cells: A Review

Recently, evidences show that cancer stem cells (CSCs) are a type of cancer cell group with self-renewal and play a huge role in tumor recurrence, metastasis, and drug resistance. Finding new treatment directions and targets for cancer prognosis and reducing mortality has become a top priority. OCT4, as a transcription factor, participates in maintaining the stem characteristics of CSCs, but the mechanism of OCT4 is often overlooked. In this review, we try to illustrate the mechanism by which OCT4 plays a role in CSCs from the perspective of genetic modification of OCT4, non-coding RNA, complexes and signaling pathways associated with OCT4. Our ultimate goal is to provide new targets for cancer treatment to prolong the survival of cancer patients.

Pioneer Transcription Factors Initiating Gene Network Changes

Pioneer transcription factors have the intrinsic biochemical ability to scan partial DNA sequence motifs that are exposed on the surface of a nucleosome and thus access silent genes that are inaccessible to other transcription factors. Pioneer factors subsequently enable other transcription factors, nucleosome remodeling complexes, and histone modifiers to engage chromatin, thereby initiating the formation of an activating or repressive regulatory sequence. Thus, pioneer factors endow the competence for fate changes in embryonic development, are essential for cellular reprogramming, and rewire gene networks in cancer cells. Recent studies with reconstituted nucleosomes in vitro and chromatin binding in vivo reveal that pioneer factors can directly perturb nucleosome structure and chromatin accessibility in different ways. This review focuses on our current understanding of the mechanisms by which pioneer factors initiate gene network changes and will ultimately contribute to our ability to control cell fates at will.

Nucleosome allostery in pioneer transcription factor binding

While recent experiments revealed that some pioneer transcription factors (TFs) can bind to their target DNA sequences inside a nucleosome, the binding dynamics of their target recognitions are poorly understood. Here we used the latest coarse-grained models and molecular dynamics simulations to study the nucleosome-binding procedure of the two pioneer TFs, Sox2 and Oct4. In the simulations for a strongly positioning nucleosome, Sox2 selected its target DNA sequence only when the target was exposed. Otherwise, Sox2 entropically bound to the dyad region nonspecifically. In contrast, Oct4 plastically bound on the nucleosome mainly in two ways. First, the two POU domains of Oct4 separately bound to the two parallel gyres of the nucleosomal DNA, supporting the previous experimental results of the partial motif recognition. Second, the POU_S domain of Oct4 favored binding on the acidic patch of histones. Then, simulating the TFs binding to a genomic nucleosome, the LIN28B nucleosome, we found that the recognition of a pseudo motif by Sox2 induced the local DNA bending and shifted the population of the rotational position of the nucleosomal DNA. The redistributed DNA phase, in turn, changed the accessibility of a distant TF binding site, which consequently affected the binding probability of a second Sox2 or Oct4. These results revealed a nucleosomal DNA-mediated allosteric mechanism, through which one TF binding event can change the global conformation, and effectively regulate the binding of another TF at distant sites. Our simulations provide insights into the binding mechanism of single and multiple TFs on the nucleosome.

Mechanisms of OCT4-SOX2 motif readout on nucleosomes

Transcription factors (TFs) regulate gene expression through chromatin where nucleosomes restrict DNA access. To study how TFs bind nucleosome-occupied motifs, we focused on the reprogramming factors OCT4 and SOX2 in mouse embryonic stem cells. We determined TF engagement throughout a nucleosome at base-pair resolution in vitro, enabling structure determination by cryo-electron microscopy at two preferred positions. Depending on motif location, OCT4 and SOX2 differentially distort nucleosomal DNA. At one position, OCT4-SOX2 removes DNA from histone H2A and histone H3; however, at an inverted motif, the TFs only induce local DNA distortions. OCT4 uses one of its two DNA-binding domains to engage DNA in both structures, reading out a partial motif. These findings explain site-specific nucleosome engagement by the pluripotency factors OCT4 and SOX2, and they reveal how TFs distort nucleosomes to access chromatinized motifs.