Investigation of binding mechanisms of nuclear proteins using confocal scanning laser microscopy and FRAP

Nuclear receptor corepressor 1 levels differentially impact the intracellular dynamics of mutant thyroid hormone receptors associated with resistance to thyroid hormone syndrome

Thyroid hormone receptor α1 (TRα1) undergoes nucleocytoplasmic shuttling and mediates gene expression in response to thyroid hormone (T3). In Resistance to Thyroid Hormone Syndrome α (RTHα), certain TRα1 mutants have higher affinity for nuclear corepressor 1 (NCoR1) and may form stable complexes that are not released in the presence of T3. Here, we examined whether NCoR1 modulates intranuclear mobility and nuclear retention of TRα1 or RTHα-associated mutants in transfected human cells, as a way of analyzing critical structural components of TRα1 and to further explore the correlation between mutations in TRα1 and aberrant intracellular trafficking. We found no significant difference in intranuclear mobility, as measured by fluorescence recovery after photobleaching, between TRα1 and select RTHα mutants, irrespective of NCoR1 expression. Nuclear-to-cytoplasmic fluorescence ratios of RTHα mutants, however, varied from TRα1 when NCoR1 was overexpressed, with a significant increase in nuclear retention for A263V and a significant decrease for A263S and R384H. In NCoR1-knockout cells, nuclear retention of A263S, A263V, P389R, A382P, C392X, and F397fs406X was significantly decreased compared to control (wild-type) cells. Luciferase reporter gene transcription mediated by TRα1 was significantly repressed by both NCoR1 overexpression and NCoR1 knockout. Most RTHα mutants showed minimal induction regardless of NCoR1 levels, but T3-mediated transcriptional activity was decreased for R384C and F397fs406X when NCoR1 was overexpressed, and also decreased for N359Y in NCoR1-knockout cells. Our results suggest a complex interaction between NCoR1 and RTHα mutants characterized by aberrant intracellular localization patterns and transcriptional activity that potentially arise from variable repressor complex stability, and may provide insight into RTHα pathogenesis on a molecular and cellular level.

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.

Fluorescence recovery after photobleaching in material and life sciences: putting theory into practice

Fluorescence recovery after photobleaching (FRAP) is a versatile tool for determining diffusion and interaction/binding properties in biological and material sciences. An understanding of the mechanisms controlling the diffusion requires a deep understanding of structure–interaction–diffusion relationships. In cell biology, for instance, this applies to the movement of proteins and lipids in the plasma membrane, cytoplasm and nucleus. In industrial applications related to pharmaceutics, foods, textiles, hygiene products and cosmetics, the diffusion of solutes and solvent molecules contributes strongly to the properties and functionality of the final product. All these systems are heterogeneous, and accurate quantification of the mass transport processes at the local level is therefore essential to the understanding of the properties of soft (bio)materials. FRAP is a commonly used fluorescence microscopy-based technique to determine local molecular transport at the micrometer scale. A brief high-intensity laser pulse is locally applied to the sample, causing substantial photobleaching of the fluorescent molecules within the illuminated area. This causes a local concentration gradient of fluorescent molecules, leading to diffusional influx of intact fluorophores from the local surroundings into the bleached area. Quantitative information on the molecular transport can be extracted from the time evolution of the fluorescence recovery in the bleached area using a suitable model. A multitude of FRAP models has been developed over the years, each based on specific assumptions. This makes it challenging for the non-specialist to decide which model is best suited for a particular application. Furthermore, there are many subtleties in performing accurate FRAP experiments. For these reasons, this review aims to provide an extensive tutorial covering the essential theoretical and practical aspects so as to enable accurate quantitative FRAP experiments for molecular transport measurements in soft (bio)materials.

Analysis of biomolecular dynamics by FRAP and computer simulation

We present a Monte Carlo simulation environment for modelling complex biological molecular interaction networks and for the design, validation, and quantitative analysis of FRAP assays to study these. The program is straightforward in its implementation and can be instructed through an intuitive script language. The simulation tool fits very well in a systems biology research setting that aims to maintain an interactive cycle of experiment-driven modelling and model-driven experimentation: the model and the experiment are in the same simulation. The full program can be obtained by request to the authors.

Nuclear proteins: finding and binding target sites in chromatin

Fluorescent protein labelling, as well as impressive progress in live cell imaging have revolutionised the view on how essential nuclear functions like gene transcription regulation and DNA repair are organised. Here, we address questions like how DNA-interacting molecules find and bind their target sequences in the vast amount of DNA. In addition, we discuss methods that have been developed for quantitative analysis of data from fluorescence recovery after photobleaching experiments (FRAP).

A nonfitting method using a spatial sine window transform for inhomogeneous effective-diffusion measurements by FRAP

Determining averaged effective diffusion constants from experimental measurements of fluorescent proteins in an inhomogeneous medium in the presence of ligand-receptor interactions poses problems of analytical tractability. Here, we introduced a nonfitting method to evaluate the averaged effective diffusion coefficient of a region of interest (which may include a whole nucleus) by mathematical processing of the entire cellular two-dimensional spatial pattern of recovered fluorescence. Spatially and temporally resolved measurements of protein transport inside cells were obtained using the fluorescence recovery after photobleaching technique. Two-dimensional images of fluorescence patterns were collected by laser-scanning confocal microscopy. The method was demonstrated by applying it to an estimation of the mobility of green fluorescent protein-tagged heterochromatin protein 1 in the nuclei of living mouse embryonic fibroblasts. This approach does not require the mathematical solution of a corresponding system of diffusion-reaction equations that is typical of conventional fluorescence recovery after photobleaching data processing, and is most useful for investigating highly inhomogeneous areas, such as cell nuclei, which contain many protein foci and chromatin domains.

A quantitative approach to analyze binding diffusion kinetics by confocal FRAP

Most of the important types of interactions that occur in cells can be characterized as binding-diffusion type processes, and can be quantified by kinetic rate constants such as diffusion coefficients (D) and binding rate constants (k(on) and k(off)). Confocal FRAP is a potentially important tool for the quantitative analysis of intracellular binding-diffusion kinetics, but how to dependably extract accurate kinetic constants from such analyses is still an open question. To this end, in this study, we developed what we believe is a new analytical model for confocal FRAP-based measurements of intracellular binding-diffusion processes, based on a closed-form equation of the FRAP formula for a spot photobleach geometry. This approach incorporates a binding diffusion model that allows for diffusion of both the unbound and bound species, and also compensates for binding diffusion that occurs during photobleaching, a critical consideration in confocal FRAP analysis. In addition, to address the problem of parametric multiplicity, we propose a scheme to reduce the number of fitting parameters in the effective diffusion subregime when D's for the bound and unbound species are known. We validate this method by measuring kinetic rate constants for the CAAX-mediated binding of Ras to membranes of the endoplasmic reticulum, obtaining binding constants of k(on) ∼ 255/s and k(off) ∼ 31/s.

Studying Smad2 intranuclear diffusion dynamics by mathematical modelling of FRAP experiments

We combine Fluorescence Recovery After Photobleaching (FRAP) experiments with mathematical modelling to study the dynamics inside the nucleus of both the TGF-β-sensitive transcriptional regulator Smad2, and Green-Fluorescent Protein (GFP). We show how combining modelling with bleaching strips of different areas allows a rigorous test of whether or not a protein is moving via diffusion as a single species. As noted recently by others, it is important to consider diffusion during the bleaching process. Neglecting it can cause serious error. Also, it is possible to use the bleaching process itself to provide an extra consistency test to the models predicting the recovery. With our method we show that the dynamics of GFP are consistent with it diffusing as a single species in a uniform environment in which flow is negligible. In contrast, the dynamics of the intracellular signal transducer Smad2 are never consistent with it moving as a single species via simple diffusion in a homogeneous environment without flow. Adding TGF-β slows down the dynamics of Smad2 but even without TGF-β, the Smad2 dynamics are influenced by one or more of: association, flow, and inhomogeneity in space of the dynamics. We suggest that the dynamics inside cells of many proteins may be poorly described by simple diffusion of a single species, and that our methodology provides a general and powerful way to test this hypothesis.

A synergic simulation-optimization approach for analyzing biomolecular dynamics in living organisms

A synergic duo simulation-optimization approach was developed and implemented to study protein-substrate dynamics and binding kinetics in living organisms. The forward problem is a system of several coupled nonlinear partial differential equations which, with a given set of kinetics and diffusion parameters, can provide not only the commonly used bleached area-averaged time series in fluorescence microscopy experiments but more informative full biomolecular/drug space-time series and can be successfully used to study dynamics of both Dirac and Gaussian fluorescence-labeled biomacromolecules in vivo. The incomplete Cholesky preconditioner was coupled with the finite difference discretization scheme and an adaptive time-stepping strategy to solve the forward problem. The proposed approach was validated with analytical as well as reference solutions and used to simulate dynamics of GFP-tagged glucocorticoid receptor (GFP-GR) in mouse cancer cell during a fluorescence recovery after photobleaching experiment. Model analysis indicates that the commonly practiced bleach spot-averaged time series is not an efficient approach to extract physiological information from the fluorescence microscopy protocols. It was recommended that experimental biophysicists should use full space-time series, resulting from experimental protocols, to study dynamics of biomacromolecules and drugs in living organisms. It was also concluded that in parameterization of biological mass transfer processes, setting the norm of the gradient of the penalty function at the solution to zero is not an efficient stopping rule to end the inverse algorithm. Theoreticians should use multi-criteria stopping rules to quantify model parameters by optimization.

Modeling capping protein FRAP and CALI experiments reveals in vivo regulation of actin dynamics

To gain insights on cellular mechanisms regulating actin polymerization, we used the Virtual Cell to model fluorescence recovery after photobleaching (FRAP) and chromophore-assisted laser inactivation (CALI) experiments on EGFP-capping protein (EGFP-CP). Modeling the FRAP kinetics demonstrated that the in vivo rate for the dissociation of CP from actin filaments is much faster (approximately 0.1 s(-1)) than that measured in vitro (0.01-0.0004 s(-1)). The CALI simulation revealed that in order to induce sustainable changes in cell morphology after CP inactivation, the cells should exhibit anticapping ability. We included the VASP protein as the anticapping agent in the modeling scheme. The model predicts that VASP affinity for barbed ends has a cooperative dependence on the concentration of VASP-barbed end complexes. This dependence produces a positive feedback that stabilizes the complexes and allows sustained growth at clustered filament tips. We analyzed the range of laser intensities that are sufficient to induce changes in cell morphology. This analysis demonstrates that FRAP experiments with EGFP-CP can be performed safely without changes in cell morphology, because, the intensity of the photobleaching beam is not high enough to produce the critical concentration of free barbed ends that will induce filament growth before diffusional replacement of EGFP-CP occurs.

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.