Recovering mixtures of fast-diffusing states from short single-particle trajectories

Epigenetic memory is governed by an effector recruitment specificity toggle in Heterochromatin Protein 1

HP1 proteins are essential for establishing and maintaining transcriptionally silent heterochromatin. They dimerize, forming a binding interface to recruit diverse chromatin-associated factors. HP1 proteins are specialized and rapidly evolve, but the extent of variation required to achieve functional specialization is unknown. To investigate how changes in amino acid sequence impacts epigenetic inheritance, we performed a targeted mutagenesis screen of the S. pombe HP1 homolog, Swi6. Substitutions within an auxiliary surface adjacent to the HP1 dimerization interface produced Swi6 variants with divergent maintenance properties. Remarkably, substitutions at a single amino acid position led to the persistent gain or loss of epigenetic inheritance. These substitutions increased Swi6 chromatin occupancy in vivo and altered Swi6-protein interactions that reprogram H3K9me maintenance. We show that relatively minor changes in Swi6 amino acid composition can lead to profound changes in epigenetic inheritance which provides a redundant mechanism to evolve novel effector specificity. .

Diffusion barriers imposed by tissue topology shape morphogen gradients

Animals use a small number of morphogens to pattern tissues, but it is unclear how evolution modulates morphogen signaling range to match tissues of varying sizes. Here, we used single molecule imaging in reconstituted morphogen gradients and in tissue explants to determine that Hedgehog diffused extra-cellularly as a monomer, and rapidly transitioned between membrane-confined and -unconfined states. Unexpectedly, the vertebrate-specific protein SCUBE1 expanded Hedgehog gradients by accelerating the transition rates between states without affecting the relative abundance of molecules in each state. This observation could not be explained under existing models of morphogen diffusion. Instead, we developed a topology-limited diffusion model in which cell-cell gaps create diffusion barriers, and morphogens can only overcome the barrier by passing through a membrane-unconfined state. Under this model, SCUBE1 promotes Hedgehog secretion and diffusion by allowing it to transiently overcome diffusion barriers. This multiscale understanding of morphogen gradient formation unified prior models and discovered novel knobs that nature can use to tune morphogen gradient sizes across tissues and organisms.

A fine kinetic balance of interactions directs transcription factor hubs to genes

Eukaryotic gene regulation relies on the binding of sequence-specific transcription factors (TFs). TFs bind chromatin transiently yet occupy their target sites by forming high-local concentration microenvironments (hubs and condensates) that increase the frequency of binding events. Despite their ubiquity, such microenvironments have been difficult to study in endogenous contexts due to technical limitations. Here, we overcome these limitations and investigate how hubs drive TF occupancy at their targets. Using a DNA binding perturbation to a hub-forming TF, Zelda, in Drosophila embryos, we find that hub properties, including the stability and frequencies of associations to targets, are key determinants of TF occupancy. Our data suggest that the targeting of these hubs is driven not just by specific DNA motif recognition, but also by a fine-tuned kinetic balance of interactions between TFs and their co-binding partners.

Disordered C-terminal domain drives spatiotemporal confinement of RNAPII to enhance search for chromatin targets

Efficient gene expression requires RNA polymerase II (RNAPII) to find chromatin targets precisely in space and time. How RNAPII manages this complex diffusive search in three-dimensional nuclear space remains largely unknown. The disordered carboxy-terminal domain (CTD) of RNAPII, which is essential for recruiting transcription-associated proteins, forms phase-separated droplets in vitro, hinting at a potential role in modulating RNAPII dynamics. In the present study, we use single-molecule tracking and spatiotemporal mapping in living yeast to show that the CTD is required for confining RNAPII diffusion within a subnuclear region enriched for active genes, but without apparent phase separation into condensates. Both Mediator and global chromatin organization are required for sustaining RNAPII confinement. Remarkably, truncating the CTD disrupts RNAPII spatial confinement, prolongs target search, diminishes chromatin binding, impairs pre-initiation complex formation and reduces transcription bursting. The present study illuminates the pivotal role of the CTD in driving spatiotemporal confinement of RNAPII for efficient gene expression.

Dynamic microenvironments shape nuclear organization and gene expression

Live imaging has revealed that the regulation of gene expression is largely driven by transient interactions. For example, many regulatory proteins bind chromatin for just seconds, and loop-like genomic contacts are rare and last only minutes. These discoveries have been difficult to reconcile with our canonical models that are predicated on stable and hierarchical interactions. Proteomic microenvironments that concentrate nuclear factors may explain how brief interactions can still mediate gene regulation by creating conditions where reactions occur more frequently. Here, we summarize new imaging technologies and recent discoveries implicating microenvironments as a potential driver of nuclear function. Finally, we propose that key properties of proteomic microenvironments, such as their size, enrichment, and lifetimes, are directly linked to regulatory function.

p300 is an obligate integrator of combinatorial transcription factor inputs

Transcription coactivators are proteins or protein complexes that mediate transcription factor (TF) function. However, they lack DNA-binding capacity, prompting the question of how they engage target loci. Three non-exclusive hypotheses have been posited: coactivators are recruited by complexing with TFs, by binding histones through epigenetic reader domains, or by partitioning into condensates through their extensive intrinsically disordered regions. Using p300 as a prototypical coactivator, we systematically mutated its annotated domains and show by single-molecule tracking in live U2OS cells that coactivator-chromatin binding depends entirely on combinatorial binding of multiple TF-interaction domains. Furthermore, we demonstrate that acetyltransferase activity opposes p300-chromatin association and that the N-terminal TF-interaction domains regulate that activity. Single TF-interaction domains are insufficient for chromatin binding and regulation of catalytic activity, implying a principle that we speculate could broadly apply to eukaryotic gene regulation: a TF must act in coordination with other TFs to recruit coactivator activity.

Extract latent features of single-particle trajectories with historical experience learning

Single-particle tracking has enabled real-time, in situ quantitative studies of complex systems. However, inferring dynamic state changes from noisy and undersampling trajectories encounters challenges. Here, we introduce a data-driven method for extracting features of subtrajectories with historical experience learning (Deep-SEES), where a single-particle tracking analysis pipeline based on a self-supervised architecture automatically searches for the latent space, allowing effective segmentation of the underlying states from noisy trajectories without prior knowledge on the particle dynamics. We validated our method on a variety of noisy simulated and experimental data. Our results showed that the method can faithfully capture both stable states and their dynamic switch. In highly random systems, our method outperformed commonly used unsupervised methods in inferring motion states, which is important for understanding nanoparticles interacting with living cell membranes, active enzymes, and liquid-liquid phase separation. Self-generating latent features of trajectories could potentially improve the understanding, estimation, and prediction of many complex systems.

Quantitative analysis of peroxisome tracks using a Hidden Markov Model

Diffusion and mobility are essential for cellular functions, as molecules are usually distributed throughout the cell and have to meet to interact and perform their function. This also involves the cytosolic migration of cellular organelles. However, observing such diffusion and interaction dynamics is challenging due to the high spatial and temporal resolution required and the accurate analysis of the diffusional tracks. The latter is especially important when identifying anomalous diffusion events, such as directed motions, which are often rare. Here, we investigate the migration modes of peroxisome organelles in the cytosol of living cells. Peroxisomes predominantly migrate randomly, but occasionally they bind to the cell's microtubular network and perform directed migration, which is difficult to quantify, and so far, accurate analysis of switching between these migration modes is missing. We set out to solve this limitation by experiments and analysis with high statistical accuracy. Specifically, we collect temporal diffusion tracks of thousands of individual peroxisomes in the HEK 293 cell line using two-dimensional spinning disc fluorescence microscopy at a high acquisition rate of 10 frames/s. We use a Hidden Markov Model with two hidden states to (1) automatically identify directed migration segments of the tracks and (2) quantify the migration properties for comparison between states and between different experimental conditions. Comparing different cellular conditions, we show that the knockout of the peroxisomal membrane protein PEX14 leads to a decrease in the directed movement due to a lowered binding probability to the microtubule. However, it does not eradicate binding, highlighting further microtubule-binding mechanisms of peroxisomes than via PEX14. In contrast, structural changes of the microtubular network explain perceived eradication of directed movement by disassembly of microtubules by Nocodazole-treatment.

Tracking live-cell single-molecule dynamics enables measurements of heterochromatinassociated protein-protein interactions

Visualizing and measuring molecular-scale interactions in living cells represents a major challenge, but recent advances in microscopy are bringing us closer to achieving this goal. Single-molecule super-resolution microscopy enables high-resolution and sensitive imaging of the positions and movement of molecules in living cells. HP1 proteins are important regulators of gene expression because they selectively bind and recognize H3K9 methylated (H3K9me) histones to form heterochromatin-associated protein complexes that silence gene expression. Here, we extended live-cell single-molecule tracking studies in fission yeast to determine how HP1 proteins interact with their binding partners in the nucleus. We measured how genetic perturbations that affect H3K9me alter the diffusive properties of HP1 proteins and each of their binding partners based on which we inferred their most likely interaction sites. Our results indicate that H3K9me promotes specific complex formation between HP1 proteins and their interactors in a spatially restricted manner, while attenuating their ability to form off-chromatin complexes. As opposed to being an inert platform or scaffold to direct HP1 binding, our studies propose a novel function for H3K9me as an active participant in enhancing HP1-associated complex formation in living cells.

Surprising Features of Nuclear Receptor Interaction Networks Revealed by Live Cell Single Molecule Imaging

Type 2 Nuclear Receptors (T2NRs) require heterodimerization with a common partner, the Retinoid X Receptor (RXR), to bind cognate DNA recognition sites in chromatin. Based on previous biochemical and over-expression studies, binding of T2NRs to chromatin is proposed to be regulated by competition for a limiting pool of the core RXR subunit. However, this mechanism has not yet been tested for endogenous proteins in live cells. Using single molecule tracking (SMT) and proximity-assisted photoactivation (PAPA), we monitored interactions between endogenously tagged retinoid X receptor (RXR) and retinoic acid receptor (RAR) in live cells. Unexpectedly, we find that higher expression of RAR, but not RXR increases heterodimerization and chromatin binding in U2OS cells. This surprising finding indicates the limiting factor is not RXR but likely its cadre of obligate dimer binding partners. SMT and PAPA thus provide a direct way to probe which components are functionally limiting within a complex TF interaction network providing new insights into mechanisms of gene regulation in vivo with implications for drug development targeting nuclear receptors.

A novel pathogenic mutation of MeCP2 impairs chromatin association independent of protein levels

Loss-of-function mutations in MECP2 cause Rett syndrome (RTT), a severe neurological disorder that mainly affects girls. Mutations in MECP2 do occur in males occasionally and typically cause severe encephalopathy and premature lethality. Recently, we identified a missense mutation (c.353G>A, p.Gly118Glu [G118E]), which has never been seen before in MECP2, in a young boy who suffered from progressive motor dysfunction and developmental delay. To determine whether this variant caused the clinical symptoms and study its functional consequences, we established two disease models, including human neurons from patient-derived iPSCs and a knock-in mouse line. G118E mutation partially reduces MeCP2 abundance and its DNA binding, and G118E mice manifest RTT-like symptoms seen in the patient, affirming the pathogenicity of this mutation. Using live-cell and single-molecule imaging, we found that G118E mutation alters MeCP2's chromatin interaction properties in live neurons independently of its effect on protein levels. Here we report the generation and characterization of RTT models of a male hypomorphic variant and reveal new insight into the mechanism by which this pathological mutation affects MeCP2's chromatin dynamics. Our ability to quantify protein dynamics in disease models lays the foundation for harnessing high-resolution single-molecule imaging as the next frontier for developing innovative therapies for RTT and other diseases.

Absolute quantification of protein number and dynamics in single cells

Quantitative characterization of protein abundance and interactions in live cells is necessary to understand and predict cellular behavior. The accurate determination of copy number for individual proteins and heterologous complexes in individual cells is critical because small changes in protein dosage, often less than two-fold, can have strong phenotypic consequences. Here, we review the merits and pitfalls of different quantitative fluorescence imaging methods for single-cell determination of protein abundance, localization, interactions, and dynamics. In particular, we discuss how scanning number and brightness (sN&B) and its variation, Raster scanning image correlation spectroscopy (RICS), exploit stochastic noise in small measurement volumes to quantify protein abundance, stoichiometry, and dynamics with high accuracy.

Single-molecule tracking to determine the abundances and stoichiometries of freely-diffusing protein complexes in living cells: Past applications and future prospects

Most biological processes in living cells rely on interactions between proteins. Live-cell compatible approaches that can quantify to what extent a given protein participates in homo- and hetero-oligomeric complexes of different size and subunit composition are therefore critical to advance our understanding of how cellular physiology is governed by these molecular interactions. Biomolecular complex formation changes the diffusion coefficient of constituent proteins, and these changes can be measured using fluorescence microscopy-based approaches, such as single-molecule tracking, fluorescence correlation spectroscopy, and fluorescence recovery after photobleaching. In this review, we focus on the use of single-molecule tracking to identify, resolve, and quantify the presence of freely-diffusing proteins and protein complexes in living cells. We compare and contrast different data analysis methods that are currently employed in the field and discuss experimental designs that can aid the interpretation of the obtained results. Comparisons of diffusion rates for different proteins and protein complexes in intracellular aqueous environments reported in the recent literature reveal a clear and systematic deviation from the Stokes-Einstein diffusion theory. While a complete and quantitative theoretical explanation of why such deviations manifest is missing, the available data suggest the possibility of weighing freely-diffusing proteins and protein complexes in living cells by measuring their diffusion coefficients. Mapping individual diffusive states to protein complexes of defined molecular weight, subunit stoichiometry, and structure promises to provide key new insights into how protein-protein interactions regulate protein conformational, translational, and rotational dynamics, and ultimately protein function.

RNA polymerase II depletion from the inactive X chromosome territory is not mediated by physical compartmentalization

Subnuclear compartmentalization has been proposed to play an important role in gene regulation by segregating active and inactive parts of the genome in distinct physical and biochemical environments. During X chromosome inactivation (XCI), the noncoding Xist RNA coats the X chromosome, triggers gene silencing and forms a dense body of heterochromatin from which the transcription machinery appears to be excluded. Phase separation has been proposed to be involved in XCI, and might explain the exclusion of the transcription machinery by preventing its diffusion into the Xist-coated territory. Here, using quantitative fluorescence microscopy and single-particle tracking, we show that RNA polymerase II (RNAPII) freely accesses the Xist territory during the initiation of XCI. Instead, the apparent depletion of RNAPII is due to the loss of its chromatin stably bound fraction. These findings indicate that initial exclusion of RNAPII from the inactive X reflects the absence of actively transcribing RNAPII, rather than a consequence of putative physical compartmentalization of the inactive X heterochromatin domain.

Dynamic switching of transcriptional regulators between two distinct low-mobility chromatin states

How chromatin dynamics relate to transcriptional activity remains poorly understood. Using single-molecule tracking, coupled with machine learning, we show that histone H2B and multiple chromatin-bound transcriptional regulators display two distinct low-mobility states. Ligand activation results in a marked increase in the propensity of steroid receptors to bind in the lowest-mobility state. Mutational analysis revealed that interactions with chromatin in the lowest-mobility state require an intact DNA binding domain and oligomerization domains. These states are not spatially separated as previously believed, but individual H2B and bound-TF molecules can dynamically switch between them on time scales of seconds. Single bound-TF molecules with different mobilities exhibit different dwell time distributions, suggesting that the mobility of TFs is intimately coupled with their binding dynamics. Together, our results identify two unique and distinct low-mobility states that appear to represent common pathways for transcription activation in mammalian cells.

p300 Is an Obligate Integrator of Combinatorial Transcription Factor Inputs

Transcription coactivators are proteins or protein complexes that mediate transcription factor (TF) function. However, they lack DNA binding capacity, prompting the question of how they engage target loci. Three non-exclusive hypotheses have been posited: coactivators are recruited by complexing with TFs, by binding histones through epigenetic reader domains, or by partitioning into phase-separated compartments through their extensive intrinsically disordered regions (IDRs). Using p300 as a prototypical coactivator, we systematically mutated its annotated domains and show by single-molecule tracking in live cells that coactivator-chromatin binding depends entirely on combinatorial binding of multiple TF-interaction domains. Furthermore, we demonstrate that acetyltransferase activity negatively impacts p300-chromatin association and that the N-terminal TF-interaction domains regulate that activity. Single TF-interaction domains are insufficient for both chromatin binding and regulation of catalytic activity, implying a principle that could broadly inform eukaryotic gene regulation: a TF must act in coordination with other TFs to recruit coactivator activity.

Single-molecule tracking (SMT): a window into live-cell transcription biochemistry

How molecules interact governs how they move. Single-molecule tracking (SMT) thus provides a unique window into the dynamic interactions of biomolecules within live cells. Using transcription regulation as a case study, we describe how SMT works, what it can tell us about molecular biology, and how it has changed our perspective on the inner workings of the nucleus. We also describe what SMT cannot yet tell us and how new technical advances seek to overcome its limitations. This ongoing progress will be imperative to address outstanding questions about how dynamic molecular machines function in live cells.

Impact of Saccharomyces cerevisiae on the Field of Single-Molecule Biophysics

Cellular functions depend on the dynamic assembly of protein regulator complexes at specific cellular locations. Single Molecule Tracking (SMT) is a method of choice for the biochemical characterization of protein dynamics in vitro and in vivo. SMT follows individual molecules in live cells and provides direct information about their behavior. SMT was successfully applied to mammalian models. However, mammalian cells provide a complex environment where protein mobility depends on numerous factors that are difficult to control experimentally. Therefore, yeast cells, which are unicellular and well-studied with a small and completely sequenced genome, provide an attractive alternative for SMT. The simplicity of organization, ease of genetic manipulation, and tolerance to gene fusions all make yeast a great model for quantifying the kinetics of major enzymes, membrane proteins, and nuclear and cellular bodies. However, very few researchers apply SMT techniques to yeast. Our goal is to promote SMT in yeast to a wider research community. Our review serves a dual purpose. We explain how SMT is conducted in yeast cells, and we discuss the latest insights from yeast SMT while putting them in perspective with SMT of higher eukaryotes.

Mechanisms governing target search and binding dynamics of hypoxia-inducible factors

Transcription factors (TFs) are classically attributed a modular construction, containing well-structured sequence-specific DNA-binding domains (DBDs) paired with disordered activation domains (ADs) responsible for protein-protein interactions targeting co-factors or the core transcription initiation machinery. However, this simple division of labor model struggles to explain why TFs with identical DNA-binding sequence specificity determined in vitro exhibit distinct binding profiles in vivo. The family of hypoxia-inducible factors (HIFs) offer a stark example: aberrantly expressed in several cancer types, HIF-1α and HIF-2α subunit isoforms recognize the same DNA motif in vitro - the hypoxia response element (HRE) - but only share a subset of their target genes in vivo, while eliciting contrasting effects on cancer development and progression under certain circumstances. To probe the mechanisms mediating isoform-specific gene regulation, we used live-cell single particle tracking (SPT) to investigate HIF nuclear dynamics and how they change upon genetic perturbation or drug treatment. We found that HIF-α subunits and their dimerization partner HIF-1β exhibit distinct diffusion and binding characteristics that are exquisitely sensitive to concentration and subunit stoichiometry. Using domain-swap variants, mutations, and a HIF-2α specific inhibitor, we found that although the DBD and dimerization domains are important, another main determinant of chromatin binding and diffusion behavior is the AD-containing intrinsically disordered region (IDR). Using Cut&Run and RNA-seq as orthogonal genomic approaches, we also confirmed IDR-dependent binding and activation of a specific subset of HIF target genes. These findings reveal a previously unappreciated role of IDRs in regulating the TF search and binding process that contribute to functional target site selectivity on chromatin.

Recovering mixtures of fast-diffusing states from short single-particle trajectories

Single-particle tracking (SPT) directly measures the dynamics of proteins in living cells and is a powerful tool to dissect molecular mechanisms of cellular regulation. Interpretation of SPT with fast-diffusing proteins in mammalian cells, however, is complicated by technical limitations imposed by fast image acquisition. These limitations include short trajectory length due to photobleaching and shallow depth of field, high localization error due to the low photon budget imposed by short integration times, and cell-to-cell variability. To address these issues, we investigated methods inspired by Bayesian nonparametrics to infer distributions of state parameters from SPT data with short trajectories, variable localization precision, and absence of prior knowledge about the number of underlying states. We discuss the advantages and disadvantages of these approaches relative to other frameworks for SPT analysis.

Detecting molecular interactions in live-cell single-molecule imaging with proximity-assisted photoactivation (PAPA)

Single-molecule imaging provides a powerful way to study biochemical processes in live cells, yet it remains challenging to track single molecules while simultaneously detecting their interactions. Here, we describe a novel property of rhodamine dyes, proximity-assisted photoactivation (PAPA), in which one fluorophore (the 'sender') can reactivate a second fluorophore (the 'receiver') from a dark state. PAPA requires proximity between the two fluorophores, yet it operates at a longer average intermolecular distance than Förster resonance energy transfer (FRET). We show that PAPA can be used in live cells both to detect protein-protein interactions and to highlight a subpopulation of labeled protein complexes in which two different labels are in proximity. In proof-of-concept experiments, PAPA detected the expected correlation between androgen receptor self-association and chromatin binding at the single-cell level. These results establish a new way in which a photophysical property of fluorophores can be harnessed to study molecular interactions in single-molecule imaging of live cells.

Enhancer-promoter interactions and transcription are largely maintained upon acute loss of CTCF, cohesin, WAPL or YY1

It remains unclear why acute depletion of CTCF (CCCTC-binding factor) and cohesin only marginally affects expression of most genes despite substantially perturbing three-dimensional (3D) genome folding at the level of domains and structural loops. To address this conundrum, we used high-resolution Micro-C and nascent transcript profiling in mouse embryonic stem cells. We find that enhancer-promoter (E-P) interactions are largely insensitive to acute (3-h) depletion of CTCF, cohesin or WAPL. YY1 has been proposed as a structural regulator of E-P loops, but acute YY1 depletion also had minimal effects on E-P loops, transcription and 3D genome folding. Strikingly, live-cell, single-molecule imaging revealed that cohesin depletion reduced transcription factor (TF) binding to chromatin. Thus, although CTCF, cohesin, WAPL or YY1 is not required for the short-term maintenance of most E-P interactions and gene expression, our results suggest that cohesin may facilitate TFs to search for and bind their targets more efficiently.

Fixation can change the appearance of phase separation in living cells

Fixing cells with paraformaldehyde (PFA) is an essential step in numerous biological techniques as it is thought to preserve a snapshot of biomolecular transactions in living cells. Fixed-cell imaging techniques such as immunofluorescence have been widely used to detect liquid-liquid phase separation (LLPS) in vivo. Here, we compared images, before and after fixation, of cells expressing intrinsically disordered proteins that are able to undergo LLPS. Surprisingly, we found that PFA fixation can both enhance and diminish putative LLPS behaviors. For specific proteins, fixation can even cause their droplet-like puncta to artificially appear in cells that do not have any detectable puncta in the live condition. Fixing cells in the presence of glycine, a molecule that modulates fixation rates, can reverse the fixation effect from enhancing to diminishing LLPS appearance. We further established a kinetic model of fixation in the context of dynamic protein-protein interactions. Simulations based on the model suggest that protein localization in fixed cells depends on an intricate balance of protein-protein interaction dynamics, the overall rate of fixation, and notably, the difference between fixation rates of different proteins. Consistent with simulations, live-cell single-molecule imaging experiments showed that a fast overall rate of fixation relative to protein-protein interaction dynamics can minimize fixation artifacts. Our work reveals that PFA fixation changes the appearance of LLPS from living cells, presents a caveat in studying LLPS using fixation-based methods, and suggests a mechanism underlying the fixation artifact.