Effect of two activators on the gating of a K2P channel.

TWIK-related potassium channel 1 (TREK1), a two-pore-domain mammalian potassium (K+) channel, regulates the resting potential across cell membranes, presenting a promising therapeutic target for neuropathy treatment. The gating of this channel converges in the conformation of the narrowest part of the pore: the selectivity filter (SF). Various hypotheses explain TREK1 gating modulation, including the dynamics of loops connecting the SF with transmembrane helices and the stability of hydrogen bond (HB) networks adjacent to the SF. Recently, two small molecules (Q6F and Q5F) were reported as activators that affect TREK1 by increasing its open probability in single-channel current measurements. Here, using molecular dynamics simulations, we investigate the effect of these ligands on the previously proposed modulation mechanisms of TREK1 gating compared to the apo channel. Our findings reveal that loop dynamics at the upper region of the SF exhibit only a weak correlation with permeation events/nonpermeation periods, whereas the HB network behind the SF appears more correlated. These nonpermeation periods arise from both distinct mechanisms: a C-type inactivation (resulting from dilation at the top of the SF), which has been described previously, and a carbonyl flipping in an SF binding site. We find that, besides the prevention of C-type inactivation in the channel, the ligands increase the probability of permeation by modulating the dynamics of the carbonyl flipping, influenced by a threonine residue at the bottom of the SF. These results offer insights for rational ligand design to optimize the gating modulation of TREK1 and related K+ channels.

Biomolecular simulations at the exascale: From drug design to organelles and beyond.

The rapid advancement in computational power available for research offers to bring not only quantitative improvements, but also qualitative changes in the field of biomolecular simulation. Here, we review the state of biomolecular dynamics simulations at the threshold to exascale resources becoming available. Both developments in parallel and distributed computing will be discussed, providing a perspective on the state of the art of both. A main focus will be on obtaining binding and conformational free energies, with an outlook to macromolecular complexes and (sub)cellular assemblies.

Selectivity filter mutations shift ion permeation mechanism in potassium channels.

Potassium (K+) channels combine high conductance with high ion selectivity. To explain this efficiency, two molecular mechanisms have been proposed. The "direct knock-on" mechanism is defined by water-free K+ permeation and formation of direct ion-ion contacts in the highly conserved selectivity filter (SF). The "soft knock-on" mechanism involves co-permeation of water and separation of K+ by water molecules. With the aim to distinguish between these mechanisms, crystal structures of the KcsA channel with mutations in two SF residues-G77 and T75-were published, where the arrangements of K+ ions and water display canonical soft knock-on configurations. These data were interpreted as evidence of the soft knock-on mechanism in wild-type channels. Here, we test this interpretation using molecular dynamics simulations of KcsA and its mutants. We show that while a strictly water-free direct knock-on permeation is observed in the wild type, conformational changes induced by these mutations lead to distinct ion permeation mechanisms, characterized by co-permeation of K+ and water. These mechanisms are characterized by reduced conductance and impaired potassium selectivity, supporting the importance of full dehydration of potassium ions for the hallmark high conductance and selectivity of K+ channels. In general, we present a case where mutations introduced at the critical points of the permeation pathway in an ion channel drastically change its permeation mechanism in a nonintuitive manner.

Interactions between selectivity filter and pore helix control filter gating in the MthK channel.

K+ channel activity can be limited by C-type inactivation, which is likely initiated in part by dissociation of K+ ions from the selectivity filter and modulated by the side chains that surround it. While crystallographic and computational studies have linked inactivation to a "collapsed" selectivity filter conformation in the KcsA channel, the structural basis for selectivity filter gating in other K+ channels is less clear. Here, we combined electrophysiological recordings with molecular dynamics simulations, to study selectivity filter gating in the model potassium channel MthK and its V55E mutant (analogous to KcsA E71) in the pore-helix. We found that MthK V55E has a lower open probability than the WT channel, due to decreased stability of the open state, as well as a lower unitary conductance. Simulations account for both of these variables on the atomistic scale, showing that ion permeation in V55E is altered by two distinct orientations of the E55 side chain. In the "vertical" orientation, in which E55 forms a hydrogen bond with D64 (as in KcsA WT channels), the filter displays reduced conductance compared to MthK WT. In contrast, in the "horizontal" orientation, K+ conductance is closer to that of MthK WT; although selectivity filter stability is lowered, resulting in more frequent inactivation. Surprisingly, inactivation in MthK WT and V55E is associated with a widening of the selectivity filter, unlike what is observed for KcsA and reminisces recent structures of inactivated channels, suggesting a conserved inactivation pathway across the potassium channel family.

Direct Detection of Bound Ammonium Ions in the Selectivity Filter of Ion Channels by Solid-State NMR.

The flow of ions across cell membranes facilitated by ion channels is an important function for all living cells. Despite the huge amount of structural data provided by crystallography, elucidating the exact interactions between the selectivity filter atoms and bound ions is challenging. Here, we detect bound 15N-labeled ammonium ions as a mimic for potassium ions in ion channels using solid-state NMR under near-native conditions. The non-selective ion channel NaK showed two ammonium peaks corresponding to its two ion binding sites, while its potassium-selective mutant NaK2K that has a signature potassium-selective selectivity filter with four ion binding sites gave rise to four ammonium peaks. Ions bound in specific ion binding sites were identified based on magnetization transfer between the ions and carbon atoms in the selectivity filters. Magnetization transfer between bound ions and water molecules revealed that only one out of four ions in the selectivity filter of NaK2K is in close contact with water, which is in agreement with the direct knock-on ion conduction mechanism where ions are conducted through the channel by means of direct interactions without water molecules in between. Interestingly, the potassium-selective ion channels investigated here (NaK2K and, additionally, KcsA-Kv1.3) showed remarkably different chemical shifts for their bound ions, despite having identical amino acid sequences and crystal structures of their selectivity filters. Molecular dynamics simulations show similar ion binding and conduction behavior between ammonium and potassium ions and identify the origin of the differences between the investigated potassium channels.

Structural plasticity of the selectivity filter in a nonselective ion channel.

The sodium potassium ion channel (NaK) is a nonselective ion channel that conducts both sodium and potassium across the cellular membrane. A new crystallographic structure of NaK reveals conformational differences in the residues that make up the selectivity filter between the four subunits that form the ion channel and the inner helix of the ion channel. The crystallographic structure also identifies a side-entry, ion-conduction pathway for Na+ permeation that is unique to NaK. NMR studies and molecular dynamics simulations confirmed the dynamical nature of the top part of the selectivity filter and the inner helix in NaK as also observed in the crystal structure. Taken together, these results indicate that the structural plasticity of the selectivity filter combined with the dynamics of the inner helix of NaK are vital for the efficient conduction of different ions through the non-selective ion channel of NaK.

The Persistent Question of Potassium Channel Permeation Mechanisms.

Potassium channels play critical roles in many physiological processes, providing a selective permeation route for K+ ions in and out of a cell, by employing a carefully designed selectivity filter, evolutionarily conserved from viruses to mammals. The structure of the selectivity filter was determined at atomic resolution by x-ray crystallography, showing a tight coordination of desolvated K+ ions by the channel. However, the molecular mechanism of K+ ions permeation through potassium channels remains unclear, with structural, functional and computational studies often providing conflicting data and interpretations. In this review, we will present the proposed mechanisms, discuss their origins, and will critically assess them against all available data. General properties shared by all potassium channels are introduced first, followed by the introduction of two main mechanisms of ion permeation: soft and direct knock-on. Then, we will discuss critical computational and experimental studies that shaped the field. We will especially focus on molecular dynamics (MD) simulations, that provided mechanistic and energetic aspects of K+ permeation, but at the same time created long-standing controversies. Further challenges and possible solutions are presented as well.

Polycation-Anionic Lipid Membrane Interactions.

Natural or synthetic polycations are used as biocides or as drug/gene carriers. Understanding the interactions between these macromolecules and cell membranes at the molecular level is therefore of great importance for the design of effective polymer biocides or biocompatible polycation-based delivery systems. Until now, details of the processes at the interface between polycations and biological systems have not been fully recognized. In this study, we consider the effect of strong polycations with quaternary ammonium groups on the properties of anionic lipid membranes that we use as a model system for protein-free cell membranes. For this purpose, we employed experimental measurements and atomic-scale molecular dynamics (MD) simulations. MD simulations reveal that the polycations are strongly hydrated in the aqueous phase and do not lose the water shell after adsorption at the bilayer surface. As a result of strong hydration, the polymer chains reside at the phospholipid headgroup and do not penetrate to the acyl chain region. The polycation adsorption involves the formation of anionic lipid-rich domains, and the density of anionic lipids in these domains depends on the length of the polycation chain. We observed the accumulation of anionic lipids only in the leaflet interacting with the polymer, which leads to the formation of compositionally asymmetric domains. Asymmetric adsorption of the polycation on only one leaflet of the anionic membrane strongly affects the membrane properties in the polycation-membrane contact areas: (i) anionic lipid accumulates in the region near the adsorbed polymer, (ii) acyl chain ordering and lipid packing are reduced, which results in a decrease in the thickness of the bilayer, and (iii) polycation-anionic membrane interactions are strongly influenced by the presence and concentration of salt. Our results provide an atomic-scale description of the interactions of polycations with anionic lipid bilayers and are fully supported by the experimental data. The outcomes are important for understanding the correlation of the structure of polycations with their activity on biomembranes.

The structure of a potassium-selective ion channel reveals a hydrophobic gate regulating ion permeation.

Protein dynamics are essential to function. One example of this is the various gating mechanisms within ion channels, which are transmembrane proteins that act as gateways into the cell. Typical ion channels switch between an open and closed state via a conformational transition which is often triggered by an external stimulus, such as ligand binding or pH and voltage differences. The atomic resolution structure of a potassium-selective ion channel named NaK2K has allowed us to observe that a hydro-phobic residue at the bottom of the selectivity filter, Phe92, appears in dual conformations. One of the two conformations of Phe92 restricts the diameter of the exit pore around the selectivity filter, limiting ion flow through the channel, while the other conformation of Phe92 provides a larger-diameter exit pore from the selectivity filter. Thus, it can be concluded that Phe92 acts as a hydro-phobic gate, regulating the flow of ions through the selectivity filter.

A ß-barrel for oil transport through lipid membranes: Dynamic NMR structures of AlkL.

The protein AlkL is known to increase permeability of the outer membrane of bacteria for hydrophobic molecules, yet the mechanism of transport has not been determined. Differing crystal and NMR structures of homologous proteins resulted in a controversy regarding the degree of structure and the role of long extracellular loops. Here we solve this controversy by determining the de novo NMR structure in near-native lipid bilayers, and by accessing structural dynamics relevant to hydrophobic substrate permeation through molecular-dynamics simulations and by characteristic NMR relaxation parameters. Dynamic lateral exit sites large enough to accommodate substrates such as carvone or octane occur through restructuring of a barrel extension formed by the extracellular loops.

Molecular mechanism of a potassium channel gating through activation gate-selectivity filter coupling.

Potassium channels are presumed to have two allosterically coupled gates, the activation gate and the selectivity filter gate, that control channel opening, closing, and inactivation. However, the molecular mechanism of how these gates regulate K+ ion flow through the channel remains poorly understood. An activation process, occurring at the selectivity filter, has been recently proposed for several potassium channels. Here, we use X-ray crystallography and extensive molecular dynamics simulations, to study ion permeation through a potassium channel MthK, for various opening levels of both gates. We find that the channel conductance is controlled at the selectivity filter, whose conformation depends on the activation gate. The crosstalk between the gates is mediated through a collective motion of channel helices, involving hydrophobic contacts between an isoleucine and a conserved threonine in the selectivity filter. We propose a gating model of selectivity filter-activated potassium channels, including pharmacologically relevant two-pore domain (K2P) and big potassium (BK) channels.

The conduction pathway of potassium channels is water free under physiological conditions.

Ion conduction through potassium channels is a fundamental process of life. On the basis of crystallographic data, it was originally proposed that potassium ions and water molecules are transported through the selectivity filter in an alternating arrangement, suggesting a "water-mediated" knock-on mechanism. Later on, this view was challenged by results from molecular dynamics simulations that revealed a "direct" knock-on mechanism where ions are in direct contact. Using solid-state nuclear magnetic resonance techniques tailored to characterize the interaction between water molecules and the ion channel, we show here that the selectivity filter of a potassium channel is free of water under physiological conditions. Our results are fully consistent with the direct knock-on mechanism of ion conduction but contradict the previously proposed water-mediated knock-on mechanism.

Conotoxin κM-RIIIJ, a tool targeting asymmetric heteromeric Kv1 channels.

The vast complexity of native heteromeric K+ channels is largely unexplored. Defining the composition and subunit arrangement of individual subunits in native heteromeric K+ channels and establishing their physiological roles is experimentally challenging. Here we systematically explored this "zone of ignorance" in molecular neuroscience. Venom components, such as peptide toxins, appear to have evolved to modulate physiologically relevant targets by discriminating among closely related native ion channel complexes. We provide proof-of-principle for this assertion by demonstrating that κM-conotoxin RIIIJ (κM-RIIIJ) from Conus radiatus precisely targets "asymmetric" Kv channels composed of three Kv1.2 subunits and one Kv1.1 or Kv1.6 subunit with 100-fold higher apparent affinity compared with homomeric Kv1.2 channels. Our study shows that dorsal root ganglion (DRG) neurons contain at least two major functional Kv1.2 channel complexes: a heteromer, for which κM-RIIIJ has high affinity, and a putative Kv1.2 homomer, toward which κM-RIIIJ is less potent. This conclusion was reached by (i) covalent linkage of members of the mammalian Shaker-related Kv1 family to Kv1.2 and systematic assessment of the potency of κM-RIIIJ block of heteromeric K+ channel-mediated currents in heterologous expression systems; (ii) molecular dynamics simulations of asymmetric Kv1 channels providing insights into the molecular basis of κM-RIIIJ selectivity and potency toward its targets; and (iii) evaluation of calcium responses of a defined population of DRG neurons to κM-RIIIJ. Our study demonstrates that bioactive molecules present in venoms provide essential pharmacological tools that systematically target specific heteromeric Kv channel complexes that operate in native tissues.

Molecular Mechanism of Polycation-Induced Pore Formation in Biomembranes.

Polycations are an attractive class of macromolecules with promising applications as drug/gene carriers and biocides. The chemical structure and concentration of a polycation determine its interaction with cellular membranes and, hence, are crucial parameters for designing efficient nontoxic polycations. However, the interaction of polycations with biomembranes at the molecular level and the corresponding free-energy landscape is not well understood. In this work, we investigate the molecular mechanism of interaction between a strong polycation substituted with alkyl moieties and zwitterionic membranes via long-time-scale all-atom molecular dynamics simulations and free-energy calculations combined with Langmuir monolayer, atomic force microscopy, and calcein-release experimental measurements. We found that the membrane activity of the polycation and its ability to induce pores in the membranes can be attributed to the polycation-induced changes in the bilayer organization, such as reduced membrane thickness, increased disorder of the acyl chains, reduced packing, and electrostatic field gradients between membrane leaflets. These changes facilitate the penetration of water into the membrane and the formation of aqueous defects/pores. The calculated free-energy profiles indicate that the polycation lowers the nucleation barrier for pore opening and the free energy for pore formation in a concentration-dependent manner. Above the critical coverage of the membrane, the polycation nucleates spontaneous pores in zwitterionic membranes. Our work demonstrates the potential of combining enhanced sampling methods in MD simulations with experiments for a quantitative description of various events in the polycation-membrane interaction cycle, such as strong adsorption on the membrane due to hydrophobic and electrostatic interactions, and pore formation.

K+ binding and proton redistribution in the E2P state of the H+, K+-ATPase.

The H+, K+-ATPase (HKA) uses ATP to pump protons into the gastric lumen against a million-fold proton concentration gradient while counter-transporting K+ from the lumen. The mechanism of release of a proton into a highly acidic stomach environment, and the subsequent binding of a K+ ion necessitates a network of protonable residues and dynamically changing protonation states in the cation binding pocket dominated by five acidic amino acid residues E343, E795, E820, D824, and D942. We perform molecular dynamics simulations of spontaneous K+ binding to all possible protonation combinations of the acidic amino acids and carry out free energy calculations to determine the optimal protonation state of the luminal-open E2P state of the pump which is ready to bind luminal K+. A dynamic pKa correlation analysis reveals the likelihood of proton transfer events within the cation binding pocket. In agreement with in-vitro measurements, we find that E795 is likely to be protonated, and that E820 is at the center of the proton transfer network in the luminal-open E2P state. The acidic residues D942 and D824 are likely to remain protonated, and the proton redistribution occurs predominantly amongst the glutamate residues exposed to the lumen. The analysis also shows that a lower number of K+ ions bind at lower pH, modeled by a higher number of protons in the cation binding pocket, in agreement with the 'transport stoichiometry variation' hypothesis.

Direct knock-on of desolvated ions governs strict ion selectivity in K+ channels.

The seeming contradiction that K+ channels conduct K+ ions at maximal throughput rates while not permeating slightly smaller Na+ ions has perplexed scientists for decades. Although numerous models have addressed selective permeation in K+ channels, the combination of conduction efficiency and ion selectivity has not yet been linked through a unified functional model. Here, we investigate the mechanism of ion selectivity through atomistic simulations totalling more than 400 µs in length, which include over 7,000 permeation events. Together with free-energy calculations, our simulations show that both rapid permeation of K+ and ion selectivity are ultimately based on a single principle: the direct knock-on of completely desolvated ions in the channels' selectivity filter. Herein, the strong interactions between multiple 'naked' ions in the four filter binding sites give rise to a natural exclusion of any competing ions. Our results are in excellent agreement with experimental selectivity data, measured ion interaction energies and recent two-dimensional infrared spectra of filter ion configurations.

In silico assessment of the conduction mechanism of the Ryanodine Receptor 1 reveals previously unknown exit pathways.

The ryanodine receptor 1 is a large calcium ion channel found in mammalian skeletal muscle. The ion channel gained a lot of attention recently, after multiple independent authors published near-atomic cryo electron microscopy data. Taking advantage of the unprecedented quality of structural data, we performed molecular dynamics simulations on the entire ion channel as well as on a reduced model. We calculated potentials of mean force for Ba2+, Ca2+, Mg2+, K+, Na+ and Cl- ions using umbrella sampling to identify the key residues involved in ion permeation. We found two main binding sites for the cations, whereas the channel is strongly repulsive for chloride ions. Furthermore, the data is consistent with the model that the receptor achieves its ion selectivity by over-affinity for divalent cations in a calcium-block-like fashion. We reproduced the experimental conductance for potassium ions in permeation simulations with applied voltage. The analysis of the permeation paths shows that ions exit the pore via multiple pathways, which we suggest to be related to the experimental observation of different subconducting states.

Correction to: The CAPOS mutation in ATP1A3 alters Na/K-ATPase function and results in auditory neuropathy which has implications for management.

The following information was inadvertently omitted in the original publication.

The CAPOS mutation in ATP1A3 alters Na/K-ATPase function and results in auditory neuropathy which has implications for management.

Cerebellar ataxia, areflexia, pes cavus, optic atrophy and sensorineural hearing impairment (CAPOS) is a rare clinically distinct syndrome caused by a single dominant missense mutation, c.2452G>A, p.Glu818Lys, in ATP1A3, encoding the neuron-specific alpha subunit of the Na+/K+-ATPase α3. Allelic mutations cause the neurological diseases rapid dystonia Parkinsonism and alternating hemiplegia of childhood, disorders which do not encompass hearing or visual impairment. We present detailed clinical phenotypic information in 18 genetically confirmed patients from 11 families (10 previously unreported) from Denmark, Sweden, UK and Germany indicating a specific type of hearing impairment-auditory neuropathy (AN). All patients were clinically suspected of CAPOS and had hearing problems. In this retrospective analysis of audiological data, we show for the first time that cochlear outer hair cell activity was preserved as shown by the presence of otoacoustic emissions and cochlear microphonic potentials, but the auditory brainstem responses were grossly abnormal, likely reflecting neural dyssynchrony. Poor speech perception was observed, especially in noise, which was beyond the hearing level obtained in the pure tone audiograms in several of the patients presented here. Molecular modelling and in vitro electrophysiological studies of the specific CAPOS mutation were performed. Heterologous expression studies of α3 with the p.Glu818Lys mutation affects sodium binding to, and release from, the sodium-specific site in the pump, the third ion-binding site. Molecular dynamics simulations confirm that the structure of the C-terminal region is affected. In conclusion, we demonstrate for the first time evidence for auditory neuropathy in CAPOS syndrome, which may reflect impaired propagation of electrical impulses along the spiral ganglion neurons. This has implications for diagnosis and patient management. Auditory neuropathy is difficult to treat with conventional hearing aids, but preliminary improvement in speech perception in some patients suggests that cochlear implantation may be effective in CAPOS patients.

Effect of Polycation Structure on Interaction with Lipid Membranes.

Interaction of polycations with lipid membranes is a very important issue in many biological and medical applications such as gene delivery or antibacterial usage. In this work, we address the influence of hydrophobic substitution of strong polycations containing quaternary ammonium groups on the polymer-zwitterionic membrane interactions. In particular, we focus on the polymer tendency to adsorb on or/and incorporate into the membrane. We used complementary experimental and computational methods to enhance our understanding of the mechanism of the polycation-membrane interactions. Polycation adsorption on liposomes was assessed using dynamic light scattering (DLS) and zeta potential measurements. The ability of the polymers to form hydrophilic pores in the membrane was evaluated using a calcein-release method. The polymer-membrane interaction at the molecular scale was explored by performing atomistic molecular dynamics (MD) simulations. Our results show that the length of the alkyl side groups plays an essential role in the polycation adhesion on the zwitterionic surface, while the degree of substitution affects the polycation ability to incorporate into the membrane. Both the experimental and computational results show that the membrane permeability can be dramatically affected by the amount of alkyl side groups attached to the polycation main chain.