Abstract
Nicotinic acetylcholine receptors (nAChRs) support the initiation and maintenance of smoking, but the long-term changes occurring in the protein complex as a result of smoking and the nicotine in tobacco are not known. Human studies and animal models have also demonstrated that increasing cholinergic tone increases behaviors related to depression, suggesting that the nAChR-associated proteome could be altered in individuals with mood disorders. We therefore immunopurified nAChRs and associated proteins for quantitative proteomic assessment of changes in protein–protein interactions of high-affinity nAChRs containing the β2 subunit (β2*-nAChRs) from either cortex of mice treated with saline or nicotine, or postmortem human temporal cortex tissue from tobacco-exposed and nonexposed individuals, with a further comparison of diagnosed mood disorder to control subjects. We observed significant effects of nicotine exposure on the β2*-nAChR-associated proteome in human and mouse cortex, particularly in the abundance of the nAChR subunits themselves, as well as putative interacting proteins that make up core components of neuronal excitability (Na/K ATPase subunits), presynaptic neurotransmitter release (syntaxins, SNAP25, synaptotagmin), and a member of a known nAChR protein chaperone family (14-3-3ζ). These findings identify candidate-signaling proteins that could mediate changes in cholinergic signaling via nicotine or tobacco use. Further analysis of identified proteins will determine whether these interactions are essential for primary function of nAChRs at presynaptic terminals. The identification of differences in the nAChR-associated proteome and downstream signaling in subjects with various mood disorders may also identify novel etiological mechanisms and reveal new treatment targets.
Significance Statement
Nicotinic acetylcholine receptors (nAChRs) are ligand-gated ion channels responsible for rapid excitatory acetylcholine signaling, which modulate numerous circuits involved in complex behaviors. In addition, nAChRs are the primary targets for the reinforcing effects of nicotine in tobacco products. Proteins associated with nAChRs are critical for receptor trafficking, localization, and function. Identifying the core elements of the nAChR-associated proteome that are conserved across human and mouse brain is essential for understanding how these receptors are regulated, and how their functions and interactions are altered by long-term exposure to nicotine.
Introduction
Individuals who smoke have significantly higher numbers of high-affinity nicotinic acetylcholine receptors (nAChRs) compared to nonsmokers, as measured in both postmortem brain studies (Breese et al., 1997) and SPECT imaging studies with a radiotracer specific for nAChRs containing the β2 subunit (β2*-nAChRs; Staley et al., 2006). This upregulation appears to be due primarily to the nicotine in tobacco, since rodent studies demonstrate that long-term exposure to nicotine alone results in a similar upregulation of β2*-nAChRs in brain (Marks et al., 1992). While this change in receptor number can alter signaling through nAChRs (Buisson and Bertrand, 2001; Govind et al., 2012), the effect of nicotine treatment and upregulation on the associated proteome is not known. Identifying the complex associated with upregulated nAChRs would be helpful in understanding both the basic biology of the nAChR system, as well as the consequences of smoking on nAChR-induced intracellular signaling. Comparison of the nAChR-associated proteome in human smokers versus rodents exposed to nicotine is also important to identify molecular changes selective for nicotine alone, rather than the ∼4000 other constituents of tobacco smoke. Further, the identification of the changes in the nAChR-associated proteome as a result of nicotine exposure could provide insights into the mechanisms underlying nAChR upregulation, as well as determining whether nicotine-mediated upregulation alters the downstream signaling of these receptors. Finally, comparison of the human and rodent nAChR-associated proteome at baseline could identify conserved signaling pathways that important for ongoing smoking behavior and other disorders associated with changes in nAChR signaling.
Among disorders associated with smoking, mood disorders (MDs) are of particular interest because of the high comorbidity reported (Lerman et al., 1996; Laje et al., 2001; Mineur and Picciotto, 2009) and because several lines of evidence implicate the dysfunction of the cholinergic system in depressive states (for review, see Mineur and Picciotto, 2010; Philip et al., 2010). Animal and human studies suggest that increased ACh activity can induce symptoms of depression, whereas decreased cholinergic activity may be associated with mania (Dilsaver, 1986; Dilsaver and Coffman, 1989). Consistent with challenge studies using a cholinesterase inhibitor to increase nAChR occupancy (Esterlis et al., 2013), imaging studies suggest that brain ACh levels may be higher in individuals with major depressive or bipolar disorder when they are actively depressed (Saricicek et al., 2012; Hannestad et al., 2013). With respect to nAChRs, decreasing the activity of β2*-nAChRs using molecular genetics, pharmacological antagonism, or partial agonism can alter behavior in tests of antidepressant efficacy (Caldarone et al., 2004; Rabenstein et al., 2006; Mineur et al., 2007, 2009). Overall nAChR signaling may be altered in individuals with depressive disorders, suggesting that it could be important to identify any potential differences in the functional regulation of β2*-nAChRs and protein–protein interactions in patients with affective disorders.
In the current study, we took a dual approach to identify changes in the nAChR-associated proteome with nicotine exposure and to determine whether any changes were also apparent in individuals with a mood disorder. We obtained postmortem tissue from the temporal cortex of human subjects and used tissue cotinine content as a determinant of tobacco use to determine whether there was any difference in the nAChR-associated proteome as a function of receptor number. (Note that we use “tobacco use” rather than “smoker” since it is not possible to discriminate between modes of tobacco use from brain metabolite levels). Among the human subjects, we evaluated cortical tissue from both healthy control subjects and matched individuals with a mood disorder diagnosis to identify any differences in receptor number or response to smoking. In addition, we exposed mice with differing levels of high-affinity nAChRs to nicotine in the long term to identify changes in the associated proteome that do or do not scale with nAChR number to determine which smoking-induced changes in human tissue can be attributed to the nicotine in tobacco. These results are the first to compare human and mouse nAChR-associated proteins, and the findings identify conserved proteins in the β2*-nAChR interactome and the effects of nicotine exposure that are valid to explore with other models.
Materials and Methods
Human temporal cortex tissue collection
Human temporal cortex samples were obtained from the Douglas Brain Bank as frozen tissue sections, ranging from 0.29 to 1.03 g wet weight of tissue. Visual inspection of each section indicated approximately equivalent amounts of gray and white matter per section. All tissue samples were numbered by the supplier, with diagnosis (major depressive, bipolar, mood disorder not otherwise specified, or no mood disorder) and cause of death as the only additional information. Experimental results were not sorted by diagnosis until after results were obtained; experimenters were blind to the groups during data collection.
Mice
Wild-type mice and mice carrying genetic deletion of the β2 nAChR subunit (Picciotto et al., 1995, 1998) or the α4 nAChR subunit (Ross et al., 2000) were used for these studies. The animals used were wild type (β2+/+), heterozygous (β2+/−), α4/β2 double-heterozygous (α4+/−β2+/−), and knock-out (β2−/−) mice. Mice were bred and maintained in the animal colonies at the University of Colorado and Yale University. Female mice were used for this study because females are more resilient to chronic nicotine administration (Hatchell and Collins, 1980) and were housed in groups of no more than five same-sex littermates per cage, on a 12 h light/dark cycle with ad libitum access to food and water. All protocols involving the use of live animals were approved by and followed the guidelines of the institutional animal care and use committees of the University of Colorado and Yale University.
Surgery and long-term nicotine treatment
Mice were fitted with jugular catheters essentially as described previously (Marks et al., 2011), with minor modifications. Briefly, mice were anesthetized with an intraperitoneal injection of ketamine/xylazine (100/8 mg/kg) and a small incision (<1 cm) was made above the right clavicle to access the superficial aspect of the jugular vein. A small catheter (inner diameter, 0.51 mm; outer diameter, 0.94 mm) made of medical-grade silastic tubing and filled with sterile saline solution was inserted into the exposed vein and anchored in position with two surgical silk sutures. The catheter was passed under the skin through a second small (<1 cm) incision made in the midthoracic region of the skin and was attached to a stainless steel post anchored under the skin with veterinary adhesive. Incisions were cleaned with 0.1% iodine in ethanol and closed with veterinary adhesive. Mice were given an intraperitoneal injection of buprenorphine (0.1 mg/kg) following surgery and were allowed to recover in warm (37°C) individual cages. For long-term infusion, mice were housed individually, and the jugular catheter was attached to a 1 ml syringe pump and infused with saline at a rate of 35 μl/h for 48 h. Vehicle control mice received sterile saline, and nicotine-treated animals received 4 mg/kg/h nicotine (free base, neutralized with HCl) for a total of 10 d. On the last day of nicotine treatment, mice were detached from the syringe pump and allowed to clear the nicotine for 2 h. Mice were killed by cervical dislocation, and the brain was rapidly removed and placed on an ice-cold surface for dissection. Each brain was dissected, and the cortex was frozen at −80°C until analysis.
[125I]-Epibatidine binding
Quantification of high-affinity nAChR binding sites in brain tissue was performed essentially as described previously (Hannestad et al., 2013) with minor modifications. Briefly, a small prism of tissue was removed from the large frozen section by dissection on an ice-cold surface with a sterile scalpel. The excised tissue was placed immediately in 0.25 ml ice-cold 0.1× phosphate buffer (PB; 12.8 mm NaCl, 0.24 mm KCl, 0.32 mm CaCl2, 0.12 mm KH2PO4, 0.12 mm MgSO4.7H2O, 2.5 mm HEPES hemi-Na, pH7.5). Samples were homogenized by hand in a glass/Teflon tissue grinder and diluted to a final volume of 0.5 ml in 0.1× PB. The tissue homogenate was centrifuged at 10,000 x g at 4°C for 20 min, and the supernatant was discarded. The tissue pellet was resuspended in 0.5 ml of 0.1× PB and centrifuged for 20 min at 10,000 × g at 4°C. This wash step was repeated a total of three times in order to ensure the removal of substances that could compete with [125]I-epibatidine binding (i.e., nicotine, ACh). For ligand binding to particulate fractions, 10 μl aliquots of tissue homogenate were added to wells of 96-well plates containing 200 pm [125I]-epibatidine as the final concentration (Whiteaker et al., 2000). Nonspecific binding was assessed in the presence of 100 μm cytisine. The protein content of each sample was determined using the method of Lowry et al. (1951).
Immunopurification of β2*-nAChRs
Purification of β2*-nAChRs with mAb295-coupled M270 Dynabeads from mouse and human cortical tissue samples was conducted essentially as described previously (McClure-Begley et al., 2013), with minor modifications. For the mouse cortex samples, dissected cortices (both hemispheres) from two mice were pooled to create one homogenate for immunoprecipitation (IP); human samples were used as is, and the volume of homogenization buffer was adjusted to equal 4 volumes based on the total wet weight of the tissue. Briefly, brain tissue sections were homogenized by hand with 17 strokes in a glass/Teflon tissue grinder with 4 volumes of ice-cold PBS (pH 8.0), and placed on ice. To each tissue homogenate, dithiobis[succinimidyl propionate] [DSP (also called Lamont's Reagent); 12 Å linker span; 10 mm in DMSO] was added to a final concentration of 1 mm to stabilize labile interactions prior to detergent extraction. Tissue homogenates were incubated on ice for 45 min with DSP prior to the addition of Tris to a final concentration of 10 mm to quench the crosslinking reaction. Samples were centrifuged at 10,000 × g at 4°C for 20 min, and the supernatant was collected and frozen at −20°C, and the pellet was washed by resuspension in TBS before being centrifuged again. The washed pellet was resuspended in 2 ml of extraction buffer (EB; in mm: 121.9 NaCl, 2.68 KCl, 10.14 Na2HPO4, 1.76 KH2PO4, 5 EDTA, 5 EGTA, 5 NaF, 0.1 Na3VO4, 1.0 PMSF, 10 μg/ml each of aprotonin, prepstatin A, leupeptin, 0.and 6% Triton X-100, pH 7.4) and incubated for 30 min with gentle rotation at room temperature. The extract was clarified by centrifugation at 10,000 × g at 4°C for 20 min, and the supernatant was collected. A 100 μl aliquot of clarified extract was taken at this point to quantify total extracted nAChR content with [125I]-epibatidine binding, as described below. The remaining extract was added to 5 ml of polystyrene culture tubes containing a volume of mAb295-Dynabeads equivalent to 50% of the input volume and incubated with gentle rotation at 4°C overnight. The following morning, sample tubes were placed on a magnetic stand to separate the mAb295 beads, and a 100 μl aliquot of the supernatant was collected to assess nAChR capture efficiency. For the quantification of nAChR immunopurification, 20 μl aliquots of input and post-IP tissue lysate were labeled with 200 pm [125I]-epibatidine overnight at 4°C, and protein precipitated by 50 μl of 40% w/v polyethylene glycol (PEG) added to each well. PEG-precipitated binding sites were captured on filter membranes as for the particulate binding studies. The depletion of binding sites was used to estimate receptor content of each sample post-elution.
Sample preparation and protein identification by liquid chromatography-tandem mass spectrometry
Dynabeads carrying mAb295-captured nAChRs and bound proteins were washed two times with 1 ml of PBS (137 mm NaCl, 27 mm KCl, 100 mm Na2HPO4, 18 mm KH2PO4, pH 7.5) containing 0.01% Tween-20, and one time with PBS prior to elution with 0.5 ml of 0.1 m NH4OH and 1 mm EDTA, pH 10.0. Elution was performed twice, and eluates were pooled prior to lyophilization in a Thermo Savant SPD1010 Speedvac. Eluted proteins were reduced, alkylated, and digested as described previously (McClure-Begley et al., 2013). Mouse cortex samples were labeled in an eight-plex iTRAQ labeling scheme, and peptide sequencing was performed on an ABSciex 5600 triple- time-of-flight mass spectrometer, with iTRAQ quantitation in ProteinPilot software, as has been described (McClure-Begley et al., 2013). For the human temporal cortex samples, the eluted protein pellet was redissolved in 30 µl of 8 m urea/0.4 m ammonium bicarbonate, pH 8.0. The addition of 3 µl of 45 mm DTT was followed by vortexing, and samples were then incubated at 37°C for 20 min. Reduced samples were cooled to room temperature and alkylated with 3 µl of 100 mm iodoacetamide for 20 min in the dark. Water (74 µl) was added to dilute the sample prior to trypsin digestion (10 µl, 0.5 mg/ml) at 37°C overnight. Samples were desalted on a Macrospin cartridge (Nest Group), and eluent was dried in a Speedvac. The pellet was resuspended in 10 µl of 70% formic acid and 40 µl of 50 mm sodium PB. Equal peptide amounts were aliquoted into glass vials for injection.
Proteomic analysis was performed using a label-free quantification (LFQ) method on an Orbitrap Elite equipped with a Waters Symmetry C18 (180 µm × 20 mm) trap column and a 1.7 µm, 75 µm × 250 mm nanoAcquity ultra high-performance liquid chromatography (UPLC) column (35°C). Trapping was performed using 99% Buffer A (99.9% water, 0.1% formic acid) and peptide separation with a linear gradient of solvent A (0.1% formic acid in water) and solvent B (0.075% formic acid in acetonitrile) over 90 min, at a flow rate of 300 nl/min. MS spectra were acquired in the Orbitrap using a 1 µm scan and a maximum injection time of 900 ms followed by three data-dependant MS/MS acquisitions in the ion trap (with precursor ions threshold of >3000). The total cycle time for both MS and MS/MS acquisition was 2.4 s. Peaks targeted for MS/MS fragmentation by collision-induced dissociation were first isolated with a 2 Da window followed by normalized collision energy of 35%. Dynamic exclusion was activated where former target ions were excluded for 30 s.
Feature extraction, chromatographic/spectral alignment, data filtering, and statistical analysis were performed using Nonlinear Dynamics Progenesis LC-MS software (www.nonlinear.com). The raw data files were imported into the program, a sample run was chosen as a reference (usually at or near the middle of all runs in a set), and all other runs were automatically aligned to that run in order to minimize retention time (RT) variability between runs. All runs were selected for the detection of peptides with an automatic detection limit for signal-to-noise ratio of 3. Features within RT ranges of 0–25 min and 110–120 min were filtered out (these are outside the range of the peptide elution profile), as were features with a charge of ≥7. A normalization factor was then calculated for each run to account for differences in sample load between injections. The experimental design was set up to group multiple injections from each run. The algorithm then calculated and tabulated the raw and normalized abundances, the maximum fold changes and the ANOVA values for each feature in the dataset. The corresponding MS and MS/MS spectra for the features were exported as a Mascot generic file (.mgf) for database searching to perform protein identification. The .mgf files created by the Progenesis LCMS were searched in-house using the Mascot Search Algorithm (version 2.2.0). The data were searched using the Swiss protein database along with the following search parameters: enzyme, trypsin; variable modifications, carbamidomethyl (Cys), oxidation (Met), and alkylated DSP crosslinker-modified lysine (C10H19N2O3S2); peptide mass tolerance, ±20 ppm; fragment mass tolerance, ±0.6 Da; charge, +7; maximum missed cleavages, 3; decoy, yes; and instrument type, ESI-TRAP. The resulting search was exported as an .xml file and read into the Progenesis LC-MS software, where search hits, with at least 95% confidence, were matched to corresponding features (precursor ions and abundances). Features/peptides with nonconflicting identification match were used to provide relative quantification of protein expression. The identification of peptides/proteins and quantitative values across the various comparisons were exported as a .csv file, and further downstream data analyses were performed. The means of normalized abundances for proteins identified in three technical replicates for each tissue sample were used in subsequent analyses.
Cotinine determination from brain tissue lysate by liquid chromatography-tandem mass spectrometry
Following detergent extraction from tissue homogenates and immunopurification of β2*-nAChR complexes with mAb295-coupled Dynabeads, the clarified lysate was frozen at −20°C. The lysate was thawed, and a 200 μl aliquot of each sample was subjected to organic extraction according to published methods (Massadeh et al., 2009) with minor modification. Briefly, the clarified tissue lysate was alkalized with 50 μl of 1 m NaOH and mixed with 2.5 ml of 1:1 dichloromethane/hexanes in a glass test tube and vortexed vigorously. The mixture was centrifuged for 10 min at 5000 × g at room temperature, and the upper organic layer was removed and placed in a clean borosilicate glass tube, then dried under a gentle stream of nitrogen at room temperature. Samples were reconstituted in 100 μl of 0.1% formic acid in water, and 10 μl of each sample was used per injection for analysis.
All liquid chromatography-tandem mass spectrometry (LC-MS/MS) identifications of cotinine were performed on an ABSciex 4000 QTrap Mass Spectrometer in positive ion mode coupled to a PerkinElmer Flexar FX-15 UPLC unit. For chromatography, solvent A was 99.9% water and 0.1% formic acid, and solvent B was 99.9% acetonitrile and 0.1% formic acid. The elution gradient was as follows: 98% solvent A, 2% solvent B from 0 to 1.8 min; 98% solvent B, 2% solvent A from 1.8 to 2.8 min; and 98% solvent A, 2% solvent B from 2.9 to 6.0 min. A selected reaction monitoring (SRM) method was developed for the quantification of cotinine present in the sample using an isotopically coded analytical standard (d-cotinine; Cerilliant Corporation). As described previously (Jacob et al., 2011), the cotinine parent ion had a mass-to-charge ratio (m/z) of 177.12 at charge +1 and yielded a specific transition with major product ions of an m/z values of 80.03 and 98.02. Standard curves showed that intensities of both major product ions were linear with respect to injected cotinine >9.5 ng/μl. Cotinine abundances for each transition from replicate injections were averaged and normalized to the total wet weight of each tissue sample in order to yield units of nanograms per gram tissue to determine nicotine use status.
Data analysis
SigmaPlot 13.0 and SPSS 23 (2015) were used for data organization and statistical analysis. All ligand binding data are reported in specific femtomoles/milligram protein. For proteins identified with LC-MS/MS, we took the sum of the normalized protein abundance for all proteins scoring above the 95% false discovery rate (FDR) cutoff from three technical replicates per sample and generated arithmetic means that were then log2 transformed to yield a continuous variable with a normal distribution (Smits et al., 2013). Linear regression analysis was performed using a Pearson’s product moment correlation coefficient (r) as an indicator of significant positive correlations. Normalization using the calculated abundance of the β2 nAChR subunit (target of the IP) in each sample as the denominator was performed, and ratios relative to control conditions were compared to assess alterations in associated proteins per unit of nAChR.
Mouse cortex iTRAQ data were processed as has been described previously (McClure-Begley et al., 2013). Briefly, mean values of the reporter ion peak areas for every protein identified above the FDR cutoff were divided by the mean value for the control condition (β2+/+, saline) for each set and then log2 transformed to yield normalized protein expression ratios. Significant effects of nicotine exposure and mouse genotype (for the mouse samples) or mood disorder diagnosis and nicotine use (human samples) were determined with ANOVA, Holm–Sidak test, or Student’s t test, where appropriate as determined by the nature of the experiment, with a 95% confidence interval. A general linear model (GLM) split by protein identification with multiple univariate comparisons was used to determine whether levels of individual human brain proteins [a normally distributed (log2 fold change) continuous response variable] that vary over multiple categorical predictor variables (nicotine, MD, both, or neither) differed in representation compared with the control condition (where the intercept is the mean of the no-nicotine exposure, no-MD diagnosis values). This was only possible for human cortical samples that were run using LFQ, in which there is a discrete interval measurement for every group. Mouse samples were run on iTRAQ in which every eight-plex experiment uses one channel as the control value, which is convoluted into the variance of the other seven reporter ion channels. Thus, the mouse GLM does not include a discrete value for the control group, since reintroducing the control group drags the regression to the origin and biases the estimates of significance.
Data are reported as mean ± SEM, with error terms specifically reported for each value.
Results
Effects of smoking and mood disorder diagnosis on nAChR number measured with [125I]-epibatidine binding
We measured high-affinity [125I]-epibatidine binding to particulate fractions in order to determine any effects of smoking status and mood disorder diagnosis on nAChR number in human temporal cortex. High levels of α4β2*-nAChRs, the subtype upregulated by concentrations of nicotine typically encountered by smokers (Teaktong et al., 2004), are found in mammalian cortex (Zoli et al., 1998). Smoking status was not included with human tissue samples, so in order to group samples appropriately, we examined the cotinine content of each sample using an SRM LC-MS/MS method with a deuterated cotinine standard. Samples were grouped according to cotinine content and α4β2*-nAChR expression, with “nicotine users” considered as those samples with >15 ng cotinine/g wet tissue weight. Cotinine levels present in tissue did not correlate linearly with [125I]-epibatidine binding sites (r 2 = 0.33), indicating that the postmortem tissue content of the major nicotine metabolite is useful principally for ordinal assignment of nicotine use (nicotine user or nonuser), and cannot be used to further discriminate “heavy/moderate/light” nicotine users. Based on these criteria, there were nine nicotine users and nine nonusers. Consistent with the high comorbidity between smoking and depression, the incidence of nicotine exposure was higher in individuals with a diagnosed mood disorder (six of nine were nicotine users) compared with individuals with no psychiatric diagnosis (three of nine were nicotine users). We observed a significant effect of cotinine level, our empirically determined tobacco use criterion, on nAChR number, as measured by epibatidine binding in homogenates (two-way ANOVA with Holm–Sidak multiple comparisons: F(1,17) = 9.794, p = 0.007). Analysis revealed no significant effect of mood disorder diagnosis on nAChR number (F(1,17) = 0.224, p = 0.673) or an interaction of nicotine × mood disorder (F(1,17) = 0.04, p = 0.836) measured by equilibrium ligand binding (Fig. 1A), which is in agreement with previous assessments of nAChR number in postmortem human cortical tissue from individuals with major depression or bipolar disorder (Saricicek et al., 2012; Hannestad et al., 2013).
Immunoprecipitation of β2*-nAChRs from human temporal cortex and identification of specific interacting proteins by tandem mass spectrometry
Following the extraction of nAChR complexes from human brain tissue, we quantified specific mAb295 immunocapture of β2*-nAChRs by measuring total [125I]-epibatidine binding sites in the solubilized tissue extract before and after incubation with mAb295-Dynabeads. In human temporal cortex samples, the recovery of β2*-nAChRs using mAb295-Dynabeads was 74.5 ± 7.6% of input across all samples, with no effect of either smoking status or mood disorder diagnosis (t test, p = 0.32, t = −1.017, df = 16). In mouse cortex samples, despite fractional differences in the percentage of upregulation induced by long-term nicotine exposure as a function of genotype, capture efficiency was 60.9 ± 3.1% of input across the samples, with no significant difference in capture efficiency between the saline control and groups receiving long-term nicotine treatment (t test, p = 0.06, t = −2.108, df = 10), indicating that the immobilized mAb295 performed equally well across the samples regardless of differences in nAChR abundance. As expected based on epibatidine binding levels in homogenates and previously published studies, we observed a significant effect of nicotine use on the amount of β2*-nAChR present in each sample (Fig. 1B); however, no significant effect of mood disorder diagnosis or nicotine × mood disorder interaction was identified (two-way ANOVA: main effect, nicotine use: F(1,17) = 18.002, p < 0.0009; main effect, mood disorder: F(1,17) = 0.236, p = 0.634; interaction: F(1,17) = 0.06, p = 0.813).
Transgenic mice lacking expression of α4 and β2 nAChR subunits have been used previously to establish a gene dose-dependent linear relationship between α4β2*-nAChRs and their putative associated proteins from whole mouse brain samples using an immunocapture protocol with iTRAQ-labeled peptides for quantitative proteomics (McClure-Begley et al., 2013). Such an approach used nAChR subunit heterozygous and knock-out mice to vary the abundance of nAChRs reliably and to examine relative changes in associated proteins with linear regressions to enhance confidence in identifying interacting proteins. In the human tissue samples used in this study, we observed significant variation in the total nAChR content of each sample, partially due to inherent variability, but also because of the impact of nicotine use. The measured abundance of the β2 nAChR subunit (target of the IP) as well as the abundance of α4 nAChR subunits measured by LFQ proteomics correlate extremely well with the calculated nAChR content of each sample obtained by the depletion of specific [125I]-epibatidine binding from tissue extracts (femtomoles vs β2-LFQ: r 2 = 0.89; femtomoles vs α4-LFQ: r 2 = 0.96; Fig. 2A,B).
We identified 118 putative interacting proteins from the abundance estimates provided by proteomic analysis of mAb295 immunopurified β2*-nAChRs from human temporal cortex. We excluded obvious contaminants (epidermal keratins, hornerin, hemoglobins, histones, cytochromes) and any protein with fewer than two unique peptides identified at a false discovery rate of 5%. With these criteria applied, we identified and quantified 78 proteins that coimmunoprecipitated with β2*-nAChRs from human postmortem temporal cortex (Table 1). Of these, 42 proteins were identified as nAChR-interacting proteins in mouse cortex, and an additional 14 were isoforms of identified mouse cortex interacting proteins. An overlapping set of 36 proteins were identified in whole mouse brain as nAChR-interacting proteins, with an additional 17 isoforms of previously identified interacting proteins (McClure-Begley et al., 2013). Further, 33 of 78 proteins identified showed statistically significant correlations of protein abundance with β2-nAChR subunit abundance (Table 2), indicating a likely specific protein–protein interaction or physical proximity at the time of receptor extraction (Table 3). Functional annotation clustering of these 33 proteins using the DAVID Bioinformatics Database analysis (Huang da et al., 2009) showed a significant enrichment for synaptic pathways and neuronal projection (enrichment score, 5.47). Therefore, stable receptor interactions that survive tissue solubilization are likely assembled β2*-nAChRs, and a subset of these nAChRs may be resident in the plasma membrane of neuronal synapses. The identification of a 14-3-3 isoform previously shown to mediate α3*nAChR trafficking to postsynaptic structures (Rosenberg et al., 2008) indicates that some β2* interactions are mediated on both sides of CNS synapses by this family of adaptor proteins. Based on the current dataset, it is possible that some of the captured β2*-nAChRs reside in the postsynaptic plasma membrane, but we failed to identify additional known postsynaptic density components or associated proteins with functional annotations placing them in that neuronal subcellular compartment, so additional experiments will be required to address variations in postsynaptic versus presynaptic β2*-nAChR-interacting proteins.
Immunoprecipitation of β2*-nAChRs from mouse cortex and identification of specific interacting proteins by tandem mass spectrometry
As for the human samples, β2 nAChR binding sites in mouse cortex samples were immunoprecipitated by mAb295-Dynabeads and subsequently quantitated by [I125]-epibatidine binding. Binding site capture normalized to total tissue weight demonstrated significant effects of nicotine treatment and genotype (Fig. 3), but not a significant interaction (two-way ANOVA; main effect, nicotine, F(1,15) = 29.05, p < 0.0009; main effect, genotype: F(3,15) = 59.314, p < 0.0009; interaction, F(3,15) = 3.846, p = 0.057).
Proteins identified and quantified from mouse cortex were analyzed, as has been described by McClure-Begley et al. (2013), using identical criteria for inclusion and exclusion as the human tissue samples. Each mouse sample was pooled from three mice (two pools per condition × eight conditions [nicotine-treated and untreated mice; wild-type, heterozygous, double-heterozygous, and knock-out mice]), and we mixed the mouse tissue such that each sample contained a total of ∼0.62 g wet weight. As expected, nAChR abundance calculated from the depletion of specific [125I]-epibatidine binding sites (as femtomoles/milligram protein) correlated very well with the iTRAQ abundance ratios of α4- and β2-nAChR subunits (r 2 = 0.77 and 0.70, respectively).
Identified proteins were considered to be significant interactors according to the overall correlation of the log2-transformed ratios of mean iTRAQ reporter ions with the control sample for each set (in these samples, the control is tissue from wild-type animals; Fig. 3, WT treated with saline). We observed a total of 184 coimmunopurified proteins from the mouse cortex samples, and, after the elimination of obvious contaminants and proteins identified in only one biological replicate, there were 94 proteins in common between the human and mouse experiments (Table 4). Of the 94 proteins identified in both biological replicates (Table 4), 5 were significantly correlated with β2-nAChR subunit abundance; statistical power was reduced compared with the human samples because individuals in the mouse study were pooled for measurement rather than used as biological replicates (Table 5). In addition, the comparative mouse nAChR-associated proteomes were collected using a peptide-labeling iTRAQ strategy, which has been shown previously to provide less comprehensive sequence coverage and fewer quantitative differences than LFQ (Latosinska et al., 2015). Of the 94 proteins identified in the nAChR interactome in mouse cortex, only 19 were not identified in the human cortex samples directly, or as a close isoform (Table 4). Further comparison to proteins associated with the β2-nAChR subunit in whole mouse brain (McClure-Begley et al., 2013) identified 37 proteins that were also immunoprecipitated from mouse and human cortex (Table 6). Comparing these proteins to those present in a common repository of contaminant proteins (Mellacheruvu et al., 2013), 11 proteins were identified <10 times in control studies, 7 were identified in <100 studies, 16 were identified in >100 control samples, and 3 were not present in the database (Table 6).
Effects of nicotine exposure and mood disorder on nAChR-associated proteins isolated from human temporal cortex
To identify changes in the β2*-nAChR-associated proteome due to nicotine exposure status or mood disorder, we compared the abundance ratios for every protein identified in each sample and performed GLM analysis to identify significant changes. Of the 33 putative interacting proteins identified following immunoprecipitation of β2*-nAChRs from human cortex (Table 3), 14 proteins differed in interaction as a result of tobacco use, as measured by a shift in the ratio of the interacting protein to nAChR number between nonexposed and tobacco-exposed individuals, with two additional protein (14-3-3ζ and ATP1A2) that showed a significant nicotine × mood disorder interaction (Fig. 4A, Table 7). None of the interactions were significantly altered by mood disorder alone (Table 7). Several of the changes in associated protein abundance appear to be affected more significantly by the combination of nicotine use and mood disorder (Fig. 4, open triangles) despite the nonsignificant interaction term calculated from the dataset as a whole. Normalizing protein abundance to β2 nAChR content within each sample identified an independent set of five proteins that were significantly altered by nicotine use (Fig. 4B, Table 8); again, we did not identify significant differences by mood disorder diagnosis.
Effects of nicotine exposure on nAChR-associated proteins isolated from mouse cortex
In the mouse iTRAQ data, GLM analysis did not identify a significant effect of nicotine on any of the nAChR-associated proteins with expression indices significantly correlated with β2-nAChR abundance. Although nicotine exposure did increase the abundance of the α4- and β2-nAChR subunits, the increase was only marginally significant (Chrna4: F(1,13) = 4.93, p = 0.06; Chrnb2: F(1,13) = 4.67, p = 0.07), likely due to the small number of pools included in the analysis. In addition, the overall n size of the experiment was reduced because the β2-nAChR subunit abundance from the control groups is used as the denominator in their ratiometric quantitation, and is not included as a correlate to avoid bias due to artificially small errors on the calculated intercept. This ratiometric quantitation also prevents normalization to β2-nAChR subunit abundance within each sample, so a depiction of change relative to “per unit” β2 is not possible. Analysis of genotype effects on iTRAQ-labeled nAChR peptides with GLM showed significance for all three subunits identified, as follows: Chrna4: F(3,13) = 16.43, p = 0.001; Chrna5: F(3,13) = 8.87, p = 0.009; Chrnb2: F(3,13) = 12.78, p = 0.003.
Discussion
In this study, we examined the effects of tobacco use and mood disorder on the expression of α4β2*-nAChRs and their associated proteins in postmortem human cortex. We also used experiments in mice to determine the effects of nicotine alone on the associated proteome. Repeated or prolonged nicotine exposure increases α4β2*-nAChR number in heterologous expression systems, animal models, and human tissue (for review, see Melroy-Greif et al., 2016). In addition, the changes in α4β2*-nAChR content after long-term exposure to nicotine are brain region and cell type specific, with cortical regions among those most reproducibly affected in both mouse and human (Marks et al., 1992, 2011; Nashmi et al., 2007). Further, a substantial body of evidence points to a cholinergic contribution to the etiology and expression of affective disorders (for review, see Mineur and Picciotto, 2010), and that β2*-nAChR occupancy by endogenous ACh is significantly increased by current or past episodes of major depressive or bipolar disorder (Saricicek et al., 2012; Hannestad et al., 2013). Studies in mouse models indicate that even acute dysregulation of cholinergic neurotransmission can elicit changes in affective states, and that long-term manipulation of acetylcholine levels leads to behavioral responses that are commonly indicative of mood disorders, specifically depression- and anxiety-like behaviors (Mineur et al., 2013). Identifying the interactions of α4β2*-nAChRs with other synaptic proteins and functional regulators could therefore identify the novel elements of neuronal plasticity underlying the acquisition of nicotine tolerance and dependence, as well as improve understanding of how cholinergic regulation of mood is accomplished in the CNS.
For this study we measured cotinine levels in the brain tissue lysate (containing all of the cytosolic and interstitial fluid components) to establish nicotine/tobacco use status for each sample. After organizing the samples according to nicotine use and mood disorder status, we observed that the levels of α4β2*-nAChRs measured by radioligand binding with [125I]-epibatidine were consistent with significant receptor upregulation due to nicotine use. The levels of α4- and β2-nAChR subunits measured by label-free quantitative proteomics were consistent with the changes observed with ligand binding. These results indicated that, when performed appropriately, quantitative immunopurification and proteomics is a valid method for evaluating changes in nAChR expression.
A different form of quantitative proteomics (iTRAQ) was used to assess changes in α4β2*-nAChRs and associated proteins from transgenic mouse tissue. iTRAQ uses amine-reactive isobaric labels to facilitate the ratiometric quantitation of peptides from multiplexed samples in a complex mixture (Wiese et al., 2007). Multiplexed iTRAQ proteomics has been used previously with nAChR subunit transgenic mice to validate the approach for the untargeted identification of putative nAChR-associated proteins (McClure-Begley et al., 2013). In this project, we used iTRAQ to identify interacting proteins that are altered in their association with β2*-nAChR due to long-term nicotine treatment. It was our goal to test the hypothesis that a core set of α4β2* interacting partners would be conserved in the cortex from mouse to human, and that the identification of these complexes could guide future studies geared toward understanding nAChR regulation and function.
We took advantage of the inherent variability in β2*-nAChR abundance across samples (due to either individual differences and/or nicotine use status in the case of the human samples, and genetic control of nAChR abundance and nicotine exposure in the mouse samples) to match linear relationships of protein abundance indices with the measured abundance of β2-nAChRs, providing confidence that the proteins coimmunopurifying with the β2*-nAChRs interact in a meaningful context. Abundance correlations would be expected to be random and nonsignificantly correlated if interactions were nonspecific. We observed a high level of overlap between human and mouse cortical samples of signaling complexes coimmunoprecipitated with β2*-nAChRs.
The functional roles of nAChRs as neuromodulators on presynaptic terminals of many different neuronal types have been described in numerous previous studies (for review, see Barik and Wonnacott, 2009). Evidence in the current study from bottom-up quantitative proteomic analysis of α4β2*-nAChR-associated proteins identified in human and mouse cortical tissue not only supports a presynaptic role for these receptors, but assigns additional proteins with known roles in neurotransmitter release and neuronal excitability that begin to describe the composition of specific signaling complexes assembled with nAChRs. Specifically, components of the presynaptic active zone and the SNARE complex (syntaxin1B, syntaxin1A, syntaxin binding protein-1, SNAP25, and synaptotagmin) were identified with significant expression correlations with β2-nAChR subunits in the human temporal cortex samples; the same core set of proteins were also identified in the mouse cortex iTRAQ analysis. Smaller sample size due to the pooling of individual mice for analyses limited statistically significant correlation to synaptophysin; however, identification across two species with differing exposure modalities increases confidence that interaction of the set of nAChR-associated proteins identified from human brain are altered as a result of nicotine exposure.
Of the 17 proteins that were identified previously with iTRAQ using stringent criteria for dose relationships in mice with varying nAChR genotype in whole mouse brain (McClure-Begley et al., 2013), 8 were identified in both mouse and human cortex in the current study (α4- and β2-nAChR subunits, CaMKIIα, spectrin α and β, myosin 10, tubulin β3, and GFAP; Tables 1, 2). Syntaxin 1B, synaptotagmin 1, syntaxin binding protein 1, and several 14-3-3 isoforms were also identified in all three datasets. These proteins, although clearly identified in the nAChR-associated proteome from whole mouse brain, did not meet the stringent criteria for variance with gene dose. One possibility is that each of these proteins is limiting in the compartment (synaptic or ER) in which it interacts with nAChRs, and that the interaction is functional and transient (during vesicle recruitment or release, or during nAChR assembly), so it does not scale with nAChR number. Of those proteins that were not identified in the current study, several proteins were represented by different isoforms or by other proteins with related functions. For example, while the current study did not identify CaMKIIγ, the β-isoform was identified in both mouse and human cortical samples, perhaps because the γ-isoform is more highly expressed in other brain regions, and the β-isoform in enriched in cortex.
Many subtypes of nAChRs have been identified as presynaptic heteroreceptors that can directly influence neurotransmitter release, and the identification of SNARE complex proteins as significant nAChR interactors in both the mouse and human tissue samples suggests a functional interaction between α4β2*-nAChRs and neurotransmitter release machinery. Minimally, this association indicates that some α4β2*-nAChRs are extremely close to sites on presynaptic terminals where neurotransmitter release is actively taking place. Such close physical association suggests that nAChRs may exert some modulatory influence on neurotransmitter release independently of action potential arrival or frequency, and highlights their potential role as a “fine-tuning” influence on activity at highly specialized synapses. It is also possible that the significant interaction of synaptic vesicle components with β2*-nAChRs indicates the presence of α4β2*-nAChR-mediated neurotransmitter release process from perisynaptic sites, similar to what has been proposed for α7 and α3* nAChRs in the chick ciliary ganglion (Coggan et al., 2005). Interestingly, the degree of interaction between syntaxin1A and β2*-nAChRs was reduced in human cortical tissue from nicotine users, raising the possibility that long-term nicotine use promotes a functional reorganization of presynaptic nAChR signaling complexes. The family of 14-3-3 isoforms have multiple functions (for review, see Yaffe, 2002), including acting as chaperones and adaptors that mediate complex protein–protein interactions, and have been shown to facilitate nAChR trafficking and assembly of α3* nAChRs at postsynaptic sites in the ciliary ganglion (Rosenberg et al., 2008). The current finding showing that 14-3-3ζ interacts with β2*-nAChRs suggests that adapter protein interactions are critical determinants of trafficking and functional regulation of nAChRs on both presynaptic and postsynaptic plasma membranes and that these interactions are likely to be context specific. Nicotine use may also induce restructuring of postsynaptic elements of cholinergic synapses, but in this study we did not reliably identify components of larger postsynaptic complexes associated with β2*-nAChRs. It is possible that in the cortex of both mice and humans, the postsynaptic β2*-nAChR interactome is more dynamic and less likely to survive the conditions used for isolating nAChRs in this study; additional studies using methods for enriching and isolating nAChRs from specific subcellular neuronal compartments is necessary to address this question.
The proteoglycan agrin, which contributes to the orchestration of nAChR clustering at the neuromuscular junction (NMJ; Wallace et al., 1991), has also been identified as an inhibitor of α3-containing Na/K ATPases (Hilgenberg et al., 2006). We identified α3 Na/K-ATPase as a β2*-nAChR-interacting protein from both human and mouse cortical tissues. Its abundance was upregulated along with the β2*-nAChR subunit in the human samples with nicotine treatment, as evidenced by the similar ratio of α3 Na/K-ATPase/β2-nAChR levels across treatment groups. This parallel increase suggests that aspects of presynaptic neuronal excitability are changed by long-term nicotine treatment beyond the alterations in excitation due to prolonged nAChR desensitization or other mechanisms of agonist-induced functional changes. An additional link between altered Na/K-ATPase activity and mood disorder is found in the Myshkin mouse model, which harbors a loss-of-function mutation in the α3Na/K-ATPase and exhibits behavioral hallmarks of the manic phase of bipolar disorder that are responsive to traditional pharmacotherapies (valproate, lithium) as well as inhibitors of downstream signaling linked to α3Na/K-ATPase (Kirshenbaum et al., 2011). Postsynaptic nAChRs at the NMJ copurify with the α2Na/K-ATPase subunit (Heiny et al., 2010; Talsma et al., 2014), and we observe significant correlation here of α2-Na/K ATPase abundance with β2*-nAChRs. Genetic risk variants for mood disorders have been found in the Na/K ATPase α2- and α3-loci (Talsma et al., 2014), highlighting the potential functional interactions between nAChRs and Na/K-ATPases in the CNS that may also impact lifetime risk for developing mood disorders. The fact that we observe α4β2*-nAChRs associated with Na/K-ATPase subunits as well as SNARE complex machinery and other mediators of synaptic activity leads us to consider that presynaptic nAChRs may exist as part of a core activity complex with substantial homology to what has been shown to exist postsynaptically at the NMJ (Heiny et al., 2010; Talsma et al., 2014). Further experiments and fine-detailed structural analysis of presynaptic nAChR scaffolds in the CNS will, it is hoped, add to this concept of conserved cholinergic regulation via novel functions of proteins previously unknown as functional regulators of nAChRs.
Another specific β2-interactor of interest for additional study is the protein TRIM21, an E3 ubiquitin ligase with roles in cell cycle progression and immune signaling (Wada and Kamitani, 2006). Muscle-type nAChR turnover is at least partially dependent on the targeted autophagy of ubiquitinated receptors (Khan et al., 2014), and the role of the ubiquitin–proteasome system in the trafficking and stability of α3* nAChRs has been shown in both cell lines (Rezvani et al., 2010) and rodent neurons (Teng et al., 2015). Different E3 ubiquitin ligases were determined to be critical mediators of the selective degradation of ubiquitinylated muscle and α3* nAChR subtypes in the previous studies (TRIM63 and CHIP, respectively); our identification of TRIM21 as a β2*-nAChR-associated E3 ubiquitin ligase in human cortical neurons indicates that future studies are warranted to determine whether it is a specific regulator of α4β2*-nAChR turnover.
A recent publication presents a hypothesis that nAChRs, including β2*-nAChR, may possess a structural motif that promotes functional interaction with Gα subunits of the heterotrimeric G-protein complex (Kabbani et al., 2013). Earlier studies using proteomic methods to interrogate β2-nAChR subunit interactions identified Gα subunits as a reliable associated protein in mouse brain tissue lysate (Kabbani et al., 2007). Here we show that β2*-nAChRs from human cortical tissue also associate with Gα subunits (particularly, the Gαo subtype), and that nicotine use significantly affects their abundance index across sample groups. Whether this change in Gαo interaction with β2*-nAChRs is functionally relevant to the physiological response to short-term or long-term nicotine use remains to be elucidated, but our results supporting the earlier discovery that a β2–Gα interaction is taking place provides considerable support for additional studies aimed at dissecting the roles of metabotropic signaling via nAChR ion channels.
Of the identified proteins, only two, 14-3-3ζ and ATP1A2, showed a significant nicotine × mood disorder interaction. As noted above, the Na/K ATPase is critically involved in neuronal energetics and at the presynaptic terminal and variants in the ATPase have been associated with mood disorders (Talsma et al., 2014). In contrast, 14-3-3 proteins are involved in chaperoning nAChRs from the endoplasmic reticulum to the cell surface, and this chaperone function is regulated by protein phosphorylation (Jeanclos et al., 2001). The significant nicotine × mood disorder interaction in the association of these proteins with β2*-nAChRs may therefore suggest that nAChR trafficking is altered differentially by nicotine exposure in individuals with mood disorders. This is an intriguing hypothesis that can be evaluated in future studies.
Our data have identified several proteins involved in the regulation of neuronal excitation and neurotransmitter exocytosis that are proximally associated with nAChRs in humans and mice, and are significantly affected by physiological conditions known to impact the expression of nAChRs themselves. While mood disorders have been demonstrated to involve significant contributions from cholinergic components, the impacts of smoking status on nAChR expression as well as the expression of interacting proteins appears to be a much more powerful perturbation. It is possible that different mood disorders (i.e., bipolar disorder vs major depression vs general anxiety disorder) are also representative of distinct entities with respect to cholinergic dysfunction that would benefit from additional exploration of discrete changes in the nAChR-associated proteome. These results indicate that the regulation of nAChRs and other membrane signaling complexes in skeletal muscle share core regulatory and functional interactions with β2*-neuronal nAChRs, albeit with different proteins that serve functionally similar roles. This study shows that proteomic methods are a powerful tool to identify proteins associated with nAChRs. Proteomic strategies were different for the human and mouse β2-nAChR-interacting protein identifications, and, despite obvious technical distinctions, strengths, and weaknesses, remarkably conserved protein complexes were found to be associated with assembled receptors. Additional proteomic studies aimed at a finer dissection of cell type-specific nAChR interactions may identify more targets that regulate and are regulated by nAChR activity to provide additional mechanistic insight into the contribution of cholinergic dysfunction in neurological disorders.
Acknowledgments
Acknowledgments: We thank Dr. Gustavo Turecki of McGill University for providing the brain samples, the W.M. Keck Biotechnology Resource Facility at the Yale School of Medicine for the SRM analysis of tissue cotinine content, and Megan Miller for critical feedback on the manuscript.
Footnotes
The authors declare no competing financial interests.
This research was funded by Grants MH-077681 and DA-14241 (to M.R.P.), K01-MH-092681 (to I.E.), DA-018343 (National Institute on Drug Abuse Proteomics Center at Yale University), and UL1-RR-024139 (Yale Clinical and Translational Science Award) from the National Institutes of Health (NIH); and by the State of Connecticut, Department of Mental Health and Addiction Services. T.D.M.-B. was supported by NIH Grant T32-MH-014276. M.J.M. and S.R.G. were supported by NIH Grants DA-003194 and DA-015663. J.M.L. was supported by NIH Grant NS-11323.
This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International, which permits unrestricted use, distribution and reproduction in any medium provided that the original work is properly attributed.