Phosphorylation sites are evolutionary checkpoints against liquid-solid transition in protein condensates

Molecular simulations of enzymatic phosphorylation of disordered proteins and their condensates

Condensation and aggregation of disordered proteins in cellular non-equilibrium environments are shaped decisively by enzymes. Enzymes called kinases phosphorylate proteins, consuming the chemical fuel ATP. Protein phosphorylation by kinases such as Casein kinase 1 delta (CK1δ) determines the interactions of neurodegeneration-linked proteins such as TDP-43. Hyperphosphorylation of TDP-43 by CK1δ may be a cytoprotective mechanism for neurons, but how CK1δ interacts with protein condensates is not known. Molecular dynamics simulations hold the promise to resolve how kinases interact with disordered proteins and their condensates, and how this shapes the phosphorylation dynamics. In practice, it is difficult to verify whether implementations of chemical-fuel driven coarse-grained simulations are thermodynamically consistent, which we address by a generally applicable and automatic Markov state modeling approach. In this work, we thus elucidate with coarse-grained simulations, drivers of how TDP-43 is phosphorylated by CK1δ and how this leads to the dissolution of TDP-43 condensates upon hyperphosphorylation.

Emerging connections: Poly(ADP-ribose), FET proteins and RNA in the regulation of DNA damage condensates

Our genome is exposed to thousands of DNA lesions every day, posing a significant threat to cellular viability. To deal with these lesions, cells have evolved sophisticated repair mechanisms collectively known as the DNA damage response. DNA double-strand breaks (DSBs) are very cytotoxic damages, and their repair requires the precise and coordinated recruitment of multiple repair factors to form nuclear foci. Recent research highlighted that these repair structures behave as biomolecular condensates, i.e. membraneless compartments with liquid-like properties. The formation of condensates is driven by weak, multivalent interactions among proteins and nucleic acids, and recent studies highlighted the roles of poly(ADP-ribose) (PAR) and RNA in regulating DSBs-related condensates. Additionally, the FET family of RNA-binding proteins (including FUS, EWS and TAF15), has emerged as a critical player in the DNA damage response, with recent evidence suggesting that FET proteins support the formation and dynamics of repair condensates. Notably, phase separation of FET proteins is implicated also in their pathological functions in cancer biology, highlighting the pervasive role of condensation. This review will provide an overview of biomolecular condensates at DSBs, focusing on the interplay among PAR and RNA in the spatiotemporal regulation of FET proteins at repair complexes. We will also discuss the role of FET condensates in cancer biology and how they are targeted for therapeutic purposes. The study of biomolecular condensates holds great promise for advancing our understanding of key cellular processes and developing novel therapeutic strategies, but requires careful consideration of potential challenges.

Semaphorin 6A phase separation sustains a histone lactylation-dependent lactate buildup in pathological angiogenesis

Ischemic retinal diseases are major causes of blindness worldwide and are characterized by pathological angiogenesis. Epigenetic alterations in response to metabolic shifts in endothelial cells (ECs) suffice to underlie excessive angiogenesis. Lactate accumulation and its subsequent histone lactylation in ECs contribute to vascular disorders. However, the regulatory mechanism of establishing and sustaining lactylation modification remains elusive. Here, we showed that lactate accumulation induced histone lactylations on H3K9 and H3K18 in neovascular ECs in the proliferative stage of oxygen-induced retinopathy. Joint CUT&Tag and scRNA-seq analyses identified Prmt5 as a target of H3K9la and H3K18la in isolated retinal ECs. EC-specific deletion of Prmt5 since the early stage of revascularization suppressed a positive feedback loop of lactate production and histone lactylation, thus inhibiting neovascular tuft formation. Mechanistically, the C-terminal intrinsically disorder region (IDR) of the transmembrane semaphorin 6A (SEMA6A) forms liquid-liquid phase separation condensates to recruit RHOA and P300, facilitating P300 phosphorylation and histone lactylation cycle. Deletion of endothelial Sema6A reduced H3K9la and H3K18la at the promoter of PRMT5 and diminished its expression. The induction of histone lactylation by SEMA6A-IDR and its pro-angiogenic effect were abrogated by deletion of Prmt5. Our study illustrates a sustainable histone lactylation machinery driven by phase separation-dependent lactyltransferase activation in dysregulated vascularization.

The Origin of the Ionic-strength Dependent Reentrant Behavior in Liquid-Liquid Phase Separation of Neutral IDPs

The effect of salt on coacervation of synthetic or biological polyelectrolytes is well-studied. However, recent experiments showed that largely neutral IDPs like FUS also undergo LLPS at physiological salt concentrations, dissolve at higher salt concentration and again phase separate at higher salt concentrations such as, [C ion ]∼3M. Here we use analytical theory and simulations to reveal the mechanism of these transitions. At low [C ion ], the ionic solution acts as a highly correlated medium conferring long-range effective attractive interactions between spatially distant FUS monomers. In this regime the ion concentration inside the condensate is higher than in the bulk solution. As [C ion ] increases, the correlation length in the ionic plasma decreases, and the condensate dissolves. Second LLPS at high [C ion ] is due to the entropy-driven crowding, and ion concentration inside the condensate is lower than in the bulk. Our study unravels a general physical mechanism of salt-dependent reentrant behavior in LLPS in neutral IDPs.

Phase Separation Clustering of Poly Ubiquitin Cargos on Ternary Mixture Lipid Membranes by Synthetically Cross-Linked Ubiquitin Binder Peptides

Ubiquitylation is involved in various physiological processes, such as signaling and vesicle trafficking, whereas ubiquitin (UB) is considered an important clinical target. The polymeric addition of UB enables cargo molecules to be recognized specifically by multivalent binding interactions with UB-binding proteins, which results in various downstream processes. Recently, protein condensate formation by ubiquitylated proteins has been reported in many independent UB processes, suggesting its potential role in governing the spatial organization of ubiquitylated cargo proteins. We created modular polymeric UB binding motifs and polymeric UB cargos by synthetic bioconjugation and protein purification. Giant unilamellar vesicles with lipid raft composition were prepared to reconstitute the polymeric UB cargo organization on the membranes. Fluorescence imaging was used to observe the outcome. The polymeric UB cargos clustered on the membranes by forming a phase separation codomain during the interaction with the multivalent UB-binding conjugate. This phase separation was valence-dependent and strongly correlated with its potent ability to form protein condensate droplets in solution. Multivalent UB binding interactions exhibited a general trend toward the formation of phase-separated condensates and the resulting condensates were either in a liquid-like or solid-like state depending on the conditions and interactions. This suggests that the polymeric UB cargos on the plasma and endosomal membranes may use codomain phase separation to assist in the clustering of UB cargos on the membranes for cargo sorting. Our findings also indicate that such phase behavior model systems can be created by a modular synthetic approach that can potentially be used to further engineer biomimetic interactions in vitro.

Expanding the molecular grammar of polar residues and arginine in FUS phase separation

A molecular grammar governing low-complexity prion-like domain phase separation (PS) has identified tyrosine and arginine as primary drivers via aromatic-aromatic and aromatic-arginine interactions. Here we show that additional residues and contacts contribute to PS, highlighting the need to include these contributions in PS theories and models. Tyrosine and arginine make important contacts beyond tyrosine-tyrosine and tyrosine-arginine, including arginine-arginine contacts. Among polar residues, glutamine contributes to PS with sequence and position specificity, contacting tyrosine, arginine and other residues, both before PS and in condensed phases. The flexibility of glycine enhances PS by allowing favorable contacts between adjacent residues and inhibits the liquid-to-solid transition. Polar residues also make sequence-specific contributions to liquid-to-solid transition, with serine positions linked to the formation of an amyloid-core structure by the FUS low-complexity domain. Hence, an extended molecular grammar expands the role of arginine and polar residues in prion-like domain protein PS and reveals the position dependence of residue contribution to solidification.

Separation of sticker-spacer energetics governs the coalescence of metastable condensates

Biological condensates often emerge as a multidroplet state and never coalesce into one large droplet within the experimental timespan. Previous work revealed that the sticker-spacer architecture of biopolymers may dynamically stabilize the multidroplet state. Here, we simulate the condensate coalescence using metadynamics approach and reveal two distinct physical mechanisms underlying the fusion of droplets. Condensates made of sticker-spacer polymers readily undergo a kinetic arrest when stickers exhibit slow exchange while fast exchanging stickers at similar levels of saturation allow merger to equilibrium states. On the other hand, condensates composed of homopolymers fuse readily until they reach a threshold density. Increase in entropy upon intercondensate mixing of chains drives the fusion of sticker-spacer chains. We map the range of mechanisms of kinetic arrest from slow sticker exchange dynamics to density mediated in terms of energetic separation of stickers and spacers. Our predictions appear to be in qualitative agreement with recent experiments probing dynamic nature of protein-RNA condensates.

Recent Progress in Modeling and Simulation of Biomolecular Crowding and Condensation Inside Cells

Macromolecular crowding in the cellular cytoplasm can potentially impact diffusion rates of proteins, their intrinsic structural stability, binding of proteins to their corresponding partners as well as biomolecular organization and phase separation. While such intracellular crowding can have a large impact on biomolecular structure and function, the molecular mechanisms and driving forces that determine the effect of crowding on dynamics and conformations of macromolecules are so far not well understood. At a molecular level, computational methods can provide a unique lens to investigate the effect of macromolecular crowding on biomolecular behavior, providing us with a resolution that is challenging to reach with experimental techniques alone. In this review, we focus on the various physics-based and data-driven computational methods developed in the past few years to investigate macromolecular crowding and intracellular protein condensation. We review recent progress in modeling and simulation of biomolecular systems of varying sizes, ranging from single protein molecules to the entire cellular cytoplasm. We further discuss the effects of macromolecular crowding on different phenomena, such as diffusion, protein-ligand binding, and mechanical and viscoelastic properties, such as surface tension of condensates. Finally, we discuss some of the outstanding challenges that we anticipate the community addressing in the next few years in order to investigate biological phenomena in model cellular environments by reproducing in vivo conditions as accurately as possible.

Sequence complexity and monomer rigidity control the morphologies and aging dynamics of protein aggregates

Understanding the biophysical basis of protein aggregation is important in biology because of the potential link to several misfolding diseases. Although experiments have shown that protein aggregates adopt a variety of morphologies, the dynamics of their formation are less well characterized. Here, we introduce a minimal model to explore the dependence of the aggregation dynamics on the structural and sequence features of the monomers. Using simulations, we demonstrate that sequence complexity (codified in terms of word entropy) and monomer rigidity profoundly influence the dynamics and morphology of the aggregates. Flexible monomers with low sequence complexity (corresponding to repeat sequences) form liquid-like droplets that exhibit ergodic behavior. Strikingly, these aggregates abruptly transition to more ordered structures, reminiscent of amyloid fibrils, when the monomer rigidity is increased. In contrast, aggregates resulting from monomers with high sequence complexity are amorphous and display nonergodic glassy dynamics. The heterogeneous dynamics of the low and high-complexity sequences follow stretched exponential kinetics, which is one of the characteristics of glassy dynamics. Importantly, at nonzero values of the bending rigidities, the aggregates age with the relaxation times that increase with the waiting time. Informed by these findings, we provide insights into aging dynamics in protein condensates and contrast the behavior with the dynamics expected in RNA repeat sequences. Our findings underscore the influence of the monomer characteristics in shaping the morphology and dynamics of protein aggregates, thus providing a foundation for deciphering the general rules governing the behavior of protein condensates.

Separation of sticker-spacer energetics governs the coalescence of metastable condensates

Biological condensates often emerge as a multi-droplet state and never coalesce into one large droplet within the experimental timespan. Previous work revealed that the sticker-spacer architecture of biopolymers may dynamically stabilize the multi-droplet state. Here, we simulate the condensate coalescence using metadynamics approach and reveal two distinct physical mechanisms underlying the fusion of droplets. Condensates made of sticker-spacer polymers readily undergo a kinetic arrest when stickers exhibit slow exchange while fast exchanging stickers at similar levels of saturation allow merger to equilibrium states. On the other hand, condensates composed of homopolymers fuse readily until they reach a threshold density. Increase in entropy upon inter-condensate mixing of chains drives the fusion of sticker-spacer chains. We map the range of mechanisms of kinetic arrest from slow sticker exchange dynamics to density mediated in terms of energetic separation of stickers and spacers. Our predictions appear to be in qualitative agreement with recent experiments probing dynamic nature of protein-RNA condensates.

IGF2BP1 phosphorylation in the disordered linkers regulates ribonucleoprotein condensate formation and RNA metabolism

The insulin-like growth factor 2 mRNA binding protein 1 (IGF2BP1) is a conserved RNA-binding protein that regulates RNA stability, localization and translation. IGF2BP1 is part of various ribonucleoprotein (RNP) condensates. However, the mechanism that regulates its assembly into condensates remains unknown. By using proteomics, we demonstrate that phosphorylation of IGF2BP1 at S181 in a disordered linker is regulated in a stress-dependent manner. Phosphomimetic mutations in two disordered linkers, S181E and Y396E, modulate RNP condensate formation by IGF2BP1 without impacting its binding affinity for RNA. Intriguingly, the S181E mutant, which lies in linker 1, impairs IGF2BP1 condensate formation in vitro and in cells, whereas a Y396E mutant in the second linker increases condensate size and dynamics. Structural approaches show that the first linker binds RNAs nonspecifically through its RGG/RG motif, an interaction weakened in the S181E mutant. Notably, linker 2 interacts with IGF2BP1's folded domains and these interactions are partially impaired in the Y396E mutant. Importantly, the phosphomimetic mutants impact IGF2BP1's interaction with RNAs and remodel the transcriptome in cells. Our data reveal how phosphorylation modulates low-affinity interaction networks in disordered linkers to regulate RNP condensate formation and RNA metabolism.

Wet Conformation of Prion-Like Domain and Intimate Correlation of Hydration and Conformational Fluctuations

Proteins with prion-like domains (PLDs) are involved in neurodegeneration-associated aggregation and are prevalent in liquid-like membrane-less organelles. These PLDs contain amyloidogenic stretches but can maintain dynamic disordered conformations, even in the condensed phase. However, the molecular mechanism underlying such intricate conformational properties of PLDs remains elusive. Here we employed molecular dynamics simulations to investigate the conformational properties of a prototypical PLD system (i.e., FUS PLD). According to our simulation results, PLD adopts a wet collapsed conformation, wherein most residues maintain sufficient hydration with the abundance of internal water. These internal water molecules can rapidly exchange between the protein interior and the bulk, enabling intensive coupling of the entire protein with its hydration environment. The dynamic exchange of water molecules is intimately correlated to the overall conformational fluctuations of PLD. Furthermore, the abundance of dynamic internal water suppresses the formation of aggregation-prone ordered structures. These results collectively elucidate the crucial role of internal water in sustaining the dynamic disordered conformation of the PLD and inhibiting its aggregation propensity.

Post-translational modification sites are present in hydrophilic cavities of alpha-synuclein, tau, FUS, and TDP-43 fibrils: A molecular dynamics study

Hydration plays a crucial role in the refolding of intrinsically disordered proteins into amyloid fibrils; however, the specific interactions between water and protein that may contribute to this process are still unknown. In our previous studies of alpha-synuclein (aSyn), we have shown that waters confined in fibril cavities are stabilizing features of this pathological fold; and that amino acids that hydrogen bond with these confined waters modulate primary and seeded aggregation. Here, we extend our aSyn molecular dynamics (MD) simulations with three new polymorphs and correlate MD trajectory information with known post-translational modifications (PTMs) and experimental data. We show that cavity residues are more evolutionarily conserved than non-cavity residues and are enriched with PTM sites. As expected, the confinement within hydrophilic cavities results in more stably hydrated amino acids. Interestingly, cavity PTM sites display the longest protein-water hydrogen bond lifetimes, three-fold greater than non-PTM cavity sites. Utilizing the deep mutational screen dataset by Newberry et al. and the Thioflavin T aggregation review by Pancoe et al. parsed using a fibril cavity/non-cavity definition, we show that hydrophobic changes to amino acids in cavities have a larger effect on fitness and aggregation rate than residues outside cavities, supporting our hypothesis that these sites are involved in the inhibition of aSyn toxic fibrillization. Finally, we expand our study to include analysis of fibril structures of tau, FUS, TDP-43, prion, and hnRNPA1; all of which contained hydrated cavities, with tau, FUS, and TDP-43 recapitulating our PTM results in aSyn fibril cavities.

Conformational fluctuations and phases in fused in sarcoma (FUS) low-complexity domain

The well-known phenomenon of phase separation in synthetic polymers and proteins has become a major topic in biophysics because it has been invoked as a mechanism of compartment formation in cells, without the need for membranes. Most of the coacervates (or condensates) are composed of Intrinsically Disordered Proteins (IDPs) or regions that are structureless, often in interaction with RNA and DNA. One of the more intriguing IDPs is the 526-residue RNA-binding protein, Fused in Sarcoma (FUS), whose monomer conformations and condensates exhibit unusual behavior that are sensitive to solution conditions. By focussing principally on the N-terminus low-complexity domain (FUS-LC comprising residues 1-214) and other truncations, we rationalize the findings of solid-state NMR experiments, which show that FUS-LC adopts a non-polymorphic fibril structure (core-1) involving residues 39-95, flanked by fuzzy coats on both the N- and C-terminal ends. An alternate structure (core-2), whose free energy is comparable to core-1, emerges only in the truncated construct (residues 110-214). Both core-1 and core-2 fibrils are stabilized by a Tyrosine ladder as well as hydrophilic interactions. The morphologies (gels, fibrils, and glass-like) adopted by FUS seem to vary greatly, depending on the experimental conditions. The effect of phosphorylation is site-specific. Simulations show that phosphorylation of residues within the fibril has a greater destabilization effect than residues that are outside the fibril region, which accords well with experiments. Many of the peculiarities associated with FUS may also be shared by other IDPs, such as TDP43 and hnRNPA2. We outline a number of problems for which there is no clear molecular explanation.

The Role of c-Abl Tyrosine Kinase in Brain and Its Pathologies

Differentiated status, low regenerative capacity and complex signaling make neuronal tissues highly susceptible to translating an imbalance in cell homeostasis into cell death. The high rate of neurodegenerative diseases in the elderly population confirms this. The multiple and divergent signaling cascades downstream of the various stress triggers challenge researchers to identify the central components of the stress-induced signaling pathways that cause neurodegeneration. Because of their critical role in cell homeostasis, kinases have emerged as one of the key regulators. Among kinases, non-receptor tyrosine kinase (Abelson kinase) c-Abl appears to be involved in both the normal development of neural tissue and the development of neurodegenerative pathologies when abnormally expressed or activated. However, exactly how c-Abl mediates the progression of neurodegeneration remains largely unexplored. Here, we summarize recent findings on the involvement of c-Abl in normal and abnormal processes in nervous tissue, focusing on neurons, astrocytes and microglial cells, with particular reference to molecular events at the interface between stress signaling, DNA damage, and metabolic regulation. Because inhibition of c-Abl has neuroprotective effects and can prevent neuronal death, we believe that an integrated view of c-Abl signaling in neurodegeneration could lead to significantly improved treatment of the disease.