Molecular design of self-coacervation phenomena in block polyampholytes

Effects of charge asymmetry on the liquid-liquid phase separation of polyampholytes and their condensate properties

Liquid-liquid phase separation (LLPS) is the mechanism underlying the formation of bio-molecular condensates which are important compartments regulating intra- and extra-cellular functions. Electrostatic interactions are some of the important driving forces of the LLPS behaviors of biomolecules. However, the understanding of the electrostatic interactions is still limited, especially in the mixtures of biomolecules with different charge patterns. Here, we focus on the electrostatic interactions in mixtures of charge-asymmetric and charge-symmetric polyampholytes and their roles in the phase separation behaviors. We build charge-asymmetric and charge-symmetric model proteins consisting of both glutamic acid (E, negatively charged) and lysine (K, positively charged), i.e. polyampholytes of E35K15 (charge asymmetric) and E25K25 (charge symmetric). Pure E25K25 can undergo LLPS. To investigate the effects of charge-asymmetric polyampholytes on the mixtures of E25K25/E35K15, we perform coarse-grained simulations to determine their phase separation. The charge-asymmetric polyampholyte E35K15 is resistant to the LLPS of the mixtures of E25K25/E35K15. The condensate density decreases with the molar fraction of E35K15 increasing to 0.4, and no LLPS occurs at the molar fraction of 0.5 and above. This can be attributed to the electrostatic repulsion between the negatively charged E35K15 polymers. We further investigate the effects of charge asymmetry on the conformations and properties of the condensates. The E35K15 polymers in the condensates exhibit a more collapsed state as the molar fraction of E35K15 increases. However, the conformation of E25K25 polymers changes slightly across different condensates. The surface tensions of condensates decline with the increase of the molar fraction of E35K15 polymers, while the diffusivity of polymers in the condensed phases is enhanced. This work elucidates the role of charge-asymmetric polyampholytes in determining the LLPS behaviours of binary mixtures of charge-symmetric and charge-asymmetric proteins as well as the properties of condensed phases.

Effects of Charge Sequence Pattern and Lysine-to-Arginine Substitution on the Structural Stability of Bioinspired Polyampholytes

A comprehensive study focusing on the combined influence of the charge sequence pattern and the type of positively charged amino acids on the formation of secondary structures in sequence-specific polyampholytes is presented. The sequences of interest consisting exclusively of ionizable amino acids (lysine, K; arginine, R; and glutamic acid, E) are (EKEK)5, (EKKE)5, (ERER)5, (ERRE)5, and (EKER)5. The stability of the secondary structure was examined at three pH values in the presence of urea and NaCl. The results presented here underscore the combined prominent effects of the charge sequence pattern and the type of positively charged monomers on secondary structure formation. Additionally, (ERRE)5 readily aggregated across a wide range of pH. In contrast, sequences with the same charge pattern, (EKKE)5, as well as the sequences with the equivalent amino acid content, (ERER)5, exhibited no aggregate formation under equivalent pH and concentration conditions.

Sequence-Tunable Phase Behavior and Intrinsic Fluorescence in Dynamically Interacting Peptides

A conceptual framework towards understanding biological condensed phases is emerging, derived from biological, biomimetic, and synthetic sequences. However, de novo peptide condensate design remains a challenge due to an incomplete understanding of the structural and interactive complexity. We designed peptide modules based on a simple repeat motif composed of tripeptide spacers (GSG, SGS, GLG) interspersed with adhesive amino acids (R/H and Y). We show, using sequence editing and a combination of computation and experiment, that n→π* interactions in GLG backbones are a dominant factor in providing sufficient backbone structure, which in turn regulates the water interface, collectively promoting liquid droplet formation. Moreover, these R(GLG)Y and H(GLG)Y condensates unexpectedly display sequence-dependent emission that is a consequence of their non-covalent network interactions, and readily observable by confocal microscopy.

Polyelectrolytes: From Seminal Works to the Influence of the Charge Sequence

We propose a selected tour of the physics of polyelectrolytes (PE) following the line initiated by de Gennes and coworkers in their seminal 1976 paper. The early works which used uniform charge distributions along the PE backbone achieved tremendous progress and set most milestones in the field. Recently, the focus has shifted to the role of the charge sequence. Revisited topics include PE complexation and polyampholytes (PA). We develop the example of a random PE in poor solvent forming pearl-necklace structures. It is shown that the pearls typically adopt very asymmetric mass and charge distributions. Individual sequences do not necessarily reflect the ensemble statistics and a rich variety of behaviors emerges (specially for PA). Pearl necklaces are dynamic structures and switch between various types of pearl-necklace structures, as described for both PE and PA.

Assessment of the partial saddle point approximation in field-theoretic polymer simulations

Field-theoretic simulations are numerical treatments of polymer field theory models that go beyond the mean-field self-consistent field theory level and have successfully captured a range of mesoscopic phenomena. Inherent in molecularly-based field theories is a "sign problem" associated with complex-valued Hamiltonian functionals. One route to field-theoretic simulations utilizes the complex Langevin (CL) method to importance sample complex-valued field configurations to bypass the sign problem. Although CL is exact in principle, it can be difficult to stabilize in strongly fluctuating systems. An alternate approach for blends or block copolymers with two segment species is to make a "partial saddle point approximation" (PSPA) in which the stiff pressure-like field is constrained to its mean-field value, eliminating the sign problem in the remaining field theory, allowing for traditional (real) sampling methods. The consequences of the PSPA are relatively unknown, and direct comparisons between the two methods are limited. Here, we quantitatively compare thermodynamic observables, order-disorder transitions, and periodic domain sizes predicted by the two approaches for a weakly compressible model of AB diblock copolymers. Using Gaussian fluctuation analysis, we validate our simulation observations, finding that the PSPA incorrectly captures trends in fluctuation corrections to certain thermodynamic observables, microdomain spacing, and location of order-disorder transitions. For incompressible models with contact interactions, we find similar discrepancies between the predictions of CL and PSPA, but these can be minimized by regularization procedures such as Morse calibration. These findings mandate caution in applying the PSPA to broader classes of soft-matter models and systems.

The effect of monomer polarizability on the stability and salt partitioning in model coacervates

Coacervation of charged polymer chains has been a topic of major interest in both polymer and biological sciences, as it is a subset of a phenomenon called liquid-liquid phase separation (LLPS). In this process the polymer-rich phase separates from the polymer-lean supernatant while still maintaining its liquid-like properties. LLPS has been shown to play a crucial role in cellular homeostasis by driving the formation of membraneless organelles. It also has the potential to be harnessed to aid in novel therapeutical applications. Recent studies have demonstrated that there is no one simple mechanism which drives LLPS, which is instead a result of the combined effect of electrostatic, dipolar, hydrophobic, and other weak interactions. Using coarse-grained polymer simulations we investigate the relatively unexplored effects of monomer polarizability and spatially varying dielectric constant on LLPS propensity, and these factors affect the properties of the resulting condensates. In order to produce spatial variations in the dielectric constant, all our simulations include explicit solvent and counterions. We demonstrate that polarizability has only a minor effect on the bulk behaviour of the condensates but plays a major role when ion partitioning and microstructure are considered. We observe that the major contribution comes from the nature of the neutral blocks as endowing them with an induced dipole changes their character from hydrophobic to hydrophilic. We hypothesize that the results of this work can aid in guiding future studies concerned with LLPS by providing a general framework and by highlighting important factors which influence LLPS.

Competition between β-Sheet and Coacervate Domains Yields Diverse Morphologies in Mixtures of Oppositely Charged Homochiral Polypeptides

Biomolecular assembly processes involving competition between specific intermolecular interactions and thermodynamic phase instability have been implicated in a number of pathological states and technological applications of biomaterials. As a model for such processes, aqueous mixtures of oppositely charged homochiral polypeptides such as poly-l-lysine and poly-l-glutamic acid have been reported to form either β-sheet-rich solid-like precipitates or liquid-like coacervate droplets depending on competing hydrogen bonding interactions. Herein, we report studies of polypeptide mixtures that reveal unexpectedly diverse morphologies ranging from partially coalescing and aggregated droplets to bulk precipitates, as well as a previously unreported re-entrant liquid-liquid phase separation at high polypeptide concentration and ionic strength. Combining our experimental results with all-atom molecular dynamics simulations of folded polypeptide complexes reveals a concentration dependence of β-sheet-rich secondary structure, whose relative composition correlates with the observed macroscale morphologies of the mixtures. These results elucidate a crucial balance of interactions that are important for controlling morphology during coacervation in these and potentially similar biologically relevant systems.

Combining particle and field-theoretic polymer models with multi-representation simulations

Particle-based and field-theoretic simulations are both widely used methods to predict the properties of polymeric materials. In general, the advantages of each method are complementary. Field-theoretic simulations are preferred for polymers with high molecular weights and can provide direct access to chemical potentials and free energies, which makes them the method-of-choice for calculating phase diagrams. The trade-off is that field-theoretic simulations sacrifice the molecular details present in particle-based simulations, such as the configurations of individual molecules and their dynamics. In this work, we describe a new approach to conduct "multi-representation" simulations that efficiently map between particle-based and field-theoretic simulations. Our approach involves the construction of formally equivalent particle-based and field-based models, which are then simulated subject to the constraint that their spatial density profiles are equal. This constraint provides the ability to directly link particle-based and field-based simulations and enables calculations that can switch between one representation to the other. By switching between particle/field representations during a simulation, we demonstrate that our approach can leverage many of the advantages of each representation while avoiding their respective limitations. Although our method is illustrated in the context of complex sphere phases in linear diblock copolymers, we anticipate that it will be useful whenever free energies, rapid equilibration, molecular configurations, and dynamic information are all simultaneously desired.

Influence of Charge Block Length on Conformation and Solution Behavior of Polyampholytes

We investigate the effect of charge block length on polyampholyte chain conformation and phase behavior using small-angle X-ray scattering (SAXS) and implicit-solvent molecular simulations. To this end, we use solid phase peptide synthesis to precision-tailor a series of polyampholytes consisting of l-glutamic acid (E) and l-lysine (K) monomers arranged in alternating blocks from 2 to 16 monomers. We observe that the polyampholytes tend to phase separate as block size increases. With addition of NaCl, phase separated polyampholytes exhibit a salting-in effect dependent on charge block length. Fourier-transform infrared (FTIR) spectroscopy reveals the presence of intramolecular hydrogen bonds that are disrupted upon the addition of NaCl, implicating both electrostatic interactions and hydrogen bonding in the phase behavior. SAXS spectra at no-added salt conditions show minimal dependence of charge block length on the radius of gyration (R_g) for soluble polyampholytes, but local chain stiffening is found to be dependent on charge block length. With increasing NaCl, consistent with electrostatic screening, all polyampholytes expand and behave as neutral or swollen chains in good solvent conditions. Molecular simulations are qualitatively consistent with experiments. Implications for understanding intracellular condensates and material design are noted.

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.

Field-Theoretic Simulation Method to Study the Liquid-Liquid Phase Separation of Polymers

Liquid-liquid phase separation (LLPS) is a process that results in the formation of a polymer-rich liquid phase coexisting with a polymer-depleted liquid phase. LLPS plays a critical role in the cell through the formation of membrane-less organelles, but it also has a number of biotechnical and biomedical applications such as drug confinement and its targeted delivery. In this chapter, we present a computational efficient methodology that uses field-theoretic simulations (FTS) with complex Langevin (CL) sampling to characterize polymer phase behavior and delineate the LLPS phase boundaries. This approach is a powerful complement to analytical and explicit-particle simulations, and it can serve to inform experimental LLPS studies. The strength of the method lies in its ability to properly sample a large ensemble of polymers in a saturated solution while including the effect of composition fluctuations on LLPS. We describe the approaches that can be used to accurately construct phase diagrams of a variety of molecularly designed polymers and illustrate the method by generating an approximation-free phase diagram for a classical symmetric diblock polyampholyte.

Multiscale Modeling of Protein-RNA Condensation in and Out of Equilibrium

A vast number of intracellular membraneless bodies also known as biomolecular condensates form through a liquid-liquid phase separation (LLPS) of biomolecules. To date, phase separation has been identified as the main driving force for a membraneless organelles such as nucleoli, Cajal bodies, stress granules, and chromatin compartments. Recently, the protein-RNA condensation is receiving increased attention, because it is closely related to the biological function of cells such as transcription, translation, and RNA metabolism. Despite the multidisciplinary efforts put forth to study the biophysical properties of protein-RNA condensates, there are many fundamental unanswered questions regarding the mechanism of formation and regulation of protein-RNA condensates in eukaryotic cells. Major challenges in studying protein-RNA condensation stem from (i) the molecular heterogeneity and conformational flexibility of RNA and protein chains and (ii) the nonequilibrium nature of transcription and cellular environment. Computer simulations, bioinformatics, and mathematical models are uniquely positioned for shedding light on the microscopic nature of protein-RNA phase separation. To this end, there is an urgent need for innovative models with the right spatiotemporal resolution for confronting the experimental observables in a comprehensive and physics-based manner. In this chapter, we will summarize the currently emerging research efforts, which employ atomistic and coarse-grained molecular models and field theoretical models to understand equilibrium and nonequilibrium aspects of protein-RNA condensation.

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.

Sequence-Controlled Secondary Structures and Stimuli Responsiveness of Bioinspired Polyampholytes

A comprehensive study focusing on the influence of the sequence charge pattern on the secondary structure preferences of annealed polyampholytes and their responsiveness to external stimuli is presented. Two sequences are designed composed entirely of ionizable amino acids (charge fraction, f = 1) and an equal number of positive and negative charges (f+ = f- = 0.5) with distinct charge patterns consisting of lysine and glutamic acid monomers. The study reveals that the sequence charge pattern has a significant influence on the secondary structure preferences of polyampholytes at physiological pH. Furthermore, it shows that external stimuli such as pH, ionic strength, and solvent dielectric constant can be used to modulate the secondary structure of the two studied sequences. The observed secondary structure transformations for the two sequences are also substantially different from those determined for uniformly charged homo-polypeptides under matching conditions.

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.

Transferable and Polarizable Coarse Grained Model for Proteins─ProMPT

The application of classical molecular dynamics (MD) simulations at atomic resolution (fine-grained level, FG), to most biomolecular processes, remains limited because of the associated computational complexity of representing all the atoms. This problem is magnified in the presence of protein-based biomolecular systems that have a very large conformational space, and MD simulations with fine-grained resolution have slow dynamics to explore this space. Current transferable coarse grained (CG) force fields in literature are either limited to only peptides with the environment encoded in an implicit form or cannot capture transitions into secondary/tertiary peptide structures from a primary sequence of amino acids. In this work, we present a transferable CG force field with an explicit representation of the environment for accurate simulations with proteins. The force field consists of a set of pseudoatoms representing different chemical groups that can be joined/associated together to create different biomolecular systems. This preserves the transferability of the force field to multiple environments and simulation conditions. We have added electronic polarization that can respond to environmental heterogeneity/fluctuations and couple it to protein's structural transitions. The nonbonded interactions are parametrized with physics-based features such as solvation and partitioning free energies determined by thermodynamic calculations and matched with experiments and/or atomistic simulations. The bonded potentials are inferred from corresponding distributions in nonredundant protein structure databases. We present validations of the CG model with simulations of well-studied aqueous protein systems with specific protein fold types─Trp-cage, Trpzip4, villin, WW-domain, and β-α-β. We also explore the applications of the force field to study aqueous aggregation of Aβ 16-22 peptides.

Direct calculation of the functional inverse of realistic interatomic potentials in field-theoretic simulations

We discuss the functional inverse problem in field-theoretic simulations for realistic pairwise potentials such as the Morse potential (widely used in particle simulations as an alternative to the 12-6 Lennard-Jones one), and we propose the following two solutions: (a) a numerical one based on direct inversion on a regular grid or deconvolution and (b) an analytical one by expressing attractive and repulsive contributions to the Morse potential as higher-order derivatives of the Dirac delta function; the resulting system of ordinary differential equations in the saddle-point approximation is solved numerically with appropriate model-consistent boundary conditions using a Newton-Raphson method. For the first time, exponential-like, physically realistic pair interactions are analytically treated and incorporated into a field-theoretic framework. The advantages and disadvantages of the two approaches are discussed in detail in connection with numerical findings from test simulations for the radial distribution function of a monatomic fluid at realistic densities providing direct evidence for the capability of the analytical method to resolve structural features down to the Angstrom scale.

Rules of Physical Mathematics Govern Intrinsically Disordered Proteins

In stark contrast to foldable proteins with a unique folded state, intrinsically disordered proteins and regions (IDPs) persist in perpetually disordered ensembles. Yet an IDP ensemble has conformational features-even when averaged-that are specific to its sequence. In fact, subtle changes in an IDP sequence can modulate its conformational features and its function. Recent advances in theoretical physics reveal a set of elegant mathematical expressions that describe the intricate relationships among IDP sequences, their ensemble conformations, and the regulation of their biological functions. These equations also describe the molecular properties of IDP sequences that predict similarities and dissimilarities in their functions and facilitate classification of sequences by function, an unmet challenge to traditional bioinformatics. These physical sequence-patterning metrics offer a promising new avenue for advancing synthetic biology at a time when multiple novel functional modes mediated by IDPs are emerging.

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.

Limitations of field-theory simulation for exploring phase separation: The role of repulsion in a lattice protein model

Field-theory simulation by the complex Langevin method offers an alternative to conventional sampling techniques for exploring the forces driving biomolecular liquid-liquid phase separation. Such simulations have recently been used to study several polyampholyte systems. Here, we formulate a field theory corresponding to the hydrophobic/polar (HP) lattice protein model, with finite same-site repulsion and nearest-neighbor attraction between HH bead pairs. By direct comparison with particle-based Monte Carlo simulations, we show that complex Langevin sampling of the field theory reproduces the thermodynamic properties of the HP model only if the same-site repulsion is not too strong. Unfortunately, the repulsion has to be taken weaker than what is needed to prevent condensed droplets from assuming an artificially compact shape. Analysis of a minimal and analytically solvable toy model hints that the sampling problems caused by repulsive interaction may stem from loss of ergodicity.

Dynamics of a single polyampholyte chain

Polymers that feature both positive and negative charges along chains, known as polyampholytes, represent a class of materials that hold promise for a new generation of energy storage devices, the design of which will require knowledge of the underlying structure and dynamics. Here, we develop a theory based on the Rouse model for the dynamic structure factor of a single polyampholyte chain in the weak coupling regime (negligible intramolecular electrostatics) or subjected to weak external electric fields (governed by linear response). Neglecting effects of small ions, we find deviations in scaling from the classic Rouse theory and make predictions for scattering experiments performed on polyampholytes. We find that, under weak coupling with arbitrarily strong fields, the dynamics are highly dependent on the charge distribution and consequently look at two representative examples-random charge densities and periodic charge densities-with different scaling properties. Under weak fields, the dynamics are largely independent of charge distribution. Finally, we investigate the influence of hydrodynamic effects and the implications of including inertial effects in the model.

Peptide Design and Self-assembly into Targeted Nanostructure and Functional Materials

Peptides have been extensively utilized to construct nanomaterials that display targeted structure through hierarchical assembly. The self-assembly of both rationally designed peptides derived from naturally occurring domains in proteins as well as intuitively or computationally designed peptides that form β-sheets and helical secondary structures have been widely successful in constructing nanoscale morphologies with well-defined 1-d, 2-d, and 3-d architectures. In this review, we discuss these successes of peptide self-assembly, especially in the context of designing hierarchical materials. In particular, we emphasize the differences in the level of peptide design as an indicator of complexity within the targeted self-assembled materials and highlight future avenues for scientific and technological advances in this field.

The Protein Folding Problem: The Role of Theory

The protein folding problem was first articulated as question of how order arose from disorder in proteins: How did the various native structures of proteins arise from interatomic driving forces encoded within their amino acid sequences, and how did they fold so fast? These matters have now been largely resolved by theory and statistical mechanics combined with experiments. There are general principles. Chain randomness is overcome by solvation-based codes. And in the needle-in-a-haystack metaphor, native states are found efficiently because protein haystacks (conformational ensembles) are funnel-shaped. Order-disorder theory has now grown to encompass a large swath of protein physical science across biology.

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.

Multiphase Coacervates Driven by Electrostatic Correlations

The liquid-liquid phase separation of a polyelectrolyte solution containing one type of negatively and two types of positively charged polymers with different charge densities is studied theoretically by random phase approximation (RPA). It is predicted that multicoacervate phases could coexist, driven purely by electrostatic correlations. The asymmetry of the linear charge density could induce an effective immiscibility between two positively charged polyelectrolytes, leading to the multiphase separation. Adding salt will induce the disappearance of the dilute phase, forming two coexisting complex phases, instead of fusion between coacervates. Raising temperature could either induce a two coexisting complex phase, or a dilute phase coexisting with a coacervate phase, depending on the bulk concentration. Our predictions are in good agreement with experiments and provide insights in the further designing of the multiphase coacervation system.

Sequence Blockiness Controls the Structure of Polyampholyte Necklaces

A scaling theory of statistical (Markov) polyampholytes is developed to understand how sequence correlations, that is, the blockiness of positive and negative charges, influences conformational behavior. An increase in the charge patchiness leads to stronger correlation attractions between oppositely charged monomers, but simultaneously, it creates a higher charge imbalance in the polyampholyte. A competition between effective short-range attractions and long-range Coulomb repulsions induces globular, pearl-necklace, or fully stretched chain conformations, depending on the average length of the block of like charges. The necklace structure and the underlying distribution of the net charge are also controlled by the sequence. Sufficiently long blocks allow for charge migration from globular beads (pearls) to strings, thereby providing a nonmonotonic change in the number of necklace beads as the blockiness increases. The sequence-dependent structure of polyampholyte necklaces is confirmed by molecular dynamics simulations. The findings presented here provide a framework for understanding the sequence-encoded conformations of synthetic polyampholytes and intrinsically disordered proteins (IDPs).

Transfer Matrix Model of pH Effects in Polymeric Complex Coacervation

Oppositely charged polyelectrolytes can undergo an associative phase separation, in a process known as polymeric complex coacervation. This phenomenon is driven by the electrostatic attraction between polyanion and polycation species, leading to the formation of a polymer-dense coacervate phase and a coexisting polymer-dilute supernatant phase. There has been significant recent interest in the physical origin and features of coacervation; yet notably, experiments often use weak polyelectrolytes the charge state of which depends on solution pH, while theoretical or computational efforts typically assume strong polyelectrolytes that remain fully charged. There have been only a few efforts to address this limitation, and thus there has been little exploration of how pH can affect complex coacervation. In this paper, we modify a transfer matrix theory of coacervation to account for acid-base equilibria, taking advantage of its ability to directly account for some local ion correlations that will affect monomer charging. We show that coacervation can stabilize the charged state of a weak polyelectrolyte via the proximity of oppositely charged monomers, and can lead to asymmetric phase diagrams where the positively and negatively charged polyelectrolytes exhibit different behaviors near the pK_a of either chain. Specifically, there is a partitioning of one of the salt species to a coacervate to maintain electroneutrality when one of the polyelectrolytes is only partially charged. This results in the depletion of the same salt species in the supernatant, and overall can suppress phase separation. We also demonstrate that, when one of the species is only partially charged, mixtures that are off-stoichiometric in volume fraction but stoichiometric in charge exhibit the greatest propensity to form coacervate phases.

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.

A Simple Explicit-Solvent Model of Polyampholyte Phase Behaviors and Its Ramifications for Dielectric Effects in Biomolecular Condensates

Biomolecular condensates such as membraneless organelles, underpinned by liquid-liquid phase separation (LLPS), are important for physiological function, with electrostatics, among other interaction types, being a prominent force in their assembly. Charge interactions of intrinsically disordered proteins (IDPs) and other biomolecules are sensitive to the aqueous dielectric environment. Because the relative permittivity of protein is significantly lower than that of water, the interior of an IDP condensate is expected to be a relatively low-dielectric regime, which aside from its possible functional effects on client molecules should facilitate stronger electrostatic interactions among the scaffold IDPs. To gain insight into this LLPS-induced dielectric heterogeneity, addressing in particular whether a low-dielectric condensed phase entails more favorable LLPS than that posited by assuming IDP electrostatic interactions are uniformly modulated by the higher dielectric constant of the pure solvent, we consider a simplified multiple-chain model of polyampholytes immersed in explicit solvents that are either polarizable or possess a permanent dipole. Notably, simulated phase behaviors of these systems exhibit only minor to moderate differences from those obtained using implicit-solvent models with a uniform relative permittivity equals to that of pure solvent. Buttressed by theoretical treatments developed here using random phase approximation and polymer field-theoretic simulations, these observations indicate a partial compensation of effects between favorable solvent-mediated interactions among the polyampholytes in the condensed phase and favorable polyampholyte-solvent interactions in the dilute phase, often netting only a minor enhancement of overall LLPS propensity from the very dielectric heterogeneity that arises from the LLPS itself. Further ramifications of this principle are discussed.

The effect of ion pairs on coacervate-driven self-assembly of block polyelectrolytes

The incorporation of oppositely charged polyelectrolytes into a block copolymer system can lead to formation of microphase separated nanostructures driven by the electrostatic complex between two oppositely charged blocks. It is a theoretical challenge to build an appropriate model to handle such coacervate-driven self-assembly, which should capture the strong electrostatic correlations for highly charged polymers. In this paper, we develop the self-consistent field theory considering the ion paring effect to predict the phase behavior of block polyelectrolytes. In our model, two types of ion pairs, the binding between two oppositely charged monomers and the binding between charged monomers and counterions, are included. Their strength of formation is controlled by two parameters K_aa and K_ac, respectively. We give a detailed analysis about how the binding strength K_ac and K_aa and salt concentration affect the self-assembled nanostructure of diblock polyelectrolyte systems. The results show that the binding between two oppositely charged blocks provides driven force for microphase separation, while the binding between charged monomers and counterions competes with the polyion pairing and thus suppresses the microphase separation. The addition of salt has a shielding effect on the charges of polymers, which is a disadvantage to microphase separation. The phase diagrams as a function of polymer concentration and salt concentration at different situations are constructed, and the influence of K_aa, K_ac, and charged block composition f_a is analyzed in depth. The obtained phase diagrams are in good agreement with currently existing experimental and theoretical results.

Surface science of cosmetic substrates, cleansing actives and formulations

The development of shampoo and cleansing formulations in cosmetics is at a crossroads due to consumer demands for better performing, more natural products and also the strong commitment of cosmetic companies to improve the sustainability of cosmetic products. In order to go beyond traditional formulations, it is of great importance to clearly establish the science behind cleansing technologies and appreciate the specificity of cleansing biological surfaces such as hair and skin. In this review, we present recent advances in our knowledge of the physicochemical properties of the hair surface from both an experimental and a theoretical point of view. We discuss the opportunities and challenges that newer, sustainable formulations bring compared to petroleum-based ingredients. The inevitable evolution towards more bio-based, eco-friendly ingredients and sustainable formulations requires a complete rethink of many well-known physicochemical principles. The pivotal role of digital sciences and modelling in the understanding and conception of new ingredients and formulations is discussed. We describe recent numerical approaches that take into account the specificities of the hair surface in terms of structuration, different methods that study the adsorption of formulation ingredients and finally the success of new data-driven approaches. We conclude with practical examples on current formulation efforts incorporating bio-surfactants, controlling foaming and searching for new rheological properties.

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.

Biomolecular Condensates: Sequence Determinants of Phase Separation, Microstructural Organization, Enzymatic Activity, and Material Properties

This perspective article highlights recent progress and emerging challenges in understanding the formation and function of membraneless organelles (MLOs). A long-term goal in the MLO field is to identify the sequence-encoded rules that dictate the formation of compositionally controlled biomolecular condensates, which cells utilize to perform a wide variety of functions. The molecular organization of the different components within a condensate can vary significantly, ranging from a homogeneous mixture to core-shell droplet structures. We provide many examples to highlight the richness of the observed behavior and potential research directions for improving our mechanistic understanding. The tunable environment within condensates can, in principle, alter enzymatic activity significantly. We examine recent examples where this was demonstrated, including applications in synthetic biology. An important question about MLOs is the role of liquid-like material properties in biological function. We discuss the need for improved quantitative characterization tools and the development of sequence-structure-dynamics relationships.

Fundamental Challenges and Outlook in Simulating Liquid-Liquid Phase Separation of Intrinsically Disordered Proteins

Intrinsically disordered proteins (IDPs) populate an ensemble of dynamic conformations, making their structural characterization by experiments challenging. Many IDPs undergo liquid-liquid phase separation into dense membraneless organelles with myriad cellular functions. Multivalent interactions in low-complexity IDPs promote the formation of these subcellular coacervates. While solution NMR, Förster resonance energy transfer (FRET), and small-angle X-ray scattering (SAXS) studies on IDPs have their own challenges, recent computational methods draw a rational trade-off to characterize the driving forces underlying phase separation. In this Perspective, we critically evaluate the scope of approximation-free field theoretic simulations, well-tempered ensemble methods, enhanced sampling techniques, coarse-grained force fields, and slab simulation approaches to offer an improved understanding of phase separation. A synergy between simulation length scale and model resolution would reduce the existing caveats and enable theories of polymer physics to elucidate finer details of liquid-liquid phase separation (LLPS). These computational advances offer promise for rigorous characterization of the IDP proteome and designing peptides with tunable material and self-assembly properties.