A computational study of barium blockades in the KcsA potassium channel based on multi-ion potential of mean force calculations and free energy perturbation

Probing Ion Configurations in the KcsA Selectivity Filter with Single-Isotope Labels and 2D IR Spectroscopy

The potassium ion (K+) configurations of the selectivity filter of the KcsA ion channel protein are investigated with two-dimensional infrared (2D IR) spectroscopy of amide I vibrations. Single 13C-18O isotope labels are used, for the first time, to selectively probe the S1/S2 or S2/S3 binding sites in the selectivity filter. These binding sites have the largest differences in ion occupancy in two competing K+ transport mechanisms: soft-knock and hard-knock. According to the former, water molecules alternate between K+ ions in the selectivity filter while the latter assumes that K+ ions occupy the adjacent sites. Molecular dynamics simulations and computational spectroscopy are employed to interpret experimental 2D IR spectra. We find that in the closed conductive state of the KcsA channel, K+ ions do not occupy adjacent binding sites. The experimental data is consistent with simulated 2D IR spectra of soft-knock ion configurations. In contrast, the simulated spectra for the hard-knock ion configurations do not reproduce the experimental results. 2D IR spectra of the hard-knock mechanism have lower frequencies, homogeneous 2D lineshapes, and multiple peaks. In contrast, ion configurations of the soft-knock model produce 2D IR spectra with a single peak at a higher frequency and inhomogeneous lineshape. We conclude that under equilibrium conditions, in the absence of transmembrane voltage, both water and K+ ions occupy the selectivity filter of the KcsA channel in the closed conductive state. The ion configuration is central to the mechanism of ion transport through potassium channels.

Acidic Polysaccharides as Green Alternatives for Barite Scale Dissolution

Barium sulfate (barite) scale poses significant challenges for processes ranging from water treatment to fossil fuel production. Here, we identify alginate (a polysaccharide derived from brown algae) as a potent, "green" alternative to commercial barite demineralizing agents. Unlike conventional treatments of inorganic scales that require caustic conditions, alginate polymers dissolve barite at near-neutral conditions. In this study, we benchmark the demineralizing efficacy of alginate against a commercial dissolver, diethylenetriaminepentaacetic acid (DTPA), using a combination of bulk dissolution assays, scanning probe microscopy, and molecular dynamics simulations. Time-resolved rates of dissolution measured in a microfluidic device show that demineralization is enhanced more than an order of magnitude under flow. In situ atomic force microscopy reveals that alginate and DTPA exhibit distinct mechanisms of surface dissolution; and surprisingly, their binary combination in alkaline media results in a synergistic cooperativity that enhances the overall rate of barite dissolution. These studies collectively demonstrate a unique approach to demineralization using an inexpensive and abundant biopolymer that enables environmentally friendly treatment of inorganic scales.

Open and Closed Structures of a Barium-Blocked Potassium Channel

Barium (Ba²⁺) is a classic permeant blocker of potassium (K⁺) channels. The "external lock-in effect" in barium block experiments, whereby the binding of external K⁺ impedes the forward translocation of the blocker, provides a powerful avenue to investigate the selectivity of the binding sites along the pore of potassium channels. Barium block experiments show that the external lock-in site is highly selective for K⁺ over Na⁺. Wild-type KcsA was crystallized in low K⁺ conditions, and the crystals were soaked in solutions containing various concentrations of barium. Structural analysis reveals open and closed gate conformations of the KcsA channel. Anomalous diffraction experiments show that Ba²⁺ primarily binds to the innermost site S4 of the selectivity filter of the open-gate conformation and also the site S2, but no binding is detected with the closed-gate conformation. Alchemical free-energy perturbation calculations indicate that the presence of a Ba²⁺ ion in the selectivity filter boosts the specificity of K⁺ binding relative to Na⁺ in the external sites S0-S2.

Selectivity of the KcsA potassium channel: Analysis and computation

Ion channels regulate the flux of ions through cell membranes and play significant roles in many physiological functions. Most of the existing literature focuses on computational approaches based on molecular dynamics simulation or numerical solution of the modified Poisson-Nernst-Planck (PNP) system. In this paper, we present an analytical and computational study of a mathematical model of the KcsA potassium channel, including the effects of ion size (Bikerman model) and solvation energy (Born model). Under equilibrium conditions, we obtain an analytical solution of our modified PNP system, which is used to explain selectivity of KcsA of various ions (K^{+}, Na^{+}, Cl^{-}, Ca^{2+}, and Ba^{2+}) due to negative permanent charges inside the filter region and the effect of ion sizes. Our results show that K^{+} is always selected over Na^{+}, as smaller Na^{+} ions have larger solvation energy. As the amount of negative charges in the filter exceeds a critical value, divalent ions (Ca^{2+} and Ba^{2+}) can enter the filter region and block the KcsA channel. For the nonequilibrium cases, due to difficulties associated with a pure analytical or numerical approach, we use a hybrid analytical-numerical method to solve the modified PNP system. Our predictions of selectivity of KcsA channels and saturation phenomenon of the current-voltage (I-V) curve agree with experimental observations.

Computational Studies of Molecular Permeation through Connexin26 Channels

A signal property of connexin channels is the ability to mediate selective diffusive movement of molecules through plasma membrane(s), but the energetics and determinants of molecular movement through these channels have yet to be understood. Different connexin channels have distinct molecular selectivities that cannot be explained simply on the basis of size or charge of the permeants. To gain insight into the forces and interactions that underlie selective molecular permeation, we investigated the energetics of two uncharged derivatized sugars, one permeable and one impermeable, through a validated connexin26 (Cx26) channel structural model, using molecular dynamics and associated analytic tools. The system is a Cx26 channel equilibrated in explicit membrane/solvent, shown by Brownian dynamics to reproduce key conductance characteristics of the native channel. The results are consistent with the known difference in permeability to each molecule. The energetic barriers extend through most of the pore length, rather than being highly localized as in ion-specific channels. There is little evidence for binding within the pore. Force decomposition reveals how, for each tested molecule, interactions with water and the Cx26 protein vary over the length of the pore and reveals a significant contribution from hydrogen bonding and interaction with K(+). The flexibility of the pore width varies along its length, and the tested molecules have differential effects on pore width as they pass through. Potential sites of interaction within the pore are defined for each molecule. The results suggest that for the tested molecules, differences in hydrogen bonding and entropic factors arising from permeant flexibility substantially contribute to the energetics of permeation. This work highlights factors involved in selective molecular permeation that differ from those that define selectivity among atomic ions.

Instantaneous ion configurations in the K+ ion channel selectivity filter revealed by 2D IR spectroscopy

Potassium channels are responsible for the selective permeation of K+ ions across cell membranes. K+ ions permeate in single file through the selectivity filter, a narrow pore lined by backbone carbonyls that compose four K+ binding sites. Here, we report on the two-dimensional infrared (2D IR) spectra of a semisynthetic KcsA channel with site-specific heavy (13C18O) isotope labels in the selectivity filter. The ultrafast time resolution of 2D IR spectroscopy provides an instantaneous snapshot of the multi-ion configurations and structural distributions that occur spontaneously in the filter. Two elongated features are resolved, revealing the statistical weighting of two structural conformations. The spectra are reproduced by molecular dynamics simulations of structures with water separating two K+ ions in the binding sites, ruling out configurations with ions occupying adjacent sites.

Multi-ion free energy landscapes underscore the microscopic mechanism of ion selectivity in the KcsA channel

Potassium (K(+)) channels are transmembrane proteins that passively and selectively allow K(+) ions to flow through them, after opening in response to an external stimulus. One of the most critical functional aspects of their function is their ability to remain very selective for K(+) over Na(+) while allowing high-throughput ion conduction at a rate close to the diffusion limit. Classically, it is assumed that the free energy difference between K(+) and Na(+) in the pore relative to the bulk solution is the critical quantity at the origin of selectivity. This is the thermodynamic view of ion selectivity. An alternative view assumes that kinetic factors play the dominant role. Recent results from a number of studies have also highlighted the great importance of the multi-ion single file on the selectivity of K(+) channels. The data indicate that having multiple K(+) ions bound simultaneously is required for selective K(+) conduction, and that a reduction in the number of bound K(+) ions destroys the multi-ion selectivity mechanism utilized by K(+) channels. In the present study, multi-ion potential of mean force molecular dynamics computations are carried out to clarify the mechanism of ion selectivity in the KcsA channel. The computations show that the multi-ion character of the permeation process is a critical element for establishing the selective ion conductivity through K(+)-channels. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.

Octa-Coordination and the Aqueous Ba(2+) Ion

The hydration structure of Ba(2+) ion is important for understanding blocking mechanisms in potassium ion channels. Here, we combine statistical mechanical theory, ab initio molecular dynamics simulations, and electronic structure methods to calculate the hydration free energy and local hydration structure of Ba(2+)(aq). The predicted hydration free energy (-304 ± 1 kcal/mol) agrees with the experimental value (-303 kcal/mol) when a maximally occupied, unimodal inner solvation shell is treated. In the local environment defined by the first shell of hydrating waters, Ba(2+) is directly and stably coordinated by eight (8) waters. Octa-coordination resembles the crystal structure of Ba(2+) and K(+) bound in potassium ion channels, but differs from the local hydration structure of K(+)(aq) determined earlier.

Ionic interactions of Ba2+ blockades in the MthK K+ channel

The movement and interaction of multiple ions passing through in single file underlie various fundamental K(+) channel properties, from the effective conduction of K(+) ions to channel blockade by Ba(2+) ions. In this study, we used single-channel electrophysiology and x-ray crystallography to probe the interactions of Ba(2+) with permeant ions within the ion conduction pathway of the MthK K(+) channel. We found that, as typical of K(+) channels, the MthK channel was blocked by Ba(2+) at the internal side, and the Ba(2+)-blocking effect was enhanced by external K(+). We also obtained crystal structures of the MthK K(+) channel pore in both Ba(2+)-Na(+) and Ba(2+)-K(+) environments. In the Ba(2+)-Na(+) environment, we found that a single Ba(2+) ion remained bound in the selectivity filter, preferably at site 2, whereas in the Ba(2+)-K(+) environment, Ba(2+) ions were predominantly distributed between sites 3 and 4. These ionic configurations are remarkably consistent with the functional studies and identify a molecular basis for Ba(2+) blockade of K(+) channels.

Modeling ion channels: past, present, and future

Ion channels are membrane-bound enzymes whose catalytic sites are ion-conducting pores that open and close (gate) in response to specific environmental stimuli. Ion channels are important contributors to cell signaling and homeostasis. Our current understanding of gating is the product of 60 plus years of voltage-clamp recording augmented by intervention in the form of environmental, chemical, and mutational perturbations. The need for good phenomenological models of gating has evolved in parallel with the sophistication of experimental technique. The goal of modeling is to develop realistic schemes that not only describe data, but also accurately reflect mechanisms of action. This review covers three areas that have contributed to the understanding of ion channels: traditional Eyring kinetic theory, molecular dynamics analysis, and statistical thermodynamics. Although the primary emphasis is on voltage-dependent channels, the methods discussed here are easily generalized to other stimuli and could be applied to any ion channel and indeed any macromolecule.

Clustering of the K+ channel GORK of Arabidopsis parallels its gating by extracellular K+

GORK is the only outward-rectifying Kv-like K(+) channel expressed in guard cells. Its activity is tightly regulated to facilitate K(+) efflux for stomatal closure and is elevated in ABA in parallel with suppression of the activity of the inward-rectifying K(+) channel KAT1. Whereas the population of KAT1 is subject to regulated traffic to and from the plasma membrane, nothing is known about GORK, its distribution and traffic in vivo. We have used transformations with fluorescently-tagged GORK to explore its characteristics in tobacco epidermis and Arabidopsis guard cells. These studies showed that GORK assembles in puncta that reversibly dissociated as a function of the external K(+) concentration. Puncta dissociation parallelled the gating dependence of GORK, the speed of response consistent with the rapidity of channel gating response to changes in the external ionic conditions. Dissociation was also suppressed by the K(+) channel blocker Ba(2+) . By contrast, confocal and protein biochemical analysis failed to uncover substantial exo- and endocytotic traffic of the channel. Gating of GORK is displaced to more positive voltages with external K(+) , a characteristic that ensures the channel facilitates only K(+) efflux regardless of the external cation concentration. GORK conductance is also enhanced by external K(+) above 1 mm. We suggest that GORK clustering in puncta is related to its gating and conductance, and reflects associated conformational changes and (de)stabilisation of the channel protein, possibly as a platform for transmission and coordination of channel gating in response to external K(+) .