A Calpha model for the transmembrane alpha helices of gap junction intercellular channels

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.

Structures of wild-type and selected CMT1X mutant connexin 32 gap junction channels and hemichannels

In myelinating Schwann cells, connection between myelin layers is mediated by gap junction channels (GJCs) formed by docked connexin 32 (Cx32) hemichannels (HCs). Mutations in Cx32 cause the X-linked Charcot-Marie-Tooth disease (CMT1X), a degenerative neuropathy without a cure. A molecular link between Cx32 dysfunction and CMT1X pathogenesis is still missing. Here, we describe the high-resolution cryo-electron cryo-myography (cryo-EM) structures of the Cx32 GJC and HC, along with two CMT1X-linked mutants, W3S and R22G. While the structures of wild-type and mutant GJCs are virtually identical, the HCs show a major difference: In the W3S and R22G mutant HCs, the amino-terminal gating helix partially occludes the pore, consistent with a diminished HC activity. Our results suggest that HC dysfunction may be involved in the pathogenesis of CMT1X.

The Multifaceted Role of Connexins in Tumor Microenvironment Initiation and Maintenance

Today's research on the processes of carcinogenesis and the vital activity of tumor tissues implies more attention be paid to constituents of the tumor microenvironment and their interactions. These interactions between cells in the tumor microenvironment can be mediated via different types of protein junctions. Connexins are one of the major contributors to intercellular communication. They form the gap junctions responsible for the transfer of ions, metabolites, peptides, miRNA, etc., between neighboring tumor cells as well as between tumor and stromal cells. Connexin hemichannels mediate purinergic signaling and bidirectional molecular transport with the extracellular environment. Additionally, connexins have been reported to localize in tumor-derived exosomes and facilitate the release of their cargo. A large body of evidence implies that the role of connexins in cancer is multifaceted. The pro- or anti-tumorigenic properties of connexins are determined by their abundance, localization, and functionality as well as their channel assembly and non-channel functions. In this review, we have summarized the data on the contribution of connexins to the formation of the tumor microenvironment and to cancer initiation and progression.

Proposed Role for Internal Lens Pressure as an Initiator of Age-Related Lens Protein Aggregation Diseases

The process that initiates lens stiffness evident in age-related lens protein aggregation diseases is thought to be mainly the result of oxidation. While oxidation is a major contributor, the exposure of lens proteins to physical stress over time increases susceptibility of lens proteins to oxidative damage, and this is believed to play a significant role in initiating these diseases. Accordingly, an overview of key physical stressors and molecular factors known to be implicated in the development of age-related lens protein aggregation diseases is presented, paying particular attention to the consequence of persistent increase in internal lens pressure.

Calmodulin-Connexin Partnership in Gap Junction Channel Regulation-Calmodulin-Cork Gating Model

In the past four decades numerous findings have indicated that gap junction channel gating is mediated by intracellular calcium concentrations ([Ca²⁺_i]) in the high nanomolar range via calmodulin (CaM). We have proposed a CaM-based gating model based on evidence for a direct CaM role in gating. This model is based on the following: CaM inhibitors and the inhibition of CaM expression to prevent chemical gating. A CaM mutant with higher Ca²⁺ sensitivity greatly increases gating sensitivity. CaM co-localizes with connexins. Connexins have high-affinity CaM-binding sites. Connexin mutants paired to wild type connexins have a higher gating sensitivity, which is eliminated by the inhibition of CaM expression. Repeated trans-junctional voltage (Vj) pulses progressively close channels by the chemical/slow gate (CaM’s N-lobe). At the single channel level, the gate closes and opens slowly with on-off fluctuations. Internally perfused crayfish axons lose gating competency but recover it by the addition of Ca-CaM to the internal perfusion solution. X-ray diffraction data demonstrate that isolated gap junctions are gated at the cytoplasmic end by a particle of the size of a CaM lobe. We have proposed two types of CaM-driven gating: “Ca-CaM-Cork” and “CaM-Cork”. In the first, the gating involves Ca²⁺-induced CaM activation. In the second, the gating occurs without a [Ca²⁺]_i rise.

Cryo-EM structure of an open conformation of a gap junction hemichannel in lipid bilayer nanodiscs

To mediate cell-to-cell communication via gap junction channels (GJCs), connexins (Cx) traffic as hexameric hemichannels to the plasma membrane, which dock end-to-end between adjacent cell membranes, thereby forming a dodecameric intercellular conduit. Hemichannels also function independently to mediate the passage of contents between the cytoplasm and extracellular space. To generate hemichannels, the mutation N176Y was introduced into the second extracellular loop of Cx26. The electron cryomicroscopy structure of the hexameric hemichannel in lipid bilayer nanodiscs displays an open pore and a 4-helix bundle transmembrane design that is nearly identical to dodecameric GJCs. In contrast to the high resolution of the transmembrane α-helices, the extracellular loops are less well resolved. The conformational flexibility of the extracellular loops may be essential to facilitate surveillance of hemichannels in apposed cells to identify compatible Cx isoforms that enable intercellular docking. Our results also provide a structural foundation for previous electrophysiologic and permeation studies of Cx hemichannels.

Aging-dependent loss of GAP junction proteins Cx46 and Cx50 in the fiber cells of human and mouse lenses accounts for the diminished coupling conductance

The homeostasis of the ocular lens is maintained by a microcirculation system propagated through gap junction channels. It is well established that the intercellular communications of the lens become deteriorative during aging. However, the molecular basis for this change in human lenses has not been well defined. Here, we present evidence to show that over 90% of Cx46 and Cx50 are lost in the fiber cells of normal human lenses aged 50 and above. From transparent to cataractous lenses, while Cx43 was upregulated, both Cx46 and Cx50 were significantly down-regulated in the lens epithelia. During aging of mouse lenses, Cx43 remained unchanged, but both Cx46 and Cx50 were significantly downregulated. Under oxidative stress treatment, mouse lenses develop in vitro cataractogenesis. Associated with this process, Cx43 was significantly upregulated, in contrast, Cx46 and Cx50 were sharply downregulated. Together, our results for the first time reveal that downregulation in Cx46 and Cx50 levels appears to be the major reason for the diminished coupling conductance, and the aging-dependent loss of Cx46 and Cx50 promotes senile cataractogenesis.

A Steric “Ball-and-Chain” Mechanism for pH-Mediated Regulation of Gap Junction Channels

Gap junction channels (GJCs) mediate intercellular communication and are gated by numerous conditions such as pH. The electron cryomicroscopy (cryo-EM) structure of Cx26 GJC at physiological pH recapitulates previous GJC structures in lipid bilayers. At pH 6.4, we identify two conformational states, one resembling the open physiological-pH structure and a closed conformation that displays six threads of density, that join to form a pore-occluding density. Crosslinking and hydrogen-deuterium exchange mass spectrometry reveal closer association between the N-terminal (NT) domains and the cytoplasmic loops (CL) at acidic pH. Previous electrophysiologic studies suggest an association between NT residue N14 and H100 near M2, which may trigger the observed movement of M2 toward M1 in our cryo-EM maps, thereby accounting for additional NT-CL crosslinks at acidic pH. We propose that these pH-induced interactions and conformational changes result in extension, ordering, and association of the acetylated NT domains to form a hexameric "ball-and-chain" gating particle.

Calmodulin-Mediated Regulation of Gap Junction Channels

Evidence that neighboring cells uncouple from each other as one dies surfaced in the late 19th century, but it took almost a century for scientists to start understanding the uncoupling mechanism (chemical gating). The role of cytosolic free calcium (Ca2+i) in cell-cell channel gating was first reported in the mid-sixties. In these studies, only micromolar [Ca2+]i were believed to affect gating-concentrations reachable only in cell death, which would discard Ca2+i as a fine modulator of cell coupling. More recently, however, numerous researchers, including us, have reported the effectiveness of nanomolar [Ca2+]i. Since connexins do not have high-affinity calcium sites, the effectiveness of nanomolar [Ca2+]i suggests the role of Ca-modulated proteins, with calmodulin (CaM) being most obvious. Indeed, in 1981 we first reported that a CaM-inhibitor prevents chemical gating. Since then, the CaM role in gating has been confirmed by studies that tested it with a variety of approaches such as treatments with CaM-inhibitors, inhibition of CaM expression, expression of CaM mutants, immunofluorescent co-localization of CaM and gap junctions, and binding of CaM to peptides mimicking connexin domains identified as CaM targets. Our gating model envisions Ca2+-CaM to directly gate the channels by acting as a plug ("Cork" gating model), and probably also by affecting connexin conformation.

A lipophilicity-based energy function for membrane-protein modelling and design

Membrane-protein design is an exciting and increasingly successful research area which has led to landmarks including the design of stable and accurate membrane-integral proteins based on coiled-coil motifs. Design of topologically more complex proteins, such as most receptors, channels, and transporters, however, demands an energy function that balances contributions from intra-protein contacts and protein-membrane interactions. Recent advances in water-soluble all-atom energy functions have increased the accuracy in structure-prediction benchmarks. The plasma membrane, however, imposes different physical constraints on protein solvation. To understand these constraints, we recently developed a high-throughput experimental screen, called dsTβL, and inferred apparent insertion energies for each amino acid at dozens of positions across the bacterial plasma membrane. Here, we express these profiles as lipophilicity energy terms in Rosetta and demonstrate that the new energy function outperforms previous ones in modelling and design benchmarks. Rosetta ab initio simulations starting from an extended chain recapitulate two-thirds of the experimentally determined structures of membrane-spanning homo-oligomers with <2.5Å root-mean-square deviation within the top-predicted five models (available online: http://tmhop.weizmann.ac.il). Furthermore, in two sequence-design benchmarks, the energy function improves discrimination of stabilizing point mutations and recapitulates natural membrane-protein sequences of known structure, thereby recommending this new energy function for membrane-protein modelling and design.

Potential of cryo-EM for high-resolution structural analysis of gap junction channels

Gap junction family proteins form conduits connecting the cytoplasm of adjacent cells, thereby enabling electrical and chemical coupling to maintain physiological homeostasis. Gap junction proteins comprise two gene families, connexins in chordates and innexins in pre-chordates. Their channel structures have been analyzed by electron or X-ray crystallography, but only a few atomic structures have been reported. Recent advances in single-particle cryo-electron microscopy (cryo-EM) will help to elucidate these structures further. Here the structural biology of gap junction channels utilizing crystallography and single-particle cryo-EM is overviewed to shed light on the functional mechanisms of cell-cell communication that are essential for multicellular organisms.

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.

Automatic Inference of Sequence from Low-Resolution Crystallographic Data

At resolutions worse than 3.5 Å, the electron density is weak or nonexistent at the locations of the side chains. Consequently, the assignment of the protein sequences to their correct positions along the backbone is a difficult problem. In this work, we propose a fully automated computational approach to assign sequence at low resolution. It is based on our surprising observation that standard reciprocal-space indicators, such as the initial unrefined R value, are sensitive enough to detect an erroneous sequence assignment of even a single backbone position. Our approach correctly determines the amino acid type for 15%, 13%, and 9% of the backbone positions in crystallographic datasets with resolutions of 4.0 Å, 4.5 Å, and 5.0 Å, respectively. We implement these findings in an application for threading a sequence onto a backbone structure. For the three resolution ranges, the application threads 83%, 81%, and 64% of the sequences exactly as in the deposited PDB structures.

Cx32 hemichannel opening by cytosolic Ca2+ is inhibited by the R220X mutation that causes Charcot-Marie-Tooth disease

Mutations of the GJB1 gene encoding connexin 32 (Cx32) cause the X-linked form of Charcot-Marie-Tooth disease (CMTX1), a demyelinating peripheral neuropathy for which there is no cure. A growing body of evidence indicates that ATP release through Cx32 hemichannels in Schwann cells could be critical for nerve myelination, but it is unknown if CMTX1 mutations alter the cytosolic Ca2+-dependent gating mechanism that controls Cx32 hemichannel opening and ATP release. The current study uncovered that loss of the C-terminus in Cx32 (R220X mutation), which causes a severe CMTX1 phenotype, inhibits hemichannel opening during a canonical IP3-mediated increase in cytosolic Ca2+ in HeLa cells. Interestingly, the gating function of R220X hemichannels was completely restored by both the intracellular and extracellular application of a peptide that mimics the Cx32 cytoplasmic loop. All-atom molecular dynamics simulations suggest that loss of the C-terminus in the mutant hemichannel triggers abnormal fluctuations of the cytoplasmic loop which are prevented by binding to the mimetic peptide. Experiments that stimulated R220X hemichannel opening by cell depolarization displayed reduced voltage sensitivity with respect to wild-type hemichannels which was explained by loss of subconductance states at the single channel level. Finally, experiments of intercellular diffusion mediated by wild-type or R220X gap junction channels revealed similar unitary permeabilities to ions, signalling molecules (cAMP) or larger solutes (Lucifer yellow). Taken together, our findings support the hypothesis that paracrine signalling alteration due to Cx32 hemichannel dysfunction underlies CMTX1 pathogenesis and suggest a candidate molecule for novel studies investigating a therapeutic approach.

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.

Connexinopathies: a structural and functional glimpse

Mutations in human connexin (Cx) genes have been related to diseases, which we termed connexinopathies. Such hereditary disorders include nonsyndromic or syndromic deafness (Cx26, Cx30), Charcot Marie Tooth disease (Cx32), occulodentodigital dysplasia and cardiopathies (Cx43), and cataracts (Cx46, Cx50). Despite the clinical phenotypes of connexinopathies have been well documented, their pathogenic molecular determinants remain elusive. The purpose of this work is to identify common/uncommon patterns in channels function among Cx mutations linked to human diseases. To this end, we compiled and discussed the effect of mutations associated to Cx26, Cx32, Cx43, and Cx50 over gap junction channels and hemichannels, highlighting the function of the structural channel domains in which mutations are located and their possible role affecting oligomerization, gating and perm/selectivity processes.

An electrostatic mechanism for Ca(2+)-mediated regulation of gap junction channels

Gap junction channels mediate intercellular signalling that is crucial in tissue development, homeostasis and pathologic states such as cardiac arrhythmias, cancer and trauma. To explore the mechanism by which Ca(2+) blocks intercellular communication during tissue injury, we determined the X-ray crystal structures of the human Cx26 gap junction channel with and without bound Ca(2+). The two structures were nearly identical, ruling out both a large-scale structural change and a local steric constriction of the pore. Ca(2+) coordination sites reside at the interfaces between adjacent subunits, near the entrance to the extracellular gap, where local, side chain conformational rearrangements enable Ca(2+)chelation. Computational analysis revealed that Ca(2+)-binding generates a positive electrostatic barrier that substantially inhibits permeation of cations such as K(+) into the pore. Our results provide structural evidence for a unique mechanism of channel regulation: ionic conduction block via an electrostatic barrier rather than steric occlusion of the channel pore.

Differential regulation of Connexin50 and Connexin46 by PI3K signaling

Gap junction channels can modify their activity in response to cell signaling pathways. Here, we demonstrate that Connexin50 (Cx50) coupling, but not Connexin46 (Cx46), increased when co-expressed with a constitutively active p110α subunit of PI3K in Xenopus oocytes. In addition, inhibition of PI3K signaling by blocking p110α, or Akt, significantly decreased gap junctional conductance in Cx50 transfected HeLa cells, with no effect on Cx46. Alterations in coupling levels were not a result of reduced Cx50 unitary conductance, suggesting that changes in the number of active channels were responsible. These data indicate that Cx50 is specifically regulated by the PI3K signaling pathway.

Gap junctional regulation of pressure, fluid force, and electrical fields in the epigenetics of cardiac morphogenesis and remodeling

Epigenetic factors of pressure load, fluid force, and electrical fields that occur during cardiac contraction affect cardiac development, morphology, function, and pathogenesis. These factors are orchestrated by intercellular communication mediated by gap junctions, which synchronize action potentials and second messengers. Misregulation of the gap junction protein connexin (Cx) alters cardiogenesis, and can be a pathogenic factor causing cardiac conduction disturbance, fatal arrhythmia, and cardiac remodeling in disease states such as hypertension and ischemia. Changes in Cx expression can occur even when the DNA sequence of the Cx gene itself is unaltered. Posttranslational modifications might reduce arrhythmogenic substrates, improve cardiac function, and promote remodeling in a diseased heart. In this review, we discuss the epigenetic features of gap junctions that regulate cardiac morphology and remodeling. We further discuss potential clinical applications of current knowledge of the structure and function of gap junctions.

Structure and closure of connexin gap junction channels

Connexin gap junctions comprise assembled channels penetrating two plasma membranes for which gating regulation is associated with a variety of factors, including voltage, pH, Ca(2+), and phosphorylation. Functional studies have established that various parts of the connexin peptides are related to channel closure and electrophysiology studies have provided several working models for channel gating. The corresponding structural models supporting these findings, however, are not sufficient because only small numbers of closed connexin structures have been reported. To fully understand the gating mechanisms, the channels should be visualized in both the open and closed states. Electron crystallography and X-ray crystallography studies recently revealed three-dimensional structures of connexin channels in a couple of states in which the main difference is the conformation of the N-terminal domain, which have helped to clarify the structure in regard to channel closure. Here the closure models for connexin gap junction channels inferred from structural and functional studies are described in the context of each domain of the connexin protein associated with gating modulation.

A history of gap junction structure: hexagonal arrays to atomic resolution

Gap junctions are specialized membrane structures that provide an intercellular pathway for the propagation and/or amplification of signaling cascades responsible for impulse propagation, cell growth, and development. Prior to the identification of the proteins that comprise gap junctions, elucidation of channel structure began with initial observations of a hexagonal nexus connecting apposed cellular membranes. Concomitant with technological advancements spanning over 50 years, atomic resolution structures are now available detailing channel architecture and the cytoplasmic domains that have helped to define mechanisms governing the regulation of gap junctions. Highlighted in this review are the seminal structural studies that have led to our current understanding of gap junction biology.

The 3.5 ångström X-ray structure of the human connexin26 gap junction channel is unlikely that of a fully open channel

The permeability of gap junction channels to metabolites, and not simply to small inorganic ions, is likely to play an important role in development, physiology as well as in etiology of several diseases. Here, we combined dual patch clamp and fluorescence imaging techniques with molecular dynamics (MD) simulations to investigate the permeation of calcein, a relatively large fluorescent tracer (MW 622 Da) through homomeric gap junction channels formed by wild type human connexin26 (hCx26wt) protomers. Our experimental data indicate that the unitary flux of calcein driven by a 125 μM concentration difference is J_pore = 226 molecule/s per channel. In the light of Eyring transition state theory adapted for the liquid phase, this value corresponds to an energy barrier of ~20 k_BT (where k_B is the Boltzmann constant and T is absolute temperature). The barrier predicted by our MD simulations, based on the 3.5 Å X–ray structural model of the hCx26wt gap junction channel, is ~45 k_BT. The main contributions to the energetics of calcein permeation originated from the interaction between the permeating molecule and the charged aminoacids lining the channel pore. Assigning a fake zero total charge to the calcein molecule yielded a value for the barrier height compatible with the experimental data. These results can be accounted for by two different (although not mutually exclusive) hypotheses: (1) the X–ray model of the hCx26wt gap junction channel is not representative of a fully open state; (2) post translational modifications affecting the hCx26wt protein in our expression system differed from the modifications undergone by the proteins in the conditions used to obtain the crystal structure. Hypothesis (1) is compatible with data indicating that, only 10% or less of the channels forming a gap junction plaque are in the open state, and therefore the averaging procedure intrinsic in the generation of the crystal structure data more closely reflects that of a closed channel. Hypothesis (2) is compatible with recent mass spectrometry data and implies that the charge of several amino acid side chains may have been altered, thus modifying substantially the permeation properties of the channels in living cells.

Permeation pathway of homomeric connexin 26 and connexin 30 channels investigated by molecular dynamics

Mutations in the genes GJB2 and GJB6 encoding human connnexin26 (hCx26) and connexin30 (hCx30), respectively, are the leading cause of non-syndromic prelingual deafness in several human populations. In this work, we exploited the high degree (77%) of sequence similarity shared by hCx26 and hCx30 to create atomistic models of homomeric hCx26 and hCx30 connexons starting from the X-ray crystallographic structure of an intercellular channel formed by hCx26 protomers at 3.5-å resolution. The equilibrium dynamics of the two protein complexes was followed for 40 ns each by Molecular Dynamics (MD) simulations. Our results indicate that, in hCx26, positively charged Lys41 residues establish a potential barrier within the fully open channel, hindering ion diffusion in the absence of an electrochemical gradient. A similar role is played, in hCx30, by negatively charged Glu49 residues. The different position and charge of these two ion sieves account for the differences in unitary conductance observed experimentally. Our results are discussed in terms of present models of voltage gating in connexin channels.

Molecular dynamics simulations of the Cx26 hemichannel: evaluation of structural models with Brownian dynamics

The recently published crystal structure of the Cx26 gap junction channel provides a unique opportunity for elucidation of the structure of the conductive connexin pore and the molecular determinants of its ion permeation properties (conductance, current-voltage [I-V] relations, and charge selectivity). However, the crystal structure was incomplete, most notably lacking the coordinates of the N-terminal methionine residue, which resides within the pore, and also lacking two cytosolic domains. To allow computational studies for comparison with the known channel properties, we completed the structure. Grand canonical Monte Carlo Brownian dynamics (GCMC/BD) simulations of the completed and the published Cx26 hemichannel crystal structure indicate that the pore is too narrow to permit significant ion flux. The GCMC/BD simulations predict marked inward current rectification and almost perfect anion selectivity, both inconsistent with known channel properties. The completed structure was refined by all-atom molecular dynamics (MD) simulations (220 ns total) in an explicit solvent and POPC membrane system. These MD simulations produced an equilibrated structure with a larger minimal pore diameter, which decreased the height of the permeation barrier formed by the N terminus. GCMC/BD simulations of the MD-equilibrated structure yielded more appropriate single-channel conductance and less anion/cation selectivity. However, the simulations much more closely matched experimentally determined I-V relations when the charge effects of specific co- and posttranslational modifications of Cx26 previously identified by mass spectrometry were incorporated. We conclude that the average equilibrated structure obtained after MD simulations more closely represents the open Cx26 hemichannel structure than does the crystal structure, and that co- and posttranslational modifications of Cx26 hemichannels are likely to play an important physiological role by defining the conductance and ion selectivity of Cx26 channels. Furthermore, the simulations and data suggest that experimentally observed heterogeneity in Cx26 I-V relations can be accounted for by variation in co- and posttranslational modifications.

Computational studies of membrane proteins: models and predictions for biological understanding

We discuss recent progresses in computational studies of membrane proteins based on physical models with parameters derived from bioinformatics analysis. We describe computational identification of membrane proteins and prediction of their topology from sequence, discovery of sequence and spatial motifs, and implications of these discoveries. The detection of evolutionary signal for understanding the substitution pattern of residues in the TM segments and for sequence alignment is also discussed. We further discuss empirical potential functions for energetics of inserting residues in the TM domain, for interactions between TM helices or strands, and their applications in predicting lipid-facing surfaces of the TM domain. Recent progresses in structure predictions of membrane proteins are also reviewed, with further discussions on calculation of ensemble properties such as melting temperature based on simplified state space model. Additional topics include prediction of oligomerization state of membrane proteins, identification of the interfaces for protein-protein interactions, and design of membrane proteins. This article is part of a Special Issue entitled: Protein Folding in Membranes.

Structure of the gap junction channel and its implications for its biological functions

Gap junctions consist of arrays of intercellular channels composed of integral membrane proteins called connexin in vertebrates. Gap junction channels regulate the passage of ions and biological molecules between adjacent cells and, therefore, are critically important in many biological activities, including development, differentiation, neural activity, and immune response. Mutations in connexin genes are associated with several human diseases, such as neurodegenerative disease, skin disease, deafness, and developmental abnormalities. The activity of gap junction channels is regulated by the membrane voltage, intracellular microenvironment, interaction with other proteins, and phosphorylation. Each connexin channel has its own property for conductance and molecular permeability. A number of studies have tried to reveal the molecular architecture of the channel pore that should confer the connexin-specific permeability/selectivity properties and molecular basis for the gating and regulation. In this review, we give an overview of structural studies and describe the structural and functional relationship of gap junction channels.

Asymmetric configurations and N-terminal rearrangements in connexin26 gap junction channels

Gap junction channels are unique in that they possess multiple mechanisms for channel closure, several of which involve the N terminus as a key component in gating, and possibly assembly. Here, we present electron crystallographic structures of a mutant human connexin26 (Cx26M34A) and an N-terminal deletion of this mutant (Cx26M34Adel2-7) at 6-Å and 10-Å resolutions, respectively. The three-dimensional map of Cx26M34A was improved by data from 60° tilt images and revealed a breakdown of the hexagonal symmetry in a connexin hemichannel, particularly in the cytoplasmic domain regions at the ends of the transmembrane helices. The Cx26M34A structure contained an asymmetric density in the channel vestibule ("plug") that was decreased in the Cx26M34Adel2-7 structure, indicating that the N terminus significantly contributes to form this plug feature. Functional analysis of the Cx26M34A channels revealed that these channels are predominantly closed, with the residual electrical conductance showing normal voltage gating. N-terminal deletion mutants with and without the M34A mutation showed no electrical activity in paired Xenopus oocytes and significantly decreased dye permeability in HeLa cells. Comparing this closed structure with the recently published X-ray structure of wild-type Cx26, which is proposed to be in an open state, revealed a radial outward shift in the transmembrane helices in the closed state, presumably to accommodate the N-terminal plug occluding the pore. Because both Cx26del2-7 and Cx26M34Adel2-7 channels are closed, the N terminus appears to have a prominent role in stabilizing the open configuration.

Design, synthesis, and properties of branch-chained maltoside detergents for stabilization and crystallization of integral membrane proteins: human connexin 26

A challenging requirement for structural studies of integral membrane proteins (IMPs) is the use of amphiphiles that replicate the hydrophobic environment of membranes. Progress has been impeded by the limited number of useful detergents and the need for a deeper understanding of their structure-activity relationships. To this end, we designed a family of detergents containing short, branched alkyl chains at the interface between the polar head and the apolar tail. This design mimics the second aliphatic chain of lipid molecules and reduces water penetration, thereby increasing the hydrophobicity within the interior of the micelle. To compare with the popular straight-chained maltoside detergents, the branch-chained beta-D-maltosides were synthesized efficiently in pure anomeric form. The branch-chained maltosides form smaller micelles by having shorter main chains, while having comparable hydrophobicity to the detergents with only straight chains. Selected branch-chained and straight-chained maltoside detergents were examined for their ability to solubilize, stabilize, and aid the crystallization of human connexin 26, an alpha-helical IMP that forms hexamers. We showed that the branch-chained maltosides with optimized micellar properties performed as well as or better than the straight-chained analogues and enabled crystallization in different space groups.

Quality assessment of protein model-structures using evolutionary conservation

Motivation: Programs that evaluate the quality of a protein structural model are important both for validating the structure determination procedure and for guiding the model-building process. Such programs are based on properties of native structures that are generally not expected for faulty models. One such property, which is rarely used for automatic structure quality assessment, is the tendency for conserved residues to be located at the structural core and for variable residues to be located at the surface.

Results: We present ConQuass, a novel quality assessment program based on the consistency between the model structure and the protein's conservation pattern. We show that it can identify problematic structural models, and that the scores it assigns to the server models in CASP8 correlate with the similarity of the models to the native structure. We also show that when the conservation information is reliable, the method's performance is comparable and complementary to that of the other single-structure quality assessment methods that participated in CASP8 and that do not use additional structural information from homologs.

Availability: A perl implementation of the method, as well as the various perl and R scripts used for the analysis are available at http://bental.tau.ac.il/ConQuass/.

Contact:nirb@tauex.tau.ac.il

Supplementary information:Supplementary data are available at Bioinformatics online.

Structures of membrane proteins

In reviewing the structures of membrane proteins determined up to the end of 2009, we present in words and pictures the most informative examples from each family. We group the structures together according to their function and architecture to provide an overview of the major principles and variations on the most common themes. The first structures, determined 20 years ago, were those of naturally abundant proteins with limited conformational variability, and each membrane protein structure determined was a major landmark. With the advent of complete genome sequences and efficient expression systems, there has been an explosion in the rate of membrane protein structure determination, with many classes represented. New structures are published every month and more than 150 unique membrane protein structures have been determined. This review analyses the reasons for this success, discusses the challenges that still lie ahead, and presents a concise summary of the key achievements with illustrated examples selected from each class.

Complex formation between calmodulin and a peptide from the intracellular loop of the gap junction protein connexin43: Molecular conformation and energetics of binding

Gap junctions are formed by a family of transmembrane proteins, connexins. Connexin43 is a widely studied member of the family, being ubiquitously expressed in a variety of tissues and a target of a large number of disease mutations. The intracellular loop of connexin43 has been shown to include a calmodulin binding domain, but detailed 3-dimensional data on the structure of the complex are not available. In this study, we used a synthetic peptide from this domain to reveal the conformation of the calmodulin-peptide complex by small angle X-ray scattering. Upon peptide binding, calmodulin lost its dumbbell shape, adopting a more globular conformation. We also studied the energetics of the interaction using calorimetry and computational methods. All our data indicate that calmodulin binds to the peptide from cx43 in the classical 'collapsed' conformation.

A description of the structural determination procedures of a gap junction channel at 3.5 A resolution

Intercellular signalling is an essential characteristic of multicellular organisms. Gap junctions, which consist of arrays of intercellular channels, permit the exchange of ions and small molecules between adjacent cells. Here, the structural determination of a gap junction channel composed of connexin 26 (Cx26) at 3.5 A resolution is described. During each step of the purification process, the protein was examined using electron microscopy and/or dynamic light scattering. Dehydration of the crystals improved the resolution limits. Phase refinement using multi-crystal averaging in conjunction with noncrystallographic symmetry averaging based on strictly determined noncrystallographic symmetry operators resulted in an electron-density map for model building. The amino-acid sequence of a protomer structure consisting of the amino-terminal helix, four transmembrane helices and two extracellular loops was assigned to the electron-density map. The amino-acid assignment was confirmed using six selenomethionine (SeMet) sites in the difference Fourier map of the SeMet derivative and three intramolecular disulfide bonds in the anomalous difference Fourier map of the native crystal.

Three residues at the interface of factor XI (FXI) monomers augment covalent dimerization of FXI

BACKGROUND

Human plasma factor XI is a homodimer, with each monomer comprising a catalytic domain and four homologous 'apple' domains. The monomers bind to each other through non-covalent bonds and through a disulfide bond between Cys321 residues in apple 4 domains.

OBJECTIVE

To identify residues essential for dimerization in the FXI monomer interface.

METHODS

Specificity-determining residues in apple 4 domains were sought by sequence alignment of FXI and prekallikrein apple domains in different species. Specific residues identified in apple 4 domains were mutagenized and expressed in baby hamster kidney (BHK) cells for evaluation of their effect on FXI dimerization, analyzed by non-reduced sodium dodecylsulfate polyacrylamide gel electrophoresis and size-exclusion chromatography.

RESULTS

Among the 19 residues of the FXI monomer interface, Leu284, Ile290 and Tyr329 were defined as specificity-determining residues. Substitutions of these residues or pairs of residues did not affect FXI synthesis and secretion from transfected BHK cells, but did impair dimerization, despite the presence of cysteine at position 321. The double mutant 284A/290A yielded predominantly a monomer, whereas all other single or double mutants yielded monomers as well as disulfide-bonded dimers.

CONCLUSIONS

The data suggest that Leu284, Ile290 and Tyr329 in the interface of FXI monomers are essential for forming non-covalently bonded dimers that facilitate formation of a disulfide-bonded stable FXI dimer.

Conformational changes in a pore-forming region underlie voltage-dependent "loop gating" of an unapposed connexin hemichannel

The structure of the pore is critical to understanding the molecular mechanisms underlying selective permeation and voltage-dependent gating of channels formed by the connexin gene family. Here, we describe a portion of the pore structure of unapposed hemichannels formed by a Cx32 chimera, Cx32*Cx43E1, in which the first extracellular loop (E1) of Cx32 is replaced with the E1 of Cx43. Cysteine substitutions of two residues, V38 and G45, located in the vicinity of the border of the first transmembrane (TM) domain (TM1) and E1 are shown to react with the thiol modification reagent, MTSEA-biotin-X, when the channel resides in the open state. Cysteine substitutions of flanking residues A40 and A43 do not react with MTSEA-biotin-X when the channel resides in the open state, but they react with dibromobimane when the unapposed hemichannels are closed by the voltage-dependent "loop-gating" mechanism. Cysteine substitutions of residues V37 and A39 do not appear to be modified in either state. Furthermore, we demonstrate that A43C channels form a high affinity Cd2+ site that locks the channel in the loop-gated closed state. Biochemical assays demonstrate that A43C can also form disulfide bonds when oocytes are cultured under conditions that favor channel closure. A40C channels are also sensitive to micromolar Cd2+ concentrations when closed by loop gating, but with substantially lower affinity than A43C. We propose that the voltage-dependent loop-gating mechanism for Cx32*Cx43E1 unapposed hemichannels involves a conformational change in the TM1/E1 region that involves a rotation of TM1 and an inward tilt of either each of the six connexin subunits or TM1 domains.

Conformational maturation and post-ER multisubunit assembly of gap junction proteins

For all previously well-characterized oligomeric integral membrane proteins, folding, multisubunit assembly, and recognition of conformationally immature molecules for degradation occurs at their organelle of synthesis. This cannot, however, be the case for the gap junction-forming protein connexin43 (Cx43), which when endogenously expressed undergoes multisubunit assembly into connexons only after its transport to the trans-Golgi network. We have developed two novel assays to assess Cx43 folding and assembly: acquisition of resistance of disulfide bonds to reduction by extracellularly added DTT and Triton X-114 detergent phase partitioning. We show that Cx43 synthesized at physiologically relevant levels undergoes a multistep conformational maturation process in which folding of connexin monomers within the ER is a prerequisite for multisubunit assembly in the TGN. Similar results were obtained with Cx32, disproving the widely reported contention that the site of endogenous beta connexin assembly is the ER. Exogenous overexpression of Cx43, Cx32, or Cx26 allows these events to take place within the ER, the first example of the TGN and ER as alternative sites for oligomeric assembly. Our findings also constitute the first biochemical evidence that defective connexin folding is a cause of the human disorder X-linked Charcot-Marie-Tooth disease.

ERp29 restricts Connexin43 oligomerization in the endoplasmic reticulum

Connexin43 (Cx43) is a gap junction protein that forms multimeric channels that enable intercellular communication through the direct transfer of signals and metabolites. Although most multimeric protein complexes form in the endoplasmic reticulum (ER), Cx43 seems to exit from the ER as monomers and subsequently oligomerizes in the Golgi complex. This suggests that one or more protein chaperones inhibit premature Cx43 oligomerization in the ER. Here, we provide evidence that an ER-localized, 29-kDa thioredoxin-family protein (ERp29) regulates Cx43 trafficking and function. Interfering with ERp29 function destabilized monomeric Cx43 oligomerization in the ER, caused increased Cx43 accumulation in the Golgi apparatus, reduced transport of Cx43 to the plasma membrane, and inhibited gap junctional communication. ERp29 also formed a specific complex with monomeric Cx43. Together, this supports a new role for ERp29 as a chaperone that helps stabilize monomeric Cx43 to enable oligomerization to occur in the Golgi apparatus.

Biological and biophysical properties of vascular connexin channels

Intercellular channels formed by connexin proteins play a pivotal role in the direct movement of ions and larger cytoplasmic solutes between vascular endothelial cells, between vascular smooth muscle cells, and between endothelial and smooth muscle cells. Multiple genetic and epigenetic factors modulate connexin expression levels and/or channel function, including cell-type-independent and cell-type-specific transcription factors, posttranslational modifications, and localized membrane targeting. Additionally, differences in protein-protein interactions, including those between connexins, significantly contribute to both vascular homeostasis and disease progression. The biophysical properties of the connexin channels identified in the vasculature, those formed by Cx37, Cx40, Cx43 and/or Cx45 proteins, are discussed in this chapter in the physiological and pathophysiological context of vessel function.

Site-directed mutagenesis reveals putative regions of protein interaction within the transmembrane domains of connexins

Through cysteine-scanning mutagenesis, the authors have compared sites within the transmembrane domains of two connexins, one from the alpha-class (Cx50) and one from the beta-class (Cx32), where amino acid substitution disrupts the function of gap junction channels. In Cx32, 11 sites resulted in no channel function, or an aberrant voltage gating phenotype referred to as "reverse gating," whereas in Cx50, 7 such sites were identified. In both connexins, the sites lie along specific faces of transmembrane helices, suggesting that these may be sites of transmembrane domain interactions. In Cx32, one broad face of the M1 transmembrane domain and a narrower, polar face of M3 were identified, including one site that was shown to come into close apposition with M4 in the closed state. In Cx50, the same face of M3 was identified, but sensitive sites in M1 differed from Cx32. Many fewer sites in M1 disrupted channel function in Cx50, and those that did were on a different helical face to the sensitive sites in Cx32. A more in depth study of two sites in M1 and M2 of Cx32 showed that side-chain length or branching are important for maintenance of normal channel behavior, consistent with this being a site of transmembrane domain interaction.

Projection structure of a N-terminal deletion mutant of connexin 26 channel with decreased central pore density

Gated gap junction channels are important cellular conduits for establishing and maintaining intercellular communication. The three-dimensional structure of a mutant human connexin 26 (Cx26M34A) by electron cryocrystallography revealed a plug-like density in the channel pore suggesting that physical blockage of the pore may be one mechanism of closure (Oshima et al. 2007, Proc Natl Acad Sci USA 104: 10034-10039). However, it remains to be determined what part of the sequence contributes to the plug. Here, we present the projection structure of an N-terminus deletion of Cx26M34A missing amino acids 2 to 7 (Cx26M34Adel2-7) crystallized in the same two-dimensional crystal form. A 10 A resolution projection map of Cx26M34Adel2-7 revealed that the plug density was dramatically reduced in comparison with that found in full-length Cx26 channel. The difference map between the deletion and full-length Cx26M34A channels strongly suggests that the N-terminus of connexin contributes to the plug for the physical closure of gap junction channels.

A fully atomistic model of the Cx32 connexon

Connexins are plasma membrane proteins that associate in hexameric complexes to form channels named connexons. Two connexons in neighboring cells may dock to form a "gap junction" channel, i.e. an intercellular conduit that permits the direct exchange of solutes between the cytoplasm of adjacent cells and thus mediate cell-cell ion and metabolic signaling. The lack of high resolution data for connexon structures has hampered so far the study of the structure-function relationships that link molecular effects of disease-causing mutations with their observed phenotypes. Here we present a combination of modeling techniques and molecular dynamics (MD) to infer side chain positions starting from low resolution structures containing only C alpha atoms. We validated this procedure on the structure of the KcsA potassium channel, which is solved at atomic resolution. We then produced a fully atomistic model of a homotypic Cx32 connexon starting from a published model of the C alpha carbons arrangement for the connexin transmembrane helices, to which we added extracellular and cytoplasmic loops. To achieve structural relaxation within a realistic environment, we used MD simulations inserted in an explicit solvent-membrane context and we subsequently checked predictions of putative side chain positions and interactions in the Cx32 connexon against a vast body of experimental reports. Our results provide new mechanistic insights into the effects of numerous spontaneous mutations and their implication in connexin-related pathologies. This model constitutes a step forward towards a structurally detailed description of the gap junction architecture and provides a structural platform to plan new biochemical and biophysical experiments aimed at elucidating the structure of connexin channels and hemichannels.

Common occurrence of internal repeat symmetry in membrane proteins

Symmetry plays significant roles in protein structure and function. Particularly, symmetric interfaces are known to act as switches for two-state conformational change. Membrane proteins often undergo two-state conformational change during the transport process of ion channels or the active/inactive transitions in receptors. Here, we provide the first comprehensive analyses of internal repeat symmetry in membrane proteins. We examined the known membrane protein structures and found that, remarkably, nearly half of them have internal repeat symmetry. Moreover, we found that the conserved cores of these internal repeats are positioned at the interface of symmetric units when they are mapped on structures. Because of the large sequence divergence that occurs between internal repeats, the inherent symmetry present in protein sequences often has only been detected after structure determination. We therefore developed a sensitive procedure to predict the internal repeat symmetry from sequence information and identified 4653 proteins that are likely to have internal repeat symmetry.

Gap junction communication influences intercellular protein distribution in the lens

Lens transparency and high refractive index presumably depend on the appropriate arrangement and distribution of lens proteins among lens fiber cells. Intercellular gap junction channels formed by alpha3 and alpha8 connexins are known to transport small molecules, ions and water, but not proteins, in the lens. Mosaic expression of green fluorescent protein (GFP) in the lens is a useful marker for monitoring macromolecule distribution between fiber cells and for constructing three-dimensional images of living lens cells. In alpha3(-/-) alpha8(-/-) double knockout (DKO) lenses, three-dimensional images of GFP-positive cells demonstrate the changes of epithelial cell surfaces and insufficient elongation of inner fiber cells. Uniform distribution of GFP between inner lens fiber cells is observed in both wild-type and alpha3(-/-) lenses. In contrast, uniform GFP distribution is slightly delayed in alpha8(-/-) lenses and is abolished in DKO lenses. Without endogenous wild-type alpha3 and alpha8 connexins, knock-in alpha3 connexin (expressed under the alpha8 gene promoter) restores the uniform distribution of GFP protein in the lens. Thus, the presence of either alpha3 or alpha8 connexins seems sufficient to support the uniform distribution of GFP between differentiated lens fiber cells. Although the mechanism that drives GFP transport between fiber cells remains unknown, this work reveals that gap junction communication plays a novel role in the regulation of intercellular protein distribution in the lens.

Purification and reconstitution of the connexin43 carboxyl terminus attached to the 4th transmembrane domain in detergent micelles

In recent years, reports have identified that many eukaryotic proteins contain disordered regions spanning greater than 30 consecutive residues in length. In particular, a number of these intrinsically disordered regions occur in the cytoplasmic segments of plasma membrane proteins. These intrinsically disordered regions play important roles in cell signaling events, as they are sites for protein-protein interactions and phosphorylation. Unfortunately, in many crystallographic studies of membrane proteins, these domains are removed because they hinder the crystallization process. Therefore, a purification procedure was developed to enable the biophysical and structural characterization of these intrinsically disordered regions while still associated with the lipid environment. The carboxyl terminal domain from the gap junction protein connexin43 attached to the 4th transmembrane domain (TM4-Cx43CT) was used as a model system (residues G178-I382). The purification was optimized for structural analysis by nuclear magnetic resonance (NMR) because this method is well suited for small membrane proteins and proteins that lack a well-structured three-dimensional fold. The TM4-Cx43CT was purified to homogeneity with a yield of approximately 6 mg/L from C41(DE3) bacterial cells, reconstituted in the anionic detergent 1-palmitoyl-2-hydroxy-sn-glycero-3-[phospho-RAC-(1-glycerol)], and analyzed by circular dichroism and NMR to demonstrate that the TM4-Cx43CT was properly folded into a functional conformation by its ability to form alpha-helical structure and associate with a known binding partner, the c-Src SH3 domain, respectively.