Reconstruction of global regulatory network from signaling to cellular functions using phosphoproteomic data.

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.

Temporal coding of insulin action through multiplexing of the AKT pathway.

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.

Analysis of cellular signal transduction from an information theoretic approach.

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.

System identification of signaling dependent gene expression with different time-scale data.

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.

Latent process genes for cell differentiation are common decoders of neurite extension length.

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 as 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 genes that are induced during the latent process (LP) 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 hours), but not phosphorylation levels of ERK, always correlated with subsequent neurite extension length (12-24 hours). 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 appear to be the common decoders for neurite extension length regardless of neurotrophic factors, and they might 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.

Comparison of hepatic responses to glucose perturbation between healthy and obese mice based on the edge type of network structures.

Interactions between various molecular species in biological phenomena give rise to numerous networks. The investigation of these networks, including their statistical and biochemical interactions, supports a deeper understanding of biological phenomena. The clustering of nodes associated with molecular species and enrichment analysis is frequently applied to examine the biological significance of such network structures. However, these methods focus on delineating the function of a node. As such, in-depth investigations of the edges, which are the connections between the nodes, are rarely explored. In the current study, we aimed to investigate the functions of the edges rather than the nodes. To accomplish this, for each network, we categorized the edges and defined the edge type based on their biological annotations. Subsequently, we used the edge type to compare the network structures of the metabolome and transcriptome in the livers of healthy (wild-type) and obese (ob/ob) mice following oral glucose administration (OGTT). The findings demonstrate that the edge type can facilitate the characterization of the state of a network structure, thereby reducing the information available through datasets containing the OGTT response in the metabolome and transcriptome.

Four features of temporal patterns characterize similarity among individuals and molecules by glucose ingestion in humans.

Oral glucose ingestion induces systemic changes of many blood metabolites related not only to glucose, but also other metabolites such as amino acids and lipids through many blood hormones. However, the detailed temporal changes in the concentrations of comprehensive metabolites and hormones over a long time by oral glucose ingestion are uncharacterized. We measured 83 metabolites and 7 hormones in 20 healthy human subjects in response to glucose ingestion. We characterized temporal patterns of blood molecules by four features: (i) the decomposability into "amplitude" and "rate" components, (ii) the similarity of temporal patterns among individuals, (iii) the relation of molecules over time among individuals, and (iv) the similarity of temporal patterns among molecules. Glucose and glucose metabolism-related hormones indicated a rapid increase, and citrulline and lipids, which indicated a rapid decrease, returned to fasting levels faster than amino acids. Compared to glucose metabolism-related molecules and lipids, amino acids showed similar temporal patterns among individuals. The four features of temporal patterns of blood molecules by oral glucose ingestion characterize the differences among individuals and among molecules.

An extensive and dynamic trans-omic network illustrating prominent regulatory mechanisms in response to insulin in the liver.

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.

Trans-omic analysis reveals obesity-associated dysregulation of inter-organ metabolic cycles between the liver and skeletal muscle.

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.

Application of information theory in systems biology.

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.

Single-Cell Information Analysis Reveals That Skeletal Muscles Incorporate Cell-to-Cell Variability as Information Not Noise.

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.

Transomics analysis reveals allosteric and gene regulation axes for altered hepatic glucose-responsive metabolism in obesity.

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.

Logical design of oral glucose ingestion pattern minimizing blood glucose in humans.

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.

Trans-omics Impact of Thymoproteasome in Cortical Thymic Epithelial Cells.

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 Vivo Decoding Mechanisms of the Temporal Patterns of Blood Insulin by the Insulin-AKT Pathway in the Liver.

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.

Increase in hepatic and decrease in peripheral insulin clearance characterize abnormal temporal patterns of serum insulin in diabetic subjects.

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.

Trans-omic Analysis Reveals Selective Responses to Induced and Basal Insulin across Signaling, Transcriptional, and Metabolic Networks.

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.

Laguerre Filter Analysis with Partial Least Square Regression Reveals a Priming Effect of ERK and CREB on c-FOS Induction.

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.

Selective control of up-regulated and down-regulated genes by temporal patterns and doses of insulin.

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.

Glucose Homeostatic Law: Insulin Clearance Predicts the Progression of Glucose Intolerance in Humans.

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.