Lateral Fenestrations in K(+)-Channels Explored Using Molecular Dynamics Simulations

GABAB Receptor: Structure, Biological Functions, and Therapy for Diseases

Gamma-aminobutyric acid (GABA) B receptors (GABABRs) that acts slowly and maintains the inhibitory tone are versatile regulators in the complex nervous behaviors and their involvement in various neuropsychiatric disorders, such as anxiety, epilepsy, pain, drug addiction, and Alzheimer's disease. Additional study advances have implied the crucial roles of GABABRs in regulating feeding-related behaviors, yet their therapeutic potential in addressing the neuropsychiatric disorders, binge eating, and feeding-related disorders remains underutilized. This general review summarized the physiological structure and functions of GABABR, explored the regulation in various psychiatric disorders, feeding behaviors, binge eating, and metabolism disorders, and fully discussed the potential of targeting GABABRs and its regulator-binding sites for the treatment of different psychiatric disorders, binge eating and even obesity. While agonists that directly bind to GABABR1 have some negative side effects, positive allosteric modulators (PAMs) that bind to GABABR2 demonstrate excellent therapeutic efficacy and tolerability and have better safety and therapeutic indexes. Moreover, phosphorylation sites of downstream GABABRs regulators may be novel therapeutic targets for psychiatric disorders, binge eating, and obesity. Further studies, clinical trials in particular, will be essential for confirming the therapeutic value of PAMs and other agents targeting the GABABR pathways in a clinical setting.

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.

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.

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.

The Inhibitory Effect of Magnolol on the Human TWIK1 Channel Is Related to G229 and T225 Sites

TWIK1 (K2P1.1/KCNK1) belongs to the potassium channels of the two-pore domain. Its current is very small and difficult to measure. In this work, we used a 100 mM NH4+ extracellular solution to increase TWIK1 current in its stable cell line expressed in HEK293. Then, the inhibition of magnolol on TWIK1 was observed via a whole-cell patch clamp experiment, and it was found that magnolol had a significant inhibitory effect on TWIK1 (IC50 = 6.21 ± 0.13 μM). By molecular docking and alanine scanning mutagenesis, the IC50 of TWIK1 mutants G229A, T225A, I140A, L223A, and S224A was 20.77 ± 3.20, 21.81 ± 7.93, 10.22 ± 1.07, 9.55 ± 1.62, and 7.43 ± 3.20 μM, respectively. Thus, we conclude that the inhibition of the TWIK1 channel by magnolol is related to G229 and T225 on the P2- pore helix.

Advances in the Understanding of Two-Pore Domain TASK Potassium Channels and Their Potential as Therapeutic Targets

TWIK-related acid-sensitive K⁺ (TASK) channels, including TASK-1, TASK-3, and TASK-5, are important members of the two-pore domain potassium (K_2P) channel family. TASK-5 is not functionally expressed in the recombinant system. TASK channels are very sensitive to changes in extracellular pH and are active during all membrane potential periods. They are similar to other K_2P channels in that they can create and use background-leaked potassium currents to stabilize resting membrane conductance and repolarize the action potential of excitable cells. TASK channels are expressed in both the nervous system and peripheral tissues, including excitable and non-excitable cells, and are widely engaged in pathophysiological phenomena, such as respiratory stimulation, pulmonary hypertension, arrhythmia, aldosterone secretion, cancers, anesthesia, neurological disorders, glucose homeostasis, and visual sensitivity. Therefore, they are important targets for innovative drug development. In this review, we emphasized the recent advances in our understanding of the biophysical properties, gating profiles, and biological roles of TASK channels. Given the different localization ranges and biologically relevant functions of TASK-1 and TASK-3 channels, the development of compounds that selectively target TASK-1 and TASK-3 channels is also summarized based on data reported in the literature.

Chanalyzer: A Computational Geometry Approach for the Analysis of Protein Channel Shape and Dynamics

Morphological analysis of protein channels is a key step for a thorough understanding of their biological function and mechanism. In this respect, molecular dynamics (MD) is a very powerful tool, enabling the description of relevant biological events at the atomic level, which might elude experimental observations, and pointing to the molecular determinants thereof. In this work, we present a computational geometry-based approach for the characterization of the shape and dynamics of biological ion channels or pores to be used in combination with MD trajectories. This technique relies on the earliest works of Edelsbrunner and on the NanoShaper software, which makes use of the alpha shape theory to build the solvent-excluded surface of a molecular system in an aqueous solution. In this framework, a channel can be simply defined as a cavity with two entrances on the opposite sides of a molecule. Morphological characterization, which includes identification of the main axis, the corresponding local radius, and the detailed description of the global shape of the cavity, is integrated with a physico-chemical description of the surface facing the pore lumen. Remarkably, the possible existence or temporary appearance of fenestrations from the channel interior towards the outer lipid matrix is also accounted for. As a test case, we applied the present approach to the analysis of an engineered protein channel, the mechanosensitive channel of large conductance.

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.

Conduction and Gating Properties of the TRAAK Channel from Molecular Dynamics Simulations with Different Force Fields

In recent years, the K2P family of potassium channels has been the subject of intense research activity. Owing to the complex function and regulation of this family of ion channels, it is common practice to complement experimental findings with the atomistic description provided by computational approaches such as molecular dynamics (MD) simulations, especially, in light of the unprecedented timescales accessible at present. However, despite recent substantial improvements, the accuracy of MD simulations is still undermined by the intrinsic limitations of force fields. Here, we systematically assessed the performance of the most popular force fields employed to study ion channels at timescales that are orders of magnitude greater than the ones accessible when these energy functions were first developed. Using 32 μs of trajectories, we investigated the dynamics of a member of the K2P ion channel family, the TRAAK channel, using two established force fields in simulations of biological systems: AMBER and CHARMM. We found that while results are comparable on the nanosecond timescales, significant inconsistencies arise at microsecond timescales.

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.

An Update on the Structure of hERG

The human voltage-sensitive K⁺ channel hERG plays a fundamental role in cardiac action potential repolarization, effectively controlling the QT interval of the electrocardiogram. Inherited loss- or gain-of-function mutations in hERG can result in dangerous “long” (LQTS) or “short” QT syndromes (SQTS), respectively, and the anomalous susceptibility of hERG to block by a diverse range of drugs underlies an acquired LQTS. A recent open channel cryo-EM structure of hERG should greatly advance understanding of the molecular basis of hERG channelopathies and drug-induced LQTS. Here we describe an update of recent research that addresses the nature of the particular gated state of hERG captured in the new structure, and the insight afforded by the structure into the molecular basis for high affinity drug block of hERG, the binding of hERG activators and the molecular basis of hERG's peculiar gating properties. Interpretation of the pharmacology of natural SQTS mutants in the context of the structure is a promising approach to understanding the molecular basis of hERG inactivation, and the structure suggests how voltage-dependent changes in the membrane domain may be transmitted to an extracellular “turret” to effect inactivation through aromatic side chain motifs that are conserved throughout the KCNH family of channels.

Elucidating the Structural Basis of the Intracellular pH Sensing Mechanism of TASK-2 K2P Channels

Two-pore domain potassium (K₂P) channels maintain the cell’s background conductance by stabilizing the resting membrane potential. They assemble as dimers possessing four transmembrane helices in each subunit. K₂P channels were crystallized in “up” and “down” states. The movements of the pore-lining transmembrane TM4 helix produce the aperture or closure of side fenestrations that connect the lipid membrane with the central cavity. When the TM4 helix is in the up-state, the fenestrations are closed, while they are open in the down-state. It is thought that the fenestration states are related to the activity of K₂P channels and the opening of the channels preferentially occurs from the up-state. TASK-2, a member of the TALK subfamily of K₂P channels, is opened by intracellular alkalization leading the deprotonation of the K245 residue at the end of the TM4 helix. This charge neutralization of K245 could be sensitive or coupled to the fenestration state. Here, we describe the relationship between the states of the intramembrane fenestrations and K245 residue in TASK-2 channel. By using molecular modeling and simulations, we show that the protonated state of K245 (K245⁺) favors the open fenestration state and, symmetrically, that the open fenestration state favors the protonated state of the lysine residue. We show that the channel can be completely blocked by Prozac, which is known to induce fenestration opening in TREK-2. K245 protonation and fenestration aperture have an additive effect on the conductance of the channel. The opening of the fenestrations with K245⁺ increases the entrance of lipids into the selectivity filter, blocking the channel. At the same time, the protonation of K245 introduces electrostatic potential energy barriers to ion entrance. We computed the free energy profiles of ion penetration into the channel in different fenestration and K245 protonation states, to show that the effects of the two transformations are summed up, leading to maximum channel blocking. Estimated rates of ion transport are in qualitative agreement with experimental results and support the hypothesis that the most important barrier for ion transport under K245⁺ and open fenestration conditions is the entrance of the ions into the channel.

Interfacial Binding Sites for Cholesterol on TRP Ion Channels

Transient receptor potential (TRP) channels are members of a large family of ion channels located in membranes rich in cholesterol, some of whose functions are affected by the cholesterol content of the membrane. Here, cholesterol binding to TRPs is studied using a docking procedure that allows the transmembrane surface of a TRP to be swept rapidly for potential binding sites at the interfaces on the two sides of the membrane. Cholesterol docking poses determined in this way match 89% of the cholesterol hemisuccinate molecules in published TRP structures when cholesterol hemisuccinate molecules unlikely to represent typical bound cholesterols are excluded. TRPs are tetrameric, with large clefts at the interfaces between subunits; cholesterol poses are located in hollows, largely within these clefts. Comparison of cholesterol poses with phospholipid binding sites suggests that binding of cholesterol to a TRP need not result in displacement of phospholipid molecules from the TRP surface.

TASK Channels Pharmacology: New Challenges in Drug Design

Rational drug design targeting ion channels is an exciting and always evolving research field. New medicinal chemistry strategies are being implemented to explore the wild chemical space and unravel the molecular basis of the ion channels modulators binding mechanisms. TASK channels belong to the two-pore domain potassium channel family and are modulated by extracellular acidosis. They are extensively distributed along the cardiovascular and central nervous systems, and their expression is up- and downregulated in different cancer types, which makes them an attractive therapeutic target. However, TASK channels remain unexplored, and drugs designed to target these channels are poorly selective. Here, we review TASK channels properties and their known blockers and activators, considering the new challenges in ion channels drug design and focusing on the implementation of computational methodologies in the drug discovery process.

Molecular determinants of chemical modulation of two-pore domain potassium channels

The K⁺ ion channels comprising the two-pore domain (K2P) family have specific biophysical roles in generating the critical regulatory K⁺ current. Ion flow through K2P channels and, implicitly, channel regulation is mediated by diverse metabolic and physical inputs such as mechanical stimulation, interaction with lipids or endogenous regulators, intra- or extracellular pH, and phosphorylation, while their function can be finely tuned by chemical compounds. In the latter category, some drug-channel interactions can lead to side effects or have clinical action, while identifying novel chemical modulators of K2Ps is an area of intense research. Due to their cellular and therapeutic importance, much attention was turned to these channels in recent years and several experimental approaches have pinpointed the molecular determinants of K2P chemical modulation. Given their unique structural features and properties, chemical modulators act on K2P channels in multiple and diverse ways. In this review, the particularities of K2P modulation by chemical compounds, such as binding modality, affinity, or position, are identified, synthesized, and linked to structural and functional properties in order to refer to how activators and blockers modify channel function and vice versa, focusing on specificity related to protein structure (and its modification) and cross-linking information among different subfamilies.

Folding and Misfolding of Human Membrane Proteins in Health and Disease: From Single Molecules to Cellular Proteostasis

Advances over the past 25 years have revealed much about how the structural properties of membranes and associated proteins are linked to the thermodynamics and kinetics of membrane protein (MP) folding. At the same time biochemical progress has outlined how cellular proteostasis networks mediate MP folding and manage misfolding in the cell. When combined with results from genomic sequencing, these studies have established paradigms for how MP folding and misfolding are linked to the molecular etiologies of a variety of diseases. This emerging framework has paved the way for the development of a new class of small molecule "pharmacological chaperones" that bind to and stabilize misfolded MP variants, some of which are now in clinical use. In this review, we comprehensively outline current perspectives on the folding and misfolding of integral MPs as well as the mechanisms of cellular MP quality control. Based on these perspectives, we highlight new opportunities for innovations that bridge our molecular understanding of the energetics of MP folding with the nuanced complexity of biological systems. Given the many linkages between MP misfolding and human disease, we also examine some of the exciting opportunities to leverage these advances to address emerging challenges in the development of therapeutics and precision medicine.

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.

The molecular basis for an allosteric inhibition of K+-flux gating in K2P channels

Two-pore-domain potassium (K_2P) channels are key regulators of many physiological and pathophysiological processes and thus emerged as promising drug targets. As for other potassium channels, there is a lack of selective blockers, since drugs preferentially bind to a conserved binding site located in the central cavity. Thus, there is a high medical need to identify novel drug-binding sites outside the conserved lipophilic central cavity and to identify new allosteric mechanisms of channel inhibition. Here, we identified a novel binding site and allosteric inhibition mechanism, disrupting the recently proposed K⁺-flux gating mechanism of K_2P channels, which results in an unusual voltage-dependent block of leak channels belonging to the TASK subfamily. The new binding site and allosteric mechanism of inhibition provide structural and mechanistic insights into the gating of TASK channels and the basis for the drug design of a new class of potent blockers targeting specific types of K_2P channels.

Microscopic view of lipids and their diverse biological functions

Biological membranes and their diverse lipid constituents play key roles in a broad spectrum of cellular and physiological processes. Characterization of membrane-associated phenomena at a microscopic level is therefore essential to our fundamental understanding of such processes. Due to the semi-fluid and dynamic nature of lipid bilayers, and their complex compositions, detailed characterization of biological membranes at an atomic scale has been refractory to experimental approaches. Computational modeling and simulation offer a highly complementary toolset with sufficient spatial and temporal resolutions to fill this gap. Here, we review recent molecular dynamics studies focusing on the diversity of lipid composition of biological membranes, or aiming at the characterization of lipid-protein interaction, with the overall goal of dissecting how lipids impact biological roles of the cellular membranes.

Inhibition of Voltage-Gated K+ Channel Kv1.5 by Antiarrhythmic Drugs

Molecular dynamics simulations are employed to determine the inhibitory mechanisms of three drugs, 5-(4-phenoxybutoxy)psoralen (PAP-1), vernakalant, and flecainide, on the voltage-gated K⁺ channel Kv1.5, a target for the treatment of cardiac arrhythmia. At neutral pH, PAP-1 is neutral, whereas the other two molecules carry one positive charge. We show that PAP-1 forms stable dimers in water, primarily through hydrophobic interactions between aromatic rings. All three molecules bind to the cavity between the Ile508 and Val512 residues from the four subunits of the channel. Once bound, the drug molecules are flexible, with the average root-mean-square fluctuation being between 2 and 3 Å, which is larger than the radius of gyration of a bulky amino acid. The presence of a monomeric PAP-1 causes the permeating K⁺ ion to dehydrate, thereby creating a significant energy barrier. In contrast, vernakalant blocks the ion permeation primarily via an electrostatic mechanism and, therefore, must be in the protonated and charged form to be effective.

Side Fenestrations Provide an "Anchor" for a Stable Binding of A1899 to the Pore of TASK-1 Potassium Channels

A1899 is a potent and selective inhibitor of the two-pore domain potassium (K_2P) channel TASK-1. It was previously reported that A1899 acts as an open-channel blocker and binds to residues of the P1 and P2 regions, the M2 and M4 segments, and the halothane response element. The recently described crystal structures of K_2P channels together with the newly identified side fenestrations indicate that residues relevant for TASK-1 inhibition are not purely facing the central cavity as initially proposed. Accordingly, the TASK-1 binding site and the mechanism of inhibition might need a re-evaluation. We have used TASK-1 homology models based on recently crystallized K_2P channels and molecular dynamics simulation to demonstrate that the highly potent TASK-1 blocker A1899 requires binding to residues located in the side fenestrations. Unexpectedly, most of the previously described residues that interfere with TASK-1 blockade by A1899 project their side chains toward the fenestration lumina, underlining the relevance of these structures for drug binding in K_2P channels. Despite its hydrophobicity, A1899 does not seem to use the fenestrations to gain access to the central cavity from the lipid bilayer. In contrast, binding of A1899 to residues of the side fenestrations might provide a physical "anchor", reflecting an energetically favorable binding mode that after pore occlusion stabilizes the closed state of the channels.