XB-ART-54876
Sci Rep
2016 Nov 03;6:36256. doi: 10.1038/srep36256.
Show Gene links
Show Anatomy links
Multi-phasic bi-directional chemotactic responses of the growth cone.
Naoki H
,
Nishiyama M
,
Togashi K
,
Igarashi Y
,
Hong K
,
Ishii S
.
???displayArticle.abstract???
The nerve growth cone is bi-directionally attracted and repelled by the same cue molecules depending on the situations, while other non-neural chemotactic cells usually show uni-directional attraction or repulsion toward their specific cue molecules. However, how the growth cone differs from other non-neural cells remains unclear. Toward this question, we developed a theory for describing chemotactic response based on a mathematical model of intracellular signaling of activator and inhibitor. Our theory was first able to clarify the conditions of attraction and repulsion, which are determined by balance between activator and inhibitor, and the conditions of uni- and bi-directional responses, which are determined by dose-response profiles of activator and inhibitor to the guidance cue. With biologically realistic sigmoidal dose-responses, our model predicted tri-phasic turning response depending on intracellular Ca2+ level, which was then experimentally confirmed by growth cone turning assays and Ca2+ imaging. Furthermore, we took a reverse-engineering analysis to identify balanced regulation between CaMKII (activator) and PP1 (inhibitor) and then the model performance was validated by reproducing turning assays with inhibitions of CaMKII and PP1. Thus, our study implies that the balance between activator and inhibitor underlies the multi-phasic bi-directional turning response of the growth cone.
???displayArticle.pubmedLink??? 27808115
???displayArticle.pmcLink??? PMC5093620
???displayArticle.link??? Sci Rep
Species referenced: Xenopus
Genes referenced: camk2g cdc42 npy4r prkg1 rac1
???attribute.lit??? ???displayArticles.show???
Figure 1. The mathematical growth cone model encompasses an activator and inhibitor system.(A) The schematic drawing depicting the one-dimensional model growth cone encountering an extracellular gradient of guidance cues. (B) In the growth cone, guidance cues, e.g., Netrin-1, lead to the production of cAMP and cGMP, which activate PKA and PKG, respectively, in turn inducing a Ca2+ increase. The intracellular Ca2+ up-regulates CaMKII and PP1, which function as an activator and inhibitor, respectively, of their effector to regulate growth cone motility (referred to as ‘X’ in the text and (D)), e.g., Rac1 and Cdc42. (C) Extracellular gradients of cues are transduced to intracellular gradients of Ca2+ (middle panel) in the same direction as the extracellular gradients, and the increase in the intracellular Ca2+ level that depends on the ratio of cAMP to cGMP regulates the turning direction of the growth cone, i.e., attraction or repulsion. The biased direction of the effector distribution is thought to be reversed between attraction and repulsion (lower panel). High and low levels of intracellular Ca2+ induce attraction and repulsion, respectively, as indicated by the dotted arrows. (D) An activator-inhibitor system of intracellular signaling in chemotaxis: Guidance signal (G) regulates activator (A) and inhibitor (I) in turn up- and down-regulating the effector X. G, A, I and X correspond to Ca2+, CaMKII, PP1 and Rac1/Cdc42, respectively, in the growth cone. (E) Following exposure to a gradient of the guidance cue G, the gradients of A and I are formed across the growth cone, thereby the gradient of X is also formed. A*, I* and X* represent the levels of A, I and X, respectively, at the center of the growth cone. ΔA, ΔI and ΔX represent the spatial differences of A, I and X, respectively, across the growth cone. (F) The turning angle of the growth cone is formulated as ΔX/X* based on the Weber-Fechner law. The balance between ΔA/A*and ΔI/I* governs the turning direction, which is illustrated by scales: attraction when ΔA/A*> ΔI / I* and repulsion when ΔA/A* < ΔI/I*. | |
Figure 2. The theoretical model is capable of producing both uni-directional and multi-phasic bi-directional turning responses.(A,C,E,G) Examples of activator and inhibitor activities in response to increasing guidance signal (intracellular Ca2+ in the case of the growth cone). Red and blue lines indicate the dose-responses of activator A and inhibitor I, respectively, to the guidance signal. The dose-response curves are given by Hill equations (Equations (8 and 9)). Parameter values are listed in Methods. (B,D,F,H) Turning responses (black lines) are calculated based on the dose-responses of activator and inhibitor, shown in (A,C,E,G) (see Methods). When the intracellular Ca2+ level is in the red and blue shaded regions, the growth cone exhibits attraction and repulsion, respectively. The turning response is given by the sum of the attractive factor (red dashed line: the first term in Equation (1)) and the repulsive factor (blue dashed line: the second term in Equation (1)). Each dashed line represents the ratio between the slope and amplitude of the dose-response curve in (B), (C), (E) or (G); e.g., the red dashed line in (B) was obtained from the slope (ΔA) and amplitude (A*) of the red solid line in (A). Thus, the theoretical model is able to produce a simple uni-directional attraction (A,B) or repulsion (C,D), bi-phasic bi-directional turning response (E,F), and tri-phasic bi-directional turning response (G,H). These complex responses are primarily due to the nonlinear responses of the activator and inhibitor to the guidance signal, typically observed in (G). | |
Figure 3. The growth cone exhibits multi-phasic bi-directional turning responses when the total concentration of [cAMP + cGMP] is changed at the fixed ratio of cAMP/cGMP.(A) Representative images of growth cones at the onset (0 hr) and after (1 hr) exposure to a gradient solution containing a mixture of the membrane permeable cyclic nucleotide analogues, [Sp-8-Br-cAMPS] and [8-Br-cGMP], at a fixed ratio of 9:1 at different concentrations (5 mM [4.5 cAMP + 0.5 cGMP]; 50 mM [45 cAMP + 5 cGMP]; and 100 mM [90 cAMP + 10 cGMP]) in the application micropipette. Superimposed traces on the right depict the total neurite trajectories examined over the 1 hr period, where the initial position of the growth cone is at the origin, and the original direction of growth is vertical. The arrows in both the representative images and trajectories mark the direction of the gradient. (B) The average growth cone turning angles in response to gradient of a mixture solution of [Sp-8-Br-cAMPS] and [8-Br-cGMP] with different concentrations (0, 1, 5, 10, 20, 50, and 100 mM) at the fixed ratio of 9:1 in the micropipette, which are indicated as “cNMP”. (C) The average growth cone turning angles in response to gradient of a solution of [Sp-8-Br-cAMPS] at different concentrations (5, 10, 20, 50, and 100 mM) in the micropipette. Positive and negative turning angles indicate attraction and repulsion, respectively. Error bars = s.e.m. (n) = number of growth cones examined. Significant differences compared to the control are indicated (*p < 0.05; **p < 0.01; Mann-Whitney U test). | |
Figure 4. Multi-phasic bi-directional turning depending on the concentration of cAMP/cGMP gradients is accompanied by a monotonic Ca2+ increase.(A) Representative images of Oregon Green 488 BAPTA-1 (OGB)-dextran fluorescent growth cones obtained at various times during exposure to either culture medium (control) or gradients of different concentrations of a 9:1 ratio of cAMP/cGMP (arrows). The fluorescence intensity is expressed as pseudo-colors (top; blue = low, red = high). (B) Sample traces (top) and summary (bottom) of the percentage changes in OGB-dextran fluorescence (ΔF/F) in the growth cones exposed to culture medium (control) or to different concentration gradients of cAMP/cGMP. OGB-dextran fluorescence was normalized to the fluorescent intensity of co-injected Texas Red–dextran in the growth cones before and after exposure to the cyclic nucleotide gradients (arrows in (A)). Cyclic nucleotide (cNMP) gradients: grey horizontal bars. Error bars = s.e.m. (n) = number of growth cones examined. Significant differences compared to the control are indicated (*p < 0.05; **p < 0.01; Mann-Whitney U test). (C) Left panel: Cumulative distribution of the average ΔF/F (%) at 10–20 min after exposure to cyclic nucleotide gradients as in (B). Right panel: Monotonic increase in the averages ΔF/F at 15–20 min that resulted from growth cone exposure to cyclic nucleotide gradients as in (B). The red curve is a regression function. | |
Figure 5. System identification of CaMKII and PP1 activities in the bi-directional turning response.(A) The schematic drawing depicting the regulation of CaMKII and PP1, which up- and down-regulate the effector X, respectively, in the growth cone model. (B) Simulation of the multi-phasic bi-directional turning response by a mathematical model where the parameters were identified by the reverse engineering approach (solid red line). Points and error bars represent the experimental results (the same as in Fig. 3B). Goodness of fitting was evaluated by the Pearson’s chi-square test (p = 0.78), indicating no significant difference. (C) Reverse-engineered dose-responses of CaMKII (Red line) and PP1 (blue line) to the total cAMP:cGMP concentration (cNMP). The dose-responses were normalized. The blue line is the sum of CaN- and Calpain-dependent PP1 (blue dashed and dotted lines). The estimated parameter values are listed in Methods. (D) Contributions of the activator (A; CaMKII) and the inhibitor (I; PP1) to growth cone turning (black line in (B)) depending on the cNMP level. Attractive and repulsive influences of the activator and inhibitor are represented by red and blue solid lines, respectively, which correspond to the first and second terms in Equation (1). The PP1-induced repulsion (blue solid line) is given by the sum of the CaN- and Calpain- dependent factors (dashed and dotted blue lines, respectively). | |
Figure 6. Bi-directional turning depends on the molecular switching between CaMKII and PP1.(A) Model simulation of growth cone turning responses during the inhibition of either the activator, CaMKII (red line) or the inhibitor, PP1 (blue line). Inhibition of CaMKII or PP1 activity was simulated by a 90% reduction of the estimated parameter for each concentration. (B) Growth cone turning responses were experimentally observed in response to a gradient solution that contained [Sp-8-Br-cAMPS] and [8-Br-cGMP] at a fixed ratio (9:1) at different concentrations in the presence of KN-93 (0.5 μM; CaMKII inhibitor, red line) or tautomycin (4 nM; PP1 inhibitor, blue line) in the assay medium. Significant differences compared to the control are indicated (*p < 0.05; **p < 0.01; Mann-Whitney U test). Goodness of fitting of model simulation to experiment was evaluated by p-value of the Pearson’s chi-square test (p < 0.01 with KN-93; p = 0.38 with tautomycin). |
References [+] :
Adler,
Logarithmic and power law input-output relations in sensory systems with fold-change detection.
2014, Pubmed
Adler, Logarithmic and power law input-output relations in sensory systems with fold-change detection. 2014, Pubmed
Bagnard, Spatial distributions of guidance molecules regulate chemorepulsion and chemoattraction of growth cones. 2000, Pubmed
Baier, Axon guidance by gradients of a target-derived component. 1992, Pubmed
Bauer, Developmental origin of patchy axonal connectivity in the neocortex: a computational model. 2014, Pubmed
Borisyuk, A developmental approach to predicting neuronal connectivity from small biological datasets: a gradient-based neuron growth model. 2014, Pubmed , Xenbase
Borisyuk, Stochasticity and functionality of neural systems: mathematical modelling of axon growth in the spinal cord of tadpole. 2008, Pubmed , Xenbase
Bradshaw, An ultrasensitive Ca2+/calmodulin-dependent protein kinase II-protein phosphatase 1 switch facilitates specificity in postsynaptic calcium signaling. 2003, Pubmed
Chua, Extending neurites sense the depth of the underlying topography during neuronal differentiation and contact guidance. 2014, Pubmed
Colamarino, The axonal chemoattractant netrin-1 is also a chemorepellent for trochlear motor axons. 1995, Pubmed
De Koninck, Sensitivity of CaM kinase II to the frequency of Ca2+ oscillations. 1998, Pubmed
Dickson, Molecular mechanisms of axon guidance. 2002, Pubmed
Dickson, Rho GTPases in growth cone guidance. 2001, Pubmed
Fisher, Quantitative analysis of cell motility and chemotaxis in Dictyostelium discoideum by using an image processing system and a novel chemotaxis chamber providing stationary chemical gradients. 1989, Pubmed
Fivaz, Robust neuronal symmetry breaking by Ras-triggered local positive feedback. 2008, Pubmed
Fleming, Ca2+/calmodulin-dependent protein kinase II regulates Tiam1 by reversible protein phosphorylation. 1999, Pubmed
Forbes, Calcium and cAMP levels interact to determine attraction versus repulsion in axon guidance. 2012, Pubmed
Gallo, Regulation of growth cone actin filaments by guidance cues. 2004, Pubmed
Glading, Epidermal growth factor activates m-calpain (calpain II), at least in part, by extracellular signal-regulated kinase-mediated phosphorylation. 2004, Pubmed
Gomez, The molecular basis for calcium-dependent axon pathfinding. 2006, Pubmed
Goodhill, Diffusion in axon guidance. 1997, Pubmed
Goodhill, Axon guidance: stretching gradients to the limit. 1998, Pubmed
Goodhill, Mathematical guidance for axons. 1998, Pubmed
Graupner, STDP in a bistable synapse model based on CaMKII and associated signaling pathways. 2007, Pubmed
Guan, Long-range Ca2+ signaling from growth cone to soma mediates reversal of neuronal migration induced by slit-2. 2007, Pubmed
Henley, Calcium mediates bidirectional growth cone turning induced by myelin-associated glycoprotein. 2004, Pubmed , Xenbase
Hentschel, Models of axon guidance and bundling during development. 1999, Pubmed
Hoeller, Chemotaxis in the absence of PIP3 gradients. 2007, Pubmed
Hong, A ligand-gated association between cytoplasmic domains of UNC5 and DCC family receptors converts netrin-induced growth cone attraction to repulsion. 1999, Pubmed , Xenbase
Hong, Calcium signalling in the guidance of nerve growth by netrin-1. 2000, Pubmed , Xenbase
Höpker, Growth-cone attraction to netrin-1 is converted to repulsion by laminin-1. 1999, Pubmed , Xenbase
Iglesias, Navigating through models of chemotaxis. 2008, Pubmed
Janetopoulos, Chemoattractant-induced phosphatidylinositol 3,4,5-trisphosphate accumulation is spatially amplified and adapts, independent of the actin cytoskeleton. 2004, Pubmed
Kaczmarek, Calpain cleaves and activates the TRPC5 channel to participate in semaphorin 3A-induced neuronal growth cone collapse. 2012, Pubmed
Keller, Traveling bands of chemotactic bacteria: a theoretical analysis. 1971, Pubmed
Keller, Distinct chemokinetic and chemotactic responses in neutrophil granulocytes. 1978, Pubmed
Kennedy, Axon guidance by diffusible chemoattractants: a gradient of netrin protein in the developing spinal cord. 2006, Pubmed
Kerokoski, Influence of phosphorylation of p35, an activator of cyclin-dependent kinase 5 (cdk5), on the proteolysis of p35. 2002, Pubmed
Kortholt, Ras activation and symmetry breaking during Dictyostelium chemotaxis. 2013, Pubmed
Krottje, A mathematical framework for modeling axon guidance. 2007, Pubmed
Lazova, Response rescaling in bacterial chemotaxis. 2011, Pubmed
Lee, Activation of CaMKII in single dendritic spines during long-term potentiation. 2009, Pubmed
Levchenko, Models of eukaryotic gradient sensing: application to chemotaxis of amoebae and neutrophils. 2002, Pubmed
Levit, Receptor methylation controls the magnitude of stimulus-response coupling in bacterial chemotaxis. 2002, Pubmed
Lisman, The molecular basis of CaMKII function in synaptic and behavioural memory. 2002, Pubmed
Lohof, Asymmetric modulation of cytosolic cAMP activity induces growth cone turning. 1992, Pubmed , Xenbase
Ma, Two complementary, local excitation, global inhibition mechanisms acting in parallel can explain the chemoattractant-induced regulation of PI(3,4,5)P3 response in dictyostelium cells. 2004, Pubmed
Matsumura, Intracellular calcium level required for calpain activation in a single myocardial cell. 2001, Pubmed
Meinhardt, Orientation of chemotactic cells and growth cones: models and mechanisms. 1999, Pubmed
Ming, cAMP-dependent growth cone guidance by netrin-1. 1997, Pubmed , Xenbase
Ming, Adaptation in the chemotactic guidance of nerve growth cones. 2002, Pubmed , Xenbase
Mortimer, Bayesian model predicts the response of axons to molecular gradients. 2009, Pubmed
Mortimer, Growth cone chemotaxis. 2008, Pubmed
Mortimer, Bayes-optimal chemotaxis. 2011, Pubmed
Mortimer, Axon guidance by growth-rate modulation. 2010, Pubmed
Murakoshi, Local, persistent activation of Rho GTPases during plasticity of single dendritic spines. 2011, Pubmed
Nakajima, Rectified directional sensing in long-range cell migration. 2014, Pubmed
Naoki, Flexible search for single-axon morphology during neuronal spontaneous polarization. 2011, Pubmed
Naoki, Stochastic control of spontaneous signal generation for gradient sensing in chemotaxis. 2008, Pubmed
Naoki, Local signaling with molecular diffusion as a decoder of Ca2+ signals in synaptic plasticity. 2005, Pubmed
Narang, Spontaneous polarization in eukaryotic gradient sensing: a mathematical model based on mutual inhibition of frontness and backness pathways. 2006, Pubmed
Nguyen, Regulation of protein phosphatase inhibitor-1 by cyclin-dependent kinase 5. 2007, Pubmed
Nishiyama, Membrane potential shifts caused by diffusible guidance signals direct growth-cone turning. 2008, Pubmed , Xenbase
Nishiyama, Cyclic AMP/GMP-dependent modulation of Ca2+ channels sets the polarity of nerve growth-cone turning. 2003, Pubmed , Xenbase
Nonaka, A multiphysical model of cell migration integrating reaction-diffusion, membrane and cytoskeleton. 2011, Pubmed
Okamoto, Switching characteristics of a model for biochemical-reaction networks describing autophosphorylation versus dephosphorylation of Ca2+/calmodulin-dependent protein kinase II. 2000, Pubmed
Piper, Endocytosis-dependent desensitization and protein synthesis-dependent resensitization in retinal growth cone adaptation. 2005, Pubmed , Xenbase
Postma, A diffusion-translocation model for gradient sensing by chemotactic cells. 2001, Pubmed
Roccasalvo, A hybrid computational model to predict chemotactic guidance of growth cones. 2015, Pubmed
Samuels, The origin of neuronal polarization: a model of axon formation. 1996, Pubmed
Sasaki, Localized Ras signaling at the leading edge regulates PI3K, cell polarity, and directional cell movement. 2004, Pubmed
Schoenwaelder, Calpain cleavage of focal adhesion proteins regulates the cytoskeletal attachment of integrin alphaIIbbeta3 (platelet glycoprotein IIb/IIIa) and the cellular retraction of fibrin clots. 1997, Pubmed
Segev, Generic modeling of chemotactic based self-wiring of neural networks. 2000, Pubmed
Sergi, Cell guidance on nanogratings: a computational model of the interplay between PC12 growth cones and nanostructures. 2013, Pubmed
Sergi, Deterministic control of mean alignment and elongation of neuron-like cells by grating geometry: a computational approach. 2015, Pubmed
Shelly, Local and long-range reciprocal regulation of cAMP and cGMP in axon/dendrite formation. 2010, Pubmed
Skoge, Cellular memory in eukaryotic chemotaxis. 2014, Pubmed
Skupsky, Distinguishing modes of eukaryotic gradient sensing. 2005, Pubmed
Song, cAMP-induced switching in turning direction of nerve growth cones. 1997, Pubmed , Xenbase
Song, Conversion of neuronal growth cone responses from repulsion to attraction by cyclic nucleotides. 1998, Pubmed , Xenbase
Song, Dictyostelium discoideum chemotaxis: threshold for directed motion. 2006, Pubmed
Stemmer, Dual calcium ion regulation of calcineurin by calmodulin and calcineurin B. 1994, Pubmed
Subramanian, A mechanistic model for eukaryotic gradient sensing: spontaneous and induced phosphoinositide polarization. 2004, Pubmed
Togashi, Cyclic GMP-gated CNG channels function in Sema3A-induced growth cone repulsion. 2008, Pubmed , Xenbase
Tojima, Attractive axon guidance involves asymmetric membrane transport and exocytosis in the growth cone. 2007, Pubmed
Tojima, Asymmetric clathrin-mediated endocytosis drives repulsive growth cone guidance. 2010, Pubmed
Tomba, Tuning the adhesive geometry of neurons: length and polarity control. 2014, Pubmed
Toriyama, A diffusion-based neurite length-sensing mechanism involved in neuronal symmetry breaking. 2010, Pubmed
Turing, The chemical basis of morphogenesis. 1953. 1990, Pubmed
Urakubo, Requirement of an allosteric kinetics of NMDA receptors for spike timing-dependent plasticity. 2008, Pubmed
Urakubo, In vitro reconstitution of a CaMKII memory switch by an NMDA receptor-derived peptide. 2014, Pubmed
Van Haastert, Sensory adaptation of Dictyostelium discoideum cells to chemotactic signals. 1983, Pubmed
Walter, A common denominator of growth cone guidance and collapse? 1990, Pubmed
Wen, A CaMKII/calcineurin switch controls the direction of Ca(2+)-dependent growth cone guidance. 2004, Pubmed , Xenbase
Xia, The role of calmodulin as a signal integrator for synaptic plasticity. 2005, Pubmed
Xiong, Cells navigate with a local-excitation, global-inhibition-biased excitable network. 2010, Pubmed
Xu, Adaptation is not required to explain the long-term response of axons to molecular gradients. 2005, Pubmed
Yamao, Distinct predictive performance of Rac1 and Cdc42 in cell migration. 2015, Pubmed
Zhabotinsky, Bistability in the Ca(2+)/calmodulin-dependent protein kinase-phosphatase system. 2000, Pubmed
Zheng, Turning of nerve growth cones induced by localized increases in intracellular calcium ions. 2000, Pubmed , Xenbase
Zheng, Turning of nerve growth cones induced by neurotransmitters. 1994, Pubmed , Xenbase
Zubler, A framework for modeling the growth and development of neurons and networks. 2009, Pubmed
van Ooyen, Using theoretical models to analyse neural development. 2011, Pubmed