Molecular mechanism of CD44 homodimerization modulated by palmitoylation and membrane environments.

The homodimerization of CD44 plays a key role in an intercellular-to-extracellular signal transduction and tumor progression. Acylated modification and specific membrane environments have been reported to mediate translocation and oligomerization of CD44; however, the underlying molecular mechanism remains elusive. In this study, extensive molecular dynamics simulations are performed to characterize the dimerization of palmitoylated CD44 variants in different bilayer environments. CD44 forms homodimer depending on the cysteines on the juxta-membrane domains, and the dimerization efficiency and packing configurations are defected by their palmitoylated modifications. In the phase-segregated (raft included) membrane, homodimerization of the palmitoylated CD44 is hardly observed, whereas PIP2 addition compensates to realize dimerization. However, the structure of CD44 homodimer formed in the phase-segregated bilayer turns susceptive and PIP2 addition allows for an extensive conformation of the cytoplasmic domain, a proposal prerequisite to access the cytoskeleton linker proteins. The results unravel a delicate competitive relationship between PIP2, lipid raft, and palmitoylation in mediating protein homodimerization, which helps to clarify the dynamic dimer conformations and involved cellular signaling of the CD44 likewise proteins.

The Molecular Mechanism of Ginsenoside Analogs Activating TMEM16A.

The calcium-activated chloride channel TMEM16A is involved in many physiological processes, and insufficient function of TMEM16A may lead to the occurrence of various diseases. Therefore, TMEM16A activators are supposed to be potentially useful for treatment of TMEM16A downregulation-inducing diseases. However, the TMEM16A activators are relatively rare, and the underlying activation mechanism of them is unclear. In the previous work, we have proved that ginsenoside Rb1 is a TMEM16A activator. In this work, we explored the activation mechanism of ginsenoside analogs on TMEM16A through analyzing the interactions between six ginsenoside analogs and TMEM16A. We identified GRg2 and GRf can directly activate TMEM16A by screening five novel ginsenosids analogs (GRb2, GRf, GRg2, GRh2, and NGR1). Isolated guinea pig ileum assay showed both GRg2 and GRf increased the amplitude and frequency of ileum contractions. We explored the molecular mechanisms of ginsenosides activating TMEM16A by combining molecular simulation with electrophysiological experiments. We proposed a TMEM16A activation process model based on the results, in which A697 on TM7 and L746 on TM8 bind to the isobutenyl of ginsenosides through hydrophobic interaction to fix the spatial location of ginsenosides. N650 on TM6 and E705 on TM7 bind to ginsenosides through electrostatic interaction, which causes the inner half of α-helix 6 to form physical contact with ginsenosides and leads to the pore opening. It should be emphasized that TMEM16A can be activated by ginsenosides only when both the above two conditions are satisfied. This is the first, to our knowledge, report of TMEM16A opening process activated by small-molecule activators. The mechanism of ginsenosides activating TMEM16A will provide important clues for TMEM16A gating mechanism and for new TMEM16A activators screening.

Molecular Mechanisms and Structural Basis of Retigabine Analogues in Regulating KCNQ2 Channel.

KCNQ2 channel is one of the important members of potassium voltage-gated channel. KCNQ2 is closely related to neuronal excitatory diseases including epilepsy and neuropathic pain, and also acts as a drug target of the anti-epileptic drug, retigabine (RTG). In the past few decades, RTG has shown strong efficacy in the treatment of refractory epilepsy but has been withdrawn from clinical use due to its multiple adverse effects in clinical phase III trials. To overcome the drawbacks of RTG, several RTG analogues have been developed with different activation potency to KCNQ2. However, the detailed molecular mechanism by which these RTG analogues regulate KCNQ2 channel remains obscure. In this study, we used molecular simulations to analyse the interaction mode between the RTG analogues and KCNQ2, and to determine their molecular mechanism of action. Our data show that the van der Waals interactions, hydrophobic interactions, hydrogen bond, halogen bond, and π-π stacking work together to maintain the binding stability of the drugs in the binding pocket. On an atomic scale, the amide group in the carbamate and the amino group in the 2-aminophenyl moiety of RTG and RL648_81 are identified as key interaction sites. Our finding provides insight into the molecular mechanism by which KCNQ2 channels are regulated by RTG analogues. It also provides direct theoretical support for optimizing design of the KCNQ2 channel openers in the future, which will help treat refractory epilepsy caused by nerve excitability.

Molecular mechanism of CaCCinh-A01 inhibiting TMEM16A channel.

TMEM16A is a calcium-activated chloride channel that is associate with several diseases, including pulmonary diseases, hypertension, diarrhea and cancer. The CaCCinh-A01 (A01) is widely recognized as an efficient blocker of TMEM16A and has been used as a tool drug to inhibit TMEM16A currents in the laboratory. A01 also has excellent pharmacokinetic properties and can be developed as a drug to target TMEM16A. However, the molecular mechanism how A01 inhibits TMEM16A is still elusive, which slows down its drug development process. Here, calculations identified that the binding pocket of A01 was located above the pore, and it was also discovered that the binding of A01 to TMEM16A not only blocked the pore but also led to its collapse. The interaction model analysis predicted that R515/K603/E623 were crucial residues for the binding between TMEM16A and A01, and the site-directed mutagenesis studies confirmed the above results. The binding mode and quantum chemical calculations showed that the carboxyl and the amide oxygen atom of A01 were the key interaction sites between TMEM16A and A01. Therefore, our study proposed the inhibitory mechanism of TMEM16A current by A01 and revealed how A01 inhibits TMEM16A at the molecular level. These findings will shed light on both the development of A01 as a potential drug for TMEM16A dysfunction-related disorders and drug screening targeting the pocket.

Recent advances in TMEM16A: Structure, function, and disease.

TMEM16A (also known as anoctamin 1, ANO1) is the molecular basis of the calcium-activated chloride channels, with ten transmembrane segments. Recently, atomic structures of the transmembrane domains of mouse TMEM16A (mTMEM16A) were determined by single-particle electron cryomicroscopy. This gives us a solid ground to discuss the electrophysiological properties and functions of TMEM16A. TMEM16A is reported to be dually regulated by Ca2+ and voltage. In addition, the dysfunction of TMEM16A has been found to be involved in many diseases including cystic fibrosis, various cancers, hypertension, and gastrointestinal motility disorders. TMEM16A is overexpressed in many cancers, including gastrointestinal stromal tumors, gastric cancer, head and neck squamous cell carcinoma (HNSCC), colon cancer, pancreatic ductal adenocarcinoma, and esophageal cancer. Furthermore, overexpression of TMEM16A is related to the occurrence, proliferation, and migration of tumor cells. To date, several studies have shown that many natural compounds and synthetic compounds have regulatory effects on TMEM16A. These small molecule compounds might be novel drugs for the treatment of diseases caused by TMEM16A dysfunction in the future. In addition, recent studies have shown that TMEM16A plays different roles in different diseases through different signal transduction pathways. This review discusses the topology, electrophysiological properties, modulators and functions of TMEM16A in mediates nociception, gastrointestinal dysfunction, hypertension, and cancer and focuses on multiple regulatory mechanisms regarding TMEM16A.

Matrine is a novel inhibitor of the TMEM16A chloride channel with antilung adenocarcinoma effects.

Calcium-activated chloride channels (CaCCs) are ion channels with key roles in physiological processes. They are abnormally expressed in various cancers, including esophageal squamous cell cancer, head and neck squamous cell carcinoma, colorectal cancer, and gastrointestinal stromal tumors. The CaCC component TMEM16A/ANO1 was recently shown to be overexpressed in lung adenocarcinoma cells and may serve as a tumorigenic protein. In this study, we determined that matrine is a potent TMEM16A inhibitor that exerts anti-lung adenocarcinoma effects. Patch clamp experiments showed that matrine inhibited TMEM16A current in a concentration-dependent manner with an IC 50 of 27.94 ± 4.78 µM. Molecular simulation and site-directed mutation experiments demonstrated that the matrine-sensitive sites of the TMEM16A channel involve the amino acids Y355, F411, and F415. Results of cell viability and wound healing assays showed that matrine significantly inhibited the proliferation and migration of LA795 cells, which exhibit high TMEM16A expression. In contrast, matrine has only weak inhibitory effect on CCD-19Lu and HeLa cells lacking TMEM16A expression. Matrine-induced effects on the proliferation and migration of LA795 cells were abrogated upon shRNA-mediated TMEM16A knockdown in LA795 cells. Finally, in vivo experiments demonstrated that matrine dramatically inhibited the growth of lung adenocarcinoma xenograft tumors in mice but did not affect mouse body weight. Collectively, these data indicate that matrine is an effective and safe TMEM16A inhibitor and that TMEM16A is the target of matrine anti-lung adenocarcinoma activity. These findings provide new insight for the development of novel treatments for lung adenocarcinoma.

Entering the spotlight: Chitosan oligosaccharides as novel activators of CaCCs/TMEM16A.

Calcium-activated chloride channels (CaCCs)/TMEM16A control diverse fundamental physiological functions, and abnormal function of TMEM16A will lead to various diseases including asthma, hypertension, gastrointestinal hypomotility and cancers. Therefore, TMEM16A as drug targets for related diseases has been increasingly concerned by researchers. In this work, COS were reported as novel natural activators of TMEM16A. It was demonstrated that COS can activate TMEM16A in a concentration dependent manner, with an EC50 of 74.5 µg/mL. Then, fluorescence experiments and inside-out patch clamp experiments were combined to confirm that COS can directly activate TMEM16A. Further, we compared the activation effects of COS monomers DP2 to DP6, with DP3 the best activator. Molecular simulation was performed to find that the binding sites between DP3 and TMEM16A are E143 and E146 in TMEM16A, and it was speculated that COS and TMEM16A may be combined by electrostatic interaction. Finally, we verified that guinea pig ileum contraction was promoted by COS and the monomers through activating TMEM16A. Collectively, COS are novel efficient natural activators of TMEM16A, with potential to be developed to treatment diseases caused by down-regulation of TMEM16A including gastrointestinal hypomotility.

Ginsenoside Rb1, a novel activator of the TMEM16A chloride channel, augments the contraction of guinea pig ileum.

Calcium-activated chloride channels (CaCCs) play important roles in many physiological processes, and the molecular basis of CaCCs has been identified as TMEM16A in many cell types. It is well established that TMEM16A is a drug target in many diseases, including cystic fibrosis, hypertension, asthma, and various tumors. Therefore, identifying potent and specific modulators of the TMEM16A channel is crucial. In this study, we identified the first natural activator of TMEM16A from traditional Chinese medicine and explored its mechanism. Our data showed that Ginsenoside Rb1 (GRb1) can activate TMEM16A directly from the intracellular side in a dose-dependent manner at an EC50 of 38.4 ± 2.14 µM. GRb1 specifically activated TMEM16A/B, but not the other previously proposed CaCC mediators such as CFTR and bestrophin. Moreover, GRb1 promoted proliferation of CHO cells stably expressing TMEM16A, in a concentration-dependent manner. Finally, we showed that GRb1 increased the amplitude and frequency of contractions in an isolated guinea pig ileum assay in vivo. In summary, GRb1 can be considered a lead compound for the development of novel drugs for the treatment of diseases caused by TMEM16A dysfunction.

Hydrocinnamic Acid Inhibits the Currents of WT and SQT3 Syndrome-Related Mutants of Kir2.1 Channel.

Gain of function in mutations, D172N and E299V, of Kir2.1 will induce type III short QT syndrome. In our previous work, we had identified that a mixture of traditional Chinese medicine, styrax, is a blocker of Kir2.1. Here, we determined a monomer, hydrocinnamic acid (HA), as the effective component from 18 compounds of styrax. Our data show that HA can inhibit the currents of Kir2.1 channel in both excised inside-out and whole-cell patch with the IC50 of 5.21 ± 1.02 and 10.08 ± 0.46 mM, respectively. The time course of HA blockage and washout are 2.3 ± 0.6 and 10.5 ± 2.6 s in the excised inside-out patch. Moreover, HA can also abolish the currents of D172N and E299V with the IC50 of 6.66 ± 0.57 and 5.81 ± 1.10 mM for D172N and E299V, respectively. Molecular docking results determine that HA binds with Kir2.1 at K182, K185, and K188, which are phosphatidylinositol 4,5-bisphosphate (PIP2) binding residues. Our results indicate that HA competes with PIP2 to bind with Kir2.1 and inhibits the currents.

Three pairs of weak interactions precisely regulate the G-loop gate of Kir2.1 channel.

Kir2.1 (also known as IRK1) plays key roles in regulation of resting membrane potential and cell excitability. To achieve its physiological roles, Kir2.1 performs a series of conformational transition, named as gating. However, the structural basis of gating is still obscure. Here, we combined site-directed mutation, two-electrode voltage clamp with molecular dynamics simulations and determined that H221 regulates the gating process of Kir2.1 by involving a weak interaction network. Our data show that the H221R mutant accelerates the rundown kinetics and decelerates the reactivation kinetics of Kir2.1. Compared with the WT channel, the H221R mutation strengthens the interaction between the CD- and G-loops (E303-R221) which stabilizes the close state of the G-loop gate and weakens the interactions between C-linker and CD-loop (R221-R189) and the adjacent G-loops (E303-R312) which destabilizes the open state of G-loop gate. Our data indicate that the three pairs of interactions (E303-H221, H221-R189 and E303-R312) precisely regulate the G-loop gate by controlling the conformation of G-loop. Proteins 2016; 84:1929-1937. © 2016 Wiley Periodicals, Inc.

Molecular simulation assisted identification of Ca(2+) binding residues in TMEM16A.

Calcium-activated chloride channels (CaCCs) play vital roles in a variety of physiological processes. Transmembrane protein 16A (TMEM16A) has been confirmed as the molecular counterpart of CaCCs which greatly pushes the molecular insights of CaCCs forward. However, the detailed mechanism of Ca(2+) binding and activating the channel is still obscure. Here, we utilized a combination of computational and electrophysiological approaches to discern the molecular mechanism by which Ca(2+) regulates the gating of TMEM16A channels. The simulation results show that the first intracellular loop serves as a Ca(2+) binding site including D439, E444 and E447. The experimental results indicate that a novel residue, E447, plays key role in Ca(2+) binding. Compared with WT TMEM16A, E447Y produces a 30-fold increase in EC50 of Ca(2+) activation and leads to a 100-fold increase in Ca(2+) concentrations that is needed to fully activate the channel. The following steered molecular dynamic (SMD) simulation data suggests that the mutations at 447 reduce the Ca(2+) dissociation energy. Our results indicated that both the electrical property and the size of the side-chain at residue 447 have significant effects on Ca(2+) dependent gating of TMEM16A.

A novel biophysical model on calcium and voltage dual dependent gating of calcium-activated chloride channel.

Ca(2+)-activated Cl(-) channels (CaCCs) are anion-selective channels and involved in physiological processes such as electrolyte/fluid secretion, smooth muscle excitability, and olfactory perception which critically depend on the Ca(2+) and voltage dual-dependent gating of channels. However, how the Ca(2+) and voltage regulate the gating of CaCCs still unclear. In this work, the authors constructed a biophysical model to illustrate the dual-dependent gating of CaCCs. For validation, we applied our model on both native CaCCs and exogenous TMEM16A which is thought to be the molecular basis of CaCCs. Our data show that the native CaCCs may share universal gating mechanism. We confirmed the assumption that by binding with the channel, Ca(2+) decreases the energy-barrier to open the channel, but not changes the voltage-sensitivity. For TMEM16A, our model indicates that the exogenous channels show different Ca(2+) dependent gating mechanism from the native ones. These results advance the understanding of intracellular Ca(2+) and membrane potential regulation in CaCCs, and shed new light on its function in aspect of physiology and pharmacology.

Combining fragment homology modeling with molecular dynamics aims at prediction of Ca²âº binding sites in CaBPs.

The family of calcium-binding proteins (CaBPs) consists of dozens of members and contributes to all aspects of the cell's function, from homeostasis to learning and memory. However, the Ca²âº-binding mechanism is still unclear for most of CaBPs. To identify the Ca²âº-binding sites of CaBPs, this study presented a computational approach which combined the fragment homology modeling with molecular dynamics simulation. For validation, we performed a two-step strategy as follows: first, the approach is used to identify the Ca²âº-binding sites of CaBPs, which have the EF-hand Ca²âº-binding site and the detailed binding mechanism. To accomplish this, eighteen crystal structures of CaBPs with 49 Ca²âº-binding sites are selected to be analyzed including calmodulin. The computational method identified 43 from 49 Ca²âº-binding sites. Second, we performed the approach to large-conductance Ca²âº-activated K⁺ (BK) channels which don't have clear Ca²âº-binding mechanism. The simulated results are consistent with the experimental data. The computational approach may shed some light on the identification of Ca²âº-binding sites in CaBPs.

Molecular dynamics simulation of TMEM16A channel: Linking structure with gating.

TMEM16A, the calcium-activated chloride channel, is broadly expressed and plays pivotal roles in diverse physiological processes. To understand the structural and functional relationships of TMEM16A, it is necessary to fully clarify the structural basis of the gating of the TMEM16A channel. Herein, we performed the protein electrostatic analysis and molecular dynamics simulation on the TMEM16A in the presence and absence of Ca2+. Data showed that the separation of TM4 and TM6 causes pore expansion, and Q646 may be a key residue for the formation of π-helix in the middle segment of TM6. Moreover, E705 was found to form a group of H-bond interactions with D554/K588/K645 below the hydrophobic gate to stabilize the closed conformation of the pore in the Ca2+-free state. Interestingly, in the Ca2+ bound state, the E705 side chain swings 100o to serve as Ca2+-binding coordination and released K645. K645 is closer to the hydrophobic gate in the calcium-bound state, which facilitates the provision of electrostatic forces for chloride ions as the ions pass through the hydrophobic gate. Our findings provide the structural-based insights to understanding the mechanisms of gating of TMEM16A.

TMEM16A Protein: Calcium-Binding Site and its Activation Mechanism.

TMEM16A mediates the calcium-activated transmembrane flow of chloride ions and a variety of physiological functions. The binding of cytoplasmic calcium ions of TMEM16A and the consequent conformational changes of it are the key issues to explore the structure-function relationship. In recent years, researchers have explored this issue through electrophysiological experiments, structure resolving, molecular dynamic simulation, and other methods. The structures of TMEM16 family members determined by cryo-Electron microscopy (cryo-EM) and X-ray crystallization provide the primary basis for the investigation of the molecular mechanism of TMEM16A. However, the binding and activation mechanism of calcium ions in TMEM16A are still unclear and controversial. This mini-review discusses four Ca2+ sensing sites of TMEM16A and analyzes activation properties of TMEM16A by them, which will help understand the structure-function relationship of TMEM16A and throw light on the molecular design targeting the TMEM16A channel.

Recent progress in structural studies on TMEM16A channel.

The calcium-activated chloride channel, also known as TMEM16A, shows both calcium and membrane potential dependent activation. The channel is expressed broadly and contributes to a variety of physiological processes, and it is expected to be a target for the treatment of diseases such as hypertension, pain, cystic fibrosis and lung cancer. A thorough understanding of the structural characteristics of TMEM16A is important to reveal its physiological and pathological roles. Recent studies have released several Cryo-EM structures of the channel, revealed the structural basis and mechanism of the gating of the channel. This review focused on the understandings of the structural basis and molecular mechanism of the gating and permeation of TMEM16A channel, which will provide important basis for the development of drugs targeting TMEM16A.

Styrax blocks inward and outward current of Kir2.1 channel.

Kir2.1 plays key roles in setting rest membrane potential and modulation of cell excitability. Mutations of Kir2.1, such as D172N or E299V, inducing gain-of-function, can cause type3 short QT syndrome (SQT3) due to the enlarged outward currents. So far, there is no clinical drug target to block the currents of Kir2.1. Here, we identified a novel blocker of Kir2.1, styrax, which is a kind of natural compound selected from traditional Chinese medicine. Our data show that styrax can abolish the inward and outward currents of Kir2.1. The IC50 of styrax for WT, D172N and E299V are 0.0113 ± 0.00075, 0.0204 ± 0.0048 and 0.0122 ± 0.0012 (in volume), respectively. The results indicate that styrax can serve as a novel blocker for Kir2.1.

Allosteric-activation mechanism of BK channel gating ring triggered by calcium ions.

Calcium ions bind at the gating ring which triggers the gating of BK channels. However, the allosteric mechanism by which Ca2+ regulates the gating of BK channels remains obscure. Here, we applied Molecular Dynamics (MD) and Targeted MD to the integrated gating ring of BK channels, and achieved the transition from the closed state to a half-open state. Our date show that the distances of the diagonal subunits increase from 41.0 Å at closed state to 45.7Å or 46.4 Å at a half-open state. It is the rotatory motion and flower-opening like motion of the gating rings which are thought to pull the bundle crossing gate to open ultimately. Compared with the 'Ca2+ bowl' at RCK2, the RCK1 Ca2+ sites make more contribution to opening the channel. The allosteric motions of the gating ring are regulated by three group of interactions. The first weakened group is thought to stabilize the close state; the second strengthened group is thought to stabilize the open state; the third group thought to lead AC region forming the CTD pore to coordinated motion, which exquisitely regulates the conformational changes during the opening of BK channels by Ca2+.

Identification of the Conformational transition pathway in PIP2 Opening Kir Channels.

The gating of Kir channels depends critically on phosphatidylinositol 4,5-bisphosphate (PIP2), but the detailed mechanism by which PIP2 regulates Kir channels remains obscure. Here, we performed a series of Targeted molecular dynamics simulations on the full-length Kir2.1 channel and, for the first time, were able to achieve the transition from the closed to the open state. Our data show that with the upward motion of the cytoplasmic domain (CTD) the structure of the C-Linker changes from a loop to a helix. The twisting of the C-linker triggers the rotation of the CTD, which induces a small downward movement of the CTD and an upward motion of the slide helix toward the membrane that pulls the inner helix gate open. At the same time, the rotation of the CTD breaks the interaction between the CD- and G-loops thus releasing the G-loop. The G-loop then bounces away from the CD-loop, which leads to the opening of the G-loop gate and the full opening of the pore. We identified a series of interaction networks, between the N-terminus, CD loop, C linker and G loop one by one, which exquisitely regulates the global conformational changes during the opening of Kir channels by PIP2.

TMEM16A/B associated CaCC: structural and functional insights.

Calcium-activated chloride channels (CaCCs) play fundamental roles in numerous physiological processes. Transmembrane proteins 16A and 16B (TMEM16A/B) were identified to be the best molecular identities of CaCCs to date. This makes molecular investigation of CaCCs become possible. This review discusses the latest findings of TMEM16A/B associated CaCCs, the calcium and voltage dual dependence，the reorganization of Ca(2+)-binding site, the mechanisms of direct or indirect activation, the structure-functional relationship, and the possible stereoscopic structure. TMEM16A and other members of the family are associated with several kinds of cancers and other chloride channelopathies. An understanding of TMEM16 associated channel proteins will shed some light on their role in oncology and in pharmacology development.