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.

The electrostatic potential inside and around α-lactalbumin: Fluctuations and mean-field models.

The electrostatic potential (EP) generated by the protein α-lactoalbumin in the presence of added salt is computed as a thermal average at a given point in space. With this aim, constant pH Monte Carlo (MC) simulations are performed within the primitive model, namely, the solvent is treated as a continuum dielectric. The study of the thermal and spatial fluctuations of the EP reveals that they are remarkably high inside the protein. The calculations indicate that fluctuations inside the protein are mainly due to the asymmetric distribution of the charge groups, while the charge fluctuations of the titratable groups play a minor role. The computed EP matches very well with the one obtained from the Poisson equation for the average charge density in spherical symmetry. The Tanford-Kirkwood multipole expansion reproduces the simulated angular-averaged potential rather accurately. Surprisingly, two of the simplest mean-field models, the linear Poisson-Boltzmann (PB) equation and Donnan potential, provide good estimations of the average EP in the effective protein surface (surface EP). The linear PB equation predicts a linear relationship between charge and surface EP, which is numerically reproduced only if the small ions within the protein are taken into account. On the other hand, the partition coefficients of the small ions inside and outside the protein predicted by Donnan theory reproduce reasonably well the simulation results.

Adsorption of flexible proteins in the 'wrong side' of the isoelectric point: Casein macropeptide as a model system.

We analyze the conditions of the adsorption of a flexible peptide onto a charged substrate in the 'wrong side' of the isoelectric point (WSIP), i.e. when surface and peptide charges have the same sign. As a model system, we focus on the casein macropeptide (CMP), both in the aglycosylated (aCMP) and fully glycosydated (gCMP) forms. We model the substrate as a uniformly charged plane while CMP is treated as a bead-and-spring model including electrostatic interactions, excluded volume effects and acid/base equilibria. Adsorption coverage, aminoacid charges and concentration profiles are computed by means of Monte Carlo simulations at fixed pH and salt concentration. We conclude that for different reasons the CMP can be adsorbed to both positively and negatively charged surfaces in the WSIP. For negatively charged surfaces, WSIP adsorption is due to the patchy distribution of charges: the peptide is attached to the surface by the positively charged end of the chain, while the repulsion of the surface for the negatively charged tail is screened by the small ions of the added salt. This effect increases with salt concentration. Conversely, a positively charged substrate induces strong charge regulation of the peptide: the acidic groups are deprotonated, and the peptide becomes negatively charged. This effect is stronger at low salt concentrations and it is more intense for gCMP than for aCMP, due to the presence of the additional sialic groups in gCMP.

Influence of macromolecular crowding on the charge regulation of intrinsically disordered proteins.

In this work we study the coupling between ionization and conformational properties of two IDPs, histatin-5 and ß-amyloid 42, in the presence of neutral and charged crowders. The latter is modeled to resemble bovine serum albumin (BSA). With this aim, semi-grand canonical Monte Carlo simulations are performed, so that the IDP charge is a dynamic property, undergoing protonation/deprotonation processes. Both ionization properties (global and specific amino acid charge and binding capacitance) and radius of gyration are analyzed in a large range of pH values and salt concentrations. Without crowder agents, the titration curve of histatin-5, a polycation, is salt-dependent while that of ß-amyloid 42, a polyampholyte, is almost unaffected. The salt concentration is found to be particularly relevant at pH values where the protein binding capacitance (directly linked with charge fluctuation) is larger. Upon addition of neutral crowders, charge regulation is observed in histatin-5, while for ß-amyloid 42 this effect is very small. The main mechanism for charge regulation is found to be the effective increase in the ionic strength due to the excluded volume. In the presence of charged crowders, a significant increase in the charge of both IDPs is observed in almost all the pH range. In this case, the IDP charge is altered not only by the increase in the effective ionic strength but also by its direct electrostatic interaction with the charged crowders.

Role of Charge Regulation and Fluctuations in the Conformational and Mechanical Properties of Weak Flexible Polyelectrolytes.

This work addresses the role of charge regulation (CR) and the associated fluctuations in the conformational and mechanical properties of weak polyelectrolytes. Due to CR, changes in the pH-value modifies the average macromolecular charge and conformational equilibria. A second effect is that, for a given average charge per site, fluctuations can alter the intensity of the interactions by means of correlation between binding sites. We investigate both effects by means of Monte Carlo simulations at constant pH-value, so that the charge is a fluctuating quantity. Once the average charge per site is available, we turn off the fluctuations by assigning the same average charge to every site. A constant charge MC simulation is then performed. We make use of a model which accounts for the main fundamental aspects of a linear flexible polyelectrolyte that is, proton binding, angle internal rotation, bond stretching and bending. Steric excluded volume and differentiated treatment for short-range and long-range interactions are also included. This model can be regarded as a kind of "minimal" in the sense that it contains a minimum number of parameters but still preserving the atomistic detail. It is shown that, if fluctuations are activated, gauche state bond probabilities increase and the persistence length decreases, so that the polymer becomes more folded. Macromolecular stretching is also analyzed in presence of CR (the charge depends on the applied force) and without CR (the charge is fixed to the value at zero force). The analysis of the low force scaling behavior concludes that Pincus exponent becomes pH-dependent. Both, with and without CR, a transition from 1/2 at high pH-values (phantom chain) to 3/5 at low pH-values (Pincus regime) is observed. Finally, the intermediate force stretching regime is investigated. It is found that CR induces a moderate influence in the force-extension curves and persistence length (which in this force regime becomes force-dependent). It is thus concluded that the effect of CR on the stretching curves is mainly due to the changes in the average charge at zero force. It is also found that, for the cases studied, the effect of steric excluded volume is almost irrelevant compared to electrostatic interactions.

Coupling of Charge Regulation and Conformational Equilibria in Linear Weak Polyelectrolytes: Treatment of Long-Range Interactions via Effective Short-Ranged and pH-Dependent Interaction Parameters.

The classical Rotational Isomeric State (RIS) model, originally proposed by Flory, has been used to rationalize a wide range of physicochemical properties of neutral polymers. However, many weak polyelectrolytes of interest are able to regulate their charge depending on the conformational state of the bonds. Recently, it has been shown that the RIS model can be coupled with the Site Binding (SB) model, for which the ionizable sites can adopt two states: protonated or deprotonated. The resulting combined scheme, the SBRIS model, allows for analyzing ionization and conformational equilibria on the same foot. In the present work, this approach is extended to include pH-dependent electrostatic Long-Range (LR) interactions, ubiquitous in weak polyelectrolytes at moderate and low ionic strengths. With this aim, the original LR interactions are taken into account by defining effective Short-Range (SR) and pH-dependent parameters, such as effective microscopic protonation constants and rotational bond energies. The new parameters are systematically calculated using variational methods. The machinery of statistical mechanics for SR interactions, including the powerful and fast transfer matrix methods, can then be applied. The resulting technique, which we will refer to as the Local Effective Interaction Parameters (LEIP) method, is illustrated with a minimal model of a flexible linear polyelectrolyte containing only one type of rotating bond. LEIP reproduces very well the pH dependence of the degree of protonation and bond probabilities obtained by semi-grand canonical Monte Carlo simulations, where LR interactions are explicitly taken into account. The reduction in the computational time in several orders of magnitude suggests that the LEIP technique could be useful in a range of areas involving linear weak polyelectrolytes, allowing direct fitting of the relevant physical parameters to the experimental quantities.

Conditional equilibrium constants in multicomponent heterogeneous adsorption: the conditional affinity spectrum.

The concept of conditional stability constant is extended to the competitive binding of small molecules to heterogeneous surfaces or macromolecules via the introduction of the conditional affinity spectrum (CAS). The CAS describes the distribution of effective binding energies experienced by one complexing agent at a fixed concentration of the rest. We show that, when the multicomponent system can be described in terms of an underlying affinity spectrum [integral equation (IE) approach], the system can always be characterized by means of a CAS. The thermodynamic properties of the CAS and its dependence on the concentration of the rest of components are discussed. In the context of metal/proton competition, analytical expressions for the mean (conditional average affinity) and the variance (conditional heterogeneity) of the CAS as functions of pH are reported and their physical interpretation discussed. Furthermore, we show that the dependence of the CAS variance on pH allows for the analytical determination of the correlation coefficient between the binding energies of the metal and the proton. Nonideal competitive adsorption isotherm and Frumkin isotherms are used to illustrate the results of this work. Finally, the possibility of using CAS when the IE approach does not apply (for instance, when multidentate binding is present) is explored.