Comparative roles of charge, π, and hydrophobic interactions in sequence-dependent phase separation of intrinsically disordered proteins.
17 Nov 2020-Proceedings of the National Academy of Sciences of the United States of America (National Academy of Sciences)-Vol. 117, Iss: 46, pp 28795-28805
TL;DR: In this article, the N-terminal intrinsically disordered region (IDR) of DEAD-box helicase Ddx4 was simulated to assess roles of electrostatic, hydrophobic, cation-π, and aromatic interactions in amino acid sequence-dependent LLPS.
Abstract: Endeavoring toward a transferable, predictive coarse-grained explicit-chain model for biomolecular condensates underlain by liquid-liquid phase separation (LLPS) of proteins, we conducted multiple-chain simulations of the N-terminal intrinsically disordered region (IDR) of DEAD-box helicase Ddx4, as a test case, to assess roles of electrostatic, hydrophobic, cation-π, and aromatic interactions in amino acid sequence-dependent LLPS. We evaluated three different residue-residue interaction schemes with a shared electrostatic potential. Neither a common hydrophobicity scheme nor one augmented with arginine/lysine-aromatic cation-π interactions consistently accounted for available experimental LLPS data on the wild-type, a charge-scrambled, a phenylalanine-to-alanine (FtoA), and an arginine-to-lysine (RtoK) mutant of Ddx4 IDR. In contrast, interactions based on contact statistics among folded globular protein structures reproduce the overall experimental trend, including that the RtoK mutant has a much diminished LLPS propensity. Consistency between simulation and experiment was also found for RtoK mutants of P-granule protein LAF-1, underscoring that, to a degree, important LLPS-driving π-related interactions are embodied in classical statistical potentials. Further elucidation is necessary, however, especially of phenylalanine's role in condensate assembly because experiments on FtoA and tyrosine-to-phenylalanine mutants suggest that LLPS-driving phenylalanine interactions are significantly weaker than posited by common statistical potentials. Protein-protein electrostatic interactions are modulated by relative permittivity, which in general depends on aqueous protein concentration. Analytical theory suggests that this dependence entails enhanced interprotein interactions in the condensed phase but more favorable protein-solvent interactions in the dilute phase. The opposing trends lead to only a modest overall impact on LLPS.
Citations
More filters
TL;DR: In this paper, the authors show that a wide range of cellular proteins, including FUS, TDP-43, Brd4, Sox2, and Annexin A11, can reenter a phase-separated regime at high salt concentrations.
Abstract: Liquid–liquid phase separation of proteins underpins the formation of membraneless compartments in living cells. Elucidating the molecular driving forces underlying protein phase transitions is therefore a key objective for understanding biological function and malfunction. Here we show that cellular proteins, which form condensates at low salt concentrations, including FUS, TDP-43, Brd4, Sox2, and Annexin A11, can reenter a phase-separated regime at high salt concentrations. By bringing together experiments and simulations, we demonstrate that this reentrant phase transition in the high-salt regime is driven by hydrophobic and non-ionic interactions, and is mechanistically distinct from the low-salt regime, where condensates are additionally stabilized by electrostatic forces. Our work thus sheds light on the cooperation of hydrophobic and non-ionic interactions as general driving forces in the condensation process, with important implications for aberrant function, druggability, and material properties of biomolecular condensates. Elucidating the molecular driving forces underlying liquid–liquid phase separation is a key objective for understanding biological function and malfunction. Here the authors show that a wide range of cellular proteins, including FUS, TDP-43, Brd4, Sox2, and Annexin A11, which form condensates at low salt concentrations, can reenter a phase-separated regime at high salt concentrations.
186 citations
TL;DR: In this article, a coarse-grained model of intrinsically disordered proteins (IDPs) with residue-level detail was developed based on an extensive set of experimental data on single-chain properties.
Abstract: Many intrinsically disordered proteins (IDPs) may undergo liquid-liquid phase separation (LLPS) and participate in the formation of membraneless organelles in the cell, thereby contributing to the regulation and compartmentalization of intracellular biochemical reactions. The phase behavior of IDPs is sequence dependent, and its investigation through molecular simulations requires protein models that combine computational efficiency with an accurate description of intramolecular and intermolecular interactions. We developed a general coarse-grained model of IDPs, with residue-level detail, based on an extensive set of experimental data on single-chain properties. Ensemble-averaged experimental observables are predicted from molecular simulations, and a data-driven parameter-learning procedure is used to identify the residue-specific model parameters that minimize the discrepancy between predictions and experiments. The model accurately reproduces the experimentally observed conformational propensities of a set of IDPs. Through two-body as well as large-scale molecular simulations, we show that the optimization of the intramolecular interactions results in improved predictions of protein self-association and LLPS.
102 citations
DOI•
24 Nov 2021TL;DR: In this article, a multiscale coarse-grained model for biomolecular phase separation is presented, which is parameterized from both atomistic simulations and bioinformatics data and accounts for the dominant role of π-π and hybrid cation-π/π-π−π interactions and the much stronger attractive contacts established by arginines than lysines.
Abstract: Various physics- and data-driven sequence-dependent protein coarse-grained models have been developed to study biomolecular phase separation and elucidate the dominant physicochemical driving forces. Here we present Mpipi, a multiscale coarse-grained model that describes almost quantitatively the change in protein critical temperatures as a function of amino acid sequence. The model is parameterized from both atomistic simulations and bioinformatics data and accounts for the dominant role of π–π and hybrid cation–π/π–π interactions and the much stronger attractive contacts established by arginines than lysines. We provide a comprehensive set of benchmarks for Mpipi and seven other residue-level coarse-grained models against experimental radii of gyration and quantitative in vitro phase diagrams, demonstrating that Mpipi predictions agree well with experiments on both fronts. Moreover, Mpipi can account for protein–RNA interactions, correctly predicts the multiphase behavior of a charge-matched poly-arginine/poly-lysine/RNA system, and recapitulates experimental liquid–liquid phase separation trends for sequence mutations on FUS, DDX4 and LAF-1 proteins. Combining bioinformatics data and atomistic simulations, this study develops a sequence-dependent coarse-grained model for biomolecular phase separation. This model achieves a quantitative agreement with experimental observations. Extensive benchmarks exemplify its performance.
82 citations
TL;DR: In this article, the hydropathy scale (HPS) coarse-grained (CG) model is used to simulate sequence-specific behavior of intrinsically disordered proteins (IDPs), including their liquid-liquid phase separation (LLPS).
Abstract: We present improvements to the hydropathy scale (HPS) coarse-grained (CG) model for simulating sequence-specific behavior of intrinsically disordered proteins (IDPs), including their liquid-liquid phase separation (LLPS). The previous model based on an atomistic hydropathy scale by Kapcha and Rossky (KR scale) is not able to capture some well-known LLPS trends such as reduced phase separation propensity upon mutations (R-to-K and Y-to-F). Here, we propose to use the Urry hydropathy scale instead, which was derived from the inverse temperature transitions in a model polypeptide with guest residues X. We introduce two free parameters to shift (Δ) and scale (µ) the overall interaction strengths for the new model (HPS-Urry) and use the experimental radius of gyration for a diverse group of IDPs to find their optimal values. Interestingly, many possible (Δ, µ) combinations can be used for typical IDPs, but the phase behavior of a low-complexity (LC) sequence FUS is only well described by one of these models, which highlights the need for a careful validation strategy based on multiple proteins. The CG HPS-Urry model should enable accurate simulations of protein LLPS and provide a microscopically detailed view of molecular interactions.
54 citations
TL;DR: This paper showed that α-synuclein (α-Syn) undergoes either spontaneous (instantaneous) or delayed liquid-liquid phase separation and subsequent liquid-to-solid phase transition, which leads to amyloid fibril formation.
Abstract: Liquid-liquid phase separation (LLPS) is a crucial phenomenon for the formation of functional membraneless organelles. However, LLPS is also responsible for protein aggregation in various neurodegenerative diseases such as amyotrophic lateral sclerosis, Alzheimer's disease, and Parkinson's disease (PD). Recently, several reports, including ours, have shown that α-synuclein (α-Syn) undergoes LLPS and a subsequent liquid-to-solid phase transition, which leads to amyloid fibril formation. However, how the environmental (and experimental) parameters modulate the α-Syn LLPS remains elusive. Here, we show that in vitro α-Syn LLPS is strongly dependent on the presence of salts, which allows charge neutralization at both terminal segments of protein and therefore promotes hydrophobic interactions supportive for LLPS. Using various purification methods and experimental conditions, we showed, depending upon conditions, α-Syn undergoes either spontaneous (instantaneous) or delayed LLPS. Furthermore, we delineate that the kinetics of liquid droplet formation (i.e., the critical concentration and critical time) is relative and can be modulated by the salt/counterion concentration, pH, presence of surface, PD-associated multivalent cations, and N-terminal acetylation, which are all known to regulate α-Syn aggregation in vitro. Together, our observations suggest that α-Syn LLPS and subsequent liquid-to-solid phase transition could be pathological, which can be triggered only under disease-associated conditions (high critical concentration and/or conditions promoting α-Syn self-assembly). This study will significantly improve our understanding of the molecular mechanisms of α-Syn LLPS and the liquid-to-solid transition.
52 citations
References
More filters
TL;DR: VMD is a molecular graphics program designed for the display and analysis of molecular assemblies, in particular biopolymers such as proteins and nucleic acids, which can simultaneously display any number of structures using a wide variety of rendering styles and coloring methods.
46,130 citations
TL;DR: In this article, the authors compared the Bernal Fowler (BF), SPC, ST2, TIPS2, TIP3P, and TIP4P potential functions for liquid water in the NPT ensemble at 25°C and 1 atm.
Abstract: Classical Monte Carlo simulations have been carried out for liquid water in the NPT ensemble at 25 °C and 1 atm using six of the simpler intermolecular potential functions for the water dimer: Bernal–Fowler (BF), SPC, ST2, TIPS2, TIP3P, and TIP4P. Comparisons are made with experimental thermodynamic and structural data including the recent neutron diffraction results of Thiessen and Narten. The computed densities and potential energies are in reasonable accord with experiment except for the original BF model, which yields an 18% overestimate of the density and poor structural results. The TIPS2 and TIP4P potentials yield oxygen–oxygen partial structure functions in good agreement with the neutron diffraction results. The accord with the experimental OH and HH partial structure functions is poorer; however, the computed results for these functions are similar for all the potential functions. Consequently, the discrepancy may be due to the correction terms needed in processing the neutron data or to an effect uniformly neglected in the computations. Comparisons are also made for self‐diffusion coefficients obtained from molecular dynamics simulations. Overall, the SPC, ST2, TIPS2, and TIP4P models give reasonable structural and thermodynamic descriptions of liquid water and they should be useful in simulations of aqueous solutions. The simplicity of the SPC, TIPS2, and TIP4P functions is also attractive from a computational standpoint.
33,683 citations
TL;DR: In this paper, a method is described to realize coupling to an external bath with constant temperature or pressure with adjustable time constants for the coupling, which can be easily extendable to other variables and to gradients, and can be applied also to polyatomic molecules involving internal constraints.
Abstract: In molecular dynamics (MD) simulations the need often arises to maintain such parameters as temperature or pressure rather than energy and volume, or to impose gradients for studying transport properties in nonequilibrium MD A method is described to realize coupling to an external bath with constant temperature or pressure with adjustable time constants for the coupling The method is easily extendable to other variables and to gradients, and can be applied also to polyatomic molecules involving internal constraints The influence of coupling time constants on dynamical variables is evaluated A leap‐frog algorithm is presented for the general case involving constraints with coupling to both a constant temperature and a constant pressure bath
25,256 citations
TL;DR: In this paper, a new Lagrangian formulation is introduced to make molecular dynamics (MD) calculations on systems under the most general externally applied, conditions of stress, which is well suited to the study of structural transformations in solids under external stress and at finite temperature.
Abstract: A new Lagrangian formulation is introduced. It can be used to make molecular dynamics (MD) calculations on systems under the most general, externally applied, conditions of stress. In this formulation the MD cell shape and size can change according to dynamical equations given by this Lagrangian. This new MD technique is well suited to the study of structural transformations in solids under external stress and at finite temperature. As an example of the use of this technique we show how a single crystal of Ni behaves under uniform uniaxial compressive and tensile loads. This work confirms some of the results of static (i.e., zero temperature) calculations reported in the literature. We also show that some results regarding the stress‐strain relation obtained by static calculations are invalid at finite temperature. We find that, under compressive loading, our model of Ni shows a bifurcation in its stress‐strain relation; this bifurcation provides a link in configuration space between cubic and hexagonal close packing. It is suggested that such a transformation could perhaps be observed experimentally under extreme conditions of shock.
13,937 citations