Phase Separation in Mixtures of Prion-Like Low Complexity Domains is Driven by the Interplay of Homotypic and Heterotypic Interactions

Direct computations of viscoelastic moduli of biomolecular condensates

In vitro facsimiles of biomolecular condensates are formed by different types of intrinsically disordered proteins including prion-like low complexity domains (PLCDs). PLCD condensates are viscoelastic materials defined by time-dependent, sequence-specific complex shear moduli. Here, we show that viscoelastic moduli can be computed directly using a generalization of the Rouse model and information regarding intra- and inter-chain contacts that is extracted from equilibrium configurations of lattice-based Metropolis Monte Carlo (MMC) simulations. The key ingredient of the generalized Rouse model is the Zimm matrix that we compute from equilibrium MMC simulations. We compute two flavors of Zimm matrices, one referred to as the single-chain model that accounts only for intra-chain contacts, and the other referred to as a collective model, that accounts for inter-chain interactions. The single-chain model systematically overestimates the storage and loss moduli, whereas the collective model reproduces the measured moduli with greater fidelity. However, in the long time, low-frequency domain, a mixture of the two models proves to be most accurate. In line with the theory of Rouse, we find that a continuous distribution of relaxation times exists in condensates. The single crossover frequency between dominantly elastic versus dominantly viscous behaviors is influenced by the totality of the relaxation modes. Hence, our analysis suggests that viscoelastic fluid-like condensates are best described as generalized Maxwell fluids. Finally, we show that the complex shear moduli can be used to solve an inverse problem to obtain distributions of relaxation times that underlie the dynamics within condensates.

Sequence-specific interactions determine viscoelasticity and aging dynamics of protein condensates

Biomolecular condensates are viscoelastic materials. Here, we report results from investigations into molecular-scale determinants of sequence-encoded and age-dependent viscoelasticity of condensates formed by prion-like low-complexity domains (PLCDs). The terminally viscous forms of PLCD condensates are Maxwell fluids. Measured viscoelastic moduli of these condensates are reproducible using a Rouse-Zimm model that accounts for the network-like organization engendered by reversible physical crosslinks among PLCDs in the dense phase. Measurements and computations show that the strengths of aromatic inter-sticker interactions determine the sequence-specific amplitudes of elastic and viscous moduli as well as the timescales over which elastic properties dominate. PLCD condensates also undergo physical aging on sequence-specific timescales. This is driven by mutations to spacer residues that weaken the metastability of terminally viscous phases. The aging of PLCD condensates is accompanied by disorder-to-order transitions, leading to the formation of non-fibrillar, beta-sheet-containing, semi-crystalline, terminally elastic, Kelvin-Voigt solids. Our results suggest that sequence grammars, which refer to the identities of stickers versus spacers in PLCDs, have evolved to afford control over the metastabilities of terminally viscous fluid phases of condensates. This selection can, in some cases, render barriers for conversion from metastable fluids to globally stable solids to be insurmountable on functionally relevant timescales.

Heterotypic interactions in the dilute phase can drive co-condensation of prion-like low-complexity domains of FET proteins and mammalian SWI/SNF complex

Prion-like domains (PLDs) are low-complexity protein sequences enriched within nucleic acid-binding proteins including those involved in transcription and RNA processing. PLDs of FUS and EWSR1 play key roles in recruiting chromatin remodeler mammalian SWI/SNF complex to oncogenic FET fusion protein condensates. Here, we show that disordered low-complexity domains of multiple SWI/SNF subunits are prion-like with a strong propensity to undergo intracellular phase separation. These PLDs engage in sequence-specific heterotypic interactions with the PLD of FUS in the dilute phase at sub-saturation conditions, leading to the formation of PLD co-condensates. In the dense phase, homotypic and heterotypic PLD interactions are highly cooperative, resulting in the co-mixing of individual PLD phases and forming spatially homogeneous co-condensates. Heterotypic PLD-mediated positive cooperativity in protein-protein interaction networks is likely to play key roles in the co-phase separation of mSWI/SNF complex with transcription factors containing homologous low-complexity domains.

Conformational Properties of Polymers at Droplet Interfaces as Model Systems for Disordered Proteins

Polymer models serve as useful tools for studying the formation and physical properties of biomolecular condensates. In recent years, the interface dividing the dense and dilute phases of condensates has been discovered to be closely related to their functionality, but the conformational preferences of the constituent proteins remain unclear. To elucidate this, we perform molecular simulations of a droplet formed by phase separation of homopolymers as a surrogate model for the prion-like low-complexity domains. By systematically analyzing the polymer conformations at different locations in the droplet, we find that the chains become compact at the droplet interface compared with the droplet interior. Further, segmental analysis revealed that the end sections of the chains are enriched at the interface to maximize conformational entropy and are more expanded than the middle sections of the chains. We find that the majority of chain segments lie tangential to the droplet surface, and only the chain ends tend to align perpendicular to the interface. These trends also hold for the natural proteins FUS LC and LAF-1 RGG, which exhibit more compact chain conformations at the interface compared to the droplet interior. Our findings provide important insights into the interfacial properties of biomolecular condensates and highlight the value of using simple polymer physics models to understand the underlying mechanisms.

Discovering chromatin dysregulation induced by protein-coding perturbations at scale

Although population-scale databases have expanded to millions of protein-coding variants, insight into variant mechanisms has not kept pace. We present PROD-ATAC, a high-throughput method for discovering the effects of protein-coding variants on chromatin. A pooled library of variants is expressed in a disease-agnostic cell line, and single-cell ATAC resolves each variant's effect on chromatin. Using PROD-ATAC, we characterized the effects of >100 oncofusions (a class of cancer-causing chimeric proteins) and controls and revealed that pioneer activity is a common feature of fusions spanning an enormous range of fusion frequencies. Further, fusion-induced dysregulation can be context-agnostic as observed mechanisms often overlapped with cancer and cell-type specific prior knowledge. We also showed that gain-of-function pioneering is common among oncofusions. This work provides a global view of fusion-induced chromatin. We uncovered convergent mechanisms among disparate oncofusions and shared modes of dysregulation across different cancers. PROD-ATAC is generalizable to any set of protein-coding variants.

Conformational Properties of Polymers at Droplet Interfaces as Model Systems for Disordered Proteins

Polymer models serve as useful tools for studying the formation and physical properties of biomolecular condensates. In recent years, the interface dividing the dense and dilute phases of condensates has been discovered to be closely related to their functionality, but the conformational preferences of the constituent proteins remain unclear. To elucidate this, we perform molecular simulations of a droplet formed by liquid‒liquid phase separation of homopolymers, as a surrogate model for the prion-like low-complexity domains. By systematically analyzing the polymer conformations at different locations in the droplet, we find that the chains become compact at the droplet interface compared to the droplet interior. Further, segmental analysis revealed that the end sections of the chains are enriched at the interface to maximize conformational entropy, and are more expanded than the middle sections of the chains. We find that the majority of chain segments lie tangential to the droplet surface and only the chain ends tend to align perpendicular to the interface. These trends also hold for the natural proteins FUC LC and LAF-1 RGG, which exhibit more compact chain conformations at the interface compared with the droplet interior. Our findings provide important insights into the interfacial properties of biomolecular condensates and highlight the value of using simple polymer physics models to understand the underlying mechanisms.