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.

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.

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.

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.

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.

K+ congeners that do not compromise Na+ activation of the Na+,K+-ATPase: hydration of the ion binding cavity likely controls ion selectivity.

The Na(+),K(+)-ATPase is essential for ionic homeostasis in animal cells. The dephosphoenzyme contains Na(+) selective inward facing sites, whereas the phosphoenzyme contains K(+) selective outward facing sites. Under normal physiological conditions, K(+) inhibits cytoplasmic Na(+) activation of the enzyme. Acetamidinium (Acet(+)) and formamidinium (Form(+)) have been shown to permeate the pump through the outward facing sites. Here, we show that these cations, unlike K(+), are unable to enter the inward facing sites in the dephosphorylated enzyme. Consistently, the organic cations exhibited little to no antagonism to cytoplasmic Na(+) activation. Na(+),K(+)-ATPase structures revealed a previously undescribed rotamer transition of the hydroxymethyl side chain of the absolutely conserved Thr(772) of the α-subunit. The side chain contributes its hydroxyl to Na(+) in site I in the E1 form and rotates to contribute its methyl group toward K(+) in the E2 form. Molecular dynamics simulations to the E1·AlF4 (-)·ADP·3Na(+) structure indicated that 1) bound organic cations differentially distorted the ion binding sites, 2) the hydroxymethyl of Thr(772) rotates to stabilize bound Form(+) through water molecules, and 3) the rotamer transition is mediated by water traffic into the ion binding cavity. Accordingly, dehydration induced by osmotic stress enhanced the interaction of the congeners with the outward facing sites and profoundly modified the organization of membrane domains of the α-subunit. These results assign a catalytic role for water in pump function, and shed light on a backbone-independent but a conformation-dependent switch between H-bond and dispersion contact as part of the catalytic mechanism of the Na(+),K(+)-ATPase.

The N Terminus of Sarcolipin Plays an Important Role in Uncoupling Sarco-endoplasmic Reticulum Ca2+-ATPase (SERCA) ATP Hydrolysis from Ca2+ Transport.

The sarcoendoplasmic reticulum Ca(2+)-ATPase (SERCA) is responsible for intracellular Ca(2+) homeostasis. SERCA activity in muscle can be regulated by phospholamban (PLB), an affinity modulator, and sarcolipin (SLN), an uncoupler. Although PLB gets dislodged from Ca(2+)-bound SERCA, SLN continues to bind SERCA throughout its kinetic cycle and promotes uncoupling of Ca(2+) transport from ATP hydrolysis. To determine the structural regions of SLN that mediate uncoupling of SERCA, we employed mutagenesis and generated chimeras of PLB and SLN. In this study we demonstrate that deletion of SLN N-terminal residues (2)ERSTQ leads to loss of the uncoupling function even though the truncated peptide can target and constitutively bind SERCA. Furthermore, molecular dynamics simulations of SLN and SERCA interaction showed a rearrangement of SERCA residues that is altered when the SLN N terminus is deleted. Interestingly, transfer of the PLB cytosolic domain to the SLN transmembrane (TM) and luminal tail causes the chimeric protein to lose SLN-like function. Further introduction of the PLB TM region into this chimera resulted in conversion to full PLB-like function. We also found that swapping PLB N and C termini with those from SLN caused the resulting chimera to acquire SLN-like function. Swapping the C terminus alone was not sufficient for this conversion. These results suggest that domains can be switched between SLN and PLB without losing the ability to regulate SERCA activity; however, the resulting chimeras acquire functions different from the parent molecules. Importantly, our studies highlight that the N termini of SLN and PLB influence their respective unique functions.

Membrane accessibility of glutathione.

Regulation of the ion pumping activity of the Na+,K+-ATPase is crucial to the survival of animal cells. Recent evidence has suggested that the activity of the enzyme could be controlled by glutathionylation of cysteine residue 45 of the ß-subunit. Crystal structures so far available indicate that this cysteine is in a transmembrane domain of the protein. Here we have analysed via fluorescence and NMR spectroscopy as well as molecular dynamics simulations whether glutathione is able to penetrate into the interior of a lipid membrane. No evidence for any penetration of glutathione into the membrane was found. Therefore, the most likely mechanism whereby the cysteine residue could become glutathionylated is via a loosening of the α-ß subunit association, creating a hydrophilic passageway between them to allow access of glutathione to the cysteine residue. By such a mechanism, glutathionylation of the protein would be expected to anchor the modified cysteine residue in a hydrophilic environment, inhibiting further motion of the ß-subunit during the enzyme's catalytic cycle and suppressing enzymatic activity, as has been experimentally observed. The results obtained, therefore, suggest a possible structural mechanism of how the Na+,K+-ATPase could be regulated by glutathione.

Molecular mechanism of Na(+),K(+)-ATPase malfunction in mutations characteristic of adrenal hypertension.

Mutations within ion-transporting proteins may severely affect their ability to traffic ions properly and thus perturb the delicate balance of ion gradients. Somatic gain-of-function mutations of the Na(+),K(+)-ATPase α1-subunit have been found in aldosterone-producing adenomas that are among the causes of hypertension. We used molecular dynamics simulations to investigate the structural consequences of these mutations, namely, Leu97 substitution by Arg (L97R), Val325 substitution by Gly (V325G), deletion of residues 93-97 (Del93-97), and deletion-substitution of residues 953-956 by Ser (EETA956S), which shows inward leak currents under physiological conditions. The first three mutations affect the structural context of the key ion-binding residue Glu327 at binding site II, which leads to the loss of the ability to bind ions correctly and to occlude the pump. The mutated residue in L97R is more hydrated, which ultimately leads to the observed proton leak. V325G mimics the structural behavior of L97R; however, it does not promote the hydration of surrounding residues. In Del93-97, a broader opening is observed because of the rearrangement of the kinked transmembrane helix 1, M1, which may explain the sodium leak measured with the mutant. The last mutant, EETA956S, opens an additional water pathway near the C-terminus, affecting the III sodium-specific binding site. The results are in excellent agreement with recent electrophysiology measurements and suggest how three mutations prevent the occlusion of the Na(+),K(+)-ATPase, with the possibility of transforming the pump into a passive ion channel, whereas the fourth mutation provides insight into the sodium binding in the E1 state.

Conformations of double-headed, triple-tailed phospholipid oxidation lipid products in model membranes.

Products of phospholipid oxidation can produce lipids with a carbonyl moiety at the end of a shortened lipid acyl tail, such as 1-palmitoyl-2-(5-oxovaleroyl)-sn-glycero-3-phosphocholine (POVPC). The carbonyl tail of POVPC can covalently bond to the free tertiary amine of a phosphatidylethanolamine lipid in a Schiff base reaction to form a conjugate lipid (SCH) with two head groups, and three acyl tails. We investigate the conformations and properties of this unique class of adduct lipids using molecular dynamics simulations, and show that their insertion into lipid bilayers of POPC increases the average cross-sectional area per lipid and decreases bilayer thickness. Significant increase in acyl tail fluidity is only observed at 25% SCH concentration. The SCH occupies a larger area per lipid than expected for a lipid with three acyl tails, owing to the interfacial location of the long spacer between the two head groups of the SCH. Schiff base formation of lipids can alter the concentration, homeostasis and localizations of phosphatidylserine and phosphatidylethanol lipids in membranes, and can therefore influence several membrane-associated processes including fusion and budding. The current work provides the first detailed structural model of this unique new class of lipids that may have important roles to play in modulating membrane properties and cell physiology.

Reinforcing the membrane-mediated mechanism of action of the anti-tuberculosis candidate drug thioridazine with molecular simulations.

Thioridazine is a well-known dopamine-antagonist drug with a wide range of pharmacological properties ranging from neuroleptic to antimicrobial and even anticancer activity. Thioridazine is a critical component of a promising multi-drug therapy against M. tuberculosis. Amongst the various proposed mechanisms of action, the cell membrane-mediated one is peculiarly tempting due to the distinctive feature of phenothiazine drug family to accumulate in selected body tissues. In this study, we employ long-scale molecular dynamics simulations to investigate the interactions of three different concentrations of thioridazine with zwitterionic and negatively charged model lipid membranes. Thioridazine partitions into the interfacial region of membranes and modifies their structural and dynamic properties, however dissimilarly so at the highest membrane-occurring concentration, that appears to be obtainable only for the negatively charged bilayer. We show that the origin of such changes is the drug induced decrease of the interfacial tension, which ultimately leads to the significant membrane expansion. Our findings support the hypothesis that the phenothiazines therapeutic activity may arise from the drug-membrane interactions, and reinforce the wider, emerging view of action of many small, bioactive compounds.

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.

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.

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.

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.

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 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.

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.

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.