Structural and Functional Consequences of Connexin 36 (Cx36) Interaction with Calmodulin

Introduction

Electrical and chemical synapses can form mixed synapses conferring synchronous activity among coupled neurons in the vertebrate nervous system (Hamzei-Sichani et al., 2012; Nagy et al., 2013; Rubio and Nagy, 2015). This enables synchronous oscillations that temporally correlate firing patterns in otherwise heterogeneous networks (Galarreta and Hestrin, 1999; Deans et al., 2001; Bennett and Zukin, 2004; Vervaeke et al., 2010; Haas et al., 2011; Postma et al., 2011). In many regions of the nervous system, synchronous rhythmic oscillations serve as a mechanism to entrain the activity of distributed neuronal networks, and promote processes related to memory formation and retrieval, sensory perception, motor control, attention and plasticity (Deans et al., 2001; Hormuzdi et al., 2004; Van Der Giessen et al., 2008; Pereda, 2014).

Neurons express connexin36 (Cx36) and form gap junctions, which are the molecular substrate of electrical synapses. Functional plasticity of neuronal gap junctions in goldfish Mauthner cells (Pereda et al., 1998, 2003, 2004; Flores et al., 2010), and in both teleost and rodent retina (Li et al., 2009; Cachope and Pereda, 2012; Kothmann et al., 2012) involves interaction of Cx36 with the calcium/calmodulin-dependent kinase II (CaMKII). In this process, CaMKII binds Cx36 in a manner analogous to binding of this kinase to the NR2B subunit of the NMDA receptor and phosphorylates the connexin protein both in situ and in vitro (Alev et al., 2008). Blocking of Cx36–CaMKII binding and phosphorylation by decoy inhibition, and deletion of CaMKII binding sites, as well as inhibition of phosphorylation sites through site-directed mutagenesis and peptide inhibition, leads to a significant loss of the “run up” activity, the term coined by us to describe the plasticity of the process (Zoidl et al., 2002; Del Corsso et al., 2012).

Interaction of Cx36 with CaMKII endows electrical synapses with functional plasticity but is in competition with binding of calmodulin (CaM), another multifunctional Ca²⁺ signaling protein (Burgoyne, 2007), to a region overlapping with the CaMKII binding motif (Burr et al., 2005). At chemical synapses both CaM and the NMDA receptor subunit NR2B synergize to lock CaMKII in an active conformation conferring calcium-dependent plasticity (Bayer et al., 2001; Lisman et al., 2002; Bhattacharyya et al., 2016). A similar mechanism involving Cx36, CaM and CaMKII appears unlikely due to the overlapping binding site. Addressing this conflict, we have proposed a model predicting sequential, competitive processing of interaction between CaM, CaMKII, and Cx36 (Alev et al., 2008). To resolve this conflict, the neuroblastoma cell line Neuro2a was used to interrogate how calcium signaling converges on Cx36 through CaM binding. Using Förster resonance energy transfer analysis (FRET), we demonstrate that formation of the Cx36–CaM binding complex is dependent on the capacity of CaM to bind calcium and strengthened when intracellular calcium concentrations rise. The formation of the Cx36–CaM complex precedes arrival of the Cx36 protein at the gap junction plaque (GJP) or is established after removal of Cx36 from the GJP. Functional dye transfer/dye recovery assays in cell pairs under physiological conditions showed that Cx36–CaM binding endows GJPs with the capacity to exchange a fluorescent dye between adjacent cells. Pharmacological blocking of CaM or mutation of the critical amino acid W277 within the core CaM binding motif of Cx36 significantly diminished this property. This observation led us to propose a revised model of calcium signaling where the localization of Cx36 in the cell affects the participants (CaM and/or CaMKII) and consequently, the outcomes arising from these interactions. This model is additionally supported by detailed structural and kinetic view of the Cx36–CaM interaction using a combination of the NMR and calorimetric approaches. In summary, this study demonstrates that priming of Cx36 by CaM represents a critical stage in the creation of plasticity at electrical synapses.

Materials and Methods

Plasmid Constructs and Site-Directed Mutagenesis

Expression vectors contained full-length rat Cx36 [NM_019281, amino acids (aa) 1–321], CaM (NM_031969, aa 1–149), and CaMKIIa (NM_012920.1, aa 1–478) in variants of pEGFP-N1 expression vectors (Clontech Laboratories Inc., Mountain View, CA, USA) (Zoidl et al., 2002; Alev et al., 2008). Full-length in-frame cloning of CaM, the calcium insensitive mutant CaM_E1234Q (Maune et al., 1992; VanScyoc et al., 2002), as well as genetically encoded organelle markers (pER-DsRed2, pGolgi/DsRed2, pMito-DsRed2) (Fliegel et al., 1989; Rizzuto et al., 1989; Watzele and Berger, 1990) was completed in two steps. The open reading frames were synthesized as gBlocks [Integrated DNA Technologies Inc. (IDT), Coralville, IA, USA] and cloned into the TA cloning vector pJet1.2 (Thermo Fisher Inc., Mississauga, ON, Canada). Then the coding regions were isolated by restriction digest and cloned in-frame into pECFP-N1 or the pDsRed2-monomer expression vector. Mutant Cx36 plasmids (G276A, W277A, R278A) were generated using the Q5^TMSite-Directed Mutagenesis Kit according to the manufacturer’s protocol [New England Biolabs Inc. (NEB), Boston, MA, USA]. Oligonucleotides were designed using the NEBaseChanger tool (NEB) and synthesized by IDT. Table 1 summarizes the oligonucleotides utilized in this study (mutations depicted in bold and underlined). All plasmid constructs used on this study were sequence verified (Eurofins MWG Operon LLC, Huntsville, AL, USA).

TABLE 1

TABLE 1. Oligonucleotides for Cx36 mutagenesis.

Cell Culture, Transient Transfection, and Western Blot

Neuroblastoma 2a (Neuro2a) cells (Olmsted et al., 1970) were cultivated in DMEM with 2 mM glutamine, 1% non-essential amino acids (NEAA), 1% penicillin and streptomycin (PS), and 10% fetal bovine serum (FBS) at 37°C in a humidified atmosphere with 5% CO₂. Cells were seeded in 24-well plates or glass-bottom dishes (MatTek Corporation, Ashland, MA, USA) and transfected with 200 ng (single transfection) or 400 ng (double transfection) endotoxin-free plasmid DNA, using the Effectene transfection protocol (Qiagen Inc., Valencia, CA, USA). For western blot, whole cell protein lysates were prepared 48 h after transfection. The 20 μg protein was separated by 10% SDS-PAGE, transferred to 0.2 μm Midi format nitrocellulose membrane and processed using the iBind^TM Western System (Bio-Rad Inc., Mississauga, ON, Canada). Primary antibodies were diluted 1:1,000 [mouse anti-GFP, Roche; rabbit anti-GFP (FL), Santa Cruz Biotechnologies, Dallas, TX, USA] and 1:20,000 (mouse anti-β-actin; Sigma-Aldrich Chemie GmbH, Munich, Germany). The secondary antibodies (LI-COR Biosciences, St. Lincoln, NE, USA) were diluted 1:20,000 (donkey anti-rabbit IRDye680LT) or 1:20,000 (goat anti-mouse IRDye800CW). Signals were detected using the Odyssey^® CLx Infrared Imaging System (LI-COR Biosciences).

Confocal Microscopy

Transfected cells were fixed with 4% paraformaldehyde for 20 min at RT, washed with PBS, and mounted with ProLong^® Antifade Mountant (Thermo Fisher Inc., Mississauga, ON, Canada) for imaging. Samples were visualized using a Zeiss LSM 700 confocal microscope with a Plan-Apochromat 63x/1.4 Oil DIC M27 objective and the ZEN 2010 program to control all hardware parameters. Images were collected by line averaging (4×) at high resolution (2048 × 2048 pixels) using single planes or z-stacks. The GJ frequency was quantified as described previously (Qu et al., 2009), by counting the number of EGFP-expressing cell pairs against the number of cell pairs having GJPs in images from randomly selected, non-overlapping visual fields collected from experimental replicates (n > 3). The GJ area was quantified using ImageJ. Here, the ImageJ free hand tool was used to draw the contours of GJs between cell pairs, followed by quantification of the GJ area using the measure tool. Mander’s overlap coefficients were calculated using the co-localization analysis if the ZEN 2010 program (Manders et al., 1993). Images were exported and further processed using ImageJ or the Imaris Scientific 3D/4D Image Processing & Analysis Software (Bitplane USA, Concord, MA, USA). Images were combined using Adobe Photoshop CS6 for presentation.

FRET Analysis

The acceptor bleach protocol was established by testing different combinations of fluorescently tagged Cx36 and CaM proteins using a Zeiss LSM 700 confocal microscope (data not shown). In the experiments shown, Neuro2a cells were transiently transfected with DsRed-monomer and CFP tagged expression vectors and processed as described above. Baseline readings were measured prior to the acceptor bleach protocol. Here, DsRed-tagged proteins were bleached using the 555 nm laser line (set to 100% intensity) and the intensity change of CFP tagged proteins was recorded using the 405 nm laser line. The experiment was completed when the bleaching of the acceptor channel reached 10% of the initial intensity. FRET efficiency was then calculated by using the FRET efficiency formula: FRET_eff = (D_post -D_pre)/D_post, where D_post is the average intensity after the bleach, and D_pre is the average intensity before the bleach (Okamoto et al., 2004) after subtracting the values of the background (noise). The FRET distance (FRET_d) between the interacting proteins was calculated using the formula: R = R_o ((1/E) - 1)^1/6, where R_o is the reference distance between the two fluorescent tags (DsRed and ECFP: 5.09) (Erickson et al., 2003), and E is the FRET efficiency. The mean values of the intensity change and distance were calculated and formulated into graphical representations. Error bars denote the standard error of the mean (SEM) of each protein–protein interaction.

Ethidium Bromide Uptake and Recovery after Photobleaching Assay

This assay is a variation of the ethidium bromide uptake assay previously reported (Prochnow et al., 2009; Kurtenbach et al., 2013). Transfected Neuro2a cells grown for 48 h in 3.5 cm MatTek cell culture dishes were incubated with 10 μM of ethidium bromide in complete growth medium for 10 min at 37°C and 5% CO_2. MatTek chambers were placed in a live-cell imaging chamber and imaged at 37°C using a Zeiss 700 confocal microscope. Cell pairs expressing Cx36–EGFP were selected and a time-lapse baseline image was recorded (defined as 100%). Next, one cell was bleached using the 555 nm laser line with the intensity set to 100% laser power. The number of iterations in the bleaching protocol was terminated when 40–50% reduction in fluorescence was achieved. In this process, a volume of 36.2 μm × 36.2 μm × 1.6 μm was bleached using a non-confocal imaging mode. The recovery of ethidium bromide fluorescence in the center plane of this volume after bleaching was measured over time in confocal imaging mode. In all experiments, three regions of interest (ROI: R0, R1, R2) were selected for analysis. R0 was placed outside the cell pair to determine background the fluorescence. ROIs were placed inside the bleached cell close (R1) and distant (R2) to the GJP. The fluorescence recovery in percent was calculated using the formula:

$$\text{Fluorescence} \text{recovery} (\%) \text{=}\frac{({\text{F}}_{\text{3} \text{min} \text{post} \text{bleach}} \text{-} {\text{F}}_{\text{bleach}})}{({\text{F}}_{\text{3} \text{min} (\text{pre-bleach})} \text{-} {\text{F}}_{\text{pre} \text{-} \text{bleach}})}$$

Pharmacology

Prior to imaging transfected cells were treated with the Ca²⁺ chelator BAPTA-AM (Thermo Fisher) (Ueda and Shah, 1992), ionomycin (Sigma-Aldrich) (Liu and Hermann, 1978), carbenoxolone (CBX; Sigma-Aldrich) (Li et al., 2001) or by applying the CaM antagonist W-7 (Sigma-Aldrich) (Weiss et al., 1982). The final concentrations in complete growth medium were 24 μM of BAPTA-AM, 10 μM of W-7, 50 μM of CBX, for 10 min or 2 μM ionomycin for 5 min at 37°C and 5% CO₂. The final solvent concentration of ethanol in ionomycin treated cells was 0.002%.

Protein Expression and Purification

A pET15 based (EMD Millipore, Etobicoke, ON, Canada) plasmid for expression of hexahistidine-tagged human calmodulin (His₆-CaM) was a gift from Prof. Mitsu Ikura (University of Toronto). From a 1 L fermentation in Luria broth, milligram quantities of >95% pure His₆-CaM were obtained by a two-step Ni-NTA affinity chromatography (Qiagen Inc., Toronto, ON, Canada) and gel filtration chromatography (Sephadex-75 HiLoad 16/60 column, GE Life Sciences) procedure. Final buffer conditions were 5 mM Tris-HCl, 100 mM NaCl, 0.05% NaN₃, pH 7.5. The CaM–Cx36 complex pursued for structural determination was expressed as a 169 aa fusion protein consisting of His₆-tagged human CaM (5–150), a short linker (GASTAAGS), and residues 276–292 of rat Cx36. A gene fragment encoding this fusion protein was synthesized by DNA2.0 (Menlo Park, CA, USA) and supplied in a T5-based expression vector (pJ441-kan). ¹⁵N, ¹³C labeled His₆-CaM–Cx36 was obtained from a 1 L fermentation in M9 media containing 1 g (99%, ¹⁵N) ammonium chloride, 3 g (99%, ¹³C-glucose) and 0.5 g (98%, ¹³C–¹⁵N) algal extract (Cambridge Isotopes, Tewksbury, MA, USA). Final buffer conditions for NMR spectroscopy were 1.5 mM protein in 5 mM Tris-HCl, 100 mM NaCl, 5 mM CaCl₂, 0.05% NaN₃, pH 7.5 with 10% (v/v) ²H₂O.

NMR Spectroscopy

Spectra were acquired at 298 K with a 700 MHz Bruker Avance spectrometer equipped with a TXI cold probe. Backbone and side chain chemical shift assignments were obtained using a suite of standard heteronuclear, triple resonance experiments (HNCACB, CBCACONH, HNCO, HNCACO, HNCA, CCH-TOCSY, HCCH-TOCSY) with non-uniform Poisson gap sparse sampling at 10–15%. Distance restraint data were obtained from ¹⁵N-edited and ¹³C-edited 3D-NOESY spectra (120 ms mixing time). All datasets were processed with NMRPipe (Delaglio et al., 1995) and analyzed with CcpNmr Analysis (Skinner et al., 2015).

Structure Determination

CcpNmr Analysis was used to integrate NOEs and calculate distance restraints based on an r^-6 function spanning 1.8–5.5 Å with a reference distance set at 3.2 Å. These restrains were used as input to CYANA (Buchner and Güntert, 2015) along with Ca²⁺ coordination distance restraints, hydrogen bonding distance restraints and torsion angle restraints predicted by the DANGLE method embedded in CCPN Analysis. Allowed backbone torsion angles and side chain rotamers were enforced by the ramaaco and rotameraco methods embedded in the CYANA suite. The final representative ensemble of 20 structures based on lowest CYANA target function was selected from an initial calculation of 400 structures. As this ensemble demonstrated good geometry and bonding, no additional refinement was pursued.

Databases

The structural ensemble was deposited in the Protein Data Bank under the accession number 2N6A.

Isothermal Titration Calorimetry (ITC)

Isothermal titration calorimetry (ITC) measurements were performed on MicroCal VP-ITC calorimeter (MicroCal Inc., Northampton, MA, USA). The recombinant H₆-CaM (affinity purified as described above) was titrated with Cx36 derived peptides reflecting the wild type (GSGWRKIKLAVRGAQAKRKSVYEIR; CanPeptide Inc., Montréal, QC, Canada) or the W227A mutant (GSGARKIKLAVRGAQAKRKSVY; BioBasic, Markham, ON, Canada). An optimal titration was achieved with 25 μM CaM in the reaction cell and 250 μM Cx36 peptide in the syringe, each in a buffer containing 50 mM NaCl, 5 mM BisTris pH 7.0, 5 mM CaCl₂. After an initial injection of 2 μL, the bulk of the titration consisted of 34 successive 8 μL injections with an equilibration delay of 300 s. Heats of dilution were determined by titrating the same peptide solutions into buffer alone. Titration profiles corrected for the heat of dilution were fitted into one-binding-site model using MicroCal Origin v5.0. All values reported were calculated based on three individual experiments.

Statistics

Statistical analysis and data presentation was performed using the IBM-SPSS and R-software packages. Results shown derive from experimental replicates with n ≥ 3. All data were analyzed for data distribution and subjected to Mann–Whitney U tests for independent samples or a paired t-test, when appropriate.

Results

Cx36–CaM at Gap Junctions in Neuro2A Cells

Previous work (Burr et al., 2005) and the sequence alignment of Cx36 orthologs representing major vertebrate phyla ranging from primitive fish to humans suggested the existence of a highly conserved potential CaM binding motif in the carboxy-terminus of Cx36 resembling a variation of the 1–8–14 subclass. This motif overlapped with a binding site for the calcium-dependent calmodulin kinase II (Figure 1A) (Alev et al., 2008).

FIGURE 1

FIGURE 1. Expression, localization, and gap junction plaque (GJP) formation of Cx36 in Neuro2a cells. (A) Sequence alignment of Cx36 orthologs from major vertebrate phyla showing a conserved CaM-binding domain overlapping with the CaMKII pseudosubstrate binding site. The 1–8–14 positions of the CaM-binding motif are highlighted. (B) Western blot analysis of transiently transfected Neuro2a cells expressing Cx36 wild-type, Cx36 mutants, CaM, and CaM_E1234Q. Expressed proteins were detected using an anti-EGFP antibody (in red), and an anti-beta-actin antibody (in green) as the control. M indicated the protein standard. (C) Rectangles indicate representative GJPs in Neuro2a cells double transfected with CaM (blue), Cx36 (red) wild-type, and mutant Cx36 proteins. Scale bar: 10 μm. (D,E) Quantification of gap junction size and frequency demonstrating that gap junction size is not affected by mutations of the core CaM binding region, but mutation W277A affects GJP frequency. Error bars box plot: maximum and minimum range with median; Mann–Whitney U test, ^∗∗∗p < 0.001.

The role of W277 in position 1 of the putative CaM-binding site was investigated after sequential alanine-scanning mutagenesis of positions 276 to 278 of the rat Cx36 protein. Western blot analysis confirmed exogenous expression of Cx36 wild-type and Cx36 mutants (G276A, W277A, and R278A) in single transfected Neuro2a cells 48 h post transfection (Figure 1B). A similar result was obtained with CaM and the Ca²⁺ insensitive mutant of CaM (CaM_E1234Q; Maune et al., 1992) generated for this study (Figure 1B). The localization of Cx36, Cx36 mutants, and CaM was investigated in double-transfected cells 48 h post transfection. The formation of GJPs was observed in cells co-expressing CaM with wild-type Cx36 and mutant proteins (Figure 1C). A quantification of the GJP size revealed no significant differences in the GJP sizes by the mutations introduced (Figure 1D) (values in μm²; WT: 1.72 ± 0.20, G276A: 1.22 ± 0.09, p = 0.05, W277A: 1.37 ± 0.11, p = 0.19, R278A: 1.13 ± 0.11, p = 0.01). The number of cells with GJP formation in transfected cells was significantly reduced for mutant W277A, but not when WT was compared to G276A and R278A (Figure 1E) (values in %; WT: 13.68 ± 2.52, G276A: 11.24 ± 2.03, p = 0.5, W277A: 3.63 ± 0.72, p = 1.9 × 10^-5, R278A: 8.69 ± 1.55, p = 0.12).

Mutation W277A Alters Subcellular Distribution and Exchange between Cell Pairs

Two types of experiments were performed to clarify the distribution of Cx36 and W277A in Neuro2A cells. When single transfected Neuro2A cells expressing Cx36–EGFP or Cx36–DsRed monomer were paired and analyzed 48 h after transfection, 3D imaging of cell pairs showed mixed GJPs (see inset in Figure 2A; Supplementary Figure S1), and limited presentation of Cx36 in mixed vesicles. Cell pairs expressing the W277A mutation showed mixed GJPs (see Supplementary Figure S2), and an extensive display of vesicular structures exchanged between cell pairs, suggesting that the mutation altered GJP stability (Figure 2A).

FIGURE 2

FIGURE 2. Colocalization of Cx36 with organelle markers indicating interaction occurring as early as at the endoplasmic reticulum. (A) Representative 3D images of Neuro2A cells expressing wild-type Cx36 or the W277A mutant. Note that cell pairs derived from single independent transfections tagging the proteins of interest with EGFP or DsRed monomer. The increased presentation of mixed red and turquoise vesicles in W277A transfected cell pairs suggests reduced GJP stability and increased protein turnover. Insets represent magnifications of the GJP region. (B) Neuro2a cells expressing wild-type Cx36–ECFP (in red) and DsRed monomer tagged organelle markers (endoplasmic reticulum, Golgi apparatus, and mitochondria, in light turquoise). The localization of Cx36 plaques is indicated with white arrows. (C) Mander’s overlap coefficients calculated with the ZEN 2010 program for Cx36 WT (n = 30) and W277A (n = 40) with the different organelle markers are shown with error bars indicating the maximum and minimum range. 3D images in (A,B) were created with Bitplane Imaris 7.6.4. Scale bars in (A,B): 10 μm. Error bars box plot: maximum and minimum range with median; Mann–Whitney U test, ^∗∗∗p < 0.01.

Next, the subcellular localization of Cx36 was visualized using Neuro2a cells double transfected with EGFP tagged Cx36 and genetically encoded cell organelle markers for endoplasmic reticulum (ER), Golgi, and mitochondria (Figure 2B, Supplementary Figure S3). Wild-type Cx36 was predominantly localized in vesicles and at GJPs. Outside of GJPs, both 2D images (data not shown) and 3D reconstructions of images from double transfected cells showed Cx36 in vesicles of variable size embedded in, or attached to ER and Golgi (Supplementary Figure S4). Co-localization analysis revealed a significant overlap of Cx36 with ER (0.695 ± 0.0413, n = 30, p = 9.837 × 10^-6) and Golgi organelles (0.651 ± 0.0355, n = 30, p = 4.209 × 10^-5). The Mander’s overlap coefficient for Cx36 and mitochondria was significantly lower (0.356 ± 0.0413, n = 30) demonstrating that exogenously expressed, tagged wild-type Cx36 is not redistributed into this organelle type (Figure 2C; Supplementary Figure S5). The distribution of Cx36 in cell pairs, together with the localization of tagged Cx36 in ER, Golgi and GJPs is consistent with the subcellular distribution of other connexin proteins. When Cx36–W277A was expressed, the Mander’s overlap coefficient was significantly higher than wild-type Cx36 in all organelles (ER: 0.814 ± 0.012, p = 0.009; Golgi: 0.811 ± 0.017, p = 1.6 × 10^-5; mitochondria: 0.812 ± 0.010, p = 2.22 × 10^-9). The distinct subcellular distribution of the W277A mutant suggests a mechanism that interferes with the transport, formation and/or stabilization of GJPs in transfected Neuro2a cells.

Cx36–CaM Binding Is Strengthened by Elevated Intracellular Calcium at Intracellular Vesicles

The interaction of Cx36–DsRed with CaM–ECFP in living cells was investigated using Fluorescence Resonance Energy Transfer (FRET) (Figure 3). The Cx36–ECFP and Cx36–DsRed assembly into mixed GJPs (Figure 2A) served as a positive FRET control. The FRET efficiency value (FRET_eff expressed in %) was 11.50 ± 0.93. The interaction of Cx36 with the Ca²⁺/calmodulin-dependent protein kinase II (CaMKII) has been shown previously (Alev et al., 2008; Del Corsso et al., 2012). Here, Cx36–ECFP and CaMKII–EYFP served as the second positive control demonstrating that intracellular calcium elevation by treatment with 2 mM extracellular calcium [Ca²⁺]_E and 2 μM ionomycin increased the FRET_eff significantly at GJPs (values in %: non-stimulated: 5.1 ± 0.4; stimulated: 8.1 ± 0.5, p = 0.0002), but not at vesicular structures (non-stimulated: 6.7 ± 0.9; stimulated: 8.1 ± 0.9, p = 0.33) (Figure 3A).

FIGURE 3

FIGURE 3. FRET analysis of Cx36 and CaM interaction. Neuro2a cells were transiently co-transfected, and FRET efficiencies (FRET_eff) determined 48 h post transfection. (A) FRET_eff of Neuro2a cells expressing active FRET control pairs. (B) DsRed and ECFP protein tags are interchangeable and with minimal impact on FRET_eff. (C,D) FRET_eff at GJPs and vesicles without and with 2 mM extracellular calcium [Ca²⁺]_E and 2 μM ionomycin in the presence of the pharmacological calmodulin blocker W-7 (10 μM), the calcium chelator BAPTA-AM (24 μM), or the Ca²⁺ insensitive CaM_E1234Q mutant. (E,F) FRET distance (FRET_dis) highlighting that calcium-dependent binding of Cx36 and CaM is significantly shortened at vesicles. Abbreviations: Iono, treatment with 2 mM [Ca²⁺]_E and 2 μM ionomycin. V, vesicles; GJP, gap junction plaque. Error bars box plot: maximum and minimum range with median; Mann–Whitney U test significance, ^∗∗∗p < 0.01, NS, not significant.

The FRET_eff for Cx36–DsRed and CaM–ECFP was comparable at GJPs and in cytoplasmic vesicles under physiological conditions without calcium stimulation (GJP: 3.23 ± 0.26; vesicles: 3.65 ± 0.54, p = 0.54) (Figure 3B). When the fluorescent protein tags (FRET pair I: Cx36–DsRed, CaM–ECFP; FRET pair II: Cx36–ECFP, CaM–DsRed) were switched (GJP: 3.35 ± 0.49, vesicles: 3.09 ± 0.40, p = 0.82) FRET_eff was not significantly different, suggesting that the protein tags had equal impact.

No significant changes to FRET_eff were observed at GJPs when intracellular calcium levels were raised (non-stimulated: 3.3 ± 0.5, stimulated: 2 mM [Ca²⁺]_E with 2 μM ionomycin; 2.8 ± 0.6, p = 0.33). Intracellular calcium reduction (24 μM BAPTA-AM; 2.8 ± 0.9, p = 0.42), or blocking of CaM with 10 μM (W-7; 2.7 ± 0.6, p = 0.44) had no impact on FRET_eff. Co-transfection with the calcium insensitive CaM_E1234Q mutant (non-stimulated: 1.1 ± 0.3; stimulated: 1.3 ± 0.4, p = 0.85) reduced FRET_eff (Figure 3C). At vesicles, raising intracellular calcium ([Ca²⁺]_I) levels significantly increased FRET_eff (non-stimulated: 3.10 ± 0.4, stimulated: 6.38 ± 0.71, p = 0.002), and led to a significant reduction after BAPTA-AM (0.30 ± 0.30, p = 0.9 × 10^-4), or W-7 application (0.88 ± 0.56, p = 3.1 × 10^-5) (Figure 3D). Similar to our measurements at the GJP, CaM_E1234Q at the vesicles showed a significantly reduced response upon [Ca²⁺]_I increase (non-stimulated: 2.18 ± 0.54, stimulated: 0.45 ± 0.27, n = 21, p = 0.005) (Figure 3D).

Next, FRET-Distance (FRET_dis) for Cx36–CFP and CaM–DsRed was compared at GJPs and the vesicles. Typically, if two interacting proteins are within 10 nm distance, FRET can be achieved (threshold indicated by the dashed line). No significant difference at the GJPs (non-stimulated: 9.11 ± 0.24, stimulated: 9.61 ± 0.33, p = 0.33), but a significant reduction in FRET_dis was observed at cytoplasmic vesicles when [Ca²⁺]_I was increased, (non-stimulated: 9.27 ± 0.22, stimulated: 8.09 ± 0.19, p = 0.002). When CaM_E1234Q was used to test the interaction with Cx36 WT-CFP, neither at the GJP (non-stimulated: 10.93 ± 0.63, stimulated: 10.13 ± 0.53, p = 0.21) nor the vesicles (non-stimulated: 9.96 ± 0.58, stimulated: 10.81 ± 0.46, p = 0.23) was significantly different (Figure 3E). All values were below 10 nm except for CaM_E1234Q, which indicated FRET was not achieved between wild-type Cx36 and CaM_E1234A. We concluded that changes in [Ca²⁺]_I are relayed to Cx36 before reaching GJPs, and that the interaction between Cx36 and CaM is critically dependent on the capacity of CaM to bind calcium.

W277 Is Critical for Calcium Dependent Cx36–CaM Binding

The role of the amino acid tryptophan in position 277 (W277) of Cx36 was investigated in calcium-dependent CaM binding together with the flanking amino acids 276 and 278. FRET_eff was determined for tagged Cx36 wild-type and Cx36 mutants in vesicles at physiological conditions, and after increasing [Ca²⁺]_I by stimulation with 2 mM [Ca²⁺]_E and 2 μM ionomycin. The FRET pair Cx36–ECFP and Cx36–DsRed served as the positive control (Figure 4). FRET_eff significantly increased in Neuro2a cells co-expressing CaM–DsRed and wild-type Cx36–ECFP (all values in %; non-stimulated: 3.65 ± 0.54, stimulated: 8.82 ± 1.19, p = 0.0003), the Cx36 mutants G276A (non-stimulated: 6.61 ± 0.63, stimulated: 9.74 ± 1.02, p = 0.02), or R278A (non-stimulated: 2.84 ± 0.29, n = 15; stimulated: 6.54 ± 1.02, n = 17, p = 0.003) (Figure 4A). In contrast, FRET_eff for W277A and CaM–DsRed were unaffected by the intracellular Ca²⁺ change (non-stimulated: 3.18 ± 0.39, n = 18; stimulated: 2.90 ± 0.69, p = 0.45).

FIGURE 4

FIGURE 4. Alanine scanning mutagenesis reveals W277 as the critical site for Cx36–CaM binding. (A) FRET_eff at vesicles in Neuro2a cells. Cells expressed fluorescently tagged Cx36–ECFP and Cx36–DsRed monomer as the control. Wild-type Cx36–DsRed monomer and Cx36 mutants (G276A, W277A, R278A) in the presence of CaM were stimulated with 2 μM of ionomycin for 10 min. (B) FRET distance with the same parameters as FRET efficiency. The dashed line indicates the FRET distance threshold at 10 nm. Error bars box plot: maximum and minimum range with median; Mann–Whitney U test, ^∗∗p < 0.05, ^∗∗∗p < 0.01.

This result was consistent with changes to the FRET distances (FRET_dis) between tagged Cx36 and CaM after stimulation with 2 mM [Ca²⁺]_E and 2 μM ionomycin. Wild-type Cx36 (all values in nm; non-stimulated: 9.15 ± 0.26, stimulated: 7.69 ± 0.22, p = 0.0003), mutant G276A (non-stimulated: 8.16 ± 0.23, stimulated: 7.48 ± 0.15, p = 0.02), or mutant R278A (non-stimulated: 9.33 ± 0.19, stimulated: 8.26 ± 0.24, p = 0.003) showed a significant decrease in FRET_dis, indicating that binding between both protein pairs increased (Figure 4B). In contrast, the FRET_dis for the Cx36 mutant W277A and CaM increased when intracellular Ca²⁺ was raised (non-stimulated: 9.22 ± 0.20, stimulated: 10.61 nm ± 0.91, p = 0.45). In summary, the FRET analysis suggested that tryptophan 277 is critical for CaM binding and that binding is sensitive to intracellular Ca²⁺ levels.

Reduced Dye-Redistribution in Neuro2a Cells Expressing W277A

The W277A mutation decreased the GJP frequency, but not the GJP size (Figures 1D,E), and caused a reduced response to [Ca²⁺]_I elevation (Figures 4A,B). To investigate impaired gap junction communication prompted by the W277A mutation, the non-invasive uptake and redistribution of ethidium bromide across GJPs under normal cell growth conditions was quantified using a method based on gap-FRAP (Abbaci et al., 2007). Pairs of Neuro2a cells expressing Cx36–ECFP or W277A–ECFP, and forming GJPs were selected for imaging (Figures 5A–D, left column). Sample traces for each experimental condition show the fluorescence in regions of interest (R1, red trace; R2, black trace) over the entire time course (Figures 5A–D, right column). Regions of interest were placed close to the GJP (R1) and at a more distant location (R2) in the same cell to capture recovery of ethidium bromide after photobleaching. Quantifications shown below (Figure 5E) demonstrated that treatment with 2 mM [Ca²⁺]_E and 2 μM ionomycin caused a significant ethidium bromide fluorescence recovery in Neuro2a cells expressing the Cx36–ECFP protein (values in %: non-stimulated: 5.82 ± 4.22, stimulated: 47.54 ± 13.56, p = 0.0045). Blocking of CaM with W-7 led to a significant reduction of fluorescence recovery when intracellular calcium levels were increased suggesting that CaM binding is critical for this response (non-stimulated: 16.60 ± 8.01, stimulated: 3.77 ± 3.41, p = 0.1427). Treatment with the gap junction blocker carbenoxolone (CBX) modulated fluorescence recovery (non-stimulated: -8.66 ± 4.65, stimulated: 8.91 ± 1.89, p = 0.0032). When Neuro2a cells expressed Cx36–W277A, no significant recovery was seen after the same treatment (non-stimulated: 8.12 ± 4.21, stimulated: -0.55 ± 8.35, p = 0.425). It was concluded that W277A expression impairs dye redistribution between cell pairs in a way similar to a gap junction blocker. Control experiments using non-transfected Neuro2a cells with and without ionomycin treatment showed limited recovery (Neuro2a without ionomycin treatment: -15.72 ± 17.52, Neuro2a with ionomycin treatment: 3.40 ± 1.23, p = 0.001). The recovery observed after CBX or ionomycin treatment likely reflects alterations to the biophysical properties of the cell membranes caused by the compounds applied.

FIGURE 5

FIGURE 5. Ethidium bromide dye uptake and redistribution after photobleaching. (A–D) Cell pairs expressing Cx36 wild-type proteins or the W277A mutant. The left column highlights examples for each treatment condition, with the cells selected for photobleaching encircled (large circles) and cell positions of regions of interest indicated (small circles; R1, R2, and R0 = background). The right column shows traces representing fluorescence over time for regions of interests. Scale bars: 10 μm. (E) Quantification of fluorescence recovery in regions of interest. Treatments: Iono (2 mM [Ca²⁺]_E and 2 μM ionomycin), W-7 (10 μM W-7, calmodulin inhibitor), CBX (50 μM carbenoxolone, gap junction blocker). Error bars display the maximum and minimum range. NS, not significant. Paired t-test, p-values ^∗∗ < 0.01, ^∗∗∗ < 0.005.

A High Resolution Solution Structure of Cx36 and CaM

Following a strategy used successfully in earlier structure determination (Ogura et al., 2012), a single-chain protein hybrid consisting of human CaM, a short glycine-rich linker, and the CaM binding region of rat Cx36 (Figure 6A) was expressed and uniformly ¹³C, ¹⁵N labeled for NMR spectroscopy. Amide ¹H, ¹⁵N HSQC spectra of this CaM–Cx36 hybrid protein (Supplementary Figure S6) and a sample of ¹⁵N-labeled CaM in a 1:1 complex with a wild-type 25 amino acid long Cx36 carboxy-terminal tail peptide were similar suggesting that the structure of the CaM–Cx36 hybrid was an accurate representation of a binary Cx36–CaM interaction.

FIGURE 6

FIGURE 6. Structural features of the Cx36–CaM complex. (A) A schematic representation of the hybrid Cx36–CaM protein used for the NMR structural determination. The N- and C-terminal lobes are colored for reference, followed by a short linker leading into a region in the carboxy-terminal cytoplasmic tail of Cx36. (B) Cartoon representation of the lowest energy model in the ensemble of twenty solutions. Four calcium ions are depicted as balls and other features colored consistently with previous panels. (C) Helical wheel representation of Cx36 bound to CaM as determined by experimental observations. W277 and V284 at anchor positions 1 and 8 are colored red and other hydrophobic amino acids colored black. (D) Sequence alignment of Cx36, smooth muscle myosin light chain kinase (smMLCK), and NMDA receptor NR1 (NMDAR NR1). Key contact sites between CaM and these ligands are made at positions 1, 8, and 14. Boxes indicate the helical amino acids in each ligand. The ligand may extend amino-terminally before anchor position 1 but for clarity, only amino acids from position 1 are shown. (E) Ligand plot (Wallace et al., 1995) of key interactions observed in the Cx36–CaM complex. Green lines denote hydrogen bonds and ionic interactions.

The Cx36–CaM hybrid protein structure was determined using 1100 distance restraints from ¹⁵N- and ¹³C-edited 3D-NOESY spectra, supplemented with a ¹²C-filtered, ¹³C-edited 2D-NOESY of a ¹³C-CaM/¹²C-Cx36 peptide complex (Table 2). The ensemble of 20 lowest energies structures is shown in Supplementary Figure S7. From the representative structure of the Cx36–CaM complex, CaM binds Cx36 in a well-known compact form (Figure 6B) with both Ca²⁺ saturated lobes enveloping a helical ligand with amphipathic features (Figure 6C). Based on the spacing of hydrophobic amino acids with W277 at anchor position 1, Cx36 resembles a 1–8–14 class ligand except that position 14 is occupied by K280 where a hydrophobic acid would typically be found. Accordingly, the Cx36–CaM protein structure compares favorably to other crystal structures of compact CaM complexes bearing 1–8–(14) ligands from smooth muscle myosin light chain kinase (Meador et al., 1992) (smMLCK; PDB:1CDL, 1.63 Å) and NMDA receptor NR1 C1 region (Ataman et al., 2007) (NMDAR NR1, PDB:2HQW, 1.81 Å) (Figure 6D). A contact map of Cx36–CaM from position 8 onward is shown in Figure 6E. The guanidino group of arginine at position 9 (R285 in Cx36) is balanced by ionic contacts from E84 in CaM. Likewise, the hydrophobic amino acid position 11 (A287 in Cx36) is supported by F19, L32, and V91 in CaM. At the penultimate contact position 14, the helix of Cx36 terminates earlier than smMLCK and NMDAR NR1 ligands. Despite this difference, three methionines in CaM (M36, M51, and M72) still contribute important hydrophobic contacts to Cx36. A lysine at position 14 in Cx36 is protruding from the CaM binding cleft toward a series of acidic amino acids in CaM (E47, D50, and E54). Taken together, a continuous network of contacts is observed between CaM and Cx36 thereby establishing Cx36 as bona fide ligand of CaM.

TABLE 2

TABLE 2. Statistics for the ensemble of 20 structures of the CaM–Cx36 hybrid protein.

Finally, ITC was performed to determine the affinity of CaM for a 25 amino acid long peptide derived from the Cx36 cytoplasmic tail (Figure 7). As predicted from the NMR structure, the Cx36 peptide bound Ca²⁺-CaM stoichiometrically with a K_d of 0.53 ± 0.05 μM (n = 3). Highlighting the role of the amino acid W277 at the position 1 of the 1–8–(14) motif, no binding was observed for the peptide bearing a W277A substitution (n = 3). Consistent with the structure and biological assays described throughout this study, peptide binding was only observed with Ca²⁺ saturated CaM.

FIGURE 7

FIGURE 7. Isothermal titration calorimetry (ITC) of Cx36 peptides and CaM. Calorimetric traces and integrated isotherms acquired at 30°C for (A) Calcium saturated CaM titrated with a Cx36 derived peptide, and (B) Calcium saturated CaM titrated with the same peptide bearing a W277A substitution. No binding was detected between CaM and W277A mutant peptide. Abbreviations: K_d, dissociation constant; n.d., not detected.

Discussion

The interaction of the neuronal gap junction protein Cx36 with CaM was investigated in Neuro2a cells and the amino acid W277 identified as a critical residue in the CaM binding motif. The binding of Cx36 to CaM was demonstrated in intracellular vesicles and at GJPs, but only the interaction at vesicles showed a calcium-dependent increase, revealed by pharmacological blocking of CaM, manipulation of intracellular calcium or overexpression of a calcium insensitive CaM protein. Mutation W277A, as well as pharmacological blocking of CaM caused a significant reduction in dye transfer between cell pairs, implicating CaM binding to Cx36 in gap junction communication. In vitro, CaM binds Cx36 with a comparable affinity to other connexins with contacts that follow a 1–8–14 motif found in other CaM-binding proteins such as smMLCK, NMDAR-NR1, and MARCKS.

The Cx36–CaM Binding Motif Is Distinct from Other Connexins

The novel details of CaM binding to rat Cx36 build on the original findings using mouse Cx36 and two perch Cx35 variants (Burr et al., 2005). The Cx36 binding motif is highly conserved across species from humans to lower vertebrates, like lamprey and zebrafish, and overlapping with a CaMKII binding site (Alev et al., 2008). In contrast, the CaM binding sites of Cx35/Cx36 are distinct from other connexins including Cx32, Cx43, Cx44, Cx46, and Cx50 in terms of ligand class (1–8 versus 1–5–10), location of the ligand binding sites, and resemble to date the only connexins activated by CaM (Zhou et al., 2007, 2009; Dodd et al., 2008; Myllykoski et al., 2009; Chen et al., 2011; Xu et al., 2012; Mruk et al., 2014).

The similarity of the Cx36 binding motifs across different species is reflected in comparable dissociation kinetics and thermodynamic characteristics found in mouse, the two perch Cx35 variants (Burr et al., 2005) and the rat Cx36 protein sequence investigated in this study. In more general terms, connexins bind CaM with dissociation kinetics and thermodynamic characteristics covering more than a 100-fold range (Dodd et al., 2008; Myllykoski et al., 2009; Zhou et al., 2009; Chen et al., 2011). Together, the specific localization of binding motifs, distinct subclasses and different dissociation kinetics suggest functional adaptations specific for each connexins protein.

The low affinity of Cx36 to CaM, at ∼500 nM, may have arisen due to the substitution of K290 at position 14 of typical 1–8–14 ligands like smMLCK, NMDAR-NR1, or MARCKS, which CaM encounters in the compact state. Furthermore, the less extensive array of hydrophobic contacts arising from the shorter side chains throughout Cx36, namely A283 and A287, cannot be compensated by larger hydrophobic amino at position 14 (K290 in Cx36). Regardless of the precise structural basis, a lower affinity may help Cx36 quickly sample CaM, CaMKII and possibly other proteins, depending upon the numerous signaling contexts it participates in throughout the life cycle of this protein.

Calmodulin Confers Calcium Signaling to Cx36 Outside the GJP

Generally, CaM binding to connexins has been investigated in the context of gap junction communication, mainly focusing on implications for channel functions in the cell membrane. On the basis of the subcellular distribution of Cx36 and mutants thereof, we have expanded this focus, showing that binding already occurs before Cx36 reaches the cell membrane forming gap junctions, or after removal from the GJP. Interestingly, the W277A mutation caused retention in vesicles and showed a reduced gap junction coupling, reflecting in part similar observations made for Cx43 (Zhou et al., 2007). Further, gap junctions made by W277A are less frequent and more prone to removal as annular plaques found distributed in cell pairs, suggesting that decorating Cx36 with CaM is critical and most likely serving multiple purposes during the life cycle of this connexin.

Ménage à Trios: Competition or Cooperation?

We have previously proposed a model predicting the sequential, cooperative processing of interaction between CaM, CaMKII, and Cx36 at the GJP (Alev et al., 2008). The main results of this study are more consistent with a sequential and competitive process, sequestering calcium-dependent actions between vesicles and GJPs. To emphasize the new findings, we propose a revised model for plasticity at electrical synapses. As highlighted in Figure 8, both CaM and CaMKII have the capacity to bind Cx36 at vesicles, with CaMKII binding at physiological calcium levels, and CaM binding only when intracellular calcium is elevated. At GJPs, CaMKII is calcium-responsive, with CaM showing no further strengthening of binding in response to calcium increase. This sequence of actions differs significantly from the activation of NMDA receptors (Bayer et al., 2001), and lead us to conclude that electrical synapses share a portfolio of interacting proteins, but mechanisms of achieving plasticity do not appear to fully recapitulate those found at chemical synapses. Instead, CaMKII binding at Cx36 at GJP seems to be the critical interaction conferring calcium-dependent plasticity at electrical synapses, while calcium-dependent sensitization of CaM appears restricted to spatially separate compartments. This could endow electrical synapses with a tunable molecular machinery using calcium as the main regulator for activation through CaMKII, while interaction of Cx36 with CaM open new avenues for linking this channel to intracellular calcium signaling pathways.

FIGURE 8

FIGURE 8. Model of calmodulin (CaM), Ca²⁺/calmodulin dependent kinase 2a (CaMKIIa) and connexin 36 (Cx36) in plasticity of the electrical synapse. Different to our original model, Cx36 can either bind calcium-activated CaM (right) or interact with CaMKIIa (left) at the vesicular level. The initial binding of CaM and Cx36 is critically dependent on elevated calcium concentrations. As the Cx36–CaM loaded vesicles traffic to the cellular membrane, the response to intracellular calcium activation shifts toward CaMKIIa interaction by calcium elevation. At the nexus, CaM serves as a kinase activator bound directly to CaMKIIa, and CaMKIIa plays an active part in modulating Cx36 nexus conductance (Illustration by Lesia Szyca, 2016).

What could be the function(s) of CaM binding to Cx36? The work by Burr et al. (2005) and our data suggest that CaM binding to Cx36 will occur when [Ca²⁺]_in levels are elevated. Such conditions typically are found during synaptic activity in neurons (Gamble and Koch, 1987). Since elevation of [Ca²⁺]_in increases Cx36–CaMKII, but not Cx36–CaM interaction at GJP, we conclude that calcium signals converge through CaM onto Cx36 before the protein reaches the GJP, ruling out a major role in plasticity.

It is reasonable to speculate that Cx36–CaM binding at vesicles may be sensitive to competition with other proteins with higher CaM affinity, suggesting that CaM binding could be orchestrated and/or sequestered by calcium affinities of CaM relative to CaMKII, or other, yet unidentified calcium binding proteins, or those activated or sensitized by calcium. Prior to its arrival at GJPs, Cx36–CaM could promote vesicular transport to the cell membrane and/or post-translational modifications affecting the oligomerization of connexins in the Golgi complex (Musil and Goodenough, 1993), as well as the assembly or disassembly of Cx36 into GJPs, potentially in a way similar to glutamate receptor trafficking to dendritic spines (Korkotian and Segal, 2007; Loo et al., 2014). Alternatively, calcium-dependent binding of CaM to Cx36 could link functional Cx36 channels in the ER directly to intracellular calcium signaling as suggested for Panx1, Panx3 or the transient receptor potential vanilloid 1 channel (TRPV1) (Vanden Abeele et al., 2006; Gallego-Sandin et al., 2009; Ishikawa et al., 2011), or by limiting calcium leakage from the ER as exemplified by CaM binding to Sec61 alpha (Erdmann et al., 2011). Addressing these consequences in future experiments will close this important knowledge gap.

Author Contributions

RS, ES, CB, CZ, and LD performed experiments. All authors analyzed data and were involved in the manuscript preparation. DS and GZ designed the study and finalized the manuscript.

Funding

This work was supported by the Natural Sciences and Engineering Research Council of Canada-Discovery Grants 238934-2011 (LD), 418142-2012 (GZ), and NIH NS092466 (DS).

Conflict of Interest Statement

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

Acknowledgment

We thank Ms. Lesia Szyca for generating the illustrations.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fnmol.2016.00120/full#supplementary-material

FIGURE S1 | Representative 3D images of Neuro2A cells expressing wild-type Cx36, tagged with EGFP, or DsRed monomer. The original projection based on the raw Z-stack before processing is shown top left. By rotating the 3D image and zooming into structural details of the Cx36 distribution it is possible to highlight the exchange of Cx36 proteins between adjacent cells connected by a GJP. Note that GJPs share a mixture of EGFP (in turquoise) and DsRed monomer (in red) tagged Cx36 proteins, whereas redistribution between cells as indicated by red and turquoise vesicles is low (Scale bar: 10 μm).

FIGURE S2 | Representative 3D images of Neuro2A cells expressing the Cx36 W277A mutant. The original projection based on the raw Z-stack before processing is shown top left. Note the mixed red and turquoise GJP and the increased redistribution of red and turquoise vesicles between adjacent cells indicating the reduced stability of the GJP and increase in protein turnover (Scale bar: 10 μm).

FIGURE S3 | Western blot detection of ER and Golgi tagged EGFP after transient expression in Neuro2a cells (in red, position indicated by lines). β-actin expression was used as internal protein standard.

FIGURE S4 | Representative 2D co-localization images of Neuro2A cells expressing wild-type EGFP tagged Cx36 and genetically encoded cell organelle markers for endoplasmic reticulum (ER), Golgi, and mitochondria Cx36. White areas indicate co-localization between Cx36 and organelle markers. (Scale bar: 10 μm).

FIGURE S5 | Representative 3D images of Neuro2a cells expressing wild-type Cx36–ECFP (in red) and DsRed monomer tagged mitochondria (in light turquoise). This image is a typical example showing the lack of co-localization between Cx36 and mitochondria (Scale bar: 10 μm).

FIGURE S6 | A ¹⁵N-edited HSQC spectum of ¹³C, ¹⁵N-labeled CaM–Cx36 hybrid protein acquired at 700 MHz. Cx36 specific amino acids are labeled red. Amino side chain resonances of asparagine and glutamine are indicated by dashed lines. An asterisk indicates an unassigned resonance. A plus sign indicates the position of a resonance below the contour limit of the plot.

FIGURE S7 | Superposition of the ensemble of 20 CaM–Cx36 NMR solution structures deposited in the Protein Data Bank (2N6A). The protein is presented as a Cα ribbon (CaM lobe 1, green; CaM lobe 2, blue; linker, gray; Cx36, orange) and the calcium ions are presented as spheres.

References

Abbaci, M., Barberi-Heyob, M., Stines, J. R., Blondel, W., Dumas, D., Guillemin, F., et al. (2007). Gap junctional intercellular communication capacity by gap-FRAP technique: a comparative study. Biotechnol. J. 2, 50–61. doi: 10.1002/biot.200600092

PubMed Abstract | CrossRef Full Text | Google Scholar

Alev, C., Urschel, S., Sonntag, S., Zoidl, G., Fort, A. G., Höher, T., et al. (2008). The neuronal connexin36 interacts with and is phosphorylated by CaMKII in a way similar to CaMKII interaction with glutamate receptors. Proc. Natl. Acad. Sci. U.S.A. 105, 20964–20969. doi: 10.1073/pnas.0805408105

PubMed Abstract | CrossRef Full Text | Google Scholar

Ataman, Z. A., Gakhar, L., Sorensen, B. R., Hell, J. W., and Shea, M. A. (2007). The NMDA receptor NR1 C1 region bound to calmodulin: structural insights into functional differences between homologous domains. Structure 15, 1603–1617. doi: 10.1016/j.str.2007.10.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Bayer, K. U., De Koninck, P., Leonard, A. S., Hell, J. W., and Schulman, H. (2001). Interaction with the NMDA receptor locks CaMKII in an active conformation. Nature 411, 801–805. doi: 10.1038/35081080

CrossRef Full Text | Google Scholar

Bennett, M. V., and Zukin, R. S. (2004). Electrical coupling and neuronal synchronization in the Mammalian brain. Neuron 41, 495–511. doi: 10.1016/S0896-6273(04)00043-1

CrossRef Full Text | Google Scholar

Bhattacharyya, M., Stratton, M. M., Going, C. C., Mcspadden, E. D., Huang, Y., Susa, A. C., et al. (2016). Molecular mechanism of activation-triggered subunit exchange in Ca(2+)/calmodulin-dependent protein kinase II. Elife 5:e13405. doi: 10.7554/eLife.13405

PubMed Abstract | CrossRef Full Text

Buchner, L., and Güntert, P. (2015). Increased reliability of nuclear magnetic resonance protein structures by consensus structure bundles. Structure 23, 425–434. doi: 10.1016/j.str.2014.11.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Burgoyne, R. D. (2007). Neuronal calcium sensor proteins: generating diversity in neuronal Ca2+ signalling. Nat. Rev. Neurosci. 8, 182–193. doi: 10.1038/nrn2093

PubMed Abstract | CrossRef Full Text | Google Scholar

Burr, G. S., Mitchell, C. K., Keflemariam, Y. J., Heidelberger, R., and O’brien, J. (2005). Calcium-dependent binding of calmodulin to neuronal gap junction proteins. Biochem. Biophys. Res. Commun. 335, 1191–1198. doi: 10.1016/j.bbrc.2005.08.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Cachope, R., and Pereda, A. E. (2012). Two independent forms of activity-dependent potentiation regulate electrical transmission at mixed synapses on the Mauthner cell. Brain Res. 1487, 173–182. doi: 10.1016/j.brainres.2012.05.059

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Y., Zhou, Y., Lin, X., Wong, H. C., Xu, Q., Jiang, J., et al. (2011). Molecular interaction and functional regulation of connexin50 gap junctions by calmodulin. Biochem. J. 435, 711–722. doi: 10.1042/BJ20101726

PubMed Abstract | CrossRef Full Text | Google Scholar

Deans, M. R., Gibson, J. R., Sellitto, C., Connors, B. W., and Paul, D. L. (2001). Synchronous activity of inhibitory networks in neocortex requires electrical synapses containing connexin36. Neuron 31, 477–485. doi: 10.1016/S0896-6273(01)00373-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Delaglio, F., Grzesiek, S., Vuister, G. W., Zhu, G., Pfeifer, J., and Bax, A. (1995). NMRPipe: a multidimensional spectral processing system based on UNIX pipes. J. Biomol. NMR 6, 277–293.

PubMed Abstract | Google Scholar

Del Corsso, C., Iglesias, R., Zoidl, G., Dermietzel, R., and Spray, D. C. (2012). Calmodulin dependent protein kinase increases conductance at gap junctions formed by the neuronal gap junction protein connexin36. Brain Res. 1487, 69–77. doi: 10.1016/j.brainres.2012.06.058

PubMed Abstract | CrossRef Full Text | Google Scholar

Dodd, R., Peracchia, C., Stolady, D., and Torok, K. (2008). Calmodulin association with connexin32-derived peptides suggests trans-domain interaction in chemical gating of gap junction channels. J. Biol. Chem. 283, 26911–26920. doi: 10.1074/jbc.M801434200

PubMed Abstract | CrossRef Full Text | Google Scholar

Erdmann, F., Schauble, N., Lang, S., Jung, M., Honigmann, A., Ahmad, M., et al. (2011). Interaction of calmodulin with Sec61alpha limits Ca2+ leakage from the endoplasmic reticulum. EMBO J. 30, 17–31. doi: 10.1038/emboj.2010.284

PubMed Abstract | CrossRef Full Text | Google Scholar

Erickson, M. G., Moon, D. L., and Yue, D. T. (2003). DsRed as a potential FRET partner with CFP and GFP. Biophys. J. 85, 599–611. doi: 10.1016/S0006-3495(03)74504-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Fliegel, L., Burns, K., Opas, M., and Michalak, M. (1989). The high-affinity calcium binding protein of sarcoplasmic reticulum. Tissue distribution, and homology with calregulin. Biochim. Biophys. Acta 982, 1–8. doi: 10.1016/0005-2736(89)90166-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Flores, C. E., Cachope, R., Nannapaneni, S., Ene, S., Nairn, A. C., and Pereda, A. E. (2010). Variability of distribution of Ca(2+)/calmodulin-dependent kinase II at mixed synapses on the mauthner cell: colocalization and association with connexin 35. J. Neurosci. 30, 9488–9499. doi: 10.1523/JNEUROSCI.4466-09.2010

PubMed Abstract | CrossRef Full Text | Google Scholar

Galarreta, M., and Hestrin, S. (1999). A network of fast-spiking cells in the neocortex connected by electrical synapses. Nature 402, 72–75. doi: 10.1038/47029

PubMed Abstract | CrossRef Full Text | Google Scholar

Gallego-Sandin, S., Rodriguez-Garcia, A., Alonso, M. T., and Garcia-Sancho, J. (2009). The endoplasmic reticulum of dorsal root ganglion neurons contains functional TRPV1 channels. J. Biol. Chem. 284, 32591–32601. doi: 10.1074/jbc.M109.019687

PubMed Abstract | CrossRef Full Text | Google Scholar

Gamble, E., and Koch, C. (1987). The dynamics of free calcium in dendritic spines in response to repetitive synaptic input. Science 236, 1311–1315. doi: 10.1126/science.3495885

PubMed Abstract | CrossRef Full Text | Google Scholar

Haas, J. S., Zavala, B., and Landisman, C. E. (2011). Activity-dependent long-term depression of electrical synapses. Science 334, 389–393. doi: 10.1126/science.1207502

PubMed Abstract | CrossRef Full Text | Google Scholar

Hamzei-Sichani, F., Davidson, K. G., Yasumura, T., Janssen, W. G., Wearne, S. L., Hof, P. R., et al. (2012). Mixed electrical-chemical synapses in adult rat hippocampus are primarily glutamatergic and coupled by connexin-36. Front. Neuroanat. 6:13. doi: 10.3389/fnana.2012.00013

PubMed Abstract | CrossRef Full Text | Google Scholar

Hormuzdi, S. G., Filippov, M. A., Mitropoulou, G., Monyer, H., and Bruzzone, R. (2004). Electrical synapses: a dynamic signaling system that shapes the activity of neuronal networks. Biochim. Biophys. Acta 1662, 113–137. doi: 10.1016/j.bbamem.2003.10.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Ishikawa, M., Iwamoto, T., Nakamura, T., Doyle, A., Fukumoto, S., and Yamada, Y. (2011). Pannexin 3 functions as an ER Ca(2+) channel, hemichannel, and gap junction to promote osteoblast differentiation. J. Cell Biol. 193, 1257–1274. doi: 10.1083/jcb.201101050

PubMed Abstract | CrossRef Full Text | Google Scholar

Korkotian, E., and Segal, M. (2007). Morphological constraints on calcium dependent glutamate receptor trafficking into individual dendritic spine. Cell Calcium 42, 41–57. doi: 10.1016/j.ceca.2006.11.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Kothmann, W. W., Trexler, E. B., Whitaker, C. M., Li, W., Massey, S. C., and O’brien, J. (2012). Nonsynaptic NMDA receptors mediate activity-dependent plasticity of gap junctional coupling in the AII amacrine cell network. J. Neurosci. 32, 6747–6759. doi: 10.1523/JNEUROSCI.5087-11.2012

PubMed Abstract | CrossRef Full Text | Google Scholar

Kurtenbach, S., Prochnow, N., Kurtenbach, S., Klooster, J., Zoidl, C., Dermietzel, R., et al. (2013). Pannexin1 channel proteins in the zebrafish retina have shared and unique properties. PLoS ONE 8:e77722. doi: 10.1371/journal.pone.0077722

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, H., Chuang, A. Z., and O’brien, J. (2009). Photoreceptor coupling is controlled by connexin 35 phosphorylation in zebrafish retina. J. Neurosci. 29, 15178–15186. doi: 10.1523/JNEUROSCI.3517-09.2009

CrossRef Full Text | Google Scholar

Li, J., Shen, H., Naus, C. C., Zhang, L., and Carlen, P. L. (2001). Upregulation of gap junction connexin 32 with epileptiform activity in the isolated mouse hippocampus. Neuroscience 105, 589–598. doi: 10.1016/S0306-4522(01)00204-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Lisman, J., Schulman, H., and Cline, H. (2002). The molecular basis of CaMKII function in synaptic and behavioural memory. Nat. Rev. Neurosci. 3, 175–190. doi: 10.1038/nrn753

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, C., and Hermann, T. E. (1978). Characterization of ionomycin as a calcium ionophore. J. Biol. Chem. 253, 5892–5894.

Google Scholar

Loo, L. S., Tang, N., Al-Haddawi, M., Dawe, G. S., and Hong, W. (2014). A role for sorting nexin 27 in AMPA receptor trafficking. Nat. Commun. 5:3176. doi: 10.1038/ncomms4176

PubMed Abstract | CrossRef Full Text | Google Scholar

Manders, E. M. M., Verbeek, F. J., and Aten, J. A. (1993). Measurement of co-localization of objects in dual-colour confocsal images. J. Microsc. 169, 375–382. doi: 10.1111/j.1365-2818.1993.tb03313.x

CrossRef Full Text | Google Scholar

Maune, J. F., Klee, C. B., and Beckingham, K. (1992). Ca2+ binding and conformational change in two series of point mutations to the individual Ca(2+)-binding sites of calmodulin. J. Biol. Chem. 267, 5286–5295.

PubMed Abstract | Google Scholar

Meador, W. E., Means, A. R., and Quiocho, F. A. (1992). Target enzyme recognition by calmodulin: 2.4 A structure of a calmodulin-peptide complex. Science 257, 1251–1255. doi: 10.1126/science.1519061

PubMed Abstract | CrossRef Full Text | Google Scholar

Mruk, K., Farley, B. M., Ritacco, A. W., and Kobertz, W. R. (2014). Calmodulation meta-analysis: predicting calmodulin binding via canonical motif clustering. J. Gen. Physiol. 144, 105–114. doi: 10.1085/jgp.201311140

PubMed Abstract | CrossRef Full Text | Google Scholar

Musil, L. S., and Goodenough, D. A. (1993). Multisubunit assembly of an integral plasma membrane channel protein, gap junction connexin43, occurs after exit from the ER. Cell 74, 1065–1077. doi: 10.1016/0092-8674(93)90728-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Myllykoski, M., Kuczera, K., and Kursula, P. (2009). Complex formation between calmodulin and a peptide from the intracellular loop of the gap junction protein connexin43: molecular conformation and energetics of binding. Biophys. Chem. 144, 130–135. doi: 10.1016/j.bpc.2009.08.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Nagy, J. I., Bautista, W., Blakley, B., and Rash, J. E. (2013). Morphologically mixed chemical-electrical synapses formed by primary afferents in rodent vestibular nuclei as revealed by immunofluorescence detection of connexin36 and vesicular glutamate transporter-1. Neuroscience 252, 468–488. doi: 10.1016/j.neuroscience.2013.07.056

PubMed Abstract | CrossRef Full Text | Google Scholar

Ogura, K., Kumeta, H., Takahasi, K., Kobashigawa, Y., Yoshida, R., Itoh, H., et al. (2012). Solution structures of yeast Saccharomyces cerevisiae calmodulin in calcium- and target peptide-bound states reveal similarities and differences to vertebrate calmodulin. Genes Cells 17, 159–172. doi: 10.1111/j.1365-2443.2012.01580.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Okamoto, K., Nagai, T., Miyawaki, A., and Hayashi, Y. (2004). Rapid and persistent modulation of actin dynamics regulates postsynaptic reorganization underlying bidirectional plasticity. Nat. Neurosci. 7, 1104–1112. doi: 10.1038/nn1311

PubMed Abstract | CrossRef Full Text | Google Scholar

Olmsted, J. B., Carlson, K., Klebe, R., Ruddle, F., and Rosenbaum, J. (1970). Isolation of microtubule protein from cultured mouse neuroblastoma cells. Proc. Natl. Acad. Sci. U.S.A. 65, 129–136.

Google Scholar

Pereda, A., O’brien, J., Nagy, J. I., Smith, M., Bukauskas, F., Davidson, K. G., et al. (2003). Short-range functional interaction between connexin35 and neighboring chemical synapses. Cell Commun. Adhes 10, 419–423. doi: 10.1080/714040462

PubMed Abstract | CrossRef Full Text | Google Scholar

Pereda, A. E. (2014). Electrical synapses and their functional interactions with chemical synapses. Nat. Rev. Neurosci. 15, 250–263. doi: 10.1038/nrn3708

PubMed Abstract | CrossRef Full Text | Google Scholar

Pereda, A. E., Bell, T. D., Chang, B. H., Czernik, A. J., Nairn, A. C., Soderling, T. R., et al. (1998). Ca2+/calmodulin-dependent kinase II mediates simultaneous enhancement of gap-junctional conductance and glutamatergic transmission. Proc. Natl. Acad. Sci. U.S.A. 95, 13272–13277. doi: 10.1073/pnas.95.22.13272

PubMed Abstract | CrossRef Full Text | Google Scholar

Pereda, A. E., Rash, J. E., Nagy, J. I., and Bennett, M. V. (2004). Dynamics of electrical transmission at club endings on the Mauthner cells. Brain Res. Brain Res. Rev. 47, 227–244. doi: 10.1016/j.brainresrev.2004.06.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Postma, F., Liu, C. H., Dietsche, C., Khan, M., Lee, H. K., Paul, D., et al. (2011). Electrical synapses formed by connexin36 regulate inhibition- and experience-dependent plasticity. Proc. Natl. Acad. Sci. U.S.A. 108, 13770–13775. doi: 10.1073/pnas.1100166108

PubMed Abstract | CrossRef Full Text | Google Scholar

Prochnow, N., Hoffmann, S., Dermietzel, R., and Zoidl, G. (2009). Replacement of a single cysteine in the fourth transmembrane region of zebrafish pannexin 1 alters hemichannel gating behavior. Exp. Brain Res. 199, 255–264. doi: 10.1007/s00221-009-1957-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Qu, C., Gardner, P., and Schrijver, I. (2009). The role of the cytoskeleton in the formation of gap junctions by Connexin 30. Exp. Cell Res. 315, 1683–1692. doi: 10.1016/j.yexcr.2009.03.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Rizzuto, R., Nakase, H., Darras, B., Francke, U., Fabrizi, G. M., Mengel, T., et al. (1989). A gene specifying subunit VIII of human cytochrome c oxidase is localized to chromosome 11 and is expressed in both muscle and non-muscle tissues. J. Biol. Chem. 264, 10595–10600.

PubMed Abstract | Google Scholar

Rubio, M. E., and Nagy, J. I. (2015). Connexin36 expression in major centers of the auditory system in the CNS of mouse and rat: Evidence for neurons forming purely electrical synapses and morphologically mixed synapses. Neuroscience 303, 604–629. doi: 10.1016/j.neuroscience.2015.07.026

PubMed Abstract | CrossRef Full Text | Google Scholar

Skinner, S. P., Goult, B. T., Fogh, R. H., Boucher, W., Stevens, T. J., Laue, E. D., et al. (2015). Structure calculation, refinement and validation using CcpNmr Analysis. Acta Crystallogr. D Biol. Crystallogr. 71, 154–161. doi: 10.1107/S1399004714026662

PubMed Abstract | CrossRef Full Text | Google Scholar

Ueda, N., and Shah, S. V. (1992). Role of intracellular calcium in hydrogen peroxide-induced renal tubular cell injury. Am. J. Physiol. 263, F214–F221.

Google Scholar

Van Der Giessen, R. S., Koekkoek, S. K., Van Dorp, S., De Gruijl, J. R., Cupido, A., Khosrovani, S., et al. (2008). Role of olivary electrical coupling in cerebellar motor learning. Neuron 58, 599–612. doi: 10.1016/j.neuron.2008.03.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Vanden Abeele, F., Bidaux, G., Gordienko, D., Beck, B., Panchin, Y. V., Baranova, A. V., et al. (2006). Functional implications of calcium permeability of the channel formed by pannexin 1. J. Cell Biol. 174, 535–546. doi: 10.1083/jcb.200601115

PubMed Abstract | CrossRef Full Text | Google Scholar

VanScyoc, W. S., Sorensen, B. R., Rusinova, E., Laws, W. R., Ross, J. B., and Shea, M. A. (2002). Calcium binding to calmodulin mutants monitored by domain-specific intrinsic phenylalanine and tyrosine fluorescence. Biophys. J. 83, 2767–2780. doi: 10.1016/S0006-3495(02)75286-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Vervaeke, K., Lorincz, A., Gleeson, P., Farinella, M., Nusser, Z., and Silver, R. A. (2010). Rapid desynchronization of an electrically coupled interneuron network with sparse excitatory synaptic input. Neuron 67, 435–451. doi: 10.1016/j.neuron.2010.06.028

PubMed Abstract | CrossRef Full Text | Google Scholar

Wallace, A. C., Laskowski, R. A., and Thornton, J. M. (1995). LIGPLOT: a program to generate schematic diagrams of protein-ligand interactions. Protein Eng. 8, 127–134. doi: 10.1093/protein/8.2.127

PubMed Abstract | CrossRef Full Text | Google Scholar

Watzele, G., and Berger, E. G. (1990). Near identity of HeLa cell galactosyltransferase with the human placental enzyme. Nucl. Acids Res. 18:7174. doi: 10.1093/nar/18.23.7174

PubMed Abstract | CrossRef Full Text | Google Scholar

Weiss, B., Prozialeck, W. C., and Wallace, T. L. (1982). Interaction of drugs with calmodulin. Biochemical, pharmacological and clinical implications. Biochem. Pharmacol. 31, 2217–2226.

PubMed Abstract | Google Scholar

Xu, Q., Kopp, R. F., Chen, Y., Yang, J. J., Roe, M. W., and Veenstra, R. D. (2012). Gating of connexin 43 gap junctions by a cytoplasmic loop calmodulin binding domain. Am. J. Physiol. Cell Physiol. 302, C1548–C1556. doi: 10.1152/ajpcell.00319.2011

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, Y., Yang, W., Lurtz, M. M., Chen, Y., Jiang, J., Huang, Y., et al. (2009). Calmodulin mediates the Ca2+-dependent regulation of Cx44 gap junctions. Biophys. J. 96, 2832–2848. doi: 10.1016/j.bpj.2008.12.3941

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, Y., Yang, W., Lurtz, M. M., Ye, Y., Huang, Y., Lee, H. W., et al. (2007). Identification of the calmodulin binding domain of connexin 43. J. Biol. Chem. 282, 35005–35017. doi: 10.1074/jbc.M707728200

PubMed Abstract | CrossRef Full Text | Google Scholar

Zoidl, G., Meier, C., Petrasch-Parwez, E., Zoidl, C., Habbes, H. W., Kremer, M., et al. (2002). Evidence for a role of the N-terminal domain in subcellular localization of the neuronal connexin36 (Cx36). J. Neurosci. Res. 69, 448–465. doi: 10.1002/jnr.10284

PubMed Abstract | CrossRef Full Text | Google Scholar