Probing the energy landscape of activation gating of the bacterial potassium channel KcsA

Transition Path Dynamics of Non-Markovian Systems across a Rough Potential Barrier

Transition paths refer to rare events in physics, chemistry, and biology where the molecules cross barriers separating stable molecular conformations. The conventional analysis of the transition path times employs a diffusive and memoryless transition over a smooth potential barrier. However, it is widely acknowledged that the free energy profile between two minima in biomolecular processes is inherently not smooth. In this article, we discuss a theoretical model with a parabolic rough potential barrier and obtain analytical results of the transition path distribution and mean transition path times by incorporating absorbing boundary conditions across the boundaries under the driving of Gaussian white noise. Further, the influence of anomalous dynamics in rough potential driven by a power-law memory kernel is analyzed by deriving a time-dependent scaled diffusion coefficient that coarse-grains the effects of roughness, and the system's dynamics is reduced to a scaled diffusion on a smooth potential. Our theoretical results are tested and validated against numerical simulations. The findings of our study show the influence of the boundary conditions, barrier height, barrier roughness, and memory effect on the transition path time distributions in a rough potential, and the validity of the scaling diffusion coefficient has been discussed.

Mechanism of external K+ sensitivity of KCNQ1 channels

KCNQ1 voltage-gated K+ channels are involved in a wide variety of fundamental physiological processes and exhibit the unique feature of being markedly inhibited by external K+. Despite the potential role of this regulatory mechanism in distinct physiological and pathological processes, its exact underpinnings are not well understood. In this study, using extensive mutagenesis, molecular dynamics simulations, and single-channel recordings, we delineate the molecular mechanism of KCNQ1 modulation by external K+. First, we demonstrate the involvement of the selectivity filter in the external K+ sensitivity of the channel. Then, we show that external K+ binds to the vacant outermost ion coordination site of the selectivity filter inducing a diminution in the unitary conductance of the channel. The larger reduction in the unitary conductance compared to whole-cell currents suggests an additional modulatory effect of external K+ on the channel. Further, we show that the external K+ sensitivity of the heteromeric KCNQ1/KCNE complexes depends on the type of associated KCNE subunits.

Roughness in the periodic potential induces absolute negative mobility in a driven Brownian ratchet

Absolute negative mobility, where particles move opposite to the direction as governed by the external load, is an anomalous transport property of a Brownian ratchet and has technological implications in mass separation and bioanalytical applications. We numerically investigated here the effect of roughness in symmetric periodic potential on the negative mobility of a driven inertial Brownian ratchet in the presence of an external load. We show that the microscopic spatial heterogeneity of the potential can generate negative mobility which would not otherwise be possible under smooth potential in the concerned parameter space. We determined the optimal condition in terms of parameter space for such anomalous behavior. Our calculations indicate that the shift of balance towards the negative velocity phase in the temporal oscillations of velocity and weakly chaotic dynamics are responsible factors for roughness-induced negative mobility. These calculations highlight a constructive role of roughness in the anomalous transport properties of Brownian ratchet.

Intermolecular functional coupling between phosphoinositides and the potassium channel KcsA

Biological membranes are composed of a wide variety of lipids. Phosphoinositides (PIPns) in the membrane inner leaflet only account for a small percentage of the total membrane lipids but modulate the functions of various membrane proteins, including ion channels, which play important roles in cell signaling. KcsA, a prototypical K⁺ channel that is small, simple, and easy to handle, has been broadly examined regarding its crystallography, in silico molecular analysis, and electrophysiology. It has been reported that KcsA activity is regulated by membrane phospholipids, such as phosphatidylglycerol. However, there has been no quantitative analysis of the correlation between direct lipid binding and the functional modification of KcsA, and it is unknown whether PIPns modulate KcsA function. Here, using contact bubble bilayer recording, we observed that the open probability of KcsA increased significantly (from about 10% to 90%) when the membrane inner leaflet contained only a small percentage of PIPns. In addition, we found an increase in the electrophysiological activity of KcsA correlated with a larger number of negative charges on PIPns. We further analyzed the affinity of the direct interaction between PIPns and KcsA using microscale thermophoresis and observed a strong correlation between direct lipid binding and the functional modification of KcsA. In conclusion, our approach was able to reconstruct the direct modification of KcsA by PIPns, and we propose that it can also be applied to elucidate the mechanism of modification of other ion channels by PIPns.

Roughness in the periodic potential enhances transport in a driven inertial ratchet

We study the effects of roughness in the asymmetric periodic potential on the transport and diffusion of an inertial Brownian particle driven by a time-periodic force in a Gaussian environment. We find that moderate roughness leads to the loss of transient anomalous diffusion, and it helps to establish normal diffusion in the weak noise limit. We uncover a contrasting effect of roughness on the transport of particles in the weak and moderate to large noise limit. In the weak noise limit, small amplitude roughness results in the increase of directed transport, whereas in the moderate to large noise limit, roughness hinders transport. The deterministic dynamics of the system reveals that the purely periodic system under smooth potential transits into a chaotic system due to the moderate roughness in the potential. Therefore our calculations demonstrate the constructive role of roughness in the transport of particles in the inertial regime.

The Role of ICL1 and H8 in Class B1 GPCRs; Implications for Receptor Activation

The first intracellular loop (ICL1) of G protein-coupled receptors (GPCRs) has received little attention, although there is evidence that, with the 8^th helix (H8), it is involved in early conformational changes following receptor activation as well as contacting the G protein β subunit. In class B1 GPCRs, the distal part of ICL1 contains a conserved R^12.48KLRCxR^2.46b motif that extends into the base of the second transmembrane helix; this is weakly conserved as a [R/H]^12.48KL[R/H] motif in class A GPCRs. In the current study, the role of ICL1 and H8 in signaling through cAMP, _iCa²⁺ and ERK1/2 has been examined in two class B1 GPCRs, using mutagenesis and molecular dynamics. Mutations throughout ICL1 can either enhance or disrupt cAMP production by CGRP at the CGRP receptor. Alanine mutagenesis identified subtle differences with regard elevation of _iCa²⁺, with the distal end of the loop being particularly sensitive. ERK1/2 activation displayed little sensitivity to ICL1 mutation. A broadly similar pattern was observed with the glucagon receptor, although there were differences in significance of individual residues. Extending the study revealed that at the CRF1 receptor, an insertion in ICL1 switched signaling bias between _iCa²⁺ and cAMP. Molecular dynamics suggested that changes in ICL1 altered the conformation of ICL2 and the H8/TM7 junction (ICL4). For H8, alanine mutagenesis showed the importance of E390^8.49b for all three signal transduction pathways, for the CGRP receptor, but mutations of other residues largely just altered ERK1/2 activation. Thus, ICL1 may modulate GPCR bias via interactions with ICL2, ICL4 and the Gβ subunit.

Computing Potential of the Mean Force Profiles for Ion Permeation Through Channelrhodopsin Chimera, C1C2

Umbrella sampling, coupled with a weighted histogram analysis method (US-WHAM), can be used to construct potentials of mean force (PMFs) for studying the complex ion permeation pathways of membrane transport proteins. Despite the widespread use of US-WHAM, obtaining a physically meaningful PMF can be challenging. Here, we provide a protocol to resolve that issue. Then, we apply that protocol to compute a meaningful PMF for sodium ion permeation through channelrhodopsin chimera, C1C2, for illustration.

Probing allosteric coupling in a constitutively open mutant of the ion channel KcsA using solid-state NMR

Transmembrane allosteric coupling is a feature of many critical biological signaling events. Here we test whether transmembrane allosteric coupling controls the potassium binding affinity of the prototypical potassium channel KcsA in the context of C-type inactivation. Activation of KcsA is initiated by proton binding to the pH gate upon an intracellular drop in pH. Numerous studies have suggested that this proton binding also prompts a conformational switch, leading to a loss of affinity for potassium ions at the selectivity filter and therefore to channel inactivation. We tested this mechanism for inactivation using a KcsA mutant (H25R/E118A) that exhibits an open pH gate across a broad range of pH values. We present solid-state NMR measurements of this open mutant at neutral pH to probe the affinity for potassium at the selectivity filter. The potassium binding affinity in the selectivity filter of this mutant, 81 mM, is about four orders of magnitude weaker than that of wild-type KcsA at neutral pH and is comparable to the value for wild-type KcsA at low pH (pH ≈ 3.5). This result strongly supports our assertion that the open pH gate allosterically affects the potassium binding affinity of the selectivity filter. In this mutant, the protonation state of a glutamate residue (E120) in the pH sensor is sensitive to potassium binding, suggesting that this mutant also has flexibility in the activation gate and is subject to transmembrane allostery.

Resonance-driven ion transport and selectivity in prokaryotic ion channels

Ion channels exhibit a remarkably high accuracy in selecting uniquely its associated type of ion. The mechanisms behind ion selectivity are not well understood. Current explanations build mainly on molecular biology and bioinformatics. Here we propose a simple physical model for ion selectivity based on the driven damped harmonic oscillator (DDHO). The driving force for this oscillator is provided by self-organizing harmonic turbulent structures in the dehydrating ionic flow through the ion channel, namely, oscillating pressure waves in one dimension, and toroidal vortices in two and three dimensions. Density fluctuations caused by these turbulences efficiently transmit their energy to aqua ions that resonate with the driving frequency. Consequently, these release their hydration shell and exit the ion channel as free ions. Existing modeling frameworks do not express the required complex spatiotemporal dynamics, hence we introduce a macroscopic continuum model for ionic dehydration and transport, based on the hydrodynamics of a dissipative ionic flow through an ion channel, subject to electrostatic and amphiphilic interactions. This model combines three classical physical fields: Navier-Stokes equations from hydrodynamics, Gauss's law from Maxwell theory, and the convection-diffusion equation from continuum physics. Numerical experiments with mixtures of chemical species of ions in various degrees of hydration indeed reveal the emergence of strong oscillations in the ionic flow that are instrumental in the efficient dehydration and cause a strong ionic jet into the cell. As such, they provide a powerful engine for the DDHO mechanism. Theoretical predictions of the modeling framework match significantly with empirical patch-clamp data. The DDHO standard response curve defines a unique resonance frequency that depends on the mass and charge of the ion. In this way, the driving oscillations act as a selection mechanism that filters out one specific ion. Application of the DDHO model to real ion data shows that this mechanism indeed clearly distinguishes between chemical species and between aqua and bare ions with a large Mahalanobis distance and high oscillator quality. The DDHO framework helps to understand how SNP mutations can cause severe genetic pathologies as they destroy the geometry of the channel and so alter the resonance peaks of the required ion type.

Conduction through a narrow inward-rectifier K+ channel pore

G-protein–gated inwardly rectifying potassium channels are important mediators of inhibitory neurotransmission. Based on microsecond-scale molecular dynamics simulations, Bernsteiner et al. propose novel gating details that may enable K⁺ flux via a direct knock-on mechanism.

Inwardly rectifying potassium (Kir) channels play a key role in controlling membrane potentials in excitable and unexcitable cells, thereby regulating a plethora of physiological processes. G-protein–gated Kir channels control heart rate and neuronal excitability via small hyperpolarizing outward K⁺ currents near the resting membrane potential. Despite recent breakthroughs in x-ray crystallography and cryo-EM, the gating and conduction mechanisms of these channels are poorly understood. MD simulations have provided unprecedented details concerning the gating and conduction mechanisms of voltage-gated K⁺ and Na⁺ channels. Here, we use multi-microsecond–timescale MD simulations based on the crystal structures of GIRK2 (Kir3.2) bound to phosphatidylinositol-4,5-bisphosphate to provide detailed insights into the channel’s gating dynamics, including insights into the behavior of the G-loop gate. The simulations also elucidate the elementary steps that underlie the movement of K⁺ ions through an inward-rectifier K⁺ channel under an applied electric field. Our simulations suggest that K⁺ permeation might occur via direct knock-on, similar to the mechanism recently shown for K_v channels.

Inactivation in the potassium channel KcsA

•

C-type inactivation in potassium channels is a nearly universal regulatory mechanism.

•

A major hypothesis states that C-type inactivation involves ion loss at the selectivity filter as an allosteric response to activation.

•

NMR is used to probe protein conformational changes in response to pH and [K⁺], demonstrating that H⁺ and K⁺ binding are allosterically coupled in KcsA.

•

The lipids are integrated parts of potassium channels in terms of structure, energetics and function.

Inactivation, the slow cessation of transmission after activation, is a general feature of potassium channels. It is essential for their function, and malfunctions in inactivation leads to numerous pathologies. The detailed mechanism for the C-type inactivation, distinct from the N-type inactivation, remains an active area of investigation. Crystallography, computational simulations, and NMR have greatly enriched our understanding of the process. Here we review the major hypotheses regarding C-type inactivation, particularly focusing on the key role played by NMR studies of the prokaryotic potassium channel KcsA, which serves as a good model for voltage gated mammalian channels.

First-passage-time distribution in a moving parabolic potential with spatial roughness

In this paper, we investigate the first-passage-time distribution (FPTD) within a time-dependent parabolic potential in the presence of roughness with two methods: the Kramers theory and a nonsingular integral equation. By spatially averaging, the rough potential is equivalent to the combination of an effective smooth potential and an effective diffusion coefficient. Based on the Kramers theory, we first obtain Kramers approximations (KAs) of FPTD for both smooth and rough potentials. As expected, KA is valid only for high barriers and small external forces, and generally applicable for high barriers in rough potentials. To overcome the shortcoming of KA, a probability asymptotic approximation (PAA) based on an integral equation is proposed, which uses the transient probability density function (PDF) of the natural boundary conditions instead of the absorbing boundary conditions. We find that PAA fits very well even for large external forces. This enables us to analytically solve the FPTD for large external forces and low barriers as a strong extension to KA. In addition, we show that in the presence of a rough potential, the PAA of FPTD is in good agreement with numerical simulations for low barrier potentials. The PAA makes it possible to investigate the first-passage problem with ultrafast varying potentials and short exiting time. Thus, KA and PAA are complementary in determining the FPTD both for various barriers and external forces. Finally, the mean first-passage time (MFPT) is studied, which illustrates that the PAA of MFPT is effective almost in the whole range of external forces, while the KA of MFPT is valid only for small external forces.

Rapid constriction of the selectivity filter underlies C-type inactivation in the KcsA potassium channel

C-type inactivation in K⁺ channels is thought to be a result of constriction of the selectivity filter. By using MD simulations, Li et al. show that rapid constriction occurs within 1–2 s when the intracellular activation gate is fully open, but not when the gate is closed or partially open.

C-type inactivation is a time-dependent process observed in many K⁺ channels whereby prolonged activation by an external stimulus leads to a reduction in ionic conduction. While C-type inactivation is thought to be a result of a constriction of the selectivity filter, the local dynamics of the process remain elusive. Here, we use molecular dynamics (MD) simulations of the KcsA channel to elucidate the nature of kinetically delayed activation/inactivation gating coupling. Microsecond-scale MD simulations based on the truncated form of the KcsA channel (C-terminal domain deleted) provide a first glimpse of the onset of C-type inactivation. We observe over multiple trajectories that the selectivity filter consistently undergoes a spontaneous and rapid (within 1–2 µs) transition to a constricted conformation when the intracellular activation gate is fully open, but remains in the conductive conformation when the activation gate is closed or partially open. Multidimensional umbrella sampling potential of mean force calculations and nonequilibrium voltage-driven simulations further confirm these observations. Electrophysiological measurements show that the truncated form of the KcsA channel inactivates faster and greater than full-length KcsA, which is consistent with truncated KcsA opening to a greater degree because of the absence of the C-terminal domain restraint. Together, these results imply that the observed kinetics underlying activation/inactivation gating reflect a rapid conductive-to-constricted transition of the selectivity filter that is allosterically controlled by the slow opening of the intracellular gate.

Effects of protein-protein interactions and ligand binding on the ion permeation in KCNQ1 potassium channel

The voltage-gated KCNQ1 potassium ion channel interacts with the type I transmembrane protein minK (KCNE1) to generate the slow delayed rectifier (I_Ks) current in the heart. Mutations in these transmembrane proteins have been linked with several heart-related issues, including long QT syndromes (LQTS), congenital atrial fibrillation, and short QT syndrome. Off-target interactions of several drugs with that of KCNQ1/KCNE1 ion channel complex have been known to cause fatal cardiac irregularities. Thus, KCNQ1/KCNE1 remains an important avenue for drug-design and discovery research. In this work, we present the structural and mechanistic details of potassium ion permeation through an open KCNQ1 structural model using the combined molecular dynamics and steered molecular dynamics simulations. We discuss the processes and key residues involved in the permeation of a potassium ion through the KCNQ1 ion channel, and how the ion permeation is affected by (i) the KCNQ1-KCNE1 interactions and (ii) the binding of chromanol 293B ligand and its derivatives into the complex. The results reveal that interactions between KCNQ1 with KCNE1 causes a pore constriction in the former, which in-turn forms small energetic barriers in the ion-permeation pathway. These findings correlate with the previous experimental reports that interactions of KCNE1 dramatically slows the activation of KCNQ1. Upon ligand-binding onto the complex, the energy-barriers along ion permeation path are more pronounced, as expected, therefore, requiring higher force in our steered-MD simulations. Nevertheless, pulling the ion when a weak blocker is bound to the channel does not necessitate high force in SMD. This indicates that our SMD simulations have been able to discern between strong and week blockers and reveal their influence on potassium ion permeation. The findings presented here will have some implications in understanding the potential off-target interactions of the drugs with the KCNQ1/KCNE1 channel that lead to cardiotoxic effects.

Conformational landscapes of membrane proteins delineated by enhanced sampling molecular dynamics simulations

The expansion of computational power, better parameterization of force fields, and the development of novel algorithms to enhance the sampling of the free energy landscapes of proteins have allowed molecular dynamics (MD) simulations to become an indispensable tool to understand the function of biomolecules. The temporal and spatial resolution of MD simulations allows for the study of a vast number of processes of interest. Here, we review the computational efforts to uncover the conformational free energy landscapes of a subset of membrane proteins: ion channels, transporters and G-protein coupled receptors. We focus on the various enhanced sampling techniques used to study these questions, how the conclusions come together to build a coherent picture, and the relationship between simulation outcomes and experimental observables.

Roughness-enhanced transport in a tilted ratchet driven by Lévy noise

The enhanced transport of particles by roughness in a tilted rough ratchet potential subject to a Lévy noise is investigated in this paper. Due to the roughness, the transport process exhibits quite different properties compared to the smooth case. We find that the roughness on the potential wall functions like a ladder to provide the convenience for particles to climb up but hinder them to slide down. The mean first passage time from one well to its right adjacent well and the mean velocity are, respectively, calculated versus the roughness, the external force, and the Lévy stability index. Our results show that the roughness is able to induce an enhancement on the mean velocity of particles and accelerate the barrier crossing process. The general conditions require a small external force and a small Lévy stability index. We find that with increasing external forces, the enhancement areas of roughness and Lévy stability index both shrink. However, for the Lévy stability index within the enhancement area, its increase will enlarge the enhancement area of roughness. On the contrary, under the same conditions we observe that for a Gaussian noise the roughness always reduces the corresponding mean velocity which is very different from the case of Lévy noise.

Transports in a rough ratchet induced by Lévy noises

We study the transport of a particle subjected to a Lévy noise in a rough ratchet potential which is constructed by superimposing a fast oscillating trigonometric function on a common ratchet background. Due to the superposition of roughness, the transport process exhibits significantly different properties under the excitation of Lévy noises compared to smooth cases. The influence of the roughness on the directional motion is explored by calculating the mean velocities with respect to the Lévy stable index α and the spatial asymmetry parameter q of the ratchet. Variations in the splitting probability have been analyzed to illustrate how roughness affects the transport. In addition, we have examined the influences of roughness on the mean first passage time to know when it accelerates or slows down the first passage process. We find that the roughness can lead to a fast reduction of the absolute value of the mean velocity for small α, however the influence is small for large α. We have illustrated that the ladder-like roughness on the potential wall increases the possibility for particles to cross the gentle side of the ratchet, which results in an increase of the splitting probability to right for the right-skewed ratchet potential. Although the roughness increases the corresponding probability, it does not accelerate the mean first passage process to the right adjacent well. Our results show that the influences of roughness on the mean first passage time are sensitive to the combination of q and α. Hence, the proper q and α can speed up the passage process, otherwise it will slow down it.

Receptor activity-modifying protein dependent and independent activation mechanisms in the coupling of calcitonin gene-related peptide and adrenomedullin receptors to Gs

Calcitonin gene-related peptide (CGRP) or adrenomedullin (AM) receptors are heteromers of the calcitonin receptor-like receptor (CLR), a class B G protein-coupled receptor, and one of three receptor activity-modifying proteins (RAMPs). How CGRP and AM activate CLR and how this process is modulated by RAMPs is unclear. We have defined how CGRP and AM induce Gs-coupling in CLR-RAMP heteromers by measuring the effect of targeted mutagenesis in the CLR transmembrane domain on cAMP production, modeling the active state conformations of CGRP and AM receptors in complex with the Gs C-terminus and conducting molecular dynamics simulations in an explicitly hydrated lipidic bilayer. The largest effects on receptor signaling were seen with H295A^5.40b, I298A^5.43b, L302A^5.47b, N305A^5.50b, L345A^6.49b and E348A^6.52b, F349A^6.53b and H374A^7.47b (class B numbering in superscript). Many of these residues are likely to form part of a group in close proximity to the peptide binding site and link to a network of hydrophilic and hydrophobic residues, which undergo rearrangements to facilitate Gs binding. Residues closer to the extracellular loops displayed more pronounced RAMP or ligand-dependent effects. Mutation of H374^7.47b to alanine increased AM potency 100-fold in the CGRP receptor. The molecular dynamics simulation showed that TM5 and TM6 pivoted around TM3. The data suggest that hydrophobic interactions are more important for CLR activation than other class B GPCRs, providing new insights into the mechanisms of activation of this class of receptor. Furthermore the data may aid in the understanding of how RAMPs modulate the signaling of other class B GPCRs.

Lévy-noise-induced transport in a rough triple-well potential

Rough energy landscape and noisy environment are two common features in many subjects, such as protein folding. Due to the wide findings of bursting or spiking phenomenon in biology science, small diffusions mixing large jumps are adopted to model the noisy environment that can be properly described by Lévy noise. We combine the Lévy noise with the rough energy landscape, modeled by a potential function superimposed by a fast oscillating function, and study the transport of a particle in a rough triple-well potential excited by Lévy noise, rather than only small perturbations. The probabilities of a particle staying in the middle well are considered under different amplitudes of roughness to find out how roughness affects the steady-state probability density function. Variations in the mean first passage time from the middle well to the right well have been investigated with respect to Lévy parameters and amplitudes of the roughness. In addition, we have examined the influences of roughness on the splitting probabilities of the first escape from the middle well. We uncover that the roughness can enhance significantly the first escape of a particle from the middle well, especially for different skewness parameters, but weak differences are found for stability index and noise intensity on the probabilities a particle staying in the middle well and splitting probability to the right.

Structural basis of control of inward rectifier Kir2 channel gating by bulk anionic phospholipids

Phospholipids are required to bind to two distinct sites on the inward rectifier potassium channel for maximal efficacy. Lee et al. show that a membrane-associating tryptophan residue in the second site can mimic the effect of phospholipid binding and cause a conformational change to reveal the primary binding site.

Inward rectifier potassium (Kir) channel activity is controlled by plasma membrane lipids. Phosphatidylinositol-4,5-bisphosphate (PIP₂) binding to a primary site is required for opening of classic inward rectifier Kir2.1 and Kir2.2 channels, but interaction of bulk anionic phospholipid (PL⁻) with a distinct second site is required for high PIP₂ sensitivity. Here we show that introduction of a lipid-partitioning tryptophan at the second site (K62W) generates high PIP₂ sensitivity, even in the absence of PL⁻. Furthermore, high-resolution x-ray crystal structures of Kir2.2[K62W], with or without added PIP₂ (2.8- and 2.0-Å resolution, respectively), reveal tight tethering of the C-terminal domain (CTD) to the transmembrane domain (TMD) in each condition. Our results suggest a refined model for phospholipid gating in which PL⁻ binding at the second site pulls the CTD toward the membrane, inducing the formation of the high-affinity primary PIP₂ site and explaining the positive allostery between PL⁻ binding and PIP₂ sensitivity.

Conformational heterogeneity in closed and open states of the KcsA potassium channel in lipid bicelles

The process of ion channel gating-opening and closing-involves local and global structural changes in the channel in response to external stimuli. Conformational changes depend on the energetic landscape that underlies the transition between closed and open states, which plays a key role in ion channel gating. For the prokaryotic, pH-gated potassium channel KcsA, closed and open states have been extensively studied using structural and functional methods, but the dynamics within each of these functional states as well as the transition between them is not as well understood. In this study, we used solution nuclear magnetic resonance (NMR) spectroscopy to investigate the conformational transitions within specific functional states of KcsA. We incorporated KcsA channels into lipid bicelles and stabilized them into a closed state by using either phosphatidylcholine lipids, known to favor the closed channel, or mutations designed to trap the channel shut by disulfide cross-linking. A distinct state, consistent with an open channel, was uncovered by the addition of cardiolipin lipids. Using selective amino acid labeling at locations within the channel that are known to move during gating, we observed at least two different slowly interconverting conformational states for both closed and open channels. The pH dependence of these conformations and the predictable disruptions to this dependence observed in mutant channels with altered pH sensing highlight the importance of conformational heterogeneity for KcsA gating.

Drug trapping in hERG K+ channels: (not) a matter of drug size?

Inhibition of hERG K+ channels by structurally diverse drugs prolongs the ventricular action potential and increases the risk of torsade de pointes arrhythmias and sudden cardiac death. The capture of drugs behind closed channel gates, so-called drug trapping, is suggested to harbor an increased pro-arrhythmic risk. In this study, the trapping mechanisms of a trapped hERG blocker propafenone and a bulky derivative (MW: 647.24 g mol-1) were studied by making use of electrophysiological measurements in combination with molecular dynamics simulations. Our study suggests that the hERG cavity is able to accommodate very bulky compounds without disturbing gate closure.

H95 Is a pH-Dependent Gate in Aquaporin 4

Aquaporin 4 (AQP4) is a transmembrane protein from the aquaporin family and is the predominant water channel in the mammalian brain. The regulation of permeability of this protein could be of potential therapeutic use to treat various forms of damage to the nervous tissue. In this work, based on data obtained from in silico and in vitro studies, a pH sensitivity that regulates the osmotic water permeability of AQP4 is demonstrated. The results indicate that AQP4 has increased water permeability at conditions of low pH in atomistic computer simulations and experiments carried out on Xenopus oocytes expressing AQP4. With molecular dynamics simulations, this effect was traced to a histidine residue (H95) located in the cytoplasmic lumen of AQP4. A mutant form of AQP4, in which H95 was replaced with an alanine (H95A), loses sensitivity to cytoplasmic pH changes in in vitro osmotic water permeability, thereby substantiating the in silico work.

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.

Allostery without a conformational change? Revisiting the paradigm

Classically, allostery induces a functional switch through a conformational change. However, lately an increasing number of studies concluded that the allostery they observe takes place through sheer dynamics. Here we explain that even if a structural comparison between the active and inactive states does not detect a conformational change, it does not mean that there is no conformational change. We list likely reasons for this lack of observation, including crystallization conditions and crystal effects; one of the states is disordered; the structural comparisons disregard the quaternary protein structure; overlooking synergy effects among allosteric effectors and graded incremental switches and too short molecular dynamics simulations. Specific functions are performed by distinct conformations; they emerge through specific interactions between conformationally selected states.

Ion permeation in K⁺ channels occurs by direct Coulomb knock-on

Potassium channels selectively conduct K(+) ions across cellular membranes with extraordinary efficiency. Their selectivity filter exhibits four binding sites with approximately equal electron density in crystal structures with high K(+) concentrations, previously thought to reflect a superposition of alternating ion- and water-occupied states. Consequently, cotranslocation of ions with water has become a widely accepted ion conduction mechanism for potassium channels. By analyzing more than 1300 permeation events from molecular dynamics simulations at physiological voltages, we observed instead that permeation occurs via ion-ion contacts between neighboring K(+) ions. Coulomb repulsion between adjacent ions is found to be the key to high-efficiency K(+) conduction. Crystallographic data are consistent with directly neighboring K(+) ions in the selectivity filter, and our model offers an intuitive explanation for the high throughput rates of K(+) channels.