Filter flexibility and distortion in a bacterial inward rectifier K+ channel: simulation studies of KirBac1.1

Mimicking classical noise in ion channels by quantum decoherence

The mechanism of selectivity in ion channels is still an open question in biology. Recent studies suggest that the selectivity filter may exhibit quantum coherence, which could help explain how ions are selected and conducted. However, environmental noise causes decoherence and loss of quantum effects. It is hoped that the effect of classical noise on ion channels can be modeled using the framework provided by quantum decoherence theory. In this paper, the behavior of the ion channel system was simulated using two models: the Spin-Boson model and the stochastic Hamiltonian model under classical noise. Additionally, using a different approach, the system's evolution was modeled as a two-level Spin-Boson model with tunneling, interacting with a bath of harmonic oscillators, based on decoherence theory. We investigated under what conditions the decoherence model approaches and deviates from the noise model. Specifically, we examined Gaussian noise and Ornstein-Uhlenbeck noise in our model. Gaussian noise shows a very good agreement with the decoherence model. By examining the results, it was found that the Spin-Boson model at a high hopping rate of potassium ions can simulate the behavior of the system in the classical noise approach for Gaussian noise.

Mutational Insight into Allosteric Regulation of Kir Channel Activity

Potassium (K⁺) channels are regulated in part by allosteric communication between the helical bundle crossing, or inner gate, and the selectivity filter, or outer gate. This network is triggered by gating stimuli. In concert, there is an allosteric network which is a conjugated set of interactions which correlate long-range structural rearrangements necessary for channel function. Inward-rectifier K⁺ (Kir) channels favor inward K⁺ conductance, are ligand-gated, and help establish resting membrane potentials. KirBac1.1 is a bacterial Kir (KirBac) channel homologous to human Kir (hKir) channels. Additionally, KirBac1.1 is gated by the anionic phospholipid ligand phosphatidylglycerol (PG). In this study, we use site-directed mutagenesis to investigate residues involved in the KirBac1.1 gating mechanism and allosteric network we previously proposed using detailed solid-state NMR (SSNMR) measurements. Using fluorescence-based K⁺ and sodium (Na⁺) flux assays, we identified channel mutants with impaired function that do not alter selectivity of the channel. In tandem, we performed coarse grain molecular dynamics simulations, observing changes in PG-KirBac1.1 interactions correlated with mutant channel activity and contacts between the two transmembrane helices and pore helix tied to this behavior. Lipid affinity is closely tied to the proximity of two tryptophan residues on neighboring subunits which lure anionic lipids to a cationic pocket formed by a cluster of arginine residues. Thus, these simulations establish a structural and functional basis for the role of each mutated site in the proposed allosteric network. The experimental and simulated data provide insight into key functional residues involved in gating and lipid allostery of K⁺ channels. Our findings also have direct implications on the physiology of hKir channels due to conservation of many of the residues identified in this work from KirBac1.1.

Start Me Up: How Can Surrounding Gangliosides Affect Sodium-Potassium ATPase Activity and Steer towards Pathological Ion Imbalance in Neurons?

Gangliosides, amphiphilic glycosphingolipids, tend to associate laterally with other membrane constituents and undergo extensive interactions with membrane proteins in cis or trans configurations. Studies of human diseases resulting from mutations in the ganglioside biosynthesis pathway and research on transgenic mice with the same mutations implicate gangliosides in the pathogenesis of epilepsy. Gangliosides are reported to affect the activity of the Na⁺/K⁺-ATPase, the ubiquitously expressed plasma membrane pump responsible for the stabilization of the resting membrane potential by hyperpolarization, firing up the action potential and ion homeostasis. Impaired Na⁺/K⁺-ATPase activity has also been hypothesized to cause seizures by several mechanisms. In this review we present different epileptic phenotypes that are caused by impaired activity of Na⁺/K⁺-ATPase or changed membrane ganglioside composition. We further discuss how gangliosides may influence Na⁺/K⁺-ATPase activity by acting as lipid sorting machinery providing the optimal stage for Na⁺/K⁺-ATPase function. By establishing a distinct lipid environment, together with other membrane lipids, gangliosides possibly modulate Na⁺/K⁺-ATPase activity and aid in “starting up” and “turning off” this vital pump. Therefore, structural changes of neuronal membranes caused by altered ganglioside composition can be a contributing factor leading to aberrant Na⁺/K⁺-ATPase activity and ion imbalance priming neurons for pathological firing.

Leveraging intrinsic flexibility to engineer enhanced enzyme catalytic activity

Dynamic motions of enzymes occurring on a broad range of timescales play a pivotal role in all steps of the reaction pathway, including substrate binding, catalysis, and product release. However, it is unknown whether structural information related to conformational flexibility can be exploited for the directed evolution of enzymes with higher catalytic activity. Here, we show that mutagenesis of residues exclusively located at flexible regions distal to the active site of Homo sapiens kynureninase (HsKYNase) resulted in the isolation of a variant (BF-HsKYNase) in which the rate of the chemical step toward kynurenine was increased by 45-fold. Mechanistic pre–steady-state kinetic analysis of the wild type and the evolved enzyme shed light on the underlying effects of distal mutations (>10 Å from the active site) on the rate-limiting step of the catalytic cycle. Hydrogen-deuterium exchange coupled to mass spectrometry and molecular dynamics simulations revealed that the amino acid substitutions in BF-HsKYNase allosterically affect the flexibility of the pyridoxal-5′-phosphate (PLP) binding pocket, thereby impacting the rate of chemistry, presumably by altering the conformational ensemble and sampling states more favorable to the catalyzed reaction.

What have molecular simulations contributed to understanding of Gram-negative bacterial cell envelopes?

Bacterial cell envelopes are compositionally complex and crowded and while highly dynamic in some areas, their molecular motion is very limited, to the point of being almost static in others. Therefore, it is no real surprise that studying them at high resolution across a range of temporal and spatial scales requires a number of different techniques. Details at atomistic to molecular scales for up to tens of microseconds are now within range for molecular dynamics simulations. Here we review how such simulations have contributed to our current understanding of the cell envelopes of Gram-negative bacteria.

Atomistic basis of opening and conduction in mammalian inward rectifier potassium (Kir2.2) channels

This paper presents the crystal structure of a forced open inward rectifier Kir2.2 channel. Molecular dynamics reveals the details of ion permeation through the open channel.

Potassium ion conduction through open potassium channels is essential to control of membrane potentials in all cells. To elucidate the open conformation and hence the mechanism of K⁺ ion conduction in the classic inward rectifier Kir2.2, we introduced a negative charge (G178D) at the crossing point of the inner helix bundle, the location of ligand-dependent gating. This “forced open” mutation generated channels that were active even in the complete absence of phosphatidylinositol-4,5-bisphosphate (PIP₂), an otherwise essential ligand for Kir channel opening. Crystal structures were obtained at a resolution of 3.6 Å without PIP₂ bound, or 2.8 Å in complex with PIP₂. The latter revealed a slight widening at the helix bundle crossing (HBC) through backbone movement. MD simulations showed that subsequent spontaneous wetting of the pore through the HBC gate region allowed K⁺ ion movement across the HBC and conduction through the channel. Further simulations reveal atomistic details of the opening process and highlight the role of pore-lining acidic residues in K⁺ conduction through Kir2 channels.

Conformational changes upon gating of KirBac1.1 into an open-activated state revealed by solid-state NMR and functional assays

Inward rectifier K⁺ (Kir) channels play an important role in reestablishing the resting membrane state of the action potential of excitable cells in humans. KirBac1.1 is a prokaryotic Kir channel with a high degree of homology to human Kir channels and can be isotopically labeled in NMR quantities for structural studies. Functional assays and NMR assignments confirm that KirBac1.1 is in a constitutively conductive state. Solid-state NMR assignments further reveal alternate conformations at key sites in the protein that are well conserved through human Kir channels, hinting at a possible allosteric network between channels. These underlying sequential and structural motifs could explain abnormal conductive properties of these channels fundamental to their native gating processes.

The conformational changes required for activation and K⁺ conduction in inward-rectifier K⁺ (Kir) channels are still debated. These structural changes are brought about by lipid binding. It is unclear how this process relates to fast gating or if the intracellular and extracellular regions of the protein are coupled. Here, we examine the structural details of KirBac1.1 reconstituted into both POPC and an activating lipid mixture of 3:2 POPC:POPG (wt/wt). KirBac1.1 is a prokaryotic Kir channel that shares homology with human Kir channels. We establish that KirBac1.1 is in a constitutively active state in POPC:POPG bilayers through the use of real-time fluorescence quenching assays and Förster resonance energy transfer (FRET) distance measurements. Multidimensional solid-state NMR (SSNMR) spectroscopy experiments reveal two different conformers within the transmembrane regions of the protein in this activating lipid environment, which are distinct from the conformation of the channel in POPC bilayers. The differences between these three distinct channel states highlight conformational changes associated with an open activation gate and suggest a unique allosteric pathway that ties the selectivity filter to the activation gate through interactions between both transmembrane helices, the turret, selectivity filter loop, and the pore helix. We also identify specific residues involved in this conformational exchange that are highly conserved among human Kir channels.

High bandwidth approaches in nanopore and ion channel recordings - A tutorial review

Transport processes through ion-channel proteins, protein pores, or solid-state nanopores are traditionally recorded with commercial patch-clamp amplifiers. The bandwidth of these systems is typically limited to 10 kHz by signal-to-noise-ratio (SNR) considerations associated with these measurement platforms. At high bandwidth, the input-referred current noise in these systems dominates, determined by the input-referred voltage noise of the transimpedance amplifier applied across the capacitance at the input of the amplifier. This capacitance arises from several sources: the parasitic capacitance of the amplifier itself; the capacitance of the lipid bilayer harboring the ion channel protein (or the membrane used to form the solid-state nanopore); and the capacitance from the interconnections between the electronics and the membrane. Here, we review state-of-the-art applications of high-bandwidth conductance recordings of both ion channels and solid-state nanopores. These approaches involve tightly integrating measurement electronics fabricated in complementary metal-oxide semiconductors (CMOS) technology with lipid bilayer or solid-state membranes. SNR improvements associated with this tight integration push the limits of measurement bandwidths, in some cases in excess of 10 MHz. Recent case studies demonstrate the utility of these approaches for DNA sequencing and ion-channel recordings. In the latter case, studies with extended bandwidth have shown the potential for providing new insights into structure-function relations of these ion-channel proteins as the temporal resolutions of functional recordings matches time scales achievable with state-of-the-art molecular dynamics simulations.

Single-channel recordings of RyR1 at microsecond resolution in CMOS-suspended membranes

Single-channel recordings are widely used to explore functional properties of ion channels. Typically, such recordings are performed at bandwidths of less than 10 kHz because of signal-to-noise considerations, limiting the temporal resolution available for studying fast gating dynamics to greater than 100 µs. Here we present experimental methods that directly integrate suspended lipid bilayers with high-bandwidth, low-noise transimpedance amplifiers based on complementary metal-oxide-semiconductor (CMOS) integrated circuits (IC) technology to achieve bandwidths in excess of 500 kHz and microsecond temporal resolution. We use this CMOS-integrated bilayer system to study the type 1 ryanodine receptor (RyR1), a Ca²⁺-activated intracellular Ca²⁺-release channel located on the sarcoplasmic reticulum. We are able to distinguish multiple closed states not evident with lower bandwidth recordings, suggesting the presence of an additional Ca²⁺ binding site, distinct from the site responsible for activation. An extended beta distribution analysis of our high-bandwidth data can be used to infer closed state flicker events as fast as 35 ns. These events are in the range of single-file ion translocations.

Dynamics of the EAG1 K+ channel selectivity filter assessed by molecular dynamics simulations

EAG1 channels belong to the KCNH family of voltage gated potassium channels. They are expressed in several brain regions and increased expression is linked to certain cancer types. Recent cryo-EM structure determination finally revealed the structure of these channels in atomic detail, allowing computational investigations. In this study, we performed molecular dynamics simulations to investigate the ion binding sites and the dynamical behavior of the selectivity filter. Our simulations suggest that sites S2 and S4 form stable ion binding sites, while ions placed at sites S1 and S3 rapidly switched to sites S2 and S4. Further, ions tended to dissociate away from S0 within less than 20 ns, due to increased filter flexibility. This was followed by water influx from the extracellular side, leading to a widening of the filter in this region, and likely non-conductive filter configurations. Simulations with the inactivation-enhancing mutant Y464A or Na⁺ ions lead to trapped water molecules behind the SF, suggesting that these simulations captured early conformational changes linked to C-type inactivation.

Multi-scaled normal mode analysis method for dynamics simulation of protein-membrane complexes: A case study of potassium channel gating motion correlations

Membrane proteins play critically important roles in many cellular activities such as ions and small molecule transportation, signal recognition, and transduction. In order to fulfill their functions, these proteins must be placed in different membrane environments and a variety of protein-lipid interactions may affect the behavior of these proteins. One of the key effects of protein-lipid interactions is their ability to change the dynamics status of membrane proteins, thus adjusting their functions. Here, we present a multi-scaled normal mode analysis (mNMA) method to study the dynamics perturbation to the membrane proteins imposed by lipid bi-layer membrane fluctuations. In mNMA, channel proteins are simulated at all-atom level while the membrane is described with a coarse-grained model. mNMA calculations clearly show that channel gating motion can tightly couple with a variety of membrane deformations, including bending and twisting. We then examined bi-channel systems where two channels were separated with different distances. From mNMA calculations, we observed both positive and negative gating correlations between two neighboring channels, and the correlation has a maximum as the channel center-to-center distance is close to 2.5 times of their diameter. This distance is larger than recently found maximum attraction distance between two proteins embedded in membrane which is 1.5 times of the protein size, indicating that membrane fluctuation might impose collective motions among proteins within a larger area. The hybrid resolution feature in mNMA provides atomic dynamics information for key components in the system without costing much computer resource. We expect it to be a conventional simulation tool for ordinary laboratories to study the dynamics of very complicated biological assemblies. The source code is available upon request to the authors.

Identification of the Conformational transition pathway in PIP2 Opening Kir Channels

The gating of Kir channels depends critically on phosphatidylinositol 4,5-bisphosphate (PIP2), but the detailed mechanism by which PIP2 regulates Kir channels remains obscure. Here, we performed a series of Targeted molecular dynamics simulations on the full-length Kir2.1 channel and, for the first time, were able to achieve the transition from the closed to the open state. Our data show that with the upward motion of the cytoplasmic domain (CTD) the structure of the C-Linker changes from a loop to a helix. The twisting of the C-linker triggers the rotation of the CTD, which induces a small downward movement of the CTD and an upward motion of the slide helix toward the membrane that pulls the inner helix gate open. At the same time, the rotation of the CTD breaks the interaction between the CD- and G-loops thus releasing the G-loop. The G-loop then bounces away from the CD-loop, which leads to the opening of the G-loop gate and the full opening of the pore. We identified a series of interaction networks, between the N-terminus, CD loop, C linker and G loop one by one, which exquisitely regulates the global conformational changes during the opening of Kir channels by PIP2.

Molecular Dynamics Simulations of KirBac1.1 Mutants Reveal Global Gating Changes of Kir Channels

Prokaryotic inwardly rectifying (KirBac) potassium channels are homologous to mammalian Kir channels. Their activity is controlled by dynamical conformational changes that regulate ion flow through a central pore. Understanding the dynamical rearrangements of Kir channels during gating requires high-resolution structure information from channels crystallized in different conformations and insight into the transition steps, which are difficult to access experimentally. In this study, we use MD simulations on wild type KirBac1.1 and an activatory mutant to investigate activation gating of KirBac channels. Full atomistic MD simulations revealed that introducing glutamate in position 143 causes significant widening at the helix bundle crossing gate, enabling water flux into the cavity. Further, global rearrangements including a twisting motion as well as local rearrangements at the subunit interface in the cytoplasmic domain were observed. These structural rearrangements are similar to recently reported KirBac3.1 crystal structures in closed and open conformation, suggesting that our simulations capture major conformational changes during KirBac1.1 opening. In addition, an important role of protein-lipid interactions during gating was observed. Slide-helix and C-linker interactions with lipids were strengthened during activation gating.

Comparative study of the energetics of ion permeation in Kv1.2 and KcsA potassium channels

Biological ion channels rely on a multi-ion transport mechanism for fast yet selective permeation of ions. The crystal structure of the KcsA potassium channel provided the first microscopic picture of this process. A similar mechanism is assumed to operate in all potassium channels, but the validity of this assumption has not been well investigated. Here, we examine the energetics of ion permeation in Shaker Kv1.2 and KcsA channels, which exemplify the six-transmembrane voltage-gated and two-transmembrane inward-rectifier channels. We study the feasibility of binding a third ion to the filter and the concerted motion of ions in the channel by constructing the potential of mean force for K(+) ions in various configurations. For both channels, we find that a pair of K(+) ions can move almost freely within the filter, but a relatively large free-energy barrier hinders the K(+) ion from stepping outside the filter. We discuss the effect of the CMAP dihedral energy correction that was recently incorporated into the CHARMM force field on ion permeation dynamics.

Cooperative nature of gating transitions in K(+) channels as seen from dynamic importance sampling calculations

The growing dataset of K(+) channel x-ray structures provides an excellent opportunity to begin a detailed molecular understanding of voltage-dependent gating. These structures, while differing in sequence, represent either a stable open or closed state. However, an understanding of the molecular details of gating will require models for the transitions and experimentally testable predictions for the gating transition. To explore these ideas, we apply dynamic importance sampling to a set of homology models for the molecular conformations of K(+) channels for four different sets of sequences and eight different states. In our results, we highlight the importance of particular residues upstream from the Pro-Val-Pro (PVP) region to the gating transition. This supports growing evidence that the PVP region is important for influencing the flexibility of the S6 helix and thus the opening of the gating domain. The results further suggest how gating on the molecular level depends on intra-subunit motions to influence the cooperative behavior of all four subunits of the K(+) channel. We hypothesize that the gating process occurs in steps: first sidechain movement, then inter-S5-S6 subunit motions, and lastly the large-scale domain rearrangements.

Proton inhibition of unitary currents of vanilloid receptors

Protons, which are released during inflammation and injury, regulate many receptors and ion channels involved in pain transduction, including capsaicin channels (transient receptor potential vanilloid receptors 1). Whereas extracellular acidification both sensitizes and directly activates the channel, it also causes concomitant reduction of the unitary current amplitudes. Here, we investigate the mechanisms and molecular basis of this inhibitory effect of protons on channel conductance. Single-channel recordings showed that the unitary current amplitudes decreased with extracellular pH in a dose-dependent manner, consistent with a model in which protons bind to a site within the channel with an apparent pKa of ∼6. The inhibition was voltage dependent, ∼65% at −60 mV and 37% at +60 mV when pH was reduced from 7.4 to 5.5. The unitary current amplitudes reached saturation at [K⁺] ≥ 1 M, and notably the maximum amplitudes did not converge with different pHs, inconsistent with a blockade model based on surface charge screening or competitive inhibition of permeating ions. Mutagenesis experiments uncovered two acidic residues critical for proton inhibition, one located at the pore entrance and the other on the pore helix. Based on homology to the KcsA structure, the two acidic residues, along with another basic residue also on the pore helix, could form a triad interacting with each other through extensive hydrogen bonds and electrostatic contacts, suggesting that protons may mediate the interactions between the selectivity filter and pore helix, thereby altering the local structure in the filter region and consequently the conductance of the channel.

Dynamics, energetics, and selectivity of the low-K+ KcsA channel structure

Potassium channels are a diverse family of integral membrane proteins through which K(+) can pass selectively. There is ongoing debate about the nature of conformational changes associated with the opening/closing and conductive/nonconductive states of potassium channels. The channels partly exert their function by varying their conductance through a mechanism known as C-type inactivation. Shortly after the activation of K(+) channels, their selectivity filter stops conducting ions at a rate that depends on various stimuli. The molecular mechanism of C-type inactivation has not been fully understood yet. However, the X-ray structure of the KcsA channel obtained in the presence of low K(+) concentration is thought to be representative of a K(+) channel in the C-type inactivated state. Here, extensive, fully atomistic molecular dynamics and free-energy simulations of the low-K(+) KcsA structure in an explicit lipid bilayer are performed to evaluate the stability of this structure and the selectivity of its binding sites. We find that the low-K(+) KcsA structure is stable on the timescale of the molecular dynamics simulations performed, and that ions preferably remain in S1 and S4. In the absence of ions, the selectivity filter evolves toward an asymmetric architecture, as already observed in other computations of the high-K(+) structure of KcsA and KirBac. The low-K(+) KcsA structure is not permeable by Na(+), K(+), or Rb(+), and the selectivity of its binding sites is different from that of the high-K(+) structure.

Many-body effects and simulations of potassium channels

The electronic polarizability of an ion or a molecule is a measure of the relative tendency of its electron cloud to be distorted from its normal shape by an electric field. On the molecular scale, in a condensed phase, any species sits in an electric field due to its neighbours, and the resulting polarization is an important contribution to the total interaction energy. Electrostatic interactions are crucial for determining the majority of chemical–physical properties of the system and electronic polarization is a fundamental component of these interactions. Thus, polarization effects should be taken into account if accurate descriptions are desired. In classical computer simulations, the forces required to drive the system are typically based on interatomic interaction potentials derived in part from electronic structure calculations or from experimental data. Owing to the difficulties in including polarization effects in classical force fields, most of them are based just on pairwise additive interaction potentials. At present, major efforts are underway to develop polarizable interaction potentials for biomolecular simulations. In this review, various ways of introducing explicit polarizability into biomolecular models and force fields are reviewed, and the progress that might be achieved in applying such methods to study potassium channels is described.

Molecular modeling and simulation studies of ion channel structures, dynamics and mechanisms

Ion channels are integral membrane proteins that enable selected ions to flow passively across membranes. Channel proteins have been the focus of computational approaches to relate their three-dimensional (3D) structure to their physiological function. We describe a number of computational tools to model ion channels. Homology modeling may be used to construct structural models of channels based on available X-ray structures. Electrostatics calculations enable an approximate evaluation of the energy profile of an ion passing through a channel. Molecular dynamics simulations and free-energy calculations provide information on the thermodynamics and kinetics of channel function.

Model development for the viral Kcv potassium channel

A computational model for the open state of the short viral Kcv potassium channel was created and tested based on homology modeling and extensive molecular-dynamics simulation in a membrane environment. Particular attention was paid to the structure of the highly flexible N-terminal region and to the protonation state of membrane-exposed lysine residues. Data from various experimental sources, NMR spectroscopy, and electrophysiology, as well as results from three-dimensional reference interaction site model integral equation theory were taken into account to select the most reasonable model among possible variants. The final model exhibits spontaneous ion transitions across the complete pore, with and without application of an external field. The nonequilibrium transport events could be induced reproducibly without abnormally large driving potential and without the need to place ions artificially at certain key positions along the transition path. The transport mechanism through the filter region corresponds to the classic view of single-file motion, which in our case is coupled to frequent exchange of ions between the innermost filter position and the cavity.

Ion channels in microbes

Studies of ion channels have for long been dominated by the animalcentric, if not anthropocentric, view of physiology. The structures and activities of ion channels had, however, evolved long before the appearance of complex multicellular organisms on earth. The diversity of ion channels existing in cellular membranes of prokaryotes is a good example. Although at first it may appear as a paradox that most of what we know about the structure of eukaryotic ion channels is based on the structure of bacterial channels, this should not be surprising given the evolutionary relatedness of all living organisms and suitability of microbial cells for structural studies of biological macromolecules in a laboratory environment. Genome sequences of the human as well as various microbial, plant, and animal organisms unambiguously established the evolutionary links, whereas crystallographic studies of the structures of major types of ion channels published over the last decade clearly demonstrated the advantage of using microbes as experimental organisms. The purpose of this review is not only to provide an account of acquired knowledge on microbial ion channels but also to show that the study of microbes and their ion channels may also hold a key to solving unresolved molecular mysteries in the future.

Intrinsic Free Energy of the Conformational Transition of the KcsA Signature Peptide from Conducting to Nonconducting State

We explore a conformational transition of the TATTVGYG signature peptide of the KcsA ion selectivity filter and its GYG to AYA mutant from the conducting α-strand state into the nonconducting pII-like state using a novel technique for multidimensional optimization of transition path ensembles and free energy calculations. We find that the wild type peptide, unlike the mutant, intrinsically favors the conducting state due to G77 backbone propensities and additional hydrophobic interaction between the V76 and Y78 side chains in water. The molecular mechanical free energy profiles in explicit water are in very good agreement with the corresponding adiabatic energies from the Generalized Born Molecular Volume (GBMV) implicit solvent model. However comparisons of the energies to higher level B3LYP/6-31G(d) Density Functional Theory calculations with Polarizable Continuum Model (PCM) suggest that the nonconducting state might be more favorable than predicted by molecular mechanics simulations. By extrapolating the single peptide results to the tetrameric channel, we propose a novel hypothesis for the ion selectivity mechanism.

Conformational changes in the selectivity filter of the open-state KcsA channel: an energy minimization study

Potassium channels switch between closed and open conformations and selectively conduct K(+) ions. There are at least two gates. The TM2 bundle at the intracellular site is the primary gate of KcsA, and rearrangements at the selectivity filter (SF) act as the second gate. The SF blocks ion flow via an inactivation process similar to C-type inactivation of voltage-gated K(+) channels. We recently generated the open-state conformation of the KcsA channel. We found no major, possibly inactivating, structural changes in the SF associated with this massive inner-pore rearrangement, which suggests that the gates might act independently. Here we energy-minimize the open state of wild-type and mutant KcsA, validating in silico structures of energy-minimized SFs by comparison with crystallographic structures, and use these data to gain insight into how mutation, ion depletion, and K(+) to Na(+) substitution influence SF conformation. Both E71 or D80 protonations/mutations and the presence/absence of protein-buried water molecule(s) modify the H-bonding network stabilizing the P-loops, spawning numerous SF conformations. We find that the inactivated state corresponds to conformations with a partially unoccupied or an entirely empty SF. These structures, involving modifications in all four P-loops, are stabilized by H-bonds between amide H and carbonyl O atoms from adjacent P-loops, which block ion passage. The inner portions of the P-loops are more rigid than the outer parts. Changes are localized to the outer binding sites, with innermost site S4 persisting in the inactivated state. Strong binding by Na(+) locally contracts the SF around Na(+), releasing ligands that do not participate in Na(+) coordination, and occluding the permeation pathway. K(+) selectivity primarily appears to arise from the inability of the SF to completely dehydrate Na(+) ions due to basic structural differences between liquid water and the "quasi-liquid" SF matrix.

Insight into the mechanism of inactivation and pH sensitivity in potassium channels from molecular dynamics simulations

Potassium (K (+)) channels can regulate ionic conduction through their pore by a mechanism, involving the selectivity filter, known as C-type inactivation. This process is rapid in the hERG K (+) channel and is fundamental to its physiological role. Although mutations within hERG are known to remove this process, a structural basis for the inactivation mechanism has yet to be characterized. Using MD simulations based on homology modeling, we observe that the carbonyl of the filter aromatic, Phe627, forming the S 0 K (+) binding site, swiftly rotates away from the conduction axis in the wild-type channel. In contrast, in well-characterized non-inactivating mutant channels, this conformational change occurs less frequently. In the non-inactivating channels, interactions with a water molecule located behind the selectivity filter are critical to the enhanced stability of the conducting state. We observe comparable conformational changes in the acid sensitive TASK-1 channel and propose a common mechanism in these channels for regulating efflux of K (+) ions through the selectivity filter.

Hydropathic analysis and comparison of KcsA and Shaker potassium channels

The similarity in structure of potassium (K(+)) channels from different families has been revealed by only recently available crystallographic 3D structural data. The hydropathic analysis presented in this work illuminates whether homologous residues perform the same functions in channels that use different gating mechanisms. We calculated and compared the hydropathic profiles of two K(+) channels, KcsA and Kv1.2 (the latter a member of the Shaker family), at their pore-forming domain. Quantitative information describing important interactions stabilizing the protein beyond obvious secondary-structure elements was extracted from the analysis and applied as a template for subsequent molecular-dynamics (MD) analyses. For example, two key groups of interactions, defining the turns that connect the transmembrane helices and responsible for the orientation of the pore helix, were identified. Our results also indicate that Asp(80) and Asp(379) play a similar role in stabilizing the P-loop of KcsA and Kv1.2, respectively, but to significantly different extents.

Backbone dynamics of alamethicin bound to lipid membranes: spin-echo electron paramagnetic resonance of TOAC-spin labels

Alamethicin F50/5 is a hydrophobic peptide that is devoid of charged residues and that induces voltage-dependent ion channels in lipid membranes. The peptide backbone is likely to be involved in the ion conduction pathway. Electron spin-echo spectroscopy of alamethicin F50/5 analogs in which a selected Aib residue (at position n = 1, 8, or 16) is replaced by the TOAC amino-acid spin label was used to study torsional dynamics of the peptide backbone in association with phosphatidylcholine bilayer membranes. Rapid librational motions of limited angular amplitude were observed at each of the three TOAC sites by recording echo-detected spectra as a function of echo delay time, 2tau. Simulation of the time-resolved spectra, combined with conventional EPR measurements of the librational amplitude, shows that torsional fluctuations of the peptide backbone take place on the subnanosecond to nanosecond timescale, with little temperature dependence. Associated fluctuations in polar fields from the peptide could facilitate ion permeation.

Negatively charged residues located near the external entrance are required for the Kir2.1 channel to function

To understand the significance of negative charges in the extracellular loops of the inwardly rectifying K(+) (Kir2.1) channel, single-point mutants (D112N, D114N, E125Q, D152E, D152K, D152N, E153D, E153K, and E153Q) and double-point mutants (D112N/D114N and D152N/E153Q) were constructed and transfected into COS-1 and HEK293 cells. All single-point mutants, except D152K, and D112N/D114N expressed functional channels. Cells transfected with the D152K and D152N/E153Q constructs did not show any inwardly rectifying K(+) currents, although fluorescence images confirmed that the channel proteins produced by D152K and D152N/E153Q were transported to the cell surface. While a tandem tetramer with one D152N subunit and three D152N/E153Q subunits, D152N-(D152N/E153Q)(3), did not express functional channels, a tandem tetramer with one E153Q subunit and three D152N/E153Q subunits, E153Q-(D152N/E153Q)(3), and that with two D152N subunits and two D152N/E153Q subunits, (D152N)(2)-(D152N/E153Q)(2), expressed channels having similar conductance and kinetics of single-channel currents to the wild-type channels. These results suggest that one negative charge of D152 or two negative charges of E153 are required for Kir2.1 channels to function. It is suggested that the contribution by D152 and E153 to the electronegative extracellular pore entrance is critical for the channel to function properly.

The protonation state of the Glu-71/Asp-80 residues in the KcsA potassium channel: a first-principles QM/MM molecular dynamics study

Although a few x-ray structures of the KcsA K(+) channel have been crystallized several issues concerning the mechanisms of the ionic permeation and the protonation state of the selectivity filter ionizable side chains are still open. Using a first-principles quantum mechanical/molecular mechanical simulation approach, we have investigated the protonation state of Glu-71 and Asp-80, two important residues located in the vicinity of the selectivity filter. Results from the dynamics show that a proton is shared between the two residues, with a slight preference for Glu-71. The proton is found to exchange on the picosecond timescale, an interesting phenomenon that cannot be observed in classical molecular dynamics. Simulations of different ionic loading states of the filter show that the probability for the proton transfer is correlated with the filter occupancy. In addition, the Glu-71/Asp-80 pair is able to modulate the potential energy profile experienced by a K(+) ion as it translates along the pore axis. These theoretical predictions, along with recent experimental results, suggest that changes of the filter structure could be associated with a shift in the Glu-Asp protonation state, which in turn would influence the ion translocation.

Molecular Dynamics Simulation of the Cytosolic Mouth in Kcv-Type Potassium Channels

The functional effect of mutations near the intracellular mouth of the short viral Kcv potassium channel was studied by molecular dynamics simulations. As a model system we used the analogously mutated and truncated KirBac1.1, a channel with known crystal structure that shares genuine local sequence motifs with Kcv. By a novel simulated annealing methodology for structural averaging, information about the structure and dynamics of the intracellular mouth was extracted and complemented by Poisson-Boltzmann and 3D-RISM (reference interaction site model) integral equation theory for the determination of the K+ free energy surface. Besides the wild-type analogue of Kcv with its experimental reference activity (truncated KirBac1.1), two variants were studied: a deletion mutant where the N-terminus is further truncated by eight amino acids, showing inactivity in the Kcv reference system, and a point mutant where the kink-forming proline at position 13 is substituted by alanine, resulting in hyperactivity. The computations reveal that the change of activity is closely related to a hydrophilic intracellular constriction formed by the C-terminal residues of the monomers. Hyperactivity of the point mutant is correlated with both sterical and electrostatic factors, while inactivity of the deletion mutant is related to a loss of specific salt bridge patterns between the C- and N-terminus at the constriction and to the consequences for ion passage barriers, as revealed by integral equation theory. The cytosolic gate, however, is probably formed by the N-terminal segment up to the proline kink and not by the constriction. The results are compared with design principles found for other channels.

Molecular dynamics simulations of inwardly rectifying (Kir) potassium channels: a comparative study

Inward rectifier potassium (Kir) channels regulate cell excitability and transport K+ ions across membranes. Homotetrameric models of three mammalian Kir channels (Kir1.1, Kir3.1, and Kir6.2) have been generated, using the KirBac3.1 transmembrane and rat Kir3.1 intracellular domain structures as templates. All three models have been explored by 10 ns molecular dynamics simulations in phospholipid bilayers. Analysis of the initial structures revealed conservation of potential lipid interaction residues (Trp/Tyr and Arg/Lys side chains near the lipid headgroup-water interfaces). Examination of the intracellular domains revealed key structural differences between Kir1.1 and Kir6.2 which may explain the difference in channel inhibition by ATP. The behavior of all three models in the MD simulations revealed that they have conformational stability similar to that seen for comparable simulations of, for example, structures derived from cryoelectron microscopy data. Local distortions of the selectivity filter were seen during the simulations, as observed in previous simulations of KirBac and in simulations and structures of KcsA. These may be related to filter gating of the channel. The intracellular hydrophobic gate does not undergo any substantial changes during the simulations and thus remains functionally closed. Analysis of lipid-protein interactions of the Kir models emphasizes the key role of the M0 (or "slide") helix which lies approximately parallel to the bilayer-water interface and forms a link between the transmembrane and intracellular domains of the channel.

Dynamics of K+ ion conduction through Kv1.2

The crystallographic structure of a potassium channel, Kv1.2, in an open state makes it feasible to simulate entire K(+) ion permeation events driven by a voltage bias and, thereby, elucidate the mechanism underlying ion conduction and selectivity of this type of channel. This Letter demonstrates that molecular dynamics simulations can provide movies of the overall conduction of K(+) ions through Kv1.2. As suggested earlier, the conduction is concerted in the selectivity filter, involving 2-3 ions residing mainly at sites identified previously by crystallography and modeling. The simulations reveal, however, the jumps of ions between these sites and identify the sequence of multi-ion configurations involved in permeation.

Unified modeling of conductance kinetics for low- and high-conductance potassium ion channels

A kinetics model is proposed for the description of ion conductance of low- and high-conductance potassium ion channels. The model describes ion permeation through the selectivity filter, which is assumed to be the only conductance determining part of the open channel. The filter occupancy can vary from zero to three ions, affecting the ion entry and exit rates. Ion motion between the binding sites inside the filter is assumed fast compared to the latter rates allowing averaging the equilibrium entry and exit rate constants over the possible ion configurations in the filter with a particular occupancy. Averaged rate constants related to a pair of adjacent occupancy states characterize a particular ion permeation mechanism. An expression for the channel conductance as a function of the symmetrical external ion concentration is derived. It comprises a sum of concentration independent conductance amplitudes for different ion permeation mechanisms weighted by the equilibrium filter occupancy probabilities. It is shown that each amplitude (i.e., maximum contribution to the channel conductance from each conductance mechanism) is proportional to an averaged exit rate constant and to quantities characterizing the effect of the applied electric field on the rate constants and the equilibrium ion distribution in the filter. The conductance expression derived provides a good description of the experimentally observed conductance-concentration curves for low-conductance (e.g., Kir2.1) and high-conductance (e.g., KcsA) potassium channels. It enables one to obtain equilibrium ion binding constants at different filter occupancies and to calculate the average number of ions in the selectivity filter for a given external ion concentration. For KcsA this number (2.0 at 200 mM) is in a good agreement with the available experimental value (2.1 at 200 mM). For the high-conductance potassium channels the net negative electrical charge around the selectivity filter increases the ion binding constants, thereby causing the larger occupancy probabilities to occur at smaller external ion concentrations compared to the low-conductance channels. This substantially increases the contributions of the two- and three-ion permeation mechanisms, with the larger conductance amplitudes leading to increased channel conductance compared to the low-conductance channels.

Base of pore loop is important for rectification, activation, permeation, and block of Kir3.1/Kir3.4

The Kir3.1/Kir3.4 channel is an inward rectifier, agonist-activated K(+) channel. The location of the binding site within the channel pore that coordinates polyamines (and is thus responsible for inward rectification) and the location of the gate that opens the channel in response to agonist activation is unclear. In this study, we show, not surprisingly, that mutation of residues at the base of the selectivity filter in the pore loop and second transmembrane domain weakens Cs(+) block and decreases selectivity (as measured by Rb(+) and spermine permeation). However, unexpectedly, the mutations also weaken inward rectification and abolish agonist activation of the channel. In the wild-type channel and 34 mutant channels, there are significant (p < 0.05) correlations among the K(D) for Cs(+) block, Rb(+) and spermine permeation, inward rectification, and agonist activation. The significance of these findings is discussed. One possible conclusion is that the selectivity filter is responsible for inward rectification and agonist activation as well as permeation and block.

Elongation of Outer Transmembrane Domain Alters Function of Miniature K+ Channel Kcv

The virus-coded channel Kcv has the typical structure of a two-transmembrane domain K(+) channel. Exceptional are its cytoplasmic domains: the C terminus basically ends inside the membrane and, hence, precludes the formation of a cytoplasmic gate by the so-called bundle crossing; the cytoplasmic N terminus is composed of only 12 amino acids. According to structural predictions, it is positioned in the membrane/aqueous interface and connected via a proline kink to the outer transmembrane domain (TM1). Here, we show that this proline kink affects channel function by determining the position of TM1 in the membrane bilayer. Extension of the hydrophobic length of TM1 by either eliminating the proline kink or introducing an alanine in TM1 augments a time- and voltage-dependent inward rectification of the channel. This suggests that the positional information of TM1 in the bilayer is transmitted to a channel gate, which is not identical with the cytoplasmic bundle crossing.

Conformational dynamics of M2 helices in KirBac channels: helix flexibility in relation to gating via molecular dynamics simulations

KirBac1.1 and 3.1 are bacterial homologues of mammalian inward rectifier K channels. We have performed extended molecular dynamics simulations (five simulations, each of >20 ns duration) of the transmembrane domain of KirBac in two membrane environments, a palmitoyl oleoyl phosphatidylcholine bilayer and an octane slab. Analysis of these simulations has focused on the conformational dynamics of the pore-lining M2 helices, which form the cytoplasmic hydrophobic gate of the channel. Principal components analysis reveals bending of M2, with a molecular hinge at the conserved glycine (Gly134 in KirBac1.1, Gly120 in KirBac3.1). More detailed analysis reveals a dimer-of-dimers type motion. The first two eigenvectors describing the motions of M2 correspond to helix kink and swivel motions. The conformational flexibility of M2 seen in these simulations correlates with differences in M2 conformation between that seen in the X-ray structures of closed channels (KcsA and KirBac) in which the helix is undistorted, and in open channels (e.g. MthK) in which the M2 helix is kinked. Thus, the simulations, albeit on a time scale substantially shorter than that required for channel gating, suggest a gating model in which the intrinsic flexibility of M2 about a molecular hinge is coupled to conformational transitions of an intracellular 'gatekeeper' domain, the latter changing conformation in response to ligand binding.

Molecular dynamics simulations of proteins in lipid bilayers

With recent advances in X-ray crystallography of membrane proteins promising many new high-resolution structures, molecular dynamics simulations will become increasingly valuable for understanding membrane protein function, as they can reveal the dynamic behavior concealed in the static structures. Dramatic increases in computational power, in synergy with more efficient computational methodologies, now allow us to carry out molecular dynamics simulations of any structurally known membrane protein in its native environment, covering timescales of up to 0.1 micros. At the frontiers of membrane protein simulations are ion channels, aquaporins, passive and active transporters, and bioenergetic proteins.

Anionic phospholipid interactions with the potassium channel KcsA: simulation studies

Molecular dynamics (MD) simulations have been used to unmask details of specific interactions of anionic phospholipids with intersubunit binding sites on the surface of the bacterial potassium channel KcsA. Crystallographic data on a diacyl glycerol fragment at this site were used to model phosphatidylethanolamine (PE), or phosphatidylglycerol (PG), or phosphatidic acid (PA) at the intersubunit binding sites. Each of these models of a KcsA-lipid complex was embedded in phosphatidyl choline bilayer and explored in a 20 ns MD simulation. H-bond analysis revealed that in terms of lipid-protein interactions PA > PG >> PE and revealed how anionic lipids (PG and PA) bind to a site provided by two key arginine residues (R(64) and R(89)) at the interface between adjacent subunits. A 27 ns simulation was performed in which KcsA (without any lipids initially modeled at the R(64)/R(89) sites) was embedded in a PE/PG bilayer. There was a progressive specific increase over the course of the simulation in the number of H-bonds of PG with KcsA. Furthermore, two specific PG binding events at R(64)/R(89) sites were observed. The phosphate oxygen atoms of bound PG formed H-bonds to the guanidinium group of R(89), whereas the terminal glycerol H-bonded to R(64). Overall, this study suggests that simulations can help identify and characterize sites for specific lipid interactions on a membrane protein surface.

A comparison between two prokaryotic potassium channels (KirBac1.1 and KcsA) in a molecular dynamics (MD) simulation study

The two potassium ion channels KirBac1.1 and KcsA are compared in a Molecular Dynamics (MD) simulation study. The location and motion of the potassium ions observed in the simulations are compared to those in the X-ray structures and previous simulations. In our simulations several of the crystallography resolved ion sites in KirBac1.1 are occupied by ions. In addition to this, two in KirBac1.1 unresolved sites where occupied by ions at sites that are in close correspondence to sites found in KcsA. There is every reason to believe that the conserved alignment of the selectivity filter in the potassium ion channel family corresponds to a very similar mechanism for ion transport across the filter. The gate residues, Phe146 in KirBac1.1 and Ala111 in KcsA acted in the simulations as effective barriers which never were passed by ions nor water molecules.

Ion conduction and selectivity in K(+) channels

Potassium (K(+)) channels are tetrameric membrane-spanning proteins that provide a selective pore for the conductance of K(+) across the cell membranes. These channels are most remarkable in their ability to discriminate K(+) from Na(+) by more than a thousandfold and conduct at a throughput rate near diffusion limit. The recent progress in the structural characterization of K(+) channel provides us with a unique opportunity to understand their function at the atomic level. With their ability to go beyond static structures, molecular dynamics simulations based on atomic models can play an important role in shaping our view of how ion channels carry out their function. The purpose of this review is to summarize the most important findings from experiments and computations and to highlight a number of fundamental mechanistic questions about ion conduction and selectivity that will require further work.

Molecular dynamics simulation approaches to K channels: conformational flexibility and physiological function

Molecular modeling and simulations enable extrapolation for the structure of bacterial potassium channels to the function of their mammalian homologues. Molecular dynamics simulations have revealed the concerted single-file motion of potassium ions and water molecules through the selectivity filter of K channels and the role of filter flexibility in ion permeation and in "fast gating." Principal components analysis of extended K channel simulations suggests that hinge-bending of pore-lining M2 (or S6) helices plays a key role in K channel gating. Based on these and other simulations, a molecular model for gating of inward rectifier K channel gating is presented.

A gate in the selectivity filter of potassium channels

The selectivity filter of potassium channels is the structural element directly responsible for the selective and rapid conduction of K+, whereas other parts of the protein are thought to function as a molecular gate that either permits or blocks the passage of ions. However, whether the selectivity filter itself also possesses the ability to play the role of a gate is an unresolved question. Using free energy molecular dynamics simulations, it is shown that the reorientation of two peptide linkages in the selectivity filter of the KcsA K+ channel can lead to a stable nonconducting conformational state. Two microscopic factors influence the transition toward such a conformational state: the occupancy of one specific cation binding site in the selectivity filter (S2), and the strength of intersubunit interactions involving the GYG signature sequence. These results suggest that such conformational transitions occurring in the selectivity filter might be related to different K+ channel gating events, including C-type (slow) inactivation.

Conformational dynamics of the ligand-binding domain of inward rectifier K channels as revealed by molecular dynamics simulations: toward an understanding of Kir channel gating

Inward rectifier (Kir) potassium channels are characterized by two transmembrane helices per subunit, plus an intracellular C-terminal domain that controls channel gating in response to changes in concentration of various ligands. Based on the crystal structure of the tetrameric C-terminal domain of Kir3.1, it is possible to build a homology model of the ATP-binding C-terminal domain of Kir6.2. Molecular dynamics simulations have been used to probe the dynamics of Kir C-terminal domains and to explore the relationship between their dynamics and possible mechanisms of channel gating. Multiple simulations, each of 10 ns duration, have been performed for Kir3.1 (crystal structure) and Kir6.2 (homology model), in both their monomeric and tetrameric forms. The Kir6.2 simulations were performed with and without bound ATP. The results of the simulations reveal comparable conformational stability for the crystal structure and the homology model. There is some decrease in conformational flexibility when comparing the monomers with the tetramers, corresponding mainly to the subunit interfaces in the tetramer. The beta-phosphate of ATP interacts with the side chain of K185 in the Kir6.2 model and simulations. The flexibility of the Kir6.2 tetramer is not changed greatly by the presence of bound ATP, other than in two loop regions. Principal components analysis of the simulated dynamics suggests loss of symmetry in both the Kir3.1 and Kir6.2 tetramers, consistent with "dimer-of-dimers" motion of subunits in C-terminal domains of the corresponding Kir channels. This is suggestive of a gating model in which a transition between exact tetrameric symmetry and dimer-of-dimers symmetry is associated with a change in transmembrane helix packing coupled to gating of the channel. Dimer-of-dimers motion of the C-terminal domain tetramer is also supported by coarse-grained (anisotropic network model) calculations. It is of interest that loss of exact rotational symmetry has also been suggested to play a role in gating in the bacterial Kir homolog, KirBac1.1, and in the nicotinic acetylcholine receptor channel.