Transfer Matrix Model of pH Effects in Polymeric Complex Coacervation

Interactions between Hyaluronic Acid and Chitosan by Isothermal Titration Calorimetry: The Effect of Ionic Strength, pH, and Polymer Molecular Weight

Hyaluronic acid (HA)/chitosan (CHI) complex coacervates have recently gained interest due to the pH-dependent ionization and semiflexibility of the polymers as well as their applicability in tissue engineering. Here, we apply isothermal titration calorimetry (ITC) to understand the apparent thermodynamics of coacervation for HA/CHI as a function of the pH, ionic strength, and chain length. We couple these ITC experiments with the knowledge of the charge states of HA and CHI from potentiometric titration to understand the mechanistic aspects of complex formation. Our data demonstrate that the driving force for the complex coacervation of HA and CHI is entropic in nature and this driving force decreased with increasing ionic strength. We also observed a decrease in the stoichiometry for ion-pairing with increasing ionic strength, which we suggest is a consequence of the changing degree of ionization for HA at higher ionic strengths. An increase in the strength of interactions with pH was hypothesized to also be a result of changes in the degree of ionization of HA, though stronger interactions were observed at the lowest pH tested, likely due to contributions from hydrogen bonding between HA and CHI.

Theory and quantitative assessment of pH-responsive polyzwitterion-polyelectrolyte complexation

We introduce a theoretical framework to describe the pH-sensitive phase behavior of polyzwitterion-polyelectrolyte complex coacervates that reasonably captures the phenomenon from recent experimental observations. The polyzwitterion is described by a combinatorial sequence of the four states in which each zwitterionic monomer can occupy: dipolar, quasi-cationic, quasi-anionic, and fully neutralized. We explore the effects of various modifiable chemical and physical properties of the polymers-such as, pKa of the pH-active charged group on the zwitterion, equilibrium constant of salt condensation on the permanently charged group on the zwitterion, degrees of polymerization, hydrophobicity (via the Flory-Huggins interaction parameter), and dipole lengths-on the window of complexation across many stoichiometric mixing ratios of polyzwitterion and polyelectrolyte. The properties that determine the net charge of the polyzwitterion have the strongest effect on the pH range in which polyzwitterion-polyelectrolyte complexation occurs. We finish with general guidance for those interested in molecular design of polyzwitterion-polyelectrolyte complex coacervates and opportunities for future investigation.

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.

Effect of Polymer Gel Elasticity on Complex Coacervate Phase Behavior

Gels are key materials in biological systems such as tissues and may control biocondensate formation and structure. To further understand the effects of elastic environments on biomacromolecular assembly, we have investigated the phase behavior and radii of complex coacervate droplets in polyacrylamide (PAM) networks as a function of gel modulus. Poly-l-lysine (PLL) and sodium hyaluronate (HA) complex coacervate phases were prepared in PAM gels with moduli varying from 0.035 to 15.0 kPa. The size of the complex coacervate droplets is reported from bright-field microscopy and confocal fluorescence microscopy. Overall, the complex coacervate droplet volume decreases inversely with the modulus. Fluorescence microscopy is also used to determine the phase behavior and concentration of fluorescently tagged HA in the complex coacervate phases as a function of ionic strength (100-270 mM). We find that the critical ionic strength and complex coacervate stability are nonmonotonic as a function of the network modulus and that the local gel concentration can be used to control phase behavior and complex coacervate droplet size scale. By understanding how elastic environments influence simple electrostatic assembly, we can further understand how biomacromolecules exist in complex, crowded, and elastic cellular environments.

Sequestration of Small Ions and Weak Acids and Bases by a Polyelectrolyte Complex Studied by Simulation and Experiment

Mixing of oppositely charged polyelectrolytes can result in phase separation into a polymer-poor supernatant and a polymer-rich polyelectrolyte complex (PEC). We present a new coarse-grained model for the Grand-reaction method that enables us to determine the composition of the coexisting phases in a broad range of pH and salt concentrations. We validate the model by comparing it to recent simulations and experimental studies, as well as our own experiments on poly(acrylic acid)/poly(allylamine hydrochloride) complexes. The simulations using our model predict that monovalent ions partition approximately equally between both phases, whereas divalent ones accumulate in the PEC phase. On a semiquantitative level, these results agree with our own experiments, as well as with other experiments and simulations in the literature. In the sequel, we use the model to study the partitioning of a weak diprotic acid at various pH values of the supernatant. Our results show that the ionization of the acid is enhanced in the PEC phase, resulting in its preferential accumulation in this phase, which monotonically increases with the pH. Currently, this effect is still waiting to be confirmed experimentally. We explore how the model parameters (particle size, charge density, permittivity, and solvent quality) affect the measured partition coefficients, showing that fine-tuning of these parameters can make the agreement with the experiments almost quantitative. Nevertheless, our results show that charge regulation in multivalent solutes can potentially be exploited in engineering the partitioning of charged molecules in PEC-based systems at various pH values.

Design Rules for the Sequestration of Viruses into Polypeptide Complex Coacervates

Encapsulation is a strategy that has been used to facilitate the delivery and increase the stability of proteins and viruses. Here, we investigate the encapsulation of viruses via complex coacervation, which is a liquid-liquid phase separation resulting from the complexation of oppositely charged polymers. In particular, we utilized polypeptide-based coacervates and explored the effects of peptide chemistry, chain length, charge patterning, and hydrophobicity to better understand the effects of the coacervating polypeptides on virus incorporation. Our study utilized two nonenveloped viruses, porcine parvovirus (PPV) and human rhinovirus (HRV). PPV has a higher charge density than HRV, and they both appear to be relatively hydrophobic. These viruses were compared to characterize how the charge, hydrophobicity, and patterning of chemistry on the surface of the virus capsid affects encapsulation. Consistent with the electrostatic nature of complex coacervation, our results suggest that electrostatic effects associated with the net charge of both the virus and polypeptide dominated the potential for incorporating the virus into a coacervate, with clustering of charges also playing a significant role. Additionally, the hydrophobicity of a virus appears to determine the degree to which increasing the hydrophobicity of the coacervating peptides can enhance virus uptake. Nonintuitive trends in uptake were observed with regard to both charge patterning and polypeptide chain length, with these parameters having a significant effect on the range of coacervate compositions over which virus incorporation was observed. These results provide insights into biophysical mechanisms, where sequence effects can control the uptake of proteins or viruses into biological condensates and provide insights for use in formulation strategies.

Electroresponse of weak polyelectrolyte brushes

End-tethered polyelectrolytes are widely used to modify substrate properties, particularly for lubrication or wetting. External stimuli, such as pH, salt concentration, or an electric field, can induce profound structural responses in weak polyelectrolyte brushes, which can be utilized to further tune substrate properties. We study the structure and electroresponsiveness of weak polyacid brushes using an inhomogeneous theory that incorporates both electrostatic and chain connectivity correlations at the Debye-Hückel level. Our calculation shows that a weak polyacid brush swells under the application of a negative applied potential, in agreement with recent experimental observation. We rationalize this behavior using a scaling argument that accounts for the effect of the surface charge. We also show that the swelling behavior has a direct influence on the differential capacitance, which can be modulated by the solvent quality, pH, and salt concentration.

Phase Transitions of Associative Biomacromolecules

Multivalent proteins and nucleic acids, collectively referred to as multivalent associative biomacromolecules, provide the driving forces for the formation and compositional regulation of biomolecular condensates. Here, we review the key concepts of phase transitions of aqueous solutions of associative biomacromolecules, specifically proteins that include folded domains and intrinsically disordered regions. The phase transitions of these systems come under the rubric of coupled associative and segregative transitions. The concepts underlying these processes are presented, and their relevance to biomolecular condensates is discussed.

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.

Driving force and pathway in polyelectrolyte complex coacervation

There is notable discrepancy between experiments and coarse-grained model studies regarding the thermodynamic driving force in polyelectrolyte complex coacervation: experiments find the free energy change to be dominated by entropy, while simulations using coarse-grained models with implicit solvent usually report a large, even dominant energetic contribution in systems with weak to intermediate electrostatic strength. Here, using coarse-grained, implicit-solvent molecular dynamics simulation combined with thermodynamic analysis, we study the potential of mean force (PMF) in the two key stages on the coacervation pathway for symmetric polyelectrolyte mixtures: polycation-polyanion complexation and polyion pair-pair condensation. We show that the temperature dependence in the dielectric constant of water gives rise to a substantial entropic contribution in the electrostatic interaction. By accounting for this electrostatic entropy, which is due to solvent reorganization, we find that under common conditions (monovalent ions, room temperature) for aqueous systems, both stages are strongly entropy-driven with negligible or even unfavorable energetic contributions, consistent with experimental results. Furthermore, for weak to intermediate electrostatic strengths, this electrostatic entropy, rather than the counterion-release entropy, is the primary entropy contribution. From the calculated PMF, we find that the supernatant phase consists predominantly of polyion pairs with vanishingly small concentration of bare polyelectrolytes, and we provide an estimate of the spinodal of the supernatant phase. Finally, we show that prior to contact, two neutral polyion pairs weakly attract each other by mutually induced polarization, providing the initial driving force for the fusion of the pairs.

Field theory description of ion association in re-entrant phase separation of polyampholytes

Phase separation of several different overall neutral polyampholyte species (with zero net charge) is studied in solution with two oppositely charged ion species that can form ion-pairs through an association reaction. A field theory description of the system, that treats polyampholyte charge sequence dependent electrostatic interactions as well as excluded volume effects, is hereby given. Interestingly, analysis of the model using random phase approximation and field theoretic simulation consistently show evidence of a re-entrant polyampholyte phase separation at high ion concentrations when there is an overall decrease of volume upon ion-association. As an illustration of the ramifications of our theoretical framework, several polyampholyte concentration vs ion concentration phase diagrams under constant temperature conditions are presented to elucidate the dependence of phase separation behavior on polyampholyte sequence charge pattern as well as ion-pair dissociation constant, volumetric effects on ion association, solvent quality, and temperature.

Quantifying charge state heterogeneity for proteins with multiple ionizable residues

Ionizable residues can release and take up protons and this has an influence on protein structure and function. The extent of protonation is linked to the overall pH of the solution and the local environments of ionizable residues. Binding or unbinding of a single proton generates a distinct charge microstate defined by a specific pattern of charges. Accordingly, the overall partition function is a sum over all charge microstates and Boltzmann weights of all conformations associated with each of the charge microstates. This ensemble-of-ensembles description recast as a q-canonical ensemble allows us to analyze and interpret potentiometric titrations that provide information regarding net charge as a function of pH. In the q-canonical ensemble, charge microstates are grouped into mesostates where each mesostate is a collection of microstates of the same net charge. Here, we show that leveraging the structure of the q-canonical ensemble allows us to decouple contributions of net proton binding and release from proton arrangement and conformational considerations. Through application of the q-canonical formalism to analyze potentiometric measurements of net charge in proteins with repetitive patterns of Lys and Glu residues, we determine the underlying mesostate pKa values and, more importantly, we estimate relative mesostate populations as a function of pH. This is a strength of using the q-canonical approach that cannot be replicated using purely site-specific analyses. Overall, our work shows how measurements of charge equilibria, decoupled from measurements of conformational equilibria, and analyzed using the framework of the q-canonical ensemble, provide protein-specific quantitative descriptions of pH-dependent populations of mesostates. This method is of direct relevance for measuring and understanding how different charge states contribute to conformational, binding, and phase equilibria of proteins.