Principles Governing the Phase Separation of Multidomain Proteins

Transient Interdomain Interactions Modulate the Monomeric Structural Ensemble and Self-Assembly of Huntingtin Exon 1

Polyglutamine (polyQ) tract length expansion (≥ 36 residues) within the N-terminal exon-1 of Huntingtin (Httex1) leads to Huntington's disease, a neurodegenerative condition marked by the presence of intranuclear Htt inclusions. Notably, the polyQ tract in Httex1 is flanked by an N-terminal coiled-coil domain -N17 (17 amino acids), which promotes the formation of soluble oligomers and brings the aggregation-prone polyQ tracts in close proximity. However, the molecular mechanisms underlying the conversion of soluble oligomers into insoluble β-rich aggregates with increasing polyQ length, remain unclear. In this study, extensive atomistic molecular dynamics (MD) simulations (aggregate time ≈0.7 milliseconds) are performed to uncover the interplay between structural transformation and domain "cross-talk" on the conformational ensemble and oligomerization of Httex1 due to polyQ expansion. Notably, MD-derived ensembles of N17-Qn-P5 monomers validated against NMR indicated that in addition to elevated α-helicity, polyQ expansion also favored transient, interdomain (N17/polyQ) interactions which resulted in the emergence of β-sheet conformations. Further, interdomain interactions modulated the stability of N17-mediated polyQ dimers and promoted a heterogeneous dimerization landscape. Finally, it is observed that the intact C-terminal proline-rich domain (PRD) promoted condensation of Httex1 through self-interactions involving its P10/P11 tracts while also interacting with N17 to suppress its α-helicity.

Multiphasic protein condensation governed by shape and valency

Liquid-liquid phase separation (LLPS) of biological macromolecules leads to the formation of various membraneless organelles. The multilayered and multiphasic form of LLPS can mediate complex cellular functions; however, the determinants of its topological features are not fully understood. Herein, we focus on synaptic proteins consisting of Ca2+/calmodulin-dependent protein kinase II (CaMKII) and its interacting partners and present a computational model that reproduces forms of LLPS, including a form of two-phase condensates, phase-in-phase (PIP) organization. The model analyses reveal that the PIP formation requires competitive binding between the proteins. The PIP forms only when CaMKII has high valency and a short linker length. Such CaMKII exhibits low surface tension, a modular structure, and slow diffusion, enabling it to stay in small biochemical domains for a long time, which is necessary for synaptic plasticity. Thus, the computational modeling reveals new structure-function relationships for CaMKII as a synaptic memory unit.

The mechanobiology of biomolecular condensates

The central goal of mechanobiology is to understand how the mechanical forces and material properties of organelles, cells, and tissues influence biological processes and functions. Since the first description of biomolecular condensates, it was hypothesized that they obtain material properties that are tuned to their functions inside cells. Thus, they represent an intriguing playground for mechanobiology. The idea that biomolecular condensates exhibit diverse and adaptive material properties highlights the need to understand how different material states respond to external forces and whether these responses are linked to their physiological roles within the cell. For example, liquids buffer and dissipate, while solids store and transmit mechanical stress, and the relaxation time of a viscoelastic material can act as a mechanical frequency filter. Hence, a liquid-solid transition of a condensate in the force transmission pathway can determine how mechanical signals are transduced within and in-between cells, affecting differentiation, neuronal network dynamics, and behavior to external stimuli. Here, we first review our current understanding of the molecular drivers and how rigidity phase transitions are set forth in the complex cellular environment. We will then summarize the technical advancements that were necessary to obtain insights into the rich and fascinating mechanobiology of condensates, and finally, we will highlight recent examples of physiological liquid-solid transitions and their connection to specific cellular functions. Our goal is to provide a comprehensive summary of the field on how cells harness and regulate condensate mechanics to achieve specific functions.

Optimal scaling of protein-water interactions coupled with targeted torsional refinements yields balanced force fields suitable for simulations of single-chain folded proteins, disordered polypeptides, and protein-protein complexes

All-atom molecular dynamics (MD) simulations based on physics-based force fields, serve as an essential complement to experiments for investigating protein structure, dynamics, and interactions. Despite significant advances in force field development, achieving a consistent balance of molecular interactions that stabilize folded proteins and protein-protein complexes while simultaneously capturing the conformational dynamics of intrinsically disordered polypeptides (IDPs), remains challenging. In this work, we systematically evaluated two current state-of-the-art force fields (i) AMBER ff03ws, and (ii) AMBER ff99SBws, by comprehensively assessing their performance on both folded domains and IDPs. By selectively scaling side chain-water interactions for uncharged residues, the refined AMBER ff03w-sc force field demonstrated improved conformational stability of folded proteins while maintaining accurate representations of IDPs. However, AMBER ff03w-sc failed to correct the discrepancies in NMR-derived ps-ns timescale backbone dynamics associated with flexible loops. Interestingly, AMBER ff99SBws retained its structural stability despite the application of upscaled interactions with water for both sidechain and backbone atoms and displayed robust agreement with NMR-derived backbone dynamics. Further, a targeted refinement of glutamine backbone torsion parameters, yielded AMBER ff99SBws-STQ', which effectively resolved discrepancies associated with glutamine α-helicity predictions. Extensive validation against small angle X-ray scattering (SAXS) and NMR chemical shifts, revealed that both refined force fields accurately reproduced chain dimensions and secondary structure propensities of disordered peptides and prion-like domains. Importantly, both force fields reliably maintained the stability of protein-protein complexes over microsecond timescales. Our systematic refinement strategies provide improved accuracy and transferability for simulating diverse protein systems, from folded domains to IDPs and protein complexes.

How a disordered linker in the Polycomb protein Polyhomeotic tunes phase separation and oligomerization

Biomolecular condensates are increasingly appreciated for their function in organizing and regulating biochemical processes in cells, including chromatin function. Condensate formation and properties are encoded in protein sequence but the mechanisms linking sequence to macroscale properties are incompletely understood. Cross species comparisons can reveal mechanisms either because they identify conserved functions or because they point to important differences. Here we use in vitro reconstitution and molecular dynamics simulations to compare Drosophila and human sequences that regulate condensate formation driven by the sterile alpha motif (SAM) oligomerization domain in the Polyhomeotic (Ph) subunit of the chromatin regulatory complex PRC1. We discover evolutionarily diverged contacts between the conserved SAM and the disordered linker that connects it to the rest of Ph. Linker-SAM interactions increase oligomerization and regulate formation and properties of reconstituted condensates. Oligomerization affects condensate dynamics but, in most cases, has little effect on their formation. Linker-SAM interactions also affect condensate formation in Drosophila and human cells, and growth in Drosophila imaginal discs. Our data show how evolutionary sequence changes in linkers connecting conserved structured domains can alter condensate properties.

Recent advances in the synthesis and application of biomolecular condensates

Biomolecular condensates (BMCs) represent a group of organized and programmed systems that participate in gene transcription, chromosome organization, cell division, tumorigenesis, and aging. However, the understanding of BMCs in terms of internal organizations and external regulations remains at an early stage. Recently, novel approaches such as synthetic biology have been used for de novo synthesis of BMCs. These synthesized BMCs (SBMCs) driven by phase separation adeptly resemble the self-assembly and dynamics of natural BMCs, offering vast potentials in basic and applied research. This review introduces recent progresses in phase separation-induced SBMCs, attempting to elaborate on the intrinsic principles and regulatory methodologies used to construct SBMCs. Furthermore, the scientific applications of SBMCs are illustrated, as indicated by the studies of chromosome structure, pathogenesis, biomanufacturing, artificial cell design, and drug delivery. The controllable SBMCs offer a powerful tool for understanding metabolic regulations, cellular organizations, and disease-associated protein aggregations, raising both opportunities and challenges in the future of biomaterial, biotechnology, and biomedicine.

Emerging regulatory mechanisms and functions of biomolecular condensates: implications for therapeutic targets

Cells orchestrate their processes through complex interactions, precisely organizing biomolecules in space and time. Recent discoveries have highlighted the crucial role of biomolecular condensates-membrane-less assemblies formed through the condensation of proteins, nucleic acids, and other molecules-in driving efficient and dynamic cellular processes. These condensates are integral to various physiological functions, such as gene expression and intracellular signal transduction, enabling rapid and finely tuned cellular responses. Their ability to regulate cellular signaling pathways is particularly significant, as it requires a careful balance between flexibility and precision. Disruption of this balance can lead to pathological conditions, including neurodegenerative diseases, cancer, and viral infections. Consequently, biomolecular condensates have emerged as promising therapeutic targets, with the potential to offer novel approaches to disease treatment. In this review, we present the recent insights into the regulatory mechanisms by which biomolecular condensates influence intracellular signaling pathways, their roles in health and disease, and potential strategies for modulating condensate dynamics as a therapeutic approach. Understanding these emerging principles may provide valuable directions for developing effective treatments targeting the aberrant behavior of biomolecular condensates in various diseases.

Spontaneous Self-Organized Order Emerging From Intrinsically Disordered Protein Polymers

Intrinsically disordered proteins (IDPs) are proteins that, despite lacking a defined 3D structure, are capable of adopting dynamic conformations. This structural adaptability allows them to play not only essential roles in crucial cellular processes, such as subcellular organization or transcriptional control, but also in coordinating the assembly of macromolecules during different stages of development. Thus, in order to artificially replicate the complex processes of morphogenesis and their dynamics, it is crucial to have materials that recapitulate the structural plasticity of IDPs. In this regard, intrinsically disordered protein polymers (IDPPs) emerge as promising materials for engineering synthetic condensates and creating hierarchically self-assembled materials. IDPPs exhibit remarkable properties for their use in biofabrication, such as functional versatility, tunable sequence order-disorder, and the ability to undergo liquid-liquid phase separation (LLPS). Recent research has focused on harnessing the intrinsic disorder of IDPPs to design complex protein architectures with tailored properties. Taking advantage of their stimuli-responsiveness and degree of disorder, researchers have developed innovative strategies to control the self-assembly of IDPPs, resulting in the creation of hierarchically organized structures and intricate morphologies. In this review, we aim to provide an overview of the latest advances in the design and application of IDPP-based materials, shedding light on the fundamental principles that control their supramolecular assembly, and discussing their application in the biomedical and nanobiotechnological fields.

Deciphering the Monomeric and Dimeric Conformational Landscapes of the Full-Length TDP-43 and the Impact of the C-Terminal Domain

The aberrant aggregation of TAR DNA-binding protein 43 kDa (TDP-43) in cells leads to the pathogenesis of multiple fatal neurodegenerative diseases. Decoding the proposed initial transition between its functional dimeric and aggregation-prone monomeric states can potentially design a viable therapeutic strategy, which is presently limited by the lack of structural detail of the full-length TDP-43. To achieve a complete understanding of such a delicate phase space, we employed a multiscale simulation approach that unearths numerous crucial features, broadly summarized in two categories: (1) state-independent features that involve inherent chain collapsibility, rugged polymorphic landscape dictated by the terminal domains, high β-sheet propensity, structural integrity preserved by backbone-based intrachain hydrogen bonds and electrostatic forces, the prominence of the C-terminal domain in the intrachain cross-domain interfaces, and equal participation of hydrophobic and hydrophilic (charged and polar) residues in cross-domain interfaces; and (2) dimerization-modulated characteristics that encompass slower collapsing dynamics, restricted polymorphic landscape, the dominance of side chains in interchain hydrogen bonds, the appearance of the N-terminal domain in the dimer interface, and the prominence of hydrophilic (specifically polar) residues in interchain homo- and cross-domain interfaces. In our work, the ill-known C-terminal domain appears as the most crucial structure-dictating domain, which preferably populates a compact conformation with a high β-sheet propensity in its isolated state stabilized by intrabackbone hydrogen bonds, and these signatures are comparatively faded in its integrated form. Validation of our simulated observables by a complementary spectroscopic approach on multiple counts ensures the robustness of the computationally predicted features of the TDP-43 aggregation landscape.

Sequence-Encoded Spatiotemporal Dependence of Viscoelasticity of Protein Condensates Using Computational Microrheology

Many biomolecular condensates act as viscoelastic complex fluids with distinct cellular functions. Deciphering the viscoelastic behavior of biomolecular condensates can provide insights into their spatiotemporal organization and physiological roles within cells. Although there is significant interest in defining the role of condensate dynamics and rheology in physiological functions, the quantification of their time-dependent viscoelastic properties is limited and is mostly done through experimental rheological methods. Here, we demonstrate that a computational passive probe microrheology technique, coupled with continuum mechanics, can accurately characterize the linear viscoelasticity of condensates formed by intrinsically disordered proteins (IDPs). Using a transferable coarse-grained protein model, we first provide a physical basis for choosing optimal values that define the attributes of the probe particle, namely, its size and interaction strength with the residues in an IDP chain. We show that the technique captures the sequence-dependent viscoelasticity of heteropolymeric IDPs that differ in either sequence charge patterning or sequence hydrophobicity. We also illustrate the technique's potential in quantifying the spatial dependence of viscoelasticity in heterogeneous IDP condensates. The computational microrheology technique has important implications for investigating the time-dependent rheology of complex biomolecular architectures, resulting in the sequence-rheology-function relationship for condensates.

Membraneless organelles in health and disease: exploring the molecular basis, physiological roles and pathological implications

Once considered unconventional cellular structures, membraneless organelles (MLOs), cellular substructures involved in biological processes or pathways under physiological conditions, have emerged as central players in cellular dynamics and function. MLOs can be formed through liquid-liquid phase separation (LLPS), resulting in the creation of condensates. From neurodegenerative disorders, cardiovascular diseases, aging, and metabolism to cancer, the influence of MLOs on human health and disease extends widely. This review discusses the underlying mechanisms of LLPS, the biophysical properties that drive MLO formation, and their implications for cellular function. We highlight recent advances in understanding how the physicochemical environment, molecular interactions, and post-translational modifications regulate LLPS and MLO dynamics. This review offers an overview of the discovery and current understanding of MLOs and biomolecular condensate in physiological conditions and diseases. This article aims to deliver the latest insights on MLOs and LLPS by analyzing current research, highlighting their critical role in cellular organization. The discussion also covers the role of membrane-associated condensates in cell signaling, including those involving T-cell receptors, stress granules linked to lysosomes, and biomolecular condensates within the Golgi apparatus. Additionally, the potential of targeting LLPS in clinical settings is explored, highlighting promising avenues for future research and therapeutic interventions.

Thermostable Nucleoid Protein Cren7 Slides Along DNA and Rapidly Dissociates From DNA While Not Inhibiting the Sliding of Other DNA-binding Protein

A nucleoid protein Cren7 compacts DNA, contributing to the living of Crenarchaeum in high temperature environment. In this study, we investigated the dynamic behavior of Cren7 on DNA and its functional relation using single-molecule fluorescence microscopy. We found two mobility modes of Cren7, sliding along DNA and pausing on it, and the rapid dissociation kinetics from DNA. The salt dependence analysis suggests a sliding with continuous contact to DNA, rather than hopping/jumping. The mutational analysis demonstrates that Cren7 slides along DNA while Trp (W26) residue interacts with the DNA. Furthermore, Cren7 does not impede the target search by a model transcription factor p53, implying no significant interference to other DNA-binding proteins on DNA. At high concentration of Cren7, the molecules form large clusters on DNA via bridging, which compacts DNA. We discuss how the dynamic behavior of Cren7 on DNA enables DNA-compaction and protein-bypass functions.

RNA binding tunes the conformational plasticity and intradomain stability of TDP-43 tandem RNA recognition motifs

TAR DNA binding protein 43 (TDP-43) is a nuclear RNA/DNA-binding protein with pivotal roles in RNA-related processes such as splicing, transcription, transport, and stability. The high binding affinity and specificity of TDP-43 toward its cognate RNA sequences (GU-rich) is mediated by highly conserved residues in its tandem RNA recognition motif (RRM) domains (aa: 104-263). Importantly, the loss of RNA binding to the tandem RRMs caused by physiological stressors and chemical modifications promotes cytoplasmic mislocalization and pathological aggregation of TDP-43. Despite the substantial implications of RNA binding in TDP-43 function and pathology, its precise effects on the intradomain stability, and conformational dynamics of the tandem RRMs is not properly understood. Here, we employed all-atom molecular dynamics (MD) simulations to assess the effect of RNA binding on the conformational landscape and intradomain stability of TDP-43 tandem RRMs. RNA limits the overall conformational space of the tandem RRMs and promotes intradomain stability through a combination of specific base stacking interactions and transient electrostatic interactions. In contrast, tandem RRMs exhibit a high intrinsic conformational plasticity in the absence of RNA, which, surprisingly, is accompanied by a tendency of RRM1 to adopt partially unfolded conformations. Overall, our simulations reveal how RNA binding dynamically tunes the structural and conformational landscape of TDP-43 tandem RRMs, contributing to physiological function and mitigating pathological aggregation.

Unraveling the Bivalent and Rapid Interactions Between a Multivalent RNA Recognition Motif and RNA: A Kinetic Approach

The kinetics of the interaction between Musashi-1 (MSI1) and RNA have been characterized using surface plasmon resonance biosensor analysis. Truncated variants of human MSI1 encompassing the two homologous RNA recognition motifs (RRM1 and RRM2) in tandem (aa 1-200), and the two RRMs in isolation (aa 1-103 and aa 104-200, respectively) were produced. The proteins were injected over sensor surfaces with immobilized RNA, varying in sequence and length, and with one or two RRM binding motifs. The interactions of the individual RRMs with all RNA variants were well described by a 1:1 interaction model. The interaction between the MSI1 variant encompassing both RRM motifs was bivalent and rapid for all RNA variants. Due to difficulties in fitting this complex data using standard procedures, we devised a new method to quantify the interactions. It revealed that two RRMs in tandem resulted in a significantly longer residence time than a single RRM. It also showed that RNA with double UAG binding motifs and potential hairpin structures forms less stable bivalent complexes with MSI1 than the single UAG motif containing linear RNA. Substituting the UAG binding motif with a CAG sequence resulted in a reduction of the affinity of the individual RRMs, but for MSI1, this reduction was strongly enhanced, demonstrating the importance of bivalency for specificity. This study has provided new insights into the interaction between MSI1 and RNA and an understanding of how individual domains contribute to the overall interaction. It provides an explanation for why many RNA-binding proteins contain dual RRMs.

Sequence-encoded Spatiotemporal Dependence of Viscoelasticity of Protein Condensates Using Computational Microrheology

Many biomolecular condensates act as viscoelastic complex fluids with distinct cellular functions. Deciphering the viscoelastic behavior of biomolecular condensates can provide insights into their spatiotemporal organization and physiological roles within cells. Though there is significant interest in defining the role of condensate dynamics and rheology in physiological functions, the quantification of their time-dependent viscoelastic properties is limited and mostly done through experimental rheological methods. Here, we demonstrate that a computational passive probe microrheology technique, coupled with continuum mechanics, can accurately characterize the linear viscoelasticity of condensates formed by intrinsically disordered proteins (IDPs). Using a transferable coarse-grained protein model, we first provide a physical basis for choosing optimal values that define the attributes of the probe particle, namely its size and interaction strength with the residues in an IDP chain. We show that the technique captures the sequence-dependent viscoelasticity of heteropolymeric IDPs that differ either in sequence charge patterning or sequence hydrophobicity. We also illustrate the technique's potential in quantifying the spatial dependence of viscoelasticity in heterogeneous IDP condensates. The computational microrheology technique has important implications for investigating the time-dependent rheology of complex biomolecular architectures, resulting in the sequence-rheology-function relationship for condensates.

Phase separation as a new form of regulation in innate immunity

Innate immunity is essential for the host against pathogens, cancer, and autoimmunity. The innate immune system encodes many sensor, adaptor, and effector proteins and relies on the assembly of higher-order signaling complexes to activate immune defense. Recent evidence demonstrates that many of the core complexes involved in innate immunity are organized as liquid-like condensates through a mechanism known as phase separation. Here, we discuss phase-separated condensates and their diverse functions. We compare the biochemical, structural, and mechanistic details of solid and liquid-like assemblies to explore the role of phase separation in innate immunity. We summarize the emerging evidence for the hypothesis that phase separation is a conserved mechanism that controls immune responses across the tree of life. The discovery of phase separation in innate immunity provides a new foundation to explain the rules that govern immune system activation and will enable the development of therapeutics to treat immune-related diseases properly.

Determinants that enable disordered protein assembly into discrete condensed phases

Cells harbour numerous mesoscale membraneless compartments that house specific biochemical processes and perform distinct cellular functions. These protein- and RNA-rich bodies are thought to form through multivalent interactions among proteins and nucleic acids, resulting in demixing via liquid-liquid phase separation. Proteins harbouring intrinsically disordered regions (IDRs) predominate in membraneless organelles. However, it is not known whether IDR sequence alone can dictate the formation of distinct condensed phases. We identified a pair of IDRs capable of forming spatially distinct condensates when expressed in cells. When reconstituted in vitro, these model proteins do not co-partition, suggesting condensation specificity is encoded directly in the polypeptide sequences. Through computational modelling and mutagenesis, we identified the amino acids and chain properties governing homotypic and heterotypic interactions that direct selective condensation. These results form the basis of physicochemical principles that may direct subcellular organization of IDRs into specific condensates and reveal an IDR code that can guide construction of orthogonal membraneless compartments.

Expanding the molecular language of protein liquid-liquid phase separation

Understanding the relationship between a polypeptide sequence and its phase separation has important implications for analysing cellular function, treating disease and designing novel biomaterials. Several sequence features have been identified as drivers for protein liquid-liquid phase separation (LLPS), schematized as a 'molecular grammar' for LLPS. Here we further probe how sequence modulates phase separation and the material properties of the resulting condensates, targeting sequence features previously overlooked in the literature. We generate sequence variants of a repeat polypeptide with either no charged residues, high net charge, no glycine residues or devoid of aromatic or arginine residues. All but one of 12 variants exhibited LLPS, albeit to different extents, despite substantial differences in composition. Furthermore, we find that all the condensates formed behaved like viscous fluids, despite large differences in their viscosities. Our results support the model of multiple interactions between diverse residue pairs-not just a handful of residues-working in tandem to drive the phase separation and dynamics of condensates.

Uncovering the molecular interactions underlying MBD2 and MBD3 phase separation

Chromatin organization controls DNA's accessibility to regulatory factors to influence gene expression. Heterochromatin, or transcriptionally silent chromatin enriched in methylated DNA and methylated histone tails, self-assembles through multivalent interactions with its associated proteins into a condensed, but dynamic state. Liquid-liquid phase separation (LLPS) of key heterochromatin regulators, such as heterochromatin protein 1 (HP1), plays an essential role in heterochromatin assembly and function. Methyl-CpG-binding protein 2 (MeCP2), the most studied member of the methyl-CpG-binding domain (MBD) family of proteins, has been recently shown to undergo LLPS in the absence and presence of methylated DNA. These studies provide a new mechanistic framework for understanding the role of methylated DNA and its readers in heterochromatin formation. However, the details of the molecular interactions by which other MBD family members undergo LLPS to mediate genome organization and transcriptional regulation are not fully understood. Here, we focus on two MBD proteins, MBD2 and MBD3, that have distinct but interdependent roles in gene regulation. Using an integrated computational and experimental approach, we uncover the homotypic and heterotypic interactions governing MBD2 and MBD3 phase separation and DNA's influence on this process. We show that despite sharing the highest sequence identity and structural homology among all the MBD protein family members, MBD2 and MBD3 exhibit differing residue patterns resulting in distinct phase separation mechanisms. Understanding the molecular underpinnings of MBD protein condensation offers insights into the higher-order, LLPS-mediated organization of heterochromatin.

Kinetic trapping organizes actin filaments within liquid-like protein droplets

Several actin-binding proteins (ABPs) phase separate to form condensates capable of curating the actin network shapes. Here, we use computational modeling to understand the principles of actin network organization within VASP condensate droplets. Our simulations reveal that the different actin shapes, namely shells, rings, and mixture states are highly dependent on the kinetics of VASP-actin interactions, suggesting that they arise from kinetic trapping. Specifically, we show that reducing the residence time of VASP on actin filaments reduces degree of bundling, thereby promoting assembly of shells rather than rings. We validate the model predictions experimentally using a VASP-mutant with decreased bundling capability. Finally, we investigate the ring opening within deformed droplets and found that the sphere-to-ellipsoid transition is favored under a wide range of filament lengths while the ellipsoid-to-rod transition is only permitted when filaments have a specific range of lengths. Our findings highlight key mechanisms of actin organization within phase-separated ABPs.

Interplay between charge distribution and DNA in shaping HP1 paralog phase separation and localization

The heterochromatin protein 1 (HP1) family is a crucial component of heterochromatin with diverse functions in gene regulation, cell cycle control, and cell differentiation. In humans, there are three paralogs, HP1α, HP1β, and HP1γ, which exhibit remarkable similarities in their domain architecture and sequence properties. Nevertheless, these paralogs display distinct behaviors in liquid-liquid phase separation (LLPS), a process linked to heterochromatin formation. Here, we employ a coarse-grained simulation framework to uncover the sequence features responsible for the observed differences in LLPS. We highlight the significance of the net charge and charge patterning along the sequence in governing paralog LLPS propensities. We also show that both highly conserved folded and less-conserved disordered domains contribute to the observed differences. Furthermore, we explore the potential co-localization of different HP1 paralogs in multicomponent assemblies and the impact of DNA on this process. Importantly, our study reveals that DNA can significantly reshape the stability of a minimal condensate formed by HP1 paralogs due to competitive interactions of HP1α with HP1β and HP1γ versus DNA. In conclusion, our work highlights the physicochemical nature of interactions that govern the distinct phase-separation behaviors of HP1 paralogs and provides a molecular framework for understanding their role in chromatin organization.

Coiled-coil domains are sufficient to drive liquid-liquid phase separation in protein models

Liquid-liquid phase separation (LLPS) is thought to be a main driving force in the formation of membraneless organelles. Examples of such organelles include the centrosome, central spindle, and stress granules. Recently, it has been shown that coiled-coil (CC) proteins, such as the centrosomal proteins pericentrin, spd-5, and centrosomin, might be capable of LLPS. CC domains have physical features that could make them the drivers of LLPS, but it is unknown if they play a direct role in the process. We developed a coarse-grained simulation framework for investigating the LLPS propensity of CC proteins, in which interactions that support LLPS arise solely from CC domains. We show, using this framework, that the physical features of CC domains are sufficient to drive LLPS of proteins. The framework is specifically designed to investigate how the number of CC domains, as well as the multimerization state of CC domains, can affect LLPS. We show that small model proteins with as few as two CC domains can phase separate. Increasing the number of CC domains up to four per protein can somewhat increase LLPS propensity. We demonstrate that trimer-forming and tetramer-forming CC domains have a dramatically higher LLPS propensity than dimer-forming coils, which shows that multimerization state has a greater effect on LLPS than the number of CC domains per protein. These data support the hypothesis of CC domains as drivers of protein LLPS, and have implications in future studies to identify the LLPS-driving regions of centrosomal and central spindle proteins.

Triblock Proteins with Weakly Dimerizing Terminal Blocks and an Intrinsically Disordered Region for Rational Design of Condensate Properties

Condensates are molecular assemblies that are formed through liquid-liquid phase separation and play important roles in many biological processes. The rational design of condensate formation and their properties is central to applications, such as biosynthetic materials, synthetic biology, and for understanding cell biology. Protein engineering is used to make a triblock structure with varying terminal blocks of folded proteins on both sides of an intrinsically disordered mid-region. Dissociation constants are determined in the range of micromolar to millimolar for a set of proteins suitable for use as terminal blocks. Varying the weak dimerization of terminal blocks leads to an adjustable tendency for condensate formation while keeping the intrinsically disordered region constant. The dissociation constants of the terminal domains correlate directly with the tendency to undergo liquid-liquid phase separation. Differences in physical properties, such as diffusion rate are not directly correlated with the strength of dimerization but can be understood from the properties and interplay of the constituent blocks. The work demonstrates the importance of weak interactions in condensate formation and shows a principle for protein design that will help in fabricating functional condensates in a predictable and rational way.

Sequence-dependent material properties of biomolecular condensates and their relation to dilute phase conformations

Material properties of phase-separated biomolecular condensates, enriched with disordered proteins, dictate many cellular functions. Contrary to the progress made in understanding the sequence-dependent phase separation of proteins, little is known about the sequence determinants of condensate material properties. Using the hydropathy scale and Martini models, we computationally decipher these relationships for charge-rich disordered protein condensates. Our computations yield dynamical, rheological, and interfacial properties of condensates that are quantitatively comparable with experimentally characterized condensates. Interestingly, we find that the material properties of model and natural proteins respond similarly to charge segregation, despite different sequence compositions. Molecular interactions within the condensates closely resemble those within the single-chain ensembles. Consequently, the material properties strongly correlate with molecular contact dynamics and single-chain structural properties. We demonstrate the potential to harness the sequence characteristics of disordered proteins for predicting and engineering the material properties of functional condensates, with insights from the dilute phase properties.

Coarse-Grained Models to Study Protein-DNA Interactions and Liquid-Liquid Phase Separation

Recent advances in coarse-grained (CG) computational models for DNA have enabled molecular-level insights into the behavior of DNA in complex multiscale systems. However, most existing CG DNA models are not compatible with CG protein models, limiting their applications for emerging topics such as protein-nucleic acid assemblies. Here, we present a new computationally efficient CG DNA model. We first use experimental data to establish the model's ability to predict various aspects of DNA behavior, including melting thermodynamics and relevant local structural properties such as the major and minor grooves. We then employ an all-atom hydropathy scale to define nonbonded interactions between protein and DNA sites, to make our DNA model compatible with an existing CG protein model (HPS-Urry), which is extensively used to study protein phase separation, and show that our new model reasonably reproduces the experimental binding affinity for a prototypical protein-DNA system. To further demonstrate the capabilities of this new model, we simulate a full nucleosome with and without histone tails, on a microsecond time scale, generating conformational ensembles and provide molecular insights into the role of histone tails in influencing the liquid-liquid phase separation (LLPS) of HP1α proteins. We find that histone tails interact favorably with DNA, influencing the conformational ensemble of the DNA and antagonizing the contacts between HP1α and DNA, thus affecting the ability of DNA to promote LLPS of HP1α. These findings shed light on the complex molecular framework that fine-tunes the phase transition properties of heterochromatin proteins and contributes to heterochromatin regulation and function. Overall, the CG DNA model presented here is suitable to facilitate micrometer-scale studies with sub-nm resolution in many biological and engineering applications and can be used to investigate protein-DNA complexes, such as nucleosomes, or LLPS of proteins with DNA, enabling a mechanistic understanding of how molecular information may be propagated at the genome level.

A complex network of interdomain interactions underlies the conformational ensemble of monomeric TDP-43 and modulates its phase behavior

TAR DNA-binding protein 43 (TDP-43) is a multidomain protein involved in the regulation of RNA metabolism, and its aggregates have been observed in neurodegenerative diseases, including amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). Numerous studies indicate TDP-43 can undergo liquid-liquid phase separation (LLPS) in vitro and is a component of biological condensates. Homo-oligomerization via the folded N-terminal domain (aa:1-77) and the conserved helical region (aa:319-341) of the disordered, C-terminal domain is found to be an important driver of TDP-43 phase separation. However, a comprehensive molecular view of TDP-43 phase separation, particularly regarding the nature of heterodomain interactions, is lacking due to the challenges associated with its stability and purification. Here, we utilize all-atom and coarse-grained (CG) molecular dynamics (MD) simulations to uncover the network of interdomain interactions implicated in TDP-43 phase separation. All-atom simulations uncovered the presence of transient, interdomain interactions involving flexible linkers, RNA-recognition motif (RRM) domains and a charged segment of disordered C-terminal domain (CTD). CG simulations indicate these inter-domain interactions which affect the conformational landscape of TDP-43 in the dilute phase are also prevalent in the condensed phase. Finally, sequence and surface charge distribution analysis coupled with all-atom simulations (at high salt) confirmed that the transient interdomain contacts are predominantly electrostatic in nature. Overall, our findings from multiscale simulations lead to a greater appreciation of the complex interaction network underlying the structural landscape and phase separation of TDP-43.

Sequence-Dependent Material Properties of Biomolecular Codensates and their Relation to Dilute Phase Conformations

Material properties of phase-separated biomolecular assemblies, enriched with disordered proteins, dictate their ability to participate in many cellular functions. Despite the significant effort dedicated to understanding how the sequence of the disordered protein drives its phase separation to form condensates, little is known about the sequence determinants of condensate material properties. Here, we computationally decipher these relationships for charged disordered proteins using model sequences comprised of glutamic acid and lysine residues as well as naturally occurring sequences of LAF1's RGG domain and DDX4's N-terminal domain. We do so by delineating how the arrangement of oppositely charged residues within these sequences influences the dynamical, rheological, and interfacial properties of the condensed phase through equilibrium and non-equilibrium molecular simulations using the hydropathy scale and Martini models. Our computations yield material properties that are quantitatively comparable with experimentally characterized condensate systems. Interestingly, we find that the material properties of both the model and natural proteins respond similarly to the segregation of charges, despite their very different sequence compositions. Condensates of the highly charge-segregated sequences exhibit slower dynamics than the uniformly charge-patterned sequences, because of their comparatively long-lived molecular contacts between oppositely charged residues. Surprisingly, the molecular interactions within the condensate are highly similar to those within a single-chain for all sequences. Consequently, the condensate material properties of charged disordered proteins are strongly correlated with their dense phase contact dynamics and their single-chain structural properties. Our findings demonstrate the potential to harness the sequence characteristics of disordered proteins for predicting and engineering the material properties of functional condensates, with insights from the dilute phase properties.

Interplay between charge distribution and DNA in shaping HP1 paralog phase separation and localization

The heterochromatin protein 1 (HP1) family is a crucial component of heterochromatin with diverse functions in gene regulation, cell cycle control, and cell differentiation. In humans, there are three paralogs, HP1α, HP1β, and HP1γ, which exhibit remarkable similarities in their domain architecture and sequence properties. Nevertheless, these paralogs display distinct behaviors in liquid-liquid phase separation (LLPS), a process linked to heterochromatin formation. Here, we employ a coarse-grained simulation framework to uncover the sequence features responsible for the observed differences in LLPS. We highlight the significance of the net charge and charge patterning along the sequence in governing paralog LLPS propensities. We also show that both highly conserved folded and less-conserved disordered domains contribute to the observed differences. Furthermore, we explore the potential co-localization of different HP1 paralogs in multicomponent assemblies and the impact of DNA on this process. Importantly, our study reveals that DNA can significantly reshape the stability of a minimal condensate formed by HP1 paralogs due to competitive interactions of HP1α with HP1β and HP1γ versus DNA. In conclusion, our work highlights the physicochemical nature of interactions that govern the distinct phase-separation behaviors of HP1 paralogs and provides a molecular framework for understanding their role in chromatin organization.

Programmable de novo designed coiled coil-mediated phase separation in mammalian cells

Membraneless liquid compartments based on phase-separating biopolymers have been observed in diverse cell types and attributed to weak multivalent interactions predominantly based on intrinsically disordered domains. The design of liquid-liquid phase separated (LLPS) condensates based on de novo designed tunable modules that interact in a well-understood, controllable manner could improve our understanding of this phenomenon and enable the introduction of new features. Here we report the construction of CC-LLPS in mammalian cells, based on designed coiled-coil (CC) dimer-forming modules, where the stability of CC pairs, their number, linkers, and sequential arrangement govern the transition between diffuse, liquid and immobile condensates and are corroborated by coarse-grained molecular simulations. Through modular design, we achieve multiple coexisting condensates, chemical regulation of LLPS, condensate fusion, formation from either one or two polypeptide components or LLPS regulation by a third polypeptide chain. These findings provide further insights into the principles underlying LLPS formation and a design platform for controlling biological processes.

Direct observation of oriented behavior of actin filaments interacting with desmin intermediate filaments

BACKGROUND

Associations between actin filaments (AFs) and intermediate filaments (IFs) are frequently observed in living cells. The crosstalk between these cytoskeletal components underpins cellular organization and dynamics; however, the molecular basis of filamentous interactions is not fully understood. Here, we describe the mode of interaction between AFs and desmin IFs (DIFs) in a reconstituted in vitro system.

METHODS

AFs (rabbit skeletal muscle) and DIFs (chicken gizzard) were labeled with fluorescent dyes. DIFs were immobilized on a heavy meromyosin (HMM)-coated collodion surface. HMM-driven AFs with ATP hydrolysis was assessed in the presence of DIFs. Images of single filaments were obtained using fluorescence microscopy. Vector changes in the trajectories of single AFs were calculated from microscopy images.

RESULTS

AF speed transiently decreased upon contact with DIF. The difference between the incoming and outgoing angles of a moving AF broadened upon contact with a DIF. A smaller incoming angle tended to result in a smaller outgoing angle in a nematic manner. The percentage of moving AFs decreased with an increasing DIF density, but the speed of the moving AFs was similar to that in the no-desmin control. An abundance of DIFs tended to exclude AFs from the HMM-coated surfaces.

CONCLUSIONS

DIFs agitate the movement of AFs with the orientation. DIFs can bind to HMMs and weaken actin-myosin interactions.

GENERAL SIGNIFICANCE

The study indicates that apart from the binding strength, the accumulation of weak interactions characteristic of filamentous structures may affect the dynamic organization of cell architecture.

Multiscale simulations reveal TDP-43 molecular-level interactions driving condensation

The RNA-binding protein TDP-43 is associated with mRNA processing and transport from the nucleus to the cytoplasm. TDP-43 localizes in the nucleus as well as accumulating in cytoplasmic condensates such as stress granules. Aggregation and formation of amyloid-like fibrils of cytoplasmic TDP-43 are hallmarks of numerous neurodegenerative diseases, most strikingly present in >90% of amyotrophic lateral sclerosis (ALS) patients. If excessive accumulation of cytoplasmic TDP-43 causes, or is caused by, neurodegeneration is presently not known. In this work, we use molecular dynamics simulations at multiple resolutions to explore TDP-43 self- and cross-interaction dynamics. A full-length molecular model of TDP-43, all 414 amino acids, was constructed from select structures of the protein functional domains (N-terminal domain, and two RNA recognition motifs, RRM1 and RRM2) and modeling of disordered connecting loops and the low complexity glycine-rich C-terminus domain. All-atom CHARMM36m simulations of single TDP-43 proteins served as guides to construct a coarse-grained Martini 3 model of TDP-43. The Martini model and a coarser implicit solvent C⍺ model, optimized for disordered proteins, were subsequently used to probe TDP-43 interactions; self-interactions from single-chain full-length TDP-43 simulations, cross-interactions from simulations with two proteins and simulations with assemblies of dozens to hundreds of proteins. Our findings illustrate the utility of different modeling scales for accessing TDP-43 molecular-level interactions and suggest that TDP-43 has numerous interaction preferences or patterns, exhibiting an overall strong, but dynamic, association and driving the formation of biomolecular condensates.

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.

A Minispidroin Guides the Molecular Design for Cellular Condensation Mechanisms in S. cerevisiae

Structural engineering of molecules for condensation is an emerging technique within synthetic biology. Liquid-liquid phase separation of biomolecules leading to condensation is a central step in the assembly of biological materials into their functional forms. Intracellular condensates can also function within cells in a regulatory manner to facilitate reaction pathways and to compartmentalize interactions. We need to develop a strong understanding of how to design molecules for condensates and how their in vivo-in vitro properties are related. The spider silk protein NT2RepCT undergoes condensation during its fiber-forming process. Using parallel in vivo and in vitro characterization, in this study, we mapped the effects of intracellular conditions for NT2RepCT and its several structural variants. We found that intracellular conditions may suppress to some extent condensation whereas molecular crowding affects both condensate properties and their formation. Intracellular characterization of protein condensation allowed experiments on pH effects and solubilization to be performed within yeast cells. The growth of intracellular NT2RepCT condensates was restricted, and Ostwald ripening was not observed in yeast cells, in contrast to earlier observations in E. coli. Our results lead the way to using intracellular condensation to screen for properties of molecular assembly. For characterizing different structural variants, intracellular functional characterization can eliminate the need for time-consuming batch purification and in vitro condensation. Therefore, we suggest that the in vivo-in vitro understanding will become useful in, e.g., high-throughput screening for molecular functions and in strategies for designing tunable intracellular condensates.

Determinants of Disordered Protein Co-Assembly Into Discrete Condensed Phases

Cells harbor numerous mesoscale membraneless compartments that house specific biochemical processes and perform distinct cellular functions. These protein and RNA-rich bodies are thought to form through multivalent interactions among proteins and nucleic acids resulting in demixing via liquid-liquid phase separation (LLPS). Proteins harboring intrinsically disordered regions (IDRs) predominate in membraneless organelles. However, it is not known whether IDR sequence alone can dictate the formation of distinct condensed phases. We identified a pair of IDRs capable of forming spatially distinct condensates when expressed in cells. When reconstituted in vitro, these model proteins do not co-partition, suggesting condensation specificity is encoded directly in the polypeptide sequences. Through computational modeling and mutagenesis, we identified the amino acids and chain properties governing homotypic and heterotypic interactions that direct selective condensation. These results form the basis of physicochemical principles that may direct subcellular organization of IDRs into specific condensates and reveal an IDR code that can guide construction of orthogonal membraneless compartments.

A synergy between site-specific and transient interactions drives the phase separation of a disordered, low-complexity domain

TAR DNA-binding protein 43 (TDP-43) is involved in key processes in RNA metabolism and is frequently implicated in many neurodegenerative diseases, including amyotrophic lateral sclerosis and frontotemporal dementia. The prion-like, disordered C-terminal domain (CTD) of TDP-43 is aggregation-prone, can undergo liquid-liquid phase separation (LLPS) in isolation, and is critical for phase separation (PS) of the full-length protein under physiological conditions. While a short conserved helical region (CR, spanning residues 319-341) promotes oligomerization and is essential for LLPS, aromatic residues in the flanking disordered regions (QN-rich, IDR1/2) are also found to play a critical role in PS and aggregation. Compared with other phase-separating proteins, TDP-43 CTD has a notably distinct sequence composition including many aliphatic residues such as methionine and leucine. Aliphatic residues were previously suggested to modulate the apparent viscosity of the resulting phases, but their direct contribution toward CTD phase separation has been relatively ignored. Using multiscale simulations coupled with in vitro saturation concentration (csat) measurements, we identified the importance of aromatic residues while also suggesting an essential role for aliphatic methionine residues in promoting single-chain compaction and LLPS. Surprisingly, NMR experiments showed that transient interactions involving phenylalanine and methionine residues in the disordered flanking regions can directly enhance site-specific, CR-mediated intermolecular association. Overall, our work highlights an underappreciated mode of biomolecular recognition, wherein both transient and site-specific hydrophobic interactions act synergistically to drive the oligomerization and phase separation of a disordered, low-complexity domain.

Protein nanocondensates: the next frontier

Over the past decade, myriads of studies have highlighted the central role of protein condensation in subcellular compartmentalization and spatiotemporal organization of biological processes. Conceptually, protein condensation stands at the highest level in protein structure hierarchy, accounting for the assembly of bodies ranging from thousands to billions of molecules and for densities ranging from dense liquids to solid materials. In size, protein condensates range from nanocondensates of hundreds of nanometers (mesoscopic clusters) to phase-separated micron-sized condensates. In this review, we focus on protein nanocondensation, a process that can occur in subsaturated solutions and can nucleate dense liquid phases, crystals, amorphous aggregates, and fibers. We discuss the nanocondensation of proteins in the light of general physical principles and examine the biophysical properties of several outstanding examples of nanocondensation. We conclude that protein nanocondensation cannot be fully explained by the conceptual framework of micron-scale biomolecular condensation. The evolution of nanocondensates through changes in density and order is currently under intense investigation, and this should lead to the development of a general theoretical framework, capable of encompassing the full range of sizes and densities found in protein condensates.

Challenges in studying the liquid-to-solid phase transitions of proteins using computer simulations

"Membraneless organelles," also referred to as biomolecular condensates, perform a variety of cellular functions and their dysregulation is implicated in cancer and neurodegeneration. In the last two decades, liquid-liquid phase separation (LLPS) of intrinsically disordered and multidomain proteins has emerged as a plausible mechanism underlying the formation of various biomolecular condensates. Further, the occurrence of liquid-to-solid transitions within liquid-like condensates may give rise to amyloid structures, implying a biophysical link between phase separation and protein aggregation. Despite significant advances, uncovering the microscopic details of liquid-to-solid phase transitions using experiments remains a considerable challenge and presents an exciting opportunity for the development of computational models which provide valuable, complementary insights into the underlying phenomenon. In this review, we first highlight recent biophysical studies which provide new insights into the molecular mechanisms underlying liquid-to-solid (fibril) phase transitions of folded, disordered and multi-domain proteins. Next, we summarize the range of computational models used to study protein aggregation and phase separation. Finally, we discuss recent computational approaches which attempt to capture the underlying physics of liquid-to-solid transitions along with their merits and shortcomings.

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.

Coiled-coil domains are sufficient to drive liquid-liquid phase separation of proteins in molecular models

Liquid-liquid phase separation (LLPS) is thought to be a main driving force in the formation of membraneless organelles. Examples of such organelles include the centrosome, central spindle, and stress granules. Recently, it has been shown that coiled-coil (CC) proteins, such as the centrosomal proteins pericentrin, spd-5, and centrosomin, might be capable of LLPS. CC domains have physical features that could make them the drivers of LLPS, but it is unknown if they play a direct role in the process. We developed a coarse-grained simulation framework for investigating the LLPS propensity of CC proteins, in which interactions which support LLPS arise solely from CC domains. We show, using this framework, that the physical features of CC domains are sufficient to drive LLPS of proteins. The framework is specifically designed to investigate how the number of CC domains, as well as multimerization state of CC domains, can affect LLPS. We show that small model proteins with as few as two CC domains can phase separate. Increasing the number of CC domains up to four per protein can somewhat increase LLPS propensity. We demonstrate that trimer-forming and tetramer-forming CC domains have a dramatically higher LLPS propensity than dimer-forming coils, which shows that multimerization state has a greater effect on LLPS than the number of CC domains per protein. These data support the hypothesis of CC domains as drivers of protein LLPS, and has implications in future studies to identify the LLPS-driving regions of centrosomal and central spindle proteins.

Liquid - liquid phase separation of tau: Driving forces, regulation, and biological implications

The past 15 years have witnessed an explosion in the studies of biomolecular condensates that are implicated in numerous biological processes and play vital roles in human health and diseases. Recent findings demonstrate that the microtubule-associated protein tau forms liquid condensates through liquid-liquid phase separation (LLPS) in in vitro experiments using purified recombinant proteins and cell-based experiments. Although in vivo studies are lacking, liquid condensates have emerged as an important assembly state of physiological and pathological tau and LLPS can regulate the function of microtubules, mediate stress granule formation, and accelerate tau amyloid aggregation. In this review, we summarize recent advances in tau LLPS, aiming to unveiling the delicate interactions driving tau LLPS. We further discuss the association of tau LLPS with physiology and disease in the context of the sophisticated regulation of tau LLPS. Deciphering the mechanisms underlying tau LLPS and the liquid-to-solid transition enables rational design of molecules that inhibit or delay the formation of tau solid species, thus providing novel targeted therapeutic strategies for tauopathies.

A coarse-grained DNA model to study protein-DNA interactions and liquid-liquid phase separation

Recent advances in coarse-grained (CG) computational models for DNA have enabled molecular-level insights into the behavior of DNA in complex multiscale systems. However, most existing CG DNA models are not compatible with CG protein models, limiting their applications for emerging topics such as protein-nucleic acid assemblies. Here, we present a new computationally efficient CG DNA model. We first use experimental data to establish the model's ability to predict various aspects of DNA behavior, including melting thermodynamics and relevant local structural properties such as the major and minor grooves. We then employ an all-atom hydropathy scale to define non-bonded interactions between protein and DNA sites, to make our DNA model compatible with an existing CG protein model (HPS-Urry), that is extensively used to study protein phase separation, and show that our new model reasonably reproduces the experimental binding affinity for a prototypical protein-DNA system. To further demonstrate the capabilities of this new model, we simulate a full nucleosome with and without histone tails, on a microsecond timescale, generating conformational ensembles and provide molecular insights into the role of histone tails in influencing the liquid-liquid phase separation (LLPS) of HP1α proteins. We find that histone tails interact favorably with DNA, influencing the conformational ensemble of the DNA and antagonizing the contacts between HP1α and DNA, thus affecting the ability of DNA to promote LLPS of HP1α. These findings shed light on the complex molecular framework that fine-tunes the phase transition properties of heterochromatin proteins and contributes to heterochromatin regulation and function. Overall, the CG DNA model presented here is suitable to facilitate micron-scale studies with sub-nm resolution in many biological and engineering applications and can be used to investigate protein-DNA complexes, such as nucleosomes, or LLPS of proteins with DNA, enabling a mechanistic understanding of how molecular information may be propagated at the genome level.

Role of Strong Localized vs Weak Distributed Interactions in Disordered Protein Phase Separation

Interaction strength and localization are critical parameters controlling the single-chain and condensed-state properties of intrinsically disordered proteins (IDPs). Here, we decipher these relationships using coarse-grained heteropolymers comprised of hydrophobic (H) and polar (P) monomers as model IDPs. We systematically vary the fraction of P monomers XP and employ two distinct particle-based models that include either strong localized attractions between only H-H pairs (HP model) or weak distributed attractions between both H-H and H-P pairs (HP+ model). To compare different sequences and models, we first carefully tune the attraction strength for all sequences to match the single-chain radius of gyration. Interestingly, we find that this procedure produces similar conformational ensembles, nonbonded potential energies, and chain-level dynamics for single chains of almost all sequences in both models, with some deviations for the HP model at large XP. However, we observe a surprisingly rich phase behavior for the sequences in both models that deviates from the expectation that similarity at the single-chain level will translate to a similar phase-separation propensity. Coexistence between dilute and dense phases is only observed up to a model-dependent XP, despite the presence of favorable interchain interactions, which we quantify using the second virial coefficient. Instead, the limited number of attractive sites (H monomers) leads to the self-assembly of finite-sized clusters of different sizes depending on XP. Our findings strongly suggest that models with distributed interactions favor the formation of liquid-like condensates over a much larger range of sequence compositions compared to models with localized interactions.

Liquid-liquid phase separation: Galectin-3 in nuclear speckles and ribonucleoprotein complexes

Nuclear speckles are subcellular structures originally characterized by punctate immunofluorescence staining of the monoclonal antibody SC35, which recognizes an epitope on SRRM2 (serine/arginine repetitive matrix protein 2) and Sfrs2, a member of the SR (serine/arginine-rich) family of splicing factors. Galectin-3 co-localizes with SC35 in nuclear speckles, which represent one group of nuclear bodies that include the nucleolus, Cajal bodies and gems, paraspeckles, etc. Although they appear to have well-delineated physical boundaries, these nuclear bodies are not membrane-bound structures but represent macromolecular assemblies arising from a phenomenon called liquid-liquid phase separation. There has been much recent interest in liquid phase condensation as a newly recognized mechanism by which a cell can organize and compartmentalize subcellular structures with distinct composition. The punctate/speckled staining of galectin-3 with SC3 demonstrates their co-localization in a phase-separated body in vivo, under conditions endogenous to the cell. The purpose of the present review is to summarize the studies that document three key features of galectin-3 for its localization in liquid phase condensates: (a) an intrinsically disordered domain; (b) oligomer formation for multivalent binding; and (c) association with RNA and ribonucleoprotein complexes.

Advances in Nucleotide Repeat Expansion Diseases: Transcription Gets in Phase

Unstable DNA repeat expansions and insertions have been found to cause more than 50 neurodevelopmental, neurodegenerative, and neuromuscular disorders. One of the main hallmarks of repeat expansion diseases is the formation of abnormal RNA or protein aggregates in the neuronal cells of affected individuals. Recent evidence indicates that alterations of the dynamic or material properties of biomolecular condensates assembled by liquid/liquid phase separation are critical for the formation of these aggregates. This is a thermodynamically-driven and reversible local phenomenon that condenses macromolecules into liquid-like compartments responsible for compartmentalizing molecules required for vital cellular processes. Disease-associated repeat expansions modulate the phase separation properties of RNAs and proteins, interfering with the composition and/or the material properties of biomolecular condensates and resulting in the formation of abnormal aggregates. Since several repeat expansions have arisen in genes encoding crucial players in transcription, this raises the hypothesis that wide gene expression dysregulation is common to multiple repeat expansion diseases. This review will cover the impact of these mutations in the formation of aberrant aggregates and how they modify gene transcription.

A Zn-dependent structural transition of SOD1 modulates its ability to undergo phase separation

The misfolding and mutation of Cu/Zn superoxide dismutase (SOD1) is commonly associated with amyotrophic lateral sclerosis (ALS). SOD1 can accumulate within stress granules (SGs), a type of membraneless organelle, which is believed to form via liquid-liquid phase separation (LLPS). Using wild-type, metal-deficient, and different ALS disease mutants of SOD1 and computer simulations, we report here that the absence of Zn leads to structural disorder within two loop regions of SOD1, triggering SOD1 LLPS and amyloid formation. The addition of exogenous Zn to either metal-free SOD1 or to the severe ALS mutation I113T leads to the stabilization of the loops and impairs SOD1 LLPS and aggregation. Moreover, partial Zn-mediated inhibition of LLPS was observed for another severe ALS mutant, G85R, which shows perturbed Zn-binding. By contrast, the ALS mutant G37R, which shows reduced Cu-binding, does not undergo LLPS. In addition, SOD1 condensates induced by Zn-depletion exhibit greater cellular toxicity than aggregates formed by prolonged incubation under aggregating conditions. Overall, our work establishes a role for Zn-dependent modulation of SOD1 conformation and LLPS properties that may contribute to amyloid formation.

Heterotypic electrostatic interactions control complex phase separation of tau and prion into multiphasic condensates and co-aggregates

Biomolecular condensates formed via phase separation of proteins and nucleic acids are thought to perform a wide range of critical cellular functions by maintaining spatiotemporal regulation and organizing intracellular biochemistry. However, aberrant phase transitions are implicated in a multitude of human diseases. Here, we demonstrate that two neuronal proteins, namely tau and prion, undergo complex coacervation driven by domain-specific electrostatic interactions to yield highly dynamic, mesoscopic liquid-like droplets. The acidic N-terminal segment of tau interacts electrostatically with the polybasic N-terminal intrinsically disordered segment of the prion protein (PrP). We employed a unique combination of time-resolved tools that encompass several orders of magnitude of timescales ranging from nanoseconds to seconds. These studies unveil an intriguing symphony of molecular events associated with the formation of heterotypic condensates comprising ephemeral, domain-specific, short-range electrostatic nanoclusters. Our results reveal that these heterotypic condensates can be tuned by RNA in a stoichiometry-dependent manner resulting in reversible, multiphasic, immiscible, and ternary condensates of different morphologies ranging from core-shell to nested droplets. This ternary system exhibits a typical three-regime phase behavior reminiscent of other membraneless organelles including nucleolar condensates. We also show that upon aging, tau:PrP droplets gradually convert into solid-like co-assemblies by sequestration of persistent intermolecular interactions. Our vibrational Raman results in conjunction with atomic force microscopy and multi-color fluorescence imaging reveal the presence of amorphous and amyloid-like co-aggregates upon maturation. Our findings provide mechanistic underpinnings of overlapping neuropathology involving tau and PrP and highlight a broader biological role of complex phase transitions in physiology and disease.

Field-Theoretic Simulation Method to Study the Liquid-Liquid Phase Separation of Polymers

Liquid-liquid phase separation (LLPS) is a process that results in the formation of a polymer-rich liquid phase coexisting with a polymer-depleted liquid phase. LLPS plays a critical role in the cell through the formation of membrane-less organelles, but it also has a number of biotechnical and biomedical applications such as drug confinement and its targeted delivery. In this chapter, we present a computational efficient methodology that uses field-theoretic simulations (FTS) with complex Langevin (CL) sampling to characterize polymer phase behavior and delineate the LLPS phase boundaries. This approach is a powerful complement to analytical and explicit-particle simulations, and it can serve to inform experimental LLPS studies. The strength of the method lies in its ability to properly sample a large ensemble of polymers in a saturated solution while including the effect of composition fluctuations on LLPS. We describe the approaches that can be used to accurately construct phase diagrams of a variety of molecularly designed polymers and illustrate the method by generating an approximation-free phase diagram for a classical symmetric diblock polyampholyte.

Molecular interactions underlying the phase separation of HP1α: role of phosphorylation, ligand and nucleic acid binding

Heterochromatin protein 1α (HP1α) is a crucial element of chromatin organization. It has been proposed that HP1α functions through liquid-liquid phase separation (LLPS), which allows it to compact chromatin into transcriptionally repressed heterochromatin regions. In vitro, HP1α can undergo phase separation upon phosphorylation of its N-terminus extension (NTE) and/or through interactions with DNA and chromatin. Here, we combine computational and experimental approaches to elucidate the molecular interactions that drive these processes. In phosphorylation-driven LLPS, HP1α can exchange intradimer hinge-NTE interactions with interdimer contacts, which also leads to a structural change from a compacted to an extended HP1α dimer conformation. This process can be enhanced by the presence of positively charged HP1α peptide ligands and disrupted by the addition of negatively charged or neutral peptides. In DNA-driven LLPS, both positively and negatively charged peptide ligands can perturb phase separation. Our findings demonstrate the importance of electrostatic interactions in HP1α LLPS where binding partners can modulate the overall charge of the droplets and screen or enhance hinge region interactions through specific and non-specific effects. Our study illuminates the complex molecular framework that can fine-tune the properties of HP1α and that can contribute to heterochromatin regulation and function.

Thermostable Proteins from HaCaT Keratinocytes Identify a Wide Breadth of Intrinsically Disordered Proteins and Candidates for Liquid-Liquid Phase Separation

Intrinsically disordered proteins (IDPs) move through an ensemble of conformations which allows multitudinous roles within a cell. Keratinocytes, the predominant cell type in mammalian epidermis, have had only a few individual proteins assessed for intrinsic disorder and its possible contribution to liquid-liquid phase separation (LLPS), especially in regard to what functions or structures these proteins provide. We took a holistic approach to keratinocyte IDPs starting with enrichment via the isolation of thermostable proteins. The keratinocyte protein involucrin, known for its resistance to heat denaturation, served as a marker. It and other thermostable proteins were identified by liquid chromatography tandem mass spectrometry and subjected to extensive bioinformatic analysis covering gene ontology, intrinsic disorder, and potential for LLPS. Numerous proteins unique to keratinocytes and other proteins with shared expression in multiple cell types were identified to have IDP traits (e.g., compositional bias, nucleic acid binding, and repeat motifs). Among keratinocyte-specific proteins, many that co-assemble with involucrin into the cell-specific structure known as the cornified envelope scored highly for intrinsic disorder and potential for LLPS. This suggests intrinsic disorder and LLPS are previously unrecognized traits for assembly of the cornified envelope, echoing the contribution of intrinsic disorder and LLPS to more widely encountered features such as stress granules and PML bodies.

Analytical Formulation and Field-Theoretic Simulation of Sequence-Specific Phase Separation of Protein-Like Heteropolymers with Short- and Long-Spatial-Range Interactions

A theory for sequence-dependent liquid-liquid phase separation (LLPS) of intrinsically disordered proteins (IDPs) in the study of biomolecular condensates is formulated by extending the random phase approximation (RPA) and field-theoretic simulation (FTS) of heteropolymers with spatially long-range Coulomb interactions to include the fundamental effects of short-range, hydrophobic-like interactions between amino acid residues. To this end, short-range effects are modeled by Yukawa interactions between multiple nonelectrostatic charges derived from an eigenvalue decomposition of pairwise residue-residue contact energies. Chain excluded volume is afforded by incompressibility constraints. A mean-field approximation leads to an effective Flory-Huggins χ parameter, which, in conjunction with RPA, accounts for the contact-interaction effects of amino acid composition and the sequence-pattern effects of long-range electrostatics in IDP LLPS, whereas FTS based on the formulation provides full sequence dependence for both short- and long-range interactions. This general approach is illustrated here by applications to variants of a natural IDP in the context of several different amino-acid interaction schemes as well as a set of different model hydrophobic-polar sequences sharing the same composition. Effectiveness of the methodology is verified by coarse-grained explicit-chain molecular dynamics simulations.