Yu Toyoshima

Decoupling of Receptor and Downstream Signals in the Akt Pathway by Its Low-Pass Filter Characteristics

The Akt pathway can serve as a signal decoupler, converting a weak receptor signal into a strong effector signal. Converting Weak into Strong Signaling pathways have complex kinetics, amplifying a signal, delaying a signal, or filtering out a signal. Fujita et al. performed kinetic analysis of the phosphorylation of three components in the epidermal growth factor receptor (EGFR) pathway—EGFR, the downstream kinase Akt, and the Akt effector ribosomal protein S6—and found that the strongest phosphorylation of the downstream effector occurred with weak, but sustained, receptor activation. Mathematical analysis indicated that Akt served as a low-pass filter to convert weak, sustained receptor signals into strong effector signals and to limit the transmission of strong, transient receptor signals. Exposure of cells to a clinically used inhibitor of the EGFR triggered a paradoxically strong phosphorylation of S6, suggesting that inhibitor data must be analyzed carefully in light of this low-pass filter characteristic. 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.

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 (GLCex). 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. GLCex was regulated by gluconeogenesis and glycolysis. These results demonstrate the selective control mechanism of glucose metabolism by temporal patterns of insulin.

Reconstruction of Insulin Signal Flow from Phosphoproteome and Metabolome Data

Cellular homeostasis is regulated by signals through multiple molecular networks that include protein phosphorylation and metabolites. However, where and when the signal flows through a network and regulates homeostasis has not been explored. We have developed a reconstruction method for the signal flow based on time-course phosphoproteome and metabolome data, using multiple databases, and have applied it to acute action of insulin, an important hormone for metabolic homeostasis. An insulin signal flows through a network, through signaling pathways that involve 13 protein kinases, 26 phosphorylated metabolic enzymes, and 35 allosteric effectors, resulting in quantitative changes in 44 metabolites. Analysis of the network reveals that insulin induces phosphorylation and activation of liver-type phosphofructokinase 1, thereby controlling a key reaction in glycolysis. We thus provide a versatile method of reconstruction of signal flow through the network using phosphoproteome and metabolome data.

Sensitivity control through attenuation of signal transfer efficiency by negative regulation of cellular signalling

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.

Automated detection and tracking of many cells by using 4D live-cell imaging data

Motivation: Automated fluorescence microscopes produce massive amounts of images observing cells, often in four dimensions of space and time. This study addresses two tasks of time-lapse imaging analyses; detection and tracking of the many imaged cells, and it is especially intended for 4D live-cell imaging of neuronal nuclei of Caenorhabditis elegans. The cells of interest appear as slightly deformed ellipsoidal forms. They are densely distributed, and move rapidly in a series of 3D images. Thus, existing tracking methods often fail because more than one tracker will follow the same target or a tracker transits from one to other of different targets during rapid moves. Results: The present method begins by performing the kernel density estimation in order to convert each 3D image into a smooth, continuous function. The cell bodies in the image are assumed to lie in the regions near the multiple local maxima of the density function. The tasks of detecting and tracking the cells are then addressed with two hill-climbing algorithms. The positions of the trackers are initialized by applying the cell-detection method to an image in the first frame. The tracking method keeps attacking them to near the local maxima in each subsequent image. To prevent the tracker from following multiple cells, we use a Markov random field (MRF) to model the spatial and temporal covariation of the cells and to maximize the image forces and the MRF-induced constraint on the trackers. The tracking procedure is demonstrated with dynamic 3D images that each contain >100 neurons of C.elegans. Availability: http://daweb.ism.ac.jp/yoshidalab/crest/ismb2014 Supplementary information: Supplementary data are available at http://daweb.ism.ac.jp/yoshidalab/crest/ismb2014 Contact: yoshidar@ism.ac.jp

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.

Accurate Automatic Detection of Densely Distributed Cell Nuclei in 3D Space

To measure the activity of neurons using whole-brain activity imaging, precise detection of each neuron or its nucleus is required. In the head region of the nematode C. elegans, the neuronal cell bodies are distributed densely in three-dimensional (3D) space. However, no existing computational methods of image analysis can separate them with sufficient accuracy. Here we propose a highly accurate segmentation method based on the curvatures of the iso-intensity surfaces. To obtain accurate positions of nuclei, we also developed a new procedure for least squares fitting with a Gaussian mixture model. Combining these methods enables accurate detection of densely distributed cell nuclei in a 3D space. The proposed method was implemented as a graphical user interface program that allows visualization and correction of the results of automatic detection. Additionally, the proposed method was applied to time-lapse 3D calcium imaging data, and most of the nuclei in the images were successfully tracked and measured.

Modulation of sensory information processing by a neuroglobin in Caenorhabditis elegans

Significance Sensory neurons encode environmental stimuli in their electrical activity and alter behavior and physiology by transmitting this information to downstream circuits. Their response properties can be characterized by tuning curves that relate stimulus parameters to neural responses. Tuning curves identify the response threshold, the stimulus features at the tuning curve peak, and high-slope regions that give maximum stimulus discrimination. Here we show that two antagonistically acting molecular oxygen sensors, a neuroglobin and a soluble guanylate cyclase, sculpt a sharp sigmoidal tuning curve in the URX oxygen sensing neurons of Caenorhabditis elegans. By combining experiments with computational modelling, we show that these changes in stimulus-encoding properties broaden C. elegans’s O2 preference. Sensory receptor neurons match their dynamic range to ecologically relevant stimulus intensities. How this tuning is achieved is poorly understood in most receptors. The roundworm Caenorhabditis elegans avoids 21% O2 and hypoxia and prefers intermediate O2 concentrations. We show how this O2 preference is sculpted by the antagonistic action of a neuroglobin and an O2-binding soluble guanylate cyclase. These putative molecular O2 sensors confer a sigmoidal O2 response curve in the URX neurons that has highest slope between 15 and 19% O2 and approaches saturation when O2 reaches 21%. In the absence of the neuroglobin, the response curve is shifted to lower O2 values and approaches saturation at 14% O2. In behavioral terms, neuroglobin signaling broadens the O2 preference of Caenorhabditis elegans while maintaining avoidance of 21% O2. A computational model of aerotaxis suggests the relationship between GLB-5–modulated URX responses and reversal behavior is sufficient to broaden O2 preference. In summary, we show that a neuroglobin can shift neural information coding leading to altered behavior. Antagonistically acting molecular sensors may represent a common mechanism to sharpen tuning of sensory neurons.

An ensemble learning approach to auto-annotation for whole-brain C. elegans imaging

Shifting from individual neuron analysis to whole-brain neural network analysis opens up new research opportunities for Caenorhabditis elegans (C. elegans). An automated data processing pipeline, including neuron detection, segmentation, tracking and annotation, will significantly improve the efficiency of analyzing whole-brain C. elegans imaging. The resulting large data sets may motivate new scientific discovery by exploiting many promising analysis tools for big data. In this study, we focus on the development of an automated annotation procedure. With only around 180 neurons in the central nervous system of a C. elegans, the annotation of each individual neuron still remains a major challenge because of the high density in space, similarity in neuron shape, unpredictable distortion of the worm’s head during motion, intrinsic variations during worm development, etc. We use an ensemble learning approach to achieve around 25% error for a test based on real experimental data. Also, we demonstrate the importance of exploring extra source of information for annotation other than the neuron positions.