# Interaction of Protein-like Nanocolloids with pH-Sensitive Polyelectrolyte Brushes

Tatiana O Popova, Ekaterina B Zhulina, Oleg V. Borisov

► To cite this version:

Tatiana O Popova, Ekaterina B Zhulina, Oleg V. Borisov. Interaction of Protein-like Nanocolloids with pH-Sensitive Polyelectrolyte Brushes. International Journal of Molecular Sciences, 2025, 26, pp.7867. 10.3390/ijms26167867. hal-05370894

HAL Id: hal-05370894

https://hal.science/hal-05370894v1

Submitted on 18 Nov 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

# Interaction of Protein-like Nanocolloids with pH-Sensitive Polyelectrolyte Brushes

Tatiana O. Popova <fnref n="1" />,<fnref n="2" /> ID, Ekaterina B. Zhulina <fnref n="2" /> ID and Oleg V. Borisov <fnref n="1" />,<fnref n="2" />,<fnref n="3" />,*

## Abstract

The self-consistent field Poisson-Boltzmann framework is applied for analysis of equilibrium partitioning of ampholytic protein-like nanocolloids between buffer solution and weak (pH-sensitive) versus strong polyelectrolyte (polyanionic) brushes with the same net charge per unit area. The position-dependent nanocolloid net charge and the insertion freeenergy profiles are derived as a function of pH and ionic strength in the solution. It is demonstrated that, similar to strong polyelectrolyte brushes, pH-sensitive brushes are capable of the uptake of nanocolloids in the vicinity of the isoelectric point, that is, when the net charge of the colloid in the buffer has either the opposite or the same sign as the ionized monomer units of the brush. At $pI \ge pK_{brush}$ and $pH \ge pI$, the particle absorption patterns by similarly (negatively) charged brushes are qualitatively similar in the cases of strong and weak polyelectrolyte brushes, but the freeenergy barrier at the brush periphery is wider for weak than for strong polyelectrolyte brushes, which may cause stronger kinetic hindrance for the nanocolloid uptake by the brush. A decrease in pH below the IEP leads to a monotonic increase in the depth of the insertion freeenergy minimum inside a strong polyelectrolyte brush, whereas for weak polyelectrolyte brushes, a more peculiar trend is predicted: due to competition between the increasing positive charge of the nanocolloid and the decreasing magnitude of the negative charge of the brush, the absorption is weakened at low pH.

**Keywords:** nanocolloids; globular proteins; polyampholytes; polyelectrolyte brushes; charge regulation; mean-field theory

## 1. Introduction

Electrostatically driven interaction and complexation of globular proteins with biological and synthetic polyelectrolytes play a prominent role in nature as well as in many technological and biomedical applications [1]. These assemblies, including immobilization of proteins, enzymes, and antibodies on solid supports, are of particular interest in the fields of biotechnology and medicine. However, when biomolecules adhere to solid surfaces, their biological function might be compromised.

As a promising alternative, immobilization of globular proteins in polyelectrolyte brushes (monolayers of charged macromolecules end-attached to planar substrates or surfaces of colloidal particles, and immersed in aqueous solutions [2–7]) has been actively explored over the last few decades [8–11] (see also ref. [12] for a comprehensive review).

<figure><img src="image_3.png" /></figure>

### check for updates

Academic Editor: Sebastian Małkowski

Received: 2 July 2025

Revised: 1 August 2025

Accepted: 7 August 2025

Published: 14 August 2025

**Citation:** Popova, T.O.; Zhulina, E.B.; Borisov, O.V. Interaction of Protein-like Nanocolloids with pH-Sensitive Polyelectrolyte Brushes. *Int. J. Mol. Sci.* **2025**, *26*, 7867. https://doi.org/10.3390/ijms26167867

**Copyright:** © 2025 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

<footnote n="1">¹ Center for Chemical Engineering, ITMO University, 197101 St. Petersburg, Russia; salamatovat170301@gmail.com</footnote>

<footnote n="2">² NRC Kurchatov Institute—PNPI—IMC, 199004 St. Petersburg, Russia</footnote>

<footnote n="3">³ Institut des Sciences Analytiques et de Physico-Chimie pour l'Environnement et les Matériaux, CNRS, Université de Pau et des Pays de l'Adour, UMR 5254, 64053 Pau, France</footnote>

<footnote>* Correspondence: oleg.borisov@univ-pau.fr</footnote>
It was found that when biomolecules are bound to polyelectrolyte brushes, their enzymatic activity is often preserved [13–16], which opens a fascinating perspective of using polyelectrolyte brushes in supported enzymatic catalysis.

The most important experimental findings that motivated and guided theoretical research on the interaction of globular proteins with polyelectrolyte (PE) brushes are (i) proteins can be absorbed by both polyanionic and polycationic brushes in the vicinity of the isoelectric point (IEP), i.e., when the charge of the protein globule is not only opposite but also of the same sign as that of the brush-forming chains. Only at $pH_b$ sufficiently far away from the IEP proteins are repelled from the similarly charged brush; and (ii) an increase in the ionic strength of the solution suppresses protein uptake by a polyelectrolyte brush on both sides of the protein IEP. These findings unambiguously indicated that protein uptake by PE brushes is electrostatically driven on both sides of the IEP.

While the interaction of proteins with an oppositely charged brush is plausibly governed by Coulomb attraction, the mechanism of protein uptake by a similarly charged brush is more delicate and remains a matter of discussion in the literature. Two possible explanations of this effect have been proposed. According to the first explanation [17–19], it is the inhomogeneity of the charge distribution on the protein surface that drives the absorption of proteins by the PE brush. Although the net charge of the protein on the “wrong side” of the IEP is of the same sign as that of the brush-forming chains, patches of charge of opposite sign on the globule surface drive the adsorption of PE chains on these patches, overcompensating for the penalty of polyelectrolyte repulsion from similarly charged patches. Evidently, this mechanism is operative provided that the distribution of anionic and cationic residues on the globule is essentially inhomogeneous.

An alternative mechanism [20,21] is based on the fact that the local pH inside the PE brush differs from that in the buffer solution due to the strong electrostatic field inside the brush. As a result, the net charge of the protein globule (controlled by the balance of ionization of pH-sensitive cationic and anionic residues) may change its sign inside the brush, so that the globule becomes charged oppositely to the brush and is driven inside by the Coulomb force. This second mechanism is more general and does not require any pre-assumption concerning the distribution of cationic and anionic residues on the globule surface. However, as it has been proven by a more accurate theoretical analysis [22–25], an inversion of the globule charge in the brush is an indispensable but not sufficient condition for the globule uptake by a similarly charged PE brush. As was demonstrated, the charge inversion leads to the appearance of a free energy barrier that may hinder protein uptake by the brush even when it is thermodynamically favorable.

Both experimental and theoretical studies reveal that the buffer $pH_b$ is the key control parameter in the interaction of globular proteins and ampholytic protein-like nanocolloids with polyelectrolyte brushes. When the PE chains are negatively charged (polyacid brush), the electrostatic driving force for globule absorption is suppressed under strongly basic conditions ($pH_b \gg pI$). It becomes operative in the vicinity of the IEP, $pH_b \ge pI$, and is further promoted upon a decrease in $pH_b < pI$. This trend is attributed to the progressively increasing net charge of the nanocolloid — from negative at $pH_b \ge pI$ to positive at $pH_b < pI$ — in the buffer solution and, with a positive offset, inside the negatively charged brush. The reverse trend is observed for positively charged (polybase) brushes.

In the theory developed earlier [22–25], it was, however, assumed that the charge on the brush-forming chains does not depend on $pH_b$, as it happens for strong (quenched) polyelectrolytes, like polystyrenesulfonate (PSS, anionic) or poly(diallyldimethylammonium chloride) (PDADMAC, cationic). As a result, the electrostatic field created by a quenched PE brush and acting on the ampholytic particle is also independent of $pH_b$. In many experiments, however, PE brushes were formed by weak (pH-sensitive) polyelectrolytes, like
poly(acrylic) acid (PAA, anionic) [26–28] or poly (2-aminoethyl methacrylate hydrochloride) (PAEMH, cationic) [28,29].

Evidently, if $pK_a$ of a pH-sensitive polyanionic/polycationic brush is significantly larger/smaller than $pI$ of the nanocolloid, the PE brush is virtually fully ionized in the vicinity of $pI$. In this case, the interactions of globular proteins with such a brush are similar to those with a strong polyelectrolyte (permanently ionized) brush. The situation may become more delicate if $pK_a \approx pI$, that is, variation in pH around the protein IEP is accompanied by a significant change in the ionization of the brush-forming PE chains. One should also bear in mind that due to cooperative intermolecular Coulomb interactions between brush-forming pH-sensitive PE chains, the effective $pK_a^{brush}$ is lower (for polyacid) or higher (for polybase) than the “bare” $pK_a$ for an individual monomer unit in the buffer and depends on the brush architecture (chain length, grafting density) and the ionic strength of the solution [30,31].

The aim of the present theoretical study is to unravel the effect of variations in buffer $pH_b$ on the interactions and uptake of protein-like ampholytic nanocolloids by brushes formed by pH-sensitive (annealing) polyanionic chains with a bare ionization constant $pK_a \approx pI$. We aim to investigate how simultaneous variations in the ionization of nanocolloid cationic and anionic residues and of monomer units in the brush-forming chains manifest themselves in the pH-controlled shape of the insertion free energy profiles and thus in the switch from absorption to depletion scenarios. To reach these goals, we implement an analytical self-consistent field Poisson-Boltzmann framework [31–33] combined with the numerical Scheutjens-Fleer calculations [34]. The latter approach must be applied under conditions of suppressed ionization of PE chains, i.e., when the free energy contributions due to the electrostatic and the excluded volume intermolecular interactions are comparable.

## 2. Results and Discussion

In Figure 1, the electrostatic potential profiles $\psi(z)$ are presented for an annealing (weak) anionic brush at different offsets $\Delta pH_{brush} = pH_b - pK_{brush}$ from the brush $pK(pK_{brush} = 5.5)$ and for a quenched (strong) PE brush with fractions $\alpha$ of permanently charged monomer units equal to the average degrees of ionization $\langle\alpha(z)\rangle$ of the weak PE brush at each respective value of $pH_b$. Four offsets were considered: $\Delta pH_{brush} = -0.4, 0, 0.4$, and 0.8 at fixed $pK_{brush} = 5.5$. Solid lines correspond to strong-stretching analytical theory predictions, while dashed lines represent numerical results obtained from the Scheutjens-Fleer self-consistent field (SF-SCF) method for the annealing brush. The dash-dot lines represent the electrostatic potential profiles for quenched PE brushes, as obtained from the strong-stretching analytical theory.

As demonstrated by Figure 1, the thickness of the annealing polyanionic brush monotonically increases upon an increase in $pH_b$, which is unambiguously attributed to an increase in the average degree of ionization $\langle\alpha(z)\rangle$ of the monomer units of the brush-forming chains that leads to an increase in the osmotic pressure inside the brush and stronger chain stretching. The local degree of ionization $\alpha(z)$ is an increasing function of the distance $z$ from the grafting surface in accordance with Equation (20). For the given choice of the brush parameters and salt concentration, the results demonstrate excellent agreement between the analytical theory (solid lines) and SF-SCF calculations (dashed lines) for the annealing brush with average degrees of ionization $\langle\alpha(z)\rangle$ between 0.1 and 0.4.

We remind readers that the analytical theory takes into account only electrostatic intermolecular interactions in the brush and assumes linear (Gaussian) elasticity of the brush-forming chains, while numerical SF-SCF modelling accounts for both electrostatic and excluded volume interactions and for the finite chains’ extensibility as well. Therefore,
certain discrepancies between the results of the analytical theory and numerical modelling arise: (i) at $<\alpha(z)> \le 0.1$, when the brush is weakly swollen and the contributions of electrostatic and excluded volume intermolecular interactions are comparable; and (ii) at $<\alpha(z)> \ge 0.4$, due to the non-linear elastic response of strongly stretched chains (underestimated in the analytical theory). Comparison of the electrostatic potential profiles generated by quenched and annealing PE brushes with the same $\alpha = <\alpha(z)>$ reveals that the quenched brush generates a stronger electrostatic field, which should promote the absorption of ampholytic nanocolloids. Remarkably, as one can see from Figure 1, PE chains in the annealing brush with $<\alpha(z)> = \alpha$ are more extended than those in the quenched PE brush (the annealing PE brush thickness is larger) because the degree of ionization of monomer units increases as a function of $z$, from the grafting surface $z=0$ towards the edge of the brush $z=H$.

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Electrostatic potential profiles for pH-sensitive (annealing) PE brush at different pH offsets $\Delta pH_{brush} = pH_b - pK_{brush}$ relative to the brush $pK_{brush} = 5.5$. Solid lines correspond to analytical theory predictions, while dashed lines represent numerical results obtained by the Scheutjens-Fleer self-consistent field (SF-SCF) method for the annealing PE brush. Dash-dot lines show electrostatic potential profiles for strong (quenched) PE brushes with average fractions of charged groups $\alpha = <\alpha(z)>$, matching those of the annealing brush at corresponding $pH_b$ values. Other parameters are chain length $N = 400$, grafting area per chain $s = 100$, and salt concentration $c_s = 0.001$ (in dimensionless units), and solvent quality corresponds to the $\theta$-condition ($\chi = 0.5$). Brush boundaries $z = H$ are indicated by arrows in the corresponding colors.</figcaption></figure>

Figure 2 illustrates similarities and differences in the nanocolloid absorption scenario by annealing and quenched PE brushes on the “wrong side” of the IEP, i.e., when the brush and the nanocolloid in the buffer are similarly (negatively, in our case) charged. In Figure 2a, the electrostatic parts $\Delta F_{ion}(z)/N_{\Sigma}k_B T$ of the nanocolloid insertion free energy into an annealing (solid lines) or quenched (dashed lines) PE brushes with a similar average fraction of charged monomer units, $\alpha = <\alpha(z)>$, are presented. The three panels, from left to right, correspond to different values of the nanocolloid isoelectric point $pI = 4.5, 5.5, 6.5$, representing nanocolloids with different fractions and ionization constants ($pK_+, pK_-$) of cationic and anionic groups. The ionization constant of the annealing brush, $pK_{brush} = 5.5$, is the same in all the panels. Curves of different colors correspond to different offsets $\delta pH_{pI} = pH_b - pI$ of the buffer $pH_b$ from the respective IEPs of the nanocolloids. For each value of $pI$, four values of $\delta pH_{pI}$ above the IEP, that is, $\delta pH_{pI} = pH_b - pI = 0, 0.3, 0.5, 0.7$, corresponding to the zero or negative charge of the nanocolloid in the buffer, are considered.
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Ionic contribution, $\Delta F_{ion}(z)/N_{Σ}k_{B}T$, to the nanocolloid insertion free energy into a weak (annealing, solid lines) or strong (quenched, dashed lines) PE brush with $\alpha = <\alpha(z)>$ (a) and position-dependent net nanocolloid charge $Q(z)/N_{Σ}k_{B}T$ and local ionization degree $\alpha(z)$ of PE chains in the annealing PE brush (b) for various values of the $pH_b$ offset above the nanocolloid's isoelectric point, defined as $\delta pH_{pI} = pH_b - pI = 0, 0.3, 0.5, 0.7$, as indicated in the inserts. The three panels (from left to right) correspond to different values of the nanocolloid isoelectric point $pI = 4.5, 5.5, 6.5$, corresponding to different fractions ($f_+$) and ionization constants of cationic and anionic groups ($pK_+, pK_-$), as indicated above the panels. The ionization constant for the annealing PE brush is $pK_{brush} = 5.5$ for all the panels. The position of the brush boundary is indicated by colored arrows in the middle panel. The brush parameters are chain length $N = 400$, grafting density $s = 100$, salt concentration $c_s = 0.001$, Flory-Huggins parameter $\chi = 0.5$ (theta-solvent conditions).</figcaption></figure>

The profiles of the position-dependent nanocolloid charge $Q(z)/N_{Σ}k_{B}T$ and local brush ionization degree $\alpha(z)$ for the same values of $pI$ and $\delta pH_{pI} = pH_b - pI$ are shown in Figure 2b, which illustrates a monotonous increase in $Q(z)$ and a concomitant decrease in $\alpha(z)$ upon approaching the grafting surface.

All three panels in Figure 2a,b demonstrate similar qualitative features in the evolution of $\Delta F_{ion}(z)$ and $Q(z)$ curves upon an increase in $pH_b$ up from the nanocolloid's IEP for both annealing and quenched PE brushes: At $pH_b = pI$, the free energy $\Delta F_{ion}(z)$ monotonously decreases upon a decrease in $z$, from zero at $z = \infty$ to a negative value attained in the edge minimum located at $z = 0$, which points to thermodynamically favorable absorption of nanocolloids in the brush. Simultaneously, the nanocolloid charge monotonously increases from zero far away from the brush to a positive value inside the brush. The edge minimum of $\Delta F_{ion}(z)$ at the grafting surface, $z = 0$, is preserved also in the vicinity of the IEP, $\delta pH_{pI} \ge 0$, when the nanocolloid is negatively charged in the buffer, $Q(z = \infty) \le 0$, but acquires a positive net charge inside the brush. Additionally, a maximum in $\Delta F_{ion}(z)$ emerges in the proximity of the brush edge. The position of the maximum coincides with the nanocolloid charge inversion point, $pH(z) = pI$, as shown in Figure 2b. The height of
the maximum increases, whereas the depth of the minimum decreases upon an increase in $\delta pH_{pI}$. Negative values of $\Delta F_{ion}(z = 0)$ indicate the presence of a thermodynamic driving force for the nanocolloid absorption in the brush, whereas positive values of $\Delta F_{ion}(z = 0)$ in the minimum indicate its metastable character, i.e., inversion of the nanocolloid charge upon insertion into the brush does not imply that the colloid absorption by the brush is thermodynamically favorable. These trends in the nanocolloid interaction with a PE brush as a function of $\delta pH_{pI}$ were previously discussed in [23,25].

The depth of the freeenergy $\Delta F_{ion}$ minimum at $z = 0$ is systematically larger for a quenched than for an annealing PE brush for the same $\alpha = \langle\alpha(z)\rangle$ and the same $\delta pH_{pI}$, which points to more efficient nanocolloid absorption by a strong than by a weak PE brush on the “wrong side” of the IEP. As one can see in Figure 2a, due to the larger thickness of the annealing PE brush than that of the quenched PE brush with the same $\alpha = \langle\alpha(z)\rangle$, the freeenergy maximum at the brush periphery is systematically wider in the annealing PE brush than in the quenched PE brush, whereas the magnitude of the maximum for given $pI$ and $pH_b$ is the same in quenched and annealing brushes, as follows from Equations (4)–(10).

The right panels in Figure 2a,b correspond to the case $pI = 6.5 > pK_{brush}$, when the brush is already strongly charged, $\langle\alpha(z)\rangle = 0.49$ at $pH_b = pI$ (i.e. at $\delta pH_{pI} = 0$), and its charge further increases upon an increase in $pH_b$, reaching $\langle\alpha(z)\rangle = 0.78$ at $\delta pH_{pI} = 0.7$. As a result of the high ionization degree, the brush-forming PE chains are stretched beyond the linear elasticity regime, and therefore, we use the numerical SF-SCF scheme to calculate the electrostatic potential in the brush. Because in a strongly ionized annealing PE brush at $pH_b \ge pK_{brush}$, the local degree of ionization $\alpha(z)$ weakly varies across the brush, the insertion freeenergy $\Delta F_{ion}(z)$ profiles at $pH_b \ge pI$ fairly coincide for quenched and annealing PE brushes with $\alpha = \langle\alpha(z)\rangle$. Remarkably, for the chosen values of $pK_+$ and $pK_-$, the nanocolloid remains very weakly charged in the vicinity of its isoelectric point due to the broad buffering region between $pK_+$ and $pK_-$, as previously discussed in [25]. As a result, the magnitude of variation of $\Delta F_{ion}(z)$ is small, and the driving force for the nanocolloid absorption is weak.

The left panels in Figure 2a,b correspond to the opposite case of $pI = 4.5 < pK_{brush}$, when the annealing PE brush is almost uncharged in the vicinity of the nanocolloid’s IEP ($\langle\alpha\rangle_{pH_b=pI} \approx 0.03$). Therefore, contributions of both ionic and excluded volume interactions between brush-forming chains are essential in determining the grafted chains’ conformations and resulting electrostatic potential profiles, which motivated us to use the numerical SF-SCF approach in this case as well. Since at $pH_b \ge pI$, the brush is only weakly charged and the magnitude of the positive charge acquired by the nanocolloid inside the brush is small (see the left panel in Figure 2b), the driving force for the nanocolloid absorption by the brush is weak. As $\delta pH_{pI}$ increases above the IEP, both the colloid and then (at even larger $pH_b$) the brush acquire negative charge, leading to mutual repulsion and eventual suppression of the driving force for absorption. Here, the difference between $\Delta F_{ion}(z)$ curves for the nanocolloid inserted into quenched or annealing PE brushes is most pronounced in the peripheral regions of the brushes, $z \le H$, with a noticeably wider free energy barrier in the latter case.

The most interesting scenario of the nanocolloid–brush interaction occurs at $pI = pK_{brush} = 5.5$, as illustrated by the middle panels in Figure 2a,b. In this case, both the nanocolloid and the brush rapidly change their charge in response to small $\delta pH_{pI}$ shifts with respect to $pI = pK_{brush}$, resulting in the largest magnitude of variation in the electrostatic free energy $\Delta F_{ion}(z)$ and the most pronounced deviations between $\Delta F_{ion}(z)$ curves for quenched and annealing PE brushes with the same $\alpha = \langle\alpha(z)\rangle$. In particular, this difference is manifested in a significantly wider freeenergy maximum in the case of the annealing PE
brush compared to the quenched one. In the $pH_b$ range considered here, the chains in the PE brushes are moderately stretched (far below the finite extensibility limit), which enabled us to derive the electrostatic potential and insertion free energy analytically (not shown results of the SF-SCF numerical modelling perfectly match the analytical results).

Figure 3 presents the electrostatic freeenergy profiles $\Delta F_{ion}(z)/N_{Σ}k_B T$ of a nanocolloid in an annealing PE brush (left panel) with a varied, $pH_b$-dependent average degree of ionization $\langle a(z) \rangle$ and a quenched PE brush (right panel) with a fixed degree of ionization $\alpha = 0.22$. The ionization constant of the monomer units in the annealing PE brush is set to $pK_{brush} = 5.5$, matching the nanocolloid's isoelectric point, $pI = 5.5$. The freeenergy profiles are shown for several values of the $pH_b$ offset $\delta pH_{pI} = -0.1, -0.3, -0.5, -0.8$, corresponding to acidic conditions ($pH_b < pI = pK_{brush}$) when the nanocolloid is positively charged, while the brush is always charged negatively.

$$pK^{-} = 5.5, pK^{+} = 5.5, f_{+} = 0.5, pI = 5.5$$

<figure><img src="image_6.png" /><figcaption>**Figure 3.** Electrostatic freeenergy $\Delta F_{ion}(z)/N_{Σ}k_B T$ profiles of a nanocolloid in a quenched PE brush with $\alpha = 0.22$ (right panel) and in an annealing PE brush (left panel) at varied offsets $\delta pH_{pI} = -0.1, -0.3, -0.5, -0.8$ (shown by the color code) from the nanocolloid IEP. For the annealing PE brush, variation in $\delta pH_{pI}$ is accompanied by a variation in the average degree of ionization $\langle a(z) \rangle$. The ionization constant of the monomer units of the annealing brush-forming chains is set to $pK_{brush} = 5.5$, matching the nanocolloid's isoelectric point ($pI = 5.5$). In the right panel, solid lines represent analytical theory results, whereas dashed lines correspond to the results obtained by SF-SCF numerical modelling. The brush parameters are $s = 100, N = 400, \chi = 0.5$, and salt concentration $c_s = 0.001$.</figcaption></figure>

Dashed lines in the right panel for the quenched PE brush represent SF-SCF numerical modelling results, which show good agreement with the results of analytical theory (solid lines). For the annealing PE brush (left panel), the freeenergy profiles were obtained by using the electrostatic potential $\psi(z)$, calculated by the SF-SCF numerical method, because the low degree of brush ionization under these conditions requires simultaneous accounting of electrostatic and excluded volume interactions between brush-forming chains.

In the case of the quenched brush, a decrease in $pH_b$ below the IEP leads to progressive protonation of the nanocolloid and an increase in its net positive charge, which enhances electrostatic attraction to the permanently negatively charged PE brush and results in an increasing depth of the freeenergy $\Delta F_{ion}(z)$ minimum inside the brush.

The effect of $pH_b$ on the insertion freeenergy $\Delta F_{ion}(z)$ profiles in the case of an annealing PE brush is more complex. As $pH_b$ decreases, the brush becomes less ionized and, as a result, produces a weaker electrostatic field. In spite of the concomitant increase in the positive charge of the nanocolloid, the balance of these two opposite trends may result in a weakening of the electrostatic driving force for absorption at sufficiently low $pH_b$. This
trend is illustrated in Figure 3 (left panel), where the depth of the free energy minimum, $|\Delta F_{\text{ion}}(z = 0)|$, increases upon a decrease in $pH_b$ from $pH_b = 5.5$ to $pH_b = 5.1$ but then decreases upon further lowering of $pH_b$.

Up to now, we have considered only the ionic contribution, $\Delta F_{\text{ion}}(z)$, to the insertion free energy, with a focus on its dependence on $pH_b$, and on defining the conditions when the electrostatic driving force for the nanocolloid absorption by the brush emerges. However, even if $\Delta F_{\text{ion}}(z)$ provides the driving force for the absorption, it can be counterbalanced by the osmotic contribution to the free energy, defined by Equation (11). The osmotic term in the free energy arises due to the difference in osmotic pressure of mobile ions in the brush interior and in the bulk of the solution, complemented by excluded volume repulsions between monomer units of the brush-forming chains. The osmotic term is always positive and produces the force expelling the nanocolloid from the brush. For a quenched PE brush, this term is $pH_b$-independent, whereas for an annealing PE brush, both the concentration of mobile ions and that of monomer units in the brush depend on its degree of ionization and thus on $pH_b$. Figure 4 shows the profiles of the total insertion free energy $\Delta F(z) = \Delta F_{\text{ion}}(z) + \Delta F_{\text{osm}}(z)$ of a nanocolloid into the annealing PE brush, which includes both electrostatic and osmotic contributions, at three different values of $pI$ and different offsets $\delta pH_{pI} = pH_b - pI$ of the buffer $pH_b$ above the IEP of the nanocolloid. The sets of parameters used are the same as in Figure 2, not volume. The ratio $\omega$ of the nanocolloid volume to the number of ionizable groups is set to $\omega = V/a^3 N_{\Sigma} = 15$.

Comparison of Figures 2 and 4 shows that accounting for the osmotic force expelling the nanocolloid from the brush leads to a systematic shift of the $\Delta F(z) = \Delta F_{\text{ion}}(z) + \Delta F_{\text{osm}}(z)$ curves with respect to the $\Delta F_{\text{ion}}(z)$ curves upwards, while the shape of the curves in the left and central panels is preserved because of the low or moderate ionization of the brush-forming chains in the vicinity of the nanocolloid IEP and the small excess osmotic pressure in the brush.

For the same values of $\delta pH_{pI}$, the addition of a positive $\Delta F_{\text{osm}}(z)$ term to the $\Delta F_{\text{ion}}(z)$ results in a decrease in the depth of the edge minimum at $z = 0$, or, at larger $\delta pH_{pI}$, to the disappearance of the minimum, while the height of the maximum increases. It is worth mentioning that the position of the total free energy maximum does not exactly match the point of the nanocolloid charge inversion.

In contrast, in the right panel, corresponding to $pI > pK_{\text{brush}}$, the brush is substantially ionized at $pH_b \approx pI$, which produces strong excess osmotic pressure inside the brush. This overcompensates the electrostatic driving force for the nanocolloid absorption even in the proximity of the IEP: the total free energy $\Delta F(z)$ monotonously increases upon approaching the grafting surface, and the nanocolloids are expelled from the brush.

In Figure 5, the profiles of the total insertion free energy $\Delta F(z) = \Delta F_{\text{ion}}(z) + \Delta F_{\text{osm}}(z)$ are shown for two values of $pI = 5.5$, 6.5 and different $pH_b$ offsets below the nanocolloid IEPs. In these cases, the nanocolloid is positively charged in the buffer, and the magnitude of its positive charge increases upon insertion into a negatively charged brush, providing an electrostatic driving force for the absorption ($\Delta F_{\text{ion}}(z)$ is negative and monotonously decreases upon approaching the grafting surface). At $pI = 6.5 > pK_{\text{brush}}$ (right panel), the high osmotic pressure in the strongly ionized brush at $pH_b \approx pI$ overperforms the electrostatic attraction and suppresses absorption. A decrease in $pH_b$ down to $\delta pH_{pI} = -1$ leads to a lowering of the freeenergy profile and the appearance of a very shallow freeenergy minimum close to the brush edge due to the increasing positive charge of the nanocolloid. The same trend is even more pronounced at $pI = pK_{\text{brush}}$, when the charge of the brush and, consequently, the osmotic pressure of counterions in the brush rapidly decrease upon lowering $pH_b$ below $pI = pK_{\text{brush}}$: a shallow minimum in the $\Delta F(z)$ curve emerges in the brush periphery even at $\delta pH_{pI} = -0.1$. Upon a further decrease in $pH_b$, this minimum
becomes deeper and is displaced towards the central region of the brush. Hence, due to the interplay of electrostatic attraction and osmotic repulsion at $pH_b$ below the IEP, absorption and accumulation of positively charged nanocolloids occurs in the peripheral and central regions of the negatively charged PE brush, while the nanocolloids are expelled from the proximity of the grafting surface, where the local osmotic pressure is higher.

$$N = 400, S = 100, C_s = 0.001, \chi = 0.5, pK_{\text{brush}} = 5.5, \omega = 15$$

<figure><img src="image_7.png" /><figcaption>**Figure 4.** Total free energy $\Delta F(z)/N_Σk_B T = (\Delta F_{\text{ion}}(z) + \Delta F_{\text{osm}}(z))/N_Σk_B T$ of the nanocolloid insertion into an annealing PE (polyanionic) brush for various values of the $pH_b$ offset above the nanocolloid's isoelectric point, $\delta pH_{pI} = pH_b - pI = 0, 0.3, 0.5, 0.7$, as indicated in the inserts. The three panels (from left to right) correspond to different values of the nanocolloid isoelectric point $pI = 4.5, 5.5, 6.5$, corresponding to different fractions, $f_+$, and ionization constants of cationic and anionic groups ($pK_+, pK_-$), as indicated above the panels. The ionization constant for the annealing PE brush is $pK_{\text{brush}} = 5.5$. The position of the brush boundary is indicated by colored arrows in the middle panel. The ratio $\omega$ of the nanocolloid volume to the number of ionizable groups is set to $\omega = V/a^3 N_Σ = 15$. The brush parameters are chain length $N = 400$, grafting density $s = 100$, salt concentration $c_s = 0.001$, and Flory-Huggins parameter $\chi = 0.5$ (theta solvent conditions).</figcaption></figure>

$$N = 400, S = 100, C_s = 0.01, \chi = 0.5, pK_{\text{brush}} = 5.5, \omega = 50$$

<figure><img src="image_8.png" /><figcaption>**Figure 5.** Total free energy $\Delta F(z)/N_Σk_B T = (\Delta F_{\text{ion}}(z) + \Delta F_{\text{osm}}(z))/N_Σk_B T$ of the nanocolloid insertion into an annealing PE (polyanionic) brush for various values of the $pH_b$ offset below the nanocolloid's isoelectric point, $\delta pH_{pI} = pH_b - pI = -0.1, -0.5, -1.0$, as indicated in the inserts. The two panels (from left to right) correspond to different values of the nanocolloid isoelectric point $pI = 5.5, 6.5$, corresponding to different fractions ($f_+$) and ionization constants of cationic and anionic groups ($pK_+, pK_-$), as indicated above the panels. The ionization constant for the annealing PE brush is $pK_{\text{brush}} = 5.5$. The ratio $\omega$ of the nanocolloid volume to the number of ionizable groups is set to $\omega = V/a^3 N_Σ = 50$. The brush parameters are: chain length $N = 400$, grafting density $s = 100$, salt concentration $c_s = 0.01$, and Flory-Huggins parameter $\chi = 0.5$ (theta-solvent conditions).</figcaption></figure>
## 3. Materials and Methods

### 3.1. Nanocolloid Charge and Insertion Free Energy

We consider a protein-like ampholytic nanocolloidal particle of volume V, comprising $N_+$ cationic and $N_-$ anionic ionizable groups on its surface, with respective acidic ionization constants $K_+$ and $K_-$. The total number of ionizable groups on the particle surface is $N_Σ = N_+ + N_-$. The fraction of cationic groups is $f_+ = N_+ / N_Σ$. The particle is immersed in an aqueous solution with a fixed $pH_b$. (Here and below, the subscript “b” refers to the bulk of the solution.) The solution contains monovalent cations and anions of added monovalent salt, with respective concentrations $c_+ = c_- = c_s$, which specify the Debye screening length in the bulk of the solution as $κ^{-1} = (8πl_Bc_s)^{-1/2}$.

The net charge of the particle in the bulk of the solution is given by

$$Q_b = \alpha_{b+}N_+ - \alpha_{b-}N_- \quad (1)$$

with

$$\alpha_{b+} = (1 + 10^{pH_b - pK_+})^{-1} \quad (2)$$

$$\alpha_{b-} = (1 + 10^{pK_- - pH_b})^{-1} \quad (3)$$

being the respective ionization degrees of the particle’s cationic and anionic residues in the bulk of the solution.

The PE chains tethered to the planar surface give rise to an electrical field with an electrostatic potential $Ψ(z)$, depending on the distance $z$ from the grafting surface, and calibrated as $Ψ(z = ∞) = 0$. The local concentration $[H^+(z)]$ of hydrogen ions in the brush follows the Boltzmann law,

$$H^+(z)] = [H_b^+] \exp(-\psi(z)) \quad (4)$$

with $\psi(z) \equiv eΨ(z) / k_B T$, which leads to the position-dependent

$$pH(z) = pH_b + \psi(z) / \ln 10 \quad (5)$$

and the position-dependent nanocolloid charge

$$Q(z) = \alpha_+(z)N_+ - \alpha_-(z)N_- \quad (6)$$

with

$$\alpha_+(z) = (1 + 10^{pH(z) - pK_+})^{-1} \quad (7)$$

$$\alpha_-(z) = (1 + 10^{pK_- - pH(z)})^{-1} \quad (8)$$

The position-dependent free energy of the nanocolloid is given by

$$\Delta F(z) = \Delta F_{ion}(z) + \Delta F_{osm}(z) \quad (9)$$

The first term in Equation (9) accounts for the ionization free energy of cationic and anionic groups in an external electrical field $\psi(z)$, which can be presented [22,25] as

$$\Delta F_{ion}(z) / k_B T = N_+ \ln \left( \frac{1 - \alpha_+(z)}{1 - \alpha_{b+}} \right) + N_- \ln \left( \frac{1 - \alpha_-(z)}{1 - \alpha_{b-}} \right) \quad (10)$$

with the reference state $\Delta F_{ion}(z = ∞) = 0$. The free energy given by Equation (10) can be either positive or negative, depending on the degrees of ionization, $\alpha_-(z)$, $\alpha_+(z)$ of anionic
and cationic groups of the colloid, as compared to their ionization degrees, $\alpha_{b+}$ and $\alpha_{b-}$, in the bulk of the solution.

The second term in Equation (9) accounts for the work performed against the excess osmotic pressure upon the nanocolloid insertion in the brush.

$$ \Delta F_{osm} = \Delta\Pi(z) \cdot V \qquad (11) $$

The differential osmotic pressure

$$ \Delta\Pi(z) = \Delta\Pi_{ion}(z) + \Delta\Pi_{pol}(z) \qquad (12) $$

comprises the ideal gas contribution due to mobile ions,

$$ \Delta\Pi_{ion}(z)/k_{B}T = 4c_{s}\sinh^{2}(\frac{\psi(z)}{2}) \qquad (13) $$

and the excluded volume contribution due to polymer chains with volume fraction $\phi(z)$,

$$ \Delta\Pi_{pol}(z)/k_{B}T = -\ln(1 - \phi(z)) - \chi\phi^{2}(z) - \phi(z) \qquad (14) $$

where $\chi$ is the Flory–Huggins parameter.

Remarkably, $\Delta\Pi_{ion}(z)$ and $\Delta F_{ion}(z)$ are non-vanishing beyond the edge of the brush due to the continuity of the electrostatic potential $\psi(z)$. Typically, the first term in Equation (12) dominates over the second one. However, the second term is important for particles embedded in a weak polyanionic brush with $pH < pK_{brush}$.

Here, we introduce the parameter $\omega$, defined as the average nanocolloid volume normalized by the volume of a polymer segment and the number of ionizable groups,

$$ \omega = \frac{V}{N_{\Sigma}a^{3}} $$

The relative magnitude of the osmotic and ionic contributions, $\frac{\Delta F_{osm}}{\Delta F_{ion}}$, depends on the parameter $\omega$, since $\Delta F_{ion} \sim N_{\Sigma}$ whereas $\Delta F_{osm} \sim V$.

## 3.2. Electrostatic Field in Weak and Strong Polyelectrolyte Brush: Analytical Poisson–Boltzmann Approach

To specify the electrostatic potential $\psi(z)$, we consider a brush of ionizable PE chains with a degree of polymerization $N$, as shown in Figure 1. The grafting density is $\sigma = a^2/s$, where $s$ is the area per chain on the grafting surface. The PE chains are assumed to be intrinsically flexible, that is, the monomer unit size $a$ coincides with the statistical segment length of the corresponding uncharged polymer. The latter is assumed to be on the order of the Bjerrum length $l_B = e^2/(\epsilon k_B T)$, which equals approximately 7 nm in water under ambient conditions.

In this section, we compare PE brushes formed by strong (quenched) and weak (annealing or pH-sensitive) anionic polyelectrolytes.

Strong PEs possess a quenched (negative) fractional charge $\alpha e$ per monomer unit (with $e$ as the elementary charge). In a weak polyacidic brush, every monomer unit is capable of acquiring an elementary negative charge via deprotonation. The fraction $\alpha$ of ionized (deprotonated) monomer units in a weak polyelectrolyte brush depends on $pH_b$ in the bulk of the solution and on the local excess electrostatic potential $\psi(z)$ in the PE brush. The maximal value of $\alpha = \alpha_b$ is attained infinitely far away from the brush, $z \to \infty$ (where $\psi(z) \to 0$), and is given by

$$ \alpha_b = (1 + 10^{pK_{brush} - pH_b})^{-1} \qquad (15) $$
Here, $K_{\text{brush}}$ is the acidic ionization constant of monomer units in the brush-forming polyacidic chains. A generalization of the theory for PE brushes of positively charged polyelectrolytes is straightforward. Below, we express all lengths in units of $a$, that is, $l_B/a \to l_B$, $s/a^2 \to s$, and $c_s$ is the dimensionless volume fraction of salt.

To calculate the electrostatic potential $\psi(z)$ created by a weak or strong polyelectrolyte brush, we employ the strong-stretching self-consistent field framework [35], which assumes the parabolic shape of the monomer self-consistent molecular potential $U(z)$ inside the brush:

$$ \frac{U(z)}{k_B T} = -\text{const} - \frac{3\pi^2}{8N^2} z^2, \quad z \le H \qquad (16) $$

Here, $z$ is the distance from the grafting surface, and $H$ is the cut-off of the polymer concentration profile (the brush thickness). The constant term in Equation (16) depends on the details of the interactions in the system. Equation (16) is valid for planar or concave brushes of monodisperse flexible chains stretched within the limits of linear (Gaussian) conformational elasticity.

In strong (quenched) PE brushes (subscript “strong”),

$$ \frac{U_{\text{strong}}(z)}{k_B T} = \alpha \psi_{\text{in}}(z), \quad z \le H \qquad (17) $$

leading (together with Equation (16)) to the explicit profile of the electrostatic potential,

$$ \psi_{\text{in}}(z) = \frac{z^2 - H^2}{H_0^2} + \psi_{\text{in}}(H), \quad z \le H \qquad (18) $$

Here, $H$ is the total brush thickness, $H_0$ is the characteristic length given by

$$ H_0 = \sqrt{\frac{8}{3\pi^2} N \alpha^{1/2} a} \qquad (19) $$

and $\psi_{\text{in}}(H)$ should be found from the condition of continuity of the electrostatic potential at the brush edge, $z = H$.

In a weak (pH-sensitive) PE brush, the ionization degree $\alpha(z)$ of a monomer unit depends on its position $z$ in the brush as

$$ \alpha(z) = (1 + 10^{pK_{\text{brush}} - pH(z)})^{-1} = [1 + \frac{1 - \alpha_b}{\alpha_b} \exp(-\psi(z))]^{-1} \qquad (20) $$

with Equation (15) specifying $\alpha_b$, and Equation (5) for $pH(z)$.

In this case, the molecular potential $U_{\text{weak}}(z)$ is related to $\alpha(z)$ in Equation (20) as [31]

$$ \frac{U_{\text{weak}}(z)}{k_B T} = \ln(1 - \alpha(z)) \qquad (21) $$

By combining Equations (16) and (21), one finds an explicit expression for the ionization profile in the weak PE brush as

$$ \alpha(z) = 1 - (1 - \alpha_H) \exp\left[\frac{\alpha_b(H^2 - z^2)}{H_{0b}^2}\right] \qquad (22) $$

with

$$ H_{0b} = \sqrt{\frac{8}{3\pi^2} N \alpha_b^{1/2}} \qquad (23) $$

and $\alpha_H \le \alpha_b$ being the degree of ionization at the edge of the brush, i.e., at $z = H$.
The value of $\alpha_H$ depends in a complex way on the following parameters of the system: $\alpha_b$ (controlled by pH, see Equation (15)), salt concentration $c_s$, chain length $N$, and grafting density $a^2/s$. As follows from Equation (22), the degree of ionization $\alpha(z)$ inside the brush monotonously increases as a function of $z$, i.e., with increasing distance from the grafting surface.

By combining Equations (20) and (22), we obtain the electrostatic potential profile, $\psi_{in}(z)$, inside a weak polyacid brush as

$$ \psi_{in}(z) = -\frac{\alpha_b(H^2 - z^2)}{H_0^2} + \ln\left\{\frac{1 - \alpha_b}{\alpha_b(1 - \alpha_H)}\left(1 - (1 - \alpha_H) \exp\left[\frac{\alpha_b(H^2 - z^2)}{H_0^2}\right]\right)\right\} \quad (24) $$

with

$$ \alpha_H = \left[1 + \frac{1 - \alpha_b}{\alpha_b} \exp(-\psi_{in}(H))\right]^{-1} \quad (25) $$

As follows from Equation (24), the electrostatic potential $\psi_{in}(z)$ inside a weak polyacidic brush follows a more complex functional dependence on $z$ than in a strong (quenched) polyelectrolyte brush.

The value of $\psi_{in}(H)$, or, equivalently, in the case of a weak polyelectrolyte brush, $\alpha(H)$, which is related to $\psi(H)$ by Equation (25), has to be specified from the condition of continuity of the electrostatic potential at the brush edge, $z = H$, i.e., $\psi_{in}(H) = \psi_{out}(H)$.

By assimilating the outer edge of the brush to a uniformly charged plane with charge

$$ \tilde{Q} = \int_{0}^{H} \rho_{in}(z) dz $$

per unit area, with

$$ \rho(z) = -\frac{1}{4\pi l_B} \frac{d^2\psi_{in}(z)}{dz^2} \quad (26) $$

being the net charge density inside the brush, the distribution of the electrostatic potential outside the brush, at $z \ge H$, can be presented as

$$ \psi_{out}(z) = -2 \ln \left[ \frac{(\kappa \tilde{\Lambda} + \sqrt{(\kappa \tilde{\Lambda})^2 + 1} - 1) + (\kappa \tilde{\Lambda} - \sqrt{(\kappa \tilde{\Lambda})^2 + 1} + 1) e^{-\kappa(z-H)}}{(\kappa \tilde{\Lambda} + \sqrt{(\kappa \tilde{\Lambda})^2 + 1} - 1) - (\kappa \tilde{\Lambda} - \sqrt{(\kappa \tilde{\Lambda})^2 + 1} + 1) e^{-\kappa(z-H)}} \right] \quad (27) $$

and

$$ \psi_{out}(H) = 2 \ln \left[ \frac{\sqrt{(\kappa \tilde{\Lambda})^2 + 1} - 1}{\kappa \tilde{\Lambda}} \right] \quad (28) $$

Here,

$$ \tilde{\Lambda} = \frac{1}{2\pi l_B |\tilde{Q}|} = \begin{cases} H_0^2/H^2 & \alpha = \text{const} \\ \alpha_H H_{0b}^2 / \alpha_b H & \alpha = \alpha(z) \end{cases} \quad (29) $$

is the Gouy-Chapman length outside strong (the upper line) and weak (the lower line) polyelectrolyte brushes.

By using the condition $\psi_{out}(H) = \psi_{in}(H)$ and substituting $\psi(H)$ given by Equation (28) into Equation (25), we obtain a closed equation for $\alpha_H$ in a weak PE brush as a function of $H$, $\kappa$, $H_0$, and $\alpha_b$ in the form

$$ \alpha_H = \left( 1 + \frac{1 - \alpha_b}{\alpha_b} \left[ \frac{\kappa \tilde{\Lambda}}{\sqrt{(\kappa \tilde{\Lambda})^2 + 1 - 1}} \right]^2 \right)^{-1} \quad (30) $$
Because mobile cations and anions inside the brush are distributed according to the Boltzmann law,

$$c_{+,in}(z) = c_s \exp(-\psi_{in}(z))$$

$$c_{-,in}(z) = c_s \exp(\psi_{in}(z))$$

the polymer concentration profile inside the brush can be found as

$$c_p(z) = \alpha^{-1}(z) [-\rho_{in}(z) - c_{-,in}(z) + c_{+,in}(z)] = \alpha^{-1}(z) [-\rho_{in}(z) - 2c_s \sinh(\psi_{in}(z))]$$

with the net charge density $\rho(z)$ inside the brush specified by Equation (26), and in the case of a strong PE brush, $\alpha(z) \equiv \alpha$.

The result is

$$c_p(z) = \frac{1}{2\pi\alpha l_B H_0^2} \left[ 1 - \kappa^2 H_0^2 \sinh \left( \frac{z^2 - H^2}{H_0^2} + 2 \ln \left[ \frac{\sqrt{(\kappa H_0^2/H)^2 + 1 - 1}}{\kappa H_0^2/H} \right] \right) \right] \quad (31)$$

for a strong PE brush, and

$$c_p(z) = \frac{\alpha_b}{2\pi l_B H_0^2} \left[ \frac{1 - (1 - \alpha_H)(1 + 2\alpha_b z^2/H_0^2) \exp[\alpha_b(H^2 - z^2)/H_0^2]}{\{1 - (1 - \alpha_H) \exp[\alpha_b(H^2 - z^2)/H_0^2]\}^3} + \right] \quad (32)$$

$$\frac{\kappa^2 H_0^2}{4\alpha_b} \left\{ \frac{\alpha_b}{1 - \alpha_b} \cdot \frac{(1 - \alpha_H) \exp[\alpha_b(H^2 - z^2)/H_0^2]}{\{1 - (1 - \alpha_H) \exp[\alpha_b(H^2 - z^2)/H_0^2]\}^2} - \frac{1 - \alpha_b}{\alpha_b} \cdot \frac{\exp[-\alpha_b(H^2 - z^2)/H_0^2]}{(1 - \alpha_H)} \right\} \quad (33)$$

for a weak PE brush.

The average degree of ionization in a weak PE brush can be calculated as

$$\langle \alpha(z) \rangle = \frac{s}{N} \int_{0}^{H} \alpha(z) c_{p}(z) dz \quad (34)$$

Finally, the brush thickness $H$ can be found from the normalization condition

$$\int_{0}^{H} c_{p}(z) dz = \frac{N}{s} \quad (35)$$

## 4. Conclusions

The interaction of ampholytic protein-like nanocolloids with polyelectrolyte brushes formed by either weak (pH-sensitive) or strong PE chains, end-tethered to a planar surface, and the equilibrium partitioning of nanocolloids between the buffer solution and the brush have been studied using an analytical self-consistent field theory that implements the Poisson-Boltzmann and strong stretching approximations. The electrostatic potential, the position-dependent particle charge, and the insertion free energy were calculated as functions of the particle and brush parameters, as well as the environmental conditions (pH and ionic strength in the solution). Under acidic conditions, when the pH-sensitive polyanionic brush is weakly charged, the electrostatic potential and polymer concentration profiles were calculated using the numerical Scheutjens-Fleer scheme.

The developed theory proves that both strong and pH-sensitive PE brushes are capable of nanocolloid uptake on the “wrong side” of its isoelectric point, that is, when the charge on the colloid particle is of the same sign as the ionized monomer units in the brush-forming chains. Such uptake may be thermodynamically favorable due to the re-ionization of pH-sensitive residues that leads to the inversion of the sign of the particle charge.
However, charge inversion requires overcoming a potential barrier upon particle insertion into the brush. This barrier may sufficiently slow down or even block the particle absorption by a similarly charged brush in the vicinity of the IEP. In the case of $pI \ge pK_{brush}$, the absorption patterns at $pH \ge pI$ are qualitatively similar for strong and weak PE brushes, but the insertion free energy barrier at the brush periphery is wider for a weak PE brush.

A decrease in the buffer $pH_b$ below the particle IEP leads to a monotonic increase in the magnitude of the positive charge on the particle, which enhances absorption by a strong (permanently negatively charged) brush. However, in the case of a pH-sensitive brush, this trend is counterbalanced by the decreasing degree of brush ionization and the concomitant decrease in the magnitude of the electrostatic potential inside the brush, which weakens and eventually suppresses absorption.

**Author Contributions:** T.O.P. investigation, article writing; E.B.Z. conceptualization, investigation; O.V.B. conceptualization, article writing. All authors have read and agreed to the published version of the manuscript.

**Funding:** Russian Science Foundation grant 23-13-00174.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** Data are contained within the article.

**Acknowledgments:** This work was supported by the Russian Science Foundation, grant 23-13-00174.

**Conflicts of Interest:** The authors declare no conflicts of interest.

## References

1. Dubin, P.; Bock, J.; Davis, R.; Schulz, D.N.; Thies, C. *Macromolecular Complexes in Chemistry and Biology*; Springer: New York, NY, USA, 1994.

2. Ruhe, J.; Ballauff, M.; Biesalski, M.; Dziezok, P.; Gruhn, F.; Johannsmann, D.; Houbenov, N.; Hugenberg, N.; Konradi, R.; Minko, S.; et al. Polyelectrolyte Brushes. *Adv. Polym. Sci.* **2004**, *165*, 79–150.

3. Chen, W.-L.; Cordero, R.; Tran H.; Ober, C.K. 50th Anniversary Perspective: Polymer Brushes: Novel Surfaces for Future Materials. *Macromolecules* **2017**, *50*, 4089–4113. [CrossRef]

4. Willott, J.D.; Murdoch, T.J.; Webber G.B.; Wanless, E.J. Physicochemical behaviour of cationic polyelectrolyte brushes. *Prog. Polym. Sci.* **2017**, *64*, 52–75. [CrossRef]

5. Toomey, R.; Tirrell, M. Functional Polymer Brushes in Aqueous Media from Self-Assembled and Surface-Initiated Polymers. *Ann. Rev. Phys. Chem.* **2008**, *59*, 493–517. [CrossRef] [PubMed]

6. Raviv, U.; Giasson, S.; Kampf, N.; Gohy, J.-F.; Jerome R.; Klein, J. Lubrication by charged polymers. *Nature* **2003**, *425*, 163–165. [CrossRef] [PubMed]

7. Ballauff, M.; Borisov, O. Polyelectrolyte brushes. *Curr. Opin. Colloid Interface Sci.* **2006**, *11*, 316–323. [CrossRef]

8. Wittemann, A.; Haupt, B.; Ballauff, M. Adsorption of proteins on spherical polyelectrolyte brushes in aqueous solution. *Phys. Chem. Chem. Phys.* **2003**, *5*, 1671–1677. [CrossRef]

9. Wittemann, A.; Ballauff, M. Interaction of proteins with linear polyelectrolytes and spherical polyelectrolyte brushes in aqueous solution. *Phys. Chem. Chem. Phys.* **2006**, *8*, 5269–5275. [CrossRef]

10. Becker, A.L.; Henzler, K.; Welsch, N.; Ballauff, M.; Borisov, O.V. Proteins and polyelectrolytes: A charge relationship. *Curr. Opinion Colloid Interface Sci.* **2012**, *17*, 90–96. [CrossRef]

11. Walkowiak, J.; Gradzielski, M.; Zauscher, S.; Ballauff, M. Interaction of Proteins with a Planar Poly(acrylic acid) Brush: Analysis by Quartz Crystal Microbalance with Dissipation Monitoring (QCM-D). *Polymers* **2021**, *13*, 122. [CrossRef]

12. Achazi, K.; Haag, R.; Ballauff, M.; Dernedde, J.; Kizhakkedathu, J.N.; Maysinger, D.; Multhaup, G. Understanding the interaction of polyelectrolyte architectures with proteins and biosystems. *Angew. Chem. Int. Ed.* **2021**, *60*, 3882–3904. [CrossRef]

13. Rosenfeldt, S.; Wittemann, A.; Ballauff, M.; Breininger, E.; Bolze, J.; Dingenouts, N. Interaction of proteins with spherical polyelectrolyte brushes in solution as studied by small-angle X-ray scattering. *Phys. Rev. E* **2004**, *70*, 061403. [CrossRef]

14. Anikin, K.; Röcker, C.; Wittemann, A.; Wiedenmann, J.; Ballauff, M.; Nienhaus, G.U. Polyelectrolyte-mediated protein adsorption: Fluorescent protein binding to individual polyelectrolyte nanospheres. *J. Phys. Chem. B* **2005**, *109*, 5418–5420. [CrossRef]
15. Hollmann, O.; Czeslik, C. Characterization of a planar poly (acrylic acid) brush as a materials coating for controlled protein immobilization. *Langmuir* **2006**, *22*, 3300–3305. [CrossRef]

16. Wittemann, A.; Ballauff, M. Secondary structure analysis of proteins embedded in spherical polyelectrolyte brushes by FT-IR spectroscopy. *Anal. Chem.* **2004**, *76*, 2813–2819. [CrossRef] [PubMed]

17. Leermakers, F.A.M.; Ballauff, M.; Borisov, O.V. On the mechanisms of interaction of globular proteins with polyelectrolyte brushes. *Langmuir* **2007**, *23*, 237–247.

18. Xu, X.; Angioletti-Uberti, S.; Lu, Y.; Dzubiella, J.; Ballauff, M. Interaction of proteins with polyelectrolytes: Comparison of theory to experiment. *Langmuir* **2018**, *35*, 5373–5391. [CrossRef]

19. Kim, S.; Sureka, H.V.; Kayitmazer, A.B.; Wang, G.; Swan, J.W.; Olsen, B.D. Effect of protein surface charge distribution on protein-polyelectrolyte complexation. *Biomacromolecules* **2020**, *21*, 3026–3037. [CrossRef]

20. Biesheuvel, P.M.; Wittemann, A. A modified box model including charged polyelectrolyte brush: Self-consistent field theory. *J. Phys. Chem. B* **2005**, *109*, 4209–4214. [CrossRef]

21. De Vos, W.M.; Leermakers, F.A.M.; De Keizer, A.; Cohen Stuart, M.A.; Klein, J.M. Field theoretical analysis of driving forces for the uptake of proteins by like-charged polyelectrolyte brushes: Effects of charge regulation and patchiness. *Langmuir* **2010**, *26*, 249–259. [CrossRef]

22. Laktionov, M.Y.; Zhulina, E.B.; Borisov, O.V. Proteins and polyampholytes interacting with polyelectrolyte brushes and microgels: The charge reversal concept revised. *Langmuir* **2021**, *37*, 2865–2873. [CrossRef]

23. Salamatova, T.O.; Laktionov, M.Y.; Zhulina, E.B.; Borisov, O.V. Polyelectrolyte Brush Interacting with Ampholytic Nanoparticles as Coarse-Grained Models of Globular Proteins: Poisson–Boltzmann Theory. *Biomacromolecules* **2023**, *24*, 2433–2446. [CrossRef] [PubMed]

24. Salamatova, T.O.; Zhulina, E.B.; Borisov, O.V. Bovine Serum Albumin interaction with polyanionic and polycationic brushes: The case theoretical study. *Int. J. Mol. Sci.* **2023**, *24*, 3395. [CrossRef] [PubMed]

25. Popova, T.O.; Zhulina, E.B.; Borisov, O.V. Polyelectrolyte brushes with protein-like nanocolloids. *Langmuir* **2024**, *40*, 1232–1246. [CrossRef] [PubMed]

26. Yoon, J.Y.; Kim, J.H.; Kim, W.S. The relationship of interaction forces in the protein adsorption onto polymeric microspheres. *Colloids Surfaces Physicochem. Eng. Asp.* **1999**, *153*, 413–419. [CrossRef]

27. Guo, X.; Ballauff, M. Spherical polyelectrolyte brushes: Comparison between annealing and quenched brushes. *Phys. Rev. E* **2001**, *164*, 051406. [CrossRef]

28. Wang, W.; Li, L.; Henzler, K.; Lu, Y.; Wang, J.; Han, H.; Tian, Y.; Wang, Y.; Zhou, Z.; Lotze, G.; et al. Protein immobilization onto cationic spherical polyelectrolyte brushes studied by small angle X-ray scattering. *Biomacromolecules* **2017**, *18*, 1574–1581. [CrossRef]

29. Wang, S.; Chen, K.; Xu, Y.; Yu, X.; Wang, W.; Li, L.; Guo, X. Protein immobilization and separation using anionic/cationic spherical polyelectrolyte brushes based on charge anisotropy. *Soft Matter* **2013**, *9*, 11276–11287. [CrossRef]

30. Israels, R.; Leermakers, F.A.M.; Fleer, G.J. On the Theory of Grafted Weak Polyacids. *Macromolecules* **1994**, *27*, 3087–3093. [CrossRef]

31. Zhulina, E.B.; Borisov, O.V. Poisson-Boltzmann Theory of pH-Sensitive (Annealing) Polyelectrolyte Brush. *Langmuir* **2011**, *27*, 10615–10633. [CrossRef]

32. Zhulina, E.B.; Borisov, O.V. Structure and Interactions of Weakly Charged Polyelectrolyte Brushes: Self-Consistent Field Theory. *J. Chem. Phys.* **1997**, *107*, 5952–5967. [CrossRef]

33. Zhulina, E.B.; Klein Wolterink, J.; Borisov, O.V. Screening Effects in Polyelectrolyte Brush: Self-Consistent Field Theory. *Macromolecules* **2000**, *33*, 4945–4953. [CrossRef]

34. Fleer, G.J.; Cohen Stuart, M.A.; Scheutjens, J.M.H.M.; Cosgrove, T.; Vincent, B. *Polymers at Interfaces*; Chapman & Hall: London, UK, 1993.

35. Semenov, A.N. Contribution to the Theory of Microphase Layering in Block-Copolymer Melts. *Sov. Phys. JETP* **1985**, *61*, 733–742.

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.