Influence of the Dielectric Constant on the Ionic Current Rectification of Bipolar Nanopores.

In this paper, we investigate how the dielectric constant, ϵ, of an electrolyte solvent influences the current rectification characteristics of bipolar nanopores. It is well recognized that bipolar nanopores with two oppositely charged regions rectify current when exposed to an alternating electric potential difference. Here, we consider dilute electrolytes with NaCl only and with a mixture of NaCl and charged nanoparticles. These systems are studied using two levels of description, all-atom explicit water molecular dynamics (MD) simulations and coarse-grained implicit solvent MD simulations. The charge density and electric potential profiles and current-voltage relationship predicted by the implicit solvent simulations with ϵ = 11.3 show good agreement with the predictions from the explicit water simulations. Under nonequilibrium conditions, the predictions of the implicit solvent simulations with a dielectric constant closer to the one of bulk water are significantly different from the predictions obtained with the explicit water model. These findings are closely aligned with experimental data on the dielectric constant of water when confined to nanometric spaces, which suggests that ϵ decreases significantly compared to its value in the bulk. Moreover, the largest electric current rectification is observed in systems containing nanoparticles when ϵ = 78.8. Using enhanced sampling, we have shown that this larger rectification arises from the presence of a significantly deeper minimum in the free energy of the system with a larger ϵ, and when a negative voltage bias is applied. Since implicit solvent models and mean-field continuum theories are often used to design Janus membranes based on bipolar nanopores, this work highlights the importance of properly accounting for the effects of confinement on the dielectric constant of the electrolyte solvent. The results presented here indicate that the dielectric constant in implicit solvent simulations may be used as an adjustable parameter to approximately account for the effects of nanometric confinement on aqueous electrolyte solvents.

Ionic Transport in Electrostatic Janus Membranes. An Explicit Solvent Molecular Dynamic Simulation.

Janus, or two-sided, charged membranes offer promise as ionic current rectifiers. In such systems, pores consisting of two regions of opposite charge can be used to generate a current from a gradient in salinity. The efficiency of nanoscale Janus pores increases dramatically as their diameter becomes smaller. However, little is known about the underlying transport processes, particularly under experimentally accessible conditions. In this work, we examine the molecular basis for rectification in Janus nanopores using an applied electric field. Molecular simulations with explicit water and ions are used to examine the structure and dynamics of all molecular species in aqueous electrolyte solutions. For several macroscopic observables, the results of such simulations are consistent with experimental observations on asymmetric membranes. Our analysis reveals a number of previously unknown features, including a pronounced local reorientation of water molecules in the pores, and a segregation of ionic species that had not been anticipated by previously reported continuum analyses of Janus pores. Using these insights, a model is proposed for ionic current rectification in which electric leakage at the pore entrance controls net transport.

A structural indicator for water built upon potential energy considerations.

We introduce a parameter-free structural indicator to classify local environments of water molecules in stable and supercooled liquid states, which reveals a clear two-peak distribution of local properties. The majority of molecules are tetrahedrally coordinated (T molecules), via low-energy hydrogen bonds. The minority component, whose relative concentration decreases with a decrease in the temperature at constant pressure, is characterized by prevalently three-coordinated molecules, giving rise to a distorted local network around them (D molecules). The inter-conversion between T and D molecules explains the increasing specific heat at constant pressure on cooling. The local structure around a T molecule resembles the one found experimentally in low-density amorphous ice (a network structure mostly composed by T molecules), while the local structure around a D molecule is reminiscent of the structural properties of high-density amorphous ice (a network structure composed by a mixture of T and D molecules).

Comparing the performance of two structural indicators for different water models while seeking for connections between structure and dynamics in the glassy regime.

In this work, we compare the performance of two structural indicators based on the degree of translational order up to the second coordination shell in three water models: SPC/E, TIP4P/2005, and TIP5P. Beyond directly contrasting their distributions for different temperatures to evidence their usefulness in estimating the fraction of structured and unstructured molecules and, when possible, their classification capability, we also correlate them with an indirect measure of structural constraint: the dynamic propensity. Furthermore, this procedure enables us to show the existence of evident correlations between structural and dynamical information. More specifically, we find that locally structured molecules display a preference for low dynamic propensity values and, more conspicuously, that locally unstructured molecules are extremely subject to high dynamic propensity. This result is particularly relevant for the supercooled regime where the establishment of firm links between the structure and dynamics has remained rather elusive since the occurrence of dynamics that vary in orders of magnitude upon supercooling usually contrast with barely noticeable overall structural changes.

Structural features of high-local-density water molecules: Insights from structure indicators based on the translational order between the first two molecular shells.

The two-liquids scenario for liquid water assumes the existence of two competing preferential local molecular structural states characterized by either low or high local density. While the former is expected to present good local order thus involving privileged structures, the latter is usually regarded as conforming a high-entropy unstructured state. A main difference in the local arrangement of such "classes" of water molecules can be inferred from the degree of translational order between the first and second molecular shells. This is so, since the low-local-density molecules present a clear gap between the first two shells while in the case of the high-local-density ones, one or more molecules from the second shell have collapsed toward the first one, thus populating the intershell region. Some structural indicators, like the widely employed local structure index and the recently introduced ζ index, have been devised precisely on the basis of this observation, being successful in detecting well-structured low-local-density molecules. However, the nature of the high-local-density state has been mainly disregarded over the years. In this work we employ molecular dynamics simulations for two water models (the extended simple point charge model and the five-site model) at the liquid and supercooled regimes combined with the inherent dynamics approach (energy minimizations of the instantaneous configurations) in order to both rationalize the detailed structural and topological information that these indicators provide and to advance in our understanding of the high-density state.

Size dependence of dynamic fluctuations in liquid and supercooled water.

We study the evolution of dynamic fluctuations averaged over different space lengths and time scales to characterize spatially and temporally heterogeneous behavior of TIP4P/2005 water in liquid and supercooled states. Analyzing a 250 000 molecules simulated system, we provide evidence of the existence, upon supercooling, of a significant enhancement of spatially localized dynamic fluctuations stemming from regions of correlated mobile molecules. We show that both the magnitude of the departure from the value expected for the system-size dependence of an uncorrelated system and the system size at which such a trivial regime is finally recovered clearly increase upon supercooling. This provides a means to estimate an upper limit to the maximum length scale of influence of the regions of correlated mobile molecules. Notably, such an upper limit grows two orders of magnitude on cooling, reaching a value corresponding to a few thousand molecules at the lowest investigated temperature.

Studies on electrostatic interactions within model nano-confined aqueous environments of different chemical nature.

We study the potential of mean force for pairs of parallel flat surfaces with attractive electrostatic interactions by employing model systems functionalized with different charged, hydrophobic and hydrophilic groups. We study the way in which the local environment (hydrophobic or hydrophilic moieties) modulates the interaction between the attractive charged groups on the plates by removing or attracting nearby water and thus screening or not the electrostatic interaction. To explicitly account for the role of the solvent and the local hydrophobicity, we also perform studies in vacuo. Additionally, the results are compared to that for non-charged plates in order to single out and rationalize the non-additivity of the different non-covalent interactions. Our simulations demonstrate that the presence of neighboring hydrophobic groups promote water removal in the vicinity of the charged groups, thus enhancing charge attraction upon self-assembly. This role of the local hydrophobicity modulating electrostatic interactions is consistent with recent qualitative descriptions in the protein binding context.

Structure and dynamics of high- and low-density water molecules in the liquid and supercooled regimes.

By combining the local structure index with potential energy minimisations we study the local environment of the water molecules for a couple of water models, TIP5P-Ew and SPC/E, in order to characterise low- and high-density "species". Both models show a similar behaviour within the supercooled regime, with two clearly distinguishable populations of unstructured and structured molecules, the fraction of the latter increasing with supercooling. Additionally, for TIP5P-Ew, we find that the structured component vanishes quickly at the normal liquid regime (above the melting temperature). Thus, while SPC/E provides a fraction of structured molecules similar to that found in X-ray experiments, we show that TIP5P-Ew underestimates such value. Moreover, unlike SPC/E, we demonstrate that TIP5P-Ew does not follow the linear dependence of the logarithm of the structured fraction with inverse temperature, as predicted by the two-order parameter model. Finally, we link structure to dynamics by showing that there exists a strong correlation between structural fluctuation and dynamics in the supercooled state with spatial correlations in both static and dynamic quantities.

Hydrophilic behavior of graphene and graphene-based materials.

Graphene and the graphene-based materials like graphite, carbon nanotubes, and fullerenes are not only usually regarded as hydrophobic but also have been widely employed as paradigms for the investigation of the behavior of water under nonpolar confinement, a question of major concern for fields ranging from biology to materials design. However, some experimental and theoretical insights seem to contradict, at least partially, such a picture. In this work, we will provide firm evidence for a neat hydrophilic nature of graphene surfaces. Our molecular dynamics studies will demonstrate that parallel graphene sheets present a strong tendency to remain fully hydrated for moderately long times (even when the equilibrium state is indeed the collapse of the plates), and thus, they are less prone to self-assembly than model hydrophobic surfaces we shall employ as control which readily undergo a hydrophobic collapse. Potential of mean force calculations will indeed make evident that the solvent exerts a repulsive contribution on the self-assembly of graphene surfaces. Moreover, we shall also quantify graphene hydrophilicity by means of the calculation of water density at two pressures and water density fluctuations. This latter study has never been performed on graphene and represents a means both to confirm and to quantify its neat hydrophilic behavior. We shall also make evident the relevance of the mildly attractive water-carbon interactions, since their artificial weakening will be shown to revert from typically hydrophilic to typically hydrophobic behavior.

Productive induced metastability in allosteric modulation of kinase function.

Allosteric modulators of kinase function are of considerable pharmacological interest as blockers or agonists of key cell-signaling pathways. They are gaining attention due to their purported higher selectivity and efficacy relative to ATP-competitive ligands. Upon binding to the target protein, allosteric inhibitors promote a conformational change that purposely facilitates or hampers ATP binding. However, allosteric binding remains a matter of contention because the binding site does not fit with a natural ligand (i.e. ATP or phosphorylation substrate) of the protein. In this study, we show that allosteric binding occurs by means of a local structural motif that promotes association with the ligand. We specifically show that allosteric modulators promote a local metastable state that is stabilized upon association. The induced conformational change generates a local enrichment of the protein in the so-called dehydrons, which are solvent-exposed backbone hydrogen bonds. These structural deficiencies that are inherently sticky are not present in the apo form and constitute a local metastable state that promotes association with the ligand. This productive induced metastability (PIM) is likely to translate into a general molecular design concept.