Differential ion dehydration energetics explains selectivity in the non-canonical lysosomal K+ channel TMEM175
Abstract
Structures of the human lysosomal K+ channel transmembrane protein 175 (TMEM175) in open and closed states revealed a novel architecture lacking the canonical K+ selectivity filter motif present in previously known K+ channel structures. A hydrophobic constriction composed of four isoleucine residues was resolved in the pore and proposed to serve as the gate in the closed state, and to confer ion selectivity in the open state. Here, we achieve higher-resolution structures of the open and closed states and employ molecular dynamics simulations to analyze the conducting properties of the putative open state, demonstrating that it is permeable to K+ and, to a lesser degree, also Na+. Both cations must dehydrate significantly to penetrate the narrow hydrophobic constriction, but ion flow is assisted by a favorable electrostatic field generated by the protein that spans the length of the pore. The balance of these opposing energetic factors explains why permeation is feasible, and why TMEM175 is selective for K+ over Na+, despite the absence of the canonical selectivity filter. Accordingly, mutagenesis experiments reveal an exquisite sensitivity of the channel to perturbations that mitigate the constriction. Together, these data reveal a novel mechanism for selective permeation of ions by TMEM175 that is unlike that of other K+ channels.
Editor's evaluation
This manuscript explores the mechanisms of permeation and selectivity in the unusual potassium-selective ion channel TMEM175, which lacks a canonical selectivity filter. The study is led by molecular dynamics simulations and free energy calculations, complemented by a cryoEM analysis and electrophysiological recordings. The authors propose a novel, single ion-based mechanism of permeation, together with a partial dehydration-driven selectivity mechanism. This study will appeal to readers interested in the structure and function of ion channels and in molecular mechanisms of ion translocation.
https://doi.org/10.7554/eLife.75122.sa0Introduction
Transmembrane protein 175 (TMEM175) is a potassium (K+)-selective cation channel that is evolutionarily distinct from all other known ion channels. In mammals, TMEM175 is expressed in lysosomes, where it is critical for maintaining lysosomal homeostasis (Cang et al., 2015). As a K+-selective channel, TMEM175 contributes to establishing the membrane potential across the lysosomal membrane (Cang et al., 2015). Loss of TMEM175 function can lead to numerous defects in the endolysosomal system including dysregulation of lysosomal pH, defects in autophagy and mitophagy, reduced glucocerebrosidase activity, and an increased susceptibility to α-synuclein cytotoxicity (Jinn et al., 2019; Jinn et al., 2017; Krohn et al., 2020). Dysregulation of lysosomal function is commonly associated with neurodegenerative disease and point mutations in TMEM175 have been identified that are strongly associated with development of Parkinson’s disease (PD) (Blauwendraat et al., 2019; Iwaki et al., 2019; Jinn et al., 2019; Jinn et al., 2017; Krohn et al., 2020; Nalls et al., 2014; Wie et al., 2021). For example, the M393T point mutation is a loss-of-function mutation associated with the increased likelihood and early onset of PD and the Q65P mutant is a gain-of-function mutation associated with a reduced likelihood of developing of PD (Wie et al., 2021).
Consistent with being evolutionarily distinct from other ion channels, TMEM175 channels from prokaryotes and humans adopt a unique architecture relative to all known ion channel structures (Brunner et al., 2020; Lee et al., 2017; Oh et al., 2020). Structures of two prokaryotic TMEM175 homologs in closed states revealed a homotetrameric channel possessing a central ion conduction pathway lined by several layers of hydrophobic residues from the first transmembrane (TM) helix of each of the four protomers (Brunner et al., 2020; Lee et al., 2017). Despite similar overall structures, analysis of the two prokaryotic TMEM175 structures led to conflicting proposals for how TMEM175 channels selectively permeate K+ ions (Brunner et al., 2020; Lee et al., 2017). However, human TMEM175 (hTMEM175) channels are homodimers, rather than homotetramers, and are ~30-fold selective for K+ over Na+, whereas prokaryotic TMEM175 channels exhibit only 2- to 4-fold selectivity for K+/Na+ (Cang et al., 2015). Because of the notable structural and functional differences between human and prokaryotic TMEM175, it remains unclear if the mechanisms gleaned from analyses of prokaryotic TMEM175 channels are generalizable for hTMEM175.
To begin to understand the mechanisms that underlie TMEM175 function in lysosomes, we determined the structure of hTMEM175 by cryo-EM and identified two conformational states (Oh et al., 2020). hTMEM175 protomers consist of two homologous repeats of six TM helices, which assemble into a dimer to form the channel. A 2-fold symmetric ion conduction pathway runs along the center of the dimer, lined by two copies of TM1 and TM7, the first helices of each repeat. Unlike most cation channels, including canonical K+ channels such as KcsA and Shaker, the pore of TMEM175 contains no clusters of hydrophilic side chains or backbone atoms that control selectivity (Doyle et al., 1998; Long et al., 2005). Rather, the most striking feature in the pore of TMEM175 is a narrow hydrophobic constriction formed by the side chains of the conserved Ile46 from TM1 and the conserved Ile271 from TM7, which we termed the isoleucine constriction (Oh et al., 2020). In one of the two cryo-EM structures, the isoleucine constriction is too narrow to accommodate K+ ions and thus we assigned it as a closed state. In the other structure, however, the isoleucine constriction is sufficiently dilated to accommodate partially hydrated K+ ions, and therefore we proposed that this structure corresponds to a conductive state of the channel. Moreover, because ions would need to be partially dehydrated to traverse the pore, it was proposed that the relative energetics of K+ and Na+ dehydration through the isoleucine constriction contributes to the ion selectivity of TMEM175 (Lee et al., 2017; Oh et al., 2020). While plausible, this proposal relies on the assumption that the putative open state is indeed conductive; however, assigning functional states to cryo-EM structural snapshots is challenging. Here, we employ high-resolution cryo-EM structure determination, molecular dynamics (MD) simulations, and electrophysiological analyses to investigate the mechanisms of ion permeation and selectivity in TMEM175.
Results
Improved cryo-EM structures of hTMEM175
Taking advantage of recent advances in cryo-EM image analysis, we reprocessed our published cryo-EM images of TMEM175 in KCl to achieve higher resolution (Punjani and Fleet, 2021; Punjani et al., 2020; Zivanov et al., 2018). The resolution of the putative open state improved from 2.64 to 2.45 Å and the closed state improved from 3.03 to 2.61 Å (Figure 1—figure supplement 1 and Appendix 1—table 1). While the improved structures are similar to those previously determined with an all-atom RMSD of 0.3 Å for the open state and 0.4 Å for the closed state, the higher-resolution reconstructions allowed us to better define the ion conduction pathways (Figure 1A–D). In both states, the ion conduction pathway is 2-fold symmetric, narrow, and largely devoid of polar contacts, particularly near the center due to the presence of several hydrophobic side chains including Ile46, Val50, Leu53, Ile271, Leu275, and Leu278 (Figure 1C–D). In the closed state, the isoleucine constriction formed by Ile46 and Ile271 is too narrow to allow ion permeation (Figure 1D). The side chains of Ile271 point toward the channel axis, narrowing the pore to a minimum radius of 0.2 Å (Figure 1D and E). In the open state, by contrast, conformational changes in the pore-lining helices TM1 and TM7 helices dilate the isoleucine constriction to a minimum radius of 1.6 Å, which appears to be sufficient to accommodate a partially dehydrated K+ ion (Figure 1C and E). In its dilated conformation, the isoleucine constriction adopts a pseudo 4-fold symmetric configuration with the side chains of both Ile46 and Ile271 adopting the mt rotamer (Figure 1—figure supplement 2).
The higher-resolution reconstructions also allowed an improved modeling of five ion-binding sites, which we call K1 through K5, and ordered water molecules in the pore of the open state. Two non-protein densities are now resolved at the K3 and K4 sites that flank the isoleucine constriction (Figure 1—figure supplement 2). We assigned these densities as K+ ions based on our previous comparison with a structure determined in the presence of Cs+ (Oh et al., 2020). The K3 and K4 ion-binding sites are positioned 3.7–4.0 Å away from the side chains of Ile46 and Ile271 and are coordinated exclusively by water molecules (Figure 1—figure supplement 2). Notably, the ions in the K3 and K4 sites are only partially hydrated. The K3-binding site on the cytoplasmic site of the constriction is coordinated by four ordered water molecules that are between 2.8 and 3.0 Å away, while the K4 site on the luminal side is coordinated by four water molecules that are between 2.9 and 3.5 Å away. The distance between the two sites is 2.8 Å and thus they are unlikely to be simultaneously occupied by ions.
For the closed state, we can now model three ion-binding sites and 25 water molecules into non-protein densities (Figure 1D). In the cytosolic region of the pore, where the channel undergoes minimal conformational changes during gating, the K1 and K2 ion-binding sites and water molecules occupy almost identical positions to those resolved in the open state (Figure 1—figure supplement 2). In contrast, there is little correspondence between the ion and water configurations at the isoleucine constriction or on the luminal side of the pore due to the conformational changes that occur to the channel during channel gating. The inward movement of Ile271 displaces the ion-binding sites on either side of the isoleucine constriction. Movements of Thr49 and Thr274, which both coordinate water molecules in the open state, result in a movement of the luminal K5 ion-binding site toward the luminal entrance of the pore. Thus, the higher-resolution structures reveal how changes in the protein structure alter the ion conduction pathway to prevent ion permeation in the closed state.
Energetics and mechanism of K+ permeation
To ascertain whether the structure of the putative open state is indeed permeable to K+, we turned to all-atom MD simulations of the channel in a model phospholipid bilayer (Figure 2—figure supplement 1). Specifically, using the enhanced-sampling methodology known as multiple-walker Metadynamics (Raiteri et al., 2006), we induced the permeation of K+ across the channel at 0 mV and symmetric 100 mM KCl and evaluated the associated potential of mean force (or free-energy profile). Importantly, these calculations used a newly formulated reaction coordinate (Materials and methods) designed to avoid a priori assumptions about the mechanism of permeation, such as the number of ions that reside within the pore at a given time.
The results from this analysis are summarized in Figures 2 and 3. We observe that K+ ions reach the isoleucine constriction readily from either side of the membrane, through a series of shallow free-energy barriers and transient-binding sites (Figure 2A). At the isoleucine constriction, however, the free energy increases steeply, peaking at about 7 kcal/mol (Figure 2A). These features recapitulate the pattern of non-protein densities resolved in the pore in the reconstruction of the putative open state; densities near Ser38 and Ala42 on the cytoplasmic side and near Val50 on the luminal side appear as metastable states in the calculated free-energy profile. Densities right below and above the central isoleucine constriction also appear as shoulders in the free-energy profile. Unlike in canonical K+ channels, however, the simulation data reveals no evidence of a multi-ion process involving concurrent occupancy of proximal-binding sites near the constriction (Figure 2—figure supplement 2). Instead, K+ ions approach and cross this constriction individually. It is apparent, too, that permeation requires a gradual but drastic depletion of the ion hydration shells. As K+ traverses the isoleucine constriction, the first two hydration shells are reduced to only ~9 molecules, down from 31 in the bulk; the first shell averages to about 4, down from 7 (Figures 2B, C ,, 3A).
Given the largely hydrophobic nature of the TMEM175 pore, it is logical that the free-energy barrier for permeation is significantly higher than those in canonical K+ channels. The finding that this barrier is only 7 kcal/mol is however intriguing, given that the cost of K+ dehydration is much larger than that (up to about 80 kcal/mol for full dehydration). This observation indicates that other factors must facilitate ion permeation in TMEM175. Indeed, further examination of the simulated trajectories using Poisson theory (Materials and methods) reveals that acidic side chains at both the luminal and cytosolic entrances of the pore generate a strongly favorable electrostatic field that spans the length of the pore (Figure 2D). Some of these residues also transiently coordinate the permeating ions as they enter and exit the pore, along with a few other carbonyl and carboxyl groups in side chains and backbone (Figure 2C, Figure 2—figure supplement 3). However, for a span of about 15 Å centered at the isoleucine constriction, direct protein contacts are minimal. That the barrier for ion conduction in TMEM175 is surmountable is thus a consequence of the degree of hydration retained by the permeating ions as well as the electrostatic field generated by the protein that attracts them to the pore interior.
Integration of the potential-of-mean-force profile obtained at 0 mV and 100 mM KCl translates into a single-channel conductance of approximately 0.23 pS (Materials and methods). At 500 mV, and assuming the structure of the channel is completely unresponsive to voltage, this conductance would in turn produce a 0.1 pA current, or a rate of K+ permeation of about 0.6 ions/μs. To assess these estimates, we calculated an independent 1-μs-long MD simulation using conventional sampling but under 500 mV (lumen positive) and 400 mM KCl. This trajectory revealed a stable pore conformation and two K+ permeation events, confirming that the open state is functionally conductive and validating the single-ion mechanism inferred from the free-energy simulations (Figure 3C and D; Figure 3—video 1 and Figure 3—video 2). K+ ions are seen to dwell at the shallow-binding site adjacent to the constriction, on the luminal side, before returning to the bulk or rapidly traversing the free-energy barrier, along with 4 water molecules in the first hydration shell (Figure 3C).
Differentials in ion dehydration relative to bulk explain selectivity for K+ over Na+
Like canonical K+ channels, TMEM175 channels are more permeable to K+ than Na+. However, in the absence of the selectivity filter seen in canonical K+ channels, it has been unclear how to explain their selectivity in structural terms. To begin to bridge this gap, we again turned to MD simulations. Specifically, we repeated the enhanced-sampling simulations carried out to examine K+ permeation after substituting the KCl solution with NaCl and examined the resulting data analogously (Figures 2 and 3). The potential-of-mean-force profile for Na+ is similar to that obtained for K+ in its overall features, including the rate-limiting barrier at the isoleucine constriction (Figure 2A) . The largest difference between the profiles is the magnitude of the barrier at the isoleucine constriction, which peaks at about 8 kcal/mol relative to bulk for Na+ compared to about 7 kcal/mol for K+. Integration of the potential-of-mean-force profile obtained for Na+ translates into a single-channel conductance of approximately 0.037 pS, that is, the calculated conductance for Na+ is 6-fold smaller than for K+. While available experimental estimates for human indicate TMEM175 a stronger preference for K+, the qualitative agreement between experimental and calculated values suggest the simulations capture the essence of the mechanism of K+ selectivity. As observed for K+, Na+ permeation proceeds via a single-ion mechanism (Figure 2—figure supplement 2, Figure 3B). Permeation is again opposed by the drastic dehydration of the ion (Figure 2B), and favored by its interaction with the electrostatic field created by the protein within the pore (Figure 2D), as well as by transient interactions with multiple residues on both cytoplasmic and luminal sides (Figure 2C). These stabilizing factors are not identical but appear comparable to those observed for K+. By contrast, the degree of dehydration required for permeation clearly differs. We observe that the depletion in the both the first and second solvation shells, relative to bulk numbers, is significantly smaller for Na+ than for K+. As will be further discussed below, this relative difference in dehydration energetics likely explains why TMEM175 is only ~30-fold more permeable to K+, even though the bulk water selectivity for Na+ is over a trillion-fold. In other words, like for K+ channels that feature a canonical selectivity filter, it can be said that the narrow pore of TMEM175 favors Na+ over K+, but not as much as bulk water. The channel is thus more permeable to K+ than Na+.
Role of the isoleucine constriction in ion selectivity and permeation
The isoleucine constriction constitutes the largest free-energy barrier to the permeation of K+ and Na+ in the pore of TMEM175 (Figure 2). This narrow constriction forces the ions to shed much of their hydration shells in order to traverse the pore. It was previously demonstrated that mutating both Ile46 and Ile271 to asparagine diminishes the selectivity of TMEM175, indicating that the hydrophobicity of the isoleucine constriction is crucial for selectivity (Lee et al., 2017). To determine if the dimensions of the constriction are also important, we used MD simulations to calculate the change in free energy that results when a K+ ion at the isoleucine constriction is exchanged with a Na+ ion, both for wild-type TMEM175 and for a I46V mutant with an expanded isoleucine constriction (Materials and methods). In both cases the exchange results in a free-energy gain (Figure 4A), consistent with the notion that Na+ is the favored species at the constriction. However, this gain is smaller than that observed for the exchange of a K+ ion in bulk water with a Na+ ion. As discussed above, it is relative favorability for K+ at the constriction compared to bulk solution that would be the expected result for a K+-selective channel (Figure 4A). Specifically, the differential relative to bulk is 1.2 kcal/mol in favor of K+, similar to the ~1 kcal/mol difference in the peak values of the free-energy curves discussed in the previous section (Figure 2A). Na+ is however favored to a larger extent in the I46V mutant, whose isoleucine constriction has a similar hydrophobicity but is expanded by the loss of a methyl group from the side chain of Ile46 in each protomer (Figure 4A); specifically, the calculations indicate that this mutant is about 2.5-fold (or 0.5 kcal/mol) less K+-selective than the wild-type channel (Figure 4A). These results again correlate with the degree of dehydration of each ion, relative to bulk values (Figure 4B).
Using whole-cell electrical recordings, we experimentally assessed the effects of expanding the size of the isoleucine constriction on selectivity by measuring the selectivity of Cs+ over Na+ for comparison with the computational analyses. We choose Cs+ rather than K+ because Cs+ blocks most canonical K+ channels, thereby reducing endogenous currents. Nevertheless, whole-cell currents recorded from cells expressing the I46A, I46V, I271A, and I271V mutants in a bi-ionic Cs+/Na+ condition were much smaller than those recorded from wild-type TMEM175 and in some cases similar to currents recorded from non-transfected cells, which hampered accurate determination of ion permeability ratios via measuring reversal potentials (Erev) (Figure 5 and Figure 5—figure supplement 3). Thus, isolation of the TMEM175-specific currents from the endogenous HEK293T and background leak was necessary to evaluate the selectivity of the TMEM175 mutants. To this end, we first performed a background subtraction with the TMEM175 inhibitor 4-aminopyridine (4-AP) (Cang et al., 2015). However, while 4-AP can effectively inhibit wild-type TMEM175, it was ineffective against the mutants and thus not suitable for isolating the TMEM175-specific currents (Figure 5—figure supplements 3 and 4). We therefore repeated the background subtraction experiments with a novel TMEM175 inhibitor that we have developed, which we call AP-6 (Figure 5 and Figure 5—figure supplements 1–3). Following AP-6 background subtraction, I46V and I271V displayed Cs+-selective currents with Erev of –13 and –16 mV, corresponding to Cs+ over Na+ permeability ratios (PCs/Na) of 1.7 and 1.9, respectively (Figure 5). Using the same background subtraction approach, we estimated the PCs/Na of wild-type TMEM175 to be 20 (Erev of –76 mV). This 10-fold reduction in PCs/Na ion permeation ratio is in line with the computational analyses of the free-energy difference in resulting from K+ to Na+ ion exchange at the isoleucine constriction of I46V and wild type. We were not able to detect any AP-6-sensitive currents from the I46A or I271A mutants either due their insensitivity toward AP-6 or inactivity of the mutant channels. These data demonstrate that increasing the size of the isoleucine constriction leads to a clear reduction in K+ selectivity. Thus, the isoleucine constriction, which is universally conserved among eukaryotic TMEM175 channels, serves as a hydrophobic selectivity filter that is exquisitely sensitive to very subtle changes in its structure.
Discussion
Permeation through selective ion channels arises from the interplay between channel-ion interactions, solvent-ion interactions, and in some cases, ion-ion interactions. The energetic cost associated with fully or partially dehydrating an ion such as K+ is so large that permeation would only very rarely occur without favorable interactions with the channel as the ion traverses the conduction pathway. In canonical K+ channels, backbone and side chain oxygen atoms provide a series of direct coordinating interactions that closely resemble those experienced by the ion in bulk water. In addition, these channels are often capable of attracting multiple ions to the pore concurrently. Thus, the cost of dehydration is almost entirely offset by channel-ion and ion-ion interactions, and the resulting free-energy barriers for permeation are small. In contrast, such interactions are lacking for much of the TMEM175 pore, and entirely absent near the isoleucine constriction. Lacking these interactions, it is unlikely that direct protein-ion interactions impart ion selectivity on hTMEM175 channels as was proposed based upon the structural analysis of a TMEM175 homolog from the prokaryote Marivirga tractuosa (Brunner et al., 2020). K+ ions crossing the pore must nevertheless shed much of their hydration shells; only ~9 of the 31 water molecules in its first two shells remain at the constriction. In addition, ions permeate TMEM175 one at a time, without direct or indirect knock-on interactions. Permeation is however strongly favored by a delocalized electrostatic field across the pore, which offsets much of the cost of ion dehydration, resulting in a free-energy barrier of ~7 kcal/mol (Figure 2). While this is a high barrier, it is not insurmountable. Indeed, considering that MD simulations using conventional force fields tend to underestimate ion conduction rates, in some cases up to 10-fold (Allen et al., 2006; Mironenko et al., 2021), we estimate the unitary conductance of hTMEM175 is 0.1–0.5 pS, which is well within the broad range of experimentally measured permeation rates for other selective cation channels such as ~100 pS for Slo K+ and 9–24 fS for CRAC channels (Tao et al., 2016; Zweifach and Lewis, 1993). Thus, despite the hydrophobicity of the isoleucine constriction, the unique architecture and amino acid make-up of the TMEM175 pore facilitates ion permeation.
Importantly, the balance between protein-ion interactions and solvent-ion interactions is dependent on the ion type, and so these opposing forces also govern ion selectivity in channels. For example, when rationalizing differences in permeability between K+ and Na+, it is useful to keep in mind that dehydration of Na+ is much more costly than that of K+; the penalty for full dehydration is ~18 kcal/mol greater, rendering it a trillion-fold less probable. However, the observed selectivity of biological K+ channels against Na+ is only 1000-fold or less, implying Na+ establishes much more favorable interactions with the selectivity filters of this kind of K+ channels, as noted in previous studies (Kim and Allen, 2011; Kopec et al., 2018). In other words, what ultimately determines the specificity of a channel is the sum of the relative cost of dehydration of the competing ions and the differential in their free energy of interaction with the channel. Remarkably, for TMEM175 these favorable protein-ion interactions are imparted by residues near the luminal and cytoplasmic entrances to the pore, quite far from the isoleucine constriction. Thus, we propose that the ion selectivity of TMEM175 primarily results from a differential in the degree of dehydration required for permeation, relative to bulk water. Incidentally, whereas recent studies have revealed differing results for simulations of ion channels performed with different force fields (Klesse et al., 2020; Ocello et al., 2020), a mechanism that relies on dehydration energetics, rather than specific ion-protein or ion-ion interactions, will be reasonably described by the force field employed in our simulations to reproduce the relative free energies of hydration (and thus dehydration) of the alkali cations (Figure 4A).
Together, our analyses provide a structural and energetic model for the unique mechanisms of ion permeation and selectivity in TMEM175 channels. By combining a narrow hydrophobic sieve with a favorable electrostatic field, hTMEM175 selectively permeates K+ ions across the lysosomal membrane, as was proposed by Lee and colleagues for the TMEM175 homolog from Chamaesiphon minutus (Lee et al., 2017). As the large conductance Slo1 K+ channel, whose K+ permeation rate is several orders of magnitude higher than TMEM175, has also been observed in lysosomes, it will be important to uncover the specific roles for TMEM175’s distinct conduction and selectivity properties in the lysosome and why mutations of TMEM175 can both increase and decrease the likelihood of developing PD (Wang et al., 2017; Zhong et al., 2016).
Materials and methods
Analysis of electron microscopic images
Request a detailed protocol7907 images of TMEM175 in KCl previously acquired in two data sets as 40-frame super-resolution movies (0.544 Å/pixel) using a Gatan K2 (Oh et al., 2020) were gain corrected, Fourier cropped by two, and aligned using whole-frame and local motion correction algorithms by Motioncor2 (Zheng et al., 2017) (1.088 Å/pixel). Whole-frame CTF parameters were determined using CTFfind 4.1.14 (Rohou and Grigorieff, 2015). Particles were automatically selected in Relion 3.0 using templates previously generated from 2D classification, resulting in 4,153,614 particles (Zivanov et al., 2018). False-positive selections and contaminants were excluded from the data using multiple rounds of heterogeneous classification in cryoSPARC v3.2 using the open and closed states, as well as several decoy classes generated from noise particles via ab initio reconstruction in cryoSPARC v3.2, resulting in a stack of 636,148 particles (Punjani et al., 2017). After Bayesian polishing in Relion and local CTF estimation and higher-order aberration correction in cryoSPARC v3.2, a consensus reconstruction was determined at resolution of 2.8 Å (Zivanov et al., 2019). A second round of Bayesian polishing in Relion 3 using a pixel size of 0.85 Å and box size of 384 yielded an improved consensus reconstruction at 2.5 Å. Iterative supervised heterogenous classification using open and closed maps low pass filtered to 6 Å as references resulted in 261,536 particles in the open state and 163,651 particles in the closed state. A 2.45 Å reconstruction of the open state particles and a 2.61 Å reconstruction of the closed state particles were obtained using non-uniform refinement in cryoSPARC v3.2 employing global and per-particle CTF correction. The reconstructions were subjected to density modification using the two unfiltered half-maps with a soft mask in Phenix (Terwilliger et al., 2019).
Model building and coordinate refinement
Request a detailed protocolThe structures of open (PDB:6WC9) and closed (PDB:6WCA) hTMEM175 were docked into the density maps in COOT and manually adjusted to fit the density (Emsley et al., 2010). Densities corresponding to TM5 and TM6 (residues 174–251) were too poorly ordered and omitted from the model. The final models are composed of residues 30–173 and 254–476. Atomic coordinates were refined against the density modified map using phenix.real_space_refinement with geometric and Ramachandran restraints maintained throughout (Adams et al., 2010). Pore radius calculations were performed using HOLE (Smart et al., 1996).
Simulation systems and general specifications
Request a detailed protocolAll simulations were calculated with NAMD 2.12 using the CHARMM36 force field for protein and lipids (Best et al., 2012; Klauda et al., 2010; Pastor and Mackerell, 2011; Phillips et al., 2005). The simulations were carried out at constant temperature (298 K) and semi-isotropic pressure (1 atm), using periodic boundary conditions and an integration time step of 2 fs. Long-range electrostatic interactions were calculated using PME, with a real-space cut-off of 12 Å. Van der Waals interactions were computed with a Lennard-Jones potential, cut-off at 12 Å with a smooth switching function taking effect at 10 Å.
The simulations were based on the high-resolution cryo-EM structure of the hTMEM175 channel (PDB 6WC9). The specific construct studied includes residues 30–165 and residues 254–476. K+ ions and water molecules originally included in the cryo-EM structure were removed. All ionizable side chains were set in their default protonation state at pH 7, except for H57, which was protonated on account of its proximity to D279, E282, and D283. The protein construct was embedded in a pre-equilibrated hydrated palmitoyl-oleoyl-phosphatidyl-choline (POPC) lipid bilayer using GRIFFIN, and enclosed in a periodic orthorhombic box of ~100 × 100 × 111 Å3 in size (Staritzbichler et al., 2011). The two resulting simulation systems contain 222 POPC lipids, 24,342 water molecules, 49 Cl- ions, and either 43 K+ or 43 Na+ ions; that is, a salt concentration of 100 mM plus counterions to neutralize the protein net charge. The simulation systems were equilibrated following a staged protocol comprising a series of restrained simulations. The protocol consists of both positional and conformational restraints, gradually weakened over 100 ns. A third simulation system was prepared with 400 mM KCl, by adding 132 K+ and 132 Cl- ions to the 100 mM KCl system (replacing water), followed by a 100 ns equilibration.
Simulation of K+ permeation under voltage
Request a detailed protocolTo simulate the flow of K+ across TMEM175 from the luminal to the cytoplasmic side, a conventional MD trajectory of 1 μs was calculated for the 400 mM KCl condition under a constant electric field perpendicular to the membrane plane, of magnitude Ez = –0.1045 kcal/(mol Å e). Given the average length of the box in this direction (Lz = 110.5 Å), the corresponding TM potential is –500 mV (Φ=Ez × Lz × 0.0434 V e-1 (kcal/mol)-1) (Roux, 2008).
Enhanced-sampling simulation of K+ and Na+ permeation
Request a detailed protocolTo induce reversible ion permeation across TMEM175 at 0 mV and 100 mM KCl or 100 mM NaCl, we used the multiple-walker Metadynamics method (Raiteri et al., 2006). In Metadynamics, a biasing potential is introduced in an MD simulation to facilitate the exploration of configurational space, as defined by one or more reaction coordinates, also known as collective variables. This biasing potential consists of a series of Gaussian functions that expands as the simulation progresses, to foster the trajectory to visit high free-energy configurations. In the multiple-walker approach, several MD simulations are carried out in parallel, each sampling a different trajectory, but with a shared biasing potential that is also collectively constructed. To study K+ permeation, we used eight walkers, each producing a 600 ns trajectory. For Na+, we also used eight walkers, 700 ns each.
In biased-sampling simulations of channel permeation, it is common for the position of one or more ions along the pore axis to be used as collective variables. This choice might be reasonable when the number of ions involved in the mechanism is known. It was not known in our case, and hence, the type of collective variable used in these simulations was newly formulated to circumvent a priori assumptions. This variable was implemented in a modified version of PLUMED 1.3 (Bonomi et al., 2009). Specifically, the variable is a quantitative descriptor of the proximity between the ions in the simulation system and a virtual center within the pore. More precisely, the variable is:
where the k index identifies each of the K+ or Na+ ions in the system and is the Cartesian Z-component of the distance between ion k and a center-of-mass defined by a group of protein atoms. C is a positive constant (set to 2 Å) that avoids numerical instabilities when and is a smoothing parameter (set to 100 Å). To successfully induce ion permeation events, we used two variables of the kind specified by Equation 1, namely and , each defined in reference to a center of atoms. Center A is defined by the Cα atoms of residues 44–47 and 268–271, while center B is defined by the Cα atoms of residues 45–48 and 269–272. Centers A and B are therefore along the pore axis, separated by 1 Å and flanking the isoleucine constriction. The Gaussian functions used to gradually construct the biasing potentials applied to and had a width 0.25 Å and were added in 4 ps intervals. To reduce systematic errors near the boundary , reflected and inverted Gaussians were added beyond this boundary for both variables as previously described (Crespo et al., 2010). In the simulations of K+ permeation, the Gaussians height was gradually raised from 0.0035 to 0.007 kcal/mol in the first 30 ns of simulation, and gradually diminished back to 0.0035 kcal/mol after 100 ns. For Na+, the Gaussians height was gradually raised from 0.0035 to 0.007 kcal/mol in the first 50 ns of simulation and diminished back to 0.0035 kcal/mol after 250 ns.
Derivation of free energies and other quantitative descriptors from biased-sampling trajectories
Request a detailed protocolA post hoc reweighting procedure was used to derive unbiased averages and histograms from the MD trajectories enhanced by the Metadynamics biasing potential (Marinelli et al., 2009). In this approach, a time average of the Metadynamics biasing potential is calculated after it becomes approximately stationary across the range of ( , ) values of interest. (Here, the last 400 ns of simulation for K+ and the last 600 ns of simulation for Na+.) This effective potential, , is then used to adjust the statistical weight of each of the simulation snapshots considered in the analysis. Specifically, the statistical weight of a given snapshot is:
where is the Boltzmann constant and is the temperature. That is, snapshots that fall in easily accessible regions of () where the accumulated Metadynamics bias is large are given a greater statistical weight, while those in unfavorable regions where the accumulated bias is less are also assigned a smaller weight. The effective potential for each simulation snapshot was derived using a 2D discretization across the space of (), wherein each snapshot is assigned to a grid point and an associated value of using a recently developed tool for free-energy analysis (Marinelli and Faraldo-Gómez, 2021).
Following this approach, it is straightforward to derive unbiased estimates of the average ion occupancy along the length of the channel pore, and thus of the corresponding free-energy landscape. Specifically, let us define this length by an axis connecting the abovementioned centers A and B (which fluctuates with the channel but is approximately perpendicular to the membrane plane). For each ion k in the simulation system, let us also define and as the projection of the ion coordinates along the pore axis and the distance to that axis, respectively, for a given snapshot . If we divide the pore axis into a series of intervals or bins, the ion occupancy of a given bin α centered in zα is:
where the i index again denotes each trajectory snapshot included in the analysis, and the sum over k is restricted to ions whose coordinate falls into bin and that are also proximal to the pore axis ( Å). It follows that the potential of mean force (or free-energy profile) associated with this occupancy distribution along zα will be:
where is a constant arbitrarily selected so that in the bulk, where the profile is nearly flat.
We followed same reweighting approach to derive unbiased estimates for other descriptors, namely average ion coordination numbers and average ion-protein electrostatic interaction energies (definitions provided below). Specifically, the mean value of either of these observables at a given position along the pore axis is:
where i again denotes a trajectory snapshot, is the descriptor of interest of ion k, in that particular snapshot, and the sum over k is again restricted to ions within the pore and within bin α.
Calculation of ion coordination numbers
Request a detailed protocolTo quantify the coordination of K+ or Na+ by either water or protein oxygen atoms as a function of the ion position along the pore axis, we evaluated the following function of the ion-oxygen distances for each simulation snapshot and each ion k in the system:
where the j index denotes all possible coordinating oxygens, and denotes their distance to ion k in snapshot . Note this function is virtually identical to a step function cut-off at distance equal to . For computational expediency, this evaluation was carried out with PLUMED 1.3 (Bonomi et al., 2009). For the evaluation of the number of water molecules in the first hydration shell, was set to 3.5 Å for K+ and to 3.2 Å for Na+, which in each case corresponds to the position of the first minimum of the ion-oxygen radial distribution functions (RDF) in bulk water (data not shown). The same values of were used for calculating the coordination with protein oxygens. Following the same criteria, the number of water molecules in the first and second shells was calculated by setting to 6.0 and 5.7 Å, respectively, as these values reflect the position of the second minimum of the RDF in each case. The two exponents of the step function in Equation 6 were set to reproduce the integral of the RDF in bulk water within .
Calculation of ion-protein electrostatic interaction energies
Request a detailed protocolTo estimate the electrostatic interaction energy between each K+ or Na+ ion and the channel, we used a continuum electrostatic model based on the linearized form of Poisson’s equation. That is, for each ion k and snapshot , we evaluated the electrostatic potential generated at the position of the ion by the atomic charges in the protein structure (), in the context of a heterogenous dielectric-constant distribution defined by the protein ( = 2 or 1), an implicit membrane surrounding the protein ( = 2), the bulk water solvent ( = 80), the water solvent inside the pore ( = 40), and the ion in question ( = 2 or 1). (Note that both charge and dielectric-constant distributions are a function of , as is the position of the ion.) The interaction between this electrostatic potential and the ion is given by:
These calculations were carried out with the PBEQ numerical solver implemented in CHARMM c44b1 (Brooks et al., 2009). The set of atomic charges employed are those in the CHARMM36 force field; the set of atomic radii used to define the dielectric boundaries derives from an existing optimization for continuum-electrostatic calculations based on the CHARMM27 force field (Best et al., 2012; Nina et al., 1997). Electrostatic potential, charge, and dielectric-constant distributions were discretized on a lattice of 150 × 150 × 150 Å3, with lattice-point spacing of 1 Å. The thickness of the implicit membrane was 28.8 Å. The low-dielectric span of pore water was defined by a cylinder of radius 18 Å and height 28.8 Å. The protein surface was defined using the REEN method, using a probe radius of 1 Å.
Calculation of ion conductance
Request a detailed protocolTo obtain an approximate estimate of the K+ and Na+ conductance of TMEM175 channel, we used the theoretical formulation proposed by Zhu and Hummer, 2012, wherein the conductance is inferred from the free-energy () and diffusion () profiles as a function of the ion position along the pore axis:
where q (=1) is the charge of the permeant ion, C is the bulk ion concentration (100 mM), and S is the cross-sectional area (314 Å2) of the cylindrical region considered to evaluate Equations 3–5. As a first approximation, we assumed the diffusion profile be flat, and used ~ 2 10–9 m2/s for K+ and ~ 1.5 10–9 m2/s for Na+, following previous simulation studies (Zhou et al., 2017; Zhu and Hummer, 2012). Note that depends linearly with , while the dependence on is exponential; thus, while it is very plausible that the diffusion constant will vary by at least a factor of 2 as the ion traverses the pore (Zhou et al., 2017), the order of magnitude of the conductance is largely set by the features of the free-energy profile.
Calculation of free energies of selectivity
Request a detailed protocolTo evaluate free-energy differences between K+ and Na+ states, whether in bulk water or at the isoleucine constriction within the TMEM175 pore, we used the free-energy perturbation method as implemented in NAMD 2.12. To ascertain the magnitude of sampling errors, all transformations were carried out in the forward (K+ to Na+) and backward (K+ to Na+) directions. The transformations were carried out in 26 steps; each step included a 50 ps equilibration, excluded from analysis, followed by 500 ps of averaging time. For bulk water, the coupling parameter lambda was varied in increments of 0.04 in the [0, 0.8] interval and in increments of 0.02 in the [0.8, 1] interval. For the protein calculations, lambda was changed in increments of 0.04 in the [0, 0.96] interval and in increments of 0.02 in the [0.96, 1] interval. To ensure that the K+ and Na+ ions involved in the transformation are in the same position, the distance between the two particles is restrained to zero using a harmonic potential of force constant 2425.58 kcal/(mol Å2). In the protein calculations, the ion(s) are confined to remain within the isoleucine constriction using a flat-bottom distance restraint, defined relative to the center-of-mass of the backbone atoms of two groups of residues, namely 44–45, 269–270 (first group) and 47–49, 272–274 (second group). These two groups also define an axis; the restraint acts on the distance between the ion and the center-of-mass, as projected on that axis. The confining potential is flat-bottomed, permitting fluctuations of ±0.1 Å in the distance, but further deviations are suppressed with a harmonic function of force constant 100 kcal/(mol Å2).
Electrophysiological analysis
Request a detailed protocolElectrophysiological recordings of TMEM175 constructs were performed in HEK293T cells (ATCC CRL-3216). HEK293T cells were recently purchased from ATCC and are negative for mycoplasma contamination. HEK293T cells cultured in DMEM supplemented with 10% FBS. To transfer cells in single dishes, cells were detached by trypsin treatment. The detached cells were transferred to poly-Lys-treated 35 mm single dishes (FluoroDish, World Precision Instruments) and incubated overnight at 37°C in fresh media. Cells in a single dish were transfected with 1.25 μg of c-term EGFP tagged hTMEM175 plasmid using 3.75 μg of PEI 25 k (Polysciences, Inc). Electrophysiological recordings were performed 48–72 hr after transfection. Prior to recording, media was replaced with a bath solution containing 145 mM Na-methanesulfonate (MS), 5 mM NaCl, 1 mM MgCl2, 1 mM CaCl2, 10 mM HEPES/Tris pH 7.4. Ten-cm-long borosilicate glasses were pulled and fire polished (Sutter instrument). Glass pipettes of resistances between 3 and 10 MΩ were filled with a pipette solution containing 150 mM Cs-MS, 5 mM MgCl2, 10 mM EGTA, 10 mM HEPES/Tris pH 7.4, and GΩ seals were formed after gentle suction. The recordings were performed in whole-cell patch clamp configuration using the following protocol: from a holding potential of 0 mV, the voltage was stepped to voltages between −100 and +100 mV, in 20 mV increments. To measure currents reduced by 4-AP or AP-6, bath solutions were perfused with a solution containing 1 mM 4-AP or 2 mM AP-6. DMSO concentration was maintained as 0.33% during the entire experiment. The currents were recorded using Axon Digidata 1550B digitizer and Clampex 10.6 (Molecular Devices, LLC) and analyzed using AxoGraph X 1.7.6 (AxoGraph Scientific). Each experiment was performed in a unique cell. To determine the IC50 of 6-AP, the dose-response curve for TMEM175 current at +100 mV in the presence of AP-6 was fit with the following equation:
where I is the normalized current, Imax and Imin are maximum and minimum normalized current, respectively. [AP-6] is the concentration of AP-6 perfused.
General chemistry
Request a detailed protocolChemical reagents and materials were purchased from Sigma-Aldrich, Thermo-Fisher, and TCI, and used without further purifications. DCM (dichloromethane), MeOH, and n-hexane for column chromatography and recrystallization were used for HPLC grade without additional purifications. Thin layer chromatography (TLC) analysis was performed for reaction monitoring on the pre-coated silica gel 60 F254 glass plates. Both starting materials and the desired product were checked by UV light (254 nm). Flash column chromatography was carried out on silica gel (400–630 mesh) to separate the target molecule.
Synthesis of compound 1 2,2'-(1,3-phenylene)bis(pyridin-4-amine) (AP-6)
Request a detailed protocolAP-6 was prepared through Suzuki-Miyaura cross-coupling reaction with palladium catalysis. In a round-bottom flask with a magnetic bar, a mixture of 2-bromopyridin-4-amine (5 mmol, 865 mg), 1,3-phenylenediboronic acid (0.75 equiv., 622 mg), K2CO3 (2 equiv., 1.38 g), and Pd(OAc)2 (7 mol%, 78.6 mg) were dissolved in H2O:EtOH solution mixture (8 mL: 32 mL), and the solution was stirred at 100°C for 24 hr under air. After completion (monitored by TLC), the reaction mixture was filtered through a Celite (after cooling to room temperature), and then the solid on the filter was washed with EtOAc. The mixture was added to brine and extracted with EtOAc for three times. The combined organic layer was dried over MgSO4, filtered, and concentrated under reduced pressure. The residue was purified by chromatography in silica gel (DCM/MeOH), and purified again with recrystallization (n-hexane/DCM) to give the desired product.
Proton nuclear magnetic resonance (1H NMR) spectra were recorded on Bruker AVANCE 500 (500 MHz). In addition, 13C{1H} NMR was measured on the same machine (125 MHz), and the spectra was fully decoupled with proton by broad band decoupling. Chemical shifts for NMR were quoted in parts per million referenced to the appropriate solvent peak (DMSO in DMSO-d6). The abbreviation codes were adopted to describe 1H NMR peak patterns; d = doublet, t = triplet, and br = broad. Coupling constants, J, were displayed in Hertz unit (Hz). Infrared (IR) spectra were recorded on Bruker Alpha FT-IR spectrometer. High-resolution mass spectra were acquired on a high-resolution Q-TOF mass spectrometer (ionization mode: ESI).
Appendix 1
Data availability
Cryo-EM maps and atomic coordinates have been deposited with the EMDB and PDB under accession codes EMD-26626 and PDB 7UNL for open TMEM175 and codes EMD-26627 and PDB 7UNM for closed TMEM175. Source data have been provided for Figure 5.
References
-
PHENIX: a comprehensive Python-based system for macromolecular structure solutionActa Crystallographica. Section D, Biological Crystallography 66:213–221.https://doi.org/10.1107/S0907444909052925
-
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
-
PLUMED: A portable plugin for free-energy calculations with molecular dynamicsComputer Physics Communications 180:1961–1972.https://doi.org/10.1016/j.cpc.2009.05.011
-
CHARMM: the biomolecular simulation programJournal of Computational Chemistry 30:1545–1614.https://doi.org/10.1002/jcc.21287
-
Metadynamics convergence law in a multidimensional systemPhysical Review. E, Statistical, Nonlinear, and Soft Matter Physics 81:055701.https://doi.org/10.1103/PhysRevE.81.055701
-
The structure of the potassium channel: molecular basis of K+ conduction and selectivityScience (New York, N.Y.) 280:69–77.https://doi.org/10.1126/science.280.5360.69
-
Features and development of CootActa Crystallographica. Section D, Biological Crystallography 66:486–501.https://doi.org/10.1107/S0907444910007493
-
p. M393T Variant as a risk factor for Parkinson DiseaseHuman Molecular Genetics 1:e136.https://doi.org/10.1093/hmg/ddz136
-
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
-
Induced Polarization in Molecular Dynamics Simulations of the 5-HT3 Receptor ChannelJournal of the American Chemical Society 142:9415–9427.https://doi.org/10.1021/jacs.0c02394
-
Genetic, Structural, and Functional Evidence Link TMEM175 to SynucleinopathiesAnnals of Neurology 87:139–153.https://doi.org/10.1002/ana.25629
-
Crystal structure of a mammalian voltage-dependent Shaker family K+ channelScience (New York, N.Y.) 309:897–903.https://doi.org/10.1126/science.1116269
-
A kinetic model of trp-cage folding from multiple biased molecular dynamics simulationsPLOS Computational Biology 5:e1000452.https://doi.org/10.1371/journal.pcbi.1000452
-
Force-Correction Analysis Method for Derivation of Multidimensional Free-Energy Landscapes from Adaptively Biased Replica SimulationsJournal of Chemical Theory and Computation 17:6775–6788.https://doi.org/10.1021/acs.jctc.1c00586
-
The Persistent Question of Potassium Channel Permeation MechanismsJournal of Molecular Biology 433:167002.https://doi.org/10.1016/j.jmb.2021.167002
-
Atomic Radii for Continuum Electrostatics Calculations Based on Molecular Dynamics Free Energy SimulationsThe Journal of Physical Chemistry B 101:5239–5248.https://doi.org/10.1021/jp970736r
-
Conduction and Gating Properties of the TRAAK Channel from Molecular Dynamics Simulations with Different Force FieldsJournal of Chemical Information and Modeling 60:6532–6543.https://doi.org/10.1021/acs.jcim.0c01179
-
Development of the CHARMM Force Field for LipidsThe Journal of Physical Chemistry Letters 2:1526–1532.https://doi.org/10.1021/jz200167q
-
Scalable molecular dynamics with NAMDJournal of Computational Chemistry 26:1781–1802.https://doi.org/10.1002/jcc.20289
-
3D variability analysis: Resolving continuous flexibility and discrete heterogeneity from single particle cryo-EMJournal of Structural Biology 213:107702.https://doi.org/10.1016/j.jsb.2021.107702
-
Efficient reconstruction of complex free energy landscapes by multiple walkers metadynamicsThe Journal of Physical Chemistry. B 110:3533–3539.https://doi.org/10.1021/jp054359r
-
CTFFIND4: Fast and accurate defocus estimation from electron micrographsJournal of Structural Biology 192:216–221.https://doi.org/10.1016/j.jsb.2015.08.008
-
HOLE: A program for the analysis of the pore dimensions of ion channel structural modelsJournal of Molecular Graphics 14:354–360.https://doi.org/10.1016/s0263-7855(97)00009-x
-
GRIFFIN: A versatile methodology for optimization of protein-lipid interfaces for membrane protein simulationsJournal of Chemical Theory and Computation 7:1167–1176.https://doi.org/10.1021/ct100576m
-
Improvement of Cryo-EM Maps by Density ModificationBiochemistry 1:845032.https://doi.org/10.1101/845032
-
A voltage-dependent K+ channel in the lysosome is required for refilling lysosomal Ca2+ storesJournal of Cell Biology 216:1715–1730.https://doi.org/10.1083/jcb.201612123
-
Theory and simulation of ion conduction in the pentameric GLIC channelJournal of Chemical Theory and Computation 8:3759–3768.https://doi.org/10.1021/ct2009279
Article and author information
Author details
Funding
National Cancer Institute (P30 CA008748)
- Richard K Hite
National Institute of General Medical Sciences (R01-GM141553)
- Richard K Hite
Josie Robertson Investigators Program
- Richard K Hite
Searle Scholars Program
- Richard K Hite
Ministry of Education (2019H1A2A1076014)
- Jooyeon Lee
National Heart, Lung, and Blood Institute (Division of Intramural Research)
- José D Faraldo-Gómez
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank the MSKCC HPC group for assistance with data processing and the members of the labs for comments on the manuscript. We also thank Dr Rahul Banerjee for early simulation studies of closed state TMEM175 not included in this manuscript. This work was supported by NIH-NCI Cancer Center Support Grant P30 CA008748, NIGMS R01-GM141553 (RKH), the Josie Robertson Investigators Program (RKH), the Searle Scholars Program (RKH), the NRF Global PhD. Fellowship program funded by the Republic of Korea Ministry of Education 2019H1A2A1076014 (JL) and the Division of Intramural Research of NHLBI-NIH (JDFG). Computational resources were in part provided by the NIH supercomputing center (Biowulf).
Copyright
This is an open-access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication.
Metrics
-
- 1,922
- views
-
- 464
- downloads
-
- 12
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Neuroscience
- Structural Biology and Molecular Biophysics
More than 20 recurrent missense gain-of-function (GOF) mutations have been identified in the sodium-activated potassium (KNa) channel gene KCNT1 in patients with severe developmental and epileptic encephalopathies (DEEs), most of which are resistant to current therapies. Defining the neuron types most vulnerable to KCNT1 GOF will advance our understanding of disease mechanisms and provide refined targets for precision therapy efforts. Here, we assessed the effects of heterozygous expression of a Kcnt1 GOF variant (Kcnt1Y777H) on KNa currents and neuronal physiology among cortical glutamatergic and GABAergic neurons in mice, including those expressing vasoactive intestinal polypeptide (VIP), somatostatin (SST), and parvalbumin (PV), to identify and model the pathogenic mechanisms of autosomal dominant KCNT1 GOF variants in DEEs. Although the Kcnt1Y777H variant had no effects on glutamatergic or VIP neuron function, it increased subthreshold KNa currents in both SST and PV neurons but with opposite effects on neuronal output; SST neurons became hypoexcitable with a higher rheobase current and lower action potential (AP) firing frequency, whereas PV neurons became hyperexcitable with a lower rheobase current and higher AP firing frequency. Further neurophysiological and computational modeling experiments showed that the differential effects of the Kcnt1Y777H variant on SST and PV neurons are not likely due to inherent differences in these neuron types, but to an increased persistent sodium current in PV, but not SST, neurons. The Kcnt1Y777H variant also increased excitatory input onto, and chemical and electrical synaptic connectivity between, SST neurons. Together, these data suggest differential pathogenic mechanisms, both direct and compensatory, contribute to disease phenotypes, and provide a salient example of how a pathogenic ion channel variant can cause opposite functional effects in closely related neuron subtypes due to interactions with other ionic conductances.
-
- Structural Biology and Molecular Biophysics
Polycystin-1 (PC1) is the protein product of the PKD1 gene whose mutation causes autosomal dominant Polycystic Kidney Disease (ADPKD). PC1 is an atypical G protein-coupled receptor (GPCR) with an autocatalytic GAIN domain that cleaves PC1 into extracellular N-terminal and membrane-embedded C-terminal (CTF) fragments. Recently, activation of PC1 CTF signaling was shown to be regulated by a stalk tethered agonist (TA), resembling the mechanism observed for adhesion GPCRs. Here, synthetic peptides of the first 9- (p9), 17- (p17), and 21-residues (p21) of the PC1 stalk TA were shown to re-activate signaling by a stalkless CTF mutant in human cell culture assays. Novel Peptide Gaussian accelerated molecular dynamics (Pep-GaMD) simulations elucidated binding conformations of p9, p17, and p21 and revealed multiple specific binding regions to the stalkless CTF. Peptide agonists binding to the TOP domain of PC1 induced close TOP-putative pore loop interactions, a characteristic feature of stalk TA-mediated PC1 CTF activation. Additional sequence coevolution analyses showed the peptide binding regions were consistent with covarying residue pairs identified between the TOP domain and the stalk TA. These insights into the structural dynamic mechanism of PC1 activation by TA peptide agonists provide an in-depth understanding that will facilitate the development of therapeutics targeting PC1 for ADPKD treatment.