Targeted modulation of protein liquid-liquid phase separation by evolution of amino-acid sequence.

Rationally and efficiently modifying the amino-acid sequence of proteins to control their ability to undergo liquid-liquid phase separation (LLPS) on demand is not only highly desirable, but can also help to elucidate which protein features are important for LLPS. Here, we propose a computational method that couples a genetic algorithm to a sequence-dependent coarse-grained protein model to evolve the amino-acid sequences of phase-separating intrinsically disordered protein regions (IDRs), and purposely enhance or inhibit their capacity to phase-separate. We validate the predicted critical solution temperatures of the mutated sequences with ABSINTH, a more accurate all-atom model. We apply the algorithm to the phase-separating IDRs of three naturally occurring proteins, namely FUS, hnRNPA1 and LAF1, as prototypes of regions that exist in cells and undergo homotypic LLPS driven by different types of intermolecular interaction, and we find that the evolution of amino-acid sequences towards enhanced LLPS is driven in these three cases, among other factors, by an increase in the average size of the amino acids. However, the direction of change in the molecular driving forces that enhance LLPS (such as hydrophobicity, aromaticity and charge) depends on the initial amino-acid sequence. Finally, we show that the evolution of amino-acid sequences to modulate LLPS is strongly coupled to the make-up of the medium (e.g. the presence or absence of RNA), which may have significant implications for our understanding of phase separation within the many-component mixtures of biological systems

Expansion of Intrinsically Disordered Proteins Increases the Range of Stability of Liquid-Liquid Phase Separation.

Proteins containing intrinsically disordered regions (IDRs) are ubiquitous within biomolecular condensates, which are liquid-like compartments within cells formed through liquid-liquid phase separation (LLPS). The sequence of amino acids of a protein encodes its phase behaviour, not only by establishing the patterning and chemical nature (e.g., hydrophobic, polar, charged) of the various binding sites that facilitate multivalent interactions, but also by dictating the protein conformational dynamics. Besides behaving as random coils, IDRs can exhibit a wide-range of structural behaviours, including conformational switching, where they transition between alternate conformational ensembles. Using Molecular Dynamics simulations of a minimal coarse-grained model for IDRs, we show that the role of protein conformation has a non-trivial effect in the liquid-liquid phase behaviour of IDRs. When an IDR transitions to a conformational ensemble enriched in disordered extended states, LLPS is enhanced. In contrast, IDRs that switch to ensembles that preferentially sample more compact and structured states show inhibited LLPS. This occurs because extended and disordered protein conformations facilitate LLPS-stabilising multivalent protein-protein interactions by reducing steric hindrance; thereby, such conformations maximize the molecular connectivity of the condensed liquid network. Extended protein configurations promote phase separation regardless of whether LLPS is driven by homotypic and/or heterotypic protein-protein interactions. This study sheds light on the link between the dynamic conformational plasticity of IDRs and their liquid-liquid phase behaviour

DNA binds to a specific site of the adhesive blood-protein von Willebrand factor guided by electrostatic interactions.

Neutrophils release their intracellular content, DNA included, into the bloodstream to form neutrophil extracellular traps (NETs) that confine and kill circulating pathogens. The mechanosensitive adhesive blood protein, von Willebrand Factor (vWF), interacts with the extracellular DNA of NETs to potentially immobilize them during inflammatory and coagulatory conditions. Here, we elucidate the previously unknown molecular mechanism governing the DNA-vWF interaction by integrating atomistic, coarse-grained, and Brownian dynamics simulations, with thermophoresis, gel electrophoresis, fluorescence correlation spectroscopy (FCS), and microfluidic experiments. We demonstrate that, independently of its nucleotide sequence, double-stranded DNA binds to a specific helix of the vWF A1 domain, via three arginines. This interaction is attenuated by increasing the ionic strength. Our FCS and microfluidic measurements also highlight the key role shear-stress has in enabling this interaction. Our simulations attribute the previously-observed platelet-recruitment reduction and heparin-size modulation, upon establishment of DNA-vWF interactions, to indirect steric hindrance and partial overlap of the binding sites, respectively. Overall, we suggest electrostatics-guiding DNA to a specific protein binding site-as the main driving force defining DNA-vWF recognition. The molecular picture of a key shear-mediated DNA-protein interaction is provided here and it constitutes the basis for understanding NETs-mediated immune and hemostatic responses

Multiscale Modelling of Biomolecular Phase Behaviour

Elucidating the physicochemical laws that govern the phase behaviour of biomolecules is a cardinal milestone to understand crucial open questions in biology ranging from the spatiotemporal organisation of the cytoplasm in the eukaryotic cell to the emergence of pathological cellular states. Understanding the governing principles of protein and nucleic acid phase behaviour implies deciphering how atomic-scale processes modulate molecular interactions, and how such interactions, in turn, dictate the large-scale thermodynamic behaviour of biomolecules. In this thesis, Molecular Dynamics and Statistical Mechanics have been used to understand the driving forces that modulate protein collective behaviour along such wide-ranging spatial scales. Simulations with a minimal coarse-grained model showed how enriching the structural ensemble of proteins in open, expanded conformations can promote their propensity to undergo liquid-liquid phase separation. This is because more expanded conformations favour inter-protein interactions which increase the molecular connectivity inside protein droplets, favouring their stability. In addition, I present a pioneering chemically specific coarse-grained force field for intrinsically disordered proteins that explicitly considers the effect of solvent in phase separation of proteins. I find our model in near-quantitative experimental agreement both with single-protein structural descriptors as well as with available thermodynamic data of phase separated-protein condensates such as water and protein content in FUS condensates. The model enables studying the structure of the solvent inside protein droplets and recovers experimentally observed salt-inhibited and salt-driven phase separation. Next, I focus on delineating the role of spontaneous beta-sheet folding – prime drivers for neurodisease-linked phase transitions – in the kinetic and transport properties of condensates as a function of the abundance and strength of such transitions. Using atomistic umbrella sampling simulations I estimated the free energy of binding when 8-residue NUP-98 sequence replicas remain disordered as opposed to when they undergo folding into a beta-sheet, finding that folding can increase the binding free energy up to an order of magnitude. Transferring such information from the atomistic to a coarse-grained model, it was established how the timescale between diffusion loss and droplet coalescence determines the shape of phase-separated protein condensates. Finally, I study how disorder-to-order transitions can give rise to single component Fused in Sacroma (FUS) multi-phase condensates. Atomistic free energy calculations were used to establish that such structural transitions can increase four-fold the interaction between FUS sequences which are prone to spontaneous beta-sheet folding. I transferred these findings to a sequence-dependent coarse-grained model which predicted that proteins that underwent folding contribute to lower the interfacial free energy of the condensates. Moreover, our dynamical minimal model informed with the findings from the atomistic and coarse-grained scales predicts single component FUS multi-phase condensates with an inhomogenous organisation where a liquid-core of disordered proteins is surrounded by a gel-like shell of the proteins which underwent disorder-to-order transitions. Taken together, this work sheds light on the molecular driving forces of protein phase behaviour that are crucial to understand cell function both in health and disease and proposes a pioneering approach that advances the realism and explanatory power of coarse-grained biomolecular simulations.EPSRSC grant agreement EP/N509620/

'RNA modulation of transport properties and stability in phase-separated condensates.

One of the key mechanisms employed by cells to control their spatiotemporal organization is the formation and dissolution of phase-separated condensates. The balance between condensate assembly and disassembly can be critically regulated by the presence of RNA. In this work, we use a chemically-accurate sequence-dependent coarse-grained model for proteins and RNA to unravel the impact of RNA in modulating the transport properties and stability of biomolecular condensates. We explore the phase behavior of several RNA-binding proteins such as FUS, hnRNPA1, and TDP-43 proteins along with that of their corresponding prion-like domains and RNA recognition motifs from absence to moderately high RNA concentration. By characterizing the phase diagram, key molecular interactions, surface tension, and transport properties of the condensates, we report a dual RNA-induced behavior: on the one hand, RNA enhances phase separation at low concentration as long as the RNA radius of gyration is comparable to that of the proteins, whereas at high concentration, it inhibits the ability of proteins to self-assemble independently of its length. On the other hand, along with the stability modulation, the viscosity of the condensates can be considerably reduced at high RNA concentration as long as the length of the RNA chains is shorter than that of the proteins. Conversely, long RNA strands increase viscosity even at high concentration, but barely modify protein self-diffusion which mainly depends on RNA concentration and on the effect RNA has on droplet density. On the whole, our work rationalizes the different routes by which RNA can regulate phase separation and condensate dynamics, as well as the subsequent aberrant rigidification implicated in the emergence of various neuropathologies and age-related diseases.EPSR

Kinetic interplay between droplet maturation and coalescence modulates shape of aged protein condensates.

Biomolecular condensates formed by the process of liquid-liquid phase separation (LLPS) play diverse roles inside cells, from spatiotemporal compartmentalisation to speeding up chemical reactions. Upon maturation, the liquid-like properties of condensates, which underpin their functions, are gradually lost, eventually giving rise to solid-like states with potential pathological implications. Enhancement of inter-protein interactions is one of the main mechanisms suggested to trigger the formation of solid-like condensates. To gain a molecular-level understanding of how the accumulation of stronger interactions among proteins inside condensates affect the kinetic and thermodynamic properties of biomolecular condensates, and their shapes over time, we develop a tailored coarse-grained model of proteins that transition from establishing weak to stronger inter-protein interactions inside condensates. Our simulations reveal that the fast accumulation of strongly binding proteins during the nucleation and growth stages of condensate formation results in aspherical solid-like condensates. In contrast, when strong inter-protein interactions appear only after the equilibrium condensate has been formed, or when they accumulate slowly over time with respect to the time needed for droplets to fuse and grow, spherical solid-like droplets emerge. By conducting atomistic potential-of-mean-force simulations of NUP-98 peptides-prone to forming inter-protein [Formula: see text]-sheets-we observe that formation of inter-peptide [Formula: see text]-sheets increases the strength of the interactions consistently with the loss of liquid-like condensate properties we observe at the coarse-grained level. Overall, our work aids in elucidating fundamental molecular, kinetic, and thermodynamic mechanisms linking the rate of change in protein interaction strength to condensate shape and maturation during ageing.Adiran Garaizar is funded by the EPRSC Doctoral Programme Training number EP/N509620/

Nucleosome plasticity is a critical element of chromatin liquid–liquid phase separation and multivalent nucleosome interactions

Abstract: Liquid–liquid phase separation (LLPS) is an important mechanism that helps explain the membraneless compartmentalization of the nucleus. Because chromatin compaction and LLPS are collective phenomena, linking their modulation to the physicochemical features of nucleosomes is challenging. Here, we develop an advanced multiscale chromatin model—integrating atomistic representations, a chemically-specific coarse-grained model, and a minimal model—to resolve individual nucleosomes within sub-Mb chromatin domains and phase-separated systems. To overcome the difficulty of sampling chromatin at high resolution, we devise a transferable enhanced-sampling Debye-length replica-exchange molecular dynamics approach. We find that nucleosome thermal fluctuations become significant at physiological salt concentrations and destabilize the 30-nm fiber. Our simulations show that nucleosome breathing favors stochastic folding of chromatin and promotes LLPS by simultaneously boosting the transient nature and heterogeneity of nucleosome–nucleosome contacts, and the effective nucleosome valency. Our work puts forward the intrinsic plasticity of nucleosomes as a key element in the liquid-like behavior of nucleosomes within chromatin, and the regulation of chromatin LLPS

Alternating one-phase and two-phase crystallization mechanisms in octahedral patchy colloids

Colloidal systems possess unique features to investigate the governing principles behind liquid-to-solid transitions. The phase diagram and crystallization landscape of colloidal particles can be finely tuned by the range, number and angular distribution of attractive interactions between the constituent particles. In this work, we present a computational study of colloidal patchy particles with high-symmetry bonding—six patches displaying octahedral symmetry—that can crystallize into distinct competing ordered phases: a cubic simple (CS) lattice, a body-centered cubic (BCC) phase, and two face-centered cubic (FCC) solids (orientationally ordered and disordered). We investigate the underlying mechanisms by which these competing crystals emerge from a disordered fluid at different pressures. Strikingly, we identify instances where the structure of the crystalline embryo corresponds to the stable solid, while in others it corresponds to a metastable crystal whose nucleation is enabled by its lower interfacial free energy with the liquid. Moreover, we find the exceptional phenomenon that, due to a subtle balance between volumetric enthalpy and interfacial free energy, the CS phase nucleates via crystalline cubic nuclei rather than through spherical clusters as the majority of crystal solids in nature. Finally, by examining growth beyond the nucleation stage, we uncover a series of alternating one-phase and two-phase crystallization mechanisms, depending on whether or not the same phase that nucleates keeps growing. Taken together, we show that an octahedral distribution of attractive sites in colloidal particles results in an extremely rich crystallization landscape where subtle differences in pressure crucially determine the crystallizing polymorph

Nucleosome plasticity is a critical element of chromatin liquid-liquid phase separation and multivalent nucleosome interactions.

Liquid-liquid phase separation (LLPS) is an important mechanism that helps explain the membraneless compartmentalization of the nucleus. Because chromatin compaction and LLPS are collective phenomena, linking their modulation to the physicochemical features of nucleosomes is challenging. Here, we develop an advanced multiscale chromatin model-integrating atomistic representations, a chemically-specific coarse-grained model, and a minimal model-to resolve individual nucleosomes within sub-Mb chromatin domains and phase-separated systems. To overcome the difficulty of sampling chromatin at high resolution, we devise a transferable enhanced-sampling Debye-length replica-exchange molecular dynamics approach. We find that nucleosome thermal fluctuations become significant at physiological salt concentrations and destabilize the 30-nm fiber. Our simulations show that nucleosome breathing favors stochastic folding of chromatin and promotes LLPS by simultaneously boosting the transient nature and heterogeneity of nucleosome-nucleosome contacts, and the effective nucleosome valency. Our work puts forward the intrinsic plasticity of nucleosomes as a key element in the liquid-like behavior of nucleosomes within chromatin, and the regulation of chromatin LLPS

Nucleosome plasticity is a critical element of chromatin liquid-liquid phase separation and multivalent nucleosome interactions.

Liquid-liquid phase separation (LLPS) is an important mechanism that helps explain the membraneless compartmentalization of the nucleus. Because chromatin compaction and LLPS are collective phenomena, linking their modulation to the physicochemical features of nucleosomes is challenging. Here, we develop an advanced multiscale chromatin model-integrating atomistic representations, a chemically-specific coarse-grained model, and a minimal model-to resolve individual nucleosomes within sub-Mb chromatin domains and phase-separated systems. To overcome the difficulty of sampling chromatin at high resolution, we devise a transferable enhanced-sampling Debye-length replica-exchange molecular dynamics approach. We find that nucleosome thermal fluctuations become significant at physiological salt concentrations and destabilize the 30-nm fiber. Our simulations show that nucleosome breathing favors stochastic folding of chromatin and promotes LLPS by simultaneously boosting the transient nature and heterogeneity of nucleosome-nucleosome contacts, and the effective nucleosome valency. Our work puts forward the intrinsic plasticity of nucleosomes as a key element in the liquid-like behavior of nucleosomes within chromatin, and the regulation of chromatin LLPS