Abstract
The ability of macrocycles, enzymes, ion channels, transporters, and DNA to differentiate among ion types is often crucial to their function. Using molecular dynamics simulations on both detailed systems and simple models, we quantify the importance of several factors which affect the ion selectivity of such molecules, including the number of coordinating ligands, their dipole moment, and their vibrational motion. The information resulting from our model systems is distilled into a series of selectivity maps that can be used to read off the relative free energy associated with binding of different ions, and to provide an estimate of the importance of the various factors. Although our maps cannot capture all elements of real systems, it is remarkable that they produce differential site-binding energies that are in line with experiment and more-detailed simulations for a variety of systems-making them useful for understanding the origins of selective binding and transport. The chemical nature of the coordinating ligands is essential for creating thermodynamic ion selectivity in flexible molecules (such as 18c6), but as the binding site becomes more rigid, the number of ligands (as in ion channels) and the reduction of thermal fluctuations (as in amino-acid transporters) can become important. In the future, our maps could aid in the determination of the local structure from binding energies and assist in the design of novel ion selective molecules.