Lattice cluster theory of associating polymers. I. Solutions of linear telechelic polymer chains

Effect of End Groups on the Cloud Point Temperature of Aqueous Solutions of Thermoresponsive Polymers: An Inside View by Flory-Huggins Theory

In an effort to gain insight into the origin of the effects of end groups on the cloud point temperature (Tcp) as a function of the polymer molar mass of thermoresponsive polymers with lower critical solution behavior in dilute aqueous solutions, we use the Flory-Huggins (FH) theory amended for end groups. The theory was applied to available experimental data sets of poly(N-isopropylacrylamide) (PNIPAM), poly(4-vinylbenzyl methoxytris(oxyethylene) ether) (PTEGSt), and poly(α-hydro-ω-(4-vinylbenzyl)tetrakis(oxyethylene) ether) (PHTrEGSt). The theory relates the variations in TcpM,ϕcp for different end groups to the effective FH χ parameter of the end groups and explains the qualitative notion that the influence of the end groups is related to the hydrophobicity/hydrophilicity of the end groups relative to that of the so called intrinsic TcpM,ϕcp response of a polymer without end groups. The limits to the applicability of the FH theory are established, and a set of possible theoretical improvements is considered. The ultimate scrutiny of the simple FH theory and suggested improved theories must await the measurement of truly thermodynamic cloud points; the available cloud points are merely estimations of the thermodynamic cloud point, for which the deviation to the true cloud point cannot be established with sufficient accuracy.

Macromolecular Crowding Is Surprisingly Unable to Deform the Structure of a Model Biomolecular Condensate

The crowded interior of a living cell makes performing experiments on simpler in vitro systems attractive. Although these reveal interesting phenomena, their biological relevance can be questionable. A topical example is the phase separation of intrinsically disordered proteins into biomolecular condensates, which is proposed to underlie the membrane-less compartmentalization of many cellular functions. How a cell reliably controls biochemical reactions in compartments open to the compositionally-varying cytoplasm is an important question for understanding cellular homeostasis. Computer simulations are often used to study the phase behavior of model biomolecular condensates, but the number of relevant parameters increases as the number of protein components increases. It is unfeasible to exhaustively simulate such models for all parameter combinations, although interesting phenomena are almost certainly hidden in their high-dimensional parameter space. Here, we have studied the phase behavior of a model biomolecular condensate in the presence of a polymeric crowding agent. We used a novel compute framework to execute dozens of simultaneous simulations spanning the protein/crowder concentration space. We then combined the results into a graphical representation for human interpretation, which provided an efficient way to search the model's high-dimensional parameter space. We found that steric repulsion from the crowder drives a near-critical system across the phase boundary, but the molecular arrangement within the resulting biomolecular condensate is rather insensitive to the crowder concentration and molecular weight. We propose that a cell may use the local cytoplasmic concentration to assist the formation of biomolecular condensates, while relying on the dense phase to reliably provide a stable, structured, fluid milieu for cellular biochemistry despite being open to its changing environment.

Model biomolecular condensates have heterogeneous structure quantitatively dependent on the interaction profile of their constituent macromolecules

Biomolecular condensates play numerous roles in cells by selectively concentrating client proteins while excluding others. These functions are likely to be sensitive to the spatial organization of the scaffold proteins forming the condensate. We use coarse-grained molecular simulations to show that model intrinsically-disordered proteins phase separate into a heterogeneous, structured fluid characterized by a well-defined length scale. The proteins are modelled as semi-flexible polymers with punctate, multifunctional binding sites in good solvent conditions. Their dense phase is highly solvated with a spatial structure that is more sensitive to the separation of the binding sites than their affinity. We introduce graph theoretic measures to quantify their heterogeneity, and find that it increases with increasing binding site number, and exhibits multi-timescale dynamics. The model proteins also swell on passing from the dilute solution to the dense phase. The simulations predict that the structure of the dense phase is modulated by the location and affinity of binding sites distant from the termini of the proteins, while sites near the termini more strongly affect its phase behaviour. The relations uncovered between the arrangement of weak interaction sites on disordered proteins and the material properties of their dense phase can be experimentally tested to give insight into the biophysical properties, pathological effects, and rational design of biomolecular condensates.

Self-assembly and glass-formation in a lattice model of telechelic polymer melts: Influence of stiffness of the sticky bonds

Telechelic polymers are chain macromolecules that may self-assemble through the association of their two mono-functional end groups (called "stickers"). A deep understanding of the relation between microscopic molecular details and the macroscopic physical properties of telechelic polymers is important in guiding the rational design of telechelic polymer materials with desired properties. The lattice cluster theory (LCT) for strongly interacting, self-assembling telechelic polymers provides a theoretical tool that enables establishing the connections between important microscopic molecular details of self-assembling polymers and their bulk thermodynamics. The original LCT for self-assembly of telechelic polymers considers a model of fully flexible linear chains [J. Dudowicz and K. F. Freed, J. Chem. Phys. 136, 064902 (2012)], while our recent work introduces a significant improvement to the LCT by including a description of chain semiflexibility for the bonds within each individual telechelic chain [W.-S. Xu and K. F. Freed, J. Chem. Phys. 143, 024901 (2015)], but the physically associative (or called "sticky") bonds between the ends of the telechelics are left as fully flexible. Motivated by the ubiquitous presence of steric constraints on the association of real telechelic polymers that impart an additional degree of bond stiffness (or rigidity), the present paper further extends the LCT to permit the sticky bonds to be semiflexible but to have a stiffness differing from that within each telechelic chain. An analytical expression for the Helmholtz free energy is provided for this model of linear telechelic polymer melts, and illustrative calculations demonstrate the significant influence of the stiffness of the sticky bonds on the self-assembly and thermodynamics of telechelic polymers. A brief discussion is also provided for the impact of self-assembly on glass-formation by combining the LCT description for this extended model of telechelic polymers with the Adam-Gibbs relation between the structural relaxation time and the configurational entropy.

Lattice model of linear telechelic polymer melts. II. Influence of chain stiffness on basic thermodynamic properties

The lattice cluster theory (LCT) for semiflexible linear telechelic melts, developed in Paper I, is applied to examine the influence of chain stiffness on the average degree of self-assembly and the basic thermodynamic properties of linear telechelic polymer melts. Our calculations imply that chain stiffness promotes self-assembly of linear telechelic polymer melts that assemble on cooling when either polymer volume fraction ϕ or temperature T is high, but opposes self-assembly when both ϕ and T are sufficiently low. This allows us to identify a boundary line in the ϕ-T plane that separates two regions of qualitatively different influence of chain stiffness on self-assembly. The enthalpy and entropy of self-assembly are usually treated as adjustable parameters in classical Flory-Huggins type theories for the equilibrium self-assembly of polymers, but they are demonstrated here to strongly depend on chain stiffness. Moreover, illustrative calculations for the dependence of the entropy density of linear telechelic polymer melts on chain stiffness demonstrate the importance of including semiflexibility within the LCT when exploring the nature of glass formation in models of linear telechelic polymer melts.

Lattice model of linear telechelic polymer melts. I. Inclusion of chain semiflexibility in the lattice cluster theory

The lattice cluster theory (LCT) for the thermodynamics of polymer systems has recently been reformulated to treat strongly interacting self-assembling polymers composed of fully flexible linear telechelic chains [J. Dudowicz and K. F. Freed, J. Chem. Phys. 136, 064902 (2012)]. Here, we further extend the LCT for linear telechelic polymer melts to include a description of chain semiflexibility, which is treated by introducing a bending energy penalty whenever a pair of consecutive bonds from a single chain lies along orthogonal directions. An analytical expression for the Helmholtz free energy is derived for the model of semiflexible linear telechelic polymer melts. The extension provides a theoretical tool for investigating the influence of chain stiffness on the thermodynamics of self-assembling telechelic polymers, and for further exploring the influence of self-assembly on glass formation in such systems.

Lattice cluster theory of associating telechelic polymers. III. Order parameter and average degree of self-assembly, transition temperature, and specific heat

The lattice cluster theory of strongly interacting, structured polymer fluids is applied to determine the thermodynamic properties of solutions of telechelic polymers that may associate through bifunctional end groups. Hence, this model represents a significant albeit natural extension of a diverse array of prior popular equilibrium polymerization models in which structureless "bead" monomers associate into chain-like clusters under equilibrium conditions. In particular, the thermodynamic description of the self-assembly of linear telechelic chains in small molecule solvents (initiated in Paper II) is systematically extended through calculations of the order parameter Φ and average degree <N> of self-assembly, the self-assembly transition temperature T(p), and the specific heat C(V) of solutions of telechelic molecules. Special focus is placed on examining how molecular and thermodynamic parameters, such as the solution composition φ, temperature T, microscopic interaction energies (ε(s) and ε), and length M of individual telechelic chains, influence the computed thermodynamic quantities that are commonly used to characterize self-assembling systems.

Lattice cluster theory of associating polymers. IV. Phase behavior of telechelic polymer solutions

The newly developed lattice cluster theory (in Paper I) for the thermodynamics of solutions of telechelic polymers is used to examine the phase behavior of these complex fluids when effective polymer-solvent interactions are unfavorable. The telechelics are modeled as linear, fully flexible, polymer chains with mono-functional stickers at the two chain ends, and these chains are assumed to self-assemble upon cooling. Phase separation is generated through the interplay of self-assembly and polymer/solvent interactions that leads to an upper critical solution temperature phase separation. The variations of the boundaries for phase stability and the critical temperature and composition are analyzed in detail as functions of the number M of united atom groups in a telechelic chain and the microscopic nearest neighbor interaction energy ε(s) driving the self-assembly. The coupling between self-assembly and unfavorable polymer/solvent interactions produces a wide variety of nontrivial patterns of phase behavior, including an enhancement of miscibility accompanying the increase of the molar mass of the telechelics under certain circumstances. Special attention is devoted to understanding this unusual trend in miscibility.