Interfacial rheology of lanthanide binding peptide surfactants at the air-water interface.

Peptide surfactants (PEPS) are studied to capture and retain rare earth elements (REEs) at air-water interfaces to enable REE separations. Peptide sequences, designed to selectively bind REEs, depend crucially on the position of ligands within their binding loop domain. These ligands form a coordination sphere that wraps and retains the cation. We study variants of lanthanide binding tags (LBTs) designed to complex strongly with Tb3+. The peptide LBT5- (with net charge -5) is known to bind Tb3+ and adsorb with more REE cations than peptide molecules, suggesting that undesired non-specific coulombic interactions occur. Rheological characterization of interfaces of LBT5- and Tb3+ solutions reveal the formation of an interfacial gel. To probe whether this gelation reflects chelation among intact adsorbed LBT5-:Tb3+ complexes or destruction of the binding loop, we study a variant, LBT3-, designed to form net neutral LBT3-:Tb3+ complexes. Solutions of LBT3- and Tb3+ form purely viscous layers in the presence of excess Tb3+, indicating that each peptide binds a single REE in an intact coordination sphere. We introduce the variant RR-LBT3- with net charge -3 and anionic ligands outside of the coordination sphere. We find that such exposed ligands promote interfacial gelation. Thus, a nuanced requirement for interfacial selectivity of PEPS is proposed: that anionic ligands outside of the coordination sphere must be avoided to prevent the non-selective recruitment of REE cations. This view is supported by simulation, including interfacial molecular dynamics simulations, and interfacial metadynamics simulations of the free energy landscape of the binding loop conformational space.

Reexpansion of charged nanoparticle assemblies in concentrated electrolytes.

Electrostatic forces in solutions are highly relevant to a variety of fields, ranging from electrochemical energy storage to biology. However, their manifestation in concentrated electrolytes is not fully understood, as exemplified by counterintuitive observations of colloidal stability and long-ranged repulsions in molten salts. Highly charged biomolecules, such as DNA, respond sensitively to ions in dilute solutions. Here, we use non-base-pairing DNA-coated nanoparticles (DNA-NP) to analyze electrostatic interactions in concentrated salt solutions. Despite their negative charge, these conjugates form colloidal crystals in solutions of sufficient divalent cation concentration. We utilize small-angle X-ray scattering (SAXS) to study such DNA-NP assemblies across the full accessible concentration ranges of aqueous CaCl2, MgCl2, and SrCl2 solutions. SAXS shows that the crystallinity and phases of the assembled structures vary with cation type. For all tested salts, the aggregates contract with added ions at low salinities and then begin expanding above a cation-dependent threshold salt concentration. Wide-angle X-ray scattering (WAXS) reveals enhanced positional correlations between ions in the solution at high salt concentrations. Complementary molecular dynamics simulations show that these ion-ion interactions reduce the favorability of dense ion configurations within the DNA brushes below that of the bulk solution. Measurements in solutions with lowered permittivity demonstrate a simultaneous increase in ion coupling and decrease in the concentration at which aggregate expansion begins, thus confirming the connection between these phenomena. Our work demonstrates that interactions between charged objects continue to evolve considerably into the high-concentration regime, where classical theories project electrostatics to be of negligible consequence.

Transformations in crystals of DNA-functionalized nanoparticles by electrolytes.

Colloidal crystals have applications in water treatments, including water purification and desalination technologies. It is, therefore, important to understand the interactions between colloids as a function of electrolyte concentration. We study the assembly of DNA-grafted gold nanoparticles immersed in concentrated electrolyte solutions. Increasing the concentration of divalent Ca2+ ions leads to the condensation of nanoparticles into face-centered-cubic (FCC) crystals at low electrolyte concentrations. As the electrolyte concentration increases, the system undergoes a phase change to body-centered-cubic (BCC) crystals. This phase change occurs as the interparticle distance decreases. Molecular dynamics analysis suggests that the interparticle interactions change from strongly repulsive to short-range attractive as the divalent-electrolyte concentration increases. A thermodynamic analysis suggests that increasing the salt concentration leads to significant dehydration of the nanoparticle environment. We conjecture that the intercolloid attractive interactions and dehydrated states favour the BCC structure. Our results gain insight into salting out of colloids such as proteins as the concentration of salt increases in the solution.

Surface polarization enhances ionic transport and correlations in electrolyte solutions nanoconfined by conductors.

Layered materials that perform mixed electron and ion transport are promising for energy harvesting, water desalination, and bioinspired functionalities. These functionalities depend on the interaction between ionic and electronic charges on the surface of materials. Here we investigate ion transport by an external electric field in an electrolyte solution confined in slit-like channels formed by two surfaces separated by distances that fit only a few water layers. We study different electrolyte solutions containing monovalent, divalent, and trivalent cations, and we consider walls made of non-polarizable surfaces and conductors. We show that considering the surface polarization of the confining surfaces can result in a significant increase in ionic conduction. The ionic conductivity is increased because the conductors' screening of electrostatic interactions enhances ionic correlations, leading to faster collective transport within the slit. While important, the change in water's dielectric constant in confinement is not enough to explain the enhancement of ion transport in polarizable slit-like channels.

Microporous water with high gas solubilities.

Liquids with permanent microporosity can absorb larger quantities of gas molecules than conventional solvents1, providing new opportunities for liquid-phase gas storage, transport and reactivity. Current approaches to designing porous liquids rely on sterically bulky solvent molecules or surface ligands and, thus, are not amenable to many important solvents, including water2-4. Here we report a generalizable thermodynamic strategy to preserve permanent microporosity and impart high gas solubilities to liquid water. Specifically, we show how the external and internal surface chemistry of microporous zeolite and metal-organic framework (MOF) nanocrystals can be tailored to promote the formation of stable dispersions in water while maintaining dry networks of micropores that are accessible to gas molecules. As a result of their permanent microporosity, these aqueous fluids can concentrate gases, including oxygen (O2) and carbon dioxide (CO2), to much higher densities than are found in typical aqueous environments. When these fluids are oxygenated, record-high capacities of O2 can be delivered to hypoxic red blood cells, highlighting one potential application of this new class of microporous liquids for physiological gas transport.

Probing the size-dependent polarizability of mesoscopic ionic clusters and their induced-dipole interactions.

Mesoscopic clusters composed of oppositely charged particles are ubiquitous in synthetic and biological soft materials. The effective interaction between these clusters is influenced by their polarizability, that is, the ability of their constituent charges to re-arrange in response to an external electrical field. Here, using coarse-grained simulations, we show that the polarizability of electrically neutral ionic clusters decreases as the number of constituent charges increases and/or their Coulombic interaction strength increases for various ion valencies, ion densities, and degrees of cluster boundary hardness. For clusters of random ionomers and their counterions, their polarizability is shown to depend on the number of polymer chains. The variation of the cluster polarizability with the cluster size indicates that throughout the assembly, the induced-dipole interactions between the clusters may be reduced substantially as they acquire more charges while maintaining zero net charge. Under certain conditions, the induced-dipole interactions may become repulsive, as inferred from our simulations with a polarizable solvent. As a result, the dipole-induced related interactions can serve as a counterbalancing force that contributes to the self-limiting aggregation of charge-containing assemblies.

A Modeling-Based Design to Engineering Protein Hydrogels with Random Copolymers.

Protein enzymes have shown great potential in numerous technological applications. However, the design of supporting materials is needed to preserve protein functionality outside their native environment. Direct enzyme-polymer self-assembly offers a promising alternative to immobilize proteins in an aqueous solution, achieving higher control of their stability and enzymatic activity in industrial applications. Herein, we propose a modeling-based design to engineering hydrogels of cytochrome P450 and of PETase with styrene/2-vinylpyridine (2VP) random copolymers. By tuning the copolymer fraction of polar groups and of charged groups via quaternization of 2VP for coassembly with cytochrome P450 and via sulfonation of styrene for coassembly with PETase, we provide quantitative guidelines to select either a protein-polymer hydrogel structure or a single-protein encapsulation. The results highlight that, regardless of the protein surface domains, the presence of polar interactions and hydration effects promote the formation of a more elongated enzyme-polymer complex, suggesting a membrane-like coassembly. On the other hand, the effectiveness of a single-protein encapsulation is reached by decreasing the fraction of polar groups and by increasing the charge fraction up to 15%. Our computational analysis demonstrates that the enzyme-polymer assemblies are first promoted by the hydrophobic interactions which lead the protein nonpolar residues to achieve the maximum coverage and to play the role of the most robust contact points. The mechanisms of coassembly are unveiled in the light of both protein and polymer physical-chemistry, providing bioconjugate phase diagrams for the optimal material design.

Insights into the Enhanced Catalytic Activity of Cytochrome c When Encapsulated in a Metal-Organic Framework.

The encapsulation of enzymes within porous materials has shown great promise, not only in protecting the enzymes from denaturation under nonbiological environments, but also, in some cases, in facilitating their enzymatic reaction rates at favorable reaction conditions. While a number of hypotheses have been developed to explain this phenomenon, the detailed structural changes of the enzymes upon encapsulation within the porous material, which are closely related to their activity, remain largely elusive. Herein, the structural change of cytochrome c (Cyt c) upon encapsulation within a hierarchical metal-organic framework, NU-1000, is investigated through a combination of experimental and computational methods, such as electron paramagnetic resonance, solid-state ultraviolet-visible spectroscopy, and all-atom explicit solvent molecular dynamics simulations. The enhanced catalytic performance of Cyt c after being encapsulated within NU-1000 is supported by the physical and in silico observations of a change around the heme ferric active center.

Water follows polar and nonpolar protein surface domains.

The conformation of water around proteins is of paramount importance, as it determines protein interactions. Although the average water properties around the surface of proteins have been provided experimentally and computationally, protein surfaces are highly heterogeneous. Therefore, it is crucial to determine the correlations of water to the local distributions of polar and nonpolar protein surface domains to understand functions such as aggregation, mutations, and delivery. By using atomistic simulations, we investigate the orientation and dynamics of water molecules next to 4 types of protein surface domains: negatively charged, positively charged, and charge-neutral polar and nonpolar amino acids. The negatively charged amino acids orient around 98% of the neighboring water dipoles toward the protein surface, and such correlation persists up to around 16 Å from the protein surface. The positively charged amino acids orient around 94% of the nearest water dipoles against the protein surface, and the correlation persists up to around 12 Å. The charge-neutral polar and nonpolar amino acids are also orienting the water neighbors in a quantitatively weaker manner. A similar trend was observed in the residence time of the nearest water neighbors. These findings hold true for 3 technically important enzymes (PETase, cytochrome P450, and organophosphorus hydrolase). Our results demonstrate that the water-amino acid degree of correlation follows the same trend as the amino acid contribution in proteins solubility, namely, the negatively charged amino acids are the most beneficial for protein solubility, then the positively charged amino acids, and finally the charge-neutral amino acids.

Self-Assembly of Charge-Containing Copolymers at the Liquid-Liquid Interface.

Quantitatively understanding the self-assembly of amphiphilic macromolecules at liquid-liquid interfaces is a fundamental scientific concern due to its relevance to a broad range of applications including bottom-up nanopatterning, protein encapsulation, oil recovery, drug delivery, and other technologies. Elucidating the mechanisms that drive assembly of amphiphilic macromolecules at liquid-liquid interfaces is challenging due to the combination of hydrophobic, hydrophilic, and Coulomb interactions, which require consideration of the dielectric mismatch, solvation effects, ionic correlations, and entropic factors. Here we investigate the self-assembly of a model block copolymer with various charge fractions at the chloroform-water interface. We analyze the adsorption and conformation of poly(styrene)-block-poly(2-vinylpyridine) (PS-b-P2VP) and of the homopolymer poly(2-vinylpyridine) (P2VP) with varying charge fraction, which is controlled via a quaternization reaction and distributed randomly along the backbone. Interfacial tension measurements show that the polymer adsorption increases only marginally at low charge fractions (<5%) but increases more significantly at higher charge fractions for the copolymer, while the corresponding randomly charged P2VP homopolymer analogues display much more sensitivity to the presence of charged groups. Molecular dynamics (MD) simulations of the experimental systems reveal that the diblock copolymer (PS-b-P2VP) interfacial activity could be mediated by the formation of a rich set of complex interfacial copolymer aggregates. Circular domains to elongated stripes are observed in the simulations at the water-chloroform interface as the charge fraction increases. These structures are shown to resemble the spherical and cylindrical helicoid structures observed in bulk chloroform as the charge fraction increases. The self-assembly of charge-containing copolymers is found to be driven by the association of the charged component in the hydrophilic block, with the hydrophobic segments extending away from the hydrophilic cores into the chloroform phase.

Hydrophobic Hydration and the Effect of NaCl Salt in the Adsorption of Hydrocarbons and Surfactants on Clathrate Hydrates.

Adsorption of functional molecules on the surface of hydrates is key in the understanding of hydrate inhibitors. We investigate the adsorption of a hydrocarbon chain, nonionic and ionic surfactants, and ions at the hydrate-aqueous interface. Our results suggest a strong connection between the water ordering around solutes in bulk and the affinity for the hydrates surface. We distinguish two types of water ordering around solutes: (i) hydrophobic hydration where water molecules form a hydrogen bond network similar to clathrate hydrates, and (ii) ionic hydration where water molecules align according to the polarity of an ionic group. The nonionic surfactant and the hydrocarbon chain induce hydrophobic hydration and are favorably adsorbed on the hydrate surface. Adsorption of ions and the ionic headgroups on the hydrate surface is not favorable because ionic hydration and the hydrogen bond structure of hydrates are incompatible. The nonionic surfactant is adsorbed by the headgroup and tail while adsorption of the ionic surfactants is not favorable through the head. Water ordering is analyzed using the hydrogen bond and tetrahedral density profiles as a function of the distance to the chemical groups. The adsorption of solutes is studied through the free energy profiles as a function of the distance to the hydrate surface. Salt lowers the melting temperature of hydrates, disrupts hydrophobic hydration, reduces the solubility of solutes in the aqueous solution, and increases the propensity of solutes to be adsorbed on hydrate surfaces. Our studies are performed by the unbiased and steered molecular dynamics simulations. The results are in line with experiments on the effect of salt and alkanes in hydrate antiagglomeration.

Enhanced Hydrate Nucleation Near the Limit of Stability.

Clathrate hydrates are crystalline structures composed of small guest molecules trapped into cages formed by hydrogen-bonded water molecules. In hydrate nucleation, water and the guest molecules may stay in a metastable fluid mixture for a long period. Metastability is broken if the concentration of the guest is above certain limit. We perform molecular dynamics (MD) simulations of supersaturated water-propane solutions close to the limit of stability. We show that hydrate nucleation can be very fast in a very narrow range of composition at moderate temperatures. Propane density fluctuations near the fluid-fluid demixing are coupled with crystallization producing en- hanced nucleation rates. This is the first report of propane-hydrate nucleation by MD simulations. We observe motifs of the crystalline structure II in line with experiments and new hydrate cages not reported in the literature. Our study relates nucleation to the fluid-fluid spinodal decomposition and demonstration that the enhanced nucleation phenomenon is more general than short range attractive interactions as suggested in nucleation of proteins.

Effects of mixed discrete surface charges on the electrical double layer.

Adsorption of surface coions and charge reversal are induced at the electrical double layer of a wall charged with positive and negative surface sites next to an electrolyte solution. While for the considered surface charge density these effects are found over a wide range of conditions, they are not observed for the typically employed surface models in equivalent conditions. Important consequences in electrophoresis experiments for different colloids with equal effective surface charge density are foreseen. This study is carried out by means of molecular dynamics simulations.

Polarity inversion of ζ-potential in concentrated colloidal dispersions.

A concentrated colloidal dispersion is studied by applying an integral equations theory to the colloidal primitive model fluid. Important effects, attributed to large size and charge and to the finite concentration of colloidal particles, are found. We observe a polarity inversion of ζ-potential for concentrated colloidal dispersions, while it is not present for a single colloidal particle at infinite dilution. An excellent qualitative agreement between our theoretical predictions and our computer simulations is observed.

Stability mechanisms for plate-like nanoparticles immersed in a macroion dispersion.

An integral equation theory and Monte Carlo simulations are applied to study a model macroion solution confined between two parallel plates immersed in a 1:1 electrolyte and the macroions' counterions. We analyze the cases in which plates are: (a) uncharged; (b) when they are like-charged to the macroions; (c) when they are oppositely charged to the macroions. For all cases a long range oscillatory behavior of the induced charge density between the plates is found (implying an overcompensation/undercompensation of the plates' charge density) and a correlation between the confined and outside fluids. The behavior of the force is discussed in terms of the macroion and ion structure inside and outside the plates. A good agreement is found between theoretical and simulation results.

Electrokinetic properties of monovalent electrolytes confined in charged nanopores: effect of geometry and ionic short-range correlations.

The electrokinetic properties (such as capillary conductance, electroviscosity, and the streaming potential) are obtained for a restricted primitive model electrolyte confined in a slitlike nanopore made up of two infinite parallel plates and in a cylindrical cavity of infinite extension. The hypernetted chain/mean spherical approximation (HNC/MSA) is used to obtain the equilibrium ionic concentration profiles inside the pores, which in turn are used to calculate the electrokinetic properties via linear hydrodynamic equations. Our results are compared with those obtained via the classical Poisson-Boltzmann (PB) theory. Important quantitative and qualitative effects, attributed to geometry and to the proper consideration of short-range correlations by HNC/MSA, are discussed.

On the regimes of charge reversal.

Charge reversal of the planar electrical double layer is studied by means of a well known integral equation theory. By a numerical analysis, a diagram is constructed with the onset points of charge reversal in the space of the fundamental variables of the system. Within this diagram, two regimes of charge reversal are identified, which are referred to as oscillatory and nonoscillatory. We found that these two regimes can be distinguished through a simple formula. Furthermore, a symmetry between electrostatic and size correlations in charge reversal is exhibited. Agreement of our results with other theories and molecular simulations data is discussed.