Challenging AQP4 druggability for NMO-IgG antibody binding using molecular dynamics and molecular interaction fields

ACT001 Relieves NMOSD Symptoms by Reducing Astrocyte Damage with an Autoimmune Antibody

Neuromyelitis optica spectrum disorder (NMOSD) is a central nervous system inflammatory demyelinating disease, the pathogenesis of which involves autoantibodies targeting the extracellular epitopes of aquaporin-4 on astrocytes. We neutralized the AQP4-IgG from NMOSD patient sera using synthesized AQP4 extracellular epitope peptides and found that the severe cytotoxicity produced by aquaporin-4 immunoglobin (AQP4-IgG) could be blocked by AQP4 extracellular mimotope peptides of Loop A and Loop C in astrocyte protection and animal models. ACT001, a natural compound derivative, has shown anti-tumor activity in various cancers. In our study, the central nervous system anti-inflammatory effect of ACT001 was investigated. The results demonstrated the superior astrocyte protection activity of ACT001 at 10 µM. Furthermore, ACT001 decreases the behavioral score in the mouse NMOSD model, which was not inferior to Methylprednisolone Sodium Succinate, the first-line therapy of NMOSD in clinical practice. In summary, our study showed that astrocytes are protected by specific peptides, or small molecular drugs, which is a new strategy for the treatment of NMOSD. It is possible for ACT001 to be a promising therapy for NMOSD.

A Novel Gating Mechanism of Aquaporin-4 Water Channel Mediated by Blast Shockwaves for Brain Edema

As the principal water channel in the brain, aquaporin-4 (AQP4) plays a vital role in brain edema, but its role in blast brain edema is unclear. On the basis of molecular simulations, we reveal the atomically detailed picture of AQP4 in response to blast shockwaves. The results show that the shockwave alone closes the AQP4 channel; however, shock-induced bubble collapse opens it. The jet from bubble collapse forcefully increases the distance between helices and the tilt angles of six helices relative to the membrane vertical direction in a very short time. The average channel size increases about 2.6 times, and the water flux rate is nearly 20 times higher than for normal states. It is responsible for abnormal water transport and a potential cause of acute blast brain edema. Additionally, the open AQP4 channel quickly returns to its normal state, which is in turn helpful for edema absorption. Thus, a novel gating mechanism for AQP4 related to the secondary structure change has been provided, which is different from the previous residue-mediated gating mechanism.

Analysis of Binding Modes of Antigen-Antibody Complexes by Molecular Mechanics Calculation

Antibodies are one of the most important protein molecules in biopharmaceutics. Due to the recent advance in technology for producing monoclonal antibodies, many structural data are available on the antigen-antibody complexes. To characterize the molecular interaction in antigen-antibody recognition, we computationally analyzed 500 complex structures by molecular mechanics calculations. The presence of Ser and Tyr is markedly large in the complementarity-determining regions (CDRs). Although Ser is abundant in CDRs, its contribution to the binding score is not large. Instead, Tyr, Asp, Glu, and Arg significantly contribute to the molecular interaction from the viewpoint of the binding score. The decomposition of the binding score suggests that the hydrophilic interaction is predominant in all CDRs compared with the hydrophobic one. The contribution of the heavy chain is larger than that of the light chain. In particular, H2 and H3 largely contribute to the binding interaction. Tyr is a main contributing residue both in H2 and H3. The positively charged residue Arg also significantly contributes to the binding score in H3, while the contribution of Lys is small. The appearance of Ser is remarkable in H2, and Asp is abundant in H3. The non-charged polar residues, Thr, Asn, and Gln, appear much in H2, compared to appearing in H3. The negatively charged residues Asp and Glu significantly contribute to the binding score in H3. The contributions of Phe and Trp are not large in spite that the aromatic residues are capable of making the π-π or CH-π interaction. Gly is commonly abundant both in H2 and H3. The average distance of the shortest direct hydrogen bond between the antigen and antibody is longer than that of the hydrogen bonds observed in the complexes between compounds and their target proteins. Therefore, the antigen-antibody interface is not so tight as the compound-target protein interface. The calculation of shape complementarity is consistent with the result of the hydrogen bonds in that the fitness of the antigen-antibody contact is not so high as that of the compound-target protein contact. There exist many water molecules at the antigen-antibody interface. These findings suggest that Tyr, Asp, Glu, and Arg are rich in H3 and work as major contributors for the interaction with the antigen. Ser, Thr, Asn, and Gln are rich in H2 and support the interaction with enhancing molecular fitness. Gly is helpful in increasing flexibility and geometrical diversity. Because the antigen-antibody binding is fundamentally hydrophilic-driven, the non-polar residues are unfavorable for mediating the contact even for the aromatic residues such as Phe and Trp.

Exploring the role of elongation Factor-Like 1 (EFL1) in Shwachman-Diamond syndrome through molecular dynamics

Shwachman-Diamond Syndrome (SDS) is an autosomal recessive disorder whose patients present mutations in two ribosome assembly proteins, the Shwachman-Bodian-Diamond Syndrome protein (SBDS) and the Elongation Factor-Like 1 (EFL1). Due to the lack of knowledge of the molecular mechanisms responsible for SDS pathogenesis, current therapy is nonspecific and focuses only at alleviating the symptoms. Building on the recent observation that EFL1 single-point mutations clinically manifest as SDS-like phenotype, we carried out comparative Molecular Dynamics (MD) simulations on three mutants, T127A, M882K and R1095Q and wild type EFL1. As supported by small angle X-ray scattering experiments, the obtained data improve the static EFL1 model resulting from the Cryo-electron microscopy and clearly show that all the mutants experience a peculiar rotation, around the hinge region, of domain IV with respect to domains I and II leading to a different conformation respect to that of wild type protein. This study supports the notion that EFL1 function is governed by an allosteric mechanism involving the concerted action of GTPase domain (domain I) and the domain IV and can help point towards new approaches to SDS treatment.Communicated by Ramaswamy H. Sarma.

Nitazoxanide inhibits paramyxovirus replication by targeting the Fusion protein folding: role of glycoprotein-specific thiol oxidoreductase ERp57

Paramyxoviridae, a large family of enveloped viruses harboring a nonsegmented negative-sense RNA genome, include important human pathogens as measles, mumps, respiratory syncytial virus (RSV), parainfluenza viruses, and henipaviruses, which cause some of the deadliest emerging zoonoses. There is no effective antiviral chemotherapy for most of these pathogens. Paramyxoviruses evolved a sophisticated membrane-fusion machine consisting of receptor-binding proteins and the fusion F-protein, critical for virus infectivity. Herein we identify the antiprotozoal/antimicrobial nitazoxanide as a potential anti-paramyxovirus drug targeting the F-protein. We show that nitazoxanide and its circulating-metabolite tizoxanide act at post-entry level by provoking Sendai virus and RSV F-protein aggregate formation, halting F-trafficking to the host plasma membrane. F-protein folding depends on ER-resident glycoprotein-specific thiol-oxidoreductase ERp57 for correct disulfide-bond architecture. We found that tizoxanide behaves as an ERp57 non-competitive inhibitor; the putative drug binding-site was located at the ERp57-b/b′ non-catalytic domains interface. ERp57-silencing mimicked thiazolide-induced F-protein alterations, suggesting an important role of this foldase in thiazolides anti-paramyxovirus activity. Nitazoxanide is used in the clinic as a safe and effective antiprotozoal/antimicrobial drug; its antiviral activity was shown in patients infected with hepatitis-C virus, rotavirus and influenza viruses. Our results now suggest that nitazoxanide may be effective also against paramyxovirus infection.

Strategies of Virtual Screening in Medicinal Chemistry

Virtual screening represents an effective computational strategy to rise-up the chances of finding new bioactive compounds by accelerating the time needed to move from an initial intuition to market. Classically, the most pursued approaches rely on ligand- and structure-based studies, the former employed when structural data information about the target is missing while the latter employed when X-ray/NMR solved or homology models are instead available for the target. The authors will focus on the most advanced techniques applied in this area. In particular, they will survey the key concepts of virtual screening by discussing how to properly select chemical libraries, how to make database curation, how to applying and- and structure-based techniques, how to wisely use post-processing methods. Emphasis will be also given to the most meaningful databases used in VS protocols. For the ease of discussion several examples will be presented.

Regulation of Aquaporin Functional Properties Mediated by the Antioxidant Effects of Natural Compounds

Some aquaporins (AQPs) have been recently demonstrated to facilitate the diffusion of hydrogen peroxide (H₂O₂) from the producing cells to the extracellular fluid, and their reactive oxygen species scavenging properties have been defined. Nevertheless, the identification of different AQPs acting as peroxiporins, their functional role in eustress and distress, and the identification of antioxidant compounds able to regulate AQP gating, remain unsolved. This study aims to investigate, in HeLa cells: (1) the expression of different AQPs; (2) the evaluation of naringenin, quercetin, (R)-aloesaponol III 8-methyl ether, marrubiin, and curcumin antioxidant profiles, via α,α-diphenyl-β-picrylhydrazyl assay; (3) the effect of the compounds on the water permeability in the presence and in the absence of oxidative stress; and (4) the effect of pre- and post-treatment with the compounds on the H₂O₂ content in heat-stressed cells. Results showed that HeLa cells expressed AQP1, 3, 8, and 11 proteins. The oxidative stress reduced the water transport, and both pre- and post-treatment with the natural compounds recovering the water permeability, with the exception of curcumin. Moreover, the pre- and post-treatment with all the compounds reduced the H₂O₂ content of heat-stressed cells. This study confirms that oxidative stress reduced water AQP-mediated permeability, reversed by some chemical antioxidant compounds. Moreover, curcumin was shown to regulate AQP gating. This suggests a novel mechanism to regulate cell signaling and survival during stress, and to manipulate key signaling pathways in cancer and degenerative diseases.

The aquaporin-4 water channel as a potential drug target in neurological disorders

INTRODUCTION

Aquaporin-4 (AQP4) is a water transporting protein expressed at the plasma membrane of astrocytes throughout the central nervous system (CNS). Analysis of AQP4 knockout mice has suggested its broad involvement in brain water balance, neuroexcitation, glial scarring, neuroinflammation, and even neurodegenerative and neuropsychiatric disorders. Broad clinical utility of AQP4 modulators has been speculated. Area covered: This review covers the biology of AQP4, evidence for its roles in normal CNS function and neurological disorders, and progress in AQP4 drug discovery. Expert opinion: Critical examination of available data reduces the lengthy potential applications list to AQP4 inhibitors for early therapy of ischemic stroke and perhaps for reduction of glial scarring following CNS injury. Major challenges in identification and clinical development of AQP4 inhibitors include the apparent poor druggability of AQPs, the many homologous AQP isoforms with broad tissue distribution and functions, technical issues with water transport assays, predicted undesired CNS and non-CNS actions, and the need for high blood-brain barrier permeation. To date, despite considerable effort, validated small-molecule AQP4 inhibitors have not been advanced. However, a biologic ('aquaporumab') is in development for neuromyelitis optica, an autoimmune inflammatory demyelinating disease where CNS pathology is initiated by binding of anti-AQP4 autoantibodies to astrocyte AQP4.

Aquaporin-Targeted Therapeutics: State-of-the-Field

Drugs targeting aquaporins have broad potential clinical applications, including cancer, obesity, edema, glaucoma, skin diseases and others. The astrocyte water channel aquaporin-4 is a particularly compelling target because of its role of brain water movement, neuroexcitation and glia scarring, and because it is the target of pathogenic autoantibodies in the neuroinflammatory demyelinating disease neuromyelitis optica . There has been considerable interest in the identification of small molecule inhibitors of aquaporins, with various candidates emerging from testing of known ion transport inhibitors, as well as compound screening and computational chemistry. However, in general, the activity of reported aquaporin inhibitors has not been confirmed on retesting, which may be due to technical problems in water transport assays used in the original identification studies, and the challenges in modulating the activity of small, compact, pore-containing membrane proteins. We review here the state of the field of aquaporin-modulating small molecules and biologics, and the challenges and opportunities in moving forward.

Novel chemotypes targeting tubulin at the colchicine binding site and unbiasing P-glycoprotein

Retrospective validation studies carried out on three benchmark databases containing a small fraction (that is 2.80%) of known tubulin binders permitted us to develop a computational platform very effective in selecting easier manageable subsets showing by far higher percentages of actives (about 25%). These studies relied on the hierarchical application of multilayer in silico screenings employing filters implying molecular shape similarity; a structure-based pharmacophore model and molecular docking campaigns. Building on this validated approach, we performed intensive prospective studies to screen a large chemical collection, including up to 3.7 millions of commercial compounds, to across an unexplored and patent space in the search of novel colchicine binding site inhibitors. Our investigation was successful in identifying a pool of 31 initial hits showing new molecular scaffolds (such as 4,5-dihydro-1H-pyrrolo[3,4-c]pyrazol-6-one and pyrazolo[1,5-a]pyrimidine). This panel of new hits resulted antiproliferative activity in the low μM range towards MCF-7 human breast cancer, HepG2 human liver cancer, HeLa human ovarian cancer and SHSY5Y human glioblastoma cell lines as well as interesting concentration-dependent inhibition of tubulin polymerization assessed through fluorescence polymerization assays. Unlike typical tubulin inhibitors, a satisfactorily low sensitivity towards P-gp was also measured in bi-directional transport studies across MDCKII-MDR1 cells for a selected subset of seven compounds.

Discovery of Novel Ligands for TNF-α and TNF Receptor-1 through Structure-Based Virtual Screening and Biological Assay

Tumor necrosis factor α (TNF-α) is overexpressed in various diseases, and it has been a validated therapeutic target for autoimmune diseases. All therapeutics currently used to target TNF-α are biomacromolecules, and limited numbers of TNF-α chemical inhibitors have been reported, which makes the identification of small-molecule alternatives an urgent need. Recent studies have mainly focused on identifying small molecules that directly bind to TNF-α or TNF receptor-1 (TNFR1), inhibit the interaction between TNF-α and TNFR1, and/or regulate related signaling pathways. In this study, we combined in silico methods with biophysical and cell-based assays to identify novel antagonists that bind to TNF-α or TNFR1. Pharmacophore model filtering and molecular docking were applied to identify potential TNF-α antagonists. In regard to TNFR1, we constructed a three-dimensional model of the TNF-α-TNFR1 complex and carried out molecular dynamics simulations to sample the conformations. The residues in TNF-α that have been reported to play important roles in the TNF-α-TNFR1 complex were removed to form a pocket for further virtual screening of TNFR1-binding ligands. We obtained 20 virtual hits and tested them using surface plasmon resonance-based assays, which resulted in one ligand that binds to TNFR1 and four ligands with different scaffolds that bind to TNF-α. T1 and R1, the two most active compounds with Kd values of 11 and 16 μM for TNF-α and TNFR1, respectively, showed activities similar to those of known antagonists. Further cell-based assays also demonstrated that T1 and R1 have similar activities compared to the known TNF-α antagonist C87. Our work has not only produced several TNF-α and TNFR1 antagonists with novel scaffolds for further structural optimization but also showcases the power of our in silico methods for TNF-α- and TNFR1-based drug discovery.

How human IgGs against myelin basic protein (MBP) recognize oligopeptides and MBP

Myelin basic protein (MBP) is a major protein of myelin-proteolipid shell of axons, and it plays an important role in pathogenesis of multiple sclerosis. In the literature, there are no data on how antibodies recognize different protein antigens including MBP. A stepwise increase in ligand complexity was used to estimate the relative contributions of virtually every amino acid residue (AA) of a specific 12-mer LSRFSWGAEGQK oligopeptide corresponding to immunodominant sequence of MBP to the light chains and to intact anti-MBP IgGs from sera of patients with multiple sclerosis. It was shown that the minimal ligands of the light chains of IgGs are many different free AAs (K_d  = 0.51-0.016 M), and each free AA interacts with the specific subsite of the light chain intended for recognition of this AA in specific LSRFSW oligopeptide. A gradual transition from Leu to LSRFSWGAEGQK leads to an increase in the affinity from 10^-1 to 2.3 × 10^-4  M because of additive interactions of the light chain with 6 AAs of this oligopeptide and then the affinity reaches plateau. The contributions of 6 various AAs to the affinity of the oligopeptide are different (K_d , M): 0.71 (S), 0.44 (R), 0.14 (F), 0.17 (S), and 0.62 (W). Affinity of nonspecific oligopeptides to the light chains of IgGs is significantly lower. Intact MBP interacts with both light and heavy chains of IgGs demonstrating 192-fold higher affinity than the specific oligopeptide. It is a first quantitative analysis of the mechanism of proteins recognition by antibodies. The thermodynamic model was constructed to describe the interactions of IgGs with MBP. The data obtained can be very useful for understanding how antibodies against many different proteins can recognize these proteins.

Comparative molecular dynamics study of neuromyelitis optica-immunoglobulin G binding to aquaporin-4 extracellular domains

Neuromyelitis optica (NMO) is an inflammatory demyelinating disease of the central nervous system in which most patients have serum autoantibodies (called NMO-IgG) that bind to astrocyte water channel aquaporin-4 (AQP4). A potential therapeutic strategy in NMO is to block the interaction of NMO-IgG with AQP4. Building on recent observation that some single-point and compound mutations of the AQP4 extracellular loop C prevent NMO-IgG binding, we carried out comparative Molecular Dynamics (MD) investigations on three AQP4 mutants, TP^137-138AA, N¹⁵³Q and V¹⁵⁰G, whose 295-ns long trajectories were compared to that of wild type human AQP4. A robust conclusion of our modeling is that loop C mutations affect the conformation of neighboring extracellular loop A, thereby interfering with NMO-IgG binding. Analysis of individual mutations suggested specific hydrogen bonding and other molecular interactions involved in AQP4-IgG binding to AQP4.

Computing osmotic permeabilities of aquaporins AQP4, AQP5, and GlpF from near-equilibrium simulations

Measuring or computing the single-channel permeability of aquaporins/aquaglyceroporins (AQPs) has long been a challenge. The measured values scatter over an order of magnitude but the corresponding Arrhenius activation energies converge in the current literature. Osmotic flux through an AQP was simulated as water current forced through the channel by kilobar hydraulic pressure or theoretically approximated as single-file diffusion. In this paper, we report large scale simulations of osmotic current under sub M gradient through three AQPs (water channels AQP4 and AQP5 and glycerol-water channel GlpF) using the mature particle mesh Ewald technique (PME) for which the established force fields have been optimized with known accuracy. These simulations were implemented with hybrid periodic boundary conditions devised to avoid the artifactitious mixing across the membrane in a regular PME simulation. The computed single-channel permeabilities at 5°C and 25°C are in agreement with recently refined experiments on GlpF. The Arrhenius activation energies extracted from our simulations for all the three AQPs agree with the in vitro measurements. The single-file diffusion approximations from our large-scale simulations are consistent with the current literature on smaller systems. From these unambiguous agreements among the in vitro and in silico studies, we observe the quantitative accuracy of the all-atom force fields of the current literature for water-channel biology. We also observe that AQP4, that is particularly rich in the central nervous system, is more efficient in water conduction and more temperature-sensitive than other water-only channels (excluding glycerol channels that also conduct water when not inhibited by glycerol).

Role of the H-bond between L53 and T56 for Aquaporin-4 epitope in Neuromyelitis Optica

Aquaporin-4 (AQP4) is the CNS water channel organized into well-ordered protein aggregates called Orthogonal Arrays of Particles (OAPs). Neuromyelitis Optica (NMO) is an autoimmune disease caused by anti-OAP autoantibodies (AQP4-IgG). Molecular Dynamics (MD) simulations have identified an H-bond between L53 and T56 as the key for AQP4 epitope and therefore of potential interest for drug design in NMO field. In the present study, we have experimentally tested this MD-prediction using the classic mutagenesis approach. We substituted T56 with V56 and tested this mutant for AQP4 aggregates and AQP4-IgG binding. gSTED super-resolution microscopy showed that the mutation does not affect AQP4 aggregate dimension; immunofluorescence and cytofluorimetric analysis demonstrated its unaltered AQP4-IgG binding, therefore invalidating the MD-prediction. We later investigated whether AQP4, expressed in Sf9 insect and HEK-293F cells, is able to correctly aggregate before and after the purification steps usually applied to obtain AQP4 crystal. The results demonstrated that AQP4-IgG recognizes AQP4 expressed in Sf9 and HEK-293F cells by immunofluorescence even though BN-PAGE analysis showed that AQP4 forms smaller aggregates when expressed in insect cells compared to mammalian cell lines. Notably, after AQP4 purification, from both insect and HEK-293F cells, no aggregates are detectable by BN-PAGE and AQP4-IgG binding is impaired in sandwich ELISA assays. All together these results indicate that 1) the MD prediction under analysis is not supported by experimental data and 2) the procedure to obtain AQP4 crystals might affect its native architecture and, as a consequence, MD simulations. In conclusion, given the complex nature of the AQP4 epitope, MD might not be the suitable for molecular medicine advances in NMO.

Effects of Mutations on Structure-Function Relationships of Matrix Metalloproteinase-1

Matrix metalloproteinase-1 (MMP-1) is one of the most widely studied enzymes involved in collagen degradation. Mutations of specific residues in the MMP-1 hemopexin-like (HPX) domain have been shown to modulate activity of the MMP-1 catalytic (CAT) domain. In order to reveal the structural and conformational effects of such mutations, a molecular dynamics (MD) study was performed of in silico mutated residues in the X-ray crystallographic structure of MMP-1 complexed with a collagen-model triple-helical peptide (THP). The results indicate an important role of the mutated residues in MMP-1 interactions with the THP and communication between the CAT and the HPX domains. Each mutation has a distinct impact on the correlated motions in the MMP-1•THP. An increased collagenase activity corresponded to the appearance of a unique anti-correlated motion and decreased correlated motions, while decreased collagenase activity corresponded both to increased and decreased anti-correlated motions.

Human Aquaporin-4 and Molecular Modeling: Historical Perspective and View to the Future

Among the different aquaporins (AQPs), human aquaporin-4 (hAQP4) has attracted the greatest interest in recent years as a new promising therapeutic target. Such a membrane protein is, in fact, involved in a multiple sclerosis-like immunopathology called Neuromyelitis Optica (NMO) and in several disorders resulting from imbalanced water homeostasis such as deafness and cerebral edema. The gap of knowledge in its functioning and dynamics at the atomistic level of detail has hindered the development of rational strategies for designing hAQP4 modulators. The application, lately, of molecular modeling has proved able to fill this gap providing a breeding ground to rationally address compounds targeting hAQP4. In this review, we give an overview of the important advances obtained in this field through the application of Molecular Dynamics (MD) and other complementary modeling techniques. The case studies presented herein are discussed with the aim of providing important clues for computational chemists and biophysicists interested in this field and looking for new challenges.

Multidisciplinary study of a new ClC-1 mutation causing myotonia congenita: a paradigm to understand and treat ion channelopathies

Myotonia congenita is an inherited disease that is characterized by impaired muscle relaxation after contraction caused by loss-of-function mutations in the skeletal muscle ClC-1 channel. We report a novel ClC-1 mutation, T335N, that is associated with a mild phenotype in 1 patient, located in the extracellular I-J loop. The purpose of this study was to provide a solid correlation between T335N dysfunction and clinical symptoms in the affected patient as well as to offer hints for drug development. Our multidisciplinary approach includes patch-clamp electrophysiology on T335N and ClC-1 wild-type channels expressed in tsA201 cells, Western blot and quantitative PCR analyses on muscle biopsies from patient and unaffected individuals, and molecular dynamics simulations using a homology model of the ClC-1 dimer. T335N channels display reduced chloride currents as a result of gating alterations rather than altered surface expression. Molecular dynamics simulations suggest that the I-J loop might be involved in conformational changes that occur at the dimer interface, thus affecting gating. Finally, the gene expression profile of T335N carrier showed a diverse expression of K⁺ channel genes, compared with control individuals, as potentially contributing to the phenotype. This experimental paradigm satisfactorily explained myotonia in the patient. Furthermore, it could be relevant to the study and therapy of any channelopathy.-Imbrici, P., Altamura, C., Camerino, G. M., Mangiatordi, G. F., Conte, E., Maggi, L., Brugnoni, R., Musaraj, K., Caloiero, R., Alberga, D., Marsano, R. M., Ricci, G., Siciliano, G., Nicolotti, O., Mora, M., Bernasconi, P., Desaphy, J.-F., Mantegazza, R., Camerino, D. C. Multidisciplinary study of a new ClC-1 mutation causing myotonia congenita: a paradigm to understand and treat ion channelopathies.

ClC-1 mutations in myotonia congenita patients: insights into molecular gating mechanisms and genotype-phenotype correlation

KEY POINTS

Loss-of-function mutations of the skeletal muscle ClC-1 channel cause myotonia congenita with variable phenotypes. Using patch clamp we show that F484L, located in the conducting pore, probably induces mild dominant myotonia by right-shifting the slow gating of ClC-1 channel, without exerting a dominant-negative effect on the wild-type (WT) subunit. Molecular dynamics simulations suggest that F484L affects the slow gate by increasing the frequency and the stability of H-bond formation between E232 in helix F and Y578 in helix R. Three other myotonic ClC-1 mutations are shown to produce distinct effects on channel function: L198P shifts the slow gate to positive potentials, V640G reduces channel activity, while L628P displays a WT-like behaviour (electrophysiology data only). Our results provide novel insight into the molecular mechanisms underlying normal and altered ClC-1 function.

ABSTRACT

Myotonia congenita is an inherited disease caused by loss-of-function mutations of the skeletal muscle ClC-1 chloride channel, characterized by impaired muscle relaxation after contraction and stiffness. In the present study, we provided an in-depth characterization of F484L, a mutation previously identified in dominant myotonia, in order to define the genotype-phenotype correlation, and to elucidate the contribution of this pore residue to the mechanisms of ClC-1 gating. Patch-clamp recordings showed that F484L reduced chloride currents at every tested potential and dramatically right-shifted the voltage dependence of slow gating, thus contributing to the mild clinical phenotype of affected heterozygote carriers. Unlike dominant mutations located at the dimer interface, no dominant-negative effect was observed when F484L mutant subunits were co-expressed with wild type. Molecular dynamics simulations further revealed that F484L affected the slow gate by increasing the frequency and stability of the H-bond formation between the pore residue E232 and the R helix residue Y578. In addition, using patch-clamp electrophysiology, we characterized three other myotonic ClC-1 mutations. We proved that the dominant L198P mutation in the channel pore also right-shifted the voltage dependence of slow gating, recapitulating mild myotonia. The recessive V640G mutant drastically reduced channel function, which probably accounts for myotonia. In contrast, the recessive L628P mutant produced currents very similar to wild type, suggesting that the occurrence of the compound truncating mutation (Q812X) or other muscle-specific mechanisms accounted for the severe symptoms observed in this family. Our results provide novel insight into the molecular mechanisms underlying normal and altered ClC-1 function.