Bases of Bacterial Sodium Channel Selectivity Among Organic Cations

Hille’s (1971) seminal study of organic cation selectivity of eukaryotic voltage-gated sodium channels showed a sharp size cut-off for ion permeation, such that no ion possessing a methyl group was permeant. Using the prokaryotic channel, NaChBac, we found some similarity and two peculiar differences in the selectivity profiles for small polyatomic cations. First, we identified a diverse group of minimally permeant cations for wildtype NaChBac, ranging in sizes from ammonium to guanidinium and tetramethylammonium; and second, for both ammonium and hydrazinium, the charge-conserving selectivity filter mutation (E191D) yielded substantial increases in relative permeability (P_X/P_Na). The relative permeabilities varied inversely with relative K_d calculated from 1D Potential of Mean Force profiles (PMFs) for the single cations traversing the channel. Several of the cations bound more strongly than Na⁺, and hence appear to act as blockers, as well as charge carriers. Consistent with experimental observations, the E191D mutation had little impact on Na⁺ binding to the selectivity filter, but disrupted the binding of ammonium and hydrazinium, consequently facilitating ion permeation across the NaChBac-like filter. We concluded that for prokaryotic sodium channels, a fine balance among filter size, binding affinity, occupancy, and flexibility seems to contribute to observed functional differences.

Elucidating factors important for monovalent cation selectivity in enzymes: E. coli β-galactosidase as a model

Many enzymes require a specific monovalent cation (M(+)), that is either Na(+) or K(+), for optimal activity. While high selectivity M(+) sites in transport proteins have been extensively studied, enzyme M(+) binding sites generally have lower selectivity and are less characterized. Here we study the M(+) binding site of the model enzyme E. coli β-galactosidase, which is about 10 fold selective for Na(+) over K(+). Combining data from X-ray crystallography and computational models, we find the electrostatic environment predominates in defining the Na(+) selectivity. In this lower selectivity site rather subtle influences on the electrostatic environment become significant, including the induced polarization effects of the M(+) on the coordinating ligands and the effect of second coordination shell residues on the charge distribution of the primary ligands. This work expands the knowledge of ion selectivity in proteins to denote novel mechanisms important for the selectivity of M(+) sites in enzymes.

Computational Studies of Glutamate Transporters

Glutamate is the major excitatory neurotransmitter in the human brain whose binding to receptors on neurons excites them while excess glutamate are removed from synapses via transporter proteins. Determination of the crystal structures of bacterial aspartate transporters has paved the way for computational investigation of their function and dynamics at the molecular level. Here, we review molecular dynamics and free energy calculation methods used in these computational studies and discuss the recent applications to glutamate transporters. The focus of the review is on the insights gained on the transport mechanism through computational methods, which otherwise is not directly accessible by experimental probes. Recent efforts to model the mammalian glutamate and other amino acid transporters, whose crystal structures have not been solved yet, are included in the review.

Molecular dynamics simulations of the mammalian glutamate transporter EAAT3

Excitatory amino acid transporters (EAATs) are membrane proteins that enable sodium-coupled uptake of glutamate and other amino acids into neurons. Crystal structures of the archaeal homolog GltPh have been recently determined both in the inward- and outward-facing conformations. Here we construct homology models for the mammalian glutamate transporter EAAT3 in both conformations and perform molecular dynamics simulations to investigate its similarities and differences from GltPh. In particular, we study the coordination of the different ligands, the gating mechanism and the location of the proton and potassium binding sites in EAAT3. We show that the protonation of the E374 residue is essential for binding of glutamate to EAAT3, otherwise glutamate becomes unstable in the binding site. The gating mechanism in the inward-facing state of EAAT3 is found to be different from that of GltPh, which is traced to the relocation of an arginine residue from the HP1 segment in GltPh to the TM8 segment in EAAT3. Finally, we perform free energy calculations to locate the potassium binding site in EAAT3, and find a high-affinity site that overlaps with the Na1 and Na3 sites in GltPh.

Molecular strategies to achieve selective conductance in NaK channel variants

A recent crystallization of several ion channels has provided strong impetus for efforts aimed at understanding the different strategies employed by nature for selective ion transport. In this work, we used two variants of the selectivity filter of NaK channel to explore molecular mechanisms that give rise to K(+)-selectivity. We computed one-dimensional (1D) and two-dimensional (2D) potentials of mean force (PMFs) for ion permeation across the channel. The results indicate that the energies for Na(+) and K(+) permeation across the selectivity filter display significant differences in positions of the binding sites and barriers. One characteristic signature of a K(+)-selective channel is the apparent preservation of the site analogous to that of S2 in KcsA. The S2-bound ion can be almost ideally dehydrated and coordinated by 6 to 8 carbonyls. In a striking contrast, the PMFs controlling transport of ions in a nonselective variant show almost identical profiles for either K(+) or Na(+) and significant involvement of water molecules in ion coordination across the entire selectivity filter. An analysis of differences in 1D PMFs for Na(+) and K(+) suggests that coordination number alone is an insufficient predictor of site selectivity, while chemical composition (ratio of carbonyls and water molecules) correlates well with preference for K(+). Multi-ion effects such as dependence of the barriers and wells for permeant ion on the type of copermeant ion were found to play a significant role in the selectivity signature of the channel as well.

Relative Free Energies for Hydration of Monovalent Ions from QM and QM/MM Simulations

Methods directly evaluating the hydration structure and thermodynamics of physiologically relevant cations (Na(+), K(+), Cl(-), etc.) have wide ranging applications in the fields of inorganic, physical, and biological chemistry. All-atom simulations based on accurate potential energy surfaces appear to offer a viable option for assessing the chemistry of ion solvation. Although MD and free energy simulations of ion solvation with classical force fields have proven their usefulness, a number of challenges still remain. One of them is the difficulty of force field benchmarking and validation against structural and thermodynamic data obtained for a condensed phase. Hybrid quantum mechanical/molecular mechanical (QM/MM) models combined with sampling algorithms have the potential to provide an accurate solvation model and to incorporate the effects from the surrounding, which is often missing in gas-phase ab initio computations. Herein, we report the results from QM/MM free energy simulations of Na(+)/K(+) and Cl(-)/Br(-) hydration where we simultaneously characterized the relative thermodynamics of ion solvation and changes in the solvation structure. The Flexible Inner Region Ensemble Separator (FIRES) method was used to impose a spatial separation between QM region and the outer sphere of solvent molecules treated with the CHARMM27 force field. FEP calculations based on QM/MM simulations utilizing the CHARMM/deMon2k interface were performed with different basis set combinations for K(+)/Na(+) and Cl(-)/Br(-) perturbations to establish the dependence of the computed free energies on the basis set level. The dependence of the computed relative free energies on the size of the QM and MM regions is discussed. The current methodology offers an accurate description of structural and thermodynamic aspects of the hydration of alkali and halide ions in neat solvents and can be used to obtain thermodynamic data on ion solvation in condensed phase along with underlying structural properties of the ion-solvent system.