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

The Molecular Mechanism of Ion Selectivity in Nanopores

Ion channels exhibit strong selectivity for specific ions over others under electrochemical potentials, such as KcsA for K+ over Na+. Based on the thermodynamic analysis, this study is focused on exploring the mechanism of ion selectivity in nanopores. It is well known that ions must lose part of their hydration layer to enter the channel. Therefore, the ion selectivity of a channel is due to the rearrangement of water molecules when entering the nanopore, which may be related to the hydrophobic interactions between ions and channels. In our recent works on hydrophobic interactions, with reference to the critical radius of solute (Rc), it was divided into initial and hydrophobic solvation processes. Additionally, the different dissolved behaviors of solutes in water are expected in various processes, such as dispersed and accumulated distributions in water. Correspondingly, as the ion approaches the nanopore, there seems to exist the "repulsive" or "attractive" forces between them. In the initial process ( Collapse

Key Words

• hydrogen bondings
• hydrophobic interactions
• ion selectivity
• nanopore
• water

Collapse

MESH Headings Collapse

Grants Collapse

Affiliation(s)

Yan-Nan Chen
Yu-Zhen Liu
Qiang Sun Key Laboratory of Orogenic Belts and Crustal Evolution, The School of Earth and Space Sciences, Ministry of Education, Peking University, Beijing 100871, China; (Y.-N.C.); (Y.-Z.L.)

Collapse

Effect of THz Waves of Different Orientations on K+ Permeation Efficiency in the KcsA Channel

Potassium (K) channels show the highest variability and most frequent alterations in expression in many tumor types, and modulation of K+ channels may represent a new window for cancer therapy. In previous work, we found that a terahertz (THz) field incident along the z-axis with a frequency of 51.87 THz increased the ion flux through K+ channels. In practice, it is difficult to ensure that the incident electromagnetic (EM) wave is strictly parallel to the direction of channel ion flow. In this paper, we found by changing the direction of the applied electric field that the EM wave of a specific frequency has the largest ion flux when the incident direction is along the ion flow, and the smallest ion flux when the incident direction is perpendicular to the ion flow, and that overall the EM wave of this frequency enhances the ion flow of the K+ channel. Changes in the direction of the applied field at a specific frequency affect the stability of the φ dihedral angle of the GLY77 residue and alter the ion permeation mechanism in the selectivity filter (SF) region, thus affecting the ion flux. Therefore, this frequency can be used to modulate K+ fluxes by THz waves to cause rapid apoptosis in potassium-overloaded tumor cells. This approach consequently represents an important tool for the treatment of cancer and is expected to be applied in practical therapy.

Regulation of Ion Permeation of the KcsA Channel by Applied Midinfrared Field

Ion transport molecules are involved in many physiological and pathological processes and are considered potential targets for cancer treatment. In the large family of ion transport molecules, potassium (K) ion channels, as surface-expressed proteins, show the highest variability and most frequent expression changes in many tumor types. The key to exploring the permeation of K⁺ through potassium channels lies in the conserved sequence TVGYG, which is common in the selectivity filter (SF) region of all potassium channels. We found that the K⁺ flux significantly increased with the help of a specific frequency terahertz electromagnetic wave (51.87 THz) in the KcsA channel using a molecular dynamics combined model through the combined simulation of the constant electric field method and ion imbalance method. This frequency has the strongest absorption peak in the infrared spectrum of -C=O groups in the SF region. With the applied electric field of 51.87 THz, the Y78 residue at the S₁ site of the SF has a smaller vibration amplitude and a more stable structure, which enables the K⁺ to bind closely with the carbonyl oxygen atoms in the SF and realize ion conduction in a more efficient direct Coulomb knock-on.

Hydrated Anions: From Clusters to Bulk Solution with Quasi-Chemical Theory

The interactions of hydrated ions with molecular and macromolecular solution and interface partners are strong on a chemical energy scale. Here, we recount the foremost ab initio theory for the evaluation of the hydration free energies of ions, namely, quasi-chemical theory (QCT). We focus on anions, particularly halides but also the hydroxide anion, because they have been outstanding challenges for all theories. For example, this work supports understanding the high selectivity for F^- over Cl^- in fluoride-selective ion channels despite the identical charge and the size similarity of these ions. QCT is built by the identification of inner-shell clusters, separate treatment of those clusters, and then the integration of those results into the broader-scale solution environment. Recent work has focused on a close comparison with mass-spectrometric measurements of ion-hydration equilibria. We delineate how ab initio molecular dynamics (AIMD) calculations on ion-hydration clusters, elementary statistical thermodynamics, and electronic structure calculations on cluster structures sampled from the AIMD calculations obtain just the free energies extracted from the cluster experiments. That theory-experiment comparison has not been attempted before the work discussed here, but the agreement is excellent with moderate computational effort. This agreement reinforces both theory and experiment and provides a numerically accurate inner-shell contribution to QCT. The inner-shell complexes involving heavier halides display strikingly asymmetric hydration clusters. Asymmetric hydration structures can be problematic for the evaluation of the QCT outer-shell contribution with the polarizable continuum model (PCM). Nevertheless, QCT provides a favorable setting for the exploitation of PCM when the inner-shell material shields the ion from the outer solution environment. For the more asymmetrically hydrated, and thus less effectively shielded, heavier halide ions clustered with waters, the PCM is less satisfactory. We therefore investigate an inverse procedure in which the inner-shell structures are sampled from readily available AIMD calculations on the bulk solutions. This inverse procedure is a remarkable improvement; our final results are in close agreement with a standard tabulation of hydration free energies, and the final composite results are independent of the coordination number on the chemical energy scale of relevance, as they should be. Finally, a comparison of anion hydration structure in clusters and bulk solutions from AIMD simulations emphasize some differences: the asymmetries of bulk solution inner-shell structures are moderated compared with clusters but are still present, and inner hydration shells fill to slightly higher average coordination numbers in bulk solution than in clusters.

Weak Cation Selectivity in HCN Channels Results From K+-Mediated Release of Na+ From Selectivity Filter Binding Sites

Hyperpolarization-activated cyclic nucleotide-gated (HCN) channels generate the pacemaker current which plays an important role in the timing of various biological processes like the heart beat. We used umbrella sampling to explore the potential of mean force for the conduction of potassium and sodium through the open HCN4 pore. Our data explain distinct functional features like low unitary conductance and weak selectivity as a result of high energetic barriers inside the selectivity filter of this channel. They exceed the 3-5 kJ/mol threshold which is presumed as maximal barrier for diffusion-limited conductance. Furthermore, simulations provide a thermodynamic explanation for the weak cation selectivity of HCN channels that contain only two ion binding sites in the selectivity filter (SF). We find that sodium ions bind more strongly to the SF than potassium and are easier released by binding of potassium than of another sodium. Hence ion transport and selectivity in HCN channels is not determined by the same mechanism as in potassium-selective channels; it rather relies on sodium as a weak blocker that can only be released by potassium.

Molecular Modeling is an Enabling Approach to Complement and Enhance Channelopathy Research

Hundreds of human membrane proteins form channels that transport necessary ions and compounds, including drugs and metabolites, yet details of their normal function or how function is altered by genetic variants to cause diseases are often unknown. Without this knowledge, researchers are less equipped to develop approaches to diagnose and treat channelopathies. High-resolution computational approaches such as molecular modeling enable researchers to investigate channelopathy protein function, facilitate detailed hypothesis generation, and produce data that is difficult to gather experimentally. Molecular modeling can be tailored to each physiologic context that a protein may act within, some of which may currently be difficult or impossible to assay experimentally. Because many genomic variants are observed in channelopathy proteins from high-throughput sequencing studies, methods with mechanistic value are needed to interpret their effects. The eminent field of structural bioinformatics integrates techniques from multiple disciplines including molecular modeling, computational chemistry, biophysics, and biochemistry, to develop mechanistic hypotheses and enhance the information available for understanding function. Molecular modeling and simulation access 3D and time-dependent information, not currently predictable from sequence. Thus, molecular modeling is valuable for increasing the resolution with which the natural function of protein channels can be investigated, and for interpreting how genomic variants alter them to produce physiologic changes that manifest as channelopathies. © 2022 American Physiological Society. Compr Physiol 12:3141-3166, 2022.

Energy landscape of the SARS-CoV-2 reveals extensive conformational heterogeneity

Cryo-electron microscopy (cryo-EM) has produced a number of structural models of the SARS-CoV-2 spike, already prompting biomedical outcomes. However, these reported models and their associated electrostatic potential maps represent an unknown admixture of conformations stemming from the underlying energy landscape of the spike protein. As with any protein, some of the spike's conformational motions are expected to be biophysically relevant, but cannot be interpreted only by static models. Using experimental cryo-EM images, we present the energy landscape of the glycosylated spike protein, and identify the diversity of low-energy conformations in the vicinity of its open (so called 1RBD-up) state. The resulting atomic refinement reveal global and local molecular rearrangements that cannot be inferred from an average 1RBD-up cryo-EM model. Here we report varied degrees of "openness" in global conformations of the 1RBD-up state, not revealed in the single-model interpretations of the density maps, together with conformations that overlap with the reported models. We discover how the glycan shield contributes to the stability of these low-energy conformations. Five out of six binding sites we analyzed, including those for engaging ACE2, therapeutic mini-proteins, linoleic acid, two different kinds of antibodies, switch conformations between their known apo- and holo-conformations, even when the global spike conformation is 1RBD-up. This apo-to-holo switching is reminiscent of a conformational preequilibrium. We found only one binding site, namely that of AB-C135 remains in apo state within all the sampled free energy-minimizing models, suggesting an induced fit mechanism for the docking of this antibody to the spike.

Ion permeation, selectivity, and electronic polarization in fluoride channels

Fluoride channels (Flucs) export toxic F^- from the cytoplasm. Crystallography and mutagenesis have identified several conserved residues crucial for fluoride transport, but the permeation mechanism at the molecular level has remained elusive. Herein, we have applied constant-pH molecular dynamics and free-energy-sampling methods to investigate fluoride permeation through a Fluc protein from Escherichia coli. We find that fluoride is facile to permeate in its charged form, i.e., F^-, by traversing through a non-bonded network. The extraordinary F^- selectivity is gained by the hydrogen-bonding capability of the central binding site and the Coulombic filter at the channel entrance. The F^- permeation rate calculated using an electronically polarizable force field is significantly more accurate compared with the experimental value than that calculated using a more standard additive force field, suggesting an essential role for electronic polarization in the F^--Fluc interactions.

Computational methods and theory for ion channel research

Ion channels are fundamental biological devices that act as gates in order to ensure selective ion transport across cellular membranes; their operation constitutes the molecular mechanism through which basic biological functions, such as nerve signal transmission and muscle contraction, are carried out. Here, we review recent results in the field of computational research on ion channels, covering theoretical advances, state-of-the-art simulation approaches, and frontline modeling techniques. We also report on few selected applications of continuum and atomistic methods to characterize the mechanisms of permeation, selectivity, and gating in biological and model channels.

Thermodynamics of ion binding and occupancy in potassium channels

Potassium channels modulate various cellular functions through efficient and selective conduction of K⁺ ions. The mechanism of ion conduction in potassium channels has recently emerged as a topic of debate. Crystal structures of potassium channels show four K⁺ ions bound to adjacent binding sites in the selectivity filter, while chemical intuition and molecular modeling suggest that the direct ion contacts are unstable. Molecular dynamics (MD) simulations have been instrumental in the study of conduction and gating mechanisms of ion channels. Based on MD simulations, two hypotheses have been proposed, in which the four-ion configuration is an artifact due to either averaged structures or low temperature in crystallographic experiments. The two hypotheses have been supported or challenged by different experiments. Here, MD simulations with polarizable force fields validated by ab initio calculations were used to investigate the ion binding thermodynamics. Contrary to previous beliefs, the four-ion configuration was predicted to be thermodynamically stable after accounting for the complex electrostatic interactions and dielectric screening. Polarization plays a critical role in the thermodynamic stabilities. As a result, the ion conduction likely operates through a simple single-vacancy and water-free mechanism. The simulations explained crystal structures, ion binding experiments and recent controversial mutagenesis experiments. This work provides a clear view of the mechanism underlying the efficient ion conduction and demonstrates the importance of polarization in ion channel simulations.

Polarization shapes the energy landscape of ion conduction in potassium channels.

Potassium and sodium ion complexes with a partial peptide of the selectivity filter in K+ channels studied by cold ion trap infrared spectroscopy: the effect of hydration

Potassium channels allow K+ to rapidly diffuse, while the selectivity filter (SF) actively blocks Na+. The presence of water in the SF during ion translocation remains under debate due to the experimental and computational challenges in characterizing the interactions between water, ions, and the SF. Our bottom-up approach has been applied to a system composed of a partial peptide of the SF (Ac-tyrosine-NHMe) with a metal ion and a single water molecule to probe these interactions. The IR photodissociation spectra of M+Ac-tyrosine-NHMe(H2O) (M = Na, K) combined with quantum chemical calculations revealed that the water molecule binding sites are ion-dependent. In addition, the ion-peptide distances are elongated significantly for the K+ complex in comparison to the Na+ complex by the addition of a single water molecule. This striking structural difference with the water molecule is discussed in relation to ion selectivity and translocation within the K+ channel.

Physics of Selective Conduction and Point Mutation in Biological Ion Channels

We introduce a statistical and linear response theory of selective conduction in biological ion channels with multiple binding sites and possible point mutation. We derive an effective grand-canonical ensemble and generalized Einstein relations for the selectivity filter, assuming strongly coordinated ionic motion, and allowing for ionic Coulomb blockade. The theory agrees well with data from the KcsA K^{+} channel and a mutant. We show that the Eisenman relations for thermodynamic selectivity follow from the condition for fast conduction and find that maximum conduction requires the binding sites to be nearly identical.

Energy Landscape of the SARS-CoV-2 Reveals Extensive Conformational Heterogeneity

Cryo-electron microscopy (cryo-EM) has produced a number of structural models of the SARS-CoV-2 spike, already prompting biomedical outcomes. However, these reported models and their associated electrostatic potential maps represent an unknown admixture of conformations stemming from the underlying energy landscape of the spike protein. As for any protein, some of the spike's conformational motions are expected to be biophysically relevant, but cannot be interpreted only by static models. Using experimental cryo-EM images, we present the energy landscape of the spike protein conformations, and identify molecular rearrangements along the most-likely conformational path in the vicinity of the open (so called 1RBD-up) state. The resulting global and local atomic refinements reveal larger movements than those expected by comparing the reported 1RBD-up and 1RBD-down cryo-EM models. Here we report greater degrees of "openness" in global conformations of the 1RBD-up state, not revealed in the single-model interpretations of the density maps, together with conformations that overlap with the reported models. We discover how the glycan shield contributes to the stability of these conformations along the minimum free-energy pathway. A local analysis of seven key binding pockets reveals that six out them, including those for engaging ACE2, therapeutic mini-proteins, linoleic acid, two different kinds of antibodies, and protein-glycan interaction sites, switch conformations between their known apo- and holo-conformations, even when the global spike conformation is 1RBD-up. This is reminiscent of a conformational pre-equilibrium. We found only one binding pocket, namely antibody AB-C135 to remain closed along the entire minimum free energy path, suggesting an induced fit mechanism for this enzyme.

Molecular Simulations of Hydrophobic Gating of Pentameric Ligand Gated Ion Channels: Insights into Water and Ions

Ion channels are proteins which form gated nanopores in biological membranes. Many channels exhibit hydrophobic gating, whereby functional closure of a pore occurs by local dewetting. The pentameric ligand gated ion channels (pLGICs) provide a biologically important example of hydrophobic gating. Molecular simulation studies comparing additive vs polarizable models indicate predictions of hydrophobic gating are robust to the model employed. However, polarizable models suggest favorable interactions of hydrophobic pore-lining regions with chloride ions, of relevance to both synthetic carriers and channel proteins. Electrowetting of a closed pLGIC hydrophobic gate requires too high a voltage to occur physiologically but may inform designs for switchable nanopores. Global analysis of ∼200 channels yields a simple heuristic for structure-based prediction of (closed) hydrophobic gates. Simulation-based analysis is shown to provide an aid to interpretation of functional states of new channel structures. These studies indicate the importance of understanding the behavior of water and ions within the nanoconfined environment presented by ion channels.

Water in Nanopores and Biological Channels: A Molecular Simulation Perspective

This Review explores the dynamic behavior of water within nanopores and biological channels in lipid bilayer membranes. We focus on molecular simulation studies, alongside selected structural and other experimental investigations. Structures of biological nanopores and channels are reviewed, emphasizing those high-resolution crystal structures, which reveal water molecules within the transmembrane pores, which can be used to aid the interpretation of simulation studies. Different levels of molecular simulations of water within nanopores are described, with a focus on molecular dynamics (MD). In particular, models of water for MD simulations are discussed in detail to provide an evaluation of their use in simulations of water in nanopores. Simulation studies of the behavior of water in idealized models of nanopores have revealed aspects of the organization and dynamics of nanoconfined water, including wetting/dewetting in narrow hydrophobic nanopores. A survey of simulation studies in a range of nonbiological nanopores is presented, including carbon nanotubes, synthetic nanopores, model peptide nanopores, track-etched nanopores in polymer membranes, and hydroxylated and functionalized nanoporous silica. These reveal a complex relationship between pore size/geometry, the nature of the pore lining, and rates of water transport. Wider nanopores with hydrophobic linings favor water flow whereas narrower hydrophobic pores may show dewetting. Simulation studies over the past decade of the behavior of water in a range of biological nanopores are described, including porins and β-barrel protein nanopores, aquaporins and related polar solute pores, and a number of different classes of ion channels. Water is shown to play a key role in proton transport in biological channels and in hydrophobic gating of ion channels. An overall picture emerges, whereby the behavior of water in a nanopore may be predicted as a function of its hydrophobicity and radius. This informs our understanding of the functions of diverse channel structures and will aid the design of novel nanopores. Thus, our current level of understanding allows for the design of a nanopore which promotes wetting over dewetting or vice versa. However, to design a novel nanopore, which enables fast, selective, and gated flow of water de novo would remain challenging, suggesting a need for further detailed simulations alongside experimental evaluation of more complex nanopore systems.

Retrieving functional pathways of biomolecules from single-particle snapshots

A primary reason for the intense interest in structural biology is the fact that knowledge of structure can elucidate macromolecular functions in living organisms. Sustained effort has resulted in an impressive arsenal of tools for determining the static structures. But under physiological conditions, macromolecules undergo continuous conformational changes, a subset of which are functionally important. Techniques for capturing the continuous conformational changes underlying function are essential for further progress. Here, we present chemically-detailed conformational movies of biological function, extracted data-analytically from experimental single-particle cryo-electron microscopy (cryo-EM) snapshots of ryanodine receptor type 1 (RyR1), a calcium-activated calcium channel engaged in the binding of ligands. The functional motions differ substantially from those inferred from static structures in the nature of conformationally active structural domains, the sequence and extent of conformational motions, and the way allosteric signals are transduced within and between domains. Our approach highlights the importance of combining experiment, advanced data analysis, and molecular simulations.

Selectivity filter modalities and rapid inactivation of the hERG1 channel

The human ether-á-go-go-related gene (hERG1) channel conducts small outward K⁺ currents that are critical for cardiomyocyte membrane repolarization. The gain-of-function mutation N629D at the outer mouth of the selectivity filter (SF) disrupts inactivation and K⁺-selective transport in hERG1, leading to arrhythmogenic phenotypes associated with long-QT syndrome. Here, we combined computational electrophysiology with Markov state model analysis to investigate how SF-level gating modalities control selective cation transport in wild-type (WT) and mutant (N629D) hERG1 variants. Starting from the recently reported cryogenic electron microscopy (cryo-EM) open-state channel structure, multiple microseconds-long molecular-dynamics (MD) trajectories were generated using different cation configurations at the filter, voltages, electrolyte concentrations, and force-field parameters. Most of the K⁺ permeation events observed in hERG1-WT simulations occurred at microsecond timescales, influenced by the spontaneous dehydration/rehydration dynamics at the filter. The SF region displayed conductive, constricted, occluded, and dilated states, in qualitative agreement with the well-documented flickering conductance of hERG1. In line with mutagenesis studies, these gating modalities resulted from dynamic interaction networks involving residues from the SF, outer-mouth vestibule, P-helices, and S5-P segments. We found that N629D mutation significantly stabilizes the SF in a state that is permeable to both K⁺ and Na⁺, which is reminiscent of the SF in the nonselective bacterial NaK channel. Increasing the external K⁺ concentration induced "WT-like" SF dynamics in N629D, in qualitative agreement with the recovery of flickering currents in experiments. Overall, our findings provide an understanding of the molecular mechanisms controlling selective transport in K⁺ channels with a nonconventional SF sequence.

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.

The permeation mechanism of organic cations through a CNG mimic channel

Several channels, ranging from TRP receptors to Gap junctions, allow the exchange of small organic solute across cell membrane. However, very little is known about the molecular mechanism of their permeation. Cyclic Nucleotide Gated (CNG) channels, despite their homology with K+ channels and in contrast with them, allow the passage of larger methylated and ethylated ammonium ions like dimethylammonium (DMA) and ethylammonium (EA). We combined electrophysiology and molecular dynamics simulations to examine how DMA interacts with the pore and permeates through it. Due to the presence of hydrophobic groups, DMA enters easily in the channel and, unlike the alkali cations, does not need to cross any barrier. We also show that while the crystal structure is consistent with the presence of a single DMA ion at full occupancy, the channel is able to conduct a sizable current of DMA ions only when two ions are present inside the channel. Moreover, the second DMA ion dramatically changes the free energy landscape, destabilizing the crystallographic binding site and lowering by almost 25 kJ/mol the binding affinity between DMA and the channel. Based on the results of the simulation the experimental electron density maps can be re-interpreted with the presence of a second ion at lower occupancy. In this mechanism the flexibility of the channel plays a key role, extending the classical multi-ion permeation paradigm in which conductance is enhanced by the plain interaction between the ions.

Ion channels and ion selectivity

Specific macromolecular transport systems, ion channels and pumps, provide the pathways to facilitate and control the passage of ions across the lipid membrane. Ion channels provide energetically favourable passage for ions to diffuse rapidly and passively according to their electrochemical potential. Selective ion channels are essential for the excitability of biological membranes: the action potential is a transient phenomenon that reflects the rapid opening and closing of voltage-dependent Na⁺-selective and K⁺-selective channels. One of the most critical functional aspects of K⁺ channels is their ability to remain highly selective for K⁺ over Na⁺ while allowing high-throughput ion conduction at a rate close to the diffusion limit. Permeation through the K⁺ channel selectivity filter is believed to proceed as a 'knockon' mechanism, in which 2-3 K⁺ ions interspersed by water molecules move in a single file. Permeation through the comparatively wider and less selective Na⁺ channels also proceeds via a loosely coupled knockon mechanism, although the ions do not need to be fully dehydrated. While simple structural concepts are often invoked to rationalize the mechanism of ion selectivity, a deeper analysis shows that subtle effects play an important role in these flexible dynamical structures.

Atomistic insight into lipid translocation by a TMEM16 scramblase

The transmembrane protein 16 (TMEM16) family of membrane proteins includes both lipid scramblases and ion channels involved in olfaction, nociception, and blood coagulation. The crystal structure of the fungal Nectria haematococca TMEM16 (nhTMEM16) scramblase suggested a putative mechanism of lipid transport, whereby polar and charged lipid headgroups move through the low-dielectric environment of the membrane by traversing a hydrophilic groove on the membrane-spanning surface of the protein. Here, we use computational methods to explore the membrane-protein interactions involved in lipid scrambling. Fast, continuum membrane-bending calculations reveal a global pattern of charged and hydrophobic surface residues that bends the membrane in a large-amplitude sinusoidal wave, resulting in bilayer thinning across the hydrophilic groove. Atomic simulations uncover two lipid headgroup-interaction sites flanking the groove. The cytoplasmic site nucleates headgroup-dipole stacking interactions that form a chain of lipid molecules that penetrate into the groove. In two instances, a cytoplasmic lipid interdigitates into this chain, crosses the bilayer, and enters the extracellular leaflet, and the reverse process happens twice as well. Continuum membrane-bending analysis carried out on homology models of mammalian homologs shows that these family members also bend the membrane-even those that lack scramblase activity. Sequence alignments show that the lipid-interaction sites are conserved in many family members but less so in those with reduced scrambling ability. Our analysis provides insight into how large-scale membrane bending and protein chemistry facilitate lipid permeation in the TMEM16 family, and we hypothesize that membrane interactions also affect ion permeation.

Potassium channels in the heart: structure, function and regulation

This paper is the outcome of the fourth UC Davis Systems Approach to Understanding Cardiac Excitation-Contraction Coupling and Arrhythmias Symposium, a biannual event that aims to bring together leading experts in subfields of cardiovascular biomedicine to focus on topics of importance to the field. The theme of the 2016 symposium was 'K⁺ Channels and Regulation'. Experts in the field contributed their experimental and mathematical modelling perspectives and discussed emerging questions, controversies and challenges on the topic of cardiac K⁺ channels. This paper summarizes the topics of formal presentations and informal discussions from the symposium on the structural basis of voltage-gated K⁺ channel function, as well as the mechanisms involved in regulation of K⁺ channel gating, expression and membrane localization. Given the critical role for K⁺ channels in determining the rate of cardiac repolarization, it is hardly surprising that essentially every aspect of K⁺ channel function is exquisitely regulated in cardiac myocytes. This regulation is complex and highly interrelated to other aspects of myocardial function. K⁺ channel regulatory mechanisms alter, and are altered by, physiological challenges, pathophysiological conditions, and pharmacological agents. An accompanying paper focuses on the integrative role of K⁺ channels in cardiac electrophysiology, i.e. how K⁺ currents shape the cardiac action potential, and how their dysfunction can lead to arrhythmias, and discusses K⁺ channel-based therapeutics. A fundamental understanding of K⁺ channel regulatory mechanisms and disease processes is fundamental to reveal new targets for human therapy.

Urea-Aromatic Stacking and Concerted Urea Transport: Conserved Mechanisms in Urea Transporters Revealed by Molecular Dynamics

Urea transporters are membrane proteins that selectively allow urea molecules to pass through. It is not clear how these transporters allow rapid conduction of urea, a polar molecule, in spite of the presence of a hydrophobic constriction lined by aromatic rings. The current study elucidates the mechanism that is responsible for this rapid conduction by performing free energy calculations on the transporter dvUT with a cumulative sampling time of about 1.3 μs. A parallel arrangement of aromatic rings in the pore enables stacking of urea with these rings, which, in turn, lowers the energy barrier for urea transport. Such interaction of the rings with urea is proposed to be a conserved mechanism across all urea-conducting proteins. The free energy landscape for the permeation of multiple urea molecules reveals an interplay between interurea interaction and the solvation state of the urea molecules. This is for the first time that multiple molecule permeation through any small molecule transporter has been modeled.