Secondary Water Pore Formation for Proton Transport in a ClC Exchanger Revealed by an Atomistic Molecular-Dynamics Simulation

Several prokaryotic ClC proteins have been demonstrated to function as exchangers that transport both chloride ions and protons simultaneously in opposite directions. However, the path of the proton through the ClC exchanger, and how the protein brings about the coupled movement of both ions are still unknown. In this work, we use an atomistic molecular dynamics (MD) simulation to demonstrate that a previously unknown secondary water pore is formed inside an Escherichia coli ClC exchanger. The secondary water pore is bifurcated from the chloride ion pathway at E148. From the systematic simulations, we determined that the glutamate residue exposed to the intracellular solution, E203, plays an important role as a trigger for the formation of the secondary water pore, and that the highly conserved tyrosine residue Y445 functions as a barrier that separates the proton from the chloride ion pathways. Based on our simulation results, we conclude that protons in the ClC exchanger are conducted via a water network through the secondary water pore, and we propose a new mechanism for the coupled transport of chloride ions and protons. It has been reported that several members of ClC proteins are not just channels that simply transport chloride ions across lipid bilayers; rather, they are exchangers that transport both the chloride ion and proton in opposite directions. However, the ion transit pathways and the mechanism of the coupled movement of these two ions have not yet been unveiled. In this article, we report a new finding (to our knowledge) of a water pore inside a prokaryotic ClC protein as revealed by computer simulation. This water pore is bifurcated from the putative chloride ion, and water molecules inside the new pore connect two glutamate residues that are known to be key residues for proton transport. On the basis of our simulation results, we conclude that the water wire that is formed inside the newly found pore acts as a proton pathway, which enables us to resolve many problems that could not be addressed by previous experimental studies.     

Residues lining the inner pore vestibule of human muscle chloride channels

Chloride channels belonging to the ClC family are ubiquitous and participate in a wide variety of physiological and pathophysiological processes. To define sequence segments in ClC channels that contribute to the formation of their ion conduction pathway, we employed a combination of site-directed mutagenesis, heterologous expression, patch clamp recordings, and chemical modification of the human muscle ClC isoform, hClC-1. We demonstrate that a highly conserved 8-amino acid motif (P3) located in the linker between transmembrane domains D2 and D3 contributes to the formation of a wide pore vestibule facing the cell interior. Similar to a previously defined pore region (P1 region), this segment functionally interacts with the corresponding segment of the contralateral subunit. The use of cysteine-specific reagents of different size revealed marked differences in the diameter of pore-forming regions implying that ClC channels exhibit a pore architecture quite similar to that of certain cation channels, in which a narrow constriction containing major structural determinants of ion selectivity is neighbored by wide vestibules on both sides of the membrane.     

Anion permeation in human ClC-4 channels

ClC-4 and ClC-5 are mammalian ClC isoforms with unique ion conduction and gating properties. Macroscopic current recordings in heterologous expression systems revealed very small currents at negative potentials, whereas a substantially larger instantaneous current amplitude and a subsequent activation were observed upon depolarization. Neither the functional basis nor the physiological impact of these channel features are currently understood. Here, we used whole-cell recordings to study pore properties of human ClC-4 channels heterologously expressed in tsA201 or HEK293 cells. Variance analysis demonstrated that the prominent rectification of the instantaneous macroscopic current amplitude is due to a voltage-dependent unitary current conductance. The single channel amplitudes are very small, i.e., 0.10 +/- 0.02 pA at +140 mV for external Cl(-) and internal I(-). Conductivity and permeability sequences were determined for various external and internal anions, and both values increase for anions with lower dehydration energies. ClC-4 exhibits pore properties that are distinct from other ClC isoforms. These differences can be explained by assuming differences in the size of the pore narrowing and the electrostatic potentials within the ion conduction pathways.     

Bead-like passage of chloride ions through ClC chloride channels

The ClC chloride channels control the ionic composition of the cytoplasm and the volume of cells, and regulate electrical excitability. Recently, it has been proposed that prokaryotic ClC channels are H+-Cl- exchange transporter. Although X-ray and molecular dynamics (MD) studies of bacterial ClC channels have investigated the filter open-close and ion permeation mechanism of channels, details have remained unclear. We performed MD simulations of ClC channels involving H+, Na+, K+, or H3O+ in the intracellular region to elucidate the open-close mechanism, and to clarify the role of H+ ion an H+-Cl- exchange transporter. Our simulations revealed that H+ and Na+ caused channel opening and the passage of Cl- ions. Na+ induced a bead-like string of Cl- -Na+-Cl--Na+-Cl- ions to form and permeate through ClC channels to the intracellular side with the widening of the channel pathway.     

Cysteine accessibility in ClC-0 supports conservation of the ClC intracellular vestibule

ClC chloride channels, which are ubiquitously expressed in mammals, have a unique double-barreled structure, in which each monomer forms its own pore. Identification of pore-lining elements is important for understanding the conduction properties and unusual gating mechanisms of these channels. Structures of prokaryotic ClC transporters do not show an open pore, and so may not accurately represent the open state of the eukaryotic ClC channels. In this study we used cysteine-scanning mutagenesis and modification (SCAM) to screen >50 residues in the intracellular vestibule of ClC-0. We identified 14 positions sensitive to the negatively charged thiol-modifying reagents sodium (2-sulfonatoethyl)methanethiosulfonate (MTSES) or sodium 4-acetamido-4'-maleimidylstilbene-2'2-disulfonic acid (AMS) and show that 11 of these alter pore properties when modified. In addition, two MTSES-sensitive residues, on different helices and in close proximity in the prokaryotic structures, can form a disulfide bond in ClC-0. When mapped onto prokaryotic structures, MTSES/AMS-sensitive residues cluster around bound chloride ions, and the correlation is even stronger in the ClC-0 homology model developed by Corry et al. (2004). These results support the hypothesis that both secondary and tertiary structures in the intracellular vestibule are conserved among ClC family members, even in regions of very low sequence similarity.     

Ion Concentration- and Voltage-Dependent Push and Pull Mechanisms of Potassium Channel Ion Conduction

The mechanism of ion conduction by potassium channels is one of the central issues in physiology. In particular, it is still unclear how the ion concentration and the membrane voltage drive ion conduction. We have investigated the dynamics of the ion conduction processes in the Kv1.2 pore domain, by molecular dynamics (MD) simulations with several different voltages and ion concentrations. By focusing on the detailed ion movements through the pore including selectivity filter (SF) and cavity, we found two major conduction mechanisms, called the III-IV-III and III-II-III mechanisms, and the balance between the ion concentration and the voltage determines the mechanism preference. In the III-IV-III mechanism, the outermost ion in the pore is pushed out by a new ion coming from the intracellular fluid, and four-ion states were transiently observed. In the III-II-III mechanism, the outermost ion is pulled out first, without pushing by incoming ions. Increases in the ion concentration and voltage accelerated ion conductions, but their mechanisms were different. The increase in the ion concentrations facilitated the III-IV-III conductions, while the higher voltages increased the III-II-III conductions, indicating that the pore domain of potassium channels permeates ions by using two different driving forces: a push by intracellular ions and a pull by voltage.     

Anion pathway and potential energy profiles along curvilinear bacterial ClC Cl- pores: electrostatic effects of charged residues

X-ray structures permit theoretical study of Cl(-) permeation along bacterial ClC Cl(-) pores. We determined the lowest energy curvilinear pathway, identified anion-coordinating amino acids, and calculated the electrostatic potential energy profiles. We find that all four bacterial ClC Cl(-) crystal structures correspond to closed states. E148 and S107 side chains form steric barriers on both sides of the crystal binding site in the StClC wild-type and EcClC wild-type crystals; both the EcClC(E148A) and EcClC(E148Q) mutants are blocked at the S107 site. We studied the effect that mutating the charge of some strongly conserved pore-lining amino acids has on the electrostatic potential energy profiles. When E148 is neutralized, it creates an electrostatic trap, binding the ion near midmembrane. This suggests a possible electrostatic mechanism for controlling anion flow: neutralize E148, displace the side chain of E148 from the pore pathway to relieve the steric barrier, then trap the anion at midmembrane, and finally either deprotonate E148 and block the pore (pore closure) or bring a second Cl(-) into the pore to promote anion flow (pore conductance). Side-chain displacement may arise by competition for the binding site between the oxygens of E148 and the anion moving down the electrostatic energy gradient. We also find that the charge state of E111 and E113 may electrostatically control anion conductance and occupancy of the binding site within the cytoplasmic pore.     

Trapping single ions inside single ion channels.

Single Ca++-activated K+ channels from rat muscle plasma membranes are inhibited by Ba++. A single Ba++ entering the channel's conduction pore induces a long-lived blocked state. This study employs Ba++ as a probe of the channel's conduction pathway to show that the channel can be forced to close with a single Ba++ ion inside the pore. A Ba++ ion inside the closed channel is trapped and cannot escape until the channel opens. The results demonstrate that in the channel's closed state, the cytoplasmic side of the conduction pore is obstructed to the passage of ions.     

Mechanism of gating and ion conductivity of a possible tetrameric pore in aquaporin-1

While substrate permeation through monomeric pores of aquaporins is well characterized, little is known about the possible tetrameric pore. AQP1 has been suggested to function as an ion channel upon cGMP activation, although this idea has been controversial. Taking a theoretical and experimental approach, we demonstrate that the current might arise through the tetrameric pore and propose a plausible mechanism for conduction and gating. In response to simulated ion permeation, immediate hydration of the putative central pore was facilitated by moderate conformational changes of pore-lining residues. cGMP is found to interact with an unusually arginine-rich, cytoplasmic loop (loop D) facilitating its outward motion, which is hypothesized to trigger the opening of a cytoplasmic gate. Physiological analyses of wild-type AQP1 and a designed mutant in which two arginines of the gating loop are replaced by alanine provide experimental support for identifying a key component of the proposed mechanism.     

Ion selectivity mechanism in a bacterial pentameric ligand-gated ion channel

The proton-gated ion channel from Gloeobacter violaceus (GLIC) is a prokaryotic homolog of the eukaryotic nicotinic acetylcholine receptor that responds to the binding of neurotransmitter acetylcholine and mediates fast signal transmission. Recent emergence of a high-resolution crystal structure of GLIC captured in a potentially open state allowed detailed, atomic-level insight into ion conduction and selectivity mechanisms in these channels. Herein, we have examined the barriers to ion conduction and origins of ion selectivity in the GLIC channel by the construction of potential-of-mean-force profiles for sodium and chloride ions inside the transmembrane region. Our calculations reveal that the GLIC channel is open for a sodium ion to transport, but presents a ∼11 kcal/mol free energy barrier for a chloride ion. Our collective findings identify three distinct contributions to the observed preference for the permeant ions. First, there is a substantial contribution due to a ring of negatively charged glutamate residues (E-2′) at the narrow intracellular end of the channel. The negative electrostatics of this region and the ability of the glutamate side chains to directly bind cations would strongly favor the passage of sodium ions while hindering translocation of chloride ions. Second, our results imply a significant hydrophobic contribution to selectivity linked to differences in the desolvation penalty for the sodium versus chloride ions in the central hydrophobic region of the pore. This hydrophobic contribution is evidenced by the large free energy barriers experienced by Cl⁻ in the middle of the pore for both GLIC and the E-2′A mutant. Finally, there is a distinct contribution arising from the overall negative electrostatics of the channel.     

Structural basis for ion conduction and gating in ClC chloride channels

Members of the ClC family of voltage-gated chloride channels are found from bacteria to mammals with a considerable degree of conservation in the membrane-inserted, pore-forming region. The crystal structures of the ClC channels of Escherichia coli and Salmonella typhimurium provide a structural framework for the entire family. The ClC channels are homodimeric proteins with an overall rhombus-like shape. Each ClC dimer has two pores each contained within a single subunit. The ClC subunit consists of two roughly repeated halves that span the membrane with opposite orientations. This antiparallel architecture defines a chloride selectivity filter within the 15-A neck of a hourglass-shaped pore. Three Cl(-) binding sites within the selectivity filter stabilize ions by interactions with alpha-helix dipoles and by chemical interactions with nitrogen atoms and hydroxyl groups of residues in the protein. The Cl(-) binding site nearest the extracellular solution can be occupied either by a Cl(-) ion or by a glutamate carboxyl group. Mutations of this glutamate residue in Torpedo ray ClC channels alter gating in electrophysiological assays. These findings reveal a form of gating in which the glutamate carboxyl group closes the pore by mimicking a Cl(-) ion.     

Ion transit pathways and gating in ClC chloride channels

ClC chloride channels possess a homodimeric structure in which each monomer contains an independent chloride ion pathway. ClC channel gating is regulated by chloride ion concentration, pH and voltage. Based on structural and physiological evidence, it has been proposed that a glutamate residue on the extracellular end of the selectivity filter acts as a fast gate. We utilized a new search algorithm that incorporates electrostatic information to explore the ion transit pathways through wild-type and mutant bacterial ClC channels. Examination of the chloride ion permeation pathways supports the importance of the glutamate residue in gating. An external chloride binding site previously postulated in physiological experiments is located near a conserved basic residue adjacent to the gate. In addition, access pathways are found for proton migration to the gate, enabling pH control at hyperpolarized membrane potentials. A chloride ion in the selectivity filter is required for the pH-dependent gating mechanism.     

Exploring gas permeability of cellular membranes and membrane channels with molecular dynamics

Aquaporins are a family of membrane proteins specialized in rapid water conduction across biological membranes. Whether these channels also conduct gas molecules and the physiological significance of this potential function have not been well understood. Here we report 140 ns of molecular dynamics simulations of membrane-embedded AQP1 and of a pure POPE bilayer addressing these questions. The permeability of AQP1 to two types of gas molecules, O2 and CO2, was investigated using two complementary methods, namely, explicit gas diffusion simulation and implicit ligand sampling. The simulations show that the central (tetrameric) pore of AQP1 can be readily used by either gas molecule to permeate the channel. The two approaches produced similar free energy profiles associated with gas permeation through the central pore: a -0.4 to -1.7 kcal/mol energy well in the middle, and a 3.6-4.6 kcal/mol energy barrier in the periplasmic vestibule. The barrier appears to be mainly due to a dense cluster of water molecules anchored in the periplasmic mouth of the central pore by four aspartate residues. Water pores show a very low permeability to O2, but may contribute to the overall permeation of CO2 due to its more hydrophilic nature. Although the central pore of AQP1 is found to be gas permeable, the pure POPE bilayer provides a much larger cross-sectional area, thus exhibiting a much lower free energy barrier for CO2 and O2 permeation. As such, gas conduction through AQP1 may only be physiologically relevant either in membranes of low gas permeability, or in cells where a major fraction of the cellular membrane is occupied by AQPs.     

The fast gating mechanism in ClC-0 channels

We investigate and then modify the hypothesis that a glutamate side chain acts as the fast gate in ClC-0 channels. We first create a putative open-state configuration of the prokaryotic ClC Cl- channel using its crystallographic structure as a basis. Then, retaining the same pore shape, the prokaryotic ClC channel is converted to ClC-0 by replacing all the nonconserved polar and charged residues. Using this open-state channel model, we carry out molecular dynamics simulations to study how the glutamate side chain can move between open and closed configurations. When the side chain extends toward the extracellular end of the channel, it presents an electrostatic barrier to Cl- conduction. However, external Cl- ions can push the side chain into a more central position where, pressed against the channel wall, it does not impede the motion of Cl- ions. Additionally, a proton from a low-pH external solution can neutralize the extended glutamate side chain, which also removes the barrier to conduction. Finally, we use Brownian dynamics simulations to demonstrate the influence of membrane potential and external Cl- concentration on channel open probability.     

Exterior site occupancy infers chloride-induced proton gating in a prokaryotic homolog of the ClC chloride channel

The ClC family of anion channels mediates the efficient, selective permeation of Cl(-) across the biological membranes of living cells under the driving force of an electrochemical gradient. In some eukaryotes, these channels are known to exhibit a unique gating mechanism, which appears to be triggered by the permeant Cl(-) anion. We infer details of this gating mechanism by studying the free energetics of Cl(-) occupancy in the pore of a prokaryotic ClC homolog. These free energetics were gleaned from 30 ns of molecular dynamics simulation on an approximately 133,000-atom system consisting of a hydrated membrane embedded StClC transporter. The binding sites for Cl(-) in the transporter were determined for the cases where the putative gating residue, Glu(148), was protonated and unprotonated. When the glutamate gate is protonated, Cl(-) favorably occupies an exterior site, S(ext), to form a queue of anions in the pore. However, when the glutamate gate is unprotonated, Cl(-) cannot occupy this site nor, consequently, pass through the pore. An additional, previously undetected, site was found in the pore near the outer membrane that exists regardless of the protonation state of Glu(148). Although this suggests that, for the prokaryotic homolog, protonation of Glu(148) may be the first step in transporting Cl(-) at the expense of H(+) transport in the opposite direction, an evolutionary argument might suggest that Cl(-) opens the ClC gate in eukaryotic channels by inducing the conserved glutamate's protonation. During an additional 20 ns free dynamics simulation, the newly discovered outermost site, S(out), and the innermost site, S(int), were seen to allow spontaneous exchange of Cl(-) ions with the bulk electrolyte while under depolarization conditions.     

Ion passage pathways and thermodynamics of the amphotericin B membrane channel

Amphotericin B is a polyene macrolide antibiotic used to treat systemic fungal infections. Amphotericin B's chemotherapeutic action requires the formation of transmembrane channels, which are known to transmit monovalent ions. We have investigated the ion passage pathways through the pore of a realistic model structure of the channel and computed the associated thermodynamic properties. Our calculations combined the free energy computations using the Poisson equation with a continuum solvent model and the molecular simulations in which solvent molecules were present explicitly. It was found that there are no substantial structural barriers to a single sodium or chloride ion passage. Thermodynamic free energy calculations showed that the path along which the ions prefer to move is off center from the channel's central axis. In accordance with experiments, Monte Carlo molecular simulations established that sodium ions can pass through the pore. When it encounters a chloride anion in the channel, the sodium cation prefers to form a solvent-bridged pair configuration with the anion.     

Nitric oxide conduction by the brain aquaporin AQP4

Involvement of aquaporins in gas conduction across the membrane and the physiological significance of this process have attracted marked attention from both experimental and theoretical studies. Previous work demonstrated that AQP1 is permeable to both CO₂ and O₂. Here we employ various simulation techniques to examine the permeability of the brain aquaporin AQP4 to NO and O₂ and to describe energetics and pathways associated with these phenomena. The energy barrier to NO and O₂ permeation through AQP4 central pore is found to be only ～3 kcal mol^?1. The results suggest that the central pore of AQP4, similar to that of AQP1, can indeed conduct gas molecules. Interestingly, despite a longer and narrower central pore, AQP4 appears to provide an energetically more favorable permeation pathway for gas molecules than AQP1, mainly due to the different orientation of its charged residues near the pore entrance. Although the low barrier against gas permeation through AQP4 indicates that it can participate in gas conduction across the cellular membrane, physiological relevance of the phenomenon remains to be established experimentally, particularly since pure lipid bilayers appear to present a more favorable pathway for gas conduction across the membrane. With an energy well of ?1.8 kcal mol^?1, the central pore of AQP4 may also act as a reservoir for NO molecules to accumulate in the membrane. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

Single-Channel Current Through Nicotinic Receptor Produced by Closure of Binding Site C-Loop

We investigated the initial coupling of agonist binding to channel gating of the nicotinic acetylcholine receptor using targeted molecular-dynamics (TMD) simulation. After TMD simulation to accelerate closure of the C-loops at the agonist binding sites, the region of the pore that passes through the cell membrane expands. To determine whether the structural changes in the pore result in ion conduction, we used a coarse-grained ion conduction simulator, Biology Boltzmann transport Monte Carlo, and applied it to two structural frames taken before and after TMD simulation. The structural model before TMD simulation represents the channel in the proposed “resting” state, whereas the model after TMD simulation represents the channel in the proposed “active” state. Under external voltage biases, the channel in the “active” state was permeable to cations. Our simulated ion conductance approaches that obtained experimentally and recapitulates several functional properties characteristic of the nicotinic acetylcholine receptor. Thus, closure of the C-loop triggers a structural change in the channel sufficient to account for the open channel current. This approach of applying Biology Boltzmann transport Monte Carlo simulation can be used to further investigate the binding to gating transduction mechanism and the structural bases for ion selection and translocation.     

Ab initio molecular orbital study of the interaction of Li, Na and K with the pore components of ion channels: Consideration of the size, structure and selectivity of the pore of the channels

Ab initio molecular orbital calculations were made for the various types of structures of the pore of the ion channel and the results were applied to the permeability model by Hille, an extension of the Eyring rate theory. In Hille's model, ion passage through the channel is regarded as a kinetic process. Accordingly, it is thought that the interaction energy between cation and ligand, the easier the passage is, due to the lower activation energy. The calculated interaction energy was in the order Li+ greater than Na+ greater than K+ for all models. The optimum size of the pore determined from the interaction energy depends on the structure of the filter. The size for the pentagon was largest, followed by the hexagon and tetragon. On the other hand, the size depends hardly at all on the kind of ligand molecules. In the case of the tetragon, the sizes for the Na and K channels were nearly the same as those estimated from the model building and inhibitor-blocking experiment. The interaction energy between the ionized carboxyl group and the cation was extremely large, clearly reflecting the experimental fact that the carboxyl group in the pore has an important role in making the passage of the cation through the channel easier by dehydrating the water molecules. By analysis of the interaction energy, it was revealed that the contribution of the electrostatic energy was predominant, although the contributions of the other effects might not be negligible. Among these effects, the value of the charge transfer energy is largest, and this is noteworthy in connection with the selective transmission of cations through the overlap of orbitals. It is concluded that the quantum-chemical indices such as interaction energy and the electronic charge calculated by the sophisticated ab initio method help to shed light on the nature of the pore of the ion channel.     

Conduction mechanisms of chloride ions in ClC-type channels

The conduction properties of ClC-0 and ClC-1 chloride channels are examined using electrostatic calculations and three-dimensional Brownian dynamics simulations. We create an open-state configuration of the prokaryotic ClC Cl(-) channel using its known crystallographic structure as a basis. Two residues that are occluding the channel are slowly pushed outward with molecular dynamics to create a continuous ion-conducting path with the minimum radius of 2.5 A. Then, retaining the same pore shape, the prokaryotic ClC channel is converted to either ClC-0 or ClC-1 by replacing all the nonconserved dipole-containing and charged amino acid residues. Employing open-state ClC-0 and ClC-1 channel models, current-voltage curves consistent with experimental measurements are obtained. We find that conduction in these pores involves three ions. We locate the binding sites, as well as pinpointing the rate-limiting steps in conduction, and make testable predictions about how the single channel current across ClC-0 and ClC-1 will vary as the ionic concentrations are increased. Finally, we demonstrate that a ClC-0 homology model created from an alternative sequence alignment fails to replicate any of the experimental observations.