Potential blockade of the human voltage-dependent anion channel by MoS2 nanoflakes.

Despite significant interest in molybdenum disulfide (MoS2) nanomaterials, particularly in biomedicine, their biological effects have been understudied. Here, we explored the effect of MoS2 nanoflakes on the ubiquitous mitochondrial porin voltage-dependent anion channel (VDAC1), using a combined computational and functional approach. All-atomic molecular dynamics simulations suggest that MoS2 nanoflakes make specific contact interactions with human VDAC1. We show that the initial contacts between hVDAC1 and the nanoflake are hydrophobic but are subsequently enhanced by a complex interplay of van der Waals (vdW), hydrophobic and electrostatic interactions in the equilibrium state. Moreover, the MoS2 nanoflake can insert into the lumen of the hVDAC1 pore. Free-energy calculations computed by the potential of mean force (PMF) verify that the blocked configuration of the MoS2-hVDAC1 complex is more energetically favorable than the non-blocked binding mode. Consistent with these predictions, we showed that MoS2 depolarizes the mitochondrial membrane potential (Ψm) and causes a decrease in the viability of mammalian tissue culture cells. These findings might shed new light on the potential biological effect of MoS2 nanomaterials.

Phosphoinositide control of membrane protein function: a frontier led by studies on ion channels.

Anionic phospholipids are critical constituents of the inner leaflet of the plasma membrane, ensuring appropriate membrane topology of transmembrane proteins. Additionally, in eukaryotes, the negatively charged phosphoinositides serve as key signals not only through their hydrolysis products but also through direct control of transmembrane protein function. Direct phosphoinositide control of the activity of ion channels and transporters has been the most convincing case of the critical importance of phospholipid-protein interactions in the functional control of membrane proteins. Furthermore, second messengers, such as [Ca(2+)]i, or posttranslational modifications, such as phosphorylation, can directly or allosterically fine-tune phospholipid-protein interactions and modulate activity. Recent advances in structure determination of membrane proteins have allowed investigators to obtain complexes of ion channels with phosphoinositides and to use computational and experimental approaches to probe the dynamic mechanisms by which lipid-protein interactions control active and inactive protein states.

A Critical Gating Switch at a Modulatory Site in Neuronal Kir3 Channels.

Inwardly rectifying potassium channels enforce tight control of resting membrane potential in excitable cells. The Kir3.2 channel, a member of the Kir3 subfamily of G-protein-activated potassium channels (GIRKs), plays several roles in the nervous system, including key responsibility in the GABAB pathway of inhibition, in pain perception pathways via opioid receptors, and is also involved in alcoholism. PKC phosphorylation acts on the channel to reduce activity, yet the mechanism is incompletely understood. Using the heterologous Xenopus oocyte system combined with molecular dynamics simulations, we show that PKC modulation of channel activity is dependent on Ser-196 in Kir3.2 such that, when this site is phosphorylated, the channel is less sensitive to PKC inhibition. This reduced inhibition is dependent on an interaction between phospho-Ser (SEP)-196 and Arg-201, reducing Arg-201 interaction with the sodium-binding site Asp-228. Neutralization of either SEP-196 or Arg-201 leads to a channel with reduced activity and increased sensitivity to PKC inhibition. This study clarifies the role of Ser-196 as an allosteric modulator of PKC inhibition and suggests that the SEP-196/Arg-201 interaction is critical for maintaining maximal channel activity. SIGNIFICANCE STATEMENT: The inwardly rectifying potassium 3.2 (Kir3.2) channel is found principally in neurons that regulate diverse brain functions, including pain perception, alcoholism, and substance addiction. Activation or inhibition of this channel leads to changes in neuronal firing and chemical message transmission. The Kir3.2 channel is subject to regulation by intracellular signals including sodium, G-proteins, ethanol, the phospholipid phosphatidylinositol bis-phosphate, and phosphorylation by protein kinases. Here, we take advantage of the recently published structure of Kir3.2 to provide an in-depth molecular view of how phosphorylation of a specific residue previously thought to be the target of PKC promotes channel gating and acts as an allosteric modulator of PKC-mediated inhibition.

Selective phosphorylation modulates the PIP2 sensitivity of the CaM-SK channel complex.

Phosphatidylinositol bisphosphate (PIP2) regulates the activities of many membrane proteins, including ion channels, through direct interactions. However, the affinity of PIP2 is so high for some channel proteins that its physiological role as a modulator has been questioned. Here we show that PIP2 is a key cofactor for activation of small conductance Ca2+-activated potassium channels (SKs) by Ca(2+)-bound calmodulin (CaM). Removal of the endogenous PIP2 inhibits SKs. The PIP2-binding site resides at the interface of CaM and the SK C terminus. We further demonstrate that the affinity of PIP2 for its target proteins can be regulated by cellular signaling. Phosphorylation of CaM T79, located adjacent to the PIP2-binding site, by casein kinase 2 reduces the affinity of PIP2 for the CaM-SK channel complex by altering the dynamic interactions among amino acid residues surrounding the PIP2-binding site. This effect of CaM phosphorylation promotes greater channel inhibition by G protein-mediated hydrolysis of PIP2.

Structural determinants of phosphatidylinositol 4,5-bisphosphate (PIP2) regulation of BK channel activity through the RCK1 Ca2+ coordination site.

Big or high conductance potassium (BK) channels are activated by voltage and intracellular calcium (Ca(2+)). Phosphatidylinositol 4,5-bisphosphate (PIP2), a ubiquitous modulator of ion channel activity, has been reported to enhance Ca(2+)-driven gating of BK channels, but a molecular understanding of this interplay or even of the PIP2 regulation of this channel's activity remains elusive. Here, we identify structural determinants in the KDRDD loop (which follows the αA helix in the RCK1 domain) to be responsible for the coupling between Ca(2+) and PIP2 in regulating BK channel activity. In the absence of Ca(2+), RCK1 structural elements limit channel activation through a decrease in the channel's PIP2 apparent affinity. This inhibitory influence of BK channel activation can be relieved by mutation of residues that (a) connect either the RCK1 Ca(2+) coordination site (Asp(367) or its flanking basic residues in the KDRDD loop) to the PIP2-interacting residues (Lys(392) and Arg(393)) found in the αB helix or (b) are involved in hydrophobic interactions between the αA and αB helix of the RCK1 domain. In the presence of Ca(2+), the RCK1-inhibitory influence of channel-PIP2 interactions and channel activity is relieved by Ca(2+) engaging Asp(367). Our results demonstrate that, along with Ca(2+) and voltage, PIP2 is a third factor critical to the integral control of BK channel activity.

PIP2 controls voltage-sensor movement and pore opening of Kv channels through the S4-S5 linker.

Voltage-gated K(+) (Kv) channels couple the movement of a voltage sensor to the channel gate(s) via a helical intracellular region, the S4-S5 linker. A number of studies link voltage sensitivity to interactions of S4 charges with membrane phospholipids in the outer leaflet of the bilayer. Although the phospholipid phosphatidylinositol-4,5-bisphosphate (PIP(2)) in the inner membrane leaflet has emerged as a universal activator of ion channels, no such role has been established for mammalian Kv channels. Here we show that PIP(2) depletion induced two kinetically distinct effects on Kv channels: an increase in voltage sensitivity and a concomitant decrease in current amplitude. These effects are reversible, exhibiting distinct molecular determinants and sensitivities to PIP(2). Gating current measurements revealed that PIP(2) constrains the movement of the sensor through interactions with the S4-S5 linker. Thus, PIP(2) controls both the movement of the voltage sensor and the stability of the open pore through interactions with the linker that connects them.

Building KCNQ1/KCNE1 channel models and probing their interactions by molecular-dynamics simulations.

The slow delayed rectifier (I(KS)) channel is composed of KCNQ1 (pore-forming) and KCNE1 (auxiliary) subunits, and functions as a repolarization reserve in the human heart. Design of I(KS)-targeting anti-arrhythmic drugs requires detailed three-dimensional structures of the KCNQ1/KCNE1 complex, a task made possible by Kv channel crystal structures (templates for KCNQ1 homology-modeling) and KCNE1 NMR structures. Our goal was to build KCNQ1/KCNE1 models and extract mechanistic information about their interactions by molecular-dynamics simulations in an explicit lipid/solvent environment. We validated our models by confirming two sets of model-generated predictions that were independent from the spatial restraints used in model-building. Detailed analysis of the molecular-dynamics trajectories revealed previously unrecognized KCNQ1/KCNE1 interactions, whose relevance in I(KS) channel function was confirmed by voltage-clamp experiments. Our models and analyses suggest three mechanisms by which KCNE1 slows KCNQ1 activation: by promoting S6 bending at the Pro hinge that closes the activation gate; by promoting a downward movement of gating charge on S4; and by establishing a network of electrostatic interactions with KCNQ1 on the extracellular surface that stabilizes the channel in a pre-open activated state. Our data also suggest how KCNE1 may affect the KCNQ1 pore conductance.

The cytosolic GH loop regulates the phosphatidylinositol 4,5-bisphosphate-induced gating kinetics of Kir2 channels.

Inwardly rectifying K(+) (Kir) channels set the resting membrane potential and regulate cellular excitability. The activity of Kir channels depends critically on the phospholipid PIP(2). The molecular mechanism by which PIP(2) regulates Kir channel gating is poorly understood. Here, we utilized a combination of computational and electrophysiological approaches to discern structural elements involved in regulating the PIP(2)-induced gating kinetics of Kir2 channels. We identify a novel role for the cytosolic GH loop. Mutations that directly or indirectly affect GH loop flexibility (e.g. V223L, E272G, D292G) increase both the on- and especially the off-gating kinetics. These effects are consistent with a model in which competing interactions between the CD and GH loops for the N terminus regulate the gating of the intracellular G loop gate.

Identification of driving genes of familial adenomatous polyposis by differential gene expression analysis and weighted gene co-expression network analysis.

BACKGROUND: Despite the advancement of new screening strategies and the advances in pharmacological therapies, the cancerization rates of familial adenomatous polyposis (FAP) are stable and even increased in the last years. Therefore, it necessitates additional research to characterize and understand the underlying mechanisms of FAP. OBJECTIVE: To determine the genes that drive the pathogenesis of familial adenomatous polyposis (FAP). METHODS: We performed on a cohort (GSE111156) gene profile, which consist of four group of gene expressions (the gene expressions of cancer, adenoma and normal tissue of duodenal cancer from patients with FAP were defined as Case N, Case A and Case C respectively, while that of adenoma tissue from patients with FAP who did not have duodenal cancer was Ctrl A). Tracking Tumor Immunophenotype (TIP) website was applied to reveal immune infiltration profile and signature genes of FAP. We merged the genes of key module (pink and midnight module) with signature genes to obtained the biomarkers related with FAP pathogenesis. The expression of these five biomarkers in FAP intratumoral region (IT) and tumor rim (TR) was detected with Quantitative Real-Time Polymerase Chain Reaction (qRT-PCR). RESULTS: In total, 220, 23 and 63 DEGs were determined in Cases C, A and N, in comparison to Ctrl A. In total, 196 and 10 DEGs were determined in Cases C and A, separately, as compared to Case N. A total of four biomarkers including CCL5, CD3G, CD2 and TLR3 were finally identified associated with pink module, while only one biomarker (KLF2) associated with midnight module was identified. All biomarkers were evidently raised in FAP IT tissues utilizing qRT-PCR. CONCLUSION: We identified five potential biomarkers for pathogenesis of FAP to understand the fundamental mechanisms of FAP progression and revealed some probable targets for the diagnosis or treatment of FAP.

Somatic mutation profiling, tumor-infiltrating leukocytes, tertiary lymphoid structures and PD-L1 protein expression in HER2-amplified colorectal cancer.

The status of human epidermal growth factor receptor 2 (HER2) for the prognosis in colorectal cancer (CRC) is controversial, and the characteristics of the somatic mutation spectrum, tumor-infiltrating leukocytes, tertiary lymphoid structures and PD-L1 protein are unknown in HER2-amplified colorectal cancer (HACC). In order to explore these characteristics along with their correlation with clinicopathological factors and prognosis in HACC. Samples of 812 CRC patients was collected. After immunohistochemistry (IHC), 59 of 812 were found to be HER2-positive, then 26 of 59 samples were further determined to be HER2 amplification by fluorescence in situ hybridization (FISH). Somatic mutation profiling of HACC was analysed using whole exome sequencing (WES). Multiplex fluorescence immunohistochemistry (mIHC) was used for tumor-infiltrating leukocytes and tertiary lymphoid structures (TLSs), while PD-L1 protein was detected by IHC. Our results indicate that the detection rates of HER2 positivity by IHC and FISH were 7.3% and 3.2% respectively, and HER2 amplification is correlated with distant tumour metastasis. The somatic mutation profiling revealed no differences between HACC and HER2-negative CRC. However, TP 53 strongly correlated with poor prognosis in HACC. Furthermore, tumor-infiltrating T cells and TLSs in the tumor immune microenvironment, as well as PD-L1 expression, were higher in HACC than in HER2-negative controls. However, none of them were associated with the prognosis of HACC. In all, HER2 amplification is correlated with distant metastasis and TP53 gene mutation may be a potential protective mechanism of HACC.

The molecular mechanism by which PIP(2) opens the intracellular G-loop gate of a Kir3.1 channel.

Inwardly rectifying potassium (Kir) channels are characterized by a long pore comprised of continuous transmembrane and cytosolic portions. A high-resolution structure of a Kir3.1 chimera revealed the presence of the cytosolic (G-loop) gate captured in the closed or open conformations. Here, we conducted molecular-dynamics simulations of these two channel states in the presence and absence of phosphatidylinositol bisphosphate (PIP(2)), a phospholipid that is known to gate Kir channels. Simulations of the closed state with PIP(2) revealed an intermediate state between the closed and open conformations involving direct transient interactions with PIP(2), as well as a network of transitional inter- and intrasubunit interactions. Key elements in the G-loop gating transition involved a PIP(2)-driven movement of the N-terminus and C-linker that removed constraining intermolecular interactions and led to CD-loop stabilization of the G-loop gate in the open state. To our knowledge, this is the first dynamic molecular view of PIP(2)-induced channel gating that is consistent with existing experimental data.

Molecular mechanism of species-dependent sweet taste toward artificial sweeteners.

The heterodimer of Tas1R2 and Tas1R3 is a broadly acting sweet taste receptor, which mediates mammalian sweet taste toward natural and artificial sweeteners and sweet-tasting proteins. Perception of sweet taste is a species-selective physiological process. For instance, artificial sweeteners aspartame and neotame taste sweet to humans, apes, and Old World monkeys but not to New World monkeys and rodents. Although specific regions determining the activation of the receptors by these sweeteners have been identified, the molecular mechanism of species-dependent sweet taste remains elusive. Using human/squirrel monkey chimeras, mutagenesis, and molecular modeling, we reveal that the different responses of mammalian species toward the artificial sweeteners aspartame and neotame are determined by the steric effect of a combination of a few residues in the ligand binding pocket. Residues S40 and D142 in the human Tas1R2, which correspond to residues T40 and E142 in the squirrel monkey Tas1R2, were found to be the critical residues for the species-dependent difference in sweet taste. In addition, human Tas1R2 residue I67, which corresponds to S67 in squirrel monkey receptor, modulates the higher affinity of neotame than of aspartame. Our studies not only shed light on the molecular mechanism of species-dependent sweet taste toward artificial sweeteners, but also provide guidance for designing novel effective artificial sweet compounds.

Functional characterization of the heterodimeric sweet taste receptor T1R2 and T1R3 from a New World monkey species (squirrel monkey) and its response to sweet-tasting proteins.

The family C G protein-coupled receptor (GPCR) T1R2 and T1R3 heterodimer functions as a broadly acting sweet taste receptor. Perception of sweet taste is a species-dependent physiological process. It has been widely reported that New World monkeys and rodents are not able to perceive some of the artificial sweeteners and sweet-tasting proteins that can be perceived by humans, apes, and Old World monkeys. Until now, only the sweet receptors of humans, mice and rats have been functionally characterized. Here we report characterization of the sweet taste receptor (T1R2/T1R3) from a species of New World primate, squirrel monkey. Our results show that the heterodimeric receptor of squirrel monkey does not respond to artificial sweeteners aspartame, neotame, cyclamate, saccharin and sweet-tasting protein monellin, but surprisingly, it does respond to thaumatin at high concentrations (>18 µM). This is the first report demonstrating that species of New World monkey can perceive some specific sweet-tasting proteins. Furthermore, the sweet receptor of squirrel monkey responses to the such sweeteners cannot be inhibited by the sweet inhibitor lactisole. We compared the response differences of the squirrel monkey and human receptors and found that the residues in T1R2 determine species-dependent sweet taste toward saccharin, while the residues in either T1R2 or T1R3 are responsible for the sweet taste difference between humans and squirrel monkeys toward monellin. Molecular models indicated that electrostatic properties of the receptors probably mediate the species-dependent response to sweet-tasting proteins.

Predicting protein interactions by Brownian dynamics simulations.

We present a newly adapted Brownian-Dynamics (BD)-based protein docking method for predicting native protein complexes. The approach includes global BD conformational sampling, compact complex selection, and local energy minimization. In order to reduce the computational costs for energy evaluations, a shell-based grid force field was developed to represent the receptor protein and solvation effects. The performance of this BD protein docking approach has been evaluated on a test set of 24 crystal protein complexes. Reproduction of experimental structures in the test set indicates the adequate conformational sampling and accurate scoring of this BD protein docking approach. Furthermore, we have developed an approach to account for the flexibility of proteins, which has been successfully applied to reproduce the experimental complex structure from the structure of two unbounded proteins. These results indicate that this adapted BD protein docking approach can be useful for the prediction of protein-protein interactions.

The prognostic and clinicopathological roles of microsatellite instability, PD-L1 expression and tumor-infiltrating leukocytes in familial adenomatous polyposis.

BACKGROUND: Microsatellite instability, programmed death-ligand 1 and tumor-infiltrating leukocytes are prognostic biomarkers in colorectal cancer but unknown toward familial adenomatous polyposis. AIM: To investigate the prognostic and clinicopathological roles of microsatellite instability, programmed death-ligand 1 and tumor-infiltrating leukocytes in familial adenomatous polyposis. METHODS: Clinical data and paraffin embedded tissues from 45 familial adenomatous polyposis patients were collected. Microsatellite instability was detected by immunohistochemistry and polymerase chain reaction. Programmed death-ligand 1 was detected by immunohistochemistry. Tumor-infiltrating leukocytes comprising CD8+ T cells, M1 and M2 tumor associated macrophages, CD56bright and CD56dim natural killer cells were analyzed using multiple fluorescence immunohistochemistry. RESULTS: Microsatellite instability high was noted in 6 samples but not associated with overall survival or progression-free survival. Programmed death-ligand 1 is negative on tumor cells but positive on tumor-infiltrating leukocytes, and positive programmed death-ligand 1 expression on tumor-infiltrating leucocytes is associated with overall survival. Low CD56bright natural killer cell infiltration was associated with longer progression-free survival and was an independent prognostic factor in FAP. CONCLUSION: For familial adenomatous polyposis, microsatellite instability high can be found but has no correlation with prognosis; programmed death-ligand 1 on tumor-infiltrating leukocytes is related with overall survival; CD56bright natural killer cell is an independent prognostic factor associating with longer progression-free survival.

Molecular Dynamics Simulation Study on Interactions of Cycloviolacin with Different Phospholipids.

Cyclotides are disulfide-rich cyclic peptides isolated from plants, which are extremely stable against thermal and proteolytic degradation, with a variety of biological activities including antibacterial, hemolytic, anti-HIV, and anti-tumor. Most of these bioactivities are related to their preference for binding to certain types of phospholipids and subsequently disrupt lipid membranes. In the present study, we use a cyclotide, cycloviolacin O2 (cyO2), as a model system to investigate its interactions with three lipid bilayers 1-palmitoyl-2-oleoylphosphatidylethanolamine (POPE), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoglycerol (POPG)-doped POPE, and 1-palmitoyl-2-oleoylphosphatidylcholine (POPC), to help understand its potential mechanism of action toward the membranes at the molecular level using molecular dynamics simulations. In our simulations, cyO2 repeatedly forms stable binding complexes with the POPE-containing bilayers, while within the same simulation time scale, it "jumps" back and forth on the surface of the POPC bilayer without a strong binding. Detailed analyses reveal that the electrostatic attraction is the main driving force for the initial bindings between cyO2 and the lipids, but with strikingly different strengths in different bilayers. For the POPE-containing bilayers, the charged residues of cyO2 attract both POPE amino and phosphate head groups favorably; meanwhile, its hydrophobic residues are deeply inserted into the lipid hydrophobic tails (core) of the membrane, thus forming stable binding complexes. In contrast, POPC lipids with three methyl groups on the amino head group create a steric hindrance when interacting with cyO2, thus resulting in a relatively difficult binding of cyO2 on POPC compared to POPE. Our current findings provide additional insights for a better understanding of how cyO2 binds to the POPE-containing membrane, which should shed light on the future cyclotide-based antibacterial agent design.

Multifaceted Regulation of Potassium-Ion Channels by Graphene Quantum Dots.

Graphene quantum dots (GQDs) are emerging as a versatile nanomaterial with numerous proposed biomedical applications. Despite the explosion in potential applications, the molecular interactions between GQDs and complex biomolecular systems, including potassium-ion (K+) channels, remain largely unknown. Here, we use molecular dynamics (MD) simulations and electrophysiology to study the interactions between GQDs and three representative K+ channels, which participate in a variety of physiological processes and are closely related to many disease states. Using MD simulations, we observed that GQDs adopt distinct contact poses with each of the three structurally distinct K+ channels. Our electrophysiological characterization of the effects of GQDs on channel currents revealed that GQDs interact with the extracellular voltage-sensing domain (VSD) of a Kv1.2 channel, augmenting current by left-shifting the voltage dependence of channel activation. In contrast, GQDs form a "lid" cluster over the extracellular mouth of inward rectifier Kir3.2, blocking the channel pore and decreasing the current in a concentration-dependent manner. Meanwhile, GQDs accumulate on the extracellular "cap domain" of K2P2 channels and have no apparent impact on the K+ flux through the channel. These results reveal a surprising multifaceted regulation of K+ channels by GQDs, which might help de novo design of nanomaterial-based channel probe openers/inhibitors that can be used to further discern channel function.

A comparative analysis of binding sites between mouse CYP2C38 and CYP2C39 based on homology modeling, molecular dynamics simulation and docking studies.

Mouse CYP2C38 and CYP2C39 are two closely related enzymes with 91.8% sequence identity. But they exhibit different substrate binding features. In this study, three-dimensional models of CYP2C38 and CYP2C39 were constructed using X-ray crystal structure of human CYP2C8 as the template based on homology modeling methods and molecular dynamics simulations. Tolbutamide as the common substrate of CYP2C38 and CYP2C39 was docked into them and positioned in their active sites with different orientation. All-trans retinoic acid (atRA) is a specific substrate for CYP2C39 and not catalyzed by CYP2C38. By comparison of active site architectures between CYP2C38 and CYP2C39, the possible reasons affecting their substrate binding were proposed. In addition, Arg241, Glu300, Leu366 and Leu476 are identified as critical residue for substrates binding.

Binding patterns and dynamics of double-stranded DNA on the phosphorene surface.

Phosphorene, a monolayer of black phosphorus, has emerged as one of the most promising two-dimensional (2D) nanomaterials for various applications in the post-graphene-discovery period due to its highly anisotropic structure and novel properties. In order to apply phosphorene in biomedical fields, it is crucial to understand how it interacts with biomolecules. Herein, we use both molecular dynamics (MD) simulations and experimental techniques to investigate the interactions of phosphorene with a dsDNA segment. Our results reveal that dsDNA can form a stable binding on the phosphorene surface through the terminal base pairs and adopt an upright orientation regardless of its initial configurations. Moreover, the binding strength of dsDNA with phosphorene is found to be mild and does not cause significant distortion in the internal structure of dsDNA. This phenomenon is attributed to the weaker dispersion interaction between dsDNA and phosphorene. Further analysis of the free energy profile calculated by the umbrella sampling technique suggests that the puckered surface morphology significantly reduces the adsorption free energy of DNA bases to phosphorene. Compared to graphene, phosphorene is found to show a milder attraction to DNA, which is confirmed by our electrophoresis experiments. We believe that these findings provide valuable insight into the molecular interactions between phosphorene and dsDNA which may prompt further investigation of phosphorene for future biomedical applications.

Molecular mechanism of phosphoinositides' specificity for the inwardly rectifying potassium channel Kir2.2.

Phosphoinositides are essential signaling lipids that play a critical role in regulating ion channels, and their dysregulation often results in fatal diseases including cardiac arrhythmia and paralysis. Despite decades of intensive research, the underlying molecular mechanism of lipid agonism and specificity remains largely unknown. Here, we present a systematic study of the binding mechanism and specificity of a native agonist, phosphatidylinositol 4,5-bisphosphate (PI(4,5)P2) and two of its variants, PI(3,4)P2 and PI(3,4,5)P3, on inwardly rectifying potassium channel Kir2.2, using molecular dynamics simulations and free energy perturbations (FEPs). Our results demonstrate that the major driving force for the PI(4,5)P2 specificity on Kir2.2 comes from the highly organized salt-bridge network formed between the charged inositol head and phosphodiester linker of PI(4,5)P2. The unsaturated arachidonic chain is also shown to contribute to the stable binding through hydrophobic interactions with nearby Kir2.2 hydrophobic residues. Consistent with previous experimental findings, our FEP results confirmed that non-native ligands, PI(3,4)P2 and PI(3,4,5)P3, show significant loss in binding affinity as a result of the substantial shift from the native binding mode and unfavorable local solvation environment. However, surprisingly, the underlying molecular pictures for the unfavorable binding of both ligands are quite distinctive: for PI(3,4)P2, it is due to a direct destabilization in the bound state, whereas for PI(3,4,5)P3, it is due to a relative stabilization in its free state. Our findings not only provide a theoretical basis for the ligand specificity, but also generate new insights into the allosteric modulation of ligand-gated ion channels.