Modulating the dynamics of NFκB and PI3K enhances the ensemble-level TNFR1 signaling mediated apoptotic response.

Cell-to-cell variability during TNFα stimulated Tumor Necrosis Factor Receptor 1 (TNFR1) signaling can lead to single-cell level pro-survival and apoptotic responses. This variability stems from the heterogeneity in signal flow through intracellular signaling entities that regulate the balance between these two phenotypes. Using systematic Boolean dynamic modeling of a TNFR1 signaling network, we demonstrate that the signal flow path variability can be modulated to enable cells favour apoptosis. We developed a computationally efficient approach "Boolean Modeling based Prediction of Steady-state probability of Phenotype Reachability (BM-ProSPR)" to accurately predict the network's ability to settle into different phenotypes. Model analysis juxtaposed with the experimental observations revealed that NFκB and PI3K transient responses guide the XIAP behaviour to coordinate the crucial dynamic cross-talk between the pro-survival and apoptotic arms at the single-cell level. Model predicted the experimental observations that ~31% apoptosis increase can be achieved by arresting Comp1 - IKK* activity which regulates the NFκB and PI3K dynamics. Arresting Comp1 - IKK* activity causes signal flow path re-wiring towards apoptosis without significantly compromising NFκB levels, which govern adequate cell survival. Priming an ensemble of cancerous cells with inhibitors targeting the specific interaction involving Comp1 and IKK* prior to TNFα exposure could enable driving them towards apoptosis.

Modulation of signaling cross-talk between pJNK and pAKT generates optimal apoptotic response.

Tumor necrosis factor alpha (TNFα) is a well-known modulator of apoptosis by maintaining a balance between proliferation and cell-death in normal cells. Cancer cells often evade apoptotic response following TNFα stimulation by altering signaling cross-talks. Thus, varying the extent of signaling cross-talk could enable optimal TNFα mediated apoptotic dynamics. Herein, we use an experimental data-driven mathematical modeling to quantitate the extent of synergistic signaling cross-talk between the intracellular entities phosphorylated JNK (pJNK) and phosphorylated AKT (pAKT) that orchestrate the phenotypic apoptosis level by modulating the activated Caspase3 dynamics. Our study reveals that this modulation is orchestrated by the distinct dynamic nature of the synergism at early and late phases. We show that this synergism in signal flow is governed by branches originating from either TNFα receptor and NFκB, which facilitates signaling through survival pathways. We demonstrate that the experimentally quantified apoptosis levels semi-quantitatively correlates with the model simulated Caspase3 transients. Interestingly, perturbing pJNK and pAKT transient dynamics fine-tunes this accumulated Caspase3 guided apoptotic response. Thus, our study offers useful insights for identifying potential targeted therapies for optimal apoptotic response.

Retroactivity induced operating regime transition in an enzymatic futile cycle.

Activated phosphorylation-dephosphorylation biochemical reaction cycles are a class of enzymatic futile cycles. A futile cycle such as a single MAPK cascade governed by two underlying enzymatic reactions permits Hyperbolic (H), Signal transducing (ST), Threshold-hyperbolic (TH) and Ultrasensitive (U) operating regimes that characterize input-output behaviour. Retroactive signalling caused by load due to sequestration of phosphorylated or unphosphorylated form of the substrate in a single enzymatic cascade without explicit feedback can introduce two-way communication, a feature not possible otherwise. We systematically characterize the operating regimes of a futile cycle subject to retroactivity in either of the substrate forms. We demonstrate that increasing retroactivity strength, which quantifies the downstream load, can trigger five possible regime transitions. Retroactivity strength is a reflection of the fraction of the substrate sequestered by its downstream target. Remarkably, the minimum required retroactivity strength to evidence any sequestration triggered regime transition demands 23% of the substrate bound to its downstream target. This minimum retroactivity strength corresponds to the transition of the dose-response curve from ST to H regime. We show that modulation of the saturation and unsaturation levels of the enzymatic reactions by retroactivity is the fundamental mechanism governing operating regime transition.

Operating regimes in a single enzymatic cascade at ensemble-level.

Single enzymatic cascade, ubiquitously found in cellular signaling networks, is a phosphorylation-dephosphorylation reaction cycle causing a transition between inactive and active states of a protein catalysed by kinase and phosphatase, respectively. Steady-state information processing ability of such a cycle (e.g., MAPK cascade) has been classified into four qualitatively different operating regimes, viz., hyperbolic (H), signal-transducing (ST), threshold-hyperbolic (TH) and ultrasensitive (U). These four regimes represent qualitatively different dose-response curves, that is, relationship between concentrations of input kinase (e.g., pMEK) and response activated protein (e.g., pERK). Regimes were identified using a deterministic model accounting for population-averaged behavior only. Operating regimes can be strongly influenced by the inherently present cell-to-cell variability in an ensemble of cells which is captured in the form of pMEK and pERK distributions using reporter-based single-cell experimentation. In this study, we show that such experimentally acquired snapshot pMEK and pERK distribution data of a single MAPK cascade can be directly used to infer the underlying operating regime even in the absence of a dose-response curve. This deduction is possible primarily due to the presence of a monotonic relationship between experimental observables RIQR, ratio of the inter-quartile range of the pERK and pMEK distribution pairs and RM, ratio of the medians of the distribution pair. We demonstrate this relationship by systematic analysis of a quasi-steady state approximated model superimposed with an input gamma distribution constrained by the stimulus strength specific pMEK distribution measured on Jurkat-T cells stimulated with PMA. As a first, we show that introduction of cell-to-cell variability only in the upstream kinase achieved by superimposition of an appropriate input pMEK distribution on the dose-response curve can predict bimodal response pERK distribution in ST regime. Implementation of the proposed method on the input-response distribution pair obtained in stimulated Jurkat-T cells revealed that while low-dosage PMA stimulation preserves the H regime observed in resting cells, high-dosage causes H to ST regime transition.

Coupling of Cellular Processes and Their Coordinated Oscillations under Continuous Light in Cyanothece sp. ATCC 51142, a Diazotrophic Unicellular Cyanobacterium.

Unicellular diazotrophic cyanobacteria such as Cyanothece sp. ATCC 51142 (henceforth Cyanothece), temporally separate the oxygen sensitive nitrogen fixation from oxygen evolving photosynthesis not only under diurnal cycles (LD) but also in continuous light (LL). However, recent reports demonstrate that the oscillations in LL occur with a shorter cycle time of ~11 h. We find that indeed, majority of the genes oscillate in LL with this cycle time. Genes that are upregulated at a particular time of day under diurnal cycle also get upregulated at an equivalent metabolic phase under LL suggesting tight coupling of various cellular events with each other and with the cell's metabolic status. A number of metabolic processes get upregulated in a coordinated fashion during the respiratory phase under LL including glycogen degradation, glycolysis, oxidative pentose phosphate pathway, and tricarboxylic acid cycle. These precede nitrogen fixation apparently to ensure sufficient energy and anoxic environment needed for the nitrogenase enzyme. Photosynthetic phase sees upregulation of photosystem II, carbonate transport, carbon concentrating mechanism, RuBisCO, glycogen synthesis and light harvesting antenna pigment biosynthesis. In Synechococcus elongates PCC 7942, a non-nitrogen fixing cyanobacteria, expression of a relatively smaller fraction of genes oscillates under LL condition with the major periodicity being 24 h. In contrast, the entire cellular machinery of Cyanothece orchestrates coordinated oscillation in anticipation of the ensuing metabolic phase in both LD and LL. These results may have important implications in understanding the timing of various cellular events and in engineering cyanobacteria for biofuel production.

Influence of mixotrophic growth on rhythmic oscillations in expression of metabolic pathways in diazotrophic cyanobacterium Cyanothece sp. ATCC 51142.

This study investigates the influence of mixotrophy on physiology and metabolism by analysis of global gene expression in unicellular diazotrophic cyanobacterium Cyanothece sp. ATCC 51142 (henceforth Cyanothece 51142). It was found that Cyanothece 51142 continues to oscillate between photosynthesis and respiration in continuous light under mixotrophy with cycle time of ∼ 13 h. Mixotrophy is marked by an extended respiratory phase compared with photoautotrophy. It can be argued that glycerol provides supplementary energy for nitrogen fixation, which is derived primarily from the glycogen reserves during photoautotrophy. The genes of NDH complex, cytochrome c oxidase and ATP synthase are significantly overexpressed in mixotrophy during the day compared to autotrophy with synchronous expression of the bidirectional hydrogenase genes possibly to maintain redox balance. However, nitrogenase complex remains exclusive to nighttime metabolism concomitantly with uptake hydrogenase. This study throws light on interrelations between metabolic pathways with implications in design of hydrogen producer strains.

Rhythmic and sustained oscillations in metabolism and gene expression of Cyanothece sp. ATCC 51142 under constant light.

Cyanobacteria, a group of photosynthetic prokaryotes, oscillate between day and night time metabolisms with concomitant oscillations in gene expression in response to light/dark cycles (LD). The oscillations in gene expression have been shown to sustain in constant light (LL) with a free running period of 24 h in a model cyanobacterium Synechococcus elongatus PCC 7942. However, equivalent oscillations in metabolism are not reported under LL in this non-nitrogen fixing cyanobacterium. Here we focus on Cyanothece sp. ATCC 51142, a unicellular, nitrogen-fixing cyanobacterium known to temporally separate the processes of oxygenic photosynthesis and oxygen-sensitive nitrogen fixation. In a recent report, metabolism of Cyanothece 51142 has been shown to oscillate between photosynthetic and respiratory phases under LL with free running periods that are temperature dependent but significantly shorter than the circadian period. Further, the oscillations shift to circadian pattern at moderate cell densities that are concomitant with slower growth rates. Here we take this understanding forward and demonstrate that the ultradian rhythm under LL sustains at much higher cell densities when grown under turbulent regimes that simulate flashing light effect. Our results suggest that the ultradian rhythm in metabolism may be needed to support higher carbon and nitrogen requirements of rapidly growing cells under LL. With a comprehensive Real time PCR based gene expression analysis we account for key regulatory interactions and demonstrate the interplay between clock genes and the genes of key metabolic pathways. Further, we observe that several genes that peak at dusk in Synechococcus peak at dawn in Cyanothece and vice versa. The circadian rhythm of this organism appears to be more robust with peaking of genes in anticipation of the ensuing photosynthetic and respiratory metabolic phases.

SHARP: genome-scale identification of gene-protein-reaction associations in cyanobacteria.

Genome scale metabolic model provides an overview of an organism's metabolic capability. These genome-specific metabolic reconstructions are based on identification of gene to protein to reaction (GPR) associations and, in turn, on homology with annotated genes from other organisms. Cyanobacteria are photosynthetic prokaryotes which have diverged appreciably from their nonphotosynthetic counterparts. They also show significant evolutionary divergence from plants, which are well studied for their photosynthetic apparatus. We argue that context-specific sequence and domain similarity can add to the repertoire of the GPR associations and significantly expand our view of the metabolic capability of cyanobacteria. We took an approach that combines the results of context-specific sequence-to-sequence similarity search with those of sequence-to-profile searches. We employ PSI-BLAST for the former, and CDD, Pfam, and COG for the latter. An optimization algorithm was devised to arrive at a weighting scheme to combine the different evidences with KEGG-annotated GPRs as training data. We present the algorithm in the form of software "Systematic, Homology-based Automated Re-annotation for Prokaryotes (SHARP)." We predicted 3,781 new GPR associations for the 10 prokaryotes considered of which eight are cyanobacteria species. These new GPR associations fall in several metabolic pathways and were used to annotate 7,718 gaps in the metabolic network. These new annotations led to discovery of several pathways that may be active and thereby providing new directions for metabolic engineering of these species for production of useful products. Metabolic model developed on such a reconstructed network is likely to give better phenotypic predictions.

Rhythm of carbon and nitrogen fixation in unicellular cyanobacteria under turbulent and highly aerobic conditions.

Nitrogen fixing cyanobacteria are being increasingly explored for nitrogenase-dependent hydrogen production. Commercial success however will depend on the ability to grow these cultures at high cell densities. Photo-limitation at high cell densities leads to hindered photoautotrophic growth while turbulent conditions, which simulate flashing light effect, can lead to oxygen toxicity to the nitrogenase enzyme. Cyanothece sp. strain ATCC 51142, a known hydrogen producer, is reported to grow and fix nitrogen under moderately oxic conditions in shake flasks. In this study, we explore the growth and nitrogen fixing potential of this organism under turbulent conditions with volumetric oxygen mass transfer coefficient (KL a) values that are up to 20-times greater than in shake flasks. In a stirred vessel, the organism grows well in turbulent regime possibly due to a simulated flashing light effect with optimal growth at Reynolds number of approximately 35,000. A respiratory burst lasting for about 4 h creates anoxic conditions intracellularly with near saturating levels of dissolved oxygen in the extracellular medium. This is concomitant with complete exhaustion of intracellular glycogen storage and upregulation of nifH and nifX, the genes encoding proteins of the nitrogenase complex. Further, the rhythmic oscillations in exhaust gas CO2 and O2 profiles synchronize faithfully with those in biochemical parameters and gene expression thereby serving as an effective online monitoring tool. These results will have important implications in potential commercial success of nitrogenase-dependent hydrogen production by cyanobacteria.

Noise propagation in two-step series MAPK cascade.

Series MAPK enzymatic cascades, ubiquitously found in signaling networks, act as signal amplifiers and play a key role in processing information during signal transduction in cells. In activated cascades, cell-to-cell variability or noise is bound to occur and thereby strongly affects the cellular response. Commonly used linearization method (LM) applied to Langevin type stochastic model of the MAPK cascade fails to accurately predict intrinsic noise propagation in the cascade. We prove this by using extensive stochastic simulations for various ranges of biochemical parameters. This failure is due to the fact that the LM ignores the nonlinear effects on the noise. However, LM provides a good estimate of the extrinsic noise propagation. We show that the correct estimate of intrinsic noise propagation in signaling networks that contain at least one enzymatic step can be obtained only through stochastic simulations. Noise propagation in the cascade depends on the underlying biochemical parameters which are often unavailable. Based on a combination of global sensitivity analysis (GSA) and stochastic simulations, we developed a systematic methodology to characterize noise propagation in the cascade. GSA predicts that noise propagation in MAPK cascade is sensitive to the total number of upstream enzyme molecules and the total number of molecules of the two substrates involved in the cascade. We argue that the general systematic approach proposed and demonstrated on MAPK cascade must accompany noise propagation studies in biological networks.

Stability analysis of the GAL regulatory network in Saccharomyces cerevisiae and Kluyveromyces lactis.

BACKGROUND: In the yeast Saccharomyces cerevisiae, interactions between galactose, Gal3p, Gal80p, and Gal4p determine the transcriptional status of the genes required for the galactose utilization. Increase in the cellular galactose concentration causes the galactose molecules to bind onto Gal3p which, via Gal80p, activates Gal4p, which induces the GAL3 and GAL80 gene transcription. Recently, a linear time-invariant multi-input multi-output (MIMO) model of this GAL regulatory network has been proposed; the inputs being galactose and Gal4p, and the outputs being the active Gal4p and galactose utilization. Unfortunately, this model assumes the cell culture to be homogeneous, although it is not so in practice. We overcome this drawback by including more biochemical reactions, and derive a quadratic ordinary differential equation (ODE) based model. RESULTS: We show that the model, referred to above, does not exhibit bistability. We establish sufficiency conditions for the domain of attraction of an equilibrium point of our ODE model for the special case of full-state feedback controller. We observe that the GAL regulatory system of Kluyveromyces lactis exhibits an aberration of monotone nonlinearity and apply the Rantzer multipliers to establish a class of stabilizing controllers for this system. CONCLUSION: Feedback in a GAL regulatory system can be used to enhance the cellular memory. We show that the system can be modeled as a quadratic nonlinear system for which the effect of feedback on the domain of attraction of the equilibrium point can be characterized using linear matrix inequality (LMI) conditions that are easily implementable in software. The benefit of this result is that a mathematically sound approach to the synthesis of full-state and partial-state feedback controllers to regulate the cellular memory is now possible, irrespective of the number of state-variables or parameters of interest.

BioPP: a tool for web-publication of biological networks.

BACKGROUND: Cellular processes depend on the function of intracellular molecular networks. The curation of the literature relevant to specific biological pathways is important for many theoretical and experimental research teams and communities. No current tool supports web publication or hosting of user-developed large scale annotated pathway diagrams. Sharing via web publication is needed to allow real-time access to the current literature pathway knowledge-base, both privately within a research team or publicly among the outside research community. Web publication also facilitates team and/or community input into the curation process while allowing centralized control of the curation and validation process. We have developed new tool to address these needs. Biological Pathway Publisher (BioPP) is a software suite for converting CellDesigner Systems Biology Markup Language (CD-SBML) formatted pathways into a web viewable format. The BioPP suite is available for private use and for depositing knowledge-bases into a newly created public repository. RESULTS: BioPP suite is a web-based application that allows pathway knowledge-bases stored in CD-SBML to be web published with an easily navigated user interface. The BioPP suite consists of four interrelated elements: a pathway publisher, an upload web-interface, a pathway repository for user-deposited knowledge-bases and a pathway navigator. Users have the option to convert their CD-SBML files to HTML for restricted use or to allow their knowledge-base to be web-accessible to the scientific community. All entities in all knowledge-bases in the repository are linked to public database entries as well as to a newly created public wiki which provides a discussion forum. CONCLUSION: BioPP tools and the public repository facilitate sharing of pathway knowledge-bases and interactive curation for research teams and scientific communities. BioPP suite is accessible at http://tsb.mssm.edu/pathwayPublisher/broadcast/