LASSI: A lattice model for simulating phase transitions of multivalent proteins

Biomolecular condensates form via phase transitions that combine phase separation or demixing and networking of key protein and RNA molecules. Proteins that drive condensate formation are either linear or branched multivalent proteins where multivalence refers to the presence of multiple protein-protein or protein-nucleic acid interaction domains or motifs within a protein. Recent work has shown that multivalent protein drivers of phase transitions are in fact biological instantiations of associative polymers. Such systems can be characterized by stickers-and-spacers architectures where stickers contribute to system-specific spatial hierarchies of directional interactions and spacers control the concentration-dependent inhomogeneities in densities of stickers around one another. The collective effects of interactions among stickers and spacers lead to the emergence of dense droplet phases wherein the stickers form percolated networks of polymers. To enable the calculation of system-specific phase diagrams of multivalent proteins, we have developed LASSI (LAttice simulations of Sticker and Spacer Interactions), which is an efficient open source computational engine for lattice-based polymer simulations built on the stickers and spacers framework. In LASSI, a specific multivalent protein architecture is mapped into a set of beads on the 3-dimensional lattice space with proper coarse-graining, and specific sticker-sticker interactions are modeled as pairwise anisotropic interactions. For efficient and broad search of the conformational ensemble, LASSI uses Monte Carlo methods, and we optimized the move set so that LASSI can handle both dilute and dense systems. Also, we developed quantitative measures to extract phase boundaries from LASSI simulations, using known and hidden collective parameters. We demonstrate the application of LASSI to two known archetypes of linear and branched multivalent proteins. The simulations recapitulate observations from experiments and importantly, they generate novel quantitative insights that augment what can be gleaned from experiments alone. We conclude with a discussion of the advantages of lattice-based approaches such as LASSI and highlight the types of systems across which this engine can be deployed, either to make predictions or to enable the design of novel condensates. Author Summary Spatial and temporal organization of molecular matter is a defining hallmark of cellular ultrastructure and recent attention has focused intensely on organization afforded by membraneless organelles, which are referred to as biomolecular condensates. These condensates form via phase transitions that combine phase separation and networking of condensate-specific protein and nucleic acid molecules. Several questions remain unanswered regarding the driving forces for condensate formation encoded in the architectures of multivalent proteins, the molecular determinants of material properties of condensates, and the determinants of compositional specificity of condensates. Building on recently recognized analogies between associative polymers and multivalent proteins, we have developed and deployed LASSI, an open source computational engine that enables the calculation of system-specific phase diagrams for multivalent proteins. LASSI relies on a priori identification of stickers and spacers within a multivalent protein and mapping the stickers onto a 3-dimensional lattice. A Monte Carlo engine that incorporates a suite of novel and established move sets enables simulations that track density inhomogeneities and changes to the extent of networking among stickers as a function of protein concentration and interaction strengths. Calculation of distribution functions and other nonconserved order parameters allow us to compute full phase diagrams for multivalent proteins modeled using a stickers-and-spacers representation on simple cubic lattices. These predictions are shown to be system-specific and allow us to rationalize experimental observations while also enabling the design of systems with bespoke phase behavior. LASSI can be deployed to study the phase behavior of multicomponent systems, which allows us to make direct contact with cellular biomolecular condensates that are in fact multicomponent systems.

[1]  Jay R. Unruh,et al.  Quantifying Nucleation In Vivo Reveals the Physical Basis of Prion-like Phase Behavior. , 2018, Molecular cell.

[2]  M. Blinov,et al.  Multivalent Molecules as Modulators of RNA Granule Size and Composition. , 2017, Biophysical journal.

[3]  D. Hernandez-Verdun,et al.  Nucleolus: the fascinating nuclear body , 2007, Histochemistry and Cell Biology.

[4]  R. Pappu,et al.  Phase Separation of Intrinsically Disordered Proteins. , 2018, Methods in enzymology.

[5]  Michael K. Rosen,et al.  Intrinsically disordered sequences enable modulation of protein phase separation through distributed tyrosine motifs , 2017, The Journal of Biological Chemistry.

[6]  R. Parker,et al.  Compositional Control of Phase-Separated Cellular Bodies , 2016, Cell.

[7]  M. M. Telo da Gama,et al.  Dynamics of network fluids. , 2017, Advances in colloid and interface science.

[8]  Gregory A Voth,et al.  A multiscale coarse-graining method for biomolecular systems. , 2005, The journal of physical chemistry. B.

[9]  Hal Daumé,et al.  A Bayesian statistics approach to multiscale coarse graining. , 2008, The Journal of chemical physics.

[10]  G. Voth The Theory of Ultra Coarse-graining , 2013 .

[11]  Hajime Tanaka,et al.  Volume-shrinking kinetics of transient gels as a consequence of dynamic interplay between phase separation and mechanical relaxation , 2018, Physical Review E.

[12]  Hue Sun Chan,et al.  Charge pattern matching as a ‘fuzzy’ mode of molecular recognition for the functional phase separations of intrinsically disordered proteins , 2017, 1707.08990.

[13]  Raoul Kopelman,et al.  Fractal Reaction Kinetics , 1988, Science.

[14]  Erin M. Langdon,et al.  RNA Controls PolyQ Protein Phase Transitions. , 2015, Molecular cell.

[15]  Diana M. Mitrea,et al.  Compositional adaptability in NPM1-SURF6 scaffolding networks enabled by dynamic switching of phase separation mechanisms , 2018, Nature Communications.

[16]  Paul J. Flory,et al.  Molecular Size Distribution in Three Dimensional Polymers. I. Gelation1 , 1941 .

[17]  R. Pappu,et al.  Profilin reduces aggregation and phase separation of huntingtin N-terminal fragments by preferentially binding to soluble monomers and oligomers , 2018, The Journal of Biological Chemistry.

[18]  B. Widom,et al.  Some Topics in the Theory of Fluids , 1963 .

[19]  T. C. B. McLeish,et al.  Polymer Physics , 2009, Encyclopedia of Complexity and Systems Science.

[20]  Daan Frenkel,et al.  Configurational bias Monte Carlo: a new sampling scheme for flexible chains , 1992 .

[21]  Josef Mikeš,et al.  Simulation of polymer network formation by the Monte Carlo method , 1982 .

[22]  A. Hyman,et al.  Liquid-liquid phase separation in biology. , 2014, Annual review of cell and developmental biology.

[23]  Lucas Pelkmans,et al.  Kinase-controlled phase transition of membraneless organelles in mitosis , 2018, Nature.

[24]  Najeeb M. Halabi,et al.  Protein Sectors: Evolutionary Units of Three-Dimensional Structure , 2009, Cell.

[25]  C. Brangwynne,et al.  The disordered P granule protein LAF-1 drives phase separation into droplets with tunable viscosity and dynamics , 2015, Proceedings of the National Academy of Sciences.

[26]  V. Nguemaha,et al.  Liquid-Liquid Phase Separation of Patchy Particles Illuminates Diverse Effects of Regulatory Components on Protein Droplet Formation , 2018, bioRxiv.

[27]  Ming-Tzo Wei,et al.  Phase behaviour of disordered proteins underlying low density and high permeability of liquid organelles. , 2017, Nature chemistry.

[28]  J. Aulbach Spatiotemporal control of light in turbid media , 2013 .

[29]  G. Fredrickson,et al.  A multi-species exchange model for fully fluctuating polymer field theory simulations. , 2014, The Journal of chemical physics.

[30]  G. Voth,et al.  Quantum theory of multiscale coarse-graining. , 2018, The Journal of chemical physics.

[31]  R. Pappu,et al.  Spontaneous driving forces give rise to protein−RNA condensates with coexisting phases and complex material properties , 2018, Proceedings of the National Academy of Sciences.

[32]  D. Resasco,et al.  Virtual Cell: computational tools for modeling in cell biology , 2012, Wiley interdisciplinary reviews. Systems biology and medicine.

[33]  D. Spector,et al.  Compartmentalization of RNA processing factors within nuclear speckles. , 2000, Journal of structural biology.

[34]  Diana M. Mitrea,et al.  Structural polymorphism in the N-terminal oligomerization domain of NPM1 , 2014, Proceedings of the National Academy of Sciences.

[35]  Christopher B. Stanley,et al.  Self-interaction of NPM1 modulates multiple mechanisms of liquid–liquid phase separation , 2018, Nature Communications.

[36]  C. Sing,et al.  Mapping the phase behavior of coacervate-driven self-assembly in diblock copolyelectrolytes. , 2019, Soft matter.

[37]  Paul J. Flory Constitution of Three-Dimensional Polymers and the Theory of Gelation , 1942 .

[38]  Jared E. Toettcher,et al.  Spatiotemporal Control of Intracellular Phase Transitions Using Light-Activated optoDroplets , 2017, Cell.

[39]  Wenwei Zheng,et al.  Sequence determinants of protein phase behavior from a coarse-grained model , 2017, bioRxiv.

[40]  Kurt Kremer,et al.  The bond fluctuation method: a new effective algorithm for the dynamics of polymers in all spatial dimensions , 1988 .

[41]  D. Theodorou Equilibration and coarse-graining methods for polymers , 2006 .

[42]  H. Chan,et al.  Coarse-grained residue-based models of disordered protein condensates: utility and limitations of simple charge pattern parameters. , 2018, Physical chemistry chemical physics : PCCP.

[43]  R. Pappu,et al.  A Molecular Grammar Governing the Driving Forces for Phase Separation of Prion-like RNA Binding Proteins , 2018, Cell.

[44]  H. Chan,et al.  A Lattice Model of Charge-Pattern-Dependent Polyampholyte Phase Separation. , 2018, The journal of physical chemistry. B.

[45]  M. Rosen,et al.  Stoichiometry controls activity of phase-separated clusters of actin signaling proteins , 2019, Science.

[46]  Ronald D. Vale,et al.  Phase separation of signaling molecules promotes T cell receptor signal transduction , 2016, Science.

[47]  Mingjie Zhang,et al.  Reconstituted Postsynaptic Density as a Molecular Platform for Understanding Synapse Formation and Plasticity , 2018, Cell.

[48]  Michael Rubinstein,et al.  Thermoreversible Gelation in Solutions of Associative Polymers. 1. Statics , 1998 .

[49]  A. Dobrynin,et al.  Solutions of associative polymers , 1997 .

[50]  R. Pappu,et al.  Injectable tissue integrating networks from recombinant polypeptides with tunable order , 2018, Nature Materials.

[51]  G. Fredrickson,et al.  Complete Phase Diagram for Liquid-Liquid Phase Separation of Intrinsically Disordered Proteins. , 2019, The journal of physical chemistry letters.

[52]  C. Brangwynne,et al.  Liquid phase condensation in cell physiology and disease , 2017, Science.

[53]  Paul S. Russo,et al.  Phase Transitions in the Assembly of MultiValent Signaling Proteins , 2016 .

[54]  W. Nowicki,et al.  Monte Carlo simulations of a polymer chain conformation. The effectiveness of local moves algorithms and estimation of entropy , 2013, Journal of Molecular Modeling.

[55]  G. Fredrickson,et al.  Recent Developments in Fully Fluctuating Field-Theoretic Simulations of Polymer Melts and Solutions. , 2016, The journal of physical chemistry. B.

[56]  S. Shanbhag,et al.  On the relationship between two popular lattice models for polymer melts. , 2008, The Journal of chemical physics.

[57]  G. Voth,et al.  Coarse-graining in interaction space: a systematic approach for replacing long-range electrostatics with short-range potentials. , 2008, The journal of physical chemistry. B.

[58]  R. Pappu,et al.  Phase separation of a yeast prion protein promotes cellular fitness , 2018, Science.

[59]  A. W. Rosenbluth,et al.  MONTE CARLO CALCULATION OF THE AVERAGE EXTENSION OF MOLECULAR CHAINS , 1955 .

[60]  C. Sing,et al.  Transfer matrix theory of polymer complex coacervation. , 2017, Soft matter.

[61]  K. Binder,et al.  Monte Carlo simulations of polymer dynamics: Recent advances , 1997 .

[62]  Nicolas L. Fawzi,et al.  Phosphorylation of the FUS low‐complexity domain disrupts phase separation, aggregation, and toxicity , 2017, The EMBO journal.

[63]  R. Pappu,et al.  Functional Implications of Intracellular Phase Transitions. , 2018, Biochemistry.

[64]  Marco Y. Hein,et al.  A Liquid-to-Solid Phase Transition of the ALS Protein FUS Accelerated by Disease Mutation , 2015, Cell.

[65]  J. Shaffer Effects of chain topology on polymer dynamics: Bulk melts , 1994 .

[66]  Rohit V Pappu,et al.  CAMELOT: A machine learning approach for coarse-grained simulations of aggregation of block-copolymeric protein sequences. , 2015, The Journal of chemical physics.

[67]  Walter H. Stockmayer,et al.  Theory of Molecular Size Distribution and Gel Formation in Branched‐Chain Polymers , 1943 .

[68]  Diana M. Mitrea,et al.  Coexisting Liquid Phases Underlie Nucleolar Subcompartments , 2016, Cell.

[69]  Christopher B. Stanley,et al.  Nucleophosmin integrates within the nucleolus via multi-modal interactions with proteins displaying R-rich linear motifs and rRNA , 2016, eLife.

[70]  David R. Liu,et al.  Sequence Determinants of Intracellular Phase Separation by Complex Coacervation of a Disordered Protein. , 2016, Molecular cell.

[71]  Jared E. Toettcher,et al.  Mapping Local and Global Liquid Phase Behavior in Living Cells Using Photo-Oligomerizable Seeds , 2019, Cell.

[72]  Benoît Roux,et al.  The Theory of Ultra-Coarse-Graining. 1. General Principles. , 2013, Journal of chemical theory and computation.

[73]  A. Hyman,et al.  The Centrosome Is a Selective Condensate that Nucleates Microtubules by Concentrating Tubulin , 2017, Cell.

[74]  Kurt Kremer,et al.  Monte Carlo simulation of lattice models for macromolecules , 1988 .

[75]  Reentrant Phase Transition Drives Dynamic Substructure Formation in Ribonucleoprotein Droplets. , 2017, Angewandte Chemie.

[76]  C. Sing,et al.  Interfacial properties of polymeric complex coacervates from simulation and theory. , 2018, The Journal of chemical physics.

[77]  P. Tomançak,et al.  RNA buffers the phase separation behavior of prion-like RNA binding proteins , 2018, Science.

[78]  S. Michnick,et al.  Probing inhomogeneous diffusion in the microenvironments of phase-separated polymers under confinement , 2018, bioRxiv.

[79]  Ying Hu,et al.  Dynamic Monte Carlo Simulation of Polymers: Cooperative Move Algorithm , 2003 .

[80]  Anthony A. Hyman,et al.  Polar Positioning of Phase-Separated Liquid Compartments in Cells Regulated by an mRNA Competition Mechanism , 2016, Cell.

[81]  Roberto Mauri,et al.  Diffusion-driven phase separation of deeply quenched mixtures , 1998 .

[82]  R. Pappu,et al.  Quantitative analysis of multilayer organization of proteins and RNA in nuclear speckles at super resolution , 2017, Journal of Cell Science.

[83]  Anthony A. Hyman,et al.  Biomolecular condensates: organizers of cellular biochemistry , 2017, Nature Reviews Molecular Cell Biology.

[84]  R. Pappu,et al.  To Mix, or To Demix, That Is the Question. , 2017, Biophysical journal.

[85]  R. Pappu,et al.  Differential solvation of intrinsically disordered linkers drives the formation of spatially organized droplets in ternary systems of linear multivalent proteins , 2018 .

[86]  Luis Kuhn Cuellar,et al.  The Eukaryotic CO2-Concentrating Organelle Is Liquid-like and Exhibits Dynamic Reorganization , 2017, Cell.

[87]  D. Frenkel,et al.  Numerical study of the phase diagram of a mixture of spherical and rodlike colloids , 1994 .

[88]  Huan‐Xiang Zhou,et al.  Fast Method for Computing Chemical Potentials and Liquid-Liquid Phase Equilibria of Macromolecular Solutions. , 2016, The journal of physical chemistry. B.

[89]  J. Wyman,et al.  Ligand-linked phase changes in a biological system: applications to sickle cell hemoglobin. , 1980, Proceedings of the National Academy of Sciences of the United States of America.

[90]  A. Hyman,et al.  An aberrant phase transition of stress granules triggered by misfolded protein and prevented by chaperone function , 2017, The EMBO journal.

[91]  Hue Sun Chan,et al.  Theories for Sequence-Dependent Phase Behaviors of Biomolecular Condensates. , 2018, Biochemistry.

[92]  Peter Tompa,et al.  Polymer physics of intracellular phase transitions , 2015, Nature Physics.

[93]  A. Semenov,et al.  Thermoreversible Gelation in Solutions of Associating Polymers. 2. Linear Dynamics , 1998 .

[94]  A. Kanagaraj,et al.  Phase Separation by Low Complexity Domains Promotes Stress Granule Assembly and Drives Pathological Fibrillization , 2015, Cell.

[95]  Lanyuan Lu,et al.  Fitting coarse-grained distribution functions through an iterative force-matching method. , 2013, The Journal of chemical physics.

[96]  Amy S. Gladfelter,et al.  mRNA structure determines specificity of a polyQ-driven phase separation , 2017, Science.

[97]  A. Hyman,et al.  Germline P Granules Are Liquid Droplets That Localize by Controlled Dissolution/Condensation , 2009, Science.

[98]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[99]  Wenwei Zheng,et al.  Relation between single-molecule properties and phase behavior of intrinsically disordered proteins , 2018, Proceedings of the National Academy of Sciences.

[100]  H. Chan,et al.  Structural and hydrodynamic properties of an intrinsically disordered region of a germ cell-specific protein on phase separation , 2017, Proceedings of the National Academy of Sciences.

[101]  Y. Sheng,et al.  Statics and dynamics of a single polymer chain confined in a tube , 2001 .

[102]  H. Chan,et al.  Sequence-Specific Polyampholyte Phase Separation in Membraneless Organelles. , 2016, Physical review letters.

[103]  C. Sing,et al.  Tuning chain interaction entropy in complex coacervation using polymer stiffness, architecture, and salt valency , 2018 .

[104]  Hong Lin,et al.  Pi-Pi contacts are an overlooked protein feature relevant to phase separation , 2018, eLife.

[105]  Monika Fuxreiter,et al.  The Structure and Dynamics of Higher-Order Assemblies: Amyloids, Signalosomes, and Granules , 2016, Cell.

[106]  P. Gennes Reptation of a Polymer Chain in the Presence of Fixed Obstacles , 1971 .

[107]  David L. Spector,et al.  Nuclear speckles: a model for nuclear organelles , 2003, Nature Reviews Molecular Cell Biology.

[108]  Sudeep Banjade,et al.  Phase transitions of multivalent proteins can promote clustering of membrane receptors , 2014, eLife.

[109]  Timothy D. Craggs,et al.  Phase Transition of a Disordered Nuage Protein Generates Environmentally Responsive Membraneless Organelles , 2015, Molecular cell.

[110]  D. Durand,et al.  Monte Carlo simulation of particle aggregation and gelation: I. Growth, structure and size distribution of the clusters , 2004, The European physical journal. E, Soft matter.

[111]  C. Sing,et al.  Designing Electrostatic Interactions via Polyelectrolyte Monomer Sequence , 2019, ACS central science.