Local anesthetic and antiepileptic drug access and binding to a bacterial voltage-gated sodium channel

Inhibition of NMDA receptors and other ion channel types by membrane-associated drugs

N-methyl-D-aspartate receptors (NMDARs) are ligand-gated ion channels present at most excitatory synapses in the brain that play essential roles in cognitive functions including learning and memory consolidation. However, NMDAR dysregulation is implicated in many nervous system disorders. Diseases that involve pathological hyperactivity of NMDARs can be treated clinically through inhibition by channel blocking drugs. NMDAR channel block can occur via two known mechanisms. First, in traditional block, charged drug molecules can enter the channel directly from the extracellular solution after NMDAR activation and channel opening. Second, uncharged molecules of channel blocking drug can enter the hydrophobic plasma membrane, and upon NMDAR activation the membrane-associated drug can transit into the channel through a fenestration within the NMDAR. This membrane-associated mechanism of action is called membrane to channel inhibition (MCI) and is not well understood despite the clinical importance of NMDAR channel blocking drugs. Intriguingly, a hydrophobic route of access for drugs is not unique to NMDARs. Our review will address inhibition of NMDARs and other ion channels by membrane-associated drugs and consider how the path of access may affect a drug's therapeutic potential.

Drugs exhibit diverse binding modes and access routes in the Nav1.5 cardiac sodium channel pore

Small molecule inhibitors of the sodium channel are common pharmacological agents used to treat a variety of cardiac and nervous system pathologies. They act on the channel via binding within the pore to directly block the sodium conduction pathway and/or modulate the channel to favor a non-conductive state. Despite their abundant clinical use, we lack specific knowledge of their protein-drug interactions and the subtle variations between different compound structures. This study investigates the binding and accessibility of nine different compounds in the pore cavity of the Nav1.5 sodium channel using enhanced sampling simulations. We find that most compounds share a common location of pore binding-near the mouth of the DII-III fenestration-associated with the high number of aromatic residues in this region. In contrast, some other compounds prefer binding within the lateral fenestrations where they compete with lipids, rather than binding in the central cavity. Overall, our simulation results suggest that the drug binding within the pore is highly promiscuous, with most drugs having multiple low-affinity binding sites. Access to the pore interior via two out of four of the hydrophobic fenestrations is favorable for the majority of compounds. Our results indicate that the polyspecific and diffuse binding of inhibitors in the pore contributes to the varied nature of their inhibitory effects and can be exploited for future drug discovery and optimization.

Reactivity of Olanzapine and Tricyclic Antidepressants on the Protective Effects of Trolox on Lipid Peroxidation Evaluated Using Fluorescence Anisotropy, Electron Paramagnetic Resonance Spectrometry, and Thermal Analysis

Multiacting receptor-targeting antipsychotics and tricyclic antidepressants stimulate various neurotransmitter receptors despite the different targets of postsynaptic receptors and presynaptic reuptake transporters. Their auxiliary and adverse effects may be caused by multiple targets or the modification of the neuronal membrane. To evaluate the membrane responses to olanzapine, imipramine, desipramine, amitriptyline, lidocaine, and dibucaine, we examined the inhibition of lipid peroxidation in egg yolk phosphatidylcholine liposomes. By contrast, their effects on membrane fluidity were measured as the suppressive contributions of the inhibitory activity of Trolox on lipid oxidation. These drugs inhibit lipid peroxidation and exclude harmful reactive oxygen species and the protective effect of Trolox. The fluorescence anisotropy of 1,6-diphenyl-1,3,5-hexatriene in saturated phospholipid liposome-containing drugs suggested that olanzapine, imipramine, and dibucaine enhanced membrane fluidity. The radical scavenging activity of 2,2-diphenylpicrylhidrazyl and galvinoxyl radicals was determined using electron paramagnetic resonance experiments, and their molecular flexibility was determined using thermograms for differential scanning calorimetry. Multiple regression analyses of the linear free energy relationship approach and comparative investigations revealed that the membranous fluidity of the liposomes, independent of the radical scavenging activity of the drugs, induced the inhibitory activity on lipid peroxidation. We discussed how these drugs act on nervous membranes and aimed to identify the relationship between uncertified functions and membranous fluidity.

Nav1.2 channel mutations preventing fast inactivation lead to SCN2A encephalopathy

SCN2A gene-related early-infantile developmental and epileptic encephalopathy (EI-DEE) is a rare and severe disorder that manifests in early infancy. SCN2A mutations affecting the fast inactivation gating mechanism can result in altered voltage dependence and incomplete inactivation of the encoded neuronal Nav1.2 channel and lead to abnormal neuronal excitability. In this study, we evaluated clinical data of seven missense Nav1.2 variants associated with DEE and performed molecular dynamics simulations, patch-clamp electrophysiology and dynamic clamp real-time neuronal modelling to elucidate the molecular and neuron-scale phenotypic consequences of the mutations. The N1662D mutation almost completely prevented fast inactivation without affecting activation. The comparison of wild-type and N1662D channel structures suggested that the ambifunctional hydrogen bond formation between residues N1662 and Q1494 is essential for fast inactivation. Fast inactivation could also be prevented with engineered Q1494A or Q1494L Nav1.2 channel variants, whereas Q1494E or Q149K variants resulted in incomplete inactivation and persistent current. Molecular dynamics simulations revealed a reduced affinity of the hydrophobic IFM-motif to its receptor site with N1662D and Q1494L variants relative to wild-type. These results demonstrate that the interactions between N1662 and Q1494 underpin the stability and the orientation of the inactivation gate and are essential for the development of fast inactivation. Six DEE-associated Nav1.2 variants, with mutations mapped to channel segments known to be implicated in fast inactivation were also evaluated. Remarkably, the L1657P variant also prevented fast inactivation and produced biophysical characteristics that were similar to those of N1662D, whereas the M1501V, M1501T, F1651C, P1658S and A1659V variants resulted in biophysical properties that were consistent with gain-of-function and enhanced action potential firing of hybrid neurons in dynamic action potential clamp experiments. Paradoxically, low densities of N1662D or L1657P currents potentiated action potential firing, whereas increased densities resulted in sustained depolarization. Our results provide novel structural insights into the molecular mechanism of Nav1.2 channel fast inactivation and inform treatment strategies for SCN2A-related EI-DEE. The contribution of non-inactivating Nav1.2 channels to neuronal excitability may constitute a distinct cellular mechanism in the pathogenesis of SCN2A-related DEE.

Voltage-Gated Ion Channels: Structure, Pharmacology and Photopharmacology

Voltage-gated ion channels are transmembrane proteins responsible for the generation and propagation of action potentials in excitable cells. Over the last decade, advancements have enabled the elucidation of crystal structures of ion channels. This progress in structural understanding, particularly in identifying the binding sites of local anesthetics, opens avenues for the design of novel compounds capable of modulating ion conduction. However, many traditional drugs lack selectivity and come with adverse side effects. The emergence of photopharmacology has provided an orthogonal way of controlling the activity of compounds, enabling the regulation of ion conduction with light. In this review, we explore the central pore region of voltage-gated sodium and potassium channels, providing insights from both structural and pharmacological perspectives. We discuss the different binding modes of synthetic compounds that can physically occlude the pore and, therefore, block ion conduction. Moreover, we examine recent advances in the photopharmacology of voltage-gated ion channels, introducing molecular approaches aimed at controlling their activity by using photosensitive drugs.

In Silico Screening Identification of Fatty Acids and Fatty Acid Derivatives with Antiseizure Activity: In Vitro and In Vivo Validation

High fat diets have been used as complementary treatments for seizure disorders for more than a century. Moreover, many fatty acids and derivatives, including the broad-spectrum antiseizure medication valproic acid, have been explored and used as pharmacological agents to treat epilepsy. In this work, we have explored the anticonvulsant potential of a large library of fatty acids and fatty acid derivatives, the LIPID MAPS Structure Database, using structure-based virtual screening to assess their ability to block the voltage-gated sodium channel 1.2 (NaV1.2), a validated target for antiseizure medications. Four of the resulting in silico hits were submitted for experimental confirmation using in vitro patch clamp experiments, and their protective role was evaluated in an acute mice seizure model, the Maximal Electroshock seizure model. These four compounds were found to protect mice against seizures. Two of them exhibited blocking effects on NaV1.2, CaV2.2, and CaV3.1.

Unplugging lateral fenestrations of NALCN reveals a hidden drug binding site within the pore region

The sodium (Na+) leak channel (NALCN) is a member of the four-domain voltage-gated cation channel family that includes the prototypical voltage-gated sodium and calcium channels (NaVs and CaVs, respectively). Unlike NaVs and CaVs, which have four lateral fenestrations that serve as routes for lipophilic compounds to enter the central cavity to modulate channel function, NALCN has bulky residues (W311, L588, M1145, and Y1436) that block these openings. Structural data suggest that occluded fenestrations underlie the pharmacological resistance of NALCN, but functional evidence is lacking. To test this hypothesis, we unplugged the fenestrations of NALCN by substituting the four aforementioned residues with alanine (AAAA) and compared the effects of NaV, CaV, and NALCN blockers on both wild-type (WT) and AAAA channels. Most compounds behaved in a similar manner on both channels, but phenytoin and 2-aminoethoxydiphenyl borate (2-APB) elicited additional, distinct responses on AAAA channels. Further experiments using single alanine mutants revealed that phenytoin and 2-APB enter the inner cavity through distinct fenestrations, implying structural specificity to their modes of access. Using a combination of computational and functional approaches, we identified amino acid residues critical for 2-APB activity, supporting the existence of drug binding site(s) within the pore region. Intrigued by the activity of 2-APB and its analogues, we tested compounds containing the diphenylmethane/amine moiety on WT channels. We identified clinically used drugs that exhibited diverse activity, thus expanding the pharmacological toolbox for NALCN. While the low potencies of active compounds reiterate the pharmacological resistance of NALCN, our findings lay the foundation for rational drug design to develop NALCN modulators with refined properties.

Unplugging lateral fenestrations of NALCN reveals a hidden drug binding site within the pore module

The sodium (Na + ) leak channel (NALCN) is a member of the four-domain voltage-gated cation channel family that includes the prototypical voltage-gated sodium and calcium channels (Na V s and Ca V s, respectively). Unlike Na V s and Ca V s, which have four lateral fenestrations that serve as routes for lipophilic compounds to enter the central cavity to modulate channel function, NALCN has bulky residues (W311, L588, M1145 and Y1436) that block these openings. Structural data suggest that oc-cluded lateral fenestrations underlie the pharmacological resistance of NALCN to lipophilic compounds, but functional evidence is lacking. To test this hypothesis, we unplugged the fenestrations of NALCN by substituting the four aforementioned resi-dues with alanine (AAAA) and compared the effects of Na V , Ca V and NALCN block-ers on both wild-type (WT) and AAAA channels. Most compounds behaved in a simi-lar manner on both channels, but phenytoin and 2-aminoethoxydiphenyl borate (2-APB) elicited additional, distinct responses on AAAA channels. Further experiments using single alanine mutants revealed that phenytoin and 2-APB enter the inner cav-ity through distinct fenestrations, implying structural specificity to their modes of ac-cess. Using a combination of computational and functional approaches, we identified amino acid residues critical for 2-APB activity, supporting the existence of drug bind-ing site(s) within the pore region. Intrigued by the activity of 2-APB and its ana-logues, we tested additional compounds containing the diphenylmethane/amine moiety on WT channels. We identified compounds from existing clinically used drugs that exhibited diverse activity, thus expanding the pharmacological toolbox for NALCN. While the low potencies of active compounds reiterate the resistance of NALCN to pharmacological targeting, our findings lay the foundation for rational drug design to develop NALCN modulators with refined properties.

Significance statement

The sodium leak channel (NALCN) is essential for survival: mutations cause life-threatening developmental disorders in humans. However, no treatment is currently available due to the resistance of NALCN to pharmacological targeting. One likely reason is that the lateral fenestrations, a common route for clinically used drugs to enter and block related ion channels, are occluded in NALCN. Using a combination of computational and functional approaches, we unplugged the fenestrations of NALCN which led us to the first molecularly defined drug binding site within the pore region. Besides that, we also identified additional NALCN modulators from existing clinically used therapeutics, thus expanding the pharmacological toolbox for this leak channel.

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.

Cannabidiol inhibits Nav channels through two distinct binding sites

Cannabidiol (CBD), a major non-psychoactive phytocannabinoid in cannabis, is an effective treatment for some forms of epilepsy and pain. At high concentrations, CBD interacts with a huge variety of proteins, but which targets are most relevant for clinical actions is still unclear. Here we show that CBD interacts with Nav1.7 channels at sub-micromolar concentrations in a state-dependent manner. Electrophysiological experiments show that CBD binds to the inactivated state of Nav1.7 channels with a dissociation constant of about 50 nM. The cryo-EM structure of CBD bound to Nav1.7 channels reveals two distinct binding sites. One is in the IV-I fenestration near the upper pore. The other binding site is directly next to the inactivated "wedged" position of the Ile/Phe/Met (IFM) motif on the short linker between repeats III and IV, which mediates fast inactivation. Consistent with producing a direct stabilization of the inactivated state, mutating residues in this binding site greatly reduced state-dependent binding of CBD. The identification of this binding site may enable design of compounds with improved properties compared to CBD itself.

Mechanosensitive pore opening of a prokaryotic voltage-gated sodium channel

Voltage-gated ion channels (VGICs) orchestrate electrical activities that drive mechanical functions in contractile tissues such as the heart and gut. In turn, contractions change membrane tension and impact ion channels. VGICs are mechanosensitive, but the mechanisms of mechanosensitivity remain poorly understood. Here, we leverage the relative simplicity of NaChBac, a prokaryotic voltage-gated sodium channel from Bacillus halodurans, to investigate mechanosensitivity. In whole-cell experiments on heterologously transfected HEK293 cells, shear stress reversibly altered the kinetic properties of NaChBac and increased its maximum current, comparably to the mechanosensitive eukaryotic sodium channel NaV1.5. In single-channel experiments, patch suction reversibly increased the open probability of a NaChBac mutant with inactivation removed. A simple kinetic mechanism featuring a mechanosensitive pore opening transition explained the overall response to force, whereas an alternative model with mechanosensitive voltage sensor activation diverged from the data. Structural analysis of NaChBac identified a large displacement of the hinged intracellular gate, and mutagenesis near the hinge diminished NaChBac mechanosensitivity, further supporting the proposed mechanism. Our results suggest that NaChBac is overall mechanosensitive due to the mechanosensitivity of a voltage-insensitive gating step associated with the pore opening. This mechanism may apply to eukaryotic VGICs, including NaV1.5.

State-Dependent Inhibition of Voltage-Gated Sodium Channels in Neuroblastoma Neuro-2A Cells by Arachidonic Acid from Halichondria okadai

Voltage-gated sodium channels (Na_v) are closely associated with epilepsy, cardiac and skeletal muscle diseases, and neuropathic pain. Several toxic compounds have been isolated from the marine sponge Halichondria okadai; however, toxic substances that modulate Na_v are yet to be identified. This study aimed to identify Na_v inhibitors from two snake venoms and H. okadai using mouse neuroblastoma Neuro-2A cells (N2A), which primarily express the specific Na_v subtype Na_v1.7, using whole-cell patch-clamp recordings. We successfully isolated arachidonic acid (AA, 1) from the hexane extract of H. okadai, and then the fatty acid-mediated modulation of Na_v in N2A was investigated in detail for the first time. Octanoic acid (2), palmitic acid (3), and oleic acid (4) showed no inhibitory activity at 100 μM, whereas AA (1), dihomo-γ-linolenic acid (DGLA, 5), and eicosapentaenoic acid (EPA, 6) showed IC₅₀ values of 6.1 ± 2.0, 58 ± 19, and 25 ± 4.0 μM, respectively (N = 4, mean ± SEM). Structure and activity relationships were investigated for the first time using two ω-3 polyunsaturated fatty acids (PUFAs), EPA (6) and eicosatetraenoic acid (ETA, 7), and two ω-6 PUFAs, AA (1) and DGLA (5), to determine their effects on a resting state, activated state, and inactivated state. Steady-state analysis showed that the half inactivation potential was largely hyperpolarized by 10 μM AA (1), while 50 μM DGLA (5), 50 μM EPA (6), and 10 μM ETA (7) led to a slight change. The percentages of the resting state block were 24 ± 1, 22 ± 1, 34 ± 4, and 38 ± 9% in the presence of AA (1), DGLA (5), EPA (6), and ETA (7), respectively, with EPA (6) and ETA (7) exhibiting a greater inhibition than both AA (1) and DGLA (5), and their inhibitions did not increase in the following depolarization pulses. None of the compounds exhibited the use-dependent block. The half recovery times from the inactivated state for the control, AA (1), DGLA (5), EPA (6), and ETA (7) were 7.67 ± 0.33, 34.3 ± 1.10, 15.5 ± 1.10, 10.7 ± 0.31, and 3.59 ± 0.18 ms, respectively, with AA (1) exhibiting a distinctively large effect. Overall, distributed binding to the resting and the inactivated states of Na_v would be significant for the inhibition of Na_v, which presumably depends on the active structure of each PUFA.

ARumenamides: A novel class of potential antiarrhythmic compounds

Background: Most therapeutics targeting cardiac voltage-gated sodium channels (Nav1.5) attenuate the sodium current (INa) conducted through the pore of the protein. Whereas these drugs may be beneficial for disease states associated with gain-of-function (GoF) in Nav1.5, few attempts have been made to therapeutically treat loss-of-function (LoF) conditions. The primary impediment to designing efficacious therapies for LoF is a tendency for drugs to occlude the Nav1.5 central pore. We hypothesized that molecular candidates with a high affinity for the fenestrations would potentially reduce pore block.

Methods and Results: Virtual docking was performed on 21 compounds, selected based on their affinity for the fenestrations in Nav1.5, which included a class of sulfonamides and carboxamides we identify as ARumenamide (AR). Six ARs, AR-051, AR-189, AR-674, AR-802, AR-807 and AR-811, were further docked against Nav1.5 built on NavAb and rNav1.5. Based on the virtual docking results, these particular ARs have a high affinity for Domain III-IV and Domain VI-I fenestrations. Upon functional characterization, a trend was observed in the effects of the six ARs on INa. An inverse correlation was established between the aromaticity of the AR’s functional moieties and compound block. Due to its aromaticity, AR-811 blocked INa the least compared with other aromatic ARs, which also decelerated fast inactivation onset. AR-674, with its aliphatic functional group, significantly suppresses INa and enhances use-dependence in Nav1.5. AR-802 and AR-811, in particular, decelerated fast inactivation kinetics in the most common Brugada Syndrome Type 1 and Long-QT Syndrome Type 3 mutant, E1784K, without affecting peak or persistent INa.

Conclusion: Our hypothesis that LoF in Nav1.5 may be therapeutically treated was supported by the discovery of ARs, which appear to preferentially block the fenestrations. ARs with aromatic functional groups as opposed to aliphatic groups efficaciously maintained Nav1.5 availability. We predict that these bulkier side groups may have a higher affinity for the hydrophobic milieu of the fenestrations, remaining there rather than in the central pore of the channel. Future refinements of AR compound structures and additional validation by molecular dynamic simulations and screening against more Brugada variants will further support their potential benefits in treating certain LoF cardiac arrhythmias.

Inhibition of NMDA receptors through a membrane-to-channel path

N-methyl-d-aspartate receptors (NMDARs) are transmembrane proteins that are activated by the neurotransmitter glutamate and are found at most excitatory vertebrate synapses. NMDAR channel blockers, an antagonist class of broad pharmacological and clinical significance, inhibit by occluding the NMDAR ion channel. A vast literature demonstrates that NMDAR channel blockers, including MK-801, phencyclidine, ketamine, and the Alzheimer’s disease drug memantine, can bind and unbind only when the NMDAR channel is open. Here we use electrophysiological recordings from transfected tsA201 cells and cultured neurons, NMDAR structural modeling, and custom-synthesized compounds to show that NMDAR channel blockers can enter the channel through two routes: the well-known hydrophilic path from extracellular solution to channel through the open channel gate, and also a hydrophobic path from plasma membrane to channel through a gated fenestration (“membrane-to-channel inhibition” (MCI)). Our demonstration that ligand-gated channels are subject to MCI, as are voltage-gated channels, highlights the broad expression of this inhibitory mechanism.

Wilcox et al. (2022) show that NMDA receptor channel blockers, some of which are clinically important drugs, can access their binding site via 2 routes: a well-known path from the extracellular solution, and another path through the plasma membrane.

Structural Advances in Voltage-Gated Sodium Channels

Voltage-gated sodium (Na_V) channels are responsible for the rapid rising-phase of action potentials in excitable cells. Over 1,000 mutations in Na_V channels are associated with human diseases including epilepsy, periodic paralysis, arrhythmias and pain disorders. Natural toxins and clinically-used small-molecule drugs bind to Na_V channels and modulate their functions. Recent advances from cryo-electron microscopy (cryo-EM) structures of Na_V channels reveal invaluable insights into the architecture, activation, fast inactivation, electromechanical coupling, ligand modulation and pharmacology of eukaryotic Na_V channels. These structural analyses not only demonstrate molecular mechanisms for Na_V channel structure and function, but also provide atomic level templates for rational development of potential subtype-selective therapeutics. In this review, we summarize recent structural advances of eukaryotic Na_V channels, highlighting the structural features of eukaryotic Na_V channels as well as distinct modulation mechanisms by a wide range of modulators from natural toxins to synthetic small-molecules.

Local anesthetics via multicomponent reactions

Local anesthetics occupy a prime position in clinical medicine as they temporarily relieve the pain by blocking the voltage-gated sodium channels. However, limited structural diversity, problems with the efficiency of syntheses and increasing toxicity, mean that alternative scaffolds with improved chemical syntheses are urgently needed. Here, we demonstrate an MCR-based approach both towards the synthesis of commercial local anesthetics and towards novel derivatives as potential anesthesia candidates via scaffold hopping. The reactions are efficient and scalable and several single-crystal structures have been obtained. In addition, our methodology has been applied to the synthesis of the antianginal drug ranolazine, via an Ugi three-component reaction. Representative derivatives from our libraries were evaluated as neuronal activity inhibitors using local field potential recordings (LFPs) in mouse hippocampal brain slices and showed very promising results. This study highlights new opportunities in drug discovery targeting local anesthetics.

Sodium Channels and Local Anesthetics-Old Friends With New Perspectives

The long history of local anesthetics (LAs) starts out in the late 19th century when the content of coca plant leaves was discovered to alleviate pain. Soon after, cocaine was established and headed off to an infamous career as a substance causing addiction. Today, LAs and related substances-in modified form-are indispensable in our clinical everyday life for pain relief during and after minor and major surgery, and dental practices. In this review, we elucidate on the interaction of modern LAs with their main target, the voltage-gated sodium channel (Navs), in the light of the recently published channel structures. Knowledge of the 3D interaction sites of the drug with the protein will allow to mechanistically substantiate the comprehensive data available on LA gating modification. In the 1970s it was suggested that LAs can enter the channel pore from the lipid phase, which was quite prospective at that time. Today we know from cryo-electron microscopy structures and mutagenesis experiments, that indeed Navs have side fenestrations facing the membrane, which are likely the entrance for LAs to induce tonic block. In this review, we will focus on the effects of LA binding on fast inactivation and use-dependent inhibition in the light of the proposed new allosteric mechanism of fast inactivation. We will elaborate on subtype and species specificity and provide insights into modelling approaches that will help identify the exact molecular binding orientation, access pathways and pharmacokinetics. With this comprehensive overview, we will provide new perspectives in the use of the drug, both clinically and as a tool for basic ion channel research.

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.

Membrane-Mediated Activity of Local Anesthetics

The activity of local anesthetics (LAs) has been attributed to the inhibition of ion channels, causing anesthesia. However, there is a growing body of research showing that LAs act on a wide range of receptors and channel proteins far beyond simple analgesia. The current concept of ligand recognition may no longer explain the multitude of protein targets influenced by LAs. We hypothesize that LAs can cause anesthesia without directly binding to the receptor proteins just by changing the physical properties of the lipid bilayer surrounding these proteins and ion channels based on LAs' amphiphilicity. It is possible that LAs act in one of the following ways: They 1) dissolve raft-like membrane microdomains, 2) impede nerve impulse propagation by lowering the lipid phase transition temperature, or 3) modulate the lateral pressure profile of the lipid bilayer. This could also explain the numerous additional effects of LAs besides anesthesia. Furthermore, the concepts of membrane-mediated activity and binding to ion channels do not have to exclude each other. If we were to consider LA as the middle part of a continuum between unspecific membrane-mediated activity on one end and highly specific ligand binding on the other end, we could describe LA as the link between the unspecific action of general anesthetics and toxins with their highly specific receptor binding. This comprehensive membrane-mediated model offers a fresh perspective to clinical and pharmaceutical research and therapeutic applications of local anesthetics. SIGNIFICANCE STATEMENT: Local anesthetics, according to the World Health Organization, belong to the most important drugs available to mankind. Their rediscovery as therapeutics and not only anesthetics marks a milestone in global pain therapy. The membrane-mediated mechanism of action proposed in this review can explain their puzzling variety of target proteins and their thus far inexplicable therapeutic effects. The new concept presented here places LAs on a continuum of structures and molecular mechanisms in between small general anesthetics and the more complex molecular toxins.

Open-state structure and pore gating mechanism of the cardiac sodium channel

The heartbeat is initiated by voltage-gated sodium channel Na_V1.5, which opens rapidly and triggers the cardiac action potential; however, the structural basis for pore opening remains unknown. Here, we blocked fast inactivation with a mutation and captured the elusive open-state structure. The fast inactivation gate moves away from its receptor, allowing asymmetric opening of pore-lining S6 segments, which bend and rotate at their intracellular ends to dilate the activation gate to ∼10 Å diameter. Molecular dynamics analyses predict physiological rates of Na⁺ conductance. The open-state pore blocker propafenone binds in a high-affinity pose, and drug-access pathways are revealed through the open activation gate and fenestrations. Comparison with mutagenesis results provides a structural map of arrhythmia mutations that target the activation and fast inactivation gates. These results give atomic-level insights into molecular events that underlie generation of the action potential, open-state drug block, and fast inactivation of cardiac sodium channels, which initiate the heartbeat.

Characterization of Compound-Specific, Concentration-Independent Biophysical Properties of Sodium Channel Inhibitor Mechanism of Action Using Automated Patch-Clamp Electrophysiology

We have developed an automated patch-clamp protocol that allows high information content screening of sodium channel inhibitor compounds. We have observed that individual compounds had their specific signature patterns of inhibition, which were manifested irrespective of the concentration. Our aim in this study was to quantify these properties. Primary biophysical data, such as onset rate, the shift of the half inactivation voltage, or the delay of recovery from inactivation, are concentration-dependent. We wanted to derive compound-specific properties, therefore, we had to neutralize the effect of concentration. This study describes how this is done, and shows how compound-specific properties reflect the mechanism of action, including binding dynamics, cooperativity, and interaction with the membrane phase. We illustrate the method using four well-known sodium channel inhibitor compounds, riluzole, lidocaine, benzocaine, and bupivacaine. Compound-specific biophysical properties may also serve as a basis for deriving parameters for kinetic modeling of drug action. We discuss how knowledge about the mechanism of action may help to predict the frequency-dependence of individual compounds, as well as their potential persistent current component selectivity. The analysis method described in this study, together with the experimental protocol described in the accompanying paper, allows screening for inhibitor compounds with specific kinetic properties, or with specific mechanisms of inhibition.

Molecular determinants of pro-arrhythmia proclivity of d- and l-sotalol via a multi-scale modeling pipeline

Drug isomers may differ in their proarrhythmia risk. An interesting example is the drug sotalol, an antiarrhythmic drug comprising d- and l- enantiomers that both block the hERG cardiac potassium channel and confer differing degrees of proarrhythmic risk. We developed a multi-scale in silico pipeline focusing on hERG channel – drug interactions and used it to probe and predict the mechanisms of pro-arrhythmia risks of the two enantiomers of sotalol. Molecular dynamics (MD) simulations predicted comparable hERG channel binding affinities for d- and l-sotalol, which were validated with electrophysiology experiments. MD derived thermodynamic and kinetic parameters were used to build multi-scale functional computational models of cardiac electrophysiology at the cell and tissue scales. Functional models were used to predict inactivated state binding affinities to recapitulate electrocardiogram (ECG) QT interval prolongation observed in clinical data. Our study demonstrates how modeling and simulation can be applied to predict drug effects from the atom to the rhythm for dl-sotalol and also increased proarrhythmia proclivity of d- vs. l-sotalol when accounting for stereospecific beta-adrenergic receptor blocking.

Local anesthetics impair the growth and self-renewal of glioblastoma stem cells by inhibiting ZDHHC15-mediated GP130 palmitoylation

BACKGROUND

A large number of preclinical studies have shown that local anesthetics have a direct inhibitory effect on tumor biological activities, including cell survival, proliferation, migration, and invasion. There are few studies on the role of local anesthetics in cancer stem cells. This study aimed to determine the possible role of local anesthetics in glioblastoma stem cell (GSC) self-renewal and the underlying molecular mechanisms.

METHODS

The effects of local anesthetics in GSCs were investigated through in vitro and in vivo assays (i.e., Cell Counting Kit 8, spheroidal formation assay, double immunofluorescence, western blot, and xenograft model). The acyl-biotin exchange method (ABE) assay was identified proteins that are S-acylated by zinc finger Asp-His-His-Cys-type palmitoyltransferase 15 (ZDHHC15). Western blot, co-immunoprecipitation, and liquid chromatograph mass spectrometer-mass spectrometry assays were used to explore the mechanisms of ZDHHC15 in effects of local anesthetics in GSCs.

RESULTS

In this study, we identified a novel mechanism through which local anesthetics can damage the malignant phenotype of glioma. We found that local anesthetics prilocaine, lidocaine, procaine, and ropivacaine can impair the survival and self-renewal of GSCs, especially the classic glioblastoma subtype. These findings suggest that local anesthetics may weaken ZDHHC15 transcripts and decrease GP130 palmitoylation levels and membrane localization, thus inhibiting the activation of IL-6/STAT3 signaling.

CONCLUSIONS

In conclusion, our work emphasizes that ZDHHC15 is a candidate therapeutic target, and local anesthetics are potential therapeutic options for glioblastoma.

The mechanism of non-blocking inhibition of sodium channels revealed by conformation-selective photolabeling

BACKGROUND AND PURPOSE

Sodium channel inhibitors can be used to treat hyperexcitability-related diseases, including epilepsies, pain syndromes, neuromuscular disorders and cardiac arrhythmias. The applicability of these drugs is limited by their nonspecific effect on physiological function. They act mainly by sodium channel block and in addition by modulation of channel kinetics. While channel block inhibits healthy and pathological tissue equally, modulation can preferentially inhibit pathological activity. An ideal drug designed to target the sodium channels of pathological tissue would act predominantly by modulation. Thus far, no such drug has been described.

EXPERIMENTAL APPROACH

Patch-clamp experiments with ultra-fast solution exchange and photolabeling-coupled electrophysiology were applied to describe the unique mechanism of riluzole on Nav1.4 sodium channels. In silico docking experiments were used to study the molecular details of binding.

KEY RESULTS

We present evidence that riluzole acts predominantly by non-blocking modulation. We propose that, being a relatively small molecule, riluzole is able to stay bound to the binding site, but nonetheless stay off the conduction pathway, by residing in one of the fenestrations. We demonstrate how this mechanism can be recognized.

CONCLUSIONS AND IMPLICATIONS

Our results identify riluzole as the prototype of this new class of sodium channel inhibitors. Drugs of this class are expected to selectively prevent hyperexcitability, while having minimal effect on cells firing at a normal rate from a normal resting potential.

Complex effects of eslicarbazepine on inhibitory micro networks in chronic experimental epilepsy

OBJECTIVE

Many antiseizure drugs (ASDs) act on voltage-dependent sodium channels, and the molecular basis of these effects is well established. In contrast, how ASDs act on the level of neuronal networks is much less understood.

METHODS

In the present study, we determined the effects of eslicarbazepine (S-Lic) on different types of inhibitory neurons, as well as inhibitory motifs. Experiments were performed in hippocampal slices from both sham-control and chronically epileptic pilocarpine-treated rats.

RESULTS

We found that S-Lic causes an unexpected reduction of feed-forward inhibition in the CA1 region at high concentrations (300 µM), but not at lower concentrations (100 µM). Concurrently, 300 but not 100 μM S-Lic significantly reduced maximal firing rates in putative feed-forward interneurons located in the CA1 stratum radiatum of sham-control and epileptic animals. In contrast, feedback inhibition was not inhibited by S-Lic. Instead, application of S-Lic, in contrast to previous data for other drugs like carbamazepine (CBZ), resulted in a lasting potentiation of feedback inhibitory post-synaptic currents (IPSCs) only in epileptic and not in sham-control animals, which persisted after washout of S-Lic. We hypothesized that this plasticity of inhibition might rely on anti-Hebbian potentiation of excitatory feedback inputs onto oriens-lacunosum moleculare (OLM) interneurons, which is dependent on Ca²⁺ -permeable α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) receptors. Indeed, we show that blocking Ca²⁺ -permeable AMPA receptors completely prevents upmodulation of feedback inhibition.

SIGNIFICANCE

These results suggest that S-Lic affects inhibitory circuits in the CA1 hippocampal region in unexpected ways. In addition, ASD actions may not be sufficiently explained by acute effects on their target channels, rather, it may be necessary to take plasticity of inhibitory circuits into account.

Structural basis for voltage-sensor trapping of the cardiac sodium channel by a deathstalker scorpion toxin

Voltage-gated sodium (Na_V) channels initiate action potentials in excitable cells, and their function is altered by potent gating-modifier toxins. The α-toxin LqhIII from the deathstalker scorpion inhibits fast inactivation of cardiac Na_V1.5 channels with IC₅₀ = 11.4 nM. Here we reveal the structure of LqhIII bound to Na_V1.5 at 3.3 Å resolution by cryo-EM. LqhIII anchors on top of voltage-sensing domain IV, wedged between the S1-S2 and S3-S4 linkers, which traps the gating charges of the S4 segment in a unique intermediate-activated state stabilized by four ion-pairs. This conformational change is propagated inward to weaken binding of the fast inactivation gate and favor opening the activation gate. However, these changes do not permit Na⁺ permeation, revealing why LqhIII slows inactivation of Na_V channels but does not open them. Our results provide important insights into the structural basis for gating-modifier toxin binding, voltage-sensor trapping, and fast inactivation of Na_V channels.

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.

Protein-ligand binding with the coarse-grained Martini model

The detailed understanding of the binding of small molecules to proteins is the key for the development of novel drugs or to increase the acceptance of substrates by enzymes. Nowadays, computer-aided design of protein–ligand binding is an important tool to accomplish this task. Current approaches typically rely on high-throughput docking essays or computationally expensive atomistic molecular dynamics simulations. Here, we present an approach to use the recently re-parametrized coarse-grained Martini model to perform unbiased millisecond sampling of protein–ligand interactions of small drug-like molecules. Remarkably, we achieve high accuracy without the need of any a priori knowledge of binding pockets or pathways. Our approach is applied to a range of systems from the well-characterized T4 lysozyme over members of the GPCR family and nuclear receptors to a variety of enzymes. The presented results open the way to high-throughput screening of ligand libraries or protein mutations using the coarse-grained Martini model.

Computer-aided design of protein-ligand binding is important for the development of novel drugs. Here authors present an approach to use the recently re-parametrized coarse-grained Martini model to perform unbiased millisecond sampling of protein-ligand binding interactions of small drug-like molecules.

A Computational Pipeline to Predict Cardiotoxicity: From the Atom to the Rhythm

RATIONALE

Drug-induced proarrhythmia is so tightly associated with prolongation of the QT interval that QT prolongation is an accepted surrogate marker for arrhythmia. But QT interval is too sensitive a marker and not selective, resulting in many useful drugs eliminated in drug discovery.

OBJECTIVE

To predict the impact of a drug from the drug chemistry on the cardiac rhythm.

METHODS AND RESULTS

In a new linkage, we connected atomistic scale information to protein, cell, and tissue scales by predicting drug-binding affinities and rates from simulation of ion channel and drug structure interactions and then used these values to model drug effects on the hERG channel. Model components were integrated into predictive models at the cell and tissue scales to expose fundamental arrhythmia vulnerability mechanisms and complex interactions underlying emergent behaviors. Human clinical data were used for model framework validation and showed excellent agreement, demonstrating feasibility of a new approach for cardiotoxicity prediction.

CONCLUSIONS

We present a multiscale model framework to predict electrotoxicity in the heart from the atom to the rhythm. Novel mechanistic insights emerged at all scales of the system, from the specific nature of proarrhythmic drug interaction with the hERG channel, to the fundamental cellular and tissue-level arrhythmia mechanisms. Applications of machine learning indicate necessary and sufficient parameters that predict arrhythmia vulnerability. We expect that the model framework may be expanded to make an impact in drug discovery, drug safety screening for a variety of compounds and targets, and in a variety of regulatory processes.

Design, synthesis and biological evaluation of N-substituted α-hydroxyimides and 1,2,3-oxathiazolidine-4-one-2,2-dioxides with anticonvulsant activity

In this investigation, we studied a family of compounds with an oxathiazolidine-4-one-2,2-dioxide skeleton and their amide synthetic precursors as new anticonvulsant drugs. The cyclic structures were synthesized using a three-step protocol that include solvent-free reactions and microwave-assisted heating. The compounds were tested in vivo through maximal electroshock seizure test in mice. All the structures showed activity at the lower doses tested (30 mg/Kg) and no signs of neurotoxicity were detected. Compound encoded as 1g displayed strong anticonvulsant effects in comparison with known anticonvulsants (ED50 = 29 mg/Kg). First approximations about the mechanisms of action of the cyclic structures were proposed by docking simulations and in vitro assays against sodium channels (patch clamp methods).

A mechanism for the activation of the mechanosensitive Piezo1 channel by the small molecule Yoda1

Mechanosensitive Piezo1 and Piezo2 channels transduce various forms of mechanical forces into cellular signals that play vital roles in many important biological processes in vertebrate organisms. Besides mechanical forces, Piezo1 is selectively activated by micromolar concentrations of the small molecule Yoda1 through an unknown mechanism. Here, using a combination of all-atom molecular dynamics simulations, calcium imaging and electrophysiology, we identify an allosteric Yoda1 binding pocket located in the putative mechanosensory domain, approximately 40 Å away from the central pore. Our simulations further indicate that the presence of the agonist correlates with increased tension-induced motions of the Yoda1-bound subunit. Our results suggest a model wherein Yoda1 acts as a molecular wedge, facilitating force-induced conformational changes, effectively lowering the channel's mechanical threshold for activation. The identification of an allosteric agonist binding site in Piezo1 channels will pave the way for the rational design of future Piezo modulators with clinical value.

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.

Characterization of Lipid-Protein Interactions and Lipid-Mediated Modulation of Membrane Protein Function through Molecular Simulation

The cellular membrane constitutes one of the most fundamental compartments of a living cell, where key processes such as selective transport of material and exchange of information between the cell and its environment are mediated by proteins that are closely associated with the membrane. The heterogeneity of lipid composition of biological membranes and the effect of lipid molecules on the structure, dynamics, and function of membrane proteins are now widely recognized. Characterization of these functionally important lipid-protein interactions with experimental techniques is however still prohibitively challenging. Molecular dynamics (MD) simulations offer a powerful complementary approach with sufficient temporal and spatial resolutions to gain atomic-level structural information and energetics on lipid-protein interactions. In this review, we aim to provide a broad survey of MD simulations focusing on exploring lipid-protein interactions and characterizing lipid-modulated protein structure and dynamics that have been successful in providing novel insight into the mechanism of membrane protein function.

A pharmacodynamic model of respiratory rate and end-tidal carbon dioxide values during anesthesia in children

It is essential to monitor the end-tidal carbon dioxide (ETCO₂) during general anesthesia and adjust the tidal volume and respiratory rate (RR). For the purpose of this study, we used a population pharmacodynamic modeling approach to establish the relationship between RR versus ETCO₂ data during general anesthesia in children, and to identify the clinical variables affecting this relationship. A prospective observational study was designed to include 51 patients (aged ≤ 12 years), including users of antiepileptic drugs (levetiracetam, valproic, or phenobarbital (n = 21)) and non-users (n = 30), scheduled to receive general anesthesia during elective surgery. When the ETCO₂ was at 40 mmHg, the RR was adjusted 1 breath per every 2 min until the ETCO₂ was 30 mmHg and recovered to 40 mmHg. Pharmacodynamic analysis using a sigmoid E_max model was performed to assess the RR-ETCO₂ relationship. As RR varied from 3 to 37 breaths per minute, the ETCO₂ changed from 40 to 30 mmHg. Hysteresis between the RR and ETCO₂ was observed and accounted for when the model was developed. The C_e50 (RR to achieve 50% of maximum decrease in ETCO₂; i.e. 35 mmHg) was 20.5 in non-users of antiepileptic drugs and 14.9 in those on antiepileptic drug medication. The values of γ (the steepness of the concentration-response relation curve) and k_eo (the first-order rate constant determining the equilibration between the RR and ETCO₂) were 7.53 and 0.467 min^-1, respectively. The C_e50 and ETCO₂ data fit to a sigmoid E_max model. In conclusion, the RR required to get the target ETCO₂ was much lower in children patients taking antiepileptic drugs than that of non-user children patients during the general anesthesia.

Structural basis for antiarrhythmic drug interactions with the human cardiac sodium channel

Voltage-gated sodium channels play a central role in cellular excitability and are key targets for drug development. Recent breakthroughs in high-resolution cryo-electron microscopy protein structure determination, Rosetta computational protein structure modeling, and multimicrosecond molecular dynamics simulations are empowering advances in structural biology to study the atomistic details of channel−drug interactions. We used Rosetta structural computational modeling and molecular dynamics simulations to study the interactions of antiarrhythmic and local anesthetic drugs with cardiac sodium channel. Our results provide crucial atomic-scale mechanistic insights into the channel–drug interactions, necessary for the rational design of novel modulators of the human cardiac sodium channel to be used for the treatment of cardiac arrhythmias.

The human voltage-gated sodium channel, hNa_V1.5, is responsible for the rapid upstroke of the cardiac action potential and is target for antiarrhythmic therapy. Despite the clinical relevance of hNa_V1.5-targeting drugs, structure-based molecular mechanisms of promising or problematic drugs have not been investigated at atomic scale to inform drug design. Here, we used Rosetta structural modeling and docking as well as molecular dynamics simulations to study the interactions of antiarrhythmic and local anesthetic drugs with hNa_V1.5. These calculations revealed several key drug binding sites formed within the pore lumen that can simultaneously accommodate up to two drug molecules. Molecular dynamics simulations identified a hydrophilic access pathway through the intracellular gate and a hydrophobic access pathway through a fenestration between DIII and DIV. Our results advance the understanding of molecular mechanisms of antiarrhythmic and local anesthetic drug interactions with hNa_V1.5 and will be useful for rational design of novel therapeutics.

Fenestrations control resting-state block of a voltage-gated sodium channel

Potency of drug action is usually determined by binding to a specific receptor site on target proteins. In contrast to this conventional paradigm, we show here that potency of local anesthetics (LAs) and antiarrhythmic drugs (AADs) that block sodium channels is controlled by fenestrations that allow drug access to the receptor site directly from the membrane phase. Voltage-gated sodium channels initiate action potentials in nerve and cardiac muscle, where their hyperactivity causes pain and cardiac arrhythmia, respectively. LAs and AADs selectively block sodium channels in rapidly firing nerve and muscle cells to relieve these conditions. The structure of the ancestral bacterial sodium channel Na_VAb, which is also blocked by LAs and AADs, revealed fenestrations connecting the lipid phase of the membrane to the central cavity of the pore. We cocrystallized lidocaine and flecainide with Na_vAb, which revealed strong drug-dependent electron density in the central cavity of the pore. Mutation of the contact residue T206 greatly reduced drug potency, confirming this site as the receptor for LAs and AADs. Strikingly, mutations of the fenestration cap residue F203 changed fenestration size and had graded effects on resting-state block by flecainide, lidocaine, and benzocaine, the potencies of which were altered from 51- to 2.6-fold in order of their molecular size. These results show that conserved fenestrations in the pores of sodium channels are crucial pharmacologically and determine the level of resting-state block by widely used drugs. Fine-tuning drug access through fenestrations provides an unexpected avenue for structure-based design of ion-channel-blocking drugs.

Molecular Dynamics Simulations of Membrane Proteins: An Overview

Simulations of membrane proteins have been rising in popularity in the past decade. Advancements in technology and force fields made it possible to simulate behavior of membrane proteins. Membrane protein simulations can now be used as supporting evidence for experimental findings, for elucidating protein mechanisms, and validating protein crystal structures. Unrelated to experimental data, these simulations can also serve to investigate larger scale processes like protein sorting, protein-membrane interactions, and more. In this review, the history as well as the state-of-the-art methodologies in membrane protein simulations will be summarized. An emphasis will be put on how to set up the system and on the current models for the different components of the simulation system. An overview of the available tools for membrane protein simulation will be given, and current limitations and prospects will also be discussed.

Role of the Interaction Motif in Maintaining the Open Gate of an Open Sodium Channel

Voltage-gated sodium channels undergo transitions between open, closed, and inactivated states, enabling regulation of the translocation of sodium ions across membranes. A recently published crystal structure of the full-length prokaryotic NavMs crystal structure in the activated open conformation has revealed the presence of a novel motif consisting of an extensive network of salt bridges involving residues in the voltage sensor, S4-S5 linker, pore, and C-terminal domains. This motif has been proposed to be responsible for maintaining an open conformation that enables ion translocation through the channel. In this study, we have used long-time molecular dynamics calculations without artificial restraints to demonstrate that the interaction network of full-length NavMs indeed prevents a rapid collapse and closure of the gate, in marked difference to earlier studies of the pore-only construct in which the gate had to be restrained to remain open. Interestingly, a frequently discussed "hydrophobic gating" mechanism at nanoscopic level is also observed in our simulations, in which the discontinuous water wire close to the gate region leads to an energetic barrier for ion conduction. In addition, we demonstrate the effects of in silico mutations of several of the key residues in the motif on the open channel's stability and functioning, correlating them with existing functional studies on this channel and homologous disease-associated mutations in human sodium channels; we also examine the effects of truncating/removing the voltage sensor and C-terminal domains in maintaining an open gate.

Propofol inhibits the voltage-gated sodium channel NaChBac at multiple sites

Voltage-gated sodium (NaV) channels are important targets of general anesthetics, including the intravenous anesthetic propofol. Electrophysiology studies on the prokaryotic NaV channel NaChBac have demonstrated that propofol promotes channel activation and accelerates activation-coupled inactivation, but the molecular mechanisms of these effects are unclear. Here, guided by computational docking and molecular dynamics simulations, we predict several propofol-binding sites in NaChBac. We then strategically place small fluorinated probes at these putative binding sites and experimentally quantify the interaction strengths with a fluorinated propofol analogue, 4-fluoropropofol. In vitro and in vivo measurements show that 4-fluoropropofol and propofol have similar effects on NaChBac function and nearly identical anesthetizing effects on tadpole mobility. Using quantitative analysis by 19F-NMR saturation transfer difference spectroscopy, we reveal strong intermolecular cross-relaxation rate constants between 4-fluoropropofol and four different regions of NaChBac, including the activation gate and selectivity filter in the pore, the voltage sensing domain, and the S4-S5 linker. Unlike volatile anesthetics, 4-fluoropropofol does not bind to the extracellular interface of the pore domain. Collectively, our results show that propofol inhibits NaChBac at multiple sites, likely with distinct modes of action. This study provides a molecular basis for understanding the net inhibitory action of propofol on NaV channels.

Location and Character of Volatile General Anesthetics Binding Sites in the Transmembrane Domain of TRPV1

It has been proposed that general anesthesia results from direct multisite interactions with multiple and diverse ion channels in the brain. An understanding of the mechanisms by which general anesthetics modulate ion channels is essential to clarify their underlying behavior and their role in reversible immobilization and amnesia. Despite the fact that volatile general anesthetics are drugs that primarily induce insensitivity to pain, they have been reported to sensitize and active the vanilloid-1 receptor, TRPV1, which is known to mediate the response of the nervous system to certain harmful stimuli and which plays a crucial role in the pain pathway. Currently, the mechanism of action of anesthetics is unknown and the precise molecular sites of interaction have not been identified. Here, using ∼2.5 μs of classical molecular dynamics simulations and metadynamics, we explore these enigmas. Binding sites are identified and the strength of the association is further characterized using alchemical free-energy calculations. Anesthetic binding/unbinding proceeds primarily through a membrane-embedded pathway, and subsequently, a complex scenario is established involving multiple binding sites featuring single or multiple occupancy states of two small volatile drugs. One of the five anesthetic binding sites reported was previously identified experimentally, and another one, importantly, is identical to that of capsaicin, one of the chemical stimuli that activate TRPV1. However, in contrast to capsaicin, isoflurane and chloroform binding free-energies render modest to no association compared to capsaicin, suggesting a different activation mechanism. Uncovering chloroform and isoflurane modulatory sites will further our understanding of the TRPV1 molecular machinery and open the possibility of developing site-specific drugs.

Predicting Structural Details of the Sodium Channel Pore Basing on Animal Toxin Studies

Eukaryotic voltage-gated sodium channels play key roles in physiology and are targets for many toxins and medically important drugs. Physiology, pharmacology, and general architecture of the channels has long been the subject of intensive research in academia and industry. In particular, animal toxins such as tetrodotoxin, saxitoxin, and conotoxins have been used as molecular probes of the channel structure. More recently, X-ray structures of potassium and prokaryotic sodium channels allowed elaborating models of the toxin-channel complexes that integrated data from biophysical, electrophysiological, and mutational studies. Atomic level cryo-EM structures of eukaryotic sodium channels, which became available in 2017, show that the selectivity filter structure and other important features of the pore domain have been correctly predicted. This validates further employments of toxins and other small molecules as sensitive probes of fine structural details of ion channels.

Electrostatic lock in the transport cycle of the multidrug resistance transporter EmrE

EmrE is a small, homodimeric membrane transporter that exploits the established electrochemical proton gradient across the Escherichia coli inner membrane to export toxic polyaromatic cations, prototypical of the wider small-multidrug resistance transporter family. While prior studies have established many fundamental aspects of the specificity and rate of substrate transport in EmrE, low resolution of available structures has hampered identification of the transport coupling mechanism. Here we present a complete, refined atomic structure of EmrE optimized against available cryo-electron microscopy (cryo-EM) data to delineate the critical interactions by which EmrE regulates its conformation during the transport process. With the model, we conduct molecular dynamics simulations of the transporter in explicit membranes to probe EmrE dynamics under different substrate loading and conformational states, representing different intermediates in the transport cycle. The refined model is stable under extended simulation. The water dynamics in simulation indicate that the hydrogen-bonding networks around a pair of solvent-exposed glutamate residues (E14) depend on the loading state of EmrE. One specific hydrogen bond from a tyrosine (Y60) on one monomer to a glutamate (E14) on the opposite monomer is especially critical, as it locks the protein conformation when the glutamate is deprotonated. The hydrogen bond provided by Y60 lowers the [Formula: see text] of one glutamate relative to the other, suggesting both glutamates should be protonated for the hydrogen bond to break and a substrate-free transition to take place. These findings establish the molecular mechanism for the coupling between proton transfer reactions and protein conformation in this proton-coupled secondary transporter.

Searching for New Leads To Treat Epilepsy: Target-Based Virtual Screening for the Discovery of Anticonvulsant Agents

The purpose of this investigation is to contribute to the development of new anticonvulsant drugs to treat patients with refractory epilepsy. We applied a virtual screening protocol that involved the search into molecular databases of new compounds and known drugs to find small molecules that interact with the open conformation of the Nav1.2 pore. As the 3D structure of human Nav1.2 is not available, we first assembled 3D models of the target, in closed and open conformations. After the virtual screening, the resulting candidates were submitted to a second virtual filter, to find compounds with better chances of being effective for the treatment of P-glycoprotein (P-gp) mediated resistant epilepsy. Again, we built a model of the 3D structure of human P-gp, and we validated the docking methodology selected to propose the best candidates, which were experimentally tested on Nav1.2 channels by patch clamp techniques and in vivo by the maximal electroshock seizure (MES) test. Patch clamp studies allowed us to corroborate that our candidates, drugs used for the treatment of other pathologies like Ciprofloxacin, Losartan, and Valsartan, exhibit inhibitory effects on Nav1.2 channel activity. Additionally, a compound synthesized in our lab, N, N'-diphenethylsulfamide, interacts with the target and also triggers significant Na1.2 channel inhibitory action. Finally, in vivo studies confirmed the anticonvulsant action of Valsartan, Ciprofloxacin, and N, N'-diphenethylsulfamide.

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.

Polyamine Modulation of Anticonvulsant Drug Response: A Potential Mechanism Contributing to Pharmacoresistance in Chronic Epilepsy

Despite the development of numerous novel anticonvulsant drugs, ∼30% of epilepsy patients remain refractory to antiepileptic drugs (AEDs). Many established and novel AEDs reduce hyperexcitability via voltage- and use-dependent inhibition of voltage-gated Na⁺ channels. For the widely used anticonvulsant carbamazepine (CBZ), use-dependent block of Na⁺ channels is significantly reduced both in experimental and human epilepsy. However, the molecular underpinnings of this potential cellular mechanism for pharmacoresistance have remained enigmatic.Here, we describe the mechanism that leads to the emergence of CBZ-resistant Na⁺ channels. We focused on the endogenous polyamine system, which powerfully modulates Na⁺ channels in a use-dependent manner. We had shown previously that the intracellular polyamine spermine is reduced in chronic epilepsy, resulting in increased persistent Na⁺ currents. Because spermine and CBZ both bind use-dependently in spatial proximity within the Na⁺ channel pore, we hypothesized that spermine loss might also be related to diminished CBZ response. Using the pilocarpine model of refractory epilepsy in male rats and whole-cell patch-clamp recordings, we first replicated the reduction of use-dependent block by CBZ in chronically epileptic animals. We then substituted intracellular spermine via the patch pipette in different concentrations. Under these conditions, we found that exogenous spermine significantly rescues use-dependent block of Na⁺ channels by CBZ. These findings indicate that an unexpected modulatory mechanism, depletion of intracellular polyamines, leads both to increased persistent Na⁺ currents and to diminished CBZ sensitivity of Na⁺ channels. These findings could lead to novel strategies for overcoming pharmacoresistant epilepsy that target the polyamine system.SIGNIFICANCE STATEMENT Pharmacoresistant epilepsy affects ∼18 million people worldwide, and intense efforts have therefore been undertaken to uncover the underlying molecular and cellular mechanisms. One of the key known candidate mechanisms of pharmacoresistance has been a loss of use-dependent Na⁺ channel block by the anticonvulsant carbamazepine (CBZ), both in human and experimental epilepsies. Despite intense scrutiny, the molecular mechanisms underlying this phenomenon have not been elucidated. We now show that a loss of intracellular spermine in chronic epilepsy is a major causative factor leading to the development of CBZ-resistant Na⁺ currents. This finding can be exploited both for the screening of anticonvulsants in expression systems, and for novel strategies to overcome pharmacoresistance that target the polyamine system.

Protonation state of inhibitors determines interaction sites within voltage-gated sodium channels

Voltage-gated sodium channels are essential for carrying electrical signals throughout the body, and mutations in these proteins are responsible for a variety of disorders, including epilepsy and pain syndromes. As such, they are the target of a number of drugs used for reducing pain or combatting arrhythmias and seizures. However, these drugs affect all sodium channel subtypes found in the body. Designing compounds to target select sodium channel subtypes will provide a new therapeutic pathway and would maximize treatment efficacy while minimizing side effects. Here, we examine the binding preferences of nine compounds known to be sodium channel pore blockers in molecular dynamics simulations. We use the approach of replica exchange solute tempering (REST) to gain a more complete understanding of the inhibitors' behavior inside the pore of NavMs, a bacterial sodium channel, and NavPas, a eukaryotic sodium channel. Using these simulations, we are able to show that both charged and neutral compounds partition into the bilayer, but neutral forms more readily cross it. We show that there are two possible binding sites for the compounds: (i) a site on helix 6, which has been previously determined by many experimental and computational studies, and (ii) an additional site, occupied by protonated compounds in which the positively charged part of the drug is attracted into the selectivity filter. Distinguishing distinct binding poses for neutral and charged compounds is essential for understanding the nature of pore block and will aid the design of subtype-selective sodium channel inhibitors.

Digging into Lipid Membrane Permeation for Cardiac Ion Channel Blocker d-Sotalol with All-Atom Simulations

Interactions of drug molecules with lipid membranes play crucial role in their accessibility of cellular targets and can be an important predictor of their therapeutic and safety profiles. Very little is known about spatial localization of various drugs in the lipid bilayers, their active form (ionization state) or translocation rates and therefore potency to bind to different sites in membrane proteins. All-atom molecular simulations may help to map drug partitioning kinetics and thermodynamics, thus providing in-depth assessment of drug lipophilicity. As a proof of principle, we evaluated extensively lipid membrane partitioning of d-sotalol, well-known blocker of a cardiac potassium channel K_v11.1 encoded by the hERG gene, with reported substantial proclivity for arrhythmogenesis. We developed the positively charged (cationic) and neutral d-sotalol models, compatible with the biomolecular CHARMM force field, and subjected them to all-atom molecular dynamics (MD) simulations of drug partitioning through hydrated lipid membranes, aiming to elucidate thermodynamics and kinetics of their translocation and thus putative propensities for hydrophobic and aqueous hERG access. We found that only a neutral form of d-sotalol accumulates in the membrane interior and can move across the bilayer within millisecond time scale, and can be relevant to a lipophilic channel access. The computed water-membrane partitioning coefficient for this form is in good agreement with experiment. There is a large energetic barrier for a cationic form of the drug, dominant in water, to cross the membrane, resulting in slow membrane translocation kinetics. However, this form of the drug can be important for an aqueous access pathway through the intracellular gate of hERG. This route will likely occur after a neutral form of a drug crosses the membrane and subsequently re-protonates. Our study serves to demonstrate a first step toward a framework for multi-scale in silico safety pharmacology, and identifies some of the challenges that lie therein.