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

Fluoride Transport and Inhibition Across CLC Transporters

The Chloride Channel (CLC) family includes proton-coupled chloride and fluoride transporters. Despite their similar protein architecture, the former exchange two chloride ions for each proton and are inhibited by fluoride, whereas the latter efficiently transport one fluoride in exchange for one proton. The combination of structural, mutagenesis, and functional experiments with molecular simulations has pinpointed several amino acid changes in the permeation pathway that capitalize on the different chemical properties of chloride and fluoride to fine-tune protein function. Here we summarize recent findings on fluoride inhibition and transport in the two prototypical members of the CLC family, the chloride/proton transporter from Escherichia coli (CLC-ec1) and the fluoride/proton transporter from Enterococcus casseliflavus (CLCF-eca).

The role of conformational change and key glutamic acid residues in the ClC-ec1 antiporter

The triple glutamine (Q) mutant (QQQ) structure of a Cl-/H+ antiporter from Escherichia coli (ClC-ec1) displaying a novel backbone arrangement has been used to challenge the long-held notion that Cl-/H+ antiporters do not operate through large conformational motions. The QQQ mutant substitutes the glutamine residue for an external glutamate E148, an internal glutamate E203, and a third glutamate E113 that hydrogen-bonds with E203. However, it is unknown if QQQ represents a physiologically relevant state, as well as how the protonation of the wild-type glutamates relates to the global dynamics. We herein apply continuous constant-pH molecular dynamics to investigate the H+-coupled dynamics of ClC-ec1. Although any large-scale conformational rearrangement upon acidification would be due to the accumulation of excess charge within the protein, protonation of the glutamates significantly impacts mainly the local structure and dynamics. Despite the fact that the extracellular pore enlarges at acidic pHs, an occluded ClC-ec1 within the active pH range of 3.5-7.5 requires a protonated E148 to facilitate extracellular Cl- release. E203 is also involved in the intracellular H+ transfer as an H+ acceptor. The water wire connection of E148 with the intracellular solution is regulated by the charge states of the E113/E203 dyad with coupled proton titration. However, the dynamics extracted from our simulations are not QQQ-like, indicating that the QQQ mutant does not represent the behavior of the wild-type ClC-ec1. These findings reinforce the necessity of having a protonatable residue at the E203 position in ClC-ec1 and suggest that a higher level of complexity exists for the intracellular H+ transfer in Cl-/H+ antiporters.

A quantitative paradigm for water-assisted proton transport through proteins and other confined spaces

Water-assisted proton transport through confined spaces influences many phenomena in biomolecular and nanomaterial systems. In such cases, the water molecules that fluctuate in the confined pathways provide the environment and the medium for the hydrated excess proton migration via Grotthuss shuttling. However, a definitive collective variable (CV) that accurately couples the hydration and the connectivity of the proton wire with the proton translocation has remained elusive. To address this important challenge-and thus to define a quantitative paradigm for facile proton transport in confined spaces-a CV is derived in this work from graph theory, which is verified to accurately describe water wire formation and breakage coupled to the proton translocation in carbon nanotubes and the Cl^-/H⁺ antiporter protein, ClC-ec1. Significant alterations in the conformations and thermodynamics of water wires are uncovered after introducing an excess proton into them. Large barriers in the proton translocation free-energy profiles are found when water wires are defined to be disconnected according to the new CV, even though the pertinent confined space is still reasonably well hydrated and-by the simple measure of the mere existence of a water structure-the proton transport would have been predicted to be facile via that oversimplified measure. In this paradigm, however, the simple presence of water is not sufficient for inferring proton translocation, since an excess proton itself is able to drive hydration, and additionally, the water molecules themselves must be adequately connected to facilitate any successful proton transport.

Simulation of pH-Dependent Conformational Transitions in Membrane Proteins: The CLC-ec1 Cl-/H+ Antiporter

Intracellular transport of chloride by members of the CLC transporter family involves a coupled exchange between a Cl⁻ anion and a proton (H⁺), which makes the transport function dependent on ambient pH. Transport activity peaks at pH 4.5 and stalls at neutral pH. However, a structure of the WT protein at acidic pH is not available, making it difficult to assess the global conformational rearrangements that support a pH-dependent gating mechanism. To enable modeling of the CLC-ec1 dimer at acidic pH, we have applied molecular dynamics simulations (MD) featuring a new force field modification scheme—termed an Equilibrium constant pH approach (ECpH). The ECpH method utilizes linear interpolation between the force field parameters of protonated and deprotonated states of titratable residues to achieve a representation of pH-dependence in a narrow range of physiological pH values. Simulations of the CLC-ec1 dimer at neutral and acidic pH comparing ECpH-MD to canonical MD, in which the pH-dependent protonation is represented by a binary scheme, substantiates the better agreement of the conformational changes and the final model with experimental data from NMR, cross-link and AFM studies, and reveals structural elements that support the gate-opening at pH 4.5, including the key glutamates Glu_in and Glu_ex.

Divergent Cl- and H+ pathways underlie transport coupling and gating in CLC exchangers and channels

The CLC family comprises H+-coupled exchangers and Cl- channels, and mutations causing their dysfunction lead to genetic disorders. The CLC exchangers, unlike canonical 'ping-pong' antiporters, simultaneously bind and translocate substrates through partially congruent pathways. How ions of opposite charge bypass each other while moving through a shared pathway remains unknown. Here, we use MD simulations, biochemical and electrophysiological measurements to identify two conserved phenylalanine residues that form an aromatic pathway whose dynamic rearrangements enable H+ movement outside the Cl- pore. These residues are important for H+ transport and voltage-dependent gating in the CLC exchangers. The aromatic pathway residues are evolutionarily conserved in CLC channels where their electrostatic properties and conformational flexibility determine gating. We propose that Cl- and H+ move through physically distinct and evolutionarily conserved routes through the CLC channels and transporters and suggest a unifying mechanism that describes the gating mechanism of both CLC subtypes.

Molecular Basis of CLC Antiporter Inhibition by Fluoride

CLC channels and transporters conduct or transport various kinds of anions, with the exception of fluoride, which acts as an effective inhibitor. Here, we performed sub-nanosecond DFT-based QM/MM simulations of the E. coli anion/proton exchanger ClC-ec1 and observed that fluoride binds incoming protons within the selectivity filter, with excess protons shared with the gating glutamate E148. Depending on E148 conformation, the competition for the proton can involve either a direct F^-/E148 interaction or the modulation of water molecules bridging the two anions. The direct interaction locks E148 in a conformation that does not allow for proton transport, and thus inhibits protein function.

Mutation of external glutamate residue reveals a new intermediate transport state and anion binding site in a CLC Cl-/H+ antiporter

The CLC family of proteins are involved in a variety of physiological processes to control cellular chloride concentration. Two distinct classes of CLC proteins, Cl^- channels and Cl^-/H⁺ antiporters, have been functionally and structurally investigated over the last several decades. Previous studies have suggested that the conformational heterogeneity of the critical glutamate residue, Glu_ex, could explain the transport cycle of CLC-type Cl^-/H⁺ antiporters. However, the presence of multiple conformations (Up, Middle, and Down) of the Glu_ex has been suggested from combined structural snapshots of 2 different CLC antiporters: CLC-ec1 from Escherichia coli and cmCLC from a thermophilic red alga, Cyanidioschyzon merolae Thus, we aimed to investigate further the heterogeneity of Glu_ex-conformations in CLC-ec1, the most deeply studied CLC antiporter, at both functional and structural levels. Here, we show that the crystal structures of the Glu_ex mutant E148D and wild-type CLC-ec1 with varying anion concentrations suggest a structural intermediate, the "Midlow" conformation. We also found that an extra anion can be located above the external Cl^--binding site in the E148D mutant when the anion concentration is high. Moreover, we observed that a carboxylate in solution can occupy either the external or central Cl^--binding site in the ungated E148A mutant using an anomalously detectable short carboxylic acid, bromoacetate. These results lend credibility to the idea that the Glu_ex can take at least 3 distinct conformational states during the transport cycle of a single CLC antiporter.

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 Simulations Reveal Key Aspects of the Proton Transport Mechanism in the ClC-ec1 Antiporter

Multiscale reactive molecular dynamics simulations are used to study proton transport through the central region of ClC-ec1, a widely studied ClC transporter that enables the stoichiometric exchange of 2 Cl(-) ions for 1 proton (H(+)). It has long been known that both Cl(-) and proton transport occur through partially congruent pathways, and that their exchange is strictly coupled. However, the nature of this coupling and the mechanism of antiporting remain topics of debate. Here multiscale simulations have been used to characterize proton transport between E203 (Glu(in)) and E148 (Glu(ex)), the internal and external intermediate proton binding sites, respectively. Free energy profiles are presented, explicitly accounting for the binding of Cl(-) along the central pathway, the dynamically coupled hydration changes of the central region, and conformational changes of Glu(in) and Glu(ex). We find that proton transport between Glu(in) and Glu(ex) is possible in both the presence and absence of Cl(-) in the central binding site, although it is facilitated by the anion presence. These results support the notion that the requisite coupling between Cl(-) and proton transport occurs elsewhere (e.g., during proton uptake or release). In addition, proton transport is explored in the E203K mutant, which maintains proton permeation despite the substitution of a basic residue for Glu(in). This collection of calculations provides for the first time, to our knowledge, a detailed picture of the proton transport mechanism in the central region of ClC-ec1 at a molecular level.

Effect of External Electric Field on Substrate Transport of a Secondary Active Transporter

Substrate transport across a membrane accomplished by a secondary active transporter (SAT) is essential to the normal physiological function of living cells. In the present research, a series of all-atom molecular dynamics (MD) simulations under different electric field (EF) strengths was performed to investigate the effect of an external EF on the substrate transport of an SAT. The results show that EF both affects the interaction between substrate and related protein's residues by changing their conformations and tunes the timeline of the transport event, which collectively reduces the height of energy barrier for substrate transport and results in the appearance of two intermediate conformations under the existence of an external EF. Our work spotlights the crucial influence of external EFs on the substrate transport of SATs and could provide a more penetrating understanding of the substrate transport mechanism of SATs.

Molecular Basis for Differential Anion Binding and Proton Coupling in the Cl(-)/H(+) Exchanger ClC-ec1

Cl–/H+ transporters of the CLC superfamily form a ubiquitous class of membrane proteins that catalyze stoichiometrically coupled exchange of Cl– and H+ across biological membranes. CLC transporters exchange H+ for halides and certain polyatomic anions, but exclude cations, F–, and larger physiological anions, such as PO43– and SO42–. Despite comparable transport rates of different anions, the H+ coupling in CLC transporters varies significantly depending on the chemical nature of the transported anion. Although the molecular mechanism of exchange remains unknown, studies on bacterial ClC-ec1 transporter revealed that Cl– binding to the central anion-binding site (Scen) is crucial for the anion-coupled H+ transport. Here, we show that Cl–, F–, NO3–, and SCN– display distinct binding coordinations at the Scen site and are hydrated in different manners. Consistent with the observation of differential bindings, ClC-ec1 exhibits markedly variable ability to support the formation of the transient water wires, which are necessary to support the connection of the two H+ transfer sites (Gluin and Gluex), in the presence of different anions. While continuous water wires are frequently observed in the presence of physiologically transported Cl–, binding of F– or NO3– leads to the formation of pseudo-water-wires that are substantially different from the wires formed with Cl–. Binding of SCN–, however, eliminates the water wires altogether. These findings provide structural details of anion binding in ClC-ec1 and reveal a putative atomic-level mechanism for the decoupling of H+ transport to the transport of anions other than Cl–.

Revealing an outward-facing open conformational state in a CLC Cl(-)/H(+) exchange transporter

CLC secondary active transporters exchange Cl^- for H⁺. Crystal structures have suggested that the conformational change from occluded to outward-facing states is unusually simple, involving only the rotation of a conserved glutamate (Glu_ex) upon its protonation. Using ¹⁹F NMR, we show that as [H⁺] is increased to protonate Glu_ex and enrich the outward-facing state, a residue ~20 Å away from Glu_ex, near the subunit interface, moves from buried to solvent-exposed. Consistent with functional relevance of this motion, constriction via inter-subunit cross-linking reduces transport. Molecular dynamics simulations indicate that the cross-link dampens extracellular gate-opening motions. In support of this model, mutations that decrease steric contact between Helix N (part of the extracellular gate) and Helix P (at the subunit interface) remove the inhibitory effect of the cross-link. Together, these results demonstrate the formation of a previously uncharacterized 'outward-facing open' state, and highlight the relevance of global structural changes in CLC function.

DOI:http://dx.doi.org/10.7554/eLife.11189.001

Cells have transporter proteins on their surface to carry molecules in and out of the cell. For example, the CLC family of transporters move two chloride ions in one direction at the same time as moving one hydrogen ion in the opposite direction.

To be able to move these ions in opposite directions, transporters have to cycle through a series of shapes in which the ions can only access alternate sides of the membrane. First, the transporter adopts an 'outward-facing' shape when the ions first bind to the transporter, then it switches into the 'occluded' shape to move the ions through the membrane. Finally, the transporter takes on the 'inward-facing' shape to release the ions on the other side of the membrane. However, structural studies of CLCs suggest that the structures of these proteins do not change much while they are moving ions, which suggests that they might work in a different way.

Khantwal, Abraham et al. have now used techniques called “nuclear magnetic resonance” and "double electron-electron resonance" to investigate how a CLC from a bacterium moves ions. The experiments suggest that when the transporter adopts the outward-facing shape, points on the protein known as Y419 and D417 shift their positions. Chemically linking two regions of the CLC prevented this movement and inhibited the transport of chloride ions across the membrane.

Khantwal, Abraham et al. then used a computer simulation to model how the protein changes shape in more detail. This model predicts that two regions of the transporter undergo major rearrangements resulting in a gate-opening motion that widens a passage to allow the chloride ions to bind to the protein. Khantwal, Abraham et al.’s findings will prompt future studies to reveal the other shapes and how CLCs transition between them.

DOI:http://dx.doi.org/10.7554/eLife.11189.002

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.

13C NMR detects conformational change in the 100-kD membrane transporter ClC-ec1

CLC transporters catalyze the exchange of Cl(-) for H(+) across cellular membranes. To do so, they must couple Cl(-) and H(+) binding and unbinding to protein conformational change. However, the sole conformational changes distinguished crystallographically are small movements of a glutamate side chain that locally gates the ion-transport pathways. Therefore, our understanding of whether and how global protein dynamics contribute to the exchange mechanism has been severely limited. To overcome the limitations of crystallography, we used solution-state (13)C-methyl NMR with labels on methionine, lysine, and engineered cysteine residues to investigate substrate (H(+)) dependent conformational change outside the restraints of crystallization. We show that methyl labels in several regions report H(+)-dependent spectral changes. We identify one of these regions as Helix R, a helix that extends from the center of the protein, where it forms the part of the inner gate to the Cl(-)-permeation pathway, to the extracellular solution. The H(+)-dependent spectral change does not occur when a label is positioned just beyond Helix R, on the unstructured C-terminus of the protein. Together, the results suggest that H(+) binding is mechanistically coupled to closing of the intracellular access-pathway for Cl(-).

Water access points and hydration pathways in CLC H+/Cl- transporters

CLC transporters catalyze transmembrane exchange of chloride for protons. Although a putative pathway for Cl(-) has been established, the pathway of H(+) translocation remains obscure. Through a highly concerted computational and experimental approach, we characterize microscopic details essential to understanding H(+)-translocation. An extended (0.4 µs) equilibrium molecular dynamics simulation of membrane-embedded, dimeric ClC-ec1, a CLC from Escherichia coli, reveals transient but frequent hydration of the central hydrophobic region by water molecules from the intracellular bulk phase via the interface between the two subunits. We characterize a portal region lined by E202, E203, and A404 as the main gateway for hydration. Supporting this mechanism, site-specific mutagenesis experiments show that ClC-ec1 ion transport rates decrease as the size of the portal residue at position 404 is increased. Beyond the portal, water wires form spontaneously and repeatedly to span the 15-Å hydrophobic region between the two known H(+) transport sites [E148 (Glu(ex)) and E203 (Glu(in))]. Our finding that the formation of these water wires requires the presence of Cl(-) explains the previously mystifying fact that Cl(-) occupancy correlates with the ability to transport protons. To further validate the idea that these water wires are central to the H(+) transport mechanism, we identified I109 as the residue that exhibits the greatest conformational coupling to water wire formation and experimentally tested the effects of mutating this residue. The results, by providing a detailed microscopic view of the dynamics of water wire formation and confirming the involvement of specific protein residues, offer a mechanism for the coupled transport of H(+) and Cl(-) ions in CLC transporters.

Surprises from an unusual CLC homolog

The chloride channel (CLC) family is distinctive in that some members are Cl(-) ion channels and others are Cl(-)/H(+) antiporters. The molecular mechanism that couples H(+) and Cl(-) transport in the antiporters remains unknown. Our characterization of a novel bacterial homolog from Citrobacter koseri, CLC-ck2, has yielded surprising discoveries about the requirements for both Cl(-) and H(+) transport in CLC proteins. First, even though CLC-ck2 lacks conserved amino acids near the Cl(-)-binding sites that are part of the CLC selectivity signature sequence, this protein catalyzes Cl(-) transport, albeit slowly. Ion selectivity in CLC-ck2 is similar to that in CLC-ec1, except that SO(4)(2-) strongly competes with Cl(-) uptake through CLC-ck2 but has no effect on CLC-ec1. Second, and even more surprisingly, CLC-ck2 is a Cl(-)/H(+) antiporter, even though it contains an isoleucine at the Glu(in) position that was previously thought to be a critical part of the H(+) pathway. CLC-ck2 is the first known antiporter that contains a nonpolar residue at this position. Introduction of a glutamate at the Glu(in) site in CLC-ck2 does not increase H(+) flux. Like other CLC antiporters, mutation of the external glutamate gate (Glu(ex)) in CLC-ck2 prevents H(+) flux. Hence, Glu(ex), but not Glu(in), is critical for H(+) permeation in CLC proteins.

Intracellular proton access in a Cl(-)/H(+) antiporter

Mutagenesis, functional analysis, and crystal structures identify a watery tunnel through which protons enter the interior of a Cl⁻/H⁺ antiport protein involved in acid resistance of enteric bacteria.

Chloride-transporting membrane proteins of the CLC family appear in two distinct mechanistic flavors: H⁺-gated Cl⁻ channels and Cl⁻/H⁺ antiporters. Transmembrane H⁺ movement is an essential feature of both types of CLC. X-ray crystal structures of CLC antiporters show the Cl⁻ ion pathway through these proteins, but the H⁺ pathway is known only inferentially by two conserved glutamate residues that act as way-stations for H⁺ in its path through the protein. The extracellular-facing H⁺ transfer glutamate becomes directly exposed to aqueous solution during the transport cycle, but the intracellular glutamate E203, Glu_in, is buried within the protein. Two regions, denoted “polar” and “interfacial,” at the intracellular surface of the bacterial antiporter CLC-ec1 are examined here as possible pathways by which intracellular aqueous protons gain access to Glu_in. Mutations at multiple residues of the polar region have little effect on antiport rates. In contrast, mutation of E202, a conserved glutamate at the protein–water boundary of the interfacial region, leads to severe slowing of the Cl⁻/H⁺ antiport rate. An X-ray crystal structure of E202Y, the most strongly inhibited of these substitutions, shows an aqueous portal leading to Glu_in physically blocked by cross-subunit interactions; moreover, this mutation has only minimal effect on a monomeric CLC variant, which necessarily lacks such interactions. The several lines of experiments presented argue that E202 acts as a water-organizer that creates a proton conduit connecting intracellular solvent with Glu_in.

Chloride-proton antiport proteins of the “CLC” superfamily are transmembrane proteins that form homodimers and are used for myriad physiological purposes, all requiring the coordinated movements of Cl⁻ anions and H⁺ cations in opposite directions across biological membranes. While the pathway for Cl⁻ ions through CLC antiporters is known, we currently have only indirect glimpses of how protons navigate their way through these membrane-embedded proteins. By combining mechanistic and structural approaches, we identify a proton-access pathway in a bacterial Cl⁻/H⁺ antiporter that allows intracellular protons to enter the protein interior and engage in the coupled antiport mechanism. We conclude that E202, a highly conserved glutamate residue, serves to organize water molecules and guide protons to the adjacent glutamate E203 (known as “Glu_in”), a critical residue for the antiport mechanism.

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.

Molecular dynamics investigation of Cl- and water transport through a eukaryotic CLC transporter

Early crystal structures of prokaryotic CLC proteins identified three Cl(-) binding sites: internal (S(int)), central (S(cen)), and external (S(ext)). A conserved external GLU (GLU(ex)) residue acts as a gate competing for S(ext). Recently, the first crystal structure of a eukaryotic transporter, CmCLC, revealed that in this transporter GLU(ex) competes instead for S(cen). Here, we use molecular dynamics simulations to investigate Cl(-) transport through CmCLC. The gating and Cl(-)/H(+) transport cycle are inferred through comparative molecular dynamics simulations with protonated and deprotonated GLU(ex) in the presence/absence of external potentials. Adaptive biasing force calculations are employed to estimate the potential of mean force profiles associated with transport of a Cl(-) ion from S(ext) to S(int), depending on the Cl(-) occupancy of other sites. Our simulations demonstrate that protonation of GLU(ex) is essential for Cl(-) transport from S(ext) to S(cen). The S(cen) site may be occupied by two Cl(-) ions simultaneously due to a high energy barrier (∼8 Kcal/mol) for a single Cl(-) ion to translocate from S(cen) to S(int). Binding two Cl(-) ions to S(cen) induces a continuous water wire from S(cen) to the extracellular solution through the side chain of the GLU(ex) gate. This may initiate deprotonation of GLU(ex), which then drives the two Cl(-) ions out of S(cen) toward the intracellular side via two putative Cl(-) transport paths. Finally, a conformational cycle is proposed that would account for the exchange stoichiometry.

Charge transport in the ClC-type chloride-proton anti-porter from Escherichia coli

The first chloride transporter identified in the superfamily of ClC chloride channels was from Escherichia coli (EClC) (Accardi, A., and Miller, C. (2004) Nature 427, 803-807). Pathways, energetics, and mechanism of proton and chloride translocation and their coupling are up to now unclear. To bridge the hydrophobic gap of proton transport, we modeled four stable buried waters into both subunits of the WT EClC structure. Together they form a "water wire" connecting Glu-203 with the chloride at the central site, which in turn connects to Glu-148, the hypothetical proton exit site. Assuming the transient production of hydrochloride in the central chloride binding site of EClC, the water wire could establish a transmembrane proton transport pathway starting from Glu-203 all the way downstream onto Glu-148. We demonstrated by electrostatic and quantum chemical computations that protonation of the central chloride is energetically feasible. We characterized all chloride occupancies and protonation states possibly relevant for the proton-chloride transport cycle in EClC and constructed a working model. Accordingly, EClC evolves through states involving up to two excess protons and between one and three chlorides, which was required to fulfill the experimentally observed 2:1 stoichiometry. We show that the Y445F and E203H mutants of EClC can operate similarly, thus explaining why they exhibit almost WT activity levels. The proposed mechanism of coupled chloride-proton transport in EClC is consistent with available experimental data and allows predictions on the importance of specific amino acids, which may be probed by mutation experiments.

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.

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.

Proton block of the CLC-5 Cl-/H+ exchanger

CLC-5 is a H(+)/Cl(-) exchanger that is expressed primarily in endosomes but can traffic to the plasma membrane in overexpression systems. Mutations altering the expression or function of CLC-5 lead to Dent's disease. Currents mediated by this transporter show extreme outward rectification and are inhibited by acidic extracellular pH. The mechanistic origins of both phenomena are currently not well understood. It has been proposed that rectification arises from the voltage dependence of a H(+) transport step, and that inhibition of CLC-5 currents by low extracellular pH is a result of a reduction in the driving force for exchange caused by a pH gradient. We show here that the pH dependence of CLC-5 currents arises from H(+) binding to a single site located halfway through the transmembrane electric field and driving the transport cycle in a less permissive direction, rather than a reduction in the driving force. We propose that protons bind to the extracellular gating glutamate E211 in CLC-5. It has been shown that CLC-5 becomes severely uncoupled when SCN(-) is the main charge carrier: H(+) transport is drastically reduced while the rate of anion movement is increased. We found that in these conditions, rectification and pH dependence are unaltered. This implies that H(+) translocation is not the main cause of rectification. We propose a simple transport cycle model that qualitatively accounts for these findings.

CLC channels and transporters: proteins with borderline personalities

Controlled chloride movement across membranes is essential for a variety of physiological processes ranging from salt homeostasis in the kidneys to acidification of cellular compartments. The CLC family is formed by two, not so distinct, sub-classes of membrane transport proteins: Cl(-) channels and H(+)/Cl(-) exchangers. All CLC's are homodimers with each monomer forming an individual Cl- permeation pathway which appears to be largely unaltered in the two CLC sub-classes. Key residues for ion binding and selectivity are also highly conserved. Most CLC's have large cytosolic carboxy-terminal domains containing two cystathionine beta-synthetase (CBS) domains. The C-termini are critical regulators of protein trafficking and directly modulate Cl- by binding intracellular ATP, H+ or oxidizing compounds. This review focuses on the recent mechanistic insights on the how the structural similarities between CLC channels and transporters translate in unexpected mechanistic analogies between these two sub-classes.

Exploring Multidimensional Free Energy Landscapes Using Time-Dependent Biases on Collective Variables

A new implementation of the adaptive biasing force (ABF) method is described. This implementation supports a wide range of collective variables and can be applied to the computation of multidimensional energy profiles. It is provided to the community as part of a code that implements several analogous methods, including metadynamics. ABF and metadynamics have not previously been tested side by side on identical systems. Here, numerical tests are carried out on processes including conformational changes in model peptides and translocation of a halide ion across a lipid membrane through a peptide nanotube. On the basis of these examples, we discuss similarities and differences between the ABF and metadynamics schemes. Both approaches provide enhanced sampling and free energy profiles in quantitative agreement with each other in different applications. The method of choice depends on the dimension of the reaction coordinate space, the height of the barriers, and the relaxation times of degrees of freedom in the orthogonal space, which are not explicitly described by the chosen collective variables.

Permeation of membranes by ribose and its diastereomers

It was recently found that ribose permeates membranes an order of magnitude faster than its diastereomers arabinose and xylose (Sacerdote, M. G.; Szostak, J. W. Proc. Natl. Acad. Sci. U.S.A. 2005, 102, 6004). On this basis it was hypothesized that differences in membrane permeability to aldopentoses provide a mechanism for preferential delivery of ribose to primitive cells for subsequent selective incorporation into nucleotides and their polymers. However, the origins of these unusually large differences have not been well understood. We address this issue in molecular dynamics simulations combined with free energy calculations. It is found that the free energy of transferring ribose from water to the bilayer is lower by 1.5-2 kcal/mol than the barrier for transferring the other two aldopentoses. The calculated and measured permeability coefficients are in excellent agreement. The sugar structures that permeate the membrane are beta-pyranoses, with a possible contribution of the alpha-anomer for arabinose. The furanoid form of ribose is not substantially involved in permeation, even though it is non-negligibly populated in aqueous solution. The differences in free energy of transfer between ribose and arabinose or xylose are attributed, at least in part, to stronger highly cooperative, intramolecular interactions between consecutive exocyclic hydroxyl groups, which are stable in nonpolar media but rare in water. Water/hexadecane partition coefficients of the sugars obtained from separate molecular dynamics simulations correlate with the calculated permeability coefficients, in qualitative agreement with the Overton rule. The relevance of our calculations to understanding the origins of life is discussed.

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.

Redox signaling across cell membranes

Generation of reactive oxygen species (ROS) by plasma membrane-localized NADPH oxidase (Nox 2) is a major mechanism of cell signaling associated with activation of the enzyme by a variety of agonists. With activation, the integral membrane flavocytochrome of Nox 2 transfers an electron from intracellular NADPH to extracellular O(2), generating superoxide anion (O(2)(*-)). The latter dismutes to H(2)O(2) which can diffuse through aquaporin channels in the plasma membrane to elicit an intracellular signaling response. O(2)(*-) also can initiate intracellular signaling by penetration of the cell membrane through anion channels (Cl(-) channel-3, ClC-3). Endosomes containing Nox2 and ClC-3 (called signaling endosomes) are composed of internalized plasma membrane and generate O(2)(*-) in the endosomal lumen to initiate signaling at intracellular sites. Thus, cellular signaling by Nox2 is dependent on the transmembrane flux of ROS. The role of this pathway has only recently been described and will require additional investigation to appreciate its physiological significance fully.

Intracellular proton regulation of ClC-0

Some CLC proteins function as passive Cl⁻ ion channels whereas others are secondary active chloride/proton antiporters. Voltage-dependent gating of the model Torpedo channel ClC-0 is modulated by intracellular and extracellular pH, possibly reflecting a mechanistic relationship with the chloride/proton coupling of CLC antiporters. We used inside-out patch clamp measurements and mutagenesis to explore the dependence of the fast gating mechanism of ClC-0 on intracellular pH and to identify the putative intracellular proton acceptor(s). Among the tested residues (S123, K129, R133, K149, E166, F214L, S224, E226, V227, C229, R305, R312, C415, H472, F418, V419, P420, and Y512) only mutants of E166, F214, and F418 qualitatively changed the pH_int dependence. No tested amino acid emerged as a valid candidate for being a pH sensor. A detailed kinetic analysis of the dependence of fast gate relaxations on pH_int and [Cl⁻]_int provided quantitative constraints on possible mechanistic models of gating. In one particular model, a proton is generated by the dissociation of a water molecule in an intrapore chloride ion binding site. The proton is delivered to the side chain of E166 leading to the opening of the channel, while the hydroxyl ion is stabilized in the internal/central anion binding site. Deuterium isotope effects confirm that proton transfer is rate limiting for fast gate opening and that channel closure depends mostly on the concentration of OH⁻ ions. The gating model is in natural agreement with the finding that only the closing rate constant, but not the opening rate constant, depends on pH_int and [Cl⁻]_int.

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.