A Modified Random Phase Approximation of Polyelectrolyte Solutions

Differential Effects of Sequence-Local versus Nonlocal Charge Patterns on Phase Separation and Conformational Dimensions of Polyampholytes as Model Intrinsically Disordered Proteins

Conformational properties of intrinsically disordered proteins (IDPs) are governed by a sequence-ensemble relationship. To differentiate the impact of sequence-local versus sequence-nonlocal features of an IDP's charge pattern on its conformational dimensions and its phase-separation propensity, the charge "blockiness" κ and the nonlocality-weighted sequence charge decoration (SCD) parameters are compared for their correlations with isolated-chain radii of gyration (Rgs) and upper critical solution temperatures (UCSTs) of polyampholytes modeled by random phase approximation, field-theoretic simulation, and coarse-grained molecular dynamics. SCD is superior to κ in predicting Rg because SCD accounts for effects of contact order, i.e., nonlocality, on dimensions of isolated chains. In contrast, κ and SCD are comparably good, though nonideal, predictors of UCST because frequencies of interchain contacts in the multiple-chain condensed phase are less sensitive to sequence positions than frequencies of intrachain contacts of an isolated chain, as reflected by κ correlating better with condensed-phase interaction energy than SCD.

Preserving condensate structure and composition by lowering sequence complexity

Biomolecular condensates play a vital role in organizing cellular chemistry. They selectively partition biomolecules, preventing unwanted cross talk and buffering against chemical noise. Intrinsically disordered proteins (IDPs) serve as primary components of these condensates due to their flexibility and ability to engage in multivalent interactions, leading to spontaneous aggregation. Theoretical advancements are critical at connecting IDP sequences with condensate emergent properties to establish the so-called molecular grammar. We proposed an extension to the stickers and spacers model, incorporating heterogeneous, nonspecific pairwise interactions between spacers alongside specific interactions among stickers. Our investigation revealed that although spacer interactions contribute to phase separation and co-condensation, their nonspecific nature leads to disorganized condensates. Specific sticker-sticker interactions drive the formation of condensates with well-defined networked structures and molecular composition. We discussed how evolutionary pressures might emerge to affect these interactions, leading to the prevalence of low-complexity domains in IDP sequences. These domains suppress spurious interactions and facilitate the formation of biologically meaningful condensates.

Theoretical treatment of complex coacervate core micelles: structure and pH-induced disassembly

Complex coacervate core micelles (C3Ms) are supramolecular soft nanostructures formed by the assembly of a block copolymer and an oppositely charged homopolymer. The coacervation of the charged segments in both macromolecules drives the formation of the core of the C3M, while the neutral block of the copolymer forms the corona. This work introduces a molecular theory (MOLT) that predicts the internal structure and stimuli-responsive properties of C3Ms and explicitly considers the chemical architecture of the polyelectrolytes, their acid-based equilibria and electrostatic and non-electrostatic interactions. In order to accurately predict complex coacervation, the correlations between charged species are incorporated into MOLT as ion-pairing processes, which are modeled using a coupled chemical equilibrium formalism. Very good agreement was observed between the experimental results in the literature and MOLT predictions for the scaling relationships that relate the dimensions of the micelle (aggregation number and sizes of the micelle and the core) to the lengths of the different blocks. MOLT was used to study the disassembly of the micelles when the solution pH is driven away from the value that guarantees the charge stoichiometry of the core. This study reveals that very sharp disassembly transitions can be obtained by tuning the length or architecture of the copolymer component, thereby suggesting potential routes to design C3Ms capable of releasing their components at very precise pH values.

Phase Behavior and Electrochemical Properties of Highly Asymmetric Redox Coacervates

This work reports the phase behavior and electrochemical properties of liquid coacervates made of ferricyanide and poly(ethylenimine). In contrast to the typical polyanion/polycation pairs used in liquid coacervates, the ferricyanide/poly(ethylenimine) system is highly asymmetric because poly(ethylenimine) has approximately 170 charges per molecule, while ferricyanide has only 3. Two types of phase diagrams were measured and fitted with a theoretical model. In the first type of diagram, the stability of the coacervate was studied in the plane given by the concentration of poly(ethylenimine) versus the concentration of ferricyanide for a fixed concentration of added monovalent salt (NaCl). The second type of diagram involved the plane given by the concentration of poly(ethylenimine) vs the concentration of the added monovalent salt for a fixed poly(ethyleneimine)/ferricyanide ratio. Interestingly, these phase diagrams displayed qualitative similarities to those of symmetric polyanion/polycation systems, suggesting that coacervates formed by a polyelectrolyte and a small multivalent ion can be treated as a specific case of polyelectrolyte coacervate. The characterization of the electrochemical properties of the coacervate revealed that the addition of monovalent salt greatly enhances charge transport, presumably by breaking ion pairs between ferricyanide and poly(ethylenimine). This finding highlights the significant influence of added salt on the transport properties of coacervates. This study provides the first comprehensive characterization of the phase behavior and transport properties of asymmetric coacervates and places these results within the broader context of the better-known symmetric polyelectrolyte coacervates.

Theory of self-coacervation in semi-dilute and concentrated zwitterionic polymer solutions

Based on the random phase approximation, we develop a molecular theory of self-coacervation in zwitterionic polymer solutions. We show that the interplay between the volume interactions of the monomeric units and electrostatic correlations of charged groups on a polymer backbone can result in liquid-liquid phase separation (self-coacervation). We analyse the behavior of the coacervate phase polymer concentration depending on the electrostatic interaction strength - the ratio of the Bjerrum length to the bond length of the chain. We establish that in a wide range of polymer concentration values - from a semi-dilute to a rather concentrated solution - the chain connectivity and excluded volume interaction of the monomeric units have an extremely weak effect on the contribution of the electrostatic interactions of the dipolar monomeric units to the total free energy. We show that for rather weak electrostatic interactions, the electrostatic correlations manifest themselves as Keesom interactions of point-like freely rotating dipoles (Keesom regime), while in the region of strong electrostatic interactions the electrostatic free energy is described by the Debye-Hückel limiting law (Debye regime). We show that for real zwitterionic coacervates the Keesom regime is realized only for sufficiently small polymer concentrations of the coacervate phase, while the Debye regime is approximately realized for rather dense coacervates. Using the mean-field variant of the density functional theory, we calculate the surface tension (surface free energy) of the "coacervate-solvent" interface as a function of the bulk polymer concentration. Obtained results can be used to estimate the parameters of the polymer chains needed for practical applications such as drug encapsulation and delivery, as well as the design of adhesive materials.

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.

Numerical Techniques for Applications of Analytical Theories to Sequence-Dependent Phase Separations of Intrinsically Disordered Proteins

Biomolecular condensates, physically underpinned to a significant extent by liquid-liquid phase separation (LLPS), are now widely recognized by numerous experimental studies to be of fundamental biological, biomedical, and biophysical importance. In the face of experimental discoveries, analytical formulations emerged as a powerful yet tractable tool in recent theoretical investigations of the role of LLPS in the assembly and dissociation of these condensates. The pertinent LLPS often involves, though not exclusively, intrinsically disordered proteins engaging in multivalent interactions that are governed by their amino acid sequences. For researchers interested in applying these theoretical methods, here we provide a practical guide to a set of computational techniques devised for extracting sequence-dependent LLPS properties from analytical formulations. The numerical procedures covered include those for the determination of spinodal and binodal phase boundaries from a general free energy function with examples based on the random phase approximation in polymer theory, construction of tie lines for multiple-component LLPS, and field-theoretic simulation of multiple-chain heteropolymeric systems using complex Langevin dynamics. Since a more accurate physical picture often requires comparing analytical theory against explicit-chain model predictions, a commonly utilized methodology for coarse-grained molecular dynamics simulations of sequence-specific LLPS is also briefly outlined.

Analytical Formulation and Field-Theoretic Simulation of Sequence-Specific Phase Separation of Protein-Like Heteropolymers with Short- and Long-Spatial-Range Interactions

A theory for sequence-dependent liquid-liquid phase separation (LLPS) of intrinsically disordered proteins (IDPs) in the study of biomolecular condensates is formulated by extending the random phase approximation (RPA) and field-theoretic simulation (FTS) of heteropolymers with spatially long-range Coulomb interactions to include the fundamental effects of short-range, hydrophobic-like interactions between amino acid residues. To this end, short-range effects are modeled by Yukawa interactions between multiple nonelectrostatic charges derived from an eigenvalue decomposition of pairwise residue-residue contact energies. Chain excluded volume is afforded by incompressibility constraints. A mean-field approximation leads to an effective Flory-Huggins χ parameter, which, in conjunction with RPA, accounts for the contact-interaction effects of amino acid composition and the sequence-pattern effects of long-range electrostatics in IDP LLPS, whereas FTS based on the formulation provides full sequence dependence for both short- and long-range interactions. This general approach is illustrated here by applications to variants of a natural IDP in the context of several different amino-acid interaction schemes as well as a set of different model hydrophobic-polar sequences sharing the same composition. Effectiveness of the methodology is verified by coarse-grained explicit-chain molecular dynamics simulations.

Phase Behavior of Ion-Containing Polymers in Polar Solvents: Predictions from a Liquid-State Theory with Local Short-Range Interactions

The thermodynamic phase behavior of charged polymers is a crucial property underlying their role in biology and various industrial applications. A complete understanding of the phase behaviors of such polymer solutions remains challenging due to the multi-component nature of the system and the delicate interplay among various factors, including the translational entropy of each component, excluded volume interactions, chain connectivity, electrostatic interactions, and other specific interactions. In this work, the phase behavior of partially charged ion-containing polymers in polar solvents is studied by further developing a liquid-state (LS) theory with local shortrange interactions. This work is based on the LS theory developed for fully-charged polyelectrolyte solutions. Specific interactions between charged groups of the polymer and counterions, between neutral segments of the polymer, and between charged segments of the polymer are incorporated into the LS theory by an extra Helmholtz free energy from the perturbed-chain statistical associating fluid theory (PC-SAFT). The influence of the sequence structure of the partially charged polymer is modeled by the number of connections between bonded segments. The effects of chain length, charge fraction, counterion valency, and specific short-range interactions are explored. A computational App for salt-free polymer solutions is developed and presented, which allows easy computation of the binodal curve and critical point by specifying values for the relevant model parameters.

A Perspective on Polythiophenes as Conformation Dependent Optical Reporters for Label-Free Bioanalytics

Poly(3-alkylthiophene) (PT)-based conjugated polyelectrolytes (CPEs) constitute an important class of responsive polymers with excellent optical properties. The electrostatic interactions between PTs and target analytes trigger complexation and concomitant conformational changes of the PT backbones that produce distinct optical responses. These conformation-induced optical responses of the PTs enable them to be utilized as reporters for detection of various analytes by employing simple UV-vis spectrophotometry or the naked eye. Numerous PTs with unique pendant groups have been synthesized to tailor their interactions with analytes such as nucleotides, ions, surfactants, proteins, and bacterial and viral pathogens. In this perspective, we discuss PT-target analyte complexation for bioanalytical applications and highlight recent advancements in point-of-care and field deployable assays. Subsequently, we highlight a few areas of critical importance for future applications of PTs as reporters, including (i) design and synthesis of specific PTs to advance the understanding of the mechanisms of interaction with target analytes, (ii) using arrays of PTs and linear discriminant analysis for selective and specific detection of target analytes, (iii) translation of conventional homogeneous solution-based assays into heterogeneous membrane-based assay formats, and finally (iv) the potential of using PT as an alternative to conjugated polymer nanoparticles and dots in bioimaging.

Extracting free energies of counterion binding to polyelectrolytes by molecular dynamics simulations

We use all-atom molecular dynamics simulations to extract ΔG^eff, the free energy of binding of potassium ions K⁺ to the partially charged polyelectrolyte poly(acrylic acid), or PAA, in dilute regimes. Upon increasing the charge fraction of PAA, the chains adopt more extended conformations, and simultaneously, potassium ions bind more strongly (i.e., with more negative ΔG^eff) to the highly charged chains to relieve electrostatic repulsions between charged monomers along the chains. We compare the simulation results with the predictions of a model that describes potassium binding to PAA chains as a reversible reaction whose binding free energy (ΔG^eff) is adjusted from its intrinsic value (ΔG) by electrostatic correlations, captured by a random phase approximation. The bare or intrinsic binding free energy ΔG, which is an input in the model, depends on the binding species and is obtained from the radial distribution function of K⁺ around the charged monomer of a singly charged, short PAA chain in dilute solutions. We find that the model yields semi-quantitative predictions for ΔG^eff and the degree of potassium binding to PAA chains, α, as a function of PAA charge fraction without using fitting parameters.

Molecularly Informed Field Theories from Bottom-up Coarse-Graining

Polymer formulations possessing mesostructures or phase coexistence are challenging to simulate using atomistic particle-explicit approaches due to the disparate time and length scales, while the predictive capability of field-based simulations is hampered by the need to specify interactions at a coarser scale (e.g., χ-parameters). To overcome the weaknesses of both, we introduce a bottom-up coarse-graining methodology that leverages all-atom molecular dynamics to molecularly inform coarser field-theoretic models. Specifically, we use relative-entropy coarse-graining to parametrize particle models that are directly and analytically transformable into statistical field theories. We demonstrate the predictive capability of this approach by reproducing experimental aqueous poly(ethylene oxide) (PEO) cloud-point curves with no parameters fit to experimental data. This synergistic approach to multiscale polymer simulations opens the door to de novo exploration of phase behavior across a wide variety of polymer solutions and melt formulations.

Subcompartmentalization of polyampholyte species in organelle-like condensates is promoted by charge-pattern mismatch and strong excluded-volume interaction

Polyampholyte field theory and explicit-chain molecular dynamics models of sequence-specific phase separation of a system with two intrinsically disordered protein (IDP) species indicate consistently that a substantial polymer excluded volume and a significant mismatch of the IDP sequence charge patterns can act in concert, but not in isolation, to demix the two IDP species upon condensation. This finding reveals an energetic-geometric interplay in a stochastic, "fuzzy" molecular recognition mechanism that may facilitate subcompartmentalization of membraneless organelles.

Microphase separation and aggregate self-assembly in brushes of oppositely charged polyelectrolytes triggered by ion pairing

This work applies a molecular theory to study the formation of lateral self-assembled aggregates in mixed brushes composed of polyanion and polycation chains. In order to overcome the well-known limitations of mean-field electrostatics to capture polyelectrolyte complexation, the formation of ion pairs between anionic and cationic groups in the polyelectrolytes is explicitly modeled in our theory as an association reaction. This feature is essential to capture the microphase separation of the mixed brush and the formation of lateral aggregates triggered by polyelectrolyte complexation. The effects of solution pH and ionic strength, surface coverage, and chain length on the morphology of the mixed brush are systematically explored. It is shown that increasing salt concentration leads to the rupture of polyelectrolyte complexes and the stabilization of the homogeneous, non-aggregated brush, providing that the formation of ion pairs between the polyelectrolytes and the salt ions in solution is explicitly accounted for by the theory. The inclusion of ion-pairing association reactions between oppositely charged polyelectrolytes within a mean-field description of electrostatics emerges from this work as a useful and simple theoretical approach to capture the formation of polyelectrolyte complexes and their responsiveness to solution ionic strength and pH.

Near-Critical Phase Behavior in Polyelectrolyte Solutions: Effect of Charge Fluctuations

The critical behaviors in polyelectrolyte (PE) solutions are studied by a renormalized Gaussian fluctuation theory. PEs are fully charged linear wormlike chains. The electrostatic interactions are considered in a continuum solvent while ignoring other interactions. The effects of temperature on the criticality are explored in a salt-free polyanion solution, and the effects of salt are explored in a symmetric mixed solution of polycation and polyanion. An unphysical phase coexistence (UPPC), in which a metastable dense phase coexists with an unstable dilute phase, always exists in the phase diagram. In the near-critical region, the UPPC could interfere with the real phase coexistence. The classical critical point is replaced by a "critical line", giving a flat top for the phase diagram. These behaviors are due to the effect of charge fluctuation in the near-critical region. Our results explain why the mean field approach overestimates the critical salt density by about 10% to experimental data in the study of coacervation between oppositely charged PEs.

Physical Principles Underlying the Complex Biology of Intracellular Phase Transitions

Many biomolecular condensates appear to form via spontaneous or driven processes that have the hallmarks of intracellular phase transitions. This suggests that a common underlying physical framework might govern the formation of functionally and compositionally unrelated biomolecular condensates. In this review, we summarize recent work that leverages a stickers-and-spacers framework adapted from the field of associative polymers for understanding how multivalent protein and RNA molecules drive phase transitions that give rise to biomolecular condensates. We discuss how the valence of stickers impacts the driving forces for condensate formation and elaborate on how stickers can be distinguished from spacers in different contexts. We touch on the impact of sticker- and spacer-mediated interactions on the rheological properties of condensates and show how the model can be mapped to known drivers of different types of biomolecular condensates.

A unified analytical theory of heteropolymers for sequence-specific phase behaviors of polyelectrolytes and polyampholytes

The physical chemistry of liquid-liquid phase separation (LLPS) of polymer solutions bears directly on the assembly of biologically functional dropletlike bodies from proteins and nucleic acids. These biomolecular condensates include certain extracellular materials and intracellular compartments that are characterized as "membraneless organelles." Analytical theories are a valuable, computationally efficient tool for addressing general principles. LLPS of neutral homopolymers is quite well described by theory, but it has been a challenge to develop general theories for the LLPS of heteropolymers involving charge-charge interactions. Here, we present a theory that combines a random-phase-approximation treatment of polymer density fluctuations and an account of intrachain conformational heterogeneity based on renormalized Kuhn lengths to provide predictions of LLPS properties as a function of pH, salt, and charge patterning along the chain sequence. Advancing beyond more limited analytical approaches, our LLPS theory is applicable to a wide variety of charged sequences ranging from highly charged polyelectrolytes to neutral or nearly neutral polyampholytes. This theory should be useful in high-throughput screening of protein and other sequences for their LLPS propensities and can serve as a basis for more comprehensive theories that incorporate nonelectrostatic interactions. Experimental ramifications of our theory are discussed.

Coarse-grained residue-based models of disordered protein condensates: utility and limitations of simple charge pattern parameters

Biomolecular condensates undergirded by phase separations of proteins and nucleic acids serve crucial biological functions. To gain physical insights into their genetic basis, we study how liquid-liquid phase separation (LLPS) of intrinsically disordered proteins (IDPs) depends on their sequence charge patterns using a continuum Langevin chain model wherein each amino acid residue is represented by a single bead. Charge patterns are characterized by the "blockiness" measure κ and the "sequence charge decoration" (SCD) parameter. Consistent with random phase approximation (RPA) theory and lattice simulations, LLPS propensity as characterized by critical temperature Tcr* increases with increasingly negative SCD for a set of sequences showing a positive correlation between κ and -SCD. Relative to RPA, the simulated sequence-dependent variation in Tcr* is often-though not always-smaller, whereas the simulated critical volume fractions are higher. However, for a set of sequences exhibiting an anti-correlation between κ and -SCD, the simulated Tcr*'s are quite insensitive to either parameter. Additionally, we find that blocky sequences that allow for strong electrostatic repulsion can lead to coexistence curves with upward concavity as stipulated by RPA, but the LLPS propensity of a strictly alternating charge sequence was likely overestimated by RPA and lattice models because interchain stabilization of this sequence requires spatial alignments that are difficult to achieve in real space. These results help delineate the utility and limitations of the charge pattern parameters and of RPA, pointing to further efforts necessary for rationalizing the newly observed subtleties.

Electrostatic correlations and the polyelectrolyte self energy

We address the effects of chain connectivity on electrostaticfluctuations in polyelectrolyte solutions using a field-theoretic, renormalizedGaussian fluctuation (RGF) theory. As in simple electrolyte solutions [Z.-G. Wang,Phys. Rev. E 81, 021501 (2010)], the RGF provides a unified theory forelectrostatic fluctuations, accounting for both dielectric and charge correlationeffects in terms of the self-energy. Unlike simple ions, the polyelectrolyte self energydepends intimately on the chain conformation, and our theory naturally provides aself-consistent determination of the response of intramolecular chain structure topolyelectrolyte and salt concentrations. The effects of the chain-conformation on theself-energy and thermodynamics are especially pronounced for flexiblepolyelectrolytes at low polymer and salt concentrations, where application of thewrong chain structure can lead to a drastic misestimation of the electrostaticcorrelations. By capturing the expected scaling behavior of chain size from dilute tosemi-dilute regimes, our theory provides improved estimates of the self energy at lowpolymer concentrations and correctly predicts the eventual N-independenceof the critical temperature and concentration of salt-free solutions of flexiblepolyelectrolytes. We show that the self energy can be interpreted in terms of aninfinite-dilution energy μ_m,0^el and a finite concentrationcorrelation correction μ^corr which tends to cancel out the formerwith increasing concentration.

An interplay of electrostatic and excluded volume interactions in the conformational behavior of a dipolar chain: theory and computer simulations

The effect of an interplay between electrostatic and excluded volume interactions on the conformational behavior of a dipolar chain has been studied theoretically and by means of molecular dynamics simulations. Every monomer unit of the dipolar chain comprises a dipole formed by a charged group of the chain and an oppositely charged counterion. The counterion is assumed to freely move around the chain but keeping the distance between oppositely charged ions (the dipole length) fixed. The novelty of the developed mean-field theory is that variations of the dipole parameters (the dipole length and the counterion size) have been accounted for in both electrostatic and excluded volume contributions to the total free energy of the dipolar chain. It has been shown that conformational transitions between swollen and collapsed states of the chain can be induced by fine-tuning the balance between electrostatic and excluded volume interactions. In particular, in low-polar media not only globule but also extended coil conformations can be realized even under strong electrostatic attraction. The results of MD simulations of a dipolar chain with variable dipolar length support theoretical conclusions.

50th Anniversary Perspective: A Perspective on Polyelectrolyte Solutions

From the beginning of life with the information-containing polymers until the present era of a plethora of water-based materials in health care industry and biotechnology, polyelectrolytes are ubiquitous with a broad range of structural and functional properties. The main attribute of polyelectrolyte solutions is that all molecules are strongly correlated both topologically and electrostatically in their neutralizing background of charged ions in highly polarizable solvent. These strong correlations and the necessary use of numerous variables in experiments on polyelectrolytes have presented immense challenges toward fundamental understanding of the various behaviors of charged polymeric systems. This Perspective presents the author's subjective summary of several conceptual advances and the remaining persistent challenges in the contexts of charge and size of polymers, structures in homogeneous solutions, thermodynamic instability and phase transitions, structural evolution with oppositely charged polymers, dynamics in polyelectrolyte solutions, kinetics of phase separation, mobility of charged macromolecules between compartments, and implications to biological systems.

Two regions of microphase separation in ion-containing polymer solutions

The phenomenon of spinodal decomposition in weakly charged polyelectrolyte solutions is studied theoretically within the random phase approximation. A novel feature of the theoretical approach is that it accounts for the effects of ionic association, i.e. ion pair and multiplet formation between counterions and ions in polymer chains, as well as the dependence of local dielectric permittivity on the polymer volume fraction Φ. The main focus is on the spinodal instability of polyelectrolyte solutions towards microscopic phase separation. It has been shown that increasing the binding energy of ions decreases the classical microphase separation region (possible at low polymer concentrations) due to the effective neutralization of the chains. A qualitatively new type of microphase separation is found in the presence of a dielectric mismatch between polymer and solvent. This new branch of microphase separation is realized at high polymer concentrations where ion association processes are the most pronounced. Typical microstructures are shown to have a period of a few nanometers like in ionomers. The driving force for the microphase formation of a new type is more favourable ion association in polymer-rich domains where ionomer-type behavior takes place. Effective attraction due to ion association promotes microscopic as well as macroscopic phase separation, even under good solvent conditions for uncharged monomer units of polymer chains. Polyelectrolyte-type behavior at low Φ and ionomer-type behavior at high Φ result in the presence of two critical points on the phase diagrams of polyelectrolyte solutions as well as two separate regions of possible microscopic structuring. Our predictions on the new type of microphase separation are supported by experimental data on polymer solutions, membranes and gels.

Explicit description of complexation between oppositely charged polyelectrolytes as an advantage of the random phase approximation over the scaling approach

Formation of single globules via 1 : 1 complexation of oppositely charged linear chains occurs prior to coacervation. F_corr is proved to be negative which is the difference between the random phase approximation (RPA) correction term and the self-energy of the chains.

Electrostatic Correlations in Polyelectrolyte Solutions

The major attribute of polyelectrolyte solutions is that all chains are strongly correlated both electrostatically and topologically. Even in very dilute solutions such that the chains are not interpenetrating, the chains are still strongly correlated. These correlations are manifest in the measured scattering intensity when such solutions are subjected to light, X-ray, and neutron radiation. The behavior of scattering intensity from polyelectrolyte solutions is qualitatively different from that of solutions of uncharged polymers. Using the technique introduced by Sir Sam Edwards, and extending the earlier work by the author on the thermodynamics of polyelectrolyte solutions, extrapolation formulas are derived for the scattering intensity from polyelectrolyte solutions. The emergence of the polyelectrolyte peak and its concentration dependence are derived. The derived theory shows that there are five regimes. Published experimental data from many laboratories are also collected into a master figure and a comparison between the present theory and experiments is presented.

Sequence-Specific Polyampholyte Phase Separation in Membraneless Organelles

Liquid-liquid phase separation of charge- and/or aromatic-enriched intrinsically disordered proteins (IDPs) is critical in the biological function of membraneless organelles. Much of the physics of this recent discovery remains to be elucidated. Here, we present a theory in the random phase approximation to account for electrostatic effects in polyampholyte phase separations, yielding predictions consistent with recent experiments on the IDP Ddx4. The theory is applicable to any charge pattern and thus provides a general analytical framework for studying sequence dependence of IDP phase separation.

Phase behavior of electrostatically complexed polyelectrolyte gels using an embedded fluctuation model

Nanostructured, responsive hydrogels formed due to electrostatic interactions have promise for applications such as drug delivery and tissue mimics. These physically cross-linked hydrogels are composed of an aqueous solution of oppositely charged triblocks with charged end-blocks and neutral, hydrophilic mid-blocks. Due to their electrostatic interactions, the end-blocks microphase separate and form physical cross-links that are bridged by the mid-blocks. The structure of this system was determined using a new, efficient embedded fluctuation (EF) model in conjunction with self-consistent field theory. The calculations using the EF model were validated against unapproximated field-theoretic simulations with complex Langevin sampling and were found consistent with small angle X-ray scattering (SAXS) measurements on an experimental system. Using both the EF model and SAXS, phase diagrams were generated as a function of end-block fraction and polymer concentration. Several structures were observed including a body-centered cubic sphere phase, a hexagonally packed cylinder phase, and a lamellar phase. Finally, the EF model was used to explore how parameters that directly relate to polymer chemistry can be tuned to modify the resulting phase diagram, which is of practical interest for the development of new hydrogels.