Quantification of Protein–Protein and Protein–DNA Interactions In Vivo, Using Fluorescence Recovery after Photobleaching

Inferences from FRAP data are model dependent: A subdiffusive analysis

Fluorescence recovery after photobleaching (FRAP) is a widely used biological experiment to study the kinetics of molecules that react and move randomly. Since the development of FRAP in the 1970s, many reaction-diffusion models have been used to interpret FRAP data. However, intracellular molecules are widely observed to move by anomalous subdiffusion instead of normal diffusion. In this article, we extend a popular reaction-diffusion model of FRAP to the case of subdiffusion modeled by a fractional diffusion equation. By analyzing this reaction-subdiffusion model, we show that FRAP data are consistent with both diffusive and subdiffusive motion in many scenarios. We illustrate this general result by fitting our model to FRAP data from glucocorticoid receptors in a cell nucleus. We further show that the assumed model of molecular motion (normal diffusion or subdiffusion) strongly impacts the biological parameter values inferred from a given experimentally observed FRAP curve. We additionally analyze our model in three simplified parameter regimes and discuss parameter identifiability for varying subdiffusion exponents.

Linker histones are fine-scale chromatin architects modulating developmental decisions in Arabidopsis

BACKGROUND

Chromatin provides a tunable platform for gene expression control. Besides the well-studied core nucleosome, H1 linker histones are abundant chromatin components with intrinsic potential to influence chromatin function. Well studied in animals, little is known about the evolution of H1 function in other eukaryotic lineages for instance plants. Notably, in the model plant Arabidopsis, while H1 is known to influence heterochromatin and DNA methylation, its contribution to transcription, molecular, and cytological chromatin organization remains elusive.

RESULTS

We provide a multi-scale functional study of Arabidopsis linker histones. We show that H1-deficient plants are viable yet show phenotypes in seed dormancy, flowering time, lateral root, and stomata formation-complemented by either or both of the major variants. H1 depletion also impairs pluripotent callus formation. Fine-scale chromatin analyses combined with transcriptome and nucleosome profiling reveal distinct roles of H1 on hetero- and euchromatin: H1 is necessary to form heterochromatic domains yet dispensable for silencing of most transposable elements; H1 depletion affects nucleosome density distribution and mobility in euchromatin, spatial arrangement of nanodomains, histone acetylation, and methylation. These drastic changes affect moderately the transcription but reveal a subset of H1-sensitive genes.

CONCLUSIONS

H1 variants have a profound impact on the molecular and spatial (nuclear) chromatin organization in Arabidopsis with distinct roles in euchromatin and heterochromatin and a dual causality on gene expression. Phenotypical analyses further suggest the novel possibility that H1-mediated chromatin organization may contribute to the epigenetic control of developmental and cellular transitions.

Using a model comparison approach to describe the assembly pathway for histone H1

Histones H1 or linker histones are highly dynamic proteins that diffuse throughout the cell nucleus and associate with chromatin (DNA and associated proteins). This binding interaction of histone H1 with the chromatin is thought to regulate chromatin organization and DNA accessibility to transcription factors and has been proven to involve a kinetic process characterized by a population that associates weakly with chromatin and rapidly dissociates and another population that resides at a binding site for up to several minutes before dissociating. When considering differences between these two classes of interactions in a mathematical model for the purpose of describing and quantifying the dynamics of histone H1, it becomes apparent that there could be several assembly pathways that explain the kinetic data obtained in living cells. In this work, we model these different pathways using systems of reaction-diffusion equations and carry out a model comparison analysis using FRAP (fluorescence recovery after photobleaching) experimental data from different histone H1 variants to determine the most feasible mechanism to explain histone H1 binding to chromatin. The analysis favors four different chromatin assembly pathways for histone H1 which share common features and provide meaningful biological information on histone H1 dynamics. We show, using perturbation analysis, that the explicit consideration of high- and low-affinity associations of histone H1 with chromatin in the favored assembly pathways improves the interpretation of histone H1 experimental FRAP data. To illustrate the results, we use one of the favored models to assess the kinetic changes of histone H1 after core histone hyperacetylation, and conclude that this post-transcriptional modification does not affect significantly the transition of histone H1 from a weakly bound state to a tightly bound state.

Molecular dynamics of histone H1

The H1 or linker histones bind dynamically to chromatin in living cells via a process that involves transient association with the nucleosome near the DNA entry/exit site followed by dissociation, translocation to a new location, and rebinding. The mean residency time of H1 on any given nucleosome is about a minute, which is much shorter than that of most core histones but considerably longer than that of most other chromatin-binding proteins, including transcription factors. Here we review recent advances in understanding the kinetic pathway of H1 binding and how it relates to linker histone structure and function. We also describe potential mechanisms by which the dynamic binding of H1 might contribute directly to the regulation of gene expression and discuss several situations for which there is experimental evidence to support these mechanisms. Finally, we review the evidence for the participation of linker histone chaperones in mediating H1 exchange.

Spatial heterogeneity of dynamics of H1 linker histone

Linker histone H1 participates in maintaining higher order chromatin structures. It is a dynamic protein that binds to DNA and exchanges rapidly with a mobile pool. Therefore, the dynamics of H1 were probed in the nuclei of intact, live cells, using an array of microscopy techniques: fluorescence recovery after photobleaching (FRAP), raster image correlation spectroscopy (RICS), fluorescence correlation spectroscopy (FCS), pair correlation functions (pCF) and fluorescence anisotropy. Combination of these techniques yielded information on H1 dynamics at small (1–100 μs: FCS, RICS, anisotropy), moderate (1–100 ms: FCS, RICS, pCF) and large (1–100 s: pCF and FRAP) time scales. These results indicate that the global movement of H1 in nuclei (at distances >1 µm) occurs at the time scale of seconds and is determined by processes other than diffusion. Moreover, a fraction of H1, which remains immobile at the time scale of tenths of seconds, is detectable. However, local (at distances <0.7 µm) H1 dynamics comprises a process occurring at a short (~3 ms) time scale and multiple processes occurring at longer (10–2,500 ms) scales. The former (fast) process (corresponding probably to H1 diffusion) is more pronounced in the nuclear regions characterized by low H1 concentration, but the latter (slow, attributable to H1 binding) in the regions of high H1 concentration. Furthermore, some regions in nuclei (possibly containing dense chromatin) may constitute barriers that impair or block movement of H1 histones within short (<1 µm) distances.

Parameter importance in FRAP acquisition and analysis: a simulation approach

Fluorescence recovery after photobleaching (FRAP) is a widespread technique used to determine intracellular reaction and diffusion parameters. In recent years, due to technical advances and an increasing number of mathematical models for analysis, there was a resurging interest in FRAP applications. However, care has to be taken when inverting parameters from such data. We study potential influences on FRAP acquisition and analysis like initial fluorescence distribution, membrane passage, and geometrical aspects. Monte Carlo simulations are employed for the investigation of reaction-diffusion processes to additionally include cases in which no analytical description is available. To assess the importance of influencing factors we apply a sensitivity method based on elementary effects providing an estimate for the global parameter space. The combination of simulations and sensitivity measure helps us to predict ranges of parameters used in acquisition and analysis for which a reliably inversion of reaction-diffusion parameters is possible. Using this approach, we show that FRAP data are highly susceptible to misinterpretation. However, by identifying the parameters of susceptibility, our analysis provides the means for taking measures to significantly improve FRAP data interpretation and analysis.

Pin1 promotes histone H1 dephosphorylation and stabilizes its binding to chromatin

The prolyl isomerase Pin1 stimulates the dephosphorylation of histone H1, stabilizing its binding to chromatin at transcriptionally active chromatin.

Histone H1 plays a crucial role in stabilizing higher order chromatin structure. Transcriptional activation, DNA replication, and chromosome condensation all require changes in chromatin structure and are correlated with the phosphorylation of histone H1. In this study, we describe a novel interaction between Pin1, a phosphorylation-specific prolyl isomerase, and phosphorylated histone H1. A sub-stoichiometric amount of Pin1 stimulated the dephosphorylation of H1 in vitro and modulated the structure of the C-terminal domain of H1 in a phosphorylation-dependent manner. Depletion of Pin1 destabilized H1 binding to chromatin only when Pin1 binding sites on H1 were present. Pin1 recruitment and localized histone H1 phosphorylation were associated with transcriptional activation independent of RNA polymerase II. We thus identify a novel form of histone H1 regulation through phosphorylation-dependent proline isomerization, which has consequences on overall H1 phosphorylation levels and the stability of H1 binding to chromatin.

Are assumptions about the model type necessary in reaction-diffusion modeling? A FRAP application

At present, fluorescence recovery after photobleaching (FRAP) data are interpreted using various types of reaction-diffusion (RD) models: the model type is usually fixed first, and corresponding model parameters are inferred subsequently. In this article, we describe what we believe to be a novel approach for RD modeling without using any assumptions of model type or parameters. To the best of our knowledge, this is the first attempt to address both model-type and parameter uncertainties in inverting FRAP data. We start from the most general RD model, which accounts for a flexible number of molecular fractions, all mobile, with different diffusion coefficients. The maximal number of possible binding partners is identified and optimal parameter sets for these models are determined in a global search of the parameter-space using the Simulated Annealing strategy. The numerical performance of the described techniques was assessed using artificial and experimental FRAP data. Our general RD model outperformed the standard RD models used previously in modeling FRAP measurements and showed that intracellular molecular mobility can only be described adequately by allowing for multiple RD processes. Therefore, it is important to search not only for the optimal parameter set but also for the optimal model type.

Core histone hyperacetylation impacts cooperative behavior and high-affinity binding of histone H1 to chromatin

Linker histones stabilize higher order chromatin structures and limit access to proteins involved in DNA-dependent processes. Core histone acetylation is thought to modulate H1 binding. In the current study, we employed kinetic modeling of H1 recovery curves obtained during fluorescence recovery after photobleaching (FRAP) experiments to determine the impact of core histone acetylation on the different variants of H1. Following brief treatments with histone deacetylase inhibitor, most variants showed no change in H1 dynamics. A change in mobility was detected only when longer treatments were used to induce high levels of histone acetylation. This hyperacetylation imparted marked changes in the dynamics of low-affinity H1 population, while conferring variant-specific changes in the mobility of H1 molecules that were strongly bound. Both the C-terminal domain (CTD) and globular domain were responsible for this differential response to TSA. Furthermore, we found that neither the CTD nor the globular domain, by themselves, undergoes a change in kinetics following hyperacetylation. This led us to conclude that hyperacetylation of core histones affects the cooperative nature of low-affinity H1 binding, with some variants undergoing a predicted decrease of almost 2 orders of magnitude.

Quantitative interpretation of binding reactions of rapidly diffusing species using fluorescence recovery after photobleaching

Fluorescence recovery after photobleaching (FRAP) measurements offer an important tool for analyzing diffusion and binding processes. Confocal scanning laser microscopes that are used in FRAP experiments bleach regions with a radially Gaussian distributed profile. Previous attempts to derive analytical expressions in the case of processes governed by fast diffusion have overlooked the characteristics of the instruments used to perform FRAP measurements and therefore led to approximating solutions. In the present paper, bleaching laser beam characteristics are incorporated into an improved model to provide a more rigorous and accurate method. The proposed model simulates binding inside bounded regions, and it leads to FRAP curves that depend on the on and off rates that can be employed to determine the rate constants. It can be used in conjunction with experimental data acquired with confocal scanning laser microscopes to investigate the biophysical properties of proteins in living cells. The model aims to improve the accuracy when determining rate constants by taking into account a more realistic scenario of the light-matter interaction.

Molecular dynamics of histone H1This paper is one of a selection of papers published in this Special Issue, entitled CSBMCB’s 51st Annual Meeting – Epigenetics and Chromatin Dynamics, and has undergone the Journal’s usual peer review process

The histone H1 family of nucleoproteins represents an important class of structural and architectural proteins that are responsible for maintaining and stabilizing higher-order chromatin structure. Essential for mammalian cell viability, they are responsible for gene-specific regulation of transcription and other DNA-dependent processes. In this review, we focus on the wealth of information gathered on the molecular kinetics of histone H1 molecules using novel imaging techniques, such as fluorescence recovery after photobleaching. These experiments have shed light on the effects of H1 phosphorylation and core histone acetylation in influencing chromatin structure and dynamics. We also delineate important concepts surrounding the C-terminal domain of H1, such as the intrinsic disorder hypothesis, and how it affects H1 function. Finally, we address the biochemical mechanisms behind low-affinity H1 binding.

Visualizing chromatin dynamics in intact cells

Chromatin and associated regulatory proteins regulate gene expression in the natural environment of the intact cell nucleus. Specific combinations of DNA-binding transcription factors and recruited coregulatory proteins alter the conformation of chromatin at promoters and enhancers of target genes to stimulate or repress transcription. The dynamic nature of the regulatory proteins active in these processes allows the cell to modulate gene expression very rapidly, an important feature in many physiological processes. Live cell imaging and photobleaching studies of fluorescently-tagged proteins reveal that many transcription factors and other chromatin-associated proteins rapidly move through the nucleoplasm. Transcription factors also transiently interact with specific regulatory sequences in chromatin, suggesting that gene activation does not require the formation of stable long-lived regulatory complexes on the chromatin. In this review we discuss how dynamic interactions allow transcriptional regulatory proteins find their targets within the nucleus, alter target chromatin structure, and modulate physiological gene expression.

Protein kinase D-mediated phosphorylation and nuclear export of sphingosine kinase 2

Sphingosine kinase (SPHK) is a key enzyme producing important messenger sphingosine 1-phosphate and is implicated in cell proliferation and suppression of apoptosis. Because the extent of agonist-induced activation of SPHK is modest, signaling via SPHK may be regulated through its localization at specific intracellular sites. Although the SPHK1 isoform has been extensively studied and characterized, the regulation of expression and function of the other isoform, SPHK2, remain largely unexplored. Here we describe an important post-translational modification, namely, phosphorylation of SPHK2 catalyzed by protein kinase D (PKD), which regulates its localization. Upon stimulation of HeLa cells by tumor promoter phorbol 12-myristate 13-acetate, a serine residue in a novel and putative nuclear export signal, identified for the first time, in SPHK2 was phosphorylated followed by SPHK2 export from the nucleus. Constitutively active PKD phosphorylated this serine residue in the nuclear export signal both in vivo and in vitro. Moreover, down-regulation of PKDs through RNA interference resulted in the attenuation of both basal and phorbol 12-myristate 13-acetate-induced phosphorylation, which was followed by the accumulation of SPHK2 in the nucleus in a manner rescued by PKD over-expression. These results indicate that PKD is a physiologically relevant enzyme for SPHK2 phosphorylation, which leads to its nuclear export for subsequent cellular signaling.

Involvement of sphingosine-1-phosphate in glutamate secretion in hippocampal neurons

Neuronal activity greatly influences the formation and stabilization of synapses. Although receptors for sphingosine-1-phosphate (S1P), a lipid mediator regulating diverse cellular processes, are abundant in the central nervous system, neuron-specific functions of S1P remain largely undefined. Here, we report two novel actions of S1P using primary hippocampal neurons as a model system: (i) as a secretagogue where S1P triggers glutamate secretion and (ii) as an enhancer where S1P potentiates depolarization-evoked glutamate secretion. Sphingosine kinase 1 (SK1), a key enzyme for S1P production, was enriched in functional puncta of hippocampal neurons. Silencing SK1 expression by small interfering RNA as well as SK1 inhibition by dimethylsphingosine resulted in a strong inhibition of depolarization-evoked glutamate secretion. Fluorescence recovery after photobleaching analysis showed translocation of SK1 from cytosol to membranes at the puncta during depolarization, which resulted in subsequent accumulation of S1P within cells. Fluorescent resonance energy transfer analysis demonstrated that the S1P(1) receptor at the puncta was activated during depolarization and that depolarization-induced S1P(1) receptor activation was inhibited in SK1-knock-down cells. Importantly, exogenously added S1P at a nanomolar concentration by itself elicited glutamate secretion from hippocampal cells even when the Na(+)-channel was blocked by tetrodotoxin, suggesting that S1P acts on presynaptic membranes. Furthermore, exogenous S1P at a picomolar level potentiated depolarization-evoked secretion in the neurons. These findings indicate that S1P, through its autocrine action, facilitates glutamate secretion in hippocampal neurons both by secretagogue and enhancer actions and may be involved in mechanisms underlying regulation of synaptic transmission.

A biological interpretation of transient anomalous subdiffusion. I. Qualitative model

Anomalous subdiffusion has been reported for two-dimensional diffusion in the plasma membrane and three-dimensional diffusion in the nucleus and cytoplasm. If a particle diffuses in a suitable infinite hierarchy of binding sites, diffusion is well known to be anomalous at all times. But if the hierarchy is finite, diffusion is anomalous at short times and normal at long times. For a prescribed set of binding sites, Monte Carlo calculations yield the anomalous diffusion exponent and the average time over which diffusion is anomalous. If even a single binding site is present, there is a very short, almost artifactual, period of anomalous subdiffusion, but a hierarchy of binding sites extends the anomalous regime considerably. As is well known, an essential requirement for anomalous subdiffusion due to binding is that the diffusing particle cannot be in thermal equilibrium with the binding sites; an equilibrated particle diffuses normally at all times. Anomalous subdiffusion due to barriers, however, still occurs at thermal equilibrium, and anomalous subdiffusion due to a combination of binding sites and barriers is reduced but not eliminated on equilibration. This physical model is translated directly into a plausible biological model testable by single-particle tracking.

In vivo assay of presynaptic microtubule cytoskeleton dynamics in Drosophila

Disrupted microtubule dynamics in neuronal synapses has been suggested as an underlying cause for several devastating neurological diseases, including Hereditary Spastic Paraplegia (HSP) and Fragile X Syndrome (FXS). However, previous studies have been restricted to indirect assays of synaptic microtubules, i.e. immunocytochemistry of microtubule-associated proteins and post-translationally modified tubulins characteristic of microtubules with different stabilities. Very little is known about synaptic microtubule dynamics in vivo, or how microtubule dynamics may be disrupted in disease states. In this study, we develop methods to analyze microtubule dynamics directly in living synaptic boutons in situ. We use fluorescence recovery after photobleaching (FRAP) of transgenic green fluorescent protein (GFP) tagged tubulin at the well-characterized Drosophila neuromuscular junction (NMJ) synapse. FRAP measurements of tubulin-GFP demonstrate biphasic recovery kinetics. Treatment with taxol to stabilize microtubules and promote microtubule assembly reduces both recovery phases. Treatment with vinblastine to disassemble microtubules increases the fast recovery phase and decreases the slow recovery phase. These data indicate that the fast recovery phase is generated by rapid diffusion of tubulin subunits and the slow phase is generated by the relatively slow turnover of microtubules. This study demonstrates that tubulin-GFP fluorescence recovery after photobleaching can be used to assay microtubule dynamics directly in living synapses.

Analysis of binding at a single spatially localized cluster of binding sites by fluorescence recovery after photobleaching

Cells contain many subcellular structures in which specialized proteins locally cluster. Binding interactions within such clusters may be analyzed in live cells using models for fluorescence recovery after photobleaching (FRAP). Here we analyze a three-dimensional FRAP model that accounts for a single spatially localized cluster of binding sites in the presence of both diffusion and impermeable boundaries. We demonstrate that models completely ignoring the spatial localization of binding yield poor estimates for the binding parameters within the binding site cluster. In contrast, we find that ignoring only the restricted axial height of the binding-site cluster is far less detrimental, thereby enabling the use of computationally less expensive models. We also identify simplified solutions to the FRAP model for limiting behaviors where either diffusion or binding dominate. We show how ignoring a role for diffusion can sometimes produce serious errors in binding parameter estimation. We illustrate application of the method by analyzing binding of a transcription factor, the glucocorticoid receptor, to a tandem array of mouse mammary tumor virus promoter sites in live cells, obtaining an estimate for an in vivo binding constant (10(-7) M), and a first approximation of an upper bound on the transcription-factor residence time at the promoter (approximately 170 ms). These FRAP analysis tools will be important for measuring key cellular binding parameters necessary for a complete and accurate description of the networks that regulate cellular behavior.

Nucleoplasmic beta-actin exists in a dynamic equilibrium between low-mobility polymeric species and rapidly diffusing populations

β-Actin, once thought to be an exclusively cytoplasmic protein, is now known to have important functions within the nucleus. Nuclear β-actin associates with and functions in chromatin remodeling complexes, ribonucleic acid polymerase complexes, and at least some ribonucleoproteins. Proteins involved in regulating actin polymerization are also found in the interphase nucleus. We define the dynamic properties of nuclear actin molecules using fluorescence recovery after photobleaching. Our results indicate that actin and actin-containing complexes are reduced in their mobility through the nucleoplasm diffusing at ∼0.5 μm² s⁻¹. We also observed that ∼20% of the total nuclear actin pool has properties of polymeric actin that turns over rapidly. This pool could be detected in endogenous nuclear actin by using fluorescent polymeric actin binding proteins and was sensitive to drugs that alter actin polymerization. Our results validate previous reports of polymeric forms of nuclear actin observed in fixed specimens and reveal that these polymeric forms are very dynamic.

The DNA binding activity of p53 displays reaction-diffusion kinetics

The tumor suppressor protein p53 plays a key role in maintaining the genomic stability of mammalian cells and preventing malignant transformation. In this study, we investigated the intracellular diffusion of a p53-GFP fusion protein using confocal fluorescence recovery after photobleaching. We show that the diffusion of p53-GFP within the nucleus is well described by a mathematical model for diffusion of particles that bind temporarily to a spatially homogeneous immobile structure with binding and release rates k1 and k2, respectively. The diffusion constant of p53-GFP was estimated to be Dp53-GFP=15.4 microm2 s-1, significantly slower than that of GFP alone, DGFP=41.6 microm2 s-1. The reaction rates of the binding and unbinding of p53-GFP were estimated as k1=0.3 s-1 and k2=0.4 s-1, respectively, values suggestive of nonspecific binding. Consistent with this finding, the diffusional mobilities of tumor-derived sequence-specific DNA binding mutants of p53 were indistinguishable from that of the wild-type protein. These data are consistent with a model in which, under steady-state conditions, p53 is latent and continuously scans DNA, requiring activation for sequence-specific DNA binding.

Dynamic relocation of poly(ADP-ribose) glycohydrolase isoforms during radiation-induced DNA damage

Poly(ADP-ribosyl)ation is a very early cellular response to DNA damage. Poly(ADP-ribose) (PAR) accumulation is transient since PAR is rapidly hydrolyzed by poly(ADP-ribose) glycohydrolase (PARG). PARG may play a prominent role in DNA damage response and repair by removing PAR from modified proteins including PARP-1. Using living cells, we provide evidence that in response to DNA damage induced by gamma-irradiation the cytoplasmic 103 kDa PARG isoform translocates into the nucleus. We further observed that the nuclear GFP-hPARG110 enzyme relocalizes to the cytoplasm in response to DNA damage. Using different GFP-PARG fusion proteins specific for the nuclear and cytoplasmic forms, we demonstrate their dynamic distribution between cytoplasm and nucleoplasm and a high mobility of major PARG isoforms by fluorescence recovery after photobleaching (FRAP). The dynamic relocation of all PARG isoforms presented in this report reveals a novel biological mechanism by which PARG could be involved in DNA damage response.

H1 family histones in the nucleus. Control of binding and localization by the C-terminal domain

H1 histones bind to DNA as they enter and exit the nucleosome. H1 histones have a tripartite structure consisting of a short N-terminal domain, a highly conserved central globular domain, and a lysine-and arginine-rich C-terminal domain. The C-terminal domain comprises approximately half of the total amino acid content of the protein, is essential for the formation of compact chromatin structures, and contains the majority of the amino acid variations that define the individual histone H1 family members. This region contains several cell cycle-regulated phosphorylation sites and is thought to function through a charge-neutralization process, neutralizing the DNA phosphate backbone to allow chromatin compaction. In this study, we use fluorescence microscopy and fluorescence recovery after photobleaching to define the behavior of the individual histone H1 subtypes in vivo. We find that there are dramatic differences in the binding affinity of the individual histone H1 subtypes in vivo and differences in their preference for euchromatin and heterochromatin. Further, we show that subtype-specific properties originate with the C terminus and that the differences in histone H1 binding are not consistent with the relatively small changes in the net charge of the C-terminal domains.