Kwon T et al. (2011), Molecular dynamics simulations of the Cx26 hemi...

XB-ART-45027

J Gen Physiol 2011 Nov 01;1385:475-93. doi: 10.1085/jgp.201110679.

Show Gene links Show Anatomy links

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

Kwon T , Harris AL , Rossi A , Bargiello TA .

???displayArticle.abstract???
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.

???displayArticle.pubmedLink??? 22006989
???displayArticle.pmcLink??? PMC3206306
???displayArticle.link??? J Gen Physiol
???displayArticle.grants??? [+]

Species referenced: Xenopus
Genes referenced: dtl gjb2 gla gnl3 mgp tbx2

???attribute.lit??? ???displayArticles.show???

	Figure 1. Structural models and pore dimensions of Cx26 hemichannels. (A) Top view of the Cx26 hemichannel corresponding to the crystal structure (Protein Data Bank accession no. 2ZW3; Maeda et al., 2009). The six protein subunits, depicted in different colors, surround a central aqueous pore. (B) Side view of the structure shown in A, illustrating the TM architecture and channel pore relative to two opposite subunits. The extracellular entrance of the channel (z = 45 Å) is at the top of the structure. (C) Pore radius of the Cx26 “crystal” hemichannel determined by HOLE plotted as a function of the z coordinate. (D) Top view of the “completed” crystal structure including the atomic coordinates of the N-terminal methionine residue (Met1), the side chains of K15, S17, and S19, as well as residues comprising the CL and the C terminus. (E) Side view of the completed crystal structure, showing the same subunits as in B. (F) Pore radius of the completed crystal structure. (G) Top view of the channel conformation that has the smallest RMSD from the average equilibrated structure, average pore dimension, and pore-lining probability. We define this structure as the “average equilibrated structure.” (H) Side view of the average equilibrated structure. (I) Pore dimensions of the average equilibrated structure. The mean pore radius is depicted by the black line. The red lines represent the dynamic variation by plotting the mean ± the standard deviation determined from the four 20-ns production-phase simulations.
	Figure 2. Skeleton drawings of selected portions of the crystal and average of all equilibrated structures to show gross changes resulting from structural equilibration by MD simulations. Superimposition of the skeleton of the crystal structure (blue) and the average of all equilibrated structures (red).
	Figure 3. Heterogeneity of I-V relations observed for human Cx26 undocked hemichannels in Xenopus oocytes obtained with outside-out patch-clamp recordings in 100 mM of symmetric KCl. Five I-V relations ranging from slight inward to slight outward are illustrated. (A) Inward/sigmoidal I-V relation obtained with a ±150-mV voltage ramp. 2 of 10 sequential current traces from one record are superimposed. Current traces from −110 to 100 mV show the open-channel I-V relation. (B) Slightly inward I-V relation obtained with a ±70-mV voltage ramp. 3 of 10 sequential traces are superimposed. (C) Linear I-V relation obtained with a ±150-mV voltage ramp. 2 of 10 sequential current traces are superimposed. Current traces are shown from −125 to 100 mV. (D) Slightly outward I-V relation obtained with a ±70-mV voltage ramp. 2 of 10 sequential current traces are superimposed. (E) Outward I-V relation obtained with a ±70-mV voltage ramp. 2 of 10 sequential current traces are superimposed. In all cases, the single-channel records display both forms of voltage-dependent gating reported for Cx26 undocked hemichannels (González et al., 2006; Sánchez et al., 2010). Loop- or slow-gating transitions that can result in full channel closure were observed at large inside negative potentials, and Vj- or fast-gating transitions that result in closure to subconductance states were observed at inside positive potentials. These events correspond to the upward current deflections from the fully open state at negative potentials and downward deflections from the fully open state at positive potentials, respectively, in the records shown. Full channel closures at negative potentials were used to subtract leak conductance in single-channel records.
	Figure 4. (A) Schematic representation of the GCMC/BD system for Cx26 hemichannels. A Cx26 hemichannel (yellow) inserted into the explicit POPC membrane system was integrated into the implicit membrane system of the GCMC/BD program to ensure a tight seal between the irregular conical surface of the Cx26 channel and the implicit membrane. The boundaries of the implicit and explicit membranes are indicated by the dotted green lines. The simulation box (66 × 66 × 120 Å) was separated by the channel/membrane into two compartments, each containing 100 mM KCl. The buffer region maintains the chemical potential in the two compartments by creating and destroying ions with a GCMC algorithm. The upper compartment (extracellular part of the channel) was defined as the ground in voltage applications. 20 replicate 450-ns simulations were performed at each of seven voltages, ±150, ±100, ±50, and 0 mV, to plot the I-V relations. Blue circles, K+; red circles, Cl−. (B) Positions of modified residues identified by Locke et al. (2009) that would alter the distribution of charge in the Cx26 channel pore, shown in a side view of two opposite subunits of the completed crystal structure. The positions of acetylated residues are colored as follows: blue, Met1; red, K15; green, K102, K103, K105, K108, K112, and K116 in CL/TM2; orange, γ-carboxyglutamated residues E42, E47, and E114.
	Figure 5. GCMC/BD simulations of the average equilibrated structure. (A) The I-V relation obtained computationally for the multi-ion case (100 mM of symmetric KCl). The black line is the total current, the red line is the Cl− current, the blue line is the K+ current, and the green line is the total current of the completed crystal structure corresponding to Fig. S5 A. (B) The pore radius of the average equilibrated structure is shown in red, and the completed crystal structure is in blue. The positions of Cα of Met1, the amide group of K41, the carbonyl of E42, and the Cα of A49 are indicated. (C) The PMF of K+ at three voltages, 150 mV (blue), 0 mV (green), and −150 mV (red), are plotted against the z coordinate. The blue, green, and red lines are the PMF determined at 150, 0, and −150 mV, respectively. (D) The PMF of Cl− determined from simulations at the same three voltages.
	Figure 6. Superimposition of experimental I-V relations (gray traces) and those computed by GCMC/BD simulations with the charge effects of the indicated protein modifications. (A) The correspondence between slight outward rectification (Fig. 3 E) and simulated channel with neutralization of Met1 and six internal lysine residues in TM2/CL. The charge neutralizations mimic acetylation. The black line is the total current, the blue line is the K+ current, the red line is the Cl− current, and the gray line is the single-channel current trace. (B) The pore radius of the simulated channel structure plotted against the z coordinate. The positions of the modified atoms are shown. The six acetylated lysine residues at the tip of TM2/CL are grouped into one position and labeled “CL-Acet 6x.” (C) The PMF of K+ at three voltages, 150 mV (blue), 0 mV (green), and −150 mV (red), plotted against the z coordinate. (D) The PMF of Cl− determined at the same three voltages. (E–H) Outwardly rectifying experimental I-V relation corresponds to simulated I-V relations in which Met1, K15, and six lysine residues in TM2/CL are neutral-mimicking acetylation. (I–L) Inward sigmoidal experimental I-V relation corresponds to simulated I-V relations in which the six lysine residues in TM2/CL are neutralized (acetylated) and the three glutamate residues are modified by γ-carboxyglutamation (Gla). (M–P) Slightly inward rectifying experimental I-V relations correspond to simulated I-V relations in which Met1, K15, and K116 are neutralized (acetylated) and the three glutamate residues are modified by γ-carboxyglutamation (Gla).

References [+] :

Allen, Atomic force microscopy of Connexin40 gap junction hemichannels reveals calcium-dependent three-dimensional molecular topography and open-closed conformations of both the extracellular and cytoplasmic faces. 2011, Pubmed
Arnesen, Proteomics analyses reveal the evolutionary conservation and divergence of N-terminal acetyltransferases from yeast and humans. 2009, Pubmed
Arnesen, Towards a functional understanding of protein N-terminal acetylation. 2011, Pubmed
Bas, Very fast prediction and rationalization of pKa values for protein-ligand complexes. 2008, Pubmed
Basu, Proteome-wide prediction of acetylation substrates. 2009, Pubmed
Berns, Calf crystallin synthesis in frog cells: the translation of lens-cell 14S RNA in oocytes. 1972, Pubmed , Xenbase
Bradshaw, N-terminal processing: the methionine aminopeptidase and N alpha-acetyl transferase families. 1998, Pubmed
Brooks, CHARMM: the biomolecular simulation program. 2009, Pubmed
Brown, Precursors of novel Gla-containing conotoxins contain a carboxy-terminal recognition site that directs gamma-carboxylation. 2005, Pubmed
Buck, Importance of the CMAP correction to the CHARMM22 protein force field: dynamics of hen lysozyme. 2006, Pubmed
Buffy, Solid-state NMR investigation of the selective perturbation of lipid bilayers by the cyclic antimicrobial peptide RTD-1. 2004, Pubmed
Caves, Locally accessible conformations of proteins: multiple molecular dynamics simulations of crambin. 1998, Pubmed
Chang, Gap junction mediated intercellular metabolite transfer in the cochlea is compromised in connexin30 null mice. 2008, Pubmed
Choudhary, Lysine acetylation targets protein complexes and co-regulates major cellular functions. 2009, Pubmed
Drag, Aminopeptidase fingerprints, an integrated approach for identification of good substrates and optimal inhibitors. 2010, Pubmed
Egwolf, Ion selectivity of alpha-hemolysin with beta-cyclodextrin adapter. II. Multi-ion effects studied with grand canonical Monte Carlo/Brownian dynamics simulations. 2010, Pubmed
Fleishman, A Calpha model for the transmembrane alpha helices of gap junction intercellular channels. 2004, Pubmed
Forte, N-terminal acetylation inhibits protein targeting to the endoplasmic reticulum. 2011, Pubmed
Frottin, The proteomics of N-terminal methionine cleavage. 2006, Pubmed
Furie, Vitamin K-dependent biosynthesis of gamma-carboxyglutamic acid. 1999, Pubmed
Goetze, Identification and functional characterization of N-terminally acetylated proteins in Drosophila melanogaster. 2009, Pubmed
Gong, Size selectivity between gap junction channels composed of different connexins. 2001, Pubmed , Xenbase
González, Species specificity of mammalian connexin-26 to form open voltage-gated hemichannels. 2006, Pubmed , Xenbase
Grossfield, Convergence of molecular dynamics simulations of membrane proteins. 2007, Pubmed
Harris, Kinetic properties of a voltage-dependent junctional conductance. 1981, Pubmed
Helbig, Profiling of N-acetylated protein termini provides in-depth insights into the N-terminal nature of the proteome. 2010, Pubmed
Hertzberg, Topology of the Mr 27,000 liver gap junction protein. Cytoplasmic localization of amino- and carboxyl termini and a hydrophilic domain which is protease-hypersensitive. 1988, Pubmed
Hoover, Canonical dynamics: Equilibrium phase-space distributions. 1985, Pubmed
Hwang, N-terminal acetylation of cellular proteins creates specific degradation signals. 2010, Pubmed
Im, A Grand Canonical Monte Carlo-Brownian dynamics algorithm for simulating ion channels. 2000, Pubmed
Ingólfsson, Gramicidin A backbone and side chain dynamics evaluated by molecular dynamics simulations and nuclear magnetic resonance experiments. I: molecular dynamics simulations. 2011, Pubmed
Ishida, PrDOS: prediction of disordered protein regions from amino acid sequence. 2007, Pubmed
Jo, CHARMM-GUI: a web-based graphical user interface for CHARMM. 2008, Pubmed
Jo, CHARMM-GUI Membrane Builder for mixed bilayers and its application to yeast membranes. 2009, Pubmed
Jo, Automated builder and database of protein/membrane complexes for molecular dynamics simulations. 2007, Pubmed
Kalmatsky, Structural studies of the N-terminus of Connexin 32 using 1H NMR spectroscopy. 2009, Pubmed
Kienker, Charge selectivity of the designed uncharged peptide ion channel Ac-(LSSLLSL)3-CONH2. 1995, Pubmed
Klauda, Update of the CHARMM all-atom additive force field for lipids: validation on six lipid types. 2010, Pubmed
Kronengold, Single-channel SCAM identifies pore-lining residues in the first extracellular loop and first transmembrane domains of Cx46 hemichannels. 2003, Pubmed , Xenbase
Kucerka, The effect of cholesterol on short- and long-chain monounsaturated lipid bilayers as determined by molecular dynamics simulations and X-ray scattering. 2008, Pubmed
Lee, Brownian dynamics simulations of ion transport through the VDAC. 2011, Pubmed
Li, Very fast empirical prediction and rationalization of protein pKa values. 2005, Pubmed
Li, Entropy localization in proteins. 2010, Pubmed
Liu, Nanomechanics of hemichannel conformations: connexin flexibility underlying channel opening and closing. 2006, Pubmed
Locke, Post-translational modifications of connexin26 revealed by mass spectrometry. 2009, Pubmed
MacKerell, All-atom empirical potential for molecular modeling and dynamics studies of proteins. 1998, Pubmed
Mackerell, Extending the treatment of backbone energetics in protein force fields: limitations of gas-phase quantum mechanics in reproducing protein conformational distributions in molecular dynamics simulations. 2004, Pubmed
Maeda, Structure of the connexin 26 gap junction channel at 3.5 A resolution. 2009, Pubmed
Müller, Conformational changes in surface structures of isolated connexin 26 gap junctions. 2002, Pubmed
Nicholson, Rat liver gap junction protein: properties and partial sequence. 1981, Pubmed
Noskov, Ion permeation through the alpha-hemolysin channel: theoretical studies based on Brownian dynamics and Poisson-Nernst-Plank electrodiffusion theory. 2004, Pubmed
Oh, Charges dispersed over the permeation pathway determine the charge selectivity and conductance of a Cx32 chimeric hemichannel. 2008, Pubmed , Xenbase
Oh, Changes in permeability caused by connexin 32 mutations underlie X-linked Charcot-Marie-Tooth disease. 1997, Pubmed , Xenbase
Oh, Molecular determinants of electrical rectification of single channel conductance in gap junctions formed by connexins 26 and 32. 1999, Pubmed
Olsson, PROPKA3: Consistent Treatment of Internal and Surface Residues in Empirical pKa Predictions. 2011, Pubmed
Oshima, Three-dimensional structure of a human connexin26 gap junction channel reveals a plug in the vestibule. 2007, Pubmed
Phillips, Scalable molecular dynamics with NAMD. 2005, Pubmed
Polevoda, A synopsis of eukaryotic Nalpha-terminal acetyltransferases: nomenclature, subunits and substrates. 2009, Pubmed
Purnick, Structure of the amino terminus of a gap junction protein. 2000, Pubmed , Xenbase
Rathmell, Biochemistry. A glucose-to-gene link. 2009, Pubmed
Roux, Theoretical and computational models of biological ion channels. 2004, Pubmed
Rui, Molecular dynamics studies of ion permeation in VDAC. 2011, Pubmed
Sadoul, The tale of protein lysine acetylation in the cytoplasm. 2011, Pubmed
Shearer, Posttranslational modifications in lens fiber connexins identified by off-line-HPLC MALDI-quadrupole time-of-flight mass spectrometry. 2008, Pubmed
Sherman, Methionine or not methionine at the beginning of a protein. 1985, Pubmed
Smart, HOLE: a program for the analysis of the pore dimensions of ion channel structural models. 1996, Pubmed
Starheim, Composition and biological significance of the human Nalpha-terminal acetyltransferases. 2009, Pubmed
Suchyna, Different ionic selectivities for connexins 26 and 32 produce rectifying gap junction channels. 1999, Pubmed , Xenbase
Sánchez, Differentially altered Ca2+ regulation and Ca2+ permeability in Cx26 hemichannels formed by the A40V and G45E mutations that cause keratitis ichthyosis deafness syndrome. 2010, Pubmed , Xenbase
Tang, Conformational changes in a pore-forming region underlie voltage-dependent "loop gating" of an unapposed connexin hemichannel. 2009, Pubmed , Xenbase
Trexler, The first extracellular loop domain is a major determinant of charge selectivity in connexin46 channels. 2000, Pubmed , Xenbase
Unger, Three-dimensional structure of a recombinant gap junction membrane channel. 1999, Pubmed
Verselis, Loop gating of connexin hemichannels involves movement of pore-lining residues in the first extracellular loop domain. 2009, Pubmed , Xenbase
Wellen, ATP-citrate lyase links cellular metabolism to histone acetylation. 2009, Pubmed
Woolf, Molecular dynamics simulation of the gramicidin channel in a phospholipid bilayer. 1994, Pubmed
Woolf, Structure, energetics, and dynamics of lipid-protein interactions: A molecular dynamics study of the gramicidin A channel in a DMPC bilayer. 1996, Pubmed
Yang, The diverse superfamily of lysine acetyltransferases and their roles in leukemia and other diseases. 2004, Pubmed
Yang, Lysine acetylation: codified crosstalk with other posttranslational modifications. 2008, Pubmed
Zhao, Distinct and gradient distributions of connexin26 and connexin30 in the cochlear sensory epithelium of guinea pigs. 2006, Pubmed
Zhao, Connexin26 is responsible for anionic molecule permeability in the cochlea for intercellular signalling and metabolic communications. 2005, Pubmed
Zhou, Identification of a pore lining segment in gap junction hemichannels. 1997, Pubmed , Xenbase
Zimmer, Topological analysis of the major protein in isolated intact rat liver gap junctions and gap junction-derived single membrane structures. 1987, Pubmed