Effects of Phosphorylation on the Structure and Backbone Dynamics of the Intrinsically Disordered Connexin43 C-terminal Domain*

Background: Phosphorylation of the connexin43 C-terminal (Cx43CT) domain regulates gap junction intercellular communication (GJIC). Results: Phosphorylation alters the α-helical propensity of the Cx43CT. Conclusion: Altering the conformational preference of the Cx43CT presents a novel mechanism for regulation of GJIC. Significance: Cx43CT residues susceptible to structural alterations are prime targets for chemical modulators of GJIC. Phosphorylation of the connexin43 C-terminal (Cx43CT) domain regulates gap junction intercellular communication. However, an understanding of the mechanisms by which phosphorylation exerts its effects is lacking. Here, we test the hypothesis that phosphorylation regulates Cx43 gap junction intercellular communication by mediating structural changes in the C-terminal domain. Circular dichroism and nuclear magnetic resonance were used to characterize the effects of phosphorylation on the secondary structure and backbone dynamics of soluble and membrane-tethered Cx43CT domains. Cx43CT phospho-mimetic isoforms, which have Asp substitutions at specific Ser/Tyr sites, revealed phosphorylation alters the α-helical content of the Cx43CT domain only when attached to the membrane. The changes in secondary structure are due to variations in the conformational preference and backbone flexibility of residues adjacent and distal to the site(s) of modification. In addition to the known direct effects of phosphorylation on molecular partner interactions, the data presented here suggest phosphorylation may also indirectly regulate binding affinity by altering the conformational preference of the Cx43CT domain.

rapid inter-conversion of an IDP between heterogeneous conformations confers several unique kinetic and thermodynamic advantages (11)(12)(13). The large hydrodynamic radius of a disordered region can increase the capture radius of the protein (14,15). "Fly casting" contributes to the high k on rates often associated with IDPs (16). The ability to adopt multiple conformations enables the molecular promiscuity typical of many IDPs (17). From a thermodynamic perspective, the entropic cost of coupled folding and binding results in high specificity, low affinity interactions, which are advantageous in cellular signaling and regulatory response (12,13). Additionally, many IDPs have a conformational preference consistent with the bound state (14). Altering the conformational preference of an IDP could modulate the favorability of a protein partner interaction (18). One cellular event regulating protein function that causes significant structural rearrangement and is disproportionately high for IDPs is phosphorylation (9,19).
The Cx43CT is differentially phosphorylated at a dozen or more residues at various times throughout the Cx43 life cycle (20 -24). Many kinases have been identified to phosphorylate Cx43 with the predominant effect being a decrease in GJIC. For example, following an ischemic event, the resulting intracellular acidification stimulates changes in the Cx43CT phosphorylation, which leads to closure of gap junction channels (25). Unfortunately, an understanding of the mechanism(s) by which Cx43 phosphorylation alters channel communication is lacking. This is best exemplified by the functional opposition of Cx43 phosphorylation at Ser-365 (PKA-mediated) and Ser-368 (PKC). Although the sites are in close proximity, phosphorylation of Ser-368 results in closure of Cx43 gap junction channels, whereas phosphorylation at Ser-365 increases GJIC and acts as a gatekeeper by blocking Ser-368 phosphorylation-induced channel closure (26). Evidence has emerged that a phosphate can directly affect binding of the Cx43CT domain with molecular partners by blocking the binding site (e.g. ZO-1 (27) and tubulin (28)) or by being directly involved in the interaction (e.g. Nedd4 (29,30)). Alternatively, phosphorylation could affect indirectly interactions between IDPs and protein partners by altering the conformational preference of the disordered region (18,31), a mechanism that has not been previously investigated for Cx43.
In this study, CD and NMR were used to characterize global and local effects of phosphorylation on the secondary structure and backbone dynamics of soluble and membrane-tethered Cx43CT domains. Asp substitutions were used to mimic phosphorylation by kinases that play an important role in the regulation of Cx43 (Table 1) (for review, see Ref. 32). The results indicate phosphorylation can alter the ␣-helical propensity of the Cx43CT at residues proximal and distal to the site(s) of modification. The changes in secondary structure are due to variations in the conformational preference and backbone flexibility of residues adjacent and distal to the site(s) of modification. Changes in the conformational preference of an IDP is one mechanism by which protein partner interactions can be regulated and it represents a novel mechanism for regulation of GJIC.

EXPERIMENTAL PROCEDURES
Protein Purification-The expression and purification of the soluble Cx43CT (Val-236 -Ile-382) and Cx43CT domain when attached to the 4th transmembrane domain (TM4-Cx43CT; Asp-197-Ile-382) in detergent micelles have been previously described (33,34). Phospho-mimetic isoforms were created by substituting sites of Ser/Tyr phosphorylation (Table 1) with Asp using site-directed mutagenesis (QuikChange Lightning Kit, Agilent Technologies). Purifications of the phospho-mimetic isoforms were done using the same protocols as the soluble and membrane-tethered wild-type (WT) proteins. All plasmid sequences were verified at the University of Nebraska Medical Center DNA Sequencing Core Facility. Protein purity was assessed using SDS-PAGE.
Nuclear Magnetic Resonance (NMR)-NMR data were acquired using a 600 MHz Varian INOVA NMR spectrometer fitted with a cryo-probe (Agilent Technologies). All NMR data for the soluble Cx43CT (200 M) and TM4-Cx43CT (800 M) were collected at 7 and 42°C, respectively. Backbone assignments of the soluble Cx43CT were obtained using the following three-dimensional experiments: HNCACB, HN(CO)CACB, HNCO, HN(CA)CO, and HNCACO. NMR data for the TM4-Cx43CT constructs were obtained in the presence of 10% TFE. Gradient-enhanced 15 N HSQC experiments used to observe backbone amide resonances were acquired with 1024 complex points in the direct dimension and 256 complex points in the indirect dimension. 15 N NOESY-HSQC data were collected with a mixing time of 50 ms. The backbone 15 N R 1 (inversion recovery) and R 2 (Carr-Purcell-Meinboom-Gill) NMR relaxation measurements were collected with 1024 and 128 complex points in 1 H and 15 N, respectively, with 32 scans per R 1 point and a recycle delay of 1 s. The R 1 longitudinal relaxation delays were 10,20,40,80,140,200, 300, 400, 600, 800, 1200, 2000, and 4000 ms. The R 2 relaxation delays were 10,30,50,70,110,150,190, and 230 ms, with a 1-ms delay between 15 N pulses. The 1 H-15 N NOE data were acquired with and without 1 H saturation. All NMR spectra were processed using NMRPipe (39) and analyzed using NMRView (40). Relaxation data were processed by rate analysis in NMRView.

Rationale for Asp-substituted Cx43CT
Isoforms-Phosphorylation of connexins has been implicated in the regulation of GJIC. However, several problems have contributed to the difficulty of elucidating the mechanism(s) by which a single kinase affects GJIC as follows: the transient nature of a particular Cx43CT phosphorylation state, the ability of many kinases to phosphorylate more than one Cx43CT residue, the ability of various kinases to phosphorylate the Cx43CT at the same time, and the inability to precisely control which residue(s) are phosphorylated (of the 27 Ser/Tyr residues in the Cx43CT, all but 8 have been observed to be phosphorylated (41,42)). Strategies employed to overcome these difficulties include the use of phosphorylation-specific Cx43CT antibodies (43) and short Cx43CT phosphopeptides (27,28). Also well appreciated is the use of Asp (or Glu) substitutions as a mimetic for phosphorylation. Asp substitutions allow for the precise control of the site(s) modified and the production of enough protein for biophysical studies. Although the size and charge of an Asp(Ϫ1) is not identical to a phosphate (Ϫ2 at pH 7.5 and ϾϪ1 at pH 5.8), in vivo and in vitro studies have confirmed Asp (or Glu) substitutions can elicit a response similar to a phosphate in Cx43 gap junction channels (26,44,45). Data were collected at either physiological (pH 7.5) or acidic (pH 5.8) conditions. Intracellular acidification, which is a major consequence of tissue ischemia during a myocardial infarction, leads to closure and degradation of gap junction channels and can be a substrate for malignant ventricular arrhythmias (46).
Electrophoretic Mobility of Asp-substituted Soluble Cx43CT Isoforms-Numerous studies have shown that differentially phosphorylated Cx43 results in multiple electrophoretic isoforms as follows: a fast migrating isoform, which includes the nonphosphorylated isoform (P0), and multiple slower migrating isoforms (P1 and P2) (20,47,48). Here, Asp substitutions were cloned into the soluble version of Cx43CT (Val-236 -Ile-382) and each phospho-mimetic isoform was analyzed by SDS-PAGE (Fig. 1). The data indicate that the wild-type (WT) Cx43CT (nonphosphorylated) and the Src, PKC, and cdc2 phospho-mimetic isoforms migrate at P0, MAPK at P1, and CK1 at P2 (Fig. 1A, top and middle). These results are consistent with the known phosphorylation state(s) that compose the P0, P1, and P2 isoforms. Thus, Asp-substituted Cx43CT isoforms mimic characteristics of endogenously phosphorylated Cx43. PKA-mediated phosphorylation at Ser-365 is typically present in the P1 isoform; however, the soluble Cx43CT S365D phospho-mimetic isoform migrated at P0 (Fig. 1A, top). In addition to Ser-365, PKA can also phosphorylate Cx43 at Ser-364, Ser-369, and Ser-373. Therefore, the migration pattern for each PKA-mediated phosphorylation site alone and together was analyzed (Fig. 1A, bottom). The data reveal a shift in migration from P0 to P1 is dependent upon PKA-mediated phosphorylation at multiple Cx43CT sites. The effect of single-site phosphorylation by each of the kinases was also tested, and none caused a shift in the migration pattern (results are summarized in Fig. 1B).
Phosphorylation-dependent shifts in electrophoretic mobility are not phenomena unique to connexins. Protein phosphorylation occurring on residues adjacent to a Pro has been reported to induce SDS-resistant changes in structure that alter electrophoretic mobility (49). The altered mobility is likely due to structural perturbations that interfere with the formation of an "ideal" random coil in the presence of SDS. Additionally, the negative charges on the phosphate group can interfere with protein-SDS interactions, thereby disrupting the charge-tomass ratio and altering the electrophoretic mobility of the protein. However, this mechanism requires a number of phosphates to be close together. For the Cx43CT domain, both mechanisms likely contribute to the electrophoretic mobility of the phospho-isoforms. Next, the possibility that phosphorylation alters the secondary structure of the soluble Cx43CT domain was investigated by evaluating the residue-specific effects of phosphorylation using NMR.
Secondary Structural Studies Using the Soluble Cx43CT Domain-The structure of the soluble Cx43CT (Ser-255-Ile-382) was determined by NMR (45). In this study, a longer construct (Val-236 -Ile-382) is used, which incorporates the other Src phosphorylation site at Tyr-247. All WT Cx43CT backbone 1 H, 15 N, and 13 C resonances (147 residues; 15 Pro) were assigned to identify specific residues and/or regions that are affected by phosphorylation. Fig. 2A is the assigned WT Cx43CT 15 N HSQC spectrum. When the WT Cx43CT 15 N HSQC spectrum is compared with spectra from all the phospho-mimetic isoforms, several residues shifted. However, these shifts are small, of similar line width, and localized next to the site of substitution (e.g. CK1 isoform; Fig. 2B). The small chemical shift differences could be the result of altered secondary structure or due to changes in the chemical environment of residues neighboring the Asp substitutions. Therefore, the CD spectrum of the WT Cx43CT was compared with each of the phospho-mimetic isoforms to determine whether the differences detected by NMR were the result of a change in secondary structure.
CD spectra of proteins have distinct MRE values for random coil, ␤-sheet, and ␣-helical structures (50). The CD spectrum of the soluble WT Cx43CT at pH 5.8 has a peak minimum near 198 nm and a small minimum at 222 nm ( Fig. 3A), characteristic of an IDP with a small amount of ␣-helical content. This observation is consistent with previous CD data and the NMR solution structure of the shorter Cx43CT (Ser-255-382) domain (7,51). Surprisingly, the CD spectra for all of the phospho-mimetic isoforms are nearly identical to the WT Cx43CT at pH 5.8 (Fig. 3A). The phospho-mimetic isoforms are also compared at pH 7.5 ( Fig. 3B) to determine whether a phosphorylation-induced structural change was contingent upon pH (46). The CD data indicate phosphorylation does not alter the secondary structure of the soluble Cx43CT. This suggests the differences in NMR chemical shifts are solely due to changes in the amino acid sequence. The possibility exists that the lack of a phosphorylation-induced structural change is due to the increased flexibility of the soluble Cx43CT domain.
Asp-substituted Membrane-tethered Cx43CT Isoforms-Our previous study determined that the Cx43CT, when attached to the 4th transmembrane domain (TM4-Cx43CT), has more native-like characteristics than the soluble Cx43CT (34). Therefore, the TM4-Cx43CT construct was used to investigate the effects of phosphorylation on the secondary structure of the Cx43CT. SDS-PAGE analysis of the soluble Cx43CT indicated multiple site modifications are needed to mimic the shifts seen with endogenously phosphorylated Cx43 (Fig. 1B). Therefore, multiple Asp substitutions were cloned into the TM4-Cx43CT to mimic complete phosphorylation by the kinases listed in Table 1, except MAPK. For the MAPK isoform (S255D,S279D,S282D), various attempts (e.g. changes in growth and expression conditions; various cell lines), including a double substitution (S279D, S282D), were made to express the phospho-mimetic isoform without success. Two MAPK single substitutions were expressed and purified, TM4-Cx43CT S279D and TM4-Cx43CT S282D . The TM4-Cx43CT S282D was chosen for characterization given that phosphorylation at Ser-282 is sufficient to induce Cx43 channel closure (52). The membrane-tethered phospho-isoforms exhibited SDS-PAGE mobility shifts consistent with the soluble Cx43CT phospho-isoforms Cx43 (Fig. 1C). Notably, the different electrophoretic mobility between MAPK S282D (Fig. 1C, TM4-Cx43CT) and S255D,S279D,S282D (Fig. 1A, soluble Cx43CT) is likely caused by the number of negative charges needed to cause a shift in the gel.
Secondary Structural Studies Using the Membrane-tethered Cx43CT Domain-In Fig. 4, the CD spectrum of the WT TM4-Cx43CT is compared with the phospho-mimetic isoforms at pH 5.8 and 7.5. MRE at 222 nm was used as a gauge for changes in ␣-helical content (50). The data indicate phosphorylation can alter both the ␣-helical content of the TM4-Cx43CT as well as the structural responsiveness to acidification ( Table 2). The CK1 and MAPK isoforms (Fig. 4, A and B) are more ␣-helical than the WT TM4-Cx43CT at pH 7.5 and exhibit an increase in ␣-helical content upon acidification. Conversely, the ␣-helical content of the Src isoform is similar to the WT TM4-Cx43CT at pH 7.5 ( Fig. 4C) but is unaffected by a decrease in pH. Likewise, the PKA-mediated and PKC isoforms show little to no change in secondary structure upon acidification (Fig. 4, D and E), but the ␣-helical content is similar to the WT TM4-Cx43CT at pH 5.8. Finally, the ␣-helical content of the cdc2 isoform at both pH values (Fig. 4F) exceeds WT TM4-Cx43CT at pH 5.8 (similar to the CK1 and MAPK isoforms), but it undergoes little to no change upon acidification (similar to the Src isoform). Including the observation that pH does not affect the secondary structure of the TM4 domain (residues Asp-197-Val-240; Fig. 3C), the data suggest phosphorylation alters the secondary structure of the membrane-tethered Cx43CT domain.
To complement the global perspective provided by the CD data, NMR was used to investigate phosphorylation-mediated changes in secondary structure on a residue-by-residue basis. Previously, we determined the NMR backbone assignment of the WT TM4-Cx43CT at pH 5.8 in the presence of 10% (v/v) TFE (53). TFE was used to stabilize innate ␣-helical content and improve NMR spectra (33). Importantly, 10% TFE did not induce ␣-helical structure in the WT TM4-Cx43CT (34) or the phospho-mimetic isoforms as assayed by CD (data not shown). Therefore, NMR 15 N HSQC spectra of all TM4-Cx43CT phospho-mimetic isoforms were collected and superimposed onto the TM4-Cx43CT WT spectrum (Fig. 5, A-F). All affected Cx43CT residues from each phospho-mimetic isoform were mapped onto the Cx43CT sequence ( Fig. 5G; Table 3). Many of the peaks that were altered in the phospho-mimetic isoforms either disappeared or were unable to be identified due to such large changes in chemical shift. Unlike the soluble Cx43CT, chemical shift changes in the TM4-Cx43CT phospho-mimetic isoforms were localized not only to regions adjacent to the site of amino acid substitution but propagated throughout the entire Cx43CT domain.
The 15 N HSQC spectra at pH 5.8 of the CK1, MAPK, and cdc2 isoforms each show upfield peak shifts (Fig. 5, A-C,  arrows). Upfield shifts were expected because these isoforms have more ␣-helical content (as determined by CD) than the WT TM4-Cx43CT. Many of the CK1, MAPK, and cdc2 isoform peaks are weak or have increased line widths (Fig. 5, A-C, circles). The increase in line width suggests these residues are in exchange between a disordered and ordered structure. The Src isoform has several peaks that shift downfield and/or increase line widths compared with the WT TM4-Cx43CT at pH 5.8 (Fig. 5D, arrows). These observations are consistent with the  AUGUST 23, 2013 • VOLUME 288 • NUMBER 34

JOURNAL OF BIOLOGICAL CHEMISTRY 24861
CD data; Src phosphorylation destabilizes ␣-helical structure to favor a disordered conformation. The 15 N HSQC for the PKC isoform is similar to the WT TM4-Cx43CT (Fig. 5E) with only a few chemical shift differences near the sites of Asp substitution. This observation is also consistent with the CD data. Conversely, the 15 N HSQC spectrum of the PKA-mediated isoform is different (Fig. 5F), even though the CD data indicate the ␣-helical content is similar to both the WT TM4-Cx43CT and PKC isoforms. Numerous 15 N HSQC peaks of the PKA-mediated isoform shift and/or change in line width. Additionally, there is an increase in the total number of resonance peaks, which is likely due to the presence of more than one stable conformation. This would explain the double band observed for the PKA-mediated isoform in the SDS-PAGE (Fig. 1C).
The phosphorylation-mediated effects could be the result of stabilizing an existing dynamic ␣-helix or disordered residues now transitioning between ␣-helical and random coiled structures (i.e. changing the conformational preference). Therefore, 15 N NOESY-HSQC experiments were performed to determine which of these two possibilities account for the changes in Cx43CT secondary structure. Chemical shift analysis of the WT TM4-Cx43CT predicted seven ␣-helical regions along the C-terminal domain (53). However, many of the expected NOEs were not present in the 15 N NOESY-HSQC spectrum indicating the ␣-helical regions are flexible and in exchange between a structured and unstructured conformation (Fig. 6, left panel). The CK1 isoform was selected as a representative phosphomimetic isoform to study because the ␣-helical content of the CK1 isoform is greater than WT and increases upon acidification. Additionally, unlike the MAPK isoform, all of the CK1 phosphorylation sites are Asp-substituted. The CK1 isoform 15 N NOESY-HSQC spectrum was collected using the identical NMR parameters, protein concentrations, and solution conditions as WT. Upon visualization of the CK1 isoform data, NOE peaks had broad line widths and weak signal intensity, when compared with the WT TM4-Cx43CT (Fig. 6, middle panel). A ϫ3 magnification of the CK1 spectrum was required to obtain the same signal intensity as the WT spectrum (Fig. 6, right  panel). If the change in CT phosphorylation state had stabilized ␣-helical structure, then additional peaks would have been detected in the CK1 isoform 15 N NOESY-HSQC. Even though the CD data indicate an increase in the overall ␣-helical content, the NOESY data suggest the TM4-Cx43CT CK1 isoform is in intermediate exchange. This observation is consistent with phosphorylation altering the fluctuation of Cx43CT residues between random coiled and ␣-helical structures.
NMR Relaxation Studies of the Membrane-tethered Cx43CT Domain-All proteins undergo time-dependent molecular motions that affect NMR relaxation rates. Molecular motions, such as protein folding, side chain rotation, molecular tumbling, and bond vibrations occur over a broad time scale (1-10 Ϫ15 s (54)). The relaxation parameters R 1 , R 2 , and 1 H-15 N NOE describe site-specific chemical exchange processes (e.g. a nucleus exchanging between environments) occurring on the

Phosphorylation-mediated Cx43CT Structural Changes
picosecond to nanosecond time scale (12,55). For a disordered protein, the folding energy landscape is relatively flat compared with a structured protein (i.e. no energetically favored conformation) resulting in greater protein picosecond to nanosecond dynamics. One of the difficulties of working with an IDP is low NMR peak dispersion, which often results in peak overlap and limits the number of peaks suitable for measurement. NMR relaxation data were collected for WT TM4-Cx43CT and each of the phospho-mimetic isoforms to determine the effect(s) of phosphorylation on the backbone dynamics of the Cx43CT domain. Relaxation measurements were obtained for 80 of the 125 assigned WT TM4-Cx43CT non-proline peaks. For the phospho-mimetic isoforms, fewer peaks were used due to additional overlap, disappearance, and/or ambiguous identification (Table 3). In total, R 1 and R 2 relaxation measurements were collected for 67 peaks of the MAPK and Src isoforms, 65 for the cdc2 isoform, 54 for the PKC isoform, 53 for the CK1 isoform, and 51 for the PKA-mediated isoform. Of note, 15 N HSQC assignment of the phospho-mimetic isoforms was determined by comparing the chemical shifts of the each of the isoforms to the previously determined WT chemical shifts (53). Unfortunately, the most interesting phospho-mimetic isoform peaks, which are those that have significant changes in chemical shift, could not be matched to a corresponding WT peak (e.g. CK1 isoform, only 4 of 23 peaks with significant changes could be unambiguously identified). Table 3 provides the phospho-mimetic isoform peaks of interest that were suitable for relaxation analysis.
The WT TM4-Cx43CT R 1 measurements show a relatively flat profile with a steep decline at the C terminus, characteristic of a disordered protein (Fig. 7A) ; Fig. 7B) also indicate fast nanosecond to picosecond time scale dynamics consistent with the TM4-Cx43CT being an IDP. R 1 and R 2 values of structured protein are ϳ2 and 4 times greater than the corresponding values of a disordered protein of similar size, respectively (56,57). 1 H-15 N NOE data are also consistent with fast dynamics (Fig. 7C, mean ϭ Ϫ0.15).
Several residues along the Cx43CT have R 1 and R 2 values Ͼ1 standard deviation (S.D.) above the mean (R 1 , Ala-311, Ala-315-Gln-317, Ala-367, and Ala-371; R 2 , Tyr-286, Gly-291, Asn-295, Tyr-301 Ala-357, and Ile-358), suggesting they are less flexible. These residues overlap with ␣-helical regions predicted from the chemical shifts of the WT TM4-Cx43CT domain (53), which encompasses the two ␣-helical regions of the soluble Cx43CT (7). However, the R 1 and R 2 values of the ␣-helical regions are still within the range of an IDP and are not suggestive of rigid secondary structure. R 1 and R 2 values Ն1 S.D. below the mean were observed from the largest stretch of consecutive disordered residues (Ala-323-Asp-339) from the WT TM4-Cx43CT predicted structure, indicating increased flexibility. Another observation was that many of the residues differ between the R 1 and R 2 data. One possible explanation is the contribution of microsecond to millisecond chemical exchange to R 2 (58).
A common approach used in analysis of NMR R 1 and R 2 relaxation data is to consider the R 2 /R 1 ratio and the R 2 R 1 product (59,60). If the sequence-specific rates of protein motion are isotropic, as is typical for a well folded protein, then the R 2 /R 1 ratio for each residue will be similar. R 2 /R 1 values that differ from the mean indicate anisotropic motions on the picosecond to nanosecond time scale. One caveat to this analysis is that R 2 relaxation rates often have a large contribution by exchange processes on the microsecond to millisecond time scale (59). The R 2 R 1 product mitigates the contribution of anisotropic motions on the picosecond to nanosecond time scale such that deviations Ͼ1 S.D. above the mean are indicative of chemical exchange on the microsecond to millisecond time scale (60). Residues within the regions Ser-282-Tyr-286, Arg-293-Cys-298, and Leu-354 -Cys-358 have R 2 /R 1 ratios Ͼ1 S.D. above the mean indicating less picosecond to nanosecond motion (Fig.  7D). However, the WT TM4-Cx43CT R 2 /R 1 ratio profile is nearly identical to the R 2 R 1 product profile (Fig. 7E). This suggests deviations Ͼ1 S.D. above the mean are due to chemical exchange on the microsecond to millisecond time scale rather than anisotropic motion.
Also examined were the relaxation profiles of the phosphomimetic isoforms. The R 2 /R 1 ratios for each of the phosphomimetic isoforms (data not shown) are nearly identical to their respective R 2 R 1 product profiles (Fig. 8), suggesting values Ͼ1 S.D. above the mean are also the result of chemical exchange processes. With the exception of the PKA isoform, all of the phospho-mimetic isoforms have two major areas with R 2 R 1 values Ͼ1 S.D. above the mean. The areas overlap with WT regions Ser-282-Tyr-286 and Leu-354 -Cys-358 but encompass slightly different residues. As an example, two regions of the CK1 isoform have increased R 2 R 1 that include more residues (Fig. 8A, Ala-276 -Thr-290 and Ser-357-Ser-362). The CK1 isoform also has residues Ͼ1 S.D. below the mean (Lys-346 -Glu-352) indicating increased picosecond to nanosecond motion in this region. The PKA-mediated isoform R 2 R 1 profile (Fig. 8F), which is the least similar to the WT profile, indicates four regions with R 2 R 1 values greater than the mean (Gly-285-Asp-292, Asn-294 -Arg-299, Ser-314 -Gln-317, and Ala-344 -Gly-350).
Finally, the 1 H-15 N NOE values for all phospho-mimetic isoforms were consistent with fast amide N-H bond vector dynamics (data not shown). Similar to the R1 and R2 data, there are slight differences in the phospho-mimetic isoform NOE profiles (e.g. CK1 isoform, residue Phe-337 shows slower dynamics than WT). Altogether, the NMR data suggest that even for the phospho-mimetic isoforms determined by CD to have the greatest ␣-helical content (e.g. MAPK, CK1, cdc2), the increased ␣-helical content is not manifesting in rigid structure. Rather more random coiled residues have now transitioned to fluctuating between random coiled and ␣-helical structure.

TABLE 2 Effects of pH and phosphorylation on TM4-Cx43CT secondary structure
Helical content was determined using DichroWeb Provencher and Glöckner method analysis using the SP175 reference. The ϭ , 1, and 2 indicate ␣-helical content relative to WT. mp means mediated phosphorylation.

Utility of the TM4-Cx43CT Construct for Investigating Phosphorylation-induced Changes in Structure-The soluble
Cx43CT has proven to be a useful model for studying the structure-function mechanisms regulating gap junctions. However, several results indicate the soluble Cx43CT may not be the best model system. For example, the electron crystallography struc-ture of the CT-truncated Cx43 mutant (terminates at Thr-263) suggested that region Ser-255-Thr-263 adopts an ␣-helical structure (4,61). NMR data for the soluble CT identified that peaks were weak from this region, suggesting these residues are in exchange between unstructured and ␣-helical conformations (62). Additionally, not all of the expected NOEs were observed in the two ␣-helical regions of the soluble CT struc-  15 N HSQC spectra of each of the phospho-mimetic isoforms (red; A-F). Phospho-mimetic isoform chemical shifts that were different from WT are mapped to the Cx43CT sequences given in G. Blue letters, residues with minor changes in chemical shift (peaks are still slightly overlapped); red letters, residues with large changes in chemical shift or peaks that disappeared; green, underlined* letters, sites of Asp substitution; black lines indicate the seven regions previously predicted to be ␣-helical. Phospho-mimetic isoforms are named based on kinase responsible for the endogenous phosphorylation. AUGUST 23, 2013 • VOLUME 288 • NUMBER 34 ture. The differences can be attributed to the constraints imposed when attached to the membrane. Therefore, expression, purification, and NMR solution conditions were optimized for a more native-like construct, the TM4-Cx43CT solubilized in detergent micelles (34,63). CD data indicated the TM4-Cx43CT has more ␣-helical content than can be attributed solely to the addition of the individual TM4 and C-terminal domains. The TM4-Cx43CT also undergoes a pronounced increase in ␣-helical content upon solvent acidification (i.e. ischemic conditions), unlike the soluble Cx43CT (Fig. 3C). The structural responsiveness to changes in pH and phosphorylation indicates the TM4-Cx43CT is a better model than the soluble Cx43CT for investigating structure-based mechanisms of Cx43 regulation.

Phosphorylation-mediated Cx43CT Structural Changes
Previously, we reported the backbone assignments and predicted secondary structure of the TM4-Cx43CT (53). The total amount of ␣-helical structure observed from CD data of the TM4-Cx43CT construct (46%) is consistent with the combined content from the seven predicted ␣-helical regions of the C-terminal domain and ␣-helical TM4 domain (34). However, not all of the NOEs were apparent in these ␣-helical regions suggesting they are flexible and in exchange between random coiled and ␣-helical conformations. Here, comparison of the TM4-Cx43CT CD spectra with a construct containing only the TM4 portion (Asp-197-Val-240) confirmed pH-dependent structural changes occurred only in the C-terminal domain; the ␣-helical content of the TM4 is unaffected by changes in pH. Therefore, we expect phosphorylation-mediated changes in secondary structure to also only occur in the C-terminal domain.
Phosphorylation Alters the Secondary Structure and Dynamics of the Cx43CT-One difficulty in understanding how Ser phosphorylation regulates protein function is the inability to predict whether phosphorylation will act to stabilize or destabilize secondary structure (31). For a disordered protein, the large negative charge of a phosphate could alter backbone hydration facilitating formation of intramolecular hydrogen bonds that stabilize secondary structure. The phosphate could also provide electrostatic (or steric) interactions that affect the formation of secondary structures. For the soluble Cx43CT domain, the effect on secondary structure correlated with the location of the phosphorylation. The PKA-mediated and PKC isoforms, which are modified at the C terminus of the Cx43CT, are both more ␣-helical than the WT TM4-Cx43CT at pH 7.5 but show little to no change in secondary structure upon acidification. The CK1 and MAPK isoforms, which have phosphorylation sites centrally located, have greater ␣-helical content than WT at pH 7.5 and become more ␣-helical upon acidification. Conversely, cdc2, which is closer to the N terminus, has more ␣-helical content at both pH 5.8 and 7.5 than WT but undergoes little to no change in ␣-helical content upon acidification. Although the use of an Asp substitution to mimic Tyr  phosphorylation is not optimal, Asp was chosen over Glu for the Src isoform to maintain consistency with the other phospho-mimetic constructs. The Src isoform is less ␣-helical than the WT at pH 5.8, as would be expected if a true phosphate were present; Tyr phosphorylation, unlike Ser phosphorylation, generally disrupts ␣-helical structure (64). Furthermore, the Aspsubstituted Tyr residues in the Src isoform elicited an opposite effect on ␣-helical content compared with Asp-substituted Ser residues in a similar location (i.e. cdc2 isoform). The 15 N HSQC spectra of the phospho-mimetic isoforms, except the PKC isoform, show the effects of phosphorylation to propagate along the C-terminal domain. The line widths of the peaks in the 15 N HSQC spectra suggest residues are in intermediate exchange between more than one conformation. Collectively, the CD and 15 N HSQC NMR data indicate phosphorylation alters the conformational preference of the C-terminal domain between random coiled and ␣-helical structure. Conformational exchange is consistent with the WT TM4-Cx43CT and CK1 isoform 15 N NOESY-HSQC data. Most of the expected NOEs were not present for residues in the ␣-helical regions, leav-ing a limited number of NMR peaks suitable for relaxation analysis. The NMR relaxation data indicate the ␣-helical structure is dynamic and in exchange with a disordered conformation. Although slight differences were observed between WT and the phospho-mimetic isoforms, the most interesting peaks to analyze would be those with significant changes (Table 3; Fig. 5G, red residues) as these residues would have the greatest variation in flexibility and/or secondary structure. NMR is similar to CD as both are averages of the ensemble of protein conformations in solution. The increased ␣-helical content in some of the TM4-Cx43CT phospho-mimetic isoforms is not due to stabilization of rigid secondary structure; phosphorylation shifts the conformational preference of random coiled residues to form dynamic ␣-helical structure.
An NMR study using a nonphosphorylated and phosphorylated Cx43CT peptide (Lys-234 -Asp-259) determined Src phosphorylation at Tyr-247 regulates the interaction between Cx43 and tubulin (28). The nonphosphorylated peptide adopted an ␣-helical conformation upon interaction with tubulin. However, phosphorylation at Tyr-247 had a negative impact on tubulin recognition. In the ␣-helical conformation, Tyr-247 would be on the hydrophobic face of the peptide involved in the tubulin interaction. Although the phosphate at Tyr-247 may act to directly inhibit tubulin binding, the data presented here suggest the phosphate may also inhibit this interaction by altering the conformational preference of these CT residues to a disordered state. This would destabilize the hydrophobic face necessary for the interaction. Changing the conformational preference of an IDP can fine-tune the thermodynamic favorability of protein partner interactions (14,18). Backbone dynamics on the picosecond to nanosecond time scale have been implicated in determining ligand specificity and affinity and have been suggested as a regulatory mechanism for IDPs in general (17,54).
How Does Phosphorylation Regulate GJIC?-Several hypotheses, which are not mutually exclusive, could explain how phosphorylation regulates GJIC at the plasma membrane as follows: 1) the negative charge of the phosphate could affect the permeability of ions through the pore via electrostatic attraction or repulsion; 2) the phosphate could stabilize favorable electrostatic interactions between the Cx43CT and the positively charged residues of the Cx43CL leading to blockage of the pore, FIGURE 8. Effects of phosphorylation on TM4-Cx43CT R 2 R 1 backbone dynamics. 15 N relaxation parameters (R 1 and R 2 ) for the phospho-mimetic isoforms were measured at 600 MHz. The R 2 R 1 values are given for each of the isoforms in A-F. Residues with chemical shift overlap or ambiguous assignments (including * sites of substitution for the cdc2 and Src isoforms) were excluded from the analysis. Vertical lines indicate sites of Cx43CT phosphorylation; mean values (solid line) and standard deviations (horizontal dotted lines) are provided. and 3) phosphorylation of the CT could directly alter the binding affinities of proteins involved in Cx43 regulation (e.g. Nedd4 (30), 14-3-3 (65), and ZO-1 (27)). Moreover, the results presented here indicate phosphorylation mediates changes in the secondary structure and flexibility of the Cx43CT domain, providing support for additional mechanisms by which phosphorylation may regulate GJIC. Additionally, CT structural changes induced by phosphorylation could 4) alter channel permeability by indirectly modulating the orientation of the transmembrane ␣-helices thereby influencing pore size and/or 5) indirectly alter the binding affinity of molecular partners involved in regulating GJIC by altering the conformational preference of the Cx43CT domain.
Conclusion-This study shows the first direct evidence of phosphorylation-mediated structural changes of Cx43. Additionally these structural changes are not due to stabilization of rigid helical structures but rather are due to a change in the conformational preference of the Cx43CT. This is a regulatory mechanism previously demonstrated for other intrinsically disordered proteins but is a novel mechanism for Cx43 regulation. Importantly, qualitative analysis of both CD and NMR relaxation was needed to identify this novel mechanism. Future studies will work toward more detailed analysis of each of the phospho-mimetic isoforms in terms of how changes in the conformational preference mediate specific protein partner interactions involved in regulating Cx43 gap junctional communication.