Experimental and computational biophysics to identify vasodilator drugs targeted at TRPV2 using agonists based on the probenecid scaffold.

TRP channels are important pharmacological targets in physiopathology. TRPV2 plays distinct roles in cardiac and neuromuscular function, immunity, and metabolism, and is associated with pathologies like muscular dystrophy and cancer. However, TRPV2 pharmacology is unspecific and scarce at best. Using in silico similarity-based chemoinformatics we obtained a set of 270 potential hits for TRPV2 categorized into families based on chemical nature and similarity. Docking the compounds on available rat TRPV2 structures allowed the clustering of drug families in specific ligand binding sites. Starting from a probenecid docking pose in the piperlongumine binding site and using a Gaussian accelerated molecular dynamics approach we have assigned a putative probenecid binding site. In parallel, we measured the EC50 of 7 probenecid derivatives on TRPV2 expressed in Pichia pastoris using a novel medium-throughput Ca2+ influx assay in yeast membranes together with an unbiased and unsupervised data analysis method. We found that 4-(piperidine-1-sulfonyl)-benzoic acid had a better EC50 than probenecid, which is one of the most specific TRPV2 agonists to date. Exploring the TRPV2-dependent anti-hypertensive potential in vivo, we found that 4-(piperidine-1-sulfonyl)-benzoic acid shows a sex-biased vasodilator effect producing larger vascular relaxations in female mice. Overall, this study expands the pharmacological toolbox for TRPV2, a widely expressed membrane protein and orphan drug target.

N-aryltetrahydroisoquinoline derivatives as HA-CD44 interaction inhibitors: Design, synthesis, computational studies, and antitumor effect.

Hyaluronic acid (HA) plays a crucial role in tumor growth and invasion through its interaction with cluster of differentiation 44 (CD44), a non-kinase transmembrane glycoprotein, among other hyaladherins. CD44 expression is elevated in many solid tumors, and its interaction with HA is associated with cancer and angiogenesis. Despite efforts to inhibit HA-CD44 interaction, there has been limited progress in the development of small molecule inhibitors. As a contribution to this endeavour, we designed and synthesized a series of N-aryltetrahydroisoquinoline derivatives based on existing crystallographic data available for CD44 and HA. Hit 2e was identified within these structures for its antiproliferative effect against two CD44+ cancer cell lines, and two new analogs (5 and 6) were then synthesized and evaluated as CD44-HA inhibitors by applying computational and cell-based CD44 binding studies. Compound 2-(3,4,5-trimethoxybenzyl)-1,2,3,4-tetrahydroisoquinolin-5-ol (5) has an EC50 value of 0.59 µM against MDA-MB-231 cells and is effective to disrupt the integrity of cancer spheroids and reduce the viability of MDA-MB-231 cells in a dose-dependent manner. These results suggest lead 5 as a promising candidate for further investigation in cancer treatment.

Redox regulation of KV7 channels through EF3 hand of calmodulin.

Neuronal KV7 channels, important regulators of cell excitability, are among the most sensitive proteins to reactive oxygen species. The S2S3 linker of the voltage sensor was reported as a site-mediating redox modulation of the channels. Recent structural insights reveal potential interactions between this linker and the Ca2+-binding loop of the third EF-hand of calmodulin (CaM), which embraces an antiparallel fork formed by the C-terminal helices A and B, constituting the calcium responsive domain (CRD). We found that precluding Ca2+ binding to the EF3 hand, but not to EF1, EF2, or EF4 hands, abolishes oxidation-induced enhancement of KV7.4 currents. Monitoring FRET (Fluorescence Resonance Energy Transfer) between helices A and B using purified CRDs tagged with fluorescent proteins, we observed that S2S3 peptides cause a reversal of the signal in the presence of Ca2+ but have no effect in the absence of this cation or if the peptide is oxidized. The capacity of loading EF3 with Ca2+ is essential for this reversal of the FRET signal, whereas the consequences of obliterating Ca2+ binding to EF1, EF2, or EF4 are negligible. Furthermore, we show that EF3 is critical for translating Ca2+ signals to reorient the AB fork. Our data are consistent with the proposal that oxidation of cysteine residues in the S2S3 loop relieves KV7 channels from a constitutive inhibition imposed by interactions between the EF3 hand of CaM which is crucial for this signaling.

Early Steps in C-Type Inactivation of the hERG Potassium Channel.

Fast C-type inactivation confers distinctive functional properties to the hERG potassium channel, and its association to inherited and acquired cardiac arrythmias makes the study of the inactivation mechanism of hERG at the atomic detail of paramount importance. At present, two models have been proposed to describe C-type inactivation in K+-channels. Experimental data and computational work on the bacterial KcsA channel support the hypothesis that C-type inactivation results from a closure of the selectivity filter that sterically impedes ion conduction. Alternatively, recent experimental structures of a mutated Shaker channel revealed a widening of the extracellular portion of the selectivity filter, which might diminish conductance by interfering with the mechanism of ion permeation. Here, we performed molecular dynamics simulations of the wild-type hERG, a non-inactivating mutant (hERG-N629D), and a mutant that inactivates faster than the wild-type channel (hERG-F627Y) to find out which and if any of the two reported C-type inactivation mechanisms applies to hERG. Closure events of the selectivity filter were not observed in any of the simulated trajectories but instead, the extracellular section of the selectivity filter deviated from the canonical conductive structure of potassium channels. The degree of widening of the potassium binding sites at the extracellular entrance of the channel was directly related to the degree of inactivation with hERG-F627Y > wild-type hERG > hERG-N629D. These findings support the hypothesis that C-type inactivation in hERG entails a widening of the extracellular entrance of the channel rather than a closure of the selectivity filter.

A Potential Route of Capsaicin to Its Binding Site in the TRPV1 Ion Channel.

Transient receptor potential (TRP) ion channels are important pharmacological targets because of their role in the perception of pain, and so, understanding their chemical regulation is essential for the development of analgesic drugs. Among the currently known TRP channel chemical agonists, capsaicin, the active compound of chili pepper, is probably the most exhaustively studied. The availability of the three-dimensional structure of the vanilloid receptor 1 (TRPV1) has fueled computational studies revealing the molecular details of capsaicin binding modes. Although this is a significant step, a comprehensible binding mechanism or pathway is invaluable for targeting TRP channels in modern pharmacology. In the present work, free-energy and enhanced sampling techniques have been used to explore a possible membrane-mediated pathway for capsaicin to enter the TRPV1 binding pocket where capsaicin accesses the protein starting at the extracellular milieu through the outer leaflet and into its binding site in the protein. The main states visited along this route have been characterized and include (i) a bound state in agreement with the binding mode "head-down, tail-up" and (ii) an alternative state corresponding to a "head-up, tail-down" binding mode. In agreement with previous reports, binding is mediated by both hydrogen bonds and van der Waals interactions, and residue Y511 is crucial for stabilizing the bound state and during the binding process. Together, these results provide a foundation to further understand TRPV channels, and they could be used to guide therapeutic design of selective inhibitors potentially leading to novel avenues for pharmacological applications targeting the TRPV1 channel.

Selective Anticancer Therapy Based on a HA-CD44 Interaction Inhibitor Loaded on Polymeric Nanoparticles.

Hyaluronic acid (HA), through its interactions with the cluster of differentiation 44 (CD44), acts as a potent modulator of the tumor microenvironment, creating a wide range of extracellular stimuli for tumor growth, angiogenesis, invasion, and metastasis. An innovative antitumor treatment strategy based on the development of a nanodevice for selective release of an inhibitor of the HA-CD44 interaction is presented. Computational analysis was performed to evaluate the interaction of the designed tetrahydroisoquinoline-ketone derivative (JE22) with CD44 binding site. Cell viability, efficiency, and selectivity of drug release under acidic conditions together with CD44 binding capacity, effect on cell migration, and apoptotic activity were successfully evaluated. Remarkably, the conjugation of this CD44 inhibitor to the nanodevice generated a reduction of the dosis required to achieve a significant therapeutic effect.

An Integrated Mass Spectrometry and Molecular Dynamics Simulations Approach Reveals the Spatial Organization Impact of Metal-Binding Sites on the Stability of Metal-Depleted Metallothionein-2 Species.

Mammalian metallothioneins (MTs) are a group of cysteine-rich proteins that bind metal ions in two α- and ß-domains and represent a major cellular Zn(II)/Cu(I) buffering system in the cell. At cellular free Zn(II) concentrations (10-11-10-9 M), MTs do not exist in fully loaded forms with seven Zn(II)-bound ions (Zn7MTs). Instead, MTs exist as partially metal-depleted species (Zn4-6MT) because their Zn(II) binding affinities are on the nano- to picomolar range comparable to the concentrations of cellular Zn(II). The mode of action of MTs remains poorly understood, and thus, the aim of this study is to characterize the mechanism of Zn(II) (un)binding to MTs, the thermodynamic properties of the Zn1-6MT2 species, and their mechanostability properties. To this end, native mass spectrometry (MS) and label-free quantitative bottom-up and top-down MS in combination with steered molecular dynamics simulations, well-tempered metadynamics (WT-MetaD), and parallel-bias WT-MetaD (amounting to 3.5 µs) were integrated to unravel the chemical coordination of Zn(II) in all Zn1-6MT2 species and to explain the differences in binding affinities of Zn(II) ions to MTs. Differences are found to be the result of the degree of water participation in MT (un)folding and the hyper-reactive character of Cys21 and Cys29 residues. The thermodynamics properties of Zn(II) (un)binding to MT2 are found to differ from those of Cd(II), justifying their distinctive roles. The potential of this integrated strategy in the investigation of numerous unexplored metalloproteins is attested by the results highlighted in the present study.

Multiple Interactions of Glucose with the Extra-Membranous Loops of GLUT1 Aid Transport.

Molecular dynamics simulations amounting to ≈8 µs demonstrate that the glucose transporter GLUT1 undergoes structural fluctuations mediated by the fluidity of the lipid bilayer and the proximity to glucose. The fluctuations of GLUT1 increase as the glucose concentration is raised. These fluctuations are more pronounced when the lipid bilayer is in the fluid compared to the gel phase. Glucose interactions are confined to the extra-membranous residues when the lipid is in the gel phase but diffuses into the transmembrane regions in the fluid phase. Proximity of glucose to GLUT1 causes asynchronous expansions of key bottlenecks at the internal and external openings of the central pore. This is accomplished only by small conformational changes at the single residue level that lower the resistance to glucose movements, thereby permitting unsteered glucose and water movements along the entire length of the pore. When glucose is near salt bridges located at the external and internal openings of the central pore, the distance separating the polar amino acid residues guarding these apertures tends to increase in both fluid and gel phases. It is evident that the multiplicity of glucose interactions, obtained with high concentrations, amplifies the structural fluctuations in GLUT1. The findings that most of the salt bridges and the bottlenecks appear to be operated by glucose proximity suggest that the main triggers to activation of transport are located within the solvent accessible linker regions in the extramembranous zones.

Ion Conduction Mechanism as a Fingerprint of Potassium Channels.

K+-channels are membrane proteins that regulate the selective conduction of potassium ions across cell membranes. Although the atomic mechanisms of K+ permeation have been extensively investigated, previous work focused on characterizing the selectivity and occupancy of the binding sites, the role of water molecules in the conduction process, or the identification of the minimum energy pathways enabling permeation. Here, we exploit molecular dynamics simulations and the analytical power of Markov state models to perform a comparative study of ion conduction in three distinct channel models. Significant differences emerged in terms of permeation mechanisms and binding site occupancy by potassium ions and/or water molecules from 100 µs cumulative trajectories. We found that, at odds with the current paradigm, each system displays a characteristic permeation mechanism, and thus, there is not a unique way by which potassium ions move through K+-channels. The high functional diversity of K+-channels can be attributed in part to the differences in conduction features that have emerged from this work. This study provides crucial information and further inspiration for wet-lab chemists designing new synthetic strategies to produce versatile artificial ion channels that emulate membrane transport for their applications in diagnosis, sensors, the next generation of water treatment technologies, etc., as the ability of synthetic channels to transport molecular ions across a bilayer in a controlled way is usually governed through the choice of metal ions, their oxidation states, or their coordination geometries.

Structure based analysis of KATP channel with a DEND syndrome mutation in murine skeletal muscle.

Developmental delay, epilepsy, and neonatal diabetes (DEND) syndrome, the most severe end of neonatal diabetes mellitus, is caused by mutation in the ATP-sensitive potassium (KATP) channel. In addition to diabetes, DEND patients present muscle weakness as one of the symptoms, and although the muscle weakness is considered to originate in the brain, the pathological effects of mutated KATP channels in skeletal muscle remain elusive. Here, we describe the local effects of the KATP channel on muscle by expressing the mutation present in the KATP channels of the DEND syndrome in the murine skeletal muscle cell line C2C12 in combination with computer simulation. The present study revealed that the DEND mutation can lead to a hyperpolarized state of the muscle cell membrane, and molecular dynamics simulations based on a recently reported high-resolution structure provide an explanation as to why the mutation reduces ATP sensitivity and reveal the changes in the local interactions between ATP molecules and the channel.

The ligand-bound state of a G protein-coupled receptor stabilizes the interaction of functional cholesterol molecules.

Cholesterol is a major component of mammalian plasma membranes that not only affects the physical properties of the lipid bilayer but also is the function of many membrane proteins including G protein-coupled receptors. The oxytocin receptor (OXTR) is involved in parturition and lactation of mammals and in their emotional and social behaviors. Cholesterol acts on OXTR as an allosteric modulator inducing a high-affinity state for orthosteric ligands through a molecular mechanism that has yet to be determined. Using the ion channel-coupled receptor technology, we developed a functional assay of cholesterol modulation of G protein-coupled receptors that is independent of intracellular signaling pathways and operational in living cells. Using this assay, we discovered a stable binding of cholesterol molecules to the receptor when it adopts an orthosteric ligand-bound state. This stable interaction preserves the cholesterol-dependent activity of the receptor in cholesterol-depleted membranes. This mechanism was confirmed using time-resolved FRET experiments on WT OXTR expressed in CHO cells. Consequently, a positive cross-regulation sequentially occurs in OXTR between cholesterol and orthosteric ligands.

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.

Critical Assessment of Common Force Fields for Molecular Dynamics Simulations of Potassium Channels.

For the last two decades, the KcsA K+ channel has served as a case study to understand how potassium channels operate at the atomic scale, and molecular dynamics simulations have contributed significantly to the current knowledge of the atomic mechanisms of conduction, inactivation, and gating in this family of membrane proteins. Currently, microsecond trajectories are becoming the new standard in the field, and consequently, it is critical to assess and compare the performance of the classical force fields ordinarily used in simulations of biological systems as well as to quantitatively assess the agreement with experimental data for trajectories of this order of magnitude. To that extent, we performed classical molecular dynamics simulations with CHARMM36 and AMBER-ff14sb force fields using atomic models based on the experimental structure of the KcsA channel in the open/conductive state, at conditions of ionic concentrations and membrane potentials resembling the ones adopted in experiments. In simulations using the CHARMM force field, the experimental open/conductive structure of the channel exhibited high conformational plasticity and fast collapse toward an occluded state. In contrast, in an identical set of simulations using the AMBER force field, no major deviations from the experimental structure were recorded. Force field development is a complex process in which many approximations are typically required and adopted. The results presented here provide additional motivation to further improve the existing models and, crucially, alert practitioners about limitations.

Affinity for the Interface Underpins Potency of Antibodies Operating In Membrane Environments.

The contribution of membrane interfacial interactions to recognition of membrane-embedded antigens by antibodies is currently unclear. This report demonstrates the optimization of this type of antibodies via chemical modification of regions near the membrane but not directly involved in the recognition of the epitope. Using the HIV-1 antibody 10E8 as a model, linear and polycyclic synthetic aromatic compounds are introduced at selected sites. Molecular dynamics simulations predict the favorable interactions of these synthetic compounds with the viral lipid membrane, where the epitope of the HIV-1 glycoprotein Env is located. Chemical modification of 10E8 with aromatic acetamides facilitates the productive and specific recognition of the native antigen, partially buried in the crowded environment of the viral membrane, resulting in a dramatic increase of its capacity to block viral infection. These observations support the harnessing of interfacial affinity through site-selective chemical modification to optimize the function of antibodies that target membrane-proximal epitopes.

Differential Stability of Aurein 1.2 Pores in Model Membranes of Two Probiotic Strains.

Aurein 1.2 is an antimicrobial peptide from the skin secretion of an Australian frog. In the previous experimental work, we reported a differential action of aurein 1.2 on two probiotic strains Lactobacillus delbrueckii subsp. bulgaricus (CIDCA 331) and Lactobacillus delbrueckii subsp. lactis (CIDCA 133). The differences found were attributed to the bilayer compositions. Cell cultures and CIDCA 331-derived liposomes showed higher susceptibility than the ones derived from the CIDCA 133 strain, leading to content leakage and structural disruption. Here, we used molecular dynamics simulations to explore these systems at the atomistic level. We hypothesize that if the antimicrobial peptides organized themselves to form a pore, it will be more stable in membranes that emulate the CIDCA 331 strain than in those of the CIDCA 133 strain. To test this hypothesis, we simulated preassembled aurein 1.2 pores embedded into bilayer models that emulate the two probiotic strains. It was found that the general behavior of the systems depends on the composition of the membrane rather than the preassemble system characteristics. Overall, it was observed that aurein 1.2 pores are more stable in the CIDCA 331 model membranes. This fact coincides with the high susceptibility of this strain against antimicrobial peptide. In contrast, in the case of the CIDCA 133 model membranes, peptides migrate to the water-lipid interphase, the pore shrinks, and the transport of water through the pore is reduced. The tendency of glycolipids to make hydrogen bonds with peptides destabilizes the pore structures. This feature is observed to a lesser extent in CIDCA 331 due to the presence of anionic lipids. Glycolipid transverse diffusion (flip-flop) between monolayers occurs in the pore surface region in all the cases considered. These findings expand our understanding of the antimicrobial peptide resistance properties of probiotic strains.

Effect of anionic lipids on ion permeation through the KcsA K+-channel.

K+-channels are responsible for the efficient and selective conduction of K+ ions across the plasma membrane. The bacterial K+-channel KcsA has historically been used to characterize various aspects of K+ conduction via computational means. The energetic barriers associated with ion translocation across the KcsA selectivity filter have been computed in various studies, leading to the proposal of two alternate mechanisms of conduction, involving or neglecting the presence of water molecules in between the permeating ions. Here, the potential of mean force of K+ permeation is evaluated for KcsA in lipid bilayers containing anionic lipids, which is known to increase the open probability of the channel. In addition, the effect of the protonation/deprotonation of residue E71, which directly interacts with the selectivity filter sequence, is assessed. Both conduction mechanisms are considered throughout. The results obtained provide novel insights into the molecular functioning of K+-channels including the inactivation process.

Mechanism of Molecular Oxygen Diffusion in a Hypoxia-Sensing Prolyl Hydroxylase Using Multiscale Simulation.

The chronic response of animals to hypoxia is mediated by the αß-heterodimeric hypoxia-inducible transcription factors (α,ß-HIFs) which upregulate the expression of sets of genes that work to ameliorate the effects of limiting dioxygen. The HIF prolyl hydroxylase domain enzymes (PHDs) are Fe(II)- and 2-oxoglutarate-dependent oxygenases that act as hypoxia-sensing components of the HIF system: prolyl-hydroxylation signals for dioxgen availability-dependent HIF-α degradation via the ubiquitin proteasome system. The unusual kinetic properties of the PHDs, in particular a high Km for dioxygen and slow reaction with dioxygen, are proposed to enable their hypoxia-sensing role. An understanding of how dioxygen is delivered to, and binds at, the active site of the PHDs is important for the development of a chemical understanding of the hypoxic response. We employed a combined multiscale approach involving classical atomistic equilibrium and nonequilibrium MD simulations combined with QM/MM trajectories to investigate dioxygen diffusion to, and binding at, the active site in the PHD2.Fe(II).2OG.HIF substrate complex; PHD2 is the most important of the three human PHDs. The transport of dioxygen to the active site is described; dioxygen transport follows a single well-defined hydrophobic tunnel, formed from both enzyme and substrate elements to reach the PHD2 active site. The results provide estimates for rate constants that define a diffusion-reaction model for dioxygen:PHD2 interactions; in combination with reported biophysical analyses they provide chemical insight into the basis of the slow reaction of PHD2 with dioxygen. They imply that the reversible binding of dioxygen is central to the hypoxia-sensing capacity of the PHDs and that different PHD HIF-α substrate combinations might have different dioxygen sensitivity profiles. The extent of HIF-α substrate prolyl hydroxylation, which signals for subsequent HIF-α degradation, may thus be a manifestation of the equilibrium between dioxygen in bulk solution and dioxygen bound to the PHD2.Fe.2OG.HIF-α substrate complex.

Insights into the Mechanisms of K+ Permeation in K+ Channels from Computer Simulations.

Ion permeation, selectivity, and the behavior of the K+ channel selectivity filter have been studied intensively in the previous two decades. The agreement among multiple approaches used to study ion flux in K+ channels suggests a consensus mechanism of ion permeation across the selectivity that has been put to the test in recent years with the proposal of an alternative way by which ions can cross the selectivity filter of K+ channels via direct Coulomb repulsion between contacting cations. Past experimental work by Zhou and MacKinnon (J. Mol. Biol. 2004, 338, 839) showed that mutation of the site S4 reduces the total occupancy of the selectivity filter to less than two ions on average by lowering the occupancy of the S2-S4 configuration without changing the S1-S3 configuration much, and this reduction of occupancy means that ion configurations different from the ones involved in the canonical mechanism are likely to be involved. At that time, calculations using complicated kinetic networks to relate occupancy to conduction did not provide deeper insight into the conduction mechanism. Here, to help solve this enigma, umbrella sampling simulations have been performed to evaluate the potential of mean force of two KcsA mutant channels where the S4 site is substituted. Our new results provide insights into the significance of threonine in this position, revealing the effect of substitution on the alternate mechanisms of conduction proposed, involving either water or vacant sites.

The Bilayer Collective Properties Govern the Interaction of an HIV-1 Antibody with the Viral Membrane.

Efficient engagement with the envelope glycoprotein membrane-proximal external region (MPER) results in robust blocking of viral infection by a class of broadly neutralizing antibodies (bnAbs) against human immunodeficiency virus (HIV). Developing an accommodation surface that engages with the viral lipid envelope appears to correlate with the neutralizing potency displayed by these bnAbs. The nature of the interactions established between the antibody and the lipid is nonetheless a matter of debate, with some authors arguing that anti-MPER specificity arises only under pathological conditions in autoantibodies endowed with stereospecific binding sites for phospholipids. However, bnAb-lipid interactions are often studied in systems that do not fully preserve the biophysical properties of lipid bilayers, and therefore, questions on binding specificity and the effect of collective membrane properties on the interaction are still open. Here, to evaluate the specificity of lipid interactions of an anti-MPER bnAb (4E10) in an intact membrane context, we determine quantitatively its association with lipid bilayers by means of scanning fluorescence correlation spectroscopy and all-atom molecular dynamic simulations. Our data support that 4E10 establishes electrostatic and hydrophobic interactions with the viral membrane surface and that the collective physical properties of the lipid bilayer influence 4E10 dynamics therein. We conclude that establishment of peripheral, nonspecific electrostatic interactions with the viral membrane through accommodation surfaces may assist high-affinity binding of HIV-1 MPER epitope at membrane interfaces. These findings highlight the importance of considering antibody-lipid interactions in the design of antibody-based anti-HIV strategies.

Modulation of the potassium channel KcsA by anionic phospholipids: Role of arginines at the non-annular lipid binding sites.

The role of arginines R64 and R89 at non-annular lipid binding sites of KcsA, on the modulation of channel activity by anionic lipids has been investigated. In wild-type (WT) KcsA reconstituted into asolectin lipid membranes, addition of phosphatidic acid (PA) drastically reduces inactivation in macroscopic current recordings. Consistent to this, PA increases current amplitude, mean open time and open probability at the single channel level. Moreover, kinetic analysis reveals that addition of PA causes longer open channel lifetimes and decreased closing rate constants. Effects akin to those of PA on WT-KcsA are observed when R64 and/or R89 are mutated to alanine, regardless of the added anionic lipids. We interpret these results as a consequence of interactions between the arginines and the anionic PA bound to the non-annular sites. NMR data shows indeed that at least R64 is involved in binding PA. Moreover, molecular dynamics (MD) simulations predict that R64, R89 and surrounding residues such as T61, mediate persistent binding of PA to the non-annular sites. Channel inactivation depends on interactions within the inactivation triad (E71-D80-W67) behind the selectivity filter. Therefore, it is expected that such interactions are affected when PA binds the arginines at the non-annular sites. In support of this, MD simulations reveal that PA binding prevents interaction between R89 and D80, which seems critical to the effectiveness of the inactivation triad. This mechanism depends on the stability of the bound lipid, favoring anionic headgroups such as that of PA, which thrive on the positive charge of the arginines.