Small ion effects on self-coacervation phenomena in block polyampholytes - PubMed
- ️Tue Jan 01 2019
Small ion effects on self-coacervation phenomena in block polyampholytes
Scott P O Danielsen et al. J Chem Phys. 2019.
Abstract
Self-coacervation is a phenomenon in which a solution of polyampholytes spontaneously phase separates into a dense liquid coacervate phase, rich in the polyampholyte, coexisting with a dilute supernatant phase. Such coacervation results in the formation of membraneless organelles in vivo and has further been applied industrially as synthetic encapsulants and coatings. It has been suggested that coacervation is primarily driven by the entropy gain from releasing counter-ions upon complexation. Using fully fluctuating field-theoretic simulations employing complex Langevin sampling and complementary molecular dynamics simulations, we have determined that the small ions contribute only weakly to the self-coacervation behavior of charge-symmetric block polyampholytes in solution. Salt partitioning between the supernatant and coacervate is also found to be negligible in the weak-binding regime at low electrostatic strengths. Asymmetries in charge distribution along the polyampholytes can cause net-charges that lead to "tadpole" configurations in dilute solution and the suppression of phase separation at low salt content. The field and particle-based simulation results are compared with analytical predictions from the random phase approximation (RPA) and postulated scaling relationships. The qualitative trends are mostly captured by the RPA, but the approximation fails at low concentration.
Figures

Schematic of the Gibbs ensemble. The volume is partitioned into two simulation cells corresponding to a dense and a dilute phase. Particle concentrations are exchanged to achieve chemical equilibrium and phase volumes are exchanged to achieve mechanical equilibrium.

Convergence of Gibbs ensemble calculations for the phase coexistence binodal curves (lB/b = 0.13, v/b3 = 0.0068, a/b = 0.8) to prior results (red dashed lines) for the diblock polyampholyte at different numbers of FTS-CL time steps per Gibbs step (symbols and connecting solid lines). Inset: expansion of the high density region to show the narrowing of error bars at different block averaging conditions.

Phase diagram showing the relatively weak effects of including explicit counter-ions on the phase diagram of a diblock polyampholyte (A50C50). Symbols are binodal points obtained from FTS-CL and lines are RPA predictions. The Gaussian smearing width is set at a/b = 0.8 and the excluded-volume parameter at v/b3 = 0.0068.

Total-density structure factor, Sn(k), and charge-weighted structure factor, Se(k), of the diblock polyampholyte (A50C50) in the dilute supernatant and coacervate as a function of the electrostatic strength, lB/b, with and without explicit counter-ions. The supernatant is at a polymer number density of ρpb3 = 7.35 × 10−7, and the coacervate is at ρpb3 = 0.0735. Symbols are obtained from FTS-CL, and lines are RPA predictions. The Gaussian smearing width is set at a/b = 0.8 and the excluded-volume parameter is set at v/b3 = 0.0068.

Phase diagram showing the modulation of critical electrostatic strength with inclusion of excess salt on the phase diagram of a diblock polyampholyte (A50C50). Symbols are binodal points obtained from FTS-CL, and lines are RPA predictions. The Gaussian smearing width is set at a/b = 0.8 and the excluded-volume parameter is set at v/b3 = 0.0068. The different salt conditions are presented as number fractions ϕsalt,±; the concentration (number density) of both polymer and salt increases along the horizontal axis.

Total-density structure factor, Sn(k), and charge-weighted structure factor, Se(k), of the diblock polyampholyte (A50C50) in the dilute supernatant and coacervate as a function of the number fraction of excess salt, ϕsalt,±. Symbols are obtained from FTS-CL, and lines are RPA predictions. The supernatant is at a polymer number density of ρpb3 = 7.35 × 10−7, and the coacervate is at ρpb3 = 0.0735. The Gaussian smearing width is set at a/b = 0.8, the excluded-volume parameter is set at v/b3 = 0.0068, and the electrostatic strength is set at lB/b = 0.065.

Phase diagram showing excess salt shifting the location of the two-phase region for the diblock polyampholyte (A50C50). The shaded region is the region between the binodal curves. Representative tie lines are shown as horizontal solid lines. Dotted (ϕsalt,± = 0.0) and dashed (ϕsalt,± = 0.3) diagonal lines show stoichiometry constraints (i.e., concentration of small ions corresponding to the polymer concentration for each salt condition.) The Gaussian smearing width is set at a/b = 0.8, the excluded-volume parameter is set at v/b3 = 0.0068, and the electrostatic strength is set at lB/b = 0.065.

Phase diagram showing charge asymmetry shrinking phase coexistence and shifting binodal boundaries to lower polyampholyte and higher salt concentrations. The shaded region is the region between the binodal curves. Representative tie lines are shown as solid lines. The Gaussian smearing width is set at a/b = 0.8, the excluded-volume parameter is set at v/b3 = 0.0068, and the electrostatic strength is set at lB/b = 0.065.

Total-density structure factor, Sn(k), and charge-weighted structure factor, Se(k), of the diblock polyampholyte in the dilute supernatant and coacervate as a function of charge asymmetry. Symbols are obtained from FTS-CL, and lines are RPA predictions. The supernatant is at a polymer number density of ρpb3 = 7.35 × 10−7, and the coacervate is at ρpb3 = 0.0735. The Gaussian smearing width is set at a/b = 0.8, the excluded-volume parameter is set at v/b3 = 0.0068, and the electrostatic strength is set at lB/b = 0.065.

Single chain structure factor, P(k), in the dilute supernatant (ρpb3 = 7.35 × 10−7) as a function of charge asymmetry. The Gaussian smearing width is set at a/b = 0.8, the excluded-volume parameter is set at v/b3 = 0.0068, and the electrostatic strength is set at lB/b = 0.065. Insets: Snapshots from the MD simulations showing globular chain conformations at zero net-charge and the emergence of “tadpole” conformations with increasing charge asymmetry.

Phase diagram showing multivalent counter-ions shrinking the region for phase coexistence. The shaded region is the region between the binodal curves. Representative tie lines are shown as solid lines. The Gaussian smearing width is set at a/b = 0.8, the excluded-volume parameter is set at v/b3 = 0.0068, and the electrostatic strength is set at lB/b = 0.065.

Total-density structure factor, Sn(k), and charge-weighted structure factor, Se(k), of the diblock polyampholyte (A50C50) in the dilute supernatant and coacervate as a function of counter-cation valency. Symbols are obtained from FTS-CL, and lines are RPA predictions. The supernatant is at a polymer number density of ρpb3 = 7.35 × 10−7, and the coacervate is at ρpb3 = 0.0735. The Gaussian smearing width is set at a/b = 0.8, the excluded-volume parameter is set at v/b3 = 0.0068, and the electrostatic strength is set at lB/b = 0.065.

Phase diagram showing the importance of charge neutrality in dictating partitioning and phase coexistence of a charge asymmetric diblock polyampholyte (A48C52) with multivalent counter-ions. The shaded region is the region between the binodal curves. Representative tie lines are shown as solid lines. The Gaussian smearing width is set at a/b = 0.8, the excluded-volume parameter is set at v/b3 = 0.0068, and the electrostatic strength is set at lB/b = 0.065.
Similar articles
-
Molecular design of self-coacervation phenomena in block polyampholytes.
Danielsen SPO, McCarty J, Shea JE, Delaney KT, Fredrickson GH. Danielsen SPO, et al. Proc Natl Acad Sci U S A. 2019 Apr 23;116(17):8224-8232. doi: 10.1073/pnas.1900435116. Epub 2019 Apr 4. Proc Natl Acad Sci U S A. 2019. PMID: 30948640 Free PMC article.
-
Theory of polyelectrolyte complexation-Complex coacervates are self-coacervates.
Delaney KT, Fredrickson GH. Delaney KT, et al. J Chem Phys. 2017 Jun 14;146(22):224902. doi: 10.1063/1.4985568. J Chem Phys. 2017. PMID: 29166038
-
Transfer Matrix Model of pH Effects in Polymeric Complex Coacervation.
Knoerdel AR, Blocher McTigue WC, Sing CE. Knoerdel AR, et al. J Phys Chem B. 2021 Aug 12;125(31):8965-8980. doi: 10.1021/acs.jpcb.1c03065. Epub 2021 Jul 30. J Phys Chem B. 2021. PMID: 34328340
-
Veis A. Veis A. Adv Colloid Interface Sci. 2011 Sep 14;167(1-2):2-11. doi: 10.1016/j.cis.2011.01.007. Epub 2011 Mar 5. Adv Colloid Interface Sci. 2011. PMID: 21377640 Free PMC article. Review.
-
Development of the modern theory of polymeric complex coacervation.
Sing CE. Sing CE. Adv Colloid Interface Sci. 2017 Jan;239:2-16. doi: 10.1016/j.cis.2016.04.004. Epub 2016 Apr 29. Adv Colloid Interface Sci. 2017. PMID: 27161661 Review.
Cited by
-
Das S, Lin YH, Vernon RM, Forman-Kay JD, Chan HS. Das S, et al. Proc Natl Acad Sci U S A. 2020 Nov 17;117(46):28795-28805. doi: 10.1073/pnas.2008122117. Epub 2020 Nov 2. Proc Natl Acad Sci U S A. 2020. PMID: 33139563 Free PMC article.
-
Huihui J, Ghosh K. Huihui J, et al. J Chem Phys. 2020 Apr 30;152(16):161102. doi: 10.1063/5.0004619. J Chem Phys. 2020. PMID: 32357776 Free PMC article.
-
Lin YH, Brady JP, Chan HS, Ghosh K. Lin YH, et al. J Chem Phys. 2020 Jan 31;152(4):045102. doi: 10.1063/1.5139661. J Chem Phys. 2020. PMID: 32007034 Free PMC article.
-
Ion Pairing and the Structure of Gel Coacervates.
Danielsen SPO, Panyukov S, Rubinstein M. Danielsen SPO, et al. Macromolecules. 2020 Nov 10;53(21):9420-9442. doi: 10.1021/acs.macromol.0c01360. Epub 2020 Oct 30. Macromolecules. 2020. PMID: 34366486 Free PMC article.
-
Complete Phase Diagram for Liquid-Liquid Phase Separation of Intrinsically Disordered Proteins.
McCarty J, Delaney KT, Danielsen SPO, Fredrickson GH, Shea JE. McCarty J, et al. J Phys Chem Lett. 2019 Apr 18;10(8):1644-1652. doi: 10.1021/acs.jpclett.9b00099. Epub 2019 Mar 27. J Phys Chem Lett. 2019. PMID: 30873835 Free PMC article.
References
-
- Kato M., Han T. N. W., Xie S. H., Shi K., Du X. L., Wu L. C., Mirzaei H., Goldsmith E. J., Longgood J., Pei J. M., Grishin N. V., Frantz D. E., Schneider J. W., Chen S., Li L., Sawaya M. R., Eisenberg D., Tycko R., and McKnight S. L., Cell 149, 753 (2012).10.1016/j.cell.2012.04.017 - DOI - PMC - PubMed