Charged patchy particle models in explicit salt: Ion distributions, electrostatic potentials, and effective interactions

Unusual Aspects of Charge Regulation in Flexible Weak Polyelectrolytes

This article reviews the state of the art of the studies on charge regulation (CR) effects in flexible weak polyelectrolytes (FWPE). The characteristic of FWPE is the strong coupling of ionization and conformational degrees of freedom. After introducing the necessary fundamental concepts, some unconventional aspects of the the physical chemistry of FWPE are discussed. These aspects are: (i) the extension of statistical mechanics techniques to include ionization equilibria and, in particular, the use of the recently proposed Site Binding-Rotational Isomeric State (SBRIS) model, which allows the calculation of ionization and conformational properties on the same foot; (ii) the recent progresses in the inclusion of proton equilibria in computer simulations; (iii) the possibility of mechanically induced CR in the stretching of FWPE; (iv) the non-trivial adsorption of FWPE on ionized surfaces with the same charge sign as the PE (the so-called "wrong side" of the isoelectric point); (v) the influence of macromolecular crowding on CR.

Interaction of Heparin with Proteins: Hydration Effects

We present a thermodynamic investigation of the interaction of heparin with lysozyme in the presence of potassium glutamate (KGlu). The binding constant K_b is measured by isothermal titration calorimetry (ITC) in a temperature range from 288 to 310 K for concentrations of KGlu between 25 and 175 mM. The free energy of binding ΔG_b derived from K_b is strongly decreasing with increasing concentration of KGlu, whereas the dependence of ΔG_b on temperature T is found to be small. The decrease of ΔG_b can be explained in terms of counterion release: Binding of lysozyme to the strong polyelectrolyte heparin liberates approximately three of the condensed counterions of heparin, thus increasing the entropy of the system. The dependence of ΔG_b on T, on the other hand, is traced back to a change of hydration of the protein and the polyelectrolyte upon complex formation. This dependence is quantitatively described by the parameter Δw that depends on T and vanishes at a characteristic temperature T₀. A comparison of the complex formation in the presence of KGlu with the one in the presence of NaCl demonstrates that the parameters related to hydration are changed considerably. The characteristic temperature T₀ in the presence of KGlu solutions is considerably smaller than that in the presence of NaCl solutions. The change of specific heat Δc_p is found to become more negative with increasing salt concentration: This finding agrees with the model-free analysis by the generalized van't Hoff equation. The entire analysis reveals a small but important change of the free energy of binding by hydration. It shows that these ion-specific Hofmeister effects can be modeled quantitatively in terms of a characteristic temperature T₀ and a parameter describing the dependence of Δc_p on salt concentration.

Self-Interactions of Two Monoclonal Antibodies: Small-Angle X-ray Scattering, Light Scattering, and Coarse-Grained Modeling

Using light scattering (LS), small-angle X-ray scattering (SAXS), and coarse-grained Monte Carlo (MC) simulations, we studied the self-interactions of two monoclonal antibodies (mAbs), PPI03 and PPI13. With LS measurements, we obtained the osmotic second virial coefficient, B₂₂, and the molecular weight, M_w, of the two mAbs, while with SAXS measurements, we studied the mAbs' self-interaction behavior in the high protein concentration regime up to 125 g/L. Through SAXS-derived coarse-grained representations of the mAbs, we performed MC simulations with either a one-protein or a two-protein model to predict B₂₂. By comparing simulation and experimental results, we validated our models and obtained insights into the mAbs' self-interaction properties, highlighting the role of both ion binding and charged patches on the mAb surfaces. Our models provide useful information about mAbs' self-interaction properties and can assist the screening of conditions driving to colloidal stability.

Toward a Suite of Coarse-Grained Models for Molecular Simulation of Monoclonal Antibodies and Therapeutic Proteins

A series of coarse-grained models for molecular simulation of proteins are considered, with emphasis on the application of predicting protein-protein self-interactions for monoclonal antibodies (MAbs). As an illustrative example and for quantitative comparison, the models are used to predict osmotic virial coefficients over a broad range of attractive and repulsive self-interactions and solution conditions for a series of MAbs where the second osmotic virial coefficient has been experimentally determined in prior work. The models are compared based on how well they can predict experimental behavior, their computational burdens, and scalability. An intermediate-resolution model is also introduced that can capture specific electrostatic interactions with improved efficiency and similar or improved accuracy when compared to the previously published models. Guidance is included for the selection of coarse-grained models more generally for capturing a balance of electrostatic, steric, and short-ranged nonelectrostatic interactions for proteins from low to high concentrations.

A modeling approach for quantitative assessment of interfacial interaction between two rough particles in colloidal systems

HYPOTHESIS AND BACKGROUND

The simulation of rough particle surface is important to understand and control the interface behavior of particles in colloidal systems. Literature analysis suggested a lack of information for an accurate model simulating the interfacial interaction between two rough particles. It is hypothesized that the total interfacial energy developed between two rough particles would depend on the surface morphologies of particles, and it could be predicted if a mathematical model to represent the interaction of two rough particles were created accurately.

EXPERIMENTS

In this study, mathematical models were developed to determine the interfacial energy created between two particles according to the XDLVO theory by considering the rippled particle theory and surface element integral (SEI) method. Three different scenarios of particle interactions were assumed in the simulation. The present study provides deep insights into particle interactions via considering aspect ratio, size, and surface roughness of two particles in colloidal systems.

FINDINGS

The assessment of the interfacial interaction revealed that an increase in the aspect ratio, surface roughness, and relative surface roughness of particles would weaken the total interaction energy generated between particles and promote particle aggregation. Increased interaction energy was predicted for the interaction of particles by increasing the particle size. The asperity ratio was more effective than the asperity number in controlling the interfacial energy between two particles. The results of this study could be used for foreseeing the interaction of rough particles, which has a significant application in particle coagulation or dispersion in colloidal systems.

Pearl-necklace assembly of human serum albumin with the poly(acrylic acid) polyelectrolyte investigated using small angle X-ray scattering (SAXS)

In this comprehensive study, the interaction of human serum albumin (HSA) with poly(acrylic acid) (PAA) was explored using small angle X-ray scattering (SAXS) combined with chromatography. The results revealed the formation of a complex between HSA macromolecules and PAA chains but solely under some specific conditions of the ionic strength and pH of the medium. In fact, this binding was found to take place only at pH close to 5 and at low ionic strength (0.15 M). Otherwise, for a higher pH and a salt concentration of 0.75 M the HSA-PAA complex tends to dissociate completely showing the reversibility of the complexation. The assessment of the influence of the HSA/PAA molar ratio on the radius of gyration of the complex suggests that 4 HSA molecules could bind to each 100 kDa PAA chain. In addition, the Porod volume evaluation for the same range of the HSA/PAA ratio confirms this assumption. Finally, an all-atom SAXS modelling study using the BUNCH program was conducted to find a compatible model that fits the HSA-PAA complex scattering data. This model allows us to portray the HSA/PAA complex as a pearl-necklace assembly with 4 HSA molecules on the 100 kDa PAA chain.

Calculation of salt-dependent free energy of binding of β-lactoglobulin homodimer formation and mechanism of dimer formation using molecular dynamics simulation and three-dimensional reference interaction site model (3D-RISM): diffuse salt ions and non-polar interactions between the monomers favor the dimer formation

There are several important phenomena in chemistry, biology, and physics where molecules (or parts of a molecule) having charges of the same sign come closer together and become stable. DNA condensation, RNA folding, colloid-colloid interactions are some of the examples of this kind. In the current work, we have investigated how β-lactoglobulin, a protein found in milk, in spite of carrying +13 charge, favors the homodimer formation in the presence of salt. We have focussed on calculating the protein-protein binding free energy in the presence of salt and identifying the thermodynamic and microscopic mechanism of the process. Estimation of binding free energy of this salt-dependent process is done by combining molecular dynamics simulation with statistical mechanical theory of three-dimensional reference interaction site model (3D-RISM). Binding free energy is evaluated from the chemical potential of the solutes as opposed to potential of mean force calculation, which gives only a constrained free energy. Our calculated values semi-quantitatively match with the experimental results. By examining the different components of binding free energy, we have found that the role of salt ions (especially of Cl^-) is to shift the equilibrium towards the dimer. Non-polar (Lennard-Jones) interactions between the monomers is also favorable to the binding free energy. However, water slightly disfavors the dimer formation. For the microscopic mechanism, heterogeneous of both Na⁺ and Cl^- near the charged residues at the binding interface and change of this charge distribution on dimer formation contribute to the stability. A fine-tuning of enthalpic and entropic effects of salt ions is found to operate at different salt concentrations. Both thermodynamic and microscopic mechanism of dimer formation gives detailed insight into the complex electrostatics of charged protein-protein binding.

Interaction of Proteins with Polyelectrolytes: Comparison of Theory to Experiment

We discuss recent investigations of the interaction of polyelectrolytes with proteins. In particular, we review our recent studies on the interaction of simple proteins such as human serum albumin (HSA) and lysozyme with linear polyelectrolytes, charged dendrimers, charged networks, and polyelectrolyte brushes. In all cases discussed here, we combined experimental work with molecular dynamics (MD) simulations and mean-field theories. In particular, isothermal titration calorimetry (ITC) has been employed to obtain the respective binding constants K_b and the Gibbs free energy of binding. MD simulations with explicit counterions but implicit water demonstrate that counterion release is the main driving force for the binding of proteins to strongly charged polyelectrolytes: patches of positive charges located on the surface of the protein become multivalent counterions of the polyelectrolyte, thereby releasing a number of counterions condensed on the polyelectrolyte. The binding Gibbs free energy due to counterion release is predicted to scale with the logarithm of the salt concentration in the system, which is verified by both simulations and experiment. In several cases, namely, for the interaction of proteins with linear polyelectrolytes and highly charged hydrophilic dendrimers, the binding constant could be calculated from simulations to very good approximation. This finding demonstrated that in these cases explicit hydration effects do not contribute to the Gibbs free energy of binding. The Gibbs free energy can also be used to predict the kinetics of protein uptake by microgels for a given system by applying dynamic density functional theory. The entire discussion demonstrates that the direct comparison of theory with experiments can lead to a full understanding of the interaction of proteins with charged polymers. Possible implications for applications, such as drug design, are discussed.

Mono-patchy zwitterionic microcolloids as building blocks for pH-controlled self-assembly

A directional molecular interaction between microcolloids can be achieved through pre-defined sites on their surface, "patches", which might make them follow each other in a controlled way and assemble into target structures of more complexity. In this article, we report the successful generation and characterization of mono-patchy melamine-formaldehyde microparticles with oppositely charged patches made of poly(methyl vinyl ether-alt-maleic acid) or polyethyleneimine via microcontact printing. The study of their self-aggregation behavior in solution shows that by change of pH, particle dimers are formed via attractive electrostatic force between the patchy and non-patchy surface of the particles, which reaches its optimum at a specific pH.

Tuning the order of colloidal monolayers: assembly of heterogeneously charged colloids close to a patterned substrate

We study the behavior of negatively charged colloids with two positively charged polar caps close to a planar patterned surface. The competition between the different anisotropic components of the particle-particle interaction is able by itself to give rise to a rich assembly scenario: colloids with charged surface patterns already form different crystalline domains when adsorbed to a homogeneously charged substrate. Here we consider substrates composed of alternating (negative/neutral, positive/neutral and positive/negative) parallel stripes and, by means of Monte Carlo simulations, we investigate the ordering of the colloids on changing the number of the stripes. We show that the additional competition between the two different lengths scales characterizing the system (i.e., the particle interaction range and the size of the stripes) gives rise to a plethora of distinct particle arrangements, where some well-defined trends can be observed. By accurately tuning the substrate charged motif it is possible to, e.g., promote specific particle arrangements, disfavor crystalline domains or induce the formation of extended, open clusters.

From discrete to continuous description of spherical surface charge distributions

The importance of electrostatic interactions in soft matter and biological systems can often be traced to non-uniform charge effects, which are commonly described using a multipole expansion of the corresponding charge distribution. The standard approach when extracting the charge distribution of a given system is to treat the constituent charges as points. This can, however, lead to an overestimation of multipole moments of high order, such as dipole, quadrupole, and higher moments. Focusing on distributions of charges located on a spherical surface - characteristic of numerous biological macromolecules, such as globular proteins and viral capsids, as well as of inverse patchy colloids - we develop a novel way of representing spherical surface charge distributions based on the von Mises-Fisher distribution. This approach takes into account the finite spatial extension of individual charges, and leads to a simple yet powerful way of describing surface charge distributions and their multipole expansions. In this manner, we analyze charge distributions and the derived multipole moments of a number of different spherical configurations of identical charges with various degrees of symmetry. We show how the number of charges, their size, and the geometry of their configuration influence the behavior and relative importance of multipole magnitudes of different order. Importantly, we clearly demonstrate how neglecting the effect of charge size leads to an overestimation of high-order multipoles. The results of our work can be applied to construct analytical models of electrostatic interactions and multipole expansion of charged particles in diverse soft matter and biological systems.

pH Dependence of Charge Multipole Moments in Proteins

Electrostatic interactions play a fundamental role in the structure and function of proteins. Due to ionizable amino acid residues present on the solvent-exposed surfaces of proteins, the protein charge is not constant but varies with the changes in the environment-most notably, the pH of the surrounding solution. We study the effects of pH on the charge of four globular proteins by expanding their surface charge distributions in terms of multipoles. The detailed representation of the charges on the proteins is in this way replaced by the magnitudes and orientations of the multipole moments of varying order. Focusing on the three lowest-order multipoles-the total charge, dipole, and quadrupole moment-we show that the value of pH influences not only their magnitudes, but more notably and importantly also the spatial orientation of their principal axes. Our findings imply important consequences for the study of protein-protein interactions and the assembly of both proteinaceous shells and patchy colloids with dissociable charge groups.

Investigating the thermal dissociation of viral capsid by lattice model

The dissociation of icosahedral viral capsids was investigated by a homogeneous and a heterogeneous lattice model. In thermal dissociation experiments with cowpea chlorotic mottle virus and probed by small-angle neutron scattering, we observed a slight shrinkage of viral capsids, which can be related to the strengthening of the hydrophobic interaction between subunits at increasing temperature. By considering the temperature dependence of hydrophobic interaction in the homogeneous lattice model, we were able to give a better estimate of the effective charge. In the heterogeneous lattice model, two sets of lattice sites represented different capsid subunits with asymmetric interaction strengths. In that case, the dissociation of capsids was found to shift from a sharp one-step transition to a gradual two-step transition by weakening the hydrophobic interaction between AB and CC subunits. We anticipate that such lattice models will shed further light on the statistical mechanics underlying virus assembly and disassembly.

Computationally Mapping pKa Shifts Due to the Presence of a Polyelectrolyte Chain around Whey Proteins

Experimental studies have shown the formation of soluble complexes in the pure repulsive Coulombic regime even when the net charges of the protein and the polyelectrolyte have the same sign ( De Kruif et al. Curr. Opin. Colloid Interface Sci. 2004 , 9 , 340 ; De Vries et al. J. Chem. Phys. 2003 , 118 , 4649 ; Grymonpre et al. Biomacromolecules 2001 , 2 , 422 ; Hattori et al. Langmuir 2000 , 16 , 9738 ). This attractive phenomenon has often been described as "complexation on the wrong side of pI". While one theory assumes the existence of "charged patches" on the protein surface from ion-dipole interactions, thus allowing a polyelectrolyte to bind to an oppositely heterogeneous charged protein region, another theoretical view considers the induced-charge interactions to be the dominant factor in these complexations. This charge regulation mechanism can be described by proton fluctuations resulting from mutual rearrangements of the distributions of the charged groups, due to perturbations of the acid-base equilibrium. Using constant-pH Monte Carlo simulations and several quantitative and visual analysis tools, we investigate the significance of each of these interactions for two whey proteins, α-lactalbumin (α-LA) and lysozyme (LYZ). Through physical chemistry parameters, free energies of interactions, and the mapping of amino acid pK_a shifts and polyelectrolyte trajectories, we show the charge regulation mechanism to be the most important contributor in protein-polyelectrolyte complexation regardless of pH, dipole moment, and protein capacitance in a low salt regime.

Electrostatics of patchy surfaces

In the study of colloidal, biological and electrochemical systems, it is customary to treat surfaces, macromolecules and electrodes as homogeneously charged. This simplified approach is proven successful in most cases, but fails to describe a wide range of heterogeneously charged surfaces commonly used in experiments. For example, recent experiments have revealed a long-range attraction between overall neutral surfaces, locally charged in a mosaic-like structure of positively and negatively charged domains ("patches"). Here, we review experimental and theoretical studies addressing the stability of heterogeneously charged surfaces, their effect on ionic profiles in solution, and the interaction between two such surfaces. We focus on electrostatics, and highlight the important new physical parameters appearing in the heterogeneous case, such as the largest patch size and inter-surface charge correlations.

Protein Immobilization onto Cationic Spherical Polyelectrolyte Brushes Studied by Small Angle X-ray Scattering

The immobilization of bovine serum albumins (BSA) onto cationic spherical polyelectrolyte brushes (SPB) consisting of a solid polystyrene (PS) core and a densely grafted poly(2-aminoethyl methacrylate hydrochloride) (PAEMH) shell was studied by small-angle X-ray scattering (SAXS). The observed dynamics of adsorption of BSA onto SPB by time-resolved SAXS can be divided into two stages. In the first stage (tens of milliseconds), the added proteins as in-between bridge instantaneously caused the aggregation of SPB. Then BSA penetrated into the brush layer driven by electrostatic attractions, and reached equilibrium in the second stage (tens of seconds). The amount of BSA immobilized onto brush layer reached the maximum when pH was increased to about 6.1 and BSA concentration to 10 g/L. The cationic SPB were confirmed to provide stronger adsorption capacity for BSA compared to anionic ones.

Interaction of Charged Patchy Protein Models with Like-Charged Polyelectrolyte Brushes

We study the adsorption of charged patchy particle models (CPPMs) on a thin film of a like-charged and dense polyelectrolyte (PE) brush (of 50 monomers per chain) by means of implicit-solvent, explicit-salt Langevin dynamics computer simulations. Our previously introduced set of CPPMs embraces well-defined one- and two-patched spherical globules, each of the same net charge and (nanometer) size, with mono- and multipole moments comparable to those of small globular proteins. We focus on electrostatic effects on the adsorption far away from the isoelectric point of typical proteins, i.e., where charge regulation plays no role. Despite the same net charge of the brush and globule, we observe large binding affinities up to tens of the thermal energy, k_BT, which are enhanced by decreasing salt concentration and increasing charge of the patch(es). Our analysis of the distance-resolved potentials of mean force together with a phenomenological description of all leading interaction contributions shows that the attraction is strongest at the brush surface, driven by multipolar, Born (self-energy), and counterion-release contributions, dominating locally over the monopolar and steric repulsions.

Limiting the valence: advancements and new perspectives on patchy colloids, soft functionalized nanoparticles and biomolecules

Artistic representation of limited valance units consisting of a soft core (in blue) and a small number of flexible bonding patches (in orange).

Effect of the surface charge distribution on the fluid phase behavior of charged colloids and proteins

A generic but simple model is presented to evaluate the effect of the heterogeneous surface charge distribution of proteins and zwitterionic nanoparticles on their thermodynamic phase behavior. By considering surface charges as continuous "patches," the rich set of surface patterns that is embedded in proteins and charged patchy particles can readily be described. This model is used to study the fluid phase separation of charged particles where the screening length is of the same order of magnitude as the particle size. In particular, two types of charged particles are studied: dipolar fluids and protein-like fluids. The former represents the simplest case of zwitterionic particles, whose charge distribution can be described by their dipole moment. The latter system corresponds to molecules/particles with complex surface charge arrangements such as those found in biomolecules. The results for both systems suggest a relation between the critical region, the strength of the interparticle interactions, and the arrangement of charged patches, where the critical temperature is strongly correlated to the magnitude of the dipole moment. Additionally, competition between attractive and repulsive charge-charge interactions seems to be related to the formation of fluctuating clusters in the dilute phase of dipolar fluids, as well as to the broadening of the binodal curve in protein-like fluids. Finally, a variety of self-assembled architectures are detected for dipolar fluids upon small changes to the charge distribution, providing the groundwork for studying the self-assembly of charged patchy particles.

Coupling and Decoupling of Rotational and Translational Diffusion of Proteins under Crowding Conditions

Molecular motion of biopolymers in vivo is known to be strongly influenced by the high concentration of organic matter inside cells, usually referred to as crowding conditions. To elucidate the effect of intermolecular interactions on Brownian motion of proteins, we performed (1)H pulsed-field gradient NMR and fluorescence correlation spectroscopy (FCS) experiments combined with small-angle X-ray scattering (SAXS) and viscosity measurements for three proteins, αB-crystalline (αBc), bovine serum albumin, and hen egg-white lysozyme (HEWL) in aqueous solution. Our results demonstrate that long-time translational diffusion quantitatively follows the expected increase of macro-viscosity upon increasing the protein concentration in all cases, while rotational diffusion as assessed by polarized FCS and previous multi-frequency (1)H NMR relaxometry experiments reveals protein-specific behavior spanning the full range between the limiting cases of full decoupling from (αBc) and full coupling to (HEWL) the macro-viscosity. SAXS was used to study the interactions between the proteins in solution, whereby it is shown that the three cases cover the range between a weakly interacting hard-sphere system (αBc) and screened Coulomb repulsion combined with short-range attraction (HEWL). Our results, as well as insights from the recent literature, suggest that the unusual rotational-translational coupling may be due to anisotropic interactions originating from hydrodynamic shape effects combined with high charge and possibly a patchy charge distribution.

Electrostatic attraction between overall neutral surfaces

Two overall neutral surfaces with positively and negatively charged domains ("patches") have been shown in recent experiments to exhibit long-range attraction when immersed in an ionic solution. Motivated by the experiments, we calculate analytically the osmotic pressure between such surfaces within the Poisson-Boltzmann framework, using a variational principle for the surface-averaged free energy. The electrostatic potential, calculated beyond the linear Debye-Hückel theory, yields an overall attraction at large intersurface separations, over a wide range of the system's controlled length scales. In particular, the attraction is stronger and occurs at smaller separations for surface patches of larger size and charge density. In this large patch limit, we find that the attraction-repulsion crossover separation is inversely proportional to the square of the patch-charge density and to the Debye screening length.

Like-charged protein-polyelectrolyte complexation driven by charge patches

We study the pair complexation of a single, highly charged polyelectrolyte (PE) chain (of 25 or 50 monomers) with like-charged patchy protein models (CPPMs) by means of implicit-solvent, explicit-salt Langevin dynamics computer simulations. Our previously introduced set of CPPMs embraces well-defined zero-, one-, and two-patched spherical globules each of the same net charge and (nanometer) size with mono- and multipole moments comparable to those of globular proteins with similar size. We observe large binding affinities between the CPPM and the like-charged PE in the tens of the thermal energy, kBT, that are favored by decreasing salt concentration and increasing charge of the patch(es). Our systematic analysis shows a clear correlation between the distance-resolved potentials of mean force, the number of ions released from the PE, and CPPM orientation effects. In particular, we find a novel two-site binding behavior for PEs in the case of two-patched CPPMs, where intermediate metastable complex structures are formed. In order to describe the salt-dependence of the binding affinity for mainly dipolar (one-patched) CPPMs, we introduce a combined counterion-release/Debye-Hückel model that quantitatively captures the essential physics of electrostatic complexation in our systems.

Tuning protein-protein interactions using cosolvents: specific effects of ionic and non-ionic additives on protein phase behavior

Cosolvents are routinely used to modulate the (thermal) stability of proteins and, hence, their interactions with proteins have been studied intensely. However, less is known about their specific effects on protein-protein interactions, which we characterize in terms of the protein phase behavior. We analyze the phase behavior of lysozyme solutions in the presence of sodium chloride (NaCl), guanidine hydrochloride (GuHCl), glycerol, and dimethyl sulfoxide (DMSO). We experimentally determined the crystallization boundary (XB) and, in combination with data on the cloud-point temperatures (CPTs), the crystallization gap. In agreement with other studies, our data indicate that the additives might affect the protein phase behavior through electrostatic screening and additive-specific contributions. At high salt concentrations, where electrostatic interactions are screened, both the CPT and the XB are found to be linear functions of the additive concentration. Their slopes quantify the additive-specific changes of the phase behavior and thus of the protein-protein interactions. While the specific effect of NaCl is to induce attractions between proteins, DMSO, glycerol and GuHCl (with increasing strength) weaken attractions and/or induce repulsions. Except for DMSO, changes of the CPT are stronger than those of the XB. Furthermore, the crystallization gap widens in the case of GuHCl and glycerol and narrows in the case of NaCl. We relate these changes to colloidal interaction models, namely square-well and patchy interactions.

Anisotropic protein-protein interactions due to ion binding

Self-association of proteins is strongly affected by long-range electrostatic interactions caused by equilibrium adsorption of small ions such as protons and multivalent metals. By affecting the molecular net charge, solution pH is thus a widely used parameter to tune stability and phase behavior of proteins. We here review recent studies where the charge distribution is perturbed not only by protons, but also by other binding ions, leading to a rich and inherently anisotropic charge distribution. Focus is on coarse grained simulation techniques, coupled to experiments of protein-protein interaction at varying salt and pH conditions. Finally, and with future bio-colloidal models in mind, we discuss the validity of coarse graining charge anisotropy using electric multipoles.