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

Electrostatic Potential of Functional Cations as a Predictor of Hydroxide Diffusion Pathways in Nanoconfined Environments of Anion Exchange Membranes

Nanoconfined anion exchange membranes (AEMs) play a vital role in emerging electrochemical technologies. The ability to control dominant hydroxide diffusion pathways is an important goal in the design of nanoconfined AEMs. Such control can shorten hydroxide transport pathways between electrodes, reduce transport resistance, and enhance device performance. In this work, we propose an electrostatic potential (ESP) approach to explore the effect of the polymer electrolyte cation spacing on hydroxide diffusion pathways from a molecular perspective. By exploring cation ESP energy surfaces and validating outcomes through prior ab initio molecular dynamics simulations of nanoconfined AEMs, we find that we can achieve control over preferred hydroxide diffusion pathways by adjusting the cation spacing. The results presented in this work provide a unique and straightforward approach to predict preferential hydroxide diffusion pathways, enabling efficient design of highly conductive nanoconfined AEM materials for electrochemical technologies.

CLC Chloride Channels and Transporters: Structure, Function, Physiology, and Disease

CLC anion transporters are found in all phyla and form a gene family of eight members in mammals. Two CLC proteins, each of which completely contains an ion translocation parthway, assemble to homo- or heteromeric dimers that sometimes require accessory β-subunits for function. CLC proteins come in two flavors: anion channels and anion/proton exchangers. Structures of these two CLC protein classes are surprisingly similar. Extensive structure-function analysis identified residues involved in ion permeation, anion-proton coupling and gating and led to attractive biophysical models. In mammals, ClC-1, -2, -Ka/-Kb are plasma membrane Cl−channels, whereas ClC-3 through ClC-7 are 2Cl−/H+-exchangers in endolysosomal membranes. Biological roles of CLCs were mostly studied in mammals, but also in plants and model organisms like yeast and Caenorhabditis elegans. CLC Cl−channels have roles in the control of electrical excitability, extra- and intracellular ion homeostasis, and transepithelial transport, whereas anion/proton exchangers influence vesicular ion composition and impinge on endocytosis and lysosomal function. The surprisingly diverse roles of CLCs are highlighted by human and mouse disorders elicited by mutations in their genes. These pathologies include neurodegeneration, leukodystrophy, mental retardation, deafness, blindness, myotonia, hyperaldosteronism, renal salt loss, proteinuria, kidney stones, male infertility, and osteopetrosis. In this review, emphasis is laid on biophysical structure-function analysis and on the cell biological and organismal roles of mammalian CLCs and their role in disease.

Chloride Ion Transport by the E. coli CLC Cl-/H+ Antiporter: A Combined Quantum-Mechanical and Molecular-Mechanical Study

We performed steered molecular dynamics (SMD) and umbrella sampling simulations of Cl^- ion migration through the transmembrane domain of a prototypical E. coli CLC Cl^-/H⁺ antiporter by employing combined quantum-mechanical (QM) and molecular-mechanical (MM) calculations. The SMD simulations revealed interesting conformational changes of the protein. While no large-amplitude motions of the protein were observed during pore opening, the side chain rotation of the protonated external gating residue Glu148 was found to be critical for full access of the channel entrance by Cl^-. Moving the anion into the external binding site (S_ext) induced small-amplitude shifting of the protein backbone at the N-terminal end of helix F. As Cl^- traveled through the pore, rigid-body swinging motions of helix R separated it from helix D. Helix R returned to its original position once Cl^- exited the channel. Population analysis based on polarized wavefunction from QM/MM calculations discovered significant (up to 20%) charge loss for Cl^- along the ion translocation pathway inside the pore. The delocalized charge was redistributed onto the pore residues, especially the functional groups containing π bonds (e.g., the Tyr445 side chain), while the charges of the H atoms coordinating Cl^- changed almost negligibly. Potentials of mean force computed from umbrella sampling at the QM/MM and MM levels both displayed barriers at the same locations near the pore entrance and exit. However, the QM/MM PMF showed higher barriers (~10 kcal/mol) than the MM PMF (~2 kcal/mol). Binding energy calculations indicated that the interactions between Cl^- and certain pore residues were overestimated by the semi-empirical PM3 Hamiltonian and underestimated by the CHARMM36 force fields, both of which were employed in the umbrella sampling simulations. In particular, CHARMM36 underestimated binding interactions for the functional groups containing π bonds, missing the stabilizations of the Cl^- ion due to electron delocalization. The results suggested that it is important to explore these quantum effects for accurate descriptions of the Cl^- transport.

Multiscale Kinetic Modeling Reveals an Ensemble of Cl-/H+ Exchange Pathways in ClC-ec1 Antiporter

Despite several years of research, the ion exchange mechanisms in chloride/proton antiporters and many other coupled transporters are not yet understood at the molecular level. Here, we present a novel approach to kinetic modeling and apply it to ion exchange in ClC-ec1. Our multiscale kinetic model is developed by (1) calculating the state-to-state rate coefficients with reactive and polarizable molecular dynamics simulations, (2) optimizing these rates in a global kinetic network, and (3) predicting new electrophysiological results. The model shows that the robust Cl:H exchange ratio (2.2:1) can indeed arise from kinetic coupling without large protein conformational changes, indicating a possible facile evolutionary connection to chloride channels. The E148 amino acid residue is shown to couple chloride and proton transport through protonation-dependent blockage of the central anion binding site and an anion-dependent pK_a value, which influences proton transport. The results demonstrate how an ensemble of different exchange pathways, as opposed to a single series of transitions, culminates in the macroscopic observables of the antiporter, such as transport rates, chloride/proton stoichiometry, and pH dependence.

Gating the glutamate gate of CLC-2 chloride channel by pore occupancy

Intracellular permeant anions, and not extracellular protons, are the predominant driver of fast gating in the hyperpolarization-activated CLC-2 chloride channel.

CLC-2 channels are dimeric double-barreled chloride channels that open in response to hyperpolarization. Hyperpolarization activates protopore gates that independently regulate the permeability of the pore in each subunit and the common gate that affects the permeability through both pores. CLC-2 channels lack classic transmembrane voltage–sensing domains; instead, their protopore gates (residing within the pore and each formed by the side chain of a glutamate residue) open under repulsion by permeant intracellular anions or protonation by extracellular H⁺. Here, we show that voltage-dependent gating of CLC-2: (a) is facilitated when permeant anions (Cl⁻, Br⁻, SCN⁻, and I⁻) are present in the cytosolic side; (b) happens with poorly permeant anions fluoride, glutamate, gluconate, and methanesulfonate present in the cytosolic side; (c) depends on pore occupancy by permeant and poorly permeant anions; (d) is strongly facilitated by multi-ion occupancy; (e) is absent under likely protonation conditions (pH_e = 5.5 or 6.5) in cells dialyzed with acetate (an impermeant anion); and (f) was the same at intracellular pH 7.3 and 4.2; and (g) is observed in both whole-cell and inside-out patches exposed to increasing [Cl⁻]_i under unlikely protonation conditions (pH_e = 10). Thus, based on our results we propose that hyperpolarization activates CLC-2 mainly by driving intracellular anions into the channel pores, and that protonation by extracellular H⁺ plays a minor role in dislodging the glutamate gate.

A tale of two CLCs: biophysical insights toward understanding ClC-5 and ClC-7 function in endosomes and lysosomes

The CLC protein family comprises both Cl(-) channels and H(+) -coupled anion transporters. The understanding of the critical role of CLC proteins in a number of physiological functions has greatly contributed to a revision of the classical paradigm that attributed to Cl(-) ions only a marginal role in human physiology. The endosomal ClC-5 and the lysosomal ClC-7 are the best characterized human CLC transporters. Their dysfunction causes Dent's disease and osteopetrosis, respectively. It had been originally proposed that they would provide a Cl(-) shunt conductance allowing efficient acidification of intracellular compartments. However, this model seems to conflict with the transport properties of these proteins and with recent physiological evidence. Currently, there is no consensus on their specific physiological role. CLC proteins present also a number of peculiar biophysical properties, such as the dimeric architecture, the co-existence of intrinsically different thermodynamic modes of transport based on similar structural principles, and the gating mechanism recently emerging for the transporters, just to name a few. This review focuses on the biophysical properties and physiological roles of ClC-5 and ClC-7.

Influences of mutations on the electrostatic binding free energies of chloride ions in Escherichia coli ClC

Mutations in ClC channel proteins may cause serious functional changes and even diseases. The function of ClC proteins mainly manifests as Cl(-) transport, which is related to the binding free energies of chloride ions. Therefore, the influence of a mutation on ClC function can be studied by investigating the mutational effect on the binding free energies of chloride ions. The present study provides quantitative and systematic investigations on the influences of residue mutations on the electrostatic binding free energies in Escherichia coli ClC (EcClC) proteins, using all-atom molecular dynamics simulations. It was found that the change of the electrostatic binding free energy decreases linearly with the increase of the residue-chloride ion distance for a mutation. This work reveals how changes in the charge of a mutated residue and in the distance between the mutated residue and the binding site govern the variations in the electrostatic binding free energies and therefore influence the transport of chloride ions and conduction in EcClC. This work would facilitate our understanding of the mutational effects on transport of chloride ions and functions of ClC proteins and provide a guideline to estimate which residue mutations will have great influences on ClC functions.

A three-state multi-ion kinetic model for conduction properties of ClC-0 chloride channel

A three-state, multiion kinetic model was proposed to enable the conduction properties of the mammalian channel ClC-0 to be well characterized. Using this rate-theory based model, the current-voltage and conductance-concentration relations were obtained. The five parameters needed were determined by fitting the data of conduction experiments of the wild-type ClC-0 and its K519C mutant. The model was then tested against available calculation and simulation data, and the energy differences between distinct chloride-occupancy states computed agreed with an independent calculation on the binding free energies solved by using the Poisson-Boltzmann equation. The average ion number of conduction and the ion passing duration calculated closely resembled the values obtained from Brownian dynamics simulations. According to the model, the decrease of conductance caused by mutating residue K519 to C519 can be attributed to the effect of K519C mutation on translocation rate constants. Our study sets up a theoretical model for ion permeation and conductance in ClC-0. It provides a starting point for experimentalists to test the three-state model, and would help in understanding the conduction mechanism of ClC-0.

Electroelastic coupling between membrane surface fluctuations and membrane-embedded charges: continuum multidielectric treatment

The coupling of electric fields and charges with membrane-water interfacial fluctuations affects membrane electroporation, ionic conductance, and voltage gating. A modified continuum model is introduced to study charge interaction with membrane-water interfacial fluctuations in multidielectric environments. By surrounding a point charge with a low dielectric sphere, the linear Poisson-Boltzmann equation is directly solved by calculating the reaction field potential via a method that eliminates singularity contributions. This allows treatment of charges located at dielectric boundaries. Two complementary mechanisms governing charge-fluctuation interactions are considered: (1) electroelastic deformation (EED), treating the membrane as an elastic slab (smectic bilayer model), and (2) electrohydrophobic solvation (EHS), accounting for water penetration into the membrane's hydrophobic core. EED often leads to large membrane thickness perturbations, far larger than those consistent with elastic model descriptions [M. B. Partenskii, G. V. Miloshevsky, and P. C. Jordan, Isr. J. Chem. 47, 385 (2007)]. We argue that a switch from EED to EHS can be energetically advantageous at intermediate perturbation amplitudes. Both perturbation mechanisms are simulated by introducing adjustable shapes optimized by the kinetic Monte Carlo reaction path following approach [G. V. Miloshevsky and P. C. Jordan, J. Chem. Phys. 122, 214901 (2005)]. The resulting energy profiles agree with those of recent atomistic molecular dynamics studies on translating a charged residue across a lipid bilayer [S. Dorairaj and T. W. Allen, Proc. Natl. Acad. Sci. U.S.A. 104, 4943 (2007)].

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.

Solvation free energies of glutamate and its metal complexes: a computer simulation study

Most biological ion channels demonstrate a high degree of selectivity for one type of ion more than others, and in many cases, how they control attaining this is still not clear. So we have studied on some metal ion compounds of glutamate. The Glutamate and its meal ion compounds (Ca(2+), Na(+), K(+) and Li(+)) were first modeled by ab initio calculations and then Monte Carlo simulation was used to calculate solvation free energies and also the complexes free energies for the related structures. The results indicated that Glutamate-Ca(2+) have more stability in water than other metal ion. Also, it was found out that the more movement in ions; less stability of the structure would result. This trend can be seen both in gas and liquid phase.

Chloride ion conduction without water coordination in the pore of ClC protein

In the present work, we have found by an atomistic molecular dynamics simulation that hydrogen atoms originating from the residues of a prokaryotic ClC protein (EcClC) stabilize the chloride ion without water molecules in the pore of ClC protein. When the chloride ion conduction is simulated by pulling a chloride ion along the pore axis, the free energy barrier for chloride ion conduction is calculated to be low (4 kcal/mol), although the chloride ion is stripped of its hydration shell as it passes through the dehydrated pore region. The calculation of the number of hydrogen atoms surrounding the chloride ion reveals that water molecules hydrating the chloride ion are replaced by polar and non-polar hydrogen atoms protruding from the protein residues. From the analysis of the pair interaction energy between the chloride ion and these hydrogen atoms, it is realized that the hydrogen atoms from the protein residues stabilize the chloride ion at the dehydrated region instead of water molecules, by which the energetic penalty for detaching water molecules from the permeating ion is compensated.

Shape-Dependent Global Deformation Modes of Large Protein Structures

Conformational changes are central to the functioning of pore-forming proteins that open and close their molecular gates in response to external stimuli such as pH, ionic strength, membrane voltage or ligand binding. Normal mode analysis (NMA) is used to identify and characterize the slowest motions in the gA, KcsA, ClC-ec1, LacY and LeuT(Aa) proteins at the onset of gating. Global deformation modes of the essentially cylindrical gA, KcsA, LacY and LeuT(Aa) biomolecules are reminiscent of global twisting, transverse and longitudinal motions in a homogeneous elastic rod. The ClC-ec1 protein executes a splaying motion in the plane perpendicular to the lipid bilayer. These global collective deformations are determined by protein shape. New methods, all-atom Monte Carlo Normal Mode Following and its simplification using a rotation-translation of protein blocks (RTB), are described and applied to gain insight into the nature of gating transitions in gA and KcsA. These studies demonstrate the severe limitations of standard NMA in characterizing the structural rearrangements associated with gating transitions. Comparison of all-atom and RTB transition pathways in gA clearly illustrates the impact of the rigid protein block approximation and the need to include all degrees of freedom and their relaxation in computational studies of protein gating. The effects of atomic level structure, pH, hydrogen bonding and charged residues on the large scale conformational changes associated with gating transitions are discussed.

Antiport mechanism for Cl(-)/H(+) in ClC-ec1 from normal-mode analysis

ClC chloride channels and transporters play major roles in cellular excitability, epithelial salt transport, volume, pH, and blood pressure regulation. One family member, ClC-ec1 from Escherichia coli, has been structurally resolved crystallographically and subjected to intensive mutagenetic, crystallographic, and electrophysiological studies. It functions as a Cl(-)/H(+) antiporter, not a Cl(-) channel; however, the molecular mechanism for Cl(-)/H(+) exchange is largely unknown. Using all-atom normal-mode analysis to explore possible mechanisms for this antiport, we propose that Cl(-)/H(+) exchange involves a conformational cycle of alternating exposure of Cl(-) and H(+) binding sites of both ClC pores to the two sides of the membrane. Both pores switch simultaneously from facing outward to facing inward, reminiscent of the standard alternating-access mechanism, which may have direct implications for eukaryotic Cl(-)/H(+) transporters and Cl(-) channels.

Discrete-state representation of ion permeation coupled to fast gating in a model of ClC chloride channels: comparison to multi-ion continuous space Brownian dynamics simulations

A discrete-state model of chloride ion motion in a ClC chloride channel is constructed, following a previously developed multi-ion continuous space model of the same system (Cheng, M. H.; Mamonov, A. B.; Dukes, J. W.; Coalson, R. D. J. Phys. Chem. B 2007, 111, 5956) that included a simplistic representation of the fast gate in this channel. The reducibility of the many-body continuous space to the eight discrete-state model considered in the present work is examined in detail by performing three-dimensional Brownian dynamics simulations of each allowed state-to-state transition in order to extract the appropriate rate constant for this process, and then inserting the pairwise rate constants thereby obtained into an appropriate set of kinetic master equations. Experimental properties of interest, including the rate of Cl(-) ion permeation through the open channel and the average rate of closing of the fast gate as a function of bulk Cl(-) ion concentrations in the intracellular and extracellular electrolyte reservoirs are computed. Good agreement is found between the results obtained via the eight discrete-state model versus the multi-ion continuous space model, thereby encouraging continued development of the discrete-state model to include more complex behaviors observed experimentally in these channels.

Proton transport pathway in the ClC Cl-/H+ antiporter

A fundamental question concerning the ClC Cl-/H+ antiporters is the nature of their proton transport (PT) pathway. We addressed this issue by using a novel computational methodology capable of describing the explicit PT dynamics in the ClC-ec1 protein. The main result is that the Glu203 residue delivers a proton from the intracellular solution to the core of ClC-ec1 via a rotation of its side chain and subsequent acid dissociation. After reorientation of the Glu203 side chain, a transient water-mediated PT pathway between Glu203 and Glu148 is established that is able to receive and translocate the proton via Grotthuss shuttling after deprotonation of Glu203. A molecular-dynamics simulation of an explicit hydrated excess proton in this pathway suggests that a negatively charged Glu148 and the central Cl- ion act together to drive H+ to the extracellular side of the membrane. This finding is consistent with the experimental result that Cl- binding to the central site facilitates the proton movement. A calculation of the PT free-energy barrier for the ClC-ec1 E203V mutant also supports the proposal that a dissociable residue is required at this position for efficient delivery of H+ to the protein interior, in agreement with recent experimental results.

Voltage-dependent and -independent titration of specific residues accounts for complex gating of a ClC chloride channel by extracellular protons

The ClC transport protein family comprises both Cl(-) ion channel and H(+)/Cl(-) and H(+)/NO(3)(-) exchanger members. Structural studies on a bacterial ClC transporter reveal a pore obstructed at its external opening by a glutamate side-chain which acts as a gate for Cl(-) passage and in addition serves as a staging post for H(+) exchange. This same conserved glutamate acts as a gate to regulate Cl(-) flow in ClC channels. The activity of ClC-2, a genuine Cl(-) channel, has a biphasic response to extracellular pH with activation by moderate acidification followed by abrupt channel closure at pH values lower than approximately 7. We have now investigated the molecular basis of this complex gating behaviour. First, we identify a sensor that couples extracellular acidification to complete closure of the channel. This is extracellularly-facing histidine 532 at the N-terminus of transmembrane helix Q whose neutralisation leads to channel closure in a cooperative manner. We go on to show that acidification-dependent activation of ClC-2 is voltage dependent and probably mediated by protonation of pore gate glutamate 207. Intracellular Cl(-) acts as a voltage-independent modulator, as though regulating the pK(a) of the protonatable residue. Our results suggest that voltage dependence of ClC-2 is given by hyperpolarisation-dependent penetration of protons from the extracellular side to neutralise the glutamate gate deep within the channel, which allows Cl(-) efflux. This is reminiscent of a partial exchanger cycle, suggesting that the ClC-2 channel evolved from its transporter counterparts.

TransPath: a computational method for locating ion transit pathways through membrane proteins

The finely tuned structures of membrane channel proteins allow selective passage of ions through the available aqueous pores. To understand channel function, it is crucial to locate the pores and study their physical and chemical properties. Here, we propose a new pore-searching algorithm (TransPath), which uses the Configurational Bias Monte Carlo (CBMC) method to generate transmembrane trajectories driven by both geometric and electrostatic features. The trajectories are binned into groups determined by a vector distance criterion. From each group, a representative trajectory is selected based on the Rosenbluth weight, and the geometrically optimal path is obtained by simulated annealing. Candidate ion pathways can then be determined by analysis of the radius and potential profiles. The proposed method and its implementation are illustrated using the bacterial KcsA potassium channel as an example. The procedure is then applied to the more complex structures of the bacterial E. coli chloride channel homolog and a homology model of the ClC-0 channel.

Conformational changes in the selectivity filter of the open-state KcsA channel: an energy minimization study

Potassium channels switch between closed and open conformations and selectively conduct K(+) ions. There are at least two gates. The TM2 bundle at the intracellular site is the primary gate of KcsA, and rearrangements at the selectivity filter (SF) act as the second gate. The SF blocks ion flow via an inactivation process similar to C-type inactivation of voltage-gated K(+) channels. We recently generated the open-state conformation of the KcsA channel. We found no major, possibly inactivating, structural changes in the SF associated with this massive inner-pore rearrangement, which suggests that the gates might act independently. Here we energy-minimize the open state of wild-type and mutant KcsA, validating in silico structures of energy-minimized SFs by comparison with crystallographic structures, and use these data to gain insight into how mutation, ion depletion, and K(+) to Na(+) substitution influence SF conformation. Both E71 or D80 protonations/mutations and the presence/absence of protein-buried water molecule(s) modify the H-bonding network stabilizing the P-loops, spawning numerous SF conformations. We find that the inactivated state corresponds to conformations with a partially unoccupied or an entirely empty SF. These structures, involving modifications in all four P-loops, are stabilized by H-bonds between amide H and carbonyl O atoms from adjacent P-loops, which block ion passage. The inner portions of the P-loops are more rigid than the outer parts. Changes are localized to the outer binding sites, with innermost site S4 persisting in the inactivated state. Strong binding by Na(+) locally contracts the SF around Na(+), releasing ligands that do not participate in Na(+) coordination, and occluding the permeation pathway. K(+) selectivity primarily appears to arise from the inability of the SF to completely dehydrate Na(+) ions due to basic structural differences between liquid water and the "quasi-liquid" SF matrix.

CLC-0 and CFTR: Chloride Channels Evolved From Transporters

CLC-0 and cystic fibrosis transmembrane conductance regulator (CFTR) Cl−channels play important roles in Cl−transport across cell membranes. These two proteins belong to, respectively, the CLC and ABC transport protein families whose members encompass both ion channels and transporters. Defective function of members in these two protein families causes various hereditary human diseases. Ion channels and transporters were traditionally viewed as distinct entities in membrane transport physiology, but recent discoveries have blurred the line between these two classes of membrane transport proteins. CLC-0 and CFTR can be considered operationally as ligand-gated channels, though binding of the activating ligands appears to be coupled to an irreversible gating cycle driven by an input of free energy. High-resolution crystallographic structures of bacterial CLC proteins and ABC transporters have led us to a better understanding of the gating properties for CLC and CFTR Cl−channels. Furthermore, the joined force between structural and functional studies of these two protein families has offered a unique opportunity to peek into the evolutionary link between ion channels and transporters. A promising byproduct of this exercise is a deeper mechanistic insight into how different transport proteins work at a fundamental level.

Proton pathways and H+/Cl- stoichiometry in bacterial chloride transporters

H+/Cl- antiport behavior has recently been observed in bacterial chloride channel homologs and eukaryotic CLC-family proteins. The detailed molecular-level mechanism driving the stoichiometric exchange is unknown. In the bacterial structure, experiments and modeling studies have identified two acidic residues, E148 and E203, as key sites along the proton pathway. The E148 residue is a major component of the fast gate, and it occupies a site crucial for both H+ and Cl- transport. E203 is located on the intracellular side of the protein; it is vital for H+, but not Cl-, transport. This suggests two independent ion transit pathways for H+ and Cl- on the intracellular side of the transporter. Previously, we utilized a new pore-searching algorithm, TransPath, to predict Cl- and H+ ion pathways in the bacterial ClC channel homolog, focusing on proton access from the extracellular solution. Here we employ the TransPath method and molecular dynamics simulations to explore H+ pathways linking E148 and E203 in the presence of Cl- ions located at the experimentally observed binding sites in the pore. A conclusion is that Cl- ions are required at both the intracellular (S(int)) and central (S(cen)) binding sites in order to create an electrostatically favorable H+ pathway linking E148 and E203; this electrostatic coupling is likely related to the observed 1H+/2Cl- stoichiometry of the antiporter. In addition, we suggest that a tyrosine residue side chain (Y445), located near the Cl- ion binding site at S(cen), is involved in proton transport between E148 and E203.

The role of a conserved lysine in chloride- and voltage-dependent ClC-0 fast gating

ClC-0 is a chloride channel whose gating is sensitive to voltage, chloride, and pH. In a previous publication, we showed that the K149C mutation causes a +70-mV shift in the voltage dependence of ClC-0 fast gating. In this paper we analyze the effects of a series of mutations at K149 on the voltage and chloride dependence of gating. By fitting our data to the previously proposed four-state model for ClC-0 fast gating, we show which steps in fast-gate opening are likely to be affected by these mutations. Computational analysis of mutant ClC-0 homology models show electrostatic contributions to chloride binding that may partially account for the effects of K149 on gating. The analysis of gating kinetics in combination with the available structural information suggests some of the structural changes likely to underpin fast-gate opening.

Permeation and gating in proteins: kinetic Monte Carlo reaction path following

We present a new Monte Carlo technique, kinetic Monte Carlo reaction path following (kMCRPF), for the computer simulation of permeation and large-scale gating transitions in protein channels. It combines ideas from Metropolis Monte Carlo (MMC) and kinetic Monte Carlo (kMC) algorithms, and is particularly suitable when a reaction coordinate is well defined. Evolution of transition proceeds on the reaction coordinate by small jumps (kMC technique) toward the nearest lowest-energy uphill or downhill states, with the jumps thermally activated (constrained MMC). This approach permits navigation among potential minima on an energy surface, finding the minimum-energy paths and determining their associated free-energy profiles. The methodological and algorithmic strategies underlying the kMCRPF method are described. We have tested it using an analytical model and applied it to study permeation through the curvilinear ClC chloride and aquaporin pores and to gating in the gramicidin A channel. These studies of permeation and gating in real proteins provide extensive procedural tests of the method.

Ion permeation dynamics in carbon nanotubes

Molecular dynamics simulations are carried out to investigate the permeation of ions and water in a membrane consisting of single wall carbon nanotubes possessing no surface charges connecting two reservoirs. Our simulations reveal that there are changes in the first hydration shell of the ions upon confinement in tubes of 0.82 or 0.90 nm effective internal diameter. Although the first minimum in the g(r) is barely changed in the nanotube compared to in the bulk solution, the hydration number of Na(+) ion is reduced by 1.0 (from 4.5 in bulk to 3.5 in the 0.90 nm tube) and the hydration number is reduced further in the 0.82 nm tube. The changes in the hydration shell of Cl(-) ion are negligible, within statistical errors. The water molecules of the first hydration shell of both ions exchange less frequently inside the tube than in the bulk solution. We compare ion trajectories for ions in the same tube under identical reservoir conditions but with different numbers of ions in the tubes. This permits investigation of changes in structure and dynamics which arise from multiple ion occupancy in a carbon nanotube possessing no surface charges. We also investigated the effects of tube flexibility. Ions enter the tubes so as to form a train of ion pairs. We find that the radial distribution profiles of Na(+) ions broaden significantly systematically with increasing number of ion pairs in the tube. The radial distribution profiles of Cl(-) ions change only slightly with increasing number of ions in the tube. Trajectories reveal that Na(+) ions do not pass each other in 0.90 nm tubes, while Cl(-) ions pass each other, as do ions of opposite charge. An ion entering the tube causes the like-charged ions preceding it in the tube to be displaced along the tube axis and positive or negative ions will exit the tube only when one or two other ions of the same charge are present in the tube. Thus, the permeation mechanism involves multiple ions and Coulomb repulsion among the ions plays an essential role.

Roles of K149, G352, and H401 in the channel functions of ClC-0: testing the predictions from theoretical calculations

The ClC family of Cl⁻ channels and transporters comprises membrane proteins ubiquitously present in species ranging from prokaryotes to mammals. The recently solved structures of the bacterial ClC proteins have provided a good model to guide the functional experiments for the eukaryotic Cl⁻ channels. Theoretical calculations based on the bacterial ClC structures have identified several residues critical for the Cl⁻ binding energy in the Cl⁻ transport pathway. It was speculated that the corresponding residues in eukaryotic Cl⁻ channels might play similar roles for the channel functions. In this study, we made a series of mutations in three such residues in eukaryotic ClC Cl⁻ channels (K149, G352, and H401 in ClC-0) and studied the functional consequences on the channel properties. A cysteine modification approach was also employed to evaluate the electrostatic effects of the charge placed at these three positions. The experimental results revealed that among the three residues tested, K149 plays the most important role in controlling both the gating and the permeation functions of ClC-0. On the other hand, mutations of H401 alter the channel conductance but not the gating properties, while mutations of G352 result in very little functional consequence. The mutation of K149 into a neutral residue leucine (K149L) shifts the activation curve and leads to flickery channel openings. The anion permeability ratios derived from bi-ionic experiments are also significantly altered in that the selectivity of Cl⁻ over other anions is decreased. Furthermore, removing the positive charge at this position reduces and increases, respectively, the accessibility of the negatively and positively charged methane thiosulfonate reagents to the pore. The control of the accessibility to charged MTS reagents and the regulation of the anion permeation support the idea that K149 exerts an electrostatic effect on the channel function, confirming the prediction from computational studies.

The ClC family of chloride channels and transporters

The ClC proteins are members of a large family of chloride transport proteins, which are involved in a variety of physiological processes. All family members share a conserved molecular architecture consisting of a complex transmembrane transport domain and a soluble regulatory domain. To date, representative structures for the two parts are available, the transmembrane domain from the structure of a bacterial homologue, the soluble domain from a eukaryotic family member. Despite the strong conservation of the structural framework, the family members show an unusually broad variety of functional behaviors, as some members work as gated chloride channels and others as secondary chloride transporters. The conservation in the structure and the functional resemblance in gating and transport mechanism suggest a strong mechanistic relationship between seemingly contradictory transport modes.

Exploring the gating mechanism in the ClC chloride channel via metadynamics

Computer simulations have been used to probe the gating mechanism in the Salmonella serovar typhimurium chloride channel (st-ClC). Specifically, the recently developed metadynamics methodology has been exploited to construct free energy surfaces as a function of the positions of either one or two chloride ions inside the pore, the position and protonation state of the key E148 residue, and the number of water molecules coordinating the translocating ions. The present calculations confirm the multi-ion mechanism in which an ion-push-ion effect lowers the main barriers to chloride ion translocation. When a second anion is taken into account, the barrier for chloride passage through the E148 narrow region is computed to be 6 kcal/mol in the wild-type channel, irrespective of the protonation state of the E148 residue, which is shown to only affect the entrance barrier. In the E148A mutant, this barrier is much lower, amounting to 3 kcal/mol. The metadynamics calculations reported herein also demonstrate that before reaching the periplasmic solution, chloride ions have to overcome an additional barrier arising from two different effects, namely the rearrangement of their solvation shell and a flip in the backbone angles of the residues E148 and G149, which reside at the end of the alphaF helix.

ClC chloride channels viewed through a transporter lens

Since its discovery, the ClC family of chloride channels has presented biophysicists with unexpected behaviours and unusual surprises. The latest of these is the realization that not only does the family feature genuine chloride channels, it also includes proton-coupled chloride transporters, which move chloride ions and protons across the membrane in opposite directions. The crystal structure of such a transporter serves as a useful platform for understanding ClC channels, and features of chloride/proton exchange-transport may provide a key for comprehending voltage-dependent gating of the channels.

Electrostatics in the cytoplasmic pore produce intrinsic inward rectification in kir2.1 channels

Inward rectifier K+ channels are important in regulating membrane excitability in many cell types. The physiological functions of these channels are related to their unique inward rectification, which has been attributed to voltage-dependent block. Here, we show that inward rectification can also be induced by neutral and positively charged residues at site 224 in the internal vestibule of tetrameric Kir2.1 channels. The order of extent of inward rectification is E224K mutant > E224G mutant > wild type in the absence of internal blockers. Mutating the glycines at the equivalent sites to lysines also rendered weak inward rectifier Kir1.1 channels more inwardly rectifying. Also, conjugating positively charged methanethiosulfonate to the cysteines at site 224 induced strong inward rectification, whereas negatively charged methanethiosulfonate alleviated inward rectification in the E224C mutant. These results suggest that charges at site 224 may control inward rectification in the Kir2.1 channel. In a D172N mutant, spermine interacting with E224 and E299 induced channel inhibition during depolarization but did not occlude the pore, further suggesting that a mechanism other than channel block is involved in the inward rectification of the Kir2.1 channel. In this and our previous studies we showed that the M2 bundle crossing and selectivity filter were not involved in the inward rectification induced by spermine interacting with E224 and E299. We propose that neutral and positively charged residues at site 224 increase a local energy barrier, which reduces K+ efflux more than K+ influx, thereby producing inward rectification.

CLC chloride channels and transporters

CLC proteins are found in cells from prokaryotes to mammals and perform functions in plasma membranes and intracellular vesicles. Several genetic human diseases and mouse models underscore their broad physiological functions in mammals. These functions range from the control of excitability to transepithelial transport, endocytotic trafficking and acidification of synaptic vesicles. The recent crystallization of bacterial CLC proteins gave surprising insights into CLC Cl(-)-channel permeation and gating and provides an excellent basis for structure-function studies. Surprisingly, the CLC from Escherichia coli functions as a Cl-/H+ exchanger, thus demonstrating the thin line separating transporters and channels.

Semimicroscopic modeling of permeation energetics in ion channels

The semimicroscopic (SMC) approach to modeling the energetics of ion permeation through biological channels provides an alternative perspective to standard molecular dynamics methods. It exploits the timescale separation between electronic and structural contributions to dielectric stabilization and accounts for electronic polarization by embedding the channel in a milieu that, on average, describes this polarization. Ions, water, and selected peptide moieties are mobile and comprise the reorganizational contribution to dielectric stabilization. The conceptual advantages and limitations of the technique are described. Methodological details are outlined, stressing three convenient electrical geometries. Practical aspects of the SMC procedure are explained, highlighting the areas ripe for further development. Finally, some specific applications are considered.

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 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.

Structure and function of clc channels

The CLC family comprises a group of integral membrane proteins whose major action is to translocate chloride (Cl-) ions across the cell membranes. Recently, the structures of CLC orthologues from two bacterial species, Salmonella typhimurium and Escherichia coli, were solved, providing the first framework for understanding the operating mechanisms of these molecules. However, most of the previous mechanistic understanding of CLC channels came from electrophysiological studies of a branch of the channel family, the muscle-type CLC channels in vertebrate species. These vertebrate CLC channels were predicted to contain two identical but independent pores, and this hypothesis was confirmed by the solved bacterial CLC structures. The opening and closing of the vertebrate CLC channels are also known to couple to the permeant ions via their binding sites in the ion-permeation pathway. The bacterial CLC structures can probably serve as a structural model to explain the gating-permeation coupling mechanism. However, the CLC-ec1 protein in E. coli was most recently shown to be a Cl- -H+ antiporter, but not an ion channel. The molecular basis to explain the difference between vertebrate and bacterial CLCs, especially the distinction between an ion channel and a transporter, remains a challenge in the structure/function studies for the CLC family.

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.

Unique Structure and Function of Chloride Transporting CLC Proteins

CLC proteins are a large structurally defined family of Cl- ion channels and H+/Cl- antiporters with nine distinct genes in mammals. The membrane-embedded part of CLC proteins bears no obvious similarity to any other class of membrane proteins, while the cytoplasmic C-terminus of most eukaryotic and some prokaryotic CLCs contains two regions with homology to cystathionine beta synthase (CBS) domains that are found in other proteins as well. Different members serve a broad range of physiological roles, including stabilization of the membrane potential, transepithelial ion transport, and vesicular acidification. Their physiological importance is underscored by the causative involvement in at least four different human genetic diseases. From functional studies of the Torpedo homologue ClC-0, a homodimeric architecture with two physically separate ion conduction pathways was anticipated and fully confirmed by solving the crystal structure of prokaryotic CLC homologues. The structure revealed a complex fold of 18 alpha-helices per subunit with at least two Cl- ions bound in the center of each protopore. A critical glutamic acid residue was identified whose side-chain seems to occupy a third Cl- ion binding site in the closed state and that moves away to allow Cl- binding. While the overall architecture and pore structure is certainly conserved from bacteria to humans, the bacterial proteins that were crystallized are actually not Cl- ion channels, but coupled H+/Cl- antiporters. These recent breakthroughs will allow us to study in further detail the structure, function, and the physiological and pathophysiological role of CLC proteins.

Brownian Dynamics Simulation for Modeling Ion Permeation Across Bionanotubes

The principles underlying Brownian dynamics (BD), its statistical consistency, and algorithms for practical implementation are outlined here. The ability to compute current flow across ion channels confers a distinct advantage to BD simulations compared to other simulation techniques. Thus, two obvious applications of BD ion channels are in calculation of the current-voltage and current-concentration curves, which can be directly compared to the physiological measurements to assess the reliability of the model and predictive power of the method. We illustrate how BD simulations are used to unravel the permeation dynamics in two biological ion channels-the KcsA K+ channel and CIC Cl- channel.

Ionic permeation free energy in gramicidin: a semimicroscopic perspective

Ion permeation through the gramicidin channel is studied using a model that circumvents two major difficulties inherent to standard simulational methods. It exploits the timescale separation between electronic and structural contributions to dielectric stabilization, accounting for the influence of electronic polarization by embedding the channel in a dielectric milieu that describes this polarization in a mean sense. The explicit mobile moieties are the ion, multipolar waters, and the carbonyls and amides of the peptide backbone. The model treats the influence of aromatic residues and the membrane dipole potential. A new electrical geometry is introduced that treats long-range electrostatics exactly and avoids problems related to periodic boundary conditions. It permits the translocating ion to make a seamless transition from nearby electrolyte to the channel interior. Other degrees of freedom (more distant bulk electrolyte and nonpolar lipid) are treated as dielectric continua. Reasonable permeation free energy profiles are obtained for potassium, rubidium, and cesium; binding wells are shallow and the central barrier is small. Estimated cationic single-channel conductances are smaller than experiment, but only by factors between 2 (rubidium) and 50 (potassium). When applied to chloride the internal barrier is large, with a corresponding miniscule single-channel conductance. The estimated relative single-channel conductances of gramicidin A, B, and C agree well with experiment.

Computer simulations of membrane proteins

Computer simulations are rapidly becoming a standard tool to study the structure and dynamics of lipids and membrane proteins. Increasing computer capacity allows unbiased simulations of lipid and membrane-active peptides. With the increasing number of high-resolution structures of membrane proteins, which also enables homology modelling of more structures, a wide range of membrane proteins can now be simulated over time spans that capture essential biological processes. Longer time scales are accessible by special computational methods. We review recent progress in simulations of membrane proteins.

A two-holed story: structural secrets about ClC proteins become unraveled?

ClC Cl(-) channels are found in almost all organisms, ranging from bacteria to mammals, in which nine Cl(-) channels belonging to the ClC family have been identified. The biophysical properties and physiological functions of ClC Cl(-) channels have been extensively reviewed. In this short review, we will focus on recent results obtained on the X-ray structure and functional properties of the prokaryotic ClC-ec1 protein and some results obtained on the role of the cytoplasmic COOH terminus of mammalian ClCs.