Muscarinic Acetylcholine Type 1 Receptor Activity Constrains Neurite Outgrowth by Inhibiting Microtubule Polymerization and Mitochondrial Trafficking in Adult Sensory Neurons

Abstract

The muscarinic acetylcholine type 1 receptor (M1R) is a metabotropic G protein-coupled receptor. Knockout of M1R or exposure to selective or specific receptor antagonists elevates neurite outgrowth in adult sensory neurons and is therapeutic in diverse models of peripheral neuropathy. We tested the hypothesis that endogenous M1R activation constrained neurite outgrowth via a negative impact on the cytoskeleton and subsequent mitochondrial trafficking. We overexpressed M1R in primary cultures of adult rat sensory neurons and cell lines and studied the physiological and molecular consequences related to regulation of cytoskeletal/mitochondrial dynamics and neurite outgrowth. In adult primary neurons, overexpression of M1R caused disruption of the tubulin, but not actin, cytoskeleton and significantly reduced neurite outgrowth. Over-expression of a M1R-DREADD mutant comparatively increased neurite outgrowth suggesting that acetylcholine released from cultured neurons interacts with M1R to suppress neurite outgrowth. M1R-dependent constraint on neurite outgrowth was removed by selective (pirenzepine) or specific (muscarinic toxin 7) M1R antagonists. M1R-dependent disruption of the cytoskeleton also diminished mitochondrial abundance and trafficking in distal neurites, a disorder that was also rescued by pirenzepine or muscarinic toxin 7. M1R activation modulated cytoskeletal dynamics through activation of the G protein (Gα13) that inhibited tubulin polymerization and thus reduced neurite outgrowth. Our study provides a novel mechanism of M1R control of Gα13 protein-dependent modulation of the tubulin cytoskeleton, mitochondrial trafficking and neurite outgrowth in axons of adult sensory neurons. This novel pathway could be harnessed to treat dying-back neuropathies since anti-muscarinic drugs are currently utilized for other clinical conditions.

Introduction

Muscarinic acetylcholine receptors constitute a sub-family of G protein-coupled receptors (GPCRs) that act as metabotropic activators of the neurotransmitter acetylcholine (ACh). Five distinct subtypes have been identified (M1-M5), based on their G-protein coupling preferences (Wess, 1996). Downstream pathways activated include phospholipase C, inositol triphosphate (IP3), cyclic adenosine monophosphate (cAMP) and altered calcium homeostasis (Eglen, 2005; Wess et al., 2007; Kruse et al., 2014). In addition, these GPCRs modulate the cytoskeleton through trimeric G protein signaling (Kapitein and Hoogenraad, 2015). For example, α and βγ subunits of heterotrimeric G proteins modulate microtubule assembly (Roychowdhury and Rasenick, 2008; Schappi et al., 2014). Activated Gα, acts as a GTPase activating protein (GAP) and increases microtubule disassembly by activating the intrinsic GTPase activity of tubulin (Roychowdhury et al., 1999).

The muscarinic acetylcholine type 1 receptor (M1R) is widely expressed in the central nervous system (CNS) (Levey, 1993; Wess et al., 2003; Jiang et al., 2014) and peripheral nervous system (PNS) (Bernardini et al., 1999; Tata et al., 2000b). Membranes isolated from hippocampus and cortex of M1R knockout (KO) mice showed a significant decrease in GTPγ-S loading to the Gα-q/11 G protein upon agonist stimulation (Felder et al., 2001). In cortical neuron cultures obtained from M1R KO mice, carbachol-stimulated phosphoinositide hydrolysis was reduced by 60% compared with wild type (Bymaster et al., 2003). In addition, phosphorylation of extracellular signal-regulated kinase (ERK) was eliminated in pyramidal neurons of hippocampal slices or cortical cultures derived from M1R KO mice (Berkeley et al., 2001; Hamilton and Nathanson, 2001).

In sympathetic neurons, ACh activation of M1R mobilizes internal Ca²⁺ stores leading to closure of M-type K⁺ channels (Kv7 subtypes) and enhancement of slow depolarization and discharge (Delmas and Brown, 2005; Brown and Passmore, 2009). In embryonic neurons, ACh modulates neurite outgrowth in a positive or negative manner based upon context (Tata et al., 2000a, 2003; Bernardini et al., 2004; Yang and Kunes, 2004). Furthermore, both adult sensory dorsal root ganglia (DRG) neurons and epidermal keratinocytes synthesize and secrete ACh (Bernardini et al., 1999; Khan et al., 2002; Nguyen et al., 2004; Grando et al., 2006; Schlereth et al., 2006; Corsetti et al., 2012). Adult rat sensory neurons of the DRG express a peripheral form of ChAT (pChAT), exhibit ChAT activity, have low AChE activity and express multiple muscarinic receptors including M1R (Bernardini et al., 1999; Tata et al., 2000b; Bellier and Kimura, 2007; Hanada et al., 2013).

We have recently reported that selective or specific antagonists of M1R elevated neurite outgrowth and augmented mitochondrial function in adult sensory neurons (Calcutt et al., 2017). These drugs also afforded protection against several different forms of peripheral neuropathy. However, the mechanism of M1R antagonist-driven neurite outgrowth and neuroprotection has not been studied in detail. Mitochondrial oxidative phosphorylation is the main mechanism providing ATP to power neuronal activities such as production of presynaptic action potentials, neurotransmitter release, postsynaptic currents and postsynaptic action potentials (Hall et al., 2012). Mitochondria are known to concentrate in regions of active signaling and high metabolic demand (Chen and Chan, 2006; Mironov, 2007; Verburg and Hollenbeck, 2008). This substantial energy demand at the nerve ending or synapse implies that neurons must have a mechanism to maintain microtubules to augment mitochondrial trafficking upon demand (Sheng and Cai, 2012; Schwarz, 2013).

In the present study we manipulated M1R expression/function in adult DRG sensory neurons and related cell lines and studied the cellular phenotypes and molecular consequences. Specifically, we tested the hypothesis that the M1R regulates the tubulin cytoskeleton, G-protein recruitment (Gα13 sub-type) and mitochondrial trafficking. We identified that excessive cholinergic signaling triggered tubulin destabilization through over-activation of Gα13 proteins. Further, we studied the ability of specific (muscarinic toxin 7: MT7) or selective (pirenzepine) M1R antagonists to ameliorate the endogenous and M1R overexpression-induced neuronal phenotypes that primarily result in a constraint on neurite outgrowth.

Materials and Methods

All animal procedures followed guidelines of University of Manitoba Animal Care Committee using Canadian Council of Animal Care rules or of the Institutional Animal Care and Use Committee at UCSD.

Cell Culture

Dorsal root ganglia from adult male Sprague-Dawley rats were dissected and dissociated using previously described methods (Akude et al., 2011; Roy Chowdhury et al., 2012; Saleh et al., 2013). All animal protocols carefully followed the Canadian Committee on Animal Care (CCAC) guidelines. Neurons were cultured in defined Hams F12 media containing 10 mM D-Glucose (N4888, Sigma, St. Louis, MO, United States) supplemented with modified Bottenstein’s N2 additives without insulin (0.1 mg/ml transferrin, 20 nM progesterone, 100 μM putrescine, 30 nM sodium selenite, 0.1 mg/ml BSA; all additives were from Sigma, St. Louis, MO, United States). In all experiments, the media was also supplemented with 0.146 g/L L-glutamine, a low dose or high dose cocktail of neurotrophic factors (Low dose = 0.1 ng/ml NGF, 1.0 ng/ml GDNF and 1 ng/ml NT-3, High dose = 1 ng/ml NGF, 10 ng/ml GDNF, 10 ng/ml NT3 – all from Promega, Madison, WI, United States), 0.1 nM insulin and 1X antibiotic antimycotic solution (A5955, Sigma). Cultures were treated with 100 nM MT7 (M-200, Alomone Labs, Jerusalem, Israel) or 1 μM pirenzepine (P7412, Sigma).

HEK293 and HTLA cells were cultivated in Dulbecco’s modified Eagle’s medium (DMEM) supplemented with 10% heat inactivated FBS. The β-arrestin null (ARRB1 and ARRB2) and Gα12/13 (GNA12 and GNA13) null HEK293 cells were obtained from the laboratory of Dr. Asuka Inoue, Tohoku University, Japan. HTLA cells were provided by Prof Bryan Roth, University of North Carolina, United States.

Cloning, Transfection and siRNA Based Gene Silencing

Total mRNA was extracted from adult rat DRGs using Trizol reagent and used for amplifying the full length M1R cDNA using the following primer sets: F: 5′-ATGAACACCTCAGTGCCCCCTGC-3′ and R: 5′-TTAGCATTGGCGGGAGGGGGTG-3′. The cDNA was cloned in the pEGFP-C1 vector (Clontech, now Takara Bio United States, Inc., Mountain View, CA, United States) in the XhoI and SacII restriction sites. In addition, the cDNA was also cloned in the pHTN Halo-Tag CMV-neo Vector (Promega, Madison, WI, United States) at the SacII and Not1 restriction sites. The plasmids were transiently transfected in freshly dissociated sensory neurons using the rat neuron nucleofection kit (VPG-1003, Amaxa, Lonza Inc., Allendale, NJ, United States) and Amaxa nucleofector-II apparatus (program 0-003) and cultured in poly L-Ornithine (P8638, Sigma) and laminin coated μ-Plate-24 well (Ibidi United States, Inc., Madison, WI, United States). The human M1-DREADD construct was obtained from Dr. Arthur Christopoulos, Monash University, Australia and sub-cloned in pEGFP-C1 vector (Abdul-Ridha et al., 2013). The rat Gα13(GNA13) was knocked down using a cocktail of three siRNAs targeted to exon 2 (AGTATCTTCCTGCTATAAGAGCC) and exon 4 (CTACATTCCGTCACAGCAAGATA and CATCAAAGACTATTTCCTAGAAT), respectively. The siRNAs were transfected in to primary sensory neurons using Amaxa transfection reagent.

Quantification of Neurite Outgrowth

The transgene transfected neurons were cultured for 48 h and then cells were fixed in 4% paraformaldehyde for 10 min and immunostained using monoclonal anti-β-tubulin III antibody. The neurons were also stained with Hoechst for nuclear staining. The neurons were imaged in an unbiased manner using a Cellomics Arrayscan-VTI high content screening (HCS) Reader (Thermo Fisher Scientific, Waltham, MA United States) and total neurite outgrowth per neuron was measured by Neuronal Profiling V4.1 software. The automated HCS reader provided a bias-free objective analysis of neurite outgrowth.

Confocal Microscopic Image Acquisition and Analysis to Determine Mitochondrial Volume and Trafficking

Mitochondrial trafficking in GFP or GFP-M1R overexpressing neurons was monitored using LSM510 confocal live cell imaging and involved co-transfection of sensory neurons with GFP/GFP-M1R and DsRed2Mito7 plasmids (Addgene plasmid #55838, a gift from Dr. Michael W. Davidson, Florida State University), respectively. The DsRed2Mito7 consists of a mitochondrial targeting sequence from subunit VIII of human cytochrome C oxidase which is placed before the Red fluorescence protein and the resultant fusion protein selectively accumulates inside the mitochondria (Van Kuilenburg et al., 1988). The transfected cells were live imaged at 10 s interval for 80 time frames (∼13 min). The time lapse images were used to generate kymographs using the ImageJ make kymograph plugin (Schneider et al., 2012). In each kymograph, the x-axis represents the position along the length of the axon and the y-axis represents time. Vertical lines indicate stationary mitochondria with no displacement during the time elapsed and diagonal lines represent moving mitochondria and their direction. Their velocity is reflected in the slope of the line. In addition we used Fiji (Schindelin et al., 2012) based MTrackJ to determine mitochondrial velocity (μm/sec) in the neurites (Meijering et al., 2012). The volume of mitochondria in the neurites was calculated by using Image J analyze particles plugin and expressed in μm per neurite length (Schindelin et al., 2012).

Western Blotting and Immune-Detection

Relative quantification of proteins was done by SDS-PAGE separation of total proteins followed by transfer to nitrocellulose membrane and immunoblotting based detection using HRP-conjugated secondary antibodies. The immunoblots were imaged in Bio-Rad Chemidoc system (Bio-Rad Laboratories Ltd., Mississauga, ON, Canada). Table 1 summarizes all the primary antibodies used in this study. The cell lysates were prepared in 1X RIPA lysis and extraction buffer (Cat No: 89900, Thermo Fisher Scientific) supplemented with 1X Halt protease and phosphatase inhibitor cocktail (Cat No: 78441, Thermo Fisher Scientific).

TABLE 1

TABLE 1. List of antibodies.

Polymerized Tubulin Quantification

The polymerized tubulin in M1R overexpressed cells (40–50% transfection efficiency) was quantified by methods described previously (Butcher et al., 2016). Briefly, the soluble fraction of tubulin was first removed by lysing the cells in a microtubule stabilizing buffer (MSB) containing 50% glycerol, 5 mM MgCl2, 0.1 mM EGTA, 0.3 mM guanosine triphosphate (grade II-S, Sigma Chemical Co.), and 10 mM sodium phosphate, pH 6.8 (Beertsen et al., 1982). Cells were harvested with a rubber scraper in MSB, homogenized, and centrifuged at 26,000 × g in a Sorvall RC2-B centrifuge (17,000 rpm in rotor SS-34 in 1.0 ml tubes; Dupont Instruments-Sorvall Biomedical Div., Dupont Co., Newtown, CT, United States) at 20°C for 30 min. Supernatants containing the soluble tubulin fraction were removed and the pellet containing polymerized tubulin and other cytoskeletal protein was assessed by immunoblotting using anti-actin and anti-β tubulin III antibodies.

Halo-M1R Pull Down and Blue-Native Polyacrylamide Gel Electrophoresis (BN-PAGE)

SH-SY5Y human neuroblastoma cells (provided by Dr. Jun-Feng Wang, University of Manitoba) were grown in dulbecco’s modified eagle medium: nutrient mixture F-12 (DMEM/F12, Thermo Fisher) supplemented with 10% fetal bovine serum (FBS, Thermo Fisher). Halo-M1R plasmid was transiently expressed in SH-SY5Y cells that were treated with 100 nM MT7 or 1 μM pirenzepine for 1 h. Halo-M1R was then pulled down using Halo-link resin by overnight incubation at 4°C (Promega Corporation, Madison, WI, United States) as per manufacturer’s instruction and the pull-down product was cleaved using TEV-protease (Promega) overnight on a constant rotating shaking platform. The cleaved fraction was resolved in SDS-PAGE and immunoblotted using anti-M1R, anti Gγ2/3/4/7 and anti-Gα12/13 antibodies. The BN-PAGE was performed as described previously (Sabbir et al., 2016). Polymeric tubulins in Gα12/13 null and native HEK293 cells were separated by BN-PAGE using microtubule stabilizing native cell lysis buffer containing 20 mM Bis-tris (pH7.0), 500 mM ε-aminocaproic acid, 20 mM NaCl, 10% Glycerol, 5 mM MgCl2, and 0.3 mM GTP (Beertsen et al., 1982).

Isoelectric Focusing

Fifty microgram of total cellular protein was precipitated by acetone and dissolved in rehydration buffer containing 8 M Urea, 2% CHAPS, 50 mM dithiothreitol (DTT) and 0.2% Bio-Lyte ampholytes pH3-10. The dissolved proteins were incubated in Zoom IPG-strip 3-10 non-linear (NL) (Thermo Fisher) for 1 h and then focused at 175 volt (V) for 15 min, 175–2000 V ramp for 45 min and 2000 V for 30 min. After the run, the strips were alkylated and resolved on 2D SDS-PAGE and immunoblotted using antibodies previously described.

Statistical Analysis

Statistical analysis was performed using Prism version 7.00 (GraphPad Software). The mean of two or more groups were compared using one-way ANOVA followed by multiple comparison tests (Siegel, 1956; Dunn, 1964). The mean of multiple experimental groups were compared with the control group by Dunnett’s post hoc multiple comparison test, whereas, the mean difference between two experimental groups were compared by Sidak’s post hoc multiple comparison test (Dunn, 1964). Comparisons between two groups were performed using Student’s t-test (unpaired). Differences were considered significant at P < 0.05. Neurite outgrowth, neurite width and mitochondrial quantification data were plotted as box and whisker plot where the end of the box represents the upper and lower quartiles and the median is marked by a horizontal line inside the box. The whiskers represent the highest and lowest values excluding the outliers. In some figures the individual data points of outliers are also indicated.

Results

MT7 and Pirenzepine Significantly Augmented Neurite Outgrowth in Cultured Primary Sensory Neurons

Pirenzepine is a selective M1R antagonist whereas MT7 is the only specific antagonist of this receptor (Birdsall et al., 1983; Max et al., 1993a). Our previous study demonstrated that both antagonists enhanced neurite outgrowth from adult sensory neurons (Calcutt et al., 2017). In order to confirm that MT7 and pirenzepine had growth promoting effects under the current conditions, we measured total neurite outgrowth from primary sensory neurons derived from adult rats and cultured in defined media containing a cocktail comprising of low or high concentrations of GDNF, NGF, and NT3 growth factors (LGF and HGF, respectively) and either 100 nM MT7 or 1 μM pirenzepine (Figures 1A–D). The growth factor concentrations in the LGF cocktail reflects concentrations of growth factors that are sub-saturating and induce small but significant increases in neurite outgrowth (Calcutt et al., 2017). The HGF cocktail contained 10-fold higher concentrations to allow us to determine whether the neuritogenic effects of MT7 and pirenzepine were also effective when growth factors were present in excess and, also to see if increased tyrosine kinase signaling masked the effect of antagonist-M1R signaling mediated growth. In the absence of muscarinic antagonists, the HGF promoted significantly (p < 0.0001) more neurite outgrowth than LGF, as assessed using unbiased automated high content imaging combined with data analysis (Figures 1A,B). Both MT7 and pirenzepine significantly increased neurite outgrowth in LGF and HGF conditions within 48 h of treatment (Figures 1A,B). In addition, binning of the entire data set of 1249 (LGF) and 1517 (HGF) neurons for their total neurite length revealed that the neuritogenic effect of pirenzepine and MT7 evenly affected neurite outgrowth of the majority of the population of neurons (Figures 1C,D). Further, pirenzepine treatment elicited significantly higher neurite outgrowth under the LGF conditions when compared to MT7, whereas under the HGF condition, the primacy was reversed (Figures 1A,B). The exact reason for the difference between the drugs in terms of synergistic effect on growth at HGF condition is not known. The difference in the chemical nature of these drugs may be responsible. MT7 is a cell impermeable 7.4 kDa protein (Krajewski et al., 2001) which binds allosterically to M1R (Max et al., 1993b; Karlsson et al., 2000) whereas pirenzepine is a cell-permeable orthosteric antagonist molecule (Caulfield and Birdsall, 1998). We used LGF conditions in subsequent experiments.

FIGURE 1

FIGURE 1. M1R antagonists, MT7 and pirenzepine, augment neurite outgrowth in primary sensory neurons and M1R overexpression inhibits neurite outgrowth. (A,B) Whiskers box (Tukey) showing total neurite outgrowth per neuron. Neurons were grown for 48 h in defined media containing LGF (A) or HGF (B) condition and treated with 100 nM MT7 or 1 [mu]M pirenzepine (PZ), respectively. N = 1249 (LGF) and 1517 (HGF), respectively. P-values were calculated by one-way ANOVA followed by post hoc multiple comparison tests. Dunnett’s multiple comparisons test was used to compare the MT7 and PZ treatment groups with the control group and Sidak’s multiple comparisons test was used to compare between the MT7 and PZ treatment groups; ∗ indicates the p-value obtained by Sidak’s multiple comparisons test. (C,D) Binning of the entire data set presented in (A,B). (E) Immunoblot showing GFP-tagged muscarinic receptors (M1R to M5R) and GFP expression in transfected adult rat DRG neurons. pEGFP-C1-(M1R-M5R) plasmids were transfected in to DRG neurons and the lysate was resolved in SDS page and subsequently immunoblotted with anti-M1R (bottom panel) and anti-GFP (top panel) antibodies. (F) Time lapse confocal images showing increasing internalization (white arrows) of the GFP-M1R following treatment with carbachol (10 [mu]M). Scale bar: 10 [mu]m. (G) Immunofluorescence images showing colocalization of 24h CCh treated GFP-M1R with endosomal marker Rab5. Scale bar: 10 [mu]m. (H) Whiskers box (Tukey) showing total neurite outgrowth per neuron, N = 634 (GFP), and N = 553 (GFP-M1R), neurons, respectively. P-value was calculated by t-test (unpaired). (I) Immunofluorescence images showing [beta]-tubulin III staining and corresponding neurite trace (red lines) images in GFP and GFP-M1R overexpressed neurons. The total neurite outgrowth measurement was performed in Cellomics ArrayScan HCS Reader using neuronal profiling software. Scale bar: 10 [mu]m.

M1R Overexpression Induced Significant Reduction in Neurite Outgrowth

The muscarinic receptor subtypes (M1–M5) show considerable heterogeneity of expression in sensory neurons(Bernardini et al., 1999; Chiu et al., 2014). In order to understand the biological function of M1R in sensory neurons, we overexpressed GFP-tagged M1R in sensory neurons and measured the impact on neurite outgrowth. M1R overexpression significantly reduced neurite outgrowth when compared to GFP-expressing neurons (Figures 1H,I). The GFP-M1R transgene-induced protein production was verified in the transfected DRG neurons by immunoblotting of the expressed recombinant proteins (GFP-M1-5R) using both anti-M1R and anti-GFP antibodies (Figure 1E). The biological functionality of the GFP-M1R recombinant protein was verified by treating GFP-M1R overexpressing DRG neurons with the broad spectrum muscarinic agonist carbachol (10 μM) followed by live confocal imaging to monitor internalization of recombinant protein. Treatment with carbachol increased presence of the recombinant protein as internalized spots that were positive for the early endosomal marker Rab5 (Zerial and McBride, 2001) (Figures 1F,G). This indicates that GFP-tagged recombinant M1R elicited agonist induced internalization response similar to that of naïve proteins. In addition, to assess functionality of the recombinant M1R in overexpressed cells, we performed a β-arrestin recruitment assay [TANGO assay: (Kroeze et al., 2015), Figure 2]. In the TANGO assay, upon activation, β-arrestin is recruited to the C-terminus of the M1R-TEV-tTA fusion protein at the TEV protease site and cleaves to release the tTA transcription factor, which after transport to the nucleus activates transcription of luciferase reporter gene (Kroeze et al., 2015) (Figure 2A). We verified expression of the M1R-TEV-tTA and β-arrestin transgenes by immunoblotting (Figure 2B). We found significant recruitment of β-arrestin due to basal M1R activation in presence of serum in the culture media in the HTLA cells (Figure 2C). This basal activation of M1R occurred in the absence of exogenously supplied agonist. However, fetal bovine serum (FBS) used in the culture media is reported to contain acetylcholine (ACh) (Lau et al., 2013) and this may have constitutively acted upon the M1R. Treatment with the muscarinic agonist, carbachol, significantly increased activation of the M1R by up to ∼7.4-fold (Figure 2C).

FIGURE 2

FIGURE 2. Constitutive basal activity of M1R and the effect of DREADD mutant on neurite outgrowth. (A) Graphical representation of the arrestin recruitment assay (TANGO) strategy, TRE-Tet response element. HA-cleavable signal sequence to promote membrane localization, FLAG-epitope tag, TEV- Tobacco Etch Virus cleavage site, V2 tail: C-terminus of the V2 vasopressin receptor (V2 tail) to promote arrestin recruitment and tTA- tetracyclin transactivator. (B) Immunoblot showing expression of M1R-TEV-tTA and [beta]-arrestin-TEV protease transgene. (C) Scatter plot showing the RLU (Relative Luminescence Units) for the drug treated HTLA cells. N = 5 independent experiments. (D) Co-expression of GFP/GFP-M1R/GFP-M1R(DREADD) with mito7-RFP in DRG neurons following 48 h of growth. Neurite outgrowth suppression was comparatively less in M1R-DREADD mutant over-expressing neurons. The white and blue arrows indicate that localization of GFP-M1R was mainly restricted in the perikaryon where as GFP was localized both in perikaryon and neurites. LSM510 confocal images acquisition parameters were same for all images. Scale bar: 100 [mu]m. (E) Whiskers box (Tukey) showing total neurite outgrowth per neuron. N = 250. p-value by one-way ANOVA followed by Dunnett’s multiple comparisons test. (F) Binning of the entire data set presented in (E).

Endogenous ACh Binds M1R to Constrain Neurite Outgrowth

Adult sensory neurons secrete ACh into the extracellular media when grown in culture (Calcutt et al., 2017). The Tango Assay revealed that ACh present in the serum supplied to the culture media may also act upon M1R and constitutively activate the receptor. We therefore hypothesized that overexpression of M1R in sensory neurons may lead to increased constitutive activation of the receptor by recruitment of trimeric G-proteins, including the Gα12/13 subtype that promotes tubulin destabilization (Roychowdhury et al., 1999). In order to test this hypothesis we transfected DRG neurons with an M1R-DREADD (designer receptors exclusively activated by designer drugs) mutant to eliminate the putative basal activity in M1R-DREADD overexpressing neurons (Figures 2D,E). The M1R-DREADD mutant contained two mutations in the conserved orthosteric site residues (Y106C and A196G in the M1R) that minimize responsiveness to ACh (Abdul-Ridha et al., 2013). Overexpression of GFP-M1R and M1R-DREADD caused low levels of neurite outgrowth (Figures 2D–F). However, neurite outgrowth in M1R-DREADD neurons was significantly higher than GFP-M1R alone over 72 h of culture (Figures 2D–F). This indicates that endogenous ACh may interact with M1R to limit neurite outgrowth.

Impact of M1R Overexpression on Mitochondrial Abundance

In an attempt to understand the mechanism of the growth inhibitory effect of M1R overexpression, we examined the abundance of mitochondria in the neurites by co-expressing DsRed2Mito7 plasmid (Figure 3A). The DsRed2Mito7 plasmid expresses DsRed protein tagged with mitochondrial targeting sequence from subunit VIII of human cytochrome C oxidase and therefore localizes specifically to mitochondria. We measured the volume of mitochondria present per unit length (μm) of neurites in the M1R and GFP overexpressing DRG neurons (Figure 3B). M1R overexpression significantly reduced the abundance of mitochondria in the neurites after 48 h of culture (Figures 3C,D). The mitochondria in both M1R and GFP over-expressing neurons accumulated mitotracker CMXRos dye, indicating that functionality of the mitochondria was not impaired in M1R overexpressed neurons despite reduced abundance (data not shown).

FIGURE 3

FIGURE 3. GFP-M1R over-expression in sensory neurons caused reduced mitochondrial abundance. (A) Representative images showing the signal of co-expression of GFP/GFP-M1R and mito7-RFP in DRG neurons. BF, bright field. (B) Immunofluorescence images showing mito7-RFP and [beta]-tubulin III stained neurites. The 4th panel image showing tracing of the mitochondrial volume used for determining the amount of mitochondria per unit length of the neurites in GFP/GFP-M1R overexpressed DRG neurons. The tracing image was created using ImageJ particle counter plugin. Scale bar: 5 [mu]m. (C) Whiskers box (min–max) showing amount of mitochondria in GFP/GFP-M1R transfected neurons. N = 101 from three independent experiments, p-value by t-test (unpaired). (D) Binning of the entire data set presented in (C).

Reduced Mitochondrial Abundance Was Associated With Impaired Cytoskeletal Structure

To determine whether reduced mitochondrial abundance was related to a defect in the actin or tubulin cytoskeletons, we immunostained DRG neurons that overexpressed GFP/GFP-M1R with phalloidin and anti-β-tubulin III antibodies. In the presence of M1R over-expression, immunofluorescent imaging of phalloidin revealed the presence of abundant and continuous actin filaments in the neurites, although the neurite tips were notably thinner compared to those of GFP expressing neurites (Figure 4A). In contrast, the β-tubulin III associated cytoskeleton appeared less abundant and discontinuous (Figure 4B). Fragmentation of the tubulin cytoskeleton was confirmed using a polymerized tubulin quantification assay in which there was significantly less polymerized tubulin in the DRG neurons overexpressing M1R compared to those that overexpressed only GFP (Figures 4C,D).

FIGURE 4

FIGURE 4. Dorsal root ganglia (DRG) neurons over-expressing M1R exhibited cytoskeletal defects. (A) Actin cytoskeleton in GFP-M1R and GFP over-expressing DRG neurons. White arrow indicates narrower neurites in GFP-M1R expressing neurons. (B) [beta]-tubulin III associated cytoskeleton in GFP-M1R and GFP over-expressing neurons. Black rectangular area is shown in a magnified view in the right panel; white arrow indicates discontinuous/continuous tubulin cytoskeleton. BF, bright field. Scale bar: 20 [mu]m. (C) Immunoblot and (D) bar graph showing relative amount of polymerized tubulin in the M1R-GFP and GFP over-expressed neurons. The data represent the mean ± SEM of three independent experiments. P = 0.0011 calculated by unpaired t-test.

Discontinuous Tubulin Cytoskeleton Was Associated With Reduced Mitochondrial Trafficking in Sensory Neurons Overexpressing M1R

In order to determine whether the discontinuous β-tubulin cytoskeleton was associated with altered mitochondrial trafficking in growing neurites, we measured velocity of mitochondrial movement in the neurites of DRG neurons that co-expressed GFP/GFP-M1R and DsRed2Mito7. Time-lapse images of the neurites at 10 s intervals were used to generate kymographs of mitochondrial trafficking (Figure 5A). The mean velocity of mitochondria was calculated as 1.9 μm/sec in GFP expressing neurons whereas in GFP-M1R expressing neurons it was significantly lower at 0.62 μm/sec (Figures 5B,C). In addition, immunostaining of neurites with a similar appearance to those used in the mitochondrial velocity measurement revealed a continuous β-tubulin-III cytoskeleton in GFP over-expressed neurons and a discontinuous β-tubulin-III cytoskeleton in M1R over-expressed neurons (Figures 5D,E).

FIGURE 5

FIGURE 5. Altered mitochondrial trafficking in M1R overexpressing neurons was rescued by antagonist treatment. (A) Representative images (top panels) of the first frame of a series of live cell time-lapse images showing the expression of DsRed2Mito7 in the mitochondria of GFP/M1R–GFP expressing neurons. The circles represent mitochondria identified by MTrackJ plugin, which then tracked migration through time in a series of time lapse images and calculated velocity of specific mitochondria. White/blue line represents neurite trace. The bottom panel represents Kymographs generated from live cell time-lapse images. The Kymograph was generated using ImageJ Kymograph plugin. The X-axis represents the physical location of mitochondria on the neurite, and the Y-axis represents the location of mitochondria in time. Streak of particles traversing the kymograph from left to right in angular lines indicates retrograde/anterograde mitochondrial motion. (B) Whisker plot showing mitochondrial velocity. DRG neurons were cultured in LGF media supplemented with 100 nM MT7 or 1 [mu]M pirenzepine for 48 h following transfection. N = 40 from three independent experiments, p-value by one-way ANOVA followed by Dunnett’s multiple comparisons test. (C) Binning of the entire data set presented in (B). (D,E) Immunofluorescence images showing mito7-RFP and [beta]-tubulin-III staining in the GFP/GFP-M1R expressing neurites. White arrows indicate continuous/discontinuous tubulin cytoskeleton in GFP/GFP-M1R expressing neurites, respectively. Scale bar: 5 [mu]m.

Muscarinic Antagonists MT7 and Pirenzepine Rescued the Cytoskeletal Defect, Aberrant Mitochondrial Distribution and Trafficking in Neurites

We investigated whether pirenzepine or MT7 could overcome M1R overexpression-induced cytoskeletal defects. M1R over-expressing neurons were maintained for 48 h and then treated for 24 h with 100 nM MT7 or 1 μM pirenzepine and total neurite outgrowth and mitochondrial velocity (Figure 5) and abundance (Figure 6) quantified. MT7 and pirenzepine treatment significantly rescued the deficits in the mitochondrial velocity in M1R overexpressed neurons, with mean mitochondrial velocities of 0.78/1.2 μm/s in MT7/pirenzepine treated neurons being significantly higher than untreated neurons (Figures 5B,C). Both MT7 and pirenzepine also caused considerable re-localization of M1R from the perikarya to the neurites, as revealed by time lapse confocal live cell images over a period of 72 h. This may indicate increased vesicular transport of internalized M1R (Figure 6A). Within the same neuron (shown in Figure 6A), following 24 h of drug treatment and upon fixation and immunostaining for β-tubulin III, there was continuity in the β-tubulin III associated cytoskeleton (right panel). MT7 or pirenzepine significantly increased total neurite outgrowth, reversed the reduced neurite caliber and increased mitochondrial abundance in DRG neurons that overexpressed GFP-M1R (Figures 6B–E).

FIGURE 6

FIGURE 6. Restoration of cytoskeleton, mitochondrial abundance and neurite outgrowth by M1R antagonists MT7 and pirenzepine treatment. (A) Time-lapse live confocal images of GFP-M1R over-expression in sensory neurons showing MT7 induced re-localization of M1R from perikaryon to neurites. The overexpressed neurons were grown in defined media for 48 h and imaged (left panel). Neurons were then treated with 100 nM MT7 for 24 h and imaged (middle panel). Right panel: The same neuron depicted left was fixed and stained for [beta]-tubulin III to show continuity of cytoskeleton. Scale bars: 50 [mu]m. (B,C) Whiskers box (Tukey) showing total neurite outgrowth per neuron (B) and average neurite width per neuron (C). p-value calculated by one-way ANOVA, N = 224, 230 and 198 cells, respectively, for (A) and N = 1250, 1268, 1280, and 1306 cells, respectively, for (B). Asterisks indicate p-value calculated by unpaired t-test. (D) Whisker box plot showing amount of mitochondria in the M1R expressing neurons treated with 100 nM MT7 or 1 [mu]M pirenzepine. N = 101, p-values were calculated by one-way ANOVA followed by Dunnett’s multiple comparisons tests. (E) Binning of the entire data set presented in (D).

Knockdown of Gα13 in Sensory Neurons Reversed the M1R Overexpression-Induced Tubulin Cytoskeleton Defect

The M1R-DREADD mutant study raised the possibility that, in normal DRG neurons, basal M1R activity resulting from binding of endogenous ACh release may destabilize the tubulin cytoskeleton through increased active G proteins. Gα12/13 proteins are known for their effect on cytoskeleton remodeling and the M1R receptor activates Gα12/13 type G proteins leading to mobilization of the small GTP-binding protein Rho through activation of Rho-GEF (RhoGTPase nucleotide exchange factor) (Haga, 2013). We, therefore, measured the relative expression of Gα12 and Gα13 proteins in sensory neurons (Figure 7A). Immunoblots revealed that DRG neurons express significantly more Gα13 compared to Gα12 (Figures 7A,B). The relative levels of Gα12 and Gα13 expression were comparable to those of the human carcinoma cell line HEK293 (Supplementary Figure S1A). We used siRNA to knockdown Gα13 in DRG neurons and used a CRISPR/Cas9 based Gα12/13 null HEK293 cell line to study the effect of cholinergic signaling through Gα13 on the tubulin cytoskeleton (Figure 7C and Supplementary Figure S1A). The HEK293 cell line has been reported to exhibit a neuronal lineage phenotype and express neuronal proteins (Stepanenko and Dmitrenko, 2015). Therefore, we considered it suitable for this study. The siRNA based knockdown of Gα13 in DRG neurons had no significant effect on neurite outgrowth (Figures 7D,E). Interestingly, knockdown of Gα13 in M1R overexpressed sensory neurons significantly reversed the suppressed neurite outgrowth (Figures 7F,G). In addition, treatment of the M1R overexpressed and Gα13 knockdown neurons with 100 nM MT7 or 1 μM pirenzepine exhibited significantly increased neurite outgrowth (Figures 7F,G).

FIGURE 7

FIGURE 7. Knockdown of G[alpha]13 reversed M1R overexpression-induced inhibition of neurite outgrowth. (A) Immunoblots showing relative expression of G[alpha]12 and G[alpha]13 proteins in cultured DRG neurons. (B) Scatter plot showing relative amount of G[alpha]12 and G[alpha]13 proteins in cultured sensory neurons. N = 5 independent experiments. p-value was calculated by unpaired t-test. (C) Immunoblots showing siRNA (cocktail of 3 siRNAs targeted to rat G[alpha]) based knockdown of G[alpha]13 protein in cultured adult rat DRG neurons. (D,F) Whisker box (Tukey) showing total neurite outgrowth per neuron. DRG neurons were transfected with a pEFGP-C1-M1R plasmid and siRNA using Amaxa nucleofection reagent and allowed to grow for 48 h. Scrambled siRNAs were used for control. In drug treatment groups, neurons were cultured in media supplemented with 100 nM MT7 or 1 [mu]M pirenzepine following transfection. The neurons were fixed after 48 h of culture, stained with [beta]-tubulin III and imaged using Cellomics ArrayScan HCS Reader. p-value by unpaired t-test or one-way ANOVA test followed by Dunnett’s multiple comparisons tests. N = 432/452 (in D) and 406/694 (in F). (E,G) Binning of the entire data set presented in (D,F).

CRISPR/Cas9 Based Gα12/13 Null HEK293 Cells Showed Abundant Tubulin Cytoskeleton but Diminished Actin Stress Fibers

RNA sequencing data in the human protein atlas (HPA) indicates that HEK293 cells express very high levels of α-tubulin compared to β-tubulin III (Uhlen et al., 2005). The transcripts per million bases (TPM) value for α-tubulin (TUBA1B isoform) and β-tubulin III has been recorded as 1923.6 and 13.7, respectively1. β-tubulin III was undetectable by immunoblotting in HEK293 cells (data not shown). We, therefore, studied α-tubulin dynamics in CRSIPR/Cas9 based Gα12/13 knockout HEK293 cells. We performed immunofluorescent labeling using phalloidin and anti-α-tubulin specific antibodies to visualize the actin and tubulin based cytoskeletal structures in Gα12/13 knockout HEK293 cells (Supplementary Figure S1E). Phalloidin staining revealed that actin stress fibers were diminished in Gα12/13 null cells, with appearance of abundant distinct punctate actin rich focal adhesion points (Supplementary Figure S1E). Further, the tubulin cytoskeleton appeared more robust and organized in Gα12/13 null cells (Supplementary Figures S1E, S2). The BN-PAGE based microtubule fractionation assay and polymerized microtubule quantitative assay showed Gα12/13 null cells exhibited significantly more polymerized tubulin than wild type HEK293 cells (Supplementary Figures S1B–D). Overexpression of M1R in Gα12/13 null cells and treatment with muscarinic agonist carbachol did not alter the tubulin networks as compared to wild type cells (Supplementary Figures S2, S3).

MT7 and Pirenzepine Modulate G Protein Interaction With the M1R

Adult DRG cultures have very low cell yields, in the range of 250,000 per culture. To enable feasible pull down of protein complexes we used human neuroblastoma SH-SY5Y cultures that allow use of millions of cells. This cell line also exhibits a cholinergic phenotype with ACh secretion, expression of muscarinic receptors and AChE shedding (Yamada et al., 2011). We examined the recruitment of trimeric G proteins to M1R by Halo-pull down assay and BN-PAGE analysis (Figures 8, 9). The Halo-pull down assay permitted a focus on the over-expressed M1R with no contamination from endogenous muscarinic receptors of mixed sub-type. The SH-SY5Y cells that overexpressed Halo-M1R were treated with drugs and Halo-tagged M1R was pulled down using Halo-linked resin. Subsequently, the pull down product was resolved in SDS-PAGE and immunoblotted with anti-M1R, anti-Gα12/13 and anti-Gγ/2/3/4/7 antibodies (Figures 8B–D). The Halo-M1R pull down fraction from the drug treated cells showed significantly elevated levels of Gγ/2/3/4/7 and Gα12/13 proteins when compared with untreated cells (Figures 8E,F).

FIGURE 8

FIGURE 8. MT7 and pirenzepine elevated sequestration of trimeric G proteins associated with M1R in SH-SY5Y cells. (A) Diagrammatic representation of the experimental strategy. (B–D) Halo-tagged M1R was expressed transiently in sensory neurons and then treated with 100 nM MT7 or 1 [mu]M PZ for 1 h. Cells were then lysed and Halo-M1R was pulled down using halo-linked resin. The halo tag was cleaved by TEV protease and the cleaved M1R associated multiprotein complex (MPC) was resolved in denaturing SDS-PAGE and immunoblotted using (B) anti-M1R, (C) anti-Gγ2/3/4/7, and (D) anti-G[alpha]12/13 antibodies. (E,F) Scatter plot showing the relative amount of G proteins (G[alpha] and Gγ, respectively) associated with M1R following drug treatment. The data represent mean ± SEM of three independent experiments. p-values (∗∗∗∗ < 0.0001) were calculated by one-way ANOVA with Dunnett’s post hoc multiple comparison test.

FIGURE 9

FIGURE 9. MT7 and pirenzepine elevated sequestration of trimeric G proteins in a protein complex associated with M1R in SH-SY5Y cells. (A,B) BN-PAGE analysis showing recruitment and sequestration of G proteins to M1R following MT7 treatment. GFP-M1R transfected cells were treated with 100 nM MT7 and incubated for 1 h. Cell lysates were then separated on 1D BN PAGE followed by 2D SDS-PAGE, (A) Control, (B) 100 nM MT7 treatment. Top panel: Coomassie stained gel piece showing 1D BN-PAGE separation of native page protein molecular weight marker. The red horizontal and vertical arrows indicate the direction of the 1D BN-PAGE and 2D SDS-PAGE, respectively. The blue rectangle shows PTMs of GFP-M1R associated with 1000 kDa and >1200 kDa MPCs. The red circle indicates native form of the GFP-M1R associated with the MPCs. The red arrow connecting the red circles indicates shift of molecular weight in MPC due to recruitment of G proteins following drug treatment. The red rectangle in the bottom panel shows the co-migration of possible interacting G proteins with the MPCs.

The 2D BN-PAGE/SDS-PAGE analysis revealed existence of 2 protein complexes at ∼900 and ∼1200 kDa equivalent molecular weights in Halo-M1R over-expressing SH-SY5Y cells (Figures 9A,B). Each complex was associated with a native ∼100 kDa Halo-tagged M1R protein and >180 kDa fractions, the latter may be derived from different PTMs of the Halo-M1R. Treatment with 100 nM MT7 caused a major shift of the ∼900 kDa protein complex to the ∼1200 kDa protein complex within 1 h of treatment, suggestive of recruitment of putative interacting proteins (Figure 9, blue and red dotted areas). In contrast, untreated cells did not show any shift in the ∼900 kDa protein complex to the ∼1200 kDa protein complex indicating less or no recruitment of interacting proteins (Figure 9A). Further, when the same blots were immunoblotted with anti-Gα12/13 antibodies, the Gα12/13 proteins appeared as spots on a vertical line corresponding to the ∼1200 kDa protein complex in the drug treated cells which suggests possible co-migration and association with M1R (Figures 9A,B, bottom panel).

Discussion

Our recent work has shown that sensory neurons derived from M1R null mice exhibit enhanced neurite outgrowth (Calcutt et al., 2017). We now demonstrate that over-expression of M1R inhibited neurite outgrowth, caused disruption of the tubulin cytoskeleton and blockade of mitochondrial trafficking in adult sensory neurons, all of which were rescued by exposure to selective or specific M1R antagonists. We then used overexpression of GFP-M1R to identify the molecular pathway components associated with specific M1R-mediated cellular phenotypes. Based on our data, we propose that ACh mediated signaling via M1R constrains neurite outgrowth via activation of Gα13 proteins, which in turn limits tubulin polymerization and mitochondrial trafficking within axons.

Overexpression of M1R in sensory neurons has biological consequences that likely arise from the presence of neuron-derived ACh in the culture environment. Cultured sensory neurons secrete endogenous ACh into the extracellular media to generate a local concentration in the range of approximately 16 μM (Calcutt et al., 2017). This far exceeds the ACh Kd of 25–35 nM measured in several regions of rat brain tissue (Kellar et al., 1985) or the Kd of 0.2–0.4 nM measured using rat brain neurons (Pavia et al., 1991; Bakker et al., 2015). In the presence of abundant extracellular ACh, overexpressed. M1R will trigger activation of Gα13 proteins that leads to the dissociation of tubulin microtubules, as seen in Figure 4. The chemogenetically modified M1R-DREADD mutant significantly reduced M1R induced growth retardation (Figure 2). Thus, ACh-driven basal activity of M1R is responsible for neurite outgrowth suppression under conditions of M1R over-expression. In addition, we propose in normal un-transfected neurons that basal endogenous M1R signaling tonically suppresses neurite outgrowth by restricting mitochondrial (and potentially vesicular) transport, thereby explaining the ability of antimuscarinic drugs to prevent/reverse this constraint.

Knockdown of Gα13 in M1R-overexpressed neurons significantly reversed the M1R-induced growth inhibitory effect and protected from tubulin destabilization in growth cones (Figure 7). Gα13 is more abundant than Gα12 in sensory neurons and M1R was linked to activation of Gα12/13 type G proteins that leads to activation of small GTP binding protein Rho through mobilization of RhoGEF (Luo et al., 2001; Siehler, 2009; Haga, 2013). Activation of Gα12/13 leads to stimulation of the Rho/Rho Kinase pathway via a subgroup of Rho guanine nucleotide exchange factors (Fukuhara et al., 2001). Involvement of Rho and ROCK (Rho-associated coiled coil forming protein kinase) in agonist-induced neurite retraction and cell rounding has been reported in N1E-115 neuroblastoma cells (Hirose et al., 1998). Dominant-negative p160-ROCK completely abolished this neurite retraction suggesting a clear link between RhoA-ROCK signaling and cytoskeleton disassembly (Hirose et al., 1998). Activation of neuronal cannabinoid receptors linked to Gα12/13 proteins triggered rapid and reversible contraction of actinomyosin cytoskeleton through a Rho-GTPase and ROCK (Roland et al., 2014). Further, the CRISPR/Cas9 based knockout of Gα12/13 proteins in the HEK293 cell line augmented the abundance of polymerized tubulin (Supplementary Figure S1). HEK293 cells have been utilized as a model system for neuronal synapse formation (Biederer and Scheiffele, 2007) as they express neuronal proteins and have neuronal cell-lineage (Shaw et al., 2002; Stepanenko and Dmitrenko, 2015). Using a NanoBIT split luciferase based RhoA biosensor that detects Gq-induced RhoA activation showed that the RhoA signal is completely lost in the Gα12/13 KO cell (Mercier et al., manuscript in revision, personal communication).

Activated G proteins regulate tubulin polymerization (Schappi et al., 2014) and tubulin binds directly to Gα or Gβγ subunits (Wang et al., 1990; Roychowdhury et al., 1999; Sarma et al., 2003). Activated GTP bound Gα promotes microtubule instability by increasing the intrinsic hydrolysis of GTP-tubulin to the less polymer stable GDP-tubulin (Roychowdhury et al., 1999, 2006). We have modelled this in Figure 10. Overexpression of Gαq in a rat pituitary cell line showed a 50% decrease in the ratio of soluble to polymerized tubulin (Ravindra et al., 1996). There is considerable cell type and isoform specificity in Gα mediated tubulin cytoskeleton dynamics (Sarma et al., 2015). We found that Gα12/13 and Gγ(2,3,4,7) were sequestered upon M1R antagonist binding (Figures 8, 9) and it is plausible that these factors also regulate tubulin polymerization in sensory neurons. In addition, some guanine nucleotide exchange factors (GEFs) for Rho GTPases, namely p115 RhoGEF (Kozasa et al., 1998), PDZ-RhoGEFs (Fukuhara et al., 1999), and LARG (Suzuki et al., 2003) can act as direct couplers of Gα12/13 proteins to small GTPases such as RhoA, Rac1, and CDC42, all of which are known to influence microtubule dynamics (Hall and Lalli, 2010). The Gα12/13-RhoGEF-RhoA pathway of GPCR has been implicated in many diseases (Siehler, 2009). Interestingly, while Gα promotes tubulin disassembly by increasing the tubulin GTPase activity, Gβγ subunits preferentially associate with GDP-bound tubulin to promote polymerization and stability of the microtubule (Roychowdhury and Rasenick, 1997; Popova and Rasenick, 2003; Roychowdhury et al., 2006). Our data from non-crosslinked halo-pull down (Figure 8) and BN-PAGE analyses (Figure 9) show augmented association and occupancy of trimeric G-proteins (Gα12/13) to M1R during MT7 or PZ treatment, which indicates elevated sequestration of these proteins on M1R. We propose that muscarinic antagonist-induced sequestration of trimeric G-proteins restricts their dissociation from the overexpressed M1R associated protein complex and thereby limits their detrimental effect on tubulin polymerization (see Figure 10). Further, we posit that the same M1R suppression pathway is occurring in normal un-transfected neurons to limit cytoskeleton formation in axons and is counteracted by these drugs. However, this does not exclude the possibility that other pathways may be involved. For example, Ca²⁺ signaling/homeostasis are known to be responsible for the maintenance of cytoskeletal integrity (Tsai et al., 2015). Cholinergic activation of M1R coupled with the Gq/11 protein generates cytosolic calcium transients via phospholipase-C signaling pathway (Langmead et al., 2008). High Ca²⁺ may act to increase the intrinsic GTP hydrolysis of tubulin and directly destabilize growing microtubule ends without changing the effective concentration of tubulin (O’Brien et al., 1997). Therefore, it is possible that excessive cholinergic signaling may imbalance intracellular Ca²⁺ homeostasis and promotes tubulin destabilization. Antagonist mediated sequestration of G-proteins may limit this response and protect Ca²⁺ induced tubulin destabilization. However, further experimentation is required to prove this hypothesis. Insertion and site-directed mutagenesis based studies have revealed potential G-protein interaction sites in the i2 and i3 loops in mAChRs (Blin et al., 1995; Liu et al., 1996; Hu et al., 2010). Overexpression of mutated M1R for disruptive Gα13 binding and subsequent reversal of cholinergic tubulin destabilization effect would provide another means of testing. The exact binding site for Gα13 protein in M1R needs to be evaluated.

FIGURE 10

FIGURE 10. Model explaining the effect of M1R overexpression and antagonism on tubulin associated cytoskeleton and mitochondrial trafficking. (Left) M1R overexpressing neurons have increased basal activity in response to secreted acetylcholine signaling through overexpressed receptor. This is turn causes recruitment of trimeric G proteins that destabilize tubulin polymers by increasing intrinsic GTPase activity of tubulin. The lack of tubulin cytoskeleton in M1R overexpressing neuron leads to decreased mitochondrial trafficking and stagnation of mitochondria in the neurites that impairs outgrowth. (Right) M1R overexpression in the presence of antagonists MT7 and pirenzepine (PZ) stabilizes tubulin polymerization. The antagonists bind to the M1R and may stabilize a specific structural ensemble, which in turn recruit trimeric G proteins. However, the antagonist mediated M1R structural ensemble may sequester the bound G proteins and makes them unavailable for exerting their effect on tubulin polymerization which in turn stabilizes microtubule cytoskeleton and promotes mitochondrial trafficking and neurite outgrowth.

Overexpression of M1R restricted the number of active mitochondria in neurites, presumably through diminished trafficking as a direct consequence of overexpressed M1R-induced disruption of the cytoskeleton. The M1R antagonists, pirenzepine and MT7, were able to protect mitochondrial transport. Pirenzepine and MT7 also enhance mitochondrial oxygen consumption rate and respiratory complex activities through activation of the AMP-activated protein kinase (AMPK)/peroxisome proliferator-activated receptor γ coactivator-1α signaling axis (Calcutt et al., 2017). Mitochondrial function was also enhanced in sensory neuron cultures derived from M1R null mice. Thus, ACh signaling through M1R can negatively regulate mitochondrial phenotype at multiple levels in the neuron that include trafficking and positioning as well as fine regulation of activity of the respiratory complexes. Optimal regulation of mitochondrial function is critical in distal regions of sensory neurons which, in human peripheral nerve, can be up to a meter from the cell body (Schwarz, 2013). The architecture of sensory neurons poses an extreme cellular environment for mitochondrial distribution and a need to supply energy to the distal endings where energy demand is high (Bernstein and Bamburg, 2003; Chowdhury et al., 2013). Any defect in mitochondrial function is likely to have a profound influence on the axon. Recent in vivo studies of mitochondrial transport along the saphenous nerve of adult mice revealed elevated anterograde transport of this organelle in axons undergoing high rates of depolarization and impulse conduction (Sajic et al., 2013). Indeed, loss of function of mitochondrial proteins such as bcl-w or mitofusin-2 results in a length-dependent dying-back sensory neuropathy (Baloh et al., 2007; Misko et al., 2010; Courchesne et al., 2011). Genetic ablation of mitochondrial transport also leads to axonal growth failure following axotomy in mice (Zhou et al., 2016). In axons, approximately 30–40% of total mitochondria are constantly engaged in saltatory motion (Lovas and Wang, 2013). It is plausible that the microtubule disruption induced by overexpression of M1R eliminated the basic framework for motor proteins to carry their mitochondrial cargo and potentially other vesicular cargos. As a result, the physical abundance of active mitochondria was diminished in the distal neuritis, which would be expected to deprive the neuronal growth cone of an essential supply of ATP, resulting in suppressed actin treadmilling and neurite outgrowth. This hypothesis, summarized in Figure 10, is further supported by reports that several microtubule-targeted chemotherapeutic agents, such as colchicine and vincristine, are known to induce a sensory neuropathy in which the distal aspect of the sensory axon gradually degenerates (Bennett et al., 2014).

Our pre-clinical studies have demonstrated that selective and specific M1R antagonists promote neurite outgrowth in adult sensory neurons in vitro (Calcutt et al., 2017). MT7 and pirenzepine were also able to prevent or reverse the distal degenerative neuropathy characteristic of diabetes, chemotherapy-induced and HIV-induced peripheral neuropathy. In the present study, we have reinforced the finding that M1R antagonists are neuronal growth promoters using an unbiased automated high throughput neurite outgrowth measurement technique in a large cohort of sensory neurons. In addition, we have highlighted a molecular mechanism that indicates that treatment with these drugs stabilizes microtubules by sequestering Gα13 proteins and promotes microtubule-based axonal transport of mitochondria, which in turn augments neurite outgrowth. Peripheral neuropathy is a major cause of human morbidity with huge associated health care costs (Gordois et al., 2003; McInnes, 2012). One particularly encouraging implication of our identification of the endogenous M1R-mediated suppression of axonal outgrowth in sensory neurons is that antimuscarinic drugs that can prevent/reverse this process have already been widely used as approved drugs for other indications (Siatkowski et al., 2008). Since small fiber axonal degeneration is an early feature of many peripheral neuropathies, the novel growth-regulating pathway we have identified could be mobilized to prevent or reverse distal neurodegeneration. Our studies support a previously unrecognized therapeutic potential for M1R antagonists in the treatment of peripheral neuropathies and unravels a novel pathway of cholinergic signaling mediated via control of microtubule dynamics through Gα13 signaling.

Author Contributions

MS conceptualized, designed and performed all the experiments, analyzed all data, prepared all figures, and wrote the manuscript. PF and NC obtained funding to support the work, aided in the design of experiments, and contributed to writing and editing of the manuscript.

Funding

This work was funded through Canadian Institutes of Health Research (CIHR) grant numbers MOP-130282 and RPA-124953 (PF) and National Institutes of Health award NS081082 (NC).

Conflict of Interest Statement

NC and PF declare that they are scientific founders of, have an equity interest in, WinSanTor Inc., a company that has licensed IP from the University of Manitoba and University of California San Diego and may potentially benefit from the research contained herein. The terms of this arrangement have been reviewed and approved by the University of California, San Diego in accordance with its conflict of interest policies. PF also serves on the Board of Directors of WinSanTor Inc.

The remaining author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We thank St. Boniface Research for additional funding support. We sincerely acknowledge Professor Arthur Christopoulos, Monash University, Australia, for providing M1R-DREADD mutant plasmid; Dr. Asuka Inoue, Tohoku University, Japan, for providing β-arrestin and Gα12/13 null cell lines; Prof. Bryan Roth, University of North Carolina, for providing HTLA cells; Dr. Jun-Feng Wang, University of Manitoba, for providing SH-SY5Y cell line. We pay tribute to the late Dr. Michael W. Davidson, Florida State University, for his lifelong contribution to science and acknowledge his contribution to the Michael Davidson fluorescent protein collection repository for public use. We thank Jennifer Chung, Santa Cruz Biotechnology, Inc. for providing numerous mouse monoclonal antibodies, Dr. Emma Van Der Westhuizen, Monash University, for sending cell lines and plasmids; Dr. Gordon Glazner, University of Manitoba, for access to the LSM510 confocal microscope.

Supplementary Material

The Supplementary Material for this article can be found online at: https://[see PDF]/articles/10.3389/fnins.2018.00402/full#supplementary-material

FIGURE S1 | Polymerized α-tubulin is augmented in Gα12/13 null HEK293 cells. (A) Immunoblots showing expression of Gα12, Gα13, and porin (VDAC1) in wild type and Gα12/13 null HEK293 cell. (B) Immunoblots showing the relative amount of polymerized tubulin in wild type and Gα12/13 null HEK293 cells. (C) In BN-PAGE based assay, polymerized tubulin was stabilized in microtubule stabilization buffer and resolved in native gradient gel. The polymerized tubulin microtubules are trapped in the native page at the top (∼1200 kDa) whereas the dimers and monomers migrated to lower molecular weight regions. The first dimension gel was denatured and the proteins were further separated in 2nd dimension SDS-PAGE and immunoblotted. The spots at ∼146 and ∼66 kDa represent α-tubulin dimer and monomers. The red dotted circle represents the polymerized tubulin. Green and blue arrows indicate direction of protein movement. (D) Scatter plot showing polymerized tubulin in wild type and Gα12/13 null HEK293 cell as revealed by polymerized tubulin and BN-Page based assays. In the BN-PAGE assay the polymerized tubulin was quantified based on the relative intensity of the monomeric tubulin. N = 3/4 independent experiments. Data represented as mean ± SEM, p-value by unpaired t-test. (E) Immunofluorescent images showing actin (phalloidin) and (α-tubulin) in HEK293 (Wild type: WT) and GNAS12/13 knockout (ΔGNAS12/13) cells. Scale bar: 20 μm.

FIGURE S2 | Actin and tubulin cytoskeleton in wild type and G12/13 KO cells. (A,B) Confocal immunofluorescent images showing F-actin (phalloidin stained, top panel, A) and tubulin (α-tubulin immunolabelled, bottom panel, B) cytoskeleton in GFP-M1R expressing wild type and G12/13 KO cells. Blue arrow indicates focal adhesions, yellow arrow indicates extended cytoplasmic processes enriched in tubulin. Scale bar: 20 μm.

FIGURE S3 | Effect of carbachol on tubulin cytoskeleton in M1R expressed wild type and G12/13 KO cells. Confocal immunofluorescence images showing the tubulin cytoskeleton. Yellow arrow indicates crest of the cell with dense packing of tubulin. White and blue arrows indicate elongated cytoplasmic processes enriched in tubulin. Blue arrow indicates localization of GFP-M1R in extended cytoplasmic processes. Z represents a particular optical slice in Z-stacked image series. Scale bar: 20 μm.

Footnote

1. ^www.proteinatlas.org

References

Abdul-Ridha, A., Lane, J. R., Sexton, P. M., Canals, M., and Christopoulos, A. (2013). Allosteric modulation of a chemogenetically modified G protein-coupled receptor. Mol. Pharmacol. 83, 521–530. doi: 10.1124/mol.112.083006

Akude, E., Zherebitskaya, E., Chowdhury, S. K., Smith, D. R., Dobrowsky, R. T., and Fernyhough, P. (2011). Diminished superoxide generation is associated with respiratory chain dysfunction and changes in the mitochondrial proteome of sensory neurons from diabetic rats. Diabetes 60, 288–297. doi: 10.2337/db10-0818

Bakker, G., Vingerhoets, W. A., Van Wieringen, J. P., De Bruin, K., Eersels, J., De Jong, J., et al. (2015). 123I-iododexetimide preferentially binds to the muscarinic receptor subtype M1 in vivo. J. Nucl. Med. 56, 317–322. doi: 10.2967/jnumed.114.147488

Baloh, R. H., Schmidt, R. E., Pestronk, A., and Milbrandt, J. (2007). Altered axonal mitochondrial transport in the pathogenesis of Charcot-Marie-Tooth disease from mitofusin 2 mutations. J. Neurosci. 27, 422–430. doi: 10.1523/JNEUROSCI.4798-06.2007

Beertsen, W., Heersche, J. N., and Aubin, J. E. (1982). Free and polymerized tubulin in cultured bone cells and Chinese hamster ovary cells: the influence of cold and hormones. J. Cell Biol. 95, 387–393. doi: 10.1083/jcb.95.2.387

Bellier, J. P., and Kimura, H. (2007). Acetylcholine synthesis by choline acetyltransferase of a peripheral type as demonstrated in adult rat dorsal root ganglion. J. Neurochem. 101, 1607–1618. doi: 10.1111/j.1471-4159.2007.04458.x

Bennett, G. J., Doyle, T., and Salvemini, D. (2014). Mitotoxicity in distal symmetrical sensory peripheral neuropathies. Nat. Rev. Neurol. 10, 326–336. doi: 10.1038/nrneurol.2014.77

Berkeley, J. L., Gomeza, J., Wess, J., Hamilton, S. E., Nathanson, N. M., and Levey, A. I. (2001). M1 muscarinic acetylcholine receptors activate extracellular signal-regulated kinase in CA1 pyramidal neurons in mouse hippocampal slices. Mol. Cell. Neurosci. 18, 512–524. doi: 10.1006/mcne.2001.1042

Bernardini, N., Levey, A. I., and Augusti-Tocco, G. (1999). Rat dorsal root ganglia express m1-m4 muscarinic receptor proteins. J. Peripher. Nerv. Syst. 4, 222–232.

Bernardini, N., Tomassy, G. S., Tata, A. M., Augusti-Tocco, G., and Biagioni, S. (2004). Detection of basal and potassium-evoked acetylcholine release from embryonic DRG explants. J. Neurochem. 88, 1533–1539. doi: 10.1046/j.1471-4159.2003.02292.x

Bernstein, B. W., and Bamburg, J. R. (2003). Actin-ATP hydrolysis is a major energy drain for neurons. J. Neurosci. 23, 1–6. doi: 10.1523/JNEUROSCI.23-01-00002.2003

Biederer, T., and Scheiffele, P. (2007). Mixed-culture assays for analyzing neuronal synapse formation. Nat. Protoc. 2, 670–676. doi: 10.1038/nprot.2007.92

Birdsall, N. J., Hulme, E. C., Stockton, J., Burgen, A. S., Berrie, C. P., Hammer, R., et al. (1983). Muscarinic receptor subclasses: evidence from binding studies. Adv. Biochem. Psychopharmacol. 37, 323–329. doi: 10.1016/0165-6147(83)90493-5

Blin, N., Yun, J., and Wess, J. (1995). Mapping of single amino acid residues required for selective activation of Gq/11 by the m3 muscarinic acetylcholine receptor. J. Biol. Chem. 270, 17741–17748. doi: 10.1074/jbc.270.30.17741

Brown, D. A., and Passmore, G. M. (2009). Neural KCNQ (Kv7) channels. Br. J. Pharmacol. 156, 1185–1195. doi: 10.1111/j.1476-5381.2009.00111.x

Butcher, A. J., Bradley, S. J., Prihandoko, R., Brooke, S. M., Mogg, A., Bourgognon, J. M., et al. (2016). An antibody biosensor establishes the activation of the M1 muscarinic acetylcholine receptor during learning and memory. J. Biol. Chem. 291, 8862–8875. doi: 10.1074/jbc.M115.681726

Bymaster, F. P., Carter, P. A., Yamada, M., Gomeza, J., Wess, J., Hamilton, S. E., et al. (2003). Role of specific muscarinic receptor subtypes in cholinergic parasympathomimetic responses, in vivo phosphoinositide hydrolysis, and pilocarpine-induced seizure activity. Eur. J. Neurosci. 17, 1403–1410. doi: 10.1046/j.1460-9568.2003.02588.x

Calcutt, N. A., Smith, D. R., Frizzi, K., Sabbir, M. G., Chowdhury, S. K., Mixcoatl-Zecuatl, T., et al. (2017). Selective antagonism of muscarinic receptors is neuroprotective in peripheral neuropathy. J. Clin. Invest. 127, 608–622. doi: 10.1172/JCI88321

Caulfield, M. P., and Birdsall, N. J. (1998). International union of pharmacology. XVII. Classification of muscarinic acetylcholine receptors. Pharmacol. Rev. 50, 279–290.

Chen, H., and Chan, D. C. (2006). Critical dependence of neurons on mitochondrial dynamics. Curr. Opin. Cell Biol. 18, 453–459. doi: 10.1016/j.ceb.2006.06.004

Chiu, I. M., Barrett, L. B., Williams, E. K., Strochlic, D. E., Lee, S., Weyer, A. D., et al. (2014). Transcriptional profiling at whole population and single cell levels reveals somatosensory neuron molecular diversity. eLife 3:e04660. doi: 10.7554/eLife.04660

Chowdhury, S. K., Smith, D. R., and Fernyhough, P. (2013). The role of aberrant mitochondrial bioenergetics in diabetic neuropathy. Neurobiol. Dis. 51, 56–65. doi: 10.1016/j.nbd.2012.03.016

Corsetti, V., Mozzetta, C., Biagioni, S., Augusti Tocco, G., and Tata, A. M. (2012). The mechanisms and possible sites of acetylcholine release during chick primary sensory neuron differentiation. Life Sci. 91, 783–788. doi: 10.1016/j.lfs.2012.08.026

Courchesne, S. L., Karch, C., Pazyra-Murphy, M. F., and Segal, R. A. (2011). Sensory neuropathy attributable to loss of Bcl-w. J. Neurosci. 31, 1624–1634. doi: 10.1523/JNEUROSCI.3347-10.2011

Delmas, P., and Brown, D. A. (2005). Pathways modulating neural KCNQ/M (Kv7) potassium channels. Nat. Rev. Neurosci. 6, 850–862. doi: 10.1038/nrn1785

Dunn, O. J. (1964). Multiple Comparisons Using Rank Sums. Technometrics 6, 241–252. doi: 10.1080/00401706.1964.10490181

Eglen, R. M. (2005). Muscarinic receptor subtype pharmacology and physiology. Prog. Med. Chem. 43, 105–136. doi: 10.1016/S0079-6468(05)43004-0

Felder, C. C., Porter, A. C., Skillman, T. L., Zhang, L., Bymaster, F. P., Nathanson, N. M., et al. (2001). Elucidating the role of muscarinic receptors in psychosis. Life Sci. 68, 2605–2613. doi: 10.1016/S0024-3205(01)01059-1

Fukuhara, S., Chikumi, H., and Gutkind, J. S. (2001). RGS-containing RhoGEFs: the missing link between transforming G proteins and Rho? Oncogene 20, 1661–1668. doi: 10.1038/sj.onc.1204182

Fukuhara, S., Murga, C., Zohar, M., Igishi, T., and Gutkind, J. S. (1999). A novel PDZ domain containing guanine nucleotide exchange factor links heterotrimeric G proteins to Rho. J. Biol. Chem. 274, 5868–5879. doi: 10.1074/jbc.274.9.5868

Gordois, A., Scuffham, P., Shearer, A., Oglesby, A., and Tobian, J. A. (2003). The health care costs of diabetic peripheral neuropathy in the US. Diabetes Care 26, 1790–1795. doi: 10.2337/diacare.26.6.1790

Grando, S. A., Pittelkow, M. R., and Schallreuter, K. U. (2006). Adrenergic and cholinergic control in the biology of epidermis: physiological and clinical significance. J. Invest. Dermatol. 126, 1948–1965. doi: 10.1038/sj.jid.5700151

Haga, T. (2013). Molecular properties of muscarinic acetylcholine receptors. Proc. Jpn. Acad. Ser. B Phys. Biol. Sci. 89, 226–256. doi: 10.2183/pjab.89.226

Hall, A., and Lalli, G. (2010). Rho and Ras GTPases in axon growth, guidance, and branching. Cold Spring Harb. Perspect. Biol. 2:a001818. doi: 10.1101/cshperspect.a001818

Hall, C. N., Klein-Flugge, M. C., Howarth, C., and Attwell, D. (2012). Oxidative phosphorylation, not glycolysis, powers presynaptic and postsynaptic mechanisms underlying brain information processing. J. Neurosci. 32, 8940–8951. doi: 10.1523/JNEUROSCI.0026-12.2012

Hamilton, S. E., and Nathanson, N. M. (2001). The M1 receptor is required for muscarinic activation of mitogen-activated protein (MAP) kinase in murine cerebral cortical neurons. J. Biol. Chem. 276, 15850–15853. doi: 10.1074/jbc.M011563200

Hanada, K., Kishimoto, S., Bellier, J. P., and Kimura, H. (2013). Peripheral choline acetyltransferase in rat skin demonstrated by immunohistochemistry. Cell Tissue Res. 351, 497–510. doi: 10.1007/s00441-012-1536-z

Hirose, M., Ishizaki, T., Watanabe, N., Uehata, M., Kranenburg, O., Moolenaar, W. H., et al. (1998). Molecular dissection of the Rho-associated protein kinase (p160ROCK)-regulated neurite remodeling in neuroblastoma N1E-115 cells. J. Cell Biol. 141, 1625–1636. doi: 10.1083/jcb.141.7.1625

Hu, J., Wang, Y., Zhang, X., Lloyd, J. R., Li, J. H., Karpiak, J., et al. (2010). Structural basis of G protein-coupled receptor-G protein interactions. Nat. Chem. Biol. 6, 541–548. doi: 10.1038/nchembio.385

Jiang, S., Li, Y., Zhang, C., Zhao, Y., Bu, G., Xu, H., et al. (2014). M1 muscarinic acetylcholine receptor in Alzheimer’s disease. Neurosci. Bull. 30, 295–307. doi: 10.1007/s12264-013-1406-z

Kapitein, L. C., and Hoogenraad, C. C. (2015). Building the Neuronal Microtubule Cytoskeleton. Neuron 87, 492–506. doi: 10.1016/j.neuron.2015.05.046

Karlsson, E., Jolkkonen, M., Mulugeta, E., Onali, P., and Adem, A. (2000). Snake toxins with high selectivity for subtypes of muscarinic acetylcholine receptors. Biochimie 82, 793–806. doi: 10.1016/S0300-9084(00)01176-7

Kellar, K. J., Martino, A. M., Hall, D. P. Jr., Schwartz, R. D., and Taylor, R. L. (1985). High-affinity binding of [3H]acetylcholine to muscarinic cholinergic receptors. J. Neurosci. 5, 1577–1582. doi: 10.1523/JNEUROSCI.05-06-01577.1985

Khan, K. M., Drescher, M. J., Hatfield, J. S., Khan, A. M., and Drescher, D. G. (2002). Muscarinic receptor subtypes are differentially distributed in the rat cochlea. Neuroscience 111, 291–302. doi: 10.1016/S0306-4522(02)00020-9

Kozasa, T., Jiang, X., Hart, M. J., Sternweis, P. M., Singer, W. D., Gilman, A. G., et al. (1998). p115 RhoGEF, a GTPase activating protein for Galpha12 and Galpha13. Science 280, 2109–2111. doi: 10.1126/science.280.5372.2109

Krajewski, J. L., Dickerson, I. M., and Potter, L. T. (2001). Site-directed mutagenesis of m1-toxin1: two amino acids responsible for stable toxin binding to M1 muscarinic receptors. Mol. Pharmacol. 60, 725–731.

Kroeze, W. K., Sassano, M. F., Huang, X. P., Lansu, K., Mccorvy, J. D., Giguere, P. M., et al. (2015). PRESTO-Tango as an open-source resource for interrogation of the druggable human GPCRome. Nat. Struct. Mol. Biol. 22, 362–369. doi: 10.1038/nsmb.3014

Kruse, A. C., Kobilka, B. K., Gautam, D., Sexton, P. M., Christopoulos, A., and Wess, J. (2014). Muscarinic acetylcholine receptors: novel opportunities for drug development. Nat. Rev. Drug Discov. 13, 549–560. doi: 10.1038/nrd4295

Langmead, C. J., Watson, J., and Reavill, C. (2008). Muscarinic acetylcholine receptors as CNS drug targets. Pharmacol. Ther. 117, 232–243. doi: 10.1016/j.pharmthera.2007.09.009

Lau, J. K., Brown, K. C., and Dasgupta, P. (2013). Measurement of Acetylcholine from Cell Lines. Bio Protoc. 3:e1007. doi: 10.21769/BioProtoc.1007

Levey, A. I. (1993). Immunological localization of m1-m5 muscarinic acetylcholine receptors in peripheral tissues and brain. Life Sci. 52, 441–448. doi: 10.1016/0024-3205(93)90300-R

Liu, J., Blin, N., Conklin, B. R., and Wess, J. (1996). Molecular mechanisms involved in muscarinic acetylcholine receptor-mediated G protein activation studied by insertion mutagenesis. J. Biol. Chem. 271, 6172–6178. doi: 10.1074/jbc.271.11.6172

Lovas, J. R., and Wang, X. (2013). The meaning of mitochondrial movement to a neuron’s life. Biochim. Biophys. Acta 1833, 184–194. doi: 10.1016/j.bbamcr.2012.04.007

Luo, W., Latchney, L. R., and Culp, D. J. (2001). G protein coupling to M1 and M3 muscarinic receptors in sublingual glands. Am. J. Physiol. Cell Physiol. 280, C884–C896. doi: 10.1152/ajpcell.2001.280.4.C884

Max, S. I., Liang, J. S., and Potter, L. T. (1993a). Purification and properties of m1-toxin, a specific antagonist of m1 muscarinic receptors. J. Neurosci. 13, 4293–4300.

Max, S. I., Liang, J. S., and Potter, L. T. (1993b). Stable allosteric binding of m1-toxin to m1 muscarinic receptors. Mol. Pharmacol. 44, 1171–1175.

McInnes, A. D. (2012). Diabetic foot disease in the United Kingdom: about time to put feet first. J. Foot Ankle Res. 5:26. doi: 10.1186/1757-1146-5-26

Meijering, E., Dzyubachyk, O., and Smal, I. (2012). Methods for cell and particle tracking. Methods Enzymol. 504, 183–200. doi: 10.1016/B978-0-12-391857-4.00009-4

Mironov, S. L. (2007). ADP regulates movements of mitochondria in neurons. Biophys. J. 92, 2944–2952. doi: 10.1529/biophysj.106.092981

Misko, A., Jiang, S., Wegorzewska, I., Milbrandt, J., and Baloh, R. H. (2010). Mitofusin 2 is necessary for transport of axonal mitochondria and interacts with the Miro/Milton complex. J. Neurosci. 30, 4232–4240. doi: 10.1523/JNEUROSCI.6248-09.2010

Nguyen, V. T., Chernyavsky, A. I., Arredondo, J., Bercovich, D., Orr-Urtreger, A., Vetter, D. E., et al. (2004). Synergistic control of keratinocyte adhesion through muscarinic and nicotinic acetylcholine receptor subtypes. Exp. Cell Res. 294, 534–549. doi: 10.1016/j.yexcr.2003.12.010

O’Brien, E. T., Salmon, E. D., and Erickson, H. P. (1997). How calcium causes microtubule depolymerization. Cell Motil. Cytoskeleton 36, 125–135. doi: 10.1002/(SICI)1097-0169(1997)36:2<125::AID-CM3>3.0.CO;2-8

Pavia, J., Marquez, E., Laukkonen, S., Martos, F., Gomez, A., and Sanchez de la Cuesta, F. (1991). M1 and M3 muscarinic receptor subtypes in rat forebrain. Methods Find. Exp. Clin. Pharmacol. 13, 653–660.

Popova, J. S., and Rasenick, M. M. (2003). G beta gamma mediates the interplay between tubulin dimers and microtubules in the modulation of Gq signaling. J. Biol. Chem. 278, 34299–34308. doi: 10.1074/jbc.M301748200

Ravindra, R., Kunapuli, S. P., Forman, L. J., Nagele, R. G., Foster, K. A., and Patel, S. A. (1996). Effect of transient overexpression of Gq alpha on soluble and polymerized tubulin pools in GH3 and AtT-20 cells. J. Cell. Biochem. 61, 392–401. doi: 10.1002/(SICI)1097-4644(19960601)61:3<392::AID-JCB6>3.0.CO;2-O

Roland, A. B., Ricobaraza, A., Carrel, D., Jordan, B. M., Rico, F., Simon, A., et al. (2014). Cannabinoid-induced actomyosin contractility shapes neuronal morphology and growth. eLife 3:e03159. doi: 10.7554/eLife.03159

Roy Chowdhury, S. K., Smith, D. R., Saleh, A., Schapansky, J., Marquez, A., Gomes, S., et al. (2012). Impaired adenosine monophosphate-activated protein kinase signalling in dorsal root ganglia neurons is linked to mitochondrial dysfunction and peripheral neuropathy in diabetes. Brain 135, 1751–1766. doi: 10.1093/brain/aws097

Roychowdhury, S., Martinez, L., Salgado, L., Das, S., and Rasenick, M. M. (2006). G protein activation is prerequisite for functional coupling between Galpha/Gbetagamma and tubulin/microtubules. Biochem. Biophys. Res. Commun. 340, 441–448. doi: 10.1016/j.bbrc.2005.12.026

Roychowdhury, S., Panda, D., Wilson, L., and Rasenick, M. M. (1999). G protein alpha subunits activate tubulin GTPase and modulate microtubule polymerization dynamics. J. Biol. Chem. 274, 13485–13490. doi: 10.1074/jbc.274.19.13485

Roychowdhury, S., and Rasenick, M. M. (1997). G protein beta1gamma2 subunits promote microtubule assembly. J. Biol. Chem. 272, 31576–31581. doi: 10.1074/jbc.272.50.31576

Roychowdhury, S., and Rasenick, M. M. (2008). Submembraneous microtubule cytoskeleton: regulation of microtubule assembly by heterotrimeric Gproteins. FEBS J. 275, 4654–4663. doi: 10.1111/j.1742-4658.2008.06614.x

Sabbir, M. G., Dillon, R., and Mowat, M. R. (2016). Dlc1 interaction with non-muscle myosin heavy chain II-A (Myh9) and Rac1 activation. Biol. Open 15, 452–460. doi: 10.1242/bio.015859

Sajic, M., Mastrolia, V., Lee, C. Y., Trigo, D., Sadeghian, M., Mosley, A. J., et al. (2013). Impulse conduction increases mitochondrial transport in adult mammalian peripheral nerves in vivo. PLoS Biol. 11:e1001754. doi: 10.1371/journal.pbio.1001754

Saleh, A., Roy Chowdhury, S. K., Smith, D. R., Balakrishnan, S., Tessler, L., Martens, C., et al. (2013). Ciliary neurotrophic factor activates NF-kappaB to enhance mitochondrial bioenergetics and prevent neuropathy in sensory neurons of streptozotocin-induced diabetic rodents. Neuropharmacology 65, 65–73. doi: 10.1016/j.neuropharm.2012.09.015

Sarma, T., Koutsouris, A., Yu, J. Z., Krbanjevic, A., Hope, T. J., and Rasenick, M. M. (2015). Activation of microtubule dynamics increases neuronal growth via the nerve growth factor (NGF)- and Galphas-mediated signaling pathways. J. Biol. Chem. 290, 10045–10056. doi: 10.1074/jbc.M114.630632

Sarma, T., Voyno-Yasenetskaya, T., Hope, T. J., and Rasenick, M. M. (2003). Heterotrimeric G-proteins associate with microtubules during differentiation in PC12 pheochromocytoma cells. FASEB J. 17, 848–859. doi: 10.1096/fj.02-0730com

Schappi, J. M., Krbanjevic, A., and Rasenick, M. M. (2014). Tubulin, actin and heterotrimeric G proteins: coordination of signaling and structure. Biochim. Biophys. Acta 1838, 674–681. doi: 10.1016/j.bbamem.2013.08.026

Schindelin, J., Arganda-Carreras, I., Frise, E., Kaynig, V., Longair, M., Pietzsch, T., et al. (2012). Fiji: an open-source platform for biological-image analysis. Nat. Methods 9, 676–682. doi: 10.1038/nmeth.2019

Schlereth, T., Birklein, F., an Haack, K., Schiffmann, S., Kilbinger, H., Kirkpatrick, C. J., et al. (2006). In vivo release of non-neuronal acetylcholine from the human skin as measured by dermal microdialysis: effect of botulinum toxin. Br. J. Pharmacol. 147, 183–187. doi: 10.1038/sj.bjp.0706451

Schneider, C. A., Rasband, W. S., and Eliceiri, K. W. (2012). NIH Image to ImageJ: 25 years of image analysis. Nat. Methods 9, 671–675. doi: 10.1038/nmeth.2089

Schwarz, T. L. (2013). Mitochondrial trafficking in neurons. Cold Spring Harb. Perspect. Biol. 5:a011304. doi: 10.1101/cshperspect.a011304

Shaw, G., Morse, S., Ararat, M., and Graham, F. L. (2002). Preferential transformation of human neuronal cells by human adenoviruses and the origin of HEK 293 cells. FASEB J. 16, 869–871. doi: 10.1096/fj.01-0995fje

Sheng, Z. H., and Cai, Q. (2012). Mitochondrial transport in neurons: impact on synaptic homeostasis and neurodegeneration. Nat. Rev. Neurosci. 13, 77–93. doi: 10.1038/nrn3156

Siatkowski, R. M., Cotter, S. A., Crockett, R. S., Miller, J. M., Novack, G. D., and Zadnik, K. (2008). Two-year multicenter, randomized, double-masked, placebo-controlled, parallel safety and efficacy study of 2% pirenzepine ophthalmic gel in children with myopia. J. AAPOS 12, 332–339. doi: 10.1016/j.jaapos.2007.10.014

Siegel, S. (ed.) (1956). Nonparametric Statistics for the Behavioral Sciences. New York, NY: McGraw-Hill.

Siehler, S. (2009). Regulation of RhoGEF proteins by G12/13-coupled receptors. Br. J. Pharmacol. 158, 41–49. doi: 10.1111/j.1476-5381.2009.00121.x

Stepanenko, A. A., and Dmitrenko, V. V. (2015). HEK293 in cell biology and cancer research: phenotype, karyotype, tumorigenicity, and stress-induced genome-phenotype evolution. Gene 569, 182–190. doi: 10.1016/j.gene.2015.05.065

Suzuki, N., Nakamura, S., Mano, H., and Kozasa, T. (2003). Galpha 12 activates Rho GTPase through tyrosine-phosphorylated leukemia-associated RhoGEF. Proc. Natl. Acad. Sci. U.S.A. 100, 733–738. doi: 10.1073/pnas.0234057100

Tata, A. M., Cursi, S., Biagioni, S., and Augusti-Tocco, G. (2003). Cholinergic modulation of neurofilament expression and neurite outgrowth in chick sensory neurons. J. Neurosci. Res. 73, 227–234. doi: 10.1002/jnr.10650

Tata, A. M., Tripiciano, A., Filippini, A., Biagioni, S., and Augusti-Tocco, G. (2000a). Muscarinic receptors modulate intracellular calcium level in chick sensory neurons. Brain Res. 866, 65–72.

Tata, A. M., Vilaro, M. T., and Mengod, G. (2000b). Muscarinic receptor subtypes expression in rat and chick dorsal root ganglia. Brain Res. Mol. Brain Res. 82, 1–10.

Tsai, F. C., Kuo, G. H., Chang, S. W., and Tsai, P. J. (2015). Ca²⁺ signaling in cytoskeletal reorganization, cell migration, and cancer metastasis. Biomed Res. Int. 2015:409245. doi: 10.1155/2015/409245

Uhlen, M., Bjorling, E., Agaton, C., Szigyarto, C. A., Amini, B., Andersen, E., et al. (2005). A human protein atlas for normal and cancer tissues based on antibody proteomics. Mol. Cell. Proteomics 4, 1920–1932. doi: 10.1074/mcp.M500279-MCP200

Van Kuilenburg, A. B., Muijsers, A. O., Demol, H., Dekker, H. L., and Van Beeumen, J. J. (1988). Human heart cytochrome c oxidase subunit VIII. Purification and determination of the complete amino acid sequence. FEBS Lett. 240, 127–132. doi: 10.1016/0014-5793(88)80353-3

Verburg, J., and Hollenbeck, P. J. (2008). Mitochondrial membrane potential in axons increases with local nerve growth factor or semaphorin signaling. J. Neurosci. 28, 8306–8315. doi: 10.1523/JNEUROSCI.2614-08.2008

Wang, N., Yan, K., and Rasenick, M. M. (1990). Tubulin binds specifically to the signal-transducing proteins, Gs alpha and Gi alpha 1. J. Biol. Chem. 265, 1239–1242.

Wess, J. (1996). Molecular biology of muscarinic acetylcholine receptors. Crit. Rev. Neurobiol. 10, 69–99. doi: 10.1615/CritRevNeurobiol.v10.i1.40

Wess, J., Duttaroy, A., Zhang, W., Gomeza, J., Cui, Y., Miyakawa, T., et al. (2003). M1-M5 muscarinic receptor knockout mice as novel tools to study the physiological roles of the muscarinic cholinergic system. Recept. Channels 9, 279–290. doi: 10.3109/10606820308262

Wess, J., Eglen, R. M., and Gautam, D. (2007). Muscarinic acetylcholine receptors: mutant mice provide new insights for drug development. Nat. Rev. Drug Discov. 6, 721–733. doi: 10.1038/nrd2379

Yamada, T., Inazu, M., Tajima, H., and Matsumiya, T. (2011). Functional expression of choline transporter-like protein 1 (CTL1) in human neuroblastoma cells and its link to acetylcholine synthesis. Neurochem. Int. 58, 354–365. doi: 10.1016/j.neuint.2010.12.011

Yang, H., and Kunes, S. (2004). Nonvesicular release of acetylcholine is required for axon targeting in the Drosophila visual system. Proc. Natl. Acad. Sci. U.S.A. 101, 15213–15218. doi: 10.1073/pnas.0308141101

Zerial, M., and McBride, H. (2001). Rab proteins as membrane organizers. Nat. Rev. Mol. Cell Biol. 2, 107–117. doi: 10.1038/35052055

Zhou, B., Yu, P., Lin, M. Y., Sun, T., Chen, Y., and Sheng, Z. H. (2016). Facilitation of axon regeneration by enhancing mitochondrial transport and rescuing energy deficits. J. Cell Biol. 214, 103–119. doi: 10.1083/jcb.201605101