Computer study of the solubilization of polymer chains in polyelectrolyte complex cores of polymeric nanoparticles in aqueous media

Changes in Ion Concentrations upon the Binding of Short Polyelectrolytes on Phospholipid Bilayers: Computer Study Addressing Interesting Physiological Consequences

This computer study was inspired by the experimental observation of Y. Qian et al. published in ACS Applied Materials and Interfaces, 2018 that the short positively charged β-peptide chains and their oligomeric analogues efficiently suppress severe medical problems caused by antimicrobial drug-resistant bacteria despite them not penetrating the bacterial membrane. Our coarse-grained molecular dynamics (dissipative particle dynamics) simulations confirm the tentative explanation of the authors of the experimental study that the potent antimicrobial activity is a result of the entropically driven release of divalent ions (mainly magnesium ions essential for the proper biological function of bacteria) into bulk solution upon the electrostatic binding of β-peptides to the bacterial membrane. The study shows that in solutions containing cations Na⁺, Ca²⁺ and Mg²⁺, and anions Cl⁻, the divalent cations preferentially concentrate close to the membrane and neutralize the negative charge. Upon the addition of positively charged oligomer chains (models of β-peptides and their analogues), the oligomers electrostatically bind to the membrane replacing divalent ions, which are released into bulk solvent. Our simulations indicate that the entropy of small ions (which controls the behavior of synthetic polyelectrolyte solutions) plays an important role in this and also in other similar biologically important systems.

DPD Modelling of the Self- and Co-Assembly of Polymers and Polyelectrolytes in Aqueous Media: Impact on Polymer Science

This review article is addressed to a broad community of polymer scientists. We outline and analyse the fundamentals of the dissipative particle dynamics (DPD) simulation method from the point of view of polymer physics and review the articles on polymer systems published in approximately the last two decades, focusing on their impact on macromolecular science. Special attention is devoted to polymer and polyelectrolyte self- and co-assembly and self-organisation and to the problems connected with the implementation of explicit electrostatics in DPD numerical machinery. Critical analysis of the results of a number of successful DPD studies of complex polymer systems published recently documents the importance and suitability of this coarse-grained method for studying polymer systems.

Temperature-responsive morphology formation of a PS-b-PI copolymer: a dissipative particle dynamics simulation study

Self-assembly responsiveness to stimuli of polystyrene-block-polyisoprene (PS-b-PI) diblock copolymer materials is explored by means of classical molecular dynamics (MD) and dissipative particle dynamics (DPD) simulations. A concerted relationship between the parameters achieved from atomistic and DPD simulations is obtained for this molecular recognition as clearly pronounced in a phase transition. Effects of temperature, model size and composition on the morphological formation were systematically investigated for the diblock copolymeric system. Structural changes resulting in the evolution of rheology as well as an equilibrium ordered structure were analyzed in terms of order parameters and radial distribution functions. From our models, various morphologies were observed including discrete clusters (sphere-liked morphology), connected clusters (gyroid-liked morphology), hexagonally packed cylinders (HEX), connected cylinders, irregular cylinders, perfect lamellae, perforated lamellae and defected lamellae. Based on this finding, a bottom-up multi-scale simulation of the PS-b-PI diblock copolymer provides a link between equilibrium copolymeric morphologies and the crucial parameters.

Solubilization of Charged Porphyrins in Interpolyelectrolyte Complexes: A Computer Study

Using coarse-grained dissipative particle dynamics (DPD) with explicit electrostatics, we performed (i) an extensive series of simulations of the electrostatic co-assembly of asymmetric oppositely charged copolymers composed of one (either positively or negatively charged) polyelectrolyte (PE) block A and one water-soluble block B and (ii) studied the solubilization of positively charged porphyrin derivatives (P+) in the interpolyelectrolyte complex (IPEC) cores of co-assembled nanoparticles. We studied the stoichiometric mixtures of 137 A10+B25 and 137 A10-B25 chains with moderately hydrophobic A blocks (DPD interaction parameter aAS=35) and hydrophilic B blocks (aBS=25) with 10 to 120 P+ added (aPS=39). The P+ interactions with other components were set to match literature information on their limited solubility and aggregation behavior. The study shows that the moderately soluble P+ molecules easily solubilize in IPEC cores, where they partly replace PE+ and electrostatically crosslink PE- blocks. As the large P+ rings are apt to aggregate, P+ molecules aggregate in IPEC cores. The aggregation, which starts at very low loadings, is promoted by increasing the number of P+ in the mixture. The positively charged copolymers repelled from the central part of IPEC core partially concentrate at the core-shell interface and partially escape into bulk solvent depending on the amount of P+ in the mixture and on their association number, AS. If AS is lower than the ensemble average ⟨AS⟩n, the copolymer chains released from IPEC preferentially concentrate at the core-shell interface, thus increasing AS, which approaches ⟨AS⟩n. If AS>⟨AS⟩n, they escape into the bulk solvent.

Langevin Dynamics Simulations of the Exchange of Complex Coacervate Core Micelles: The Role of Nonelectrostatic Attraction and Polyelectrolyte Length

Complex coacervate core micelles (C3Ms) are promising encapsulators for a wide variety of (bio)molecules. To protect and stabilize their cargo, it is essential to control their exchange dynamics. Yet, to date, little is known about the kinetic stability of C3Ms and the dynamic equilibrium of molecular building blocks with micellar species. Here we study the C3M exchange during the initial micellization by using Langevin dynamics simulations. In this way, we show that charge neutral heterocomplexes consisting of multiple building blocks are the primary mediator for exchange. In addition, we show that the kinetic stability of the C3Ms can be tuned not only by the electrostatic interaction but also by the nonelectrostatic attraction between the polyelectrolytes, the polyelectrolyte length ratio, and the overall polyelectrolyte length. These insights offer new rational design guides to aid the development of new C3M encapsulation strategies.