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

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.

DNA-Mediated Formation of Phase-Separated Coacervates of the Nucleic Acid-Binding Domain of TAR DNA-Binding Protein (TDP-43) Prevents Its Amyloid-Like Misfolding

Sequestration of protein molecules and nucleic acids to stress granules is one of the most promising strategies that cells employ to protect themselves from stress. In vitro, studies suggest that the nucleic acid-binding domain of TDP-43 (TDP-43tRRM) undergoes amyloid-like aggregation to β-sheet-rich structures in low pH stress. In contrast, we observed that the TDP-43tRRM undergoes complex coacervation in the presence of ssDNA to a dense and light phase, preventing its amyloid-like aggregation. The soluble light phase consists of monomeric native-like TDP-43tRRM. The microscopic data suggest that the dense phase consists of spherical coacervates with limited internal dynamics. We performed multiparametric analysis by employing various biophysical techniques and found that complex coacervation depends on the concentration and ratio of the participating biomolecules and is driven by multivalent interactions. The modulation of these forces due to environmental conditions or disease mutations regulates the extent of coacervation, and the weakening of interactions between TDP-43tRRM and ssDNA leads to amyloid-like aggregation of TDP-43tRRM. Our results highlight a competition among the native state, amyloid-like aggregates, and complex coacervates tuned by various environmental factors. Together, our results illuminate an alternate function of TDP-43tRRM in response to pH stress in the presence of the ssDNA.

Pathological C-terminal phosphomimetic substitutions alter the mechanism of liquid-liquid phase separation of TDP-43 low complexity domain

C-terminally phosphorylated TAR DNA-binding protein of 43 kDa (TDP-43) marks the proteinaceous inclusions that characterize a number of age-related neurodegenerative diseases, including amyotrophic lateral sclerosis, frontotemporal lobar degeneration and Alzheimer's disease. TDP-43 phosphorylation at S403/S404 and (especially) at S409/S410 is, in fact, accepted as a biomarker of proteinopathy. These residues are located within the low complexity domain (LCD), which also drives the protein's liquid-liquid phase separation (LLPS). The impact of phosphorylation at these LCD sites on phase separation of the protein is a topic of great interest, as these post-translational modifications and LLPS are both implicated in proteinopathies. Here, we employed a combination of experimental and simulation-based approaches to explore this question on a phosphomimetic model of the TDP-43 LCD. Our turbidity and fluorescence microscopy data show that phosphomimetic Ser-to-Asp substitutions at residues S403, S404, S409 and S410 alter the LLPS behavior of TDP-43 LCD. In particular, unlike the LLPS of unmodified protein, LLPS of the phosphomimetic variants displays a biphasic dependence on salt concentration. Through coarse-grained modeling, we find that this biphasic salt dependence is derived from an altered mechanism of phase separation, in which LLPS-driving short-range intermolecular hydrophobic interactions are modulated by long-range attractive electrostatic interactions. Overall, this in vitro and in silico study provides a physiochemical foundation for understanding the impact of pathologically relevant C-terminal phosphorylation on the LLPS of TDP-43 in a more complex cellular environment.

Conformational Enigma of TDP-43 Misfolding in Neurodegenerative Disorders

Misfolding and aggregation of the protein remain some of the most common phenomena observed in neurodegeneration. While there exist multiple neurodegenerative disorders characterized by accumulation of distinct proteins, what remains particularly interesting is the ability of these proteins to undergo a conformational change to form aggregates. TDP-43 is one such nucleic acid binding protein whose misfolding is associated with many neurogenerative diseases including amyotrophic lateral sclerosis (ALS) and fronto-temporal lobar degeneration (FTLD). TDP-43 protein assumes several different conformations and oligomeric states under the diseased condition. In this review, we explore the intrinsic relationship between the conformational variability of TDP-43 protein, with a particular focus on the RRM domains, and its propensity to undergo aggregation. We further emphasize the probable mechanism behind the formation of these conformations and suggest a potential diagnostic and therapeutic strategy in the context of these conformational states of the protein.

Rational Tuning of the Concentration-independent Enrichment of Prion-like Domains in Stress Granules

Stress granules (SGs) are large ribonucleoprotein assemblies that form in response to acute stress in eukaryotes. SG formation is thought to be initiated by liquid-liquid phase separation (LLPS) of key proteins and RNA. These molecules serve as a scaffold for recruitment of client molecules. LLPS of scaffold proteins in vitro is highly concentration-dependent, yet biomolecular condensates in vivo contain hundreds of unique proteins, most of which are thought to be clients rather than scaffolds. Many proteins that localize to SGs contain low-complexity, prion-like domains (PrLDs) that have been implicated in LLPS and SG recruitment. The degree of enrichment of proteins in biomolecular condensates such as SGs can vary widely, but the underlying basis for these differences is not fully understood. Here, we develop a toolkit of model PrLDs to examine the factors that govern efficiency of PrLD recruitment to stress granules. Recruitment was highly sensitive to amino acid composition: enrichment in SGs could be tuned through subtle changes in hydrophobicity. By contrast, SG recruitment was largely insensitive to PrLD concentration at both a population level and single-cell level. These observations point to a model wherein PrLDs are enriched in SGs through either simple solvation effects or interactions that are effectively non-saturable even at high expression levels.

Differential Effects of Sequence-Local versus Nonlocal Charge Patterns on Phase Separation and Conformational Dimensions of Polyampholytes as Model Intrinsically Disordered Proteins

Conformational properties of intrinsically disordered proteins (IDPs) are governed by a sequence-ensemble relationship. To differentiate the impact of sequence-local versus sequence-nonlocal features of an IDP's charge pattern on its conformational dimensions and its phase-separation propensity, the charge "blockiness" κ and the nonlocality-weighted sequence charge decoration (SCD) parameters are compared for their correlations with isolated-chain radii of gyration (Rgs) and upper critical solution temperatures (UCSTs) of polyampholytes modeled by random phase approximation, field-theoretic simulation, and coarse-grained molecular dynamics. SCD is superior to κ in predicting Rg because SCD accounts for effects of contact order, i.e., nonlocality, on dimensions of isolated chains. In contrast, κ and SCD are comparably good, though nonideal, predictors of UCST because frequencies of interchain contacts in the multiple-chain condensed phase are less sensitive to sequence positions than frequencies of intrachain contacts of an isolated chain, as reflected by κ correlating better with condensed-phase interaction energy than SCD.

Multiscale simulations reveal the driving forces of p53C phase separation accelerated by oncogenic mutations

Liquid-Liquid phase separation (LLPS) of p53 to form liquid condensates has been implicated in cellular functions and dysfunctions. The p53 condensates may serve as amyloid fibril precursors to initiate p53 aggregation, which is associated with oncogenic gain-of-function and various human cancers. M237I and R249S mutations located in p53 core domain (p53C) have been detected respectively in glioblastomas and hepatocellular carcinoma. Interestingly, these p53C mutants can also undergo LLPS and liquid-to-solid phase transition, which are faster than wild type p53C. However, the underlying molecular basis governing the accelerated LLPS and liquid-to-solid transition of p53C remain poorly understood. Herein, we explore the M237I/R249S mutation-induced structural alterations and phase separation behavior of p53C by employing multiscale molecular dynamics simulations. All-atom simulations revealed conformational disruptions in the zinc-binding domain of the M237I mutant and in both loop3 and zinc-binding domain of the R249S mutant. The two mutations enhance hydrophobic exposure of those regions and attenuate intramolecular interactions, which may hasten the LLPS and aggregation of p53C. Martini 3 coarse-grained simulations demonstrated spontaneous phase separation of p53C and accelerated effects of M237I/R249S mutations on the phase separation of p53C. Importantly, we find that the regions with enhanced intermolecular interactions observed in coarse-grained simulations coincide with the disrupted regions with weakened intramolecular interactions observed in all-atom simulations, indicating that M237I/R249S mutation-induced local structural disruptions expedite the LLPS of p53C. This study unveils the molecular mechanisms underlying the two cancer-associated mutation-accelerated LLPS and aggregation of p53C, providing avenues for anticancer therapy by targeting the phase separation process.

Context-dependent, fuzzy protein interactions: Towards sequence-based insights

Predicting protein interactions in the cellular environment still remains a challenge in the AlphaFold era. Protein interactions, similarly to their structures, sample a continuum from ordered to disordered states, with specific partners in many bound configurations. A multiplicity of binding modes (MBM) enables transition between these states under different cellular conditions. This review focuses on how the cellular environment affects protein interactions, highlighting the molecular mechanisms, biophysical origin, and sequence-based principles of context-dependent, fuzzy interactions. It summarises experimental and computational approaches to address the challenge of interaction heterogeneity and its contribution to a wide range of biological functions. These insights will help in understanding complex cellular processes, involving conversions between protein assembly states, such as from liquid-like droplet state to the amyloid state.

Fundamental Aspects of Phase-Separated Biomolecular Condensates

Biomolecular condensates, formed through phase separation, are upending our understanding in much of molecular, cell, and developmental biology. There is an urgent need to elucidate the physicochemical foundations of the behaviors and properties of biomolecular condensates. Here we aim to fill this need by writing a comprehensive, critical, and accessible review on the fundamental aspects of phase-separated biomolecular condensates. We introduce the relevant theoretical background, present the theoretical basis for the computation and experimental measurement of condensate properties, and give mechanistic interpretations of condensate behaviors and properties in terms of interactions at the molecular and residue levels.

Structural details of helix-mediated TDP-43 C-terminal domain multimerization

The primarily disordered C-terminal domain (CTD) of TAR DNA binding protein-43 (TDP-43), a key nuclear protein in RNA metabolism, forms neuronal inclusions in several neurodegenerative diseases. A conserved region (CR, spanning residues 319-341) in CTD forms transient helix-helix contacts important for its higher-order oligomerization and function that are disrupted by ALS-associated mutations. However, the structural details of CR assembly and the explanation for several ALS-associated variants' impact on phase separation and function remain unclear due to challenges in analyzing the dynamic association of TDP-43 CTD using traditional structural biology approaches. By employing an integrative approach, combining biophysical experiments, biochemical assays, AlphaFold2-Multimer (AF2-Multimer), and atomistic simulations, we generated structural models of helical oligomerization of TDP-43 CR. Using NMR, we first established that the native state of TDP-43 CR under physiological conditions is α-helical. Next, alanine scanning mutagenesis revealed that while hydrophobic residues in the CR are important for CR assembly, phase separation and TDP-43 nuclear retention function, polar residues down regulate these processes. Finally, pairing AF2-Multimer modeling with AAMD simulations indicated that dynamic, oligomeric assemblies of TDP-43 that are stabilized by a methionine-rich core with specific contributions from a tryptophan/leucine pair. In conclusion, our results advance the structural understanding of the mechanisms driving TDP-43 function and provide a window into the initial stages of its conversion into pathogenic aggregates.

Structure and interactions of prion-like domains in transcription factor Efg1 phase separation

Candida albicans, a prominent member of the human microbiome, can make an opportunistic switch from commensal coexistence to pathogenicity accompanied by an epigenetic shift between the white and opaque cell states. This transcriptional switch is under precise regulation by a set of transcription factors (TFs), with Enhanced Filamentous Growth Protein 1 (Efg1) playing a central role. Previous research has emphasized the importance of Efg1's prion-like domain (PrLD) and the protein's ability to undergo phase separation for the white-to-opaque transition of C. albicans. However, the underlying molecular mechanisms of Efg1 phase separation have remained underexplored. In this study, we delved into the biophysical basis of Efg1 phase separation, revealing the significant contribution of both N-terminal (N) and C-terminal (C) PrLDs. Through NMR structural analysis, we found that Efg1 N-PrLD and C-PrLD are mostly disordered but have prominent partial α-helical secondary structures in both domains. NMR titration experiments suggest that the partially helical structures in N-PrLD act as hubs for self-interaction as well as Efg1 interaction with RNA. Using condensed-phase NMR spectroscopy, we uncovered diverse amino acid interactions underlying Efg1 phase separation. Particularly, we highlight the indispensable role of tyrosine residues within the transient α-helical structures of PrLDs particularly in the N-PrLD compared to the C-PrLD in stabilizing phase separation. Our study provides evidence that the transient α-helical structure is present in the phase-separated state and highlights the particular importance of aromatic residues within these structures for phase separation. Together, these results enhance the understanding of C. albicans transcription factor interactions that lead to virulence and provide a crucial foundation for potential antifungal therapies targeting the transcriptional switch.

Transient interdomain interactions modulate the monomeric structural ensemble and oligomerization landscape of Huntingtin Exon 1

Polyglutamine expansion (≥ 36 residues) within the N-terminal exon-1 of Huntingtin (Httex1) leads to Huntington's disease, a neurogenerative condition marked by the presence of intranuclear Htt inclusions. Notably, the polyglutamine tract in Httex1 is flanked by an N-terminal coiled-coil domain - N17 (17 amino acids), which undergoes self-association to promote the formation of soluble Httex1 oligomers and brings the aggregation-prone polyQ tracts in close spatial proximity. However, the mechanisms underlying the subsequent conversion of soluble oligomers into insoluble β-rich aggregates with increasing polyQ length, remain unclear. Current knowledge suggests that expansion of the polyQ tract increases its helicity, and this favors its oligomerization and aggregation. In addition, studies utilizing conformation-specific antibodies and a stable coiled-coil heterotetrametric system fused to polyQ indicate that domain "cross-talk" (i.e., interdomain interactions) may be necessary to efficiently promote the emergence of toxic conformations (in monomers and oligomers) and fibrillar aggregation. Here, we performed extensive atomistic molecular dynamics (MD) simulations (aggregate time ∼ 0.7 ms) of N17-polyQ fragments to uncover the interplay between structural transformation and domain "cross-talk" on the monomeric structural ensemble and oligomerization landscape of Httex1. Our simulation ensembles of N17-polyQ monomers validated against 13 C NMR chemical shifts indicated that in addition to elevated α-helicity, polyQ expansion also favors transient, interdomain (N17-polyQ) interactions which result in the emergence of β-conformations. Further, interdomain interactions decreased the overall stability of N17-mediated dimers by counteracting the stabilizing effect of increased α-helicity and promoted a heterogenous oligomerization landscape on the sub-microsecond timescale. Overall, our study uncovers the significance of domain "cross-talk" in modulating the monomeric conformational ensemble and oligomerization landscape of Httex1 to favor the formation of amyloid aggregates.

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.

Co-Assembly of Cancer Drugs with Cyclo-HH Peptides: Insights from Simulations and Experiments

Peptide-based nanomaterials can serve as promising drug delivery agents, facilitating the release of active pharmaceutical ingredients while reducing the risk of adverse reactions. We previously demonstrated that Cyclo-Histidine-Histidine (Cyclo-HH), co-assembled with cancer drug Epirubicin, zinc, and nitrate ions, can constitute an attractive drug delivery system, combining drug self-encapsulation, enhanced fluorescence, and the ability to transport the drug into cells. Here, we investigated both computationally and experimentally whether Cyclo-HH could co-assemble, in the presence of zinc and nitrate ions, with other cancer drugs with different physicochemical properties. Our studies indicated that Methotrexate, in addition to Epirubicin and its epimer Doxorubicin, and to a lesser extent Mitomycin-C and 5-Fluorouracil, have the capacity to co-assemble with Cyclo-HH, zinc, and nitrate ions, while a significantly lower propensity was observed for Cisplatin. Epirubicin, Doxorubicin, and Methorexate showed improved drug encapsulation and drug release properties, compared to Mitomycin-C and 5-Fluorouracil. We demonstrated the biocompatibility of the co-assembled systems, as well as their ability to intracellularly release the drugs, particularly for Epirubicin, Doxorubicin, and Methorexate. Zinc and nitrate were shown to be important in the co-assembly, coordinating with drugs and/or Cyclo-HH, thereby enabling drug-peptide as well as drug-drug interactions in successfully formed nanocarriers. The insights could be used in the future design of advanced cancer therapeutic systems with improved properties.

Pathological C-terminal phosphomimetic substitutions alter the mechanism of liquid-liquid phase separation of TDP-43 low complexity domain

C-terminally phosphorylated TAR DNA-binding protein of 43 kDa (TDP-43) marks the proteinaceous inclusions that characterize a number of age-related neurodegenerative diseases, including amyotrophic lateral sclerosis, frontotemporal lobar degeneration and Alzheimer's disease. TDP-43 phosphorylation at S403/S404, and especially at S409/S410, is in fact accepted as a biomarker of proteinopathy. These residues are located within the low complexity domain (LCD), which also drives the protein's liquid-liquid phase separation (LLPS). The impact of phosphorylation at these LCD sites on phase separation of the protein is a topic of great interest, as these post-translational modifications and LLPS are both implicated in proteinopathies. Here, we employed a combination of experimental and simulation-based approaches to explore this question on a phosphomimetic model of the TDP-43 LCD. Our turbidity and fluorescence microscopy data show that Ser-to-Asp substitutions at residues S403, S404, S409 and S410 alter the LLPS behavior of TDP-43 LCD. In particular, in contrast to the unmodified protein, the phosphomimetic variants display a biphasic dependence on salt concentration. Through coarse-grained modeling, we find that this biphasic salt dependence is derived from an altered mechanism of phase separation, in which LLPS-driving short-range intermolecular hydrophobic interactions are modulated by long-range attractive electrostatic interactions. Overall, this in vitro and in silico study provides a physiochemical foundation for understanding the impact of pathologically-relevant C-terminal phosphorylation on the LLPS of the TDP-43 in a more complex cellular environment.

Expanding the molecular grammar of polar residues and arginine in FUS prion-like domain phase separation and aggregation

A molecular grammar governing low-complexity prion-like domains phase separation (PS) has been proposed based on mutagenesis experiments that identified tyrosine and arginine as primary drivers of phase separation via aromatic-aromatic and aromatic-arginine interactions. Here we show that additional residues make direct favorable contacts that contribute to phase separation, highlighting the need to account for these contributions in PS theories and models. We find that tyrosine and arginine make important contacts beyond only tyrosine-tyrosine and tyrosine-arginine, including arginine-arginine contacts. Among polar residues, glutamine in particular contributes to phase separation with sequence/position-specificity, making contacts with both tyrosine and arginine as well as other residues, both before phase separation and in condensed phases. For glycine, its flexibility, not its small solvation volume, favors phase separation by allowing favorable contacts between other residues and inhibits the liquid-to-solid (LST) transition. Polar residue types also make sequence-specific contributions to aggregation that go beyond simple rules, which for serine positions is linked to formation of an amyloid-core structure by the FUS low-complexity domain. Hence, here we propose a revised molecular grammar expanding the role of arginine and polar residues in prion-like domain protein phase separation and aggregation.

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.

Toward Accurate Simulation of Coupling between Protein Secondary Structure and Phase Separation

Intrinsically disordered proteins (IDPs) frequently mediate phase separation that underlies the formation of a biomolecular condensate. Together with theory and experiment, efficient coarse-grained (CG) simulations have been instrumental in understanding the sequence-specific phase separation of IDPs. However, the widely used Cα-only models are limited in capturing the peptide nature of IDPs, particularly backbone-mediated interactions and effects of secondary structures, in phase separation. Here, we describe a hybrid resolution (HyRes) protein model toward a more accurate description of the backbone and transient secondary structures in phase separation. With an atomistic backbone and coarse-grained side chains, HyRes can semiquantitatively capture the residue helical propensity and overall chain dimension of monomeric IDPs. Using GY-23 as a model system, we show that HyRes is efficient enough for the direct simulation of spontaneous phase separation and, at the same time, appears accurate enough to resolve the effects of single His to Lys mutations. HyRes simulations also successfully predict increased β-structure formation in the condensate, consistent with available experimental CD data. We further utilize HyRes to study the phase separation of TPD-43, where several disease-related mutants in the conserved region (CR) have been shown to affect residual helicities and modulate the phase separation propensity as measured by the saturation concentration. The simulations successfully recapitulate the effect of these mutants on the helicity and phase separation propensity of TDP-43 CR. Analyses reveal that the balance between backbone and side chain-mediated interactions, but not helicity itself, actually determines phase separation propensity. These results support that HyRes represents an effective protein model for molecular simulation of IDP phase separation and will help to elucidate the coupling between transient secondary structures and phase separation.

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.

Structure and position-specific interactions of prion-like domains in transcription factor Efg1 phase separation

Candida albicans, a prominent member of the human microbiome, can make an opportunistic switch from commensal coexistence to pathogenicity accompanied by an epigenetic shift between the white and opaque cell states. This transcriptional switch is under precise regulation by a set of transcription factors (TFs), with Enhanced Filamentous Growth Protein 1 (Efg1) playing a central role. Previous research has emphasized the importance of Egf1's prion-like domain (PrLD) and the protein's ability to undergo phase separation for the white-to-opaque transition of C. albicans. However, the underlying molecular mechanisms of Efg1 phase separation have remained underexplored. In this study, we delved into the biophysical basis of Efg1 phase separation, revealing the significant contribution of both N-terminal (N) and C-terminal (C) PrLDs. Through NMR structural analysis, we found that Efg1 N-PrLD and C-PrLD are mostly disordered though have prominent partial α-helical secondary structures in both domains. NMR titration experiments suggest that the partially helical structures in N-PrLD act as hubs for self-interaction as well as Efg1 interaction with RNA. Using condensed-phase NMR spectroscopy, we uncovered diverse amino acid interactions underlying Efg1 phase separation. Particularly, we highlight the indispensable role of tyrosine residues within the transient α-helical structures of PrLDs particularly in the N-PrLD compared to the C-PrLD in stabilizing phase separation. Our study provides evidence that the transient α-helical structure is present in the phase separated state and highlights the particular importance of aromatic residues within these structures for phase separation. Together, these results enhance the understanding of C. albicans TF interactions that lead to virulence and provide a crucial foundation for potential antifungal therapies targeting the transcriptional switch.

Accurate Simulation of Coupling between Protein Secondary Structure and Liquid-Liquid Phase Separation

Intrinsically disordered proteins (IDPs) frequently mediate liquid-liquid phase separation (LLPS) that underlies the formation of membraneless organelles. Together with theory and experiment, efficient coarse-grained (CG) simulations have been instrumental in understanding sequence-specific phase separation of IDPs. However, the widely-used Cα-only models are severely limited in capturing the peptide nature of IDPs, including backbone-mediated interactions and effects of secondary structures, in LLPS. Here, we describe a hybrid resolution (HyRes) protein model for accurate description of the backbone and transient secondary structures in LLPS. With an atomistic backbone and coarse-grained side chains, HyRes accurately predicts the residue helical propensity and chain dimension of monomeric IDPs. Using GY-23 as a model system, we show that HyRes is efficient enough for direct simulation of spontaneous phase separation, and at the same time accurate enough to resolve the effects of single mutations. HyRes simulations also successfully predict increased beta-sheet formation in the condensate, consistent with available experimental data. We further utilize HyRes to study the phase separation of TPD-43, where several disease-related mutants in the conserved region (CR) have been shown to affect residual helicities and modulate LLPS propensity. The simulations successfully recapitulate the effect of these mutants on the helicity and LLPS propensity of TDP-43 CR. Analyses reveal that the balance between backbone and sidechain-mediated interactions, but not helicity itself, actually determines LLPS propensity. We believe that the HyRes model represents an important advance in the molecular simulation of LLPS and will help elucidate the coupling between IDP transient secondary structures and phase separation.