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

Computational Tools Unravel Putative Sterol Binding Sites in the Lysosomal NPC1 Protein

Two proteins have been linked as the critical components in the molecular mechanisms involved in the Niemann Pick type C (NPC) disease: NPC1, a 140 kDa polytopic membrane-bound protein, and the smaller (132 residues), water-soluble NPC2 protein. NPC1 is believed to act in tandem with NPC2, transferring cholesterol and other sterols out of the LE/Lys compartments. Mutations in either NPC1 or NPC2 can lead to an accumulation of cholesterol and lipids in the LE/Lys, the primary phenotype of the NPC disease, but approximately 95% of identified disease-causing mutations have been mapped to the membrane-bound NPC1 protein. Here, we investigate the full length, membrane-bound NPC1 protein computationally using a combination of molecular modeling, docking, and molecular dynamics (MD) simulations. An analysis of titratable amino acid side chains, several buried in protein pockets, reveals several nonstandard protonation states for the low-pH scenario (pH 5) that is realized in the lysosome. Together with the location of these buried amino acids, docking studies have identified putative lipid binding domains that are in close proximity to amino acids that, when mutated, are connected to NPC1 loss-of-function. Using energy analyses and MD simulations, we analyze these domains as potential cholesterol binding sites and propose the possibility of multiple sterol binding pockets enabling the intramolecular transport of sterol molecules to the transmembrane domain.

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.

Two Cl Ions and a Glu Compete for a Helix Cage in the CLC Proton/Cl- Antiporter

The ubiquitously expressed CLC chloride transporters are involved in a great variety of physiological functions. The CLC protein fold is shared by Cl^- channels and 2Cl^-:1H⁺ antiporters. The antiporters pump three charges per cycle across the membrane with two Cl ions moving in the opposite direction of one proton. Multiconformational continuum electrostatics was used to calculate the coupled thermodynamics of the protonation of the extracellular-facing gating Glu (E_x) and Cl^- binding to the external (S_x) and central (S_c) sites in CLC-ec1, the Escherichia coli exchanger. S_x, S_c, and E_x are buried within the protein where the intersection of two helix N-termini creates a region with a strong, localized positive potential for anion binding. Our chemical potential titrations describe the thermodynamic linkage for binding the Cl^- to each site and protons to E_x. We find that the 2Cl^-:1H⁺ binding stoichiometry is a result of Cl^- binding to S_x requiring H⁺ binding to E_x, whereas Cl^- binding to S_c does not lead to proton uptake. When S_x binds a Cl^-, the protonated E_x moves upward, out of the positive helix cage. The increasing E_x proton affinity on binding the first Cl^- reduces the cost of binding the second Cl^- at either S_x or S_c. Despite the repulsion among the anions, the lowest energy states have two anions bound in the helix cage. The state with no Cl^- is not favored electrostatically, but relies on E_x blocking S_x and on the central residues Y445 and S107 blocking S_c.

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.

pK(A) in proteins solving the Poisson-Boltzmann equation with finite elements

Knowledge on pK(A) values is an eminent factor to understand the function of proteins in living systems. We present a novel approach demonstrating that the finite element (FE) method of solving the linearized Poisson-Boltzmann equation (lPBE) can successfully be used to compute pK(A) values in proteins with high accuracy as a possible replacement to finite difference (FD) method. For this purpose, we implemented the software molecular Finite Element Solver (mFES) in the framework of the Karlsberg+ program to compute pK(A) values. This work focuses on a comparison between pK(A) computations obtained with the well-established FD method and with the new developed FE method mFES, solving the lPBE using protein crystal structures without conformational changes. Accurate and coarse model systems are set up with mFES using a similar number of unknowns compared with the FD method. Our FE method delivers results for computations of pK(A) values and interaction energies of titratable groups, which are comparable in accuracy. We introduce different thermodynamic cycles to evaluate pK(A) values and we show for the FE method how different parameters influence the accuracy of computed pK(A) values.

Identification of a proton-chloride antiporter (EriC) by Himar1 transposon mutagenesis in Lactobacillus reuteri and its role in histamine production

The gut microbiome may modulate intestinal immunity by luminal conversion of dietary amino acids to biologically active signals. The model probiotic organism Lactobacillus reuteri ATCC PTA 6475 is indigenous to the human microbiome, and converts the amino acid L-histidine to the biogenic amine, histamine. Histamine suppresses tumor necrosis factor (TNF) production by human myeloid cells and is a product of L-histidine decarboxylation, which is a proton-facilitated reaction. A transposon mutagenesis strategy was developed based on a single-plasmid nisin-inducible Himar1 transposase/transposon delivery system for L. reuteri. A highly conserved proton-chloride antiporter gene (eriC), a gene widely present in the gut microbiome was discovered by Himar1 transposon (Tn)-mutagenesis presented in this study. Genetic inactivation of eriC by transposon insertion and genetic recombineering resulted in reduced ability of L. reuteri to inhibit TNF production by activated human myeloid cells, diminished histamine production by the bacteria and downregulated expression of histidine decarboxylase cluster genes compared to those of WT 6475. EriC belongs to a large family of ion transporters that includes chloride channels and proton-chloride antiporters and may facilitate the availability of protons for the decarboxylation reaction, resulting in histamine production by L. reuteri. This report leverages the tools of bacterial genetics for probiotic gene discovery. The findings highlight the widely conserved nature of ion transporters in bacteria and how ion transporters are coupled with amino acid decarboxylation and contribute to microbiome-mediated immunomodulation.

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.

Understanding selectin counter-receptor binding from electrostatic energy computations and experimental binding studies

Higher organisms defend themselves against invading micro-organisms and harmful substances with their immune system. Key players of the immune system are the white blood cells (WBC), which in case of infection move in an extravasation process from blood vessels toward infected tissue promoting inflammation. This process starts with the attachment of the WBC to the blood vessel wall, mediated by protein pair interactions of selectins and counter-receptors (C-R). Individual selectin C-R binding is weak and varies only moderately between the three selectin types. Multivalency enhances such small differences, rendering selectin-binding type specific. In this work, we study selectin C-R binding, the initial step of extravasation. We performed electrostatic energy computations based on the crystal structure of one selectin type co-crystallized with the ligating part of the C-R. The agreement with measured free energies of binding is satisfactory. Additionally, we modeled selectin mutant structures in order to explain differences in binding of the different selectin types. To verify our modeling procedures, surface plasmon resonance data were measured for several mutants and compared with computed binding affinities. Binding affinities computed with soaked rather than co-crystallized selectin C-R structures do not agree with measured data. Hence, these structures are inappropriate to describe the binding mode. The analysis of selectin/C-R binding unravels the role played by individual molecular components in the binding event. This opens new avenues to prevent immune system malfunction, designing drugs that can control inflammatory processes by moderating selectin C-R binding.

Evolution of the genetic code by incorporation of amino acids that improved or changed protein function

Fifty years have passed since the genetic code was deciphered, but how the genetic code came into being has not been satisfactorily addressed. It is now widely accepted that the earliest genetic code did not encode all 20 amino acids found in the universal genetic code as some amino acids have complex biosynthetic pathways and likely were not available from the environment. Therefore, the genetic code evolved as pathways for synthesis of new amino acids became available. One hypothesis proposes that early in the evolution of the genetic code four amino acids-valine, alanine, aspartic acid, and glycine-were coded by GNC codons (N = any base) with the remaining codons being nonsense codons. The other sixteen amino acids were subsequently added to the genetic code by changing nonsense codons into sense codons for these amino acids. Improvement in protein function is presumed to be the driving force behind the evolution of the code, but how improved function was achieved by adding amino acids has not been examined. Based on an analysis of amino acid function in proteins, an evolutionary mechanism for expansion of the genetic code is described in which individual coded amino acids were replaced by new amino acids that used nonsense codons differing by one base change from the sense codons previously used. The improved or altered protein function afforded by the changes in amino acid function provided the selective advantage underlying the expansion of the genetic code. Analysis of amino acid properties and functions explains why amino acids are found in their respective positions in the genetic code.

Understanding photosynthetic light-harvesting: a bottom up theoretical approach

We discuss a bottom up approach for modeling photosynthetic light-harvesting. Methods are reviewed for a full structure-based parameterization of the Hamiltonian of pigment-protein complexes (PPCs). These parameters comprise (i) the local transition energies of the pigments in their binding sites in the protein, the site energies; (ii) the couplings between optical transitions of the pigments, the excitonic couplings; and (iii) the spectral density characterizing the dynamic modulation of pigment transition energies and excitonic couplings by protein vibrations. Starting with quantum mechanics perturbation theory, we provide a microscopic foundation for the standard PPC Hamiltonian and relate the expressions obtained for its matrix elements to quantities that can be calculated with classical molecular mechanics/electrostatics approaches including the whole PPC in atomic detail and using charge and transition densities obtained with quantum chemical calculations on the isolated building blocks of the PPC. In the second part of this perspective, the Hamiltonian is utilized to describe the quantum dynamics of excitons. Situations are discussed that differ in the relative strength of excitonic and exciton-vibrational coupling. The predictive power of the approaches is demonstrated in application to different PPCs, and challenges for future work are outlined.

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.

The coupled proton transport in the ClC-ec1 Cl(-)/H(+) antiporter

Using a reactive molecular dynamics simulation methodology, the free energy barrier for water-mediated proton transport between the two proton gating residues Glu(203) and Glu(148) in the ClC-ec1 antiporter, including the Grotthuss mechanism of proton hopping, was calculated. Three different chloride-binding states, with 1), both the central and internal Cl(-), 2), the central Cl(-) only, and 3), the internal Cl(-) only, were considered and the coupling to the H(+) transport studied. The results show that both the central and internal Cl(-) are essential for the proton transport from Glu(203) to Glu(148) to have a favorite free energy driving force. The rotation of the Glu(148) side chain was also found to be independent of the internal chloride binding state. These results emphasize the importance of the 2:1 stoichiometry of this well-studied Cl(-)/H(+) antiporter.