Insights into the Atomistic Mechanisms of Phosphorylation in Disrupting Liquid-Liquid Phase Separation and Aggregation of the FUS Low-Complexity Domain

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.

BEMM-GEN: A Toolkit for Generating a Biomolecular Environment-Mimicking Model for Molecular Dynamics Simulation

Understanding the influence of the cellular environment on protein conformations is crucial for elucidating protein functions within living cells. In studies using molecular dynamics (MD) simulation, carbon nanotubes and hydrophobic cages have been widely used to emulate the cellular environment inside specific large biomolecules such as ribosome tunnels and chaperones. However, recent studies suggest that these uniform hydrophobic models may not adequately capture the environmental effects inside each biomolecule. Based on these facts, it is necessary to generate spherical and cylindrical models with varied chemical properties corresponding to the components within target biomolecules. We developed a biomolecular environment-mimicking model generator (BEMM-GEN) that generates spherical and cylindrical models with user-specified chemical properties and allows the integration of arbitrary protein conformations into the generated models. BEMM-GEN provides model and protein complex structures, along with the corresponding parameter files for MD simulation (AMBER and GROMACS), and users immediately run their MD simulation based on the generated input files. BEMM-GEN can be freely downloaded and installed via a Python package manager (pip install BEMM-gen). The source code files and a user manual for operation are provided on GitHub (https://github.com/y4suda/BEMM-GEN).

Influence of ALS-linked M337V mutation on the conformational ensembles of TDP-43321-340 peptide monomer and dimer

The transactive response (TAR) DNA/RNA-binding protein 43 (TDP-43) can self-assemble into both functional stress granules via liquid-liquid phase separation (LLPS) and pathogenic amyloid fibrillary aggregates that are closely linked to amyotrophic lateral sclerosis. Previous experimental studies reported that the low complexity domain (LCD) of TDP-43 plays an essential role in the LLPS and aggregation of the full-length protein, and it alone can also undergo LLPS to form liquid droplets mainly via intermolecular interactions in the 321-340 region. And the ALS-associated M337V mutation impairs LCD's LLPS and facilitates liquid-solid phase transition. However, the underlying atomistic mechanism is not well understood. Herein, as a first step to understand the M337V-caused LLPS disruption of TDP-43 LCD mediated by the 321-340 region and the fibrillization enhancement, we investigated the conformational properties of monomer/dimer of TDP-43321-340 peptide and its M337V mutant by performing extensive all-atom explicit-solvent replica exchange molecular dynamic simulations. Our simulations demonstrate that M337V mutation alters the residue regions with high helix/β-structure propensities and thus affects the conformational ensembles of both monomer and dimer. M337V mutation inhibits helix formation in the N-terminal Ala-rich region and the C-terminal mutation site region, while facilitating their long β-sheet formation, albeit with a minor impact on the average probability of both helix structure and β-structure. Further analysis of dimer system shows that M337V mutation disrupts inter-molecular helix-helix interactions and W334-W334 π-π stacking interactions which were reported to be important for the LLPS of TDP-43 LCD, whereas enhances the overall peptide residue-residue interactions and weakens peptide-water interactions, which is conducive to peptide fibrillization. This study provides mechanistic insights into the M337V-mutation-induced impairment of phase separation and facilitation of fibril formation of TDP-43 LCD.

Ribosome Tunnel Environment Drives the Formation of α-Helix during Cotranslational Folding

Protein conformations in cells are not solely determined by amino acid sequences; they also depend on cellular environments. For instance, the ribosome tunnel induces its specific α-helix formation during cotranslational folding. Owing to the link between these temporally α-helix and biological functions, the mechanism of α-helix formation inside the ribosome tunnel has been previously explored. Consequently, the conformational restrictions of the tunnel were considered one of the driving forces of α-helix formation. Conversely, the ribosomal tunnel environment, including its chemical properties, appears to influence the α-helix formation. However, a comprehensive analysis of the ribosome tunnel environment's impact on the α-helix formation has not been conducted yet due to challenges in experimentally controlling it. Therefore, as a new computational approach, we proposed a ribosome environment-mimicking model (REMM) based on the radius and components of the experimentally determined ribosome tunnel structures. Using REMM, we assessed the impact of the ribosome tunnel environment on α-helix formation. Herein, we employed carbon nanotubes (CNT) as a reference model alongside REMM because CNT reproduce conformational restrictions rather than the ribosome tunnel environment. Quantitatively, the ability to reproduce the α-helix of nascent peptides in the experimental structure was compared between the CNT and REMM using enhanced all-atom molecular dynamics simulations. Consequently, the REMM more accurately reproduced the α-helix of the nascent peptides than the CNT, highlighting the significance of the ribosome tunnel environment in α-helix formation. Additionally, we analyzed the properties of the peptide inside each model to reveal the mechanism of ribosome tunnel-specific α-helix formation. Consequently, we revealed that the chemical diversities of the tunnel are essential for the formation of backbone-to-backbone hydrogen bonds in the peptides. In conclusion, the ribosome tunnel environment, with the diverse chemical properties, drives its specific α-helix formation. By proposing REMM, we newly provide the technical basis for investigating the protein conformations in various cellular environments.

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.

Coarse-Grained Modeling of Liquid-Liquid Phase Separation in Cells: Challenges and Opportunities

Liquid-liquid phase separation (LLPS) within cells gives rise to membraneless organelles, which play pivotal roles in numerous cellular functions. A comprehensive understanding of the functional aspects of intrinsically disordered protein (IDP) condensates necessitates elucidating their inherent structures and establishing correlations with biological functions. Coarse-grained (CG) molecular dynamics (MD) simulations present a promising avenue for gaining insights into LLPS mechanisms of biomacromolecules. Essential to this endeavor is the development of tailored CG force fields for MD simulations, incorporating the full spectrum of biomolecules involved in the formation of condensates and accounting for real-time biochemical reactions coupled to the LLPS. Moreover, developing accurate theoretical frameworks and establishing links between condensate structure and its function are imperative for a thorough comprehension of LLPS of biological systems.

Molecular Insights into the Differential Effects of Acetylation on the Aggregation of Tau Microtubule-Binding Repeats

Aggregation of tau protein into intracellular fibrillary inclusions is characterized as the hallmark of tauopathies, including Alzheimer's disease and chronic traumatic encephalopathy. The microtubule-binding (MTB) domain of tau, containing either three or four repeats with sequence similarities, plays an important role in determining tau's aggregation. Previous studies have reported that abnormal acetylation of lysine residues displays a distinct effect on the formation of pathological tau aggregates. However, the underlying molecular mechanism remains mostly elusive. In this study, we performed extensive replica exchange molecular dynamics (REMD) simulations of 144 μs in total to systematically investigate the dimerization of four tau MTB repeats and explore the impacts of Lys280 (K280) or Lys321 (K321) acetylation on the conformational ensembles of the R2 or R3 dimer. Our results show that R3 is the most prone to aggregation among the four repeats, followed by R2 and R4, while R1 displays the weakest aggregation propensity with a disordered structure. Acetylation of K280 could promote the aggregation of R2 peptides by increasing the formation of β-sheet structures and strengthening the interchain interaction. However, K321 acetylation decreases the β-sheet content of the R3 dimer, reduces the ability of R3 peptides to form long β-strands, and promotes the stable helix structure formation. The salt bridge and Y310-Y310 π-π stacking interactions of the R3 dimer are greatly weakened by K321 acetylation, resulting in the inhibition of dimerization. This study uncovers the structural ensembles of tau MTB repeats and provides mechanistic insights into the influences of acetylation on tau aggregation, which may deepen the understanding of the pathogenesis of tauopathies.

Crosstalk between protein post-translational modifications and phase separation

The phenomenon of phase separation is quite common in cells, and it is involved in multiple processes of life activities. However, the current research on the correlation between protein modifications and phase separation and the interference with the tendency of phase separation has some limitations. Here we focus on several post-translational modifications of proteins, including protein phosphorylation modification at multiple sites, methylation modification, acetylation modification, ubiquitination modification, SUMOylation modification, etc., which regulate the formation of phase separation and the stability of phase separation structure through multivalent interactions. This regulatory role is closely related to the development of neurodegenerative diseases, tumors, viral infections, and other diseases, and also plays essential functions in environmental stress, DNA damage repair, transcriptional regulation, signal transduction, and cell homeostasis of living organisms, which provides an idea to explore the interaction between novel protein post-translational modifications and phase separation. Video Abstract.

Molecular mechanisms involved in the destabilization of two types of R3-R4 tau fibrils associated with chronic traumatic encephalopathy by Fisetin

Chronic traumatic encephalopathy is a neurodegenerative tauopathy pathologically characterized by fibrillary tau aggregates in the depth of sulci. Clearing fibrous tau aggregates is considered a promising strategy in the treatment of CTE. Fisetin (FS), a natural polyphenolic small molecule, was confirmed to disassociate the tau filaments in vitro. However, the molecular mechanisms of FS in destabilizing the CTE-related R3-R4 tau fibrils remain largely unknown. In this study, we compared the atomic-level structural differences of the two types of CTE-related R3-R4 tau fibrils and explored the influence and molecular mechanisms of FS on the two types of fibrils by conducting multiple molecular dynamics (MD) simulations. The results reveal that the type 1 fibril displays higher structural stability than the type 2 fibril, with a lower root-mean-square-fluctuation value and higher β-sheet structure probability. FS can destabilize both types of fibrils by decreasing the β-sheet structure content, interrupting the mainchain H-bond network, and increasing the solvent accessible surface area and β7-β8 angle of the fibrils. H-bonding, π-π stacking and cation-π are the common interactions driving FS molecules binding on the two types of fibrils, while the hydrophobic interaction occurs only in the type 2 fibril. Due to the relatively short simulation time, our study captures the early molecular mechanisms. However, it does provide beneficial information for the design of drugs to prevent or treat CTE.

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.

Machine Learning Subtle Conformational Change due to Phosphorylation in Intrinsically Disordered Proteins

Phosphorylation of intrinsically disordered proteins/regions (IDPs/IDRs) has a profound effect in biological functions such as cell signaling, protein folding or unfolding, and long-range allosteric effects. However, here we focus on two IDPs, namely 83-residue IDR transcription factor Ash1 and 92-residue long N-terminal region of CDK inhibitor Sic1 protein, found in Saccharomyces cerevisiae, for which experimental measurements of average conformational properties, namely, radius of gyration and structure factor, indicate negligible changes upon phosphorylation. Here, we show that a judicious dissection of conformational ensemble via combination of unsupervised machine learning and extensive molecular dynamics (MD) trajectories can highlight key differences and similarities among the phosphorylated and wild-type IDP. In particular, we develop Markov state model (MSM) using the latent-space dimensions of an autoencoder, trained using multi-microsecond long MD simulation trajectories. Examination of structural changes among the states, prior to and upon phosphorylation, captured several similarities and differences in their backbone contact maps, secondary structure, and torsion angles. Hydrogen bonding analysis revealed that phosphorylation not only increases the number of hydrogen bonds but also switches the pattern of hydrogen bonding between the backbone and side chain atoms with the phosphorylated residues. We also observe that although phosphorylation introduces salt bridges, there is a loss of the cation-π interaction. Phosphorylation also improved the probability for long-range hydrophobic contacts and also enhanced interaction with water molecules and improved the local structure of water as evident from the geometric order parameters. The observations on these machine-learnt states gave important insights, as it would otherwise be difficult to determine experimentally which is important, if we were to understand the role of phosphorylation of IDPs in their biological functions.

Benchmarking of force fields to characterize the intrinsically disordered R2-FUS-LC region

Intrinsically Disordered Proteins (IDPs) play crucial roles in numerous diseases like Alzheimer's and ALS by forming irreversible amyloid fibrils. The effectiveness of force fields (FFs) developed for globular proteins and their modified versions for IDPs varies depending on the specific protein. This study assesses 13 FFs, including AMBER and CHARMM, by simulating the R2 region of the FUS-LC domain (R2-FUS-LC region), an IDP implicated in ALS. Due to the flexibility of the region, we show that utilizing multiple measures, which evaluate the local and global conformations, and combining them together into a final score are important for a comprehensive evaluation of force fields. The results suggest c36m2021s3p with mTIP3p water model is the most balanced FF, capable of generating various conformations compatible with known ones. In addition, the mTIP3P water model is computationally more efficient than those of top-ranked AMBER FFs with four-site water models. The evaluation also reveals that AMBER FFs tend to generate more compact conformations compared to CHARMM FFs but also more non-native contacts. The top-ranking AMBER and CHARMM FFs can reproduce intra-peptide contacts but underperform for inter-peptide contacts, indicating there is room for improvement.

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.

E46K Mutation of α-Synuclein Preorganizes the Intramolecular Interactions Crucial for Aggregation

Aggregation of α-synuclein is central to the pathogenesis of Parkinson's disease. The most toxic familial mutation E46K accelerates the aggregation process by an unknown mechanism. Herein, we provide a clue by investigating the influence of E46K on monomeric α-synuclein and its relation to aggregation with molecular dynamics simulations. The E46K mutation suppresses β-sheet structures in the N-terminus while promoting those at the key fibrillization region named NACore. Even though WT and E46K monomers share conserved intramolecular interactions with fibrils, E46K abolishes intramolecular contacts within the N-terminus which are present in the WT monomer but absent in fibrils. Network analysis identifies residues 36-53 as the interaction core of the WT monomer. Upon mutation, residues 36-46 are expelled to water due to aggravated electrostatic repulsion in the 43KTKK46 segment. Instead, NACore (residues 68-78) becomes the interaction hub and connects preceding residues 47-56 and the C-terminus. Consequently, residues 47-95 which belong to the fibril core form more compact β-sheets. Overall, the interaction network of E46K is more like fibrils than WT, stabilizing the fibril-like conformations. Our work provides mechanistic insights into the faster aggregation of the E46K mutant. It implies a close link between monomeric conformations and fibrils, which would spur the development of therapeutic strategies.

Dissecting the Effect of ALS Mutation G335D on the Early Aggregation of the TDP-43 Amyloidogenic Core Peptide: Helix-to-β-Sheet Transition and Conformational Shift

The aggregation of TAR DNA-binding protein of 43 kDa (TDP-43) into fibrillary deposits is associated with amyotrophic lateral sclerosis (ALS). The 311-360 fragment of TDP-43 (TDP-43311-360), the amyloidogenic core region, can spontaneously aggregate into fibrils, and the ALS-associated mutation G335D has an enhanced effect on TDP-43311-360 fibrillization. However, the molecular mechanism underlying G335D-enhanced aggregation at atomic level remains largely unknown. By utilizing all-atom molecular dynamics (MD) and replica exchange with solute tempering 2 (REST2) simulations, we investigated influences of G335D on the dimerization (the first step of aggregation) and conformational ensemble of the TDP-43311-360 peptide. Our simulations show that G335D mutation increases inter-peptide interactions, especially inter-peptide hydrogen-bonding interactions in which the mutant site has a relatively large contribution, and enhances the dimerization of TDP-43311-360 peptides. The α-helix regions in the NMR-resolved conformation of the TDP-43311-360 monomer (321-330 and 335-343) play an essential role in the formation of the dimer. G335D mutation induces helix unfolding and promotes α-to-β conversion. G335D mutation alters the conformational distribution of TDP-43311-360 dimers and causes population shift from helix-rich to β-sheet-rich conformations, which facilitates the fibrillization of the TDP-43311-360 peptide. Our MD and REST2 simulation results suggest that the 321-330 region is of paramount importance to α-to-β transition and could be the initiation site for TDP-43311-360 fibrillization. Our work reveals the mechanism underlying the enhanced aggregation propensity of the G335D TDP-43311-360 peptide, which provides atomistic insights into the G335D mutation-caused pathogenicity of TDP-43 protein.

Dissecting the Molecular Mechanisms of the Co-Aggregation of Aβ40 and Aβ42 Peptides: A REMD Simulation Study

The aggregation of amyloid-β protein (Aβ) into oligomers and amyloid fibrils is closely related to Alzheimer's disease (AD). Aβ40 and Aβ42, as two most prominent isoforms of Aβ peptides, can cross-interact with each other and form co-aggregates, which affect the progression of the disease. However, the molecular determinants underlying Aβ40 and Aβ42 cross-interaction and the structural details of their co-oligomers remain elusive. Herein, we performed all-atom explicit-solvent replica exchange molecular dynamics simulations on Aβ40-Aβ42 heterogeneous and Aβ40/Aβ42 homogeneous dimer systems to dissect the co-aggregation mechanisms of the two isoforms. Our results show that the interpeptide main-chain interaction of Aβ40-Aβ42 is stronger than that of Aβ40-Aβ40 and Aβ42-Aβ42. The positions of hotspot residues in heterodimers and homodimers display high similarity, implying similar molecular recognition sites for both cross-interaction and self-interaction. Contact maps of Aβ40-Aβ42 heterodimers reveal that residue pairs crucial for cross-interaction are mostly located in the C-terminal hydrophobic regions of Aβ40 and Aβ42 peptides. Conformational analysis shows that Aβ40 and Aβ42 monomers can co-assemble into β-sheet-rich heterodimers with shorter β-sheets than those in homodimers, which is decremental to monomer addition. Similar molecular recognition sites and β-sheet distribution of Aβ40 and Aβ42 peptides are observed in heterodimers and homodimers, which may provide the molecular basis for the two isoforms' co-aggregation and cross-seeding. Our work dissects the co-aggregation mechanisms of Aβ40 and Aβ42 peptides at the atomic level, which will help for in-depth understanding of the cross-talk between the two Aβ isoforms and the pathogenesis of AD.

Dissecting how ALS-associated D290V mutation enhances pathogenic aggregation of hnRNPA2286-291 peptides: Dynamics and conformational ensembles

The aggregation of RNA binding proteins, including hnRNPA1/2, TDP-43 and FUS, is heavily implicated in causing or increasing disease risk for a series of neurodegenerative diseases such as amyotrophic lateral sclerosis (ALS). A recent experimental study demonstrated that an ALS-related D290V mutation in the low complexity domain (LCD) of hnRNPA2 can enhance the aggregation propensity of wild type (WT) hnRNPA2_286-291 peptide. However, the underlying molecular mechanisms remain elusive. Herein, we investigated effects of D290V mutation on aggregation dynamics of hnRNPA2_286-291 peptide and the conformational ensemble of hnRNPA2_286-291 oligomers by performing all-atom molecular dynamic and replica-exchange molecular dynamic simulations. Our simulations demonstrate that D290V mutation greatly reduces the dynamics of hnRNPA2_286-291 peptide and that D290V oligomers possess higher compactness and β-sheet content than WT, indicative of mutation-enhanced aggregation capability. Specifically, D290V mutation strengthens inter-peptide hydrophobic, main-chain hydrogen bonding and side-chain aromatic stacking interactions. Those interactions collectively lead to the enhancement of aggregation capability of hnRNPA2_286-291 peptides. Overall, our study provides insights into the dynamics and thermodynamic mechanisms underlying D290V-induced disease-causing aggregation of hnRNPA2_286-291, which could contribute to better understanding of the transitions from reversible condensates to irreversible pathogenic aggregates of hnRNPA2 LCD in ALS-related diseases.

Glutathionylation on RNA-binding proteins: a regulator of liquid‒liquid phase separation in the pathogenesis of amyotrophic lateral sclerosis

RNA-binding proteins (RBPs) containing low-sequence complexity domains mediate the formation of cellular condensates and membrane-less organelles with biological functions via liquid‒liquid phase separation (LLPS). However, the abnormal phase transition of these proteins induces the formation of insoluble aggregates. Aggregates are pathological hallmarks of neurodegenerative diseases, such as amyotrophic lateral sclerosis (ALS). The molecular mechanisms underlying aggregate formation by ALS-associated RPBs remain largely unknown. This review highlights emerging studies on various posttranslational modifications (PTMs) related to protein aggregation. We begin with the introduction of several ALS-associated RBPs that form aggregates induced by phase separation. In addition, we highlight our recent discovery of a new PTM involved in the phase transition during the pathogenesis of fused-in-sarcoma (FUS)-associated ALS. We suggest a molecular mechanism through which LLPS mediates glutathionylation in FUS-linked ALS. This review aims to provide a detailed overview of the key molecular mechanisms of LLPS-mediated aggregate formation by PTMs, which will help further the understanding of the pathogenesis and development of ALS therapeutics.

Effects of Mass Change on Liquid–Liquid Phase Separation of the RNA-Binding Protein Fused in Sarcoma

In recent years, many experimental and theoretical studies of protein liquid–liquid phase separation (LLPS) have shown its important role in the processes of physiology and pathology. However, there is a lack of definite information on the regulation mechanism of LLPS in vital activities. Recently, we found that the intrinsically disordered proteins with the insertion/deletion of a non-interacting peptide segment or upon isotope replacement could form droplets, and the LLPS states are different from the proteins without those. We believed that there is an opportunity to decipher the LLPS mechanism with the mass change perspective. To investigate the effect of molecular mass on LLPS, we developed a coarse-grained model with different bead masses, including mass 1.0, mass 1.1, mass 1.2, mass 1.3, and mass 1.5 in atomic units or with the insertion of a non-interacting peptide (10 aa) and performed molecular dynamic simulations. Consequently, we found that the mass increase promotes the LLPS stability, which is based on decreasing the z motion rate and increasing the density and the inter-chain interaction of droplets. This insight into LLPS by mass change paves the way for the regulation and relevant diseases on LLPS.

ALS-Linked A315T and A315E Mutations Enhance β-Barrel Formation of the TDP-43307-319 Hexamer: A REST2 Simulation Study

Pathogenic mutations of transactivation response element DNA-binding protein 43 (TDP-43) are closely linked with amyotrophic lateral sclerosis (ALS). It was recently reported that two ALS-linked familial mutants A315T and A315E of TDP-43_307-319 peptides can self-assemble into oligomers including tetramers, hexamers, and octamers, among which hexamers were suggested to form the β-barrel structure. However, due to the transient nature of oligomers, their conformational properties and the atomic mechanisms underlying the β-barrel formation remain largely elusive. Herein, we investigated the hexameric conformational distributions of the wild-type (WT) TDP-43_307-319 fragment and its A315T and A315E mutants by performing all-atom explicit-solvent replica exchange with solute tempering 2 simulations. Our simulations reveal that each peptide can self-assemble into diverse conformations including ordered β-barrels, bilayer β-sheets and/or monolayer β-sheets, and disordered complexes. A315T and A315E mutants display higher propensity to form β-barrel structures than the WT, which provides atomic explanation for their enhanced neurotoxicity reported previously. Detailed interaction analysis shows that A315T and A315E mutations increase inter-molecular interactions. Also, the β-barrel structures formed by the three different peptides are stabilized by distinct inter-peptide side-chain hydrogen bonding, hydrophobic, and aromatic stacking interactions. This study demonstrates the enhanced β-barrel formation of the TDP-43_307-319 hexamer by the pathogenic A315T and A315E mutations and reveals the underlying molecular determinants, which may be helpful for in-depth understanding of the ALS-mutation-induced neurotoxicity of TDP-43 protein.

Atomistic Insights into the Inhibitory Mechanism of Tyrosine Phosphorylation against the Aggregation of Human Tau Fragment PHF6

Abnormal aggregation of the microtubule-associated protein tau into intracellular fibrillary inclusions is characterized as the hallmark of tauopathies, including Alzheimer's disease and chronic traumatic encephalopathy. The hexapeptide ³⁰⁶VQIVYK³¹¹ (PHF6) of R3 plays an important role in the aggregation of tau. Recent experimental studies reported that phosphorylation of residue tyrosine 310 (Y310) could decrease the propensity of PHF6 to form fibrils and inhibit tau aggregation. However, the underlying inhibitory mechanism is not well understood. In this work, we systematically investigated the influences of phosphorylation on the conformational ensembles and oligomerization dynamics of PHF6 by performing extensive all-atom molecular dynamics (MD) simulations. Our replica exchange MD simulations demonstrate that Y310 phosphorylation could effectively suppress the formation of β-structure and shift PHF6 oligomers toward coil-rich aggregates. The interaction analyses show that hydrogen bonding and hydrophobic interactions among PHF6 peptides, as well as Y310-Y310 π-π stacking and I308-Y310 CH-π interactions, are weakened by phosphorylation. Additional microsecond MD simulations show that Y310 phosphorylation could inhibit the oligomerization of PHF6 by preventing the formation of large β-sheet oligomers and multi-layer β-sheet aggregates. This study provides mechanistic insights into the phosphorylation-inhibited tau aggregation, which may be helpful for the in-depth understanding of the pathogenesis of tauopathies.

Atomistic Insights into A315E Mutation-Enhanced Pathogenicity of TDP-43 Core Fibrils

The aggregation of TAR DNA-binding protein of 43 kDa (TDP-43) into fibrillary deposits is implicated in amyotrophic lateral sclerosis (ALS), and some hereditary mutations localized in the low complexity domain (LCD) facilitate the formation of pathogenic TDP-43 fibrils. A recent cryo-EM study reported the atomic-level structures of the A315E TDP-43 LCD (residues 288-319, TDP-43_288-319) core fibril in which the protofilaments have R-shaped structures and hypothesized that A315E U-shaped protofilaments can readily convert to R-shaped protofilaments compared to the wild-type (WT) ones. There are no atomic structures of WT protofilaments available yet. Herein, we performed extensive all-atom explicit-solvent molecular dynamics simulations on A315E and WT protofilaments starting from both the cryo-EM-determined R-shaped and our constructed U-shaped structures. Our simulations show that WT protofilaments also adopt the R-shaped structures but are less stable than their A315E counterparts. Except for R293-E315 salt bridges, N312-F316 hydrophobic interactions and F316-F316 π-π stacking interactions are also crucial for the stabilization of the neck region of the R-shaped A315E protofilaments. The loss of R293-E315 salt bridges and the weakened interactions of N312-F316 and F316-F316 result in the reduced stability of the R-shaped WT protofilaments. Simulations starting from U-shaped folds reveal that A315E protofilaments can spontaneously convert to the cryo-EM-derived R-shaped protofilaments, whereas WT protofilaments convert to R-shape-like structures with remodeled neck regions. The R-shape-like WT protofilaments could act as intermediate states slowing down the U-to-R transition. This study reveals that A315E mutation can not only enhance the structural stability of the R-shaped TDP-43_288-319 protofilaments but also promote the U-to-R transition, which provides atomistic insights into the A315E mutation-enhanced TDP-43 pathogenicity in ALS.