Natural killer cell immune synapse formation and cytotoxicity are controlled by tension of the target interface.

Natural killer (NK) cells can kill infected or transformed cells via a lytic immune synapse. Diseased cells may exhibit altered mechanical properties but how this impacts NK cell responsiveness is unknown. We report that human NK cells were stimulated more effectively to secrete granzymes A and B, FasL (also known as FasLG), granulysin and IFNγ, by stiff (142 kPa) compared to soft (1 kPa) planar substrates. To create surrogate spherical targets of defined stiffness, sodium alginate was used to synthesise soft (9 kPa), medium (34 kPa) or stiff (254 kPa) cell-sized beads, coated with antibodies against activating receptor NKp30 (also known as NCR3) and the integrin LFA-1 (also known as ITGAL). Against stiff beads, NK cells showed increased degranulation. Polarisation of the microtubule-organising centre and lytic granules were impaired against soft targets, which instead resulted in the formation of unstable kinapses. Thus, by varying target stiffness to characterise the mechanosensitivity of immune synapses, we identify soft targets as a blind spot in NK cell recognition. This article has an associated First Person interview with the co-first authors of the paper.

Multiplexing information flow through dynamic signalling systems.

We consider how a signalling system can act as an information hub by multiplexing information arising from multiple signals. We formally define multiplexing, mathematically characterise which systems can multiplex and how well they can do it. While the results of this paper are theoretical, to motivate the idea of multiplexing, we provide experimental evidence that tentatively suggests that the NF-κB transcription factor can multiplex information about changes in multiple signals. We believe that our theoretical results may resolve the apparent paradox of how a system like NF-κB that regulates cell fate and inflammatory signalling in response to diverse stimuli can appear to have the low information carrying capacity suggested by recent studies on scalar signals. In carrying out our study, we introduce new methods for the analysis of large, nonlinear stochastic dynamic models, and develop computational algorithms that facilitate the calculation of fundamental constructs of information theory such as Kullback-Leibler divergences and sensitivity matrices, and link these methods to a new theory about multiplexing information. We show that many current models such as those of the NF-κB system cannot multiplex effectively and provide models that overcome this limitation using post-transcriptional modifications.

Disentangling juxtacrine from paracrine signalling in dynamic tissue.

Prolactin is a major hormone product of the pituitary gland, the central endocrine regulator. Despite its physiological importance, the cell-level mechanisms of prolactin production are not well understood. Having significantly improved the resolution of real-time-single-cell-GFP-imaging, the authors recently revealed that prolactin gene transcription is highly dynamic and stochastic yet shows space-time coordination in an intact tissue slice. However, it still remains an open question as to what kind of cellular communication mediates the observed space-time organization. To determine the type of interaction between cells we developed a statistical model. The degree of similarity between two expression time series was studied in terms of two distance measures, Euclidean and geodesic, the latter being a network-theoretic distance defined to be the minimal number of edges between nodes, and this was used to discriminate between juxtacrine from paracrine signalling. The analysis presented here suggests that juxtacrine signalling dominates. To further determine whether the coupling is coordinating transcription or post-transcriptional activities we used stochastic switch modelling to infer the transcriptional profiles of cells and estimated their similarity measures to deduce that their spatial cellular coordination involves coupling of transcription via juxtacrine signalling. We developed a computational model that involves an inter-cell juxtacrine coupling, yielding simulation results that show space-time coordination in the transcription level that is in agreement with the above analysis. The developed model is expected to serve as the prototype for the further study of tissue-level organised gene expression for epigenetically regulated genes, such as prolactin.

Quantitative analysis reveals crosstalk mechanisms of heat shock-induced attenuation of NF-κB signaling at the single cell level.

Elevated temperature induces the heat shock (HS) response, which modulates cell proliferation, apoptosis, the immune and inflammatory responses. However, specific mechanisms linking the HS response pathways to major cellular signaling systems are not fully understood. Here we used integrated computational and experimental approaches to quantitatively analyze the crosstalk mechanisms between the HS-response and a master regulator of inflammation, cell proliferation, and apoptosis the Nuclear Factor κB (NF-κB) system. We found that populations of human osteosarcoma cells, exposed to a clinically relevant 43°C HS had an attenuated NF-κB p65 response to Tumor Necrosis Factor α (TNFα) treatment. The degree of inhibition of the NF-κB response depended on the HS exposure time. Mathematical modeling of single cells indicated that individual crosstalk mechanisms differentially encode HS-mediated NF-κB responses while being consistent with the observed population-level responses. In particular "all-or-nothing" encoding mechanisms were involved in the HS-dependent regulation of the IKK activity and IκBα phosphorylation, while others involving transport were "analogue". In order to discriminate between these mechanisms, we used live-cell imaging of nuclear translocations of the NF-κB p65 subunit. The single cell responses exhibited "all-or-nothing" encoding. While most cells did not respond to TNFα stimulation after a 60 min HS, 27% showed responses similar to those not receiving HS. We further demonstrated experimentally and theoretically that the predicted inhibition of IKK activity was consistent with the observed HS-dependent depletion of the IKKα and IKKß subunits in whole cell lysates. However, a combination of "all-or-nothing" crosstalk mechanisms was required to completely recapitulate the single cell data. We postulate therefore that the heterogeneity of the single cell responses might be explained by the cell-intrinsic variability of HS-modulated IKK signaling. In summary, we show that high temperature modulates NF-κB responses in single cells in a complex and unintuitive manner, which needs to be considered in hyperthermia-based treatment strategies.

Glucocorticoid receptor regulates accurate chromosome segregation and is associated with malignancy.

The glucocorticoid receptor (GR) is a member of the nuclear receptor superfamily, which controls programs regulating cell proliferation, differentiation, and apoptosis. We have identified an unexpected role for GR in mitosis. We discovered that specifically modified GR species accumulate at the mitotic spindle during mitosis in a distribution that overlaps with Aurora kinases. We found that Aurora A was required to mediate mitosis-driven GR phosphorylation, but not recruitment of GR to the spindle. GR was necessary for mitotic progression, with increased time to complete mitosis, frequency of mitotic aberrations, and death in mitosis observed following GR knockdown. Complementation studies revealed an essential role for the GR ligand-binding domain, but no clear requirement for ligand binding in regulating chromosome segregation. The GR N-terminal domain, and specifically phosphosites S203 and S211, were not required. Reduced GR expression results in a cell cycle phenotype, with isolated cells from mouse and human subjects showing changes in chromosome content over prolonged passage. Furthermore, GR haploinsufficient mice have an increased incidence of tumor formation, and, strikingly, these tumors are further depleted for GR, implying additional GR loss as a consequence of cell transformation. We identified reduced GR expression in a panel of human liver, lung, prostate, colon, and breast cancers. We therefore reveal an unexpected role for the GR in promoting accurate chromosome segregation during mitosis, which is causally linked to tumorigenesis, making GR an authentic tumor suppressor gene.

p63 is required beside p53 for PERP-mediated apoptosis in uveal melanoma.

BACKGROUND: PERP (p53 apoptosis effector related to PMP-22), a transcriptional target of p53, is downregulated and contributes to the impairment of apoptosis in uveal melanoma (UM). Intriguingly, PERP is not induced in UM despite functional p53. p63, located on chromosome 3, which is characteristically altered in high-risk UM, can transactivate PERP. Here, we determine the functional role of p63 expression in the initiation of p53/PERP-mediated apoptosis in UM. METHODS: PERP expression was monitored by quantitative PCR (qPCR) and immunoblotting in UM cell lines treated with DNA-damaging agents. The functional role of p63 was assessed by transient expression of p63-turbo GFP (p63-tGFP) in the apoptosis- resistant, 3q-deficient OCM-1 cells. Expression and localisation of p63, PERP and p53, and induction of apoptosis were characterised by qPCR, immunoblotting and live cell confocal microscopy. RESULTS: PERP expression was significantly downregulated in all UM cell lines. DNA-damaging treatments failed to induce apoptosis and activate PERP in OCM-1 cells, which displayed non-functional levels of p63. Expression of p63-tGFP induced apoptosis with marked increase in PERP expression and associated p53 accumulation. CONCLUSIONS: Lack of p63 contributes to reduced PERP levels and impaired p53-mediated apoptosis in UM. p63 expression is required for PERP-mediated apoptosis in UM.

A stochastic transcriptional switch model for single cell imaging data.

Gene expression is made up of inherently stochastic processes within single cells and can be modeled through stochastic reaction networks (SRNs). In particular, SRNs capture the features of intrinsic variability arising from intracellular biochemical processes. We extend current models for gene expression to allow the transcriptional process within an SRN to follow a random step or switch function which may be estimated using reversible jump Markov chain Monte Carlo (MCMC). This stochastic switch model provides a generic framework to capture many different dynamic features observed in single cell gene expression. Inference for such SRNs is challenging due to the intractability of the transition densities. We derive a model-specific birth-death approximation and study its use for inference in comparison with the linear noise approximation where both approximations are considered within the unifying framework of state-space models. The methodology is applied to synthetic as well as experimental single cell imaging data measuring expression of the human prolactin gene in pituitary cells.

Measurement of single-cell dynamics.

Populations of cells are almost always heterogeneous in function and fate. To understand the plasticity of cells, it is vital to measure quantitatively and dynamically the molecular processes that underlie cell-fate decisions in single cells. Early events in cell signalling often occur within seconds of the stimulus, whereas intracellular signalling processes and transcriptional changes can take minutes or hours. By contrast, cell-fate decisions, such as whether a cell divides, differentiates or dies, can take many hours or days. Multiparameter experimental and computational methods that integrate quantitative measurement and mathematical simulation of these noisy and complex processes are required to understand the highly dynamic mechanisms that control cell plasticity and fate.

Tight control of hypoxia-inducible factor-α transient dynamics is essential for cell survival in hypoxia.

Intracellular signaling involving hypoxia-inducible factor (HIF) controls the adaptive responses to hypoxia. There is a growing body of evidence demonstrating that intracellular signals encode temporal information. Thus, the dynamics of protein levels, as well as protein quantity and/or localization, impacts on cell fate. We hypothesized that such temporal encoding has a role in HIF signaling and cell fate decisions triggered by hypoxic conditions. Using live cell imaging in a controlled oxygen environment, we observed transient 3-h pulses of HIF-1α and -2α expression under continuous hypoxia. We postulated that the well described prolyl hydroxylase (PHD) oxygen sensors and HIF negative feedback regulators could be the origin of the pulsatile HIF dynamics. We used iterative mathematical modeling and experimental analysis to scrutinize which parameter of the PHD feedback could control HIF timing and we probed for the functional redundancy between the three main PHD proteins. We identified PHD2 as the main PHD responsible for HIF peak duration. We then demonstrated that this has important consequences, because the transient nature of the HIF pulse prevents cell death by avoiding transcription of p53-dependent pro-apoptotic genes. We have further shown the importance of considering HIF dynamics for coupling mathematical models by using a described HIF-p53 mathematical model. Our results indicate that the tight control of HIF transient dynamics has important functional consequences on the cross-talk with key signaling pathways controlling cell survival, which is likely to impact on HIF targeting strategies for hypoxia-associated diseases such as tumor progression and ischemia.

A method of 'speed coefficients' for biochemical model reduction applied to the NF-κB system.

The relationship between components of biochemical network and the resulting dynamics of the overall system is a key focus of computational biology. However, as these networks and resulting mathematical models are inherently complex and non-linear, the understanding of this relationship becomes challenging. Among many approaches, model reduction methods provide an avenue to extract components responsible for the key dynamical features of the system. Unfortunately, these approaches often require intuition to apply. In this manuscript we propose a practical algorithm for the reduction of biochemical reaction systems using fast-slow asymptotics. This method allows the ranking of system variables according to how quickly they approach their momentary steady state, thus selecting the fastest for a steady state approximation. We applied this method to derive models of the Nuclear Factor kappa B network, a key regulator of the immune response that exhibits oscillatory dynamics. Analyses with respect to two specific solutions, which corresponded to different experimental conditions identified different components of the system that were responsible for the respective dynamics. This is an important demonstration of how reduction methods that provide approximations around a specific steady state, could be utilised in order to gain a better understanding of network topology in a broader context.

Dynamic analysis of stochastic transcription cycles.

In individual mammalian cells the expression of some genes such as prolactin is highly variable over time and has been suggested to occur in stochastic pulses. To investigate the origins of this behavior and to understand its functional relevance, we quantitatively analyzed this variability using new mathematical tools that allowed us to reconstruct dynamic transcription rates of different reporter genes controlled by identical promoters in the same living cell. Quantitative microscopic analysis of two reporter genes, firefly luciferase and destabilized EGFP, was used to analyze the dynamics of prolactin promoter-directed gene expression in living individual clonal and primary pituitary cells over periods of up to 25 h. We quantified the time-dependence and cyclicity of the transcription pulses and estimated the length and variation of active and inactive transcription phases. We showed an average cycle period of approximately 11 h and demonstrated that while the measured time distribution of active phases agreed with commonly accepted models of transcription, the inactive phases were differently distributed and showed strong memory, with a refractory period of transcriptional inactivation close to 3 h. Cycles in transcription occurred at two distinct prolactin-promoter controlled reporter genes in the same individual clonal or primary cells. However, the timing of the cycles was independent and out-of-phase. For the first time, we have analyzed transcription dynamics from two equivalent loci in real-time in single cells. In unstimulated conditions, cells showed independent transcription dynamics at each locus. A key result from these analyses was the evidence for a minimum refractory period in the inactive-phase of transcription. The response to acute signals and the result of manipulation of histone acetylation was consistent with the hypothesis that this refractory period corresponded to a phase of chromatin remodeling which significantly increased the cyclicity. Stochastically timed bursts of transcription in an apparently random subset of cells in a tissue may thus produce an overall coordinated but heterogeneous phenotype capable of acute responses to stimuli.

Cellular heterogeneity in TNF/TNFR1 signalling: live cell imaging of cell fate decisions in single cells.

Cellular responses to TNF are inherently heterogeneous within an isogenic cell population and across different cell types. TNF promotes cell survival by activating pro-inflammatory NF-κB and MAPK signalling pathways but may also trigger apoptosis and necroptosis. Following TNF stimulation, the fate of individual cells is governed by the balance of pro-survival and pro-apoptotic signalling pathways. To elucidate the molecular mechanisms driving heterogenous responses to TNF, quantifying TNF/TNFR1 signalling at the single-cell level is crucial. Fluorescence live-cell imaging techniques offer real-time, dynamic insights into molecular processes in single cells, allowing for detection of rapid and transient changes, as well as identification of subpopulations, that are likely to be missed with traditional endpoint assays. Whilst fluorescence live-cell imaging has been employed extensively to investigate TNF-induced inflammation and TNF-induced cell death, it has been underutilised in studying the role of TNF/TNFR1 signalling pathway crosstalk in guiding cell-fate decisions in single cells. Here, we outline the various opportunities for pathway crosstalk during TNF/TNFR1 signalling and how these interactions may govern heterogenous responses to TNF. We also advocate for the use of live-cell imaging techniques to elucidate the molecular processes driving cell-to-cell variability in single cells. Understanding and overcoming cellular heterogeneity in response to TNF and modulators of the TNF/TNFR1 signalling pathway could lead to the development of targeted therapies for various diseases associated with aberrant TNF/TNFR1 signalling, such as rheumatoid arthritis, metabolic syndrome, and cancer.

Unregulated actin polymerization by WASp causes defects of mitosis and cytokinesis in X-linked neutropenia.

Specific mutations in the human gene encoding the Wiskott-Aldrich syndrome protein (WASp) that compromise normal auto-inhibition of WASp result in unregulated activation of the actin-related protein 2/3 complex and increased actin polymerizing activity. These activating mutations are associated with an X-linked form of neutropenia with an intrinsic failure of myelopoiesis and an increase in the incidence of cytogenetic abnormalities. To study the underlying mechanisms, active mutant WASp(I294T) was expressed by gene transfer. This caused enhanced and delocalized actin polymerization throughout the cell, decreased proliferation, and increased apoptosis. Cells became binucleated, suggesting a failure of cytokinesis, and micronuclei were formed, indicative of genomic instability. Live cell imaging demonstrated a delay in mitosis from prometaphase to anaphase and confirmed that multinucleation was a result of aborted cytokinesis. During mitosis, filamentous actin was abnormally localized around the spindle and chromosomes throughout their alignment and separation, and it accumulated within the cleavage furrow around the spindle midzone. These findings reveal a novel mechanism for inhibition of myelopoiesis through defective mitosis and cytokinesis due to hyperactivation and mislocalization of actin polymerization.

Pulsatile patterns of pituitary hormone gene expression change during development.

Important questions in biology have emerged recently concerning the timing of transcription in living cells. Studies on clonal cell lines have shown that transcription is often pulsatile and stochastic, with implications for cellular differentiation. Currently, information regarding transcriptional activity at cellular resolution within a physiological context remains limited. To investigate single-cell transcriptional activity in real-time in living tissue we used bioluminescence imaging of pituitary tissue from transgenic rats in which luciferase gene expression is driven by a pituitary hormone gene promoter. We studied fetal and neonatal pituitary tissue to assess whether dynamic patterns of transcription change during tissue development. We show that gene expression in single cells is highly pulsatile at the time endocrine cells first appear but becomes stabilised as the tissue develops in early neonatal life. This stabilised transcription pattern might depend upon tissue architecture or paracrine signalling, as isolated cells, generated from enzymatic dispersion of the tissue, display pulsatile luminescence. Nascent cells in embryonic tissue also showed coordinated transcription activity over short distances further indicating that cellular context is important for transcription activity. Overall, our data show that cells alter their patterns of gene expression according to their context and developmental stage, with important implications for cellular differentiation.

Population robustness arising from cellular heterogeneity.

Heterogeneity between individual cells is a common feature of dynamic cellular processes, including signaling, transcription, and cell fate; yet the overall tissue level physiological phenotype needs to be carefully controlled to avoid fluctuations. Here we show that in the NF-kappaB signaling system, the precise timing of a dual-delayed negative feedback motif [involving stochastic transcription of inhibitor kappaB (IkappaB)-alpha and -epsilon] is optimized to induce heterogeneous timing of NF-kappaB oscillations between individual cells. We suggest that this dual-delayed negative feedback motif enables NF-kappaB signaling to generate robust single cell oscillations by reducing sensitivity to key parameter perturbations. Simultaneously, enhanced cell heterogeneity may represent a mechanism that controls the overall coordination and stability of cell population responses by decreasing temporal fluctuations of paracrine signaling. It has often been thought that dynamic biological systems may have evolved to maximize robustness through cell-to-cell coordination and homogeneity. Our analyses suggest in contrast, that this cellular variation might be advantageous and subject to evolutionary selection. Alternative types of therapy could perhaps be designed to modulate this cellular heterogeneity.

Overexpression of IκB⍺ modulates NF-κB activation of inflammatory target gene expression.

Cells respond to inflammatory stimuli such as cytokines by activation of the nuclear factor-κB (NF-κB) signalling pathway, resulting in oscillatory translocation of the transcription factor p65 between nucleus and cytoplasm in some cell types. We investigate the relationship between p65 and inhibitor-κB⍺ (IκBα) protein levels and dynamic properties of the system, and how this interaction impacts on the expression of key inflammatory genes. Using bacterial artificial chromosomes, we developed new cell models of IκB⍺-eGFP protein overexpression in a pseudo-native genomic context. We find that cells with high levels of the negative regulator IκBα remain responsive to inflammatory stimuli and maintain dynamics for both p65 and IκBα. In contrast, canonical target gene expression is dramatically reduced by overexpression of IκBα, but can be partially rescued by overexpression of p65. Treatment with leptomycin B to promote nuclear accumulation of IκB⍺ also suppresses canonical target gene expression, suggesting a mechanism in which nuclear IκB⍺ accumulation prevents productive p65 interaction with promoter binding sites. This causes reduced target promoter binding and gene transcription, which we validate by chromatin immunoprecipitation and in primary cells. Overall, we show how inflammatory gene transcription is modulated by the expression levels of both IκB⍺ and p65. This results in an anti-inflammatory effect on transcription, demonstrating a broad mechanism to modulate the strength of inflammatory response.

Dynamic organisation of prolactin gene expression in living pituitary tissue.

Gene expression in living cells is highly dynamic, but temporal patterns of gene expression in intact tissues are largely unknown. The mammalian pituitary gland comprises several intermingled cell types, organised as interdigitated networks that interact functionally to generate co-ordinated hormone secretion. Live-cell imaging was used to quantify patterns of reporter gene expression in dispersed lactotrophic cells or intact pituitary tissue from bacterial artificial chromosome (BAC) transgenic rats in which a large prolactin genomic fragment directed expression of luciferase or destabilised enhanced green fluorescent protein (d2EGFP). Prolactin promoter activity in transgenic pituitaries varied with time across different regions of the gland. Although amplitude of transcriptional responses differed, all regions of the gland displayed similar overall patterns of reporter gene expression over a 50-hour period, implying overall co-ordination of cellular behaviour. By contrast, enzymatically dispersed pituitary cell cultures showed unsynchronised fluctuations of promoter activity amongst different cells, suggesting that transcriptional patterns were constrained by tissue architecture. Short-term, high resolution, single cell analyses in prolactin-d2EGFP transgenic pituitary slice preparations showed varying transcriptional patterns with little correlation between adjacent cells. Together, these data suggest that pituitary tissue comprises a series of cell ensembles, which individually display a variety of patterns of short-term stochastic behaviour, but together yield long-range and long-term coordinated behaviour.

Physiological levels of TNFalpha stimulation induce stochastic dynamics of NF-kappaB responses in single living cells.

Nuclear factor kappa B (NF-kappaB) signalling is activated by cellular stress and inflammation and regulates cytokine expression. We applied single-cell imaging to investigate dynamic responses to different doses of tumour necrosis factor alpha (TNFalpha). Lower doses activated fewer cells and those responding showed an increasingly variable delay in the initial NF-kappaB nuclear translocation and associated IkappaBalpha degradation. Robust 100 minute nuclear:cytoplasmic NF-kappaB oscillations were observed over a wide range of TNFalpha concentrations. The result is supported by computational analyses, which identified a limit cycle in the system with a stable 100 minute period over a range of stimuli, and indicated no co-operativity in the pathway activation. These results suggest that a stochastic threshold controls functional all-or-nothing responses in individual cells. Deterministic and stochastic models simulated the experimentally observed activation threshold and gave rise to new predictions about the structure of the system and open the way for better mechanistic understanding of physiological TNFalpha activation of inflammatory responses in cells and tissues.

A systematic survey of the response of a model NF-κB signalling pathway to TNFα stimulation.

White's lab established that strong, continuous stimulation with tumour necrosis factor-α (TNFα) can induce sustained oscillations in the subcellular localisation of the transcription factor nuclear factor κB (NF-κB). But the intensity of the TNFα signal varies substantially, from picomolar in the blood plasma of healthy organisms to nanomolar in diseased states. We report on a systematic survey using computational bifurcation theory to explore the relationship between the intensity of TNFα stimulation and the existence of sustained NF-κB oscillations. Using a deterministic model developed by Ashall et al. in 2009, we find that the system's responses to TNFα are characterised by a supercritical Hopf bifurcation point: above a critical intensity of TNFα the system exhibits sustained oscillations in NF-kB localisation. For TNFα below this critical value, damped oscillations are observed. This picture depends, however, on the values of the model's other parameters. When the values of certain reaction rates are altered the response of the signalling pathway to TNFα stimulation changes: in addition to the sustained oscillations induced by high-dose stimulation, a second oscillatory regime appears at much lower doses. Finally, we define scores to quantify the sensitivity of the dynamics of the system to variation in its parameters and use these scores to establish that the qualitative dynamics are most sensitive to the details of NF-κB mediated gene transcription.

Apoptotic priming is defined by the dynamic exchange of Bcl-2 proteins between mitochondria and cytosol.

Apoptosis is regulated by interactions between the BH3-only and multi-domain Bcl-2 family proteins. These interactions are integrated on the outer mitochondrial membrane (OMM) where they set the threshold for apoptosis, known as mitochondrial priming. However, how mitochondrial priming is controlled at the level of single cells remains unclear. Retrotranslocation of Bcl-XL has been proposed as one mechanism, removing pro-apoptotic Bcl-2 proteins from the OMM, thus reducing priming. Contrary to this view, we now show that Bcl-XL retrotranslocation is inhibited by binding to its BH3-only partners, resulting in accumulation of these protein complexes on mitochondria. We find that Bcl-XL retrotranslocation dynamics are tightly coupled to mitochondrial priming. Quantifying these dynamics indicates the heterogeneity in priming between cells within a population and predicts how they subsequently respond to a pro-apoptotic signal.