Conformation and Ionization Behavior of Charge-Regulating Polyelectrolyte Brushes in a Poor Solvent

Unexpected Two-Stage Swelling of Weak Polyelectrolyte Brushes with Divalent Counterions

We use particle-based, coarse-grained simulations to study the influence of divalent counterions on a weak polyelectrolyte brush. Our simulations show a profound influence of even small concentrations of divalent salt on the titration behavior of the brush, which is shown to be a combined effect of electrostatic interactions and the Donnan effect. Furthermore, we examine the partitioning of mono- and divalent counterions into the brush. We demonstrate the preferred uptake of divalent ions by the brush, which is further enhanced by electrostatic correlation effects. Finally, our simulations reveal a hitherto unobserved two-stage swelling of the brush as a function of the pH in the presence of divalent salt. This phenomenon arises as a consequence of charge regulation and ion partitioning.

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.

Impact of Charge Regulation on Self-Assembly of Zwitterionic Nanoparticles

Zwitterionic modification of colloids with weak acids and bases represents a promising strategy in creating functional materials with tunable properties and modeling the self-organization of charged proteins. However, accurate incorporation of the dynamic dissociation or association of ionization groups known as charge regulation (CR) is often intractable in theoretical and computational investigations since charge redistribution and configuration need to be evolved self-consistently. Using hybrid Monte Carlo and molecular dynamics simulations, we demonstrate that a dilute suspension of overall charge-neutral zwitterionic Janus nanoparticles shows a conformational transition from an open assembly of string or bundle to compact cluster along with the variation in pH. The behavior under CR is qualitatively different from the commonly employed constant charge condition where the transition is absent. The CR-induced clustering is due to the inhomogeneous and fluctuating charges localized near the equatorial boundary of the Janus particle. These features are enhanced particularly at low salt concentration and high electrostatic coupling strength. Our results indicate the critical role of charge regulation in the spatial self-organization of zwitterionic nanoparticles.