Ionic size effects on the Poisson-Boltzmann theory

Adsorption of polyelectrolytes in the presence of varying dielectric discontinuity between solution and substrate

We examine the interactions between polyelectrolytes (PEs) and uncharged substrates under conditions corresponding to a dielectric discontinuity between the aqueous solution and the substrate. To this end, we vary the relevant system characteristics, in particular the substrate dielectric constant ɛs under different salt conditions. We employ coarse-grained molecular dynamics simulations with rodlike PEs in salt solutions with explicit ions and implicit water solvent with dielectric constant ɛw = 80. As expected, at low salt concentrations, PEs are repelled from the substrates with ɛs < ɛw but are attracted to substrates with a high dielectric constant due to image charges. This attraction considerably weakens for high salt and multivalent counterions due to enhanced screening. Furthermore, for monovalent salt, screening enhances adsorption for weakly charged PEs, but weakens it for strongly charged ones. Meanwhile, multivalent counterions have little effect on weakly charged PEs, but prevent adsorption of highly charged PEs, even at low salt concentrations. We also find that correlation-induced charge inversion of a PE is enhanced close to the low dielectric constant substrates, but suppressed when the dielectric constant is high. To explore the possibility of a PE monolayer formation, we examine the interaction of a pair of like-charged PEs aligned parallel to a high dielectric constant substrate with ɛs = 8000. Our main conclusion is that monolayer formation is possible only for weakly charged PEs at high salt concentrations of both monovalent and multivalent counterions. Finally, we also consider the energetics of a PE approaching the substrate perpendicular to it, in analogy to polymer translocation. Our results highlight the complex interplay between electrostatic and steric interactions and contribute to a deeper understanding of PE-substrate interactions and adsorption at substrate interfaces with varying dielectric discontinuities from solution, ubiquitous in biointerfaces, PE coating applications, and designing adsorption setups.

Interactions between Rigid Polyelectrolytes Mediated by Ordering and Orientation of Multivalent Nonspherical Ions in Salt Solutions

Multivalent ions in solutions with polyelectrolytes (PEs) induce electrostatic correlations that can drastically change ion distributions around the PEs and their mutual interactions. Using coarse-grained molecular dynamics simulations, we show how in addition to valency, ion shape and concentration can be harnessed as tools to control rigid like-charged PE-PE interactions. We demonstrate a correlation between the orientational ordering of aspherical ions and how they mediate the effective PE-PE attraction induced by multivalency. The interaction type, strength, and range can thus be externally controlled in ionic solutions. Our results can be used as generic guidelines to tune the self-assembly of like-charged polyelectrolytes by variation of the characteristics of the ions.

Cumulative Millisecond-Long Sampling for a Comprehensive Energetic Evaluation of Aqueous Ionic Liquid Effects on Amino Acid Interactions

The interactions of amino acid side-chains confer diverse energetic contributions and physical properties to a protein's stability and function. Various computational tools estimate the effect of changing a given amino acid on the protein's stability based on parametrized (free) energy functions. When parametrized for the prediction of protein stability in water, such energy functions can lead to suboptimal results for other solvents, such as ionic liquids (IL), aqueous ionic liquids (aIL), or salt solutions. However, to our knowledge, no comprehensive data are available describing the energetic effects of aIL on intramolecular protein interactions. Here, we present the most comprehensive set of potential of mean force (PMF) profiles of pairwise protein-residue interactions to date, covering 50 relevant interactions in water, the two biotechnologically relevant aIL [BMIM/Cl] and [BMIM/TfO], and [Na/Cl]. These results are based on a cumulated simulation time of >1 ms. aIL and salt ions can weaken, but also strengthen, specific residue interactions by more than 3 kcal mol^-1, depending on the residue pair, residue-residue configuration, participating ions, and concentration, necessitating considering such interactions specifically. These changes originate from a complex interplay of competitive or cooperative noncovalent ion-residue interactions, changes in solvent structural dynamics, or unspecific charge screening effects and occur at the contact distance but also at larger, solvent-separated distances. This data provide explanations at the atomistic and energetic levels for complex IL effects on protein stability and should help improve the prediction accuracies of computational tools that estimate protein stability based on (free) energy functions.

Modified Poisson–Boltzmann theory for polyelectrolytes in monovalent salt solutions with finite-size ions

We present a soft-potential-enhanced Poisson–Boltzmann (SPB) theory to efficiently capture ion distributions and electrostatic potential around rodlike charged macromolecules. The SPB model is calibrated with a coarse-grained particle-based model for polyelectrolytes (PEs) in monovalent salt solutions as well as compared to a full atomistic molecular dynamics simulation with the explicit solvent. We demonstrate that our modification enables the SPB theory to accurately predict monovalent ion distributions around a rodlike PE in a wide range of ion and charge distribution conditions in the weak-coupling regime. These include excess salt concentrations up to 1M and ion sizes ranging from small ions, such as Na+ or Cl−, to softer and larger ions with a size comparable to the PE diameter. The work provides a simple way to implement an enhancement that effectively captures the influence of ion size and species into the PB theory in the context of PEs in aqueous salt solutions.

Assembly-based pathways of crystallization

Solution crystallization of materials ranging from simple salts to complex supramolecular assemblies has long been viewed through the lens of classical nucleation and growth theories in which monomeric building blocks...

A self-consistent Ornstein-Zernike jellium for highly charged colloids (microgels) in suspensions with added salt

This work discusses a jellium scheme, built within the framework of the multicomponent Ornstein-Zernike (OZ) equation, which is capable of describing the collective structure of suspensions of highly charged colloids with added salt, even in the presence of finite-size multivalent microions. This approach uses a suitable approximation to decouple the microion-microion correlations from the macroion-microion profiles, which in combination with the methodology from the dressed ion theory (DIT) gives a full account of the electrostatic effective potential among the colloids. The main advantages of the present contribution reside in its ability to manage the short-range potentials and non-linear correlations among the microions, as well as its realistic characterization of the ionic clouds surrounding each macroion. The structure factors predicted by this jellium scheme are contrasted with previously reported experimental results for microgel suspensions with monovalent salts (2019Phys. Rev. E100032602), thus validating its high accuracy in these situations. The present theoretical analysis is then extended to microgel suspensions with multivalent salts, which reveals the prominent influence of the counterion valence on the makeup of the effective potentials. Although the induced differences may be difficult to identify through the mesoscopic structure, our results suggest that the microgel collapsing transition may be used to enhance these distinct effects, thus giving a feasible experimental probe for these phenomena.

Chemodynamic features of nanoparticles: Application to understanding the dynamic life cycle of SARS-CoV-2 in aerosols and aqueous biointerfacial zones

We review concepts involved in describing the chemodynamic features of nanoparticles and apply the framework to gain physicochemical insights into interactions between SARS-CoV-2 virions and airborne particulate matter (PM). Our analysis is highly pertinent given that the World Health Organisation acknowledges that SARS-CoV-2 may be transmitted by respiratory droplets, and the US Center for Disease Control and Prevention recognises that airborne transmission of SARS-CoV-2 can occur. In our theoretical treatment, the virion is assimilated to a core-shell nanoparticle, and contributions of various interaction energies to the virion-PM association (electrostatic, hydrophobic, London-van der Waals, etc.) are generically included. We review the limited available literature on the physicochemical features of the SARS-CoV-2 virion and identify knowledge gaps. Despite the lack of quantitative data, our conceptual framework qualitatively predicts that virion-PM entities are largely able to maintain equilibrium on the timescale of their diffusion towards the host cell surface. Comparison of the relevant mass transport coefficients reveals that virion biointernalization demand by alveolar host cells may be greater than the diffusive supply. Under such conditions both the free and PM-sorbed virions may contribute to the transmitted dose. This result points to the potential for PM to serve as a shuttle for delivery of virions to host cell targets. Thus, our critical review reveals that the chemodynamics of virion-PM interactions may play a crucial role in the transmission of COVID-19, and provides a sound basis for explaining reported correlations between episodes of air pollution and outbreaks of COVID-19.

Osmotic stress and pore nucleation in charged biological nanoshells and capsids

A model system is proposed to investigate the chemical equilibrium and mechanical stability of biological spherical-like nanoshells in contact with an aqueous solution with added dissociated electrolyte of a given concentration. The ionic chemical equilibrium across the permeable shell is investigated in the framework of an accurate Density Functional Theory (DFT) that incorporates electrostatic and hardcore correlations beyond the traditional mean-field (e.g., Poisson-Boltzmann) limit. The accuracy of the theory is tested by a direct comparison with Monte Carlo (MC) simulations. A simple analytical expression is then deduced which clearly highlights the entropic, electrostatic, and self-energy contributions to the osmotic stress over the shell in terms of the calculated ionic profiles. By invoking a continuum mean-field elastic approach to account for the shell surface stress upon osmotic stretching, the mechanical equilibrium properties of the shell under a wide variety of ionic strengths and surface charges are investigated. The model is further coupled to a continuum mechanical approach similar in structure to a Classical Nucleation Theory (CNT) to address the question of mechanical stability of the shells against a pore nucleation. This allows us to construct a phase diagram which delimits the mechanical stability of capsids for different ionic strengths and shell surface charges.

Connecting energetics to dynamics in particle growth by oriented attachment using real-time observations

The interplay between crystal and solvent structure, interparticle forces and ensemble particle response dynamics governs the process of crystallization by oriented attachment (OA), yet a quantitative understanding is lacking. Using ZnO as a model system, we combine in situ TEM observations of single particle and ensemble assembly dynamics with simulations of interparticle forces and responses to relate experimentally derived interparticle potentials to the underlying interactions. We show that OA is driven by forces and torques due to a combination of electrostatic ion-solvent correlations and dipolar interactions that act at separations well beyond 5 nm. Importantly, coalignment is achieved before particles reach separations at which strong attractions drive the final jump to contact. The observed barrier to attachment is negligible, while dissipative factors in the quasi-2D confinement of the TEM fluid cell lead to abnormal diffusivities with timescales for rotation much less than for translation, thus enabling OA to dominate.