Mechanisms of anion permeation within ion channels and nanopores remain poorly understood. Recent cryo-electron microscopy structures of the human bestrophin 1 chloride channel (hBest1) provide an opportunity to evaluate ion interactions predicted by molecular dynamics (MD) simulations against experimental observations. We implement the fully polarizable forcefield AMOEBA in MD simulations of open and partially-open states of the hBest1. The AMOEBA forcefield models multipole moments up to the quadrupole; therefore, it captures induced dipole and anion-π interactions. By including polarization we demonstrate the key role that aromatic residues play in ion permeation and the functional advantages of pore asymmetry within the highly conserved hydrophobic neck of the pore. We establish that these only arise when electronic polarization is included in the molecular models. We also show that Cl⁻ permeation in this region can be achieved through hydrophobic solvation concomitant with partial ion dehydration, which is compensated for by the formation of contacts with the edge of the phenylalanine ring. Furthermore, we demonstrate how polarizable simulations can help determine the identity of ion-like densities within high-resolution cryo-EM structures. Crucially, neglecting polarization in simulation of these systems results in the localization of Cl⁻ at positions that do not correspond with their experimentally resolved location. Overall, our results demonstrate the importance of including electronic polarization in realistic and physically accurate models of biological systems.
Keywords: Bestrophin channels; Molecular dynamics; anion permeation; anion-π; chloride selectivity.