Next Article in Journal
Synergistic Enhancement of Carrier Migration by SnO2/ZnO@GO Heterojunction for Rapid Degradation of RhB
Previous Article in Journal
Integrating Transcriptome and Chemical Analyses to Provide Insights into Biosynthesis of Terpenoids and Flavonoids in the Medicinal Industrial Crop Andrographis paniculate and Its Antiviral Medicinal Parts
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Molecular Mechanism of Ion Selectivity in Nanopores

Key Laboratory of Orogenic Belts and Crustal Evolution, The School of Earth and Space Sciences, Ministry of Education, Peking University, Beijing 100871, China
*
Author to whom correspondence should be addressed.
Molecules 2024, 29(4), 853; https://doi.org/10.3390/molecules29040853
Submission received: 9 November 2023 / Revised: 3 February 2024 / Accepted: 12 February 2024 / Published: 14 February 2024
(This article belongs to the Section Chemical Biology)

Abstract

:
Ion channels exhibit strong selectivity for specific ions over others under electrochemical potentials, such as KcsA for K+ over Na+. Based on the thermodynamic analysis, this study is focused on exploring the mechanism of ion selectivity in nanopores. It is well known that ions must lose part of their hydration layer to enter the channel. Therefore, the ion selectivity of a channel is due to the rearrangement of water molecules when entering the nanopore, which may be related to the hydrophobic interactions between ions and channels. In our recent works on hydrophobic interactions, with reference to the critical radius of solute (Rc), it was divided into initial and hydrophobic solvation processes. Additionally, the different dissolved behaviors of solutes in water are expected in various processes, such as dispersed and accumulated distributions in water. Correspondingly, as the ion approaches the nanopore, there seems to exist the “repulsive” or “attractive” forces between them. In the initial process (<Rc), the energy barrier related to “repulsive” force may be expected as ions enter the channel. Regarding the ion selectivity of nanopores, this may be due to the energy barrier between the ion and channel, which is closely related to the ion size and pore radius. Additionally, these may be demonstrated by the calculated potential mean forces (PMFs) using molecular dynamics (MD) simulations.

1. Introduction

Ion channels with nanoscale pores are specialized membrane proteins. They play a vital role in allowing specific ion movements across the cell membrane, which is known as ion selectivity and is a fundamental property defining ion channel function. From this, they perform key functions in cell behavior, such as electrical excitability, regulation of cell volume, regulation of cytoplasmic ion concentration and pH, and hormone and neurotransmitter secretion [1]. Additionally, due to analogies with biological membranes, much attention has been paid to the selectively permeable membranes with sub-nm pores because of their potential applications, such as in water filtration, molecular separation, and the desalination of seawater [2,3,4]. Therefore, it is important to investigate the mechanism of ion selectivity in nanopores.
In fact, the current knowledge on such a selective process was laid down by MacKinnon’s [5,6,7] group in the reports on the structure of KcsA, a potassium channel from Streptomyces lividans. Potassium (K+) channels are tetrameric membrane-spanning proteins that provide a selective pore for the conductance of K+ across the cell membranes. This structure revealed that, in the narrow pore region forming the selectivity filter, the backbone carbonyl oxygen atoms (as well as threonine hydroxyls) of the TVGYG sequence form a series of cage-like K+ binding sites, termed S1–S4, and a non-cage-forming S0 site. Under physiological conditions, the channels are K+-selective protein pores in cell membranes, which display up to a 1000-fold preference for K+ over Na+ [8,9], despite the fact that the two ions are both spherical in nature, have the same charge, and differ in atomic radii by only 0.38 Å. Therefore, the selectivity filter is the functional unit that allows K+ channels to distinguish potassium (K+) and sodium (Na+) ions.
Many experimental and theoretical works have been conducted on the ion selectivity of the KcsA channel [10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31]. Additionally, the selectively permeable graphene nanopores have also been attracting much attention to unravel the mechanism of ion selectivity in filters [32,33]. These show that, under transmembrane voltage bias, a nanopore containing four carbonyl groups to mimic the selectivity filter of the KcsA K+ channel preferentially conducts K+ over Na+ [32]. From the free energy calculations, these indicate that, compared to K+, the smaller radius of Na+ leads to a significantly higher free energy barrier in the nanopore of a certain size [32]. Additionally, it indicates that the ion selectivity of such biomimetic graphene nanopores can be simply controlled by the size of the nanopore [33], and the graphene nanopores with a distance of about 3.9 Å between two neighboring oxygen atoms present excellent ion selectivity for Na+ and K+ ions.
So far, various models have been proposed to understand the mechanism of ion selectivity of channels. The classical “snug-fit” mechanism [34] suggested that the permeating ions would need to fit snugly in a rigid selectivity filter, which would allow the coordinating residues to solvate the preferred ions and provide a parsimonious explanation for the exclusion of smaller ions, such as Na+ in potassium channels. In the field strength model of Eisenman [35,36], it is suggested that ion selectivity may be regulated by the strength of the interactions between the ions and the residues. From Noskov et al.’s study [37], as residues are packed as tightly as they are in the selectivity filter, it is necessary to consider not only the ion–residue interactions but also the residue–residue interactions. Therefore, they proposed a “carbonyl-repulsion” mechanism whereby the increased repulsion upon coordinating the smaller Na+ leads to K+/Na+ selectivity in K+ channels. In addition, the “over-coordination hypothesis” is proposed by several groups [38,39,40,41]. This means that each cation is coordinated by more residues than it would experience in bulk water. Therefore, ion selectivity is due to altering the number of residues, not their field strength. From the above, further study is necessary to understand the mechanism of ion selectivity in nanopores.
In fact, as K+ passes through the KcsA filter, the ion must be almost completely dehydrated within the structural confinement of a narrow pore, which undoubtedly leads to the changes in hydrogen bondings of water. The ion selectivity of nanopores may be related to the rearrangement of water molecules when the ion enters the nanopore. In combination with our recent studies [42,43,44,45,46,47,48,49,50,51], this is closely related to hydrophobic interactions, which may be reasonably ascribed to the structural competition between interfacial and bulk water.
Based on thermodynamic analysis and MD simulations, this work is devoted to investigating the molecular mechanism of ion selectivity in nanopores. In our recent studies [47,48,49,50,51], hydration free energy was derived and applied to study the origin of hydrophobic effects. Hydration free energy is dependent on the solute size. With reference to the critical radius of solute (Rc), this is divided into initial and hydrophobic solvation processes, corresponding to the different dissolving behaviors of solutes in water, such as dispersed and accumulated distributions in solution [47,48,49,50,51]. Correspondingly, as the ion approaches the nanopore, “repulsive” or “attractive” forces may be expected between them. In the initial process (<Rc), the energy barrier is expected as ions enter the channel. The ion selectivity in nanopores may be related to the energy barrier when the ion enters the channel. It is found that the energy barrier may be dependent on the ion size and pore radius. Additionally, the barrier may also be modulated by the direct interactions between the ions and nanopores.

2. Results and Discussion

2.1. Thermodynamics of Ion Selectivity

In thermodynamics, Gibbs free energy (ΔG) can be utilized to investigate whether a process is likely to occur. As the solutes are embedded into water, this includes the solute–solute, the solute–water, and the water–water interactions. The total Gibbs energy may be expressed as
Δ G T o t a l = Δ G S o l u t e s o l u t e + Δ G S o l u t e w a t e r + Δ G W a t e r w a t e r
In fact, this may be applied to investigate the ion selectivity of channels as the ions pass through the filters.
For the selectivity filter of the KcsA channel, the ion selectivity represents a competition of K+ and Na+ for a binding site in the filter. The selectivity for K+ over Na+ is due to the difference in free energies of the ions in the pore, which departs from the corresponding difference in the bulk solution. It follows that the ion selectivity of binding is fundamentally governed by differences in relative free energies [31], and the ion selectivity can be stated from the point of view of thermodynamic equilibrium as
Δ Δ G K + N a + = Δ G s i t e K + N a + Δ G b u l k K + N a +
Thermodynamically, it can be understood to allow the preferential permeation of K+ over Na+ by providing an unfavorable environment for Na+ in the presence of K+.
Regarding the KcsA K+ channel, there are mainly O-atom quaternion rings with three sizes: 3.12/3.33 Å, 3.89 Å, and 4.22 Å. When the ions permeate the nanopore channels, one unambiguous observation based on the crystallographic X-ray structures of K+ channels is that the K+ ions must become essentially dehydrated to pass through the narrow selectivity filter [6] and directly interact with the filter. Therefore, as the ions pass through the nanopore channels, this undoubtedly leads to the changes in water structure.
As the ion permeates the nanopore, the ion selectivity of the nanopore may be closely related to not only the changes in hydrogen bondings in water, but also the direct interactions between the ion and channel. Therefore, the ion selectivity of the nanopore may be related to various interactions, such as those due to the solute–solute, solute–water, and water–water interactions. In fact, before the ions are affected by the interactions between ions and the channel, it is necessary for the ions to approach the filter. This may be due to the rearrangement of water molecules as ions approach the nanopore, which is related to ΔGWater-water and ΔGSolute-water. Therefore, it is important to investigate the structure of liquid water, and the effects of dissolved solute on the water structure.
Numerous experimental and theoretical works have been conducted to investigate the structure of liquid water. To date, a variety of structural models for water have been proposed, broadly categorized into two groups: (a) mixture models and (b) distorted hydrogen bond or continuum models [52]. The mixture model postulates the concurrent existence of two distinct types of structures. In the latter model, this implies that water consists of a random, three-dimensional network of hydrogen bonds with a wide distribution of O-H∙∙∙O hydrogen bond angles and distances. However, unlike in the mixture models, the water networks cannot be “broken” or separated into distinct molecular species.
The OH vibrations are sensitive to hydrogen bondings and are employed to study the structure of water. Based on the dependence of OH vibrations on water clusters ((H2O)n) [42,43,44], it is found that, as three-dimensional hydrogen bondings emerge (n ≥ 6), OH vibrations are predominantly influenced by the local hydrogen bondings (within the first shell) of water molecules, and the effects of hydrogen bonding beyond the first shell on OH vibrations are weak. Hence, different OH vibrations may reasonably be attributed to the OH vibrations engaged into various local hydrogen-bonded networks of a water molecule.
Under ambient conditions, the Raman OH stretching bands of water can be fitted into five sub-bands. These may be attributed to OH vibrations involved in various local hydrogen bonds, including DDAA (double donor-double acceptor, tetrahedral hydrogen bonding), DDA (double donor-single acceptor), DAA (single donor-double acceptor), and DA (single donor-single acceptor) hydrogen bondings, respectively. In ambient water, a water molecule interacts with neighboring water molecules (the first shell) through various local hydrogen-bonded networks [44]. It is different from the mixture and continuum structural models of water. Additionally, under ambient conditions, the local hydrogen-bonded networks of a water molecule may be influenced by factors such as temperature, pressure, dissolved salt, and confined environments.
When a solute is embedded in water, an interface appears between the solute and water, which undoubtedly affects the structure of liquid water. As the OH vibration is primarily influenced by the local hydrogen bondings of a water molecule, the dissolved solute mainly affects the structure of interfacial water (the topmost water layer at the solute/water interface). In theory, vibrational sum frequency generation (SFG) spectroscopy is a technique selectively applicable to interfaces and is widely used to examine the structure of the air/water interface. Based on the SFG study [45] on the air/water interface, it is observed that, in contrast to bulk water, no DDAA (tetrahedral) hydrogen bondings are present in interfacial water.
The solute mainly affects the hydrogen bondings of interfacial water. The formation of the solute/water interface is related to the loss of DDAA hydrogen bonding in interfacial water layer. Once the ratio of the interfacial water layer to the volume is established, it can be employed to calculate the Gibbs free energy between solute and water (interfacial water). Therefore, ΔGSolute-water can reasonably be expressed as
Δ G S o l u t e w a t e r = n H B · Δ G D D A A · R I n t e r f a c i a l     w a t e r
where RInterfacial water represents the ratio of interfacial to bulk water, ∆GDDAA is the Gibbs energy of DDAA hydrogen bondings, and nHB means the average hydrogen bondings of a water molecule. For DDAA hydrogen bonidngs, nHB is equal to 2.
The hydration free energy refers to the change in Gibbs energy when an ion (or molecule) is transferred from a vacuum (or the gas phase) to a solvent. Thermodynamically, the process of inserting a hard sphere into water is equivalent to the formation of an empty spherical cavity in water. After the solute is regarded as an ideal sphere, the hydration free energy may be expressed as [47]
Δ G H y d r a t i o n = Δ G W a t e r w a t e r + Δ G S o l u t e w a t e r = Δ G W a t e r w a t e r + 8   ·   Δ G D D A A   ·   r H 2 O R
where rH2O is the average radius of an H2O molecule, and R represents the radius of the solute. At ambient conditions, the Gibbs energy of DDAA hydrogen bondings is determined to be −2.66 kJ/mol [47].
In principle, the lower the hydration free energy, the more thermodynamically stable the system becomes. The hydration free energy is the sum of Gibbs energy of pure water (ΔGWater-water) and that of interfacial water (ΔGSolute-water) related to the dissolved solute. Therefore, the structural transition is expected when ΔGWater-water equals ΔGSolute-water (Figure 1).
Δ G W a t e r w a t e r = Δ G S o l u t e w a t e r             R c = 8 · Δ G D D A A · r H 2 O Δ G W a t e r w a t e r
where Rc represents the critical radius of solute. At 293 K and 0.1 MPa, Rc is calculated to be 6.5 Å. With increasing the solute size (or concentrations), it is divided into the initial and hydrophobic solvation processes [47,48]. Additionally, the Gibbs free energy of interfacial water (ΔGSolute-water) is inversely proportional to the radius of solute (or the ratio of surface area to volume of solute). As the solute size (or solute concentrations) increases, different dissolved behaviors can be expected for various processes (Figure 1) [47,48].
In the initial solvation process (∆GSolute-water < ∆GWater-water, both of which are negative), the hydration free energy is dominated by the Gibbs energy of interfacial water (ΔGSolute-water). To achieve greater thermodynamic stability, the solutes tend to be dispersed in water to maximize ΔGSolute-water. Consequently, water molecules are expected to exist between solutes. In other words, it appears that there are “repulsive” forces between the solutes as they are brought closer together. In the hydrophobic solvation process (∆GSolute-water > ∆GWater-water), the hydration free energy is primarily governed by the Gibbs energy of water (ΔGWater-water). To maximize the Gibbs energy of water, the dissolved solutes are expected to be aggregated to minimize the interfacial water. It seems that there are “attractive” forces between the dissolved solutes.
Hydrophobic effects are generally defined as the tendency of non-polar molecules (or molecular surfaces) to aggregate in an aqueous solution. In our recent investigations [47,48] on hydrophobic interactions, they were reasonably characterized as the tendency to minimize the ratio of the surface area to the volume of the solutes in order to maximize the hydrogen bondings of water. Of course, this is due to the dissolved solute mainly affecting the structure of interfacial water, and the hydrogen bondings of interfacial water are weaker than bulk water. Additionally, the origin of hydrophobic interactions arises from the structural competition between interfacial and bulk water [47,48].
The Gibbs free energy between solute and water (ΔGSolute-water) is closely related to the water molecules of the interfacial water layer. When two solutes are embedded into water, the interfacial water of two solutes may be inversely proportional to the distance between them. Therefore, ΔGSolute-water is reasonably expressed as
Δ G S o l u t e w a t e r 1 / R S o l u t e s o l u t e
where RSolute-solute means the distance between solutes. As two identical spheres are dissolved in ambient water, the Rc is 3.25 Å [48]. As the distance between solutes is larger or less than Rc, hydrophobic interactions between solutes may be regarded as the “attractive” or the “repulsive” forces [48].
To investigate the ion selectivity of nanopores, the test solute is a sphere, and the target filter is the membrane embedded with a nanopore (Figure 2). As the distance between the ion and filter is decreased, the hydration free energy may reasonably be expressed as
Δ G H y d r a t i o n = Δ G W a t e r w a t e r + Δ G I o n f i l t e r d x d y
where ΔGIon-filter is the Gibbs energy of interfacial water arising from the ion and filter, which is inversely proportional to the distance between them, ∝1/dIon-filter, where dIon-filter is the distance between the ion and filter (Figure 2).
In fact, as the nanopore is filled with a circle with corresponding radius R, the filter becomes the membrane with the same size as the filter (Figure 2). Therefore, as the ion approaches the filter, the Gibbs free energy (ΔGIon-filter) related to the ion and filter equals the difference in Gibbs energy between the ion and membrane (ΔGIon-membrane), and that between the ion and filled circle with radius R (ΔGIon-circle) (Figure 2), which may be expressed as
Δ G I o n f i l t e r = Δ G I o n m e m b r a n e Δ G I o n c i r c l e
As the ion enters the nanopore, the interfacial water molecules related to the ion and nanopore are changed into bulk water, which may be related to the ΔGIon-circle.
In comparison with the membrane, the test ion is treated as a point after the size of the ion is ignored. Additionally, the ion is restrained to perpendicularly move to the center of the membrane. Therefore, the Gibbs energy between the ion and membrane (ΔGIon-membrane) may mathematically be determined as (Figure 2)
Δ G I o n m e m b r a n e 1 / d I o n m e m b r a n e = π · w 2 + d 2 d
where w is the radius of membrane, and d means the distance between the centers of the ion and membrane. In addition, ΔGIon-membrane is higher than ΔGWater-water (ΔGIon-membrane >ΔGWater-water), and the hydrophobic interactions between the test ion and membrane substrate can reasonably be regarded as the “attractive” force [48].
Additionally, Equation (6) may also be used to determine the Gibbs free energy between the ion and the circle (ΔGIon-circle). In this work, the nanopore (or circle) is simply regarded as a sphere with the radius RPore. Of course, ΔGIon-circle is also dependent on the distance between them. This may be reasonably expressed as
Δ G I o n c i r c l e 1 / d I o n c i r c l e = 1 / R P o r e + r I o n + s e p
where dIon-circle is the distance between the centers of the ion and the circle (or nanopore), RPore is the pore radius, rIon is the ion size, and sep is the separation between the ion and the circle (Figure 2). Therefore, when the ion approaches the nanopore, the Gibbs energy of interfacial water may be related to the pore radius (RPore), ion size (rIon), and the separation between them. However, different from the ΔGIon-membrane, the hydrophobic interactions between the ion and the circle (ΔGIon-circle) may be the “attractive” or “repulsive” forces, which may be closely related to the ion size and pore radius.
As the ion and filter are embedded into water, hydration free energy is the sum of the Gibbs energy of water and interfacial water, and the latter is composed of the interactions between the ion and membrane substrate (ΔGIon-membrane), and those between the ion and the circle with radius R (ΔGIon-circle). From Equation (10), it can be seen that a decrease in the dIon-circle may increase the ΔGIon-circle. In our recent works on hydrophobic interactions [47,48], dc (or Rc) was expected to be
Δ G I o n m e m b r a n e Δ G I o n c i r c l e     ( d = d c ) = Δ G W a t e r w a t e r C o r r e s p o n d i n g   R P o r e , c   a n d   r I o n , c
From the equation, dc is expected as the ion and channel come into contact (or with the separation between them being zero, sep = 0). Therefore, as an ion enters the channel, the dc may be closely related to the ion size (rIon) and the radius of the channel (RPore). Regarding an ion with a specific size rIon, the corresponding pore radius (RPore,c) may be expected for the ion, RPore,c for rIon,c = rIon, or vice versa. Based on the discussion on hydrophobic interactions, with reference to dc, hydrophobic interactions between the ion and filter may be either “repulsive” or “attractive” forces when the ion approaches the filter.
Regarding the specific ion with radius rIon, as the radius of the pore is larger than RPore,c (RPore > RPore,c and rIon = rIon,c), it is associated with the hydrophobic process. Therefore, as the ion approaches the filter, an “attractive” force may be expected between them, and no energy barrier is expected (Figure 3). In other words, the ion is expected to penetrate the pore of the filter, similar to the ion diffusion in aqueous solutions.
Additionally, as the radius of the nanopore is less than RPore,c (RPore < RPore,c and rIon = rIon,c), it is related to the initial solvation process. In the initial solvation process, the “repulsive” force between solutes is expected as they pushed together in water. In other words, as the ion approaches the nanopore, the energy barrier may be expected between the ion and the filter (Figure 3).
From the above, as the ion approaches the nanopore, the energy barrier between the ion and the channel may be related to the difference in Gibbs energy between ΔGIon-circle and ΔGIon-circle (d = dc). This may be expressed as
Δ G E n e r g y   b a r r i e r = Δ G I o n c i r c l e Δ G I o n c i r c l e     ( d = d c ) 1 / R P o r e + r I o n 1 / R P o r e , c + r I o n , c
Therefore, when an ion approaches the nanopore, the energy barrier between them is related to the ion size (rIon) and pore radius (RPore).
Additionally, regarding the specific pore with the radius being RPore, the energy barrier may also be expected as the ion size is less than rIon,c, RPore = RPore,c and rIon < rIon,c. For the filter with specific pore radius (RPore = RPore,c), it can be derived that the energy barrier may be dependent on the ion size (rIon). In comparison with rIon,c, the lower energy barrier may be expected for the larger ion, and the higher energy barrier is related to the smaller ion (Figure 4). This may be utilized to understand the ion selectivity of the nanopore.
When the ions permeate the KcsA K+ channel, they become essentially dehydrated to pass through the narrow selectivity filter. Therefore, as the ion enters the nanopore, direct interactions, such as van der Waals interactions, are expected between the ion and nanopore. This means that the energy barrier may also be affected by the intermolecular interactions between the ion and nanopore. In other words, the ion selectivity of the nanopore may be modulated through changing the molecular polarity of the channel. In fact, this can be demonstrated with many theoretical simulations [31,32].
As ion size is less than rIon,c (rIon < rIon,c for RPore = RPore,c), an energy barrier is expected between the ion and filter, and the ion may be suppressed to enter the nanopore of the filter. In other words, the external driving force, such as electrical potential or pressure, is necessarily added to the ion to surmount the energy barrier and pass through the nanopore. Of course, this is in accordance with the experimental studies on ion channel of the filter, which means that the electrical transmembrane potential is necessary inside and outside of the membrane to drive the ion through the channel.
From the above, the energy barrier between the ion and nanopore is related to ion size and pore radius. Regarding the nanopore with specific radius RPore (RPore = RPore,c), a lower energy barrier is expected for large ions, and a higher energy barrier is related to small ions. Therefore, compared to small ions, the weaker external force may be enough to make large ions pass through the nanopore. This means that the ion selectivity of nanopores may also be related to the strength of external forces.
When the ion approaches the nanopore, “repulsive” or “attractive” forces may be expected between the ion and nanopore. From the above, dc is related to not only the ion size (rIon) but also the radius of nanopore (RPore). Regarding the nanopore with specific RPore (RPore = RPore,c), the corresponding critical ion size (rIon,c) may be expected. As ion size is less than rIon,c, the energy barrier between the ion and nanopore is expected, rIon < rIon,c and RPore = RPore,c. With reference to rIon,c, the lower energy barrier may be expected for large ions. Additionally, as the ion passes through the filter, the direct interactions between the ion and channel, such as van der Waals interactions, become important, thereby influencing the energy barrier. In addition, the external driving force is necessarily added to the ion so that it can pass through the filter. This may be utilized to understand the ion selectivity of nanopores.

2.2. MD Simulations

From the above discussion, the ion selectivity of nanopores is closely related to the energy barrier as the ion enters the nanopore of the filter. Furthermore, the energy barrier between the ions and nanopore is related to both ion size and the radius of the pore. To understand the dependence of ion selectivity on the ion size, different ions (K+ and Na+) are forced to pass through the specific nanopore with a radius of 3.7 Å. In addition, to investigate the effects of nanopore pore radius on ion selectivity, MD simulations are also conducted on the systems in which a K+ ion enters the nanopores with different radii, such as 1.8 Å, 3.7 Å, and 5.6 Å.
In this work, the PMFs are calculated through ABF methods as a K+ ion perpendicularly approaching the nanopore with a radius of 3.7 Å (Figure 5). With a decrease in the distance between the K+ ions and nanopores, various minima can be found on the calculated PMFs (Figure 5). The minimum at 7.5 Å is found in the PMF, corresponding to a double water molecular layer located between the ion and nanopore (Figure 6a). Another minimum is located at 4.7 Å, attributed to the solvent-separated PMF, which means that only a water molecular layer can be found between the K+ ion and the graphene (Figure 6b).
From the calculated PMFs at 300 K, with decreasing distance between the K+ ion and graphene embedded with nanopores, energy barriers can be found between the neighboring minima (Figure 5). Of course, they may be related to the expulsion of a single water layer in the confined volume as the solutes moved closer together. Additionally, an obvious energy barrier can be found on the calculated PMFs as the distance between the K+ ions and nanopores is less than 5.0 Å, in which only one water layer can be found between them (Figure 6b), especially the distance less than 2.4 Å where K+ lies at the graphene surface (Figure 6c). Due to the obvious energy barrier between K+ ion and graphene, the ion is suppressed (or forbidden) from entering the nanopore of graphene (Figure 6d).
In addition, as the ion enters the nanopore, no water molecules can be found between the K+ ion and the nanopore of graphene, and the ion comes into contact with the graphene. This indicates that the energy barrier between the K+ ion and the nanopore may be closely related to the expulsion of interfacial water when the ion enters the nanopore. In other words, the above energy barrier may be related to the water molecules changing from interfacial to bulk water as the K+ ion approaches the nanopore of graphene. In combination with our recent studies [47,48], it may be concluded that this is closely related to hydrophobic interactions between the ion and nanopore.
Based on the above discussion on hydrophobic interactions, the dissolved K+ ion and graphene may be dispersed in water in the initial solvation process. To be more thermodynamically stable, water molecules can be expected between the K+ ion and graphene. In other words, it seems that there is a “repulsive” force as the K+ ion enters the filter, in which no water molecules can be found between them. Therefore, when the K+ ion enters the nanopore, the energy barrier may be due to the expulsion of interfacial water between the K+ ion and the graphene, and the K+ ion comes into contact with graphene in solutions. Of course, this is also in agreement with the above thermodynamic analysis on the energy barrier.
In our Raman spectroscopic studies [42,43,44], the Raman OH vibrations are mainly dependent on the local hydrogen-bonded network of a water molecule as three-dimensional hydrogen bondings appear in the water. Therefore, in ambient conditions, the dissolved solute mainly affects the structure of interfacial water. Of course, this is demonstrated by other studies on the structure and dynamics of water around ions using neutron and X-ray diffraction [53], X-ray absorption spectroscopy [54], femtosecond time-resolved infrared (fs-IR) vibrational spectroscopy [55,56,57], and optical Kerr effect spectroscopy [58]. These mean that the effects of ions on water may be largely limited to the first solvation shell.
The dissolved solutes primarily affect the structure of interfacial water. When they are dissolved in water, they can be divided into interfacial and bulk water. To investigate the structural reorganization between interfacial and bulk water, they are calculated as the K+ ion permeates the nanopore (Figure 7). With the decrease in the distance between the ion and nanopore, there is a decrease in the molecular number for interfacial water (Figure 7a), but an increase in the number for bulk water (Figure 7b). When the K+ ion enters the nanopore, the interfacial water associated with the ion and nanopore may be expelled into bulk water. Additionally, this is also related to the appearance of the energy barrier between the ion and nanopore.
To investigate the changes in hydrogen bondings when the K+ ion is forced to enter the nanopore of graphene, the hydrogen bondings of water are also determined. In this study, the geometrical definition of hydrogen bonding is utilized to determine the hydrogen bondings in water [59]. According to the geometrical definition, hydrogen bonding is considered to exist between two neighboring water molecules if the oxygen–oxygen distance (rOO) and ∠OOH angle between two water molecules are less than 3.5 Å and 30°, respectively.
From the MD simulations, the hydrogen bonding number (average number of hydrogen bonds per water molecule, nHB) can be calculated (Figure 8). In comparison with bulk water, due to the truncations of hydrogen bonding at the solute/water interface, the hydrogen bonding of interfacial water is less than that of bulk water (Figure 8). Therefore, the origin of hydrophobic interactions is due to the competition between the hydrogen bondings of bulk and interfacial water.
Regarding the nanopore with specific radius RPore, the energy barrier between the ion and nanopore is expected as the ion size is less than rIon,c, rIon < rIon,c and RPore = RPore,c. Additionally, the above energy barrier may be dependent on the ion size. To understand the dependence of the energy barrier on the ion size, MD simulations are carried out on different ions (K+ and Na+) passing through the specific nanopore with a radius of 3.7 Å. From the calculated PMFs, the obvious barrier can also be found as the distance between the Na+ ion and the filter is less than 3.1 Å (Figure 5). However, compared to the Na+ ion, a lower energy barrier can be expected as the K+ ion enters the nanopore.
From the above discussion, the origin of the energy barrier between the ion and nanopore is closely related to the expulsion of the interfacial water layer between them as the ion enters the nanopore. This is in accordance with the fact that the ion will lose the hydration shell and become essentially dehydrated as the ion enters the pore. The dissolved ion mainly affects the structure of the interfacial water, and the coordination number within the first hydration shell is related to the ion size. Regarding Na+, the number of water molecules in the first hydration layer is 5.4, which is less than the corresponding number of K+ (8.0). This may be applied to understand the dependence of the energy barrier on ion size.
Due to the energy barrier, when the ion enters the nanopore, the ion may be suppressed or forbidden from entering the channel. Therefore, the external driving force, such as the transmembrane electrical voltage or pressure, is necessarily added to the ion so that it may surmount the energy barrier and pass through the filter. Regarding the channel with a specific radius pore RPore, the energy barrier is dependent on the ion size. Under the same electrical voltage, it can be expected that the nanopore preferentially allows the specific ion with the lower barrier to pass through it. In comparison with small ions, a lower energy barrier is expected for larger ions. Therefore, compared to the smaller Na+ ions, the nanopore prefers the larger K+ ions to pass through the channel. This may be utilized to understand the ion selectivity of nanopores, such as KcsA for K+ over Na+.
In this study, a few MD simulations were conducted to investigate the effects of electrical voltage on ion selectivity (Supplementary Materials). In the simulations, various electrical voltages were added to the aqueous solution, which was composed of 10 KCl and 10 NaCl molecules. It is necessary to provide enough voltage (2.0 kcal·mol−1·Å−1·e−1) so that the K+ ions may firstly pass through the filter with a pore radius of 3.7 Å (Figure 9). With increasing voltage, more K+ ions may pass through the filter. However, some Na+ ions were found to penetrate the nanopore at high voltages. As the voltage is 3.0 kcal·mol−1·Å−1·e−1, the filter gave the suitable ion selectivity of K+/(K+ + Na+) (Figure 9). Regarding the specific ion, it is important to design the corresponding nanopore for the ion. Additionally, it is also necessary to choose the reasonable voltage to drive the specific ion to pass through the corresponding nanopore of the filter.
In addition, the energy barrier between the ion and nanopore may be affected by not only the ion size but also the pore radius of the nanopore (RPore < RPore,c and rIon = rIon,c). In this study, to investigate the dependence of ion selectivity on the radius of the nanopore, MD simulations were also conducted on the systems in which the K+ ion approaches the nanopore with different radii, namely 1.8 Å, 3.7 Å, and 5.6 Å. From the calculated PMFs (Figure 10), the energy barrier can be found as the K+ ion enters the nanopore with the different radii. However, with an increase in the radius of the nanopore, this decreases the energy barrier between the K+ ion and nanopore. This is in agreement with the above discussion on the effects of pore radius on the energy barrier.
From the calculated PMFs (Figure 10), a decrease in pore radius of the nanopore may lead to the increase in an energy barrier as the K+ ion enters the nanopore. Therefore, the strong transmembrane electrical voltage or driving pressure is necessarily added to the ion so that it can surmount the energy barrier and pass through the small nanopore. In fact, due to the strong energy barrier when K+ approaches the nanopore with a pore radius of 1.8 Å, the ion may be forbidden from entering the small nanopore even under high electrical voltage. Of course, the filter with this small pore radius is not suitably applied to separate the ion from others.
Additionally, an increase in the pore radius of the nanopore leads to a decrease in the energy barrier as K+ permeates the nanopore (Figure 10). In fact, no obvious energy barrier is found as K+ enters the nanopore with a pore radius of 5.6 Å. Therefore, the K+ ion may be expected to approach and freely enter the nanopore. In other words, this indicates that, with an increase in the pore radius, the hydrophobic interactions between the K+ ion and the nanopore may change from “repulsive” to “attractive” forces. Additionally, this also means that the nanopore with a pore radius of 5.6 Å may lose the ion selectivity for K+.
From the above discussion, it can be stated that the ion selectivity of the nanopore is associated with the energy barrier between the ion and the nanopore, arising from the transformation of water molecules from interfacial to bulk water as the ion enters the nanopore. Therefore, if the interfacial water of the ion is still preserved when it passes through the channel, it may be expected that no energy barrier can be found. In fact, it can be found that the interfacial water of the K+ ion can still be found as K+ permeates the nanopore with a pore radius of 5.6 Å (Figure 11). In addition, this also indicates that the effects of the dissolved ions on water structure may be mainly limited to the first hydration shell of the ion.
The ion selectivity of the nanopore is due to the difference in the energy barrier as various ions pass through the channel. Thermodynamically, the nanopore may prefer the specific ion with a lower barrier to permeate the filter. From this study, it can be derived that the energy barrier is closely related to the initial solvation process in which the Gibbs free energy of interfacial water is lower than that of bulk water. As the ions pass through the nanopore, they may become dehydrated. Therefore, the energy barrier is due to the water molecules transforming from interfacial to bulk water, which is dependent on not only the ion size but also the pore radius of the nanopore. Additionally, as the ion enters the channel, the energy barrier may also be affected by the direct interactions between the ion and nanopore. Due to the energy barrier between the ions and nanopore, the driving forces, such as the electrical voltage, may be necessarily added to the ions so that they may surmount the energy barrier and pass through the nanopore (Figure 12).

3. Method

3.1. Simulated Systems

The MD simulations were conducted using the program NAMD 2.12 [60]. The simulations were carried out on different systems, each containing a graphene filter embedded with a nanopore. The ion selectivity of the nanopore is related to the ion size and pore radius. After considering the pore radius of the KcsA channel, three nanopores with different radii were constructed: 1.8 Å, 3.7 Å, and 5.6 Å.
The empirical Chemistry at Harvard Macromolecular Mechanics (CHARMM) force field [61] was utilized to describe interatomic interactions. The water molecules were simulated using the intermolecular three-point potential (TIP3P) water model [62]. Non-bonded van der Waals interactions were smoothly switched to zero between 10 and 12 Å. The electrostatic interactions were simulated using the particle mesh Ewald (PME) algorithm.
In this study, the MD simulations were carried out in the NVT ensemble. The simulated temperature was 300 K, employing moderately damped Langevin dynamics. The simulated box was 40 Å × 40 Å × 60 Å. In addition, periodic boundary conditions were utilized in the three directions of Cartesian spaces. In the work, PMFs were determined through the adaptive biasing force (ABF) calculations [63,64,65,66,67,68]. The total time for ABF simulations was 80 ns. Additionally, a few simulations were also performed to study the ion selectivity of nanopores under various electrical voltages (Supplementary Materials). For MD simulations under electrical voltages, the total simulated time was 84 ns. Additionally, the simulated results were analyzed through VMD package [69].

3.2. ABF Calculations

In this work, PMFs were calculated to investigate the mechanism of ion selectivity of nanopores. These were carried out using NAMD with ABF [63,64,65,66,67,68] extensions integrated in the Collective Variables (CVs) module [70].
The ABF method is based on the thermodynamic integration of average force acting on coordinates, which is a combination of probability density and constraint force methods. In the ABF method, a biasing force opposing the actual force arising from system components is periodically applied to the reaction coordinate to generate what is effectively a random walk along the reaction coordinate (purely diffusive dynamics). For the ABF method, the free energy along a transition coordinate can be seen as a potential arising from the average force acting along the coordinate.
In this study, the target nanopore is fixed, and the test cation (K+ or Na+) is restrained to perpendicularly move to the center of graphene nanopore along Z axis. In the ABF calculations, the distances between the ions and nanopores are regarded as CVs, which are broken down into several consecutive windows, each with 2.5–2.7 Å wide.

4. Conclusions

Based on thermodynamics analysis and MD simulations, this work is devoted to investigating the mechanism of ion selectivity in nanopores. From the work, the following conclusions are derived:
(1) The ion selectivity of the channel is related to the rearrangement of water molecules as the ion enters the channel. In combination with our recent studies, the ion selectivity of the nanopore may be related to hydrophobic interactions between the ion and nanopore.
(2) In our recent works on hydrophobic interactions, with reference to Rc, it was divided into initial and hydrophobic solvation processes. In the hydrophobic process (>Rc), the ions are expected to freely pass through the nanopore (rIon = rIon,c, RPore > RPore,c) or are forbidden from passing through the channel (rIon > rIon,c, RPore = RPore,c). It seems that there are the “attractive” forces between the ions and nanopore.
(3) In the initial process (<Rc), the energy barrier can be expected as the ion enters the channel (rIon<rIon,c, RPore = RPore,c), which is related to the ion selectivity of the nanopore. It is found that the energy barrier is related to the ion size and pore radius. Additionally, the ion selectivity of the nanopore is also affected by the external forces added to the ions to overcome the energy barrier.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules29040853/s1, A few MD simulations were performed to study the ion selectivity of the nanopore under various electrical voltages. Some simulated trajectories are provided in the Supplementary Data.

Author Contributions

Writing—review and editing, Y.-N.C.; writing—review and editing, Y.-Z.L.; conceptualization, formal analysis, investigation, resources, data curation, writing—original draft preparation, supervision, writing—review and editing, Q.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article and Supplementary Materials.

Acknowledgments

The editor and reviewers are greatly appreciated for providing good suggestions to revise the paper.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Hille, B. Ionic Channels of Excitable Membranes, 3rd ed.; Sinauer Associates Inc.: Sunderland, UK, 2001. [Google Scholar]
  2. Konatham, D.; Yu, J.; Ho, T.A.; Striolo, A. Simulation insights for graphene-based water desalination membranes. Langmuir 2013, 29, 11884–11897. [Google Scholar] [CrossRef]
  3. Wang, Y.H.; He, Z.J.; Gupta, K.M.; Shi, Q.; Lu, R.F. Molecular dynamics study on water desalination through functionalized nanoporous grapheme. Carbon 2017, 116, 120–127. [Google Scholar] [CrossRef]
  4. Nguyen, C.T.; Beskok, A. Charged nanoporous graphene membranes for water desalination. Phys. Chem. Chem. Phys. 2019, 21, 9483–9494. [Google Scholar] [CrossRef]
  5. Mackinnon, R. Potassium channels and the atomic basis of selective ion conduction (Nobel lecture). Angew. Chem. Int. Ed. 2004, 43, 4265–4277. [Google Scholar] [CrossRef] [PubMed]
  6. Zhou, Y.; Morais-Cabral, J.H.; Kaufman, A.; MacKinnon, R. Chemistry of ion coordination and hydration revealed by a K+ channel-Fab complex at 2.0 A resolution. Nature 2001, 414, 43–48. [Google Scholar] [CrossRef]
  7. Doyle, D.A.; Cabral, J.M.; Pfuetzner, R.A.; Kuo, A.; Gulbis, J.M.; Cohen, S.L.; Chait, B.T.; MacKinnon, R. The structure of the potassium channel: Molecular basis of K+ conduction and selectivity. Science 1998, 280, 69–77. [Google Scholar] [CrossRef]
  8. Heginbotham, L.; MacKinnon, R. Conduction properties of the cloned Shaker channel. Biophys. J. 1993, 65, 2089–2096. [Google Scholar] [CrossRef] [PubMed]
  9. LeMasurier, M.; Heginbotham, L.; Miller, C. KscA: It’s an ion channel. J. Gen. Physiol. 2001, 118, 303–314. [Google Scholar] [CrossRef]
  10. Medovoy, D.; Perozo, E.; Roux, B. Multi-ion free energy landscapes underscore the microscopic mechanism of ion selectivity in the KcsA channel. BBA Biomembranes 2016, 1858, 1722–1732. [Google Scholar] [CrossRef]
  11. Dudev, T.; Lim, C. Determinants of K+ vs Na+ selectivity in potassium channels. J. Am. Chem. Soc. 2009, 131, 8092–8101. [Google Scholar] [CrossRef] [PubMed]
  12. Dixit, P.D.; Merchant, S.; Asthagiri, D. Ion selectivity in the KcsA potassium channel from the perspective of the ion binding site. Biophys. J. 2009, 96, 2138–2145. [Google Scholar] [CrossRef]
  13. Bostick, D.L.; Arora, K.; Brooks III, C.L. K+/Na+ selectivity in toy cation binding site models is determined by the ‘host’. Biophys. J. 2009, 96, 3887–3896. [Google Scholar] [CrossRef] [PubMed]
  14. Bostick, D.L.; Brooks III, C.L. Selectivity in K+ channels is due to topological control of the permeant ion’s coordinated state. Proc. Natl. Acad. Sci. USA 2007, 104, 9260–9265. [Google Scholar] [CrossRef]
  15. Yu, H.; Yu, S.; Roux, N.B. Two mechanisms of ion selectivity in protein binding sites. Proc. Natl. Acad. Sci. USA 2010, 107, 20329–20334. [Google Scholar] [CrossRef] [PubMed]
  16. Fowler, P.W.; Tai, K.; Sansom, M.S.P. The selectivity of K+ Ion channels: Testing the hypotheses. Biophys. J. 2008, 95, 5062–5072. [Google Scholar] [CrossRef]
  17. Lockless, S.W.; Zhou, M.; MacKinnon, R. Structural and thermodynamic properties of selective ion binding in a K+ channel. PLoS Biol. 2007, 5, 1079–1088. [Google Scholar] [CrossRef] [PubMed]
  18. Furini, S.; Domene, C. Selectivity and permeation of alkali metal ions in K+-channels. J. Mol. Biol. 2011, 409, 867–878. [Google Scholar] [CrossRef] [PubMed]
  19. Liu, S.; Bian, X.; Lockless, S.W. Preferential binding of K+ ions in the selectivity filter at equilibrium explains high selectivity of K+ channels. J. Gen. Physiol. 2012, 140, 671–679. [Google Scholar] [CrossRef]
  20. Valiyaveetil, F.I.; Leonetti, M.; Muir, T.W.; MacKinnon, R. Ion selectivity in a semisynthetic K+ channel locked in the conductive conformation. Science 2006, 314, 1004–1007. [Google Scholar] [CrossRef]
  21. Thompson, A.N.; Kim, I.; Panosian, T.D.; Iverson, T.M.; Allen, T.W.; Nimigean, C.M. Mechanism of potassium-channel selectivity revealed by Na+ and Li+ binding sites within the KcsA pore. Nat. Struct. Mol. Biol. 2009, 16, 1317–1324. [Google Scholar] [CrossRef]
  22. Egwolf, B.; Roux, B. Ion selectivity of the KcsA channel: A perspective from multi-ion free energy landscapes. J. Mol. Biol. 2010, 401, 831–842. [Google Scholar] [CrossRef] [PubMed]
  23. Noskov, S.Y.; Roux, B. Ion selectivity in potassium channels. Biophys. Chem. 2006, 124, 279–291. [Google Scholar] [CrossRef] [PubMed]
  24. Andersen, O.S. Perspectives on: Ion selectivity. J. Gen. Physiol. 2011, 137, 5393–5395. [Google Scholar] [CrossRef] [PubMed]
  25. Dixit, P.D.; Asthagiri, D. Thermodynamics of ion selectivity in the KcsA K+ channel. J. Gen. Physiol. 2011, 137, 427–433. [Google Scholar] [CrossRef] [PubMed]
  26. Naranjo, D.; Moldenhauer, H.; Pincuntureo, M.; Díaz-Franulic, I. Pore size matters for potassium channel conductance. J. Gen. Physiol. 2016, 148, 277–291. [Google Scholar] [CrossRef]
  27. Horn, R.; Roux, B.; Åqvist, J. Permeation redux: Thermodynamics and kinetics of ion movement through potassium channels. Biophys. J. 2014, 106, 1859–1863. [Google Scholar] [CrossRef]
  28. Kim, I.; Allen, T.W. On the selective ion binding hypothesis for potassium channels. Proc. Natl. Acad. Sci. USA 2011, 108, 17963–17968. [Google Scholar] [CrossRef]
  29. Beckstein, O.; Tai, K.; Sansom, M.S.P. Not ions alone: Barriers to ion permeation in nanopores and channels. J. Am. Chem. Soc. 2004, 126, 14694–14695. [Google Scholar] [CrossRef]
  30. Roux, B.; Bernèche, S.; Egwolf, B.; Lev, B.; Noskov, S.Y.; Rowley, C.N.; Yu, H. Ion selectivity in channels and transporters. J. Gen. Physiol. 2011, 137, 415–426. [Google Scholar] [CrossRef]
  31. Roux, B. Ion channels and ion selectivity essays in biochemistry. Essays Biochem. 2017, 61, 201–209. [Google Scholar]
  32. He, Z.J.; Zhou, J.; Lu, X.H.; Corry, B. Bioinspired graphene nanopores with voltage-tunable ion selectivity for Na+ and K+. ACS Nano 2013, 7, 10148–10157. [Google Scholar] [CrossRef] [PubMed]
  33. Kang, Y.; Zhang, Z.; Shi, H.; Zhang, J.; Liang, L.; Wang, Q.; Ågren, H.; Tu, Y. Na+ and K+ ion selectivity by size-controlled biomimetic graphene nanopores. Nanoscale 2014, 6, 10666–10672. [Google Scholar] [CrossRef] [PubMed]
  34. Bezanilla, F.; Armstrong, C.M. Negative conductance caused by entry of sodium and cesium ions into the potassium channels of squid axons. J. Gen. Physiol. 1972, 60, 588–608. [Google Scholar] [CrossRef]
  35. Eisenman, G.; Horn, R. Ionic selectivity revisited: The role of kinetic and equlibrium processes in ion permeation through channels. J. Membr. Biol. 1983, 76, 197–225. [Google Scholar] [CrossRef] [PubMed]
  36. Eisenman, G. Cation selective electrodes and their mode of operation. Biophys. J. 1962, 2, 259–323. [Google Scholar] [CrossRef]
  37. Noskov, S.Y.; Bernèche, S.; Roux, B. Control of ion selectivity in potassium channels by electrostatic and dynamic properties of carbonyl ligands. Nature 2004, 431, 830–834. [Google Scholar] [CrossRef]
  38. Kraszewski, S.; Boiteux, C.; Ramseyer, C.; Girardet, C. Determination of the charge profile in the KcsA selectivity filter using ab initio calculations and molecular dynamics simulations. Phys. Chem. Chem. Phys. 2009, 11, 8606–8613. [Google Scholar] [CrossRef]
  39. Varma, S.; Rempe, S.B. Tuning ion coordination architectures to enable selective partitioning. Biophys. J. 2007, 93, 1093–1099. [Google Scholar] [CrossRef]
  40. Thomas, M.; Jayatilaka, D.; Corry, B. The predominant role or coordination number in potassium channel selectivity. Biophys. J. 2007, 93, 2635–2643. [Google Scholar] [CrossRef]
  41. Varma, S.; Sabo, D.; Rempe, S.B. K+/Na+ selectivity in K channels and valinomycin: Over-coordination versus cavity-size constraints. J. Mol. Biol. 2008, 376, 13–22. [Google Scholar] [CrossRef]
  42. Sun, Q. The Raman OH stretching bands of liquid water. Vib. Spectrosc. 2009, 51, 213–217. [Google Scholar] [CrossRef]
  43. Sun, Q. Raman spectroscopic study of the effects of dissolved NaCl on water structure. Vib. Spectrosc. 2012, 62, 110–114. [Google Scholar] [CrossRef]
  44. Sun, Q. Local statistical interpretation for water structure. Chem. Phys. Lett. 2013, 568, 90–94. [Google Scholar] [CrossRef]
  45. Sun, Q.; Guo, Y. Vibrational sum frequency generation spectroscopy of the air/water interface. J. Mol. Liq. 2016, 213, 28–32. [Google Scholar] [CrossRef]
  46. Sun, Q. The Effects of dissolved hydrophobic and hydrophilic groups on water structure. J. Solut. Chem. 2020, 49, 1473–1484. [Google Scholar] [CrossRef]
  47. Sun, Q. The physical origin of hydrophobic effects. Chem. Phys. Lett. 2017, 672, 21–25. [Google Scholar] [CrossRef]
  48. Sun, Q.; Su, X.W.; Cheng, C.B. The dependence of hydrophobic interactions on the solute size. Chem. Phys. 2019, 516, 199–205. [Google Scholar] [CrossRef]
  49. Sun, Q.; Zhang, M.X.; Cui, S. The structural origin of hydration repulsive force. Chem. Phys. Lett. 2019, 714, 30–36. [Google Scholar] [CrossRef]
  50. Sun, Q.; Wang, W.; Cui, S. Directional nature of hydrophobic interactions: Implications for the mechanism of molecular recognition. Chem. Phys. 2021, 547, 111200. [Google Scholar] [CrossRef]
  51. Sun, Q. The hydrophobic effects: Our current understanding. Molecules 2022, 27, 7009. [Google Scholar] [CrossRef]
  52. Stanley, H.E.; Teixeira, J. Interpretation of the unusual behavior of H2O and D2O at low temperatures: Tests of a percolation model. J. Chem. Phys. 1980, 73, 3404–3422. [Google Scholar] [CrossRef]
  53. Collins, K.D.; Neilson, G.W.; Enderby, J.E. Ions in water: Characterizing the forces that control chemical processes and biological structure. Biophys. Chem. 2007, 128, 95–104. [Google Scholar] [CrossRef] [PubMed]
  54. Cappa, C.D.; Smith, J.D.; Messer, B.M.; Cohen, R.C.; Saykally, R.J. Effects of cations on the hydrogen bond network of liquid water:  New results from X-ray absorption spectroscopy of liquid microjets. J. Phys. Chem. B 2006, 110, 5301–5309. [Google Scholar] [CrossRef] [PubMed]
  55. Omta, A.W.; Kropman, M.F.; Woutersen, S.; Bakker, H.J. Negligible effect of ions on the hydrogen-bond structure in liquid water. Science 2003, 301, 347–349. [Google Scholar] [CrossRef] [PubMed]
  56. Park, S.; Fayer, M.D. Hydrogen bond dynamics in aqueous NaBr solutions. Proc. Natl. Acad. Sci. USA 2007, 104, 16731–16738. [Google Scholar] [CrossRef] [PubMed]
  57. Moilanen, D.E.; Wong, D.; Rosenfeld, D.E.; Fenn, E.E.; Fayer, M.D. Ion–water hydrogen-bond switching observed with 2D IR vibrational echo chemical exchange spectroscopy. Proc. Natl. Acad. Sci. USA 2009, 106, 375–380. [Google Scholar] [CrossRef] [PubMed]
  58. Turton, D.A.; Hunger, J.; Hefter, G.; Buchner, R.; Wynne, K. Glasslike behavior in aqueous electrolyte solutions. J. Chem. Phys. 2008, 128, 161102. [Google Scholar] [CrossRef]
  59. Kumar, R.; Schmidt, J.R.; Skinner, J.L. Hydrogen bonding definitions and dynamics in liquid water. J. Chem. Phys. 2007, 126, 204107. [Google Scholar] [CrossRef]
  60. Phillips, J.C.; Braun, R.; Wang, W.; Gumbart, J.; Tajkhorshid, E.; Villa, E.; Chipot, C.; Skeel, R.D.; Kale, L.; Schulten, K. Scalable molecular dynamics with NAMD. J. Comput. Chem. 2005, 26, 1781–1802. [Google Scholar] [CrossRef]
  61. MacKerell, A.D.; Bashford, D.; Bellott, M.; Dunbrack, R.L.; Evanseck, J.D.; Field, M.J.; Fischer, S.; Gao, J.; Guo, H.; Ha, S.; et al. All-atom empirical potential for molecular modeling and dynamics studies of proteins. J. Phys. Chem. B 1998, 102, 3586–3616. [Google Scholar] [CrossRef]
  62. Price, D.J.; Brooks III, C.L. A modified TIP3P water potential for simulation with Ewald summation. J. Chem. Phys. 2004, 121, 10096–10103. [Google Scholar] [CrossRef]
  63. Darve, E.; Pohorille, A. Calculating free energies using average force. J. Chem. Phys. 2001, 115, 9169–9183. [Google Scholar] [CrossRef]
  64. Rodriguez-Gomeza, D.; Darveb, E.; Pohorille, A. Assessing the efficiency of free energy calculation methods. J. Chem. Phys. 2004, 120, 3563. [Google Scholar] [CrossRef]
  65. Darve, E.; Rodríguez-Gómez, D.; Pohorille, A. Adaptive biasing force method for scalar and vector free energy calculations. J. Chem. Phys. 2008, 128, 144120. [Google Scholar] [CrossRef]
  66. Hénin, J.; Chipot, C. Overcoming free energy barriers using unconstrained molecular dynamics simulations. J. Chem. Phys. 2004, 121, 2904–2914. [Google Scholar] [CrossRef] [PubMed]
  67. Hénin, J.; Fiorin, G.; Chipot, C.; Klein, M.L. Exploring multidimensional free energy landscapes using time-dependent biases on collective variables. J. Chem. Theory Comput. 2010, 6, 35–47. [Google Scholar] [CrossRef] [PubMed]
  68. Comer, J.; Gumbart, J.C.; Hénin, J.; Leliévre, T.; Pohorille, A.; Chipot, C. The adaptive biasing force method: Everything you always wanted to know but were afraid to ask? J. Phys. Chem. B 2015, 119, 1129–1151. [Google Scholar] [CrossRef] [PubMed]
  69. Humphrey, W.; Dalke, A.; Schulten, K. VMD-Visual molecular dynamics. J. Mol. Graphics 1996, 14, 33–38. [Google Scholar] [CrossRef] [PubMed]
  70. Fiorin, G.; Klein, M.L.; Hénin, J. Using collective variables to drive molecular dynamics simulations. Mol. Phys. 2013, 111, 3345–3362. [Google Scholar] [CrossRef]
Figure 1. Hydration free energy at 293 K and 0.1 MPa. With reference to Rc, it can be divided into initial and hydrophobic solvation processes. The Rc is drawn as a dashed line. (I) Initial solvation process. (II) Hydrophobic solvation process.
Figure 1. Hydration free energy at 293 K and 0.1 MPa. With reference to Rc, it can be divided into initial and hydrophobic solvation processes. The Rc is drawn as a dashed line. (I) Initial solvation process. (II) Hydrophobic solvation process.
Molecules 29 00853 g001
Figure 2. (a) To investigate the ion selectivity of nanopore, the ion moves perpendicularly to the filter embedded with a nanopore with radius R. Both side and top views are shown. (b) In thermodynamics, the Gibbs free energy (ΔGIon-filter) related to the ion and filter equals the difference in Gibbs energy between the ion and membrane (ΔGIon-membrane), and that between the ion and filled circle with radius R (ΔGIon-circle). Additionally, d is the distance between the centers of ion and filter, and sep means the separation between the ion and nanopore.
Figure 2. (a) To investigate the ion selectivity of nanopore, the ion moves perpendicularly to the filter embedded with a nanopore with radius R. Both side and top views are shown. (b) In thermodynamics, the Gibbs free energy (ΔGIon-filter) related to the ion and filter equals the difference in Gibbs energy between the ion and membrane (ΔGIon-membrane), and that between the ion and filled circle with radius R (ΔGIon-circle). Additionally, d is the distance between the centers of ion and filter, and sep means the separation between the ion and nanopore.
Molecules 29 00853 g002
Figure 3. (a) Regarding the hydrophobic interactions between ion and filter, the critical distance (dc) can be expected, which is related to the pore radius (RPore,c) and ion size (rIon,c). (b) For the specific ion, the “attractive” force between the ion and filter is expected as the pore radius is larger than RPore,c. (c) As the specific ion approaches the nanopore, the “repulsive” force is expected as the pore radius is less than RPore,c. It seems that there is an energy barrier between the ion and nanopore.
Figure 3. (a) Regarding the hydrophobic interactions between ion and filter, the critical distance (dc) can be expected, which is related to the pore radius (RPore,c) and ion size (rIon,c). (b) For the specific ion, the “attractive” force between the ion and filter is expected as the pore radius is larger than RPore,c. (c) As the specific ion approaches the nanopore, the “repulsive” force is expected as the pore radius is less than RPore,c. It seems that there is an energy barrier between the ion and nanopore.
Molecules 29 00853 g003
Figure 4. (a) As the ion approaches the nanopore, dc can be expected. (b) The ion is banned from passing through the nanopore when the ion is larger than rIon,c. (c,d) When the ion enters the channel with a specific radius, the energy barrier is expected when the ion size is less than rIon,c. For the specific nanopore with radius R, the energy barrier is related to the ion size. A lower energy barrier is expected for a larger ion.
Figure 4. (a) As the ion approaches the nanopore, dc can be expected. (b) The ion is banned from passing through the nanopore when the ion is larger than rIon,c. (c,d) When the ion enters the channel with a specific radius, the energy barrier is expected when the ion size is less than rIon,c. For the specific nanopore with radius R, the energy barrier is related to the ion size. A lower energy barrier is expected for a larger ion.
Molecules 29 00853 g004
Figure 5. The PMFs as K+ and Na+ permeate the nanopore with a radius of 3.7 Å. During K+ approaches the nanopore, various minima are found at 7.5 Å, 4.7 Å, and 2.4 Å, which are labeled as a, b, and c. Compared with Na+, a lower energy barrier is found as K+ enters the nanopore. Due to the difference in energy barrier, this leads to the ion selectivity for K+ over Na+.
Figure 5. The PMFs as K+ and Na+ permeate the nanopore with a radius of 3.7 Å. During K+ approaches the nanopore, various minima are found at 7.5 Å, 4.7 Å, and 2.4 Å, which are labeled as a, b, and c. Compared with Na+, a lower energy barrier is found as K+ enters the nanopore. Due to the difference in energy barrier, this leads to the ion selectivity for K+ over Na+.
Molecules 29 00853 g005
Figure 6. Snapshots of water molecules between K+ and the nanopore at various distances shown in Figure 5: (a) 7.5 Å, (b) 4.7 Å, and (c) 2.4 Å. (d) The nanopore with a pore radius of 3.7 Å is also shown.
Figure 6. Snapshots of water molecules between K+ and the nanopore at various distances shown in Figure 5: (a) 7.5 Å, (b) 4.7 Å, and (c) 2.4 Å. (d) The nanopore with a pore radius of 3.7 Å is also shown.
Molecules 29 00853 g006
Figure 7. The changes in interfacial and bulk water as K+ passes through the nanopore. The letters a, b, and c mean the various distances between K+ and nanopore as shown in Figure 5.
Figure 7. The changes in interfacial and bulk water as K+ passes through the nanopore. The letters a, b, and c mean the various distances between K+ and nanopore as shown in Figure 5.
Molecules 29 00853 g007
Figure 8. The hydrogen bondings number per water molecule in: (a) interfacial water and bulk water, and (b) total water as K+ permeates the nanopore.
Figure 8. The hydrogen bondings number per water molecule in: (a) interfacial water and bulk water, and (b) total water as K+ permeates the nanopore.
Molecules 29 00853 g008
Figure 9. The effects of electrical voltage on the ion selectivity of a nanopore with a pore radius of 3.7 Å. The inlet shows the ion selectivity for K+/(K+ + Na+).
Figure 9. The effects of electrical voltage on the ion selectivity of a nanopore with a pore radius of 3.7 Å. The inlet shows the ion selectivity for K+/(K+ + Na+).
Molecules 29 00853 g009
Figure 10. The PMFs as the K+ ion approaches the nanopores with different radii: 1.8 Å, 3.7 Å, and 5.6 Å.
Figure 10. The PMFs as the K+ ion approaches the nanopores with different radii: 1.8 Å, 3.7 Å, and 5.6 Å.
Molecules 29 00853 g010
Figure 11. The nanopore with a pore radius of 1.8 Å (a) and 5.6 Å (b). (c,d) No obvious obstacle is found as K+ passes through the nanopore with a radius of 5.6 Å; this is because the interfacial water of K+ is preserved when it enters the nanopore.
Figure 11. The nanopore with a pore radius of 1.8 Å (a) and 5.6 Å (b). (c,d) No obvious obstacle is found as K+ passes through the nanopore with a radius of 5.6 Å; this is because the interfacial water of K+ is preserved when it enters the nanopore.
Molecules 29 00853 g011
Figure 12. The mechanism of ion selectivity of nanopore. (a) In combination with our recent studies, the critical pore radius (RPore,c) and ion size (rIon,c) are expected during the ion approaches the nanopore. (b) During the ion enters the filter, energy barrier is expected as ion size being less than rIon,c. (c) External force is necessarily added on the ion so that it enters the nanopore.
Figure 12. The mechanism of ion selectivity of nanopore. (a) In combination with our recent studies, the critical pore radius (RPore,c) and ion size (rIon,c) are expected during the ion approaches the nanopore. (b) During the ion enters the filter, energy barrier is expected as ion size being less than rIon,c. (c) External force is necessarily added on the ion so that it enters the nanopore.
Molecules 29 00853 g012
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Chen, Y.-N.; Liu, Y.-Z.; Sun, Q. The Molecular Mechanism of Ion Selectivity in Nanopores. Molecules 2024, 29, 853. https://doi.org/10.3390/molecules29040853

AMA Style

Chen Y-N, Liu Y-Z, Sun Q. The Molecular Mechanism of Ion Selectivity in Nanopores. Molecules. 2024; 29(4):853. https://doi.org/10.3390/molecules29040853

Chicago/Turabian Style

Chen, Yan-Nan, Yu-Zhen Liu, and Qiang Sun. 2024. "The Molecular Mechanism of Ion Selectivity in Nanopores" Molecules 29, no. 4: 853. https://doi.org/10.3390/molecules29040853

APA Style

Chen, Y. -N., Liu, Y. -Z., & Sun, Q. (2024). The Molecular Mechanism of Ion Selectivity in Nanopores. Molecules, 29(4), 853. https://doi.org/10.3390/molecules29040853

Article Metrics

Back to TopTop