pyMBE: The Python-based molecule builder for ESPResSo.

We present the Python-based Molecule Builder for ESPResSo (pyMBE), an open source software application to design custom coarse-grained (CG) models, as well as pre-defined models of polyelectrolytes, peptides, and globular proteins in the Extensible Simulation Package for Research on Soft Matter (ESPResSo). The Python interface of ESPResSo offers a flexible framework, capable of building custom CG models from scratch. As a downside, building CG models from scratch is prone to mistakes, especially for newcomers in the field of CG modeling, or for molecules with complex architectures. The pyMBE module builds CG models in ESPResSo using a hierarchical bottom-up approach, providing a robust tool to automate the setup of CG models and helping new users prevent common mistakes. ESPResSo features the constant pH (cpH) and grand-reaction (G-RxMC) methods, which have been designed to study chemical reaction equilibria in macromolecular systems with many reactive species. However, setting up these methods for systems, which contain several types of reactive groups, is an error-prone task, especially for beginners. The pyMBE module enables the automatic setup of cpH and G-RxMC simulations in ESPResSo, lowering the barrier for newcomers and opening the door to investigate complex systems not studied with these methods yet. To demonstrate some of the applications of pyMBE, we showcase several case studies where we successfully reproduce previously published simulations of charge-regulating peptides and globular proteins in bulk solution and weak polyelectrolytes in dialysis. The pyMBE module is publicly available as a GitHub repository (https://github.com/pyMBE-dev/pyMBE), which includes its source code and various sample and test scripts, including the ones that we used to generate the data presented in this article.

Charge Regulation Triggers Condensation of Short Oligopeptides to Polyelectrolytes.

Electrostatic interactions between charged macromolecules are ubiquitous in biological systems, and they are important also in materials design. Attraction between oppositely charged molecules is often interpreted as if the molecules had a fixed charge, which is not affected by their interaction. Less commonly, charge regulation is invoked to interpret such interactions, i.e., a change of the charge state in response to a change of the local environment. Although some theoretical and simulation studies suggest that charge regulation plays an important role in intermolecular interactions, experimental evidence supporting such a view is very scarce. In the current study, we used a model system, composed of a long polyanion interacting with cationic oligolysines, containing up to 8 lysine residues. We showed using both simulations and experiments that while these lysines are only weakly charged in the absence of the polyanion, they charge up and condense on the polycations if the pH is close to the pKa of the lysine side chains. We show that the lysines coexist in two distinct populations within the same solution: (1) practically nonionized and free in solution; (2) highly ionized and condensed on the polyanion. Using this model system, we demonstrate under what conditions charge regulation plays a significant role in the interactions of oppositely charged macromolecules and generalize our findings beyond the specific system used here.

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.

Reply to the 'Comment on "Simulations of ionization equilibria in weak polyelectrolyte solutions and gels"' by J. Landsgesell, L. Nová, O. Rud, F. Uhlík, D. Sean, P. Hebbeker, C. Holm and P. Kosovan, Soft Matter, 2019, 15, 1155-1185.

Levin and Bakhshandeh suggested in their comment that (1), we stated in our recent review that pH-pKA is a universal parameter for titrating systems, that (2), we omitted to mention in our review the broken symmetry of the constant pH algorithm, and that (3), a constant pH simulation must include a grand-canonical exchange of ions with the reservoir. As a reply to (1), we point out that Levin and Bakhshandeh misquoted and hence invalidated our original statement. We therefore explain in detail under which circumstances pH-pKA can be a universal parameter, and also demonstrate why their numerical example is not in contradiction to our statement. Moreover, the fact that pH-pKA is not a universal parameter for titrating systems is well known in the pertinent literature. Regarding (2), we admit that the symmetry-breaking feature of the constant pH algorithm has escaped our attention at the time of writing the review. We added some clarifying remarks to this behavior. Concerning (3), we point out that the grand-canonical coupling and the resultant Donnan potential are not features of single-phase systems, but are essential for two-phase systems, as was shown in a recent paper by some of us, see J. Landsgesell et al., Macromolecules, 2020, 53, 3007-3020.

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.

Anionically Functionalized Glycogen Encapsulates Melittin by Multivalent Interaction.

We developed acid-functionalized glycogen conjugates as supramolecular carriers for efficient encapsulation and inhibition of a model cationic peptide melittin─the main component of honeybee venom. For this purpose, we synthesized and characterized a set of glycogens, functionalized to various degrees by several different acid groups. These conjugates encapsulate melittin up to a certain threshold amount, beyond which they precipitate. Computer simulations showed that sufficiently functionalized conjugates electrostatically attract melittin, resulting in its efficient encapsulation in a broad pH range around the physiological pH. Hemolytic assays confirmed in vitro that the effective inhibition of melittin's hemolytic activity occurs for highly functionalized samples, whereas no inhibition is observed when using low-functionalized conjugates. It can be concluded that functional glycogens are promising carriers for cationic molecular cargos or antidotes against animal venoms under conditions, in which suitable properties such as biodegradability and biocompatibility are crucial.

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.

Macromolecular diffusion in crowded media beyond the hard-sphere model.

The effect of macromolecular crowding on diffusion beyond the hard-core sphere model is studied. A new coarse-grained model is presented, the Chain Entanglement Softened Potential (CESP) model, which takes into account the macromolecular flexibility and chain entanglement. The CESP model uses a shoulder-shaped interaction potential that is implemented in the Brownian Dynamics (BD) computations. The interaction potential contains only one parameter associated with the chain entanglement energetic cost (Ur). The hydrodynamic interactions are included in the BD computations via Tokuyama mean-field equations. The model is used to analyze the diffusion of a streptavidin protein among different sized dextran obstacles. For this system, Ur is obtained by fitting the streptavidin experimental long-time diffusion coefficient Dlongversus the macromolecular concentration for D50 (indicating their molecular weight in kg mol-1) dextran obstacles. The obtained Dlong values show better quantitative agreement with experiments than those obtained with hard-core spheres. Moreover, once parametrized, the CESP model is also able to quantitatively predict Dlong and the anomalous exponent (α) for streptavidin diffusion among D10, D400 and D700 dextran obstacles. Dlong, the short-time diffusion coefficient (Dshort) and α are obtained from the BD simulations by using a new empirical expression, able to describe the full temporal evolution of the diffusion coefficient.

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.