Charge Transfer and Polarization for Chloride Ions Bound in ClC Transport Proteins: Natural Bond Orbital and Energy Decomposition Analyses

Fluoride Ion Binding and Translocation in the CLCF Fluoride/Proton Antiporter: Molecular Insights from Combined Quantum-Mechanical/Molecular-Mechanical Modeling

CLCF fluoride/proton antiporters move fluoride ions out of bacterial cells, leading to fluoride resistance in these bacteria. However, many details about their operating mechanisms remain unclear. Here, we report a combined quantum-mechanical/molecular-mechanical (QM/MM) study of a CLCF homologue from Enterococci casseliflavus (Eca), in accord with the previously proposed windmill mechanism. Our multiscale modeling sheds light on two critical steps in the transport cycle: (i) the external gating residue E118 pushing a fluoride in the external binding site into the extracellular vestibule and (ii) an incoming fluoride reconquering the external binding site by forcing out E118. Both steps feature competitions for the external binding site between the negatively charged carboxylate of E118 and the fluoride. Remarkably, the displaced E118 by fluoride accepts a proton from the nearby R117, initiating the next transport cycle. We also demonstrate the importance of accurate quantum descriptions of fluoride solvation. Our results provide clues to the mysterious E318 residue near the central binding site, suggesting that the transport activities are unlikely to be disrupted by the glutamate interacting with a well-solvated fluoride at the central binding site. This differs significantly from the structurally similar CLC chloride/proton antiporters, where a fluoride trapped deep in the hydrophobic pore causes the transporter to be locked down. A free-energy barrier of 10-15 kcal/mol was estimated via umbrella sampling for a fluoride ion traveling through the pore to repopulate the external binding site.

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

Anion Pathways in the NarK Nitrate/Nitrite Exchanger

NarK nitrate/nitrite antiporter imports nitrate (a mineral form of the essential element nitrogen) into the cell and exports nitrite (a metabolite that can be toxic in high concentrations) out of the cell. However, many details about its operational mechanism remain poorly understood. In this work, we performed steered molecular dynamics simulations of anion translocations and quantum-chemistry model calculations of the binding sites to study the wild-type NarK protein and its R89K mutant. Our results shed light on the importance of the two strictly conserved binding-site arginine residues (R89 and R305) and two glycine-rich signature motifs (G164-M176 and G408-F419) in anion movement through the pore. We also observe conformational changes of the protein during anion migration. For the R89K mutant, our quantum calculations reveal a competition for a proton between the anion (especially nitrite) and lysine, which can potentially slow down or even trap the anion in the pore. Our findings provide a possible explanation for the striking experimental finding that the arginine-to-lysine mutation, despite preserving the charge, impedes or abolishes anion transport in such mutants of NarK and other similar nitrate/nitrite exchangers.

Molecular Simulations of Aqueous Electrolytes: Role of Explicit Inclusion of Charge Transfer into Force Fields

We describe here simulations of aqueous salt solutions that are performed using an explicit charge transfer force field. The emphasis of the discussion is on the calculation of a dynamical property of the solutions: self-diffusion of water. While force fields that are based on pairwise additive potentials or on potentials with explicit inclusion of polarization or with scaled charges can provide at best a qualitative agreement with experiments, force fields with explicit inclusion of charge transfer can produce quantitative agreement with experiment for NaCl and KCl solutions. We argue that a force field with explicit charge transfer contains new physics absent in the previously used force fields described in recent reviews of molecular simulations of aqueous electrolytes.

Comparison of the DFT-SAPT and Canonical EDA Schemes for the Energy Decomposition of Various Types of Noncovalent Interactions

Interaction energies computed with density functional theory can be divided into physically meaningful components by symmetry-adapted perturbation theory (DFT-SAPT) or the canonical energy decomposition analysis (EDA). In this work, the decomposition results obtained by these schemes were compared for more than 200 hydrogen-, halogen-, and pnicogen-bonded, dispersion-bound, and mixed complexes to investigate their similarity in the evaluation of the nature of noncovalent interactions. BLYP functional with D3(BJ) correction was used for the EDA scheme, whereas asymptotically corrected PBE0 functional for DFT-SAPT provided some of the best combinations for description of noncovalent interactions. Both schemes provide similar results concerning total interaction energies and insight into the individual energy components. For most complexes, the dominant energetic term was identified equally by both decomposition schemes. Because the canonical EDA is computationally less demanding than the DFT-SAPT, the former can be especially used in cases where the systems investigated are very large.

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.

Energy decomposition analysis approaches and their evaluation on prototypical protein–drug interaction patterns

The partitioning of the interaction energy into chemical components such as electrostatics, polarization, and charge transfer is possible with energy decomposition analysis approaches. We review and evaluate these for biomolecular applications.