Orientational Pathways during Protein Translocation through Polymer-Modified Nanopores.

Protein translocation through nanopores holds significant promise for applications in biotechnology, biomolecular analysis, and medicine. However, the interpretation of signals generated by the translocation of the protein remains challenging. In this way, it is crucial to gain a comprehensive understanding on how macromolecules translocate through a nanopore and to identify what are the critical parameters that govern the process. In this study, we investigate the interplay between protein charge regulation, orientation, and nanopore surface modifications using a theoretical framework that allows us to explicitly take into account the acid-base reactions of the titrable amino acids in the proteins and in the polyelectrolytes grafted to the nanopore surface. Our goal is to thoroughly characterize the translocation process of different proteins (GFP, ß-lactoglobulin, lysozyme, and RNase) through nanopores modified with weak polyacids. Our calculations show that the charge regulation mechanism exerts a profound effect on the translocation process. The pH-dependent interactions between proteins and charged polymers within the nanopore lead to diverse free energy landscapes with barriers, wells, and flat regions dictating translocation efficiency. Comparison of different proteins allows us to identify the significance of protein isoelectric point, size, and morphology in the translocation behavior. Taking advantage of these insights, we propose pH-responsive nanopores that can load proteins at one pH and release them at another, offering opportunities for controlled protein delivery, separation, and sensing applications.

Theoretical treatment of complex coacervate core micelles: structure and pH-induced disassembly.

Complex coacervate core micelles (C3Ms) are supramolecular soft nanostructures formed by the assembly of a block copolymer and an oppositely charged homopolymer. The coacervation of the charged segments in both macromolecules drives the formation of the core of the C3M, while the neutral block of the copolymer forms the corona. This work introduces a molecular theory (MOLT) that predicts the internal structure and stimuli-responsive properties of C3Ms and explicitly considers the chemical architecture of the polyelectrolytes, their acid-based equilibria and electrostatic and non-electrostatic interactions. In order to accurately predict complex coacervation, the correlations between charged species are incorporated into MOLT as ion-pairing processes, which are modeled using a coupled chemical equilibrium formalism. Very good agreement was observed between the experimental results in the literature and MOLT predictions for the scaling relationships that relate the dimensions of the micelle (aggregation number and sizes of the micelle and the core) to the lengths of the different blocks. MOLT was used to study the disassembly of the micelles when the solution pH is driven away from the value that guarantees the charge stoichiometry of the core. This study reveals that very sharp disassembly transitions can be obtained by tuning the length or architecture of the copolymer component, thereby suggesting potential routes to design C3Ms capable of releasing their components at very precise pH values.

Water Self-Dissociation is Insensitive to Nanoscale Environments.

Nanoconfinement effects on water dissociation and reactivity remain controversial, despite their importance to understand the aqueous chemistry at interfaces, pores, or aerosols. The pKw in confined environments has been assessed from experiments and simulations in a few specific cases, leading to dissimilar conclusions. Here, with the use of carefully designed ab initio simulations, we demonstrate that the energetics of bulk water dissociation is conserved intact to unexpectedly small length-scales, down to aggregates of only a dozen molecules or pores of widths below 2 nm. The reason is that most of the free-energy involved in water autoionization comes from breaking the O-H covalent bond, which has a comparable barrier in the bulk liquid, in a small droplet of nanometer size, or in a nanopore in the absence of strong interfacial interactions. Thus, dissociation free-energy profiles in nanoscopic aggregates or in 2D slabs of 1 nm width reproduce the behavior corresponding to the bulk liquid, regardless of whether the corresponding nanophase is delimited by a solid or a gas interface. The present work provides a definite and fundamental description of the mechanism and thermodynamics of water dissociation at different scales with broader implications on reactivity and self-ionization at the air-liquid interface.

Revisiting the Mechanisms of Charge Transport in Solutions of Redox-Active Molecules Using Computer Simulations: When and Why Do Analytical Theories Fail?

Understanding charge transport is essential for the development of energy-storage applications. This work introduces a new theoretical methodology to model diffusive charge transport in solutions of redox-active molecules by combining Langevin dynamics for the spatial degrees of freedom and a master-equation formalism to describe the electron-hopping events between redox molecules. The model is used to analyze the effects of the concentration of the redox molecules and the strength of the intermolecular interactions on the charge-transport mechanism. In the past, the rate of charge transport has been modeled with the analytical Dahms-Ruff equation; however, this is a mean-field equation, whose range of validity has not been tested with less approximate theories. We show that the Dahms-Ruff equation fails to quantitatively predict the diffusion coefficient for charge transport for large concentrations of the redox species and high bimolecular electron-transfer rates, i.e., the most relevant conditions for energy-storage applications. Under these conditions, the diffusion coefficient for charge transport obtained from simulations is larger than that predicted from the Dahm-Ruff equation because of the formation of transient clusters of redox molecules. Also, intermolecular interactions, which are not taken into account by the Dahms-Ruff equation, play a central role in the charge transport of redox species. We show that the apparent diffusion coefficient experiences a maximum with respect to the strength of the intermolecular attractions. This maximum is traced back to the formation of clusters and their two opposite effects on the diffusion coefficient: electron hopping is fast within a cluster but inefficient between neighboring clusters.

Molecular Theory: A Tool for Predicting the Outcome of Self-Assembly of Polymers, Nanoparticles, Amphiphiles, and Other Soft Materials.

The supramolecular organization of soft materials, such as colloids, polymers, and amphiphiles, results from a subtle balance of weak intermolecular interactions and entropic forces. This competition can drive the self-organization of soft materials at the nano-/mesoscale. Modeling soft-matter self-assembly requires, therefore, considering a complex interplay of forces at the relevant length scales without sacrificing the molecular details that define the chemical identity of the system. This mini-review focuses on the application of a tool known as molecular theory to study self-assembly in different types of soft materials. This tool is based on extremizing an approximate free energy functional of the system, and, therefore, it provides a direct, computationally affordable estimation of the stability of different self-assembled morphologies. Moreover, the molecular theory explicitly incorporates structural details of the chemical species in the system, accounts for their conformational degrees of freedom, and explicitly includes their chemical equilibria. This mini-review introduces the general ideas behind the theoretical formalism and discusses its advantages and limitations compared with other theoretical tools commonly used to study self-assembled soft materials. Recent application examples are discussed: the self-patterning of polyelectrolyte brushes on planar and curved surfaces, the formation of nanoparticle (NP) superlattices, and the self-organization of amphiphiles into micelles of different shapes. Finally, prospective methodological improvements and extensions (also relevant for related theoretical tools) are analyzed.

Mechanisms of Enzymatic Transduction in Nanochannel Biosensors.

The immobilization of enzymes in solid-state nanochannels is a new avenue for the design of biosensors with outstanding selectivity and sensitivity. This work reports the first theoretical model of an enzymatic nanochannel biosensor. The model is applied to the system previously experimentally studied by Lin, et al. (Anal. Chem. 2014, 86, 10546): a hourglass nanochannel modified by glucose oxidase for the detection of glucose. Our predictions are in good agreement with experimental observations as a function of the applied potential, pH and glucose concentration. The sensing mechanism results from the combination of three processes: i) the establishment of a steady-state proton concentration gradient due to a reaction-diffusion mechanism, ii) the effect of that gradient on the charge of the adsorbed enzymes and native surface groups, and iii) the effect of the resulting surface charge on the ionic current. Strategies to improve the sensor performance based on this mechanism are identified and discussed.

Nanopore gates via reversible crosslinking of polymer brushes: a theoretical study.

Polymer-brush-modified nanopores are synthetic structures inspired by the gated transport exhibited by their biological counterparts. This work theoretically analyzes how the reversible crosslinking of a polymer network by soluble species can be used to control transport through nanochannels and pores. The study was performed with a molecular theory that allows inhomogeneities in the three spatial dimensions and explicitly takes into account the size, shape and conformations of all molecular species, considers the intermolecular interactions between the polymers and the soluble crosslinkers and includes the presence of a translocating particle inside the pore. It is shown than increasing the concentration of the soluble crosslinkers in bulk solution leads to a gradual increase of its number within the pore until a critical bulk concentration is reached. At the critical concentration, the number of crosslinkers inside the pore increases abruptly. For long chains, this sudden transition triggers the collapse of the polymer brush to the center of the nanopore. The resulting structure increases the free-energy barrier that a translocating particle has to surmount to go across the pore and modifies the route of translocation from the axis of the pore to its walls. On the other hand, for short polymer chains the crosslinkers trigger the collapse of the brush to the pore walls, which reduces the translocation barrier.

Electrochemically Generated Nanobubbles: Invariance of the Current with Respect to Electrode Size and Potential.

Gas-producing electrochemical reactions are key to energy conversion and generation technologies. Bubble formation dramatically decreases gas-production rates on nanoelectrodes, by confining the reaction to the electrode boundary. This results in the collapse of the current to a stationary value independent of the potential. Startlingly, these residual currents also appear to be insensitive to the nanoelectrode diameter in the 5 to 500 nm range. These results are counterintuitive, as it may be expected that the current be proportional to the circumference of the electrode, i.e., the length of the three-phase line where the reaction occurs. Here, we use molecular simulations and a kinetic model to elucidate the origin of current insensitivity with respect to the potential and establish its relationship to the size of nanoelectrodes. We provide critical insights for the design and operation of nanoscale electrochemical devices and demonstrate that nanoelectrode arrays maximize conversion rates compared to macroscopic electrodes with same total area.

Mechanisms of Nucleation and Stationary States of Electrochemically Generated Nanobubbles.

Gas evolving reactions are ubiquitous in the operation of electrochemical devices. Recent studies of individual gas bubbles on nanoelectrodes have resulted in unprecedented control and insights on their formation. The experiments, however, lack the spatial resolution to elucidate the molecular pathway of nucleation of nanobubbles and their stationary size and shape. Here we use molecular simulations with an algorithm that mimics the electrochemical formation of gas, to investigate the mechanisms of nucleation of gas bubbles on nanoelectrodes, and characterize their stationary states. The simulations reproduce the experimental currents in the induction and stationary stages, and indicate that surface nanobubbles nucleate through a classical mechanism. We identify three distinct regimes for bubble nucleation, depending on the binding free energy per area of bubble to the electrode, Δγbind. If Δγbind is negative, the nucleation is heterogeneous and the nanobubble remains bound to the electrode, resulting in a low-current stationary state. For very negative Δγ, the bubble fully wets the electrode, forming a one-layer-thick micropancake that nucleates without supersaturation. On the other hand, when Δγbind > 0 the nanobubble nucleates homogeneously close to the electrode, but never attaches to it. We conclude that all surface nanobubbles must nucleate heterogeneously. The simulations reveal that the size and contact angle of stationary nanobubbles increase with the reaction driving force, although their residual current is invariant. The myriad of driven nonequilibrium stationary states with the same rate of production of gas, but distinct bubble properties, suggests that these dissipative systems have attractors that control the stationary current.

Pore condensation and freezing is responsible for ice formation below water saturation for porous particles.

Ice nucleation in the atmosphere influences cloud properties, altering precipitation and the radiative balance, ultimately regulating Earth's climate. An accepted ice nucleation pathway, known as deposition nucleation, assumes a direct transition of water from the vapor to the ice phase, without an intermediate liquid phase. However, studies have shown that nucleation occurs through a liquid phase in porous particles with narrow cracks or surface imperfections where the condensation of liquid below water saturation can occur, questioning the validity of deposition nucleation. We show that deposition nucleation cannot explain the strongly enhanced ice nucleation efficiency of porous compared with nonporous particles at temperatures below -40 °C and the absence of ice nucleation below water saturation at -35 °C. Using classical nucleation theory (CNT) and molecular dynamics simulations (MDS), we show that a network of closely spaced pores is necessary to overcome the barrier for macroscopic ice-crystal growth from narrow cylindrical pores. In the absence of pores, CNT predicts that the nucleation barrier is insurmountable, consistent with the absence of ice formation in MDS. Our results confirm that pore condensation and freezing (PCF), i.e., a mechanism of ice formation that proceeds via liquid water condensation in pores, is a dominant pathway for atmospheric ice nucleation below water saturation. We conclude that the ice nucleation activity of particles in the cirrus regime is determined by the porosity and wettability of pores. PCF represents a mechanism by which porous particles like dust could impact cloud radiative forcing and, thus, the climate via ice cloud formation.

Stability and Vapor Pressure of Aqueous Aggregates and Aerosols Containing a Monovalent Ion.

The incidence of charged particles on the nucleation and the stability of aqueous aggregates and aerosols was reported more than a century ago. Many studies have been conducted ever since to characterize the stability, structure, and nucleation barrier of ion-water droplets. Most of these studies have focused on the free-energy surface as a function of cluster size, with an emphasis on the role of ionic charge and radius. This knowledge is fundamental to go beyond the rudimentary ion-induced classical nucleation theory. In the present article, we address this problem from a different perspective, by computing the vapor pressures of (H2O)nLi+ and (H2O)nCl- aggregates using molecular simulations. Our calculations shed light on the structure, the critical size, the range of stability, and the role of ion-water interactions in aqueous clusters. Moreover, they allow one to assess the accuracy of the classical thermodynamic model, highlighting its strengths and weaknesses.

High-Resolution Coarse-Grained Model of Hydrated Anion-Exchange Membranes that Accounts for Hydrophobic and Ionic Interactions through Short-Ranged Potentials.

Molecular simulations provide a versatile tool to study the structure, anion conductivity, and stability of anion-exchange membrane (AEM) materials and can provide a fundamental understanding of the relation between structure and property of membranes that is key for their use in fuel cells and other applications. The quest for large spatial and temporal scales required to model the multiscale structure and transport processes in the polymer electrolyte membranes, however, cannot be met with fully atomistic models, and the available coarse-grained (CG) models suffer from several challenges associated with their low-resolution. Here, we develop a high-resolution CG force field for hydrated polyphenylene oxide/trimethylamine chloride (PPO/TMACl) membranes compatible with the mW water model using a hierarchical parametrization approach based on Uncertainty Quantification and reference atomistic simulations modeled with the Generalized Amber Force Field (GAFF) and TIP4P/2005 water. The parametrization weighs multiple properties, including coordination numbers, radial distribution functions (RDFs), self-diffusion coefficients of water and ions, relative vapor pressure of water in the solution, hydration enthalpy of the tetramethylammonium chloride (TMACl) salt, and cohesive energy of its aqueous solutions. We analyze the interdependence between properties and address how to compromise between the accuracies of the properties to achieve an overall best representability. Our optimized CG model FFcomp quantitatively reproduces the diffusivities and RDFs of the reference atomistic model and qualitatively reproduces the experimental relative vapor pressure of water in solutions of tetramethylammonium chloride. These properties are of utmost relevance for the design and operation of fuel cell membranes. To our knowledge, this is the first CG model that includes explicitly each water and ion and accounts for hydrophobic, ionic, and intramolecular interactions explicitly parametrized to reproduce multiple properties of interest for hydrated polyelectrolyte membranes. The CG model of hydrated PPO/TMACl water is about 100 times faster than the reference atomistic GAFF-TIP4P/2005 model. The strategy implemented here can be used in the parametrization of CG models for other substances, such as biomolecular systems and membranes for desalination, water purification, and redox flow batteries. We anticipate that the large spatial and temporal simulations made possible by the CG model will advance the quest for anion-exchange membranes with improved transport and mechanical properties.

Vapor Pressure of Aqueous Solutions of Electrolytes Reproduced with Coarse-Grained Models without Electrostatics.

The vapor pressure of water is a key property in a large class of applications from the design of membranes for fuel cells and separations to the prediction of the mixing state of atmospheric aerosols. Molecular simulations have been used to compute vapor pressures, and a few studies on liquid mixtures and solutions have been reported on the basis of the Gibbs Ensemble Monte Carlo method in combination with atomistic force fields. These simulations are costly, making them impractical for the prediction of the vapor pressure of complex materials. The goal of the present work is twofold: (1) to demonstrate the use of the grand canonical screening approach ( Factorovich , M. H. J. Chem. Phys. 2014 , 140 , 064111 ) to compute the vapor pressure of solutions and to extend the methodology for the treatment of systems without a liquid-vapor interface and (2) to investigate the ability of computationally efficient high-resolution coarse-grained models based on the mW monatomic water potential and ions described exclusively with short-range interactions to reproduce the relative vapor pressure of aqueous solutions. We find that coarse-grained models of LiCl and NaCl solutions faithfully reproduce the experimental relative pressures up to high salt concentrations, despite the inability of these models to predict cohesive energies of the solutions or the salts. A thermodynamic analysis reveals that the coarse-grained models achieve the experimental activity coefficients of water in solution through a compensation of severely underestimated hydration and vaporization free energies of the salts. Our results suggest that coarse-grained models developed to replicate the hydration structure and the effective ion-ion attraction in solution may lead to this compensation. Moreover, they suggest an avenue for the design of coarse-grained models that accurately reproduce the activity coefficients of solutions.