Charge regulation of colloidal particles in aqueous solutions

On the Validity of Constant pH Simulations

Constant pH (cpH) simulations are now a standard tool for investigating charge regulation in coarse-grained models of polyelectrolytes and colloidal systems. Originally developed for studying solutions with implicit ions, extending this method to systems with explicit ions or solvents presents several challenges. Ensuring proper charge neutrality within the simulation cell requires performing titration moves in sync with the insertion or deletion of ions, a crucial aspect often overlooked in the literature. Contrary to the prevailing views, cpH simulations are inherently grand-canonical, meaning that the controlled pH is that of the reservoir. The presence of the Donnan potential between the implicit reservoir and the simulation cell introduces significant differences between titration curves calculated for open and closed systems; the pH of an isolated (closed) system is different from the pH of the reservoir for the same protonation state of the polyelectrolyte. To underscore this point, in this paper, we will compare the titration curves calculated using the usual cpH algorithm with those from the exact canonical simulation algorithm. In the latter case, titration moves adhere to the correct detailed balance condition, and pH is calculated using the recently introduced surface Widom insertion algorithm. Our findings reveal a very significant difference between the titration isotherms obtained using the standard cpH algorithm and the canonical titration algorithm, emphasizing the importance of using the correct simulation approach when studying charge regulation of polyelectrolytes, proteins, and colloidal particles.

Canonical titration simulations

We present a Monte Carlo approach for performing titration simulations in the canonical ensemble. The standard constant pH (cpH) simulation methods are intrinsically grand canonical, allowing us to study the protonation state of molecules only as a function of pH in the reservoir. Due to the Donnan potential between a system and an (implicit) reservoir of a semi-grand canonical simulation, the pH of the reservoir can be significantly different from that of an isolated system, for an identical protonation state. The new titration method avoids this difficulty by using the canonical reactive Monte Carlo algorithm to calculate the protonation state of macromolecules as a function of the total number of protons present inside the simulation cell. The pH of an equilibrated system is then calculated using a new surface insertion Widom algorithm, which bypasses the difficulties associated with the bulk Widom particle insertion for intermediate and high pH values. To properly treat the long range Coulomb force, we use the Ewald summation method, showing the importance of the Bethe potential for calculating the pH of canonical systems.

Electrostatic interactions between charge regulated spherical macroions

We study the interaction between two charge regulating spherical macroions with dielectric interior and dissociable surface groups immersed in a monovalent electrolyte solution. The charge dissociation is modelled via the Frumkin-Fowler-Guggenheim isotherm, which allows for multiple adsorption equilibrium states. The interactions are derived from the solutions of the mean-field Poisson-Boltzmann type theory with charge regulation boundary conditions. For a range of conditions we find symmetry breaking transitions from symmetric to asymmetric charge distribution exhibiting annealed charge patchiness, which results in like-charge attraction even in a univalent electrolyte-thus fundamentally modifying the nature of electrostatic interactions in charge-stabilized colloidal suspensions.

Charge fluctuations in charge-regulated systems: dependence on statistical ensemble

We investigate charge regulation of nanoparticles in concentrated suspensions, focusing on the effect of different statistical ensembles. We find that the choice of ensemble does not affect the mean charge of nanoparticles, but significantly alters the magnitude of its fluctuation. Specifically, we compared the behaviors of colloidal charge fluctuations in the semi-grand canonical and canonical ensembles and identified significant differences between the two. The choice of ensemble-whether the system is isolated or is in contact with a reservoir of acid and salt-will, therefore, affect the Kirkwood-Shumaker fluctuation-induced force inside concentrated suspensions. Our results emphasize the importance of selecting an appropriate ensemble that accurately reflects the experimental conditions when studying fluctuation-induced forces between polyelectrolytes, proteins, and colloidal particles in concentrated suspensions.

Theory of Charge Regulation of Colloidal Particles in Electrolyte Solutions

We present a theory that enables us to (i) calculate the effective surface charge of colloidal particles and (ii) efficiently obtain titration curves for different salt concentrations. The theory accounts for the shift of pH of solution due to the presence of 1:1 electrolyte. It also accounts self-consistently for the electrostatic potential produced by the deprotonated surface groups. To examine the accuracy of the theory, we have performed extensive reactive Monte Carlo simulations, which show excellent agreement between theory and simulations without any adjustable parameters.

Brownian dynamics simulations of hard rods in external fields and with contact interactions

We propose a simulation method for Brownian dynamics of hard rods in one dimension for arbitrary continuous external force fields. It is an event-driven procedure based on the fragmentation and mergers of clusters formed by particles in contact. It allows one to treat particle interactions in addition to the hard-sphere exclusion as long as the corresponding interaction forces are continuous functions of the particle coordinates. We furthermore develop a treatment of sticky hard spheres as described by Baxter's contact interaction potential.

A method for quantifying how the activity of an enzyme is affected by the net charge of its nearest crowded neighbor

The electrostatic effects of protein crowding have not been systematically explored. Rather, protein crowding is generally studied with co‐solvents or crowders that are electrostatically neutral, with no methods to measure how the net charge (Z) of a crowder affects protein function. For example, can the activity of an enzyme be affected electrostatically by the net charge of its neighbor in crowded milieu? This paper reports a method for crowding proteins of different net charge to an enzyme via semi‐random chemical crosslinking. As a proof of concept, RNase A was crowded (at distances ≤ the Debye length) via crosslinking to different heme proteins with Z = +8.50 ± 0.04, Z = +6.39 ± 0.12, or Z = −10.30 ± 1.32. Crosslinking did not disrupt the structure of proteins, according to amide H/D exchange, and did not inhibit RNase A activity. For RNase A, we found that the electrostatic environment of each crowded neighbor had significant effects on rates of RNA hydrolysis. Crowding with cationic cytochrome c led to increases in activity, while crowding with anionic “supercharged” cytochrome c or myoglobin diminished activity. Surprisingly, electrostatic crowding effects were amplified at high ionic strength (I = 0.201 M) and attenuated at low ionic strength (I = 0.011 M). This salt dependence might be caused by a unique set of electric double layers at the dimer interspace (maximum distance of 8 Å, which cannot accommodate four layers). This new method of crowding via crosslinking can be used to search for electrostatic effects in protein crowding.

Applicability of the linearized Poisson-Boltzmann theory to contact angle problems and application to the carbon dioxide-brine-solid systems

In colloidal science and bioelectrostatics, the linear Poisson Boltzmann equation (LPBE) has been used extensively for the calculation of potential and surface charge density. Its fundamental assumption rests on the premises of low surface potential. In the geological sequestration of carbon dioxide in saline aquifers, very low pH conditions coupled with adsorption induced reduction of surface charge density result in low pH conditions that fit into the LPB theory. In this work, the Gouy–Chapman model of the electrical double layer has been employed in addition to the LPBE theory to develop a contact angle model that is a second-degree polynomial in pH. Our model contains the point of zero charge pH of solid surface. To render the model applicable to heterogeneous surfaces, we have further developed a model for the effective value of the point of zero charge pH. The point of zero charge pH model when integrated into our model enabled us to determine the point of zero charge pH of sandstone, quartz and mica using literature based experimental data. In this regard, a literature based thermodynamic model was used to calculate carbon dioxide solubility and pH of aqueous solution. Values of point of zero charge pH determined in this paper agree with reported ones. The novelty of our work stems from the fact that we have used the LPB theory in the context of interfacial science completely different from the classical approach, where the focus is on interparticle electrostatics involving colloidal stabilization.

Accelerated simulation method for charge regulation effects

The net charge of solvated entities, ranging from polyelectrolytes and biomolecules to charged nanoparticles and membranes, depends on the local dissociation equilibrium of individual ionizable groups. Incorporation of this phenomenon, charge regulation (CR), in theoretical and computational models requires dynamic, configuration-dependent recalculation of surface charges and is therefore typically approximated by assuming constant net charge on particles. Various computational methods exist that address this. We present an alternative, particularly efficient CR Monte Carlo method (CR-MC), which explicitly models the redistribution of individual charges and accurately samples the correct grand-canonical charge distribution. In addition, we provide an open-source implementation in the large-scale Atomic/Molecular Massively Parallel Simulator molecular dynamics (MD) simulation package, resulting in a hybrid MD/CR-MC simulation method. This implementation is designed to handle a wide range of implicit-solvent systems that model discreet ionizable groups or surface sites. The computational cost of the method scales linearly with the number of ionizable groups, thereby allowing accurate simulations of systems containing thousands of individual ionizable sites. By matter of illustration, we use the CR-MC method to quantify the effects of CR on the nature of the polyelectrolyte coil-globule transition and on the effective interaction between oppositely charged nanoparticles.

Reactive Monte Carlo simulations for charge regulation of colloidal particles

We use a reactive Monte Carlo simulation method and the primitive model of electrolyte to study acid-base equilibrium that controls charge regulation in colloidal systems. The simulations are performed in a semi-grand canonical ensemble in which colloidal suspension is in contact with a reservoir of salt and strong acid. The interior of colloidal particles is modeled as a low dielectric medium, different from the surrounding water. The effective colloidal charge is calculated for different numbers of surface acidic groups, pH, salt concentrations, and types of electrolyte. In the case of potassium chloride, the titration curves are compared with the experimental measurements obtained using potentiometric titration. A good agreement is found between simulations and experiments. In the case of lithium chloride, the specific ionic adsorption is taken into account through the partial dehydration of lithium ion.

Anion Specific Effects at Negatively Charged Interfaces: Influence of Cl-, Br-, I-, and SCN- on the Interactions of Na+ with the Carboxylic Acid Moiety

Unlike counterion interactions with charged interfaces, the influence of co-ions is only scarcely reported in the literature. In this work, the effect of SCN- and the halide co-ions in the interactions of Na+ with carboxylic acid Langmuir monolayers is investigated by using vibrational sum frequency spectroscopy. At 1 M concentrations in the subphase, the identity of the anion is shown to have a remarkable influence on the charging behavior and degree of deprotonation of the monolayer, with ions ordering in the sequence I- > SCN- > Cl- ≈ Br-. The same trend is observed at both pH 6 and pH 9 when the monolayer is intrinsically more charged. Spectroscopic evidence is found for both the presence of I- and SCN- in the interfacial region at levels close to their detection limits. The results contradict electrostatic theories on charged interfaces where co-ions are not expected to play any significant role. The higher propensity for the large polarizable anions to deprotonate the monolayer is explained in terms of their ability to modify the cations affinity toward the carboxylic acid groups present at the surface.

Measuring how two proteins affect each other's net charge in a crowded environment

Theory predicts that the net charge (Z) of a protein can be altered by the net charge of a neighboring protein as the two approach one another below the Debye length. This type of charge regulation suggests that a protein's charge and perhaps function might be affected by neighboring proteins without direct binding. Charge regulation during protein crowding has never been directly measured due to analytical challenges. Here, we show that lysine specific protein crosslinkers (NHS ester-Staudinger pairs) can be used to mimic crowding by linking two non-interacting proteins at a maximal distance of ~7.9 Å. The net charge of the regioisomeric dimers and preceding monomers can then be determined with lysine-acyl "protein charge ladders" and capillary electrophoresis. As a proof of concept, we covalently linked myoglobin (Z_monomer  = -0.43 ± 0.01) and α-lactalbumin (Z_monomer  = -4.63 ± 0.05). Amide hydrogen/deuterium exchange and circular dichroism spectroscopy demonstrated that crosslinking did not significantly alter the structure of either protein or result in direct binding (thus mimicking crowding). Ultimately, capillary electrophoretic analysis of the dimeric charge ladder detected a change in charge of ΔZ = -0.04 ± 0.09 upon crowding by this pair (Z_dimer  = -5.10 ± 0.07). These small values of ΔZ are not necessarily general to protein crowding (qualitatively or quantitatively) but will vary per protein size, charge, and solvent conditions.

Lattice-gas model of a charge regulated planar surface

In this work, we consider a lattice-gas model of charge regulation with electrostatic interactions within the Debye-Hückel level of approximation. In addition to long-range electrostatic interactions, the model incorporates the nearest-neighbor interactions for representing non-electrostatic forces between adsorbed ions. The Frumkin-Fowler-Guggenheim isotherm obtained from the mean-field analysis accurately reproduces the simulation data points.

Debye-Hückel Free Energy of an Electric Double Layer with Discrete Charges Located at a Dielectric Interface

Poisson-Boltzmann theory provides an established framework to calculate properties and free energies of an electric double layer, especially for simple geometries and interfaces that carry continuous charge densities. At sufficiently small length scales, however, the discreteness of the surface charges cannot be neglected. We consider a planar dielectric interface that separates a salt-containing aqueous phase from a medium of low dielectric constant and carries discrete surface charges of fixed density. Within the linear Debye-Hückel limit of Poisson-Boltzmann theory, we calculate the surface potential inside a Wigner-Seitz cell that is produced by all surface charges outside the cell using a Fourier-Bessel series and a Hankel transformation. From the surface potential, we obtain the Debye-Hückel free energy of the electric double layer, which we compare with the corresponding expression in the continuum limit. Differences arise for sufficiently small charge densities, where we show that the dominating interaction is dipolar, arising from the dipoles formed by the surface charges and associated counterions. This interaction propagates through the medium of a low dielectric constant and alters the continuum power of two dependence of the free energy on the surface charge density to a power of 2.5 law.