The role of time delay in adaptive cellular negative feedback systems.

Adaptation in cellular systems is often mediated by negative feedbacks, which usually come with certain time delays causing several characteristic response patterns including an overdamped response, damped or sustained oscillations. Here, we analyse generic two-dimensional delay differential equations with delayed negative feedback describing the dynamics of biochemical adaptive signal-response networks. We derive explicit thresholds and boundaries showing how time delay determines characteristic response patterns of these networks. Applying our theoretical analyses to concrete data we show that adaptation to osmotic stress in yeast is optimal in the sense of minimizing adaptation time without causing oscillatory behaviour, i.e., a critically damped response. In addition, our framework demonstrates that a slight increase of time delay in the NF-κB system might induce a switch from damped to sustained oscillatory behaviour. Thus, we demonstrate how delay differential equations can be used to explicitly study the delay in biochemical negative feedback systems. Our analysis also provides insight into how time delay may tune biological signal-response patterns and control the systems behaviour.

FoCo: a simple and robust quantification algorithm of nuclear foci.

BACKGROUND: The number of γH2AX foci per nucleus is an accepted measure of the number of DNA double-strand breaks in single cells. One of the experimental techniques for γH2AX detection in cultured cells is immunofluorescent labelling of γH2AX and nuclei followed by microscopy imaging and analysis. RESULTS: In this study, we present the algorithm FoCo for reliable and robust automatic nuclear foci counting in single cell images. FoCo has the following advantages with respect to other software packages: i) the ability to reliably quantify even densely distributed foci, e.g., on images of cells subjected to radiation doses up to 10 Gy, ii) robustness of foci quantification in the sense of suppressing out-of-focus background signal, and iii) its simplicity. FoCo requires only 5 parameters that have to be adjusted by the user. CONCLUSIONS: FoCo is an open-source user-friendly software with GUI for individual foci counting, which is able to produce reliable and robust foci quantifications even for low signal/noise ratios and densely distributed foci.

To senesce or not to senesce: how primary human fibroblasts decide their cell fate after DNA damage.

Excessive DNA damage can induce an irreversible cell cycle arrest, called senescence, which is generally perceived as an important tumour-suppressor mechanism. However, it is unclear how cells decide whether to senesce or not after DNA damage. By combining experimental data with a parameterized mathematical model we elucidate this cell fate decision at the G1-S transition. Our model provides a quantitative and conceptually new understanding of how human fibroblasts decide whether DNA damage is beyond repair and senesce. Model and data imply that the G1-S transition is regulated by a bistable hysteresis switch with respect to Cdk2 activity, which in turn is controlled by the Cdk2/p21 ratio rather than cyclin abundance. We experimentally confirm the resulting predictions that to induce senescence i) in healthy cells both high initial and elevated background DNA damage are necessary and sufficient, and ii) in already damaged cells much lower additional DNA damage is sufficient. Our study provides a mechanistic explanation of a) how noise in protein abundances allows cells to overcome the G1-S arrest even with substantial DNA damage, potentially leading to neoplasia, and b) how accumulating DNA damage with age increasingly sensitizes cells for senescence.

Development of a robust DNA damage model including persistent telomere-associated damage with application to secondary cancer risk assessment.

Mathematical modelling has been instrumental to understand kinetics of radiation-induced DNA damage repair and associated secondary cancer risk. The widely accepted two-lesion kinetic (TLK) model assumes two kinds of double strand breaks, simple and complex ones, with different repair rates. Recently, persistent DNA damage associated with telomeres was reported as a new kind of DNA damage. We therefore extended existing versions of the TLK model by new categories of DNA damage and re-evaluated those models using extensive data. We subjected different versions of the TLK model to a rigorous model discrimination approach. This enabled us to robustly select a best approximating parsimonious model that can both recapitulate and predict transient and persistent DNA damage after ionizing radiation. Models and data argue for i) nonlinear dose-damage relationships, and ii) negligible saturation of repair kinetics even for high doses. Additionally, we show that simulated radiation-induced persistent telomere-associated DNA damage foci (TAF) can be used to predict excess relative risk (ERR) of developing secondary leukemia after fractionated radiotherapy. We suggest that TAF may serve as an additional measure to predict cancer risk after radiotherapy using high dose rates. This may improve predicting risk-dose dependency of ionizing radiation especially for long-term therapies.

Nested autoinhibitory feedbacks alter the resistance of homeostatic adaptive biochemical networks.

Negative feedback control is a ubiquitous feature of biochemical systems, as is time delay between a signal and its response. Negative feedback in conjunction with time delay can lead to oscillations. In a cellular context, it might be beneficial to mitigate oscillatory behaviour to avoid recurring stress situations. This can be achieved by increasing the distance between the parameters of the system and certain thresholds, beyond which oscillations occur. This distance has been termed resistance. Here, we prove that in a generic three-dimensional negative feedback system the resistance of the system is modified by nested autoinhibitory feedbacks. Our system features negative feedbacks through both input-inhibition as well as output-activation, a signalling component with mass conservation and perfect adaptation. We show that these features render the system applicable to biological data, exemplified by the high osmolarity glycerol system in yeast and the mammalian p53 system. Output-activation is better supported by data than input-inhibition and also shows distinguished properties with respect to the system's stimulus. Our general approach might be useful in designing synthetic systems in which oscillations can be tuned by synthetic autoinhibitory feedbacks.