Evaluating inositol phospholipid interactions with inward rectifier potassium channels and characterising their role in disease

Molecular Dynamics Simulation of Kir6.2 Variants Reveals Potential Association with Diabetes Mellitus

Diabetes mellitus (DM) represents a problem for the healthcare system worldwide. DM has very serious complications such as blindness, kidney failure, and cardiovascular disease. In addition to the very bad socioeconomic impacts, it influences patients and their families and communities. The global costs of DM and its complications are huge and expected to rise by the year 2030. DM is caused by genetic and environmental risk factors. Genetic testing will aid in early diagnosis and identification of susceptible individuals or populations using ATP-sensitive potassium (KATP) channels present in different tissues such as the pancreas, myocardium, myocytes, and nervous tissues. The channels respond to different concentrations of blood sugar, stimulation by hormones, or ischemic conditions. In pancreatic cells, they regulate the secretion of insulin and glucagon. Mutations in the KCNJ11 gene that encodes the Kir6.2 protein (a major constituent of KATP channels) were reported to be associated with Type 2 DM, neonatal diabetes mellitus (NDM), and maturity-onset diabetes of the young (MODY). Kir6.2 harbors binding sites for ATP and phosphatidylinositol 4,5-diphosphate (PIP2). The ATP inhibits the KATP channel, while the (PIP2) activates it. A Kir6.2 mutation at tyrosine330 (Y330) was demonstrated to reduce ATP inhibition and predisposes to NDM. In this study, we examined the effect of mutations on the Kir6.2 structure using bioinformatics tools and molecular dynamic simulations (SIFT, PolyPhen, SNAP2, PANTHER, PhD&SNP, SNP&Go, I-Mutant, MuPro, MutPred, ConSurf, HOPE, and GROMACS). Our results indicated that M199R, R201H, R206H, and Y330H mutations influence Kir6.2 structure and function and therefore may cause DM. We conclude that MD simulations are useful techniques to predict the effects of mutations on protein structure. In addition, the M199R, R201H, R206H, and Y330H variant in the Kir6.2 protein may be associated with DM. These results require further verification in protein-protein interactions, Kir6.2 function, and case-control studies.

Electro-metabolic signaling

Precise matching of energy substrate delivery to local metabolic needs is essential for the health and function of all tissues. Here, we outline a mechanistic framework for understanding this critical process, which we refer to as electro-metabolic signaling (EMS). All tissues exhibit changes in metabolism over varying spatiotemporal scales and have widely varying energetic needs and reserves. We propose that across tissues, common signatures of elevated metabolism or increases in energy substrate usage that exceed key local thresholds rapidly engage mechanisms that generate hyperpolarizing electrical signals in capillaries that then relax contractile elements throughout the vasculature to quickly adjust blood flow to meet changing needs. The attendant increase in energy substrate delivery serves to meet local metabolic requirements and thus avoids a mismatch in supply and demand and prevents metabolic stress. We discuss in detail key examples of EMS that our laboratories have discovered in the brain and the heart, and we outline potential further EMS mechanisms operating in tissues such as skeletal muscle, pancreas, and kidney. We suggest that the energy imbalance evoked by EMS uncoupling may be central to cellular dysfunction from which the hallmarks of aging and metabolic diseases emerge and may lead to generalized organ failure states-such as diverse flavors of heart failure and dementia. Understanding and manipulating EMS may be key to preventing or reversing these dysfunctions.

Alchemical Free Energy Calculations on Membrane-Associated Proteins

Membrane proteins have diverse functions within cells and are well-established drug targets. The advances in membrane protein structural biology have revealed drug and lipid binding sites on membrane proteins, while computational methods such as molecular simulations can resolve the thermodynamic basis of these interactions. Particularly, alchemical free energy calculations have shown promise in the calculation of reliable and reproducible binding free energies of protein-ligand and protein-lipid complexes in membrane-associated systems. In this review, we present an overview of representative alchemical free energy studies on G-protein-coupled receptors, ion channels, transporters as well as protein-lipid interactions, with emphasis on best practices and critical aspects of running these simulations. Additionally, we analyze challenges and successes when running alchemical free energy calculations on membrane-associated proteins. Finally, we highlight the value of alchemical free energy calculations calculations in drug discovery and their applicability in the pharmaceutical industry.

Rapid Characterization of the Functional and Pharmacological Consequences of Cantú Syndrome KATP Channel Mutations in Intact Cells

Gain-of-function of KATP channels, resulting from mutations in either KCNJ8 (encoding inward rectifier sub-family 6 [Kir6.1]) or ABCC9 (encoding sulphonylurea receptor [SUR2]), cause Cantú syndrome (CS), a channelopathy characterized by excess hair growth, coarse facial appearance, cardiomegaly, and lymphedema. Here, we established a pipeline for rapid analysis of CS mutation consequences in Landing pad HEK 293 cell lines stably expressing wild type (WT) and mutant human Kir6.1 and SUR2B. Thallium-influx and cell membrane potential, reported by fluorescent Tl-sensitive Fluozin-2 and voltage-sensitive bis-(1,3-dibutylbarbituric acid)trimethine oxonol (DiBAC4(3)) dyes, respectively, were used to assess channel activity. In the Tl-influx assay, CS-associated Kir6.1 mutations increased sensitivity to the ATP-sensitive potassium (KATP) channel activator, pinacidil, but there was strikingly little effect of pinacidil for any SUR2B mutations, reflecting unexpected differences in the molecular mechanisms of Kir6.1 versus SUR2B mutations. Compared with the Tl-influx assay, the DiBAC4(3) assay presents more significant signal changes in response to subtle KATP channel activity changes, and all CS mutants (both Kir6.1 and SUR2B), but not WT channels, caused marked hyperpolarization, demonstrating that all mutants were activated under ambient conditions in intact cells. Most SUR2 CS mutations were markedly inhibited by <100 nM glibenclamide, but sensitivity to inhibition by glibenclamide, repaglinide, and PNU37883A was markedly reduced for Kir6.1 CS mutations. Understanding functional consequences of mutations can help with disease diagnosis and treatment. The analysis pipeline we have developed has the potential to rapidly identify mutational consequences, aiding future CS diagnosis, drug discovery, and individualization of treatment. SIGNIFICANCE STATEMENT: We have developed new fluorescence-based assays of channel activities and drug sensitivities of Cantú syndrome (CS) mutations in human Kir6.1/SUR2B-dependent KATP channels, showing that Kir6.1 mutations increase sensitivity to potassium channel openers, while SUR2B mutations markedly reduce K channel opener (KCO) sensitivity. However, both Kir6.1 and SUR2B CS mutations are both more hyperpolarized than WT cells under basal conditions, confirming pathophysiologically relevant gain-of-function, validating DiBAC4(3) fluorescence to characterize hyperpolarization induced by KATP channel activity under basal, non KCO-activated conditions.

Nutrigenomics of inward rectifier potassium channels

Inwardly rectifying potassium (Kir) channels play a key role in maintaining the resting membrane potential and supporting potassium homeostasis. There are many variants of Kir channels, which are usually tetramers in which the main subunit has two trans-membrane helices attached to two N- and C-terminal cytoplasmic tails with a pore-forming loop in between that contains the selectivity filter. These channels have domains that are strongly modulated by molecules present in nutrients found in different diets, such as phosphoinositols, polyamines and Mg²⁺. These molecules can impact these channels directly or indirectly, either allosterically by modulation of enzymes or via the regulation of channel expression. A particular type of these channels is coupled to cell metabolism and inhibited by ATP (KATP channels, essential for insulin release and for the pathogenesis of metabolic diseases like diabetes mellitus). Genomic changes in Kir channels have a significant impact on metabolism, such as conditioning the nutrients and electrolytes that an individual can take. Thus, the nutrigenomics of ion channels is an important emerging field in which we are attempting to understand how nutrients and diets can affect the activity and expression of ion channels and how genomic changes in such channels may be the basis for pathological conditions that limit nutrition and electrolyte intake. In this contribution we briefly review Kir channels, discuss their nutrigenomics, characterize how different components in the diet affect their function and expression, and suggest how their genomic changes lead to pathological phenotypes that affect diet and electrolyte intake.

Mutational Insight into Allosteric Regulation of Kir Channel Activity

Potassium (K⁺) channels are regulated in part by allosteric communication between the helical bundle crossing, or inner gate, and the selectivity filter, or outer gate. This network is triggered by gating stimuli. In concert, there is an allosteric network which is a conjugated set of interactions which correlate long-range structural rearrangements necessary for channel function. Inward-rectifier K⁺ (Kir) channels favor inward K⁺ conductance, are ligand-gated, and help establish resting membrane potentials. KirBac1.1 is a bacterial Kir (KirBac) channel homologous to human Kir (hKir) channels. Additionally, KirBac1.1 is gated by the anionic phospholipid ligand phosphatidylglycerol (PG). In this study, we use site-directed mutagenesis to investigate residues involved in the KirBac1.1 gating mechanism and allosteric network we previously proposed using detailed solid-state NMR (SSNMR) measurements. Using fluorescence-based K⁺ and sodium (Na⁺) flux assays, we identified channel mutants with impaired function that do not alter selectivity of the channel. In tandem, we performed coarse grain molecular dynamics simulations, observing changes in PG-KirBac1.1 interactions correlated with mutant channel activity and contacts between the two transmembrane helices and pore helix tied to this behavior. Lipid affinity is closely tied to the proximity of two tryptophan residues on neighboring subunits which lure anionic lipids to a cationic pocket formed by a cluster of arginine residues. Thus, these simulations establish a structural and functional basis for the role of each mutated site in the proposed allosteric network. The experimental and simulated data provide insight into key functional residues involved in gating and lipid allostery of K⁺ channels. Our findings also have direct implications on the physiology of hKir channels due to conservation of many of the residues identified in this work from KirBac1.1.

The dynamic interplay of PIP2 and ATP in the regulation of the KATP channel

ATP-sensitive potassium (KATP ) channels couple the intracellular ATP concentration to insulin secretion. KATP channel activity is inhibited by ATP binding to the Kir6.2 tetramer and activated by phosphatidylinositol 4,5-bisphosphate (PIP2 ). Here, we use molecular dynamics simulation, electrophysiology and fluorescence spectroscopy to show that ATP and PIP2 occupy different binding pockets that share a single amino acid residue, K39. When both ligands are present, simulations suggest that K39 shows a greater preference to co-ordinate with PIP2 than with ATP. They also predict that a neonatal diabetes mutation at K39 (K39R) increases the number of hydrogen bonds formed between K39 and PIP2 , potentially accounting for the reduced ATP inhibition observed in electrophysiological experiments. Our work suggests that PIP2 and ATP interact allosterically to regulate KATP channel activity. KEY POINTS: The KATP channel is activated by the binding of phosphatidylinositol 4,5-bisphosphate (PIP2 ) lipids and inactivated by the binding of ATP. K39 has the potential to bind to both PIP2 and ATP. A mutation to this residue (K39R) results in neonatal diabetes. This study uses patch-clamp fluorometry, electrophysiology and molecular dynamics simulation. We show that PIP2 competes with ATP for K39, and this reduces channel inhibition by ATP. We show that K39R increases channel affinity to PIP2 by increasing the number of hydrogen bonds with PIP2 , when compared with the wild-type K39. This therefore decreases KATP channel inhibition by ATP.

Transition between conformational states of the TREK-1 K2P channel promoted by interaction with PIP2

Members of the TREK family of two-pore domain (K2P) potassium channels are highly sensitive to regulation by membrane lipids, including phosphatidylinositol-4,5-bisphosphate (PIP₂). Previous studies have demonstrated that PIP₂ increases TREK1 channel activity, however, the mechanistic understanding of the conformational transitions induced by PIP₂ remain unclear. Here, we used coarse-grained molecular dynamics (CG-MD) and atomistic MD simulations to model the PIP₂ binding site on both the up and down state conformations of TREK-1. We also calculated the free energy of PIP₂ binding relative to other anionic phospholipids in both conformational states using potential of mean force (PMF) and free energy perturbation (FEP) calculations. Our results identify state-dependent binding of PIP₂ to sites involving the proximal C-terminus and we show that PIP₂ promotes a conformational transition from a down state towards an intermediate that resembles the up state. These results are consistent with functional data for PIP₂ regulation and together provide evidence for a structural mechanism of TREK-1 channel activation by phosphoinositides.

Specific interactions of peripheral membrane proteins with lipids: what can molecular simulations show us?

Peripheral membrane proteins (PMPs) can reversibly and specifically bind to biological membranes to carry out functions such as cell signalling, enzymatic activity, or membrane remodelling. Structures of these proteins and of their lipid-binding domains are typically solved in a soluble form, sometimes with a lipid or lipid headgroup at the binding site. To provide a detailed molecular view of PMP interactions with the membrane, computational methods such as molecular dynamics (MD) simulations can be applied. Here, we outline recent attempts to characterise these binding interactions, focusing on both intracellular proteins, such as phosphatidylinositol phosphate (PIP)-binding domains, and extracellular proteins such as glycolipid-binding bacterial exotoxins. We compare methods used to identify and analyse lipid-binding sites from simulation data and highlight recent work characterising the energetics of these interactions using free energy calculations. We describe how improvements in methodologies and computing power will help MD simulations to continue to contribute to this field in the future.

Entropy in the Molecular Recognition of Membrane Protein-Lipid Interactions

Understanding the molecular driving forces that underlie membrane protein-lipid interactions requires the characterization of their binding thermodynamics. Here, we employ variable-temperature native mass spectrometry to determine the thermodynamics of lipid binding events to the human G-protein-gated inward rectifier potassium channel, Kir3.2. The channel displays distinct thermodynamic strategies to engage phosphatidylinositol (PI) and phosphorylated forms thereof. The addition of a 4'-phosphate to PI results in an increase in favorable entropy. PI with two or more phosphates exhibits more complex binding, where lipids appear to bind two nonidentical sites on Kir3.2. Remarkably, the interaction of 4,5-bisphosphate PI with Kir3.2 is solely driven by a large, favorable change in entropy. Installment of a 3'-phosphate to PI(4,5)P₂ results in an altered thermodynamic strategy. The acyl chain of the lipid has a marked impact on binding thermodynamics and, in some cases, enthalpy becomes favorable.

Relative Affinities of Protein-Cholesterol Interactions from Equilibrium Molecular Dynamics Simulations

Specific interactions of lipids with membrane proteins contribute to protein stability and function. Multiple lipid interactions surrounding a membrane protein are often identified in molecular dynamics (MD) simulations and are, increasingly, resolved in cryo-electron microscopy (cryo-EM) densities. Determining the relative importance of specific interaction sites is aided by determination of lipid binding affinities using experimental or simulation methods. Here, we develop a method for determining protein-lipid binding affinities from equilibrium coarse-grained MD simulations using binding saturation curves, designed to mimic experimental protocols. We apply this method to directly obtain affinities for cholesterol binding to multiple sites on a range of membrane proteins and compare our results with free energies obtained from density-based equilibrium methods and with potential of mean force calculations, getting good agreement with respect to the ranking of affinities for different sites. Thus, our binding saturation method provides a robust, high-throughput alternative for determining the relative consequence of individual sites seen in, e.g., cryo-EM derived membrane protein structures surrounded by an array of ancillary lipid densities.

From Bench to Biomolecular Simulation: Phospholipid Modulation of Potassium Channels

Potassium (K+) ion channels are crucial in numerous cellular processes as they hyperpolarise a cell through K+ conductance, returning a cell to its resting potential. K+ channel mutations result in multiple clinical complications such as arrhythmia, neonatal diabetes and migraines. Since 1995, the regulation of K+ channels by phospholipids has been heavily studied using a range of interdisciplinary methods such as cellular electrophysiology, structural biology and computational modelling. As a result, K+ channels are model proteins for the analysis of protein-lipid interactions. In this review, we will focus on the roles of lipids in the regulation of K+ channels, and how atomic-level structures, along with experimental techniques and molecular simulations, have helped guide our understanding of the importance of phospholipid interactions.

Simulating PIP2-Induced Gating Transitions in Kir6.2 Channels

ATP-sensitive potassium (KATP) channels consist of an inwardly rectifying K+ channel (Kir6.2) pore, to which four ATP-sensitive sulfonylurea receptor (SUR) domains are attached, thereby coupling K+ permeation directly to the metabolic state of the cell. Dysfunction is linked to neonatal diabetes and other diseases. K+ flux through these channels is controlled by conformational changes in the helix bundle region, which acts as a physical barrier for K+ permeation. In addition, the G-loop, located in the cytoplasmic domain, and the selectivity filter might contribute to gating, as suggested by different disease-causing mutations. Gating of Kir channels is regulated by different ligands, like Gβγ, H+, Na+, adenosine nucleotides, and the signaling lipid phosphatidyl-inositol 4,5-bisphosphate (PIP2), which is an essential activator for all eukaryotic Kir family members. Although molecular determinants of PIP2 activation of KATP channels have been investigated in functional studies, structural information of the binding site is still lacking as PIP2 could not be resolved in Kir6.2 cryo-EM structures. In this study, we used Molecular Dynamics (MD) simulations to examine the dynamics of residues associated with gating in Kir6.2. By combining this structural information with functional data, we investigated the mechanism underlying Kir6.2 channel regulation by PIP2.

Identification and assessment of cardiolipin interactions with E. coli inner membrane proteins

Integral membrane proteins are localized and/or regulated by lipids present in the surrounding bilayer. While bacteria have relatively simple membranes, there is ample evidence that many bacterial proteins bind to specific lipids, especially the anionic lipid cardiolipin. Here, we apply molecular dynamics simulations to assess lipid binding to 42 different Escherichia coli inner membrane proteins. Our data reveal an asymmetry between the membrane leaflets, with increased anionic lipid binding to the inner leaflet regions of the proteins, particularly for cardiolipin. From our simulations, we identify >700 independent cardiolipin binding sites, allowing us to identify the molecular basis of a prototypical cardiolipin binding site, which we validate against structures of bacterial proteins bound to cardiolipin. This allows us to construct a set of metrics for defining a high-affinity cardiolipin binding site on bacterial membrane proteins, paving the way for a heuristic approach to defining other protein-lipid interactions.