A DNA methylation clock associated with age-related illnesses and mortality is accelerated in men with combat PTSD.

DNA methylation patterns at specific cytosine-phosphate-guanine (CpG) sites predictably change with age and can be used to derive "epigenetic age", an indicator of biological age, as opposed to merely chronological age. A relatively new estimator, called "DNAm GrimAge", is notable for its superior predictive ability in older populations regarding numerous age-related metrics like time-to-death, time-to-coronary heart disease, and time-to-cancer. PTSD is associated with premature mortality and frequently has comorbid physical illnesses suggestive of accelerated biological aging. This is the first study to assess DNAm GrimAge in PTSD patients. We investigated the acceleration of GrimAge relative to chronological age, denoted "AgeAccelGrim" in combat trauma-exposed male veterans with and without PTSD using cross-sectional and longitudinal data from two independent well-characterized veteran cohorts. In both cohorts, AgeAccelGrim was significantly higher in the PTSD group compared to the control group (N = 162, 1.26 vs -0.57, p = 0.001 and N = 53, 0.93 vs -1.60 Years, p = 0.008), suggesting accelerated biological aging in both cohorts with PTSD. In 3-year follow-up study of individuals initially diagnosed with PTSD (N = 26), changes in PTSD symptom severity were correlated with AgeAccelGrim changes (r = 0.39, p = 0.049). In addition, the loss of CD28 cell surface markers on CD8 + T cells, an indicator of T-cell senescence/exhaustion that is associated with biological aging, was positively correlated with AgeAccelGrim, suggesting an immunological contribution to the accelerated biological aging. Overall, our findings delineate cellular correlates of biological aging in combat-related PTSD, which may help explain the increased medical morbidity and mortality seen in this disease.

Dual-Color Single-Cell Imaging of the Suprachiasmatic Nucleus Reveals a Circadian Role in Network Synchrony.

The suprachiasmatic nucleus (SCN) acts as a master pacemaker driving circadian behavior and physiology. Although the SCN is small, it is composed of many cell types, making it difficult to study the roles of particular cells. Here we develop bioluminescent circadian reporter mice that are Cre dependent, allowing the circadian properties of genetically defined populations of cells to be studied in real time. Using a Color-Switch PER2::LUCIFERASE reporter that switches from red PER2::LUCIFERASE to green PER2::LUCIFERASE upon Cre recombination, we assess circadian rhythms in two of the major classes of peptidergic neurons in the SCN: AVP (arginine vasopressin) and VIP (vasoactive intestinal polypeptide). Surprisingly, we find that circadian function in AVP neurons, not VIP neurons, is essential for autonomous network synchrony of the SCN and stability of circadian rhythmicity.

Role of enhanced glucocorticoid receptor sensitivity in inflammation in PTSD: insights from computational model for circadian-neuroendocrine-immune interactions.

Although glucocorticoid resistance contributes to increased inflammation, individuals with posttraumatic stress disorder (PTSD) exhibit increased glucocorticoid receptor (GR) sensitivity along with increased inflammation. It is not clear how inflammation coexists with a hyperresponsive hypothalamic-pituitary-adrenal (HPA) axis. To understand this better, we developed and analyzed an integrated mathematical model for the HPA axis and the immune system. We performed mathematical simulations for a dexamethasone (DEX) suppression test and IC50-dexamethasone for cytokine suppression by varying model parameters. The model analysis suggests that increasing the steepness of the dose-response curve for GR activity may reduce anti-inflammatory effects of GRs at the ambient glucocorticoid levels, thereby increasing proinflammatory response. The adaptive response of proinflammatory cytokine-mediated stimulatory effects on the HPA axis is reduced due to dominance of the GR-mediated negative feedback on the HPA axis. To verify these hypotheses, we analyzed the clinical data on neuroendocrine variables and cytokines obtained from war-zone veterans with and without PTSD. We observed significant group differences for cortisol and ACTH suppression tests, proinflammatory cytokines TNFα and IL6, high-sensitivity C-reactive protein, promoter methylation of GR gene, and IC50-DEX for lysozyme suppression. Causal inference modeling revealed significant associations between cortisol suppression and post-DEX cortisol decline, promoter methylation of human GR gene exon 1F (NR3C1-1F), IC50-DEX, and proinflammatory cytokines. We noted significant mediation effects of NR3C1-1F promoter methylation on inflammatory cytokines through changes in GR sensitivity. Our findings suggest that increased GR sensitivity may contribute to increased inflammation; therefore, interventions to restore GR sensitivity may normalize inflammation in PTSD.

A Randomized, Placebo-Controlled Double-Blind Trial of a Closed-Loop Glucagon System for Postbariatric Hypoglycemia.

BACKGROUND: Postbariatric hypoglycemia (PBH) can threaten safety and reduce quality of life. Current therapies are incompletely effective. METHODS: Patients with PBH were enrolled in a double-blind, placebo-controlled, crossover trial to evaluate a closed-loop glucose-responsive automated glucagon delivery system designed to reduce severe hypoglycemia. A hypoglycemia detection and mitigation algorithm was embedded in the artificial pancreas system connected to a continuous glucose monitor (CGM, Dexcom) driving a patch infusion pump (Insulet) filled with liquid investigational glucagon (Xeris) or placebo (vehicle). Sensor/plasma glucose responses to mixed meal were assessed during 2 study visits. The system delivered up to 2 doses of study drug (300/150 µg glucagon or equal-volume vehicle) if triggered by the algorithm. Rescue dextrose was given for plasma glucose <55 mg/dL or neuroglycopenia. RESULTS: Twelve participants (11 females/1 male, age 52 ± 2, 8 ± 1 years postsurgery, mean ± SEM) completed all visits. Predictive hypoglycemia alerts prompted automated drug delivery postmeal, when sensor glucose was 114 ± 7 vs 121 ± 5 mg/dL (P = .39). Seven participants required rescue glucose after vehicle but not glucagon (P = .008). Five participants had severe hypoglycemia (<55 mg/dL) after vehicle but not glucagon (P = .03). Nadir plasma glucose was higher with glucagon vs vehicle (67 ± 3 vs 59 ± 2 mg/dL, P = .004). Plasma glucagon rose after glucagon delivery (1231 ± 187 vs 16 ± 1 pg/mL at 30 minutes, P = .001). No rebound hyperglycemia occurred. Transient infusion site discomfort was reported with both glucagon (n = 11/12) and vehicle (n = 10/12). No other adverse events were observed. CONCLUSION: A CGM-guided closed-loop rescue system can detect imminent hypoglycemia and deliver glucagon, reducing severe hypoglycemia in PBH. CLINICAL TRIALS REGISTRATION: NCT03255629.

Metabolomic analysis of male combat veterans with post traumatic stress disorder.

Posttraumatic stress disorder (PTSD) is associated with impaired major domains of psychology and behavior. Individuals with PTSD also have increased co-morbidity with several serious medical conditions, including autoimmune diseases, cardiovascular disease, and diabetes, raising the possibility that systemic pathology associated with PTSD might be identified by metabolomic analysis of blood. We sought to identify metabolites that are altered in male combat veterans with PTSD. In this case-control study, we compared metabolomic profiles from age-matched male combat trauma-exposed veterans from the Iraq and Afghanistan conflicts with PTSD (n = 52) and without PTSD (n = 51) ('Discovery group'). An additional group of 31 PTSD-positive and 31 PTSD-negative male combat-exposed veterans was used for validation of these findings ('Test group'). Plasma metabolite profiles were measured in all subjects using ultrahigh performance liquid chromatography/tandem mass spectrometry and gas chromatography/mass spectrometry. We identified key differences between PTSD subjects and controls in pathways related to glycolysis and fatty acid uptake and metabolism in the initial 'Discovery group', consistent with mitochondrial alterations or dysfunction, which were also confirmed in the 'Test group'. Other pathways related to urea cycle and amino acid metabolism were different between PTSD subjects and controls in the 'Discovery' but not in the smaller 'Test' group. These metabolic differences were not explained by comorbid major depression, body mass index, blood glucose, hemoglobin A1c, smoking, or use of analgesics, antidepressants, statins, or anti-inflammatories. These data show replicable, wide-ranging changes in the metabolic profile of combat-exposed males with PTSD, with a suggestion of mitochondrial alterations or dysfunction, that may contribute to the behavioral and somatic phenotypes associated with this disease.

Entrainment of Circadian Rhythms Depends on Firing Rates and Neuropeptide Release of VIP SCN Neurons.

The mammalian suprachiasmatic nucleus (SCN) functions as a master circadian pacemaker, integrating environmental input to align physiological and behavioral rhythms to local time cues. Approximately 10% of SCN neurons express vasoactive intestinal polypeptide (VIP); however, it is unknown how firing activity of VIP neurons releases VIP to entrain circadian rhythms. To identify physiologically relevant firing patterns, we optically tagged VIP neurons and characterized spontaneous firing over 3 days. VIP neurons had circadian rhythms in firing rate and exhibited two classes of instantaneous firing activity. We next tested whether physiologically relevant firing affected circadian rhythms through VIP release. We found that VIP neuron stimulation with high, but not low, frequencies shifted gene expression rhythms in vitro through VIP signaling. In vivo, high-frequency VIP neuron activation rapidly entrained circadian locomotor rhythms. Thus, increases in VIP neuronal firing frequency release VIP and entrain molecular and behavioral circadian rhythms. VIDEO ABSTRACT.

Design and Clinical Evaluation of a Novel Low-Glucose Prediction Algorithm with Mini-Dose Stable Glucagon Delivery in Post-Bariatric Hypoglycemia.

BACKGROUND: Postbariatric hypoglycemia (PBH) is a complication of bariatric surgery with limited therapeutic options. We developed an event-based system to predict and detect hypoglycemia based on continuous glucose monitor (CGM) data and recommend delivery of minidose liquid glucagon. METHODS: We performed an iterative development clinical study employing a novel glucagon delivery system: a Dexcom CGM connected to a Windows tablet running a hypoglycemia prediction algorithm and an Omnipod pump filled with an investigational stable liquid glucagon formulation. Meal tolerance testing was performed in seven participants with PBH and history of neuroglycopenia. Glucagon was administered when hypoglycemia was predicted. Primary outcome measures included the safety and feasibility of this system to predict and prevent severe hypoglycemia. Secondary outcomes included hypoglycemia prediction by the prediction algorithm, minimization of time below hypoglycemia threshold using glucagon, and prevention of rebound hyperglycemia. RESULTS: The hypoglycemia prediction algorithm alerted for impending hypoglycemia in the postmeal state, prompting delivery of glucagon (150 µg). After observations of initial incomplete efficacy to prevent hypoglycemia in the first two participants, system modifications were implemented: addition of PBH-specific detection algorithm, increased glucagon dose (300 µg), and a second glucagon dose if needed. These modifications, together with rescue carbohydrates provided to some participants, contributed to progressive improvements in glucose time above the hypoglycemia threshold (75 mg/dL). CONCLUSIONS: Preliminary results indicate that our event-based automatic monitoring algorithm successfully predicted likely hypoglycemia. Minidose glucagon therapy was well tolerated, without prolonged or severe hypoglycemia, and without rebound hyperglycemia.

Ontogeny of Circadian Rhythms and Synchrony in the Suprachiasmatic Nucleus.

In mammals, the suprachiasmatic nucleus (SCN) of the hypothalamus coordinates daily rhythms including sleep-wake, hormone release, and gene expression. The cells of the SCN must synchronize to each other to drive these circadian rhythms in the rest of the body. The ontogeny of circadian cycling and intercellular coupling in the SCN remains poorly understood. Recent in vitro studies have recorded circadian rhythms from the whole embryonic SCN. Here, we tracked the onset and precision of rhythms in PERIOD2 (PER2), a clock protein, within the SCN isolated from embryonic and postnatal mice of undetermined sex. We found that a few SCN cells developed circadian periodicity in PER2 by 14.5 d after mating (E14.5) with no evidence for daily cycling on E13.5. On E15.5, the fraction of competent oscillators increased dramatically corresponding with stabilization of their circadian periods. The cells of the SCN harvested at E15.5 expressed sustained, synchronous daily rhythms. By postnatal day 2 (P2), SCN oscillators displayed the daily, dorsal-ventral phase wave in clock gene expression typical of the adult SCN. Strikingly, vasoactive intestinal polypeptide (VIP), a neuropeptide critical for synchrony in the adult SCN, and its receptor, VPAC2R, reached detectable levels after birth and after the onset of circadian synchrony. Antagonists of GABA or VIP signaling or action potentials did not disrupt circadian synchrony in the E15.5 SCN. We conclude that endogenous daily rhythms in the fetal SCN begin with few noisy oscillators on E14.5, followed by widespread oscillations that rapidly synchronize on E15.5 by an unknown mechanism.SIGNIFICANCE STATEMENT We recorded the onset of PER2 circadian oscillations during embryonic development in the mouse SCN. When isolated at E13.5, the anlagen of the SCN expresses high, arrhythmic PER2. In contrast, a few cells show noisy circadian rhythms in the isolated E14.5 SCN and most show reliable, self-sustained, synchronized rhythms in the E15.5 SCN. Strikingly, this synchrony at E15.5 appears before expression of VIP or its receptor and persists in the presence of blockers of VIP, GABA or neuronal firing. Finally, the dorsal-ventral phase wave of PER2 typical of the adult SCN appears ∼P2, indicating that multiple signals may mediate circadian synchrony during the ontogeny of the SCN.

Quantity and accessibility for specific targeting of receptors in tumours.

Synaphic (ligand-directed) targeting of drugs is an important potential new approach to drug delivery, particularly in oncology. Considerable success with this approach has been achieved in the treatment of blood-borne cancers, but the advances with solid tumours have been modest. Here, we have studied the number and availability for ligand binding of the receptors for two targeting ligands. The results show that both paucity of total receptors and their poor availability are major bottlenecks in drug targeting. A tumour-penetrating peptide greatly increases the availability of receptors by promoting transport of the drug to the extravascular tumour tissue, but the number of available receptors still remains low, severely limiting the utility of the approach. Our results emphasize the importance of using drugs with high specific activity to avoid exceeding receptor capacity because any excess drug conjugate would lose the targeting advantage. The mathematical models we describe make it possible to focus on those aspects of the targeting mechanism that are most likely to have a substantial effect on the overall efficacy of the targeting.

A neuropeptide speeds circadian entrainment by reducing intercellular synchrony.

Shift work or transmeridian travel can desynchronize the body's circadian rhythms from local light-dark cycles. The mammalian suprachiasmatic nucleus (SCN) generates and entrains daily rhythms in physiology and behavior. Paradoxically, we found that vasoactive intestinal polypeptide (VIP), a neuropeptide implicated in synchrony among SCN cells, can also desynchronize them. The degree and duration of desynchronization among SCN neurons depended on both the phase and the dose of VIP. A model of the SCN consisting of coupled stochastic cells predicted both the phase- and the dose-dependent response to VIP and that the transient phase desynchronization, or "phase tumbling", could arise from intrinsic, stochastic noise in small populations of key molecules (notably, Period mRNA near its daily minimum). The model also predicted that phase tumbling following brief VIP treatment would accelerate entrainment to shifted environmental cycles. We tested this using a prepulse of VIP during the day before a shift in either a light cycle in vivo or a temperature cycle in vitro. Although VIP during the day does not shift circadian rhythms, the VIP pretreatment approximately halved the time required for mice to reentrain to an 8-h shifted light schedule and for SCN cultures to reentrain to a 10-h shifted temperature cycle. We conclude that VIP below 100 nM synchronizes SCN cells and above 100 nM reduces synchrony in the SCN. We show that exploiting these mechanisms that transiently reduce cellular synchrony before a large shift in the schedule of daily environmental cues has the potential to reduce jet lag.

Modeling the intra- and extracellular cytokine signaling pathway under heat stroke in the liver.

Heat stroke (HS) is a life-threatening illness induced by prolonged exposure to a hot environment that causes central nervous system abnormalities and severe hyperthermia. Current data suggest that the pathophysiological responses to heat stroke may not only be due to the immediate effects of heat exposure per se but also the result of a systemic inflammatory response syndrome (SIRS). The observation that pro- (e.g., IL-1) and anti-inflammatory (e.g., IL-10) cytokines are elevated concomitantly during recovery suggests a complex network of interactions involved in the manifestation of heat-induced SIRS. In this study, we measured a set of circulating cytokine/soluble cytokine receptor proteins and liver cytokine and receptor mRNA accumulation in wild-type and tumor necrosis factor (TNF) receptor knockout mice to assess the effect of neutralization of TNF signaling on the SIRS following HS. Using a systems approach, we developed a computational model describing dynamic changes (intra- and extracellular events) in the cytokine signaling pathways in response to HS that was fitted to novel genomic (liver mRNA accumulation) and proteomic (circulating cytokines and receptors) data using global optimization. The model allows integration of relevant biological knowledge and formulation of new hypotheses regarding the molecular mechanisms behind the complex etiology of HS that may serve as future therapeutic targets. Moreover, using our unique modeling framework, we explored cytokine signaling pathways with three in silico experiments (e.g. by simulating different heat insult scenarios and responses in cytokine knockout strains in silico).

Collagen cross-linking using rose bengal and green light to increase corneal stiffness.

PURPOSE: Photochemical cross-linking of corneal collagen is an evolving treatment for keratoconus and other ectatic disorders. We evaluated collagen cross-linking by rose bengal plus green light (RGX) in rabbit eyes and investigated factors important for clinical application. METHODS: Rose bengal (RB, 0.1%) was applied to deepithelialized corneas of enucleated rabbit eyes for 2 minutes. The diffusion distance of RB into the stroma was measured by fluorescence microscopy on frozen sections. RB-stained corneas were exposed to green (532-nm) light for 3.3 to 9.9 minutes (50-150 J/cm(2)). Changes in the absorption spectrum during the irradiation were recorded. Corneal stiffness was measured by uniaxial tensiometry. The spatial distribution of the stromal elastic modulus was assessed by Brillouin microscopy. Viable keratocytes were counted on H&E-stained sections 24 hours posttreatment. RESULTS: RB penetrated approximately 100 µm into the corneal stroma and absorbed >90% of the incident green light. RGX (150 J/cm(2)) increased stromal stiffness by 3.8-fold. The elastic modulus increased in the anterior approximately 120 µm of stroma. RB was partially photobleached during the 2-minute irradiation, but reapplication of RB blocked light transmission by >70%. Spectral measurements suggested that RGX initiated cross-linking by an oxygen-dependent mechanism. RGX did not decrease keratocyte viability. CONCLUSIONS: RGX significantly increases cornea stiffness in a rapid treatment (≅12 minutes total time), does not cause toxicity to keratocytes and may be used to stiffen corneas thinner than 400 µm. Thus, RGX may provide an attractive approach to inhibit progression of keratoconus and other ectatic disorders.

Identification of small molecule activators of cryptochrome.

Impairment of the circadian clock has been associated with numerous disorders, including metabolic disease. Although small molecules that modulate clock function might offer therapeutic approaches to such diseases, only a few compounds have been identified that selectively target core clock proteins. From an unbiased cell-based circadian phenotypic screen, we identified KL001, a small molecule that specifically interacts with cryptochrome (CRY). KL001 prevented ubiquitin-dependent degradation of CRY, resulting in lengthening of the circadian period. In combination with mathematical modeling, our studies using KL001 revealed that CRY1 and CRY2 share a similar functional role in the period regulation. Furthermore, KL001-mediated CRY stabilization inhibited glucagon-induced gluconeogenesis in primary hepatocytes. KL001 thus provides a tool to study the regulation of CRY-dependent physiology and aid development of clock-based therapeutics of diabetes.

Core module biomarker identification with network exploration for breast cancer metastasis.

BACKGROUND: In a complex disease, the expression of many genes can be significantly altered, leading to the appearance of a differentially expressed "disease module". Some of these genes directly correspond to the disease phenotype, (i.e. "driver" genes), while others represent closely-related first-degree neighbours in gene interaction space. The remaining genes consist of further removed "passenger" genes, which are often not directly related to the original cause of the disease. For prognostic and diagnostic purposes, it is crucial to be able to separate the group of "driver" genes and their first-degree neighbours, (i.e. "core module") from the general "disease module". RESULTS: We have developed COMBINER: COre Module Biomarker Identification with Network ExploRation. COMBINER is a novel pathway-based approach for selecting highly reproducible discriminative biomarkers. We applied COMBINER to three benchmark breast cancer datasets for identifying prognostic biomarkers. COMBINER-derived biomarkers exhibited 10-fold higher reproducibility than other methods, with up to 30-fold greater enrichment for known cancer-related genes, and 4-fold enrichment for known breast cancer susceptible genes. More than 50% and 40% of the resulting biomarkers were cancer and breast cancer specific, respectively. The identified modules were overlaid onto a map of intracellular pathways that comprehensively highlighted the hallmarks of cancer. Furthermore, we constructed a global regulatory network intertwining several functional clusters and uncovered 13 confident "driver" genes of breast cancer metastasis. CONCLUSIONS: COMBINER can efficiently and robustly identify disease core module genes and construct their associated regulatory network. In the same way, it is potentially applicable in the characterization of any disease that can be probed with microarrays.

Confidence from uncertainty--a multi-target drug screening method from robust control theory.

BACKGROUND: Robustness is a recognized feature of biological systems that evolved as a defence to environmental variability. Complex diseases such as diabetes, cancer, bacterial and viral infections, exploit the same mechanisms that allow for robust behaviour in healthy conditions to ensure their own continuance. Single drug therapies, while generally potent regulators of their specific protein/gene targets, often fail to counter the robustness of the disease in question. Multi-drug therapies offer a powerful means to restore disrupted biological networks, by targeting the subsystem of interest while preventing the diseased network from reconciling through available, redundant mechanisms. Modelling techniques are needed to manage the high number of combinatorial possibilities arising in multi-drug therapeutic design, and identify synergistic targets that are robust to system uncertainty. RESULTS: We present the application of a method from robust control theory, Structured Singular Value or µ- analysis, to identify highly effective multi-drug therapies by using robustness in the face of uncertainty as a new means of target discrimination. We illustrate the method by means of a case study of a negative feedback network motif subject to parametric uncertainty. CONCLUSIONS: The paper contributes to the development of effective methods for drug screening in the context of network modelling affected by parametric uncertainty. The results have wide applicability for the analysis of different sources of uncertainty like noise experienced in the data, neglected dynamics, or intrinsic biological variability.

Fathead minnow steroidogenesis: in silico analyses reveals tradeoffs between nominal target efficacy and robustness to cross-talk.

BACKGROUND: Interpreting proteomic and genomic data is a major challenge in predictive ecotoxicology that can be addressed by a systems biology approach. Mathematical modeling provides an organizational platform to consolidate protein dynamics with possible genomic regulation. Here, a model of ovarian steroidogenesis in the fathead minnow, Pimephales promelas, (FHM) is developed to evaluate possible transcriptional regulation of steroid production observed in microarray studies. RESULTS: The model was developed from literature sources, integrating key signaling components (G-protein and PKA activation) with their ensuing effect on steroid production. The model properly predicted trajectory behavior of estradiol and testosterone when fish were exposed to fadrozole, a specific aromatase inhibitor, but failed to predict the steroid hormone behavior occurring one week post-exposure as well as the increase in steroid levels when the stressor was removed. In vivo microarray data implicated three modes of regulation which may account for over-production of steroids during a depuration phase (when the stressor is removed): P450 enzyme up-regulation, inhibin down-regulation, and luteinizing hormone receptor up-regulation. Simulation studies and sensitivity analysis were used to evaluate each case as possible source of compensation to endocrine stress. CONCLUSIONS: Simulation studies of the testosterone and estradiol response to regulation observed in microarray data supported the hypothesis that the FHM steroidogenesis network compensated for endocrine stress by modulating the sensitivity of the ovarian network to global cues coming from the hypothalamus and pituitary. Model predictions of luteinizing hormone receptor regulation were consistent with depuration and in vitro data. These results challenge the traditional approach to network elucidation in systems biology. Generally, the most sensitive interactions in a network are targeted for further elucidation but microarray evidence shows that homeostatic regulation of the steroidogenic network is likely maintained by a mildly sensitive interaction. We hypothesize that effective network elucidation must consider both the sensitivity of the target as well as the target's robustness to biological noise (in this case, to cross-talk) when identifying possible points of regulation.

Circadian phase resetting via single and multiple control targets.

Circadian entrainment is necessary for rhythmic physiological functions to be appropriately timed over the 24-hour day. Disruption of circadian rhythms has been associated with sleep and neuro-behavioral impairments as well as cancer. To date, light is widely accepted to be the most powerful circadian synchronizer, motivating its use as a key control input for phase resetting. Through sensitivity analysis, we identify additional control targets whose individual and simultaneous manipulation (via a model predictive control algorithm) out-perform the open-loop light-based phase recovery dynamics by nearly 3-fold. We further demonstrate the robustness of phase resetting by synchronizing short- and long-period mutant phenotypes to the 24-hour environment; the control algorithm is robust in the presence of model mismatch. These studies prove the efficacy and immediate application of model predictive control in experimental studies and medicine. In particular, maintaining proper circadian regulation may significantly decrease the chance of acquiring chronic illness.

Identifying fragilities in biochemical networks: robust performance analysis of Fas signaling-induced apoptosis.

Proper control of apoptotic signaling is critical to immune response and development in multicellular organisms. Two tools from control engineering are applied to a mathematical model of Fas ligand signaling-induced apoptosis. Structured singular value analysis determines the volume in parameter space within which the system parameters may exist and still maintain efficacious signaling, but is limited to linear behaviors. Sensitivity analysis can be applied to nonlinear systems but is difficult to relate to performance criteria. Thus, structured singular value analysis is used to quantify performance during apoptosis rejection, ensuring that the system remains sensitive but not overly so to apoptotic stimuli. Sensitivity analysis is applied when the system has switched to the death-inducing, apoptotic steady state to determine parameters significant to maintaining the bistability. The analyses reveal that the magnitude of the death signal is fragile to perturbations in degradation parameters (failures in the ubiquitin/proteasome mechanism) while the timing of signal expression can be tuned by manipulating local parameters. Simultaneous parameter uncertainty highlights apoptotic fragility to disturbances in the ubiquitin/proteasome system. Sensitivity analysis reveals that the robust signaling characteristics of the apoptotic network is due to network architecture, and the apoptotic signaling threshold is best manipulated by interactions upstream of the apoptosome.

A molecular model for intercellular synchronization in the mammalian circadian clock.

The mechanisms and consequences of synchrony among heterogeneous oscillators are poorly understood in biological systems. We present a multicellular, molecular model of the mammalian circadian clock that incorporates recent data implicating the neurotransmitter vasoactive intestinal polypeptide (VIP) as the key synchronizing agent. The model postulates that synchrony arises among circadian neurons because they release VIP rhythmically on a daily basis and in response to ambient light. Two basic cell types, intrinsically rhythmic pacemakers and damped oscillators, are assumed to arise from a distribution of Period gene transcription rates. Postsynaptic neurons show time-of-day dependent responses to VIP binding through a signaling cascade that activates Period mRNA transcription. The heterogeneous cell ensemble model self-synchronizes, entrains to ambient light-dark cycles, and desynchronizes in constant bright light or upon removal of VIP signaling. The degree of synchronicity observed depends on cell-specific features (e.g., mean and variability of parameters within the rhythm-generating loop), in addition to the more commonly studied effect of intercellular coupling strength. These simulations closely replicate experimental data and predict that heterogeneous oscillations (e.g., sustained, damped, and arrhythmic) arise from small differences in the molecular parameters between cells, that damped oscillators participate in entrainment and synchrony of the ensemble of cells, and that constant light desynchronizes oscillators by maximizing VIP release.