Calcium-gated potassium channel blockade via membrane-facing fenestrations.

Quaternary ammonium blockers were previously shown to bind in the pore to block both open and closed conformations of large-conductance calcium-activated potassium (BK and MthK) channels. Because blocker entry was assumed through the intracellular entryway (bundle crossing), closed-pore access suggested that the gate was not at the bundle crossing. Structures of closed MthK, a Methanobacterium thermoautotrophicum homolog of BK channels, revealed a tightly constricted intracellular gate, leading us to investigate the membrane-facing fenestrations as alternative pathways for blocker access directly from the membrane. Atomistic free energy simulations showed that intracellular blockers indeed access the pore through the fenestrations, and a mutant channel with narrower fenestrations displayed no closed-state TPeA block at concentrations that blocked the wild-type channel. Apo BK channels display similar fenestrations, suggesting that blockers may use them as access paths into closed channels. Thus, membrane fenestrations represent a non-canonical pathway for selective targeting of specific channel conformations, opening novel ways to selectively drug BK channels.

Ball-and-chain inactivation in a calcium-gated potassium channel.

Inactivation is the process by which ion channels terminate ion flux through their pores while the opening stimulus is still present1. In neurons, inactivation of both sodium and potassium channels is crucial for the generation of action potentials and regulation of firing frequency1,2. A cytoplasmic domain of either the channel or an accessory subunit is thought to plug the open pore to inactivate the channel via a 'ball-and-chain' mechanism3-7. Here we use cryo-electron microscopy to identify the molecular gating mechanism in calcium-activated potassium channels by obtaining structures of the MthK channel from Methanobacterium thermoautotrophicum-a purely calcium-gated and inactivating channel-in a lipid environment. In the absence of Ca2+, we obtained a single structure in a closed state, which was shown by atomistic simulations to be highly flexible in lipid bilayers at ambient temperature, with large rocking motions of the gating ring and bending of pore-lining helices. In Ca2+-bound conditions, we obtained several structures, including multiple open-inactivated conformations, further indication of a highly dynamic protein. These different channel conformations are distinguished by rocking of the gating rings with respect to the transmembrane region, indicating symmetry breakage across the channel. Furthermore, in all conformations displaying open channel pores, the N terminus of one subunit of the channel tetramer sticks into the pore and plugs it, with free energy simulations showing that this is a strong interaction. Deletion of this N terminus leads to functionally non-inactivating channels and structures of open states without a pore plug, indicating that this previously unresolved N-terminal peptide is responsible for a ball-and-chain inactivation mechanism.

Atomistic Simulations of Membrane Ion Channel Conduction, Gating, and Modulation.

Membrane ion channels are the fundamental electrical components in the nervous system. Recent developments in X-ray crystallography and cryo-EM microscopy have revealed what these proteins look like in atomic detail but do not tell us how they function. Molecular dynamics simulations have progressed to the point that we can now simulate realistic molecular assemblies to produce quantitative calculations of the thermodynamic and kinetic quantities that control function. In this review, we summarize the state of atomistic simulation methods for ion channels to understand their conduction, activation, and drug modulation mechanisms. We are at a crossroads in atomistic simulation, where long time scale observation can provide unbiased exploration of mechanisms, supplemented by biased free energy methodologies. We illustrate the use of these approaches to describe ion conduction and selectivity in voltage-gated sodium and acid-sensing ion channels. Studies of channel gating present a significant challenge, as activation occurs on longer time scales. Enhanced sampling approaches can ensure convergence on minimum free energy pathways for activation, as illustrated here for pentameric ligand-gated ion channels that are principal to nervous system function and the actions of general anesthetics. We also examine recent studies of local anesthetic and antiepileptic drug binding to a sodium channel, revealing sites and pathways that may offer new targets for drug development. Modern simulations thus offer a range of molecular-level insights into ion channel function and modulation as a learning platform for mechanistic discovery and drug development.

Selective ion permeation involves complexation with carboxylates and lysine in a model human sodium channel.

Bacterial and human voltage-gated sodium channels (Navs) exhibit similar cation selectivity, despite their distinct EEEE and DEKA selectivity filter signature sequences. Recent high-resolution structures for bacterial Navs have allowed us to learn about ion conduction mechanisms in these simpler homo-tetrameric channels, but our understanding of the function of their mammalian counterparts remains limited. To probe these conduction mechanisms, a model of the human Nav1.2 channel has been constructed by grafting residues of its selectivity filter and external vestibular region onto the bacterial NavRh channel with atomic-resolution structure. Multi-µs fully atomistic simulations capture long time-scale ion and protein movements associated with the permeation of Na+ and K+ ions, and their differences. We observe a Na+ ion knock-on conduction mechanism facilitated by low energy multi-carboxylate/multi-Na+ complexes, akin to the bacterial channels. These complexes involve both the DEKA and vestibular EEDD rings, acting to draw multiple Na+ into the selectivity filter and promote permeation. When the DEKA ring lysine is protonated, we observe that its ammonium group is actively participating in Na+ permeation, presuming the role of another ion. It participates in the formation of a stable complex involving carboxylates that collectively bind both Na+ and the Lys ammonium group in a high-field strength site, permitting pass-by translocation of Na+. In contrast, multiple K+ ion complexes with the DEKA and EEDD rings are disfavored by up to 8.3 kcal/mol, with the K+-lysine-carboxylate complex non-existent. As a result, lysine acts as an electrostatic plug that partially blocks the flow of K+ ions, which must instead wait for isomerization of lysine downward to clear the path for K+ passage. These distinct mechanisms give us insight into the nature of ion conduction and selectivity in human Nav channels, while uncovering high field strength carboxylate binding complexes that define the more general phenomenon of Na+-selective ion transport in nature.

Potassium dependent structural changes in the selectivity filter of HERG potassium channels.

The fine tuning of biological electrical signaling is mediated by variations in the rates of opening and closing of gates that control ion flux through different ion channels. Human ether-a-go-go related gene (HERG) potassium channels have uniquely rapid inactivation kinetics which are critical to the role they play in regulating cardiac electrical activity. Here, we exploit the K+ sensitivity of HERG inactivation to determine structures of both a conductive and non-conductive selectivity filter structure of HERG. The conductive state has a canonical cylindrical shaped selectivity filter. The non-conductive state is characterized by flipping of the selectivity filter valine backbone carbonyls to point away from the central axis. The side chain of S620 on the pore helix plays a central role in this process, by coordinating distinct sets of interactions in the conductive, non-conductive, and transition states. Our model represents a distinct mechanism by which ion channels fine tune their activity and could explain the uniquely rapid inactivation kinetics of HERG.

Determinants of ion selectivity in ASIC1a- and ASIC2a-containing acid-sensing ion channels.

Trimeric acid-sensing ion channels (ASICs) contribute to neuronal signaling by converting extracellular acidification into excitatory sodium currents. Previous work with homomeric ASIC1a implicates conserved leucine (L7') and consecutive glycine-alanine-serine (GAS belt) residues near the middle, and conserved negatively charged (E18') residues at the bottom of the pore in ion permeation and/or selectivity. However, a conserved mechanism of ion selectivity throughout the ASIC family has not been established. We therefore explored the molecular determinants of ion selectivity in heteromeric ASIC1a/ASIC2a and homomeric ASIC2a channels using site-directed mutagenesis, electrophysiology, and molecular dynamics free energy simulations. Similar to ASIC1a, E18' residues create an energetic preference for sodium ions at the lower end of the pore in ASIC2a-containing channels. However, and in contrast to ASIC1a homomers, ion permeation through ASIC2a-containing channels is not determined by L7' side chains in the upper part of the channel. This may be, in part, due to ASIC2a-specific negatively charged residues (E59 and E62) that lower the energy of ions in the upper pore, thus making the GAS belt more important for selectivity. This is confirmed by experiments showing that the L7'A mutation has no effect in ASIC2a, in contrast to ASIC1a, where it eliminated selectivity. ASIC2a triple mutants eliminating both L7' and upper charges did not lead to large changes in selectivity, suggesting a different role for L7' in ASIC2a compared with ASIC1a channels. In contrast, we observed measurable changes in ion selectivity in ASIC2a-containing channels with GAS belt mutations. Our results suggest that ion conduction and selectivity in the upper part of the ASIC pore may differ between subtypes, whereas the essential role of E18' in ion selectivity is conserved. Furthermore, we demonstrate that heteromeric channels containing mutations in only one of two ASIC subtypes provide a means of functionally testing mutations that render homomeric channels nonfunctional.

Comparison of permeation mechanisms in sodium-selective ion channels.

Voltage-gated sodium channels are the molecular components of electrical signaling in the body, yet the molecular origins of Na+-selective transport remain obscured by diverse protein chemistries within this family of ion channels. In particular, bacterial and mammalian sodium channels are known to exhibit similar relative ion permeabilities for Na+ over K+ ions, despite their distinct signature EEEE and DEKA sequences. Atomic-level molecular dynamics simulations using high-resolution bacterial channel structures and mammalian channel models have begun to describe how these sequences lead to analogous high field strength ion binding sites that drive Na+ conduction. Similar complexes have also been identified in unrelated acid sensing ion channels involving glutamate and aspartate side chains that control their selectivity. These studies suggest the possibility of a common origin for Na+ selective binding and transport.

A model of the effects of cancer cell motility and cellular adhesion properties on tumour-immune dynamics.

We present a three-dimensional model simulating the dynamics of an anti-cancer T-cell response against a small, avascular, early-stage tumour. Interactions at the tumour site are accounted for using an agent-based model (ABM), while immune cell dynamics in the lymph node are modelled as a system of delay differential equations (DDEs). We combine these separate approaches into a two-compartment hybrid ABM-DDE system to capture the T-cell response against the tumour. In the ABM at the tumour site, movement of tumour cells is modelled using effective physical forces with a specific focus on cell-to-cell adhesion properties and varying levels of tumour cell motility, thus taking into account the ability of cancer cells to spread and form clusters. We consider the effectiveness of the immune response over a range of parameters pertaining to tumour cell motility, cell-to-cell adhesion strength and growth rate. We also investigate the dependence of outcomes on the distribution of tumour cells. Low tumour cell motility is generally a good indicator for successful tumour eradication before relapse, while high motility leads, almost invariably, to relapse and tumour escape. In general, the effect of cell-to-cell adhesion on prognosis is dependent on the level of tumour cell motility, with an often unpredictable cross influence between adhesion and motility, which can lead to counterintuitive effects. In terms of overall tumour shape and structure, the spatial distribution of cancer cells in clusters of various sizes has shown to be strongly related to the likelihood of extinction.

A selectivity filter at the intracellular end of the acid-sensing ion channel pore.

Increased extracellular proton concentrations during neurotransmission are converted to excitatory sodium influx by acid-sensing ion channels (ASICs). 10-fold sodium/potassium selectivity in ASICs has long been attributed to a central constriction in the channel pore, but experimental verification is lacking due to the sensitivity of this structure to conventional manipulations. Here, we explored the basis for ion selectivity by incorporating unnatural amino acids into the channel, engineering channel stoichiometry and performing free energy simulations. We observed no preference for sodium at the "GAS belt" in the central constriction. Instead, we identified a band of glutamate and aspartate side chains at the lower end of the pore that enables preferential sodium conduction.