Computational Studies of Molecular Permeation through Connexin26 Channels

An Ala/Glu difference in E1 of Cx26 and Cx30 contributes to their differential anionic permeabilities

Two closely related connexins, Cx26 and Cx30, share widespread expression in the cochlear cellular networks. Gap junction channels formed by these connexins have been shown to have different permeability profiles, with Cx30 showing a strongly reduced preference for anionic tracers. The pore-forming segment of the first extracellular loop, E1, identified by computational studies of the Cx26 crystal structure to form a parahelix and a narrowed region of the pore, differs at a single residue at position 49. Cx26 contains an Ala and Cx30, a charged Glu at this position, and cysteine scanning in hemichannels identified this position to be pore-lining. To assess whether the Ala/Glu difference affects permeability, we modeled and quantified Lucifer Yellow transfer between HeLa cell pairs expressing WT Cx26 and Cx30 and variants that reciprocally substituted Glu and Ala at position 49. Cx26(A49E) and Cx30(E49A) substitutions essentially reversed the Lucifer Yellow permeability profile when accounting for junctional conductance. Moreover, by using a calcein efflux assay in single cells, we observed a similar reduced anionic preference in undocked Cx30 hemichannels and a reversal with reciprocal Ala/Glu substitutions. Thus, our data indicate that Cx26 and Cx30 gap junction channels and undocked hemichannels retain similar permeability characteristics and that a single residue difference in their E1 domains can largely account for their differential permeabilities to anionic tracers. The higher anionic permeability of Cx26 compared with Cx30 suggests that these connexins may serve distinct signaling functions in the cochlea, perhaps reflected in the vastly higher prevalence of Cx26 mutations in human deafness.

Molecular permeation through large pore channels: computational approaches and insights

Computational methods such as molecular dynamics (MD) have illuminated how single-atom ions permeate membrane channels and how selectivity among them is achieved. Much less is understood about molecular permeation through eukaryotic channels that mediate the flux of small molecules (e.g. connexins, pannexins, LRRC8s, CALHMs). Here we describe computational methods that have been profitably employed to explore the movements of molecules through wide pores, revealing mechanistic insights, guiding experiments, and suggesting testable hypotheses. This review illustrates MD techniques such as voltage-driven flux, potential of mean force, and mean first-passage-time calculations, as applied to molecular permeation through wide pores. These techniques have enabled detailed and quantitative modeling of molecular interactions and movement of permeants at the atomic level. We highlight novel contributors to the transit of molecules through these wide pathways. In particular, the flexibility and anisotropic nature of permeant molecules, coupled with the dynamics of pore-lining residues, lead to bespoke permeation dynamics. As more eukaryotic large-pore channel structures and functional data become available, these insights and approaches will be important for understanding the physical principles underlying molecular permeation and as guides for experimental design.

Large-pore connexin hemichannels function like molecule transporters independent of ion conduction

Connexin hemichannels were identified as the first members of the eukaryotic large-pore channel family that mediate permeation of both atomic ions and small molecules between the intracellular and extracellular environments. The conventional view is that their pore is a large passive conduit through which both ions and molecules diffuse in a similar manner. In stark contrast to this notion, we demonstrate that the permeation of ions and of molecules in connexin hemichannels can be uncoupled and differentially regulated. We find that human connexin mutations that produce pathologies and were previously thought to be loss-of-function mutations due to the lack of ionic currents are still capable of mediating the passive transport of molecules with kinetics close to those of wild-type channels. This molecular transport displays saturability in the micromolar range, selectivity, and competitive inhibition, properties that are tuned by specific interactions between the permeating molecules and the N-terminal domain that lies within the pore-a general feature of large-pore channels. We propose that connexin hemichannels and, likely, other large-pore channels, are hybrid channel/transporter-like proteins that might switch between these two modes to promote selective ion conduction or autocrine/paracrine molecular signaling in health and disease processes.

Connexin Gap Junction Channels and Hemichannels: Insights from High-Resolution Structures

Connexins (Cxs) are a family of integral membrane proteins, which function as both hexameric hemichannels (HCs) and dodecameric gap junction channels (GJCs), behaving as conduits for the electrical and molecular communication between cells and between cells and the extracellular environment, respectively. Their proper functioning is crucial for many processes, including development, physiology, and response to disease and trauma. Abnormal GJC and HC communication can lead to numerous pathological states including inflammation, skin diseases, deafness, nervous system disorders, and cardiac arrhythmias. Over the last 15 years, high-resolution X-ray and electron cryomicroscopy (cryoEM) structures for seven Cx isoforms have revealed conservation in the four-helix transmembrane (TM) bundle of each subunit; an αβ fold in the disulfide-bonded extracellular loops and inter-subunit hydrogen bonding across the extracellular gap that mediates end-to-end docking to form a tight seal between hexamers in the GJC. Tissue injury is associated with cellular Ca2+ overload. Surprisingly, the binding of 12 Ca2+ ions in the Cx26 GJC results in a novel electrostatic gating mechanism that blocks cation permeation. In contrast, acidic pH during tissue injury elicits association of the N-terminal (NT) domains that sterically blocks the pore in a "ball-and-chain" fashion. The NT domains under physiologic conditions display multiple conformational states, stabilized by protein-protein and protein-lipid interactions, which may relate to gating mechanisms. The cryoEM maps also revealed putative lipid densities within the pore, intercalated among transmembrane α-helices and between protomers, the functions of which are unknown. For the future, time-resolved cryoEM of isolated Cx channels as well as cryotomography of GJCs and HCs in cells and tissues will yield a deeper insight into the mechanisms for channel regulation. The cytoplasmic loop (CL) and C-terminal (CT) domains are divergent in sequence and length, are likely involved in channel regulation, but are not visualized in the high-resolution X-ray and cryoEM maps presumably due to conformational flexibility. We expect that the integrated use of synergistic physicochemical, spectroscopic, biophysical, and computational methods will reveal conformational dynamics relevant to functional states. We anticipate that such a wealth of results under different pathologic conditions will accelerate drug discovery related to Cx channel modulation.

Molecular Dynamics Simulation of Permeation Through Connexin Channels

Molecular dynamics (MD) simulations are a collection of computational tools that can be used to trace intermolecular interactions at the sub-nanometer level. They offer possibilities that are often unavailable to experimental methods, making MD an ideal complementary technique for the understanding a plethora of biological processes. Thanks to significant efforts by many groups of developers around the world, setting up and running MD simulations has become progressively simpler. However, simulating ionic permeation through membrane channels still presents significant caveats.MD simulations of connexin (Cx) hemichannels (HCs) are particularly problematic because HCs create wide pores in the plasma membrane, and the lateral sizes of the extracellular and intracellular regions are quite different. In this chapter, we provide a detailed instruction to perform MD simulations aimed at computationally modeling the permeation of inorganic ions and larger molecules through Cx HCs.

The Hydrophobic Residues in Amino Terminal Domains of Cx46 and Cx50 Are Important for Their Gap Junction Channel Ion Permeation and Gating

Lens gap junctions (GJs) formed by Cx46 and Cx50 are important to keep lens transparency. Functional studies on Cx46 and Cx50 GJs showed that the V_j-gating, single channel conductance (γ_j), gating polarity, and/or channel open stability could be modified by the charged residues in the amino terminal (NT) domain. The role of hydrophobic residues in the NT on GJ properties is not clear. Crystal and cryo-EM GJ structures have been resolved, but the NT domain structure has either not been resolved or has showed very different orientations depending on the component connexins and possibly other experimental conditions, making it difficult to understand the structural basis of the NT in V_j-gating and γ_j. Here, we generated missense variants in Cx46 and Cx50 NT domains and studied their properties by recombinant expression and dual whole-cell patch clamp experiments on connexin-deficient N2A cells. The NT variants (Cx46 L10I, N13E, A14V, Q15N, and Cx50 I10L, E13N, V14A, N15Q) were all able to form functional GJs with similar coupling%, except Cx46 N13E, which showed a significantly reduced coupling%. The GJs of Cx46 N13E, A14V and Cx50 E13N, N15Q showed a reduced coupling conductance. V_j-gating of all the variant GJs were similar to the corresponding wild-type GJs except Cx46 L10I. The γ_j of Cx46 N13E, A14V, Cx50 E13N, and N15Q GJs was reduced to 51%, 82%, 87%, and 74%, respectively, as compared to their wild-type γ_js. Structural models of Cx46 L10I and A14V predicted steric clashes between these residues and the TM2 residues, which might be partially responsible for our observed changes in GJ properties. To verify the importance of hydrophobic interactions, we generated a variant, Cx50 S89T, which also shows a steric clash and failed to form a functional GJ. Our experimental results and structure models indicate that hydrophobic interactions between the NT and TM2 domain are important for their V_j-gating, γ_j, and channel open stability in these and possibly other GJs.

A novel voltage-clamp/dye uptake assay reveals saturable transport of molecules through CALHM1 and connexin channels

Large-pore channels permeable to small molecules such as ATP, in addition to atomic ions, are emerging as important regulators in health and disease. Nonetheless, their mechanisms of molecular permeation and selectivity remain mostly unexplored. Combining fluorescence microscopy and electrophysiology, we developed a novel technique that allows kinetic analysis of molecular permeation through connexin and CALHM1 channels in Xenopus oocytes rendered translucent. Using this methodology, we found that (1) molecular flux through these channels saturates at low micromolar concentrations, (2) kinetic parameters of molecular transport are sensitive to modulators of channel gating, (3) molecular transport and ionic currents can be differentially affected by mutation and gating, and (4) N-terminal regions of these channels control transport kinetics and permselectivity. Our methodology allows analysis of how human disease-causing mutations affect kinetic properties and permselectivity of molecular signaling and enables the study of molecular mechanisms, including selectivity and saturability, of molecular transport in other large-pore channels.

Free energy and kinetics of cAMP permeation through connexin26 via applied voltage and milestoning

The connexin family is a diverse group of highly regulated wide-pore channels permeable to biological signaling molecules. Despite the critical roles of connexins in mediating selective molecular signaling in health and disease, the basis of molecular permeation through these pores remains unclear. Here, we report the thermodynamics and kinetics of binding and transport of a second messenger, adenosine-3',5'-cyclophosphate (cAMP), through a connexin26 hemichannel (Cx26). First, inward and outward fluxes of cAMP molecules solvated in KCl solution were obtained from 4 μs of ± 200 mV simulations. These fluxes data yielded a single-channel permeability of cAMP and cAMP/K⁺ permeability ratio consistent with experimentally measured values. The results from voltage simulations were then compared with the potential of mean force (PMF) and the mean first passage times (MFPTs) of a single cAMP without voltage, obtained from a total of 16.5 μs of Voronoi-tessellated Markovian milestoning simulations. Both the voltage simulations and the milestoning simulations revealed two cAMP-binding sites, for which the binding constants K_D and dissociation rates k_off were computed from PMF and MFPTs. The protein dipole inside the pore produces an asymmetric PMF, reflected in unequal cAMP MFPTs in each direction once within the pore. The free energy profiles under opposite voltages were derived from the milestoning PMF and revealed the interplay between voltage and channel polarity on the total free energy. In addition, we show how these factors influence the cAMP dipole vector during permeation, and how cAMP affects the local and nonlocal pore diameter in a position-dependent manner.

Connexin 43 Mutations Lead to Increased Hemichannel Functionality in Skin Disease

Gap junctional channels are specialized components of the cellular membrane that allow the intercellular passage of small metabolites, ions, and second messengers to maintain homeostasis. They are comprised of members of the connexin gene family that encode a wide array of proteins that are expressed in nearly every tissue type. Cx43 is perceived to be the most broadly expressed connexin in humans, with several genetic skin diseases being linked to Cx43 mutations specifically. These mutations, in large, produce a gain of functional hemichannels that contribute to the phenotypes of Erythrokeratoderma Variabilis et Progressiva (EKVP), Palmoplantar Keratodemra Congenital Alopecia-1 (PPKCA1), and others that produce large conductance and increased permselectivity in otherwise quiescent structures. Gaining functional hemichannels can have adverse effects in the skin, inducing apoptosis via Ca²⁺ overload or increased ATP permeability. Here, we review the link between Cx43 and skin disease. We aim to provide insight into the mechanisms regulating the normal and pathophysiological gating of these essential proteins, as well as address current therapeutic strategies. We also demonstrate that transient transfection of neuro-2a (N2a) cells with mutant Cx43 cDNA resulted in increased hemichannel activity compared to wild-type Cx43 and untransfected cells, which is consistent with other studies in the current literature.

The connexin26 human mutation N14K disrupts cytosolic intersubunit interactions and promotes channel opening

The N14K mutation in the N-terminal domain of connexin26 produces gain-of-function hemichannels. Valdez Capuccino et al. reveal that the mutation disrupts intersubunit interactions between the N terminus and the second transmembrane domain–cytoplasmic loop transition.

A group of human mutations within the N-terminal (NT) domain of connexin 26 (Cx26) hemichannels produce aberrant channel activity, which gives rise to deafness and skin disorders, including keratitis-ichthyosis-deafness (KID) syndrome. Structural and functional studies indicate that the NT of connexin hemichannels is folded into the pore, where it plays important roles in permeability and gating. In this study, we explore the molecular basis by which N14K, an NT KID mutant, promotes gain of function. In macroscopic and single-channel recordings, we find that the N14K mutant favors the open conformation of hemichannels, shifts calcium and voltage sensitivity, and slows deactivation kinetics. Multiple copies of MD simulations of WT and N14K hemichannels, followed by the Kolmogorov–Smirnov significance test (KS test) of the distributions of interaction energies, reveal that the N14K mutation significantly disrupts pairwise interactions that occur in WT hemichannels between residue K15 of one subunit and residue E101 of the adjacent subunit (E101 being located at the transition between transmembrane segment 2 [TM2] and the cytoplasmic loop [CL]). Double mutant cycle analysis supports coupling between the NT and the TM2/CL transition in WT hemichannels, which is disrupted in N14K mutant hemichannels. KS tests of the α carbon correlation coefficients calculated over MD trajectories suggest that the effects of the N14K mutation are not confined to the K15–E101 pairs but extend to essentially all pairwise residue correlations between the NT and TM2/CL interface. Together, our data indicate that the N14K mutation increases hemichannel open probability by disrupting interactions between the NT and the TM2/CL region of the adjacent connexin subunit. This suggests that NT–TM2/CL interactions facilitate Cx26 hemichannel closure.

Electrical coupling and its channels

As the physiology of synapses began to be explored in the 1950s, it became clear that electrical communication between neurons could not always be explained by chemical transmission. Instead, careful studies pointed to a direct intercellular pathway of current flow and to the anatomical structure that was (eventually) called the gap junction. The mechanism of intercellular current flow was simple compared with chemical transmission, but the consequences of electrical signaling in excitable tissues were not. With the recognition that channels were a means of passive ion movement across membranes, the character and behavior of gap junction channels came under scrutiny. It became evident that these gated channels mediated intercellular transfer of small molecules as well as atomic ions, thereby mediating chemical, as well as electrical, signaling. Members of the responsible protein family in vertebrates-connexins-were cloned and their channels studied by many of the increasingly biophysical techniques that were being applied to other channels. As described here, much of the evolution of the field, from electrical coupling to channel structure-function, has appeared in the pages of the Journal of General Physiology.

Structure of native lens connexin 46/50 intercellular channels by cryo-EM

Gap junctions establish direct pathways for cell-to-cell communication through the assembly of twelve connexin subunits that form intercellular channels connecting neighbouring cells. Co-assembly of different connexin isoforms produces channels with unique properties and enables communication across cell types. Here we used single-particle cryo-electron microscopy to investigate the structural basis of connexin co-assembly in native lens gap junction channels composed of connexin 46 and connexin 50 (Cx46/50). We provide the first comparative analysis to connexin 26 (Cx26), which-together with computational studies-elucidates key energetic features governing gap junction permselectivity. Cx46/50 adopts an open-state conformation that is distinct from the Cx26 crystal structure, yet it appears to be stabilized by a conserved set of hydrophobic anchoring residues. 'Hot spots' of genetic mutations linked to hereditary cataract formation map to the core structural-functional elements identified in Cx46/50, suggesting explanations for many of the disease-causing effects.

Alterations at Arg76 of human connexin 46, a residue associated with cataract formation, cause loss of gap junction formation but preserve hemichannel function

The connexins are members of a family of integral membrane proteins that form gap junction channels between apposed cells and/or hemichannels across the plasma membranes. The importance of the arginine at position 76 (Arg76) in the structure and/or function of connexin 46 (Cx46) is highlighted by its conservation across the entire connexin family and the occurrence of pathogenic mutations at this (or the corresponding homologous) residue in a number of human diseases. Two mutations at Arg76 in Cx46 are associated with cataracts in humans, highlighting the importance of this residue. We examined the expression levels and macroscopic and single-channel properties of human Cx46 and compared them with those for two pathogenic mutants, namely R76H and R76G. To gain further insight into the role of charge at this position, we generated two additional nonnaturally occurring mutants, R76K (charge conserving) and R76E (charge inverting). We found that, when expressed exogenously in Neuro2a cells, all four mutants formed membrane hemichannels, inducing membrane permeability at levels comparable to those recorded in cells expressing the wild-type Cx46. In contrast, the number of gap-junction plaques and the magnitude of junctional coupling were reduced by all four mutations. To gain further insight into the role of Arg76 in the function of Cx46, we performed homology modeling of Cx46 and in silico mutagenesis of Arg76 to Gly, His, or Glu. Our studies suggest that the loss of interprotomeric interactions has a significant effect on the extracellular domain conformation and dynamics, thus affecting the hemichannel docking required for formation of cell-cell channels.

Substitution impact of highly conserved arginine residue at position 75 in GJB1 gene in association with X-linked Charcot-Marie-tooth disease: A computational study

X-linked Charcot-Marie-Tooth type 1 X (CMTX1) disease is a subtype of Charcot-Marie-Tooth (CMT), which is mainly caused by mutations in the GJB1 gene. It is also known as connexin 32 (Cx32) that leads to Schwann cell abnormalities and peripheral neuropathy. CMTX1 is considered as the second most common form of CMT disease. The aim of this study is to computationally predict the potential impact of different single amino acid substitutions at position 75 of Cx32, from arginine (R) to proline (P), glutamine (Q) and tryptophan (W). This position is known to be highly conserved among the family of connexin. To understand the structural and functional changes due to these single amino acid substitutions, we employed a homology-modeling technique to build the three-dimensional structure models for the native and mutant proteins. The protein structures were further embedded into a POPC lipid bilayer, inserted into a water box, and subjected to molecular dynamics simulation for 50 ns. Our results show that the mutants R75P, R75Q and R75W display variable structural conformation and dynamic behavior compared to the native protein. Our data proves useful in predicting the potential pathogenicity of the mutant proteins and is expected to serve as a platform for drug discovery for patients with CMT.

Defining the factors that affect solute permeation of gap junction channels

This review focuses on the biophysical properties and structure of the pore and vestibule of homotypic gap junction channels as they relate to channel permeability and selectivity. Gap junction channels are unique in their sole role to connect the cytoplasm of two adjacent cells. In general, these channels are considered to be poorly selective, possess open probabilities approximating unity, and exhibit mean open times ranging from milliseconds to seconds. These properties suggest that such channels can function as delivery pathways from cell to cell for solutes that are significantly larger than monovalent ions. We have taken quantitative data from published works concerning unitary conductance, ion flux, and permeability for homotypic connexin 43 (Cx43), Cx40, Cx26, Cx50, and Cx37, and performed a comparative analysis of conductance and/or ion/solute flux versus diffusion coefficient. The analysis of monovalent cation flux portrays the pore as equivalent to an aqueous space where hydrogen bonding and weak interactions with binding sites dominate. For larger solutes, size, shape and charge are also significant components in determining the permeation rate. This article is part of a Special Issue entitled: Gap Junction Proteins edited by Jean Claude Herve.

A structural and functional comparison of gap junction channels composed of connexins and innexins

Methods such as electron microscopy and electrophysiology led to the understanding that gap junctions were dense arrays of channels connecting the intracellular environments within almost all animal tissues. The characteristics of gap junctions were remarkably similar in preparations from phylogenetically diverse animals such as cnidarians and chordates. Although few studies directly compared them, minor differences were noted between gap junctions of vertebrates and invertebrates. For instance, a slightly wider gap was noted between cells of invertebrates and the spacing between invertebrate channels was generally greater. Connexins were identified as the structural component of vertebrate junctions in the 1980s and innexins as the structural component of pre-chordate junctions in the 1990s. Despite a lack of similarity in gene sequence, connexins and innexins are remarkably similar. Innexins and connexins have the same membrane topology and form intercellular channels that play a variety of tissue- and temporally specific roles. Both protein types oligomerize to form large aqueous channels that allow the passage of ions and small metabolites and are regulated by factors such as pH, calcium, and voltage. Much more is currently known about the structure, function, and structure-function relationships of connexins. However, the innexin field is expanding. Greater knowledge of innexin channels will permit more detailed comparisons with their connexin-based counterparts, and provide insight into the ubiquitous yet specific roles of gap junctions. © 2016 Wiley Periodicals, Inc. Develop Neurobiol 77: 522-547, 2017.

Gap junction structure: unraveled, but not fully revealed

Gap junction channels facilitate the intercellular exchange of ions and small molecules, a process that is critical for the function of many different kinds of cells and tissues. Recent crystal structures of channels formed by one connexin isoform (connexin26) have been determined, and they have been subjected to molecular modeling. These studies have provided high-resolution models to gain insights into the mechanisms of channel conductance, molecular permeability, and gating. The models share similarities, but there are some differences in the conclusions reached by these studies. Many unanswered questions remain to allow an atomic-level understanding of intercellular communication mediated by connexin26. Because some domains of the connexin polypeptides are highly conserved (like the transmembrane regions), it is likely that some features of the connexin26 structure will apply to other members of the family of gap junction proteins. However, determination of high-resolution structures and modeling of other connexin channels will be required to account for the diverse biophysical properties and regulation conferred by the differences in their sequences.

Accessing gap-junction channel structure-function relationships through molecular modeling and simulations

Background

Gap junction channels (GJCs) are massive protein channels connecting the cytoplasm of adjacent cells. These channels allow intercellular transfer of molecules up to ~1 kDa, including water, ions and other metabolites. Unveiling structure-function relationships coded into the molecular architecture of these channels is necessary to gain insight on their vast biological function including electrical synapse, inflammation, development and tissular homeostasis. From early works, computational methods have been critical to analyze and interpret experimental observations. Upon the availability of crystallographic structures, molecular modeling and simulations have become a valuable tool to assess structure-function relationships in GJCs. Modeling different connexin isoforms, simulating the transport process, and exploring molecular variants, have provided new hypotheses and out-of-the-box approaches to the study of these important channels.

Methods

Here, we review foundational structural studies and recent developments on GJCs using molecular modeling and simulation techniques, highlighting the methods and the cross-talk with experimental evidence.

Results and discussion

By comparing results obtained by molecular modeling and simulations techniques with structural and functional information obtained from both recent literature and structural databases, we provide a critical assesment of structure-function relationships that can be obtained from the junction between theoretical and experimental evidence.

Mechanism of gating by calcium in connexin hemichannels

Aberrant opening of nonjunctional connexin hemichannels at the plasma membrane is associated with many diseases, including ischemia and muscular dystrophy. Proper control of hemichannel opening is essential to maintain cell viability and is achieved by physiological levels of extracellular Ca²⁺, which drastically reduce hemichannel activity. Here we examined the role of conserved charged residues that form electrostatic networks near the extracellular entrance of the connexin pore, a region thought to be involved in gating rearrangements of hemichannels. Molecular dynamics simulations indicate discrete sites for Ca²⁺ interaction and consequent disruption of salt bridges in the open hemichannels. Experimentally, we found that disruption of these salt bridges by mutations facilitates hemichannel closing. Two negatively charged residues in these networks are putative Ca²⁺ binding sites, forming a Ca²⁺-gating ring near the extracellular entrance of the pore. Accessibility studies showed that this Ca²⁺-bound gating ring does not prevent access of ions or small molecules to positions deeper into the pore, indicating that the physical gate is below the Ca²⁺-gating ring. We conclude that intra- and intersubunit electrostatic networks at the extracellular entrance of the hemichannel pore play critical roles in hemichannel gating reactions and are tightly controlled by extracellular Ca²⁺ Our findings provide a general mechanism for Ca²⁺ gating among different connexin hemichannel isoforms.