Search

Over recent years, new light has been shed on aspects of information processing in cells. The quantification of information, as described by Shannon’s information theory, is a basic and powerful tool that can be applied to various fields, such as communication, statistics, and computer science, as well as to information processing within cells. It has also been used to infer the network structure of molecular species. However, the difficulty of obtaining sufficient sample sizes and the computational burden associated with the high-dimensional data often encountered in biology can result in bottlenecks in the application of information theory to systems biology. This article provides an overview of the application of information theory to systems biology, discussing the associated bottlenecks and reviewing recent work.

SUMMARY The thymic function to produce self-protective and self-tolerant T cells is chiefly mediated by cortical thymic epithelial cells (cTECs) and medullary TECs (mTECs). Recent studies including single-cell transcriptomic analyses have highlighted a rich diversity in functional mTEC subpopulations. Because of their limited cellularity, however, the biochemical characterization of TECs, including the proteomic profiling of cTECs and mTECs, has remained unestablished. Utilizing genetically modified mice that carry enlarged but functional thymuses, here we show a combination of proteomic and transcriptomic profiles for cTECs and mTECs, which identified signature molecules that characterize a developmental and functional contrast between cTECs and mTECs. Our results reveal a highly specific impact of the thymoproteasome on proteasome subunit composition in cTECs and provide an integrated trans-omics platform for further exploration of thymus biology., In Brief Ohigashi et al. show that the use of cyclin D1-transgenic mice allows quantitative proteomic analysis of cortical and medullary thymic epithelial cells (TECs). Results provide a trans-omics platform for further exploration of TEC biology and reveal the specific impact of the thymoproteasome on proteasome subunit composition in cortical TECs., Graphical Abstract

Summary An effective combination of multi-omic datasets can enhance our understanding of complex biological phenomena. To build a context-dependent network with multiple omic layers, i.e., a trans-omic network, we perform phosphoproteomics, transcriptomics, proteomics, and metabolomics of murine liver for 4 h after insulin administration and integrate the resulting time series. Structural characteristics and dynamic nature of the network are analyzed to elucidate the impact of insulin. Early and prominent changes in protein phosphorylation and persistent and asynchronous changes in mRNA and protein levels through non-transcriptional mechanisms indicate enhanced crosstalk between phosphorylation-mediated signaling and protein expression regulation. Metabolic response shows different temporal regulation with transient increases at early time points across categories and enhanced response in the amino acid and nucleotide categories at later time points as a result of process convergence. This extensive and dynamic view of the trans-omic network elucidates prominent regulatory mechanisms that drive insulin responses through intricate interlayer coordination.

Summary Systemic metabolic homeostasis is regulated by inter-organ metabolic cycles involving multiple organs. Obesity impairs inter-organ metabolic cycles, resulting in metabolic diseases. The systemic landscape of dysregulated inter-organ metabolic cycles in obesity has yet to be explored. Here, we measured the transcriptome, proteome, and metabolome in the liver and skeletal muscle and the metabolome in blood of fasted wild-type and leptin-deficient obese (ob/ob) mice, identifying components with differential abundance and differential regulation in ob/ob mice. By constructing and evaluating the trans-omic network controlling the differences in metabolic reactions between fasted wild-type and ob/ob mice, we provided potential mechanisms of the obesity-associated dysfunctions of metabolic cycles between liver and skeletal muscle involving glucose-alanine, glucose-lactate, and ketone bodies. Our study revealed obesity-associated systemic pathological mechanisms of dysfunction of inter-organ metabolic cycles., Graphical abstract, Highlights • Multi-omic data in liver and skeletal muscle of WT and ob/ob mice were measured • We developed the trans-omic network of differentially regulated metabolic reactions • Dysregulation of inter-organ metabolic cycles associated with obesity was revealed, Biological Sciences ; Endocrinology ; Systems Biology ; Omics ; Proteomics ; Metabolomics ; Transcriptomic

Impaired glucose tolerance associated with obesity causes postprandial hyperglycemia and can lead to type 2 diabetes. To study the differences in liver metabolism in healthy and obese states, we constructed and analyzed transomics glucose-responsive metabolic networks with layers for metabolites, expression data for metabolic enzyme genes, transcription factors, and insulin signaling proteins from the livers of healthy and obese mice. We integrated multiomics time course data from wild-type and leptin-deficient obese (ob/ob) mice after orally administered glucose. In wild-type mice, metabolic reactions were rapidly regulated within 10 min of oral glucose administration by glucose-responsive metabolites, which functioned as allosteric regulators and substrates of metabolic enzymes, and by Akt-induced changes in the expression of glucose-responsive genes encoding metabolic enzymes. In ob/ob mice, the majority of rapid regulation by glucose-responsive metabolites was absent. Instead, glucose administration produced slow changes in the expression of carbohydrate, lipid, and amino acid metabolic enzyme-encoding genes to alter metabolic reactions on a time scale of hours. Few regulatory events occurred in both healthy and obese mice. Thus, our transomics network analysis revealed that regulation of glucose-responsive liver metabolism is mediated through different mechanisms in healthy and obese states. Rapid changes in allosteric regulators and substrates and in gene expression dominate the healthy state, whereas slow changes in gene expression dominate the obese state.

Cellular signaling regulates various cellular functions via protein phosphorylation. Phosphoproteomic data potentially include information for a global regulatory network from signaling to cellular functions, but a procedure to reconstruct this network using such data has yet to be established. In this paper, we provide a procedure to reconstruct a global regulatory network from signaling to cellular functions from phosphoproteomic data by integrating prior knowledge of cellular functions and inference of the kinase-substrate relationships (KSRs). We used phosphoproteomic data from insulin-stimulated Fao hepatoma cells and identified protein phosphorylation regulated by insulin specifically over-represented in cellular functions in the KEGG database. We inferred kinases for protein phosphorylation by KSRs, and connected the kinases in the insulin signaling layer to the phosphorylated proteins in the cellular functions, revealing that the insulin signal is selectively transmitted via the Pi3k-Akt and Erk signaling pathways to cellular adhesions and RNA maturation, respectively. Thus, we provide a method to reconstruct global regulatory network from signaling to cellular functions based on phosphoproteomic data.

Summary The concentrations of insulin selectively regulate multiple cellular functions. To understand how insulin concentrations are interpreted by cells, we constructed a trans-omic network of insulin action in FAO hepatoma cells using transcriptomic data, western blotting analysis of signaling proteins, and metabolomic data. By integrating sensitivity into the trans-omic network, we identified the selective trans-omic networks stimulated by high and low doses of insulin, denoted as induced and basal insulin signals, respectively. The induced insulin signal was selectively transmitted through the pathway involving Erk to an increase in the expression of immediate-early and upregulated genes, whereas the basal insulin signal was selectively transmitted through a pathway involving Akt and an increase of Foxo phosphorylation and a reduction of downregulated gene expression. We validated the selective trans-omic network in vivo by analysis of the insulin-clamped rat liver. This integrated analysis enabled molecular insight into how liver cells interpret physiological insulin signals to regulate cellular functions., Graphical Abstract, Highlights • We constructed a trans-omic network of insulin action using multi-omic data • The trans-omic network integrates phosphorylation, transcription, and metabolism • We divided induced- and basal-insulin-dependent trans-omic network by sensitivity • Induced and basal insulin triggered selective trans-omic network in the liver, Systems Biology; Omics; Metabolomics; Transcriptomics

Impaired glucose tolerance associated with obesity causes postprandial hyperglycemia and can lead to type 2 diabetes. To study the differences in liver metabolism in the healthy and obese states, we constructed and analyzed trans-omic glucose-responsive metabolic networks with layers for metabolites, expression data for metabolic enzyme genes, transcription factors, and insulin signaling proteins from the livers of healthy and obese mice. We integrated multi-omic time-course data from wild-type (WT) and leptin-deficient obese (ob/ob) mice after orally administered glucose. In WT mice, metabolic reactions were rapidly regulated (within 10 minutes of oral glucose administration) primarily by glucose-responsive metabolites, especially ATP and NADP+, which functioned as allosteric regulators and substrates of metabolic enzymes, and by Akt-dependent glucose-responsive genes encoding metabolic enzymes. Inob/obmice, most rapid regulation by glucose-responsive metabolites was absent; instead, glucose administration produced slow changes in the expression of metabolic enzyme-encoding genes to alter metabolic reactions in a time scale of hours. Few common regulatory events occurred in both the healthy and obese mice. Thus, our trans-omic network analysis revealed regulation of liver metabolism in response to glucose is mediated through different mechanisms in the healthy and obese states: Rapid changes in allosteric regulators and substrates and in gene expression dominate the healthy state, and slow transcriptional regulation dominates the obese state.One Sentence SummaryRapid changes in regulatory metabolites and gene expression dominate the healthy state, and slow transcriptional regulation dominates the obese state.

Cells transmit information about extracellular stimulation through signaling pathways to control cellular function. A signaling pathway can be regarded as a communication channel. In the analysis of channels of cell populations, intercellular variation is considered noise. However, intercellular variation enables individual cells to encode different information. Therefore, at the single-cell level, each cell can be regarded as an independent channel. Thus, we propose that responses of cells of the same type in tissues, such as the fibers in a skeletal muscle, should be regarded as a multiple-cell channel composed of single-cell channels, in which intercellular variation contains information. Here, we applied electrical pulses to individual myotubes from cultured C2C12 cells or dissociated skeletal muscle fibers and measured Ca2+ responses or contraction, respectively, to estimate information capacity in a biological system. For each muscle cell system, we found that a single-cell channel transmitted more information than did a cell-population channel, indicating that the cellular response is consistent with each cell (low intracellular variation) but different among individual cells (high intercellular variation). As cell number and thus the number of single-cell channels increased, a multiple-cell channel transmitted more information by incorporating the differences among individual cells. Thus, a tissue with small intracellular and large intercellular variations has the capacity to distinguish differences in stimulation intensity to precisely control physiological function.One Sentence SummarySmall intracellular and large intercellular variations increase information transmission for precise control of tissue function.

Signal transduction processes the information of various cellular functions, including cell proliferation, differentiation, and death. The information for controlling cell fate is transmitted by concentrations of cellular signaling molecules. However, how much information is transmitted in signaling pathways has thus far not been investigated. Shannon's information theory paves the way to quantitatively analyze information transmission in signaling pathways. The theory has recently been applied to signal transduction, and mutual information of signal transduction has been determined to be a measure of information transmission. We review this work and provide an overview of how signal transduction transmits informational input and exerts biological output.

Cell-to-cell variability in signal transduction in biological systems is often considered noise. However, intercellular variation (i.e., cell-to-cell variability) has the potential to enable individual cells to encode different information. Here, we show that intercellular variation increases information transmission of skeletal muscle. We analyze the responses of multiple cultured myotubes or isolated skeletal muscle fibers as a multiple-cell channel composed of single-cell channels. We find that the multiple-cell channel, which incorporates intercellular variation as information, not noise, transmitted more information in the presence of intercellular variation than in the absence according to the "response diversity effect," increasing in the gradualness of dose response by summing the cell-to-cell variable dose responses. We quantify the information transmission of human facial muscle contraction during intraoperative neurophysiological monitoring and find that information transmission of muscle contraction is comparable to that of a multiple-cell channel. Thus, our data indicate that intercellular variation can increase the information capacity of tissues.

Excessive increase in blood glucose level after eating increases the risk of macroangiopathy, and a method for not increasing the postprandial blood glucose level is desired. However, a logical design method of the dietary ingestion pattern controlling the postprandial blood glucose level has not yet been established. We constructed a mathematical model of blood glucose control by oral glucose ingestion in three healthy human subjects, and predicted that intermittent ingestion 30 min apart was the optimal glucose ingestion patterns that minimized the peak value of blood glucose level. We confirmed with subjects that this intermittent pattern consistently decreased the peak value of blood glucose level. We also predicted insulin minimization pattern, and found that the intermittent ingestion 30 min apart was optimal, which is similar to that of glucose minimization pattern. Taken together, these results suggest that the glucose minimization is achieved by suppressing the peak value of insulin concentration, rather than by enhancing insulin concentration. This approach could be applied to design optimal dietary ingestion patterns., The key points of this study are three-fold: The first point is the physiological impact. Intuitively, an ingestion pattern that minimizes the maximum blood glucose level is expected to be a slow, continuous ingestion. We found that the optimal minimization pattern was an intermittent pattern with 30 min intervals. The second point is the methodology. We constructed the mathematical model as a forward problem, and in turn, predicted input pattern by control the output pattern of interest as an inverse problem. The third point is the experiment design used to generate the mathematical model; dense time course data for six ingestion patterns combining of 3 doses of glucose level and 2 durations of ingestion.

Insulin plays a central role in glucose homeostasis, and impairment of insulin action causes glucose intolerance and leads to type 2 diabetes mellitus (T2DM). A decrease in the transient peak and sustained increase of circulating insulin following an infusion of glucose accompany T2DM pathogenesis. However, the mechanism underlying this abnormal temporal pattern of circulating insulin concentration remains unknown. Here we show that changes in opposite direction of hepatic and peripheral insulin clearance characterize this abnormal temporal pattern of circulating insulin concentration observed in T2DM. We developed a mathematical model using a hyperglycemic and hyperinsulinemic-euglycemic clamp in 111 subjects, including healthy normoglycemic and diabetic subjects. The hepatic and peripheral insulin clearance significantly increase and decrease, respectively, from healthy to borderline type and T2DM. The increased hepatic insulin clearance reduces the amplitude of circulating insulin concentration, whereas the decreased peripheral insulin clearance changes the temporal patterns of circulating insulin concentration from transient to sustained. These results provide further insight into the pathogenesis of T2DM, and thus may contribute to develop better treatment of this condition., Glucose homeostasis: roles of insulin clearance Type 2 diabetes mellitus (T2DM) is one of the fastest growing public health problems, characterized by chronic hyperglycemia with the failure of glucose homeostasis. Evaluating alteration in biological functions regulating circulating glucose concentration is complicated due to the mutual relation between circulating glucose and insulin. A team led by Wataru Ogawa at Kobe University designed clinical experiments for breaking such feedback relations, and a team led by Shinya Kuroda at University of Tokyo developed mathematical models for specifically quantifying the functions from the clinical data. The estimated model parameters revealed the significant increase in hepatic and decrease in peripheral insulin clearance, which occur before and after insulin delivery into systemic circulation, respectively, from healthy to T2DM subjects. Model analysis suggested these insulin clearances centrally regulate the dynamics of circulating insulin concentration in the glucose-insulin regulatory system.

Summary Cells respond to various extracellular stimuli through a limited number of signaling pathways. One strategy to process such stimuli is to code the information into the temporal patterns of molecules. Although we showed that insulin selectively regulated molecules depending on its temporal patterns using Fao cells, the in vivo mechanism remains unknown. Here, we show how the insulin-AKT pathway processes the information encoded into the temporal patterns of blood insulin. We performed hyperinsulinemic-euglycemic clamp experiments and found that, in the liver, all temporal patterns of insulin are encoded into the insulin receptor, and downstream molecules selectively decode them through AKT. S6K selectively decodes the additional secretion information. G6Pase interprets the basal secretion information through FoxO1, while GSK3β decodes all secretion pattern information. Mathematical modeling revealed the mechanism via differences in network structures and from sensitivity and time constants. Given that almost all hormones exhibit distinct temporal patterns, temporal coding may be a general principle of system homeostasis by hormones.

This paper deals with an input/output signal recovery problem for nonlinear multiple-input single-output autoregressive exogenous (ARX) models with input nonlinearity, which are used in data-driven systems biology. A matrix rank minimization approach is applied, and a new signal recovery algorithm for nonlinear ARX models is provided. The proposed algorithm recovers output signals and nonlinear-transformed input signals on a linear subspace using some assumptions about nonlinear functions and does not require the exact knowledge of nonlinear functions. Numerical examples using experimental data of signal transduction of cellular systems show the efficiency of the proposed algorithm.

Cells decode information of signaling activation at a scale of tens of minutes by downstream gene expression with a scale of hours to days, leading to cell fate decisions such as cell differentiation. However, no system identification method with such different time scales exists. Here we used compressed sensing technology and developed a system identification method using data of different time scales by recovering signals of missing time points. We measured phosphorylation of ERK and CREB, immediate early gene expression products, and mRNAs of decoder genes for neurite elongation in PC12 cell differentiation and performed system identification, revealing the input–output relationships between signaling and gene expression with sensitivity such as graded or switch-like response and with time delay and gain, representing signal transfer efficiency. We predicted and validated the identified system using pharmacological perturbation. Thus, we provide a versatile method for system identification using data with different time scales., Author summary The key points of this study are two-fold: The first point is the decoding mechanism for cell differentiation. We previously demonstrated the encoding mechanism of cell fate decision information by transient and sustained ERK activation in PC12 cells, and also identified the decoding genes essential for cell differentiation in PC12 cells, including Metrnl, Dclk1, and Serpinb1a, denoted as LP (latent process) genes, which are the decoders of neurite length information. Importantly, the expression levels of the LP genes, but not the phosphorylation level of ERK, correlate with neurite length. Thus, the decoding mechanism of signaling activities by LP gene expression is a key issue for understanding the mechanism of cell differentiation. Here we identified a selective NGF- and PACAP-signaling decoding system by LP gene expression for neurite extension by developing a system identification method. The second point is the modeling. Cells decode information of signaling activation at a scale of tens of minutes by downstream gene expression with a scale of hours to days, leading to cell fate decisions such as cell differentiation. However, no system identification method with such different time scales exists. Here we developed a signal recovery technique in the field of compressed sensing originally developed for image analysis to biological sparse data of different time scales of signaling and gene expression.

SUMMARYThe concentration and temporal pattern of insulin selectively regulate multiple cellular functions. To understand how insulin dynamics are interpreted by cells, we constructed a trans-omic network of insulin action in FAO hepatoma cells from three networks—a phosphorylation-dependent cellular functions regulatory network using phosphoproteomic data, a transcriptional regulatory network using phosphoproteomic and transcriptomic data, and a metabolism regulatory network using phosphoproteomic and metabolomic data. With the trans-omic regulatory network, we identified selective regulatory networks that mediate differential responses to insulin. Akt and Erk, hub molecules of insulin signaling, encode information of a wide dynamic range of dose and time of insulin. Down-regulated genes and metabolites in glycolysis had high sensitivity to insulin (fasting insulin signal); up-regulated genes and dicarboxylic acids in the TCA cycle had low sensitivity (fed insulin signal). This integrated analysis enables molecular insight into how cells interpret physiologically fed and fasting insulin signals.HighlightsWe constructed a trans-omic network of insulin action using multi-omic data.The trans-omic network integrates phosphorylation, transcription, and metabolism.We classified signaling, transcriptome, and metabolome by sensitivity to insulin.We identified fed and fasting insulin signal flow across the trans-omic network.

SUMMARYCells decode information of signaling activation at a scale of tens of minutes by downstream gene expression with a scale of hours to days, leading to cell fate decisions such as cell differentiation. However, no system identification method with such different time scales exists. Here we used compressed sensing technology and developed a system identification method using data of different time scales by recovering signals of missing time points. We measured phosphorylation of ERK and CREB, immediate early gene expression products, and mRNAs of decoder genes for neurite elongation in PC12 cell differentiation and performed system identification, revealing the input–output relationships between signaling and gene expression with sensitivity such as graded or switch-like response and with time delay and gain, representing signal transfer efficiency. We predicted and validated the identified system using pharmacological perturbation. Thus, we provide a versatile method for system identification using data with different time scales.HighlightsWe developed a system identification method using compressed sensing.This method allowed us to find a pathway using data of different time scales.We identified a selective signaling-decoding system by gene expression.We validated the identified system by pharmacological perturbation.eTOC BlurbWe describe a system identification method of molecular networks with different time-scale data using a signal recovery technique in compressed sensing.

UTokyo Research掲載「細胞内のビッグデータから大規模ネットワークの再構築に成功」URI: http://www.u-tokyo.ac.jp/ja/utokyo-research/research-news/reconstruction-of-molecular-network-from-cellular-big-data/, UTokyo Research "Reconstruction of molecular network from cellular big data" URI: http://www.u-tokyo.ac.jp/en/utokyo-research/research-news/reconstruction-of-molecular-network-from-cellular-big-data/

Secretion of insulin transiently increases after eating, resulting in a high circulating concentration. Fasting limits insulin secretion, resulting in a low concentration of insulin in the circulation. We analyzed transcriptional responses to different temporal patterns and doses of insulin in the hepatoma FAO cells and identified 13 up-regulated and 16 down-regulated insulin-responsive genes (IRGs). The up-regulated IRGs responded more rapidly than did the down-regulated IRGs to transient stepwise or pulsatile increases in insulin concentration, whereas the down-regulated IRGs were repressed at lower concentrations of insulin than those required to stimulate the up-regulated IRGs. Mathematical modeling of the insulin response as two stages-(i) insulin signaling to transcription and (ii)transcription and mRNA stability-indicated that the first stage was the more rapid stage for the down-regulated IRGs, whereas the second stage of transcription was the more rapid stage for the up-regulated IRGs. A subset of the IRGs that were up-regulated or down-regulated in the FAO cells was similarly regulated in the livers of rats injected with a single dose of insulin. Thus, not only can cells respond to insulin but they can also interpret the intensity and pattern of signal to produce distinct transcriptional responses. These results provide insight that may be useful in treating obesity and type 2 diabetes associated with aberrant insulin production or tissue responsiveness.

One of the unique characteristics of cellular signaling pathways is that a common signaling pathway can selectively regulate multiple cellular functions of a hormone; however, this selective downstream control through a common signaling pathway is poorly understood. Here we show that the insulin-dependent AKT pathway uses temporal patterns multiplexing for selective regulation of downstream molecules. Pulse and sustained insulin stimulations were simultaneously encoded into transient and sustained AKT phosphorylation, respectively. The downstream molecules, including ribosomal protein S6 kinase (S6K), glucose-6-phosphatase (G6Pase), and glycogen synthase kinase-3β (GSK3β) selectively decoded transient, sustained, and both transient and sustained AKT phosphorylation, respectively. Selective downstream decoding is mediated by the molecules' network structures and kinetics. Our results demonstrate that the AKT pathway can multiplex distinct patterns of blood insulin, such as pulse-like additional and sustained-like basal secretions, and the downstream molecules selectively decode secretion patterns of insulin.

Signaling networks are made up of limited numbers of molecules and yet can code information that controls different cellular states through temporal patterns and a combination of signaling molecules. In this study, we used a data-driven modeling approach, the Laguerre filter with partial least square regression, to describe how temporal and combinatorial patterns of signaling molecules are decoded by their downstream targets. The Laguerre filter is a time series model used to represent a nonlinear system based on Volterra series expansion. Furthermore, with this approach, each component of the Volterra series expansion is expanded by Laguerre basis functions. We combined two approaches, application of a Laguerre filter and partial least squares (PLS) regression, and applied the combined approach to analysis of a signal transduction network. We applied the Laguerre filter with PLS regression to identify input and output (IO) relationships between MAP kinases and the products of immediate early genes (IEGs). We found that Laguerre filter with PLS regression performs better than Laguerre filter with ordinary regression for the reproduction of a time series of IEGs. Analysis of the nonlinear characteristics extracted using the Laguerre filter revealed a priming effect of ERK and CREB on c-FOS induction. Specifically, we found that the effects of a first pulse of ERK enhance the subsequent effects on c-FOS induction of treatment with a second pulse of ERK, a finding consistent with prior molecular biological knowledge. The variable importance of projections and output loadings in PLS regression predicted the upstream dependency of each IEG. Thus, a Laguerre filter with partial least square regression approach appears to be a powerful method to find the processing mechanism of temporal patterns and combination of signaling molecules by their downstream gene expression.

The demand for extracting rules from high dimensional real world data is increasing in various fields. However, the possible redundancy of such data sometimes makes it difficult to obtain a good generalization ability for novel samples. To resolve this problem, we provide a scheme that reduces the effective dimensions of data by pruning redundant components for bicategorical classification based on the Bayesian framework. First, the potential of the proposed method is confirmed in ideal situations using the replica method. Unfortunately, performing the scheme exactly is computationally difficult. So, we next develop a tractable approximation algorithm, which turns out to offer nearly optimal performance in ideal cases when the system size is large. Finally, the efficacy of the developed classifier is experimentally examined for a real world problem of colon cancer classification, which shows that the developed method can be practically useful., 13 pages, 6 figures

In biological data, it is often the case that objects are described in two or more representations. In order to perform classification based on such data, we have to combine them in a certain way. In the context of kernel machines, this task amounts to mix several kernel matrices into one. In this paper, we present two ways to mix kernel matrices, where the mixing weights are optimized to minimize the cross validation error. In bacteria classification and gene function prediction experiments, our methods significantly outperformed single kernel classifiers in most cases.

Homeostatic control of blood glucose is regulated by a complex feedback loop between glucose and insulin, of which failure leads to diabetes mellitus. However, physiological and pathological nature of the feedback loop is not fully understood. We made a mathematical model of the feedback loop between glucose and insulin using time course of blood glucose and insulin during consecutive hyperglycemic and hyperinsulinemic-euglycemic clamps in 113 subjects with variety of glucose tolerance including normal glucose tolerance (NGT), impaired glucose tolerance (IGT) and type 2 diabetes mellitus (T2DM). We analyzed the correlation of the parameters in the model with the progression of glucose intolerance and the conserved relationship between parameters. The model parameters of insulin sensitivity and insulin secretion significantly declined from NGT to IGT, and from IGT to T2DM, respectively, consistent with previous clinical observations. Importantly, insulin clearance, an insulin degradation rate, significantly declined from NGT, IGT to T2DM along the progression of glucose intolerance in the mathematical model. Insulin clearance was positively correlated with a product of insulin sensitivity and secretion assessed by the clamp analysis or determined with the mathematical model. Insulin clearance was correlated negatively with postprandial glucose at 2h after oral glucose tolerance test. We also inferred a square-law between the rate constant of insulin clearance and a product of rate constants of insulin sensitivity and secretion in the model, which is also conserved among NGT, IGT and T2DM subjects. Insulin clearance shows a conserved relationship with the capacity of glucose disposal among the NGT, IGT and T2DM subjects. The decrease of insulin clearance predicts the progression of glucose intolerance.

Simple Signals? Cells process information about themselves and their surroundings through biochemical signaling pathways. Uda et al. (p. 558 8) used a recently developed cytometric method to quantitate signaling through biochemical pathways in individual rat pheochromocytoma cells responding to growth factors. The signaling pathways studied provided about 1 bit of information, or only enough for a binary (on or off) decision. In spite of the simplicity, the results showed interactions between pathways with shared components. In some cases, information carried between inputs and intermediate outputs was less than that between the input and more “downstream” outputs, indicating that information was carried through multiple paths. Similarly, in the presence of pharmacological inhibitors of one pathway, others were able to compensate to allow robust transfer of information. Thus, in spite of noise and variation in signal intensities in individual cells, robust transfer of information from the growth factors was achieved.

Cells use common signaling molecules for the selective control of downstream gene expression and cell-fate decisions. The relationship between signaling molecules and downstream gene expression and cellular phenotypes is a multiple-input and multiple-output (MIMO) system and is difficult to understand due to its complexity. For example, it has been reported that, in PC12 cells, different types of growth factors activate MAP kinases (MAPKs) including ERK, JNK, and p38, and CREB, for selective protein expression of immediate early genes (IEGs) such as c-FOS, c-JUN, EGR1, JUNB, and FOSB, leading to cell differentiation, proliferation and cell death; however, how multiple-inputs such as MAPKs and CREB regulate multiple-outputs such as expression of the IEGs and cellular phenotypes remains unclear. To address this issue, we employed a statistical method called partial least squares (PLS) regression, which involves a reduction of the dimensionality of the inputs and outputs into latent variables and a linear regression between these latent variables. We measured 1,200 data points for MAPKs and CREB as the inputs and 1,900 data points for IEGs and cellular phenotypes as the outputs, and we constructed the PLS model from these data. The PLS model highlighted the complexity of the MIMO system and growth factor-specific input-output relationships of cell-fate decisions in PC12 cells. Furthermore, to reduce the complexity, we applied a backward elimination method to the PLS regression, in which 60 input variables were reduced to 5 variables, including the phosphorylation of ERK at 10 min, CREB at 5 min and 60 min, AKT at 5 min and JNK at 30 min. The simple PLS model with only 5 input variables demonstrated a predictive ability comparable to that of the full PLS model. The 5 input variables effectively extracted the growth factor-specific simple relationships within the MIMO system in cell-fate decisions in PC12 cells.

Sensitivity is one of the hallmarks of biological and pharmacological responses. However, the principle of controlling sensitivity remains unclear. Here we theoretically analyse a simple biochemical reaction and find that the signal transfer efficiency of the transient peak amplitude attenuates depending on the strength of negative regulation. We experimentally find that many signalling pathways in various cell lines, including the Akt and ERK pathways, can be approximated by simple biochemical reactions and that the same property of the attenuation of signal transfer efficiency was observed for such pathways. Because of this property, a downstream molecule should show higher sensitivity to an activator and lower sensitivity to an inhibitor than an upstream molecule. Indeed, we experimentally verify that S6, which lies downstream of Akt, shows lower sensitivity to an epidermal growth factor receptor inhibitor than Akt. Thus, cells can control downstream sensitivity through the attenuation of signal transfer efficiency by changing the expression level of negative regulators.

A latent process involving signal transduction and gene expression is needed as a preparation step for cellular function. We previously found that nerve growth factor (NGF)-induced cell differentiation has a latent process, which is dependent on ERK activity and gene expression and required for subsequent neurite extension. A latent process can be considered a preparation step that decodes extracellular stimulus information into cellular functions; however, molecular mechanisms of this process remain unknown. We identified Metrnl, Dclk1, and Serpinb1a as latent process (LP) genes that are induced during the latent process with distinct temporal expression profiles and are required for subsequent neurite extension in PC12 cells. The LP genes showed distinct dependency on the duration of ERK activity, and they were also induced during the latent process of PACAP- and forskolin-induced cell differentiation. Regardless of neurotrophic factors, expression levels of the LP genes during the latent process (0–12 h), but not phosphorylation levels of ERK, always correlated with subsequent neurite extension length (12–24 h). Overexpression of all LP genes together, but not of each gene separately, enhanced NGF-induced neurite extension. The LP gene products showed distinct spatial localization. Thus, the LP genes appeared to be the common decoders for neurite extension length regardless of neurotrophic factors, and they may function in distinct temporal and spatial manners during the latent process. Our findings provide molecular insight into the physiological meaning of the latent process as the preparation step for decoding information for future phenotypic change.

In cellular signal transduction, the information in an external stimulus is encoded in temporal patterns in the activities of signaling molecules; for example, pulses of a stimulus may produce an increasing response or may produce pulsatile responses in the signaling molecules. Here, we show how the Akt pathway, which is involved in cell growth, specifically transmits temporal information contained in upstream signals to downstream effectors. We modeled the epidermal growth factor (EGF)-dependent Akt pathway in PC12 cells on the basis of experimental results. We obtained counterintuitive results indicating that the sizes of the peak amplitudes of receptor and downstream effector phosphorylation were decoupled; weak, sustained EGF receptor (EGFR) phosphorylation, rather than strong, transient phosphorylation, strongly induced phosphorylation of the ribosomal protein S6, a molecule downstream of Akt. Using frequency response analysis, we found that a three-component Akt pathway exhibited the property of a low-pass filter and that this property could explain decoupling of the peak amplitudes of receptor phosphorylation and that of downstream effectors. Furthermore, we found that lapatinib, an EGFR inhibitor used as an anticancer drug, converted strong, transient Akt phosphorylation into weak, sustained Akt phosphorylation, and, because of the low-pass filter characteristics of the Akt pathway, this led to stronger S6 phosphorylation than occurred in the absence of the inhibitor. Thus, an EGFR inhibitor can potentially act as a downstream activator of some effectors.