Abstract
Mechanisms of anion permeation within ion channels and nanopores remain poorly understood. Recent cryo-electron microscopy structures of the human bestrophin 1 Cl(-) channel (hBest1) provide an opportunity to evaluate ion interactions predicted by molecular dynamics (MD) simulations against experimental observations. Here, we implement the fully polarizable force field AMOEBA in MD simulations on different conformations of hBest1. This force field models multipole moments up to the quadrupole. Using this approach, we model key biophysical properties of the channel that can only be simulated when electronic polarization is included in the molecular models and show that Cl(-) permeation through the neck of the pore is achieved through hydrophobic solvation concomitant with partial ion dehydration. Furthermore, we demonstrate how such polarizable simulations can help determine the identity of ion-like densities within high-resolution cryo-EM structures and demonstrate that neglecting polarization places Cl(-) at positions that do not correspond to their experimentally resolved location. Overall, our results demonstrate the importance of including electronic polarization in realistic and physically accurate models of biological systems, especially channels and pores that selectively permeate anions.