Quantitative metabolic fluxes regulated by trans-omic networks.

Cells change their metabolism in response to internal and external conditions by regulating the trans-omic network, which is a global biochemical network with multiple omic layers. Metabolic flux is a direct measure of the activity of a metabolic reaction that provides valuable information for understanding complex trans-omic networks. Over the past decades, techniques to determine metabolic fluxes, including 13C-metabolic flux analysis (13C-MFA), flux balance analysis (FBA), and kinetic modeling, have been developed. Recent studies that acquire quantitative metabolic flux and multi-omic data have greatly advanced the quantitative understanding and prediction of metabolism-centric trans-omic networks. In this review, we present an overview of 13C-MFA, FBA, and kinetic modeling as the main techniques to determine quantitative metabolic fluxes, and discuss their advantages and disadvantages. We also introduce case studies with the aim of understanding complex metabolism-centric trans-omic networks based on the determination of metabolic fluxes.

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.

NMDAR-Mediated Ca2+ Increase Shows Robust Information Transfer in Dendritic Spines.

A dendritic spine is a small structure on the dendrites of a neuron that processes input timing information from other neurons. Tens of thousands of spines are present on a neuron. Why are spines so small and many? We addressed this issue by using the stochastic simulation model of N-methyl-D-aspartate receptor (NMDAR)-mediated Ca2+ increase. NMDAR-mediated Ca2+ increase codes the input timing information between prespiking and postspiking. We examined how much the input timing information is encoded by Ca2+ increase against presynaptic fluctuation. We found that the input timing information encoded in the cell volume (103µm3) largely decreases against the presynaptic fluctuation, whereas that in the spine volume (10-1µm3) slightly decreases. Therefore, the input timing information encoded in the spine volume is more robust against presynaptic fluctuation than that in the cell volume. We further examined the mechanism of the robust information transfer in the spine volume. We demonstrated that the condition for the robustness is that the stochastic NMDAR-mediated Ca2+ increase (intrinsic noise) becomes much larger than the presynaptic fluctuation (extrinsic noise). When the presynaptic fluctuation is large, the condition is satisfied in the spine volume but not in the cell volume. Moreover, we compared the input timing information encoded in many small spines with that encoded in a single large spine. We found that the input timing information encoded in many small spines are larger than that in a single large spine when presynaptic fluctuation is large because of their robustness. Thus, robustness is a functional reason why dendritic spines are so small and many.

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.

Automatic Quantitative Segmentation of Myotubes Reveals Single-cell Dynamics of S6 Kinase Activation.

Automatic cell segmentation is a powerful method for quantifying signaling dynamics at single-cell resolution in live cell fluorescence imaging. Segmentation methods for mononuclear and round shape cells have been developed extensively. However, a segmentation method for elongated polynuclear cells, such as differentiated C2C12 myotubes, has yet to be developed. In addition, myotubes are surrounded by undifferentiated reserve cells, making it difficult to identify background regions and subsequent quantification. Here we developed an automatic quantitative segmentation method for myotubes using watershed segmentation of summed binary images and a two-component Gaussian mixture model. We used time-lapse fluorescence images of differentiated C2C12 cells stably expressing Eevee-S6K, a fluorescence resonance energy transfer (FRET) biosensor of S6 kinase (S6K). Summation of binary images enhanced the contrast between myotubes and reserve cells, permitting detection of a myotube and a myotube center. Using a myotube center instead of a nucleus, individual myotubes could be detected automatically by watershed segmentation. In addition, a background correction using the two-component Gaussian mixture model permitted automatic signal intensity quantification in individual myotubes. Thus, we provide an automatic quantitative segmentation method by combining automatic myotube detection and background correction. Furthermore, this method allowed us to quantify S6K activity in individual myotubes, demonstrating that some of the temporal properties of S6K activity such as peak time and half-life of adaptation show different dose-dependent changes of insulin between cell population and individuals.Key words: time lapse images, cell segmentation, fluorescence resonance energy transfer, C2C12, myotube.

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.

Small-Volume Effect Enables Robust, Sensitive, and Efficient Information Transfer in the Spine.

Why is the spine of a neuron so small that it can contain only small numbers of molecules and reactions inevitably become stochastic? We previously showed that, despite such noisy conditions, the spine exhibits robust, sensitive, and efficient features of information transfer using the probability of Ca2+ increase; however, the mechanisms are unknown. In this study, we show that the small volume effect enables robust, sensitive, and efficient information transfer in the spine volume, but not in the cell volume. In the spine volume, the intrinsic noise in reactions becomes larger than the extrinsic noise of input, resulting in robust information transfer despite input fluctuation. In the spine volume, stochasticity makes the Ca2+ increase occur with a lower intensity of input, causing higher sensitivity to lower intensity of input. The volume-dependency of information transfer increases its efficiency in the spine volume. Thus, we propose that the small-volume effect is the functional reason why the spine has to be so small.

In vitro reconstitution of a CaMKII memory switch by an NMDA receptor-derived peptide.

Ca(2+)/Calmodulin-dependent protein kinase II (CaMKII) has been shown to play a major role in establishing memories through complex molecular interactions including phosphorylation of multiple synaptic targets. However, it is still controversial whether CaMKII itself serves as a molecular memory because of a lack of direct evidence. Here, we show that a single holoenzyme of CaMKII per se serves as an erasable molecular memory switch. We reconstituted Ca(2+)/Calmodulin-dependent CaMKII autophosphorylation in the presence of protein phosphatase 1 in vitro, and found that CaMKII phosphorylation shows a switch-like response with history dependence (hysteresis) only in the presence of an N-methyl-D-aspartate receptor-derived peptide. This hysteresis is Ca(2+) and protein phosphatase 1 concentration-dependent, indicating that the CaMKII memory switch is not simply caused by an N-methyl-D-aspartate receptor-derived peptide lock of CaMKII in an active conformation. Mutation of a phosphorylation site of the peptide shifted the Ca(2+) range of hysteresis. These functions may be crucial for induction and maintenance of long-term synaptic plasticity at hippocampal synapses.

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.

Physical hydrogels composed of polyampholytes demonstrate high toughness and viscoelasticity.

Hydrogels attract great attention as biomaterials as a result of their soft and wet nature, similar to that of biological tissues. Recent inventions of several tough hydrogels show their potential as structural biomaterials, such as cartilage. Any given application, however, requires a combination of mechanical properties including stiffness, strength, toughness, damping, fatigue resistance and self-healing, along with biocompatibility. This combination is rarely realized. Here, we report that polyampholytes, polymers bearing randomly dispersed cationic and anionic repeat groups, form tough and viscoelastic hydrogels with multiple mechanical properties. The randomness makes ionic bonds of a wide distribution of strength. The strong bonds serve as permanent crosslinks, imparting elasticity, whereas the weak bonds reversibly break and re-form, dissipating energy. These physical hydrogels of supramolecular structure can be tuned to change multiple mechanical properties over wide ranges by using diverse ionic combinations. This polyampholyte approach is synthetically simple and dramatically increases the choice of tough hydrogels for applications.

The selective control of glycolysis, gluconeogenesis and glycogenesis by temporal insulin patterns.

Insulin governs systemic glucose metabolism, including glycolysis, gluconeogenesis and glycogenesis, through temporal change and absolute concentration. However, how insulin-signalling pathway selectively regulates glycolysis, gluconeogenesis and glycogenesis remains to be elucidated. To address this issue, we experimentally measured metabolites in glucose metabolism in response to insulin. Step stimulation of insulin induced transient response of glycolysis and glycogenesis, and sustained response of gluconeogenesis and extracellular glucose concentration (GLC(ex)). Based on the experimental results, we constructed a simple computational model that characterises response of insulin-signalling-dependent glucose metabolism. The model revealed that the network motifs of glycolysis and glycogenesis pathways constitute a feedforward (FF) with substrate depletion and incoherent feedforward loop (iFFL), respectively, enabling glycolysis and glycogenesis responsive to temporal changes of insulin rather than its absolute concentration. In contrast, the network motifs of gluconeogenesis pathway constituted a FF inhibition, enabling gluconeogenesis responsive to absolute concentration of insulin regardless of its temporal patterns. GLC(ex) was regulated by gluconeogenesis and glycolysis. These results demonstrate the selective control mechanism of glucose metabolism by temporal patterns of insulin.

Mathematical Models of the Effect of Glucagon on Glycemia in Individuals With Type 2 Diabetes Treated With Dapagliflozin.

Context: Sodium-glucose cotransporter 2 (SGLT2) inhibitors lower blood glucose levels by promoting urinary glucose excretion, but their overall effects on hormonal and metabolic status remain unclear. Objective: We here investigated the roles of insulin and glucagon in the regulation of glycemia in individuals treated with an SGLT2 inhibitor using mathematical model analysis. Methods: Hyperinsulinemic-euglycemic clamp and oral glucose tolerance tests were performed in 68 individuals with type 2 diabetes treated with the SGLT2 inhibitor dapagliflozin. Data previously obtained from such tests in 120 subjects with various levels of glucose tolerance and not treated with an SGLT2 inhibitor were examined as a control. Mathematical models of the feedback loops connecting glucose and insulin (GI model) or glucose, insulin, and glucagon (GIG model) were generated. Results: Analysis with the GI model revealed that the disposition index/clearance, which is defined as the product of insulin sensitivity and insulin secretion divided by the square of insulin clearance and represents the glucose-handling ability of insulin, was significantly correlated with glycemia in subjects not taking an SGLT2 inhibitor but not in those taking dapagliflozin. Analysis with the GIG model revealed that a metric defined as the product of glucagon sensitivity and glucagon secretion divided by glucagon clearance (designated production index/clearance) was significantly correlated with blood glucose level in subjects treated with dapagliflozin. Conclusion: Treatment with an SGLT2 inhibitor alters the relation between insulin effect and blood glucose concentration, and glucagon effect may account for variation in glycemia among individuals treated with such drugs.

High-Performance Genetically Encoded Green Fluorescent Biosensors for Intracellular l-Lactate.

l-Lactate is a monocarboxylate produced during the process of cellular glycolysis and has long generally been considered a waste product. However, studies in recent decades have provided new perspectives on the physiological roles of l-lactate as a major energy substrate and a signaling molecule. To enable further investigations of the physiological roles of l-lactate, we have developed a series of high-performance (ΔF/F = 15 to 30 in vitro), intensiometric, genetically encoded green fluorescent protein (GFP)-based intracellular l-lactate biosensors with a range of affinities. We evaluated these biosensors in cultured cells and demonstrated their application in an ex vivo preparation of Drosophila brain tissue. Using these biosensors, we were able to detect glycolytic oscillations, which we analyzed and mathematically modeled.

Trans-omic analysis reveals opposite metabolic dysregulation between feeding and fasting in liver associated with obesity.

Dysregulation of liver metabolism associated with obesity during feeding and fasting leads to the breakdown of metabolic homeostasis. However, the underlying mechanism remains unknown. Here, we measured multi-omics data in the liver of wild-type and leptin-deficient obese (ob/ob) mice at ad libitum feeding and constructed a differential regulatory trans-omic network of metabolic reactions. We compared the trans-omic network at feeding with that at 16 h fasting constructed in our previous study. Intermediate metabolites in glycolytic and nucleotide metabolism decreased in ob/ob mice at feeding but increased at fasting. Allosteric regulation reversely shifted between feeding and fasting, generally showing activation at feeding while inhibition at fasting in ob/ob mice. Transcriptional regulation was similar between feeding and fasting, generally showing inhibiting transcription factor regulations and activating enzyme protein regulations in ob/ob mice. The opposite metabolic dysregulation between feeding and fasting characterizes breakdown of metabolic homeostasis associated with obesity.

Raf kinase inhibitory protein is required for cerebellar long-term synaptic depression by mediating PKC-dependent MAPK activation.

It was demonstrated previously that a positive feedback loop, including protein kinase C (PKC) and mitogen-activated protein kinase (MAPK), is required for the gradual expression of cerebellar long-term depression (LTD). PKC and MAPK are mutually activated in this loop. MAPK-dependent PKC activation is likely to be mediated by phospholipase A2. On the other hand, it is not clear how PKC activates MAPK. Therefore, the entire picture of this loop was not fully understood. We here test the hypothesis that this loop is completed by the PKC substrate, Raf kinase inhibitory protein (RKIP). To test this hypothesis, we used a mutant form of RKIP that is not phosphorylated by PKC and thus constitutively inhibits Raf-1 and MEK, upstream kinases of MAPK. When this RKIP mutant was introduced into Purkinje cells of mouse cerebellar slices through patch-clamp electrodes, LTD was blocked, while wild-type (WT) RKIP had no effect on LTD. Physiological epistasis experiments demonstrated that RKIP works downstream of PKC and upstream of MAPK during LTD induction. Furthermore, biochemical analyses demonstrated that endogenous RKIP dissociates from Raf-1 and MEK during LTD induction in a PKC-dependent manner, suggesting that RKIP binding-dependent inhibition of Raf-1 and MEK is removed upon LTD induction. We therefore conclude that PKC-dependent regulation of RKIP leads to MAPK activation, with RKIP completing the positive feedback loop that is required for LTD.

Features extracted using tensor decomposition reflect the biological features of the temporal patterns of human blood multimodal metabolome.

High-throughput omics technologies have enabled the profiling of entire biological systems. For the biological interpretation of such omics data, two analyses, hypothesis- and data-driven analyses including tensor decomposition, have been used. Both analyses have their own advantages and disadvantages and are mutually complementary; however, a direct comparison of these two analyses for omics data is poorly examined.We applied tensor decomposition (TD) to a dataset representing changes in the concentrations of 562 blood molecules at 14 time points in 20 healthy human subjects after ingestion of 75 g oral glucose. We characterized each molecule by individual dependence (constant or variable) and time dependence (later peak or early peak). Three of the four features extracted by TD were characterized by our previous hypothesis-driven study, indicating that TD can extract some of the same features obtained by hypothesis-driven analysis in a non-biased manner. In contrast to the years taken for our previous hypothesis-driven analysis, the data-driven analysis in this study took days, indicating that TD can extract biological features in a non-biased manner without the time-consuming process of hypothesis generation.

DI/cle, a Measure Consisting of Insulin Sensitivity, Secretion, and Clearance, Captures Diabetic States.

CONTEXT: Insulin clearance is implicated in regulation of glucose homeostasis independently of insulin sensitivity and insulin secretion. OBJECTIVE: To understand the relation between blood glucose and insulin sensitivity, secretion, and clearance. METHODS: We performed a hyperglycemic clamp, a hyperinsulinemic-euglycemic clamp, and an oral glucose tolerance test (OGTT) in 47, 16, and 49 subjects with normal glucose tolerance (NGT), impaired glucose tolerance (IGT), and type 2 diabetes mellitus (T2DM), respectively. Mathematical analyses were retrospectively performed on this dataset. RESULTS: The disposition index (DI), defined as the product of insulin sensitivity and secretion, showed a weak correlation with blood glucose levels, especially in IGT (r = 0.04; 95% CI, -0.63 to 0.44). However, an equation relating DI, insulin clearance, and blood glucose levels was well conserved regardless of the extent of glucose intolerance. As a measure of the effect of insulin, we developed an index, designated disposition index/clearance, (DI/cle) that is based on this equation and corresponds to DI divided by the square of insulin clearance. DI/cle was not impaired in IGT compared with NGT, possibly as a result of a decrease in insulin clearance in response to a reduction in DI, whereas it was impaired in T2DM relative to IGT. Moreover, DI/cle estimated from a hyperinsulinemic-euglycemic clamp, OGTT, or a fasting blood test were significantly correlated with that estimated from 2 clamp tests (r = 0.52; 95% CI, 0.37 to 0.64, r = 0.43; 95% CI, 0.24 to 0.58, r = 0.54; 95% CI, 0.38 to 0.68, respectively). CONCLUSION: DI/cle can serve as a new indicator for the trajectory of changes in glucose tolerance.

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.

Optogenetic decoding of Akt2-regulated metabolic signaling pathways in skeletal muscle cells using transomics analysis.

Insulin regulates various cellular metabolic processes by activating specific isoforms of the Akt family of kinases. Here, we elucidated metabolic pathways that are regulated in an Akt2-dependent manner. We constructed a transomics network by quantifying phosphorylated Akt substrates, metabolites, and transcripts in C2C12 skeletal muscle cells with acute, optogenetically induced activation of Akt2. We found that Akt2-specific activation predominantly affected Akt substrate phosphorylation and metabolite regulation rather than transcript regulation. The transomics network revealed that Akt2 regulated the lower glycolysis pathway and nucleotide metabolism and cooperated with Akt2-independent signaling to promote the rate-limiting steps in these processes, such as the first step of glycolysis, glucose uptake, and the activation of the pyrimidine metabolic enzyme CAD. Together, our findings reveal the mechanism of Akt2-dependent metabolic pathway regulation, paving the way for Akt2-targeting therapeutics in diabetes and metabolic disorders.