C-type inactivation of potassium channels fine-tunes the electrical signaling in excitable cells through an internal timing mechanism that is mediated by a hydrogen bond network in the channels' selectively filter. Previously, we used nonsense suppression to highlight the role of the conserved Trp434-Asp447 indole hydrogen bond in Shaker potassium channels with a non-hydrogen bonding homologue of tryptophan, Ind (Pless et al., 2013). Here, molecular dynamics simulations indicate that the Trp434Ind hydrogen bonding partner, Asp447, unexpectedly 'flips out' towards the extracellular environment, allowing water to penetrate the space behind the selectivity filter while simultaneously reducing the local negative electrostatic charge. Additionally, a protein engineering approach is presented whereby split intein sequences are flanked by endoplasmic reticulum retention/retrieval motifs (ERret) are incorporated into the N- or C- termini of Shaker monomers or within sodium channels two-domain fragments. This system enabled stoichiometric control of Shaker monomers and the encoding of multiple amino acids within a channel tetramer.https://doi.org/10.7554/eLife.18976.001
C-type inactivation, also termed ‘slow-inactivation’ due to its delayed kinetics, in voltage-gated potassium channels describes a phenomenon by which ionic conductance through the channel pore is decreased in the presence of continued depolarizing stimulus (Choi et al., 1991; Hoshi et al., 1991). In a physiological setting, this mechanism represents a form of activity-dependent macromolecular memory that can alter action-potential firing rates (Aldrich et al., 1979; Roeper et al., 1997). The basis for C-type inactivation has been investigated with functional, structural and computational approaches, which have produced a consensus mechanism whereby channels become non-conducting due to local perturbations in the selectivity filter (Yellen et al., 1994; Baukrowitz and Yellen, 1995; Kiss and Korn, 1998; Perozo et al., 1998; Roux and MacKinnon, 1999; Bernèche and Roux, 2001; Cordero-Morales et al., 2007; Panyi and Deutsch, 2007; Cuello et al., 2010; Ostmeyer et al., 2013; Thomson et al., 2014). However, the nature of the conformational change and role of specific residues remain an area of active experimentation and debate (Devaraneni et al., 2013; Hoshi and Armstrong, 2013).
We previously identified an intra- and inter-subunit hydrogen bond network that controls the rate of entry into the non-conducting pore conformation (Pless et al., 2013). Our study employed nonsense suppression in the Xenopus laevis oocyte expression system to genetically encode a variety of synthetic noncanonical amino acids (ncAAs). In the context of the Xenopus laevis oocyte, the cRNA of the target protein contains an introduced stop codon (often amber, TAG) that is subsequently suppressed by a co-injected orthogonal tRNA that has been misacylated in vitro with the ncAA. This approach was used to encode a tryptophan homologue, 2-amino-3-indol-1-yl-propionic acid (referred to as Ind) that fails to participate in hydrogen bonding at position 434, a site previously identified as being important for channel inactivation. Our original efforts failed to determine the functional role of individual hydrogen bonds with Ind because subunit concatenation in the tetrameric construct resulted in a significant loss of expressing channels, an outcome that prevented the application of nonsense suppression (Pless et al., 2013). Reduced channel expression is a common challenge of the nonsense suppression method where amino acid incorporation varies greatly from position to position within a single reading frame and between channel families. Furthermore, given the central role of the Trp434-Asp447 hydrogen bond in the mechanism of C-type inactivation, it was not possible to simultaneously encode Ind within each of the four subunits of the channel because this resulted in a non-conducting phenotype which lacked the relatively high expression levels needed to measure gating currents which originate from the movement of the channel voltage-sensor in the transmembrane electric field (Perozo et al., 1993). Thus, we concluded that while one can control the stoichiometry of individual subunit composition in a concatenated construct, their low expression impeded the use of nonsense suppression in combination with this approach. In our original publication, this shortcoming was partially overcome by combining WT cRNA with Ind-tRNA and Trp434TAG cRNA within an individual oocyte. This experimental approach resulted in predicted accelerated, but complex, inactivation kinetics consistent with a mixed population of channels with 1, 2 or 3 Ind434-containing monomers (see Figure 2F, Pless et al., 2013, eLife). In the present study, we examined the mechanism of Ind-induced inactivation by generating a new approach to control subunit stoichiometry to test individual hydrogen bond contributions in Shaker channel C-type inactivation and more generally, for the encoding of multiple ncAA within a single channel or an ion channel complex in the application of nonsense suppression.
Targeted replacement of Trp434 with the isosteric ncAA Ind was used previously to highlight the general importance of the indole hydrogen bond of this highly conserved Trp residue. The Ind amino acid has seen limited use in structure-function studies for hydrogen bond testing (Lacroix et al., 2012; Pless et al., 2013, 2014; Kim et al., 2015; Zhang et al., 2015) and none thus far in structural or computational studies that would inform on its tolerance once encoded within a protein. For instance, it is not known if the Ind substitution would specifically limit hydrogen bonding at the indole nitrogen atom (the purpose for which it was designed) or could it also produce a new side-chain orientation that itself alters protein function. A computational approach based on molecular dynamics (MD) simulations was therefore employed to examine the local dynamics of the Ind side-chain and its relation to the intra-subunit hydrogen bond accepting side-chain Asp447. These MD simulations reveal key conformational and dynamic differences between wild type (WT) and Trp434Ind of the Shaker channel selectivity filter in the conductive state (Long et al., 2005) (Figure 1A). In WT Shaker, a hydrogen bond is formed between Trp434 and Asp447, restraining the two residues close to each other during a large fraction of a 500 ns trajectory (Figure 1C). In the absence of the Asp447-Trp434 hydrogen bond in the Trp434Ind Shaker, there is an increased distance between these two residues in all of the four subunits (Figure 1C), which is associated with a flipping of the Asp447 side-chain toward the extracellular bulk solution (Figure 1B). See also Supplemental Videos 1 and 2 for wildtype W434 and W434Ind channels which depict 150 ns vignettes of the full molecular dynamics simulations.
The breaking of the Asp447-Trp434 hydrogen bond in the Trp434Ind mutant has two major effects in the vicinity of the selectivity filter. First, the formation of a hydrogen bond between Asp447-Trp434 in the wild type blocks the exchange of water molecules between extracellular bulk and the water binding pocket behind the selectivity filter, as shown in the water occupancy map in Figure 1D. In the absence of this hydrogen bond in Trp434Ind, the increased spatial separation of Asp447 and Trp434Ind promotes the formation of a continuous connection between the bulk solution and the water molecules in the buried pocket. The easier accessibility of the buried water-binding pocket probably accelerates C-type inactivation of Shaker. Second, the flipping of Asp447 leads to an additional effect to the selectivity filter. Considering Asp447 in each subunit carries a negative charge, the reorientation of the Asp447 side-chain from facing inside to outside (Figure 1C) perturbs the local electrostatic potential at the outer end of the selectivity filter, altering the binding affinity of potassium ions in multiple binding sites, and is likely to shift the relative stability between conductive and nonconductive states.
An engineered protein approach was designed to circumvent issues raised with expressing Shaker monomer concatemers: complementary split inteins (Wu et al., 1998; Muralidharan and Muir, 2006), were incorporated into the N- or C- termini of Shaker monomers (Figure 2A). Split inteins drive a post-translational event involving precise excision of the intein sequence and formation of a normal peptide bond between the flanking sequences (N- and C-exteins) (Saleh and Perler, 2006). We reasoned that Shaker monomers would form split intein-dependent dimers, and Shaker channels would consist of dimers of dimers. In order to exclude nonconforming Shaker monomers from the surface membrane, we encoded endoplasmic reticulum retention/retrieval (ERret) motifs which flank the split inteins (Figure 2A) that would be removed following the intein ligation process (Figure 2B). We used the Kir6.2 ERret motif because it works in a variety of eukaryotic cell types, has strong retention properties, and functions independently from its position within a protein (Zerangue et al., 1999). Oocytes injected with only one of the monomer cRNAs did not show measurable voltage-dependent K+ selective currents (Figure 2A, lower panels). This is presumably due to the retention of the monomers within the ER as a result of the appended ERret motifs. In contrast, co-expression of Shaker monomers with complementary split inteins resulted in robust outward voltage-dependent potassium currents due to the nature of the split intein protein ligation mechanism. The formation of tetrameric dimer-of-dimer complexes with no ERret motif allows trafficking of Shaker channels from the ER to the plasma membrane. The level of expressed potassium current was stable over a large range of injected Shaker monomer cRNA concentrations (Figure 2B) suggesting a limiting cellular component in channel assembly and/or trafficking. This method allows for control of subunit stoichiometry and the possibility for encoding multiple Ind side-chains within a single channel tetramer by placing the introduced TAG mutation within one (or both) of the two complementary fragments that form a dimer.
In order to test cooperativity of the hydrogen bond at position 434, the ERret-intein approach was used to express Shaker channels with two Ind residues at position 434 per channel (Figure 3). In experiments where both WT and Trp434TAG complimentary ERret-intein cRNAs and Trp-tRNA were injected, robust WT-like currents were observed, demonstrating successful Trp incorporation (Figure 3A, left). Conversely, when Ind-tRNA was expressed within the C-terminal Shaker monomer of each dimer, C-type inactivation was significantly accelerated, consistent with the breakage of two of the four hydrogen bonds around the central permeation pathway (Figure 3A, right) while voltage dependence of activation was not altered (Figure 3B, Table 1). At +20 mV, WT Shaker channels had an inactivation time constant of 3900 +/- 0.2 milliseconds (first-order exponential) compared to 529 +/-141 and 142 +/-14 milliseconds when one or two hydrogen bonds were broken through Ind incorporation, respectively (Figure 3C).
In order to examine the general utility of the ERret-intein sequences, they were constructed within two-domain hemi-channels of the skeletal muscle isoform of the voltage-gated sodium channel NaV1.4. Previously, it has been shown that both sodium and calcium channels are amenable to functional complementation via expression of ‘split’ hemi-channels consisting of DI-DIII and DIV or DI-II and DIII-DIV, respectively (Stühmer et al., 1989; Ahern et al., 2001). This result suggests that sodium and calcium channels need not be tethered on a single polypeptide to form a functional tetrameric channel. We therefore reasoned that co-assembly of sodium channel domains I-II and III-IV would likely bring ERret-intein sequences on the C-terminus of DI-II and the N-terminus of DIII-IV in close proximity and promote intein dependent cleavage of the ERret sequences and formation of a peptide bond between domains II and III (Figure 4A). To this end, cRNA for Nav1.4 DI-II and DIII-IV ERret-intein containing hemi-channel fragments were expressed in Xenopus oocyte and subsequent currents were recorded by TEVC. As with the Shaker intein constructs, expression of either individual fragment did not result in expressed ionic currents (Figure 4D; triangles and squares). However, co-expression of the split NaV1.4 halves resulted in robust inward voltage-gated sodium currents (Figure 4C– D, diamonds) similar to WT full-length channels (Figure 4D, circles). Notably, full-length sodium channels produced from the ERret-inteins had biophysical properties that were indistinguishable from WT sodium channels, p>0.01, Table 1. To examine the general versatility of the approach, we expressed the NaV1.4hemi-channel construct cDNA by transient transfection into HEK293 cells and examined the ligation efficiency by Western blot using C-terminal HA epitopes that are on both N- and C-channel fragments (Figure 4B). In the case of the DI-II-HA construct, this epitope is cleaved out during the in-cell ligation reaction but the full-length channel can still be identified by the remaining C-terminal HA epitope on the DIII-IV protein (see methods for protein sequence and construct design). The Western blot data show that each hemi-channel is robustly expressed, producing bands of ~150 kDa (Figure 4B, lane 2) and ~110 kDa (Figure 4B, lane 3) for the DI-II-HA and DIII-DIV-HA truncated constructs, respectively. Co-expression of the channel halves resulted in the formation of a ~250 kDa band (lanes 4 and 5) similar to WT NaV1.4 full-length channels (lane 1). The efficiency of protein ligation was similar to that of intein ligated Ca2+ channel hemi-channel fragments where also, interestingly, the fragment corresponding to the DIII-IV construct was limiting (Subramanyam et al., 2013). To rule out the possibility that the ligation reaction occurred during sample preparation, i.e after the fragments had been solubilized, cell lysates from NaV1.4 DI-II or DIII-IV expressing HEK cells were mixed and then incubated for 10, 30 or 60 min, Figure 4, Figure 1—figure supplement 1. No higher band corresponding with the full-length NaV1.4 was observed at these time points. Thus, the biochemical data demonstrate that the NaV1.4 DI-II and DIII-IV ERret-intein constructs are competent in HEK cells where they are ligated to produce full-length sodium channels.
Finally, we aimed to see if the split-intein ERret strategy could assist with nonsense suppression efficiency when attempting to simultaneously suppress two stop sites within a single cRNA transcript. The data in Figure 3 with Shaker ERret-inteins show that this approach can be used to produce tetrameric K+ channel complexes containing two ncAA, but these are within the dimer of dimers K+ channel structural context (i.e. not within the same polypeptide). In contrast, double suppression within the same reading frame poses a particular challenge. Specifically, the ability of an acylated tRNA in the oocyte expression system to suppress stop sites has variable efficiency that is compounded when two stop sites are present, resulting in the expression of full-length channels that are far below that of WT channels, and often immeasurable using electrophysiological methods. Thus, suppression at the second site down-stream from the first compounds the overall efficiency of producing a full-length peptide. We hypothesized that if the two sites are suppressed independently, within the individual Nav1.4 I-II and III-IV hemi-channels for instance, and then subsequently ligated via the ERret-intein sequences that the reduced expression would not be compounded but only limited by the site with the lower suppression efficiency. To test this possibility directly, stop (TAG) sites were engineered into the DI-II construct at Tyr401 near the outer vestibule of the selectivity filter and within the DIII-IV hemi-channel at Phe1304, a site within the DIII-IV linker. When expressed as hemi-channels, the simultaneous dual suppression with Phe-tRNA at Tyr401 and Phe1304 yielded macroscopic sodium currents with gating similar to wildtype NaV1.4 demonstrating that the split ERret-intein method enables dual suppression, Table 1, Figure 5—figure supplement 1. The incorporation was specific to encoding of a Phe residue, because co-injecting the NaV1.4 TAG cRNA with a tRNA lacking an acylated amino acid failed to produce sodium current, Figure 5.
Taken together, the results presented here advance the specific molecular understanding of C-type inactivation and introduces a novel approach which employs the use ncAAs with split intein protein engineering. MD simulations indicate that the Trp434Ind substitution alters the orientation of the hydrogen bond partner, Asp447. Notably, the ‘flipped out’ orientation of Asp447 occurs in WT channels, but does so briefly and returns to the hydrogen-bonded state. Thus, the simulations suggest that the Ind substitution increases the likelihood of the ‘flipped out’ Asp447 conformation, which disrupts the local environment of the selectivity filter. This orientation of Asp447 allows water to fill the space behind the selectivity filter and simultaneously alters the local electrostatics properties, two outcomes that likely impact the occupancy of potassium ions in the filter. One possibility is that the reorientation of Asp447 may represent an event in the mechanism of C-type inactivation as the channel transitions from the open/conductive (O/C) to a C-type inactivated conformation (O/I), see Panyi and Deutsch, 2007. Further, the breakage of two versus one of the four Asp447-Trp434 hydrogen bonds further accelerates C-type inactivation, possibly due to the loss of an additional hydrogen bond and ultimately results in a concomitant destabilization of the open pore. However, the data do not inform on the details of movements required to break the hydrogen bonds or the nature of the final non-conducting state, an event that occurs downstream from hydrogen bond breakage.
We previously attempted to understand domain-specific contributions of hydrogen bonding and C-type inactivation with a Shaker concatemer with all four subunits tethered into a single peptide (Yang et al., 1997). However, while this construct tolerated a single Trp434Phe mutation by site-directed mutagenesis (Pless et. al, 2013, eLife) encoding an ncAA was not feasible as the concatemer displayed a significant decrease in overall expression, a common effect of protein concatenation. Here we describe a protein engineering approach for altering side-chain chemistry prior to concatenation that produces macroscopic potassium currents from Shaker channels containing two Ind amino acids. Split intein sequences have been used to ligate a variety of protein fragments (Shah and Muir, 2014) and the ERret motif used in this study has been shown to work in several types of eukaryotic cells including yeast (Zerangue et al., 1999). Therefore, this approach could be used to examine other proteins in other cellular environments for studying the effects of dominant negative mutations, heteromers, etc. Western blot data of expressed Nav1.4 ERret-intein hemi-channels demonstrates the effectiveness of the approach outside of the Xenopus laevis oocyte expression system and also suggests the ERret-intein can be transferred to other channel types aside from potassium channels. Encoding of multiple ncAA within an individual protein complex remains a significant challenge. Where success is reported, it stems primarily from improving the efficiency of the suppression process and tRNA delivery (Neumann et al., 2010; Wan et al., 2010; Johnson et al., 2011; Chatterjee et al., 2013; Xiao et al., 2013; Lammers et al., 2014). In contrast, the approach described here changes the nature of the dual suppression problem by rescuing each transcript independently prior to protein ligation thus abrogating the loss of expression due to the compounded efficiency of multiple suppression events. Further, the sodium channel data obtained with dual TAG sites suggests that the approach will be useful for suppression of multiple stop sites, a possibility that will be valuable for encoded FRET, once such probes become available.
Initial coordinates for the molecular dynamic simulations were modeled based on the crystal structure of Kv1.2/Kv2.1 Voltage-gated potassium channel chimera (PDB ID:2R9R). Considering the high sequence identity (87%) of the pore domain between Kv1.2 and Shaker, several non-conserved residues were mutated from Kv1.2 to Shaker, and all of the conserved residue side-chains still maintain the same conformation as the crystal structure. Residues were assigned their standard protonation state at pH7.0. Only the pore domain of the channel was embedded in a bilayer of POPC/POPG (3:1) lipids and solvated in 0.2 M KCl using the web service CHARMM-GUI [1-3] (total number of atoms ~ 44,500). The CHARMM36 force field (MacKerell et al., 1998, 2004; Klauda et al., 2010; Best et al., 2012) was used for proteins and lipids, and explicit water was described with TIP3P model (Jorgensen et al., 1983). For mutant Trp434Ind, the unnatural amino acid was substituted using CHARMM-GUI PDB Reader (Jo et al., 2014), and the force field was obtained from CHARMM General force field (CGenFF) (Vanommeslaeghe et al., 2010). All-atom simulations were performed in constant NPT conditions at a temperature of 310 K and a constant pressure of 1.0 atm, after initial minimization and equilibrations. All the simulations were performed under periodic boundary conditions with a time step of 2 fs. Throughout the simulations, bond distances involving hydrogen atoms were fixed using the SHAKE algorithm (Ryckaert et al., 1977). Short-range non-bonded interactions were calculated using a cutoff distance of 12 Å, and long-range electrostatic interactions were calculated using the particle mesh Ewald (PME) method (Darden et al., 1993). Molecular dynamic simulations were carried out using NAMD 2.10 (Phillips et al., 2005) or OpenMM 6.2 (Eastman et al., 2013). For both wild type and mutant Trp434Ind, MD simulations were performed for 500 ns.
Shaker IR (Inactivation Removed by deletion of amino acids 6–46) cDNA in pBSTA was used as the parent clone for all cDNA plasmids. Plasmids containing the DnaE split inteins (DnaE-n and DnaE-c) were a generous gift from Henry M. Colecraft (Columbia University, New York, NY) and were previoulsy published in Subramanyam et al., 2013. Two ERret-intein Shaker plasmids were generated. For the first cDNA plasmid, the stop codon was removed and a glycine rich linker (GSAGSAAGSGEFAS) was added to the c-terminus of the Shaker coding region using standard PCR. The DnaE-n sequence (bold) followed by a flexible linker (PINQPANTHERPR), Kir6.2 ER retention/retrieval motif (italics) and canonical KKXX ER retention/retrieval motif (underlined) was generated and ligated into the afformentioned plasmid:
(CLSYETEILTVEYGLLPIGKIVEKRIECTVYSVDNNGNIYTQPVAQWHDRGEQEVFEYCLEDGSLIRATKDHKFMTVDGQMLPIDEIFERELDLMRVDNLPNPINQPANTHERPRLLDALTLASSRGPLRKRSVAVAKAKPKFSISPDSLSPRKKFQstop). For the companion Shaker cDNA plasmid, the start codon was removed and a flexible linker (ARSMEAS) was added to the n-terminus of the Shaker coding region using standard PCR. The Kir6.2 ER retention/retrieval motif (italics) and DnaE-c sequence (bold) was generated and ligated 5’ inframe with the flexible linker using standard molecular biology techniques:
(MLLDALTLASSRGPLRKRSVAVAKAKPKFSISPDSLSPRIKIATRKYLGKQNVYDIGVERDHNFALKNGFIASNCFN). Two ERret-intein rNav1.4 plasmids were generated, Nav1.4 Domain I-II- Intein ERret and Nav1.4 ERret Intein III-IV in pcDNA3.1(+) Zeocin (Promega, Madison WI). For the front half Nav1.4 Domain I-II- Intein ERret plasmid, the DnaE-n intein sequence (bold) was ligated inframe following rNav1.4 Glu985. The intein sequence was followed by tandem human influenza hemagglutinin (HA) epitope sequences (underlined) that were used as both flexible linkers and for Western blot analysis. The HA linker was followed by the Kir6.2 ER retention/retrieval motif (italics) and canonical KKXX ER retention/retrieval motif (italics underlined) (CLSYETEILTVEYGLLPIGKIVEKRIECTVYSVDNNGNIYTQPVAQWHDRGEQEVFEYCLEDGSLIRATKDHKFMTVDGQMLPIDEIFERELDLMRVDNLPNYPYDVPDYAYPYDVPDYLLDALTLASSRGPLRKRSVAVAKAKPKFSISPDSLSPRKKFQstop). For the back half Nav1.4 ERret Intein III-IV plasmid, the ERret Dna-c intein (bold) and tandem HA epitope (underlined) was ligated 5 ‘and inframe of rNav1.4 Glu986 (MLLDALTLASSRGPLRKRSVAVAKAKPKFSISPDSLSYPYDVPDYAYPYDVPDYIKIATRKYLGKQNVYDIGVERDHNFALKNGFIASNCFN).
WT and Trp434TAG Shaker and all Nav1.4 cRNAs were transcribed from a pBSTA vector using the mMessage mMachine T7 Kit (Thermo Fisher Scientific, Waltham, MA). Purification of the cRNA from the transcription reaction was conducted on columns from the RNeasy Mini Kit (Qiagen, Hilden, Germany). Concentration was determined by absorbance measurements at 260 nm and quality was confirmed on a 1% agarose gel (RNase-free).
A modified (G73) version of Tetrahymena thermophila tRNA, THG73 (Saks et al., 1996), lacking the last two nucleotides was transcribed in vitro using CellScript T7-Scribe Standard RNA IVT Kit (CELLSCRIPT). DNA template 20 µg of annealed oligonucleotides coding for the tRNA preceded by a T7 promoter were used (Forward: ATTCGTAATACGACTCACTATAGGTTCTATAGTATAGCGGTTAGTACTGGGGACTCTAAATCCCTTGACCTGGGTTCGAATCCCAGTAGGACCGC; Reverse:
GCGGTCCTACTGGGATTCGAACCCAGGTCAAGGGATTTAGAGTCCCCAGTACTAACCGCTATACTATAGAACCTATAGTGAGTCGTATTACGAAT; Integrated DNA Technologies, Coralville, IA). The total reaction volume was adjusted to 100 µl and the kit reagents were added in the following amounts: 10 µl of 10X T7-Scribe transcription buffer, 7.5 µl of each nucleotide (100 mM stocks), 10 µl of 100 mM Dithiothreitol, 2.5 µl ScriptGuard RNase Inhibitor, 10 µl T7-Scribe enzyme solution. After the reaction was incubated for 4–5 hr at 37°C, the DNA template was digested with 5 µl DNase (1 U/µl) provided with the kit for 30–60 min. The tRNA was extracted from the reaction with acidic phenol chlorophorm (5:1, pH 4.5) and precipitated with ethanol. The precipitates tRNA was pelleted, washed, dried and resuspended in 100 µl DEPC-treated water and further purified with Chroma Spin-30 columns (Clontech, Mountain View, CA). The procedure yielded roughly 100 µl of ~5 µg/µl tRNA that was stored in aliquots at −80°C. Prior to the ligation reaction of the 73-mer THG73 to the ncAA-pCA conjugate (or just the pCA to yield a complete, nonacylated tRNA), the tRNA was folded in 10 mM HEPES (pH 7.4) by heating at 94°C for 3 min and subsequent gradual cool-down to ~10°C. 25 µg of folded tRNA in 30 µl of 10 mM HEPES (pH 7.4) were mixed with 28 µl of DEPC-treated water, 8 µl of 3 mM ncAA-pCA (or pCA) dimethylsulphoxide stock, 8 µl of 10X T4 RNA Ligase 1 buffer, 1 µl 10 mM ATP and 5 µl of T4 RNA Ligase 1 (New England Biolabs. Ipswich, MA) then incubated at 4°C for 2 hr. The (misacylated) tRNA was extracted from the samples with acidic phenol chlorophorm (5:1, pH 4.5) and precipitated with ethanol. The precipitated tRNA pellets were washed, dried in a Speedvac and stored at −80°C.
Xenopus laevis oocytes (stage V and VI) were purchased from Ecocyte (Austin, TX). Prior to injection each tRNA pellet was resuspended in 1.5 µl of ice-cold 3 mM sodium acetate and pelleted at 21,000 xg, 4°C for 25 min. For two electrode voltage clamp experiments cRNA amounts are indicated in each figure. After injection, oocytes were kept in OR-3 (50% Leibovitz’s medium, 250 mg/l gentamycin, 1 mM L-glutamine, 10 mM HEPES pH = 7.6) at 18°C for ~24–48 hr.
TEVC was performed as described previously (Pless et al., 2011). In brief, voltage-dependent potassium currents were recorded in in standard Ringer (in mM: 116 NaCl, 2 KCl, 1 MgCl2, 0.5 CaCl2, 5 HEPES, pH 7.4) using an OC-725C voltage clamp amplifier (Warner Instruments, Hamden, CT). Glass microelectrodes backfilled with 3 M KCl had resistances of 0.5–3 MΩ. Data were filtered at 1 kHz and digitized at 10 kHz using a Digidata 1322 A (Molecular Devices, Sunnyvale, CA) controlled by the pClamp 9.2 software. Shaker currents were elicited by +10 mV voltage steps from a holding potential of −80 mV to +40 mV and Nav1.4 currents were elicited by +5 mV voltage steps from a holding potential of −100 mV from −80 mV to +40 mV. Clampfit 9.2 software was used for current analysis. Numbers of oocytes, current analysis and statistical significance are indicated in the appropriate figure legends or the main text. All values are presented as mean ± SEM. To determine statistical significance Student’s t-test (two-tailed distribution; two-sample equal variance) was performed. The threshold for significance was p=0.01.
Nav1.4 cDNA plasmids were transfected into HEK293 cells 12–18 hr after plating at 70% confluency using standard calcium-phoshate reagents. 18–24 hr post transfection the HEK293 cells were washed twice with PBS and scraped in PBS supplemented with PMSF (0.1 mM final concentration)and Benzamidine (0.75 µM final concentration) and spun at 4000 g. The cell pellet was resuspended in solubilization buffer (50 mM Tris-HCl pH 7.4, 250 mM NaCl, 1% Triton, 0.5 µg/ml Pepstatin, 2.5 µg/ml Aprotinin, 2.5 µg/ml Leupeptin, 0.1 mM PMSF, 0.75 mM Benzamidine) on ice and Vortex Genied for 15mins at 4 C. Non-solubilized protein was removed by centrifugation at 20,000 g. Equal cell-lysate was loaded on a 3-15% gradient SDS-page in the presence of 1% 2 mercaptoethanol, separated at 55 V O/N and transferred to 0.45 µM LF PVDF (Bio-Rad, CA, USA). PVDF was immunoblotted using anti-Human influenza hemagglutinin (HA) antibody (1:2000; Biolegend, CA) in 2% NF milk and imaged on LI-COR Odyssey Imaging System (LI-COR, NE, USA).
Mechanism of frequency-dependent broadening of molluscan neurone soma spikesThe Journal of Physiology 291:531–544.https://doi.org/10.1113/jphysiol.1979.sp012829
Optimization of the additive CHARMM all-atom protein force field targeting improved sampling of the backbone φ, ψ and side-chain χ(1) and χ(2) dihedral anglesJournal of Chemical Theory and Computation 8:3257–3273.https://doi.org/10.1021/ct300400x
Molecular driving forces determining potassium channel slow inactivationNature Structural & Molecular Biology 14:1062–1069.https://doi.org/10.1038/nsmb1309
Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systemsThe Journal of Chemical Physics 98:10089–10092.https://doi.org/10.1063/1.464397
OpenMM 4: a reusable, extensible, hardware independent library for high performance molecular simulationJournal of Chemical Theory and Computation 9:461–469.https://doi.org/10.1021/ct300857j
C-type inactivation of voltage-gated K+ channels: pore constriction or dilation?The Journal of General Physiology 141:151–160.https://doi.org/10.1085/jgp.201210888
CHARMM-GUI PDB manipulator for advanced modeling and simulations of proteins containing nonstandard residuesAdvances in Protein Chemistry and Structural Biology 96:235–265.https://doi.org/10.1016/bs.apcsb.2014.06.002
RF1 knockout allows ribosomal incorporation of unnatural amino acids at multiple sitesNature Chemical Biology 7:779–786.https://doi.org/10.1038/nchembio.657
Comparison of simple potential functions for simulating liquid waterThe Journal of Chemical Physics 79:926–935.https://doi.org/10.1063/1.445869
Atomic basis for therapeutic activation of neuronal potassium channelsNature Communications 6:8116.https://doi.org/10.1038/ncomms9116
Modulation of C-type inactivation by K+ at the potassium channel selectivity filterBiophysical Journal 74:1840–1849.https://doi.org/10.1016/S0006-3495(98)77894-4
Update of the CHARMM all-atom additive force field for lipids: validation on six lipid typesThe Journal of Physical Chemistry B 114:7830–7843.https://doi.org/10.1021/jp101759q
All-atom empirical potential for molecular modeling and dynamics studies of proteinsThe Journal of Physical Chemistry. B 102:3586–3616.https://doi.org/10.1021/jp973084f
Extending the treatment of backbone energetics in protein force fields: limitations of gas-phase quantum mechanics in reproducing protein conformational distributions in molecular dynamics simulationsJournal of Computational Chemistry 25:1400–1415.https://doi.org/10.1002/jcc.20065
Probing the cavity of the slow inactivated conformation of shaker potassium channelsThe Journal of General Physiology 129:403–418.https://doi.org/10.1085/jgp.200709758
Three-dimensional architecture and gating mechanism of a K+ channel studied by EPR spectroscopyNature Structural Biology 5:459–469.https://doi.org/10.1038/nsb0698-459
Asymmetric functional contributions of acidic and aromatic side chains in sodium channel voltage-sensor domainsThe Journal of General Physiology 143:645–656.https://doi.org/10.1085/jgp.201311036
Numerical integration of the cartesian equations of motion of a system with constraints: molecular dynamics of n-alkanesJournal of Computational Physics 23:327–341.https://doi.org/10.1016/0021-9991(77)90098-5
An engineered Tetrahymena tRNAGln for in vivo incorporation of unnatural amino acids into proteins by nonsense suppressionThe Journal of Biological Chemistry 271:23169–23175.https://doi.org/10.1074/jbc.271.38.23169
CHARMM general force field: a force field for drug-like molecules compatible with the CHARMM all-atom additive biological force fieldsJournal of Computational Chemistry 31:671–690.https://doi.org/10.1002/jcc.21367
A facile system for genetic incorporation of two different noncanonical amino acids into one protein in Escherichia coliAngewandte Chemie International Edition 49:3211–3214.https://doi.org/10.1002/anie.201000465
Genetic incorporation of multiple unnatural Amino acids into proteins in mammalian cellsAngewandte Chemie International Edition 52:14080–14083.https://doi.org/10.1002/anie.201308137
How does the W434F mutation block current in Shaker potassium channels?The Journal of General Physiology 109:779–789.https://doi.org/10.1085/jgp.109.6.779
A conserved residue cluster that governs kinetics of ATP-dependent gating of Kir6.2 potassium channelsThe Journal of Biological Chemistry 290:15450–15461.https://doi.org/10.1074/jbc.M114.631960
Richard AldrichReviewing Editor; The University of Texas at Austin, United States
In the interests of transparency, eLife includes the editorial decision letter and accompanying author responses. A lightly edited version of the letter sent to the authors after peer review is shown, indicating the most substantive concerns; minor comments are not usually included.
Thank you for submitting your article "Atomic mutagenesis in potassium channels with engineered stoichiometry" for consideration by eLife. Your article has been reviewed by two peer reviewers, including Baron Chanda (Reviewer #2), and the evaluation has been overseen by Richard Aldrich as the Senior Editor and Reviewing Editor.
The reviewers have discussed the reviews with one another and the Reviewing Editor has drafted this decision to help you prepare a revised submission.
The paper of Lueck et al. is a technically sound and very creative report that attempts to address some fundamental questions about the atomic details underlying C-type inactivation gating in K+ channels. For many years it has been known that a single substitution at the position 443 of the Shaker K+ channel (W434F) yields a protein non-functional, presumably deeply C-type inactivated. C-type inactivation gating is widely accepted to be the functional collapse of the channel's selectivity filter. Structural changes at the channel's selectivity filter have been proposed to produce its structural collapse causing C-type inactivation and ceasing in this way ion conduction. The Lueck et al. paper deal with the identification of molecular interactions that are crucial in determining the conductive conformation of the Shaker K+ channel selectivity filter. Particularly interesting to keep the Shaker channel conductive is the highly conserved indole hydrogen bond between the pair Trp434-Asp447. It has been shown before by Dr. Ahern group that the incorporation of a non-hydrogen bonding homologue of tryptophan at position 434 obliterate ion conduction, presumably by constitutively C-type inactivating the Shaker channel, which is similar to the phenotype displayed by the mutation W343F. In this sequel of the previous paper by Pless, Galpin et al. 2013, the authors make two very provocative and significant observations that notably contribute to cement our understanding about C-type inactivation gating in K+ channels. First, they propose, based in MD simulation studies, that the Asp447 flips out toward the channel external vestibule and gives space for water molecules to enter a "crack or crevice" behind the channel selectivity filter. This event (the hydration of this crack) and/or a change in the electrostatic potential around the potassium binding sites of the channel's selectivity filter, trigger the functional collapse of latter. Second, the authors in a very innovative way present a new strategy to dissect the contribution of each subunit within the homotetramer channel to the C-type inactivation process. These combined strategies and the ideas derived from this experimental approach I think will be well received among the members of the ion channel community.
The following revisions will be required for acceptance:
1) It would have been very powerful and conclusive to perform an SDS PAGE analysis documenting the tetrameric "dimer of dimer complexes" obtained after removal of the ERret motif and putative recovering of trafficking to the oocyte membranes.
2) A more complete statistical analysis is required. Figure 2 is missing a more detailed characterization of the ERret-intein cRNAs constructs to conclude they show "WT-like currents". Please provide the n, statistics and standard deviation of the current expressed by this construct as well as the GV curves to make sure the gating machinery after all these manipulations remain intact. The same applied to Figure 3.https://doi.org/10.7554/eLife.18976.013
[…] The following revisions will be required for acceptance:
1) It would have been very powerful and conclusive to perform an SDS PAGE analysis documenting the tetrameric "dimer of dimer complexes" obtained after removal of the ERret motif and putative recovering of trafficking to the oocyte membranes.
We agree completely with this sentiment. The resubmitted manuscript now contains new Western blot data demonstrating the high fidelity of in-cell protein ligation using newly designed sodium hemi-channels (Figure 4B). The reason for using sodium channels was two-fold. First, to further demonstrate the versatility of the approach in terms of channel type, second, the expression system used for these experiments was HEK cells as opposed to Xenopus oocytes. Thus the data support the notion that the ERret/intein sequences are generally tolerated between channel and cell types.
2) A more complete statistical analysis is required. Figure 2 is missing a more detailed characterization of the ERret-intein cRNAs constructs to conclude they show "WT-like currents". Please provide the n, statistics and standard deviation of the current expressed by this construct as well as the G/V curves to make sure the gating machinery after all these manipulations remain intact. The same applied to Figure 3.
We now include a table (Table 1) summarizing the potassium and sodium channel electrophysiological data (GV curves and Boltzmann fits) with pairwise statistical analysis and sample sizes.https://doi.org/10.7554/eLife.18976.014
- Jason D Galpin
- Christopher A Ahern
- Daniel T Infield
- Jason D Galpin
- Christopher A Ahern
- Jing Li
- Benoît Roux
- Christopher A Ahern
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
CAA is supported by NIH/NIGMS (GM106569) and the American Heart Association Established Investigator (A22180002). BR is supported by NIH/NIGMS (GM GM062342). This collaborative work was supported by the Membrane Protein Structural Dynamics Consortium, which is funded by NIH/NIGMS (GM087519). We thank Dr. Stephan Pless comments on the manuscript. The authors declare no conflict of interest.
- Richard Aldrich, The University of Texas at Austin, United States
© 2016, Lueck et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.