Structure of a bacterial voltage-gated sodium channel pore reveals mechanisms of opening and closing

The chemistry of electrical signaling in sodium channels from bacteria and beyond

Electrical signaling is essential for all fast processes in biology, but its molecular mechanisms have been uncertain. This review article focuses on studies of bacterial sodium channels in order to home in on the essential molecular and chemical mechanisms underlying transmembrane ion conductance and voltage-dependent gating without the overlay of complex protein interactions and regulatory mechanisms in mammalian sodium channels. This minimalist approach has yielded a nearly complete picture of sodium channel function at the atomic level that are mostly conserved in mammalian sodium channels, including sodium selectivity and conductance, voltage sensing and activation, electromechanical coupling to pore opening and closing, slow inactivation, and pathogenic dysfunction in a debilitating channelopathy. Future studies of nature's simplest sodium channels may continue to yield key insights into the fundamental molecular and chemical principles of their function and further elucidate the chemical basis of electrical signaling.

Functional effects of drugs and toxins interacting with NaV1.4

NaV1.4 is a voltage-gated sodium channel subtype that is predominantly expressed in skeletal muscle cells. It is essential for producing action potentials and stimulating muscle contraction, and mutations in NaV1.4 can cause various muscle disorders. The discovery of the cryo-EM structure of NaV1.4 in complex with β1 has opened new possibilities for designing drugs and toxins that target NaV1.4. In this review, we summarize the current understanding of channelopathies, the binding sites and functions of chemicals including medicine and toxins that interact with NaV1.4. These substances could be considered novel candidate compounds or tools to develop more potent and selective drugs targeting NaV1.4. Therefore, studying NaV1.4 pharmacology is both theoretically and practically meaningful.

Water, Protons, and the Gating of Voltage-Gated Potassium Channels

Ion channels are ubiquitous throughout all forms of life. Potassium channels are even found in viruses. Every cell must communicate with its surroundings, so all cells have them, and excitable cells, in particular, especially nerve cells, depend on the behavior of these channels. Every channel must be open at the appropriate time, and only then, so that each channel opens in response to the stimulus that tells that channel to open. One set of channels, including those in nerve cells, responds to voltage. There is a standard model for the gating of these channels that has a section of the protein moving in response to the voltage. However, there is evidence that protons are moving, rather than protein. Water is critical as part of the gating process, although it is hard to see how this works in the standard model. Here, we review the extensive evidence of the importance of the role of water and protons in gating these channels. Our principal example, but by no means the only example, will be the Kv1.2 channel. Evidence comes from the effects of D2O, from mutations in the voltage sensing domain, as well as in the linker between that domain and the gate, and at the gate itself. There is additional evidence from computations, especially quantum calculations. Structural evidence comes from X-ray studies. The hydration of ions is critical in the transfer of ions in constricted spaces, such as the gate region and the pore of a channel; we will see how the structure of the hydrated ion fits with the structure of the channel. In addition, there is macroscopic evidence from osmotic experiments and streaming current measurements. The combined evidence is discussed in the context of a model that emphasizes the role of protons and water in gating these channels.

Microbial membrane transport proteins and their biotechnological applications

Because of the hydrophobic nature of the membrane lipid bilayer, the majority of the hydrophilic solutes require special transportation mechanisms for passing through the cell membrane. Integral membrane transport proteins (MTPs), which belong to the Major Intrinsic Protein Family, facilitate the transport of these solutes across cell membranes. MTPs including aquaporins and carrier proteins are transmembrane proteins spanning across the cell membrane. The easy handling of microorganisms enabled the discovery of a remarkable number of transport proteins specific to different substances. It has been realized that these transporters have very important roles in the survival of microorganisms, their pathogenesis, and antimicrobial resistance. Astonishing features related to the solute specificity of these proteins have led to the acceleration of the research on the discovery of their properties and the development of innovative products in which these unique properties are used or imitated. Studies on microbial MTPs range from the discovery and characterization of a novel transporter protein to the mining and screening of them in a large transporter library for particular functions, from simulations and modeling of specific transporters to the preparation of biomimetic synthetic materials for different purposes such as biosensors or filtration membranes. This review presents recent discoveries on microbial membrane transport proteins and focuses especially on formate nitrite transport proteins and aquaporins, and advances in their biotechnological applications.

Anticonvulsant Classes and Possible Mechanism of Actions

Epilepsy is considered one of the most common neurological disorders worldwide; it needs long-term or life-long treatment. Despite the presence of several novel antiepileptic drugs, approximately 30% patients still suffer from drug-resistant epilepsy. Subsequently, searching for new anticonvulsants with lower toxicity and better efficacy is still in paramount demand. Using target-based studies in the discovery of novel antiepileptics is uncommon owing to the insufficient information on the molecular pathway of epilepsy and complex mode of action for most of known antiepileptic drugs. In this review, we investigated the properties of anticonvulsants, types of epileptic seizures, and mechanism of action for anticonvulsants.

Binding of carboxylate and water to monovalent cations

The interactions of carboxylate anions with water and cations are important for a wide variety of systems, both biological and synthetic. To gain insight on properties of the local complexes, we apply density functional theory, to treat the complex electrostatic interactions, and investigate mixtures with varied numbers of carboxylate anions (acetate) and waters binding to monovalent cations, Li+, Na+ and K+. The optimal structure with overall lowest free energy contains two acetates and two waters such that the cation is four-fold coordinated, similar to structures found earlier for pure water or pure carboxylate ligands. More generally, the complexes with two acetates have the lowest free energy. In transitioning from the overall optimal state, exchanging an acetate for water has a lower free energy barrier than exchanging water for an acetate. In most cases, the carboxylates are monodentate and in the first solvation shell. As water is added to the system, hydrogen bonding between waters and carboxylate O atoms further stabilizes monodentate structures. These structures, which have strong electrostatic interactions that involve hydrogen bonds of varying strength, are significantly polarized, with ChelpG partial charges that vary substantially as the bonding geometry varies. Overall, these results emphasize the increasing importance of water as a component of binding sites as the number of ligands increases, thus affecting the preferential solvation of specific metal ions and clarifying Hofmeister effects. Finally, structural analysis correlated with free energy analysis supports the idea that binding to more than the preferred number of carboxylates under architectural constraints are a key to ion transport.

Determining the structure of the bacterial voltage-gated sodium channel NaChBac embedded in liposomes by cryo electron tomography and subtomogram averaging

Voltage-gated sodium channels shape action potentials that propagate signals along cells. When the membrane potential reaches a certain threshold, the channels open and allow sodium ions to flow through the membrane depolarizing it, followed by the deactivation of the channels. Opening and closing of the channels is important for cellular signalling and regulates various physiological processes in muscles, heart and brain. Mechanistic insights into the voltage-gated channels are difficult to achieve as the proteins are typically extracted from membranes for structural analysis which results in the loss of the transmembrane potential that regulates their activity. Here, we report the structural analysis of a bacterial voltage-gated sodium channel, NaChBac, reconstituted in liposomes under an electrochemical gradient by cryo electron tomography and subtomogram averaging. We show that the small channel, most of the residues of which are embedded in the membrane, can be localized using a genetically fused GFP. GFP can aid the initial alignment to an average resulting in a correct structure, but does not help for the final refinement. At a moderate resolution of ˜16 Å the structure of NaChBac in an unrestricted membrane bilayer is 10% wider than the structure of the purified protein previously solved in nanodiscs, suggesting the potential movement of the peripheral voltage-sensing domains. Our study explores the limits of structural analysis of membrane proteins in membranes.

Ion channel selectivity through ion-modulated changes of selectivity filter pKa values

In bacterial voltage-gated sodium channels, the passage of ions through the pore is controlled by a selectivity filter (SF) composed of four glutamate residues. The mechanism of selectivity has been the subject of intense research, with suggested mechanisms based on steric effects, and ion-triggered conformational change. Here, we propose an alternative mechanism based on ion-triggered shifts in pKa values of SF glutamates. We study the NavMs channel for which the open channel structure is available. Our free-energy calculations based on molecular dynamics simulations suggest that pKa values of the four glutamates are higher in solution of K+ ions than in solution of Na+ ions. Higher pKa in the presence of K+ stems primarily from the higher population of dunked conformations of the protonated Glu sidechain, which exhibit a higher pKa shift. Since pKa values are close to the physiological pH, this results in predominant population of the fully deprotonated state of glutamates in Na+ solution, while protonated states are predominantly populated in K+ solution. Through molecular dynamics simulations we calculate that the deprotonated state is the most conductive, the singly protonated state is less conductive, and the doubly protonated state has significantly reduced conductance. Thus, we propose that a significant component of selectivity is achieved through ion-triggered shifts in the protonation state, which favors more conductive states for Na+ ions and less conductive states for K+ ions. This mechanism also suggests a strong pH dependence of selectivity, which has been experimentally observed in structurally similar NaChBac channels.

New aryl and acylsulfonamides as state-dependent inhibitors of Nav1.3 voltage-gated sodium channel

Voltage-gated sodium channels (Na_vs) play an essential role in neurotransmission, and their dysfunction is often a cause of various neurological disorders. The Na_v1.3 isoform is found in the CNS and upregulated after injury in the periphery, but its role in human physiology has not yet been fully elucidated. Reports suggest that selective Na_v1.3 inhibitors could be used as novel therapeutics to treat pain or neurodevelopmental disorders. Few selective inhibitors of this channel are known in the literature. In this work, we report the discovery of a new series of aryl and acylsulfonamides as state-dependent inhibitors of Na_v1.3 channels. Using a ligand-based 3D similarity search and subsequent hit optimization, we identified and prepared a series of 47 novel compounds and tested them on Na_v1.3, Na_v1.5, and a selected subset also on Na_v1.7 channels in a QPatch patch-clamp electrophysiology assay. Eight compounds had an IC₅₀ value of less than 1 μM against the Na_v1.3 channel inactivated state, with one compound displaying an IC₅₀ value of 20 nM, whereas activity against the inactivated state of the Na_v1.5 channel and Na_v1.7 channel was approximately 20-fold weaker. None of the compounds showed use-dependent inhibition of the cardiac isoform Na_v1.5 at a concentration of 30 μM. Further selectivity testing of the most promising hits was measured using the two-electrode voltage-clamp method against the closed state of the Na_v1.1-Na_v1.8 channels, and compound 15b displayed small, yet selective, effects against the Na_v1.3 channel, with no activity against the other isoforms. Additional selectivity testing of promising hits against the inactivated state of the Na_v1.3, Na_v1.7, and Na_v1.8 channels revealed several compounds with robust and selective activity against the inactivated state of the Na_v1.3 channel among the three isoforms tested. Moreover, the compounds were not cytotoxic at a concentration of 50 μM, as demonstrated by the assay in human HepG2 cells (hepatocellular carcinoma cells). The novel state-dependent inhibitors of Na_v1.3 discovered in this work provide a valuable tool to better evaluate this channel as a potential drug target.

Molecular Dynamics Simulations of Ion Permeation in Human Voltage-Gated Sodium Channels

The recent determination of cryo-EM structures of voltage-gated sodium (Na_v) channels has revealed many details of these proteins. However, knowledge of ionic permeation through the Na_v pore remains limited. In this work, we performed atomistic molecular dynamics (MD) simulations to study the structural features of various neuronal Na_v channels based on homology modeling of the cryo-EM structure of the human Na_v1.4 channel and, in addition, on the recently resolved configuration for Na_v1.2. In particular, single Na⁺ permeation events during standard MD runs suggest that the ion resides in the inner part of the Na_v selectivity filter (SF). On-the-fly free energy parametrization (OTFP) temperature-accelerated molecular dynamics (TAMD) was also used to calculate two-dimensional free energy surfaces (FESs) related to single/double Na⁺ translocation through the SF of the homology-based Na_v1.2 model and the cryo-EM Na_v1.2 structure, with different realizations of the DEKA filter domain. These additional simulations revealed distinct mechanisms for single and double Na⁺ permeation through the wild-type SF, which has a charged lysine in the DEKA ring. Moreover, the configurations of the ions in the SF corresponding to the metastable states of the FESs are specific for each SF motif. Overall, the description of these mechanisms gives us new insights into ion conduction in human Na_v cryo-EM-based and cryo-EM configurations that could advance understanding of these systems and how they differ from potassium and bacterial Na_v channels.

Inhibition of voltage-gated sodium ion channel by corannulene and computational inversion blockage underlying mechanisms

Corannulene (Cor), a special carbon material, evidenced strong protein binding capacity which regulating lysozyme crystallization and controlling reactive oxygen species (ROS) generation. Ion channel protein play role in regulating ion channel functions to affect physiological functions. However, the interaction between Cor and ion channel protein have not been studied. In this study, PEG/Cor nanoparticles (PEG/Cor Nps) were prepared by mPEG-DSPE. The PEG/Cor Nps localized in cytoplasm and produced cytotoxicity at high concentration. Whole cell patch clamp examined ion channel functions after incubate PEG/Cor Nps with PC-12 cell. we found that PEG/Cor Nps inhibited voltage-gated Na⁺ ion channels in a dose- and time-dependent manner but not act on voltage-gated K⁺ ion channels. The potential mechanisms were revealed by all-atom molecular dynamic (MD) simulations. The results showed that PEG/Cor Nps block the pore of sodium ion channel protein due to dose- and time-dependent accumulation. Besides, the orientation angle (θ) configuration of PEG/Cor Nps will be inverted with the accumulation to generate two blocking mechanisms. Different from other carbon nanomaterials, the blockage mechanism of PEG/Cor Nps provides novel insights into the mechanisms of interaction between carbon nanomaterials and protein.

An α-π transition in S6 shapes the conformational cycle of the bacterial sodium channel NavAb

Voltage-gated sodium channels play an important role in electrical signaling in excitable cells. In response to changes in membrane potential, they cycle between nonconducting and conducting conformations. With recent advances in structural biology, structures of sodium channels have been captured in several distinct conformations, which are thought to represent different functional states. However, it has been difficult to capture the intrinsically transient open state. We recently showed that a proposed open state of the bacterial sodium channel NavMs was not conductive and that a conformational change involving a transition to a π-helix in the pore-lining S6 helix converted this structure into a conducting state. However, the relevance of this structural feature in other sodium channels, and its implications for the broader gating cycle, remained unclear. Here, we propose a comparable open state of another class of bacterial channel from Aliarcobacter butzleri (NavAb) with characteristic pore hydration, ion permeation, and drug binding properties. Furthermore, we show that a π-helix transition can lead to pore opening and that such a conformational change blocks fenestrations in the inner helix bundle. We also discover that a region in the C-terminal domain can undergo a disordering transition proposed to be important for pore opening. These results support a role for a π-helix transition in the opening of NavAb, enabling new proposals for the structural annotation and drug modulation mechanisms in this important sodium channel model.

Simulation and Machine Learning Methods for Ion-Channel Structure Determination, Mechanistic Studies and Drug Design

Ion channels are expressed in almost all living cells, controlling the in-and-out communications, making them ideal drug targets, especially for central nervous system diseases. However, owing to their dynamic nature and the presence of a membrane environment, ion channels remain difficult targets for the past decades. Recent advancement in cryo-electron microscopy and computational methods has shed light on this issue. An explosion in high-resolution ion channel structures paved way for structure-based rational drug design and the state-of-the-art simulation and machine learning techniques dramatically improved the efficiency and effectiveness of computer-aided drug design. Here we present an overview of how simulation and machine learning-based methods fundamentally changed the ion channel-related drug design at different levels, as well as the emerging trends in the field.

Supramolecular Mechanosensitive Potassium Channel Formed by Fluorinated Amphiphilic Cyclophane

Inspired by mechanosensitive potassium channels found in nature, we developed a fluorinated amphiphilic cyclophane composed of fluorinated rigid aromatic units connected via flexible hydrophilic octa(ethylene glycol) chains. Microscopic and emission spectroscopic studies revealed that the cyclophane could be incorporated into the hydrophobic layer of the lipid bilayer membranes and self-assembled to form a supramolecular transmembrane ion channel. Current recording measurements using cyclophane-containing planer lipid bilayer membranes successfully demonstrated an efficient transmembrane ion transport. We also demonstrated that the ion transport property was sensitive to the mechanical forces applied to the membranes. In addition, ion transport assays using pH-sensitive fluorescence dye revealed that the supramolecular channel possesses potassium ion selectivity. We also performed all-atom hybrid quantum-mechanical/molecular mechanical simulations to assess the channel structures at atomic resolution and the mechanism of selective potassium ion transport. This research demonstrated the first example of a synthetic mechanosensitive potassium channel, which would open a new door to sensing and manipulating biologically important processes and purification of key materials in industries.

Quaternary structure independent folding of voltage-gated ion channel pore domain subunits

Every voltage-gated ion channel (VGIC) has a pore domain (PD) made from four subunits, each comprising an antiparallel transmembrane helix pair bridged by a loop. The extent to which PD subunit structure requires quaternary interactions is unclear. Here, we present crystal structures of a set of bacterial voltage-gated sodium channel (BacNa_V) 'pore only' proteins that reveal a surprising collection of non-canonical quaternary arrangements in which the PD tertiary structure is maintained. This context-independent structural robustness, supported by molecular dynamics simulations, indicates that VGIC-PD tertiary structure is independent of quaternary interactions. This fold occurs throughout the VGIC superfamily and in diverse transmembrane and soluble proteins. Strikingly, characterization of PD subunit-binding Fabs indicates that non-canonical quaternary PD conformations can occur in full-length VGICs. Together, our data demonstrate that the VGIC-PD is an autonomously folded unit. This property has implications for VGIC biogenesis, understanding functional states, de novo channel design, and VGIC structural origins.

Different acyl-CoA:diacylglycerol acyltransferases vary widely in function, and a targeted amino acid substitution enhances oil accumulation

Triacylglycerols (TAGs) are the major component of plant storage lipids such as oils. Acyl-CoA:diacylglycerol acyltransferase (DGAT) catalyzes the final step of the Kennedy pathway, and is mainly responsible for plant oil accumulation. We previously found that the activity of Vernonia DGAT1 was distinctively higher than that of Arabidopsis and soybean DGAT1 in a yeast microsome assay. In this study, the DGAT1 cDNAs of Arabidopsis, Vernonia, soybean, and castor bean were introduced into Arabidopsis. All Vernonia DGAT1-expressing lines showed a significantly higher oil content (49% mean increase compared with the wild-type) followed by soybean and castor bean. Most Arabidopsis DGAT1-overexpressing lines did not show a significant increase. In addition to these four DGAT1 genes, sunflower, Jatropha, and sesame DGAT1 genes were introduced into a TAG biosynthesis-defective yeast mutant. In the yeast expression culture, DGAT1s from Arabidopsis, castor bean, and soybean only slightly increased the TAG content; however, DGAT1s from Vernonia, sunflower, Jatropha, and sesame increased TAG content >10-fold more than the former three DGAT1s. Three amino acid residues were characteristically common in the latter four DGAT1s. Using soybean DGAT1, these amino acid substitutions were created by site-directed mutagenesis and substantially increased the TAG content.

Druggability of Voltage-Gated Sodium Channels-Exploring Old and New Drug Receptor Sites

Voltage-gated ion channels are important drug targets because they play crucial physiological roles in both excitable and non-excitable cells. About 15% of clinical drugs used for treating human diseases target ion channels. However, most of these drugs do not provide sufficient specificity to a single subtype of the channels and their off-target side effects can be serious and sometimes fatal. Recent advancements in imaging techniques have enabled us for the first time to visualize unique and hidden parts of voltage-gated sodium channels in different structural conformations, and to develop drugs that further target a selected functional state in each channel subtype with the potential for high precision and low toxicity. In this review we describe the druggability of voltage-gated sodium channels in distinct functional states, which could potentially be used to selectively target the channels. We review classical drug receptors in the channels that have recently been structurally characterized by cryo-electron microscopy with natural neurotoxins and clinical drugs. We further examine recent drug discoveries for voltage-gated sodium channels and discuss opportunities to use distinct, state-dependent receptor sites in the voltage sensors as unique drug targets. Finally, we explore potential new receptor sites that are currently unknown for sodium channels but may be valuable for future drug discovery. The advancement presented here will help pave the way for drug development that selectively targets voltage-gated sodium channels.

P-Loop Channels: Experimental Structures, and Physics-Based and Neural Networks-Based Models

The superfamily of P-loop channels includes potassium, sodium, and calcium channels, as well as TRP channels and ionotropic glutamate receptors. A rapidly increasing number of crystal and cryo-EM structures have revealed conserved and variable elements of the channel structures. Intriguing differences are seen in transmembrane helices of channels, which may include π-helical bulges. The bulges reorient residues in the helices and thus strongly affect their intersegment contacts and patterns of ligand-sensing residues. Comparison of the experimental structures suggests that some π-bulges are dynamic: they may appear and disappear upon channel gating and ligand binding. The AlphaFold2 models represent a recent breakthrough in the computational prediction of protein structures. We compared some crystal and cryo-EM structures of P-loop channels with respective AlphaFold2 models. Folding of the regions, which are resolved experimentally, is generally similar to that predicted in the AlphaFold2 models. The models also reproduce some subtle but significant differences between various P-loop channels. However, patterns of π-bulges do not necessarily coincide in the experimental and AlphaFold2 structures. Given the importance of dynamic π-bulges, further studies involving experimental and theoretical approaches are necessary to understand the cause of the discrepancy.

Fenestropathy of Voltage-Gated Sodium Channels

Voltage-gated sodium channels (Nav) are responsible for the initiation and propagation of action potentials in excitable cells. From pain to heartbeat, these integral membrane proteins are the ignition stations for every sensation and action in human bodies. They are large (>200 kDa, 24 transmembrane helices) multi-domain proteins that couple changes in membrane voltage to the gating cycle of the sodium-selective pore. Nav mutations lead to a multitude of diseases - including chronic pain, cardiac arrhythmia, muscle illnesses, and seizure disorders - and a wide variety of currently used therapeutics block Nav. Despite this, the mechanisms of action of Nav blocking drugs are only modestly understood at this time and many questions remain to be answered regarding their state- and voltage-dependence, as well as the role of the hydrophobic membrane access pathways, or fenestrations, in drug ingress or egress. Nav fenestrations, which are pathways that connect the plasma membrane to the central cavity in the pore domain, were discovered through functional studies more than 40 years ago and once thought to be simple pathways. A variety of recent genetic, structural, and pharmacological data, however, shows that these fenestrations are actually key functional regions of Nav that modulate drug binding, lipid binding, and influence gating behaviors. We discovered that some of the disease mutations that cause arrhythmias alter amino acid residues that line the fenestrations of Nav1.5. This indicates that fenestrations may play a critical role in channel's gating, and that individual genetic variation may also influence drug access through the fenestrations for resting/inactivated state block. In this review, we will discuss the channelopathies associated with these fenestrations, which we collectively name "Fenestropathy," and how changes in the fenestrations associated with the opening of the intracellular gate could modulate the state-dependent ingress and egress of drugs binding in the central cavity of voltage gated sodium channels.

Characterizing fenestration size in sodium channel subtypes and their accessibility to inhibitors

Voltage-gated sodium channels (Nav) underlie the electrical activity of nerve and muscle cells. Humans have nine different subtypes of these channels, which are the target of small-molecule inhibitors commonly used to treat a range of conditions. Structural studies have identified four lateral fenestrations within the Nav pore module that have been shown to influence Nav pore blocker access during resting-state inhibition. However, the structural differences among the nine subtypes are still unclear. In particular, the dimensions of the four individual fenestrations across the Nav subtypes and their differential accessibility to pore blockers is yet to be characterized. To address this, we applied classical molecular dynamics simulations to study the recently published structures of Nav1.1, Nav1.2, Nav1.4, Nav1.5, and Nav1.7. Although there is significant variability in the bottleneck sizes of the Nav fenestrations, the subtypes follow a common pattern, with wider DI-II and DIII-IV fenestrations, a more restricted DII-III fenestration, and the most restricted DI-IV fenestration. We further identify the key bottleneck residues in each fenestration and show that the motions of aromatic residue sidechains govern the bottleneck radii. Well-tempered metadynamics simulations of Nav1.4 and Nav1.5 in the presence of the pore blocker lidocaine also support the DI-II fenestration being the most likely access route for drugs. Our computational results provide a foundation for future in vitro experiments examining the route of drug access to sodium channels. Understanding the fenestrations and their accessibility to drugs is critical for future analyses of diseases mutations across different sodium channel subtypes, with the potential to inform pharmacological development of resting-state inhibitors and subtype-selective drug design.

The T1-tetramerisation domain of Kv1.2 rescues expression and preserves function of a truncated NaChBac sodium channel

Cytoplasmic domains frequently promote functional assembly of multimeric ion channels. To investigate structural determinants of this process, we generated the ‘T1‐chimera’ construct of the NaChBac sodium channel by truncating its C‐terminal domain and splicing the T1‐tetramerisation domain of the Kv1.2 channel to the N terminus. Purified T1‐chimera channels were tetrameric, conducted Na⁺ when reconstituted into proteoliposomes, and were functionally blocked by the drug mibefradil. Both the T1‐chimera and full‐length NaChBac had comparable expression levels in the membrane, whereas a NaChBac mutant lacking a cytoplasmic domain had greatly reduced membrane expression. Our findings support a model whereby bringing the transmembrane regions into close proximity enables their tetramerisation. This phenomenon is found with other channels, and thus, our findings substantiate this as a common assembly mechanism.

An open state of a voltage-gated sodium channel involving a π-helix and conserved pore-facing asparagine

Voltage-gated sodium (Nav) channels play critical roles in propagating action potentials and otherwise manipulating ionic gradients in excitable cells. These channels open in response to membrane depolarization, selectively permeating sodium ions until rapidly inactivating. Structural characterization of the gating cycle in this channel family has proved challenging, particularly due to the transient nature of the open state. A structure from the bacterium Magnetococcus marinus Nav (NavMs) was initially proposed to be open, based on its pore diameter and voltage-sensor conformation. However, the functional annotation of this model, and the structural details of the open state, remain disputed. In this work, we used molecular modeling and simulations to test possible open-state models of NavMs. The full-length experimental structure, termed here the α-model, was consistently dehydrated at the activation gate, indicating an inability to conduct ions. Based on a spontaneous transition observed in extended simulations, and sequence/structure comparison to other Nav channels, we built an alternative π-model featuring a helix transition and the rotation of a conserved asparagine residue into the activation gate. Pore hydration, ion permeation, and state-dependent drug binding in this model were consistent with an open functional state. This work thus offers both a functional annotation of the full-length NavMs structure and a detailed model for a stable Nav open state, with potential conservation in diverse ion-channel families.

Structural modeling of Nav1.5 pore domain in closed state

The voltage-dependent cardiac sodium channel plays a key role in cardiac excitability and conduction and it is the drug target of medically important. However, its atomic- resolution structure is still lack. Here, we report a modeled structure of Na_v1.5 pore domain in closed state. The structure was constructed by Rosetta-membrane homology modeling method based on the template of eukaryotic Na_v channel Na_vPaS and selected by energy and direct coupling analysis (DCA). Moreover, this structure was optimized through molecular dynamical simulation in the lipid membrane bilayer. Finally, to validate the constructed model, the binding energy and binding sites of closed-state local anesthetics (LAs) in the modeled structure were computed by the MM-GBSA method and the results are in agreement with experiments. The modeled structure of Na_v1.5 pore domain in closed state may be useful to explore molecular mechanism of a state-dependent drug binding and helpful for new drug development.

Function-Related Dynamics in Multi-Spanning Helical Membrane Proteins Revealed by Solution NMR

A primary biological function of multi-spanning membrane proteins is to transfer information and/or materials through a membrane by changing their conformations. Therefore, particular dynamics of the membrane proteins are tightly associated with their function. The semi-atomic resolution dynamics information revealed by NMR is able to discriminate function-related dynamics from random fluctuations. This review will discuss several studies in which quantitative dynamics information by solution NMR has contributed to revealing the structural basis of the function of multi-spanning membrane proteins, such as ion channels, GPCRs, and transporters.

Determinants of conductance of a bacterial voltage-gated sodium channel

Through molecular dynamics (MD) and free energy simulations in electric fields, we examine the factors influencing conductance of bacterial voltage-gated sodium channel Na_vMs. The channel utilizes four glutamic acid residues in the selectivity filter (SF). Previously, we have shown, through constant pH and free energy calculations of pKa values, that fully deprotonated, singly protonated, and doubly protonated states are all feasible at physiological pH, depending on how many ions are bound in the SF. With 173 MD simulations of 450 or 500 ns and additional free energy simulations, we determine that the conductance is highest for the deprotonated state and decreases with each additional proton bound. We also determine that the pKa value of the four glutamic residues for the transition between deprotonated and singly protonated states is close to the physiological pH and that there is a small voltage dependence. The pKa value and conductance trends are in agreement with experimental work on bacterial Na_v channels, which show a decrease in maximal conductance with lowering of pH, with pKa in the physiological range. We examine binding sites for Na⁺ in the SF, compare with previous work, and note a dependence on starting structures. We find that narrowing of the gate backbone to values lower than the crystal structure's backbone radius reduces the conductance, whereas increasing the gate radius further does not affect the conductance. Simulations with some amount of negatively charged lipids as opposed to purely neutral lipids increases the conductance, as do simulations at higher voltages.

In vitro discovery of novel prokaryotic ion channel candidates for antiarrhythmic gene therapy

Sudden cardiac death continues to have a devastating impact on public health prompting the continued efforts to develop more effective therapies for cardiac arrhythmias. Among different approaches to normalize function of ion channels and prevent arrhythmogenic remodeling of tissue substrate, cardiac cell and gene therapies are emerging as promising strategies to restore and maintain normal heart rhythm. Specifically, the ability to genetically enhance electrical excitability of diseased hearts through voltage-gated sodium channel (VGSC) gene transfer could improve velocity of action potential conduction and act to stop reentrant circuits underlying sustained arrhythmias. For this purpose, prokaryotic VGSC genes are promising therapeutic candidates due to their small size (<1kb) and potential to be effectively packaged in adeno-associated viral (AAV) vectors and delivered to cardiomyocytes for stable, long-term expression. This article describes a versatile method to discover and characterize novel prokaryotic ion channels for use in gene and cell therapies for heart disease including cardiac arrhythmias. Detailed protocols are provided for: (1) identification of potential ion channel candidates from large genomic databases, (2) candidate screening and characterization using site-directed mutagenesis and engineered human excitable cell system and, (3) candidate validation using electrophysiological techniques and an in vitro model of impaired cardiac impulse conduction.

Structural Pharmacology of Voltage-Gated Sodium Channels

Voltage-gated sodium (Na_V) channels initiate and propagate action potentials in excitable tissues to mediate key physiological processes including heart contraction and nervous system function. Accordingly, Na_V channels are major targets for drugs, toxins and disease-causing mutations. Recent breakthroughs in cryo-electron microscopy have led to the visualization of human Na_V1.1, Na_V1.2, Na_V1.4, Na_V1.5 and Na_V1.7 channel subtypes at high-resolution. These landmark studies have greatly advanced our structural understanding of channel architecture, ion selectivity, voltage-sensing, electromechanical coupling, fast inactivation, and the molecular basis underlying Na_V channelopathies. Na_V channel structures have also been increasingly determined in complex with toxin and small molecule modulators that target either the pore module or voltage sensor domains. These structural studies have provided new insights into the mechanisms of pharmacological action and opportunities for subtype-selective Na_V channel drug design. This review will highlight the structural pharmacology of human Na_V channels as well as the potential use of engineered and chimeric channels in future drug discovery efforts.

Structural basis of cytoplasmic NaV1.5 and NaV1.4 regulation

Voltage-gated sodium channels (NaVs) are membrane proteins responsible for the rapid upstroke of the action potential in excitable cells. There are nine human voltage-sensitive NaV1 isoforms that, in addition to their sequence differences, differ in tissue distribution and specific function. This review focuses on isoforms NaV1.4 and NaV1.5, which are primarily expressed in skeletal and cardiac muscle cells, respectively. The determination of the structures of several eukaryotic NaVs by single-particle cryo-electron microscopy (cryo-EM) has brought new perspective to the study of the channels. Alignment of the cryo-EM structure of the transmembrane channel pore with x-ray crystallographic structures of the cytoplasmic domains illustrates the complementary nature of the techniques and highlights the intricate cellular mechanisms that modulate these channels. Here, we review structural insights into the cytoplasmic C-terminal regulation of NaV1.4 and NaV1.5 with special attention to Ca2+ sensing by calmodulin, implications for disease, and putative channel dimerization.

Changes in Ion Selectivity Following the Asymmetrical Addition of Charge to the Selectivity Filter of Bacterial Sodium Channels

Voltage-gated sodium channels (NaVs) play fundamental roles in eukaryotes, but their exceptional size hinders their structural resolution. Bacterial NaVs are simplified homologues of their eukaryotic counterparts, but their use as models of eukaryotic Na⁺ channels is limited by their homotetrameric structure at odds with the asymmetric Selectivity Filter (SF) of eukaryotic NaVs. This work aims at mimicking the SF of eukaryotic NaVs by engineering radial asymmetry into the SF of bacterial channels. This goal was pursued with two approaches: the co-expression of different monomers of the NaChBac bacterial channel to induce the random assembly of heterotetramers, and the concatenation of four bacterial monomers to form a concatemer that can be targeted by site-specific mutagenesis. Patch-clamp measurements and Molecular Dynamics simulations showed that an additional gating charge in the SF leads to a significant increase in Na⁺ and a modest increase in the Ca²⁺ conductance in the NavMs concatemer in agreement with the behavior of the population of random heterotetramers with the highest proportion of channels with charge -5e. We thus showed that charge, despite being important, is not the only determinant of conduction and selectivity, and we created new tools extending the use of bacterial channels as models of eukaryotic counterparts.

The conformational cycle of a prototypical voltage-gated sodium channel

Electrical signaling was a dramatic development in evolution, allowing complex single-cell organisms like Paramecium to coordinate movement and early metazoans like worms and jellyfish to send regulatory signals rapidly over increasing distances. But how are electrical signals generated in biology? In fact, voltage-gated sodium channels conduct sodium currents that initiate electrical signals in all kingdoms of life, from bacteria to man. They are responsible for generating the action potential in vertebrate nerve and muscle, neuroendocrine cells, and other cell types1,2. Because of the high level of conservation of their core structure, it is likely that their fundamental mechanisms of action are conserved as well. Here we describe the complete cycle of conformational changes that a bacterial sodium channel undergoes as it transitions from resting to activated/open and inactivated/closed states, based on high-resolution structural studies of a single sodium channel. We further relate this conformational cycle of the ancestral sodium channel to the function of its vertebrate orthologs. The strong conservation of amino acid sequence and three-dimensional structure suggests that this model, at a fundamental level, is relevant for both prokaryotic and eukaryotic sodium channels, as well as voltage-gated calcium and potassium channels.

Changes in Membrane Protein Structural Biology

Membrane proteins are essential components of many biochemical processes and are important pharmaceutical targets. Membrane protein structural biology provides the molecular rationale for these biochemical process as well as being a highly useful tool for drug discovery. Unfortunately, membrane protein structural biology is a difficult area of study due to low protein yields and high levels of instability especially when membrane proteins are removed from their native environments. Despite this instability, membrane protein structural biology has made great leaps over the last fifteen years. Today, the landscape is almost unrecognisable. The numbers of available atomic resolution structures have increased 10-fold though advances in crystallography and more recently by cryo-electron microscopy. These advances in structural biology were achieved through the efforts of many researchers around the world as well as initiatives such as the Membrane Protein Laboratory (MPL) at Diamond Light Source. The MPL has helped, provided access to and contributed to advances in protein production, sample preparation and data collection. Together, these advances have enabled higher resolution structures, from less material, at a greater rate, from a more diverse range of membrane protein targets. Despite this success, significant challenges remain. Here, we review the progress made and highlight current and future challenges that will be overcome.

Membrane protein crystallography in the era of modern structural biology

The aim of structural biology has been always the study of biological macromolecules structures and their mechanistic behaviour at molecular level. To achieve its goal, multiple biophysical methods and approaches have become part of the structural biology toolbox. Considered as one of the pillars of structural biology, X-ray crystallography has been the most successful method for solving three-dimensional protein structures at atomic level to date. It is however limited by the success in obtaining well-ordered protein crystals that diffract at high resolution. This is especially true for challenging targets such as membrane proteins (MPs). Understanding structure-function relationships of MPs at the biochemical level is vital for medicine and drug discovery as they play critical roles in many cellular processes. Though difficult, structure determination of MPs by X-ray crystallography has significantly improved in the last two decades, mainly due to many relevant technological and methodological developments. Today, numerous MP crystal structures have been solved, revealing many of their mechanisms of action. Yet the field of structural biology has also been through significant technological breakthroughs in recent years, particularly in the fields of single particle electron microscopy (cryo-EM) and X-ray free electron lasers (XFELs). Here we summarise the most important advancements in the field of MP crystallography and the significance of these developments in the present era of modern structural biology.

Conservation and divergence in NaChBac and NaV1.7 pharmacology reveals novel drug interaction mechanisms

Voltage-gated Na⁺ (Na_V) channels regulate homeostasis in bacteria and control membrane electrical excitability in mammals. Compared to their mammalian counterparts, bacterial Na_V channels possess a simpler, fourfold symmetric structure and have facilitated studies of the structural basis of channel gating. However, the pharmacology of bacterial Na_V remains largely unexplored. Here we systematically screened 39 Na_V modulators on a bacterial channel (NaChBac) and characterized a selection of compounds on NaChBac and a mammalian channel (human Na_V1.7). We found that while many compounds interact with both channels, they exhibit distinct functional effects. For example, the local anesthetics ambroxol and lidocaine block both Na_V1.7 and NaChBac but affect activation and inactivation of the two channels to different extents. The voltage-sensing domain targeting toxin BDS-I increases Na_V1.7 but decreases NaChBac peak currents. The pore binding toxins aconitine and veratridine block peak currents of Na_V1.7 and shift activation (aconitine) and inactivation (veratridine) respectively. In NaChBac, they block the peak current by binding to the pore residue F224. Nonetheless, aconitine has no effect on activation or inactivation, while veratridine only modulates activation of NaChBac. The conservation and divergence in the pharmacology of bacterial and mammalian Na_V channels provide insights into the molecular basis of channel gating and will facilitate organism-specific drug discovery.

Employing NaChBac for cryo-EM analysis of toxin action on voltage-gated Na+ channels in nanodisc

NaChBac, the first bacterial voltage-gated Na⁺ (Na_v) channel to be characterized, has been the prokaryotic prototype for studying the structure-function relationship of Na_v channels. Discovered nearly two decades ago, the structure of NaChBac has not been determined. Here we present the single particle electron cryomicroscopy (cryo-EM) analysis of NaChBac in both detergent micelles and nanodiscs. Under both conditions, the conformation of NaChBac is nearly identical to that of the potentially inactivated Na_vAb. Determining the structure of NaChBac in nanodiscs enabled us to examine gating modifier toxins (GMTs) of Na_v channels in lipid bilayers. To study GMTs in mammalian Na_v channels, we generated a chimera in which the extracellular fragment of the S3 and S4 segments in the second voltage-sensing domain from Na_v1.7 replaced the corresponding sequence in NaChBac. Cryo-EM structures of the nanodisc-embedded chimera alone and in complex with HuwenToxin IV (HWTX-IV) were determined to 3.5 and 3.2 Å resolutions, respectively. Compared to the structure of HWTX-IV-bound human Na_v1.7, which was obtained at an overall resolution of 3.2 Å, the local resolution of the toxin has been improved from ∼6 to ∼4 Å. This resolution enabled visualization of toxin docking. NaChBac can thus serve as a convenient surrogate for structural studies of the interactions between GMTs and Na_v channels in a membrane environment.

Hydrophobic Drug/Toxin Binding Sites in Voltage-Dependent K+ and Na+ Channels

In the Na_v channel family the lipophilic drugs/toxins binding sites and the presence of fenestrations in the channel pore wall are well defined and categorized. No such classification exists in the much larger K_v channel family, although certain lipophilic compounds seem to deviate from binding to well-known hydrophilic binding sites. By mapping different compound binding sites onto 3D structures of Kv channels, there appear to be three distinct lipid-exposed binding sites preserved in K_v channels: the front and back side of the pore domain, and S2-S3/S3-S4 clefts. One or a combination of these sites is most likely the orthologous equivalent of neurotoxin site 5 in Na_v channels. This review describes the different lipophilic binding sites and location of pore wall fenestrations within the K_v channel family and compares it to the knowledge of Na_v channels.

Simulating ion channel activation mechanisms using swarms of trajectories

Atomic-level studies of protein activity represent a significant challenge as a result of the complexity of conformational changes occurring on wide-ranging timescales, often greatly exceeding that of even the longest simulations. A prime example is the elucidation of protein allosteric mechanisms, where localized perturbations transmit throughout a large macromolecule to generate a response signal. For example, the conversion of chemical to electrical signals during synaptic neurotransmission in the brain is achieved by specialized membrane proteins called pentameric ligand-gated ion channels. Here, the binding of a neurotransmitter results in a global conformational change to open an ion-conducting pore across the nerve cell membrane. X-ray crystallography has produced static structures of the open and closed states of the proton-gated GLIC pentameric ligand-gated ion channel protein, allowing for atomistic simulations that can uncover changes related to activation. We discuss a range of enhanced sampling approaches that could be used to explore activation mechanisms. In particular, we describe recent application of an atomistic string method, based on Roux's "swarms of trajectories" approach, to elucidate the sequence and interdependence of conformational changes during activation. We illustrate how this can be combined with transition analysis and Brownian dynamics to extract thermodynamic and kinetic information, leading to understanding of what controls ion channel function. © 2019 Wiley Periodicals, Inc.

Veratridine: A Janus-Faced Modulator of Voltage-Gated Sodium Ion Channels

Voltage-gated sodium ion channels (Na_Vs) are integral to both neuronal and muscular signaling and are a primary target for a number of proteinaceous and small molecule toxins. Included among these neurotoxins is veratridine (VTD), a C-nor-D homosteroidal alkaloid from the seeds of members of the Veratrum genus. VTD binds to Na_V within the pore region, causing a hyperpolarizing shift in the activation threshold in addition to reducing peak current. We have characterized the activity of VTD against heterologously expressed rat Na_V1.4 and have demonstrated that VTD acts on the channel as either an agonist or antagonist depending on the nature of the electrophysiological stimulation protocol. Structure-activity studies with VTD and VTD derivatives against Na_V mutants show that the functional duality of VTD can be decoupled. These findings suggest that the dichotomous activity of VTD may derive from two distinct, use-dependent binding orientations of the toxin.

Role of Voltage-Gated Sodium Channels in the Mechanism of Ether-Induced Unconsciousness

Despite continuous clinical use for more than 170 years, the mechanism of general anesthetics has not been completely characterized. In this review, we focus on the role of voltage-gated sodium channels in the sedative-hypnotic actions of halogenated ethers, describing the history of anesthetic mechanisms research, the basic neurobiology and pharmacology of voltage-gated sodium channels, and the evidence for a mechanistic interaction between halogenated ethers and sodium channels in the induction of unconsciousness. We conclude with a more integrative perspective of how voltage-gated sodium channels might provide a critical link between molecular actions of the halogenated ethers and the more distributed network-level effects associated with the anesthetized state across species.

Voltage-gated Sodium Channels and Blockers: An Overview and Where Will They Go?

Voltage-gated sodium (Nav) channels are critical players in the generation and propagation of action potentials by triggering membrane depolarization. Mutations in Nav channels are associated with a variety of channelopathies, which makes them relevant targets for pharmaceutical intervention. So far, the cryoelectron microscopic structure of the human Nav1.2, Nav1.4, and Nav1.7 has been reported, which sheds light on the molecular basis of functional mechanism of Nav channels and provides a path toward structure-based drug discovery. In this review, we focus on the recent advances in the structure, molecular mechanism and modulation of Nav channels, and state updated sodium channel blockers for the treatment of pathophysiology disorders and briefly discuss where the blockers may be developed in the future.

Valproic acid interactions with the NavMs voltage-gated sodium channel

Valproic acid (VPA) is an anticonvulsant drug that is also used to treat migraines and bipolar disorder. Its proposed biological targets include human voltage-gated sodium channels, among other membrane proteins. We used the prokaryotic NavMs sodium channel, which has been shown to be a good exemplar for drug binding to human sodium channels, to examine the structural and functional interactions of VPA. Thermal melt synchrotron radiation circular dichroism spectroscopic binding studies of the full-length NavMs channel (which includes both pore and voltage sensor domains), and a pore-only construct, undertaken in the presence and absence of VPA, indicated that the drug binds to and destabilizes the channel, but not the pore-only construct. This is in contrast to other antiepileptic compounds that have previously been shown to bind in the central hydrophobic core of the pore region of the channel, and that tend to increase the thermal stability of both pore-only constructs and full-length channels. Molecular docking studies also indicated that the VPA binding site is associated with the voltage sensor, rather than the hydrophobic cavity of the pore domain. Electrophysiological studies show that VPA influences the block and inactivation rates of the NavMs channel, although with lower efficacy than classical channel-blocking compounds. It thus appears that, while VPA is capable of binding to these voltage-gated sodium channels, it has a very different mode and site of action than other anticonvulsant compounds.

X-ray crystallography of TRP channels

Transient receptor potential (TRP) ion channels are molecular sensors of a large variety of stimuli including temperature, mechanical stress, voltage, small molecules including capsaicin and menthol, and lipids such as phosphatidylinositol 4,5-bisphosphate (PIP₂). Since the same TRP channels may respond to different physical and chemical stimuli, they can serve as signal integrators. Many TRP channels are calcium permeable and contribute to Ca²⁺ homeostasis and signaling. Although the TRP channel family was discovered decades ago, only recently have the structures of many of these channels been solved, largely by cryo-electron microscopy (cryo-EM). Complimentary to cryo-EM, X-ray crystallography provides unique tools to unambiguously identify specific atoms and can be used to study ion binding in channel pores. In this review we describe crystallographic studies of the TRP channel TRPV6. The methodology used in these studies may serve as a template for future structural analyses of different types of TRP and other ion channels.

Resting-State Structure and Gating Mechanism of a Voltage-Gated Sodium Channel

Voltage-gated sodium (NaV) channels initiate action potentials in nerve, muscle, and other electrically excitable cells. The structural basis of voltage gating is uncertain because the resting state exists only at deeply negative membrane potentials. To stabilize the resting conformation, we inserted voltage-shifting mutations and introduced a disulfide crosslink in the VS of the ancestral bacterial sodium channel NaVAb. Here, we present a cryo-EM structure of the resting state and a complete voltage-dependent gating mechanism. The S4 segment of the VS is drawn intracellularly, with three gating charges passing through the transmembrane electric field. This movement forms an elbow connecting S4 to the S4-S5 linker, tightens the collar around the S6 activation gate, and prevents its opening. Our structure supports the classical "sliding helix" mechanism of voltage sensing and provides a complete gating mechanism for voltage sensor function, pore opening, and activation-gate closure based on high-resolution structures of a single sodium channel protein.

Voltage-gated sodium channels: structures, functions, and molecular modeling

Voltage-gated sodium channels (VGSCs), formed by 24 transmembrane segments arranged into four domains, have a key role in the initiation and propagation of electrical signaling in excitable cells. VGSCs are involved in a variety of diseases, including epilepsy, cardiac arrhythmias, and neuropathic pain, and therefore have been regarded as appealing therapeutic targets for the development of anticonvulsant, antiarrhythmic, and local anesthetic drugs. In this review, we discuss recent advances in understanding the structures and biological functions of VGSCs. In addition, we systematically summarize eight pharmacologically distinct ligand-binding sites in VGSCs and representative isoform-selective VGSC modulators in clinical trials. Finally, we review studies on molecular modeling and computer-aided drug design (CADD) for VGSCs to help understanding of biological processes involving VGSCs.

Covalent linkage of bacterial voltage-gated sodium channels

Background

Bacterial sodium channels are important models for understanding ion permeation and selectivity. However, their homotetrameric structure limits their use as models for understanding the more complex eukaryotic voltage-gated sodium channels (which have a pseudo-heterotetrameric structure formed from an oligomer composed of four domains). To bridge this gap we attempted to synthesise oligomers made from four covalently linked bacterial sodium channel monomers and thus resembling their eukaryotic counterparts.

Results

Western blot analyses revealed NaChBac oligomers to be inherently unstable whereas intact expression of NavMs oligomers was possible. Immunodectection using confocal microscopy and electrophysiological characterisation of NavMs tetramers confirmed plasma membrane localisation and equivalent functionality with wild type NavMs channels when expressed in human embryonic kidney cells.

Conclusion

This study has generated new tools for the investigation of eukaryotic channels. The successful covalent linkage of four bacterial Nav channel monomers should permit the introduction of radial asymmetry into the structure of bacterial Nav channels and enable the known structures of these channels to be used to gain unique insights into structure-function relationships of their eukaryotic counterparts.

Electronic supplementary material

The online version of this article (10.1186/s13628-019-0049-5) contains supplementary material, which is available to authorized users.

Emerging Diversity in Lipid-Protein Interactions

Membrane lipids interact with proteins in a variety of ways, ranging from providing a stable membrane environment for proteins to being embedded in to detailed roles in complicated and well-regulated protein functions. Experimental and computational advances are converging in a rapidly expanding research area of lipid-protein interactions. Experimentally, the database of high-resolution membrane protein structures is growing, as are capabilities to identify the complex lipid composition of different membranes, to probe the challenging time and length scales of lipid-protein interactions, and to link lipid-protein interactions to protein function in a variety of proteins. Computationally, more accurate membrane models and more powerful computers now enable a detailed look at lipid-protein interactions and increasing overlap with experimental observations for validation and joint interpretation of simulation and experiment. Here we review papers that use computational approaches to study detailed lipid-protein interactions, together with brief experimental and physiological contexts, aiming at comprehensive coverage of simulation papers in the last five years. Overall, a complex picture of lipid-protein interactions emerges, through a range of mechanisms including modulation of the physical properties of the lipid environment, detailed chemical interactions between lipids and proteins, and key functional roles of very specific lipids binding to well-defined binding sites on proteins. Computationally, despite important limitations, molecular dynamics simulations with current computer power and theoretical models are now in an excellent position to answer detailed questions about lipid-protein interactions.