Potential of mean force and transient states in polyelectrolyte pair complexation

Polyelectrolyte complexation of two oppositely charged symmetric polymers: A minimal theory

Interplay of Coulomb interaction energy, free ion entropy, and conformational elasticity is a fascinating aspect in polyelectrolytes (PEs). We develop a theory for complexation of two oppositely charged PEs, a process known to be the precursor to the formation of complex coacervates in PE solutions, to explore the underlying thermodynamics of complex formation, at low salts. The theory considers general degrees of solvent polarity and dielectricity within an implicit solvent model, incorporating a varying Coulomb strength. Explicit calculation of the free energy of complexation and its components indicates that the entropy of free counterions and salt ions and the Coulomb enthalpy of bound ion-pairs dictate the equilibrium of PE complexation. This helps decouple the self-consistent dependency of charge and size of the uncomplexed parts of the polyions, derive an analytical expression for charge, and evaluate the free energy components as functions of chain overlap. Complexation is observed to be driven by enthalpy gain at low Coulomb strengths, driven by entropy gain of released counterions but opposed by enthalpy loss due to reduction of ion-pairs at moderate Coulomb strengths, and progressively less favorable due to enthalpy loss at even higher Coulomb strengths. The total free energy of the system is found to decrease linearly with an overlap of chains. Thermodynamic predictions from our model are in good quantitative agreement with simulations in literature.

Affinity of disordered protein complexes is modulated by entropy-energy reinforcement

Intrinsically disordered proteins (IDPs), which are very common and essential to many biological activities, sometimes function via interaction with another IDP and form a fuzzy complex, which can be highly stable. It is unclear what the biophysical forces are that govern their thermodynamics and specificity, which are essential for de novo fuzzy complex design. Here, we explored the fuzzy complex formed between ProTα and H1, which are oppositely charged IDPs, by swapping the charges between them, generating variants that have either greater polyampholytic or polyelectrolytic nature as well as different charge patterns. Charge swapping and shuffling dramatically change the affinity of the fuzzy complex, which is contributed to by both enthalpy and entropy, where the latter is dominated by counterion release.

The association between two intrinsically disordered proteins (IDPs) may produce a fuzzy complex characterized by a high binding affinity, similar to that found in the ultrastable complexes formed between two well-structured proteins. Here, using coarse-grained simulations, we quantified the biophysical forces driving the formation of such fuzzy complexes. We found that the high-affinity complex formed between the highly and oppositely charged H1 and ProTα proteins is sensitive to electrostatic interactions. We investigated 52 variants of the complex by swapping charges between the two oppositely charged proteins to produce sequences whose negatively or positively charged residue content was more homogeneous or heterogenous (i.e., polyelectrolytic or polyampholytic, having higher or lower absolute net charges, respectively) than the wild type. We also changed the distributions of oppositely charged residues within each participating sequence to produce variants in which the charges were segregated or well mixed. Both types of changes significantly affect binding affinity in fuzzy complexes, which is governed by both enthalpy and entropy. The formation of H1–ProTa is supported by an increase in configurational entropy and by entropy due to counterion release. The latter can be twice as large as the former, illustrating the dominance of counterion entropy in modulating the binding thermodynamics. Complexes formed between proteins with greater absolute net charges are more stable, both enthalpically and entropically, indicating that enthalpy and entropy have a mutually reinforcing effect. The sensitivity of the thermodynamics of the complex to net charge and the charge pattern within each of the binding constituents may provide a means to achieve binding specificity between IDPs.

Coarse-Grained Modeling of Ion-Containing Polymers

Ion-containing polymers have continued to be an important research focus for several decades due to their use as an electrolyte in energy storage and conversion devices. Elucidation of connections between the mesoscopic structure and multiscale dynamics of the ions and solvent remains incompletely understood. Coarse-grained modeling provides an efficient approach for exploring the structural and dynamical properties of these soft materials. The unique physicochemical properties of such polymers are of broad interest. In this review, we summarize the current development and understanding of the structure-property relationship of ion-containing polymers and provide insights into the design of such materials determined from coarse-grained modeling and simulations accompanying significant advances in experimental strategies. We specifically concentrate on three types of ion-containing polymers: proton exchange membranes (PEMs), anion exchange membranes (AEMs), and polymerized ionic liquids (polyILs). We posit that insight into the similarities and differences in these materials will lead to guidance in the rational design of high-performance novel materials with improved properties for various power source technologies.

Probing the protein corona around charged macromolecules: interpretation of isothermal titration calorimetry by binding models and computer simulations

Isothermal titration calorimetry (ITC) is a widely used tool to experimentally probe the heat signal of the formation of the protein corona around macromolecules or nanoparticles. If an appropriate binding model is applied to the ITC data, the heat of binding and the binding stoichiometry as well as the binding affinity per protein can be quantified and interpreted. However, the binding of the protein to the macromolecule is governed by complex microscopic interactions. In particular, due to the steric and electrostatic protein–protein interactions within the corona as well as cooperative, charge renormalization effects of the total complex, the application of standard (e.g., Langmuir) binding models is questionable and the development of more appropriate binding models is very challenging. Here, we discuss recent developments in the interpretation of the Langmuir model applied to ITC data of protein corona formation, exemplified for the well-defined case of lysozyme coating highly charged dendritic polyglycerol sulfate (dPGS), and demonstrate that meaningful data can be extracted from the fits if properly analyzed. As we show, this is particular useful for the interpretation of ITC data by molecular computer simulations where binding affinities can be calculated but it is often not clear how to consistently compare them with the ITC data. Moreover, we discuss the connection of Langmuir models to continuum binding models (where no discrete binding sites have to be assumed) and their possible extensions toward the inclusion of leading order cooperative electrostatic effects.

Recent progress in the science of complex coacervation

Complex coacervation is an associative, liquid-liquid phase separation that can occur in solutions of oppositely-charged macromolecular species, such as proteins, polymers, and colloids. This process results in a coacervate phase, which is a dense mix of the oppositely-charged components, and a supernatant phase, which is primarily devoid of these same species. First observed almost a century ago, coacervates have since found relevance in a wide range of applications; they are used in personal care and food products, cutting edge biotechnology, and as a motif for materials design and self-assembly. There has recently been a renaissance in our understanding of this important class of material phenomena, bringing the science of coacervation to the forefront of polymer and colloid science, biophysics, and industrial materials design. In this review, we describe the emergence of a number of these new research directions, specifically in the context of polymer-polymer complex coacervates, which are inspired by a number of key physical and chemical insights and driven by a diverse range of experimental, theoretical, and computational approaches.

Role of Associative Charging in the Entropy-Energy Balance of Polyelectrolyte Complexes

Polyelectrolytes may be classified into two primary categories (strong and weak) depending on how their charge state responds to the local environment. Both of these find use in many applications, including drug delivery, gene therapy, layer-by-layer films, and fabrication of ion filtration membranes. The mechanism of polyelectrolyte complexation is, however, still not completely understood, though experimental investigations suggest that entropy gain due to release of counterions is the key driving force for strong polyelectrolyte complexation. Here we perform a comprehensive thermodynamic investigation through coarse-grained molecular simulations permitting us to calculate the free energy of complex formation. Importantly, our expanded-ensemble methods permit the explicit separation of energetic and entropic contributions to the free energy. Our investigations indicate that entropic contributions indeed dominate the free energy of complex formation for strong polyelectrolytes, but are less important than energetic contributions when weak electrostatic coupling or weak polyelectrolytes are present. Our results provide a new view of the free energy of polyelectrolyte complex formation driven by polymer association, which should also arise in systems with large charge spacings or bulky counterions, both of which act to weaken ion-polymer binding.

The Mechanism for siRNA Transmembrane Assisted by PMAL

The capacity of silencing genes makes small interfering RNA (siRNA) appealing for curing fatal diseases. However, the naked siRNA is vulnerable to and degraded by endogenous enzymes and is too large and too negatively charged to cross cellular membranes. An effective siRNA carrier, PMAL (poly(maleic anhydride-alt-1-decene) substituted with 3-(dimethylamino) propylamine), has been demonstrated to be able to assist siRNA transmembrane by both experiments and molecular simulation. In the present work, the mechanism of siRNA transmembrane assisted by PMAL was studied using steered molecular dynamics simulations based on the martini coarse-grained model. Here two pulling rates, i.e., 10⁻⁶ and 10⁻⁵ nm·ps⁻¹, were chosen to imitate the passive and active transport of siRNA, respectively. Potential of mean force (PMF) and interactions among siRNA, PMAL, and lipid bilayer membrane were calculated to describe the energy change during siRNA transmembrane processes at various conditions. It is shown that PMAL-assisted siRNA delivery is in the mode of passive transport. The PMAL can help siRNA insert into lipid bilayer membrane by lowering the energy barrier caused by siRNA and lipid bilayer membrane. PMAL prefers to remain in the lipid bilayer membrane and release siRNA. The above simulations establish a molecular insight of the interaction between siRNA and PMAL and are helpful for the design and applications of new carriers for siRNA delivery.

Counterion-Release Entropy Governs the Inhibition of Serum Proteins by Polyelectrolyte Drugs

Dendritic polyelectrolytes constitute high potential drugs and carrier systems for biomedical purposes. Still, their biomolecular interaction modes, in particular those determining the binding affinity to proteins, have not been rationalized. We study the interaction of the drug candidate dendritic polyglycerol sulfate (dPGS) with serum proteins using isothermal titration calorimetry (ITC) interpreted and complemented with molecular computer simulations. Lysozyme is first studied as a well-defined model protein to verify theoretical concepts, which are then applied to the important cell adhesion protein family of selectins. We demonstrate that the driving force of the strong complexation, leading to a distinct protein corona, originates mainly from the release of only a few condensed counterions from the dPGS upon binding. The binding constant shows a surprisingly weak dependence on dPGS size (and bare charge) which can be understood by colloidal charge-renormalization effects and by the fact that the magnitude of the dominating counterion-release mechanism almost exclusively depends on the interfacial charge structure of the protein-specific binding patch. Our findings explain the high selectivity of P- and L-selectins over E-selectin for dPGS to act as a highly anti-inflammatory drug. The entire analysis demonstrates that the interaction of proteins with charged polymeric drugs can be predicted by simulations with unprecedented accuracy. Thus, our results open new perspectives for the rational design of charged polymeric drugs and carrier systems.

Understanding the differences of the ligand binding/unbinding pathways between phosphorylated and non-phosphorylated ARH1 using molecular dynamics simulations

ADP-ribosylhydrolases (ARH1, ARH2 and ARH3) are a family of enzymes to catalyze ADP-ribosylation, a reversible and covalent post-translational modification (PTM). There are four phosphorylated sites (Tyr-4, Tyr-19, Tyr-20, and Tyr-205) in ARH1. To explore the structural changes and functional impact induced by phosphorylation, molecular dynamics (MD) simulations and steered molecular dynamics (SMD) simulations were performed for the phosphorylated and non-phosphorylated ARH1 with the ligands. MD simulations results indicate that: (1) Glu-25 is more frequently in the α helix group in the phosphorylated state with the adenosine-5-diphosphate-ribosylarginine (ADP-RA) complex (51.56%) than that of the non-phosphorylated state(2.12%); (2) Ser-124 and Ser-264 become less flexible in the phosphorylated state with ADP-RA complex, which helps two residues form hydrogen bonds with ADP-RA; and (3) Tyr-211 is also less flexible in the phosphorylated state with ADP-RA complex, which helps stabilize the cation-π interaction of Y211-R119. All these changes facilitate ADP-RA to bind ARH1. In addition, according to the crystal structure of adenosine-5-diphosphate-ribose (ADP-ribose) in complex with non-phosphorylated and phosphorylated ARH1, the possible unbinding pathways of ADP-ribose from non-phosphorylated and phosphorylated ARH1 were explored respectively using SMD simulations. Our results show that phosphorylated ARH1 has more ordered structures than the non-phosphorylated type.