Minimum criteria for DNA damage-induced phase advances in circadian rhythms

A coupled model between circadian, cell-cycle, and redox rhythms reveals their regulation of oxidative stress

Most organisms possess three biological oscillators, circadian clock, cell cycle, and redox rhythm, which are autonomous but interact each other. However, whether their interactions and autonomy are beneficial for organisms remains unclear. Here, we modeled a coupled oscillator system where each oscillator affected the phase of the other oscillators. We found that multiple types of coupling prevent a high H2O2 level in cells at M phase. Consequently, we hypothesized a high H2O2 sensitivity at the M phase and found that moderate coupling reduced cell damage due to oxidative stress by generating appropriate phase relationships between three rhythms, whereas strong coupling resulted in an elevated cell damage by increasing the average H2O2 level and disrupted the cell cycle. Furthermore, the multicellularity model revealed that phase variations among cells confer flexibility in synchronization with environments at the expense of adaptability to the optimal environment. Thus, both autonomy and synchrony among the oscillators are important for coordinating their phase relationships to minimize oxidative stress, and couplings balance them depending on environments.

Bilateral Feedback in Oscillator Model Is Required to Explain the Coupling Dynamics of Hes1 with the Cell Cycle

Biological processes are governed by the expression of proteins, and for some proteins, their level of expression can fluctuate periodically over time (i.e., they oscillate). Many oscillatory proteins (e.g., cell cycle proteins and those from the HES family of transcription factors) are connected in complex ways, often within large networks. This complexity can be elucidated by developing intuitive mathematical models that describe the underlying critical aspects of the relationships between these processes. Here, we provide a mathematical explanation of a recently discovered biological phenomenon: the phasic position of the gene Hes1’s oscillatory expression at the beginning of the cell cycle of an individual human breast cancer stem cell can have a predictive value on how long that cell will take to complete a cell cycle. We use a two-component model of coupled oscillators to represent Hes1 and the cell cycle in the same cell with minimal assumptions. Inputting only the initial phase angles, we show that this model is capable of predicting the dynamic mitosis to mitosis behaviour of Hes1 and predicting cell cycle length patterns as found in real-world experimental data. Moreover, we discover that bidirectional coupling between Hes1 and the cell cycle is critical within the system for the data to be reproduced and that nonfixed asymmetry in the interactions between the oscillators is required. The phase dynamics we present here capture the complex interplay between Hes1 and the cell cycle, helping to explain nongenetic cell cycle variability, which has critical implications in cancer treatment contexts.

The Circadian Gene Per1 Plays an Important Role in Radiation-Induced Apoptosis and DNA Damage in Glioma

Objective:

Period1 (PER1), a core circadian gene, not only modulates circadian rhythm but may also play an important role in other biological processes, including pathways involved in the proliferation and apoptosis of tumor cells. In this study, we investigated the mechanism by which the downregulated expression of PER1 promotes the apoptosis of wild-type P53 human glioma U343 cells exposed to X-rays.

Methods:

U343 cells were exposed to 6 mV 10 Gy X-ray irradiation after infection with an shRNA lentivirus to reduce the expression of PER1 and were analyzed by SCGE analysis, flow cytometry, qRT-PCR, and western blotting.

Result:

SCGE analysis revealed that compared with the controls, U343 cells expressing low levels of PER1 showed minor DNA damage when exposed to X-ray irradiation (P<0.05), and the flow cytometry assay showed lower death rates (P<0.05). RT-PCR and western blot analysis both revealed decreased expression of CHK2 and P53, which regulate DNA damage and repair via the CHK2-P53 pathway, and decreased expression of C-MYC, which is related to cell apoptosis.

Conclusion:

Our research suggests that PER1 may play an important role in tumor radiotherapy, which is attributable to enhanced chk2-P53 signaling and proapoptotic processes. These findings provide a new target for the clinical treatment of glioma and a reliable basis for postradiation therapy and gene therapy for glioma and other cancers.

Systems Chronotherapeutics

Chronotherapeutics aim at treating illnesses according to the endogenous biologic rhythms, which moderate xenobiotic metabolism and cellular drug response. The molecular clocks present in individual cells involve approximately fifteen clock genes interconnected in regulatory feedback loops. They are coordinated by the suprachiasmatic nuclei, a hypothalamic pacemaker, which also adjusts the circadian rhythms to environmental cycles. As a result, many mechanisms of diseases and drug effects are controlled by the circadian timing system. Thus, the tolerability of nearly 500 medications varies by up to fivefold according to circadian scheduling, both in experimental models and/or patients. Moreover, treatment itself disrupted, maintained, or improved the circadian timing system as a function of drug timing. Improved patient outcomes on circadian-based treatments (chronotherapy) have been demonstrated in randomized clinical trials, especially for cancer and inflammatory diseases. However, recent technological advances have highlighted large interpatient differences in circadian functions resulting in significant variability in chronotherapy response. Such findings advocate for the advancement of personalized chronotherapeutics through interdisciplinary systems approaches. Thus, the combination of mathematical, statistical, technological, experimental, and clinical expertise is now shaping the development of dedicated devices and diagnostic and delivery algorithms enabling treatment individualization. In particular, multiscale systems chronopharmacology approaches currently combine mathematical modeling based on cellular and whole-body physiology to preclinical and clinical investigations toward the design of patient-tailored chronotherapies. We review recent systems research works aiming to the individualization of disease treatment, with emphasis on both cancer management and circadian timing system–resetting strategies for improving chronic disease control and patient outcomes.

Circadian rhythm influences genome-wide transcriptional responses to (131)I in a tissue-specific manner in mice

Background

Circadian variation of gene expression is often neglected when ionizing radiation-induced effects are studied, whether in animal models or in cell culture. This study characterized diurnal variation of genome-wide transcriptional regulation and responses of potential biomarkers and signature genes in normal mouse tissues at 24 h after i.v. administration of ¹³¹I.

Methods

Female BALB/c nude mice were i.v. injected with 90 kBq ¹³¹I at 9:00 a.m., 12:00 p.m., or 3:00 p.m. and killed after 24 h (n = 4/group). Paired control groups were mock-treated (n = 3–4/group). The kidneys, liver, lungs, spleen, and thyroid were excised, snap-frozen, and stored at −80 °C until extraction of total RNA. RNA microarray technology was used for genome-wide expression analysis. Enriched biological processes were categorized after cellular function. Signature genes for ionizing radiation and thyroid hormone-induced responses were taken from the literature. Absorbed dose was estimated using the Medical Internal Radiation Dose (MIRD) formalism.

Results

The thyroid received an absorbed dose of 5.9 Gy and non-thyroid tissues received 0.75–2.2 mGy over 24 h. A distinct peak in the total number of significantly regulated transcripts was observed at 9:00 a.m. in the thyroid, but 3 h later in the kidney cortex, kidney medulla, and liver. Transcriptional regulation in the lungs and spleen was marginal. Associated cellular functions generally varied in quality and response strength between morning, noon, and afternoon. In the thyroid, 25 genes were significantly regulated at all investigated times of day, and 24 thereof showed a distinct pattern of pronounced down-regulation at 9:00 a.m. and comparatively weak up-regulation at later times. Eleven of these genes belonged to the species-specific kallikrein subfamily Klk1b. Responses in signature genes for thyroid hormone-induced responses were more frequent than for ionizing radiation, and trends persisted irrespective of time of day.

Conclusion

Diurnal variation of genome-wide transcriptional responses to 90 kBq ¹³¹I was demonstrated for the thyroid, kidney cortex and medulla, and liver, whereas variation was only marginal in the lungs and spleen. Overall, significant detection of potential biomarkers and signature genes was validated at each time of day, although direction of regulation and fold-change differed between morning, noon, and afternoon. These findings suggest that circadian rhythm should be considered in radiation research and that biological and analytical endpoints should be validated for circadian robustness.

Electronic supplementary material

The online version of this article (doi:10.1186/s13550-015-0150-y) contains supplementary material, which is available to authorized users.

Predicted Role of NAD Utilization in the Control of Circadian Rhythms during DNA Damage Response

The circadian clock is a set of regulatory steps that oscillate with a period of approximately 24 hours influencing many biological processes. These oscillations are robust to external stresses, and in the case of genotoxic stress (i.e. DNA damage), the circadian clock responds through phase shifting with primarily phase advancements. The effect of DNA damage on the circadian clock and the mechanism through which this effect operates remains to be thoroughly investigated. Here we build an in silico model to examine damage-induced circadian phase shifts by investigating a possible mechanism linking circadian rhythms to metabolism. The proposed model involves two DNA damage response proteins, SIRT1 and PARP1, that are each consumers of nicotinamide adenine dinucleotide (NAD), a metabolite involved in oxidation-reduction reactions and in ATP synthesis. This model builds on two key findings: 1) that SIRT1 (a protein deacetylase) is involved in both the positive (i.e. transcriptional activation) and negative (i.e. transcriptional repression) arms of the circadian regulation and 2) that PARP1 is a major consumer of NAD during the DNA damage response. In our simulations, we observe that increased PARP1 activity may be able to trigger SIRT1-induced circadian phase advancements by decreasing SIRT1 activity through competition for NAD supplies. We show how this competitive inhibition may operate through protein acetylation in conjunction with phosphorylation, consistent with reported observations. These findings suggest a possible mechanism through which multiple perturbations, each dominant during different points of the circadian cycle, may result in the phase advancement of the circadian clock seen during DNA damage.

Many physiological processes are regulated by the circadian clock, and we are continuing to learn about the role of the circadian clock in disease. Research in recent years has begun to shed light on the feedback mechanisms that exist between circadian regulation and other processes, including metabolism and the response to DNA damage. A challenge has been to understand the dynamic nature of the protein interactions of these processes, which often involve protein modification as a means of communicating cellular states, such as damaged DNA. Here we have devised a model that simulates an alteration of the circadian clock that is observed during DNA damage response. A novel aspect of this model is the inclusion of SIRT1, a protein that regulates core circadian proteins through modification and helps to repress gene expression. SIRT1 is dependent on a metabolite regulated by the circadian clock and is depleted during DNA damage. In conjunction with a second form of protein modification, our results suggest that multiple forms of protein modification may contribute to the experimentally observed alterations to circadian function.

Circadian rhythms synchronize mitosis in Neurospora crassa

The cell cycle and the circadian clock communicate with each other, resulting in circadian-gated cell division cycles. Alterations in this network may lead to diseases such as cancer. Therefore, it is critical to identify molecular components that connect these two oscillators. However, molecular mechanisms between the clock and the cell cycle remain largely unknown. A model filamentous fungus, Neurospora crassa, is a multinucleate system used to elucidate molecular mechanisms of circadian rhythms, but not used to investigate the molecular coupling between these two oscillators. In this report, we show that a conserved coupling between the circadian clock and the cell cycle exists via serine/threonine protein kinase-29 (STK-29), the Neurospora homolog of mammalian WEE1 kinase. Based on this finding, we established a mathematical model that predicts circadian oscillations of cell cycle components and circadian clock-dependent synchronized nuclear divisions. We experimentally demonstrate that G1 and G2 cyclins, CLN-1 and CLB-1, respectively, oscillate in a circadian manner with bioluminescence reporters. The oscillations of clb-1 and stk-29 gene expression are abolished in a circadian arrhythmic frq(ko) mutant. Additionally, we show the light-induced phase shifts of a core circadian component, frq, as well as the gene expression of the cell cycle components clb-1 and stk-29, which may alter the timing of divisions. We then used a histone hH1-GFP reporter to observe nuclear divisions over time, and show that a large number of nuclear divisions occur in the evening. Our findings demonstrate the circadian clock-dependent molecular dynamics of cell cycle components that result in synchronized nuclear divisions in Neurospora.

Heat shock partially dissociates the overlapping modules of the yeast protein-protein interaction network: a systems level model of adaptation

Network analysis became a powerful tool giving new insights to the understanding of cellular behavior. Heat shock, the archetype of stress responses, is a well-characterized and simple model of cellular dynamics. S. cerevisiae is an appropriate model organism, since both its protein-protein interaction network (interactome) and stress response at the gene expression level have been well characterized. However, the analysis of the reorganization of the yeast interactome during stress has not been investigated yet. We calculated the changes of the interaction-weights of the yeast interactome from the changes of mRNA expression levels upon heat shock. The major finding of our study is that heat shock induced a significant decrease in both the overlaps and connections of yeast interactome modules. In agreement with this the weighted diameter of the yeast interactome had a 4.9-fold increase in heat shock. Several key proteins of the heat shock response became centers of heat shock-induced local communities, as well as bridges providing a residual connection of modules after heat shock. The observed changes resemble to a 'stratus-cumulus' type transition of the interactome structure, since the unstressed yeast interactome had a globally connected organization, similar to that of stratus clouds, whereas the heat shocked interactome had a multifocal organization, similar to that of cumulus clouds. Our results showed that heat shock induces a partial disintegration of the global organization of the yeast interactome. This change may be rather general occurring in many types of stresses. Moreover, other complex systems, such as single proteins, social networks and ecosystems may also decrease their inter-modular links, thus develop more compact modules, and display a partial disintegration of their global structure in the initial phase of crisis. Thus, our work may provide a model of a general, system-level adaptation mechanism to environmental changes.

Multiscale complexity in the mammalian circadian clock

The field of systems biology studies how the interactions among individual components (e.g. genes and proteins) yield interesting and complex behavior. The circadian (daily) timekeeping system in mammals is an ideal system to study complexity because of its many biological scales (from genes to animal behavior). A wealth of data at each of these scales has recently been discovered. Within each scale, modeling can advance our understanding of challenging problems that arise in studying mammalian timekeeping. However, future work must focus on bridging the multiple spatial and temporal scales in the modeling of SCN network. Here we review recent advances, and then delve into a few areas that are promising research directions. We also discuss the flavor of modeling needed (simple or detailed) as well as new techniques that are needed to meet the challenges in modeling data across scales.

BioModels.net Web Services, a free and integrated toolkit for computational modelling software

Exchanging and sharing scientific results are essential for researchers in the field of computational modelling. BioModels.net defines agreed-upon standards for model curation. A fundamental one, MIRIAM (Minimum Information Requested in the Annotation of Models), standardises the annotation and curation process of quantitative models in biology. To support this standard, MIRIAM Resources maintains a set of standard data types for annotating models, and provides services for manipulating these annotations. Furthermore, BioModels.net creates controlled vocabularies, such as SBO (Systems Biology Ontology) which strictly indexes, defines and links terms used in Systems Biology. Finally, BioModels Database provides a free, centralised, publicly accessible database for storing, searching and retrieving curated and annotated computational models. Each resource provides a web interface to submit, search, retrieve and display its data. In addition, the BioModels.net team provides a set of Web Services which allows the community to programmatically access the resources. A user is then able to perform remote queries, such as retrieving a model and resolving all its MIRIAM Annotations, as well as getting the details about the associated SBO terms. These web services use established standards. Communications rely on SOAP (Simple Object Access Protocol) messages and the available queries are described in a WSDL (Web Services Description Language) file. Several libraries are provided in order to simplify the development of client software. BioModels.net Web Services make one step further for the researchers to simulate and understand the entirety of a biological system, by allowing them to retrieve biological models in their own tool, combine queries in workflows and efficiently analyse models.