Study of the solid-liquid-vapour phase equilibria of flexible chain molecules using Wertheim's thermodynamic perturbation theory

RNA length has a non-trivial effect in the stability of biomolecular condensates formed by RNA-binding proteins

Biomolecular condensates formed via liquid-liquid phase separation (LLPS) play a crucial role in the spatiotemporal organization of the cell material. Nucleic acids can act as critical modulators in the stability of these protein condensates. To unveil the role of RNA length in regulating the stability of RNA binding protein (RBP) condensates, we present a multiscale computational strategy that exploits the advantages of a sequence-dependent coarse-grained representation of proteins and a minimal coarse-grained model wherein proteins are described as patchy colloids. We find that for a constant nucleotide/protein ratio, the protein fused in sarcoma (FUS), which can phase separate on its own-i.e., via homotypic interactions-only exhibits a mild dependency on the RNA strand length. In contrast, the 25-repeat proline-arginine peptide (PR25), which does not undergo LLPS on its own at physiological conditions but instead exhibits complex coacervation with RNA-i.e., via heterotypic interactions-shows a strong dependence on the length of the RNA strands. Our minimal patchy particle simulations suggest that the strikingly different effect of RNA length on homotypic LLPS versus RBP-RNA complex coacervation is general. Phase separation is RNA-length dependent whenever the relative contribution of heterotypic interactions sustaining LLPS is comparable or higher than those stemming from protein homotypic interactions. Taken together, our results contribute to illuminate the intricate physicochemical mechanisms that influence the stability of RBP condensates through RNA inclusion.

'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.

'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.

Salt dependent phase behavior of intrinsically disordered proteins from a coarse-grained model with explicit water and ions

Multivalent proteins and nucleic acids can self-assemble into biomolecular condensates that contribute to compartmentalize the cell interior. Computer simulations offer a unique view to elucidate the mechanisms and key intermolecular interactions behind the dynamic formation and dissolution of these condensates. In this work, we present a novel approach to include explicit water and salt in sequence-dependent coarse-grained (CG) models for proteins and RNA, enabling the study of biomolecular condensate formation in a salt-dependent manner. Our framework combines a reparameterized version of the HPS protein force field with the monoatomic mW water model and the mW-ion potential for NaCl. We show how our CG model qualitatively captures the experimental radius of the gyration trend of a subset of intrinsically disordered proteins and reproduces the experimental protein concentration and water percentage of the human fused in sarcoma (FUS) low-complexity-domain droplets at physiological salt concentration. Moreover, we perform seeding simulations as a function of salt concentration for two antagonist systems: the engineered peptide PR₂₅ and poly-uridine/poly-arginine mixtures, finding good agreement with their reported in vitro phase behavior with salt concentration in both cases. Taken together, our work represents a step forward towards extending sequence-dependent CG models to include water and salt, and to consider their key role in biomolecular condensate self-assembly.

Valency and Binding Affinity Variations Can Regulate the Multilayered Organization of Protein Condensates with Many Components

Biomolecular condensates, which assemble via the process of liquid-liquid phase separation (LLPS), are multicomponent compartments found ubiquitously inside cells. Experiments and simulations have shown that biomolecular condensates with many components can exhibit multilayered organizations. Using a minimal coarse-grained model for interacting multivalent proteins, we investigate the thermodynamic parameters governing the formation of multilayered condensates through changes in protein valency and binding affinity. We focus on multicomponent condensates formed by scaffold proteins (high-valency proteins that can phase separate on their own via homotypic interactions) and clients (proteins recruited to condensates via heterotypic scaffold-client interactions). We demonstrate that higher valency species are sequestered to the center of the multicomponent condensates, while lower valency proteins cluster towards the condensate interface. Such multilayered condensate architecture maximizes the density of LLPS-stabilizing molecular interactions, while simultaneously reducing the surface tension of the condensates. In addition, multilayered condensates exhibit rapid exchanges of low valency proteins in and out, while keeping higher valency proteins-the key biomolecules involved in condensate nucleation-mostly within. We also demonstrate how modulating the binding affinities among the different proteins in a multicomponent condensate can significantly transform its multilayered structure, and even trigger fission of a condensate into multiple droplets with different compositions.

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.

A perturbed-chain equation of state based on Wertheim TPT for the fully flexible LJ chains in the fluid and solid phases

In the framework of thermodynamic perturbation theory (TPT), a new perturbed-chain equation of state (EOS) is presented for a fully flexible Lennard-Jones (LJ) chain system. The EOS is the sum of repulsive and perturbation contributions. The reference term of the EOS is derived based on first- and second-order TPT of Wertheim for the chains interacting with each other through the Weeks-Chandler-Anderson potential model. In order to derive the perturbation term, we have used the radial distribution function of the hard-chain system with a chain range of m = 2-10 and packing fraction range of η = 0.10-0.72, which cover the entire density range from vapor to solid phases. The performance of the EOS is tested against simulation data of the compressibility factor, residual internal energy, and phase equilibrium. A close agreement was observed across all cases. The EOS has three pure component parameters and is able to describe the global vapor-liquid-solid phase diagram of the LJ chain.

Equations of state for the fully flexible WCA chains in the fluid and solid phases based on Wertheims-TPT2

Based on Wertheim's second order thermodynamic perturbation theory (TPT2), equations of state (EOSs) are presented for the fluid and solid phases of tangent, freely jointed spheres. It is considered that the spheres interact with each other through the Weeks-Chandler-Anderson (WCA) potential. The developed TPT2 EOS is the sum of a monomeric reference term and a perturbation contribution due to bonding. MC NVT simulations are performed to determine the structural properties of the reference system in the reduced temperature range of 0.6 ≤ T* ≤ 4.0 and the packing fraction range of 0.1 ≤ η ≤ 0.72. Mathematical functions are fitted to the simulation results of the reference system and employed in the framework of Wertheim's theory to develop TPT2 EOSs for the fluid and solid phases. The extended EOSs are compared to the MC NPT simulation results of the compressibility factor and internal energy of the fully flexible chain systems. Simulations are performed for the WCA chain system for chain lengths of up to 15 at T* = 1.0, 1.5, 2.0, 3.0. Across all the reduced temperatures, the agreement between the results of the TPT2 EOS and MC simulations is remarkable. Overall Average Absolute Relative Percent Deviation at T* = 1.0 for the compressibility factor in the entire chain lengths we covered is 0.51 and 0.77 for the solid and fluid phases, respectively. Similar features are observed in the case of residual internal energy.

Structure and phase behavior of a confined nanodroplet composed of the flexible chain molecules

A polymer density functional theory has been employed for investigating the structure and phase behaviors of the chain polymer, which is modelled as the tangentially connected sphere chain with an attractive interaction, inside the nanosized pores. The excess free energy of the chain polymer has been approximated as the modified fundamental measure-theory for the hard spheres, the Wertheim's first-order perturbation for the chain connectivity, and the mean-field approximation for the van der Waals contribution. For the value of the chemical potential corresponding to a stable liquid phase in the bulk system and a metastable vapor phase, the flexible chain molecules undergo the liquid-vapor transition as the pore size is reduced; the vapor is the stable phase at small volume, whereas the liquid is the stable phase at large volume. The wide liquid-vapor coexistence curve, which explains the wide range of metastable liquid-vapor states, is observed at low temperature. The increase of temperature and decrease of pore size result in a narrowing of liquid-vapor coexistence curves. The increase of chain length leads to a shift of the liquid-vapor coexistence curve towards lower values of chemical potential. The coexistence curves for the confined phase diagram are contained within the corresponding bulk liquid-vapor coexistence curve. The equilibrium capillary phase transition occurs at a higher chemical potential than in the bulk phase.

Shear viscosity of the Lennard-Jones chain fluid in its gaseous, supercritical, and liquid states

Extensive nonequilibrium molecular dynamics (NEMD) simulations have been carried out in order to estimate the Newtonian shear viscosity of a fluid composed of short chains (up to 16 segments) of jointed spheres [Lennard-Jones chain (LJC)] over a large range of thermodynamic conditions. Using the NEMD results, it is shown that the zero-density contribution decreases with the chain length for a given temperature and is simply proportional to N-1/2 , where N is the number of spheres composing the chain. In addition, it has been noticed that the residual shear viscosity is proportional to the chain length. Then, using these results, a relation is proposed to correlate the shear viscosity of the LJC fluid using the LJ fluid (the monomer) as a reference. It is shown that this correlation is able to provide an excellent estimation of the LJC fluid viscosity compared to NEMD results for N<or=16 over the domain 0<or=rho <or=1.1 and 0.7<or=T <or=6 . Finally, it is shown that the LJC model is unambiguously more efficient than a simple LJ approximation when applied to estimate the shear viscosity of n -butane, if only the sphere or segment diameter is used as an adjustable parameter in both models.

Computer simulation study of the global phase behavior of linear rigid Lennard-Jones chain molecules: Comparison with flexible models

The global phase behavior (i.e., vapor-liquid and fluid-solid equilibria) of rigid linear Lennard-Jones (LJ) chain molecules is studied. The phase diagrams for three-center and five-center rigid model molecules are obtained by computer simulation. The segment-segment bond lengths are L = sigma, so that models of tangent monomers are considered in this study. The vapor-liquid equilibrium conditions are obtained using the Gibbs ensemble Monte Carlo method and by performing isobaric-isothermal NPT calculations at zero pressure. The phase envelopes and critical conditions are compared with those of flexible LJ molecules of tangent segments. An increase in the critical temperature of linear rigid chains with respect to their flexible counterparts is observed. In the limit of infinitely long chains the critical temperature of linear rigid LJ chains of tangent segments seems to be higher than that of flexible LJ chains. The solid-fluid equilibrium is obtained by Gibbs-Duhem integration, and by performing NPT simulations at zero pressure. A stabilization of the solid phase, an increase in the triple-point temperature, and a widening of the transition region are observed for linear rigid chains when compared to flexible chains with the same number of segments. The triple-point temperature of linear rigid LJ chains increases dramatically with chain length. The results of this work suggest that the fluid-vapor transition could be metastable with respect to the fluid-solid transition for chains with more than six LJ monomer units.

An accurate density functional theory for the vapor-liquid interface of associating chain molecules based on the statistical associating fluid theory for potentials of variable range

A Helmholtz free energy density functional is developed to describe the vapor-liquid interface of associating chain molecules. The functional is based on the statistical associating fluid theory with attractive potentials of variable range (SAFT-VR) for the homogenous fluid [A. Gil-Villegas, A. Galindo, P. J. Whitehead, S. J. Mills, G. Jackson, and A. N. Burgess, J. Chem. Phys. 106, 4168 (1997)]. A standard perturbative density functional theory (DFT) is constructed by partitioning the free energy density into a reference term (which incorporates all of the short-range interactions, and is treated locally) and an attractive perturbation (which incorporates the long-range dispersion interactions). In our previous work [F. J. Blas, E. Martin del Rio, E. de Miguel, and G. Jackson, Mol. Phys. 99, 1851 (2001); G. J. Gloor, F. J. Blas, E. Martin del Rio, E. de Miguel, and G. Jackson, Fluid Phase Equil. 194, 521 (2002)] we used a mean-field version of the theory (SAFT-HS) in which the pair correlations were neglected in the attractive term. This provides only a qualitative description of the vapor-liquid interface, due to the inadequate mean-field treatment of the vapor-liquid equilibria. Two different approaches are used to include the correlations in the attractive term: in the first, the free energy of the homogeneous fluid is partitioned such that the effect of correlations are incorporated in the local reference term; in the second, a density averaged correlation function is incorporated into the perturbative term in a similar way to that proposed by Toxvaerd [S. Toxvaerd, J. Chem. Phys. 64, 2863 (1976)]. The latter is found to provide the most accurate description of the vapor-liquid surface tension on comparison with new simulation data for a square-well fluid of variable range. The SAFT-VR DFT is used to examine the effect of molecular chain length and association on the surface tension. Different association schemes (dimerization, straight and branched chain formation, and network structures) are examined separately. The surface tension of the associating fluid is found to be bounded between the nonassociating and fully associated limits (both of which correspond to equivalent nonassociating systems). The temperature dependence of the surface tension is found to depend strongly on the balance between the strength and range of the association, and on the particular association scheme. In the case of a system with a strong but very localized association interaction, the surface tension exhibits the characteristic "s shaped" behavior with temperature observed in fluids such as water and alkanols. The various types of curves observed in real substances can be reproduced by the theory. It is very gratifying that a DFT based on SAFT-VR free energy can provide an accurate quantitative description of the surface tension of both the model and experimental systems.