Effective pore size and radius of capture for K(+) ions in K-channels.

Reconciling protein functional data with crystal structure is arduous because rare conformations or crystallization artifacts occur. Here we present a tool to validate the dimensions of open pore structures of potassium-selective ion channels. We used freely available algorithms to calculate the molecular contour of the pore to determine the effective internal pore radius (r(E)) in several K-channel crystal structures. r(E) was operationally defined as the radius of the biggest sphere able to enter the pore from the cytosolic side. We obtained consistent r(E) estimates for MthK and Kv1.2/2.1 structures, with r(E) = 5.3-5.9 Å and r(E) = 4.5-5.2 Å, respectively. We compared these structural estimates with functional assessments of the internal mouth radii of capture (r(C)) for two electrophysiological counterparts, the large conductance calcium activated K-channel (r(C) = 2.2 Å) and the Shaker Kv-channel (r(C) = 0.8 Å), for MthK and Kv1.2/2.1 structures, respectively. Calculating the difference between r(E) and r(C), produced consistent size radii of 3.1-3.7 Å and 3.6-4.4 Å for hydrated K(+) ions. These hydrated K(+) estimates harmonize with others obtained with diverse experimental and theoretical methods. Thus, these findings validate MthK and the Kv1.2/2.1 structures as templates for open BK and Kv-channels, respectively.

K-channel crystal structures reveal that they are machines optimized for efficient and selective K+ ion transport. Thanks to the groundbreaking work from the Mackinnon lab, we know the structure of several K-channels in detail1. A great deal of effort has been invested in building a conceptual scaffold that makes functional sense of these K-channel structures. In the pore domain of K-channels, this scaffold works pretty well for understanding toxin binding and K+ selectivity. This is partly because the external vestibule and the selectivity filters are experimentally more accessible, and also because of the low variance in atom’s space coordinates among different K-channels across diverse crystallization conditions23456. Thus, the structure of the external vestibule together with that of the selectivity filter enjoy a solid functional reputation. By contrast, the internal vestibule of the pore seems to be much less well-defined. On one hand, the dimensions of the cytosolic aspects of the pore structure differ from channel to channel in the various crystal structures3578. On the other, the internal vestibule seems to be the flexible part of the protein where the voltage controlled gate is located. For example, the structures of KcsA (PDB:1K4C) and Slo2.2 (PDB:5A6E) seem to correspond to channels crystallized in the closed conformation because their pore’s internal opening size is smaller than that of a hydrated K+ 59. Meanwhile, the structures of the bacterial MthK (PDB: 4HYO) and of the mammalian Kv1.2/2.1 paddle chimera (PDB: 2R9R) appear to be those of open channels, with the MthK internal vestibule being ~10 Å wider than that of Kv1.2/2.1371011.

The size difference in the pore internal vestibule in the structure of the MthK vs. the Kv1.2/2.1 chimera structure was somehow electrophysiologically corroborated in their functional counterparts, the mouse large conductance calcium and voltage gated K+-channel (BK) and the Drosophila voltage gated Shaker K+-channel, respectively. Cysteine substitution scanning accessibility experiments1213 in Shaker and BK channels, and also side-chain volume changes in residues located at the internal entrance of BK14, suggested a pore several angstrom wider than Shaker’s. Nevertheless in diffusional determinations of their radii of capture, BK is only ~1.4 Å wider (see below)1516. To what extent do these differences represent structurally different K-channels, or do they just reflect diverse conformations? Are they the result of rarely visited conformational states or are they caused by crystallization artifacts? These questions are especially important when, based on sequence homology or functional properties such as single channel conductance and pharmacology, we use the atomic coordinates of one channel as a structural template for a distantly related one, as is the case when using the structure of MthK to model the BK channel.

Here, we present a simple methodology to estimate the functional dimensions of the internal entrance to the pore in K+-channel structural models. Briefly: for several K+-channel crystal structures we first estimated their effective pore size17 and then compared these values with functional estimations of their diffusional radius of capture obtained from diffusion limited unitary currents1516181920. To validate the methodology, the difference between these two estimates must equal robust approximations of the hydrodynamic radius of K+ ions.

The Effective Pore Size

The effective dimension for permeation in a K+-channel pore is determined by the effective size of the K+ ion, which, in turn, is determined by the average number of water molecules forming its hydration shell. At one extreme, if K+ ion coming from the cytosolic side were a point charge, having no physical dimension and no hydration shell, and the protein were a rigid structure, the effective sectional area to pass should be delimited by the Van der Waals envelope of the pore entrance21. But if instead K+ were a solid spherical body of finite dimensions, it will not pass across surface tortuosities smaller than its own size, resulting in a less rugose, and decreased, sectional area available for permeation. In this other extreme, if the hydrated K+-complex is larger than the pore entrance, it will not enter. In this latter case, from the ion’s perspective, the pore is closed. Thus, if we assume the ion to be a rigid spherical probe of variable size moving into a circular pore, we define the effective pore radius (rE) as the radius of the largest sphere that is able to pass through the pore17.

To estimate rE on potassium channel structures we used freely available algorithms (SURF, MSMS, and HOLLOW). These routines calculate the molecular surface of proteins left by the contact of rigid spherical test probes of variable radius. In SURF and MSMS the probe is rolled on top of the Van der Waals surface of the channel structure212223 (Fig. 1 in21). When the probe contacts simultaneously more than one atom, these contacts produce a series of disconnected patches that are filled with concave surfaces having the probe´s curvature. Thus, the resulting envelope is a molecular surface that is less detailed as the probe grows bigger (See for example first column in Fig. 1). On the other hand, HOLLOW consisted in filling the protein holes and voids with virtual overlapping oxygen atoms placed at fixed intervals defined in a 0.2–0.5 Å cubic grid24. The molecular surface is then defined by the envelope of the virtual oxygen atoms that are not contacted by the probe (See for example second column in Fig. 1).

Figure 1

The effective opening radius in two K-channel structures.

Shown are ~5 Å slabs of the molecular surface of the large conductance bacterial channel MthK (PDB: 4HYO) and the chimeric Kv1.2/2.1 voltage gated K-channel (PDB: 2R9R). The molecular surfaces were calculated as indicated in the main text. The probe radius (in Å) is shown above each representation. The top row shows the molecular surfaces for a 4 Å-radius probe (about the size of a hydrated K+). The middle row (with underlined radii values) shows the size of the largest spheres able to pass through the pore (rE). The lower row shows that probes 0.1 Å bigger cannot pass. With SURF, these probes leave a dimple or a bump at the pore entrance, depending on whether the rolling trajectory began outside or inside the cavity, respectively (radius = 5.7 Å and 4.6 Å, for MthK and Kv1.2/2.1, respectively). Meanwhile, with HOLLOW, the non-permeating probes (radius = 5.9 Å and 5.3 Å, for MthK and Kv1.2/2.1, respectively) leave the cavity full of virtual O atoms (red spheres). Two opposite pore-helices and K+ ions in selectivity filter (in green) are shown for reference. The T1, the β-subunit, and the voltage sensing domains are not represented in 2R9R for clarity. Figure prepared with VMD (http://www.ks.uiuc.edu/Research/vmd/).

Examples of the molecular surfaces left by SURF and HOLLOW (Fig. 1) show that smaller probes go inside the cavity defining a detailed molecular surface of the pore walls, whereas bigger probes leave poorer surface detail. Probes big enough will not be able to enter, transforming the pore into another dimple in the protein surface. When this condition is satisfied (Bottom row), we state a 0.1 Å smaller radius as the biggest probe able to enter, defining the effective pore radius, rE17. This 0.1 Å dissimilarity makes a dramatic difference in the probe’s ability to enter into the inner vestibule (compare the middle with the bottom row).

We measured rE in all crystallographic structures available today of four different K+-channels: Kv1.2/2.1 paddle-chimera, KvAP, MthK and KcsA channels (Table 1). For Kv 1.2/2.1, SURF, MSMS and HOLLOW gave “average” values of rE = 4.5 Å, rE = 4.5 Å, and rE = 5.2 Å, respectively. For KvAP, rE values were less regular but, in general agreed with those of Kv1.2/2.1. Also, for MthK, rE = 5.3–5.6 Å, rE = 5.3–5.5 Å, and rE = 5.7–6.0 Å, for SURF, MSMS, and HOLLOW, respectively. For KcsA rE values were much smaller, as expected for a closed channel.

Table 1

Effective pore radii (r E) for several K+ channel crystals.

Protein	PDB code	Crystal conditions & resolution	Sphere radius (Å)						Ref.
			hollow		Surf		msms
			rE	rK	rE	rK	rE	rK
Kv1.2/2.1 Paddle chimera	2R9R	Complex with lipids at 2.4 Å	5.2	4.4	4.5	3.7	4.4	3.6	7
	4JTA	Complex with ChTx at 2.5 Å	5.2	4.4	4.5	3.7	4.6	3.8	2
KvAP	2A0L	3.9 Å	4.4	3.6	4.1	3.3	4.2	3.4	36
	1ORQ	3.2 Å	4.9	4.1	5.1	4.3	5.5	4.7	37
MthK	4HYO	1.65 Å	5.9	3.7	5.6	3.4	5.5	3.3	10
	3LDC	1.45 Å	5.7	3.5	5.3	3.1	5.3	3.1	38
	1LNQ	3.3 Å***	5.8	3.6	5.6	3.4	5.5	3.3	3
KcsA	1K4C	2.0 Å	2.1		<1.4		<1.4		4
	1BL8	3.2 Å	<1.4		<1.4		<1.4		5

Surface calculation algorithms were used to define the cast envelope left by a spherical probe of varying radius222324. A probe able to enter the pore from the channel’s cytosolic side created a continuous surface connecting the cytosolic face with the pore walls. However, when the probe cannot longer enter, it casts a bump in the place of the pore. Using this criterion, rE was defined as the largest sphere able to enter the pore. rK was calculated from Eq. (3). For Kv1.2/2.1 and KvAP we used rC = 0.8 Å16, while rC = 2.2 Å was used for MthK channels15. ChTx stand for charybdotoxin. ***The biological assembly of this structure is not fully symmetric at the pore entrance.

The Radius of Capture

We compared rE with functional estimates of radius of capture (rC) in the two counterparts of the MthK and Kv1.2/2.1 structures: the BK, and the Shaker channels, respectively1516 (To our knowledge, these are the only rC estimates available for K+-channels). rC is an assessment of the dimension of the pore’s entrance obtained from single channel current measurements in a regime in which diffusion of K+ ions into the pore entrance is the rate limiting step for passage through the channel. If the transmembrane voltage and/or the recording solution viscosity are high enough, as the applied voltage increases the unitary ionic current departs significantly from the Ohms law; approaching a plateau. This asymptotic saturation amplitude represents the diffusion limited rate with which ions approach the channel entrance. Thus, the bigger the pore entrance, or the higher the ion concentration, or the larger the diffusion coefficient, the higher is the limiting asymptotic current. If the pore is assumed to be a hemispheric sink into which approaching ions vanish, it is possible from diffusional collisions theory to infer its dimension from the limited unitary current, iDL, according to1516181920.

where z is the permeant ion valence, eo the elementary charge, c the bulk ion concentration, and D its diffusion coefficient. If the permeant ion is a point charge, rC represents the average radius of the pore’s entrance1518. However, because ions have finite size and pores usually aren’t circular, shape and size must be considered. For simplicity, let’s assume that ions are solid spheres reaching a circular pore with radius rO, then,

where rK is the ion´s spherical radius151820. For point charges, rO = rC. However, for a given, experimentally obtained rC, estimates of rO will grows linearly with rK, because the number of effective collisions will be indistinguishable from those produced by point charges having identical radial deviations from the central trajectory16.

To reach asymptotic currents values at experimentally possible voltages (between 100–350 mV), 2 molar sucrose was added to the cytosolic recording solution to increase viscosity by ~7-fold, and values of rC = 2.2 Å and rC = 0.8 Å, were obtained for BK and Shaker internal entrances, respectively1516. According to Eq. 2, by knowing rK, we could obtain rO, their vestibule average radius; nevertheless, the hydrated ion size is not well defined, in part because hydration is a fuzzy arrangement of water molecules bound with dissimilar energies and lifetimes151625.

To compare structural (rE) with functional (rC) estimates of the internal pore vestibule, here we make rO = rE, then:

To be meaningful, Eq. 3 must provide consistent estimates for rK. In fact, Table 1 show that this was true across all structural models and channel types, regardless of the surface calculation algorithm used. Although, HOLLOW gave consistently larger estimates, rK fell mostly in the 3.5–4.0 Å interval, providing a robust validation for pore openings dimensions. Strikingly, the values for rK are well in agreement with the hydrated size of K+ ions estimates obtained using other very different experimental and theoretical methods2526272829. A 3.5–4 Å radius sphere occupies a volume of 180–270 Å3 which is comparable to the anti-prismatic geometric arrangement of waters co-crystalized “in flagrante” with K+ in the KcsA structure (~300 Å3)430, and in agreement with several studies, suggest a coordination number of 8 or less25262729. Of course, our intention is not to provide another estimate of the hydrodynamic size of hydrated K+ but to show the robustness of Eq. 3 in comparing structural and functional data.

We should bear in mind that the rC estimates are restricted to the hydrated path to the pore entrance, and may contain two sources of error: on the one hand, an overestimation of rC arises with the presence of the gating ring in BK and the tetramerization domain in Shaker. In typical recording solutions, both contribute modestly to single channel conductance (up to 30%) by increasing the local K+ concentration in corresponding proportion93132. On the other hand, a narrowing of the hydrated path promoted by the hydrophobic nature of the channel’s inner cavity residues may cause an underestimation of the sectional access to the pore; nevertheless at very positive voltages the water content in the cavity matches the cavity’s volume33. In this regard, it is interesting to note that our estimates of rE for Kv1.2/2.1 chimera are in the 4.4–5.2 Å range, just enough to allow K+ hydration to occur in this confined space as in bulk solution, according to molecular dynamics calculations on Shaker34. In narrower pores, hydration number decreases and the energetic cost of putting ions inside may be too high to make permeation possible34. In contrast, a wider cavity as MthK with rE = 5.3–5.9 Å, may provide the incoming hydrated K+ with a second, loosely attached hydration shell that may reduce the ion’s energy and/or the friction of its rigid hydration cage163435.

The rC estimates shown here stem from studies in liquid phase, at room temperature, away from equilibrium, and with fully operative friction forces151620. By contrast, K+ channel structures are obtained at low temperature, in solid phase, and in equilibrium. Thus, the agreement between these functional and structural data is satisfying. They gave consistent estimations of the hydrodynamic radius of K+ ions in physiological solutions and, while we cannot ascertain the degree of opening of KvAP structures, we can propose that KcsA structures are closed, whereas Kv1.2/2.1, KvAP, and MthK are bona fide open channels having different pore dimensions, possibly explaining differences in unitary conductance16.

Additional Information

How to cite this article: Moldenhauer, H. et al. Effective pore size and radius of capture for K+ ions in K-channels. Sci. Rep. 6, 19893; doi: 10.1038/srep19893 (2016).